Sample Preparation for Fungal Community Analysis by High-Throughput Sequencing of Barcode Amplicons.

Fungal species participate in vast numbers of processes in the landscape around us. However, their cryptic mycelial growth, inside various substrates and in highly diverse species assemblages, has been a major obstacle to thorough analysis of fungal communities, hampering exhaustive description of the fungal kingdom. Technological developments allowing rapid, high-throughput sequencing of mixed communities from many samples at once are currently having a tremendous impact in fungal community ecology. Universal DNA extraction followed by amplification and sequencing of fungal species-level barcodes such as the nuclear internal transcribed spacer (ITS) region now enables identification and relative quantification of fungal community members across well-replicated experimental settings.Here, we present the sample preparation procedure presently used in our laboratory for fungal community analysis by high-throughput sequencing of amplified ITS2 markers. We focus on the procedure optimized for studies of total fungal communities in humus-rich soils, wood, and litter. However, this procedure can be applied to other sample types and markers. We focus on the laboratory-based part of sample preparation, i.e., the procedure from the point where samples enter the laboratory until amplicons are submitted for sequencing. Our procedure comprises four main parts: (1) universal DNA extraction, (2) optimization of PCR conditions, (3) production of tagged ITS amplicons, and (4) preparation of the multiplexed amplicon pool to be sequenced. The presented procedure is independent of the specific high-throughput sequencing technology used, which makes it highly versatile.

Comparative analyses of the Hymenoscyphus fraxineus and Hymenoscyphus albidus genomes reveals potentially adaptive differences in secondary metabolite and transposable element repertoires.

BACKGROUND: The dieback epidemic decimating common ash (Fraxinus excelsior) in Europe is caused by the invasive fungus Hymenoscyphus fraxineus. In this study we analyzed the genomes of H. fraxineus and H. albidus, its native but, now essentially displaced, non-pathogenic sister species, and compared them with several other members of Helotiales. The focus of the analyses was to identify signals in the genome that may explain the rapid establishment of H. fraxineus and displacement of H. albidus. RESULTS: The genomes of H. fraxineus and H. albidus showed a high level of synteny and identity. The assembly of H. fraxineus is 13 Mb longer than that of H. albidus', most of this difference can be attributed to higher dispersed repeat content (i.e. transposable elements [TEs]) in H. fraxineus. In general, TE families in H. fraxineus showed more signals of repeat-induced point mutations (RIP) than in H. albidus, especially in Long-terminal repeat (LTR)/Copia and LTR/Gypsy elements. Comparing gene family expansions and 1:1 orthologs, relatively few genes show signs of positive selection between species. However, several of those did appeared to be associated with secondary metabolite genes families, including gene families containing two of the genes in the H. fraxineus-specific, hymenosetin biosynthetic gene cluster (BGC). CONCLUSION: The genomes of H. fraxineus and H. albidus show a high degree of synteny, and are rich in both TEs and BGCs, but the genomic signatures also indicated that H. albidus may be less well equipped to adapt and maintain its ecological niche in a rapidly changing environment.

A tipping point in carbon storage when forest expands into tundra is related to mycorrhizal recycling of nitrogen.

Tundra ecosystems are global belowground sinks for atmospheric CO2 . Ongoing warming-induced encroachment by shrubs and trees risks turning this sink into a CO2 source, resulting in a positive feedback on climate warming. To advance mechanistic understanding of how shifts in mycorrhizal types affect long-term carbon (C) and nitrogen (N) stocks, we studied small-scale soil depth profiles of fungal communities and C-N dynamics across a subarctic-alpine forest-heath vegetation gradient. Belowground organic stocks decreased abruptly at the transition from heath to forest, linked to the presence of certain tree-associated ectomycorrhizal fungi that contribute to decomposition when mining N from organic matter. In contrast, ericoid mycorrhizal plants and fungi were associated with organic matter accumulation and slow decomposition. If climatic controls on arctic-alpine forest lines are relaxed, increased decomposition will likely outbalance increased plant productivity, decreasing the overall C sink capacity of displaced tundra.

The conifer root rot pathogens Heterobasidion irregulare and Heterobasidion occidentale employ different strategies to infect Norway spruce.

Heterobasidion irregulare and H. occidentale are two closely related conifer root rot pathogens in the H. annosum sensu lato (s.l.) species complex. The two species H. irregulare and H. occidentale have different host preference with pine and non-pine tree species favored, respectively. The comparison of transcriptomes of H. irregulare and H. occidentale growing in Norway spruce bark, a susceptible host non-native to North America, showed large differences in gene expression. Heterobasidion irregulare induced more genes involved in detoxification of host compounds and in production of secondary metabolites, while the transcriptome induced in H. occidentale was more oriented towards carbohydrate degradation. Along with their separated evolutionary history, the difference might be driven by their host preferences as indicated by the differentially expressed genes enriched in particular Gene Ontology terms.

