Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.     

Investigating Microbial Eukaryotic Diversity from a Global Census: Insights from a Comparison of Pyrotag and Full-Length Sequences of 18S rRNA Genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.     

Synthetic Statistical Approach Reveals a High Degree of Richness of Microbial Eukaryotes in an Anoxic Water Column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.     

Low-diversity fungal assemblage in an Antarctic Dry Valleys soil

The McMurdo Dry Valleys of Antarctica present extreme environmental challenges. Life is restricted to patchy occurrence of lichens, mosses and invertebrates, plus soil microbial communities. Fungi have been described in lichen symbioses but relatively little is known about the occurrence of free-living soil fungi in the Dry Valleys. A challenge in estimating fungal species richness has been the extent to which estimates based on either cultivation or environmental DNA reflect the active assemblage in cold-arid soils. Here, we describe analysis for inland Dry Valleys soil of environmental DNA and RNA (cDNA) to infer total and putative metabolically active assemblages, respectively, plus cultivation approaches using a variety of laboratory growth conditions. Environmental sequences indicated a highly restricted assemblage of just seven phylotypes that affiliated phylogenetically within two known genera, Helicodendron and Zalerion, plus previously unidentified fungal phylotypes. None of the commonly encountered molds and mitosporic genera recorded from maritime Antarctic locations were encountered. A striking difference was observed in the frequency of recovery for phylotypes between libraries. This suggests that both species richness and beta diversity estimates based on DNA libraries have the potential to misinform putatively active assemblages. Cultivation yielded a cold-tolerant Zalerion strain that affiliated with DNA and RNA library clones, and a psychrotrophic yeast (Debaryomyces hansenii), which was not detected using either culture-independent approach.     

Synthetic statistical approach reveals a high degree of richness of microbial eukaryotes in an anoxic water column

Molecular surveys suggest that communities of microbial eukaryotes are remarkably rich, because even large clone libraries seem to capture only a minority of species. This provides a qualitative picture of protistan richness but does not measure its real extent either locally or globally. Statistical analysis can estimate a community's richness, but the specific methods used to date are not always well grounded in statistical theory. Here we study a large protistan molecular survey from an anoxic water column in the Cariaco Basin (Caribbean Sea). We group individual 18S rRNA gene sequences into operational taxonomic units (OTUs) using different cutoff values for sequence similarity (99 to 50%) and systematically apply parametric models and nonparametric estimators to the OTU frequency data to estimate the total protistan diversity. The parametric models provided statistically sound estimates of protistan richness, with biologically meaningful standard errors, maximal data usage, and extensive model diagnostics and were preferable to the available nonparametric tools. Our clone library exceeded 700 clones but still covered only a minority of species and less than half of the larger protistan clades. Our estimates of total protistan richness portray the target community as very rich at all OTU levels, with hundreds of different populations apparently co-occurring in the small (3-liter) volume of our sample, as well as dozens of clades of the highest taxonomic order. These estimates are among the first for microbial eukaryotes that are obtained using state-of-the-art statistical methods and can serve as benchmark numbers for the local diversity of protists.     

Molecular screening of free-living microbial eukaryotes: diversity and distribution using a meta-analysis

Numerous environmental gene library studies have shown that eukaryote microbial diversity is much greater than expected. Molecular surveys of several 'extreme' and some more anthropomorphically commonplace environments have revealed many previously unsampled micro-eukaryotic lineages. However, it cannot be assumed that all of the sequences recovered from these studies are derived from real organisms, and for those that are, many questions remain about their distribution and ecology. Integrating all available sequence data from these studies reveals patterns of distribution, diversity and evolutionary relationships that are not accessible from independent analyses of the individual surveys and enables us to review the wider implications of such studies.     

Amplification by PCR artificially reduces the proportion of the rare biosphere in microbial communities

The microbial world has been shown to hold an unimaginable diversity. The use of rRNA genes and PCR amplification to assess microbial community structure and diversity present biases that need to be analyzed in order to understand the risks involved in those estimates. Herein, we show that PCR amplification of specific sequence targets within a community depends on the fractions that those sequences represent to the total DNA template. Using quantitative, real-time, multiplex PCR and specific Taqman probes, the amplification of 16S rRNA genes from four bacterial species within a laboratory community were monitored. Results indicate that the relative amplification efficiency for each bacterial species is a nonlinear function of the fraction that each of those taxa represent within a community or multispecies DNA template. Consequently, the low-proportion taxa in a community are under-represented during PCR-based surveys and a large number of sequences might need to be processed to detect some of the bacterial taxa within the 'rare biosphere'. The structure of microbial communities from PCR-based surveys is clearly biased against low abundant taxa which are required to decipher the complete extent of microbial diversity in nature.     

南海陆坡沉积物细菌丰度预测

