Distinct RNA-dependent RNA polymerases are required for RNAi triggered by double-stranded RNA versus truncated transgenes in Paramecium tetraurelia

Biochemical and nanotechnological approaches to combat phytoparasitic nematodes

The foundation of most food production systems underpinning global food security is the careful management of soil resources. Embedded in the concept of soil health is the impact of diverse soil-borne pests and pathogens, and phytoparasitic nematodes represent a particular challenge. Root-knot nematodes and cyst nematodes are severe threats to agriculture, accounting for annual yield losses of US$157 billion. The control of soil-borne phytoparasitic nematodes conventionally relies on the use of chemical nematicides, which can have adverse effects on the environment and human health due to their persistence in soil, plants, and water. Nematode-resistant plants offer a promising alternative, but genetic resistance is species-dependent, limited to a few crops, and breeding and deploying resistant cultivars often takes years. Novel approaches for the control of phytoparasitic nematodes are therefore required, those that specifically target these parasites in the ground whilst minimizing the impact on the environment, agricultural ecosystems, and human health. In addition to the development of next-generation, environmentally safer nematicides, promising biochemical strategies include the combination of RNA interference (RNAi) with nanomaterials that ensure the targeted delivery and controlled release of double-stranded RNA. Genome sequencing has identified more than 75 genes in root knot and cyst nematodes that have been targeted with RNAi so far. But despite encouraging results, the delivery of dsRNA to nematodes in the soil remains inefficient. In this review article, we describe the state-of-the-art RNAi approaches targeting phytoparasitic nematodes and consider the potential benefits of nanotechnology to improve dsRNA delivery.

Challenges and Opportunities Arising from Host-Botrytis cinerea Interactions to Outline Novel and Sustainable Control Strategies: The Key Role of RNA Interference

The necrotrophic plant pathogenic fungus Botrytis cinerea (Pers., 1794), the causative agent of gray mold disease, causes significant losses in agricultural production. Control of this fungal pathogen is quite difficult due to its wide host range and environmental persistence. Currently, the management of the disease is still mainly based on chemicals, which can have harmful effects not only on the environment and on human health but also because they favor the development of strains resistant to fungicides. The flexibility and plasticity of B. cinerea in challenging plant defense mechanisms and its ability to evolve strategies to escape chemicals require the development of new control strategies for successful disease management. In this review, some aspects of the host-pathogen interactions from which novel and sustainable control strategies could be developed (e.g., signaling pathways, molecules involved in plant immune mechanisms, hormones, post-transcriptional gene silencing) were analyzed. New biotechnological tools based on the use of RNA interference (RNAi) are emerging in the crop protection scenario as versatile, sustainable, effective, and environmentally friendly alternatives to the use of chemicals. RNAi-based fungicides are expected to be approved soon, although they will face several challenges before reaching the market.

Role of Argonaute proteins in RNAi pathway in Plutella xylostella: A review

Argonaute (Ago) proteins act as key elements in RNA interference (RNAi) pathway, orchestrating the intricate machinery of gene regulation within eukaryotic cells. Within the RNAi pathway, small RNA molecules, including microRNA (miRNA), small interfering RNA (siRNA), and PIWI-interacting RNA (piRNA), collaborate with Ago family member proteins such as Ago1, Ago2, and Ago3 to form the RNA-induced silencing complex (RISC). This RISC complex, in turn, either cleaves the target mRNA or inhibits the process of protein translation. The precise contributions of Ago proteins have been well-established in numerous animals and plants, although they still remain unclear in some insect species. This review aims to shed light on the specific roles played by Ago proteins within the RNAi mechanism in a destructive lepidopteran pest, the diamondback moth (Plutella xylostella). Furthermore, we explore the potential of double-stranded RNA (dsRNA)-mediated RNAi as a robust genetic tool in pest management strategies. Through an in-depth examination of Ago proteins and dsRNA-mediated RNAi, this review seeks to contribute to our understanding of innovative approaches for controlling this pest and potentially other insect species of agricultural significance.

Paramecium Genetics, Genomics, and Evolution

The ciliate genus Paramecium served as one of the first model systems in microbial eukaryotic genetics, contributing much to the early understanding of phenomena as diverse as genome rearrangement, cryptic speciation, cytoplasmic inheritance, and endosymbiosis, as well as more recently to the evolution of mating types, introns, and roles of small RNAs in DNA processing. Substantial progress has recently been made in the area of comparative and population genomics. Paramecium species combine some of the lowest known mutation rates with some of the largest known effective populations, along with likely very high recombination rates, thereby harboring a population-genetic environment that promotes an exceptionally efficient capacity for selection. As a consequence, the genomes are extraordinarily streamlined, with very small intergenic regions combined with small numbers of tiny introns. The subject of the bulk of Paramecium research, the ancient Paramecium aurelia species complex, is descended from two whole-genome duplication events that retain high degrees of synteny, thereby providing an exceptional platform for studying the fates of duplicate genes. Despite having a common ancestor dating to several hundred million years ago, the known descendant species are morphologically indistinguishable, raising significant questions about the common view that gene duplications lead to the origins of evolutionary novelties.

Genome-Wide Identification and Characterization of Major RNAi Genes Highlighting Their Associated Factors in Cowpea (Vigna unguiculata (L.) Walp.)

In different regions of the world, cowpea (Vigna unguiculata (L.) Walp.) is an important vegetable and an excellent source of protein. It lessens the malnutrition of the underprivileged in developing nations and has some positive effects on health, such as a reduction in the prevalence of cancer and cardiovascular disease. However, occasionally, certain biotic and abiotic stresses caused a sharp fall in cowpea yield. Major RNA interference (RNAi) genes like Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) are essential for the synthesis of their associated factors like domain, small RNAs (sRNAs), transcription factors, micro-RNAs, and cis-acting factors that shield plants from biotic and abiotic stresses. In this study, applying BLASTP search and phylogenetic tree analysis with reference to the Arabidopsis RNAi (AtRNAi) genes, we discovered 28 VuRNAi genes, including 7 VuDCL, 14 VuAGO, and 7 VuRDR genes in cowpea. We looked at the domains, motifs, gene structures, chromosomal locations, subcellular locations, gene ontology (GO) terms, and regulatory factors (transcription factors, micro-RNAs, and cis-acting elements (CAEs)) to characterize the VuRNAi genes and proteins in cowpea in response to stresses. Predicted VuDCL1, VuDCL2(a, b), VuAGO7, VuAGO10, and VuRDR6 genes might have an impact on cowpea growth, development of the vegetative and flowering stages, and antiviral defense. The VuRNAi gene regulatory features miR395 and miR396 might contribute to grain quality improvement, immunity boosting, and pathogen infection resistance under salinity and drought conditions. Predicted CAEs from the VuRNAi genes might play a role in plant growth and development, improving grain quality and production and protecting plants from biotic and abiotic stresses. Therefore, our study provides crucial information about the functional roles of VuRNAi genes and their associated components, which would aid in the development of future cowpeas that are more resilient to biotic and abiotic stress. The manuscript is available as a preprint at this link: doi:10.1101/2023.02.15.528631v1.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory element analyses in sunflower (Helianthus annuus)

