Scratching the surface of the rare biosphere with ribosomal sequence tag primers

Increasingly large datasets of 16S rRNA gene sequences reveal new information about the extent of microbial diversity and the surprising extent of the rare biosphere. Currently, many of the largest datasets are represented by short and variable ribosomal sequence tags (RSTs) that are limited in their ability to accurately assign sequences to broad-scale phylogenetic trees. In this study, we selected 30 rare RSTs from existing sequence datasets and designed primers to amplify c. 1400 bases of the 16S rRNA gene to determine whether these sequences were represented by existing databases or if they might reveal new lineages within the Bacteria. Approximately one-third of the RST primers successfully amplified longer portions of these low-abundance 16S rRNA genes in a specific manner. Subsequent phylogenetic analysis demonstrated that most of these sequences were (1) distantly related to existing cultivated microorganisms and (2) closely related to uncultivated clone sequences that were recently deposited in GenBank. The presence of so many recently collected 16S rRNA gene reference sequences in existing databases suggests that progress is being made quickly towards a microbial census, one which has begun scratching the surface of the 'rare biosphere'.     

Detection of Desulfotomaculum in an Italian rice paddy soil by 16S ribosomal nucleic acid analyses

Two specific primers were developed for the amplification of 16S rRNA genes of Desulfotomaculum lineage 1 to detect members of the genus Desulfotomaculum in rice field soil. The combination of both primers in PCR allowed the specific amplification and cloning of ten 16S rDNA sequences of this group from rice paddy soil DNA extracts. The phylogenetic analysis showed that these sequences formed a deeply branching cluster within Desulfotomaculum lineage 1, together with two sequences from the database and two sequences from a hydrocarbon-contaminated aquifer. Dissimilarity values to validly described species, including recently isolated strains of Desulfotomaculum from rice paddy microcosms, were higher than 12%. Within the new cluster the cloned sequences formed three separate groups which were each represented by at least two sequences with identities of >/=99% while one sequence represented an additional group. The sequences should represent sulfate-reducing organisms because they clearly fell into the physiologically coherent group of Gram-positive sulfate reducers. The relative abundance of bacteria of the Desulfotomaculum lineage 1 in rice paddy soil and root samples was estimated with rRNA dot blot hybridizations of extracted RNA. The relative RNA content of Desulfotomaculum lineage 1 was 0.55% in the bulk soil and 1% in the rice root samples, respectively, of the total 16S rRNA content (probe Eub338). Hybridization of rRNA with a probe targeting the new cluster represented by the cloned sequences confirmed the high abundance of 16S rRNA sequences from this cluster in the rice paddy field samples. Another hybridization probe detecting Desulfotomaculum acetoxidans and two closely related Desulfotomaculum isolates from rice paddy soil indicated that these bacteria were less abundant.     

Sequence Heterogeneity between the Two Genes Encoding 16s Rrna from the Halophilic Archaebacterium Haloarcula Marismortui

The halophilic archaebacterium, Haloarcula marismortui, contains two nonadjacent ribosomal RNA operons, designated rrnA and rrnB, in its genome. The 16S rRNA genes within these operons are 1472 nucleotides in length and differ by nucleotide substitutions at 74 positions. The substitutions are not uniformly distributed but rather are localized within three domains of 16S rRNA; more than two-thirds of the differences occur within the domain bounded by nucleotides 508 and 823. This domain is known to be important for P site binding of aminoacylated tRNA and for 30-50S subunit association. Using S1 nuclease protection, it has been shown that the 16S rRNAs transcribed from both operons are equally represented in the functional 70S ribosome population. Comparison of these two H. marismortui sequences to the 16S gene sequences from related halophilic genera suggests that (i) in diverging genera, mutational differences in 16S gene sequences are not clustered but rather are more generally distributed throughout the length of the 16S sequence, and (ii) the rrnB sequence, particularly within the 508-823 domain, is more different from the out group sequences than is the rrnA sequence. Several possible explanations for the evolutionary origin and maintenance of this sequence heterogeneity within 16S rRNA of H. marismortui are discussed.     

