Phylogenetic analyses of ribosomal DNA-containing bacterioplankton genome fragments from a 4000 m vertical profile in the North Pacific Subtropical Gyre

High-throughput identification of rRNA gene-containing clones in large insert metagenomic libraries is difficult, because of the high background of host ribosomal RNA (rRNA) and rRNA genes. To address this challenge, a membrane hybridization method was developed to identify all bacterial small subunit rRNA-containing fosmid clones of microbial community DNA from seven different depths in the North Pacific Subtropical Gyre. Out of 101,376 clones screened, 751 rDNA-containing clones were identified that grouped in ∼60 different clades. Several rare sequences only remotely related to known groups were detected, including a Wolbachia -related sequence containing a putative intron or intervening sequence, as well as seven sequences from Order Myxococcales not previously detected in pelagic habitats. Stratified, depth-specific population structure was evident within both cultured and uncultured lineages. Conversely, some eurybathyal members of the genera Alcanivorax and Rhizobium shared identical small subunit ribosomal DNA sequences that were distributed from surface waters to the 4000 m depth. Comparison with similar analyses in Monterey Bay microbial communities revealed previously recognized, as well as some distinctive, depth-stratified partitioning that distinguished coastal from open ocean bacterioplankton populations. While some bias was evident in fosmid clone recovery in a few particular lineages, the overall phylogenetic group recovery and distributions were consistent with previous studies, as well as with direct shotgun sequence data from the same source DNA.     

Archaeal phylogeny based on ribosomal proteins

Until recently, phylogenetic analyses of Archaea have mainly been based on ribosomal RNA (rRNA) sequence comparisons, leading to the distinction of the two major archaeal phyla: the Euryarchaeota and the Crenarchaeota. Here, thanks to the recent sequencing of several archaeal genomes, we have constructed a phylogeny based on the fusion of the sequences of the 53 ribosomal proteins present in most of the archaeal species. This phylogeny was remarkably congruent with the rRNA phylogeny, suggesting that both reflected the actual phylogeny of the domain Archaea even if some nodes remained unresolved. In both cases, the branches leading to hyperthermophilic species were short, suggesting that the evolutionary rate of their genes has been slowed down by structural constraints related to environmental adaptation. In addition, to estimate the impact of lateral gene transfer (LGT) on our tree reconstruction, we used a new method that revealed that 8 genes out of the 53 ribosomal proteins used in our study were likely affected by LGT. This strongly suggested that a core of 45 nontransferred ribosomal protein genes existed in Archaea that can be tentatively used to infer the phylogeny of this domain. Interestingly, the tree obtained using only the eight ribosomal proteins likely affected by LGT was not very different from the consensus tree, indicating that LGT mainly brought random phylogenetic noise. The major difference involves organisms living in similar environments, suggesting that LGTs are mainly directed by the physical proximity of the organisms rather than by their phylogenetic proximity.     

Phylogenetic screening of a bacterial,metagenomic library using homing endonuclease restriction and marker insertion

Metagenomics provides access to the uncultured majority of the microbial world. The approaches employed in this field have, however, had limited success in linking functional genes to the taxonomic or phylogenetic origin of the organism they belong to. Here we present an efficient strategy to recover environmental DNA fragments that contain phylogenetic marker genes from metagenomic libraries. Our method involves the cleavage of 23S ribsosmal RNA (rRNA) genes within pooled library clones by the homing endonuclease I-CeuI followed by the insertion and selection of an antibiotic resistance cassette. This approach was applied to screen a library of 6500 fosmid clones derived from the microbial community associated with the sponge Cymbastela concentrica. Several fosmid clones were recovered after the screen and detailed phylogenetic and taxonomic assignment based on the rRNA gene showed that they belong to previously unknown organisms. In addition, compositional features of these fosmid clones were used to classify and taxonomically assign a dataset of environmental shotgun sequences. Our approach represents a valuable tool for the analysis of rapidly increasing, environmental DNA sequencing information.     

