Identifying viral integration sites using SeqMap 2.0

Retroviral integration has been implicated in several biomedical applications, including identification of cancer-associated genes and malignant transformation in gene therapy clinical trials. We introduce an efficient and scalable method for fast identification of viral vector integration sites from long read high-throughput sequencing. Individual sequence reads are masked to remove non-genomic sequence, aligned to the host genome and assembled into contiguous fragments used to pinpoint the position of integration. AVAILABILITY AND IMPLEMENTATION: The method is implemented in a publicly accessible web server platform, SeqMap 2.0, containing analysis tools and both private and shared lab workspaces that facilitate collaboration among researchers. Available at http://seqmap.compbio.iupui.edu/.     

The avian retroviral integration protein cleaves the terminal sequences of linear viral DNA at the in vivo sites of integration.

The purified integration protein (IN) of avian myeloblastosis virus is shown to nick double-stranded oligodeoxynucleotide substrates that mimic the ends of the linear form of viral DNA. In the presence of Mg2+, nicks are created 2 nucleotides from the 3' OH ends of both the U5 plus strand and the U3 minus strand. Similar cleavage is observed in the presence of Mn2+ but only when the extent of the reaction is limited. Neither the complementary strands nor sequences representing the termini of human immunodeficiency virus type 1 DNA were cleaved at analogous positions. Analysis of a series of substrates containing U5 base substitutions has defined the sequence requirements for site-selective nicking; nucleotides near the cleavage site are most critical for activity. The minimum substrate size required to demonstrate significant activity corresponds to the nearly perfect 15-base terminal inverted repeat. This in vitro activity of IN thus produces viral DNA ends that are joined to host DNA in vivo and corresponds to an expected early step in the integrative recombination reaction. These results provide the first enzymatic support using purified retroviral proteins for a linear DNA precursor to the integrated provirus.     

Analyzing the number of common integration sites of viral vectors--new methods and computer programs

Vectors based on γ-retroviruses or lentiviruses have been shown to stably express therapeutical transgenes and effectively cure different hematological diseases. Molecular follow up of the insertional repertoire of gene corrected cells in patients and preclinical animal models revealed different integration preferences in the host genome including clusters of integrations in small genomic areas (CIS; common integrations sites). In the majority, these CIS were found in or near genes, with the potential to influence the clonal fate of the affected cell. To determine whether the observed degree of clustering is statistically compatible with an assumed standard model of spatial distribution of integrants, we have developed various methods and computer programs for γ-retroviral and lentiviral integration site distribution. In particular, we have devised and implemented mathematical and statistical approaches for comparing two experimental samples with different numbers of integration sites with respect to the propensity to form CIS as well as for the analysis of coincidences of integration sites obtained from different blood compartments. The programs and statistical tools described here are available as workspaces in R code and allow the fast detection of excessive clustering of integration sites from any retrovirally transduced sample and thus contribute to the assessment of potential treatment-related risks in preclinical and clinical retroviral gene therapy studies.     

Moloney murine leukemia virus integration protein produced in yeast binds specifically to viral att sites.

The integration protein (IN) of Moloney murine leukemia virus (MuLV), purified after being produced in yeast cells, has been analyzed for its ability to bind its putative viral substrates, the att sites. An electrophoretic mobility shift assay revealed that the Moloney MuLV IN protein binds synthetic oligonucleotides containing att sequences, with specificity towards its cognate (MuLV) sequences. The terminal 13 base pairs, which are identical at both ends of viral DNA, are sufficient for binding if present at the ends of oligonucleotide duplexes in the same orientation as in linear viral DNA. However, only weak binding was observed when the same sequences were positioned within a substrate in a manner simulating att junctions in circular viral DNA with two long terminal repeats. Binding to att sites in oligonucleotides simulating linear viral DNA was dependent on the presence of the highly conserved CA residues preceding the site for 3' processing (an IN-dependent reaction that removes two nucleotides from the 3' ends of linear viral DNA); mutation of CA to TG abolished binding, and a CA to TA change reduced affinity by at least 20-fold. Removal of either the terminal two base pairs from both ends of the oligonucleotide duplex or the terminal two nucleotides from the 3' ends of each strand did not affect binding. The removal of three 3' terminal nucleotides, however, abolished binding, suggesting an essential role for the A residue immediately upstream of the 3' processing site in the binding reaction. These results help define the sequence requirements for att site recognition by IN, explain the conservation of the subterminal CA dinucleotide, and provide a simple assay for sequence-specific IN activity.     

