1. Cell Biology
  2. Genetics and Genomics
Download icon

Endogenous RNA interference is driven by copy number

  1. Cristina Cruz
  2. Jonathan Houseley  Is a corresponding author
  1. The Babraham Institute, United Kingdom
Research Article
  • Cited 15
  • Views 1,753
  • Annotations
Cite this article as: eLife 2014;3:e01581 doi: 10.7554/eLife.01581


A plethora of non-protein coding RNAs are produced throughout eukaryotic genomes, many of which are transcribed antisense to protein-coding genes and could potentially instigate RNA interference (RNAi) responses. Here we have used a synthetic RNAi system to show that gene copy number is a key factor controlling RNAi for transcripts from endogenous loci, since transcripts from multi-copy loci form double stranded RNA more efficiently than transcripts from equivalently expressed single-copy loci. Selectivity towards transcripts from high-copy DNA is therefore an emergent property of a minimal RNAi system. The ability of RNAi to selectively degrade transcripts from high-copy loci would allow suppression of newly emerging transposable elements, but such a surveillance system requires transcription. We show that low-level genome-wide pervasive transcription is sufficient to instigate RNAi, and propose that pervasive transcription is part of a defense mechanism capable of directing a sequence-independent RNAi response against transposable elements amplifying within the genome.


eLife digest

Genes contain the codes that are needed to make the proteins used by cells. This code is transcribed to make a messenger RNA molecule that is then translated to make a protein. However, other types of RNA called non-coding RNA molecules can disrupt this process by binding to messenger RNA molecules, with matching sequences, before translation begins. This phenomenon, which is known as RNA interference, involves enzymes called Dicer and Argonaute.

Many cells contain large numbers of non-coding RNA molecules—so called because they are not translated to produce proteins—and many of these are capable of starting the process of RNA interference. However, most do not, and the reasons for this are not understood. Now, work by Cruz and Houseley has provided new insight into this phenomenon by showing that it is related to the number of copies of the gene encoding such RNAs in the genome.

Yeast cells normally do not have the genes for RNA interference, but Cruz and Houseley used genetically engineered yeast cells containing Dicer and Argonaute. Although most of the messenger RNA molecules in these cells showed no change, the expression of some genes with high ‘copy numbers’ was reduced. Further experiments that involved adding more and more copies of other genes showed that RNA interference could selectively target messenger RNA molecules produced from genes with an increased copy number—particularly if the copies of the genes were clustered in one location in the genome.

RNA interference is also used to defend against DNA sequences that invade and multiply within a genome, such as viruses and other ‘genetic parasites’. As such, the effect observed by Cruz and Houseley could explain why entire genomes are often continuously copied to RNA at low levels. This activity would allow the monitoring of the genome for the invasion of any genetic parasites that had multiplied to high numbers. Following on from this work, the next challenge will be to understand how gene copy number and location are balanced to achieve a selective RNA interference system.



Over the past decade, our understanding of the complexity of the eukaryotic transcriptome has been revolutionized. Genome-wide sequencing studies in many organisms have revealed that protein-coding mRNAs are augmented by a multitude of non-protein coding RNAs (ncRNAs), many produced from regions of the genome traditionally considered to be transcriptionally silent (The ENCODE Project Consortium, 2012; Bertone et al., 2004; Cheng et al., 2005; David et al., 2006; Birney et al., 2007). Functional data for the vast majority of ncRNAs are currently lacking, with only a few examples characterized in any detail; however, the diversity of mechanisms by which these act suggests that ncRNAs have a rich and varied biology that is largely still to be sampled.

Long ncRNAs which overlap protein-coding genes have the potential to modulate the expression of their cognate coding RNA. Early characterized examples in yeast were thought to work by directly disrupting transcription factor or polymerase binding to the promoter of the coding RNA (Martens et al., 2004; Hongay et al., 2006); however, more recent data implicate specific chromatin structure changes in repression (Gelfand et al., 2011; Hainer et al., 2011), and many other cases of ncRNAs that alter chromatin modifications have been described (Camblong et al., 2007; Berretta et al., 2008; Houseley et al., 2008; Pinskaya et al., 2009; van Werven et al., 2012). Chromatin modifications are not necessarily repressive, and ncRNAs that enhance expression of their overlapping coding gene have also been described (Uhler et al., 2007; Hirota et al., 2008). In these examples, chromatin modifications are deposited during transcription, and therefore the act of transcription rather than the ncRNA itself is important. This is not always the case, and in higher eukaryotes multiple cis-acting ncRNAs have also been characterized, particularly as functional agents in imprinting. For example, Air and Kcnq1ot1 act in cis to deposit repressive chromatin marks and DNA methylation, but these ncRNAs interact with chromatin modifiers and allele specificity is achieved by restriction of the ncRNA to the vicinity of the transcription site, although the importance of the transcription itself remains controversial (Nagano et al., 2008; Pandey et al., 2008; Redrup et al., 2009; Latos et al., 2012).

Genomes are also replete with low abundance and unstable RNA. The vast majority of ncRNAs in budding yeast are unstable (Neil et al., 2009; van Dijk et al., 2011), limiting the potential action of the RNAs themselves, although the transcription of such RNAs can still alter gene expression (reviewed in Houseley, 2012). Such unstable RNAs are also widespread in higher eukaryotes, probably with similar functional roles (Chekanova et al., 2007; Preker et al., 2008). More mysterious is the pervasive transcription that permeates eukaryotic genomes; the ENCODE project found that almost all the human genome is transcribed at some point, but the products of this transcription are vanishingly rare (Cheng et al., 2005; Birney et al., 2007; Kapranov et al., 2007; Goodman et al., 2012). It appears that regions of the genome which are not actively transcribed for other reasons undergo pervasive transcription; however, it is not known whether this pervasive transcription simply represents transcriptional noise or whether the transcription or RNAs themselves have important but as yet undiscovered functions.

Systems in which a ncRNA is transcribed antisense to a sense protein-coding RNA are common and have strong regulatory potential (Figure 1A) (Derrien et al., 2012; Carninci et al., 2005; Xu et al., 2009). It has been suggested that, since antisense ncRNAs are perfectly complementary to their cognate mRNA, the two species could form double stranded RNA (dsRNA) that would be a substrate for the RNA interference system (RNAi). During a basic RNAi response, dsRNA is cleaved by the endonuclease Dicer into short interfering RNA (siRNA), of which one strand is then loaded onto an Argonaute protein. The Argonaute–siRNA complex can anneal to complementary sequences in target RNAs, which are then cleaved by the endonuclease activity of Argonaute. RNAi was originally discovered in Caenorhabditis elegans and rapidly linked to the phenomenon of post-transcriptional gene silencing in plants; however, almost all eukaryotes contain Dicer and Argonaute orthologues and therefore have some form of RNAi system (Hamilton and Baulcombe, 1999; Fire, 1998; Hannon, 2002). RNAi probably evolved to protect cells against dsRNA viruses, a role which is maintained in plants, insects, and lower eukaryotes (Ding, 2010) and has recently been described in mammalian cells (Li et al., 2013; Maillard et al., 2013). RNAi also forms a potent defense against transposons, and high-copy transposon-derived sequences are excellent targets for RNAi, giving rise to copious siRNAs in most eukaryotes including mammals (Yang and Kazazian, 2006; Slotkin and Martienssen, 2007; Babiarz et al., 2008). In addition to degrading transposon-derived and viral RNA, siRNAs can mediate transcriptional repression of target RNAs through chromatin modifications and DNA methylation, although this activity is seemingly much stronger in lower eukaryotes and plants than in mammals (Martienssen et al., 2005; Lejeune et al., 2010; Zhang and Zhu, 2011). However, the source of the dsRNA that is processed into siRNA is not always obvious, nor is the mechanism by which cells differentiate host and transposon-derived sequences. siRNA-mediated repression is complemented by PIWI-interacting RNAs (piRNAs), which bind to the Argonaut-related PIWI-domain proteins and enforce transposon repression in the germline of eukaryotes from worms to mammals (Siomi et al., 2011). piRNAs are derived from specific genomic clusters, but it is unclear how the transcripts from these clusters are selected for processing into primary piRNAs and many of the processing enzymes remain to be identified (reviewed in Ishizu et al., 2012).

Frequency of annotated antisense non-protein coding RNAs (ncRNAs) and effects on mRNA abundance.

(A) Schematic of an example sense mRNA-antisense (ncRNA) system. (B) Number of annotated open reading frames (ORFs) with antisense transcripts. Positions of CUTs, SUTs, and XUTs were collated with expressed ORFs (Xu et al., 2009; van Dijk et al., 2011), SUTs were later re-classified as XUTs were removed. Overlaps between ORFs expressed in glucose media (total 5171, Xu et al., 2009) and other RNAs were calculated and summed for increasing minimum overlaps of 50–500 bp. ORF–ORF overlaps and ORF–ncRNA overlaps were analyzed separately as ORF–ORF overlaps are consistently smaller. Detailed figures are given in Table 1. (C) Abundance of short interfering RNAs (siRNAs) in RNA interference (RNAi)+ strain produced from expressed ORFs with and without an annotated overlapping antisense ncRNA, based on read counts from published high-throughput sequencing data (Drinnenberg et al., 2011). Minimum antisense overlap with ORF was set at 250 bp; only ORFs with >100 reads in the wild-type poly(A)+ library were assessed to remove noise. Stated p value calculated by Student’s t test. (D) Abundance of mRNA in RNAi+ cells relative to wild-type; data source and categories as in C, differences were not significant.


Various studies have looked for endogenous sense-antisense RNA pairs that instigate RNAi responses. Efficient generation of siRNAs from endogenous sense-antisense systems (endo-siRNA) has been observed in plants under stress (Borsani et al., 2005; Katiyar-Agarwal et al., 2006), and mammalian oocytes generate endo-siRNAs that can mediate mRNA knockdown (Tam et al., 2008; Watanabe et al., 2008). However, although endo-siRNAs have been detected outside the germline in mammals, they are surprisingly under-represented where sense and antisense RNAs are co-expressed (Faghihi and Wahlestedt, 2006; Okamura et al., 2008; Carlile et al., 2009), and overall there is a positive correlation between antisense and sense RNA expression in mammalian genomes, which is inconsistent with RNAi (Derrien et al., 2012; Katayama et al., 2005). This raises questions about whether endogenous sense-antisense systems do in fact form dsRNA in vivo and, if so, whether all dsRNA is equivalently accessible to Dicer.

