Secondary Logo

Journal Logo

Basic and Translational Science

Antiretroviral Therapy Down-Regulates Innate Antiviral Response Genes in Patients With AIDS in Sub-Saharan Africa

Boulware, David R MD, MPH, DTM&H*†; Meya, David B MBChB, Mmed; Bergemann, Tracy L PhD§; Williams, Darlisha MPH*‖; Vlasova-St. Louis, Irina A MD, PhD†‖; Rhein, Josh MD*†; Staddon, Jack MD, PhD*†‖; Kambugu, Andrew MBChB, Mmed; Janoff, Edward N MD; Bohjanen, Paul R MD, PhD*†‖

Author Information
JAIDS Journal of Acquired Immune Deficiency Syndromes: December 1, 2010 - Volume 55 - Issue 4 - p 428-438
doi: 10.1097/QAI.0b013e3181ef4963

Abstract

INTRODUCTION

HIV infection elicits innate and adaptive host immune responses that contribute to control of the infection. However, the antiviral response is insufficient to prevent establishment of chronic infection, characterized by ongoing viral replication and viremia, destruction of lymphocytes, scarring of lymphatic tissues, and immune activation.1 Although antiretroviral therapy (ART) interrupts HIV-mediated immune destruction and increases the number and function of CD4+ T cells, the fundamental molecular mechanisms underlying HIV-associated immune compromise and subsequent ART-associated immune reconstitution remain largely unknown.

Microarray-based gene expression analyses have been used to identify pathologic pathways in cancer pathophysiology, to choose therapeutic targets, to predict treatment responses, and to identify novel pathways to categorize the pathophysiology of infectious diseases, including HIV/AIDS.2-6 Among 9 in vivo studies involving HIV-1 and microarrays between 2000 and 2006,7 all but one were cross-sectional, and the sole longitudinal study followed only 2 subjects.8 Subsequently, Li et al2,9 identified ∼200 ART-responsive genes in lymph nodes obtained by serial biopsies over 1 month in 5 HIV-infected individuals receiving ART and follow on experiments in 24 subjects. The ART-responsive genes included those encoding regulators of innate defense, cell trafficking, and tissue repair.2 More recently, Vahey et al3 examined gene expression changes in primates infected with simian immunodeficiency virus after 2 weeks of ART, and Jacquelin et al10 explored the longitudinal gene expression changes in African green monkeys. In 2010, Rotger et al6 reported cross-sectional genome-wide mRNA correlates of viral control. The initial studies of the effects of ART on gene expression in HIV-infected persons have been limited in size and duration, and none included longitudinal analyses in persons with AIDS.

We sought to use gene expression microarrays to characterize ART-induced changes in gene expression in peripheral blood of HIV-infected humans as a foundation to ultimately understand the role of ART in immune reconstitution. The use of peripheral blood in patients with successful responses to ART reveals the range and patterns of HIV-associated gene dysregulation and provides a normative framework for subsequent comparative studies in patients who fail therapy or develop ART-related complications, such as drug toxicity or immune reconstitution inflammatory syndrome (IRIS). This reverse challenge experiment can characterize the host response to HIV infection by quantifying the intraperson changes in gene expression occurring when HIV is blocked with successful ART.

We conducted gene expression microarray analysis of whole blood samples from a prospective cohort of 10 Ugandans with AIDS before and 2, 4, 8, and 24 weeks after initiation of ART. We identified approximately 200 transcripts whose expression was consistently changed after starting ART, approximately half of which are novel ART-responsive transcripts. Concomitant with the ART-associated drop in plasma HIV RNA, we observed consistent declines in the expression of genes regulating innate antiviral responses [eg, APOBEC3G, TRIM6, interferon (IFN)-signaling], pattern recognition receptors (DDX58), immune activation (LAMP3, NFκβ-pathway), cell apoptosis (IRF7, MX1, TREX1), cell proliferation, and numerous genes of unknown function. Our results suggest that HIV infection induces antiviral host responses, including innate antiviral responses, which are subsequently turned off when the virus is blocked with successful ART.

METHODS

Research Subjects and HIV Care

We prospectively monitored 24 ART-naive subjects with AIDS (<200 CD4+/μL), but no known active opportunistic infections, for one year after initiating ART. Exclusion criteria included anemia (hemoglobin < 8 g/dL), hepatitis (liver transaminases >2.5 × upper limit of normal), and renal insufficiency (creatinine >1.5 × upper limit of normal). Subjects were enrolled between May and July 2006 at the Infectious Disease Institute in Kampala, Uganda, and seen every 2 weeks for 12 weeks, then monthly through 52 weeks. At each clinic visit, subjects were screened for new infections, IRIS, malignancies, treatment failure, and ART compliance via pill counts. Antiretroviral regimens were stavudine, lamivudine, and nevirapine, or zidovudine, lamivudine, and efavirenz. Plasma HIV RNA and CD4+ T cells were measured every 12 weeks. Informed consent was obtained from subjects, and ethics approval was obtained from the University of Minnesota, Makerere University, and Uganda National Council of Science and Technology. Ten of 17 subjects with immunological recovery (CD4+ increase >50 cells/μL), viral suppression at 24 weeks (<400 copies/mL), and no new opportunistic infections or adverse events over 1 year of monitoring were chosen for microarray gene expression analyses.

Sample Collection and Gene Expression Analysis

