RNA interference of two somatic targets show stark differences in their requirements for MUT-16 and RDE-10.

(A) Overview of RNA interference in somatic cells. Double-stranded RNA (dsRNA, blue) enters the cell through the importer SID-1 (import, teal), after which it is processed by the dsRNA-binding protein RDE-4 and the endonuclease Dicer into 1° short interfering RNAs (siRNAs) that are bound by the primary Argonaute RDE-1 (1° processing, pink). mRNA transcripts (green) recognized by these 1° siRNAs are modified after cleavage by the 3′ addition of UG repeats (pUG RNA) and act as templates for the amplification of 2° siRNAs aided by the intrinsically disordered protein MUT-16, the Maelstrom-domain containing protein RDE-10, and RNA-dependent RNA polymerases (2° amplification, orange). These 2° siRNAs can bind secondary Argonaute(s) (e.g., NRDE-3), which can then recognize additional complementary targets (2° recognition) and cause gene silencing. See text for details. (B) Gene schematics depicting the mutant alleles found in a genetic screen (black) and/or created using genome editing (blue). Black boxes indicate exons and red dots indicate locations of mutations. Allele names (e.g., jam139) and expected amino acid change in the corresponding proteins (e.g., mutation of a tyrosine codon to a stop codon, Y294*) are indicated. See Fig. S1 for details of genetic screen. (C and D) Response to bli-1 or unc-22 RNAi in different mutants. For each mutant, the fraction of animals that showed bli-1 silencing or unc-22 silencing (fraction silenced) and the numbers of animals scored (n) are shown. Asterisks indicate p < 0.05 for each comparison (brackets) using Wilson’s estimates with continuity correction and error bars represent 95% confidence interval. (C) Of five isolates with a mutation in mut-16, four (jam138, jam140, jam141, and jam247) failed to silence bli-1 (blue) but retained unc-22 silencing (orange). The other mutant failed to silence both genes and additionally had a mutation in rde-10 (mut-16(jam139) rde-10(jam248)). (D) Mutants created using genome editing recapitulated the selective silencing of unc-22 in mut-16(-) single mutants (mut-16(jam148)) and the failure to silence both genes in mut-16(-) rde-10(-) double mutants (mut-16(jam148) rde-10(jam206)). Using genome editing to recreate the jam248 mutation, which is expected to make a mutant protein (RDE-10(S228F)) that disrupts the Maelstrom domain (see Fig. S1), resulted in animals (rde-10(jam196)) that showed unc-22 silencing but not bli-1 silencing. (E) Selective requirement for a regulator could reflect two underlying mechanisms of RNA silencing: (1) Two parallel pathways (left, A/B path vs C/D path) that are differentially used for different target genes; or (2) One intersecting network (right, A/B/C/D network) with quantitative contributions by all regulators along with different thresholds for each target gene.

Gene-specific requirements and complex redundancy can arise from a single RNA regulatory network.

(A) Schematic (as in Fig. 1) depicting nrde-3 alleles. (B) Feeding RNAi of bli-1 and unc-22. Fractions silenced, numbers scored, comparisons, asterisks, and error bars are as in Fig. 1. Single mutants lacking NRDE-3 (nrde-3(jam205)) fail to silence bli-1 but not unc-22. Double mutants fail to silence both targets. (C and D) Mutual constraints among parameters required for a single RNA regulatory network to support experimental results. (C, left) Model for a single network of interactors that regulate all RNAi targets in somatic cells. All targets require import (SID-1) and processing (RDE-4 and RDE-1) of dsRNA. Branching after 1° siRNA processing results in four distinct paths (Nm, Nr, Om, Or) that together contribute to gene silencing, which could occur through co-transcriptional gene silencing (TGS) and/or post-transcriptional gene silencing (PTGS) mechanisms. (C, right) Representation of simulation workflow. First, random values between 0 and 2 were drawn for each of the four variables (Nm, Nr, Om, Or). Second, constraints were added based on the experimental results in Fig. 2B and Fig. 1D. Third, allowed values that satisfied all experimental conditions were culled. Of 100,000 sets of random values simulated (0 to 2 for Nm, Nr, Om, Or and 0 to 100 for the ratio of thresholds Tbli-1/Tunc-22), 145 were consistent with all observed responses to RNAi. These allowed numbers reveal the domain of parameter values that support the observed range of gene silencing outcomes using feeding RNAi. (D, left) The contribution of NRDE-3 via MUT-16 (Nm) versus that via RDE-10 (Nr) for different ratios of thresholds for bli-1 versus unc-22 silencing (Tbli-1/Tunc-22) are shown. (D, center and right) The relative contributions to silencing that require MUT-16 (Nm + Om, D, center) or NRDE-3 (Nm + Nr, D, right) do not frequently take extreme values and both support a low value for the ratio of thresholds (Tbli-1/Tunc-22 < ∼5.5 despite allowed values of up to 100). (E) Feeding RNAi of unc-22 assayed as in Fig. 1, but using aged plates, resulting in weaker RNAi. Animals that lack MUT-16 (mut-16(jam148)) have the most severe defect, followed by animals lacking RDE-10 (rde-10(jam206)), which is followed by animals lacking NRDE-3 (nrde-3(jam205)). (F) Feeding RNAi of unc-54 or dpy-7. Fractions silenced, numbers scored, comparisons, asterisks, and error bars are as in Fig. 1. Silencing of unc-54 showed a partial dependency on NRDE-3, while silencing of dpy-7 showed a partial dependency on MUT-16, suggesting that the quantitative requirement for a regulator can differ depending on the target.

