Use of 16S rRNA and rpoB Genes as Molecular Markers for Microbial Ecology Studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such “heterogeneity hot spots” occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Use of 16S rRNA and rpoB genes as molecular markers for microbial ecology studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such "heterogeneity hot spots" occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

Intragenomic Heterogeneity of 16S rRNA Genes Causes Overestimation of Prokaryotic Diversity

Ever since Carl Woese introduced the use of 16S rRNA genes for determining the phylogenetic relationships of prokaryotes, this method has been regarded as the “gold standard” in both microbial phylogeny and ecology studies. However, intragenomic heterogeneity within 16S rRNA genes has been reported in many investigations and is believed to bias the estimation of prokaryotic diversity. In the current study, 2,013 completely sequenced genomes of bacteria and archaea were analyzed and intragenomic heterogeneity was found in 952 genomes (585 species), with 87.5% of the divergence detected being below the 1% level. In particular, some extremophiles (thermophiles and halophiles) were found to harbor highly divergent 16S rRNA genes. Overestimation caused by 16S rRNA gene intragenomic heterogeneity was evaluated at different levels using the full-length and partial 16S rRNA genes usually chosen as targets for pyrosequencing. The result indicates that, at the unique level, full-length 16S rRNA genes can produce an overestimation of as much as 123.7%, while at the 3% level, an overestimation of 12.9% for the V6 region may be introduced. Further analysis showed that intragenomic heterogeneity tends to concentrate in specific positions, with the V1 and V6 regions suffering the most intragenomic heterogeneity and the V4 and V5 regions suffering the least intragenomic heterogeneity in bacteria. This is the most up-to-date overview of the diversity of 16S rRNA genes within prokaryotic genomes. It not only provides general guidance on how much overestimation can be introduced when applying 16S rRNA gene-based methods, due to its intragenomic heterogeneity, but also recommends that, for bacteria, this overestimation be minimized using primers targeting the V4 and V5 regions.     

Igniting taxonomic curiosity: The amazing story of Amazonocrinis with the description of a new genus Ahomia gen. nov. and novel species of Ahomia,Amazonocrinis, and Dendronalium from the biodiversity-rich northeast region of India

Five cyanobacterial strains exhibiting Nostoc-like morphology were sampled from the biodiversity hotspots of the northeast region of India and characterized using a polyphasic approach. Molecular and phylogenetic analysis using the 16S rRNA gene indicated that the strains belonged to the genera Amazonocrinis and Dendronalium. In the present investigation, the 16S rRNA gene phylogeny clearly demarcated two separate clades of Amazonocrinis. The strain MEG8-PS clustered along with Amazonocrinis nigriterrae CENA67, which is the type strain of the genus. The other three strains ASM11-PS, RAN-4C-PS, and NP-KLS-5A-PS clustered in a different clade that was phylogenetically distinct from the Amazonocrinis sensu stricto clade. Interestingly, while the 16S rRNA gene phylogeny exhibited two separate clusters, the 16S–23S ITS region analysis did not provide strong support for the phylogenetic observation. Subsequent analyses raised questions regarding the resolving power of the 16S–23S ITS region at the genera level and the associated complexities in cyanobacterial taxonomy. Through this study, we describe a novel genus Ahomia to accommodate the members clustering outside the Amazonocrinis sensu stricto clade. In addition, we describe five novel species, Ahomia kamrupensis, Ahomia purpurea, Ahomia soli, Amazonocrinis meghalayensis, and Dendronalium spirale, in accordance with the International Code of Nomenclature for algae, fungi, and plants (ICN). Apart from further enriching the genera Amazonocrinis and Dendronalium, the current study helps to resolve the taxonomic complexities revolving around the genus Amazonocrinis and aims to attract researchers to the continued exploration of the tropical and subtropical cyanobacteria for interesting taxa and lineages.     

Unusual intragenomic and interspecific variability of <Emphasis Type="Italic">Geobacillus</Emphasis> 16S-23S rRNA internal transcribed spacers

