Diverse Evolutionary Trajectories for Small RNA Biogenesis Genes in the Oomycete Genus Phytophthora

A novel and ubiquitous miRNA-involved regulatory module ensures precise phosphorylation of RNA polymerase II and proper transcription

Proper transcription orchestrated by RNA polymerase II (RNPII) is crucial for cellular development, which is rely on the phosphorylation state of RNPII's carboxyl-terminal domain (CTD). Sporangia, developed from mycelia, are essential for the destructive oomycetes Phytophthora, remarkable transcriptional changes are observed during the morphological transition. However, how these changes are rapidly triggered and their relationship with the versatile RNPII-CTD phosphorylation remain enigmatic. Herein, we found that Phytophthora capsici undergone an elevation of Ser5-phosphorylation in its uncanonical heptapeptide repeats of RNPII-CTD during sporangia development, which subsequently changed the chromosomal occupation of RNPII and primarily activated transcription of certain genes. A cyclin-dependent kinase, PcCDK7, was highly induced and phosphorylated RNPII-CTD during this morphological transition. Mechanistically, a novel DCL1-dependent microRNA, pcamiR1, was found to be a feedback modulator for the precise phosphorylation of RNPII-CTD by complexing with PcAGO1 and regulating the accumulation of PcCDK7. Moreover, this study revealed that the pcamiR1-CDK7-RNPII regulatory module is evolutionarily conserved and the impairment of the balance between pcamiR1 and PcCDK7 could efficiently reduce growth and virulence of P. capsici. Collectively, this study uncovers a novel and evolutionary conserved mechanism of transcription regulation which could facilitate correct development and identifies pcamiR1 as a promising target for disease control.

Mutations in PpAGO3 Lead to Enhanced Virulence of Phytophthora parasitica by Activation of 25-26 nt sRNA-Associated Effector Genes

Oomycetes represent a unique group of plant pathogens that are destructive to a wide range of crops and natural ecosystems. Phytophthora species possess active small RNA (sRNA) silencing pathways, but little is known about the biological roles of sRNAs and associated factors in pathogenicity. Here we show that an AGO gene, PpAGO3, plays a major role in the regulation of effector genes hence the pathogenicity of Phytophthora parasitica. PpAGO3 was unique among five predicted AGO genes in P. parasitica, showing strong mycelium stage-specific expression. Using the CRISPR-Cas9 technology, we generated PpAGO3^ΔRGG1-3 mutants that carried a deletion of 1, 2, or 3 copies of the N-terminal RGG motif (QRGGYD) but failed to obtain complete knockout mutants, which suggests its vital role in P. parasitica. These mutants showed increased pathogenicity on both Nicotiana benthamiana and Arabidopsis thaliana plants. Transcriptome and sRNA sequencing of PpAGO3^ΔRGG1 and PpAGO3^ΔRGG3 showed that these mutants were differentially accumulated with 25–26 nt sRNAs associated with 70 predicted cytoplasmic effector genes compared to the wild-type, of which 13 exhibited inverse correlation between gene expression and 25–26 nt sRNA accumulation. Transient overexpression of the upregulated RXLR effector genes, PPTG_01869 and PPTG_15425 identified in the mutants PpAGO3^ΔRGG1 and PpAGO3^ΔRGG3, strongly enhanced N. benthamiana susceptibility to P. parasitica. Our results suggest that PpAGO3 functions together with 25–26 nt sRNAs to confer dynamic expression regulation of effector genes in P. parasitica, thereby contributing to infection and pathogenicity of the pathogen.

An Arabidopsis downy mildew non-RxLR effector suppresses induced plant cell death to promote biotroph infection

Our understanding of obligate biotrophic pathogens is limited by lack of knowledge concerning the molecular function of virulence factors. We established Arabidopsis host-induced gene silencing (HIGS) to explore gene functions of Hyaloperonospora arabidopsidis, including CYSTEINE-RICH PROTEIN (HaCR)1, a potential secreted effector gene of this obligate biotrophic pathogen. HaCR1 HIGS resulted in H. arabidopsidis-induced local plant cell death and reduced pathogen reproduction. We functionally characterized HaCR1 by ectopic expression in Nicotiana benthamiana. HaCR1 was capable of inhibiting effector-triggered plant cell death. Consistent with this, HaCR1 expression in N. benthamiana led to stronger disease symptoms caused by the hemibiotrophic oomycete pathogen Phytophthora capsici, but reduced disease symptoms caused by the necrotrophic fungal pathogen Botrytis cinerea. Expressing HaCR1 in transgenic Arabidopsis confirmed higher susceptibility to H. arabidopsidis and to the bacterial hemibiotrophic pathogen Pseudomonas syringae. Increased H. arabidopsidis infection was in accordance with reduced PATHOGENESIS RELATED (PR)1 induction. Expression of full-length HaCR1 was required for its function, which was lost if the signal peptide was deleted, suggesting its site of action in the plant apoplast. This study provides phytopathological and molecular evidence for the importance of this widespread, but largely unexplored class of non-RxLR effectors in biotrophic oomycetes.

Identification and characterization of Dicer-like genes in leaf rust pathogen (Puccinia triticina) of wheat