分别利用参数模型和无参数估计法预测南海陆坡沉积物柱MD05-2896中的细菌丰度.基于非培养的PCR-RFLP的16SrRNA基因分子技术,扩增了沉积物柱中的细菌16S rRNA基因序列,并构建16S rRNA基因文库.系统发育分析表明16S rRNA基因文库中,大多数序列属于17个已知的“门”.分别以99％、97％、90％和80％序列一致性作为分类单元分界点,将16SrRNA基因序列组群为分类单元.使用逆高斯分布模型、对数正态分布模型、负二项式分布模型、帕雷托分布模型、双指数分布模型以及ACE、ACE-1等估计方法预测不同分类单元分类水平下的细菌丰度.结果表明在“种”级分类水平上,负二项式分布为最优估计模型,估计细菌丰度为244±10(SE).不过,受实验条件的限制,该估计值可能偏低.     

Molecular evidence for widespread occurrence of Foraminifera in soils

Environmental SSU rDNA‐based surveys are contributing to the dramatic revision of eukaryotic high‐level diversity and phylogeny as the number of sequence data increases. This ongoing revolution gives the opportunity to test for the presence of some eukaryotic taxa in environments where they have not been found using classical microscopic observations. Here, we test whether the foraminifera, a group of single‐celled eukaryotes, considered generally as typical for the marine ecosystems are present in soil. We performed foraminiferal‐specific nested PCR on 20 soil DNA samples collected in contrasted environments. Unexpectedly, we found that the majority of the samples contain foraminiferal SSU rDNA sequences. In total, we obtained 49 sequences from 17 localities. Phylogenetic analysis clusters them in four groups branching among the radiation of early foraminiferal lineages. Three of these groups also include sequences originated from previous freshwater surveys, suggesting that there were up to four independent colonization events of terrestrial and/or freshwater ecosystems by ancestral foraminifera. As shown by our data, foraminifera are a widespread and diverse component of soil microbial communities. Yet, identification of terrestrial foraminiferal species and understanding of their ecological role represent an exciting challenge for future research.     

The generation and maintenance of diversity in microbial communities

Microorganisms play a central role in the regulation of ecosystem processes, and they comprise the vast majority of species on Earth. With recent developments in molecular methods, it has become tractable to quantify the extent of microbial diversity in natural environments. Here we examine this revolution in our understanding of microbial diversity, and we explore the factors that contribute to the seemingly astounding numbers of microbial taxa found within individual environmental samples. We conducted a meta-analysis of bacterial richness estimates from a variety of ecosystems. Nearly all environments contained hundreds to thousands of bacterial taxa, and richness levels increased with the number of individuals in a sample, a pattern consistent with those reported for nonmicrobial taxa. A cursory comparison might suggest that bacterial richness far exceeds the richness levels typically observed for plant and animal taxa. However, the apparent diversity of bacterial communities is influenced by phylogenetic breadth and allometric scaling issues. When these features are taken into consideration, the levels of microbial diversity may appear less astounding. Although the fields of ecology and biogeography have traditionally ignored microorganisms, there are no longer valid excuses for neglecting microorganisms in surveys of biodiversity. Many of the concepts developed to explain plant and animal diversity patterns can also be applied to microorganisms once we reconcile the scale of our analyses to the scale of the organisms being observed. Furthermore, knowledge from microbial systems may provide insight into the mechanisms that generate and maintain species richness in nonmicrobial systems.     

Effects of OTU Clustering and PCR Artifacts on Microbial Diversity Estimates

Next-generation sequencing has increased the coverage of microbial diversity surveys by orders of magnitude, but differentiating artifacts from rare environmental sequences remains a challenge. Clustering 16S rRNA sequences into operational taxonomic units (OTUs) organizes sequence data into groups of 97 % identity, helping to reduce data volumes and avoid analyzing sequencing artifacts by grouping them with real sequences. Here, we analyze sequence abundance distributions across environmental samples and show that 16S rRNA sequences of >99 % identity can represent functionally distinct microorganisms, rendering OTU clustering problematic when the goal is an accurate analysis of organism distribution. Strict postsequencing quality control (QC) filters eliminated the most prevalent artifacts without clustering. Further experiments proved that DNA polymerase errors in polymerase chain reaction (PCR) generate a significant number of substitution errors, most of which pass QC filters. Based on our findings, we recommend minimizing the number of PCR cycles in DNA library preparation and applying strict postsequencing QC filters to reduce the most prevalent artifacts while maintaining a high level of accuracy in diversity estimates. We further recommend correlating rare and abundant sequences across environmental samples, rather than clustering into OTUs, to identify remaining sequence artifacts without losing the resolution afforded by high-throughput sequencing.     

Mitochondrial cox1 sequence data reliably uncover patterns of insect diversity but suffer from high lineage-idiosyncratic error rates

Background

