Genomic Analysis Suggests that mRNA Destabilization by the Microprocessor Is Specialized for the Auto-Regulation of Dgcr8

Abstract

Background

The Microprocessor, containing the RNA binding protein Dgcr8 and RNase III enzyme Drosha, is responsible for processing primary microRNAs to precursor microRNAs. The Microprocessor regulates its own levels by cleaving hairpins in the 5′UTR and coding region of the Dgcr8 mRNA, thereby destabilizing the mature transcript.

Methodology/Principal Findings

To determine whether the Microprocessor has a broader role in directly regulating other coding mRNA levels, we integrated results from expression profiling and ultra high-throughput deep sequencing of small RNAs. Expression analysis of mRNAs in wild-type, Dgcr8 knockout, and Dicer knockout mouse embryonic stem (ES) cells uncovered mRNAs that were specifically upregulated in the Dgcr8 null background. A number of these transcripts had evolutionarily conserved predicted hairpin targets for the Microprocessor. However, analysis of deep sequencing data of 18 to 200nt small RNAs in mouse ES, HeLa, and HepG2 indicates that exonic sequence reads that map in a pattern consistent with Microprocessor activity are unique to Dgcr8.

Conclusion/Significance

We conclude that the Microprocessor's role in directly destabilizing coding mRNAs is likely specifically targeted to Dgcr8 itself, suggesting a specialized cellular mechanism for gene auto-regulation.

About the Authors:

Archana Shenoy

Affiliation: The Eli and Edythe Broad Center of Regeneration Medicine and Stem Cell Research, Center for Reproductive Sciences, and Department of Urology, University of California San Francisco, San Francisco, California, United States of America

Robert Blelloch

* E-mail: [email protected]

Affiliation: The Eli and Edythe Broad Center of Regeneration Medicine and Stem Cell Research, Center for Reproductive Sciences, and Department of Urology, University of California San Francisco, San Francisco, California, United States of America

Introduction

MicroRNA maturation involves two processing steps [1]. First, a long primary miRNA (pri-miRNA) is cleaved by the Microprocessor, containing the RNA binding protein Dgcr8 and the RNAseIII enzyme Drosha, to produce a 60–75 nucleotide hairpin precursor miRNA (pre-miRNA) in the nucleus [2], [3], [4], [5], [6]. The pre-miRNA is translocated to the cytoplasm where it is cleaved to a miRNA duplex (∼19–25 nt in length) by the RNAseIII enzyme Dicer [7]. A single strand of the duplex enters the RNA induced silencing complex (RISC) with the help of another RNA binding protein, TRBP [8], [9]. Dicer has roles outside of the maturation of canonical miRNAs. For example, in mouse ES cells, Dicer processes other subclasses of miRNAs including mirtrons and short hairpin RNAs as well as endogenous siRNAs [10]. Similarly, Dicer processes endogenous siRNAs in mouse oocytes [11], [12]. Consistent with these additional roles of Dicer, Dgcr8 knockout (KO) ES cells have less severe phenotypes than Dicer knockout ES cells [13].

The Microprocessor was recently shown to have an additional role in directly destabilizing a mRNA target. Specifically, it can cleave hairpins in the 5′UTR and coding region of the Dgcr8 mRNA, which in turns destabilizes the mature transcript [14], [15], [16]. This negative feedback loop on Dgcr8 suggests the importance of tight homeostatic control of the Microprocessor in normal cellular function. The finding that the Microprocessor can directly influence Dgcr8 mRNA levels raises the possibility that this mechanism may affect many other mRNAs.

To further test whether there is a broader role of the Microprocessor in the direct regulation of mRNAs, we evaluated the mRNA and small non-coding RNA profiles of wild-type, Dgcr8 KO and Dicer KO cells as well as a recently published data set of small RNAs less than 200 nucleotides from human Hela and HepG2 cell lines [17]. While many mRNAs were differentially expressed between Dgcr8 and Dicer KO ES cells, there was no evidence for Microprocessor based processing of these mRNAs, with the striking exception of Dgcr8 itself. Similarly, analysis of the Hela and HepG2 data sets identified many sequence reads from the Dgcr8 hairpins showing a pattern consistent with Microprocessor activity, but none from any other predicted hairpins within spliced mRNAs. These findings suggest that the Microprocessor's role in directly regulating mRNA levels is specific to auto-regulation of Dgcr8, highlighting the importance of this negative feedback regulation of Microprocessor levels.

Results

mRNAs regulated by a direct Microprocessor cleavage mechanism should be upregulated in cells deficient for the Microprocessor, but not in Dicer deficient cells. Therefore, we evaluated coding mRNA profiling data from wild-type, Dgcr8 KO and Dicer KO mouse ES cells. Normalized mRNA levels in Dgcr8 KO and Dicer KO cells were compared to wild-type ES cells (Figure 1). Most mRNAs that were upregulated or downregulated were similarly altered in both mutants. However, similar to previous studies [14], [18], we found multiple mRNAs whose expression were specifically altered in cells that lacked Dgcr8. Using a false discovery rate of 5%, there were 778 transcripts there were upregulated in Dgcr8 KO cells relative to both wild-type and Dicer KO. There were 843 transcripts that were downregulated.

