Review and re-analysis of domain-specific 16S primers

The Polymerase Chain Reaction (PCR) has facilitated the detection of unculturable microorganisms in virtually any environmental source and has thus been used extensively in the assessment of environmental microbial diversity. This technique relies on the assumption that the gene sequences present in the environment are complementary to the "universal" primers used in their amplification. The recent discovery of new taxa with 16S rDNA sequences not complementary to standard universal primers suggests that current 16S rDNA libraries are not representative of true prokaryotic biodiversity. Here we re-assess the specificity of commonly used 16S rRNA gene primers and present these data in tabular form designed as a tool to aid simple analysis, selection and implementation. In addition, we present two new primer pairs specifically designed for effective "universal" Archaeal 16S rDNA sequence amplification. These primers are found to amplify sequences from Crenarchaeote and Euryarchaeote type strains and environmental DNA.     

Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies

16S ribosomal RNA gene (rDNA) amplicon analysis remains the standard approach for the cultivation-independent investigation of microbial diversity. The accuracy of these analyses depends strongly on the choice of primers. The overall coverage and phylum spectrum of 175 primers and 512 primer pairs were evaluated in silico with respect to the SILVA 16S/18S rDNA non-redundant reference dataset (SSURef 108 NR). Based on this evaluation a selection of ‘best available’ primer pairs for Bacteria and Archaea for three amplicon size classes (100–400, 400–1000, ≥1000 bp) is provided. The most promising bacterial primer pair (S-D-Bact-0341-b-S-17/S-D-Bact-0785-a-A-21), with an amplicon size of 464 bp, was experimentally evaluated by comparing the taxonomic distribution of the 16S rDNA amplicons with 16S rDNA fragments from directly sequenced metagenomes. The results of this study may be used as a guideline for selecting primer pairs with the best overall coverage and phylum spectrum for specific applications, therefore reducing the bias in PCR-based microbial diversity studies.     

Cross-Kingdom Amplification Using Bacteria-Specific Primers: Complications for Studies of Coral Microbial Ecology

PCR amplification of pure bacterial DNA is vital to the study of bacterial interactions with corals. Commonly used Bacteria-specific primers 8F and 27F paired with the universal primer 1492R amplify both eukaryotic and prokaryotic rRNA genes. An alternative primer set, 63F/1542R, is suggested to resolve this problem.     

Design and Experimental Application of a Novel Non-Degenerate Universal Primer Set that Amplifies Prokaryotic 16S rRNA Genes with a Low Possibility to Amplify Eukaryotic rRNA Genes

The deep sequencing of 16S rRNA genes amplified by universal primers has revolutionized our understanding of microbial communities by allowing the characterization of the diversity of the uncultured majority. However, some universal primers also amplify eukaryotic rRNA genes, leading to a decrease in the efficiency of sequencing of prokaryotic 16S rRNA genes with possible mischaracterization of the diversity in the microbial community. In this study, we compared 16S rRNA gene sequences from genome-sequenced strains and identified candidates for non-degenerate universal primers that could be used for the amplification of prokaryotic 16S rRNA genes. The 50 identified candidates were investigated to calculate their coverage for prokaryotic and eukaryotic rRNA genes, including those from uncultured taxa and eukaryotic organelles, and a novel universal primer set, 342F-806R, covering many prokaryotic, but not eukaryotic, rRNA genes was identified. This primer set was validated by the amplification of 16S rRNA genes from a soil metagenomic sample and subsequent pyrosequencing using the Roche 454 platform. The same sample was also used for pyrosequencing of the amplicons by employing a commonly used primer set, 338F-533R, and for shotgun metagenomic sequencing using the Illumina platform. Our comparison of the taxonomic compositions inferred by the three sequencing experiments indicated that the non-degenerate 342F-806R primer set can characterize the taxonomic composition of the microbial community without substantial bias, and is highly expected to be applicable to the analysis of a wide variety of microbial communities.     

Novel primers for 16S rRNA-based archaeal community analyses in environmental samples

