Transcriptome and metabolite profiling of the infection cycle of Zymoseptoria tritici on wheat reveals a biphasic interaction with plant immunity involving differential pathogen chromosomal contributions and a variation on the hemibiotrophic lifestyle definition

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

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.

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.

Infection-specific transcriptional patterns of the maize pathogen Cochliobolus heterostrophus unravel genes involved in asexual development and virulence

Southern corn leaf blight (SCLB) caused by Cochliobolus heterostrophus is a destructive disease that threatens global maize (Zea mays) production. Despite many studies being conducted, very little is known about molecular processes employed by the pathogen during infection. There is a need to understand the fungal arms strategy and identify novel functional genes as targets for fungicide development. Transcriptome analysis based on RNA sequencing was carried out across conidia germination and host infection by C. heterostrophus. The present study revealed major changes in C. heterostrophus gene expression during host infection. Several differentially expressed genes (DEGs) induced during C. heterostrophus infection could be involved in the biosynthesis of secondary metabolites, peroxisome, energy metabolism, amino acid degradation and oxidative phosphorylation. In addition, histone acetyltransferase, secreted proteins, peroxisomal proteins, NADPH oxidase and transcription factors were selected for further functional validation. Here, we demonstrated that histone acetyltransferases (Hat2 and Rtt109), secreted proteins (Cel61A and Mep1), peroxisomal proteins (Pex11A and Pex14), NADPH oxidases (NoxA, NoxD and NoxR) and transcription factors (Crz1 and MtfA) play essential roles in C. heterostrophus conidiation, stress adaption and virulence. Taken together, our study revealed major changes in gene expression associated with C. heterostrophus infection and identified a diverse repertoire of genes critical for successful infection.

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.

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.

Comparison of the impact of two key fungal signalling pathways on Zymoseptoria tritici infection reveals divergent contribution to invasive growth through distinct regulation of infection-associated genes

The lifecycle of Zymoseptoria tritici requires a carefully regulated asymptomatic phase within the wheat leaf following penetration of the mesophyll via stomata. Here we compare the roles in this process of two key fungal signalling pathways, mutants of which were identified through forward genetics due to their avirulence on wheat. Whole-genome resequencing of avirulent Z. tritici T-DNA transformants identified disruptive mutations in ZtBCK1 from the kinase cascade of the cell wall integrity (CWI) pathway, and the adenylate cyclase gene ZtCYR1. Targeted deletion of these genes abolished the pathogenicity of the fungus and led to similar in vitro phenotypes to those associated with disruption of putative downstream kinases, both supporting previous studies and confirming the importance of these pathways in virulence. RNA sequencing was used to investigate the effect of ZtBCK1 and ZtCYR1 deletion on gene expression in both the pathogen and host during infection. ZtBCK1 was found to be required for the adaptation to the host environment, controlling expression of infection-associated secreted proteins, including known virulence factors. Meanwhile, ZtCYR1 is implicated in controlling the switch to necrotrophy, regulating expression of effectors associated with this transition. This represents the first study to compare the influence of CWI and cAMP signalling on in planta transcription of a fungal plant pathogen, providing insights into their differential regulation of candidate effectors during invasive growth.

Seeking the interspecies crosswalk for filamentous microbe effectors

Both pathogenic and symbiotic microorganisms modulate the immune response and physiology of their host to establish a suitable niche. Key players in mediating colonization outcome are microbial effector proteins that act either inside (cytoplasmic) or outside (apoplastic) the plant cells and modify the abundance or activity of host macromolecules. We compile novel insights into the much-disputed processes of effector secretion and translocation of filamentous organisms, namely fungi and oomycetes. We report how recent studies that focus on unconventional secretion and effector structure challenge the long-standing image of effectors as conventionally secreted proteins that are translocated with the aid of primary amino acid sequence motifs. Furthermore, we emphasize the potential of diverse, unbiased, state-of-the-art proteomics approaches in the holistic characterization of fungal and oomycete effectomes.

Biology, taxonomy, genetics, and management of Zymoseptoria tritici: the causal agent of wheat leaf blotch

Septoria tritici blotch or Septoria leaf blotch has been used for long time, but leaf blotch is a correct disease name. Moreover, Lb resistant gene is the correct name, but, not Stb gene. It has sexual and asexual parts on the mycelia, known as heterothallic fungi. Its pathogenic diversity ranged from 40% to 93% and has produced a wide variety of AvrLb6 haplotypes. M. graminicola has a plasmogamy and karyogamy sexual process. The pathogen can use macropycnidiospores, micropycnidiospores, and pycnidia vegetative growths for infection and overwintering. Synthetic M3, Kavkaz-K4500, Synthetic 6×, and TE9111 wheat genotypes have horizontal resistance. Avirulence (Avr) genes in Z. tritici and their matching wheat (R) genes indicate gene for gene mechanisms of resistance. Twenty-two R genes (vertical resistance) have been identified. In both horizontal and vertical resistance, different Lb genes have been broken down due to new Z.tritici virulent gene and currently Lb19 resistant gene is being recommended. Mixing of resistant and susceptible cultivars is also the most effective management strategy. Moreover, different cultural practices and biological control have been proposed. Lastly, different fungicides are also available. However, in developing countries cultivar mixture, isolates diversity, biological control, and epidemic studies have been greatly missed.

Polygenic pathogen networks influence transcriptional plasticity in the Arabidopsis-Botrytis pathosystem

Bidirectional flow of information shapes the outcome of the host-pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens' genetic variation and both organisms' transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host's genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.

Defence-related metabolic changes in wheat (Triticum aestivum L.) seedlings in response to infection by Puccinia graminis f. sp. tritici

