Personalized Cancer Monitoring Assay for the Detection of ctDNA in Patients with Solid Tumors

BACKGROUND

Highly sensitive molecular assays have been developed to detect plasma-based circulating tumor DNA (ctDNA), and emerging evidence suggests their clinical utility for monitoring minimal residual disease and recurrent disease, providing prognostic information, and monitoring therapy responses in patients with solid tumors. The Invitae Personalized Cancer Monitoring™ assay uses a patient-specific, tumor-informed variant signature identified through whole exome sequencing to detect ctDNA in peripheral blood of patients with solid tumors.

METHODS

The assay's tumor whole exome sequencing and ctDNA detection components were analytically validated using 250 unique human specimens and nine commercial reference samples that generated 1349 whole exome sequencing and cell-free DNA (cfDNA)-derived libraries. A comparison of tumor and germline whole exome sequencing was used to identify patient-specific tumor variant signatures and generate patient-specific panels, followed by targeted next-generation sequencing of plasma-derived cfDNA using the patient-specific panels with anchored multiplex polymerase chain reaction chemistry leveraging unique molecular identifiers.

RESULTS

Whole exome sequencing resulted in overall sensitivity of 99.8% and specificity of > 99.9%. Patient-specific panels were successfully designed for all 63 samples (100%) with ≥ 20% tumor content and 24 (80%) of 30 samples with ≥ 10% tumor content. Limit of blank studies using 30 histologically normal, formalin-fixed paraffin-embedded specimens resulted in 100% expected panel design failure. The ctDNA detection component demonstrated specificity of > 99.9% and sensitivity of 96.3% for a combination of 10 ng of cfDNA input, 0.008% allele frequency, 50 variants on the patient-specific panels, and a baseline threshold. Limit of detection ranged from 0.008% allele frequency when utilizing 60 ng of cfDNA input with 18-50 variants in the patient-specific panels (> 99.9% sensitivity) with a baseline threshold, to 0.05% allele frequency when using 10 ng of cfDNA input with an 18-variant panel with a monitoring threshold (> 99.9% sensitivity).

CONCLUSIONS

The Invitae Personalized Cancer Monitoring assay, featuring a flexible patient-specific panel design with 18-50 variants, demonstrated high sensitivity and specificity for detecting ctDNA at variant allele frequencies as low as 0.008%. This assay may support patient prognostic stratification, provide real-time data on therapy responses, and enable early detection of residual/recurrent disease.

Detection of Biallelic Loss of DNA Repair Genes in Formalin-Fixed, Paraffin-Embedded Tumor Samples Using a Novel Tumor-Only Sequencing Panel

Patient selection for synthetic lethal-based cancer therapy may be improved by assessment of gene-specific loss of heterozygosity (LOH) and biallelic loss of function (LOF). This report describes SyNthetic lethal Interactions for Precision Diagnostics (SNiPDx), a targeted next-generation sequencing (NGS) panel for detection of LOH and biallelic LOF alterations in 26 target genes focused on DNA damage response pathways, in tumor-only formalin-fixed, paraffin-embedded (FFPE) samples. NGS was performed across all exons of these 26 genes and encompassed a total of 7632 genome-wide single-nucleotide polymorphisms on genomic DNA from 80 FFPE solid tumor samples. The Fraction and Allele-Specific Copy Number Estimates from Tumor Sequencing algorithm was optimized to assess tumor purity and copy number based on heterozygous single-nucleotide polymorphisms. SNiPDx demonstrated high sensitivity (95%) and specificity (91%) for LOH detection compared with whole genome sequencing. Positive agreement with local NGS-based testing in the detection of genetic alterations was 95%. SNiPDx detected 93% of biallelic ATM LOF mutations, 100% of ATM single-nucleotide variants and small insertions/deletions, and 100% of all ATM LOH status events identified by orthogonal NGS-based testing. SNiPDx is a novel, clinically feasible test for analysis of allelic status in FFPE tumor samples, which demonstrated high accuracy when compared with other NGS-based approaches in clinical use.

Feasibility of blood testing combined with PET-CT to screen for cancer and guide intervention

Cancer treatments are often more successful when the disease is detected early. We evaluated the feasibility and safety of multicancer blood testing coupled with positron emission tomography-computed tomography (PET-CT) imaging to detect cancer in a prospective, interventional study of 10,006 women not previously known to have cancer. Positive blood tests were independently confirmed by a diagnostic PET-CT, which also localized the cancer. Twenty-six cancers were detected by blood testing. Of these, 15 underwent PET-CT imaging and nine (60%) were surgically excised. Twenty-four additional cancers were detected by standard-of-care screening and 46 by neither approach. One percent of participants underwent PET-CT imaging based on false-positive blood tests, and 0.22% underwent a futile invasive diagnostic procedure. These data demonstrate that multicancer blood testing combined with PET-CT can be safely incorporated into routine clinical care, in some cases leading to surgery with intent to cure.

Apolipoprotein L-1 renal risk variants form active channels at the plasma membrane driving cytotoxicity