RNA interference (RNAi) regulates a variety of eukaryotic gene expressions that are engaged in response to stress, growth, and the conservation of genomic stability during developmental phases. It is also intimately connected to the post-transcriptional gene silencing (PTGS) process and chromatin modification levels. The entire process of RNA interference (RNAi) pathway gene families mediates RNA silencing. The main factors of RNA silencing are the Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) gene families. To the best of our knowledge, genome-wide identification of RNAi gene families like DCL, AGO, and RDR in sunflower (Helianthus annuus) has not yet been studied despite being discovered in some species. So, the goal of this study is to find the RNAi gene families like DCL, AGO, and RDR in sunflower based on bioinformatics approaches. Therefore, we accomplished an inclusive in silico investigation for genome-wide identification of RNAi pathway gene families DCL, AGO, and RDR through bioinformatics approaches such as (sequence homogeneity, phylogenetic relationship, gene structure, chromosomal localization, PPIs, GO, sub-cellular localization). In this study, we have identified five DCL (HaDCLs), fifteen AGO (HaAGOs), and ten RDR (HaRDRs) in the sunflower genome database corresponding to the RNAi genes of model plant Arabidopsis thaliana based on genome-wide analysis and a phylogenetic method. The analysis of the gene structure that contains exon-intron numbers, conserved domain, and motif composition analyses for all HaDCL, HaAGO, and HaRDR gene families indicated almost homogeneity among the same gene family. The protein-protein interaction (PPI) network analysis illustrated that there exists interconnection among identified three gene families. The analysis of the Gene Ontology (GO) enrichment showed that the detected genes directly contribute to the RNA gene-silencing and were involved in crucial pathways. It was observed that the cis-acting regulatory components connected to the identified genes were shown to be responsive to hormone, light, stress, and other functions. That was found in HaDCL, HaAGO, and HaRDR genes associated with the development and growth of plants. Finally, we are able to provide some essential information about the components of sunflower RNA silencing through our genome-wide comparison and integrated bioinformatics analysis, which open the door for further research into the functional mechanisms of the identified genes and their regulatory elements.

Species-Specific Duplication of Surface Antigen Genes in Paramecium

Paramecium is a free-living ciliate that undergoes antigenic variation and still the functions of these variable surface antigen coats in this non-pathogenic ciliate remain elusive. Only a few surface antigen genes have been described, mainly in the two model species P. tetraurelia strain 51 and P. primaurelia strain 156. Given the lack of suitable sequence data to allow for phylogenetics and deeper sequence comparisons, we screened the genomes of six different Paramecium species for serotype genes and isolated 548 candidates. Our approach identified the subfamilies of the isogenes of individual serotypes that were mostly represented by intrachromosomal gene duplicates. These showed different duplication levels, and chromosome synteny suggested rather young duplication events after the emergence of the P. aurelia species complex, indicating a rapid evolution of surface antigen genes. We were able to identify the different subfamilies of the surface antigen genes with internal tandem repeats, which showed consensus motifs across species. The individual isogene families showed additional consensus motifs, indicating that the selection pressure holds individual amino acids constant in these repeats. This may be a hint of the receptor function of these antigens rather than a presentation of random epitopes, generating the variability of these surface molecules.

Paramecium epigenetics in development and proliferation

The term epigenetics is used for any layer of genetic information aside from the DNA base-sequence information. Mammalian epigenetic research increased our understanding of chromatin dynamics in terms of cytosine methylation and histone modification during differentiation, aging, and disease. Instead, ciliate epigenetics focused more on small RNA-mediated effects. On the one hand, these do concern the transport of RNA from parental to daughter nuclei, representing a regulated transfer of epigenetic information across generations. On the other hand, studies of Paramecium, Tetrahymena, Oxytricha, and Stylonychia revealed an almost unique function of transgenerational RNA. Rather than solely controlling chromatin dynamics, they control sexual progeny's DNA content quantitatively and qualitatively. Thus epigenetics seems to control genetics, at least genetics of the vegetative macronucleus. This combination offers ciliates, in particular, an epigenetically controlled genetic variability. This review summarizes the epigenetic mechanisms that contribute to macronuclear heterogeneity and relates these to nuclear dimorphism. This system's adaptive and evolutionary possibilities raise the critical question of whether such a system is limited to unicellular organisms or binuclear cells. We discuss here the relevance of ciliate genetics and epigenetics to multicellular organisms.

Emergent RNA-RNA interactions can promote stability in a facultative phototrophic endosymbiosis

Stable endosymbiosis between eukaryotic microbes has driven the evolution of further cellular complexity. Yet the mechanisms that can act to stabilize an emergent eukaryote–eukaryote endosymbiosis are unclear. Using the model facultative endosymbiotic system, Paramecium bursaria, we demonstrate that endosymbiont–host RNA–RNA interactions can drive a cost to host growth upon endosymbiont digestion. These RNA–RNA interactions are facilitated by the host RNA-interference system. For endosymbiont messenger RNA sharing a high level of sequence identity with host transcripts, this process can result in host gene knockdown. We propose that these endosymbiont–host RNA–RNA interactions—“RNA-interference collisions”—represent an emergent mechanism to sanction the host for breakdown of the endosymbiosis, promoting the stability of the facultative endosymbiotic interaction.

Eukaryote–eukaryote endosymbiosis was responsible for the spread of chloroplast (plastid) organelles. Stability is required for the metabolic and genetic integration that drives the establishment of new organelles, yet the mechanisms that act to stabilize emergent endosymbioses—between two fundamentally selfish biological organisms—are unclear. Theory suggests that enforcement mechanisms, which punish misbehavior, may act to stabilize such interactions by resolving conflict. However, how such mechanisms can emerge in a facultative endosymbiosis has yet to be explored. Here, we propose that endosymbiont–host RNA–RNA interactions, arising from digestion of the endosymbiont population, can result in a cost to host growth for breakdown of the endosymbiosis. Using the model facultative endosymbiosis between Paramecium bursaria and Chlorella spp., we demonstrate that this mechanism is dependent on the host RNA-interference (RNAi) system. We reveal through small RNA (sRNA) sequencing that endosymbiont-derived messenger RNA (mRNA) released upon endosymbiont digestion can be processed by the host RNAi system into 23-nt sRNA. We predict multiple regions of shared sequence identity between endosymbiont and host mRNA, and demonstrate through delivery of synthetic endosymbiont sRNA that exposure to these regions can knock down expression of complementary host genes, resulting in a cost to host growth. This process of host gene knockdown in response to endosymbiont-derived RNA processing by host RNAi factors, which we term “RNAi collisions,” represents a mechanism that can promote stability in a facultative eukaryote–eukaryote endosymbiosis. Specifically, by imposing a cost for breakdown of the endosymbiosis, endosymbiont–host RNA–RNA interactions may drive maintenance of the symbiosis across fluctuating ecological conditions.

