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

Physiological and molecular responses of a resistant and susceptible wheat cultivar to the fungal wheat pathogen Zymoseptoria tritici

Zymoseptoria tritici is the causal agent of Septoria tritici blotch (STB), one of the most economically destructive wheat foliar diseases. In this study, we explore the physiological and molecular changes elicited in two wheat cultivars with divergent responses (Taichung 29 = susceptible, and Shafir = resistant) upon infection by Z. tritici. Our aim is to uncover novel insights into the intricate mechanisms that govern wheat defense against Z. tritici infection. Our quantitative histopathological study showed that H2O2 accumulated in the resistant cultivar to a higher degree compared to the susceptible cultivar at the biotrophic and switching phase. Additionally, we combined qPCR with a targeted quantitative HPLC technique to evaluate the expression profiles of 13 defense-related genes and profile the polyphenolic compounds induced differentially in the STB susceptible and resistant cultivar. Our finding indicated that five out of 13 genes were strongly up-regulated in the resistant cultivar compared with that of the susceptible one at eight days post-inoculation (dpi), corresponding to the transition phase present in the infection process of Z. tritici. Finally, our targeted HPLC analysis demonstrated that the traced phenolic compounds were highly elevated in the susceptible cultivar infected by Z. tritici compared with that of the resistant cultivar. In conclusion, our comprehensive analysis unveils a robust defense response in the resistant wheat cultivar Shafir, characterized by heightened H2O2 accumulation, significant up-regulation of key defense-related genes during the transition phase, and a distinct profile of polyphenolic compounds, shedding light on the intricate mechanisms contributing to its resistance against Z. tritici, thereby providing valuable insights for the development of more resilient wheat varieties.

An array of Zymoseptoria tritici effectors suppress plant immune responses

Zymoseptoria tritici is the most economically significant fungal pathogen of wheat in Europe. However, despite the importance of this pathogen, the molecular interactions between pathogen and host during infection are not well understood. Herein, we describe the use of two libraries of cloned Z. tritici effectors that were screened to identify effector candidates with putative pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI)-suppressing activity. The effectors from each library were transiently expressed in Nicotiana benthamiana, and expressing leaves were treated with bacterial or fungal PAMPs to assess the effectors' ability to suppress reactive oxygen species (ROS) production. From these screens, numerous effectors were identified with PTI-suppressing activity. In addition, some effectors were able to suppress cell death responses induced by other Z. tritici secreted proteins. We used structural prediction tools to predict the putative structures of all of the Z. tritici effectors and used these predictions to examine whether there was enrichment of specific structural signatures among the PTI-suppressing effectors. From among the libraries, multiple members of the killer protein-like 4 (KP4) and killer protein-like 6 (KP6) effector families were identified as PTI suppressors. This observation is intriguing, as these protein families were previously associated with antimicrobial activity rather than virulence or host manipulation. This data provides mechanistic insight into immune suppression by Z. tritici during infection and suggests that, similar to biotrophic pathogens, this fungus relies on a battery of secreted effectors to suppress host immunity during early phases of colonization.

Functional characterization of extracellular and intracellular catalase-peroxidases involved in virulence of the fungal wheat pathogen Zymoseptoria tritici

