Finished genome of the fungal wheat pathogen Mycosphaerella graminicola reveals dispensome structure, chromosome plasticity, and stealth pathogenesis

Genome evolution in intracellular parasites: Microsporidia and Apicomplexa

Microsporidia and Apicomplexa are eukaryotic, single-celled, intracellular parasites with huge public health and economic importance. Typically, these parasites are studied separately, emphasizing their uniqueness and diversity. In this review, we explore the huge amount of genomic data that has recently become available for the two groups. We compare and contrast their genome evolution and discuss how their transitions to intracellular life may have shaped it. In particular, we explore genome reduction and compaction, genome expansion and ploidy, gene shuffling and rearrangements, and the evolution of centromeres and telomeres.

Defense Pathways of Wheat Plants Inoculated with Zymoseptoria tritici under NaCl Stress Conditions: An Overview

Due to being sessile, plants develop a broad range of defense pathways when they face abiotic or biotic stress factors. Although plants are subjected to more than one type of stress at a time in nature, the combined effects of either multiple stresses of one kind (abiotic or biotic) or more kinds (abiotic and biotic) have now been realized in agricultural lands due to increases in global warming and environmental pollution, along with population increases. Soil-borne pathogens, or pathogens infecting aerial parts, can have devastating effects on plants when combined with other stressors. Obtaining yields or crops from sensitive or moderately resistant plants could be impossible, and it could be very difficult from resistant plants. The mechanisms of combined stress in many plants have previously been studied and elucidated. Recent studies proposed new defense pathways and mechanisms through signaling cascades. In light of these mechanisms, it is now time to develop appropriate strategies for crop protection under multiple stress conditions. This may involve using disease-resistant or stress-tolerant plant varieties, implementing proper irrigation and drainage practices, and improving soil quality. However, generation of both stress-tolerant and disease-resistant crop plants is of crucial importance. The establishment of a database and understanding of the defense mechanisms under combined stress conditions would be meaningful for the development of resistant and tolerant plants. It is clear that leaf pathogens show great tolerance to salinity stress and result in pathogenicity in crop plants. We noticed that regulation of the stomata through biochemical applications and some effort with the upregulation of the minor gene expressions indirectly involved with the defense mechanisms could be a great way to increase the defense metabolites without interfering with quality parameters. In this review, we selected wheat as a model plant and Zymoseptoria tritici as a model leaf pathogen to evaluate the defense mechanisms under saline conditions through physiological, biochemical, and molecular pathways and suggested various ways to generate tolerant and resistant cereal plants.

The Egyptian wheat cultivar Gemmeiza-12 is a source of resistance against the fungus Zymoseptoria tritici

BACKGROUND

Wheat is one of the world's most important cereal crops. However, the fungal pathogen Zymoseptoria tritici can cause disease epidemics, leading to reduced yields. With climate change and development of new agricultural areas with suitable environments, Z. tritici may advance into geographical areas previously unaffected by this pathogen. It is currently unknown how Egyptian wheat will perform in the face of this incoming threat. This project aimed to assess the resistance of Egyptian wheat germplasm to Z. tritici, to identify cultivars with high levels of resistance and characterise the mechanism(s) of resistance present in these cultivars.

RESULTS

Eighteen Egyptian wheat cultivars were screened against two Z. tritici model isolates and exhibited a wide spectrum of responses. This ranged from resistance to complete susceptibility to one or both isolates tested. The most highly resistant cultivars from the initial screen were then tested under two environmental conditions against modern UK field isolates. Disease levels under UK-like conditions were higher, however, symptom development on the cultivar Gemmeiza-12 was noticeably slower than on other Egyptian wheats. The robustness of the resistance shown by Gemmeiza-12 was confirmed in experiments mimicking Egyptian environmental conditions, where degree of Z. tritici infection was lower. The Kompetitive allele-specific PCR (KASP) diagnostic assay suggested the presence of an Stb6 resistant allele in several Egyptian wheats including Gemmeiza-12. Infection assays using the IPO323 WT and IPO323ΔAvrStb6 mutant confirmed the presence of Stb6 in several Egyptian cultivars including Gemmeiza-12. Confocal fluorescence microscopy demonstrated that growth of the IPO323 strain is blocked at the point of stomatal penetration on Gemmeiza-12, consistent with previous reports of Stb gene mediated resistance. In addition to this R-gene mediated resistance, IPO323 spores showed lower adherence to leaves of Gemmeiza-12 compared to UK wheat varieties, suggesting other aspects of leaf physiology may also contribute to the resistance phenotype of this cultivar.

CONCLUSION

These results indicate that Gemmeiza-12 will be useful in future breeding programs where improved resistance to Z. tritici is a priority.

Guanosine-specific single-stranded ribonuclease effectors of a phytopathogenic fungus potentiate host immune responses

Plants activate immunity upon recognition of pathogen-associated molecular patterns. Although phytopathogens have evolved a set of effector proteins to counteract plant immunity, some effectors are perceived by hosts and induce immune responses. Here, we show that two secreted ribonuclease effectors, SRN1 and SRN2, encoded in a phytopathogenic fungus, Colletotrichum orbiculare, induce cell death in a signal peptide- and catalytic residue-dependent manner, when transiently expressed in Nicotiana benthamiana. The pervasive presence of SRN genes across Colletotrichum species suggested the conserved roles. Using a transient gene expression system in cucumber (Cucumis sativus), an original host of C. orbiculare, we show that SRN1 and SRN2 potentiate host pattern-triggered immunity responses. Consistent with this, C. orbiculare SRN1 and SRN2 deletion mutants exhibited increased virulence on the host. In vitro analysis revealed that SRN1 specifically cleaves single-stranded RNAs at guanosine, leaving a 3'-end phosphate. Importantly, the potentiation of C. sativus responses by SRN1 and SRN2, present in the apoplast, depends on ribonuclease catalytic residues. We propose that the pathogen-derived apoplastic guanosine-specific single-stranded endoribonucleases lead to immunity potentiation in plants.

Preventing multi-resistance: New insights for managing fungal adaptation

Sustainable crop protection is vital for food security, yet it is under threat due to the adaptation of a diverse and evolving pathogen population. Resistance can be managed by maximising the diversity of selection pressure through dose variation and the spatial and temporal combination of active ingredients. This study explores the interplay between operational drivers for maximising the sustainability of management strategies in relation to the resistance status of fungal populations. We applied an experimental evolution approach to three artificial populations of Zymoseptoria tritici, an economically significant wheat pathogen, each differing in initial resistance status. Our findings reveal that diversified selection pressure curtails the selection of resistance in naïve populations and those with low frequencies of single resistance. Increasing the number of modes of action most effectively delays resistance development, surpassing the increase in the number of fungicides, fungicide choice based on resistance risk, and temporal variation in fungicide exposure. However, this approach favours generalism in the evolved populations. The prior presence of multiple resistant isolates and their subsequent selection in populations override the effects of diversity in management strategies, thereby invalidating any universal ranking. Therefore, the initial resistance composition must be specifically considered in sustainable resistance management to address real-world field situations.

Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Recent reactivation of a pathogenicity-associated transposable element is associated with major chromosomal rearrangements in a fungal wheat pathogen