Species diversity of magnetotactic bacteria from the Ol’khovka River, Russia

A fraction of magnetotactic bacteria was isolated by magnetic separation from the water and silt samples collected from the Ol’khovka River (Kislovodsk, Russia). A 16S rRNA clone library was obtained from the total DNA of the fraction by PCR amplification and molecular cloning. Phylogenetic analysis of 67 16S rRNA gene sequences of randomly selected clones demonstrated that two phylotypes of magnetotactic bacteria were present in the library: the first phylotype consisted of 42 sequences and the second one included only one sequence. The remaining 24 sequences belonged to non-magnetotactic bacteria. According to the results of phylogenetic analysis, both phylotypes were magnetotactic cocci; the predominant sequences were almost identical to the 16S rRNA sequence of the freshwater coccus TB24 (X81185.1) identified earlier among the magnetotactic bacteria isolated from Lake Chiemsee (Bavaria). The phylotype represented by a single sequence formed a separate branch in the dendrogram, with 97% similarity between its sequence and that of TB24. The discovered phylotypes formed with the sequences of uncultured freshwater magnetotactic cocci a separate branch within the class Alphaproteobacteria and presumably belonged to a separate family within the recently described order Magnetococcales. Despite the fact that phylogenetic analysis of the 16S rRNA clone library did not reveal any phylotypes of magnetotactic spirilla, after the secondary enrichment of the fraction of magnetotactic bacteria using the “race track” technique, a new strain of magnetotactic spirilla, Magnetospirillum SO-1, was isolated. The closest relative of strain SO-1 was the previously described magnetotactic spirillum Magnetospirillum magneticum AMB-1.     

Legionella species diversity in an acidic biofilm community in Yellowstone National Park

Legionella species are frequently detected in aquatic environments, but their occurrence in extreme, acidic, geothermal habitats has not been explored with cultivation-independent methods. We investigated a predominately eukaryotic algal mat community in a pH 2.7 geothermal stream in Yellowstone National Park for the presence of Legionella and potential host amoebae. Our analyses, using PCR amplification with Legionella-specific primers targeting 16S rRNA genes, detected four known Legionella species, as well as Legionella sequences from species that are not represented in sequence databases, in mat samples and cultivated isolates. The nonrandom occurrence of sequences detected at lower (30 degrees C) and higher (35 to 38 degrees C) temperatures suggests that natural thermal gradients in the stream influence Legionella species distributions in this mat community. We detected only one sequence, Legionella micdadei, from cultivated isolates. We cultured and sequenced partial 18S rRNA gene regions from two potential hosts, Acanthamoeba and Euglena species.     

Use of Subtractive Hybridization To Design Habitat-Based Oligonucleotide Probes for Investigation of Natural Bacterial Communities

We describe a rapid oligonucleotide probe design strategy based on subtractive hybridization which yields probes for 16S rRNA or rRNA genes of individual members of microbial communities that are specific within the context of those communities. This strategy circumvents the need to sequence many similar or identical clones of dominant members of a community. Radioactively labeled subfragments of a cloned 16S rRNA gene sequence for which a probe is required (target) were hybridized with biotinylated total 16S ribosomal DNA (rDNA) amplified from the microbial community, and the hybrids formed were subsequently discarded. The remaining enriched fragments were used to screen a library consisting of cloned subfragments of the target sequence by colony hybridization in order to identify the variable regions of the 16S rRNA gene with the required specificity. The sequencing of random clones in one 16S rDNA library demonstrated that only those clones with 100% sequence identity with the probe fragment were detected by it. Moreover, sequencing of other, randomly selected, probe-positive clones revealed 100% sequence identity with the probe. Probes developed in this way tended to correspond to more variable regions of the 16S rRNA if the target sequences were similar to the sequences of other clones in the library and to less variable regions if the target sequences were phylogenetically isolated within the clone library. Although the absolute specificity of the latter probes, as assessed by comparison with available database sequences, was lower than the absolute specificity of the probes from the more variable regions, they were specific within the context of the environmental samples from which they were derived.     

