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

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

Dissecting the molecular interactions between wheat and the fungal pathogen Zymoseptoria tritici

The Dothideomycete fungus Zymoseptoria tritici (previously known as Mycosphaerella graminicola and Septoria tritici) is the causative agent of Septoria tritici leaf blotch (STB) disease of wheat (Triticum aestivum L.). In Europe, STB is the most economically damaging disease of wheat, with an estimated ∼€1 billion per year in fungicide expenditure directed toward its control. Here, an overview of our current understanding of the molecular events that occur during Z. tritici infection of wheat leaves is presented. On the host side, this includes the contribution of (1) the pathogen-associated molecular pattern-triggered immunity (PTI) layer of plant defense, and (2) major Stb loci for resistance against Z. tritici. On the pathogen side of the interaction, we consolidate evidence from recent bioinformatic, transcriptomic and proteomic studies that begin to explain the contribution of Z. tritici effector proteins to the biphasic lifestyle of the fungus. This includes the discovery of chitin-binding proteins in the Z. tritici secretome, which contribute to evasion of immune surveillance by this pathogen, and the possible existence of ‘necrotrophic’ effectors from Z. tritici, which may actively stimulate host recognition in a manner similar to related necrotrophic fungal pathogens. We finish by speculating on how some of these recent fundamental discoveries might be harnessed to help improve resistance to STB in the world’s second largest food crop

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

sRNA profiling combined with gene function analysis reveals a lack of evidence for cross-kingdom RNAi in the wheat – Zymoseptoria tritici pathosystem

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonization and disease development. Here we characterized RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assessed their contribution to pathogenesis. Computational analysis of fungal sRNA and host mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict their mRNA targets in wheat. 389 in planta-induced sRNA loci were identified. sRNAs generated from some of these loci were predicted to target wheat mRNAs including those potentially involved in pathogen defense. However, molecular approaches failed to validate targeting of selected wheat mRNAs by fungal sRNAs. Mutant strains of Z. tritici carrying deletions of genes encoding key components of RNAi such as Dicer-like (DCL) and Argonaute (AGO) proteins were generated, and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results did suggest the existence of non-canonical DCL-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA targeting essential fungal genes applied in vitro or generated from an RNA virus vector in planta in a procedure known as HIGS (Host-Induced Gene Silencing) was ineffective in preventing Z. tritici growth or disease. We also demonstrated that Z. tritici is incapable of dsRNA uptake. Collectively, our data suggest that RNAi approaches for gene function analyses in this fungal species and potentially also as a control measure may not be as effective as has been demonstrated for some other plant pathogenic fungi

Analysis of small RNAs in the Zymoseptoria tritici – wheat interaction

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonisation and disease development in plants. Here we characterise RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assess their contribution to pathogenesis. Computational analysis of fungal sRNA and wheat mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict putative mRNA targets in wheat. In total, 389 in planta-induced sRNA loci were identified in Z. tritici. sRNAs generated from some of these loci were predicted to target wheat mRNAs, including some that have previously been implicated in defence against pathogens. However, biochemical approaches were unable to successfully validate targeting of selected wheat mRNAs by fungal sRNAs. Z. tritici gene deletion strains deficient for key RNAi components were generated and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results do point to the existence of non-canonical Dicer-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA applied in vitro or generated from an RNA virus vector in planta was ineffective at triggering gene silencing or reducing growth of Z. tritici. We conclude that neither in vitro nor in planta RNAi approaches are likely to be useful for gene function analyses or as a viable control measure for this pathogen

Analysis of small RNA silencing in Zymoseptoria tritici – wheat interactions

Cross-kingdom small RNA (sRNA) silencing has recently emerged as a mechanism facilitating fungal colonization and disease development. Here we characterized RNAi pathways in Zymoseptoria tritici, a major fungal pathogen of wheat, and assessed their contribution to pathogenesis. Computational analysis of fungal sRNA and host mRNA sequencing datasets was used to define the global sRNA populations in Z. tritici and predict their mRNA targets in wheat. 389 in planta-induced sRNA loci were identified. sRNAs generated from some of these loci were predicted to target wheat mRNAs including those potentially involved in pathogen defense. However, molecular approaches failed to validate targeting of selected wheat mRNAs by fungal sRNAs. Mutant strains of Z. tritici carrying deletions of genes encoding key components of RNAi such as Dicer-like (DCL) and Argounate (AGO) proteins were generated, and virulence bioassays suggested that these are dispensable for full infection of wheat. Nonetheless, our results did suggest the existence of non-canonical DCL-independent pathway(s) for sRNA biogenesis in Z. tritici. dsRNA targeting essential fungal genes applied in vitro or generated from an RNA virus vector in planta in a procedure known as HIGS (Host-Induced Gene Silencing) was ineffective in preventing Z. tritici growth or disease. We also demonstrated that Z. tritici is incapable of dsRNA uptake. Collectively, our data suggest that RNAi approaches for gene function analyses in this fungal species and potentially also as a control measure may not be as effective as has been demonstrated for some other plant pathogenic fungi

