Development of an Analysis Pipeline Characterizing Multiple Hypervariable Regions of 16S rRNA Using Mock Samples

Objectives

There is much speculation on which hypervariable region provides the highest bacterial specificity in 16S rRNA sequencing. The optimum solution to prevent bias and to obtain a comprehensive view of complex bacterial communities would be to sequence the entire 16S rRNA gene; however, this is not possible with second generation standard library design and short-read next-generation sequencing technology.

Methods

This paper examines a new process using seven hypervariable or V regions of the 16S rRNA (six amplicons: V2, V3, V4, V6-7, V8, and V9) processed simultaneously on the Ion Torrent Personal Genome Machine (Life Technologies, Grand Island, NY). Four mock samples were amplified using the 16S Ion Metagenomics Kit^™ (Life Technologies) and their sequencing data is subjected to a novel analytical pipeline.

Results

Results are presented at family and genus level. The Kullback-Leibler divergence (D_KL), a measure of the departure of the computed from the nominal bacterial distribution in the mock samples, was used to infer which region performed best at the family and genus levels. Three different hypervariable regions, V2, V4, and V6-7, produced the lowest divergence compared to the known mock sample. The V9 region gave the highest (worst) average D_KL while the V4 gave the lowest (best) average D_KL. In addition to having a high D_KL, the V9 region in both the forward and reverse directions performed the worst finding only 17% and 53% of the known family level and 12% and 47% of the genus level bacteria, while results from the forward and reverse V4 region identified all 17 family level bacteria.

Conclusions

The results of our analysis have shown that our sequencing methods using 6 hypervariable regions of the 16S rRNA and subsequent analysis is valid. This method also allowed for the assessment of how well each of the variable regions might perform simultaneously. Our findings will provide the basis for future work intended to assess microbial abundance at different time points throughout a clinical protocol.     

Taxonomic Precision of Different Hypervariable Regions of 16S rRNA Gene and Annotation Methods for Functional Bacterial Groups in Biological Wastewater Treatment

High throughput sequencing of 16S rRNA gene leads us into a deeper understanding on bacterial diversity for complex environmental samples, but introduces blurring due to the relatively low taxonomic capability of short read. For wastewater treatment plant, only those functional bacterial genera categorized as nutrient remediators, bulk/foaming species, and potential pathogens are significant to biological wastewater treatment and environmental impacts. Precise taxonomic assignment of these bacteria at least at genus level is important for microbial ecological research and routine wastewater treatment monitoring. Therefore, the focus of this study was to evaluate the taxonomic precisions of different ribosomal RNA (rRNA) gene hypervariable regions generated from a mix activated sludge sample. In addition, three commonly used classification methods including RDP Classifier, BLAST-based best-hit annotation, and the lowest common ancestor annotation by MEGAN were evaluated by comparing their consistency. Under an unsupervised way, analysis of consistency among different classification methods suggests there are no hypervariable regions with good taxonomic coverage for all genera. Taxonomic assignment based on certain regions of the 16S rRNA genes, e.g. the V1&V2 regions – provide fairly consistent taxonomic assignment for a relatively wide range of genera. Hence, it is recommended to use these regions for studying functional groups in activated sludge. Moreover, the inconsistency among methods also demonstrated that a specific method might not be suitable for identification of some bacterial genera using certain 16S rRNA gene regions. As a general rule, drawing conclusions based only on one sequencing region and one classification method should be avoided due to the potential false negative results.     

High-Resolution Taxonomic Profiling of the Subgingival Microbiome for Biomarker Discovery and Periodontitis Diagnosis

