A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Phylogenomics of eukaryotes: Impact of missing data on large alignments

Resolving the relationships between Metazoa and other eukaryotic groups as well as between metazoan phyla is central to the understanding of the origin and evolution of animals. The current view is based on limited data sets, either a single gene with many species (e.g., ribosomal RNA) or many genes but with only a few species. Because a reliable phylogenetic inference simultaneously requires numerous genes and numerous species, we assembled a very large data set containing 129 orthologous proteins (similar to30,000 aligned amino acid positions) for 36 eukaryotic species. Included in the alignments are data from the choanoflagellate Monosiga ovata, obtained through the sequencing of about 1,000 cDNAs. We provide conclusive support for choanoflagellates as the closest relative of animals and for fungi as the second closest. The monophyly of Plantae and chromalveolates was recovered but without strong statistical support. Within animals, in contrast to the monophyly of Coelomata observed in several recent large-scale analyses, we recovered a paraphyletic Coelamata, with nematodes and platyhelminths nested within. To include a diverse sample of organisms, data from EST projects were used for several species, resulting in a large amount of missing data in our alignment (about 25%). By using different approaches, we verify that the inferred phylogeny is not sensitive to these missing data. Therefore, this large data set provides a reliable phylogenetic framework for studying eukaryotic and animal evolution and will be easily extendable when large amounts of sequence information become available from a broader taxonomic range.

A novel method for the transport and analysis of genetic material from polyps and zooxanthellae of scleractinian corals

We have developed a new simple method for transport, storage, and analysis of genetic material from the corals Agaricia agaricites, Dendrogyra cylindrica, Eusmilia ancora, Meandrina meandrites, Montastrea annularis, Porites astreoides, Porites furcata, Porites porites, and Siderastrea siderea at room temperature. All species yielded sufficient DNA from a single FTA(R) card (19 mug-43 ng) for subsequent PCR amplification of both coral and zooxanthellar DNA. The D1 and D2 variable region of the large Subunit rRNA gene (LSUrDNA) was amplified from the DNA of P. furcata and S. siderea by PCR. Electrophoresis yielded two major DNA bands: an 800-base pair (bp) DNA, which represented the coral ribosomal RNA (rRNA) gene, and a 600-bp DNA, which represented the zooxanthellar srRNA gene. Extraction of DNA from the bands yielded between 290 mug total DNA (S. siderea coral DNA) and 9 mug total DNA (P. furcata zooxanthellar DNA). The ability to transport and store genetic material from scleractinian corals without resort to laboratory facilities in the field allows for the molecular Study of a far wider range and variety of coral sites than have been studied to date. (C) 2003 Elsevier Science B.V. All rights reserved.

Development of a PCR-based assay for rapid and reliable identification of pathogenic Fusaria

Identification of Fusarium species has always been difficult due to confusing phenotypic classification systems. We have developed a fluorescent-based polymerase chain reaction assay that allows for rapid and reliable identification of five toxigenic and pathogenic Fusarium species. The species includes Fusarium avenaceum, F. culmorum, F. equiseti, F. oxysporum and F. sambucinum. The method is based on the PCR amplification of species-specific DNA fragments using fluorescent oligonucleotide primers, which were designed based on sequence divergence within the internal transcribed spacer region of nuclear ribosomal DNA. Besides providing an accurate, reliable, and quick diagnosis of these Fusaria, another advantage with this method is that it reduces the potential for exposure to carcinogenic chemicals as it substitutes the use of fluorescent dyes in place of ethidium, bromide. Apart from its multidisciplinary importance and usefulness, it also obviates the need for gel electrophoresis. (C) 2002 Published by Elsevier Science B.V. on behalf of the Federation of European Microbiological Societies.

Biodiversity of human faecal bacteria isolated from phytic acid enriched chemostat fermenters