Transposable elements (TEs) are key drivers of genomic variation contributing to recent adaptation in most species. Yet, the evolutionary origins and insertion dynamics within species remain poorly understood. We recapitulate the spread of the pathogenicity-associated Styx element across five species that last diverged ∼11 000 years ago. We show that the element likely originated in the Zymoseptoria fungal pathogen genus and underwent multiple independent reactivation events. Using a global 900-genome panel of the wheat pathogen Zymoseptoria tritici, we assess Styx copy number variation and identify renewed transposition activity in Oceania and South America. We show that the element can mobilize to create additional Styx copies in a four-generation pedigree. Importantly, we find that new copies of the element are not affected by genomic defenses suggesting minimal control against the element. Styx copies are preferentially located in recombination breakpoints and likely triggered multiple types of large chromosomal rearrangements. Taken together, we establish the origin, diversification and reactivation of a highly active TE with likely major consequences for chromosomal integrity and the expression of disease.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Use of chitin:DNA ratio to assess growth form in fungal cells

BACKGROUND

Dimorphism, the ability to switch between a 'yeast-like' and a hyphal growth form, is an important feature of certain fungi, including important plant and human pathogens. The switch to hyphal growth is often associated with virulence, pathogenicity, biofilm formation and stress resistance. Thus, the ability to accurately and efficiently measure fungal growth form is key to research into these fungi, especially for discovery of potential drug targets. To date, fungal growth form has been assessed microscopically, a process that is both labour intensive and costly.

RESULTS

Here, we unite quantification of the chitin in fungal cell walls and the DNA in nuclei to produce a methodology that allows fungal cell shape to be estimated by calculation of the ratio between cell wall quantity and number of nuclei present in a sample of fungus or infected host tissue. Using the wheat pathogen Zymoseptoria tritici as a test case, with confirmation in the distantly related Fusarium oxysporum, we demonstrate a close, linear relationship between the chitin:DNA ratio and the average polarity index (length/width) of fungal cells. We show the utility of the method for estimating growth form in infected wheat leaves, differentiating between the timing of germination in two different Z. tritici isolates using this ratio. We also show that the method is robust to the occurrence of thick-walled chlamydospores, which show a chitin:DNA ratio that is distinct from either 'yeast-like' blastospores or hyphae.

CONCLUSIONS

The chitin:DNA ratio provides a simple methodology for determining fungal growth form in bulk tissue samples, reducing the need for labour-intensive microscopic studies requiring specific staining or GFP-tags to visualise the fungus within host tissues. It is applicable to a range of dimorphic fungi under various experimental conditions.

Distribution and conservation of simple sequence repeats in plant pathogenic species of Zymoseptoria and development of genomic resources for its orphaned species

To better understand the structure and evolution of the genomes of four plant pathogenic species of Zymoseptoria, we analyzed the occurrence, relative abundance (RA), and density (RD) of simple sequence repeats (SSRs) in their whole genome and transcriptome sequences. In this study, SSRs are defined as repeats of more than 12 bases in length. The genome and transcriptome sequences of Zymoseptoria ardabiliae show the highest RA (201.1 and 129.9) and RD (3229.4 and 1928.2) of SSRs, while those of Zymoseptoria pseudotritici show the lowest RA (167.2 and 118.5) and RD (2482.2 and 1687.0). The majority of SSRs in the genomic and transcriptome sequences of species were trinucleotide SSRs, while dinucleotide SSRs were the least common. The most common trinucleotide motifs in the transcriptomic sequences across all species were those that encoded the amino acid arginine. As per our motif conservation study, Zymoseptoria tritici (12.4%) possessed the most unique motifs, while Z. pseudotritici (3.9%) had the fewest. Overall, only 38.1% of the motifs were found to be conserved among the species. Gene enrichment studies reveal that three of the species, Z. ardabiliae, Zymoseptoria brevis, and Z. pseudotritici, have SSRs in their genes related to cellular metabolism, while the remaining Z. tritici harbors SSRs in genes related to DNA synthesis and gene expression. In an effort to improve the genetic resources for the orphan species of pathogenic Zymoseptoria, a total of 73,134 primers were created. The genomic resources developed in this study could help with analyses of genetic relatedness within the population and the development of species-specific markers.

Dimethylpolysulfides production as the major mechanism behind wheat fungal pathogen biocontrol, by Arthrobacter and Microbacterium actinomycetes

IMPORTANCE

As the management of wheat fungal diseases becomes increasingly challenging, the use of bacterial agents with biocontrol potential against the two major wheat phytopathogens, Fusarium graminearum and Zymoseptoria tritici, may prove to be an interesting alternative to conventional pest management. Here, we have shown that dimethylpolysulfide volatiles are ubiquitously and predominantly produced by wheat-associated Microbacterium and Arthrobacter actinomycetes, displaying antifungal activity against both pathogens. By limiting pathogen growth and DON virulence factor production, the use of such DMPS-producing strains as soil biocontrol inoculants could limit the supply of pathogen inocula in soil and plant residues, providing an attractive alternative to dimethyldisulfide fumigant, which has many non-targeted toxicities. Notably, this study demonstrates the importance of bacterial volatile organic compound uptake by inhibited F. graminearum, providing new insights for the study of volatiles-mediated toxicity mechanisms within bacteria-fungus signaling crosstalk.

A pangenome-guided manually curated library of transposable elements for Zymoseptoria tritici

OBJECTIVES

High-quality species-specific transposable element (TE) libraries are required for studies to elucidate the evolutionary dynamics of TEs and gain an understanding of their impacts on host genomes. Such high-quality TE resources are severely lacking for species in the fungal kingdom. To facilitate future studies on the putative role of TEs in rapid adaptation observed in the fungal wheat pathogen Zymoseptoria tritici, we produced a manually curated TE library. This was generated by detecting TEs in 19 reference genome assemblies representing the global diversity of the species supplemented by multiple sister species genomes. Improvements over previous TE libraries have been made on TE boundary resolution, detection of ORFs, TE domains, terminal inverted repeats, and class-specific motifs.

DATA DESCRIPTION

A TE consensus library for Z. tritici formatted for use with RepeatMasker. This data is relevant to other researchers investigating TE-host evolutionary dynamics in Z. tritici or who are interested in comparative studies of the fungal kingdom. Further, this TE library can be used to improve gene annotation. Finally, this TE library increases the number of manually curated TE datasets, providing resources to further our understanding of TE diversity.

Combined reference-free and multi-reference based GWAS uncover cryptic variation underlying rapid adaptation in a fungal plant pathogen

