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

Alterations in plant sugar metabolism: signatory of pathogen attack

This review summarizes the current understanding, future challenges and ongoing quest on sugar metabolic alterations that influence the outcome of plant-pathogen interactions. Intricate cellular and molecular events occur during plant-pathogen interactions. They cause major metabolic perturbations in the host and alterations in sugar metabolism play a pivotal role in governing the outcome of various kinds of plant-pathogen interactions. Sugar metabolizing enzymes and transporters of both host and pathogen origin get differentially regulated during the interactions. Both plant and pathogen compete for utilizing the host sugar metabolic machinery and in turn promote resistant or susceptible responses. However, the kind of sugar metabolism alteration that is beneficial for the host or pathogen is yet to be properly understood. Recently developed tools and methodologies are facilitating research to understand the intricate dynamics of sugar metabolism during the interactions. The present review elaborates current understanding, future challenges and ongoing quest on sugar metabolism, mobilization and regulation during various plant-pathogen interactions.

Dissection of Ramularia Leaf Spot Disease by Integrated Analysis of Barley and Ramularia collo-cygni Transcriptome Responses

Ramularia leaf spot disease (RLS), caused by the ascomycete fungus Ramularia collo-cygni, has emerged as a major economic disease of barley. No substantial resistance has been identified, so far, among barley genotypes and, based on the epidemiology of the disease, a quantitative genetic determinacy of RLS has been suggested. The relative contributions of barley and R. collo-cygni genetics to disease infection and epidemiology are practically unknown. Here, we present an integrated genome-wide analysis of host and pathogen transcriptome landscapes identified in a sensitive barley cultivar following infection by an aggressive R. collo-cygni isolate. We compared transcriptional responses in the infected and noninfected leaf samples in order to identify which molecular events are associated with RLS symptom development. We found a large proportion of R. collo-cygni genes to be expressed in planta and that many were also closely associated with the infection stage. The transition from surface to apoplastic colonization was associated with downregulation of cell wall-degrading genes and upregulation of nutrient uptake and resistance to oxidative stresses. Interestingly, the production of secondary metabolites was dynamically regulated within the fungus, indicating that R. collo-cygni produces a diverse panel of toxic compounds according to the infection stage. A defense response against R. collo-cygni was identified in barley at the early, asymptomatic infection and colonization stages. We found activation of ethylene signaling, jasmonic acid signaling, and phenylpropanoid and flavonoid pathways to be highly induced, indicative of a classical response to necrotrophic pathogens. Disease development was found to be associated with gene expression patterns similar to those found at the onset of leaf senescence, when nutrients, possibly, are used by the infecting fungus. These analyses, combining both barley and R. collo-cygni transcript profiles, demonstrate the activation of complex transcriptional programs in both organisms.

Highly flexible infection programs in a specialized wheat pathogen

Many filamentous plant pathogens exhibit high levels of genomic variability, yet the impact of this variation on host-pathogen interactions is largely unknown. We have addressed host specialization in the wheat pathogen Zymoseptoria tritici. Our study builds on comparative analyses of infection and gene expression phenotypes of three isolates and reveals the extent to which genomic variation translates into phenotypic variation. The isolates exhibit genetic and genomic variation but are similarly virulent. By combining confocal microscopy, disease monitoring, staining of ROS, and comparative transcriptome analyses, we conducted a detailed comparison of the infection processes of these isolates in a susceptible wheat cultivar. We characterized four core infection stages: establishment, biotrophic growth, lifestyle transition, and necrotrophic growth and asexual reproduction that are shared by the three isolates. However, we demonstrate differentiated temporal and spatial infection development and significant differences in the expression profiles of the three isolates during the infection stages. More than 20% of the genes were differentially expressed and these genes were located significantly closer to transposable elements, suggesting an impact of epigenetic regulation. Further, differentially expressed genes were enriched in effector candidates suggesting that isolate-specific strategies for manipulating host defenses are present in Z. tritici. We demonstrate that individuals of a host-specialized pathogen have highly differentiated infection programs characterized by flexible infection development and functional redundancy. This illustrates how high genetic diversity in pathogen populations results in highly differentiated infection phenotypes, which fact needs to be acknowledged to understand host-pathogen interactions and pathogen evolution.

Phosphopantetheinyl transferase (Ppt)-mediated biosynthesis of lysine, but not siderophores or DHN melanin, is required for virulence of Zymoseptoria tritici on wheat

Zymoseptoria tritici is the causal agent of Septoria tritici blotch (STB) disease of wheat. Z. tritici is an apoplastic fungal pathogen, which does not penetrate plant cells at any stage of infection, and has a long initial period of symptomless leaf colonisation. During this phase it is unclear to what extent the fungus can access host plant nutrients or communicate with plant cells. Several important primary and secondary metabolite pathways in fungi are regulated by the post-translational activator phosphopantetheinyl transferase (Ppt) which provides an essential co-factor for lysine biosynthesis and the activities of non-ribosomal peptide synthases (NRPS) and polyketide synthases (PKS). To investigate the relative importance of lysine biosynthesis, NRPS-based siderophore production and PKS-based DHN melanin biosynthesis, we generated deletion mutants of ZtPpt. The ∆ZtPpt strains were auxotrophic for lysine and iron, non-melanised and non-pathogenic on wheat. Deletion of the three target genes likely affected by ZtPpt loss of function (Aar- lysine; Nrps1-siderophore and Pks1- melanin), highlighted that lysine auxotrophy was the main contributing factor for loss of virulence, with no reduction caused by loss of siderophore production or melanisation. This reveals Ppt, and the lysine biosynthesis pathway, as potential targets for fungicides effective against Z. tritici.

Transcriptomic and proteomic analyses of mulberry (Morus atropurpurea) fruit response to Ciboria carunculoides

The aim of this work was to gain insights into the molecular mechanisms and dynamics of the mulberry (Morus atropurpurea) fruit response to Ciboria carunculoides infection. A transcriptomic and proteomic study was carried out based on RNA sequencing and isobaric tags for relative and absolute quantification analysis, respectively. These data were then validated using quantitative real-time PCR and multiple reaction monitoring assays. Comparative analyses revealed that 9.0% of the transcriptome and 20.8% of the proteome were differentially regulated after C. carunculoides infection at the early stage (stage 1) and middle stage (stage 2), but correlation analysis revealed that only 145 genes were differentially regulated at both the transcriptome and proteome levels. The combined transcriptome and proteome analysis showed that plant hormone signal transduction, calcium-mediated defense signaling, transcription factors, and secondary metabolites were stimulated, whereas photosynthesis and cellular growth-related metabolism were suppressed after C. carunculoides infection. These finding provide theoretical foundation for disease resistance breeding of C. carunculoides. BIOLOGICAL SIGNIFICANCE: Ciboria carunculoides is a major fungal pathogen that infects mulberry fruit, leading to extensive damage and productivity loss. Despite this major impact, the mulberry fruit response to C. carunculoides infection has yet to be characterized. This study provides the first system-wide datasets with which to examine changes in the transcriptome and proteome after C. carunculoides infection in mulberry fruit. The results showed that plant hormone signal transduction, calcium-mediated defense signaling, and other pathways were stimulated, whereas photosynthesis and cellular growth-related metabolism were suppressed by C. carunculoides. These results will lead to a better understanding of the molecular mechanisms triggered in mulberry fruit in response to C. carunculoides infection and will provide new molecular targets for regulating defense responses to fungal pathogens in berry fruits.

