Methods for integration site distribution analyses in animal cell genomes

Non-Random Pattern of Integration for Epstein-Barr Virus with Preference for Gene-Poor Genomic Chromosomal Regions into the Genome of Burkitt Lymphoma Cell Lines

Background: Epstein-Barr virus (EBV) is an oncogenic virus found in about 95% of endemic Burkitt lymphoma (BL) cases. In latently infected cells, EBV DNA is mostly maintained in episomal form, but it can also be integrated into the host genome, or both forms can coexist in the infected cells. Methods: In this study, we mapped the chromosomal integration sites of EBV (EBV-IS) into the genome of 21 EBV+ BL cell lines (BL-CL) using metaphase fluorescence in situ hybridization (FISH). The data were used to investigate the EBV-IS distribution pattern in BL-CL, its relation to the genome instability, and to assess its association to common fragile sites and episomes. Results: We detected a total of 459 EBV-IS integrated into multiple genome localizations with a preference for gene-poor chromosomes. We did not observe any preferential affinity of EBV to integrate into common and rare fragile sites or enrichment of EBV-IS at the chromosomal breakpoints of the BL-CL analyzed here, as other DNA viruses do. Conclusions: We identified a non-random integration pattern into 13 cytobands, of which eight overlap with the EBV-IS in EBV-transformed lymphoblastoid cell lines and with a preference for gene- and CpGs-poor G-positive cytobands. Moreover, it has been demonstrated that the episomal form of EBV interacts in a non-random manner with gene-poor and AT-rich regions in EBV+ cell lines, which may explain the observed affinity for G-positive cytobands in the EBV integration process. Our results provide new insights into the patterns of EBV integration in BL-CL at the chromosomal level, revealing an unexpected connection between the episomal and integrated forms of EBV.

A computational framework to assess genome-wide distribution of polymorphic human endogenous retrovirus-K In human populations

Human Endogenous Retrovirus type K (HERV-K) is the only HERV known to be insertionally polymorphic; not all individuals have a retrovirus at a specific genomic location. It is possible that HERV-Ks contribute to human disease because people differ in both number and genomic location of these retroviruses. Indeed viral transcripts, proteins, and antibody against HERV-K are detected in cancers, auto-immune, and neurodegenerative diseases. However, attempts to link a polymorphic HERV-K with any disease have been frustrated in part because population prevalence of HERV-K provirus at each polymorphic site is lacking and it is challenging to identify closely related elements such as HERV-K from short read sequence data. We present an integrated and computationally robust approach that uses whole genome short read data to determine the occupation status at all sites reported to contain a HERV-K provirus. Our method estimates the proportion of fixed length genomic sequence (k-mers) from whole genome sequence data matching a reference set of k-mers unique to each HERV-K locus and applies mixture model-based clustering of these values to account for low depth sequence data. Our analysis of 1000 Genomes Project Data (KGP) reveals numerous differences among the five KGP super-populations in the prevalence of individual and co-occurring HERV-K proviruses; we provide a visualization tool to easily depict the proportion of the KGP populations with any combination of polymorphic HERV-K provirus. Further, because HERV-K is insertionally polymorphic, the genome burden of known polymorphic HERV-K is variable in humans; this burden is lowest in East Asian (EAS) individuals. Our study identifies population-specific sequence variation for HERV-K proviruses at several loci. We expect these resources will advance research on HERV-K contributions to human diseases.

Human Endogenous Retrovirus type K (HERV-K) is the youngest of retrovirus families in the human genome and is the only group of endogenous retroviruses that has polymorphic members; a locus containing a HERV-K can be occupied in one individual but empty in others. HERV-Ks could contribute to disease risk or pathogenesis but linking one of the known polymorphic HERV-K to a specific disease has been difficult. We develop an easy to use method that reveals the considerable variation existing among global populations in the prevalence of individual and co-occurring polymorphic HERV-K, and in the number of HERV-K that any individual has in their genome. Our study provides a reference of diversity for the currently known polymorphic HERV-K in global populations and tools needed to determine the profile of all known polymorphic HERV-K in the genome of any patient population.

The integration of a macrophage-adapted live vaccine strain of equine infectious anaemia virus (EIAV) in the horse genome

