Evolution of the EKA family of powdery mildew avirulence-effector genes from the ORF 1 of a LINE retrotransposon

Role of pathogen's effectors in understanding host-pathogen interaction

Pathogens can pose challenges to plant growth and development at various stages of their life cycle. Two interconnected defense strategies prevent the growth of pathogens in plants, i.e., molecular patterns triggered immunity (PTI) and pathogenic effector-triggered immunity (ETI) that often provides resistance when PTI no longer functions as a result of pathogenic effectors. Plants may trigger an ETI defense response by directly or indirectly detecting pathogen effectors via their resistance proteins. A typical resistance protein is a nucleotide-binding receptor with leucine-rich sequences (NLRs) that undergo structural changes as they recognize their effectors and form associations with other NLRs. As a result of dimerization or oligomerization, downstream components activate "helper" NLRs, resulting in a response to ETI. It was thought that ETI is highly dependent on PTI. However, recent studies have found that ETI and PTI have symbiotic crosstalk, and both work together to create a robust system of plant defense. In this article, we have summarized the recent advances in understanding the plant's early immune response, its components, and how they cooperate in innate defense mechanisms. Moreover, we have provided the current perspective on engineering strategies for crop protection based on up-to-date knowledge.

Genetic co-option into plant-filamentous pathogen interactions

Plants are engaged in a coevolutionary arms race with their pathogens that drives rapid diversification and specialization of genes involved in resistance and virulence. However, some major innovations in plant-pathogen interactions, such as molecular decoys, trans-kingdom RNA interference, two-speed genomes, and receptor networks, evolved through the expansion of the functional landscape of genes. This is a typical outcome of genetic co-option, the evolutionary process by which available genes are recruited into new biological functions. Co-option into plant-pathogen interactions emerges generally from (i) cis-regulatory variation, (ii) horizontal gene transfer (HGT), (iii) mutations altering molecular promiscuity, and (iv) rewiring of gene networks and protein complexes. Understanding these molecular mechanisms is key for the functional and predictive biology of plant-pathogen interactions.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

Comparative proteomic analysis reveals molecular differences between incompatible and compatible interaction of Erysiphe pisi in garden pea

Comparative proteome analysis of Erysiphe pisi-infected pea genotypes; JI-2480 carrying er2 resistant gene and Arkel, the susceptible genotype by liquid chromatography- mass spectrometry (LCMS/MS QTOF) at 72 h post inoculation (hpi) revealed several differentially abundant proteins (DAPs) of both the host and the pathogen. The functional annotation of proteins through gene enrichment and KEGG pathway analyses revealed strong up-regulation of pathogenesis related protein NPR1, proteins related to defense, transportation and signal transduction, hypersensitive response, cell wall modifications, phenylpropanoid and metabolic pathways in J-72. Significant abundance of membrane-related polypeptides, kinase domains and small GTPase signal transduction-related proteins suggested their major role in plant defense. The abundance of cellular antioxidant protein, catalase and its isozyme along with calreticulin-1 and 2 in J-72 confirmed their intervention in maintaining a redox balance in powdery mildew defense. High abundance levels of Glycolysis-related proteins indicated it as a major pathway for energy source during fungal growth. The majority of pathogenicity and virulence genes were downregulated in J-72 compared to A-72, while four EKA (Effectors homologues to Avk1 and Avra10) like avirulence proteins were significantly upregulated in incompatible interaction suggesting their role in eliciting hypersensitive response in pea against E. pisi.

Characterization of Two Satellite DNA Families in the Genome of the Oomycete Plant Pathogen Phytophthora parasitica

Satellite DNA is a class of repetitive sequences that are organized in long arrays of tandemly repeated units in most eukaryotes. Long considered as selfish DNA, satellite sequences are now proposed to contribute to genome integrity. Despite their potential impact on the architecture and evolution of the genome, satellite DNAs have not been investigated in oomycetes due to the paucity of genomic data and the difficulty of assembling highly conserved satellite arrays. Yet gaining knowledge on the structure and evolution of genomes of oomycete pathogens is crucial to understanding the mechanisms underlying adaptation to their environment and to proposing efficient disease control strategies. A de novo assembly of the genome of Phytophthora parasitica, an important oomycete plant pathogen, led to the identification of several families of tandemly repeated sequences varying in size, copy number, and sequence conservation. Among them, two abundant families, designated as PpSat1 and PpSat2, displayed typical features of satellite DNA and were collectively designated as PpSat. These two satellite families differ by their length, sequence, organization, genomic environment, and evolutionary dynamics. PpSat1, but not PpSat2, presented homologs among oomycetes. This observation, as well as the characterization of transcripts of PpSat families, suggested that these satellite DNA families likely play a conserved role within this important group of pathogens.

