Effector discovery in the fungal wheat pathogen Zymoseptoria tritici

Fungal plant pathogens, such as Zymoseptoria tritici (formerly known as Mycosphaerella graminicola), secrete repertoires of effectors to facilitate infection or trigger host defence mechanisms. The discovery and functional characterization of effectors provides valuable knowledge that can contribute to the design of new and effective disease management strategies. Here, we combined bioinformatics approaches with expression profiling during pathogenesis to identify candidate effectors of Z. tritici. In addition, a genetic approach was conducted to map quantitative trait loci (QTLs) carrying putative effectors, enabling the validation of both complementary strategies for effector discovery. In planta expression profiling revealed that candidate effectors were up‐regulated in successive waves corresponding to consecutive stages of pathogenesis, contrary to candidates identified by QTL mapping that were, overall, expressed at low levels. Functional analyses of two top candidate effectors (SSP15 and SSP18) showed their dispensability for Z. tritici pathogenesis. These analyses reveal that generally adopted criteria, such as protein size, cysteine residues and expression during pathogenesis, may preclude an unbiased effector discovery. Indeed, genetic mapping of genomic regions involved in specificity render alternative effector candidates that do not match the aforementioned criteria, but should nevertheless be considered as promising new leads for effectors that are crucial for the Z. tritici–wheat pathosystem.     

ATP prevents Woronin bodies from sealing septal pores in unwounded cells of the fungus Zymoseptoria tritici

Septa of filamentous ascomycetes are perforated by septal pores that allow communication between individual hyphal compartments. Upon injury, septal pores are plugged rapidly by Woronin bodies (WBs), thereby preventing extensive cytoplasmic bleeding. The mechanism by which WBs translocate into the pore is not known, but it has been suggested that wound‐induced cytoplasmic bleeding “flushes” WBs into the septal opening. Alternatively, contraction of septum‐associated tethering proteins may pull WBs into the septal pore. Here, we investigate WB dynamics in the wheat pathogen Zymoseptoria tritici. Ultrastructural studies showed that 3.4 ± 0.2 WBs reside on each side of a septum and that single WBs of 128.5 ± 3.6 nm in diameter seal the septal pore (41 ± 1.5 nm). Live cell imaging of green fluorescent ZtHex1, a major protein in WBs, and the integral plasma membrane protein ZtSso1 confirms WB translocation into the septal pore. This was associated with the occasional formation of a plasma membrane “balloon,” extruding into the dead cell, suggesting that the plasma membrane rapidly seals the wounded septal pore wound. Minor amounts of fluorescent ZtHex1‐enhanced green fluorescent protein (eGFP) appeared associated with the “ballooning” plasma membrane, indicating that cytoplasmic ZtHex1‐eGFP is recruited to the extending plasma membrane. Surprisingly, in ~15% of all cases, WBs moved from the ruptured cell into the septal pore. This translocation against the cytoplasmic flow suggests that an active mechanism drives WB plugging. Indeed, treatment of unwounded and intact cells with the respiration inhibitor carbonyl cyanide m‐chlorophenyl hydrazone induced WB translocation into the pores. Moreover, carbonyl cyanide m‐chlorophenyl hydrazone treatment recruited cytoplasmic ZtHex1‐eGFP to the lateral plasma membrane of the cells. Thus, keeping the WBs out of the septal pores, in Z. tritici, is an ATP‐dependent process.     

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.     

Comprehensive proteomic analysis of Penicillium verrucosum

Mass spectrometric identification of proteins in species lacking validated sequence information is a major problem in veterinary science. In the present study, we used ochratoxin A producing Penicillium verrucosum to identify and quantitatively analyze proteins of an organism with yet no protein information available. The work presented here aimed to provide a comprehensive protein identification of P. verrucosum using shotgun proteomics. We were able to identify 3631 proteins in an “ab initio” translated database from DNA sequences of P. verrucosum. Additionally, a sequential window acquisition of all theoretical fragment‐ion spectra analysis was done to find differentially regulated proteins at two different time points of the growth curve. We compared the proteins at the beginning (day 3) and at the end of the log phase (day 12).     

Sources of resistance and susceptibility to Septoria tritici blotch of wheat