Understanding how pathogens defend themselves against host defence mechanisms, such as hydrogen peroxide (H2O2) production, is crucial for comprehending fungal infections. H2O2 poses a significant threat to invading fungi due to its potent oxidizing properties. Our research focuses on the hemibiotrophic fungal wheat pathogen Zymoseptoria tritici, enabling us to investigate host-pathogen interactions. We examined two catalase-peroxidase (CP) genes, ZtCpx1 and ZtCpx2, to elucidate how Z. tritici deals with host-generated H2O2 during infection. Our analysis revealed that ZtCpx1 was up-regulated during biotrophic growth and asexual spore formation in vitro, while ZtCpx2 showed increased expression during the transition from biotrophic to necrotrophic growth and in-vitro vegetative growth. Deleting ZtCpx1 increased the mutant's sensitivity to exogenously added H2O2 and significantly reduced virulence, as evidenced by decreased Septoria tritici blotch symptom severity and fungal biomass production. Reintroducing the wild-type ZtCpx1 allele with its native promoter into the mutant strain restored the observed phenotypes. While ZtCpx2 was not essential for full virulence, the ZtCpx2 mutants exhibited reduced fungal biomass development during the transition from biotrophic to necrotrophic growth. Moreover, both CP genes act synergistically, as the double knock-out mutant displayed a more pronounced reduced virulence compared to ΔZtCpx1. Microscopic analysis using fluorescent proteins revealed that ZtCpx1 was localized in the peroxisome, indicating its potential role in managing host-generated reactive oxygen species during infection. In conclusion, our research sheds light on the crucial roles of CP genes ZtCpx1 and ZtCpx2 in the defence mechanism of Z. tritici against host-generated hydrogen peroxide.

Highly Repetitive Genome of Coniella granati (syn. Pilidiella granati), the Causal Agent of Pomegranate Fruit Rot, Encodes a Minimalistic Proteome with a Streamlined Arsenal of Effector Proteins

This study describes the first genome sequence and analysis of Coniella granati, a fungal pathogen with a broad host range, which is responsible for postharvest crown rot, shoot blight, and canker diseases in pomegranates. C. granati is a geographically widespread pathogen which has been reported across Europe, Asia, the Americas, and Africa. Our analysis revealed a 46.8 Mb genome with features characteristic of hemibiotrophic fungi. Approximately one third of its genome was compartmentalised within 'AT-rich' regions exhibiting a low GC content (30 to 45%). These regions primarily comprised transposable elements that are repeated at a high frequency and interspersed throughout the genome. Transcriptome-supported gene annotation of the C. granati genome revealed a streamlined proteome, mirroring similar observations in other pathogens with a latent phase. The genome encoded a relatively compact set of 9568 protein-coding genes with a remarkable 95% having assigned functional annotations. Despite this streamlined nature, a set of 40 cysteine-rich candidate secreted effector-like proteins (CSEPs) was predicted as well as a gene cluster involved in the synthesis of a pomegranate-associated toxin. These potential virulence factors were predominantly located near repeat-rich and AT-rich regions, suggesting that the pathogen evades host defences through Repeat-Induced Point mutation (RIP)-mediated pseudogenisation. Furthermore, 23 of these CSEPs exhibited homology to known effector and pathogenicity genes found in other hemibiotrophic pathogens. The study establishes a foundational resource for the study of the genetic makeup of C. granati, paving the way for future research on its pathogenicity mechanisms and the development of targeted control strategies to safeguard pomegranate production.

Comparative genomics of the extremophile Cryomyces antarcticus and other psychrophilic Dothideomycetes

Over a billion years of fungal evolution has enabled representatives of this kingdom to populate almost all parts of planet Earth and to adapt to some of its most uninhabitable environments including extremes of temperature, salinity, pH, water, light, or other sources of radiation. Cryomyces antarcticus is an endolithic fungus that inhabits rock outcrops in Antarctica. It survives extremes of cold, humidity and solar radiation in one of the least habitable environments on Earth. This fungus is unusual because it produces heavily melanized, meristematic growth and is thought to be haploid and asexual. Due to its growth in the most extreme environment, it has been suggested as an organism that could survive on Mars. However, the mechanisms it uses to achieve its extremophilic nature are not known. Comparative genomics can provide clues to the processes underlying biological diversity, evolution, and adaptation. This effort has been greatly facilitated by the 1000 Fungal Genomes project and the JGI MycoCosm portal where sequenced genomes have been assembled into phylogenetic and ecological groups representing different projects, lifestyles, ecologies, and evolutionary histories. Comparative genomics within and between these groups provides insights into fungal adaptations, for example to extreme environmental conditions. Here, we analyze two Cryomyces genomes in the context of additional psychrophilic fungi, as well as non-psychrophilic fungi with diverse lifestyles selected from the MycoCosm database. This analysis identifies families of genes that are expanded and contracted in Cryomyces and other psychrophiles and may explain their extremophilic lifestyle. Higher GC contents of genes and of bases in the third positions of codons may help to stabilize DNA under extreme conditions. Numerous smaller contigs in C. antarcticus suggest the presence of an alternative haplotype that could indicate the sequenced isolate is diploid or dikaryotic. These analyses provide a first step to unraveling the secrets of the extreme lifestyle of C. antarcticus.