Variable genome evolution in fungi after transposon-mediated amplification of a housekeeping gene

Background

Transposable elements (TEs) can be key drivers of evolution, but the mechanisms and scope of how they impact gene and genome function are largely unknown. Previous analyses revealed that TE-mediated gene amplifications can have variable effects on fungal genomes, from inactivation of function to production of multiple active copies. For example, a DNA methyltransferase gene in the wheat pathogen Zymoseptoria tritici (synonym Mycosphaerella graminicola) was amplified to tens of copies, all of which were inactivated by Repeat-Induced Point mutation (RIP) including the original, resulting in loss of cytosine methylation. In another wheat pathogen, Pyrenophora tritici-repentis, a histone H3 gene was amplified to tens of copies with little evidence of RIP, leading to many potentially active copies. To further test the effects of transposon-aided gene amplifications on genome evolution and architecture, the repetitive fraction of the significantly expanded genome of the banana pathogen, Pseudocercospora fijiensis, was analyzed in greater detail.

Results

These analyses identified a housekeeping gene, histone H3, which was captured and amplified to hundreds of copies by a hAT DNA transposon, all of which were inactivated by RIP, except for the original. In P. fijiensis the original H3 gene probably was not protected from RIP, but most likely was maintained intact due to strong purifying selection. Comparative analyses revealed that a similar event occurred in five additional genomes representing the fungal genera Cercospora, Pseudocercospora and Sphaerulina.

Conclusions

These results indicate that the interplay of TEs and RIP can result in different and unpredictable fates of amplified genes, with variable effects on gene and genome evolution.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0177-0) contains supplementary material, which is available to authorized users.

Small RNA discovery in the interaction between barley and the powdery mildew pathogen

Background

Plants encounter pathogenic and non-pathogenic microorganisms on a nearly constant basis. Small RNAs such as siRNAs and miRNAs/milRNAs influence pathogen virulence and host defense responses. We exploited the biotrophic interaction between the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh), and its diploid host plant, barley (Hordeum vulgare) to explore fungal and plant sRNAs expressed during Bgh infection of barley leaf epidermal cells.

Results

RNA was isolated from four fast-neutron immune-signaling mutants and their progenitor over a time course representing key stages of Bgh infection, including appressorium formation, penetration of epidermal cells, and development of haustorial feeding structures. The Cereal Introduction (CI) 16151 progenitor carries the resistance allele Mla6, while Bgh isolate 5874 harbors the AVR_a6 avirulence effector, resulting in an incompatible interaction. Parallel Analysis of RNA Ends (PARE) was used to verify sRNAs with likely transcript targets in both barley and Bgh. Bgh sRNAs are predicted to regulate effectors, metabolic genes, and translation-related genes. Barley sRNAs are predicted to influence the accumulation of transcripts that encode auxin response factors, NAC transcription factors, homeodomain transcription factors, and several splicing factors. We also identified phasing small interfering RNAs (phasiRNAs) in barley that overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs).

Conclusions

These data suggest that Bgh sRNAs regulate gene expression in metabolism, translation-related, and pathogen effectors. PARE-validated targets of predicted Bgh milRNAs include both EKA (effectors homologous to AVR_k1 and AVR_a10) and CSEP (candidate secreted effector protein) families. We also identified barley phasiRNAs and miRNAs in response to Bgh infection. These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots. In addition, we identified conserved miRNAs, novel miRNA candidates, and barley genome mapped sRNAs that have PARE validated transcript targets in barley. The miRNA target transcripts are enriched in transcription factors, signaling-related proteins, and photosynthesis-related proteins. Together these results suggest both barley and Bgh control metabolism and infection-related responses via the specific accumulation and targeting of genes via sRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5947-z) contains supplementary material, which is available to authorized users.

Multiple pairs of allelic MLA immune receptor-powdery mildew AVRA effectors argue for a direct recognition mechanism