The aim of this study was to evaluate the inter-and intraspecific as well as intragenomic variability of Geobacillus 16S–23S rRNA internal transcribed spacers without tRNA genes and to compare these sequences with sequences bearing tRNA genes. In this study the structural analysis was performed in a unique way because the length and the sequence of the structural blocks were adjusted to fit the structure of 16S–23S rRNA internal transcribed spacers of five different Geobacillus species. Our study demonstrated the mosaic-like structure of 16S–23S rRNA internal transcribed spacers in Geobacillus. Some characteristics of these spacers of geobacilli were not previously reported for other bacteria: unusually short conserved sequence in the 5′ end region, some identical conserved blocks in both 5′ and 3′ regions of 16S–23S rRNA internal transcribed spacers, the same sequence blocks in both 16S–23S and 23S–5S rRNA intergenic spacers. Our study demonstrated quite uniform arrangement of the sequence blocks in Geobacillus thermodenitrificans. This species diverged early in the phylogenetic tree of the genus Geobacillus. For the phylogenetically recent species Geobacillus kaustophilus and Geobacillus lituanicus the low inter-and intraspecific, but high intragenomic variability, as a consequence of recent phylogenetic events, was established.     

Phylogenetic analysis of cultivation‐resistant terrestrial cyanobacteria with massive sheaths (Stigonema spp. and Petalonema alatum,Nostocales, Cyanobacteria) using single‐cell and filament sequencing of environmental samples

Molecular assessment of a large portion of traditional cyanobacterial taxa has been hindered by the failure to isolate and grow them in culture. In this study, we developed an optimized protocol for single cell/filament isolation and 16S rRNA gene sequencing of terrestrial cyanobacteria with large mucilaginous sheaths, and applied it to determine the phylogenetic position of typical members of the genera Petalonema and Stigonema. A methodology based on a glass‐capillary isolation technique and a semi‐nested PCR protocol enabled reliable sequencing of the 16S rRNA gene from all samples analyzed. Ten samples covering seven species of Stigonema from Europe, North and Central America, and Hawaii, and the type species of Petalonema from Slovakia were sequenced. Contrary to some previous studies, which proposed a relationship with heteropolar nostocalean cyanobacteria, Petalonema appeared to belong to the family Scytonemataceae. Analysis of Stigonema specimens recovered a unique coherent phylogenetic cluster, substantially broadening our knowledge of the molecular diversity within this genus. Neither the uni‐ to biseriate species nor the multiseriate species formed monophyletic subclusters within the genus. Typical multiseriate species of Stigonema clustered in a phylogenetic branch derived from uni‐ to biseriate S. ocellatum Thuret ex Bornet & Flahault in our analysis, suggesting that species with more complex thalli may have evolved from the more simple ones. We propose the technique tested in this study as a promising tool for a future revision of the molecular taxonomy in cyanobacteria.     

Intra-Genomic Heterogeneity in 16S rRNA Genes in Strictly Anaerobic Clinical Isolates from Periodontal Abscesses

Background

Members of the genera Prevotella, Veillonella and Fusobacterium are the predominant culturable obligate anaerobic bacteria isolated from periodontal abscesses. When determining the cumulative number of clinical anaerobic isolates from periodontal abscesses, ambiguous or overlapping signals were frequently encountered in 16S rRNA gene sequencing chromatograms, resulting in ambiguous identifications. With the exception of the genus Veillonella, the high intra-chromosomal heterogeneity of rrs genes has not been reported.

Methods

The 16S rRNA genes of 138 clinical, strictly anaerobic isolates and one reference strain were directly sequenced, and the chromatograms were carefully examined. Gene cloning was performed for 22 typical isolates with doublet sequencing signals for the 16S rRNA genes, and four copies of the rrs-ITS genes of 9 Prevotella intermedia isolates were separately amplified by PCR, sequenced and compared. Five conserved housekeeping genes, hsp60, recA, dnaJ, gyrB1 and rpoB from 89 clinical isolates of Prevotella were also amplified by PCR and sequenced for identification and phylogenetic analysis along with 18 Prevotella reference strains.

Results

Heterogeneity of 16S rRNA genes was apparent in clinical, strictly anaerobic oral bacteria, particularly in the genera Prevotella and Veillonella. One hundred out of 138 anaerobic strains (72%) had intragenomic nucleotide polymorphisms (SNPs) in multiple locations, and 13 strains (9.4%) had intragenomic insertions or deletions in the 16S rRNA gene. In the genera Prevotella and Veillonella, 75% (67/89) and 100% (19/19) of the strains had SNPs in the 16S rRNA gene, respectively. Gene cloning and separate amplifications of four copies of the rrs-ITS genes confirmed that 2 to 4 heterogeneous 16S rRNA copies existed.

Conclusion

Sequence alignment of five housekeeping genes revealed that intra-species nucleotide similarities were very high in the genera Prevotella, ranging from 94.3–100%. However, the inter-species similarities were relatively low, ranging from 68.7–97.9%. The housekeeping genes rpoB and gyrB1 were demonstrated to be alternative classification markers to the species level based on intra- and inter-species comparisons, whereas based on phylogenetic tree rpoB proved to be reliable phylogenetic marker for the genus Prevotella.     