Oncogenes

Genetic elements responsible for key steps in the conversion of normal cells to malignancy, initially identified in oncogenic viruses, have normal cellular origins. These normal cell genes, protooncogenes, are highly conserved in an evolutionary context, and this fact, along with other data on expression and product localization, are hypothesized to play an important role in growth and differentiation. The number of such genes is limited to perhaps two dozen and they can be subdivided into families, thus implying fine-tuning functions. The conversion of protooncogenes to oncogenes can be based on mutation or selective stimulation based on chromosomal alterations, e.g., translocations, loss of control elements, or insertion of active promoters. A class of control elements, designated antioncogenes, has been tentatively identified in certain familial malignancies. Data are accumulating on the chromosomal localization of both onco- and antioncogenes, justifying prospects of defining cancer at the molecular level.     

Physical map of infectious baboon type C viral DNA and sites of integration in infected cells.

Three species of unintegrated viral DNAs were found in permissive cells infected with baboon type C virus. The major species was a 9.0-kilobase (kb) linear DNA that was infectious. A restriction endonuclease map of this DNA was constructed and oriented with respect to the viral RNA. The linear DNA had a 0.6-kb sequence repeated at each terminus. These terminal repeat sequences were required for infectivity of the viral DNA. The minor species of the unintegrated viral DNAs were covalently closed circles of 9.0 and 8.4 kb. The smaller circle was in two- to threefold excess over the larger circle. The difference appeared to be that the smaller circle lacked one of the two 0.6-kb repeat sequences found in the larger circle. Restriction endonuclease maps of the integrated viral DNAs were constructed, and the sequences on both viral DNA and cellular DNA that are involved in integration were determined. The integrated viral DNA map was identical to that of the unintegrated infectious 9.0-kb linear DNA. Therefore, a specific site in the terminal repeat sequence of the viral DNA was used to integrate with the host cell DNA. The sizes of the cellular DNA fragments were different from clone to clone but stable with cell passage. Therefore, many sites in the cell DNA can recombine with the viral DNA.     

Non-homologous integration of transforming vectors in the fungus Podospora anserina: sequences of junctions at the integration sites

Transformation of the ura5-6 mutant strain of Podospora anserina with a recombinant vector carrying the ura5+ gene often results in the integration of the transforming plasmid by non-homologous recombination outside of the genomic ura5 locus. To investigate the mechanism of such integration, we rescued the integrated plasmid from three transformants. In two cases, the rescued plasmid was highly altered compared with the original transforming vector. We cloned the junctions between plasmidic DNA and genomic DNA of the transformants and determined their nucleotide sequences. It was found that there was little homology between plasmidic and genomic DNA sequences. Moreover, in all cases deletions of plasmid sequences at the integration site had occurred. These rearrangements can be explained by the formation of multimeric plasmids prior to integration.     

DNA sequence at the integration sites of the insertion element IS1.

We have detected two independent occurrences of insertion mutations in the lacl gene of E. Coli, and have used small plasmids carrying the l gene to purify large amounts of DNA containing these insertions. Analyses with restriction endonucleases and DNA sequencing techniques establish that both insertions involve the previously characterized element IS1. In each case, the integration of IS1 into the l gene DNA is associated with a directly repeated sequence of 9 nucleotides appearing at each end of the insertion element. Since one of these sequences was present in the wild-type gene, the second sequence either preexisted in the IS1 before integration, or else was generated by the process of insertion itself. The 9 base repeat is different in both cases. We discuss the relevance of these findings to the mechanism of integration of transposable elements.     

Mapping of recombinant retrovirus integration sites that cause expression of the viral genome in murine embryonal carcinoma cells