The tight integration of the RNAi system into the physiology of most eukaryotic cells makes it very difficult to disentangle direct and indirect effects of mutating RNAi components (reviewed in Ketting, 2011). To elucidate factors important for the induction of RNAi by endogenous sense-antisense systems, we therefore used a recently described synthetic system in which RNAi is reconstituted in Saccharomyces cerevisiae by the introduction of Argonaute and Dicer from the related yeast S. castellii (Drinnenberg et al., 2009). S. cerevisiae is highly unusual in lacking an endogenous RNAi system, allowing maintenance of the symbiotic dsRNA Killer virus (Drinnenberg et al., 2011). The reconstituted system is functional, since RNAi+ S. cerevisiae efficiently degrades exogenous hairpin RNAs and endogenous Ty retrotransposon transcripts; however, no clear mRNA expression changes are detectable in these cells (Drinnenberg et al., 2009, 2011).


Antisense ncRNAs exist for many S. cerevisiae genes; combining published datasets we found that 15–30% of yeast open reading frames (ORFs) have an annotated antisense ncRNA depending on the minimum size of overlap considered (Xu et al., 2009; van Dijk et al., 2011), not including overlapping convergent ORFs (Figure 1B and Table 1). In RNAi+ cells, these ORFs produce more siRNA than ORFs lacking an antisense (Figure 1C), showing that they are transcribed into dsRNA that is targeted by the RNAi machinery. However, these ORFs do not show reduced mRNA levels in RNAi+ cells consistent with published data, suggesting that insufficient siRNAs are produced to elicit a detectable mRNA knockdown (Derrien et al., 2012; Katayama et al., 2005; Drinnenberg et al., 2011) (Figure 1D).

Table 1

Stability of antisense ncRNAs

Overlap size (bp)Overlap typeTotals
50700106657552215961448 (91%)
100449100854347515071367 (91%)
15021696750842514231306 (92%)
20013693147840313581249 (92%)
2509689343438012871181 (92%)
3006981239136311891085 (91%)
350627593353511086990 (91%)
40054694292334992899 (91%)
45051637244322904811 (90%)
50048591204302828741 (89%)
  1. Number of ORFs overlapping with ORFs and various classes of ncRNAs, with various minimum size cut-offs for the overlapping region.

  2. Totals are given for ORFs overlapping with ncRNAs and with unstable ncRNAs, including a percentage of overlapping ncRNAs that are unstable.

  3. ncRNA: non-protein coding RNA; ORF: optical reading frame; XUT: Xrn1-sensitive unstable transcript (degraded in cytoplasm); CUT: cryptic unstable transcript (degraded by nuclear exosome); SUT: stable unannotated transcript (not known to be degraded).

We first asked whether any endo-siRNA pairs are degraded by RNAi in this reconstituted system. RNAi+ cells produce abundant siRNAs from sub-telomeric Y′ elements and from the ribosomal DNA (rDNA) intergenic spacers (Drinnenberg et al., 2011 and Figure 2—figure supplement 1) and, despite transcriptional repression by the histone deacetylase Sir2, both regions transcribe sense and antisense ncRNAs that could hybridize to form dsRNA (Aparicio et al., 1991; Yamada et al., 1998; Kobayashi and Ganley, 2005; Houseley et al., 2007) (Figure 2A,D). Northern blots revealed full-length ncRNAs from both strands of the Y′ elements in wild-type cells (Figure 2B lanes 1,5 marked with arrows); these were largely absent in the RNAi+ strain, being replaced by heterogeneous degradation products and readily detectable siRNAs (Figure 2B,C). Despite weak transcriptional repression in this genetic background, ncRNAs and siRNAs were more abundant in sir2Δ mutants reinforcing the precursor–product relationship (Figure 2B,C). Equivalent results were seen for the rDNA intergenic spacer region (Figure 2D–F). These data show that endo-siRNA pairs can form RNAi substrates and undergo efficient degradation by a minimal RNAi system.

Figure 2 with 1 supplement see all
High-copy endogenous sense-antisense pairs instigate efficient RNA interference (RNAi).

(A) Schematic diagram of sub-telomeric Y′ elements. (B) Northern analysis of Y′ element non-protein coding RNAs (ncRNAs) comparing wild-type and RNAi+ strains in SIR2 and sir2Δ backgrounds. 18S ribosomal RNA is shown as a loading control. Arrows indicate full-length RNA species. (C) Northern analysis of Y′ element-derived short interfering RNAs (siRNAs) from cells in B, tRNAs are shown as a loading control. (DF) Equivalent analysis of rDNA intergenic spacer ncRNAs.


One distinguishing feature of these regions is high copy number; to determine whether copy number amplification can drive RNAi, we examined MAL32 which has a clearly defined antisense RNA that is co-expressed with the sense mRNA (Figure 3—figure supplement 1). MAL32 is effectively present at two copies in the haploid genome as the orthologous gene MAL12 shares 99.5% nucleotide identity, reducing potential transcriptional repression. As for approximately 90% of yeast antisense RNAs (Table 1), MAL32 antisense RNA is highly unstable and is only detectable in strains lacking the nuclear exosome co-factor Trf4 (reviewed in Houseley and Tollervey, 2009) (Figure 3—figure supplement 1), but endogenous MAL32 mRNA was not down-regulated in the RNAi+ strain even in trf4Δ cells (Figure 3B and Figure 3—figure supplement 2A). When expressed from a high-copy plasmid, however, MAL32 mRNA was significantly down-regulated by RNAi (Figure 3C) with concurrent production of siRNA (Figure 3D). The MAL32 antisense RNA was only detected in these experiments as a smear of degradation products and was not noticeably altered by RNAi, probably because nuclear degradation acts faster than RNAi on this substrate. To confirm that the knockdown of MAL32 mRNA was not an indirect effect of the strain background or an undirected Argonaute cleavage, we reconstituted the RNAi system in the BY4741 background using separate plasmids expressing Dicer and Argonaute. dsRNA from the MAL32 locus was detectable in these cells and was removed by Dicer; however, a significant knockdown of the mRNA was only observed in cells expressing both Dicer and Argonaute, confirming that the knockdown represents a genuine RNAi response (Figure 3—figure supplement 3).

Figure 3 with 6 supplements see all
Copy number amplification of coding genes can instigate RNA interference (RNAi).

(A) Schematic of MAL32 locus. (B) MAL32 mRNA and antisense non-protein coding RNA (ncRNA) comparing wild-type and RNAi+ strains in wild-type and trf4Δ backgrounds; cells grown on YP raffinose (extended image and quantification shown in Figure 3—figure supplement 2A). (C) mRNA and antisense ncRNA from MAL32 cloned onto a high-copy plasmid in wild-type and RNAi+ strains. Lanes 3,4 show empty vector control. Antisense panel shows degradation products, no full-length antisense is detectable due to Trf4 activity. (D) Short interfering RNA (siRNA) analysis of cells from C. (E) Schematic of GAL4 locus. (F) GAL4 mRNA and antisense ncRNA in wild-type and RNAi+ strains; cells grown in YP galactose (extended image and quantification shown in Figure 3—figure supplement 2B). (G) mRNA and antisense ncRNA from GAL4 locus cloned onto a high-copy plasmid in wild-type and RNAi + strains. Lanes 5,6 show empty vector, signal is from genomic GAL4; note that cells used here are diploids to mitigate defects in galactose response (see ‘Materials and methods’). Lanes 3,4 show a previously described GAL4 antisense mutant (Geisler et al., 2012); this removes detectable antisense RNA for genomic GAL4, but the mutant sequence still expresses an antisense ncRNA when cloned on the high-copy plasmid (see Figure 3—figure supplement 4). (H) siRNA analysis of cells in G. For quantification, n = 4 biological replicates, error bars represent ± 1se, *p<0.05, ***p<0.01 by Student’s t test, y axes in arbitrary units.


We then tested GAL4, which is a single-copy gene with a co-expressed antisense that is degraded in the cytoplasm by Xrn1 (Geisler et al., 2012) (Figure 3E). Cells lacking Xrn1 show increased levels of antisense and, unexpectedly, sense RNA, but, as for MAL32, we did not detect a significant decrease in full-length RNA in RNAi+ xrn1Δ cells (Figure 3F and Figure 3—figure supplement 2B). However, amplification of the locus by cloning on a high-copy plasmid leads to significant degradation of the sense RNA along with the antisense RNA by RNAi (Figure 3G,H). Surprisingly, this occurred even in a known mutant that lacks antisense expression (Geisler et al., 2012), but 5′ RACE (Rapid Amplification of 5′ Complementary DNA Ends) experiments revealed that antisense RNA is still produced by this mutant when expressed from a high-copy plasmid, even if it is too heterogeneous for detection by northern blot (Figure 3—figure supplement 4). Taken together, these experiments on MAL32 and GAL4 demonstrate that increasing gene copy number can make the products of a normal gene susceptible to RNAi.

One potential confounding factor in these experiments is the copy number of the high-copy plasmids; if the copy number drops in RNAi+ cells, this would provide a trivial explanation for the observed knockdowns. However, Southern blotting revealed that RNAi+ cells contained approximately twofold more plasmid than the controls, which would tend to decrease rather than increase the apparent effect of RNAi (Figure 3—figure supplement 5A,B). It is likely that RNAi degrades the mRNA for the plasmid-encoded selectable marker and 2µ maintenance genes (2µ genes produce copious siRNA, Figure 3—figure supplement 5C), and the plasmid copy number rises to compensate for this. We also wanted to use a different method to confirm the northern blot results. We therefore lysed wild-type and RNAi+ cells containing the MAL32 and GAL4 plasmids, precipitated dsRNA using a specific monoclonal antibody (Schonborn et al., 1991; Gullerova and Proudfoot, 2012), and assayed total and dsRNA levels by quantitative RT-PCR. MAL32 mRNA was knocked down approximately 75%, as observed by northern blot analysis, while dsRNA was reduced 11-fold, consistent with specific removal of dsRNA by Dicer. GAL4 mRNA knockdown was measured at 80% in this assay, again with an 11-fold reduction in dsRNA in the RNAi+ strain (Figure 3—figure supplement 6).