Stem rust caused by the pathogen Puccinia graminis f. sp. tritici is a destructive fungal disease-causing major grain yield losses in wheat. Therefore, understanding the plant defence regulation and function in response to the pathogen attack is required. As such, an untargeted LC-MS-based metabolomics approach was employed as a tool to dissect and understand the biochemical responses of Koonap (resistant) and Morocco (susceptible) wheat varieties infected with two different races of P. graminis (2SA88 [TTKSF] and 2SA107 [PTKST]). Data was generated from the infected and non-infected control plants harvested at 14- and 21- days post-inoculation (dpi), with 3 biological replicates per sample under a controlled environment. Chemo-metric tools such as principal component analysis (PCA), orthogonal projection to latent structures-discriminant analysis (OPLS-DA) were used to highlight the metabolic changes using LC-MS data of the methanolic extracts generated from the two wheat varieties. Molecular networking in Global Natural Product Social (GNPS) was further used to analyse biological networks between the perturbed metabolites. PCA and OPLS-DA analysis showed cluster separations between the varieties, infection races and the time-points. Distinct biochemical changes were also observed between the races and time-points. Metabolites were identified and classified using base peak intensities (BPI) and single ion extracted chromatograms from samples, and the most affected metabolites included flavonoids, carboxylic acids and alkaloids. Network analysis also showed high expression of metabolites from thiamine and glyoxylate, such as flavonoid glycosides, suggesting multi-faceted defence response strategy by understudied wheat varieties towards P. graminis pathogen infection. Overall, the study provided the insights of the biochemical changes in the expression of wheat metabolites in response to stem rust infection.

Genome-wide characterization of the VQ genes in Triticeae and their functionalization driven by polyploidization and gene duplication events in wheat

Valine-glutamine motif-containing (VQ) proteins are transcriptional cofactors widely involved in plant growth, development, and response to various stresses. Although the VQ family has been genome-wide identified in some species, but the knowledge regarding duplication-driven functionalization of VQ genes among evolutionarily related species is still lacking. Here, 952 VQ genes have been identified from 16 species, emphasizing seven Triticeae species including the bread wheat. Comprehensive phylogenetic and syntenic analyses allow us to establish the orthologous relationship of VQ genes from rice (Oryza sativa) to bread wheat (Triticum aestivum). The evolutionary analysis revealed that whole-genome duplication (WGD) drives the expansion of OsVQs, while TaVQs expansion is associated with a recent burst of gene duplication (RBGD). We also analyzed the motif composition and molecular properties of TaVQ proteins, enriched biological functions, and expression patterns of TaVQs. We demonstrate that WGD-derived TaVQs have become divergent in both protein motif composition and expression pattern, while RBGD-derived TaVQs tend to adopt specific expression patterns, suggesting their functionalization in certain biological processes or in response to specific stresses. Furthermore, some RBGD-derived TaVQs are found to be associated with salt tolerance. Several of the identified salt-related TaVQ proteins were located in the cytoplasm and nucleus and their salt-responsive expression patterns were validated by qPCR analysis. Yeast-based functional experiments confirmed that TaVQ27 may be a new regulator to salt response and regulation. Overall, this study lays the foundation for further functional validation of VQ family members within the Triticeae species.

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.

The transcriptional landscape of plant infection by the rice blast fungus Magnaporthe oryzae reveals distinct families of temporally co-regulated and structurally conserved effectors

The rice blast fungus Magnaporthe oryzae causes a devastating disease that threatens global rice (Oryza sativa) production. Despite intense study, the biology of plant tissue invasion during blast disease remains poorly understood. Here we report a high-resolution transcriptional profiling study of the entire plant-associated development of the blast fungus. Our analysis revealed major temporal changes in fungal gene expression during plant infection. Pathogen gene expression could be classified into 10 modules of temporally co-expressed genes, providing evidence for the induction of pronounced shifts in primary and secondary metabolism, cell signaling, and transcriptional regulation. A set of 863 genes encoding secreted proteins are differentially expressed at specific stages of infection, and 546 genes named MEP (Magnaportheeffector protein) genes were predicted to encode effectors. Computational prediction of structurally related MEPs, including the MAX effector family, revealed their temporal co-regulation in the same co-expression modules. We characterized 32 MEP genes and demonstrate that Mep effectors are predominantly targeted to the cytoplasm of rice cells via the biotrophic interfacial complex and use a common unconventional secretory pathway. Taken together, our study reveals major changes in gene expression associated with blast disease and identifies a diverse repertoire of effectors critical for successful infection.

Genetics and breeding for resistance against four leaf spot diseases in wheat (Triticum aestivum L.)

In wheat, major yield losses are caused by a variety of diseases including rusts, spike diseases, leaf spot and root diseases. The genetics of resistance against all these diseases have been studied in great detail and utilized for breeding resistant cultivars. The resistance against leaf spot diseases caused by each individual necrotroph/hemi-biotroph involves a complex system involving resistance (R) genes, sensitivity (S) genes, small secreted protein (SSP) genes and quantitative resistance loci (QRLs). This review deals with resistance for the following four-leaf spot diseases: (i) Septoria nodorum blotch (SNB) caused by Parastagonospora nodorum; (ii) Tan spot (TS) caused by Pyrenophora tritici-repentis; (iii) Spot blotch (SB) caused by Bipolaris sorokiniana and (iv) Septoria tritici blotch (STB) caused by Zymoseptoria tritici.

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.

Deciphering immune responses primed by a bacterial lipopeptide in wheat towards Zymoseptoria tritici