Genome-wide identification of DCL, AGO and RDR gene families and their associated functional regulatory elements analyses in banana (Musa acuminata)

RNA silencing is mediated through RNA interference (RNAi) pathway gene families, i.e., Dicer-Like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RDR) and their cis-acting regulatory elements. The RNAi pathway is also directly connected with the post-transcriptional gene silencing (PTGS) mechanism, and the pathway controls eukaryotic gene regulation during growth, development, and stress response. Nevertheless, genome-wide identification of RNAi pathway gene families such as DCL, AGO, and RDR and their regulatory network analyses related to transcription factors have not been studied in many fruit crop species, including banana (Musa acuminata). In this study, we studied in silico genome-wide identification and characterization of DCL, AGO, and RDR genes in bananas thoroughly via integrated bioinformatics approaches. A genome-wide analysis identified 3 MaDCL, 13 MaAGO, and 5 MaRDR candidate genes based on multiple sequence alignment and phylogenetic tree related to the RNAi pathway in banana genomes. These genes correspond to the Arabidopsis thaliana RNAi silencing genes. The analysis of the conserved domain, motif, and gene structure (exon-intron numbers) for MaDCL, MaAGO, and MaRDR genes showed higher homogeneity within the same gene family. The Gene Ontology (GO) enrichment analysis exhibited that the identified RNAi genes could be involved in RNA silencing and associated metabolic pathways. A number of important transcription factors (TFs), e.g., ERF, Dof, C2H2, TCP, GATA and MIKC_MADS families, were identified by network and sub-network analyses between TFs and candidate RNAi gene families. Furthermore, the cis-acting regulatory elements related to light-responsive (LR), stress-responsive (SR), hormone-responsive (HR), and other activities (OT) functions were identified in candidate MaDCL, MaAGO, and MaRDR genes. These genome-wide analyses of these RNAi gene families provide valuable information related to RNA silencing, which would shed light on further characterization of RNAi genes, their regulatory elements, and functional roles, which might be helpful for banana improvement in the breeding program.

Characterization of the RNA-interference pathway as a tool for reverse genetic analysis in the nascent phototrophic endosymbiosis, Paramecium bursaria

Endosymbiosis was fundamental for the evolution of eukaryotic complexity. Endosymbiotic interactions can be dissected through forward- and reverse-genetic experiments, such as RNA-interference (RNAi). However, distinguishing small (s)RNA pathways in a eukaryote-eukaryote endosymbiotic interaction is challenging. Here, we investigate the repertoire of RNAi pathway protein-encoding genes in the model nascent endosymbiotic system, Paramecium bursaria-Chlorella spp. Using comparative genomics and transcriptomics supported by phylogenetics, we identify essential proteome components of the small interfering (si)RNA, scan (scn)RNA and internal eliminated sequence (ies)RNA pathways. Our analyses reveal that copies of these components have been retained throughout successive whole genome duplication (WGD) events in the Paramecium clade. We validate feeding-induced siRNA-based RNAi in P. bursaria via knock-down of the splicing factor, u2af1, which we show to be crucial to host growth. Finally, using simultaneous knock-down 'paradox' controls to rescue the effect of u2af1 knock-down, we demonstrate that feeding-induced RNAi in P. bursaria is dependent upon a core pathway of host-encoded Dcr1, Piwi and Pds1 components. Our experiments confirm the presence of a functional, host-derived RNAi pathway in P. bursaria that generates 23-nt siRNA, validating the use of the P. bursaria-Chlorella spp. system to investigate the genetic basis of a nascent endosymbiosis.

In silico identification and characterization of AGO, DCL and RDR gene families and their associated regulatory elements in sweet orange (Citrus sinensis L.)

RNA interference (RNAi) plays key roles in post-transcriptional and chromatin modification levels as well as regulates various eukaryotic gene expressions which are involved in stress responses, development and maintenance of genome integrity during developmental stages. The whole mechanism of RNAi pathway is directly involved with the gene-silencing process by the interaction of Dicer-Like (DCL), Argonaute (AGO) and RNA-dependent RNA polymerase (RDR) gene families and their regulatory elements. However, these RNAi gene families and their sub-cellular locations, functional pathways and regulatory components were not extensively investigated in the case of economically and nutritionally important fruit plant sweet orange (Citrus sinensis L.). Therefore, in silico characterization, gene diversity and regulatory factor analysis of RNA silencing genes in C. sinensis were conducted by using the integrated bioinformatics approaches. Genome-wide comparison analysis based on phylogenetic tree approach detected 4 CsDCL, 8 CsAGO and 4 CsRDR as RNAi candidate genes in C. sinensis corresponding to the RNAi genes of model plant Arabidopsis thaliana. The domain and motif composition and gene structure analyses for all three gene families exhibited almost homogeneity within the same group members. The Gene Ontology enrichment analysis clearly indicated that the predicted genes have direct involvement into the gene-silencing and other important pathways. The key regulatory transcription factors (TFs) MYB, Dof, ERF, NAC, MIKC_MADS, WRKY and bZIP were identified by their interaction network analysis with the predicted genes. The cis-acting regulatory elements associated with the predicted genes were detected as responsive to light, stress and hormone functions. Furthermore, the expressed sequence tag (EST) analysis showed that these RNAi candidate genes were highly expressed in fruit and leaves indicating their organ specific functions. Our genome-wide comparison and integrated bioinformatics analyses provided some necessary information about sweet orange RNA silencing components that would pave a ground for further investigation of functional mechanism of the predicted genes and their regulatory factors.

dsRNA Uptake in Plant Pests and Pathogens: Insights into RNAi-Based Insect and Fungal Control Technology

Efforts to develop more environmentally friendly alternatives to traditional broad-spectrum pesticides in agriculture have recently turned to RNA interference (RNAi) technology. With the built-in, sequence-specific knockdown of gene targets following delivery of double-stranded RNA (dsRNA), RNAi offers the promise of controlling pests and pathogens without adversely affecting non-target species. Significant advances in the efficacy of this technology have been observed in a wide range of species, including many insect pests and fungal pathogens. Two different dsRNA application methods are being developed. First, host induced gene silencing (HIGS) harnesses dsRNA production through the thoughtful and precise engineering of transgenic plants and second, spray induced gene silencing (SIGS) that uses surface applications of a topically applied dsRNA molecule. Regardless of the dsRNA delivery method, one aspect that is critical to the success of RNAi is the ability of the target organism to internalize the dsRNA and take advantage of the host RNAi cellular machinery. The efficiency of dsRNA uptake mechanisms varies across species, and in some uptake is negligible, rendering them effectively resistant to this new generation of control technologies. If RNAi-based methods of control are to be used widely, it is critically important to understand the mechanisms underpinning dsRNA uptake. Understanding dsRNA uptake mechanisms will also provide insight into the design and formulation of dsRNAs for improved delivery and provide clues into the development of potential host resistance to these technologies.