Murine embryonal carcinoma (EC) cells do not normally express Moloney murine leukemia virus genes. Earlier, rare EC cell lines were isolated that expressed proviral neomycin resistance (neo) gene. This expression was dependent on cellular enhancer or promoter sequences that flank the proviral integration site. Four such integration sites, designated as Mint (for Moloney murine leukemia virus integration and expression sites in EC cells), have been mapped on mouse chromosomes. Minta, Mintb, Mintc and Mintd are unlinked and mapped on different chromosomes (Chr), Chr 10, Chr 1, Chr 5 and the X Chr, respectively. None of these loci appear to be linded to any known Mo-MuLV proviral integration sites previously mapped. These enhancer and promoter loci may represent a new set of genes active in undifferentiated embryonic cells.The terminology for Mint genes has been discussed with and approved by the International Committee on Standard Genetic Nomenclature for Mice (M.T. Davisson, The Jackson Laboratory, Bar Harbor, Me.).     

HIV整合位点的生物信息学分析

已有研究证明,HIV整合位点的选择,与宿主基因组功能和染色体结构特性之间存在着紧密的联系,但是具体的选择机制还不明确。使用聚类分析、特征提取、分类分析等生物信息学方法,对HIV整合位点序列进行分析和研究,挖掘HIV整合位点序列之间的关系,探索HIV整合位点的选择规律。通过实验和计算,从HIV整合位点集合中提取出了含有6个特征向量的向量集,该向量集与大部分整合位点的特征向量具有较高的相关性,从而提示了HIV整合位点选择中的规律性,即符合向量集的宿主DNA序列可为HIV的整合位点。研究结果为进一步揭示HIV整合位点的选择机制提供了可供参考的依据。     

Transgene integration in aspen: structures of integration sites and mechanism of T-DNA integration

To obtain insight into the mechanism of transferred DNA (T-DNA) integration in a long-lived tree system, we analysed 30 transgenic aspen lines. In total, 27 right T-DNA/plant junctions, 20 left T-DNA/plant junctions, and 10 target insertions from control plants were obtained. At the right end, the T-DNA was conserved up to the cleavage site in 18 transgenic lines (67%), and the right border repeat was deleted in nine junctions. Nucleotides from the left border repeat were present in 19 transgenic lines out of 20 cases analysed. However, only four (20%) of the left border ends were conserved to the processing end, indicating that the T-DNA left and right ends are treated mechanistically differently during the T-DNA integration process. Comparison of the genomic target sites prior to integration to the T-DNA revealed that the T-DNA inserted into the plant genome without any notable deletion of genomic sequence in three out of 10 transgenic lines analysed. However, deletions of DNA ranging in length from a few nucleotides to more than 500 bp were observed in other transgenic lines. Filler DNAs of up to 235 bp were observed on left and/or right junctions of six transgenic lines, which in most cases originated from the nearby host genomic sequence or from the T-DNA. Short sequence similarities between recombining strands near break points, in particular for the left T-DNA end, were observed in most of the lines analysed. These results confirm the well-accepted T-DNA integration model based on single-stranded annealing followed by ligation of the right border which is preserved by the VirD2 protein. However, a second category of T-DNA integration was also identified in nine transgenic lines, in which the right border of the T-DNA was partly truncated. Such integration events are described via a model for the repair of genomic double-strand breaks in somatic plant cells based on synthesis-dependent strand-annealing. This report in a long-lived tree system provides major insight into the mechanism of transgene integration.     

Retrovirus integration and chromatin structure: Moloney murine leukemia proviral integration sites map near DNase I-hypersensitive sites.

The chromatin conformation of mouse genome regions containing Moloney murine leukemia proviral intergration sites in two Mov mouse strains and randomly selected integration sites in virus-infected mouse 3T3 fibroblasts was analyzed. All integrations have occurred into chromosomal regions containing several DNase-hypersensitive sites, and invariably the proviral integration sites map within a few hundred base pairs of a DNase-hypersensitive site. The probability that this close association between proviral integration sites and DNase-hypersensitive sites was due to chance was calculated to be extremely low (2 X 10(-4]. Because the proviral integrations analyzed were not selected for an altered phenotype, our results suggest that DNase-hypersensitive regions are preferred targets for retrovirus integration.