MAT – gene structure and mating behavior of Hymenoscyphus fraxineus and Hymenoscyphus albidus

The Native Hymenoscyphus albidus and the Invasive Hymenoscyphus fraxineus Are Similar in Their Necrotrophic Growth Phase in Ash Leaves

The populations of European ash and its harmless fungal associate Hymenoscyphus albidus are in decline owing to ash dieback caused by the invasive Hymenoscyphus fraxineus, a fungus that in its native range in Asia is a harmless leaf endophyte of local ash species. To clarify the behavior of H. albidus and its spatial and temporal niche overlap with the invasive relative, we used light microscopy, fungal species-specific qPCR assays, and PacBio long-read amplicon sequencing of the ITS1-5.8S-ITS2 region to examine fungal growth and species composition in attached leaves of European ash. The plant material was collected from a healthy stand in central Norway, where ash saplings in late autumn showed leaflet vein necrosis like that commonly related to H. fraxineus. For reference, leaflet samples were analyzed from stands with epidemic level of ash dieback in southeastern Norway and Estonia. While H. albidus was predominant in the necrotic veins in the healthy stand, H. fraxineus was predominant in the diseased stands. Otherwise, endophytes with pathogenic potential in the genera Venturia (anamorph Fusicladium), Mycosphaerella (anamorph Ramularia), and Phoma, and basidiomycetous yeasts formed the core leaflet mycobiome both in the healthy and diseased stands. In necrotic leaf areas with high levels of either H. albidus or H. fraxineus DNA, one common feature was the high colonization of sclerenchyma and phloem, a region from which the ascomata of both species arise. Our data suggest that H. albidus can induce necrosis in ash leaves, but that owing to low infection pressure, this first takes place in tissues weakened by autumn senescence, 1-2 months later in the season than what is characteristic of H. fraxineus at an epidemic phase of ash dieback. The most striking difference between these fungi would appear to be the high fecundity of H. fraxineus. The adaptation to a host that is phylogenetically closely related to European ash, a tree species with high occurrence frequency in Europe, and the presence of environmental conditions favorable to H. fraxineus life cycle completion in most years may enable the build-up of high infection pressure and challenge of leaf defense prior to autumn senescence.

Population Genomics Reveals Population Structure and Mating-Type Loci in Marssonina brunnea

Marssonina brunnea is an important fungal pathogen of poplar trees. We collected 32 M. brunnea f.sp. multigermtubi (MbMu) and three M. brunnea f.sp. monogermtubi (MbMo) isolates from four poplar species in three Chinese regions and performed genome resequencing. An annotation of SNPs of MbMu indicated that the SNPs potentially have a functional effect on 69.2% of the predicted genes. Using the SNP dataset of nonredundant isolates, a structure and principal component analysis revealed that MbMu and MbMo belong to two genetically distinct populations. By contrast, subpopulation structures could not be found among MbMu isolates. A neighbor-net analysis and a homoplasy index test provided evidence of recombination among MbMu isolates. The short distance (109–174 bp) of linkage disequilibrium half-decay supported the presence of a high level of recombination in the MbMu population. The genetic architectures of the MAT loci of MbMu and MbMo were revealed by searching genome assemblies or by homology-based cloning, and a BLAST search verified each isolate carrying one of the two opposite MAT loci. This study revealed that the MbMu population contains a wide range of functional variants, shows high-frequency recombination, and exhibits heterothallic mating systems, indicating high evolutionary potential and a resultant threat to poplar plantations.

Genomic landscape of a relict fir-associated fungus reveals rapid convergent adaptation towards endophytism

Comparative and pan-genomic analyses of the endophytic fungus Pezicula neosporulosa (Helotiales, Ascomycota) from needles of the relict fir, Abies beshanzuensis, showed expansions of carbohydrate metabolism and secondary metabolite biosynthetic genes characteristic for unrelated plant-beneficial helotialean, such as dark septate endophytes and ericoid mycorrhizal fungi. The current species within the relatively young Pliocene genus Pezicula are predominantly saprotrophic, while P. neosporulosa lacks such features. To understand the genomic background of this putatively convergent evolution, we performed population analyses of 77 P. neosporulosa isolates. This revealed a mosaic structure of a dozen non-recombining and highly genetically polymorphic subpopulations with a unique mating system structure. We found that one idiomorph of a probably duplicated mat1-2 gene was found in putatively heterothallic isolates, while the other co-occurred with mat1-1 locus suggesting homothallic reproduction for these strains. Moreover, 24 and 81 genes implicated in plant cell-wall degradation and secondary metabolite biosynthesis, respectively, showed signatures of the balancing selection. These findings highlight the evolutionary pattern of the two gene families for allowing the fungus a rapid adaptation towards endophytism and facilitating diverse symbiotic interactions.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07837-2.