Microbial pathogens often harbor substantial functional diversity driven by structural genetic variation. Rapid adaptation from such standing variation threatens global food security and human health. Genome-wide association studies (GWAS) provide a powerful approach to identify genetic variants underlying recent pathogen adaptation. However, the reliance on single reference genomes and single nucleotide polymorphisms (SNPs) obscures the true extent of adaptive genetic variation. Here, we show quantitatively how a combination of multiple reference genomes and reference-free approaches captures substantially more relevant genetic variation compared to single reference mapping. We performed reference-genome based association mapping across 19 reference-quality genomes covering the diversity of the species. We contrasted the results with a reference-free (i.e., k-mer) approach using raw whole-genome sequencing data in a panel of 145 strains collected across the global distribution range of the fungal wheat pathogen Zymoseptoria tritici. We mapped the genetic architecture of 49 life history traits including virulence, reproduction and growth in multiple stressful environments. The inclusion of additional reference genome SNP datasets provides a nearly linear increase in additional loci mapped through GWAS. Variants detected through the k-mer approach explained a higher proportion of phenotypic variation than a reference genome-based approach and revealed functionally confirmed loci that classic GWAS approaches failed to map. The power of GWAS in microbial pathogens can be significantly enhanced by comprehensively capturing structural genetic variation. Our approach is generalizable to a large number of species and will uncover novel mechanisms driving rapid adaptation of pathogens.

The expression landscape and pangenome of long non-coding RNA in the fungal wheat pathogen Zymoseptoria tritici

Long non-coding RNAs (lncRNAs) are regulatory molecules interacting in a wide array of biological processes. lncRNAs in fungal pathogens can be responsive to stress and play roles in regulating growth and nutrient acquisition. Recent evidence suggests that lncRNAs may also play roles in virulence, such as regulating pathogenicity-associated enzymes and on-host reproductive cycles. Despite the importance of lncRNAs, only a few model fungi have well-documented inventories of lncRNA. In this study, we apply a recent computational pipeline to predict high-confidence lncRNA candidates in Zymoseptoria tritici, an important global pathogen of wheat impacting global food production. We analyse genomic features of lncRNAs and the most likely associated processes through analyses of expression over a host infection cycle. We find that lncRNAs are frequently expressed during early infection, before the switch to necrotrophic growth. They are mostly located in facultative heterochromatic regions, which are known to contain many genes associated with pathogenicity. Furthermore, we find that lncRNAs are frequently co-expressed with genes that may be involved in responding to host defence signals, such as oxidative stress. Finally, we assess pangenome features of lncRNAs using four additional reference-quality genomes. We find evidence that the repertoire of expressed lncRNAs varies substantially between individuals, even though lncRNA loci tend to be shared at the genomic level. Overall, this study provides a repertoire and putative functions of lncRNAs in Z. tritici enabling future molecular genetics and functional analyses in an important pathogen.

Chromosome-scale assembly uncovers genomic compartmentation of Fusarium oxysporum f. sp. albedinis, the causal agent of Bayoud disease in date palm

Date palm (Phoenixdactylifera) is the most significant crop across North Africa and the Middle East. However, the crop faces a severe threat from Bayoud disease caused by the fungal pathogen Fusarium oxysporum f. sp. albedinis (FOA). FOA is a soil-borne fungus that infects the roots and vascular system of date palms, leading to widespread destruction of date palm plantations in North Africa over the last century. This is considered the most devastating pathogen of oasis agriculture in North Africa and responsible for loss of 13 million trees in Algeria and Morocco alone. In this study, we present a chromosome-scale high-quality genome assembly of the virulent isolate Foa 44, which provides valuable insights into understanding the genetic basis of Bayoud disease. The genome assembly consists of 11 chromosomes and 40 unplaced contigs, totalling 65,971,825 base pairs in size. It exhibits a GC ratio of 47.77% and a TE (transposable element) content of 17.30%. Through prediction and annotation, we identified 20,416 protein-coding genes. By combining gene and repeat densities analysis with alignment to Fusarium oxysporum f. sp. lycopersici (FOL) 4287 isolate genome sequence, we determined the core and lineage-specific compartments in Foa 44, shedding light on the genome structure of this pathogen. Furthermore, a phylogenomic analysis based on the 3,292 BUSCOs core genome revealed a distinct clade of FOA isolates within the Fusarium oxysporum species complex (FOSC). Notably, the genealogies of the five identified Secreted In Xylem (SIX) genes (1, 6, 9, 11 and 14) in FOA displayed a polyphyletic pattern, suggesting a horizontal inheritance of these effectors. These findings provide a valuable genomics toolbox for further research aimed at combatting the serious biotic constraints posed by FOA to date palm. This will pave the way for a deeper understanding of Bayoud disease and facilitate the development of effective diagnostic tools and control measures.

Epiphytic proliferation of Zymoseptoria tritici isolates on resistant wheat leaves

The wheat pathogen Zymoseptoria tritici is capable of a long period of pre-invasive epiphytic growth. Studies have shown that virulent isolates vary in the extent, duration and growth form of this epiphytic growth, and the fungus has been observed to undergo behaviours such as asexual reproduction by budding and vegetative fusion of hyphae on the leaf surface. This epiphytic colonisation has been investigated very little during interactions in which an isolate of Z. tritici is unable to colonise the apoplast, as occurs during avirulence. However, avirulent isolates have been seen to undergo sexual crosses in the absense of leaf penetration, and it is widely accepted that the main point of distinction between virulent and avirulent isolates occurs at the point of attempted leaf penetration or attempted apoplastic growth, which fails in the avirulent case. In this work, we describe extensive epiphytic growth in three isolates which are unable or have very limited ability to invade the leaf, and show that growth form is as variable as for fully virulent isolates. We demonstrate that during certain interactions, Z. tritici isolates rarely invade the leaf and form pycnidia, but induce necrosis. These isolates are able to achieve higher epiphytic biomass than fully virulent isolates during asymptomatic growth, and may undergo very extensive asexual reproduction on the leaf surface. These findings have implications for open questions such as whether and how Z. tritici obtains nutrients on the leaf surface and the nature of its interaction with wheat defences.

In silico characterization, molecular docking, and dynamic simulation of a novel fungal cell-death suppressing effector, MoRlpA as potential cathepsin B-like cysteine protease inhibitor during rice blast infection

The blast fungus Magnaporthe oryzae is one of the most notorious pathogens affecting rice production worldwide. The cereal killer employs a special class of small secreted proteins called effectors to manipulate and perturb the host metabolism. In turn, the host plants trigger effector-triggered immunity (ETI) via localized cell death and hypersensitive response (HR). We have identified and characterized a novel secreted effector MoRlpA from M. oryzae by extensive in silico methods. The localization studies suggested that it is exclusively secreted in the host apoplasts. Interestingly, MoRlpA interacts with a protease, cathepsin B from rice with highest affinity. The 3D structural models of both the proteins were generated. Cathepsin B-like cysteine proteases are usually involved in programmed cell death (PCD) and autophagy in plants which lead to generation of HR upon infection. Our results suggest that MoRlpA interacts with rice cathepsin B-like cysteine protease and demolish the host counter-attack by suppressing cell death and HR during an active blast infection. This was further validated by molecular docking and molecular dynamic simulation analyses. The important residues involved in the rice-blast pathogen interactions were deciphered. Overall, this research highlights stable interactions between MoRlpA-OsCathB during rice blast pathogenesis and providing an insight into how this novel RlpA protease inhibitor-cum-effector modulates the host's apoplast to invade the host tissues and establish a successful infection. Thus, this research will help to develop potential fungicide to block the binding region of MoRlpA target so that the cryptic pathogen would be recognized by the host. HIGHLIGHTSFor the first time, a novel secreted effector protein, MoRlpA has been identified and characterised from M. oryzae in silicoMoRlpA contains a rare lipoprotein A-like DPBB domain which is often an enzymatic domain in other systemsMoRlpA as an apoplastic effector interacts with the rice protease OsCathB to suppress the cell death and hypersensitive response during rice blast infectionThe three-dimensional structures of both the MoRlpA and OsCathB proteins were predictedMoRlpA-OsCathB interactions were analysed by molecular docking and molecular dynamic simulation studiesCommunicated by Ramaswamy H. Sarma.

