- Open Access
Novel principles of gamma-retroviral insertional transcription activation in murine leukemia virus-induced end-stage tumors
Retrovirology volume 11, Article number: 36 (2014)
Insertional mutagenesis screens of retrovirus-induced mouse tumors have proven valuable in human cancer research and for understanding adverse effects of retroviral-based gene therapies. In previous studies, the assignment of mouse genes to individual retroviral integration sites has been based on close proximity and expression patterns of annotated genes at target positions in the genome. We here employed next-generation RNA sequencing to map retroviral-mouse chimeric junctions genome-wide, and to identify local patterns of transcription activation in T-lymphomas induced by the murine leukemia gamma-retrovirus SL3-3. Moreover, to determine epigenetic integration preferences underlying long-range gene activation by retroviruses, the colocalization propensity with common epigenetic enhancer markers (H3K4Me1 and H3K27Ac) of 6,117 integrations derived from end-stage tumors of more than 2,000 mice was examined.
We detected several novel mechanisms of retroviral insertional mutagenesis: bidirectional activation of mouse transcripts on opposite sides of a provirus including transcription of unannotated mouse sequence; sense/antisense-type activation of genes located on opposite DNA strands; tandem-type activation of distal genes that are positioned adjacently on the same DNA strand; activation of genes that are not the direct integration targets; combination-type insertional mutagenesis, in which enhancer activation, alternative chimeric splicing and retroviral promoter insertion are induced by a single retrovirus. We also show that irrespective of the distance to transcription start sites, the far majority of retroviruses in end-stage tumors colocalize with H3K4Me1 and H3K27Ac-enriched regions in murine lymphoid tissues.
We expose novel retrovirus-induced host transcription activation patterns that reach beyond a single and nearest annotated gene target. Awareness of this previously undescribed layer of complexity may prove important for elucidation of adverse effects in retroviral-based gene therapies. We also show that wild-type gamma-retroviruses are frequently positioned at enhancers, suggesting that integration into regulatory regions is specific and also subject to positive selection for sustaining long-range gene activation in end-stage tumors. Altogether, this study should prove useful for extrapolating adverse outcomes of retroviral vector therapies, and for understanding fundamental cellular regulatory principles and retroviral biology.
The murine leukemia viruses (MLVs) are slowly-transforming gamma-retroviruses that induce tumors of hematopoietic origin upon injection into susceptible newborn mice [1, 2]. MLVs are simple retroviruses that contain only the gag, pol and env genes flanked at both ends by major regulatory regions, the long terminal repeats (LTRs) composed of U3, R and U5. A hallmark of retroviral replication is the stable integration of the ~9 kb genome into host chromosomes as a provirus which makes gamma-retroviral and lentiviral-based vectors attractive gene delivery vehicles for some therapeutic applications in human gene therapy ([3–6] for recent reviews). MLVs alter mouse gene expression via insertional mutagenesis (IM), which induces transcriptional or post-transcriptional deregulation of affected genes. The major determinants of retroviral IM are the LTRs, which contain viral promoter and enhancer elements in U3. Insertion of the provirus upstream from the first exon or in early introns may induce transcription of the target mouse gene driven by the retroviral promoter and this is known as promoter insertion. Transcripts may also undergo alternative splicing resulting from the use of either proviral or cryptic intronic splice sites, while integration in the 3′-end of a gene may result in truncation of the native transcript due to usage of a proviral poly (A) signal present in the LTR [3–7]. MLV enhancer activation mutagenesis is induced by enhancer elements present in U3, which augment transcription from cellular promoters over long distances by recruitment of transcription factors. IM screens based on retroviral and transposon mouse models have proven an effective approach in identifying human cancer genes [3–6, 8, 9].
IM has received increased attention due to the occurrence of adverse events following gamma-retroviral vector-based gene therapy to correct the X-linked severe combined immunodeficiency disease (X-SCID) in which a lymphoproliferative disorder was induced in patients by activation of the LMO2 oncogene [10–13]. In another clinical trial to correct the X-linked chronic granulomatous disease (X-CGD) immunodeficiency, gamma-retroviral vector integration caused activation of MDS1-EVI1, resulting in myelodysplastic syndrome (MDS) with transition to acute myeloid leukemia (AML) and development of monosomy 7 [14, 15]. Therefore, understanding fundamental principles of retrovirus-induced tumorigenesis including integration site preferences and mechanisms of host sequence deregulation is vital for the improvement of therapy safety and assessment of possible genetic disruptions in the course of treatment.
While earlier studies suggested that MLVs preferentially target DNase I hypersensitive regions, transcription start sites (TSSs) and CpG islands [16, 17], several in vitro studies and studies of cultured primary cells have shown that gamma-retroviral vectors favor integration into nucleosomal DNA [18–24]. Chromatin associated with histone methylations and acetylations including H3K4Me3 and H3K27Ac/H3K4Me1 which are common promoter and enhancer markers, respectively [25–28], are major targets of gamma-retroviral vector integration [22–24].
The bromodomain and extraterminal (BET) proteins mediate MLV integration at TSSs by tethering the integrase to acetylated H3 and H4 tails , and expression of a BET fusion protein containing the chromatin binding domain of the lentiviral integrase cofactor LEDGF/p75 results in retargeting of MLV to match the integration profile of HIV . While promoters and promoter-proximal enhancers are functionally similar [31, 32], promoter-distal enhancers may enhance transcription from thousands of bp away. One example is the distal limb bud enhancer of mouse Shh which is positioned 1 Mb from the Shh promoter [33, 34]. The mammalian genome is demarcated into coregulated enhancer and promoter units known as EPUs , and it has been shown that proviruses positioned in c-Myb upstream elements that coincide with enhancers in the EPU of this gene, establish physical contact with the c-Myb promoter through DNA looping . Taken together, this suggests that the three-dimensional structure of the genome influences promoter-distal integration-mediated long-range gene activation. It is currently not known if integration outside TSSs is mediated by BET proteins or other factors.
In the past, MLV integrations were identified using anchored-type PCR methods and Sanger sequencing, including e.g. a special 5′-end primer called a splinkerette, or inverse PCR [36–38]. Next-generation sequencing (NGS) has proven effective in the analysis of retroviral biology in diverse contexts, including HIV-1 infection and gamma-retroviral and lentiviral vector systems, and is commonly coupled with anchored-type PCR to determine positions of retroviral integration in DNA [39–46].
In this study we subjected four NMRI mouse tumors induced upon infection with the rapid T-lymphomagenic MLV SL3-3 wild-type strain to strand-specific and paired-end RNA sequencing (RNA-seq) to determine retroviral integrations genome-wide, and to identify local patterns of retroviral IM. We then used a dataset of 6,117 SL3-3 integrations derived from lymphoid tumors of more than 2,000 mice to determine the colocalization propensity of proviruses with ENCODE immunoprecipitation with sequencing (ChIP-seq) enhancer and promoter markers in vivo. Our study unravels novel mechanisms of retroviral IM involving bidirectional and tandem-type activation patterns, as well as more complex patterns including activation of major unannotated transcripts and combination-type activation where transcription is induced by promoter insertion, chimeric alternative splicing and enhancer activation by a single provirus. We also show that the majority of proviruses in tumors are located at ChIP-seq H3K4Me1/H3K27Ac-enriched positions irrespective of their distance to TSSs, suggesting that insertion into cellular regulatory regions is highly specific and subject to positive selection during tumorigenesis for sustaining long-range gene activation.
Chimeric sequencing reads expose retroviral integration sites
In this study RNA-seq was used without specific enrichment to map integrations, and simultaneously determine transcript expression levels at sites of integration. We made whole-transcriptome libraries of four thymic tumors induced by wild-type SL3-3 in mice of an inbred NMRI strain. The four tumors are referred to as 324, 327, 359 and 410. We obtained ~290 million 101-base reads exceeding an average of ~70 million single reads per tumor (online available sequencing data). Integrations were mapped by analysis of 14 paired-end sequence signatures that expose genetic structural alterations in tumors (Additional file 1: Figure S1). We manually examined the evidence for each retroviral-mouse chimeric position and assigned 92 integrations supported by chimeric fusions, and 44 regions for which fusions were not directly covered in sequencing (Figure 1A). Table 1 provides an overview of the RNA-seq and other integration datasets used in this study. The complete lists of integrations from RNA-seq and the splinkerette-based PCR screen of NMRI mice in Table 1 are provided in Additional files 2 and 3, respectively.
We observed a proportional increase in the number of integrations with increasing sequencing depth, and with the exception of tumor 359, the fraction of chimeric read pairs aligning in a direction from mouse to viral sequence, and vice versa, were nearly evenly distributed (Figure 1A). This indicates that both mouse and retroviral promoters contribute to chimeric transcription in tumors. In tumor 359, the nearly one-sided distribution of chimeric read pairs (from mouse to viral sequence) resulted from integrations into growth-factor independence 1 (Gfi1), where chimeric transcription is enhanced by high levels of mouse transcription (described below).
Comparison of RefSeq gene annotations showed significant overlap (P = 2.03E-30, hypergeometric probability) (Figure 1B) of genes assigned from integrations supported by chimeric fusions to those assigned from the conventional integration datasets in Table 1. Since these datasets comprise only 5,183 or 17% of RefSeq annotations gene assignment from the corresponding RNA-seq integrations is nonrandom. The distribution of the integrations with respect to the closest genes is shown in Figure 1C. Genes assigned from RNA-seq integrations that were not supported by chimeric fusions also showed significant overlap (P = 1.87e-08, hypergeometric probability) with genes assigned from the conventional screens. However, the potential risk of misassignment is higher since the exact locations of proviruses are not known (the annotations are available in Additional file 2).
We confirmed 59 of 92 (64%) integrations supported by chimeric fusions in DNA analyses using PCR (described in Methods). Since the majority of integrations that were not confirmed obtained a minimum coverage in sequencing (Figure 1E), it is possible that an unknown fraction did not amplify due to low copy numbers. In total, 22 out of 33 integrations that were not confirmed by DNA analyses could be assigned to a previously tagged gene (Figure 1D and E), while two out of the remaining 11 integrations were possibly misassigned in Figure 1B when considering their proximity to previous integrations (Figure 1D, purple arrows). Assuming that integrations which were neither confirmed by DNA analyses, nor assigned to a previously tagged gene are false-positives, the estimated error-rate is <12%. It should be noted that integrations supported by chimeric fusions showed a distribution at enhancer peak midpoints comparable to the BALB/c and NMRI integration datasets in Table 1 (described below).
We identified several integrations in RTCGD-assigned common integration sites (CISs) showing that multiple CISs are targeted in tumors: Arf6, Ccnd3, Chd9, Coro1a, Frat1, Gfi1/Evi5, Hivep1, Hsp90b1, Ikzf1, Kis2, Mef2c, Mir17, Myb, Mycn, Pvt1, Rasgrp1 and Thada (Additional file 2). While the four tumors harbored integrations in Gfi1 and ecotropic viral-integration Site 5 (Evi5) (the Gfi1/Evi5 locus), tumor 359 contained six integrations at this locus, four of which were positioned in sense in the 3′-UTR of Gfi1 (Figure 1E). The expression of this gene was comparably high in tumors (sequencing data available online). The differential coverage at chimeric junctions in the Gfi1 3′-UTR in tumor 359 most likely represents the expansion of distinct cellular subpopulations. Tumor 327 contained six interspersed integrations at the Gfi-1/Evi-5 locus, while tumor 410 contained two integrations at this locus. Tumor 324 contained only a single integration at the Gfi1/Evi5 locus however it contained two integrations at the Myc/Pvt1 locus. We have previously shown the occurrence of more than one integration at the same position in the same tumor, indicating that insertion at such positions trigger the onset of tumorigenesis and/or are facilitated by earlier stage or parallel mutations in individual leukemias  (discussed below).
Novel principles of retroviral insertional mutagenesis
In the following sections we describe five novel mechanisms of IM that are shown in their genomic context in Figure 2. The activation mechanisms are distinguished by diverse patterns of mouse transcription deregulation including: bidirectional activation of mouse transcription on opposite sides of a provirus (Figure 2A); tandem-type activation of distal genes that are positioned adjacently on the same DNA strand (Figure 2B); sense/antisense-type activation of genes located on opposite strands (Figure 2C); activation of genes that are not direct targets of retroviral insertion (Figure 2D); combination-type IM where enhancer activation, alternative chimeric splicing and retroviral promoter insertion are used by a single provirus to alter the expression pattern of a mouse gene (Figure 2E). To confirm common activation patterns in other tumors harboring integrations at the same positions (Table 2) we used quantitative real-time PCR (qPCR) (Figure 3) and/or rapid amplification of cDNA ends (RACE) (described below). Fragments per kb of exon per million fragments mapped (FPKM) values for transcripts at the loci in Figure 2 as well as additional loci (described below) for individual tumors are shown in Additional file 4.
Bidirectional activation of Klf7and unknown RNA transcripts
Integration at the Klf7 locus induces a complex bidirectional-type activation pattern where major deregulation of mouse sequence is induced on both sides of the provirus (Figure 2A). This results in transcriptional activation of a large 165 kb region that is not annotated (dashed line), upregulation of Klf7 as well as local activation of mouse sequence at the provirus (red arrow). The Klf7 locus contains two integration clusters one of which is positioned ~35 kb downstream of Klf7 while the other is located in the first intron of the gene. Klf7 was found consistently overexpressed in tumors harboring proviral insertions (Figure 3A) however activation of the 165 kb region was restricted to two tumors (324 and 2110) that contained in sense integrations in the downstream cluster (Figure 4A). Therefore, activation of the 165 kb region appears to be dependent on the orientation of proviruses in the downstream cluster.
To determine whether the tumors shared common transcripts in the unannotated region, they were subjected to 5′-RACE which resulted in the identification of several unknown transcripts initiating at alternative start sites (Figure 4A). The longest transcript spanned ~90 kb within the activated region, while one transcript resulted from five splicing events. We did not identify any similarity hits in common gene or RNA databases (data not shown), suggesting that the transcripts are either random products or perhaps the remains of a gene that has lost sequence similarity to any known gene.
In principle, bidirectional activation could be thought to originate from the transcriptional activities of both LTRs of an intact provirus. However, based on RNA-seq chimeric transcription was only detected at the 3′-LTR/mouse junction, proceeding in a direction opposite to that of Klf7 and the unannotated 165 kb region (Figure 4A). Therefore, transcription is induced locally by retroviral promoter insertion, while Klf7 and the 165-kb region appear to be activated by enhancer activation mutagenesis.
Concurrent transcription from both LTRs was observed at other loci, including e.g. the Tmem30b/Prkch locus where antisense transcription (previously described by us ) and retroviral promoter insertion by the 5′-LTR and 3′-LTR, respectively, was detected by RNA-seq and 5′-RACE (Additional file 1: Figure S2). Integration at this locus also results in bidirectional mouse sequence activation, including transcriptional activation of unannotated sequence. Other loci for which unannotated transcripts were detected are described below.
Integrations at Ccr9 tandem-activate Lztfl1 and Slc6a20a
We observed retrovirus-induced activation of independent transcripts at other loci including Slc6a20a/Lztfl1/Ccr9 and Syn2/AK038749 (described below). Integration at Ccr9 induces tandem-activation of the distal genes Slc6a20a and Lztfl1 which are located adjacently in the same orientation (Figure 2B), as well as increased transcription of opposite polarity (red arrow). The Slc6a20a/Lztfl1/Ccr9 locus contains five integrations clustered in a region immediately upstream of Ccr9. Therefore, based on closest distance, Ccr9 would be expected to be a target of retroviral IM. However, we did not find this gene to be systematically deregulated in tumors suggesting that Ccr9 is not activated by proviruses (Figure 3B). SLC6A20 encodes an amino acid transporter, and is duplicated in mice, i.e. paralogues Slc6a20a and Slc6a20b. Slc6a20b was not deregulated in any of the tumors harboring integrations at Ccr9 (data not shown).
Therefore integration at Ccr9 activates Slc6a20a and Lztfl1, exclusively, in three independent tumors analyzed. We speculate if the tandem-type activation pattern is induced by integration into a shared regulatory region through which proviruses can deregulate both Slc6a20a and Lztfl1. The regulation of genes through shared enhancers has been described previously and appears to be a regulatory principle that is applicable genome-wide [25, 56, 57].
Sense/antisense activation of Syn2and non-coding AK038749
The sense/antisense activation mechanism at Syn2/AK038749 induces expression of the brain genes Syn2 and non-coding AK038749 located in opposite orientation (Figures 2C and 3C). In human glioblastoma multiforme, SYN2 and TIMP4, a metalloproteinase-encoding gene located within an intron of SYN2, are subjects to reciprocal deregulation . We did not detect Timp4 expression in any tumor by qPCR analysis using three different primer pairs (data not shown). We confirmed expression of AK038749 by 5′-RACE (Figure 4B). Interestingly, the integrations at Syn2/AK038749 are located in an enhancer region marked by the histone marker H3K4Me1 that is coregulated with the Syn2 promoter, marked by H3K4Me3 and H3K4Me1, in mouse cerebellum and cortex  (Figure 4B). Since the AK038749 and Syn2 transcription initiation sites are positioned only ~200 bp from each other, it is possible that proviruses activate a bidirectional mouse promoter. In fact, promoters of many coding genes transcribe non-coding RNAs in the opposite direction .
We have previously reported on the activation of neuronal neurogranin (Nrgn) in T-cell lymphomas induced by SL3-3 integration in the Esam/Vsig2/Nrgn/Siae/Spa17 locus where the Nrgn expression level in lymphoid tumors corresponded to brain levels . Therefore, the deregulation pattern suggests that enhancer-promoter associations in mouse lymphomas, which are normally active in brain, are reinstated following retroviral integration, resulting in activation of oppositely located brain genes at this locus.
Integrations at the Col4a5/Irs4 locus activate Irs4only
The activation pattern at the Col4a5/Irs4 locus shows resemblance to deregulations detected at the Ccr9 locus. However, in this case one gene (Col4a5), functions as a hotspot for retroviral integrations that induce expression of another distal gene (Irs4) (Figure 2D), without affecting the expression of Col4a5 itself in any tumor analyzed (Figure 3D, and data not shown). The Col4a5/Irs4 locus contains in the excess of 30 integrations delimited by two major clusters that intersect at Irs4. Col4a5 contains 10 of these sites, one of which is located more than 220 kb from the Irs4 promoter. The finding that IRS4 but not COL4A5 activation in T-ALL involves translocation of both genes to the T-cell receptor beta locus  is particularly interesting considering that insertion of the retroviral enhancer mimics such oncogenic rearrangements. We detected an activation pattern at the Wwox locus fairly similar to the patterns observed at the Klf7, Col4a5/Irs4 and Ccr9 loci (Additional file 1: Figure S3). Integration into Wwox induces increased expression of Wwox as well as activation of distal Maf and unannotated mouse sequence of opposite polarity both of which are located outside the actual integration target.
Celf2is activated by a combination-type mechanism of insertional mutagenesis
The deregulation shown in Figure 2E results from a combination-type activation mechanism where several modes of IM are employed by a provirus to alter the expression of Celf2, whereby a non-chimeric transcript variant initiating from a far upstream position is activated, including also increased transcription in the provirus-containing intron (red arrow).
Remapping of chimeric sequencing reads against the Celf2 reference sequence containing the integrated provirus made it possible to link upstream and downstream Celf2 sequence to the provirus (Figure 4D). RNA splicing is mediated by the viral env splice acceptor site, and involves three splice donor (SD) sites: the canonical SD located ~39 kb upstream from the insertion site; intronic iSD1 and iSD2, located ~38 and ~2 kb from the integration site, respectively. We confirmed splicing between the canonical SD site and the env splice site by RT-PCR and Sanger sequencing (data not shown). In the 3′-end of Celf2, splicing is mediated by iSD3 and iSD4, located 15 bp and 180 bp from the integration site, respectively. The coverage by sequencing reads at the 3′-LTR/mouse junction indicated that Celf2 is also activated by retroviral promoter insertion. Transcription initiation at the canonical U3/R position was confirmed by 5′- RACE as shown in Figure 4D.
We wanted to confirm the activation pattern in other tumors containing proviruses in the Celf2 integration cluster (Figure 4C). Due to inadequate tumor tissue the analysis included one tumor (128) in which the provirus is integrated in the opposite orientation. In this tumor only the 3′-end transcript levels were increased, showing that no transcript variant is expressed (Figure 3E). The implications of alternatively spliced variants of the CELF/Bruno-like family members is not fully understood, however differential expression of Celf2 isoforms has been related to separate tissues, as well as fetal versus adult developmental stages . Since, the tumors originate from similarly-aged adult mice (described in Methods) the non-chimeric variant initiating at a far upstream position in 410 is most likely induced by retroviral enhancer-activation and not a developmentally regulated difference. The lack of 5′-end activation in tumor 128 is possibly attributable to mouse strain or tissue differences and/or to the different position and orientation of the provirus in this tumor.
Long-range IM is sustained by integration into enhancers in end-stage tumors
Gamma-retroviral vectors show a strong propensity for insertion into nucleosomal DNA, and in this study the integration-enrichment at enhancers of promoter-distal integration clusters was examined to determine epigenetic features underlying possible long-range gene activation in end-stage tumors. We exploited publically available ENCODE ChIP-seq data from Ren’s laboratory  (Additional file 1: Figure S4) and major splinkerette-based PCR integration datasets derived from tumors of more than 2,000 BALB/c and NMRI mice (Table 1). To account for possible bias introduced by genomic regions containing high numbers of integrations a supplementary reduced dataset was also analyzed which excluded highly tagged regions (described in Methods).
We found that the integrations were distributed in clusters peaking at immediate (~1,000 bp), intermediate (~8 kb) and distal (>10 kb) positions relative to UCSC TSSs (Table 3 and Additional file 1: Figure S5). Thus the majority of proviruses are not located in the immediate vicinity of TSSs. The enhancer-colocalization analyses are summarized in Table 3. The far majority of integrations colocalize with regions enriched in H3K4Me1 and/or H3K27Ac irrespective of their distance to TSSs, while the fraction of integrations colocalizing with the promoter marker H3K4Me3 decreases as expected.
We determined integration-enrichments of promoter-distal integrations at ChIP-seq peaks (Figure 5), and observed strong enrichments at H3K4Me1, H3K27Ac and H3K4Me3 regions of lymphoid tissues compared to random (P < 0.001 for all enrichments shown, described in Methods). The chromosomal distribution of integrations colocalizing directly in Figure 5 is shown in Additional file 1: Figure S6. At a distance of 5,000 bp from the juncture of ChIP-seq peaks the sequential decrease in colocalizing integrations was nearly proportional to random for markers including H3K4Me1 and the insulator CTCF (data not shown), as a result of the larger amount of the genome these markers comprise when extended to 5,000 bp (Figure 5). Considering direct overlaps, the integration enrichments reached as high as ~50 and ~30-fold (H3K27Ac) for the complete and reduced integration sets, respectively, showing that proviruses target enhancers directly. In total 2,757 of 3,329 (83%) promoter-distal integrations were positioned at enhancer markers of the lymphoid tissues, while the number for the reduced dataset was 803 of 1,175 (68%) (Table 3). The distribution of promoter-distal integrations relative to specific ChIP-seq enhancer markers is shown in Additional file 1: Figure S7.
We observed consistently lowered enrichment of integrations at enhancers of brain indicating that retroviral integration targets enhancers tissue-specifically, and therefore the enrichments in most cases reflect integration into enhancers that are common to brain and the lymphoid tissues (data not shown). The distribution of promoter-distal integrations with respect to the midpoints of H3K4Me1 and H3K27Ac ChIP-seq peaks shows that the proviruses form a dense cluster peaking at ~800 bp in case of lymphoid enhancers (i.e. spleen and thymus), while the distribution is markedly distorted in case of brain enhancers (i.e. cortex and cerebellum) (Figure 6).
Considering overlapping integrations in Table 3 at intermediate and distal positions, ~12-17% of the integrations overlap both H3K4Me1/H3K27Ac and H3K4Me3 markers. Enrichment of both enhancer and promoter markers may be caused by physical interaction between these regulatory elements as predicted by DNA looping, whereby the histone modification is assigned to the interacting region [63–65]. It could also represent regions at unannotated TSSs, where proviruses are located at promoters and promoter-proximal enhancers resulting in a colocalization pattern similar to the promoter-immediate integrations in Table 3.
In this study we used RNA-seq and ENCODE ChIP-seq data to determine positions of integration, local patterns of mouse sequence deregulation and integration-enhancer colocalization in end-stage tumors induced by the SL3-3 retrovirus. We mapped 136 integrations in only four tumors. 59 out of 92 (64%) integrations supported by chimeric fusions were confirmed by DNA analyses (Figure 1D and E), while the majority of these (~90%) showed an expected distribution close to enhancer midpoints (Figure 6). While the number of integrations in only four tumors may seem high comparably high numbers of integrations in transposon and retroviral tumor models have been reported previously [45, 66]: e.g. deep sequencing analysis of mouse mammary tumor virus (MMTV)-induced tumors, revealed an average of 27 integrations per tumor, a number that was significantly higher when including integrations supported by only a single chimeric sequencing read .
We cannot completely exclude that some insertions that had low coverage in RNA-seq may represent endogenous false-positives even-though they were unique for each tumor. We anticipate that factors such as viral mutation and/or recombination during tumorigenesis may complicate validation of genuine integration events, and would therefore require optimization of PCR conditions and cloning of bands that differ from an expected size. Mutant viruses were indeed detected in this study by Sanger sequencing including several 72-bp repeat variants of SL3-3 which are known to accumulate in end-stage tumors [67, 68] (data not shown). MLVs also frequently recombine with endogenous sequences to generate replication-competent derivative viruses which add to the complexity of MLV-induced tumors [69, 70].
The majority of integrations could not be correlated with local transcriptional activation by Cufflinks differential expression analysis (online available sequencing data) using known gene annotations though infrequent activation of mouse (unannotated) sequence was observed at several positions as described throughout this manuscript. We have previously shown that integration into major CISs, including Myc and Rasgrp1, does not consistently induce measureable changes in transcript levels, indicating that deregulations may be imposed post-transcriptionally, in some cases . This could happen following integration within a gene where premature polyadenylation induced by the retroviral poly (A) signal results in a truncated transcript and a protein with possibly altered function rather than increased gene expression, reviewed in [4, 5]. In addition, searching for disrupted target regions is also complicated by the fact that deregulations may be imposed over long distances. Alternatively, subpopulations of cells that constitute a minor fraction of clonal tumors may become out-averaged in RNA-seq, thereby masking events of retrovirus-induced transcription activation.
We found multiple integrations in the same loci, and in the same tumors, including Gfi-1/Evi-5 and Myc/Pvt1. The probability of finding multiple integrations in the same region in one tumor by chance is exceedingly small. Therefore, integrations into such positions are either pivotal for tumorigenesis or promoted in a context-dependent manner. An example is LMO2 and IL2RG cooperation in leukemia cases among X-SCID patients where insertional mutagenesis of LMO2 results in an increased growth advantage in the presence of IL2RG[71, 72]. Therefore, earlier stage or parallel mutations may predispose progenitor cells for gaining an even further growth advantage following insertion into certain loci. Integration into CISs is commonly considered important as such integrations should mark disruptions that triggered the onset of tumorigenesis in progenitor cells , and in this study numerous integrations in RTCGD-assigned CISs were detected e.g. Ccnd3, Frat1, Gfi1/Evi5, Myb and Myc/Pvt1 (Additional file 2). Cancer genomes frequently acquire mutations that alter disease progression and become subject to purifying selection in late stages of tumorigenesis [74–76] suggesting that integrations in subsidiary loci may contribute to progression rather than onset.
The human orthologs of the genes described in this study are implicated in cancer (Figure 2 and Additional file 1: Figures S2-S3 which show deregulations of the Prkch and Wwox loci, respectively). Elevated levels of the transcription factor KLF7 is associated with minimal residual disease and relapse following chemotherapeutic treatment of childhood acute lymphoblastic leukemia (ALL) and KLF7 also promote early T cell survival [77, 78]. SYN2 encodes a neuronal phosphoprotein (Synapsin 2) and is deregulated in human glioblastoma multiforme and breast cancer [58, 79]. As described already, IRS4 is activated in T-ALL following COL4A5/IRS4 translocation involving the T-cell receptor beta locus . The CUGBP-ETR3-like factors and Bruno-like (or CELF/Bruno-like) family of RNA-binding proteins regulate RNA splicing, translation and mRNA stability. The family comprises six members including CELF2, which functions as a tumor suppressor in colon cancer . TMEM30B and PRKCH which are located adjacently on human chromosome 14 have been implicated in brain and breast cancer, respectively [81, 82]. SLC6A20, LZTFL1 and CCR9 are located on 3p21.3 in a region that is often found eliminated in tumors, referred to as CER1 (commonly eliminated region 1). Chromosome 3 abnormalities have been proposed to mediate tumor formation due to loss of putative tumor suppressor genes [83–85]. WWOX also functions as a tumor suppressor in cervical cancer where it induces apoptosis and inhibits proliferation . Although such insertional activation of tumor suppressors appears at least multifaceted, tumorigenesis may be augmented by overexpression of tumor suppressors which themselves may also exhibit oncogene properties context-dependently [87, 88].
It has become widely established that retroviral vector integration favors nucleosomal over naked DNA, including the epigenetic promoter and enhancer markers H3K4Me3 and H3K4Me/H3K27Ac, respectively [18–24]. Importantly, the studies by De Rijck et al. and Sharma et al. have exposed a fundamental mechanism for integration of MLVs at TSSs mediated by BET proteins [29, 30]. Therefore, BET proteins appear to form the MLV counterpart of HIV LEDGF/p75. These studies are particularly interesting with respect to understanding fundamental integration site selection patterns, as well as the underlying chromatin structure that promotes host transcription activation by an integrated retrovirus. While the majority of MLVs do not target TSSs it is currently not known whether BET proteins perform a universal function in directing proviruses to distal regions.
Short-term models provide limited information, if any, on the selection pressure during tumorigenesis, and in our study integration-enhancer colocalizations in end-stage tumors from more than 2,000 mice was examined. We found that ~45% of the integrations were concentrated in clusters at immediate and intermediate distances from TSSs, while the remaining integrations were dispersed throughout an extended promoter-distal cluster (Table 3 and Additional file 1: Figure S5). We found that proviruses located promoter-distally showed a strong propensity for integration at enhancers (83%) comparable to integrations positioned in closer vicinity to TSSs (~90%), and this tendency was also confirmed in a reduced integration dataset where colocalization is not biased by highly tagged genomic regions (Figure 5 and Table 3). Considering that ~50% of ENCODE H3K4Me1 and H3K27Ac ChIP-seq sequence of spleen and thymus (together comprising ~181 Mb in total) are located promoter distally  (i.e. beyond 10 kb from UCSC TSSs ), the proportion of colocalizing integrations in promoter-distal regions, as well as the dense clustering of integrations at enhancer midpoints (Figure 6) accentuates that enhancers are indeed major targets of MLV integration in tumors. Our results are comparable to the studies by De Ravin et al.  and LaFave et al. , where it was found that ~87% of vector integrations in CD34+ cells overlap H3K4Me1-enriched ChIP-seq peaks throughout the genome , and that integration is driven by promoters and strong enhancers in human HepG2 and K562 cells .
We found that a considerable fraction of proviruses did not directly target enhancers rather they appeared positionally offset or scattered randomly in the genome (Figures 5 and 6). For the reduced dataset offset integrations accounted for 32% of all promoter-distal integrations, and 22% in case of intermediately positioned proviruses (Table 3). The integration pattern is likely influenced by the chromatin state of individual end-stage tumors [90, 91], and therefore it is difficult to estimate an absolute fraction of integrations that target enhancers which optimally requires the sequencing of hundreds of tumors in this particular case. Considering the looping model in which a promoter-distal enhancer is brought into proximity of a promoter, perhaps positional offsets or displacements of retroviral integrations reflect local structural features of the genome whereby the LTR becomes favorably positioned for interaction with a host promoter. Taking into account the size of a provirus or a retroviral vector, it seems plausible that some degree of local structural reorganization of the genome should accompany integration. In the study by Zhang et al.  the authors observed positional expansion, comprising several kb, of enhancers (H3K4Me1) in tumor cells containing proviruses at the c-Myb locus. It would be interesting to address at a large-scale how the local chromatin environment changes in response to retroviral integration in tumors, as well as determine BET protein binding sites genome-wide. This should show if BET tethering comprises a mechanism for integration outside TSSs, and possibly also account for integrations offset from enhancers.
In its basic form, a retroviral vector contains a therapeutic gene in place of the gag, pol, and env genes, and is delivered in the form of a replication-defective virus particle. Although vector-based treatments have proven highly effective in human clinical trials IM constitutes a major safety concern due to the development of leukemias in a minority of patients following treatment of SCID and CGD [10–15] (described above). In this study we have exposed several novel principles of gamma-retroviral-induced deregulations which altogether share a prominent complexity that reaches beyond any previously described disruption (Figure 2 and Additional file 1: Figures S2-S3).
We have shown several examples where more than a single gene is subject to activation from the same or opposite DNA strands including both proximal genes and genes whose initiation sites are positioned distantly from each other (the AK038749/Syn2 and Slc6a20a/Lztfl1/Ccr9 loci). We have shown that unannotated RNAs are transcribed at sites of integration as well as in more distal regions (the Klf7, Tmem30b/Prkch and Wwox loci), and that the provirus may activate transcription outside the actual gene in which it is integrated (the Col4a5/Irs4 and Wwox loci). We have also shown occurrences of combination-type activation patterns where retroviral promoter insertion and enhancer activation mutagenesis are employed by a provirus to alter the expression of a single gene (Celf2) as well as genomic regions outside genes (the Klf7 and Tmem30b/Prkch loci).
It may seem counterintuitive that retrovirus-induced deregulations should in general be restricted to a single mode of mutagenesis however there is a remarkable absence of studies reporting the use of combination-type activation mechanisms. One reason could be that it is difficult without prior knowledge to dissect deregulations involving chimeric and non-chimeric transcripts, respectively, to the overall expression pattern at a target locus considering also (1) transcription from the non-infected allele, and (2) unannotated transcripts. We believe that combination-type deregulation may represent a more general mechanism exposed by whole-transcriptome RNA-seq in our study.
The promoters of many coding genes transcribe non-coding RNAs in the opposite direction, and are therefore bidirectional per se, and in cancers, non-coding RNAs specific for certain malignant phenotype as well as pseudogenes are commonly expressed [93, 94]. Moreover, regulation through enhancer-promoter interaction is not strictly pairwise, rather enhancers may be shared between separate promoters to intricately coregulate the expression of more than one gene [25, 56]. Enhancers may contribute to the establishment of a intra and interchromosomal three-dimensional regulatory networks , and they may also be positioned in genes located adjacently to the genes that they regulate . Noncoding intergenic transcription may regulate nearby protein-coding genes and an L1 LINE at AZU1 provides an example where expression of this non-LTR retrotransposon correlates with the expression of multiple surrounding genes in addition to AZU1. Considering the human globin locus, an ERV-9 LTR-element modulates long-range transcription factor occupancies at multiple cis-linked genes thereby coordinating gene switching during hematopoiesis, and it also activates intergenic RNAs at low levels as a result of transient DNA looping with multiple intergenic sites . Therefore, at least for retrotransposons transcriptional activities are sustained which present a complexity comparable to the transcription patterns in MLV-induced tumors.
Taken together the deregulations described by us comply with emerging principles of complex higher-order genome regulation and show that MLVs have evolved to hijack such routes to activate multiple regions resulting in complex and long-ranging deregulations that are difficult to evaluate using conventional methods. Our analyses of integration-enhancer colocalizations in tumors also strongly support such models. Therefore, we anticipate that retrovirus-induced deregulations of equivalent or similar complexity are broadly applicable to other loci throughout the genome.
The analysis of MLV-induced mouse tumors using RNA-seq has revealed novel mechanisms of retroviral insertional mutagenesis resulting in deregulations that reach beyond a single and nearest annotated gene target. Awareness of this previously undescribed layer of complexity regarding host sequence activations may prove important for elucidating adverse effects in retroviral-based gene therapies. We have also shown that wild-type gamma-retroviruses are positioned at enhancers of lymphoid tumors irrespective of their distance to TSSs, showing that insertion into regulatory regions is highly specific and also subject to positive selection during tumorigenesis. This suggests a mechanism whereby the provirus exploits the higher-order genome regulatome for sustaining long-range deregulations in tumors. This study should prove useful for extrapolating adverse outcomes of retroviral vector therapies, and for understanding fundamental cellular regulatory principles.
Availability of supporting data
Sequences are available from the NCBI short read archive (accession no. SRP041565). Primer sequences not provided in Additional file 5 are available on request.
Mouse infection and splinkerette-based PCR mapping of integrations
Inbred BALB/c and NMRI mice were infected with the rapid lymphomagenic MLV SL3-3 strain as described previously [48–50]. Upon disease or appearance of tumors (in 60 to 70 days) the mice were sacrificed and spleen and thymus organs eviscerated and kept frozen at −80°C. The approx. size of thymic tumors was ~1.5 cm in the longest dimension while that of splenic tumors was ~3 cm. DNA was extracted using the DNeasy Tissue kit (Qiagen) and integrations determined using an automated high-throughput splinkerette-PCR method [36, 38]. We have previously published results from the screen of ~2,000 BALB/c mice (Table 1). The NMRI dataset in Additional file 3 contains integrations from several cohorts of mice infected with SL3-3. The total number of mice was 120, including a cohort of 22 mice from which the four tumors subjected to RNA-seq originate.
Total RNA was extracted from thymus and spleen by the guanidinium thiocyanate-phenol-chloroform method (TRIzol, Invitrogen) and purified on spin columns using the RNeasy Mini Kit (Qiagen), including on-column DNase I treatment (Qiagen), and stored at −80°C. The concentration and purity was routinely determined on a NanoDrop (Thermo Scientific), and RNA integrity was assessed by microfluidics technology on the Agilent BioAnalyzer 2100 (Agilent).
cDNA library preparation and sequencing
Total RNA from four SL3-3-induced NMRI mouse tumors was depleted of rRNA by hybridization using Ribo-Zero rRNA Removal Kit for Human/Mouse/Rat (Epicentre). The samples were then concentrated with RNA Clean & Concentrator (Zymo Research), and cDNA libraries prepared with the ScriptSeq RNA-Seq Library Preparation Kit (Epicentre) following the manufacturer’s guidelines. The libraries incorporated barcodes for multiplexed sequencing using the RNA-Seq Barcode Primers for Illumina (Epicentre). The resulting cDNA libraries were amplified by limited 12-cycle PCR and size-fractionated using high-percentage agarose gel electrophoresis. The size of the sequenced libraries was 400–450 bp including adaptors. The libraries were selectively quantified on a LightCycler 480 II Real-Time PCR System (Roche), using the KAPA Library Quantification Assay (Kapa Biosystems). The libraries were pooled in equal amounts on one paired-end flow cell lane using the cBot cluster generation process (Illumina), and sequenced on an Illumina HiSeq 2000 producing 2x101-base sequencing reads.
RNA-seq-based integration mapping
To filter chimeric read pairs a custom reference was built manually using BLAT  based on the content of SL3-3 LTR-like sequence in the NCBIM37/mm9 mouse reference assembly (Additional file 1: Figure S8). In brief, SL3-3 LTR [GenBank:AF169256] sequence windows were aligned to the mouse genome in 5 bp increments to identify mouse sequence windows showing 90% or higher identity by sequence. The custom reference included also SL3-3 sequences [GenBank:X00862 and GenBank:AF169256]. Subsequently, 25 bp sequencing reads were mapped against the custom reference with the Burrows-Wheeler Aligner (BWA)  followed by sequestration of read pairs in which only one of the mates mapped. To identify integrations full length reads were quality trimmed and aligned separately with BLAT against the NCBIM37/mm9 reference assembly modified to contain SL3-3 on a separate chromosome. The chimeric pairs produce either one of 14 paired-end sequence signatures, based on U3, U5, SD or no chimeric fusion coverage (Additional file 1: Figure S1). We manually examined the evidence for each of the integrations and integrations common to the samples were considered endogenous false-positives based on the improbability of insertion at the exact same position in several samples. Integrations supported by a chimeric fusion contained intact sequence, at one or both LTRs. In any case the minimum requirement was 40 bp of unambiguously mapped sequence. The integrations are listed in Additional file 2.
Analysis of ENCODE data
ChIP-seq data  were obtained from ENCODE (Additional file 1: Figure S4), while SL3-3 integrations were obtained from the screens of NMRI and BALB/c mice shown in Table 1. We determined the colocalization of integrations in end-stage tumors with enhancers using complete (n = 6,117) and reduced (n = 2,127) integration datasets. The analyses were performed using BEDTools . As control size-matched random datasets were used in iterative simulations (n = 1,000) for each analysis shown in Figure 5. In brief, the ChIP-seq datasets were annotated according to UCSC TSS annotations, excluding unmappable sequence. The median peak length and sequence content (the size in bp) of each dataset within the searched genome was then used for picking random intervals without replacement from NCBIM37/mm9. We defined direct overlaps based on extension of ChIP-seq peaks by 1,250 bp. In each case intersection was performed using 1 bp integration coordinates. The determination of empirical p-values based on random sampling was performed as previously described . The reduced SL3-3 integration dataset was assembled by clustering integrations in bins of 2,500 bp across the genome. Clusters larger than this (containing more than two integration positions) were excluded. The position of integrations for clusters that contained more than one integration site was defined as the cluster midpoint. Out of the 2,127 coordinates in the reduced integration dataset 1,981 (93%) represented individual integration sites.
Sequence alignment, visualization and expression analysis
FASTQ files were processed including adapter clipping and quality trimming, and aligned to the NCBI37/mm9 mouse assembly with Bowtie/TopHat . The visualizations in Figures 2 and 4, as well as Additional file 1: Figures S2-S3 are based on forward read alignments for which the sequence coverage of each sample was computed strand-specifically with BEDTools  and scaled using Cufflinks library metrics . The images were generated using the UCSC Genome Browser and post-processed using standard vector imaging tools. The remapping of Celf2 in Figure 4E was performed using Celf2 (chr2:6,453,742-7,029,527) and the SL3-3 proviral sequence (built using [GenBank:X00862, GenBank:AF169256]) inserted at the predetermined integration site (Additional file 2) including also the TSD. The Cufflinks transcript expression values for the loci shown in Figure 2, as well as (Additional file 1: Figures S2-S3) are available in Additional file 4.
DNA analysis, RT-PCR, qPCR, RACE and cloning
DNA and RNA were copurified with TRIzol (Invitrogen) followed by DNA back-extraction and clean-up using QIAamp spin columns (Qiagen). To confirm integrations (Figure 1E) 40 ng of genomic DNA was subjected to 40 cycles of PCR (DreamTaq, Thermo Scientific) using three or more primer pairs designed to yield products of size differences discernible by gel electrophoresis (Additional file 1: Figure S9). In case a distinct band pattern did not appear select products were Sanger sequenced (described below). For RT-PCR and the qPCR-results shown in Figure 3 primer sequences from OriGene were used with exceptions. 500 ng of RNA was used for cDNA synthesis with the qScript cDNA SuperMix (Quanta Biosciences). PCR-reactions (DreamTaq, Thermo Scientific) were performed in a 2720 Thermo Cycler (Applied Biosystems) using standard cycling conditions, i.e. annealing at Tm -5°C. For qPCR measurements a Stratagene Mx3000P cycler was used (Applied Biosciences). Samples were measured in duplicates using PerfeCta SYBR Green FastMix (Quanta Biosciences) in 10 μL reactions, and beta-actin as reference. As controls, tumors were used without known integrations at the loci shown in Figure 2 and included for tumors not subjected to RNA-seq a minimum of three samples from the same mouse strain (BALB/c or NMRI) and tissue (thymus or spleen). Above-threshold values were defined as Ct values greater than 30 cycles of amplification. RACE was carried out using the SMARTer RACE cDNA amplification Kit (Clontech) in nested reactions. cDNA synthesis for 5′-RACE included random priming, and was performed according to manufacturer’s guidelines. For Sanger sequencing PCR products were cloned in the pCR4-TOPO vector using the TOPO TA Cloning Kit for Sequencing (Invitrogen). The primers used in this study are shown in Additional file 5.
Murine leukemia virus
Long terminal repeat
Common integration site
Chromatin Immunoprecipitation with sequencing
Enhancer and promoter unit
Transcriptional start site
Retrovirus and Transposon tagged Cancer Gene Database
Quantitative real-time PCR
Rapid amplification of cDNA ends
Fragments per kilobase of transcript per million mapped reads
Target site duplication.
Hays EF, Levy JA: Differences in lymphomagenic properties of AKR mouse retroviruses. Virology. 1984, 138 (1): 49-57. 10.1016/0042-6822(84)90146-6.
Pedersen FS, Crowther RL, Tenney DY, Reimold AM, Haseltine WA: Novel leukaemogenic retroviruses isolated from cell line derived from spontaneous AKR tumour. Nature. 1981, 292 (5819): 167-170. 10.1038/292167a0.
Knight S, Collins M, Takeuchi Y: Insertional mutagenesis by retroviral vectors: current concepts and methods of analysis. Curr Gene Ther. 2013, 13 (3): 211-227. 10.2174/1566523211313030006.
Nowrouzi A, Glimm H, Von Kalle C, Schmidt M: Retroviral vectors: post entry events and genomic alterations. Viruses. 2011, 3 (5): 429-455.
Cavazza A, Moiani A, Mavilio F: Mechanisms of retroviral integration and mutagenesis. Hum Gene Ther. 2013, 24 (2): 119-131. 10.1089/hum.2012.203.
Gabriel R, Schmidt M, Von Kalle C: Integration of retroviral vectors. Curr Opin Immunol. 2012, 24 (5): 592-597. 10.1016/j.coi.2012.08.006.
Dabrowska MJ, Dybkaer K, Johnsen HE, Wang B, Wabl M, Pedersen FS: Loss of MicroRNA targets in the 3′ untranslated region as a mechanism of retroviral insertional activation of growth factor independence 1. J Virol. 2009, 83 (16): 8051-8061. 10.1128/JVI.00427-09.
Uren AG, Kool J, Berns A, Van Lohuizen M: Retroviral insertional mutagenesis: past, present and future. Oncogene. 2005, 24 (52): 7656-7672. 10.1038/sj.onc.1209043.
Mikkers H, Berns A: Retroviral insertional mutagenesis: tagging cancer pathways. Adv Cancer Res. 2003, 88: 53-99.
Gaspar HB, Parsley KL, Howe S, King D, Gilmour KC, Sinclair J, Brouns G, Schmidt M, Von Kalle C, Barington T, Jakobsen MA, Christensen HO, Al Ghonaium A, White HN, Smith JL, Levinsky RJ, Ali RR, Kinnon C, Thrasher AJ: Gene therapy of X-linked severe combined immunodeficiency by use of a pseudotyped gammaretroviral vector. Lancet. 2004, 364 (9452): 2181-2187. 10.1016/S0140-6736(04)17590-9.
Hacein-Bey-Abina S, Le Deist F, Carlier F, Bouneaud C, Hue C, De Villartay JP, Thrasher AJ, Wulffraat N, Sorensen R, Dupuis-Girod S, Fischer A, Davies EG, Kuis W, Leiva L, Cavazzana-Calvo M: Sustained correction of X-linked severe combined immunodeficiency by ex vivo gene therapy. N Engl J Med. 2002, 346 (16): 1185-1193. 10.1056/NEJMoa012616.
Cavazzana-Calvo M, Hacein-Bey S, de Saint BG, Gross F, Yvon E, Nusbaum P, Selz F, Hue C, Certain S, Casanova JL, Bousso P, Deist FL, Fischer A: Gene therapy of human severe combined immunodeficiency (SCID)-X1 disease. Science. 2000, 288 (5466): 669-672. 10.1126/science.288.5466.669.
Howe SJ, Mansour MR, Schwarzwaelder K, Bartholomae C, Hubank M, Kempski H, Brugman MH, Pike-Overzet K, Chatters SJ, de Ridder D, Gilmour KC, Adams S, Thornhill SI, Parsley KL, Staal FJ, Gale RE, Linch DC, Bayford J, Brown L, Quaye M, Kinnon C, Ancliff P, Webb DK, Schmidt M, von Kalle C, Gaspar HB, Thrasher AJ: Insertional mutagenesis combined with acquired somatic mutations causes leukemogenesis following gene therapy of SCID-X1 patients. J Clin Invest. 2008, 118 (9): 3143-3150. 10.1172/JCI35798.
Stein S, Ott MG, Schultze-Strasser S, Jauch A, Burwinkel B, Kinner A, Schmidt M, Kramer A, Schwable J, Glimm H, Koehl U, Preiss C, Ball C, Martin H, Gohring G, Schwarzwaelder K, Hofmann WK, Karakaya K, Tchatchou S, Yang R, Reinecke P, Kuhlcke K, Schlegelberger B, Thrasher AJ, Hoelzer D, Seger R, von Kalle C, Grez M: Genomic instability and myelodysplasia with monosomy 7 consequent to EVI1 activation after gene therapy for chronic granulomatous disease. Nature medicine. 2010, 16 (2): 198-204. 10.1038/nm.2088.
Ott MG, Schmidt M, Schwarzwaelder K, Stein S, Siler U, Koehl U, Glimm H, Kuhlcke K, Schilz A, Kunkel H, Naundorf S, Brinkmann A, Deichmann A, Fischer M, Ball C, Pilz I, Dunbar C, Du Y, Jenkins NA, Copeland NG, Luthi U, Hassan M, Thrasher AJ, Hoelzer D, von Kalle C, Seger R, Grez M: Correction of X-linked chronic granulomatous disease by gene therapy, augmented by insertional activation of MDS1-EVI1, PRDM16 or SETBP1. Nature medicine. 2006, 12 (4): 401-409. 10.1038/nm1393.
Vijaya S, Steffen DL, Robinson HL: Acceptor sites for retroviral integrations map near DNase I-hypersensitive sites in chromatin. J Virol. 1986, 60 (2): 683-692.
Wu X, Li Y, Crise B, Burgess SM: Transcription start regions in the human genome are favored targets for MLV integration. Science. 2003, 300 (5626): 1749-1751. 10.1126/science.1083413.
Cattoglio C, Pellin D, Rizzi E, Maruggi G, Corti G, Miselli F, Sartori D, Guffanti A, Di Serio C, Ambrosi A, De Bellis G, Mavilio F: High-definition mapping of retroviral integration sites identifies active regulatory elements in human multipotent hematopoietic progenitors. Blood. 2010, 116 (25): 5507-5517. 10.1182/blood-2010-05-283523.
Roth SL, Malani N, Bushman FD: Gammaretroviral Integration into Nucleosomal Target DNA In Vivo. J Virol. 2011, 85 (14): 7393-7401. 10.1128/JVI.00635-11.
Biasco L, Ambrosi A, Pellin D, Bartholomae C, Brigida I, Roncarolo MG, Di Serio C, Von Kalle C, Schmidt M, Aiuti A: Integration profile of retroviral vector in gene therapy treated patients is cell-specific according to gene expression and chromatin conformation of target cell. EMBO Mol Med. 2011, 3 (2): 89-101. 10.1002/emmm.201000108.
Moiani A, Miccio A, Rizzi E, Severgnini M, Pellin D, Suerth JD, Baum C, De Bellis G, Mavilio F: Deletion of the LTR enhancer/promoter has no impact on the integration profile of MLV vectors in human hematopoietic progenitors. PLoS One. 2013, 8 (1): e55721-10.1371/journal.pone.0055721.
Lafave MC, Varshney GK, Gildea DE, Wolfsberg TG, Baxevanis AD, Burgess SM: MLV integration site selection is driven by strong enhancers and active promoters. Nucleic Acids Res. 2014, 42 (7): 4257-4269. 10.1093/nar/gkt1399.
Santoni FA, Hartley O, Luban J: Deciphering the code for retroviral integration target site selection. PLoS Comput Biol. 2010, 6 (11): e1001008-10.1371/journal.pcbi.1001008.
De Ravin SS, Su L, Theobald N, Choi U, Macpherson JL, Poidinger M, Symonds G, Pond SM, Ferris AL, Hughes SH, Malech HL, Wu X: Enhancers are major targets for murine leukemia virus vector integration. J Virol. 2014, 88 (8): 4504-4513. 10.1128/JVI.00011-14.
Shen Y, Yue F, McCleary DF, Ye Z, Edsall L, Kuan S, Wagner U, Dixon J, Lee L, Lobanenkov VV, Ren B: A map of the cis-regulatory sequences in the mouse genome. Nature. 2012, 488 (7409): 116-120. 10.1038/nature11243.
Creyghton MP, Cheng AW, Welstead GG, Kooistra T, Carey BW, Steine EJ, Hanna J, Lodato MA, Frampton GM, Sharp PA, Boyer LA, Young RA, Jaenisch R: Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc Natl Acad Sci U S A. 2010, 107 (50): 21931-21936. 10.1073/pnas.1016071107.
Zhang B, Day DS, Ho JW, Song L, Cao J, Christodoulou D, Seidman JG, Crawford GE, Park PJ, Pu WT: A dynamic H3K27ac signature identifies VEGFA-stimulated endothelial enhancers and requires EP300 activity. Genome Res. 2013, 23 (6): 917-927. 10.1101/gr.149674.112.
Tian Y, Jia Z, Wang J, Huang Z, Tang J, Zheng Y, Tang Y, Wang Q, Tian Z, Yang D, Zhang Y, Fu X, Song J, Liu S, van Velkinburgh JC, Wu Y, Ni B: Global mapping of H3K4me1 and H3K4me3 reveals the chromatin state-based cell type-specific gene regulation in human Treg cells. PLoS One. 2011, 6 (11): e27770-10.1371/journal.pone.0027770.
Sharma A, Larue RC, Plumb MR, Malani N, Male F, Slaughter A, Kessl JJ, Shkriabai N, Coward E, Aiyer SS, Green PL, Wu L, Roth MJ, Bushman FD, Kvaratskhelia M: BET proteins promote efficient murine leukemia virus integration at transcription start sites. Proc Natl Acad Sci U S A. 2013, 110 (29): 12036-12041. 10.1073/pnas.1307157110.
De Rijck J, de Kogel C, Demeulemeester J, Vets S, El Ashkar S, Malani N, Bushman FD, Landuyt B, Husson SJ, Busschots K, Gijsbers R, Debyser Z: The BET family of proteins targets moloney murine leukemia virus integration near transcription start sites. Cell reports. 2013, 5 (4): 886-894. 10.1016/j.celrep.2013.09.040.
Taher L, Smith R, Kim M, Ahituv N, Ovcharenko I: Sequence signatures extracted from proximal promoters can be used to predict distal enhancers. Genome Biol. 2013, 14 (10): R117-10.1186/gb-2013-14-10-r117.
Maston GA, Evans SK, Green MR: Transcriptional regulatory elements in the human genome. Annu Rev Genomics Hum Genet. 2006, 7: 29-59. 10.1146/annurev.genom.7.080505.115623.
Lettice LA, Heaney SJ, Purdie LA, Li L, De Beer P, Oostra BA, Goode D, Elgar G, Hill RE, De Graaff E: A long-range Shh enhancer regulates expression in the developing limb and fin and is associated with preaxial polydactyly. Hum Mol Genet. 2003, 12 (14): 1725-1735. 10.1093/hmg/ddg180.
Sagai T, Hosoya M, Mizushina Y, Tamura M, Shiroishi T: Elimination of a long-range cis-regulatory module causes complete loss of limb-specific Shh expression and truncation of the mouse limb. Development. 2005, 132 (4): 797-803. 10.1242/dev.01613.
Zhang J, Markus J, Bies J, Paul T, Wolff L: Three Murine Leukemia Virus Integration Regions within 100 Kilobases Upstream of c-myb Are Proximal to the 5′ Regulatory Region of the Gene through DNA Looping. J Virol. 2012, 86 (19): 10524-10532. 10.1128/JVI.01077-12.
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 (1): 153-159. 10.1038/ng950.
Uren AG, Mikkers H, Kool J, van der Weyden L, Lund AH, Wilson CH, Rance R, Jonkers J, van Lohuizen M, Berns A, Adams DJ: A high-throughput splinkerette-PCR method for the isolation and sequencing of retroviral insertion sites. Nat Protoc. 2009, 4 (5): 789-798. 10.1038/nprot.2009.64.
Devon RS, Porteous DJ, Brookes AJ: Splinkerettes–improved vectorettes for greater efficiency in PCR walking. Nucleic Acids Res. 1995, 23 (9): 1644-1645. 10.1093/nar/23.9.1644.
Ustek D, Sirma S, Gumus E, Arikan M, Cakiris A, Abaci N, Mathew J, Emrence Z, Azakli H, Cosan F, Cakar A, Parlak M, Kursun O: A genome-wide analysis of lentivector integration sites using targeted sequence capture and next generation sequencing technology. Infect Genet Evol. 2012, 12 (7): 1349-1354. 10.1016/j.meegid.2012.05.001.
Ciuffi A, Barr SD: Identification of HIV integration sites in infected host genomic DNA. Methods. 2011, 53 (1): 39-46. 10.1016/j.ymeth.2010.04.004.
Arens A, Appelt J-U, Bartholomae CC, Gabriel R, Paruzynski A, Gustafson D, Cartier N, Aubourg P, Deichmann A, Glimm H, von Kalle C, Schmidt M: Bioinformatic clonality analysis of next-generation sequencing-derived viral vector integration sites. Human Gene Therapy Methods. April 2012, 23 (2): 111-118.
Schopman NC, Willemsen M, Liu YP, Bradley T, Van Kampen A, Baas F, Berkhout B, Haasnoot J: Deep sequencing of virus-infected cells reveals HIV-encoded small RNAs. Nucleic Acids Res. 2012, 40 (1): 414-427. 10.1093/nar/gkr719.
Lefebvre G, Desfarges S, Uyttebroeck F, Munoz M, Beerenwinkel N, Rougemont J, Telenti A, Ciuffi A: Analysis of HIV-1 expression level and sense of transcription by high-throughput sequencing of the infected cell. J Virol. 2011, 85 (13): 6205-6211. 10.1128/JVI.00252-11.
Cesana D, Sgualdino J, Rudilosso L, Merella S, Naldini L, Montini E: Whole transcriptome characterization of aberrant splicing events induced by lentiviral vector integrations. J Clin Invest. 2012, 122 (5): 1667-1676. 10.1172/JCI62189.
Koudijs MJ, Klijn C, van der Weyden L, Kool J, ten Hoeve J, Sie D, Prasetyanti PR, Schut E, Kas S, Whipp T, Cuppen E, Wessels L, Adams DJ, Jonkers J: High-throughput semiquantitative analysis of insertional mutations in heterogeneous tumors. Genome Res. 2011, 21 (12): 2181-2189. 10.1101/gr.112763.110.
Chang ST, Sova P, Peng X, Weiss J, Law GL, Palermo RE, Katze MG: Next-generation sequencing reveals HIV-1-mediated suppression of T cell activation and RNA processing and regulation of noncoding RNA expression in a CD4+ T cell line. mBio. 2011, 2: e00134-11.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L: Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012, 7 (3): 562-578. 10.1038/nprot.2012.016.
Beck-Engeser GB, Lum AM, Huppi K, Caplen NJ, Wang BB, Wabl M: Pvt1-encoded microRNAs in oncogenesis. Retrovirology. 2008, 5: 4-10.1186/1742-4690-5-4.
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 (49): 18680-18684. 10.1073/pnas.0609030103.
Lum AM, Wang BB, Li L, Channa N, Bartha G, Wabl M: Retroviral activation of the mir-106a microRNA cistron in T lymphoma. Retrovirology. 2007, 4: 5-10.1186/1742-4690-4-5.
Akagi K, Suzuki T, Stephens RM, Jenkins NA, Copeland NG: RTCGD: retroviral tagged cancer gene database. Nucleic Acids Res. 2004, 32 (Database issue): D523-D527.
Dabrowska MJ, Ejegod D, Lassen LB, Johnsen HE, Wabl M, Pedersen FS, Dybkaer K: Gene expression profiling of murine T-cell lymphoblastic lymphoma identifies deregulation of S-phase initiating genes. Leuk Res. 2013, 37 (10): 1383-1390. 10.1016/j.leukres.2013.04.012.
Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26 (6): 841-842. 10.1093/bioinformatics/btq033.
Rasmussen MH, Ballarin-Gonzalez B, Liu J, Lassen LB, Fuchtbauer A, Fuchtbauer EM, Nielsen AL, Pedersen FS: Antisense transcription in gammaretroviruses as a mechanism of insertional activation of host genes. J Virol. 2010, 84 (8): 3780-3788. 10.1128/JVI.02088-09.
Kikuchi R, Yagi S, Kusuhara H, Imai S, Sugiyama Y, Shiota K: Genome-wide analysis of epigenetic signatures for kidney-specific transporters. Kidney Int. 2010, 78 (6): 569-577. 10.1038/ki.2010.176.
Eun B, Sampley ML, Good AL, Gebert CM, Pfeifer K: Promoter cross-talk via a shared enhancer explains paternally biased expression of Nctc1 at the Igf2/H19/Nctc1 imprinted locus. Nucleic Acids Res. 2013, 41 (2): 817-826. 10.1093/nar/gks1182.
Sanyal A, Lajoie BR, Jain G, Dekker J: The long-range interaction landscape of gene promoters. Nature. 2012, 489 (7414): 109-113. 10.1038/nature11279.
Dong H, Luo L, Hong S, Siu H, Xiao Y, Jin L, Chen R, Xiong M: Integrated analysis of mutations, miRNA and mRNA expression in glioblastoma. BMC Syst Biol. 2010, 4: 163-10.1186/1752-0509-4-163.
Wei W, Pelechano V, Jarvelin AI, Steinmetz LM: Functional consequences of bidirectional promoters. Trends Genet. 2011, 27 (7): 267-276. 10.1016/j.tig.2011.04.002.
Nielsen AA, Kjartansdottir KR, Rasmussen MH, Sorensen AB, Wang B, Wabl M, Pedersen FS: Activation of the brain-specific neurogranin gene in murine T-cell lymphomas by proviral insertional mutagenesis. Gene. 2009, 442 (1–2): 55-62.
Karrman K, Kjeldsen E, Lassen C, Isaksson M, Davidsson J, Andersson A, Hasle H, Fioretos T, Johansson B: The t(X;7)(q22;q34) in paediatric T-cell acute lymphoblastic leukaemia results in overexpression of the insulin receptor substrate 4 gene through illegitimate recombination with the T-cell receptor beta locus. Br J Haematol. 2009, 144 (4): 546-551. 10.1111/j.1365-2141.2008.07453.x.
Li D, Bachinski LL, Roberts R: Genomic organization and isoform-specific tissue expression of human NAPOR (CUGBP2) as a candidate gene for familial arrhythmogenic right ventricular dysplasia. Genomics. 2001, 74 (3): 396-401. 10.1006/geno.2001.6558.
Hon GC, Hawkins RD, Ren B: Predictive chromatin signatures in the mammalian genome. Hum Mol Genet. 2009, 18 (R2): R195-R201. 10.1093/hmg/ddp409.
Bulger M, Groudine M: Looping versus linking: toward a model for long-distance gene activation. Genes Dev. 1999, 13 (19): 2465-2477. 10.1101/gad.13.19.2465.
Robertson AG, Bilenky M, Tam A, Zhao Y, Zeng T, Thiessen N, Cezard T, Fejes AP, Wederell ED, Cullum R, Euskirchen G, Krzywinski M, Birol I, Snyder M, Hoodless PA, Hirst M, Marra MA, Jones SJ: Genome-wide relationship between histone H3 lysine 4 mono- and tri-methylation and transcription factor binding. Genome Res. 2008, 18 (12): 1906-1917. 10.1101/gr.078519.108.
Keng VW, Villanueva A, Chiang DY, Dupuy AJ, Ryan BJ, Matise I, Silverstein KA, Sarver A, Starr TK, Akagi K, Tessarollo L, Collier LS, Powers S, Lowe SW, Jenkins NA, Copeland NG, Llovet JM, Largaespada DA: A conditional transposon-based insertional mutagenesis screen for genes associated with mouse hepatocellular carcinoma. Nat Biotechnol. 2009, 27 (3): 264-274. 10.1038/nbt.1526.
Martiney MJ, Rulli K, Beaty R, Levy LS, Lenz J: Selection of reversions and suppressors of a mutation in the CBF binding site of a lymphomagenic retrovirus. J Virol. 1999, 73 (9): 7599-7606.
Morrison HL, Soni B, Lenz J: Long terminal repeat enhancer core sequences in proviruses adjacent to c-myc in T-cell lymphomas induced by a murine retrovirus. J Virol. 1995, 69 (1): 446-455.
Evans LH, Cloyd MW: Friend and Moloney murine leukemia viruses specifically recombine with different endogenous retroviral sequences to generate mink cell focus-forming viruses. Proc Natl Acad Sci U S A. 1985, 82 (2): 459-463. 10.1073/pnas.82.2.459.
Cuypers HT, Selten G, Quint W, Zijlstra M, Maandag ER, Boelens W, Van Wezenbeek P, Melief C, Berns A: Murine leukemia virus-induced T-cell lymphomagenesis: integration of proviruses in a distinct chromosomal region. Cell. 1984, 37 (1): 141-150. 10.1016/0092-8674(84)90309-X.
Dave UP, Akagi K, Tripathi R, Cleveland SM, Thompson MA, Yi M, Stephens R, Downing JR, Jenkins NA, Copeland NG: Murine leukemias with retroviral insertions at Lmo2 are predictive of the leukemias induced in SCID-X1 patients following retroviral gene therapy. PLoS Genet. 2009, 5 (5): e1000491-10.1371/journal.pgen.1000491.
Hacein-Bey-Abina S, Von Kalle C, Schmidt M, McCormack MP, Wulffraat N, Leboulch P, Lim A, Osborne CS, Pawliuk R, Morillon E, Sorensen R, Forster A, Fraser P, Cohen JI, de Saint BG, Alexander I, Wintergerst U, Frebourg T, Aurias A, Stoppa-Lyonnet D, Romana S, Radford-Weiss I, Gross F, Valensi F, Delabesse E, Macintyre E, Sigaux F, Soulier J, Leiva LE, Wissler M, et al: LMO2-associated clonal T cell proliferation in two patients after gene therapy for SCID-X1. Science. 2003, 302 (5644): 415-419. 10.1126/science.1088547.
Abel U, Deichmann A, Nowrouzi A, Gabriel R, Bartholomae CC, Glimm H, Von Kalle C, Schmidt M: Analyzing the number of common integration sites of viral vectors–new methods and computer programs. PLoS One. 2011, 6 (10): e24247-10.1371/journal.pone.0024247.
Mullighan CG, Phillips LA, Su X, Ma J, Miller CB, Shurtleff SA, Downing JR: Genomic analysis of the clonal origins of relapsed acute lymphoblastic leukemia. Science. 2008, 322 (5906): 1377-1380. 10.1126/science.1164266.
Stratton MR, Campbell PJ, Futreal PA: The cancer genome. Nature. 2009, 458 (7239): 719-724. 10.1038/nature07943.
Torkamani A, Schork NJ: Identification of rare cancer driver mutations by network reconstruction. Genome Res. 2009, 19 (9): 1570-1578. 10.1101/gr.092833.109.
Flotho C, Coustan-Smith E, Pei D, Cheng C, Song G, Pui C-H, Downing JR, Campana D: A set of genes that regulate cell proliferation predicts treatment outcome in childhood acute lymphoblastic leukemia. Blood. 2007, 110 (4): 1271-1277. 10.1182/blood-2007-01-068478.
Schuettpelz LG, Gopalan PK, Giuste FO, Romine MP, van Os R, Link DC: Kruppel-like factor 7 overexpression suppresses hematopoietic stem and progenitor cell function. Blood. 2012, 120 (15): 2981-2989. 10.1182/blood-2012-02-409839.
Park KS, Youn HJ, Jung SH: A study for expression and biological function of N-myc downstream regulated gene 2 in breast cancer. J Breast Cancer. 2007, 10 (3): 180-192. 10.4048/jbc.2007.10.3.180.
Ramalingam S, Ramamoorthy P, Subramaniam D, Anant S: Reduced expression of RNA binding protein CELF2, a putative tumor suppressor gene in colon cancer. Immunogastroenterology. 2012, 1 (1): 27-33. 10.7178/ig.1.1.7.
Masso-Welch PA, Winston JS, Edge S, Darcy KM, Asch H, Vaughan MM, Ip MM: Altered expression and localization of PKC eta in human breast tumors. Breast Cancer Res Treat. 2001, 68 (3): 211-223. 10.1023/A:1012265703669.
Perez-Magan E, Campos-Martin Y, Mur P, Fiano C, Ribalta T, Garcia JF, Rey JA, Rodriguez de Lope A, Mollejo M, Melendez B: Genetic alterations associated with progression and recurrence in meningiomas. J Neuropathol Exp Neurol. 2012, 71 (10): 882-893. 10.1097/NEN.0b013e31826bf704.
Yang Y, Kiss H, Kost-Alimova M, Kedra D, Fransson I, Seroussi E, Li J, Szeles A, Kholodnyuk I, Imreh MP, Fodor K, Hadlaczky G, Klein G, Dumanski JP, Imreh S: A 1-Mb PAC contig spanning the common eliminated region 1 (CER1) in microcell hybrid-derived SCID tumors. Genomics. 1999, 62 (2): 147-155. 10.1006/geno.1999.5952.
Kholodnyuk I, Kost-Alimova M, Kashuba V, Gizatulin R, Szeles A, Stanbridge EJ, Zabarovsky ER, Klein G, Imreh S: A 3p21.3 region is preferentially eliminated from human chromosome 3/mouse microcell hybrids during tumor growth in SCID mice. Genes Chromosomes Cancer. 1997, 18 (3): 200-211. 10.1002/(SICI)1098-2264(199703)18:3<200::AID-GCC6>3.0.CO;2-5.
Ingvarsson S: Tumor suppressor genes on human chromosome 3 and cancer pathogenesis. Cancer Genomics Proteomi. 2005, 2 (4): 247-253.
Qu J, Lu W, Li B, Lu C, Wan X: WWOX induces apoptosis and inhibits proliferation in cervical cancer and cell lines. Int J Mol Cell Med. 2013, 31 (5): 1139-1147.
Messina S, Frati L, Leonetti C, Zuchegna C, Di Zazzo E, Calogero A, Porcellini A: Dual-specificity phosphatase DUSP6 has tumor-promoting properties in human glioblastomas. Oncogene. 2011, 30 (35): 3813-3820. 10.1038/onc.2011.99.
Hong L, Li X, Jin H, Yan L, Wu K, Ding J, Zhao Y, Guo W, Fan D: Up-regulation of tumor suppressor genes might promote the malignant phenotype of cancer cells. Med Hypotheses. 2007, 69 (6): 1379-10.1016/j.mehy.2007.03.029.
Karolchik D, Hinrichs AS, Kent WJ: The UCSC genome browser. Curr Protoc Hum Genet. 2011, Chapter 18: Unit18 16-editorial board, Jonathan L Haines [et al.]
Gonzalez-Perez A, Jene-Sanz A, Lopez-Bigas N: The mutational landscape of chromatin regulatory factors across 4,623 tumor samples. Genome Biol. 2013, 14 (9): r106-10.1186/gb-2013-14-9-r106.
Zeitz MJ, Ay F, Heidmann JD, Lerner PL, Noble WS, Steelman BN, Hoffman AR: Genomic interaction profiles in breast cancer reveal altered chromatin architecture. PLoS One. 2013, 8 (9): e73974-10.1371/journal.pone.0073974.
Orekhova AS, Rubtsov PM: Bidirectional promoters in the transcription of mammalian genomes. Biochemistry. 2013, 78 (4): 335-341.
Kalyana-Sundaram S, Kumar-Sinha C, Shankar S, Robinson DR, Wu YM, Cao X, Asangani IA, Kothari V, Prensner JR, Lonigro RJ, Iyer MK, Barrette T, Shanmugam A, Dhanasekaran SM, Palanisamy N, Chinnaiyan AM: Expressed pseudogenes in the transcriptional landscape of human cancers. Cell. 2012, 149 (7): 1622-1634. 10.1016/j.cell.2012.04.041.
Cheetham SW, Gruhl F, Mattick JS, Dinger ME: Long noncoding RNAs and the genetics of cancer. Br J Cancer. 2013, 108 (12): 2419-2425. 10.1038/bjc.2013.233.
Ong C-T, Corces VG: Enhancer function: new insights into the regulation of tissue-specific gene expression. Nat Rev Genet. 2011, 12 (4): 283-293.
Bulger M, Groudine M: Functional and mechanistic diversity of distal transcription enhancers. Cell. 2011, 144 (3): 327-339. 10.1016/j.cell.2011.01.024.
Faulkner GJ, Kimura Y, Daub CO, Wani S, Plessy C, Irvine KM, Schroder K, Cloonan N, Steptoe AL, Lassmann T, Waki K, Hornig N, Arakawa T, Takahashi H, Kawai J, Forrest AR, Suzuki H, Hayashizaki Y, Hume DA, Orlando V, Grimmond SM, Carninci P: The regulated retrotransposon transcriptome of mammalian cells. Nature genetics. 2009, 41 (5): 563-571. 10.1038/ng.368.
Pi W, Zhu X, Wu M, Wang Y, Fulzele S, Eroglu A, Ling J, Tuan D: Long-range function of an intergenic retrotransposon. Proc Natl Acad Sci U S A. 2010, 107 (29): 12992-12997. 10.1073/pnas.1004139107.
Kent WJ: BLAT–the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664. 10.1101/gr.229202. Article published online before March 2002.
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.
North BV, Curtis D, Sham PC: A note on the calculation of empirical P values from Monte Carlo procedures. Am J Hum Genet. 2002, 71 (2): 439-441. 10.1086/341527.
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-1111. 10.1093/bioinformatics/btp120.
This work was supported by grants from The Danish Cancer Society and The Danish Research Council for Health and Disease to FSP, and by the National Institutes of Health Grant R01 AI041570 to MW. We are grateful for the enthusiastic discussions with Mads Heilskov Rasmussen during the project, and for the technical assistance provided by Lone Højgaard Nielsen.
The authors declare that they have no competing interest.
MS designed the deep-sequencing experiments, analyzed data and drafted the manuscript. MW carried out the splinkerette-based-PCR determination of MLV integrations. IRR and MS performed qPCR and RACE experiments. FSP conceived the study and participated in its design and direction. All authors read and consented to the contents of the manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Paired-end RNA-seq signatures expose retroviral integration sites. Figure S2. The Tmem30b/Prkch locus is deregulated by a bidirectional-type activation mechanism. Figure S3. Integration in Wwox induces overexpression of distal Maf and activation of unannotated transcription outside Wwox. Figure S4. ChIP-seq datasets from ENCODE. Figure S5. Integrations in end-stage tumors form clusters at immediate, intermediate and distal positions from TSSs. Figure S6. Chromosomal distributions of promoter-distal integrations in the complete and reduced integration datasets. Figure S7. Distribution of colocalizing integrations with respect to H3K4Me1 and H3K27Ac ChIP-seq peaks from spleen and thymus. Figure S8. Content of SL3-3 LTR-like sequence in the mouse genome assembly (NCBIM37/mm9). Figure S9. PCR confirmation of integrations identified in RNA sequencing. (PPTX 4 MB)
Additional file 2: Integrations identified using RNA-seq. In columns 2–3 the integrations for each tumor are numbered according to Figures 1C to E. The positions and orientations (sense, S and antisense, AS) of the integrations with respect to the reference genome, and the presence (+) or absence (−) of chimeric fusions are shown in columns 4–8. In cases the chimeric fusion point was not covered in sequencing, the integration positions correspond to the 5′-end of the murine read mates. In columns 9–12 the distance of integrations to the nearest RefSeq annotation are shown as well as the positions of proviruses (exon, intron, or outside). Columns 13–15 show if the genes have been tagged (+) or not (−) in the RTCGD or in screens of BALB/c and NMRI mice (described in the main text). The distance of the RNA-seq integrations to the nearest integrations from these screens is indicated in column 16. The four tumors (324 through 410) had also been analyzed in the NMRI mouse screen. Seven out of 13 integrations identified in the latter screen were detected with deep sequencing (Additional files 2 and 3, compared). The analyses were performed on different tumor sections. (XLSX 36 KB)
Additional file 3: Integrations identified in NMRI mice using splinkerette-based PCR. The dataset contains integrations from several cohorts of NMRI mice infected with SL3-3 including a cohort of 20 mice from which the four mouse tumors (324 through 410) subjected to RNA-seq originate (described in Methods). The integrations are sorted in order of decreasing number of tags at a locus. The first column shows tumors used in the present study for RNA-seq and/or qPCR and RACE. The layout column shows the orientations of integrations with respect to the genes: e.g. at chr5:108,167,153 the provirus is located 6,660 bp downstream from Evi5 in the opposite (antisense) orientation relative to this gene. (XLSX 52 KB)
Additional file 4: RNA-seq transcript expression values. This table lists Cufflinks fragments per kilobase of transcript per million mapped reads (FPKM)  values using the UCSC mouse transcriptome prediction track for the deregulated loci described in the manuscript. The lower (conf_lo) and upper (conf_hi) bounds of the 95% confidence interval of transcript abundances is also indicated. The complete RNA-seq data is available online. (XLSX 19 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Sokol, M., Wabl, M., Ruiz, I.R. et al. Novel principles of gamma-retroviral insertional transcription activation in murine leukemia virus-induced end-stage tumors. Retrovirology 11, 36 (2014). https://doi.org/10.1186/1742-4690-11-36
- Insertional mutagenesis
- Deep sequencing
- Chromatin immunoprecipication with sequencing (ChIP-seq)
- Retroviral integration sites
- RNA sequencing (RNA-seq)