Life on a leaf: the epiphyte to pathogen continuum and interplay in the phyllosphere

Epiphytic microbes are those that live for some or all of their life cycle on the surface of plant leaves. Leaf surfaces are a topologically complex, physicochemically heterogeneous habitat that is home to extensive, mixed communities of resident and transient inhabitants from all three domains of life. In this review, we discuss the origins of leaf surface microbes and how different biotic and abiotic factors shape their communities. We discuss the leaf surface as a habitat and microbial adaptations which allow some species to thrive there, with particular emphasis on microbes that occupy the continuum between epiphytic specialists and phytopathogens, groups which have considerable overlap in terms of adapting to the leaf surface and between which a single virulence determinant can move a microbial strain. Finally, we discuss the recent findings that the wheat pathogenic fungus Zymoseptoria tritici spends a considerable amount of time on the leaf surface, and ask what insights other epiphytic organisms might provide into this pathogen, as well as how Z. tritici might serve as a model system for investigating plant-microbe-microbe interactions on the leaf surface.

Stomatal penetration: the cornerstone of plant resistance to the fungal pathogen Zymoseptoria tritici

BACKGROUND

Septoria tritici blotch (STB), caused by the foliar fungus Zymoseptoria tritici, is one of the most damaging disease of wheat in Europe. Genetic resistance against this fungus relies on different types of resistance from non-host resistance (NHR) and host species specific resistance (HSSR) to host resistance mediated by quantitative trait loci (QTLs) or major resistance genes (Stb). Characterizing the diversity of theses resistances is of great importance for breeding wheat cultivars with efficient and durable resistance. While the functional mechanisms underlying these resistance types are not well understood, increasing piece of evidence suggest that fungus stomatal penetration and early establishment in the apoplast are both crucial for the outcome of some interactions between Z. tritici and plants. To validate and extend these previous observations, we conducted quantitative comparative phenotypical and cytological analyses of the infection process corresponding to 22 different interactions between plant species and Z. tritici isolates. These interactions included four major bread wheat Stb genes, four bread wheat accessions with contrasting quantitative resistance, two species resistant to Z. tritici isolates from bread wheat (HSSR) and four plant species resistant to all Z. tritici isolates (NHR).

RESULTS

Infiltration of Z. tritici spores into plant leaves allowed the partial bypass of all bread wheat resistances and durum wheat resistance, but not resistances from other plants species. Quantitative comparative cytological analysis showed that in the non-grass plant Nicotiana benthamiana, Z. tritici was stopped before stomatal penetration. By contrast, in all resistant grass plants, Z. tritici was stopped, at least partly, during stomatal penetration. The intensity of this early plant control process varied depending on resistance types, quantitative resistances being the least effective. These analyses also demonstrated that Stb-mediated resistances, HSSR and NHR, but not quantitative resistances, relied on the strong growth inhibition of the few Z. tritici penetrating hyphae at their entry point in the sub-stomatal cavity.

CONCLUSIONS

In addition to furnishing a robust quantitative cytological assessment system, our study uncovered three stopping patterns of Z. tritici by plant resistances. Stomatal resistance was found important for most resistances to Z. tritici, independently of its type (Stb, HSSR, NHR). These results provided a basis for the functional analysis of wheat resistance to Z. tritici and its improvement.

What are the 100 most cited fungal genera?

