Evolution of virulence in fungal plant pathogens: exploiting fungal genomics to control plant disease

Large-scale genomic analyses with machine learning uncover predictive patterns associated with fungal phytopathogenic lifestyles and traits

Invasive plant pathogenic fungi have a global impact, with devastating economic and environmental effects on crops and forests. Biosurveillance, a critical component of threat mitigation, requires risk prediction based on fungal lifestyles and traits. Recent studies have revealed distinct genomic patterns associated with specific groups of plant pathogenic fungi. We sought to establish whether these phytopathogenic genomic patterns hold across diverse taxonomic and ecological groups from the Ascomycota and Basidiomycota, and furthermore, if those patterns can be used in a predictive capacity for biosurveillance. Using a supervised machine learning approach that integrates phylogenetic and genomic data, we analyzed 387 fungal genomes to test a proof-of-concept for the use of genomic signatures in predicting fungal phytopathogenic lifestyles and traits during biosurveillance activities. Our machine learning feature sets were derived from genome annotation data of carbohydrate-active enzymes (CAZymes), peptidases, secondary metabolite clusters (SMCs), transporters, and transcription factors. We found that machine learning could successfully predict fungal lifestyles and traits across taxonomic groups, with the best predictive performance coming from feature sets comprising CAZyme, peptidase, and SMC data. While phylogeny was an important component in most predictions, the inclusion of genomic data improved prediction performance for every lifestyle and trait tested. Plant pathogenicity was one of the best-predicted traits, showing the promise of predictive genomics for biosurveillance applications. Furthermore, our machine learning approach revealed expansions in the number of genes from specific CAZyme and peptidase families in the genomes of plant pathogens compared to non-phytopathogenic genomes (saprotrophs, endo- and ectomycorrhizal fungi). Such genomic feature profiles give insight into the evolution of fungal phytopathogenicity and could be useful to predict the risks of unknown fungi in future biosurveillance activities.

Methionine biosynthesis enzyme MoMet2 is required for rice blast fungus pathogenicity by promoting virulence gene expression via reducing 5mC modification

The emergence of fungicide resistance severely threatens crop production by limiting the availability and application of established fungicides. Therefore, it is urgent to identify new fungicidal targets for controlling plant diseases. Here, we characterized the function of a conserved homoserine O-acetyltransferase (HOA) from the rice blast fungus Magnaporthe oryzae that could serve as the candidate antifungal target. Deletion of the MoMET2 and MoCYS2 genes encoding HOAs perturbed the biosynthesis of methionine and S-adenyl methionine, a methyl group donor for epigenetic modifications, and severely attenuated the development and virulence of M. oryzae. The ∆Momet2 mutant is significantly increased in 5-methylcytosine (5mC) modification that represses the expression of genes required for pathogenicity, including MoGLIK and MoCDH-CYT. We further showed that host-induced gene silencing (HIGS) targeting MoMET2 and MoCYS2 effectively controls rice blasts. Our studies revealed the importance of HOA in the development and virulence of M. oryzae, which suggests the potential feasibility of HOA as new targets for novel anti-rice blast measurements.

Virulence Profiles and Genome-Wide Association Study for Ascochyta lentis Isolates Collected from Australian Lentil-Growing Regions

Ascochyta lentis, the causal organism of Ascochyta blight (AB) of lentil (Lens culinaris), has been shown to produce an avirulence effector protein that mediates AB resistance in certain lentil cultivars. The two known forms of the effector protein were identified from a biparental mapping population between isolates that have reciprocal virulence on 'PBA Hurricane XT' and 'Nipper'. The effector AlAvr1-1 was described for the PBA Hurricane XT-avirulent isolate P94-24 and AlAvr1-2 characterized in the PBA Hurricane XT-virulent isolate AlKewell. Here, we performed a genome-wide association study to identify other loci associated with AB for a differential set of lentil cultivars from a diverse panel of isolates collected in the Australian lentil-growing regions from 2013 to 2020. The chromosome 3 AlAvr1 locus was strongly associated with the PBA Hurricane XT, 'Indianhead', and Nipper disease responses, but one other genomic region on chromosome 11 was also associated with the Nipper disease trait. Our results corroborate earlier work that identified the AlAvr1 locus for field-collected isolates that span the period before release and after widespread adoption of PBA Hurricane XT. A multiplex PCR assay was developed to differentiate the genes AlAvr1-1 and AlAvr1-2 to predict PBA Hurricane XT avirulence and pathotype designation in the diversity panel. Increasing numbers of the PBA Hurricane XT-virulent pathotype 2 isolates across that time indicate strong selection for isolates with the AlAvr1-2 allele. Furthermore, one other region of the A. lentis genome may contribute to the pathogen-host interaction for lentil AB.