Background: Myo-inositol hexaphosphate (IP6) or phytic acid is found mostly in cereals and legumes and is thought to possess anti-carcinogenic properties. Aim: To isolate and identify faecal bacteria capable of phytic acid metabolism and to assess the effectiveness of prebiotics (dietary oligosaccharides, metabolised by selective colonic bacteria) in preserving the integrity of phytic acid. Methods: Faecal samples from three volunteers were used in continuous culture experiments under varying conditions of pH, substrate concentration and dilution rates, seventy three different isolates cultured at steady state were then screened for phytic acid metabolism and identified through partial sequencing of their 16S rRNA genes (16S ribosomal ribonucleic acid). Utilisation of phytic acid was also assessed in a continuous culture system enriched with prebiotic fructooligosaccharides (FOS). Results: Bacteroides spp., Clostridium spp. and facultatively anaerobic bacteria generally appeared to maintain viable counts in the presence of phytic acid. Bifidobacterium spp. and Lactobacillus spp. appeared less able to maintain viable counts in the presence of phytic acid. These results were confirmed by an increase in viable counts of Bacteroides spp., Clostridium spp. and a decrease in viable counts of Bifidobacterium spp. and Lactobacillus spp. once phytic acid was introduced to a FOS enriched continuous culture. Conclusions: The phytate metabolising biodiversity from the human large intestine does not appear to encompass major bacterial genera associated with beneficial or benign health effects (e.g. Lactobacillus spp. and Bifidobacterium spp).

Reclassification of Bacteroides putredinis (Weinberg et al., 1937) in a New Genus Alistipes gen. nov., as Alistipes putredinis comb. nov., and description of Alistipes finegoldii sp nov., from human sources

During studies on the bacteriology of appendicitis in children, we often isolated from inflamed and non-inflamed tissue samples, an unusual bile-resistant pigment-producing strictly anaerobic gram-negative rod. Phenotypically this organism resembles members of Bacteroides fragilis group of species, as it is resistant to bile and exhibits a special-potency-disk pattern (resistance to vancomycin, kanamycin and colistin) typical for the B. fragilis group. However, the production of brown pigment on media containing haemolysed blood and a cellular fatty acid composition dominated by iso-C15:0, suggests that the organism most closely resembles species of the genus Porphyromonas. However, the unidentified organism differs from porphyromonads by being bile-resistant and by not producing butyrate as a metabolic end-product. Comparative 16S ribosomal RNA gene sequencing studies show the unidentified organism represents a distinct sub-line, associated with but distinct from, the miss-classified species Bacteroides putredinis. The clustering of the unidentified bacterium with Bacteroides putredinis was statistically significant, but they displayed >4% sequence divergence with each other. Chromosomal DNA-DNA pairing studies further confirmed the separateness of the unidentified bacterium and Bacteroides putredinis. Based on phenotypic and phylogenetic considerations, it is proposed that Bacteroides putredinis and the unidentified bacterium from human sources be classified in a new genus Alistipes, as Alistipes putredinis comb. nov. and Alistipes finegoldii sp. nov., respectively. The type strain of Alistipes finegoldii is CCUG 46020(T) (= AHN2437(T)).

Variation in 16S-23S rRNA intergenic spacer regions in Photobacterium damselae: a mosaic-like structure

Phenotypically, Photobacterium damselae subsp. piscicida and P. damselae subsp. damselae are easily distinguished. However, their 16S rRNA gene sequences are identical, and attempts to discriminate these two subspecies by molecular tools are hampered by their high level of DNA-DNA similarity. The 16S-23S rRNA internal transcribed spacers (ITS) were sequenced in two strains of Photobacterium damselae subsp. piscicida and two strains of P. damselae subsp. damselae to determine the level of molecular diversity in this DNA region. A total of 17 different ITS variants, ranging from 803 to 296 bp were found, some of which were subspecies or strain specific. The largest ITS contained four tRNA genes (tDNAs) coding for tRNA(Glu(UUC)), tRNA(LyS(UUU)), tRNA(Val(UAC)), and tRNA(Ala(GGC)). Five amplicons contained tRNA(Glu(UUC)) combined with two additional tRNA genes, including tRNA(Lys(UUU)), tRNA(Val(UAC)), or tRNA(Ala(UGC)). Five amplicons contained tRNA(Ile(GAU)) and tRNA(Ala(UGC)). Two amplicons contained tRNA(Glu(UUC)) and tRNA(Val(UGC)). Two different isoacceptor tRNA(Ala) genes (GGC and UGC anticodons) were found. The five smallest amplicons contained no tRNA genes. The tRNA-gene combinations tRNA(Glu(UUC)) -tRNA(Val(UAC)) -tRNA(Ala(UGC)) and tRNA(Glu(UUC)) -tRNA(Ala(UGC)) have not been previously reported in bacterial ITS regions. The number of copies of the ribosomal operon (rrn) in the P. damselae chromosome ranged from at least 9 to 12. For ITS variants coexisting in two strains of different subspecies or in strains of the same subspecies, nucleotide substitution percentages ranged from 0 to 2%. The main source of variation between ITS variants was due to different combinations of DNA sequence blocks, constituting a mosaic-like structure.