DSYB catalyses the key step of dimethylsulfoniopropionate biosynthesis in many phytoplankton

Dimethylsulfoniopropionate (DMSP) is a globally important organosulfur molecule and the major precursor for dimethyl sulfide. These compounds are important info-chemicals, key nutrients for marine microorganisms, and are involved in global sulfur cycling, atmospheric chemistry and cloud formation1,2,3. DMSP production was thought to be confined to eukaryotes, but heterotrophic bacteria can also produce DMSP through the pathway used by most phytoplankton4, and the DsyB enzyme catalysing the key step of this pathway in bacteria was recently identified5. However, eukaryotic phytoplankton probably produce most of Earth’s DMSP, yet no DMSP biosynthesis genes have been identified in any such organisms. Here we identify functional dsyB homologues, termed DSYB, in many phytoplankton and corals. DSYB is a methylthiohydroxybutryate methyltransferase enzyme localized in the chloroplasts and mitochondria of the haptophyte Prymnesium parvum, and stable isotope tracking experiments support these organelles as sites of DMSP synthesis. DSYB transcription levels increased with DMSP concentrations in different phytoplankton and were indicative of intracellular DMSP. Identification of the eukaryotic DSYB sequences, along with bacterial dsyB, provides the first molecular tools to predict the relative contributions of eukaryotes and prokaryotes to global DMSP production. Furthermore, evolutionary analysis suggests that eukaryotic DSYB originated in bacteria and was passed to eukaryotes early in their evolution

Why Are Outcomes Different for Registry Patients Enrolled Prospectively and Retrospectively? Insights from the Global Anticoagulant Registry in the FIELD-Atrial Fibrillation (GARFIELD-AF).

Background: Retrospective and prospective observational studies are designed to reflect real-world evidence on clinical practice, but can yield conflicting results. The GARFIELD-AF Registry includes both methods of enrolment and allows analysis of differences in patient characteristics and outcomes that may result. Methods and Results: Patients with atrial fibrillation (AF) and ≥1 risk factor for stroke at diagnosis of AF were recruited either retrospectively (n = 5069) or prospectively (n = 5501) from 19 countries and then followed prospectively. The retrospectively enrolled cohort comprised patients with established AF (for a least 6, and up to 24 months before enrolment), who were identified retrospectively (and baseline and partial follow-up data were collected from the emedical records) and then followed prospectively between 0-18 months (such that the total time of follow-up was 24 months; data collection Dec-2009 and Oct-2010). In the prospectively enrolled cohort, patients with newly diagnosed AF (≤6 weeks after diagnosis) were recruited between Mar-2010 and Oct-2011 and were followed for 24 months after enrolment. Differences between the cohorts were observed in clinical characteristics, including type of AF, stroke prevention strategies, and event rates. More patients in the retrospectively identified cohort received vitamin K antagonists (62.1% vs. 53.2%) and fewer received non-vitamin K oral anticoagulants (1.8% vs . 4.2%). All-cause mortality rates per 100 person-years during the prospective follow-up (starting the first study visit up to 1 year) were significantly lower in the retrospective than prospectively identified cohort (3.04 [95% CI 2.51 to 3.67] vs . 4.05 [95% CI 3.53 to 4.63]; p = 0.016). Conclusions: Interpretations of data from registries that aim to evaluate the characteristics and outcomes of patients with AF must take account of differences in registry design and the impact of recall bias and survivorship bias that is incurred with retrospective enrolment. Clinical Trial Registration: - URL: http://www.clinicaltrials.gov . Unique identifier for GARFIELD-AF (NCT01090362)

Improved risk stratification of patients with atrial fibrillation: an integrated GARFIELD-AF tool for the prediction of mortality, stroke and bleed in patients with and without anticoagulation.