Plant immunity induction with natural biocontrol compounds is a valuable and promising ecofriendly tool that fits with sustainable agriculture and healthy food. Despite the agroeconomic significance of wheat, the mechanisms underlying its induced defense responses remain obscure. We reveal here, using combined transcriptomic, metabolomic and cytologic approach, that the lipopeptide mycosubtilin from the beneficial bacterium Bacillus subtilis, protects wheat against Zymoseptoria tritici through a dual mode of action (direct and indirect) and that the indirect one relies mainly on the priming rather than on the elicitation of plant defense-related mechanisms. Indeed, the molecule primes the expression of 80 genes associated with sixteen functional groups during the early stages of infection, as well as the accumulation of several flavonoids during the period preceding the fungal switch to the necrotrophic phase. Moreover, genes involved in abscisic acid (ABA) biosynthesis and ABA-associated signaling pathways are regulated, suggesting a role of this phytohormone in the indirect activity of mycosubtilin. The priming-based bioactivity of mycosubtilin against a biotic stress could result from an interaction of the molecule with leaf cell plasma membranes that may mimic an abiotic stress stimulus in wheat leaves. This study provides new insights into induced immunity in wheat and opens new perspectives for the use of mycosubtilin as a biocontrol compound against Z. tritici.

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.

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.

Fungal Secondary Metabolites and Small RNAs Enhance Pathogenicity during Plant-Fungal Pathogen Interactions

Fungal plant pathogens use proteinaceous effectors as well as newly identified secondary metabolites (SMs) and small non-coding RNA (sRNA) effectors to manipulate the host plant's defense system via diverse plant cell compartments, distinct organelles, and many host genes. However, most molecular studies of plant-fungal interactions have focused on secreted effector proteins without exploring the possibly equivalent functions performed by fungal (SMs) and sRNAs, which are collectively known as "non-proteinaceous effectors". Fungal SMs have been shown to be generated throughout the plant colonization process, particularly in the early biotrophic stages of infection. The fungal repertoire of non-proteinaceous effectors has been broadened by the discovery of fungal sRNAs that specifically target plant genes involved in resistance and defense responses. Many RNAs, particularly sRNAs involved in gene silencing, have been shown to transmit bidirectionally between fungal pathogens and their hosts. However, there are no clear functional approaches to study the role of these SM and sRNA effectors. Undoubtedly, fungal SM and sRNA effectors are now a treasured land to seek. Therefore, understanding the role of fungal SM and sRNA effectors may provide insights into the infection process and identification of the interacting host genes that are targeted by these effectors. This review discusses the role of fungal SMs and sRNAs during plant-fungal interactions. It will also focus on the translocation of sRNA effectors across kingdoms, the application of cross-kingdom RNA interference in managing plant diseases and the tools that can be used to predict and study these non-proteinaceous effectors.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

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.

Stb6 mediates stomatal immunity, photosynthetic functionality, and the antioxidant system during the Zymoseptoria tritici-wheat interaction

This study offers new perspectives on the biochemical and physiological changes that occur in wheat following a gene-for-gene interaction with the fungal pathogen Zymoseptoria tritici. The Z. tritici isolate IPO323, carries AvrStb6, while ΔAvrStb6#33, lacks AvrStb6. The wheat cultivar (cv.) Shafir, bears the corresponding resistance gene Stb6. Inoculation of cv. Shafir with these isolates results in two contrasted phenotypes, offering a unique opportunity to study the immune response caused by the recognition of AvrStb6 by Stb6. We employed a variety of methodologies to dissect the physiological and biochemical events altered in cv. Shafir, as a result of the AvrStb6-Stb6 interaction. Comparative analysis of stomatal conductance demonstrated that AvrStb6-Stb6 mediates transient stomatal closures to restrict the penetration of Zymoseptoria tritici. Tracking photosynthetic functionality through chlorophyll fluorescence imaging analysis demonstrated that AvrStb6-Stb6 retains the functionality of photosynthesis apparatus by promoting Non-Photochemical Quenching (NPQ). Furthermore, the PlantCV image analysis tool was used to compare the H₂O₂ accumulation and incidence of cell death (2, 4, 8, 12, 16, and 21 dpi), over Z. tritici infection. Finally, our research shows that the AvrStb6-Stb6 interaction coordinates the expression and activity of antioxidant enzymes, both enzymatic and non-enzymatic, to counteract oxidative stress. In conclusion, the Stb6-AvrStb6 interaction in the Z. tritici-wheat pathosystem triggers transient stomatal closure and maintains photosynthesis while regulating oxidative stress.

Induction of distinct plant cell death programs by secreted proteins from the wheat pathogen Zymoseptoria tritici

Cell death processes in eukaryotes shape normal development and responses to the environment. For plant-microbe interactions, initiation of host cell death plays an important role in determining disease outcomes. Cell death pathways are frequently initiated following detection of pathogen-derived molecules which can lead to resistance or susceptibility to disease depending on pathogen lifestyle. We previously identified several small secreted proteins (SSPs) from the wheat-infecting fungus Zymoseptoria tritici that induce rapid cell death in Nicotiana benthamiana following Agrobacterium-mediated delivery and expression (agroinfiltration). Here we investigated whether the execution of host cells was mechanistically similar in response to different Z. tritici SSPs. Using RNA sequencing, we found that transient expression of four Z. tritici SSPs led to massive transcriptional reprogramming within 48 h of agroinfiltration. We observed that distinct host gene expression profiles were induced dependent on whether cell death occurs in a cell surface immune receptor-dependent or -independent manner. These gene expression profiles involved differential transcriptional networks mediated by WRKY, NAC and MYB transcription factors. In addition, differential expression of genes belonging to different classes of receptor-like proteins and receptor-like kinases was observed. These data suggest that different Z. tritici SSPs trigger differential transcriptional reprogramming in plant cells.

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.

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.

Scaling-up to understand tree-pathogen interactions: A steep, tough climb or a walk in the park?

Plants have proficient tools that allow them to survive interactions with pathogens. Upon attack, they respond with specific countermeasures, which are controlled by the immune system. However, defences can fail and this failure exposes plants to fast-spreading devastation. Trees face similar challenges to other plants and their immune system allows them to mount defences against pathogens. However, their slow growth, longevity, woodiness, and size can make trees a challenging system to study. Here, we review scientific successes in plant systems, highlight the key challenges and describe the enormous opportunities for pathology research in trees. We discuss the benefits that scaling-up our understanding on tree-pathogen interactions can provide in the fight against plant pathogenic threats.