An association genetics analysis was conducted to investigate the genetics of resistance to Septoria tritici blotch, caused by the fungus Zymoseptoria tritici (alternatively Mycosphaerella graminicola), in cultivars and breeding lines of wheat (Triticum aestivum) used in the UK between 1860 and 2000. The population was tested with Diversity Array Technology (DArT) and simple‐sequence repeat (SSR or microsatellite) markers. The lines formed a single population with no evidence for subdivision, because there were several common ancestors of large parts of the pedigree. Quantitative trait loci (QTLs) controlling Septoria resistance were postulated on 11 chromosomes, but 38% of variation was not explained by the identified QTLs. Calculation of best linear unbiased predictions (BLUPs) identified lineages of spring and winter wheat carrying different alleles for resistance and susceptibility. Abundant variation in Septoria resistance may be exploited by crossing well‐adapted cultivars in different lineages to achieve transgressive segregation and thus breed for potentially durable quantitative resistance, whereas phenotypic selection for polygenic quantitative resistance should be effective in breeding cultivars with increased resistance. The most potent allele reducing susceptibility to Septoria, on chromosome arm 6AL, was associated with reduced leaf size. Genes which increase susceptibility to Septoria may have been introduced inadvertently into UK wheat breeding programmes from cultivars used to increase yield, rust resistance and eyespot resistance between the 1950s and 1980s. This indicates the need to consider trade‐offs in plant breeding when numerous traits are important and to be cautious about the use of non‐adapted germplasm.     

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.     

Variation in aggressiveness components of Zymoseptoria tritici populations in Iran

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, has emerged and evolved as a pathogen on wheat in the Fertile Crescent. Iran is located in this region, and study on the ancient pathogen populations in the country can improve our understanding on adaptive potential of aggressiveness, the role of local adaptation in shaping population structure and the involvement of selection and genetic drift in favouring aggressiveness adaptation to environments. To this aim, three aggressiveness components including days until first lesion (DUFL), days until first pycnidia (DUFP) and percentage of leaf area covered by pycnidia, in five populations of the pathogen collected from different provinces of Iran were compared in greenhouse conditions. All populations except for Golestan showed high within‐population diversity for the examined traits. No difference in aggressiveness components was found between fungal collections from Khuzestan and Fars; however, significant variation was evident among the populations originated from other provinces. Comparisons of estimated Q_ST values to F_ST indicated that genetic differentiation in pycnidial coverage has been the result of selection imposed by different variables; however, the divergence found for DUFL and DUFP has been achieved by genetic drift. The possible mechanisms involved in aggressiveness diversity of the pathogen populations and the impact of these findings on breeding programs for quantitative resistance are discussed.     

2‐DE proteomic approach to the Botrytis cinerea secretome induced with different carbon sources and plant‐based elicitors

Botrytis cinerea is a phytopathogenic fungus infecting a number of crops (tomatoes, grapes and strawberries), which has been adopted as a model system in molecular phytopathology. B. cinerea uses a wide variety of infection strategies, which are mediated by a set of genes/proteins called pathogenicity/virulence factors. Many of these factors have been described as secreted proteins, and thus the study of this sub‐proteome, the secretome, under changing circumstances can help us to understand the roles of these factors, possibly revealing new loci for the fight against the pathogen. A 2‐DE, MALDI TOF/TOF‐based approach has been developed to establish the proteins secreted to culture media supplemented with different carbon sources and plant‐based elicitors (in this study: glucose, cellulose, starch, pectin and tomato cell walls). Secreted proteins were obtained from the culture media by deoxycholate‐trichloroacetic acid/phenol extraction, and 76 spots were identified, yielding 95 positive hits that correspond to 56 unique proteins, including several known virulence factors (i.e. pectin methyl esterases, xylanases and proteases). The observed increases in secretion of proteins with established virulence‐related functions indicate that this in vitro‐induction/proteome‐mining approach is a promising strategy for discovering new pathogenicity factors and dissecting infection mechanisms in a discrete fashion.     

Proteomic analysis of the phytopathogenic fungus Botrytis cinerea during cellulose degradation