The mycoparasitic fungus Clonostachys rosea responds with both common and specific gene expression during interspecific interactions with fungal prey.

Clonostachys rosea is a necrotrophic mycoparasitic fungus, used for biological control of plant pathogenic fungi. A better understanding of the underlying mechanisms resulting in successful biocontrol is important for knowledge-based improvements of the application and use of biocontrol in agricultural production systems. Transcriptomic analyses revealed that C. rosea responded with both common and specific gene expression during interactions with the fungal prey species Botrytis cinerea and Fusarium graminearum. Genes predicted to encode proteins involved in membrane transport, biosynthesis of secondary metabolites and carbohydrate-active enzymes were induced during the mycoparasitic attack. Predicted major facilitator superfamily (MFS) transporters constituted 54% of the induced genes, and detailed phylogenetic and evolutionary analyses showed that a majority of these genes belonged to MFS gene families evolving under selection for increased paralog numbers, with predicted functions in drug resistance and transport of carbohydrates and small organic compounds. Sequence analysis of MFS transporters from family 2.A.1.3.65 identified rapidly evolving loop regions forming the entry to the transport tunnel, indicating changes in substrate specificity as a target for selection. Deletion of the MFS transporter gene mfs464 resulted in mutants with increased growth inhibitory activity against F. graminearum, providing evidence for a function in interspecific fungal interactions. In summary, we show that the mycoparasite C. rosea can distinguish between fungal prey species and modulate its transcriptomic responses accordingly. Gene expression data emphasize the importance of secondary metabolites in mycoparasitic interactions.

Draft genome of the brown-rot fungus Fomitopsis pinicola GR9-4.

Basidiomycete brown-rot fungi have a huge importance for wood decomposition and thus the global carbon cycle. Here, we present the genome sequence of Fomitopsis pinicola GR9-4 which represent different F. pinicola clade than the previously sequenced North American isolate FP-58527 SS1. The genome was sequenced by using a paired-end sequence library of Illumina and a 2.5k and 5k mate-pair library (ABI SOLiD). The final assembly adds up to a size of 45 Mb (including gaps between contigs), with a GC-content of 56%. The gene prediction resulted in 13,888 gene models. The genome sequence will be used as a basis for understanding population genomics, genome-wide association studies and wood decay mechanisms of this brown-rot fungus.

Sample Preparation for Fungal Community Analysis by High-Throughput Sequencing of Barcode Amplicons.

Fungal species participate in vast numbers of processes in the landscape around us. However, their often cryptic growth, inside various substrates and in highly diverse species assemblages, has been a major obstacle to thorough analysis of fungal communities, hampering exhaustive description of the fungal kingdom. Recent technological developments allowing rapid, high-throughput sequencing of mixed communities from many samples at once are currently having a tremendous impact in fungal community ecology. Universal DNA extraction followed by amplification and sequencing of fungal species-level barcodes such as the nuclear internal transcribed spacer (ITS) region now enable identification and relative quantification of fungal community members across well-replicated experimental settings. Here, we present the sample preparation procedure presently used in our laboratory for fungal community analysis by high-throughput sequencing of amplified ITS2 markers. We focus on the procedure optimized for studies of total fungal communities in humus-rich soils, wood, and litter. However, this procedure can be applied to other sample types and markers. We focus on the laboratory-based part of sample preparation, that is, the procedure from the point where samples enter the laboratory until amplicons are submitted for sequencing. Our procedure comprises four main parts: (1) universal DNA extraction, (2) optimization of PCR conditions, (3) production of tagged ITS amplicons, and (4) preparation of the multiplexed amplicon mix to be sequenced. The presented procedure is independent of the specific high-throughput sequencing technology used, which makes it highly versatile.

Genus-Specific Primers for Study of Fusarium Communities in Field Samples.

Fusarium is a large and diverse genus of fungi of great agricultural and economic importance, containing many plant pathogens and mycotoxin producers. To date, high-throughput sequencing of Fusarium communities has been limited by the lack of genus-specific primers targeting regions with high discriminatory power at the species level. In the present study, we evaluated two Fusarium-specific primer pairs targeting translation elongation factor 1 (TEF1). We also present the new primer pair Fa+7/Ra+6. Mock Fusarium communities reflecting phylogenetic diversity were used to evaluate the accuracy of the primers in reflecting the relative abundance of the species. TEF1 amplicons were subjected to 454 high-throughput sequencing to characterize Fusarium communities. Field samples from soil and wheat kernels were included to test the method on more-complex material. For kernel samples, a single PCR was sufficient, while for soil samples, nested PCR was necessary. The newly developed primer pairs Fa+7/Ra+6 and Fa/Ra accurately reflected Fusarium species composition in mock DNA communities. In field samples, 47 Fusarium operational taxonomic units were identified, with the highest Fusarium diversity in soil. The Fusarium community in soil was dominated by members of the Fusarium incarnatum-Fusarium equiseti species complex, contradicting findings in previous studies. The method was successfully applied to analyze Fusarium communities in soil and plant material and can facilitate further studies of Fusarium ecology.