Characterization of the 23S and 5S rRNA genes and 23S-5S intergenic spacer region (ITS-2) of Photobacterium damselae

The 23S ribosomal RNA (rRNA) gene has been sequenced in strains of the fish pathogens Photobacterium damselae subsp. damselae (ATCC 33539) and subsp. piscicida (ATCC 29690), showing that 3 nucleotide positions are clearly different between subspecies. In addition, the 5S rRNA gene plus the intergenic spacer region between the 23S and 5S rRNA genes (ITS-2) were amplified, cloned and sequenced for the 2 reference strains as well as the field isolates RG91 (subsp. damselae) and DI21 (subsp. piscicida). A 100% similarity was found for the consensus 5S rRNA gene sequence in the 2 subspecies, although some microheterogeneity was detected as inter-cistronic variability within the same chromosome. Sequence analysis of the spacer region between the 23S and 5S rRNA genes revealed 2 conserved and 3 variable nucleotide sequence blocks, and 4 different modular organizations were found. The ITS-2 spacer region exhibited both inter-subspecies and inter-cistronic polymorphism, with a mosaic-like structure. The EMBL accession numbers for the 23S, 5S and ITS-2 sequences are: P. damselae subsp. piscicida 5S gene (AJ274379), P. damselae subsp. damselae 23S gene (Y18520), subsp. piscicida 23S gene (Y17901), R damselae subsp. piscicida ITS-2 (AJ250695, AJ250696), P. damselae subsp. damselae ITS-2 (AJ250697, AJ250698).

Molecular inventory of faecal microflora in patients with Crohn's disease

Intestinal microbial community is involved in the pathogenesis of Crohn's disease, but knowledge of its potential abnormalities has been limited by the impossibility to grow many dominant intestinal bacteria. Using sequence analysis of randomly cloned bacterial 16S ribosomal DNA, the dominant faecal species from four Crolin's disease patients and four controls were compared. Whereas marked inter-individual differences were observed in the faecal microflora of patients, three remained distantly related to controls on the basis of their operational taxonomic unit composition. Bacteroides vidgatus and closely related organisms represented the only molecular species shared by all patients and exhibited an unusually high rate of occurrence. Escherichia coli clones were isolated only in two patients with ileocolonic Crohn's disease. Moreover, numerous clones belonged to phylogenetic groups or species that are commonly not dominant in the faecal microflora of healthy subjects: Pectinatus, Sutterella, Verritcomicrobium, Fusobacterium, Clostridium disporicum, clostridium glycolicum, Clostridium ramosum, Clostridium innocuum and Clostridium perfringens. (C) 2004 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.

Subdoligranulum variabile gen. nov., sp nov from human feces

