Previous Article | Next Article ![]()
Journal of Virology, September 2004, p. 9458-9473, Vol. 78, No. 17
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.17.9458-9473.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
HIV and AIDS Malignancy Branch, National Cancer Institute, National Institutes of Health, Bethesda, Maryland
Received 24 February 2004/ Accepted 20 May 2004
|
|
|---|
|
|
|---|
Several studies have described changes in certain selected cellular genes due to the expression of HIV proteins (20, 60, 61). Other studies have described changes in cellular gene expression due to acute infection following the addition of infectious virus to cells (16, 29, 67), and additional studies have described changes in host cell gene expression that accompany infections by other viruses, offering interesting insights into viral pathogenesis (12, 33, 37, 51, 63, 72). While these studies provide insights into the extensive effects of specific viral proteins or acute HIV infection (reviewed in reference 4), the global gene expression changes that may accompany reactivation of cells latently infected with HIV and completion of the lytic phase of viral replication have not been described.
Current antiretroviral agents can control HIV replication but cannot eliminate latently infected cells. Models incorporating higher-order decay kinetics (30, 31, 53) suggest that therapy must last some 60 years to eradicate the virus (23, 49). The predictions of these models may reflect the slow loss of virus from "reservoirs." Latent HIV infection can exist in many such reservoirs, such as macrophages and resting memory CD4+ T cells (14, 46), which appear to "archive" HIV species (55).
Several lines of evidence show that treating host cells carrying HIV provirus with different agents can end latency (24, 26, 40, 47, 66, 70). For example, HIV-infected cells can be activated in vitro by inducing agents such as phorbol myristyl acetate (PMA), which can induce lytic replication in chronically infected cell lines. However, many agents that affect latent HIV infection are toxic or may not eliminate latent reservoirs in patients (13). Since viral reactivation is necessary for targeting by antiviral drugs (9), identification of new approaches that eject HIV from latency when used together with effective antiviral therapy might lead to the reduction of latent HIV reservoirs in an infected patient.
We studied the pattern of host cell gene expression when cells latently infected with HIV were induced to undergo the rest of the viral replication cycle. We found that cellular gene expression patterns are altered in latently infected cells even before they are induced to undergo lytic infection. Following the emergence of HIV from the latent state and during the subsequent completion of the lytic infection cycle, the host cell exhibited carefully ordered changes in the expression of a subset of its genes, which shadowed the well-known, ordered changes in the pattern of viral gene expression characteristic of the HIV replication cycle. Reasoning that the altered expression of some of the host cell genes prior to induction of lytic replication might be related to the establishment and maintenance of latency in these cells, we treated the chronically infected cells with agents targeting the gene products of certain differentially expressed genes and found that some of these agents could activate lytic HIV replication. The differentially expressed cellular genes may thus offer novel cellular targets for agents that activate latent virus into a lytic replication cycle and may provide potential new therapeutic approaches to eliminate latent viral reservoirs.
|
|
|---|
Cells were induced by addition of 20 ng of phorbol myristyl acetate (PMA; Sigma, St. Louis, Mo.) per ml for 1 h, after which the cells were washed with phosphate-buffered saline. Cells were harvested by centrifugation at 1,000 rpm for 10 min at 0.5, 3, 6, 8, 12, 18, 24, 48, 72, and 96 h after induction. HIV-infected and uninfected cells maintained and harvested in parallel with the PMA-treated cells but not induced with PMA were also harvested. For the time course experiment, 3'-azido-3'-deoxythymidine (AZT; Sigma) was not added to the ACH-2 or A3.01 cells in order to keep conditions as close to those of an acute infection as allowed by the experimental model. Harvested cells were washed thrice with ice-cold phosphate-buffered saline to remove the medium; cell pellets were snap frozen with an ethanol-dry ice mixture and stored at 80°C for subsequent RNA extraction. Three independent time course experiments (biological replicates) were performed.
In order to compare the expression profiles of chronically infected cell lines, ACH-2, U1, and J1.1 and their uninfected parental cell lines, A3.01, U-937, and Jurkat cells, respectively, were grown under identical conditions but in the presence of AZT (250 nM) in growth medium, and the cells were harvested as described above. In these studies, no inducing agent was used in either chronically infected or uninfected parental cell lines so as to study changes in cellular gene expression cells latently infected with HIV and uninfected cells.
Flow cytometry. To confirm viral replication following PMA induction, we measured the accumulation of intracellular p24 over a period of 48 h by measuring cell populations labeled with anti-p24 fluorescein isothiocyanate-labeled antibody by flow cytometry. Cells (ACH-2 and A3.01) were washed twice with ice-cold PBS and suspended in 50 µl of ice-cold permeabilization buffer (BD Biosciences, San Jose, Calif.), and incubated at 4°C in the dark for 30 min. The cells were fixed with the CytoFix/CytoPerm kit (BD Biosciences), and 5 µl of fluorescein isothiocyanate-labeled p24 antibody KC57 (Beckman Coulter) was added to detect intracellular p24. A3.01 samples labeled with fluorescein isothiocyanate-labeled p24 antibody served as controls to ensure that the parental cell line did not show any p24 accumulation over that in samples labeled with fluorescein isothiocyanate-labeled mouse immunoglobulin G1 (IgG1; Immunotech, Hialeah, Fla.), which was used as the isotype control. Following incubation on ice for 30 min in the dark, the cells were washed thrice with permeabilization buffer, resuspended in 300 µl of permeabilization buffer, and analyzed with a Becton Dickinson FACSCAN instrument (BD Biosciences) in conjunction with CellQuest software (BD Biosciences) for flow cytometric analysis.
Total RNA extraction. Total RNA was extracted with RNEasy Midiprep kits per the manufacturer's protocol (Qiagen, Valencia, Calif.). RNA concentrations and purity were measured by spectrophotometry, and RNA quality (absence of RNA degradation) was assessed by gel electrophoresis. The RNA concentration was adjusted to the levels required for subsequent microarray experiment protocols by concentration in a SpeedVac (Savant Instruments, Holbrook, Calif.). RNA samples (6 to 7 µg/µl) were stored in 100 µl of Tris-EDTA buffer at 80°C.
Real-time RT-PCR quantitation of viral RNA and cellular RNA. Quantitation of HIV viral mRNA was carried out by real-time PCR with an ABI Prism 7000 instrument (Applied Biosystems, Foster City, Calif.). A housekeeping gene, glyceraldehyde phosphate dehydrogenase (GAPDH), whose expression was unchanged in ACH-2 cells before and after PMA treatment was used as a normalization control. RNA from the samples was subjected to DNase treatment to remove contaminating DNA, and the DNase was inactivated with the DNase Free kit (Amersham Biosciences, Piscataway, N.J.) according to the manufacturer's protocols; 2 µg of RNA was reverse transcribed with the Taqman reverse transcription (RT) kit from Applied Biosystems per the manufacturer's specifications. Briefly, the reaction mixture (50 µl) was incubated at 65°C for 5 min, followed by 37°C for 45 min and 94°C for 5 min, and then cooled on ice, and 1/40th aliquots of the corresponding samples were used in a real-time PCR with Taqman probes labeled with reporter dye label FAM and quencher dye label TAMRA at the 5' and 3' ends, respectively.
Primer probe pairs were designed with PrimerExpress (Applied Biosystems). The reactions were carried out in triplicate for each time point, and the changes observed were normalized to expression of GAPDH for each time point. Melting-curve analysis and an assessment of the efficiency of the PCRs over a given concentration range were performed to determine if any errors occurred during the reactions, for example, primer dimer formation or poor priming effects. Real-time PCRs were set up as described above with primers and probes specific for early (multiply spliced mRNA) and late (unspliced mRNA). The sequences of the 5' and 3' primers used to quantitate early mRNA were 5'-CGAAGAGCTCATCAGAACAGTCA-3' and 5'-TTGGGAGGTGGGTCTGCTT-3'. The sequence of the labeled probe was 5'-CTTCTCTATCAAAGCAGACCCACCTCC-3', and it overlapped the splice site of the HIV-1 Rev sequence. The sequence detection primers for unspliced, late RNA were SK38 and SK39 from the HIV-1 Gene Amplimer kit (Applied Biosystems). A TAMRA-labeled probe identical in sequence to SK19 (Applied Biosystems) was used for real-time PCR quantitation of the late viral RNA species. Standards from the kit were diluted to calculate the copy number of virus based on gag mRNA concentrations.
Real-time RT-PCR analysis was also carried out for selected cellular genes with gene-specific primer-probe pairs and Taqman detection primers. Differences in mRNA expression in uninduced ACH-2 samples and the corresponding A3.01 samples were determined with the protocol described for quantitation of viral mRNA. Real-time RT-PCR quantitation was performed for genes PSMC5, p44s10 (proteasome subunits) Egr1 (early growth response 1), HDAC1 (histone deacetylase 1), NK4, EIF4, and SFRS3 to confirm that these genes were differentially expressed in the latently infected ACH-2 cells compared to the uninfected A3.01 parental cells. Primer-probe pairs specific for each gene were designed with PrimerExpress (Applied Biosystems). The sequences of the detection primers and probes for each gene are available as supplemental material (Table S1).
Microarray studies. Total RNAs obtained from induced chronically infected and the corresponding uninfected parental cells were used for microarray experiments. For each time point, RNA from the induced chronically infected ACH-2 cells and RNA from the corresponding induced, uninfected A3.01 cells were compared to minimize effects due to PMA induction. Microarrays were obtained from the National Cancer Institute Microarray Facility, Advanced Technology Center (Gaithersburg, Md.). The microarrays (Hs. UniGem2) contained 10,395 cDNA spots on each glass slide. The cDNAs were selected for spotting on the slides based on their known or probable involvement in oncogenesis, signal transduction, apoptosis, immune function, inflammatory pathways, cellular transport, transcription, protein translation, and other important cellular functions. A number of expressed sequence tags (ESTs) from unknown genes homologous to known genes and cDNAs encoding housekeeping genes were also included in these gene sets.
For each time point, 50 µg of total RNA from PMA-induced ACH-2 cells and 70 µg of total RNA from PMA-induced A3.01 cells was labeled with indocarbocyanine (Cy3)-dUTP and indodicarbocyanine (Cy5)-dUTP, respectively, as previously described (34, 60). Larger amounts of RNA were used for Cy5 labeling to minimize the disparities in dye incorporation. Each sample of RNA from PMA-induced, infected cells from a particular time point was compared with RNA from the corresponding PMA-induced, uninfected cells from the same time point for subsequent hybridization to the same array to ensure accurate comparisons and to eliminate interarray variability. The labeled cDNAs were then combined and purified with MicroCon YM-30 (Millipore, Bedford, Mass.) spin column filters to remove any unincorporated nucleotides; 8 to 10 µg each of Cot-1 DNA (Boehringer Mannheim, Indianapolis, Ind.), Saccharomyces cerevisiae tRNA (Sigma), and poly(A) (Amersham Biosciences) were added to the reaction mixture and heated at 100°C for 1 min. Hybridization of the labeled cDNA to the microarray was carried out at 65°C overnight, followed by washes with 1x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), 0.2x SSC, and 0.05x SSC. The slides were dried by centrifugation at 1,000 rpm for 3 min and then scanned as described below. RNA samples from three identical but independently conducted time course experiments were tested. Microarray experiments were performed at least twice for each time point (technical replicates) of each experiment. We also compared AZT-treated ACH-2 cells to untreated ACH-2 cells to determine whether any differences in gene expression might be due solely to AZT.
Microarray experiments for each chronically infected cell line (J1.1, U1, and ACH-2) were also conducted as per the protocol described above. Samples from eight independent experiments per cell line were used for studying the gene expression patterns in chronically infected cell lines. To compensate for dye labeling bias that might be due to differences in Cy5 and Cy3 labeling efficiency and preferential dye incorporation by some mRNA species, RNA from the same samples labeled with Cy5 (70 µg of RNA) and Cy3 (50 µg of RNA) were cohybridized to the same array and scanned, data were analyzed for all the cell lines studied with identical filtering and statistical tests, and genes showing dye incorporation bias were eliminated from further analysis as described below.
Microarray scanning and data analysis. The slides were scanned with an Axon GenePix 4000 scanner (Axon Instruments, Union City, Calif.). The photomultiplier tube values were adjusted to obtain equivalent intensities at both wavelengths used, 635 and 532 nm for the indodicarbocyanine and indocarbocyanine channels, respectively. Image analysis was performed with GenePix analysis software (Axon Instruments), and data analysis was performed with the MicroArray Database (mAdb) system hosted by the Center for Information Technology and Center for Cancer Research at NIH (http://nciarray.nci.nih.gov). Each array was normalized with Lowess normalization (71). Normalization across arrays over the time course was not feasible because no particular time point could be established as a median for cellular gene expression data, since viral gene expression and associated host cell gene expression were different for different time points.
For each time point, there were at least six data sets (two technical replicates times three biological replicates). For each time period (0.5 to 8 h, 12 to 24 h, and 48 to 96 h postinduction), there were at least 18 datasets (at least three time points per time period). Only arrays that met initial spot size and intensity criteria and whose normalization ratio (between the two signal intensities) was close to 1 (0.85 to 1.15) were analyzed. Filtering criteria were as follows. For each spot, signal intensity must be at least twice that of the background intensity; each gene must have values in at least 70% of the arrays; and each array must have values for at least 70% of the gene spots. Genes that showed dye labeling bias in a particular cell line after normalization were excluded from that gene set prior to further analysis. This was determined with a one-sample t test on mean log ratios for replicate arrays with the same sample labeled with both Cy3 and Cy5.
The raw data from these experiments have been deposited into the National Center for Biotechnology Informations Gene Expression Omnibus database (www.ncbi.nlm.nih.gov/geo/) under accession numbers GSE1441 and GSE1443.
Statistical analysis. Comparison of the expression profiles of infected versus uninfected (both induced by PMA) cell lines at a given time point was performed with univariate parametric and multivariate permutation tests based on the one-sample random variance t statistic in BRB-ArrayTools (http://linus.nci.nih.gov/BRB-ArrayTools) (62). Since RNAs from infected and uninfected cell lines corresponding to the same time point were paired and cohybridized on the same array, interarray sources of variation were minimized and differential expression could be detected by a statistically significant nonzero mean log ratio in biologically independent replicates. All biological replicates that passed the filtering criteria described above were used in the analyses. Technical replicates were averaged. The random-variance model enabled variance information to be shared across genes without assuming that all genes had the same variance (69). For comparison of expression of latently infected versus uninfected cell lines, significance was based on P < 0.001 for a parametric one-sample random variance t test. For evaluation of differential expression between infected and uninfected cell lines at fixed times after induction, a multivariate permutation test based on the one-sample random variance t statistic was used in which the proportion of false discoveries was limited to 0.10 with 90% confidence (36, 62).
Hierarchical clustering analyses of the resulting data sets were done with the mAdb system as well as the Cluster and TreeView software programs (Stanford University, Stanford, Calif.). Since a gene exhibiting statistically significant differences in expression need not alter the physiology of a cell in a biologically meaningful way, an additional pathway analysis of the genes that showed significant differential expression was performed with analysis tools provided by the NIH mAdb database (http://nciarray.nci.nih.gov) and querying the database of the Cancer Genome Anatomy Project (CGAP) (http://cgap.nci.nih.gov/) with pathway information provided by the Kyoto Encyclopedia of Genes and Genomes (KEGG) (www.genome.ad.jp/kegg/) and Biocarta (www.biocarta.com) pathway databases. Gene ontology summary analyses, which allow grouping of genes based on their molecular function, were also performed to ensure that the changes observed in our studies were not due to general dysregulation of all genes studied. Where possible, the observed/expected ratio (O/E) for each functional class of genes was determined to ensure that observed changes in the number of genes differentially expressed within a gene class were greater than that expected by chance (O/E > 1). This provided stringent thresholds for pathway and functional classification of the differentially expressed genes.
Latency reactivation studies. Cells (ACH-2, J1.1, and U1) were seeded at a concentration of 2 x 105 cells/ml in 24-well plates in a volume of 1 ml. Clastolactacystin-beta-lactone, a proteasome inhibitor; resveratrol, an Egr1 activator; or trichostatin, a histone deacetylase inhibitor (Biomol Research Laboratories, Plymouth Meeting, Pa.) was dissolved in sterile dimethyl sulfoxide and further diluted with medium to obtain the desired final concentrations. The final concentration of dimethyl sulfoxide in contact with the cells was never greater than 0.001% at any dose tested. AZT (250 nM) was added to the chronically infected cells in order to inhibit p24 production that may be caused by low levels of actively replicating virus present along with the chronically infected cells and to ensure that any increases in p24 expression would be attributable to activation of latent provirus and not due to subsequent amplification via additional rounds of viral replication. Cells were incubated with different concentrations of either clastolactacystin-beta-lactone, resveratrol, or trichostatin at 37°C, and 200-µl samples of cell supernatant were collected at 24 h after treatment. Cells incubated with tumor necrosis factor alpha (0.5 µg/ml) in the presence of AZT served as a positive control. Cells treated with AZT alone served as a negative control. Cells not treated with AZT were also examined.
Samples were mixed with lysing buffer (10% Triton X-100; Sigma) to inactivate virus and diluted fivefold with sample diluent (1% bovine serum albumin, 0.2% Tween 20 in RPMI 1640). p24 expression was assayed by enzyme-linked immunosorbent assay with HIV-1 p24 antigen capture kits (AIDS Vaccine Program, Frederick, Md.) per the manufacturer's specifications. Briefly, plates were washed with plate wash buffer, and samples were added in duplicate wells. The samples (100 µl) were incubated for 2 h at 37°C. The plates were washed, and rabbit anti-HIV p24 antibody was added at a 1:400 dilution. Following incubation for 1 h, the plates were washed, and goat anti-rabbit IgG peroxidase-labeled antibody at a 1:300 dilution was added. The plates were incubated for 1 h at 37°C, followed by washing and addition of a two-component substrate. The substrate solution consisted of equal volumes of tetramethylbenzedine peroxidase substrate and peroxidase solution B (Kirkegaard and Perry Laboratories, Gaithersburg, Md.). Samples were incubated for 30 min at room temperature, and reactions were stopped by addition of 1 N hydrochloric acid solution. The absorbance was measured at 450 nm with a SpectraMax 250 spectrophotometer (Molecular Devices Corporation, Sunnyvale, Calif.). The samples were assayed in duplicate, and experiments were performed at least thrice with independent cell samples.
|
|
|---|
We treated the ACH-2 chronically infected cell line with PMA to trigger the initiation and completion of the lytic replication cycle. Virus production was confirmed by flow cytometry for the late HIV protein p24 (Fig. 1). p24 production was low in the absence of PMA (8.2%) in ACH-2 cells not treated with AZT. A3.01, the parental uninfected cell line, did not show any p24-specific staining over that of the isotype control (data not shown). At 6 h postinduction, 62% of the ACH-2 cells showed p24 production, and from 12 h up to 96 h postinduction, nearly all cells were positive for p24 production, indicating that a lytic HIV infection was under way in essentially all the cells in the culture (Fig. 1). By 48 h postinduction, flow cytometry analysis showed high levels of p24 production but with increased cell death due to cytopathic effects (data not shown).
![]() View larger version (27K): [in a new window] |
FIG. 1. Flow cytometric analysis of chronically infected ACH-2 cells before and after induction. Uninduced cells and cells from serial time points were fixed and permeabilized for intracellular p24 labeling. As an isotype control, cell samples were labeled with mouse IgG1. For each sample, 100,000 events were collected. Each sample histogram labeled for p24 (red) is overlaid with the control histogram labeled for the isotype control (green). (A) Uninduced ACH-2 cells, showing minimal p24 accumulation with 8.2% of cells infected. (B) ACH-2 cells at 0.5 h postinduction (p.i.) with 7.4% of cells positive for p24. (C) ACH-2 cells at 6 h postinduction, with 61.6% of cells infected. (D, E, and F) ACH-2 cells at 12, 18, and 24 h postinduction, respectively, showing complete infection. Flow cytometric analysis was performed on all batches of cells to ensure active replication of HIV following induction with PMA. Data from one induction experiment are shown. The data indicate that viral replication occurs in an ordered manner postinduction and that complete infection of cells is achieved within 12 h postinduction of chronically infected ACH-2 cells.
|
Viral mRNA expression by real-time RT-PCR. In order to confirm that production of early and late viral mRNAs was under way in lytically induced cells and to determine the relative proportions of early (multiply spliced; tat, rev, and nef) and late (unspliced; gag, pol, and env) HIV-1 mRNAs over the time course, we used real-time RT-PCR to quantitate the message classes. In our system, multiply spliced mRNA expression increased after lytic induction and showed a maximum 40-fold increase over uninduced ACH-2 cells at 8 h postinduction. Unspliced mRNA concentrations showed a gradual increase, with a 148-fold increase at 18 h postinduction (Fig. 2). The results indicate that viral RNA expression in our lytically induced cells followed known kinetic expression patterns (6, 35). There was a clear distinction in the peak expression of the early and late viral mRNAs, indicating that viral RNA expression followed a discrete temporal pattern in these cells.
![]() View larger version (18K): [in a new window] |
FIG. 2. Levels of expression of multiply spliced (MS) and unspliced (US) HIV-1 mRNA pre- and postinduction of chronically infected ACH-2 cells. Real-time RT-PCRs were carried out with Taqman probes specific for early (multiply spliced) and late (unspliced) transcripts of HIV-1 and tagged with FAM and TAMRA fluorescent dyes at the 5' and 3' ends, respectively. Reactions were performed in triplicate for each time point as described in Materials and Methods, and average values are shown. Maximal changes in mRNA levels for early transcripts (multiply spliced) were observed at 8 h postinduction. The changes for late transcripts (unspliced) showed a maximal increase at 18 h postinduction.
|
Data from the technical and biological replicates for each time point were normalized for each array as described in Materials and Methods. Statistical analyses to identify genes that were significantly differentially expressed were performed with univariate and multivariate random-variance one-sample t tests with all the replicates that passed the filtering criteria. Based on the statistical analyses of the genes showing altered expression, 131 genes showed altered expression even prior to induction; 1,740 of all the analyzable spots showed statistically significant (P < 0.001) altered gene expression at some point, either before induction or over the entire period of the lytic replication cycle (Fig. 3).
![]() View larger version (56K): [in a new window] |
FIG. 3. Hierarchical clustering of differentially expressed cellular genes before and after induction of chronically infected ACH-2 cells. The figure shows the hierarchical clustering of the cellular genes that showed significant differential expression (P < 0.001) across the time course (before induction up to 96 h postinduction) following reactivation of chronically infected ACH-2 cells as described in the text. Genes shown in red showed upregulation, those in green were downregulated, while those that did not show any change with respect to a normalized matched control are shown in black. The gray areas indicate missing data for the given gene and time point. The magnified panels indicate selected kinetic profiles that were seen before and following induction into active viral replication. (A) Upregulation of selected genes observed before induction; (B) upregulation of genes immediately following induction; (C) genes that are upregulated prior to induction and downregulated 12 to 24 h postinduction; (D) genes that are upregulated in the early stage following reactivation but downregulated in the intermediate stage; (E) genes that are downregulated before induction but upregulated in the intermediate stage followed by downregulation in the late stage (48 to 96 h postinduction).
|
The changes in gene expression observed during the late period (48 to 96 h postinduction) cannot be solely associated with the process of lytic replication because during this time period the cells showed cytopathic effects. However, the data are presented in order to provide the most complete possible view of the cellular environment following lytic replication. Some genes showed differential expression only during the early stages of lytic replication and returned to the levels observed in the uninfected cells, while certain other genes were initially unaltered and showed differential expression at later time points.
Certain groups of genes showed similar patterns of regulation following induction (Fig. 3). In the early time period, 1,334 genes were differentially expressed. In the intermediate time period, 756 genes were differentially expressed. The late time period (48 to 96 h postinduction) showed the least change, with 566 genes exhibiting significantly altered expression (P < 0.001). Many of the genes that were differentially expressed in the early time period also showed either a similar or the opposite trend in their expression patterns during the other time periods; hence, some genes were included in the analysis of both time periods.
A number of discrete patterns of gene regulation were observed. Several cellular genes showed distinct temporal expression patterns during the lytic replication cycle, an expected finding, but more interestingly, a smaller number of genes appeared to be differentially expressed in the latently infected ACH-2 cells compared to their parental uninfected cells even before induction of the lytic cycle. A total of 131 genes showed a significant change (P < 0.001) in their expression prior to induction. They included genes encoding transcription factors, components of proteasomes, factors that control immune function, apoptosis, and other functional classes. For gene classes that were annotated in the gene ontology database (www.geneontology.org) (5), the O/E ratio for the number of genes within a functional class that were differentially expressed was set at greater than 1 (O/E > 1) to eliminate functional classes where the number of differentially expressed genes was not greater than that randomly expected. However, not all genes that were significantly differentially expressed are annotated in the database. Extensive literature studies for functional significance and classification were conducted in such cases to ensure that important classes of genes were still included in the analyses. An abbreviated listing of the genes grouped according to known functions that were differentially expressed before induction is given in Table 1.
|
View this table: [in a new window] |
TABLE 1. Functionally related genes that were differentially expressed prior to induction in chronically infected ACH-2 cellsa
|
Real-time RT-PCR analysis was performed on a selected set of cellular genes with gene-specific sequence detection primers and probes to determine if the microarray results correlated with the actual normalized differences in the corresponding ACH-2 and A3.01 samples prior to induction. For selected genes that were found to be statistically and/or biologically significant, differential expression was confirmed by RT-PCR quantitation (supplemental material, Table S3).
Trends seen in pathway profiles during lytic replication. While a number of genes that had altered expression during lytic replication could be classified based on their functionality, an alternative method was also used to determine if the genes that showed differential expression grouped into known cellular pathways. Genes encoding components of several distinct pathways were regulated in a coordinated fashion during lytic viral replication. Of the 1,740 genes that were differentially expressed over the lytic infection cycle in ACH-2 cells, 697 were assigned to various known pathways with the CGAP pathway databases, used in conjunction with the MicroArray database mAdb data analysis tools (http://nciarray.nci.nih.gov). Based on the information in the database, these genes were found to lie in a total of 385 known pathways. The remainder of the genes could not be classified into any known pathways listed in the CGAP database.
The principal pathways affected over the time course and classified based on the number of genes involved in a particular pathway that are differentially expressed are shown in Fig. 4 (and in the supplemental material, Table S2). We observed that a number of pathways involved in signaling, cell cycle, and transcription showed maximum changes even prior to induction. During the early stage (0.5 to 8 h postinduction), a number of metabolic pathways and signaling pathways showed similar patterns. The intermediate stage (12 to 24 h postinduction) showed maximal changes in pathways involved in immune response modulation and cell survival. These patterns correlated well with the levels of early (multiply spliced) and late (unspliced) HIV protein gene expression and known viral replication effects on the host cell (32, 44, 59).
![]() View larger version (40K): [in a new window] |
FIG. 4. Trends seen in pathways that show differential expression before and after induction of chronically infected ACH-2 cells. Pathway profiles observed prior to induction and following reactivation of ACH-2 cells with PMA over a period of 96 h. The figure shows the number of genes in each pathway that were differentially expressed in a particular pathway. (A) Pathways that were maximally altered prior to induction. (B) Pathways that showed maximum change during the early phase of the lytic cycle (0.5 to 8 h postinduction). (C) Pathways that showed maximal change during from 12 to 24 h postinduction. Most pathways did not show any change during the period from 48 to 96 h postinduction. The groups above are a selected representation of the various pathways that changed differentially prior to induction and/or over the time course studied. Classification of the altered genes into various pathways was performed with the CGAP pathway databases.
|
Experiments were performed on eight independent cell samples for each cell line, and similar parameters were applied for filter criteria, gene selection, and statistical analysis as with the ACH-2 cell line, as described in Materials and Methods. In these sets of experiments, 24 arrays were analyzed, and 8,902 of the 10,395 gene spots passed the selection criteria. Statistical analysis of expression ratios of AZT-treated ACH-2 cells and untreated ACH-2 cells over their respective A3.01 controls did not show any significant differences in gene expression profiles, indicating that changes in gene expression were not due to low levels of actively replicating viral population (data not shown). Genes that had shown differential dye incorporation in our dye labeling bias experiments were excluded from each cell line data (43 genes in ACH-2, 18 genes in Jurkat, and 22 genes in U937 cells). Upon analyzing the resulting data sets, we found that 131 genes were differentially expressed in ACH-2, 65 genes were differentially expressed in J1.1, and 155 genes were differentially expressed in U1 cells compared to their uninfected AZT-treated parental cell lines. While stringent statistical thresholds were used for genes that were differentially expressed in each cell line, we reasoned that if a gene showed up as significantly differentially expressed (P < 0.001) in at least one cell line, then the expression data for that gene in the other cell lines may be important even if not found statistically significant for that cell line. We found that ACH-2 and J1.1 were more similar in their profiles, which may be due to their common T-cell lineage, compared to U1, which is a promonocytic cell line (Fig. 5).
![]() View larger version (32K): [in a new window] |
FIG. 5. Hierarchical clustering of genes that showed differential expression across three chronically infected cell lines prior to induction. Hierarchical clustering of differentially expressed genes that showed a significant change in expression (P < 0.001) in the chronically infected cell lines ACH-2, U1, and J1.1. Genes shown in red were upregulated, those in green exhibited downregulation, and black indicates normal expression. Gray areas indicate missing values. Many genes were altered similarly across the cell lines. Each cell line also showed some unique patterns of cellular expression. Data are the averages of values from eight independent samples per cell line. The magnified portions of the cluster highlight some of the patterns of gene expression across the cell lines. (A) Genes that were upregulated in all three cell lines; (B) genes that were downregulated in all three cell lines; (C) genes that were upregulated in ACH-2 and J1.1 and downregulated in U1 cells; (D) genes which showed no significant similarity in their expression in the three cell lines.
|
Reactivation of chronically infected cells by targeting specific gene classes. We reasoned that the differential expression of certain cellular genes in latently infected cells might be involved in maintaining the virus in latency and so hypothesized that treating latently infected cells with agents targeting the products of differentially expressed genes may force the virus out of latency. This may be important because reactivation of proviral HIV into lytic infection by targeting cellular factors, notably including cellular factors differentially expressed in latent, chronically infected cell lines, may provide new approaches to reduce or eliminate latent viral reservoirs. To test our hypothesis, we focused on the proteasome class of genes and Egr1 due to the availability of specific agents. We also found other differentially expressed genes such as the histone deacetylases, HDAC1 and HDAC2, which may represent good targets for reactivation of viral replication. Treatment of latently infected cells with agents that alter histone acetylation has already been shown to trigger lytic replication in latently infected cells (21, 40, 66).
We observed that the proteasome class of genes showed increased expression in ACH-2 cells even prior to induction but not in J1.1 and U1 cells. Hence, we sought to determine whether inhibition of proteasomes would induce latent provirus into lytic replication. We studied the effects of different concentrations of a proteasome inhibitor, clastolactacystin-beta-lactone (19), on p24 protein levels in chronically infected ACH-2 cells. A previous study has shown that proteasome inhibitors are capable of increasing the efficiency of an acute HIV infection (58). We therefore studied the effect of the proteasome inhibitor in the presence of 250 nM AZT (50% inhibitory concentration, 50 nM) to ensure that any increases in p24 levels were not due to increased efficiency of infection by actively replicating virus present in chronically infected cells or to subsequent rounds of viral replication. ACH-2 cells treated with AZT alone showed slightly lower p24 amounts than cells that were not treated with AZT, consistent with the observation that the chronically infected cells supported a low but detectable amount of ongoing viral replication.
Tumor necrosis factor alpha, a known activator of latently infected cells, caused a 20- to 24-fold increase in p24 expression in ACH-2 cells treated with AZT, and 5 µM clastolactacystin-beta-lactone caused a 10- to 15-fold increase in p24 expression compared to AZT-treated ACH-2 cells, indicating that clastolactacystin-beta-lactone activated latent provirus into lytic replication. The increases in p24 expression observed in cells treated with clastolactacystin-beta-lactone were dose dependent, indicating a drug-related response (Fig. 6A). This effect was not observed in similarly treated J1.1 and U1 cells (data not shown), which is consistent with the finding that the expression of a large number of proteasomal subunits was not upregulated in these chronically infected cells. We also tested PS-341, a highly specific proteasome inhibitor that binds to the beta-5 subunit of the proteasome and is approved for clinical treatment of multiple myeloma (1, 3). However, due to the extreme cytotoxicity of PS-341 even at a low concentration (100 nM), we were unable to observe any changes in p24 expression under our experimental conditions (1 to 1,000 nM) in latently infected cells (data not shown).
![]() View larger version (16K): [in a new window] |
FIG. 6. Effects of specific agents on HIV p24 production in chronically infected ACH-2 cells. Different concentrations of (A) the proteasome inhibitor clastolactacystin-beta-lactone, (B) resveratrol, an Egr1 activator, or (C) trichostatin, a histone deacetylase inhibitor were tested in chronically infected ACH-2 cells treated with 250 nM AZT. Samples were collected 24 h after addition of agent, and p24 concentrations were determined by enzyme-linked immunosorbent assay. p24 production in cells treated with tumor necrosis factor alpha (TNF ) was used as a positive control in ACH-2 cells compared to control (AZT-treated) cells. p24 production in untreated cells (no AZT) was also determined. Experiments were performed in triplicate, and values are representative of three independent experiments. Microarray data for the specific genes targeted are shown for each agent tested.
|
Considerable prior work has been done on the influence of chromatin structure and histone acetylation status on HIV gene expression (21, 40, 66). These studies postulated that agents that can alter the histone acetylation state of the host cell would activate lytic HIV replication, and this was indeed found to be the case. However, there has been no indication that the cellular machinery that maintains histones in the deacetylated state might be upregulated in cells latently infected with HIV. In our expression data, we observed that HDAC1 and HDAC2 were upregulated more than twofold prior to induction. Also, the gene encoding YY1, known to repress the HIV-1 long terminal repeat through histone deacetylase recruitment (17) was also significantly upregulated prior to induction in ACH-2 cells. We confirmed that treating cells latently infected with HIV with trichostatin, an agent that targets this class of differentially expressed gene products, activates lytic HIV replication (Fig. 6C). Studies with agents that target specific cellular factors not only provide new approaches to cause viral reactivation, but also validate the biological relevance of these differentially expressed gene classes.
|
|
|---|
We studied the responses of 9,127 cellular genes following induction of a lytic replication cycle in cells chronically infected with HIV. Of these, 1,740 genes showed statistically significant differential expression at least once over the period of the lytic cycle. We classified the differentially expressed genes based on function and into known pathways by querying pathway databases hosted by CGAP. Classification of differentially expressed genes into known pathways in addition to functional classification provides valuable information regarding pathways that show differential expression in functionally unrelated genes but change similarly across the lytic replication cycle.
Our studies indicated that, following induction of cells latently infected with HIV, clear, temporally ordered changes in the cellular gene expression pattern accompanied the completion of the viral lytic cycle. Some of these changes may result from incidental changes extraneous to the processes of viral replication, for example, nonspecific cytopathic effects; some may reflect specific modifications to host cell physiology produced by the targeted action of viral gene products aimed at enhancing the ability of the host cell to support viral replication; and some may represent attempts by the host cell to blunt the effects of viral replication. Several genes previously implicated as having important roles in viral replication were differentially expressed, suggesting that at least some of the observed changes may be important for HIV pathophysiology. Other studies of acute infection have also shown some of these effects (for example, differential expression of genes encoding transcription factors and factors involved in apoptosis and signaling), indicating that the changes that accompany lytic replication are not too different from those associated with acute infection (29, 67).
While a number of gene classes did show similar trends during our lytic replication study, which correlated well with findings from other groups studying acute HIV infection, a direct and quantitative comparison of all gene classes may not be accurate due to the different controls employed in studies conducted by other groups, which may cause a skewed interpretation of the results. Also, since most other published studies looked at acute infection in a small number of cells against a background of a large population of uninfected cells, many gene classes that showed up as differentially expressed in our study may have failed to pass the filtering criteria due to high background noise in those studies.
In our microarray experiments, genes encoding histone deacetylases were upregulated prior to induction, correlating with studies that histone deacetylase inhibitors end viral latency. Genes involved in the mitogen-activated protein kinase signaling pathway are also upregulated during lytic replication, which corroborates studies that show that mitogen-activated protein kinase inhibitors suppress lytic replication (70). However, due to the innate complexity of many signaling pathways and the presence of feedback mechanisms and because various posttranslational regulatory mechanisms may be involved, a direct correlation may not always be present. The identification of differentially expressed gene classes may thus provide only a helpful starting point for further detailed research into the mechanisms that influence viral latency and the initiation of lytic replication. Assessing the precise significance of the many changes in host cell gene expression pattern that accompany the completion of the lytic infection cycle will necessarily require much additional work. However, the identification of many cellular genes that are differentially expressed during HIV lytic replication may offer additional targets for agents aimed at blocking HIV replication.
While changes in cellular gene expression during the viral lytic cycle were expected, a more important finding was that a significant number of genes showed altered expression even prior to induction. These changes could result either from the direct action of low levels of certain HIV gene products expressed in latency or, more likely, because changes in the cellular pattern of gene expression were selected for as the chronically infected cell lines were isolated and cloned. However, regardless of the mechanisms through which the alterations in cellular gene expression developed, the differential expression of these cellular genes may play a role in HIV latency and the maintenance of latency. While the differential expression of cellular genes in the latently infected cells may not be directly caused by the latent provirus, it is reasonable to postulate that the products of at least some of these differentially expressed genes may collaborate in the maintenance of the latent state, even though those genes were differentially expressed in the chronically infected cell lines as a consequence of the development and cloning of the chronically infected cell lines. Treating latently infected cells with agents targeting the differentially expressed cellular genes may help eject HIV from the latent state, thus providing novel approaches for the reduction of latent viral reservoirs with agents that target the differentially expressed cellular genes.
The effects seen prior to induction are probably not due merely to low levels of HIV replication in the latently infected cell population. Many genes that were downregulated during the lytic phase showed increased expression prior to induction and vice versa, for example, some transcription factors and DEAD-box proteins. Our flow cytometry studies showed only a low level (8%) of active replication prior to induction. In microarray experiments, signals originating from such a small fraction of the cells are likely to go undetected. Consistent with this, our microarray data for AZT-treated, latently infected ACH-2 cells yielded the same gene sets as obtained for untreated ACH-2 cells, indicating that changes in expression ratios were due mainly to the latently infected population and not to the low level of cells harboring actively replicating virus. Also, most of the cellular genes that were differentially expressed prior to induction also showed a gradual decrease in the level of up- or downregulation as the lytic phase progressed. It is therefore unlikely that most of the changes in gene expression observed prior to induction in chronically infected cell lines are due to actively replicating virus present in only a small fraction of cells. Hence, the effects on cellular gene expression seen prior to induction are most likely associated with the state of latent infection in these cell lines.
Based on our hypothesis that targeting the products of genes that were differentially expressed during viral latency would lead to viral reactivation, we showed that specific targeting of certain genes and gene classes can indeed force the virus out of latency. The mechanisms contributing to the observed effects may not be clear in all cases. For example, the mechanisms by which proteasome inhibitors end viral latency are not immediately obvious, though some available data may suggest an explanation. Studies have shown that HIV TAR sequence (present in the 5' untranslated region of all HIV mRNAs) serves as the substrate for proteasomal degradation (28). Also, Nef, an early protein which is present in the host cell even before proviral integration, is processed by proteasomes, which provides many of the cytotoxic T-lymphocyte epitopes of Nef (45). Studies have also shown that Tat has an inhibitory effect on proteasome function by physically binding to 20S proteasome, thereby reducing viral protein degradation (59).
We propose that in chronically infected ACH-2 cells, the higher expression of proteasomes may lead to increased degradation of HIV mRNA. This effect, along with insufficient levels of full-length Tat and Nef proteins, would tend to inhibit HIV lytic replication. In such a case, addition of a proteasome inhibitor may release the suppressive effect of proteasomal degradation, potentiate the function of Tat by binding to the active sites on the proteasome, and thus allow accumulation of sufficient early viral proteins to trigger viral reactivation. Proteasome inhibition also inhibits NF-
B activity, which should decrease viral replication. However, it must be noted that even during active viral replication, use of a proteasome inhibitor increases viral replication efficiency (58). Binding of Tat to proteasomes has been shown to affect the peptidase activity (32). This may explain why we did not observe viral reactivation with PS-341 (2), which specifically inhibits the chymotryptic activity of proteasomes (3). A broader range of proteasomal inhibition may be required for viral reactivation. Clastolactacystin-beta-lactone, which caused viral reactivation, affects all three catalytic activities of the proteasome. The broad range of effects that the different active sites of proteasomes have on different pathways will necessarily make understanding the crucial mechanisms underlying these phenomena much more difficult.
Egr1 was downregulated during viral latency in ACH-2 cells. Resveratrol, an Egr1 activator, caused viral reactivation. Increased levels of Egr1 lead to growth arrest (G1/S and G2/M) (52), similar to that observed during viral replication. Downregulation of Egr1 may lead to decreased Egr1 levels, preventing growth arrest and creating conditions unfavorable for viral replication. The effect observed with resveratrol may be a combination of manipulating the latent cellular environment and mimicking the effects of active HIV replication on the cell cycle, providing favorable conditions for viral replication. Resveratrol has been shown to not affect the levels of NF-
B in other cell lines studied (52), suggesting that mechanisms other than NF-
B modulation may contribute to viral latency and reactivation. Our studies with specific agents that target cellular factors encoded by genes that are differentially expressed during viral latency not only provide new approaches to target viral latency but also provide support for the biological relevance of the genes that are differentially expressed.
Our studies in latently infected cell lines identified a common set of genes that were differentially expressed in all three cell lines studied. Notable among them was the downregulation of Cdc42 and Lyn, which are necessary for viral replication and maturation. Cdc42 and Lyn have been shown to be critical for activation of Nef-associated kinase (PAK2) and for binding of the PXXP motif, respectively (43, 57). Their downregulation in cell lines of different origins indicates that they may be critical for the cellular maintenance of viral latency. We were unable to test the effects of targeting the gene products of genes that were commonly differentially expressed in all cell lines because no specific agents for Cdc42 or Lyn are currently available.
Interestingly, three genes associated with acute myeloid leukemia, MNDA, CEBPalpha, and Meis1, were also downregulated in all the latently infected cell lines, similar to cdc42 and Lyn, indicating a new set of genes previously not related to HIV pathogenesis that may be involved in latency maintenance. MNDA also has been shown to bind with LANA (latency-associated nuclear antigen), a human herpesvirus 8 protein, upon viral reactivation by alpha interferon (27). There may be interactions between HIV proteins and MNDA too, which may be inhibited in latently infected cells due to downregulation of MNDA, thus preventing viral replication.
A small set of 20 genes were also upregulated in all three cell lines. While none of these has been previously linked to HIV infection or latency, some of them have been shown to interact with other viruses. For example, CDT1 inhibition is observed in active cytomegalovirus infection, which has been proposed as a mechanism by which cellular dysregulation due to cytomegalovirus occurs (7, 68). In our studies, CDT1 expression was upregulated in all three cell lines, suggesting that the cell may be maintaining latency by modulating certain cellular replication check points. Further research into the mechanisms of latency maintenance and reactivation and future therapeutic advances may be facilitated by the development of agents that target the genes that are similarly differentially expressed in these cell lines.
The existence of long-lasting HIV reservoirs is the principal barrier preventing the eradication of HIV infection (8, 49). Since there are currently no known phenotypic differences between latently infected and uninfected cells (10), gene expression profiling may provide crucial evidence of differences in the cellular environment between the two populations. Our data also suggest that specific targeting of genes showing differential expression in latently infected cells may force viral reactivation.
In summary, our studies show that a number of cellular genes and pathways are altered in chronically infected cell lines during viral latency. Following induction of a lytic replication, a number of genes and pathways are altered in a temporally ordered manner in the presence of actively replicating HIV. Some of these genes have been implicated in viral pathogenesis previously. However, a number of the differentially expressed genes have not been previously associated with HIV replication and/or viral latency. Our studies show that specific targeting of cellular gene products based only on carefully chosen gene expression data can be used to force HIV out of latency. In a clinical setting, triggering lytic replication in latently infected cells while maintaining good control of viral replication with antiretroviral drugs may potentially eliminate viral reservoirs. While the compounds used in our studies to eject virus from latency may not have immediate clinical utility, other agents targeting host cell gene products differentially expressed in latently infected cells may prove clinically useful.
Our studies provide a set of genes that were similarly differentially expressed in all the latently infected cell lines, which may yield a greater understanding of the functions of certain known cofactors and expand our knowledge of the factors involved in latency maintenance. Our studies also indicate new classes of cellular genes that may constitute good targets for agents that inhibit lytic viral replication. Such agents targeting host cell functions may be particularly helpful in blocking HIV replication by virus that is resistant to existing agents that target viral gene products. Our work thus may suggest new approaches to inhibit viral replication and new approaches to activate lytic infection, which may help decrease or eliminate latent viral reservoirs.
We are very grateful to Richard Simon for advice on statistical analyses of microarray data. We thank Mike Lu and Ann Gordon for assistance. We thank John Coffin, Kuan-Teh Jeang, Paul Meltzer, and Richard Koup for critical reading of the manuscript.
Supplemental material for this article may be found at http://jvi.asm.org. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»