Retroviral insertions in the VISION database identify molecular pathways in mouse lymphoid leukemia and lymphoma

AKXD recombinant inbred (RI) strains develop a variety of leukemias and lymphomas due to somatically acquired insertions of retroviral DNA into the genome of hematopoetic cells that can mutate cellular proto-oncogenes and tumor suppressor genes. We generated a new set of tumors from nine AKXD RI strains selected for their propensity to develop B-cell tumors, the most common type of human hematopoietic cancers. We employed a PCR technique called viral insertion site amplification (VISA) to rapidly isolate genomic sequence at the site of provirus insertion. Here we describe 550 VISA sequence tags (VSTs) that identify 74 common insertion sites (CISs), of which 21 have not been identified previously. Several suspected proto-oncogenes and tumor suppressor genes lie near CISs, providing supportive evidence for their roles in cancer. Furthermore, numerous previously uncharacterized genes lie near CISs, providing a pool of candidate disease genes for future research. Pathway analysis of candidate genes identified several signaling pathways as common and powerful routes to blood cancer, including Notch, E-protein, NFκB, and Ras signaling. Misregulation of several Notch signaling genes was confirmed by quantitative RT-PCR. Our data suggest that analyses of insertional mutagenesis on a single genetic background are biased toward the identification of cooperating mutations. This tumor collection represents the most comprehensive study of the genetics of B-cell leukemia and lymphoma development in mice. We have deposited the VST sequences, CISs in a genome viewer, histopathology, and molecular tumor typing data in a public web database called VISION (Viral Insertion Sites Identifying Oncogenes), which is located at . Keith C. Weiser and Bin Liu are authors that contributed equally to this work.     

Sequence and spacing requirements of a retrovirus integration site

Following infection, retroviruses insert a DNA copy of their RNA genome into the host cell genome. This integrative recombination reaction occurs at specific sites on the viral DNA: inverted repeat sequences near the termini of the linear DNA form of the viral genome. We have described elsewhere the generation and analysis of deletion mutations at one of the inverted repeat sequences in Moloney murine leukemia virus. We describe here the effects of insertion mutations made at this locus. Our results show that substantial sequence changes at the site of recombination can be tolerated, and that the spacing between the cleavage sites on the viral DNA can be expanded as well as contracted while still allowing efficient viral integration. After several rounds of virus replication, each of the insertion mutants gave rise to pseudorevertants with new alterations at the integration site.     

Molecular mechanisms of retroviral integration site selection

Retroviral replication proceeds through an obligate integrated DNA provirus, making retroviral vectors attractive vehicles for human gene-therapy. Though most of the host cell genome is available for integration, the process of integration site selection is not random. Retroviruses differ in their choice of chromatin-associated features and also prefer particular nucleotide sequences at the point of insertion. Lentiviruses including HIV-1 preferentially integrate within the bodies of active genes, whereas the prototypical gammaretrovirus Moloney murine leukemia virus (MoMLV) favors strong enhancers and active gene promoter regions. Integration is catalyzed by the viral integrase protein, and recent research has demonstrated that HIV-1 and MoMLV targeting preferences are in large part guided by integrase-interacting host factors (LEDGF/p75 for HIV-1 and BET proteins for MoMLV) that tether viral intasomes to chromatin. In each case, the selectivity of epigenetic marks on histones recognized by the protein tether helps to determine the integration distribution. In contrast, nucleotide preferences at integration sites seem to be governed by the ability for the integrase protein to locally bend the DNA duplex for pairwise insertion of the viral DNA ends. We discuss approaches to alter integration site selection that could potentially improve the safety of retroviral vectors in the clinic.     

The distribution of transgene insertion sites in barley determined by physical and genetic mapping

The exact site of transgene insertion into a plant host genome is one feature of the genetic transformation process that cannot, at present, be controlled and is often poorly understood. The site of transgene insertion may have implications for transgene stability and for potential unintended effects of the transgene on plant metabolism. To increase our understanding of transgene insertion sites in barley, a detailed analysis of transgene integration in independently derived transgenic barley lines was carried out. Fluorescence in situ hybridization (FISH) was used to physically map 23 transgene integration sites from 19 independent barley lines. Genetic mapping further confirmed the location of the transgenes in 11 of these lines. Transgene integration sites were present only on five of the seven barley chromosomes. The pattern of transgene integration appeared to be nonrandom and there was evidence of clustering of independent transgene insertion events within the barley genome. In addition, barley genomic regions flanking the transgene insertion site were isolated for seven independent lines. The data from the transgene flanking regions indicated that transgene insertions were preferentially located in gene-rich areas of the genome. These results are discussed in relation to the structure of the barley genome.     

Sint1, a common integration site in SL3-3-induced T-cell lymphomas, harbors a putative proto-oncogene with homology to the septin gene family