A universal method for the rapid isolation of all known classes of functional silencing small RNAs

Diverse classes of silencing small (s)RNAs operate via ARGONAUTE-family proteins within RNA-induced-silencing-complexes (RISCs). Here, we have streamlined various embodiments of a Q-sepharose-based RISC-purification method that relies on conserved biochemical properties of all ARGONAUTEs. We show, in multiple benchmarking assays, that the resulting 15-min benchtop extraction procedure allows simultaneous purification of all known classes of RISC-associated sRNAs without prior knowledge of the samples-intrinsic ARGONAUTE repertoires. Optimized under a user-friendly format, the method – coined ‘TraPR’ for Trans-kingdom, rapid, affordable Purification of RISCs – operates irrespectively of the organism, tissue, cell type or bio-fluid of interest, and scales to minute amounts of input material. The method is highly suited for direct profiling of silencing sRNAs, with TraPR-generated sequencing libraries outperforming those obtained via gold-standard procedures that require immunoprecipitations and/or lengthy polyacrylamide gel-selection. TraPR considerably improves the quality and consistency of silencing sRNA sample preparation including from notoriously difficult-to-handle tissues/bio-fluids such as starchy storage roots or mammalian plasma, and regardless of RNA contaminants or RNA degradation status of samples.

Feeding exogenous dsRNA interferes with endogenous sRNA accumulation in Paramecium

Supply of exogenous dsRNA (exo-dsRNA), either by injection or by feeding, is a fast and powerful alternative to classical knockout studies. The biotechnical potential of feeding techniques is evident from the numerous studies focusing on oral administration of dsRNA to control pests and viral infection in crops/animal farming. We aimed to dissect the direct and indirect effects of exo-dsRNA feeding on the endogenous short interfering RNA (endo-siRNA) populations of the free-living ciliate Paramecium. We introduced dsRNA fragments against Dicer1 (DCR1), involved in RNA interference (RNAi) against exo- and few endo-siRNAs, and an RNAi unrelated gene, ND169. Any feeding, even the control dsRNA, diminishes genome wide the accumulation of endo-siRNAs and mRNAs. This cannot be explained by direct off-target effects and suggests mechanistic overlaps of the exo- and endo-RNAi mechanisms. Nevertheless, we observe a stronger down-regulation of mRNAs in DCR1 feeding compared with ND169 knockdown. This is likely due to the direct involvement of DCR1 in endo-siRNA accumulation. We further observed a cis-regulatory effect on mRNAs that overlap with phased endo-siRNAs. This interference of exo-dsRNA with endo-siRNAs warrants further investigations into secondary effects in target species/consumers, risk assessment of dsRNA feeding applications, and environmental pollution with dsRNA.

Exogenous RNAi mechanisms contribute to transcriptome adaptation by phased siRNA clusters in Paramecium

Extensive research has characterized distinct exogenous RNAi pathways interfering in gene expression during vegetative growth of the unicellular model ciliate Paramecium. However, role of RNAi in endogenous transcriptome regulation, and environmental adaptation is unknown. Here, we describe the first genome-wide profiling of endogenous sRNAs in context of different transcriptomic states (serotypes). We developed a pipeline to identify, and characterize 2602 siRNA producing clusters (SRCs). Our data show no evidence that SRCs produce miRNAs, and in contrast to other species, no preference for strand specificity of siRNAs. Interestingly, most SRCs overlap coding genes and a separate group show siRNA phasing along the entire open reading frame, suggesting that the mRNA transcript serves as a source for siRNAs. Integrative analysis of siRNA abundance and gene expression levels revealed surprisingly that mRNA and siRNA show negative as well as positive associations. Two RNA-dependent RNA Polymerase mutants, RDR1 and RDR2, show a drastic loss of siRNAs especially in phased SRCs accompanied with increased mRNA levels. Importantly, most SRCs depend on both RDRs, reminiscent to primary siRNAs in the RNAi against exogenous RNA, indicating mechanistic overlaps between exogenous and endogenous RNAi contributing to flexible transcriptome adaptation.

Functional lability of RNA-dependent RNA polymerases in animals

RNA interference (RNAi) requires RNA-dependent RNA polymerases (RdRPs) in many eukaryotes, and RNAi amplification constitutes the only known function for eukaryotic RdRPs. Yet in animals, classical model organisms can elicit RNAi without possessing RdRPs, and only nematode RNAi was shown to require RdRPs. Here we show that RdRP genes are much more common in animals than previously thought, even in insects, where they had been assumed not to exist. RdRP genes were present in the ancestors of numerous clades, and they were subsequently lost at a high frequency. In order to probe the function of RdRPs in a deuterostome (the cephalochordate Branchiostoma lanceolatum), we performed high-throughput analyses of small RNAs from various Branchiostoma developmental stages. Our results show that Branchiostoma RdRPs do not appear to participate in RNAi: we did not detect any candidate small RNA population exhibiting classical siRNA length or sequence features. Our results show that RdRPs have been independently lost in dozens of animal clades, and even in a clade where they have been conserved (cephalochordates) their function in RNAi amplification is not preserved. Such a dramatic functional variability reveals an unexpected plasticity in RNA silencing pathways.

RNA interference (RNAi) is a conserved gene regulation system in eukaryotes. In non-animal eukaryotes, it necessitates RNA-dependent RNA polymerases (“RdRPs”). Among animals, only nematodes appear to require RdRPs for RNAi. Yet additional animal clades have RdRPs and it is assumed that they participate in RNAi. Here, we find that RdRPs are much more common in animals than previously thought, but their genes were independently lost in many lineages. Focusing on a species with RdRP genes (a cephalochordate), we found that it does not use them for RNAi. While RNAi is the only known function for eukaryotic RdRPs, our results suggest additional roles. Eukaryotic RdRPs thus have a complex evolutionary history in animals, with frequent independent losses and apparent functional diversification.

Secondary amplification of siRNA machinery limits the application of spray-induced gene silencing

