<title>Background</title> <p>RNA interference (RNAi) was first described in the nematode <italic>Caenorhabditis elegans </italic>[##REF##9486653##1##,##REF##9860998##2##], when it was found that the injection of double-stranded RNA (dsRNA) into a hermaphrodite worm resulted in the degradation of endogenous mRNA corresponding in sequence to the injected dsRNA. This resulted in a loss of function phenotype for the target gene. Potent RNAi effects can also be achieved in <italic>C. elegans </italic>by soaking the nematodes in a dsRNA solution [##REF##9841401##3##] or by feeding the nematodes on bacteria engineered to express dsRNA of the target gene [##REF##9804418##4##]. The potential for RNAi by feeding to be used in a large scale manner for functional genomics was demonstrated by the construction of a feeding library corresponding to 86% of the 19,000 predicted genes in <italic>C. elegans </italic>[##REF##12828945##5##]. This library has already been screened for genes involved in a wide range of biological processes including development [##REF##11483502##6##,##REF##16816419##7##], ageing [##REF##12845331##8##, ####REF##15998808##9##, ##REF##17411345##10####17411345##10##], cell biology [##REF##15272318##11##,##REF##17189861##12##] and neurobiology [##REF##17213328##13##].</p> <p>RNAi effects have now been detected in all eukaryote groups; phylogenetic analyses suggest that the main components of the RNAi machinery can be traced back to the last common ancestor of the eukaryotes [##REF##16691418##14##]. Gene silencing by RNAi is the outcome of a complex co-ordinated process. Many of the essential proteins involved in this process have been identified. RNAi is initiated by the cleavage of exogenous dsRNA into short interfering RNAs (siRNAs) by a dsRNA-specific RNAse enzyme Dicer [##REF##11201747##15##]. These double-stranded ~21 nucleotide siRNAs interact with, and activate, a multiprotein RNA-induced silencing complex (RISC) and the activated RISC then recognises and cleaves RNA molecules homologous to the interfering siRNA [##REF##10749213##16##,##REF##11701122##17##]. Gene silencing via RNAi can also be triggered by endogenous double stranded micro RNAs (miRNAs) which regulate gene expression at the post-transcriptional level in metazoans [##REF##11081512##18##] and plants [##REF##16669754##19##]. RNAi-related mechanisms are also involved in heterochromatin formation [##REF##12193640##20##,##REF##16009128##21##].</p> <p>Although the core features of the RNAi response are conserved among eukaryotes, several taxon-specific variations occur. Taxa differ in their ability to take up dsRNA, to transmit siRNA systemically to other tissues and in the magnitude and duration of the RNAi response. While various cell types can efficiently take up dsRNA in <italic>Drosophila melanogaster </italic>[##REF##10823906##22##], the response is cell autonomous and the flies are unable to produce a systemic response to the localized introduction of dsRNA [##REF##12592004##23##]. Naturally occurring systemic transmission of siRNAs also appears to be absent from vertebrates. Other organisms including <italic>C. elegans </italic>[##REF##9486653##1##], the insect parasitic nematode <italic>Heterorhabditis bacteriophora </italic>[##REF##17803822##24##], the planarian worm <italic>Schmidtea mediterranea </italic>[##REF##10220416##25##], the coelenterate <italic>Hydra </italic>[##REF##10491269##26##], insect species including a grasshopper [##REF##16202143##27##] and cockroach [##REF##16403020##28##] and several plant species [##REF##12181491##29##, ####REF##14555775##30##, ##REF##16242131##31##, ##REF##17287115##32####17287115##32##] have the ability to generate a systemic RNAi response. In <italic>C. elegans </italic>two genes <italic>sid-1 </italic>and <italic>sid-2 </italic>are required for the uptake of dsRNA from the intestine (sid = <underline>s</underline>ystemic RNA<underline>i</underline>-<underline>d</underline>eficient) [##REF##11834782##33##,##REF##17563372##34##]. SID-1 is a multispan transmembrane protein which appears to function as a channel for the diffusion of dsRNA [##REF##11834782##33##,##REF##12970568##35##], and expression of nematode SID-1 in <italic>Drosophila </italic>cells which lack a <italic>sid-1 </italic>homolog resulted in the rapid import of dsRNA into the cells [##REF##12970568##35##]. Although <italic>sid-2 </italic>mutants fail to show a systemic RNAi response to ingested dsDNA, the injection of dsRNA into the pseudocoel of these mutants generated a strong systemic RNAi response in the treated worms [##REF##17563372##34##]. SID-2 encodes a single-pass transmembrane protein which is expressed in the membranes of the intestine lumen and although its mode of action has yet to be determined, SID-2 is an essential requirement for environmental RNAi in <italic>C. elegans </italic>[##REF##17563372##34##].</p> <p>Because of the complexity of the dsRNAi pathway, the loss of a specific RNAi component protein or the occurrence of mutations may lead to the lack of a systemic RNAi response. Exogenous dsRNAi function has been lost from specific lineages such as <italic>Saccharomyces cerevisiae </italic>[##REF##16691418##14##] and <italic>Trypanosoma cruzi </italic>[##REF##14698430##36##]. Lack of both a systemic and a cell-autonomous RNAi response to injected dsRNA was also observed in the free living nematodes <italic>Oscheius tipulae </italic>[##REF##12586702##37##,##UREF##0##38##] and <italic>Pristionchus pacificus </italic>[##REF##15155582##39##,##REF##18050490##40##]. Testing of strains within the <italic>Caenorhabditis </italic>clade for susceptibility to RNAi by different delivery methods indicated that most strains are capable of a systemic RNAi response if dsRNA is delivered via microinjection, however six of the eight <italic>Caenorhabditis </italic>species tested, including <italic>C. briggsae</italic>, were resistant to dsRNAi via soaking or feeding [##REF##17563372##34##,##UREF##1##41##]. When SID-2 from <italic>C. elegans </italic>was expressed in <italic>C. briggsae </italic>it conferred sensitivity to RNAi by soaking [##REF##17563372##34##]. A naturally occurring strain of <italic>C. elegans </italic>from Hawaii was found to be resistant to RNAi by feeding or injection of germline-expressed genes and this defect was the result of mutations in the Argonaute gene <italic>ppw-1 </italic>[##REF##12225671##42##]. The minimum defining feature of a RISC complex is that it contains one Argonaute family member and a small-RNA guide strand [##REF##17173028##43##]. The Argonaute gene family contains 20 members in <italic>C. elegans </italic>[##REF##17029813##44##], which suggests that several distinct RISC complexes may be formed in this organism. Resistance or partial resistance to RNAi has also been demonstrated in animal parasitic nematode species such as <italic>Haemonchus contortus </italic>[##REF##16469321##45##] and <italic>Ostertagia ostertagi </italic>[##REF##16879764##46##]. This resistance to RNAi may be due to difficulties in developing suitable <italic>in vitro </italic>culture conditions and delivery methods for these parasites, but the available data suggest that the RNAi pathway may not be fully functional in some animal parasitic nematodes [##REF##17276139##47##]. By contrast, several plant parasitic nematode species are susceptible to RNAi, including the nematodes <italic>Heterodera glycines </italic>and <italic>Globodera pallida </italic>[##REF##12182331##48##] and <italic>Meloidogyne incognita </italic>[##REF##16255249##49##]. Viney and Thompson [##REF##18028931##50##] propose that the limited success of RNAi in animal parasitic nematodes to date might be circumvented if, instead of the soaking protocols currently in use, the dsRNA were delivered to the worms' pseudocoel. Alternatively it may be necessary to genetically engineer animal parasitic nematodes to express the <italic>sid-1 </italic>and <italic>sid-2 </italic>genes from <italic>C. elegans </italic>[##REF##18028931##50##].</p> <p>Native <italic>Drosophila </italic>S2 cells lack a SID-1 homolog, however these cells become competent to take up dsRNA when exposed to brief serum starvation in the presence of dsRNA, followed by incubation in medium containing serum [##REF##10823906##22##,##REF##11752672##51##]. Systemic RNAi effects in <italic>C. elegans </italic>body wall muscles are also enhanced by starvation [##REF##11834782##33##]. Observations in <italic>C. elegans </italic>and <italic>Drosophila </italic>suggest that environmental or physiological signals may be important in the uptake and spread of dsRNA. Using an RNAi approach Saleh <italic>et al. </italic>[##REF##16862146##52##] have shown that dsRNA is taken up in serum starved <italic>Drosophila </italic>S2 cells by receptor mediated endocytosis and these authors demonstrated that a functional intracellular vesicle transport system is also required for systemic RNAi in <italic>C. elegans</italic>. It has been hypothesized that an RNAi response to exogenous dsRNA evolved as a defence against infectious transposons and viruses [##REF##12016302##53##,##REF##16369574##54##]. Many viruses encode viral suppressors of RNA silencing (VSRs) to counteract the hosts' defence against viral encoded dsRNAs [##REF##16768647##55##]. VSRs target several of the key steps in the RNA silencing pathway and VSRs encoded by viruses of different families often share no homology at the primary amino acid sequence level [##REF##16768647##55##]. Pharmacological evidence suggests that in serum starved <italic>Drosophila </italic>S2 cells dsRNA uptake is mediated by pattern recognition receptors of the scavenger-receptor family [##REF##16862146##52##], a family of receptors with roles in the innate immune response [##REF##11790542##56##]. The ability to mount a systemic RNAi response to ingested dsRNA is absent in many <italic>Caenorhabditis </italic>species and it is also poorly developed in animal parasitic nematodes, and both systemic and cell autonomous RNAi are absent in the nematodes <italic>P. pacificus </italic>and <italic>O. tipulae</italic>. In contrast <italic>C. elegans, H. bacteriophora </italic>and several species of plant parasitic nematodes are susceptible to RNAi from environmentally supplied dsRNA. This phenotypic variability among nematodes may result from differing niche-specific selective pressures affecting environmental RNAi and the innate immune response.</p> <p>We have previously shown that the free living bacterial feeding nematode <italic>Panagrolaimus superbus </italic>is a suitable model for studies on the molecular basis of anhydrobiosis because of its remarkable ability to survive extreme desiccation [##REF##15939782##57##]. As a bacterial feeder, <italic>Panagrolaimus </italic>can be readily cultured in the laboratory with <italic>E. coli </italic>as a food source, using methods developed for <italic>C. elegans </italic>[##REF##4366476##58##]. We therefore investigated whether <italic>Panagrolaimus </italic>is susceptible to RNAi by feeding. Our results show that potent and specific RNAi effects can be achieved in <italic>Panagrolaimus </italic>by feeding it on <italic>E. coli </italic>expressing dsRNA for genes with essential conserved functions. This demonstration that <italic>Panagrolaimus </italic>is amenable to RNAi by feeding will allow the development of high throughput methods of RNAi screening for <italic>P. superbus </italic>and further enhance the utility of this nematode as a model system for the study of the molecular biology of anhydrobiosis, and as a possible satellite model nematode for comparative and functional genomics.</p>
<title>Methods</title> <title>Source and culturing of <italic>Panagrolaimus sp</italic>. and <italic>Caenorhabditis elegans</italic></title> <p>The nematode strains <italic>C. elegans </italic>N2, and <italic>Panagrolaimus </italic>sp. PS1159 were obtained from the <italic>Caenorhabditis </italic>Genetics Centre (CGC). <italic>P. superbus </italic>(strain DF5050) was obtained from Prof. Bjorn Sohlenius, Swedish Museum of Natural History, Stockholm. The nematodes were cultured on NGM agar plates [##REF##15939782##57##].</p> <title><italic>E. coli </italic>feeding clones</title> <p>Two <italic>E. coli </italic>RNAi feeding clones from the <italic>C. elegans </italic>chromosome I feeding library [##REF##11099033##59##] were kindly provided by Dr. Julie Ahringer, Department of Genetics, University of Cambridge. These clones contained inserts from the embryonic lethal genes <italic>Ce-ran-4 </italic>(R05D11.3) and <italic>Ce-lmn-1 </italic>(DY3.2). Bacterial feeding strains were constructed for the <italic>Panagrolaimus </italic>genes using the L4440 plasmid vector [##REF##9804418##4##] and <italic>Escherichia coli </italic>HT115(DE3) cells as described by Kamath and Ahringer [##REF##12828945##5##]. The <italic>Panagrolaimus </italic>genes were selected from a panel of expressed sequence tags (ESTs) which we have generated for <italic>P. superbus </italic>and <italic>Panagrolaimus </italic>sp. PS1159 (Tyson <italic>et al. </italic>unpublished). The target genes were selected based on their gene identity as determined by BLAST analysis [##UREF##7##86##] and the RNAi phenotypes of their orthologs in <italic>C. elegans</italic>, as listed in WormBase [##UREF##8##87##]. The accession numbers of the <italic>Panagrolaimus </italic>ESTs and their cDNA fragment sizes are as follows: <italic>ef1b</italic>, EU368942 (625 bp); <italic>rps-2</italic>, EU368941 (649 bp); <italic>ifb-1</italic>, EU368943 (360 bp) and <italic>act-2</italic>, EU368944 (391 bp). Gene fragments of the <italic>ef1b </italic>and <italic>rps-2 </italic>genes for insertion into the L4440 plasmid were prepared by PCR using the gene specific primer pairs EF1 Forward/EF1 Reverse and RPS Forward/RPS Reverse (The sequences of these PCR primers are listed in Additional File ##SUPPL##0##1##). Gene fragments corresponding to the entire EST sequences of the <italic>ifb-1 </italic>and <italic>act-2 </italic>genes were prepared from the plasmid name using M13 F and R primers and cloned into the L4440 vector. A separate control experiment (data not presented) showed that an RNAi feeding clone containing dsRNA for the multiple cloning site of the pDNR-LIB plasmid had no observable phenotype when fed to <italic>P. superbus </italic>or <italic>Panagrolaimus </italic>sp. PS1159.</p> <title>RNAi feeding protocol</title> <p>NGM agar plates (9 cm) containing 1 mM iso-propyl-β-D-thiogalactopyranoside (IPTG), 50 ug/ml ampicillin (Amp) and 12.5 ug/ml tetracycline (Tet) were prepared. The plates were allowed to dry inverted for ~4 days at room temperature. They were then were spread with cells from an overnight <italic>E. coli </italic>culture (strain HTT1(DE3)) containing the empty L4440 plasmid (control) or with <italic>E. coli </italic>HTT1(DE3) transformed with an L4440 plasmid containing the target gene fragment inserted into the multiple cloning site, and flanked on both sides by an IPTG-inducible T7 RNA polymerase promoter [##REF##9804418##4##]. The plates were grown at 20°C overnight to allow a lawn of <italic>E. coli </italic>to form and to begin induction and expression of dsRNA. Axenised eggs were prepared from gravid female nematodes by bleaching them with hypochlorite [##UREF##9##88##]. The axenised eggs were rinsed three times in S Basal buffer [##REF##4366476##58##], then aliquotted onto NGM agar feeding plates. The plates were incubated at 20°C and monitored daily for the presence of visible RNAi phenotypes in the treated nematodes. To analyse embryonic lethal RNAi phenotypes young L4 stage worms were rinsed from the treatment and control plates using sterile S Basal buffer and approximately 5 worms were transferred to each of 5–8 replicate 3 cm NGM/IPTG/Amp/Tet agar plates containing a lawn of either control <italic>E. coli </italic>or <italic>E. coli </italic>expressing dsRNA for the target gene. These worms were allowed to grow to adulthood and lay eggs for 24 hours and they were then transferred to a second set of replicate plates. The worms were again allowed to lay eggs for 24 hours before being transferred to a third set of plates. The worms were allowed to lay eggs on the third set of plates until laying ceased and the adult worms were then removed from the plates. The number of dead embryos on each of the three sets of plates was determined by microscopic observation over a number of days for both the control and experimental treatments. Standard deviation and standard errors were calculated using the individual counts of surviving nematodes and not on proportional data, but in the histograms survival data are presented as percentage survival.</p> <p>To observe possible embryonic lethal RNAi phenotypes synchronous eggs from control and RNAi treated cultures were rinsed from the agar surface using S Basal buffer and centrifuged at 3,000 g for five minutes. The concentrated eggs were transferred to a glass slide and gently covered with a glass cover slip. Embryos were observed using an Olympus EX51 microscope fitted with a U-CMAD3 imaging system and an Olympus U-RFL-T camera.</p> <title>RT-PCR analysis of gene expression levels</title> <p>Individual adult female nematodes whose progeny displayed an embryonic lethal RNAi phenotype were used for RT-PCR. RNA was isolated from single nematodes using the Absolutely RNA Microprep Kit (Stratagene), which also included a DNAse treatment prior to the reverse transcription step to remove any genomic DNA (gDNA) contamination. As a control, PCR was attempted on each sample after DNAse treatment and prior to RT-PCR to ensure all gDNA had been removed. PCR samples were then analysed by spectrophotometry and agarose gel electrophoresis and samples showing no gDNA amplification were chosen for RT-PCR experiments. Each worm RNA preparation yielded 30 μl of extracted RNA. Six μl of the extracted RNA solution was used for each RT-PCR reaction. RT-PCR was carried out using the QuanTitect SYBR Green RT-PCR kit (Qiagen). The upstream test forward primers used for RT-PCR analysis targeted a region not included in the dsRNA feeding construct. The primers used were the EF1-F-TEST and the EF1 reverse primers for the <italic>ef1b </italic>gene and the RPS-F-TEST and the RPS reverse primers for the <italic>rps</italic>-2 gene. The sequences of all PCR primers are listed in Additional File ##SUPPL##0##1##. RT-PCR was also carried using the D3 expansion region of the ribosomal RNA gene as an internal control gene for each single worm RNA preparation. Amplification of the rRNA D3 domain was carried out as described previously [##REF##15939782##57##]. The RT-PCR conditions were as follows: 94°C for 5 min at the start of the reaction was followed by 30 minutes at 50°C to allow reverse transcription of the RNA. The PCR cycling conditions were as follows: 94°C, 1 min; 55°C, 1.5 min; 72°C, 2 min. To obtain PCR product during the linear phase of the PCR reaction a 10 μl aliquot from the PCR reaction mixture was sampled at 20, 25, 30 and 35 cycles for each gene used. The products of the RT-PCR reactions were visualised on 1% agarose gels with 10 μl of each RT-PCR product loaded per gel. Image analysis of the RT-PCR gels was carried out and band intensities were measured in machine units using the ImageQuant™ TL software package from GE Healthcare.</p> <title>Phylogenetic analyses</title> <p>Nucleotide sequences for the 18S small subunit rRNA nuclear gene of 34 taxa were downloaded from the NCBI database [##UREF##10##89##]. The accession numbers for these sequences are given in Additional File ##SUPPL##1##2##. The sequences were aligned using the ClustalW algorithm as implemented in the MEGA software package (version 3.1) [##REF##15260895##90##]. The alignment was manually edited by removing poorly aligned blocks at the 5' and 3' ends and ambiguously aligned internal segments. The final alignment contained 1,275 aligned characters. Phylogenetic trees were constructed using the Neighbor-Joining method [##REF##3447015##91##] and Minimum Evolution method [##UREF##11##92##] as implemented by MEGA version 3.1. The tree topologies recovered trees by both methods were very similar, differing only in the placement of the <italic>Oscheius </italic>clade (see Figure ##FIG##5##6## and Additional File ##SUPPL##2##3##).</p> <title>Statistical analyses</title> <p>All statistical analyses were carried out using Minitab™ version 14.1 (Minitab Inc., State College, Pennsylvania, USA). The 2-tailed T-Test was used to test for difference between populations. The null hypothesis that population means are all equal was tested at the <italic>P </italic>= 0.05 level of significance. T-Tests were carried out using the individual counts of surviving nematodes and not on proportional data.</p>
<title>Results and Discussion</title> <title><italic>Panagrolaimus </italic>fed with <italic>E. coli </italic>expressing dsRNA for <italic>C. elegans </italic>embryonic lethal genes give rise to progeny with highly penetrant embryonic lethal phenotypes</title> <p>We tested two species of <italic>Panagrolaimus </italic>for RNAi effects using <italic>E. coli </italic>clones from a <italic>C. elegans </italic>feeding library [##REF##11099033##59##]. The selected clones contained inserts from the genes <italic>Ce-lmn-1 </italic>(DY3.2) and <italic>Ce-ran-4 </italic>(R05D11.3). The functions of these genes are conserved among eukaryotes and they yield severe embryonic lethal RNAi phenotypes in <italic>C. elegans </italic>[##REF##11099033##59##]. The <italic>C. elegans </italic>genome contains a single lamin gene, <italic>lmn-1 </italic>which encodes a nuclear lamin protein [##REF##11071918##60##]. Previous RNAi studies have shown that <italic>C. elegans lmn-1 (RNAi) </italic>embryos have defects in chromatin organization, cell cycle progression and chromosome segregation [##REF##11071918##60##] and that embryogenesis is halted before the embryo is enclosed by epidermal tissue [##REF##11483502##6##]. The <italic>Ce-lmn-1 (RNAi) Panagrolaimus </italic>embryos show similar defects to those described for <italic>C. elegans </italic>and the terminal arrest phenotype is similar in both nematode genera (Figure ##FIG##0##1A##, ##FIG##0##1D##, ##FIG##0##1G##). <italic>Ce-ran-4 </italic>encodes an ortholog of nuclear transport factor 2 [##REF##10889207##61##] which is responsible for importing the Ras-related nuclear protein RanGDP into the nucleus [##REF##9822603##62##]. Ran is an essential component in the transport of many proteins and nucleic acids between the nucleus and the cytoplasm [##REF##9582265##63##]. RNAi of <italic>ran-4 </italic>in <italic>C. elegans </italic>resulted in an embryonic lethal phenotype in which the embryos arrested as a multicellular mass (&lt;200 cells) and the nuclear membranes were indistinct [##REF##11483502##6##,##REF##10889207##61##]. The terminal arrest phenotypes of <italic>Ce-ran-4 (RNAi) </italic>embryos in <italic>Panagrolaimus </italic>were also similar to those of <italic>C. elegans </italic>(Figure ##FIG##0##1B##, ##FIG##0##1E##, ##FIG##0##1H##). Although embryonic lethal RNAi phenotypes were observed for both species of <italic>Panagrolaimus</italic>, the RNAi effect showed greater penetrance in <italic>P. superbus </italic>than in <italic>Panagrolaimus </italic>sp. 1159 (Figure ##FIG##1##2##). From Figure ##FIG##1##2## it can also be seen that the penetrance of the heterologous RNAi effect observed in <italic>P. superbus </italic>for both dsRNA constructs was comparable to penetrance of the RNAi phenotypes observed in <italic>C. elegans </italic>which had been treated with homologous dsRNA constructs. All three species tested (<italic>C. elegans</italic>, <italic>Panagrolaimus </italic>sp. PS1159 and <italic>P. superbus</italic>) showed significant differences in the percentage dead embryos between the experimental and control populations for both target genes tested (Two sample T-Test, P = 0.00).</p> <title>Embryonic lethal RNAi phenotypes obtained by feeding <italic>Panagrolaimus </italic>with <italic>E. coli </italic>expressing dsRNA for <italic>Panagrolaimus </italic>genes</title> <p><italic>E. coli </italic>feeding strains were constructed using sequences for two genes from <italic>Panagrolaimus </italic>sp. PS1159 whose orthologs have embryonic lethal RNAi phenotypes in <italic>C. elegans</italic>. These were the <italic>elongation factor 1-beta (ef1b) </italic>gene and the ribosomal protein small subunit S2 (<italic>rps-2</italic>) gene. The <italic>Ce-rps-2 </italic>gene demonstrated both a sterile maternal RNAi phenotype [##REF##12529635##64##] and an embryonic lethal RNAi phenotype in <italic>C. elegans </italic>[##REF##12176360##65##]. The <italic>Panagrolaimus rps-2 </italic>gene also yielded an embryonic lethal RNAi phenotype in both <italic>Panagrolaimus </italic>species tested (Figure ##FIG##2##3A##), with a significant difference in percentage dead embryos between the experimental and control populations in both cases (Two sample T-Test, P = 0.00). Here also the heterologous RNAi effect in <italic>P. superbus </italic>(80% embryonic lethality) was stronger than the homologous RNAi effect in <italic>Panagrolaimus </italic>sp. PS1159 (38% embryonic lethality). The average number of eggs laid per female nematode was compared for the control and RNAi treated <italic>Panagrolaimus </italic>populations and a Two sample T-Test (p = 0.62) did not reveal any significant difference in fecundity between the populations. The dsRNAi embryonic lethal phenotype for the <italic>Panagrolaimus </italic>sp. PS1159 <italic>ef1b </italic>gene also showed a stronger heterologous RNAi effect in <italic>P. superbus </italic>(53% embryonic lethality) than the homologous RNAi effect in <italic>Panagrolaimus </italic>sp. PS1159 (39% embryonic lethality), however there was a significant difference in the percentage of dead embryos between the experimental and control populations in both nematode species (Two sample T-Test, P = 0.00).</p> <title>Exposure of <italic>Panagrolaimus </italic>sp. PS1159 to homologous dsRNA by feeding results in a reduction of mRNA transcript levels for the target genes</title> <p>Single nematode RT-PCR (Figure ##FIG##3##4##) showed that a reduction in mRNA transcript levels occurred for the target genes in the RNAi treated <italic>Panagrolaimus </italic>sp. PS1159 nematodes. For the <italic>ef1b </italic>gene the reduction in transcript level was 30% that of the controls while for the <italic>rps-2 </italic>gene the mRNA transcript levels were reduced to 55% of the controls (Figure ##FIG##3##4C, D##). There was no significant difference between the treated and control worms in the mRNA transcript levels for the D3 segment of the nuclear large subunit rRNA gene which was used as an internal control. This indicates that the dsRNAi effect in <italic>Panagrolaimus </italic>is specific for the targeted genes.</p> <title>Exposing <italic>Panagrolaimus superbus </italic>to homologous dsRNA for genes encoding contractile proteins results in potent and highly penetrant locomotor and growth defects</title> <p><italic>E. coli </italic>feeding strains were constructed for two genes from <italic>Panagrolaimus suprbus </italic>whose <italic>C. elegans </italic>orthologs encode contractile proteins. These genes were <italic>Psu-ifb-1 </italic>and <italic>Psu-act-2</italic>. Each of these genes has multiple RNAi phenotypes in <italic>C. elegans</italic>. The <italic>Ce-ifb-1 </italic>gene encodes two isoforms of a cytoplasmic intermediate filament protein that are expressed in the hypodermal cells of <italic>C. elegans</italic>, where they form attachment structures which connect the muscles to the cuticle [##REF##1860880##66##]. These attachment structures transmit the contractile forces generated by the body muscles across the hypodermis to the cuticle. In <italic>C. elegans </italic>both IFB-1 isoforms have roles in epidermal morphogenesis and muscle attachment [##REF##14975728##67##]. <italic>C. elegans </italic>larvae exposed to <italic>ifb-1 </italic>RNAi by feeding displayed a paralysis phenotype in which the hypodermis becomes detached from the cuticle [##REF##14529618##68##]. Exposure of <italic>P. superbus </italic>L1 larvae to <italic>Psu-ifb-1 </italic>dsRNA by feeding resulted in a larval arrest and paralysis phenotype (Figure ##FIG##4##5A##). The phenotype was highly penetrant and the treated larvae were severely paralysed. Microscopic observations revealed waves or ripples of muscle contraction moving from the anterior to the posterior of the immobile worms, indicating that the contractile forces generated by the muscles were not being transmitted to the cuticle.</p> <p>The <italic>C. elegans </italic>genome encodes five highly conserved actin genes. ACT-1, -2, -3 and -4 all share 99% amino acid identity, while ACT-5, the most divergent, shares 93% identity with the other actin isoforms [##REF##6527380##69##]. The <italic>act-1, -2, -3 </italic>and <italic>-4 </italic>genes function primarily in muscle and myofilament containing cells in <italic>C. elegans </italic>[##REF##6527380##69##,##REF##8406009##70##], while <italic>act-5 </italic>expression is restricted to cells which contain microvilli [##REF##15872090##71##]. A variety of mutant and RNAi phenotypes have been associated with impared actin gene function in <italic>C. elegans</italic>, including defects in embryonic cell division [##REF##16407404##72##], defects in muscle contractile function [##REF##6098683##73##] including pharyngeal function [##REF##8462849##74##], and growth defects [##REF##15872090##71##]. <italic>P. superbus </italic>larvae grown on <italic>E. coli </italic>expressing <italic>Psu-act-2 </italic>dsRNA displayed severe larval lethal, larval arrest, and paralysed phenotypes (Figure ##FIG##4##5B##). The treated larvae were unable to grow beyond the first or second larval stages. The worms were immobile, although pharyngeal function appeared to be retained.</p> <title><italic>Panagrolaimus</italic>: a potential satellite model nematode for functional genomics</title> <p>Comparative studies show that extensive genetic diversity occurs within the Phylum Nematoda [##UREF##2##75##,##REF##14624247##76##]. Consequently it is not always possible to extrapolate from the model nematode <italic>C. elegans </italic>to other members of the Phylum Nematoda; studies in additional representative nematode species are often necessary. Several additional species of plant and animal parasitic nematodes have been selected as model organisms. Comparative studies are also carried out on other species in the <italic>Caenorhabditis </italic>clade [##REF##14624247##76##,##REF##17240335##77##] and on the free-living bacterial feeding nematodes <italic>Oeschius tipulae </italic>[##UREF##0##38##,##UREF##3##78##] and <italic>Pristionchus pacificus </italic>[##REF##18050490##40##,##REF##16700067##79##]. Neither <italic>O. tipulae </italic>nor <italic>P. pacificus </italic>(both members of the suborder Rhabditina) display gene silencing in response to exogenous dsRNA. However injection of antisense morpholino DNA oligonucleotides into both of these satellite models results in a partial inactivation of the corresponding gene product [##REF##12586702##37##] or in phenotypes with low penetrance [##REF##15155582##39##]. Morpholino antisense oligonucleotides block mRNA translation when complementary to the 5' untranslated region of the target gene or to the first 25 bases downstream of the start codon, and they achieve their effect by preventing ribosomal binding [##REF##11884031##80##]. The use of morpholinos presents difficulties for high-throughput gene knockout experiments because of the need for specific sequence information at the 5' end of the gene, the requirement for a means of delivery of the oligonucleotides to the target cells, and the lack of a systemic response.</p> <p>We have previously shown that <italic>Panagrolaimus </italic>is a suitable model for studies on the molecular basis of anhydrobiosis [##REF##15939782##57##]. Panagrolaimid nematodes occupy a diversity of niches ranging from polar, temperate and semi-arid soils to terrestrial mosses. Thus many <italic>Panagrolaimus </italic>species are adapted to tolerate freezing and desiccation [##UREF##4##81##, ####UREF##5##82##, ##REF##12477892##83####12477892##83##]. Our current data show that when fed on <italic>E. coli </italic>expressing dsRNA targeting <italic>Panagrolaimus </italic>genes the nematode mounts a potent and specific systemic RNAi response. This makes <italic>Panagrolaimus </italic>an ideal system with which to carry out high-throughput RNAi screens for studies on the molecular basis of anhydrobiosis and cryobiosis. <italic>Panagrolaimus </italic>is a member of the suborder Tylenchina of the Order Rhabditida and thus it is more distantly related to <italic>C. elegans </italic>than <italic>Oeschius </italic>or <italic>Pristionchus </italic>(Figure ##FIG##5##6##). The Suborder Tylenchina contains three major lineages [##UREF##6##84##]: the Panagrolaimorpha which includes the free living bacterial feeding panagrolaimids as well as the zooparasitic strongyloids and steinernematids; the Cephalobomorpha, bacterial feeding nematodes that, like the panagrolaimaids, have a widespread distribution and are capable of surviving in extremely dry and cold environments and the Tylenchomorpha, a group of plant parasites and fungal feeders which frequently has anhydrobiotic encysted infective stages (e.g. <italic>Globodera</italic>, <italic>Meliodogyne</italic>), as well as species which can enter anhydrobiosis throughout most or all of their life cycle (e.g. <italic>Aphelencus avenae</italic>). Within the genus <italic>Panagrolaimus </italic>several changes in reproductive mode have occurred. Examples of gonochoristic (male female), hermaphroditic (self-fertilization), and parthenogenetic (development of unfertilized eggs) reproduction have been reported [##REF##10772792##85##]. Thus <italic>Panagrolaimus </italic>is an excellent model system to study the evolutionary origins and developmental consequences of parthenogenetic reproduction. Being closely related to the zooparasitic strongyloids and steinernematids and the plant parasitic tylenchid nematodes, <italic>Panagrolaimus </italic>could also be used in comparative genomic approaches to study the evolution of parasitic lifestyles.</p>
<title>Results and Discussion</title> <title><italic>Panagrolaimus </italic>fed with <italic>E. coli </italic>expressing dsRNA for <italic>C. elegans </italic>embryonic lethal genes give rise to progeny with highly penetrant embryonic lethal phenotypes</title> <p>We tested two species of <italic>Panagrolaimus </italic>for RNAi effects using <italic>E. coli </italic>clones from a <italic>C. elegans </italic>feeding library [##REF##11099033##59##]. The selected clones contained inserts from the genes <italic>Ce-lmn-1 </italic>(DY3.2) and <italic>Ce-ran-4 </italic>(R05D11.3). The functions of these genes are conserved among eukaryotes and they yield severe embryonic lethal RNAi phenotypes in <italic>C. elegans </italic>[##REF##11099033##59##]. The <italic>C. elegans </italic>genome contains a single lamin gene, <italic>lmn-1 </italic>which encodes a nuclear lamin protein [##REF##11071918##60##]. Previous RNAi studies have shown that <italic>C. elegans lmn-1 (RNAi) </italic>embryos have defects in chromatin organization, cell cycle progression and chromosome segregation [##REF##11071918##60##] and that embryogenesis is halted before the embryo is enclosed by epidermal tissue [##REF##11483502##6##]. The <italic>Ce-lmn-1 (RNAi) Panagrolaimus </italic>embryos show similar defects to those described for <italic>C. elegans </italic>and the terminal arrest phenotype is similar in both nematode genera (Figure ##FIG##0##1A##, ##FIG##0##1D##, ##FIG##0##1G##). <italic>Ce-ran-4 </italic>encodes an ortholog of nuclear transport factor 2 [##REF##10889207##61##] which is responsible for importing the Ras-related nuclear protein RanGDP into the nucleus [##REF##9822603##62##]. Ran is an essential component in the transport of many proteins and nucleic acids between the nucleus and the cytoplasm [##REF##9582265##63##]. RNAi of <italic>ran-4 </italic>in <italic>C. elegans </italic>resulted in an embryonic lethal phenotype in which the embryos arrested as a multicellular mass (&lt;200 cells) and the nuclear membranes were indistinct [##REF##11483502##6##,##REF##10889207##61##]. The terminal arrest phenotypes of <italic>Ce-ran-4 (RNAi) </italic>embryos in <italic>Panagrolaimus </italic>were also similar to those of <italic>C. elegans </italic>(Figure ##FIG##0##1B##, ##FIG##0##1E##, ##FIG##0##1H##). Although embryonic lethal RNAi phenotypes were observed for both species of <italic>Panagrolaimus</italic>, the RNAi effect showed greater penetrance in <italic>P. superbus </italic>than in <italic>Panagrolaimus </italic>sp. 1159 (Figure ##FIG##1##2##). From Figure ##FIG##1##2## it can also be seen that the penetrance of the heterologous RNAi effect observed in <italic>P. superbus </italic>for both dsRNA constructs was comparable to penetrance of the RNAi phenotypes observed in <italic>C. elegans </italic>which had been treated with homologous dsRNA constructs. All three species tested (<italic>C. elegans</italic>, <italic>Panagrolaimus </italic>sp. PS1159 and <italic>P. superbus</italic>) showed significant differences in the percentage dead embryos between the experimental and control populations for both target genes tested (Two sample T-Test, P = 0.00).</p> <title>Embryonic lethal RNAi phenotypes obtained by feeding <italic>Panagrolaimus </italic>with <italic>E. coli </italic>expressing dsRNA for <italic>Panagrolaimus </italic>genes</title> <p><italic>E. coli </italic>feeding strains were constructed using sequences for two genes from <italic>Panagrolaimus </italic>sp. PS1159 whose orthologs have embryonic lethal RNAi phenotypes in <italic>C. elegans</italic>. These were the <italic>elongation factor 1-beta (ef1b) </italic>gene and the ribosomal protein small subunit S2 (<italic>rps-2</italic>) gene. The <italic>Ce-rps-2 </italic>gene demonstrated both a sterile maternal RNAi phenotype [##REF##12529635##64##] and an embryonic lethal RNAi phenotype in <italic>C. elegans </italic>[##REF##12176360##65##]. The <italic>Panagrolaimus rps-2 </italic>gene also yielded an embryonic lethal RNAi phenotype in both <italic>Panagrolaimus </italic>species tested (Figure ##FIG##2##3A##), with a significant difference in percentage dead embryos between the experimental and control populations in both cases (Two sample T-Test, P = 0.00). Here also the heterologous RNAi effect in <italic>P. superbus </italic>(80% embryonic lethality) was stronger than the homologous RNAi effect in <italic>Panagrolaimus </italic>sp. PS1159 (38% embryonic lethality). The average number of eggs laid per female nematode was compared for the control and RNAi treated <italic>Panagrolaimus </italic>populations and a Two sample T-Test (p = 0.62) did not reveal any significant difference in fecundity between the populations. The dsRNAi embryonic lethal phenotype for the <italic>Panagrolaimus </italic>sp. PS1159 <italic>ef1b </italic>gene also showed a stronger heterologous RNAi effect in <italic>P. superbus </italic>(53% embryonic lethality) than the homologous RNAi effect in <italic>Panagrolaimus </italic>sp. PS1159 (39% embryonic lethality), however there was a significant difference in the percentage of dead embryos between the experimental and control populations in both nematode species (Two sample T-Test, P = 0.00).</p> <title>Exposure of <italic>Panagrolaimus </italic>sp. PS1159 to homologous dsRNA by feeding results in a reduction of mRNA transcript levels for the target genes</title> <p>Single nematode RT-PCR (Figure ##FIG##3##4##) showed that a reduction in mRNA transcript levels occurred for the target genes in the RNAi treated <italic>Panagrolaimus </italic>sp. PS1159 nematodes. For the <italic>ef1b </italic>gene the reduction in transcript level was 30% that of the controls while for the <italic>rps-2 </italic>gene the mRNA transcript levels were reduced to 55% of the controls (Figure ##FIG##3##4C, D##). There was no significant difference between the treated and control worms in the mRNA transcript levels for the D3 segment of the nuclear large subunit rRNA gene which was used as an internal control. This indicates that the dsRNAi effect in <italic>Panagrolaimus </italic>is specific for the targeted genes.</p> <title>Exposing <italic>Panagrolaimus superbus </italic>to homologous dsRNA for genes encoding contractile proteins results in potent and highly penetrant locomotor and growth defects</title> <p><italic>E. coli </italic>feeding strains were constructed for two genes from <italic>Panagrolaimus suprbus </italic>whose <italic>C. elegans </italic>orthologs encode contractile proteins. These genes were <italic>Psu-ifb-1 </italic>and <italic>Psu-act-2</italic>. Each of these genes has multiple RNAi phenotypes in <italic>C. elegans</italic>. The <italic>Ce-ifb-1 </italic>gene encodes two isoforms of a cytoplasmic intermediate filament protein that are expressed in the hypodermal cells of <italic>C. elegans</italic>, where they form attachment structures which connect the muscles to the cuticle [##REF##1860880##66##]. These attachment structures transmit the contractile forces generated by the body muscles across the hypodermis to the cuticle. In <italic>C. elegans </italic>both IFB-1 isoforms have roles in epidermal morphogenesis and muscle attachment [##REF##14975728##67##]. <italic>C. elegans </italic>larvae exposed to <italic>ifb-1 </italic>RNAi by feeding displayed a paralysis phenotype in which the hypodermis becomes detached from the cuticle [##REF##14529618##68##]. Exposure of <italic>P. superbus </italic>L1 larvae to <italic>Psu-ifb-1 </italic>dsRNA by feeding resulted in a larval arrest and paralysis phenotype (Figure ##FIG##4##5A##). The phenotype was highly penetrant and the treated larvae were severely paralysed. Microscopic observations revealed waves or ripples of muscle contraction moving from the anterior to the posterior of the immobile worms, indicating that the contractile forces generated by the muscles were not being transmitted to the cuticle.</p> <p>The <italic>C. elegans </italic>genome encodes five highly conserved actin genes. ACT-1, -2, -3 and -4 all share 99% amino acid identity, while ACT-5, the most divergent, shares 93% identity with the other actin isoforms [##REF##6527380##69##]. The <italic>act-1, -2, -3 </italic>and <italic>-4 </italic>genes function primarily in muscle and myofilament containing cells in <italic>C. elegans </italic>[##REF##6527380##69##,##REF##8406009##70##], while <italic>act-5 </italic>expression is restricted to cells which contain microvilli [##REF##15872090##71##]. A variety of mutant and RNAi phenotypes have been associated with impared actin gene function in <italic>C. elegans</italic>, including defects in embryonic cell division [##REF##16407404##72##], defects in muscle contractile function [##REF##6098683##73##] including pharyngeal function [##REF##8462849##74##], and growth defects [##REF##15872090##71##]. <italic>P. superbus </italic>larvae grown on <italic>E. coli </italic>expressing <italic>Psu-act-2 </italic>dsRNA displayed severe larval lethal, larval arrest, and paralysed phenotypes (Figure ##FIG##4##5B##). The treated larvae were unable to grow beyond the first or second larval stages. The worms were immobile, although pharyngeal function appeared to be retained.</p> <title><italic>Panagrolaimus</italic>: a potential satellite model nematode for functional genomics</title> <p>Comparative studies show that extensive genetic diversity occurs within the Phylum Nematoda [##UREF##2##75##,##REF##14624247##76##]. Consequently it is not always possible to extrapolate from the model nematode <italic>C. elegans </italic>to other members of the Phylum Nematoda; studies in additional representative nematode species are often necessary. Several additional species of plant and animal parasitic nematodes have been selected as model organisms. Comparative studies are also carried out on other species in the <italic>Caenorhabditis </italic>clade [##REF##14624247##76##,##REF##17240335##77##] and on the free-living bacterial feeding nematodes <italic>Oeschius tipulae </italic>[##UREF##0##38##,##UREF##3##78##] and <italic>Pristionchus pacificus </italic>[##REF##18050490##40##,##REF##16700067##79##]. Neither <italic>O. tipulae </italic>nor <italic>P. pacificus </italic>(both members of the suborder Rhabditina) display gene silencing in response to exogenous dsRNA. However injection of antisense morpholino DNA oligonucleotides into both of these satellite models results in a partial inactivation of the corresponding gene product [##REF##12586702##37##] or in phenotypes with low penetrance [##REF##15155582##39##]. Morpholino antisense oligonucleotides block mRNA translation when complementary to the 5' untranslated region of the target gene or to the first 25 bases downstream of the start codon, and they achieve their effect by preventing ribosomal binding [##REF##11884031##80##]. The use of morpholinos presents difficulties for high-throughput gene knockout experiments because of the need for specific sequence information at the 5' end of the gene, the requirement for a means of delivery of the oligonucleotides to the target cells, and the lack of a systemic response.</p> <p>We have previously shown that <italic>Panagrolaimus </italic>is a suitable model for studies on the molecular basis of anhydrobiosis [##REF##15939782##57##]. Panagrolaimid nematodes occupy a diversity of niches ranging from polar, temperate and semi-arid soils to terrestrial mosses. Thus many <italic>Panagrolaimus </italic>species are adapted to tolerate freezing and desiccation [##UREF##4##81##, ####UREF##5##82##, ##REF##12477892##83####12477892##83##]. Our current data show that when fed on <italic>E. coli </italic>expressing dsRNA targeting <italic>Panagrolaimus </italic>genes the nematode mounts a potent and specific systemic RNAi response. This makes <italic>Panagrolaimus </italic>an ideal system with which to carry out high-throughput RNAi screens for studies on the molecular basis of anhydrobiosis and cryobiosis. <italic>Panagrolaimus </italic>is a member of the suborder Tylenchina of the Order Rhabditida and thus it is more distantly related to <italic>C. elegans </italic>than <italic>Oeschius </italic>or <italic>Pristionchus </italic>(Figure ##FIG##5##6##). The Suborder Tylenchina contains three major lineages [##UREF##6##84##]: the Panagrolaimorpha which includes the free living bacterial feeding panagrolaimids as well as the zooparasitic strongyloids and steinernematids; the Cephalobomorpha, bacterial feeding nematodes that, like the panagrolaimaids, have a widespread distribution and are capable of surviving in extremely dry and cold environments and the Tylenchomorpha, a group of plant parasites and fungal feeders which frequently has anhydrobiotic encysted infective stages (e.g. <italic>Globodera</italic>, <italic>Meliodogyne</italic>), as well as species which can enter anhydrobiosis throughout most or all of their life cycle (e.g. <italic>Aphelencus avenae</italic>). Within the genus <italic>Panagrolaimus </italic>several changes in reproductive mode have occurred. Examples of gonochoristic (male female), hermaphroditic (self-fertilization), and parthenogenetic (development of unfertilized eggs) reproduction have been reported [##REF##10772792##85##]. Thus <italic>Panagrolaimus </italic>is an excellent model system to study the evolutionary origins and developmental consequences of parthenogenetic reproduction. Being closely related to the zooparasitic strongyloids and steinernematids and the plant parasitic tylenchid nematodes, <italic>Panagrolaimus </italic>could also be used in comparative genomic approaches to study the evolution of parasitic lifestyles.</p>
<title>Conclusion</title> <p>Our results show that potent and specific RNAi effects can be achieved in two species of <italic>Panagrolaimus </italic>by feeding them on <italic>E. coli </italic>expressing dsRNA for genes with essential conserved functions. RNAi was first described in the nematode <italic>C. elegans</italic>, however some nematode species are refractory to RNAi, including <italic>O. tipulae </italic>[##REF##12586702##37##,##UREF##0##38##], <italic>P. pacificus </italic>[##REF##15155582##39##], several species of <italic>Caenorhabdtis </italic>[##REF##17563372##34##,##UREF##1##41##] and some animal parasitic nematodes such as <italic>H. contortus </italic>[##REF##16469321##45##] and <italic>O. ostertagi </italic>[##REF##16879764##46##]. In contrast many species of plant parasitic nematodes have potent RNAi responses to exogenous RNAi [##REF##12182331##48##,##REF##16255249##49##]. Here we show that two species of <italic>Panagrolaimus </italic>are amenable to RNAi by feeding. <italic>Panagrolaimus </italic>is a soil dwelling bacterial feeding nematode, as are <italic>C. elegans</italic>, <italic>P. pacificus </italic>and <italic>O. tipulae</italic>, but as a member of the suborder Tylenchina <italic>Pangarolaimus </italic>is more distantly related to <italic>C. elegans </italic>than <italic>Oeschius </italic>or <italic>Pristionchus</italic>. Thus our data identify another nematode infraorder which is amenable to RNAi and provide additional information on the diversity of RNAi phenotypes in nematodes.</p> <p>Panagrolaimaid nematodes occupy a diversity of niches ranging from polar, temperate and semi-arid soils to terrestrial mosses. Thus many <italic>Panagrolaimus </italic>species are adapted to tolerate freezing and desiccation [##REF##15939782##57##,##UREF##4##81##, ####UREF##5##82##, ##REF##12477892##83####12477892##83##]. Within the genus <italic>Panagrolaimus </italic>several changes in reproductive mode have occurred. Thus <italic>Panagrolaimus </italic>is an excellent model system to study the evolutionary origins and developmental consequences of parthenogenetic reproduction. Our demonstration that <italic>Panagrolaimus </italic>is amenable to RNAi by feeding will allow the development of high throughput methods of RNAi screening for <italic>P. superbus</italic>. This greatly enhances the utility of this nematode as a model system for the study of the molecular biology of anhydrobiosis and cryobiosis and as a possible satellite model nematode for comparative and functional genomics.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Gene silencing by RNA interference (RNAi) is a powerful tool for functional genomics. Although RNAi was first described in <italic>Caenorhabditis elegans</italic>, several nematode species are unable to mount an RNAi response when exposed to exogenous double stranded RNA (dsRNA). These include the satellite model organisms <italic>Pristionchus pacificus </italic>and <italic>Oscheius tipulae</italic>. Available data also suggest that the RNAi pathway targeting exogenous dsRNA may not be fully functional in some animal parasitic nematodes. The genus <italic>Panagrolaimus </italic>contains bacterial feeding nematodes which occupy a diversity of niches ranging from polar, temperate and semi-arid soils to terrestrial mosses. Thus many <italic>Panagrolaimus </italic>species are adapted to tolerate freezing and desiccation and are excellent systems to study the molecular basis of environmental stress tolerance. We investigated whether <italic>Panagrolaimus </italic>is susceptible to RNAi to determine whether this nematode could be used in large scale RNAi studies in functional genomics.</p> <title>Results</title> <p>We studied two species: <italic>Panagrolaimus </italic>sp. PS1159 and <italic>Panagrolaimus superbus</italic>. Both nematode species displayed embryonic lethal RNAi phenotypes following ingestion of <italic>Escherichia coli </italic>expressing dsRNA for the <italic>C. elegans </italic>embryonic lethal genes <italic>Ce-lmn-1 </italic>and <italic>Ce-ran-4</italic>. Embryonic lethal RNAi phenotypes were also obtained in both species upon ingestion of dsRNA for the <italic>Panagrolaimus </italic>genes <italic>ef1b </italic>and <italic>rps-2</italic>. Single nematode RT-PCR showed that a significant reduction in mRNA transcript levels occurred for the target <italic>ef1b </italic>and <italic>rps-2 </italic>genes in RNAi treated <italic>Panagrolaimus </italic>sp. 1159 nematodes. Visible RNAi phenotypes were also observed when <italic>P. superbus </italic>was exposed to dsRNA for structural genes encoding contractile proteins. All RNAi phenotypes were highly penetrant, particularly in <italic>P. superbus</italic>.</p> <title>Conclusion</title> <p>This demonstration that <italic>Panagrolaimus </italic>is amenable to RNAi by feeding will allow the development of high throughput methods of RNAi screening for <italic>P. superbus</italic>. This greatly enhances the utility of this nematode as a model system for the study of the molecular biology of anhydrobiosis and cryobiosis and as a possible satellite model nematode for comparative and functional genomics. Our data also identify another nematode infraorder which is amenable to RNAi and provide additional information on the diversity of RNAi phenotypes in nematodes.</p>
<title>Authors' contributions</title> <p>AJS and TT did the RNAi experiments. ID selected and tested the <italic>C. elegans </italic>feeding constructs. ID and AJS carried out the microscopy. AJS did the RT-PCR experiments. JB selected and prepared the <italic>Panagrolaimus </italic>sp. PS1159 cDNA sequences. AMB did the phylogenetic reconstruction. AMB, AJS, TT and ID designed the study and AMB coordinated the study. AJS and AMB drafted the manuscript. All authors read and approved the final manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was funded by a Science Foundation Ireland Investigator Programme Grant and by the Leverhulme Trust. We thank Dr. Julie Ahringer who kindly sent us the <italic>Caenorhabditis elegans </italic>feeding strains for the <italic>Ce-ran-4 </italic>and <italic>Ce-lmn-1 </italic>genes.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>The RNAi phenotypes obtained among the progeny of <italic>Panagrolaimus </italic>females fed on <italic>E. coli </italic>clones expressing dsRNA for the <italic>Caenorhabaditis elegans </italic>for the embryonic lethal genes <italic>Ce-lmn-1 </italic>and <italic>Ce-ran-4</italic></bold>. (A) Embryo from a <italic>P. superbus </italic>female fed on <italic>Ce-lmn-1</italic>; (B) Embryo from a <italic>P. superbus </italic>female fed on <italic>Ce-ran-4</italic>; (C) Control embryo from <italic>P. superbus </italic>fed on <italic>E. coli </italic>containing the empty L4440 plasmid vector; (D) Embryo from a <italic>Panagrolaimus </italic>sp. PS1159 female fed on <italic>Ce-lmn-1</italic>; (E) Embryo from a <italic>Panagrolaimus </italic>sp. PS1159 female fed on <italic>Ce-ran-4</italic>; (F) Control embryo from <italic>Panagrolaimus </italic>sp. PS1159 fed on <italic>E. coli </italic>containing the empty L4440 plasmid vector; (G) Embryo from a <italic>C. elegans </italic>female fed on <italic>Ce-lmn-1</italic>; (H) Embryo from a <italic> C. elegans </italic>female fed on <italic>Ce-ran-4</italic>; (I) Control embryo from <italic>C. elegans </italic>fed on <italic>E. coli </italic>containing the empty L4440 plasmid vector. The scale bar represents 30 μM.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>The penetrance of the dsRNAi embryonic lethal phenotypes observed among the progeny of <italic>C. elegans</italic>, <italic>P. superbus </italic>and <italic>Panagrolaimus </italic>sp. PS1159 following ingestion of <italic>E. coli </italic>expressing dsRNA for the embryonic lethal genes <italic>Ce-lmn-1 </italic>and <italic>Ce-ran-4</italic></bold>. (A) The percentage of embryonic lethal RNAi phenotypes observed in progeny of nematodes exposed to <italic>Ce-ran-4</italic>; (B) The percentage of embryonic lethal RNAi phenotypes observed in progeny of nematodes exposed to <italic>Ce-lmn-1</italic>. For both genes tested the control nematodes were fed on <italic>E. coli </italic>containing the empty L4440 plasmid vector. Data are presented as the mean and standard error of the mean of a minimum of 5 biological replicates, each replicate count containing a minimum of 50 embryos. ***<italic>P </italic>&lt; 0.001; **<italic>P </italic>&lt; 0.01 as compared to controls using a 2-tailed T test.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>The percentage of embryonic lethal RNAi phenotypes obtained by feeding <italic>P. superbus </italic>and <italic>Panagrolaimus </italic>sp. PS1159 on <italic>E. coli </italic>feeding clones expressing dsRNA for the <italic>ef1b </italic>and <italic>rps-2 </italic>genes from <italic>Panagrolaimus </italic>sp. PS1159</bold>. (A) The percentage of embryonic lethal progeny observed among the progeny of female nematodes exposed to dsRNA for <italic>P-rps-2</italic>; (B) The percentage of embryonic lethal progeny observed in nematodes exposed to dsRNA for <italic>P-ef1b</italic>. For both genes tested the control nematodes were fed on <italic>E. coli </italic>containing the empty L4440 plasmid vector. Data are presented as the mean and standard error of the mean of a minimum of 5 biological replicates, each replicate containing a minimum of 20 embryos ***<italic>P </italic>&lt; 0.001; **<italic>P </italic>&lt; 0.01 as compared to controls using a 2-tailed T test.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Single nematode RT-PCR to assess transcript levels for the <italic>P-ef1b </italic>and <italic>P-rps-2 </italic>genes in control and RNAi treated <italic>Panagrolaimus </italic>sp. PS1159 nematodes.</bold> RNA was isolated from individual nematodes whose progeny showed a positive RNAi phenotype and from individual control nematodes fed on <italic>E. coli </italic>containing the empty L4440 plasmid vector. Transcript levels were determined by RT-PCR for the target gene and for the <italic>rRNA D3 </italic>gene (which served as an internal control). (A) Transcript levels for the <italic>P-ef1b </italic>and <italic>rRNA D3 </italic>genes in individual <italic>Panagrolaimus </italic>sp. PS1159 adult females fed on <italic>E. coli </italic>expressing <italic>P-ef1b </italic>dsRNA (each of the three replicates corresponds to a separate nematode); (B) Transcript levels for the <italic>P-rps-2 </italic>and <italic>rRNA D3 </italic>genes in individual <italic>Panagrolaimus </italic>sp. PS1159 females fed on <italic>E. coli </italic>expressing <italic>P-rps-2 </italic>dsRNA (each of the three replicates corresponds to a separate nematode); (C, D) Quantification of the RT-PCR products for the <italic>P-ef1b </italic>and <italic>P-rps-2 </italic>genes in control and RNAi treated <italic>Panagrolaimus </italic>sp. PS1159 nematodes. The band intensities are given in ImageQuant™ machine units (imu).</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Larval arrest phenotypes obtained by plating <italic>P. superbus </italic>eggs onto agar plates containing a lawn of <italic>E. coli </italic>expressing dsRNA for <italic>P. superbus </italic>genes encoding contractile proteins</bold>. (A) dsRNAi phenotype for <italic>Psu-ifb-1 </italic>(cytoplasmic intermediate filament protein) and (B) dsRNAi phenotype for <italic>Psu-act-2 </italic>(actin-2). These images were taken 9 days post hatching. By this time the control embryos which had been plated onto agar plates containing a lawn of <italic>E. coli </italic>containing the empty L4440 plasmid vector had reached adulthood (C). The scale bar represents 100 μM.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>The molecular phylogeny of selected members of the Order Rhabditida showing the phylogenetic relationships between <italic>Caenorhabditis elegans </italic>and other rhabditid and tylenchid nematodes</bold>. Lineages in blue correspond to free-living nematodes, plant parasite lineages are green and animal parasite lineages are red. This Minimum Evolution tree [92] was constructed using 18S small subunit rRNA gene sequences. Bootstrap values (10,000 replicates) are placed next to the nodes.</p></caption></fig>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Table 1</bold>. Sequences of the PCR primers used for single nematode RT-PCR.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p><bold>Table 2</bold>. Accession Numbers of the small subunit ribosomal RNA sequences used to construct a molecular phylogeny of selected members of the Order Rhabditida.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p><bold>Additional figure</bold>. The molecular phylogeny of selected members of the Order Rhabditida showing the phylogenetic relationships between Caenorhabditis elegans and other rhabditid and tylenchid nematodes.  Lineages in blue correspond to free-living nematodes, plant parasite lineages are green and animal parasite lineages are red.  This Neighbor-Joining tree [91] was constructed using 18S small subunit rRNA gene sequences.  Bootstrap values (10,000 replicates) are placed next to the nodes.</p></caption></supplementary-material>
<graphic xlink:href="1471-2199-9-58-1"/> <graphic xlink:href="1471-2199-9-58-2"/> <graphic xlink:href="1471-2199-9-58-3"/> <graphic xlink:href="1471-2199-9-58-4"/> <graphic xlink:href="1471-2199-9-58-5"/> <graphic xlink:href="1471-2199-9-58-6"/>
<media xlink:href="1471-2199-9-58-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2199-9-58-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2199-9-58-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"collab": ["F\u00e9lix MA"], "surname": ["The C elegans Research Community"], "italic": ["Oscheius tipulae"], "source": ["WormBook"], "year": ["2006"]}, {"surname": ["Descotte", "Montgomery"], "given-names": ["V", "MK"], "article-title": [" Differing susceptibilities to RNAi within the "], "italic": ["Caenorhabditis clade"], "source": ["Dev Biol"], "year": ["2003"], "volume": ["259"], "fpage": ["466"]}, {"surname": ["Meldal", "Debenham", "De Ley", "De Ley", "Vanfleteren", "Vierstraete", "Bert", "Borgonie", "Moens", "Tyler", "Austen", "Blaxter", "Rogers", "Lambshead"], "given-names": ["BH", "NJ", "P", "IT", "JR", "AR", "W", "G", "T", "PA", "MC", "ML", "AD", "PJ"], "article-title": ["An improved molecular phylogeny of the Nematoda with special emphasis on marine taxa"], "source": ["Mol Phylo Evol"], "year": ["2006"], "volume": ["42"], "fpage": ["622"], "lpage": ["636"]}, {"collab": ["F\u00e9lix MA", "Delattre M", "Dichtel ML"], "article-title": ["Comparative developmental studies using Oscheius/Dolichorhabditis sp. CEW1 (Rhabditidae)."], "source": ["Nematology "], "year": ["2000"], "volume": ["2"], "fpage": ["88"], "lpage": ["89"]}, {"surname": ["Ricci", "Pagani"], "given-names": ["C", "M"], "article-title": ["Desiccation of "], "italic": ["Panagrolaimus rigidus"], "source": ["Hydrobiologia"], "year": ["1997"], "volume": ["347"], "fpage": ["1"], "lpage": ["13"], "pub-id": ["10.1023/A:1002979522816"]}, {"surname": ["Aroian", "Carta", "Kaloshian", "Sternberg"], "given-names": ["RV", "L", "I", "PW"], "article-title": ["A free-living "], "italic": ["Panagrolaimus"], "source": ["J Nematol"], "year": ["1993"], "volume": ["25"], "fpage": ["500"], "lpage": ["502"]}, {"collab": ["De Ley P", "Blaxter ML"], "surname": ["Lee D"], "article-title": [" Systematic position and phylogeny"], "source": ["The Biology of Nematodes"], "year": ["2002"], "publisher-name": [" Harwood Academic Publishers"], "fpage": ["1\u201330"]}, {"article-title": ["Blast Home"]}, {"article-title": ["Wormbase"]}, {"collab": ["Lewis JA", "Fleming JT"], "surname": ["Epstein HF, Shakes DC"], "article-title": ["Basic Culture Methods"], "source": ["Caenorhabditis elegans: Modern Biological Analysis of an Organism"], "year": ["1995"], "volume": [" 48"], "publisher-name": [" Methods in Cell Biology"], "fpage": ["3"], "lpage": ["29"]}, {"article-title": ["National Center for Biotechnology Information"]}, {"surname": ["Rzhetsky", "Nei"], "given-names": ["A", "M"], "article-title": ["A simple method for estimating and testing minimum-evolution trees"], "source": ["Mol Biol Evol"], "year": ["1992"], "volume": ["9"], "fpage": ["945"], "lpage": ["967"]}]
BMC Mol Biol. 2008 Jun 19; 9:58