Recently evolved alleles of Apolipoprotein L-1 (APOL1) provide increased protection against African trypanosome parasites while also significantly increasing the risk of developing kidney disease in humans. APOL1 protects against trypanosome infections by forming ion channels within the parasite, causing lysis. While the correlation to kidney disease is robust, there is little consensus concerning the underlying disease mechanism. We show in human cells that the APOL1 renal risk variants have a population of active channels at the plasma membrane, which results in an influx of both Na⁺ and Ca²⁺. We propose a model wherein APOL1 channel activity is the upstream event causing cell death, and that the activate-state, plasma membrane-localized channel represents the ideal drug target to combat APOL1-mediated kidney disease.

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.

INSPIIRED: Quantification and Visualization Tools for Analyzing Integration Site Distributions

Analysis of sites of newly integrated DNA in cellular genomes is important to several fields, but methods for analyzing and visualizing these datasets are still under development. Here, we describe tools for data analysis and visualization that take as input integration site data from our INSPIIRED pipeline. Paired-end sequencing allows inference of the numbers of transduced cells as well as the distributions of integration sites in target genomes. We present interactive heatmaps that allow comparison of distributions of integration sites to genomic features and that support numerous user-defined statistical tests. To summarize integration site data from human gene therapy samples, we developed a reproducible report format that catalogs sample population structure, longitudinal dynamics, and integration frequency near cancer-associated genes. We also introduce a novel summary statistic, the UC50 (unique cell progenitors contributing the most expanded 50% of progeny cell clones), which provides a single number summarizing possible clonal expansion. Using these tools, we characterize ongoing longitudinal characterization of a patient from the first trial to treat severe combined immunodeficiency-X1 (SCID-X1), showing successful reconstitution for 15 years accompanied by persistence of a cell clone with an integration site near the cancer-associated gene CCND2. Software is available at https://github.com/BushmanLab/INSPIIRED.

Structural and sequencing analysis of local target DNA recognition by MLV integrase

Target-site selection by retroviral integrase (IN) proteins profoundly affects viral pathogenesis. We describe the solution nuclear magnetic resonance structure of the Moloney murine leukemia virus IN (M-MLV) C-terminal domain (CTD) and a structural homology model of the catalytic core domain (CCD). In solution, the isolated MLV IN CTD adopts an SH3 domain fold flanked by a C-terminal unstructured tail. We generated a concordant MLV IN CCD structural model using SWISS-MODEL, MMM-tree and I-TASSER. Using the X-ray crystal structure of the prototype foamy virus IN target capture complex together with our MLV domain structures, residues within the CCD α2 helical region and the CTD β1-β2 loop were predicted to bind target DNA. The role of these residues was analyzed in vivo through point mutants and motif interchanges. Viable viruses with substitutions at the IN CCD α2 helical region and the CTD β1-β2 loop were tested for effects on integration target site selection. Next-generation sequencing and analysis of integration target sequences indicate that the CCD α2 helical region, in particular P187, interacts with the sequences distal to the scissile bonds whereas the CTD β1-β2 loop binds to residues proximal to it. These findings validate our structural model and disclose IN-DNA interactions relevant to target site selection.

Outcomes following gene therapy in patients with severe Wiskott-Aldrich syndrome

IMPORTANCE

Wiskott-Aldrich syndrome is a rare primary immunodeficiency associated with severe microthrombocytopenia. Partially HLA antigen-matched allogeneic hematopoietic stem cell (HSC) transplantation is often curative but is associated with significant comorbidity.

OBJECTIVE

To assess the outcomes and safety of autologous HSC gene therapy in Wiskott-Aldrich syndrome.

DESIGN, SETTING, AND PARTICIPANTS

Gene-corrected autologous HSCs were infused in 7 consecutive patients with severe Wiskott-Aldrich syndrome lacking HLA antigen-matched related or unrelated HSC donors (age range, 0.8-15.5 years; mean, 7 years) following myeloablative conditioning. Patients were enrolled in France and England and treated between December 2010 and January 2014. Follow-up of patients in this intermediate analysis ranged from 9 to 42 months.

INTERVENTION

A single infusion of gene-modified CD34+ cells with an advanced lentiviral vector.

MAIN OUTCOMES AND MEASURES

Primary outcomes were improvement at 24 months in eczema, frequency and severity of infections, bleeding tendency, and autoimmunity and reduction in disease-related days of hospitalization. Secondary outcomes were improvement in immunological and hematological characteristics and evidence of safety through vector integration analysis.

RESULTS

Six of the 7 patients were alive at the time of last follow-up (mean and median follow-up, 28 months and 27 months, respectively) and showed sustained clinical benefit. One patient died 7 months after treatment of preexisting drug-resistant herpes virus infection. Eczema and susceptibility to infections resolved in all 6 patients. Autoimmunity improved in 5 of 5 patients. No severe bleeding episodes were recorded after treatment, and at last follow-up, all 6 surviving patients were free of blood product support and thrombopoietic agonists. Hospitalization days were reduced from a median of 25 days during the 2 years before treatment to a median of 0 days during the 2 years after treatment. All 6 surviving patients exhibited high-level, stable engraftment of functionally corrected lymphoid cells. The degree of myeloid cell engraftment and of platelet reconstitution correlated with the dose of gene-corrected cells administered. No evidence of vector-related toxicity was observed clinically or by molecular analysis.