Cloning and study of inserted sequences of gene 28S in Drosophila melanogaster ribosomal RNA

Cloning of fragments of ribosomal genes containing insertions in the 28S RNA gene has been reported earlier. Subcloning of DNA fragments corresponding to insertion sequences and their hybridization with DNA, RNA and polytene chromosomes from different flies is described. Type 1 insertions (containing BamI sites) are highly heterogeneous in length and sequence even in homozygotes. Type 2 insertions (with EcoRI sites) are rather homogeneous. Two types of insertions are represented in the D. melanogaster genome by 50 and 30 copies, respectively. Restriction fragments with insertions significantly differ in DNA from embryos and larvae. D. simulans and D. virilis also contain the sequences of both types of insertions, though in fewer number of copies. Type 1 insertions seem to be poorly transcribed, and type 2 insertions are not transcribed at all. Among 2000 recombinant clones screened a number of DI plasmids hybridizing to isolated insertions were obtained. Six of them were mapped with restriction endonucleases and hybridized with insertion fragments. rRNA and polytene chromosomes. All of these DI plasmids hybridize with the nucleoli, one with the chromocenter and one with the 79F 3L site. In LI9, not coding for rRNA, the sequences, corresponding to two types on insertions are located only a few kilobases apart. D17a does not encode for rRNA, but hybridizes in situ only with the nucleoli.     

Molecular analysis of a haplosporidian parasite from cultured New Zealand abalone Haliotis iris

In the Austral summer and autumn of 2000 and 2001, mortalities of black-footed abalone Haliotis iris (Martyn, 1784) occurred in a commercial facility in New Zealand. Histological analyses suggested that infection by a haplosporidian parasite was responsible. To confirm identification as a haplosporidian and to help determine if this parasite represented a new, undescribed species, DNA was extracted from infected host tissues scored as positive for infection by histological examination. Small-subunit rRNA (SSU rRNA) gene sequences from both the host abalone and a parasitic organism were amplified by PCR and characterized. Although the sequence for this parasite was novel, not matching any known SSU rRNA gene sequences, phylogenetic analyses strongly supported grouping this parasite with the haplosporidians. Parsimony analyses placed the parasite at the base of the phylum Haplosporidia, ancestral to Urosporidium crescens and the Haplosporidium, Bonamia, and Minchinia species. Sequencing of multiple parasite DNA clones revealed a single polymorphic site in the haplosporidian SSU rRNA gene sequence.     

Sequence ordinations: a multivariate analysis approach to analysing large sequence data sets

Ordination is a powerful method for analysing complex data setsbut has been largely ignored in sequence analysis. This papershows how to use principal coordinates analysis to find low–dimensionalrepresentations of distance matrices derived from aligned setsof sequences. The method takes a matrix of Euclidean distancesbetween all pairs of sequence and finds a coordinate space wherethe distances are exactly preserved The main problem is to finda measure of distance between aligned sequences that is Euclidean.The simplest distance function is the square root of the percentagedifference (as measured by identities) between two sequences,where one ignores any positions in the alignment where thereis a gap in any sequence. If one does not ignore positions witha gap, the distances cannot be guaranteed to be Euclidean butthe deleterious effects are trivial. Two examples of using themethod are shown. A set of 226 aligned globins were analysedand the resulting ordination very successfully represents theknown patterns of relationship between the sequences. In theother example, a set of 610 aligned 5S rRNA sequences were analysed.Sequence ordinations complement phylogenetic analyses. Theyshould not be viewed as a complete alternative.     

The nucleotide sequences of 5.8-S ribosomal RNA from Xenopus laevis and Xenopus borealis

