The ITS region as a target for characterization of fungal communities using emerging sequencing technologies

Molecular Identification and Antifungal Susceptibility of Clinical Isolates of Sporothrix schenckii Complex in Medellin, Colombia

BACKGROUND

Sporotrichosis is a subcutaneous mycosis that affects humans and other animals. Infection prevails in tropical and subtropical countries. Until a few years ago, it was considered that two varieties of Sporothrix schenckii caused this mycosis, but by applying molecular taxonomic markers, it has been demonstrated that there are several cryptic species within S. schenckii complex which varies in susceptibility, virulence, and geographic distribution.

OBJECTIVE

This study aimed to identify the clinical isolates of Sporothrix spp. from patients with sporotrichosis in Medellin, Colombia, using two markers and to evaluate the in vitro susceptibility to itraconazole.

METHODS

Thirty-four clinical isolates of Sporothrix spp. from Colombia, three from Mexico, and one from Guatemala were identified through sequencing of the noncoding region ITS-1 + 5.8SDNAr + ITS-2 and of the fragment containing exons 3 and 4 of the β-tubulin gene. Clinical isolate sequences were compared with GenBank reference sequences using the BLASTN tool, and then, phylogenetic analysis was performed. Besides, the in vitro susceptibility to itraconazole was evaluated by determining the minimum inhibitory concentrations according to the CLSI M38-A2 method.

RESULTS

Clinical isolates were identified by morphology as Sporothrix spp. Using the molecular markers, ITS and β-tubulin, isolates were identified as S. schenckii sensu stricto (25) and Sporothrix globosa (13). Susceptibility to itraconazole was variable among clinical isolates.

CONCLUSION

This is the first scientific publication that identifies species that cause sporotrichosis in Colombia, along with the antifungal susceptibility to itraconazole.

Small-Scale Soil Microbial Community Heterogeneity Linked to Landform Historical Events on King George Island, Maritime Antarctica

Although research on microbial biogeography has made great progress in the past decade, distributions of terrestrial microbial communities in extreme environments such as Antarctica are not well understood. In addition, knowledge of whether and how historical contingencies affect microbial distributions at small spatial scales is lacking. Here, we analyzed soil-borne microbial (bacterial, archaeal, and fungal) communities in 12 quadrat plots around the Fildes Region of King George Island, maritime Antarctica, and the communities were divided into two groups according to the soil elemental compositions and environmental attributes of Holocene raised beach and Tertiary volcanic stratigraphy. Prokaryotic communities of the two groups were well separated; the prokaryotic data were primarily correlated with soil elemental compositions and were secondly correlated with environmental attributes (e.g., soil pH, total organic carbon, NO3 -, and vegetation coverage; Pearson test, r = 0.59 vs. 0.52, both P < 0.01). The relatively high abundance of P, S, Cl, and Br in Group 1 (Holocene raised beach site) was likely due to landform uplift. Lithophile-elements (Si, Al, Ca, Sr, Ti, V, and Fe) correlated with prokaryotic communities in Group 2 may have originated from weathering of Tertiary volcanic rock. No significant correlations were found between the fungal community distribution and both the soil elemental composition and environmental attributes in this study; however, Monte Carlo tests revealed that elements Sr and Ti, soil pH, sampling altitude, and moss and lichen species numbers had significant impacts on fungal communities. The elements and nutrients accumulated during the formation of different landforms influenced the development of soils, plant growth, and microbial communities, and this resulted in small-scale spatially heterogeneous biological distributions. These findings provide new evidence that geological evolutionary processes in the Fildes Region were crucial to its microbial community development, and the results highlight that microbial distribution patterns are the legacies of historical events at this small spatial scale. Based on this study, the ice-free regions in maritime Antarctica represent suitable research sites for studying the influence of geomorphological features on microbial distributions, and we envision the possibility of a site-specific landform assignment through the analysis of the soil prokaryotic community structure.

Effects of phenolic acids from ginseng rhizosphere on soil fungi structure, richness and diversity in consecutive monoculturing of ginseng

Ginseng yield and quality are seriously compromised by consecutive monoculturing in northeastern China. The imbalance of soil fungi communities and autotoxicity of ginseng are the major factors in consecutive monoculturing ginseng crops. Soil fungal communities were identified using Illumina MiSeq sequencing, applied to soils that consecutively cultured ginseng (CCG) for six years and new forest soil (NFS), or receiving application of phenolic acids (PAs). The CCG field received five treatments with five different phenolic acids, including gallic acid (GA), salicylic acid (SA), 3-phenylpropionic acid (3-PA), benzoic acid (BA) and cinnamic acid (CA), which were detected from ginseng rhizosphere in consecutive cropping soil. Fungal richness, fungi diversity, community composition, relative taxon abundances, root rot disease, and growth rate were compared among the different treatments. 579 fungal operational taxonomic units at 97% ITS sequence identity were found among 201,617 sequence reads derived from 18 separate soil samples. Members of the phylum Ascomycota dominated the soil fungal communities, and putative pathogens, such as Fusarium, Gibberella and Nectriaceae_unclassified which may include the abundant sexual morph of Cylindrocarpon destructans, showed higher relative abundances in the CCG fields. Compared to the CCG and NFS fields, PAs (except CA) enhanced the fungi richness and decreased fungi diversity. Cluster analysis indicated that the PAs (except CA) changed the fungi structure in a uniform way. PAs stimulate root rot disease and enhance disease severity, restricting plant growth. The results suggest that the PAs (except CA) may enhance the fungi richness, decrease the fungi diversity and changed the fungi structure to increase fungal pathogen loads, which could explain the declined yield and quality of ginseng in consecutively monocultured ginseng crops.

Database establishment for the secondary fungal DNA barcode translational elongation factor 1α (TEF1α) 1

