Skip to main content

Advertisement

Retroviral activation of the mir-106a microRNA cistron in T lymphoma

Abstract

Retroviral insertion into a host genome is a powerful tool not only for the discovery of cancer genes, but also for the discovery of potential oncogenic noncoding RNAs. In a large-scale mouse T lymphocyte tumor screen we found a high density of integrations upstream of the mir-106a microRNA cistron. In tumors containing an integration, the primary transcript encoding the mir-106a cistron was overexpressed five to 20-fold compared with that of control tumors; concomitantly, the mature mir-106a and mir-363 microRNAs were highly overexpressed as well. These findings suggest the mir-106a cistron plays an important role in T cell tumorigenesis.

Findings

Retroviral insertions into the genome of a host can induce tumor formation by altering gene expression or function. Integration of a retrovirus near a gene can induce overexpression of the gene through the viral promoter or enhancer, while insertion of a retrovirus into a gene can cause both activation and inactivation. If the affected genes are proto-oncogenes or tumor suppressor genes, the insertion events may lead to tumor formation [1]. Consequently, retroviral mutagenesis has been used to search entire genomes for genes involved in cancer development [24], including oncogenic microRNAs (miRNAs) [5]. MiRNAs are short (~22 bp) noncoding RNAs that are implicated in gene regulation and cancer [610]. In a large-scale retroviral insertion mutagenesis screen, we used the murine leukemia virus (MLV) strain SL3-3, which causes T lymphomas [11], and identified several miRNAs that are potentially involved in tumorigenesis. We previously demonstrated that a group of these retroviral insertions induces overexpression of the oncogenic mmu-mir-17 miRNA cistron in mouse tumors [5]. Here we build on our validation of the retrovirus insertional mutagenesis method to identify oncogenic miRNA and present another potentially oncogenic miRNA cistron, mmu-mir-106a. In this screen, male BALB/c mice were treated with ethyl-nitroso-urea (ENU) and bred to normal female mice. ENU treatment was conducted to increase the recovery of tumor suppressors in the F1 progeny through mutagenesis of the paternal allele. Newborn offspring mice were then injected with MLV strain SL3-3. After becoming moribund due to tumor development, mice were euthanized and thymus and spleen tissues were collected and stored at -80°C. Locations of the SL3-3 provirus integration sites were identified as previously described using a splinkerette based PCR method [3] that amplifies genomic DNA flanking the 5' LTR of the virus.

We identified 6234 integration sites in 2199 tumors; of these tumors, 76 sites were located on chromosome X upstream of a miRNA cluster containing mmu-mir-106a, mmu-mir-20b, mmu-mir-19b-2, mmu-mir-92-2, and mmu-mir-363. The locations of the integrations ranged from 1.5 kb to 22 kb upstream of the miRNA cluster (Fig. 1), with proviral inserts in both sense and anti-sense orientations with respect to the primary RNA transcript encoding the miRNA cistron. The Mouse Retroviral Tagged Cancer Gene Database [12], which compiles retroviral insertions into the genomic DNA from various non-T cell derived mouse tumors, also lists 10 integrations located upstream of the mmu-mir-106a cluster. Furthermore, Hwang et al. found that EST AI464896, which maps to the same location as mmu-mir-363, was overexpressed in tumors with proviral MLV integrations into this region [13]. The radiation leukemia virus (RadLV) also frequently integrates at this locus and a group of five differentially spliced noncoding RNAs known as Kis2 (GenBank Accession numbers AY940614-AY940618) are overexpressed in these tumors [14]. Because the Kis2 transcripts lie directly upstream of the mir-106a miRNA cluster (mmu-mir-106a overlaps these transcripts by four bases), they likely are part of the primary transcripts containing the miRNA cluster.

Figure 1
figure1

Map of integration sites upstream of the mmu-mir-106a cistron. A map of the SL3-3 retroviral sites upstream of the mir-106a cluster is displayed using the UCSC genome website browser (February 2006 version of the mm8 genome assembly). Insertion sites are depicted as vertical handlebars. Tumors assayed by quantitative PCR are numbered and noted in black text.