We used Affymetrix microarrays for our gene expression studies performed on peripheral whole blood because this approach has been extensively validated. The correlation between Affymetrix microarrays and real-time polymerase chain reaction is very high (r = 0.87 to 0.92),11,12 and the reproducibility of microarrays actually exceeds that of polymerase chain reaction, even for common applications such as HIV-1 viral load quantification.13,14 We used peripheral whole blood for gene expression studies, rather than isolated peripheral blood mononuclear cells because in vitro artifacts can be induced during cellular isolation, ex vivo handling, and processing.15-18 We avoided these in vitro artifacts by drawing blood into PAXGene tubes, which immediately lyses cells and stabilizes RNA.15,18

RNA for gene expression analysis was isolated from peripheral whole blood at the time of commencing ART (baseline) and at 2, 4, 8, and 24 weeks of ART. Peripheral whole blood (2.5 mL) was collected into PaxGene tubes (PreAnalytiX, Hombrechtikon, Switzerland), incubated at room temperature for 2 hours, and frozen at −80°C. RNA was extracted using PreAnalytix PAXGene Blood RNA Kits and treated with Dnase I to remove DNA. Each subject's RNA was processed on the same day. The quality of extracted RNA was checked with a NanoDrop spectrophotometer (Thermo Fisher Scientific, Wilmington, DE). cDNA, prepared from 1 μg of RNA, was used to synthesize biotinylated-labeled RNA using the MessageAmp II aRNA Amplification Kit (Ambion, Austin, TX). Fragmented, labeled RNA (15 μg) was hybridized onto Human Genome U133_Plus2.0 arrays (Affymetrix, Santa Clara, CA) as specified by the manufacturer and scanned with an Affymetrix GeneChip 3000 Scanner. The U133_Plus2.0 array has 54,120 probe sets for 38,572 characterized genes.

Statistics

Feature intensities for each microarray were condensed into single intensity values for each probe set using Robust Multichip Average normalization. Robust Multichip Average performs background adjustment, quantile normalization, and summarization from the perfect-match probes.19 Postnormalization, we used a linear mixed effects model to test for changes in the expression series over time.

The regression tested all time points of 2, 4, 8, and 24 weeks versus baseline. The model was adjusted for the difference in log2CD4+ counts from baseline to 24 weeks (as a covariate) to assure that changes in expression were not related solely to CD4 change. The levels of other hematologic cell types (eg, total lymphocytes, CD8+ T cells, neutrophils, monocytes, eosinophils, basophils) did not change over time in the cohort (minimum P = 0.11). A person-specific random effect allowed correlation between time points within person. The vector of coefficients for time points derived from the regression model was tested for statistical significance. For each coefficient, we test whether the difference was significantly different from zero using the appropriate partial t test for fixed effects. In addition to tests for each time point, we also performed an overall test for a change at all time points.

Because there were a large number of tests (4 time points + 1 overall test × 54,675 probe sets = 273,375 statistical tests in total), P values need to be adjusted for multiple comparisons. We employed the Benjamini-Hochberg correction to control the false discovery rate for the 273,375 tests.20 Of the 273,375 P values, 66,212 were statistically significant with a range of adjusted P values from 1.7 × 10−6 to 0.05 (unadjusted P values of I.1 × 10−10 to 0.007). All presented P values were adjusted for multiple comparisons, unless otherwise stated.

We then filtered the list of transcripts based on (1) consistent direction of expression change (up or down-regulated) at all time points, (2) genes which were statistically significant at all 4, 8, 24 weeks (unadjusted P < 0.05) time points, and (3) statistically significant overall in the regression model (adjusted P < 0.05).

Validation

To validate changes in gene transcripts, serial serum samples collected at baseline, 2,4,8,12, and 24 weeks were analyzed for changes in the levels of 27 cytokines/chemokines (Bio-Rad, Hercules, CA) per the manufacturer's protocol via a Luminex 100 system (Austin, TX). The cytokine levels were Iog2 transformed and longitudinal analysis was performed mirroring the statistical analysis of gene expression using the R nlme package. We modeled changes in cytokines by linear regression versus weeks of ART with and without a covariate of absolute CD4 count and individual slope estimates for each time point after baseline. An intraperson random effects variable was utilized for intraperson comparison of repeated measures.

Identification of Gene Networks and Functional Pathways

Ingenuity Pathways Analysis (www.ingenuity.com) was used to visualize gene regulation networks, gene-gene interactions, functional pathways, and biological processes influenced by ART. The Ingenuity Pathways Knowledge Base is the largest curated database of previously published findings on mammalian biology from the public literature. Ontology results were cross validated with the DAVID public database.21

To identify transcripts related to T-cell activation, we performed a comparative analysis between transcripts that were down-regulated with ART and a previously published gene expression dataset of human T-lymphocytes stimulated in vitro with anti-CD3+/anti-CD28+ using Affymetrix U95 and U133A microarrays.22

RESULTS

