Figures
Abstract
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.
Author Summary
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.
Citation: Dalzell JJ, McVeigh P, Warnock ND, Mitreva M, Bird DM, Abad P, et al. (2011) RNAi Effector Diversity in Nematodes. PLoS Negl Trop Dis 5(6): e1176. https://doi.org/10.1371/journal.pntd.0001176
Editor: Matty Knight, Biomedical Research Institute, United States of America
Received: March 3, 2011; Accepted: April 4, 2011; Published: June 7, 2011
Copyright: © 2011 Dalzell et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was funded in part by The Department of Education and Learning for Northern Ireland, The Department of Agriculture and Rural Development for Northern Ireland and The National Institutes of Health (NIH R01 AI49162, to AGM, TAD). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
RNA interference (RNAi) is a reverse genetics technique which permits the ablation of mRNA by introduction of complementary double-stranded RNA (dsRNA), through cellular mechanisms common to most eukaryotes (for review, see [1]) and provides a functional genomics platform in a range of organisms, including those intractable to traditional genetic manipulations. One such group of organisms are the parasitic nematodes for which there have been recent expansions in transcriptomic and genomic datasets [2]–[5].
Several groups have attempted to apply the RNAi protocols pioneered in Caenorhabditis elegans to parasitic nematodes. Significant progress has been made in plant-parasitic nematodes (PPNs) in which RNAi is an established experimental technique [6]–[8], and may have utility for parasite control in plants genetically engineered to express PPN-transcript-specific dsRNA [9], [10]. In contrast, RNAi experiments in animal- and human-parasitic nematodes have had variable levels of success (for reviews, see [11]–[13]). Of note are experiments reporting inefficient or inconsistent transcript knockdown, highlighted by successful silencing of only 3 of 8 Ostertagia ostertagi genes [14] and 2 of 11 Haemonchus contortus genes [15]. In H. contortus, one feature of successful RNAi appears to be the location of target gene expression, since genes predicted to be expressed in environmentally-exposed tissues are more readily silenced [16]. RNAi difficulties have also been seen in Heligmosomoides polygyrus [17] and the non-parasitic species Pristionchus pacificus and Oscheius sp1 CEW1 [18]–[20]. Notably, inter-species differences are apparent even within the genus Caenorhabditis, where C. briggsae (unlike C. elegans) is unable to take up dsRNA from the environment, due to a SID-2 which displays aberrant RNAi functionality [21].
Hypotheses to explain RNAi difficulties in parasitic nematodes have been reported, and include: (i) the lack of appropriate in vitro culture systems for parasitic nematodes [15]; (ii) inappropriate methods of dsRNA delivery, i.e. delivered externally, where microinjection directly into the worm is more effective in C. elegans [13]; (iii) differences in RNAi effector protein functionality [13], [15]; and (iv) differences in the complement of RNAi effectors between nematodes [12], [13], [15], [17]. The latter hypothesis has been confirmed for the apicomplexan Plasmodium spp. (the causative agents of malaria), which are refractory to RNAi due to deficiencies in key pathway components [22]–[24]. Here, we test this hypothesis in nematodes by investigating the complement of RNAi pathway proteins in selected nematode datasets. Using 77 C. elegans RNAi pathway proteins as query sequences, we performed BLAST trawls of nematode-derived genomic and transcriptomic resources. Our searches focused on high-quality sequence datasets, including the draft genomes of Trichinella spiralis (Clade I/clade 2; here and throughout, we utilize clade delineations of both Blaxter et al. (denoted clades I–V [25]) and Holterman et al. (denoted clades 1–12 [26]), Ascaris suum (Clade III/clade 8), Brugia malayi (Clade III/clade 8), Meloidogyne incognita (Clade IV/clade 12), Meloidogyne hapla (Clade IV/clade 12), Caenorhabditis brenneri, Caenorhabditis briggsae, Caenorhabditis japonica, Caenorhabditis remanei (Clade V/clade 9), Haemonchus contortus (Clade V/clade 9), and Pristionchus pacificus (Clade V/clade 9) as well as the transcriptomes of Oesophagostomum dentatum (Clade V/clade 9) and Ancylostoma caninum (Clade V/clade 9). We find that the RNAi effector complements of these species, whilst quantitatively different are qualitatively similar with regard to the presence of functional groupings, yielding no major inter-species differences except that all were notably less diverse than in Caenorhabditis spp. These data suggest that variable susceptibilities to RNAi amongst parasitic nematodes cannot be adequately explained by differences in RNAi effector complement between such species.
Materials and Methods
Reciprocal BLAST Methodology
Seventy-seven C. elegans proteins known to be involved in core aspects of RNAi were identified from literature (Figure 1). These proteins were separated into five core functional groups; namely, small RNA biosynthesis, dsRNA uptake and spreading, AGOs and RISC, RNAi inhibitors, and nuclear effectors. Protein sequences were retrieved from WormBase (www.wormbase.org; release WS206) and used as search strings in a series of primary translated nucleotide (tBLASTn) and protein BLASTs (BLASTp) [27] against genome and transcriptome databases described below. All primary BLAST hits returning with a bitscore ≥40 and an expect value ≤0.01 were manually translated to amino acid sequence in six reading frames (www.expasy.ch/tools/dna.html), and analysed for identity and domain structure by BLASTp (through NCBI's Conserved Domain Database service) and InterProScan (www.ebi.ac.uk/Tools/InterProScan). The appropriate reading frame in each case (usually that with the largest uninterrupted open reading frame [ORF], however this was determined empirically on a case by case basis) was then subjected to reciprocal tBLASTn and BLASTp against the C. elegans non-redundant nucleotide and protein databases on the NCBI BLAST server (http://www.ncbi.nlm.nih.gov/BLAST), using default settings. The identity of the top-scoring reciprocal BLAST hit was accepted as identity of the relevant primary hit, as long as that identity was also supported by domain structure analysis (see Datasets S1, S2, S3, S4, S5). In the case of H. contortus, primary tBLASTn searches were performed and the separate high scoring return sequences were concatenated into a single sequence (to facilitate reciprocation) and used as reciprocal tBLASTn and BLASTp searches against C. elegans, as before.
(1) Exogenously applied double-stranded RNA (dsRNA) and small interfering RNA (exo-siRNA) are thought to enter cells via SID (Systemic RNA Interference Defective) proteins SID-1/RSD-8 and SID-2. (2) Endogenous RNAi-based pathways begin in the nucleus; micro-interfering RNA (miRNA) synthesis begins with transcription of hairpin-looped primary miRNA (pri-miRNA) transcripts from intergenic, intronic or antisense regions. pri-miRNAs are processed by the DRSH-1/PASH-1 complex to pre-miRNA, which are exported from the nucleus by exportin proteins XPO-1, -2 and -3. Endogenous siRNAs (endo-siRNAs) are also produced from genomic regions, and exported by XPO-1, 2, and -3. (3) Both pre-miRNAs and exogenously applied dsRNA molecules are bound and cleaved by the dicer complex, which consists of the RNAse III-like nuclease DCR-1, the dsRNA-binding proteins RDE-1 and -4, the helicases DRH-1 and DRH-3/EKL-3, the RNA-dependent RNA-polymerase (RdRP) RRF-1, and the uncharacterized protein, AIN-1. Dicer cleaves dsRNA to produce siRNA molecules, and pre-miRNA to mature miRNA, both of which are substrates for the RNA-induced silencing complex (RISC). (4) Both siRNAs and miRNAs are the focus of a battery of inhibitors, which allow down-regulation of the RNAi response. (5) The RISC complex incorporates a single strand of miRNA or siRNA (termed the guide strand), and binds a complementary mRNA strand, eliciting gene silencing by either mRNA destruction or translational repression (6). The central catalytic component of RISC is an argonaute (AGO) protein, allied with the nuclease TSN-1, the RNA-binding protein VIG-1, and AIN-1. (7) The RNAi response may be amplified by the action of the RdRPs RRF-1 and -2, SMG-5, RDE-2/MUT-8 and MUT-7, which produce a population of single-stranded RNAs bearing N-terminal tri-phosphates from a target mRNA template. (8) These secondary siRNAs interact with Secondary-siRNA-specific AGOs (SAGO-1 and -2), terminating in down-regulation of target transcript. Secondary siRNAs can also spread between cells through RSD-2, -3 and -6, resulting in intercellular spread of the RNAi effect (9), and can be imported into the nucleus by NRDE-3, which elicits transcriptional silencing of nascent RNA transcripts as part of nuclear RISC (nucRISC) (10). siRNAs may also control aspects of nuclear RNAi, including histone methylation, chromatin formation and chromosome segregation (11). Dashed lines indicate miRNA-based pathways, solid lines indicate siRNA-based pathways.
Databases
The M. incognita (http://www.inra.fr/meloidogyne_incognita/genomic_resources) and B. malayi (http://blast.jcvi.org/er-blast/index.cgi?project=bma1) genomes were searched using BLASTp to predicted protein sets, in addition to tBLASTn against available contig assembly, unplaced reads and associated ESTs [3], [5]. The M. hapla genome was searched using BLASTp against public release 4 (HapPep4: www.hapla.org) of the hand annotated and experimentally-validated M. hapla protein set [28], in addition to tBLASTn against the 10× contig assembly [4]. The H. contortus genome was searched using tBLASTn against the supercontig 26/08/09 database (http://www.sanger.ac.uk/cgi-bin/blast/submitblast/h_contortus). A. suum, A. caninum, T. spiralis and O. dentatum primary BLASTp and tBLASTn searches were performed using the datasets generated at Washington University, St Louis (available at www.nematode.net, [29]), as above; reciprocal BLAST searches against C. elegans datasets were then performed as before. Using the core eukaryotic genes as a reference [30], we estimated that 93% of the A. caninum [31]; and 87% of the O. dentatum transcriptome is identified, making these two dataset comparable to the full proteomes predicted from the genomes of the other species included in this study. C. brenneri, C. briggsae, C. japonica and C. remanei datasets were accessed through WormBase. Searches were also performed against publically-available nematode expressed sequence tags (ESTs) available through GenBank (www.ncbi.nlm.nih.gov), using methods as described above.
Results and Discussion
In the absence of tractable methods for transgenesis or genetic manipulation, RNAi represents one of the few molecular genetics tools that can currently be applied to parasitic nematodes. However, reports documenting difficulties in the application of RNAi to some nematodes [11]–[13], [16], [17] suggest diminished potential for RNAi-based gene function and target validation studies in some species. We have employed primary sequence similarity-based methodology to identify putative orthologs of C. elegans RNAi pathway (Figure 1) proteins in a dataset of 13 nematode genomes/transcriptomes, as a means of investigating the inter-species conservation of RNAi effectors which might inform the wider utility of RNAi in parasitic nematodes. We selected these datasets in the first instance as those represented by predicted protein sets, which were most consistent with our primary protein similarity-based search methodology. Latterly, we extended our searches to include the publically-available genome construct of H. contortus (which at the time of searching lacked adequate gene predictions), due to the central importance of this species to the question of RNAi's applicability to animal-parasitic nematodes. While we recognise that we have omitted several other publically-available nematode genome datasets from our analyses, we considered that they did not meet our criteria for inclusion, as described above. Whilst several authors refer to the apparent presence/absence of a few RNAi effector proteins in single nematode species, one small scale study examined the occurrence of 18 such proteins across H. contortus, B. malayi and the flatworm parasite Schistosoma mansoni [17]. The putative orthologs that we identified are summarised in Tables 1–5, with corresponding protein sequences described in Datasets S1, S2, S3, S4, S5. While we addressed publically-available nematode ESTs in our searches, these contributed little to our analyses due to the fragmentary nature of their coverage of RNAi effector proteins (see Table S1).
C. elegans displays an expanded repertoire of RNAi effectors relative to other nematodes
Perhaps the most striking observation is that each of the parasite species considered here possessed only a fraction of our original search set of 77 C. elegans RNAi proteins (Table 6), with all displaying a greatly contracted suite of RNAi effector proteins; of the original 77 C. elegans search strings, H. contortus returned 46, A. suum 44, A. caninum 40, O. dentatum 38, P. pacificus 36, B. malayi 35, M. hapla 28, M. incognita 27, and T. spiralis 22. This reduction in diversity (which could suggest either that: (i) orthologs of the C. elegans proteins are absent from the species in question; (ii) they have diverged to a degree that is unrecognisable on a primary sequence level, or (iii) our datasets possess significant areas of inadequate coverage such that additional RNAi effector genes await discovery in these species) was observed across all of the functional groupings in our dataset, but was most pronounced within the proteins responsible for uptake/spread of dsRNA. In contrast, the other Caenorhabditid species possessed an RNAi effector complement much closer to that of C. elegans; C. briggsae 65, C. remanei 65, C. brenneri 63, and C. japonica 60 (Table 6). However, both parasitic and free-living species returned only a subset of putative AGO orthologs relative to C. elegans. AGO analysis presented a significant challenge within our sequence similarity searches, due in part to significant areas of sequence similarity between functionally disparate C. elegans proteins. In many cases our BLAST analysis presented a clustering of multiple distinct AGOs around an individual C. elegans ortholog. Additionally, in some examples we could identify putative AGO orthologs which reciprocated to non-cleavage competent C. elegans proteins, but which encoded catalytic residues consistent with cleavage-competency themselves [32]. Clearly, using gross sequence similarity as an identification tool for AGOs underestimates functional diversity (data not shown), and as a result, we considered that an in depth analysis of AGO family diversity was beyond the scope of this study. This did not represent an issue for the analysis of other RNAi pathway protein families.
Small RNA biosynthesis
Small RNA-based genetic regulatory pathways are ubiquitous in eukaryotes, and represent a set of proteins with conserved function and structure in evolutionarily distant organisms. As such, our analysis of proteins that perform nuclear biosynthesis, nuclear export and cytoplasmic processing of small RNAs such as miRNAs (Figure 1; for recent review, see [1]) should provide a positive control measure for both our approach, and sequence data quality. These core proteins were well conserved within our dataset (Table 1; Dataset S1) - transcripts encoding many of the proteins required for siRNA and miRNA processing, including RNase III enzymes (drosha, DRSH-1; pasha, PASH-1; dicer, DCR-1), RNA helicases (dicer-related helicases DRH-1 and -3), and exportins (XPO-1 and -3) are highly conserved across the genomic and transcriptomic datasets considered here, although orthologs of the dsRNA-binding protein and dicer-complex cofactor, RDE-4, were notably absent from all of the parasites except B. malayi and A. caninum.
dsRNA uptake and spreading
Our dataset recognizes five C. elegans genes putatively responsible for dsRNA uptake and spread, identified from mutant screens for defects in systemic RNAi (the RNAi spreading defective mutants rsd-2, -3 and -6, and the systemic RNAi defective mutants sid-1 and -2). Much interest has centered on SIDs as core determinants of dsRNA uptake/spreading mechanisms. These transmembrane proteins were first described in C. elegans as mediators of systemic and environmental RNAi due to their role in transmembrane transport of dsRNA [21], [33]. Putative SID orthologs have since been described in disparate organisms including mammalian cells [34], trematode flatworms [35], crustaceans [36] and insects [37], [38] (although Drosophila melanogaster does not possess known SID orthologs, heterologous expression of C. elegans SID-1 sensitizes Drosophila cells to RNAi by soaking [39]). Similarly, expression of SID-1 in C. elegans neurons reverses the neuronal intractability of this species [40]. The role of SID-2 in environmental RNAi has been demonstrated by functional expression of C. elegans SID-2 in C. briggsae, a transformation which confers susceptibility to environmental RNAi in this species [21]. Given the importance of SID-1 and -2 to functional RNAi in C. elegans, it is surprising that these proteins are so poorly conserved in other nematodes, where putative SID-1 orthologs were identified in H. contortus and O. dentatum only (Table 2) and sid-2 was not identified outside the Caenorhabditis genus. Similarly poor conservation was observed with RSD-2 (not identified) and RSD-6 (seen only in P. pacificus). RSD-3 is the sole perfectly conserved spreading protein in our dataset, occurring in all 13 species (see Table 2; Dataset S2). Evidence from C. elegans implicates RSD-3 in intercellular spread since rsd-3 null mutants are able to take up dsRNA from the gut lumen, but are unable to distribute this dsRNA into the germline [41]. Despite lacking identifiable orthologs of SID-1, and -2, as well as RSD-2 and -6, plant-parasitic Meloidogyne and Globodera spp. display systemic RNAi following soaking in dsRNA/siRNA [7], [8], [42], [43], suggesting that alternative uptake proteins (e.g. fed mutants; see [44]), or mechanisms are involved, perhaps similar to the receptor-mediated endocytotic dsRNA uptake process seen in insect gut cells [45]. Intriguingly, our own unpublished data demonstrate a phenomenon of well conserved miRNA target transcript up-regulation in response to dsRNA/siRNA soaking of M. incognita, G. pallida and A. suum, possibly in response to a ubiquitous saturation of RNAi pathway effectors shared between exogenous (dsRNA/siRNA) and endogenous (miRNA) small RNA pathways, which could indicate that uptake is not limiting for these nematodes ([46]; unpublished observations). Additionally, we cannot discount the possibility that poorly-characterised morphological differences, such as cuticle permeability, better enable dsRNA uptake or propagation in PPNs relative to other parasite species.
Secondary siRNA amplification
In C. elegans, plants [47], and Neurospora [48], the RNAi effect is greatly amplified by the action of RNA-dependent RNA polymerases (RdRPs), which produce a population of secondary siRNAs from the target mRNA template [41], [49]–[52]. Further examples of RdRP-catalyzed amplification mechanisms have recently been reported in Paramecium tetraurelia, where multiple RdRPs appear to exist [53], and in Drosophila, where a non-canonical RdRP has been identified [54]. The most well-conserved RdRP in our dataset is EGO-1 (Enhancer of Glp-One [glp-1]), which appears in seven species (Table 2). RRF-3 (RNA-dependent RNA polymerase family member 3), which coordinates complex and ill-understood interactions between RNAi inhibition and amplification of the secondary siRNA response is reasonably well conserved, with RRF-1 less so. EGO-1 is an RdRP with core functions in transcription of “WAGO” (worm-specific AGO [55])-interacting 22G-RNAs responsible for silencing events involved in genome surveillance [56], [57] and with additional roles in germline development [58], heterochromatin assembly [59], [60], holocentric chromosome segregation [61], and P-granule function [62]. In light of these core roles, the inter-species conservation of EGO-1 is unsurprising. RRF-3, which is also reasonably well-conserved, was traditionally referred to as an inhibitory RdRP [63], although through recent work has been implicated in the production of secondary 26G-RNAs which seed a two-step process of secondary amplification against endogenous targets (endo-siRNAs) [57], [64], [65]. It is also believed that nonsense-mediated decay (NMD) proteins SMG-2 (Suppressor with Morphological effects on Genitalia 2), -5 and -6 may play a role in the induction and maintenance of secondary amplification [66], a hypothesis supported by analysis of smg null mutants which are defective for RNAi initiation [67]. SMG-2 and -6 are perfectly conserved across the genomes and transcriptomes considered here, while SMG-5 is not well conserved (see Table 2; Dataset S2). Conservation of EGO-1 suggests that all of the nematode species examined here are capable of some degree of secondary RNAi amplification, consistent with previous observations of the potency of RNAi in PPNs, where soaking in as little as 0.1 µg/ml dsRNA was capable of eliciting significant and consistent knockdown of transcripts in Globodera pallida and M. incognita second stage juveniles (J2s) [8].
AGOs and RISC
C. elegans possesses at least 27 distinct AGOs (including pseudogenes C06A1.4 and C14B1.7) [32], which constitute the central effectors of the RNA-induced silencing complex (RISC), conferring both function and specificity to RISC. All of the nematodes in our dataset possessed multiple distinct AGOs (Table 3). A subset of well-conserved AGOs (defined according to closest C. elegans BLAST match) included the miRNA-interacting AGO, ALG-1 (Argonaute [Plant]-Like Gene), as well as several endo-siRNA-interacting AGOs including the 26G-RNA-interacting ALG-4 [68], and the 22G-RNA-interacting WAGOs, R06C7.1 and F58G1.1 [55]. Some members of the PIWI-clade of AGOs, such as PRG-1 (Piwi-Related Gene 1), PRG-2, ERGO-1 (Endogenous Rnai deficient arGOnaute 1) and the AGO/PIWI-clade secondary AGOs SAGO-1 and SAGO-2, are not well conserved. Surprisingly, RDE-1, which is believed to be the main AGO involved in silencing events triggered by exogenous dsRNA in C. elegans, was only identified in the animal parasitic nematodes A. suum, H. contortus and A. caninum. Thus the AGOs known in C. elegans to be responsible for endogenous regulation of gene expression are well conserved, while the AGOs responsible for executing RNAi triggered by exogenous dsRNA are not. However, as previously stated, our identification strategy does not account for the possibility that other uncharacterized AGOs exist in each nematode species, performing roles comparable to those AGOs which we could not identify. A further four C. elegans AGOs (M03D4.7; T23D8.7; ZK218.8, NRDE-3) did not appear to be present within our parasite dataset. The AGO NRDE-3, is responsible for nuclear translocation of RNAi triggers in C. elegans, and is involved in processes which lead to heritability of gene silencing events. As NRDE-3 is completely absent from the parasite datasets considered here, this may indicate that silencing events cannot be passed between generations of parasitic nematodes. Our data suggest that most nematodes have smaller AGO complements than C. elegans, although the impact this has on functional diversity is unknown. The contracted complement of AGOs identified in the parasite species relative to C. elegans is consistent with their propensity for gene loss [69]. This could indicate redundancy in the function of individual AGOs within C. elegans, or conversely a reduced functionality within the parasites considered here. Interestingly, ERGO-1 is involved in the function of endogenous siRNA populations within C. elegans [57], [65] but is poorly conserved perhaps indicating a differential small RNA population dynamic between species. Again, the poor conservation of such proteins in RNAi-competent plant-parasitic species would seem to suggest that such deficiencies need not undermine RNAi functionality.
In addition to the catalytic AGO protein, RISCs also comprise several protein co-factors, including multiple dsRNA-binding proteins and exonucleases which are thought to pass from elements of the biosynthetic machinery (Figure 1), although these co-factors are in fact quite poorly characterized, even in C. elegans. Our analysis reveals that TSN-1 (Tudor Staphylococcal Nuclease 1), which is a common component of RISC in C. elegans, Drosophila and mammalian cells [70], is well conserved across the species considered here (Table 3; Dataset S3). The ALG interacting protein AIN-1, responsible for targeting miRNA-bound ALGs to P-bodies [71], [72], is also reasonably well-conserved, being present in seven species. VIG-1, the C. elegans ortholog of Drosophila VASA intronic gene which regulates transition between larval and adult cellular fates though interaction with the let-7 miRNA [73], was identified in five of our eight species.
RNAi inhibitors
Proteins with RNAi-inhibiting function were first characterized in C. elegans, leading to the identification of RNAi-hypersensitive null mutant strains of RRF-3 [63] and ERI-1 [74]. Only two RNAi inhibitor orthologs, the DEDDh-like 3′-5′ siRNA exonuclease ERI-1 and the miRNA 5′-3′ exonuclease XRN-2 (XRN RiboNuclease related 2), are fully conserved across our genomic and transcriptomic datasets (Table 4; Dataset S4). Sporadically-conserved inhibitors included the adenosine deaminases ADR-1 and -2 [75], and LIN-15b, while orthologs of ERI-3, -5 and -6/7 [76] were not identified outside Caenorhabditis spp.
Nuclear effectors
The RNAi pathway affects a number of poorly understood nuclear silencing mechanisms. We found that an uncharacterized nuclear effector, EKL-1 (Enhancer of KSR-1 Lethality 1 [KSR-1 is a Ras-ERK signaling scaffold protein] [77]) was the most highly conserved between species (Table 5; Dataset S5). Other chromatin-associated proteins, helicases and methylation factors are conserved to varying degrees, however MES-3 (Maternal Effect Sterile 3), RDE-2 (RNAi Defective 2), EKL-5 and MUT-16 were only found in Caenorhabditis spp.
Conclusions
In spite of the contrasting experimental evidence from published studies, our data indicate that diverse nematode species possess the machinery required to facilitate an RNAi response. Our inability to culture many animal parasitic nematodes under in vitro conditions may represent one of the main reasons why RNAi is difficult to perform in these species. Certainly, where RNAi has been most successful in nematodes it has been in species/life-stages amenable to laboratory culture, e.g. free living species such as C. elegans or free-living stages of parasites such as PPN J2 larvae, and more recently in vivo in mosquito-stage Brugia [78], although some readily-cultured species seem refractory to RNAi [12]. Additionally, given that small non-coding RNAs are heavily involved in various cellular stress responses [79], it may be that adverse culture conditions lead to their increased expression, resulting in saturation of available RISC proteins, which would interfere with the organism's ability to direct an RNAi response to an exogenous trigger. If such saturation events varied between cells and/or tissues, then this could account for differing knockdown susceptibilities between some genes. Further, we have little information on differences in RNAi effector protein expression level or localization between species and/or life-stages, which might account for the observed variability. Other possible explanations for RNAi disparities include factors for which we have limited information, such as uncharacterized morphological differences between species (e.g. permeability of the cuticle to nucleic acids), or allelic diversity in discrete worm populations which may affect RNAi susceptibility in a similar fashion to drug susceptibility/resistance. In summary, our data do not support inter-species disparities in RNAi effector protein complements as an explanation for differences in RNAi competencies. Whilst the Caenorhabditid spp. encode significantly more RNAi pathway effectors than the other nematodes considered here, qualitative similarities in functional groupings across species with variable RNAi susceptibilities validate our conclusion.
Supporting Information
Dataset S1.
Nematode proteins associated with biosynthesis and nuclear export of small RNA; domains and sequence data. (*, putative stop codon)
https://doi.org/10.1371/journal.pntd.0001176.s001
(DOC)
Dataset S2.
Nematode proteins putatively responsible for secondary amplification, uptake and intercellular spread of siRNA; domains and sequence data. (*, putative stop codon)
https://doi.org/10.1371/journal.pntd.0001176.s002
(DOC)
Dataset S3.
Components of nematode RNA-induced silencing complex (RISC); domains and sequence data. (*, putative stop codon)
https://doi.org/10.1371/journal.pntd.0001176.s003
(DOC)
Dataset S4.
Nematode RNA interference (RNAi) inhibitor proteins; domains and sequence data. (*, putative stop codon)
https://doi.org/10.1371/journal.pntd.0001176.s004
(DOC)
Dataset S5.
Nematode chromatin modifiers, histone methylation factors, and other nuclear effectors; domains and sequence data. (*, putative stop codon)
https://doi.org/10.1371/journal.pntd.0001176.s005
(DOC)
Table S1.
Distribution of RNAi pathway components in nematode expressed-sequence tag, genome and transcriptome databases. Select tab at page bottom to view RISC-associated proteins, or all other proteins. Key to colours described below. Species were assigned to Clades according to nematode molecular phylogeny as described by [25], [26].
https://doi.org/10.1371/journal.pntd.0001176.s006
(XLS)
Acknowledgments
Thanks to Matt Berriman and John Gilleard for granting access to unpublished H. contortus RNAseq data. We would also like to acknowledge the contribution of the laboratories responsible for generating publically-available nematode sequence datasets, without whom this work would have been impossible.
Author Contributions
Conceived and designed the experiments: JJD PM NDW AGM. Performed the experiments: JJD PM NDW. Analyzed the data: JJD PM NDW AGM. Contributed reagents/materials/analysis tools: MM DMB PA. Wrote the paper: JJD PM NDW CCF TAD AM NJM AGM.
References
- 1. Liu Q, Paroo Z (2010) Biochemical principles of small RNA pathways. Annu Rev Biochem 79: 295–319.
- 2. Parkinson J, Mitreva M, Hall N, Blaxter M, McCarter JP (2003) 400000 Nematode ESTs on the net. Trends Parasitol 19: 283–286.
- 3. Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, et al. (2007) Draft genome of the filarial nematode parasite Brugia malayi. Science 317: 1756–1760.
- 4. Opperman CH, Bird DM, Williamson VM, Rokhsar DS, Burke M, et al. (2008) Sequence and genetic map of Meloidogyne hapla: A compact nematode genome for plant parasitism. Proc Natl Acad Sci USA 105: 14802–14807.
- 5. Abad P, Gouzy J, Aury JM, Castagnone-Sereno P, Danchin EG, et al. (2008) Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita. Nat Biotechnol 26: 909–915.
- 6. Urwin PE, Lilley CJ, Atkinson HJ (2002) Ingestion of double-stranded RNA by pre-parasitic juvenile cyst nematodes leads to RNA interference. Mol Plant Microbe Interact 15: 747–752.
- 7. Kimber MJ, McKinney S, McMaster S, Day TA, Fleming CC, et al. (2007) Flp gene disruption in a parasitic nematode reveals motor dysfunction and unusual neuronal sensitivity to RNA interference. FASEB J 21: 1233–1243.
- 8. Dalzell JJ, McMaster S, Fleming CC, Maule AG (2010) Short interfering RNA-mediated gene silencing in Globodera pallida and Meloidogyne incognita infective stage juveniles. Int J Parasitol 40: 91–100.
- 9. Huang G, Allen R, Davis EL, Baum TJ, Hussey RS (2006) Engineering broad root-knot resistance in transgenic plants by RNAi silencing of a conserved and essential root-knot nematode parasitism gene. Proc Natl Acad Sci USA 103: 14301–14306.
- 10. Yadav BC, Veluthambi K, Subramaniam K (2006) Host-generated double stranded RNA induces RNAi in plant-parasitic nematodes and protects the host from infection. Mol Biochem Parasitol 148: 219–222.
- 11. Geldhof P, Visser A, Clark D, Saunders G, Britton C, et al. (2007) RNA interference in parasitic helminths: Current situation, potential pitfalls and future prospects. Parasitology 134: 609–619.
- 12. Knox DP, Geldhof P, Visser A, Britton C (2007) RNA interference in parasitic nematodes of animals: A reality check? Trends Parasitol 23: 105–107.
- 13. Viney ME, Thompson FJ (2008) Two hypotheses to explain why RNA interference does not work in animal parasitic nematodes. Int J Parasitol 38: 43–47.
- 14. Visser A, Geldhof P, deMaere V, Knox DP, Vercruysse J, et al. (2006) Efficacy and specificity of RNA interference in larval life-stages of Ostertagia ostertagi. Parasitology 133: 777–783.
- 15. Geldhof P, Murray L, Couthier A, Gilleard JS, McLauchlan G, et al. (2006) Testing the efficacy of RNA interference in Haemonchus contortus. Int J Parasitol 36: 801–810.
- 16. Samarasinghe B, Knox DP, Britton C (2010) Factors affecting susceptibility to RNA interference in Haemonchus contortus and in vivo silencing of an H11 aminopeptidase gene. Int J Parasitol 41: 51–59.
- 17. Lendner M, Doligalska M, Lucius R, Hartmann S (2008) Attempts to establish RNA interference in the parasitic nematode Heligmosomoides polygyrus. Mol Biochem Parasitol 161: 21–31.
- 18. Louvet-Vallee S, Kolotuev I, Podbilewicz B, Felix MA (2003) Control of vulval competence and centering in the nematode Oscheius sp. 1 CEW1. Genetics 163: 133–146.
- 19. Pires da Silva A, Sommer RJ (2004) Conservation of the global sex determination gene tra-1 in distantly related nematodes. Gene Dev 18: 1198–1208.
- 20.
Pires da Silva A, WormBook, ed (2006) Pristionchus pacificus genetic protocols. The C. elegans Research Community, WormBook. doi/10.1895/wormbook.1.114.1.
- 21. Winston WM, Sutherlin M, Wright AJ, Feinberg EH, Hunter CP (2007) Caenorhabditis elegans SID-2 is required for environmental RNA interference. Proc Natl Acad Sci USA 104: 10565–10570.
- 22. Blackman MJ (2003) RNAi in protozoan parasites: What hope for the apicomplexa? Protist 154: 177–180.
- 23. Brown AE, Catteruccia F (2006) Toward silencing the burden of malaria: Progress and prospects for RNAi-based approaches. BioTechniques Apr: Suppl38–44.
- 24. Baum J, Papenfuss AT, Mair GR, Janse CJ, Vlachou D, et al. (2009) Molecular genetics and comparative genomics reveal RNAi is not functional in malaria parasites. Nucleic Acids Res 37: 3788–3798.
- 25. Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, et al. (1998) A molecular evolutionary framework for the phylum Nematoda. Nature 392: 71–75.
- 26. Holterman M, van der Wurff A, van den Elsen S, van Megen H, Bongers T, et al. (2006) Phylum-wide analysis of SSU rDNA reveals deep phylogenetic relationships among nematodes and accelerated evolution toward crown Clades. Mol Biol Evol 23: 1792–1800.
- 27. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403–410.
- 28. Mbeunkui F, Scholl EH, Opperman CH, Goshe MB, Bird DMcK (2010) Proteomic and bioinformatic analysis of the root-knot nematode Meloidogyne hapla: The basis for plant parasitism. J Proteome Res 9: 5370–5381.
- 29. Martin J, Abubucker S, Wylie T, Yin Y, Wang Z, et al. (2009) Nematode.net update 2008: improvements enabling more efficient data mining and comparative nematode genomics. Nucleic Acids Res 37: D571–578.
- 30. Parra G, Bradnam K, Ning Z, Keane T, Korf I (2008) Assessing the gene space in draft genomes. Nucleic Acids Res 37: 289–297.
- 31. Wang Z, Abubucker S, Martin J, Wilson RK, Hawdon J, et al. (2010) Characterising Ancylostoma caninum transcriptome and exploring nematode parasitic adaptation. BMC Genomics 11: 307.
- 32. Yigit E (2006) Analysis of the C. elegans argonaute family reveals that distinct argonautes act sequentially during RNAi. Cell 127: 747–757.
- 33. Winston WM, Molodowitch C, Hunter CP (2002) Systemic RNAi in C. elegans requires the putative transmembrane protein SID-1. Science 295: 2456–2459.
- 34. Duxbury MS, Ashley SW, Whang EE (2005) RNA interference: A mammalian SID-1 homologue enhances siRNA uptake and gene silencing efficacy in human cells. Biochem Biophys Res Commun 331: 459–463.
- 35. Krautz-Peterson G, Bhardwaj R, Faghiri Z, Tararam CA, Skelly PJ (2010) RNA interference in schistosomes: Machinery and methodology. Parasitology 137: 485–495.
- 36. Labreuche Y, Veloso A, de la Vega E, Gross PS, Chapman RW, et al. (2010) Non-specific activation of antiviral immunity and induction of RNA interference may engage the same pathway in the pacific white leg shrimp Litopenaeus vannamei. Dev Comp Immunol 34: 1209–1218.
- 37. Tomoyasu Y, Miller SC, Tomita S, Schoppmeier M, Grossman D, et al. (2008) Exploring systemic RNA interference in insects: A genome-wide survey for RNAi genes in Tribolium. Genome Biol 9: R10.
- 38. Xu W, Han Z (2008) Cloning and phylogenetic analysis of sid-1-like genes from aphids. J Insect Sci 8: 1–6.
- 39. Feinberg EH, Hunter CP (2003) Transport of dsRNA into cells by the transmembrane protein SID-1. Science 301: 1545–1547.
- 40. Calixto A, Chelur D, Topalidou I, Chen X, Chalfie M (2010) Enhanced neuronal RNAi in C. elegans using SID-1. Nat Methods 7: 554–559.
- 41. Tijsterman M, May RC, Simmer F, Okihara KL, Plasterk RH (2004) Genes required for systemic RNA interference in Caenorhabditis elegans. Curr Biol 14: 111–116.
- 42. Rosso MN, Jones JT, Abad P (2009) RNAi and functional genomics in plant parasitic nematodes. Annu Rev Phytopathol 47: 207–232.
- 43. Dalzell JJ, Warnock ND, Stevenson MA, Mousley A, Fleming CC, et al. (2010) Short interfering RNA-mediated knockdown of drosha and pasha in undifferentiated Meloidogyne incognita eggs leads to irregular growth and embryonic lethality. Int J Parasitol 40: 1303–1310.
- 44. Britton A, Murray L (2006) Using Caenorhabditis elegans for functional analysis of genes of parasitic nematodes. Int J Parasitol 36: 651–659.
- 45. Saleh MC, van Rij RP, Hekele A, Gillis A, Foley E, et al. (2006) The endocytic pathway mediates cell entry of dsRNA to induce RNAi silencing. Nat Cell Biol 8: 793–802.
- 46. Dalzell JJ, McMaster S, Johnston MJ, Kerr R, Fleming CC, et al. (2009) Non-nematode-derived double-stranded RNAs induce profound phenotypic changes in Meloidogyne incognita and Globodera pallida infective juveniles. Int J Parasitol 39: 1503–1516.
- 47. Mourrain P, Béclin C, Elmayan T, Feuerbach F, Gordon C, et al. (2000) Arabidopsis SGS2 and SGS3 genes are required for posttranscriptional gene silencing and natural virus resistance. Cell 101: 533–542.
- 48. Cogoni C, Macino G (1999) Gene silencing in Neurospora crassa requires a protein homologous to RNA-dependent RNA polymerase. Nature 399: 166–169.
- 49. Sijen T, Fleenor J, Simmer F, Thijssen KL, Parrish S, et al. (2001) On the role of RNA amplification in dsRNA-triggered gene silencing. Cell 107: 465–476.
- 50. Aoki K, Moriguchi H, Yoshioka T, Okawa K, Tabara H (2007) In vitro analyses of the production and activity of secondary small interfering RNAs in C. elegans. EMBO J 26: 5007–5019.
- 51. Pak J, Fire A (2007) Distinct populations of primary and secondary effectors during RNAi in C. elegans. Science 315: 241–244.
- 52. Sijen T, Steiner FA, Thijssen KL, Plasterk RH (2007) Secondary siRNAs result from unprimed RNA synthesis and form a distinct class. Science 315: 244–247.
- 53. Marker S, Le Mouel A, Meyer E, Simon M (2010) Distinct RNA-dependent RNA polymerases are required for RNAi triggered by double-stranded RNA versus truncated transgenes in Paramecium tetraurelia. Nucleic Acids Res 38: 4092–4107.
- 54. Lipardi C, Paterson BM (2009) Identification of an RNA-dependent RNA polymerase in Drosophila involved in RNAi and transposon suppression. Proc Natl Acad Sci USA 106: 15645–15650.
- 55. Gu W, Shirayama M, Conte D Jr, Vasale J, Batista PJ, et al. (2009) Distinct argonaute-mediated 22G-RNA pathways direct genome surveillance in the C. elegans germline. Mol Cell 36: 231–244.
- 56. Smardon A, Spoerke JM, Stacey SC, Klein ME, Mackin N, et al. (2000) EGO-1 is related to RNA-directed RNA polymerase and functions in germ-line development and RNA interference in C. elegans. Curr Biol 10: 169–178.
- 57. Vasale JJ, Gu W, Thivierge C, Batista PJ, Claycomb JM, et al. (2010) Sequential rounds of RNA-dependent RNA transcription drive endogenous small-RNA biogenesis in the ERGO-1/Argonaute pathway. Proc Natl Acad Sci USA 107: 3582–3587.
- 58. Vought VE, Ohmachi M, Lee MH, Maine EM (2005) EGO-1, a putative RNA-directed RNA polymerase, promotes germline proliferation in parallel with GLP-1/notch signaling and regulates the spatial organization of nuclear pore complexes and germline P granules in Caenorhabditis elegans. Genetics 170: 1121–1132.
- 59. Maine EM, Hauth J, Ratliff T, Vought VE, She X, et al. (2005) EGO-1, a putative RNA-dependent RNA polymerase, is required for heterochromatin assembly on unpaired DNA during C. elegans meiosis. Curr Biol 15: 1972–1978.
- 60. She X, Xu X, Fedotov A, Kelly WG, Maine EM (2009) Regulation of heterochromatin assembly on unpaired chromosomes during Caenorhabditis elegans meiosis by components of a small RNA-mediated pathway. PLoS Genet 5: e1000624.
- 61. Claycomb JM, Batista PJ, Pang KM, Gu W, Vasale JJ, et al. (2009) The argonaute CSR-1 and its 22G-RNA cofactors are required for holocentric chromosome segregation. Cell 139: 123–134.
- 62. Updike DL, Strome S (2009) A genomewide RNAi screen for genes that affect the stability, distribution and function of P granules in Caenorhabditis elegans. Genetics 183: 1397–1419.
- 63. Simmer F, Tijsterman M, Parrish S, Koushika SP, Nonet ML, et al. (2002) Loss of the putative RNA-directed RNA polymerase RRF-3 makes C. elegans hypersensitive to RNAi. Curr Biol 12: 1317–1319.
- 64. Han T, Manoharan AP, Harkins TT, Bouffard P, Fitzpatrick C, et al. (2009) 26G endo-siRNAs regulate spermatogenic and zygotic gene expression in Caenorhabditis elegans. Proc Natl Acad Sci USA 106: 18674–18679.
- 65. Gent JI, Lamm AT, Pavelec DM, Maniar JM, Parameswaran P, et al. (2010) Distinct phases of siRNA synthesis in an endogenous RNAi pathway in C. elegans soma. Mol Cell 37: 679–689.
- 66. Mango SE (2001) Stop making nonSense: the C. elegans smg genes. Trends Genet 17: 646–653.
- 67. Kim JK, Gabel HW, Kamath RS, Tewari M, Pasquinelli A, et al. (2005) Functional genomic analysis of RNA interference in C. elegans. Science 308: 1164–1167.
- 68. Conine CC, Batista PJ, Gu W, Claycomb JM, Chaves DA, et al. (2010) Argonautes ALG-3 and ALG-4 are required for spermatogenesis-specific 26G-RNAs and thermotolerant sperm in Caenorhabditis elegans. Proc Natl Acad Sci USA 107: 3588–3593.
- 69. Wasmuth J, Schmid R, Hedley A, Blaxter M (2008) On the extent and origins of genetic novelty in the phylum Nematoda. Plos Negl Trop Dis 2: e258.
- 70. Caudy AA, Ketting RF, Hammond SM, Denli AM, Bathoorn AM, et al. (2003) A micrococcal nuclease homologue in RNAi effector complexes. Nature 425: 411–414.
- 71. Ding L, Spencer A, Morita K, Han M (2005) The developmental timing regulator AIN-1 interacts with miRISCs and may target the argonaute protein ALG-1 to cytoplasmic P bodies in C. elegans. Mol Cell 19: 437–447.
- 72. Zhang L, Ding L, Cheung TH, Dong MQ, Chen J, et al. (2007) Systematic identification of C. elegans miRISC proteins, miRNAs, and mRNA targets by their interactions with GW182 proteins AIN-1 and AIN-2. Mol Cell 28: 598–613.
- 73. Chan SP, Ramaswamy G, Choi EY, Slack FJ (2008) Identification of specific let-7 microRNA binding complexes in Caenorhabditis elegans. RNA 14: 2104–2114.
- 74. Kennedy S, Wang D, Ruvkun G (2004) A conserved siRNA-degrading RNase negatively regulates RNA interference in C. elegans. Nature 427: 645–649.
- 75. Knight SW, Bass BL (2002) The role of RNA editing by ADARs in RNAi. Mol Cell 10: 809–817.
- 76. Pavelec DM, Lachowiec J, Duchaine TF, Smith HE, Kennedy S (2009) Requirement for the ERI/DICER complex in endogenous RNA interference and sperm development in Caenorhabditis elegans. Genetics 183: 1283–1295.
- 77. Rocheleau CE, Cullison K, Huang K, Bernstein Y, Spilker AC, et al. (2008) The Caenorhabditis elegans ekl (enhancer of ksr-1 lethality) genes include putative components of a germline small RNA pathway. Genetics 178: 1431–1443.
- 78. Song C, Gallup JM, Day TA, Bartholomay LC, Kimber MJ (2010) Development of an In Vivo RNAi Protocol to Investigate Gene Function in the Filarial Nematode, Brugia malayi. PLoS Pathogens 6: e1001239.
- 79. Shukla LI, Chinnusamy V, Sunkar R (2008) The role of microRNAs and other endogenous small RNAs in plant stress responses. Biochim Biophys Acta Gene Regul Mech 1779: 743.748