To determine whether the retroviral integrations in this region affected the expression of the mir-106a cistron, we used quantitative PCR (qPCR) to measure expression levels of the primary transcript (Kis2) and the mature miRNAs (mmu-mir-106a and mmu-mir-363) in tumors containing mir-106a cistron integrations as well as in control tumors lacking such integrations. To measure primary transcript (Kis2) expression levels, a probe and primer set was designed to AY940616, which is a common exon to three of the alternatively spliced forms of Kis2. The probe and primers for AY940616 were as follows: 5'-TGTGTCCCTGAAGTTTATTGGTGT-3', 5'-GGGTCACGAGCTCCCTCC-3', and 5'-[6-FAM]-CCCCCATCAACACAAACATTCCATCA-[3BHQ1]-3'. MiRNAs and low molecular weight RNAs were isolated from frozen mouse tumor tissue using the Purelink miRNA Isolation Kit (Invitrogen). Large fraction RNAs were then purified by eluting the high molecular weight RNA bound to the first column (used for the miRNA purification). cDNA was generated from total RNA by reverse transcription with random hexamers using the SuperScript First-Strand Synthesis System for RT-PCR (Invitrogen). qPCR runs were conducted on the MX3000P (Stratagene). All qPCR reactions were run in triplicate. As controls, tumors not containing integrations near the mmu-mir-106a-363 cluster were also assayed. Beta-actin was used as the endogenous reference gene (Mouse ACTB 20× VIC-MGB probe set, Applied Biosystems) and control tumor 1 was used as the calibrator sample in the calculation of 2-ΔΔCt values (relative expression). All relative expression values were normalized such that the average of the tumor controls was set to 1.

Representative tumors with integration sites spanning the upstream region of mir-106a were measured for expression of the miRNA primary transcript (Fig. 1 and Table 1). In 16 of the 21 tumors assayed, expression of AY940616 was elevated five to 20 fold as compared to the average expression of tumors with no integrations at this locus (Fig. 2A). This confirms the previous report that proviral integrations in this region can increase expression of the Kis2 locus [14].

Figure 2
figure2

Expression of the primary transcript and mature species of the mmu-mir-106a cistron. Quantitative PCR data for tumors with integrations upstream of the mir-106a cistron. (A) Relative expression of AY940616 in tumors containing integration sites near the mir-106a miRNA cluster. Control tumors contain integration sites at locations in the genome other than the mir-106a region. "N" is cDNA generated from normal mouse spleen RNA (Ambion). (B) Relative expression of the mature species of mmu-mir-106a and mmu-mir-363 in tumors containing integration sites near the mir-106a cistron. Tumors are numbered as in Figure 1.

Table 1 Assayed mmu-mir-106a cistron integrations

The mature species of mmu-mir-106a and mmu-mir-363 were then measured by RT-qPCR using a stem-loop RT primer specific for each miRNA [15]. Accordingly, 50 ng of each tumor miRNA preparation was reverse transcribed with the SuperScript First-Strand Synthesis System for RT-PCR using the following stem loop RT primers (50 nM final concentration) 5'-GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTACCTG-3'(mmu-mir-106a) and 5'-GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGACTTACAG-3' (mmu-mir-363). The reverse transcription reactions were diluted 1:200 and 5 μl of these dilutions were used in the 25 μl qPCR reactions. The annealing step was 50°C for 60s. The qPCR probes and primers were as follows: mmu-mir-106a: 5'-CGGCAAAGTGCTAACAGT-3', 5'-GTGCAGGGTCCGAGGT-3', 5'- [6-FAM]- CACTGGATACGACTACCTGC- [BHQ1]-3'; and mmu-mir-363: 5'-TGCGGATTGCACGGTATC-3', 5'-GTGCAGGGTCCGAGGT-3', 5'- [6-FAM]- CACTGGATACGACTTACAGATG- [BHQ1]-3'. Synthetic RNA oligos (IDT) were used to generate a calibration curve for each miRNA: 5'-CAAAGUGCUAACAGUGCAGGUA-3' (mmu-mir-106a) and 5'-AUUGCACGGUAUCCAUCUGUAA-3'(mmu-mir-363). Amplification efficiencies of the calibration curves for mmu-mir-106a and mmu-mir-363 were respectively 67% and 69%. Concentrations of the mature species were calculated using the calibration curves and then normalized by the average of the control tumors, to calculate relative expression levels.