The oral microbiome plays a key role for caries, periodontitis, and systemic diseases. A method for rapid, high-resolution, robust taxonomic profiling of subgingival bacterial communities for early detection of periodontitis biomarkers would therefore be a useful tool for individualized medicine. Here, we used Illumina sequencing of the V1-V2 and V5-V6 hypervariable regions of the 16S rRNA gene. A sample stratification pipeline was developed in a pilot study of 19 individuals, 9 of whom had been diagnosed with chronic periodontitis. Five hundred twenty-three operational taxonomic units (OTUs) were obtained from the V1-V2 region and 432 from the V5-V6 region. Key periodontal pathogens like Porphyromonas gingivalis, Treponema denticola, and Tannerella forsythia could be identified at the species level with both primer sets. Principal coordinate analysis identified two outliers that were consistently independent of the hypervariable region and method of DNA extraction used. The linear discriminant analysis (LDA) effect size algorithm (LEfSe) identified 80 OTU-level biomarkers of periodontitis and 17 of health. Health- and periodontitis-related clusters of OTUs were identified using a connectivity analysis, and the results confirmed previous studies with several thousands of samples. A machine learning algorithm was developed which was trained on all but one sample and then predicted the diagnosis of the left-out sample (jackknife method). Using a combination of the 10 best biomarkers, 15 of 17 samples were correctly diagnosed. Training the algorithm on time-resolved community profiles might provide a highly sensitive tool to detect the onset of periodontitis.     

Evaluation of different partial 16S rRNA gene sequence regions for phylogenetic analysis of microbiomes

Operational taxonomic units (OTUs) are conventionally defined at a phylogenetic distance (0.03—species, 0.05—genus, 0.10—family) based on full-length 16S rRNA gene sequences. However, partial sequences (700 bp or shorter) have been used in most studies. This discord may affect analysis of diversity and species richness because sequence divergence is not distributed evenly along the 16S rRNA gene. In this study, we compared a set each of bacterial and archaeal 16S rRNA gene sequences of nearly full length with multiple sets of different partial 16S rRNA gene sequences derived therefrom (approximately 440-700 bp), at conventional and alternative distance levels. Our objective was to identify partial sequence region(s) and distance level(s) that allow more accurate phylogenetic analysis of partial 16S rRNA genes. Our results showed that no partial sequence region could estimate OTU richness or define OTUs as reliably as nearly full-length genes. However, the V1-V4 regions can provide more accurate estimates than others. For analysis of archaea, we recommend the V1-V3 and the V4-V7 regions and clustering of species-level OTUs at 0.03 and 0.02 distances, respectively. For analysis of bacteria, the V1-V3 and the V1-V4 regions should be targeted, with species-level OTUs being clustered at 0.04 distance in both cases.     

Genetic Identification of Members of the Genus Corynebacterium at Genus and Species Levels with 16S rDNA-Targeted Probes

16S rRNA gene-targeted probes were designed for the identification of corynebacteria at the genus and species levels. The genus-specific probe hybridized all clinically important members of the genus Corynebacterium and could distinguish them from other coryneform bacteria and phylogenetically related high G + C% gram-positive bacteria, including Actinomyces, Rhodococcus, Gordona, Nocardia, Streptomyces, Brevibacterium and Mycobacterium. The species-specific probes for C. jeikeium and C. diphtheriae could differentiate these two species from other members of this genus. The probes were used to select corynebacteria among gram-positive clinical isolates which had been tentatively identified as corynebacteria by biochemical tests. We screened 59 strains with the genus-specific probe; 51 strains hybridized to the genus-specific probe, 8 did not. Of the 51 strains that hybridized to the genus-specific probe, 1 hybridized to the C. diphtheriae species probe and 13 hybridized to the C. jeikeium species probe. The 8 strains that did not hybridize to the genus probe were further characterized by analyzing cell wall diaminopimelic acid and partial 16S rRNA sequencing. The results indicated that these strains were distributed in the genera Arthrobacter and Brevibacterium.     

Bacterial community structure and location in Stilton cheese

The microbial diversity occurring in Stilton cheese was evaluated by 16S ribosomal DNA analysis with PCR-denaturing gradient gel electrophoresis. DNA templates for PCR experiments were directly extracted from the cheese as well as bulk cells harvested from a variety of viable-count media. The variable V3 and V4-V5 regions of the 16S genes were analyzed. Closest relatives of Lactococcus lactis, Enterococcus faecalis, Lactobacillus plantarum, Lactobacillus curvatus, Leuconostoc mesenteroides, Staphylococcus equorum, and Staphylococcus sp. were identified by sequencing of the DGGE fragments. Fluorescently labeled oligonucleotide probes were developed to detect Lactococcus lactis, Lactobacillus plantarum, and Leuconostoc mesenteroides in fluorescence in situ hybridization (FISH) experiments, and their specificity for the species occurring in the community of Stilton cheese was checked in FISH experiments carried out with reference cultures. The combined use of these probes and the bacterial probe Eub338 in FISH experiments on Stilton cheese sections allowed the assessment of the spatial distribution of the different microbial species in the dairy matrix. Microbial colonies of bacteria showed a differential location in the different parts of the cheese examined: the core, the veins, and the crust. Lactococci were found in the internal part of the veins as mixed colonies and as single colonies within the core. Lactobacillus plantarum was detected only underneath the surface, while Leuconostoc microcolonies were homogeneously distributed in all parts observed. The combined molecular approach is shown to be useful to simultaneously describe the structure and location of the bacterial flora in cheese. The differential distribution of species found suggests specific ecological reasons for the establishment of sites of actual microbial growth in the cheese, with implications of significance in understanding the ecology of food systems and with the aim of achieving optimization of the fermentation technologies as well as preservation of traditional products.     