Nucleotide-binding domain and leucine-rich repeat (NLR)-containing proteins in plants and animals mediate intracellular pathogen sensing. Plant NLRs typically detect strain-specific pathogen effectors and trigger immune responses often linked to localized host cell death. The barley Mla disease resistance locus has undergone extensive functional diversification in the host population and encodes numerous allelic NLRs each detecting a matching isolate-specific avirulence effector (AVRA) of the fungal pathogen Blumeria graminis f. sp. hordei (Bgh). We report here the isolation of Bgh AVRa7, AVRa9, AVRa10, and AVRa22, which encode small secreted proteins recognized by allelic MLA7, MLA9, MLA10, and MLA22 receptors, respectively. These effectors are sequence-unrelated, except for allelic AVRa10 and AVRa22 that are co-maintained in pathogen populations in the form of a balanced polymorphism. Contrary to numerous examples of indirect recognition of bacterial effectors by plant NLRs, co-expression experiments with matching Mla-AVRa pairs indicate direct detection of the sequence-unrelated fungal effectors by MLA receptors.

LINE-1 ORF-1p enhances the transcription factor activity of pregnenolone X receptor and promotes sorafenib resistance in hepatocellular carcinoma cells

Background

LINE-1 ORF-1p is encoded by the human pro-oncogene LINE-1. Our previous work showed that LINE-1 ORF-1p could enhance the resistance of hepatocellular carcinoma (HCC) cells to antitumor agents. However, the mechanisms involved in LINE-1 ORF-1p-mediated drug resistance remain largely unknown.

Materials and methods

The endogenous mRNA level of LINE-1 ORF-1p in clinical HCC specimens was examined using quantitative PCR (qPCR). The prognosis of HCC patients was assessed using time to progression and overall survival. The transcription factor activity of pregnenolone X receptor (PXR) was examined using luciferase gene reporter assays, qPCR, chromatin immunoprecipitation assays and cellular subfraction assays. Protein interaction between LINE-1 ORF-1p and PXR was detected by co-immunoprecipitation. The effect of LINE-1 ORF-1p on sorafenib resistance in HCC cells was studied using in vitro and in vivo models.

Results

A high level of LINE-1 ORF-1p in clinical specimens was related to poor prognosis in patients who received sorafenib treatment. LINE-1 ORF-1p increased the transcription factor activity of PXR by interacting with PXR and enhancing its cytoplasmic/nuclear translocation, and recruiting PXR to its downstream gene promoter, in turn enhancing the expression of the sorafenib resistance-related genes, CYP3A4 and mdr-1. LINE-1 ORF-1p enhanced the resistance to and clearance of sorafenib in HCC cells.

Conclusion

LINE-1 ORF-1p enhances the transcription factor activation of PXR and promotes the clearance of and resistance to sorafenib in HCC cells.

Comparative transcriptome profiling of Blumeria graminis f. sp. tritici during compatible and incompatible interactions with sister wheat lines carrying and lacking Pm40

Blumeria graminis f. sp. tritici (Bgt) is an obligate biotrophic fungus that causes wheat powdery mildew, which is a devastating disease in wheat. However, little is known about the pathogenesis of this fungus, and differences in the pathogenesis of the same pathogen at various resistance levels in hosts have not been determined. In the present study, leaf tissues of both Pm40-expressing hexaploid wheat line L658 and its Pm40-deficient sister line L958 were harvested at 0 (without inoculation), 6, 12, 24, 48 and 72 hours post-inoculation (hpi) with Bgt race 15 and then subjected to RNA sequencing (RNA-seq). In addition, we also observed changes in fungal growth morphology at the aforementioned time points. There was a high correlation between percentage of reads mapped to the Bgt reference genome and biomass of the fungus within the leaf tissue during the growth process. The percentage of mapped reads of Bgt in compatible interactions was significantly higher (at the p<0.05 level) than that of reads in incompatible interactions from 24 to 72 hpi. Further functional annotations indicated that expression levels of genes encoding H+-transporting ATPase, putative secreted effector proteins (PSEPs) and heat shock proteins (HSPs) were significantly up-regulated in compatible interactions compared with these levels in incompatible interactions, particularly at 72 hpi. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis suggested that genes involved in the endocytosis pathway were also enriched in compatible interactions. Overall, genes encoding H+-transporting ATPase, PSEPs and HSPs possibly played crucial roles in successfully establishing the pathogenesis of compatible interactions during late stages of inoculation. The study results also indicated that endocytosis is likely to play a potential role in Bgt in establishing compatible interactions.

