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.     

Divergence thresholds and divergent biodiversity estimates: can metabarcoding reliably describe zooplankton communities?

DNA metabarcoding is a promising method for describing communities and estimating biodiversity. This approach uses high‐throughput sequencing of targeted markers to identify species in a complex sample. By convention, sequences are clustered at a predefined sequence divergence threshold (often 3%) into operational taxonomic units (OTUs) that serve as a proxy for species. However, variable levels of interspecific marker variation across taxonomic groups make clustering sequences from a phylogenetically diverse dataset into OTUs at a uniform threshold problematic. In this study, we use mock zooplankton communities to evaluate the accuracy of species richness estimates when following conventional protocols to cluster hypervariable sequences of the V4 region of the small subunit ribosomal RNA gene (18S) into OTUs. By including individually tagged single specimens and “populations” of various species in our communities, we examine the impact of intra‐ and interspecific diversity on OTU clustering. Communities consisting of single individuals per species generated a correspondence of 59–84% between OTU number and species richness at a 3% divergence threshold. However, when multiple individuals per species were included, the correspondence between OTU number and species richness dropped to 31–63%. Our results suggest that intraspecific variation in this marker can often exceed 3%, such that a single species does not always correspond to one OTU. We advocate the need to apply group‐specific divergence thresholds when analyzing complex and taxonomically diverse communities, but also encourage the development of additional filtering steps that allow identification of artifactual rRNA gene sequences or pseudogenes that may generate spurious OTUs.     

Analysis of the gut microbiota in the old order amish and its relation to the metabolic syndrome

Obesity has been linked to the human gut microbiota; however, the contribution of gut bacterial species to the obese phenotype remains controversial because of conflicting results from studies in different populations. To explore the possible dysbiosis of gut microbiota in obesity and its metabolic complications, we studied men and women over a range of body mass indices from the Old Order Amish sect, a culturally homogeneous Caucasian population of Central European ancestry. We characterized the gut microbiota in 310 subjects by deep pyrosequencing of bar-coded PCR amplicons from the V1-V3 region of the 16S rRNA gene. Three communities of interacting bacteria were identified in the gut microbiota, analogous to previously identified gut enterotypes. Neither BMI nor any metabolic syndrome trait was associated with a particular gut community. Network analysis identified twenty-two bacterial species and four OTUs that were either positively or inversely correlated with metabolic syndrome traits, suggesting that certain members of the gut microbiota may play a role in these metabolic derangements.     

Exploration of bacterial community classes in major human habitats

Background

Determining bacterial abundance variation is the first step in understanding bacterial similarity between individuals. Categorization of bacterial communities into groups or community classes is the subsequent step in describing microbial distribution based on abundance patterns. Here, we present an analysis of the groupings of bacterial communities in stool, nasal, skin, vaginal and oral habitats in a healthy cohort of 236 subjects from the Human Microbiome Project.

Results

We identify distinct community group patterns in the anterior nares, four skin sites, and vagina at the genus level. We also confirm three enterotypes previously identified in stools. We identify two clusters with low silhouette values in most oral sites, in which bacterial communities are more homogeneous. Subjects sharing a community class in one habitat do not necessarily share a community class in another, except in the three vaginal sites and the symmetric habitats of the left and right retroauricular creases. Demographic factors, including gender, age, and ethnicity, significantly influence community composition in several habitats. Community classes in the vagina, retroauricular crease and stool are stable over approximately 200 days.

Conclusion

The community composition, association of demographic factors with community classes, and demonstration of community stability deepen our understanding of the variability and dynamics of human microbiomes. This also has significant implications for experimental designs that seek microbial correlations with clinical phenotypes.     

PCR biases distort bacterial and archaeal community structure in pyrosequencing datasets