Target region selection is a critical determinant of community fingerprints generated by 16S pyrosequencing

Pyrosequencing of 16S rRNA genes allows for in-depth characterization of complex microbial communities. Although it is known that primer selection can influence the profile of a community generated by sequencing, the extent and severity of this bias on deep-sequencing methodologies is not well elucidated. We tested the hypothesis that the hypervariable region targeted for sequencing and primer degeneracy play important roles in influencing the composition of 16S pyrotag communities. Subgingival plaque from deep sites of current smokers with chronic periodontitis was analyzed using Sanger sequencing and pyrosequencing using 4 primer pairs. Greater numbers of species were detected by pyrosequencing than by Sanger sequencing. Rare taxa constituted nearly 6% of each pyrotag community and less than 1% of the Sanger sequencing community. However, the different target regions selected for pyrosequencing did not demonstrate a significant difference in the number of rare and abundant taxa detected. The genera Prevotella, Fusobacterium, Streptococcus, Granulicatella, Bacteroides, Porphyromonas and Treponema were abundant when the V1-V3 region was targeted, while Streptococcus, Treponema, Prevotella, Eubacterium, Porphyromonas, Campylobacter and Enterococcus predominated in the community generated by V4-V6 primers, and the most numerous genera in the V7-V9 community were Veillonella, Streptococcus, Eubacterium, Enterococcus, Treponema, Catonella and Selenomonas. Targeting the V4-V6 region failed to detect the genus Fusobacterium, while the taxa Selenomonas, TM7 and Mycoplasma were not detected by the V7-V9 primer pairs. The communities generated by degenerate and non-degenerate primers did not demonstrate significant differences. Averaging the community fingerprints generated by V1-V3 and V7-V9 primers provided results similar to Sanger sequencing, while allowing a significantly greater depth of coverage than is possible with Sanger sequencing. It is therefore important to use primers targeted to these two regions of the 16S rRNA gene in all deep-sequencing efforts to obtain representational characterization of complex microbial communities.     

Effect of primers hybridizing to different evolutionarily conserved regions of the small-subunit rRNA gene in PCR-based microbial community analyses and genetic profiling

Genetic profiling techniques of microbial communities based on PCR-amplified signature genes, such as denaturing gradient gel electrophoresis or single-strand-conformation polymorphism (SSCP) analysis, are normally done with PCR products of less than 500-bp. The most common target for diversity analysis, the small-subunit rRNA genes, however, are larger, and thus, only partial sequences can be analyzed. Here, we compared the results obtained by PCR targeting different variable (V) regions (V2 and V3, V4 and V5, and V6 to V8) of the bacterial 16S rRNA gene with primers hybridizing to evolutionarily conserved flanking regions. SSCP analysis of single-stranded PCR products generated from 13 different bacterial species showed fewer bands with products containing V4-V5 (average, 1.7 bands per organism) than with V2-V3 (2.2 bands) and V6-V8 (2.3 bands). We found that the additional bands (>1 per organism) were caused by intraspecies operon heterogeneities or by more than one conformation of the same sequence. Community profiles, generated by PCR-SSCP from bacterial-cell consortia extracted from rhizospheres of field-grown maize (Zea mays), were analyzed by cloning and sequencing of the dominant bands. A total of 48 sequences could be attributed to 34 different strains from 10 taxonomical groups. Independent of the primer pairs, we found proteobacteria (alpha, beta, and gamma subgroups) and members of the genus Paenibacillus (low G+C gram-positive) to be the dominant organisms. Other groups, however, were only detected with single primer pairs. This study gives an example of how much the selection of different variable regions combined with different specificities of the flanking "universal" primers can affect a PCR-based microbial community analysis.     