Next generation sequencing technologies for in depth analyses of complex microbial communities rely on rational primer design based on up-to-date reference databases. Most of the 16S rRNA-gene based analyses of environmental Archaea community composition use PCR primers developed from small data sets several years ago, making an update long overdue. Here we present a new set of archaeal primers targeting the 16S rRNA gene designed from 8500 aligned archaeal sequences in the SILVA database. The primers 340F-1000R showed a high archaeal specificity (< 1% bacteria amplification) covering 93 and 97% of available sequences for Crenarchaeota and Euryarchaeota respectively. In silico tests of the primers revealed at least 38% higher coverage for Archaea compared to other commonly used primers. Empirical tests with clone libraries confirmed the high specificity of the primer pair to Archaea in three biomes: surface waters in the Arctic Ocean, the pelagic zone of a temperate lake and a methanogenic bioreactor. The clone libraries featured both Euryarchaeota and Crenarchaeota in variable proportions and revealed dramatic differences in the archaeal community composition and minimal phylogenetic overlap between samples.     

Diversity of Methane-Cycling Archaea in Hydrothermal Sediment Investigated by General and Group-Specific PCR Primers

The zonation of anaerobic methane-cycling Archaea in hydrothermal sediment of Guaymas Basin was studied by general primer pairs (mcrI, ME1/ME2, mcrIRD) targeting the alpha subunit of methyl coenzyme M reductase gene (mcrA) and by new group-specific mcrA and 16S rRNA gene primer pairs. The mcrIRD primer pair outperformed the other general mcrA primer pairs in detection sensitivity and phylogenetic coverage. Methanotrophic ANME-1 Archaea were the only group detected with group-specific primers only. The detection of 14 mcrA lineages surpasses the diversity previously found in this location. Most phylotypes have high sequence similarities to hydrogenotrophs, methylotrophs, and anaerobic methanotrophs previously detected at Guaymas Basin or at hydrothermal vents, cold seeps, and oil reservoirs worldwide. Additionally, five mcrA phylotypes belonging to newly defined lineages are detected. Two of these belong to deeply branching new orders, while the others are new species or genera of Methanopyraceae and Methermicoccaceae. Downcore diversity decreases from all groups detected in the upper 6 cm (∼2 to 40°C, sulfate measurable to 4 cm) to only two groups below 6 cm (>40°C). Despite the presence of hyperthermophilic genera (Methanopyrus, Methanocaldococcus) in cooler surface strata, no genes were detected below 10 cm (≥60°C). While mcrA-based and 16S rRNA gene-based community compositions are generally congruent, the deeply branching mcrA cannot be assigned to specific 16S rRNA gene lineages. Our study indicates that even among well-studied metabolic groups and in previously characterized model environments, major evolutionary branches are overlooked. Detecting these groups by improved molecular biological methods is a crucial first step toward understanding their roles in nature.     

Anaerobic Microbial Communities in Lake Pavin, a Unique Meromictic Lake in France

The Bacteria and Archaea from the meromictic Lake Pavin were analyzed in samples collected along a vertical profile in the anoxic monimolimnion and were compared to those in samples from the oxic mixolimnion. Nine targeted 16S rRNA oligonucleotide probes were used to assess the distribution of Bacteria and Archaea and to investigate the in situ occurrence of sulfate-reducing bacteria and methane-producing Archaea involved in the terminal steps of the anaerobic degradation of organic material. The diversity of the complex microbial communities was assessed from the 16S rRNA polymorphisms present in terminal restriction fragment (TRF) depth patterns. The densities of the microbial community increased in the anoxic layer, and Archaea detected with probe ARCH915 represented the largest microbial group in the water column, with a mean Archaea/Eubacteria ratio of 1.5. Terminal restriction fragment length polymorphism (TRFLP) analysis revealed an elevated archaeal and bacterial phylotype richness in anoxic bottom-water samples. The structure of the Archaea community remained rather homogeneous, while TRFLP patterns for the eubacterial community revealed a heterogeneous distribution of eubacterial TRFs.     

Phylogenetic Diversity Analysis of Subterranean Hot Springs in Iceland