Extraordinary Genome Instability and Widespread Chromosome Rearrangements During Vegetative Growth

The haploid genome of the pathogenic fungus Zymoseptoria tritici is contained on "core" and "accessory" chromosomes. While 13 core chromosomes are found in all strains, as many as eight accessory chromosomes show presence/absence variation and rearrangements among field isolates. The factors influencing these presence/absence polymorphisms are so far unknown. We investigated chromosome stability using experimental evolution, karyotyping, and genome sequencing. We report extremely high and variable rates of accessory chromosome loss during mitotic propagation in vitro and in planta Spontaneous chromosome loss was observed in 2 to >50% of cells during 4 weeks of incubation. Similar rates of chromosome loss in the closely related Zymoseptoria ardabiliae suggest that this extreme chromosome dynamic is a conserved phenomenon in the genus. Elevating the incubation temperature greatly increases instability of accessory and even core chromosomes, causing severe rearrangements involving telomere fusion and chromosome breakage. Chromosome losses do not affect the fitness of Zymoseptoria tritici in vitro, but some lead to increased virulence, suggesting an adaptive role of this extraordinary chromosome instability.

CaC3H14 encoding a tandem CCCH zinc finger protein is directly targeted by CaWRKY40 and positively regulates the response of pepper to inoculation by Ralstonia solanacearum

Tandem CCCH zinc finger (TZnF) proteins have been implicated in plant defence, but their role in pepper (Capsicum annuum) is unclear. In the present study, the role of CaC3H14, a pepper TZnF protein, in the immune response of pepper plants to Ralstonia solanacearum infection was characterized. When fused to the green fluorescent protein, CaC3H14 was localized exclusively to the nuclei in leaf cells of Nicotiana benthamiana plants transiently overexpressing CaC3H14. Transcript abundance of CaC3H14 was up-regulated by inoculation with R. solanacearum. Virus-induced silencing of CaC3H14 increased the susceptibility of the plants to R. solanacearum and down-regulated the genes associated with the hypersensitive response (HR), specifically HIR1 and salicylic acid (SA)-dependent PR1a. By contrast, silencing resulted in the up-regulation of jasmonic acid (JA)-dependent DEF1 and ethylene (ET) biosynthesis-associated ACO1. Transient overexpression of CaC3H14 in pepper triggered an intensive HR, indicated by cell death and hydrogen peroxide (H₂ O₂ ) accumulation, up-regulated PR1a and down-regulated DEF1 and ACO1. Ectopic overexpression of CaC3H14 in tobacco plants significantly decreased the susceptibility of tobacco plants to R. solanacearum. It also up-regulated HR-associated HSR515, immunity-associated GST1 and the SA-dependent marker genes NPR1 and PR2, but down-regulated JA-dependent PR1b and ET-dependent EFE26. The CaC3H14 promoter and was bound and its transcription was up-regulated by CaWRKY40. Collectively, these results indicate that CaC3H14 is transcriptionally targeted by CaWRKY40, is a modulator of the antagonistic interaction between SA and JA/ET signalling, and enhances the defence response of pepper plants to infection by R. solanacearum.

Comparative proteomic analysis of maize (Zea mays L.) seedlings under rice black-streaked dwarf virus infection

BACKGROUND

Maize rough dwarf disease (MRDD) is a severe disease that has been occurring frequently in southern China and many other Asian countries. MRDD is caused by the infection of Rice black streaked dwarf virus (RBSDV) and leads to significant economic losses in maize production. To well understand the destructive effects of RBSDV infection on maize growth, comparative proteomic analyses of maize seedlings under RBSDV infection was performed using an integrated approach involving LC-MS/MS and Tandem Mass Tag (TMT) labeling.

RESULTS

In total, 7615 maize proteins, 6319 of which were quantified. A total of 116 differentially accumulated proteins (DAPs) were identified, including 35 up- and 81 down-regulated proteins under the RBSDV infection. Enrichment analysis showed that the DAPs were most strongly associated with cyanoamino acid metabolism, protein processing in ER, and ribosome-related pathways. Two sulfur metabolism-related proteins were significantly reduced, indicating that sulfur may participate in the resistance against RBSDV infection. Furthermore, 15 DAPs involved in six metabolic pathways were identified in maize under the RBSDV infection.

CONCLUSIONS

Our data revealed that the responses of maize to RBSDV infection were controlled by various metabolic pathways.

Genome-Based Discovery of Polyketide-Derived Secondary Metabolism Pathways in the Barley Pathogen Ramularia collo-cygni

Ramularia collo-cygni causes Ramularia leaf spot (RLS) disease of barley. The fungus develops asymptomatically within its host until late in the growing season, when necrotic lesions become visible on upper leaves. Fungal secondary metabolites (SM) have been proposed as important factors in RLS lesion formation but the biosynthetic pathways involved remain largely unknown. Mining the R. collo-cygni genome revealed the presence of 10 polyketide synthases (PKS), 10 nonribosomal peptide synthetases (NRPS), and 3 hybrid PKS-NRPS (HPS) identified within clusters of genes with predicted functions associated with secondary metabolism. SM core genes along with their predicted transcriptional regulators exhibited transcriptional coexpression during infection of barley plants. Moreover, their expression peaked during early stages of host colonization and preceded or overlapped with the appearance of disease symptoms, suggesting that SM may manipulate the host to promote colonization or protect R. collo-cygni from competing organisms. Accordingly, R. collo-cygni inhibited the growth of several fungi in vitro, indicating that it synthesized and excreted antifungal agents. Taken together, these findings demonstrate that the R. collo-cygni genome contains the genetic architecture to synthesize a wide range of SM and suggests that coexpression of PKS and HPS is associated with competitive colonization of the host and early symptom development.

The Genome Biology of Effector Gene Evolution in Filamentous Plant Pathogens

Filamentous pathogens, including fungi and oomycetes, pose major threats to global food security. Crop pathogens cause damage by secreting effectors that manipulate the host to the pathogen's advantage. Genes encoding such effectors are among the most rapidly evolving genes in pathogen genomes. Here, we review how the major characteristics of the emergence, function, and regulation of effector genes are tightly linked to the genomic compartments where these genes are located in pathogen genomes. The presence of repetitive elements in these compartments is associated with elevated rates of point mutations and sequence rearrangements with a major impact on effector diversification. The expression of many effectors converges on an epigenetic control mediated by the presence of repetitive elements. Population genomics analyses showed that rapidly evolving pathogens show high rates of turnover at effector loci and display a mosaic in effector presence-absence polymorphism among strains. We conclude that effective pathogen containment strategies require a thorough understanding of the effector genome biology and the pathogen's potential for rapid adaptation.

Exploring the Genetic Regulation of Asexual Sporulation in Zymoseptoria tritici