We studied 10 Ugandan subjects with AIDS, randomly selected from among 17 subjects without new opportunistic infections and with appropriate immunologic recovery and viral suppression by 24 weeks of ART (see Supplemental Digital Content 1, part 1A, https://links.lww.com/QAI/A97). Clinical characteristics, including the low baseline CD4 counts (median 30 cells/μL; mean 72 cells/μL), are summarized in Table 1. All subjects achieved virologic suppression (<400 copies/mL) by 24 weeks of ART, and the median CD4 count increased to 102 cells per microliter (IQR: 91-173 cells/μL).

T1-3
TABLE 1:
HIV-Related Parameters in 10 Subjects Selected for Gene Expression Analysis

We characterized ART-induced changes in gene expression after 4, 8, or 24 weeks of ART. The complete dataset is posted at http://bohjanenlab.umn.edu/. We found that more genes were down-regulated than up-regulated, and that the majority of genes were either not expressed or not affected by ART (see Supplemental Digital Content 1, part B, https://links.lww.com/QAI/A97). We found 239 probesets with consistently changed expression at each individual time point (unadjusted P < 0.05) and over the entire time course based on a linear regression model (adjusted P < 0.05, adjusting for 273,375 multiple comparisons). These 239 probesets represented 208 unique genes, of which 160 were down-regulated and 48 were up-regulated (see Supplemental Digital Content 2, https://links.lww.com/QAI/A98). Even after adjusting for increases in CD4+ T-cell numbers over time, the list remained unchanged. Examples of the down-regulated transcripts are shown in Table 2. Most of the down-regulated genes showed progressively decreased expression over time during ART (Fig. 1, see Supplemental Digital Content 3, https://links.lww.com/QAI/A99).

T2-3
TABLE 2:
Downregulation of 53 Genes at all Time Points With HIV Therapy
F1-3
FIGURE 1:
ART-associated down-regulation of gene expression, A, This heat map presents the average fold changes in gene expression from baseline of the 53 genes which were consistently down-regulated with ART. These genes were discovered by linear regression (adjusted P < 0.001) and further filtered by (1) being down-regulated in direction at all time points after starting ART and (2) being testing for statistical significance for the change from baseline at individual 8 and 24 week time points by separate paired t tests (adjusted P < 0.05). Decreased gene expression is represented by the intensity of blue color. To the right of the panel, the functions and pathways are indicated for genes with known relationship to: interferon signaling, apoptosis, immune response, viral immune response, and T-cell activation. Changes by individual subject are provided online (see Supplemental Digital Content 3, https://links.lww.com/QAI/A99) and representative examples of individual genes displayed online as line graphs (see Supplemental Digital Content 5, https://links.lww.com/QAI/A101); B, Displays the top canonical (mechanistic) pathways down-regulated at 8-weeks of ART. The significance of the association between the data set and the canonical pathway was measured in 2 ways: First, the Fischer exact test was used to calculate a P value determining the probability that the association between the genes in the dataset and the canonical pathway is explained by chance alone (Blue bars, primary y axis). Second, a ratio of the number of genes from the data set that map to the pathway divided by the total number of genes that map to the canonical pathway is displayed as orange line (secondary y axis).

Biological Processes

Among the 208 genes with differential regulation, 137 were mapped to a known biological process by gene ontology (108 down-regulated; 29 up-regulated) and 71 had unknown biological process. The known functions that were down-regulated included response to virus (n = 16, P = 10−12), response to other (foreign) organisms (n = 17, P = 10−9), responses to biotic stimulus (n = 20, P = 10−9), multiorganism processes (n = 17, P = 10−7), responses to stimulus (n = 49, P = 10−6), immune response (n = 24, P = 10−6), immune system process (n = 26, P = 10−4), and defense response (n = 16, P = .007) (Table 3). Among the 160 consistently down-regulated genes, 29 (14%) were known to have HIV protein interactions. Up-regulated genes did not map to consistent biologic processes, had largely unknown functions, and none had known HTV-protein interactions.

T3-3
TABLE 3:
Functional Categories of Genes With Altered Expression After HIV Therapy in Persons With AIDS

The annotated functions of the down-regulated genes suggest that many are involved in immune responses (both innate and antiviral), cellular proliferation, and/or immune activation. The specific down-regulated transcripts encoded several important components of innate antiviral immune responses, including host restriction factors such as APOBEC3A, TRIM6, and TRIM38 (P < 3 × 104) (see Supplemental Digital Content 2, https://links.lww.com/QAI/A98)23,24 and chemokines such as CXCL9 and CXCL10. Paradoxically, transcripts that encode mediators of apoptosis (XAF1) and cell proliferation pathways were down-regulated during ART, highlighting that HIV pathogenesis is characterized by cellular activation and proliferation and cell death. Notably, 33% (52 of 160) of these down-regulated genes currently do not have known gene biologic ontology biological process assignments, and 31% (49 of 160) have unknown molecular function, demonstrating that numerous genes of unknown function are up-regulated by HIV infection.

Interferons play prominent roles in innate and antiviral responses. Indeed, annotated pathway analyses revealed that transcripts encoding multiple components of the canonical (mechanistic) pathway for type-I interferon signaling pathways were coordinately down-regulated at all points with ARV-associated viral suppression (Fig. lB) (Fig. 2) (P < 0.001) Changes in the expression of transcripts encoding components of other pathways involved with innate control of viral replication were seen at later time points. For example, transcripts encoding components of the protein ubiquitination pathway (eg, TRIM5, CUL5, and UBR2) and APOBEC3G were also coordinately down-regulated, but only after 8 weeks. In this context, decreased protein ubiquitination/degradation likely reflects decreased cell death after virologic suppression.25 Overall, we identified numerous novel ART-responsive genes that encode important components of the host response to HIV infection. That these changes occurred in the context of ART-induced suppression of HIV replication suggests that these processes were up-regulated upon exposure to HIV and subsequently down-regulated when HIV was blocked with ART. Thus, ART serves as a “reverse challenge” to identify HTV-induced responses.

F2-3
FIGURE 2:
Coordinate down-regulation of genes involved in interferon signaling and innate immune responses. The biological network of genes within the type-I interferon-signaling pathway is shown. Genes colored in blue were down-regulated with ART, with the intensity of blue reflecting the degree of down-regulation.

Relationship With T-Cell Activation

The majority (71%) of the genes down-regulated in response to ART are regulated during T-cell activation. Of the 208 transcripts influenced by ART, the expression of 154 transcripts has been studied over time after ex vivo stimulation of primary human T-cells with anti-CD3 and anti-CD28 antibodies using an older generation Affymetrix microarray.14 Of the 154 transcripts, 110 showed alteration in expression compared with resting/unstimulated T cells (P < 0.05) and 84 transcripts exhibited very early changes in gene expression within 30 minutes of stimulation (62 up-regulated, 22 down-regulated with a fold change of >1.5) (see Supplemental Digital Content 2, https://links.lww.com/QAI/A98). These results suggest that many of the ART-responsive genes are regulated during T-cell activation.

Protein Levels of ART-Responsive Cytokines

For some of the ART-responsive genes encoding secreted proteins, we measured protein levels in serum to determine if changes in RNA expression correlated with changes in protein expression in order to provide validation for our results. We serially measured protein levels of 28 serum cytokines in the same 10 patients. Serum protein levels for all of the cytokines tested are provided (see Supplemental Digital Content 4, https://links.lww.com/QAI/A100). Consistent with the down-regulation in gene expression, we identified decreasing serum levels of the chemokines CXCL9 (coefficient estimates at week 8 = −1.348, week 24 = −2.752, P < 0.0001) and CXCL10 (coefficient estimates at week 8 = −2.327, week 24 = −2.891, P < 0.0001) over the first 24 weeks of ART. Thus, serum levels of CXCL10 (also called IP-10) decreased by 80% (95% CI: 33% to 118%), whereas gene expression decreased by 101%. Levels of IL-12, IL-1ra, IL-6, TNF-α, and TRAIL showed neither a statistically significant decline with ART by microarray nor by serum protein concentrations over the first 24 weeks of ART (P > 0.2). Overall, our data showed that ART-associated changes in protein levels of serum cytokines and chemokines correlated well with changes in RNA expression.

DISCUSSION

AIDS is caused by profound T-cell loss and immunosuppression in the context of dramatic chronic immune activation. HIV infection activates antiviral immune responses that combat the virus, whereas at the same time, the virus promotes cell death through direct killing of infected cells and by indirect apoptotic mechanisms. We propose that the genes that were up-regulated in response to the virus during uncontrolled viremia are subsequently down-regulated when viral replication is blocked with ART. Thus, evaluating the changes in gene expression in peripheral blood in response to ART can provide insight into HIV pathogenesis by identifying gene expression pathways that were turned on in response to HIV infection.

We performed microarray gene expression analysis over time in 10 patients with AIDS from Sub-Saharan Africa who were successfully treated with ART. This study represents the largest human prospective longitudinal characterization of ART-related gene expression to date. We identified 160 genes down-regulated after ART, including 87 that are novel HIV/ART-responsive genes, not previously reported in in vitro or in vivo studies.2,4,6,9,27-29,38 The major functions of the ART-responsive genes include (1) innate anti-viral immune responses, including interferon responses, (2) immune activation and cellular proliferation/apoptosis, and (3) unknown function. The diversity of ART-responsive genes reveals the breadth of the host responses to HIV and the ability of ART to reverse these effects.

Many of the novel down-regulated ART-responsive genes that we identified encode regulators of innate immune responses, including innate antiviral responses. For example, transcripts encoding APOBEC3A and APOBEC3G were down-regulated after ART. APOBEC3G, which inhibits HIV replication through its cytidine deaminase activity,39 was down-regulated only at the 8-week time point, whereas the related APOBEC3A, which has no known anti-HIV activity,40 was down-regulated at all time points. HIV has developed mechanisms to counteract effective antiviral restriction by APOBEC proteins. For example, the HIV Vif protein binds directly to APOBEC3G and inhibits its function by targeting APOBEC3G for degradation via the ubiquitin-proteasome pathway.41-43 Interestingly, we found expression of genes encoding components of this pathway,44 including CUL5, RBX1, SKP1, and UBR2 were down-regulated after 8 weeks of ART. Because the ubiquitin-proteasome pathway promotes the degradation of APOBEC3G by Vif,41-43,45 the induction of these genes could represent an antiviral response. Our results suggest that genes controlling APOBEC and the protein degradation pathway that targets Vif are coordinately induced in response to HIV infection, because the expression of these genes was abrogated after ART. These responses likely occur as a host response to try to limit viral replication.

Expression of the TRIM5 restriction factor gene was also down-regulated at the 8-week time point. In many species, TRIM5 prevents retroviral infection by binding to retroviruses and recruits the proteasome to degrade viral proteins.23,24,46 Although human TRIM5 has only weak activity against HIV, TRIM5 is effective at preventing other retroviruses from infecting human cells.23,24,47,48 Another restriction factor down-regulated was TRIM22. In the context of HIV, TRIM22 is induced by interferons to disrupt viral replication by preventing the HIV structural protein Gag from trafficking to the cell surface.49 TRIM22 had high pre-ART levels of expression in peripheral blood (>99th percentile), and was down-regulated at 4 weeks, 8 weeks, and overall. The induction of viral restriction genes such as TRIM genes after HIV infection represents additional innate antiviral responses that attempt to limit viral replication, and these responses are subsequently down-regulated with ART.

Numerous other genes encoding components of innate responses were also down-regulated. For example, 16 down-regulated genes included pattern recognition receptors. ADAR, a gene up-regulated in response to poliovirus dsRNA and simian immunodeficiency virus in African green monkeys,10,50 which may modulate the posttranscriptional regulation of HIV-1 env gene expression,51 was down-regulated. Additional PRR genes involved in innate responses to dsRNA that were down-regulated included DDX58 (also known as RIG-1) and thereafter IRF7.5,9,10,52 Although DDX58 has not been previously associated with HIV infection, DDX58 is involved in antiviral innate immunity.53,54 We also found the transcript encoding IFIH1, a pattern recognition receptor functioning as an IFN-induced helicase, was down-regulated. IFIH1 has been associated with increased HIV viral mRNA expression in vitro,33 suggesting that HIV may usurp this antiviral response to promote viral replication. Because viral replication is virtually stopped by effective ART, understandably, many of these HIV-induced antiviral and innate defense pathways are turned off with ART-mediated viral suppression. These PRR genes are most likely induced as part of a host warning system triggered by HIV to allow other foreign material to be rapidly recognized.

Genes encoding multiple components of interferon-response pathways were also down-regulated after ART, suggesting that these pathways were induced by HIV infection and are components of the innate antiviral response to HIV (Fig. 3). As circulating HIV virions decline after ART, decreased expression of type-1 interferon signaling genes, including IFI27, ISG15, ISG20 and IFI44L, occurs, likely due to lack of antigenic challenge once the viral stimulus is eliminated. ISG20 is an IFN-induced3′→5′ exonuclease specific for ssRNA that is up-regulated with HIV-infection in vitro.55ISG15 has known antiviral properties in the late stage of HIV-1 assembly by inhibiting the ubiquitination of HIV Gag during endosomal trafficking.56,57 These interferon signaling genes are known to be activated in B-lymphocytes in response to HIV viremia.4 Once CD4+ cells are infected by HIV, interferon responses induce apoptosis of infected cells as an antiviral response facilitating the elimination of infected cells and preventing both viral replication and infection of new cells.55 The upregulation of interferon and interferon response pathways may help limit viral replication, but the associated apoptosis probably contributes to immune destruction and immunosuppression.

F3-3
FIGURE 3:
Innate antiviral immune responses down-regulated with HIV therapy. Known innate anti-viral immune response genes which are down-regulated with ART are shown in blue and HIV proteins are shown in green. After viral entry, TRIM5 binds to the viral capsule and inhibits infection, although it only a weak inhibitor of HIV infection. DDX58 primarily binds to dsRNA viruses but may also bind to dsDNA in the cytoplasm to activate downstream IRF7. Post HIV RNA transcription, the HIV protein Vif (virion infectivity factor) binds to APOBEC3G targeting it for destruction by the ubiquitin-proteasome pathway via known interactions with CUL5 and RBX1 and recognition of the E3-ligase complex by UBR2. Additional partial innate defenses include TRIM22 which blocks the trafficking of HIV Gag protein and ISG16 which interferes with viral budding. RT, reverse transcription.

Because HIV pathogenesis is characterized by a state of chronic immune activation, it is not surprising that numerous transcripts related to immune activation were prominently down-regulated with ART. Of 160 down-regulated genes, 110 are T-cell activation genes that we have shown previously to be up-regulated after ex vivo stimulation of primary human T-cells with anti-CD3 and anti-CD28 antibodies.22 Several transcripts encoding well-known markers of cellular activation, including LAMPS (CD63), were down-regulated over time after initiation of ART,58,59 demonstrating that ART produces an overall decrease in immune activation. We also found that immune activation genes involved in cellular proliferation/apoptosis pathways such as the NF-κβ pathway and XAF1 were down-regulated. NF-κβ is a multi-subunit transcription factor that regulates the expression of numerous genes involved in immune activation, including interleukin-2 and TNF-a.60 NF-κβ acts as a transcription factor activating HIV transcription and facilitating viral replication.44 In addition, NF-κβ is involved in stress responses, cell death, and virus-mediated apoptosis.60 Another down-regulated apoptosis regulator was XAF1, an IFN-stimulated gene that mediates type-1 IFN-induced sensitization to TRAIL, thereafter causing apoptosis.61 Chronic immune activation and the associated activation-induced apoptosis is a major contributor to immune cellular depletion in HIV-infected persons. Our results suggest that although peripheral CD4+ T-cell counts increase after ART, the overall expression of proliferation-related and apoptosis-related transcripts in peripheral blood decreases. Thus, the contributions of cellular proliferation, activation and apoptosis to promoting viral replication and causing immunodeficiency are abrogated during ART.

A large percentage of ART-responsive genes identified (33%) have unknown biological processes and functions. The finding that ART causes down-regulation of these genes suggests that they had been induced by HIV and may therefore play roles in the host response to HIV. It is possible, however, that the expression of some of these genes may have changed as a direct consequence of exposure of human cells to antiretroviral agents, rather than as a consequence of ART-mediated viral suppression. These genes with unknown function may be targets for future studies to explore their function and role in HIV pathogenesis. A large number of genes of the genes we identified to be down-regulated after ART (and presumably up-regulated during HIV infection) are newly described in this context or have unknown function, highlighting that we are only beginning to unravel the complex nature of immune regulation during HIV infection.

Prior recognition of ART-responsive genes in human subjects has been limited and our results significantly extend our knowledge about the influence of ART on gene expression, whereas some of our results confirm results from previous studies. For example, of our 160 genes that were down-regulated with ART, 43 overlap with genes that up-regulated in CD4+ cells in association with increasing HIV viral load,6 suggesting that our hypothesis that these genes that were down-regulated in response to ART, had previously been up-regulated in response the HIV.

Our results also confirm the previous report from Li et al2 showing that IFN-related genes, such as OAS, STAT, IRF7, CCR1, CCR5, CXCL9, CXCL10, and MX1 genes, were down-regulated in lymph node biopsy specimens after ART, and these results are strikingly similar to our study. In our prospective, longitudinal study, however, we identified many heretofore unrecognized ART-responsive genes in humans. This greater range of response was revealed in part because we utilized the current Affymetrix microarray which probed the entire 38,572-gene genome, whereas, Li, et al,2 examined the expression of only 9838 genes with the earlier generation U95A microarray. Indeed, 35% of ART-responsive genes we identified were not represented on the earlier microarray. In contrast to our results, Li et al2 also found down-regulation of the transcripts encoding the cytokines IFN-γ, MIP-1β, and TRAIL in lymphatic tissue of persons with less advanced HIV; findings we did not observe in peripheral blood mRNA or serum cytokines. We sampled blood from 10 persons with AIDS in Sub-Saharan Africa (median CD4+ 30 cells/μL) over multiple standardized time points for 24 weeks, whereas, Li et al2 tested lymphatic tissues from 5 US subjects with less advanced HIV infection (mean CD4+ of 415 ± 226/μL) and evaluated changes at 1 month only. The different results likely derive from differences in sample type, stage of infection, and timing of ART.

In conclusion, we characterized the host response to HIV through a reverse challenge approach by examining gene expression changes associated with ART-induced viral suppression in persons with AIDS. We identified numerous transcripts encoding important components of innate antiviral pathways, including antiviral restriction factors and pattern recognition receptors, which were not previously known to be regulated by HIV infection. Our results confirmed existing paradigms suggesting that chronic immune activation, cellular proliferation, and apoptosis play roles in HIV pathogenesis, and we extended knowledge in the field by defining gene expression pathways and components that regulate these responses during HIV infection. Our results also provide information about the normal response to ART in AIDS patients that may be useful for identifying abnormal responses to ART in persons who develop adverse events such as IRIS, drug reactions, or new AIDS-related opportunistic infections. For example, AIDS patients with these adverse events on ART might not show the characteristic gene expression response to ART, but rather may express abnormal gene expression signatures of immune activation or inflammation that could be used to identify or diagnose adverse events.

ACKNOWLEDGMENTS

The authors thank Ashley Haase, Quinsheng Li, and Anne Marie Weber-Main for critically reading this article.

REFERENCES

1. Schacker TW, Nguyen PL, Beilman GJ, et al. Collagen deposition in HIV-1 infected lymphatic tissues and T cell homeostasis. J Clin Invest. 2002;110:1133-1139.
2. Li Q, Schacker T, Carlis J, et al. Functional genomic analysis of the response of HIV-1-infected lymphatic tissue to antiretroviral therapy. J Infect Dis. 2004;189:572-582.
3. Vahey MT, Ockenhouse CF, Wang Z, et al. Impact of antiretroviral treatment on gene expression in peripheral blood mononuclear cells from SIVmac251-infected macaques. J Infect Dis. 2007;196:384-393.
4. Moir S, Malaspina A, Pickeral OK, et al. Decreased survival of B cells of HIV-viremic patients mediated by altered expression of receptors of the TNF superfamily. J Exp Med. 2004;200:587-599.
5. Montano M, Rarick M, Sebastiani P, et al. Gene-expression profiling of HIV-1 infection and perinatal transmission in Botswana. Genes Immun. 2006;7:298-309.
6. Rotger M, Dang KK, Fellay J, et al. Genome-wide mRNA expression correlates of viral control in CD4+ T-cells from HIV-1-infected individuals. PLoS Pathog. 2010;6:el000781.
7. Giri MS, Nebozhyn M, Showe L, et al. Microarray data on gene modulation by HIV-1 in immune cells: 2000-2006. J Leukoc Biol. 2006;80:1031-1043.
8. Guadalupe M, Reay E, Sankaran S, et al. Severe CD4+ T-cell depletion in gut lymphoid tissue during primary human immunodeficiency virus type 1 infection and substantial delay in restoration following highly active antiretroviral therapy. J Virol. 2003;77:11708-11717.
9. Li Q, Smith AJ, Schacker TW, et al. Microarray analysis of lymphatic tissue reveals stage-specific, gene expression signatures in HIV-1 infection. J Immunol. 2009;183:1975-1982.
10. Jacquelin B, Mayau V, Targat B, et al. Nonpathogenic SIV infection of African green monkeys induces a strong but rapidly controlled type IIFN response. J Clin Invest. 2009;119:3544-3555.
11. Consortium M, Shi L, Reid LH, et al. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006;24:1151-1161.
12. Chen JJ, Hsueh HM, Delongchamp RR, et al. Reproducibility of microarray data: a further analysis of microarray quality control (MAQC) data. BMC Bioinformatics. 2007;8:412.
13. Spanish Society of Infectious Diseases and Clinical Microbiology, Orta Mira N, Guna Serrano Mdel R, Latorre Martínez JC, et al. Analysis of the results of the SEIMC External Quality Control Program for HIV-1 and HCV viral loads, 2007 [in Spanish]. Enferm Infecc Microbiol Clin. 2008;26(Suppl 13):8-13.
14. Brambilla D, Jennings C, Aldrovandi G, et al. Multicenter evaluation of use of dried blood and plasma spot specimens in quantitative assays for human immunodeficiency virus RNA: measurement, precision, and RNA stability. J Clin Microbiol. 2003;41:1888-1893.
15. Baechler EC, Batliwalla FM, Karypis G, et al. Expression levels for many genes in human peripheral blood cells are highly sensitive to ex vivo incubation. Genes Immun. 2004;5:347-353.
16. Debey S, Schoenbeck U, Hellmich M, et al. Comparison of different isolation techniques prior gene expression profiling of blood derived cells: impact on physiological responses, on overall expression and the role of different cell types. Pharmacogenomics J. 2004;4:193-207.
17. Feezor RJ, Baker HV, Mindrinos M, et al. Whole blood and leukocyte RNA isolation for gene expression analyses. Physiol Genomics. 2004;19:247-254.
18. Chai V, Vassilakos A, Lee Y, et al. Optimization of the PAXgene blood RNA extraction system for gene expression analysis of clinical samples. J Clin Lab Anal. 2005;19:182-188.
19. Bolstad BM, Irizarry RA, Astrand M, et al. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185-193.
20. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289-300.
21. Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44-57.
22. Raghavan A, Dhalla M, Bakheet T, et al. Patterns of coordinate down-regulation of ARE-containing transcripts following immune cell activation. Genomics. 2004;84:1002-1013.
23. Huthoff H, Towers GJ. Restriction of retroviral replication by APOBEC3 G/F and TRIM5alpha. Trends Microbiol. 2008;16:612-619.
24. Langelier CR, Sandrin V, Eckert DM, et al. Biochemical characterization of a recombinant TRIM5alpha protein that restricts human immunodeficiency virus type 1 replication. J Virol. 2008;82:11682-11694.
25. Glickman MH, Ciechanover A. The ubiquitin-proteasome proteolytic pathway: destruction for the sake of construction. Physiol Rev. 2002;82:373-428.
26. Imbeault M, Ouellet M, Tremblay MJ. Microarray study reveals that HIV-1 induces rapid type-I interferon-dependent p53 mRNA up-regulation in human primary CD4+ T cells. Retrovirology. 2009;6:5.
27. Rempel H, Calosing C, Sun B, et al. Sialoadhesin expressed on IFN-induced monocytes binds HIV-1 and enhances infectivity. PLoS ONE. 2008;3:el967.
28. van der Kuyl AC, van den Burg R, Zorgdrager F, et al. Sialoadhesin (CD169) expression in CD14+ cells is upregulated early after HIV-1 infection and increases during disease progression. PLoS ONE. 2007;2:e257.
29. Silverman RH, Sengupta DN. Translational regulation by HIV leader RNA, TAT, and interferon-inducible enzymes. J Exp Pathol. 1990;5:69-77.
30. Vollbrecht T, Brackmann H, Henrich N, et al. Impact of changes in antigen level on CD38/PD-1 co-expression on HlV-specific CD8 T cells in chronic, untreated HIV-1 infection. J Med Virol. 2010;82:358-370.
31. Coetzee LM, Tay SS, Lawrie D, et al. From research tool to routine test: CD38 monitoring in HIV patients. Cytometry B Clin Cytom. 2009;76:375-384.
32. Tuaillon E, Al Tabaa Y, Baillat V, et al. Close association of CD8+/CD38 bright with HIV-1 replication and complex relationship with CD4+ T-cell count. Cytometry B Clin Cytom. 2009;76B:249-260.
33. Cocude C, Truong MJ, Billaut-Mulot O, et al. A novel cellular RNA helicase, RH116, differentially regulates cell growth, programmed cell death and human immunodeficiency virus type 1 replication. J Gen Virol. 2003;84(Pt 12):3215-3225.
34. Kalia V, Sarkar S, Gupta P, et al. Rational site-directed mutations of the LLP-1 and LLP-2 lentivirus lytic peptide domains in the intracytoplasmic tail of human immunodeficiency virus type 1 gp41 indicate common functions in cell-cell fusion but distinct roles in virion envelope incorporation. J Virol. 2003;77:3634-3646.
35. Kreuzer KA, Dayer JM, Rockstroh JK, et al. The IL-1 system in HIV infection: peripheral concentrations of IL-1 beta, IL-1 receptor antagonist and soluble IL-1 receptor type II. Clin Exp Immunol. 1997;109:54-58.
36. Yadav A, Fitzgerald P, Sajadi MM, et al. Increased expression of suppressor of cytokine signaling-1 (SOCS-1): a mechanism for dysregulated T helper-1 responses in HIV-1 disease. Virology. 2009;385:126-133.
37. Yan N, Cherepanov P, Daigle JE, et al. The SET complex acts as a barrier to autointegration of HIV-1. PLoS Pathog. 2009;5:el000327.
    38. Zauli G, Milani D, Mirandola P, et al. HIV-1 Tat protein down-regulates CREB transcription factor expression in PC12 neuronal cells through a phosphatidylinositol 3-kinase/AKT/cyclic nucleoside phosphodiesterase pathway. FASEB J. 2001;15:483-491.
    39. Mangeat B, Turelli P, Caron G, et al. Broad antiretroviral defence by human APOBEC3G through lethal editing of nascent reverse transcripts. Nature. 2003;424:99-103.
    40. Goila-Gaur R, Khan M, Miyagi E, et al. Targeting APOBEC3A to the viral nucleoprotein complex confers antiviral activity. Retrovirology. 2007;4:61.
    41. Yu X, Yu Y, Liu B, et al. Induction of APOBEC3G Ubiquitination and Degradation by an HIV-1 Vif-Cul5-SCF Complex. Science. 2003;302:1056-1060.
    42. Mehle A, Goncalves J, Santa-Marta M, et al. Phosphorylation of a novel SOCS-box regulates assembly of the HIV-1 Vif-Cul5 complex that promotes APOBEC3G degradation. Genes Dev. 2004;18:2861-2866.
    43. Xiao Z, Ehrlich E, Yu Y, et al. Assembly of HIV-1 Vif-Cul5 E3 ubiquitin ligase through a novel zinc-binding domain-stabilized hydrophobic interface in Vif. Virology. 2006;349:290-299.
    44. Ockenhouse CF, Bernstein WB, Wang Z, et al. Functional genomic relationships in HIV-1 disease revealed by gene-expression profiling of primary human peripheral blood mononuclear cells. J Infect Dis. 2005;191:2064-2074.
    45. Zhang W, Chen G, Niewiadomska AM, et al. Distinct determinants in HIV-1 Vif and human APOBEC3 proteins are required for the suppression of diverse host anti-viral proteins. PLoS ONE. 2008;3:e3963.
    46. Sawyer SL, Emerman M, Malik HS. Discordant evolution of the adjacent antiretro viral genes TRIM22 and TRIM5 in mammals. PLoS Pathog. 2007;3:el97.
    47. Kaiser SM, Malik HS, Emerman M. Restriction of an extinct retrovirus by the human TRIM5alpha antiviral protein. Science. 2007;316:1756-1758.
    48. Stremlau M, Owens CM, Perron MJ, et al. The cytoplasmic body component TRIM5alpha restricts HIV-1 infection in Old World monkeys. Nature. 2004;427:848-853.
    49. Barr SD, Smiley JR, Bushman FD. The Interferon Response Inhibits HIV Particle Production by Induction of TRIM22. PLoS Pathog. 2008;4:el000007.
    50. Morrissey LM, Kirkegaard K. Regulation of a double-stranded RNA modification activity in human cells. Mol Cell Biol. 1991;11:3719-3725.
    51. Phuphuakrat A, Kraiwong R, Boonarkart C, et al. Double-stranded RNA adenosine deaminases enhance expression of human immunodeficiency virus type 1 proteins. J Virol. 2008;82:10864-10872.
    52. Harenberg A, Guillaume F, Ryan EJ, et al. Gene profiling analysis of ALVAC infected human monocyte derived dendritic cells. Vaccine. 2008;26:5004-5013.
    53. Gack MU, Shin YC, Joo CH, et al. TRIM25 RING-finger E3 ubiquitin ligase is essential for RIG-I-mediated antiviral activity. Nature. 2007;446:916-920.
    54. Biacchesi S, LeBerre M, Lamoureux A, et al. MAVS plays a major role in induction of the fish innate immune response against RNA and DNA viruses. J. Virol. 2009:83;7815-7827.
    55. Espert L, Degols G, Lin YL, et al. Interferon-induced exonuclease ISG20 exhibits an antiviral activity against human immunodeficiency virus type 1. J Gen Virol. 2005;86(Pt 8):2221-2229.
    56. Okumura A, Lu G, Pitha-Rowe I, et al. Innate antiviral response targets HIV-1 release by the induction of ubiquitin-like protein ISG15. Proc Natl Acad Sci USA. 2006;103:1440-1445.
    57. Harty RN, Pitha PM, Okumura A. Antiviral activity of innate immune protein ISG15. J Innate Immun. 2009;l:397-404.
    58. Chen H, Dziuba N, Friedrich B, et al. A critical role for CD63 in HIV replication and infection of macrophages and cell lines. Virology. 2008;379:191-196.
    59. Kaufmann GR, Zaunders JJ, Cunningham P, et al. Rapid restoration of CD4 T cell subsets in subjects receiving antiretroviral therapy during primary HIV-1 infection. AIDS. 2000;14:2643-2651.
    60. Bren GD, Trushin SA, Whitman J, et al. HIV gpl20 induces, NF-kappaB dependent, HIV replication that requires procaspase 8. PLoS ONE. 2009;4:e4875.
    61. Micali OC, Cheung HH, Plenchette S, et al. Silencing of the XAF1 gene by promoter hypermethylation in cancer cells and reactivation to TRAIL-sensitization by IFN-beta. BMC Cancer. 2007;7:52.
    Keywords:

    adult; gene expression profiling/BL; HAART; HIV; AIDS; oligonucleotide array sequence analysis

    Supplemental Digital Content

    © 2010 Lippincott Williams & Wilkins, Inc.