The wheat-Septoria conflict: a new front opening up?

Septoria tritici blotch resistance gene Stb15 encodes a lectin receptor-like kinase

Septoria tritici blotch (STB), caused by the Dothideomycete fungus Zymoseptoria tritici, is one of the most damaging diseases of bread wheat (Triticum aestivum)1 and the target of costly fungicide applications2. In line with the fungus's apoplastic lifestyle, STB resistance genes isolated to date encode receptor-like kinases (RLKs) including a wall-associated kinase (Stb6) and a cysteine-rich kinase (Stb16q)3,4. Here we used genome-wide association studies on a diverse panel of 300 whole-genome shotgun-sequenced wheat landraces (WatSeq consortium5) to identify a 99-kb region containing six candidates for the Stb15 resistance gene. Mutagenesis and transgenesis confirmed a gene encoding an intronless G-type lectin RLK as Stb15. The characterization of Stb15 exemplifies the unexpected diversity of RLKs conferring Z. tritici resistance in wheat.

Exploitation of the genetic potential of Thinopyrum and Agropyron genera to protect wheat from diseases and environmental stresses

Common wheat is one of the most important food crops in the world. Grain harvests can be increased by reducing losses from diseases and environmental stresses. The tertiary gene pool, including Thinopyrum spp., is a valuable resource for increasing genetic diversity and wheat resistance to fungal diseases and abiotic stresses. Distant hybridization between wheat and Thinopyrum spp. began in the 1920s in Russia, and later continued in different countries. The main results were obtained using the species Th. ponticum and Th. intermedium. Additionally, introgression material was created based on Th. elongatum, Th. bessarabicum, Th. junceiforme, Agropyron cristatum. The results of introgression for resistance to diseases (leaf, stem, and stripe rusts; powdery mildew; Fusarium head blight; and Septoria blotch) and abiotic stresses (drought, extreme temperatures, and salinity) to wheat was reviewed. Approaches to improving the agronomic properties of introgression breeding material (the use of irradiation, ph-mutants and compensating Robertsonian translocations) were described. The experience of long-term use in the world of a number of genes from the tertiary gene pool in protecting wheat from leaf and stem rust was observed. Th. ponticum is a nonhost for Puccinia triticina (Ptr) and P. graminis f. sp. tritici (Pgt) and suppresses the development of rust fungi on the plant surface. Wheat samples with the tall wheatgrass genes Lr19, Lr38, Sr24, Sr25 and Sr26 showed defence mechanisms similar to nonhosts resistance. Their influence led to disruption of the development of surface infection structures and fungal death when trying to penetrate the stomata (prehaustorial resistance or stomatal immunity). Obviously, a change in the chemical properties of fungal surface structures of races virulent to Lr19, Lr24, Sr24, Sr25, and Sr26 leads to a decrease in their adaptability to the environment. This possibly determined the durable resistance of cultivars to leaf and stem rusts in different regions. Alien genes with a similar effect are of interest for breeding cultivars with durable resistance to rust diseases and engineering crops with the help of molecular technologies.

The Zymoseptoria tritici Avirulence Factor AvrStb6 Accumulates in Hyphae Close to Stomata and Triggers a Wheat Defense Response Hindering Fungal Penetration

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is one of Europe's most damaging wheat pathogens, causing significant economic losses. Genetic resistance is a common strategy to control the disease, Stb6 being a resistance gene used for more than 100 years in Europe. This study investigates the molecular mechanisms underlying Stb6-mediated resistance. Utilizing confocal microscopy imaging, we determined that Z. tritici epiphytic hyphae mainly accumulate the corresponding avirulence factor AvrStb6 in close proximity to stomata. Consequently, the progression of AvrStb6-expressing avirulent strains is hampered during penetration. The fungal growth inhibition co-occurs with a transcriptional reprogramming in wheat characterized by an induction of immune responses, genes involved in stomatal regulation, and cell wall-related genes. Overall, we shed light on the gene-for-gene resistance mechanisms in the wheat-Z. tritici pathosystem at the cytological and transcriptomic level, and our results highlight that stomatal penetration is a critical process for pathogenicity and resistance. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

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.

A potential alternative to fungicides using actives-free (meth)acrylate polymers for protection of wheat crops from fungal attachment and infection

Fungicidal compounds are actives widely used for crop protection from fungal infection, but they can also kill beneficial organisms, enter the food chain and promote resistant pathogen strains from overuse. Here we report the first field crop trial of homopolymer materials that prevent fungal attachment, showing successful crop protection via an actives-free approach. In the trial, formulations containing two candidate polymers were applied to young wheat plants that were subject to natural infection with the wheat pathogen Zymoseptoria tritici. A formulation containing one of the candidate polymers, poly(di(ethylene glycol) ethyl ether acrylate) (abbreviated DEGEEA), produced a significant reduction (26%) in infection of the crop by Z. tritici, delivering protection against fungal infection that compared favourably with three different commercially established fungicide programmes tested in parallel. Furthermore, the sprayed polymers did not negatively affect wheat growth. The two lead polymer candidates were initially identified by bio-performance testing using in vitro microplate- and leaf-based assays and were taken forward successfully into a programme to optimize and scale-up their synthesis and compound them into a spray formulation. Therefore, the positive field trial outcome has also established the validity of the smaller-scale, laboratory-based bioassay data and scale-up methodologies used. Because fungal attachment to plant surfaces is a first step in many crop infections, this non-eluting polymer: (i) now offers significant potential to deliver protection against fungal attack, while (ii) addressing the fourth and aligning with the eleventh principles of green chemistry by using chemical products designed to preserve efficacy of function while reducing toxicity. A future focus should be to develop the material properties for this and other applications including other fungal pathogens.

Unlocking the story of resistance to Zymoseptoria tritici in Tunisian old durum wheat germplasm based on population structure analysis

BACKGROUND

Septoria tritici blotch (STB) remains a significant obstacle to durum wheat cultivation on a global scale. This disease remains a challenge for farmers, researchers, and breeders, who are collectively dedicated to reduce its damage and improve wheat resistance. Tunisian durum wheat landraces have been recognized as valuable genetic ressources that exhibit resistance to biotic and abiotic stresses and therefore play a crucial role in breeding program aimed at creating new wheat varieties resistant to fungal diseases as STB, as well as adapted to climate change constraints.

RESULTS

A total of 366 local durum wheat accessions were assessed for resistance to two virulent Tunisian isolates of Zymoseptoria tritici Tun06 and TM220 under field conditions. Population structure analysis of the durum wheat accessions, performed with 286 polymorphic SNPs (PIC > 0.3) covering the entire genome, identified three genetic subpopulations (GS1, GS2 and GS3) with 22% of admixed genotypes. Interestingly, all of the resistant genotypes were among GS2 or admixed with GS2.

CONCLUSIONS

This study revealed the population structure and the genetic distribution of the resistance to Z. tritici in the Tunisian durum wheat landraces. Accessions grouping pattern reflected the geographical origins of the landraces. We suggested that GS2 accessions were mostly derived from eastern Mediterranean populations, unlike GS1 and GS3 that originated from the west. Resistant GS2 accessions belonged to landraces Taganrog, Sbei glabre, Richi, Mekki, Badri, Jneh Khotifa and Azizi. Furthermore, we suggested that admixture contributed to transmit STB resistance from GS2 resistant landraces to initially susceptible landraces such as Mahmoudi (GS1), but also resulted in the loss of resistance in the case of GS2 susceptible Azizi and Jneh Khotifa accessions.

Baseline sensitivity of European Zymoseptoria tritici populations to the complex III respiration inhibitor fenpicoxamid

BACKGROUND