[Figure omitted. See PDF.]

Figure 1. Transcripts differentially regulated in Dgcr8 KO relative to WT and Dicer KO ES cells.

The sets of genes differentially up- and down- regulated in Dgcr8 KO relative to Dicer KO and WT ES cells were determined based on a cutoff of FDR <5%. Data are represented as a mean of 3 biological replicates of WT, Dgcr8 KO and Dicer KO arrays. Transcripts positive for EvoFold hairpin predictions and transcripts with 5 or more small RNAs mapping to their exons are shown (see legend). Arrow points to Dgcr8 expression levels, which, as expected, is down in Dgcr8 KO (exon 3 deletion results in premature termination codon and, hence, non-sense mediated RNA decay [33].

https://doi.org/10.1371/journal.pone.0006971.g001

If genes specifically upregulated in Dgcr8 KO cells are normally cleaved by the Microprocessor, there should be hairpin substrates for the complex within these mRNAs. Therefore, we searched for evolutionary conserved hairpins within these mRNAs using predictions generated by the EvoFold algorithm [16]. The 5′UTR hairpin in Dgcr8 was first identified by this method. EvoFold predictions are grouped based on their location in CDS, 5′UTR, 3′UTR, intron and intergenic regions. We determined mouse genome coordinates for EvoFold hairpins in CDS, 5′UTR and 3′UTR regions (see Methods), mapped them to the coding mRNA database, and compared the relative expression levels of all positive hits in Dgcr8 KO, Dicer KO, and wild-type ES cells (Figure 1). A total of 824 out of 23805 (3.5%) coding mRNAs contained predicted hairpins. Of these 824, 43 mRNAs were specifically upregulated in Dgcr8 KO cells, while 24 mRNAs were specifically downregulated in the Dgcr8 KO cells. Therefore, there was a subset of genes specifically upregulated in Dgcr8 KO cells that contain predicted hairpins and hence could be direct targets of the Microprocessor.

If hairpins within the Dgcr8 KO- upregulated gene set are indeed cleaved by the Microprocessor, we hypothesized that there would be Dgcr8-dependent small RNAs that map to these hairpins. Therefore, we evaluated ultra-high throughput deep sequencing data representing small RNAs ranging from 18-32 nucleotides from the wild-type, Dgcr8 KO and Dicer KO ES cells. As expected, multiple sequence reads mapped to the EvoFold predicted 5′UTR and coding region hairpins of Dgcr8 mRNA in WT cells (Figure 2). None of the reads mapping to the coding region hairpin were found in either Dgcr8 or Dicer KO libraries confirming their Dgcr8- and Dicer-dependence (Figure 2B). Interestingly, two sequence reads mapping to the 5′UTR hairpin were found in the Dicer KO library (Figure 2A). One of these reads mapped just 5′ to the hairpin. Such Dgcr8-dependent, Dicer-independent reads have been previously observed at miRNA loci in Drosophila and mouse small RNA sequencing studies and appear to be a 5′ remnant of Drosha cleavage that is further degraded by an unknown 5′-3′ exonuclease [10], [19]. The remaining read that was uncovered in the Dicer KO library had a 5′ end that did not map to the 5′ or 3′ end of the hairpin suggesting that it was a degradation product of the full length hairpin. Analysis of all EvoFold-predicted hairpins in the Dgcr8 KO-upregulated set of coding mRNAs failed to identify a single other hairpin with corresponding small RNAs.

[Figure omitted. See PDF.]

Figure 2. The distribution of reads across hairpins in the first exon of Dgcr8 in mES cells.

The location of each small RNA read relative to the exon is represented by a grey bar and was generated using the custom tracks feature on the UCSC genome browser. For each RNA species, the number of reads that were obtained with that sequence is indicated at the left. The predicted secondary structure is represented below the genomic sequence. Genomic coordinates are based on UCSC Known Genes annotations (mm8) (A) 5′UTR hairpin (B) CDS hairpin. Small RNA reads in WT cells are represented by black bars and small RNA reads in Dicer KO cells are represented by a grey bar.

https://doi.org/10.1371/journal.pone.0006971.g002

Analysis of only EvoFold predicted loci could miss poorly conserved hairpins. Therefore, to extend the analysis, sequencing reads from WT ES cells were mapped to all exons of the transcripts whose expression was altered in Dgcr8 KO versus WT and Dicer KO cells. 7 out of the 778 Dgcr8 KO- upregulated transcripts and 15 out of the 844 downregulated transcripts had at least 5 small RNA reads that overlapped with their exons (Figure 1). As Microprocessor activity is predicted to destabilize the mRNAs, we looked more closely at the 7 transcripts upregulated in Dgcr8 KO cells. The small RNAs that mapped within exonic regions of these annotated transcripts fell into two groups based on their distribution. Three had multiple small RNAs with a similar 5′ or 3′ end, consistent with specific endonuclease cleavage (Figure 3A and Figure S1, S2). The remaining five (two from the same transcript, Arrdc-3) had small RNAs mapping across the exon without shared 5′ or 3′ ends consistent with degradation (Figure 3B and Figure S3, S4, S5, S6). All of these small RNAs were present in the Dgcr8 null background (Figure 3A–B and Figure S1–S6). Hence, they are not products of Microprocessor cleavage.

[Figure omitted. See PDF.]

Figure 3. Representative examples of read distribution in exons with >5 reads in WT cells.

The location of unique small RNA reads from WT (black bars), Dgcr8 KO (dark grey bars) and Dicer KO (grey bars) are represented. For each RNA species, the number of reads that were obtained with that sequence is indicated at the left. Genomic coordinates are based on UCSC Known Genes annotations (mm8) (A) Example showing reads that are localized to a small window consistent with Microprocessor cleavage but are not Dgcr8-dependent (B) Reads are distributed across the range of the exon and most likely represent degradation products.

https://doi.org/10.1371/journal.pone.0006971.g003

A small number of annotated miRNAs map to exonic regions of coding genes (∼37 in mice) [20]. Therefore, analogous to Dgcr8, the host genes for these miRNAs might be expected to be downregulated by Microprocessor-induced cleavage. Upon examination of the exonic miRNAs, we found only 10 to fully lie within annotated exons (Table S1). We were able to find small RNA reads to three of these exonic miRNAs (mmu-miR-21, mmu-miR-671, mmu-miR-147). However, the mRNA levels of the host genes of these three miRNAs were not altered in the Dgcr8 and Dicer KO ES cells. Therefore, production of these miRNAs does not appear to influence the overall levels of the annotated host mRNAs. Together, these detailed analyses of both mRNA expression profiling and small RNA sequencing data from ES cells failed to uncover any genes other than Dgcr8 that are directly destabilized by the Microprocessor.

It is possible that 18–32 nucleotide small RNA sequencing missed Microprocessor-cleaved exonic hairpins that are sequestered and/or are not processed by Dicer. Microprocessor miRNAs are typically 60–75 nucleotides in length. Therefore, to directly identify these hairpins, we analyzed ultra high-throughput sequencing data sequence sets produced from small RNAs less than 200 nucleotides in length from Hela and HepG2 cells [17]. The forty small RNA libraries generated in the study were derived from whole cell, cytoplasmic and nuclear fractions, as well as from cells following enzymatic treatments that enrich for either mono-, di-, tri-phosphate modified or 5′ capped RNAs. Sequence reads from all forty libraries were mapped to exonic EvoFold hairpins. The largest number of hits, 184, mapped to the Dgcr8 5′UTR hairpin and 4 mapped to the coding region hairpin (Figure 4). Most of these reads had a uniform 5′ end consistent with Microprocessor cleavage. There was an additional read just 5′ to the hairpin, a likely remnant of the Microprocessor cleavage, similar to that seen in the ES cell small RNA libraries (Figure 2). A large number (166 out of 184) of the 5′ UTR reads were derived from nuclear libraries, consistent with previous work showing that the cleaved 5′UTR hairpin is confined to the nuclear fraction [14]. When mapping reads from the libraries to known pre-miRNA hairpins, many reads extend beyond the known mature miRNA into the loop region of the hairpin (Figure S7), thereby confirming that these libraries contain hairpin products of the Microprocessor cleavage. These findings show that the analysis of the Hela and HepG2 small RNA data sets should identify other hairpins that are cleaved by the Microprocessor even if they are not further processed.

[Figure omitted. See PDF.]

Figure 4. Read distribution across hairpins in the first exon of Dgcr8 in <200nt small RNA sequencing data from HeLa and HepG2 cells.

Small RNA locations are presented as in Figure 2. Genomic coordinates are based on UCSC Known Genes annotations (hg18) (A) 5′UTR hairpin (B) CDS hairpin.

https://doi.org/10.1371/journal.pone.0006971.g004

In order to identify any other potential mRNA substrates, we next mapped the HeLa and HepG2 datasets to all UTR and CDS EvoFold loci. There were 106 additional EvoFold hairpins containing overlapping small RNAs, although the number of reads mapping to any one of these hairpins was much less than seen for Dgcr8 (Table S2). Only four of these hairpins had at least 5 sequence reads. Furthermore, none of the small RNA reads in these hairpins mapped in a manner consistent with Microprocessor cleavage. That is, they had heterogeneous 5′ and 3′ends and/or the ends went beyond the extremes of the hairpins (Figure 5A–D). For example, the second highest-ranking hairpin, which mapped to the gene RPS3, had 14 reads. However, unlike the reads mapping to the Dgcr8 hairpins, they did not have a defined 5′ end, but instead mapped across the locus, more consistent with RNA degradation than Microprocessor cleavage. Therefore, analysis of small RNAs less than 200 nucleotides failed to identify any Evofold loci within exons other than Dgcr8 that are cleaved in a Microprocessor-like fashion.

[Figure omitted. See PDF.]

Figure 5. Read distribution across hairpins positive for >5 small RNA reads in HeLa cell <200 nt small RNA sequencing data.

Small RNA locations are presented as in Figure 2. Genomic coordinates are based on UCSC Known Genes annotations (hg18) Reads across hairpins in (A) RPS3 (B) HIST1H4C (C) RHOB and (D) RPS8.

https://doi.org/10.1371/journal.pone.0006971.g005

Again, limiting the analysis to Evofold predicted hairpins would miss non-conserved hairpins. Therefore, we mapped small RNAs from HeLa and HepG2 libraries to exons of transcripts upregulated over 2-fold with siRNA-mediated knockdown of both Drosha and Dgcr8 relative to siGFP. Expression information was extracted from recently published microarray data in HeLa cells (see Methods) [14]. As expected, Dgcr8, which was upregulated in the Drosha knockdown sample, had 188 small RNAs mapping to the first exon. Upon examining protein-coding genes upregulated in both Drosha and Dgcr8 knockdown samples, 31 transcripts had > = 10 small RNAs mapping to at least one exon (45 exons total, Table S3). Notably, 15 out of the 31 were genes that encode ribosomal protein subunits, which are highly abundant in cells. Out of the 31, 11 transcripts had small RNA reads distributed over the exon, as would be expected for degradation products. The remaining 20 transcripts had small RNA reads clustering in small window(s) within exons. However, further examination of the regions in these 20 transcripts using RNAfold did not reveal the presence of any good hairpin structures, in contrast to the Dgcr8 small RNA mapping-regions. In summary, analysis of the ultra high-throughput sequence reads of RNAs less than 200 nucleotides, like the ES cell small RNA dataset, showed that a role of the Microprocessor in direct mRNA regulation is likely limited to Dgcr8.

Discussion

Our findings show a focused role for the Microprocessor in destabilizing coding mRNAs by the direct cleavage and destabilization of spliced transcripts. Indeed, we only find evidence for the cleavage and destabilizaton of Dgcr8. Similar to previous reports, our mRNA profiling analysis of wild-type, Dgcr8, and Dicer deficient cells identifies many mRNAs that are specifically upregulated with the loss of Dgcr8 [18]. The presence of such mRNAs would be consistent with Microprocessor regulation of coding mRNAs through direct cleavage and may be a broadly used mechanism of mRNA regulation. However, closer analysis of these mRNAs and evaluation of ultra-high throughput deep sequencing for small RNAs either in the 18–32 or <200 nucleotide range failed to identify any additional mRNAs that are regulated by such a mechanism. We cannot exclude the possibility that rare examples of Microprocessor-mediated destabilization of mRNAs may be found in specific cellular contexts or at levels too low to be identified using current deep sequencing technology. However, in this study, we examined data from cell lines representing three different tissues: ES (inner cell mass of the blastocyst), Hela (kidney), and HepG2 (liver). Furthermore, in all deep sequencing datasets examined, we find numerous reads to the Dgcr8 hairpins but are unable to find a single additional similar candidate, suggesting that any additional examples would be extremely rare.

The absence of Dgcr8-dependent upregulation of the host transcripts carrying annotated exonic miRNAs is worth noting (Table S1). A number of these annotated exonic miRNAs were not present in our small RNA libraries even though the host gene is clearly expressed. This finding may be the result of mis-annotation of these sequences as miRNAs or that processing of the hairpins is somehow suppressed in ES cells. One example of an annotated exonic miRNA that is present in large numbers in ES cells is miR-21. Its host gene, Tmem49, is not upregulated in Dgcr8 or Dicer knockout ES cells. Possible explanations include: 1) there are alternative transcripts responsible for miR-21 production either from an alternative promoter or an alternative splicing event or 2) only a small subset of the Tmem49 transcripts is processed by the Microprocessor to produce the pre-miR-21 hairpin.