During studies on the microflora of human feces we have isolated a strictly anaerobic, non-spore-forming, Gram-negative staining organism which exhibits a somewhat variable coccus-shaped morphology. Comparative 16S ribosomal RNA gene sequencing studies show the unidentified organism is phylogenetically a member of the Clostridium leptum supra-generic rRNA cluster and displays a close affinity to some rDNA clones derived from human and pig feces. The nearest named relatives of the unidentified isolate corresponded to Faecalibacterium prausnitzii (formerly Fusobacterium prausnitzii) displaying a 16S rRNA sequence divergence of approximately 9%, with Anaerofilum agile and A. pentosovorans the next closest relatives of the unidentified bacterium (sequence divergence approximately 10%). Based on phenotypic and phylogenetic considerations, it is proposed that the unusual coccoid-shaped organism be classified as a new genus and species, Subdoligranulum variabile. The type strain of S. variabile is BI 114(T) (= CCUG 47106(T) = DSM 15176(T)). (C) 2004 Elsevier Ltd. All rights reserved.

Characterization of some Actinomyces-like isolates from human clinical sources: description of Varibaculum cambriensis gen. nov., sp nov

Fifteen strains of an anaerobic, catalase-negative, gram-positive diphtheroid-shaped bacterium recovered from human sources were characterized by phenotypic and molecular chemical and molecular genetic methods. The unidentified bacterium showed some resemblance to Actinomyces species and related taxa, but biochemical testing, polyacrylamide gel electrophoresis analysis of whole-cell proteins, and amplified 16S ribosomal DNA restriction analysis indicated the strains were distinct from all currently named Actinomyces species and related taxa. Comparative 16S rRNA gene sequencing studies showed that the bacterium represents a hitherto-unknown phylogenetic line that is related to but distinct from Actinomyces, Actinobaculum, Arcanobacterium, and Mobiluncus. We propose, on the basis of phenotypic and phylogenetic evidence, that the unknown bacterium from human clinical specimens should be classified as a new genus and species, Varibaculum cambriensis gen. nov., sp. nov. The type strain of Varibaculum cambriensis sp. nov. is CCUG 44998(T) = CIP 107344(T).

Gene transcription in Daphnia magna: effects of acute exposure to a carbamate insecticide and an acetanilide herbicide.

Daphnia magna is a key invertebrate in the freshwater environment and is used widely as a model in ecotoxicological measurements and risk assessment. Understanding the genomic responses of D. magna to chemical challenges will be of value to regulatory authorities worldwide. Here we exposed D. magna to the insecticide methomyl and the herbicide propanil to compare phenotypic effects with changes in mRNA expression levels. Both pesticides are found in drainage ditches and surface water bodies standing adjacent to crops. Methomyl, a carbamate insecticide widely used in agriculture, inhibits acetylcholinesterase, a key enzyme in nerve transmission. Propanil, an acetanilide herbicide, is used to control grass and broad-leaf weeds. The phenotypic effects of single doses of each chemical were evaluated using a standard immobilisation assay. Immobilisation was linked to global mRNA expression levels using the previously estimated 48h-EC(1)s, followed by hybridization to a cDNA microarray with more than 13,000 redundant cDNA clones representing >5000 unique genes. Following exposure to methomyl and propanil, differential expression was found for 624 and 551 cDNAs, respectively (one-way ANOVA with Bonferroni correction, Pribosomal proteins, transcription factors). Methomyl induced the transcription of genes involved in specific processes such as ion homeostasis and xenobiotic metabolism. Propanil highly promoted haemoglobin synthesis and up-regulated genes specifically related to defence mechanisms (e.g., innate immunity response systems) and neuronal pathways. Pesticide-specific toxic responses were found but there is little evidence for transcriptional responses purely restricted to genes associated with the pesticide target site or mechanism of toxicity.

Examination of faecal Bifidobacterium populations in breast- and formula-fed infants during the first 18 months of life.

