Amplification of soil fungal community DNA using the ITS86F and ITS4 primers

Internal transcribed spacer (ITS) 86F and ITS4 and the ITS1-F and ITS86R primer pairs were tested to specifically amplify fungal community DNA extracted from soil. Libraries were constructed from PCR-amplified fragments, sequenced and compared against sequences deposited in GenBank. The results confirmed that the ITS86F and ITS4 primer pair was selectively specific for the Ascomycetes, Basidiomycetes and Zygomycetes fungal clades. Amplified products generated by the ITS1F and ITS86R primer pair also aligned with sequences from a range of species within the Ascomycete and Basidiomycete clades but not from the Zygomycete. Both primer sets demonstrated fungal specificity and appear to be well suited for rapid PCR-based (fingerprinting) analysis of environmental fungal community DNA. This is the first reported use and assessment of the ITS86F and ITS4 and the ITS1-F and ITS86R primer pairs in amplifying fungal community DNA from soil.     

Comparison of ITS1 and ITS2 rDNA in 454 sequencing of hyperdiverse fungal communities

Next-generation amplicon sequencing is a powerful tool in ecological studies of fungi. Technological development suggests that short fragment high-throughput techniques, e.g. Illumina, will gain importance in fungal community analyses. Thus there is a need for short (250 bp) and informative molecular identifiers. Here we compared ITS1 vs. ITS2 rDNA using empirical data from a study of hyperdiverse leaf-associated fungal communities. Our results suggest that ITS2 may be more variable and recovers more of the molecular diversity. We confirm an earlier in silico study showing that ITS1 and ITS2 yielded somewhat different taxonomic community compositions when blasted against public databases. However, we demonstrate that both ITS1 and ITS2 reveal similar patterns in community structure when analyzed in a community ecology context.     

Fungal diversity in sheep (Ovis aries) and cattle (Bos taurus) feces assessed by comparison of 18S, 28S and ITS ribosomal regions

To explore the fungal diversity in ruminant feces for bioenergy, libraries based on internal transcribed spacer (ITS), 18S rRNA, and 28S rRNA regions were constructed, respectively. Although the libraries were constructed from the same DNA extracts, the fungal taxa analyses based on these libraries are different. The ITS and 28S libraries comprised higher proportions of fungal clones than 18S libraries, and the ITS libraries converged into the lower diversities. The ITS libraries could be used to analyze the fungal community. The 18S libraries were suitable for the fungi and protozoa community. However, the 28S are suitable for analysis of Ascomycota fungi. The major fungal taxa in cattle feces analyzed by ITS, 18S, and 28S libraries are similar to those of sheep feces, respectively. The fungal taxa detected by the ITS library comprised only 20 % fungal taxa detected by the three libraries. The 18S library comprised 30 % fungal taxa; the 28S library comprised about 50 % fungal taxa. The results indicated that primer sets toward different DNA regions lead to the difference in structures of fungal community. So the selection of primer sets may influence the fungal communities, and libraries based on single primer sets may underestimate the fungal diversity. The comparison of ITS, 18S, and 28S libraries could fid more diverse fungi than that based on only one library.     

The rpb2 gene represents a viable alternative molecular marker for the analysis of environmental fungal communities

Although the commonly used internal transcribed spacer region of rDNA (ITS) is well suited for taxonomic identification of fungi, the information on the relative abundance of taxa and diversity is negatively affected by the multicopy nature of rDNA and the existence of ITS paralogues. Moreover, due to high variability, ITS sequences cannot be used for phylogenetic analyses of unrelated taxa. The part of single‐copy gene encoding the second largest subunit of RNA polymerase II (rpb2) was thus compared with first spacer of ITS as an alternative marker for the analysis of fungal communities in spruce forest topsoil, and their applicability was tested on a comprehensive mock community. In soil, rpb2 exhibited broad taxonomic coverage of the entire fungal tree of life including basal fungal lineages. The gene exhibited sufficient variation for the use in phylogenetic analyses and taxonomic assignments, although it amplifies also paralogues. The fungal taxon spectra obtained with rbp2 region and ITS1 corresponded, but sequence abundance differed widely, especially in the basal lineages. The proportions of OTU counts and read counts of major fungal groups were close to the reality when rpb2 was used as a molecular marker while they were strongly biased towards the Basidiomycota when using the ITS primers ITS1/ITS4. Although the taxonomic placement of rbp2 sequences is currently more difficult than that of the ITS sequences, its discriminative power, quantitative representation of community composition and suitability for phylogenetic analyses represent significant advantages.     

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.     