Increasing gene copy number also increases RNA production. To separate the contributions of RNA abundance and copy number, we analyzed existing genome-wide data (Hobson et al., 2012; Drinnenberg et al., 2011). If siRNA formation depends only on precursor RNA abundance, a positive correlation should be observed between total RNA abundance and siRNA abundance, and there should be no difference between distributions of single-copy loci and multi-copy loci. We observed little evidence for such a positive correlation (Figure 4—figure supplement 1); however, plots of siRNA versus total RNA abundance for single-copy and multi-copy loci showed strikingly different distributions, with multi-copy loci clearly biased towards higher siRNA production (Figure 4A). To quantify this difference, loci were segregated into eight groups of increasing total RNA abundance and siRNA abundance was assessed for single-copy and multi-copy loci in each group (Figure 4B). siRNA production was significantly higher from multi-copy loci than from single-copy loci in all except the lowest category of RNA abundance. This result was robust to changes in the threshold between low and high copy, and was still observed in a comparison of low to medium copy number, showing that high-copy Ty1 retrotransposons were not distorting the analysis (Figure 4—figure supplement 2). A normalization step is required in these analyses to deal with mapping of sequence reads to multi-copy loci (discussed in detail in ‘Materials and methods’); however, the same differences were observed with no normalization or a different normalization scheme (Figure 4—figure supplement 3). These surprising results show that multi-copy loci produce more siRNA than single-copy loci with equivalent RNA abundance. If this observation is real, the siRNA:total RNA ratio should be predictive of copy number, an important test since this comparison requires no copy number normalization. As predicted, the 1% of genome with the highest siRNA:total RNA ratio is massively enriched for multi-copy loci (Figure 4C), and when this ratio was plotted across a chromosome, an obvious correlation was observed between regions of high-copy number and regions with high siRNA:total RNA abundance (Figure 4D). These analyses clearly demonstrate that selectivity towards the products of multi-copy loci is an emergent property of a minimal RNAi system.

Figure 4 with 3 supplements see all
Multi-copy loci are preferentially targeted by RNA interference (RNAi).

(A) Short interfering RNA (siRNA) (Drinnenberg et al., 2011) and total RNA (Hobson et al., 2012) abundance for loci with copy number <2 (left, single-copy) or ≥2 (right, multi-copy). (B) Quantification of data from A binned into categories of increasing total RNA level, with p values for pairwise comparisons of siRNA abundance in single-copy and multi-copy datasets using the Wilcoxon Rank Sum test. (C) Copy number distribution of the 1% of loci with the highest siRNA:total RNA ratio compared with other loci; difference is significant by Wilcoxon Rank Sum test, p<2.2 × 10−16, loci scoring below noise threshold (0–2 category in B) were removed. n values for tests in B and C are given in Table 2. (D) Comparison of copy number with siRNA:total RNA ratio across chromosome I.


We then directly tested the effect of copy number at the MAL32 locus. We constructed strains in which MAL32 sense and antisense RNAs were expressed at similar levels from multi-copy or single-copy loci by over-expressing single-copy sense and antisense (Figure 5A). In this system, both sense and antisense RNAs were produced at higher levels from the single-copy system (Figure 5B compare lanes 1 and 3) but more siRNAs were produced from the multi-copy system (Figure 5C compare lanes 2 and 4). The over-expression of both RNAs from the single-copy MAL32 locus led to the production of easily detectable siRNA, as would be expected; however, this result directly demonstrates that gene copy number influences the formation of siRNA above and beyond the effect of total RNA abundance. The increased siRNA production in these cells is most likely due to enhanced dsRNA formation in the multi-copy system. To confirm this, we quantified MAL32 RNA in wild-type cells that is resistant to the single-strand specific nuclease RNase A, and observed significantly more RNase-resistant material in cells expressing MAL32 from the multi-copy system than the single-copy system (Figure 5D). This experiment shows that a multi-copy locus produces more dsRNA than an equivalently expressed single-copy locus in wild-type cells without the RNAi system, explaining the increased siRNA formation in RNAi+ cells.

Single gene analysis of copy number effect on RNA interference (RNAi).

(A) Schematic of single-copy and multi-copy MAL32 system. (B) Northern analysis of MAL32 RNA from single-copy and multi-copy systems. All visible species are included in antisense quantification. (C) MAL32 short interfering RNA (siRNA) abundance from cells in B. (D) RNase A sensitivity of MAL32 in wild-type cells expressing multi-copy and single-copy MAL32. Cells were lysed on ice, treated with RNase A as indicated and analyzed by northern blot. 25S and L-A (a double stranded RNA) are shown as controls for loading and RNase specificity. n = 3 biological replicates, error bars ±1se, *p<0.05, ***p<0.01 by Student’s t test, y axes in arbitrary units.


The influence of copy number suggested that dsRNA formation and potentially siRNA production may occur in the nucleus. We initially tested this by immunofluorescence against Dicer and dsRNA in mixed populations of wild-type and RNAi+ cells (Figure 6—figure supplement 1). Dicer was present in cytoplasmic foci, but also showed a diffuse cytoplasmic staining (compare the indicated wild-type and RNAi+ cells) and, under wide-field imaging, appeared to be present in the nucleus. However, super-resolution images of the same cells showed nuclear exclusion of Dicer; therefore, if Dicer is present in the nucleus, it is at low levels. dsRNA staining in these cells revealed many cytoplasmic foci, presumably Killer virus dsRNAs that are known to be incompletely cleared by RNAi (Drinnenberg et al., 2011), but did not show unambiguous nuclear staining.

As an alternative, we asked whether the spatial configuration of gene copies within the nucleus could affect siRNA formation; such an effect would provide strong evidence for the formation of dsRNA in the nucleus. In systems which undergo efficient RNAi such as the rDNA and 2µ plasmids, all gene copies are clustered together in a small sub-nuclear volume. To test the importance of this clustering, we used the sense-antisense system at the TRP1 locus which produces detectable siRNA even when present in the genome in only two tandem copies (data not shown). We generated strains with three tandem copies of TRP1 on a single plasmid (Clustered, Cls) or three unlinked copies (Dispersed, Dsp) (Figure 6A), and expressed Dicer without Argonaute to allow siRNA formation but minimize the effect of RNAi on total RNA levels. Quantification of total RNA showed that both systems produced similar amounts of sense and antisense RNA molecules (Figure 6B), although this experiment was complicated by read-through transcripts of antisense TRP1 from the clustered system (Figure 6B lanes 1,2), a behavior that was not prevented by insertion of transcriptional stop cassettes between the repeats. Nonetheless, the clustered system produced fourfold more TRP1 siRNA than the dispersed system (Figure 6C), showing that close nuclear juxtaposition of transcriptional loci enhances dsRNA formation.

Figure 6 with 1 supplement see all
Clustered loci show higher efficiency of short interfering RNA (siRNA) formation.

(A) Comparison of the systems used. Three copies of TRP1 were placed either in tandem on a single-copy plasmid (Clustered, Cls) or a single copy was left in the genome at the TRP1 locus and two further copies were placed on different single-copy plasmids (Dispersed, Dsp). Dicer was expressed from a single-copy plasmid. (B) Sense and antisense RNA expression in clustered and dispersed systems with and without Dicer. Quantification shows that Dicer alone has little effect on total RNA levels. Read-through species visible in lanes 1,2 are included in the quantification; values have been normalized for the different number of probe binding sites in the read-through RNAs. In the absence of this normalization (i.e., counting the number of binding sites rather than the number of molecules), the clustered antisense is approximately twofold more abundant than the dispersed antisense, which is insufficient to explain the difference in siRNA formation. (C) siRNA produced from TRP1 in clustered and dispersed systems. n = 3 biological replicates, error bars ±1se, ***p<0.01 by Student’s t test, y axes in arbitrary units.


While testing the effects of copy number amplification on siRNA production, we noticed that even low abundance sense-antisense ncRNA pairs (selected from a published dataset, Xu et al., 2009) underwent efficient RNAi when amplified to high copy number. For both the SUT176 and SUT430 systems (Figure 7A), the sense and antisense RNAs are barely detectable by northern blot and are clearly not targeted for degradation by RNAi (Figure 7B). However, after cloning on high-copy plasmids, the full-length RNAs became highly susceptible to RNAi and produced copious siRNAs (Figure 7C,D). This raised the interesting prospect that low abundance pervasive transcription would be sufficient to trigger efficient RNAi responses from sequences that undergo copy number amplification.

RNA interference (RNAi) against transcripts from amplified low-expression systems.

(A) Schematic diagrams of SUT176 and SUT430 loci. (B) Northern analysis of SUT176 and SUT430 transcripts from single-copy genomic loci in wild-type and RNAi+ cells. Ty1 RNA is a positive control for RNAi, ACT1 is a loading control. (C) Analysis of SUT176 and SUT430 non-protein coding RNAs (ncRNAs) expressed from high-copy plasmids in wild-type and RNAi+ cells. Amplified regions are indicated in A. (D) Short interfering RNA analysis of cells in C.


Clear examples of pervasive transcription are not well defined in any organism because, by definition, the products of pervasive transcription are almost undetectable. We therefore chose to examine the GAL gene cluster (Figure 8A), which is tightly transcriptionally repressed in cells grown in glucose. Under these conditions, antisense ncRNAs are produced from the GAL10 ORF with a known steady-state abundance of one RNA molecule per 14 cells (Houseley et al., 2008; Pinskaya et al., 2009) (arrows in Figure 8B lane 1). Transcription of these ncRNAs is abrogated in a previously described Reb1 binding site mutant (RBSΔ), leaving almost no detectable RNAs from this locus (Figure 8B lane 3). For reasons that remain unclear, the GAL cluster is slightly de-repressed in the RNAi+ strain (Figure 8B compare lanes 1,2); nonetheless, the RBSΔ RNAi+ strain (Figure 8B lane 4) only produces very low level heterogeneous transcripts from GAL10, suggesting that it forms a good model of pervasive transcription. Cloning either wild-type or RBSΔ GAL clusters onto high-copy plasmids substantially increased the levels of detectable ncRNA as expected (Figure 8B lanes 5,7), and these ncRNAs were processed into easily detectable siRNAs (Figure 8C lanes 6,8). Therefore, ncRNAs produced at the level of pervasive transcription are sufficient to mediate extensive siRNA production when the copy number of the transcribing locus is increased.