The Zymoseptoria tritici white collar-1 gene, ZtWco-1, is required for development and virulence on wheat

The fungus Zymoseptoria tritici causes Septoria Tritici Blotch (STB), which is one of the most devastating diseases of wheat in Europe. There are currently no fully durable methods of control against Z. tritici, so novel strategies are urgently required. One of the ways in which fungi are able to respond to their surrounding environment is through the use of photoreceptor proteins which detect light signals. Although previous evidence suggests that Z. tritici can detect light, no photoreceptor genes have been characterised in this pathogen. This study characterises ZtWco-1, a predicted photoreceptor gene in Z. tritici. The ZtWco-1 gene is a putative homolog to the blue light photoreceptor from Neurospora crassa, wc-1. Z. tritici mutants with deletions in ZtWco-1 have defects in hyphal branching, melanisation and virulence on wheat. In addition, we identify the putative circadian clock gene ZtFrq in Z. tritici. This study provides evidence for the genetic regulation of light detection in Z. tritici and it open avenues for future research into whether this pathogen has a circadian clock.

Bioinspired Rhamnolipid Protects Wheat Against Zymoseptoria tritici Through Mainly Direct Antifungal Activity and Without Major Impact on Leaf Physiology

Rhamnolipids (RLs), glycolipids biosynthesized by the Pseudomonas and Burkholderia genera, are known to display various activities against a wide range of pathogens. Most previous studies on RLs focused on their direct antimicrobial activity, while only a few reports described the mechanisms by which RLs induce resistance against phytopathogens and the related fitness cost on plant physiology. Here, we combined transcriptomic and metabolomic approaches to unravel the mechanisms underlying RL-induced resistance in wheat against the hemibiotrophic fungus Zymoseptoria tritici, a major pathogen of this crop. Investigations were carried out by treating wheat plants with a bioinspired synthetic mono-RL with a 12-carbon fatty acid tail, dodecanoyl α/β-L-rhamnopyranoside (Rh-Est-C12), under both infectious and non-infectious conditions to examine its potential wheat defense-eliciting and priming bioactivities. Whereas, Rh-Est-C12 conferred to wheat a significant protection against Z. tritici (41% disease severity reduction), only a slight effect of this RL on wheat leaf gene expression and metabolite accumulation was observed. A subset of 24 differentially expressed genes (DEGs) and 11 differentially accumulated metabolites (DAMs) was scored in elicitation modalities 2, 5, and 15 days post-treatment (dpt), and 25 DEGs and 17 DAMs were recorded in priming modalities 5 and 15 dpt. Most changes were down-regulations, and only a few DEGs and DAMs associated with resistance to pathogens were identified. Nevertheless, a transient early regulation in gene expression was highlighted at 2 dpt (e.g., genes involved in signaling, transcription, translation, cell-wall structure, and function), suggesting a perception of the RL by the plant upon treatment. Further in vitro and in planta bioassays showed that Rh-Est-C12 displays a significant direct antimicrobial activity toward Z. tritici. Taken together, our results suggest that Rh-Est-C12 confers protection to wheat against Z. tritici through direct antifungal activity and, to a lesser extent, by induction of plant defenses without causing major alterations in plant metabolism. This study provides new insights into the modes of action of RLs on the wheat-Z. tritici pathosystem and highlights the potential interest in Rh-Est-C12, a low-fitness cost molecule, to control this pathogen.

Regulatory Modules of Metabolites and Protein Phosphorylation in Arabidopsis Genotypes With Altered Sucrose Allocation

Multi-omics data sets are increasingly being used for the interpretation of cellular processes in response to environmental cues. Especially, the posttranslational modification of proteins by phosphorylation is an important regulatory process affecting protein activity and/or localization, which, in turn, can have effects on metabolic processes and metabolite levels. Despite this importance, relationships between protein phosphorylation status and metabolite abundance remain largely underexplored. Here, we used a phosphoproteomics-metabolomics data set collected at the end of day and night in shoots and roots of Arabidopsis to propose regulatory relationships between protein phosphorylation and accumulation or allocation of metabolites. For this purpose, we introduced a novel, robust co-expression measure suited to the structure of our data sets, and we used this measure to construct metabolite-phosphopeptide networks. These networks were compared between wild type and plants with perturbations in key processes of sugar metabolism, namely, sucrose export (sweet11/12 mutant) and starch synthesis (pgm mutant). The phosphopeptide-metabolite network turned out to be highly sensitive to perturbations in sugar metabolism. Specifically, KING1, the regulatory subunit of SnRK1, was identified as a primary candidate connecting protein phosphorylation status with metabolism. We additionally identified strong changes in the fatty acid network of the sweet11/12 mutant, potentially resulting from a combination of fatty acid signaling and metabolic overflow reactions in response to high internal sucrose concentrations. Our results further suggest novel protein-metabolite relationships as candidates for future targeted research.

A Laminarin-Based Formulation Protects Wheat Against Zymoseptoria tritici via Direct Antifungal Activity and Elicitation of Host Defense-Related Genes

The present study aimed to evaluate the potential of the laminarin-based formulation Vacciplant to protect and induce resistance in wheat against Zymoseptoria tritici, a major pathogen on this crop. Under greenhouse conditions, a single foliar spraying of the product 2 days before inoculation with Z. tritici reduced disease severity and pycnidium density by 42 and 45%, respectively. Vacciplant exhibited a direct antifungal activity on Z. tritici conidial germination both in vitro and in planta. Moreover, it reduced in planta substomatal colonization as well as pycnidium formation on treated leaves. Molecular investigations revealed that Vacciplant elicits but did not prime the expression of several wheat genes related to defense pathways, including phenylpropanoids (phenylalanine ammonia-lyase and chalcone synthase), octadecanoids (lipoxygenase and allene oxide synthase), and pathogenesis-related proteins (β-1,3-endoglucanase and chitinase). By contrast, it did not modulate the expression of oxalate oxidase gene involved in the reactive oxygen species metabolism. Ultrahigh-performance liquid chromatography-mass spectrometry analysis indicated limited changes in leaf metabolome after product application in both noninoculated and inoculated conditions, suggesting a low metabolic cost associated with induction of plant resistance. This study provides evidence that the laminarin-based formulation confers protection to wheat against Z. tritici through direct antifungal activity and elicitation of plant defense-associated genes.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds—namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