As 16S rRNA gene targeted massively parallel sequencing has become a common tool for microbial diversity investigations, numerous advances have been made to minimize the influence of sequencing and chimeric PCR artifacts through rigorous quality control measures. However, there has been little effort towards understanding the effect of multi-template PCR biases on microbial community structure. In this study, we used three bacterial and three archaeal mock communities consisting of, respectively, 33 bacterial and 24 archaeal 16S rRNA gene sequences combined in different proportions to compare the influences of (1) sequencing depth, (2) sequencing artifacts (sequencing errors and chimeric PCR artifacts), and (3) biases in multi-template PCR, towards the interpretation of community structure in pyrosequencing datasets. We also assessed the influence of each of these three variables on α- and β-diversity metrics that rely on the number of OTUs alone (richness) and those that include both membership and the relative abundance of detected OTUs (diversity). As part of this study, we redesigned bacterial and archaeal primer sets that target the V3-V5 region of the 16S rRNA gene, along with multiplexing barcodes, to permit simultaneous sequencing of PCR products from the two domains. We conclude that the benefits of deeper sequencing efforts extend beyond greater OTU detection and result in higher precision in β-diversity analyses by reducing the variability between replicate libraries, despite the presence of more sequencing artifacts. Additionally, spurious OTUs resulting from sequencing errors have a significant impact on richness or shared-richness based α- and β-diversity metrics, whereas metrics that utilize community structure (including both richness and relative abundance of OTUs) are minimally affected by spurious OTUs. However, the greatest obstacle towards accurately evaluating community structure are the errors in estimated mean relative abundance of each detected OTU due to biases associated with multi-template PCR reactions.     

Global patterns of diversity and community structure in marine bacterioplankton

Because of their small size, great abundance and easy dispersal, it is often assumed that marine planktonic microorganisms have a ubiquitous distribution that prevents any structured assembly into local communities. To challenge this view, marine bacterioplankton communities from coastal waters at nine locations distributed world-wide were examined through the use of comprehensive clone libraries of 16S ribosomal RNA genes, used as operational taxonomic units (OTU). Our survey and analyses show that there were marked differences in the composition and richness of OTUs between locations. Remarkably, the global marine bacterioplankton community showed a high degree of endemism, and conversely included few cosmopolitan OTUs. Our data were consistent with a latitudinal gradient of OTU richness. We observed a positive relationship between the relative OTU abundances and their range of occupation, i.e. cosmopolitans had the largest population sizes. Although OTU richness differed among locations, the distributions of the major taxonomic groups represented in the communities were analogous, and all local communities were similarly structured and dominated by a few OTUs showing variable taxonomic affiliations. The observed patterns of OTU richness indicate that similar evolutionary and ecological processes structured the communities. We conclude that marine bacterioplankton share many of the biogeographical and macroecological features of macroscopic organisms. The general processes behind those patterns are likely to be comparable across taxa and major global biomes.     

Implications of pyrosequencing error correction for biological data interpretation

There has been a rapid proliferation of approaches for processing and manipulating second generation DNA sequence data. However, users are often left with uncertainties about how the choice of processing methods may impact biological interpretation of data. In this report, we probe differences in output between two different processing pipelines: a de-noising approach using the AmpliconNoise algorithm for error correction, and a standard approach using quality filtering and preclustering to reduce error. There was a large overlap in reads culled by each method, although AmpliconNoise removed a greater net number of reads. Most OTUs produced by one method had a clearly corresponding partner in the other. Although each method resulted in OTUs consisting entirely of reads that were culled by the other method, there were many more such OTUs formed in the standard pipeline. Total OTU richness was reduced by AmpliconNoise processing, but per-sample OTU richness, diversity and evenness were increased. Increases in per-sample richness and diversity may be a result of AmpliconNoise processing producing a more even OTU rank-abundance distribution. Because communities were randomly subsampled to equalize sample size across communities, and because rare sequence variants are less likely to be selected during subsampling, fewer OTUs were lost from individual communities when subsampling AmpliconNoise-processed data. In contrast to taxon-based diversity estimates, phylogenetic diversity was reduced even on a per-sample basis by de-noising, and samples switched widely in diversity rankings. This work illustrates the significant impacts of processing pipelines on the biological interpretations that can be made from pyrosequencing surveys. This study provides important cautions for analyses of contemporary data, for requisite data archiving (processed vs. non-processed data), and for drawing comparisons among studies performed using distinct data processing pipelines.     

基于高通量测序的鱼菜共生池塘与普通池塘微生物群落结构比较