Fenpicoxamid is a recently developed fungicide belonging to the quinone inside inhibitor (QiI) group. This is the first fungicide within this group to be active against the Zymoseptoria tritici, which causes Septoria tritici blotch on wheat. The occurrence of pre-existing resistance mechanisms was monitored, using sensitivity assays and Illumina sequencing, in Z. tritici populations sampled in multiple European countries before the introduction of fenpicoxamid.

RESULTS

Although differences in sensitivity to all three fungicides tested (fenpicoxamid, fentin chloride and pyraclostrobin) existed between the isolate collections, no alterations associated with QiI resistance were detected. Among the isolates, a range in sensitivity to fenpicoxamid was observed (ratio between most sensitive/least sensitive = 53.1), with differences between the most extreme isolates when tested in planta following limited fenpicoxamid treatment. Sensitivity assays using fentin chloride suggest some of the observed differences in fenpicoxamid sensitivity are associated with multi-drug resistance. Detailed monitoring of the wider European population using Illumina-based partial sequencing of the Z. tritici also only detected the presence of G143A, with differences in frequencies of this alteration observed across the region.

CONCLUSIONS

This study provides a baseline sensitivity for European Z. tritici populations to fenpicoxamid. Target-site resistance appears to be limited or non-existing in European Z. tritici populations prior to the introduction of fenpicoxamid. Non-target site resistance mechanisms exist, but their impact in the field is predicted to be limited. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

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.

Spatio-temporal distribution of DMI and SDHI fungicide resistance of Zymoseptoria tritici throughout Europe based on frequencies of key target-site alterations

BACKGROUND

Over the past decade, demethylation inhibitor (DMI) and succinate dehydrogenase inhibitor (SDHI) fungicides have been extensively used to control to septoria tritici blotch, caused by Zymoseptoria tritici on wheat. This has led to the development and selection of alterations in the target-site enzymes (CYP51 and SDH, respectively).

RESULTS

Taking advantage of newly and previously developed qPCR assays, the frequency of key alterations associated with DMI (CYP51-S524T) and SDHI (SDHC-T79N/I, C-N86S and C-H152R) resistance was assessed in Z. tritici-infected wheat leaf samples collected from commercial crops (n = 140) across 14 European countries prior to fungicide application in the spring of 2019. This revealed the presence of a West to East gradient in the frequencies of the most common key alterations conferring azole (S524T) and SDHI resistance (T79N and N86S), with the highest frequencies measured in Ireland and Great Britain. These observations were corroborated by sequencing (CYP51 and SDH subunits) and sensitivity phenotyping (prothioconazole-desthio and fluxapyroxad) of Z. tritici isolates collected from a selection of field samples. Additional sampling made at the end of the 2019 season confirmed the continued increase in frequency of the targeted alterations. Investigations on historical leaf DNA samples originating from different European countries revealed that the frequency of all key alterations (except C-T79I) has been gradually increasing over the past decade.

CONCLUSION

Whilst these alterations are quickly becoming dominant in Ireland and Great Britain, scope still exists to delay their selection throughout the wider European population, emphasizing the need for the implementation of fungicide antiresistance measures. © 2021 Society of Chemical Industry.

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.

Development of a PacBio Long-Read Sequencing Assay for High Throughput Detection of Fungicide Resistance in Zymoseptoria tritici

Fungicide resistance has become a challenging problem in management of Septoria tritici blotch (STB), caused by Zymoseptoria tritici, the most destructive disease of winter wheat throughout western and northern Europe. To ensure the continued effectiveness of those fungicides currently used, it is essential to monitor the development and spread of such resistance in field populations of the pathogen. Since resistance to the key families of fungicides used for STB control (demethyalation inhibitors or azoles, succinate dehydrogenase inhibitors or SDHIs and Quinone outside Inhibitors or QoIs) is conferred through target-site mutations, the potential exists to monitor resistance through the molecular detection of alterations in the target site genes. As more efficient fungicides were developed and applied, the pathogen has continuously adapted through accumulating multiple target-site alterations. In order to accurately monitor these changes in field populations, it is therefore becoming increasingly important to completely sequence the targeted genes. Here we report the development of a PacBio assay that facilitates the multiplex amplification and long-read sequencing of the target gene(s) for the azole (CYP51), SDHI (Sdh B, C, and D), and QoI (cytochrome b) fungicides. The assay was developed and optimised using three Irish Z. tritici collections established in spring 2017, which capture the range of fungicide resistance present in modern European populations of Z. tritici. The sequences obtained through the PacBio assay were validated using traditional Sanger sequencing and in vitro sensitivity screenings. To further exploit the long-read and high throughput potential of PacBio sequencing, an additional nine housekeeping genes (act, BTUB, cal, cyp, EF1, GAPDH, hsp80-1, PKC, TFC1) were sequenced and used to provide comprehensive Z. tritici strain genotyping.

Mixed infections alter transmission potential in a fungal plant pathogen

Infections by more than one strain of a pathogen predominate under natural conditions. Mixed infections can have significant, though often unpredictable, consequences for overall virulence, pathogen transmission and evolution. However, effects of mixed infection on disease development in plants often remain unclear and the critical factors that determine the outcome of mixed infections remain unknown. The fungus Zymoseptoria tritici forms genetically diverse infections in wheat fields. Here, for a range of pathogen traits, we experimentally decompose the infection process to determine how the outcomes and consequences of mixed infections are mechanistically realized. Different strains of Z. tritici grow in close proximity and compete in the wheat apoplast, resulting in reductions in growth of individual strains and in pathogen reproduction. We observed different outcomes of competition at different stages of the infection. Overall, more virulent strains had higher competitive ability during host colonization, and less virulent strains had higher transmission potential. We showed that within‐host competition can have a major effect on infection dynamics and pathogen population structure in a pathogen and host genotype‐specific manner. Consequently, mixed infections likely have a major effect on the development of septoria tritici blotch epidemics and the evolution of virulence in Z. tritici.

Biochemical characterisation of class III biotin protein ligases from Botrytis cinerea and Zymoseptoria tritici

Biotin protein ligase (BPL) is an essential enzyme in all kingdoms of life, making it a potential target for novel anti-infective agents. Whilst bacteria and archaea have simple BPL structures (class I and II), the homologues from certain eukaryotes such as mammals, insects and yeast (class III) have evolved a more complex structure with a large extension on the N-terminus of the protein in addition to the conserved catalytic domain. The absence of atomic resolution structures of any class III BPL hinders structural and functional analysis of these enzymes. Here, two new class III BPLs from agriculturally important moulds Botrytis cinerea and Zymoseptoria tritici were characterised alongside the homologue from the prototypical yeast Saccharomyces cerevisiae. Circular dichroism and ion mobility-mass spectrometry analysis revealed conservation of the overall tertiary and secondary structures of all three BPLs, corresponding with the high sequence similarity. Subtle structural differences were implied by the different thermal stabilities of the enzymes and their varied Michaelis constants for their interactions with ligands biotin, MgATP, and biotin-accepting substrates from different species. The three BPLs displayed different preferences for fungal versus bacterial protein substrates, providing further evidence that class III BPLs have a 'substrate validation' activity for selecting only appropriate proteins for biotinylation. Selective, potent inhibition of these three BPLs was demonstrated despite sequence and structural homology. This highlights the potential for targeting BPL for novel, selective antifungal therapies against B. cinerea, Z. tritici and other fungal species.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

Insights into the resistance of a synthetically-derived wheat to Septoria tritici blotch disease: less is more

BACKGROUND