With new or emerging fungal infections, human and animal fungal pathogens are a growing threat worldwide. Current diagnostic tools are slow, non-specific at the species and subspecies levels, and require specific morphological expertise to accurately identify pathogens from pure cultures. DNA barcodes are easily amplified, universal, short species-specific DNA sequences, which enable rapid identification by comparison with a well-curated reference sequence collection. The primary fungal DNA barcode, ITS region, was introduced in 2012 and is now routinely used in diagnostic laboratories. However, the ITS region only accurately identifies around 75% of all medically relevant fungal species, which has prompted the development of a secondary barcode to increase the resolution power and suitability of DNA barcoding for fungal disease diagnostics. The translational elongation factor 1α (TEF1α) was selected in 2015 as a secondary fungal DNA barcode, but it has not been implemented into practice, due to the absence of a reference database. Here, we have established a quality-controlled reference database for the secondary barcode that together with the ISHAM-ITS database, forms the ISHAM barcode database, available online at http://its.mycologylab.org/ . We encourage the mycology community for active contributions.

Evaluation of the ribosomal DNA internal transcribed spacer (ITS), specifically ITS1 and ITS2, for the analysis of fungal diversity by deep sequencing

The nuclear ribosomal DNA internal transcribed spacer (ITS) has been widely used to assess the fungal composition in different environments by deep sequencing. To evaluate the ITS in the analysis of fungal diversity, comparisons of the clustering and taxonomy generated by sequencing with different portions of the whole fragment were conducted in this study. For a total of 83,120 full-length ITS sequences obtained from the UNITE database, it was found that, on average, ITS1 varied more than ITS2 within the kingdom Fungi; this variation included length and GC content variations and polymorphisms, with some polymorphisms specific to particular fungal groups. The taxonomic accuracy for ITS was higher than that for ITS1 or ITS2. The commonly used operational taxonomic unit (OTU) for evaluating fungal diversity and richness assigned several species to a single OTU even with clustering at 99.00% sequence similarity. The clustering and taxonomic capacities did not differ between ITS1 and ITS2. However, the OTU commonality between ITS1 and ITS2 was very low. To test this observation further, 219,741 pyrosequencing reads, including 39,840 full-length ITS sequences, were obtained from 10 soil samples and were clustered into OTUs. The pyrosequencing results agreed with the results of the in silico analysis. ITS1 might overestimate the fungal diversity and richness. Analyses using ITS, ITS1 and ITS2 yielded several different taxa, and the taxonomic preferences for ITS and ITS2 were similar. The results demonstrated that ITS2 alone might be a more suitable marker for revealing the operational taxonomic richness and taxonomy specifics of fungal communities when the full-length ITS is not available.

ITSxpress: Software to rapidly trim internally transcribed spacer sequences with quality scores for marker gene analysis