RNA-Seq and Gene Regulatory Network Analyses Uncover Candidate Genes in the Early Defense to Two Hemibiotrophic Colletorichum spp. in Strawberry

Two hemibiotrophic pathogens, Colletotrichum acutatum (Ca) and C. gloeosporioides (Cg), cause anthracnose fruit rot and anthracnose crown rot in strawberry (Fragaria × ananassa Duchesne), respectively. Both Ca and Cg can initially infect through a brief biotrophic phase, which is associated with the production of intracellular primary hyphae that can infect host cells without causing cell death and establishing hemibiotrophic infection (HBI) or quiescent (latent infections) in leaf tissues. The Ca and Cg HBI in nurseries and subsequent distribution of asymptomatic infected transplants to fruit production fields is the major source of anthracnose epidemics in North Carolina. In the absence of complete resistance, strawberry varieties with good fruit quality showing rate-reducing resistance have frequently been used as a source of resistance to Ca and Cg. However, the molecular mechanisms underlying the rate-reducing resistance or susceptibility to Ca and Cg are still unknown. We performed comparative transcriptome analyses to examine how rate-reducing resistant genotype NCS 10-147 and susceptible genotype ‘Chandler’ respond to Ca and Cg and identify molecular events between 0 and 48 h after the pathogen-inoculated and mock-inoculated leaf tissues. Although plant response to both Ca and Cg at the same timepoint was not similar, more genes in the resistant interaction were upregulated at 24 hpi with Ca compared with those at 48 hpi. In contrast, a few genes were upregulated in the resistant interaction at 48 hpi with Cg. Resistance response to both Ca and Cg was associated with upregulation of MLP-like protein 44, LRR receptor-like serine/threonine-protein kinase, and auxin signaling pathway, whereas susceptibility was linked to modulation of the phenylpropanoid pathway. Gene regulatory network inference analysis revealed candidate transcription factors (TFs) such as GATA5 and MYB-10, and their downstream targets were upregulated in resistant interactions. Our results provide valuable insights into transcriptional changes during resistant and susceptible interactions, which can further facilitate assessing candidate genes necessary for resistance to two hemibiotrophic Colletotrichum spp. in strawberry.

Unravelling consensus genomic regions conferring leaf rust resistance in wheat via meta-QTL analysis

Leaf rust, caused by the fungus Puccinia triticina Erikss (Pt), is a destructive disease affecting wheat (Triticum aestivum L.) and a threat to food security. Developing resistant cultivars represents a useful method of disease control, and thus, understanding the genetic basis for leaf rust resistance is required. To this end, a comprehensive bibliographic search for leaf rust resistance quantitative trait loci (QTL) was performed, and 393 QTL were collected from 50 QTL mapping studies. Afterward, a consensus map with a total length of 4,567 cM consisting of different types of markers (simple sequence repeat [SSR], diversity arrays technology [DArT], chip-based single-nucleotide polymorphism [SNP] markers, and SNP markers from genotyping-by-sequencing) was used for QTL projection, and meta-QTL (MQTL) analysis was performed on 320 QTL. A total of 75 MQTL were discovered and refined to 15 high-confidence MQTL (hcmQTL). The candidate genes discovered within the hcmQTL interval were then checked for differential expression using data from three transcriptome studies, resulting in 92 differentially expressed genes (DEGs). The expression of these genes in various leaf tissues during wheat development was explored. This study provides insight into leaf rust resistance in wheat and thereby provides an avenue for developing resistant cultivars by incorporating the most important hcmQTL.

Chickpea Roots Undergoing Colonisation by Phytophthora medicaginis Exhibit Opposing Jasmonic Acid and Salicylic Acid Accumulation and Signalling Profiles to Leaf Hemibiotrophic Models

Hemibiotrophic pathogens cause significant losses within agriculture, threatening the sustainability of food systems globally. These microbes colonise plant tissues in three phases: a biotrophic phase followed by a biotrophic-to-necrotrophic switch phase and ending with necrotrophy. Each of these phases is characterized by both common and discrete host transcriptional responses. Plant hormones play an important role in these phases, with foliar models showing that salicylic acid accumulates during the biotrophic phase and jasmonic acid/ethylene responses occur during the necrotrophic phase. The appropriateness of this model to plant roots has been challenged in recent years. The need to understand root responses to hemibiotrophic pathogens of agronomic importance necessitates further research. In this study, using the root hemibiotroph Phytophthora medicaginis, we define the duration of each phase of pathogenesis in Cicer arietinum (chickpea) roots. Using transcriptional profiling, we demonstrate that susceptible chickpea roots display some similarities in response to disease progression as previously documented in leaf plant–pathogen hemibiotrophic interactions. However, our transcriptomic results also show that chickpea roots do not conform to the phytohormone responses typically found in leaf colonisation by hemibiotrophs. We found that quantified levels of salicylic acid concentrations in root tissues decreased significantly during biotrophy while jasmonic acid concentrations were significantly induced. This study demonstrated that a wider spectrum of plant species should be investigated in the future to understand the physiological changes in plants during colonisation by soil-borne hemibiotrophic pathogens before we can better manage these economically important microbes.

The Parastagonospora nodorum necrotrophic effector SnTox5 targets the wheat gene Snn5 and facilitates entry into the leaf mesophyll