Little is known about the initial, symptomless (latent) phase of the devastating wheat disease Septoria tritici blotch. However, speculations as to its impact on fungal success and disease severity in the field have suggested that a long latent phase is beneficial to the host and can reduce inoculum build up in the field over a growing season. The winter wheat cultivar Stigg is derived from a synthetic hexaploid wheat and contains introgressions from wild tetraploid wheat Triticum turgidum subsp. dicoccoides, which contribute to cv. Stigg's exceptional STB resistance, hallmarked by a long latent phase. We compared the early transcriptomic response to Zymoseptoria tritici of cv. Stigg to a susceptible wheat cultivar, to elucidate the mechanisms of and differences in pathogen recognition and disease response in these two hosts.

RESULTS

The STB-susceptible cultivar Longbow responds to Z. tritici infection with a stress response, including activation of hormone-responsive transcription factors, post translational modifications, and response to oxidative stress. The activation of key genes associated with these pathways in cv. Longbow was independently observed in a second susceptible wheat cultivar based on an independent gene expression study. By comparison, cv. Stigg is apathetic in response to STB, and appears to fail to activate a range of defence pathways that cv. Longbow employs. Stigg also displays some evidence of sub-genome bias in its response to Z. tritici infection, whereas the susceptible cv. Longbow shows even distribution of Z. tritici responsive genes across the three wheat sub-genomes.

CONCLUSIONS

We identify a suite of disease response genes that are involved in early pathogen response in susceptible wheat cultivars that may ultimately lead to susceptibility. In comparison, we hypothesise that rather than an active defence response to stave off disease progression, cv. Stigg's defence strategy is molecular lethargy, or a lower-amplitude of pathogen recognition that may stem from cv. Stigg's wild wheat-derived ancestry. Overall, we present insights into cv. Stigg's exceptional resistance to STB, and present key biological processes for further characterisation in this pathosystem.

Quantification of In Planta Zymoseptoria tritici Progression Through Different Infection Phases and Related Association with Components of Aggressiveness

In planta growth of Zymoseptoria tritici, causal agent of Septoria tritici blotch of wheat, during the infection process has remained an understudied topic due to the long symptomless latent period before the emergence of fruiting bodies. In this study, we attempted to understand the relationship between in planta growth of Z. tritici relative to the primary components of aggressiveness, i.e., latent period and pycnidia coverage in regard to contrasting host resistance. We tested isolates collected from Ireland against the susceptible cultivar Gallant and cultivar Stigg, which has strong partial resistance. A clear isolate-host interaction effect (F = 3.018; P = 0.005, and F = 6.008; P < 0.001) for latent period and pycnidia coverage, respectively, was identified. Furthermore, during the early infection phase of latency from 5 to 11 days postinoculation (dpi), in planta growth rate of fungal biomass was significantly (F = 30.06; P < 0.001) more affected by host resistance than isolate specificity (F = 1.27; P = 0.27), indicating the importance of host resistance in the early infection phase. In planta Z. tritici growth rates in cultivar Gallant spiked between 11 and 16 dpi followed by a continuous fall onward, whereas in cultivar Stigg it was slowly progressive in nature. From correlation and regression analysis, we found that the in planta growth rate preceding the average latent period of cultivar Gallant has more influence on latency duration and pycnidia production. Likewise, correlation between component of aggressiveness and in planta growth rate of pathogen supports our understanding of aggressiveness to be driven by the pathogen's multiplication capacity within host tissue.

Wheat pathogen Zymoseptoria tritici N-myristoyltransferase inhibitors: on-target antifungal activity and an unusual metabolic defense mechanism

Zymoseptoria tritici is the causative agent of Septoria tritici blotch (STB), which costs billions of dollars annually to major wheat-producing countries in terms of both fungicide use and crop loss. Agricultural pathogenic fungi have acquired resistance to most commercially available fungicide classes, and the rate of discovery and development of new fungicides has stalled, demanding new approaches and insights. Here we investigate a potential mechanism of targeting an important wheat pathogen Z. tritici via inhibition of N-myristoyltransferase (NMT). We characterize Z. tritici NMT biochemically for the first time, profile the in vivo Z. tritici myristoylated proteome and identify and validate the first Z. tritici NMT inhibitors. Proteomic investigation of the downstream effects of NMT inhibition identified an unusual and novel mechanism of defense against chemical toxicity in Z. tritici through the application of comparative bioinformatics to deconvolute function from the previously largely unannotated Z. tritici proteome. Research into novel fungicidal modes-of-action is essential to satisfy an urgent unmet need for novel fungicide targets, and we anticipate that this study will serve as a useful proteomics and bioinformatics resource for researchers studying Z. tritici.

Wheat Encodes Small, Secreted Proteins That Contribute to Resistance to Septoria Tritici Blotch

During plant-pathogen interactions, pathogens secrete many rapidly evolving, small secreted proteins (SSPs) that can modify plant defense and permit pathogens to colonize plant tissue. The fungal pathogen Zymoseptoria tritici is the causal agent of Septoria tritici blotch (STB), one of the most important foliar diseases of wheat, globally. Z. tritici is a strictly apoplastic pathogen that can secrete numerous proteins into the apoplast of wheat leaves to promote infection. We sought to determine if, during STB infection, wheat also secretes small proteins into the apoplast to mediate the recognition of pathogen proteins and/or induce defense responses. To explore this, we developed an SSP-discovery pipeline to identify small, secreted proteins from wheat genomic data. Using this pipeline, we identified 6,998 SSPs, representing 2.3% of all proteins encoded by the wheat genome. We then mined a microarray dataset, detailing a resistant and susceptible host response to STB, and identified 141 Z. tritici- responsive SSPs, representing 4.7% of all proteins encoded by Z. tritici - responsive genes. We demonstrate that a subset of these SSPs have a functional signal peptide and can interact with Z. tritici SSPs. Transiently silencing two of these wheat SSPs using virus-induced gene silencing (VIGS) shows an increase in susceptibility to STB, confirming their role in defense against Z. tritici.

Fungicide Sensitivity Shifting of Zymoseptoria tritici in the Finnish-Baltic Region and a Novel Insertion in the MFS1 Promoter

Septoria tritici blotch (STB) is caused by the ascomycete Zymoseptoria tritici and one of the predominating diseases in wheat (Triticum aestivum) in Europe. The control of STB is highly reliant on frequent fungicide applications. The primary objective of this study was to assess sensitivity levels of Z. tritici to different fungicide groups. The fungicides included in this study were epoxiconazole, prothioconazole-desthio, tebuconazole, and fluxapyroxad. A panel of 63 isolates from Estonia, Latvia, and Lithuania, and 10 isolates from Finland were tested. Fungicide sensitivity testing was carried out as a bioassay analyzing single pycnidium isolates on different fungicide concentrations. The average EC₅₀ value in Baltic countries and Finland to epoxiconazole was high ranging from 1.04 to 2.19 ppm. For prothioconazole-desthio and tebuconazole, EC₅₀ varied from 0.01 to 0.24 ppm, and 1.25 to 18.23 ppm, respectively. The average EC₅₀ value for fluxapyroxad varied from 0.07 to 0.33 ppm. To explain the range of sensitivity, the samples were analyzed for CYP51 and Sdh mutations, as well as cytb G143A, CYP51 overexpression, and multidrug resistance (MDR). Frequencies of ZtCYP51 mutations D134G, V136A/C, A379G, I381V, and S524T in the Finnish-Baltic region were lower than in other European countries, but have increased compared to previous years. The frequency of cytb G143A conferring strobilurin resistance also augmented to 50-70% in the Z. tritici populations from Estonia, Finland, Latvia, and Lithuania. No Sdh mutations were found in this study, and neither strains of MDR phenotypes. However, we found a strain harboring a previously unknown transposon insertion in the promoter of the MFS1 gene, involved in drug efflux and multi-drug resistance. This new insert, however, does not confer an MDR phenotype to the strain.