鱼菜共生(Aquaponics)作为一种可持续、循环型、零排放的新型的复合耕作体系, 因具有良好的生态环境效应和食品安全保障, 成为解决农业生态危机的有效途径。研究采用Illumina高通量测序平台对西北地区典型池塘鱼菜共生养殖模式下, 养殖水环境与根系细菌微生物16S rRNA基因(V3—V5区)进行高通量测序。比较群落结构和微生物多样性表明, 细菌总有效序列为56444条, 细菌物种平均注释OTU数目为945条。在鱼菜共生模式下的优势门类为变形菌门(Proteobacteria)、拟杆菌门(Bacteroidetes)、厚壁菌门(Firmicutes)、蓝菌门(Cyanobacteria)。根系微生物群落中占优势地位的分别为拟杆菌门(Bacteroidetes)和变形菌门(Proteobacteria)。OTU聚类分析结果显示6个样品共有OTU数目为165, 其中新疆水生野生动物救护中心鱼菜共生池塘水样(JH.KX)特有OTU数量最高(137), 其空心菜根系水样(JH.R)最低(30), 而乌鲁木齐市米东区长山子镇水源地养殖基地空心菜根系水样(MD.R)和新疆水生野生动物救护中心对照池塘水样(JH.C)数量相同(85)。通过对各OTU进行注释, 发现了根瘤菌、鱼类致病菌、耐寒细菌、降解有机物细菌、降低重金属及毒素危害的细菌。研究为鱼菜共生生态养殖模式下微生物群落结构及其生态调控机理提供了参考依据。     

A metagenomic approach to characterization of the vaginal microbiome signature in pregnancy