Zymoseptoria tritici is the causal agent of septoria tritici blotch, a devastating fungal disease of wheat which can cause up to 40% yield loss. One of the ways in which Z. tritici spreads in the field is via rain splash-dispersed asexual pycnidiospores, however there is currently limited understanding of the genetic mechanisms governing the development of these propagules. In order to explore whether the existing models for conidiation in ascomycete fungi apply to Z. tritici, homologs to the well-characterized Aspergillus nidulans genes abacus (abaA), bristle (brlA), fluffy B (flbB), fluffy C (flbC), and stunted (stuA) were identified and knocked-out by Agrobacterium-mediated transformation. Although deletion of the ZtAbaA, ZtBrlA1, and ZtFlbB genes had no apparent effect on Z. tritici asexual sporulation or on pathogenicity, deletion of ZtFlbC or ZtBrlA2 resulted in mutants with reduced pycnidiospore production compared to the parental IPO323 strain. Deletion of ZtStuA gave non-pigmented mutants with altered vegetative growth and eliminated asexual sporulation and pathogenicity. These findings suggest that the well-established A. nidulans model of asexual sporulation is only partially applicable to Z. tritici, and that this pathogen likely uses additional, as yet uncharacterized genes to control asexual sporulation.

Exploring the biological roles of Dothideomycetes ABC proteins: Leads from their phylogenetic relationships with functionally-characterized Ascomycetes homologs

BACKGROUND

The ATP-binding cassette (ABC) superfamily is one of the largest, ubiquitous and diverse protein families in nature. Categorized into nine subfamilies, its members are important to most organisms including fungi, where they play varied roles in fundamental cellular processes, plant pathogenesis or fungicide tolerance. However, these proteins are not yet well-understood in the class Dothideomycetes, which includes several phytopathogens that infect a wide range of food crops including wheat, barley and maize and cause major economic losses.

RESULTS

We analyzed the genomes of 14 Dothideomycetes fungi (Test set) and seven well-known Ascomycetes fungi (Model set- that possessed gene expression/ functional analysis data about the ABC genes) and predicted 578 and 338 ABC proteins from each set respectively. These proteins were classified into subfamilies A to I, which revealed the distribution of the subfamily members across the Dothideomycetes and Ascomycetes genomes. Phylogenetic analysis of Dothideomycetes ABC proteins indicated evolutionary relationships among the subfamilies within this class. Further, phylogenetic relationships among the ABC proteins from the Model and the Test fungi within each subfamily were analyzed, which aided in classifying these proteins into subgroups. We compiled and curated functional and gene expression information from the previous literature for 118 ABC genes and mapped them on the phylogenetic trees, which suggested possible roles in pathogenesis and/or fungicide tolerance for the newly identified Dothideomycetes ABC proteins.

CONCLUSIONS

The present analysis is one of the firsts to extensively analyze ABC proteins from Dothideomycetes fungi. Their phylogenetic analysis and annotating the clades with functional information indicated a subset of Dothideomycetes ABC genes that could be considered for experimental validation for their roles in plant pathogenesis and/or fungicide tolerance.

A fungal avirulence factor encoded in a highly plastic genomic region triggers partial resistance to septoria tritici blotch

Cultivar-strain specificity in the wheat-Zymoseptoria tritici pathosystem determines the infection outcome and is controlled by resistance genes on the host side, many of which have been identified. On the pathogen side, however, the molecular determinants of specificity remain largely unknown. We used genetic mapping, targeted gene disruption and allele swapping to characterise the recognition of the new avirulence factor Avr3D1. We then combined population genetic and comparative genomic analyses to characterise the evolutionary trajectory of Avr3D1. Avr3D1 is specifically recognised by wheat cultivars harbouring the Stb7 resistance gene, triggering a strong defence response without preventing pathogen infection and reproduction. Avr3D1 resides in a cluster of putative effector genes located in a genome region populated by independent transposable element insertions. The gene was present in all 132 investigated strains and is highly polymorphic, with 30 different protein variants identified. We demonstrated that specific amino acid substitutions in Avr3D1 led to evasion of recognition. These results demonstrate that quantitative resistance and gene-for-gene interactions are not mutually exclusive. Localising avirulence genes in highly plastic genomic regions probably facilitates accelerated evolution that enables escape from recognition by resistance proteins.

Beyond the one-way ANOVA for 'omics data

BACKGROUND

With ever increasing accessibility to high throughput technologies, more complex treatment structures can be assessed in a variety of 'omics applications. This adds an extra layer of complexity to the analysis and interpretation, in particular when inferential univariate methods are applied en masse. It is well-known that mass univariate testing suffers from multiplicity issues and although this has been well documented for simple comparative tests, few approaches have focussed on more complex explanatory structures.

RESULTS

Two frameworks are introduced incorporating corrections for multiplicity whilst maintaining appropriate structure in the explanatory variables. Within this paradigm, a choice has to be made as to whether multiplicity corrections should be applied to the saturated model, putting emphasis on controlling the rate of false positives, or to the predictive model, where emphasis is on model selection. This choice has implications for both the ranking and selection of the response variables identified as differentially expressed. The theoretical difference is demonstrated between the two approaches along with an empirical study of lipid composition in Arabidopsis under differing levels of salt stress.

CONCLUSIONS

Multiplicity corrections have an inherent weakness when the full explanatory structure is not properly incorporated. Although a unifying 'single best' recommendation is not provided, two reasonable alternatives are provided and the applicability of these approaches is discussed for different scenarios where the aims of analysis will differ. The key result is that the point at which multiplicity is incorporated into the analysis will fundamentally change the interpretation of the results, and the choice of approach should therefore be driven by the specific aims of the experiment.

Origin and Spread of Spliceosomal Introns: Insights from the Fungal Clade Zymoseptoria

Spliceosomal introns are a key feature of eukaryote genome architecture and have been proposed to originate from selfish group II introns from an endosymbiotic bacterium, that is, the ancestor of mitochondria. However, the mechanisms underlying the wide spread of spliceosomal introns across eukaryotic genomes have been obscure. In this study, we characterize the dynamic evolution of spliceosomal introns in the fungal genus Zymoseptoria at different evolutionary scales, that is, within a genome, among conspecific strains within species, and between different species. Within the genome, spliceosomal introns can proliferate in unrelated genes and intergenic regions. Among conspecific strains, spliceosomal introns undergo rapid turnover (gains and losses) and frequent sequence exchange between geographically distinct strains. Furthermore, spliceosomal introns could undergo introgression between distinct species, which can further promote intron invasion and proliferation. The dynamic invasion and proliferation processes of spliceosomal introns resemble the life cycles of mobile selfish (group I/II) introns, and these intron movements, at least in part, account for the dramatic processes of intron gain and intron loss during eukaryotic evolution.

Proximal Phenotyping and Machine Learning Methods to Identify Septoria Tritici Blotch Disease Symptoms in Wheat

Phenotyping with proximal sensors allow high-precision measurements of plant traits both in the controlled conditions and in the field. In this work, using machine learning, an integrated analysis was done from the data obtained from spectroradiometer, infrared thermometer, and chlorophyll fluorescence measurements to identify most predictive proxy measurements for studying Septoria tritici blotch (STB) disease of wheat. The random forest (RF) models for chlorosis and necrosis identified photosystem II quantum yield (QY) and vegetative indices (VIs) associated with the biochemical composition of leaves as the top predictive variables for identifying disease symptoms. The RF model for chlorosis was validated with a validation set (R2: 0.80) and in an independent test set (R2: 0.55). Based on the results, it can be concluded that the proxy measurements for photosystem II, chlorophyll content, carotenoid, and anthocyanin levels and leaf surface temperature can be successfully used to detect STB. Further validation of these results in the field will enable application of these predictive variables for detection of STB in the field.

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.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