Figure 8 with 1 supplement see all
RNA interference (RNAi) against pervasive transcripts from the repressed GAL cluster.

(A) Schematic representation of the GAL cluster. (B) Non-strand-specific northern blot of non-protein coding RNAs (ncRNAs) produced from the GAL locus present at single-copy (lanes 1–4) or high-copy (lanes 5–8), showing wild-type and RBSΔ mutant. Arrow indicates GAL10 antisense RNA. Strand-specific northern blots for the same RNA are shown in Figure 8—figure supplement 1. (C) GAL10 short interfering RNA (siRNA) from the same cells as in B. (D) Expression of GAL10 mRNA from a single-copy genomic locus under the control of a Cu2+-responsive promoter in wild-type and RNAi+ strains carrying an empty vector (lanes 1,2), high-copy wild-type GAL cluster (lanes 3,4), or high-copy RBSΔ GAL cluster (lanes 5,6). n = 3 biological replicates, error bars ±1se, ***p<0.01 by Student’s t test, y axes in arbitrary units.


The siRNAs produced from the high-copy GAL10 locus are clearly sufficient to degrade the GAL10 ncRNAs in the RNAi+ background (Figure 8B compare lanes 5,7 with lanes 6,8); however, a classical RNAi response should be able to degrade RNA expressed from a separate locus. To test this we introduced the high-copy GAL cluster plasmids into a strain in which the single-copy genomic GAL10 ORF is expressed at high levels from a Cu2+-dependent promoter, allowing expression of the GAL10 mRNA from the single-copy locus while the GAL clusters present on the high-copy plasmids remain fully repressed. As observed for the GAL10 ncRNAs, the GAL10 mRNA was expressed at higher levels in the RNAi+ strain than in the wild-type (Figure 8D compare lanes 1,2) but, nonetheless, both wild-type and RBSΔ high-copy GAL cluster plasmids caused highly significant >50% knockdowns of the GAL10 mRNA compared with the empty vector control (Figure 8D lanes 2,4,6). This was not an indirect effect of the high-copy GAL clusters alone as, in the wild-type background, GAL10 mRNA levels were slightly increased by the presence of the GAL plasmids (Figure 8D lanes 1,3,5). These data demonstrate that pervasive transcription of a high-copy locus is sufficient to instigate an effective RNAi response that can mediate the degradation of a target mRNA in trans.


The ability of the RNAi system to selectively target the products of high-copy sequences such as transposons provides a remarkably efficient genome defense, as well as an effective way to differentiate heterochromatic regions, which are often repetitive, from gene-rich euchromatin. Here we have demonstrated that RNAi has an innate preference for the products of high-copy sequences, probably because RNA from high-copy sequences forms dsRNA more efficiently. It has long been known that cells can recognize and silence high-copy DNA, which would form a basic defense against uncontrolled amplification of transposable elements (reviewed in Hsieh and Fire, 2000). This of course requires a mechanism to count copy number, or at least differentiate high- and low-copy regions, which has remained mysterious. Our data show that RNAi provides such a mechanism by selectively targeting the products of high-copy loci.

The production of siRNA from high-copy DNA, which would be the basis of such a counting mechanism, absolutely requires that all DNA is transcribed; if this does not occur, transposable elements that remain transcriptionally silent would be invisible to the system. Pervasive transcription, the general background of very low level RNA produced across the genome, ensures that the vast majority of the genome is transcribed, and therefore that no region remains completely silent (Cheng et al., 2005; Birney et al., 2007; Kapranov et al., 2007; Goodman et al., 2012). The extent to which mammalian genomes are pervasively transcribed has been controversial; however, many of the questions revolve around whether the pervasive transcripts represent defined functional products or whether much of the detected RNA represents random transcriptional noise (van Bakel et al., 2010; Clark et al., 2011). For a general surveillance function, it does not really matter whether pervasive transcription is formed of many discrete transcripts or occurs at random since either process should be sufficient to generate dsRNA. If a large proportion of pervasive transcription does represent random noise, this would be actively advantageous; random transcription would be sequence independent, and therefore transposable elements could not become fully silent by mutating individual promoter sequences. We suggest that the primary function of pervasive transcription lies in ensuring the whole genome is transcribed to allow identification and suppression of transposable elements; this does not conflict with the idea that some proportion of these transcripts may have additional functions.

We propose that hybridization kinetics explains the dependence of RNAi on copy number (shown in Figure 9). The rate of formation of dsRNA from single stranded sense and antisense RNA is proportional to the concentration of each strand of RNA, and so is inversely proportional to the square of the reaction volume. Technically, the reaction volume is the non-excluded volume of the cell; however, this assumes a uniform distribution of RNA throughout the cell. In reality the RNA is far from evenly distributed, so some small volumes may have very high concentrations of RNA and, within these volumes, the rate of dsRNA formation will be dramatically higher than in the bulk of the cell. Single-copy loci cannot simultaneously transcribe sense and antisense RNA (Hobson et al., 2012) so, although such a locus can give rise to a mixed population of sense and antisense RNA in the cytoplasm over time, in the vicinity of the transcription site only one sense of RNA should ever be present, assuming efficient RNA export. Annealing of sense and antisense RNA must therefore occur in the relatively large non-excluded cytoplasmic volume of the cell, which will be inefficient except for very highly expressed RNA. In contrast, at multi-copy loci the simultaneous production of sense and antisense RNA from many closely spaced sites can lead to high concentrations of sense and antisense RNA around the transcription site, causing efficient duplex formation and hence efficient RNAi.

Proposed mechanism for RNA interference (RNAi) on high-copy loci.

The rate of double stranded RNA (dsRNA) formation, the required first step in RNAi, is highly dependent on local RNA concentration. Single-copy loci cannot simultaneously transcribe sense and antisense RNA, allowing RNA export to occur before the strands meet and requiring hybridization to occur in the relatively large cytoplasmic volume. In contrast, sense and antisense RNA can be simultaneously transcribed from different parts of a multi-copy locus and, therefore, if the copies are closely juxtaposed in the genome or in 3D space, the local concentration of sense and antisense RNA around the transcription sites should be high, favoring dsRNA formation.


This mechanism predicts that dsRNA formation should occur in the nucleus, but we were not able to detect Dicer in the nucleus by immunofluorescence. This reflects the situation in higher eukaryotes where Dicer is largely cytoplasmic, but recent experiments in Drosophila and mammalian cells have detected small quantities of nuclear Dicer, particularly associated with chromatin (Sinkkonen et al., 2010; Cernilogar et al., 2011; Gullerova et al., 2011; Doyle et al., 2013). Low but functional levels of Dicer may therefore be present in the nucleus of RNAi+ cells and able to generate siRNA. Alternatively, the dsRNA may be exported and processed in the cytoplasm. To our knowledge, there is no clear evidence for or against export of dsRNA by normal pathways; certainly these would not be too large or too structured to pass through nuclear pores compared, for example, with pre-ribosomes.

One notable prediction of this mechanism is that clustering of multi-copy transcription sites would be a particularly efficient way to increase the local density of sense and antisense RNA. All the systems we have described in this paper are clustered: rDNA repeats are arranged in tandem, telomeres are known to cluster together at various points in the cell cycle (Gotta et al., 1996), and high-copy 2µ plasmids exist in a discrete focus that is vital for copy number maintenance (Velmurugan et al., 2000; Wong et al., 2002). The comparison of siRNA formation from clustered and dispersed TRP1 loci provides experimental evidence for this effect since, for a given quantity of sense and antisense RNA, the clustered system produces more siRNA. Although the clustered system also produces read-through transcripts, these would not have a higher hybridization rate than the non-read-through RNAs as the hybridization rate depends on the frequency of collisions between molecules. Intriguingly, Tf2 retrotransposons in Schizosaccharomyces pombe are clustered by the action of centromere protein B (CENP-B), which also silences these elements through histone deacetylation (Cam et al., 2008). This clustering would allow cells to produce siRNA against the Tf2 elements through pervasive transcription, although multiple mechanisms silence Tf2 retrotransposons (Yamanaka et al., 2012), providing an extra defense against retrotransposon activation. Similarly, gypsy retrotransposons in Drosophila are known to cluster (Gerasimova et al., 2000), which may again facilitate siRNA production. Hence, the clustering of transposable elements by factors such as CENP-B would facilitate their recognition by RNAi and allow for selective RNA degradation.

Mammalian germline cells are replete with small RNAs including endogenous siRNA (Watanabe et al., 2006, 2008; Tam et al., 2008; Song et al., 2011), and siRNAs in sperm and oocytes show a pronounced bias towards high-copy sequences that would be effectively explained by the selectivity of the RNAi system towards the products of high-copy loci (Watanabe et al., 2006; Song et al., 2011). However, it remains unclear how some dsRNA precursors of siRNAs are generated, particularly for retrotransposons that are primarily expressed only on the sense strand. We suggest that pervasive transcription would provide sufficient antisense RNA for this role, just as we observed for high-copy GAL cluster sequences in yeast.

In comparison to the germline, the response of mammalian somatic cells to dsRNA is distinctly muted. dsRNA could be processed into siRNA, be altered by RNA editing (Hogg et al., 2011), or could activate the interferon response leading to apoptosis (Gantier and Williams, 2007). However, transgenic mice expressing a hairpin dsRNA construct produce minimal siRNAs, little edited RNA, and show no phenotype that might indicate cell death (Nejepinska et al., 2012a). Nonetheless, siRNAs produced from LINE-1 retrotransposons have been detected in cell culture (Yang and Kazazian, 2006), and high-copy transfected plasmids expressing a sense-antisense RNA pair do produce detectable siRNA in HEK293 cells (Nejepinska et al., 2012b). This shows that a basic siRNA response with an apparent bias towards high-copy sequences is functional in mammalian somatic cells.