The demand for scientific biodiversity data is increasing, but taxonomic expertise is often limited or not available. DNA sequencing is a potential remedy to overcome this taxonomic impediment. Mitochondrial DNA is most commonly used, e.g., for species identification (“DNA barcoding”). Here, we present the first study in arthropods based on a near-complete species sampling of a family-level taxon from the entire Australian region. We aimed to assess how reliably mtDNA data can capture species diversity when many sister species pairs are included. Then, we contrasted phylogenetic subsampling with the hitherto more commonly applied geographical subsampling, where sister species are not necessarily captured.

Methodology/Principal Findings

We sequenced 800 bp cox1 for 1,439 individuals including 260 Australian species (78% species coverage). We used clustering with thresholds of 1 to 10% and general mixed Yule Coalescent (GMYC) analysis for the estimation of species richness. The performance metrics used were taxonomic accuracy and agreement between the morphological and molecular species richness estimation. Clustering (at the 3% level) and GMYC reliably estimated species diversity for single or multiple geographic regions, with an error for larger clades of lower than 10%, thus outperforming parataxonomy. However, the rates of error were higher for some individual genera, with values of up to 45% when very recent species formed nonmonophyletic clusters. Taxonomic accuracy was always lower, with error rates above 20% and a larger variation at the genus level (0 to 70%). Sørensen similarity indices calculated for morphospecies, 3% clusters and GMYC entities for different pairs of localities was consistent among methods and showed expected decrease over distance.

Conclusion/Significance

Cox1 sequence data are a powerful tool for large-scale species richness estimation, with a great potential for use in ecology and β-diversity studies and for setting conservation priorities. However, error rates can be high in individual lineages.     

HOST SWITCHING PROMOTES DIVERSITY IN HOST‐SPECIALIZED MYCOPARASITIC FUNGI: UNCOUPLED EVOLUTION IN THE BIATOROPSIS‐USNEA SYSTEM

Fungal mycoparasitism—fungi parasitizing other fungi—is a common lifestyle in some basal lineages of the basidiomycetes, particularly within the Tremellales. Relatively nonaggressive mycoparasitic fungi of this group are in general highly host specific, suggesting cospeciation as a plausible speciation mode in these associations. Species delimitation in the Tremellales is often challenging because morphological characters are scant. Host specificity is therefore a great aid to discriminate between species but appropriate species delimitation methods that account for actual diversity are needed to identify both specialist and generalist taxa and avoid inflating or underestimating diversity. We use the Biatoropsis‐Usnea system to study factors inducing parasite diversification. We employ morphological, ecological, and molecular data—methods including genealogical concordance phylogenetic species recognition (GCPSR) and the general mixed Yule‐coalescent (GMYC) model—to assess the diversity of fungi currently assigned to Biatoropsis usnearum. The degree of cospeciation in this association is assessed with two cophylogeny analysis tools (ParaFit and Jane 4.0). Biatoropsis constitutes a species complex formed by at least seven different independent lineages and host switching is a prominent force driving speciation, particularly in host specialists. Combining ITS and nLSU is recommended as barcode system in tremellalean fungi.     

Unexpected cryptic species diversity of parasites of the family Xenidae (Strepsiptera) with a constant diversification rate over time

Parasitism is one of the most successful and ancient strategies. Due to the specialized lifestyle of parasites, they are usually affected by reductions and changes in their body plan in comparison with nonparasitic sister groups. Extreme environmental conditions may impose restraints on behavioural or physiological adaptations to a specific host and limit morphological changes associated with speciation. Such morphological homogeneity has led to the diversity of parasites being underestimated in morphological studies. By contrast, the species concept has dramatically changed in many parasitic groups during recent decades of study using DNA sequence data. Here we tested the phenomenon of cryptic species diversity in the twisted‐wing parasite family Xenidae (Strepsiptera) using nuclear and mitochondrial DNA sequence data for a broad sample of Xenidae. We used three quantitative methods of species delimitation from the molecular phylogenetic data – one distance‐based (ABGD) and two tree‐based (GMYC, bPTP). We found 77–96 putative species in our data and suggested the number of Xenidae species to be more diverse than expected. We identified 67 hosts to species level and almost half of them were not previously known as hosts of Xenidae. The mean number of host species per putative species varied between 1.39 and 1.55. The constant rate in net diversification can be explained by the flexibility of this parasitic group, represented by their ability to colonize new host lineages combined with passive long‐range dispersal by hosts.     

Emerging diversity within chrysophytes, choanoflagellates and bicosoecids based on molecular surveys

