HIV latency and integration site placement in five cell-based models
- Scott Sherrill-Mix1Email author,
- Mary K Lewinski2,
- Marylinda Famiglietti3,
- Alberto Bosque4,
- Nirav Malani1,
- Karen E Ocwieja1,
- Charles C Berry5,
- David Looney2, 6,
- Liang Shan7,
- Luis M Agosto9,
- Matthew J Pace8,
- Robert F Siliciano7,
- Una O’Doherty8,
- John Guatelli2, 6,
- Vicente Planelles4 and
- Frederic D Bushman1
© Sherrill-Mix et al.; licensee BioMed Central Ltd. 2013
Received: 26 April 2013
Accepted: 12 August 2013
Published: 16 August 2013
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.
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.
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.
Highly active antiretroviral therapy (HAART) can suppress HIV-1 replication in infected patients, but the ability of HIV to persist as an inducible reservoir of latent proviruses [1–3] obstructs eradication of the virus and functional cure . These latent proviruses are long lived [5, 6] and relatively invisible to the immune system [2, 7]. The potential for even a single virus to restart infection despite successful antiviral therapy means that it may be necessary to eliminate all latent proviruses to eradicate HIV from an infected person.
After integration, a positive feedback loop of Tat transactivation appears to partition proviral gene activity into either of two stable states [8–10]—abundant Tat driving high proviral expression or little Tat leading to quiescent latency. Similar to the positional effect variegation observed in fruit fly chromosomal rearrangements [11, 12], studies on cell clones with single integrations show that differing integration sites can have large differences in proviral expression [13–15]. These data suggest that integration site location, along with the cellular environment [15–18], influences the balance between latency and proviral expression.
Associations between latency and genomic features have also been reported in collections of integration sites from cell culture models although the consistency of these effects across model systems and their relationships to latency in patients remains uncertain. Lewinski et al.  reported that proviruses integrated in gene deserts, alphoid repeats and highly expressed genes are more likely to have low expression. Shan et al.  reported an association between latency and integration in the same transcriptional orientation as host genes. Pace et al.  found that silent and expressed provirus integration sites differed in the abundance and expression levels of nearby genes, GC content, CpG islands and alphoid repeats. In model systems with defined integration sites, Lenasi et al.  reported decreased and Han et al.  reported increased viral transcription when the provirus is downstream of a highly expressed host gene.
Cell-based models of latency are important for many aspects of HIV research, including screening small molecules that can reverse latency and potentially allow eradication [24, 25]. Location-driven differences in expression are preserved even after DNA methyltransferase and histone deacetylase inhibitor treatments , which suggests that integration location has the potential to confound “shock and kill” anti-latency treatments [26, 27]. A greater understanding of the effects of integration site location on latency could thus affect antiretroviral development.
To search for features of integration site associated with latency, we generated a set of inducible and expressed integration sites using a primary central memory CD4 + T cell model of latency [28, 29], collected four previously reported integration site datasets and modeled the effects of genomic features near the integration site on the expression status of these proviruses. Although some genomic features associated with latency in individual models, no feature was consistently associated with proviral expression across all five cell culture models. However, closely neighboring proviruses within the same cellular model shared the same latency status much more often than expected by chance suggesting that chromosomal position of integration affects latency but that the mechanism remains unclear or differs between cell culture models. Thus these data help inform the design of experiments in HIV eradication research.
HIV-1 integration datasets from in vitro models of latency
Time of harvest after infection
Generation of expressed vs. silent/inducible
Silent/inducible unique sites
Expressed unique sites
HIV vector pEV731(LTR-Tat-IRES-GFP)
TNF α, GFP expression
Bcl-2 transduced CD4 +
Primary CD4 +T cells (Bcl-2transduced)
HIV NL4-3- Δ6-drEGFP(inactivated gag, vif,vpr, vpu, nef and envreplaced by GFP)
3 days + 3-4weeks + 3 days
Anti-CD3, anti-CD28 antibodies, GFP expression
Active CD4 +
Primary activeCD4 + T cells
High vs. low Gag
Resting CD4 +
Primary restingCD4 + T cells
High vs. low Gag
Central Memory CD4 +
Primary centralmemory CD4 +T cells
HIV NL4-3 ΔNef GFP
2 days/9 days
Anti-CD3, anti-CD28 antibodies, GFP expression
Genomic data available for comparison to HIV integration sites
T cell expression
% GC, HIV PWM score, distance to centrosome, chromosomal position
DNA, LINE, Low_complexity, LTR, Other, RC, RNA, rRNA, Satellite, scRNA, Simple_repeat, SINE, snRNA, srpRNA, tRNA, alphoid
H2AK5ac, H2AK9ac, H2BK120ac, H2BK12ac, H2BK20ac, H2BK5ac, H3K14ac, H3K18ac, H3K23ac, H3K27ac, H3K36ac, H3K4ac, H3K9ac, H4K12ac, H4K16ac, H4K5ac, H4K8ac, H4K91ac
Histone methylation and other proteins
CTCF, H2AZ, H2BK5me1, H3K27me1, H3K27me2, H3K27me3, H3K36me1, H3K36me3, H3K4me1, H3K4me2, H3K4me3, H3K79me1, H3K79me2, H3K79me3, H3K9me1, H3K9me2, H3K9me3, H3R2me1, H3R2me2, H4K20me1, H4K20me3, H4R3me2, PolII
HATs and HDACs
Resting-HDAC1, Resting-HDAC2, Resting-HDAC3, Resting-HDAC6, Resting-p300, Resting-CBP, Resting-MOF, Resting-PCAF, Resting-Tip60, Active-HDAC6, Active-Tip60
In gene, in gene (same strand), gene count, distance to nearest gene, in exon, in intron
To investigate whether a combination of variables may affect latency, we fit a lasso-regularized logistic regression, as implemented in the R package glmnet , to predict latency using the genomic variables. The relationship between silent/inducible status and each genomic variable was allowed to vary between models by including the interaction of genomic features with dummy variables indicating cellular model. The λ smoothing parameter of the lasso regression was optimized by finding the λ with lowest classification error in 480-fold cross validation and finding the simplest model with misclassification error within one standard error.
When each dataset was fit individually with leave-one-out cross validation, improvements in cross-validated misclassification error were only observed in the Active CD4 + (5.8% decrease in misclassification error, standard error: 2.1) and Jurkat (6.7% decrease in misclassification error, standard error: 3.5) samples (Figure 2B-F). There was no overlap in variables selected for the Active CD4 + and Jurkat samples.
Finding little global association between latency and genomic features, we investigated whether predictors of latency reported previously by single studies were consistently associated with latency across studies.
Model systems with defined integration sites show upstream transcription can interfere with viral transcription  and that cellular transcription in the same orientation may interfere with viral transcription  or increase viral transcription  and in opposite orientations may decrease transcription . In integration site studies, integration outside genes appears to increase latency  but high transcription of nearby host cell genes may cause increased latency [19, 20]. In addition, Tat or other viral proteins may affect cellular transcription [41, 42].
Lewinski et al.  also reported decreased latency near CpG islands and reasoned this was tied to the increased latency in gene deserts. In the Resting CD4 + sample, silent sites were on average further from CpG islands than expressed sites (Bonferroni corrected Welch’s two sample T test, p=0.006), but there was no significant relationship between silent/inducible status and log distance to CpG island after Bonferroni correction if the integration site’s location inside or outside of a gene was accounted for first (analysis of deviance).
Alphoid repeats are repetitive DNA sequences found largely in the heterochromatin of centromeres . Integration near heterochromatic alphoid repeats has been reported to associate with latency [14, 19, 21]. Looking only at uniquely mapping sites, there was no statistically significant association between latency and location inside an alphoid repeat in pooled or individual samples (Fisher’s exact test).
Histone marks or chromatin remodeling, especially involving the key “Nuc-1” histone near the transcription start site in the viral LTR, appear to affect viral expression [15, 45, 46]. Based on this effect, histone deacetylase inhibitors have been developed as potential HIV treatments and show some promise in disrupting latency . In these genome-wide datasets, we do not have information on the state of individual LTR nucleosomes. However, repressive chromatin does seem to spread to nearby locations if not blocked by insulators [11, 12] and the state of neighboring chromatin could affect proviral transcription independently of provirus-associated histones.
Although the appearance of several significantly associated acetylation marks might suggest acetylation exerts a considerable effect on the expression of a provirus, there are strong correlations among these marks, so their effects may not be independent. To account for the correlations between these variables, we performed a principal component analysis (PCA) to convert the correlated acetylation marks into a series of uncorrelated principal components that capture much of the variance within a few components. Here, the first principal component explained 59% of the variance and the first ten components 84%. Several of these principal components again displayed significant associations with latency in the Active CD4 +, Resting CD4 + and Central Memory CD4 + samples but no significant correlations in the Bcl-2 transduced CD4 + or Jurkat samples (Figure 7B). A logistic regression of expression status on the first ten principal components and sample did not reduce misclassification error from a base model including only sample in 480-fold cross validation (base model misclassification error: 36.4%, PCA model: 36.5%). This suggests that acetylation of neighboring chromatin does not exert strong effects on latency in all samples.
Here we compared the latency status of HIV-1 proviruses in five model systems with the genomic features surrounding their integration sites. Surprisingly, no relationships between genomic features near the integration location and latency achieved significance in all models. Proviruses from the same cellular model integrated in nearby positions did share the same latency status much more often than predicted by chance, indicating the existence of local features influencing latency, but these were not consistent among models. This suggests that whatever features are affecting latency are highly local and model-specific, and that we may not have access to all relevant chromosomal features (e.g. [47–50]).
In addition to differences in experimental conditions, methodological issues have the potential to obscure patterns. Examples include multiply infected cells, inactivated viruses and inaccurate assessment of HIV gene activity—each of these are discussed below.
A latent provirus integrated into the same cell as an expressed provirus will be erroneously sorted as expressed, potentially confounding analysis. A low multiplicity of infection (MOI) will help to avoid this problem, but there is still the potential for a significant proportion of the cells studied to contain multiple integrations. This problem arises because although cells with multiple integrations form a small proportion of total cells, most of the total are cells lacking an integrated provirus and thus are excluded by experimental design. For example, assuming integrations are Poisson distributed with an MOI of 0.1 (1 integration per 10 cells), 90.5% of cells will not contain a provirus, 9% of cells will contain one proviral integration and 0.5% of cells will contain multiple integrations. The cells without an integration are not amplified by HIV-targeted PCR leaving only 9.5% of the total cells. Of these cells actually under study, 4.9% will contain multiple integrations. Thus the signal from expressed proviruses may be muted by the presence of latent proviruses in the expressed population.
The replication cycle of HIV is error prone, and a significant proportion of virions contain mutated genomes . In studies that do not check for inducibility, mutant proviruses integrated in regions of the genome otherwise favorable to proviral expression can be sorted into the latent pool due to mutational inactivation. This problem of inactivated provirus is worse when latent provirus are rare and exacerbated further when looking at latency in the cells of HIV patients due to selective enrichment of inactivated proviruses incapable of spreading infection . Here, the effects of mutation are minimized in the datasets that required inducible viral expression (Jurkat, Bcl-2 transduced CD4 +, Central Memory CD4 +) but may be a confounder in the two datasets that were sorted based on lack of viral expression only (Active CD4 +, Resting CD4 +).
Inaccurate staining or leaky markers may also result in misclassification of proviruses. False positives and false negatives will result in incorrectly sorted latent and expressed integrations. For example, if 5% of cells not containing Gag are labeled as Gag+ and there are an equal amount of latent and expressed integration sites, then 4.8% of integrations labeled expressed will actually be latent. If a category is rare, false staining has even greater potential to cause error. For example, if only 5% of sites are latent and a Gag stain has a false negative rate of 5%, then we would expect 48.7% of sites classified as latent to actually be mislabeled expressed integrations.
Attempts to induce latent proviruses in patients have so far focused on using histone deacetylase inhibitors, raising interest in associations with histone acetylation in these data. An important caveat in results from these genome-wide data is that histone modification near the integrated provirus may not be representative of modification within the provirus at the key “Nuc-1” nucleosome of the transcription start site , though local correlations in chromatin states are well established from studies of position effect variegation [11, 12]. We found that some histone acetylation marks were significantly associated with viral expression in some but not all samples (Figures 1 and 7). This lack of association may be due to a lack of power in these studies, but the confidence intervals suggest that any correlations between acetylations and latency are unlikely to be strong. These weak correlations raise the possibility that there are populations of latent proviruses that are not associated with acetylation and may not be inducible by histone deacetylase inhibitors.
This study highlights that the choice of model system can have a large effect on measurements of latency. Further studies are needed to determine which in vitro models best reflect latency in vivo. Different cell culture models may report genuinely different mechanisms of latency. While we did see some relationship between histone acetylation and latency, paralleling a recent clinical trial of SAHA , associations with histone acetylation did not explain a large fraction of the difference between latent and expressed proviruses in any of the five models. One possible explanation is that there may be multiple mechanisms that maintain proviruses in a latent state. To be successful, shock-and-kill treatments must induce and destroy all latent proviruses to eliminate HIV from an infected individual, raising the question of whether multiple simultaneous inducing treatments will be necessary.
Availability of supporting data
Sequence reads from the Central Memory CD4 + sample reported here, the Resting CD4 + and Active CD4 + data reported by Pace et al. , the Bcl-2 transduced CD4 + data reported by Shan et al.  and reprocessed data originally reported by Lewinski et al.  are available at the Sequence Read Archive under accession number SRP028573.
Naive CD4 + T cells were purified by negative selection from peripheral blood mononuclear cells. The cells were activated with anti-CD3 and anti-CD28 (+TGF-beta, anti-IL-12, and anti-IL-4) to generate “non-polarized” cells (the in vitro equivalent of central memory T cells). Five days after isolation, cells were infected with an NL4-3-based virus with GFP in place of Nef and the LAI envelope (X4) provided in trans at a concentration of 500 ng of p24 as measured by ELISA per million cells. Based on previous experience with this model, this amount of p24 should produce an MOI of approximately 0.15. Cells were cultured in the presence of IL-2. Two days post-infection, cells were sorted for GFP+; this active population expresses GFP even when treated with flavopiridol, although for this study they were not treated. The inducible population was the set of GFP negative cells from the initial sort that, 9 days post-infection, were activated with anti-CD3 and anti-CD28 and sorted for GFP production.
Genomic DNA from the inducible and expressed populations was digested with MseI, ligated to an adapter, and amplified by ligation-mediated PCR essentially as in Wu et al.  and Mitchell et al.  except that the nested PCR primers included sequence for the Ion Torrent P1 adapter and adapter A sequence with a 5 base barcode sequence specific to the inducible or expressed conditions. Amplicons were sequenced using an Ion Torrent Personal Genome Machine (PGM) according to manufacturer’s instructions using an Ion 316 chip and the Ion PGM 200 Sequencing kit (Life Technologies). The sequence reads were sorted into samples by barcode. All reads were required to match the expected 5′ sequence with a Levenshtein edit distance less than 3 from the expected barcode, 5′ primer and HIV long terminal repeat (LTR). The 5′ primer and HIV sequence, along with the 3′ primer if present, were trimmed from the read. Sequences with less than 24 bases remaining or containing any eight base window with an average quality less than 15 were discarded. Duplicate reads and reads forming an exact substring of a longer read were removed.
Previously published data
We collected integration sites from three previously reported studies (Table 1), for a total of four expressed versus silent/inducible pairs of samples. These studies used primary CD4 + T cells or Jurkat cells infected with HIV or HIV-derived constructs as cell culture models of latency. Flow cytometry allowed cells expressing viral encoded proteins to be sorted from non-expressing cells. In two of the studies, these non-expressing populations were stimulated to ensure that the provirus could be aroused from latency. Specific differences in protocol between the study sets are summarized below.
Lewinski et al.  infected Jurkat cells with a VSV-G pseudotyped, GFP-expressing pEV731 HIV construct (LTR-Tat-IRES-GFP)  at an MOI of 0.1. The cells were sorted into GFP+ and GFP- two to four days after infection. GFP+ cells were sorted again two weeks after infection and cells that were again GFP+ were collected for integration site sequencing. GFP- cells were sorted for GFP negativity twice more then stimulated with TNFalpha. Cells that were GFP+ after stimulation were collected for integration site sequencing. DNA was digested with MseI or a combination of NheI, SpeI and XbaI, ligated to adapters for nested PCR, amplified and sequenced by Sanger capillary electrophoresis.
Bcl-2 transduced CD4 +
Shan et al.  transduced CD4 + T cells with Bcl-2, costimulated with bound anti-CD3 and soluble anti-CD28 antibodies, interleukin-2 and T cell growth factor and then infected with X4-pseudotyped GFP-expressing NL4-3- Δ6-drEGFP construct  at an MOI of less than 0.1. DNA was extracted, digested with PstI and circularized . HIV-human junctions were amplified by reverse PCR and sequenced using Sanger capillary electrophoresis.
Active CD4 + & Resting CD4 +
Pace et al.  spinoculated CD4 + T cells with HIV NL4-3 at an MOI of 0.1. After 96 hours, the cells were stained for intracellular Gag CD25, CD69 and HLA-DR and sorted into four subpopulations based on activation state and Gag expression; activated Gag-, activated Gag+, resting Gag- and resting Gag+. The ability of the viruses to reactivate was not tested although previous studies have shown that the majority are likely inducible . Genomic DNA was extracted and digested with restriction enzymes MseI and Tsp509 and ligated to adapters. Proviral LTR-host genome junctions were sequenced by 454 pyrosequencing after nested PCR.
All datasets were processed using the hiReadsProcessor R package . Adaptor trimmed reads were aligned to UCSC freeze hg19 using BLAT . Genomic alignments were scored and required to start within the first three bases of a read with 98% identity. Alignments for a given read with a BLAT score less than the maximum score for that read were discarded. Reads giving rise to multiple best scoring genomic alignments were excluded, while reads with a single best hit were dereplicated and converged if within 5 bp of each other. The Bcl-2 transduced CD4 + sample was sequenced from U3 in the 5′ HIV LTR while the other samples were sequenced from U5 in the 3′ LTR. To account for the 5 base duplication of host DNA caused by HIV integration, the chromosomal coordinates of the Bcl-2 transduced CD4 + sample were adjusted by ±4 bases.
To allow for alignment difficulties in the analysis of genomic repeats, reads with multiple best scoring alignments, along with the single best hit reads used above, were included in the repeat analyses. If any best scoring alignment for a read fell within a repeat, then that read was considered to map to that repeat.
A total of 140 whole genome features for CD4 + T-cells were gathered from data sources indicated in Table 2. For features encoded as peaks or hotspots, the log of the distance of each integration site to the nearest border was used for modeling. Integration sites from HIV 89.6 infection in primary CD4 + T cells (unpublished data) were used to count nearby integrations and determine a ±20 bp position weight matrix for integration targets. Illumina RNA-Seq from active CD4 + cells (unpublished data) was used to estimate raw cellular expression and fragments per kilobase of transcript per million mapped reads for genes as calculated by Cufflinks . For sequence-based data like RNA-Seq and ChIP-Seq, the number of reads aligned within a ± 50, 500, 5,000 50,000 and 500,000 bp windows of each integration site were counted and log transformed. In addition, chromatin state classifications derived from a hidden Markov model based on histone marks and a few binding factors  were included as binary variables. All data from previous genomic freezes were converted to hg19 using liftover .
All statistical analysis was performed in R 2.15.2 . The analyses are described in a reproducible report (Additional file 1). The annotated integration site data necessary to perform the analyses (Additional files 2 and 3) and the compilable code (Additional file 4) to generate this reproducible report are provided as supplemental information. The new Central Memory CD4 + data set was analyzed as in Berry et al.  (Additional file 5). The integration patterns appeared similar to previously reported HIV integration site datasets .
We would like to thank Werner Witke for assistance with IonTorrent sequencing. This work was supported in part by NIH grants R01 AI 052845-11 to FDB, R21AI 096993 and K02AI078766 to UO’D, 5T32HG000046 to SS-M, AI087508 to VP and R01AI038201 to JG, the Penn Genome Frontiers Institute, the University of Pennsylvania Center for AIDS Research (CFAR) P30 AI 045008, the University of California, San Diego, CFAR P30 AI036214 and the Collaboratory of AIDS Researchers for Eradication U19AI096113 (R.F.S., V.P., project leaders).
- Chun TW, Finzi D, Margolick J, Chadwick K, Schwartz D, Siliciano RF: In vivo fate of HIV-1-infected T cells: quantitative analysis of the transition to stable latency. Nat Med. 1995, 1 (12): 1284-1290. 10.1038/nm1295-1284. http://www.nature.com/nm/journal/v1/n12/full/nm1295-1284.html,View ArticlePubMedGoogle Scholar
- Chun TW, Carruth L, Finzi D, Shen X, DiGiuseppe JA, Taylor H, Hermankova M, Chadwick K, Margolick J, Quinn TC, Kuo YH, Brookmeyer R, Zeiger MA, Barditch-Crovo P, Siliciano RF: Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature. 1997, 387 (6629): 183-188. 10.1038/387183a0. http://dx.doi.org/10.1038/387183a0,View ArticlePubMedGoogle Scholar
- Davey RT, Bhat N, Yoder C, Chun TW, Metcalf JA, Dewar R, Natarajan V, Lempicki RA, Adelsberger JW, Miller KD, Kovacs JA, Polis MA, Walker RE, Falloon J, Masur H, Gee D, Baseler M, Dimitrov DS, Fauci AS, Lane HC: HIV-1 and T cell dynamics after interruption of highly active antiretroviral therapy (HAART) in patients with a history of sustained viral suppression. Proc Natl Acad Sci U S A. 1999, 96 (26): 15109-15114. 10.1073/pnas.96.26.15109. http://dx.doi.org/10.1073/pnas.96.26.15109,PubMed CentralView ArticlePubMedGoogle Scholar
- Richman DD, Margolis DM, Delaney M, Greene WC, Hazuda D, Pomerantz RJ: The challenge of finding a cure for HIV infection. Science. 2009, 323 (5919): 1304-1307. 10.1126/science.1165706. http://dx.doi.org/10.1126/science.1165706,View ArticlePubMedGoogle Scholar
- Finzi D, Blankson J, Siliciano JD, Margolick JB, Chadwick K, Pierson T, Smith K, Lisziewicz J, Lori F, Flexner C, Quinn TC, Chaisson RE, Rosenberg E, Walker B, Gange S, Gallant J, Siliciano RF: Latent infection of CD4+ T cells provides a mechanism for lifelong persistence of HIV-1, even in patients on effective combination therapy. Nat Med. 1999, 5 (5): 512-517. 10.1038/8394. http://dx.doi.org/10.1038/8394,View ArticlePubMedGoogle Scholar
- Siliciano JD, Kajdas J, Finzi D, Quinn TC, Chadwick K, Margolick JB, Kovacs C, Gange SJ, Siliciano RF: Long-term follow-up studies confirm the stability of the latent reservoir for HIV-1 in resting CD4+ T cells. Nat Med. 2003, 9 (6): 727-728. 10.1038/nm880. http://dx.doi.org/10.1038/nm880,View ArticlePubMedGoogle Scholar
- Finzi D, Hermankova M, Pierson T, Carruth LM, Buck C, Chaisson RE, Quinn TC, Chadwick K, Margolick J, Brookmeyer R, Gallant J, Markowitz M, Ho DD, Richman DD, Siliciano RF: Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. Science. 1997, 278 (5341): 1295-1300. 10.1126/science.278.5341.1295. http://dx.doi.org/10.1126/science.278.5341.1295,View ArticlePubMedGoogle Scholar
- Weinberger LS, Dar RD, Simpson ML: Transient-mediated fate determination in a transcriptional circuit of HIV. Nat Genet. 2008, 40 (4): 466-470. 10.1038/ng.116. http://dx.doi.org/10.1038/ng.116,View ArticlePubMedGoogle Scholar
- Singh A, Razooky B, Cox CD, Simpson ML, Weinberger LS: Transcriptional bursting from the HIV-1 promoter is a significant source of stochastic noise in HIV-1 gene expression. Biophys J. 2010, 98 (8): L32-L34. 10.1016/j.bpj.2010.03.001. http://dx.doi.org/10.1016/j.bpj.2010.03.001,PubMed CentralView ArticlePubMedGoogle Scholar
- Razooky BS, Weinberger LS: Mapping the architecture of the HIV-1 Tat circuit: A decision-making circuit that lacks bistability and exploits stochastic noise. Methods. 2011, 53: 68-77. 10.1016/j.ymeth.2010.12.006. http://dx.doi.org/10.1016/j.ymeth.2010.12.006,PubMed CentralView ArticlePubMedGoogle Scholar
- Muller HJ: Types of visible variations induced by X-rays in Drosophila. J Genet. 1930, 22: 299-334. 10.1007/BF02984195. http://www.ias.ac.in/j_archive/jgenet/22/vol22contents.html,View ArticleGoogle Scholar
- Gaszner M, Felsenfeld G: Insulators: exploiting transcriptional and epigenetic mechanisms. Nat Rev Genet. 2006, 7 (9): 703-713. 10.1038/nrg1925. http://dx.doi.org/10.1038/nrg1925,View ArticlePubMedGoogle Scholar
- Jordan A, Defechereux P, Verdin E: The site of HIV-1 integration in the human genome determines basal transcriptional activity and response to Tat transactivation. EMBO J. 2001, 20 (7): 1726-1738. 10.1093/emboj/20.7.1726. http://dx.doi.org/10.1093/emboj/20.7.1726,PubMed CentralView ArticlePubMedGoogle Scholar
- Jordan A, Bisgrove D, Verdin E: HIV reproducibly establishes a latent infection after acute infection of T cells in vitro. EMBO J. 2003, 22 (8): 1868-1877. 10.1093/emboj/cdg188. http://dx.doi.org/10.1093/emboj/cdg188,PubMed CentralView ArticlePubMedGoogle Scholar
- Pearson R, Kim YK, Hokello J, Lassen K, Friedman J, Tyagi M, Karn J: Epigenetic silencing of human immunodeficiency virus (HIV) transcription by formation of restrictive chromatin structures at the viral long terminal repeat drives the progressive entry of HIV into latency. J Virol. 2008, 82 (24): 12291-12303. 10.1128/JVI.01383-08. http://dx.doi.org/10.1128/JVI.01383-08,PubMed CentralView ArticlePubMedGoogle Scholar
- Romerio F, Gabriel MN, Margolis DM: Repression of human immunodeficiency virus type 1 through the novel cooperation of human factors YY1 and LSF. J Virol. 1997, 71 (12): 9375-9382. http://jvi.asm.org/content/71/12/9375.long,PubMed CentralPubMedGoogle Scholar
- Coull JJ, Romerio F, Sun JM, Volker JL, Galvin KM, Davie JR, Shi Y, Hansen U, Margolis DM: The human factors YY1 and LSF repress the human immunodeficiency virus type 1 long terminal repeat via recruitment of histone deacetylase 1. J Virol. 2000, 74 (15): 6790-6799. 10.1128/JVI.74.15.6790-6799.2000. http://jvi.asm.org/content/74/15/6790,PubMed CentralView ArticlePubMedGoogle Scholar
- He G, Margolis DM: Counterregulation of chromatin deacetylation and histone deacetylase occupancy at the integrated promoter of human immunodeficiency virus type 1 (HIV-1) by the HIV-1 repressor YY1 and HIV-1 activator Tat. Mol Cell Biol. 2002, 22 (9): 2965-2973. 10.1128/MCB.22.9.2965-2973.2002. http://dx.doi.org/10.1128/?MCB.22.9.2965-2973.2002,PubMed CentralView ArticlePubMedGoogle Scholar
- Lewinski MK, Bisgrove D, Shinn P, Chen H, Hoffmann C, Hannenhalli S, Verdin E, Berry CC, Ecker JR, Bushman FD: Genome-wide analysis of chromosomal features repressing human immunodeficiency virus transcription. J Virol. 2005, 79 (11): 6610-6619. 10.1128/JVI.79.11.6610-6619.2005. http://dx.doi.org/10.1128/JVI.79.11.6610-6619.2005,PubMed CentralView ArticlePubMedGoogle Scholar
- Shan L, Yang HC, Rabi SA, Bravo HC, Shroff NS, Irizarry RA, Zhang H, Margolick JB, Siliciano JD, Siliciano RF: Influence of host gene transcription level and orientation on HIV-1 latency in a primary-cell model. J Virol. 2011, 85 (11): 5384-5393. 10.1128/JVI.02536-10. http://dx.doi.org/10.1128/JVI.02536-10,PubMed CentralView ArticlePubMedGoogle Scholar
- Pace MJ, Graf EH, Agosto LM, Mexas AM, Male F, Brady T, Bushman FD, O’Doherty U: Directly infected resting CD4+ T cells can produce HIV Gag without spreading infection in a model of HIV latency. PLoS Pathog. 2012, 8 (7): e1002818-10.1371/journal.ppat.1002818. http://dx.doi.org/10.1371/journal.ppat.1002818,PubMed CentralView ArticlePubMedGoogle Scholar
- Lenasi T, Contreras X, Peterlin BM: Transcriptional interference antagonizes proviral gene expression to promote HIV latency. Cell Host Microbe. 2008, 4 (2): 123-133. 10.1016/j.chom.2008.05.016. http://dx.doi.org/10.1016/j.chom.2008.05.016,PubMed CentralView ArticlePubMedGoogle Scholar
- Han Y, Lin YB, An W, Xu J, Yang HC, O’Connell K, Dordai D, Boeke JD, Siliciano JD, Siliciano RF: Orientation-dependent regulation of integrated HIV-1 expression by host gene transcriptional readthrough. Cell Host Microbe. 2008, 4 (2): 134-146. 10.1016/j.chom.2008.06.008. http://dx.doi.org/10.1016/j.chom.2008.06.008,PubMed CentralView ArticlePubMedGoogle Scholar
- Shan L, Deng K, Shroff NS, Durand CM, Rabi SA, Yang HC, Zhang H, Margolick JB, Blankson JN, Siliciano RF: Stimulation of HIV-1-specific cytolytic T lymphocytes facilitates elimination of latent viral reservoir after virus reactivation. Immunity. 2012, 36 (3): 491-501. 10.1016/j.immuni.2012.01.014. http://dx.doi.org/10.1016/j.immuni.2012.01.014,PubMed CentralView ArticlePubMedGoogle Scholar
- Boehm D, Calvanese V, Dar RD, Xing S, Schroeder S, Martins L, Aull K, Li PC, Planelles V, Bradner JE, Zhou MM, Siliciano RF, Weinberger L, Verdin E, Ott M: BET bromodomain-targeting compounds reactivate HIV from latency via a Tat-independent mechanism. Cell Cycle. 2013, 12 (3): 452-462. 10.4161/cc.23309. http://dx.doi.org/10.4161/cc.23309,PubMed CentralView ArticlePubMedGoogle Scholar
- Savarino A, Mai A, Norelli S, Daker SE, Valente S, Rotili D, Altucci L, Palamara AT, Garaci E: “Shock and kill” effects of class I-selective histone deacetylase inhibitors in combination with the glutathione synthesis inhibitor buthionine sulfoximine in cell line models for HIV-1 quiescence. Retrovirology. 2009, 6: 52-10.1186/1742-4690-6-52. http://dx.doi.org/10.1186/1742-4690-6-52,PubMed CentralView ArticlePubMedGoogle Scholar
- Archin NM, Liberty AL, Kashuba AD, Choudhary SK, Kuruc JD, Crooks AM, Parker DC, Anderson EM, Kearney MF, Strain MC, Richman DD, Hudgens MG, Bosch RJ, Coffin JM, Eron JJ, Hazuda DJ, Margolis DM: Administration of vorinostat disrupts HIV-1 latency in patients on antiretroviral therapy. Nature. 2012, 487 (7408): 482-485. 10.1038/nature11286. http://dx.doi.org/10.1038/nature11286,PubMed CentralView ArticlePubMedGoogle Scholar
- Bosque A, Planelles V: Induction of HIV-1 latency and reactivation in primary memory CD4+ T cells. Blood. 2009, 113: 58-65. 10.1182/blood-2008-07-168393. http://dx.doi.org/10.1182/blood-2008-07-168393,PubMed CentralView ArticlePubMedGoogle Scholar
- Bosque A, Planelles V: Studies of HIV-1 latency in an ex vivo model that uses primary central memory T cells. Methods. 2011, 53: 54-61. 10.1016/j.ymeth.2010.10.002. http://dx.doi.org/10.1016/j.ymeth.2010.10.002,PubMed CentralView ArticlePubMedGoogle Scholar
- Plesa G, Dai J, Baytop C, Riley JL, June CH, O’Doherty U: Addition of deoxynucleosides enhances human immunodeficiency virus type 1 integration and 2LTR formation in resting CD4+ T cells. J Virol. 2007, 81 (24): 13938-13942. 10.1128/JVI.01745-07. http://dx.doi.org/10.1128/JVI.01745-07,PubMed CentralView ArticlePubMedGoogle Scholar
- Rosenbloom KR, Sloan CA, Malladi VS, Dreszer TR, Learned K, Kirkup VM, Wong MC, Maddren M, Fang R, Heitner SG, Lee BT, Barber GP, Harte RA, Diekhans M, Long JC, Wilder SP, Zweig AS, Karolchik D, Kuhn RM, Haussler D, Kent WJ: ENCODE Data in the UCSC Genome Browser: year 5 update. Nucleic Acids Res. 2013, 41 (D1): D56—D63-http://dx.doi.org/10.1093/nar/gks1172,PubMed CentralView ArticlePubMedGoogle Scholar
- Han J, Park SG, Bae JB, Choi J, Lyu JM, Park SH, Kim HS, Kim YJ, Kim S, Kim TY: The characteristics of genome-wide DNA methylation in naïve CD4+ T cells of patients with psoriasis or atopic dermatitis. Biochem Biophys Res Commun. 2012, 422: 157-163. 10.1016/j.bbrc.2012.04.128. http://dx.doi.org/10.1016/j.bbrc.2012.04.128,View ArticlePubMedGoogle Scholar
- Meyer LR, Zweig AS, Hinrichs AS, Karolchik D, Kuhn RM, Wong M, Sloan CA, Rosenbloom KR, Roe G, Rhead B, Raney BJ, Pohl A, Malladi VS, Li CH, Lee BT, Learned K, Kirkup V, Hsu F, Heitner S, Harte RA, Haeussler M, Guruvadoo L, Goldman M, Giardine BM, Fujita PA, Dreszer TR, Diekhans M, Cline MS, Clawson H, Barber GP, Haussler D, Kent WJ: The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013, 41 (D1): D64—D69-http://dx.doi.org/10.1093/nar/gks1048,PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Z, Zang C, Rosenfeld JA, Schones DE, Barski A, Cuddapah S, Cui K, Roh TY, Peng W, Zhang MQ, Zhao K: Combinatorial patterns of histone acetylations and methylations in the human genome. Nat Genet. 2008, 40 (7): 897-903. 10.1038/ng.154. http://dx.doi.org/10.1038/ng.154,PubMed CentralView ArticlePubMedGoogle Scholar
- Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009. http://dx.doi.org/10.1016/j.cell.2007.05.009,View ArticlePubMedGoogle Scholar
- Wang Z, Zang C, Cui K, Schones DE, Barski A, Peng W, Zhao K: Genome-wide mapping of HATs and HDACs reveals distinct functions in active and inactive genes. Cell. 2009, 138 (5): 1019-1031. 10.1016/j.cell.2009.06.049. http://dx.doi.org/10.1016/j.cell.2009.06.049,PubMed CentralView ArticlePubMedGoogle Scholar
- Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, Wei G, Zhao K: Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132 (5): 887-898. 10.1016/j.cell.2008.02.022. http://dx.doi.org/10.1016/j.cell.2008.02.022,View ArticlePubMedGoogle Scholar
- Hsu F, Kent WJ, Clawson H, Kuhn RM, Diekhans M, Haussler D: The UCSC Known Genes. Bioinformatics. 2006, 22 (9): 1036-1046. 10.1093/bioinformatics/btl048. http://dx.doi.org/10.1093/bioinformatics/btl048,View ArticlePubMedGoogle Scholar
- Friedman J, Hastie T, Tibshirani R: Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010, 33: 1-22. http://www.jstatsoft.org/v33/i01/,PubMed CentralView ArticlePubMedGoogle Scholar
- Greger IH, Demarchi F, Giacca M, Proudfoot NJ: Transcriptional interference perturbs the binding of Sp1 to the HIV-1 promoter. Nucleic Acids Res. 1998, 26 (5): 1294-1301. 10.1093/nar/26.5.1294. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC147389/,PubMed CentralView ArticlePubMedGoogle Scholar
- De Marco A, Biancotto C, Knezevich A, Maiuri P, Vardabasso C, Marcello A: Intragenic transcriptional cis-activation of the human immunodeficiency virus 1 does not result in allele-specific inhibition of the endogenous gene. Retrovirology. 2008, 5: 98-10.1186/1742-4690-5-98. http://dx.doi.org/10.1186/1742-4690-5-98,PubMed CentralView ArticlePubMedGoogle Scholar
- Chang ST, Sova P, Peng X, Weiss J, Law GL, Palermo RE, Katze MG: Next-generation sequencing reveals HIV-1-mediated suppression of T cell activation and RNA processing and regulation of noncoding RNA expression in a CD4+ T cell line. MBio. 2011, 2 (5): http://dx.doi.org/10.1128/mBio.00134-11,Google Scholar
- Waye JS, Willard HF: Nucleotide sequence heterogeneity of alpha satellite repetitive DNA: a survey of alphoid sequences from different human chromosomes. Nucleic Acids Res. 1987, 15 (18): 7549-7569. 10.1093/nar/15.18.7549.PubMed CentralView ArticlePubMedGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110 (1-4): 462-467. 10.1159/000084979. http://dx.doi.org/10.1159/000084979,View ArticlePubMedGoogle Scholar
- Verdin E, Paras P, Van, Lint C: Chromatin disruption in the promoter of human immunodeficiency virus type 1 during transcriptional activation. EMBO J. 1993, 12 (8): 3249-3259. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC413592/,PubMed CentralPubMedGoogle Scholar
- Van Lint C, Emiliani S, Ott M, Verdin E: Transcriptional activation and chromatin remodeling of the HIV-1 promoter in response to histone acetylation. EMBO J. 1996, 15 (5): 1112-1120. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC450009/,PubMed CentralPubMedGoogle Scholar
- Lassen KG, Ramyar KX, Bailey JR, Zhou Y, Siliciano RF: Nuclear retention of multiply spliced HIV-1 RNA in resting CD4+ T cells. PLoS Pathog. 2006, 2 (7): e68-10.1371/journal.ppat.0020068. http://dx.doi.org/10.1371/journal.ppat.0020068,PubMed CentralView ArticlePubMedGoogle Scholar
- Dieudonné M, Maiuri P, Biancotto C, Knezevich A, Kula A, Lusic M, Marcello A: Transcriptional competence of the integrated HIV-1 provirus at the nuclear periphery. EMBO J. 2009, 28 (15): 2231-2243. 10.1038/emboj.2009.141. http://dx.doi.org/10.1038/emboj.2009.141,PubMed CentralView ArticlePubMedGoogle Scholar
- Siliciano RF, Greene WC: HIV Latency. Cold Spring Harb Perspect Med. 2011, 1: a007096-10.1101/cshperspect.a007096. http://dx.doi.org/10.1101/cshperspect.a007096,PubMed CentralView ArticlePubMedGoogle Scholar
- Lusic M, Marini B, Ali H, Lucic B, Luzzati R, Giacca M: Proximity to PML nuclear bodies regulates HIV-1 latency in CD4+ T cells. Cell Host Microbe. 2013, 13 (6): 665-677. 10.1016/j.chom.2013.05.006. http://dx.doi.org/10.1016/j.chom.2013.05.006,View ArticlePubMedGoogle Scholar
- Mansky LM, Temin HM: Lower in vivo mutation rate of human immunodeficiency virus type 1 than that predicted from the fidelity of purified reverse transcriptase. J Virol. 1995, 69 (8): 5087-5094. http://jvi.asm.org/content/69/8/5087.short,PubMed CentralPubMedGoogle Scholar
- Wu X, Li Y, Crise B, Burgess SM: Transcription start regions in the human genome are favored targets for MLV integration. Science. 2003, 300 (5626): 1749-1751. 10.1126/science.1083413. http://dx.doi.org/10.1126/science.1083413,View ArticlePubMedGoogle Scholar
- Mitchell RS, Beitzel BF, Schroder ARW, Shinn P, Chen H, Berry CC, Ecker JR, Bushman FD: Retroviral DNA integration: ASLV, HIV, and MLV show distinct target site preferences. PLoS Biol. 2004, 2 (8): e234-10.1371/journal.pbio.0020234. http://dx.doi.org/10.1371/journal.pbio.0020234,PubMed CentralView ArticlePubMedGoogle Scholar
- Mochizuki H, Schwartz JP, Tanaka K, Brady RO, Reiser J: High-titer human immunodeficiency virus type 1-based vector systems for gene delivery into nondividing cells. J Virol. 1998, 72 (11): 8873-8883. http://jvi.asm.org/content/72/11/8873.abstract,PubMed CentralPubMedGoogle Scholar
- Han Y, Lassen K, Monie D, Sedaghat AR, Shimoji S, Liu X, Pierson TC, Margolick JB, Siliciano RF, Siliciano JD: Resting CD4+ T cells from human immunodeficiency virus type 1 (HIV-1)-infected individuals carry integrated HIV-1 genomes within actively transcribed host genes. J Virol. 2004, 78 (12): 6122-6133. 10.1128/JVI.78.12.6122-6133.2004. http://dx.doi.org/10.1128/JVI.78.12.6122-6133.2004,PubMed CentralView ArticlePubMedGoogle Scholar
- Malani N: hiReadsProcessor R package. http://github.com/malnirav/hiReadsProcessor,
- Kent WJ: BLAT–the BLAST-like alignment tool. Genome Res. 2002, 12 (4): 656-664. http://dx.doi.org/10.1101/gr.229202,PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-515. 10.1038/nbt.1621. http://dx.doi.org/10.1038/nbt.1621,PubMed CentralView ArticlePubMedGoogle Scholar
- Ernst J, Kellis M: Discovery and characterization of chromatin states for systematic annotation of the human genome. Nat Biotechnol. 2010, 28 (8): 817-825. 10.1038/nbt.1662. http://dx.doi.org/10.1038/nbt.1662,PubMed CentralView ArticlePubMedGoogle Scholar
- Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, Diekhans M, Furey TS, Harte RA, Hsu F, Hillman-Jackson J, Kuhn RM, Pedersen JS, Pohl A, Raney BJ, Rosenbloom KR, Siepel A, Smith KE, Sugnet CW, Sultan-Qurraie A, Thomas DJ, Trumbower H, Weber RJ, Weirauch M, Zweig AS, Haussler D, Kent WJ: The UCSC genome browser database: update 2006. Nucleic Acids Res. 2006, 34 (Database issue): D590-D598. http://dx.doi.org/10.1093/nar/gkj144,PubMed CentralView ArticlePubMedGoogle Scholar
- R Core Team: R: A Language and Environment for Statistical Computing. 2012, Vienna: R Foundation for Statistical ComputingGoogle Scholar
- Berry C, Hannenhalli S, Leipzig J, Bushman FD: Selection of target sites for mobile DNA integration in the human genome. PLoS Comput Biol. 2006, 2 (11): e157-10.1371/journal.pcbi.0020157. http://dx.doi.org/10.1371/journal.pcbi.0020157,PubMed CentralView ArticlePubMedGoogle Scholar
- Wang GP, Ciuffi A, Leipzig J, Berry CC, Bushman FD: HIV integration site selection: analysis by massively parallel pyrosequencing reveals association with epigenetic modifications. Genome Res. 2007, 17 (8): 1186-1194. 10.1101/gr.6286907. http://dx.doi.org/10.1101/gr.6286907,PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.