Integration is an important feature of retroviruses and retrovirus-based therapeutic transfection vectors. The non-primate lentivirus equine infectious anaemia virus (EIAV) primarily targets macrophages/monocytes in vivo. Investigation of the integration features of EIAV_DLV121 strains, which are adapted to donkey monocyte-derived macrophages (MDMs), is of great interest. In this study, we analysed the integration features of EIAV_DLV121 in equine MDMs during in vitro infection. Our previously published integration sites (IS) for EIAV_FDDV13 in fetal equine dermal (FED) cells were also analysed in parallel as references. Sequencing of the host genomic regions flanking the viral IS showed that reference sequence (RefSeq) genes were preferentially targeted for integration by EIAV_DLV121. Introns, AT-rich regions, long interspersed nuclear elements (LINEs) and DNA transposons were also predominantly biased toward viral insertion, which is consistent with EIAV_FDDV13 integration into the horse genome in FED cells. In addition, the most significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, specifically gag junctions for EIAV_DLV121 and tight junctions for EIAV_FDDV13, are regulators of metabolic function, which is consistent with the common bioprocesses, specifically cell cycle and chromosome/DNA organization, identified by gene ontology (GO) analysis. Our results demonstrate that EIAV integration occurs in regions that harbour structural and topological features of local chromatin in both macrophages and fibroblasts. Our data on EIAV will facilitate further understanding of lentivirus infection and the development of safer and more effective gene therapy vectors.

Ub-ISAP: a streamlined UNIX pipeline for mining unique viral vector integration sites from next generation sequencing data

BACKGROUND

The analysis of viral vector genomic integration sites is an important component in assessing the safety and efficiency of patient treatment using gene therapy. Alongside this clinical application, integration site identification is a key step in the genetic mapping of viral elements in mutagenesis screens that aim to elucidate gene function.

RESULTS

We have developed a UNIX-based vector integration site analysis pipeline (Ub-ISAP) that utilises a UNIX-based workflow for automated integration site identification and annotation of both single and paired-end sequencing reads. Reads that contain viral sequences of interest are selected and aligned to the host genome, and unique integration sites are then classified as transcription start site-proximal, intragenic or intergenic.

CONCLUSION

Ub-ISAP provides a reliable and efficient pipeline to generate large datasets for assessing the safety and efficiency of integrating vectors in clinical settings, with broader applications in cancer research. Ub-ISAP is available as an open source software package at https://sourceforge.net/projects/ub-isap/ .

GCN2 phosphorylates HIV-1 integrase and decreases HIV-1 replication by limiting viral integration

GCN2 is a serine/threonine kinase involved in cellular stress response related to amino acid starvation. Previously, we showed that GCN2 interacts with HIV-1 integrase and is activated during HIV-1 infection. Herein, we identified HIV-1 integrase as a previously unknown substrate of GCN2 in vitro with a major site of phosphorylation at residue S255 located in the C-terminal domain of HIV-1 integrase. The underlying mechanism was investigated and it appeared that the integrase active site was required in order for GCN2 to target the integrase residue S255. Moreover, various integrases from other retroviruses (e.g. MLV, ASV) were also recognized as a substrate by GCN2. In cells, HIV-1 lentiviral particles harboring mutation at integrase position 255 were affected in their replication. Preventing phosphorylation resulted in an increase in infectivity that correlated with an increase in viral DNA integration. Infectivity of MLV was also higher in cells knocked-out for GCN2 suggesting a conserved mechanism to control viral replication. Altogether, our data suggest that GCN2 may constitute a general guardian of genome stability by regulating foreign DNA integration and as such be part of the antiviral armamentarium of the cell.

Integration site selection by retroviruses and transposable elements in eukaryotes

Transposable elements and retroviruses are found in most genomes, can be pathogenic and are widely used as gene-delivery and functional genomics tools. Exploring whether these genetic elements target specific genomic sites for integration and how this preference is achieved is crucial to our understanding of genome evolution, somatic genome plasticity in cancer and ageing, host-parasite interactions and genome engineering applications. High-throughput profiling of integration sites by next-generation sequencing, combined with large-scale genomic data mining and cellular or biochemical approaches, has revealed that the insertions are usually non-random. The DNA sequence, chromatin and nuclear context, and cellular proteins cooperate in guiding integration in eukaryotic genomes, leading to a remarkable diversity of insertion site distribution and evolutionary strategies.

Chromatin states shape insertion profiles of the piggyBac, Tol2 and Sleeping Beauty transposons and murine leukemia virus

DNA transposons and retroviruses are versatile tools in functional genomics and gene therapy. To facilitate their application, we conducted a genome-wide insertion site profiling of the piggyBac (PB), Tol2 and Sleeping Beauty (SB) transposons and the murine leukemia virus (MLV) in mouse embryonic stem cells (ESCs). PB and MLV preferred highly expressed genes, whereas Tol2 and SB preferred weakly expressed genes. However, correlations with DNase I hypersensitive sites were different for all vectors, indicating that chromatin accessibility is not the sole determinant. Therefore, we analysed various chromatin states. PB and MLV highly correlated with Cohesin, Mediator and ESC-specific transcription factors. Notably, CTCF sites were correlated with PB but not with MLV, suggesting MLV prefers smaller promoter-enhancer loops, whereas PB insertion encompasses larger chromatin loops termed topologically associating domains. Tol2 also correlated with Cohesin and CTCF. However, correlations with ESC-specific transcription factors were weaker, suggesting that Tol2 prefers transcriptionally weak chromatin loops. Consistently, Tol2 insertions were associated with bivalent histone modifications characteristic of silent and inducible loci. SB showed minimum preference to all chromatin states, suggesting the least adverse effect on adjacent genes. These results will be useful for vector selection for various applications.