A barley powdery mildew fungus non-autonomous retrotransposon encodes a peptide that supports penetration success on barley

Pathogens overcome plant immunity by means of secreted effectors. Host effector targets often act in pathogen defense, but might also support fungal accommodation or nutrition. The barley ROP GTPase HvRACB is involved in accommodation of fungal haustoria of the powdery mildew fungus Blumeria graminis f.sp. hordei (Bgh) in barley epidermal cells. We found that HvRACB interacts with the ROP-interactive peptide 1 (ROPIP1) that is encoded on the active non-long terminal repeat retroelement Eg-R1 of Bgh. Overexpression of ROPIP1 in barley epidermal cells and host-induced post-transcriptional gene silencing (HIGS) of ROPIP1 suggested that ROPIP1 is involved in virulence of Bgh. Bimolecular fluorescence complementation and co-localization supported that ROPIP1 can interact with activated HvRACB in planta. We show that ROPIP1 is expressed by Bgh on barley and translocated into the cytoplasm of infected barley cells. ROPIP1 is recruited to microtubules upon co-expression of MICROTUBULE ASSOCIATED ROP GTPase ACTIVATING PROTEIN (HvMAGAP1) and can destabilize cortical microtubules. The data suggest that Bgh ROPIP targets HvRACB and manipulates host cell microtubule organization for facilitated host cell entry. This points to a possible neo-functionalization of retroelement-derived transcripts for the evolution of a pathogen virulence effector.

Novel jack-in-the-box effector of the barley powdery mildew pathogen?

This article comments on: Nottensteiner M, Zechmann B, McCollum C, Hückelhoven R. 2018. A barley powdery mildew fungus non-autonomous retrotransposon encodes a peptide that supports penetration success on barley. Journal of Experimental Botany 69, 3745–3758.

Signatures of host specialization and a recent transposable element burst in the dynamic one-speed genome of the fungal barley powdery mildew pathogen

Background

Powdery mildews are biotrophic pathogenic fungi infecting a number of economically important plants. The grass powdery mildew, Blumeria graminis, has become a model organism to study host specialization of obligate biotrophic fungal pathogens. We resolved the large-scale genomic architecture of B. graminis forma specialis hordei (Bgh) to explore the potential influence of its genome organization on the co-evolutionary process with its host plant, barley (Hordeum vulgare).

Results

The near-chromosome level assemblies of the Bgh reference isolate DH14 and one of the most diversified isolates, RACE1, enabled a comparative analysis of these haploid genomes, which are highly enriched with transposable elements (TEs). We found largely retained genome synteny and gene repertoires, yet detected copy number variation (CNV) of secretion signal peptide-containing protein-coding genes (SPs) and locally disrupted synteny blocks. Genes coding for sequence-related SPs are often locally clustered, but neither the SPs nor the TEs reside preferentially in genomic regions with unique features. Extended comparative analysis with different host-specific B. graminis formae speciales revealed the existence of a core suite of SPs, but also isolate-specific SP sets as well as congruence of SP CNV and phylogenetic relationship. We further detected evidence for a recent, lineage-specific expansion of TEs in the Bgh genome.

Conclusions

The characteristics of the Bgh genome (largely retained synteny, CNV of SP genes, recently proliferated TEs and a lack of significant compartmentalization) are consistent with a “one-speed” genome that differs in its architecture and (co-)evolutionary pattern from the “two-speed” genomes reported for several other filamentous phytopathogens.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4750-6) contains supplementary material, which is available to authorized users.

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.

AvrPm2 encodes an RNase-like avirulence effector which is conserved in the two different specialized forms of wheat and rye powdery mildew fungus

There is a large diversity of genetically defined resistance genes in bread wheat against the powdery mildew pathogen Blumeria graminis (B. g.) f. sp. tritici. Many confer race-specific resistance to this pathogen, but until now only the mildew avirulence gene AvrPm3a2/f2 that is recognized by Pm3a/f was known molecularly. We performed map-based cloning and genome-wide association studies to isolate a candidate for the mildew avirulence gene AvrPm2. We then used transient expression assays in Nicotiana benthamiana to demonstrate specific and strong recognition of AvrPm2 by Pm2. The virulent AvrPm2 allele arose from a conserved 12 kb deletion, while there is no protein sequence diversity in the gene pool of avirulent B. g. tritici isolates. We found one polymorphic AvrPm2 allele in B. g. triticale and one orthologue in B. g. secalis and both are recognized by Pm2. AvrPm2 belongs to a small gene family encoding structurally conserved RNase-like effectors, including Avra13 from B. g. hordei, the cognate Avr of the barley resistance gene Mla13. These results demonstrate the conservation of functional avirulence genes in two cereal powdery mildews specialized on different hosts, thus providing a possible explanation for successful introgression of resistance genes from rye or other grass relatives to wheat.

