Transcriptome analysis of primary monocytes from HIV-positive patients with differential responses to antiretroviral therapy
© Wu et al.; licensee BioMed Central Ltd. 2013
Received: 25 October 2013
Accepted: 18 December 2013
Published: 27 December 2013
Despite the significant contributions of monocytes to HIV persistence, the HIV-monocyte interaction remains elusive. For patients on antiretroviral therapy, previous studies observed a virological suppression rate of >70% and suggested complete viral suppression as the primary goal. Although some studies have reported genetic dysregulations associated with HIV disease progression, research on ex vivo-derived monocytic transcriptomes from HIV+ patients with differential responses to therapy is limited. This study investigated the monocytic transcriptome distinctions between patients with sustained virus suppression and those with virological failure during highly active antiretroviral therapy (HAART).
Genome-wide transcriptomes of primary monocytes from five HIV+ patients on HAART who sustainably controlled HIV to below detection level (BDL), five HIV+ patients on HAART who consecutively experienced viremia, and four healthy HIV sero-negative controls were analyzed using Illumina microarray. Pairwise comparisons were performed to identify differentially expressed genes followed by quantitative PCR validation. Gene set enrichment analysis was used to check the consistency of our dataset with previous studies, as well as to detect the global dysregulations of the biological pathways in monocytes between viremic patients and BDLs.
Pairwise comparisons including viremic patients versus controls, BDL versus controls, and viremic patients versus BDLs identified 473, 76, and 59 differentially expressed genes (fold change > 2 and FDR < 0.05), respectively. The reliability of our dataset was confirmed by gene set enrichment analysis showing that 6 out of 10 published gene lists were significantly enriched (FDR < 0.01) in at least one of the three pairwise comparisons. In the comparison of viremic patients versus BDLs, gene set enrichment analysis revealed that the pathways characterizing the primary functions of monocytes including antigen processing and presentation, FcγR mediated phagocytosis, and chemokine signaling were significantly up-regulated in viremic patients.
This study revealed the first transcriptome distinctions in monocytes between viremic patients and BDLs on HAART. Our results reflected the outcome balanced between the subversion of the monocyte transcriptome by HIV and the compensatory effect adapted by host cells. The up-regulation of antigen presentation pathway in viremic patients particularly highlighted the role of the interface between innate and adaptive immunity in HIV disease progression.
KeywordsHIV Monocyte Transcriptome HAART Virological failure Antigen presentation
Monocytes, a key cell type in innate and adaptive immunity, have a propensity to differentiate into macrophages or dendritic cells [1, 2]. This differentiation ability, along with the activities of antigen presentation, migration, chemotaxis, and phagocytosis , enables them to play crucial roles in HIV pathogenesis. Though less permissive to HIV infection than T cells and macrophages [4, 5], monocytes could be infected by HIV [6, 7] and infectious virus can be isolated from circulating monocytes in untreated patients and highly active antiretroviral therapy (HAART) responders [4, 8]. By harboring and trafficking HIV into various tissue compartments through differentiating into tissue macrophages or dendritic cells, monocytes serve as important viral reservoirs [9, 10]. During therapy, monocytes can maintain HIV replication throughout HAART as antiretroviral drugs may not block viral replication in monocytes as efficiently as they do in CD4+ T cells . Moreover, undifferentiated monocytic precursor cells (such as CD34+ progenitor cells) infected by HIV may pass the virus to progeny monocytes and keep on renewing the viral pool in peripheral blood monocytes [11, 12].
Despite the significant contributions of monocytes to HIV persistence, the underlying pathogenic mechanism is not fully understood. To better understand HIV pathogenesis at the genomic level, genome-wide transcriptomic studies of monocytes and monocyte-derived macrophages (MDM) have been carried out. For example, studies using monocytes/MDM infected by HIV in vitro revealed the key areas of monocyte dysfunctions related to inflammation , cytokine networks , cell cycle , cytoskeleton , and signaling pathways [17, 18]. Other studies using ex vivo-derived monocytes identified an anti-apoptosis gene signature in viremic patients , a mixed phenotype with both increased and decreased pro-inflammatory features in patients with high viral load [20, 21], and a novel candidate gene NAMPT correlating with the viral load in therapy-naïve patients . Recently, by comparing the monocyte transcriptomes from HIV+ progressors and therapy naïve non-progressors, we have shown the systematic alteration of the interrelated pathways such as Toll-like receptor (TLR) signaling and cytokine-cytokine receptor interaction in viremic patients . Although these studies have provided large datasets to facilitate our understanding, current knowledge on the dysregulations of monocytic transcriptome during HIV disease progression remains far from complete. In particular, none of the previous studies has looked into the global dysregulations of the biological pathways in monocytes from patients with sustained virus suppression versus those with virological failure during HAART, as we previously did on T cell subsets . Regarding the virological suppression rate by HAART, the study on the patients receiving HAART for 12 months in Nigeria has observed a virological suppression rate of 76.7% versus a virological failure rate of 23.4% , whereas another study on the UK cohort has reported that 73.5% of the patients initiating HAART achieved complete virological suppression within 6 months .
Clinical profiles of study patients
Viral load (copies/ml)
Total months on HAART
RNA integrity number
Cluster analysis and identification of differentially expressed genes
Differentially expressed genes in the comparison of VIR versus BDL
Entrez gene name
Defensin, alpha 1
Defensin, alpha 3
Host cell factor C1 regulator 1 (XPO1 dependent)
Chromosome 20 open reading frame 27
E4F transcription factor 1
Transmembrane and ubiquitin-like domain containing 2
Zinc finger protein 467
Solute carrier family 43, member 2
Chromosome 16 open reading frame 58
Differentially expressed in FDCP 6 homolog (mouse)
Tyrosylprotein sulfotransferase 2
Discs, large (Drosophila) homolog-associated protein 4
Lymphocyte-specific protein 1
Mitogen-activated protein kinase kinase 3
OAF homolog (Drosophila)
Mahogunin ring finger 1, E3 ubiquitin protein ligase
PML-RARA regulated adaptor molecule 1
Zinc finger protein, FOG family member 1
Complement component 5a receptor 1
E1A binding protein p300
Cleft lip and palate associated transmembrane protein 1
GEM interacting protein
Solute carrier family 22, member 18
SH3KBP1 binding protein 1
Neutrophil cytosolic factor 1B pseudogene
Lymphotoxin beta receptor (TNFR superfamily, member 3)
Family with sequence similarity 89, member B
Glycosylphosphatidylinositol anchor attachment 1
Aldehyde dehydrogenase 3 family, member B1
Protein phosphatase 2, regulatory subunit A, alpha
Trafficking protein particle complex 9
Syntaxin binding protein 2
Complement component 1, q subcomponent, B chain
Actinin, alpha 4
Josephin domain containing 2
Rho GDP dissociation inhibitor (GDI) alpha
T-cell, immune regulator 1, ATPase, H + transporting, lysosomal V0 subunit A3
Integrin, alpha L (antigen CD11A (p180), lymphocyte function-associated antigen 1; alpha polypeptide)
Ankyrin repeat and BTB (POZ) domain containing 1
Family with sequence similarity 108, member A3, pseudogene
Paroxysmal nonkinesigenic dyskinesia
Zinc finger protein 580
Mitogen-activated protein kinase kinase kinase 11
Jun B proto-oncogene
Spleen focus forming virus (SFFV) proviral integration oncogene spi1
Lectin, galactoside-binding, soluble, 9
Basigin (Ok blood group)
Tweety homolog 3 (Drosophila)
Uncoupling protein 2 (mitochondrial, proton carrier)
Ring finger protein 19B
Major histocompatibility complex, class II, DR beta 5
qPCR validation and the comparison to previous microarray studies
qPCR confirmation of differentially expressed genes
The GSEA of our dataset compared with the published gene lists derived from monocyte/MDM transcriptomes modulated by HIV
Gene set name
Pathways/biological functions DEGs involved in HIV infection/disease progression
Gene set size
VIRvs BDL FDR
VIRvs CTR FDR
BDLvs CTR FDR
Chemokine and cytokine
in vitro gp120-treated vs mock-treated MDM
Cell cycle regulators, translation, cell signaling, TNF, MAPK
in vitro HIV-infected vs mock-treated MDM
Interferon stimulated genes, host defense genes
in vitro HIV-infected vs mock-treated MDM
Signal transduction, transcription, cell cycle and apoptosis, adhesion molecules and receptors, chemokines and cytokines, proteases and protease inhibitors, metabolism
in vitro HIV-infected vs mock-treated MDM
Signaling components, transcription factors, cytokines, apoptotic and anti-apoptotic factors, growth factors, anti-HIV infection genes
in vitro HIV-infected vs mock-infected U937 human promonocytes
Type I interferon responses, NF-κB, mitogen-activated protein kinase, Jun signaling pathways, general immune activation, immune down-regulation, protein degradation, protein secretion, and apoptosis
ex vivo correlations between changes in gene expression values and changes in monocyte cytokine levels in HIV+ patients on and off therapy
Apoptosis-related gene signatures, TNF-α signaling, CD40L/CD40 signaling, MAPK signaling, p53 modulation
ex vivo HIV+ patients vs healthy controls
Van_HIV_Serostatus _associated _validated
Apoptosis, cell cycle, transcriptional regulation, immune response, protein trafficking, lipid metabolism
ex vivo HIV+ patients vs healthy controls
TLR2-agonist stimulated gene signature TNF (NFκB), p53 and MAPK networks
ex vivo and in vitro HIV+ patients vs healthy controls stimulated vs non-stimulated controls
ERK/MAPK, TNF/IL6 (NFκB) and p53 gene networks, apoptosis-related gene signatures
ex vivo HIV+ patients vs healthy controls
Pathway analysis by GSEA
As we aimed to identify monocyte dysfunction in relation to HIV disease progression, our subsequent pathway analysis focused on the comparison between VIR and BDL groups. GSEA was attempted using the KEGG pathways including 186 gene sets. Unlike DEG approach, GSEA uses the whole gene expression dataset to identify enriched pathways and is better able to detect small coordinated changes in gene expression in the context of gene set [31, 32].
The ten immune-related pathways significantly up-regulated in viremic patients versus BDLs
Gene set name
Core enrichment genes contributing to the pathway enrichment
HLA-DRB5, HLA-DRB1, LTBR, ITGB7, CD40, TNFSF13, HLA-DOA, CCR9
HLA-DRB5, HLA-DRB1, CD4, HLA-C, IFI30, NFYC, TAPBP, HSPA1B, CD74, HLA-DOA, TAP2, HLA-E, PDIA3, PSME1, HSP90AB1, RFX5, CALR, KIR2DS5, HLA-F
ACTN4, EP300, TCF7L2, PTPN6, IQGAP1, CSNK2B, SRC, FYN, RAC1, CSNK2A2, WASF2, PVRL4, PTPN1, BAIAP2, CTNNB1, PVRL2, ACTB
HLA-DRB5, HLA-DRB1, FLT3LG, CD4, TNF, CD37, IL1B, CD9, CD36, EPOR, ANPEP, CD1A, CD19, ITGA2B, FCER2, CD1D, FCGR1A
HLA-DRB5, HLA-DRB1, ITGAL, ITGB7, CD40, CD4, HLA-C, SELL, PECAM1, SPN, HLA-DOA, ICAM2, HLA-E, ICAM3, CLDN3, PVRL2, JAM3, HLA-F, CD28, CLDN18, PVR, CD99, ICOS, PDCD1LG2, CD80, CLDN5, NEO1, PDCD1, ITGAM, PTPRF, HLA-DRA, CNTN2, HLA-B, CD58, NLGN3, CTLA4, ITGA4, SIGLEC1
ARSA, ACP2, SGSH, CTSZ, MCOLN1, CTSA, CTSD, SORT1, FUCA1, MAN2B1, GBA, TCIRG1, NAPSA, IDUA, NAGPA, M6PR, ATP6V0D1, NEU1, ACP5, GGA1, LAPTM5, CD63, AP1S1, AP3B1, AP3D1, NAGA, CLN3, IDS, TPP1
YWHAE, ARHGDIA, RPS6KA4, SH2B2, MAP2K2, MAPKAPK2, MAPK7, SORT1, CSK, RELA, SHC1, PIK3CG, KIDINS220, CRK, IRAK1, BAX, SH2B3, RAC1, MAP2K5, AKT1, CALM1, MAP3K3, MAPK9, ABL1, RPS6KA3, SH2B1, IRS1, RAF1, HRAS, MAPK14
MARCKSL1, NCF1, PIP5K1C, VASP, HCK, PIK3CG, DNM2, CRK, SYK, ARPC4, ARPC1A, SPHK1, RAC1, ARPC5L, LYN, AKT1, WASF2, VAV1, ARPC1B, LOC653888, FCGR1A, RAF1, PRKCA, GSN, DNM1, PAK1, VAV2, CFL1, AKT2
PLCB2, NCF1, GNB2, ELMO1, CX3CR1, GRK6, PPBP, GNAI2, CSK, GNB4, RELA, SHC1, HCK, ARRB1, PIK3CG, CRK, RASGRP2, CCR1, GNGT2, GNG5, GNG3, STAT5B, RAC1, LYN, AKT1, CCL2, JAK2, VAV1, GNG8, CCR2, CCR9, PRKACA, CCR7, XCR1, TIAM2, GNG12, RAF1, ARRB2, HRAS, CXCR5, PTK2B, PAK1, CCR3, VAV2
C1QB, C1QC, C5AR1, SERPING1, C1QA, C2, KLKB1, C3AR1, THBD, SERPINA1, C4A, SERPINA5
Out of the ten immune-related pathways (Table 5), two were involved in cell differentiation and development (hematopoietic cell lineage and neurotrophin signaling), two were associated with transendothelial migration (cell adhesion molecules and adherens junctions), and the remaining six covered all the major aspects of innate immunity, including chemokine signaling, IgA production, complement cascade, phagocytosis, lysosome, and antigen presentation. All the core enrichment genes contributing to the up-regulation of the immune-related pathways are listed in Table 5. Since the pathways of antigen presentation, phagocytosis, and chemokine signaling characterize the primary functions of monocytes, we further inspected their core enrichment genes.
Up-regulation of antigen presentation pathway in VIR versus BDL
Up-regulation of phagocytosis pathway in VIR versus BDL
Up-regulation of chemokine signaling pathway in VIR versus BDL
Promoter motif analysis by GSEA
Promoter motif analysis by GSEA used gene sets that contained gene members sharing the same transcription factor binding site defined in the TRANSFAC database . This analysis can identify the coordinated changes of the genes under the control of a certain transcriptional regulator. The GSEA revealed that 37 and 42 gene sets were significantly up- and down- regulated in the VIR group (FDR < 0.05) compared to the BDL group (Additional file 4). While the majority of the significantly altered gene sets contained regulatory motifs matching for the annotated transcription factors, 12 gene sets contained motifs which did not match any known transcription factor. The most significantly up- and down- regulated gene sets in the VIR group (both FDR = 0.03) had gene members sharing the regulatory motif for transcription factors NFκB and MYC, respectively.
Our study has provided the first snapshot into the transcriptome distinctions of the primary monocytes between HIV+ patients on HAART who consecutively experienced viremia and HIV+ patients on HAART who sustainably controlled HIV to below detection level (Table 1). The main objective of the study was to identify genomic signatures associated with HIV disease progression. The relevance of the identified DEGs was initially confirmed by querying HIV interaction database using DAVID (Table 2; Additional file 2). The comparison between our dataset and previous microarray studies on monocyte/MDM transcriptomes further confirmed the reliability of our data (Table 4). The previously identified gene sets, which were highly significant (FDR < 0.01) in our dataset, reflected the altered biological functions including cytokine networks, cell cycle, signaling pathways, metabolism, immune responses, and transcriptional regulation. These biological themes were also represented by the gene sets which only reached the relaxed significance threshold (FDR < 0.25). The failure of these gene sets to achieve higher statistical cut-offs could possibly be attributed to certain biological and/or technical variations across multiple microarray studies. For example, differences in the cohort such as age, gender, and ethnicity could contribute to the biological variations, whereas microarray inter-platform differences could contribute to the technical variations. Compared to DEG approach which focused on discrete genes, GSEA enabled a more comprehensive detection of genes contributing to the enrichment of the pathways correlated with disease progression. Therefore, the subsequent discussion focused on the pathways significantly associated with progressive phase of HIV disease revealed by GSEA. Since the immune-related pathways have the most direct relevance to HIV disease and this aspect has been under constant investigations, our subsequent discussions shall center on the immune-related pathways including antigen presentation, phagocytosis, and chemokine signaling.
At the transcriptome level, we observed an overall up-regulation of antigen processing and presentation pathway in the VIR group compared to the BDL group (FDR = 0.016; Figure 3; Table 5), which was manifested by the up-regulation of both endogenous (MHC-I) and exogenous (MHC-II) signaling branches. Five molecules from MHC-II pathway were detected as core enrichment genes including IFI30, CD74, HLA-DOA, HLA-DRB1, and HLA-DRB5 (Table 5). The coordinated up-regulation of HLA-DR was consistent with the previous studies reporting higher expression levels of surface HLA-DR on monocytes from HIV+ individuals compared to sero-negative controls [34–36]. As for the MHC-I pathway, the majority of the core enrichment genes were reported to have interactions with HIV, such as HSP90AB1  and PA28 (PSME1) . Moreover, 3 MHC-I molecules including HLA-C, HLA-E, and HLA-F were also detected as core enrichment genes with the coordinated up-regulation in the viremic patients. The importance of HLA-C was highlighted in a very recent study showing that increased HLA-C expression was correlated with increased likelihood of cytotoxic T lymphocyte (CTL) responses and frequency of viral escape mutation . In addition, the substantially increased expression of MHC-I including HLA-C has been reported in monocytes in HIV progressors compared to HAART-suppressed patients, which was suggested to contribute to the generation of the dysfunctional naive CD8-low T cells that emerged during disease progression . Taken together, our observation that the up-regulation of HLA-C expression on monocytes was associated with plasma viremia could possibly reflect a compensatory effect imposed by impaired CTL responses. However, this compensatory effect may be subverted by HIV, resulting in further CD8+ T cell dysfunction as proposed by Favre et al. .
As antigen presentation pathway is right at the interface between innate and adaptive immunity, its significant up-regulation in the viremic patients provides direct evidence that adaptive immune response is perturbed through the components of innate immunity during HIV disease progression. This interference is systematic throughout the process of antigen presentation as demonstrated by the core enrichment genes covering not only the aforementioned MHC molecules, but also the genes associated with antigen digestion, loading and transportation (Table 5). In fact, the dysregulation of antigen presentation pathway is not the only pathway altered at the interface. Other innate immunity pathways crucial for the development of adaptive immunity have also been associated with HIV disease progression by this and previous studies [23, 24, 41]. For example, this study detected that complement cascade regulating both B and T cell responses  was significantly up-regulated in the viremic patients (Table 5). Consistently, complement pathway was also found to be up-regulated in the viremic patients versus BDLs by our previous study on primary CD4+ T cells . In addition, a recent study has reported that monocytes and complement system contributed to the tuberculosis-associated immune reconstitution inflammatory syndrome in HIV-TB co-infected patients . Altogether, these studies have demonstrated the importance of complement components in HIV disease progression. Recently, another key component of innate immunity, Toll-like receptor signaling (TLR) pathway, has been found to be significantly down-regulated in monocytes from viremic patients versus long-term non-progressors . The different direction of the changes may reflect the various aspects of HIV-host interactions that contribute to disease progression, such as HIV persistence and impairment of T cell functions. Despite this difference, all of the aforementioned studies point towards the adaptive immunity being perturbed at the interface where innate and adaptive immunity interact during HIV disease progression.
Consistent with the previous reports on phagocytosis dysfunction in monocytes upon HIV infection [43–45], we observed the significant up-regulation of FcγR-mediated phagocytosis pathway (HSA04666) in the VIR group versus the BDL group (FDR = 0.029; Figure 4). The expression of FCGR1A (FcγRI/CD64), the sole high affinity receptor for monomeric IgG, was coordinately increased along with other core enrichment genes in the viremic patients. This observation was consistent with the previous study reporting the increased expression of FcγRI on monocytes from acute HIV infection patients compared to chronically-infected individuals and healthy controls . In addition, the previous study suggested that monocytes might expand and/or up-regulate the expression of this high affinity FcγR in response to the burst of viral replication. This host response may also account for the increased FCGR1A expression in the viremic patients on HAART, as detected in this study.
The perturbation of phagocytosis pathway was further manifested by the up-regulation of the genes encoding for kinases in early signaling events such as SRC kinases (HCK and LYN) and SYK, as well as the genes encoding for proteins involved in oxidative burst (SPHK1/NADPH oxidase) and cytoskeleton rearrangement (Table 5). The majority of the core enrichment genes, which spread along the signaling branches of PI3K – AKT, PIP5K – VASP/WASP/WAVE – ARP2/3, and VAV – RAC – PAK1 – CFL1, have also been implicated previously in HIV pathogenesis. A recent study using gene knockdown and pharmacological inhibition of SRC and AKT identified them as mediators of HIV-induced inhibition of autophagocytosis in bystander macrophages/monocytic cells . The elevated AKT activity was also detected in HIV-infected macrophages, and PI3K/AKT inhibitors were suggested as a novel therapy for interfering with the establishment of HIV reservoirs . In the signaling branch pointing towards ARP2/3, WAVE2 was involved in the activation of the actin polymerization nucleator ARP2/3 complex, which played a role in the migration of the viral core components toward the cell nucleus and the efficient infection of HIV in cell lines [49, 50]. In addition, HIV Nef was found to activate the signaling branch of VAV – RAC – PAK (p21), which markedly enhanced the NADPH oxidase response . Nef also mediated CFL1 phosphorylation, which was crucial for maintaining actin homeostasis . Taken together, the up-regulated expression of the genes involved in various kinase cascades and the gene encoding for FcγRI may reflect the HIV interference with the host genetic machinery on the one hand, and the compensatory processes adapted by the host on the other hand during virological failure.
Chemokine signaling pathway (HSA04062) was also significantly up-regulated in the VIR group versus the BDL group (FDR = 0.038), which was manifested by the systematic up-regulation of the genes encoding for chemokine ligand 2 (CCL2/MCP-1) and chemokine receptors including CCR1, CCR2, CCR3, CCR7, CCR9, CX3CR1, CXCR5, and XCR1 (Table 5; Figure 5). The up-regulation of MCP-1 in viremic patients was also reported by Pulliam et al. , which was suggested to be involved in enhancing HIV production and spread . In addition, it was shown that the sequential activation of SRC, MAPKs, and PI3K – AKT – NFκB pathways resulted in the increased expression of MCP-1 . Given the fact that the majority of these genes overlapped with the core enrichment genes we detected, it is thus plausible to hypothesize that the up-regulation of MAPK and NFκB braches signaling through Gα subunit may be related to the coordinated up-regulation of MCP-1, which could contribute to HIV spread. This potential link could be one of the mechanisms underlying the association of the up-regulation of chemokine signaling pathway with virological failure during HAART.
In the comparison of VIR versus BDL, the promoter motif analysis yielded a list of significantly up- and down- regulated gene sets with members of each gene set containing the same binding site for a certain transcription factor (Additional file 4). The transcription factors implicated by these significant gene sets included some of the well-known regulators, such as SP (SP1 and SP3), NFκB (RELA binding motif), AP1, AP2, and CREB1, which have been demonstrated to play crucial roles in HIV transcription regulation . Interestingly, the most significantly up-regulated gene set in the VIR group (FDR = 0.03) contained the binding motif for RELA, which was also identified as the core enrichment gene in chemokine signaling pathway in the VIR group. This overlap not only demonstrated the consistency between pathway and promoter motif analysis, but also implicated the complexity of the transcriptional regulation underlying HIV-monocyte interaction. In addition to the aforementioned transcription factors, some of the significant gene sets contained the binding motifs without matching any known transcription factor, which may suggest potential novel regulators that warrant future investigations.
A few limitations of this study should be noted. First, this study used a cross-sectional design, which could not provide dynamic findings from a longitudinal perspective. Secondly, this study used a relatively small sample size appropriate for the pilot investigation and future studies using larger sample size are thus warranted to further confirm the results. Finally, although the GSEA identified a panel of significantly altered gene sets with high relevance to HIV disease progression during therapy, the predictive nature of GSEA as the common limitation of statistical tools should be noted and interpretations should be made with caution. To overcome this limitation, future biological experiments should be conducted to directly confirm and further explore these findings.
This study has revealed the first transcriptome distinctions in monocytes between HIV+ patients on HAART who consecutively experienced plasma viremia and HIV+ patients on HAART who sustainably controlled plasma viremia to below detection level. Compared to the BDL group, the pathways characterizing the primary functions of monocytes including antigen processing and presentation, FcγR mediated phagocytosis, and chemokine signaling, were significantly up-regulated in the VIR group. Our results reflected the outcome balanced between the subversion of monocyte transcriptome by HIV and the compensatory effects adapted by host cells. Furthermore, the altered pathways of antigen presentation and complement cascade highlighted that HIV manipulated adaptive immune response via innate immunity components at the interface of innate-adaptive immunity during disease progression. These data offered new comparative insights into the perturbed genetic networks of ex vivo-derived monocytes subverted by HIV during disease progression on therapy. In-depth functional studies on the regulation of these pathways and the corresponding core enrichment genes along with proteomic analysis may further confirm our findings and provide detailed molecular mechanisms underlying HIV-monocyte interaction.
Patient profiles and collection protocol
Five HIV+ patients on HAART who sustainably controlled HIV to below detection level (viral load < 40 copies of HIV RNA/ml), five HIV+ patients on HAART who consecutively experienced viremia, and four healthy HIV sero-negative controls were studied (Table 1). Patients with viral load < 1,000 were chosen to represent patients on HAART with virological failure because the previous study has demonstrated that the transcriptome profiling of these patients was more homogenous . It has also been shown that monocyte transcriptomes from patients with viral load < 1,000 and those with higher viral load exhibited similar dysregulated pathways during disease progression. Patients in the VIR and BDL groups were on HAART for > 52 weeks. The BDL and VIR groups had a broad range of CD4+ T cell counts (100–796 cells/μl) as our objective was to find unbiased differences of gene expression profiling of monocytes between patients on HAART with sustained virus suppression and virological failure regardless of the degree of T cell decline. This grouping criterion by viral load has been successfully used by the previous studies [20, 23, 24]. These patients received two NRTIs (zidovudine, lamivudine, stavudine, emtricitabine, tenofovir) in association with one or two protease inhibitors (darunavir, ritonavir, indinavir, saquinavir, atazanavir). Ten patients were from the HIV clinic at Westmead Hospital and the four healthy controls were from the Australian Red Cross Blood Service in Sydney. This study was approved by the Sydney West Area Health Services Research Ethics Committee, and all blood samples were collected after individual informed written consent.
Purification of CD14+ monocytes and RNA isolation
A single blood sample (10–20 ml in EDTA) was obtained from each patient. After separation of plasma, primary PBMCs were isolated immediately after obtaining blood samples by Ficoll-gradient centrifugation and purified. This aspect was strictly followed in our experiments because of previously described lower RNA yields and possible changes in gene expression profiles upon storage of blood . CD14+ monocytes were then obtained by positive isolation with antibody-conjugated magnetic beads according to the manufacturer’s instructions (Miltenyi Biotech, Germany) with a purity > 98.6% as verified by flow cytometry. Binding of antibody to CD14 does not trigger signal transduction and a previous study has clearly demonstrated that CD14 positive selection does not alter cellular transcriptome by comparing gene expression profiles in parallel using either positive or negative selection . Total RNA was isolated from purified cells using RNeasy Mini kit (Qiagen Pty Ltd., Clifton Hill, Victoria, Australia) with an integrated step of on-column DNase treatment.
cRNA preparation, microarray hybridization and scanning
RNA quality was checked by Agilent Bioanalyzer and RNA Integrity Scores were higher than 9 for all the samples (Table 1). cRNA amplification and labeling with biotin were performed using Illumina TotalPrep RNA amplification kit (Ambion, Inc., Austin, USA) with 250 ng total RNA as input material. cRNA yields were quantified with Agilent Bioanalyzer and 750 ng cRNAs were hybridized to Illumina HumanHT-12 v3 Expression BeadChips (Illumina, Inc., San Diego, USA). Each chip contains 12 arrays and each array contains >48,000 gene transcripts, of which, 46,000 are derived from human genes in the National Center for Biotechnology Information (NCBI) Reference Sequence (RefSeq) and UniGene databases. All reagents and equipment used for hybridization were purchased from Illumina, Inc. According to the manufacturer’s protocol, cRNAs were hybridized to arrays for 16 hours at 58°C before being washed and stained with streptavidin-Cy3. Then the beadchips were centrifuged to dry and scanned on the Illumina BeadArray Reader confocal scanner. To minimize the batch effect, the microarray chips were all processed at the single site using the same platform with the identical setting of the parameters by the same experimenter. The microarray dataset has been submitted to GEO (Accession Number GSE52900).
Analysis of differentially expressed genes
The quality of the entire data set was assessed by box plot and density plot of bead intensities, density plot of coefficient of variance, pairwise MAplot, pairwise plot with microarray correlation, cluster dendrogram, and non-metric multidimensional scaling using R/Bioconductor and the lumi package . Based on the quality assessment, all 14 samples were deemed suitable for further analysis. Data normalization was performed using log2 transform and a robust spline normalization (RSN) implemented in the lumi package for R/Bioconductor [56, 57]. Cluster analysis of gene expression profiling was carried out using dist and hclust functions from R stats package. Euclidean distance and complete linkage were used for distance metric and linkage criterion, respectively. To reduce false positives of differentially expressed genes, genes below detectable limit (based on a detection p value cut-off 0.01) were removed from the dataset. A linear model fit in conjunction with an empirical Bayes statistics was used to identify candidate DEGs . P values were corrected for multiple testing using FDR adjustment implemented in lumi package. Pairwise comparisons for the 3 groups were carried out and candidate DEGs with fold change >2 and FDR <0.05 were identified for each of the comparisons. The list of the DEGs derived from the comparison of VIR versus BDL was uploaded to DAVID for the detection of the DEGs showing overlap with the genes in HIV interaction database at NCBI [59, 60].
Gene set enrichment analysis
GSEA was used for the comparison of our dataset with the published DEG lists from the previous studies (in vivo and ex vivo since 2002) [15, 16, 19, 21, 22, 27–30], the investigation of global dysregulations of the biological pathways, and the promoter motif analysis. For the comparison with previous studies, 10 DEG lists were used from the studies on monocyte/MDM transcriptomes modulated by HIV (Table 4). For the pathway investigation, the gene sets were from MsigDB , catalog C2 functional sets, subcatalog KEGG pathways, which included 186 gene sets from pathway databases. For the promoter motif analysis, C3 motif gene sets which contained gene members sharing the same transcription factor binding site were used. This collection included 615 gene sets and each of them was annotated by a TRANSFAC record [33, 61].
Instead of focusing on discrete DEGs, we analyzed the entire transcriptome data with GSEA to identify genes coordinately regulated in predefined gene sets from various biological pathways . For each group comparison, GSEA was performed using the normalized data of entire 48,803 transcripts (GSEA version 2.07, Broad Institute http://www.broad.mit.edu/gsea). First, a ranked list was obtained by ranking all genes according to the correlation between their expression and the group distinction using the metric signal to noise ratio. Then the association between a given gene set and the group was measured by the non-parametric running sum statistic termed the enrichment score (ES), which was calculated by walking down the ranked list (increasing ES when encountering a gene in the given gene set and decreasing ES when encountering a gene not in the gene set). To estimate the statistical significance of the ES, a nominal p value was calculated by permuting the genes 1,000 times. To adjust for multiple hypothesis testing, the maximum ES was normalized to account for the gene set size (NES) and the false discovery rate (FDR) corresponding to each NES was also calculated. Along with the pathway enrichment results, the details report for each significant pathway was simultaneously generated. This report listed the details of each gene member in columns, one of which indicated whether this gene member was “core enrichment gene” or not. The core enrichment genes account for the enrichment signal of the pathway and the inspection of them can reveal a biologically important subset within the pathway .
Real-time quantitative PCR
Ten genes and 14 pairs of group comparison were selected for validation based on the coverage of different levels and directions of fold change, different ranges of FDR values, and/or biological significance. Purified total cellular RNA was used for reverse transcription with oligo d (T) and Superscript III followed by RNase H treatment (Invitrogen Life Technologies). The cDNA was then subjected to qPCR in a 96-well format in triplicate reactions with defined primers and SYBR Green (Invitrogen Life Technologies). The qPCR reactions were carried out for the extended cohort consisting of 10 viremic patients, 10 BDLs, and 9 healthy controls (including all the original samples used in the microarray) using Mx3005P™ QPCR System (Stratagene). The mean expressions of the tested genes in each group were obtained and the housekeeping gene GAPDH was used as an internal control and the normalizer for all data. The fold change was calculated by the relative quantitation method 2-(ddCt). Primer sequences for each transcript are available from the authors upon request.
This work was funded by the grant from the ACH2. JQW received an NHMRC postdoctoral training fellowship (GNT1016870). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- Randolph GJ, Inaba K, Robbiani DF, Steinman RM, Muller WA: Differentiation of phagocytic monocytes into lymph node dendritic cells in vivo. Immunity 1999, 11: 753-761. 10.1016/S1074-7613(00)80149-1PubMedView ArticleGoogle Scholar
- Geissmann F: The origin of dendritic cells. Nat Immunol 2007, 8: 558-560. 10.1038/ni0607-558PubMedView ArticleGoogle Scholar
- Sasse TR, Wu JQ, Li Z, Saksena NK: Monocytes and their role in human immunodeficiency virus pathogenesis. Am J Infect Dis 2012, 8: 92-105. 10.3844/ajidsp.2012.92.105View ArticleGoogle Scholar
- McElrath MJ, Pruett JE, Cohn ZA: Mononuclear phagocytes of blood and bone marrow: comparative roles as viral reservoirs in human immunodeficiency virus type 1 infections. Proc Natl Acad Sci USA 1989, 86: 675-679. 10.1073/pnas.86.2.675PubMedPubMed CentralView ArticleGoogle Scholar
- Lewin SR, Kirihara J, Sonza S, Irving L, Mills J, Crowe SM: HIV-1 DNA and mRNA concentrations are similar in peripheral blood monocytes and alveolar macrophages in HIV-1-infected individuals. AIDS 1998, 12: 719-727. 10.1097/00002030-199807000-00008PubMedView ArticleGoogle Scholar
- Gendelman HE, Orenstein JM, Martin MA, Ferrua C, Mitra R, Phipps T, Wahl LA, Lane HC, Fauci AS, Burke DS, et al.: Efficient isolation and propagation of human immunodeficiency virus on recombinant colony-stimulating factor 1-treated monocytes. J Exp Med 1988, 167: 1428-1441. 10.1084/jem.167.4.1428PubMedView ArticleGoogle Scholar
- Collman R, Hassan NF, Walker R, Godfrey B, Cutilli J, Hastings JC, Friedman H, Douglas SD, Nathanson N: Infection of monocyte-derived macrophages with human immunodeficiency virus type 1 (HIV-1). Monocyte-tropic and lymphocyte-tropic strains of HIV-1 show distinctive patterns of replication in a panel of cell types. J Exp Med 1989, 170: 1149-1163. 10.1084/jem.170.4.1149PubMedView ArticleGoogle Scholar
- Crowe SM, Sonza S: HIV-1 can be recovered from a variety of cells including peripheral blood monocytes of patients receiving highly active antiretroviral therapy: a further obstacle to eradication. J Leukoc Biol 2000, 68: 345-350.PubMedGoogle Scholar
- Fulcher JA, Hwangbo Y, Zioni R, Nickle D, Lin X, Heath L, Mullins JI, Corey L, Zhu T: Compartmentalization of human immunodeficiency virus type 1 between blood monocytes and CD4+ T cells during infection. J Virol 2004, 78: 7883-7893. 10.1128/JVI.78.15.7883-7893.2004PubMedPubMed CentralView ArticleGoogle Scholar
- Zhu T: HIV-1 in peripheral blood monocytes: an underrated viral source. J Antimicrob Chemother 2002, 50: 309-311. 10.1093/jac/dkf143PubMedView ArticleGoogle Scholar
- Alexaki A, Wigdahl B: HIV-1 infection of bone marrow hematopoietic progenitor cells and their role in trafficking and viral dissemination. PLoS Pathog 2008, 4: e1000215. 10.1371/journal.ppat.1000215PubMedPubMed CentralView ArticleGoogle Scholar
- Alexaki A, Liu Y, Wigdahl B: Cellular reservoirs of HIV-1 and their role in viral persistence. Curr HIV Res 2008, 6: 388-400. 10.2174/157016208785861195PubMedPubMed CentralView ArticleGoogle Scholar
- 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: 1031-1043. 10.1189/jlb.0306157PubMedView ArticleGoogle Scholar
- Izmailova E, Bertley FM, Huang Q, Makori N, Miller CJ, Young RA, Aldovini A: HIV-1 Tat reprograms immature dendritic cells to express chemoattractants for activated T cells and macrophages. Nat Med 2003, 9: 191-197. 10.1038/nm822PubMedView ArticleGoogle Scholar
- Coberley CR, Kohler JJ, Brown JN, Oshier JT, Baker HV, Popp MP, Sleasman JW, Goodenow MM: Impact on genetic networks in human macrophages by a CCR5 strain of human immunodeficiency virus type 1. J Virol 2004, 78: 11477-11486. 10.1128/JVI.78.21.11477-11486.2004PubMedPubMed CentralView ArticleGoogle Scholar
- Cicala C, Arthos J, Selig SM, Dennis G Jr, Hosack DA, Van Ryk D, Spangler ML, Steenbeke TD, Khazanie P, Gupta N, et al.: HIV envelope induces a cascade of cell signals in non-proliferating target cells that favor virus replication. Proc Natl Acad Sci USA 2002, 99: 9380-9385. 10.1073/pnas.142287999PubMedPubMed CentralView ArticleGoogle Scholar
- Kilareski EM, Shah S, Nonnemacher MR, Wigdahl B: Regulation of HIV-1 transcription in cells of the monocyte-macrophage lineage. Retrovirology 2009, 6: 118. 10.1186/1742-4690-6-118PubMedPubMed CentralView ArticleGoogle Scholar
- Herbein G, Gras G, Khan KA, Abbas W: Macrophage signaling in HIV-1 infection. Retrovirology 2010, 7: 34. 10.1186/1742-4690-7-34PubMedPubMed CentralView ArticleGoogle Scholar
- Giri MS, Nebozyhn M, Raymond A, Gekonge B, Hancock A, Creer S, Nicols C, Yousef M, Foulkes AS, Mounzer K, et al.: Circulating monocytes in HIV-1-infected viremic subjects exhibit an antiapoptosis gene signature and virus- and host-mediated apoptosis resistance. J Immunol 2009, 182: 4459-4470. 10.4049/jimmunol.0801450PubMedPubMed CentralView ArticleGoogle Scholar
- Pulliam L, Sun B, Rempel H: Invasive chronic inflammatory monocyte phenotype in subjects with high HIV-1 viral load. J Neuroimmunol 2004, 157: 93-98. 10.1016/j.jneuroim.2004.08.039PubMedView ArticleGoogle Scholar
- Tilton JC, Johnson AJ, Luskin MR, Manion MM, Yang J, Adelsberger JW, Lempicki RA, Hallahan CW, McLaughlin M, Mican JM, et al.: Diminished production of monocyte proinflammatory cytokines during human immunodeficiency virus viremia is mediated by type I interferons. J Virol 2006, 80: 11486-11497. 10.1128/JVI.00324-06PubMedPubMed CentralView ArticleGoogle Scholar
- Van den Bergh R, Florence E, Vlieghe E, Boonefaes T, Grooten J, Houthuys E, Tran HT, Gali Y, De Baetselier P, Vanham G, Raes G: Transcriptome analysis of monocyte-HIV interactions. Retrovirology 2010, 7: 53. 10.1186/1742-4690-7-53PubMedPubMed CentralView ArticleGoogle Scholar
- Wu JQ, Sasse TR, Wolkenstein G, Conceicao V, Saksena MM, Soedjono M, Perera SS, Wang B, Dwyer DE, Saksena NK: Transcriptome analysis of primary monocytes shows global down-regulation of genetic networks in HIV viremic patients versus long-term non-progressors. Virology 2013, 435: 308-319. 10.1016/j.virol.2012.10.026PubMedView ArticleGoogle Scholar
- Wu JQ, Dwyer DE, Dyer WB, Yang YH, Wang B, Saksena NK: Genome-wide analysis of primary CD4+ and CD8+ T cell transcriptomes shows evidence for a network of enriched pathways associated with HIV disease. Retrovirology 2011, 8: 18. 10.1186/1742-4690-8-18PubMedPubMed CentralView ArticleGoogle Scholar
- Anude CJ, Eze E, Onyegbutulem HC, Charurat M, Etiebet MA, Ajayi S, Dakum P, Akinwande O, Beyrer C, Abimiku A, Blattner W: Immuno-virologic outcomes and immuno-virologic discordance among adults alive and on anti-retroviral therapy at 12 months in Nigeria. BMC Infect Dis 2013, 13: 113. 10.1186/1471-2334-13-113PubMedPubMed CentralView ArticleGoogle Scholar
- Street E, Curtis H, Sabin CA, Monteiro EF, Johnson MA: British HIV Association (BHIVA) national cohort outcomes audit of patients commencing antiretrovirals from naive. HIV Med 2009, 10: 337-342. 10.1111/j.1468-1293.2009.00692.xPubMedView ArticleGoogle Scholar
- Gekonge B, Giri MS, Kossenkov AV, Nebozyhn M, Yousef M, Mounzer K, Showe L, Montaner LJ: Constitutive gene expression in monocytes from chronic HIV-1 infection overlaps with acute Toll-like receptor induced monocyte activation profiles. PLoS ONE 2012, 7: e41153. 10.1371/journal.pone.0041153PubMedPubMed CentralView ArticleGoogle Scholar
- Vazquez N, Greenwell-Wild T, Marinos NJ, Swaim WD, Nares S, Ott DE, Schubert U, Henklein P, Orenstein JM, Sporn MB, Wahl SM: Human immunodeficiency virus type 1-induced macrophage gene expression includes the p21 gene, a target for viral regulation. J Virol 2005, 79: 4479-4491. 10.1128/JVI.79.7.4479-4491.2005PubMedPubMed CentralView ArticleGoogle Scholar
- Wen W, Chen S, Cao Y, Zhu Y, Yamamoto Y: HIV-1 infection initiates changes in the expression of a wide array of genes in U937 promonocytes and HUT78 T cells. Virus Res 2005, 113: 26-35. 10.1016/j.virusres.2005.04.002PubMedView ArticleGoogle Scholar
- Woelk CH, Ottones F, Plotkin CR, Du P, Royer CD, Rought SE, Lozach J, Sasik R, Kornbluth RS, Richman DD, Corbeil J: Interferon gene expression following HIV type 1 infection of monocyte-derived macrophages. AIDS Res Hum Retroviruses 2004, 20: 1210-1222. 10.1089/aid.2004.20.1210PubMedView ArticleGoogle Scholar
- Kim SY, Volsky DJ: PAGE: parametric analysis of gene set enrichment. BMC Bioinforma 2005, 6: 144. 10.1186/1471-2105-6-144View ArticleGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 2005, 102: 15545-15550. 10.1073/pnas.0506580102PubMedPubMed CentralView ArticleGoogle Scholar
- Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, et al.: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res 2003, 31: 374-378. 10.1093/nar/gkg108PubMedPubMed CentralView ArticleGoogle Scholar
- Abel PM, McSharry C, Galloway E, Ross C, Severn A, Toner G, Gruer L, Wilkinson PC: Heterogeneity of peripheral blood monocyte populations in human immunodeficiency virus-1 seropositive patients. FEMS Microbiol Immunol 1992, 5: 317-323.PubMedView ArticleGoogle Scholar
- Koethe SM, Carrigan DR, Turner PA: Increased density of HLA-DR antigen on monocytes of patients infected with the human immunodeficiency virus. J Med Virol 1989, 29: 82-87. 10.1002/jmv.1890290203PubMedView ArticleGoogle Scholar
- Locher C, Vanham G, Kestens L, Kruger M, Ceuppens JL, Vingerhoets J, Gigase P: Expression patterns of Fc gamma receptors, HLA-DR and selected adhesion molecules on monocytes from normal and HIV-infected individuals. Clin Exp Immunol 1994, 98: 115-122.PubMedPubMed CentralView ArticleGoogle Scholar
- Joshi P, Sloan B, Torbett BE, Stoddart CA: Heat shock protein 90AB1 and hyperthermia rescue infectivity of HIV with defective cores. Virology 2013, 436: 162-172. 10.1016/j.virol.2012.11.005PubMedPubMed CentralView ArticleGoogle Scholar
- Huang X, Seifert U, Salzmann U, Henklein P, Preissner R, Henke W, Sijts AJ, Kloetzel PM, Dubiel W: The RTP site shared by the HIV-1 Tat protein and the 11S regulator subunit alpha is crucial for their effects on proteasome function including antigen processing. J Mol Biol 2002, 323: 771-782. 10.1016/S0022-2836(02)00998-1PubMedView ArticleGoogle Scholar
- Apps R, Qi Y, Carlson JM, Chen H, Gao X, Thomas R, Yuki Y, Del Prete GQ, Goulder P, Brumme ZL, et al.: Influence of HLA-C expression level on HIV control. Science 2013, 340: 87-91. 10.1126/science.1232685PubMedPubMed CentralView ArticleGoogle Scholar
- Favre D, Stoddart CA, Emu B, Hoh R, Martin JN, Hecht FM, Deeks SG, McCune JM: HIV disease progression correlates with the generation of dysfunctional naive CD8(low) T cells. Blood 2011, 117: 2189-2199. 10.1182/blood-2010-06-288035PubMedPubMed CentralView ArticleGoogle Scholar
- Tran HT, Van den Bergh R, Loembe MM, Worodria W, Mayanja-Kizza H, Colebunders R, Mascart F, Stordeur P, Kestens L, De Baetselier P, Raes G: Modulation of the complement system in monocytes contributes to tuberculosis-associated immune reconstitution inflammatory syndrome. AIDS 2013, 27: 1725-1734. 10.1097/QAD.0b013e328361648bPubMedView ArticleGoogle Scholar
- Carroll MC: The complement system in regulation of adaptive immunity. Nat Immunol 2004, 5: 981-986. 10.1038/ni1113PubMedView ArticleGoogle Scholar
- Baldwin GC, Fleischmann J, Chung Y, Koyanagi Y, Chen IS, Golde DW: Human immunodeficiency virus causes mononuclear phagocyte dysfunction. Proc Natl Acad Sci USA 1990, 87: 3933-3937. 10.1073/pnas.87.10.3933PubMedPubMed CentralView ArticleGoogle Scholar
- Collini P, Noursadeghi M, Sabroe I, Miller RF, Dockrell DH: Monocyte and macrophage dysfunction as a cause of HIV-1 induced dysfunction of innate immunity. Curr Mol Med 2010, 10: 727-740. 10.2174/156652410793384141PubMedView ArticleGoogle Scholar
- Michailidis C, Giannopoulos G, Vigklis V, Armenis K, Tsakris A, Gargalianos P: Impaired phagocytosis among patients infected by the human immunodeficiency virus: implication for a role of highly active anti-retroviral therapy. Clin Exp Immunol 2012, 167: 499-504. 10.1111/j.1365-2249.2011.04526.xPubMedPubMed CentralView ArticleGoogle Scholar
- Dugast AS, Tonelli A, Berger CT, Ackerman ME, Sciaranghella G, Liu Q, Sips M, Toth I, Piechocka-Trocha A, Ghebremichael M, Alter G: Decreased Fc receptor expression on innate immune cells is associated with impaired antibody-mediated cellular phagocytic activity in chronically HIV-1 infected individuals. Virology 2011, 415: 160-167. 10.1016/j.virol.2011.03.012PubMedPubMed CentralView ArticleGoogle Scholar
- Van Grol J, Subauste C, Andrade RM, Fujinaga K, Nelson J, Subauste CS: HIV-1 inhibits autophagy in bystander macrophage/monocytic cells through Src-Akt and STAT3. PLoS ONE 2010, 5: e11733. 10.1371/journal.pone.0011733PubMedPubMed CentralView ArticleGoogle Scholar
- Chugh P, Bradel-Tretheway B, Monteiro-Filho CM, Planelles V, Maggirwar SB, Dewhurst S, Kim B: Akt inhibitors as an HIV-1 infected macrophage-specific anti-viral therapy. Retrovirology 2008, 5: 11. 10.1186/1742-4690-5-11PubMedPubMed CentralView ArticleGoogle Scholar
- Sirois M, Robitaille L, Sasik R, Estaquier J, Fortin J, Corbeil J: R5 and ×4 HIV viruses differentially modulate host gene expression in resting CD4+ T cells. AIDS Res Hum Retroviruses 2008, 24: 485-493. 10.1089/aid.2007.0120PubMedView ArticleGoogle Scholar
- Komano J, Miyauchi K, Matsuda Z, Yamamoto N: Inhibiting the Arp2/3 complex limits infection of both intracellular mature vaccinia virus and primate lentiviruses. Mol Biol Cell 2004, 15: 5197-5207. 10.1091/mbc.E04-04-0279PubMedPubMed CentralView ArticleGoogle Scholar
- Vilhardt F, Plastre O, Sawada M, Suzuki K, Wiznerowicz M, Kiyokawa E, Trono D, Krause KH: The HIV-1 Nef protein and phagocyte NADPH oxidase activation. J Biol Chem 2002, 277: 42136-42143. 10.1074/jbc.M200862200PubMedView ArticleGoogle Scholar
- Ishaq M, Lin BR, Bosche M, Zheng X, Yang J, Huang D, Lempicki RA, Natarajan V: LIM kinase 1 - dependent cofilin 1 pathway and actin dynamics mediate nuclear retinoid receptor function in T lymphocytes. BMC Mol Biol 2011, 12: 41. 10.1186/1471-2199-12-41PubMedPubMed CentralView ArticleGoogle Scholar
- Vicenzi E, Alfano M, Ghezzi S, Gatti A, Veglia F, Lazzarin A, Sozzani S, Mantovani A, Poli G: Divergent regulation of HIV-1 replication in PBMC of infected individuals by CC chemokines: suppression by RANTES, MIP-1alpha, and MCP-3, and enhancement by MCP-1. J Leukoc Biol 2000, 68: 405-412.PubMedGoogle Scholar
- Yao H, Yang Y, Kim KJ, Bethel-Brown C, Gong N, Funa K, Gendelman HE, Su TP, Wang JQ, Buch S: Molecular mechanisms involving sigma receptor-mediated induction of MCP-1: implication for increased monocyte transmigration. Blood 2010, 115: 4951-4962. 10.1182/blood-2010-01-266221PubMedPubMed CentralView ArticleGoogle Scholar
- Lyons PA, Koukoulaki M, Hatton A, Doggett K, Woffendin HB, Chaudhry AN, Smith KG: Microarray analysis of human leucocyte subsets: the advantages of positive selection and rapid purification. BMC Genomics 2007, 8: 64. 10.1186/1471-2164-8-64PubMedPubMed CentralView ArticleGoogle Scholar
- Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics 2008, 24: 1547-1548. 10.1093/bioinformatics/btn224PubMedView ArticleGoogle Scholar
- Lin SM, Du P, Huber W, Kibbe WA: Model-based variance-stabilizing transformation for Illumina microarray data. Nucleic Acids Res 2008, 36: e11.PubMedPubMed CentralView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 2004, 3: Article3.PubMedGoogle Scholar
- Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: database for annotation, visualization, and integrated discovery. Genome Biol 2003, 4: P3. 10.1186/gb-2003-4-5-p3PubMedView ArticleGoogle Scholar
- Fu W, Sanders-Beer BE, Katz KS, Maglott DR, Pruitt KD, Ptak RG: Human immunodeficiency virus type 1, human protein interaction database at NCBI. Nucleic Acids Res 2009, 37: D417-D422. 10.1093/nar/gkn708PubMedPubMed CentralView ArticleGoogle Scholar
- Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3’ UTRs by comparison of several mammals. Nature 2005, 434: 338-345. 10.1038/nature03441PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.