Phylogeographic and Phylogenomic Structure of the Quarantine Plant Pathogen Colletotrichum liriopes, Including New Reports in the United States

Global agricultural trade has accelerated the emergence and re-emergence of new plant pathogens. In the United States, the fungal pathogen Colletotrichum liriopes is still considered a foreign quarantine pathogen that affects ornamental plants (i.e., Liriope spp.). Even though this species has been reported in East Asia on various asparagaceous hosts, its first and only report in the United States was in 2018. However, that study used only ITS nrDNA for identification, and no available culture or voucher specimen was maintained. The main objective of the present study was to determine the geographic and host distribution of specimens identified as C. liriopes. To accomplish this, new and existing isolates, sequences, and genomes obtained from various hosts and geographic locations (i.e., China, Colombia, Mexico, and the United States) were compared with the ex-type of C. liriopes. Multilocus phylogenetic (ITS, Tub2, GAPDH, CHS-1, and HIS3), phylogenomic, and splits tree analyses revealed that all the studied isolates/sequences form a well-supported clade with little intraspecific variation. Morphological characterizations support these findings. The minimum spanning network, low nucleotide diversity, and negative Tajima's D from both multilocus and genomic data suggest that there was a recent movement/invasion of a few East Asian genotypes to other countries where the ornamental plants are produced (e.g., South America) and subsequently to the importing countries, such as the United States. The study reveals that the geographic and host distribution of C. liriopes sensu stricto is expanded to the United States (i.e., at least Maryland, Mississippi, and Tennessee) and on various hosts in addition to Asparagaceae and Orchidaceae. The present study produces fundamental knowledge that can be used in efforts to reduce costs or losses from agricultural trade and to expand our understanding of pathogen movement.

Horizontally Transferred DNA in the Genome of the Fungus Pyricularia oryzae is Associated With Repressive Histone Modifications

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

H4K20me3 is important for Ash1-mediated H3K36me3 and transcriptional silencing in facultative heterochromatin in a fungal pathogen

Facultative heterochromatin controls development and differentiation in many eukaryotes. In metazoans, plants, and many filamentous fungi, facultative heterochromatin is characterized by transcriptional repression and enrichment with nucleosomes that are trimethylated at histone H3 lysine 27 (H3K27me3). While loss of H3K27me3 results in derepression of transcriptional gene silencing in many species, additional up- and downstream layers of regulation are necessary to mediate control of transcription in chromosome regions enriched with H3K27me3. Here, we investigated the effects of one histone mark on histone H4, namely H4K20me3, in the fungus Zymoseptoria tritici, a globally important pathogen of wheat. Deletion of kmt5, the gene encoding the sole methyltransferase responsible for H4K20 methylation, resulted in global derepression of transcription, especially in regions of facultative heterochromatin. Derepression in the absence of H4K20me3 not only affected known genes but also a large number of novel, previously undetected transcripts generated from regions of facultative heterochromatin on accessory chromosomes. Transcriptional activation in kmt5 deletion strains was accompanied by a complete loss of Ash1-mediated H3K36me3 and chromatin reorganization affecting H3K27me3 and H3K4me2 distribution in regions of facultative heterochromatin. Strains with H4K20L, M or Q mutations in the single histone H4 gene of Z. tritici recapitulated these chromatin changes, suggesting that H4K20me3 is important for Ash1-mediated H3K36me3. The ∆kmt5 mutants we obtained were more sensitive to genotoxic stressors than wild type and both, ∆kmt5 and ∆ash1, showed greatly increased rates of accessory chromosome loss. Taken together, our results provide insights into an unsuspected mechanism involved in the assembly and maintenance of facultative heterochromatin.

A secreted protease-like protein in Zymoseptoria tritici is responsible for avirulence on Stb9 resistance gene in wheat

Zymoseptoria tritici is the fungal pathogen responsible for Septoria tritici blotch on wheat. Disease outcome in this pathosystem is partly determined by isolate-specific resistance, where wheat resistance genes recognize specific fungal factors triggering an immune response. Despite the large number of known wheat resistance genes, fungal molecular determinants involved in such cultivar-specific resistance remain largely unknown. We identified the avirulence factor AvrStb9 using association mapping and functional validation approaches. Pathotyping AvrStb9 transgenic strains on Stb9 cultivars, near isogenic lines and wheat mapping populations, showed that AvrStb9 interacts with Stb9 resistance gene, triggering an immune response. AvrStb9 encodes an unusually large avirulence gene with a predicted secretion signal and a protease domain. It belongs to a S41 protease family conserved across different filamentous fungi in the Ascomycota class and may constitute a core effector. AvrStb9 is also conserved among a global Z. tritici population and carries multiple amino acid substitutions caused by strong positive diversifying selection. These results demonstrate the contribution of an 'atypical' conserved effector protein to fungal avirulence and the role of sequence diversification in the escape of host recognition, adding to our understanding of host-pathogen interactions and the evolutionary processes underlying pathogen adaptation.

Quantitative and qualitative plant-pathogen interactions call upon similar pathogenicity genes with a spectrum of effects

Septoria leaf blotch is a foliar wheat disease controlled by a combination of plant genetic resistances and fungicides use. R-gene-based qualitative resistance durability is limited due to gene-for-gene interactions with fungal avirulence (Avr) genes. Quantitative resistance is considered more durable but the mechanisms involved are not well documented. We hypothesize that genes involved in quantitative and qualitative plant-pathogen interactions are similar. A bi-parental population of Zymoseptoria tritici was inoculated on wheat cultivar 'Renan' and a linkage analysis performed to map QTL. Three pathogenicity QTL, Qzt-I05-1, Qzt-I05-6 and Qzt-I07-13, were mapped on chromosomes 1, 6 and 13 in Z. tritici, and a candidate pathogenicity gene on chromosome 6 was selected based on its effector-like characteristics. The candidate gene was cloned by Agrobacterium tumefaciens-mediated transformation, and a pathology test assessed the effect of the mutant strains on 'Renan'. This gene was demonstrated to be involved in quantitative pathogenicity. By cloning a newly annotated quantitative-effect gene in Z. tritici that is effector-like, we demonstrated that genes underlying pathogenicity QTL can be similar to Avr genes. This opens up the previously probed possibility that 'gene-for-gene' underlies not only qualitative but also quantitative plant-pathogen interactions in this pathosystem.

Natural variation in Avr3D1 from Zymoseptoria sp. contributes to quantitative gene-for-gene resistance and to host specificity

