Evolution, selection and isolation: a genomic view of speciation in fungal plant pathogens

Multiple Horizontal Mini-chromosome Transfers Drive Genome Evolution of Clonal Blast Fungus Lineages

Crop disease pandemics are often driven by asexually reproducing clonal lineages of plant pathogens that reproduce asexually. How these clonal pathogens continuously adapt to their hosts despite harboring limited genetic variation, and in absence of sexual recombination remains elusive. Here, we reveal multiple instances of horizontal chromosome transfer within pandemic clonal lineages of the blast fungus Magnaporthe (Syn. Pyricularia) oryzae. We identified a horizontally transferred 1.2Mb accessory mini-chromosome which is remarkably conserved between M. oryzae isolates from both the rice blast fungus lineage and the lineage infecting Indian goosegrass (Eleusine indica), a wild grass that often grows in the proximity of cultivated cereal crops. Furthermore, we show that this mini-chromosome was horizontally acquired by clonal rice blast isolates through at least nine distinct transfer events over the past three centuries. These findings establish horizontal mini-chromosome transfer as a mechanism facilitating genetic exchange among different host-associated blast fungus lineages. We propose that blast fungus populations infecting wild grasses act as genetic reservoirs that drive genome evolution of pandemic clonal lineages that afflict cereal crops.

The Sordariomycetes: an expanding resource with Big Data for mining in evolutionary genomics and transcriptomics

Advances in genomics and transcriptomics accompanying the rapid accumulation of omics data have provided new tools that have transformed and expanded the traditional concepts of model fungi. Evolutionary genomics and transcriptomics have flourished with the use of classical and newer fungal models that facilitate the study of diverse topics encompassing fungal biology and development. Technological advances have also created the opportunity to obtain and mine large datasets. One such continuously growing dataset is that of the Sordariomycetes, which exhibit a richness of species, ecological diversity, economic importance, and a profound research history on amenable models. Currently, 3,574 species of this class have been sequenced, comprising nearly one-third of the available ascomycete genomes. Among these genomes, multiple representatives of the model genera Fusarium, Neurospora, and Trichoderma are present. In this review, we examine recently published studies and data on the Sordariomycetes that have contributed novel insights to the field of fungal evolution via integrative analyses of the genetic, pathogenic, and other biological characteristics of the fungi. Some of these studies applied ancestral state analysis of gene expression among divergent lineages to infer regulatory network models, identify key genetic elements in fungal sexual development, and investigate the regulation of conidial germination and secondary metabolism. Such multispecies investigations address challenges in the study of fungal evolutionary genomics derived from studies that are often based on limited model genomes and that primarily focus on the aspects of biology driven by knowledge drawn from a few model species. Rapidly accumulating information and expanding capabilities for systems biological analysis of Big Data are setting the stage for the expansion of the concept of model systems from unitary taxonomic species/genera to inclusive clusters of well-studied models that can facilitate both the in-depth study of specific lineages and also investigation of trait diversity across lineages. The Sordariomycetes class, in particular, offers abundant omics data and a large and active global research community. As such, the Sordariomycetes can form a core omics clade, providing a blueprint for the expansion of our knowledge of evolution at the genomic scale in the exciting era of Big Data and artificial intelligence, and serving as a reference for the future analysis of different taxonomic levels within the fungal kingdom.

Multi-parental fungal mapping population study to detect genomic regions associated with Pyrenophora teres f. teres virulence

In recent years multi-parental mapping populations (MPPs) have been widely adopted in many crops to detect quantitative trait loci (QTLs) as this method can compensate for the limitations of QTL analyses using bi-parental mapping populations. Here we report the first multi-parental nested association mapping (MP-NAM) population study used to detect genomic regions associated with host-pathogenic interactions. MP-NAM QTL analyses were conducted on 399 Pyrenophora teres f. teres individuals using biallelic, cross-specific and parental QTL effect models. A bi-parental QTL mapping study was also conducted to compare the power of QTL detection between bi-parental and MP-NAM populations. Using MP-NAM with 399 individuals detected a maximum of eight QTLs with a single QTL effect model whilst only a maximum of five QTLs were detected with an individual bi-parental mapping population of 100 individuals. When reducing the number of isolates in the MP-NAM to 200 individuals the number of QTLs detected remained the same for the MP-NAM population. This study confirms that MPPs such as MP-NAM populations can be successfully used in detecting QTLs in haploid fungal pathogens and that the power of QTL detection with MPPs is greater than with bi-parental mapping populations.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Enforcement of Postzygotic Species Boundaries in the Fungal Kingdom

Understanding the molecular basis of speciation is a primary goal in evolutionary biology. The formation of the postzygotic reproductive isolation that causes hybrid dysfunction, thereby reducing gene flow between diverging populations, is crucial for speciation. Using various advanced approaches, including chromosome replacement, hybrid introgression and transcriptomics, population genomics, and experimental evolution, scientists have revealed multiple mechanisms involved in postzygotic barriers in the fungal kingdom. These results illuminate both unique and general features of fungal speciation. Our review summarizes experiments on fungi exploring how Dobzhansky-Muller incompatibility, killer meiotic drive, chromosome rearrangements, and antirecombination contribute to postzygotic reproductive isolation. We also discuss possible evolutionary forces underlying different reproductive isolation mechanisms and the potential roles of the evolutionary arms race under the Red Queen hypothesis and epigenetic divergence in speciation.

Telomere-to-telomere genome sequence of the model mould pathogen Aspergillus fumigatus

The pathogenic fungus Aspergillus fumigatus is a major etiological agent of fungal invasive and chronic diseases affecting tens of millions of individuals worldwide. Draft genome sequences of two clinical isolates (Af293 and A1163) are commonly used as reference genomes for analyses of clinical and environmental strains. However, the reference sequences lack coverage of centromeres, an accurate sequence for ribosomal repeats, and a comprehensive annotation of chromosomal rearrangements such as translocations and inversions. Here, we used PacBio Single Molecule Real-Time (SMRT), Oxford Nanopore and Illumina HiSeq sequencing for de novo genome assembly and polishing of two laboratory reference strains of A. fumigatus, CEA10 (parental isolate of A1163) and its descendant A1160. We generated full length chromosome assemblies and a comprehensive telomere-to-telomere coverage for CEA10 and near complete assembly of A1160 including ribosomal repeats and the sequences of centromeres, which we discovered to be composed of long transposon elements. We envision these high-quality reference genomes will become fundamental resources to study A. fumigatus biology, pathogenicity and virulence, and to discover more effective treatments against diseases caused by this fungus.