A very specific role for the Microprocessor in destabilizing Dgcr8 and hence providing a negative feedback on Microprocessor levels itself suggests that homeostatic control of microRNA processing is central to normal cellular physiology. This is consistent with recent findings showing that much regulation is occurring at the level of Microprocessor activity. For example hnRNAP, Lin28, and KSRP have been suggested to regulate Microprocessor activity on specific miRNAs [21], [22], [23], [24]. Furthermore, SMAD signaling alters the processing of pri- to pre-miR-21 [25]. A carefully controlled balance between the levels of the Microprocessor and these regulators are likely important for proper physiologic function.

Dgcr8 levels are differentially regulated during development and in cancers. Interestingly, some cancers have decreased, while other cancers have increased levels of Dgcr8 [26], [27]. Similarly, Dicer levels and/or activity appear to be altered in cancers [26], [27], [28], [29]. A direct role for changes in processing activity in cancer is supported by a mouse model of lung cancer where heterozygous loss of Dicer promotes tumor progression [30]. Together, these findings suggest that the biogenesis of miRNAs is not simply a passive process, but rather a tightly controlled one. Therefore, it will be important to determine in greater detail how the level and the activity of the biogenesis machinery influence the molecular constitution of cells.

Materials and Methods

Solexa sequencing data for Dgcr8 KO, Dicer KO and WT cells were previously published[10]. Information about exonic miRNAs and host genes was extracted from the CoGemir database [20].