CONCLUSIONS AND RELEVANCE

This study demonstrated the feasibility of the use of gene therapy in patients with Wiskott-Aldrich syndrome. Controlled trials with larger numbers of patients are necessary to assess long-term outcomes and safety.

A modified γ-retrovirus vector for X-linked severe combined immunodeficiency

BACKGROUND

In previous clinical trials involving children with X-linked severe combined immunodeficiency (SCID-X1), a Moloney murine leukemia virus-based γ-retrovirus vector expressing interleukin-2 receptor γ-chain (γc) complementary DNA successfully restored immunity in most patients but resulted in vector-induced leukemia through enhancer-mediated mutagenesis in 25% of patients. We assessed the efficacy and safety of a self-inactivating retrovirus for the treatment of SCID-X1.

METHODS

We enrolled nine boys with SCID-X1 in parallel trials in Europe and the United States to evaluate treatment with a self-inactivating (SIN) γ-retrovirus vector containing deletions in viral enhancer sequences expressing γc (SIN-γc).

RESULTS

All patients received bone marrow-derived CD34+ cells transduced with the SIN-γc vector, without preparative conditioning. After 12.1 to 38.7 months of follow-up, eight of the nine children were still alive. One patient died from an overwhelming adenoviral infection before reconstitution with genetically modified T cells. Of the remaining eight patients, seven had recovery of peripheral-blood T cells that were functional and led to resolution of infections. The patients remained healthy thereafter. The kinetics of CD3+ T-cell recovery was not significantly different from that observed in previous trials. Assessment of insertion sites in peripheral blood from patients in the current trial as compared with those in previous trials revealed significantly less clustering of insertion sites within LMO2, MECOM, and other lymphoid proto-oncogenes in our patients.

CONCLUSIONS

This modified γ-retrovirus vector was found to retain efficacy in the treatment of SCID-X1. The long-term effect of this therapy on leukemogenesis remains unknown. (Funded by the National Institutes of Health and others; ClinicalTrials.gov numbers, NCT01410019, NCT01175239, and NCT01129544.).

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.

A new class of multimerization selective inhibitors of HIV-1 integrase

The quinoline-based allosteric HIV-1 integrase (IN) inhibitors (ALLINIs) are promising candidates for clinically useful antiviral agents. Studies using these compounds have highlighted the role of IN in both early and late stages of virus replication. However, dissecting the exact mechanism of action of the quinoline-based ALLINIs has been complicated by the multifunctional nature of these inhibitors because they both inhibit IN binding with its cofactor LEDGF/p75 and promote aberrant IN multimerization with similar potencies in vitro. Here we report design of small molecules that allowed us to probe the role of HIV-1 IN multimerization independently from IN-LEDGF/p75 interactions in infected cells. We altered the rigid quinoline moiety in ALLINIs and designed pyridine-based molecules with a rotatable single bond to allow these compounds to bridge between interacting IN subunits optimally and promote oligomerization. The most potent pyridine-based inhibitor, KF116, potently (EC₅₀ of 0.024 µM) blocked HIV-1 replication by inducing aberrant IN multimerization in virus particles, whereas it was not effective when added to target cells. Furthermore, KF116 inhibited the HIV-1 IN variant with the A128T substitution, which confers resistance to the majority of quinoline-based ALLINIs. A genome-wide HIV-1 integration site analysis demonstrated that addition of KF116 to target or producer cells did not affect LEDGF/p75-dependent HIV-1 integration in host chromosomes, indicating that this compound is not detectably inhibiting IN-LEDGF/p75 binding. These findings delineate the significance of correctly ordered IN structure for HIV-1 particle morphogenesis and demonstrate feasibility of exploiting IN multimerization as a therapeutic target. Furthermore, pyridine-based compounds present a novel class of multimerization selective IN inhibitors as investigational probes for HIV-1 molecular biology.

The administration of highly active-antiretroviral therapy (HAART) has changed what was once a terminal disease into a manageable chronic infection. The success of HAART is manifested by reduced mortality and morbidity of HIV-1 infected patients. However, evolution of HIV-1 strains resistant to current therapies is a major clinical problem in the fight against AIDS. Therefore, new inhibitors with novel mechanisms of action are needed. One such mechanism is to target multimerization of HIV-1 integrase. In the present study, we report the design of pyridine-based small molecules that contain a rotatable single bond to allow optimal bridging between interacting integrase subunits. As a result, pyridine-based compounds stabilized interacting IN subunits and promoted aberrant, higher order integrase multimerization. The most potent compound, KF116, potently inhibited HIV-1 replication by interfering with proper maturation of HIV-1 particles, whereas KF116 at therapeutically relevant (submicromolar) concentrations had no detectable effects on LEDGF/p75 mediated HIV-1 integration. Our findings highlight HIV-1 integrase multimerization as a plausible therapeutic target and offer a path for designing improved inhibitors for potential clinical use.