In recent years, a substantial amount of data on aquatic protists has been obtained from culture-independent molecular approaches, unveiling a large diversity and the existence of new lineages. However, sequences affiliated with minor groups (in terms of clonal abundance) have often been under-analyzed, and this hides a potentially relevant source of phylogenetic information. Here we have searched public databases for 18S rDNA sequences of chrysophytes, choanoflagellates and bicosoecids retrieved from molecular surveys of protists. These three groups are often considered to account for most of the heterotrophic flagellates, an important functional component in microbial food webs. They represented a significant fraction of clones in freshwater studies, whereas their relative clonal abundance was low in marine studies. The novelty displayed by this dataset was notable. Most environmental sequences were distant to sequences of cultured organisms, indicating a significant bias in the representation of taxa in culture. Moreover, they were often distant to sequences from other molecular surveys, suggesting an insufficient sequencing effort to characterize the in situ diversity of these groups. Phylogenetic trees with complete sequences present the most accurate representation of the diversity of these groups, with the emergence of several new clades formed exclusively by environmental sequences. Exhaustive data mining in sequence databases allowed the identification of new diversity hidden inside chrysophytes, choanoflagellates and bicosoecids.     

Long‐term development of species richness in a central European beech (Fagus sylvatica) forest affected by windthrow—Support for the intermediate disturbance hypothesis?

On the basis of long‐term surveys of permanent plots and traps, we examined the communities of saproxylic beetles, fungi, herbs, and trees on an untreated 22 ha large beech forest windthrow and asked whether the results lend support to the intermediate disturbance hypothesis (IDH). We studied species richness and the similarity of community composition. Additionally, we grouped species by their frequency trend over time to successional model types to examine whether, corresponding to the IDH, the diversity of these groups explained peak richness at intermediate intervals after the disturbance. In line with the IDH, species richness showed a hump‐backed temporal course for alpha and gamma diversity. We found evidence for a linear succession directly after the disturbance. This, however, did not continue, and in all species groups, a partial recovery of the initial community was observed. In the case of fungi, herbs, and trees, but not for saproxylic beetles, alpha diversity was driven by the diversity of the successional model types. Our results underline that the mechanisms driving species richness after disturbances are more complex than the IDH suggests and that these mechanisms vary with species group. We assumed that, besides competition, legacy effects, facilitation, habitat heterogeneity, and random saturation of the species pool are important. In case of trees and herbs, we found indications for strong legacy and competition effects. For fungi and beetles, substrate heterogeneity and microclimate were assumed to be important. We concluded that disturbances contribute to increasing species richness not only by reducing the effectiveness of competitors but also by increasing the amount and diversity of resources, as well as their rate of change over time.     

Performance of the COX1 gene as a marker for the study of metabolically active Pezizomycotina and Agaricomycetes fungal communities from the analysis of soil RNA

In temperate forest soils, filamentous ectomycorrhizal and saprotrophic fungi affiliated to the Agaricomycetes and Pezizomycotina contribute to key biological processes. The diversity of soil fungal communities is usually estimated by studying molecular markers such as nuclear ribosomal gene regions amplified from soil-extracted DNA. However, this approach only reveals the presence of the corresponding genomic DNA in the soil sample and may not reflect the diversity of the metabolically active species. To circumvent this problem, we investigated the performance of the mitochondrial cytochrome c oxidase 1 (COX1)-encoding gene as a fungal molecular marker for environmental RNA-based studies. We designed PCR primers to specifically amplify Agaricomycetes and Pezizomycotina COX1 partial sequences and amplified them from both soil DNA and reverse-transcribed soil RNA. As a control, we also amplified the nuclear internal transcribed spacer ribosomal region from soil DNA. Fungal COX1 sequences were readily amplified from soil-extracted nucleic acids and were not significantly contaminated by nontarget sequences. We show that the relative abundance of fungal taxonomic groups differed between the different sequence data sets, with for example ascomycete COX1 sequences being more abundant among sequences amplified from soil DNA than from soil cDNAs.     

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.     

Phylogenetic diversity of macromycetes and woody plants along an elevational gradient in Eastern Mexico

Phylogenetic information provides insight into the ecological and evolutionary processes that organize species assemblages. We compared patterns of phylogenetic diversity among macromycete and woody plant communities along a steep elevational gradient in eastern Mexico to better understand the evolutionary processes that structure their communities. Macrofungi and trees were counted and identified in eight sites from 100 to 3500 m asl, and sequence data retrieved from GenBank for the same or closely related species were used to reconstruct their phylogenies. Patterns of species richness and phylogenetic diversity were similar for both macrofungi and trees, but macromycete richness and diversity peaked at mid‐elevations, whereas woody plant richness and diversity did not show significant trends with elevation. Phylogenetic similarity among sites was low for both groups and decreased as elevational distance between sites increased. Macromycete communities displayed phylogenetic overdispersion at low elevations and phylogenetic clustering at high elevations; the latter is consistent with environmental filtering at high elevation sites. Woody plants generally exhibited phylogenetic clustering, consistent with the potential importance of environmental filtering throughout the elevational gradient.