Spray-induced gene silencing (SIGS) is an innovative strategy for crop protection. However, the mechanism of SIGS is not known. Here, we first demonstrate that secondary small interfering RNA (siRNA) amplification limits the application of SIGS. A myosin5 gene (Myo5) was chosen as the target of SIGS in an agronomically important pathogen-Fusarium asiaticum. Five segments corresponding to the different regions of the Myo5 gene were found to efficiently silence Myo5, resulting in cell wall defects, life cycle disruption and virulence reduction. Myo5-8 (one of the Myo5 segments) induced sequence-specific RNA interference (RNAi) activity in F. asiaticum, F. graminearum, F. tricinctum and F. oxysporum, but not in other fungi, in vitro. Remarkably, the silencing of Myo5 lasted for only 9 h unless the double-stranded RNA (dsRNA) was continuously supplied, because F. asiaticum is unable to maintain siRNA amplification. After spraying on plants, dsRNAs were more efficiently taken up via the wounded surface. The antifungal activity of dsRNAs taken up by plant cells was higher and longer lasting than that dried onto the plant surface. In contrast with dsRNAs in fungi, dsRNAs in plant cells could efficiently turn into substantial siRNAs via secondary amplification machinery. Our findings provide new implications to develop SIGS as a mainstream disease control strategy against Fusarium and other fungi.

Two Sets of Piwi Proteins Are Involved in Distinct sRNA Pathways Leading to Elimination of Germline-Specific DNA

Piwi proteins and piRNAs protect eukaryotic germlines against the spread of transposons. During development in the ciliate Paramecium, two Piwi-dependent sRNA classes are involved in the elimination of transposons and transposon-derived DNA: scan RNAs (scnRNAs), associated with Ptiwi01 and Ptiwi09, and iesRNAs, whose binding partners we now identify as Ptiwi10 and Ptiwi11. scnRNAs derive from the maternal genome and initiate DNA elimination during development, whereas iesRNAs continue DNA targeting until the removal process is complete. Here, we show that scnRNAs and iesRNAs are processed by distinct Dicer-like proteins and bind Piwi proteins in a mutually exclusive manner, suggesting separate biogenesis pathways. We also demonstrate that the PTIWI10 gene is transcribed from the developing nucleus and that its transcription depends on prior DNA excision, suggesting a mechanism of gene expression control triggered by the removal of short DNA segments interrupting the gene.

•

Identification of two Piwi proteins (Ptiwi10/11) associated with iesRNAs

•

Piwi proteins bind Dicer-produced sRNAs and remove passenger strands

•

Ptiwi10 is expressed from the new somatic macronucleus

•

DNA elimination activates gene transcription

Environmental Temperature Controls Accumulation of Transacting siRNAs Involved in Heterochromatin Formation

Genes or alleles can interact by small RNAs in a homology dependent manner meaning that short interfering (siRNAs) can act in trans at the chromatin level producing stable and heritable silencing phenotypes. Because of the puzzling data on endogenous paramutations, their impact contributing to adaptive evolution in a Lamarckian manner remains unknown. An increasing number of studies characterizes the underlying siRNA accumulation pathways using transgene experiments. Also in the ciliate Paramecium tetraurelia, we induce trans silencing on the chromatin level by injection of truncated transgenes. Here, we characterize the efficiency of this mechanism at different temperatures showing that silencing of the endogenous genes is temperature dependent. Analyzing different transgene constructs at different copy numbers, we dissected whether silencing efficiency is due to varying precursor RNAs or siRNA accumulation. Our data shows that silencing efficiency correlates with more efficient accumulation of primary siRNAs at higher temperatures rather than higher expression of precursor RNAs. Due to higher primary levels, secondary siRNAs also show temperature dependency and interestingly increase their relative proportion to primary siRNAs. Our data shows that efficient trans silencing on the chromatin level in P. tetraurelia depends on environmental parameters, thus being an important epigenetic factor limiting regulatory effects of siRNAs.

Characterization of the Small RNA Transcriptome of the Marine Coccolithophorid, Emiliania huxleyi

Small RNAs (smRNAs) control a variety of cellular processes by silencing target genes at the transcriptional or post-transcription level. While extensively studied in plants, relatively little is known about smRNAs and their targets in marine phytoplankton, such as Emiliania huxleyi (E. huxleyi). Deep sequencing was performed of smRNAs extracted at different time points as E. huxleyi cells transition from logarithmic to stationary phase growth in batch culture. Computational analyses predicted 18 E. huxleyi specific miRNAs. The 18 miRNA candidates and their precursors vary in length (18–24 nt and 71–252 nt, respectively), genome copy number (3–1,459), and the number of genes targeted (2–107). Stem-loop real time reverse transcriptase (RT) PCR was used to validate miRNA expression which varied by nearly three orders of magnitude when growth slows and cells enter stationary phase. Stem-loop RT PCR was also used to examine the expression profiles of miRNA in calcifying and non-calcifying cultures, and a small subset was found to be differentially expressed when nutrients become limiting and calcification is enhanced. In addition to miRNAs, endogenous small RNAs such as ra-siRNAs, ta-siRNAs, nat-siRNAs, and piwiRNAs were predicted along with the machinery for the biogenesis and processing of si-RNAs. This study is the first genome-wide investigation smRNAs pathways in E. huxleyi. Results provide new insights into the importance of smRNAs in regulating aspects of physiological growth and adaptation in marine phytoplankton and further challenge the notion that smRNAs evolved with multicellularity, expanding our perspective of these ancient regulatory pathways.

Two sets of RNAi components are required for heterochromatin formation in trans triggered by truncated transgenes

Across kingdoms, RNA interference (RNAi) has been shown to control gene expression at the transcriptional- or the post-transcriptional level. Here, we describe a mechanism which involves both aspects: truncated transgenes, which fail to produce intact mRNA, induce siRNA accumulation and silencing of homologous loci in trans in the ciliate Paramecium. We show that silencing is achieved by co-transcriptional silencing, associated with repressive histone marks at the endogenous gene. This is accompanied by secondary siRNA accumulation, strictly limited to the open reading frame of the remote locus. Our data shows that in this mechanism, heterochromatic marks depend on a variety of RNAi components. These include RDR3 and PTIWI14 as well as a second set of components, which are also involved in post-transcriptional silencing: RDR2, PTIWI13, DCR1 and CID2. Our data indicates differential processing of nascent un-spliced and long, spliced transcripts thus suggesting a hitherto-unrecognized functional interaction between post-transcriptional and co-transcriptional RNAi. Both sets of RNAi components are required for efficient trans-acting RNAi at the chromatin level and our data indicates similar mechanisms contributing to genome wide regulation of gene expression by epigenetic mechanisms.

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