Altering murine leukemia virus integration through disruption of the integrase and BET protein family interaction

We report alterations to the murine leukemia virus (MLV) integrase (IN) protein that successfully result in decreasing its integration frequency at transcription start sites and CpG islands, thereby reducing the potential for insertional activation. The host bromo and extraterminal (BET) proteins Brd2, 3 and 4 interact with the MLV IN protein primarily through the BET protein ET domain. Using solution NMR, protein interaction studies, and next generation sequencing, we show that the C-terminal tail peptide region of MLV IN is important for the interaction with BET proteins and that disruption of this interaction through truncation mutations affects the global targeting profile of MLV vectors. The use of the unstructured tails of gammaretroviral INs to direct association with complexes at active promoters parallels that used by histones and RNA polymerase II. Viruses bearing MLV IN C-terminal truncations can provide new avenues to improve the safety profile of gammaretroviral vectors for human gene therapy.

Bimodal high-affinity association of Brd4 with murine leukemia virus integrase and mononucleosomes

The importance of understanding the molecular mechanisms of murine leukemia virus (MLV) integration into host chromatin is highlighted by the development of MLV-based vectors for human gene-therapy. We have recently identified BET proteins (Brd2, 3 and 4) as the main cellular binding partners of MLV integrase (IN) and demonstrated their significance for effective MLV integration at transcription start sites. Here we show that recombinant Brd4, a representative of the three BET proteins, establishes complementary high-affinity interactions with MLV IN and mononucleosomes (MNs). Brd4(1–720) but not its N- or C-terminal fragments effectively stimulate MLV IN strand transfer activities in vitro. Mass spectrometry- and NMR-based approaches have enabled us to map key interacting interfaces between the C-terminal domain of BRD4 and the C-terminal tail of MLV IN. Additionally, the N-terminal fragment of Brd4 binds to both DNA and acetylated histone peptides, allowing it to bind tightly to MNs. Comparative analyses of the distributions of various histone marks along chromatin revealed significant positive correlations between H3- and H4-acetylated histones, BET protein-binding sites and MLV-integration sites. Our findings reveal a bimodal mechanism for BET protein-mediated MLV integration into select chromatin locations.

The major cellular sterol regulatory pathway is required for Andes virus infection

The Bunyaviridae comprise a large family of RNA viruses with worldwide distribution and includes the pathogenic New World hantavirus, Andes virus (ANDV). Host factors needed for hantavirus entry remain largely enigmatic and therapeutics are unavailable. To identify cellular requirements for ANDV infection, we performed two parallel genetic screens. Analysis of a large library of insertionally mutagenized human haploid cells and a siRNA genomic screen converged on components (SREBP-2, SCAP, S1P and S2P) of the sterol regulatory pathway as critically important for infection by ANDV. The significance of this pathway was confirmed using functionally deficient cells, TALEN-mediated gene disruption, RNA interference and pharmacologic inhibition. Disruption of sterol regulatory complex function impaired ANDV internalization without affecting virus binding. Pharmacologic manipulation of cholesterol levels demonstrated that ANDV entry is sensitive to changes in cellular cholesterol and raises the possibility that clinically approved regulators of sterol synthesis may prove useful for combating ANDV infection.

As obligate, intracellular parasites viruses are dependent upon the host cell for numerous factors and processes. However, for many important viruses few of the required host factors have been identified. Hantaviruses are rodent-borne viruses that are associated with severe human disease. Transmission to humans occurs sporadically with a recent notable example in Yosemite National park. In the present study, we utilized two independent genetic strategies to discover cellular factors needed for replication of the highly pathogenic hantavirus Andes virus. We found that four genes, encoding components of a complex involved in regulation of cholesterol synthesis and uptake, were critical for Andes virus infection. Drugs that inhibit an enzyme in this complex or that reduce cellular cholesterol levels effectively blocked Andes virus infection, suggesting new ways for combating this pathogenic virus.

Comparing DNA integration site clusters with scan statistics

MOTIVATION

Gene therapy with retroviral vectors can induce adverse effects when those vectors integrate in sensitive genomic regions. Retroviral vectors are preferred that target sensitive regions less frequently, motivating the search for localized clusters of integration sites and comparison of the clusters formed by integration of different vectors. Scan statistics allow the discovery of spatial differences in clustering and calculation of false discovery rates providing statistical methods for comparing retroviral vectors.

RESULTS

A scan statistic for comparing two vectors using multiple window widths is proposed with software to detect clustering differentials and compute false discovery rates. Application to several sets of experimentally determined HIV integration sites demonstrates the software. Simulated datasets of various sizes and signal strengths are used to determine the power to discover clusters and evaluate a convenient lower bound. This provides a toolkit for planning evaluations of new gene therapy vectors.

AVAILABILITY AND IMPLEMENTATION

The geneRxCluster R package containing a simple tutorial and usage hints is available from http://www.bioconductor.org.

Recombinant adeno-associated virus integration sites in murine liver after ornithine transcarbamylase gene correction

