Fungal cross-talk: an integrated approach to study distance communication.

Despite the interest on fungi as eukaryotic model systems, the molecular mechanisms regulating the fungal non-self-recognition at a distance have not been studied so far. This paper investigates the molecular mechanisms regulating the cross-talk at a distance between two filamentous fungi, Trichoderma gamsii and Fusarium graminearum which establish a mycoparasitic interaction where T. gamsii and F. graminearum play the roles of mycoparasite and prey, respectively. In the present work, we use an integrated approach involving dual culture tests, comparative genomics and transcriptomics to investigate the fungal interaction before contact ('sensing phase'). Dual culture tests demonstrate that growth rate of F. graminearum accelerates in presence of T. gamsii at the sensing phase. T. gamsii up-regulates the expression of a ferric reductase involved in iron acquisition, while F. graminearum up-regulates the expression of genes coding for transmembrane transporters and killer toxins. At the same time, T. gamsii decreases the level of extracellular interaction by down-regulating genes coding for hydrolytic enzymes acting on fungal cell wall (chitinases). Given the importance of fungi as eukaryotic model systems and the ever-increasing genomic resources available, the integrated approach hereby presented can be applied to other interactions to deepen the knowledge on fungal communication at a distance.

Transcriptional responses in developing lesions of European common ash (Fraxinus excelsior) reveal genes responding to infection by Hymenoscyphus fraxineus.

BACKGROUND: With the expanding ash dieback epidemic that has spread across the European continent, an improved functional understanding of the disease development in afflicted hosts is needed. The study investigated whether differences in necrosis extension between common ash (Fraxinus excelsior) trees with different levels of susceptibility to the fungus Hymenoscyphus fraxineus are associated with, and can be explained by, the differences in gene expression patterns. We inoculated seemingly healthy branches of each of two resistant and susceptible ash genotypes with H. fraxineus grown in a common garden. RESULTS: Ten months after the inoculation, the length of necrosis on the resistant genotypes were shorter than on the susceptible genotypes. RNA sequencing of bark samples collected at the border of necrotic lesions and from healthy tissues distal to the lesion revealed relatively limited differences in gene expression patterns between susceptible and resistant genotypes. At the necrosis front, only 138 transcripts were differentially expressed between the genotype categories while 1082 were differentially expressed in distal, non-symptomatic tissues. Among these differentially expressed genes, several genes in the mevalonate (MVA) and iridoid pathways were found to be co-regulated, possibly indicating increased fluxes through these pathways in response to H. fraxineus. Comparison of transcriptional responses of symptomatic and non-symptomatic ash in a controlled greenhouse experiment revealed a relatively small set of genes that were differentially and concordantly expressed in both studies. This gene-set included the rate-limiting enzyme in the MVA pathway and a number of transcription factors. Furthermore, several of the concordantly expressed candidate genes show significant similarity to genes encoding players in the abscisic acid- or Jasmonate-signalling pathways. CONCLUSIONS: A set of candidate genes, concordantly expressed between field and greenhouse experiments, was identified. The candidates are associated with hormone signalling and specialized metabolite biosynthesis pathways indicating the involvement of these pathways in the response of the host to infection by H. fraxineus.

Optimized metabarcoding with Pacific biosciences enables semi-quantitative analysis of fungal communities.

Recent studies have questioned the use of high-throughput sequencing of the nuclear ribosomal internal transcribed spacer (ITS) region to derive a semi-quantitative representation of fungal community composition. However, comprehensive studies that quantify biases occurring during PCR and sequencing of ITS amplicons are still lacking. We used artificially assembled communities consisting of 10 ITS-like fragments of varying lengths and guanine-cytosine (GC) contents to evaluate and quantify biases during PCR and sequencing with Illumina MiSeq, PacBio RS II and PacBio Sequel I technologies. Fragment length variation was the main source of bias in observed community composition relative to the template, with longer fragments generally being under-represented for all sequencing platforms. This bias was three times higher for Illumina MiSeq than for PacBio RS II and Sequel I. All 10 fragments in the artificial community were recovered when sequenced with PacBio technologies, whereas the three longest fragments (> 447 bases) were lost when sequenced with Illumina MiSeq. Fragment length bias also increased linearly with increasing number of PCR cycles but could be mitigated by optimization of the PCR setup. No significant biases related to GC content were observed. Despite lower sequencing output, PacBio sequencing was better able to reflect the community composition of the template than Illumina MiSeq sequencing.