Geothermal energy has been harnessed and used for domestic heating in Iceland. In wells that are typically drilled to a depth of 1,500 to 2,000 m, the temperature of the source water is 50 to 130°C. The bottoms of the boreholes can therefore be regarded as subterranean hot springs and provide a unique opportunity to study the subterranean biosphere. Large volumes of geothermal fluid from five wells and a mixture of geothermal water from 50 geothermal wells (hot tap water) were sampled and concentrated through a 0.2-μm-pore-size filter. Cells were observed in wells RG-39 (91.4°C) and MG-18 (71.8°C) and in hot tap water (76°C), but no cells were detected in wells SN-4, SN-5 (95 to 117°C), and RV-5 (130°C). Archaea and Bacteria were detected by whole-cell fluorescent in situ hybridization. DNAs were extracted from the biomass, and small-subunit rRNA genes (16S rDNAs) were amplified by PCR using primers specific for the Archaea and Bacteria domains. The PCR products were cloned and sequenced. The sequence analysis showed 11 new operational taxonomic units (OTUs) out of 14, 3 of which were affiliated with known surface OTUs. Samples from RG-39 and hot tap water were inoculated into enrichment media and incubated at 65 and 85°C. Growth was observed only in media based on geothermal water. 16S rDNA analysis showed enrichments dominated with Desulfurococcales relatives. Two strains belonging to Desulfurococcus mobilis and to the Thermus/Deinococcus group were isolated from borehole RG-39. The results indicate that subsurface volcanic zones are an environment that provides a rich subsurface for novel thermophiles.     

Phylogenetic Characterization and In Situ Detection of a Cytophaga-Flexibacter-Bacteroides Phylogroup Bacterium in Tuber borchii Vittad. Ectomycorrhizal Mycelium

Mycorrhizal ascomycetous fungi are obligate ectosymbionts that colonize the roots of gymnosperms and angiosperms. In this paper we describe a straightforward approach in which a combination of morphological and molecular methods was used to survey the presence of potentially endo- and epiphytic bacteria associated with the ascomycetous ectomycorrhizal fungus Tuber borchii Vittad. Universal eubacterial primers specific for the 5′ and 3′ ends of the 16S rRNA gene (16S rDNA) were used for PCR amplification, direct sequencing, and phylogenetic analyses. The 16S rDNA was amplified directly from four pure cultures of T. borchii Vittad. mycelium. A nearly full-length sequence of the gene coding for the prokaryotic small-subunit rRNA was obtained from each T. borchii mycelium studied. The 16S rDNA sequences were almost identical (98 to 99% similarity), and phylogenetic analysis placed them in a single unique rRNA branch belonging to the Cytophaga-Flexibacter-Bacteroides (CFB) phylogroup which had not been described previously. In situ detection of the CFB bacterium in the hyphal tissue of the fungus T. borchii was carried out by using 16S rRNA-targeted oligonucleotide probes for the eubacterial domain and the Cytophaga-Flexibacter phylum, as well as a probe specifically designed for the detection of this mycelium-associated bacterium. Fluorescent in situ hybridization showed that all three of the probes used bound to the mycelium tissue. This study provides the first direct visual evidence of a not-yet-cultured CFB bacterium associated with a mycorrhizal fungus of the genus Tuber.     

Advantage of Using Inosine at the 3′ Termini of 16S rRNA Gene Universal Primers for the Study of Microbial Diversity

To overcome the shortcomings of universal 16S rRNA gene primers 8F and 907R when studying the diversity of complex microbial communities, the 3′ termini of both primers were replaced with inosine. A comparison of the clone libraries derived using both primer sets showed seven bacterial phyla amplified by the altered primer set (8F-I/907R-I) whereas the original set amplified sequences belonging almost exclusively to Proteobacteria (95.8%). Sequences belonging to Firmicutes (42.6%) and Thermotogae (9.3%) were more abundant in a library obtained by using 8F-I/907R-I at a PCR annealing temperature of 54°C, while Proteobacteria sequences were more frequent (62.7%) in a library obtained at 50°C, somewhat resembling the result obtained using the original primer set. The increased diversity revealed by using primers 8F-I/907R-I confirms the usefulness of primers with inosine at the 3′ termini in studying the microbial diversity of environmental samples.     

Phylogeny of Dissimilatory Sulfite Reductases Supports an Early Origin of Sulfate Respiration