Phylogenetic Studies of <Emphasis Type="Italic">Gordonia</Emphasis> Species Based on <Emphasis Type="Italic">gyrB</Emphasis> and <Emphasis Type="Italic">secA1</Emphasis> Gene Analyses

Phylogenetic relations within the genus Gordonia were analyzed using partial gyrB and secA1 gene sequences of 23 type species in comparison with those of 16S rRNA gene. The gyrB and secA1 phylogenies showed agreement with that constructed using 16S rRNA gene sequences. The degrees of divergence of the gyrB and secA1 genes were approximately 3.4 and 1.7 times greater, respectively, than that of 16S rRNA gene. The gyrB gene showed more discriminatory power than either the secA1 or 16S rRNA gene, facilitating clear differentiation of any two Gordonia species using gyrB gene analysis. Our data indicate that gyrB and secA1 gene sequences are useful as markers for phylogenetic study and identification at the species level of the genus Gordonia.     

An Evaluation of Four Phylogenetic Markers in <Emphasis Type="Italic">Nostoc</Emphasis>: Implications for Cyanobacterial Phylogenetic Studies at the Intrageneric Level

The success of some phylogenetic markers in cyanobacteria owes to the design of cyanobacteria-specific primers, but a few studies have directly investigated the evolution “behavior” of the loci. In this study, we performed a case study in Nostoc to evaluate rpoC1, hetR, rbcLX, and 16S rRNA–tRNA^Ile–tRNA^Ala-23S rRNA internal transcribed spacer (ITS) as phylogenetic markers. The results indicated that the gene trees of these loci are not congruent with the phylogeny based on 16S rRNA gene. The mechanisms contributing to the incongruence include randomized variation and recombination. As the results suggested, one should be careful to choose the molecular markers for phylogenetic reconstruction at the intrageneric level in cyanobacteria.     

Synechococcus: 3 billion years of global dominance

Cyanobacteria are among the most important primary producers on the Earth. However, the evolutionary forces driving cyanobacterial species diversity remain largely enigmatic due to both their distinction from macro‐organisms and an undersampling of sequenced genomes. Thus, we present a new genome of a Synechococcus‐like cyanobacterium from a novel evolutionary lineage. Further, we analyse all existing 16S rRNA sequences and genomes of Synechococcus‐like cyanobacteria. Chronograms showed extremely polyphyletic relationships in Synechococcus, which has not been observed in any other cyanobacteria. Moreover, most Synechococcus lineages bifurcated after the Great Oxidation Event, including the most abundant marine picoplankton lineage. Quantification of horizontal gene transfer among 70 cyanobacterial genomes revealed significant differences among studied genomes. Horizontal gene transfer levels were not correlated with ecology, genome size or phenotype, but were correlated with the age of divergence. All findings were synthetized into a novel model of cyanobacterial evolution, characterized by serial convergence of the features, that is multicellularity and ecology.     

Reptodigitus Chapmanii (Nostocales,Hapalosiphonaceae) Gen. Nov.: A Unique Nostocalean (Cyanobacteria) Genus Based on a Polyphasic Approach1

The Nostocales is a monophyletic, heterocytous lineage of cyanobacteria capable of akinete production and division in multiple planes, depending upon family-level clade. While present in a variety of ecosystems, the diversity of the Nostocales has been poorly elucidated. Due to environmentally -induced phenotypic plasticity, morphology alone is often insufficient to determine the true phylogenetic placement of these taxa. In order to bridge this gap, taxonomists now employ the polyphasic approach, combining methods such as morphological analysis, phylogenetic analysis based on DNA sequence and genetic identity based on ribosomal genes, and secondary structure of the 16S-23S ITS and 16S rRNA gene sequences, as well as ecological characterization. Using this combined approach, a new genus and species (Reptodigitus chapmanii gen. et sp. nov.) isolated from the St. Johns River (Jacksonville, Florida, USA) within the Nostocales is herein described. Phylogenetic analyses place this taxon within the Hapalosiphonaceae, sister to the clade containing Fischerella, Hapalosiphon, and Westiellopsis. The 16S-23S ITS secondary folding structure analysis also supports the erection of this new genus.     

Detection of Seven Major Evolutionary Lineages in Cyanobacteria Based on the 16S rRNA Gene Sequence Analysis with New Sequences of Five Marine Synechococcus Strains