Changes in fungal communities along a boreal forest soil fertility gradient.

Boreal forests harbour diverse fungal communities with decisive roles in decomposition and plant nutrition. Although changes in boreal plant communities along gradients in soil acidity and nitrogen (N) availability are well described, less is known about how fungal taxonomic and functional groups respond to soil fertility factors. We analysed fungal communities in humus and litter from 25 Swedish old-growth forests, ranging from N-rich Picea abies stands to acidic and N-poor Pinus sylvestris stands. 454-pyrosequencing of ITS2 amplicons was used to analyse community composition, and biomass was estimated by ergosterol analysis. Fungal community composition was significantly related to soil fertility at the levels of species, genera/orders and functional groups. Ascomycetes dominated in less fertile forests, whereas basidiomycetes increased in abundance in more fertile forests, both in litter and humus. The relative abundance of mycorrhizal fungi in the humus layer remained high even in the most fertile soils. Tolerance to acidity and nitrogen deficiency seems to be of greater importance than plant carbon (C) allocation patterns in determining responses of fungal communities to soil fertility, in old-growth boreal forests.

Population genomic analyses reveal possible drivers of population divergence.

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. ), Saccharomyces uvarum (Almeida et al. ), Coccidioides spp. (Neafsey et al. ) and Cryptococcus gatti (Engelthaler et al. ). In this issue of Molecular Ecology, Branco et al. () apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.

Intronic and plasmid-derived regions contribute to the large mitochondrial genome sizes of Agaricomycetes.

Sizes of mitochondrial genomes vary extensively between fungal species although they typically contain a conserved set of core genes. We have characterised the mitochondrial genome of the conifer root rot pathogen Heterobasidion irregulare and compared the size, gene content and structure of 20 Basidiomycete mitochondrial genomes. The mitochondrial genome of H. irregulare was 114, 193 bp and contained a core set of 15 protein coding genes, two rRNA genes and 26 tRNA genes. In addition, we found six non-conserved open reading frames (ORFs) and four putative plasmid genes clustered in three separate regions together with 24 introns and 14 intronic homing endonuclease genes, unequally spread across seven of the core genes. The size differences among the 20 Basidiomycetes can largely be explained by length variation of intergenic regions and introns. The Agaricomycetes contained the nine largest mitochondrial genomes in the Basidiomycete group and Agaricomycete genomes are significantly (p < 0.001) larger than the other Basidiomycetes. A feature of the Agaricomycete mitochondrial genomes in this study was the simultaneous occurrence of putative plasmid genes and non-conserved ORFs, with Cantharellus cibarius as only exception, where no non-conserved ORF was identified. This indicates a mitochondrial plasmid origin of the non-conserved ORFs or increased mitochondrial genome dynamics of species harbouring mitochondrial plasmids. We hypothesise that two independent factors are the driving forces for large mitochondrial genomes: the homing endonuclease genes in introns and integration of plasmid DNA.

Partial genetic turnover in neandertals: continuity in the East and population replacement in the West.

Remarkably little is known about the population-level processes leading up to the extinction of the neandertal. To examine this, we use mitochondrial DNA sequences from 13 neandertal individuals, including a novel sequence from northern Spain, to examine neandertal demographic history. Our analyses indicate that recent western European neandertals (<48 kyr) constitute a tightly defined group with low mitochondrial genetic variation in comparison with both eastern and older (>48 kyr) European neandertals. Using control region sequences, Bayesian demographic simulations provide higher support for a model of population fragmentation followed by separate demographic trajectories in subpopulations over a null model of a single stable population. The most parsimonious explanation for these results is that of a population turnover in western Europe during early Marine Isotope Stage 3, predating the arrival of anatomically modern humans in the region.

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex.

Introduction: Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species. Methods: In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers. Results and Discussion: Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

A pile of pipelines: An overview of the bioinformatics software for metabarcoding data analyses.

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects.

Comparative genomics highlights the importance of drug efflux transporters during evolution of mycoparasitism in Clonostachys subgenus Bionectria (Fungi, Ascomycota, Hypocreales).

