Temporal transcriptional response to latency reversing agents identifies specific factors regulating HIV-1 viral transcriptional switch
© Venkatachari et al. 2015
Received: 7 September 2015
Accepted: 25 September 2015
Published: 6 October 2015
Latent HIV-1 reservoirs are identified as one of the major challenges to achieve HIV-1 cure. Currently available strategies are associated with wide variability in outcomes both in patients and CD4+ T cell models. This underlines the critical need to develop innovative strategies to predict and recognize ways that could result in better reactivation and eventual elimination of latent HIV-1 reservoirs.
Results and discussion
In this study, we combined genome wide transcriptome datasets post activation with Systems Biology approach (Signaling and Dynamic Regulatory Events Miner, SDREM analyses) to reconstruct a dynamic signaling and regulatory network involved in reactivation mediated by specific activators using a latent cell line. This approach identified several critical regulators for each treatment, which were confirmed in follow-up validation studies using small molecule inhibitors. Results indicate that signaling pathways involving JNK and related factors as predicted by SDREM are essential for virus reactivation by suberoylanilide hydroxamic acid. ERK1/2 and NF-κB pathways have the foremost role in reactivation with prostratin and TNF-α, respectively. JAK-STAT pathway has a central role in HIV-1 transcription. Additional evaluation, using other latent J-Lat cell clones and primary T cell model, also confirmed that many of the cellular factors associated with latency reversing agents are similar, though minor differences are identified. JAK-STAT and NF-κB related pathways are critical for reversal of HIV-1 latency in primary resting T cells.
These results validate our combinatorial approach to predict the regulatory cellular factors and pathways responsible for HIV-1 reactivation in latent HIV-1 harboring cell line models. JAK-STAT have a role in reversal of latency in all the HIV-1 latency models tested, including primary CD4+ T cells, with additional cellular pathways such as NF-κB, JNK and ERK 1/2 that may have complementary role in reversal of HIV-1 latency.
KeywordsHIV-1 Latency Transcripts SDREM Gene expression Transcription factors
Development of combined antiretroviral therapy (cART) has made it possible to treat and control HIV-1 infection in infected individuals. Despite these major strides in the management of HIV-1 infection, HIV-1 cure still remains a challenge. Virus persists indefinitely even when the patient is treated with cART. Studies have identified resting CD4+ T cells as the main reservoir for the latent virus [1–5]. Discontinuation or interruption of treatment in well-controlled subjects results in increasing virus titer and recurrence of clinical symptoms associated with HIV/AIDS. This is due to the ability of the virus to establish a state of latent infection in a small number of cells where the virus remains transcriptionally silent for long periods of time [6–8]. Currently available antiretroviral agents are effective against actively replicating virus and prevent new infection; however, they fail to eliminate the latent reservoirs. Hence there is an urgent need to develop novel therapies that can eliminate the latent viral reservoirs. It is established that CD4+ resting memory T cells are the latent cells in vivo [9, 10], that are non permissive for viral gene expression. However, the mechanism(s) and the factors involved in this process are not fully defined. Though resting CD4+ T cells are the main reservoir for the latent virus, additional reservoirs include other cell types of T cell linage, including naïve T cells; macrophages; and cells of the central nervous systems [11–15]. Various methods have been proposed to activate and kill the latently infected cells; however, these strategies have resulted in variable outcomes in clinical trials, indicating that a better understanding of the factors and/or mechanisms regulating reactivation of latent virus is needed.
Studies have shown that viral latency is facilitated in part by transcription factors such as NF-κB, AP-1, NFAT and histone modification as well as DNA methylation in HIV-1 LTR in different HIV latency cell models and in primary cells [16–21]. However, compounds/reactivators targeting histone acetylation (SAHA, suberoylanilide hydroxamic acid), NF-κB activation (prostratin) or T cell activation (α-CD3/α-CD28) do not consistently reactivate the latently infected cells in patients, suggesting that HIV-1 latency is controlled by multiple host factors and regulatory proteins. Similarly the ability of latency reversing agents to reactivate latent HIV-1 in various in vitro T cell based HIV-1 latency models are also highly variable, suggesting that multiple factors and mechanisms are involved in different cell types [22–25].
In an effort to identify the factors and their regulatory pathways involved in reactivation mediated by a specific activator, we generated and integrated genome wide transcriptome datasets using robust Systems Biology approaches. As a “proof of concept”, ACH-2 cells, a T cell line harboring a single copy of latent HIV-1 provirus, was activated with SAHA, prostratin or TNF-α and transcriptome analyses was performed at multiple time points to identify host cellular transcripts and upstream transcription factors (TFs) that regulated activation. ACH-2 cells have a mutation in TAR region affecting the Tat mediated transactivation, this helps to identify the factors involved in early stages of HIV-1 latency reversal in the absence of viral proteins including Tat, though absence of effective Tat-TAR axis can promote latency status of infected cells. Our methods were used to reconstruct dynamic signaling and regulatory networks for these activators. The reconstructed model identified several key regulators for each treatment, and the specific set of genes they control either directly or indirectly. Furthermore, specific cellular factors and signaling pathways identified by SDREM analyses are validated using small molecule inhibitors. Our combinatorial approach was able to identify the key pathways involved in latent HIV-1 reactivation. JAK-STAT appears to have a critical role in HIV-1 transcription independent of the method of reactivation. Additionally, regulatory pathways involving JNK is required for virus reactivation by SAHA in ACH-2 cells, whereas, ERK1/2 and NF-κB, pathways are involved with reactivation with prostratin and TNF-α, respectively. Evaluation using another HIV-1 latent cell line, J-Lat as well as human primary resting T cells also confirmed that many of the factors associated with latency reversing agents are similar in different cell lines, though differences are also noted. Here we have demonstrated that temporal transcriptional analysis assessed by SDREM is a useful tool to identify the key regulators and/or pathways involved in latent HIV-1 reactivation.
Cell culture and reagents
ACH-2, A3.01 and J-Lat cell lines were obtained through the NIH AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH (ACH-2 cell line from Dr. Thomas Folks [26, 27]; J-Lat Full Length Clones from Dr. Eric Verdin [28, 29]) and were maintained in RPMI containing 10 % FBS, 1 % l-glutamine and 1 % penicillin–streptomycin (GIBCO). Multiple J-Lat cell lines derived from the same parent cell line, Jurkat T cell line, were included in our study. These cell lines contain either a full length GFP reporter virus (HIV-1 ΔNΔE-GFP) (3 cell clones—FL8.4, FL9.2 and FL10.6) or LTR–Tat-IRES-GFP (5 cell clones—TG82, TGA1, TGA2, TGA7 and TGH2). SAHA (suberoylanilide hydroxamic acid), prostratin, TNF-α and PHA-M (Phytohemagglutanin-M) were obtained from Sigma-Aldrich and R&D Systems. IΚK2 inhibitor V and SB203580 were obtained from CalBiochem. FK506 (Tacrolimus), Cyclosporin A and SB600125 were purchased from Abcam Biochemicals. Rottlerin, U0126 and WP 1066 were obtained from Sigma-Aldrich, Cell Signaling Technology and Enzo Life Sciences, respectively. In time kinetic experiments, the cells were treated with a compound and cells were collected at multiple time points from the single pool to reduce variability between time points within an experiment (paired design, repeated measures).
Purification of primary resting CD4+ T cells
Peripheral blood mononuclear cells were isolated from HIV-1 seronegative whole blood by Ficoll-Paque density gradient centrifugation (GE Healthcare). CD4+ T cells were purified by magnetic bead negative selection using a CD4+ T cell isolation kit (Miltenyi Biotec) and resting CD4+ T cells (rCD4+) were purified by magnetic bead negative selection using anti-CD25, anti-CD69 and anti-HLA-DR antibodies (Miltenyi Biotec).
Primary T cell treatment, infection and stimulation
Freshly isolated rCD4+ T cells were cultured at a density of 1–2 × 106 cells/mL in RPMI-1640 supplemented with 10 % FBS and pencillin/streptomycin/l-glutamine (100 Units/mL; 100 µg/mL; 0.292 mg/mL respectively, Life technologies) along with the chemokine CCL19 (100 nM, R&D systems), as described previously . After being in culture for 2 days, rCD4+ T cells were infected with HIV-1LAI (MOI = 1, titered on GHOST cells) for two to 3 h at 37 °C. The cells were then washed twice with media to remove any free virus and replaced back at 37 °C in fresh media. Two and four days post-infection, rIL-2 (10 U/mL, Roche) and EFV (300 nM, NIH AIDS Reagent Repository) were added to each well to inhibit any new rounds of infection that may occur. Seven days post-infection, the cells were washed again with media and plated in a 96-well plate at 100,000 cells/well. Cells were either stimulated with anti-CD3/CD28 (3 beads/cell, Life technologies) or 5 µM prostratin, and one of five inhibitors: 10 µM AZD 6244, 25 µM IΚK2-V, 50 µM SP600125, 5 µM WP1066 or 50 µM SB203580 were included at either 4 h prior to or 4 h post stimulation. Four days post-stimulation, half of the media was replaced with fresh media containing rIL-2 (10 U/mL, Roche) and EFV (300 nM).
Surface staining of the cells was performed with CD2, CD3, CD4, CD28 antibodies or isotype controls as described . For the detection of intracellular p24, fixation and permeabilization were carried out using the CytoFix–CytoPerm kit (BD Biosciences, Mountainview, CA, USA) and intracellular p24 staining was performed at room temperature for 1 h using 1 µl of anti-p24-FITC antibody (Coulter, Miami, FL, USA; clone KC47) per 106 cells, followed by two washes in Perm-Wash buffer, and finally resuspended in FACS buffer. Samples were analyzed using Fortessa (BD Biosciences) with 20,000-gated events acquired for each sample, and the results were analyzed using FlowJo software (Tree Star, Inc., OR, USA).
Real-time RT-PCR analysis
Total RNA extracted using RNeasy Mini Kit (Qiagen, CA, USA) was used to quantitate RPLPO, Multiply Spliced RNA (MS-RNA) and HIV-1 gag–pol transcripts by real-time PCR as described before [32, 33]. Briefly, a two-step RT-PCR was performed as follows: RNA was reverse transcribed using Taqman Reverse Transcription Reagents (Applied Biosystems, Foster City, CA, USA); Real-time PCR was carried out in triplicate using primer/probe sets specific for MS-RNA, HIV-1 gag-pol and ribosomal large protein (RPLPO, Applied Biosystems). The comparative CT method was used to determine the relative level of MS-RNA and HIV-1 gag transcripts by normalizing to the RPLPO control transcript.
For HIV-1 RNA quantification from primary cell culture, culture supernatant samples were spun at 16,100×g for 70 min to pellet HIV-1 virions. HIV-1 RNA was extracted from the virions using the RNeasy PLUS Mini Kit per the manufacturers protocol (Qiagen). To quantify total HIV-1 RNA in the culture supernatant, the extracted HIV-1 RNA samples were first converted into cDNA followed by real-time PCR using the protocols previously described  with few modification (AffinityScript Multiple Temperature RT (Agilent technologies) was used instead of Superscript II RT). The primers and probe used to quantify HIV-1 RNA were used as described previously . High copy number HIV-1 RNA transcripts were serially diluted to use as a RNA standard also as previously described .
Transcriptome profiling and data analysis
Illumina HT-12 V4 array bead chips (Illumina, Inc., San Diego, CA, USA) were used for whole genome transcriptome analysis for mRNA profiling after different treatment of ACH-2 cells. Each array targets about 47,231 probes that include 28,688 well-characterized or annotated coding transcripts along with 11,121 coding transcripts with provisional annotation and remaining being non-coding transcripts and splice variants. RNA samples (1 μg) were labeled using the ‘TotalPrep RNA’ labeling kit (Ambion), reverse transcribed to cDNA; cRNA was synthesized from cDNA with labeling and hybridized onto array bead chips overnight on rocker and scanned on ‘iScan system’, according to the manufacturer’s protocols as well as standardized protocols developed by the Genomics and Proteomics Core Laboratories at the University of Pittsburgh. Datasets will be deposited in NCBI gene expression and hybridization array data repository GEO database. The data were analyzed using GenomeStudio to identify the differentially regulated gene transcripts. The data were normalized by rank invariant method and no background subtraction was included, additionally, the missing samples were excluded. For calculating differential expression, the Illumina custom model was included along with multiple testing corrections using Benjamini and Hochberg False Discovery Rate, which is a standard methodology recommended by GenomeStudio to compare paired data . The differential score is a transformation of the p value that provides directionality to the p-value based on the difference between the average signal at time point zero versus different time points. The formula used for calculating Differential score = 10 × (Mean signal intensity at given time point (μt) − Mean Signal intensity at time point 0 (μt0)) × Log10p. A Differential score of ±13, corresponding to p < 0.05 was considered as the cut-off to identify significantly regulated transcripts.
Gene set enrichment analysis (GSEA)
To identify the biological process/function associated at virus replication at initial virus reactivation and later productive stage, the transcriptome data was analyzed using GSEA/MSigDB (version 4.0) (http://www.broadinstitute.org/gsea/msigdb/annotate.jsp) [37, 38]. First, a list of genes (regulated by more than twofolds, with p-value <0.05) was obtained for the time point in each treatment corresponding to virus reactivation and gag production/virus release (Multiple probes for the same gene was integrated together and analyzed at gene level). The identified genes were then analyzed using GSEA, with an FDR q-value below 0.05. This represents genes coordinately regulated in predefined gene sets from various biological pathways.
Signaling and dynamic regulatory events miner (SDREM)
To reconstruct signaling and regulatory networks activated following different treatments, we used SDREM as described [39, 40]. For the regulatory part, SDREM integrates condition specific time series gene expression data with global protein-DNA interaction data to identify bifurcation events in a time series (places where the expression of previously co-expressed set of genes diverges)–and the transcription factors (TFs) controlling these split events. While some TFs are transcriptionally activated, others are only activated post translationally via signaling networks. To identify and explain these TFs, the second part of SDREM links sources (host proteins that directly interact with the virus, in this case human proteins that were determined to physically interact with HIV-1 viral proteins ) to the TFs determined to regulate the regulatory network. This part of SDREM uses experimentally derived protein–protein interaction (PPI) from Biogrid and protein modification data  to infer such pathways—while imposing the constraint that the direction of PPI in the inferred pathways is consistent. These two parts (regulatory and signaling reconstruction) iterate a fixed number of times until the final network is obtained. See Gitter and Bar-Joseph  for complete details.
ACH-2 cells or A3.01 cells were treated with SAHA, prostratin, or TNF-α, and at the indicated time point, the cells were washed twice with PBS and lysed in RIPA buffer containing 50 mM Tris (pH 7.5), 150 mM NaCl, 1 % Triton X-100, 1 mM sodium orthovanadate, 10 mM sodium fluoride, 1.0 mM phenylmethylsulfonyl fluoride, 0.05 % deoxycholate, 10 % sodium dodecyl sulfate, aprotinin (0.07 trypsin inhibitor unit/ml), and the protease inhibitors leupeptin, chymostatin, and pepstatin (1 μg/ml; Sigma). Cell lysates were clarified by centrifugation, and total cell lysates (30 μg) were separated on a 12 % sodium dodecyl sulfate–polyacrylamide gel (SDS-PAGE) electrophoresis gel, transferred, and immunoblotted with anti-NT5C3 (santa cruz biotechnology) or anti- tubulin antibodies. The blots were developed using an ECL kit (Amersham Biosciences, Piscataway, NJ, USA).
Gag p24 ELISA
Following different treatment of ACH-2 and J-Lat cells, supernatants were collected and analyzed for the amount of HIV-1 p24 Gag. Gag (p24), was measured by using enzyme-linked immunosorbent assay (ELISA) kit (Zeptometrix, Buffalo, NY, USA) according to the manufacturer’s protocol.
The results from p24 ELISA, comparative reactivation studies, and inhibitor experiments are expressed as mean ± standard deviation. The data are analyzed using the Student’s t test for paired samples.
Multiple complementary factors are involved in reversal of latent HIV-1
To assess whether these activators could have an additive effect in reversal of latent virus, ACH-2 cells were treated with a combination of submaximal concentrations of these compounds to reactivate latent virus (Fig. 1c). Results indicate that combining SAHA and prostratin resulted in an increase in percentage of cells expressing p24 Gag. Treatment with SAHA or prostratin alone, showed 57–65 and 32–37 % of p24 positive cells, respectively, whereas a combination of these two activators resulted in 87–92 % of p24 positive cells. These results suggest an additive effect, implying that independent cellular factors are involved in SAHA and prostratin mediated latent HIV-1 reactivation and that these factors could complement each other to potentiate the effect. A similar increase in percentage of p24 Gag positive cells was observed when SAHA was combined with TNF-α, an increase from 23–28 to 83–87 % p24 positive cells. Similarly, when prostratin and TNF-α were combined together, the p24 positive cells increased to 46-62 %, a lesser extent than the additive effect, suggesting that prostratin and TNF-α may share few common cellular factors and pathways in virus reactivation. Combining SAHA, prostratin and TNF-α together resulted in maximum reactivation (93–97 % p24 positive cells), similar to that observed in PMA treated positive control. Collectively, these results suggest that HIV-1 reactivation from latency is mediated by multiple pathways involving different cellular factors that are complementary as well as common depending on the context of the activator used.
Kinetics of latent HIV-1 reactivation in ACH-2 cells indicates that transcriptional reactivation of viral RNA is an early event
To understand the time kinetics required for HIV-1 transcriptional reactivation in these cells, multiply spliced viral transcripts and unspliced gag RNA was assessed over time (Fig. 2b). Total RNA was extracted from part of the cells collected at the same time points to assess MS-RNA and unspliced Gag RNA (Fig. 2b, c) and virus released in the supernatant was measured by p24 ELISA (Fig. 2d), as described in materials and methods. Results indicate that prostratin initiated virus transcription within 2–4 h and that the peak expression of MS-RNA was reached at 6 h. Six hours post treatment there was a gradual decrease in the multiply spliced RNA and a corresponding increase in full-length gag transcripts occurred (Fig. 2b, c). The lesser fold change in cellular gag transcripts can be attributed to the packaging of full length gag transcripts in virus particles and their release in the supernatant, as measured by an increase in p24 virus Gag (Fig. 2d), resulting in sinusoidal curves of gag transcripts. In case of SAHA, we observed an increase in multiply spliced viral transcripts around 4–6 h, which was to a lesser extent than the fold change observed in prostratin. In case of gag transcripts, the increase was also to a lesser extent in comparison to prostratin occurring around 10 h. This correlates with the detection of p24 Gag by flow cytometry. The release of virus particles in SAHA as measured by an increase in p24 Gag protein in the supernatant was observed around 10 h. TNF-α resulted in an increase in multiply spliced transcripts by 6 h though the increase was very modest in multiple experiments; similarly, the increase in gag transcripts was observed around 10 h. Virus released in supernatant in TNF-α treatment of ACH-2 cells was similar to SAHA, where an increase in p24 was detected by 10 h and continued to increase over time. These results suggest that HIV-1 reactivation at transcription level in ACH-2 cells is an early event occurring within the initial 4–6 h post reactivation. Additionally, we noted that reactivation of latent virus in ACH-2 cells is a sequential process similar to acute HIV-1 infection in T cells, where the initial transcription initiation is followed by transcription of multiply spliced viral RNA transcripts, and a subsequent shift to full length gag transcripts and the final stage of viral RNA packaging and virus release.
Transcriptome analysis to identify specific factors involved in virus reactivation
Differentially regulated cellular transcripts in ACH-2 at indicated time points following treatment with SAHA, prostratin or TNF-α
A single copy of HIV-1 in ACH-2 is integrated in the intron between exons 5 and 6 in NT5C3A, variant 1 on chromosome 7 . Analysis of transcriptional changes in cellular genes around this region clearly identifies that NT5C3 as the only host gene that is consistently induced with the reactivation of the virus in ACH-2 cells to 7.4- 18.5 fold. None of the closely associated gene transcripts namely RP9, BBS9, FKBP9 or RP9P changes more than 1.4 to 0.9 fold with reversal of HIV-1 latency (Additional file 2: Table S1, Additional file 3: Fig. S2A). Time kinetics of changes in NT5C3 transcripts identifies that transcripts are induced as early as 2 h following treatment of ACH-2 cells with prostratin and TNF-α, and the levels continue to increase over time. With SAHA an increase in NT5C3 transcript is noticed at 4 h and this increase remains progressive and sustained, finally resulting in highest fold change of 18.5 in comparison to prostratin (12.6 fold change) and TNF-α (7.4 fold change) (Additional file 3: Fig. S2B). Induction of NT5C3 transcripts as detected by RNA hybridization in Illumina HT-12 chip was also detected at the protein level by western blot using NT5C3 specific antibody (Additional file 3: Fig. S2C). It is observed that the NT5C3 protein is not present in the media control, whereas detected upon treatment with prostratin, SAHA or TNF-α. Comparison of changes in protein level over time with the parent cell line, A3.01 cells, indicates that NT5C3 is expressed in A3.01 cells even in the absence of latency reactivating agents, Additional file 3: Fig. S2C, lane 5, suggests that, the integration of HIV-1 virus and its latent state in ACH-2 cells is associated with suppression of host cellular gene, NT5C3. With the reactivation of latent HIV-1 with either prostratin, SAHA or TNF-α is associated with reversal of this inhibition of NT5C3 expression.
Gene expression data sets were assessed to identify the pathways or major functions that are targeted by the reactivators using GSEA as described . Results presented in Fig. 3b focused on time points corresponding to the expression of viral transcripts and Gag synthesis/release indicate that the differentially expressed genes prior and/or during the time point represent primarily genes belong to transcription factors, protein kinases, cell differentiation markers and cancer related genes as well as cytokines and growth factors.
Transcriptional regulators and pathways involved in the reactivation of HIV-1
List of top 30 cellular factors involved in SAHA, prostratin or TNF-α induced latent HIV-1 reactivation in ACH-2 cells
Similar analyses were performed for prostratin induced reactivation in ACH-2 cells (Fig. 4c, d). Results indicate that TFs—HNF1A, POU2F1, SPI1, ETS1, LEF1, POU2F1, CEBP, SMAD, TCF7L2, RELA, GATA, STATs, BRCA1, ZEB1, FOXC1 and MZF1 are regulating the transcriptional changes observed at 4 h post prostratin treatment (Fig. 4c). It is also noted that the expression levels of ETS1 and LEF1 are increased and the expression of GATA3, STAT5A and CEBPB is reduced at this time point. As prostratin is an activator of Protein Kinase C (PKC), PKC was included as a source and the regulatory pathways and factors that were predicted to mediate the changes observed in the regulatory models are identified in Fig. 4d. The top 30 factors predicted by SDREM as critical for prostratin effects in ACH-2 cells are also presented in Table 2. Results indicate that MAPK1, SMAD2/4, SRC, MAX, GSK3B, and LEF1 are identified as unique factors involved in prostratin induced signaling and predict PRKCQ as the main driving factor responsible for the changes observed in cellular transcripts during early time duration of 0–4 h following prostratin treatment.
Unlike prostratin, TNF-α exhibited delayed reactivation kinetics in ACH-2 cells. SDREM analysis of time kinetic data obtained from TNF-α treated cells presented in Fig. 4e, suggests that TFs- VDR, RUNX1, SP1 and SP3 are contributing to the differential regulation of cellular transcripts observed at 2 h. TFs related to NF-κB—NF-κB1, NF-κB2, RELA; JAK-STAT related factors—STAT4, STAT5A, STAT6; CEBPA, CEBPE, GATA3, CUX1, ETV4, ETS1, FKBP4, CD40, ZEB1 are identified to be responsible for changes observed in the cellular transcripts at 2 and 4 h. The regulatory pathways and factors likely responsible for changes in TFs originating from TNFRSF1A and TNFRSF1B are indicated in Fig. 4f and Table 2. TFs, Jun, Myc, RelA, STAT1/3, TP53, NR3C1, CEBPB, and CREBBP, which are known to bind to HIV-1 LTR were present in all three treatments. Together these results suggest that specific transcription factors and their associated cellular signaling pathways are involved in HIV-1 reactivation in the ACH-2 cell line. It should be noted that SDREM identifies the upstream regulatory factors as the most probable factors responsible for changes observed in the transcripts at the identified time point and hence changes (activation/inhibition) in these upstream factors occur prior to the time points when the changes in transcripts were evaluated.
Validating the predicted cellular factors and TFs in latent HIV-1 reactivation in ACH-2 cells by specific inhibitors
In the case of prostratin, IKK2 inhibitor V, SP600125, U0126, AZD6244 and Rottlerin have a significant effect in preventing virus reactivation, suggesting that prostratin uses multiple alternate pathways involving NF-κB, JNK, ERK1/2 and PKC to induce HIV-1 reactivation (Fig. 4c, d). These results correlate with the SDREM prediction for prostratin suggesting a major role for ERK1/2. Similarly, SDREM analysis was accurate in identifying factors involved in TNF-α induced HIV-1 reactivation, which involves NF-κB, JNK and PKC, but not NFAT. p38 inhibitor did not increase prostratin or TNF-α mediated HIV-1 reactivation. The percentage of cells expressing p24 was significantly below the background expression observed with DMSO, suggesting that the JAK → STAT pathway is involved in regulating HIV-1 LTR activity. Irrespective of reactivation agents, STATs were also identified as a common factor in SDREM analysis for all these activating agents. Furthermore, it is important to note that no cellular toxicity was observed in these cells with these inhibitors, with the exception of WP1066 and Rottlerin, which showed 40–60 % reduction in cell viability (Fig. 5b); therefore only live cells were gated in flow cytometry based assessment of p24 positive cells. Additionally, dose dependent inhibitory effects can be observed with WP1066 and Rottlerin, and consistent inhibition was noted at lower concentration, when the associated cellular toxicity was minimal (Additional file 4: Fig. S3).
ACH-2 cells were incubated with the inhibitors in the absence of reactivating agents, to evaluate whether these inhibitors could potentially modulate HIV-1 transcription. Results indicate that only the p38 inhibitor–SB203580 and JAK-STAT inhibitor–WP1066 significantly altered the HIV-1 transcription (Fig. 5c). WP1066 reduces the basal level HIV-1 transcription in ACH-2 cells suggesting a role for STAT in HIV-1 transcription. SB203580 increased HIV-1 transcription from 15 to 42 % at 16 h post treatment, suggesting that p38 inhibitor SB203580 could directly reactivate latent HIV-1 in ACH-2 cells.
Validation of signaling pathways and regulatory networks in J-Lat latent cells
Next, we validated the specific TFs and regulatory networks in J-Lat based latency cell lines using inhibitors as described in our methods. Results indicate that pathways are similar in J-Lat and ACH-2 cells, with minor differences (Fig. 6b). JAK-STAT inhibitor, WP1066 and PKC inhibitor—Rottlerin have major effects in inhibiting reactivation of latent HIV-1 in J-Lat clones (close to background GFP level), suggesting that these signaling pathways are critical for HIV-1 LTR activity (Additional file 5: Fig. S4). ERK1/2 inhibitor U0126 and AZD6244 have a significant effect in inhibiting HIV-1 reactivation specifically in prostratin (~74–82 %), suggesting that ERK1/2 may be the major pathway involved downstream of PKC in prostratin mediated latent HIV-1 reactivation in J-Lat clones. Similarly NF-κB is the critical factor for TNF-α mediated reactivation in J-Lat cells (57–47 % inhibition with IΚK2 inhibitor V in J-Lat FL10.6 cells, 67–83 % inhibition in J-Lat TGA1 cells and 90–94 % inhibition in J-Lat TGA2 cell line). However, JNK exhibits a central role in activating HIV-1 LTR with SAHA, prostratin or TNF-α in J-Lat cells. Interestingly, it is also noted that IΚK2 inhibitor V inhibited SAHA mediated reactivation in J-Lat cells suggesting a role for NF-κB (~62–83 %). This correlates with the relatively reduced ability of SP600125 to inhibit the SAHA mediated reactivation in comparison to reactivation observed ACH-2 cells (60–70 % inhibition in J-Lat cells versus 85–95 % inhibition in ACH-2 cells). This suggests that when JNK mediated signaling is blocked in J-Lat, NF-κB related pathways could still lead to reactivation of latent HIV-1.
Time dependent inhibition of signaling pathways and regulatory networks in J-Lat cells identify unique role for JAK- STAT early in latent virus reactivation
With ERK1/2 inhibitor, we observed a more gradual loss of inhibition over time (Fig. 7b, d, f, h) when prostratin and PHA-M are included to reactivated the virus. The inhibition is >95 % when AZD 6244 was included 4 h prior to activation, this inhibition gradually decreased to 40–60 % at 8 h post stimulation. JNK inhibitor, SP600125, also demonstrated consistent inhibition around 40–60 % with both the cell lines with all the tested reactivating agents but for prostratin in J-Lat TGA1cells (Fig. 7f), when included 4 h prior to stimulation. This inhibitory effect reduced progressively with time over to 40 % especially in J-Lat cell line FL10.6 (Fig. 7a–d) but did not change drastically till 8 h in J-Lat cell line TGA1. Rottlerin showed consistent inhibition (~80–90 %) over time till 8 h in both the cell lines when stimulated with SAHA, TNF-α and PHA-M, but with prostratin, the inhibitory effect of rottlerin progressively reduces to ~20 % at 8 h post stimulation (Fig. 7b, f). NF-κB inhibitor was able to consistently inhibit virus reactivation (~60–90 %) that remains stable over time in both the cells with all the tested reactivators but for in J-Lat cells TGA1 with prostratin (Fig. 7f). These results suggest that JAK-STAT has a role in early stage of reactivation, in the initial 2–4 h, and inhibiting JAK-STAT pathway after this interval does not inhibit latent virus reactivation with all the tested activators, though with delayed kinetics for TNF-α and PHA-M. And this early role of JAK-STAT pathway is essential for viral reactivation which cannot be complemented by other signal or factors. With other regulatory factors like NF-κB, PKC, JNK and ERK1/2, the activation of these factors are essential at all-time points for effective reactivation, though at later time points, the inhibition of these activated factors can complement each other, which leads to gradual loss of inhibitory activity.
Role of signaling pathways and regulatory networks in latent virus reactivation in primary T cells
Inability to consistently reactivate the latent viral reservoir is presently the main impediment in the “shock and kill” strategy to achieve HIV-1 cure. Numerous therapeutic agents and methods have been proposed to reactivate the latent virus but the results vary widely between individual patients, suggesting that multiple factors are involved in determining the reactivation of latent HIV-1. This variability in outcomes is true for in vitro primary CD4+ T cell based models as well as established T cell lines, where different activators are effective to reactivate latent virus to variable extent. Hence there is an urgent need to develop innovative approaches to predict and identify strategies that could result in more effective reactivation and eventual viral elimination. Here in this study, using ACH-2, a T cell line (that has a single latent copy of replication competent provirus), we describe the application of a robust Systems Biology approach to predict factors that are potentially involved in latent HIV-1 reactivation. Changes in transcripts of NT5C3 gene in ACH-2 cells suggest that onset of HIV latency in these cell lines is associated with silencing of the host gene where the virus integrated and can be induced upon virus reactivation, suggesting methods by which a retroviral genetic element can regulate host cellular expression. Additionally these changes indicate that HIV latency and reactivation in ACH-2 cells is associated with local chromosomal changes that are specific to the region of virus integration, which could be due changes in the nucleosomal assembly and/or DNA methylation.
Using Systems Biology method, we integrate several genome wide transcriptome datasets to reconstruct a dynamic signaling and regulatory network for each reactivating agent. This helped to specifically identify the critical factors for individual treatments along with their time of activation and the specific set of genes regulated by these critical factors. We validated our approach using specific inhibitors to confirm the cellular factors that are critical for latent HIV-1 reactivation. The latent HIV-1 virus in ACH-2 cells is integrated in chromosome 7 and has a point mutation in TAR- Tat axis C37T. It has been previously identified that the TAR mutation renders the virus less responsive to Tat and suggests that Tat-TAR axis may have a role in HIV-1 latency in patients . The mutation in TAR region helps us to understand the cellular pathways involved in HIV-1 latency reversal in the early phase in the absence of viral protein Tat. Additional reports also suggest a role for DNA methylation in LTR region of HIV-1 virus identified in ACH-2 as the major contributing factor for latency and reactivation from latency . Therefore, results from ACH-2 cells were further confirmed in J-Lat cell clones and primary T cell HIV latency model that do not contain any mutations and similar results were observed, thus confirming the validity of the cell lines as latent cell model. SDREM reconstructs accurate models of the cross talk between signaling pathways and transcriptional regulatory networks within cells, which are essential to understand complex host cellular response programs . This computational method combines condition-specific time series expression data with general protein interaction data to reconstruct dynamic and causal networks [50, 51]. Similar analyses were performed in H1N1 influenza infection to identify strain specific targets of Influenza infection .
Analyses of time kinetic transcriptome data corresponding to the initial period of reactivation extending up to 6 h in ACH-2 cells post treatment with SAHA, clearly identified Jun as the key regulatory factor that is responsible for the increase in cellular transcripts observed in the initial 2 h. SDREM also predicts NFATC1-4, CEBPE, TBP and ATF2 as potential factors modulating cellular transcript levels. These factors are known to bind to HIV-1 LTR and have a role in HIV-1 transcription. However, interpretation of predicted factors associated with a decrease in cellular transcript level is complex, since it can mean that either association or dissociation of these predicted TFs may lead to suppression of cellular transcripts). Some of the other top 30 factors regulated by SAHA include—STATs, TP53, RB1, NF-κB related factors, BRCA1, NR3C1, CTNNB1, SMADs and Myc (Table 2). Though these factors are part of the predicted regulatory network, validating the role of these factors is difficult, as targeting specific cellular factors can result in cells using alternative compensatory networks, and multiple TFs can be regulated by a single upstream regulatory pathway, hence the recommended approach includes targeting specific signaling or regulatory pathways rather than TFs. The SDREM prediction for SAHA highly correlates with the results in validation experiments, where Jun is identified as the critical cellular factor and validation experiments using JNK inhibitor SP600125 inhibits SAHA induced reactivation of HIV-1 in ACH-2 cells, and J-Lat cell lines supporting the validity of these analyses.
Multiple factors, including ETS1, CEBP, STATs, RELA and GATA 2/3, which are known to bind to HIV-1 LTR, are predicted by SDREM as active in ACH-2 cells within the initial 4 h following prostratin treatment. The predicted regulators that are unique for prostratin include MAPK1, SMADs2/4 and others. This prediction correlates well with the results obtained in validation experiments, where we observed specific inhibition of prostratin induced HIV-1 reactivation in the presence of ERK1/2, NF-κB pathway and JNK inhibitors. The SDREM prediction is also accurate for TNF-α treatment, which identified NF-κB related factors, STATs, CEBP, ETS1 and GATA3, and the validation results correlated with the inhibitory effect observed with IΚK-2 inhibitor V and SP600125. Additionally, it can be noted that STATs are common in the network for all the three treatments and correlate with the inhibitory effect of WP1066, though STAT may also have a role in regulation of basal HIV-1 LTR transcription, as we observed a decrease in p24 positive cells with WP1066 in the absence of latency reactivating agents. Similar results describing a central role for STAT were also reported in acute HIV-1 infection and in reactivation of latent HIV infection by using STAT inhibitors Ruxolitinib and Tofacitinib . Our time kinetics data in J-Lat cell lines further refines the role of STAT and restricts the requirement of JAK-STAT pathway to the initial immediate early stage of latent HIV reactivation. Other factors that are predicted in common for all three treatments include BRCA1, CEBPB, JUN, MYC, NR3C1, RELA, SMAD3, TP53 and UBC. Previous results suggest that Myc may have a role proviral latency by recruiting Histone Deacetylase 1 to the HIV-1 promoter , and we also found that c-Myc levels are reduced in the three treatments tested in ACH-2 cells, which may have a role in HIV-1 reactivation. The role of other factors identified as common or unique needs to be evaluated in patients.
The minor differences in degree of response and the cellular pathway involved are expected in different cell types. Previously it has been demonstrated that the latent HIV-1 virus in the human monocytic cell line U1 can be induced by cytokines, including TNF-α, which is sensitive to U0126 and has been suggested to involve cooperative interaction of AP-1 and NF-κB at the HIV-1 LTR . Our results though suggests that AP-1 and NF-κB have a critical role in TNF-α induced reactivation of latent virus in ACH-2 and J-Lat cells, but this does not require MAPK as U0126 was unable to inhibit the TNF-α mediated reactivation of latent HIV-1. This could be explained as the availability of alternate signaling pathways in different cell types, which might compensate for missing factors. Also, in the primary T cell model developed by Bosque et al. , it was shown that Cyclosporin A and p38 MAPK inhibitor inhibited αCD3/αCD28 or PHA-M mediated reactivation of latent HIV-1, but in ACH-2 cells, NFAT and p38-MAPK did not seem to have a role in latent HIV-1 reactivation, though NFAT seems to have a minor role in PHA-M induced latent virus reactivation in J-Lat clones. In our primary resting T cell model, we see a major role for JAK-STAT and NFκB in reversal of latency mediated by αCD3/αCD28 or prostratin, with variable role for JNK and ERK1/2 which were found to have a minor role in some donors but not all and the inhibition was not statistically significant. SB203580, p38 MAPK inhibitor, did not show any inhibition in our primary T cell model, suggesting no role for p38 in latent virus reactivation. Similar variations were observed when reactivity of different in vitro primary T cell models and cell lines were compared . SAHA was not tested in our primary cell model, as it is ineffective and does not reactivate latent HIV-1 as has been observed in in vitro studies using resting CD4+ T cells isolated from aviremic infected individuals [55–57]. Though ACH-2 cells have good dynamic range in their response to LRAs (Latency Reversing Agents), it can also be noted that the background of non-latent virus containing cells in ACH-2 population is around ~13–20 %. This background could potentially undermine the magnitude of fold change observed in cellular transcripts induced by LRAs at the population level. A twofold change in transcripts in the 75–80 % of cells will be observed as 1.5–1.6 fold change at the population. However, SDREM considers the changes observed in global transcripts over time that are independent of the magnitude of fold change to predict the upstream regulating factors, hence the reduced fold changes in transcripts does not affect the analyses.
Surprisingly, this temporal requirement of JAK-STAT was not observed at 4 h with prostratin or αCD3/αCD28 in primary T cell model suggesting that the absence could be due to delayed kinetics of the JAK-STAT dependent phase of virus reactivation. This is true in J-Lat cells when they were reactivated with TNF-α and PHA-M, where JAK-STAT dependent phase extended till 6 h post activation. Furthermore, within an individual, different latent reservoirs might respond differently based on the availability of cellular signaling pathways and other factors influencing the signal transduction. Given the short duration time kinetics data required and the small amount of cells required to obtain transcriptome data, SDREM has potential to be adapted to characterize cellular signaling pathways—to identify defects in pathways of HIV-1 subjects, and to identify alternate pathways that can bypass the defects to reactivate the HIV-1 transcription with minimal adverse effects. Also, the SDREM approach will help in assay development—to compare functional signaling pathways in cell models that can closely simulate conditions present in cells obtained from latent reservoir of HIV-1 patients.
Our results support the wide variations observed in outcomes both in patients and CD4+ T cell models, suggesting that cells derived from HIV-1 infected patients may respond differently to various activators and that a single strategy may not be optimal for all individuals, though JAK-STAT and NF-κB seems to be the common regulatory factors essential for latent HIV reactivation. JAK-STAT has an essential role in immediate early phase of HIV-1 latency reversal, where activation of STAT is critical for reversal of HIV-1 latency, though the duration of this immediate early phase varies between latent HIV-1 cell models tested and the LRAs used. Other cellular signaling pathways including NFκB, NFAT, JNK, and ERK1/2 are also involved in latent HIV-1 reversal and could play a complement role with each other that is essential during reactivation. SDREM based prediction of cellular signaling pathways will help to optimize the management of HIV-1 patients to achieve the best outcome and cure.
human immunodeficiency virus type 1
latency reversing agents
combined antiretroviral therapy
multiply spliced RNA
signaling and dynamic regulatory events miner
suberoylanilide hydroxamic acid
tumor necrosis factor-α
protein kinase C
phorbol 12-myristate 13-acetate
ribosomal large protein
histone acetyl transferase
heterochromatin protein 1
- RNA POLII:
RNA polymerase II
Conceived and designed the experiments: VA, NJV, ZBJ, NSC, JMZ, SJ. Performed the experiments: NJV, JMZ. Analyzed the data: VA, NJV, ZBJ, NSC, JMZ, SJ, AEM, AC. Wrote the paper: VA, NJV, ZBJ, NSC, JMZ, SJ. All authors read and approved the final manuscript.
This work was in part supported by R01 MH087247 to VA from NIMH, NIH and U01 AI-35041 based ARRA supplements to VA from NIAID, NIH. NSC was supported by P50GM082251.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Siliciano RF, Greene WC. HIV latency. Cold Spring Harb Perspect Med. 2011;1:a007096.PubMed CentralView ArticlePubMedGoogle Scholar
- Finzi D, Blankson J, Siliciano JD, Margolick JB, Chadwick K, Pierson T, Smith K, Lisziewicz J, Lori F, Flexner C, et al. 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:512–7.View ArticlePubMedGoogle Scholar
- Marcello A. Latency: the hidden HIV-1 challenge. Retrovirology. 2006;3:7.PubMed CentralView ArticlePubMedGoogle Scholar
- Chan CN, Dietrich I, Hosie MJ, Willett BJ. Recent developments in human immunodeficiency virus-1 latency research. J Gen Virol. 2013;94:917–32.PubMed CentralView ArticlePubMedGoogle Scholar
- Persaud D, Pierson T, Ruff C, Finzi D, Chadwick KR, Margolick JB, Ruff A, Hutton N, Ray S, Siliciano RF. A stable latent reservoir for HIV-1 in resting CD4(+) T lymphocytes in infected children. J Clin Invest. 2000;105:995–1003.PubMed CentralView ArticlePubMedGoogle Scholar
- Choudhary SK, Margolis DM. Curing HIV: pharmacologic approaches to target HIV-1 latency. Annu Rev Pharmacol Toxicol. 2011;51:397–418.PubMed CentralView ArticlePubMedGoogle Scholar
- Chun TW, Fauci AS. HIV reservoirs: pathogenesis and obstacles to viral eradication and cure. AIDS. 2012;26:1261–8.View ArticlePubMedGoogle Scholar
- Chun TW, Fauci AS. Latent reservoirs of HIV: obstacles to the eradication of virus. Proc Natl Acad Sci USA. 1999;96:10958–61.PubMed CentralView ArticlePubMedGoogle Scholar
- Siliciano JD, Siliciano RF. The latent reservoir for HIV-1 in resting CD4+ T cells: a barrier to cure. Curr Opin HIV AIDS. 2006;1:121–8.PubMedGoogle 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:6122–33.PubMed CentralView ArticlePubMedGoogle Scholar
- Chun TW, Carruth L, Finzi D, Shen X, DiGiuseppe JA, Taylor H, Hermankova M, Chadwick K, Margolick J, Quinn TC, et al. Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature. 1997;387:183–8.View ArticlePubMedGoogle Scholar
- Thompson KA, Cherry CL, Bell JE, McLean CA. Brain cell reservoirs of latent virus in presymptomatic HIV-infected individuals. Am J Pathol. 2011;179:1623–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Wong JK, Hezareh M, Gunthard HF, Havlir DV, Ignacio CC, Spina CA, Richman DD. Recovery of replication-competent HIV despite prolonged suppression of plasma viremia. Science. 1997;278:1291–5.View ArticlePubMedGoogle Scholar
- 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:1284–90.View ArticlePubMedGoogle Scholar
- Lambotte O, Taoufik Y, de Goer MG, Wallon C, Goujard C, Delfraissy JF. Detection of infectious HIV in circulating monocytes from patients on prolonged highly active antiretroviral therapy. J Acquir Immune Defic Syndr. 2000;23:114–9.View ArticlePubMedGoogle Scholar
- Matalon S, Rasmussen TA, Dinarello CA. Histone deacetylase inhibitors for purging HIV-1 from the latent reservoir. Mol Med. 2011;17:466–72.PubMed CentralView ArticlePubMedGoogle Scholar
- Chan JK, Greene WC. NF-kappaB/Rel: agonist and antagonist roles in HIV-1 latency. Curr Opin HIV AIDS. 2011;6:12–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Williams SA, Chen LF, Kwon H, Fenard D, Bisgrove D, Verdin E, Greene WC. Prostratin antagonizes HIV latency by activating NF-kappaB. J Biol Chem. 2004;279:42008–17.View ArticlePubMedGoogle Scholar
- Jiang G, Espeseth A, Hazuda DJ, Margolis DM. c-Myc and Sp1 contribute to proviral latency by recruiting histone deacetylase 1 to the human immunodeficiency virus type 1 promoter. J Virol. 2007;81:10914–23.PubMed CentralView ArticlePubMedGoogle Scholar
- Marban C, Suzanne S, Dequiedt F, de Walque S, Redel L, Van Lint C, Aunis D, Rohr O. Recruitment of chromatin-modifying enzymes by CTIP2 promotes HIV-1 transcriptional silencing. EMBO J. 2007;26:412–23.PubMed CentralView ArticlePubMedGoogle Scholar
- Blazkova J, Murray D, Justement JS, Funk EK, Nelson AK, Moir S, Chun TW, Fauci AS. Paucity of HIV DNA methylation in latently infected, resting CD4+ T cells from infected individuals receiving antiretroviral therapy. J Virol. 2012;86:5390–2.PubMed CentralView ArticlePubMedGoogle Scholar
- Spina CA, Anderson J, Archin NM, Bosque A, Chan J, Famiglietti M, Greene WC, Kashuba A, Lewin SR, Margolis DM, et al. An in-depth comparison of latent HIV-1 reactivation in multiple cell model systems and resting CD4+ T cells from aviremic patients. PLoS Pathog. 2013;9:e1003834.PubMed CentralView ArticlePubMedGoogle Scholar
- Ruelas DS, Greene WC. An integrated overview of HIV-1 latency. Cell. 2013;155:519–29.PubMed CentralView ArticlePubMedGoogle Scholar
- Krishnan V, Zeichner SL. Host cell gene expression during human immunodeficiency virus type 1 latency and reactivation and effects of targeting genes that are differentially expressed in viral latency. J Virol. 2004;78:9458–73.PubMed CentralView ArticlePubMedGoogle Scholar
- Olivares I, Ballester A, Lombardia L, Dominguez O, Lopez-Galindez C. Human immunodeficiency virus type 1 chronic infection is associated with different gene expression in MT-4, H9 and U937 cell lines. Virus Res. 2009;139:22–31.View ArticlePubMedGoogle Scholar
- Clouse KA, Powell D, Washington I, Poli G, Strebel K, Farrar W, Barstad P, Kovacs J, Fauci AS, Folks TM. Monokine regulation of human immunodeficiency virus-1 expression in a chronically infected human T cell clone. J Immunol. 1989;142:431–8.PubMedGoogle Scholar
- Folks TM, Clouse KA, Justement J, Rabson A, Duh E, Kehrl JH, Fauci AS. Tumor necrosis factor alpha induces expression of human immunodeficiency virus in a chronically infected T-cell clone. Proc Natl Acad Sci USA. 1989;86:2365–8.PubMed CentralView 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:1726–38.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:1868–77.PubMed CentralView ArticlePubMedGoogle Scholar
- Saleh S, Solomon A, Wightman F, Xhilaga M, Cameron PU, Lewin SR. CCR7 ligands CCL19 and CCL21 increase permissiveness of resting memory CD4+ T cells to HIV-1 infection: a novel model of HIV-1 latency. Blood. 2007;110:4161–4.View ArticlePubMedGoogle Scholar
- Venkatachari NJ, Majumder B, Ayyavoo V. Human immunodeficiency virus (HIV) type 1 Vpr induces differential regulation of T cell costimulatory molecules: direct effect of Vpr on T cell activation and immune function. Virology. 2007;358:347–56.View ArticlePubMedGoogle Scholar
- Saleh S, Wightman F, Ramanayake S, Alexander M, Kumar N, Khoury G, Pereira C, Purcell D, Cameron PU, Lewin SR. Expression and reactivation of HIV in a chemokine induced model of HIV latency in primary resting CD4+ T cells. Retrovirology. 2011;8:80.PubMed CentralView ArticlePubMedGoogle Scholar
- Venkatachari NJ, Alber S, Watkins SC, Ayyavoo V. HIV-1 infection of DC: evidence for the acquisition of virus particles from infected T cells by antigen uptake mechanism. PLoS One. 2009;4:e7470.PubMed CentralView ArticlePubMedGoogle Scholar
- Palmer S, Wiegand AP, Maldarelli F, Bazmi H, Mican JM, Polis M, Dewar RL, Planta A, Liu S, Metcalf JA, et al. New real-time reverse transcriptase-initiated PCR assay with single-copy sensitivity for human immunodeficiency virus type 1 RNA in plasma. J Clin Microbiol. 2003;41:4531–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Cillo AR, Krishnan A, Mitsuyasu RT, McMahon DK, Li S, Rossi JJ, Zaia JA, Mellors JW. Plasma viremia and cellular HIV-1 DNA persist despite autologous hematopoietic stem cell transplantation for HIV-related lymphoma. J Acquir Immune Defic Syndr. 2013;63:438–41.PubMed CentralView ArticlePubMedGoogle Scholar
- Illumina: GenomeStudio Gene Expression Module v1.0 User Guide (11319121). 2015:103.
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50.PubMed CentralView ArticlePubMedGoogle Scholar
- Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34:267–73.View ArticlePubMedGoogle Scholar
- Gitter A, Bar-Joseph Z. Identifying proteins controlling key disease signaling pathways. Bioinformatics. 2013;29:i227–36.PubMed CentralView ArticlePubMedGoogle Scholar
- Schulz MH, Pandit KV, Lino Cardenas CL, Ambalavanan N, Kaminski N, Bar-Joseph Z. Reconstructing dynamic microRNA-regulated interaction networks. Proc Natl Acad Sci USA. 2013;110:15686–91.PubMed CentralView ArticlePubMedGoogle Scholar
- Navratil V, de Chassey B, Meyniel L, Delmotte S, Gautier C, Andre P, Lotteau V, Rabourdin-Combe C. VirHostNet: a knowledge base for the management and the analysis of proteome-wide virus-host interaction networks. Nucleic Acids Res. 2009;37:D661–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Keshava Prasad TS, Goel R, Kandasamy K, Keerthikumar S, Kumar S, Mathivanan S, Telikicherla D, Raju R, Shafreen B, Venugopal A, et al. Human Protein Reference Database–2009 update. Nucleic Acids Res. 2009;37:D767–72.PubMed CentralView ArticlePubMedGoogle Scholar
- Pace MJ, Agosto L, Graf EH, O’Doherty U. HIV reservoirs and latency models. Virology. 2011;411:344–54.PubMed CentralView ArticlePubMedGoogle Scholar
- Hakre S, Chavez L, Shirakawa K, Verdin E. HIV latency: experimental systems and molecular models. FEMS Microbiol Rev. 2012;36:706–16.PubMed CentralView ArticlePubMedGoogle Scholar
- Finzi D, Hermankova M, Pierson T, Carruth LM, Buck C, Chaisson RE, Quinn TC, Chadwick K, Margolick J, Brookmeyer R, et al. Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. Science. 1997;278:1295–300.View ArticlePubMedGoogle Scholar
- Chun TW, Davey RT Jr, Ostrowski M, Shawn Justement J, Engel D, Mullins JI, Fauci AS. Relationship between pre-existing viral reservoirs and the re-emergence of plasma viremia after discontinuation of highly active anti-retroviral therapy. Nat Med. 2000;6:757–61.View ArticlePubMedGoogle Scholar
- Ishida T, Hamano A, Koiwa T, Watanabe T. 5′ long terminal repeat (LTR)-selective methylation of latently infected HIV-1 provirus that is demethylated by reactivation signals. Retrovirology. 2006;3:69.PubMed CentralView ArticlePubMedGoogle Scholar
- Emiliani S, Van Lint C, Fischle W, Paras P Jr, Ott M, Brady J, Verdin E. A point mutation in the HIV-1 Tat responsive element is associated with postintegration latency. Proc Natl Acad Sci USA. 1996;93:6377–81.PubMed CentralView ArticlePubMedGoogle Scholar
- Gitter A, Carmi M, Barkai N, Bar-Joseph Z. Linking the signaling cascades and dynamic regulatory networks controlling stress responses. Genome Res. 2013;23:365–76.PubMed CentralView ArticlePubMedGoogle Scholar
- Schulz MH, Devanny WE, Gitter A, Zhong S, Ernst J, Bar-Joseph Z. DREM 2.0: Improved reconstruction of dynamic regulatory networks from time-series expression data. BMC Syst Biol. 2012;6:104.PubMed CentralView ArticlePubMedGoogle Scholar
- Bar-Joseph Z, Gitter A, Simon I. Studying and modelling dynamic biological processes using time-series gene expression data. Nat Rev Genet. 2012;13:552–64.View ArticlePubMedGoogle Scholar
- Gavegnano C, Detorio M, Montero C, Bosque A, Planelles V, Schinazi RF. Ruxolitinib and tofacitinib are potent and selective inhibitors of HIV-1 replication and virus reactivation in vitro. Antimicrob Agents Chemother. 2014;58:1977–86.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang X, Chen Y, Gabuzda D. ERK MAP kinase links cytokine signals to activation of latent HIV-1 infection by stimulating a cooperative interaction of AP-1 and NF-kappaB. J Biol Chem. 1999;274:27981–8.View ArticlePubMedGoogle Scholar
- Bosque A, Planelles V. Induction of HIV-1 latency and reactivation in primary memory CD4 + T cells. Blood. 2009;113:58–65.PubMed CentralView ArticlePubMedGoogle Scholar
- Cillo AR, Sobolewski MD, Bosch RJ, Fyne E, Piatak M Jr, Coffin JM, Mellors JW. Quantification of HIV-1 latency reversal in resting CD4+ T cells from patients on suppressive antiretroviral therapy. Proc Natl Acad Sci USA. 2014;111:7078–83.PubMed CentralView ArticlePubMedGoogle Scholar
- Bullen CK, Laird GM, Durand CM, Siliciano JD, Siliciano RF. New ex vivo approaches distinguish effective and ineffective single agents for reversing HIV-1 latency in vivo. Nat Med. 2014;20:425–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Blazkova J, Chun TW, Belay BW, Murray D, Justement JS, Funk EK, Nelson A, Hallahan CW, Moir S, Wender PA, Fauci AS. Effect of histone deacetylase inhibitors on HIV production in latently infected, resting CD4(+) T cells from infected individuals receiving effective antiretroviral therapy. J Infect Dis. 2012;206:765–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Grivennikov SI, Karin M. Dangerous liaisons: STAT3 and NF-kappaB collaboration and crosstalk in cancer. Cytokine Growth Factor Rev. 2010;21:11–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Crotti A, Chiara GD, Ghezzi S, Lupo R, Jeeninga RE, Liboi E, Lievens PM, Vicenzi E, Bovolenta C, Berkhout B, Poli G. Heterogeneity of signal transducer and activator of transcription binding sites in the long-terminal repeats of distinct HIV-1 subtypes. Open Virol J. 2007;1:26–32.PubMed CentralView ArticlePubMedGoogle Scholar