Parastagonospora nodorum is an economically important necrotrophic fungal pathogen of wheat. Parastagonospora nodorum secretes necrotrophic effectors that target wheat susceptibility genes to induce programmed cell death (PCD). In this study, we cloned and functionally validated SnTox5 and characterized its role in pathogenesis. We used whole genome sequencing, genome-wide association study (GWAS) mapping, CRISPR-Cas9-based gene disruption, gain-of-function transformation, quantitative trait locus (QTL) analysis, haplotype and isoform analysis, protein modeling, quantitative PCR, and laser confocal microscopy to validate SnTox5 and functionally characterize SnTox5. SnTox5 is a mature 16.26 kDa protein with high structural similarity to SnTox3. Wild-type and mutant P. nodorum strains and wheat genotypes of SnTox5 and Snn5, respectively, were used to show that SnTox5 not only targets Snn5 to induce PCD but also facilitates the colonization of the mesophyll layer even in the absence of Snn5. Here we show that SnTox5 facilitates the efficient colonization of the mesophyll tissue and elicits PCD specific to host lines carrying Snn5. The homology to SnTox3 and the ability of SnTox5 to facilitate the colonizing of the mesophyll also suggest a role in the suppression of host defense before PCD induction.

Resistance of the Wheat Cultivar ‘Renan’ to Septoria Leaf Blotch Explained by a Combination of Strain Specific and Strain Non-Specific QTL Mapped on an Ultra-Dense Genetic Map

Quantitative resistance is considered more durable than qualitative resistance as it does not involve major resistance genes that can be easily overcome by pathogen populations, but rather a combination of genes with a lower individual effect. This durability means that quantitative resistance could be an interesting tool for breeding crops that would not systematically require phytosanitary products. Quantitative resistance has yet to reveal all of its intricacies. Here, we delve into the case of the wheat/Septoria tritici blotch (STB) pathosystem. Using a population resulting from a cross between French cultivar Renan, generally resistant to STB, and Chinese Spring, a cultivar susceptible to the disease, we built an ultra-dense genetic map that carries 148,820 single nucleotide polymorphism (SNP) markers. Phenotyping the interaction was done with two different Zymoseptoria tritici strains with contrasted pathogenicities on Renan. A linkage analysis led to the detection of three quantitative trait loci (QTL) related to resistance in Renan. These QTL, on chromosomes 7B, 1D, and 5D, present with an interesting diversity as that on 7B was detected with both fungal strains, while those on 1D and 5D were strain-specific. The resistance on 7B was located in the region of Stb8 and the resistance on 1D colocalized with Stb19. However, the resistance on 5D was new, so further designated Stb20q. Several wall-associated kinases (WAK), nucleotide-binding and leucine-rich repeats (NB-LRR) type, and kinase domain carrying genes were present in the QTL regions, and some of them were expressed during the infection. These results advocate for a role of Stb genes in quantitative resistance and for resistance in the wheat/STB pathosystem being as a whole quantitative and polygenic.

Two Novel Dimorphism-Related Virulence Factors of Zymoseptoria tritici Identified Using Agrobacterium-Mediated Insertional Mutagenesis

Diseases caused by dimorphic phytopathogenic and systemic dimorphic fungi have markedly increased in prevalence in the last decades, and understanding the morphogenic transition to the virulent state might yield novel means of controlling dimorphic fungi. The dimorphic fungus Z. tritici causes significant economic impact on wheat production, and yet the regulation of the dimorphic switch, a key first step in successful plant colonization, is still largely unexplored in this fungus. The fungus is amenable to suppression by fungicides at this switch point, and the identification of the factors controlling the dimorphic switch provides a potential source of novel targets to control Septoria tritici blotch (STB). Inhibition of the dimorphic switch can potentially prevent penetration and avoid any damage to the host plant. The aim of the current work was to unveil genetic determinants of the dimorphic transition in Z. tritici by using a forward genetics strategy. Using this approach, we unveiled two novel factors involved in the switch to the pathogenic state and used reverse genetics and complementation to confirm the role of the novel virulence factors and further gained insight into the role of these genes, using transcriptome analysis via RNA-Seq. The transcriptomes generated potentially contain key determinants of the dimorphic transition.

Plant PhysioSpace: a robust tool to compare stress response across plant species

Generalization of transcriptomics results can be achieved by comparison across experiments. This generalization is based on integration of interrelated transcriptomics studies into a compendium. Such a focus on the bigger picture enables both characterizations of the fate of an organism and distinction between generic and specific responses. Numerous methods for analyzing transcriptomics datasets exist. Yet, most of these methods focus on gene-wise dimension reduction to obtain marker genes and gene sets for, for example, pathway analysis. Relying only on isolated biological modules might result in missing important confounders and relevant contexts. We developed a method called Plant PhysioSpace, which enables researchers to compute experimental conditions across species and platforms without a priori reducing the reference information to specific gene sets. Plant PhysioSpace extracts physiologically relevant signatures from a reference dataset (i.e. a collection of public datasets) by integrating and transforming heterogeneous reference gene expression data into a set of physiology-specific patterns. New experimental data can be mapped to these patterns, resulting in similarity scores between the acquired data and the extracted compendium. Because of its robustness against platform bias and noise, Plant PhysioSpace can function as an inter-species or cross-platform similarity measure. We have demonstrated its success in translating stress responses between different species and platforms, including single-cell technologies. We have also implemented two R packages, one software and one data package, and a Shiny web application to facilitate access to our method and precomputed models.

Remarkable recent changes in the genetic diversity of the avirulence gene AvrStb6 in global populations of the wheat pathogen Zymoseptoria tritici