The murine retrovirus SL3-3 is a potent inducer of T-cell lymphomas when inoculated into susceptible newborn mice. Previously, DNAs from twenty SL3-3-induced tumors were screened by PCR for provirus integration sites. Two out of 20 tumors demonstrated clonal provirus insertion into a common region. This region has now been isolated and characterized. The region, named SL3-3 integration site 1 (Sint1), maps to the distal end of mouse chromosome 11, corresponding to human chromosome 17q25, and may be identical to a mouse mammary tumor virus integration site in a T-cell lymphoma, Pad3. Two overlapping genomic lambda clones spanning about 35 kb were isolated and used as a starting point for a search for genes in the neighborhood of the virus integration sites. A genomic fragment was used as a hybridization probe to isolate a 3-kb cDNA clone, the expression of which was upregulated in one of two tumors harboring a provirus in Sint1. The cDNA clone is predicted to encode a protein which shows 97.0% identity to a human septin-like protein encoded by a gene which has been found as a fusion partner gene of MLL in an acute myeloid leukemia with a t(11;17)(q23;q25). Together these findings raise the possibility that a proto-oncogene belonging to the septin family, and located about 15 kb upstream of the provirus integration sites, is involved in murine leukemia virus-induced T-cell lymphomagenesis.     

Long-range mapping of Mis-2, a common provirus integration site identified in murine leukemia virus-induced thymomas and located 160 kilobase pairs downstream of Myb.

The nondefective Moloney murine leukemia virus (MuLV) induces clonal or oligoclonal T-cell tumors in mice or rats. The proviruses of these nondefective MuLVs have been shown to act as insertion mutagens most frequently activating an adjacent cellular gene involved in cell growth control. Mutations by provirus insertions, recognized as common provirus integration sites, have been instrumental in identifying novel cellular genes involved in tumor formation. We have searched for new common provirus integration sites in Moloney MuLV-induced thymomas. Using cellular sequences flanking a provirus cloned from one of these tumors, we found one region, designated Mis-2, which was the target of provirus integration in a low (3%) percentage of these tumors. Mis-2 was mapped on mouse chromosome 10, approximately 160 kbp downstream of myb. The Mis-2 region may contain a novel gene involved in tumor development.     

Evidence for the involvement of pim-2, a new common proviral insertion site, in progression of lymphomas.

We have compared proviral integrations near (putative) proto-oncogenes in Moloney murine leukemia virus-induced primary and transplanted T cell lymphomas. We previously found proviruses integrated near c-myc, pim-1, and N-myc in primary tumors (Selten et al., 1984; Van Lohuizen et al., 1989a; Van Lohuizen et al., 1989b). We have now identified an additional common proviral integration site, called pim-2, that carries somatically acquired proviruses in the majority of transplanted tumors. In primary tumors integration near pim-2 is usually undetectable or present in only a minor fraction of the tumor cells. This subpopulation selectively grows out upon transplantation. Insertion near pim-2 is a relatively late event in tumorigenesis and is often preceded by proviral insertions in other common insertion sites, yielding tumor clones which carry proviruses in up to three different common insertion sites within the same cell (c-myc, pim-1 and pim-2). The data suggest that pim-2 plays an important role in tumor progression.     

Integration Profiling of Gene Function With Dense Maps of Transposon Integration

Understanding how complex networks of genes integrate to produce dividing cells is an important goal that is limited by the difficulty in defining the function of individual genes. Current resources for the systematic identification of gene function such as siRNA libraries and collections of deletion strains are costly and organism specific. We describe here integration profiling, a novel approach to identify the function of eukaryotic genes based upon dense maps of transposon integration. As a proof of concept, we used the transposon Hermes to generate a library of 360,513 insertions in the genome of Schizosaccharomyces pombe. On average, we obtained one insertion for every 29 bp of the genome. Hermes integrated more often into nucleosome free sites and 33% of the insertions occurred in ORFs. We found that ORFs with low integration densities successfully identified the genes that are essential for cell division. Importantly, the nonessential ORFs with intermediate levels of insertion correlated with the nonessential genes that have functions required for colonies to reach full size. This finding indicates that integration profiles can measure the contribution of nonessential genes to cell division. While integration profiling succeeded in identifying genes necessary for propagation, it also has the potential to identify genes important for many other functions such as DNA repair, stress response, and meiosis.     

Discovery of transgene insertion sites by high throughput sequencing of mate pair libraries

Background

Transgenesis by random integration of a transgene into the genome of a zygote has become a reliable and powerful method for the creation of new mouse strains that express exogenous genes, including human disease genes, tissue specific reporter genes or genes that allow for tissue specific recombination. Nearly 6,500 transgenic alleles have been created by random integration in embryos over the last 30 years, but for the vast majority of these strains, the transgene insertion sites remain uncharacterized.

Results

To obtain a complete understanding of how insertion sites might contribute to phenotypic outcomes, to more cost effectively manage transgenic strains, and to fully understand mechanisms of instability in transgene expression, we’ve developed methodology and a scoring scheme for transgene insertion site discovery using high throughput sequencing data.

Conclusions

