Disrupting MLV integrase:BET protein interaction biases integration into quiescent chromatin and delays but does not eliminate tumor activation in a MYC/Runx2 mouse model
Authors:
Lorenz Loyola aff001; Vasudevan Achuthan aff002; Kathryn Gilroy aff003; Gillian Borland aff004; Anna Kilbey aff004; Nancy Mackay aff004; Margaret Bell aff005; Jodie Hay aff004; Sriram Aiyer aff001; Dylan Fingerman aff001; Rodrigo A. Villanueva aff001; Ewan Cameron aff005; Christine A. Kozak aff006; Alan N. Engelman aff002; James Neil aff004; Monica J. Roth aff001
Authors place of work:
Rutgers-Robert Wood Johnson Medical School, Dept of Pharmacology, Piscataway, New Jersey, United States of America
aff001; Department of Cancer Immunology and Virology, Dana-Farber Cancer Institute, Boston, Massachusetts, United States of America
aff002; Beatson Institute for Cancer Research, Institute of Cancer Sciences, University of Glasgow, Glasgow, United Kingdom
aff003; MRC Univ. of Glasgow Centre for Virus Research, College of Medicine, Veterinary Medicine and Life Sciences, University of Glasgow, Glasgow, United Kingdom
aff004; Univ. of Glasgow School of Veterinary Medicine, Department of Veterinary Pathology Bearsden, United Kingdom
aff005; NIH, NIAID, Bethesda, Maryland, United States of America
aff006; Harvard Medical School, Department of Medicine, Boston, Massachusetts, United States of America
aff007
Published in the journal:
Disrupting MLV integrase:BET protein interaction biases integration into quiescent chromatin and delays but does not eliminate tumor activation in a MYC/Runx2 mouse model. PLoS Pathog 15(12): e32767. doi:10.1371/journal.ppat.1008154
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1008154
Summary
Murine leukemia virus (MLV) integrase (IN) lacking the C-terminal tail peptide (TP) loses its interaction with the host bromodomain and extraterminal (BET) proteins and displays decreased integration at promoter/enhancers and transcriptional start sites/CpG islands. MLV lacking the IN TP via an altered open reading frame was used to infect tumorigenesis mouse model (MYC/Runx2) animals to observe integration patterns and phenotypic effects, but viral passage resulted in the restoration of the IN TP through small deletions. Mice subsequently infected with an MLV IN lacking the TP coding sequence (TP-) showed an improved median survival by 15 days compared to wild type (WT) MLV infection. Recombination with polytropic endogenous retrovirus (ERV), Pmv20, was identified in seven mice displaying both fast and slow tumorigenesis, highlighting the strong selection within the mouse to maintain the full-length IN protein. Mapping the genomic locations of MLV in tumors from an infected mouse with no observed recombination with ERVs, TP-16, showed fewer integrations at TSS and CpG islands, compared to integrations observed in WT tumors. However, this mouse succumbed to the tumor in relatively rapid fashion (34 days). Analysis of the top copy number integrants in the TP-16 tumor revealed their proximity to known MLV common insertion sites genes while maintaining the MLV IN TP- genotype. Furthermore, integration mapping in K562 cells revealed an insertion preference of MLV IN TP- within chromatin profile states associated with weakly transcribed heterochromatin with fewer integrations at histone marks associated with BET proteins (H3K4me1/2/3, and H3K27Ac). While MLV IN TP- showed a decreased overall rate of tumorigenesis compared to WT virus in the MYC/Runx2 model, MLV integration still occurred at regions associated with oncogenic driver genes independently from the influence of BET proteins, either stochastically or through trans-complementation by functional endogenous Gag-Pol protein.
Keywords:
Chromatin – DNA methylation – Mammalian genomics – Protein interactions – Genomic libraries – Polymerase chain reaction – Mouse models – carcinogenesis
Introduction
Integration is an essential step for retroviral replication and pathogenesis (for review [1]). The integration substrate generated through reverse transcription is linear viral DNA containing a copy of the long terminal repeat (LTR) at each end. Integration is driven by the viral integrase (IN), structurally defined by distinct functional domains and regions, including the N-terminal region (NTR) containing the N-terminal extension (NED) and HHCC zinc-binding domain (NTD) [2], the catalytic core domain (CCD) [3], and the C-terminal domain (CTD) containing an unstructured tail peptide (TP) [4, 5]. Integration proceeds via two distinct IN catalytic activities; 3’ processing and strand transfer. During processing of the gammaretrovirus murine leukemia virus (MLV) LTR, a TT dinucleotide is hydrolyzed, exposing CAOH 3’ ends. IN then uses the 3’ hydroxyl groups to cut the chromosomal DNA in a staggered fashion, which joins the viral DNA 3’ ends to the 5’ phosphates of the host DNA cut. Repair of the gapped integration intermediate, which contains unjoined viral DNA 5’ ends with protruding 5’-AA single stranded (ss) DNA tails, yields a 4 bp duplication of host chromosomal DNA flanking the integrated MLV provirus [6, 7].
Host protein interactions are key determinants of integration preferences of retroviruses [8–11]. MLV displays preferential integration near regions of high transcriptional activity, such as promoter and enhancer regions and transcription start sites (TSSs) [12, 13]. The interaction of the MLV IN with the host bromo- and extraterminal (BET) domain proteins influences target-site selection. BET proteins interact with various chromatin remodeler proteins [14], therefore guiding MLV integration to highly transcriptionally active chromatin regions. MLV IN has an unstructured TP within the CTD that interacts with the extraterminal (ET) domain of BET proteins [4, 5]. Removal of this interaction by substitution or truncation of the TP reduces preferential integration at TSS and CpG islands, which can redistribute the integration profile and decrease the oncogenic effects of MLV insertional mutagenesis [4, 15].
In mice, Moloney MLV (M-MLV) is a non-acute retrovirus and thus insertional activation of proto-oncogenes at identified common insertions sites (CISs) is the predominant mechanism of oncogenesis [1, 16], requiring a long-latency period varying between 4–12 months [16]. A transgenic mouse overexpressing two MLV CIS genes, MYC and Runx2, from CD2 promoters exhibits early onset lymphomagenesis through a synergistic mechanism that is proposed to neutralize p53 activation [17, 18]. Infection of MYC/Runx2 mice with WT MLV reduced animal survival by 10 days [19–21]. Additionally, neonatal infection of this mouse model with M-MLV WT virus accelerated tumorigenesis and increased clonal complexity through various insertional mutagenesis sites [20]. Analysis of these integration sites through next-generation sequencing and subsequent comparison with reference genomes and ChIP-seq data sets identified a panel of MLV CIS that accelerated the oncogenic process [20]. MYC/Runx2 mice are accordingly an established model to study the relationship between M-MLV integration at predetermined CISs and tumorigenesis [20, 21].
Murine gammaretroviruses are classified based on their exogenous versus endogenous localization and receptor usage [22, 23]. Inbred strains of mice harbor endogenous type C MLVs, which fall into three general classes depending on their receptor usage, and thus their host and tissue specificities. These classes include the ecotropic viruses, limited to rodents (mCAT1 receptor), xenotropic viruses (excluded from infection of inbred mice [22]; Xpr1 phosphate exporter receptor), and polytropic/mixed polytropic viruses [24], infecting mouse and nonrodent species [22, 23] using Xpr1 as their receptor [25]. Although many endogenous ecotropic and some xenotropic viruses can form infectious particles, the endogenous polytropic MLVs (P-MLVs) do not produce replication competent viruses [22]. However, these sequences are an abundant source for recombination, when challenged with alternative defective or replication competent viruses. The generation of such recombinants frequently results in viruses with improved virulence and the exchange of the viral env gene [26–28]. Of significance to this study, C57BL/10 mice express xenotropic MLV from the Bxv1 locus, which can be a source of viral proteins as well as genetic material [29].
MLV-based vectors were used in initial human gene therapy trials [30] as well as selected CAR-T cell therapies [31]. In multiple clinical trials, including X-linked SCID [32–34], X-linked chronic granulomatous disease [35] and Wiskott-Aldrich Syndrome [36], but not ADA deficiencies [30], insertional mutagenesis resulted in the outgrowth of oligoclonal populations due to trans-activation of proto-oncogenes [37]. Subsequent approaches involving self-inactivating (SIN) vectors [38] or lentiviral vectors have been used to address these outcomes [30]. Alternatively, addressing the integration target-site bias of gammaretroviruses to integrate preferentially at promoter/enhancer regions by altering or eliminating their interaction with host BET protein could alter the oncogenic potential of these vectors [4, 39].
In this study, the genotoxicity of replication-competent M-MLV lacking the IN region required for interaction with the host BET proteins (IN TP-) was directly examined using the MYC/Runx2 transgenic mouse model [20, 21]. Analysis of tumor progression within MYC/Runx2 mice infected neonatally with MLV that maintain the IN TP- genotype as well as the integration preferences of such mutants in human K562 cells highlights the preferential integration into quiescent states and the strong selective pressure on MLV to maintain the IN tail peptide, through either internal deletions or recombination with endogenous retroviruses. The impact of our findings on the development of MLV-based vectors for human gene therapy is discussed.
Results
Previously, in vitro studies indicated that the integration bias of MLV IN lacking the tail peptide towards CpG islands and TSS was reduced due to the loss of interaction with host BET proteins [4]. In this study, the effects of virus lacking the BET interaction domain within the IN TP were examined using the mouse MYC/Runx2 model. This model is advantageous because lymphomas form within 36 days, and lymphomagenesis is accelerated by a further 10 days following insertional mutagenesis by M-MLV at known CIS.
Generation and characterization of replication competent M-MLV lacking the IN TP
The MLV IN TP region and the env coding sequence overlap in alternative open reading frames. In order to analyze the effects of an M-MLV bearing IN C-terminal truncations on the MYC/Runx2 mouse model, replication competent viruses that terminated the IN protein without altering the expression of the ecotropic M-MLV Env were required.
Studies utilizing MLV IN-XN
Initial experiments were performed using the previously reported pNCA-C IN-XN construct [40], which truncated the C-terminal 23 aa of IN while maintaining a viral titer within 2 fold of WT M-MLV in 293T cells [4]. As illustrated in Fig 1A, the IN-XN construct introduces a stop codon within IN, upstream of the env coding region, resulting in a frame-shift of the sequence encoding the IN C-terminal region. Survival of MYC/Runx2 mice infected with WT M-MLV or MLV IN-XN in comparison to an uninfected control was monitored by Date of Death (DoD) over a 115-day period following neonatal injection. The survival curves of WT (n = 30) and IN-XN (n = 40) infected mice were statistically indistinguishable (P = 0.089) using the Log-rank test but both showed significant differences using the same test to the uninfected control (****P<0.0001, **P<0.0021 respectively) (Fig 1F). Median survival times of WT and IN-XN infected mice were 35 and 43 days, respectively, and both succumbed to tumors significantly faster than the uninfected control (median survival 54 days). To investigate further, IN-XN virus isolated from the tumors of three infected mice was introduced into 293mCAT human cells that express the mouse ecotropic receptor [41], which facilitated the isolation of infectious virus in the absence of endogenous mouse viruses. Remarkably, the viruses transferred to 293mCAT cells from two independent mice (XN3; DoD 30d and XN35; DoD 70d) harbored deletions of 20 and 5 bases, respectively (Fig 1B). These deletions counteracted the IN stop codon and restored the C-terminus, encoded in an alternate reading frame, onto the IN protein (Fig 1C). Functionally, such deletions alter the spacing between the IN C-terminal domain (CTD) SH3 fold and the region of the C-terminus that becomes structured upon binding to the host Brd ET domain (Fig 1C) [41]. DNA from the XN2 mouse (DoD = 41d) maintained the IN-XN genotype. From this we noted that reversion to the WT TP sequence did not necessarily correlate with early-onset DoD.
Studies utilizing MLV IN TP-
In order to address the complication of small deletions restoring the TP to the IN protein, a second construct that codon optimized the env reading frame, pNCA-C IN TP-, was generated. This construct design eliminated the coding potential of the IN TP within the env overlapping region and incorporated multiple stop codons into the non-env reading frames, thus eliminating the potential to restore the WT IN sequence through simple deletion (Fig 1A). Single-round infection of D17 cells with MLV IN TP- confirmed that the viral titer was equivalent to MLV IN-XN; in this assay, a 10-fold decrease compared with the wild-type (WT) MLV was observed (Fig 1D). However, highlighting the infectious potential of both MLV IN-XN and TP- viruses, their titers were 3-orders of magnitude greater than the D184N IN mutant virus that carried the Asp>Asn substitution in the enzyme active site [41]. We were also interested in examining the effects of three amino acid substitutions in the IN CCD (E266A, L268A, and Y269A) that were previously reported to disrupt MLV IN interaction with Brd2 proteins using co-immunoprecipitation experiments [42]. However, the viral titers of MLV bearing IN E266A, IN L268A, and IN Y269A were equivalent to the catalytically inactive IN D184N mutant (Fig 1D). Thus, tissue culture and animal studies with these viruses were not possible.
Viral passage of IN-XN and TP- viruses in D17/pJET cells, which express the ecotropic mCAT receptor [43], displayed similar infection time courses, with viral capsid proteins (CA) detected in the medium on day 5. This is a delay of two days compared to the WT MLV (Fig 1E). Competency of the IN TP- virus was further confirmed by western blot analysis of the viral CA, IN, Env proteins released from D17/pJET producer cell line (S1 Fig), with both Env (74 kDa) and CA (30 kDa) detected in the media of cells infected with either MLV IN TP- or WT MLV virus. The truncation of the IN TP was stable after passage in tissue culture (S1 Fig, lane 1), and the predicted decreased molecular weight compared with WT IN was maintained (S1 Fig, lane 2). As expected, no viral proteins (Env, CA or IN) were detected for the replication-defective pNCA-C IN D184N virus (S1 Fig, lane 3).
Having verified that the MLV IN TP- virus was replication competent and that the truncation was stable, survival experiments were performed in MYC/Runx2 mice to determine if this optimized construct affected tumorigenesis (Fig 1G and Fig 2). As expected [19–21] (Fig 1F), mice infected with WT MLV exhibited significantly poorer survival than the uninfected controls (Log-rank test: ****P<0.0001; (Fig 1G)), while the lifespan of MLV IN TP- infected mice was extended compared to WT MLV (P = 0.006) (Fig 1G). Interestingly, the survival curve of MLV IN TP- infected mice (n = 23) showed no significant differences (P = 0.26) from that of the uninfected mice (n = 25), in contrast to the observed curve for the mice infected with MLV IN-XN virus (Fig 1F). Of note, the survival curve for IN TP- infected mice showed a biphasic trend with one third of the mice developing tumors comparatively early. Additionally, all of the IN TP- mice were deceased by day 72, while some uninfected mice survived until day 115. The median survival times for untreated mice as well as those infected with WT MLV or MLV IN TP- were 53, 35, and 50 days, respectively. These observed phenotypic differences suggest that infection with MLV IN TP- delayed tumorigenesis compared to WT MLV in this mouse model.
Analysis of integration site preferences in MYC/Runx2 mice infected with WT or TP- MLV
Previous studies showed that truncation of IN TP decreased preferential integration at TSS and CpG islands by >50% in tissue culture [4, 15], and experiments were performed to determine if this was also the case in the MYC/Runx2 mice. Fig 2 outlines the overall MYC/Runx2 mouse model (Fig 2A) and workflow for analysis of MLV integration sites (Fig 2B), which were bioinformatically mapped from genomic DNA samples from MYC/Runx2 mouse tumors following ligation-mediated PCR (LM-PCR) and next generation sequencing (NGS) (Table 1). Initially, libraries were generated from tumors from 4 representative mice infected with WT MLV (WT6, 8, 10, and 12), 5 representative mice infected with MLV IN TP- (TP-4, 6, 7, 9, and 16), as well as one non-infected control (NC) mouse (Table 1). Integration sites were initially mapped to +/- 1 kb of TSSs and CpG islands, as described previously [44]. TSS and CpG island-proximal integration frequencies were 11.8–13.4% and 11.9–14.3%, respectively, for WT tumors. Surprisingly, tumors from 4 of the 5 MYC/Runx2 mice infected with MLV IN TP- (TP-4, 6, 7, and 9) showed similar preferential integration at TSSs and CpG islands to WT MLV (10.2–11.9% and 11.3–13.4% respectively). In contrast, MLV integrations from TP-16 tumor at TSSs and CpG islands were markedly lower compared to other tumors from IN TP- and WT MLV-infected mice. MLV integration sites proximal to TSSs and CpG islands from all tumor samples were statistically different compared to the NC control (Table in S1 Table). MLV integrations from TP-16 were statistically different from WT6, even though each animal had the same DoD (Day 34, Table in S2 Table) (Fisher’s test, P< 0.001, Table in S1 Table). Furthermore, the percentage of MLV TP-16 integrations at TSSs (7.3%) and CpG islands (6.3%) was respectively1.7- and 2.0-fold lower than the mean from all WT MLV tumors at these markers (TSSs = 12.4% ± 0.35 and CpG = 12.9% ± 0.57). While the integration profiles of MLV TP-16 are consistent with the loss of the IN TP region, insertions from tumors TP-4, 6, 7, and 9 were more similar to the WT virus. Unlike the IN-XN construct, the IN TP construct cannot revert to functional IN through deletions in the viral genome; however, recombination with endogenous retroviruses could have restored the IN TP in TP-4, 6, 7, and 9 tumors.
Detection of M-MLV recombination with ERVs
To verify the structure of the virus that persisted in TP-16 tumors and to investigate the cause of WT-like behavior of TP-4, 6, 7, and 9 tumors, the samples were analyzed for recombination with ERVs (Fig 2A). Recombinants were detected using PCR with primer pairs that included a primer for each of two different classes of known ERVs (polytropic and xenotropic) and an M-MLV primer (S2 Fig, Table in S3 Table, [24]). Amphotropic Env primers were included in the analysis as a control for amplification of laboratory constructs [45].
Analysis of recombinant MLV through infection of 293mCAT cells
To facilitate this analysis, infectious virus isolated from tumor cells was used to infect human 293mCAT cells that express the mouse ecotropic virus receptor. Env recombinants with polytropic and xenotropic MLV would also be infectious on this cell line. Transferring the virus to 293mCAT cells eliminates the potential for background amplification products from ERVs in mouse cells. PCR analyses performed on viruses derived from TP-4, 6, 7, 9, 16 tumors using primers specified in Table in S3 Table are shown in Fig 3 (corresponding DNA gels in S3 Fig). Viral PCR products were detected for TP-6, 7, and 9 using the RT_universal_fwd primer and two independent reverse polytropic specific primers (Polytropic_JS4_rev and Polytropic_JS5_rev) (S3A and S3B Fig). The absence of mouse contamination in the 293mCAT DNA was confirmed with negative PCR results for intracisternal particle A (IAP) and mouse mitochondrial cyclooxygenase-2 (mCOX2) DNA sequences (S3C and S3D Fig). Only TP-9 yielded a PCR product using the RT_universal_fwd primer paired with either the Amphotropic_rev or the Xenotropic_JS10_rev, but the quality and quantity of this product was insufficient for subsequent analysis.
For these studies, it is of interest to define the crossover junctions within recombinants, as these directly address the potential restoration of the IN TP and receptor recognition of the subsequent virus. Thus, 3’ and 5’ sequence junctions of M-MLV/polytopic ERV recombinants were determined. Polytropic ERV (P-ERV) DNA segments from TP-6, 7, 9 had close homology to the P-ERV Pmv20 [46, 47]. The 5’ junction points for the recombinants were found to be within the IN region, with those of TP-6 and TP-9 within the IN CCD (Fig 3A and 3B, region 1, green) and that of TP-7 within the IN CTD (Fig 3A and 3B, region 2, blue). The 3’ junction points for all corresponded to the C-terminus end of the surface (SU) subunit of Env glycoprotein. In summary, the segments that recombined into the M-MLV TP- virus contain the WT IN TP sequence, resulting in the restoration of expression of full-length IN during viral spread in the mice.
Analysis of recombinant MLV in mouse tumors
Having shown that recombination had occurred in some of the mice infected with MLV IN TP, tumors from other mice with early DoD were investigated to see if recombination events had also occurred. Tumor DNA samples from MLV IN TP-4, 12, 13, 15, 16, 17, 18, and 19 (Table in S2 Table) were thus analyzed by nested PCR. The primary PCR product was generated using a forward primer that hybridized to all MLV classes (RT_universal_fwd) and a reverse M-MLV primer (MLV_LTR_U3_rev) that amplifies 3’ M-MLV DNA sequences including IN and Env. For the second round, polytropic recombinants were detected using a Polytropic_JS5_fwd primer and an M-MLV reverse primer (7791_rev) (Fig 3A, bottom). Recombinant viruses were detected in tumor samples TP-4, 15, 18 and 19, and were successfully sequenced to identify the 3’ junction points (Fig 3B). However, the 5’ junction breakpoints were not identified in these recombinants. As above, the polytropic segment of these recombinants was homologous to Pmv20. The 3’ breakpoints for the TP- recombinants varied throughout the C-terminus of the Env SU (Fig 3B (region 3; bottom)). The 3’ breakpoints of TP-9 and TP-18 could not be distinguished in this analysis.
The recombinants identified within the IN (Region 1) had crossover junctions distinct from that of the MCF247 and C58v2 recombinant MLVs, with TP-6 having partial overlap with the previously identified crossover region of M965 MLV (Fig 3B, black boxes) [46]. Similarly, within Region 3, the crossover identified in TP-7 partially overlapped that previously described for PTV-1 [46]. Notably, the presence of recombination did not correlate with tumorigenesis in the MYC/Runx2 mouse. For example, TP-16 with a DoD of 34 days had no detectable recombinants in any of its DNA samples, while TP-4, 6 and 9 all had recombinant virus and DoDs of 50, 57 and 63 days, respectively (Table in S2 Table).
The integration profile of MLV IN TP-16 tumor is distinct from that of WT MLV
Analysis of the MLV IN TP- genome within TP-16 tumor DNA did not detect recombination with ERVs and indicated decreased viral integrations at TSSs and CpG islands (Table 1). We therefore sought to determine detailed profiles of WT and TP- MLV integration in tumors that arose with similar kinetics, as evidenced by the same 34 day DoD of WT6 and TP-16 animals (Table in S2 Table; Fig 4). As expected, WT MLV integrated symmetrically around TSSs (black, Fig 4A left) compared to the NC sample (orange, Fig 4A right), with the majority of the integrants located within 1 kb of TSSs. TP-16 distribution of integrants around TSSs was asymmetric and more dispersed compared to WT6 (Fig 4A blue, center), with comparatively increased integration events located ±8 kb from TSSs. Comparison of the MLV integration profile in these tumors with Brd4 binding sites (data taken from ENCODE ID GSM1262345, murine AML MLL-AF9/NrasG12D cells) indicates that MLV IN TP-16 integrations at Brd4 sites was decreased by approximately 20% compared to MLV WT6 (Fig 4B). Additionally, approximately 50% fewer TP-16 integrations were observed +/- 1 kb from the histone modifications H3K4me1 and H3K4me3 associated with active chromatin when compared to WT6 (ENCODE IDs ENCSR000CCI and ENCSR000CCJ, respectively), and were at a levels that were comparable with the NC (Fig 4C). Statistical analysis of the WT and TP-16 integrants with respect to Brd4, H3K4me1 and H3K4me 3 is presented in Table in S1 Table. H3K4me1 and H3K4me3 are considered MLV supermarkers [48], with H3K4me3 associated with nucleosome-bound BET protein [49].
In addition to the histone marks studied above, BET proteins associate significantly with active enhancer features, specifically acetylated histone tails. In the human lymphoma cell line Ly1 DLBCL, 79.1% of H3K27Ac sites overlap with Brd4 and 92.2% of chromatin bound Brd4 is at regions of the H3K27Ac active enhancer mark [50]. MLV integrations are also reported to be highly enriched at H3K27Ac sites [5, 12, 42, 51]. 31.4% of MLV integration sites from WT6 tumors were within H3K27Ac sites (ENCODE ID ENCFF001KYG) (Fig 4D), whereas only 12.8% of the MLV IN TP-16 integration sites overlapped with H3K27Ac sites. Furthermore, all of the WT samples analyzed had similar overlap with H3K27Ac sites (31.2% for WT8, 35.3% for WT10, 37.1% for WT12) (Table in S4 Table). The distribution of the distance between retroviral integration sites (RISs) and the nearest H3K27Ac peak was significantly broader for TP-16 than for the cumulative RISs from all WT tumors (Wilcoxon Rank Sum, p<2.2 x 10−16; Fig 4E).
Viruses extracted from MLV IN TP-16 tumor were used to infect 293mCAT cells and integration site analysis was performed after 14 days. The integration profile in 293mCAT for MLV IN TP-16 paralleled that of MLV IN-XN (Table 2). Overall, these integration profiles indicate that the viral population from the TP-16 tumor maintained the MLV IN TP- genotype within the mouse and that the integrants are distributed further away from promoter and active enhancers compared to WT MLV.
Integration sites within the MLV IN TP-16 tumors localize to known MLV integration common insertion sites (CISs)
Integration site copy number is representative of predominant integrants that clonally expand during tumorigenesis [20]. Table in S5 Table summarizes the key MLV integrants present in the IN TP-16 tumor. The top 25 CISs for WT MLV in the MYC/Runx2 mouse have previously been identified [20]. Within the IN TP-16 tumor, integrants localized at 11 of the 25 CIS (S5 Table). Significantly, of the top ten copy number integrants, seven were within these previously identified loci (Table in S5 Table). The viral integrants from the TP-16 tumors at three of these CIS genes, Mapk13, Ccnd1, Hdac6, is shown schematically in Fig 5A & 5B. Insertions at all three genes showed an orientation bias for the genomic (+) strand (mm10). For Ccnd1 and Rasgrp1, this orientation bias is antisense to the host gene transcription, which is consistent with enhancer insertional activation reported previously [52–54]. Interestingly, the Rasgrp1 cluster of integrations mapped >91 kb upstream of the Rasgrp1 promoter, whereas for Ccnd1, the high copy number insert was located 791 bp from the promoter (black triangle, Fig 5A). It is striking that three independent insertions in Hdac6 were highly abundant, and all were within the gene: the first disrupted exon 3 and corresponded to the N-terminus of the protein, while two abundant insertions localized within introns spanning the gene (introns 5 and 28) (Table in S5 Table).
Integrants at CISs encode the IN TP- genotype
It was of considerable interest to verify that the MLV inserted into these three genes maintained the IN TP- genotype. Based on the known insertion sites, a nested PCR was developed to amplify the IN region of pol with specific host sequences (Fig 5C and Table in S5 Table). Sequencing of the resulting PCR products verified that all four insertions at Hdac6 (2 insertions), Ccnd1, and Rasgrp1 maintained the parental IN TP- mutation.
Mapping integrations of MLV IN TP- in K562 cells for 15 chromatin states and histone modifications
To further investigate the pattern of MLV IN TP- integration, a single round of infection of leukemia cell line K562 with this virus was performed. As expected, these cells recapitulated the decreased integration percentages at TSSs and CpG islands that was previously observed in 293mCAT cells for the IN-XN construct (Tables 2 and 3, Fig 6A). As previously reported [13], the use of chromatin states provides a different approach to understand the genome landscape. In this approach, clusters of chromatin marks are used to define functional active states of chromosomes, which are specific for each cell line. For K562, 15 chromosome states have been utilized to analyze MLV integration [13], and 85% of the integrants were shown to map to strong enhancer and active promoter regions. The overlap between MLV WT and TP- integration sites and the components of the 15-chromatin state model was investigated in K562 cells. As shown in Fig 6B, 73.8% of MLV WT integrations mapped to the same three highest states that were identified previously [13], which were annotated as active promoter (state 1) and strong enhancer (states 4 and 5). In contrast, MLV IN TP- integrations displayed a divergent preference, with the top two chromatin states being heterochromatin (state 13, (21.6%)) and weakly transcribed region (state 11, 19.4%) (Fig 6B, bottom). Loss of the IN TP decreased, but did not eliminate, integrations at active promoters (state 1) and enhancers (strong enhancers, states 4 and 5, and weak enhancers, state 7), which cumulatively accounted for 35% of the IN TP- integrants and are drastically reduced compared to WT MLV (82%).
The median chromosomal coverage in K562 for chromatin states 11 and 13 were reported to be 11.3 and 71.4% [55], respectively, which correlates with the observed integrations from the coverage of the computer-generated random integration control (RIC) within these states (Fig 6C). For state 11, corresponding to weakly transcribed regions, loss of the IN TP increased integration frequency 4.8-fold compared to WT (Fig 6C; 19.4% for IN TP-; 4.01% for WT MLV). Significantly, this frequency is 2-fold above the RIC, indicating a bias for integration into these weakly transcribed regions (Table S1 Table). Similarly, integration into heterochromatin (state 13), increased ~ 5-fold in MLV IN TP- (21.63%) compared to WT MLV (4.06%). Thus, in the absence of the IN TP, MLV integration preference for heterochromatin and weakly transcribed regions increases, while integration at active promoters and enhancers was diminished.
Decreased integration at BET-associated epigenetic marks for IN TP- MLV in K562 cells
Chromosomal states are defined in part through profiling combinations of known histone modifications and the occupancy of various cis-regulatory elements by known protein factors [55]. Directed by the bromodomains, acetylated histone modifications are important in determining BET proteins interactions, however additional marks including H3K4me2/3 are elevated in Brd-bound nucleosomes [49]. MLV integrations strongly associate with H3K4me1, H3K4me2, H3K4me3, H3K27ac, H2Az and H3K9ac modified chromatin [5, 12, 13, 15, 44, 51]. Fig 6D and S4 Fig show the proportion of WT and IN TP- integrants found in regions marked by the various epigenetic marks and proteins associated with integration site preference in K562 cells. ChIP-seq data for all chromatin modifications was obtained from the ENCODE consortium, and enrichment of viral integration over RIC was computed for each modification. Significantly, the largest decrease in enrichment upon loss of the IN TP was associated with Brd4 binding sites (17.5-fold for WT vs 7.1-fold for MLV IN TP-) (Fig 6D and S4 Fig). The WT virus was most highly enriched at the epigenetic modifications H3K4me3, H3K4me2, H3K27ac and H3K9ac, which is consistent with previous reports [13]. Interestingly, loss of IN TP resulted in a >45% decrease in fold enrichment at H3K4me3, H3K4me2, and H3K27ac sites, which is consistent with the loss of association with BET proteins.
Discussion
MYC/Runx2 transgenic mouse model and MLV-induced tumorigenesis
The discovery that the MLV IN protein interacts with host BET proteins, directing integration into sites of active expression marked by acetylated chromatin, gave rise to the question of whether the loss of the MLV BET binding site would affect the pathogenesis of the virus. In this study, we used MYC/Runx2 mice, a transgenic model displaying rapid tumor formation that is further accelerated by MLV infection [18, 19, 21]. These mice were infected with MLV lacking the IN tail peptide. Globally, the infection time course for the IN TP- series showed a biphasic DoD curve, with one-third of the mice developing tumors early, as with WT infection, and two-thirds of the mice developing tumors later, paralleling the uninfected MYC/Runx2 control mice. We extensively characterized IN TP-16 samples, which displayed early tumorigenesis and showed no detectable signs of recombination with endogenous viruses. Characterization of TP-16 tumor viral integration sites indicated a profile dispersed away from TSSs and decreased association with H3K27ac marks (Fig 4A and 4C). The IN mutation was preserved within the tumor, and the characteristic target-site profile for IN lacking the TP was maintained following 293mCAT cell infection with IN TP-16 tumor-derived virus. However, integrations into CISs known to accelerate tumor formation were also present, suggesting that stochastic integration into an open chromatin hotspot can still provide the positive selection required for tumor outgrowth. It is noteworthy that 7/10 of the highest MLV TP-16 copy number integration sites were within known CISs (Table in S5 Table).
Selective pressure to maintain the IN tail peptide in mice
Our results highlight the strong pressure within the in vivo model to maintain the presence of the IN-TP. Interestingly, the TP region overlaps with the Env signal peptide, yet the codon bias within 4070A MLV is maintained towards the IN reading frame [56]. Sequence conservation between gammaretroviruses identified the conserved sequence W(X7)PLK(I/L)R within the TP [4]. From our studies, the initial viruses (MLV IN-XN) emerging from the tumors indicated that selective pressures restored the TP through deletion mutations within the virus coding region, removing a small putative linker region within the IN C-terminus to restore the TP in an alternative reading frame. When this option was removed, through optimizing codon usage for the IN/Env region towards the env sequence, thereby destroying the coding potential of the TP, circulating virus restored the TP through recombination with endogenous polytropic viruses, with cross-over junctions within the IN and Env coding regions.
Identification of a linker region between the IN CTD SH3 and TP
The deletion mutations that restored the BET interacting domain onto the IN-XN construct identify a linker region between the IN CTD SH3 fold (ending at KAADPG [3, 40, 57] and the IN BET interaction domain (initiating at IN W390 [4, 15]). After passage in mice, the spacing between the IN SH3 and TP domains is reduced from nine amino acids (GGGPSSRLT) to 3 amino acids (GRK), two of which are not normally encoded by the IN protein. This indicates that the length of this linker region, as well as the composition of the amino acids between the IN SH3 fold and the BET binding site can be substantially altered while still maintaining virus viability in vivo [57].
Recombination with ERVs and temporal effects of tumorigenesis
Recombination with the endogenous polytropic Pmv20 virus was detected in multiple tumors (IN TP-4, 6, 7, 9, 15, 18, and 19). While this recombination restored the IN TP, the presence of recombinant virus did not correlate with mouse DoDs; IN TP-16 did not recombine with polytropic virus (DoD 34 days), whereas IN TP-4, 6, 7, 9, 15, 18, and 19, which did undergo recombination, had DoDs of 50, 57, 37, 63, 33 and 36 days, respectively. Our analysis cannot determine the point in the infection time course at which recombination occurred, and it is possible that the recombination event for IN TP-9, for example, which had a DoD of 63 days, occurred late in tumor development and thus did not have a marked effect on tumor progression. It is interesting to note that following infection with IN TP-, no mice survived beyond day 70, whereas some NC mice survived beyond 110 days, suggesting that the presence of a recombinant virus may have affected long-term survival.
Sequence analysis of viral/host DNA junctions
Interestingly, when the 3’ LTR junctions of 3 of the top 10 RISs (Mapk13, Ccnd1, and Hdac6) were sequenced as a PCR population, heterogeneity was observed at positions +1/+2 downstream of the MLV LTR TCTTTCA terminus. Viral integration involves the cleavage of the LTR terminal dinucleotides, exposing the conserved 3’ CA for strand transfer and generating a 5’ ss AA tail at the 5’ ends of the viral DNA substrate. The observed sequence heterogeneity corresponded with either the predicted host DNA sequence or nucleotides from repair of the 5’ ss viral DNA tail. Heterogeneity at the 5’ LTR-host junction was not observed within the Illumina sequence reads at these genes nor in the 5’LTR analysis of the Hdac6 gene. The mechanism for this heterogeneity requires further analysis and could reflect simple sequencing compression or a unique repair mechanisms at the viral/host DNA junction in the MYC/Runx2 mice model, in which p53 activity is reported to be suppressed [18, 19].
Mechanisms of retroviral enhancer activation
For retroviruses, enhancer activation usually occurs upstream of the gene in the antisense orientation, or downstream in the sense orientation [54]. Indeed, this was the orientation bias observed for integrations into both Ccnd1 and Rasgrp1 in TP-16 tumor DNA. For Ccnd1, the top integrant was in the integration cluster near the promoter (Fig 5A). Similarly, previous studies of WT MLV in MYC/Runx2 tumors indicated a cluster of insertions at the 5’ end of the Ccnd1 gene, predominantly upstream of the coding sequence [20]. The proximal cluster initiates overexpression via retroviral enhancer elements [58]. Ccnd1 has an important role in cell cycle regulation, and overexpression induces the formation of different cancer types [59–61]. For Rasgrp1, the integrant analyzed was 91,758 bases upstream of the Rasgrp1 promoter, and oriented with the viral promoter in the opposite orientation. Activation of Rasgrp1 would most likely occur through an enhancer activation event and sequence analysis indicated this integrant maintained the IN TP- sequence. It was initially surprising to see the RIS orientation bias >91 kb upstream of Rasgrp1 promoter. For human T cell leukemia virus-1 (HTLV-1), long-range interactions between target gene promoters and viral enhancers are facilitated through chromatin looping utilizing the host zinc finger binding protein CTCF [62, 63]. The HTLV-1 provirus contains a CTCF nucleotide-binding motif that has been shown to mediate clone-specific deregulation of host transcription from distances up to 300 kb [62], and CTCF-mediated cis contacts within the host genome can be as far as 1.4 Mb [62, 63]. Although M-MLV does not encode a known CTCF binding motif, CTCF binding sites have been identified at the promoter region of Rasgrp1 and ~10 kb downstream from the integration site (ENCODE reference ENCFF310MUQ). CTCF-mediated transcription varies depending on the cell type. Validation of CTCF binding would require circular chromosome conformation capture (3C or ChIA-PET) analysis from the tumors, which is not available for this study. However, the presence of these CTCF binding sites provides a potential mechanism for the MLV enhancer to interact with Rasgrp1 promoters that are distant from each other, thereby driving overexpression of Rasgrp1 concomitant with tumorigenesis in these mice [64–66]. For Ccnd1, CTCF-dependent long-range loops have been identified that reposition distal clusters of retroviral insertions, driving gene activation [58].
Biased integration orientation within the HDAC6 gene
The Hdac6 gene also displayed biased integration in an orientation opposite to that of transcription, however these integrants were within the Hdac6 gene. For IN TP-16, the three most abundant integrants in the library mapped within exon 3 proximal to the gene 5’ end and within introns at the center and toward the 3’ end of the gene. The results imply a loss of function through oncogenic selection, however the mechanism cannot be determined. Hdac6 is reported to interact with Runx2 [67, 68] as well as being involved in multiple cellular processes, including organization of the immune synapse, cell migration, protein degradation, and viral infections [69].
IN domains involved in BET protein recognition
The absence of the IN TP reduced the integration bias for strong enhancers and active promoters and increased preferences for heterochromatin and weakly transcribed regions. The secondary preference towards active regulatory elements (state 1, 4 and 5) was maintained, which corresponds to median genome coverage of only 2.5%. The observed integration of MLV TP- at active promoter/enhancers could be the result of additional IN sequences interacting with BET proteins [42]. Although our studies indicated that the interaction with BET domains occurs predominantly through the IN TP, others have mapped three amino acid changes within the IN CCD, specifically E266A, L268A, and Y269A, which resulted in the loss of Brd2 interaction as assessed by coimmunoprecipitation [42]. We have analyzed these mutants in the context of single-round infection, which indicated that all three mutant viruses had negligible titers, akin to the IN catalytic mutant virus D184N (Fig 1D). Molecular modeling indicates that these positions are located in the predicted MLV CCD α6 helix (residues 264–270) and are close to the CCD dimer interface (S5 Fig). The model shows the close proximity of residues Y269 and L268 between monomers and thus alanine substitutions of L268 and Y269 could disrupt the multimerization of the CCD. Thus the effects of L268A and Y269A changes on IN binding to Brd2 were plausibly indirect.
Chromatin marks associated with MLV integration and transcriptionally quiescent regions
Heterochromatin and weakly transcribed regions in the 15-state model share the absence of chromatin marks H3K4me1/2/3, H3K27me3, H3K27ac, H3K9ac and CTCF, with the weakly transcribed state 11 containing low levels of H3K36me3 and H4K20me1 [55]. A more defined model with 25 and 50 chromatin states has recently been described [70], which makes further subdivisions based on an expanded set of chromatin marks. In these, the heterochromatin state is distinguished from the quiescent state by the presence of the H3K9me3 mark. While K562 ENCODE data is not available for the full expanded set of chromatin marks, we observed no significant fold enrichment of WT MLV or IN TP- integrations at H3K9me3. Therefore, most of the MLV IN TP- integrations in heterochromatin would reasonably be categorized as targeting the quiescent chromatin state. Similar to the heterochromatin state in 15-chromatin scale, the quiescent state is defined by the large absence of any histone modifications and is transcriptionally inactive with low annotated non-coding and coding transcripts [70, 71].
Recognition nucleosomes by MLV IN
This is the first study to define where MLV integration is directed in the absence of the BET protein interaction, beyond noting a decreased bias away from active promoters and enhancers [12, 13, 72]. Two models could explain the integration preference into regions with limited modified histones observed for State 11. In the first model, MLV IN may display an innate recognition of unmodified histone tails. In the second model, modified histones may present a steric hindrance for IN binding. Both MLV and prototype foamy virus (PFV) encode an N-terminal extension domain (NED) [2, 8] though PFV IN does not encode a homologous TP. Two loop regions within the PFV IN CCD-CCD dimer interface interact with the H2A-H2B heterodimer, specifically with the C-terminal helix of H2B and N-terminal of H2A [73–75]. Binding to the nucleosome results in a 7Å deformation of the target DNA and ultimately drives viral integration into heterochromatin regions, Lamin A/B1 rich-regions, and intergenic regions [73–76]. Transposing the loop regions of the PFV IN CCD that are implicated in the H2B interaction onto MLV IN resulted in viral titers equivalent to those of the catalytically inactive MLV D184N implying that these surface exposed loops in MLV IN are integral to IN stability or sites of secondary interactions. H2A and H2B are the most diverse histones, which along with distinct post-translational patterns, contribute to the complexity and variability of H2A-H2B dimers [77]. It is possible that the MLV IN CCD-CCD dimer may interact with a specific variant of H2A-H2B heterodimer with distinct modifications independent of the BET proteins.
Potential of phenotypic complementation of IN by ERVs
In the mouse model, we cannot exclude a role for phenotypic complementation of IN function through co-packaging of TP- Gag-Pol with full-length Gag-Pol expressed by ERVs in tumor cells, specifically transcomplementation by endogenous Gag-Pol proteins from either the Emv2 [22, 78] or the xenotropic MLV Bxv-1 loci. C57BL/10 mice have been documented to express xenotropic MLV at low levels in vivo, which can be induced in tissue culture with IdU [29]. Thus, although the xenotropic Env would exclude infection by these endogenous viruses, complementation of the IN protein in trans cannot be eliminated. This could result in the low level of bias towards the TSS integration observed within the MYC/Runx2 tumors in the absence of IN TP region recombination events. This issue of transcomplementation by xenotropic or ecotropic MLV would not be expected to be problematic were IN TP- virus to be used in non-murine cells, including human K562 cells. However, each species has their own ERVs, which are much less likely to contribute to the target-site selection of MLV-based vectors.
MLV IN TP- vectors for gene therapy
The overall goal of these experiments was to determine whether MLV IN TP- might represent a safer gene delivery vector as compared to MLV harboring WT IN. Importantly, the uniform acceleration of tumorigenesis by MLV in the MYC/Runx2 mouse was not observed with IN lacking the BET interaction domain at its C-terminus. However, a major limitation of this mouse system is the strong selective pressure to maintain IN-TP function, either directly by recombination or indirectly via transcomplementation. Alternative animal models or assay systems [79, 80] would be beneficial to assess the full potential of MLV IN TP- as a vector in the absence of endogenous elements influencing integration preferences. Our results indicate that removing the TP was insufficient to redirect all integrations away from active promoters and strong enhancers, or to eliminate the stochastic events that can select for oncogenic activation. Ultimately, a modified vector that combines SIN LTRs to eliminate strong viral enhancers [10, 81] with insertions/replacement of the IN TP- to redirect integration to less active regions [39, 82] could decrease vector genotoxicity and overcome current limitations for clinical applications.
Materials and methods
Ethics statement
All animal experiments were performed under the EU directive 2010 and UK Animals (Scientific Procedures) Act 1986. This specific study was undertaken on project license number P3C97B34C (awarded to E Cameron) ethically approved by the University of Glasgow (Establishment license number XC2FD842E) Animal Welfare & Ethics Review Board and the UK Home Office.
Cell lines
Human embryonic kidney 293T cells, 293mCAT cells (expressing mCAT receptor) [83] and canine osteosarcoma D17/pJET cells (expresses the mCAT-1 receptor) [43] were maintained in Dulbecco's modified Eagle medium (DMEM; Gibco #11965) supplemented with 10% (vol/vol) heat-inactivated fetal bovine serum (Atlanta Biologicals # S1245OH) and 1x antibiotic-antimycotic (100 units/mL of penicillin, 100 μg/mL of streptomycin, and 0.25 μg/mL Amphotericin B) (Gibco #15240). The parental 293, 293T and D17 cells were obtained from ATCC. Human K562 cells were acquired from ATCC (CCL-243) and maintained in Iscove’s modified Dulbecco’s medium (IMDM, # 12440079) supplemented with 10% (vol/vol) heat-inactivated fetal bovine serum (Atlanta Biologicals # S1245OH).
Plasmid and vector construction
The replication-competent M-MLV proviral construct pNCA-C [84] and pNCA-C IN-XN (previously named in6215a [40]), bearing a 23-aa truncation of the IN tail peptide (TP) of the C-terminal domain (CTD) was previously described [40]. To generate a codon-optimized pNCA-C-TP-, a 137 bp gene block (IDT) was chemically synthesized and amplified using primers NCACXN_ScaI6330_rev and NCACXN_NotI6220_fwd. Overlapping PCR of this fragment with a ScaI-ClaI fragment from pNCA-C (generated using primers NCAC_8290_rev and NCACXN_6327_fwd) resulted in a NotI-ClaI fragment, which was exchanged into NotI/ClaI digested pNCA-C IN-XN. Generation of the pNCA-C IN D184N was previously described [41]. Sequences of all oligonucleotides are provided in Table in S3 Table.
The three residues implicated for BET protein binding (MLV IN E266, L268, Y269) were substituted to alanine using overlapping PCR with KOD polymerase [42]. PCR of the first fragment was amplified with primer 102510NdeIINteinIN forward and point mutant specific reverse primer (E266A_rev, L268A_rev, Y269A_rev) and the second fragment was amplified using the point mutant specific forward primer (E266A_fwd, L268A_fwd, Y269A_fwd) and the 102510XhoIInteinIN1-407 reverse. The overlapping PCR fragment was introduced into pNCA-C using the HindIII and PmlI sites. All mutations were verified using Sanger DNA sequencing.
DEAE-dextran transient transfection of proviral DNA clones in D17pJET cells
Transient expression of the pNCA-C based proviral constructs was performed as previously described [3, 85] using 500 ng pNCA-C based plasmids. Tissue culture supernatant was monitored for viral spread using enzyme-linked immunosorbent assay (ELISA) against MLV p30 [86]. Cultures were maintained for at least 14 days prior to analysis.
LacZ viral titer assay and single round infection in D17 cells
2 x 106 293T cells were transfected with 0.8 μg pMD2.G (Addgene) expressing the vesicular stomatitis virus glycoprotein (VSV-G), 0.8 μg pRT43.2Tnlsβ-gal [87] a retroviral packaging vector expressing lacZ and 0.8 μg M-MLV viral genome with wild-type (pNCA-C) or MLV TP- (pNCA-C TP-) using Fugene 6 (Promega #E2691) overnight as directed by the manufacturer [88]. Viral supernatant was collected, filtered through a 0.45 μm syringe and viral particles were quantified using ELISA against the MLV p30 [86]. 1x105 cells D17 cells on 3.5-cm gridded plates were infected with media containing 10 ng of CA in 2 ml DMEM in the presence of 8 μg/ml polybrene. Medium was replaced with fresh DMEM after 24 h of infection. Cells were stained for LacZ expression as previously described [45].
Western blot
Viruses were collected from D17/pJET viral producer cell lines. For CA and Env western blot, 2 ml of viral supernatant was spun at 15,000 x g for 30 min and the viral pellet was resuspended in 20 μl phosphate buffer saline (PBS). Samples were run on a 10% SDS-PAG and transferred to polyvinyl difluoride (PVDF) membranes using Bio-Rad Trans-Blot Turbo Transfer System. Immunoblots were developed using goat anti-p30 (CA) (1:2000, 81S-263) and goat anti-Env (1:1000, 80S-019) (Quality Biotech) with bovine anti-goat HRP (80S-035-180) as secondary antibody (1:10,000). For IN western blot, 10 ml of viral supernatant was pelleted at 15,000 x g for 30 min and the proteins were visualized using 1:1 mix (1:1000) of antiserum from Rabbit 3 and Rabbit 4, Bleed 5 with goat anti rabbit HRP (Pierce #31460) as secondary antibody (1: 5,000) [89].
Infection of MLV into MYC/Runx2 mice
The MYC/Runx2 transgenic mice are on a C57BL6 x CBA/Ca background. Infection and maintenance of the mice was as described [21]. For WT and mutant MLV, viruses were obtained from 293mCAT cells to avoid recombination with endogenous viruses prior to infection. Briefly, mice were inoculated intraperitoneally with virus isolated from tissue culture supernatant (105 TCID50) within 24 h of birth. Date of Death (DoD) was monitored over a 115-day period. For each mouse, a tissue fragment was thawed from liquid N2, minced and incubated in medium at 37°C for 2–3 hrs. The medium was spun at 1,200 rpm and then filtered (0.45μm) before adding to 293mCAT cells. Cells were cultured for at least 7 days before harvest for DNA isolation. Amplification of the integrated MLV genomes from the 293mCAT cells was performed using primers 4924 and 7791, previously named 3807 and 6320, respectively [45] (Table in S3 Table). PCR products were cloned using TA cloning. Individual colonies from mice XN-2, 3, and 35 were selected and sequenced for presence of the TP coding region.
Detection of mouse DNA
Genomic DNA isolates from 293mCAT cells were analyzed for mouse DNA contamination by examining for mouse intracisternal particle A (IAP) and mouse mitochondrial cyclooxygenase-2 (mCOX2). The primers used were mouse_IAP_fwd (ATAATCTGCGCATGAGCCAAGG) and Mouse_IAP_rev (AGGAAGAACACCACAGACCAG) [90], and primers used for COX2 were Mouse_mt__COX2_fwd (5′ TTC TAC CAG CTG TAA TCC TTA 3′) and Mouse_mt_COX2__rev (5′ GTT TTA GGT CGT TTG TTG GGA T 3′) [91]. PCR analysis was performed using KOD HotStart.
Analysis of M-MLV recombination with endogenous retroviruses
Detection of recombination with ERVs utilized primers recognizing M-MLV, polytropic, xenotropic ERV and the RT_universal_primer (Table in S3 Table) [92]. PCR analysis on genomic DNA of 293mCAT virus infected cells was performed using a combination of RT universal primer (5’ CCTACTCCGAAGACCCCTCGA-3’) and primers specific for polytropic and xenotropic ERVs (Table in S2 Table) using KOD Hotstart polymerase (Millipore, 71086) according to suggested parameters. PCR products from the reaction with RT_Universal_primer and Polytropic_JS5_rev on 293mCAT infected cells from TP-6, 7 and 9 mice were cloned into pCR4-TOPO vector using the TOPO TA kit following the protocol provided by the manufacturer (Invitrogen, K4575-40). Recombinant plasmids were sequenced using the T3/T7 sequencing primers from the manufacturer, and 4981_fwd and MLV_IN_T159A_fwd to determine the 5’ recombination junction. PCR with Polytropic_JS5_fwd and 7791_reverse primers determined the 3’ recombination junction for the same 293mCAT samples. PCR analysis of genomic DNA from mice tumor or thymus samples required a nested PCR. First round PCR used the RT_universal_fwd primer and the MLV_LTR_U3_rev primer. Second round of PCR used a primer pair of 7791_reverse primer and the Polytropic_JS5_ fwd primer.
Single round of infection of K562 cells
Transfection of 293Lenti-X cells using Mirus TransIT-Lenti transfection reagent with WT MLV and MLV TP−along with pMD2.G (Addgene) generated viruses for single round infection [93]. Viruses were quantitated using ELISA as previously described [86]. K562 cells (5x105 cells) were prepared a day prior for infection in 6-well plates. For WT MLV, 500 ng of p30 was added to one well of K562 cells and for MLV IN TP-, 5000 ng was added. The plate was spinoculated at 1,500 g for 1 h followed by incubated for 4 h at 37°C [94]. Supernatant was removed and the cells were grown for 24 h. The cells were collected for genomic DNA extraction (see below).
Integration target-site analysis
Genomic DNA from infected mouse tumor, thymus, 293mCAT, and K562 cells was extracted (QIAGEN #6941) and used to generate libraries for MLV integration sites. Control uninfected libraries were generated from C57BL/6J genomic DNA. Protocols for library preparation were adapted from [93]. Genomic DNA sample (5 μg) were subjected to two rounds of sonication with the following parameters: duty cycle: 5%; intensity: 3; cycles per burst: 200; time: 80 sec). Purification of DNA for the next generation library protocol used MinElute Reaction Cleanup kit (Qiagen 28204). The sonicated DNA was purified and ends of DNA fragments were repaired using End-It DNA End-Repair Kit (ER0720) as described in manufacturer’s protocol and purified after repaired DNA ends were A-tailed using Klenow Fragment (M0212S). All kits were used as described by the manufacturer’s protocol. Linker short and long strands (Table in S6 Table) were annealed by heating to 90°C and slowly cooling to room temperature in steps of 1°C per min. The annealed linkers were ligated with assigned genomic DNA sample with 3000 U of T4 ligase (M0202M) overnight at 12°C and purified. The first round of PCR used a LTR specific primer and linker specific primer (Table in S6 Table) with adapter sequence and primer binding sequence adapted from [93]. For the mouse DNA samples, the MLV_LTR_U3_rev primer (5’- GCGTTACTTAAGCTAGCTTGCCAAACCTAC-3) was used [54]. For 293mCAT and K562 cells, MLV_LTR_U5 primer (5’-CCTTGGGAGGGTCTCCTCTGAGT-3’) was used. Four PCR reactions of 100 ng DNA each were setup for each genomic DNA sample with PCR KOD Hotstart polymerase (Millipore, 71086) under these parameters: One cycle: 98°C for 2 min; 30 cycles: 98°C for 15 sec, 60°C for 30 sec, 70°C for 45 sec. The reactions were pooled and purified. The second round of PCR used a second round LTR specific primer and the same linker specific primers. These second round LTR specific primers encode a 6-nucleotide index or barcode sequence compatible for NGS, an adapter sequence for DNA clustering and a sequencing primer binding site (Table in S7 Table). Reaction parameters for the second round of PCR paralleled those of the first PCR round. All reactions were pooled and purified. Libraries were analyzed for quality and sequenced using the Illumina MiSeq system at the Molecular Biology Core Facilities at the Dana-Farber Cancer Institute. The sequences reported in this paper have been deposited in the National Center for Biotechnology Information Sequence Read Archive (Bioproject id # PRJNA548288).
Bioinformatics
Bioinformatic analyses of integration sites were performed as described in [95]. LTR and linker sequences were cropped from 150 bp paired end reads using custom Python scripts, and the cropped reads were mapped to the reference genome (mm10 for mouse samples and hg19 for libraries from human cell lines) using HISAT2 [96]. Results were then filtered to retain high-quality alignments using SAMtools [97] and unique (deduplicate) integration sites were extracted and formatted to the browser extensible data (BED) format using custom Python scripts. Copy number of MLV integration sites from the tumor samples was calculated post filtering of high-quality alignments and prior to deduplication using custom R scripts. Copy number was calculated as described [98]. Briefly, copy number of the integrants was defined as the number of sequences having the same integration site but different breakpoints in the host DNA. Only breakpoints that were >3nt apart were counted as independent events.
Integration sites obtained from the tumor of non-infected control (NC) mouse was considered as background amplification of endogenous retroviruses obtained using this pipeline, and hence integration sites from all tumor samples overlapping with the sites from NC were computationally removed. BEDtools software suite [99] was then used to correlate unique integration sites proximal to genomic annotations such as TSS and CpG islands obtained from the University of California Santa Cruz (UCSC) database (http://genome.ucsc.edu/cgi-bin/hgTables)). Fraction of integration sites enriched at chromatin associated with Brd4 binding sites and various histone modifications (Table in S8 Table) was also computed using BEDtools suite. For ChipSeq datasets based on prior genomic versions, the coordinates were converted to the genome build used in the study (mm10 for mouse samples and hg19 for human samples) using the liftover utility from the UCSC database. Distance of integration sites from TSS, CpG islands, and Brd4 binding sites were calculated using BEDtools and histograms comparing the obtained distribution of integrants were plotted using ggplot2 [100].
Genomic annotations showing the chromatin state segmentation of K562 (wgEncodeEH000790) defined by HMM from ENCODE/BROAD was downloaded from the UCSC genome browser. Custom R scripts were written to segregate the individual chromosome state definitions from the master file and BEDtools was used to correlate integration sites within each chromosome states.
Statistical analyses
Single factor ANOVA test was performed as described in [101] to confirm the significant changes within the experiment (P <0.05) and Fisher’s exact test was used from computing the statistical comparisons (Table in S1 Table). Other utilized tests were as described in text.
Amplification of MLV integrants at targeted loci
TP-16 integration-specific primers were designed based on the genomic location of integrants mapped using NGS (Table in S3 Table). First round PCR used forward primer at the mouse genomic sequence upstream to the integrant and MLV_LTR_U3 reverse primer using PrimeStar GXL DNA polymerase (Takara R050A) or RT_universal_fwd primer and the TP-16 specific reverse primer at the mouse genomic sequence downstream from the integrant using KOD HotStart polymerase. Second round PCR used the same TP-16 specific mouse genomic primer and either primers MLV 20R_reverse or NCAC 6327_reverse or NCAC 5166_forward. For some integrants (Hdac6 intron 28 and near Ccnd1), additional single linear amplification PCR to amplify the first round PCR was included [102]. Products were sequenced using either LTR_outside_fwd, E266K_fwd, or gene specific primers.
Analysis of retroviral insertions overlapping H3K27ac peaks
H3K27ac ChIP-seq narrow peak files from C57BL/6 mouse thymus were obtained from ENCODE (project code ENCFF001KYC). Overlap between the H3K27ac peaks and retroviral insertion sites for each sample was assessed using the R Bioconductor packages ChIPpeakAnno and GenomicRanges. Distances between H3K27ac peaks and retroviral insertion sites were mapped using the distanceToNearest() function from the R Bioconductor package GRanges, then density plots were produced using the ggpubr and ggplot2 packages. The statistical significance of differences between distances was assessed using Wilcoxon Rank Sum tests.
Supporting information
S1 Fig [tif]
Western blot analysis of MLV -associated proteins.
S2 Fig [erv]
Scheme of the MLV genome and primers.
S3 Fig [a]
PCR for MLV recombinants from mice tumors and detection of mouse DNA contamination in viral infected human cell line 293mCAT.
S4 Fig [tif]
Percentage and fold enrichment WT MLV and MLV TP within different histone modifications and Brd4 identified by ChipSeq data for K562 cells.
S5 Fig [a]
MLV IN CCD and prototype foamy virus (PFV) IN CCD.
S1 Table [docx]
Fisher’s test for statistical comparison of integration profile.
S2 Table [docx]
mice infected with WT MLV and MLV IN TP viruses.
S3 Table [docx]
List of oligonucleotide primers.
S4 Table [docx]
Percent overlap of RISs from mouse tumors and H3K27Ac peaks.
S5 Table [docx]
Comparison of known CISs with MLV IN TP16 tumor integrants.
S6 Table [docx]
Linker sequence and linker specific primers.
S7 Table [docx]
MLV LTR specific primers for second round PCR.
S8 Table [docx]
Genomic annotations and ChipSeq datasets used in the study.
Zdroje
1. Fan H, Johnson C. Insertional oncogenesis by non-acute retroviruses: implications for gene therapy. Viruses. 2011;3(4):398–422. doi: 10.3390/v3040398 21994739
2. Guan R, Aiyer S, Cote ML, Xiao R, Jiang M, Acton TB, et al. X-ray crystal structure of the N-terminal region of Moloney murine leukemia virus integrase and its implications for viral DNA recognition. Proteins: Structure, Function, and Bioinformatics. 2017;85(4):647–56.
3. Aiyer S, Rossi P, Malani N, Schneider WM, Chandar A, Bushman FD, et al. Structural and sequencing analysis of local target DNA recognition by MLV integrase. Nucleic Acids Res. 2015;43(11):5647–63. doi: 10.1093/nar/gkv410 25969444
4. Aiyer S, Swapna GV, Malani N, Aramini JM, Schneider WM, Plumb MR, et al. Altering murine leukemia virus integration through disruption of the integrase and BET protein family interaction. Nucleic Acids Res. 2014;42(9):5917–28. doi: 10.1093/nar/gku175 24623816
5. De Rijck J, de Kogel C, Demeulemeester J, Vets S, El Ashkar S, Malani N, et al. The BET family of proteins targets Moloney murine leukemia virus integration near transcription start sites. Cell Rep. 2013;5(4):886–94. doi: 10.1016/j.celrep.2013.09.040 24183673
6. Roth MJ, Schwartzberg PL, Goff SP. Structure of the termini of DNA intermediates in the integration of retroviral DNA: dependence on IN function and terminal DNA sequence. Cell. 1989;58(1):47–54. doi: 10.1016/0092-8674(89)90401-7 2546673
7. Craigie R, Fujiwara T, Bushman F. The IN protein of Moloney murine leukemia virus processes the viral DNA ends and accomplishes their integration in vitro. Cell. 1990;62(4):829–37. doi: 10.1016/0092-8674(90)90126-y 2167180
8. Engelman A, Cherepanov P. Retroviral integrase structure and DNA recombination mechanism. Microbiology Spectrum. 2014;2(6):1011–33.
9. Kvaratskhelia M, Sharma A, Larue RC, Serrao E, Engelman A. Molecular mechanisms of retroviral integration site selection. Nucleic Acids Res. 2014;42(16):10209–25. doi: 10.1093/nar/gku769 25147212
10. Maetzig T, Galla M, Baum C, Schambach A. Gammaretroviral vectors: biology, technology and application. Viruses. 2011;3(6):677–713. doi: 10.3390/v3060677 21994751
11. Serrao E, Ballandras-Colas A, Cherepanov P, Maertens GN, Engelman AN. Key determinants of target DNA recognition by retroviral intasomes. Retrovirology. 2015;12:39. doi: 10.1186/s12977-015-0167-3 25924943
12. De Ravin SS, Su L, Theobald N, Choi U, Macpherson JL, Poidinger M, et al. Enhancers are major targets for murine leukemia virus vector integration. J Virol. 2014;88(8):4504–13. doi: 10.1128/JVI.00011-14 24501411
13. LaFave MC, Varshney GK, Gildea DE, Wolfsberg TG, Baxevanis AD, Burgess SM. MLV integration site selection is driven by strong enhancers and active promoters. Nucleic Acids Res. 2014;42(7):4257–69. doi: 10.1093/nar/gkt1399 24464997
14. Hsu SC, Blobel GA. The role of Bromodomain and Extraterminal Motif (BET) proteins in chromatin structure. Cold Spring Harb Symp Quant Biol. 2017;82:37–43. doi: 10.1101/sqb.2017.82.033829 29196562
15. El Ashkar S, De Rijck J, Demeulemeester J, Vets S, Madlala P, Cermakova K, et al. BET-independent MLV-based vectors target away from promoters and regulatory elements. Mol Ther Nucleic Acids. 2014;3:e179. doi: 10.1038/mtna.2014.33 25072693
16. Steffen D. Proviruses are adjacent to c-myc in some murine leukemia virus-induced lymphomas. Proc Natl Acad Sci U S A. 1984;81(7):2097–101. doi: 10.1073/pnas.81.7.2097 6326104
17. Hay J, Gilroy K, Huser C, Kilbey A, Mcdonald A, MacCallum A, et al. Collaboration of MYC and RUNX2 in lymphoma simulates T-cell receptor signaling and attenuates p53 pathway activity. J Cell Biochem 2019, 120(10):18332–18345. doi: 10.1002/jcb.29143 31257681
18. Neil JC, Gilroy K, Borland G, Hay J, Terry A, Kilbey A. The RUNX Genes as conditional oncogenes: Insights from retroviral targeting and mouse models. In: Groner Y, Ito Y, Liu P, Neil JC, Speck NA, van Wijnen A, editors. RUNX Proteins in Development and Cancer. Singapore: Springer Singapore; 2017. p. 247–64.
19. Blyth K, Vaillant F, Hanlon L, Mackay N, Bell M, Jenkins A, et al. Runx2 and MYC collaborate in lymphoma development by suppressing apoptotic and growth arrest pathways in vivo. Cancer Res. 2006;66(4):2195–201. doi: 10.1158/0008-5472.CAN-05-3558 16489021
20. Huser CA, Gilroy KL, de Ridder J, Kilbey A, Borland G, Mackay N, et al. Insertional mutagenesis and deep profiling reveals gene hierarchies and a Myc/p53-dependent bottleneck in lymphomagenesis. PLoS Genet. 2014;10(2):e1004167. doi: 10.1371/journal.pgen.1004167 24586197
21. Stewart M, Mackay N, Hanlon L, Blyth K, Scobie L, Cameron E, et al. Insertional mutagenesis reveals progression genes and checkpoints in MYC/Runx2 lymphomas. Cancer Res. 2007;67(11):5126–33. doi: 10.1158/0008-5472.CAN-07-0433 17545590
22. Kozak CA. Origins of the endogenous and infectious laboratory mouse gammaretroviruses. Viruses. 2014;7(1):1–26. doi: 10.3390/v7010001 25549291
23. Greenwood AD, Ishida Y, O'Brien SP, Roca AL, Eiden MV. Transmission, evolution, and endogenization: Lessons learned from recent retroviral invasions. Microbiol Mol Biol Rev. 2018;82(1):e00044–17. doi: 10.1128/MMBR.00044-17 29237726
24. Stoye JP, Coffin JM. The four classes of endogenous murine leukemia virus: structural relationships and potential for recombination. J Virol. 1987;61(9):2659–69. 3039159
25. Giovannini D, Touhami J, Charnet P, Sitbon M, Battini JL. Inorganic phosphate export by the retrovirus receptor XPR1 in metazoans. Cell Rep. 2013;3(6):1866–73. doi: 10.1016/j.celrep.2013.05.035 23791524
26. Fan H. Leukemogenesis by Moloney murine leukemia virus: a multistep process. Trends Microbiol. 1997;5(2):74–82. doi: 10.1016/S0966-842X(96)10076-7 9108934
27. Stoye JP, Moroni C, Coffin JM. Virological events leading to spontaneous AKR thymomas. J Virol. 1991;65(3):1273–85. 1847454
28. Ott D, Friedrich R, Rein A. Sequence analysis of Amphotropic and lOAl murine leukemia viruses: Close relationship to Mink Cell Focus-Inducing Viruses. J Virol. 1990;64(2):757–66. 2153240
29. Kozak C, Rowe WP. Genetic mapping of xenotropic leukemia virus-inducing loci in two mouse strains. Science. 1978;199(4336):1448. doi: 10.1126/science.204014 204014
30. Kohn DB. Historical perspective on the current renaissance for hematopoietic stem cell gene therapy. Hematol Oncol Clin North Am. 2017;31(5):721–35. doi: 10.1016/j.hoc.2017.06.006 28895843
31. Milone MC, Bhoj VG. The pharmacology of T cell therapies. Mol Ther—Methods Clin Dev. 2018;8:210–21. doi: 10.1016/j.omtm.2018.01.010 29552577
32. Hacein-Bey-Abina S Von Kalle C, Schmidt M, McCormack MP, Wulffraat N, Leboulch P, et al. LMO2-associated clonal T cell proliferation in two patients after gene therapy for SCID-X1. Science. 2003;302((5644)):415–9. doi: 10.1126/science.1088547 14564000
33. Hacein-Bey-Abina S, Garrigue A, Wang GP, Soulier J, Lim A, Morillon E, et al. Insertional oncogenesis in 4 patients after retrovirus-mediated gene therapy of SCID-X1. J Clin Invest. 2008;118(9):3132–42. doi: 10.1172/JCI35700 18688285
34. Hacein-Bey-Abina S, Hauer J, Lim A, Picard C, Wang GP, Berry CC, et al. Efficacy of gene therapy for X-linked severe combined immunodeficiency. N Engl J Med. 2010;363(4):355–64. doi: 10.1056/NEJMoa1000164 20660403
35. Stein S, Ott MG, Schultze-Strasser S, Jauch A, Burwinkel B, Kinner A, et al. Genomic instability and myelodysplasia with monosomy 7 consequent to EVI1 activation after gene therapy for chronic granulomatous disease. Nat Med. 2010;16(2):198–204. doi: 10.1038/nm.2088 20098431
36. Braun CJ, Boztug K, Paruzynski A, Witzel M, Schwarzer A, Rothe M, et al. Gene therapy for Wiskott-Aldrich syndrome—long-term efficacy and genotoxicity. Sci Transl Med. 2014;6(227):227ra33.
37. Wu C, Dunbar CE. Stem cell gene therapy: the risks of insertional mutagenesis and approaches to minimize genotoxicity. Front Med. 2011;5(4):356–71. doi: 10.1007/s11684-011-0159-1 22198747
38. Hacein-Bey-Abina S, Pai SY, Gaspar HB, Armant M, Berry CC, Blanche S, et al. A modified gamma-retrovirus vector for X-linked severe combined immunodeficiency. N Engl J Med. 2014;371(15):1407–17. doi: 10.1056/NEJMoa1404588 25295500
39. El Ashkar S, Van Looveren D, Schenk F, Vranckx LS, Demeulemeester J, De Rijck J, et al. Engineering next-generation BET-independent MLV vectors for safer gene therapy. Mol Ther Nucleic Acids. 2017;7:231–45. doi: 10.1016/j.omtn.2017.04.002 28624199
40. Roth MJ. Mutational analysis of the carboxyl terminus of the Moloney murine leukemia virus integration protein. J Virol. 1991;65:2141–5. 2002557
41. Schneider WM, Wu D-t, Amin V, Aiyer S, Roth MJ. MuLV IN mutants responsive to HDAC inhibitors enhance transcription from unintegrated retroviral DNA. Virology. 2012;426(2):188–96. doi: 10.1016/j.virol.2012.01.034 22365328
42. Gupta SS, Maetzig T, Maertens GN, Sharif A, Rothe M, Weidner-Glunde M, et al. Bromo- and extraterminal domain chromatin regulators serve as cofactors for murine leukemia virus integration. J Virol. 2013;87(23):12721–36. doi: 10.1128/JVI.01942-13 24049186
43. O'Reilly L, Roth MJ. Second-site changes affect viability of amphotropic/ecotropic chimeric enveloped murine leukemia viruses. J Virol. 2000;74((2)):899–913. doi: 10.1128/jvi.74.2.899-913.2000 10623753
44. Moiani A, Miccio A, Rizzi E, Severgnini M, Pellin D, Suerth JD, et al. Deletion of the LTR enhancer/promoter has no impact on the integration profile of MLV vectors in human hematopoietic progenitors. PLoS One. 2013;8(1):e55721. doi: 10.1371/journal.pone.0055721 23383272
45. Peredo C, O'Reilly L, Gray K, Roth MJ. Characterization of chimeras between the ecotropic Moloney murine leukemia virus and the amphotropic 4070A envelope proteins. J Virol. 1996;70(5):3142–52. 8627794
46. Bamunusinghe D, Liu Q, Plishka R, Dolan MA, Skorski M, Oler AJ, et al. Recombinant origins of pathogenic and nonpathogenic mouse gammaretroviruses with polytropic host range. J Virol. 2017;91(21):e00855–17. doi: 10.1128/JVI.00855-17 28794032
47. Jern P, Stoye JP, Coffin JM. Role of APOBEC3 in genetic diversity among endogenous murine leukemia viruses. PLoS Genet. 2007;3(10):2014–22. doi: 10.1371/journal.pgen.0030183 17967065
48. Santoni FA, Hartley O, Luban J. Deciphering the code for retroviral integration target site selection. PLoS Comput Biol. 2010;6(11):e1001008. doi: 10.1371/journal.pcbi.1001008 21124862
49. LeRoy G, Chepelev I, DiMaggio PA, Blanco MA, Zee BM, Zhao K, et al. Proteogenomic characterization and mapping of nucleosomes decoded by Brd and HP1 proteins. Genome Biol. 2012;13(8):R68. doi: 10.1186/gb-2012-13-8-r68 22897906
50. Chapuy B, McKeown MR, Lin CY, Monti S, Roemer MG, Qi J, et al. Discovery and characterization of super-enhancer-associated dependencies in diffuse large B cell lymphoma. Cancer Cell. 2013;24(6):777–90. doi: 10.1016/j.ccr.2013.11.003 24332044
51. Gilroy KL, Terry A, Naseer A, de Ridder J, Allahyar A, Wang W, et al. Gamma-retrovirus integration marks cell type-specific cancer genes: A novel profiling tool in cancer genomics. PLoS One. 2016;11(4):e0154070. doi: 10.1371/journal.pone.0154070 27097319
52. Ranzani M, Annunziato S, Adams DJ, Montini E. Cancer gene discovery: exploiting insertional mutagenesis. Mol Cancer Res. 2013;11(10):1141–58. doi: 10.1158/1541-7786.MCR-13-0244 23928056
53. Touw IP, Erkeland SJ. Retroviral insertion mutagenesis in mice as a comparative oncogenomics tool to identify disease genes in human leukemia. Mol Ther. 2007;15(1):13–9. doi: 10.1038/sj.mt.6300040 17164770
54. Uren AG, Kool J, Berns A, van Lohuizen M. Retroviral insertional mutagenesis: past, present and future. Oncogene. 2005;24(52):7656–72. doi: 10.1038/sj.onc.1209043 16299527
55. Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011;473(7345):43–9. doi: 10.1038/nature09906 21441907
56. Ey P, Freeman N, Bela B, Li P, McInnes J. Nucleotide sequence of the murine leukemia virus amphotropic strain 4070A integrase (IN) coding region and comparative structural analysis of the inferred polypeptide. Arch Virol. 1997;142:1757–70. doi: 10.1007/s007050050195 9672635
57. Puglia J, Wang T, Smith-Snyder C, Cote M, Scher M, Pelletier J, et al. Revealing domain structure through linker-scanning analysis of the murine leukemia virus (MuLV) RNase H and MuLV and human immunodeficiency virus type 1 integrase proteins. J Virol. 2006;80(19):9497–510. doi: 10.1128/JVI.00856-06 16973554
58. Pattison JM, Wright JB, Cole MD. Retroviruses hijack chromatin loops to drive oncogene expression and highlight the chromatin architecture around proto-oncogenic loci. PLoS One. 2015;10(3):e0120256. doi: 10.1371/journal.pone.0120256 25799187
59. Arita K, Maeda-Kasugai Y, Ohshima K, Tsuzuki S, Suguro-Katayama M, Karube K, et al. Generation of mouse models of lymphoid neoplasm using retroviral gene transduction of in vitro-induced germinal center B and T cells. Exp Hematol. 2013;41(8):731–41 e9. doi: 10.1016/j.exphem.2013.04.001 23583576
60. Junk DJ, Cipriano R, Stampfer M, Jackson MW. Constitutive CCND1/CDK2 activity substitutes for p53 loss, or MYC or oncogenic RAS expression in the transformation of human mammary epithelial cells. PLoS One. 2013;8(2):e53776. doi: 10.1371/journal.pone.0053776 23390492
61. Nakagawa M, Tsuzuki S, Honma K, Taguchi O, Seto M. Synergistic effect of Bcl2, Myc and Ccnd1 transforms mouse primary B cells into malignant cells. Haematologica. 2011;96(9):1318–26. doi: 10.3324/haematol.2011.041053 21606168
62. Melamed A, Yaguchi H, Miura M, Witkover A, Fitzgerald TW, Birney E, et al. The human leukemia virus HTLV-1 alters the structure and transcription of host chromatin in cis. Elife. 2018;7: e36245. doi: 10.7554/eLife.36245 29941091
63. Satou Y, Miyazato P, Ishihara K, Yaguchi H, Melamed A, Miura M, et al. The retrovirus HTLV-1 inserts an ectopic CTCF-binding site into the human genome. Proc Natl Acad Sci U S A. 2016;113(11):3054–9. doi: 10.1073/pnas.1423199113 26929370
64. Ong CT, Corces VG. Enhancer function: new insights into the regulation of tissue-specific gene expression. Nat Rev Genet. 2011;12(4):283–93. doi: 10.1038/nrg2957 21358745
65. Ong CT, Corces VG. CTCF: an architectural protein bridging genome topology and function. Nat Rev Genet. 2014;15(4):234–46. doi: 10.1038/nrg3663 24614316
66. Tang Z, Luo OJ, Li X, Zheng M, Zhu JJ, Szalaj P, et al. CTCF-mediated human 3D genome architecture reveals chromatin topology for transcription. Cell. 2015;163(7):1611–27. doi: 10.1016/j.cell.2015.11.024 26686651
67. Ozaki T, Wu D, Sugimoto H, Nagase H, Nakagawara A. Runt-related transcription factor 2 (RUNX2) inhibits p53-dependent apoptosis through the collaboration with HDAC6 in response to DNA damage. Cell Death Dis. 2013;4(4):e610.
68. Westendorf JJ, Zaidi SK, Cascino JE, Kahler R, van Wijnen AJ, Lian JB, et al. Runx2 (Cbfa1, AML-3) interacts with histone deacetylase 6 and represses the p21(CIP1/WAF1) promoter. Mol Cell Biol. 2002;22(22):7982–92. doi: 10.1128/MCB.22.22.7982-7992.2002 12391164
69. Valenzuela-Fernandez A, Cabrero JR, Serrador JM, Sanchez-Madrid F. HDAC6: a key regulator of cytoskeleton, cell migration and cell-cell interactions. Trends Cell Biol. 2008;18(6):291–7. doi: 10.1016/j.tcb.2008.04.003 18472263
70. Ernst J, Kellis M. Large-scale imputation of epigenomic datasets for systematic annotation of diverse human tissues. Nat Biotechnol. 2015;33(4):364–76. doi: 10.1038/nbt.3157 25690853
71. Hoffman MM, Ernst J, Wilder SP, Kundaje A, Harris RS, Libbrecht M, et al. Integrative annotation of chromatin elements from ENCODE data. Nucleic Acids Res. 2013;41(2):827–41. doi: 10.1093/nar/gks1284 23221638
72. Poletti V, Mavilio F. Interactions between retroviruses and the host cell genome. Mol Ther Methods Clin Dev. 2018;8:31–41. doi: 10.1016/j.omtm.2017.10.001 29159201
73. Maskell DP, Renault L, Serrao E, Lesbats P, Matadeen R, Hare S, et al. Structural basis for retroviral integration into nucleosomes. Nature. 2015;523(7560):366–9. doi: 10.1038/nature14495 26061770
74. McGinty RK, Tan S. Recognition of the nucleosome by chromatin factors and enzymes. Current Opinion in Structural Biology. 2016;37:54–61. doi: 10.1016/j.sbi.2015.11.014 26764865
75. Wilson MD, Renault L, Maskell DP, Ghoneim M, Pye VE, Nans A, et al. Retroviral integration into nucleosomes through DNA looping and sliding along the histone octamer. Nat Commun. 2019;10(1):4138. doi: 10.1038/s41467-019-11557-3
76. Lesbats P, Serrao E, Maskell DP, Pye VE, O'Reilly N, Lindemann D, et al. Structural basis for spumavirus GAG tethering to chromatin. Proc Natl Acad Sci U S A. 2017;114(21):5509–14. doi: 10.1073/pnas.1621159114 28490494
77. Shaytan AK, Landsman D, Panchenko AR. Nucleosome adaptability conferred by sequence and structural variations in histone H2A-H2B dimers. Curr Opin Struct Biol. 2015;32:48–57. doi: 10.1016/j.sbi.2015.02.004 25731851
78. King SR, Berson BJ, Risser R. Mechanism of interaction between endogenous ecotropic murine leukemia viruses in (Balb/C X C57bl/6) hybrid-cells. Virology. 1988;162(1):1–11. doi: 10.1016/0042-6822(88)90388-1 2447699
79. Biasco L, Rothe M, Buning H, Schambach A. Analyzing the genotoxicity of retroviral vectors in hematopoietic cell gene therapy. Mol Ther Methods Clin Dev. 2018;8:21–30. doi: 10.1016/j.omtm.2017.10.002 29159200
80. Rothe M, Schambach A, Biasco L. Safety of gene therapy: new insights to a puzzling case. Curr Gene Ther. 2014;14(6):429–36. doi: 10.2174/1566523214666140918110905 25245088
81. Cavazza A, Cocchiarella F, Bartholomae C, Schmidt M, Pincelli C, Larcher F, et al. Self-inactivating MLV vectors have a reduced genotoxic profile in human epidermal keratinocytes. Gene Ther. 2013;20(9):949–57. doi: 10.1038/gt.2013.18 23615186
82. Nam JS, Lee JE, Lee KH, Yang Y, Kim SH, Bae GU, et al. Shifting retroviral Vector integrations away from transcriptional start sites via DNA-binding protein domain insertion into Integrase. Mol Ther Methods Clin Dev. 2019;12:58–70. doi: 10.1016/j.omtm.2018.11.001 30534579
83. Schneider WM, Brzezinski JD, Aiyer S, Malani N, Gyuricza M, Bushman FD, et al. Viral DNA tethering domains complement replication-defective mutations in the p12 protein of MuLV Gag. Proc Natl Acad Sci U S A. 2013;110(23):9487–92. doi: 10.1073/pnas.1221736110 23661057
84. Felkner RH, Roth MJ. Mutational analysis of the N-linked glycosylation sites of the SU envelope protein of Moloney murine leukemia virus. J Virol. 1992;66(7):4258–64. 1318404
85. McCutchan JH, Pagano JS. Enhancement of the infectivity of simian virus 40 deoxyribonucleic acid with diethylaminoethyl-dextran. J Natl Cancer Inst. 1968;41:351–7. 4299537
86. Wu DT, Aiyer S, Villanueva RA, Roth MJ. Development of an enzyme-linked immunosorbent assay based on the murine leukemia virus p30 capsid protein. J Virol Methods. 2013;193(2):332–6. doi: 10.1016/j.jviromet.2013.06.020 23810854
87. Ting YT, Wilson CA, Farrell KB, Chaudry GJ, Eiden MV. Simian sarcoma-associated virus fails to infect Chinese hamster cells despite the presence of functional gibbon ape leukemia virus receptors. J Virol. 1998;72(12):9453–8. 9811678
88. Valdivieso-Torres L, Sarangi A, Whidby J, Marcotrigiano J, Roth MJ. Role of cysteines in stabilizing the randomized receptor binding domains within feline leukemia virus Envelope proteins. J Virol. 2015;90(6):2971–80. doi: 10.1128/JVI.02544-15 26719270
89. Tanese N, Roth MJ, Goff SP. Analysis of retroviral pol gene products with antisera raised against fusion proteins produced in Escherichia coli. J Virol. 1986;59:328–40. 2426463
90. Uphoff CC, Lange S, Denkmann SA, Garritsen HS, Drexler HG. Prevalence and characterization of murine leukemia virus contamination in human cell lines. PLoS One. 2015;10(4):e0125622. doi: 10.1371/journal.pone.0125622 25927683
91. Zheng H, Jia H, Shankar A, Heneine W, Switzer WM. Detection of murine leukemia virus or mouse DNA in commercial RT-PCR reagents and human DNAs. PLoS One. 2011;6(12):e29050. doi: 10.1371/journal.pone.0029050 22205995
92. Stoye JP, Coffin SM. Polymorphism of murine endogenous proviruses revealed by using virus class-specific oligonucleotide probes. J Virol. 1988;62:168–75. 2824845
93. Serrao E, Cherepanov P, Engelman AN. Amplification, next-generation sequencing, and genomic DNA mapping of retroviral integration sites. J Vis Exp. 2016(109) doi: 10.3791/53840 27023428
94. Wang GG, Calvo KR, Pasillas MP, Sykes DB, Hacker H, Kamps MP. Quantitative production of macrophages or neutrophils ex vivo using conditional Hoxb8. Nature Methods. 2006;3(4):287–93. doi: 10.1038/nmeth865 16554834
95. Achuthan V, Perreira JM, Sowd GA, Puray-Chavez M, McDougall WM, Paulucci-Holthauzen A, et al. Capsid-CPSF6 interaction licenses nuclear HIV-1 trafficking to sites of viral DNA integration. Cell Host Microbe. 2018;24(3):392–404 e8. doi: 10.1016/j.chom.2018.08.002 30173955
96. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60. doi: 10.1038/nmeth.3317 25751142
97. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9. doi: 10.1093/bioinformatics/btp352 19505943
98. Maldarelli F, Wu X, Su L, Simonetti FR, Shao W, Hill S, et al. HIV latency. Specific HIV integration sites are linked to clonal expansion and persistence of infected cells. Science. 2014;345(6193):179–83. doi: 10.1126/science.1254194 24968937
99. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2. doi: 10.1093/bioinformatics/btq033 20110278
100. Wickham H. ggplot2. New York: Springer-Verlag; 2016.
101. Achuthan V, DeStefano JJ. Alternative divalent cations (Zn(2)(+), Co(2)(+), and Mn(2)(+)) are not mutagenic at conditions optimal for HIV-1 reverse transcriptase activity. BMC Biochem. 2015;16:12. doi: 10.1186/s12858-015-0041-x 25934642
102. Jia X, Lin X, Chen J. Linear and exponential TAIL-PCR: a method for efficient and quick amplification of flanking sequences adjacent to Tn5 transposon insertion sites. AMB Express. 2017;7(1):195. doi: 10.1186/s13568-017-0495-x 29098449
103. Gouet P, Robert X. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014;42(W1):W320–W4.
104. Pei J, Kim B-H, Grishin NV. PROMALS3D: a tool for multiple protein sequence and structure alignments. Nucleic Acids Res. 2008;36(7):2295–300. doi: 10.1093/nar/gkn072 18287115
105. Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–12. doi: 10.1002/jcc.20084 15264254
Štítky
Hygiena a epidemiologie Infekční lékařství LaboratořČlánek vyšel v časopise
PLOS Pathogens
2019 Číslo 12
- Farmakovigilanční studie perorálních antivirotik indikovaných v léčbě COVID-19
- 10 bodů k očkování proti COVID-19: stanovisko České společnosti alergologie a klinické imunologie ČLS JEP
Nejčtenější v tomto čísle
- Coxiella burnetii Type 4B Secretion System-dependent manipulation of endolysosomal maturation is required for bacterial growth
- IL-22 produced by type 3 innate lymphoid cells (ILC3s) reduces the mortality of type 2 diabetes mellitus (T2DM) mice infected with Mycobacterium tuberculosis
- The pandemic Escherichia coli sequence type 131 strain is acquired even in the absence of antibiotic exposure
- A role of hypoxia-inducible factor 1 alpha in Mouse Gammaherpesvirus 68 (MHV68) lytic replication and reactivation from latency