Detection of a Conspecific Mycovirus in Two Closely Related Native and Introduced Fungal Hosts and Evidence for Interspecific Virus Transmission

Hymenoscyphus albidus is a native fungus in Europe where it behaves as a harmless decomposer of leaves of common ash. Its close relative Hymenoscyphus fraxineus was introduced into Europe from Asia and currently threatens ash (Fraxinus sp.) stands all across the continent causing ash dieback. H. fraxineus isolates from Europe were previously shown to harbor a mycovirus named Hymenoscyphus fraxineus Mitovirus 1 (HfMV1). In the present study, we describe a conspecific mycovirus that we detected in H. albidus. HfMV1 was consistently identified in H. albidus isolates (mean prevalence: 49.3%) which were collected in the sampling areas before the arrival of ash dieback. HfMV1 strains in both fungal hosts contain a single ORF of identical length (717 AA) for which a mean pairwise identity of 94.5% was revealed. The occurrence of a conspecific mitovirus in H. albidus and H. fraxineus is most likely the result of parallel virus evolution in the two fungal hosts. HfMV1 sequences from H. albidus showed a higher nucleotide diversity and a higher number of mutations compared to those from H. fraxineus, probably due to a bottleneck caused by the introduction of H. fraxineus in Europe. Our data also points to multiple interspecific virus transfers from H. albidus to H. fraxineus, which could have contributed to the intraspecific virus diversity found in H. fraxineus.

Propagule Pressure Build-Up by the Invasive Hymenoscyphus fraxineus Following Its Introduction to an Ash Forest Inhabited by the Native Hymenoscyphus albidus

Dieback of European ash, caused by the ascomycete Hymenoscyphus fraxineus originating from Asia, has rapidly spread across Europe, and is threatening this keystone tree at a continental scale. High propagule pressure is characteristic to invasive species. Consistently, the enormous production of windborne ascospores by H. fraxineus in an ash forest with epidemic level of disease obviously facilitates its invasiveness and long distance spread. To understand the rate of build-up of propagule pressure by this pathogen following its local introduction, during 2011-2017 we monitored its sporulation at a newly infested ash stand in south-western Norway characterized with mild winters and cool summers. We also monitored the propagule pressure by Hymenoscyphus albidus, a non-pathogenic native species that competes for the same sporulation niche with H. fraxineus. During the monitoring period, crown condition of ash trees had impaired, and 20% of the dominant trees were severely damaged in 2017. H. fraxineus showed an exponential increase in spore production between 2012 and 2015, followed by drastic decline in 2016 and 2017. During 2011-2013, the two Hymenoscyphus species showed similar sporulation level, but thereafter spores of H. albidus were no longer detected. The data suggest that following local introduction, the population of H. fraxineus reaches rapidly an exponential growth stage if the local weather conditions are favorable for ascomata maturation across years. In the North Atlantic climate, summer temperatures critically influence the pathogen infection pressure, warm summers allowing the population to grow according to its biotic potential, whereas cold summers can cause a drastic decline in propagule pressure.

The ash dieback invasion of Europe was founded by two genetically divergent individuals

Accelerating international trade and climate change make pathogen spread an increasing concern. Hymenoscyphus fraxineus, the causal agent of ash dieback, is a fungal pathogen that has been moving across continents and hosts from Asian to European ash. Most European common ash trees (Fraxinus excelsior) are highly susceptible to H. fraxineus, although a minority (~5%) have partial resistance to dieback. Here, we assemble and annotate a H. fraxineus draft genome which approaches chromosome scale. Pathogen genetic diversity across Europe and in Japan, reveals a strong bottleneck in Europe, though a signal of adaptive diversity remains in key host interaction genes. We find that the European population was founded by two divergent haploid individuals. Divergence between these haplotypes represents the ancestral polymorphism within a large source population. Subsequent introduction from this source would greatly increase adaptive potential of the pathogen. Thus, further introgression of H. fraxineus into Europe represents a potential threat and Europe-wide biological security measures are needed to manage this disease.

The homothallic mating-type locus of the conifer needle endophyte Phialocephala scopiformis DAOMC 229536 (order Helotiales)