The fungus Aspergillus fumigatus causes invasive and chronic diseases worldwide. Here, Bowyer et al. use long-read and short-read sequencing to generate complete chromosome assemblies and telomere-to-telomere coverage for two isolates, thus providing high-quality reference genomes as fundamental resources to study this pathogen.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities

Many species of fungal plant pathogens coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. The rice blast fungus Pyricularia oryzae is a widespread model plant pathogen displaying population subdivision. However, most studies of natural variation in P. oryzae have been limited in genomic or geographic resolution, and host adaptation is the only factor that has been investigated extensively as a contributor to population subdivision. In an effort to complement previous studies, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus covering a broad geographical range. Using single-nucleotide polymorphism genotyping data for 886 isolates sampled from 152 sites in 51 countries, we showed that population subdivision of P. oryzae in one recombining and three clonal lineages with broad distributions persisted with deeper sampling. We also extended previous findings by showing further population subdivision of the recombining lineage into one international and three Asian clusters, and by providing evidence that the three clonal lineages of P. oryzae were found in areas with different prevailing environmental conditions, indicating niche separation. Pathogenicity tests and bioinformatic analyses using an extended set of isolates and rice varieties indicated that partial specialization to rice subgroups contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong additional barriers to gene flow between clonal lineages. Our results demonstrate that the spread of a fungal pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct, widely distributed, lineages.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

Predicting global numbers of teleomorphic ascomycetes

Sexual reproduction is the basic way to form high genetic diversity and it is beneficial in evolution and speciation of fungi. The global diversity of teleomorphic species in Ascomycota has not been estimated. This paper estimates the species number for sexual ascomycetes based on five different estimation approaches, viz. by numbers of described fungi, by fungus:substrate ratio, by ecological distribution, by meta-DNA barcoding or culture-independent studies and by previous estimates of species in Ascomycota. The assumptions were made with the currently most accepted, “2.2–3.8 million” species estimate and results of previous studies concluding that 90% of the described ascomycetes reproduce sexually. The Catalogue of Life, Species Fungorum and published research were used for data procurement. The average value of teleomorphic species in Ascomycota from all methods is 1.86 million, ranging from 1.37 to 2.56 million. However, only around 83,000 teleomorphic species have been described in Ascomycota and deposited in data repositories. The ratio between described teleomorphic ascomycetes to predicted teleomorphic ascomycetes is 1:22. Therefore, where are the undiscovered teleomorphic ascomycetes? The undescribed species are no doubt to be found in biodiversity hot spots, poorly-studied areas and species complexes. Other poorly studied niches include extremophiles, lichenicolous fungi, human pathogens, marine fungi, and fungicolous fungi. Undescribed species are present in unexamined collections in specimen repositories or incompletely described earlier species. Nomenclatural issues, such as the use of separate names for teleomorph and anamorphs, synonyms, conspecific names, illegitimate and invalid names also affect the number of described species. Interspecies introgression results in new species, while species numbers are reduced by extinctions.

The First Mitochondrial Genome of Ciborinia camelliae and Its Position in the Sclerotiniaceae Family

Ciborinia camelliae is the causal agent of camellia flower blight (CFB). It is a hemibiotrophic pathogen, inoperculate Discomycete of the family Sclerotiniaceae. It shows host and organ specificity infecting only flowers of species belonging to the genus Camellia, causing serious damage to the ornamental component of the plant. In this work, the first mitochondrial genome of Ciborinia camellia is reported. The mitogenome was obtained by combining Illumina short read and Nanopore long read technology. To resolve repetitive elements, specific primers were designed and used for Sanger sequencing. The manually curated mitochondrial DNA (mtDNA) of the Italian strain DSM 112729 is a circular sequence of 114,660 bp, with 29.6% of GC content. It contains two ribosomal RNA genes, 33 transfer RNAs, one RNase P gene, and 62 protein-coding genes. The latter include one gene coding for a ribosomal protein (rps3) and the 14 typical proteins involved in the oxidative metabolism. Moreover, a partial mtDNA assembled from a contig list was obtained from the deposited genome assembly of a New Zealand strain of C. camelliae. The present study contributes to understanding the mitogenome arrangement and the evolution of this phytopathogenic fungus in comparison to other Sclerotiniaceae species and confirms the usefulness of mitochondrial analysis to define phylogenetic positioning of this newly sequenced species.

"Jumping Jack": Genomic Microsatellites Underscore the Distinctiveness of Closely Related Pseudoperonospora cubensis and Pseudoperonospora humuli and Provide New Insights Into Their Evolutionary Past

Downy mildews caused by obligate biotrophic oomycetes result in severe crop losses worldwide. Among these pathogens, Pseudoperonospora cubensis and P. humuli, two closely related oomycetes, adversely affect cucurbits and hop, respectively. Discordant hypotheses concerning their taxonomic relationships have been proposed based on host-pathogen interactions and specificity evidence and gene sequences of a few individuals, but population genetics evidence supporting these scenarios is missing. Furthermore, nuclear and mitochondrial regions of both pathogens have been analyzed using microsatellites and phylogenetically informative molecular markers, but extensive comparative population genetics research has not been done. Here, we genotyped 138 current and historical herbarium specimens of those two taxa using microsatellites (SSRs). Our goals were to assess genetic diversity and spatial distribution, to infer the evolutionary history of P. cubensis and P. humuli, and to visualize genome-scale organizational relationship between both pathogens. High genetic diversity, modest gene flow, and presence of population structure, particularly in P. cubensis, were observed. When tested for cross-amplification, 20 out of 27 P. cubensis-derived gSSRs cross-amplified DNA of P. humuli individuals, but few amplified DNA of downy mildew pathogens from related genera. Collectively, our analyses provided a definite argument for the hypothesis that both pathogens are distinct species, and suggested further speciation in the P. cubensis complex.

Genomic rearrangements generate hypervariable mini-chromosomes in host-specific isolates of the blast fungus

Supernumerary mini-chromosomes–a unique type of genomic structural variation–have been implicated in the emergence of virulence traits in plant pathogenic fungi. However, the mechanisms that facilitate the emergence and maintenance of mini-chromosomes across fungi remain poorly understood. In the blast fungus Magnaporthe oryzae (Syn. Pyricularia oryzae), mini-chromosomes have been first described in the early 1990s but, until very recently, have been overlooked in genomic studies. Here we investigated structural variation in four isolates of the blast fungus M. oryzae from different grass hosts and analyzed the sequences of mini-chromosomes in the rice, foxtail millet and goosegrass isolates. The mini-chromosomes of these isolates turned out to be highly diverse with distinct sequence composition. They are enriched in repetitive elements and have lower gene density than core-chromosomes. We identified several virulence-related genes in the mini-chromosome of the rice isolate, including the virulence-related polyketide synthase Ace1 and two variants of the effector gene AVR-Pik. Macrosynteny analyses around these loci revealed structural rearrangements, including inter-chromosomal translocations between core- and mini-chromosomes. Our findings provide evidence that mini-chromosomes emerge from structural rearrangements and segmental duplication of core-chromosomes and might contribute to adaptive evolution of the blast fungus.