Allelic barley MLA immune receptors recognize sequence-unrelated avirulence effectors of the powdery mildew pathogen

Disease-resistance genes encoding intracellular nucleotide-binding domain and leucine-rich repeat proteins (NLRs) are key components of the plant innate immune system and typically detect the presence of isolate-specific avirulence (AVR) effectors from pathogens. NLR genes define the fastest-evolving gene family of flowering plants and are often arranged in gene clusters containing multiple paralogs, contributing to copy number and allele-specific NLR variation within a host species. Barley mildew resistance locus a (Mla) has been subject to extensive functional diversification, resulting in allelic resistance specificities each recognizing a cognate, but largely unidentified, AVR_a gene of the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). We applied a transcriptome-wide association study among 17 Bgh isolates containing different AVR_a genes and identified AVR_a1 and AVR_a13, encoding candidate-secreted effectors recognized by Mla1 and Mla13 alleles, respectively. Transient expression of the effector genes in barley leaves or protoplasts was sufficient to trigger Mla1 or Mla13 allele-specific cell death, a hallmark of NLR receptor-mediated immunity. AVR_a1 and AVR_a13 are phylogenetically unrelated, demonstrating that certain allelic MLA receptors evolved to recognize sequence-unrelated effectors. They are ancient effectors because corresponding loci are present in wheat powdery mildew. AVR_A1 recognition by barley MLA1 is retained in transgenic Arabidopsis, indicating that AVR_A1 directly binds MLA1 or that its recognition involves an evolutionarily conserved host target of AVR_A1 Furthermore, analysis of transcriptome-wide sequence variation among the Bgh isolates provides evidence for Bgh population structure that is partially linked to geographic isolation.

Avirulence Genes in Cereal Powdery Mildews: The Gene-for-Gene Hypothesis 2.0

The gene-for-gene hypothesis states that for each gene controlling resistance in the host, there is a corresponding, specific gene controlling avirulence in the pathogen. Allelic series of the cereal mildew resistance genes Pm3 and Mla provide an excellent system for genetic and molecular analysis of resistance specificity. Despite this opportunity for molecular research, avirulence genes in mildews remain underexplored. Earlier work in barley powdery mildew (B.g. hordei) has shown that the reaction to some Mla resistance alleles is controlled by multiple genes. Similarly, several genes are involved in the specific interaction of wheat mildew (B.g. tritici) with the Pm3 allelic series. We found that two mildew genes control avirulence on Pm3f: one gene is involved in recognition by the resistance protein as demonstrated by functional studies in wheat and the heterologous host Nicotiana benthamiana. A second gene is a suppressor, and resistance is only observed in mildew genotypes combining the inactive suppressor and the recognized Avr. We propose that such suppressor/avirulence gene combinations provide the basis of specificity in mildews. Depending on the particular gene combinations in a mildew race, different genes will be genetically identified as the "avirulence" gene. Additionally, the observation of two LINE retrotransposon-encoded avirulence genes in B.g. hordei further suggests that the control of avirulence in mildew is more complex than a canonical gene-for-gene interaction. To fully understand the mildew-cereal interactions, more knowledge on avirulence determinants is needed and we propose ways how this can be achieved based on recent advances in the field.

Mildew-Omics: How Global Analyses Aid the Understanding of Life and Evolution of Powdery Mildews

The common powdery mildew plant diseases are caused by ascomycete fungi of the order Erysiphales. Their characteristic life style as obligate biotrophs renders functional analyses in these species challenging, mainly because of experimental constraints to genetic manipulation. Global large-scale ("-omics") approaches are thus particularly valuable and insightful for the characterisation of the life and evolution of powdery mildews. Here we review the knowledge obtained so far from genomic, transcriptomic and proteomic studies in these fungi. We consider current limitations and challenges regarding these surveys and provide an outlook on desired future investigations on the basis of the various -omics technologies.