Fifteen tumors with integrations in this region were assayed by qPCR for the mature species of mmu-mir-106a and mmu-mir-363. Approximately 70% of these tumors had increased expression levels of mmu-mir-106a by two to six fold, and of mmu-mir-363 by four to 12 fold over the average expression of tumors with no integrations in this region (Fig. 2B). The mature miRNA expression difference between tumors with integrations in this region and the tumor controls was statistically significant [p < 0.00001 (mmu-mir-106a) and p < 0.0001 (mmu-mir-363)] by a two sample unequal variance Student's t test. From these data we conclude that retroviral integrations in the Kis2 region cause overexpression not only of the primary RNA, but also of the mature species of the mir-106a cluster. This, in turn, suggests that the miRNA cluster can drive the development of T lymphomas. Although there is a possibility that these integrations also may affect the expression of other oncogenes and tumor suppressors in this region, our data clearly indicates a majority of these integrations induce the expression of the mir-106a cluster.

As the mir-106a cistron is a homolog of the oncogenic mir-17 cistron [16], it is not unexpected that mir-106a would also be involved tumorigenesis. Indeed, in human solid tumors, mir-106a expression is increased in colon, pancreas, and prostate tumors; and mir-92-2 expression is increased in pancreas, prostate, and stomach tumors [17]. Given the sequence similarity between the mir-17 and mir-106a cistrons, it is likely that these clusters have overlapping gene targets. In humans, the mir-106a cistron contains several paralogs to members of the mir-17 cistron including mir-17, mir-19b-1, and mir-92-1 [16], which are implicated in cancer development: overexpression of the mir-17 cluster accelerates lymphoma formation from cells of mice overexpressing c-Myc [7]. The mir-17 cluster is also overexpressed in human lung cancer [18]. However, in breast cancer cells, mir-17-5p expression is decreased; there it acts as a translational repressor of the oncogene AIB1 (amplified in breast cancer 1) [19], and in this context may formally act as a tumor suppressor.

It is well established that tumorigenesis is the result of accumulating several cooperating mutations that drive relentless proliferation and aid in metastases. Viral insertional mutagenesis, though perhaps not providing all the mutations necessary for a full-blown tumor, follows this multistep scenario. Although in general the superinfection barrier largely prevents multiple proviral integrations within the same cell, re-infection does happen over time. Because it is a rare event, such cells are selected over the others only when these integrations also give a growth advantage. As a consequence, in general, most viral insertions ("co-mutations") in a single tumor are thought to be causative in its formation. With the caveats of potential passenger genes and potential oligoclonality of tumors, co-mutation analysis may be a powerful way to find cooperating signaling pathways in tumorigenesis.

We detected multiple insertion sites in all of the tumor samples we assayed from the mir-106a cluster. Genes near common co-integration sites for these tumors include Ahi1, Evi5, and Gfi1, candidates previously appearing in retroviral screens [12], as well as PVT1, a noncoding RNA frequently amplified with myc [20]. A summary of all integration sites in the assayed tumors is listed in Table 2.

Table 2 Summary of integrations in tumors assayed for the mmu-mir-106a cistron

Through retroviral insertion in the mouse, we have discovered another potentially oncogenic microRNA cluster, mir-106a-363. Retroviral insertion caused significant overexpression of this microRNA cluster indicating its role in tumor development. This study further demonstrates the power of retrovirus insertion as a tool to discover new oncogenic noncoding RNAs.