A Method for Studying Protistan Diversity Using Massively Parallel Sequencing of V9 Hypervariable Regions of Small-Subunit Ribosomal RNA Genes

Background

Massively parallel pyrosequencing of amplicons from the V6 hypervariable regions of small-subunit (SSU) ribosomal RNA (rRNA) genes is commonly used to assess diversity and richness in bacterial and archaeal populations. Recent advances in pyrosequencing technology provide read lengths of up to 240 nucleotides. Amplicon pyrosequencing can now be applied to longer variable regions of the SSU rRNA gene including the V9 region in eukaryotes.

Methodology/Principal Findings

We present a protocol for the amplicon pyrosequencing of V9 regions for eukaryotic environmental samples for biodiversity inventories and species richness estimation. The International Census of Marine Microbes (ICoMM) and the Microbial Inventory Research Across Diverse Aquatic Long Term Ecological Research Sites (MIRADA-LTERs) projects are already employing this protocol for tag sequencing of eukaryotic samples in a wide diversity of both marine and freshwater environments.

Conclusions/Significance

Massively parallel pyrosequencing of eukaryotic V9 hypervariable regions of SSU rRNA genes provides a means of estimating species richness from deeply-sampled populations and for discovering novel species from the environment.     

Succession of microbial communities during hot composting as detected by PCR-single-strand-conformation polymorphism-based genetic profiles of small-subunit rRNA genes

A cultivation-independent technique for genetic profiling of PCR-amplified small-subunit rRNA genes (SSU rDNA) was chosen to characterize the diversity and succession of microbial communities during composting of an organic agricultural substrate. PCR amplifications were performed with DNA directly extracted from compost samples and with primers targeting either (i) the V4-V5 region of eubacterial 16S rRNA genes, (ii) the V3 region in the 16S rRNA genes of actinomycetes, or (iii) the V8-V9 region of fungal 18S rRNA genes. Homologous PCR products were converted to single-stranded DNA molecules by exonuclease digestion and were subsequently electrophoretically separated by their single-strand-conformation polymorphism (SSCP). Genetic profiles obtained by this technique showed a succession and increasing diversity of microbial populations with all primers. A total of 19 single products were isolated from the profiles by PCR reamplification and cloning. DNA sequencing of these molecular isolates showed similarities in the range of 92.3 to 100% to known gram-positive bacteria with a low or high G+C DNA content and to the SSU rDNA of gamma-Proteobacteria. The amplified 18S rRNA gene sequences were related to the respective gene regions of Candida krusei and Candida tropicalis. Specific molecular isolates could be attributed to different composting stages. The diversity of cultivated bacteria isolated from samples taken at the end of the composting process was low. A total of 290 isolates were related to only 6 different species. Two or three of these species were also detectable in the SSCP community profiles. Our study indicates that community SSCP profiles can be highly useful for the monitoring of bacterial diversity and community successions in a biotechnologically relevant process.     

Rapid identification of bacteria on the basis of polymerase chain reaction-amplified ribosomal DNA spacer polymorphisms.

To facilitate genus and species level identification of a broad range of bacteria without the requirement of presumptive identification, we have developed a unified set of primers and polymerase chain reaction conditions to amplify spacer regions between the 16S and 23S genes in the prokaryotic rRNA genetic loci. Spacer regions within these loci show a significant level of length and sequence polymorphism across both genus and species lines. A generic pair of priming sequences was selected for the amplification of these polymorphisms from highly conserved sequences in the 16S and 23S genes occurring adjacent to these polymorphic regions. This single set of primers and reaction conditions was used for the amplification of 16S-23S spacer regions for over 300 strains of bacteria belonging to eight genera and 28 species or serotypes, including Listeria, Staphylococcus, and Salmonella species and additional species related to these pathogenic organisms. When the spacer amplification products were resolved by electrophoresis, the resulting patterns could be used to distinguish all of the species of bacteria within the test group. Unique elements in the amplification product patterns generally clustered at the species level, although some genus-specific characteristics were also observed. On the basis of the results obtained with our test group of 300 bacterial strains, amplification of the 16S-23S ribosomal spacer region is a suitable process for generating a data base for use in a polymerase chain reaction-based identification method, which can be comprehensively applied to the bacterial kingdom.     