Similar to other molecular approaches to transgene insertion site discovery, high-throughput sequencing of standard paired-end libraries is hindered by low signal to noise ratios. This problem is exacerbated when the transgene consists of sequences that are also present in the host genome. We’ve found that high throughput sequencing data from mate-pair libraries are more informative when compared to data from standard paired end libraries. We also show examples of the genomic regions that harbor transgenes, which have in common a preponderance of repetitive sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-367) contains supplementary material, which is available to authorized users.     

Sequence analysis of HIV-1 insertion sites in peripheral blood lymphocytes.

An essential component of the HIV-1 life cycle involves insertion in the genome of an infected cell. The site of HIV-1 integration has the potential to disrupt a gene and perturb a normal cellular function. To begin to address whether disease pathogenesis may correlate with the site of insertion, flanking cellular sequences at these HIV integrated regions were directly amplified from peripheral blood mononuclear cells DNA from a broad range of infected individuals using an inverse polymerase chain reaction strategy. Amplified flanking regions were sequenced and examined for similarity to the nucleic acid database. In this group of analyzed samples, the HIV-1 provirus was inserted within non-coding regions throughout the genome of the infected host, in which 7/14 sites were positioned in close proximity to different Alu repetitive elements while 2/14 sites were located within intron sequences. Insertions were also detected at sites without a specific gene designation but not within short tandem repetitive sequences, telomeres or centromeric repeat regions. Altogether, it is expected that this approach will yield new information on sites of integration by HIV-1 that may be associated with the pathogenic manifestations of disease progression.     

Retargeting transposon insertions by the adeno-associated virus Rep protein

The Sleeping Beauty (SB), piggyBac (PB) and Tol2 transposons are promising instruments for genome engineering. Integration site profiling of SB, PB and Tol2 in human cells showed that PB and Tol2 insertions were enriched in genes, whereas SB insertions were randomly distributed. We aimed to introduce a bias into the target site selection properties of the transposon systems by taking advantage of the locus-specific integration system of adeno-associated virus (AAV). The AAV Rep protein binds to Rep recognition sequences (RRSs) in the human genome, and mediates viral integration into nearby sites. A series of fusion constructs consisting of the N-terminal DNA-binding domain of Rep and the transposases or the N57 domain of SB were generated. A plasmid-based transposition assay showed that Rep/SB yielded a 15-fold enrichment of transposition at a particular site near a targeted RRS. Genome-wide insertion site analysis indicated that an approach based on interactions between the SB transposase and Rep/N57 enriched transgene insertions at RRSs. We also provide evidence of biased insertion of the PB and Tol2 transposons. This study provides a comparative insight into target site selection properties of transposons, as well as proof-of-principle for targeted chromosomal transposition by composite protein-protein and protein-DNA interactions.     

A high throughput method for genome-wide analysis of retroviral integration

Retroviral and lentiviral vectors integrate their DNA into the host cell genome leading to stable transgene expression. Integration preferentially occurs in the proximity of active genes, and may in some case disturb their activity, with adverse toxic consequences. To efficiently analyze high numbers of lentiviral insertion sites in the DNA of transduced cells, we developed an improved high-throughput method called vector integration tag analysis (VITA). VITA is based on the identification of Genomic Tags associated to the insertion sites, which are used as signatures of the integration events. We use the capacity of MmeI to cleave DNA at a defined distance of its recognition site, in order to generate 21 bp long tags from libraries of junction fragments between vector and cellular DNA. The length of the tags is sufficient in most cases, to identify without ambiguity an unique position in the human genome. Concatenation, cloning and sequencing of the tags allow to obtain information about 20–25 insertion sites in a single sequencing reaction. As a validation of this method, we have characterized 1349 different lentiviral vector insertion sites in transduced HeLa cells, from only 487 sequencing reactions, with a background of <2% false positive tags.     

Weak palindromic consensus sequences are a common feature found at the integration target sites of many retroviruses

Integration into the host genome is one of the hallmarks of the retroviral life cycle and is catalyzed by virus-encoded integrases. While integrase has strict sequence requirements for the viral DNA ends, target site sequences have been shown to be very diverse. We carefully examined a large number of integration target site sequences from several retroviruses, including human immunodeficiency virus type 1, simian immunodeficiency virus, murine leukemia virus, and avian sarcoma-leukosis virus, and found that a statistical palindromic consensus, centered on the virus-specific duplicated target site sequence, was a common feature at integration target sites for these retroviruses.     

Genetic alterations by human papillomaviruses in oncogenesis.

The integration sites in the cellular genome of human papillomavirus are located in chromosomal regions always associated with oncogenes or other known tumor phenotypes. Two regions, 8q24 and 12q13, are common to several cases of cervical carcinoma and can have integrated more than one type of papillomavirus DNA. These two chromosomal regions contain several genes implicated in oncogenesis. These observations strongly imply that viral integration sites of DNA tumor viruses can be used as the access point to chromosomal regions where genes implicated in the tumor phenotype are located, a situation similar to that of non-transforming retroviruses.