The genomes of plant pathogens often exhibit an architecture that facilitates high rates of dynamic rearrangements and genetic diversification in virulence associated regions. These regions, which tend to be gene sparse and repeat rich, are thought to serve as a cradle for adaptive evolution. Supernumerary chromosomes, i.e. chromosomes that are only present in some but not all individuals of a species, are a special type of structural variation that have been observed in plants, animals, and fungi. Here we identified and studied supernumerary mini-chromosomes in the blast fungus Magnaporthe oryzae, a pathogen that causes some of the most destructive plant diseases. We found that rice, foxtail millet and goosegrass isolates of this pathogen contain mini-chromosomes with distinct sequence composition. All mini-chromosomes are rich in repetitive genetic elements and have lower gene densities than core-chromosomes. Further, we identified virulence-related genes on the mini-chromosome of the rice isolate. We observed large-scale genomic rearrangements around these loci, indicative of a role of mini-chromosomes in facilitating genome dynamics. Taken together, our results indicate that mini-chromosomes contribute to genome rearrangements and possibly adaptive evolution of the blast fungus.

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

Meiotic recombination in the offspring of Microbotryum hybrids and its impact on pathogenicity

BACKGROUND

Hybridization is a central mechanism in evolution, producing new species or introducing important genetic variation into existing species. In plant-pathogenic fungi, adaptation and specialization to exploit a host species are key determinants of evolutionary success. Here, we performed experimental crosses between the two pathogenic Microbotryum species, M. lychnidis-dioicae and M. silenes-acaulis that are specialized to different hosts. The resulting offspring were analyzed on phenotypic and genomic levels to describe genomic characteristics of hybrid offspring and genetic factors likely involved in host-specialization.

RESULTS

Genomic analyses of interspecific fungal hybrids revealed that individuals were most viable if the majority of loci were inherited from one species. Interestingly, species-specific loci were strictly controlled by the species' origin of the mating type locus. Moreover we detected signs of crossing over and chromosome duplications in the genomes of the analyzed hybrids. In Microbotryum, mitochondrial DNA was found to be uniparentally inherited from the a2 mating type. Genome comparison revealed that most gene families are shared and the majority of genes are conserved between the two species, indicating very similar biological features, including infection and pathogenicity processes. Moreover, we detected 211 candidate genes that were retained under host-driven selection of backcrossed lines. These genes and might therefore either play a crucial role in host specialization or be linked to genes that are essential for specialization.

CONCLUSION

The combination of genome analyses with experimental selection and hybridization is a promising way to investigate host-pathogen interactions. This study manifests genetic factors of host specialization that are required for successful biotrophic infection of the post-zygotic stage, but also demonstrates the strong influence of intra-genomic conflicts or instabilities on the viability of hybrids in the haploid host-independent stage.

Population Analyses Reveal Two Host-Adapted Clades of Pseudoperonospora cubensis, the Causal Agent of Cucurbit Downy Mildew, on Commercial and Wild Cucurbits

Pseudoperonospora cubensis, the causal agent of cucurbit downy mildew, is an airborne, obligate oomycete pathogen that re-emerged in 2004 and causes foliar disease and yield losses in all major cucurbit crops in the United States. Approximately 60 species in the family Cucurbitaceae have been reported as hosts of P. cubensis. Commercial hosts including cucumber, cantaloupe, pumpkin, squash, and watermelon are grown in North Carolina and many host species occur in the wild as weeds. Little is known about the contribution of wild cucurbits to the yearly epidemic; thus, this study aimed to determine the role of commercial and wild cucurbits in the structuring of P. cubensis populations in North Carolina, a region with high pathogen diversity. Ten microsatellite markers were used to analyze 385 isolates from six commercial and four wild cucurbits from three locations representing different growing regions across North Carolina. Population analyses revealed that wild and commercial cucurbits are hosts of P. cubensis in the United States, that host is the main factor structuring P. cubensis populations, and that P. cubensis has two distinct, host-adapted clades at the cucurbit species level, with clade 1 showing random mating and evidence of recombination and clade 2 showing nonrandom mating and no evidence of recombination. Our findings have implications for disease management because clade-specific factors such as host susceptibility and inoculum availability of each clade by region may influence P. cubensis outbreaks in different commercial cucurbits, timing of fungicide applications, and phenotyping for breeding efforts.

Genome comparisons suggest an association between Ceratocystis host adaptations and effector clusters in unique transposable element families

Ceratocystis fimbriata is a host specific fungal pathogen of sweet potato (Ipomoea batatas). The closely related species, C. manginecans, is an important pathogen of trees (e.g. Acacia mangium and Mangifera indica) but has never been isolated from tuber crops. The genetic factors that determine the host range and host specificity of these species have not been determined. The aim of this study was to compare the genomes of C. fimbriata and C. manginecans in order to identify species-specific genetic differences that could be associated with host specificity. This included whole-genome alignments as well as comparisons of gene content and transposable elements (TEs). The genomes of the two species were found to be very similar, sharing similar catalogues of CAZymes, peptidases and lipases. However, the genomes of the two species also varied, harbouring species-specific genes (e.g. small secreted effectors, nutrient processing proteins and stress response proteins). A portion of the TEs identified (17%) had a unique distribution in each species. Transposable elements appeared to have played a prominent role in the divergence of the two species because they were strongly associated with chromosomal translocations and inversions as well as with unique genomic regions containing species-specific genes. Two large effector clusters, with unique TEs in each species, were identified. These effectors displayed non-synonymous mutations and deletions, conserved within a species, and could serve as mutational hot-spots for the development of host specificity in the two species.

Threats to global food security from emerging fungal and oomycete crop pathogens

Emerging fungal and oomycete pathogens infect staple calorie crops and economically important commodity crops, thereby posing a significant risk to global food security. Our current agricultural systems - with emphasis on intensive monoculture practices - and globalized markets drive the emergence and spread of new pathogens and problematic traits, such as fungicide resistance. Climate change further promotes the emergence of pathogens on new crops and in new places. Here we review the factors affecting the introduction and spread of pathogens and current disease control strategies, illustrating these with the historic example of the Irish potato famine and contemporary examples of soybean rust, wheat blast and blotch, banana wilt and cassava root rot. Our Review looks to the future, summarizing what we see as the main challenges and knowledge gaps, and highlighting the direction that research must take to face the challenge of emerging crop pathogens.