Nucleotide sequences of the 5S ribosomal RNA genes and their adjacent regions in Schizosaccharomyces pombe.

The organization of 5S ribosomal RNA (rRNA) genes in the genome of Schizosaccharomyces pombe has been investigated by restriction and hybridization analyses. The 5S rRNA genes were not linked to the other three species of rRNA genes which formed a repeating unit of 6.9 megadaltons, but located in other regions surrounded by heterogeneous sequences. The 5S rRNA gene organization in S. pombe is therefore different from those in other yeasts; Saccharomyces cerevisiae and Torulopsis utilis. Four restriction segments of different sizes each containing a single 5S rRNA gene were cloned on a bacterial plasmid, and the sequences in and around the RNA coding regions were determined. In the RNA coding regions, the sequences in four clones were identical with an exception that one residue has been substituted in one clone. In the flanking regions, the sequences were extremely rich in the AT-content and highly heterogeneous. The sequences were also markedly different from those in the corresponding regions of the other two yeasts. THe presence of T-clusters in the regions immediately after the RNA coding sequences was only notable homology among the four clones and the other two yeasts.     

Phylogenetic Screening of Ribosomal RNA Gene-Containing Clones in Bacterial Artificial Chromosome (BAC) Libraries from Different Depths in Monterey Bay

Marine picoplankton are central mediators of many oceanic biogeochemical processes, but much of their biology and ecology remains ill defined. One approach to better defining these environmentally significant microbes involves the acquisition of genomic data that can provide information about genome content, metabolic capabilities, and population variability in picoplankton assemblages. Previously, we constructed and phylogenetically screened a Bacterial Artificial Chromosome (BAC) library from surface water picoplankton of Monterey Bay. To further describe niche partitioning, metabolic variability, and population structure in coastal picoplankton populations, we constructed and compared several picoplankton BAC libraries recovered from different depths in Monterey Bay. To facilitate library screening, a rapid technique was developed (ITS-LH-PCR) to identify and quantify ribosomal RNA (rRNA) gene-containing BAC clones in BAC libraries. The approach exploited natural length variations in the internal transcribed spacer (ITS) located between SSU and LSU rRNA genes, as well as the presence and location of tRNA-alanine coding genes within the ITS. The correspondence between ITS-LH-PCR fragment sizes and 16S rRNA gene phylogenies facilitated rapid identification of rRNA genes in BAC clones without requiring direct DNA sequencing. Using this approach, 35 phylogenetic groups (previously identified by cultivation or PCR-based rRNA gene surveys) were detected and quantified among the BAC clones. Since the probability of recovering chimeric rRNA gene sequences in large insert BAC clones was low, we used these sequences to identify potentially chimeric sequences from previous PCR amplified clones deposited in public databases. Full-length SSU rRNA gene sequences from picoplankton BAC libraries, cultivated bacterioplankton, and nonchimeric RNA genes were then used to refine phylogenetic analyses of planktonic marine gamma Proteobacteria, Roseobacter, and Rhodospirillales species.(M.T. Suzuki and C.M. Preston) These authors contributed equally to this paper.     

Phylogenetic analysis of ribosomal RNA operons from uncultivated coastal marine bacterioplankton

Analyses of small subunit ribosomal RNA genes (SSU rDNAs) have significantly influenced our understanding of the composition of aquatic microbial assemblages. Unfortunately, SSU rDNA sequences often do not have sufficient resolving power to differentiate closely related species. To address this general problem for uncultivated bacterioplankton taxa, we analysed and compared sequences of polymerase chain reaction (PCR)-generated and bacterial artificial chromosome (BAC)-derived clones that contained most of the SSU rDNAs, the internal transcribed spacer (ITS) and the large subunit ribosomal RNA gene (LSU rDNA). The phylogenetic representation in the rRNA operon PCR library was similar to that reported previously in coastal bacterioplankton SSU rDNA libraries. We observed good concordance between the phylogenetic relationships among coastal bacterioplankton inferred from SSU or LSU rDNA sequences. ITS sequences confirmed the close intragroup relationships among members of the SAR11, SAR116 and SAR86 clades that were predicted by SSU and LSU rDNA sequence analyses. We also found strong support for homologous recombination between the ITS regions of operons from the SAR11 clade.     