INSPIIRED: A Pipeline for Quantitative Analysis of Sites of New DNA Integration in Cellular Genomes

Integration of new DNA into cellular genomes mediates replication of retroviruses and transposons; integration reactions have also been adapted for use in human gene therapy. Tracking the distributions of integration sites is important to characterize populations of transduced cells and to monitor potential outgrow of pathogenic cell clones. Here, we describe a pipeline for quantitative analysis of integration site distributions named INSPIIRED (integration site pipeline for paired-end reads). We describe optimized biochemical steps for site isolation using Illumina paired-end sequencing, including new technology for suppressing recovery of unwanted contaminants, then software for alignment, quality control, and management of integration site sequences. During library preparation, DNAs are broken by sonication, so that after ligation-mediated PCR the number of ligation junction sites can be used to infer abundance of gene-modified cells. We generated integration sites of known positions in silico, and we describe optimization of sample processing parameters refined by comparison to truth. We also present a novel graph-theory-based method for quantifying integration sites in repeated sequences, and we characterize the consequences using synthetic and experimental data. In an accompanying paper, we describe an additional set of statistical tools for data analysis and visualization. Software is available at https://github.com/BushmanLab/INSPIIRED.

Amplification, Next-generation Sequencing, and Genomic DNA Mapping of Retroviral Integration Sites

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of diverse libraries of retroviral integration sites using ligation-mediated PCR (LM-PCR) amplification and next-generation sequencing (NGS), (2) mapping the genomic location of each virus-host junction using BEDTools, and (3) analyzing the data for statistical relevance. Genomic DNA extracted from infected cells is fragmented by digestion with restriction enzymes or by sonication. After suitable DNA end-repair, double-stranded linkers are ligated onto the DNA ends, and semi-nested PCR is conducted using primers complementary to both the long terminal repeat (LTR) end of the virus and the ligated linker DNA. The PCR primers carry sequences required for DNA clustering during NGS, negating the requirement for separate adapter ligation. Quality control (QC) is conducted to assess DNA fragment size distribution and adapter DNA incorporation prior to NGS. Sequence output files are filtered for LTR-containing reads, and the sequences defining the LTR and the linker are cropped away. Trimmed host cell sequences are mapped to a reference genome using BLAT and are filtered for minimally 97% identity to a unique point in the reference genome. Unique integration sites are scrutinized for adjacent nucleotide (nt) sequence and distribution relative to various genomic features. Using this protocol, integration site libraries of high complexity can be constructed from genomic DNA in three days. The entire protocol that encompasses exogenous viral infection of susceptible tissue culture cells to integration site analysis can therefore be conducted in approximately one to two weeks. Recent applications of this technology pertain to longitudinal analysis of integration sites from HIV-infected patients.

Host Factors in Retroviral Integration and the Selection of Integration Target Sites

In order to replicate, a retrovirus must integrate a DNA copy of the viral RNA genome into a chromosome of the host cell. The study of retroviral integration has advanced considerably in the past few years. Here we focus on host factor interactions and the linked area of integration targeting. Genome-wide screens for cellular factors affecting HIV replication have identified a series of host cell proteins that may mediate subcellular trafficking for preintegration complexes, nuclear import, and integration target site selection. The cell transcriptional co-activator protein LEDGF/p75 has been identified as a tethering factor important for HIV integration, and recently, BET proteins (Brd2, 4, and 4) have been identified as tethering factors for the gammaretroviruses. A new class of HIV inhibitors has been developed targeting the HIV-1 IN-LEDGF binding site, though surprisingly these inhibitors appear to block assembly late during replication and do not act at the integration step. Going forward, genome-wide studies of HIV-host interactions offer many new starting points to investigate HIV replication and identify potential new inhibitor targets.

Characterization of Equine Infectious Anemia Virus Integration in the Horse Genome