Various strains of the mycoparasitic fungal species Clonostachys rosea are used commercially as biological control agents for the control of fungal plant diseases in agricultural crop production. Further improvements of the use and efficacy of C. rosea in biocontrol require a mechanistic understanding of the factors that determines the outcome of the interaction between C. rosea and plant pathogenic fungi. Here, we determined the genome sequences of 11 Clonostachys strains, representing five species in Clonostachys subgenus Bionectria, and performed a comparative genomic analysis with the aim to identify gene families evolving under selection for gene gains or losses. Several gene families predicted to encode proteins involved in biosynthesis of secondary metabolites, including polyketide synthases, nonribosomal peptide syntethases and cytochrome P450s, evolved under selection for gene gains (p ≤ .05) in the Bionectria subgenus lineage. This was accompanied with gene copy number increases (p ≤ .05) in ATP-binding cassette (ABC) transporters and major facilitator superfamily (MFS) transporters predicted to contribute to drug efflux. Most Clonostachys species were also characterized by high numbers of auxiliary activity (AA) family 9 lytic polysaccharide monooxygenases, AA3 glucose-methanol-choline oxidoreductases and additional carbohydrate-active enzyme gene families with putative activity (or binding) towards xylan and rhamnose/pectin substrates. Particular features of the C. rosea genome included expansions (p ≤ .05) of the ABC-B4 multidrug resistance transporters, the ABC-C5 multidrug resistance-related transporters and the 2.A.1.3 drug:H + antiporter-2 MFS drug resistance transporters. The ABC-G1 pleiotropic drug resistance transporter gene abcG6 in C. rosea was induced (p ≤ .009) by exposure to the antifungal Fusarium mycotoxin zearalenone (1121-fold) and various fungicides. Deletion of abcG6 resulted in mutants with reduced (p < .001) growth rates on media containing the fungicides boscalid, fenhexamid and iprodione. Our results emphasize the role of biosynthesis of, and protection against, secondary metabolites in Clonostachys subgenus Bionectria.

Production of ectomycorrhizal mycelium peaks during canopy closure in Norway spruce forests.

*Here, species composition and biomass production of actively growing ectomycorrhizal (EM) mycelia were studied over the rotation period of managed Norway spruce (Picea abies) stands in south-western Sweden. *The EM mycelia were collected using ingrowth mesh bags incubated in the forest soil during one growing season. Fungal biomass was estimated by ergosterol analysis and the EM species were identified by 454 sequencing of internal transcribed spacer (ITS) amplicons. Nutrient availability and the fungal biomass in soil samples were also estimated. *Biomass production peaked in young stands (10-30 yr old) before the first thinning phase. Tylospora fibrillosa dominated the EM community, especially in these young stands, where it constituted 80% of the EM amplicons derived from the mesh bags. Species richness increased in older stands. *The establishment of EM mycelial networks in young Norway spruce stands requires large amounts of carbon, while much less is needed to sustain the EM community in older stands. The variation in EM biomass production over the rotation period has implications for carbon sequestration rates in forest soils.

Ancient DNA reveals traces of Iberian Neolithic and Bronze Age lineages in modern Iberian horses.

Multiple geographical regions have been proposed for the domestication of Equus caballus. It has been suggested, based on zooarchaeological and genetic analyses that wild horses from the Iberian Peninsula were involved in the process, and the overrepresentation of mitochondrial D1 cluster in modern Iberian horses supports this suggestion. To test this hypothesis, we analysed mitochondrial DNA from 22 ancient Iberian horse remains belonging to the Neolithic, the Bronze Age and the Middle Ages, against previously published sequences. Only the medieval Iberian sequence appeared in the D1 group. Neolithic and Bronze Age sequences grouped in other clusters, one of which (Lusitano group C) is exclusively represented by modern horses of Iberian origin. Moreover, Bronze Age Iberian sequences displayed the lowest nucleotide diversity values when compared with modern horses, ancient wild horses and other ancient domesticates using nonparametric bootstrapping analyses. We conclude that the excessive clustering of Bronze Age horses in the Lusitano group C, the observed nucleotide diversity and the local continuity from wild Neolithic Iberian to modern Iberian horses, could be explained by the use of local wild mares during an early Iberian domestication or restocking event, whereas the D1 group probably was introduced into Iberia in later historical times.