Although molecular phylogenetic studies of cyanobacteria on the basis of the 16S rRNA gene sequence have been reported, the topologies were unstable, especially in the inner branchings. Our analysis of 16S rRNA gene phylogeny by the maximum-likelihood and neighbor-joining methods combined with rate homogeneous and heterogeneous models revealed seven major evolutionary lineages of the cyanobacteria, including prochlorophycean organisms. These seven lineages are always stable on any combination of these methods and models, fundamentally corresponding to phylogenetic relationships based on other genes, e.g., psbA, rbcL, rnpB, rpoC, and tufA. Moreover, although known genotypic and phenotypic characters sometimes appear paralleled in independent lineages, many characters are not contradictory within each group. Therefore we propose seven evolutionary groups as a working hypothesis for successive taxonomic reconstruction. New 16S rRNA sequences of five unicellular cyanobacterial strains, PCC 7001, PCC 7003, PCC 73109, PCC 7117, and PCC 7335 of Synechococcus sp., were determined in this study. Although all these strains have been assigned to ``marine clusters B and C,' they were separated into three lineages. This suggests that the organisms classified in the genus Synechococcus evolved diversely and should be reclassified in several independent taxonomic units. Moreover, Synechococcus strains and filamentous cyanobacteria make a monophyletic group supported by a comparatively high statistical confidence value (80 to 100%) in each of the two independent lineages; therefore, these monophylies probably reflect the convergent evolution of a multicellular organization. Received: 3 September 1998 / Accepted: 30 November 1998     

POLYPHASIC CHARACTERIZATION OF DOLICHOSPERMUM SPP. AND SPHAEROSPERMOPSIS SPP. (NOSTOCALES,CYANOBACTERIA): MORPHOLOGY, 16S rRNA GENE SEQUENCES AND FATTY ACID AND SECONDARY METABOLITE PROFILES1

The genera Dolichospermum (Ralfs ex Bornet et Flahault) Wacklin, L. Hoffm. et Komárek and Sphaerospermopsis Zapomělová, Jezberová, Hrouzek, Hisem, K. ?eháková et Komárk.‐Legn. represent a highly diversified group of planktonic cyanobacteria that have been recently separated from the traditional genus Anabaena Bory ex Bornet et Flahault. In this study, morphological diversity, phylogeny of the 16S rRNA gene, production of fatty acids, and secondary metabolite profiles were evaluated in 33 strains of 14 morphospecies isolated from the Czech Republic. Clustering of the strains based on 16S rRNA gene sequences corresponded to wider groups of species in terms of morphology. The overall secondary metabolite and fatty acid profiles, however, were not correlated to each other and neither were they correlated to the 16S rRNA phylogeny nor the morphology of the strains. Nevertheless, a minor part of the detected secondary metabolites (19% of all compounds) was present only in close relatives and can be thus considered as autapomorphic features.     

Reconstruction of phylogenetic relationships within Grapsidae (Crustacea: Brachyura) and comparison of trans-isthmian versus amphi-atlantic gene flow based on mtDNA

Following taxonomic revisions in recent years, the originally large family Grapsidae MacLeay, 1838 has become a relatively small and morphologically homogeneous family in terms of adult and larval morphology. Most available molecular studies including more than one genus of the family have also suggested monophyly of the corresponding taxa. However, no single phylogenetic study has ever included all constituent genera of the Grapsidae. In the current study, a molecular phylogeny based on sequences of the mitochondrial 16S rRNA gene from all eight grapsid genera and 34 species is presented and suggests that up to four genera are not monophyletic. This is mainly due to the polyphyletic nature of the genus Pachygrapsus which can be found in six different lineages of the phylogeny, suggesting that the genus currently does not represent a single evolutionary lineage and is in need of taxonomic revision. Amphi-atlantic and trans-isthmian species pairs or populations in four genera are compared and reveal relatively constant and pronounced divergences across the Panama Isthmus as opposed to moderate divergences across the Atlantic Ocean, thereby suggesting occurrence of gene flow across the Atlantic Ocean during the past three million years.     

Evolution of gene orders in genomes of cyanobacteria

Genomes of 23 strains of cyanobacteria were comparatively analyzed using quantitative methods of estimation of gene order similarity. It has been found that reconstructions of phylogenesis of cyanobacteria based on the comparison of the orders of genes in chromosomes and nucleotide sequences appear to be similar. This confirms the applicability of quantitative measures of similarity of gene orders for phylogenetic reconstructions. In the evolution of marine unicellular planktonic cyanobacteria, genome rearrangements are fixed with a low rate (about 3% of gene order changes per 1% of 16S rRNA changes), whereas in other groups of cyanobacteria the gene order can change several times more rapidly. The gene orders in genomes of cyanobacteria and chloroplasts preserve a considerable degree of similarity. The closest relatives of chloroplasts among the analyzed cyanobacteria are likely to be strains from hot springs belonging to the genus Synechococcus. Comparative analysis of gene orders and nucleotide sequences strongly suggests that Synechococcus strains from different environments (sea, fresh waters, hot springs) are not related and belong to evolutionally distant lines.     