Recombinant adeno-associated viruses (rAAVs) have been tested in humans and other large mammals without adverse events. However, one study of mucopolysaccharidosis VII correction in mice showed repeated integration of rAAV in cells from hepatocellular carcinoma (HCC) in the Dlk1-Dio3 locus, suggesting possible insertional mutagenesis. In contrast, another study found no association of rAAV integration with HCC, raising questions about the generality of associations between liver transformation and integration at Dlk1-Dio3. Here we report that in rAAV-treated ornithine transcarbamylase (Otc)-deficient mice, four examples of integration sites in Dlk1-Dio3 could be detected in specimens from liver nodule/tumors, confirming previous studies of rAAV integration in the Dlk1-Dio3 locus in the setting of another murine model of metabolic disease. In one case, the integrated vector was verified to be present at about one copy per cell, consistent with clonal expansion. Another verified integration site in liver nodule/tumor tissue near the Tax1bp1 gene was also detected at about one copy per cell. The Dlk1-Dio3 region has also been implicated in human HCC and so warrants careful monitoring in ongoing human clinical trials with rAAV vectors.

The BET family of proteins targets moloney murine leukemia virus integration near transcription start sites

A hallmark of retroviral replication is integration of the viral genome into host cell DNA. This characteristic makes retrovirus-based vectors attractive delivery vehicles for gene therapy. However, adverse events in gene therapeutic trials, caused by activation of proto-oncogenes due to murine leukemia virus (MLV)-derived vector integration, hamper their application. Here, we show that bromodomain and extraterminal (BET) proteins (BRD2, BRD3, and BRD4) and MLV integrase specifically interact and colocalize within the nucleus of the cell. Inhibition of the BET proteins' chromatin interaction via specific bromodomain inhibitors blocks MLV virus replication at the integration step. MLV integration site distribution parallels the chromatin binding profile of BET proteins, and expression of an artificial fusion protein of the BET integrase binding domain with the chromatin interaction domain of the lentiviral targeting factor LEDGF/p75 retargets MLV integration away from transcription start sites and into the body of actively transcribed genes, conforming to the HIV integration pattern. Together, these data validate BET proteins as MLV integration targeting factors.

HIV latency and integration site placement in five cell-based models

Background

HIV infection can be treated effectively with antiretroviral agents, but the persistence of a latent reservoir of integrated proviruses prevents eradication of HIV from infected individuals. The chromosomal environment of integrated proviruses has been proposed to influence HIV latency, but the determinants of transcriptional repression have not been fully clarified, and it is unclear whether the same molecular mechanisms drive latency in different cell culture models.

Results

Here we compare data from five different in vitro models of latency based on primary human T cells or a T cell line. Cells were infected in vitro and separated into fractions containing proviruses that were either expressed or silent/inducible, and integration site populations sequenced from each. We compared the locations of 6,252 expressed proviruses to those of 6,184 silent/inducible proviruses with respect to 140 forms of genomic annotation, many analyzed over chromosomal intervals of multiple lengths. A regularized logistic regression model linking proviral expression status to genomic features revealed no predictors of latency that performed better than chance, though several genomic features were significantly associated with proviral expression in individual models. Proviruses in the same chromosomal region did tend to share the same expressed or silent/inducible status if they were from the same cell culture model, but not if they were from different models.

Conclusions

The silent/inducible phenotype appears to be associated with chromosomal position, but the molecular basis is not fully clarified and may differ among in vitro models of latency.

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.

Comparative analysis of the recently discovered hAT transposon TcBuster in human cells

Background

Transposons are useful tools for creating transgenic organisms, insertional mutagenesis, and genome engineering. TcBuster, a novel hAT-family transposon system derived from the red flour beetle Tribolium castaneum, was shown to be highly active in previous studies in insect embryoes.

Methodology/Principal Findings

We tested TcBuster for its activity in human embryonic kidney 293 (HEK-293) cells. Excision footprints obtained from HEK-293 cells contained small insertions and deletions consistent with a hAT-type repair mechanism of hairpin formation and non-homologous end-joining. Genome-wide analysis of 23,417 piggyBac, 30,303 Sleeping Beauty, and 27,985 TcBuster integrations in HEK-293 cells revealed a uniquely different integration pattern when compared to other transposon systems with regards to genomic elements. TcBuster experimental conditions were optimized to assay TcBuster activity in HEK-293 cells by colony assay selection for a neomycin-containing transposon. Increasing transposon plasmid increased the number of colonies, whereas gene transfer activity dependent on codon-optimized transposase plasmid peaked at 100 ng with decreased colonies at the highest doses of transposase DNA. Expression of the related human proteins Buster1, Buster3, and SCAND3 in HEK-293 cells did not result in genomic integration of the TcBuster transposon. TcBuster, Tol2, and piggyBac were compared directly at different ratios of transposon to transposase and found to be approximately comparable while having their own ratio preferences.

Conclusions/Significance

TcBuster was found to be highly active in mammalian HEK-293 cells and represents a promising tool for mammalian genome engineering.

HRP-2 determines HIV-1 integration site selection in LEDGF/p75 depleted cells

Background