Gene regulation by small RNA pathways is ubiquitous among eukaryotes, but little is known about small RNA pathways in the Stramenopile kingdom. Phytophthora, a genus of filamentous oomycetes, contains many devastating plant pathogens, causing multibillion-dollar damage to crops, ornamental plants, and natural environments. The genomes of several oomycetes including Phytophthora species such as the soybean pathogen P. sojae, have been sequenced, allowing evolutionary analysis of small RNA-processing enzymes. This study examined the evolutionary origins of the oomycete small RNA-related genes Dicer-like (DCL), and RNA-dependent RNA polymerase (RDR) through broad phylogenetic analyses of the key domains. Two Dicer gene homologs, DCL1 and DCL2, and one RDR homolog were cloned and analyzed from P. sojae. Gene expression analysis revealed only minor changes in transcript levels among different life stages. Oomycete DCL1 homologs clustered with animal and plant Dicer homologs in evolutionary trees, whereas oomycete DCL2 homologs clustered basally to the tree along with Drosha homologs. Phylogenetic analysis of the RDR homologs confirmed a previous study that suggested the last common eukaryote ancestor possessed three RDR homologs, which were selectively retained or lost in later lineages. Our analysis clarifies the position of some Unikont and Chromalveolate RDR lineages within the tree, including oomycete homologs. Finally, we analyzed alterations in the domain structure of oomycete Dicer and RDR homologs, specifically focusing on the proposed domain transfer of the DEAD-box helicase domain from Dicer to RDR. Implications of the oomycete domain structure are discussed, and possible roles of the two oomycete Dicer homologs are proposed.

Epigenetic regulation of serotype expression antagonizes transcriptome dynamics in Paramecium tetraurelia

Phenotypic variation of a single genotype is achieved by alterations in gene expression patterns. Regulation of such alterations depends on their time scale, where short-time adaptations differ from permanently established gene expression patterns maintained by epigenetic mechanisms. In the ciliate Paramecium, serotypes were described for an epigenetically controlled gene expression pattern of an individual multigene family. Paradoxically, individual serotypes can be triggered in Paramecium by alternating environments but are then stabilized by epigenetic mechanisms, thus raising the question to which extend their expression follows environmental stimuli. To characterize environmental adaptation in the context of epigenetically controlled serotype expression, we used RNA-seq to characterize transcriptomes of serotype pure cultures. The resulting vegetative transcriptome resource is first analysed for genes involved in the adaptive response to the altered environment. Secondly, we identified groups of genes that do not follow the adaptive response but show co-regulation with the epigenetically controlled serotype system, suggesting that their gene expression pattern becomes manifested by similar mechanisms. In our experimental set-up, serotype expression and the entire group of co-regulated genes were stable among environmental changes and only heat-shock genes altered expression of these gene groups. The data suggest that the maintenance of these gene expression patterns in a lineage represents epigenetically controlled robustness counteracting short-time adaptation processes.

The Complexity of Posttranscriptional Small RNA Regulatory Networks Revealed by In Silico Analysis of Gossypium arboreum L. Leaf, Flower and Boll Small Regulatory RNAs

MicroRNAs (miRNAs) and secondary small interfering RNAs (principally phased siRNAs or trans-acting siRNAs) are two distinct subfamilies of small RNAs (sRNAs) that are emerging as key regulators of posttranscriptional gene expression in plants. Both miRNAs and secondary-siRNAs (sec-siRNAs) are processed from longer RNA precursors by DICER-LIKE proteins (DCLs). Gossypium arboreum L., also known as tree cotton or Asian cotton, is a diploid, possibly ancestral relative of tetraploid Gossypium hirsutum L., the predominant type of commercially grown cotton worldwide known as upland cotton. To understand the biological significance of these gene regulators in G. arboreum, a bioinformatics analysis was performed on G. arboreum small RNAs produced from G. arboreum leaf, flower, and boll tissues. Consequently, 263 miRNAs derived from 353 precursors, including 155 conserved miRNAs (cs-miRNAs) and 108 novel lineage-specific miRNAs (ls-miRNAs). Along with miRNAs, 2,033 miRNA variants (isomiRNAs) were identified as well. Those isomiRNAs with variation at the 3'-miRNA end were expressed at the highest levels, compared to other types of variants. In addition, 755 pha-siRNAs derived 319 pha-siRNA gene transcripts (PGTs) were identified, and the potential pha-siRNA initiators were predicted. Also, 2,251 non-phased siRNAs were found as well, of which 1,088 appeared to be produced by so-called cis- or trans-cleavage of the PGTs observed at positions differing from pha-siRNAs. Of those sRNAs, 148 miRNAs/isomiRNAs and 274 phased/non-phased siRNAs were differentially expressed in one or more pairs of tissues examined. Target analysis revealed that target genes for both miRNAs and pha-siRNAs are involved a broad range of metabolic and enzymatic activities. We demonstrate that secondary siRNA production could result from initial cleavage of precursors by both miRNAs or isomiRNAs, and that subsequently produced phased and unphased siRNAs could result that also serve as triggers of a second round of both cis- and trans-cleavage of additional siRNAs, leading to the formation of complex sRNA regulatory networks mediating posttranscriptional gene silencing. Results from this study extended our knowledge on G. arboreum sRNAs and their biological importance, which would facilitate future studies on regulatory mechanism of tissue development in cotton and other plant species.

Bioinformatics analysis of small RNAs in pima (Gossypium barbadense L.)

Small RNAs (sRNAs) are ~20 to 24 nucleotide single-stranded RNAs that play crucial roles in regulation of gene expression. In plants, sRNAs are classified into microRNAs (miRNAs), repeat-associated siRNAs (ra-siRNAs), phased siRNAs (pha-siRNAs), cis and trans natural antisense transcript siRNAs (cis- and trans-nat siRNAs). Pima (Gossypium barbadense L.) is one of the most economically important fiber crops, producing the best and longest spinnable fiber. Although some miRNAs are profiled in Pima, little is known about siRNAs, the largest subclass of plant sRNAs. In order to profile these gene regulators in Pima, a comprehensive analysis of sRNAs was conducted by mining publicly available sRNA data, leading to identification of 678 miRNAs, 3,559,126 ra-siRNAs, 627 pha-siRNAs, 136,600 cis-nat siRNAs and 79,994 trans-nat siRNAs. The 678 miRNAs, belonging to 98 conserved and 402 lineage-specific families, were produced from 2,138 precursors, of which 297 arose from introns, exons, or intron/UTR-exon junctions of protein-coding genes. Ra-siRNAs were produced from various repeat loci, while most (97%) were yielded from retrotransposons, especially LTRs (long terminal repeats). The genes encoding auxin-signaling-related proteins, NBS-LRRs and transcription factors were major sources of pha-siRNAs, while two conserved TAS3 homologs were found as well. Most cis-NATs in Pima overlapped in enclosed and convergent orientations, while a few hybridized in divergent and coincided orientations. Most cis- and trans-nat siRNAs were produced from overlapping regions. Additionally, characteristics of length and the 5’-first nucleotide of each sRNA class were analyzed as well. Results in this study created a valuable molecular resource that would facilitate studies on mechanism of controlling gene expression.

Primary and secondary siRNA synthesis triggered by RNAs from food bacteria in the ciliate Paramecium tetraurelia

