|Home | About | Journals | Submit | Contact Us | Français|
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 molecular mechanisms of cell cycle exit are poorly understood. Studies on lymphocytes at cell cycle exit after growth factor deprivation have predominantly focused on the initiation of apoptosis. We aimed to study gene expression profile of primary and immortalised IL-2-dependent human T cells forced to exit the cell cycle by growth factor withdrawal, before apoptosis could be evidenced.
By the Affymetrix microarrays HG-U133 2.0 Plus, 53 genes were distinguished as differentially expressed before and soon after IL-2 deprivation. Among those, PIM1, BCL2, IL-8, HBEGF, DUSP6, OSM, CISH, SOCS2, SOCS3, LIF and IL13 were down-regulated and RPS24, SQSTM1, TMEM1, LRRC8D, ECOP, YY1AP1, C1orf63, ASAH1, SLC25A46 and MIA3 were up-regulated. Genes linked to transcription, cell cycle, cell growth, proliferation and differentiation, cell adhesion, and immune functions were found to be overrepresented within the set of the differentially expressed genes.
Cell cycle exit of the growth factor-deprived T lymphocytes is characterised by a signature of differentially expressed genes. A coordinate repression of a set of genes known to be induced during T cell activation is observed. However, growth arrest following exit from the cell cycle is actively controlled by several up-regulated genes that enforce the non-dividing state. The identification of genes involved in cell cycle exit and quiescence provides new hints for further studies on the molecular mechanisms regulating the non-dividing state of a cell, the mechanisms closely related to cancer development and to many biological processes.
Cell growth, proliferation and arrest mechanisms are crucial for development and homeostasis. Extensive studies have examined the mechanisms of cell activation and cell cycle progression, which involve sequential effects of growth factors and their receptors and activation of intracellular signal transduction pathways, transcription factors, cyclins and cyclin-dependent kinases . Studies on lymphocytes withdrawn from the cell cycle have focused on the initiation of apoptosis following growth factor deprivation [2-8] or on the cell-to-cell contact mechanisms that prevent cells from apopotosis [9-11]. Still, the molecular mechanisms of cell cycle exit remain poorly understood. With the lack of mitogenic signals, and in response to antiproliferative signals such as mitogen withdrawal, growth factor starvation, contact inhibition or DNA damage, a proliferating cell arrests to become quiescent (a reversible, nondividing state) or senescent, or to undergo apoptotic death. Most somatic cells in an adult body remain in a postmitotic G0 phase. However, during tissue regeneration and repair, wound healing and immune response, cells re-enter the cell cycle and subsequently withdraw from proliferation.
Studies in lymphocytes and fibroblasts have shown that resting cells are not in a passive state resulting simply from the lack of stimulation. Pajalunga et al.  summarized cell cycle exit as "a shift in the balance between positive and negative regulators of proliferation in favour of the latter". Indeed, resting cells were found to express sets of up-regulated and down-regulated genes that maintained an active state of non-division [1,12,13,15,16]. A group of genes required for cell cycle exit and the maintenance of cell quiescence in human fibroblasts following serum deprivation has been recently identified . Comparative studies of transcriptional profiles of resting and stimulated lymphocytes [1,18] have shown that T lymphocyte quiescence is actively maintained by products of a set of genes, including TOB and KLF (LKLF, GKLF and BKLF) and FOXO families of transcription factors. Other candidate genes involved in lymphocyte quiescence, including TSC-22 and Dyrk1, have been identified in a mouse model among genes highly expressed in resting lymphocytes and down-regulated after T or B cell activation .
The present study on a human T cell model may provide new insights into the molecular changes involved in cell proliferation control, cell cycle exit and cellular quiescence. Cellular quiescence represents an important safeguard against tumorigenesis; therefore the identification of quiescence-controlling genes and a better understanding of the relevant molecular pathways may open new possibilities in cancer diagnosis and treatment. Thus, in our study, proliferating interleukin (IL)-2-dependent human T cells were forced to exit cell cycle by growth factor withdrawal, and their gene expression profiles were examined.
To gain new insights into cell cycle exit of T lymphocytes, we performed gene expression profiling in human IL-2-dependent T lymphoblasts deprived of the growth factor. Following IL-2 withdrawal, IL-2-dependent, proliferating T lymphoblasts cease to divide and undergo apoptosis . Therefore, transcription profiles of T cells were analysed soon after IL-2 withdrawal (8 hours), before apoptotic changes could be observed. In addition, to minimize effects associated with the heterogeneity of primary T cell populations, we also studied clonal populations of immortalised IL-2-dependent T lymphoblasts and looked for the common changes in the two cellular models. We believe that this approach also limited the "non-specific" gene expression changes in the primary IL-2-deprived T lymphoblasts, i.e. not related to IL-2 withdrawal, but rather to the possible paracrine signalling in polyclonal populations by other cytokines characterised by a redundant activity.
Five cell lines were subjected to IL-2 deprivation and microarray analysis, three primary, IL-2-dependent T lymphoblast lines derived from three different donors (denoted "6", "43" and "j") and two immortalised T cell lines (denoted line 5 and S9), examined in three biological replicates each. Altogether, 18 microarray hybridisations were performed.
In the first step, unsupervised Principal Component Analysis was performed to assess the reliability of the experiments. Biological replicates were shown to cluster together. The gene expression profile of all the analysed cell lines was shown to be affected by IL-2 withdrawal, but the major variability of the gene expression profile in the analysed dataset was found between primary and immortalised cell lines (unpublished data). Differences between primary and immortalised cells were much greater than those resulting from IL-2 withdrawal. Therefore, our analysis was carried out in a block design, and changes in gene expression following IL-2 withdrawal were assessed separately in primary and immortalised cell lines. By means of the parametric non-corrected t-test we selected 158 differentially expressed genes at p < 0.001 (see Additional File 1: Table A1 for genes differentiating between samples before and after IL-2 withdrawal, selected by a non-paired analysis). The probability of a random selection of such a number of genes was low (0.03 in a global test) which confirmed the reliability of the observed differences. A false discovery rate for the selected genes ranged from 0.003 to 0.12. The hierarchical clustering of the 158 differentially expressed probe-sets obtained in this analysis resulted in an ideal discrimination between the samples before and after IL-2 withdrawal, with a rather uniform pattern of expression (Figure (Figure11).
Having confirmed the statistical significance of differences related to IL-2 withdrawal, we have selected the most consistently changed genes in all five different cell lines analysed, by employing a paired analysis (of the corresponding cells before and after IL-2 withdrawal). The univariate non-corrected permutation test revealed 166 probe-sets changed significantly at P < 0.001 (while in a global test, a borderline significance of P = 0.0625 was obtained). A hundred and eight of these genes were repressed after growth factor withdrawal, and the remaining 58 were up-regulated (see Additional File 2: Table A2 for genes differentiating between samples before and after IL-2 withdrawal, selected by a paired sample analysis). Out of those, DUSP6, OSM, CISH, SOCS2, LIF and IL13 were found to be the most prominently down-regulated (more than 20-fold), and the expression of EMB, SOX4, RPS24, WNT5B, SLC1A4, SQSTM1, SUMF1, TMEM1, STK38, LAMP2 and C17orf42 was found to be the most increased (2-fold to 3.5-fold change).
We carried out an external validation of the results obtained in the microarray analysis by a qRT-PCR. For validation we used the same cell lines and we chose 8 genes from the list of 158 most significantly changed transcripts, and a further 4 genes, known for their interesting function, with lower significance of differences. Most of the qRT-PCR results confirmed the microarray findings (Table (Table1):1): 7 of 8 top genes and 3 of 4 less significant transcripts were confirmed in a qRT-PCR experiment, usually with a similar fold-change.
To select transcripts consistently changed upon IL-2 withdrawal, we combined the two lists of 166 and 158 significant genes, obtained in the paired and non-paired analyses, respectively. We identified a set of 53 genes comprised of 13 up-regulated and 40 down-regulated genes, which we designated a "T lymphocyte cell cycle exit signature" (Table (Table22).
The up-regulated genes identified after IL-2 deprivation in both primary and immortalised cells included SQSTM1, ECOP, YY1AP1, RPS24, TMEM1, LRRC8D, C1orf63, ASAH1, SLC25A46 and MIA3, some of which have been associated with cellular growth control.
The SQSTM1 protein, previously known as p62 protein, interacts with the ubiquitinated proteins to mediate their clearance and is an important scaffold molecule in the RANK-NF-kappaB signalling pathway . In addition, while the p62 protein is also a ligand for Lck kinase of the Src protein-tyrosine kinase family that is involved in T-cell receptor-dependent activation , Köller et al.  stress the role of Lck in transient T cell unresponsiveness, mediated by a selective IL-2 deficiency. Thus, SQSTM1 may play an important role in cell growth, especially because it may facilitate cell survival through signalling cascades for example those that result in Akt activation, by interaction with protein kinase C, as shown in human neuronal cells .
Interestingly, the product of the ECOP gene, which was found here to be up-regulated, has also been previously shown to be involved in NF-kappa B-related regulation of cellular growth. ECOP (EGFR-Co-amplified and over-expressed protein) has been implicated as a key regulator in the NF-kappaB signalling, and it has been postulated that high-level, amplification-mediated ECOP expression, such as that occurring in tumours with amplified EGFR, could contribute to resistance to apoptosis .
The identified genes may as well be associated with other cell proliferation pathways. For example, a product of the YYAP1 gene, called YY1-associated protein, a co-activator of the YY1 (Yin Yang 1) transcription factor (YYAP1), had been initially identified as a long splice variant of HCC (hepatocellular carcinoma)-specific protein encoded by the HCCA2 gene (HCC-associated gene 2). An over-expression of HCCA2 has been found to lead to a cell cycle arrest at G0/G1 phase and to an inhibition of cell proliferation .
For the other up-regulated genes distinguished in this study, we could not identify links to the known mechanisms of cell cycle arrest, but some of the genes have been previously associated with different cancers. An over-expression of ASAH1 has been found in different malignant tumours, including prostate cancer, head and neck squamous cell cancers and T-cell large granular lymphocyte (LGL) leukaemia [26-28]. An inhibition of acid ceramidase, the product of ASAH1, the over-expression of which in cancer cells has been implicated in drug resistance, is suggested as an efficient and promising novel treatment strategy [29,30]. The ribosomal protein gene RPS24 has been associated with hepatocellular carcinoma . The MIA3 (TANGO) protein has been identified as a tumour suppressor in malignant melanoma and in colon and hepatocellular carcinomas [31,32].
Out of genes encoding for the suggested crucial regulators of T lymphocyte quiescence, such as TOB and members of the KLF and FOXO transcription factors' families, as well as TSC-22 and Dyrk1, which had previously been found to be overexpressed in resting T cells [1,18,33-35], we only found the FOXO3 gene to be slightly, but significantly up-regulated after IL-2 withdrawal. This discrepancy might be explained by differences between the different cell models employed (mouse vs. human) and by different experimental settings. On the other hand, the activity of these transcription factors may be regulated posttranscriptionally, so that their mRNA levels may not necessarily reflect their activity. Thus, the expression of some of the known FOXO target genes was examined and found to be affected following IL-2 withdrawal. IL7R, the FOXO1 target gene was up to 4-fold overexpressed in the primary T lymphocytes, but not in immortalised T cells. CDKN1B (p27), another FOXO target in human, mouse and C. elegans models [1,34-37], encoding a cell-cycle regulatory protein, presented slightly, but significantly increased expression. Similarly, Pink1, a recently identified FOXO target, the product of which helps to protect lymphocytes from apoptosis after growth factor deprivation , was among genes of significantly elevated expression in our model of cell-cycle exit. Studies of Murphy  on C. elegans model pointed at metallothionein genes as possible FOXO targets. In the present study, several metallothionein genes (MT1H, MT1X, MT1G, MT2A, MT1M, MT1E, MT1F) were overexpressed (2.6 to 3.8 times) following IL-2 withdrawal. In addition, we show here, that Cyclin D expression, reported to be downregulated by FOXO and KLF [1,35], decreased after IL-2 withdrawal. However, some genes reported by other research groups as human FOXO targets, such as p130 gene, Cyclin G2 , or IGF1R, BCL-2BP, or Sestrin 1 found to be FOXO3a-inducible genes in mouse T cells , as well as most of the KLF targets summarised by Tsachanis  and Yusuf , were unaffected by IL-2 deprivation in our experimental model.
A set of genes identified as significantly down-regulated 8 hours following IL-2 withdrawal included PIM1, BCL2, IL-8, HBEGF, DUSP6, OSM, CISH, SOCS2, LIF and IL13. It is notable that some of the identified genes, namely PIM1, BCL2, DUSP6, CISH and SOCS2 were previously demonstrated to be the IL-2 target genes or the regulators of IL-2 signalling [39-41]. In addition, it is clear that PIM1, CISH and suppressors of cytokine signalling (SOCS) are implicated in the biological actions of IL-2 , and that IL-2 was found to regulate mRNA levels of SOCS2, CISH, DUSP5 and DUSP6 . It has also been demonstrated that OSM can quickly up-regulate CISH . The coordinated repression of OSM and CISH that we showed in T cells following IL-2 withdrawal is consistent with these observations.
To perform a functional analysis of the up- and down-regulated genes, the 158 genes selected in the non-paired comparison between IL-2-deprived and non-deprived lymphocytes were first analysed by Gene Ontology (GO) Categories (see Additional File 3. Table A3 for gene ontology analysis of genes selected in a non-paired sample comparison between IL-2-deprived and non-deprived lymphocytes). The number of genes changed in each category was compared with the expected quantity of genes and the greatest enrichment was found in two molecular function (MF) classes, GO:0016868, intramolecular transferase activity, phosphotransferases; (observed/expected ratio = 43.9) and GO:0005126, hematopoietin/interferon-class (D200-domain) cytokine receptor binding, (observed/expected ratio = 26.3) and two biological process (BP) classes, GO:0007259, JAK-STAT cascade (observed/expected ratio = 16.7) and GO:0009968, negative regulation of signal transduction, (observed/expected ratio = 6.8). Secondly, a gene set analysis was performed by comparing the overall significance of gene groups defined by GO categories. Out of the 2348 GO classes 317 were found to be significantly affected by IL-2 withdrawal, as shown by at least one of the three tests used to assess the significance of the differences. The selected GO categories are shown in Table Table3.3. The LS/KS permutation tests, which find gene sets that have more differentially expressed genes among the classes than expected by chance, identified 206 significant gene-sets. Efron-Tibshirani's test, which uses 'maxmean' statistics to identify differentially expressed gene-sets, found 231 significant gene-sets (under 200 permutations). Eleven out of the 215 cellular component categories, 69 of 545 molecular function (MF) categories and 237 of 1588 biological process (BP) categories were significant. Thirty-nine GO categories with p < 0.005 in all of the three applied gene group significance tests were identified (see Additional File 4: Table A4 for overall significance of expression changes following IL-2 deprivation in gene groups defined by GO categories).
The same analytical procedure applied to the database of biological pathways (Biocarta, human database) resulted in 54 significantly up-regulated gene sets. Selected pathways affected by IL-2 withdrawal are shown in Table Table4.4. Both analyses, GO and Biocarta, confirmed the consistent up-regulation of the genes belonging to the JAK-STAT pathway.
Our data demonstrate for the first time that IL-2 withdrawal induces a coordinate repression of the same set of genes that have been found to be induced during T cell activation as IL-2 targets [44,45]. A similar "symmetry" has been previously shown in human fibroblasts, where approximately half of the genes of the early response to serum stimulation were correspondingly repressed at the cell cycle exit after mitogen deprivation .
A few previous studies have shown that some mRNAs that are expressed in resting T and B lymphocytes become repressed following cell activation, and it was suggested that cell quiescence is under an active transcriptional control . Similarly, as shown recently, exit from the cell cycle of human fibroblasts is under control of a "quiescence program" dependent upon a set of genes that enforce the non-dividing state, and ensure the reversibility of the cell cycle arrest [15,17]. Our data support the idea of cell cycle arrest as an active state, controlled by some up-regulated genes.
While in vitro T lymphocytes deprived of a growth factor would eventually undergo apoptosis, in vivo a small proportion of antigen-activated cells would exit to G0 and become quiescent. Cellular quiescence is thought to be an indispensable state for the maintenance of lymphocyte homeostasis following immune response, and therefore it is an important barrier against tumorigenesis [46-48]. Some autoimmune and chronic inflammatory disorders involving an excessive lymphocyte proliferation, were found to be associated with increased risks of lymphoma [44,45]. Thus, it is very likely that a deregulated expression of genes responsible for cellular quiescence may contribute to the development of some lymphoid malignancies.
To assess whether T cells adopt quiescent, memory cell phenotype upon IL-2 withdrawal in our model, we compared the gene expression patterns obtained with those described by other researchers as characterising memory vs. naïve vs. effector T cells. Two molecules known to facilitate T cell homing to lymphoid tissues, SELL (CD62L) and CCR7 are expressed predominantly by naïve and memory T cells [49-52], but CCR7 expression seems not to be affected by growth factor withdrawal . In our model of cell cycle exit the SELL gene was up-regulated in primary T cell lines, while CCR7, was unchanged in primary T cells and slightly down-regulated in immortalised T cells after IL-2 deprivation. Out of other genes identified by Holmes et al.  as up-regulated in memory T cells, T cells deprived of IL-2 overexpressed SMAD3, TRIM66 and NAP1L3. There was also an increased expression of IL7R, thought to be expressed in memory  and naïve [50,52] T cells, but this was found not significant by the qRT-PCR (Table (Table1)1) when analysed for all T cells, possibly because this was truth only in primary T cells. The CD58 (LFA3) gene, encoding an adhesion molecule known to be expressed on activated and memory T cells , was indeed highly expressed in actively growing cells and consistently up-regulated after IL-2 withdrawal in our experimental system.
Holmes et al.  reported Granzyme B and KLRD1 expression to be induced in activated T cells. Our experiments showed that growth factor-deprived T cells presented decreased Granzyme B expression, while KLRD1 expression was slightly upregulated, but only in primary T cells. A hierarchical clustering of some genes outlined by Holmes et al.  is shown in an Additional file 5: Figure A1.
Haining et al.  demonstrated a gene expression profile common for CD8 and CD4 memory T cells. We showed that some of the genes of this profile were significantly changed at cell cycle exit; there was an increased expression of S100A4, S100A11, ANXA1, ANXA2, CRIP1 and CASP1, which, according to Haining et al. , are up-regulated in memory T cells, and a reduced expression of Granzyme A, downregulated in memory cells .
In summary, T cells withdrawn from cell cycle gained some crucial characteristics of memory T cells. Nevertheless, the cell cycle exit signature that we describe here is different from memory T cell gene expression signatures published so far.
Further studies will reveal whether there is a common mechanism of quiescence for different cell types, such as the "common quiescence program" in fibroblasts stimulated to exit to G0 by different stimuli . Interestingly, it has been found recently that a group of genes required for cell cycle exit in human fibroblasts following serum deprivation, were coordinately repressed in many types of human cancers, and repression of these genes predicted an increased risk of cancer progression and death in breast cancer patients . Our studies on a model of human primary and immortalised IL-2-dependent T lymphoblasts extended these observations and brought new insights into the understanding of molecular events important in cell cycle exit, deregulation of which may relate to cancer development.
Cell cycle exit of growth factor-deprived primary and immortalised T lymphocytes is characterised by a gene expression signature, comprising 13 up-regulated and 40 down-regulated genes. Some of these genes are implicated in the mechanisms of cell growth control and show deregulated expression in different tumours. The gene set identified in the non-paired analysis was enriched in transcription, cell cycle, cell growth, proliferation and differentiation, cell adhesion, and immune function genes. Cells at cell cycle exit seem to be under an active transcriptional control. Upon IL-2 withdrawal, some genes become repressed, including those known to be induced during T-cell activation, while others become up-regulated. The identification of genes involved in cell cycle exit and quiescence provides new hints for further studies on the molecular mechanisms regulating the non-dividing state of a cell, the mechanisms closely related to cancer development and to many biological processes.
Primary, IL-2-dependent T lymphoblast cell lines ("6", "43" and "j") were generated and propagated as previously described . Briefly, peripheral blood mononuclear cells (PBMC) derived from three healthy donors ("6", "43" and "j") were suspended in a standard medium (RPMI 1640, 10–12% FCS, 50 μg/ml gentamycin [Sigma]), activated for 24 h with 20 μg/ml of wheat germ agglutinin (Pharmacia) and subsequently cultured in the standard medium supplemented with 20 U/ml of recombinant (r) IL-2 (R&D). Flow cytometry cell surface analysis showed that the three primary T lymphoblast cell populations were composed of T cells exclusively, as they expressed all the pan-T markers, and a subset of the cells also expressed CD4 or CD8 surface molecules (unpublished data).
Viable cells were assessed every 3–7 days by the trypan blue exclusion test and resuspended to the concentration of 0.5–2.0 × 105 cells/ml in a fresh standard medium supplemented with rIL-2. The number of population doublings (PD) was calculated according to the formula:
The cells "6", "43" and "j" were collected after 4–5, 7–8 and 6–7 PD, respectively. The samples of the three primary cell lines were analysed as biological replicates.
Spontaneously immortalised IL-2-dependent T cell lines were derived from normal spleen (line 5) and from PBMC derived from a Nijmegen Breakage Syndrome patient (line S9), as previously described [56,57]. Samples of the immortalised cell lines were collected in three biological replicates each.
The differential gene expression analysis was performed in cell samples collected directly from the IL-2-containing cultures (samples: 6/0, 43/0, j/0, line 5/0 and S9/0) and after 8-hrs incubation of the lymphoblasts in a standard medium without IL-2, following triple thorough washing of the cells in a standard medium (samples: 6/8, 43/8, j/8, line 5/8 and S9/8). In the cultures deprived of IL-2 for 8 hours, no apoptosis was observed in the primary T lymphoblast lines, and in the immortalised cell lines a negligible proportion of cells had the characteristics of early apoptotic cells, as detected by flow cytometry following Annexin-V staining with the use of the Annexin-V-FLUOS Staining Kit of ROCHE (Table (Table5).5). But S-phase depletion and G1/S block were observed, as evaluated by flow cytometry with the use of the BD Cycletest Plus (Becton Dickinson) (unpublished data).
Total RNA was isolated using Nucleospin RNA II kit (Macherey-Nagel, Germany). RNA quality was assessed using the Agilent 2100 Bioanalyzer and RNA 6000 Nano Chip Kit (Agilent Technologies).
The Affymetrix GeneChip Human Genome U133 Plus 2.0 arrays, which measure the expression level of over 47,000 transcripts and variants, including 38,500 well-characterised human genes, were used. cDNA synthesis was carried out from 5 μg of RNA (One-Cycle cDNA Synthesis Kit, Affymetrix). After purification (GeneChip Sample Cleanup Module), 7 μl of double-stranded cDNA were used for biotinylated cRNA synthesis (IVT Labeling Kit, Affymetrix). Labeled cRNA was purified using a GeneChip Sample Cleanup Module, fragmented and hybridised with genome array. Washing, staining with streptavidin-phycoerythrin conjugate and scanning of the arrays in Affymetrx GeneChip 3000 scanner were performed as recommended by the Affymetrix Gene Expression Analysis Technical Manual.
All 18 arrays were normalised by the GC-RMA algorithm using a Bioconductor gcrma package. We filtered out all genes with small amplitude of differences (less than 10% of samples with at least 1.5-fold difference in either direction from the median gene expression among all samples). The filtered dataset comprised of 18770 probe-sets. The full dataset has been deposited in the Gene Expression Omnibus repository (accession no GSE13909).
Differentially expressed genes were identified by a random-variance t-test, with the statistical significance threshold set at p < 0.001. A stringent significance cut-off was used to limit the number of false positive genes obtained (at this threshold no more than 19 false positives are expected in each analysis). A global test of the differences in expression profiles between the classes was also performed by sample permutation. Simultaneously, we assessed a False Discovery Rate for each probeset using the Benjamini-Hochberg method.
To examine the effects of IL-2 withdrawal, samples were compared either according to the type of the cell line (primary/immortalised) in a block design or in a paired sample t-test analysis (for this analysis the replicates of the same sample were averaged).
For cluster analysis, genes identified in the non-paired analysis were used to obtain a heatmap plot with dendrograms, using Ward's hierarchical cluster algorithm (Bioconductor).
Gene Ontology classes for the resulting lists of genes were compared to the number of genes expected, considering the number of genes on the microarray in each category. GO classes and parent classes with at least 5 observations in the selected subset and with an observed vs expected ratio of at least 2 were analysed.
The statistical significance of global gene expression changes in gene ontology categories and BioCarta pathways was assessed by LS/KS permutation tests and the Efron-Tibshirani's GSA maxmean test. All data analyses were performed using the BRB-ArrayTools (developed by Dr. Richard Simon and BRB Array Tools development team) or by Bioconductor packages.
Based on the microarray results, the following genes were selected for validation by the real-time RT-PCR technique: DUSP6, ASAH1, HBEGF, MT1F, MT1X, IL7R, IL8, CLK4, C3orf59, SOCS2, OSM and CISH. Specific primers were designed using Primer Express v.2.0 software (Applied Biosystems). Detailed primer sequences are available on request.
The cDNA was synthesised with the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems) using 1 μg of total RNA as a template, and Oligo-dT-Primers (Applied Biosystems), according to the manufacturer's instructions. The real-time (= quantitative, q) RT-PCR reactions were performed with Power SYBR Green PCR Master Mix (Applied Biosystems) using the 7000 Real-Time PCR system (Applied Biosystems). Reaction mixtures (25 μl) were subjected to a two-step program: 95°C for 10 minutes followed by 40 cycles of 95°C for 15 s and 60°C for 1 minute. The data were analysed using ABI PRISM 7000 SDS Software. The relative fold difference in gene expression was calculated by the ΔΔCt method using the ACTB gene as a normalisation control. The statistical significance of the fold change was assessed by the ANOVA test (Statistica package).
BP: biological process; GO: gene ontology; IL: interleukin; MF: molecular function; PBMC: peripheral blood mononuclear cells; r: recombinant.
MCh conceived of the study and participated in its design and coordination, carried out RNA isolation, analysed data and interpreted the study, drafted, wrote and edited the manuscript. JKS conceived of the study and participated in its design and coordination, carried out cell cultures, analysed data and interpreted the study, drafted and wrote the manuscript. MG carried out real-time PCR verification of gene expression, participated in drafting the manuscript. MO-W performed microarray gene expression profiling. MJ performed the biostatistical analyses, participated in drafting the manuscript. AP participated in the biostatistical analyses. BJ participated in the design of the study. JS conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.
Table A1. Genes differentiating between samples before and after IL-2 withdrawal identified by a non-paired sample T-test. Genes selected at the nominal univariate comparison level of 0.001.
Table A2. Genes differentiating between samples before and after IL-2 withdrawal identified by a paired sample T-test. Genes selected at the nominal univariate comparison level of 0.001.
Table A3. Gene Ontology (GO) analysis of genes selected in a non-paired sample comparison between IL-2-deprived and non-deprived lymphocytes. The number of genes changed in each category was compared with the number of expected occurrences. Only GO classes and parent classes with at least 5 observations in the selected subset and with an "observed vs. expected" ratio of at least 2 were shown.
Table A4. Assessment of overall significance of expression changes following IL-2 deprivation in gene groups defined by Gene Ontology categories. Three independent tests: LS, KS and Efron-Tibshirani GSA were applied to select significantly affected gene classes.
Figure A1. Hierarchical clustering analysis of selected memory T cell genes. The analysis based on the expression values of some of the genes identified by Holmes et. al  to characterise memory T cells. Samples before IL-2 withdrawal are in the light-blue columns, samples after IL-2 withdrawal in the navy-blue columns. Primary samples are in the yellow column, immortalised samples are in the green one.
We thank Prof. Leendert Looijenga, PhD, for the critical reading of the manuscript. The study was supported by the grant no 2 P05A 150 29 (to JKS, MCh, MG, and JS) and in part by the grant 2 P05A 022 30 (to BJ) of the Polish Ministry of Science and Higher Education.