1. The nucleotide sequence of 5.8-S rRNA from Xenopus laevis is given; it differs by a C in equilibrium U transition at position 140 from the 5.8-S rRNA of Xenopus borealis. 2. The sequence contains two completely modified and two partially modified residues. 3. Three different 5' nucleotides are found: pU-C-G (0.4) pC-G (0.2) and pG (0.4). 4. The 3' terminus is C not U as in all other 5.8-S sequences so far determined. 5. The X. laevis sequence differs from the mammalian and turtle sequences by five and six residue changes respectively. 6. A ribonuclease-resistant hairpin loop is a principle feature of secondary structure models proposed for this molecule. 7. Sequence heterogeneity may occur at one position at a very low level (approximately 0.01) in X. laevis 5.8-S rRNA, while none was detected in X. borealis or HeLa cell 5.8-S rRNA.     

A comparison of the small ribosomal RNA genes from the mitochondrial DNA of the great apes and humans: sequence, structure, evolution, and phylogenetic implications

Restriction endonuclease fragments produced by EcoRI/AvaI or KpnI digestion and containing the small (12S) ribosomal RNA (rRNA) genes from the mitochondrial DNAs (mtDNAs) of the common chimpanzee, pygmy chimpanzee, gorilla, and orangutan were inserted into the plasmids pBR322 or pADD1. After species verification the inserted fragments were digested with SauIIIA, subcloned into M13mp7 vectors, and sequenced. The small rRNA gene sequences were compared with each other and with the published human sequence (Anderson et al. 1981). Substitutions were detected at 118 of the 955 nucleotide positions compared. Pairwise, the sequence differences ranged from 1% (between the chimpanzee species) to 9% (comparisons involving the orangutan); the proportion that were transitions ranged from 87% to 100%. Deletions and/or additions were noted at seven locations. With respect to evolutionary sequence lability, kinetic analysis indicated the presence of at least two classes of nucleotide positions; the more labile class occurs in sequences thought to form self-complementary duplexes (stems) in the mature rRNA. The high frequency of compensating substitutions, which maintain base-pairing within these sequences, corroborates their inferred structure. Phylogenetic inferences drawn from the sequence comparisons support the notion of an approximately equidistant relationship among chimpanzees, gorilla, and man, with the orangutan much less closely related. However, inference from a shared deletion suggests that the gorilla and the chimpanzees may be more closely related to one another than they are to man.     

Common 5S rRNA Variants Are Likely To Be Accepted in Many Sequence Contexts

Over evolutionary time RNA sequences which are successfully fixed in a population are selected from among those that satisfy the structural and chemical requirements imposed by the function of the RNA. These sequences together comprise the structure space of the RNA. In principle, a comprehensive understanding of RNA structure and function would make it possible to enumerate which specific RNA sequences belong to a particular structure space and which do not. We are using bacterial 5S rRNA as a model system to attempt to identify principles that can be used to predict which sequences do or do not belong to the 5S rRNA structure space. One promising idea is the very intuitive notion that frequently seen sequence changes in an aligned data set of naturally occurring 5S rRNAs would be widely accepted in many other 5S rRNA sequence contexts. To test this hypothesis, we first developed well-defined operational definitions for a Vibrio region of the 5S rRNA structure space and what is meant by a highly variable position. Fourteen sequence variants (10 point changes and 4 base-pair changes) were identified in this way, which, by the hypothesis, would be expected to incorporate successfully in any of the known sequences in the Vibrio region. All 14 of these changes were constructed and separately introduced into the Vibrio proteolyticus 5S rRNA sequence where they are not normally found. Each variant was evaluated for its ability to function as a valid 5S rRNA in an E. coli cellular context. It was found that 93% (13/14) of the variants tested are likely valid 5S rRNAs in this context. In addition, seven variants were constructed that, although present in the Vibrio region, did not meet the stringent criteria for a highly variable position. In this case, 86% (6/7) are likely valid. As a control we also examined seven variants that are seldom or never seen in the Vibrio region of 5S rRNA sequence space. In this case only two of seven were found to be potentially valid. The results demonstrate that changes that occur multiple times in a local region of RNA sequence space in fact usually will be accepted in any sequence context in that same local region.     

On molecular taxonomy: what is in a name?

