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

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 [2][3][4], including oncogenic microRNAs (miRNAs) [5]. MiRNAs are short (~22 bp) noncoding RNAs that are implicated in gene regulation and cancer [6][7][8][9][10]. 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.
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 miR-NAs (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'-TGTGTC- MiR-NAs 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].
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'-GTCGTATCCAGTGCAGGGTC-CGAGGTATTCGCACTGGATACGACTACCTG-3'(mmumir-106a) and 5'-GTCGTATCCAGTGCAGGGTCCGAGG-TATTCGCACTGGATACGACTTACAG-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: 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 Map of integration sites upstream of the mmu-mir-106a cistron Figure 1 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. mmu-mir-363 mmu-mir-92-2 mmu-mir-19b-2 mmu-mir-20b mmu-mir-106a Expression of the primary transcript and mature species of the mmu-mir-106a cistron  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 forma-tion 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.
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.