Microorganisms that use sulfate as a terminal electron acceptor for anaerobic respiration play a central role in the global sulfur cycle. Here, we report the results of comparative sequence analysis of dissimilatory sulfite reductase (DSR) genes from closely and distantly related sulfate-reducing organisms to infer the evolutionary history of DSR. A 1.9-kb DNA region encoding most of the α and β subunits of DSR could be recovered only from organisms capable of dissimilatory sulfate reduction with a PCR primer set targeting highly conserved regions in these genes. All DNA sequences obtained were highly similar to one another (49 to 89% identity), and their inferred evolutionary relationships were nearly identical to those inferred on the basis of 16S rRNA. We conclude that the high similarity of bacterial and archaeal DSRs reflects their common origin from a conserved DSR. This ancestral DSR was either present before the split between the domains Bacteria, Archaea, and Eucarya or laterally transferred between Bacteria and Archaea soon after domain divergence. Thus, if the physiological role of the DSR was constant over time, then early ancestors of Bacteria and Archaea already possessed a key enzyme of sulfate and sulfite respiration.     

Population Structure and Phylogenetic Characterization of Marine Benthic Archaea in Deep-Sea Sediments

During the past few years Archaea have been recognized as a widespread and significant component of marine picoplankton assemblages and, more recently, the presence of novel archaeal phylogenetic lineages has been reported in coastal marine benthic environments. We investigated the relative abundance, vertical distribution, phylogenetic composition, and spatial variability of Archaea in deep-sea sediments collected from several stations in the Atlantic Ocean. Quantitative oligonucleotide hybridization experiments indicated that the relative abundance of archaeal 16S rRNA in deep-sea sediments (1500 m deep) ranged from about 2.5 to 8% of the total prokaryotic rRNA. Clone libraries of PCR-amplified archaeal rRNA genes (rDNA) were constructed from 10 depth intervals obtained from sediment cores collected at depths of 1,500, 2,600, and 4,500 m. Phylogenetic analysis of rDNA sequences revealed the presence of a complex archaeal population structure, whose members could be grouped into discrete phylogenetic lineages within the two kingdoms, Crenarchaeota and Euryarchaeota. Comparative denaturing gradient gel electrophoresis profile analysis of archaeal 16S rDNA V3 fragments revealed a significant depth-related variability in the composition of the archaeal population.     

Does peptide-nucleic acid (PNA) clamping of host plant DNA benefit ITS1 amplicon-based characterization of the fungal endophyte community?

Fungal endophyte community amplicon sequencing can lose a significant number of informative reads due to host-plant co-amplification. Blocking of plant-specific sequences with peptide nucleic acid (PNA) clamps has been shown to improve metrics of detected microbial diversity in studies targeting 16S and 18S regions of rRNA genes. However, PNA clamping has not been applied to the plant ITS region of rRNA gene – a widely accepted fungal marker. By applying PNA clamping technique to ITS amplicon sequencing of the endophytic fungal community of elderberry this study shows that PNA clamping significantly reduces host-plant co-amplification with the universal ITS1/ITS4 primer set. However, PNA clamping in combination with the discriminatory ITS1F/ITS2 primer set did not improve the metrics of fungal endophyte community ITS amplicon Illumina sequencing. This study shows that PNA clamping does not add practical benefit to taxonomic profiling of plant-associated fungal communities if the primers are already specific enough to exclude amplification of host DNA.     

Influence of DNA Extraction Method, 16S rRNA Targeted Hypervariable Regions,and Sample Origin on Microbial Diversity Detected by 454 Pyrosequencing in Marine Chemosynthetic Ecosystems

Next-generation sequencing (NGS) opens up exciting possibilities for improving our knowledge of environmental microbial diversity, allowing rapid and cost-effective identification of both cultivated and uncultivated microorganisms. However, library preparation, sequencing, and analysis of the results can provide inaccurate representations of the studied community compositions. Therefore, all these steps need to be taken into account carefully. Here we evaluated the effects of DNA extraction methods, targeted 16S rRNA hypervariable regions, and sample origins on the diverse microbes detected by 454 pyrosequencing in marine cold seep and hydrothermal vent sediments. To assign the reads with enough taxonomic precision, we built a database with about 2,500 sequences from Archaea and Bacteria from deep-sea marine sediments, affiliated according to reference publications in the field. Thanks to statistical and diversity analyses as well as inference of operational taxonomic unit (OTU) networks, we show that (i) while DNA extraction methods do not seem to affect the results for some samples, they can lead to dramatic changes for others; and (ii) the choice of amplification and sequencing primers also considerably affects the microbial community detected in the samples. Thereby, very different proportions of pyrosequencing reads were obtained for some microbial lineages, such as the archaeal ANME-1, ANME-2c, and MBG-D and deltaproteobacterial subgroups. This work clearly indicates that the results from sequencing-based analyses, such as pyrosequencing, should be interpreted very carefully. Therefore, the combination of NGS with complementary approaches, such as fluorescence in situ hybridization (FISH)/catalyzed reporter deposition (CARD)-FISH or quantitative PCR (Q-PCR), would be desirable to gain a more comprehensive picture of environmental microbial communities.     