Gene sequences of small portions of the genome are often used for premature detailed taxonomic changes, neglecting polyphasic taxonomy, which should also consider phenotypical characteristics. Three examples are given: (i) Recently, members of the genera Eperythrozoon and Haemobartonella have been moved, correctly so, from the Rickettsiales to the Mycoplasmatales, but were assigned to the genus Mycoplasma, mostly on the basis of 16S rRNA sequence analysis. Not only is the 16S rRNA sequence similarity between 'classical' Mycoplasma and these species of Eperythrozoon and Haemobartonella less than that between some other well-recognised bacterial genera, but their biological differences amply justify their classification in different genera of the Mycoplasmatales. Furthermore, the move creates considerable confusion, as it necessitates new names for some species, with more confusion likely to come when the 16S rRNA sequences of the type species of Eperythrozoon, a name which has priority over Mycoplasma, will be analysed. (ii) In the Rickettsiales, members of the genera Anaplasma, Ehrlichia, Cowdria, Neorickettsia and Wolhbachia are so closely related phylogenetically on the basis of 16S rRNA sequences, and for some also of groESL operon sequences, that they have recently been fused, correctly so, into one family, the Anaplasmataceae, while the tribes Ehrlichieae and Wolbachieae have been abolished. Sequence diversity within the 'classical' genus Ehrlichia has led to classifying E. phagocytophila (including E. equi and the agent of human granulocytic ehrlichiosis), E. platys and E. bovis in the genus Anaplasma, while others have been retained in Ehrlichia, which also includes Cowdria ruminantium. E. sennetsu and E. risticii have been transferred to the genus Neorickettsia. 16S rRNA and GroEL sequences of 'classical' Anaplasma and some members of 'classical' Ehrlichia do show a close relationship, but differences in citrate synthase gene sequences, the GC content of this gene, and sequences of the gene encoding the beta-subunit of RNA polymerase, not to speak of the phenotypical differences, do not justify the fusion into one genus. Because of the phylogenetical diversity in Ehrlichia it is recommended that a new genus name be created for the E. phagocytophila genogroup (and E. platys and E. bovis). (iii) One of the conclusions of studies on the phylogeny of ticks of the subfamilies Rhipicephalinae and Hyalomminae, based on nucleotide sequences from 12S rRNA, cytochrome c oxidase I, the internal transcribed spacer 2, 18S rRNA, as well as morphological characters, is that Boophilus should be considered as a subgenus of Rhipicephalus. While Boophilus and Rhipicephalus are undoubtedly close, the obviously important morphological and biological differences between the genera Rhipicephalus and Boophilus are thus overruled by similarities in the sequences of a number of genes and this leads to considerable confusion. Polyphasic taxonomy amply justifies maintaining Boophilus as a separate genus, phylogenetically near to Rhipicephalus. This note is a plea for a cautious and balanced approach to taxonomy, taking into account molecular genotypical information, as far as is possible from different genes, as well as phenotypical characteristics.     

The first information on Schistocephalus nemachili (Cestoda: Ligulidae) from the Siberian groundling in the Lake Baikal region

Morphological data on Schistocephalus nemachili plerocercoid (a parasite of Nemachilis babatulus) and a nucleotide sequence of its 18S rRNA are presented. Samples were collected from the groundling inhabiting the Isumrudnoe lake (karst lake, belonging to the Lena River basin, 1200 m altitude). Length of plerocercoids 37-70 mm (mean 52.9 mm), width 1.6-7.3 (3.8), number of segments in srobile--229-348 (290). The sequence data on the 5'-terminal major domain of small rRNA of S. nemachili were obtained by mean of polymerase chain reaction and direct sequencing, they displayed 574 base positions. The nucleotide sequence was compared with known Plathelminthes sequences from the EMBL data base to compute data matrix under Juke and Cantor's model (1969). Distance matrix was used for obtaining a similarity dendrogram by Saitou and Nei's method (1987). Comparative analysis showed S. nemachili nucleotide sequence clustering with the cestode sequence (E. granulosus) within the dendrogram containing four trematodes, one cestode and one turbellarian sequences.     