Lens epithelium–derived growth factor (LEDGF/p75) is a cellular co-factor of HIV-1 integrase (IN) that tethers the viral pre-integration complex to the host cell chromatin and determines the genome wide integration site distribution pattern of HIV-1. Recently, we demonstrated that HIV-1 replication was reduced in LEDGF/p75 knockout (KO) cells. LEDGF/p75 KO significantly altered the integration site preference of HIV-1, but the pattern remained distinct from a computationally generated matched random control set (MRC), suggesting the presence of an alternative tethering factor. We previously identified Hepatoma-derived growth factor related protein 2 (HRP-2) as a factor mediating LEDGF/p75-independent HIV-1 replication. However, the role of HRP-2 in HIV-1 integration site selection was not addressed.

Findings

We studied the HIV-1 integration site distribution in the presence and absence of LEDGF/p75 and/or HRP-2, and in LEDGF/p75-depleted cells that overexpress HRP-2. We show that HRP-2 functions as a co-factor of HIV-1 IN in LEDGF/p75-depleted cells. Endogenous HRP-2 only weakly supported HIV-1 replication in LEDGF/p75 depleted cells. However, HRP-2 overexpression rescued HIV-1 replication and restored integration in RefSeq genes to wild-type levels. Additional HRP-2 KD in LEDGF/p75-depleted cells reduces integration frequency in transcription units and shifts the integration distribution towards random.

Conclusions

We demonstrate that HRP-2 overexpression can compensate for the absence of LEDGF/p75 and indicate that the residual bias in integration targeting observed in the absence of LEDGF/p75 can be ascribed to HRP-2. Knockdown of HRP-2 upon LEDGF/p75 depletion results in a more random HIV-1 integration pattern. These data therefore reinforce the understanding that LEDGF/p75 is the dominant HIV-1 IN co-factor.

Role of the PWWP domain of lens epithelium-derived growth factor (LEDGF)/p75 cofactor in lentiviral integration targeting

LEDGF/p75 is a chromatin-interacting, cellular cofactor of HIV integrase that dictates lentiviral integration site preference. In this study we determined the role of the PWWP domain of LEDGF/p75 in tethering and targeting of the lentiviral pre-integration complex, employing potent knockdown cell lines allowing analysis in the absence of endogenous LEDGF/p75. Deletion of the PWWP domain resulted in a diffuse subnuclear distribution pattern, loss of interaction with condensed chromatin, and failure to rescue proviral integration, integration site distribution, and productive virus replication. Substitution of the PWWP domain of LEDGF/p75 with that of hepatoma-derived growth factor or HDGF-related protein-2 rescued viral replication and lentiviral integration site distribution in LEDGF/p75-depleted cells. Replacing all chromatin binding elements of LEDGF/p75 with full-length hepatoma-derived growth factor resulted in more integration in genes combined with a preference for CpG islands. In addition, we showed that any PWWP domain targets SMYD1-like sequences. Analysis of integration preferences of lentiviral vectors for epigenetic marks indicates that the PWWP domain is critical for interactions specifying the relationship of integration sites to regions enriched in specific histone post-translational modifications.

A gene-rich, transcriptionally active environment and the pre-deposition of repressive marks are predictive of susceptibility to KRAB/KAP1-mediated silencing

Background

KRAB-ZFPs (Krüppel-associated box domain-zinc finger proteins) are vertebrate-restricted transcriptional repressors encoded in the hundreds by the mouse and human genomes. They act via an essential cofactor, KAP1, which recruits effectors responsible for the formation of facultative heterochromatin. We have recently shown that KRAB/KAP1 can mediate long-range transcriptional repression through heterochromatin spreading, but also demonstrated that this process is at times countered by endogenous influences.

Method

To investigate this issue further we used an ectopic KRAB-based repressor. This system allowed us to tether KRAB/KAP1 to hundreds of euchromatic sites within genes, and to record its impact on gene expression. We then correlated this KRAB/KAP1-mediated transcriptional effect to pre-existing genomic and chromatin structures to identify specific characteristics making a gene susceptible to repression.

Results

We found that genes that were susceptible to KRAB/KAP1-mediated silencing carried higher levels of repressive histone marks both at the promoter and over the transcribed region than genes that were insensitive. In parallel, we found a high enrichment in euchromatic marks within both the close and more distant environment of these genes.

Conclusion

Together, these data indicate that high levels of gene activity in the genomic environment and the pre-deposition of repressive histone marks within a gene increase its susceptibility to KRAB/KAP1-mediated repression.

Abstract 2715: Host genomic environment determines HTLV-1 clone size in vivo