Tramesan Elicits Durum Wheat Defense against the Septoria Disease Complex

The Septoria Leaf Blotch Complex (SLBC), caused by the two ascomycetes Zymoseptoria tritici and Parastagonospora nodorum, can reduce wheat global yearly yield by up to 50%. In the last decade, SLBC incidence has increased in Italy; notably, durum wheat has proven to be more susceptible than common wheat. Field fungicide treatment can efficiently control these pathogens, but it leads to the emergence of resistant strains and adversely affects human and animal health and the environment. Our previous studies indicated that active compounds produced by Trametes versicolor can restrict the growth of mycotoxigenic fungi and the biosynthesis of their secondary metabolites (e.g., mycotoxins). Specifically, we identified Tramesan: a 23 kDa α-heteropolysaccharide secreted by T. versicolor that acts as a pro-antioxidant molecule in animal cells, fungi, and plants. Foliar-spray of Tramesan (3.3 μM) on SLBC-susceptible durum wheat cultivars, before inoculation of causal agents of Stagonospora Nodorum Blotch (SNB) and Septoria Tritici Blotch (STB), significantly decreased disease incidence both in controlled conditions (SNB: -99%, STB: -75%) and field assays (SNB: -25%, STB: -30%). We conducted these tests were conducted under controlled conditions as well as in field. We showed that Tramesan increased the levels of jasmonic acid (JA), a plant defense-related hormone. Tramesan also increased the early expression (24 hours after inoculation - hai) of plant defense genes such as PR4 for SNB infected plants, and RBOH, PR1, and PR9 for STB infected plants. These results suggest that Tramesan protects wheat by eliciting plant defenses, since it has no direct fungicidal activity. In field experiments, the yield of durum wheat plants treated with Tramesan was similar to that of healthy untreated plots. These results encourage the use of Tramesan to protect durum wheat against SLBC.

Small RNA Bidirectional Crosstalk During the Interaction Between Wheat and Zymoseptoria tritici

Cross-kingdom RNA interference (RNAi) has been shown to play important roles during plant-pathogen interactions, and both plants and pathogens can use small RNAs (sRNAs) to silence genes in each other. This bidirectional cross-kingdom RNAi was still unexplored in the wheat-Zymoseptoria tritici pathosystem. Here, we performed a detailed analysis of the sRNA bidirectional crosstalk between wheat and Z. tritici. Using a combination of small RNA sequencing (sRNA-seq) and microRNA sequencing (mRNA-seq), we were able to identify known and novel sRNAs and study their expression and their action on putative targets in both wheat and Z. tritici. We predicted the target genes of all the sRNAs in either wheat or Z. tritici transcriptome and used degradome analysis to validate the cleavage of these gene transcripts. We could not find any clear evidence of a cross-kingdom RNAi acting by mRNA cleavage in this pathosystem. We also found that the fungal sRNA enrichment was lower in planta than during in vitro growth, probably due to the lower expression of the only Dicer gene of the fungus during plant infection. Our results support the recent finding that Z. tritici sRNAs cannot play important roles during wheat infection. However, we also found that the fungal infection induced wheat sRNAs regulating the expression of specific wheat genes, including auxin-related genes, as an immune response. These results indicate a role of sRNAs in the regulation of wheat defenses during Z. tritici infection. Our findings contribute to improve our understanding of the interactions between wheat and Z. tritici.

In-Field Detection and Quantification of Septoria Tritici Blotch in Diverse Wheat Germplasm Using Spectral-Temporal Features

Hyperspectral remote sensing holds the potential to detect and quantify crop diseases in a rapid and non-invasive manner. Such tools could greatly benefit resistance breeding, but their adoption is hampered by i) a lack of specificity to disease-related effects and ii) insufficient robustness to variation in reflectance caused by genotypic diversity and varying environmental conditions, which are fundamental elements of resistance breeding. We hypothesized that relying exclusively on temporal changes in canopy reflectance during pathogenesis may allow to specifically detect and quantify crop diseases while minimizing the confounding effects of genotype and environment. To test this hypothesis, we collected time-resolved canopy hyperspectral reflectance data for 18 diverse genotypes on infected and disease-free plots and engineered spectral-temporal features representing this hypothesis. Our results confirm the lack of specificity and robustness of disease assessments based on reflectance spectra at individual time points. We show that changes in spectral reflectance over time are indicative of the presence and severity of Septoria tritici blotch (STB) infections. Furthermore, the proposed time-integrated approach facilitated the delineation of disease from physiological senescence, which is pivotal for efficient selection of STB-resistant material under field conditions. A validation of models based on spectral-temporal features on a diverse panel of 330 wheat genotypes offered evidence for the robustness of the proposed method. This study demonstrates the potential of time-resolved canopy reflectance measurements for robust assessments of foliar diseases in the context of resistance breeding.

A novel starch-binding laccase from the wheat pathogen Zymoseptoria tritici highlights the functional diversity of ascomycete laccases

BACKGROUND

Laccases are multicopper oxidases, which are assigned into auxiliary activity family 1 (AA1) in the CAZy database. These enzymes, catalyzing the oxidation of phenolic and nonphenolic substrates coupled to reduction of O2 to H2O, are increasingly attractive as eco-friendly oxidation biocatalysts. Basidiomycota laccases are well characterized due to their potential in de-lignification of lignocellulose. By contrast, insight into the biochemical diversity of Ascomycota counterparts from saprophytes and plant pathogens is scarce.

RESULTS

Here, we report the properties of the laccase from the major wheat pathogen Zymoseptoria tritici (ZtrLac1A), distinguished from common plant fungal pathogens by an apoplastic infection strategy. We demonstrate that ZtrLac1A is appended to a functional starch-binding module and displays an activity signature disfavoring relatively apolar phenolic redox mediators as compared to the related biochemically characterized laccases. By contrast, the redox potential of ZtrLac1A (370 mV vs. SHE) is similar to ascomycetes counterparts. The atypical specificity is consistent with distinctive sequence substitutions and insertions in loops flanking the T1 site and the enzyme C-terminus compared to characterized laccases.

CONCLUSIONS

ZtrLac1A is the first reported modular laccase appended to a functional starch-specific carbohydrate binding module of family 20 (CBM20). The distinct specificity profile of ZtrLac1A correlates to structural differences in the active site region compared to previously described ascomycetes homologues. These differences are also highlighted by the clustering of the sequence of ZtrLac1A in a distinct clade populated predominantly by plant pathogens in the phylogenetic tree of AA1 laccases. The possible role of these laccases in vivo merits further investigations. These findings expand our toolbox of laccases for green oxidation and highlight the binding functionality of CBM-appended laccases as versatile immobilization tags.

Sexual reproduction of Zymoseptoria tritici on durum wheat in Tunisia revealed by presence of airborne inoculum, fruiting bodies and high levels of genetic diversity

Septoria tritici blotch (STB) caused by the heterothallic ascomycete Zymoseptoria tritici is currently one of the most devastating diseases of wheat worldwide. The extent of sexual reproduction of this pathogen is well documented on bread wheat, but not on durum wheat. The objective of the present study was to quantify the occurrence of Z. tritici sexual reproduction on durum wheat in the Tunisian environment. The assessment was undertaken using a triple approach combining fruiting body assessment, ascospore trapping and population genetic analyses. The results highlighted the formation of pseudothecia on leaves and stubble from the autumn until the end of the growing season. Likewise, qPCR monitoring highlighted a constant release of Z. tritici airborne inoculum during the wheat-growing season, with a peak of production at the end of the season. Genetic investigations using microsatellites revealed high levels of gene and genotypic diversities, an equal distribution of mating types, and a lack of genetic clustering within and between growing seasons. Taken together, these findings indicate that Z. tritici undergoes sexual reproduction on durum wheat in Tunisia at least to the same extent than on bread wheat in Western Europe, and that the dry and warm climate does not affect the mating process of the fungus. Frequent occurrence of sexual reproduction is a valuable knowledge to take into account in STB control strategies on durum wheat.

