A DCL3 dicing code within Pol IV-RDR2 transcripts diversifies the siRNA pool guiding RNA-directed DNA methylation
Abstract
In plants, selfish genetic elements, including retrotransposons and DNA viruses, are transcriptionally silenced by RNA-directed DNA methylation. Guiding the process are short interfering RNAs (siRNAs) cut by DICER-LIKE 3 (DCL3) from double-stranded precursors of ~30 bp that are synthesized by NUCLEAR RNA POLYMERASE IV (Pol IV) and RNA-DEPENDENT RNA POLYMERASE 2 (RDR2). We show that Pol IV’s choice of initiating nucleotide, RDR2’s initiation 1–2 nt internal to Pol IV transcript ends and RDR2’s terminal transferase activity collectively yield a code that influences which precursor end is diced and whether 24 or 23 nt siRNAs are produced. By diversifying the size, sequence, and strand specificity of siRNAs derived from a given precursor, alternative patterns of DCL3 dicing allow for maximal siRNA coverage at methylated target loci.
Editor's evaluation
The paper is of interest to RNA biologists, especially to those who study small RNAs. The findings deepen our understanding of the rules of DCL3 dicing and explain how 23-nt and 24-nt siRNAs in the RdDM pathway are produced.
https://doi.org/10.7554/eLife.73260.sa0Introduction
In eukaryotes, short interfering RNAs (siRNAs) are used to suppress DNA transcription or mRNA translation, thereby playing important roles in gene regulation (Shabalina and Koonin, 2008; Ipsaro and Joshua-Tor, 2015; Martienssen and Moazed, 2015). Plant siRNAs range in size from 21 to 24 nt, with DCL3-dependent 23 and 24 nt siRNAs (Xie et al., 2004; Henderson et al., 2006) accounting for ~90% of the total siRNA pool (Kasschau et al., 2007; Zhang et al., 2007; Mosher et al., 2008). The 24 nt siRNAs stably associate with one of several Argonaute family proteins, primarily ARGONAUTE 4 (AGO4) (Zilberman et al., 2003; Qi et al., 2006), and guide resulting complexes to target loci via base-pairing interactions with long noncoding RNAs synthesized by multisubunit NUCLEAR RNA POLYMERASE V (Figure 1A; Wierzbicki et al., 2008; Wierzbicki et al., 2009). Protein-protein interactions between AGO4 and the C-terminal domain of the Pol V largest subunit, or Pol V-associated protein SPT5L, also contribute to AGO4 localization at target loci (El-Shami et al., 2007; He et al., 2009; Lahmy et al., 2016). Subsequent recruitment of the de novo DNA methyltransferase, DRM2 (Cao and Jacobsen, 2002; Zhong et al., 2014) and other chromatin modifying enzymes then leads to the establishment of repressive chromatin environments that inhibit promoter-dependent transcription by DNA-dependent RNA Polymerases I, II, or III (Matzke and Mosher, 2014; Wendte and Pikaard, 2017). In this way, RNA-directed DNA methylation (RdDM) facilitates transcriptional silencing at thousands of loci throughout plant genomes.
-
Figure 1—source data 1
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig1-data1-v2.pdf
-
Figure 1—source data 2
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig1-data2-v2.pdf
-
Figure 1—source data 3
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig1-data3-v2.pdf
Biogenesis of siRNAs involved in RdDM begins with DNA transcription by Pol IV, a 12-subunit DNA-dependent RNA polymerase that evolved as a specialized form of Pol II (Ream et al., 2009). Pol IV associates with RDR2 (Law et al., 2011; Haag et al., 2012) via direct physical interaction (Mishra et al., 2021) to form a multi-functional enzyme complex. The transcription reactions of Pol IV and RDR2 are tightly coupled (Singh et al., 2019) with Pol IV’s inability to maintain processivity through double-stranded DNA regions proposed to results in Pol IV arrest and backtracking, coincident with DNA template-nontemplate strand reannealing and extrusion of the nascent transcript’s 3′ end for engagement by RDR2 (Fukudome et al., 2021). Resulting double-stranded RNAs (dsRNAs), ranging in size from ~25 to 40 bp in vivo (Blevins et al., 2015; Zhai et al., 2015; Yang et al., 2016; Ye et al., 2016) are then diced by DCL3. Biogenesis of 23 and 24 nt siRNAs can be recapitulated in vitro upon incubating purified Pol IV, RDR2, and DCL3 with single-stranded bacteriophage M13 template DNA and nucleoside triphosphates, indicating that no other activities are needed (Singh et al., 2019).
Although DCL3 generates both 23 and 24 nt siRNAs, immunoprecipitated AGO4 associates almost entirely with 24 nt siRNAs (Mi et al., 2008; Havecker et al., 2010), making the significance of 23 nt siRNAs unclear. We have proposed that when paired with 24 nt siRNAs, the 23 nt RNAs help specify the stable AGO4-association of the 24 nt strands (Singh et al., 2019). In this study, we investigated the rules of DCL3 dicing, revealing how Pol IV and RDR2-encoded cues program DCL3 dicing patterns and account for the biogenesis of either 23 or 24 nt siRNAs. We show that DCL3 preferentially binds dsRNAs with 3′ overhangs and that overhangs are present at both ends of Pol IV-RDR2 dsRNA transcripts to facilitate dicing from either end, with two distinct activities of RDR2 accounting for these 3′ overhangs. The choice of nucleotide incorporated at the 5′ terminus of Pol IV transcripts also affects DCL3 activity, and which end of the dsRNA is diced. Our evidence indicates that DCL3 measures only one strand of its substrate dsRNAs, with RNase III domain B cutting the measured strand to produce 24 nt siRNAs and domain A cutting the non-measured strand to produce either 23 or 24 nt siRNAs. Collectively, our experiments reveal a code comprised of sequence and structural features intrinsic to paired Pol IV-RDR2 transcripts that programs alternative DCL3 dicing patterns, thereby diversifying the siRNA pool that guides RNA-directed DNA methylation to complementary target loci.
Results
DCL3 preferentially dices double-stranded RNAs with 3′ overhangs
An untemplated nucleotide is often present at the 3′ ends of RDR2 transcripts and is enriched among 23 nt, but not 24 nt siRNAs (Singh and Pikaard, 2019). RDR2’s terminal transferase activity can account for the addition of these untemplated 3′ terminal nucleotides (Blevins et al., 2015). These observations led to the model in Figure 1B in which DCL3 interacts with base-paired Pol IV-RDR2-transcripts, measures and cuts 24 nt from the 5′ terminus of the Pol IV strand (shown as the top strand in Figure 1B and throughout this paper), and makes a second cut of the RDR2 strand that is offset by 2 nt from the Pol IV-strand cut site (Singh et al., 2019; Singh and Pikaard, 2019). Because of the untemplated nucleotide added to the 3′ end of the RDR2 strand, the predicted dicing product is an asymmetric 24/23 nt siRNA duplex (Singh et al., 2019).
To test whether DCL3 can, in fact, carry out the hypothetical reaction of Figure 1B, recombinant DCL3 produced in insect cells was purified to near-homogeneity (Figure 1—figure supplement 1 -panel A) and incubated with a dsRNA generated by annealing 37 and 38 nt RNA strands (see Supplementary file 1a for RNA strand sequences). This dsRNA has a 3′ overhang of 1 nt on the left side, as drawn in Figure 1C, to mimic the overhang attributable to nucleotide addition by RDR2’s terminal transferase activity. DCL3 cuts the 37/38 nt dsRNA to yield 23 and 24 nt products (Figure 1C, lanes 1 and 5) as predicted by the model. Similarly, DCL3 cuts a dsRNA with a 2 nt overhang, formed by annealing 37 and 39 nt RNAs, generating two 24 nt RNAs (Figure 1C, lanes 2 and 6). These results are consistent with the hypothesis that DCL3 measures and cuts 24 nt from the 5′ end of the top strand, and makes an offset cut of the bottom strand to leave a 2 nt 3′ overhang, as expected from prior studies (Zhang et al., 2004; Takeshita et al., 2007; Nagano et al., 2014). The results also suggest that the length of the bottom strand can vary, indicating that its 3′ end is not anchored at a fixed position within DCL3.
Given a dsRNA substrate with two blunt ends, DCL3 generates dicing products of diverse sizes, including 21, 22, 24, and 25 nt RNAs (Figure 1C, compare lanes 3 and 7). Gel mobility depends on RNA sequence, thus RNAs of the same length, but from opposite strands of the dsRNA, differ in mobility, as is apparent in Figure 1C. Importantly, 21, 22, and 25 nt dicing products are not observed at significant levels among siRNAs produced in vitro by Pol IV, RDR2, and DCL3 (Singh et al., 2019), nor have RNAs of these sizes been attributed to DCL3 activity based on analysis of dcl3 mutants in vivo (Xie et al., 2004; Henderson et al., 2006). These considerations, combined with the results of Figure 1C, indicate that dsRNAs with 3′ overhangs, not blunt ends, best account for the 23 and 24 nt siRNAs generated by DCL3.
To more definitively test whether DCL3 prefers 3′ overhangs over blunt ends, we performed DCL3 cleavage assays using dsRNAs that have a 3′ overhang on the left side and a blunt end on the right side, with either the top strand or bottom strand end-labeled with 32P (Figure 1D). If the 3′ overhanging end is preferred by DCL3, a 24 nt 32P-labeled dicing product is expected from the 32P-labeled top strand and a 16 nt product is expected from the 32P-labeled bottom strand. This is observed, with dsRNA substrates that have 3′ overhangs of 1 or 2 nt yielding nearly identical products (Figure 1D, compare lanes 1–4 to lanes 5–8) and reaction time-courses (Figure 1—figure supplement 1 – (B)). In contrast, if DCL3 is incubated with a dsRNA with two blunt ends, DCL3 dices from both ends, generating 16 and 24 nt products from both strands (Figure 1E). Collectively, the results of Figure 1 show that DCL3 will cut dsRNAs with 3′ overhanging ends or blunt ends but prefers 3′ overhanging ends.
DCL3 measures 24 nt from the recessed 5′ end of dsRNAs with 3′ overhangs
Experiments testing the dicing of dsRNAs formed by annealing RNA strands ranging in size from 22 to 25 nt yielded additional evidence that DCL3 measures only one strand of its dsRNA substrates (Figure 2). Using a duplex whose top strand was 24 nt and whose bottom strand was 25 nt, thus generating a 1 nt 3′ overhang on the left side, the 24 nt top strand was not cut by DCL3, but the bottom strand was trimmed by 2 nt to yield a 23 nt product (Figure 2A). We next tested a duplex with two blunt ends, formed by annealing two 24 nt RNAs (Figure 2B; note that these 24 nt RNAs have different gel mobilities). Dicing by DCL3 occurred from the left or right sides, due to the absence of a 3′ overhang to bias dicing to one side. In each case, one 24 nt strand was uncut whereas the paired strand was trimmed to 22 nt RNA (Figure 2B, lane 3; note that the alternative 22 nt RNAs also have different gel mobilities).
An informative dicing pattern was observed for a duplex formed by annealing a 23 nt top strand with a 24 nt bottom strand and having a 1 nt 3′ overhang on the left side. Dicing of this substrate resulted in the 24 nt bottom strand being cut, but trimmed by only 1 nt (Figure 2B, lanes 4 and 5). This indicates that DCL3’s measurement of an interval equivalent to 24 nt is a consequence of interaction with the recessed 5′ end of the top strand. The bottom strand was cut at a position offset by 2 nt from where the 24th nucleotide would have been present in the top strand had the top strand been long enough.
Next, we annealed a 24 nt top strand and a 23 nt bottom strand to generate a duplex with a 1 nt 3′ overhang on the left side and a 2 nt 3′ overhang on the right side (Figure 2C), which mimics the dicing product proposed in Figure 1B (and Singh et al., 2019) and verified in Figure 1C. DCL3 does not cut either stand of this duplex RNA (Figure 2C, compare lanes 2 and 3). This suggests that the asymmetric 24/23 nt dsRNA duplex fits into the enzyme with the RNase domains aligned with the pre-existing ends of the RNA strands. As a result, no cutting occurs. A 24/24 nt duplex having 2 nt 3′ overhangs at both ends was also not cut by DCL3 (Figure 2C, compare lanes 4 and 5). These results indicate that DCL3 dicing products are not substrates for further dicing.
Unlike 24/23 and 24/24 dicing products that are not diced further, duplexes formed by annealing 23 and 22 nt RNA strands, in various permutations, yielded a ladder of digestion products as short as 17 nt (Figure 2D). This suggests that dsRNAs whose strands are shorter than typical DCL3 dicing products do not fit into the enzyme in a way that fixes their position, allowing for cutting at variable positions. The results also raise the possibility that DCL3 could play a role in the cleavage and turnover of dsRNAs smaller than its dicing products in vivo.
Strand-cutting specificities of DCL3’s RNase III domains
To determine how DCL3 is oriented on its dsRNA substrates, we mutated RNase III domains A and B by converting glutamates E1146 and/or E1329 to glutamines (Figure 3A). We then tested the cutting of a dsRNA formed by annealing 26 and 27 nt RNAs, with a 1 nt 3′ overhang on the left side (Figure 3B). DCL3 with wild-type, unmutated RNase domains cut both strands of the dsRNA to yield 24 and 23 nt siRNAs (compare lanes 5 and 6). Mutating both RNase III domains abolished dicing (compare lanes 5 and 7). Mutating RNase III domain A, while leaving domain B unchanged, resulted in top strand (26 nt) cutting to 24 nt, but no bottom strand (27 nt) cutting (lane 8). Conversely, mutating domain B, but leaving domain A unchanged, allowed bottom strand cutting, to 23 nt (lane 9). Repeating these experiments using a dsRNA substrate with a 2 nt overhang, instead of a 1 nt overhang, yielded equivalent results, the only difference being the generation of a 24 nt product from the bottom strand due to the additional overhanging nucleotide (Figure 3—figure supplement 1). Collectively, these results indicate that the PAZ domain binds the recessed 5′ end of the top RNA strand, opposite the 3′ overhang of the bottom strand, and cuts that strand using RNase III domain B at a distance equivalent to 24 nt. The paired strand is cut by RNase III domain A, with the 3′ overhanging end of this strand not anchored at a fixed distance, allowing its diced products to vary in length (Figure 3C).
DCL3 substrate recognition is influenced by 5′ terminal nucleotide and phosphorylation status
Pol IV transcripts tend to have a purine at their 5′ ends (as denoted in Figure 1B), with A being more prevalent than G (Singh et al., 2019). However, pyrimidines can also serve as the initiating nucleotide, with C used more frequently than U (Singh et al., 2019). We thus tested whether the 5′-terminal nucleotide affects DCL3 dicing using dsRNA substrates with a 1 nt 3′ overhang on the left side, a blunt-end on the right side, and a top strand that begins with A, G, C, or U (Figure 4A). This experiment revealed that DCL3 cleaves dsRNA substrates whose top strands begin with 5′ A or U (lanes 5–10) more efficiently than those beginning with C or G (lanes 11–16), in agreement with prior results comparing cell-free lysates of wild-type and dcl3 mutant plants (Nagano et al., 2014). The fact that A or U allows for similar dicing efficiency suggests that an AU or UA pair is preferred over a GC or CG pair at the precursor’s terminus.
-
Figure 4—source data 1
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data1-v2.pdf
-
Figure 4—source data 2
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data2-v2.pdf
-
Figure 4—source data 3
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data3-v2.pdf
-
Figure 4—source data 4
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data4-v2.xlsx
-
Figure 4—source data 5
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data5-v2.pdf
-
Figure 4—source data 6
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data6-v2.pdf
-
Figure 4—source data 7
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data7-v2.pdf
-
Figure 4—source data 8
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data8-v2.xlsx
-
Figure 4—source data 9
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig4-data9-v2.pdf
We next tested whether DCL3 is affected by the nucleotide at the 3′ overhanging position of the bottom strand (Figure 4B). No significant difference in cleavage efficiency was observed for bottom strands having A, U, C, or G as the overhanging nucleotide, indicating that although the presence of 3′ overhang is important for dicing (Figure 1D), the identity of the overhanging nucleotide is unimportant.
The 5′ terminal nucleotide of a nascent transcript is expected to possess a 5′ triphosphate group, and biochemical evidence indicates that this is true for RDR2 transcripts (Singh et al., 2019). However, Pol IV-dependent transcripts generated in vivo or in vitro can be cloned via ligation reactions that require a 5′ monophosphate, without prior enzymatic treatments (Blevins et al., 2015; Li et al., 2015; Zhai et al., 2015; Singh et al., 2019), suggesting the possibility of an intrinsic or associated pyrophosphatase activity. We tested whether monophosphate or triphosphate groups at the 5′ end of a dsRNA substrate affect DCL3 dicing. For this experiment, dsRNA substrates had 1 nt 3′ overhangs (uridines in each case) and 5′ adenosines at each end (Figure 4C), making both ends similarly conducive to dicing. The top strand in all cases was 5′ end-labeled with a 32P monophosphate. The 5′ end of the bottom strand had either a triphosphate, a monophosphate, or a hydroxyl group. Dicing initiated by measuring 24 nt from the 5′ end of the top strand yields a 24 nt labeled product, whereas dicing measured 24 nt from the 5′ end of the bottom strand yields a 16 nt labeled product. The ratio of 24 versus 16 nt products thus provides a way to assess the relative affinity of DCL3 for dicing the alternative ends of the substrate dsRNA. This experiment revealed that if the bottom strand has a 5′ hydroxyl group, a strong 24 nt signal and weak 16 nt signal (barely above background levels) are observed (Figure 4C, lane 4). This indicates that DCL3 preferentially engaged the left side of the dsRNA substrate, which has a 5′ monophosphate, and disfavored the right side, which has a 5′ hydroxyl. However, if the bottom strand also has a 5′ monophosphate (lane 5), or a triphosphate (lane 6), more right-side dicing occurs suggesting that monophosphates or triphosphates are similarly conducive to DCL3 engagement, but with monophosphates slightly preferred.
DCL3 dicing does not require ATP
Drosophila melanogaster Dicer-2 displays ATP-dependent activity for substrates with blunt-ends, but not overhanging ends (Sinha et al., 2018). Moreover, experiments comparing cell-free lysates of wild-type versus dcl3 mutant Arabidopsis plants suggested that binding, but not hydrolysis, of ATP is needed for DCL3 activity (Nagano et al., 2014). We thus tested whether dicing by purified DCL3 is ATP-dependent using dsRNA substrates formed by annealing a 37 nt top strand RNA, 5′ end-labeled with 32P, to either an unlabeled 37 nt bottom strand, thereby generating a dsRNA with two blunt ends (Figure 5A), or to an unlabeled 38 nt bottom strand, thereby generating a dsRNA with a 1 nt overhang on the left side and a blunt end on the right side (Figure 5B). We also tested a dsRNAs with a 2 nt 3′ overhang (Figure 5—figure supplement 1). To preclude the potential involvement of ATP pre-associated with DCL3, we incubated FLAG-DCL3 with apyrase during affinity purification to hydrolyze any contaminating ATP to AMP and inorganic phosphate, as demonstrated in control reactions (see Figure 5—figure supplement 2). We also used a 1:5 molar ratio of DCL3:dsRNA substrate such that complete dicing of the dsRNAs would require each DCL3 protein to carry out multiple dicing reactions.
In reactions containing 0 mM ATP, 5 mM ATP, or 5 mM ATP-γ-S, a non-hydrolyzable form of ATP, no requirement of ATP for DCL3 activity was observed, as monitored by the generation of labeled 24 nt dicing products (lanes 2–10 of Figure 5A, B; see also Figure 5—figure supplement 1). We conclude that Arabidopsis DCL3 does not require ATP binding or hydrolysis for dicing.
Overhangs at each end of dsRNAs substrates explain 24 and 23 nt siRNA biogenesis from either RNA strand
The model in Figure 1B, supported by the experiments of Figures 1—3, accounts for how DCL3 can produce a diced duplex consisting of a 23 nt siRNA derived from the RDR2-transcribed strand and a 24 nt siRNA derived from the Pol IV strand. However, the model does not account for our prior RNA-seq data revealing that 23 and 24 nt siRNAs are generated from both strands (Singh et al., 2019). This led us to deduce a new model (Figure 6A) in which DCL3’s ability to dice from either the left or right side of a duplex and DCL3’s penchant for 3′ overhangs were considered. Scenario 1 of Figure 6A is the scenario of Figure 1B, with the left-side 3′ overhang generated by RDR2’s addition of an untemplated nucleotide (N). In scenarios 2 and 3, we hypothesized that 3′ overhangs on the right side could occur if RDR2 initiates transcription 1 or 2 nt internal to the 3′ end of the Pol IV transcript.
-
Figure 6—source data 1
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig6-data1-v2.pdf
-
Figure 6—source data 2
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig6-data2-v2.pdf
-
Figure 6—source data 3
- https://cdn.elifesciences.org/articles/73260/elife-73260-fig6-data3-v2.pdf
To experimentally test for single-stranded overhangs at the ends of Pol IV and RDR2 transcripts, we conducted transcription reactions as in Singh et al., 2019, using 32P to specifically label either the Pol IV or RDR2 strands of their dsRNA products. We then treated the transcripts with S1 nuclease, which digests single-stranded nucleic acids (Figure 6B). For these experiments, Pol IV transcription was initiated using an RNA primer hybridized to a T-less (lacking thymidines) DNA template (see diagram in Figure 6B and Supplementary file 1b for oligonucleotide sequences). Downstream of the primer used to initiate transcription, a non-template DNA oligonucleotide is annealed to the template DNA to induce Pol IV arrest and RDR2 engagement of the Pol IV transcript’s 3′ end, enabling synthesis of the RDR2-transcribed strand (Singh et al., 2019). By 5′ end-labeling the RNA primer, only first-strand RNAs synthesized by Pol IV are labeled and detected. Conversely, by using an unlabeled RNA primer and including alpha-labeled 32P-ATP in the reactions, second-strand RNAs synthesized by RDR2 are specifically labeled. Note that because adenosines are not present in first-strand RNAs generated by Pol IV transcription of a T-less DNA template, first strands are not labeled.
Treatment of transcription reaction products with increasing amounts of S1 nuclease resulted in a progressive 1–2 nt shortening of labeled Pol IV transcripts (Figure 6B, lanes 1–4). Because the 32P label is at the 5′ end of these Pol IV transcripts, S1 trimming must be occurring at their 3′ ends. Body-labeled RDR2 transcripts were shortened by 1 nt upon digestion with S1 nuclease (Figure 6B, lanes 5–8), a result consistent with the removal of the single untemplated nucleotide added to the 3′ end of RDR2 transcripts. Collectively, the data of Figure 6B indicate that the 3′ ends of both the Pol IV transcripts and the RDR2 transcripts overhang the paired strand, fitting the predictions of the model in Figure 6A.
Because Pol IV transcripts within Pol IV-RDR2 dsRNAs are trimmed 1–2 nt by S1 nuclease at their 3′ ends, this suggests that RDR2 initiates 1 or 2 nt internal to the ends of Pol IV transcripts. As a test of this hypothesis, we used a 5′ end-labeled synthetic RNA as a template for recombinant RDR2 (Mishra et al., 2021) and subjected resulting dsRNA products to S1 nuclease digestion (Figure 6C). S1 trimmed the 5′-labeled template strand of the RNA duplexes by 1–2 nt (lanes 1 and 2). As controls, we synthesized RNA strands with perfect complementarity to the template strand (lanes 3 and 4) or were shorter by 1 (lanes 5 and 6) or 2 nt (lanes 7 and 8) at their 5′ ends, annealed these RNAs to the 5′ end-labeled template RNA strand, and subjected resulting dsRNAs to S1 nuclease digestion. S1 digestion patterns for synthetic dsRNAs with 1 or 2 nt overhangs most closely resembled the digestion patterns observed for RDR2 transcription products (compare lanes 2, 6, and 8). We conclude that RDR2 initiates transcription 1 or 2o nt internal to the 3′ ends of the Pol IV transcripts it uses as templates, consistent with the results of Figure 6B, scenarios 2 and 3 of Figure 6A, and the results of a different assay (Fukudome et al., 2021).
Labeling Pol IV and RDR2 transcripts in different ways and observing how siRNAs were derived from the labeled precursor strands revealed that 24 and 23 nt siRNAs come from opposite ends of Pol IV and RDR2-transcribed RNA strands (Figure 6D). By initiating Pol IV-RDR2 transcription with an RNA primer labeled with 32P on its 5′ end, only the Pol IV strands of resulting dsRNAs are labeled (as in Figure 6B). Upon dicing of these dsRNAs into siRNAs, only labeled 24 nt siRNA products are detected (Figure 6D, lane 2). Note that a labeled band of ~25 nt is also apparent in lanes 1 and 2, but this is a Pol IV- and DCL3-independent RNA that results from RDR2 transcription of the 16 nt RNA primer, generating an initial transcript that folds back on itself into a partial stem-loop structure, followed by further elongation upon transcription of the single-stranded portion of the stem.
Unlike the dicing results obtained using end-labeled Pol IV strands, in which only labeled 24 nt siRNAs were detected following dicing, body-labeling of Pol IV or RDR2 strands yields labeled 24 and 23 siRNAs upon dicing (Figure 6D, lanes 4 and 6).
In a third labeling strategy, we initiated Pol IV transcripts using a dephosphorylated primer and unlabeled nucleotide triphosphates to generate Pol IV-RDR2 dsRNAs. Resulting dsRNAs were then diced by DCL3 and dicing products were subjected to capping reactions using vaccinia virus capping enzyme and 32P-GTP to label siRNAs that retain the triphosphorylated 5′ end of the RDR2 strand. In this experiment, only 24 nt siRNAs were capped with 32P-GTP (Figure 6D, lane 10), suggesting that only 24 nt siRNAs come from the 5′ ends of RDR2 transcripts.
Collectively, the experiments of Figure 6D reveal that the siRNAs derived from the 5′ ends of Pol IV transcripts and the 5′ ends of RDR2 transcripts are 24 nt, consistent with DCL3 measuring and cutting 24 nt from the 5′ end of either strand. We thus deduce that the 23 nt siRNA dicing products, which are observed only when Pol IV or RDR2 transcripts are body-labeled (lanes 4 and 6), come from the 3′ ends of Pol IV and RDR2 transcripts.
Discussion
In vivo, highly abundant DCL3-dependent 24 and 23 nt siRNAs map to genomic loci subjected to RNA-directed DNA methylation, forming a swarm of overlapping siRNAs that make relationships between siRNAs of opposite polarity, or alternative size, very difficult to discern. Likewise, determining which siRNA precursor strands are synthesized by Pol IV and which are synthesized by RDR2 is problematic in vivo, as the enzymes work as a complex and are co-dependent (Blevins et al., 2015; Singh et al., 2019). Our in vitro biochemical experiments circumvent these uncertainties. Several important insights came from recapitulating siRNA biogenesis by incubating purified Pol IV, RDR2, and DCL3 with single-stranded bacteriophage M13 as the source of template DNA (Singh et al., 2019). These experiments showed definitively that Pol IV acts first, transcribing the single-stranded template DNA to generate RNA transcripts whose 5′ to 3′ polarity is opposite that of the DNA. RDR2 acts second, using the Pol IV transcripts as templates to generate complementary RNAs that are oriented in the same 5′ to 3′ polarity as the DNA template strand, and thus could not have been directly transcribed from the DNA. This ability to unambiguously discriminate the Pol IV transcripts from the RDR2 transcripts, based on their polarity relative to the template DNA, also allowed strand-specificities of DCL3-diced siRNAs to be assigned. These studies showed that 24 and 23 nt siRNAs are generated from both the Pol IV and the RDR2-transcribed strands of siRNA precursors, with 24 nt siRNAs outnumbering 23 nt siRNAs by several fold (Singh et al., 2019). Moreover, these studies showed that an untemplated nucleotide at the 3′ terminus is a characteristic of RDR2 transcripts, but not Pol IV transcripts. The untemplated 3′-terminal nucleotide persists in siRNAs derived from RDR2 transcripts and, importantly, is selectively enriched among 23 nt siRNAs but not 24 nt siRNAs (Singh et al., 2019). However, 23 nt siRNAs also come from the 3′ ends of Pol IV transcripts (Figure 6D), helping explain why not all 23 nt siRNAs possess an untemplated nucleotide (Wang et al., 2016; Singh et al., 2019).
The occurrence of an untemplated nucleotide at the end of 23 nt siRNAs derived from RDR2 transcripts led to the model in Figure 1B. In this model, diced 23 nt siRNAs derived from the 3′ ends of RDR2 transcripts are paired with 24 nt siRNAs derived from the 5′ ends of Pol IV transcripts. However, the model did not account for how 24 nt siRNAs can come from the RDR2 strand, explain why RDR2-strand 24 nt siRNAs outnumber 23 nt siRNAs, or account for 23 nt siRNAs derived from the Pol IV strand (Singh et al., 2019). Our current study provides the answers to these questions. A key finding is that RDR2 does not initiate second strand RNA synthesis precisely at the 3′ end of a Pol IV transcript, but instead initiates 1–2 nt internal to the transcript, as has also shown by a biochemical study of RDR2 transcription (Fukudome et al., 2021). As a result, a 3′ overhang of 1–2 nt is generated upon RDR2 initiation. A second key finding is that the 3′ overhangs, attributable to RDR2’s mode of initiation or its terminal transferase activity, generate preferred DCL3 substrates. Pol IV-RDR2 dsRNAs are typically shorter than 40 bp, thus they can be diced only once (Blevins et al., 2015; Zhai et al., 2015), either from the left side or right side. Our experiments reveal that this left/right choice is influenced by the 5′ terminal nucleotide at the recessed end of an overhang, with adenosines or uridines favored over guanosines or cytosines. Importantly, Pol IV transcription tends to begin with a purine, with adenosine used somewhat more frequently than guanosine (Singh et al., 2019). Collectively, our experiments suggest that if Pol IV initiates with an adenosine, left-side dicing is favored, with DCL3 generating a diced duplex consisting of a 24 nt siRNA from the 5′ end of the Pol IV strand paired with a 23 nt siRNA from the 3′ end of the RDR2 strand (Figure 7, dicing scenario 1). By contrast, if Pol IV transcription begins with a G, we propose that this encourages right-side dicing, especially if the RDR2 strand has a 5′ A or U that is optimal for DCL3 engagement. Indeed, RNA-seq data have shown that RDR2 transcripts initiate most frequently with A or U (Singh et al., 2019). Moreover, consensus sequences of in vivo precursor 5′ and 3′ ends best matched consensus sequences for siRNA 5′ and 3′ ends by invoking left-side dicing for precursors that initiated with A and right-side dicing for precursors that initiated with G (Blevins et al., 2015). Our current study shows that it is this right-side dicing that accounts for 24 nt siRNAs derived from the RDR2 strand and 23 nt siRNAs derived from the Pol IV strand. There are two ways to generate 24 nt siRNAs from each strand and only one way to generate 23 nt siRNAs from each strand (Figures 6A and 7), which fits with the fact that 24 nt siRNAs are more abundant than 23 nt siRNAs, both in vivo and in vitro (Blevins et al., 2015; Singh et al., 2019).
Our study also provides new insight into the minimum length of Pol IV or RDR2 transcripts that can serve as siRNA precursors or siRNAs. The experiments of Figure 2 indicate that Pol IV transcripts that are 24 nt or longer are sufficient to give rise to DCL3-dependent siRNAs. Moreover, if Pol IV were to generate a 24 nt transcript and RDR2 were to initiate 2 nt internal to this transcript and add an extra untemplated nucleotide to the 3′ end of its transcript, the result would be a 24/23 nt dsRNA duplex. This duplex would not need to be diced by DCL3 because it is already indistinguishable from a dicing product (see Figure 2C). Importantly, in dcl2 dcl3 dcl4 triple mutants, we detected Pol IV-dependent 24 nt RNAs at ~7% of the level found in wild-type plants (Blevins et al., 2015). This may account for low levels of RNA-directed DNA methylation that persist in dicer mutants, suggesting an alternative explanation to putative dicer-independent pathways that have been proposed to be distinct from the canonical RNA-directed DNA methylation pathway yet remain unidentified (Yang et al., 2016; Ye et al., 2016).
Another prior finding in need of reconsideration is the observation that 24 nt siRNAs associated with AGO4 tend to begin with a 5′ adenosine (Mi et al., 2008), a finding that has been interpreted as evidence that AGO4 actively selects siRNAs that begin with an adenosine. We have shown that Pol IV and RDR2 transcripts most frequently begin with adenosine (Singh et al., 2019) and our current study shows that DCL3 preferentially dices substrates with A at their 5′ ends. Thus Pol IV, RDR2, and DCL3 activities could collectively account for much of the 5′ A-enrichment among siRNAs that become loaded into AGO4. Consistent with this interpretation, 24 nt RNAs that have a 5′ A, U, C, or G are assimilated into AGO4 with similar efficiency in vitro (Wang and Pikaard, unpublished).
Our results illustrate how siRNAs of different lengths can be produced by a Dicer endonuclease that measures a fixed length of RNA, an apparent paradox. This fixed length presumably corresponds to the distance from a phosphorylated or triphosphorylated 5′ terminal nucleotide, bound by the PAZ domain, as in other Dicers (Park et al., 2011; Tian et al., 2014; Sinha et al., 2018), to the RNase III domain B cleavage site of DCL3, a length equivalent to 24 nt. The siRNAs of alternative size (23 nt) come from the opposite strand of RNA. The 5′ end of the non-measured strand is generated by the RNase III domain A cleavage reaction, but the 3′ end of this strand is apparently not anchored at a set distance, unlike human dicer (Park et al., 2011), and can overhang the 5′ end of the measured strand by 1 nt, to yield a 23 nt siRNA, or 2 nt to yield a 24 nt siRNA (see Figures 6A and 7). A structural study of DCL3 bound to dsRNA, published while our study was in review, supports this interpretation (Wang et al., 2021).
Our study adds to our knowledge of how individual reactions within the RNA-directed DNA methylation pathway specify what will happen in the next step. Pol IV and RDR2 physically interact (Mishra et al., 2021) and their activities are tightly coupled such that Pol IV transcripts are channeled directly to RDR2 (Singh et al., 2019) by a proposed mechanism involving Pol IV backtracking (Fukudome et al., 2021). Lack of Pol IV-interaction explains why the other five RNA-dependent RNA polymerases in Arabidopsis thaliana are not redundant with RDR2. Our current study now shows that idiosyncrasies of Pol IV and RDR2 transcription, including the choice of nucleotide used to initiate transcription, the alternative start site positions of RDR2 internal to Pol IV transcripts, and the non-templated addition of a nucleotide at the 3′ end of RDR2 transcripts, collectively inform DCL3 how to dice the dsRNAs at the next step of the pathway. This results in production of siRNAs that can be 23 or 24 nt, derived from either strand. Experimental evidence indicates that the specification of 23 and 24 nt siRNAs has functional significance, by specifying that the 24 nt RNA of a 23/24 nt RNA duplex becomes the siRNA stably associated with AGO4 (Wang and Pikaard, unpublished). Thus, by generating 24 and 23 nt siRNAs from either strand and either end of the Pol IV-RDR2 dsRNAs, DCL3 can diversify the set of siRNAs that can be derived from short dsRNA precursors, enabling siRNA base-pairing with Pol V transcripts transcribed from either strand of target loci DNA to maximize DNA methylation.
Materials and methods
Overexpression and purification of recombinant DCL3
Request a detailed protocolRecombinant FLAG epitope-tagged DCL3 was expressed and purified as previously described (Singh et al., 2019). Briefly, High Five insect cells (Thermo Fisher Scientific) grown in Express Five serum-free medium (Thermo Fisher Scientific) were used to produce DCL3 using baculovirus mediated protein expression. The High Five cells that had been infected with bavculovirus at a multiplicity of infection (MOI) of 1.5 were collected, pelleted at 500×g for 5 min, flash-frozen in liquid nitrogen, and stored at –80°C. To purify DCL3, High Five cell pellets (~5 ml) were thawed on ice for 10–15 min and lysed using 45 ml hypertonic lysis buffer (50 mM HEPES-KOH (pH 7.5), 400 mM NaCl, 5 mM MgSO4, 10% glycerol, 2 mM DTT, 1 mM PMSF, 1% protease inhibitor cocktail, and 0.01% IPEGAL CA-630 (Sigma-Aldrich)). The lysate was centrifuged at 50,000×g for 30 min at 4°C. The supernatant was incubated with 0.5 ml anti-FLAG M2 resin (Sigma-Aldrich) for 2 hr at 4°C. Beads were collected by centrifugation for 3 min at 300×g at 4°C and washed 3× as follows: 1× with 50 ml lysis buffer, 1× with 50 ml lysis buffer +5 mM ATP·Mg2+ (Sigma-Aldrich), 1× with 50 ml low salt wash buffer (50 mM HEPES-KOH (pH 7.5), 150 mM NaCl, 5 mM MgSO4, 10% glycerol, 1 mM DTT, and 0.01% IPEGAL CA-630). DCL3-bound FLAG resin was transferred to a small gravity-flow column, washed with 10 ml low salt buffer, and eluted 5× with 1 volume elution buffer (0.5 mg/ml 3× FLAG peptide, 50 mM HEPES-KOH (pH 7.5), 150 mM NaCl, 5 mM MgSO4, 10% glycerol, 1 mM DTT, and 0.01% IPEGAL CA-630). Eluted fractions were pooled and concentrated using a Centricon filter (EMD Millipore) with a 30 kDa cutoff size and analyzed by electrophoresis on a 4%–20% SDS-PAGE gel and Coomassie Blue staining. Recombinant DCL3 was stored at –20°C in storage buffer (50 mM HEPES-KOH (pH 7.5), 150 mM NaCl, 5 mM MgSO4, 45% glycerol, 1 mM DTT, and 0.01% IPEGAL CA-630).
To generate point mutations in the two RNase III domains, the wild-type DCL3 cDNA construct in pUC57 were mutated using site-directed, ligase independent mutagenesis. Resulting constructs were confirmed by sequencing. DCL3 sequences were then sub-cloned into the pFastBac HT B vector (Thermo Fisher Scientific) and bacmids and baculovirus were produced as previously described (Singh et al., 2019). Overexpression and purification of active-site mutants were performed for the wild-type version of the protein.
For the ATP-dependency test of Figure 5, DCL3 was expressed and purified as described above except that 5 mM ATP·Mg2+ was omitted from the second wash. Following the third wash, DCL3-bound anti-FLAG resin was transferred into a small gravity-flow column and washed with 4 ml 1× Apyrase reaction buffer (NEB). The column was then capped and the resin was resuspended in 0.5 ml 1× Apyrase reaction buffer (NEB) containing 10 units of Apyrase protein (NEB). The DCL3-bound FLAG resin was then incubated with Apyrase for 30 min at 24°C. The column was then uncapped and allowed to drain. To remove any residual Ca2+, the resin was then washed with 4 ml EGTA buffer (150 mM NaCl, 50 mM HEPES-KOH (pH 7.5), 20 mM EGTA, 10% glycerol, 1 mM DTT, and 0.01% IPEGAL CA-630). Following the EGTA wash, the DCL3-bound FLAG resin was washed with 10 ml low salt buffer (50 mM HEPES-KOH (pH 7.5), 150 mM NaCl, 5 mM MgSO4, 10% glycerol, 1 mM DTT, and 0.01% IPEGAL CA-630), and eluted with five successive additions of 1 volume of elution buffer (0.5 mg/ml 3× FLAG peptide, 50 mM HEPES-KOH (pH 7.5), 150 mM NaCl, 5 mM MgSO4, 10% glycerol, 1 mM DTT, and 0.01% IPEGAL CA-630). Eluted fractions containing DCL3 were concentrated, analyzed, and stored as described above.
Synthetic nucleic aAcids used in dicing assays
Request a detailed protocolRNA oligonucleotides used in this study were purchased from Integrated DNA Technologies, Inc, except for the 5′ triphosphorylated oligo used in Figure 4C, which was purchased from BioSynthesis, Inc Oligonucleotides used for dicing assays are listed in Supplementary file 1a.
RNA oligonucleotides were gel-purified using 15% denaturing polyacrylamide gel electrophoresis. As necessary, RNAs were monophosphorylated on their 5′ ends using T4 polynucleotide kinase (NEB) and either 3 mM unlabeled ATP (Sigma-Aldrich) or 25 µCi of [γ32P]-ATP, 6000 Ci/mmol (PerkinElmer). Following PNK treatment, reactions were passed through Performa spin columns (EdgeBio) at 1000×g for 3 min. For non-radioactive dicing reactions, equimolar amounts of RNA oligos were mixed in annealing buffer (30 mM HEPES-KOH (pH 7.6), 100 mM potassium acetate), whereas for radioactive dicing reactions a 10% excess of unlabeled oligos were mixed with 32P-labeled oligos in annealing buffer. RNAs were incubated for 5 min in an 85°C water bath and annealed by allowing the water bath to cool to room temperature.
DCL3 dicing assays
Request a detailed protocolDouble-stranded RNA substrates were diced in 40 µl reactions containing either 25 nM (32P-labeled) or 50 nM (nonradioactive) dsRNA substrate, 25 nM DCL3, 50 mM HEPES-KOH (pH 7.5), 150 mM NaCl, 10 mM MgSO4, 10% glycerol, 1 mM DTT, 0.01% IPEGAL CA-630, and 0.4 U/µl RiboLock RNase Inhibitor (Thermo Fisher Scientific) for 30–60 min at room temperature and stopped by incubation at 72°C for 5 min. Reactions were then passed through Performa spin columns (EdgeBio) by centrifuging at 1000×g for 3 min and adjusted to 0.3 M sodium acetate (pH 5.2). 15 µg GlycoBlue(Thermo Fisher Scientific) was added and RNAs were precipitated with 3 volumes of isopropanol at –20°C overnight. Precipitated RNAs were pelleted by centrifugation at 16,000×g for 30 min, washed with 1 ml 70% ethanol, and resuspended in 10 µl 2× RNA Loading Dye (NEB). Resuspended RNAs were incubated at 72°C for 5 min and resolved on 15% polyacrylamide 7 M urea gels. For radioactive assays, gels were transferred to filter paper, vacuum dried, and subjected to phosphorimaging using a Typhoon scanner (GE Healthcare). For non-radioactive assays, gels were stained with SYBR Gold (Invitrogen) in 0.5× TBE for 30 min at room temperature.
For time-course assays, reactions were carried out as described above except 5 nM DCL3 was used and in 80–110 µl reaction volumes, depending on the number of time points assayed. To collect and stop reactions at each time point, 25 µl of all reactions were transferred simultaneously using a multichannel pipette to tubes containing 2 µl 0.5 M EDTA (pH 8.0).
Pol IV-RDR2 transcription assays
Request a detailed protocolTo detect single-stranded regions of dsRNAs synthesized by Pol IV and RDR2, in vitro transcription reactions were performed as previously described (Singh et al., 2019). Oligonucleotides used for in vitro transcription reactions are listed in Supplementary file 1b. Resulting Pol IV-RDR2 dsRNAs were then subjected to the indicated amounts of S1 nuclease digestion for 10 min at 37°C using the manufacturer supplied S1 nuclease digestion buffer. Following this, RNAs were precipitated with 3 volumes of isopropanol, 1/10th volume of 3 M sodium acetate (pH 5.2), 20 µg of GlycoBlue (Thermo Fisher Scientific), and overnight incubation at –20°C. Precipitated RNAs were pelleted by centrifugation 16,000×g for 30 min, washed with 1 ml 70% ethanol, and resuspended in 2× RNA Loading Dye (NEB). Resuspended RNAs were incubated at 72°C for 5 min and loaded onto a 15% polyacrylamide 7 M urea gel. Following electrophoresis, gels were transferred to filter paper, vacuum dried, and subjected to autoradiography using BioMax XAR film (Kodak).
Dicing and capping of siRNAs generated from in vitro Pol IV-RDR2-DCL3 reactions were performed as previously described (Singh et al., 2019).
RDR2 transcription assays
Request a detailed protocolRecombinant RDR2 was expressed and purified as previously described (Blevins et al., 2015).
For RDR2 transcription reactions, recombinant RDR2 (280 nM) was mixed with 5′ 32P-labeled template strand RNA (50 nM) in a 30 µl reaction containing 25 mM HEPES-KOH (pH 7.9), 20 mM ammonium acetate, 50 mM NaCl, 2 mM MgCl2, 0.4 U/µl RNase Inhibitor (NEB), 0.1 mM EDTA, 0.01% Triton X-100, 3% PEG8000, and 0.1 mM each of GTP, CTP, UTP, and ATP. Reactions were incubated at room temperature for 2 hr. Then, the reaction products were split into 10 µl aliquots and S1 nuclease 10× reaction buffer (Promega) was added to 1× concentration. Following this, S1 nuclease (Promega) or an equivalent amount of 1× S1 nuclease reaction buffer was added. After incubation at 37°C for 15 min, reactions were stopped by addition of 50 µl of Proteinase K solution (100 mM Tris-HCl (pH 7.9), 250 mM NaCl, 1 mM MgCl2, 1% SDS, 0.4 mg/ml Proteinase K RNA grade (Invitrogen), 0.3 mg/ml GlycoBlue (Thermo Fisher Scientific)), and incubated for an additional 30 min at 37°C. RNAs were then precipitated by adding 180 µl 100% ethanol and incubating overnight at –20°C. Precipitated RNAs were pelleted by centrifugation at 16,000×g for 20 min, washed with 300 µl 70% ethanol, and resuspended in 10 µl 2× TBE-Urea Sample Buffer (Invitrogen). Resuspended RNAs were incubated at 75°C for 5 min, snap-cooled on ice, and resolved on 15% polyacrylamide 7 M urea gels. Gels were transferred to filter paper, vacuum dried, and subjected to phosphorimaging using a Typhoon scanner (GE Healthcare).
Quantification and statistical analysis
Request a detailed protocolDCL3 assays were quantified using Image Lab version 6.0.1 software (Bio-Rad Laboratories). Substrate and diced RNA bands were boxed and diced band signal was calculated as the % of total signal. In Figure 4, means for triplicate reactions are plotted, with error bars showing the standard error of the mean.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting file; Source Data files have been provided for all figures.
References
-
From guide to target: molecular insights into eukaryotic RNA-interference machineryNature Structural & Molecular Biology 22:20–28.https://doi.org/10.1038/nsmb.2931
-
Evidence for ARGONAUTE4-DNA interactions in RNA-directed DNA methylation in plantsGenes & Development 30:2565–2570.https://doi.org/10.1101/gad.289553.116
-
RNAi and heterochromatin assemblyCold Spring Harbor Perspectives in Biology 7:a019323.https://doi.org/10.1101/cshperspect.a019323
-
RNA-directed DNA methylation: an epigenetic pathway of increasing complexityNature Reviews. Genetics 15:394–408.https://doi.org/10.1038/nrg3683
-
Distinct substrate specificities of Arabidopsis DCL3 and DCL4Nucleic Acids Research 42:1845–1856.https://doi.org/10.1093/nar/gkt1077
-
Origins and evolution of eukaryotic RNA interferenceTrends in Ecology & Evolution 23:578–587.https://doi.org/10.1016/j.tree.2008.06.005
-
Reconstitution of siRNA Biogenesis In Vitro: Novel Reaction Mechanisms and RNA Channeling in the RNA-Directed DNA Methylation PathwayQuantitative Biology (Beijing, China) 1:195–201.https://doi.org/10.1101/sqb.2019.84.039842
-
Dicer uses distinct modules for recognizing dsRNA terminiScience (New York, N.Y.) 359:329–334.https://doi.org/10.1126/science.aaq0921
-
Homodimeric structure and double-stranded RNA cleavage activity of the C-terminal RNase III domain of human dicerJournal of Molecular Biology 374:106–120.https://doi.org/10.1016/j.jmb.2007.08.069
-
Genome-wide analysis of single non-templated nucleotides in plant endogenous siRNAs and miRNAsNucleic Acids Research 44:7395–7405.https://doi.org/10.1093/nar/gkw457
-
Mechanism of siRNA production by a plant Dicer-RNA complex in dicing-competent conformationScience (New York, N.Y.) 374:1152–1157.https://doi.org/10.1126/science.abl4546
-
The RNAs of RNA-directed DNA methylationBiochimica et Biophysica Acta. Gene Regulatory Mechanisms 1860:140–148.https://doi.org/10.1016/j.bbagrm.2016.08.004
-
RNA polymerase V transcription guides ARGONAUTE4 to chromatinNature Genetics 41:630–634.https://doi.org/10.1038/ng.365
-
ARGONAUTE4 control of locus-specific siRNA accumulation and DNA and histone methylationScience (New York, N.Y.) 299:716–719.https://doi.org/10.1126/science.1079695
Article and author information
Author details
Funding
National Institutes of Health (GM077590)
- Craig S Pikaard
Howard Hughes Medical Institute (Investigator:Pikaard)
- Akihito Fukudome
- Vibhor Mishra
- Feng Wang
- Craig S Pikaard
Indiana University (Carlos O. Miller Graduate Student Fellowship)
- Andrew Loffer
- Jasleen Singh
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank the Drosophila Genome Resource Center at Indiana University for use of their insect cell culture facilities. This research was supported by NIH Grant GM077590 and funds to CSP as an Investigator of the Howard Hughes Medical Institute. AL and JS were supported, in part, by Carlos O Miller graduate fellowships at Indiana University.
Copyright
© 2022, Loffer et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,219
- views
-
- 208
- downloads
-
- 26
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Chromosomes and Gene Expression
The mRNA 5'-cap structure removal by the decapping enzyme DCP2 is a critical step in gene regulation. While DCP2 is the catalytic subunit in the decapping complex, its activity is strongly enhanced by multiple factors, particularly DCP1, which is the major activator in yeast. However, the precise role of DCP1 in metazoans has yet to be fully elucidated. Moreover, in humans, the specific biological functions of the two DCP1 paralogs, DCP1a and DCP1b, remain largely unknown. To investigate the role of human DCP1, we generated cell lines that were deficient in DCP1a, DCP1b, or both to evaluate the importance of DCP1 in the decapping machinery. Our results highlight the importance of human DCP1 in decapping process and show that the EVH1 domain of DCP1 enhances the mRNA-binding affinity of DCP2. Transcriptome and metabolome analyses outline the distinct functions of DCP1a and DCP1b in human cells, regulating specific endogenous mRNA targets and biological processes. Overall, our findings provide insights into the molecular mechanism of human DCP1 in mRNA decapping and shed light on the distinct functions of its paralogs.
-
- Chromosomes and Gene Expression
- Computational and Systems Biology
Genes are often regulated by multiple enhancers. It is poorly understood how the individual enhancer activities are combined to control promoter activity. Anecdotal evidence has shown that enhancers can combine sub-additively, additively, synergistically, or redundantly. However, it is not clear which of these modes are more frequent in mammalian genomes. Here, we systematically tested how pairs of enhancers activate promoters using a three-way combinatorial reporter assay in mouse embryonic stem cells. By assaying about 69,000 enhancer-enhancer-promoter combinations we found that enhancer pairs generally combine near-additively. This behaviour was conserved across seven developmental promoters tested. Surprisingly, these promoters scale the enhancer signals in a non-linear manner that depends on promoter strength. A housekeeping promoter showed an overall different response to enhancer pairs, and a smaller dynamic range. Thus, our data indicate that enhancers mostly act additively, but promoters transform their collective effect non-linearly.