HTLV-1 persists by driving clonal proliferation of infected T-lymphocytes. A high proviral load predisposes to the inflammatory and malignant diseases associated with HTLV-1. Yet the reasons for the remarkable variation within and between individuals in the size of HTLV-1-infected clones remain unknown. We devised a high-throughput protocol to map the genomic location and quantify the abundance of >91,000 distinct integration sites from 61 HTLV-1+ individuals and >2,100 sites from in vitro infection. We show that a typical HTLV-1-infected host carries between 500 and 5000 distinct integrated proviral clones. We demonstrate that negative selection dominates during chronic infection, favouring survival of proviruses integrated in transcriptionally silenced DNA: this selection is significantly stronger in asymptomatic carriers. The high proviral load characteristic of HTLV-1-associated inflammatory disease results from a larger number of integrated proviruses and not as previously believed from a difference in clonality. The abundance of surviving integrated proviruses is determined by genomic features of the flanking DNA and the transcriptional orientation of the provirus relative to the nearest host gene. HTLV-1 clonal expansion in vivo is favoured by orientation of the provirus in the same sense as the nearest host gene.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 2715. doi:10.1158/1538-7445.AM2011-2715

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.

HIV integration targeting: a pathway involving Transportin-3 and the nuclear pore protein RanBP2

Genome-wide siRNA screens have identified host cell factors important for efficient HIV infection, among which are nuclear pore proteins such as RanBP2/Nup358 and the karyopherin Transportin-3/TNPO3. Analysis of the roles of these proteins in the HIV replication cycle suggested that correct trafficking through the pore may facilitate the subsequent integration step. Here we present data for coupling between these steps by demonstrating that depletion of Transportin-3 or RanBP2 altered the terminal step in early HIV replication, the selection of chromosomal sites for integration. We found that depletion of Transportin-3 and RanBP2 altered integration targeting for HIV. These knockdowns reduced HIV integration frequency in gene-dense regions and near gene-associated features, a pattern that differed from that reported for depletion of the HIV integrase binding cofactor Psip1/Ledgf/p75. MLV integration was not affected by the Transportin-3 knockdown. Using siRNA knockdowns and integration targeting analysis, we also implicated several additional nuclear proteins in proper target site selection. To map viral determinants of integration targeting, we analyzed a chimeric HIV derivative containing MLV gag, and found that the gag replacement phenocopied the Transportin-3 and RanBP2 knockdowns. Thus, our data support a model in which Gag-dependent engagement of the proper transport and nuclear pore machinery mediate trafficking of HIV complexes to sites of integration.

HIV continues to be responsible for approximately two million deaths worldwide each year. As part of the viral replication cycle, the viral cDNA is transported through the nuclear pore into the nucleus where it integrates into the host cell genome. HIV integrates non-randomly, likely choosing integration sites within the host chromosomes that best enable the viral genes to be expressed and, ultimately, progeny virus to be produced. HIV uses host factors to guide its selection of integration sites. Here we demonstrate that components of the nuclear trafficking and nuclear pore machinery are required for HIV to achieve its normal pattern of integration sites. This finding suggests that passage of the virus through the nuclear pore into the nucleus is coupled to downstream integration events and enables the virus to achieve its final position within the host genome. Our study provides new insights into two important steps of the HIV replication cycle and suggests possible new targets for anti-retroviral drugs.

High-definition mapping of retroviral integration sites defines the fate of allogeneic T cells after donor lymphocyte infusion

The infusion of donor lymphocytes transduced with a retroviral vector expressing the HSV-TK suicide gene in patients undergoing hematopoietic stem cell transplantation for leukemia/lymphoma promotes immune reconstitution and prevents infections and graft-versus-host disease. Analysis of the clonal dynamics of genetically modified lymphocytes in vivo is of crucial importance to understand the potential genotoxic risk of this therapeutic approach. We used linear amplification-mediated PCR and pyrosequencing to build a genome-wide, high-definition map of retroviral integration sites in the genome of peripheral blood T cells from two different donors and used gene expression profiling and bioinformatics to associate integration clusters to transcriptional activity and to genetic and epigenetic features of the T cell genome. Comparison with matched random controls and with integrations obtained from CD34⁺ hematopoietic stem/progenitor cells showed that integration clusters occur within chromatin regions bearing epigenetic marks associated with active promoters and regulatory elements in a cell-specific fashion. Analysis of integration sites in T cells obtained ex vivo two months after infusion showed no evidence of integration-related clonal expansion or dominance, but rather loss of cells harboring integration events interfering with RNA post-transcriptional processing. The study shows that high-definition maps of retroviral integration sites are a powerful tool to analyze the fate of genetically modified T cells in patients and the biological consequences of retroviral transduction.

Genomic safe harbors permit high β-globin transgene expression in thalassemia induced pluripotent stem cells

Realizing the therapeutic potential of human induced pluripotent stem (iPS) cells will require robust, precise and safe strategies for genetic modification, as cell therapies that rely on randomly integrated transgenes pose oncogenic risks. Here we describe a strategy to genetically modify human iPS cells at 'safe harbor' sites in the genome, which fulfill five criteria based on their position relative to contiguous coding genes, microRNAs and ultraconserved regions. We demonstrate that ∼10% of integrations of a lentivirally encoded β-globin transgene in β-thalassemia-patient iPS cell clones meet our safe harbor criteria and permit high-level β-globin expression upon erythroid differentiation without perturbation of neighboring gene expression. This approach, combining bioinformatics and functional analyses, should be broadly applicable to introducing therapeutic or suicide genes into patient-specific iPS cells for use in cell therapy.