The internally transcribed spacer (ITS) region between the small subunit ribosomal RNA gene and large subunit ribosomal RNA gene is a widely used phylogenetic marker for fungi and other taxa. The eukaryotic ITS contains the conserved 5.8S rRNA and is divided into the ITS1 and ITS2 hypervariable regions. These regions are variable in length and are amplified using primers complementary to the conserved regions of their flanking genes. Previous work has shown that removing the conserved regions results in more accurate taxonomic classification. An existing software program, ITSx, is capable of trimming FASTA sequences by matching hidden Markov model profiles to the ends of the conserved genes using the software suite HMMER. ITSxpress was developed to extend this technique from marker gene studies using Operational Taxonomic Units (OTU's) to studies using exact sequence variants; a method used by the software packages Dada2, Deblur, QIIME 2, and Unoise. The sequence variant approach uses the quality scores of each read to identify sequences that are statistically likely to represent real sequences. ITSxpress enables this by processing FASTQ rather than FASTA files. The software also speeds up the trimming of reads by a factor of 14-23 times on a 4-core computer by temporarily clustering highly similar sequences that are common in amplicon data and utilizing optimized parameters for Hmmsearch. ITSxpress is available as a QIIME 2 plugin and a stand-alone application installable from the Python package index, Bioconda, and Github.

A highly efficient electrophoretic method for discrimination between two Neoscytalidium species using a specific fungal internal transcribed spacer (ITS) fragment

Neoscytalidium (or N.) dimidiatum and N. novaehollandiae are two aggressive plant pathogenic species that affect several agricultural crops. Early detection and identification of these fungi are of critical importance to bring about the effective minimization to the threat they pose to the infected plants. Herein, two species of Neoscytalidium were rapidly discriminated by utilizing the rRNA internal transcribed (ITS4-5.8S-ITS5) PCR primers. A total of 100 isolates of Neoscytalidium species, which were isolated from Iraqi canker-infected fig trees, were included in this study. Two discrete electrophoretic PCR bands were observed in Neoscytalidium isolates-A-variants were about 546 bp, while B-variants were about 993 bp in length. The comprehensive phylogenetic analysis of both DNA variants revealed that A-variants resided between N. novaehollandiae and N. hyalinum, while B-variants were closely related to N. dimidiatum. Furthermore, the highly specific re-constructed tree of both electrophoretic variants demonstrated that B-variants share a high similarity with N. novaehollandiae. Additionally, the secondary structures for both variants were predicted computationally to reveal the structural patterns that each variant follows. In conclusion, a small rRNA locus comprising 22 nucleotides that differs in the two variants is potentially responsible for this species-specific classification. The main divergence in the amplified loci led to the classification of these fungal variants into two main species, namely N. dimidiatum and N. novaehollandiae, demonstrating that the amplification by ITS4-ITS5 rRNA fragment is a beneficial strategy that can be employed for the assessment of Neoscytalidium diversity in the natural ecosystems.

Trends in bacterial and fungal communities in ant nests observed with Terminal-Restriction Fragment Length Polymorphism (T-RFLP) and Next Generation Sequencing (NGS) techniques-validity and compatibility in ecological studies

Microbes are ubiquitous and often occur in functionally and taxonomically complex communities. Unveiling these community dynamics is one of the main challenges of microbial research. Combining a robust, cost effective and widely used method such as Terminal Restriction Fragment Length Polymorphism (T-RFLP) with a Next Generation Sequencing (NGS) method (Illumina MiSeq), offers a solid alternative for comprehensive assessment of microbial communities. Here, these two methods were combined in a study of complex bacterial and fungal communities in the nest mounds of the ant Formica exsecta, with the aim to assess the degree to which these methods can be used to complement each other. The results show that these methodologies capture similar spatiotemporal variations, as well as corresponding functional and taxonomical detail, of the microbial communities in a challenging medium consisting of soil, decomposing plant litter and an insect inhabitant. Both methods are suitable for the analysis of complex environmental microbial communities, but when combined, they complement each other well and can provide even more robust results. T-RFLP can be trusted to show similar general community patterns as Illumina MiSeq and remains a good option if resources for NGS methods are lacking.

Diversity of Fungal Communities in Heshang Cave of Central China Revealed by Mycobiome-Sequencing

Deciphering of the mycobiome in pristine karst caves has been impeded by constraints of remote locations, inaccessibility to specimens and technical limitations, which greatly restricted in-depth understanding of mycobiomes in subterranean ecosystem. Here, mycobiomes of Heshang Cave in south-western karst region of China were investigated by Illumina HiSeq sequencing of fungal rRNA-ITS1 gene across different habitats. In total 793,502 ITS1 reads and 2,179 OTUs from 778 Mb reads after stringent quality control (Q30) and 453 genera, 72 orders and 19 classes within 6 phyla were detected. Ascomycota (42% OTUs) dominated across the five habitats. Shannon-Wiener index varied from 1.25 to 7.62 and community richness was highest in drip waters, followed by weathered rocks, bat guanos, sediments, and air samples. Mycobiomes displayed specificity to five habitats and more distinct OTUs were found in weathered rocks (12%) and drip waters (9%). In contrast, only 6.60% core OTUs were shared by five habitats. Notably, weathered rocks possessed more indicator groups and were revealed for the first time to be dominated by Sordariomycetes (43%). The community richness of air mycobiomes increased from cave entrance to the innermost part and dominated by the indicator groups of Penicillium mallochii (>30%) and P. herquei (>9%). Our work represents the largest attempt to date to a systematical investigation of oligotrophic solution-cave-associated mycobiomes in China. Our discovery of high diversity of mycobiomes in Heshang Cave also suggests that eukaryotic microorganisms may play a crucial role in subsurface environments.

Molecular Identification of Endophytic Fungi from Banana Leaves (Musa spp.)

Endophytic fungi are part of microbial community found in various types of plant tissues including the leave, and display a range of symbiotic interactions with the plant host. In this study, endophytic fungi isolated from banana leaves were identified using ITS (Internal Transcribed Spacer region) sequences of which 10 genera comprising 17 species were molecularly identified. Endophytic fungal species identified were Nigrospora oryzae, Nigrospora sphaerica, Colletotrichum gloeosporioides, Colletotrichum siamense, Fusarium equiseti, Fusarium chlamydosporum, Phoma sorghina, Pestalotiopsis oxyanthi, Pestalotiopsis theae, Pestalotiopsis eugeniae, Penicillium steckii, Penicillium purpurogenum, Bipolaris papendorfii, Bipolaris sp., Lasidiodiplodia theobromae, Cochliobolus intermedius dan Aspergillus niger. The present study showed that several endophytic fungal genera/species are common plant pathogen and there is a possibility that these endophytes can become pathogenic. Some of the fungal endophyte might be mutualist or saprophyte.

ITS2 metabarcoding analysis complements lichen mycobiome diversity data

Lichen thalli harbor complex fungal communities (mycobiomes) of species with divergent trophic and ecological strategies. The complexity and diversity of lichen mycobiomes are still largely unknown, despite surveys combining culture-based methods and high-throughput sequencing (HTS). The results of such surveys are strongly influenced by the barcode locus chosen, its sensitivity in discriminating taxa, and the depth to which public sequence repositories cover the phylogenetic spectrum of fungi. Here, we use HTS of the internal transcribed spacer 2 (ITS2) to assess the taxonomic composition and diversity of a well-characterized, alpine rock lichen community that includes thalli symptomatically infected by lichenicolous fungi as well as asymptomatic thalli. Taxa belonging to the order Chaetothyriales are the major components of the observed lichen mycobiomes. We predict sequences representative of lichenicolous fungi characterized morphologically and assess their asymptomatic presence in lichen thalli. We demonstrated the limitations of metabarcoding in fungi and show how the estimation of species diversity widely differs when ITS1 or ITS2 are used as barcode, and particularly biases the detection of Basidiomycota. The complementary analysis of both ITS1 and ITS2 loci is therefore required to reliably estimate the diversity of lichen mycobiomes.

Molecular survey of basidiomycetes and divergence time estimation: An Indian perspective

This study outlines the biodiversity of mushrooms of India. It reveals the molecular biodiversity and divergence time estimation of basidiomycetes from Gujarat, India. A total of 267 mushrooms were collected from 10 locations across the state. 225 ITS sequences were generated belonging to 105 species, 59 genera and 29 families. Phylogenetic analysis of Agaricaceae reveals monophyletic clade of Podaxis differentiating it from Coprinus. Further, the ancient nature of Podaxis supports the hypothesis that gasteroid forms evolved from secotioid forms. Members of Polyporaceae appeared polyphyletic. Further, our results of a close phylogenetic relationship between Trametes and Lenziteslead us to propose that the genera Trametes may by enlarged to include Lenzites. The tricholomatoid clade shows a clear demarcation for Entolomataceae. However, Lyophyllaceae and Tricholomataceae could not be distinguished clearly. Distribution studies of the mushrooms showed omnipresence of Ganoderma and Schizophyllum. Further, divergence time estimation shows that Dacrymycetes evolved in the Neoproterozoic Era and Hymenochaetales diverged from Agaricomycetes during the Silurian period.

The species of Coleosporium (Pucciniales) on Solidago in North America

Species of Coleosporium (Pucciniales) are rust fungi that typically alternate between pines and angiosperms. In North America, species of Coleosporium often infect Solidago (goldenrods), although their taxonomy on these hosts is unresolved. Joseph. C. Arthur and George B. Cummins regarded these as a single species, Coleosporium solidaginis (fide Arthur) or C. asterum (fide Cummins), but later inoculation studies demonstrated the presence of more than one species, distinguishable by their aecial hosts. A more recent taxonomic study of Coleosporium found that specimens on Solidago identified as C. asterum in North America were not conspecific with the type, which is from Japan, prompting the present study. Herein, we conducted a systematic study on ca. 60 collections of Coleosporium infecting species of Asteraceae from North America using regions of ribosomal DNA and morphology of teliospores and basidia. Our data indicate at least three species of Coleosporium occur on Solidago in North America, C. solidaginis, C. montanum comb. nov., which is proposed for the taxon that has commonly been identified as C. asterum in North America, and C. delicatulum, all of which can be differentiated by morphology of their basidia. In addition, the challenges of marker selection for molecular barcoding of rust fungi is discussed.

DNA metabarcoding uncovers fungal diversity of mixed airborne samples in Italy

Fungal spores and mycelium fragments are particles which become and remain airborne and have been subjects of aerobiological studies. The presence and the abundance of taxa in aerobiological samples can be very variable and impaired by changeable climatic conditions. Because many fungi produce mycotoxins and both their mycelium fragments and spores are potential allergens, monitoring the presence of these taxa is of key importance. So far data on exposure and sensitization to fungal allergens are mainly based on the assessment of few, easily identifiable taxa and focused only on certain environments. The microscopic method used to analyze aerobiological samples and the inconspicuous fungal characters do not allow a in depth taxonomical identification. Here, we present a first assessment of fungal diversity from airborne samples using a DNA metabarcoding analysis. The nuclear ITS2 region was selected as barcode to catch fungal diversity in mixed airborne samples gathered during two weeks in four sites of North-Eastern and Central Italy. We assessed the taxonomic composition and diversity within and among the sampled sites and compared the molecular data with those obtained by traditional microscopy. The molecular analyses provide a tenfold more comprehensive determination of the taxa than the traditional morphological inspections. Our results prove that the metabarcoding analysis is a promising approach to increases quality and sensitivity of the aerobiological monitoring. The laboratory and bioinformatic workflow implemented here is now suitable for routine, high-throughput, regional analyses of airborne fungi.

The Intestinal Mycobiota in Wild Zebrafish Comprises Mainly Dothideomycetes While Saccharomycetes Predominate in Their Laboratory-Reared Counterparts

As an integral part of the resident microbial community of fish intestinal tract, the mycobiota is expected to play important roles in health and disease resistance of the host. The composition of the diverse fungal communities, which colonize the intestine, is greatly influenced by the host, their diet and geographic origin. Studies of fungal communities are rare and the majority of previous studies have relied on culture-based methods. In particular, fungal communities in fish are also poorly characterized. The aim of this study was to provide an in-depth overview of the intestinal mycobiota in a model fish species (zebrafish, Danio rerio) and to determine differences in fungal composition between wild and captive specimens. We have profiled the intestinal mycobiota of wild-caught (Sharavati River, India), laboratory-reared (Bodø, Norway) and wild-caught-laboratory-kept (Uttara, India) zebrafish by sequencing the fungal internal transcribed spacer 2 region on the Illumina MiSeq platform. Wild fish were exposed to variable environmental factors, whereas both laboratory groups were kept in controlled conditions. There were also differences in husbandry practices at Bodø and Uttara, particularly diet. Zebrafish from Bodø were reared in the laboratory for over 10 generations, while wild-caught-laboratory-kept fish from Uttara were housed in the laboratory for only 2 months before sample collection. The intestine of zebrafish contained members of more than 15 fungal classes belonging to the phyla Ascomycota, Basidiomycota, and Zygomycota. Fungal species richness and diversity distinguished the wild-caught and laboratory-reared zebrafish communities. Wild-caught zebrafish-associated mycobiota comprised mainly Dothideomycetes in contrast to their Saccharomycetes-dominated laboratory-reared counterparts. The predominant Saccharomycetes in laboratory-reared fish belonged to the saprotrophic guild. Another characteristic feature of laboratory-reared fish was the significantly higher abundance of Cryptococcus (Tremellomycetes) compared to wild fish. This pioneer study has shed light into the differences in the intestinal fungal communities of wild-caught and laboratory-reared zebrafish and the baseline data generated will enrich our knowledge on fish mycobiota.

Metabarcoding analysis of eukaryotic microbiota in the gut of HIV-infected patients

Research on the relationship between changes in the gut microbiota and human disease, including AIDS, is a growing field. However, studies on the eukaryotic component of the intestinal microbiota have just begun and have not yet been conducted in HIV-infected patients. Moreover, eukaryotic community profiling is influenced by the use of different methodologies at each step of culture-independent techniques. Herein, initially, four DNA extraction protocols were compared to test the efficiency of each method in recovering eukaryotic DNA from fecal samples. Our results revealed that recovering eukaryotic components from fecal samples differs significantly among DNA extraction methods. Subsequently, the composition of the intestinal eukaryotic microbiota was evaluated in HIV-infected patients and healthy volunteers through clone sequencing, high-throughput sequencing of nuclear ribosomal internal transcribed spacers 1 (ITS1) and 2 (ITS2) amplicons and real-time PCRs. Our results revealed that not only richness (Chao-1 index) and alpha diversity (Shannon diversity) differ between HIV-infected patients and healthy volunteers, depending on the molecular strategy used, but also the global eukaryotic community composition, with little overlapping taxa found between techniques. Moreover, our results based on cloning libraries and ITS1/ITS2 metabarcoding sequencing showed significant differences in fungal composition between HIV-infected patients and healthy volunteers, but without distinct clusters separating the two groups. Malassezia restricta was significantly more prevalent in fecal samples of HIV-infected patients, according to cloning libraries, whereas operational taxonomic units (OTUs) belonging to Candida albicans and Candida tropicalis were significantly more abundant in fecal samples of HIV-infected patients compared to healthy subjects in both ITS subregions. Finally, real-time PCR showed the presence of Microsporidia, Giardia lamblia, Blastocystis and Hymenolepis diminuta in different proportions in fecal samples from HIV patients as compared to healthy individuals. Our work revealed that the use of different sequencing approaches can impact the perceived eukaryotic diversity results of the human gut. We also provide a more comprehensive view of the eukaryotic community in the gut of HIV-infected patients through the complementarity of the different molecular techniques used. Combining these various methodologies may provide a gold standard for a more complete characterization of the eukaryotic microbiome in future studies.

Analysis of the influence of living environment and age on vaginal fungal microbiome in giant pandas (Ailuropoda melanoleuca) by high throughput sequencing

A recent study has described the normal vaginal bacterial community in giant pandas, but there is a lack of knowledge of the fungal community residing in the vagina of giant pandas. In order to comprehensively understand the vaginal fungal microbial diversity and abundance in giant pandas, high throughput sequencing was used to analyse the ITS1 region, based on thirteen samples taken from the pandas' vaginas, which were grouped by sampling points and age. The results showed that the most abundant phyla were Basidiomycota (73.37%), followed by Ascomycota (20.04%), Zygomycota (5.23%), Glomeromycota (0.014%) and Chytridiomycota (0.006%). At the genus level, Guehomyces (37.92%) was the most abundant, followed by Cladosporium (9.072%), Trichosporon (6.2%) and Mucor (4.97%). Furthermore, Candida only accounted for a low percentage of the vaginal fungal community. With the saturation of rarefaction curves and fungal diversity indices, the samples from Dujiangyan and Chungking Safari Park (DC group) showed a higher fungal species richness and diversity than other living environments. Shannon diversity indices showed significant difference between group WL (Wolong nature reserve) and DC (P < .05). Additionally, a higher diversity was found in ten to fifteen years old (Group 2) than other groups. Group 2 and Group 3 displayed significant differences in the diversities of their vaginal fungal communities (P < .05). These data that has been collected from this research will be helpful for further study to improve the reproductive status of giant pandas.

Estimating and mitigating amplification bias in qualitative and quantitative arthropod metabarcoding

Amplicon based metabarcoding promises rapid and cost-efficient analyses of species composition. However, it is disputed whether abundance estimates can be derived from metabarcoding due to taxon specific PCR amplification biases. PCR-free approaches have been suggested to mitigate this problem, but come with considerable increases in workload and cost. Here, we analyze multilocus datasets of diverse arthropod communities, to evaluate whether amplification bias can be countered by (1) targeting loci with highly degenerate primers or conserved priming sites, (2) increasing PCR template concentration, (3) reducing PCR cycle number or (4) avoiding locus specific amplification by directly sequencing genomic DNA. Amplification bias is reduced considerably by degenerate primers or targeting amplicons with conserved priming sites. Surprisingly, a reduction of PCR cycles did not have a strong effect on amplification bias. The association of taxon abundance and read count was actually less predictable with fewer cycles. Even a complete exclusion of locus specific amplification did not exclude bias. Copy number variation of the target loci may be another explanation for read abundance differences between taxa, which would affect amplicon based and PCR free methods alike. As read abundance biases are taxon specific and predictable, the application of correction factors allows abundance estimates.

Culturomics and Amplicon-based Metagenomic Approaches for the Study of Fungal Population in Human Gut Microbiota

Herein, the mycobiota was characterized in fecal samples from sick patients and healthy subjects, collected from different geographical locations and using both culturomics and amplicon-based metagenomics approaches. Using the culturomics approach, a total of 17,800 fungal colonies were isolated from 14 fecal samples, and resulted in the isolation of 41 fungal species, of which 10 species had not been previously reported in the human gut. Deep sequencing of fungal-directed ITS1 and ITS2 amplicons led to the detection of a total of 142 OTUs and 173 OTUs from the ITS1 and ITS2 regions, respectively. Ascomycota composed the largest fraction of the total OTUs analyzed (78.9% and 68.2% of the OTUs from the ITS1 and ITS2 regions, respectively), followed by Basidiomycota (16.9% and 30.1% of the OTUs from the ITS1 and ITS2 regions, respectively). Interestingly, the results demonstrate that the ITS1/ITS2 amplicon sequencing provides different information about gut fungal communities compared to culturomics, though both approaches complete each other in assessing fungal diversity in fecal samples. We also report higher fungal diversity and abundance in patients compared to healthy subjects. In conclusion, combining both culturomic and amplicon-based metagenomic approaches may be a novel strategy towards analyzing fungal compositions in the human gut.

Forgotten fungi-the gut mycobiome in human health and disease

The human body is home to a complex and diverse microbial ecosystem that plays a central role in host health. This includes a diversity of fungal species that is collectively referred to as our 'mycobiome'. Although research into the mycobiome is still in its infancy, its potential role in human disease is increasingly recognised. Here we review the existing literature available on the human mycobiota with an emphasis on the gut mycobiome, including how fungi interact with the human host and other microbes. In doing so, we provide a comprehensive critique of the methodologies available to research the human mycobiota as well as highlighting the latest research findings from mycological surveys of different groups of interest including infants, obese and inflammatory bowel disease cohorts. This in turn provides new insights and directions for future studies in this burgeoning research area.

Foliar fungi of Betula pendula: impact of tree species mixtures and assessment methods

Foliar fungi of silver birch (Betula pendula) in an experimental Finnish forest were investigated across a gradient of tree species richness using molecular high-throughput sequencing and visual macroscopic assessment. We hypothesized that the molecular approach detects more fungal taxa than visual assessment, and that there is a relationship among the most common fungal taxa detected by both techniques. Furthermore, we hypothesized that the fungal community composition, diversity, and distribution patterns are affected by changes in tree diversity. Sequencing revealed greater diversity of fungi on birch leaves than the visual assessment method. One species showed a linear relationship between the methods. Species-specific variation in fungal community composition could be partially explained by tree diversity, though overall fungal diversity was not affected by tree diversity. Analysis of specific fungal taxa indicated tree diversity effects at the local neighbourhood scale, where the proportion of birch among neighbouring trees varied, but not at the plot scale. In conclusion, both methods may be used to determine tree diversity effects on the foliar fungal community. However, high-throughput sequencing provided higher resolution of the fungal community, while the visual macroscopic assessment detected functionally active fungal species.

The importance of molecular markers and primer design when characterizing biodiversity from environmental DNA

Environmental DNA (eDNA) comprises DNA fragments that have been shed into the environment by organisms, and which can be extracted from environmental samples such as water or soil. Characterization of eDNA can allow researchers to infer the presence or absence of species from a particular site without the need to locate and identify individuals, and therefore may provide an extremely valuable tool for quantifying biodiversity. However, as is often the case with relatively new protocols, methodological challenges remain. A number of earlier reviews have discussed these challenges, but none have provided extensive treatment of the critical decisions surrounding molecular markers and primer development for use in eDNA assays. This review discusses a number of options and approaches that can be used when determining which primers and gene regions are most appropriate for either targeted species detection or metabarcoding macro-organisms from eDNA. The latter represents a new field that is growing rapidly, and which has the potential to revolutionize future assessments of community and ecosystem diversity.

Vertical Structure of Phyllosphere Fungal Communities in a Tropical Forest in Thailand Uncovered by High-Throughput Sequencing

Phyllosphere fungi harbor a tremendous species diversity and play important ecological roles. However, little is known about their distribution patterns within forest ecosystems. We examined how species diversity and community composition of phyllosphere fungi change along a vertical structure in a tropical forest in Thailand. Fungal communities in 144 leaf samples from 19 vertical layers (1.28-34.4 m above ground) of 73 plant individuals (27 species) were investigated by metabarcoding analysis using Ion Torrent sequencing. In total, 1,524 fungal operational taxonomic units (OTUs) were detected among 890,710 reads obtained from the 144 leaf samples. Taxonomically diverse fungi belonging to as many as 24 orders of Ascomycota and 21 orders of Basidiomycota were detected, most of which inhabited limited parts of the lowest layers closest to the forest floor. Species diversity of phyllosphere fungi was the highest in the lowest layers closest to the forest floor, decreased with increasing height, and lowest in the canopy; 742 and 55 fungal OTUs were detected at the lowest and highest layer, respectively. On the layers close to the forest floor, phyllosphere fungal communities were mainly composed of low frequency OTUs and largely differentiated among plant individuals. Conversely, in the canopy, fungal communities consisted of similar OTUs across plant individuals, and as many as 86.1%-92.7% of the OTUs found in the canopy (≥22 m above ground) were also distributed in the lower layers. Overall, our study showed the variability of phyllosphere fungal communities along the vertical gradient of plant vegetation and environmental conditions, suggesting the significance of biotic and abiotic variation for the species diversity of phyllosphere fungi.

Optimization of Xylanase Production through Response Surface Methodology by Fusarium sp. BVKT R2 Isolated from Forest Soil and Its Application in Saccharification

Xylanses are hydrolytic enzymes with wide applications in several industries like biofuels, paper and pulp, deinking, food, and feed. The present study was aimed at hitting at high yield xylanase producing fungi from natural resources. Two highest xylanase producing fungal isolates-Q12 and L1 were picked from collection of 450 fungal cultures for the utilization of xylan. These fungal isolates-Q12 and L1 were identified basing on ITS gene sequencing analysis as Fusarium sp. BVKT R2 (KT119615) and Fusarium strain BRR R6 (KT119619), respectively with construction of phylogenetic trees. Fusarium sp. BVKT R2 was further optimized for maximum xylanase production and the interaction effects between variables on production of xylanase were studied through response surface methodology. The optimal conditions for maximal production of xylanase were sorbitol 1.5%, yeast extract 1.5%, pH of 5.0, Temperature of 32.5°C, and agitation of 175 rpm. Under optimal conditions, the yields of xylanase production by Fusarium sp. BVKT R2 was as high as 4560 U/ml in SmF. Incubation of different lignocellulosic biomasses with crude enzyme of Fusarium sp. BVKT R2 at 37°C for 72 h could achieve about 45% saccharification. The results suggest that Fusarium sp. BVKT R2 has potential applications in saccharification process of biomass.

Mycofier: a new machine learning-based classifier for fungal ITS sequences

BACKGROUND

The taxonomic and phylogenetic classification based on sequence analysis of the ITS1 genomic region has become a crucial component of fungal ecology and diversity studies. Nowadays, there is no accurate alignment-free classification tool for fungal ITS1 sequences for large environmental surveys. This study describes the development of a machine learning-based classifier for the taxonomical assignment of fungal ITS1 sequences at the genus level.

RESULTS

A fungal ITS1 sequence database was built using curated data. Training and test sets were generated from it. A Naïve Bayesian classifier was built using features from the primary sequence with an accuracy of 87 % in the classification at the genus level.

CONCLUSIONS

The final model was based on a Naïve Bayes algorithm using ITS1 sequences from 510 fungal genera. This classifier, denoted as Mycofier, provides similar classification accuracy compared to BLASTN, but the database used for the classification contains curated data and the tool, independent of alignment, is more efficient and contributes to the field, given the lack of an accurate classification tool for large data from fungal ITS1 sequences. The software and source code for Mycofier are freely available at https://github.com/ldelgado-serrano/mycofier.git .

Metagenomic Analysis of Fungal Diversity on Strawberry Plants and the Effect of Management Practices on the Fungal Community Structure of Aerial Organs

An amplicon metagenomic approach based on the ITS2 region of fungal rDNA was used to identify the composition of fungal communities associated with different strawberry organs (leaves, flowers, immature and mature fruits), grown on a farm using management practices that entailed the routine use of various chemical pesticides. ITS2 sequences clustered into 316 OTUs and Ascomycota was the dominant phyla (95.6%) followed by Basidiomycota (3.9%). Strawberry plants supported a high diversity of microbial organisms, but two genera, Botrytis and Cladosporium, were the most abundant, representing 70–99% of the relative abundance (RA) of all detected sequences. According to alpha and beta diversity analyses, strawberry organs displayed significantly different fungal communities with leaves having the most diverse fungal community, followed by flowers, and fruit. The interruption of chemical treatments for one month resulted in a significant modification in the structure of the fungal community of leaves and flowers while immature and mature fruit were not significantly affected. Several plant pathogens of other plant species, that would not be intuitively expected to be present on strawberry plants such as Erysiphe, were detected, while some common strawberry pathogens, such as Rhizoctonia, were less evident or absent.

The defensive role of foliar endophytic fungi for a South American tree

Fungal endophytes colonize living internal plant tissues without causing any visible symptoms of disease. Endophytic fungi associated with healthy leaves may play an important role in the protection of hosts against herbivores and pathogens. In this study, the diversity of foliar endophytic fungi (FEF) of the southern temperate tree Embothrium coccineum (Proteaceae), as well as their role in plant protection in nature was determined. Fungal endophytes were isolated from 40 asymptomatic leaves by the culture method for molecular identification of the 18S rRNA gene. A relationship between FEF frequency and plant protection was evaluated in juveniles of E. coccineum Fungal endophyte frequency was estimated using real-time PCR analyses to determine endophyte DNA content per plant. A total of 178 fungal isolates were identified, with sequence data revealing 34 different operational taxonomic units (OTUs). A few common taxa dominated the fungal endophyte community, whereas most taxa qualified as rare. A significant positive correlation between plant protection (evaluated in terms of percentage of leaf damage) and FEF frequency was found. Furthermore, in vitro confrontation assays indicated that FEF were able to inhibit the growth of fungal pathogens. The data showed a relatively high diversity of fungal endophytes associated with leaves of E. coccineum, and suggest a positive relationship between fungal endophyte frequencies in leaves and host protection in nature.

Highlights of DNA Barcoding in identification of salient microorganisms like fungi

Fungi, the second largest kingdom of eukaryotic life, are diverse and widespread. Fungi play a distinctive role in the production of different products on industrial scale, like fungal enzymes, antibiotics, fermented foods, etc., to give storage stability and improved health to meet major global challenges. To utilize algae perfectly for human needs, and to pave the way for getting a healthy relationship with fungi, it is important to identify them in a quick and robust manner with molecular-based identification system. So, there is a technique that aims to provide a well-organized method for species level identifications and to contribute powerfully to taxonomic and biodiversity research is DNA Barcoding. DNA Barcoding is generally achieved by the retrieval of a short DNA sequence - the 'barcode' - from a standard part of the genome and that barcode is then compared with a library of reference barcode sequences derived from individuals of known identity for identification.

Determining threshold values for barcoding fungi: lessons from Cortinarius (Basidiomycota), a highly diverse and widespread ectomycorrhizal genus

Different distance-based threshold selection approaches were used to assess and compare use of the internal transcribed spacer (ITS) region to distinguish among 901 Cortinarius species represented by >3000 collections. Sources of error associated with genetic markers and selection approaches were explored and evaluated using MOTUs from genus and lineage based-alignments. Our study indicates that 1%-2% more species can be distinguished by using the full-length ITS barcode as compared to either the ITS1 or ITS2 regions alone. Optimal threshold values for different picking approaches and genetic marker lengths inferred from a subset of species containing major lineages ranged from 97.0% to 99.5% sequence similarity using clustering optimization and UNITE SH, and from 1% to 2% sequence dissimilarity with CROP. Errors for the optimal cutoff ranged from 0% to 70%, and these can be reduced to a maximum of 22% when excluding species lacking a barcode gap. A threshold value of 99% is suitable for distinguishing species in the majority of lineages in the genus using the entire ITS region but only 90% of the species could be identified using just the ITS1 or ITS2 region. Prior identification of species, lacking barcode gaps and their subsequent separate analyses, maximized the accuracy of threshold approaches.

Ribosomal DNA and Plastid Markers Used to Sample Fungal and Plant Communities from Wetland Soils Reveals Complementary Biotas

Though the use of metagenomic methods to sample below-ground fungal communities is common, the use of similar methods to sample plants from their underground structures is not. In this study we use high throughput sequencing of the ribulose-bisphosphate carboxylase large subunit (rbcL) plastid marker to study the plant community as well as the internal transcribed spacer and large subunit ribosomal DNA (rDNA) markers to investigate the fungal community from two wetland sites. Observed community richness and composition varied by marker. The two rDNA markers detected complementary sets of fungal taxa and total fungal composition clustered according to primer rather than by site. The composition of the most abundant plants, however, clustered according to sites as expected. We suggest that future studies consider using multiple genetic markers, ideally generated from different primer sets, to detect a more taxonomically diverse suite of taxa compared with what can be detected by any single marker alone. Conclusions drawn from the presence of even the most frequently observed taxa should be made with caution without corroborating lines of evidence.

Diversity of Marine-Derived Fungal Cultures Exposed by DNA Barcodes: The Algorithm Matters

Marine fungi are an understudied group of eukaryotic microorganisms characterized by unresolved genealogies and unstable classification. Whereas DNA barcoding via the nuclear ribosomal internal transcribed spacer (ITS) provides a robust and rapid tool for fungal species delineation, accurate classification of fungi is often arduous given the large number of partial or unknown barcodes and misidentified isolates deposited in public databases. This situation is perpetuated by a paucity of cultivable fungal strains available for phylogenetic research linked to these data sets. We analyze ITS barcodes produced from a subsample (290) of 1781 cultured isolates of marine-derived fungi in the Bioresources Library located at the Australian Institute of Marine Science (AIMS). Our analysis revealed high levels of under-explored fungal diversity. The majority of isolates were ascomycetes including representatives of the subclasses Eurotiomycetidae, Hypocreomycetidae, Sordariomycetidae, Pleosporomycetidae, Dothideomycetidae, Xylariomycetidae and Saccharomycetidae. The phylum Basidiomycota was represented by isolates affiliated with the genera Tritirachium and Tilletiopsis. BLAST searches revealed 26 unknown OTUs and 50 isolates corresponding to previously uncultured, unidentified fungal clones. This study makes a significant addition to the availability of barcoded, culturable marine-derived fungi for detailed future genomic and physiological studies. We also demonstrate the influence of commonly used alignment algorithms and genetic distance measures on the accuracy and comparability of estimating Operational Taxonomic Units (OTUs) by the automatic barcode gap finder (ABGD) method. Large scale biodiversity screening programs that combine datasets using algorithmic OTU delineation pipelines need to ensure compatible algorithms have been used because the algorithm matters.

Assessing fungal community structure from mineral surfaces in Kartchner Caverns using multiplexed 454 pyrosequencing

Research on the distribution and structure of fungal communities in caves is lacking. Kartchner Caverns is a wet and mineralogically diverse carbonate cave located in an escarpment of Mississippian Escabrosa limestone in the Whetstone Mountains, Arizona, USA. Fungal diversity from speleothem and rock wall surfaces was examined with 454 FLX Titanium sequencing technology using the Internal Transcribed Spacer 1 as a fungal barcode marker. Fungal diversity was estimated and compared between speleothem and rock wall surfaces, and its variation with distance from the natural entrance of the cave was quantified. Effects of environmental factors and nutrient concentrations in speleothem drip water at different sample sites on fungal diversity were also examined. Sequencing revealed 2,219 fungal operational taxonomic units (OTUs) at the 95% similarity level. Speleothems supported a higher fungal richness and diversity than rock walls. However, community membership and the taxonomic distribution of fungal OTUs at the class level did not differ significantly between speleothems and rock walls. Both OTU richness and diversity decreased significantly with increasing distance from the natural cave entrance. Community membership and taxonomic distribution of fungal OTUs also differed significantly between the sampling sites closest to the entrance and those furthest away. There was no significant effect of temperature, CO2 concentration, or drip water nutrient concentration on fungal community structure on either speleothems or rock walls. Together, these results suggest that proximity to the natural entrance is a critical factor in determining fungal community structure on mineral surfaces in Kartchner Caverns.

Metabarcoding Analysis of Fungal Diversity in the Phyllosphere and Carposphere of Olive (Olea europaea)

The fungal diversity associated with leaves, flowers and fruits of olive (Olea europaea) was investigated in different phenological stages (May, June, October and December) using an implemented metabarcoding approach. It consisted of the 454 pyrosequencing of the fungal ITS2 region and the subsequent phylogenetic analysis of relevant genera along with validated reference sequences. Most sequences were identified up to the species level or were associated with a restricted number of related taxa enabling supported speculations regarding their biological role. Analyses revealed a rich fungal community with 195 different OTUs. Ascomycota was the dominating phyla representing 93.6% of the total number of detected sequences followed by unidentified fungi (3.6%) and Basidiomycota (2.8%). A higher level of diversity was revealed for leaves compared to flowers and fruits. Among plant pathogens the genus Colletotrichum represented by three species (C. godetiae syn. C. clavatum, C. acutatum s.s and C. karstii) was the most abundant on ripe fruits but it was also detected in other organs. Pseudocercospora cladosporioides was detected with a high frequency in all leaf samples and to a less extent in ripe fruits. A much lower relative frequency was revealed for Spilocaea oleagina and for other putative pathogens including Fusarium spp., Neofusicoccum spp., and Alternaria spp. Among non-pathogen taxa, Aureobasidium pullulans, the species complex of Cladosporium cladosporioides and Devriesia spp. were the most represented. This study highlights the existence of a complex fungal consortium including both phytopathogenic and potentially antagonistic microorganisms that can have a significant impact on olive productions.

DNA barcoding of fungi causing infections in humans and animals

Correct species identification is becoming increasingly important in clinical diagnostics. Till now, many mycological laboratories rely on conventional phenotypic identification. But this is slow and strongly operator-dependent. Therefore, to improve the quality of pathogen identification, rapid, reliable, and objective identification methods are essential. One of the most encouraging approaches is molecular barcoding using the internal transcribed spacer (ITS) of the rDNA, which is rapid, easily achievable, accurate, and applicable directly from clinical specimens. It relies on the comparison of a single ITS sequence with a curated reference database. The International Society for Human and Animal Mycology (ISHAM) working group for DNA barcoding has recently established such a database, focusing on the majority of human and animal pathogenic fungi (ISHAM-ITS, freely accessible at http://www.isham.org/ or directly from http://its.mycologylab.org). For some fungi the use of secondary barcodes may be necessary.

The use of genus-specific amplicon pyrosequencing to assess phytophthora species diversity using eDNA from soil and water in Northern Spain

Phytophthora is one of the most important and aggressive plant pathogenic genera in agriculture and forestry. Early detection and identification of its pathways of infection and spread are of high importance to minimize the threat they pose to natural ecosystems. eDNA was extracted from soil and water from forests and plantations in the north of Spain. Phytophthora-specific primers were adapted for use in high-throughput Sequencing (HTS). Primers were tested in a control reaction containing eight Phytophthora species and applied to water and soil eDNA samples from northern Spain. Different score coverage threshold values were tested for optimal Phytophthora species separation in a custom-curated database and in the control reaction. Clustering at 99% was the optimal criteria to separate most of the Phytophthora species. Multiple Molecular Operational Taxonomic Units (MOTUs) corresponding to 36 distinct Phytophthora species were amplified in the environmental samples. Pyrosequencing of amplicons from soil samples revealed low Phytophthora diversity (13 species) in comparison with the 35 species detected in water samples. Thirteen of the MOTUs detected in rivers and streams showed no close match to sequences in international sequence databases, revealing that eDNA pyrosequencing is a useful strategy to assess Phytophthora species diversity in natural ecosystems.

Complementary amplicon-based genomic approaches for the study of fungal communities in humans

Recent studies highlight the importance of intestinal fungal microbiota in the development of human disease. Infants, in particular, are an important population in which to study intestinal microbiomes because microbial community structure and dynamics during this formative window of life have the potential to influence host immunity and metabolism. When compared to bacteria, much less is known about the early development of human fungal communities, owing partly to their lower abundance and the relative lack of established molecular and taxonomic tools for their study. Herein, we describe the development, validation, and use of complementary amplicon-based genomic strategies to characterize infant fungal communities and provide quantitative information about Candida, an important fungal genus with respect to intestinal colonization and human disease. Fungal communities were characterized from 11 infant fecal samples using primers that target the internal transcribed spacer (ITS) 2 locus, a region that provides taxonomic discrimination of medically relevant fungi. Each sample yielded an average of 27,553 fungal sequences and Candida albicans was the most abundant species identified by sequencing and quantitative PCR (qPCR). Low numbers of Candida krusei and Candida parapsilosis sequences were observed in several samples, but their presence was detected by species-specific qPCR in only one sample, highlighting a challenge inherent in the study of low-abundance organisms. Overall, the sequencing results revealed that infant fecal samples had fungal diversity comparable to that of bacterial communities in similar-aged infants, which correlated with the relative abundance of C. albicans. We conclude that targeted sequencing of fungal ITS2 amplicons in conjunction with qPCR analyses of specific fungi provides an informative picture of fungal community structure in the human intestinal tract. Our data suggests that the infant intestine harbors diverse fungal species and is consistent with prior culture-based analyses showing that the predominant fungus in the infant intestine is C. albicans.