Septoria tritici blotch (STB), caused by the fungus Zymoseptoria tritici, is one of the most economically important diseases of wheat. Recently, both factors of a gene-for-gene interaction between Z. tritici and wheat, the wheat receptor-like kinase Stb6 and the Z. tritici secreted effector protein AvrStb6, have been identified. Previous analyses revealed a high diversity of AvrStb6 haplotypes present in earlier Z. tritici isolate collections, with up to c.18% of analysed isolates possessing the avirulence isoform of AvrStb6 identical to that originally identified in the reference isolate IPO323. With Stb6 present in many commercial wheat cultivars globally, we aimed to assess potential changes in AvrStb6 genetic diversity and the incidence of haplotypes allowing evasion of Stb6-mediated resistance in more recent Z. tritici populations. Here we show, using targeted resequencing of AvrStb6, that this gene is universally present in field isolates sampled from major wheat-growing regions of the world in 2013-2017. However, in contrast to the data from previous AvrStb6 population studies, we report a complete absence of the originally described avirulence isoform of AvrStb6 amongst modern Z. tritici isolates. Moreover, a remarkably small number of haplotypes, each encoding AvrStb6 protein isoforms conditioning virulence on Stb6-containing wheat, were found to predominate among modern Z. tritici isolates. A single virulence isoform of AvrStb6 was found to be particularly abundant throughout the global population. These findings indicate that, despite the ability of Z. tritici to sexually reproduce on resistant hosts, AvrStb6 avirulence haplotypes tend to be eliminated in subsequent populations.

In Planta Transcriptome and Proteome Profiles of Spongospora subterranea in Resistant and Susceptible Host Environments Illuminates Regulatory Principles Underlying Host-Pathogen Interaction

Simple Summary

Infections of potato tubers and roots by Spongospora subterranea result in powdery scab and root diseases. Losses due to infections with S. subterranea are substantial in most potato-growing regions of the world with no fully effective treatments available. Understanding the gene regulation of pathogens in their host is dependent on multidimensional datasets. In this study, we profiled the transcriptome and proteome of S. subterranea within the susceptible and resistant host. Enzyme activity and nucleic acid metabolism appear to be important to the virulence of S. subterranea. Our results provide a good resource for future functional studies of powdery scab and might be useful in S. subterranea inoculum management.

Abstract

Spongospora subterranea is an obligate biotrophic pathogen, causing substantial economic loss to potato industries globally. Currently, there are no fully effective management strategies for the control of potato diseases caused by S. subterranea. To further our understanding of S. subterranea biology during infection, we characterized the transcriptome and proteome of the pathogen during the invasion of roots of a susceptible and a resistant potato cultivar. A total of 7650 transcripts from S. subterranea were identified in the transcriptome analysis in which 1377 transcripts were differentially expressed between two cultivars. In proteome analysis, we identified 117 proteins with 42 proteins significantly changed in comparisons between resistant and susceptible cultivars. The functional annotation of transcriptome data indicated that the gene ontology terms related to the transportation and actin processes were induced in the resistant cultivar. The downregulation of enzyme activity and nucleic acid metabolism in the resistant cultivar suggests a probable influence of these processes in the virulence of S. subterranea. The protein analysis results indicated that the majority of differentially expressed proteins were related to the metabolic processes and transporter activity. The present study provides a comprehensive molecular insight into the multiple layers of gene regulation that contribute to S. subterranea infection and development in planta and illuminates the role of host immunity in affecting pathogen responses.

Characterization of Effector-Target Interactions in Necrotrophic Pathosystems Reveals Trends and Variation in Host Manipulation

Great strides have been made in defining the details of the plant defense response involving biotrophic fungal and bacterial pathogens. The groundwork for the current model was laid by H.H. Flor and others who defined the gene-for-gene hypothesis, which is now known to involve effector-triggered immunity (ETI). PAMP-triggered immunity (PTI) is also a highly effective response to most pathogens because of the recognition of common pathogen molecules by pattern recognition receptors. In this article, we consider the three pathogens that make up the foliar disease complex of wheat, Zymoseptoria tritici, Pyrenophora tritici-repentis, and Parastagonospora nodorum, to review the means by which necrotrophic pathogens circumvent, or outright hijack, the ETI and PTI pathways to cause disease.

Reprogramming of the wheat transcriptome in response to infection with Claviceps purpurea, the causal agent of ergot

BACKGROUND

Ergot, caused by the fungal pathogen Claviceps purpurea, infects the female flowers of a range of cereal crops, including wheat. To understand the interaction between C. purpurea and hexaploid wheat we undertook an extensive examination of the reprogramming of the wheat transcriptome in response to C. purpurea infection through floral tissues (i.e. the stigma, transmitting and base ovule tissues of the ovary) and over time.

RESULTS

C. purpurea hyphae were observed to have grown into and down the stigma at 24 h (H) after inoculation. By 48H hyphae had grown through the transmitting tissue into the base, while by 72H hyphae had surrounded the ovule. By 5 days (D) the ovule had been replaced by fungal tissue. Differential gene expression was first observed at 1H in the stigma tissue. Many of the wheat genes differentially transcribed in response to C. purpurea infection were associated with plant hormones and included the ethylene (ET), auxin, cytokinin, gibberellic acid (GA), salicylic acid and jasmonic acid (JA) biosynthetic and signaling pathways. Hormone-associated genes were first detected in the stigma and base tissues at 24H, but not in the transmitting tissue. Genes associated with GA and JA pathways were seen in the stigma at 24H, while JA and ET-associated genes were identified in the base at 24H. In addition, several defence-related genes were differential expressed in response to C. purpurea infection, including antifungal proteins, endocytosis/exocytosis-related proteins, NBS-LRR class proteins, genes involved in programmed cell death, receptor protein kinases and transcription factors. Of particular interest was the identification of differential expression of wheat genes in the base tissue well before the appearance of fungal hyphae, suggesting that a mobile signal, either pathogen or plant-derived, is delivered to the base prior to colonisation.

CONCLUSIONS

Multiple host hormone biosynthesis and signalling pathways were significantly perturbed from an early stage in the wheat - C. purpurea interaction. Differential gene expression at the base of the ovary, ahead of arrival of the pathogen, indicated the potential presence of a long-distance signal modifying host gene expression.

