Small RNAs (sRNAs) are important non-coding RNA regulators, playing key roles in developmental regulation, transposon suppression, environmental response, host–pathogen interaction and other diverse biological processes. However, their roles in oomycetes are poorly understood. Here, we performed sRNA sequencing and RNA sequencing of Phytophthora parasitica at stages of vegetative growth and infection of Arabidopsis roots to examine diversity and function of sRNAs in P. parasitica, a model hemibiotrophic oomycete plant pathogen. Our results indicate that there are two distinct types of sRNA-generating loci in P. parasitica genome, giving rise to clusters of 25–26 nt and 21 nt sRNAs, respectively, with no significant strand-biases. The 25–26 nt sRNA loci lie predominantly in gene-sparse and repeat-rich regions, and overlap with over 7000 endogenous gene loci. These overlapped genes are typically P. parasitica species-specific, with no homologies to the sister species P. infestans. They include approximately 40% RXLR effector genes, 50% CRN effector genes and some elicitor genes. The transcripts of most of these genes could not be detected at both the vegetative mycelium and infection stages as revealed by RNA sequencing, indicating that the 25–26 nt sRNAs are associated with efficient silencing of these genes. The 21 nt sRNA loci typically overlap with the exon regions of highly expressed genes, suggesting that the biogenesis of the 21 nt sRNAs may be dependent on the level of gene transcription and that these sRNAs do not mediate efficient silencing of homologous genes. Analyses of the published P. infestans sRNA and mRNA sequencing data consistently show that the 25–26 nt sRNAs, but not the 21 nt sRNAs, may mediate efficient gene silencing in Phytophthora.

SUMMARY Oomycetes are diploid eukaryotic microorganisms that seriously threaten sustainable crop production. MicroRNAs (miRNAs) and corresponding natural antisense transcripts (NATs) are important regulators of multiple biological processes. However, little is known about their roles in plant immunity against oomycete pathogens. In this study, we report the identification and functional characterization of miR398b and its cis‐NAT, the core‐2/I‐branching beta‐1,6‐N‐acetylglucosaminyltransferase gene (AtC2GnT), in plant immunity. Gain‐ and loss‐of‐function assays revealed that miR398b mediates Arabidopsis thaliana susceptibility to Phytophthora parasitica by targeting Cu/Zn‐Superoxidase Dismutase1 (CSD1) and CSD2, leading to suppressed expression of CSD1 and CSD2 and decreased plant disease resistance. We further showed that AtC2GnT transcripts could inhibit the miR398b–CSDs module via inhibition of pri‐miR398b expression, leading to elevated plant resistance to P. parasitica. Furthermore, quantitative reverse transcription PCR, RNA ligase‐mediated 5′‐amplification of cDNA ends (RLM‐5′ RACE), and transient expression assays indicated that miR398b suppresses the expression of AtC2GnT. We generated AtC2GnT‐silenced A. thaliana plants by CRISPR/Cas9 or RNA interference methods, and the Nicotiana benthamiana NbC2GnT‐silenced plants by virus‐induced gene silencing. Pathogenicity assays showed that the C2GnT‐silenced plants were more susceptible, while AtC2GnT‐overexpressing plants exhibited elevated resistance to P. parasitica. AtC2GnT encodes a Golgi‐localized protein, and transient expression of AtC2GnT enhanced N. benthamiana resistance to Phytophthora pathogens. Taken together, our results revealed a positive role of AtC2GnT and a negative regulatory loop formed by miR398b and AtC2GnT in regulating plant resistance to P. parasitica.

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.