Materials and methods

Yeast strains, plasmids and culture conditions

Request a detailed protocol

Yeast deletion strains (Supplementary file 1) were created by standard methods using the oligonucleotides in Supplementary file 1. Plasmids are described in Supplementary file 1 with construction details. Cells were grown on rich media (2% peptone, 1% yeast extract, 2% sugar) or synthetic media (0.69% yeast nitrogen base with ammonium sulfate, amino acids, 2% sugar) for plasmid assays. GAL10 mRNA was induced with 20 µM CuSO4 in Figure 8D. Media components were purchased from Formedium. Cells were grown to mid-log (OD 0.4–0.6) at 30°C for most experiments or at 25°C for experiments involving trf4Δ mutants. The W303 background strain used here has defects in galactose induction in synthetic media, so strains in Figure 3G,H were diploids of W303xBY4741 that show a normal galactose response.

RNA extraction and northern analysis

Request a detailed protocol

RNA was extracted by three procedures described below. High molecular weight RNA was prepared using the hot phenol method for all experiments except Figures 7B, 3B, 3F and Figure 3—figure supplements 1,2 where guanidine thiocyanate (GTC)-phenol preparations were used. 5–10 µg glyoxylated RNA was resolved on 1.2% gels as described (Sambrook and Russell, 2001), transferred to Hybond N+ membrane (GE) and hybridized with probes listed in Supplementary file 1 using either Church Hyb (Sambrook and Russell, 2001) or UltraHyb (Life Technologies). RNA probes were hybridized at 65°C and washed at 65°C using 0.1× SSC, 0.1% SDS, DNA probes in Church Hyb were hybridized at 65°C and washed at 65°C with 0.5× SSC, 0.1% SDS, DNA probes in UltraHyb were hybridized at 42°C and washed at 55°C using 0.2× SSC, 0.1% SDS. Small RNA enriched fractions were isolated using the mirVANA kit (Ambion). 4–10 µg small RNA was separated on 15% polyacrylamide/8 M urea gels containing 20 mM MOPS or 1× TBE, transferred in 20 mM MOPS or 0.5× TBE to Hybond N membrane (GE) and chemically cross-linked as described (Pall and Hamilton, 2008). We observed no difference in cross-linking efficiency between MOPS and TBE gels, but resolution of TBE gels was superior in our hands. siRNAs were detected using random primed probes (Supplementary file 1) hybridized in UltraHyb Oligo (Life Technologies) at 42°C and washed with 2× SSC, 0.5% SDS at 42°C, U6 control oligonucleotide was labeled using T4 polynucleotide kinase and hybridized in Church Hyb under the same conditions.

Hot phenol RNA preparation

Request a detailed protocol

10 × 107 cells in 15 ml tubes were re-suspended in 600 µl TES (10 mM Tris pH 7.5, 10 mM EDTA, 0.5% SDS) and 600 µl phenol pH 7. The mixture was incubated at 65°C for 20 min with 30 s vortexing every 5 min, before briefly chilling on ice. Samples were centrifuged for 5 min and the upper phase extracted. This phase (5–600 µl) was extracted twice with phenol:chloroform (5:1) and once with chloroform before precipitation with 50 µl 3 M sodium acetate (NaOAc) pH 5.2 and 1.1 ml ethanol. The pellet was washed with 70% ethanol and re-suspended in 30 µl water.

GTC-phenol RNA preparation

Request a detailed protocol

2 × 107 cells were lysed by 5 min vortexing at 4°C with 50 µl glass beads and 40 µl GTC-phenol (2.1 M GTC, 26.5 mM Na citrate pH7, 5.3 mM EDTA, 76 mM β-mercaptoethanol, 1.06% N-lauryl sarcosine, 50% phenol pH7). 600 µl GTC-phenol was added, mixed, and samples were heated at 65°C for 10 min, then placed on ice for 10 min. 160 µl 100 mM NaOAc pH 5.2 and 300 µl chloroform:isoamyl alcohol (24:1) were added, samples were vortexed and centrifuged at top speed for 5 min. The upper phase was extracted, precipitated with 1 ml ethanol, washed with 70% ethanol and the pellet re-suspended in 6 µl water. 3 µl RNA was analyzed per lane.

Small RNA purification

Request a detailed protocol

Small RNAs were isolated using a mirVANA kit (Life Technologies) with minor modifications. 35 × 107 cells were thoroughly re-suspended in 100 µl lysis/binding buffer, 200 µl glass beads were added, and the samples were vortexed for 5 min at 4°C. 500 µl lysis/binding buffer were added and the samples were mixed before proceeding with the isolation as per the manufacturer’s instructions, starting with addition of the miRNA homogenate additive. After isolation the samples were generally re-precipitated and re-suspended in 20 µl water.

RNase A treatment

Request a detailed protocol

20 × 107 cells were harvested and split into two aliquots, then re-suspended in 600 µl 10 mM Tris pH 7.5, 10 mM EDTA on ice. Cells were lysed with glass beads (10 cycles of 30 s vortex, 60 s on ice), and 5 µg RNase A was added to one aliquot followed by 30 min incubation on ice. After centrifugation for 10 min at 4500×g, 600 µl lysate was extracted, SDS added to 0.5%, and RNA extracted by the hot phenol method as above. RNase A treated samples were re-suspended in 12 µl water.

DNA isolation and Southern blotting

Request a detailed protocol

Cells from 2 ml saturated culture were washed with 50 mM EDTA, then spheroplasted with 250 µl 0.34 U/ml lyticase (Sigma L4025) in 1.2 M sorbitol, 50 mM EDTA, 10 mM DTT. After centrifuging at 1000×g, the cells were gently resuspended in 400 µl of 0.3% SDS, 50 mM EDTA, 100 µg/ml RNase A and incubated at 37°C for 30 min. 4 µl of 20 mg/ml proteinase K was added and the samples were mixed by inversion and heated to 65°C for 30 min. 160 µl 5 M potassium acetate was added after cooling to room temperature, the samples were mixed by inversion and then chilled on ice for 30 min. After 10 min centrifuging at 10,000×g, the supernatant was poured into a new tube containing 500 µl phenol:chloroform pH 8 and the samples were mixed on a wheel for 15 min. The samples were centrifuged for 10 min at 10,000×g and the upper phase was extracted using cut tips and precipitated with 400 µl isopropanol. Pellets were washed with 70% ethanol, air-dried and left overnight at 4°C to dissolve in 30 µl TE. After gentle mixing, 10 µl of each sample was digested with 40 U of EcoRV, ethanol precipitated and separated on a 25 cm 1% TBE gel at 90 V overnight. The gel was washed in 0.25 N HCl for 15 min, 0.5 N NaOH for 45 min, and twice in 1.5 M NaCl 0.5 M Tris pH 7.5 for 20 min before transfer to HyBond N+ membrane in 6× SSC. The membrane was probed for URA3 and GAL7 in Church Hyb at 65°C and washed with 0.5× SSC, 0.1% SDS at 65°C.


Request a detailed protocol

Cells were grown to OD 0.5 in YPD and the cultures mixed as required, 4 ml per coverslip. The cells were fixed with 440 µl 37% formaldehyde (Merck, microscopy grade) for 15 min at room temperature, then centrifuged for 2 min at 4600 rpm. The cells were washed three times with 1 ml of buffer B (0.1 M sodium phosphate pH 7.5, 1.2 M sorbitol), then re-suspended in 100 µl buffer B containing 3 µl 17 U/µl lyticase (Sigma L2524) and 10 mM DTT for 15 min. The cells were centrifuged for 2 min at 1000×g, then washed with 1 ml buffer B. The cells were re-suspended in 40 µl buffer B, applied to a poly-L-lysine coated coverslip (Zeiss 18 × 18 × 0.170 ± 0.005 mm) and left for 20 min before washing twice with buffer B. Coverslips were treated with −20°C methanol for 6 min, then dipped in −20°C acetone for 10 s, followed by two washes with PBS. Coverslips were blocked for 30 min with 5% milk 0.3% Triton-X100 in PBS, washed with PBS, then incubated overnight at 4°C with primary antibodies in 50 µl 1% BSA 0.3% Triton-X100 in PBS. Coverslips were washed three times with PBS and incubated for 30 min at room temperature with secondary antibodies 1:1000 in same buffer as primaries. After washing three times with PBS, the samples were dehydrated with 70%, 90%, and 100% ethanol and mounted in Pro-long Gold with DAPI (Life Technologies). Antibodies were rabbit anti-GFP (Life Technologies A11122) at 1:500 and mouse anti-dsRNA J2 (ESC 10010200) at 1:1000. Images were acquired using a Nikon N-SIM microscope comprising a Nikon Ti-E microscope, Nikon 100× 1.49 NA lens, Nikon SIM illumination module, and Andor iXon 897 EM-CCD camera. SIM data were acquired in ‘3D-SIM’ mode using five phases and three rotations. DAPI, Alexa Fluor 488, and Alexa Fluor 594 dyes were excited using 405, 488, and 561 nm laser light, respectively. Super-resolution images were reconstructed using Nikon Elements software. Equivalent wide-field images were reconstructed in FIJI (ImageJ, NIH) by summing the phase shifts from one grid rotation.

Imaging and analysis

Request a detailed protocol

Gels and phosphor screens were imaged using FLA 3000 (Fuji) or FLA 7000 (GE) systems. Quantification was performed using AIDA (Fuji) or ImageQuant (GE), depending on the scanner used. Images were prepared for publication with ImageJ by smoothing and minimal contrast enhancement. Images from the FLA3000 had a Gamma Shift of 3 applied.


Request a detailed protocol

5′ RACE was performed with an ExactSTART Eukaryotic mRNA 5′- and 3′-RACE Kit (Epicentre) as per manufacturer’s instructions, except that reverse transcription was primed from random hexamers.

RNA immunoprecipitation

Request a detailed protocol