Genetic Dissection of Resistance to the Three Fungal Plant Pathogens Blumeria graminis, Zymoseptoria tritici, and Pyrenophora tritici-repentis Using a Multiparental Winter Wheat Population

Bread wheat (Triticum aestivum L.) is one of the world’s most important crop species. The development of new varieties resistant to multiple pathogens is an ongoing task in wheat breeding, especially in times of increasing demand for sustainable agricultural practices. Despite this, little is known about the relations between various fungal disease resistances at the genetic level, and the possible consequences for wheat breeding strategies. As a first step to fill this gap, we analyzed the genetic relations of resistance to the three fungal diseases – powdery mildew (PM), septoria tritici blotch (STB), and tan spot (TS) – using a winter wheat multiparent advanced generation intercross population. Six, seven, and nine QTL for resistance to PM, STB, and TS, respectively, were genetically mapped. Additionally, 15 QTL were identified for the three agro-morphological traits plant height, ear emergence time, and leaf angle distribution. Our results suggest that resistance to STB and TS on chromosome 2B is conferred by the same genetic region. Furthermore, we identified two genetic regions on chromosome 1AS and 7AL, which are associated with all three diseases, but not always in a synchronal manner. Based on our results, we conclude that parallel marker-assisted breeding for resistance to the fungal diseases PM, STB, and TS appears feasible. Knowledge of the genetic co-localization of alleles with contrasting effects for different diseases, such as on chromosome 7AL, allows the trade-offs of selection of these regions to be better understood, and ultimately determined at the genic level.

Genome-wide Association Mapping and Prediction of Adult Stage Septoria tritici Blotch Infection in European Winter Wheat via High-Density Marker Arrays

blotch (STB) caused by the fungus is a devastating foliar disease of wheat ( L.) that can lead to substantial yield losses. Quantitative genetic resistance has been proposed as a durable strategy for STB control. In this study, we dissected the genetic basis of STB infection in 371 European wheat varieties based on 35k and 90k single nucleotide polymorphism marker arrays. The phenotypic data analyses suggested that large genetic variance exists for STB infection with a broad-sense heritability of 0.78. Genome-wide association studies (GWAS) propose the highly quantitative nature of STB infection with potential associations on chromosomes 1A, 1B, 2D, 4A, 5A, 6A, 6D, 7A, and 7B. Increased marker density in GWAS by combining markers from both arrays helped to detect additional markers explaining increased genotypic variance. Linkage disequilibrium analyses revealed genes with a possible role in disease resistance. The potential of genomic prediction (GP) assessed via two models accounting for additive effects and additive plus epistatic interactions among the loci suggested the possibility of genomic selection for improved STB resistance. Genomic prediction results also indicated that the higher-order epistatic interactions are not abundant and that both marker platforms are equally suitable for GP of STB infection. Our results provide further understanding of the quantitative genetic nature of STB infection, serve as a resource for marker-assisted breeding, and highlight the potential of genomic selection for improved STB resistance.

Genetic Structure of Zymoseptoria tritici in Northern France at Region, Field, Plant, and Leaf Layer Scales

Population genetic structure of the worldwide-distributed wheat pathogen Zymoseptoria tritici has been extensively studied at large geographical scales, but to a much less extent at small or local spatial scales. A total of 627 single-conidial fungal isolates were sampled from several locations in northern France (Hauts-de-France Region) to assess fungal genetic structure at region, field, plant, and leaf layer scales, using highly polymorphic microsatellite markers and mating type idiomorphs. Important and overall similar levels of both gene and genotype diversities (gene diversity values of ≥0.44 and haplotype frequencies of ≥94%) were found at all the examined scales. Such rates of diversity are likely due to an active sexual recombination in the investigated areas, as revealed by equal proportions of the two mating types scored in all sampled populations. Interestingly, a rare occurrence of clones among lesions from a same leaf, as well as among leaves from different plant leaf layers (e.g., upper versus lower leaves), was highlighted, indicating that ascospores contribute much more than expected to Z. tritici epidemics, compared with pycnidiospores. Population structure and analyses of molecular variance revealed significant genetic differentiation at the regional scale (G_ST = 0.23) and, as expected, not at the other more local scales (G_ST ≤ 0.01). Further analyses using Bayesian and unweighted neighbor-joining statistical methods detected six genetic clusters within the regional population, overall distributed according to the locations from which the isolates were sampled. Neither clear directional relative migration linked to the geographical distribution of the locations, nor isolation by distance, were observed. Separate evolutionary trajectories caused by selection and adaptations to habitat heterogeneity could be the main forces shaping such structuration. This study provides new insights into the epidemiology and the genetic structure of Z. tritici at small local and, for the first time, at single plant and leaf layer scales. Such findings would be helpful in implementing effective control strategies.

Hyperspectral Canopy Sensing of Wheat Septoria Tritici Blotch Disease

Producing quantitative and reliable measures of crop disease is essential for resistance breeding, but is challenging and time consuming using traditional phenotyping methods. Hyperspectral remote sensing has shown potential for the detection of plant diseases, but its utility for phenotyping large and diverse populations of plants under field conditions requires further evaluation. In this study, we collected canopy hyperspectral data from 335 wheat varieties using a spectroradiometer, and we investigated the use of canopy reflectance for detecting the Septoria tritici blotch (STB) disease and for quantifying the severity of infection. Canopy- and leaf-level infection metrics of STB based on traditional visual assessments and automated analyses of leaf images were used as ground truth data. Results showed (i) that canopy reflectance and the selected spectral indices show promise for quantifying STB infections, and (ii) that the normalized difference water index (NDWI) showed the best performance in detecting STB compared to other spectral indices. Moreover, partial least squares (PLS) regression models allowed for an improvement in the prediction of STB metrics. The PLS discriminant analysis (PLSDA) model calibrated based on the spectral data of four reference varieties was able to discriminate between the diseased and healthy canopies among the 335 varieties with an accuracy of 93% (Kappa = 0.60). Finally, the PLSDA model predictions allowed for the identification of wheat genotypes that are potentially more susceptible to STB, which was confirmed by the STB visual assessment. This study demonstrates the great potential of using canopy hyperspectral remote sensing to improve foliar disease assessment and to facilitate plant breeding for disease resistance.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

Comparative Transcriptomics Reveals How Wheat Responds to Infection by Zymoseptoria tritici

The fungus Zymoseptoria tritici causes septoria tritici blotch (STB) on wheat, an important disease globally and the most damaging wheat disease in Europe. Despite the global significance of STB, the molecular basis of wheat defense against Z. tritici is poorly understood. Here, we use a comparative transcriptomic study to investigate how wheat responds to infection by four distinct strains of Z. tritici. We examined the response of wheat across the entire infection cycle, identifying both shared responses to the four strains and strain-specific responses. We found that the early asymptomatic phase is characterized by strong upregulation of genes encoding receptor-like kinases and pathogenesis-related proteins, indicating the onset of a defense response. We also identified genes that were differentially expressed among the four fungal strains, including genes related to defense. Genes involved in senescence were induced during both the asymptomatic phase and at late stages of infection, suggesting manipulation of senescence processes by both the plant and the pathogen. Our findings illustrate the need, when identifying important genes affecting disease resistance in plants, to include multiple pathogen strains.