Successful host colonization by plant pathogens requires the circumvention of host defense responses, frequently through sequence modifications in secreted pathogen proteins known as avirulence factors (Avrs). Although Avr sequences are often polymorphic, the contribution of these polymorphisms to virulence diversity in natural pathogen populations remains largely unexplored. We used molecular genetic tools to determine how natural sequence polymorphisms of the avirulence factor Avr3D1 in the wheat pathogen Zymoseptoria tritici contributed to adaptive changes in virulence. We showed that there is a continuous distribution in the magnitude of resistance triggered by different Avr3D1 isoforms and demonstrated that natural variation in an Avr gene can lead to a quantitative resistance phenotype. We further showed that homologues of Avr3D1 in two nonpathogenic sister species of Z. tritici are recognized by some wheat cultivars, suggesting that Avr-R gene-for-gene interactions can contribute to nonhost resistance. We suggest that the mechanisms underlying host range, qualitative resistance, and quantitative resistance are not exclusive.

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

Effector Identification in Plant Pathogens

Effectors play a central role in determining the outcome of plant-pathogen interactions. As key virulence proteins, effectors are collectively indispensable for disease development. By understanding the virulence mechanisms of effectors, fundamental knowledge of microbial pathogenesis and disease resistance have been revealed. Effectors are also considered double-edged swords because some of them activate immunity in disease resistant plants after being recognized by specific immune receptors, which evolved to monitor pathogen presence or activity. Characterization of effector recognition by their cognate immune receptors and the downstream immune signaling pathways is instrumental in implementing resistance. Over the past decades, substantial research effort has focused on effector biology, especially concerning their interactions with virulence targets or immune receptors in plant cells. A foundation of this research is robust identification of the effector repertoire from a given pathogen, which depends heavily on bioinformatic prediction. In this review, we summarize methodologies that have been used for effector mining in various microbial pathogens which use different effector delivery mechanisms. We also discuss current limitations and provide perspectives on how recently developed analytic tools and technologies may facilitate effector identification and hence generation of a more complete vision of host-pathogen interactions. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

A highly multiplexed assay to monitor pathogenicity, fungicide resistance and gene flow in the fungal wheat pathogen Zymoseptoria tritici

Crop pathogens pose severe risks to global food production due to the rapid rise of resistance to pesticides and host resistance breakdowns. Predicting future risks requires monitoring tools to identify changes in the genetic composition of pathogen populations. Here we report the design of a microfluidics-based amplicon sequencing assay to multiplex 798 loci targeting virulence and fungicide resistance genes, and randomly selected genome-wide markers for the fungal pathogen Zymoseptoria tritici. The fungus causes one of the most devastating diseases on wheat showing rapid adaptation to fungicides and host resistance. We optimized the primer design by integrating polymorphism data from 632 genomes of the same species. To test the performance of the assay, we genotyped 192 samples in two replicates. Analysis of the short-read sequence data generated by the assay showed a fairly stable success rate across samples to amplify a large number of loci. The performance was consistent between samples originating from pure genomic DNA as well as material extracted directly from infected wheat leaves. In samples with mixed genotypes, we found that the assay recovers variations in allele frequencies. We also explored the potential of the amplicon assay to recover transposable element insertion polymorphism relevant for fungicide resistance. As a proof-of-concept, we show that the assay recovers the pathogen population structure across French wheat fields. Genomic monitoring of crop pathogens contributes to more sustainable crop protection and yields.

Combined pangenomics and transcriptomics reveals core and redundant virulence processes in a rapidly evolving fungal plant pathogen

BACKGROUND

Studying genomic variation in rapidly evolving pathogens potentially enables identification of genes supporting their "core biology", being present, functional and expressed by all strains or "flexible biology", varying between strains. Genes supporting flexible biology may be considered to be "accessory", whilst the "core" gene set is likely to be important for common features of a pathogen species biology, including virulence on all host genotypes. The wheat-pathogenic fungus Zymoseptoria tritici represents one of the most rapidly evolving threats to global food security and was the focus of this study.

RESULTS

We constructed a pangenome of 18 European field isolates, with 12 also subjected to RNAseq transcription profiling during infection. Combining this data, we predicted a "core" gene set comprising 9807 sequences which were (1) present in all isolates, (2) lacking inactivating polymorphisms and (3) expressed by all isolates. A large accessory genome, consisting of 45% of the total genes, was also defined. We classified genetic and genomic polymorphism at both chromosomal and individual gene scales. Proteins required for essential functions including virulence had lower-than average sequence variability amongst core genes. Both core and accessory genomes encoded many small, secreted candidate effector proteins that likely interact with plant immunity. Viral vector-mediated transient in planta overexpression of 88 candidates failed to identify any which induced leaf necrosis characteristic of disease. However, functional complementation of a non-pathogenic deletion mutant lacking five core genes demonstrated that full virulence was restored by re-introduction of the single gene exhibiting least sequence polymorphism and highest expression.

CONCLUSIONS

These data support the combined use of pangenomics and transcriptomics for defining genes which represent core, and potentially exploitable, weaknesses in rapidly evolving pathogens.

Genome sequencing and annotation of Cercospora sesami, a fungal pathogen causing leaf spot to Sesamum indicum

Cercospora sesami is a plant pathogen that causes leaf spot disease in sesame plants worldwide. In this study, genome sequence assembly of C. sesami isolate Cers 52-10 (MCC 9069) was generated using native paired-end and mate-pair DNA sequencing based on the Illumina HiSeq 2500 platform. The genome assembly of C. sesami is 34.3 Mb in size with an N50 of 26,222 bp and an average GC content of 53.02%. A total number of 10,872 genes were predicted in this study, out of which 9,712 genes were functionally annotated. Genes assigned to carbohydrate-active enzyme classes were also identified during the study. A total of 80 putative effector candidates were predicted and functionally annotated. The C. sesami genome sequence is available at DDBJ/ENA/GenBank, and other associated information is submitted to Mendeley's data.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03468-4.

Disentangling the complex gene interaction networks between rice and the blast fungus identifies a new pathogen effector

Studies focused solely on single organisms can fail to identify the networks underlying host-pathogen gene-for-gene interactions. Here, we integrate genetic analyses of rice (Oryza sativa, host) and rice blast fungus (Magnaporthe oryzae, pathogen) and uncover a new pathogen recognition specificity of the rice nucleotide-binding domain and leucine-rich repeat protein (NLR) immune receptor Pik, which mediates resistance to M. oryzae expressing the avirulence effector gene AVR-Pik. Rice Piks-1, encoded by an allele of Pik-1, recognizes a previously unidentified effector encoded by the M. oryzae avirulence gene AVR-Mgk1, which is found on a mini-chromosome. AVR-Mgk1 has no sequence similarity to known AVR-Pik effectors and is prone to deletion from the mini-chromosome mediated by repeated Inago2 retrotransposon sequences. AVR-Mgk1 is detected by Piks-1 and by other Pik-1 alleles known to recognize AVR-Pik effectors; recognition is mediated by AVR-Mgk1 binding to the integrated heavy metal-associated (HMA) domain of Piks-1 and other Pik-1 alleles. Our findings highlight how complex gene-for-gene interaction networks can be disentangled by applying forward genetics approaches simultaneously to the host and pathogen. We demonstrate dynamic coevolution between an NLR integrated domain and multiple families of effector proteins.

