Abstract Background Human central memory CD4 T cells are characterized by their capacity of proliferation and differentiation into effector memory CD4 T cells. Homeostasis of central memory CD4 T cells is considered a key factor sustaining the asymptomatic stage of Human Immunodeficiency Virus type 1 (HIV-1) infection, while progression to acquired immunodeficiency syndrome is imputed to central memory CD4 T cells homeostatic failure. We investigated if central memory CD4 T cells from patients with HIV-1 infection have a gene expression profile impeding proliferation and survival, despite their activated state. Methods Using gene expression microarrays, we analyzed mRNA expression patterns in naive, central memory, and effector memory CD4 T cells from healthy controls, and naive and central memory CD4 T cells from patients with HIV-1 infection. Differentially expressed genes, defined by Log[2] Fold Change (FC) ≥ |0.5| and Log (odds) > 0, were used in pathway enrichment analyses. Results Central memory CD4 T cells from patients and controls showed comparable expression of differentiation-related genes, ruling out an effector-like differentiation of central memory CD4 T cells in HIV infection. However, 210 genes were differentially expressed in central memory CD4 T cells from patients compared with those from controls. Expression of 75 of these genes was validated by semi quantitative RT-PCR, and independently reproduced enrichment results from this gene expression signature. The results of functional enrichment analysis indicated movement to cell cycle phases G1 and S (increased CCNE1, MKI67, IL12RB2, ADAM9, decreased FGF9, etc.), but also arrest in G2/M (increased CHK1, RBBP8, KIF11, etc.). Unexpectedly, the results also suggested decreased apoptosis (increased CSTA, NFKBIA, decreased RNASEL, etc.). Results also suggested increased IL-1β, IFN-γ, TNF, and RANTES (CCR5) activity upstream of the central memory CD4 T cells signature, consistent with the demonstrated milieu in HIV infection. Conclusions Our findings support a model where progressive loss of central memory CD4 T cells in chronic HIV-1 infection is driven by increased cell cycle entry followed by mitotic arrest, leading to a non-apoptotic death pathway without actual proliferation, possibly contributing to increased turnover. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3308-8) contains supplementary material, which is available to authorized users. Keywords: HIV, Immunologic Memory, Cell Cycle, Cell Death, CD4-Positive T-Lymphocytes, Transcriptome, Homeostasis Background Acute HIV infection depletes mucosal CD4 T cells, mainly effector memory (T[EM]) cells, rapidly and profoundly [[45]1–[46]3]. The ensuing chronic phase is largely asymptomatic, even though mucosal tissues are not replenished with T[EM] cells [[47]4]. Simian immunodeficiency virus (SIV) infection of rhesus macaques (an animal model of human HIV disease) shows that opportunistic control infection in the chronic phase is mediated by remnant mucosal T[EM] cells supplied by the differentiation of central memory (T[CM]) cells in lymph nodes [[48]5, [49]6]. Additionally, human T[CM] cells are also characterized by their capacity of proliferation and differentiation into T[EM] cells [[50]7, [51]8]. Thus, homeostasis of T[CM] cells is considered a key factor sustaining the asymptomatic stage of HIV infection, while progression to acquired immunodeficiency syndrome is attributed to homeostatic failure of T[CM] cells [[52]5, [53]6, [54]9–[55]12]. It is unclear how this homeostatic equilibrium is lost during chronic infection. CD4 T cell maturation subpopulations (T[N], T[CM], and T[EM]) [[56]7] are differentially affected by HIV infection [[57]13, [58]14]; with T[EM] cells being HIV’s main target [[59]15]. T[CM] cells can be infected in a lower proportion by HIV, which has led to propose that direct virion-mediated cytopathicity could gradually eliminate them, leading to poor homeostatic activity [[60]6]. Nevertheless, direct cytopathicity by HIV [[61]16] cannot completely explain CD4 T cell depletion during chronic infection [[62]17–[63]20], which suggests the participation of indirect pathogenic mechanisms, particularly chronic activation [[64]12, [65]21]. Additionally, CD4 T cells from patients with HIV could be intrinsically altered, as suggested by the limited proportion of HIV-infected patients recovering their pre-infection CD4 T cells counts under virus-controlling antiretroviral therapy [[66]22]. In this regard, we have found intrinsic dysfunctions in activated T[CM] cells from HIV-infected patients, as a lowered IL-2 response and CD40L induction after T cell receptor (TCR)-mediated stimulation [[67]23, [68]24], which could decrease their proliferative, differentiation, and survival capacities. In order to determine if circulating T[CM] cells from HIV-infected patients have a transcriptome consistent with activation, but simultaneously with altered capacities to divide and survive, we compared the ex-vivo messenger Ribonucleic acid (mRNA) whole-genome expression patterns of CD4 T naive (T[N]) and T[CM] cells from HIV^+ patients with T[N], T[CM], and T[EM] cells from healthy controls. We found a T[CM] cell signature in HIV-1 infection suggesting that the loss of this subpopulation may be driven by increased cell cycle entry followed by mitotic arrest possibly leading to cell death in a non-senescent or effector-like state. Methods Participants This study was approved by the boards of Instituto Nacional de Enfermedades Respiratorias Ismael Cosío Villegas (reference number B29-11), and Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán (reference number 1403). All patients signed written informed consent according with the Helsinki Protocol. Blood samples were obtained from 9 HIV¯ controls, and 6 HIV^+ patients. Patients had median 480 CD4 T cells/μL blood (range 330–757), and median 121 563 HIV-ribonucleic acid (RNA) copies/mL-blood (23 883–41 2584). Among them, patients providing T[CM] cells had viral loads of 23 883, 81 834 and 107 732 HIV RNA copies/mL-blood, and CD4 T cell counts of 439, 473 and 491 CD4 T cells/μL blood, respectively. Relative telomere length was determined in samples from ten additional HIV¯ controls, and ten additional HIV^+patients with median 628 CD4 T cells/ μL-blood (194–1 128) and median 485 882 HIV-RNA copies/mL-blood (3 870–3 500 000). Patients were antiretroviral therapy-naive, free of opportunistic infections and malignancies, and were not taking any immunomodulatory drugs. Isolation of CD4 T cell subpopulations Peripheral blood mononuclear cells (PBMCs) were purified from 50 to 60 mL of peripheral blood by sedimentation on Lymphoprep (Fresenius Kabi Norge, Oslo, Norway). CD4 T[N] (CD45RA^+ CCR7^+), T[CM] (CD45RA¯ CCR7^+) and T[EM] (CD45RA¯ CCR7¯) cells were purified from PBMCs using immunomagnetic beads (Miltenyi Biotec, Bergisch Gladbach, Germany). Subpopulation purity was determined according to the expression of CD4, CD45RA and CCR7, using anti-CD4-APC-Cy7, anti-CD45RA-APC (BD Biosciences, San José, CA, USA), and anti-CCR7-PE (Miltenyi Biotec) fluorochrome-conjugated antibodies (See Additional file [69]1). Cells were analyzed in a FACSCanto II flow cytometer (BD Biosciences). Cells with purity >90% were used. Membrane CD38 was detected with an anti-CD38-biotin (Miltenyi Biotec) plus streptavidin PerCp-Cy5.5 (Biolegend, San Diego, CA, USA). RNA extraction and microarray analysis Total RNA was obtained from three T[N], three T[CM], and three T[EM] CD4 T cell samples from healthy controls, and three T[N] and three T[CM]CD4 T cell samples from HIV^+ patients, using RNeasy Mini Kit (Qiagen, Venlo, Netherlands). Each RNA sample proceeded from a different subject. Scarcity of patients’ T[EM] cells precluded obtaining sufficient RNA. Microarray gene expression analysis used equimolar concentrations of total RNA from T cell subpopulations. Complementary deoxyribonucleic acid (cDNA) synthesis, amplification, and gene expression profiling were performed according to the manufacturer’s instructions (Affymetrix WT Sense Target labeling assay manual, California, USA). Labeled DNA was added to hybridization cocktail and injected into the array (GeneChip Human Gene 1.0 ST Array, Affymetrix). Washing and staining steps were performed in the GeneChip Fluidics Station 450 (Affymetrix). Probe arrays were scanned using a GeneChip Scanner 3000 7G (Affymetrix). Data were deposited in GEO, series record [70]GSE73968. Background correction and normalization were performed with Robust Multiarray Average Method (RMA) [[71]25] using Bioconductor package [[72]26] of R [[73]27]. A Principal component analysis (PCA) of normalized signals from all genes in each microarray was performed using R [[74]27]. Modeling gene expression was performed using linear models of Limma package [[75]28]. The B-statistic was used as significant measure to define differentially expressed genes. This statistic is computed as the posterior odds of differential expression. It is reformulated in terms of a moderated t-statistic in which posterior residual standard deviations are used in place of ordinary standard deviations. Essentially, the B-statistic compromises between individual gene variance estimates and a single variance estimate for all genes. The probabilities are transformed to a scale that goes from –Inf to Inf using log odds. The B-statistic is analogous to the adjusted p-value, which addresses statistical significance for multiple comparisons. Here, genes with Log 2 Fold Change (FC) ≥ |0.5| and Log odds > 0 were considered as differentially expressed. Limma statistics such as adjusted p-value and the B statistic can be seen in Additional file [76]2. FDR Benjamini Hochberg multiple testing correction [[77]29] was applied to control the number of false positives. Both B statistic and adjusted p-value showed consistency across differentially expressed genes. Unsupervised 2-way hierarchical clustering analysis of gene expression data was performed using Euclidian distance and average linkage with gplots [[78]30] of R [[79]27]. Venn diagrams were made with Venny 2.0.2 [[80]31]. Functional enrichment analyses were performed with Data Base for Annotation, Visualization and Integrated Discovery (DAVID) [[81]32, [82]33], Gen Set Enrichment Analysis (GSEA) [[83]34] and Ingenuity Pathway Analysis (IPA, QIAGEN Redwood City, CA, USA). DAVID uses a Fisher Exact test in order to determine gene-enrichment in annotation terms. A gene set is enriched when the proportion of genes in a list that falls into an annotation term differs from the background model. The EASE score is a modified Fisher exact p-value. Basically, if n is the number of genes in the list that falls into a given annotation term, n-1 is used to compute the p-value [[84]32, [85]33]. Gene set enrichment methods also implement strategies for addressing the issue of multiple testing hypotheses. GSEA uses a ranking procedure to produce a gene list from the full expression matrix. This is done by computing an Enrichment Score (ES(S)). It controls the ratio of false positives to the total number of gene sets attaining a fixed level of significance using FDR [[86]34]. IPA assesses enrichment (i. e. biological functions that could be increased or decreased given the observed gene expression patterns) using a Fisher exact p-value. Additionally, it computes a Z score that allows inferring upstream transcriptional regulators and expectable enriched functions, based on statistical significance by comparing the match between observed and predicted up/down regulation patterns. The null model is referred as activation Z-score [[87]35]. Predicted regulation patterns are based on previously reported causal relationships between relevant genes and functions [[88]35]. Semi-quantitative real-time PCR We used B2M, GAPDH, POLR2A, and TBP as reference genes to normalize expression. RNA proceeded from the samples used for microarray analysis. cDNA was synthesized from ~100 ng total RNA with Transcriptor First Strand cDNA Synthesis Kit (Roche Applied Science, Mannheim, Germany), using random hexamers and performing one cycle of 10 min 25 °C; 30 min 55 °C, and 5 min 85 °C. cDNA was stored at −20 °C until use. PCR amplifications were performed by high-throughput gene expression analysis using DNA binding dye Evagreen (SsoFast MasterMix, Biorad, California, USA) for product detection, and specific primers for each gene (DELTAgene Assays, Fluidigm Corporation, California, USA). Specific target pre-amplification of each cDNA and a cleanup step were performed as described elsewhere [[89]36]. We performed semiquantitative RT-PCR using Fast Gene Expression Analysis with EvaGreen (Biorad), following the Biomark System Protocol (Fluidigm Corporation, California, USA). Assay mixes (100 μM of each pair of primers, 2X Assay Loading Reagent, and 1X TE buffer), sample mixes (pre-amplified cDNA, 2X SsoFast MasterMix (BioRad), and 20X DNA Binding Dye Sample Loading Reagent (Fluidigm), were loaded into a 96.96 Dynamic Array (Fluidigm), using the IFC Controller HX (Fluidigm), and were then transferred to a BioMark HD device (Fluidigm) for the PCR cycles (40 min 70 °C, 30s 60 °C; 60s 95 °C, then 30 cycles of 5 s 96 °C, 60s 60 °C). Melting curves were determined at the 60 to 95 °C rise, with a temperature change rate of 1 °C/3 s. Ct values were obtained with Fluidigm Real-Time PCR Analysis Version 4.1.3 software (Fluidigm). Only Ct values <30 and amplicons with only 1 melting curve were used. Geometric means of four reference genes were used to normalize expression data [[90]37]. Relative expression was calculated as ΔΔCt. Expression of each gene was determined with six technical replicates per sample. Normality was verified using Kolmogorov-Smirnov test, which quantifies the distance between the empirical distribution of the sample and the cumulative distribution of the reference distribution, which in this case is assumed to be normal. Group differences were analyzed with Student’s t test. Data management and statistics were done with Reshape [[91]38] and fBasics [[92]39] packages of R [[93]27]. Relative telomere length Telomere PNA kit/FITC (Dako, São Paulo, Brazil) was used following the manufacturer’s instructions, including thymocytes from 6-week old mice as reference for normalization. Briefly, samples were prepared by mixing 10^6 mouse thymocytes and 10^6 T[CM] cells. The mixture was distributed into four tubes. 150 μl of FITC-labeled peptide nucleic acid (PNA) probe solution was added into two tubes while 150 μl of unlabeled PNA probe solution was added into the other two. Samples were hybridized in a pre-warmed heating block (TB2 Thermoblock, Biometra, Göttingen, Germany) set at 82 °C, 10 min, and left overnight at room temperature. Samples were washed twice. Between washing steps, samples were heated to 40 °C in a pre-warmed TB2 Thermoblock (Biometra) for 10 min. Samples were resuspended in 250 μL of DNA staining solution (1X), and stored overnight at 4 °C, away of light. Then, samples were analyzed by flow cytometry in a FACSCanto II (BD Biosciences). Results T[CM] cells from HIV^+ patients are not more differentiated and are not more senescent Unsupervised principal component analysis of normalized whole genome expression data segregated samples of each maturation subpopulation, and further separated samples originating from persons with HIV and samples from controls (Fig. [94]1a). Thus, phenotype-based classification of differentiation subpopulations [[95]7, [96]8] reliably reflected distinct gene expression programs, as previously reported [[97]40–[98]42], which were altered by HIV infection. Fig. 1. Fig. 1 [99]Open in a new tab Gene expression does not support greater differentiation or senescence of T[CM] cells from HIV^+ patients. a Principal component analysis of the entire microarray data set of each subpopulations from HIV^+ patients (red triangles and red circles) and controls (blue triangles, blue circles and blue squares). The first three principal components are shown, accounting for 55% of variance in a three dimensional plot. b Heat map resulting from hierarchical clustering of genes related with normal differentiation (pairwise comparisons between not infected subpopulations). Each row represents a differentially expressed gene. Each column represents each independent sample. The unsupervised two-way hierarchical clustering is shown as a dendrogram for genes (left), and a dendrogram for samples (top). In the upper dendrogram (samples) the independent resulting nodes, each one corresponding to a maturation subpopulation, is encircled in green. c Sequential downregulation of selected naïve-associated genes, and sequential upregulation of selected effector-associated genes when samples are arranged according to the linear differentiation model. Data are represented as means ± 1 SEM of three donors (blue) and three patients (red). TBX21 expression difference between T[CM] and T[CM] HIV was analyzed with Student’s t test. d Relative telomere length of central memory CD4 T cells (T[CM]) from HIV^+ patients (red triangles) and controls (blue circles), Student’s t-test was used to compare groups. We were unable to obtain sufficient RNA from T[EM] cells from patients due to their small number We asked if differential gene expression by patients’ and controls’ T[CM] cells reflected greater differentiation of patients’ cells (towards effector stages) [[100]23, [101]24]. Using the criteria defined in methods (Log[2]FC ≥ |0.5| and Log (odds) > 0) we looked in the whole transcriptome for all differentially expressed genes in the following pair-wise comparisons of CD4 T cell subpopulations from controls: T[CM] vs. T[N], T[EM] vs. T[CM,] and T[EM] vs. T[N] (arrows a, b and c in Fig. [102]2a). The resulting 1858 differentially expressed genes are subsequently referred to as differentiation-related genes (corresponding to subpopulations in distinct stages of differentiation). We performed an unsupervised 2-way hierarchical clustering analysis of these 1858 differentiation-related genes (Fig. [103]1b, and Additional file [104]3). T[N] and T[CM] cells from patients grouped with their counterparts from controls (Fig. [105]1b). Samples of a same subpopulation were assigned to a same node (green circles 1, 2, and 3 on Fig. [106]1b), regardless of their HIV status. The expression of the differentiation-related genes progressively decreased or increased in the order of linear differentiation (T[N] → T[CM] → T[EM] ), agreeing with previous reports [[107]40–[108]43] (Fig. [109]1c). For instance, LEF1, ACTN1, FOXP1, IL6ST and CERS6 reportedly undergoing down-regulation in naive T cells after antigen recognition and differentiation [[110]44–[111]48], along with TAF4B, appeared progressively down regulated when samples were ordered according to the linear model of peripheral differentiation (Fig. [112]1c, left panel). These changes agree with previous reports [[113]42]. Conversely, differentiation and effector function-associated transcripts, like EOMES, TBX21 (t-bet), PRDM1 (Blimp-1) [[114]49, [115]50], GZMA and PRF1 [[116]51, [117]52], were gradually increased in the same order (Fig. [118]1c right panel). A same pattern was followed by the expression of KLRG1,an indicator of replicative senescence [[119]53, [120]54] (Fig. [121]1c). TBX21 (t-bet) was the only gene with increased expression in T[CM] cells from patients, compared with controls (p = 0.003), which, along with the increased expression of IL12R e IL18R, suggests a Th1-skewed response driven by HIV infection. A Th1-skewed response was also predicted by Ingenuity Canonical Pathway analysis (See Additional file [122]4). Thus, T[CM] cells from patients did not seem to be more differentiated than their counterparts from controls, but appeared polarized to Th1 functions. Fig. 2. Fig. 2 [123]Open in a new tab Unique T[CM] cell signature in HIV infection. Differential expression was defined as Log[2] of fold change (Log [FC]) ≥ |0.5|, and Log (odds) > 0. a Pairwise comparisons of samples of CD4 T cells subpopulations from HIV^+ and HIV¯ groups indicated by arrows a and d. Number of genes differentially expressed in each comparison are shown. Blue circles, controls’ samples; red squares, HIV^+ patients’ samples. b Venn diagram of sets of differentially expressed genes. Each pairwise comparison is depicted by a colored oval. The number of differentially expressed genes found in more than one comparison appear in the intersections. c Heat map displaying a two-way unsupervised hierarchical clustering of 210 differentially expressed gens distinguishing HIV^+ patients’ T[CM] cells (red bar) and controls’ T[CM]cells (blue bar), grouped in dendrograms. Each column represents an independent sample (biological replica) of each subpopulation, numbered 1 to 3. Each row corresponds to a differentially expressed gene We then asked if patients’ T[CM] cells had a longer replicative history, which would entail a shortening of telomeres. We did not find any difference in relative telomere length between T[CM] cells from patients and controls (p = 0.737, Fig. [124]1d), agreeing with KLRG1 expression [[125]53, [126]54], and suggesting that they are not in a more senescent state. T[CM] gene expression signature in HIV infection Having ruled out a greater differentiation of patients’ T[CM] cells, we investigated if the gene expression signature of these cells revealed a functional state that could explain loss of homeostatic capacity. Using the criteria defined in methods (Log[2]FC ≥ |0.5| and Log (odds) > 0), we looked for genes that were differentially expressed by T[CM] cells from HIV^+ patient and T[CM] from controls. We found a total of 210 differentially expressed genes. We refer to this 210- gene list as the gene expression signature of T[CM] cells in HIV infection (See Fig. [127]2a arrow d, b red oval and Additional file [128]2). This gene expression signature was obtained from the transcriptome independently of the list of 1858 differentiation-related genes. Among these 210 differentially expressed genes, 137 were absent in all other pairwise comparison (Fig. [129]2a, b). Hierarchical clustering analysis showed clear and consistent differences in the relative expression of these 210 genes between patients and controls (Fig. [130]2c). Of note, biological replicates were very homogeneous. We analyzed the HIV T[CM] signature with the enrichment analysis tools IPA, GSEA, and DAVID. These different analyses consistently yielded four general functional categories that were modified in T[CM] cells from patients: cell cycle, DNA damage and repair, apoptosis, and immune responses (Table [131]1). Notably, the 137 genes uniquely distinguishing T[CM] cells from patients and controls (Fig. [132]2b red oval) sufficed to yield the same four functional categories when analyzed with DAVID and IPA. This suggests that the enriched functions largely depend on the T[CM] cell signature. GSEA rendered a larger set of altered immune functions, likely because it uses data from the entire microarray, and because it detects more modest changes when the members of a function or pathway show a strong correlation [[133]34] (Table [134]1). GSEA identified Toll-like receptors (TLR), type I interferons, IL1, and NLRs signaling, plus NFκB activation, all of them related to an inflammatory milieu. Independently, Ingenuity Upstream Regulator Analysis [[135]35] assigned the greatest z-scores and the most significant p-values to the activity of IL-1B, TNF, NFκB complex, and CCL5, as possible upstream molecules eliciting the expression changes constituting the T[CM] signature (See Additional file [136]5, upstream analysis). Cell cycle, DNA damage and repair, and apoptosis (greatly related functions [[137]55, [138]56]) appeared consistently in the output of all enrichment analysis tools (Table [139]1). There were several functional categories closely related with G0/G1/S transition and G2/M checkpoints in the output of GSEA analysis. IPA, which weighs its predictions, displayed increased proliferation and cell survival, decreased apoptosis and decreased cell death. In contrast, IPA’s output simultaneously indicated an increase in cytostasis, movement to interphase, and a decrease in mitosis (Table [140]1). Table 1. Enriched categories of functions according to TCM gene expression signature in HIV infection General category DAVID GSEA IPA EASE < 0.05 Number of genes FDR < 0.05, p < 0.001 p < 0.01 Number of genes Prediction sense Cell cycle Cell cycle 15 DNA replication Proliferation of tumor cell lines 31 Positive Cell division 10 Cell cycle Proliferation of cells 53 Positive mitosis 9 Mitotic M/G1 Cytostasis 6 Positive G1/S transition Cytostasis od tumor cell lines 5 Positive Cell cycle check points Interphase of tumor cell lines 11 Positive Cyclin E associated event during G1/S transition Cell survival 28 Positive Assembly of pre-replicative complex Interphase 13 Positive G0 and early G1 Mitosis 10 Negative G2/M check points DNA damage or repair p53 signaling pathway 5 p53 dependent G1 DNA damage response ATM signal pathway 3 Cell cycle checkpoints 5 Apoptosis Apoptosis 13 Apoptosis of tumor cell lines 32 Negative Cell death of cancer cells 6 Negative Apoptosis of cervical cancer cell lines 10 Negative Immune responses Toll endogenous pathway Synthesis of reactive oxygen species 7 Positive IL1 signaling IFN-alpha/beta signaling Chemokine receptors bind chemokines NOD like receptors signaling NFKB activation by IKKS complex Myd88 cascade TLR4 signaling IL12 pathway FOXO Pathway [141]Open in a new tab Enriched categories of functions according to differential expression of 210 genes in T[CM] cells from HIV^+ patients and controls. DAVID and IPA tools show in a column the number of genes supporting each prediction. EASE Score is the P-Value of a modified Fisher Exact test of the significance of gene enrichment in a gen-set. FDR: False discovery rate. In IPA, the sign indicates if the function would be up-regulated (positive) or down-regulated (negative) We re-analyzed mRNA expression by RT-PCR of 91 genes of the HIV T[CM] signature that were associated with enriched functions, and B2M, GAPDH, POLR2A, and TBP as reference genes. This analysis validated 75 genes (82%) of the signature (See Additional file [142]6). Five genes failed amplification, and 11 were not differentially expressed when assessed by RT-PCR. The expression of reference genes did not differ between samples (See Additional file [143]7). An analysis with IPA using only the 75 validated genes yielded the same enriched functions and pathways as microarray data (Fig. [144]3a, b, c, d, e, and Additional file [145]8). While some genes were related with more than one function, many were related exclusively with one function (Fig. [146]3f), supporting an unambiguous prediction. Finally, analyzing the T[CM] cells used for microarray and RT-PCR analyses, we found that surface expression of the CD38 protein was more frequent among T[CM] cells from patients than among those from controls, consistent with mRNA results (Fig. [147]3g). Fig. 3. Fig. 3 [148]Open in a new tab Cell cycle impairments are predicted using RT-PCR-validated genes from the signature. Each graph (a to e) represents the group of validated genes by RT-PCR associated with an increased enriched function (red title and positive z-score) or decreased enriched function (blue title and negative z-core). f Venn diagram depicting the differentially expressed genes within sets corresponding to three IPA predictions. Intersections correspond to genes appearing in more than one prediction. g Frequency of surface expression of CD38 on T[CM] cells from two groups, corresponding to gene expression results A model of T[CM] cell death in HIV infection Since the predictions of increased proliferation and increased cytostasis were incompatible, and the prediction of reduced apoptosis did not agree with previous evidence [[149]13, [150]57–[151]60], we took into account that enrichment tools base their predictions on a broad set of previous findings, ranging from very particular to very general ones. Accordingly, we investigated if the predictions were based on more demarcated processes, and if these processes were compatible. With this purpose, we reviewed the references supporting