50 ml of cells at 0.7 × 107 cells/ml were harvested, washed, and frozen on nitrogen. Cells were thawed, washed twice in 1 ml lysis buffer (50 mM HEPES pH7.5, 50 mM KCl, 5 mM MgCl2, 1 mM DTT, 1× complete protease inhibitors (Roche)) and transferred to 2 ml tubes, then re-suspended in 300 µl lysis buffer. 300 µl zirconium beads were added and cells were lysed by vortexing 10 × 30 s with 30 s on ice between cycles. The lysate was clarified by centrifuging twice for 5 min at 14,000 rpm, a 12 µl aliquot was removed for total RNA and the remaining lysate was split in half and 2.5 µl mouse anti-dsRNA J2 (ESC 10010200) added to one aliquot. Antibody was bound for 2 hr at 4°C, then 20 µl GammaBind beads (GE) pre-blocked overnight with 1% BSA were added and incubation continued on a wheel for 2 hr at 4°C. The beads were washed 5× for 10 min with 1 ml wash buffer (10 mM Tris pH 7.5, 120 mM NaCl, 5 mM MgCl2, 0.1% NP-40, 1 mM DTT). RNA was extracted from the beads and total lysate samples using Tri-reagent (Sigma) according to the manufacturer’s instructions. 1 µg total lysate and whole immunoprecipitation samples were treated with RQ1 DNase (Promega), purified by phenol:chloroform and ethanol precipitation, then reverse transcribed from random hexamers using Superscript III (Life Technologies). Quantitative PCR was performed using Maxima SYBR qPCR mix (Fermentas).


Antisense analysis

Request a detailed protocol

Locations of XUTs (van Dijk et al., 2011) were merged with CUT and SUT locations (Xu et al., 2009), along with expression validated ORFs (Xu et al., 2009), and overlap between ORFs and other features calculated using an R script.

Table 2

n values for statistical tests

4B0–2 single-copy19,8920–2 multi-copy1128
2–4 single-copy36,4742–4 multi-copy1111
4–6 single-copy34,9214–6 multi-copy1177
6–8 single-copy16,2086–8 multi-copy1356
8–10 single-copy40608–10 multi-copy1424
10–12 single-copy92210–12 multi-copy1505
12–14 single-copy27012–14 multi-copy887
14–16 single-copy5614–16 multi-copy199
4CBulk 1–2 copies85,495Top 1% 1–2 copies252
Bulk 2–4 copies1419Top 1% 2–4 copies95
Bulk 4–8 copies544Top 1% 4–8 copies98
Bulk 8–16 copies981Top 1% 8–16 copies238
Bulk 16–32 copies1978Top 1% 16–32 copies212
Bulk >32 copies695Top 1% >32 copies25
4—Supplement 1A0–2 single-copy19,8920–2 multi-copy836
2–4 single-copy36,4742–4 multi-copy906
4–6 single-copy34,9214–6 multi-copy724
6–8 single-copy16,2086–8 multi-copy387
8–10 single-copy40608–10 multi-copy189
10–12 single-copy92210–12 multi-copy167
12–14 single-copy27012–14 multi-copy128
14–16 single-copy5614–16 multi-copy29
4—Supplement 1B0–2 low-copy20,7280–2 high-copy292
2–4 low-copy37,3802–4 high-copy205
4–6 low-copy35,6454–6 high-copy453
6–8 low-copy16,5956–8 high-copy969
8–10 low-copy42498–10 high-copy1235
10–12 low-copy108910–12 high-copy1338
12–14 low-copy39812–14 high-copy759
14–16 low-copy8514–16 high-copy170
4-Supplement 2A0–2 single-copy20,0300–2 multi-copy1999
2–4 single-copy36,4942–4 multi-copy2104
4–6 single-copy34,8834–6 multi-copy2098
6–8 single-copy16,1786–8 multi-copy1668
8–10 single-copy40788–10 multi-copy644
10–12 single-copy87310–12 multi-copy204
12–14 single-copy21512–14 multi-copy66
14–16 single-copy5214–16 multi-copy10
4-Supplement 2B0–2 single-copy19,8920–2 multi-copy1128
2–4 single-copy36,4742–4 multi-copy1111
4–6 single-copy34,9214–6 multi-copy1177
6–8 single-copy16,2086–8 multi-copy1356
8–10 single-copy40608–10 multi-copy1424
10–12 single-copy92210–12 multi-copy1505
12–14 single-copy27012–14 multi-copy887
14–16 single-copy5614–16 multi-copy199

siRNA analysis

Request a detailed protocol