Identification and characterization of metagenomic fragments from tidal flat sediment

Phylogenetic surveys based on cultivation-independent methods have revealed that tidal flat sediments are environments with extensive microbial diversity. Since most of prokaryotes in nature cannot be easily cultivated under general laboratory conditions, our knowledge on prokaryotic dwellers in tidal flat sediment is mainly based on the analysis of metagenomes. Microbial community analysis based on the 16S rRNA gene and other phylogenetic markers has been widely used to provide important information on the role of microorganisms, but it is basically an indirect means, compared with direct sequencing of metagenomic DNAs. In this study, we applied a sequence-based metagenomic approach to characterize uncultivated prokaryotes from tidal flat sediment. Two large-insert genomic libraries based on fosmid were constructed from tidal flat metagenomic DNA. A survey based on end-sequencing of selected fosmid clones resulted in the identification of clones containing 274 bacterial and 16 archaeal homologs in which majority were of proteobacterial origins. Two fosmid clones containing large metagenomic DNAs were completely sequenced using the shotgun method. Both DNA inserts contained more than 20 genes encoding putative proteins which implied their ecological roles in tidal flat sediment. Phylogenetic analyses of evolutionary conserved proteins indicate that these clones are not closely related to known prokaryotes whose genome sequence is known, and genes in tidal flat may be subjected to extensive lateral gene transfer, notably between domains Bacteria and Archaea. This is the first report demonstrating that direct sequencing of metagenomic gene library is useful in underpinning the genetic makeup and functional roles of prokaryotes in tidal flat sediments.     

Organisation of the S10, spc and alpha ribosomal protein gene clusters in prokaryotic genomes

Although it is well known that there is no long range colinearity in gene order in bacterial genomes, it is thought that there are several regions that are under strong structural constraints during evolution, in which gene order is extremely conserved. One such region is the str locus, containing the S10-spc-alpha operons. These operons contain genes coding for ribosomal proteins and for a number of housekeeping genes. We compared the organisation of these gene clusters in 111 sequenced prokaryotic genomes (99 bacterial and 12 archaeal genomes). We also compared the organisation to the phylogeny based on 16S ribosomal RNA gene sequences and the sequences of the ribosomal proteins L22, L16 and S14. Our data indicate that there is much variation in gene order and content in these gene clusters, both in bacterial as well as in archaeal genomes. Our data indicate that differential gene loss has occurred on multiple occasions during evolution. We also noted several discrepancies between phylogenetic trees based on 16S rRNA gene sequences and sequences of ribosomal proteins L16, L22 and S14, suggesting that horizontal gene transfer did play a significant role in the evolution of the S10-spc-alpha gene clusters.     

Phylogenetic and metagenomic analysis of Verrucomicrobia in former agricultural grassland soil

The bacterial phylum Verrucomicrobia has a widespread distribution, and is known to be one of the most common and diverse phyla in soil habitats. However, members of this phylum have typically been recalcitrant to cultivation methods, hampering the study of this presumably important bacterial group. In this study, we examine the phylogenetic diversity of the Verrucomicrobia in a former agricultural field and gain access to genomic information via a metagenomic approach. We examined Verrucomicrobia -like 16S rRNA gene sequences recovered from general bacterial and phylum-specific libraries, revealing a dominance of subdivisions 1 and 2. A PCR-based screening method was developed to identify inserts containing verrucomicrobial 16S rRNA genes within a large-insert metagenomic library, and on screening of 28 800 clones, four fosmids were identified as containing verrucomicrobial genomic DNA. Full-length sequencing of fosmid inserts and gene annotation identified a total of 98 ORFs, representing a range of functions. No conservation of gene order was observed adjacent to the ribosomal operons. Fosmid inserts were further analyzed for tetranucleotide frequencies to identify remnants of past horizontal gene transfer events. The metagenomic approach utilized proved to be suitable for the recovery of verrucomicrobial genomic DNA, thereby providing a window into the genomes of members of this important, yet poorly characterized, bacterial phylum.     