Novel Fungal Pathogenicity and Leaf Defense Strategies Are Revealed by Simultaneous Transcriptome Analysis of Colletotrichum fructicola and Strawberry Infected by This Fungus

Colletotrichum fructicola, which is part of the C. gloeosporioides species complex, can cause anthracnose diseases in strawberries worldwide. However, the molecular interactions between C. fructicola and strawberry are largely unknown. A deep RNA-sequencing approach was applied to gain insights into the pathogenicity mechanisms of C. fructicola and the defense response of strawberry plants at different stages of infection. The transcriptome data showed stage-specific transcription accompanied by a step-by-step strawberry defense response and the evasion of this defense system by fungus. Fungal genes involved in plant cell wall degradation, secondary metabolism, and detoxification were up-regulated at different stage of infection. Most importantly, C. fructicola infection was accompanied by a large number of highly expressed effectors. Four new identified effectors function in the suppression of Bax-mediated programmed cell death. Strawberry utilizes pathogen-associated molecular patterns (PAMP)-triggered immunity and effector-triggered immunity to prevent C. fructicola invasion, followed by the initiation of downstream innate immunity. The up-regulation of genes related to salicylic acid provided evidence that salicylic acid signaling may serve as the core defense signaling mechanism, while jasmonic acid and ethylene pathways were largely inhibited by C. fructicola. The necrotrophic stage displayed a significant up-regulation of genes involved in reactive oxygen species activation. Collectively, the transcriptomic data of both C. fructicola and strawberry shows that even though plants build a multilayered defense against infection, C. fructicola employs a series of escape or antagonizing mechanisms to successfully infect host cells.

A Review of the Interactions between Wheat and Wheat Pathogens: Zymoseptoria tritici, Fusarium spp. and Parastagonospora nodorum

Zymoseptoria tritici is a hemibiotrophic pathogen which causes Septoria leaf blotch in wheat. The pathogenesis of the disease consists of a biotrophic phase and a necrotrophic phase. The pathogen infects the host plant by suppressing its immune response in the first stage of infection. Hemibiotrophic pathogens of the genus Fusarium cause Fusarium head blight, and the necrotrophic Parastagonosporanodorum is responsible for Septoria nodorum blotch in wheat. Cell wall-degrading enzymes in plants promote infections by necrotrophic and hemibiotrophic pathogens, and trichothecenes, secondary fungal metabolites, facilitate infections caused by fungi of the genus Fusarium. There are no sources of complete resistance to the above pathogens in wheat. Defense mechanisms in wheat are controlled by many genes encoding resistance traits. In the wheat genome, the characteristic features of loci responsible for resistance to pathogenic infections indicate that at least several dozen genes encode resistance to pathogens. The molecular interactions between wheat and Z. tritici, P. nodorum and Fusarium spp. pathogens have been insufficiently investigated. Most studies focus on the mechanisms by which the hemibiotrophic Z. tritici suppresses immune responses in plants and the role of mycotoxins and effector proteins in infections caused by P. nodorum and Fusarium spp. fungi. Trichothecene glycosylation and effector proteins, which are involved in defense responses in wheat, have been described at the molecular level. Recent advances in molecular biology have produced interesting findings which should be further elucidated in studies of molecular interactions between wheat and fungal pathogens. The Clustered Regularly-Interspaced Short Palindromic Repeats/ CRISPR associated (CRISPR/Cas) system can be used to introduce targeted mutations into the wheat genome and confer resistance to selected fungal diseases. Host-induced gene silencing and spray-induced gene silencing are also useful tools for analyzing wheat-pathogens interactions which can be used to develop new strategies for controlling fungal diseases.

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.

Distinct Trajectories of Massive Recent Gene Gains and Losses in Populations of a Microbial Eukaryotic Pathogen

Differences in gene content are a significant source of variability within species and have an impact on phenotypic traits. However, little is known about the mechanisms responsible for the most recent gene gains and losses. We screened the genomes of 123 worldwide isolates of the major pathogen of wheat Zymoseptoria tritici for robust evidence of gene copy number variation. Based on orthology relationships in three closely related fungi, we identified 599 gene gains and 1,024 gene losses that have not yet reached fixation within the focal species. Our analyses of gene gains and losses segregating in populations showed that gene copy number variation arose preferentially in subtelomeres and in proximity to transposable elements. Recently lost genes were enriched in virulence factors and secondary metabolite gene clusters. In contrast, recently gained genes encoded mostly secreted protein lacking a conserved domain. We analyzed the frequency spectrum at loci segregating a gene presence–absence polymorphism in four worldwide populations. Recent gene losses showed a significant excess in low-frequency variants compared with genome-wide single nucleotide polymorphism, which is indicative of strong negative selection against gene losses. Recent gene gains were either under weak negative selection or neutral. We found evidence for strong divergent selection among populations at individual loci segregating a gene presence–absence polymorphism. Hence, gene gains and losses likely contributed to local adaptation. Our study shows that microbial eukaryotes harbor extensive copy number variation within populations and that functional differences among recently gained and lost genes led to distinct evolutionary trajectories.

The role of reactive oxygen in the development of Ramularia leaf spot disease in barley seedlings

BACKGROUND AND AIMS

Ramularia collo-cygni is an ascomycete fungus that colonizes barley primarily as a benign endophyte, although this interaction can become pathogenic, causing the disease Ramularia leaf spot (RLS). Factors, particularly reactive oxygen species, that resulted in the transition of the fungus from endophyte to necrotrophic parasite and the development of disease symptoms were investigated.

METHODS

Disease development in artificially inoculated seedlings of barley varieties varying in partial resistance to RLS was related to exposure to abiotic stress prior to inoculation. Histochemical and molecular analysis determined the effect of R. collo-cygni colonization on accumulation of reactive oxygen species and antioxidant gene expression. Development of RLS on barley lines defective in antioxidant enzymes and with altered redox status or non-functional chloroplasts was compared with the accumulation of fungal biomass to determine how these factors affect disease symptom expression.

KEY RESULTS

Exposure to abiotic stress increased symptom development in all susceptible and most partially resistant barley varieties, in association with greater hydrogen peroxide (H2O2) levels in leaves. Decreased activity of the antioxidant enzymes superoxide dismutase and catalase in transgenic and mutant plants had no effect on the disease transition, whereas manipulation of H2O2 levels during asymptomatic growth of the fungus increased disease symptoms in most susceptible varieties but not in partially resistant plants. Barley mutants that undergo rapid loss of green leaf area when infected by R. collo-cygni or albino mutants with non-functional chloroplasts showed reduced development of RLS symptoms.

CONCLUSIONS

These results imply that in seedlings the pathogenic transition of the normally endophytic fungus R. collo-cygni does not result from senescence as such, but rather is promoted by factors that result in changes to host reactive oxygen species. Barley varieties vary in the extent to which these factors promote RLS disease.

Stress and sexual reproduction affect the dynamics of the wheat pathogen effector AvrStb6 and strobilurin resistance