Gene-specific requirements for NRDE-3 can be bypassed in two ways.

(A) Gene schematics (as in Fig. 1) of rde-4 and eri-1. (B and C) Loss of ERI-1 can bypass the NRDE-3 requirements for silencing bli-1 but not the other requirements for silencing bli-1 or unc-22. Feeding RNAi targeting bli-1 (B) or unc-22 (C) with fractions silenced, numbers scored, comparisons, asterisks, and error bars as in Fig. 1. (B) Loss of ERI-1 (mg366, jam263, and jam264 alleles) can compensate for the role of NRDE-3 (tm1116 and jam205 alleles) but not of RDE-10 (jam206 allele) or MUT-16 (jam148, jam265, jam266, jam267, and jam268 alleles) in bli-1 silencing. See Table S3 for additional information. (C) Silencing of unc-22 is not restored by loss of ERI-1 (mg366, jam260, jam261, and jam262 alleles) in mutants that also lack any two of mut-16 (jam148 and jam240 alleles), rde-10 (jam196 and jam206 alleles), or nrde-3 (jam205 allele). See Table S3 for additional information. (D to E) Overexpression of RDE-4 in the hypodermis can bypass the requirement for NRDE-3 in bli-1 silencing. (D) Minimal amounts of RDE-4 are sufficient for somatic silencing. (top) Schematic depicting generation of male progeny with paternal inheritance of a single-copy transgene (Si[…]) that expresses rde-4(+) or rde-1(+) under the control of the mex-5 promoter (mex-5p) in the germline (green) of rde-4(-) or rde-1(-) animals, respectively (germline– and intestine-enriched RDE, based on rescue of RNAi in rde-1(-) animals (39)). (bottom) Male cross progeny with the transgene were scored after feeding only F1 animals, showing that unlike animals with germline– and intestine-enriched RDE-1, animals with similarly enriched RDE-4 can rescue both unc-22 and bli-1 silencing. Thus, small amounts of RDE-4 potentially mis-expressed in the hypodermis or a non-autonomous effect of RDE-4 from the germline or intestine is sufficient for silencing in the muscle and hypodermis. † indicates p<0.05 when compared to either wild type or the rde-4(-) mutant and other symbols are as in (B). (E) Silencing of bli-1 is restored in nrde-3(tm1116); rde-4(ne301) double mutants when rde-4(+) is overexpressed in the hypodermis (Si[nas-9p::rde-4(+)]). (F) Silencing of bli-1 is not restored in animals lacking MUT-16 or RDE-10, despite the overexpression of RDE-4 in the hypodermis and/or the loss of ERI-1. (G) Summary depicting the bypass of NRDE-3 when ERI-1 is eliminated and/or RDE-4 is overexpressed. The increase in dsRNA processing increases the contributions of NRDE-3-independent paths to silencing.

A quantitative model allows exploration of parameters for RNAi in C. elegans.

(A) RNAi can be modeled to understand (from left) the responding populations of animals, the outline of mechanism, the process highlighting RNA species, the biochemical sub-steps (e.g., 2° siRNA production using pUG RNA, right). This study focuses on the overall RNA-mediated processes that accompany RNAi. (B) Schematic (left) and ordinary differential equations (right) describing the production and turnover of different RNA species. Rate constant for 1° siRNA processing from dsRNA (k1), binding constant for 1° siRNA binding target mRNA (k2), rate constant for pUG RNA production (k3), rate constant for 2° siRNA production (k4), binding constant for 2° siRNAs binding mRNA (k5) or pre-mRNA (k6), rate constant for export and splicing of transcripts from the nucleus (k7), rate constant for repression of transcription (k8), and rate constant for new transcript production (k9). Other terms are described in the figure. (C) Relative changes in the concentrations of each RNA ([RNA] vs. time in a.u.; dsRNA, mRNA, pre-mRNA, pUG RNA, 1° siRNA, and 2° siRNA) for an example set of parameters (all turnover rates = 0.05, k1 = 1, k2 = 0.01, k3 = 1, k4 = 0.05*lm = 0.5, k5 = 0.01, k6 = 0.01, k7 = 0.1, k8 = 0.05, k9 = 7.5) are illustrated. A reduction to 10% of initial mRNA concentration is designated as the threshold for detecting a defect upon knockdown ([m]kd), the time needed to reach the threshold (tkd) and the time for which mRNA levels remain below the threshold (tkd) are also indicated. (D) Relationship between the duration of knockdown and the time to knockdown (tkd and tkd are as in (C)). (E) Relationship between mRNA concentration and 2° siRNA accumulation. The minimum mRNA concentrations and maximum 2° siRNA concentrations reached for different transcripts with two different binding constants of 2° siRNAs binding to mRNA (k5 = 0.1, red and k5 = 0.01, blue) are plotted. Also see Fig. S5. (F) Impact of doubling transcription on transcripts with different knockdown parameters. Each transcript is colored based on its initial duration of knockdown (tkd, blue to red gradient) before a 2-fold increase in the rate constant for transcription (k9) and the resultant fractional change in the duration of knockdown tkd is plotted against that in the time to knockdown tkd. (G) Genes with higher turnover are harder to knockdown. Response of mRNAs and their respective 1° siRNA with the same steady-state concentrations but with different rates of mRNA turnover (solid lines: Tm = 0, large dashes: Tm = 0.05, small dashes: Tm = 0.5) upon addition of 10 molecules of dsRNA are shown. (inset) Relationship of the minimum concentration of mRNA ([mRNA]min) with its Tm in response to a fixed amount of dsRNA.