Microbial Diversity in Maras Salterns, a Hypersaline Environment in the Peruvian Andes

Maras salterns are located 3,380 m above sea level in the Peruvian Andes. These salterns consist of more than 3,000 little ponds which are not interconnected and act as crystallizers where salt precipitates. These ponds are fed by hypersaline spring water rich in sodium and chloride. The microbiota inhabiting these salterns was examined by fluorescence in situ hybridization (FISH), 16S rRNA gene clone library analysis, and cultivation techniques. The total counts per milliliter in the ponds were around 2 × 10⁶ to 3 × 10⁶ cells/ml, while the spring water contained less than 100 cells/ml and did not yield any detectable FISH signal. The microbiota inhabiting the ponds was dominated (80 to 86% of the total counts) by Archaea, while Bacteria accounted for 10 to 13% of the 4′,6′-diamidino-2-phenylindole (DAPI) counts. A total of 239 16S rRNA gene clones were analyzed (132 Archaea clones and 107 Bacteria clones). According to the clone libraries, the archaeal assemblage was dominated by microorganisms related to the cosmopolitan square archaeon “Haloquadra walsbyi,” although a substantial number of the sequences in the libraries (31% of the 16S rRNA gene archaeal clones) were related to Halobacterium sp., which is not normally found in clone libraries from solar salterns. All the bacterial clones were closely related to each other and to the γ-proteobacterium “Pseudomonas halophila” DSM 3050. FISH analysis with a probe specific for this bacterial assemblage revealed that it accounted for 69 to 76% of the total bacterial counts detected with a Bacteria-specific probe. When pond water was used to inoculate solid media containing 25% total salts, both extremely halophilic Archaea and Bacteria were isolated. Archaeal isolates were not related to the isolates in clone libraries, although several bacterial isolates were very closely related to the “P. halophila” cluster found in the libraries. As observed for other hypersaline environments, extremely halophilic bacteria that had ecological relevance seemed to be easier to culture than their archaeal counterparts.     

Vertical Distribution of Archaea and Bacteria in a Meromictic Lake as Determined by Fluorescent In Situ Hybridization

The prokaryotic cells distribution in the water column of the coastal saline meromictic Lake Faro (Messina, Italy) was investigated by microscopic counting techniques. Water samples were collected at a central station from the surface to the bottom, when waters were characterized by a marked stratification. A “red-water” layer, caused by a dense growth of photosynthetic sulfur bacteria, was present at a depth of 15 m, defining a transition area between oxic (mixolimnion) and anoxic (monimolimnion) layers. Fluorescently labeled 16S rRNA oligonucleotide, group-specific probes were used to determine the abundance of Bacteria and Archaea, and their subgroups, Green Sulfur Bacteria (GSB), Sulfate Reducing Bacteria (SRB), Cyanobacteria and Chromatium okenii, and Crenarchaeota and Euryarchaeota, as key elements of the microbial community. Bacteria decreased from surface to bottom, while Archaea increased with depth and reached the maximum value at 30 m, where they outnumbered the Bacteria. Bacteria and picophytoplankton prevailed in the mixolimnion. At the chemocline high numbers of prokaryotic cells were present, mainly represented by Cyanobacteria, Chromatium okenii and Euryarchaeota. GSB, SRB, and Crenarchaeota prevailed below the chemocline. Although Archaea constitute a minor fraction of microbial community, they could represent active contributors to the meromictic Lake Faro ecosystem.     

Use of a Hierarchical Oligonucleotide Primer Extension Approach for Multiplexed Relative Abundance Analysis of Methanogens in Anaerobic Digestion Systems