In various organisms, an efficient RNAi response can be triggered by feeding cells with bacteria producing double-stranded RNA (dsRNA) against an endogenous gene. However, the detailed mechanisms and natural functions of this pathway are not well understood in most cases. Here, we studied siRNA biogenesis from exogenous RNA and its genetic overlap with endogenous RNAi in the ciliate Paramecium tetraurelia by high-throughput sequencing. Using wild-type and mutant strains deficient for dsRNA feeding we found that high levels of primary siRNAs of both strands are processed from the ingested dsRNA trigger by the Dicer Dcr1, the RNA-dependent RNA polymerases Rdr1 and Rdr2 and other factors. We further show that this induces the synthesis of secondary siRNAs spreading along the entire endogenous mRNA, demonstrating the occurrence of both 3′-to-5′ and 5′-to-3′ transitivity for the first time in the SAR clade of eukaryotes (Stramenopiles, Alveolates, Rhizaria). Secondary siRNAs depend on Rdr2 and show a strong antisense bias; they are produced at much lower levels than primary siRNAs and hardly contribute to RNAi efficiency. We further provide evidence that the Paramecium RNAi machinery also processes single-stranded RNAs from its bacterial food, broadening the possible natural functions of exogenously induced RNAi in this organism.

Genomic characterization of variable surface antigens reveals a telomere position effect as a prerequisite for RNA interference-mediated silencing in Paramecium tetraurelia

Antigenic or phenotypic variation is a widespread phenomenon of expression of variable surface protein coats on eukaryotic microbes. To clarify the mechanism behind mutually exclusive gene expression, we characterized the genetic properties of the surface antigen multigene family in the ciliate Paramecium tetraurelia and the epigenetic factors controlling expression and silencing. Genome analysis indicated that the multigene family consists of intrachromosomal and subtelomeric genes; both classes apparently derive from different gene duplication events: whole-genome and intrachromosomal duplication. Expression analysis provides evidence for telomere position effects, because only subtelomeric genes follow mutually exclusive transcription. Microarray analysis of cultures deficient in Rdr3, an RNA-dependent RNA polymerase, in comparison to serotype-pure wild-type cultures, shows cotranscription of a subset of subtelomeric genes, indicating that the telomere position effect is due to a selective occurrence of Rdr3-mediated silencing in subtelomeric regions. We present a model of surface antigen evolution by intrachromosomal gene duplication involving the maintenance of positive selection of structurally relevant regions. Further analysis of chromosome heterogeneity shows that alternative telomere addition regions clearly affect transcription of closely related genes. Consequently, chromosome fragmentation appears to be of crucial importance for surface antigen expression and evolution. Our data suggest that RNAi-mediated control of this genetic network by trans-acting RNAs allows rapid epigenetic adaptation by phenotypic variation in combination with long-term genetic adaptation by Darwinian evolution of antigen genes.

IMPORTANCE

Alternating surface protein structures have been described for almost all eukaryotic microbes, and a broad variety of functions have been described, such as virulence factors, adhesion molecules, and molecular camouflage. Mechanisms controlling gene expression of variable surface proteins therefore represent a powerful tool for rapid phenotypic variation across kingdoms in pathogenic as well as free-living eukaryotic microbes. However, the epigenetic mechanisms controlling synchronous expression and silencing of individual genes are hardly understood. Using the ciliate Paramecium tetraurelia as a (epi)genetic model, we showed that a subtelomeric gene position effect is associated with the selective occurrence of RNAi-mediated silencing of silent surface protein genes, suggesting small interfering RNA (siRNA)-mediated epigenetic cross talks between silent and active surface antigen genes. Our integrated genomic and molecular approach discloses the correlation between gene position effects and siRNA-mediated trans-silencing, thus providing two new parameters for regulation of mutually exclusive gene expression and the genomic organization of variant gene families.

Genome-wide analysis of genetic and epigenetic control of programmed DNA deletion

During the development of the somatic genome from the Paramecium germline genome the bulk of the copies of ∼45 000 unique, internal eliminated sequences (IESs) are deleted. IES targeting is facilitated by two small RNA (sRNA) classes: scnRNAs, which relay epigenetic information from the parental nucleus to the developing nucleus, and iesRNAs, which are produced and used in the developing nucleus. Why only certain IESs require sRNAs for their removal has been enigmatic. By analyzing the silencing effects of three genes: PGM (responsible for DNA excision), DCL2/3 (scnRNA production) and DCL5 (iesRNA production), we identify key properties required for IES elimination. Based on these results, we propose that, depending on the exact combination of their lengths and end bases, some IESs are less efficiently recognized or excised and have a greater requirement for targeting by scnRNAs and iesRNAs. We suggest that the variation in IES retention following silencing of DCL2/3 is not primarily due to scnRNA density, which is comparatively uniform relative to IES retention, but rather the genetic properties of IESs. Taken together, our analyses demonstrate that in Paramecium the underlying genetic properties of developmentally deleted DNA sequences are essential in determining the sensitivity of these sequences to epigenetic control.

A forward genetic screen reveals essential and non-essential RNAi factors in Paramecium tetraurelia

In most eukaryotes, small RNA-mediated gene silencing pathways form complex interacting networks. In the ciliate Paramecium tetraurelia, at least two RNA interference (RNAi) mechanisms coexist, involving distinct but overlapping sets of protein factors and producing different types of short interfering RNAs (siRNAs). One is specifically triggered by high-copy transgenes, and the other by feeding cells with double-stranded RNA (dsRNA)-producing bacteria. In this study, we designed a forward genetic screen for mutants deficient in dsRNA-induced silencing, and a powerful method to identify the relevant mutations by whole-genome sequencing. We present a set of 47 mutant alleles for five genes, revealing two previously unknown RNAi factors: a novel Paramecium-specific protein (Pds1) and a Cid1-like nucleotidyl transferase. Analyses of allelic diversity distinguish non-essential and essential genes and suggest that the screen is saturated for non-essential, single-copy genes. We show that non-essential genes are specifically involved in dsRNA-induced RNAi while essential ones are also involved in transgene-induced RNAi. One of the latter, the RNA-dependent RNA polymerase RDR2, is further shown to be required for all known types of siRNAs, as well as for sexual reproduction. These results open the way for the dissection of the genetic complexity, interconnection, mechanisms and natural functions of RNAi pathways in P. tetraurelia.

Genome-defence small RNAs exapted for epigenetic mating-type inheritance

In the ciliate Paramecium, transposable elements and their single-copy remnants are deleted during the development of somatic macronuclei from germline micronuclei, at each sexual generation. Deletions are targeted by scnRNAs, small RNAs produced from the germ line during meiosis that first scan the maternal macronuclear genome to identify missing sequences, and then allow the zygotic macronucleus to reproduce the same deletions. Here we show that this process accounts for the maternal inheritance of mating types in Paramecium tetraurelia, a long-standing problem in epigenetics. Mating type E depends on expression of the transmembrane protein mtA, and the default type O is determined during development by scnRNA-dependent excision of the mtA promoter. In the sibling species Paramecium septaurelia, mating type O is determined by coding-sequence deletions in a different gene, mtB, which is specifically required for mtA expression. These independently evolved mechanisms suggest frequent exaptation of the scnRNA pathway to regulate cellular genes and mediate transgenerational epigenetic inheritance of essential phenotypic polymorphisms.