Host resistance and fungicide treatments are cornerstones of plant-disease control. Here, we show that these treatments allow sex and modulate parenthood in the fungal wheat pathogen Zymoseptoria tritici. We demonstrate that the Z. tritici-wheat interaction complies with the gene-for-gene model by identifying the effector AvrStb6, which is recognized by the wheat resistance protein Stb6. Recognition triggers host resistance, thus implying removal of avirulent strains from pathogen populations. However, Z. tritici crosses on wheat show that sex occurs even with an avirulent parent, and avirulence alleles are thereby retained in subsequent populations. Crossing fungicide-sensitive and fungicide-resistant isolates under fungicide pressure results in a rapid increase in resistance-allele frequency. Isolates under selection always act as male donors, and thus disease control modulates parenthood. Modeling these observations for agricultural and natural environments reveals extended durability of host resistance and rapid emergence of fungicide resistance. Therefore, fungal sex has major implications for disease control.

The Biotrophic Development of Ustilago maydis Studied by RNA-Seq Analysis

The maize smut fungus Ustilago maydis is a model organism for elucidating host colonization strategies of biotrophic fungi. Here, we performed an in depth transcriptional profiling of the entire plant-associated development of U. maydis wild-type strains. In our analysis, we focused on fungal metabolism, nutritional strategies, secreted effectors, and regulatory networks. Secreted proteins were enriched in three distinct expression modules corresponding to stages on the plant surface, establishment of biotrophy, and induction of tumors. These modules are likely the key determinants for U. maydis virulence. With respect to nutrient utilization, we observed that expression of several nutrient transporters was tied to these virulence modules rather than being controlled by nutrient availability. We show that oligopeptide transporters likely involved in nitrogen assimilation are important virulence factors. By measuring the intramodular connectivity of transcription factors, we identified the potential drivers for the virulence modules. While known components of the b-mating type cascade emerged as inducers for the plant surface and biotrophy module, we identified a set of yet uncharacterized transcription factors as likely responsible for expression of the tumor module. We demonstrate a crucial role for leaf tumor formation and effector gene expression for one of these transcription factors.

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.

Characterization of an antimicrobial and phytotoxic ribonuclease secreted by the fungal wheat pathogen Zymoseptoria tritici

The fungus Zymoseptoria tritici is the causal agent of Septoria Tritici Blotch (STB) disease of wheat leaves. Zymoseptoria tritici secretes many functionally uncharacterized effector proteins during infection. Here, we characterized a secreted ribonuclease (Zt6) with an unusual biphasic expression pattern. Transient expression systems were used to characterize Zt6, and mutants thereof, in both host and non-host plants. Cell-free protein expression systems monitored the impact of Zt6 protein on functional ribosomes, and in vitro assays of cells treated with recombinant Zt6 determined toxicity against bacteria, yeasts and filamentous fungi. We demonstrated that Zt6 is a functional ribonuclease and that phytotoxicity is dependent on both the presence of a 22-amino-acid N-terminal 'loop' region and its catalytic activity. Zt6 selectively cleaves both plant and animal rRNA species, and is toxic to wheat, tobacco, bacterial and yeast cells, but not to Z. tritici itself. Zt6 is the first Z. tritici effector demonstrated to have a likely dual functionality. The expression pattern of Zt6 and potent toxicity towards microorganisms suggest that, although it may contribute to the execution of wheat cell death, it is also likely to have an important secondary function in antimicrobial competition and niche protection.

The Importance of Experimental Design, Quality Assurance, and Control in Plant Metabolomics Experiments

The output of metabolomics relies to a great extent upon the methods and instrumentation to identify, quantify, and access spatial information on as many metabolites as possible. However, the most modern machines and sophisticated tools for data analysis cannot compensate for inappropriate harvesting and/or sample preparation procedures that modify metabolic composition and can lead to erroneous interpretation of results. In addition, plant metabolism has a remarkable degree of complexity, and the number of identified compounds easily surpasses the number of samples in metabolomics analyses, increasing false discovery risk. These aspects pose a large challenge when carrying out plant metabolomics experiments. In this chapter, we address the importance of a proper experimental design taking into consideration preventable complications and unavoidable factors to achieve success in metabolomics analysis. We also focus on quality control and standardized procedures during the metabolomics workflow.

Non-parent of Origin Expression of Numerous Effector Genes Indicates a Role of Gene Regulation in Host Adaption of the Hybrid Triticale Powdery Mildew Pathogen

Powdery mildew is an important disease of cereals. It is caused by one species, Blumeria graminis, which is divided into formae speciales each of which is highly specialized to one host. Recently, a new form capable of growing on triticale (B.g. triticale) has emerged through hybridization between wheat and rye mildews (B.g. tritici and B.g. secalis, respectively). In this work, we used RNA sequencing to study the molecular basis of host adaptation in B.g. triticale. We analyzed gene expression in three B.g. tritici isolates, two B.g. secalis isolates and two B.g. triticale isolates and identified a core set of putative effector genes that are highly expressed in all formae speciales. We also found that the genes differentially expressed between isolates of the same form as well as between different formae speciales were enriched in putative effectors. Their coding genes belong to several families including some which contain known members of mildew avirulence (Avr) and suppressor (Svr) genes. Based on these findings we propose that effectors play an important role in host adaptation that is mechanistically based on Avr-Resistance gene-Svr interactions. We also found that gene expression in the B.g. triticale hybrid is mostly conserved with the parent-of-origin, but some genes inherited from B.g. tritici showed a B.g. secalis-like expression. Finally, we identified 11 unambiguous cases of putative effector genes with hybrid-specific, non-parent of origin gene expression, and we propose that they are possible determinants of host specialization in triticale mildew. These data suggest that altered expression of multiple effector genes, in particular Avr and Svr related factors, might play a role in mildew host adaptation based on hybridization.

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.

Biotrophy-necrotrophy switch in pathogen evoke differential response in resistant and susceptible sesame involving multiple signaling pathways at different phases

Infection stages of charcoal rot fungus Macrophomina phaseolina in sesame revealed for the first time a transition from biotrophy via BNS (biotrophy-to-necrotrophy switch) to necrotrophy as confirmed by transcriptional studies. Microscopy using normal and GFP-expressing pathogen showed typical constricted thick intercellular bitrophic hyphae which gave rise to thin intracellular necrotrophic hyphae during BNS and this stage was delayed in a resistant host. Results also show that as the pathogen switched its strategy of infection, the host tailored its defense strategy to meet the changing situation. Less ROS accumulation, upregulation of ROS signaling genes and higher antioxidant enzyme activities post BNS resulted in resistance. There was greater accumulation of secondary metabolites and upregulation of secondary metabolite-related genes after BNS. A total of twenty genes functioning in different aspects of plant defense that were monitored over a time course during the changing infection phases showed a coordinated response. Experiments using phytohormone priming and phytohormone inhibitors showed that resistance resulted from activation of JA-ET signaling pathway. Most importantly this defense response was more prompt in the resistant than the susceptible host indicating that a resistant host makes different choices from a susceptible host during infection which ultimately influences the severity of the disease.

Forward Genetics Approach Reveals Host Genotype-Dependent Importance of Accessory Chromosomes in the Fungal Wheat Pathogen Zymoseptoria tritici