Genome Analysis of the Broad Host Range Necrotroph Nalanthamala psidii Highlights Genes Associated With Virulence

Guava wilt disease is caused by the fungus Nalanthamala psidii. The wilt disease results in large-scale destruction of orchards in South Africa, Taiwan, and several Southeast Asian countries. De novo assembly, annotation, and in-depth analysis of the N. psidii genome were carried out to facilitate the identification of characteristics associated with pathogenicity and pathogen evolution. The predicted secretome revealed a range of CAZymes, proteases, lipases and peroxidases associated with plant cell wall degradation, nutrient acquisition, and disease development. Further analysis of the N. psidii carbohydrate-active enzyme profile exposed the broad-spectrum necrotrophic lifestyle of the pathogen, which was corroborated by the identification of putative effectors and secondary metabolites with the potential to induce tissue necrosis and cell surface-dependent immune responses. Putative regulatory proteins including transcription factors and kinases were identified in addition to transporters potentially involved in the secretion of secondary metabolites. Transporters identified included important ABC and MFS transporters involved in the efflux of fungicides. Analysis of the repetitive landscape and the detection of mechanisms linked to reproduction such as het and mating genes rendered insights into the biological complexity and evolutionary potential of N. psidii as guava pathogen. Hence, the assembly and annotation of the N. psidii genome provided a valuable platform to explore the pathogenic potential and necrotrophic lifestyle of the guava wilt pathogen.

The Wheat Stem Rust (Puccinia graminis f. sp. tritici) Population from Washington Contains the Most Virulent Isolates Reported on Barley

A diverse sexual population of wheat stem rust, Puccinia graminis f. sp. tritici, exists in the Pacific Northwest region of the United States because of the natural presence of Mahonia spp. that serves as alternate hosts to complete its sexual life cycle. The region appears to be a center of stem rust diversity in North America where novel virulence gene combinations can emerge that could overcome deployed barley and wheat stem rust resistances. A total of 100 single pustule isolates derived from stem rust samples collected from barley in Eastern Washington during the 2019 growing season were assayed for virulence on the two known effective barley stem rust resistance genes/loci, Rpg1 and the rpg4/5-mediated resistance locus (RMRL) at the seedling stage. Interestingly, 99% of the P. graminis f. sp. tritici isolates assayed were virulent on barley variety Morex carrying the Rpg1 gene, and 62% of the isolates were virulent on the variety Golden Promise transformant (H228.2c) that carries a single-copy insertion of the Rpg1 gene from Morex and is more resistant than Morex to many Rpg1 avirulent isolates. Also, 16% of the isolates were virulent on the near isogenic line HQ-1, which carries the RMRL introgression from the barley line Q21861 in the susceptible Harrington background. Alarmingly, 10% of the isolates were virulent on barley line Q21861, which contains both Rpg1 and RMRL. Thus, we report on the first P. graminis f. sp. tritici isolates worldwide with virulence on both Rpg1 and RMRL when stacked together, representing the most virulent P. graminis f. sp. tritici isolates reported on barley.

Identification of Putative Virulence Genes by DNA Methylation Studies in the Cereal Pathogen Fusarium graminearum