不同因素对六盘山地区甘肃鼢鼠肠道细菌多样性的影响

【目的】研究甘肃鼢鼠肠道细菌的群落结构和多样性,探讨不同环境因素对六盘山地区甘肃鼢鼠肠道细菌多样性的影响。【方法】通过Illumina高通量测序技术,对甘肃鼢鼠的36个盲肠样品进行16S rRNA V3–V4区的高通量测序,分析了肠道细菌多样性、丰度和群落结构,探讨地域、性别和季节等因素对甘肃鼢鼠肠道细菌多样性的影响。【结果】甘肃鼢鼠肠道细菌群落主要包括3个门,其中Firmicutes门占主导地位,其次是Bacteroidetes和Proteobacteri。在属水平,优势菌属分别为Oscillospira、Ruminococcus、Coprococcus和Desulfovibrio等。不同县(区)样品中,彭阳县、隆德县和泾源县三个县的甘肃鼢鼠肠道细菌菌群结构相似度较高,海原县甘肃鼢鼠肠道细菌群落结构在组内相似度高,与其他县(区)相似度低;雌性甘肃鼢鼠肠道细菌群落结构相似性高,而雄性细菌群落结构在样品间差异较大;甘肃鼢鼠肠道细菌菌群多样性秋季显著高于春季,细菌群落结构相似度秋季高于春季。【结论】不同地域、性别和季节因素对甘肃鼢鼠肠道细菌群落结构产生显著的影响,甘肃鼢鼠肠道细菌群落结构和多样性的变化对基于食源的季节性变化具有积极的响应。     

Analysis of the fecal microflora of human subjects consuming a probiotic product containing Lactobacillus rhamnosus DR20

The composition of the fecal microflora of 10 healthy subjects was monitored before (6-month control period), during (6-month test period), and after (3-month posttest period) the administration of a milk product containing Lactobacillus rhamnosus DR20 (daily dose, 1.6 x 10(9) lactobacilli). Monthly fecal samples were examined by a variety of methods, including bacteriological culture analysis, fluorescent in situ hybridization with group-specific DNA probes, denaturing gradient gel electrophoresis of the V2-V3 region of 16S rRNA genes amplified by PCR, gas-liquid chromatography, and bacterial enzyme activity analysis. The composition of the Lactobacillus population of each subject was analyzed by pulsed-field gel electrophoresis of bacterial DNA digests in order to differentiate between DR20 and other strains present in the samples. Representative isolates of lactobacilli were identified to the species level by sequencing the V2-V3 region of their 16S rRNA genes and comparing the sequences obtained (BLAST search) to sequences in the GenBank database. DR20 was detected in the feces of all of the subjects during the test period, but at different frequencies. The presence of DR20 among the numerically predominant strains was related to the presence or absence of a stable indigenous population of lactobacilli during the control period. Strain DR20 did not persist at levels of >10(2) cells per g in the feces of most of the subjects after consumption of the product ceased; the only exception was one subject in which this strain was detected for 2 months during the posttest period. We concluded that consumption of the DR20-containing milk product transiently altered the Lactobacillus and enterococcal contents of the feces of the majority of consumers without markedly affecting biochemical or other bacteriological factors.     

Denaturing gradient gel electrophoresis analysis with different primers of subgingival bacterial communities under mechanical debridement

DGGE of 16S rDNA is one of the most frequently used methods to study microbial communities. In this study, the DGGE profiles of different 16S rDNA regions of the periodontal pathogens Porphyromonas gingivalis, Fusobacterium nucleatum, and Prevotella nigrescens were investigated. The results suggested that V3-V5 and V6-V8 fragments may be suitable for community analysis of subgingival bacteria. Further analysis of subgingival samples with V3-V5 and V6-V8 regions as target fragments suggested that, in chronic periodontitis, re-colonization by periodontal bacteria with a population very similar to the baseline may occur by 6 weeks after mechanical debridement.     

PCR-DGGE analysis for the identification of microbial populations from Argentinean dry fermented sausages