The fungal wheat pathogen Zymoseptoria tritici possesses a large complement of accessory chromosomes showing presence/absence polymorphism among isolates. These chromosomes encode hundreds of genes; however, their functional role and why the chromosomes have been maintained over long evolutionary times are so far not known. In this study, we addressed the functional relevance of eight accessory chromosomes in reference isolate IPO323. We induced chromosome losses by inhibiting the β-tubulin assembly during mitosis using carbendazim and generated several independent isogenic strains, each lacking one of the accessory chromosomes. We confirmed chromosome losses by electrophoretic karyotyping and whole-genome sequencing. To assess the importance of the individual chromosomes during host infection, we performed in planta assays comparing disease development results in wild-type and chromosome mutant strains. Loss of the accessory chromosomes 14, 16, 18, 19, and 21 resulted in increased virulence on wheat cultivar Runal but not on cultivars Obelisk, Titlis, and Riband. Moreover, some accessory chromosomes affected the switch from biotrophy to necrotrophy as strains lacking accessory chromosomes 14, 18, 19, and 21 showed a significantly earlier onset of necrosis than the wild type on the Runal cultivar. In general, we observed that the timing of the lifestyle switch affects the fitness of Z. tritici. Taking the results together, this study was the first to use a forward-genetics approach to demonstrate a cultivar-dependent functional relevance of the accessory chromosomes of Z. tritici during host infection.

Zymoseptoria tritici is a prominent fungal pathogen of wheat of worldwide distribution. This fungus shows a remarkable genome organization, with a large number of chromosomes that are present in only some isolates and therefore considered to be “accessory” chromosomes. To date, the function of these accessory chromosomes in Z. tritici has been unknown, although their maintenance in the species over evolutionary times suggests a functional relevance. Here we deleted whole accessory chromosomes to test the effect of these chromosomes on host specificity and virulence of the fungus. We show for the first time that some accessory chromosomes of Z. tritici affect the fitness of the fungus during host infection in a cultivar-dependent manner. These results show that the accessory chromosomes encode host-specific virulence determinants having a negative effect on fitness. Understanding the population dynamic of the accessory chromosomes and the molecular interaction of pathogen and plant traits is crucial to improve wheat-breeding strategies.

A conserved fungal glycosyltransferase facilitates pathogenesis of plants by enabling hyphal growth on solid surfaces

Pathogenic fungi must extend filamentous hyphae across solid surfaces to cause diseases of plants. However, the full inventory of genes which support this is incomplete and many may be currently concealed due to their essentiality for the hyphal growth form. During a random T-DNA mutagenesis screen performed on the pleomorphic wheat (Triticum aestivum) pathogen Zymoseptoria tritici, we acquired a mutant unable to extend hyphae specifically when on solid surfaces. In contrast "yeast-like" growth, and all other growth forms, were unaffected. The inability to extend surface hyphae resulted in a complete loss of virulence on plants. The affected gene encoded a predicted type 2 glycosyltransferase (ZtGT2). Analysis of >800 genomes from taxonomically diverse fungi highlighted a generally widespread, but discontinuous, distribution of ZtGT2 orthologues, and a complete absence of any similar proteins in non-filamentous ascomycete yeasts. Deletion mutants of the ZtGT2 orthologue in the taxonomically un-related fungus Fusarium graminearum were also severely impaired in hyphal growth and non-pathogenic on wheat ears. ZtGT2 expression increased during filamentous growth and electron microscopy on deletion mutants (ΔZtGT2) suggested the protein functions to maintain the outermost surface of the fungal cell wall. Despite this, adhesion to leaf surfaces was unaffected in ΔZtGT2 mutants and global RNAseq-based gene expression profiling highlighted that surface-sensing and protein secretion was also largely unaffected. However, ΔZtGT2 mutants constitutively overexpressed several transmembrane and secreted proteins, including an important LysM-domain chitin-binding virulence effector, Zt3LysM. ZtGT2 likely functions in the synthesis of a currently unknown, potentially minor but widespread, extracellular or outer cell wall polysaccharide which plays a key role in facilitating many interactions between plants and fungi by enabling hyphal growth on solid matrices.

Lifestyle, gene gain and loss, and transcriptional remodeling cause divergence in the transcriptomes of Phytophthora infestans and Pythium ultimum during potato tuber colonization

BACKGROUND

How pathogen genomes evolve to support distinct lifestyles is not well-understood. The oomycete Phytophthora infestans, the potato blight agent, is a largely biotrophic pathogen that feeds from living host cells, which become necrotic only late in infection. The related oomycete Pythium ultimum grows saprophytically in soil and as a necrotroph in plants, causing massive tissue destruction. To learn what distinguishes their lifestyles, we compared their gene contents and expression patterns in media and a shared host, potato tuber.

RESULTS

Genes related to pathogenesis varied in temporal expression pattern, mRNA level, and family size between the species. A family's aggregate expression during infection was not proportional to size due to transcriptional remodeling and pseudogenization. Ph. infestans had more stage-specific genes, while Py. ultimum tended towards more constitutive expression. Ph. infestans expressed more genes encoding secreted cell wall-degrading enzymes, but other categories such as secreted proteases and ABC transporters had higher transcript levels in Py. ultimum. Species-specific genes were identified including new Pythium genes, perforins, which may disrupt plant membranes. Genome-wide ortholog analyses identified substantial diversified expression, which correlated with sequence divergence. Pseudogenization was associated with gene family expansion, especially in gene clusters.

CONCLUSION

This first large-scale analysis of transcriptional divergence within oomycetes revealed major shifts in genome composition and expression, including subfunctionalization within gene families. Biotrophy and necrotrophy seem determined by species-specific genes and the varied expression of shared pathogenicity factors, which may be useful targets for crop protection.

Omics-based hybrid prediction in maize

Complementing genomic data with other "omics" predictors can increase the probability of success for predicting the best hybrid combinations using complex agronomic traits. Accurate prediction of traits with complex genetic architecture is crucial for selecting superior candidates in animal and plant breeding and for guiding decisions in personalized medicine. Whole-genome prediction has revolutionized these areas but has inherent limitations in incorporating intricate epistatic interactions. Downstream "omics" data are expected to integrate interactions within and between different biological strata and provide the opportunity to improve trait prediction. Yet, predicting traits from parents to progeny has not been addressed by a combination of "omics" data. Here, we evaluate several "omics" predictors-genomic, transcriptomic and metabolic data-measured on parent lines at early developmental stages and demonstrate that the integration of transcriptomic with genomic data leads to higher success rates in the correct prediction of untested hybrid combinations in maize. Despite the high predictive ability of genomic data, transcriptomic data alone outperformed them and other predictors for the most complex heterotic trait, dry matter yield. An eQTL analysis revealed that transcriptomic data integrate genomic information from both, adjacent and distant sites relative to the expressed genes. Together, these findings suggest that downstream predictors capture physiological epistasis that is transmitted from parents to their hybrid offspring. We conclude that the use of downstream "omics" data in prediction can exploit important information beyond structural genomics for leveraging the efficiency of hybrid breeding.

A role for random, humidity-dependent epiphytic growth prior to invasion of wheat by Zymoseptoria tritici