Human immunodeficiency virus (HIV)-1 has a unique integration profile in the human genome relative to murine and avian retroviruses. Equine infectious anemia virus (EIAV) is another well-studied lentivirus that can also be used as a promising retro-transfection vector, but its integration into its native host has not been characterized. In this study, we mapped 477 integration sites of the EIAV strain EIAV_FDDV13 in fetal equine dermal (FED) cells during in vitro infection. Published integration sites of EIAV and HIV-1 in the human genome were also analyzed as references. Our results demonstrated that EIAV_FDDV13 tended to integrate into genes and AT-rich regions, and it avoided integrating into transcription start sites (TSS), which is consistent with EIAV and HIV-1 integration in the human genome. Notably, the integration of EIAV_FDDV13 favored long interspersed elements (LINEs) and DNA transposons in the horse genome, whereas the integration of HIV-1 favored short interspersed elements (SINEs) in the human genome. The chromosomal environment near LINEs or DNA transposons potentially influences viral transcription and may be related to the unique EIAV latency states in equids. The data on EIAV integration in its natural host will facilitate studies on lentiviral infection and lentivirus-based therapeutic vectors.

Coherence analysis discriminates between retroviral integration patterns in CD34(+) cells transduced under differing clinical trial conditions

Unequivocal demonstration of the therapeutic utility of γ-retroviral vectors for gene therapy applications targeting the hematopoietic system was accompanied by instances of insertional mutagenesis. These events stimulated the ongoing development of putatively safer integrating vector systems and analysis methods to characterize and compare integration site (IS) biosafety profiles. Continuing advances in next-generation sequencing technologies are driving the generation of ever-more complex IS datasets. Available bioinformatic tools to compare such datasets focus on the association of integration sites (ISs) with selected genomic and epigenetic features, and the choice of these features determines the ability to discriminate between datasets. We describe the scalable application of point-process coherence analysis (CA) to compare patterns produced by vector ISs across genomic intervals, uncoupled from association with genomic features. To explore the utility of CA in the context of an unresolved question, we asked whether the differing transduction conditions used in the initial Paris and London SCID-X1 gene therapy trials result in divergent genome-wide integration profiles. We tested a transduction carried out under each condition, and showed that CA could indeed resolve differences in IS distributions. Existence of these differences was confirmed by the application of established methods to compare integration datasets.

Bioinformatics and HIV latency

Despite effective treatment, HIV is not completely eliminated from the infected organism because of the existence of viral reservoirs. A major reservoir consists of infected resting CD4+ T cells, mostly of memory type, that persist over time due to the stable proviral insertion and a long cellular lifespan. Resting cells do not produce viral particles and are protected from viral-induced cytotoxicity or immune killing. However, these latently infected cells can be reactivated by stochastic events or by external stimuli. The present review focuses on novel genome-wide technologies applied to the study of integration, transcriptome, and proteome characteristics and their recent contribution to the understanding of HIV latency.

BET-independent MLV-based Vectors Target Away From Promoters and Regulatory Elements

Stable integration in the host genome renders murine leukemia virus (MLV)-derived vectors attractive tools for gene therapy. Adverse events in otherwise successful clinical trials caused by proto-oncogene activation due to vector integration hamper their application. MLV and MLV-based vectors integrate near strong enhancers, active promoters, and transcription start sites (TSS) through specific interaction of MLV integrase (IN) with the bromodomain and extra-terminal (BET) family of proteins, accounting for insertional mutagenesis. We identified a BET-interaction motif in the C-terminal tail of MLV IN conserved among gammaretroviruses. By deletion of this motif or a single point mutation (IN_W390A), BET-independent MLV (BinMLV) were engineered. BinMLV vectors carrying IN_W390A integrate at wild-type efficiency, with an integration profile that no longer correlates with BET chromatin distribution nor with the traditional markers of MLV integration. In particular, BinMLV vector integration associated less with oncogene TSS compared to the MLV vectors currently used in clinical trials. Together, these findings open perspectives to increase the biosafety of gammaretroviral vectors for gene therapy.

Allosteric inhibition of human immunodeficiency virus integrase: late block during viral replication and abnormal multimerization involving specific protein domains

HIV-1 replication in the presence of antiviral agents results in evolution of drug-resistant variants, motivating the search for additional drug classes. Here we report studies of GSK1264, which was identified as a compound that disrupts the interaction between HIV-1 integrase (IN) and the cellular factor lens epithelium-derived growth factor (LEDGF)/p75. GSK1264 displayed potent antiviral activity and was found to bind at the site occupied by LEDGF/p75 on IN by x-ray crystallography. Assays of HIV replication in the presence of GSK1264 showed only modest inhibition of the early infection steps and little effect on integration targeting, which is guided by the LEDGF/p75-IN interaction. In contrast, inhibition of late replication steps was more potent. Particle production was normal, but particles showed reduced infectivity. GSK1264 promoted aggregation of IN and preformed LEDGF/p75-IN complexes, suggesting a mechanism of inhibition. LEDGF/p75 was not displaced from IN during aggregation, indicating trapping of LEDGF/p75 in aggregates. Aggregation assays with truncated IN variants revealed that a construct with catalytic and C-terminal domains of IN only formed an open polymer associated with efficient drug-induced aggregation. These data suggest that the allosteric inhibitors of IN are promising antiviral agents and provide new information on their mechanism of action.