While current major national research efforts (i.e., the NIH Human Microbiome Project) will enable comprehensive metagenomic characterization of the adult human microbiota, how and when these diverse microbial communities take up residence in the host and during reproductive life are unexplored at a population level. Because microbial abundance and diversity might differ in pregnancy, we sought to generate comparative metagenomic signatures across gestational age strata. DNA was isolated from the vagina (introitus, posterior fornix, midvagina) and the V5V3 region of bacterial 16S rRNA genes were sequenced (454FLX Titanium platform). Sixty-eight samples from 24 healthy gravidae (18 to 40 confirmed weeks) were compared with 301 non-pregnant controls (60 subjects). Generated sequence data were quality filtered, taxonomically binned, normalized, and organized by phylogeny and into operational taxonomic units (OTU); principal coordinates analysis (PCoA) of the resultant beta diversity measures were used for visualization and analysis in association with sample clinical metadata. Altogether, 1.4 gigabytes of data containing >2.5 million reads (averaging 6,837 sequences/sample of 493 nt in length) were generated for computational analyses. Although gravidae were not excluded by virtue of a posterior fornix pH >4.5 at the time of screening, unique vaginal microbiome signature encompassing several specific OTUs and higher-level clades was nevertheless observed and confirmed using a combination of phylogenetic, non-phylogenetic, supervised, and unsupervised approaches. Both overall diversity and richness were reduced in pregnancy, with dominance of Lactobacillus species (L. iners crispatus, jensenii and johnsonii, and the orders Lactobacillales (and Lactobacillaceae family), Clostridiales, Bacteroidales, and Actinomycetales. This intergroup comparison using rigorous standardized sampling protocols and analytical methodologies provides robust initial evidence that the vaginal microbial 16S rRNA gene catalogue uniquely differs in pregnancy, with variance of taxa across vaginal subsite and gestational age.     

Helminth Colonization Is Associated with Increased Diversity of the Gut Microbiota

Soil-transmitted helminths colonize more than 1.5 billion people worldwide, yet little is known about how they interact with bacterial communities in the gut microbiota. Differences in the gut microbiota between individuals living in developed and developing countries may be partly due to the presence of helminths, since they predominantly infect individuals from developing countries, such as the indigenous communities in Malaysia we examine in this work. We compared the composition and diversity of bacterial communities from the fecal microbiota of 51 people from two villages in Malaysia, of which 36 (70.6%) were infected by helminths. The 16S rRNA V4 region was sequenced at an average of nineteen thousand sequences per samples. Helminth-colonized individuals had greater species richness and number of observed OTUs with enrichment of Paraprevotellaceae, especially with Trichuris infection. We developed a new approach of combining centered log-ratio (clr) transformation for OTU relative abundances with sparse Partial Least Squares Discriminant Analysis (sPLS-DA) to enable more robust predictions of OTU interrelationships. These results suggest that helminths may have an impact on the diversity, bacterial community structure and function of the gut microbiota.     

Denaturing gradient gel electrophoresis analyses of the vertical distribution and diversity of Vibrio spp. populations in the Cariaco Basin

The Cariaco system is the second largest permanently anoxic marine water body in the world. Its water column is characterized by a pronounced vertical layering of microbial communities. The goal of our study was to investigate the vertical distribution and diversity of Vibrio spp. present in the Cariaco Basin waters using denaturing gradient gel electrophoresis (DGGE) of PCR-amplified 16S rDNA fragments. Representatives of the Vibrio genus were detected by nested and direct PCR in seawater at 10 depths. Sequence analyses of 55 DGGE bands revealed that only 11 different operational taxonomic units (OTU) are identified as Vibrio species. Between one and five OTUs were detected at each depth and the most common OTUs were OTU 1 and OTU 2, which phylogenetically clustered with Vibrio chagasii and Vibrio fortis, respectively. OTUs 3 and 4 were only found in the anoxic zone and were identified as Vibrio orientalis and Vibrio neptunius, respectively. Several Vibrio species detected are potentially pathogenic to human, prawns and corals such as Vibrio parahaemolyticus, Vibrio fischeri and Vibrio shilonii. In the Cariaco Basin, different Vibrio species were found to be specific to specific depths strata, suggesting that this genus is a natural component of the microbial communities in this marine redox environment.     

Variation in Mycorrhizal Associations with Tulasnelloid Fungi among Populations of Five Dactylorhiza Species

Background

Orchid species rely on mycorrhizal symbioses with fungi to complete their life cycle. Although there is mounting evidence that orchids can associate with several fungi from different clades or families, less is known about the actual geographic distribution of these fungi and how they are distributed across different orchid species within a genus.

Methodology/Principal Findings

We investigated among-population variation in mycorrhizal associations in five species of the genus Dactylorhiza (D. fuchsii, D. incarnata, D. maculata, D. majalis and D. praetermissa) using culture-independent detection and identification techniques enabling simultaneous detection of multiple fungi in a single individual. Mycorrhizal specificity, determined as the number of fungal operational taxonomic units (OTUs), and phylogenetic diversity of fungi were compared between species, whereas discriminant analysis was used to compare mycorrhizal spectra across populations and species. Based on a 95% cut-off value in internal transcribed spacer (ITS) sequence similarity, a total of ten OTUs was identified belonging to three different clades within the Tulasnellaceae. Most OTUs were found in two or more Dactylorhiza species, and some of them were common and widespread, occurring in more than 50% of all sampled populations. Each orchid species associated with at least five different OTUs, whereas most individuals also associated with two or more fungal OTUs at the same time. Phylogenetic diversity, corrected for species richness, was not significantly different between species, confirming the generality of the observed orchid mycorrhizal associations.

Conclusions/Significance

We found that the investigated species of the genus Dactylorhiza associated with a wide range of fungal OTUs from the Tulasnellaceae, some of which were widespread and common. These findings challenge the idea that orchid rarity is related to mycorrhizal specificity and fungal distribution.     

Fungal disease and temperature alter skin microbiome structure in an experimental salamander system

Pathogens compete with host microbiomes for space and resources. Their shared environment impacts pathogen–microbiome–host interactions, which can lead to variation in disease outcome. The skin microbiome of red‐backed salamanders (Plethodon cinereus) can reduce infection by the pathogen Batrachochytrium dendrobatidis (Bd) at moderate infection loads, with high species richness and high abundance of competitors as putative mechanisms. However, it is unclear if the skin microbiome can reduce epizootic Bd loads across temperatures. We conducted a laboratory experiment to quantify skin microbiome and host responses (P. cinereus: n = 87) to Bd at mimicked epizootic loads across temperatures (13, 17 and 21°C). We quantified skin microbiomes using 16S rRNA gene metabarcoding and identified operational taxonomic units (OTUs) taxonomically similar to culturable bacteria known to kill Bd (anti‐Bd OTUs). Prior to pathogen exposure, temperature changed the microbiome (OTU richness decreased by 12% and the abundance of anti‐Bd OTUs increased by 18% per degree increase in temperature), but these changes were not predictive of disease outcome. After exposure, Bd changed the microbiome (OTU richness decreased by 0.1% and the abundance of anti‐Bd OTUs increased by 0.2% per 1% increase in Bd load) and caused high host mortality across temperatures (35/45: 78%). Temperature indirectly impacted microbiome change and mortality through its direct effect on pathogen load. We did not find support for the microbiome impacting Bd load or host survival. Our research reveals complex host, pathogen, microbiome and environmental interactions to demonstrate that during epizootic events the microbiome will be unlikely to reduce pathogen invasion, even for putatively Bd‐resistant species.     

Patterns of community composition and diversity in latent fungi of living Quercus serrata trunks across a range of oak wilt prevalence and climate variables in Japan

Given that forest dieback due to emerging pests is increasing under global warming, understanding the relationships between pests, climate, and forest biodiversity is an urgent priority. In Japan, mass attacks of an ambrosia beetle, vectoring a pathogenic fungus, cause oak wilt outbreaks in recent decades. Here, the associations of oak wilt and climate with wood-inhabiting fungal communities in apparently healthy Quercus serrata trunks were investigated using DNA metabarcoding in seven sites along a climatic gradient in Japan. Amplicon sequencing of the fungal internal transcribed spacer 1 region generated 1,339,958 sequence reads containing 879 fungal operational taxonomic units (OTUs) in 234 wood samples. OTU compositions were significantly different between sites with and without oak wilt. OTU richness increased with temperature and precipitation at sites where oak wilt was present, but this relationship was not observed at sites without oak wilt, possibly due to interaction between oak wilt and climate.     

Locality or habitat? Exploring predictors of biodiversity in Amazonia

Amazonia is an environmentally heterogeneous and biologically megadiverse region, and its biodiversity varies considerably over space. However, existing knowledge on Amazonian biodiversity and its environmental determinants stems almost exclusively from studies of macroscopic above‐ground organisms, notably vertebrates and trees. In contrast, diversity patterns of most other organisms remain elusive, although some of them, for instance microorganisms, constitute the overwhelming majority of taxa in any given location, both in terms of diversity and abundance. Here, we use DNA metabarcoding to estimate prokaryote and eukaryote diversity in environmental soil and litter samples from 39 survey plots in a longitudinal transect across Brazilian Amazonia using 16S and 18S gene sequences, respectively. We characterize richness and community composition based on operational taxonomic units (OTUs) and test their correlation with longitude and habitat. We find that prokaryote and eukaryote OTU richness and community composition differ significantly among localities and habitats, and that prokaryotes are more strongly structured by locality and habitat type than eukaryotes. Our results 1) provide a first large‐scale mapping of Amazonian soil biodiversity, suggesting that OTU richness patterns might follow substantially different patterns from those observed for macro‐organisms; and 2) indicate that locality and habitat factors interact in determining OTU richness patterns and community composition. This study shows the potential of DNA metabarcoding in unveiling Amazonia's outstanding diversity, despite the lack of complete reference sequence databases for the organisms sequenced.     