Bifidobacteria in the infant faecal microbiota have been the focus of much interest, especially during the exclusive milk-feeding period and in relation to the fortification of infant formulae to better mimic breast milk. However, longitudinal studies examining the diversity and dynamics of the Bifidobacterium population of infants are lacking, particularly in relation to the effects of weaning. Using a polyphasic strategy, the Bifidobacterium populations of breast- and formula-fed infants were examined during the first 18 months of life. Bifidobacterium-specific denaturing gradient gel electrophoresis demonstrated that breast-fed infants harboured greater diversity than formula-fed infants and the diversity of the infants' Bifidobacterium populations increased with weaning. Twenty-seven distinctive banding profiles were observed from ∼1100 infant isolates using ribosomal intergenic spacer analysis, 14 biotypes of which were confirmed to be members of the genus Bifidobacterium. Two profiles (H, Bifidobacterium longum subsp. infantis; and I, Bifidobacterium bifidum) were common culturable biotypes, seen in 9/10 infants, while profile E (Bifidobacterium breve) was common among breast-fed infants. Overall, inter- and intra-individual differences were observed in the Bifidobacterium populations of infants between 1 and 18 months of age, although weaning was associated with increased diversity of the infant Bifidobacterium populations. Breast-fed infants generally harboured a more complex Bifidobacterium microbiota than formula-fed infants.

Regulation of the cardiomyocyte transcriptome vs translatome by endothelin-1 and insulin: translational regulation of 5' terminal oligopyrimidine tract (TOP) mRNAs by insulin

Background: Changes in cellular phenotype result from underlying changes in mRNA transcription and translation. Endothelin-1 stimulates cardiomyocyte hypertrophy with associated changes in mRNA/protein expression and an increase in the rate of protein synthesis. Insulin also increases the rate of translation but does not promote overt cardiomyocyte hypertrophy. One mechanism of translational regulation is through 5' terminal oligopyrimidine tracts (TOPs) that, in response to growth stimuli, promote mRNA recruitment to polysomes for increased translation. TOP mRNAs include those encoding ribosomal proteins, but the full panoply remains to be established. Here, we used microarrays to compare the effects of endothelin-1 and insulin on the global transcriptome of neonatal rat cardiomyocytes, and on mRNA recruitment to polysomes (i.e. the translatome). Results: Globally, endothelin-1 and insulin (1 h) promoted >1.5-fold significant (false discovery rate < 0.05) changes in expression of 341 and 38 RNAs, respectively. For these transcripts with this level of change there was little evidence of translational regulation. However, 1336 and 712 RNAs had >1.25-fold significant changes in expression in total and/or polysomal RNA induced by endothelin-1 or insulin, respectively, of which ~35% of endothelin-1-responsive and ~56% of insulin-responsive transcripts were translationally regulated. Of mRNAs for established proteins recruited to polysomes in response to insulin, 49 were known TOP mRNAs with a further 15 probable/possible TOP mRNAs, but 49 had no identifiable TOP sequences or other consistent features in the 5' untranslated region. Conclusions: Endothelin-1, rather than insulin, substantially affects global transcript expression to promote cardiomyocyte hypertrophy. Effects on RNA recruitment to polysomes are subtle, with differential effects of endothelin-1 and insulin on specific transcripts. Furthermore, although insulin promotes recruitment of TOP mRNAs to cardiomyocyte polysomes, not all recruited mRNAs are TOP mRNAs.

Phylogenetic systematics of Sternbergia (Amaryllidaceae) based on plastid and ITS sequence data

The phylogenetics of Sternbergia (Amaryllidaceae) were studied using DNA sequences of the plastid ndhF and matK genes and nuclear internal transcribed spacer (ITS) ribosomal region for 38, 37 and 32 ingroup and outgroup accessions, respectively. All members of Sternbergia were represented by at least one accession, except S. minoica and S. schubertii, with additional taxa from Narcissus and Pancratium serving as principal outgroups. Sternbergia was resolved and supported as sister to Narcissus and composed of two primary subclades: S. colchiciflora sister to S. vernalis, S. candida and S. clusiana, with this clade in turn sister to S. lutea and its allies in both Bayesian and bootstrap analyses. A clear relationship between the two vernal flowering members of the genus was recovered, supporting the hypothesis of a single origin of vernal flowering in Sternbergia. However, in the S. lutea complex, the DNA markers examined did not offer sufficient resolving power to separate taxa, providing some support for the idea that S. sicula and S. greuteriana are conspecific with S. lutea