Zymoseptoria tritici causes Septoria leaf blotch of wheat. The prevailing paradigm of the Z. tritici-wheat interaction assumes fungal ingress through stomata within 24-48h, followed by days of symptomless infection. This is extrapolated from studies testing the mode of fungal ingress under optimal infection conditions. Here, we explicitly assess the timing of entry, using GFP-tagged Z. tritici. We show that early entry is comparatively rare, and extended epiphytic growth possible. We test the hypotheses that our data diverge from earlier studies due to: i. random ingress of Z. tritici into the leaf, with some early entry events; ii. previous reliance upon fungal stains, combined with poor attachment of Z. tritici to the leaf, leading to increased likelihood of observing internal versus external growth, compared to using GFP; iii. use of exceptionally high humidity to promote entry in previous studies. We combine computer simulation of leaf-surface growth with thousands of in planta observations to demonstrate that while spores germinate rapidly on the leaf, over 95% of fungi remain epiphytic, growing randomly over the leaf for ten days or more. We show that epiphytic fungi are easily detached from leaves by rinsing and that humidity promotes epiphytic growth, increasing infection rates. Together, these results explain why epiphytic growth has been dismissed and early ingress assumed. The prolonged epiphytic phase should inform studies of pathogenicity and virulence mutants, disease control strategies, and interpretation of the observed low in planta growth, metabolic quiescence and evasion of plant defences by Zymoseptoria during symptomless infection.

In silico prediction and characterization of secondary metabolite biosynthetic gene clusters in the wheat pathogen Zymoseptoria tritici

BACKGROUND

Fungal pathogens of plants produce diverse repertoires of secondary metabolites, which have functions ranging from iron acquisition, defense against immune perturbation, to toxic assaults on the host. The wheat pathogen Zymoseptoria tritici causes Septoria tritici blotch, a foliar disease which is a significant threat to global food security. Currently, there is limited knowledge of the secondary metabolite arsenal produced by Z. tritici, which significantly restricts mechanistic understanding of infection. In this study, we analyzed the genome of Z. tritici isolate IP0323 to identify putative secondary metabolite biosynthetic gene clusters, and used comparative genomics to predict their encoded products.

RESULTS

We identified 32 putative secondary metabolite clusters. These were physically enriched at subtelomeric regions, which may facilitate diversification of cognate products by rapid gene rearrangement or mutations. Comparative genomics revealed a four gene cluster with significant similarity to the ferrichrome-A biosynthetic locus of the maize pathogen Ustilago maydis, suggesting this siderophore is deployed by Z. tritici to acquire iron. The Z. tritici genome also contains several isoprenoid biosynthetic gene clusters, including one with high similarity to a carotenoid/opsin producing locus in several fungi. Furthermore, we identify putative phytotoxin biosynthetic clusters, suggesting Z. tritici can produce an epipolythiodioxopiperazine, and a polyketide and non-ribosomal peptide with predicted structural similarities to fumonisin and the Alternaria alternata AM-toxin, respectively. Interrogation of an existing transcriptional dataset suggests stage specific deployment of numerous predicted loci during infection, indicating an important role of these secondary metabolites in Z. tritici disease.

CONCLUSIONS

We were able to assign putative biosynthetic products to numerous clusters based on conservation amongst other fungi. However, analysis of the majority of secondary metabolite loci did not enable prediction of a cluster product, and consequently the capacity of these loci to play as yet undetermined roles in disease or other stages of the Z. tritici lifecycle is significant. These data will drive future experimentation for determining the role of these clusters and cognate secondary metabolite products in Z. tritici virulence, and may lead to discovery of novel bioactive molecules.

Adaptation to the Host Environment by Plant-Pathogenic Fungi

Many fungi can live both saprophytically and as endophyte or pathogen inside a living plant. In both environments, complex organic polymers are used as sources of nutrients. Propagation inside a living host also requires the ability to respond to immune responses of the host. We review current knowledge of how plant-pathogenic fungi do this. First, we look at how fungi change their global gene expression upon recognition of the host environment, leading to secretion of effectors, enzymes, and secondary metabolites; changes in metabolism; and defense against toxic compounds. Second, we look at what is known about the various cues that enable fungi to sense the presence of living plant cells. Finally, we review literature on transcription factors that participate in gene expression in planta or are suspected to be involved in that process because they are required for the ability to cause disease.

A fungal wheat pathogen evolved host specialization by extensive chromosomal rearrangements

Fungal pathogens can rapidly evolve virulence towards resistant crops in agricultural ecosystems. Gains in virulence are often mediated by the mutation or deletion of a gene encoding a protein recognized by the plant immune system. However, the loci and the mechanisms of genome evolution enabling rapid virulence evolution are poorly understood. We performed genome-wide association mapping on a global collection of 106 strains of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe, to identify polymorphisms linked to virulence on two wheat varieties. We found 25 distinct genomic loci associated with reproductive success of the pathogen. However, no locus was shared between the host genotypes, suggesting host specialization. The main locus associated with virulence encoded a highly expressed, small secreted protein. Population genomic analyses showed that the gain in virulence was explained by a segregating gene deletion polymorphism. The deletion was likely adaptive by preventing detection of the encoded protein. Comparative genomics of closely related species showed that the locus emerged de novo since speciation. A large cluster of transposable elements in direct proximity to the locus generated extensive rearrangements leading to multiple independent gene losses. Our study demonstrates that rapid turnover in the chromosomal structure of a pathogen can drive host specialization.

An integrated RNAseq-1H NMR metabolomics approach to understand soybean primary metabolism regulation in response to Rhizoctonia foliar blight disease

BACKGROUND

Rhizoctonia solani AG1-IA is a devastating phytopathogen causing Rhizoctonia foliar blight (RFB) of soybean worldwide with yield losses reaching 60%. Plant defense mechanisms are complex and information from different metabolic pathways is required to thoroughly understand plant defense regulation and function. Combining information from different "omics" levels such as transcriptomics, metabolomics, and proteomics is required to gain insights into plant metabolism and its regulation. As such, we studied fluctuations in soybean metabolism in response to R. solani infection at early and late disease stages using an integrated transcriptomics-metabolomics approach, focusing on the regulation of soybean primary metabolism and oxidative stress tolerance.

RESULTS

Transcriptomics (RNAseq) and metabolomics (1H NMR) data were analyzed individually and by integration using bidirectional orthogonal projections to latent structures (O2PLS) to reveal possible links between the metabolome and transcriptome during early and late infection stages. O2PLS analysis detected 516 significant transcripts, double that reported in the univariate analysis, and more significant metabolites than detected in partial least squares discriminant analysis. Strong separation of treatments based on integration of the metabolomes and transcriptomes of the analyzed soybean leaves was revealed, similar trends as those seen in analyses done on individual datasets, validating the integration method being applied. Strong fluctuations of soybean primary metabolism occurred in glycolysis, the TCA cycle, photosynthesis and photosynthates in response to R. solani infection. Data were validated using quantitative real-time PCR on a set of specific markers as well as randomly selected genes. Significant increases in transcript and metabolite levels involved in redox reactions and ROS signaling, such as peroxidases, thiamine, tocopherol, proline, L-alanine and GABA were also recorded. Levels of ethanol increased 24 h post-infection in soybean leaves, and alcohol dehydrogenase (ADH) loss-of-function mutants of Arabidopsis thaliana had higher necrosis than wild type plants.

CONCLUSIONS