Sequencing data for mRNA and siRNA fractions in RNAi strains (GEO accession GSE31300) (Drinnenberg et al., 2011) and W303 total RNA (GEO accession GSE38383) (Hobson et al., 2012) were quality and adapter trimmed using Trim Galore (v0.2.3; default options; http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) and aligned to the yeast genome (build SGD1.01) using Bowtie (Langmead et al., 2009) (v0.12.7; default settings plus ‘--best’) allowing non-unique sequences to be assigned at random. For expressed gene analysis (Figure 1), reads overlapping each ORF were binned and only ORFs with >100 reads were used. For siRNA profiles (Figure 2—figure supplement 1), reads were binned over 50 bp intervals using SeqMonk (http://www.bioinformatics.babraham.ac.uk/projects/seqmonk/). For siRNA versus total RNA expression (Figure 4), read counts were quantified in consecutive 100 bp bins across the genome using SeqMonk, bins with >10,0000 total RNA reads were excluded as were bins derived from 2µ sequence which is single copy in the genome sequence but high copy in reality, and a pseudocount of one read was added to total and siRNA read counts for each bin. Total RNA levels were multiplied by copy number to correct for the division of reads amongst copies that occurs during mapping, or alternate normalization was applied in Figure 4—figure supplement 2 (see below for the reasoning underlying this copy number normalization methodology). The copy number for each bin was determined by splitting the complete genomic sequence into overlapping 20 bp segments at 1 bp intervals and re-mapping to the genome with reads allowed to align to all perfectly matching sequences, producing a measure of the number of genomic sequences matching each 100 bp bin. Little difference was seen if one mismatch was allowed (data not shown).

Read count normalization for multi-copy sequences

Request a detailed protocol

Multi-copy loci are problematic for standard high-throughput sequencing mapping pipelines and are commonly discarded. Reads mapping to a non-unique genome sequence are usually assigned at random to one copy in the genome, therefore the total reads are divided evenly amongst the copies and the apparent abundance of RNA matching each copy is effectively divided by the copy number. In order to assess total RNA abundance (as in Figure 4A,B), we multiplied all total RNA read counts by the copy number to obtain the real total RNA abundance. However, we decided that the siRNA abundance should be analyzed per producing locus because each copy in the genome was analyzed separately for comparison with single-copy loci (i.e., we quantified how many siRNA reads an individual copy of a multi-copy locus produced, not how many the combined copies produced). We therefore did not multiply the siRNA read counts by the copy number. Such normalizations clearly have the potential to introduce systematic biases, and we therefore repeated the analysis in Figure 4 either with no copy number normalization or with both total and siRNA read counts multiplied by copy number (Figure 4—figure supplement 2). We found that although the distributions changed somewhat, the majority of total RNA abundance categories had higher siRNA levels for multi-copy loci irrespective of the copy number normalization applied. We note that the alternative copy number analysis in Figure 4C,D did not require any normalization; indeed, any copy number normalization would simply cancel out in the calculation of the siRNA:total RNA ratio, therefore any systematic bias that might be introduced by copy number normalization would have no effect.


  1. 1
  2. 2
  3. 3
  4. 4
  5. 5
    Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project
    1. E Birney
    2. JA Stamatoyannopoulos
    3. A Dutta
    4. R Guigo
    5. TR Gingeras
    6. EH Margulies
    7. Z Weng
    8. M Snyder
    9. ET Dermitzakis
    10. RE Thurman
    11. MS Kuehn
    12. CM Taylor
    13. S Neph
    14. CM Koch
    15. S Asthana
    16. A Malhotra
    17. I Adzhubei
    18. JA Greenbaum
    19. RM Andrews
    20. P Flicek
    21. PJ Boyle
    22. H Cao
    23. NP Carter
    24. GK Clelland
    25. S Davis
    26. N Day
    27. P Dhami
    28. SC Dillon
    29. MO Dorschner
    30. H Fiegler
    31. PG Giresi
    32. J Goldy
    33. M Hawrylycz
    34. A Haydock
    35. R Humbert
    36. KD James
    37. BE Johnson
    38. EM Johnson
    39. TT Frum
    40. ER Rosenzweig
    41. N Karnani
    42. K Lee
    43. GC Lefebvre
    44. PA Navas
    45. F Neri
    46. SC Parker
    47. PJ Sabo
    48. R Sandstrom
    49. A Shafer
    50. D Vetrie
    51. M Weaver
    52. S Wilcox
    53. M Yu
    54. FS Collins
    55. J Dekker
    56. JD Lieb
    57. TD Tullius
    58. GE Crawford
    59. S Sunyaev
    60. WS Noble
    61. I Dunham
    62. F Denoeud
    63. A Reymond
    64. P Kapranov
    65. J Rozowsky
    66. D Zheng
    67. R Castelo
    68. A Frankish
    69. J Harrow
    70. S Ghosh
    71. A Sandelin
    72. IL Hofacker
    73. R Baertsch
    74. D Keefe
    75. S Dike
    76. J Cheng
    77. HA Hirsch
    78. EA Sekinger
    79. J Lagarde
    80. JF Abril
    81. A Shahab
    82. C Flamm
    83. C Fried
    84. J Hackermüller
    85. J Hertel
    86. M Lindemeyer
    87. K Missal
    88. A Tanzer
    89. S Washietl
    90. J Korbel
    91. O Emanuelsson
    92. JS Pedersen
    93. N Holroyd
    94. R Taylor
    95. D Swarbreck
    96. N Matthews
    97. MC Dickson
    98. DJ Thomas
    99. MT Weirauch
    100. J Gilbert
    101. J Drenkow
    102. I Bell
    103. X Zhao
    104. KG Srinivasan
    105. WK Sung
    106. HS Ooi
    107. KP Chiu
    108. S Foissac
    109. T Alioto
    110. M Brent
    111. L Pachter
    112. ML Tress
    113. A Valencia
    114. SW Choo
    115. CY Choo
    116. C Ucla
    117. C Manzano
    118. C Wyss
    119. E Cheung
    120. TG Clark
    121. JB Brown
    122. M Ganesh
    123. S Patel
    124. H Tammana
    125. J Chrast
    126. CN Henrichsen
    127. C Kai
    128. J Kawai
    129. U Nagalakshmi
    130. J Wu
    131. Z Lian
    132. J Lian
    133. P Newburger
    134. X Zhang
    135. P Bickel
    136. JS Mattick
    137. P Carninci
    138. Y Hayashizaki
    139. S Weissman
    140. T Hubbard
    141. RM Myers
    142. J Rogers
    143. PF Stadler
    144. TM Lowe
    145. CL Wei
    146. Y Ruan
    147. K Struhl
    148. M Gerstein
    149. SE Antonarakis
    150. Y Fu
    151. ED Green
    152. U Karaöz
    153. A Siepel
    154. J Taylor
    155. LA Liefer
    156. KA Wetterstrand
    157. PJ Good
    158. EA Feingold
    159. MS Guyer
    160. GM Cooper
    161. G Asimenos
    162. CN Dewey
    163. M Hou
    164. S Nikolaev
    165. JI Montoya-Burgos
    166. A Löytynoja
    167. S Whelan
    168. F Pardi
    169. T Massingham
    170. H Huang
    171. NR Zhang
    172. I Holmes
    173. JC Mullikin
    174. A Ureta-Vidal
    175. B Paten
    176. M Seringhaus
    177. D Church
    178. K Rosenbloom
    179. WJ Kent
    180. EA Stone
    181. S Batzoglou
    182. N Goldman
    183. RC Hardison
    184. D Haussler
    185. W Miller
    186. A Sidow
    187. ND Trinklein
    188. ZD Zhang
    189. L Barrera
    190. R Stuart
    191. DC King
    192. A Ameur
    193. S Enroth
    194. MC Bieda
    195. J Kim
    196. AA Bhinge
    197. N Jiang
    198. J Liu
    199. F Yao
    200. VB Vega
    201. CW Lee
    202. P Ng
    203. A Shahab
    204. A Yang
    205. Z Moqtaderi
    206. Z Zhu
    207. X Xu
    208. S Squazzo
    209. MJ Oberley
    210. D Inman
    211. MA Singer
    212. TA Richmond
    213. KJ Munn
    214. A Rada-Iglesias
    215. O Wallerman
    216. J Komorowski
    217. JC Fowler
    218. P Couttet
    219. AW Bruce
    220. OM Dovey
    221. PD Ellis
    222. CF Langford
    223. DA Nix
    224. G Euskirchen
    225. S Hartman
    226. AE Urban
    227. P Kraus
    228. S Van Calcar
    229. N Heintzman
    230. TH Kim
    231. K Wang
    232. C Qu
    233. G Hon
    234. R Luna
    235. CK Glass
    236. MG Rosenfeld
    237. SF Aldred
    238. SJ Cooper
    239. A Halees
    240. JM Lin
    241. HP Shulha
    242. X Zhang
    243. M Xu
    244. JN Haidar
    245. Y Yu
    246. Y Ruan
    247. VR Iyer
    248. RD Green
    249. C Wadelius
    250. PJ Farnham
    251. B Ren
    252. RA Harte
    253. AS Hinrichs
    254. H Trumbower
    255. H Clawson
    256. J Hillman-Jackson
    257. AS Zweig
    258. K Smith
    259. A Thakkapallayil
    260. G Barber
    261. RM Kuhn
    262. D Karolchik
    263. L Armengol
    264. CP Bird
    265. PI de Bakker
    266. AD Kern
    267. N Lopez-Bigas
    268. JD Martin
    269. BE Stranger
    270. A Woodroffe
    271. E Davydov
    272. A Dimas
    273. E Eyras
    274. IB Hallgrímsdóttir
    275. J Huppert
    276. MC Zody
    277. GR Abecasis
    278. X Estivill
    279. GG Bouffard
    280. X Guan
    281. NF Hansen
    282. JR Idol
    283. VV Maduro
    284. B Maskeri
    285. JC McDowell
    286. M Park
    287. PJ Thomas
    288. AC Young
    289. RW Blakesley
    290. DM Muzny
    291. E Sodergren
    292. DA Wheeler
    293. KC Worley
    294. H Jiang
    295. GM Weinstock
    296. RA Gibbs
    297. T Graves
    298. R Fulton
    299. ER Mardis
    300. RK Wilson
    301. M Clamp
    302. J Cuff
    303. S Gnerre
    304. DB Jaffe
    305. JL Chang
    306. K Lindblad-Toh
    307. ES Lander
    308. M Koriabine
    309. M Nefedov
    310. K Osoegawa
    311. Y Yoshinaga
    312. B Zhu
    313. PJ de Jong
    Nature 447:799–816.
  6. 6
  7. 7
  8. 8
  9. 9
  10. 10
  11. 11
    The transcriptional landscape of the mammalian genome
    1. P Carninci
    2. T Kasukawa
    3. S Katayama
    4. J Gough
    5. MC Frith
    6. N Maeda
    7. R Oyama
    8. T Ravasi
    9. B Lenhard
    10. C Wells
    11. R Kodzius
    12. K Shimokawa
    13. VB Bajic
    14. SE Brenner
    15. S Batalov
    16. AR Forrest
    17. M Zavolan
    18. MJ Davis
    19. LG Wilming
    20. V Aidinis
    21. JE Allen
    22. A Ambesi-Impiombato
    23. R Apweiler
    24. RN Aturaliya
    25. TL Bailey
    26. M Bansal
    27. L Baxter
    28. KW Beisel
    29. T Bersano
    30. H Bono
    31. AM Chalk
    32. KP Chiu
    33. V Choudhary
    34. A Christoffels
    35. DR Clutterbuck
    36. ML Crowe
    37. E Dalla
    38. BP Dalrymple
    39. B de Bono
    40. G Della Gatta
    41. D di Bernardo
    42. T Down
    43. P Engstrom
    44. M Fagiolini
    45. G Faulkner
    46. CF Fletcher
    47. T Fukushima
    48. M Furuno
    49. S Futaki
    50. M Gariboldi
    51. P Georgii-Hemming
    52. TR Gingeras
    53. T Gojobori
    54. RE Green
    55. S Gustincich
    56. M Harbers
    57. Y Hayashi
    58. TK Hensch
    59. N Hirokawa
    60. D Hill
    61. L Huminiecki
    62. M Iacono
    63. K Ikeo
    64. A Iwama
    65. T Ishikawa
    66. M Jakt
    67. A Kanapin
    68. M Katoh
    69. Y Kawasawa
    70. J Kelso
    71. H Kitamura
    72. H Kitano
    73. G Kollias
    74. SP Krishnan
    75. A Kruger
    76. SK Kummerfeld
    77. IV Kurochkin
    78. LF Lareau
    79. D Lazarevic
    80. L Lipovich
    81. J Liu
    82. S Liuni
    83. S McWilliam
    84. M Madan Babu
    85. M Madera
    86. L Marchionni
    87. H Matsuda
    88. S Matsuzawa
    89. H Miki
    90. F Mignone
    91. S Miyake
    92. K Morris
    93. S Mottagui-Tabar
    94. N Mulder
    95. N Nakano
    96. H Nakauchi
    97. P Ng
    98. R Nilsson
    99. S Nishiguchi
    100. S Nishikawa
    101. F Nori
    102. O Ohara
    103. Y Okazaki
    104. V Orlando
    105. KC Pang
    106. WJ Pavan
    107. G Pavesi
    108. G Pesole
    109. N Petrovsky
    110. S Piazza
    111. J Reed
    112. JF Reid
    113. BZ Ring
    114. M Ringwald
    115. B Rost
    116. Y Ruan
    117. SL Salzberg
    118. A Sandelin
    119. C Schneider
    120. C Schönbach
    121. K Sekiguchi
    122. CA Semple
    123. S Seno
    124. L Sessa
    125. Y Sheng
    126. Y Shibata
    127. H Shimada
    128. K Shimada
    129. D Silva
    130. B Sinclair
    131. S Sperling
    132. E Stupka
    133. K Sugiura
    134. R Sultana
    135. Y Takenaka
    136. K Taki
    137. K Tammoja
    138. SL Tan
    139. S Tang
    140. MS Taylor
    141. J Tegner
    142. SA Teichmann
    143. HR Ueda
    144. E van Nimwegen
    145. R Verardo
    146. CL Wei
    147. K Yagi
    148. H Yamanishi
    149. E Zabarovsky
    150. S Zhu
    151. A Zimmer
    152. W Hide
    153. C Bult
    154. SM Grimmond
    155. RD Teasdale
    156. ET Liu
    157. V Brusic
    158. J Quackenbush
    159. C Wahlestedt
    160. JS Mattick
    161. DA Hume
    162. C Kai
    163. D Sasaki
    164. Y Tomaru
    165. S Fukuda
    166. M Kanamori-Katayama
    167. M Suzuki
    168. J Aoki
    169. T Arakawa
    170. J Iida
    171. K Imamura
    172. M Itoh
    173. T Kato
    174. H Kawaji
    175. N Kawagashira
    176. T Kawashima
    177. M Kojima
    178. S Kondo
    179. H Konno
    180. K Nakano
    181. N Ninomiya
    182. T Nishio
    183. M Okada
    184. C Plessy
    185. K Shibata
    186. T Shiraki
    187. S Suzuki
    188. M Tagami
    189. K Waki
    190. A Watahiki
    191. Y Okamura-Oho
    192. H Suzuki
    193. J Kawai
    194. Y Hayashizaki
    195. FANTOM Consortium
    196. Genome Network Project Core Group
    Science 309:1559–1563.
  12. 12
  13. 13
  14. 14
  15. 15
  16. 16
  17. 17
  18. 18
  19. 19
  20. 20
  21. 21
  22. 22
  23. 23
  24. 24
  25. 25
  26. 26
  27. 27
  28. 28
  29. 29
  30. 30
  31. 31
  32. 32
  33. 33
  34. 34
  35. 35
  36. 36
  37. 37
  38. 38
  39. 39
  40. 40
  41. 41
  42. 42
  43. 43
  44. 44
  45. 45
  46. 46
  47. 47
  48. 48
  49. 49
  50. 50
  51. 51
  52. 52
  53. 53
  54. 54
  55. 55
  56. 56
  57. 57
  58. 58
  59. 59
  60. 60
  61. 61
  62. 62
  63. 63
  64. 64
  65. 65
  66. 66
  67. 67
  68. 68
  69. 69
    Molecular cloning: a laboratory manual (3rd edition)
    1. J Sambrook
    2. DW Russell
    Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press.
  70. 70
  71. 71
  72. 72
  73. 73
  74. 74
    Male germ cells express abundant endogenous siRNAs
    1. R Song
    2. GW Hennig
    3. Q Wu
    4. C Jose
    5. H Zheng
    6. W Yan
    Proceedings of the National Academy of Sciences of the United States of America 108:13159–13164.
  75. 75
  76. 76
  77. 77
    A role for noncoding transcription in activation of the yeast PHO5 gene
    1. JP Uhler
    2. C Hertel
    3. JQ Svejstrup
    Proceedings of the National Academy of Sciences of the United States of America 104:8011–8016.
  78. 78
  79. 79
  80. 80
  81. 81
  82. 82
  83. 83
  84. 84
  85. 85
  86. 86
  87. 87
  88. 88
  89. 89

Decision letter

  1. Nick J Proudfoot
    Reviewing Editor; University of Oxford, United Kingdom

eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.

Thank you for sending your work entitled “Endogenous RNA interference is driven by copy number” for consideration at eLife. Your article has been favorably evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors.

The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.

This study employs the so-called RNAi+ S. cerevisiae strain generated in the Bartel lab (Drinnenberg et al 2009) that is transformed with S. Castellii Dicer and Ago genes. This original study suggested that this RNAi+ yeast would be a good model system to explore RNAi gene silencing mechanisms. It is clear that this is the case as this new study makes excellent use of this strain. What is shown here is that endogenous convergent transcription (especially mRNA sense and ncRNA antisense arrangement) generates siRNAs in RNAi+ cells but this does not correlate with reduced mRNA abundance (Figure 1). Also repetitive sequences (Y′ RNA and rRNA) generate abundant siRNAs and show significant RNAi dependent RNA degradation. Some of these data are however predicted from the original Drinnenberg et al study. New results are presented in Figure 3 showing in two test gene cases (MAL32 and GAL4) that forced sense antisense expression of these genes on multicopy transformed plasmids generate both siRNAs and also show clear gene silencing. These results lead to the key observation that it isn't the overall level of expression (sense antisense) that is critical but rather it is the fact that expression comes from multicopy genes (rather than one highly expressed gene). This point is initially made from bioinformatic analysis of the RNAi+ transcriptiome (Drinneneberg et al) where multicopy genes generate more siRNA than single copy rather than simply correlating with totals RNA levels (Figure 4). Figures 5 and 6 then validate the bioinformatics very clearly. The integrated single copy MAL32 expressing high level sense and antisense transcription due to the use of strong convergent arranged promoters does not generate as much siRNA as multicopy low expressed MAL32. Also sense antisense ncRNAs expressed on multicopy plasmids generate significant amounts of siRNA. Finally some evidence was provided that siRNAs from high-copy GAL10 constructs (making sense and antisense RNAs) can generate siRNA that silence another GAL10 gene in trans.

Overall the three reviewers of this paper while aware of the somewhat artificial nature of the RNAi+ S. cerevisiae strain consider that the experimental benefit of this system has been put to excellent use and has allowed significant new insight in RNAi mechanism to be obtained. The key additional experiments that we feel are necessary to make a revised manuscript acceptable for publication by eLife are as follows.

1) The actual copy number of the 2µ plasmids transformed into the RNAi+ strain should be directly measured by Southern blotting