Intra-specific and intra-sporocarp ITS variation of ectomycorrhizal fungi as assessed by rDNA sequencing of sporocarps and pooled ectomycorrhizal roots from a <Emphasis Type="Italic">Quercus</Emphasis> woodland

The Internal Transcribed Spacer (ITS) regions of ribosomal DNA are widely used as markers for phylogenetic analyses and environmental sampling from a variety of organisms including fungi, plants, and animals. In theory, concerted evolution homogenizes multicopy genes so that little or no variation exists within populations or individuals. However, contrary to theory, ITS variation has been confirmed in populations and individuals from a diverse range of eukaryotes. The presence of intraspecific and intra-individual variation in multicopy genes has important implications for ecological and phylogenetic studies, yet relatively little is known about natural variation of these genes, particularly at the community level. In this study, we examined intraspecific and intra-sporocarp ITS variation by DNA sequencing from sporocarps and pooled roots from 68 species of ectomycorrhizal fungi collected at a single site in a Quercus woodland. We detected ITS variation in 27 species, roughly 40% of the taxa examined. Although intraspecific ITS variation was generally low (0.16–2.85%, mean = 0.74%), it was widespread within this fungal community. We detected ITS variation in both sporocarps and ectomycorrhizal roots, and variation was present within species of Ascomycota and Basidiomycota, two distantly related lineages within the Fungi. We discuss the implications of such widespread ITS variability with special reference to DNA-based environmental sampling from diverse fungal communities. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Soil chemical properties,rather than landscape context,influence woodland fungal communities along an urban‐rural gradient

With the expansion of cities around the world there is a growing interest in the factors that influence biodiversity and ecosystem processes in urban areas. Fungi are exceptionally diverse and play key roles in ecosystem function, yet despite predictions of negative impacts due to urbanization, fungi have been generally overlooked in urban ecological studies. We surveyed fungi in 16 remnant river red gum (Eucalyptus camaldulensis: Myrtaceae) woodlands along a gradient of 4–35 km from the city of Melbourne (south‐east Australia). Using both sporocarp surveys and terminal restriction fragment length polymorphism (T‐RFLP; primer pair ITS1‐F‐ITS4), we examined relationships between fungal community composition, landscape context (i.e. urbanization) and soil physicochemical properties. Community compositions from sporocarp data were significantly correlated with those from T‐RFLP data, largely because of correlations with ectomycorrhizal sporocarps (Spearman rank correlation coefficients ρ 0.31–0.42) rather than saprotrophic fungi (ρ 0.18–0.21). Principal components analysis of soil properties and non‐metric multidimensional scaling ordinations of fungal community composition showed no clear separation of sites according to urbanization, and there were no significant correlations between fungal community composition and urbanization. However, fungal community composition was significantly correlated with soil chemical properties (ρ 0.41–0.55). These data suggest that site‐scale soil properties, and associated effects of past and current land management activities, were more important in determining fungal community composition than the landscape‐level influences of urbanization.     

Exploring the accuracy of amplicon‐based internal transcribed spacer markers for a fungal community

With the continual improvement in high‐throughput sequencing technology and constant updates to fungal reference databases, the use of amplicon‐based DNA markers as a tool to reveal fungal diversity and composition in various ecosystems has become feasible. However, both primer selection and the experimental procedure require meticulous verification. Here, we computationally and experimentally evaluated the accuracy and specificity of three widely used or newly designed internal transcribed spacer (ITS) primer sets (ITS1F/ITS2, gITS7/ITS4 and 5.8S‐Fun/ITS4‐Fun). In silico evaluation revealed that primer coverage varied at different taxonomic levels due to differences in degeneracy and the location of primer sets. Using even and staggered mock community standards, we identified different proportions of chimeric and mismatch reads generated by different primer sets, as well as great variation in species abundances, suggesting that primer selection would affect the results of amplicon‐based metabarcoding studies. Choosing proofreading and high‐fidelity polymerase (KAPA HiFi) could significantly reduce the percentage of chimeric and mismatch sequences, further reducing inflation of operational taxonomic units. Moreover, for two types of environmental fungal communities, plant endophytic and soil fungi, it was demonstrated that the three primer sets could not reach a consensus on fungal community composition or diversity, and that primer selection, not experimental treatment, determines observed soil fungal community diversity and composition. Future DNA marker surveys should pay greater attention to potential primer effects and improve the experimental scheme to increase credibility and accuracy.     

New insights into the fungal community from the raw genomic sequence data of fig wasp Ceratosolen solmsi