The 5' spreading of small RNAs in Dictyostelium discoideum depends on the RNA-dependent RNA polymerase RrpC and on the dicer-related nuclease DrnB

RNA interference (RNAi) is a gene-regulatory mechanism in eukarya that is based on the presence of double stranded RNA and that can act on both, the transcription or post-transcriptional level. In many species, RNA-dependent RNA polymerases (RdRPs) are required for RNAi. To study the function of the three RdRPs in the amoeba Dictyostelium discoideum, we have deleted the encoding genes rrpA, rrpB and rrpC in all possible combinations. In these strains, expression of either antisense or hairpin RNA constructs against the transgene lacZ resulted in a 50% reduced β-Galactosidase activity. Northern blots surprisingly revealed unchanged lacZ mRNA levels, indicative of post-transcriptional regulation. Only in rrpC knock out strains, low levels of β-gal small interfering RNAs (siRNAs) could be detected in antisense RNA expressing strains. In contrast to this, and at considerably higher levels, all hairpin RNA expressing strains featured β-gal siRNAs. Spreading of the silencing signal to mRNA sequences 5′ of the original hairpin trigger was observed in all strains, except when the rrpC gene or that of the dicer-related nuclease DrnB was deleted. Thus, our data indicate that transitivity of an RNA silencing signal exists in D. discoideum and that it requires the two enzymes RrpC and DrnB.

Dynamic chromatin remodelling of ciliate macronuclear DNA as determined by an optimized chromatin immunoprecipitation (ChIP) method for Paramecium tetraurelia

We report the detailed evaluation of crucial parameters for chromatin immunoprecipitation (ChIP) of macronuclear DNA in the unicellular eukaryote Paramecium tetraurelia. Optimized parameters include crosslinking conditions, chromatin sonication and antibody titration thus providing a detailed protocol for successful ChIP in P. tetraurelia. As this ciliate is bacterivorous and RNAi by feeding represents a powerful tool for analysis of gene function, we moreover determined the effects of ingested nucleic acids by food bacteria. Feasibility of our protocol is demonstrated by characterisation of chromatin remodelling at promoters of cytosolic HSP70 isoforms during transcriptional activation under heat shock conditions by analyzing RNA abundance, nucleosome occupancy and levels of H3 lysine 9 acetylation.

RNAi effector diversity in nematodes

While RNA interference (RNAi) has been deployed to facilitate gene function studies in diverse helminths, parasitic nematodes appear variably susceptible. To test if this is due to inter-species differences in RNAi effector complements, we performed a primary sequence similarity survey for orthologs of 77 Caenorhabditis elegans RNAi pathway proteins in 13 nematode species for which genomic or transcriptomic datasets were available, with all outputs subjected to domain-structure verification. Our dataset spanned transcriptomes of Ancylostoma caninum and Oesophagostomum dentatum, and genomes of Trichinella spiralis, Ascaris suum, Brugia malayi, Haemonchus contortus, Meloidogyne hapla, Meloidogyne incognita and Pristionchus pacificus, as well as the Caenorhabditis species C. brenneri, C. briggsae, C. japonica and C. remanei, and revealed that: (i) Most of the C. elegans proteins responsible for uptake and spread of exogenously applied double stranded (ds)RNA are absent from parasitic species, including RNAi-competent plant-nematodes; (ii) The Argonautes (AGOs) responsible for gene expression regulation in C. elegans are broadly conserved, unlike those recruited during the induction of RNAi by exogenous dsRNA; (iii) Secondary Argonautes (SAGOs) are poorly conserved, and the nuclear AGO NRDE-3 was not identified in any parasite; (iv) All five Caenorhabditis spp. possess an expanded RNAi effector repertoire relative to the parasitic nematodes, consistent with the propensity for gene loss in nematode parasites; (v) In spite of the quantitative differences in RNAi effector complements across nematode species, all displayed qualitatively similar coverage of functional protein groups. In summary, we could not identify RNAi effector deficiencies that associate with reduced susceptibility in parasitic nematodes. Indeed, similarities in the RNAi effector complements of RNAi refractory and competent nematode parasites support the broad applicability of this research genetic tool in nematodes.

Many organisms regulate gene expression through an RNA interference (RNAi) pathway, first characterized in the nematode Caenorhabditis elegans. This pathway can be triggered experimentally using double-stranded (ds)RNA to selected gene targets, thereby allowing researchers to ‘silence’ individual genes and so investigate their function. It is hoped that this technology will facilitate gene silencing in important parasitic nematodes that impose a considerable health and economic burden on mankind. Unfortunately, differences in RNAi susceptibility have been observed between species. Here we investigated the possibility that differences in the complement of effector proteins involved in the RNAi pathway are responsible for these differences in susceptibility. Our data revealed that most facets of the RNAi pathway are well represented across parasitic nematodes, although there were fewer pathway proteins in other nematodes compared to C. elegans. In contrast, the proteins responsible for uptake and spread of dsRNA are not well represented in parasitic nematodes. However, the importance of these differences is undermined by our observation that the protein complements in all the parasites were qualitatively similar, regardless of RNAi-susceptibility. Clearly, differences in the RNAi pathway of parasitic nematodes do not explain the variations in susceptibility to experimental RNAi.

Functional specialization of Piwi proteins in Paramecium tetraurelia from post-transcriptional gene silencing to genome remodelling

Proteins of the Argonaute family are small RNA carriers that guide regulatory complexes to their targets. The family comprises two major subclades. Members of the Ago subclade, which are present in most eukaryotic phyla, bind different classes of small RNAs and regulate gene expression at both transcriptional and post-transcriptional levels. Piwi subclade members appear to have been lost in plants and fungi and were mostly studied in metazoa, where they bind piRNAs and have essential roles in sexual reproduction. Their presence in ciliates, unicellular organisms harbouring both germline micronuclei and somatic macronuclei, offers an interesting perspective on the evolution of their functions. Here, we report phylogenetic and functional analyses of the 15 Piwi genes from Paramecium tetraurelia. We show that four constitutively expressed proteins are involved in siRNA pathways that mediate gene silencing throughout the life cycle. Two other proteins, specifically expressed during meiosis, are required for accumulation of scnRNAs during sexual reproduction and for programmed genome rearrangements during development of the somatic macronucleus. Our results indicate that Paramecium Piwi proteins have evolved to perform both vegetative and sexual functions through mechanisms ranging from post-transcriptional mRNA cleavage to epigenetic regulation of genome rearrangements.