2) The Castellii Dcr1 gene should be knocked out to show loss of the high-copy number silencing observed for MAL32 and GAL4.

3) The location of dsRNA and Dicer in the transformed RNAi+ cells should be directly visualised (J2 antibody for dsRNA and Tagged Dcr1/anti Dcr1) to show either chromatin or cytolasmic association. See Schönborn, J. et al. NAR 19, 2993-3000, 1991; Gullerova and Proudfoot NSMB 19, 1193-1201 2013 for J2 antibody.

4) Better quantitation of RNAi effect on mRNA levels and level of dsRNA is required. This can be achieved by RTqPCR and for dsRNA level of J2 dsRNA immunoprecipitation.


Author response

1) The actual copy number of the 2µ plasmids transformed into the RNAi+ strain should be directly measured by Southern blotting

We have included a southern for 2µ copy number as Figure 3–figure supplement 5. We have also shown a northern blot for the MAL32 high-copy plasmid in which the plasmid signals are clearly visible. Both analyses show that the copy number does not decrease in the RNAi+ strain – it actually increases, an observation we discuss in the text.

2) The Castellii Dcr1 gene should be knocked out to show loss of the high-copy number silencing observed for MAL32 and GAL4.

We decided to combine this with changing the genetic background of the RNAi+ strain. We cloned Dcr1 and Ago1 onto single copy plasmids and introduced these into BY4741 cells (the original background is W303), along with the high-copy plasmids. The result for MAL32 is very clear, showing that both genes are required for knockdown (Figure 3–figure supplement 3). However, in this genetic background GAL4 is only expressed very weakly from the high-copy plasmid; the RNA forms dsRNA efficiently as expected and this is removed by Dicer, although this does not lead to a clear knockdown. The GAL4 knockdown was the weakest we observed in the original experiments, and RNAi knockdowns in this new strain are not as strong as in the original strain (60% in new strain vs. 80% in the original strain for the MAL32 mRNA). Combined with the confounding effect of the increased plasmid copy number in RNAi+ cells, the GAL4 knockdown becomes too weak to test the importance of Dicer. We do not think that this compromises the conclusions of the paper.

3) The location of dsRNA and Dicer in the transformed RNAi+ cells should be directly visualised (J2 antibody for dsRNA and Tagged Dcr1/anti Dcr1) to show either chromatin or cytolasmic association. See Schönborn, J. et al. NAR 19, 2993-3000, 1991; Gullerova and Proudfoot NSMB 19, 1193-1201 2013 for J2 antibody.

We have included this data as Figure 6–figure supplement 1. We could not convincingly detect Dicer in nuclei, and the dsRNA signal is dominated by strong cytoplasmic signals, most likely from Killer virus RNAs. Although RNAi should technically clear Killer virus, the analysis performed by the Bartel lab shows that this clearance is far from complete. We have used a different approach to improve confidence in nuclear formation of dsRNA; if dsRNA forms in the nucleus then the juxtaposition of transcribing copies should be important, and we demonstrate that this is true for the TRP1 gene in the new Figure 6. We have referred to both of these experiments in the Discussion.

4) Better quantitation of RNAi effect on mRNA levels and level of dsRNA is required. This can be achieved by RTqPCR and for dsRNA level of J2 dsRNA immunoprecipitation.

We have added a number of additional results here. Firstly, we have added quantification for the effect of RNAi on endogenous MAL32 and GAL4 mRNA along with full length northern blot images that show degradation patterns (Figure 3–figure supplement 2). We have also performed the J2 immunoprecipitations (Figure 3–figure supplement 6) – these show a clear loss dsRNA in the RNAi strain, and confirm the Northern blot results for mRNA knockdown.


Article and author information

Author details

  1. Cristina Cruz

    Epigenetics Programme, The Babraham Institute, Cambridge, United Kingdom
    CC, Conception and design, Acquisition of data, Analysis and interpretation of data
    Competing interests
    The authors declare that no competing interests exist.
  2. Jonathan Houseley

    Epigenetics Programme, The Babraham Institute, Cambridge, United Kingdom
    JH, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    Competing interests
    The authors declare that no competing interests exist.


Wellcome Trust (088335)

  • Cristina Cruz
  • Jonathan Houseley

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.


We thank Simon Andrews, Felix Krueger, Laura Biggins and Anne Segonds-Pichon for Bioinformatics support, Simon Walker for microscopy support, David Bartel for strains, Robin Allshire, Wolf Reik and Sarah Elderkin for comments, and Alex Murray and Tim Hore for critical reading. This work was supported by the Wellcome Trust [grant number 088335].

Reviewing Editor

  1. Nick J Proudfoot, University of Oxford, United Kingdom

Publication history

  1. Received: September 23, 2013
  2. Accepted: December 29, 2013
  3. Version of Record published: February 11, 2014 (version 1)


© 2014, Cruz and Houseley

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.


  • 1,753
    Page views
  • 178
  • 15

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Download citations (links to download the citations from this article in formats compatible with various reference manager tools)

Open citations (links to open the citations from this article in various online reference manager services)

Further reading

    1. Cell Biology
    2. Structural Biology and Molecular Biophysics
    Salah A Baker et al.
    Research Article Updated

    Interstitial cells of Cajal (ICC) generate pacemaker activity responsible for phasic contractions in colonic segmentation and peristalsis. ICC along the submucosal border (ICC-SM) contribute to mixing and more complex patterns of colonic motility. We show the complex patterns of Ca2+ signaling in ICC-SM and the relationship between ICC-SM Ca2+ transients and activation of smooth muscle cells (SMCs) using optogenetic tools. ICC-SM displayed rhythmic firing of Ca2+transients ~ 15 cpm and paced adjacent SMCs. The majority of spontaneous activity occurred in regular Ca2+ transients clusters (CTCs) that propagated through the network. CTCs were organized and dependent upon Ca2+ entry through voltage-dependent Ca2+ conductances, L- and T-type Ca2+ channels. Removal of Ca2+ from the external solution abolished CTCs. Ca2+ release mechanisms reduced the duration and amplitude of Ca2+ transients but did not block CTCs. These data reveal how colonic pacemaker ICC-SM exhibit complex Ca2+-firing patterns and drive smooth muscle activity and overall colonic contractions.