Natural variation of root lesion nematode antagonism in the biocontrol fungus Clonostachys rosea and identification of biocontrol factors through genome-wide association mapping.

Biological control is a promising approach to reduce plant diseases caused by nematodes to ensure high productivity in agricultural production. Large-scale analyses of genetic variation in fungal species used for biocontrol can generate knowledge regarding interaction mechanisms that can improve efficacy of biocontrol applications. In this study, we performed a genome-wide association study (GWAS) for in vitro antagonism against the root lesion nematode Pratylenchus penetrans in 53 previously genome re-sequenced strains of the biocontrol fungus Clonostachys rosea. Nematode mortality in C. rosea potato dextrose broth (PDB) culture filtrates was highly variable and showed continuous variation (p < .001) between strains, indicating a polygenic inheritance. Twenty-one strains produced culture filtrates with higher (p ≤ .05) nematode mortality compared with the PDB control treatment, while ten strains lowered (p ≤ .05) the mortality. The difference in in vitro antagonism against P. penetrans correlated with antagonism against the soybean cyst nematode Heterodera glycines, indicating lack of host specificity in C. rosea. An empirical Bayesian multiple hypothesis testing approach identified 279 single nucleotide polymorphism markers significantly (local false sign rate < 10-10) associated with the trait. Genes present in the genomic regions associated with nematicidal activity included several membrane transporters, a chitinase and genes encoding proteins predicted to biosynthesize secondary metabolites. Gene deletion strains of the predicted nonribosomal peptide synthetase genes nps4 and nps5 were generated and showed increased (p ≤ .001) fungal growth and conidiation rates compared to the wild type. Deletion strains also exhibited reduced (p < .001) nematicidal activity and reduced (p ≤ .05) biocontrol efficacy against nematode root disease and against fusarium foot rot on wheat. In summary, we show that the GWAS approach can be used to identify biocontrol factors in C. rosea, specifically the putative nonribosomal peptide synthetases NPS4 and NPS5.

Estimating the Fitness Effect of Deleterious Mutations During the Two Phases of the Life Cycle: A New Method Applied to the Root-Rot Fungus Heterobasidion parviporum.

Many eukaryote species, including taxa such as fungi or algae, have a lifecycle with substantial haploid and diploid phases. A recent theoretical model predicts that such haploid-diploid lifecycles are stable over long evolutionary time scales when segregating deleterious mutations have stronger effects in homozygous diploids than in haploids and when they are partially recessive in heterozygous diploids. The model predicts that effective dominance-a measure that accounts for these two effects-should be close to 0.5 in these species. It also predicts that diploids should have higher fitness than haploids on average. However, an appropriate statistical framework to conjointly investigate these predictions is currently lacking. In this study, we derive a new quantitative genetic model to test these predictions using fitness data of two haploid parents and their diploid offspring, and genome-wide genetic distance between haploid parents. We apply this model to the root-rot basidiomycete fungus Heterobasidion parviporum-a species where the heterokaryotic (equivalent to the diploid) phase is longer than the homokaryotic (haploid) phase. We measured two fitness-related traits (mycelium growth rate and the ability to degrade wood) in both homokaryons and heterokaryons, and we used whole-genome sequencing to estimate nuclear genetic distance between parents. Possibly due to a lack of power, we did not find that deleterious mutations were recessive or more deleterious when expressed during the heterokaryotic phase. Using this model to compare effective dominance among haploid-diploid species where the relative importance of the two phases varies should help better understand the evolution of haploid-diploid life cycles.

Out in the Cold: Identification of Genomic Regions Associated With Cold Tolerance in the Biocontrol Fungus Clonostachys rosea Through Genome-Wide Association Mapping.