Extensive chromosomal rearrangements and rapid evolution of novel effector superfamilies contribute to host adaptation and speciation in the basal ascomycetous fungi

The basal ascomycetes in genus Taphrina have strict host specificity and coevolution with their host plants, making them appealing models for studying the genomic basis of ecological divergence and host adaption. We therefore performed genome sequencing and comparative genomics of different Taphrina species with distinct host ranges to reveal their evolution. We identified frequent chromosomal rearrangements and highly dynamic lineage-specific (LS) genomic regions in Taphrina genomes. The LS regions occur at the flanking regions of chromosomal breakpoints, and are greatly enriched for DNA repeats, non-core genes, and in planta up-regulated genes. Furthermore, we identified hundreds of candidate secreted effector proteins (CSEPs) that are commonly organized in gene clusters that form distinct AT-rich isochore-like regions. Nearly half of the CSEPs constitute two novel superfamilies with modular structures unique to Taphrina. These CSEPs are commonly up-regulated during infection, enriched in the LS regions, evolved faster, and underwent extensive gene gain and loss in different species. In addition to displaying signatures of positive selection, functional characterization of selected CSEP genes confirmed their roles in suppression of plant defence responses. Overall, our results showed that extensive chromosomal rearrangements and rapidly evolving CSEP superfamilies play important roles in speciation and host adaptation in the early-branching ascomycetous fungi.

Population Genomic Analyses Reveal Connectivity via Human-Mediated Transport across Populus Plantations in North America and an Undescribed Subpopulation of Sphaerulina musiva

Domestication of plant species has affected the evolutionary dynamics of plant pathogens in agriculture and forestry. A model system for studying the consequences of plant domestication on the evolution of an emergent plant disease is the fungal pathogen Sphaerulina musiva. This ascomycete causes leaf spot and stem canker disease of Populus spp. and their hybrids. A population genomics approach was used to determine the degree of population structure and evidence for selection on the North American population of S. musiva. In total, 122 samples of the fungus were genotyped identifying 120,016 single-nucleotide polymorphisms after quality filtering. In North America, S. musiva has low to moderate degrees of differentiation among locations. Three main genetic clusters were detected: southeastern United States, midwestern United States and Canada, and a new British Columbia cluster (BC2). Population genomics suggest that BC2 is a novel genetic cluster from central British Columbia, clearly differentiated from previously reported S. musiva from coastal British Columbia, and the product of a single migration event. Phenotypic measurements from greenhouse experiments indicate lower aggressiveness of BC2 on Populus trichocarpa. In summary, S. musiva has geographic structure across broad regions indicative of gene flow among clusters. The interconnectedness of the North American S. musiva populations across large geographic distances further supports the hypothesis of anthropogenic-facilitated transport of the pathogen.

New insights into the evolution of host specificity of three Penicillium species and the pathogenicity of P. Italicum involving the infection of Valencia orange (Citrus sinensis)

Blue and green molds, the common phenotypes of post-harvest diseases in fruits, are mainly caused by Penicillium fungal species, including P. italicum, P. digitatum, and P. expansum. We sequenced and assembled the genome of a P. italicum strain, which contains 31,034,623 bp with 361 scaffolds and 627 contigs. The mechanisms underlying the evolution of host specificity among the analyzed Penicillium species were associated with the expansion of protein families, genome restructuring, horizontal gene transfer, and positive selection pressure. A dual-transcriptome analysis following the infection of Valencia orange (Citrus sinensis) by P. italicum resulted in the annotation of 9,307 P. italicum genes and 24,591 Valencia orange genes. The pathogenicity of P. italicum may be due to the activation of effectors, including 51 small secreted cysteine-rich proteins, 110 carbohydrate-active enzymes, and 12 G protein-coupled receptors. Additionally, 211 metabolites related to the interactions between P. italicum and Valencia orange were identified by gas chromatography-time of flight mass spectrography, three of which were further confirmed by ultra-high performance liquid chromatography triple quadrupole mass spectrometry. A metabolomics analysis indicated that P. italicum pathogenicity is associated with the sphingolipid and salicylic acid signaling pathways. Moreover, a correlation analysis between the metabolite contents and gene expression levels suggested that P. italicum induces carbohydrate metabolism in Valencia orange fruits as part of its infection strategy. This study provides useful information regarding the genomic determinants that drive the evolution of host specificity in Penicillium species and clarifies the host-plant specificity during the infection of Valencia orange by P. italicum.

IMPORTANCE

P. italicum GL_Gan1, a local strain in Guangzhou, China, was sequenced. Comparison of the genome of P. italicum GL_Gan1 with other pathogenic Penicillium species, P. digitatum and P. expansum, revealed that the expansion of protein families, genome restructuring, HGT, and positive selection pressure were related to the host range expansion of the analyzed Penicillium species. Moreover, gene gains or losses might be associated with the speciation of these Penicillium species. In addition, the molecular basis of host-plant specificity during the infection of Valencia orange (Citrus sinensis) by P. italicum was also elucidated by transcriptomic and metabolomics analysis. The data presented herein may be useful for further elucidating the molecular basis of the evolution of host specificity of Penicillium species and for illustrating the host-plant specificity during the infection of Valencia orange by P. italicum.

Genome Sequencing and Transcriptome Analysis of the Hop Downy Mildew Pathogen Pseudoperonospora humuli Reveal Species-Specific Genes for Molecular Detection

Pseudoperonospora humuli is an obligate oomycete pathogen of hop (Humulus lupulus) that causes downy mildew, an important disease in most production regions in the Northern Hemisphere. The pathogen can cause a systemic infection in hop, overwinter in the root system, and infect propagation material. Substantial yield loss may occur owing to P. humuli infection of strobiles (seed cones), shoots, and cone-bearing branches. Fungicide application and cultural practices are the primary methods to manage hop downy mildew. However, effective, sustainable, and cost-effective management of downy mildew can be improved by developing early detection systems to inform on disease risk and timely fungicide application. However, no species-specific diagnostic assays or genomic resources are available for P. humuli. The genome of the P. humuli OR502AA isolate was partially sequenced using Illumina technology and assembled with ABySS. The assembly had a minimum scaffold length of 500 bp and an N50 (median scaffold length of the assembled genome) of 19.2 kbp. A total number of 18,656 genes were identified using MAKER standard gene predictions. Additionally, transcriptome assemblies were generated using RNA-seq and Trinity for seven additional P. humuli isolates. Bioinformatics analyses of next generation sequencing reads of P. humuli and P. cubensis (a closely related sister species) identified 242 candidate species-specific P. humuli genes that could be used as diagnostic molecular markers. These candidate genes were validated using polymerase chain reaction against a diverse collection of isolates from P. humuli, P. cubensis, and other oomycetes. Overall, four diagnostic markers were found to be uniquely present in P. humuli. These candidate markers identified through comparative genomics can be used for pathogen diagnostics in propagation material, such as rhizomes and vegetative cuttings, or adapted for biosurveillance of airborne sporangia, an important source of inoculum in hop downy mildew epidemics.