Background

To date, biologists have discovered a large amount of valuable information from assembled genomes, but the abundant microbial data that is hidden in the raw genomic sequence data of plants and animals is usually ignored. In this study, the richness and composition of fungal community were determined in the raw genomic sequence data of Ceratosolen solmsi (RGSD-CS).

Results

To avoid the interference from sequences of C. solmsi, the unmapped raw data (about 17.1%) was obtained by excluding the assembled genome of C. solmsi from RGSD-CS. Comparing two fungal reference datasets, internal transcribed spacer (ITS) and large ribosomal subunit (LSU) of rRNA, the ITS dataset discovered a more diverse fungal community and was therefore selected as the reference dataset for evaluating the fungal community based on the unmapped raw data. The threshold of 95% sequence identity revealed many more matched fungal reads and fungal richness in the unmapped raw data than those by identities above 95%. Based on the threshold of 95% sequence identity, the fungal community of RGSD-CS was primarily composed of Saccharomycetes (88.4%) and two other classes (Agaricomycetes and Sordariomycetes, 8.3% in total). Compared with the fungal community of other reported fig wasps, Agaricomycetes and Eurotiomycetes were found to be unique to C. solmsi. In addition, the ratio of total fungal reads to RGSD-CS was estimated to be at least 4.8 × 10⁻³, which indicated that a large amount of fungal data was contained in RGSD-CS. However, rarefaction measure indicated that a deeper sequencing coverage with RGSD-CS was required to discover the entire fungal community of C. solmsi.

Conclusion

This study investigated the richness and composition of fungal community in RGSD-CS and provided new insights into the efficient study of microbial diversity using raw genomic sequence data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0370-3) contains supplementary material, which is available to authorized users.     

Diversity and composition of ectomycorrhizal community on seedling roots: the role of host preference and soil origin

As the main source of inocula, ectomycorrhizal (ECM) fungal propagules are critical for root colonization and seedling survival in deforested areas. It is essential to know factors that may affect the diversity and composition of ECM fungal community on roots of seedlings planted in deforest areas during reforestation. We quantitatively evaluated the effect of host plant and soil origin on ECM fungal propagule community structure established on roots of Castanopsis fargesii, Lithocarpus harlandii, Pinus armandii, and Pinus massoniana growing in soils from local natural forests and from sites deforested by clear-cut logging in the 1950s and 1960s. ECM root tips were sampled in April, July, and October of 2006, and ECM fungal communities were determined using ECM root morphotyping, internal transcribed spacer (ITS)-RFLP, and ITS sequencing. A total of 36 ECM fungal species were observed in our study, and species richness varied with host species and soil origin. Decreased colonization rates were found in all host species except for L. harlandii, and reduced species richness was found in all host species except for P. armandii in soil from the deforested site, which implied the great changes in ECM fungal community composition. Our results showed that 33.3% variance in ECM fungal community composition could be explained by host plant species and 4.6% by soil origin. Results of indicator species analysis demonstrated that 14 out of 19 common ECM fungal species showed significant preference to host plant species, suggesting that the host preference of ECM fungi was one of the most important mechanisms in structuring ECM fungal community. Accordingly, the host plant species should be taken into account in the reforestation of deforested areas due to the strong and commonly existed host preference of ECM fungi.     

Comparing activated sludge fungal community population diversity using denaturing gradient gel electrophoresis and terminal restriction fragment length polymorphism

We compared the relative values of denaturing gradient gel electrophoresis (DGGE) and terminal restriction fragment length polymorphism (T-RFLP) for profiling fungal communities in wastewater treatment plants using both ITS and 18S rRNA gene fragments as phylogenetic markers. A similar number of fungal ribotypes was obtained with both methods for the same treatment plant when the ITS primer set was used, while a greater number of ribotypes was obtained with T-RFLP compared to DGGE with the 18S rRNA primer set. Non-metric multi-dimensional scaling of presence/absence data and analysis of similarity showed that both methods could distinguish between the different plant communities at a statistically significant level (p < 0.05), regardless of which phylogenetic marker was used. The data suggest that both methods can be used preferably together to profile activated sludge fungal communities. A comparison of profiles generated with both these phylogenetic markers based on the number of ribotypes/bands, suggests that the 18S rRNA region is more discriminatory than the ITS region. Detected differences in fungal community compositions between plants probably reflect differences in their influent compositions and operational parameters.     

Initial fungal colonizer affects mass loss and fungal community development in Picea abies logs 6 yr after inoculation