Lentiviral and targeted cellular barcoding reveals ongoing clonal dynamics of cell lines in vitro and in vivo

Background

Cell lines are often regarded as clonal, even though this simplifies what is known about mutagenesis, transformation and other processes that destabilize them over time. Monitoring these clonal dynamics is important for multiple areas of biomedical research, including stem cell and cancer biology. Tracking the contributions of individual cells to large populations, however, has been constrained by limitations in sensitivity and complexity.

Results

We utilize cellular barcoding methods to simultaneously track the clonal contributions of tens of thousands of cells. We demonstrate that even with optimal culturing conditions, common cell lines including HeLa, K562 and HEK-293 T exhibit ongoing clonal dynamics. Starting a population with a single clone diminishes but does not eradicate this phenomenon. Next, we compare lentiviral and zinc-finger nuclease barcode insertion approaches, finding that the zinc-finger nuclease protocol surprisingly results in reduced clonal diversity. We also document the expected reduction in clonal complexity when cells are challenged with genotoxic stress. Finally, we demonstrate that xenografts maintain clonal diversity to a greater extent than in vitro culturing of the human non-small-cell lung cancer cell line HCC827.

Conclusions

We demonstrate the feasibility of tracking and quantifying the clonal dynamics of entire cell populations within multiple cultured cell lines. Our results suggest that cell heterogeneity should be considered in the design and interpretation of in vitro culture experiments. Aside from clonal cell lines, we propose that cellular barcoding could prove valuable in modeling the clonal behavior of heterogeneous cell populations over time, including tumor populations treated with chemotherapeutic agents.

Integrase residues that determine nucleotide preferences at sites of HIV-1 integration: implications for the mechanism of target DNA binding

Retroviruses favor target-DNA (tDNA) distortion and particular bases at sites of integration, but the mechanism underlying HIV-1 selectivity is unknown. Crystal structures revealed a network of prototype foamy virus (PFV) integrase residues that distort tDNA: Ala188 and Arg329 interact with tDNA bases, while Arg362 contacts the phosphodiester backbone. HIV-1 integrase residues Ser119, Arg231, and Lys258 were identified here as analogs of PFV integrase residues Ala188, Arg329 and Arg362, respectively. Thirteen integrase mutations were analyzed for effects on integrase activity in vitro and during virus infection, yielding a total of 1610 unique HIV-1 integration sites. Purine (R)/pyrimidine (Y) dinucleotide sequence analysis revealed HIV-1 prefers the tDNA signature (0)RYXRY(4), which accordingly favors overlapping flexible dinucleotides at the center of the integration site. Consistent with roles for Arg231 and Lys258 in sequence specific and non-specific binding, respectively, the R231E mutation altered integration site nucleotide preferences while K258E had no effect. S119A and S119T integrase mutations significantly altered base preferences at positions −3 and 7 from the site of viral DNA joining. The S119A preference moreover mimicked wild-type PFV selectivity at these positions. We conclude that HIV-1 IN residue Ser119 and PFV IN residue Ala188 contact analogous tDNA bases to effect virus integration.

A large-scale zebrafish gene knockout resource for the genome-wide study of gene function

With the completion of the zebrafish genome sequencing project, it becomes possible to analyze the function of zebrafish genes in a systematic way. The first step in such an analysis is to inactivate each protein-coding gene by targeted or random mutation. Here we describe a streamlined pipeline using proviral insertions coupled with high-throughput sequencing and mapping technologies to widely mutagenize genes in the zebrafish genome. We also report the first 6144 mutagenized and archived F₁'s predicted to carry up to 3776 mutations in annotated genes. Using in vitro fertilization, we have rescued and characterized ∼0.5% of the predicted mutations, showing mutation efficacy and a variety of phenotypes relevant to both developmental processes and human genetic diseases. Mutagenized fish lines are being made freely available to the public through the Zebrafish International Resource Center. These fish lines establish an important milestone for zebrafish genetics research and should greatly facilitate systematic functional studies of the vertebrate genome.

Quantitation of HIV DNA integration: effects of differential integration site distributions on Alu-PCR assays