As a proof-of-concept, this study offers novel insights into the biological correlations and identification of candidate genes and metabolites that can be used in soybean breeding for resistance to R. solani AG1-IA infection. Additionally, these findings imply that alcohol and its associated gene product ADH may have important roles in plant resistance to R. solani AG1-IA causing foliar blight.

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.

Early molecular signatures of responses of wheat to Zymoseptoria tritici in compatible and incompatible interactions

Zymoseptoria tritici, the causal agent of septoria tritici blotch, a serious foliar disease of wheat, is a necrotrophic pathogen that undergoes a long latent period. Emergence of insensitivity to fungicides, and pesticide reduction policies, mean there is a pressing need to understand septoria and control it through greater varietal resistance. Stb6 and Stb15, the most common qualitative resistance genes in modern wheat cultivars, determine specific resistance to avirulent fungal genotypes following a gene-for-gene relationship. This study investigated compatible and incompatible interactions of wheat with Z. tritici using eight combinations of cultivars and isolates, with the aim of identifying molecular responses that could be used as markers for disease resistance during the early, symptomless phase of colonization. The accumulation of TaMPK3 was estimated using western blotting, and the expression of genes implicated in gene-for-gene interactions of plants with a wide range of other pathogens was measured by qRT-PCR during the presymptomatic stages of infection. Production of TaMPK3 and expression of most of the genes responded to inoculation with Z. tritici but varied considerably between experimental replicates. However, there was no significant difference between compatible and incompatible interactions in any of the responses tested. These results demonstrate that the molecular biology of the gene-for-gene interaction between wheat and Zymoseptoria is unlike that in many other plant diseases, indicate that environmental conditions may strongly influence early responses of wheat to infection by Z. tritici, and emphasize the importance of including both compatible and incompatible interactions when investigating the biology of this complex pathosystem.

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.

The complete genome sequence of the phytopathogenic fungus Sclerotinia sclerotiorum reveals insights into the genome architecture of broad host range pathogens

Sclerotinia sclerotiorum is a phytopathogenic fungus with over 400 hosts including numerous economically important cultivated species. This contrasts many economically destructive pathogens that only exhibit a single or very few hosts. Many plant pathogens exhibit a “two-speed” genome. So described because their genomes contain alternating gene rich, repeat sparse and gene poor, repeat-rich regions. In fungi, the repeat-rich regions may be subjected to a process termed repeat-induced point mutation (RIP). Both repeat activity and RIP are thought to play a significant role in evolution of secreted virulence proteins, termed effectors. We present a complete genome sequence of S. sclerotiorum generated using Single Molecule Real-Time Sequencing technology with highly accurate annotations produced using an extensive RNA sequencing data set. We identified 70 effector candidates and have highlighted their in planta expression profiles. Furthermore, we characterized the genome architecture of S. sclerotiorum in comparison to plant pathogens that exhibit “two-speed” genomes. We show that there is a significant association between positions of secreted proteins and regions with a high RIP index in S. sclerotiorum but we did not detect a correlation between secreted protein proportion and GC content. Neither did we detect a negative correlation between CDS content and secreted protein proportion across the S. sclerotiorum genome. We conclude that S. sclerotiorum exhibits subtle signatures of enhanced mutation of secreted proteins in specific genomic compartments as a result of transposition and RIP activity. However, these signatures are not observable at the whole-genome scale.

Alterations in rice chloroplast integrity, photosynthesis and metabolome associated with pathogenesis of Rhizoctonia solani

Sheath blight disease is caused by a necrotrophic fungal pathogen Rhizoctonia solani and it continues to be a challenge for sustainable rice cultivation. In this study, we adopted a multi-pronged approach to understand the intricacies of rice undergoing susceptible interactions with R. solani. Extensive anatomical alteration, chloroplast localized ROS, deformed chloroplast ultrastructure along with decreased photosynthetic efficiency were observed in infected tissue. GC-MS based metabolite profiling revealed accumulation of glycolysis and TCA cycle intermediates, suggesting enhanced respiration. Several aromatic and aliphatic amino acids along with phenylpropanoid intermediates were also accumulated, suggesting induction of secondary metabolism during pathogenesis. Furthermore, alterations in carbon metabolism along with perturbation of hormonal signalling were highlighted in this study. The gene expression analysis including RNAseq profiling reinforced observed metabolic alterations in the infected tissues. In conclusion, the present study unravels key events associated during susceptible rice-R. solani interactions and identifies metabolites and transcripts that are accumulated in infected tissues.

Systems genetics reveals a transcriptional network associated with susceptibility in the maize-grey leaf spot pathosystem

We used a systems genetics approach to elucidate the molecular mechanisms of the responses of maize to grey leaf spot (GLS) disease caused by Cercospora zeina, a threat to maize production globally. Expression analysis of earleaf samples in a subtropical maize recombinant inbred line population (CML444 × SC Malawi) subjected in the field to C. zeina infection allowed detection of 20 206 expression quantitative trait loci (eQTLs). Four trans-eQTL hotspots coincided with GLS disease QTLs mapped in the same field experiment. Co-expression network analysis identified three expression modules correlated with GLS disease scores. The module (GY-s) most highly correlated with susceptibility (r = 0.71; 179 genes) was enriched for the glyoxylate pathway, lipid metabolism, diterpenoid biosynthesis and responses to pathogen molecules such as chitin. The GY-s module was enriched for genes with trans-eQTLs in hotspots on chromosomes 9 and 10, which also coincided with phenotypic QTLs for susceptibility to GLS. This transcriptional network has significant overlap with the GLS susceptibility response of maize line B73, and may reflect pathogen manipulation for nutrient acquisition and/or unsuccessful defence responses, such as kauralexin production by the diterpenoid biosynthesis pathway. The co-expression module that correlated best with resistance (TQ-r; 1498 genes) was enriched for genes with trans-eQTLs in hotspots coinciding with GLS resistance QTLs on chromosome 9. Jasmonate responses were implicated in resistance to GLS through co-expression of COI1 and enrichment of genes with the Gene Ontology term 'cullin-RING ubiquitin ligase complex' in the TQ-r module. Consistent with this, JAZ repressor expression was highly correlated with the severity of GLS disease in the GY-s susceptibility network.

Using Population and Comparative Genomics to Understand the Genetic Basis of Effector-Driven Fungal Pathogen Evolution

Epidemics caused by fungal plant pathogens pose a major threat to agro-ecosystems and impact global food security. High-throughput sequencing enabled major advances in understanding how pathogens cause disease on crops. Hundreds of fungal genomes are now available and analyzing these genomes highlighted the key role of effector genes in disease. Effectors are small secreted proteins that enhance infection by manipulating host metabolism. Fungal genomes carry 100s of putative effector genes, but the lack of homology among effector genes, even for closely related species, challenges evolutionary and functional analyses. Furthermore, effector genes are often found in rapidly evolving chromosome compartments which are difficult to assemble. We review how population and comparative genomics toolsets can be combined to address these challenges. We highlight studies that associated genome-scale polymorphisms with pathogen lifestyles and adaptation to different environments. We show how genome-wide association studies can be used to identify effectors and other pathogenicity-related genes underlying rapid adaptation. We also discuss how the compartmentalization of fungal genomes into core and accessory regions shapes the evolution of effector genes. We argue that an understanding of genome evolution provides important insight into the trajectory of host-pathogen co-evolution.