Picea abies logs were inoculated with Resinicium bicolor, Fomitopsis pinicola or left un-inoculated and placed in an old-growth boreal forest. Mass loss and fungal community data were collected after 6 yr to test whether simplification of the fungal community via inoculation affects mass loss and fungal community development. Three techniques were used to survey communities: (1) observation of fruiting structures; (2) culturing on media; and (3) cloning and sequencing of ITS rDNA. Fruit body surveys detected the smallest number of species (18, 3.8 per log), DNA-based methods detected the most species (72, 31.7 per log), and culturing detected an intermediate number (23, 7.2 per log). Initial colonizer affected community development and inoculation with F. pinicola led to significantly greater mass loss. Relationships among fungal community composition, community richness and mass loss are complex and further work is needed to determine whether simplification of fungal communities affects carbon sequestration in forests.     

Soil fungal community changes in response to long-term fire cessation and N fertilization in tallgrass prairie

In grasslands, fire management and fertilization are established drivers of plant community change, but associated soil fungal responses are less well defined. We predicted that soil fungal communities would change seasonally, that decades of fire cessation and nitrogen (N) fertilization would alter fungal distributions, and that plant and fungal community change would be correlated. Surface soils were sampled monthly for 1 y from a 30-y fire by fertilization experiment to evaluate fungal community dynamics and assess correlation with plant community heterogeneity. ITS gene community composition was seasonally stable, excepting increased arbuscular mycorrhizal fungal summer abundance in the burned, fertilized treatment. Long-term treatments affected soil fungal and plant communities, with correlated heterogeneity patterns. Despite woody encroachment in the fire cessation treatment, soil fungal communities did not resemble those of forests. This study provides evidence supporting the strength of feedbacks between fungal and plant community change in response to long-term grassland fire and N management changes.     

Forest Age and Plant Species Composition Determine the Soil Fungal Community Composition in a Chinese Subtropical Forest

Fungal diversity and community composition are mainly related to soil and vegetation factors. However, the relative contribution of the different drivers remains largely unexplored, especially in subtropical forest ecosystems. We studied the fungal diversity and community composition of soils sampled from 12 comparative study plots representing three forest age classes (Young: 10–40 yrs; Medium: 40–80 yrs; Old: ≥80 yrs) in Gutianshan National Nature Reserve in South-eastern China. Soil fungal communities were assessed employing ITS rDNA pyrotag sequencing. Members of Basidiomycota and Ascomycota dominated the fungal community, with 22 putative ectomycorrhizal fungal families, where Russulaceae and Thelephoraceae were the most abundant taxa. Analysis of similarity showed that the fungal community composition significantly differed among the three forest age classes. Forest age class, elevation of the study plots, and soil organic carbon (SOC) were the most important factors shaping the fungal community composition. We found a significant correlation between plant and fungal communities at different taxonomic and functional group levels, including a strong relationship between ectomycorrhizal fungal and non-ectomycorrhizal plant communities. Our results suggest that in subtropical forests, plant species community composition is the main driver of the soil fungal diversity and community composition.     

Towards a unified paradigm for sequence‐based identification of fungi