Characterization and identification of productivity-associated rhizobacteria in wheat

The rhizosphere is populated by a numerous and diverse array of rhizobacteria, and many impact productivity in largely unknown ways. Here we characterize the rhizobacterial community in a wheat variety categorized according to shoot biomass using 16S rRNA pyrosequencing abundance data. Plants were grown in homogenized field soil under greenhouse conditions, and DNA was extracted and pyrosequenced, resulting in 29,007 quality sequences. Operational taxonomic units (OTUs) that were significantly associated with biomass productivity were identified using an exact test adjusted for the false-discovery rate. The productivity deviation expressed as a percentage of the total mean square for regression (PMSR) was determined for each OTU. Out of 719 OTUs, 42 showed significant positive associations and 39 showed significant negative associations (q value, ≤0.05). OTUs with the greatest net positive associations, by genus, were as follows: Duganella, OTU 43 and OTU 3; Janthinobacterium, OTU 278; Pseudomonas, OTU 588; and Cellvibrio, OTU 1847. Those with negative associations were as follows: Bacteria, OTU 273; Chryseobacterium, OTU 508; Proteobacteria, OTU 249; and Enterobacter, OTU 357. Shoot biomass productivity was strongly correlated with the balance between the overall abundances of positive- and negative-productivity-associated OTUs. High-productivity rhizospheres contained 9.2 significant positives for every negatively associated rhizobacterium, while low-productivity rhizospheres showed 2.3 significant negatives for every positively associated rhizobacterium. Overall rhizobacterial community diversity as measured by the Chao1, Shannon, and Simpson indexes was nonlinearly related to productivity, closely fitting a wavelike cubic equation. We conclude that shoot biomass productivity is strongly related to the ratio of positive- to negative-productivity-associated rhizobacteria in the rhizosphere. This study identifies significant OTUs composing the productive and unproductive rhizobacterial communities.     

Polymerase matters: non-proofreading enzymes inflate fungal community richness estimates by up to 15 %

