Previous Article | Next Article ![]()
Journal of Virology, January 2003, p. 1392-1402, Vol. 77, No. 2
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.2.1392-1402.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Svetlana A. Mikheeva,
Gemma C. O'Keeffe, Michael G. Katze, Roger E. Bumgarner, Gary K. Geiss,
and James I. Mullins
Department of Microbiology, University of Washington School of Medicine, Seattle, Washington 98195-8070
Received 29 July 2002/ Accepted 11 October 2002
|
|
|---|
4,600 cellular RNA transcripts were assessed in CD4+-T-cell lines at different times after infection with human immunodeficiency virus type 1 strain BRU (HIV-1BRU) using DNA microarrays. We found that several classes of genes were inhibited by HIV-1BRU infection, consistent with the G2 arrest of HIV-1-infected cells induced by Vpr. These included genes involved in cell division and transcription, a family of DEAD-box proteins (RNA helicases), and all genes involved in translation and splicing. However, the overall level of cell activation and signaling was increased in infected cells, consistent with strong virus production. These included a subgroup of transcription factors, including EGR1 and JUN, suggesting they play a specific role in the HIV-1 life cycle. Some regulatory changes were cell line specific; however, the majority, including enzymes involved in cholesterol biosynthesis, of changes were regulated in most infected cell lines. Compendium analysis comparing gene expression profiles of our HIV-1 infection experiments to those of cells exposed to heat shock, interferon, or influenza A virus indicated that HIV-1 infection largely induced specific changes rather than simply activating stress response or cytokine response pathways. Thus, microarray analysis confirmed several known HIV-1 host cell interactions and permitted identification of specific cellular pathways not previously implicated in HIV-1 infection. Continuing analyses are expected to suggest strategies for impacting HIV-1 replication in vivo by targeting these pathways. |
|
|---|
Recently, the study of pathogen-host interactions has been greatly aided by gene expression analysis using high-density microarrays. Microarrays consist of ordered sets of cDNAs or oligonucleotides, each representing a single gene, immobilized as spots on solid substrates. The level of expression of each gene is then investigated by hybridization with probes derived from mRNA isolated from cell populations of interest. In this way, expression levels of thousands of genes can be measured simultaneously, generating a global picture of the transcriptional changes effected by pathogens such as bacteria, viruses, and parasites (reviewed in references 2, 23, and 32).
Previous microarray studies on HIV-1 have used infection of cell lines at low multiplicity of infection (MOI) (MOI, 0.1 to 0.5) to study changes up to 3 days after infection (8, 15, 39) or heterogeneous target cell populations (45). These experimental protocols result in nonsynchronous infection and in selection of certain cells resistant to infection, as several rounds of infection are needed to derive a completely infected population. This reduces sensitivity and specificity for the detection of gene expression changes in infected cells. The goal of our study was to analyze global changes that occur in a near-synchronous infection cycle of HIV-1 to increase our chances of finding relevant changes induced by HIV-1 infection. To this end, we used (i) infection at a high MOI with sampling within the first 24 h postinfection, in an effort to study a single replication cycle within synchronously infected cells; (ii) cell lines, to minimize heterogeneous host cell responses to infection and to allow precise enumeration of infected cells; and (iii) several different T-cell lines, to study variability and consistency of changes induced by viral infection. This approach has allowed us to identify cellular pathways not previously implicated in HIV-1 infection.
|
|
|---|
HIV-1BRU was obtained from J.-M. Bechet and L. Montagnier through the ARP. Virus stocks were made by exposing CCRF-CEM cells (107/ml) to virus at an MOI of 0.01. After 1 h, cells were washed and then cultured at 106/ml. After 3 days, culture medium was replaced, and cells were cultured for an additional 24 h before harvesting the virus-containing supernatant. Infectivity was determined on the same cells and was generally
107 50% tissue culture infective doses/ml. Supernatants were also obtained from mock infections, with cells treated as described above but without adding infectious virus. Heat-inactivated virus was prepared by heating for 2 h at 56°C in a water bath with frequent mixing. All mock and virus preparations were kept on ice until added to cells (within 30 min).
Infections. Infection of T-cell lines with HIV-1BRU was performed at an MOI of 2 at 37°C (Table 1). After 1 h of exposure to virus, cells were washed with phosphate-buffered saline and cultured at 106 cells per ml in RPMI with 10% fetal bovine serum. Mock cultures were treated the same way, substituting virus-containing supernatant with non-virus-containing tissue culture supernatant. Cells were harvested for RNA extraction at the indicated times by washing twice with phosphate-buffered saline and resuspending in 4 M guanidinium thiocyanate-25 mM sodium citrate-0.5% Sarkosyl-0.1 M 2-mercaptoethanol.
|
View this table: [in a new window] |
TABLE 1. CD4+-T-cell line infection experimentsa with HIV-1BRU
|
Hybridization probes. Isolation of total RNA was performed using a modified guanidinium thiocyanate method (7). Oligotex (Qiagen, Valencia, Calif.) was used to isolate mRNA according to the manufacturer's protocol. Fluorescence-labeled cDNA probes were generated and hybridized to the slides as described by Geiss et al. (14).
Microarray construction and analysis. The cDNA microarrays used in this study were generated in-house, and a description of the cDNA clones and construction has been published (14). Each slide contained 9,216 spots, corresponding to 4,608 genes spotted in duplicate. Included in this number was a set of 384 selected control genes spotted on each slide. These included HIV-1 genes used as positive controls, nonhuman genes used as negative controls, and a variety of genes selected from the Research Genetics 15K human gene set (Invitrogen, Carlsbad, Calif.).
Data analysis and selection of differentially expressed clones. For each of the infection conditions (Table 1), duplicate slides were hybridized with probes generated from the same RNAs but with the fluorescent labels reversed to control for dye-specific effects (15). Intensity values in Cy3 and Cy5 channels were extracted from each image with customized software (Spot-On Image, developed by R. E. Bumgarner and Erick Hammersmark). The resulting data files were entered into the Expression Array Manager (EAM) (developed by Jeff Furlong). EAM includes a database for storing microarray data, software tools for uploading data to the database, and software for uploading data to the Rosetta Resolver gene expression data analysis system (purchased from Rosetta Inpharmatics, Kirkland, Wash.). Spot-On data files were combined with information regarding slide and sample preparation and hybridization images. Resolver was then used to analyze hybridization intensity for each spot, to calculate ratios and standard deviations of replicate experiments after elimination of uninformative data, and to select statistically significant, differentially expressed genes. This included normalization for nonlinearities in the Cy3/Cy5 ratio as a function of intensity. The lists of genes selected in all experiments were combined to allow clustering using agglomerative clustering with complete linkage as the heuristic criterion to compute the between-cluster similarity measures (18). In complete linkage, the updated distance (or dissimilarity) of the merged cluster (A,B) is equal to the longest distance over all pairs of objects between clusters A and B.
RT-PCR quantification. Differential gene expression for a sample of eight genes was verified using reverse transcription (RT)-PCR with gene-specific primer pairs (Table 2). These genes were chosen based on either a substantial change in their expression levels by microarray analysis or their being regulated by the same transcription factor (SREBP [5]). Total RNA from a 24-h HIV-1BRU infection of CEMfh cells was used for quantification by multiplex RT-PCR using QuantumRNA Classic 18S internal standards (Ambion, Austin, Tex.) according to the manufacturer's protocol. PCR products were quantified on an Agilent (Palo Alto, Calif.) model 2100 Bioanalyzer DNA chip with internal standards as controls for variable amplification and loading of PCR products.
|
View this table: [in a new window] |
TABLE 2. Gene-specific primers for RT-PCR
|
|
|
|---|
![]() View larger version (43K): [in a new window] |
FIG. 1. Infection of CD4+-T-cell line CCRF-CEM with HIV-1BRU. (A and B) Cells were exposed to infectious virus at a multiplicity of infection of 2, heat-inactivated virus, or mock-infected cell supernatant. At 1 and 24 h postexposure, cells were fixed and permeabilized for intracellular p24 staining. Only cells exposed to the infectious virus stock showed detectable levels of p24 at 1 h postinfection (A) and increased p24 expression after 24 h (B). (C) Expression of HIV-1 RNA transcripts during infection. The graph depicts results for 17 HIV-1BRU infection experiments: a time course using CCRF-CEM cells and after 24 h of infection of CEMfh, CEMss, Jurkat, and SupT1 cells. Twelve PCR products representing different regions of the HIV-1BRU genome were spotted on each slide; nucleotide positions relative to the HIV-1BRU genome were 317 to 1281 (gag), 853 to 2471 (gag), 2179 to 3097 (pol), 3015 to 4363 (pol), 4337 to 5598 (pol), 5538 to 6974 (env), 6583 to 7404 (env), 7115 to 8631 (env), 8276 to 9225 (nef), 5402 to 5666 (tat), 7951 to 8272 (rev), 4619 to 5257 (vif), and 5124 to 5446 (vpr). Plotted are hybridization intensities of each of the 12 spots in the treated cells. HI (-) refers to mock infection versus heat-inactivated BRU, and HI (+) refers to heat-inactivated HIV-1BRU versus infectious HIV-1BRU.
|
Twelve control spots on each array contained PCR products spanning the length of the HIV-1 genome in segments. None of these control spots showed expression above background levels after hybridization with probes from mock-infected or heat-treated cells. The expression of HIV-1 gene segments in the HIV-1BRU-infected cells is shown in Fig. 1C. HIV-1 RNA was evident 1 h after exposure to virus, reflecting input virus bound to cells, and as expected, the abundance of these sequences decreased at 4 and 8 h postinfection as the genomic viral RNA was reverse transcribed into cDNA and then degraded. At 12 h postinfection the first increase in viral transcripts due to productive infection was observed, with maximum levels achieved by 24 h, and the intensity for each segment correlated with the length of the probe. RNA expression kinetics were consistent with the intracellular p24 expression kinetics measured by FACS (compare Fig. 1).
Several criteria were examined to select the set of statistically significant, differentially expressed genes (Fig. 2). P values chosen to define each experiment's set of regulated genes greatly influenced the numbers of genes selected, with close to 1,000 genes selected using a P of
0.05 and fewer than 300 selected using a P of
0.001 for the 24-h time points (Fig. 2A). The majority of changes observed at each time point were less than twofold for both down-regulation (Fig. 2B) and up-regulation (Fig. 2C). Changes were consistent across duplicate experiments for genes down- or up-regulated by both small or larger amounts, as the correlation coefficient between duplicate experiments was similar whether the experiments included all differentially regulated genes or just those regulated by at least 1.5- or 2-fold (Fig. 2D). However, decreasing the P cutoff did increase the correlation coefficient for duplicate experiments (Fig. 2D). On the basis of these analyses, we primarily used P values to select regulated genes, irrespective of the degree of change.
![]() View larger version (54K): [in a new window] |
FIG. 2. Genes found to be differentially regulated under different selection criteria. (A) Effect of signature P on number of genes selected. (B) Number of down-regulated genes in each experiment using a signature P of 0.005. Depicted are genes down-regulated by at least 2-fold (white bars) or 1.5-fold (grey bars), and all differentially regulated genes (black bars). (C) Number of up-regulated genes in each experiment using a signature P of 0.005 and depicted as in panel B. (D) Effect of signature P and relative change on the correlation coefficient between duplicate experiments.
|
0.005 is shown in Fig. 2B and C, respectively. Consistent with viral gene expression, few cellular genes were found to be differentially regulated as early as 1 h postinfection (5 genes), while the majority of changes were observed by 24 h postinfection (281 genes). The lists of differentially expressed genes selected at each time point were combined to allow clustering of genes over the entire time course. The experiments with heat-inactivated virus (HI- and HI+) were included as specificity controls. Because every time point was tested at least two times, only genes differentially regulated (P
0.005) in at least 2 of the 13 CCRF-CEM experiments were included in the following analyses.
Kinetics and cell-type-specific expression of host cell genes.
409 host cell genes fulfilled our requirements for differential expression (see supplemental data S1 and S2 for hierarchical cluster, gene names, IMAGE clone numbers and mean ratios) and were divided into groups based on their known function (Fig. 3). The majority of changes in gene expression were observed at 24 h. These included several genes previously implicated in HIV-1 infection, such as EGR1 (11), CDK4 (30), CD69 (3), CD71 (33, 37), CD87 (43), and JUN (6, 25), providing support for the validity of our other findings. Overall, genes involved in signaling and communication were stimulated, while genes involved in cell division, transcription, translation, and splicing, were inhibited by HIV-1BRU infection (Fig. 3). Many receptors and kinases were up-regulated, including CD69, T-cell receptor
, LDLR, mitogen-activated protein 2 kinase 1, and DYRK2. Positive cell cycle regulators were down-regulated and negative cell cycle regulators were up-regulated, reflecting G2 arrest of HIV-1 infected cells. Five members of the SR family of splicing factors were down-regulated, perhaps to facilitate the production of unspliced and single-spliced viral RNA messages needed at the later stages of the virus life cycle. Likewise, the expression of seven DEAD-box proteins (DDX, RNA helicases), several elongation initiation factors, and heterogeneous nuclear ribonucleoproteins was down-regulated. Though most transcription factors were down-regulated (e.g., MYC), a small group of transcription factors including EGR1, JUN, TNRC3, DRPLA, ARP, SMARCD2, SP140, GATA3, RELB, ELF1, and TAL1 were up-regulated by HIV-1 infection, implicating these factors in facilitating HIV-1 production. The majority of enzymes were down-regulated, including strong repression of MAT2A, ODC1, UNG, and ADSL. Interestingly, the genes encoding enzymes in the cholesterol biosynthesis pathway with sterol-responsive elements in their promoters (5), were all consistently up-regulated following HIV-1BRU infection (van 't Wout et al., submitted for publication).
![]() View larger version (66K): [in a new window] |
FIG. 3. Gene groups up- or down-regulated by HIV-1BRU infection of CCRF-CEM cells. (A) Up-regulated genes; (B) down-regulated genes. Genes were divided into seven broad categories based on their function (middle pie charts), with each category depicted by a separate pattern and numbers reflecting percentages of genes in that category compared to total number of genes regulated. Genes involved in cell signaling and communication were mostly up-regulated; this category was further divided into seven classes (pie charts at left). Genes involved in gene and protein expression were mostly down-regulated; this category was further divided into eight classes (pie charts at right).
|
0.005) in at least one cell line 24 h postinfection (see supplemental data S3 for common and anti-correlated responses between genes regulated in each cell line). Figure 4A shows the 319 genes that were twofold differentially expressed (P
0.005) in at least one cell line, illustrating both the common and cell-specific regulation observed. Responses in CEM, CEMss, and CEMfh cells were more similar to each other than to those in Jurkat or SupT1 cells, corresponding to the different origins of the various cell lines (see supplemental data S3 for correlation coefficients between the 24 h infections of each cell line). Of the 1,484 regulated genes, 71 genes were significantly differentially expressed in all eight experiments (46 down- and 25 up-regulated), shown in Fig. 4B and C (P
0.10).
![]() View larger version (67K): [in a new window] |
FIG. 4. Gene clusters impacted by HIV-1BRU infection of five CD4+-T-cell lines. Log10 ratios are depicted in green (down-regulated) or red (up-regulated). The magnitude of the regulation is illustrated by the intensity of the color. HI, heat-inactivated BRU versus BRU. (A) Cell-type-specific responses to HIV-1 infection. Three hundred nineteen genes exhibited >2-fold change in expression in at least one of the eight HIV-1-infected cell line experiments (P 0.005). (B and C) Homogeneous responses to HIV-1 infection. Seventy-one genes were regulated in all cell lines with greater than 90% confidence (P 0.10); 46 genes were down-regulated (B), and 25 genes were up-regulated (C).
|
![]() View larger version (24K): [in a new window] |
FIG. 5. Verification of RNA expression and comparison to changes in protein levels. Array results were compared to semiquantitative RT-PCR for eight genes (A) and for three membrane antigens by FACS (B). (A) Total RNA samples from a 24-h infection of CEMfh cells with HIV-1BRU or mock infection were tested with primers for 18S internal control and gene-specific primers for eight different genes at two different cycles. The graph depicts results obtained using microarrays on the x axis and those obtained by RT-PCR on the y axis. (B) Five CD4+-T-cell lines were infected with HIV-1BRU or mock infected with culture supernatant from uninfected cells. After 24 h, cells were washed and stained with specific purified monoclonal antibodies or isotype controls. The level of regulation corresponded with the level of p24 expression in each cell line (data not shown). The graph depicts results obtained using microarrays on the x axis and those obtained by FACS on the y axis.
|
Specificity of gene expression changes induced by HIV-1BRU infection.
To assess the specificity of the changes in gene expression induced by HIV-1 infection, we performed a compendium analysis, comparing genes impacted by HIV-1 and influenza A virus infection, as well as by interferon and heat shock treatment of cells. Figure 6 shows a two-dimensional cluster analysis of genes versus experiments depicting 271 genes that were twofold regulated (P
0.005) in at least 2 of the 34 experiments analyzed. In general, genes impacted by each of the four infections or treatments clustered separately, indicating that they were differentially regulated in response to these agents. However, a set of genes was coordinately regulated in the interferon and influenza virus experiments, reflecting the interferon response in influenza virus-infected cells. Comparing the heat shock, interferon, and influenza virus- and HIV-1-induced expression profiles, it is clear that genes regulated in one group were either not regulated or regulated in the opposite direction than in most of the other experiment groups. This indicates that HIV-1 infection largely induces changes specific to HIV-1 infection, rather than simply activation of stress-response or cytokine response pathways.
![]() View larger version (83K): [in a new window] |
FIG. 6. Impact of HIV infection compared to influenza A virus infection, heat shock, and interferon treatment. Two-dimensional cluster analysis of 271 genes in 34 array experiments comparing treatment of A549 cells with influenza A virus (16); HeLa, HFH, or HH2 cells with interferon (Geiss et al., unpublished observations); heat shock in CD4+-T-cell lines (van 't Wout et al., unpublished observations); and HIV-1 infection in CD4+-T-cell lines. Ratios were expressed as described in Fig. 4. Highlighted branches depict genes specifically regulated in HIV-1 infection experiments.
|
|
|
|---|
Transcription of several membrane-bound receptors was regulated by HIV-1 infection, notably CD53 and CD69 (both up-regulated) and CD71 (down-regulated). Down-regulation of the transferrin receptor (CD71) has been shown to involve Nef (33, 37) and may be related to observations of iron accumulation in patients (4, 38). Iron metabolism in HIV-1-infected cells could also be affected by the strong down-regulation of the transcription factor MYC we observed (49).
Interestingly, FACS analysis of membrane-bound receptor proteins revealed down-regulation of one of the three proteins tested (CD53) despite increases in the levels of RNA (Fig. 5). CD53 has been implicated in mediating CD2 signaling pathways and thus is involved in immune responses (1). Moreover, patients with genetic CD53 deficiency show increased susceptibility to viral infections (29), suggesting that HIV-1 might interfere with this potentially antiviral pathway. The viral protein Nef has been shown to down-regulate several membrane bound receptors, including CD4 and major histocompatibility complex molecules (31). CD53 may be another target for Nef-mediated down-regulation of surface expression because of its close association with CD4 (10). Thus, the increased mRNA expression levels of CD53 measured on arrays (and confirmed by RT-PCR) may reflect a cellular feedback loop sensing low surface expression.
Infected cells had decreased levels of proliferation or cycling, as cell cycle repressors were up-regulated and inducers were down-regulated. This could reflect the impact of the viral protein Vpr, which has been shown to induce G2 arrest (19, 21, 34, 36). Genes involved in gene and protein expression, including transcription factors, DNA helicases, and splice factors, were mostly down-regulated. A small group of transcription factors was up-regulated, suggesting a specific role in increasing HIV-1 production. Indeed, some of these factors have been previously implicated in HIV-1 infection (6, 11, 12, 25, 51).
Although many genes seem to be specific for the particular system used, a small group of genes has been consistently found to be regulated in most studies of the impact of HIV-1 infection on host cell gene expression. This group includes thymosin beta 10, prothymosin alpha, and T-cell receptor alpha chain (8, 15, 45). The majority of metabolic enzymes regulated in this study were decreased in HIV-1-infected cells. However, a group of enzymes known to be involved in cholesterol biosynthesis and regulated by the sterol responsive element binding protein 2 (SREBP2) were up-regulated (van 't Wout et al., submitted for publication) in all cell lines under study.
The changes induced by HIV-1 infection were most pronounced at 24 h postinfection. By this time, clear differences were observed between different CD4+-T-cell lines infected with HIV-1. Target cells have long been known to play an important role in HIV infection and pathogenesis, although the specific interactions and the importance of these different interactions to the disease process are unknown. Species barriers to infection are partly explained by the absence or mutation of host proteins essential for the viral life cycle in nonpermissive cells (44, 47). Similarly, certain permissive cells lack antiviral pathways that inhibit HIV-1 infection and can be infected with mutant strains lacking the viral protein Vif, thought to inhibit an antiviral pathway (27, 42). CD4+ T cells have been classified as permissive (e.g., CEMss and SupT1) and nonpermissive (e.g., CEM and PBMC) to infection with viruses lacking Vif (27, 42). Indeed, our gene expression studies found several groups of genes regulated by HIV-1 infection in only specific cell lines. Studying the differences between host-cell responses in different target cells provides an opportunity to identify host proteins essential for HIV-1 infection.
Gene expression profiling using microarrays is increasingly being used to study host-pathogen interactions. By determining the pattern of gene expression at different times, it is possible to elucidate which host genes are impacted, and how, over the course of infection. Identification of genes that are differentially regulated and the characterization of their functions provide a promising window on the understanding of pathogenicity and design of novel therapeutic agents.
This work was supported by grants from the American Foundation for AIDS Research (A.B.V.W.) and from the U.S. Public Health Service (J.I.M.).
Present address: Department of Internal Medicine, Division of Infectious Diseases, University of Texas Southwestern Medical Center, Dallas, TX 75390-9113. ![]()
Present address: Department of Pediatrics, University of Washington School of Medicine, Seattle, WA 91895-6320. ![]()
Present address: Rexagen Corporation, Seattle, WA 98103. ![]()
|
|
|---|
B activation: possible correlation with cell surface Tat-binding proteins. AIDS Res. Hum. Retrovir. 15:1209-1218.[CrossRef][Medline]
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»