In this study, we established a rapid multiplex method to detect the relative abundances of amplified 16S rRNA genes from known cultivatable methanogens at hierarchical specificities in anaerobic digestion systems treating industrial wastewater and sewage sludge. The method was based on the hierarchical oligonucleotide primer extension (HOPE) technique and combined with a set of 27 primers designed to target the total archaeal populations and methanogens from 22 genera within 4 taxonomic orders. After optimization for their specificities and detection sensitivity under the conditions of multiple single-nucleotide primer extension reactions, the HOPE approach was applied to analyze the methanogens in 19 consortium samples from 7 anaerobic treatment systems (i.e., 513 reactions). Among the samples, the methanogen populations detected with order-level primers accounted for >77.2% of the PCR-amplified 16S rRNA genes detected using an Archaea-specific primer. The archaeal communities typically consisted of 2 to 7 known methanogen genera within the Methanobacteriales, Methanomicrobiales, and Methanosarcinales and displayed population dynamic and spatial distributions in anaerobic reactor operations. Principal component analysis of the HOPE data further showed that the methanogen communities could be clustered into 3 distinctive groups, in accordance with the distribution of the Methanosaeta, Methanolinea, and Methanomethylovorans, respectively. This finding suggested that in addition to acetotrophic and hydrogenotrophic methanogens, the methylotrophic methanogens might play a key role in the anaerobic treatment of industrial wastewater. Overall, the results demonstrated that the HOPE approach is a specific, rapid, and multiplexing platform to determine the relative abundances of targeted methanogens in PCR-amplified 16S rRNA gene products.     

Quantitative multiplexing analysis of PCR-amplified ribosomal RNA genes by hierarchical oligonucleotide primer extension reaction

A method, termed hierarchical oligonucleotide primer extension (HOPE), is developed for quantitative, multiplexing detection of DNA targets present in PCR-amplified community 16S rRNA genes. It involves strand extension reaction and multiple oligonucleotide primers modified with different lengths of polyA at the 5′ end and targeting 16S rRNA genes at different phylogenetic specificities. On annealing to the targets, these primers are extended with a single fluorescently labeled dideoxynucleoside triphosphate or a dye-terminator. Using a DNA autosequencer, these extended primers are separated and identified by size and dye color, and quantified and normalized based on the fluorescence intensities and internal size standards. Using a primer-to-target ratio >1000, constant primer extension efficiencies can be obtained with individual primers to establish a ‘calibration factor’ between individual primers and a universal or domain-specific primer, providing the relative abundance of targeted rRNA genes with respect to total rRNA genes. HOPE up to 10-plexing is demonstrated to correctly identify 20 different bacterial strains, and quantify different Bacteroides spp. in 16S rRNA gene amplicons from different model bacteria mixtures and the influent and effluent of a wastewater treatment plant. Single mismatch discrimination with detection sensitivity of a target down to 0.01–0.05% of total DNA template is achieved.     

16S/18S/ITS扩增子高通量测序引物的生物信息学评估和改进

【背景】在过去的十几年里,基于核糖体RNA基因的扩增子测序技术被广泛用于各种生态系统中微生物群落的多样性检测。扩增子测序的使用极大地促进了土壤、水体、空气等环境中微生物生态的相关研究。【目的】随着高通量测序技术的不断发展和参考数据库的不断更新,针对不同的环境样本的引物选择和改进仍然需要更深入的校验。【方法】本文收集了目前在微生物群落研究中被广泛采用的标记基因扩增通用引物,包括16S rRNA基因扩增常用的8对通用引物和2对古菌引物、9对真菌转录间隔区(internal transcribed spacer,ITS)基因扩增引物,以及18S rRNA基因扩增的4对真核微生物通用引物和1对真菌特异性引物。这些引物中包括了地球微生物组计划(Earth Microbiome Project,EMP)推荐的2对16S rRNA基因扩增引物、1对ITS1基因扩增引物和1对18S rRNA基因扩增引物。采用最近更新的标准数据库对这些引物进行了覆盖度和特异性评价。【结果】EMP推荐的引物依然具有较高的覆盖度,而其他引物在覆盖度及对特定环境或类群的特异性上也各有特点。此外,最近有研究对这些通用引物进行了一些改进,而我们也发现,一个碱基的变化都可能会导致评价结果或扩增产物发生明显变化,简并碱基的引入既可以覆盖更多的物种,但同时也会在一定程度上降低关注物种的特异性。【结论】研究结果为微生态研究中标记基因的引物选择提供了一个广泛的指导,但在关注具体科学问题时,引物的选择仍需数据指导与实验尝试。