The ascomycete Botrytis cinerea is a phytopathogenic fungus infecting and causing significant yield losses in a number of crops. Moreover, in the last few years, B. cinerea has been adopted as an important model system in molecular phytopathology. In spite of these contributions, the molecular basis of the infection cycle remains unclear. Proteomic approaches have revealed significant information about the infective cycle of several pathogens, including B. cinerea. The main aim of this study is to make available a proteomic database containing a significant number of identified proteins from B. cinerea. In brief, three independent B. cinerea cultures supplemented with carboxymethylcellulose were used, and the extracted proteins were independently separated by 2‐D PAGE to obtain the proteome map from B. cinerea. Two hundred and sixty‐seven spots were selected for MALDI TOF/TOF MS analysis, resulting in 303 positive identifications, mostly representing unannotated proteins. Identified proteins were then classified into categories using the PANTHER classification system ( www.pantherdb.org ), showing the relevance of protein metabolism and modification process and oxidoreductase activity. Since cellulose is one of the major components of the plant cell wall, many of the identified proteins may have a crucial role in the pathogenicity process. In brief, this proteomic map of B. cinerea will be a useful basis for exploring the proteins involved in the infection cycle, which will in turn provide new targets for crop diagnosis and focused fungicide design.     

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.     

A proteomic approach to identifying proteins differentially expressed in conidia and mycelium of the entomopathogenic fungus Metarhizium acridum

Metarhizium spp. is an important worldwide group of entomopathogenic fungi used as an interesting alternative to chemical insecticides in programs of agricultural pest and disease vector control. Metarhizium conidia are important in fungal propagation and also are responsible for host infection. Despite their importance, several aspects of conidial biology, including their proteome, are still unknown. We have established conidial and mycelial proteome reference maps for Metarhizium acridum using two-dimensional gel electrophoresis (2-DE) and matrix-assisted laser desorption/ionization-time of flight-mass spectrometry (MALDI-TOF MS). In all, 1130±102 and 1200±97 protein spots were detected in ungerminated conidia and fast-growing mycelia, respectively. Comparison of the two protein-expression profiles reveled that only 35% of the protein spots were common to both developmental stages. Out of 94 2-DE protein spots (65 from conidia, 25 from mycelia and two common to both) analyzed using mass spectrometry, seven proteins from conidia, 15 from mycelia and one common to both stages were identified. The identified protein spots exclusive to conidia contained sequences similar to known fungal stress-protector proteins (such as heat shock proteins (HSP) and 6-phosphogluconate dehydrogenase) plus the fungal allergen Alt a 7, actin and the enzyme cobalamin-independent methionine synthase. The identified protein spots exclusive to mycelia included proteins involved in several cell housekeeping biological processes. Three proteins (HSP 90, 6-phosphogluconate dehydrogenase and allergen Alt a 7) were present in spots in conidial and mycelial gels, but they differed in their locations on the two gels.     

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.     

High-pH reversed-phase chromatography with fraction concatenation for 2D proteomic analysis

Orthogonal high-resolution separations are critical for attaining improved analytical dynamic range and protein coverage in proteomic measurements. High-pH reversed-phase liquid chromatography (RPLC), followed by fraction concatenation, affords better peptide analysis than conventional strong cation-exchange chromatography applied for 2D proteomic analysis. For example, concatenated high-pH RPLC increased identification of peptides (by 1.8-fold) and proteins (by 1.6-fold) in shotgun proteomics analyses of a digested human protein sample. Additional advantages of high-pH RPLC with fraction concatenation include improved protein sequence coverage, simplified sample processing and reduced sample losses, making this an attractive alternative to strong cation-exchange chromatography in conjunction with second-dimension low-pH RPLC for 2D proteomics analyses.     

Comprehensive proteomic analysis of human bile

Bile serves diverse functions from metabolism to transport. In addition to acids and salts, bile is composed of proteins secreted or shed by the hepatobiliary system. Although there have been previous efforts to catalog biliary proteins, an in-depth analysis of the bile proteome has not yet been reported. We carried out fractionation of non-cancerous bile samples using a multipronged approach (SDS-PAGE, SCX and OFFGEL) followed by MS analysis on an LTQ-Orbitrap Velos mass spectrometer using high resolution at both MS and MS/MS levels. We identified 2552 proteins - the largest number of proteins reported in human bile till date. To our knowledge, there are no previous studies employing high-resolution MS reporting a more detailed catalog of any body fluid proteome in a single study. We propose that extensive fractionation coupled to high-resolution MS can be used as a standard methodology for in-depth characterization of any body fluid. This catalog should serve as a baseline for the future studies aimed at discovering biomarkers from bile in gallbladder, hepatic, and biliary cancers.     

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.