- Research article
- Open Access
- Open Peer Review
Genome-wide search for the genes accountable for the induced resistance to HIV-1 infection in activated CD4+ T cells: apparent transcriptional signatures, co-expression networks and possible cellular processes
BMC Medical Genomicsvolume 6, Article number: 15 (2013)
Upon co-stimulation with CD3/CD28 antibodies, activated CD4 + T cells were found to lose their susceptibility to HIV-1 infection, exhibiting an induced resistant phenotype. This rather unexpected phenomenon has been repeatedly confirmed but the underlying cell and molecular mechanisms are still unknown.
We first replicated the reported system using the specified Dynal beads with PHA/IL-2-stimulated and un-stimulated cells as controls. Genome-wide expression and analysis were then performed by using Agilent whole genome microarrays and established bioinformatics tools.
We showed that following CD3/CD28 co-stimulation, a homogeneous population emerged with uniform expression of activation markers CD25 and CD69 as well as a memory marker CD45RO at high levels. These cells differentially expressed 7,824 genes when compared with the controls on microarrays. Series-Cluster analysis identified 6 distinct expression profiles containing 1,345 genes as the representative signatures in the permissive and resistant cells. Of them, 245 (101 potentially permissive and 144 potentially resistant) were significant in gene ontology categories related to immune response, cell adhesion and metabolism. Co-expression networks analysis identified 137 “key regulatory” genes (84 potentially permissive and 53 potentially resistant), holding hub positions in the gene interactions. By mapping these genes on KEGG pathways, the predominance of actin cytoskeleton functions, proteasomes, and cell cycle arrest in induced resistance emerged. We also revealed an entire set of previously unreported novel genes for further mining and functional validation.
This initial microarray study will stimulate renewed interest in exploring this system and open new avenues for research into HIV-1 susceptibility and its reversal in target cells, serving as a foundation for the development of novel therapeutic and clinical treatments.
A comprehensive picture of the host factors putatively supporting HIV-1 (human immunodeficiency virus, type 1) replication in cells has emerged from recent siRNA studies [1–4] and meta-analysis . Genome-wide landscapes of host genes and proteins involved in HIV-1 infection and disease progression have also been established in gene array studies [6–9] and novel proteomic approaches [10–12]. In parallel, genome-wide association studies (GWAS) have revealed a set of inheritable genetic variations in large populations related to susceptibility to HIV-1 infection [13–15].
Against these exciting developments, there is still no global view of the host cellular factors that render the target cells resistance to infection, in spite of a few well-studied restriction factors [16–19]. A recent genome-wide screening for novel restriction factors  further highlighted this awareness and our interest in establishing a more holistic pictures of the host determinants working against HIV-1 susceptibility.
A crucial starting point for studies of this type is the consideration of sample sources, cell types, and experimental settings. In this regard, CD4 + T cells are the first choice; they are the major cell type amongst all the susceptible targets and reservoirs of HIV-1 infection . Apart from their intrinsic susceptibility, CD4 + T cells’ activation in vivo (during the natural courses of HIV-1 infection) and in vitro (typically with PHA/IL-2 stimulation) is generally recognized as an absolute prerequisite for the virus to replicate productively . However, rather unexpectedly, Levine et al. (1996) found that activation by co-stimulation with CD28 led to a complete loss of susceptibility to HIV-1 infection in these cells. This phenomenon was subsequently confirmed by several independent groups [23–26]. In spite of the fact that several further studies attempted to explore the potential of using thus-stimulated cells for the treatment of SIV/HIV infection in monkey models and clinical settings, progress in elucidating the underlying cellular and molecular mechanisms seems to have halted since 2002 for unknown reasons.
In the present study, in order to investigate the possible mechanisms of the observed reversal of HIV-1 susceptibility in these activated CD4 + T cells, we first replicated the reported experimental settings  and then performed genome-wide expression analysis using Agilent microarrays.
Isolation and stimulation of CD4+T cells
Buffy coats were obtained from healthy donors (Kunming Blood Station) and peripheral blood mononuclear cells (PBMCs) were isolated by Ficool-Hypaque (TBD Sciences) gradient centrifugation. Resting CD4 + T cells were then purified by magnetic negative selection (CD4 + T cell Isolation Kit II, Miltenyi Biotec) and only the aliquots of cells with purity > 95% as determined by flow cytometry were used for further analysis. This study was reviewed and approved by the internal review board of the Kunming Institute of Zoology, Chinese Academy of Sciences (approval ID: RTYX20090910-1, approval date: 2009-09-10). All donors provided written informed consent for participation in this study.
Cells were stimulated as previously reported [27, 28]. Briefly, freshly isolated CD4 + T cells were resuspended in RPMI 1640 medium (Gibco) supplemented with 10% heat inactivated fetal bovine serum (Gibco) and 20 mM HEPES (Amresco) and seeded in 6-well plates either at an initial density of 2 × 106 cells/well with human recombinant IL-2 (100 U/ml) and PHA (5 μg/ml) or 0.5 × 106 cells/well with polystyrene beads coated anti-CD3/CD28 antibodies (Dynal beads CD3/CD28 T Cell Expander, Dynal) at a bead to cell ratio of 3:1. Half media were changed every 2 days in CD3/CD28 costimualted CD4 + T cells and every 3 days in PHA/IL-2 stimulated cells. Cells were then cultured at 37°C in a humidified incubator with 5% CO2 for 6 days. Cells and the derived RNA samples and data sets were labeled “P”, “R”, and “B”, respectively, according to the PHA/IL-2 stimulated, un-stimulated resting and beads-stimulated settings.
Aliquots of cells to be analyzed were washed with FACS buffer (PBS supplemented with 1% bovine serum albumin (BSA)) and stained for 40 minutes in the dark at 4°C with anti-CD45RO-FITC, anti-CD25-PE, anti-CXCR4, fluorescent dye 5-(and −6)-carboxyfluorescein diacetate succinimidyl ester (CFSE) (Sigma-Aldrich), anti-Ki67(Abcam), anti-CCR5 (Biolegend) and anti-CD69-PE-Cy-5 (BD Biosciences) and appropriate isotype controls. After washing with FACS (Fluorescence Activated Cell Sorter) buffer, cells were fixed with 4% paraformaldehyde and analyzed on a FACS Calibur using Cell Quest and FlowJo 7.6.1.
Gene expression profiling with Agilent microarrays
Total RNA was extracted from stimulated or un-stimulated CD4 + T cells by TRIzol (Invitrogen) followed by a purification using RNeasy columns (Qiagen) according to the manufacturer’s protocols. The amount and quality of RNA preparations were evaluated on an Agilent 2100 Bioanalyzer with RNA6000 Nano Reagents and Supplies (Agilent).
Quality-checked RNAs were then transcribed with the First-Strand cDNA Synthesis Kit (Agilent) and their expression data obtained using Agilent 4 × 44 K Human Whole-Genome 60-mer oligonucleotide microarrays according to the protocols by the manufacturer. The original microarray data from this study are available at the NCBI GEO database (http://www.ncbi.nlm.nih.gov/geo/) under the accession number GSE34252.
Analysis of microarray data
The normalized ratio of the gene expression signals was log2 transformed and hierarchical clustering was performed with average linkage. The clustered heatmap was visualized using Treeview. The RVM (Random variance model) f-test was applied to filter differentially expressed genes for the different situations. After the significance analysis and FDR (false discovery rate) analysis, we selected the differentially expressed genes according to the p-value and FDR threshold set at p < 0.01 and FDR <0.01 [29–31]. And the fold changes of any two groups are more than 2.
Series-Cluster analysis was performed to identify the global trends and model profiles of expression according to signal density under the “P”, “R” and “B” conditions and in the P-R-B sequence. Fisher’s exact test and the multiple comparison test [32, 33] were applied to identify the model profiles with probability significantly higher than expected as random.
Gene ontology (GO) analysis  was performed to facilitate elucidating the biological implications of unique genes in the significant or representative profiles. GO analysis was used to find the main function of the genes having the same expression trend according to the Gene Ontology. Fisher’s exact test and χ 2 test were applied to identify the significant GO categories and FDR was used to correct the p-values.
Gene co-expression networks analysis  were performed to track the interactions among the differentially expressed genes, according to their normalized signal intensity in 6 representative profiles. Pearson correlation was applied to each pair of genes and the significantly correlated pairs were used to construct the network . To locate the core regulatory genes in the networks, k-core scoring was introduced to simplify graph topology analysis [37, 38]. A k-core of a given gene indicates its hub or nodal status with connection to “k” other genes in a network [37, 38]. Accordingly, the genes with largest k-core scores and highest degrees of connection were identified as “key regulatory genes” in a network  and those unique to each network were selected as “marker genes” for the cellular status of intrinsic susceptibility or induced resistance to HIV-1 infection.
Cellular characterization of stimulated and un-stimulated CD4+ T cells
Following the previously published details of the initial experiment (Levine et al. 1996) and other independent studies [23–26], we re-established de novo and verified the experimental systems (Figure 1). After stimulation for 6 days with the CD3/CD28 antibodies coated beads (labeled “B” for beads), purified CD4 + T cells became enlarged and highly proliferative, forming a large amount of big cell colonies (Figure 1A). When compared with the un-stimulated (labeled “R” for resting) or PHA/IL-2-stimulated (labeled “P” for PHA) cells at the functional level, the induced resistance to HIV-1 infection in these co-stimulated cells [23–28] was also readily replicated in our earlier experiments (described as a brief report in a home journal ). We monitored the cell viability and proliferation first by Trypan blue exclusion test which showed more than 95% viable cells for all the experiments (data not shown) and then by Ki67 and CFSE staining (Additional file 1: Figure S1 and Additional file 2: Figure S2). We then proceeded to analyze the surface phenotypes of the generated cells and found that a homogeneous population of cells emerged following the co-stimulation (Figure 1B). In contrast to the controls, these cells were uniformly and simultaneously stained positive for both the general and early activation markers (CD25 and CD69, respectively) at high levels (Figure 1B). Another important marker for activated and memory T cells, CD45RO, was also uniformly expressed following co-stimulation in these “B” cells but in a much lower percentage than in “R” or “P” cells (Figure 1B).
Global expression profiles and signatures
Using Agilent 4 × 44 K Human Whole-Genome Microarrays, we found that the overall gene expression patterns were clearly and sharply different among the “R”, “P” and “B” cells. As shown in Figure 2 (left) with their clustered heatmaps and hierarchical patterns, in total 7,824 genes represented 8,128 transcripts (out of the 41,000 test probes on the chips) were differentially expressed according to the RVM (Random variance model) algorithm (p-value < 0.05, FDR < 0.05) (Additional file 3: Table S1).
In subsequent Series-Cluster analysis, we identified 16 possible profiles (Figure 2 middle) (Additional file 4: Table S2), which represent the overall expression patterns. Of these, 10 (#1, 5, 6, 7, 8, 9, 10, 11, 12 and 15, Figure 2 middle) showed the sharp difference between “B” and “P” cells. More specifically, genes down-regulated in “P” cells are shown in the “Pdown” (#1) profile and those up-regulated in these cells were grouped in the “Pup” (#9) profile (Figure 2 right). Similarly, genes up- or down-regulated only in “B” are represented in “Bup” (#8) and “Bdown” (#15) profiles. We merged three profiles (#5, 6 and 7) that contain the genes down-regulated in “P” but up-regulated in “B” into a new “PdownBup” profile, in order to show their special importance in the reduced resistance as they contain the genes which were down regulated in the susceptible “P” cells and at the same up-regulated in the resistant “B” cells. In the same way, we merged profiles #10, 11 and 12 into a new “PupBdown” profile to show the genes that were up-regulated in “P” cells but at the same time down-regulated in the “B” cells (Figure 2 right).
These 6 representative profiles contain 1,345 differentially expressed genes in “B” and “P” cells (Figure 3 left) and give the apparent expression signatures either unique to “P” and “B” cells or polarized expressed in “P” and “B” cells at two opposite extremities. The remaining 6,479 genes differentially expressed in the other 6 profiles (#2, 3, 4, 13, 14 and 16, Figure 2 middle) were unique to the “R” cells, up- or down-regulated when compared with “P” and “B” cells as the baseline. Due to the nature of the “R” cells that are seemingly susceptive to the infection but do not support active viral replication, the biological importance of these genes in HIV-1 infection are more complicated to dissect and interpret. We kept aside these 6,479 “R” cell-unique genes for future studies.
Functional categories and significant ontologies of the differentially expressed genes in “P” and “B” cells
It is immediately noticeable (Figure 3 left) that most of the 1,345 differentially expressed genes in the 6 representative profiles falls in the “B” cell-related expression, with 471 genes in the “Bdown” profile and 465 in “Bup”, respectively, representing the down-regulated and up-regulated un-permissive genes in the resistant “B” cells (Figure 3 left, Additional file 5: Table S3). In the “Bdown” profile, 222 genes with known functions gave a rather complicated ontology category but, after applying Fisher’s exact test, χ 2 test, and FDR analysis, a clearer picture emerged. Statistically, 43 genes were significantly involved in 4 major functional categories (Figure 3 middle), showing: (a) response to stimulus (immune response, cell surface receptor linked signal transduction and positive regulation of natural killer cell mediated cytotoxicity); (b) metabolic process; (c) multicellular organismal development; and (d) cellular process are predominantly influenced by the co-stimulation.
Other part of exactly the same 4 major categories was found significantly influenced in the “Pup” and “PupBdown” profiles, with part of one more category (localization) being also up-regulated in the permissive “P” cells or down-regulated in the un-permissive “B” cells (Figure 3 middle). Combined GO analysis of the genes in the “Pup”, “PupBdown” and Bdown” profiles showed one more category (cell proliferation) was also significantly influenced (Figure 3 right). Altogether, these 6 significantly affected functional categories contain 101 genes (Figure 3 right) representing the whole set of the “permissive genes” that are associated positively with the highly susceptible status or negatively with the induced resistance of the cells. As shown in Figure 4A, these permissive genes are clearly hierarchical and highly enriched in the categories related to defense, immune/inflammatory response and signal transduction under the general heading of response to stimuli .
More genes (352) with known functions were identified in the “Bup” profile with 117 being significantly involved in 7 major GO categories, showing that (a) cellular process (cell adhesion, cell differentiation and regulation of transcription), (b) multicellular organismal development, (c) metabolic process, (d) response to stimulus, (e) localization, (f) circulatory system process, and (g) others are predominantly up-regulated following the co-stimulations. Exactly the same 7 categories were also up-regulated in the “PdownBup” profile, although the numbers of the genes in each category varied (Figure 3 middle right). These 7 categories contain all the 4 significant ones in the “Pdown” profile and one more category (cell fate determination) was also identified. Altogether, the 8 categories contain 144 genes representing the whole set of the “resistant genes” that are associated positively with the induced resistance or negatively with the highly susceptible status of the cells (Figure 3 right). However, they are dispersed in more diverse categories (Figure 4B) than the permissive ones (Figure 4A), indicating that the resistant genes have wider influence in cellular functions and may hold the crucial check points for HIV-1 infection in “multicellular organismal development”, “cellular process” and “metabolic process”. One of the striking examples is that 30 genes involved in the process of cell adhesion were found up-regulated in “B” cells (profiles “Bup” and “PdownBup”) (Figure 4B) (Additional file 6: Table S4).
Direct interactions of the genes with known functions with HIV-1 proteins
When subjected to further screening, 22 of the 144 potential “resistant genes” and 23 of the 101 potential “permissive genes” were identified for their direct interactions with HIV-1 proteins (Figure 5A & Additional file 7: Table S5), according to the data in the “HIV-1 Human Protein Interaction Network” (http://www.ncbi.nlm.nih.gov/RefSeq/HIVInteractions/). Their significant functional categories are shown in Figure 5B.
Specifically, 3 (Env, Nef, and Tat) of HIV-1 proteins attracted most of the host (no matter permissive or resistant) factors in this set of data, suggesting that host factors may have the capacity to counterbalance with the permissive genes by interacting with the same early genes in HIV-1 infection. This is coincided with the initial report that the novel anti-HIV effect induced by co-stimulation happens around viral entry and before integration . However, in terms of host gene functions, no overlap was identified between the resistant and permissive genes as the former were mapped on the categories of cell adhesion, multicellular organismal development and cell differentiation, while the latter on immune response, cell surface receptor signal transduction, cellular defense response and chemotaxis. Fewer but prominent genes, e.g., CD59, IL-2, IL-15, INDO and FN1 were found to interact with Gag and Vpr, although we did not find any genes in all our profiles capable of interacting with HIV-1 Vif or Rev.
Co-expressed genes and their networks in the permissive “P” and un-permissive “B” cells
Given that HIV-1 tends to interact with “key” host proteins, such as bottlenecks and hubs in gene interaction networks [42, 43], we used the genes in the representative profiles and constructed 6 gene co-expression networks (Figure 6 left). We then applied the “k-core” scores (see Methods) to identify those with highest networking degrees as the “key regulatory” genes that may play pivotal roles in gene interactions and regulations. Thus, 53 genes in profiles “Pup”, “Bdown” and “PupBdown”, which tentatively contain the “permissive genes” were chosen as “key regulatory” genes and, among them, 29 genes with known functions (underlined in Figure 6 middle) are involved in immune response, signal transduction and so on (Figure 6 right). Similarly, 84 genes in profiles “Pdown”, “Bup” and “PdownBup”, which tentatively contain the “resistant genes”, stood out with the highest k-core scores and, out of them, 36 genes (underlined in Figure 6 middle) are involved in cell adhesion, multicellular organismal development and so on (Figure 6 right).
KEGG pathways and schematic overview of HIV-host gene interactions in the resistant “B” cells
We then mapped all the genes with known functions in the significant categories (Figure 3) and all the “key regulatory” genes (Figure 6) on KEGG pathways from the Kyoto Encyclopedia of Genes and Genomes database. A clear picture emerged in the resistant “B” cells showing (Figure 7) that (a) the up-regulated genes in “Bup” and “PdownBup” profiles were enriched in the cellular processes (Figure 8) including: filamentous actin (F-actin), tight junction, actomyosin assembly contraction, proteasomes, proteolysis, lysosomes, degradation, and Na+ Ca2+ exchange, and (b) the down-regulated genes in “Bdown” and “PupBdown” profiles were enriched in the process of actin polymeration, apoptosis, cell cycle checkpoint, ER to Golgi transport and terminally misfolded.
Verification of the differential gene expression by FACS
We then performed FACS analysis to examine the expression levels of two HIV-1 co-receptors, CCR5 and CXCR4. As shown in (Figure 9), CCR5 was up-regulated in the permissive “P” cells but down-regulated in the resistant “B” cells, while CXCR4 remained down-regulated in both “P” and “B” cells. These data not only validated our microarray analysis—albeit only a small part of the whole data set—but also were consistent with the initial findings that the induced resistance was only effective to CCR5-dependent strains of the HIV-1 virus [24, 44, 45]. This also confirmed our previous report  showing that, following co-stimulation, the levels of CCR5 expression varied to some degrees in the 12 individuals tested, but the general trend was clearly down-regulated. Conceivably, the down-regulation of the critical co-receptor CCR5 accounts, at least in part, for the induced resistance in the co-stimulated cells.
Clues for novels genes accounting for HIV-1 susceptibility and resistance
It should be noted that the aforementioned co-expression network analysis was performed using all the genes with and without known functions. Strikingly, among the total 137 “key regulatory” genes, only 9 genes have been previously reported to interact with HIV-1 (Table 1) while the remaining majority (128 genes) as of yet have no confirmation of their potentials to interact with HIV-1 genes or proteins. When going back to the 1,345 differentially expressed genes in the 6 representative profiles (Figure 2 left), around a third of them (513 genes) do not yet have any assigned functions. In view of their prominence in the representative profiles and co-expression networks, these unknown genes constitute a new source for future studies on their roles in biological processes and potential involvement in determining HIV-1 susceptibility and resistance.
In the present study, we re-established an experimental system for cellular and molecular studies on the induced resistance to HIV-1 infection in CD4 + T cells as first reported by Levine et al. . This cellular system includes two methods most commonly used to activate and expand CD4 + T cells in vitro. Incubation with PHA/IL-2 induces proliferating CD4 + T cells which are highly susceptible to HIV-1 infection and highly permissive for the subsequent viral replication, whereas co-stimulation with CD3/CD28 reverses the intrinsic susceptibility in these cells and renders them with resistance to HIV-1 infection and un-permissiveness for viral replication. Thus, these two extremely polarized statuses of activated CD4 + T cells, with the un-stimulated resting cells (susceptible to the infection but not permissive for the rapid viral replication) as control in the middle of the “spectrum of permissiveness”, constitute a highly insightful system in the search for host factors responsible for HIV-1 susceptibility and its reversal in CD4 + T cells.
We were attracted more by the biological significance and potentials of the induced resistance in the “B” cells. We first replicated this in vitro system and verified its usefulness (Figure 1 and ) at the cellular level. Considering the controversy around the source cells, i.e., which subsets (naïve or memory CD4 + T cells) could be induced to generate the resistance [24, 26, 46], we adopted a strategy using total CD4 + T cell population to start our study without prior cell sorting before co-stimulations. This gave us the chance to address the issue of which subset of CD4 + T cells was generated to be responsible for the induced resistance and led us to the first set of our findings, which shows a homogenous cell population emerged from the co-stimulation. These cells expressed three activation markers (CD25, CD69 and CD45RO) at high levels at the same time. Although the expression of these classical activation markers may not ultimately account for the induced resistance, they served well as indicators for the highly activated population and, more relevant to the following analysis, the homogeneity of the induced cells gave us the confidence for the subsequent whole-genome wide search for the “key regulatory genes” and the “core marker genes” in the status switch to resistance.
Apart from the above novel findings in defining the cell population emerged during the reversal of susceptibility, the present study is the first whole-genome-wide analysis for the genes that account for the induced resistance. Once again, the system allowed us to pin down the potential permissive and resistant genes since susceptible “P” and the resistant “B” cells are polarized at the two extremities of the “permissive spectrum”. As a result, the 1,345 genes in the 6 specific profiles for the “B” and “P” cells could be truly the representatives accounting for the intrinsic susceptibility to HIV-1 infection and its induced reversal in the target cells. This simple dichotomy classification further allowed us to analyze the differentially expressed genes for their gene ontology and co-expression with ease and confidence.
We have not included the popular pathway and transcript factor analysis in this initial genome-wide search but, by focusing on the gene categories and co-expression, managed to identify the major cellular processes involved the reversed susceptibility in the major HIV-1 target cells. Prominently, these include actin cytoskeleton system, protein degradation and cell cycle arrest. First, several lines of evidence highlighted that actin cytoskeleton were regarded as a barrier and the hijack of actin cytoskeleton facilitates entry of HIV into its target cells [47, 48]. In our study, 8 genes (ACTN2, PTPRF, SORBS1, PAR3, PVRL2, MYL9, MYH10 and PAK2) were found up-regulated in the resistant “B” cells (“Bup” profile in Figure 8A). These genes mediate assembly and contraction of actin and actomyosin and formation of tight junctions. Whether they formed a barrier against HIV-1 infection warrant further functional verifications. To the contrary, genes associated with actin polymerization (FN1 and ITGAX, Figure 8A), which have been shown to promote viral binding and entry  and is also necessary in chemotaxis and cytokinesis  was sharply down regulated in the “B” cells. Given that chemotaxis may actually serve to fuel the infection response by recruiting susceptible, activated CD4 + T cells to the virus, ultimately aiding viral dissemination , we therefore hypothesized that co-stimulation could also block the tracks for HIV-1 release and dissemination by reducing actin polymerization and chemotactic response.
Second, genes involved in ubiquitin-dependent degradation by the proteasome (WNT3, FZD1, FZD2, FZD3, FBXW11, RAD23B, DBAHV2, HLPA4L, IRAK3, IL1RAP, PROC and PROS1), which is the third step of endoplasmic reticulum associated protein degradation, were also predominantly in “B” cells (Figure 8B). This is consistent with previously report that the ubiquitin-proteasome was negatively associated with HIV-1 replication by acting to destroy incoming viral complexes at the early steps [50, 51]. We also found that ATPVOA4 and IGF2R, which are important for lysosome-mediated degradation of endocytosed proteins [52, 53], were up-regulated in “B” cells (Figure 8B). Therefore, co-stimulation could enhance the cellular defense to HIV-1 infection by up-regulating these lysosome- and proteasome-related genes .
Third, it is also noticeable that 8 unique genes (Table 2) involved in cell cycle arrest were up-regulated in “B” cells. Apparently, these genes may have helped for the co-stimulated cells to reach the highly proliferative status. However, their role in HIV-1/AIDS infection and disease progression need further detailed studies since modifications to the cell cycles was observed in both long-term non-progressors (LTNPs) and progressive patients [55, 56]. Anyhow, this line of evidence reflects that host target cells may have intrinsic mechanisms, which, when activated, would halt viral infection at the specific point of the cell cycle.
The other two major cellular programs which were profoundly influenced by the co-stimulation were metabolism and apoptosis (Figure 4B). This seems inconsistent with previous studies showing that (a) numerous metabolism-associated genes were down-regulated in LTNPs  and (b) apoptosis is widely accepted as a mechanism for T cell depletion in vitro and in vivo[3, 57, 58]. Further 5 genes (CA2, MED28, PLOD2, SLC2A3 and TMTC1), which were up-regulated as probably resistant genes in our “B” cells, were proposed as “host dependent factors (HDF)” in recent small-interfering RNA-knockdown screens [1–4, 20]. More surprisingly perhaps, in our set of potentially resistant genes, we did not find any overlap with those reported to be associated with LTNPs , viral latency [59, 60], or even the “resistance genes” identified by virus-host interaction network analysis using public data sets  (data not shown). However, these apparent discrepancies do not necessarily discredit our studies but, to the contrary, may favorably underscore the uniqueness and importance of the highly controlled experimental system we adopted and call for further studies down the line.
Finally, we chose to examine the expression of CCR5 and CXCR4, the 2 major HIV-1 co-receptors, on the cell surface by FACS analysis as the first part of our ongoing efforts to verify our array data and analysis. To our satisfaction, their expression in this set of experiments and our previous screening  consistently corroborated our bioinformatics analysis and helped in explaining the fact that the induced resistance by the co-stimulation is a CCR5-dependent phenomenon [24, 44, 45]. Obviously, large-scale validation and functional studies at both mRNA and protein levels are needed to corroborate the role of the key regulatory genes (e.g., ITGAX, VEGFA, FN1, CCND1, CA12, LAMB3, MYL9, etc.) and their pathways as the major determinants for the induced resistance. We hope the present study serves as a fresh call for a renewed interest in the experimental system as a useful model for the search of novel host resistant factors.
In summary, we replicated a simple yet powerful cellular system with CD3/CD28 co-stimulation and confirmed its usefulness in studying the induced resistance to HIV-1 in CD4 + T cells. This initial microarray study, although still descriptive and correlative in nature, allowed us a chance to glean valuable new insights into this phenomenon. Based on the overall expression patterns and apparent signatures of the differentially expressed genes, we managed to pin down 245 (101 potentially permissive and 144 potentially resistant) significant functional genes and 137 (84 potentially permissive and 53 potentially resistant) “key regulatory” genes involved in the reversal of target cell susceptibility to HIV-1 infection. Moreover, we showed that this system has the potential as a rich source for the search of novel genes accounting for the intrinsic susceptibility to HIV-1 and its reversal in the major target cells. We hope our findings will stimulate renewed interest in investigating the mechanisms underlying this phenomenon. We surmise that future studies of this kind will help in furthering our understanding of HIV-1 infection, which may ultimately then lead to the development of novel biomarkers and therapeutics.
Brass AL, Dykxhoorn DM, Benita Y, Yan N, Engelman A, Xavier RJ, Lieberman J, Elledge SJ: Identification of host proteins required for HIV infection through a functional genomic screen. Science. 2008, 319 (5865): 921-926. 10.1126/science.1152725.
Konig R, Zhou YY, Elleder D, Diamond TL, Bonamy GMC, Irelan JT, Chiang CY, Tu BP, De Jesus PD, Lilley CE: Global analysis of host-pathogen interactions that regulate early-stage HIV-1 replication. Cell. 2008, 135 (1): 49-60. 10.1016/j.cell.2008.07.032.
Zhou HL, Xu M, Huang Q, Gates AT, Zhang XHD, Castle JC, Stec E, Ferrer M, Strulovici B, Hazuda DJ: Genome-Scale RNAi Screen for Host Factors Required for HIV Replication. Cell Host Microbe. 2008, 4 (5): 495-504. 10.1016/j.chom.2008.10.004.
Yeung ML, Houzet L, Yedavalli VSRK, Jeang KT: A genome-wide short hairpin RNA screening of jurkat T-cells for human proteins contributing to productive HIV-1 replication. J Biol Chem. 2009, 284 (29): 19463-19473. 10.1074/jbc.M109.010033.
Bushman FD, Malani N, Fernandes J, D’Orso I, Cagney G, Diamond TL, Zhou HL, Hazuda DJ, Espeseth AS, Konig R: Host Cell Factors in HIV Replication: Meta-Analysis of Genome-Wide Studies. PLoS Pathog. 2009, 5 (5): e1000437-10.1371/journal.ppat.1000437.
Mehla R, Ayyavoo V: Gene Array Studies in HIV-1 Infection. Curr HIV/AIDS Rep. 2012, 9: 34-43. 10.1007/s11904-011-0100-x.
Giri MS, Nebozhyn M, Showe L, Montaner LJ: Microarray data on gene modulation by HIV-1 in immune cells: 2000–2006. J Leukoc Biol. 2006, 80 (5): 1031-1043. 10.1189/jlb.0306157.
Li QS, Smith AJ, Schacker TW, Carlis JV, Duan L, Reilly CS, Haase AT: Microarray Analysis of Lymphatic Tissue Reveals Stage-Specific, Gene Expression Signatures in HIV-1 Infection. J Immunol. 2009, 183 (3): 1975-1982. 10.4049/jimmunol.0803222.
Smith AJ, Li QS, Wietgrefe SW, Schacker TW, Reilly CS, Haase AT: Host Genes Associated with HIV-1 Replication in Lymphatic Tissue. J Immunol. 2010, 185 (9): 5417-5424. 10.4049/jimmunol.1002197.
Jager S, Cimermancic P, Gulbahce N, Johnson JR, McGovern KE, Clarke SC, Shales M, Mercenne G, Pache L, Li K: Global landscape of HIV-human protein complexes. Nature. 2012, 481 (7381): 365-370.
Zhang L, Zhang X, Ma Q, Zhou H: Host Proteome Research in HIV Infection. Genomics Proteomics Bioinformatics. 2010, 8 (1): 1-9. 10.1016/S1672-0229(10)60001-0.
Kontijevskis A, Wikberg JES, Komorowski J: Computational proteomics analysis of HIV‒1 protease interactome. Proteins Struct Function Bioinformatics. 2007, 68 (1): 305-312. 10.1002/prot.21415.
Limou S, Le Clerc S, Coulonges C, Carpentier W, Dina C, Delaneau O, Labib T, Taing L, Sladek R, Deveau C: Genomewide Association Study of an AIDS-Nonprogression Cohort Emphasizes the Role Played by HLA Genes (ANRS Genomewide Association Study 02). J Infect Dis. 2009, 199 (3): 419-426. 10.1086/596067.
Le Clerc S, Limou S, Coulonges C, Carpentier W, Dina C, Taing L, Delaneau O, Labib T, Sladek R, Deveau C: Genomewide Association Study of a Rapid Progression Cohort Identifies New Susceptibility Alleles for AIDS (ANRS Genomewide Association Study 03). J Infect Dis. 2009, 200 (8): 1194-1201. 10.1086/605892.
Guergnon J, Theodorou I: What did we learn on host’s genetics by studying large cohorts of HIV-1-infected patients in the genome-wide association era?. Curr Opin HIV AIDS. 2011, 6 (4): 290-296. 10.1097/COH.0b013e3283478449.
Sheehy AM, Gaddis NC, Choi JD, Malim MH: Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein. Nature. 2002, 418 (6898): 646-650. 10.1038/nature00939.
Stremlau M, Owens CM, Perron MJ, Kiessling M, Autissier P, Sodroski J: The cytoplasmic body component TRIM5α restricts HIV-1 infection in Old World monkeys. Nature. 2004, 427 (6977): 848-853. 10.1038/nature02343.
Neil SJD, Zang T, Bieniasz PD: Tetherin inhibits retrovirus release and is antagonized by HIV-1 Vpu. Nature. 2008, 451 (7177): 425-430. 10.1038/nature06553.
Tyagi M, Kashanchi F: New and novel intrinsic host repressive factors against HIV-1: PAF1 complex, HERC5 and others. Retrovirology. 2012, 9 (1): 19-10.1186/1742-4690-9-19.
Liu L, Oliveira NMM, Cheney KM, Pade C, Dreja H, Bergin AMH, Borgdorff V, Beach DH, Bishop CL, Dittmar MT: A whole genome screen for HIV restriction factors. Retrovirology. 2011, 8: 94-10.1186/1742-4690-8-94.
Levy JA: Pathogenesis of human-immunodeficiency-virus infection. Microbiol Rev. 1993, 57 (1): 183-289.
Ho M, Armstrong J, McMahon D, Pazin G, Huang XL, Rinaldo C, Whiteside T, Tripoli C, Levine G, Moody D: A phase-1 study of adoptive transfer of autologous CD8+ T-lymphocytes in patients with acquired-immunodeficiency-syndrome (AIDS)-related complex or AIDS. Blood. 1993, 81 (8): 2093-2101.
Roederer M, Raju PA, Mitra DK, Herzenberg LA: HIV does not replicate in naive CD4 T cells stimulated with CD3/CD28. J Clin Invest. 1997, 99 (7): 1555-10.1172/JCI119318.
Mengozzi M, Malipatlolla M, De Rosa SC, Herzenberg LA, Roederer M: Naive CD4 T cells inhibit CD28-costimulated R5 HIV replication in memory CD4 T cells. Proc Natl Acad Sci. 2001, 98 (20): 11644-10.1073/pnas.211205098.
Barker E, Bossart KN, Levy JA: Differential effects of CD28 costimulation on HIV production by CD4+ cells. J Immunol. 1998, 161 (11): 6223.
Spina CA, Prince HE, Richman DD: Preferential replication of HIV-1 in the CD45RO memory cell subset of primary CD4 lymphocytes in vitro. J Clin Invest. 1997, 99 (7): 1774-1785. 10.1172/JCI119342.
Levine BL, Mosca JD, Riley JL, Carroll RG, Vahey MT, Jagodzinski LL, Wagner KF, Mayers DL, Burke DS, Weislow OS: Antiviral effect and ex vivo CD4+ T cell proliferation in HIV-positive patients as a result of CD28 costimulation. Science. 1996, 272 (5270): 1939-10.1126/science.272.5270.1939.
Carroll RG, Riley JL, Levine BL, Feng Y, Kaushal S, Ritchey DW, Bernstein W, Weislow OS, Brown CR, Berger EA: Differential regulation of HIV-1 fusion cofactor expression by CD28 costimulation of CD4+ T cells. Science. 1997, 276 (5310): 273-276. 10.1126/science.276.5310.273.
Wright GW, Simon RM: A random variance model for detection of differential gene expression in small microarray experiments. Bioinformatics. 2003, 19 (18): 2448-10.1093/bioinformatics/btg345.
Yang H, Crawford N, Lukes L, Finney R, Lancaster M, Hunter KW: Metastasis predictive signature profiles pre-exist in normal tissues. Clin Exper Meta. 2005, 22 (7): 593-603. 10.1007/s10585-005-6244-6.
Clarke R, Ressom HW, Wang A, Xuan J, Liu MC, Gehan EA, Wang Y: The properties of high-dimensional data spaces: implications for exploring gene and protein expression data. Nat Rev Cancer. 2008, 8 (1): 37-49. 10.1038/nrc2294.
Ramoni MF, Sebastiani P, Kohane IS: Cluster analysis of gene expression dynamics. Proc Natl Acad Sci. 2002, 99 (14): 9121-10.1073/pnas.132656399.
Miller LD, Long PM, Wong L, Mukherjee S, McShane LM, Liu ET: Optimal gene expression analysis by microarrays. Cancer Cell. 2002, 2 (5): 353-361. 10.1016/S1535-6108(02)00181-2.
Harris MA, Clark JI, Ireland A, Lomax J, Ashburner M, Collins R, Eilbeck K, Lewis S, Mungall C, Richter J: The Gene Ontology (GO) project in 2006. Nucleic Acids Res. 2006, 34: D322-D326.
Pujana MA, Han JD, Starita LM, Stevens KN, Tewari M, Ahn JS, Rennert G, Moreno V, Kirchhoff T, Gold B: Network modeling links breast cancer susceptibility and centrosome dysfunction. Nat Genet. 2007, 39 (11): 1338-1349. 10.1038/ng.2007.2.
Prieto C, Risueno A, Fontanillo C, De las Rivas J: Human gene coexpression landscape: confident network derived from tissue transcriptomic profiles. PLoS One. 2008, 3 (12): e3911-10.1371/journal.pone.0003911.
Barabasi AL, Oltvai ZN: Network biology: understanding the cell’s functional organization. Nat rev. 2004, 5 (2): 101-113. 10.1038/nrg1272.
Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297 (5586): 1551-1555. 10.1126/science.1073374.
Chen F, Zhu HH, Zhou LF, Li J, Zhao LY, Wu SS, Wang J, Liu W, Chen Z: Genes related to the very early stage of ConA-induced fulminant hepatitis: a gene-chip-based study in a mouse model. BMC Genomics. 2010, 11 (1): 240-10.1186/1471-2164-11-240.
Chen L, Xu W-w, Han M-j, Guo Y, Wang J, Cui W-q, Jia M-h, Ma Y-l, Lu L, Zhang H-T: CD3/CD28 costimulation-induced resistance to HIV-1 and cell subsets analysis. Int J Immunol. 2011, 34 (5): 297-301.
Zheng Q, Wang XJ: GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008, 36 (suppl 2): W358-W363.
Jonathan D, John P, David R: The biological context of HIV-1 host interactions reveals subtle insights into a system hijack. BMC Syst Biol. 2010, 4: 80-10.1186/1752-0509-4-80.
Jeong H, Mason SP, Barabasi AL, Oltvai ZN: Lethality and centrality in protein networks. Nature. 2001, 411 (6833): 41-42. 10.1038/35075138.
Oswald-Richter K, Grill SM, Leelawong M, Tseng M, Kalams SA, Hulgan T, Haas DW, Unutmaz D: Identification of a CCR5-expressing T cell subset that is resistant to R5-tropic HIV infection. PLoS Pathog. 2007, 3 (4): e58-10.1371/journal.ppat.0030058.
Creson JR, Lin AA, Li Q, Broad DF, Roberts MR, Anderson SJ: The mode and duration of anti-CD28 costimulation determine resistance to infection by macrophage-tropic strains of human immunodeficiency virus type 1 in vitro. J Virol. 1999, 73 (11): 9337-9347.
Riley JL, Levine BL, Craighead N, Francomano T, Kim D, Carroll RG, June CH: Native and memory CD4 T cells differ in their susceptibilities to human immunodeficiency virus type 1 infection following CD28 costimulation: Implications for transmission and pathogenesis. J Virol. 1998, 72 (10): 8273-8280.
Naghavi MH, Goff SP: Retroviral proteins that interact with the host cell cytoskeleton. Curr Opin Immunol. 2007, 19 (4): 402-407. 10.1016/j.coi.2007.07.003.
Liu Y, Belkina NV, Shaw S: HIV Infection of T Cells: Actin-in and Actin-out. Sci Signal. 2009, 2 (66): pe23-10.1126/scisignal.266pe23.
Foley JF, Yu CR, Solow R, Yacobucci M, Peden KWC, Farber JM: Roles for CXC chemokine ligands 10 and 11 in recruiting CD4+ T cells to HIV-1-infected monocyte-derived macrophages, dendritic cells, and lymph nodes. J Immunol. 2005, 174 (8): 4892.
Butler SL, Johnson EP, Bushman FD: Human immunodeficiency virus cDNA metabolism: notable stability of two-long terminal repeat circles. J Virol. 2002, 76 (8): 3739-3747. 10.1128/JVI.76.8.3739-3747.2002.
Schwartz O, Maréchal V, Friguet B, Arenzana-Seisdedos F, Heard JM: Antiviral activity of the proteasome on incoming human immunodeficiency virus type 1. J Virol. 1998, 72 (5): 3845-3850.
Baumeister W, Walz J, Proteolysis C: The proteasome: paradigm review of a self-compartmentalizing protease. Cell. 1998, 92: 367-380. 10.1016/S0092-8674(00)80929-0.
Coux O, Tanaka K, Goldberg AL: Structure and functions of the 20S and 26S proteasomes. Annu Rev Biochem. 1996, 65 (1): 801-847. 10.1146/annurev.bi.65.070196.004101.
Wei BL, Denton PW, O’Neill E, Luo T, Foster JL, Garcia JV: Inhibition of lysosome and proteasome function enhances human immunodeficiency virus type 1 infection. J Virol. 2005, 79 (9): 5705-5712. 10.1128/JVI.79.9.5705-5712.2005.
Sankaran S, Guadalupe M, Reay E, George MD, Flamm J, Prindiville T, Dandekar S: Gut mucosal T cell responses and gene expression correlate with protection against disease in long-term HIV-1-infected nonprogressors. Proc Natl Acad Sci USA. 2005, 102 (28): 9860-9865. 10.1073/pnas.0503463102.
Salgado M, Lopez-Romero P, Callejas S, Lopez M, Labarga P, Dopazo A, Soriano V, Rodes B: Characterization of host genetic expression patterns in HIV-infected individuals with divergent disease progression. Virology. 2011, 411 (1): 103-112. 10.1016/j.virol.2010.12.037.
Varbanov M, Espert L, Biard-Piechaczyk M: Mechanisms of CD4 T-cell depletion triggered by HIV-1 viral proteins. AIDS Rev. 2006, 8 (4): 221-236.
MacPherson JI, Dickerson JE, Pinney JW, Robertson DL: Patterns of HIV-1 protein interaction identify perturbed host-cellular subsystems. PLoS Comput Biol. 2010, 6 (7): e1000863-10.1371/journal.pcbi.1000863.
Krishnan V, Zeichner SL: Host cell gene expression during human immunodeficiency virus type I latency and reactivation and effects of targeting genes that are differentially expressed in viral latency. J Virol. 2004, 78 (17): 9458-9473. 10.1128/JVI.78.17.9458-9473.2004.
Chun TW, Justement JS, Lempicki RA, Yang J, Dennis G, Hallahan CW, Sanford C, Pandya P, Liu S, McLaughlin M: Gene expression and viral prodution in latently infected, resting CD4+ T cells in viremic versus aviremic HIV-infected individuals. Proc Natl Acad Sci. 2003, 100 (4): 1908-1913. 10.1073/pnas.0437640100.
Huang T, Xu Z, Chen L, Cai YD, Kong X: Computational Analysis of HIV-1 Resistance Based on Gene Expression Profiles and the Virus-Host Interaction Network. PLoS One. 2011, 6 (3): e17291-10.1371/journal.pone.0017291.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/6/15/prepub
We thank the National Engineering Center for Biochip in Shanghai for the RNA quality control and microarray hybridization, and Genminix Informatics Ltd., Co and Novel Bioinformatics Ltd., Co for their technical assistance in the bioinformatic analysis. This work was supported by the “100 Talents” Program of the Chinese Academy of Sciences (No. A0820), National Natural Science Foundation of China (NSFC, Nos. 31070817 & 30771951), and Applied Fundamental Research Program of Yunnan Province (Nos. 2008CC001 & 2010CD101).
The authors declared that they have no competing interests.
HTZ conceived and designed the study. WWX, MJH & LC performed the experiments. WWX & DC analyzed the data. YG contributed reagents/materials. WWX, DQL & HTZ interpreted the data and wrote the paper. AW revised the manuscript. All authors read and approved the final manuscript.
Wen-Wen Xu, Miao-Jun Han contributed equally to this work.
Electronic supplementary material
Additional file 2: Figure S2: Comparative Ki67 staining among “R”, “P” and “B” cells. Cells of “R” were stained with anti-Ki67 and cells of “P” and “B” were stained with anti-Ki67 on day 3 and day 6 and it also showed significant T cell proliferation during the culture period. Red line: isotype control; green line: Ki67. (TIFF 555 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.