Transposable Elements Adaptive Role in Genome Plasticity, Pathogenicity and Evolution in Fungal Phytopathogens

Transposable elements (TEs) are agents of genetic variability in phytopathogens as they are a source of adaptive evolution through genome diversification. Although many studies have uncovered information on TEs, the exact mechanism behind TE-induced changes within the genome remains poorly understood. Furthermore, convergent trends towards bigger genomes, emergence of novel genes and gain or loss of genes implicate a TE-regulated genome plasticity of fungal phytopathogens. TEs are able to alter gene expression by revamping the cis-regulatory elements or recruiting epigenetic control. Recent findings show that TEs recruit epigenetic control on the expression of effector genes as part of the coordinated infection strategy. In addition to genome plasticity and diversity, fungal pathogenicity is an area of economic concern. A survey of TE distribution suggests that their proximity to pathogenicity genes TEs may act as sites for emergence of novel pathogenicity factors via nucleotide changes and expansion or reduction of the gene family. Through a systematic survey of literature, we were able to conclude that the role of TEs in fungi is wide: ranging from genome plasticity, pathogenicity to adaptive behavior in evolution. This review also identifies the gaps in knowledge that requires further elucidation for a better understanding of TEs' contribution to genome architecture and versatility.

Global distribution of mating types shows limited opportunities for mating across populations of fungi causing boxwood blight disease

Boxwood blight is a disease threat to natural and managed landscapes worldwide. To determine mating potential of the fungi responsible for the disease, Calonectria pseudonaviculata and C. henricotiae, we characterized their mating-type (MAT) loci. Genomes of C. henricotiae, C. pseudonaviculata and two other Calonectria species (C. leucothoes, C. naviculata) were sequenced and used to design PCR tests for mating-type from 268 isolates collected from four continents. All four Calonectria species have a MAT locus that is structurally consistent with the organization found in heterothallic ascomycetes, with just one idiomorph per individual isolate. Mating type was subdivided by species: all C. henricotiae isolates possessed the MAT1-1 idiomorph, whereas all C. pseudonaviculata isolates possessed the MAT1-2 idiomorph. To determine the potential for divergence at the MAT1 locus to present a barrier to interspecific hybridization, evolutionary analysis was conducted. Phylogenomic estimates showed that C. henricotiae and C. pseudonaviculata diverged approximately 2.1 Mya. However, syntenic comparisons, phylogenetic analyses, and estimates of nucleotide divergence across the MAT1 locus and proximal genes identified minimal divergence in this region of the genome. These results show that in North America and parts of Europe, where only C. pseudonaviculata resides, mating is constrained by the absence of MAT1-1. In regions of Europe where C. henricotiae and C. pseudonaviculata currently share the same host and geographic range, it remains to be determined whether or not these two recently diverged species are able to overcome species barriers to mate.

QTL mapping of mycelial growth and aggressiveness to distinct hosts in Ceratocystis pathogens

Some species of Ceratocystis display strong host specificity, such as C. fimbriata sensu stricto that is restricted to sweet potato (Ipomoea batatas) as host. In contrast, the closely related C. manginecans, infects Acacia mangium and Mangifera indica but is not pathogenic to I. batatas. Despite the economic importance of these fungi, knowledge regarding the genetic factors that influence their pathogenicity and host specificity is limited. A recent inheritance study, based on an interspecific cross between C. fimbriata and C. manginecans and the resultant 70 F₁ progeny, confirmed that traits such as mycelial growth rate, spore production and aggressiveness on A. mangium and I. batatas are regulated by multiple genes. In the present study, a quantitative trait locus (QTL) analysis was performed to determine the genomic loci associated with these traits. All 70 progeny isolates were genotyped with SNP markers and a linkage map was constructed. The map contained 467 SNPs, distributed across nine linkage groups, with a total length of 1203 cm. Using the progeny genotypes and phenotypes, one QTL was identified on the linkage map for mycelial growth rate, one for aggressiveness to A. mangium and two for aggressiveness to I. batatas (P < 0.05). Two candidate genes, likely associated with mycelial growth rate, were identified in the QTL region. The three QTLs associated with aggressiveness to different hosts contained candidate genes involved in protein processing, detoxification and regions with effector genes and high transposable element density. The results provide a foundation for studies considering the function of genes regulating various quantitative traits in Ceratocystis.

Comparative genome analysis indicates high evolutionary potential of pathogenicity genes in Colletotrichum tanaceti

Colletotrichum tanaceti is an emerging foliar fungal pathogen of commercially grown pyrethrum (Tanacetum cinerariifolium). Despite being reported consistently from field surveys in Australia, the molecular basis of pathogenicity of C. tanaceti on pyrethrum is unknown. Herein, the genome of C. tanaceti (isolate BRIP57314) was assembled de novo and annotated using transcriptomic evidence. The inferred putative pathogenicity gene suite of C. tanaceti comprised a large array of genes encoding secreted effectors, proteases, CAZymes and secondary metabolites. Comparative analysis of its putative pathogenicity gene profiles with those of closely related species suggested that C. tanaceti likely has additional hosts to pyrethrum. The genome of C. tanaceti had a high repeat content and repetitive elements were located significantly closer to genes inferred to influence pathogenicity than other genes. These repeats are likely to have accelerated mutational and transposition rates in the genome, resulting in a rapid evolution of certain CAZyme families in this species. The C. tanaceti genome showed strong signals of Repeat Induced Point (RIP) mutation which likely caused its bipartite nature consisting of distinct gene-sparse, repeat and A-T rich regions. Pathogenicity genes within these RIP affected regions were likely to have a higher evolutionary rate than the rest of the genome. This "two-speed" genome phenomenon in certain Colletotrichum spp. was hypothesized to have caused the clustering of species based on the pathogenicity genes, to deviate from taxonomic relationships. The large repertoire of pathogenicity factors that potentially evolve rapidly due to the plasticity of the genome, indicated that C. tanaceti has a high evolutionary potential. Therefore, C. tanaceti poses a high-risk to the pyrethrum industry. Knowledge of the evolution and diversity of the putative pathogenicity genes will facilitate future research in disease management of C. tanaceti and other Colletotrichum spp.

Advances in understanding obligate biotrophy in rust fungi