Culturable Actinobacteria from the Marine Sponge Hymeniacidon perleve: Isolation and Phylogenetic Diversity by 16S rRNA gene-RFLP Analysis

A total of 106 actinobacteria associated with the marine sponge Hymeniacidon perleve collected from the Yellow Sea, China were isolated using eight different media. The number of species and genera of actinobacteria recovered from the different media varied significantly, underlining the importance of optimizing the isolation conditions. The phylogenetic diversity of the actinobacteria isolates was assessed using 16S rRNA gene amplification–restriction fragment length polymorphism (RFLP) analysis of the 106 strains with different morphologies. The RFLP fingerprinting of selected strains by HhaI-digestion of the 16S rRNA genes resulted in 11 different patterns. The HhaI-RFLP analysis gave good resolution for the identification of the actinobacteria isolates at the genus level. A phylogenetic analysis using 16S rRNA gene sequences revealed that the isolates belonged to seven genera of culturable actinobacteria including Actinoalloteichus, Micromonospora, Nocardia, Nocardiopsis, Pseudonocardia, Rhodococcus, and Streptomyces. The dominant genus was Streptomyces, which represented 74% of the isolates. Three of the strains identified are candidates for new species.     

Gene copy number variation and its significance in cyanobacterial phylogeny

ABSTRACT: BACKGROUND: In eukaryotes, variation in gene copy numbers is often associated with deleterious effects, but may also have positive effects. For prokaryotes, studies on gene copy number variation are rare. Previous studies have suggested that high numbers of rRNA gene copies can be advantageous in environments with changing resource availability, but further association of gene copies and phenotypic traits are not documented. We used one of the morphologically most diverse prokaryotic phyla to test whether numbers of gene copies are associated with levels of cell differentiation. RESULTS: We implemented a search algorithm that identified 44 genes with highly conserved copies across 22 fully sequenced cyanobacterial taxa. For two very basal cyanobacterial species, Gloeobacter violaceus and a thermophilic Synechococcus species, distinct phylogenetic positions previously found were supported by identical protein coding gene copy numbers. Furthermore, we found that increased ribosomal gene copy numbers showed a strong correlation to cyanobacteria capable of terminal cell differentiation. Additionally, we detected extremely low variation of 16S rRNA sequence copies within the cyanobacteria. We compared our results for 16S rRNA to three other eubacterial phyla (Chroroflexi, Spirochaetes and Bacteroidetes). Based on Bayesian phylogenetic inference and the comparisons of genetic istances, we could confirm that cyanobacterial 16S rRNA paralogs and orthologs show significantly stronger conservation than found in other eubacterial phyla. Conclusions: A higher number of ribosomal operons could potentially provide an advantage to terminally differentiated cyanobacteria. Furthermore, we suggest that 16S rRNA gene copies in cyanobacteria are homogenized by both concerted evolution and purifying selection. In addition, the small ribosomal subunit in cyanobacteria appears to evolve at extraordinary slow evolutionary rates, an observation that has been made previously for morphological characteristics of cyanobacteria.     

Benthic cyanobacterial bloom impacts the reefs of South Florida (Broward County, USA)

Benthic cyanobacteria of the genus Lyngbya can form prominent mats and blooms in tropical and subtropical coral reef and seagrass habitats worldwide. A Lyngbya bloom on the reef tract offshore of Broward County, Florida, was first noted in 2002, and although it is seasonally variable in its distribution and abundance, it has persisted and spread over the past 3 years. In this study, the most abundant species of Lyngbya found in the blooms have been identified and compared to other species of Lyngbya by morphological and molecular methods. The most common species of Lyngbya is consistent with the properties of Lyngbya confervoides C. Agardh. The 16S ribosomal DNA sequence shares 88–92% identity with other known Lyngbya sequences, suggesting that this bloom consists primarily of a new, previously unsequenced species of Lyngbya. The second most common Lyngbya in the bloom is consistent with Lyngbya polychroa. This persistent bloom is a concern because it smothers octocorals and other invertebrates and negatively impacts these southeastern Florida reefs.