In many studies of HIV replication, it is useful to quantify the number of HIV proviruses in cells against a background of unintegrated forms of the HIV DNA. A popular method for doing so involves quantitative PCR using one primer complementary to the HIV long terminal repeat (LTR), and a second primer complementary to a cellular Alu repeat, so that PCR product only forms from templates where a provirus is integrated in the human genome near an Alu repeat. However, several recent studies have identified conditions that alter distributions of HIV integration sites relative to genes. Because Alu repeats are enriched in gene rich regions, this raises the question of whether altered integration site distributions might confound provirus abundance measurements using the Alu-PCR method. Here modified versions of the HIV tethering protein LEDGF/p75 were used to retarget HIV integration outside of transcription units, and show that this has a negligible effect on Alu-PCR quantitation of proviral abundance. Thus altered integration targeting, at least to the degree achieved here, is not a major concern when using the Alu-PCR assay.

rDNA-directed integration by an HIV-1 integrase--I-PpoI fusion protein

Integrating viral vectors are efficient gene transfer tools, but their integration patterns have been associated with genotoxicity and oncogenicity. The recent development of highly specific designer nucleases has enabled target DNA modification and site-specific gene insertion at desired genomic loci. However, a lack of consensus exists regarding a perfect genomic safe harbour (GSH) that would allow transgenes to be stably and reliably expressed without adversely affecting endogenous gene structure and function. Ribosomal DNA (rDNA) has many advantages as a GSH, but efficient means to target integration to this locus are currently lacking. We tested whether lentivirus vector integration can be directed to rDNA by using fusion proteins consisting of the Human Immunodeficiency Virus 1 (HIV-1) integrase (IN) and the homing endonuclease I-PpoI, which has natural cleavage sites in the rDNA. A point mutation (N119A) was introduced into I-PpoI to abolish unwanted DNA cleavage by the endonuclease. The vector-incorporated IN-I-PpoIN119A fusion protein targeted integration into rDNA significantly more than unmodified lentivirus vectors, with an efficiency of 2.7%. Our findings show that IN-fusion proteins can be used to modify the integration pattern of lentivirus vectors, and to package site-specific DNA-recognizing proteins into vectors to obtain safer transgene integration.

Decade-long safety and function of retroviral-modified chimeric antigen receptor T cells

The success of adoptive T cell gene transfer for treatment of cancer and HIV is predicated on generating a response that is both durable and safe. We report long-term results from three clinical trials to evaluate gammaretroviral vector-engineered T cells for HIV. The vector encoded a chimeric antigen receptor (CAR) composed of CD4 linked to the CD3ζ signaling chain (CD4ζ). CAR T cells were detected in 98% of samples tested for at least 11 years after infusion at frequencies that exceeded average T cell levels after most vaccine approaches. The CD4ζ transgene retained expression and function. There was no evidence of vector-induced immortalization of cells; integration site distributions showed no evidence of persistent clonal expansion or enrichment for integration sites near genes implicated in growth control or transformation. The CD4ζ T cells had stable levels of engraftment, with decay half-lives that exceeded 16 years, in marked contrast to previous trials testing engineered T cells. These findings indicate that host immunosuppression before T cell transfer is not required to achieve long-term persistence of gene-modified T cells. Further, our results emphasize the safety of T cells modified by retroviral gene transfer in clinical application, as measured in >500 patient-years of follow-up. Thus, previous safety issues with integrating viral vectors are hematopoietic stem cell or transgene intrinsic, and not a general feature of retroviral vectors. Engineered T cells are a promising form of synthetic biology for long-term delivery of protein-based therapeutics. These results provide a framework to guide the therapy of a wide spectrum of human diseases.

Integrase-independent HIV-1 infection is augmented under conditions of DNA damage and produces a viral reservoir

HIV-1 possesses a viral protein, integrase (IN), which is necessary for its efficient integration in target cells. However, it has been reported that an IN-defective HIV strain is still capable of integration. Here, we assessed the ability of wild type (WT) HIV-1 to establish infection in the presence of IN inhibitors. We observed a low, yet clear infection of inhibitor-incubated cells infected with WT HIV which was identical to cells infected with IN-deficient HIV, D64A. Furthermore, the IN-independent integration could be enhanced by the pretreatment of cells with DNA-damaging agents suggesting that integration is mediated by a DNA repair system. Moreover, significantly faster viral replication kinetics with augmented viral DNA integration was observed after infection in irradiated cells treated with IN inhibitor compared to nonirradiated cells. Altogether, our results suggest that HIV DNA has integration potential in the presence of an IN inhibitor and may serve as a virus reservoir.

Safe harbours for the integration of new DNA in the human genome

Interactions between newly integrated DNA and the host genome limit the reliability and safety of transgene integration for therapeutic cell engineering and other applications. Although targeted gene delivery has made considerable progress, the question of where to insert foreign sequences in the human genome to maximize safety and efficacy has received little attention. In this Opinion article, we discuss 'genomic safe harbours' - chromosomal locations where therapeutic transgenes can integrate and function in a predictable manner without perturbing endogenous gene activity and promoting cancer.

Gammaretroviral integration into nucleosomal target DNA in vivo