Whole genome sequencing of Clarireedia aff. paspali reveals potential pathogenesis factors in Clarireedia species, causal agents of dollar spot in turfgrass

Dollar spot is one of the most damaging diseases in turfgrass, reducing its quality and playability. Two species, Clarireedia monteithiana and C. jacksonii (formerly Sclerotinia homoeocarpa) have been reported so far in the United States To study the Clarireedia genome, two isolates H2 and H3, sampled from seashore paspalum in Hawaii in 2019 were sequenced via Illumina paired-end sequencing by synthesis technology and PacBio SMRT sequencing. Both isolates were identified as C. aff. paspali, a novel species in the United States Using short and long reads, C. aff. paspali H3 contained 193 contigs with 48.6 Mbp and presented the most completed assembly and annotation among Clarireedia species. Out of the 13,428 protein models from AUGUSTUS, 349 cytoplasmic effectors and 13 apoplastic effectors were identified by EffectorP. To further decipher Clarireedia pathogenicity, C. aff. paspali genomes (H2 and H3), as well as available C. jacksonii (LWC-10 and HRI11), C. monteithiana (DRR09 and RB-19) genomes were screened for fifty-four pathogenesis determinants, previously identified in S. sclerotiorum. Seventeen orthologs of pathogenicity genes have been identified in Clarireedia species involved in oxalic acid production (pac1, nox1), mitogen-activated protein kinase cascade (pka1, smk3, ste12), appressorium formation (caf1, pks13, ams2, rgb1, rhs1) and glycolytic pathway (gpd). Within these genes, 366 species-specific SNPs were recorded between Clarireedia species; twenty-eight were non-synonymous and non-conservative. The predicted protein structure of six of these genes showed superimposition of the models among Clarireedia spp. The genomic variations revealed here could potentially lead to differences in pathogenesis and other physiological functions among Clarireedia species.

New plant immunity elicitors from a sugar beet byproduct protect wheat against Zymoseptoria tritici

The current worldwide context promoting agroecology and green agriculture require the discovery of new ecofriendly and sustainable plant protection tools. Plant resistance inducers, called also elicitors, are one of the most promising alternatives fitting with such requirements. We produced here a set of 30 molecules from pyroglutamic acid, bio-sourced from sugar beet byproducts, and examined for their biological activity on the major agro-economically pathosystem wheat-Zymoseptoria tritici. Foliar application of the molecules provided significant protection rates (up to 63% disease severity reduction) for 16 among them. Structure-activity relationship analysis highlighted the importance of all chemical groups of the pharmacophore in the bioactivity of the molecules. Further investigations using in vitro and in planta antifungal bioassays as well as plant molecular biomarkers revealed that the activity of the molecules did not rely on direct biocide activity towards the pathogen, but rather on the activation of plant defense mechanisms dependent on lipoxygenase, phenylalanine ammonia-lyase, peroxidase, and pathogenesis-related protein pathways. This study reports a new family of bio-sourced resistance inducers and provides new insights into the valorization of agro-resources to develop the sustainable agriculture of tomorrow.

Prediction of effector protein structures from fungal phytopathogens enables evolutionary analyses

Elucidating the similarity and diversity of pathogen effectors is critical to understand their evolution across fungal phytopathogens. However, rapid divergence that diminishes sequence similarities between putatively homologous effectors has largely concealed the roots of effector evolution. Here we modelled the structures of 26,653 secreted proteins from 14 agriculturally important fungal phytopathogens, six non-pathogenic fungi and one oomycete with AlphaFold 2. With 18,000 successfully predicted folds, we performed structure-guided comparative analyses on two aspects of effector evolution: uniquely expanded sequence-unrelated structurally similar (SUSS) effector families and common folds present across the fungal species. Extreme expansion of lineage-specific SUSS effector families was found only in several obligate biotrophs, Blumeria graminis and Puccinia graminis. The highly expanded effector families were the source of conserved sequence motifs, such as the Y/F/WxC motif. We identified new classes of SUSS effector families that include known virulence factors, such as AvrSr35, AvrSr50 and Tin2. Structural comparisons revealed that the expanded structural folds further diversify through domain duplications and fusion with disordered stretches. Putatively sub- and neo-functionalized SUSS effectors could reconverge on regulation, expanding the functional pools of effectors in the pathogen infection cycle. We also found evidence that many effector families could have originated from ancestral folds conserved across fungi. Collectively, our study highlights diverse effector evolution mechanisms and supports divergent evolution as a major force in driving SUSS effector evolution from ancestral proteins.

Fungal plant pathogen "mutagenomics" reveals tagged and untagged mutations in Zymoseptoria tritici and identifies SSK2 as key morphogenesis and stress-responsive virulence factor

"Mutagenomics" is the combination of random mutagenesis, phenotypic screening, and whole-genome re-sequencing to uncover all tagged and untagged mutations linked with phenotypic changes in an organism. In this study, we performed a mutagenomics screen on the wheat pathogenic fungus Zymoseptoria tritici for altered morphogenetic switching and stress sensitivity phenotypes using Agrobacterium-mediated "random" T-DNA mutagenesis (ATMT). Biological screening identified four mutants which were strongly reduced in virulence on wheat. Whole genome re-sequencing defined the positions of the T-DNA insertion events and revealed several unlinked mutations potentially affecting gene functions. Remarkably, two independent reduced virulence mutant strains, with similarly altered stress sensitivities and aberrant hyphal growth phenotypes, were found to have a distinct loss of function mutations in the ZtSSK2 MAPKKK gene. One mutant strain had a direct T-DNA insertion affecting the predicted protein's N-terminus, while the other possessed an unlinked frameshift mutation towards the C-terminus. We used genetic complementation to restore both strains' wild-type (WT) function (virulence, morphogenesis, and stress response). We demonstrated that ZtSSK2 has a non-redundant function with ZtSTE11 in virulence through the biochemical activation of the stress-activated HOG1 MAPK pathway. Moreover, we present data suggesting that SSK2 has a unique role in activating this pathway in response to specific stresses. Finally, dual RNAseq-based transcriptome profiling of WT and SSK2 mutant strains revealed many HOG1-dependent transcriptional changes in the fungus during early infection and suggested that the host response does not discriminate between WT and mutant strains during this early phase. Together these data define new genes implicated in the virulence of the pathogen and emphasise the importance of a whole genome sequencing step in mutagenomic discovery pipelines.

Cytospora paraplurivora sp. nov. isolated from orchards with fruit tree decline syndrome in Ontario, Canada