Oxidation of Aliphatic Aldehydes by Tetraethylammonium Chlorochromate: A Kinetic and Mechanistic Study

The oxidation of six aliphatic aldehydes by tetraethylammonium chlorochromate (TEACC) in dimethyl sulfoxide (DMSO) leads to the formation of the corresponding carboxylic acids. The reaction is first order each in TEACC and the aldehyde. The reaction is catalysed by hydrogen ions. The hydrogen–ion dependence has the form: kobs = a + b[H+]. The oxidation of deuteriated acetaldehyde, MeCDO, exhibited a substantial primary kinetic isotope effect ( kH/ kD = 5.83 at 298 K). The oxidation of acetaldehyde has been studied in 19 different organic solvents. The solvent effect has been analysed using Taft's and Swain's multiparametric equations. The rate constants correlate well with Taft's σ* values; reaction constants being negative. A mechanism involving transfer of hydride ion has been suggested.

Methods for integration site distribution analyses in animal cell genomes

The question of where retroviral DNA becomes integrated in chromosomes is important for understanding (i) the mechanisms of viral growth, (ii) devising new anti-retroviral therapy, (iii) understanding how genomes evolve, and (iv) developing safer methods for gene therapy. With the completion of genome sequences for many organisms, it has become possible to study integration targeting by cloning and sequencing large numbers of host-virus DNA junctions, then mapping the host DNA segments back onto the genomic sequence. This allows statistical analysis of the distribution of integration sites relative to the myriad types of genomic features that are also being mapped onto the sequence scaffold. Here we present methods for recovering and analyzing integration site sequences.

Integration target site selection by a resurrected human endogenous retrovirus

At least 8% of the human genome was formed by integration of retroviral DNA sequences. Here we analyze the forces directing the accumulation of human endogenous retroviruses (HERVs) by comparing de novo HERV integration targeting with the distribution of fixed HERV elements in the human genome. All known genomic HERVs are inactive due to mutation, but we were able to study integration targeting using a reconstituted consensus HERV-K (designated HERV-K(Con)). We found that HERV-K(Con) integrated preferentially in transcription units, in gene-rich regions, and near features associated with active transcription units and associated regulatory regions. In contrast, genomic HERV-K proviruses are found preferentially outside transcription units. The minority of genomic HERVKs present inside transcription units are in opposite transcriptional orientation relative to the host gene, the orientation predicted to be minimally disruptive to host mRNA synthesis, but de novo HERV-K(Con) integration within transcription units showed no orientation bias. We also found that the youngest HERV-K elements in the human genome showed a distribution intermediate between de novo HERV-K(Con) integration sites and older fixed HERV-Ks. These findings indicate that accumulation of HERVs in the human germline is a two-step process: integration targeting biases direct initial accumulation, then purifying selection leads to loss of proviruses disrupting gene function.

Analysis of lentiviral vector integration in HIV+ study subjects receiving autologous infusions of gene modified CD4+ T cells

Lentiviral vector-based gene therapy has been used to target the human immunodeficiency virus (HIV) using an antisense env payload. We have analyzed lentiviral-vector integration sites from three treated individuals. We compared integration sites from the ex vivo vector-transduced CD4+ cell products to sites from cells recovered at several times after infusion. Integration sites were analyzed using 454 pyrosequencing, yielding a total of 7,782 unique integration sites from the ex vivo product and 237 unique sites from cells recovered after infusion. Integrated vector copies in both data sets were found to be strongly enriched within active genes and near epigenetic marks associated with active transcription units. Analysis of integration relative to nucleosome structure on target DNA indicated favoring of integration in outward facing DNA major grooves on the nucleosome surface. There was no indication that growth of transduced cells after infusion resulted in enrichment for integration sites near proto-oncogene 5'-ends or within tumor suppressor genes. Thus, this first look at the longitudinal evolution of cells transduced with a lentiviral vector after infusion of gene modified CD4+ cells provided no evidence for abnormal expansions of cells due to vector-mediated insertional activation of proto-oncogenes.

Bacteriophage Mu integration in yeast and mammalian genomes

Genomic parasites have evolved distinctive lifestyles to optimize replication in the context of the genomes they inhabit. Here, we introduced new DNA into eukaryotic cells using bacteriophage Mu DNA transposition complexes, termed 'transpososomes'. Following electroporation of transpososomes and selection for marker gene expression, efficient integration was verified in yeast, mouse and human genomes. Although Mu has evolved in prokaryotes, strong biases were seen in the target site distributions in eukaryotic genomes, and these biases differed between yeast and mammals. In Saccharomyces cerevisiae transposons accumulated outside of genes, consistent with selection against gene disruption. In mouse and human cells, transposons accumulated within genes, which previous work suggests is a favorable location for efficient expression of selectable markers. Naturally occurring transposons and viruses in yeast and mammals show related, but more extreme, targeting biases, suggesting that they are responding to the same pressures. These data help clarify the constraints exerted by genome structure on genomic parasites, and illustrate the wide utility of the Mu transpososome technology for gene transfer in eukaryotic cells.