OBJECTIVES: To provide an accurate, web-based tool for stratifying patients with atrial fibrillation to facilitate decisions on the potential benefits/risks of anticoagulation, based on mortality, stroke and bleeding risks. DESIGN: The new tool was developed, using stepwise regression, for all and then applied to lower risk patients. C-statistics were compared with CHA2DS2-VASc using 30-fold cross-validation to control for overfitting. External validation was undertaken in an independent dataset, Outcome Registry for Better Informed Treatment of Atrial Fibrillation (ORBIT-AF). PARTICIPANTS: Data from 39 898 patients enrolled in the prospective GARFIELD-AF registry provided the basis for deriving and validating an integrated risk tool to predict stroke risk, mortality and bleeding risk. RESULTS: The discriminatory value of the GARFIELD-AF risk model was superior to CHA2DS2-VASc for patients with or without anticoagulation. C-statistics (95% CI) for all-cause mortality, ischaemic stroke/systemic embolism and haemorrhagic stroke/major bleeding (treated patients) were: 0.77 (0.76 to 0.78), 0.69 (0.67 to 0.71) and 0.66 (0.62 to 0.69), respectively, for the GARFIELD-AF risk models, and 0.66 (0.64-0.67), 0.64 (0.61-0.66) and 0.64 (0.61-0.68), respectively, for CHA2DS2-VASc (or HAS-BLED for bleeding). In very low to low risk patients (CHA2DS2-VASc 0 or 1 (men) and 1 or 2 (women)), the CHA2DS2-VASc and HAS-BLED (for bleeding) scores offered weak discriminatory value for mortality, stroke/systemic embolism and major bleeding. C-statistics for the GARFIELD-AF risk tool were 0.69 (0.64 to 0.75), 0.65 (0.56 to 0.73) and 0.60 (0.47 to 0.73) for each end point, respectively, versus 0.50 (0.45 to 0.55), 0.59 (0.50 to 0.67) and 0.55 (0.53 to 0.56) for CHA2DS2-VASc (or HAS-BLED for bleeding). Upon validation in the ORBIT-AF population, C-statistics showed that the GARFIELD-AF risk tool was effective for predicting 1-year all-cause mortality using the full and simplified model for all-cause mortality: C-statistics 0.75 (0.73 to 0.77) and 0.75 (0.73 to 0.77), respectively, and for predicting for any stroke or systemic embolism over 1 year, C-statistics 0.68 (0.62 to 0.74). CONCLUSIONS: Performance of the GARFIELD-AF risk tool was superior to CHA2DS2-VASc in predicting stroke and mortality and superior to HAS-BLED for bleeding, overall and in lower risk patients. The GARFIELD-AF tool has the potential for incorporation in routine electronic systems, and for the first time, permits simultaneous evaluation of ischaemic stroke, mortality and bleeding risks. CLINICAL TRIAL REGISTRATION: URL: http://www.clinicaltrials.gov. Unique identifier for GARFIELD-AF (NCT01090362) and for ORBIT-AF (NCT01165710)

Rapid transcriptional plasticity of duplicated gene clusters enables a clonally reproducing aphid to colonise diverse plant species

Background: The prevailing paradigm of host-parasite evolution is that arms races lead to increasing specialisation via genetic adaptation. Insect herbivores are no exception and the majority have evolved to colonise a small number of closely related host species. Remarkably, the green peach aphid, Myzus persicae, colonises plant species across 40 families and single M. persicae clonal lineages can colonise distantly related plants. This remarkable ability makes M. persicae a highly destructive pest of many important crop species. Results: To investigate the exceptional phenotypic plasticity of M. persicae, we sequenced the M. persicae genome and assessed how one clonal lineage responds to host plant species of different families. We show that genetically identical individuals are able to colonise distantly related host species through the differential regulation of genes belonging to aphid-expanded gene families. Multigene clusters collectively upregulate in single aphids within two days upon host switch. Furthermore, we demonstrate the functional significance of this rapid transcriptional change using RNA interference (RNAi)-mediated knock-down of genes belonging to the cathepsin B gene family. Knock-down of cathepsin B genes reduced aphid fitness, but only on the host that induced upregulation of these genes. Conclusions: Previous research has focused on the role of genetic adaptation of parasites to their hosts. Here we show that the generalist aphid pest M. persicae is able to colonise diverse host plant species in the absence of genetic specialisation. This is achieved through rapid transcriptional plasticity of genes that have duplicated during aphid evolution