Transcriptional Responses Associated with Virulence and Defence in the Interaction between Heterobasidion annosum s.s. and Norway Spruce.

Heterobasidion annosum sensu lato is a serious pathogen causing root and stem rot to conifers in the northern hemisphere and rendering the timber defective for sawing and pulping. In this study we applied next-generation sequencing to i) identify transcriptional responses unique to Heterobasidion-inoculated Norway spruce and ii) investigate the H. annosum transcripts to identify putative virulence factors. To address these objectives we wounded or inoculated 30-year-old Norway spruce clones with H. annosum and 454-sequenced the transcriptome of the interaction at 0, 5 and 15 days post inoculation. The 491,860 high-quality reads were de novo assembled and the relative expression was analysed. Overall, very few H. annosum transcripts were represented in our dataset. Three delta-12 fatty acid desaturase transcripts and one Clavaminate synthase-like transcript, both associated with virulence in other pathosystems, were found among the significantly induced transcripts. The analysis of the Norway spruce transcriptional responses produced a handful of differentially expressed transcripts. Most of these transcripts originated from genes known to respond to H. annosum. However, three genes that had not previously been reported to respond to H. annosum showed specific induction to inoculation: an oxophytodienoic acid-reductase (OPR), a beta-glucosidase and a germin-like protein (GLP2) gene. Even in a small data set like ours, five novel highly expressed Norway spruce transcripts without significant alignment to any previously annotated protein in Genbank but present in the P. abies (v1.0) gene catalogue were identified. Their expression pattern suggests a role in defence. Therefore a more complete survey of the transcriptional responses in the interactions between Norway spruce and its major pathogen H. annosum would probably provide a better understanding of gymnosperm defence than accumulated until now.

Insights on the evolution of mycoparasitism from the genome of Clonostachys rosea.

Clonostachys rosea is a mycoparasitic fungus that can control several important plant diseases. Here, we report on the genome sequencing of C. rosea and a comparative genome analysis, in order to resolve the phylogenetic placement of C. rosea and to study the evolution of mycoparasitism as a fungal lifestyle. The genome of C. rosea is estimated to 58.3 Mb, and contains 14,268 predicted genes. A phylogenomic analysis shows that C. rosea clusters as sister taxon to plant pathogenic Fusarium species, with mycoparasitic/saprotrophic Trichoderma species in an ancestral position. A comparative analysis of gene family evolution reveals several distinct differences between the included mycoparasites. Clonostachys rosea contains significantly more ATP-binding cassette (ABC) transporters, polyketide synthases, cytochrome P450 monooxygenases, pectin lyases, glucose-methanol-choline oxidoreductases, and lytic polysaccharide monooxygenases compared with other fungi in the Hypocreales. Interestingly, the increase of ABC transporter gene number in C. rosea is associated with phylogenetic subgroups B (multidrug resistance proteins) and G (pleiotropic drug resistance transporters), whereas an increase in subgroup C (multidrug resistance-associated proteins) is evident in Trichoderma virens. In contrast with mycoparasitic Trichoderma species, C. rosea contains very few chitinases. Expression of six group B and group G ABC transporter genes was induced in C. rosea during exposure to the Fusarium mycotoxin zearalenone, the fungicide Boscalid or metabolites from the biocontrol bacterium Pseudomonas chlororaphis. The data suggest that tolerance toward secondary metabolites is a prominent feature in the biology of C. rosea.

Finding the founder of Stockholm - a kinship study based on Y-chromosomal, autosomal and mitochondrial DNA.

Historical records claim that Birger Magnusson (died 1266), famous regent of Sweden and the founder of Stockholm, was buried in Varnhem Abbey in Västergötland. After being lost for centuries, his putative grave was rediscovered during restoration work in the 1920s. Morphological analyses of the three individuals in the grave concluded that the older male, the female and the younger male found in the grave were likely to be Birger, his second wife Mechtild of Holstein and his son Erik from a previous marriage. More recent evaluations of the data from the 1920s seriously questioned these conclusions, ultimately leading to the reopening and reexamination of the grave in 2002. Ancient DNA-analyses were performed to investigate if the relationship between the three individuals matched what we would expect if the individuals were Birger, Erik and Mechtild. We used pyrosequencing of Y-chromosomal and autosomal SNPs and compared the results with haplogroup frequencies of modern Swedes to investigate paternal relations. Possible maternal kinship was investigated by deep FLX-sequencing of overlapping mtDNA amplicons. The authenticity of the sequences was examined using data from independent extractions, massive clonal data, the c-statistics, and real-time quantitative data. We show that the males carry the same Y-chromosomal haplogroup and thus we cannot reject a father-son type of relation. Further, as shown by the mtDNA analyses, none of the individuals are maternally related. We conclude that the graves indeed belong to Birger, Erik and Mechtild, or to three individuals with the exact same kind of biological relatedness.