The Evolutionary History of the Genus Acanthamoeba and the Identification of Eight New 18S rRNA Gene Sequence Types

ABSTRACT The 18S rRNA gene ( Rns ) phylogeny of Acanthamoeba is being investigated as a basis for improvements in the nomenclature and taxonomy of the genus. We previously analyzed Rns sequences from 18 isolates from morphological groups 2 and 3 and found that they fell into four distinct evolutionary lineages we called sequence types T1-T4. Here, we analyzed sequences from 53 isolates representing 16 species and including 35 new strains. Eight additional lineages (sequence types T5-T12) were identified. Four of the 12 sequence types included strains from more than one nominal species. Thus, sequence types could be equated with species in some cases or with complexes of closely related species in others. The largest complex, sequence type T4, which contained six closely related nominal species, included 24 of 25 keratitis isolates. Rns sequence variation was insufficient for full phylogenetic resolution of branching orders within this complex, but the mixing of species observed at terminal nodes confirmed that traditional classification of isolates has been inconsistent. One solution to this problem would be to equate sequence types and single species. Alternatively, additional molecular information will be required to reliably differentiate species within the complexes. Three sequence types of morphological group 1 species represented the earliest divergence in the history of the genus and, based on their genetic distinctiveness, are candidates for reclassification as one or more novel genera.     

How close is close: 16S rRNA sequence identity may not be sufficient to guarantee species identity.

16S rRNA (genes coding for rRNA) sequence comparisons were conducted with the following three psychrophilic strains: Bacillus globisporus W25T (T = type strain) and Bacillus psychrophilus W16AT, and W5. These strains exhibited more than 99.5% sequence identity and within experimental uncertainty could be regarded as identical. Their close taxonomic relationship was further documented by phenotypic similarities. In contrast, previously published DNA-DNA hybridization results have convincingly established that these strains do not belong to the same species if current standards are used. These results emphasize the important point that effective identity of 16S rRNA sequences is not necessarily a sufficient criterion to guarantee species identity. Thus, although 16S rRNA sequences can be used routinely to distinguish and establish relationships between genera and well-resolved species, very recently diverged species may not be recognizable.     

毛头鬼伞中一条新28S rRNA的发现及其同源性分析

本研究从担子菌毛头鬼伞(Coprinus comatus)菌丝中分离获得一条新的28S rRNA序列,序列长度为906bp(GenBank accession No.GU568178)。该序列是我们前期在从毛头鬼伞中克隆一种烟草花叶病毒(TMV)的抗性蛋白基因y3时意外获得的一条非目的条带。将此获得的序列通过NCBI的BLAST,以及与其同源序列进行Clustal w和MEGA聚类分析,证实该序列是28S rRNA,同时还发现毛头鬼伞的系统进化关系比较离散。此外,在这一新28S rRNA与TMV的抗性蛋白基因y3之间发现有两个同源区段有可能是PCR扩增y3基因时出现非目的条带的原因。在这两个同源区段中,其一区段与克隆y3基因时所用的PCR引物之一有较高的相似性,另一区段也是一般PCR引物的类似物。本研究中新28S rRNA序列的获得是PCR扩增中出现非目的条带的新例,该序列的发现及聚类分析的结果有助于真菌基因组学研究及真菌生物分子分类系统的建立。     

Evidence for genetic divergence in ribosomal RNA genes in mycobacteria

DNA was isolated from Mycobacterium phlei and from M. smegmatis. Each DNA sample was restricted with endonucleases, the fragments were separated by agarose gel electrophoresis and transferred to nitrocellulose film. Fragments of DNA containing rRNA sequences were identified by means of 125I-labelled rRNA of M. phlei or of M. smegmatis. The distributions of restriction endonuclease sites within the rRNA gene(s) and flanking sequences were found to be characteristic for each of the two species. Hybridizations with heterologous probes indicate that although M. phlei rRNA and M. smegmatis rRNA share regions of sequence homology, they are probably not identical in primary structure. The results suggest that the rRNA genes might prove to be useful taxonomic markers for mycobacteria.