Some of the earliest studies of retroviral integration targeting reported that sites of gammaretroviral DNA integration were positively correlated with DNase I-hypersensitive sites in chromatin. This led to the suggestion that open chromatin was favorable for integration. More recent deep sequencing experiments confirmed that gammaretroviral integration sites and DNase I cleavage sites are associated in genome-wide surveys. Paradoxically, in vitro studies of integration show that nucleosomal DNA is actually favored over naked DNA, raising the question of whether integration target DNA in chromosomes is wrapped in nucleosomes or nucleosome free. In this study we examined gammaretroviral integration by infecting primary human CD4(+) T lymphocytes with a murine leukemia virus (MLV)-based retroviral vector or xenotropic murine leukemia virus-related virus (XMRV), and isolated 32,585 unique integration sites using ligation-mediated PCR and 454 pyrosequencing. CD4(+) T lymphocytes were chosen for study because of the particularly dense genome-wide mapping of chromatin features available for comparison. Analysis relative to predicted nucleosome positions showed that gammaretroviruses direct integration into outward-facing major grooves on nucleosome-wrapped DNA, similar to the integration pattern of HIV. Also, a suite of histone modifications correlated with gene activity are positively associated with integration by both MLV and XMRV. Thus, we conclude that favored integration near DNase I-hypersensitive sites does not imply that integration takes place exclusively in nucleosome-free regions.

Gammaretroviral integration into nucleosomal target DNA in vivo

Some of the earliest studies of retroviral integration targeting reported that sites of gammaretroviral DNA integration were positively correlated with DNase I-hypersensitive sites in chromatin. This led to the suggestion that open chromatin was favorable for integration. More recent deep sequencing experiments confirmed that gammaretroviral integration sites and DNase I cleavage sites are associated in genome-wide surveys. Paradoxically, in vitro studies of integration show that nucleosomal DNA is actually favored over naked DNA, raising the question of whether integration target DNA in chromosomes is wrapped in nucleosomes or nucleosome free. In this study we examined gammaretroviral integration by infecting primary human CD4(+) T lymphocytes with a murine leukemia virus (MLV)-based retroviral vector or xenotropic murine leukemia virus-related virus (XMRV), and isolated 32,585 unique integration sites using ligation-mediated PCR and 454 pyrosequencing. CD4(+) T lymphocytes were chosen for study because of the particularly dense genome-wide mapping of chromatin features available for comparison. Analysis relative to predicted nucleosome positions showed that gammaretroviruses direct integration into outward-facing major grooves on nucleosome-wrapped DNA, similar to the integration pattern of HIV. Also, a suite of histone modifications correlated with gene activity are positively associated with integration by both MLV and XMRV. Thus, we conclude that favored integration near DNase I-hypersensitive sites does not imply that integration takes place exclusively in nucleosome-free regions.

A method to sequence and quantify DNA integration for monitoring outcome in gene therapy

Human genetic diseases have been successfully corrected by integration of functional copies of the defective genes into human cells, but in some cases integration of therapeutic vectors has activated proto-oncogenes and contributed to leukemia. For this reason, extensive efforts have focused on analyzing integration site populations from patient samples, but the most commonly used methods for recovering newly integrated DNA suffer from severe recovery biases. Here, we show that a new method based on phage Mu transposition in vitro allows convenient and consistent recovery of integration site sequences in a form that can be analyzed directly using DNA barcoding and pyrosequencing. The method also allows simple estimation of the relative abundance of gene-modified cells from human gene therapy subjects, which has previously been lacking but is crucial for detecting expansion of cell clones that may be a prelude to adverse events.

Deciphering the code for retroviral integration target site selection

Upon cell invasion, retroviruses generate a DNA copy of their RNA genome and integrate retroviral cDNA within host chromosomal DNA. Integration occurs throughout the host cell genome, but target site selection is not random. Each subgroup of retrovirus is distinguished from the others by attraction to particular features on chromosomes. Despite extensive efforts to identify host factors that interact with retrovirion components or chromosome features predictive of integration, little is known about how integration sites are selected. We attempted to identify markers predictive of retroviral integration by exploiting Precision-Recall methods for extracting information from highly skewed datasets to derive robust and discriminating measures of association. ChIPSeq datasets for more than 60 factors were compared with 14 retroviral integration datasets. When compared with MLV, PERV or XMRV integration sites, strong association was observed with STAT1, acetylation of H3 and H4 at several positions, and methylation of H2AZ, H3K4, and K9. By combining peaks from ChIPSeq datasets, a supermarker was identified that localized within 2 kB of 75% of MLV proviruses and detected differences in integration preferences among different cell types. The supermarker predicted the likelihood of integration within specific chromosomal regions in a cell-type specific manner, yielding probabilities for integration into proto-oncogene LMO2 identical to experimentally determined values. The supermarker thus identifies chromosomal features highly favored for retroviral integration, provides clues to the mechanism by which retrovirus integration sites are selected, and offers a tool for predicting cell-type specific proto-oncogene activation by retroviruses.