Changing cis-regulatory regions or open reading frames can alter the observed dependence on NRDE-3 for silencing by dsRNA.

(A) Schematics (as in Fig. 1) of the hybrid genes created to test the role of 3′ and 5′ cis-regulatory sequences in the response to dsRNA (see Methods for details). (B) The need for NRDE-3 to cause silencing using bli-1-dsRNA is reduced when the 3′ cis-regulatory sequences of bli-1 are replaced with that of unc-22 (bli-1(jam307)). Fractions silenced, numbers scored, comparisons, asterisks, and error bars are as in Fig. 1 (C) Representative images (DIC, top; GFP fluorescence, bottom) of the head region of bli-1p::gfp::cye-1 3′cis animals either in response to a control vector (pL4440, left) or in response to gfp-dsRNA (right)). The presence or absence of gfp in the head region was used to score the fraction silenced in (D). Scale bar = 100µm. (D) Use of the bli-1 promoter alone does not confer the genetic requirements for silencing bli-1. Silencing of GFP expressed from the bli-1p::gfp::cye-1 3′cis transgene (bli-1p::gfp(jam297)) by gfp-dsRNA was measured in a wild-type background, in animals lacking NRDE-3 (nrde-3(jam205)), in animals lacking RDE-10 (rde-10(jam206)), in animals lacking MUT-16 (mut-16(jam148)), and in animals lacking both MUT-16 and RDE-10 (mut-16(jam148) rde-10(jam196)). Fractions silenced, numbers scored, comparisons, asterisks, and error bars are as in Fig. 1.

Animals recover from a pulse of RNAi and production of pUG RNAs is restricted despite continuous exposure to dsRNA.

(A) Response to a pulse of feeding RNAi. (top) Schematic of assay. Animals were exposed to unc-22 RNAi for one hour and then returned to OP50 plates. (bottom) A separate cohort of animals was scored for silencing after each subsequent 24hr period. Fractions silenced, numbers scored, comparisons, asterisks, and error bars are as in Fig. 1. A weak Unc-22 defect indicates animals that were nearly completely still except for a slight twitch of the head or of the tail. (B) pUG RNA production in response to continuous exposure to unc-22 dsRNA. (top) Schematic depicting the PCR primers used for the RT-PCR to detect pUG RNAs. Two sets of primers (0kb, purple; 1kb, blue) positioned 5′ of the unc-22-dsRNA (orange), combined with an RT primer that contains nine CA repeats were used. (bottom) Populations of small RNAs that might be detected by pUG-PCR as above. Population I would suggest that amplified 2° siRNAs in addition to 1° siRNAs are capable of guiding target mRNA cleavage and poly-UG addition. Population II would suggest only 1° siRNAs can initiate pUG RNA production. (C) (left) Distribution of DNA amplified from pUG RNAs. Lanes with PCR products amplified from total RNA of animals fed unc-22 dsRNA (unc-22) or L4440 dsRNA (control) isolated from three biological replicates each (1-3), or a no-template control (NTC) with no added RNA (-) are shown. Different bands are detected for each primer set (0kb, top vs. 1kb, bottom). A gene with a poly-UG sequence encoded in the genome (gsa-1) serves as a loading control. (right) Schematic of siRNA production in somatic cells. Successive rounds of amplified small RNAs would map progressively closer to the 5′ end of the target transcript. Since we did not detect poly-UG addition upstream of the region homologous to the dsRNA trigger, it is unlikely that 3° siRNAs are being produced.

Overview of findings.

(left) Amplification of small RNAs occurs through an intersecting network of regulators (orange arrows) and in response to somatic RNAi, the addition of poly-UG repeats (green) is restricted to a ‘pUG zone’ (red) that is homologous to the dsRNA trigger. (right) In response to a pulse of dsRNA (red), levels of the target mRNA (blue) recover, suggesting that the RNA silencing intermediates (1° siRNAs, purple; pUG RNAs, green; 2° siRNAs, brown) undergo turnover.