Puccinia triticina (P. triticina) is one of the most devastating fungal pathogens of wheat which causes significant annual yield loss to the crop. Understanding the gene regulatory mechanism of the biotrophic pathogen is one of the important aspects of host-pathogen interaction studies. Dicer-like genes are considered as important mediators of RNAi-based gene regulation. In this study, we report the presence of three Dicer-like genes (Pt-DCL1, Pt-DCL2, Pt-DCL3) in P. triticina genome identified through computational and biological analyses. Quantitative real-time PCR studies revealed an increase in the expression of these genes in germinating spore stages. Heterologous expression combined with mass spectrometry analysis of Pt-DCL2 confirmed the presence of a canonical Dicer-like gene in P. triticina. Phylogenetic analysis of the Pt-DCLs with the Dicer-like proteins from other organisms showed a distinct cluster of rust pathogens from the order Pucciniales. The results indicated a species-specific duplication of Dicer-like genes within the wheat rust pathogens. This study, for the first time, reports the presence of Dicer-dependent RNAi pathway in P. triticina that may play a role in gene regulatory mechanism of the pathogen during its development. Our study serves as a vital source of information for further RNAi-based molecular studies for better understanding and management of the wheat leaf rust disease.

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.

Oomycete small RNAs bind to the plant RNA-induced silencing complex for virulence

The exchange of small RNAs (sRNAs) between hosts and pathogens can lead to gene silencing in the recipient organism, a mechanism termed cross-kingdom RNAi (ck-RNAi). While fungal sRNAs promoting virulence are established, the significance of ck-RNAi in distinct plant pathogens is not clear. Here, we describe that sRNAs of the pathogen Hyaloperonospora arabidopsidis, which represents the kingdom of oomycetes and is phylogenetically distant from fungi, employ the host plant’s Argonaute (AGO)/RNA-induced silencing complex for virulence. To demonstrate H. arabidopsidis sRNA (HpasRNA) functionality in ck-RNAi, we designed a novel CRISPR endoribonuclease Csy4/GUS reporter that enabled in situ visualization of HpasRNA-induced target suppression in Arabidopsis. The significant role of HpasRNAs together with AtAGO1 in virulence was revealed in plant atago1 mutants and by transgenic Arabidopsis expressing a short-tandem-target-mimic to block HpasRNAs, that both exhibited enhanced resistance. HpasRNA-targeted plant genes contributed to host immunity, as Arabidopsis gene knockout mutants displayed quantitatively enhanced susceptibility.

Expansion and Divergence of Argonaute Genes in the Oomycete Genus Phytophthora

Modulation of gene expression through RNA interference is well conserved in eukaryotes and is involved in many cellular processes. In the oomycete Phytophthora, research on the small RNA machinery and function has started to reveal potential roles in the pathogen, but much is still unknown. We examined Argonaute (AGO) homologs within oomycete genome sequences, especially among Phytophthora species, to gain a clearer understanding of the evolution of this well-conserved protein family. We identified AGO homologs across many representative oomycete and stramenopile species, and annotated representative homologs in P. sojae. Furthermore, we demonstrate variable transcript levels of all identified AGO homologs in comparison to previously identified Dicer-like (DCL) and RNA-dependent RNA polymerase (RDR) homologs. Our phylogenetic analysis further refines the relationship of the AGO homologs in oomycetes and identifies a conserved tandem duplication of AGO homologs in a subset of Phytophthora species.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

Phytophthora megakarya and P. palmivora, closely related causal agents of cacao black pod rot, underwent increases in genome sizes and gene numbers by different mechanisms

Phytophthora megakarya (Pmeg) and Phytophthora palmivora (Ppal) are closely related species causing cacao black pod rot. Although Ppal is a cosmopolitan pathogen, cacao is the only known host of economic importance for Pmeg. Pmeg is more virulent on cacao than Ppal. We sequenced and compared the Pmeg and Ppal genomes and identified virulence-related putative gene models (PGeneM) that may be responsible for their differences in host specificities and virulence. Pmeg and Ppal have estimated genome sizes of 126.88 and 151.23 Mb and PGeneM numbers of 42,036 and 44,327, respectively. The evolutionary histories of Pmeg and Ppal appear quite different. Postspeciation, Ppal underwent whole-genome duplication whereas Pmeg has undergone selective increases in PGeneM numbers, likely through accelerated transposable element-driven duplications. Many PGeneMs in both species failed to match transcripts and may represent pseudogenes or cryptic genetic reservoirs. Pmeg appears to have amplified specific gene families, some of which are virulence-related. Analysis of mycelium, zoospore, and in planta transcriptome expression profiles using neural network self-organizing map analysis generated 24 multivariate and nonlinear self-organizing map classes. Many members of the RxLR, necrosis-inducing phytophthora protein, and pectinase genes families were specifically induced in planta . Pmeg displays a diverse virulence-related gene complement similar in size to and potentially of greater diversity than Ppal but it remains likely that the specific functions of the genes determine each species’ unique characteristics as pathogens.

The tRNA-Derived Small RNAs Regulate Gene Expression through Triggering Sequence-Specific Degradation of Target Transcripts in the Oomycete Pathogen Phytophthora sojae

Along with the well-studied microRNA (miRNA) and small interfering RNA (siRNA) is a new class of transfer RNA-derived small RNA (tsRNA), which has recently been detected in multiple organisms and is implicated in gene regulation. However, while miRNAs and siRNAs are known to repress gene expression through sequence-specific RNA cleavage or translational repression, how tsRNAs regulate gene expression remains unclear. Here we report the identification and functional characterization of tsRNAs in the oomycete pathogen Phytophthora sojae. We show that multiple tRNAs are processed into abundant tsRNAs, which accumulate in a similar developmental stage-specific manner and are negatively correlated with the expression of predicted target genes. Degradome sequencing and 5' RLM RACE experiments indicate tsRNAs can trigger degradation of target transcripts. Transient expression assays using GUS sensor constructs confirmed the requirement of sequence complementarity in tsRNA-mediated RNA degradation in P. sojae. Our results show that the tsRNA are a class of functional endogenous sRNAs and suggest that tsRNA regulate gene expression through inducing sequence-specific degradation of target RNAs in oomycetes.