|
|
|
RESEARCH COMMUNICATION
1 Broad Institute of Massachussetts Institute of Technology and Harvard University, Cambridge, Massachusetts 02141, USA; 2 Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA; 3 European Molecular Biology Laboratory, 69117 Heidelberg, Germany; 4 Department of Biology, Howard Hughes Medical Institute and Whitehead Institute for Biomedical Research, Massachusetts Institute of Technology Cambridge, Massachusetts 02139, USA; 5 Watson School of Biological Sciences and Howard Hughes Medical Institute, Cold Spring Harbor Laboratory, Cold Spring Harbor, New York 11724, USA
| Abstract |
|---|
|
|
|---|
22-nucleotide RNAs that are processed from characteristic precursor hairpins and pair to sites in messages of protein-coding genes to direct post-transcriptional repression. Here, we report that the miRNA iab-4 locus in the Drosophila Hox cluster is transcribed convergently from both DNA strands, giving rise to two distinct functional miRNAs. Both sense and antisense miRNA products target neighboring Hox genes via highly conserved sites, leading to homeotic transformations when ectopically expressed. We also report sense/antisense miRNAs in mouse and find antisense transcripts close to many miRNAs in both flies and mammals, suggesting that additional sense/antisense pairs exist.
[Keywords: Drosophila; miR-iab-4; Hox; antisense miRNAs]]
Received September 6, 2007; revised version accepted November 2, 2007.
22-nucleotide (nt) RNAs that regulate gene expression post-transcriptionally (Bartel 2004| Results and Discussion |
|---|
|
|
|---|
|
To assess the possible biological roles of the two iab-4 miRNAs, we examined fly genes for potential target sites by searching for conserved matches to the seed region of the miRNAs (Lewis et al. 2005
). We found highly conserved target sites for miR-iab-4AS in the 3' untranslated regions (UTRs) of several Hox genes that are proximal to the iab-4 locus and are expressed in the neighboring more anterior embryonic segments: abd-A, Ubx, and Antennapedia (Antp) have four, five, and two seed sites, respectively, most of which are conserved across 12 Drosophila species that diverged 40 million years ago (Fig. 2A; Supplemental Fig. S2; Drosophila 12 Genomes Consortium 2007
; Stark et al. 2007a
). More than two highly conserved sites for one miRNA is exceptional for fly 3' UTRs, placing these messages among the most confidently predicted miRNA targets and suggesting that they might be particularly responsive to the presence of the miRNA. The strong predicted targeting of proximal Hox genes was reminiscent of previously characterized miR-iab-4 targeting of Ubx in flies and miR-196 targeting of HoxB8 in vertebrates (Mansfield et al. 2004
; Yekta et al. 2004
; Hornstein et al. 2005
; Ronshaugen et al. 2005
).
|
If miR-iab-4AS were able to potently down-regulate Ubx in the fly, its misexpression should result in a Ubx loss-of-function phenotype, a line of reasoning that has often been used to study the functions and regulatory relationships of Hox genes. Ubx is expressed throughout the haltere imaginal disc, where it represses wing-specific genes and specifies haltere identity (Weatherbee et al. 1998
). When we expressed mir-iab-4AS in the haltere imaginal disc under bx-Gal4 control, a clear homeotic transformation of halteres to wings was observed (Fig. 3). The halteres developed sense organs characteristic of the wing margin and their size increased severalfold, features typical of transformation to wing (Weatherbee et al. 1998
). Consistent with the increased number of miR-iab4AS target sites, the transformation was stronger than that reported for expression of iab-4 (Ronshaugen et al. 2005
), for which we confirmed changes in morphology but did not find wing-like growth (Fig. 3D).
|
The expression patterns and regulatory connections between Hox genes and the two iab-4 miRNAs show an intriguing pattern in which the miRNAs appear to reinforce Hox gene-mediated transcriptional regulation (Fig. 4A). In particular, miR-iab-4AS would reinforce the posterior expression boundary of abd-A, Ubx, and Antp, supporting their transcriptional repression by Abd-B. mir-iab-4 appears to support abd-A- and Abd-B-mediated repression of Ubx, reinforcing the abd-A/Ubx expression domains and the posterior boundary of Ubx expression. Furthermore, both iab-4 miRNAs have conserved target sites in Antp, which is also repressed by Abd-B, abd-A, and Ubx. The iab-4 miRNAs thus appear to support the established regulatory hierarchy among Hox transcription factors, which exhibits "posterior prevalence," in that more posterior Hox genes repress more anterior ones and are dominant in specifying segment identity (for reviews, see McGinnis and Krumlauf 1992
; Pearson et al. 2005
). Interestingly, Abd-B and mir-iab-4AS are expressed in the same segments, and the majority of cis-regulatory elements controlling Abd-B expression are located 3' of Abd-B (Boulet et al. 1991
). This places them near the inferred transcription start of mir-iab-4AS, where they potentially direct the coexpression of these genes. Similarly, abd-A and mir-iab-4 may be coregulated as both are transcribed divergently, potentially under the control of shared upstream elements.
|
The genomic arrangement of two miRNAs that are expressed from the same locus but on different strands might provide a simple and efficient means to create nonoverlapping miRNA expression domains (Fig. 4B). Such sense/antisense miRNAs could restrict each others transcription, either by direct transcriptional interference, as shown for overlapping convergently transcribed genes (Shearwin et al. 2005
; Hongay et al. 2006
), or post-transcriptionally, possibly via RNA–RNA duplexes formed by the complementary transcripts. Sense/antisense miRNAs would usually differ at their 5' ends and thereby target distinct sets of genes, which might help define and establish sharp boundaries between expression domains. Coupled with feedback loops or coregulation of miRNAs and genes in cis or trans, this arrangement could provide a powerful regulatory switch. The iab-4 miRNAs might be a special case of tight regulatory integration in which miRNAs and proximal genes appear coregulated transcriptionally in cis and repress each other both transcriptionally and post-transcriptionally.
It is perhaps surprising that no antisense miRNA had been found previously, even though, for example, the intriguing expression pattern of the iab-4 transcripts had been reported nearly two decades ago (Cumberledge et al. 1990
; Bae et al. 2002
), and iab-4 lies in one of the most extensively studied regions of the Drosophila genome. The frequent occurrence of antisense transcripts (Yelin et al. 2003
; Katayama et al. 2005
) suggests that more antisense miRNAs might exist. Indeed, up to 13% of known Drosophila, 20% of mouse, and 31% of human miRNAs are located in introns of host genes transcribed on the opposite strand or are within 50 nt of antisense ESTs or cDNAs (Supplemental Table S1). These include an antisense transcript overlapping human mir-196 (see also Mainguy et al. 2007
). However, because of the contribution of noncanonical base pairs, particularly G:U pairs that become less favorable A:C in the antisense strand, many miRNA antisense transcripts will not fold into hairpin structures suitable for miRNA biogenesis, which explains the propensity of miRNA gene predictions to identify the correct strand (Lim et al. 2003
). Nonetheless, in a recent prediction effort, 22 sequences reverse-complementary to known Drosophila miRNAs showed scores seemingly compatible with miRNA processing (Stark et al. 2007b
). Deep sequencing of small RNA libraries from Drosophila confirmed the processing of small RNAs from four of these high-scoring antisense candidates (Ruby et al. 2007
), and the ovary/testes libraries used here showed antisense reads for an additional Drosophila miRNA (mir-312) (see Supplemental Tables S2, S3). In addition, using high-throughput sequencing of small RNA libraries from mice, we found sequencing reads that uniquely matched the mouse genome in loci antisense to 10 annotated mouse miRNAs. Eight of the inferred antisense miRNAs were supported by multiple independent reads, and two of them had reads from both the mature miRNA and the star sequence (Supplemental Table S2). These results suggest that sense/antisense miRNAs could be more generally employed in diverse contexts and in species as divergent as flies and mammals.
| Materials and methods |
|---|
|
|
|---|
3' UTRs were amplified from Drosophila melanogaster genomic DNA and cloned in pCR2.1 for site-directed mutagenesis. The following primer pairs were used to amplify the indicated 3' UTR: abd-A (tc tagaGCGGTCAGCAAAGTCAACTC; gtcgacATGGATGGGTTCTCGT TGCAG), Ubx (tctagaATCCTTAGATCCTTAGATCCTTAG; ctcgag ATGGTTTGAATTTCCACTGA), and Abd-B (tctagaGCCACCACCT GAACCTTAG; aactcgagCGGAGTAATGCGAAGTAATTG). QuickChange multisite-directed mutagenesis was used to mutate all miR-iab-4AS seed sites from ATACGT to ATAGGT, per the manufacturers directions (Stratagene). Wild-type and mutated 3' UTRs were subcloned into pCJ40 between SacI and NotI sites to make Renilla luciferase reporters. Plasmid pCJ71 contains the abd-A wild-type 3' UTR, pCJ72 contains the Ubx wild-type 3' UTR, pCJ74 contains the Abd-B wild-type 3' UTR, pCJ75 contains the abd-A mutated 3' UTR, and pCJ76 contains the Ubx mutated 3' UTR fused to Renilla luciferase. The control let-7 expression vector was obtained by amplifying let-7 from genomic DNA with primers 474 base pairs (bp) upstream of and 310 bp downstream from the let-7 hairpin and cloning it into pMT-puro. To express miR-iab-4 and miR-iab-4AS, a 430-bp genomic fragment containing the miR-iab-4 hairpin was cloned, in either direction, downstream from the tubulin promoter as described in Stark et al. (2005)
. For the UAS-miR-iab-4 and UAS-miR-iab-4AS constructs, the same 430-bp genomic fragment containing the miR-iab-4 hairpin was cloned downstream from pUAST-DSred2 (Stark et al. 2003
) in either direction.
Reporter assays
For the luciferase assays, 2 ng of p2129 (firefly luciferase), 4 ng of Renilla reporter, 48 ng of miRNA expression plasmid, and 48 ng of p2032 (GFP) were cotransfected with 0.3 µL Fugene HD per well of a 96-well plate. Twenty-four hours after transfection, expression of Renilla luciferase was induced by addition of 500 µM CuSO4 to the culture media. Twenty-four hours after induction, reporter activity was measured with the Dual-Glo luciferase kit (Promega), per the manufacturers instructions on a Tecan Safire II plate reader.
The ratio of Renilla:firefly luciferase activity was measured for each well. To calculate fold repression, the ratio of Renilla:firefly for reporters cotransfected with let-7 was set to 1. The Wilcoxon rank-sum test was used to assess the significance of changes in fold repression of wild-type reporters compared with mutant reporters. Geometric means from 16 transfections representing four replicates of four independent transfections are shown. Error bars represent the fourth highest and lowest values of each set.
Drosophila strains
UAS-miR-iab-4 and UAS-miR-iab-4AS flies were generated by injection of the corresponding plasmids into w1118 embryos. bxMS1096-GAL4 flies were obtained from the Bloomington Stock Center.
In situ hybridization and protein stainings
Double in situ hybridization for the miRNA primary transcripts was performed as described in Stark et al. (2005)
. Probes were generated using PCR on genomic DNA with primers TCAGAGCATGCAGAGACAT AAAG, TTGTAGATTGAAATCGGACACG for iab-4 sense and ATTT TACTGGGTGTCTGGGAAAG, TAGAAACTGAGACGGAGAAGCAG for iab-4 antisense. Protein stainings were performed as described in Patel (1994)
. Antibodies used were mouse anti-Ubx (1:30), mouse anti-abd-A (1:5), and HRP-conjugated goat anti-mouse (Dianova, 1:3000).
RT–PCRs
Total RNA was isolated using Trizol (Invitrogen), treated with RQI DNase (Promega), and used for strand-specific cDNA synthesis with SuperScript III (Invitrogen). Primers for cDNA synthesis were CATATAA CAAAGTGCTACGTG (iab-4 sense) and CTTTATCTGCATTTG GATCCG (iab-4 antisense). Both primers were used for subsequent amplification.
Small library sequencing
Drosophila small RNAs were cloned from adult ovaries and testes as described previously (Brennecke et al. 2007
) and sequenced using Solexa sequencing. A total of 657,251 sequencing reads uniquely matched known Drosophila miRNAs (Rfam release 9.2), and the 69 miRNAs with unique matches had 1011 matches on average (Stark et al. 2007b
). Two miRNAs had unique matches to the antisense hairpin (Supplemental Tables S2, S3). Mouse small RNAs were cloned from wild-type and c-kit mutant ovaries (Supplemental Table S4; G. Hannon, pers. comm.) and from Comma-Dbgeo cells, a murine mammary epithelial cell line (Ibarra et al. 2007
), and were sequenced using Solexa sequencing. A total of 4,217,883 reads uniquely matched known mouse miRNAs (Rfam release 9.2), and the 286 miRNAs with unique reads showed 256 reads on average. Sequencing reads matching to the plus and minus strand of known mouse miRNAs with antisense reads are listed in Supplemental Table S3.
Multiple sequence alignments and target site prediction
The multiple sequence alignments for the indicated Hox 3' UTRs were obtained from the University of California at Santa Cruz (UCSC) genome browser (Kent et al. 2002
) and were slightly manually adjusted. We predicted target sites according to Lewis et al. (2005)
by searching for 3' UTR seed sites (reverse-complementary to miRNA positions 2–8 or matching to "A" + reverse complement of miRNA positions 2–7).
Antisense transcripts near known miRNAs
To assess the fraction of Drosophila, human, and mouse miRNAs that are also putatively transcribed on both strands and might give rise to antisense miRNAs, we determined the number of miRNAs that are near known transcripts on the opposite strand. We obtained the coordinates of all introns of protein-coding genes and all mapped ESTs or cDNAs for the three species from the UCSC genome browser (Kent et al. 2002
). We intersected them with the miRNA coordinates from Rfam (release 9.2; Griffiths-Jones et al. 2006
), requiring miRNAs and transcripts to be on opposite strands and at a distance of at most 50 nt. For each miRNA, we recorded the number of antisense transcripts and their identifiers. Note that some of the transcripts might have been mapped to more than one place in the genome, such that the intersection represents an upper estimate based on the currently known transcripts.
| Acknowledgments |
|---|
|
|
|---|
| Footnotes |
|---|
7 Present address: Temasek Life Sciences Laboratory, The National University of Singapore, Singapore 117604. ![]()
E-MAIL alex.stark{at}mit.edu; FAX (617) 253-7512. ![]()
9 E-MAIL manoli{at}mit.edu; FAX (617) 253-7512. ![]()
Supplemental material is available at http://www.genesdev.org.
Article is online at http://www.genesdev.org/cgi/doi/10.1101/gad.1613108
| References |
|---|
|
|
|---|
Bartel, D.P. 2004. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 116: 281–297.[CrossRef][Medline]
Bender, W. 2008. MicroRNAs in the Drosophila bithorax complex. Genes & Dev. (this issue) doi: 10.1101/gad.1614208.
Boulet, A.M., Lloyd, A., and Sakonju, S. 1991. Molecular definition of the morphogenetic and regulatory functions and the cis-regulatory elements of the Drosophila Abd-B homeotic gene. Development 111: 393–405.[Abstract]
Brennecke, J., Stark, A., Russell, R.B., and Cohen, S.M. 2005. Principles of microRNA-target recognition. PLoS Biol. 3: e85. doi: 10.1371/journal.pbio.0030085.[CrossRef][Medline]
Brennecke, J., Aravin, A.A., Stark, A., Dus, M., Kellis, M., Sachidanandam, R., and Hannon, G.J. 2007. Discrete small RNA-generating loci as master regulators of transposon activity in Drosophila. Cell 128: 1089–1103.[CrossRef][Medline]
Cumberledge, S., Zaratzian, A., and Sakonju, S. 1990. Characterization of two RNAs transcribed from the cis-regulatory region of the abd-A domain within the Drosophila bithorax complex. Proc. Natl. Acad. Sci. 87: 3259–3263.
Drosophila 12 Genomes Consortium. 2007. Evolution of genes and genomes on the Drosophila phylogeny. Nature 450: 203–218.[CrossRef][Medline]
Duboule, D. 1998. Vertebrate hox gene regulation: Clustering and/or colinearity? Curr. Opin. Genet. Dev. 8: 514–518.[CrossRef][Medline]
Farh, K.K., Grimson, A., Jan, C., Lewis, B.P., Johnston, W.K., Lim, L.P., Burge, C.B., and Bartel, D.P. 2005. The widespread impact of mammalian microRNAs on mRNA repression and evolution. Science 310: 1817–1821.
Griffiths-Jones, S., Grocock, R.J., van Dongen, S., Bateman, A., and Enright, A.J. 2006. miRBase: MicroRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 34 (Database issue): D140–D144. doi: 10.1093/nar/gkj112.
Hongay, C.F., Grisafi, P.L., Galitski, T., and Fink, G.R. 2006. Antisense transcription controls cell fate in Saccharomyces cerevisiae. Cell 127: 735–745.[CrossRef][Medline]
Hornstein, E., Mansfield, J.H., Yekta, S., Hu, J.K., Harfe, B.D., McManus, M.T., Baskerville, S., Bartel, D.P., and Tabin, C.J. 2005. The microRNA miR-196 acts upstream of Hoxb8 and Shh in limb development. Nature 438: 671–674.[CrossRef][Medline]
Ibarra, I., Erlich, Y., Muthuswamy, S.K., Sachidanandam, R., and Hannon, G.J. 2007. A microRNA fingerprint of mammary epithelial stem cells. Genes & Dev. 21: 3238–3243.
Karch, F., Bender, W., and Weiffenbach, B. 1990. abdA expression in Drosophila embryos. Genes & Dev. 4: 1573–1587.
Katayama, S., Tomaru, Y., Kasukawa, T., Waki, K., Nakanishi, M., Nakamura, M., Nishida, H., Yap, C.C., Suzuki, M., Kawai, J., et al. 2005. Antisense transcription in the mammalian transcriptome. Science 309: 1564–1566.
Kent, W.J., Sugnet, C.W., Furey, T.S., Roskin, K.M., Pringle, T.H., Zahler, A.M., and Haussler, D. 2002. The human genome browser at UCSC. Genome Res. 12: 996–1006.
Khvorova, A., Reynolds, A., and Jayasena, S.D. 2003. Functional siRNAs and miRNAs exhibit strand bias. Cell 115: 209–216.[CrossRef][Medline]
Lewis, E.B. 1978. A gene complex controlling segmentation in Drosophila. Nature 276: 565–570.[CrossRef][Medline]
Lewis, B.P., Burge, C.B., and Bartel, D.P. 2005. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120: 15–20.[CrossRef][Medline]
Lim, L.P., Lau, N.C., Weinstein, E.G., Abdelhakim, A., Yekta, S., Rhoades, M.W., Burge, C.B., and Bartel, D.P. 2003. The microRNAs of Caenorhabditis elegans. Genes & Dev. 17: 991–1008.
Macias, A., Casanova, J., and Morata, G. 1990. Expression and regulation of the abd-A gene of Drosophila. Development 110: 1197–1207.
Mainguy, G., Koster, J., Woltering, J., Jansen, H., and Durston, A. 2007. Extensive polycistronism and antisense transcription in the mammalian hox clusters. PLoS ONE 2: e356. doi: 10.1371/journal.pone.0000356.[CrossRef]
Mansfield, J.H., Harfe, B.D., Nissen, R., Obenauer, J., Srineel, J., Chaudhuri, A., Farzan-Kashani, R., Zuker, M., Pasquinelli, A.E., Ruvkun, G., et al. 2004. MicroRNA-responsive sensor transgenes uncover Hox-like and other developmentally regulated patterns of vertebrate microRNA expression. Nat. Genet. 36: 1079–1083.[CrossRef][Medline]
McGinnis, W. and Krumlauf, R. 1992. Homeobox genes and axial patterning. Cell 68: 283–302.[CrossRef][Medline]
Patel, N.H. 1994. Imaging neuronal subsets and other cell types in whole-mount Drosophila embryos and larvae using antibody probes. Methods Cell Biol. 44: 445–487.[Medline]
Pearson, J.C., Lemons, D., and McGinnis, W. 2005. Modulating Hox gene functions during animal body patterning. Nat. Rev. Genet. 6: 893–904.[CrossRef][Medline]
Ronshaugen, M., Biemar, F., Piel, J., Levine, M., and Lai, E.C. 2005. The Drosophila microRNA iab-4 causes a dominant homeotic transformation of halteres to wings. Genes & Dev. 19: 2947–2952.
Ruby, J.G., Stark, A., Johnston, W.K., Kellis, M., Bartel, D.P., and Lai, E.C. 2007. Evolution, biogenesis, expression, and target predictions of a substantially expanded set of Drosophila microRNAs. Genome Res. doi: 10.1101/gr.6597907.
Schwarz, D.S., Hutvagner, G., Du, T., Xu, Z., Aronin, N., and Zamore, P.D. 2003. Asymmetry in the assembly of the RNAi enzyme complex. Cell 115: 199–208.[CrossRef][Medline]
Shearwin, K.E., Callen, B.P., and Egan, J.B. 2005. Transcriptional interference—A crash course. Trends Genet. 21: 339–345.[CrossRef][Medline]
Stark, A., Brennecke, J., Russell, R.B., and Cohen, S.M. 2003. Identification of Drosophila microRNA targets. PLoS Biol. 1: E60. doi: 10.1371/journal.pbio.0000060.[CrossRef][Medline]
Stark, A., Brennecke, J., Bushati, N., Russell, R.B., and Cohen, S.M. 2005. Animal microRNAs confer robustness to gene expression and have a significant impact on 3'UTR evolution. Cell 123: 1133–1146.[CrossRef][Medline]
Stark, A., Lin, M.F., Kheradpour, P., Pedersen, J.S., Parts, L., Carlson, J.W., Crosby, M.A., Rasmussen, M.D., Roy, S., Deoras, A.N., et al. 2007a. Discovery of functional elements in 12 Drosophila genomes using evolutionary signatures. Nature 450: 219–232.[CrossRef][Medline]
Stark, A., Kheradpour, P., Parts, L., Brennecke, J., Hodges, E., Hannon, G.J., and Kellis, M. 2007b. Systematic discovery and characterization of fly microRNAs using 12 Drosophila genomes. Genome Res. doi: 10.1101/gr.6593807.
Weatherbee, S.D., Halder, G., Kim, J., Hudson, A., and Carroll, S. 1998. Ultrabithorax regulates genes at several levels of the wing-patterning hierarchy to shape the development of the Drosophila haltere. Genes & Dev. 12: 1474–1482.
Yekta, S., Shih, I.H., and Bartel, D.P. 2004. MicroRNA-directed cleavage of HOXB8 mRNA. Science 304: 594–596.
Yelin, R., Dahary, D., Sorek, R., Levanon, E.Y., Goldstein, O., Shoshan, A., Diber, A., Biton, S., Tamir, Y., Khosravi, R., et al. 2003. Widespread occurrence of antisense transcription in the human genome. Nat. Biotechnol. 21: 379–386.[CrossRef][Medline]
Yoder, J.H. and Carroll, S.B. 2006. The evolution of abdominal reduction and the recent origin of distinct Abdominal-B transcript classes in Diptera. Evol. Dev. 8: 241–251.[CrossRef][Medline]
Zuker, M. 2003. Mfold Web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 31: 3406–3415.
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:
![]() |
W. Bender MicroRNAs in the Drosophila bithorax complex Genes & Dev., January 1, 2008; 22(1): 14 - 19. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Tyler, K. Okamura, W.-J. Chung, J. W. Hagen, E. Berezikov, G. J. Hannon, and E. C. Lai Functionally distinct regulatory RNAs generated by bidirectional transcription and processing of microRNA loci Genes & Dev., January 1, 2008; 22(1): 26 - 36. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||