The identification of a transposon affecting the asexual reproduction of the wheat pathogen Zymoseptoria tritici

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is a fungal wheat pathogen that causes significant global yield losses. Within Z. tritici populations, quantitative differences in virulence among different isolates are commonly observed; however, the genetic components that underpin these differences remain elusive. In this study, intraspecific comparative transcriptomic analysis was used to identify candidate genes that contribute to differences in virulence on the wheat cultivar WW2449. This led to the identification of a multicopy gene that was not expressed in the high-virulence isolate when compared to the medium- and low-virulence isolates. Further investigation suggested this gene resides in a 7.9-kb transposon. Subsequent long-read sequencing of the isolates used in the transcriptomic analysis confirmed that this gene did reside in an active Class II transposon, which is composed of four genes named REP9-1 to -4. Silencing and overexpression of REP9-1 in two distinct genetic backgrounds demonstrated that its expression alone reduces the number of pycnidia produced by Z. tritici during infection. The REP9-1 gene identified within a Class II transposon is the first discovery of a gene in a transposable element that influences the virulence of Z. tritici. This discovery adds further complexity to genetic loci that contribute to quantitative virulence in this important pathogen.

Differential Physiological Prerequisites and Gene Expression Profiles of Conidial Anastomosis Tube and Germ Tube Formation in Colletotrichum gloeosporioides

The conidia of a hemibiotrophic fungus, Colletotrichum gloeosporioides, can conventionally form a germ tube (GT) and develop into a fungal colony. Under certain conditions, they tend to get connected through a conidial anastomosis tube (CAT) to share the nutrients. CAT fusion is believed to be responsible for the generation of genetic variations in few asexual fungi, which appears problematic for effective fungal disease management. The physiological and molecular requirements underlying the GT formation versus CAT fusion remained underexplored. In the present study, we have deciphered the physiological prerequisites for GT formation versus CAT fusion in C. gloeosporioides. GT formation occurred at a high frequency in the presence of nutrients, while CAT fusion was found to be higher in the absence of nutrients. Younger conidia were found to form GT efficiently, while older conidia preferentially formed CAT. Whole transcriptome analysis of GT and CAT revealed highly differential gene expression profiles, wherein 11,050 and 9786 genes were differentially expressed during GT formation and CAT fusion, respectively. A total of 1567 effector candidates were identified; out of them, 102 and 100 were uniquely expressed during GT formation and CAT fusion, respectively. Genes coding for cell wall degrading enzymes, germination, hyphal growth, host-fungus interaction, and virulence were highly upregulated during GT formation. Meanwhile, genes involved in stress response, cell wall remodeling, membrane transport, cytoskeleton, cell cycle, and cell rescue were highly upregulated during CAT fusion. To conclude, the GT formation and CAT fusion were found to be mutually exclusive processes, requiring differential physiological conditions and sets of DEGs in C. gloeosporioides. This study will help in understanding the basic CAT biology in emerging fungal model species of the genus Colletotrichum.

Three LysM effectors of Zymoseptoria tritici collectively disarm chitin-triggered plant immunity

Chitin is a major structural component of fungal cell walls and acts as a microbe-associated molecular pattern (MAMP) that, on recognition by a plant host, triggers the activation of immune responses. To avoid the activation of these responses, the Septoria tritici blotch (STB) pathogen of wheat, Zymoseptoria tritici, secretes LysM effector proteins. Previously, the LysM effectors Mg1LysM and Mg3LysM were shown to protect fungal hyphae against host chitinases. Furthermore, Mg3LysM, but not Mg1LysM, was shown to suppress chitin-induced reactive oxygen species (ROS) production. Whereas initially a third LysM effector gene was disregarded as a presumed pseudogene, we now provide functional data to show that this gene also encodes a LysM effector, named Mgx1LysM, that is functional during wheat colonization. While Mg3LysM confers a major contribution to Z. tritici virulence, Mgx1LysM and Mg1LysM contribute to Z. tritici virulence with smaller effects. All three LysM effectors display partial functional redundancy. We furthermore demonstrate that Mgx1LysM binds chitin, suppresses the chitin-induced ROS burst, and is able to protect fungal hyphae against chitinase hydrolysis. Finally, we demonstrate that Mgx1LysM is able to undergo chitin-induced polymerization. Collectively, our data show that Z. tritici utilizes three LysM effectors to disarm chitin-triggered wheat immunity.

Saccharin Provides Protection and Activates Defense Mechanisms in Wheat Against the Hemibiotrophic Pathogen Zymoseptoria tritici

Plant resistance inducers are among the most promising alternatives to develop sustainable crop protection. Here, we examined the ability of saccharin, a metabolite derived from probenazole, to protect wheat against Zymoseptoria tritici, the most frequently occurring and damaging foliar pathogen on this crop. The experiments were performed in the greenhouse by treating seedlings of the wheat cultivar 'Alixan' with 15 mM of saccharin 2 days before challenge inoculation with the Z. tritici pathogenic strain T02596. Foliar application of saccharin resulted in 77% lower disease severity than in nontreated control plants. In vitro and in planta assays showed that saccharin did not exhibit any direct antifungal effect on spore germination or hyphal growth. Molecular investigations from 2 to 7 days posttreatment (dpt) revealed that saccharin treatment upregulates the expression of genes encoding for lipoxygenase (LOX) at all sampled time points and pathogenesis-related protein 1 (PR1) at 7 dpt, in both noninfectious and infectious contexts, as well as peroxidase (POX2) in noninfectious conditions. However, saccharin did not induce significant change in the expression of PAL gene encoding for phenylalanine ammonia-lyase. Our findings report for the first time the potential of saccharin to confer protection in wheat against Z. tritici through an elicitation and priming of LOX and PR gene-related defense pathways. Additional investigations would provide a better deciphering of defense mechanisms activated by this molecule in wheat against Z. tritici.