A new species of Cytospora was isolated from cankered wood of Prunus spp. during a survey of orchards exhibiting symptoms of fruit tree decline syndrome in southern Ontario, Canada. We found isolates that are morphologically similar to species in the Cytosporaceae family, which is characterized by single or labyrinthine locules, filamentous conidiophores or clavate to elongate obovoid asci and allantoid, hyaline conidia. Multi-gene phylogenetic analysis of ITS, LSU, act and tef1- α showed that the isolates form a distinct clade, sister to Cytospora plurivora. Morphologically, our isolates showed differences in the length of conidia and culture characteristics compared to C. plurivora, suggesting the establishment of a new species. The species is described as Cytospora paraplurivora sp. nov. and placed in the family Cytosporaceae of Diaporthales. Additionally, we sequenced, assembled and characterized the genome of the representative isolate for this new species. The phylogenomic analysis confirms the species order and family level classification. C. paraplurivora sp. nov. has the potential to severely affect stone fruits production, causing cankers and dieback in stressed trees, and eventually leads to tree decline. Pathogenicity tests show that the species is pathogenic to Prunus persica var. persica.

Complete Genome Sequences of Septoria linicola: A Resource for Studying a Damaging Flax Pathogen

Fungal genus Septoria causes diseases in a wide range of plants. Here, we report the first genome sequences of two strains of Septoria linicola, the causal agent of the pasmo disease of flax (Linum usitatissimum). The genome of the first strain, SE15195, was fully assembled in 16 chromosomes, while 35 unitigs were obtained for a second strain, SE14017. Structural annotations predicted 13,096 and 13,085 protein-encoding genes and transposable elements content of 19.0 and 18.1% of the genome for SE15195 and SE14017, respectively. The four smaller chromosomes 13 to 16 show genomics features of potential accessory chromosomes. The assembly of these two genomes is a new resource for studying S. linicola and improving management of pasmo. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Heterochromatin in the fungal plant pathogen, Zymoseptoria tritici: Control of transposable elements, genome plasticity and virulence

Heterochromatin is a repressive chromatin state that plays key roles in the functional organisation of eukaryotic genomes. In fungal plant pathogens, effector genes that are required for host colonization tend to be associated with heterochromatic regions of the genome that are enriched with transposable elements. It has been proposed that the heterochromatin environment silences effector genes in the absence of host and dynamic chromatin remodelling facilitates their expression during infection. Here we discuss this model in the context of the key wheat pathogen, Zymoseptoria tritici. We cover progress in understanding the deposition and recognition of heterochromatic histone post translational modifications in Z. tritici and the role that heterochromatin plays in control of genome plasticity and virulence.

Distinct roles for different autophagy-associated genes in the virulence of the fungal wheat pathogen Zymoseptoria tritici

The fungal wheat pathogen Zymoseptoria tritici causes major crop losses as the causal agent of the disease Septoria tritici blotch. The infection cycle of Z. tritici displays two distinct phases, beginning with an extended symptomless phase of 1-2 weeks, before the fungus induces host cell death and tissue collapse in the leaf. Recent evidence suggests that the fungus uses little host-derived nutrition during asymptomatic colonisation, raising questions as to the sources of energy required for this initial growth phase. Autophagy is crucial for the pathogenicity of other fungal plant pathogens through its roles in supporting cellular differentiation and growth under starvation. Here we characterised the contributions of the autophagy genes ZtATG1 and ZtATG8 to the development and virulence of Z. tritici. Deletion of ZtATG1 led to inhibition of autophagy but had no impact on starvation-induced hyphal differentiation or virulence, suggesting that autophagy is not required for Z. tritici pathogenicity. Contrastingly, ZtATG8 deletion delayed the transition to necrotrophic growth, despite having no influence on filamentous growth under starvation, pointing to an autophagy-independent role of ZtATG8 during Z. tritici infection. To our knowledge, this study represents the first to find autophagy not to contribute to the virulence of a fungal plant pathogen, and reveals novel roles for different autophagy-associated proteins in Z. tritici.

Multi-site fungicides suppress banana Panama disease, caused by Fusarium oxysporum f. sp. cubense Tropical Race 4

Global banana production is currently challenged by Panama disease, caused by Fusarium oxysporum f.sp. cubense Tropical Race 4 (FocTR4). There are no effective fungicide-based strategies to control this soil-borne pathogen. This could be due to insensitivity of the pathogen to fungicides and/or soil application per se. Here, we test the effect of 12 single-site and 9 multi-site fungicides against FocTR4 and Foc Race1 (FocR1) in quantitative colony growth, and cell survival assays in purified FocTR4 macroconidia, microconidia and chlamydospores. We demonstrate that these FocTR4 morphotypes all cause Panama disease in bananas. These experiments reveal innate resistance of FocTR4 to all single-site fungicides, with neither azoles, nor succinate dehydrogenase inhibitors (SDHIs), strobilurins or benzimidazoles killing these spore forms. We show in fungicide-treated hyphae that this innate resistance occurs in a subpopulation of "persister" cells and is not genetically inherited. FocTR4 persisters respond to 3 μg ml-1 azoles or 1000 μg ml-1 strobilurins or SDHIs by strong up-regulation of genes encoding target enzymes (up to 660-fold), genes for putative efflux pumps and transporters (up to 230-fold) and xenobiotic detoxification enzymes (up to 200-fold). Comparison of gene expression in FocTR4 and Zymoseptoria tritici, grown under identical conditions, reveals that this response is only observed in FocTR4. In contrast, FocTR4 shows little innate resistance to most multi-site fungicides. However, quantitative virulence assays, in soil-grown bananas, reveals that only captan (20 μg ml-1) and all lipophilic cations (200 μg ml-1) suppress Panama disease effectively. These fungicides could help protect bananas from future yield losses by FocTR4.

Genome-wide association mapping reveals genes underlying population-level metabolome diversity in a fungal crop pathogen

Background

Fungi produce a wide range of specialized metabolites (SMs) involved in biotic interactions. Pathways for the production of SMs are often encoded in clusters of tightly arranged genes identified as biosynthetic gene clusters. Such gene clusters can undergo horizontal gene transfers between species and rapid evolutionary change within species. The acquisition, rearrangement, and deletion of gene clusters can generate significant metabolome diversity. However, the genetic basis underlying variation in SM production remains poorly understood.

Results

Here, we analyzed the metabolite production of a large population of the fungal pathogen of wheat, Zymoseptoria tritici. The pathogen causes major yield losses and shows variation in gene clusters. We performed untargeted ultra-high performance liquid chromatography-high resolution mass spectrometry to profile the metabolite diversity among 102 isolates of the same species. We found substantial variation in the abundance of the detected metabolites among isolates. Integrating whole-genome sequencing data, we performed metabolite genome-wide association mapping to identify loci underlying variation in metabolite production (i.e., metabolite-GWAS). We found that significantly associated SNPs reside mostly in coding and gene regulatory regions. Associated genes encode mainly transport and catalytic activities. The metabolite-GWAS identified also a polymorphism in the 3′UTR region of a virulence gene related to metabolite production and showing expression variation.

Conclusions

Taken together, our study provides a significant resource to unravel polymorphism underlying metabolome diversity within a species. Integrating metabolome screens should be feasible for a range of different plant pathogens and help prioritize molecular studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01422-z.

Multi-Layer Biosensor for Pre-Symptomatic Detection of Puccinia strifformis, the Causal Agent of Yellow Rust

