Chromatin structure characteristics of pre-miRNA genomic sequences
© Zhu et al; licensee BioMed Central Ltd. 2011
Received: 30 October 2010
Accepted: 25 June 2011
Published: 25 June 2011
MicroRNAs (miRNAs) are non-coding RNAs with important roles in regulating gene expression. Recent studies indicate that transcription and cleavage of miRNA are coupled, and that chromatin structure may influence miRNA transcription. However, little is known about the relationship between the chromatin structure and cleavage of pre-miRNA from pri-miRNA.
By analysis of genome-wide nucleosome positioning data sets from human and Caenorhabditis elegans (C. elegans), we found an enrichment of positioned nucleosome on pre-miRNA genomic sequences, which is highly correlated with GC content within pre-miRNA. In addition, obvious enrichments of three histone modifications (H2BK5me1, H3K36me3 and H4K20me1) as well as RNA Polymerase II (RNAPII) were observed on pre-miRNA genomic sequences corresponding to the active-promoter miRNAs and expressed miRNAs.
Our results revealed the chromatin structure characteristics of pre-miRNA genomic sequences, and implied potential mechanisms that can recognize these characteristics, thus improving pre-miRNA cleavage.
MicroRNAs (miRNAs) are a class of approximate 22-nt non-coding RNAs capable of binding to the 3' UTR of their target genes, resulting in transcriptional destabilization and/or post-transcriptional inhibition [1, 2]. Typically, miRNA biogenesis initiates at transcription by RNA Polymerase II (RNAPII) , and produces a primary transcript (pri-miRNA). A canonical pri-miRNA is composed of a stem, a terminal loop, and long flanking sequences. Next, pri-miRNA is cleaved by a class 2 RNase III named Drosha (and its cofactor DGCR8) to yield a precursor miRNA (pre-miRNA) with a characteristic hairpin structure. The Drosha-DGCR8 complex is termed as Microprocessor .
Eukaryotic genomes are packaged together with nucleosomes making up chromatin. This organization embodies two distinct structures: sequence and chromatin structure. At the aspect of sequence, the pri-miRNAs are quite diverse. For example, different sites of stem in pri-miRNA may have different numbers of unpaired nucleotides, and the length of terminal loop by bioinformatic prediction could range from 3-nt to more than 10-nt [5–7]. It is also worth noting that no common sequence motif associated with human pri-miRNAs has been discovered . Thus, a confusing issue is introduced: how can these diverse RNA structures be identified by the same complex and precisely cleaved to derive pre-miRNA. Recent studies indicate that the double-stranded stem and the unpaired flanking regions are crucial for the recognition and cleavage by Microprocessor [8–10]. However, the factors which contribute to the recognition of pre-miRNA genomic sequences are still far from being well-understood. Duan et al. indicated that the Drosha processing can be hindered by a single nucleotide polymorphism emerging in the stem of pri-miRNAs , whereas by studying the tumor-associated mutation, Diederichs et al. found that the Drosha processing remains intact even if the secondary structure of pri-miRNA is altered by many sequence variations found in human cancers . These facts highlight the flexibility of the pre-miRNA recognition.
Chromatin structure, mainly concerned with how and where the DNA is packaged around nucleosomes, is affected by chemical modifications of histone proteins, including acetylation, methylation, phosphorylation and ubiquitination . Saito et al. showed that some of these histone modifications can alter the miRNA expression level . Also, Scott et al. found that treatment with HDAC inhibitor (stabilizing histone acetylation) can result in a rapid change in miRNA expression . Recently, researchers employed transcriptionally active features to systematically identify promoters for known miRNA genes [16–18]. Exploiting these predicted miRNA promoters, Barski et al. found that histone modifications create a chromatin environment that can modulate transcriptional activation of miRNAs .
Traditionally, the sequence and chromatin structure are thought to function separately, sequence (primary transcript) at the RNA level, and chromatin structure at the DNA level. However, accumulating evidence suggests that processes occurring at these two levels are coupled. Recent research suggests that Drosha cleavage of pre-miRNAs from pri-miRNAs occurs co-transcriptionally, for both independently transcribed and intron-located miRNAs . Also, Jan et al. observed that the pri-miRNAs retained at transcription sites are processed more efficiently into pre-miRNAs than the released pri-miRNAs, and the efficiently processed pri-miRNAs are enriched in chromatin-associated regions . The results suggest that pri-miRNA processing is influenced by the chromatin structure when coupled with transcription. Recently, genome-wide data for nucleosome positioning and several epigenetic histone modifications have been released [21–23]. These data enable us to assess the distribution of nucleosome and histone modifications along pre-miRNA genomic sequences as well as their relationships to other genomic features. A pioneering study by Ozsolak et al. discovered an enrichment of nucleosome on pre-miRNA genomic sequence , suggesting that the pre-miRNA has specific chromatin characteristics. In this study, a series of statistical analyses about chromatin structure and related genomic features were performed, including nucleosome occupancy, specific epigenetic histone modifications, GC content, and RNAPII occupancy along pre-miRNA genomic sequences. Through these analyses, we found that the nucleosomes are enriched on all pre-miRNA genomic sequences corresponding to miRNAs experiencing different phrases within biogenesis (potentially transcribed, untranscribed, expressed and unexpressed miRNAs), and miRNAs located in different positions relative to protein-coding genes (intronic, intergenic, and exonic miRNAs). In addition, three histone modifications (H2BK5me1, H3K36me3 and H4K20me1) as well as RNAPII are obviously enriched along the pre-miRNA genomic sequences corresponding to potentially transcribed miRNAs and expressed miRNAs. The factors, which contribute to pre-miRNA recognition, are previously understood to only lie in the primary transcripts (RNA level), but our results suggest that they may also reside in the chromatin structure (DNA level), thus providing a new way to better appreciate the connection between chromatin structure, sequence and pre-miRNA cleavage.
Nucleosome enrichment on pre-miRNA genomic sequence
We derived the mature miRNA expression profile and putative miRNA promoters of CD4+ T cells from Ref. . These promoters are predicted employing characteristic transcriptionally active features for genes (RNAPII, H3K4me3 and H2A.Z), which are indirect reflection of transcription. Based on predicted miRNA promoters and mature miRNA expression profiles, we derived four groups of miRNAs, including 230 active-promoter miRNAs, 243 inactive-promoter miRNAs, 206 expressed miRNAs and 329 unexpressed miRNAs (group definition details are seen in Methods). It is noted that overlaps were allowed between the categories to address the different steps during the miRNA biogenesis. The group annotation for each miRNA is presented in Additional file 2. The nucleosome occupancy is assessed in these four groups, respectively. As shown in Figure 1C, we discovered higher nucleosome occupancy on pre-miRNA genomic sequences corresponding to miRNAs of all three groups than their flanking regions (Wilcoxon signed-rank test p-value < 2.2e-12 for active-promoter, expressed and unexpressed miRNAs, and p-value = 3.761e-12 for inactive-promoter miRNAs), and found that the enrichments of nucleosome occupancy in active-promoter miRNAs and expressed miRNAs are significantly higher than the inactive-promoter miRNAs and unexpressed miRNAs (Wilcoxon signed-rank test p-value = 1.676e-06 for comparison between active-promoter and inactive-promoter, and p-value = 0.0001272 for expressed and unexpressed). In addition, we divided pre-miRNAs into five groups according to their corresponding mature miRNA expression levels: the pre-miRNAs are treated as being with "Lowest expression", whose expression levels are zero (no mature miRNAs found), and the remaining pre-miRNAs are divided into four groups with equal size. As shown in Figure 1D, nucleosome occupancies are almost uniformly distributed in each of the five groups. Furthermore, we evaluated the correlation between mature miRNA expression levels and corresponding nucleosome occupancies using Kendall's rank correlation test (tau = 0.1134140, p-value = 0.0004294). These facts suggest that there is no significant dependency between nucleosome occupancy and mature miRNA expression levels, consistently with the previous study in three human cancer cell lines . A detailed list of miRNAs with regards to their expression profiles and nucleosome occupancies was presented in Additional file 2.
Histone modification enrichment on human pre-miRNA genomic sequence
Effect of GC content on nucleosome occupancy
To further investigate the correlation between DNA sequence and nucleosome occupancy, we divided pre-miRNAs into five equal bins according to the GC content in a 400-nt window surrounding the center of pre-miRNA genomic sequence, and examined the correlation between GC content and nucleosome occupancy (the details are in the method section). As shown in Figure 3C, nucleosome occupancy was increased with the elevated GC content on pre-miRNA genomic sequence, except the group with the highest GC content (60%~82%). This is similar to the observation on exons in human CD4+ T cells . The facts above suggest that GC composition could be one of the factors that influence nucleosome occupancy on pre-miRNA genomic sequence.
RNA Polymerase II enrichment on human pre-miRNA genomic sequence
Our analyses revealed that nucleosome occupancy is significantly higher in regions of the genome encoding pre-miRNAs of human CD4+ T cells with regards to their surrounding regions. The nucleosome enrichments are observed on all pre-miRNA genomic sequences corresponding to miRNAs experiencing different phrases within biogenesis (potentially transcribed, untranscribed, expressed and unexpressed miRNAs), and miRNAs located in different positions relative to protein-coding genes (intronic, intergenic, and exonic miRNAs). The similar enrichment of nucleosome is discovered in C. elegans, implying that this phenomenon is evolutionarily conserved. In addition, high distributions of three histone modifications (H2BK5me1, H3K36me3 and H4K20me1) as well as RNAPII were observed on pre-miRNA genomic sequences corresponding to potentially transcribed miRNAs and expressed miRNAs.
Nucleosome formation potentially depends significantly on the GC content of DNA sequence [29–31]. Consistently, Ozsolak et al. observed that GC content within pre-miRNA genomic sequence is significantly higher than the surrounding regions . Moreover, we found that the high distributions emerge in all pre-miRNA genomic sequences located in different positions relative to protein-coding genes (intronic, intergenic, and exonic pre-miRNAs). In addition, we found a high correlation between GC content and nucleosome occupancy within pre-miRNA genomic sequence. These facts suggest GC content could be one of the potential factors that influence nucleosome occupancy in pre-miRNA genomic sequence.
Our results revealed chromatin structure characteristics of pre-miRNA genomic sequences. Furthermore, molecular mechanisms should be required to recognize the characteristics, whereas some of the recognition mechanisms just reside in the cleavage mechanisms, since the recognition of pre-miRNA for organism denotes cleavage of pre-miRNA to some extent. There are two possible mechanisms for chromatin structure to improve miRNA cleavage. First, the presence of a positioned nucleosome within the pre-miRNA genomic sequence can act as 'speed bump' to reduce the rate of RNAPII elongation, thus improving pre-miRNA cleavage. Our analyses revealed that nucleosome occupancy is significantly higher in regions of the genome corresponding to the pre-miRNAs, also, we found that there exists enrichment of RNAPII within pre-miRNA, consistently with previous results in several specific miRNAs of Hela cell line . This enrichment could result from the reduction of RNAPII elongation rate caused by nucleosome bound on pre-miRNA. Furthermore, the reduction of RNAPII elongation rate has been shown to promote splicing . Also, due to the rate reduction, the nascent transcript could retain at DNA template longer. Pawlicki et al. indicated that the retention at chromatin could facilitate pri-miRNA processing . Our results showed that the nucleosome occupancies in the genomic sequences of active-promoter miRNAs and expressed miRNAs are significantly higher than the inactive-promoter miRNAs and unexpressed miRNAs. These results suggest that nucleosomes bound on pre-miRNA genomic sequence could slow the elongation rate of RNAPII, thereby improving pre-miRNA cleavage. Similar phenomenon has been reported in some recent studies in splicing site selection for protein-coding genes [25, 26, 40]. In those studies, nucleosome enrichment is found within exons, compared to introns, and this enrichment within exons can slow the elongation rate of RNAPII to promote splicing. In our studies, we further made a comparison between the nucleosome occupancies within pre-miRNA and exons, and the corresponding occupancy patterns are similar. This comparison strengthens our hypothesis. However, we did not find a significant correlation between nucleosome occupancy and RNAPII level. This is probably because RNAPII level within pre-miRNA genomic sequence is determined not only by nucleosome occupancy, but also by the transcription level of corresponding miRNA. Thus, their relationship can be evaluated, only when the transcription influence is excluded. Similarly, the correlation was not observed between nucleosome occupancy in pre-miRNA genomic sequence and corresponding mature miRNA expression level, consistently with the previous study . This is probably because miRNAs are regulated at multiple levels through their biogenesis . Particularly, pre-miRNA is processed by RNase III endonuclease Dicer to generate the mature miRNA. This process takes place in cytoplasm, which is not influenced by the chromatin structure. More detailed analyses and relevant dataset publication are required to address the relationship between nucleosome occupancy, RNAPII and miRNA expression levels.
Second, some specific histone modifications on pre-miRNA genomic sequence may contribute to pre-miRNA cleavage. Our results showed that there exist enrichments of three specific histone modifications (H2BK5me1, H3K36me3 and H4K20me1) on pre-miRNA genomic sequences corresponding to active-promoter miRNAs and expressed miRNAs. Consistently, several studies have demonstrated that these three histone modifications are the robust transcriptionally active characteristics [22, 28, 42–44], whereas miRNA cleavage recently is indicated to take place co-transcriptionally [19, 20]. Therefore, we hypothesize that these specific features may influence the pri-miRNA cleavage. In addition, numerous Microprocessor-associated proteins were identified as splicing factors [45, 46]. Kornblihtt indicates that diverse splicing factors associate with a carboxy terminal domain (CTD) of RNAPII , whereas phosphorylated CTD of RNAPII can recruit Histone-lysine N-methyltransferase SETD2, which modulates the histone modification H3K36me3 . Thus, H3K36me3 may interact with Microprocessor via the mediation of CTD, thus influencing the pre-miRNA cleavage. Moreover, Morlando et al. found an enrichment of Drosha in pre-miRNA, and this enrichment correlates with an enrichment of RNAPII . Considering the cross-talk between H3K36me3 and RNAPII mentioned above, this fact also suggests that specific histone modifications may contribute to pre-miRNA cleavage. Recent study indicated that H3K36me3 can interact with Polypyrimidine tract-binding protein (PTBP) , and moreover, the change of PTBP levels can differentially influence specific miRNA expression . These facts further validated our hypothesis.
Recently, several studies have been carried out which suggest that exons carry enrichments of nucleosome positioning and specific histone modifications [25, 26, 40, 51], implying that chromatin structure can also modulate exon selection in the process of mRNA maturation. In these studies, similar mechanisms to our study are pointed out. Some similarities can be found between mRNA maturation and miRNA biogenesis. 1) At the DNA level, both of transcriptions rely on RNAPII. 2) At the RNA level, pre-mRNA processing takes place in a large complex known as Spliceosome , whereas Spliceosome and Microprocessor share numerous key components [45, 53, 54]. 3) Recent researches demonstrate that pre-miRNAs are also generated co-transcriptionally, similar to protein-coding genes [19, 20]. These similarities also suggest that nucleosome bound on pre-miRNA and exons could share the similar mechanism for splicing.
Some groups indicated that introns carry low nucleosome occupancy, histone modification, and GC content [25, 26, 40]. However, our results showed that miRNA-coding intron has an enrichment of nucleosome, and this nucleosome carries specific histone modifications, such as H2BK5me1, H3K36me3 and H4K20me1. In addition, we found that GC content is much higher within intronic pre-miRNAs than the flanking regions in the corresponding introns. These differences may be critical to the biogenesis of intronic pre-miRNAs, since both of the two possible mechanisms of pre-miRNA cleavage mentioned above depend on nucleosome occupancy and specific histone modification. Moreover, many studies suggested that mature miRNAs and mRNAs derive from a common nascent transcript [19, 39, 55]. Consequently, the cleavage of pre-miRNA is possible to influence exon splicing. Given that intronic pre-miRNA is cleaved by Drosha with the help of chromatin structure, the exposed host intron transcript will be rapidly degraded by RNA exonucleases, thus enhancing the exon splicing efficiency . Thus, due to the nucleosome enrichment within intronic pre-miRNA, the corresponding splicing signals may refrain from disruption.
A pioneering study by Ozsolak et al. discovered the nucleosome enrichment on pre-miRNA genomic sequence , suggesting that pre-miRNA has specific chromatin structure characteristics. However, the characteristics reside in the DNA level, whereas the cleavage of pre-miRNA occurs in the RNA level. Thus, a confusing problem is raised, how these chromatin structure characteristics can be recognized. Recent study indicated these two levels are coupled for pre-miRNA cleavage . Inspired by the facts above, a series of analyses was performed in this paper. Enrichments of three specific histone modifications as well as RNAPII were also observed on pre-miRNA genomic sequences corresponding to active-promoter miRNAs and expressed miRNAs. These facts together with relevant analyses implied an RNAPII-mediated mechanism that recognizes the chromatin structure characteristics, and links up DNA and RNA levels, thus improving pre-miRNA cleavage.
MiRNAs are non-coding RNAs with important roles in regulating gene expression. Traditionally, the sequence and chromatin structure are thought to function separately, sequence at RNA level, and chromatin structure at the DNA level. However, accumulating evidence suggests that processes occurring at these two levels are coupled. Recent research suggests that Drosha cleavage of pre-miRNAs from pri-miRNAs occurs co-transcriptionally. However, little is known about the relationship between the chromatin structure and cleavage of pre-miRNA from pri-miRNA. By analysis of genome-wide nucleosome positioning data sets from human and C. elegans, we found an enrichment of positioned nucleosome in pre-miRNA, which is highly correlated with GC content within pre-miRNA. In addition, by analysis of genome-wide histone modification data sets from human, we discovered obvious enrichments of three histone modifications (H2BK5me1, H3K36me3 and H4K20me1) on pre-miRNA genomic sequences corresponding active-promoter miRNAs and expressed miRNAs. Moreover, we found enrichment of RNAPII in pre-miRNA genomic sequence. Our results revealed the chromatin structure characteristics of pre-miRNA genomic sequences, and imply potential mechanisms to recognize the characteristics, thus improving pre-miRNA cleavage. There are two possible explanations. First, the presence of a positioned nucleosome within the pre-miRNA can act as 'speed bump' to reduce the rate of RNAPII elongation to improve pre-miRNA cleavage. Second, some specific histone modifications on pre-miRNA may contribute to pre-miRNA cleavage by recruiting Microprocessor. The crucial factors for pre-miRNA recognition are previously understood to only lie in the sequence (RNA level), but our results suggest that they may also reside in the chromatin structure (DNA level), thus providing a new way to better appreciate the connection between chromatin structure, sequence and pre-miRNA cleavage.
Pre-miRNA annotation and miRNA expression profile
We downloaded the annotation of 715 human pre-miRNAs from miRBase  and human refGene from the UCSC genome browser  (both relative to the hg18 version of the human genome). By comparing the positions of pre-miRNAs obtained from miRBase with the positions of exons obtained from refGene, we obtained intronic, intergenic and exonic pre-miRNAs. In addition, we derived the mature miRNA expression profile and putative miRNA promoters of CD4+ T cells from Ref. . These promoters are predicted employing characteristic transcriptionally active features (RNAPII, H3K4me3 and H2A.Z obtained from ), which are the indirect reflection of miRNA transcription. According to mature miRNA expression and miRNA promoter prediction, we obtained four groups of miRNAs, expressed miRNAs, unexpressed miRNAs, active-promoter miRNAs and inactive-promoter miRNAs. miRNAs are expressed, if reads by deep sequencing corresponding to miRNAs are detected; otherwise, miRNAs are unexpressed. Active promoters of miRNAs are identified in 250,000-nt upstream of pre-miRNAs, if all three enrichments of RNAPII, H3K4me3 and H2A.Z are observed, and enrichment of RNAPII co-localizes with the one of H2A.Z and/or H3K4me3, or else, enrichment of H3K4me3 co-localizes with the one of H2A.Z in the upstream regions of the miRNA genomic sequences; otherwise, if miRNAs are unexpressed, and corresponding active promoters are not identified, they are treated as being with inactive promoters. The promoter prediction and enrichment identification of RNAPII, H3K4me3 and H2A.Z refer to .
Nucleosome and histone modification data
We downloaded three datasets of nucleosome positioning and histone modifications from [21–23] (relative to the hg18 version of the human genome). We also obtained genomic sequence of C. elegans from the UCSC table browser. Adjusted coverage of nucleosome occupancy in C. elegans was obtained from  (relative to the ce6 version of the C. elegans genome).
We aligned the centers of every pre-miRNA genomic sequences, and defined the center as position zero, downstream positions as plus, and upstream positions as minus. We calculated the scores by assigning each 10-nt sliding window the average of the scores mapping to 140-nt surrounding this window.
To assess whether the enrichment of specific histone modification within pre-miRNA is significant compared to its flanking regions, we made a comparison between two kinds of segments:
2. 2,000-nt windows in the flanking regions of pre-miRNA genomic sequences. We chose different regions for intronic, intergenic and exonic miRNAs, respectively. For intronic miRNA, select every 2,000 nucleotides as a window, from the boundary of the center window above to the splicing sites of the host intron. For intergenic and exonic miRNAs, select two 2,000-nt windows away from the two boundaries of the center window, respectively. (refer to Figure 6B for details)
We made comparison between the nucleosome occupancies within 400-nt windows surrounding the centers of these two segments. The statistical significance of the difference was evaluated by the Wilcoxon signed-rank test (refer to Figure 6C for details). Also, we calculated the mean nucleosome occupancy within a 400-nt window surrounding the center of each miRNA, which is presented in Additional file 2. Based on this dataset, we respectively employed Kendall's rank correlation test to evaluate the correlation between nucleosome occupancies on pre-miRNA genomic sequences and the corresponding mature miRNA expression profiles, and utilized Wilcoxon signed-rank test to make comparison between nucleosome occupancies on pre-miRNA genomic sequences corresponding to miRNAs experiencing different phrases in biogenesis (potentially transcribed, untranscribed, expressed and unexpressed miRNAs) and located in different positions relative to protein-coding genes (intronic, intergenic and exonic miRNAs).
Comparison between nucleosome occupancy on pre-miRNAs, exons and introns
To make comparison between nucleosome occupancy on pre-miRNA, exon and intron, we summed up the nucleosome occupancy on every pre-miRNA, exon and intron, respectively. Since the lengths of pre-miRNA, exon and intron are different, we divided the sum of nucleosome occupancy of every pre-miRNA, exon or intron by its length, respectively, to obtain the corresponding distributions of nucleosome occupancy for each nucleotide. The statistical significance of the difference between nucleosome occupancies on pre-miRNA, exon and intron was evaluated by the Wilcoxon signed-rank test.
where N is the number of nucleotide G or C in S, and T is the total number of nucleotides in S. We assigned score of position i the GC content of a 100-nt window surrounding this position. The calculation of GC dinucleotides, AT and AT dinucleotides is similar to that of GC content.
To observe the correlation between GC content and nucleosome occupancy, we carried out the following work. 1) Selected the same group of segments utilized for assessing the difference of histone modifications between pre-miRNA and its surrounding regions (which is mentioned in the "Nucleosome and histone modification data" subsection). 2) Calculated the GC-scores within a 400-nt window surrounding the centers of segments. 3) According to the means of the GC-scores, divided each group of segments into five equal bins. 4) Calculated the nucleosome occupancy in each bin.
We used the R language (R Development Core Team 2008) to perform all analyses, using standard methods and packages.
RNA Polymerase II
monomethylation of H2BK5
trimethylation of H3K4
trimethylation of H3K36
monomethylation of H4K20
carboxy terminal domain
Polypyrimidine tract-binding protein.
This work is supported by the National Natural Science Foundation of China (60973078, 60901075), the Hi-Tech Research and Development Program of China (2007AA02Z302) and the Natural Science Foundation of Heilongjiang Province of China (LC2009C35).
- Lai EC: Micro RNAs are complementary to 3' UTR sequence motifs that mediate negative post-transcriptional regulation. Nat Genet. 2002, 30 (4): 363-364. 10.1038/ng865.PubMedView Article
- Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002.PubMed CentralPubMedView Article
- Lee Y, Kim M, Han J, Yeom KH, Lee S, Baek SH, Kim VN: MicroRNA genes are transcribed by RNA polymerase II. EMBO J. 2004, 23 (20): 4051-4060. 10.1038/sj.emboj.7600385.PubMed CentralPubMedView Article
- Seitz H, Zamore PD: Rethinking the microprocessor. Cell. 2006, 125 (5): 827-829. 10.1016/j.cell.2006.05.018.PubMedView Article
- Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T: Identification of novel genes coding for small expressed RNAs. Science. 2001, 294 (5543): 853-858. 10.1126/science.1064921.PubMedView Article
- Lau NC, Lim LP, Weinstein EG, Bartel DP: An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science. 2001, 294 (5543): 858-862. 10.1126/science.1065062.PubMedView Article
- Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001, 294 (5543): 862-864. 10.1126/science.1065329.PubMedView Article
- Han J, Lee Y, Yeom KH, Nam JW, Heo I, Rhee JK, Sohn SY, Cho Y, Zhang BT, Kim VN: Molecular basis for the recognition of primary microRNAs by the Drosha-DGCR8 complex. Cell. 2006, 125 (5): 887-901. 10.1016/j.cell.2006.03.043.PubMedView Article
- Zeng Y, Cullen BR: Sequence requirements for micro RNA processing and function in human cells. RNA. 2003, 9 (1): 112-123. 10.1261/rna.2780503.PubMed CentralPubMedView Article
- Zeng Y, Cullen BR: Efficient processing of primary microRNA hairpins by Drosha requires flanking nonstructured RNA sequences. J Biol Chem. 2005, 280 (30): 27595-27603. 10.1074/jbc.M504714200.PubMedView Article
- Duan R, Pak C, Jin P: Single nucleotide polymorphism associated with mature miR-125a alters the processing of pri-miRNA. Hum Mol Genet. 2007, 16 (9): 1124-1131. 10.1093/hmg/ddm062.PubMedView Article
- Diederichs S, Haber DA: Sequence variations of microRNAs in human cancer: alterations in predicted secondary structure do not affect processing. Cancer Res. 2006, 66 (12): 6097-6104. 10.1158/0008-5472.CAN-06-0537.PubMedView Article
- Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128 (4): 693-705. 10.1016/j.cell.2007.02.005.PubMedView Article
- Saito Y, Liang G, Egger G, Friedman JM, Chuang JC, Coetzee GA, Jones PA: Specific activation of microRNA-127 with downregulation of the proto-oncogene BCL6 by chromatin-modifying drugs in human cancer cells. Cancer Cell. 2006, 9 (6): 435-443. 10.1016/j.ccr.2006.04.020.PubMedView Article
- Scott GK, Mattie MD, Berger CE, Benz SC, Benz CC: Rapid alteration of microRNA levels by histone deacetylase inhibition. Cancer Res. 2006, 66 (3): 1277-1281. 10.1158/0008-5472.CAN-05-3632.PubMedView Article
- Ozsolak F, Poling LL, Wang Z, Liu H, Liu XS, Roeder RG, Zhang X, Song JS, Fisher DE: Chromatin structure analyses identify miRNA promoters. Genes Dev. 2008, 22 (22): 3172-3183. 10.1101/gad.1706508.PubMed CentralPubMedView Article
- Wang G, Wang Y, Shen C, Huang YW, Huang K, Huang TH, Nephew KP, Li L, Liu Y: RNA polymerase II binding patterns reveal genomic regions involved in microRNA gene regulation. PLoS One. 2010, 5 (11): e13798-10.1371/journal.pone.0013798.PubMed CentralPubMedView Article
- Barski A, Jothi R, Cuddapah S, Cui K, Roh TY, Schones DE, Zhao K: Chromatin poises miRNA- and protein-coding genes for expression. Genome Res. 2009, 19 (10): 1742-1751. 10.1101/gr.090951.109.PubMed CentralPubMedView Article
- Morlando M, Ballarino M, Gromak N, Pagano F, Bozzoni I, Proudfoot NJ: Primary microRNA transcripts are processed co-transcriptionally. Nat Struct Mol Biol. 2008, 15 (9): 902-909. 10.1038/nsmb.1475.PubMedView Article
- Pawlicki JM, Steitz JA: Primary microRNA transcript retention at sites of transcription leads to enhanced microRNA production. Journal of Cell Biology. 2008, 182 (1): 61-76. 10.1083/jcb.200803111.PubMed CentralPubMedView Article
- Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, Wei G, Zhao K: Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132 (5): 887-898. 10.1016/j.cell.2008.02.022.PubMedView Article
- Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009.PubMedView Article
- Wang Z, Zang C, Rosenfeld JA, Schones DE, Barski A, Cuddapah S, Cui K, Roh TY, Peng W, Zhang MQ, Zhao K: Combinatorial patterns of histone acetylations and methylations in the human genome. Nat Genet. 2008, 40 (7): 897-903. 10.1038/ng.154.PubMed CentralPubMedView Article
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36 Database: D154-158.
- Schwartz S, Meshorer E, Ast G: Chromatin organization marks exon-intron structure. Nat Struct Mol Biol. 2009, 16 (9): 990-995. 10.1038/nsmb.1659.PubMedView Article
- Tilgner H, Nikolaou C, Althammer S, Sammeth M, Beato M, Valcarcel J, Guigo R: Nucleosome positioning as a determinant of exon recognition. Nature Structural & Molecular Biology. 2009, 16 (9): 996-U124. 10.1038/nsmb.1658.View Article
- Valouev A, Ichikawa J, Tonthat T, Stuart J, Ranade S, Peckham H, Zeng K, Malek JA, Costa G, McKernan K, Sidow A, Fire A, Johnson SM: A high-resolution, nucleosome position map of C. elegans reveals a lack of universal sequence-dictated positioning. Genome Res. 2008, 18 (7): 1051-1063. 10.1101/gr.076463.108.PubMed CentralPubMedView Article
- Vakoc CR, Sachdeva MM, Wang HX, Blobel GA: Profile of histone lysine methylation across transcribed mammalian chromatin. Molecular and Cellular Biology. 2006, 26 (24): 9185-9195. 10.1128/MCB.01529-06.PubMed CentralPubMedView Article
- Segal E, Fondufe-Mittendorf Y, Chen L, Thastrom A, Field Y, Moore IK, Wang JP, Widom J: A genomic code for nucleosome positioning. Nature. 2006, 442 (7104): 772-778. 10.1038/nature04979.PubMed CentralPubMedView Article
- Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, Segal E: The DNA-encoded nucleosome organization of a eukaryotic genome. Nature. 2009, 458 (7236): 362-366. 10.1038/nature07667.PubMed CentralPubMedView Article
- Peckham HE, Thurman RE, Fu Y, Stamatoyannopoulos JA, Noble WS, Struhl K, Weng Z: Nucleosome positioning signals in genomic DNA. Genome Res. 2007, 17 (8): 1170-1177. 10.1101/gr.6101007.PubMed CentralPubMedView Article
- Johnson SM, Tan FJ, McCullough HL, Riordan DP, Fire AZ: Flexibility and constraint in the nucleosome core landscape of Caenorhabditis elegans chromatin. Genome Res. 2006, 16 (12): 1505-1516. 10.1101/gr.5560806.PubMed CentralPubMedView Article
- Denli AM, Tops BB, Plasterk RH, Ketting RF, Hannon GJ: Processing of primary microRNAs by the Microprocessor complex. Nature. 2004, 432 (7014): 231-235. 10.1038/nature03049.PubMedView Article
- McGhee JD, Felsenfeld G: Nucleosome structure. Annu Rev Biochem. 1980, 49: 1115-1156. 10.1146/annurev.bi.49.070180.005343.PubMedView Article
- Allemand E, Batsche E, Muchardt C: Splicing, transcription, and chromatin: a menage a trois. Curr Opin Genet Dev. 2008, 18 (2): 145-151. 10.1016/j.gde.2008.01.006.PubMedView Article
- de la Mata M, Alonso CR, Kadener S, Fededa JP, Blaustein M, Pelisch F, Cramer P, Bentley D, Kornblihtt AR: A slow RNA polymerase II affects alternative splicing in vivo. Mol Cell. 2003, 12 (2): 525-532. 10.1016/j.molcel.2003.08.001.PubMedView Article
- Kornblihtt AR: Coupling transcription and alternative splicing. Adv Exp Med Biol. 2007, 623: 175-189. 10.1007/978-0-387-77374-2_11.PubMedView Article
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30 (1): 207-210. 10.1093/nar/30.1.207.PubMed CentralPubMedView Article
- Kim YK, Kim VN: Processing of intronic microRNAs. EMBO J. 2007, 26 (3): 775-783. 10.1038/sj.emboj.7601512.PubMed CentralPubMedView Article
- Andersson R, Enroth S, Rada-Iglesias A, Wadelius C, Komorowski J: Nucleosomes are well positioned in exons and carry characteristic histone modifications. Genome Res. 2009, 19 (10): 1732-1741. 10.1101/gr.092353.109.PubMed CentralPubMedView Article
- Winter J, Jung S, Keller S, Gregory RI, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009, 11 (3): 228-234. 10.1038/ncb0309-228.PubMedView Article
- Bell O, Wirbelauer C, Hild M, Scharf AN, Schwaiger M, MacAlpine DM, Zilbermann F, van Leeuwen F, Bell SP, Imhof A, Garza D, Peters AH, Schubeler D: Localized H3K36 methylation states define histone H4K16 acetylation during transcriptional elongation in Drosophila. EMBO J. 2007, 26 (24): 4974-4984. 10.1038/sj.emboj.7601926.PubMed CentralPubMedView Article
- Du HN, Fingerman IM, Briggs SD: Histone H3 K36 methylation is mediated by a trans-histone methylation pathway involving an interaction between Set2 and histone H4. Genes & Development. 2008, 22 (20): 2786-2798. 10.1101/gad.1700008.View Article
- Talasz H, Lindner HH, Sarg B, Helliger W: Histone H4-lysine 20 monomethylation is increased in promoter and coding regions of active genes and correlates with hyperacetylation. Journal of Biological Chemistry. 2005, 280 (46): 38814-38822. 10.1074/jbc.M505563200.PubMedView Article
- Shiohama A, Sasaki T, Noda S, Minoshima S, Shimizu N: Nucleolar localization of DGCR8 and identification of eleven DGCR8-associated proteins. Exp Cell Res. 2007, 313 (20): 4196-4207. 10.1016/j.yexcr.2007.07.020.PubMedView Article
- Wen X, Tannukit S, Paine ML: TFIP11 interacts with mDEAH9, an RNA helicase involved in spliceosome disassembly. Int J Mol Sci. 2008, 9 (11): 2105-2113. 10.3390/ijms9112105.PubMed CentralPubMedView Article
- Kornblihtt AR: Chromatin, transcript elongation and alternative splicing. Nat Struct Mol Biol. 2006, 13 (1): 5-7. 10.1038/nsmb0106-5.PubMedView Article
- Edmunds JW, Mahadevan LC, Clayton AL: Dynamic histone H3 methylation during gene induction: HYPB/Setd2 mediates all H3K36 trimethylation. EMBO J. 2008, 27 (2): 406-420. 10.1038/sj.emboj.7601967.PubMed CentralPubMedView Article
- Luco RF, Pan Q, Tominaga K, Blencowe BJ, Pereira-Smith OM, Misteli T: Regulation of Alternative Splicing by Histone Modifications. Science. 2010, 327 (5968): 996-1000. 10.1126/science.1184208.PubMed CentralPubMedView Article
- Xu M, Hecht NB: Polypyrimidine tract-binding protein 2 binds to selective, intronic messenger RNA and microRNA targets in the mouse testis. Biol Reprod. 2011, 84 (3): 435-439. 10.1095/biolreprod.110.087114.PubMedView Article
- Kolasinska-Zwierz P, Down T, Latorre I, Liu T, Liu XS, Ahringer J: Differential chromatin marking of introns and expressed exons by H3K36me3. Nat Genet. 2009, 41 (3): 376-381. 10.1038/ng.322.PubMed CentralPubMedView Article
- Nilsen TW: The spliceosome: the most complex macromolecular machine in the cell?. Bioessays. 2003, 25 (12): 1147-1149. 10.1002/bies.10394.PubMedView Article
- Gregory RI, Yan KP, Amuthan G, Chendrimada T, Doratotaj B, Cooch N, Shiekhattar R: The Microprocessor complex mediates the genesis of microRNAs. Nature. 2004, 432 (7014): 235-240. 10.1038/nature03120.PubMedView Article
- Gregory RI, Chendrimada TP, Shiekhattar R: MicroRNA biogenesis: isolation and characterization of the microprocessor complex. Methods Mol Biol. 2006, 342: 33-47.PubMed
- Smalheiser NR: EST analyses predict the existence of a population of chimeric microRNA precursor-mRNA transcripts expressed in normal human and mouse tissues. Genome Biol. 2003, 4 (7): 403-10.1186/gb-2003-4-7-403.PubMed CentralPubMedView Article
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Research. 2002, 12 (6): 996-1006.PubMed CentralPubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.