When HIV-1, murine leukemia virus (MLV), or other retroviruses infect a cell, the virus generates a DNA copy of the viral RNA genome and ligates the cDNA within host chromosomal DNA. This integration reaction occurs at sites throughout the host cell genome, but little is known about how integration sites are selected. We attempted to identify markers predictive of retroviral integration by comparing the genome-wide binding sites for more than 60 factors with 14 retroviral integration datasets. We borrowed Precision-Recall methods from the Information Retrieval field for extracting information from highly skewed datasets such as these. For MLV and other gammaretroviruses, strong association was observed with STAT1, acetylation of H3 and H4 at several positions, and methylation of H2AZ, H3K4, and K9. We generated a supermarker by combining high scoring markers. The supermarker localized within 2 kB of 75% of MLV proviruses and predicted the likelihood of integration within specific chromosomal regions in a cell-type specific manner. This study identified chromosomal features highly favored for retroviral integration. It also provides clues to the mechanism by which retrovirus integration sites are selected, and offers a tool for predicting cell-type specific proto-oncogene activation by retroviruses.

Assessing the potential for AAV vector genotoxicity in a murine model

Gene transfer using adeno-associated virus (AAV) vectors has great potential for treating human disease. Recently, questions have arisen about the safety of AAV vectors, specifically, whether integration of vector DNA in transduced cell genomes promotes tumor formation. This study addresses these questions with high-dose liver-directed AAV-mediated gene transfer in the adult mouse as a model (80 AAV-injected mice and 52 controls). After 18 months of follow-up, AAV-injected mice did not show a significantly higher rate of hepatocellular carcinoma compared with controls. Tumors in mice treated with AAV vectors did not have significantly different amounts of vector DNA compared with adjacent normal tissue. A novel high-throughput method for identifying AAV vector integration sites was developed and used to clone 1029 integrants. Integration patterns in tumor tissue and adjacent normal tissue were similar to each other, showing preferences for active genes, cytosine-phosphate-guanosine islands, and guanosine/cytosine-rich regions. [corrected] Gene expression data showed that genes near integration sites did not show significant changes in expression patterns compared with genes more distal to integration sites. No integration events were identified as causing increased oncogene expression. Thus, we did not find evidence that AAV vectors cause insertional activation of oncogenes and subsequent tumor formation.

Identification of HIV integration sites in infected host genomic DNA

The integration of the Human Immunodeficiency Virus (HIV) genetic information into the host genome is fundamental for its replication and long-term persistence in the host. Isolating and characterizing the integration sites can be useful for obtaining data such as identifying the specific genomic location of integration or understanding the forces dictating HIV integration site selection. The methods outlined in this article describe a highly efficient and precise technique for identifying HIV integration sites in the host genome on a small scale using molecular cloning techniques and standard sequencing or on a massive scale using 454 pyrosequencing.

Retroviral integration site selection

The stable insertion of a copy of their genome into the host cell genome is an essential step of the life cycle of retroviruses. The site of viral DNA integration, mediated by the viral-encoded integrase enzyme, has important consequences for both the virus and the host cell. The analysis of retroviral integration site distribution was facilitated by the availability of the human genome sequence, revealing the non-random feature of integration site selection and identifying different favored and disfavored genomic locations for individual retroviruses. This review will summarize the current knowledge about retroviral differences in their integration site preferences as well as the mechanisms involved in this process.

Mechanistic and pharmacological analyses of HIV-1 integration

Significant advances have transpired in the human immunodeficiency virus type 1 (HIV-1) integration field in recent years. Considering its essential nature, integrase has long been a target of interest for antiviral drug development. The most significant advance was the approval of the Merck compound raltegravir, the first licensed integrase inhibitor, in October 2007. Another milestone was the identification and characterization of specific nucleoprotein complexes that mediate integrase 3' processing and DNA strand transfer activities in vitro. Genome-wide distribution analyses have furthermore revealed that different retroviruses differentially target distinctive regions of chromatin during integration. For examples, lentiviruses favor actively transcribed genes whereas gammaretroviruses such as Moloney murine leukemia virus prefer transcriptional start sites. Though the underlying mechanisms are unknown for most retroviruses, the lentiviral preference is in large part guided through the interaction with the integrase binding protein lens epithelium-derived growth factor (LEDGF)/p75. Experimental methods that formed the foundations for each of these advances, as well as other techniques topical to the study of HIV-1 integration, are described in this issue of Methods.