Contents Summary 1190 I. Introduction 1190 II. Rust fungi: a diverse and serious threat to agriculture 1191 III. The different facets of rust life cycles and unresolved questions about their evolution 1191 IV. The biology of rust infection 1192 V. Rusts in the genomics era: the ever-expanding list of candidate effector genes 1195 VI. Functional characterization of rust effectors 1197 VII. Putting rusts to sleep: Pucciniales research outlooks 1201 Acknowledgements 1202 References 1202 SUMMARY: Rust fungi (Pucciniales) are the largest group of plant pathogens and represent one of the most devastating threats to agricultural crops worldwide. Despite the economic importance of these highly specialized pathogens, many aspects of their biology remain obscure, largely because rust fungi are obligate biotrophs. The rise of genomics and advances in high-throughput sequencing technology have presented new options for identifying candidate effector genes involved in pathogenicity mechanisms of rust fungi. Transcriptome analysis and integrated bioinformatics tools have led to the identification of key genetic determinants of host susceptibility to infection by rusts. Thousands of genes encoding secreted proteins highly expressed during host infection have been reported for different rust species, which represents significant potential towards understanding rust effector function. Recent high-throughput in planta expression screen approaches (effectoromics) have pushed the field ahead even further towards predicting high-priority effectors and identifying avirulence genes. These new insights into rust effector biology promise to inform future research and spur the development of effective and sustainable strategies for managing rust diseases.

Asexual reproduction and growth rate: independent and plastic life history traits in Neurospora crassa

Trade-offs among traits influencing fitness are predicted by life history theory because resources allocated to one function are unavailable to another. Here we examine the relationship between two such traits, asexual reproduction and growth rate, in the filamentous fungus Neurospora crassa, where shared genetic and physiological factors and a source-sink energetic relationship between growth and reproduction may constrain the evolution of these traits. To test growth-reproduction relationships in this species, we independently selected on mycelial growth rate or asexual spore production in a heterogeneous lab-derived population and evaluated the response of the non-selected traits. Combined with phenotypes for the 20 wild strains used to produce the heterogeneous population and the genome-wide genotypes of 468 strains, these data show that growth and reproduction are highly plastic in N. crassa and do not trade off either among wild strains or after laboratory selection in two environments. Rather, we find no predictable growth-reproduction relationship in the environments tested, indicating an effective absence of genetic constraint between these traits. Our results suggest that growth rate and asexual reproduction may not respond predictably to environmental change and suggest that reliance on a single trait as a proxy for fitness in fungal studies may be inadvisable.

Mating-type loci of Ustilago esculenta are essential for mating and development

Ustilago esculenta is closely related to the smut fungus Ustilago maydis and, in an endophytic-like life in the plant Zizania latifolia, only infects host stems and causes swollen stems to form edible galls called Jiaobai in China. In order to study its different modes of invasion and sites of symptom development from other smut fungi at the molecular level, we first characterized the a and b mating-type loci of U. esculenta. The a loci contained three a mating-type alleles, encoding two pheromones and one pheromone receptor per allele. The pheromone/receptor system controlled the conjugation formation, the initial step of mating, in which each pheromone was specific for recognition by only one mating partner. In addition, there are at least three b alleles identified in U. esculenta, encoding two subunits of heterodimeric homeodomain transcription factors bE and bW, responsible for hyphal growth and invasiveness. Hyphal formation, elongation and invasion after mating of two compatible partners occurred, only when a heterodimer complex was formed by the bE and bW proteins derived from different alleles. We also demonstrated that even with only one paired pheromone-pheromone receptor, the active b locus heterodimer triggered hyphal growth and infection.

Gene exchange between two divergent species of the fungal human pathogen, Coccidioides

The fungal genus Coccidioides is composed of two species, Coccidioides immitis and Coccidioides posadasii. These two species are the causal agents of coccidioidomycosis, a pulmonary disease also known as valley fever. The two species are thought to have shared genetic material due to gene exchange in spite of their long divergence. To quantify the magnitude of shared ancestry between them, we analyzed the genomes of a population sample from each species. Next, we inferred what is the expected size of shared haplotypes that might be inherited from the last common ancestor of the two species and find a cutoff to find what haplotypes have conclusively been exchanged between species. Finally, we precisely identified the breakpoints of the haplotypes that have crossed the species boundary and measure the allele frequency of each introgression in this sample. We find that introgressions are not uniformly distributed across the genome. Most, but not all, of the introgressions segregate at low frequency. Our results show that divergent species can share alleles, that species boundaries can be porous, and highlight the need for a systematic exploration of gene exchange in fungal species.

Comparative Genomics of Aspergillus flavus S and L Morphotypes Yield Insights into Niche Adaptation

Aspergillus flavus, the primary causal agent for aflatoxin contamination on crops, consists of isolates with two distinct morphologies: isolates of the S morphotype produce numerous small sclerotia and lower numbers of conidia while isolates of the L morphotype produce fewer large sclerotia and abundant conidia. The morphotypes also differ in aflatoxin production with S isolates consistently producing high concentrations of aflatoxin, whereas L isolates range from atoxigenic to highly toxigenic. The production of abundant sclerotia by the S morphotype suggests adaptation for long-term survival in the soil, whereas the production of abundant conidia by the L morphotype suggests adaptation for aerial dispersal to the phyllosphere. To identify genomic changes that support differential niche adaption, the sequences of three S and three L morphotype isolates were compared. Differences in genome structure and gene content were identified between the morphotypes. A >530 kb inversion between the morphotypes affect a secondary metabolite gene cluster and a cutinase gene. The morphotypes also differed in proteins predicted to be involved in carbon/nitrogen metabolism, iron acquisition, antimicrobial defense, and evasion of host immunity. The S morphotype genomes contained more intact secondary metabolite clusters indicating there is higher selection pressure to maintain secondary metabolism in the soil and that it is not limited to aflatoxin production. The L morphotype genomes were enriched in amino acid transporters, suggesting efficient nitrogen transport may be critical in the nutrient limited phyllosphere. These findings indicate the genomes of the two morphotypes differ beyond developmental genes and have diverged as they adapted to their respective niches.

High intraspecific genome diversity in the model arbuscular mycorrhizal symbiont Rhizophagus irregularis

Arbuscular mycorrhizal fungi (AMF) are known to improve plant fitness through the establishment of mycorrhizal symbioses. Genetic and phenotypic variations among closely related AMF isolates can significantly affect plant growth, but the genomic changes underlying this variability are unclear. To address this issue, we improved the genome assembly and gene annotation of the model strain Rhizophagus irregularis DAOM197198, and compared its gene content with five isolates of R. irregularis sampled in the same field. All isolates harbor striking genome variations, with large numbers of isolate-specific genes, gene family expansions, and evidence of interisolate genetic exchange. The observed variability affects all gene ontology terms and PFAM protein domains, as well as putative mycorrhiza-induced small secreted effector-like proteins and other symbiosis differentially expressed genes. High variability is also found in active transposable elements. Overall, these findings indicate a substantial divergence in the functioning capacity of isolates harvested from the same field, and thus their genetic potential for adaptation to biotic and abiotic changes. Our data also provide a first glimpse into the genome diversity that resides within natural populations of these symbionts, and open avenues for future analyses of plant-AMF interactions that link AMF genome variation with plant phenotype and fitness.