Evolutionary analyses of the avirulence effector AvrStb6 in global populations of Zymoseptoria tritici identify candidate amino acids involved in recognition

We analysed the population genetic diversity of AvrStb6, the first avirulence gene cloned from the wheat pathogen Zymoseptoria tritici, using 142 Z. tritici strains sampled from four wheat fields growing on three continents. Although AvrStb6 was located in a recombination hotspot, it was found in every strain, with 71 polymorphic sites that produced 41 distinct DNA haplotypes encoding 30 AvrStb6 protein isoforms. An AvrStb6 homologue was found in the closest known relative, Z. pseudotritici, but not in three other closely related Zymoseptoria species, indicating that this gene has emerged in Zymoseptoria quite recently. Two AvrStb6 homologues with nucleotide similarities greater than 70% were identified on chromosome 10 in all Z. tritici isolates, suggesting that AvrStb6 belongs to a multigene family of candidate effectors that has expanded recently through gene duplication. The AvrStb6 sequences exhibited strong evidence for non-neutral evolution, including a large number of non-synonymous mutations, with significant positive diversifying selection operating on nine of the 82 codons. It appears that balancing selection is operating across the entire gene in natural field populations. There was also evidence for co-evolving codons within the gene that may reflect compensatory mutations associated with the evasion of recognition by Stb6. Intragenic recombination also appears to have affected the diversity of AvrStb6.

Pangenome analyses of the wheat pathogen Zymoseptoria tritici reveal the structural basis of a highly plastic eukaryotic genome

Background

Structural variation contributes substantially to polymorphism within species. Chromosomal rearrangements that impact genes can lead to functional variation among individuals and influence the expression of phenotypic traits. Genomes of fungal pathogens show substantial chromosomal polymorphism that can drive virulence evolution on host plants. Assessing the adaptive significance of structural variation is challenging, because most studies rely on inferences based on a single reference genome sequence.

Results

We constructed and analyzed the pangenome of Zymoseptoria tritici, a major pathogen of wheat that evolved host specialization by chromosomal rearrangements and gene deletions. We used single-molecule real-time sequencing and high-density genetic maps to assemble multiple genomes. We annotated the gene space based on transcriptomics data that covered the infection life cycle of each strain. Based on a total of five telomere-to-telomere genomes, we constructed a pangenome for the species and identified a core set of 9149 genes. However, an additional 6600 genes were exclusive to a subset of the isolates. The substantial accessory genome encoded on average fewer expressed genes but a larger fraction of the candidate effector genes that may interact with the host during infection. We expanded our analyses of the pangenome to a worldwide collection of 123 isolates of the same species. We confirmed that accessory genes were indeed more likely to show deletion polymorphisms and loss-of-function mutations compared to core genes.

Conclusions

The pangenome construction of a highly polymorphic eukaryotic pathogen showed that a single reference genome significantly underestimates the gene space of a species. The substantial accessory genome provides a cradle for adaptive evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0457-4) contains supplementary material, which is available to authorized users.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

Comparison of Models and Whole-Genome Profiling Approaches for Genomic-Enabled Prediction of Septoria Tritici Blotch, Stagonospora Nodorum Blotch, and Tan Spot Resistance in Wheat

The leaf spotting diseases in wheat that include Septoria tritici blotch (STB) caused by , Stagonospora nodorum blotch (SNB) caused by , and tan spot (TS) caused by pose challenges to breeding programs in selecting for resistance. A promising approach that could enable selection prior to phenotyping is genomic selection that uses genome-wide markers to estimate breeding values (BVs) for quantitative traits. To evaluate this approach for seedling and/or adult plant resistance (APR) to STB, SNB, and TS, we compared the predictive ability of least-squares (LS) approach with genomic-enabled prediction models including genomic best linear unbiased predictor (GBLUP), Bayesian ridge regression (BRR), Bayes A (BA), Bayes B (BB), Bayes Cπ (BC), Bayesian least absolute shrinkage and selection operator (BL), and reproducing kernel Hilbert spaces markers (RKHS-M), a pedigree-based model (RKHS-P) and RKHS markers and pedigree (RKHS-MP). We observed that LS gave the lowest prediction accuracies and RKHS-MP, the highest. The genomic-enabled prediction models and RKHS-P gave similar accuracies. The increase in accuracy using genomic prediction models over LS was 48%. The mean genomic prediction accuracies were 0.45 for STB (APR), 0.55 for SNB (seedling), 0.66 for TS (seedling) and 0.48 for TS (APR). We also compared markers from two whole-genome profiling approaches: genotyping by sequencing (GBS) and diversity arrays technology sequencing (DArTseq) for prediction. While, GBS markers performed slightly better than DArTseq, combining markers from the two approaches did not improve accuracies. We conclude that implementing GS in breeding for these diseases would help to achieve higher accuracies and rapid gains from selection.

A small secreted protein in Zymoseptoria tritici is responsible for avirulence on wheat cultivars carrying the Stb6 resistance gene

Zymoseptoria tritici is the causal agent of Septoria tritici blotch, a major pathogen of wheat globally and the most damaging pathogen of wheat in Europe. A gene-for-gene (GFG) interaction between Z. tritici and wheat cultivars carrying the Stb6 resistance gene has been postulated for many years, but the genes have not been identified. We identified AvrStb6 by combining quantitative trait locus mapping in a cross between two Swiss strains with a genome-wide association study using a natural population of c. 100 strains from France. We functionally validated AvrStb6 using ectopic transformations. AvrStb6 encodes a small, cysteine-rich, secreted protein that produces an avirulence phenotype on wheat cultivars carrying the Stb6 resistance gene. We found 16 nonsynonymous single nucleotide polymorphisms among the tested strains, indicating that AvrStb6 is evolving very rapidly. AvrStb6 is located in a highly polymorphic subtelomeric region and is surrounded by transposable elements, which may facilitate its rapid evolution to overcome Stb6 resistance. AvrStb6 is the first avirulence gene to be functionally validated in Z. tritici, contributing to our understanding of avirulence in apoplastic pathogens and the mechanisms underlying GFG interactions between Z. tritici and wheat.

Comparative Transcriptome Analyses in Zymoseptoria tritici Reveal Significant Differences in Gene Expression Among Strains During Plant Infection

Zymoseptoria tritici is an ascomycete fungus that causes Septoria tritici blotch, a globally distributed foliar disease on wheat. Z. tritici populations are highly polymorphic and exhibit significant quantitative variation for virulence. Despite its importance, the genes responsible for quantitative virulence in this pathogen remain largely unknown. We investigated the expression profiles of four Z. tritici strains differing in virulence in an experiment conducted under uniform environmental conditions. Transcriptomes were compared at four different infection stages to characterize the regulation of gene families thought to be involved in virulence and to identify new virulence factors. The major components of the fungal infection transcriptome showed consistent expression profiles across strains. However, strain-specific regulation was observed for many genes, including some encoding putative virulence factors. We postulate that strain-specific regulation of virulence factors can determine the outcome of Z. tritici infections. We show that differences in gene expression may be major determinants of virulence variation among Z. tritici strains, adding to the already known contributions to virulence variation based on differences in gene sequence and gene presence/absence polymorphisms.

Linear Correlation Analysis of Zymoseptoria tritici Aggressiveness with In Vitro Growth Rate