Phylogeny-function analysis of (meta)genomic libraries: screening for expression of ribosomal RNA genes by large-insert library fluorescent in situ hybridization (LIL-FISH)

We assessed the utility of fluorescent in situ hybridization (FISH) in the screening of clone libraries of (meta)genomic or environmental DNA for the presence and expression of bacterial ribosomal RNA (rRNA) genes. To establish proof-of-principle, we constructed a fosmid-based library in Escherichia coli of large-sized genomic DNA fragments of the mycophagous soil bacterium Collimonas fungivorans, and hybridized 768 library clones with the Collimonas-specific fluorescent probe CTE998-1015. Critical to the success of this approach (which we refer to as large-insert library FISH or LIL-FISH) was the ability to induce fosmid copy number, the exponential growth status of library clones in the FISH assay and the use of a simple pooling strategy to reduce the number of hybridizations. Twelve out of 768 E. coli clones were suspected to harbour and express Collimonas 16S rRNA genes based on their hybridization to CTE998-1015. This was confirmed by the finding that all 12 clones were also identified in an independent polymerase chain reaction-based screening of the same 768 clones using a primer set for the specific detection of Collimonas 16S ribosomal DNA (rDNA). Fosmids isolated from these clones were grouped by restriction analysis into two distinct contigs, confirming that C. fungivorans harbours at least two 16S rRNA genes. For one contig, representing 1-2% of the genome, the nucleotide sequence was determined, providing us with a narrow but informative view of Collimonas genome structure and content.     

PCR/DGGE and 16S rRNA gene library analysis of the colonic microbiota of HLA-B27/beta2-microglobulin transgenic rats