Rare taxa overwhelm metabarcoding data generated using next-generation sequencing (NGS). Low frequency Operational Taxonomic Units (OTUs) may be artifacts generated by PCR-amplification errors resulting from polymerase mispairing. We analyzed two Internal Transcribed Spacer 2 (ITS2) MiSeq libraries generated with proofreading (ThermoScientific Phusion^®) and non-proofreading (ThermoScientific Phire^®) polymerases from the same MiSeq reaction, the same samples, using the same DNA tags, and with two different clustering methods to evaluate the effect of polymerase and clustering tool choices on the estimates of richness, diversity and community composition. Our data show that, while the overall communities are comparable, OTU richness is exaggerated by the use of the non-proofreading polymerase–up to 15 % depending on the clustering method, and on the threshold of low frequency OTU removal. The overestimation of richness also consistently led to underestimation of community evenness, a result of increase in the low frequency OTUs. Stringent thresholds of eliminating the rare reads remedy this issue; exclusion of reads that occurred ≤10 times reduced overestimated OTU numbers to <0.3 %. As a result of these findings, we strongly recommend the use of proofreading polymerases to improve the data integrity as well as the use of stringent culling thresholds for rare sequences to minimize overestimation of community richness.     

Dynamic bacterial communities on reverse-osmosis membranes in a full-scale desalination plant

To better understand biofouling of seawater reverse osmosis (SWRO) membranes, bacterial diversity was characterized in the intake water, in subsequently pretreated water and on SWRO membranes from a full-scale desalination plant (FSDP) during a 9 month period. 16S rRNA gene fingerprinting and sequencing revealed that bacterial communities in the water samples and on the SWRO membranes were very different. For the different sampling dates, the bacterial diversity of the active and the total bacterial fractions of the water samples remained relatively stable over the sampling period whereas the bacterial community structure on the four SWRO membrane samples was significantly different. The richness and evenness of the SWRO membrane bacterial communities increased with usage time with an increase in the Shannon diversity index of 2.2 to 3.7. In the oldest SWRO membrane (330 days), no single operational taxonomic unit (OTU) dominated and the majority of the OTUs fell into the Alphaproteobacteria or the Planctomycetes. In striking contrast, a Betaproteobacteria OTU affiliated to the genus Ideonella was dominant and exclusively found in the membrane used for the shortest time (10 days). This suggests that bacteria belonging to this genus could be one of the primary colonizers of the SWRO membrane. Knowledge of the dominant bacterial species on SWRO membranes and their dynamics should help guide culture studies for physiological characterization of biofilm forming species.     

Ironing out the wrinkles in the rare biosphere through improved OTU clustering

Deep sequencing of PCR amplicon libraries facilitates the detection of low‐abundance populations in environmental DNA surveys of complex microbial communities. At the same time, deep sequencing can lead to overestimates of microbial diversity through the generation of low‐frequency, error‐prone reads. Even with sequencing error rates below 0.005 per nucleotide position, the common method of generating operational taxonomic units (OTUs) by multiple sequence alignment and complete‐linkage clustering significantly increases the number of predicted OTUs and inflates richness estimates. We show that a 2% single‐linkage preclustering methodology followed by an average‐linkage clustering based on pairwise alignments more accurately predicts expected OTUs in both single and pooled template preparations of known taxonomic composition. This new clustering method can reduce the OTU richness in environmental samples by as much as 30–60% but does not reduce the fraction of OTUs in long‐tailed rank abundance curves that defines the rare biosphere.     

Many broadly-shared mycobionts characterize mycorrhizal interactions of two coexisting epiphytic orchids in a high elevation tropical forest

Interspecific interactions play an important role in community assembly. A basic ecological question is whether interactions are specialized (one to one) or generalized (many to many). Specialization of interactions should ideally be assessed across several populations because species could be specialists at a particular site but generalists when several sites are considered. Mycorrhizal interactions are fundamental for orchid life and distribution, but their level of specialization is still under debate. To understand the extent to which epiphytic orchids are specialists in their mycorrhizal interactions, we studied the richness and phylogenetic structure of mycobionts across different sites, and the similarity in the mycobiont composition between coexisting orchid species. We sequenced the nrDNA ITS2 region and explored the mycobiont communities associated with two epiphytic orchids, Epidendrum marsupiale and Cyrtochilum pardinum, at two elevations within two sites in Ecuador. We found 108 OTUs belonging to Serendipitaceae (66), Ceratobasidiaceae (22), Atractiellales (11) and Tulasnellaceae (9). Orchids at the highest elevations hosted the highest OTU richness. The two orchid species shared a high percentage of mycobionts between all sites. No phylogenetic structure within orchid mycorrhizal communities was found at any sites or elevations. Our results indicate that the studied orchids are generalists and share a broad group of mycobionts (16 OTUs) with no apparent niche segregation within or between sites.