DNA methylation mediates organisms’ adaptations to environmental changes in a wide range of species. We investigated if a such a strategy is also adopted by Fusarium graminearum in regulating virulence toward its natural hosts. A virulent strain of this fungus was consecutively sub-cultured for 50 times (once a week) on potato dextrose agar. To assess the effect of subculturing on virulence, wheat seedlings and heads (cv. A416) were inoculated with subcultures (SC) 1, 23, and 50. SC50 was also used to re-infect (three times) wheat heads (SC50×3) to restore virulence. In vitro conidia production, colonies growth and secondary metabolites production were also determined for SC1, SC23, SC50, and SC50×3. Seedling stem base and head assays revealed a virulence decline of all subcultures, whereas virulence was restored in SC50×3. The same trend was observed in conidia production. The DNA isolated from SC50 and SC50×3 was subject to a methylation content-sensitive enzyme and double-digest, restriction-site-associated DNA technique (ddRAD-MCSeEd). DNA methylation analysis indicated 1024 genes, whose methylation levels changed in response to the inoculation on a healthy host after subculturing. Several of these genes are already known to be involved in virulence by functional analysis. These results demonstrate that the physiological shifts following sub-culturing have an impact on genomic DNA methylation levels and suggest that the ddRAD-MCSeEd approach can be an important tool for detecting genes potentially related to fungal virulence.

DIM5/KMT1 controls fungal insect pathogenicity and genome stability by methylation of histone H3K4, H3K9 and H3K36

Mono-, di- and tri-methylation of histone H3 Lys 9, Lys 4, and Lys 36 (H3K_me1/me2/me3) required for mediation of DNA-based cellular events in eukaryotes usually rely upon the activities of histone lysine methyltransferases (KMTs) classified to the KMT1, KMT2, and KMT3 families, respectively. Here, an H3K9-specific DIM5/KMT1 orthologue, which lacks a C-terminal post-SET domain and localizes mainly in nucleus, is reported to have both conserved and noncanonical roles in methylating the H3 core lysines in Beauveria bassiana, an insect-pathogenic fungus serving as a main source of wide-spectrum fungal insecticides. Disruption of dim5 led to abolishment of H3K9me3 and marked attenuation of H3K4me1/me2, H3K9me1/me2 and H3K36me2. Consequently, the Δdim5 mutant lost the whole insect pathogenicity through normal cuticle infection, and was compromised severely in virulence through cuticle-bypassing infection (hemocoel injection) and also in a series of cellular events critical for the fungal virulence and lifecycle in vivo and in vitro, including reduced hyphal growth, blocked conidiation, impeded proliferation in vivo, altered carbohydrate epitopes, disturbed cell cycle, reduced biosynthesis and secretion of cuticle-degrading enzymes, and increased sensitivities to various stresses. Among 1,201 dysregulated genes (up/down ratio: 712:489) associated with those phenotypic changes, 92 (up/down ratio: 59:33) encode transcription factors and proteins or enzymes involved in posttranslational modifications, implying that the DIM5-methylated H3 core lysines could act as preferential marks of those transcription-active genes crucial for global gene regulation. These findings uncover a novel scenario of DIM5 and its indispensability for insect-pathogenic lifestyle and genome stability of B. bassiana.

The potential use of the Penicillium chrysogenum antifungal protein PAF, the designed variant PAFopt and its γ-core peptide Pγopt in plant protection

The prevention of enormous crop losses caused by pesticide-resistant fungi is a serious challenge in agriculture. Application of alternative fungicides, such as antifungal proteins and peptides, provides a promising basis to overcome this problem; however, their direct use in fields suffers limitations, such as high cost of production, low stability, narrow antifungal spectrum and toxicity on plant or mammalian cells. Recently, we demonstrated that a Penicillium chrysogenum-based expression system provides a feasible tool for economic production of P. chrysogenum antifungal protein (PAF) and a rational designed variant (PAFopt ), in which the evolutionary conserved γ-core motif was modified to increase antifungal activity. In the present study, we report for the first time that γ-core modulation influences the antifungal spectrum and efficacy of PAF against important plant pathogenic ascomycetes, and the synthetic γ-core peptide Pγopt , a derivative of PAFopt , is antifungal active against these pathogens in vitro. Finally, we proved the protective potential of PAF against Botrytis cinerea infection in tomato plant leaves. The lack of any toxic effects on mammalian cells and plant seedlings, as well as the high tolerance to harsh environmental conditions and proteolytic degradation further strengthen our concept for applicability of these proteins and peptide in agriculture.