Zymoseptoria tritici is a globally distributed plant-pathogenic fungus causing Septoria tritici blotch of wheat. In this study, the in vitro growth rates and aggressiveness of 141 genetically distinct isolates sampled from four wheat fields on three continents were assessed to determine the association of these two ecological parameters. Aggressiveness was assessed on two spring wheat cultivars (‘Toronit’ and ‘Greina’) in a greenhouse using percentages of leaf area covered by lesions and pycnidia. We found a positive correlation between aggressiveness of pathogen strains on the two cultivars, consistent with a quantitative and host-nonspecific interaction in this pathosystem. We also found a positive correlation between aggressiveness and average growth rate at two temperatures, suggesting that in vitro pathogen growth rate may make a significant contribution to pathogen aggressiveness.

How Knowledge of Pathogen Population Biology Informs Management of Septoria Tritici Blotch

Zymoseptoria tritici (previously Mycosphaerella graminicola) causes Septoria tritici blotch (STB) on wheat. The population biology of Z. tritici has been exceptionally well characterized as a result of intensive studies conducted over nearly 30 years. These studies provided important insights into the biology, epidemiology and evolutionary history of Z. tritici that will prove useful for management of STB. The well-documented, rapid adaptation of Z. tritici populations to fungicide applications and deployment of wheat cultivars carrying both major gene and quantitative resistance reflects the high evolutionary potential predicted by the large effective population size, high degree of gene flow and high levels of recombination found in field populations of Z. tritici globally. QST studies that assessed the global diversity for several important quantitative traits confirmed the adaptive potential of field populations and laid the groundwork for quantitative trait loci (QTL) mapping studies. QTL mapping elucidated the genetic architecture of each trait and led to identification of candidate genes affecting fungicide resistance, thermal adaptation, virulence, and host specialization. The insights that emerged through these analyses of Z. tritici population biology can now be used to generate actionable disease management strategies aimed at sustainably reducing losses due to STB. The high evolutionary potential found in field populations of Z. tritici requires deployment of a corresponding dynamically diverse set of control measures that integrate cultural, chemical, biological and resistance breeding strategies. In this review, we describe and prioritize STB control strategies based on current knowledge of Z. tritici population biology and propose a future research agenda oriented toward long-term STB management.

Comparative transcriptomic analyses of Zymoseptoria tritici strains show complex lifestyle transitions and intraspecific variability in transcription profiles

Zymoseptoria tritici causes Septoria tritici blotch (STB) on wheat. The disease interaction is characterized by clearly defined temporal phases of infection, ultimately resulting in the death of host tissue. Zymoseptoria tritici is a highly polymorphic species with significant intraspecific variation in virulence profiles. We generated a deep transcriptomic sequencing dataset spanning the entire time course of an infection using a previously uncharacterized, highly virulent Z. tritici strain isolated from a Swiss wheat field. We found that seven clusters of gene transcription profiles explained the progression of the infection. The earliest highly up-regulated genes included chloroperoxidases, which may help the fungus cope with plant defences. The onset of necrotrophy was characterized by a concerted up-regulation of proteases, plant cell wall-degrading enzymes and lipases. Functions related to nutrition and growth characterized late necrotrophy and the transition to saprotrophic growth on dead plant tissue. We found that the peak up-regulation of genes essential for mating coincided with the necrotrophic phase. We performed an intraspecies comparative transcriptomics analysis using a comparable time course infection experiment of the genome reference isolate IPO323. Major components of the fungal infection transcriptome were conserved between the two strains. However, individual small, secreted proteins, proteases and cell wall-degrading enzymes showed strongly differentiated transcriptional profiles between isolates. Our analyses illustrate that successful STB infections involve complex transcriptomic remodelling to up-regulate distinct gene functions. Heterogeneity in transcriptomes among isolates may explain some of the considerable variation in virulence and host specialization found within the species.

An Improved Method for Measuring Quantitative Resistance to the Wheat Pathogen Zymoseptoria tritici Using High-Throughput Automated Image Analysis

Zymoseptoria tritici causes Septoria tritici blotch (STB) on wheat. An improved method of quantifying STB symptoms was developed based on automated analysis of diseased leaf images made using a flatbed scanner. Naturally infected leaves (n = 949) sampled from fungicide-treated field plots comprising 39 wheat cultivars grown in Switzerland and 9 recombinant inbred lines (RIL) grown in Oregon were included in these analyses. Measures of quantitative resistance were percent leaf area covered by lesions, pycnidia size and gray value, and pycnidia density per leaf and lesion. These measures were obtained automatically with a batch-processing macro utilizing the image-processing software ImageJ. All phenotypes in both locations showed a continuous distribution, as expected for a quantitative trait. The trait distributions at both sites were largely overlapping even though the field and host environments were quite different. Cultivars and RILs could be assigned to two or more statistically different groups for each measured phenotype. Traditional visual assessments of field resistance were highly correlated with quantitative resistance measures based on image analysis for the Oregon RILs. These results show that automated image analysis provides a promising tool for assessing quantitative resistance to Z. tritici under field conditions.

A new class of fatty acid allene oxide formed by the DOX-P450 fusion proteins of human and plant pathogenic fungi, C. immitis and Z. tritici

Linoleate dioxygenase-cytochrome P450 (DOX-CYP) fusion enzymes are common in pathogenic fungi. The DOX domains form hydroperoxy metabolites of 18:2n-6, which can be transformed by the CYP domains to 1,2- or 1,4-diols, epoxy alcohols, or to allene oxides. We have characterized two novel allene oxide synthases (AOSs), namely, recombinant 8R-DOX-AOS of Coccidioides immitis (causing valley fever) and 8S-DOX-AOS of Zymoseptoria tritici (causing septoria tritici blotch of wheat). The 8R-DOX-AOS oxidized 18:2n-6 sequentially to 8R-hydroperoxy-9Z,12Z-octadecadienoic acid (8R-HPODE) and to an allene oxide, 8R(9)-epoxy-9,12Z-octadecadienoic acid, as judged from the accumulation of the α-ketol, 8S-hydroxy-9-oxo-12Z-octadecenoic acid. The 8S-DOX-AOS of Z. tritici transformed 18:2n-6 sequentially to 8S-HPODE and to an α-ketol, 8R-hydroxy-9-oxo-12Z-octadecenoic acid, likely formed by hydrolysis of 8S(9)-epoxy-9,12Z-octadecadienoic acid. The 8S-DOX-AOS oxidized [8R-(2)H]18:2n-6 to 8S-HPODE with retention of the (2)H-label, suggesting suprafacial hydrogen abstraction and oxygenation in contrast to 8R-DOX-AOS. Both enzymes oxidized 18:1n-9 and 18:3n-3 to α-ketols, but the catalysis of the 8R- and 8S-AOS domains differed. 8R-DOX-AOS transformed 9R-HPODE to epoxy alcohols, but 8S-DOX-AOS converted 9S-HPODE to an α-ketol (9-hydroxy-10-oxo-12Z-octadecenoic acid) and epoxy alcohols in a ratio of ∼1:2. Whereas all fatty acid allene oxides described so far have a conjugated diene impinging on the epoxide, the allene oxides formed by 8-DOX-AOS are unconjugated.

QTL mapping of temperature sensitivity reveals candidate genes for thermal adaptation and growth morphology in the plant pathogenic fungus Zymoseptoria tritici