Microarray analysis

Microarray experiments on the wild-type, Dgcr8 KO, and Dicer KO cells were performed by the Gladstone Genomic Core Facility using the Affymetrix 1.0 mouse gene ST arrays with 3 biological replicates per genotype (wild-type (v6.5), Dgcr8 knockout, Dicer knockout ES cells) [10], [31]. Dgcr8 and Dicer knockout ES cell derivation and culture has been previously described [10], [13]. Protocol used for preparation of RNA and hybridization for microarray has been previously described [31]. Array data was normalized using the robust multi-array average (RMA) algorithm. Normalized data has been deposited at GEO (#GSE16923). Genes upregulated and downregulated in Dgcr8 KO relative to WT and Dicer KO were determined by FDR analysis using the SAM software package from Stanford. (http://www-stat.stanford.edu/~tibs/SAM/). Specifically, two sets of genes were determined: 1) Genes upregulated in Dgcr8 KO relative to Dicer KO and 2) Genes upregulated in Dgcr8 KO relative to wild-type. Overlapping transcripts between these two sets of genes were assigned to the group upregulated in Dgcr8 KO relative to both Dicer KO and WT. Genes downregulated in Dgcr8 KO relative to Dicer and WT were determined using the same approach. For the analysis of overlap between small RNAs or predicted hairpins and protein coding mRNAs, we excluded Affy transcripts annotated only as miRNAs, transcripts mapping to the mitochondrial genome, chromosome Y and transcripts missing gene ID annotations.

HeLa cell microarray data was previously published [14]. siGFP, siDrosha and siDgcr8 expression data was averaged for 24 and 48 hr timepoints for each Affy ID, which resulted in 4 biological samples/gene. AffyIDs upregulated at least 2 fold (n = 1195) in both siDrosha and siDgcr8 relative to siGFP were analyzed further.

Mapping small RNA reads to exons

Small RNA reads from the Solexa sequencing dataset were first mapped to the genome (mouse, version mm8) using Eland. Uniquely mapping small RNA reads were mapped to exons by examining overlap between genomic coordinates of a small RNA read and each exon. Any small RNA overlapping with beginning and end of an exon as well as lying within in an exon was included as a positive hit. Exon information was determined using annotations from the UCSC Known Genes and Ensemble databases (mouse, version mm8) and all transcripts were collapsed to match to Affy ID annotations[32].

For analysis of data from HeLa and HepG2 cells, small RNAs from all libraries were first mapped to the genome (hg18) using Eland. Sequence length of HeLa cell libraries ranged from 15 to 26 nt. Sequence length of HepG2 cell libraries ranged from 15 to 36 nt. Genomic coordinates of the small RNAs were then mapped to exons of transcripts upregulated with siDrosha and siDgcr8 relative to siGFP. Exon information was determined using RefSeq annotations, which were matched to Affy IDs. The positive hits were further filtered manually of snoRNAs. The remaining exons were ranked based on the number of small RNA reads and exons containing >10 small RNAs were analyzed further using custom tracks at the UCSC genome browser. For exons with small RNA reads localized to a small window, sequences surrounding the small RNA reads were extracted and fold predictions were generated using RNAfold.

Mapping Small RNA reads to miRNA hairpins

Genomic locations of miRNA hairpins were extracted from miRBase and converted to the hg18 assembly using the liftover tool. Genomic coordinates of small RNA sequences (25 to 36 nt) from HeLa and HepG2 cells were mapped to miRNA hairpin locations.

EvoFold Analysis

Lists of long CDS and 5′UTR hairpins and their location in the human genome (mapping based on May 2004 release) were downloaded from the EvoFold database (available online at:http://www.cbse.ucsc.edu/~jsp/EvoFold/) [16]. The genomic coordinates were converted to the mouse genome (version mm8) using the LiftOver tool at the UCSC genome browser. Predicted hairpins were then mapped to mouse exons from UCSC known genes and Ensemble data sets and matched to the corresponding Affy IDs. Small RNAs were mapped to the hairpins using genomic coordinates using the same approach used when mapping small RNAs to exons.

For analysis of data from HeLa and HepG2 cells, EvoFold UTR and CDS hairpin coordinates were converted to this version using the Liftover tool at the UCSC genome browser [32]. Small RNAs were directly mapped to the hairpins as described earlier.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

https://doi.org/10.1371/journal.pone.0006971.s001

(0.66 MB TIF)

Table S2.

https://doi.org/10.1371/journal.pone.0006971.s002

(3.62 MB TIF)

Table S3.

https://doi.org/10.1371/journal.pone.0006971.s003

(1.56 MB TIF)

Figure S1.

Distribution of small RNA reads from Dgcr8 KO and WT libraries across the Hn1 exon. Small RNA locations are presented as in Figure 2 (WT reads = black bars, Dgcr8 KO reads = grey bars). Genomic coordinates are based on UCSC Known Genes annotations (mm8).

https://doi.org/10.1371/journal.pone.0006971.s004

(0.66 MB TIF)

Figure S2.

Distribution of small RNA reads from Dgcr8 KO and WT libraries across the Atbf1 exon, presented as in Figure S1.

https://doi.org/10.1371/journal.pone.0006971.s005

(0.31 MB TIF)

Figure S3.

Distribution of small RNA reads from Dgcr8 KO and WT libraries across the Adam23 exon, presented as in Figure S1.

https://doi.org/10.1371/journal.pone.0006971.s006

(0.26 MB TIF)

Figure S4.

Distribution of small RNA reads from Dgcr8 KO and WT libraries across the Zfp462 exon, presented as in Figure S1.

https://doi.org/10.1371/journal.pone.0006971.s007

(0.43 MB TIF)

Figure S5.

Distribution of small RNA reads from Dgcr8 KO and WT libraries across the Arrdc3 exon, presented as in Figure S1.

https://doi.org/10.1371/journal.pone.0006971.s008

(0.45 MB TIF)

Figure S6.

Distribution of small RNA reads from Dgcr8 KO and WT libraries across the Arrdc3 exon, presented as in Figure S1.

https://doi.org/10.1371/journal.pone.0006971.s009

(0.55 MB TIF)

Figure S7.

Distribution of 25 to 36 nt sequence reads from HeLa and HepG2 cell <200nt libraries across pre-miRNA hairpins. Locations are presented as in Figure 2. Genomic coordinates are based on mirbase annotations.

https://doi.org/10.1371/journal.pone.0006971.s010

(0.71 MB TIF)

Acknowledgments

We thank Michelle Cleary, Miguel Ramalho-Santos, and members of the Blelloch Lab for critical reading of the manuscript.

Author Contributions

Conceived and designed the experiments: AS RB. Performed the experiments: AS. Analyzed the data: AS RB. Wrote the paper: AS RB.

Citation: Shenoy A, Blelloch R (2009) Genomic Analysis Suggests that mRNA Destabilization by the Microprocessor Is Specialized for the Auto-Regulation of Dgcr8. PLoS ONE4(9): e6971. https://doi.org/10.1371/journal.pone.0006971

References

1. Babiarz JE, Blelloch R (2009) Small RNAs – their biogenesis, regulation and function in embryonic stem cells. JE BabiarzR. Blelloch2009Small RNAs – their biogenesis, regulation and function in embryonic stem cells.In StemBook (www.stembook.org/node/583). In StemBook (www.stembook.org/node/583).

2. Han J, Lee Y, Yeom KH, Kim YK, Jin H, et al. (2004) The Drosha-DGCR8 complex in primary microRNA processing. Genes Dev 18: 3016–3027.J. HanY. LeeKH YeomYK KimH. Jin2004The Drosha-DGCR8 complex in primary microRNA processing.Genes Dev1830163027

3. Lee Y, Ahn C, Han J, Choi H, Kim J, et al. (2003) The nuclear RNase III Drosha initiates microRNA processing. Nature 425: 415–419.Y. LeeC. AhnJ. HanH. ChoiJ. Kim2003The nuclear RNase III Drosha initiates microRNA processing.Nature425415419

4. Denli AM, Tops BB, Plasterk RH, Ketting RF, Hannon GJ (2004) Processing of primary microRNAs by the Microprocessor complex. Nature 432: 231–235.AM DenliBB TopsRH PlasterkRF KettingGJ Hannon2004Processing of primary microRNAs by the Microprocessor complex.Nature432231235

5. Gregory RI, Yan KP, Amuthan G, Chendrimada T, Doratotaj B, et al. (2004) The Microprocessor complex mediates the genesis of microRNAs. Nature 432: 235–240.RI GregoryKP YanG. AmuthanT. ChendrimadaB. Doratotaj2004The Microprocessor complex mediates the genesis of microRNAs.Nature432235240

6. Landthaler M, Yalcin A, Tuschl T (2004) The human DiGeorge syndrome critical region gene 8 and Its D. melanogaster homolog are required for miRNA biogenesis. Curr Biol 14: 2162–2167.M. LandthalerA. YalcinT. Tuschl2004The human DiGeorge syndrome critical region gene 8 and Its D. melanogaster homolog are required for miRNA biogenesis.Curr Biol1421622167

7. Bernstein E, Caudy AA, Hammond SM, Hannon GJ (2001) Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature 409: 363–366.E. BernsteinAA CaudySM HammondGJ Hannon2001Role for a bidentate ribonuclease in the initiation step of RNA interference.Nature409363366

8. Chendrimada TP, Gregory RI, Kumaraswamy E, Norman J, Cooch N, et al. (2005) TRBP recruits the Dicer complex to Ago2 for microRNA processing and gene silencing. Nature 436: 740–744.TP ChendrimadaRI GregoryE. KumaraswamyJ. NormanN. Cooch2005TRBP recruits the Dicer complex to Ago2 for microRNA processing and gene silencing.Nature436740744

9. Haase AD, Jaskiewicz L, Zhang H, Laine S, Sack R, et al. (2005) TRBP, a regulator of cellular PKR and HIV-1 virus expression, interacts with Dicer and functions in RNA silencing. EMBO Rep 6: 961–967.AD HaaseL. JaskiewiczH. ZhangS. LaineR. Sack2005TRBP, a regulator of cellular PKR and HIV-1 virus expression, interacts with Dicer and functions in RNA silencing.EMBO Rep6961967

10. Babiarz JE, Ruby JG, Wang Y, Bartel DP, Blelloch R (2008) Mouse ES cells express endogenous shRNAs, siRNAs, and other Microprocessor-independent, Dicer-dependent small RNAs. Genes Dev 22: 2773–2785.JE BabiarzJG RubyY. WangDP BartelR. Blelloch2008Mouse ES cells express endogenous shRNAs, siRNAs, and other Microprocessor-independent, Dicer-dependent small RNAs.Genes Dev2227732785

11. Tam OH, Aravin AA, Stein P, Girard A, Murchison EP, et al. (2008) Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes. Nature 453: 534–538.OH TamAA AravinP. SteinA. GirardEP Murchison2008Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes.Nature453534538

12. Watanabe T, Totoki Y, Toyoda A, Kaneda M, Kuramochi-Miyagawa S, et al. (2008) Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes. Nature 453: 539–543.T. WatanabeY. TotokiA. ToyodaM. KanedaS. Kuramochi-Miyagawa2008Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes.Nature453539543

13. Wang Y, Medvid R, Melton C, Jaenisch R, Blelloch R (2007) DGCR8 is essential for microRNA biogenesis and silencing of embryonic stem cell self-renewal. Nat Genet 39: 380–385.Y. WangR. MedvidC. MeltonR. JaenischR. Blelloch2007DGCR8 is essential for microRNA biogenesis and silencing of embryonic stem cell self-renewal.Nat Genet39380385

14. Han J, Pedersen JS, Kwon SC, Belair CD, Kim YK, et al. (2009) Posttranscriptional crossregulation between Drosha and DGCR8. Cell 136: 75–84.J. HanJS PedersenSC KwonCD BelairYK Kim2009Posttranscriptional crossregulation between Drosha and DGCR8.Cell1367584

15. Triboulet R, Chang HM, Lapierre RJ, Gregory RI (2009) Post-transcriptional control of DGCR8 expression by the Microprocessor. RNA 15: 1005–1011.R. TribouletHM ChangRJ LapierreRI Gregory2009Post-transcriptional control of DGCR8 expression by the Microprocessor.RNA1510051011

16. Pedersen JS, Bejerano G, Siepel A, Rosenbloom K, Lindblad-Toh K, et al. (2006) Identification and classification of conserved RNA secondary structures in the human genome. PLoS Comput Biol 2: e33.JS PedersenG. BejeranoA. SiepelK. RosenbloomK. Lindblad-Toh2006Identification and classification of conserved RNA secondary structures in the human genome.PLoS Comput Biol2e33

17. ENCODE (2009) Post-transcriptional processing generates a diversity of 5′-modified long and short RNAs. Nature 457: 1028–1032.ENCODE2009Post-transcriptional processing generates a diversity of 5′-modified long and short RNAs.Nature45710281032

18. Kadener S, Rodriguez J, Abruzzi KC, Khodor YL, Sugino K, et al. (2009) Genome-wide identification of targets of the drosha-pasha/DGCR8 complex. RNA. S. KadenerJ. RodriguezKC AbruzziYL KhodorK. Sugino2009Genome-wide identification of targets of the drosha-pasha/DGCR8 complex.RNA

19. Ruby JG, Stark A, Johnston WK, Kellis M, Bartel DP, et al. (2007) Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res 17: 1850–1864.JG RubyA. StarkWK JohnstonM. KellisDP Bartel2007Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs.Genome Res1718501864

20. Maselli V, Di Bernardo D, Banfi S (2008) CoGemiR: a comparative genomics microRNA database. BMC Genomics 9: 457.V. MaselliD. Di BernardoS. Banfi2008CoGemiR: a comparative genomics microRNA database.BMC Genomics9457

21. Michlewski G, Guil S, Semple CA, Caceres JF (2008) Posttranscriptional regulation of miRNAs harboring conserved terminal loops. Mol Cell 32: 383–393.G. MichlewskiS. GuilCA SempleJF Caceres2008Posttranscriptional regulation of miRNAs harboring conserved terminal loops.Mol Cell32383393

22. Viswanathan SR, Daley GQ, Gregory RI (2008) Selective blockade of microRNA processing by Lin28. Science 320: 97–100.SR ViswanathanGQ DaleyRI Gregory2008Selective blockade of microRNA processing by Lin28.Science32097100

23. Newman MA, Thomson JM, Hammond SM (2008) Lin-28 interaction with the Let-7 precursor loop mediates regulated microRNA processing. RNA 14: 1539–1549.MA NewmanJM ThomsonSM Hammond2008Lin-28 interaction with the Let-7 precursor loop mediates regulated microRNA processing.RNA1415391549

24. Trabucchi M, Briata P, Garcia-Mayoral M, Haase AD, Filipowicz W, et al. (2009) The RNA-binding protein KSRP promotes the biogenesis of a subset of microRNAs. Nature. M. TrabucchiP. BriataM. Garcia-MayoralAD HaaseW. Filipowicz2009The RNA-binding protein KSRP promotes the biogenesis of a subset of microRNAs.Nature

25. Davis BN, Hilyard AC, Lagna G, Hata A (2008) SMAD proteins control DROSHA-mediated microRNA maturation. Nature 454: 56–61.BN DavisAC HilyardG. LagnaA. Hata2008SMAD proteins control DROSHA-mediated microRNA maturation.Nature4545661

26. Ambs S, Prueitt RL, Yi M, Hudson RS, Howe TM, et al. (2008) Genomic profiling of microRNA and messenger RNA reveals deregulated microRNA expression in prostate cancer. Cancer Res 68: 6162–6170.S. AmbsRL PrueittM. YiRS HudsonTM Howe2008Genomic profiling of microRNA and messenger RNA reveals deregulated microRNA expression in prostate cancer.Cancer Res6861626170

27. Merritt WM, Lin YG, Han LY, Kamat AA, Spannuth WA, et al. (2008) Dicer, Drosha, and outcomes in patients with ovarian cancer. N Engl J Med 359: 2641–2650.WM MerrittYG LinLY HanAA KamatWA Spannuth2008Dicer, Drosha, and outcomes in patients with ovarian cancer.N Engl J Med35926412650

28. Chiosea S, Jelezcova E, Chandran U, Acquafondata M, McHale T, et al. (2006) Up-regulation of dicer, a component of the MicroRNA machinery, in prostate adenocarcinoma. Am J Pathol 169: 1812–1820.S. ChioseaE. JelezcovaU. ChandranM. AcquafondataT. McHale2006Up-regulation of dicer, a component of the MicroRNA machinery, in prostate adenocarcinoma.Am J Pathol16918121820

29. Melo SA, Ropero S, Moutinho C, Aaltonen LA, Yamamoto H, et al. (2009) A TARBP2 mutation in human cancer impairs microRNA processing and DICER1 function. Nat Genet 41: 365–370.SA MeloS. RoperoC. MoutinhoLA AaltonenH. Yamamoto2009A TARBP2 mutation in human cancer impairs microRNA processing and DICER1 function.Nat Genet41365370

30. Kumar MS, Lu J, Mercer KL, Golub TR, Jacks T (2007) Impaired microRNA processing enhances cellular transformation and tumorigenesis. Nat Genet 39: 673–677.MS KumarJ. LuKL MercerTR GolubT. Jacks2007Impaired microRNA processing enhances cellular transformation and tumorigenesis.Nat Genet39673677

31. Wang Y, Baskerville S, Shenoy A, Babiarz JE, Baehner L, et al. (2008) Embryonic stem cell-specific microRNAs regulate the G1-S transition and promote rapid proliferation. Nat Genet 40: 1478–1483.Y. WangS. BaskervilleA. ShenoyJE BabiarzL. Baehner2008Embryonic stem cell-specific microRNAs regulate the G1-S transition and promote rapid proliferation.Nat Genet4014781483

32. Karolchik D, Kuhn RM, Baertsch R, Barber GP, Clawson H, et al. (2008) The UCSC Genome Browser Database: 2008 update. Nucleic Acids Res 36: D773–779.D. KarolchikRM KuhnR. BaertschGP BarberH. Clawson2008The UCSC Genome Browser Database: 2008 update.Nucleic Acids Res36D773779

33. Rebbapragada I, Lykke-Andersen J (2009) Execution of nonsense-mediated mRNA decay: what defines a substrate? Curr Opin Cell Biol 21: 394–402.I. RebbapragadaJ. Lykke-Andersen2009Execution of nonsense-mediated mRNA decay: what defines a substrate?Curr Opin Cell Biol21394402