Phytophthora betacei, a new species within Phytophthora clade 1c causing late blight on Solanum betaceum in Colombia

Over the past few years, symptoms akin to late blight disease have been reported on a variety of crop plants in South America. Despite the economic importance of these crops, the causal agents of the diseases belonging to the genus Phytophthora have not been completely characterized. In this study, a new Phytophthora species was described in Colombia from tree tomato (Solanum betaceum), a semi-domesticated fruit grown in northern South America. Comprehensive phylogenetic, morphological, population genetic analyses, and infection assays to characterize this new species, were conducted. All data support the description of the new species, Phytophthora betacei sp. nov. Phylogenetic analyses suggest that this new species belongs to clade 1c of the genus Phytophthora and is a close relative of the potato late blight pathogen, P. infestans. Furthermore, it appeared as the sister group of the P. andina strains collected from wild Solanaceae (clonal lineage EC-2). Analyses of morphological and physiological characters as well as host specificity showed high support for the differentiation of these species. Based on these results, a complete description of the new species is provided and the species boundaries within Phytophthora clade 1c in northern South America are discussed.

Comparative whole-genome analysis reveals artificial selection effects on Ustilago esculenta genome

Ustilago esculenta, infects Zizania latifolia, and induced host stem swollen to be a popular vegetable called Jiaobai in China. It is the long-standing artificial selection that maximizes the occurrence of favourable Jiaobai, and thus maintaining the plant-fungi interaction and modulating the fungus evolving from plant pathogen to entophyte. In this study, whole genome of U. esculenta was sequenced and transcriptomes of the fungi and its host were analysed. The 20.2 Mb U. esculenta draft genome of 6,654 predicted genes including mating, primary metabolism, secreted proteins, shared a high similarity to related Smut fungi. But U. esculenta prefers RNA silencing not repeat-induced point in defence and has more introns per gene, indicating relatively slow evolution rate. The fungus also lacks some genes in amino acid biosynthesis pathway which were filled by up-regulated host genes and developed distinct amino acid response mechanism to balance the infection-resistance interaction. Besides, U. esculenta lost some surface sensors, important virulence factors and host range-related effectors to maintain the economic endophytic life. The elucidation of the U. esculenta genomic information as well as expression profiles can not only contribute to more comprehensive insights into the molecular mechanism underlying artificial selection but also into smut fungi-host interactions.

Population genomic footprints of host adaptation, introgression and recombination in coffee leaf rust

Coffee leaf rust, caused by Hemileia vastatrix (Hv), represents the biggest threat to coffee production worldwide and ranks amongst the most serious fungal diseases in history. Despite a recent series of outbreaks and emergence of hypervirulent strains, the population evolutionary history and potential of this pathogen remain poorly understood. To address this issue, we used restriction site-associated DNA sequencing (RADseq) to generate ∼19 000 single nucleotide polymorphisms (SNPs) across a worldwide collection of 37 Hv samples. Contrary to the long-standing idea that Hv represents a genetically unstructured and cosmopolitan species, our results reveal the existence of a cryptic species complex with marked host tropism. Using phylogenetic and pathological data, we show that one of these lineages (C3) infects almost exclusively the most economically valuable coffee species (tetraploids that include Coffea arabica and interspecific hybrids), whereas the other lineages (C1 and C2) are severely maladapted to these hosts, but successfully infect diploid coffee species. Population dynamic analyses suggest that the C3 group may be a recent 'domesticated' lineage that emerged via host shift from diploid coffee hosts. We also found evidence of recombination occurring within this group, which could explain the high pace of pathotype emergence despite the low genetic variation. Moreover, genomic footprints of introgression between the C3 and C2 groups were discovered and raise the possibility that virulence factors may be quickly exchanged between groups with different pathogenic abilities. This work advances our understanding of the evolutionary strategies used by plant pathogens in agro-ecosystems with direct and far-reaching implications for disease control.

Recent admixture between species of the fungal pathogen Histoplasma

Hybridization between species of pathogens has the potential to speed evolution of virulence by providing the raw material for adaptation through introgression or by assembling new combinations of virulence traits. Fungal diseases are a source high morbidity, and remain difficult to treat. Yet the frequency of hybridization between fungal species has rarely been explored, and the functional role of introgressed alleles remains largely unknown. Histoplasma mississippiense and H. ohiense are sympatric throughout their range in North America and have distinct virulence strategies, making them an ideal system to examine the role introgression may play in fungal pathogens. We identified introgressed tracts in the genomes of a sample of H. mississippiense and H. ohiense isolates. We found strong evidence in each species for recent admixture, but introgressed alleles were present at low frequencies, suggesting that they were deleterious. Consistent with this, coding and regulatory sequences were strongly depleted within introgressed regions, whereas intergenic regions were enriched, indicating that functional introgressed alleles were frequently deleterious in their new genomic context. Surprisingly, we found only two isolates with substantial admixture: the H. mississippiense and H. ohiense genomic reference strains, WU24 and G217B, respectively. Our results show that recent admixture has occurred, that it is frequently deleterious and that conclusions based on studies of the H. mississippiense and H. ohiense type strains should be revisited with more representative samples from the genus.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

Comparative Genomics of Smut Pathogens: Insights From Orphans and Positively Selected Genes Into Host Specialization

Host specialization is a key evolutionary process for the diversification and emergence of new pathogens. However, the molecular determinants of host range are poorly understood. Smut fungi are biotrophic pathogens that have distinct and narrow host ranges based on largely unknown genetic determinants. Hence, we aimed to expand comparative genomics analyses of smut fungi by including more species infecting different hosts and to define orphans and positively selected genes to gain further insights into the genetics basis of host specialization. We analyzed nine lineages of smut fungi isolated from eight crop and non-crop hosts: maize, barley, sugarcane, wheat, oats, Zizania latifolia (Manchurian rice), Echinochloa colona (a wild grass), and Persicaria sp. (a wild dicot plant). We assembled two new genomes: Ustilago hordei (strain Uhor01) isolated from oats and U. tritici (strain CBS 119.19) isolated from wheat. The smut genomes were of small sizes, ranging from 18.38 to 24.63 Mb. U. hordei species experienced genome expansions due to the proliferation of transposable elements and the amount of these elements varied among the two strains. Phylogenetic analysis confirmed that Ustilago is not a monophyletic genus and, furthermore, detected misclassification of the U. tritici specimen. The comparison between smut pathogens of crop and non-crop hosts did not reveal distinct signatures, suggesting that host domestication did not play a dominant role in shaping the evolution of smuts. We found that host specialization in smut fungi likely has a complex genetic basis: different functional categories were enriched in orphans and lineage-specific selected genes. The diversification and gain/loss of effector genes are probably the most important determinants of host specificity.