The nuclear ribosomal internal transcribed spacer (ITS) region is the formal fungal barcode and in most cases the marker of choice for the exploration of fungal diversity in environmental samples. Two problems are particularly acute in the pursuit of satisfactory taxonomic assignment of newly generated ITS sequences: (i) the lack of an inclusive, reliable public reference data set and (ii) the lack of means to refer to fungal species, for which no Latin name is available in a standardized stable way. Here, we report on progress in these regards through further development of the UNITE database ( http://unite.ut.ee ) for molecular identification of fungi. All fungal species represented by at least two ITS sequences in the international nucleotide sequence databases are now given a unique, stable name of the accession number type (e.g. Hymenoscyphus pseudoalbidus|GU586904|SH133781.05FU), and their taxonomic and ecological annotations were corrected as far as possible through a distributed, third‐party annotation effort. We introduce the term ‘species hypothesis’ (SH) for the taxa discovered in clustering on different similarity thresholds (97–99%). An automatically or manually designated sequence is chosen to represent each such SH. These reference sequences are released ( http://unite.ut.ee/repository.php ) for use by the scientific community in, for example, local sequence similarity searches and in the QIIME pipeline. The system and the data will be updated automatically as the number of public fungal ITS sequences grows. We invite everybody in the position to improve the annotation or metadata associated with their particular fungal lineages of expertise to do so through the new Web‐based sequence management system in UNITE.     

Changes in ectomycorrhizal community structure on two containerized oak hosts across an experimental hydrologic gradient

Shifts in ectomycorrhizal (ECM) community structure were examined across an experimental hydrologic gradient on containerized seedlings of two oak species, Quercus montana and Quercus palustris, inoculated from a homogenate of roots from mature oak trees. At the end of one growing season, seedlings were harvested, roots were sorted by morphotype, and proportional colonization of each type was determined. DNA was subsequently extracted from individual root tips for polymerase chain reaction, restriction fragment length polymorphism, and rDNA sequencing of the ITS1/5.8S/ITS2 region to determine identities of fungal morphotypes. Twelve distinct molecular types were identified. Analysis of similarity showed that ECM fungal assemblages shifted significantly in composition across the soil moisture gradient. Taxa within the genus Tuber and the family Thelephoraceae were largely responsible for the changes in fungal assemblages. There were also significant differences in ECM community assemblages between the two oak host species. These results demonstrate that the structure of ECM fungal communities depends on both the abiotic and biotic environments and can shift with changes in soil moisture as well as host plant, even within the same genus.     

Fungal community diversity and soil health in intensive potato cropping systems of the east Po valley, northern Italy

An ecological approach was used to investigate the relationship between diversity of soil fungal communities and soil‐borne pathogen inoculum in a potato growing area of northern Italy affected by yield decline. The study was performed in 14 sites with the same tillage management practices: 10 named ‘potato sites’, that for many years had been intensely cultivated with potatoes, and 4 named ‘rotation sites’, subject to a 4‐year rotation without potatoes or any recurrent crop for many years. Fungal communities were recorded using conventional (soil fungi by plate count and endophytic fungi as infection frequency on pot‐grown potato plant roots in soil samples) and molecular approaches [Basidiomycetes and Ascomycetes with specific and denaturing gradient gel electrophoresis (DGGE) analysis]. Diversity of fungal communities in potato sites was significantly lower than that in rotation sites. In addition, fungal communities in rotation sites showed lower Berger–Parker dominance than those in the potato sites, suggesting that rotation sites had a higher diversity as well as a better fungal community balance than potato sites. The ANalysis Of SIMilarity test of soil fungi and root endophytic fungi revealed that the two cropping systems differed significantly for species composition. Root endophytic fungal communities showed a greater ability to colonise potato roots in soil samples from potato sites than those from rotation sites. Moreover, the majority of endophytic root fungal community species in potato sites belonged to the potato root rot complex and storage disease (Colletotrichum coccodes, Fusarium solani and Fusarium oxysporum), while those in rotation sites were mainly ubiquitous or saprobic fungi. Soil rDNA analyses showed that Ascomycetes were much more frequent than Basidiomycetes in all the soils examined. DGGE analysis, with the Ascomycete‐specific primer (ITS1F/ITS4A), did not reveal distinctions between the communities found at the potato and rotation sites, although the same analysis showed differences between the communities of Basidiomycetes (specific primer ITS1F/ITS4B). These findings showed that recurrent potato cropping affected diversity and composition of soil fungal communities and induced a shift in specialisation of the endophytic fungi towards potato.     

Influence of Elevated CO2 on the Fungal Community in a Coastal Scrub Oak Forest Soil Investigated with Terminal-Restriction Fragment Length Polymorphism Analysis

Sixteen open-top chambers (diameter, 3.66 m) were established in a scrub oak habitat in central Florida where vegetation was removed in a planned burn prior to chamber installation. Eight control chambers have been continuously exposed to ambient air and eight have been continuously exposed to elevated CO₂ at twice-ambient concentration (~700 ppm) for 5 years. Soil cores were collected from each chamber to examine the influence of elevated atmospheric CO₂ on the fungal community in different soil fractions. Each soil sample was physically fractionated into bulk soil, rhizosphere soil, and roots for separate analyses. Changes in relative fungal biomass were estimated by the ergosterol technique. In the bulk soil and root fractions, a significantly increased level of ergosterol was detected in the elevated CO₂ treatments relative to ambient controls. Fungal community composition was determined by terminal-restriction fragment length polymorphism (T-RFLP) analysis of the internal transcribed spacer (ITS) region. The specificities of different ITS primer sets were evaluated against plant and fungal species isolated from the experimental site. Changes in community composition were assessed by principal component analyses of T-RFLP profiles resolved by capillary electrophoresis. Fungal species richness, defined by the total number of terminal restriction fragments, was not significantly affected by either CO₂ treatment or soil fraction.