The global diversity of fungi has been estimated between 2 to 11 million species, of which only about 155 000 have been named. Most fungi are invisible to the unaided eye, but they represent a major component of biodiversity on our planet, and play essential ecological roles, supporting life as we know it. Although approximately 20 000 fungal genera are presently recognised, the ecology of most remains undetermined. Despite all this diversity, the mycological community actively researches some fungal genera more commonly than others. This poses an interesting question: why have some fungal genera impacted mycology and related fields more than others? To address this issue, we conducted a bibliometric analysis to identify the top 100 most cited fungal genera. A thorough database search of the Web of Science, Google Scholar, and PubMed was performed to establish which genera are most cited. The most cited 10 genera are Saccharomyces, Candida, Aspergillus, Fusarium, Penicillium, Trichoderma, Botrytis, Pichia, Cryptococcus and Alternaria. Case studies are presented for the 100 most cited genera with general background, notes on their ecology and economic significance and important research advances. This paper provides a historic overview of scientific research of these genera and the prospect for further research. Citation: Bhunjun CS, Chen YJ, Phukhamsakda C, Boekhout T, Groenewald JZ, McKenzie EHC, Francisco EC, Frisvad JC, Groenewald M, Hurdeal VG, Luangsa-ard J, Perrone G, Visagie CM, Bai FY, Błaszkowski J, Braun U, de Souza FA, de Queiroz MB, Dutta AK, Gonkhom D, Goto BT, Guarnaccia V, Hagen F, Houbraken J, Lachance MA, Li JJ, Luo KY, Magurno F, Mongkolsamrit S, Robert V, Roy N, Tibpromma S, Wanasinghe DN, Wang DQ, Wei DP, Zhao CL, Aiphuk W, Ajayi-Oyetunde O, Arantes TD, Araujo JC, Begerow D, Bakhshi M, Barbosa RN, Behrens FH, Bensch K, Bezerra JDP, Bilański P, Bradley CA, Bubner B, Burgess TI, Buyck B, Čadež N, Cai L, Calaça FJS, Campbell LJ, Chaverri P, Chen YY, Chethana KWT, Coetzee B, Costa MM, Chen Q, Custódio FA, Dai YC, Damm U, de Azevedo Santiago ALCM, De Miccolis Angelini RM, Dijksterhuis J, Dissanayake AJ, Doilom M, Dong W, Alvarez-Duarte E, Fischer M, Gajanayake AJ, Gené J, Gomdola D, Gomes AAM, Hausner G, He MQ, Hou L, Iturrieta-González I, Jami F, Jankowiak R, Jayawardena RS, Kandemir H, Kiss L, Kobmoo N, Kowalski T, Landi L, Lin CG, Liu JK, Liu XB, Loizides M, Luangharn T, Maharachchikumbura SSN, Makhathini Mkhwanazi GJ, Manawasinghe IS, Marin-Felix Y, McTaggart AR, Moreau PA, Morozova OV, Mostert L, Osiewacz HD, Pem D, Phookamsak R, Pollastro S, Pordel A, Poyntner C, Phillips AJL, Phonemany M, Promputtha I, Rathnayaka AR, Rodrigues AM, Romanazzi G, Rothmann L, Salgado-Salazar C, Sandoval-Denis M, Saupe SJ, Scholler M, Scott P, Shivas RG, Silar P, Souza-Motta CM, Silva-Filho AGS, Spies CFJ, Stchigel AM, Sterflinger K, Summerbell RC, Svetasheva TY, Takamatsu S, Theelen B, Theodoro RC, Thines M, Thongklang N, Torres R, Turchetti B, van den Brule T, Wang XW, Wartchow F, Welti S, Wijesinghe SN, Wu F, Xu R, Yang ZL, Yilmaz N, Yurkov A, Zhao L, Zhao RL, Zhou N, Hyde KD, Crous PW (2024). What are the 100 most cited fungal genera? Studies in Mycology 108: 1-411. doi: 10.3114/sim.2024.108.01.

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.

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.

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.

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.