The Pattern and Function of DNA Methylation in Fungal Plant Pathogens

To successfully infect plants and trigger disease, fungal plant pathogens use various strategies that are dependent on characteristics of their biology and genomes. Although pathogenic fungi are different from animals and plants in the genomic heritability, sequence feature, and epigenetic modification, an increasing number of phytopathogenic fungi have been demonstrated to share DNA methyltransferases (MTases) responsible for DNA methylation with animals and plants. Fungal plant pathogens predominantly possess four types of DNA MTase homologs, including DIM-2, DNMT1, DNMT5, and RID. Numerous studies have indicated that DNA methylation in phytopathogenic fungi mainly distributes in transposable elements (TEs), gene promoter regions, and the repetitive DNA sequences. As an important and heritable epigenetic modification, DNA methylation is associated with silencing of gene expression and transposon, and it is responsible for a wide range of biological phenomena in fungi. This review highlights the relevant reports and insights into the important roles of DNA methylation in the modulation of development, pathogenicity, and secondary metabolism of fungal plant pathogens. Recent evidences prove that there are massive links between DNA and histone methylation in fungi, and they commonly regulate fungal development and mycotoxin biosynthesis.

Colletotrichum: species complexes, lifestyle, and peculiarities of some sources of genetic variability

The genus Colletotrichum comprises species with different lifestyles but is mainly known for phytopathogenic species that infect crops of agronomic relevance causing considerable losses. The fungi of the genus Colletotrichum are distributed in species complexes and within each complex some species have particularities regarding their lifestyle. The most commonly found and described lifestyles in Colletotrichum are endophytic and hemibiotrophic phytopathogenic. Several of these phytopathogenic species show wide genetic variability, which makes long-term maintenance of resistance in plants difficult. Different mechanisms may play an important role in the emergence of genetic variants but are not yet fully understood in this genus. These mechanisms include heterokaryosis, a parasexual cycle, sexual cycle, transposable element activity, and repeat-induced point mutations. This review provides an overview of the genus Colletotrichum, the species complexes described so far and the most common lifestyles in the genus, with a special emphasis on the mechanisms that may be responsible, at least in part, for the emergence of new genotypes under field conditions.

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.

The emergence of the multi-species NIP1 effector in Rhynchosporium was accompanied by high rates of gene duplications and losses

Plant pathogens secrete effector proteins to manipulate the host and facilitate infection. Cognate hosts trigger strong defence responses upon detection of these effectors. Consequently, pathogens and hosts undergo rapid coevolutionary arms races driven by adaptive evolution of effectors and receptors. Because of their high rate of turnover, most effectors are thought to be species-specific and the evolutionary trajectories are poorly understood. Here, we investigate the necrosis-inducing protein 1 (NIP1) effector in the multihost pathogen genus Rhynchosporium. We retraced the evolutionary history of the NIP1 locus using whole-genome assemblies of 146 strains covering four closely related species. NIP1 orthologues were present in all species but the locus consistently segregated presence-absence polymorphisms suggesting long-term balancing selection. We also identified previously unknown paralogues of NIP1 that were shared among multiple species and showed substantial copy-number variation within R. commune. The NIP1A paralogue was under significant positive selection suggesting that NIP1A is the dominant effector variant coevolving with host immune receptors. Consistent with this prediction, we found that copy number variation at NIP1A had a stronger effect on virulence than NIP1B. Our analyses unravelled the origins and diversification mechanisms of a pathogen effector family shedding light on how pathogens gain adaptive genetic variation.

Pyramiding for Resistance Durability: Theory and Practice