Different thermal environments impose strong, differential selection on populations, leading to local adaptation, but the genetic basis of thermal adaptation is poorly understood. We used quantitative trait locus (QTL) mapping in the fungal wheat pathogen Zymoseptoria tritici to study the genetic architecture of thermal adaptation and identify candidate genes. Four wild-type strains originating from the same thermal environment were crossed to generate two mapping populations with 263 (cross 1) and 261 (cross 2) progeny. Restriction site-associated DNA sequencing was used to genotype 9745 (cross 1) and 7333 (cross 2) single-nucleotide polymorphism markers segregating within the mapping population. Temperature sensitivity was assessed using digital image analysis of colonies growing at two different temperatures. We identified four QTLs for temperature sensitivity, with unique QTLs found in each cross. One QTL had a logarithm of odds score >11 and contained only six candidate genes, including PBS2, encoding a mitogen-activated protein kinase kinase associated with low temperature tolerance in Saccharomyces cerevisiae. This and other QTLs showed evidence for pleiotropy among growth rate, melanization and growth morphology, suggesting that many traits can be correlated with thermal adaptation in fungi. Higher temperatures were highly correlated with a shift to filamentous growth among the progeny in both crosses. We show that thermal adaptation has a complex genetic architecture, with natural populations of Z. tritici harboring significant genetic variation for this trait. We conclude that Z. tritici populations have the potential to adapt rapidly to climate change and expand into new climatic zones.

The impact of Septoria tritici Blotch disease on wheat: An EU perspective

•

Zymospetoria tritici is a threat to wheat production in the EU Z.t.’s plastic genome increases the potential severity of this threat in the future.

•

Climate change may also affect the risk from Z.t.

•

We estimate the spore numbers produced by Z.t. during each infection cycle.

•

We calculate 1) the economic value of wheat in the three main EU producers 2) the cost of and economic return for fungicide treatment of wheat vs Z.t.

Zymoseptoria tritici is the causal agent of one of the European Union’s most devastating foliar diseases of wheat: Septoria tritici Blotch (STB). It is also a notable pathogen of wheat grown in temperate climates throughout the world. In this commentary, we highlight the importance of STB on wheat in the EU. To better understand STB, it is necessary to consider the host crop, the fungal pathogen and their shared environment. Here, we consider the fungus per se and its interaction with its host and then focus on a more agricultural overview of the impact STB on wheat. We consider the climatic and weather factors which influence its spread and severity, allude to the agricultural practices which may mitigate or enhance its impact on crop yields, and evaluate the economic importance of wheat as a food and animal feed crop in the UK and EU. Finally, we estimate the cost of STB disease to EU agriculture.

Genetics of resistance to Zymoseptoria tritici and applications to wheat breeding

This paper reviews current knowledge about genes for resistance to Septoria tritici blotch (STB) of wheat, caused by Zymoseptoria tritici (formerly Mycosphaerella graminicola). These genes can be placed into two classes, although a few may have characteristics of both classes. Qualitative resistance is controlled by genes which control large fractions of genetic variation, 21 of which have been discovered and mapped so far. Most of them have been shown to be genotype-specific, being effective against the minority of Z. tritici isolates which are avirulent, and Stb6 has been shown to control a gene-for-gene relationship. Most qualitative resistances are unlikely to be durable and some formerly effective genes have been overcome by the evolution of pathogen virulence. Quantitative resistance is generally controlled by genes with small-to-moderate effects on STB. They have generally weaker specificity than qualitative genes and have provided more durable resistance. 89 genome regions carrying quantitative trait loci (QTL) or meta-QTL have been identified to date. Some QTL have been mapped at or near loci of qualitative genes, especially Stb6, which is present in several sources of resistance. Another gene of particular interest is Stb16q, which has been effective against all Z. tritici isolates tested so far. In addition to resistance, the susceptibility of wheat cultivars to STB can also be reduced by disease escape traits, some of which may be undesirable in breeding. The fundamental requirements for breeding for STB-resistance are genetic diversity for resistance in wheat germplasm and a field trial site at which STB epidemics occur regularly and effective selection can be conducted for resistance combined with other desirable traits. If these are in place, knowledge of resistance genes can be applied to improving control of STB.

Potential and limits of whole genome prediction of resistance to Fusarium head blight and Septoria tritici blotch in a vast Central European elite winter wheat population

Fusarium head blight and Septoria tritici blotch resistances are complex traits and can be improved efficiently by genomic selection modeling main and epistatic effects. Enhancing the resistance against Fusarium head blight (FHB) and Septoria tritici blotch (STB) is of central importance for a sustainable wheat production. Our study is based on a large experimental data set of 2325 inbred lines genotyped with 12,642 SNP markers and phenotyped in multi-environmental trials for FHB and STB resistance as well as for plant height. Our objectives were to (1) investigate the impact of plant height on FHB and STB severity, (2) examine the potential of marker-assisted selection, and (3) study the prediction ability of genomic selection modeling main and epistatic effects. We observed low correlations between plant height and FHB (r = -0.15; P < 0.05) as well as STB severity (r = -0.17; P < 0.05) suggesting negligible morphological resistances. Cross-validation in combination with association mapping revealed absence of large effect QTL impeding an efficient pyramiding of different resistance loci through marker-assisted selection. The prediction ability of genomic selection was high amounting to 0.6 for FHB and 0.5 for STB resistance. Therefore, genomic selection is a promising tool to improve FHB and STB resistance in wheat.

Adult-plant resistance to Septoria tritici blotch in hexaploid spring wheat

New QTL for Septoria tritici blotch detected in hexapoid spring wheat under field conditions across diverse environments. Septoria tritici blotch caused by the ascomycete fungus Zymoseptoria tritici presents a serious and consistent challenge to global wheat production. In particular the augmented use of soil management practices that leave large amounts of wheat stubble on the soil surface and global warming increases the chance of Septoria tritici blotch epidemics to emerge more frequently including in developing countries. Two recombinant inbred line populations developed from a cross between the susceptible Moroccan spring bread wheat variety 'NASMA' and the CIMMYT resistant lines, 'IAS20*5/H567.71' and 'RPB709.71/COC' were used to study the genetics and map adult-plant resistance to Septoria tritici blotch under field conditions in different environments. Resistance to Septoria tritici blotch in both populations was quantitative and overall, five across environment consistent resistance loci on chromosomes 1BS, 3AL, 5AL and 7AS were detected in the two populations. The QTL on chromosome 1BS and 7AS are likely to be allelic with the known Septoria tritici blotch genes Stb3 and Stb11. All identified QTL were additive and explained between 4 and 27 % of the phenotypic variation. Epistatic interaction was not observed. Low cost KASP assays were developed as flanking markers for all five QTL that will facilitate molecular breeding. Our study represents the first mapping effort under field conditions utilizing two spring bread wheat resistant sources evaluated over multiple environments.

The Impact of Recombination Hotspots on Genome Evolution of a Fungal Plant Pathogen

Recombination has an impact on genome evolution by maintaining chromosomal integrity, affecting the efficacy of selection, and increasing genetic variability in populations. Recombination rates are a key determinant of the coevolutionary dynamics between hosts and their pathogens. Historic recombination events created devastating new pathogens, but the impact of ongoing recombination in sexual pathogens is poorly understood. Many fungal pathogens of plants undergo regular sexual cycles, and sex is considered to be a major factor contributing to virulence. We generated a recombination map at kilobase-scale resolution for the haploid plant pathogenic fungus Zymoseptoria tritici. To account for intraspecific variation in recombination rates, we constructed genetic maps from two independent crosses. We localized a total of 10,287 crossover events in 441 progeny and found that recombination rates were highly heterogeneous within and among chromosomes. Recombination rates on large chromosomes were inversely correlated with chromosome length. Short accessory chromosomes often lacked evidence for crossovers between parental chromosomes. Recombination was concentrated in narrow hotspots that were preferentially located close to telomeres. Hotspots were only partially conserved between the two crosses, suggesting that hotspots are short-lived and may vary according to genomic background. Genes located in hotspot regions were enriched in genes encoding secreted proteins. Population resequencing showed that chromosomal regions with high recombination rates were strongly correlated with regions of low linkage disequilibrium. Hence, genes in pathogen recombination hotspots are likely to evolve faster in natural populations and may represent a greater threat to the host.