AIMS: To determine the phylogenetic composition of the colonic microbiota of transgenic (TG) HLA-B27 rats using 16S ribosomal RNA (rRNA) gene sequences obtained from denaturing gradient gel electrophoresis (DGGE) gels and sequences from a 16S rRNA gene library. METHODS AND RESULTS: Colonic microbiota of TG and nontransgenic (NT) rats harboured by 10-week-old and 6-month-old animals was screened using PCR/DGGE. Six months old TG rats had marked inflammation of the colon compared with 10-week-old TG and NT rats. The DGGE profiles of rats with inflamed colon were similar from rat to rat (Dice's Similarity Coefficient proximal colon 73%, distal colon 83%) whereas profiles from animals without inflammation were dissimilar (52-64%). Identifications of bacterial origins of 16S rRNA gene sequences obtained from DGGE gels (200 bp) and from 16S rRNA clones (450 bp) of the colonic microbiota of diseased rats gave sequences most closely phylogenetically affiliated with uncultured or unknown bacteria. CONCLUSIONS: PCR/DGGE was shown to be an effective method to compare the colonic microbiota composition of TG and NT rats relative to the progression of inflammatory disease. Sequencing of 16S rRNA gene fragments from DGGE gels or 16S rRNA gene clones from a random library showed that uncultured or unknown bacteria were most commonly detected by both methods. It can be concluded that it would be better in future studies to search for the antigens produced by the gut microbiota against which the dysfunctional immune system reacts rather than seek phylogenetic associations. SIGNIFICANCE AND IMPACT OF THE STUDY: PCR/DGGE can be used as a rapid initial screening method to compare the composition of bacterial communities of initially unknown composition that are associated with the development of intestinal disease.     

Metagenomic assessment of a sulfur-oxidizing enrichment culture derived from marine sediment

The biological oxidation of reduced sulfur compounds is a critically important process in global sulfur biogeochemistry. In this study, we enriched from marine sediments under denitrifying conditions, chemolithotrophic sulfur oxidizers that could oxidize a variety of reduced sulfur compounds: thiosulfate, tetrathionate, sulfide, and polysulfide. Two major phylotypes of 16S rRNA gene (>99% identity in each phylotype) were detected in this enrichment culture. In order to characterize sulfide oxidation, we sequenced and characterized one fosmid clone (43.6 kb) containing the group I sulfide-quinone reductase (sqr) gene. Interestingly, four putative rhodanese genes were found in this clone. Furthermore, comparative alignment with the closest genome of Thiomicrospira crunogena XCL2 revealed that three homologous genes were located within the vicinity of the sqr gene. Fosmid clones harboring carbon fixation (cbbL and cbbM) and denitrification (narG) genes were screened, and the phylogeny of the functional genes was analyzed. Along with the comparison between the sqr-containing fosmid clones and the relevant -proteobacteria, our phylogenetic study based on the 16S rRNA gene and carbon fixation genes suggest the prevalence of chemolithotrophic -proteobacteria in the denitrifying cultures. The findings of this study imply that a combination of cultivation and metagenomic approaches might provide us with a glimpse into the characteristics of sulfur oxidizers in marine sediments.     

Ribosomal and protein coding gene based multigene phylogeny on the family Streptomycetaceae

The phylogenetic relationship among the three genera of the family Streptomycetaceae was examined using the small and large subunit ribosomal RNA genes, and the gyrB, rpoB, trpB, atpD and recA genes. The total stretches of the analyzed ribosomal genes were 4.2kb, and those of five protein coding genes were 4.5 kb. The resultant phylogenetic trees confirmed that each genus formed an independent clade in the majority of cases. The G+C contents of rRNA genes were 56.9-58.9 mol%, and those of protein coding genes were 65.4-72.4 mol%, the latter being closer to those of the genomic DNAs. The average nucleotide sequence identity between the organisms were 94.1-96.4% for rRNA genes and 85.7-90.6% for protein coding genes, thus indicating that protein coding genes can give higher resolution than rRNA genes. In addition, the protein coding gene trees were more stable than the rRNA gene trees, supported by higher bootstrap values and other treeing algorithms. Moreover, the genome data of six Streptomyces species indicated that many protein coding genes exhibited higher correlations with genome relatedness. The combined gene sequences were also shown to give a better resolution with higher stability than any single genes, though not necessarily more correlated with genome relatedness. It is evident from this study that the rRNA gene based phylogeny can be misleading, and also that protein coding genes have a number of advantages over the rRNA genes as the phylogenetic markers including a high correlation with the genome relatedness.     

Evidence for existence of "mesotogas," members of the order Thermotogales adapted to low-temperature environments

All cultivated isolates of the bacterial order Thermotogales are either thermophiles or hyperthermophiles, but Thermotogales 16S rRNA gene sequences have been detected in many mesophilic anaerobic and microaerophilic environments, particularly within communities involved in the remediation of pollutants. Here we provide metagenomic evidence for the existence of Thermotogales lineages, which we informally call "mesotoga," that are adapted to growth at lower temperatures. Two fosmid clones containing mesotoga DNA, originating from a low-temperature enrichment culture that degrades a polychlorinated biphenyl congener, were sequenced. Phylogenetic analysis clearly puts this bacterial lineage within the Thermotogales order, with the rRNA gene trees and 21 of 58 open reading frames strongly supporting this relationship. An analysis of protein sequence composition showed that mesotoga proteins are adapted to function at lower temperatures than are their identifiable homologs from thermophilic and hyperthermophilic members of the order Thermotogales, supporting the notion that this bacterium lives and grows optimally at lower temperatures. The phylogenetic analysis suggests that the mesotoga lineage from which our fosmids derive has used both the acquisition of genes from its neighbors and the modification of existing thermophilic sequences to adapt to a mesophilic lifestyle.     

Widespread distribution of archaeal reverse gyrase in thermophilic bacteria suggests a complex history of vertical inheritance and lateral gene transfers

Reverse gyrase, an enzyme of uncertain funtion, is present in all hyperthermophilic archaea and bacteria. Previous phylogenetic studies have suggested that the gene for reverse gyrase has an archaeal origin and was transferred laterally (LGT) to the ancestors of the two bacterial hyperthermophilic phyla, Thermotogales and Aquificales. Here, we performed an in-depth analysis of the evolutionary history of reverse gyrase in light of genomic progress. We found genes coding for reverse gyrase in the genomes of several thermophilic bacteria that belong to phyla other than Aquificales and Thermotogales. Several of these bacteria are not, strictly speaking, hyperthermophiles because their reported optimal growth temperatures are below 80 degrees C. Furthermore, we detected a reverse gyrase gene in the sequence of the large plasmid of Thermus thermophilus strain HB8, suggesting a possible mechanism of transfer to the T. thermophilus strain HB8 involving plasmids and transposases. The archaeal part of the reverse gyrase tree is congruent with recent phylogenies of the archaeal domain based on ribosomal proteins or RNA polymerase subunits. Although poorly resolved, the complete reverse gyrase phylogeny suggests an ancient acquisition of the gene by bacteria via one or two LGT events, followed by its secondary distribution by LGT within bacteria. Finally, several genes of archaeal origin located in proximity to the reverse gyrase gene in bacterial genomes have bacterial homologues mostly in thermophiles or hyperthermophiles, raising the possibility that they were co-transferred with the reverse gyrase gene. Our new analysis of the reverse gyrase history strengthens the hypothesis that the acquisition of reverse gyrase may have been a crucial evolutionary step in the adaptation of bacteria to high-temperature environments. However, it also questions the role of this enzyme in thermophilic bacteria and the selective advantage its presence could provide.     

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.     

Evaluation of ribosomal RNA gene restriction patterns for the classification of Bacillus species and related genera

AIMS: To identify Bacillus species and related genera by fingerprinting based on ribosomal RNA gene restriction patterns; to compare ribosomal RNA gene restriction patterns-based phylogenetic trees with trees based on 16S rRNA gene sequences; to evaluate the usefulness of ribosomal RNA gene restriction patterns as a taxonomic tool for the classification of Bacillus species and related genera. METHODS AND RESULTS: Seventy-eight bacterial species which include 42 Bacillus species, 31 species from five newly created Bacillus-related genera, and five species from five phenotypically related genera were tested. A total of 77 distinct 16S rRNA gene hybridization banding patterns were obtained. The dendrogram resulting from UPGMA analysis showed three distinct main genetic clusters at the 75% banding pattern similarity. A total of 77 distinct 23S and 5S rRNA genes hybridization banding patterns were obtained, and the dendrogram showed four distinct genetic clusters at the 75% banding pattern similarity. A third dendrogram was constructed using a combination of the data from the 16S rRNA gene fingerprinting and the 23S and 5S rRNA genes fingerprinting. It revealed three distinct main phylogenetic clusters at the 75% banding pattern similarity. CONCLUSIONS: The Bacillus species along with the species from related genera were identified successfully and differentiated by ribosomal RNA gene restriction patterns, and most were distributed with no apparent order in various clusters on each of the three dendrograms. SIGNIFICANCE AND IMPACT OF THE STUDY: Our data indicate that ribosomal RNA gene restriction patterns can be used to reconstruct the phylogeny of the Bacillus species and derived-genera that approximates, but does not duplicate, phylogenies based on 16S rRNA gene sequences.