The yellow rust of wheat (caused by Puccinia striiformis f. sp. tritici) is a devastating fungal infection that is responsible for significant wheat yield losses. The main challenge with the detection of this disease is that it can only be visually detected on the leaf surface between 7 and 10 days after infection, and by this point, counter measures such as the use of fungicides are generally less effective. The hypothesis of this study is to develop and use a compact electrochemical-based biosensor for the early detection of P. striiformis, thus enabling fast countermeasures to be taken. The biosensor that was developed consists of three layers. The first layer mimics the wheat leaf surface morphology. The second layer consists of a sucrose/agar mixture that acts as a substrate and contains a wheat-derived terpene volatile organic compound that stimulates the germination and growth of the spores of the yellow rust pathogen P. s. f. sp. tritici. The third layer consists of a nonenzymatic glucose sensor that produces a signal once invertase is produced by P. striiformis, which comes into contact with the second layer, thereby converting sucrose to glucose. The results show the proof that this innovative biosensor can enable the detection of yellow rust spores in 72 h.

Zymoseptoria tritici white-collar complex integrates light, temperature and plant cues to initiate dimorphism and pathogenesis

Transitioning from spores to hyphae is pivotal to host invasion by the plant pathogenic fungus Zymoseptoria tritici. This dimorphic switch can be initiated by high temperature in vitro (~27 °C); however, such a condition may induce cellular heat stress, questioning its relevance to field infections. Here, we study the regulation of the dimorphic switch by temperature and other factors. Climate data from wheat-growing areas indicate that the pathogen sporadically experiences high temperatures such as 27 °C during summer months. However, using a fluorescent dimorphic switch reporter (FDR1) in four wild-type strains, we show that dimorphic switching already initiates at 15-18 °C, and is enhanced by wheat leaf surface compounds. Transcriptomics reveals 1261 genes that are up- or down-regulated in hyphae of all strains. These pan-strain core dimorphism genes (PCDGs) encode known effectors, dimorphism and transcription factors, and light-responsive proteins (velvet factors, opsins, putative blue light receptors). An FDR1-based genetic screen reveals a crucial role for the white-collar complex (WCC) in dimorphism and virulence, mediated by control of PCDG expression. Thus, WCC integrates light with biotic and abiotic cues to orchestrate Z. tritici infection.

Inactivation of a candidate effector gene of Zymoseptoria tritici affects its sporulation

BACKGROUND

Wheat is one of the most important staple crops produced worldwide. Its susceptibility to plant diseases reduces its production significantly. One of the most important diseases of wheat is septoria tritici blotch, a devastating disease observed in fields with wet and temperate conditions. Z. tritici secretes effector proteins to influence the host's defense mechanisms, as is typical of plant pathogens. In this investigation, we evaluated the pathogenicity of some Zymoseptoria tritici effector candidate genes having a signal peptide for secretion with no known function.

METHODS AND RESULTS

Three genes named Mycgr3G104383, Mycgr3G104444 and Mycgr3G105826 were knocked out separately through homologous recombination, generating Z. tritici IPO323 mutants lacking the functional copy of the corresponding genes. While KO1 and KO3 mutants did not show any significant differences during phenotypic and virulence investigations, the KO2 mutant generated exclusively macropycnidiospores in artificial media, different from wild-type IPO323 which produce only micropycidiospores. The mycelial growth capability of KO2 was also severely attenuated in all of the investigated growth conditions. These changes were observed independent of growth media and growth temperatures, implying that changes were genetic and inherited through generations. Virulence of knockout mutants in wheat leaves was observed to be similar to the wild-type IPO323.

CONCLUSION

Understanding the biology of Z. tritici and its interactions with wheat will reveal new strategies to fight septoria tritici blotch, enabling breeding wheat cultivars resistant to a broader spectrum of Z. tritici strains. Furthermore, gene knockout via homologous recombination proved to be a powerful tool for discovering novel gene functions.

Current Status and Future Perspectives of Genomics Research in the Rust Fungi

Rust fungi in Pucciniales have caused destructive plant epidemics, have become more aggressive with new virulence, rapidly adapt to new environments, and continually threaten global agriculture. With the rapid advancement of genome sequencing technologies and data analysis tools, genomics research on many of the devastating rust fungi has generated unprecedented insights into various aspects of rust biology. In this review, we first present a summary of the main findings in the genomics of rust fungi related to variations in genome size and gene composition between and within species. Then we show how the genomics of rust fungi has promoted our understanding of the pathogen virulence and population dynamics. Even with great progress, many questions still need to be answered. Therefore, we introduce important perspectives with emphasis on the genome evolution and host adaptation of rust fungi. We believe that the comparative genomics and population genomics of rust fungi will provide a further understanding of the rapid evolution of virulence and will contribute to monitoring the population dynamics for disease management.

Characterization of Host-Specific Genes from Pine- and Grass-Associated Species of the Fusarium fujikuroi Species Complex

The Fusarium fujikuroi species complex (FFSC) includes socioeconomically important pathogens that cause disease for numerous crops and synthesize a variety of secondary metabolites that can contaminate feedstocks and food. Here, we used comparative genomics to elucidate processes underlying the ability of pine-associated and grass-associated FFSC species to colonize tissues of their respective plant hosts. We characterized the identity, possible functions, evolutionary origins, and chromosomal positions of the host-range-associated genes encoded by the two groups of fungi. The 72 and 47 genes identified as unique to the respective genome groups were potentially involved in diverse processes, ranging from transcription, regulation, and substrate transport through to virulence/pathogenicity. Most genes arose early during the evolution of Fusarium/FFSC and were only subsequently retained in some lineages, while some had origins outside Fusarium. Although differences in the densities of these genes were especially noticeable on the conditionally dispensable chromosome of F. temperatum (representing the grass-associates) and F. circinatum (representing the pine-associates), the host-range-associated genes tended to be located towards the subtelomeric regions of chromosomes. Taken together, these results demonstrate that multiple mechanisms drive the emergence of genes in the grass- and pine-associated FFSC taxa examined. It also highlighted the diversity of the molecular processes potentially underlying niche-specificity in these and other Fusarium species.

Non-Mendelian transmission of accessory chromosomes in fungi

Non-Mendelian transmission has been reported for various genetic elements, ranging from small transposons to entire chromosomes. One prime example of such a transmission pattern are B chromosomes in plants and animals. Accessory chromosomes in fungi are similar to B chromosomes in showing presence/absence polymorphism and being non-essential. How these chromosomes are transmitted during meiosis is however poorly understood—despite their often high impact on the fitness of the host. For several fungal organisms, a non-Mendelian transmission or a mechanistically unique meiotic drive of accessory chromosomes have been reported. In this review, we provide an overview of the possible mechanisms that can cause the non-Mendelian transmission or meiotic drives of fungal accessory chromosomes. We compare processes responsible for the non-Mendelian transmission of accessory chromosomes for different fungal eukaryotes and discuss the structural traits of fungal accessory chromosomes affecting their meiotic transmission. We conclude that research on fungal accessory chromosomes, due to their small size, ease of sequencing, and epigenetic profiling, can complement the study of B chromosomes in deciphering factors that influence and regulate the non-Mendelian transmission of entire chromosomes.