Different PCR-DGGE protocols were evaluated to monitor fermentation process and to investigate bacterial communities developed in two artisanal Argentinean fermented sausages. Bacterial universal primers frequently used in PCR-denaturing gradient gel electrophoresis (DGGE) were evaluated. Lactic acid bacteria (LAB) and staphylococci species isolated from Tucumán sausages were used to determine the experimental conditions for PCR amplification and DGGE differentiation. Total microbial DNA extracted directly from both fermented sausages was subjected to DGGE analysis. PCR-DGGE results were different for each set of primers used. Primers Bact-0124f(GC)-Uni-0515r and V1f(GC)-V1r showed to be efficient to differentiate LAB and Staphylococcus cultures while the set V3f(GC)-Uni-0515r allowed to demonstrate the succession of different Lactobacillus and Staphylococcus species during ripening process. An intense band corresponding to Lactobacillus sakei was observed to be present in both samples. Staphylococcus saprophyticus was only observed in Tucumán sausage while a band identified as Brochothrix thermophacta was detected in Córdoba sausage. PCR-DGGE analysis of different 16S rDNA amplicons was able to discriminate between LAB and Gram-positive, coagulase-negative cocci, resulting an effective tool to establish the microbiota developed in artisanal dry sausages.     

The Relation between Oral Candida Load and Bacterial Microbiome Profiles in Dutch Older Adults

Currently there are no evidence-based ecological measures for prevention of overgrowth and subsequent infection by fungi in the oral cavity. The aim of this study was to increase our knowledge on fungal-bacterial ecological interactions. Salivary Candida abundance of 82 Dutch adults aged 58-80 years was established relative to the bacterial load by quantitative PCR analysis of the Internal Transcribed (ITS) region (Candida) and 16S rDNA gene (bacteria). The salivary microbiome was assessed using barcoded pyrosequencing of the bacterial hypervariable regions V5-V7 of 16S rDNA. Sequencing data was preprocessed by denoising and chimera removal, clustered in Operational Taxonomic Units (OTUs) and assigned to taxonomy. Both OTU-based (PCA, diversity statistics) and phylogeny-based analyses (UniFrac, PCoA) were performed. Saliva of Dutch older adults contained 0-4 × 10(8) CFU/mL Candida with a median Candida load of 0.06%. With increased Candida load the diversity of the salivary microbiome decreased significantly (p<0.001). Increase in the Candida load correlated positively with class Bacilli, and negatively with class Fusobacteria, Flavobacteria, and Bacteroidia. Microbiomes with high Candida load were less diverse and had a distinct microbial composition towards dominance by saccharolytic and acidogenic bacteria - streptococci. The control of the acidification of the oral environment may be a potential preventive measure for Candida outgrowth that should be evaluated in longitudinal clinical intervention trials.     

A core human microbiome as viewed through 16S rRNA sequence clusters

We explore the microbiota of 18 body sites in over 200 individuals using sequences amplified V1-V3 and the V3-V5 small subunit ribosomal RNA (16S) hypervariable regions as part of the NIH Common Fund Human Microbiome Project. The body sites with the greatest number of core OTUs, defined as OTUs shared amongst 95% or more of the individuals, were the oral sites (saliva, tongue, cheek, gums, and throat) followed by the nose, stool, and skin, while the vaginal sites had the fewest number of OTUs shared across subjects. We found that commonalities between samples based on taxonomy could sometimes belie variability at the sub-genus OTU level. This was particularly apparent in the mouth where a given genus can be present in many different oral sites, but the sub-genus OTUs show very distinct site selection, and in the vaginal sites, which are consistently dominated by the Lactobacillus genus but have distinctly different sub-genus V1-V3 OTU populations across subjects. Different body sites show approximately a ten-fold difference in estimated microbial richness, with stool samples having the highest estimated richness, followed by the mouth, throat and gums, then by the skin, nasal and vaginal sites. Richness as measured by the V1-V3 primers was consistently higher than richness measured by V3-V5. We also show that when such a large cohort is analyzed at the genus level, most subjects fit the stool "enterotype" profile, but other subjects are intermediate, blurring the distinction between the enterotypes. When analyzed at the finer-scale, OTU level, there was little or no segregation into stool enterotypes, but in the vagina distinct biotypes were apparent. Finally, we note that even OTUs present in nearly every subject, or that dominate in some samples, showed orders of magnitude variation in relative abundance emphasizing the highly variable nature across individuals.     