We describe the complete mating-type (MAT) locus for Phialocephala scopiformis Canadian Collection of Fungal Cultures (DAOMC) 229536 - a basal lineage within Vibrisseaceae. This strain is of interest due to its ability to produce the important antiinsectan rugulosin. We also provide some of the first insights into the genome structure and gene inventory of nonclavicipitalean endophytes. Sequence was obtained through shotgun sequencing of the entire P. scopiformis genome, and the MAT locus was then determined by comparing this genomic sequence to known MAT loci within the Phialocephala fortinii s.l.-Acephala applanata species complex. We also tested the relative levels of sequence conservation for MAT genes within Vibrisseaceae (n = 10), as well as within the Helotiales (n = 27). Our results: (1) show a homothallic gene arrangement for P. scopiformis [MAT1-1-1, MAT1-2-1, and MAT1-1-3 genes are present], (2) increase the genomic survey of homothallism within Vibrisseaceae, (3) confirm that P. scopiformis contains a unique S-adenosyl-l-methionine-dependent methyltransferase (SAM-Mtase) gene proximal to its MAT locus, while also lacking a cytoskeleton assembly control protein (sla2) gene, and (4) indicate that MAT1-1-1 is the more useful molecular marker amongst the MAT genes for phylogenetic reconstructions aimed at tracking evolutionary shifts in reproductive strategy and/or MAT loci gene composition within the Helotiales.

Organization and evolution of mating-type genes in three Stagonosporopsis species causing gummy stem blight of cucurbits and leaf spot and dry rot of papaya

Population divergence and speciation of closely related lineages can result from reproductive differences leading to genetic isolation. An increasing number of fungal diseases of plants and animals have been determined to be caused by morphologically indistinguishable species that are genetically distinct, thereby representing cryptic species. We were interested in identifying if mating systems among three Stagonosporopsis species (S. citrulli, S. cucurbitacearum, and S. caricae) causing gummy stem blight (GSB) of cucurbits or leaf spot and dry rot of papaya differed, possibly underlying species divergence. Additionally, we were interested in identifying evolutionary pressures acting on the genes controlling mating in these fungi. The mating-type loci (MAT1) of three isolates from each of the three species were identified in draft genome sequences. For the three species, MAT1 was structurally identical and contained both mating-type genes necessary for sexual reproduction, which suggests that all three species are homothallic. However, both MAT1-1-1 and MAT1-2-1 were divergent among species showing rapid evolution with a much greater number of amino acid-changing substitutions detected for the reproductive genes compared with genes flanking MAT1. Positive selection was detected in MAT1-2-1, especially in the highly conserved high mobility group (MATA_HMG-box) domain. Thus, the mating-type genes are rapidly evolving in GSB fungi, but a difference in mating systems among the three species does not underlie their divergence.

A role for the asexual spores in infection of Fraxinus excelsior by the ash-dieback fungus Hymenoscyphus fraxineus

The invasive pathogen, ash dieback fungus Hymenoscyphus fraxineus, is spreading rapidly across Europe. It shows high levels of outcrossing and limited population structure, even at the epidemic front. The anamorphic (asexual) form produces prolific conidia, thought to function solely as spermatia (male gametes), facilitating gene flow between sympatric strains. Here, we show that conidia are capable of germination on ash leaves and in vitro, and can infect seedlings via leaves or soil. In leaves, germlings form structures resembling fruiting bodies. Additionally, H. fraxineus colonises ash debris and grows in soil in the absence of ash tissues. We propose an amended life-cycle in which wind-dispersed, insect-vectored or water-spread conidia infect ash and may sporulate in planta, as well as in forest debris. This amplifies inoculum levels of different strains in ash stands. In combination with their function as spermatia, conidia thus act to maximise gene flow between sympatric strains, including those originally present at low inoculum. Such mixing increases evolutionary potential, as well as enhancing the likelihood of gene introgression from closely-related strains or assimilation of further genetic diversity from parental Asian populations. This scenario increases the adaptability of H. fraxineus to new climates and, indeed, onto new host species.

Mating type genes in the genus Neofusicoccum: Mating strategies and usefulness in species delimitation

The genus Neofusicoccum includes species with wide geographical and plant host distribution, some of them of economic importance. The genus currently comprises 27 species that are difficult to identify based on morphological features alone. Thus, species differentiation is based on phylogenetic species recognition using multigene genealogies. In this study, we characterised the mating type genes of Neofusicoccum species. Specific primers were designed to amplify and sequence MAT genes in several species and a PCR-based mating type diagnostic assay was developed. Homothallism was the predominant mating strategy among the species tested. Furthermore, the potential of mating type gene sequences for species delimitation was evaluated. Phylogenetic analyses were performed on both MAT genes and compared with multigene genealogies using sequences of the ribosomal internal transcribed spacer region, translation elongation factor 1-alpha and beta-tubulin. Phylogenies based on mating type genes could discriminate between the species analysed and are in concordance with the results obtained with the more conventional multilocus phylogenetic analysis approach. Thus, MAT genes represent a powerful tool to delimit cryptic species in the genus Neofusicoccum.