Durable disease resistance is a key component of global food security, and combining resistance genes into "pyramids" is an important way to increase durability of resistance. The mechanisms by which pyramids impart durability are not well known. The traditional view of resistance pyramids considers the use of major resistance gene (R-gene) combinations deployed against pathogens that are primarily asexual. Interestingly, published examples of the successful use of pyramids in the traditional sense are rare. In contrast, most published descriptions of durable pyramids in practice are for cereal rusts, and tend to indicate an association between durability and cultivars combining major R-genes with incompletely expressed, adult plant resistance genes. Pyramids have been investigated experimentally for a diversity of pathogens, and many reduce disease levels below that of the single best gene. Resistance gene combinations have been identified through phenotypic reactions, molecular markers, and challenge against effector genes. As resistance genes do not express equally in all genetic backgrounds, however, a combination of genetic information and phenotypic analyses provide the ideal scenario for testing of putative pyramids. Not all resistance genes contribute equally to pyramids, and approaches have been suggested to identify the best genes and combinations of genes for inclusion. Combining multiple resistance genes into a single plant genotype quickly is a challenge that is being addressed through alternative breeding approaches, as well as through genomics tools such as resistance gene cassettes and gene editing. Experimental and modeling tests of pyramid durability are in their infancy, but have promise to help direct future studies of pyramids. Several areas for further work on resistance gene pyramids are suggested.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Fungal Plant Pathogenesis Mediated by Effectors

The interactions between fungi and plants encompass a spectrum of ecologies ranging from saprotrophy (growth on dead plant material) through pathogenesis (growth of the fungus accompanied by disease on the plant) to symbiosis (growth of the fungus with growth enhancement of the plant). We consider pathogenesis in this article and the key roles played by a range of pathogen-encoded molecules that have collectively become known as effectors.

Pterostilbene Is a Potential Candidate for Control of Blackleg in Canola

Two stilbenes, resveratrol and pterostilbene, exhibit antifungal activity against Leptosphaeria maculans, the fungal pathogen responsible for blackleg (stem canker) in canola (Brassica napus). In vitro studies on the effect of these stilbenes on L. maculans mycelial growth and conidia germination showed that pterostilbene is a potent fungicide and sporicide, but resveratrol only exerted minor inhibition on L. maculans. Cell viability of hyphae cultures was markedly reduced by pterostilbene and SYTOX green staining showed that cell membrane integrity was compromised. We demonstrate that pterostilbene exerts fungicidal activity across 10 different L. maculans isolates and the compound confers protection to the blackleg-susceptible canola cv. Westar seedlings. The potential of pterostilbene as a control agent against blackleg in canola is discussed.

Comparative genome analysis and genome evolution of members of the magnaporthaceae family of fungi

BACKGROUND

Magnaporthaceae, a family of ascomycetes, includes three fungi of great economic importance that cause disease in cereal and turf grasses: Magnaporthe oryzae (rice blast), Gaeumannomyces graminis var. tritici (take-all disease), and Magnaporthe poae (summer patch disease). Recently, the sequenced and assembled genomes for these three fungi were reported. Here, the genomes were compared for orthologous genes in order to identified genes that are unique to the Magnaporthaceae family of fungi. In addition, ortholog clustering was used to identify a core proteome for the Magnaporthaceae, which was examined for diversifying and purifying selection and evidence of two-speed genome evolution.

RESULTS

A genome-scale comparative study was conducted across 74 fungal genomes to identify clusters of orthologous genes unique to the three Magnaporthaceae species as well as species specific genes. We found 1149 clusters that were unique to the Magnaporthaceae family of fungi with 295 of those containing genes from all three species. Gene clusters involved in metabolic and enzymatic activities were highly represented in the Magnaporthaceae specific clusters. Also highly represented in the Magnaporthaceae specific clusters as well as in the species specific genes were transcriptional regulators. In addition, we examined the relationship between gene evolution and distance to repetitive elements found in the genome. No correlations between diversifying or purifying selection and distance to repetitive elements or an increased rate of evolution in secreted and small secreted proteins were observed.

CONCLUSIONS

Taken together, these data show that at the genome level, there is no evidence to suggest multi-speed genome evolution or that proximity to repetitive elements play a role in diversification of genes.