Effect of Primers Hybridizing to Different Evolutionarily Conserved Regions of the Small-Subunit rRNA Gene in PCR-Based Microbial Community Analyses and Genetic Profiling

Genetic profiling techniques of microbial communities based on PCR-amplified signature genes, such as denaturing gradient gel electrophoresis or single-strand-conformation polymorphism (SSCP) analysis, are normally done with PCR products of less than 500-bp. The most common target for diversity analysis, the small-subunit rRNA genes, however, are larger, and thus, only partial sequences can be analyzed. Here, we compared the results obtained by PCR targeting different variable (V) regions (V2 and V3, V4 and V5, and V6 to V8) of the bacterial 16S rRNA gene with primers hybridizing to evolutionarily conserved flanking regions. SSCP analysis of single-stranded PCR products generated from 13 different bacterial species showed fewer bands with products containing V4-V5 (average, 1.7 bands per organism) than with V2-V3 (2.2 bands) and V6-V8 (2.3 bands). We found that the additional bands (>1 per organism) were caused by intraspecies operon heterogeneities or by more than one conformation of the same sequence. Community profiles, generated by PCR-SSCP from bacterial-cell consortia extracted from rhizospheres of field-grown maize (Zea mays), were analyzed by cloning and sequencing of the dominant bands. A total of 48 sequences could be attributed to 34 different strains from 10 taxonomical groups. Independent of the primer pairs, we found proteobacteria (α, β, and γ subgroups) and members of the genus Paenibacillus (low G+C gram-positive) to be the dominant organisms. Other groups, however, were only detected with single primer pairs. This study gives an example of how much the selection of different variable regions combined with different specificities of the flanking “universal” primers can affect a PCR-based microbial community analysis.     

Assessing SSU rDNA Barcodes in Foraminifera: A Case Study using Bolivina quadrata

The Small Subunit Ribosomal RNA gene (SSU rDNA) is a widely used tool to reconstruct phylogenetic relationships among foraminiferal species. Recently, the highly variable regions of this gene have been proposed as DNA barcodes to identify foraminiferal species. However, the resolution of these barcodes has not been well established, yet. In this study, we evaluate four SSU rDNA hypervariable regions (37/f, 41/f, 43/e, and 45/e) as DNA barcodes to distinguish among species of the genus Bolivina, with particular emphasis on Bolivina quadrata for which ten new sequences ( KY468817 – KY468826 ) were obtained during this study. Our analyses show that a single SSU rDNA hypervariable sequence is insufficient to resolve all Bolivina species and that some regions (37/f and 41/f) are more useful than others (43/e and 45/e) to distinguish among closely related species. In addition, polymorphism analyses reveal a high degree of variability. In the context of barcoding studies, these results emphasize the need to assess the range of intraspecific variability of DNA barcodes prior to their application to identify foraminiferal species in environmental samples; our results also highlight the possibility that a longer SSU rDNA region might be required to distinguish among species belonging to the same taxonomic group (i.e. genus).     

Geographical location determines the population structure in phyllosphere microbial communities of a salt-excreting desert tree

The leaf surfaces of Tamarix, a salt-secreting desert tree, harbor a diverse community of microbial epiphytes. This ecosystem presents a unique combination of ecological characteristics and imposes a set of extreme stress conditions. The composition of the microbial community along ecological gradients was studied from analyses of microbial richness and diversity in the phyllosphere of three Tamarix species in the Mediterranean and Dead Sea regions in Israel and in two locations in the United States. Over 200,000 sequences of the 16S V6 and 18S V9 hypervariable regions revealed a diverse community, with 788 bacterial and 64 eukaryotic genera but only one archaeal genus. Both geographic location and tree species were determinants of microbial community structures, with the former being more dominant. Tree leaves of all three species in the Mediterranean region were dominated by Halomonas and Halobacteria, whereas trees from the Dead Sea area were dominated by Actinomycetales and Bacillales. Our findings demonstrate that microbial phyllosphere communities on different Tamarix species are highly similar in the same locale, whereas trees of the same species that grow in different climatic regions host distinct microbial communities.