There is an increasing importance for using biocontrol agents in combating plant diseases sustainably and in the long term. As large scale genomic sequencing becomes economically viable, the impact of single nucleotide polymorphisms (SNPs) on biocontrol-associated phenotypes can be easily studied across entire genomes of fungal populations. Here, we improved a previously reported genome assembly of the biocontrol fungus Clonostachys rosea strain IK726 using the PacBio sequencing platform, which resulted in a total genome size of 70.7 Mbp and 21,246 predicted genes. We further performed whole-genome re-sequencing of 52 additional C. rosea strains isolated globally using Illumina sequencing technology, in order to perform genome-wide association studies in conditions relevant for biocontrol activity. One such condition is the ability to grow at lower temperatures commonly encountered in cryic or frigid soils in temperate regions, as these will be prevalent for protecting growing crops in temperate climates. Growth rates at 10°C on potato dextrose agar of the 53 sequenced strains of C. rosea were measured and ranged between 0.066 and 0.413 mm/day. Performing a genome wide association study, a total of 1,478 SNP markers were significantly associated with the trait and located in 227 scaffolds, within or close to (< 1000 bp distance) 265 different genes. The predicted gene products included several chaperone proteins, membrane transporters, lipases, and proteins involved in chitin metabolism with possible roles in cold tolerance. The data reported in this study provides a foundation for future investigations into the genetic basis for cold tolerance in fungi, with important implications for biocontrol.

A genome-wide association study identifies genomic regions for virulence in the non-model organism Heterobasidion annosum s.s.

The dense single nucleotide polymorphisms (SNP) panels needed for genome wide association (GWA) studies have hitherto been expensive to establish and use on non-model organisms. To overcome this, we used a next generation sequencing approach to both establish SNPs and to determine genotypes. We conducted a GWA study on a fungal species, analysing the virulence of Heterobasidion annosum s.s., a necrotrophic pathogen, on its hosts Picea abies and Pinus sylvestris. From a set of 33,018 single nucleotide polymorphisms (SNP) in 23 haploid isolates, twelve SNP markers distributed on seven contigs were associated with virulence (P<0.0001). Four of the contigs harbour known virulence genes from other fungal pathogens and the remaining three harbour novel candidate genes. Two contigs link closely to virulence regions recognized previously by QTL mapping in the congeneric hybrid H. irregulare × H. occidentale. Our study demonstrates the efficiency of GWA studies for dissecting important complex traits of small populations of non-model haploid organisms with small genomes.

A 2nd generation linkage map of Heterobasidion annosum s.l. based on in silico anchoring of AFLP markers.

In this study, we present a 2(nd) generation genetic linkage map of a cross between the North American species Heterobasidion irregulare and H. occidentale, based on the alignment of the previously published 1(st) generation map to the parental genomes. We anchored 216 of the original 308 AFLP markers to their respective restriction sites using an in silico-approach. The map resolution was improved by adding 146 sequence-tagged microsatellite markers and 39 sequenced gene markers. The new markers confirmed the positions of the anchored AFLP markers, fused the original 39 linkage groups together into 17, and fully expanded 12 of these to single groups covering entire chromosomes. Map coverage of the genome increased from 55.3% to 92.8%, with 96.3% of 430 markers collinearly aligned with the genome sequence. The anchored map also improved the H. irregulare assembly considerably. It identified several errors in scaffold arrangements and assisted in reducing the total number of major scaffolds from 18 to 15. This denser, more comprehensive map allowed sequence-based mapping of three intersterility loci and one mating type locus. This demonstrates the possibility to utilize an in silico procedure to convert anonymous markers into sequence-tagged ones, as well as the power of a sequence-anchored linkage map and its usefulness in the assembly of a whole genome sequence.

New primers to amplify the fungal ITS2 region--evaluation by 454-sequencing of artificial and natural communities.

With recent methodological advances, molecular markers are increasingly used for semi-quantitative analyses of fungal communities. The aim to preserve quantitative relationships between genotypes through PCR places new demands on primers to accurately match target sites and provide short amplicons. The internal transcribed spacer (ITS) region of the ribosome encoding genes is a commonly used marker for many fungal groups. Here, we describe three new primers - fITS7, gITS7 and fITS9, which may be used to amplify the fungal ITS2 region by targeting sites in the 5.8S encoding gene. We evaluated the primers and compared their performance with the commonly used ITS1f primer by 454-sequencing of both artificially assembled templates and field samples. When the entire ITS region was amplified using the ITS1f/ITS4 primer combination, we found strong bias against species with longer amplicons. This problem could be overcome by using the new primers, which produce shorter amplicons and better preserve the quantitative composition of the template. In addition, the new primers yielded more diverse amplicon communities than the ITS1f primer.