Population genomic analysis of the rice blast fungus reveals specific events associated with expansion of three main clades

We examined the genomes of 100 isolates of Magnaporthe oryzae (Pyricularia oryzae), the causal agent of rice blast disease. We grouped current field populations of M. oryzae into three major globally distributed groups. A genetically diverse group, clade 1, which may represent a group of closely related lineages, contains isolates of both mating types. Two well-separated clades, clades 2 and 3, appear to have arisen as clonal lineages distinct from the genetically diverse clade. Examination of genes involved in mating pathways identified clade-specific diversification of several genes with orthologs involved in mating behavior in other fungi. All isolates within each clonal lineage are of the same mating type. Clade 2 is distinguished by a unique deletion allele of a gene encoding a small cysteine-rich protein that we determined to be a virulence factor. Clade 3 isolates have a small deletion within the MFA2 pheromone precursor gene, and this allele is shared with an unusual group of isolates we placed within clade 1 that contain AVR1-CO39 alleles. These markers could be used for rapid screening of isolates and suggest specific events in evolution that shaped these populations. Our findings are consistent with the view that M. oryzae populations in Asia generate diversity through recombination and may have served as the source of the clades 2 and 3 isolates that comprise a large fraction of the global population.

Inheritance of phenotypic traits in the progeny of a Ceratocystis interspecific cross

Ceratocystis fimbriata is a fungal plant pathogen that causes black rot on Ipomoea batatas. Based on inoculation studies on numerous tree species, the pathogen is known to be host specific. The closely related species, Ceratocystis manginecans, causes severe wilt on a broad range of tree hosts, including Mangifera indica, Acacia mangium and other leguminous tree species. The genetic factors underlying the pathogenicity and host specificity of Ceratocystis species have rarely been investigated. In this study, an F₁ population of 70 recombinant progeny from a cross between C. fimbriata and C. manginecans was generated and the inheritance of various phenotypic traits was investigated. Results showed that colony colour, growth rate, asexual spore production and aggressiveness to I. batatas and A. mangium are all quantitative traits with high levels of heritability. However, conidia production and aggressiveness appeared to be regulated by a small number of genes. No correlation could be found between aggressiveness and other phenotypic traits, suggesting that these are inherited independently. This is the first study to consider genetic inheritance of pathogenicity and host specificity in Ceratocystis species and the results will contribute, in future, to the identification of quantitative trait loci and candidate genes associated with the traits investigated.

How to boost marine fungal research: A first step towards a multidisciplinary approach by combining molecular fungal ecology and natural products chemistry

Marine fungi have attracted attention in recent years due to increased appreciation of their functional role in ecosystems and as important sources of new natural products. The concomitant development of various "omic" technologies has boosted fungal research in the fields of biodiversity, physiological ecology and natural product biosynthesis. Each of these research areas has its own research agenda, scientific language and quality standards, which have so far hindered an interdisciplinary exchange. Inter- and transdisciplinary interactions are, however, vital for: (i) a detailed understanding of the ecological role of marine fungi, (ii) unlocking their hidden potential for natural product discovery, and (iii) designing access routes for biotechnological production. In this review and opinion paper, we describe the two different "worlds" of marine fungal natural product chemists and marine fungal molecular ecologists. The individual scientific approaches and tools employed are summarised and explained, and enriched with a first common glossary. We propose a strategy to find a multidisciplinary approach towards a comprehensive view on marine fungi and their chemical potential.

Emerging Fungal Threats to Plants and Animals Challenge Agriculture and Ecosystem Resilience

While fungi can make positive contributions to ecosystems and agro-ecosystems, for example, in mycorrhizal associations, they can also have devastating impacts as pathogens of plants and animals. In undisturbed ecosystems, most such negative interactions will be limited through the coevolution of fungi with their hosts. In this article, we explore what happens when pathogenic fungi spread beyond their natural ecological range and become invasive on naïve hosts in new ecosystems. We will see that such invasive pathogens have been problematic to humans and their domesticated plant and animal species throughout history, and we will discuss some of the most pressing fungal threats of today.

Careful with That Axe, Gene, Genome Perturbation after a PEG-Mediated Protoplast Transformation in Fusarium verticillioides

Fusarium verticillioides causes ear rot disease in maize and its contamination with fumonisins, mycotoxins harmful for humans and livestock. Lipids, and their oxidized forms, may drive the fate of this disease. In a previous study, we have explored the role of oxylipins in this interaction by deleting by standard transformation procedures a linoleate diol synthase-coding gene, lds1, in F. verticillioides. A profound phenotypic diversity in the mutants generated has prompted us to investigate more deeply the whole genome of two lds1-deleted strains. Bioinformatics analyses pinpoint significant differences in the genome sequences emerged between the wild type and the lds1-mutants further than those trivially attributable to the deletion of the lds1 locus, such as single nucleotide polymorphisms, small deletion/insertion polymorphisms and structural variations. Results suggest that the effect of a (theoretically) punctual transformation event might have enhanced the natural mechanisms of genomic variability and that transformation practices, commonly used in the reverse genetics of fungi, may potentially be responsible for unexpected, stochastic and henceforth off-target rearrangements throughout the genome.

Co-occurrence and hybridization of anther-smut pathogens specialized on Dianthus hosts

Host specialization has important consequences for the diversification and ecological interactions of obligate pathogens. The anther-smut disease of natural plant populations, caused by Microbotryum fungi, has been characterized by specialized host-pathogen interactions, which contribute in part to the isolation among these numerous fungal species. This study investigated the molecular variation of Microbotryum pathogens within the geographic and host-specific distributions on wild Dianthus species in southern European Alps. In contrast to prior studies on this pathogen genus, a range of overlapping host specificities was observed for four delineated Microbotryum lineages on Dianthus hosts, and their frequent co-occurrence within single-host populations was quantified at local and regional scales. In addition to potential consequences for direct pathogen competition, the sympatry of Microbotryum lineages led to hybridization between them in many populations, and these admixed genotypes suffered significant meiotic sterility. Therefore, this investigation of the anther-smut fungi reveals how variation in the degrees of host specificity can have major implications for ecological interactions and genetic integrity of differentiated pathogen lineages.