References

  1. 1.

    Uren AG, Kool J, Berns A, van Lohuizen M: Retroviral insertional mutagenesis: past, present and future. Oncogene. 2005, 24: 7656-7672. 10.1038/sj.onc.1209043.

  2. 2.

    Suzuki T, Shen H, Akagi K, Morse HC, Malley JD, Naiman DQ, Jenkins NA, Copeland NG: New genes involved in cancer identified by retroviral tagging. Nat Genet. 2002, 32: 166-174. 10.1038/ng949.

  3. 3.

    Mikkers H, Allen J, Knipscheer P, Romeijn L, Hart A, Vink E, Berns A: High-throughput retroviral tagging to identify components of specific signaling pathways in cancer. Nat Genet. 2002, 32: 153-159. 10.1038/ng950.

  4. 4.

    Li J, Shen H, Himmel KL, Dupuy AJ, Largaespada DA, Nakamura T, Shaughnessy JD, Jenkins NA, Copeland NG: Leukaemia disease genes: large-scale cloning and pathway predictions. Nat Genet. 1999, 23: 348-353. 10.1038/15531.

  5. 5.

    Wang CL, Wang BB, Bartha G, Li L, Channa N, Klinger M, Killeen N, Wabl M: Activation of an oncogenic microRNA cistron by provirus integration. Proc Natl Acad Sci U S A. 2006, 103: 18680-18684. 10.1073/pnas.0609030103.

  6. 6.

    Shamovsky I, Ivannikov M, Kandel ES, Gershon D, Nudler E: RNA-mediated response to heat shock in mammalian cells. Nature. 2006, 440: 556-560. 10.1038/nature04518.

  7. 7.

    He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435: 828-833. 10.1038/nature03552.

  8. 8.

    Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, Shimizu M, Rattan S, Bullrich F, Negrini M, Croce CM: Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A. 2004, 101: 2999-3004. 10.1073/pnas.0307323101.

  9. 9.

    O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435: 839-843. 10.1038/nature03677.

  10. 10.

    Costinean S, Zanesi N, Pekarsky Y, Tili E, Volinia S, Heerema N, Croce CM: Pre-B cell proliferation and lymphoblastic leukemia/high-grade lymphoma in E(mu)-miR155 transgenic mice. Proc Natl Acad Sci U S A. 2006, 103: 7024-7029. 10.1073/pnas.0602266103.

  11. 11.

    Lenz J, Crowther R, Klimenko S, Haseltine W: Molecular cloning of a highly leukemogenic, ecotropic retrovirus from an AKR mouse. J Virol. 1982, 43: 943-951.

  12. 12.

    Akagi K, Suzuki T, Stephens RM, Jenkins NA, Copeland NG: RTCGD: retroviral tagged cancer gene database. Nucleic Acids Res. 2004, 32: D523-7. 10.1093/nar/gkh013.

  13. 13.

    Hwang HC, Martins CP, Bronkhorst Y, Randel E, Berns A, Fero M, Clurman BE: Identification of oncogenes collaborating with p27Kip1 loss by insertional mutagenesis and high-throughput insertion site analysis. Proc Natl Acad Sci U S A. 2002, 99: 11293-11298. 10.1073/pnas.162356099.

  14. 14.

    Landais S, Quantin R, Rassart E: Radiation leukemia virus common integration at the Kis2 locus: simultaneous overexpression of a novel noncoding RNA and of the proximal Phf6 gene. J Virol. 2005, 79: 11443-11456. 10.1128/JVI.79.17.11443-11456.2005.

  15. 15.

    Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-10.1093/nar/gni178.

  16. 16.

    Tanzer A, Stadler PF: Molecular evolution of a microRNA cluster. J Mol Biol. 2004, 339: 327-335. 10.1016/j.jmb.2004.03.065.

  17. 17.

    Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006, 103: 2257-2261. 10.1073/pnas.0510565103.

  18. 18.

    Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005, 65: 9628-9632. 10.1158/0008-5472.CAN-05-2352.

  19. 19.

    Hossain A, Kuo MT, Saunders GF: Mir-17-5p Regulates Breast Cancer Cell Proliferation by Inhibiting Translation of AIB1 mRNA. Mol Cell Biol. 2006, 26: 8191-8201. 10.1128/MCB.00242-06.

  20. 20.

    Shtivelman E, Bishop JM: The PVT gene frequently amplifies with MYC in tumor cells. Mol Cell Biol. 1989, 9: 1148-1154.

Download references

Acknowledgements

This work was supported by NIH grant CA100266 to MW, and Synergenics, LLC. We thank Dr. Clifford Wang for his technical advice and for his comments on the manuscript.

Author information

Correspondence to Matthias Wabl.

Additional information

Competing interests

The authors declare a financial interest in Picobella, LLC.

Authors' contributions

AML carried out the RNA isolation, quantitative PCR, expression data analysis, and drafted the manuscript. GB, LL, NC, and BBW carried out the tag recovery and identification. BBW and MW planned and directed the execution of the retroviral screen, the design of the study, and the writing of the manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article

Keywords

  • Murine Leukemia Virus
  • Primary Transcript
  • miRNA Cluster
  • Mature Species
  • microRNA Cluster