- Open Access
T cell epigenetic remodeling and accelerated epigenetic aging are linked to long-term immune alterations in childhood cancer survivors
Clinical Epigenetics volume 10, Article number: 138 (2018)
Cancer treatments have substantially improved childhood cancer survival but are accompanied by long-term complications, notably chronic inflammatory diseases. We hypothesize that cancer treatments could lead to long-term epigenetic changes in immune cells, resulting in increased prevalence of inflammatory diseases in cancer survivors.
To test this hypothesis, we established the epigenetic and transcriptomic profiles of immune cells from 44 childhood cancer survivors (CCS, > 16 years old) on full remission (> 5 years) who had received chemotherapy alone or in combination with total body irradiation (TBI) and hematopoietic stem cell transplant (HSCT). We found that more than 10 years post-treatment, CCS treated with TBI/HSCT showed an altered DNA methylation signature in T cell, particularly at genes controlling immune and inflammatory processes and oxidative stress. DNA methylation remodeling in T cell was partially associated with chronic expression changes of nearby genes, increased frequency of type 1 cytokine-producing T cell, elevated systemic levels of these cytokines, and over-activation of related signaling pathways. Survivors exposed to TBI/HSCT were further characterized by an Epigenetic-Aging-Signature of T cell consistent with accelerated epigenetic aging. To investigate the potential contribution of irradiation to these changes, we established two cell culture models. We identified that radiation partially recapitulated the immune changes observed in survivors through a bystander effect that could be mediated by circulating factors.
Cancer treatments, in particular TBI/HSCT, are associated with long-term immune disturbances. We propose that epigenetic remodeling of immune cells following cancer therapy augments inflammatory- and age-related diseases, including metabolic complications, in childhood cancer survivors.
Improvement in the cure rate of childhood cancer, and in particular acute lymphoblastic leukemia, has led to a significant increase in the number of long-term survivors . However, more than 95% of childhood cancer survivors (CCS) develop chronic health conditions by the age of 45 . The most frequently reported complications include cardiac, endocrine, or neurological conditions , chronic diseases that all share an inflammatory component. While a combination of cancer treatment, lifestyle and environmental factors could contribute to the development of these chronic conditions [4, 5], we and others have identified total body irradiation (TBI) and hematopoietic stem cells transplant (HSCT) as major risk factors for the development of metabolic, cardiovascular, and other health conditions in childhood cancer survivors [6, 7]. When compared to chemotherapy alone, exposure to radiation resulted in a higher cumulative incidence of several chronic health conditions 20 years post-diagnosis, with TBI being the mode of radiation causing the highest incidence of severe conditions . We further reported in a cohort of 248 CCS that the prevalence of metabolic complications was more than doubled at a median of 12.7 years after remission . More recently, we have established in a mouse model that irradiation in the absence of HSCT led to altered muscle and adipose lineage commitment and metabolic dysfunction later in life .
Immune disturbances contribute to cardiovascular, neurological, and metabolic conditions, notably in early development of these diseases [9, 10]. Radiation-induced immune response, including macrophage activation, neutrophil, and lymphocyte recruitment, leads to the production of pro-inflammatory mediators to support anti-tumor activity and clearance of apoptotic cells [11, 12], with a similar response having been described in response to chemotherapy . Although this inflammatory process is usually confined to the acute phase of the radiation response , long-term immune changes could also persist following irradiation, chemotherapy, or HSCT [15, 16]. As there is still limited evidence to support the existence of long-term chronic inflammation, the mechanisms responsible for long-term chronic immune disturbances in cancer survivors and their potential consequences on survivors’ health still remain unknown.
We hypothesized that childhood cancer survivors, in particular when exposed to TBI/HSCT, present long-lasting immune alterations that could contribute to diseases later in life. Our results show that exposure to TBI/HSCT is associated with persistent chronic inflammation, a remodeling of immune cells epigenome and altered immune functions. We further establish that T cell from childhood cancer survivors presents an epigenetic signature consistent with accelerated aging. We propose that persistent alterations of the immune system in childhood cancer survivors contributes to higher risk of secondary diseases, including cardio-metabolic diseases.
To study the long-term effect of cancer treatment on immune functions, we recruited childhood cancer survivors who had received chemotherapy alone or in combination with TBI/HSCT in the course of their treatment, with clinical characteristics and cancer history presented in Table 1. The majority of our population was composed of survivors of acute lymphoblastic or myeloid leukemia, as well as Hodgkin’s and non-Hodgkin’s lymphoma, and all survivors received chemotherapy. Our two groups were age-matched, with no difference in age at inclusion, at diagnosis, or in number of years on remission (Table 1). Survivors in the TBI/HSCT group tended to be shorter and lighter (p = 0.05 and p = 0.07 respectively).
TBI/HSCT is associated with long-term epigenetic remodeling in T cell
Reports have suggested that inflammation could persist years after cancer treatments , potentially through a long-term “memory” of the treatment, implying epigenetic mechanisms. Thus, we assessed the DNA methylation profile of immune cells. CD4+ (p = 0.05) and CD8+ T cell (p = 0.07) tended to show lower global DNA methylation levels in CCS treated with TBI/HSCT, while monocytes tended to show higher methylation levels (p = 0.06), with no difference in other populations (Fig. 1a). We then tested if these changes in methylation were gene-specific and performed reduced representation bisulfite sequencing (RRBS) in CD4+ T cell. We found differential methylation at cytosines which corresponded to 419 different neighboring genes in the TBI/HSCT group (Fig. 1b, Additional file 1: Table S1). Out of the 419 differentially methylated regions (DMRs), 9 were identified in enhancers and 47 in promoters (Additional file 1: Table S1). Differential methylation was also found at microRNA (miRNA) loci (Additional file 1: Table S1), including miR-942 and miR-378c, which have been implicated in acute response to irradiation, or miR-124 which is linked to T cell polarization and inflammation [18, 19]. Using GO analysis on our 419 coding genes containing differentially methylated cytosines, we identified enrichments for GO terms linked to oxidative stress (15/77 GO terms found in our data, Fig. 1c, Additional file 1: Table S2) and genes related to immune processes (21/77 terms) including positive and negative regulation of cytokines production (Fig. 1c, Additional file 1: Table S2). Transcriptomic analysis of CD4+ T cell by RNA-seq returned a small subset of genes differentially expressed (29, Fig. 1d, Additional file 1: Table S3a). Although we did not detect differential methylation at the proximity of the differentially expressed genes, except for adhesion G protein-coupled receptor G1 (also known as GPR56), we found enrichment for GO terms related to inflammation, similar to the DMR-associated genes (Fig. 1e, Additional file 1: Table S3B). Collectively, these results support that TBI/HSCT stably alters the epigenome of T cell.
Cancer survivors treated with TBI/HSCT show long-term altered immune function
To determine the consequences of the remodeling of the T cell epigenome, we analyzed T cell function (Fig. 2a). In mitogen-stimulated CD4+ and CD8+ T cell from the TBI/HSCT group, we found a greater frequency of interferon (IFN)-γ but not interleukin (IL)-4-producing cells (Fig. 2b, Additional file 1: Figure S1a). Intracellular signaling pathways, including mTOR (mammalian target of rapamycin) and MAPK (mitogen-activated protein kinases), are critical regulators of T cell polarization and cytokines production (Fig. 2a). Quantification of p38, ribosomal protein S6 kinase1 (S6 k1), c-Jun N-terminal kinase (JNK) and Akt phosphorylation in CD4+ and CD8+ T cell showed a marked increase for all four proteins in both groups upon mitogen stimulation (Fig. 2c, Additional file 1: Figure S1b). However, the TBI/HSCT group showed over-activation of p38 and S6 k1 under resting conditions in CD4+ (Fig. 2c) and CD8+ cells (Additional file 1: Figure S1b), potentially linked to the increased frequency of type 1 cytokine-producing T cell. These functional changes in T cell from the TBI/HSCT group were associated with elevated systemic levels of IFN-γ and tumor necrosis factor (TNF)-α (p < 0.05, Fig. 2d), suggesting a link between radiation, HSCT, and long-term inflammation.
Acute exposure to radiation triggers a transient stress response resulting in pro-inflammatory cytokines production . Thus, we tested the contribution of direct irradiation to T cell polarization and MAPK/mTOR activation using an in vitro model (Fig. 3a). Direct irradiation did not affect p-p38 or pS6 k1 (p > 0.05, Fig. 3b), or pJNK, pAkt-Ser473, and the frequency of IL-4/IFN-γ-producing cells in Jurkat cells (Additional file 1: Figure S2a-b). We then investigated the indirect effect of irradiation by exposing mitogen-stimulated PBMCs collected from healthy donors to the conditioned media obtained from irradiated fibroblasts (Fig. 3c). We found that conditioned media from fibroblasts irradiated with a dose of 3Gy increased phosphorylation of p38 and S6 K1 in CD4+ cells (Fig. 3d) but had no effect on pJNK or pAkt-Ser473 (Additional file 1: Figure S2c). Exposure to conditioned media from fibroblasts irradiated with a dose of 3Gy also had no effect on the frequency of CD4+ cells or CD8+ cells producing IFN-γ (p = 0.14, Fig. 3d) or IL-4 (p > 0.05, Fig. 3d, Additional file 1: Figure S3a-b for CD8+ cells). These results suggest that irradiation may exert its effect on T cell through secreted factors. To identify proteins secreted by irradiated cells, we analyzed cell culture supernatants from irradiated (1) fibroblasts or (2) fully differentiated adipocytes at 24 h and 7 days post-exposure using mass spectrometry (Fig. 4a). Most detected proteins were secreted in both irradiated or non-irradiated conditions (Fig. 4b), with some proteins uniquely detected in the irradiated or non-irradiated conditions, although this did not reach statistical significance (Additional file 1: Table S4). We next analyzed plasma samples from CCS and detected more than 50 proteins with a distinct expression pattern (Fig. 4c). GO analysis identified a specific enrichment for pathways related to cell stress, immune processes and inflammation (Fig. 4d). Of interest, the alpha-1-antitrypsin protein (encoded by the SERPINA1 gene) was decreased in the TBI/HSCT group. Given the role of alpha-1 antitrypsin in the inhibition of pro-inflammatory cytokines secretion , decreased levels of this protein could participate in the immune changes in CCS treated with TBI/HSCT. Our cluster analysis suggested that the expression pattern of these proteins was group-specific. Collectively, our data suggest that secreted factors participate in a long-term alteration of the immune system in survivors exposed to irradiation and hematopoietic stem cell transplant.
Childhood cancer survivors are characterized by accelerated epigenetic aging
Given that T cell epigenetic remodeling specifically affected genes involved in inflammation and oxidative stress, two processes associated with aging, we investigated the contribution of epigenetic factors in cancer treatment-induced accelerated aging. Using the Epigenetic-Aging-Signature method , we found that both groups of survivors had increased predicted biological age (Fig. 1f), with CCS treated with TBI/HSCT showing a higher predicted biological age (p = 0.01, Fig. 1f). Moreover, the chronological age was correlated to the predicted age in non-irradiated CCS (n = 8, r = 0.85, p < 0.01, Fig. 1f) but not in the TBI/HSCT group (n = 7, p > 0.05, Fig. 1f). Since all survivors treated with TBI also received HSCT, our results could indicate that the predicted biological age reflect the age of the donor as previously suggested . However, we did not find any difference between the donors’ and the recipients’ chronological age (Fig. 1f). In addition, the inclusion of the donors’ age in the model did not change the correlation between chronological and predicted age (Fig. 1f), suggesting that TBI/HSCT could be responsible for the accelerated epigenetic aging in T cell. When we compared the methylation levels on three specific CpG sites used to establish the Epigenetic-Aging-Signature, we found a hypermethylation only on PDE4C in the TBI/HSCT group (p > 0.05 for ASPA and ITGA, p < 0.01 for PDE4C, Additional file 1: Figure S1c). Consistent with the Epigenetic-Aging-Signature, transcriptomic analysis of T cells showed that more than 10% of differentially expressed genes belonged to the “innate and adaptive immunity, cytokine, and chemokine” cluster of genes, which constitutes genes previously reported to contribute to the aging signature of peripheral blood (Additional file 1: Table S3A) . The present results support a role for cancer treatment in the accelerated aging of T cell in cancer survivors.
Multiple reports have suggested an increased prevalence of inflammatory diseases in cancer survivors. We thus investigated the presence of cardio-metabolic risk factors in our cohort. Consistent with previous studies [3, 7], CCS treated with TBI/HSCT had impaired metabolic profile, including increased triglycerides and reduced high-density lipoprotein cholesterol (HDL-C) levels (p < 0.05, Table 1). Importantly, these changes were not due to a difference in age, age at diagnosis, body mass index (BMI), or waist-to-hip ratio (p > 0.05, Table 1). The level of insulin resistance was also higher in the TBI/HSCT group (p = 0.04, Table 1), while fasting blood glucose and insulin levels also tended to be higher in that group (p = 0.1 and p = 0.07 respectively; Table 1). These results confirm that CCS, in particular those who received TBI/HSCT, are at increased cardio-metabolic risks with a potential link with chronic inflammation.
Our results suggest that cancer treatment stably remodels the T cell epigenome at specific genes controlling inflammation and oxidative stress in CCS. Using various markers of immune cell function and signaling pathways regulating pro-inflammatory cytokines production, we establish a link between T cell epigenetic remodeling and low-grade inflammation in survivors, in particular those who received TBI/HSCT. We identified that TBI/HSCT is associated with accelerated aging and propose a role for epigenetic changes in the elevated risk of age- and inflammation-related diseases long after treatment in CCS .
We found that CCS treated with TBI/HSCT are characterized by systemic inflammation. The secretion of cytokines by immune cells is regulated by epigenetic mechanisms . Our global DNA methylation analysis of T cell from survivors supports a dramatic remodeling of the epigenome after TBI/HSCT. Radiation-induced long-term epigenetic changes have been reported in cultured human keratinocytes and rodents [25, 26] and could mediate the imprinting of a long-term “memory” of the radiation insult, potentially affecting immune cell function. HSCT has also been associated with an inflammatory response, although limited information regarding its long-term effect has been reported . Nonetheless, it can be hypothesized that even sub-clinical levels of graft-versus-host disease could have contributed to the immune changes reported here. We have previously shown that the epigenome of lymphocytes is linked to environmental influences, notably nutritional [27, 28]. Epigenetic mechanisms including changes in DNA methylation, histone modifications, or even miRNAs are well-described as modulators of immune cell differentiation and function [29, 30]. For example, differentiation and polarized activation of T helper or cytotoxic cells are controlled by global DNA methylation at particular gene loci [31, 32]. Although we did not observe specific changes at previously reported gene loci, we identified several differentially methylated genes involved in critical pathways controlling T cell differentiation or cytokine production. It cannot be ruled out that by investigating epigenetic changes in global CD4+ T cell, we might have missed DNA methylation changes specific to effector, memory, or regulatory T cell [31, 32]. In addition, these T cell subtypes present different epigenetic signature [31, 32]. Thus, our results could also suggest changes in the distribution of T cell subsets.
We found a greater frequency of type 1 cytokine-producing T cell in irradiated survivors, with these cells showing higher activation of p38 MAP kinase pathway and mTOR Complex 1 (mTORC1) signaling (S6 k1), potentially through a bystander effect. Direct and indirect exposure to ionizing radiation activate MAPKs (p38/JNK) and the PI3K (phosphoinositide 3-kinase)/Akt/mTOR pathways [12, 33], and both p38 and mTORC1 are critical for Th1 polarization and type 1 cytokines production [34, 35]. The over-activation of both p38 and mTORC1 is consistent with the greater frequency of Th1 cells and the higher levels of pro-inflammatory cytokines in survivors who received TBI/HSCT. The response to radiation involves the mobilization of immune cells and pro-inflammatory cytokine production, although these immune changes usually resolve within hours to days upon exposure . Similarly, HSCT can be associated with an inflammatory response, though the long-term consequences remain unclear . Delayed or chronic immune activation after radiation has been reported in atomic bomb survivors more than 50 years after exposure . Impaired immune function has also been reported up to 10 years following radiation in long-term cancer survivors , although in both cases, the molecular mechanisms were not identified. Our results suggest that such immune changes could result from an indirect effect of irradiation, potentially due to circulating factors. This bystander effect has been previously suggested as playing a key role in the adaptations following irradiation . Although we could not identify specific factors driving the epigenetic remodeling and altered immune function, it is likely that such factors might contribute to these phenomena . We can speculate that the observed changes in the T cell methylome and T cell function could be involved in the perpetuation of the cancer treatment-induced pro-inflammatory programming. Limitations associated with our in vitro model to investigate the direct effect of irradiation, in particular the absence of the regulator of the Akt/mTOR pathway Phosphatase and TENsin homolog (PTEN) in Jurkat cells, forces us to caution when interpreting our results using this cell system. Further work is needed using other cell models to fully clarify the influence of irradiation on the epigenome.
Inflammation and oxidative stress are associated with several processes involved in aging . Recently, changes in both the epigenome and the transcriptome of circulating leukocytes have been linked to the biological (rather than chronological) age of an individual [22, 39]. A marked remodeling of the epigenome encompassing global DNA hypomethylation and gene-specific hypermethylation at particular loci has been reported during aging , consistent with the hypomethylation we observed in T helper cells at the global DNA level. At the gene level, particular loci or even CpG sites have been identified as predictive markers of biological age  and allowed us to identify an increased predicted biological age in cancer survivors. Several components of the radiation response can contribute to accelerated aging or inflammatory diseases, including low-grade systemic inflammation, activation of MAPK/mTOR signaling, epigenetic remodeling, or oxidative stress [12, 33, 38, 40, 41]. In the present study, we identified the presence of such disturbances in cancer survivors more than 10 years post-treatment, which potentially explains the accelerated aging reported in this population . We and others have previously shown that both reactive oxygen species and pro-inflammatory cytokines can affect DNA methylation [42, 43]. We can speculate that the initial inflammatory response and oxidative burst in response to treatment might have contributed to the epigenetic remodeling of immune cells, although the timing and nature of the events that might have led to this “long-term memory” are still to be unraveled. These changes in the T cell methylome might be involved in the perpetuation of the pro-inflammatory programming contributing to the early onset of cardio-metabolic complications in cancer survivors.
Improvements in childhood cancer survival are associated with higher prevalence of treatment-associated diseases later in life, notably of age-related conditions including cardiometabolic diseases. We identified low-grade inflammation and altered immune cell function in survivors treated with TBI/HSCT. This was associated with a dramatic remodeling of the T cell epigenetic signature and features of accelerated epigenetic aging. These immune changes could account for the increased risk of inflammatory and age-associated diseases later in life. A greater understanding of the mechanisms by which cancer treatments affect the epigenome of immune cells may help to decrease long-term complications in cancer survivors.
We recruited 44 childhood cancer survivors, diagnosed between 0 and 18 years, in full remission for more than 5 years and 16 years old and older at the time of enrolment. Survivors who had experienced graft-versus-host disease were not included in this study. This project was approved by the Human Research Ethics Committee (UNSW Sydney, HREC-10/017) and performed in accordance with Helsinki Declaration procedures. The nature of the project was explained and participants provided informed consent. Participants had received chemotherapy only (non-IRR) or a combination of chemotherapy with TBI and HSCT (TBI/HSCT). Anthropometric data and medical history were recorded and BMI calculated.
Fasting blood samples were collected and stored at − 80 °C. Serum insulin (Immulite 2000, Diagnostic Products Corporation, CA, USA), glucose levels, HDL-C, low-density lipoprotein cholesterol (LDL-C), and triglycerides were measured (LX20 analyzer, Beckman Coulter, CA, USA), and the insulin resistance index was calculated (Homeostasis Model Assessment, HOMA-IR) . Systemic cytokine levels including IL-2, IL-4, IL-5, IL-10, IL-12, IL-13, IFN-γ, TNF-α, and granulocyte-monocyte colony stimulating factor were measured (Bio-Plex Pro Human Cytokine Th1/Th2 Assay, Bio-Rad Laboratories Inc., Australia, Intra-assay CV: 7–12% and Inter-assay CV: 6–10%). Results were reported only for cytokines that were detectable in more than 30% of the participants.
Peripheral blood mononuclear cells (PBMCs) were isolated from whole blood and cryopreserved . PBMCs were thawed, left to recover overnight, and mitogen-stimulated for 6 h (2 μg/ml phorbol 12-myristate 13-acetate/ionomycin, Sigma-Aldrich, Australia), with the addition of monensin (BioLegend, CA, USA) after 2 h. PBMCs were then fixed (2% paraformaldehyde), stained for dead cells (ethidium monoazide bromide, Sigma-Aldrich), and permeabilized in ice-cold methanol. Cells were then stained using CD4-V500, CD-PE-Cy7, IL-4-PE, IFN-γ-Pacific blue (BD Biosciences, Australia), and phosphorylated pAkt-Ser473-AF647 (BD Biosciences), pS6k1-AF488, pJNK-AF488, and phosphorylated p38-AF647 (Cell Signaling Technology, MA, USA). PBMCs were analyzed on FACS CantoII (BD Biosciences), with data post-compensated and analyzed using FlowJo (Tree Star Inc., OR, USA). Results were expressed as percentage of IL-4+/IFN-γ+ cells or normalized median fluorescence intensity for pAkt, pS6k1, pJNK, and p-p38 in CD4+ and CD8+ cells.
In vitro and co-culture models
Jurkat cells (Clone E6–1, ATCC, VA, USA) were irradiated with doses of 0, 1, 3, or 6Gy (Xrad 320, Precision X-ray Inc., CT, USA) and left to recover for 24 h. At days 1, 14, and 28, post-irradiation cells were mitogen-stimulated and analyzed by flow cytometry as described above. Adipocytes (3 T3-L1, CL-173, ATCC) were irradiated with doses of 0, 1, 3, and 6Gy and left to recover for 24 h prior to the collection of culture media. PBMCs from six young healthy donors were mitogen-stimulated in conditioned media from irradiated adipocytes and analyzed by flow cytometry as described above.
Protein sample preparation and mass spectrometry
Plasma samples (14ul) from childhood cancer survivors were depleted from the 14 most abundant proteins (Sigma Seppro IgY14 Spin columns). Depleted plasma sample and conditioned media from cultured fibroblasts and fibroblasts differentiated into adipocytes (only 0/3Gy and at 24 h/7 days post-irradiation) were processed as previously described with slight modifications (see “Additional file 1”) . A survey scan m/z 350–1750 was acquired in the Orbitrap and lock mass enabled with data-dependent tandem MS analysis performed using a top-speed approach (cycle time of 2 s). MS2 spectra were fragmented by HCD activation mode and the ion trap was selected as the mass analyzer. Peak lists were generated using Mascot Daemon and submitted to Mascot (version 2.5.1, Matrix Science). Label-free quantitation was carried out using MaxQuant (version 18.104.22.168) and Perseus (version 22.214.171.124) and pathway analysis using G-Profiler.
Global DNA methylation
Global DNA methylation was measured in PBMCs, using surface markers (CD3-PE, CD4, or CD21-APC, CD8, or CD14-PerCP, BD Biosciences) and anti-5-methylcytidine antibody (AbD Serotec, Bio-Rad) conjugated to AF488 (Zenon AF488 Mouse IgG1 labeling kit, Life Technologies, Australia) as previously described .
PBMCs were sorted on Influx Cell Sorter (BD Biosciences) using specific surface markers (CD4 and CD19/CD20-BV421, CD8-BV510, CD3-FITC, CD56-PE-CF594, CD45-PerCP-Cy5.5, and CD14-APC, BD Biosciences). Cells were then lysed before DNA/RNA extraction (AllPrep DNA/RNA/miRNA Universal kit on QIAcube, Qiagen, Australia). CD4+ T cell DNA was processed using reduced representation bisulfite sequencing with minor changes . Genomic DNA was fragmented by overnight incubation with Msp1 restriction enzyme and adenylated with dATP and Klenow DNA polymerase (3′ to 5′ exo minus) before ligation to TruSeq adapters (Illumina, CA, USA). Twelve ligated samples were pooled and processed for bisulfite conversion using EZ DNA methylation kit (Zymo Research, CA, USA). Libraries were enriched by PCR, purified with AMPure beads, and sequenced on a HiSeq (50 bp single-end sequencing, Illumina, Danish National High-Throughput DNA Sequencing Centre, Denmark). Preprocessed reads were aligned to hg38 with Bismark , and differential methylation was analyzed with MethylKit . All regions were annotated with the Bioconductor package ChIPseeker . CpGs with more than 25% methylation change and q value < 0.01 were included for further analysis. Gene Ontology (GO) analysis was performed using a hypergeometric test on individual CpGs and scatter plots generated with REVIGO . T cell-specific enhancers were downloaded from SlideBase [51, 52] and lifted to hg38 using the UCSC liftOver tool , before counting overlaps with DMRs.
DNA methylation was measured on sorted CD8+ T cell at three specific CpG sites (ASPA, ITGA2B, and PDE4C) by bisulfite pyrosequencing. Beta values were then used in a multivariate model to estimate the predicted biological age as previously described .
Whole transcriptome sequencing
RNA from CD4+ T cell was sequenced following standard protocol (Truseq Stranded total RNA protocol, Illumina). After depletion of ribosomal RNA and fragmentation, cDNA synthesis was performed on the first strand, followed by second strand synthesis. cDNA was cleaned using AMPure beads, adenylated and then ligated to adapters. Following an additional beads cleanup, DNA fragments were enriched by PCR and subjected to a final beads cleanup. Libraries were sequenced on a HiSeq (100 bp single-end sequencing, Illumina, Danish National High-Throughput DNA Sequencing Centre). Reads were aligned to hg38 using Rsubread aligner, and the number of reads aligning with Gencode (version 23) was counted with featureCounts . Genes with less than 1 RPKM in half the samples were excluded. Differential expression was calculated with edgeR  using glmFit/glmLRT functions and tagwise dispersion, with correction for multiple testing using the Benjamini–Hochberg procedure. GO analysis was performed using GOrilla , and scatter plots were generated with REVIGO .
Data are expressed as mean +/− standard deviation and analyzed using SPSS Statistics (IBM, NY, USA). Normality of the distribution was tested using the Skewness and Kurtosis tests and anthropometric data, biological markers, cytokines levels, T cell polarization, and global DNA methylation or predicted age were compared using Student t test or Mann-Whitney U test. A chi-square test was used to compare the frequencies between the two groups. The effect of irradiation/mitogen stimulation on intracellular signaling was tested using two-way ANOVA for repeated measures or a combination of Friedman’s ANOVA and Mann-Whitney U test. Correlation levels were tested using the Pearson or Spearman tests.
Coleman M, Forman D, Bryant H, Butler J, Rachet B, Maringe C, Nur U, Tracey E, Coory M, Hatcher J, et al. Cancer survival in Australia, Canada, Denmark, Norway, Sweden, and the UK, 1995–2007 (the International Cancer Benchmarking Partnership): an analysis of population-based cancer registry data. Lancet. 2011;377:127–38.
Hudson MM, Ness KK, Gurney JG, Mulrooney DA, Chemaitilly W, Krull KR, Green DM, Armstrong GT, Nottage KA, Jones KE, et al. Clinical ascertainment of health outcomes among adults treated for childhood cancer. JAMA. 2013;309:2371–81.
Diller L, Chow EJ, Gurney JG, Hudson MM, Kadin-Lottick NS, Kawashima TI, Leisenring WM, Meacham LR, Mertens AC, Mulrooney DA, et al. Chronic disease in the Childhood Cancer Survivor Study cohort: a review of published findings. J Clin Oncol. 2009;27:2339–55.
Ehrhardt MJ, Mulrooney DA. Metabolic syndrome in adult survivors of childhood cancer: the intersection of oncology, endocrinology, and cardiology. Lancet Diabetes Endocrinol. 2015;3:494–6.
Smith WA, Li C, Nottage KA, Mulrooney DA, Armstrong GT, Lanctot JQ, Chemaitilly W, Laver JH, Srivastava DK, Robison LL, et al. Lifestyle and metabolic syndrome in adult survivors of childhood cancer: a report from the St. Jude Lifetime Cohort Study. Cancer. 2014;120:2742–50.
Meacham LR, Sklar CA, Li S, Liu Q, Gimpel N, Yasui Y, Whitton JA, Stovall M, Robison LL, Oeffinger KC. Diabetes mellitus in long-term survivors of childhood cancer. Increased risk associated with radiation therapy: a report for the childhood cancer survivor study. Arch Intern Med. 2009;169:1381–8.
Neville KA, Cohn RJ, Steinbeck KS, Johnston K, Walker JL. Hyperinsulinemia, impaired glucose tolerance, and diabetes mellitus in survivors of childhood cancer: prevalence and risk factors. J Clin Endocrinol Metab. 2006;91:4401–7.
Nylander V, Ingerslev LR, Andersen E, Fabre O, Garde C, Rasmussen M, Citirikkaya K, Bæk J, Christensen GL, Aznar M, et al. Ionizing radiation potentiates high-fat diet-induced insulin resistance and reprograms skeletal muscle and adipose progenitor cells. Diabetes. 2016;65:3573–84.
Geng S, Chen K, Yuan R, Peng L, Maitra U, Diao N, Chen C, Zhang Y, Hu Y, Qi C-F, et al. The persistence of low-grade inflammatory monocytes contributes to aggravated atherosclerosis. Nat Commun. 2016;7:13436.
Winer DA, Winer S, Shen L, Wadia PP, Yantha J, Paltser G, Tsui H, Wu P, Davidson MG, Alonso MN, et al. B cells promote insulin resistance through modulation of T cells and production of pathogenic IgG antibodies. Nat Med. 2011;17:610–7.
Coates PJ, Lorimore SA, Wright EG. Damaging and protective cell signalling in the untargeted effects of ionizing radiation. Mutat Res. 2004;568:5–20.
Formenti SC, Demaria S. Systemic effects of local radiotherapy. Lancet Oncol. 2009;10:718–26.
Apetoh L, Ghiringhelli F, Tesniere A, Obeid M, Ortiz C, Criollo A, Mignot G, Maiuri MC, Ullrich E, Saulnier P, et al. Toll-like receptor 4–dependent contribution of the immune system to anticancer chemotherapy and radiotherapy. Nat Med. 2007;13:1050–9.
Stoecklein VM, Osuka A, Ishikawa S, Lederer MR, Wanke-Jellinek L, Lederer JA. Radiation exposure induces inflammasome pathway activation in immune cells. J Immunol. 2015;194:1178–89.
Fuks Z, Strober S, Bobrove A, Sasazuki T, McMichael A, Kaplan H. Long term effects of radiation of T and B lymphocytes in peripheral blood of patients with Hodgkin's disease. J Clin Invest. 1976;58:803.
Neven B, Leroy S, Decaluwe H, Le Deist F, Picard C, Moshous D, Mahlaoui N, Debre M, Casanova JL, Dal Cortivo L, et al. Long-term outcome after hematopoietic stem cell transplantation of a single-center cohort of 90 patients with severe combined immunodeficiency. Blood. 2009;113:4114–24.
Chow EJ, Simmons JH, Roth CL, Baker KS, Hoffmeister PA, Sanders JE, Friedman DL. Increased cardiometabolic traits in pediatric survivors of acute lymphoblastic leukemia treated with total body irradiation. Biol Blood Marrow Transplant. 2010;16:1674–81.
Chaudhry MA, Omaruddin RA, Brumbaugh CD, Tariq MA, Pourmand N. Identification of radiation-induced microRNA transcriptome by next-generation massively parallel sequencing. Jo Radiat Res. 2013;54:808–22.
Jiang S, Li C, McRae G, Lykken E, Sevilla J, Liu S-Q, Wan Y, Li Q-J. MeCP2 reinforces STAT3 signaling and the generation of effector CD4+ T cells by promoting miR-124-mediated suppression of SOCS5. Sci Signal. 2014;7:ra25.
Pott GB, Chan ED, Dinarello CA, Shapiro L. α-1-Antitrypsin is an endogenous inhibitor of proinflammatory cytokine production in whole blood. J Leukoc Biol. 2009;85:886–95.
Weidner CI, Ziegler P, Hahn M, Brummendorf TH, Ho AD, Dreger P, Wagner W. Epigenetic aging upon allogeneic transplantation: the hematopoietic niche does not affect age-associated DNA methylation. Leukemia. 2015;29:985–8.
Peters MJ, Joehanes R, Pilling LC, Schurmann C, Conneely KN, Powell J, Reinmaa E, Sutphin GL, Zhernakova A, Schramm K, et al. The transcriptional landscape of age in human peripheral blood. Nat Commun. 2015;6:8570.
Ness KK, Krull KR, Jones KE, Mulrooney DA, Armstrong GT, Green DM, Chemaitilly W, Smith WA, Wilson CL, Sklar CA, et al. Physiologic frailty as a sign of accelerated aging among adult survivors of childhood cancer: a report from the St Jude Lifetime Cohort Study. J Clin Oncol. 2013;31:4496–503.
Raghuraman S, Donkin I, Versteyhe S, Barrès R, Simar D. The emerging role of epigenetics in inflammation and immunometabolism. Trends Endocrinol Metab. 2016;27:782–95.
Koturbash I, Boyko A, Rodriguez-Juarez R, McDonald RJ, Tryndyak VP, Kovalchuk I, Pogribny IP, Kovalchuk O. Role of epigenetic effectors in maintenance of the long-term persistent bystander effect in spleen in vivo. Carcinogenesis. 2007;28:1831–8.
Kaup S, Kaup S, Grandjean V, Grandjean V, Mukherjee R, Mukherjee R, Kapoor A, Kapoor A, Keyes E, Keyes E, et al. Radiation-induced genomic instability is associated with DNA methylation changes in cultured human keratinocytes. Mutat Res. 2006;597:87–97.
Jacobsen MJ, Mentzel CMJ, Olesen AS, Huby T, Jorgensen CB, Barrès R, Fredholm M, Simar D. Altered methylation profile of lymphocytes is concordant with perturbation of lipids metabolism and inflammatory response in obesity. J Diabetes Res. 2016;2016:8539057.
Simar D, Versteyhe S, Donkin I, Liu J, Hesson L, Nylander V, Fossum A, Barrès R. DNA methylation is altered in B and NK lymphocytes in obese and type 2 diabetic human. Metabolism. 2014;63:1188–97.
Busslinger M, Tarakhovsky A. Epigenetic control of immunity. Cold Spring Harb Perspect Biol. 2014;6:a019307-a.
Kirchner H, Nylen C, Laber S, Barres R, Yan J, Krook A, Zierath JR, Naslund E. Altered promoter methylation of PDK4, IL1 B, IL6, and TNF after Roux-en Y gastric bypass. Surg Obes Relat Dis. 2014;10:671–8.
Mullen AC, Hutchins AS, High FA, Lee HW, Sykes KJ, Chodosh LA, Reiner SL. Hlx is induced by and genetically interacts with T-bet to promote heritable TH1 gene induction. Nat Immunol. 2002;3:652–8.
Lee DU, Agarwal S, Rao A. Th2 lineage commitment and efficient IL-4 production involves extended demethylation of the IL-4 gene. Immunity. 2002;16:649–60.
Multhoff G, Radons J. Radiation, inflammation, and immune responses in cancer. Front Oncol. 2012;2:58.
Waickman AT, Powell JD. mTOR, metabolism, and the regulation of T-cell differentiation and function. Immunol Rev. 2012;249:43–58.
Rincón M, Davis RJ. Regulation of the immune response by stress-activated protein kinases. Immunol Rev. 2009;228:212–24.
Shan Y-X, Jin S-Z, Liu X-D, Liu Y, Liu S-Z. Ionizing radiation stimulates secretion of pro-inflammatory cytokines: dose–response relationship, mechanisms and implications. Radiat Environ Biophys. 2006;46:21–9.
Kusunoki Y, Yamaoka M, Kubo Y, Hayashi T, Kasagi F, Douple EB, Nakachi K. T-cell immunosenescence and inflammatory response in atomic bomb survivors. Radiat Res. 2010;174:870–6.
López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. The hallmarks of aging. Cell. 2013;153:1194–217.
Weidner CI, Lin Q, Koch CM, Eisele L, Beier F, Ziegler P, Bauerschlag DO, ckel K-HJ, Erbel R, hleisen TWM, et al. Aging of blood can be tracked by DNA methylation changes at just three CpG sites. Genome Biol. 2014;15:R24.
Datta K, Suman S, Fornace J, Albert J. Radiation persistently promoted oxidative stress, activated mTOR via PI3K/Akt, and downregulated autophagy pathway in mouse intestine. Int J Biochem Cell Biol. 2014;57:167–76.
Kuzmina NS, Lapteva NS, Rubanovich AV. Hypermethylation of gene promoters in peripheral blood leukocytes in humans long term after radiation exposure. Environ Res. 2016;146:10–7.
Barrès R, Osler ME, Yan J, Rune A, Fritz T, Caidahl K, Krook A, Zierath JR. Non-CpG methylation of the PGC-1alpha promoter through DNMT3B controls mitochondrial density. Cell Metab. 2009;10:189–98.
Li Y, Gorelik G, Strickland FM, Richardson BC. Oxidative Stress, T cell DNA methylation, and lupus. Arthritis Rheumatol. 2014;66:1574–82.
Wallace TM, Levy JC, Matthews DR. Use and abuse of HOMA modeling. Diabetes Care. 2004;27:1487–95.
Unnikrishnan A, Guan YF, Huang Y, Beck D, Thoms JAI, Peirs S, Knezevic K, Ma S, de Walle IV, de Jong I, et al. A quantitative proteomics approach identifies ETV6 and IKZF1 as new regulators of an ERG-driven transcriptional network. Nucleic Acids Res. 2016;44:10644–61.
Gu H, Smith ZD, Bock C, Boyle P, Gnirke A, Meissner A. Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat Protoc. 2011;6:468–81.
Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27:1571–2.
Akalin A, Kormaksson M, Li S, Garrett-Bakelman FE, Figueroa ME, Melnick A, Mason CE. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13:R87.
Yu G, Wang L-G, He Q-Y. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics. 2015;31:2382–3.
Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.
Andersson R, Gebhard C, Miguel-Escalada I, Hoof I, Bornholdt J, Boyd M, Chen Y, Zhao X, Schmidl C, Suzuki T, et al. An atlas of active enhancers across human cell types and tissues. Nature. 2014;507:455–61.
Ienasescu H, Li K, Andersson R, Vitezic M, Rennie S, Chen Y, Vitting-Seerup K, Lagoni E, Boyd M, Bornholdt J, et al. On-the-fly selection of cell-specific enhancers, genes, miRNAs and proteins across the human body using SlideBase. Database. 2016;2016:baw144.
Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, Diekhans M, Furey TS, Harte RA, Hsu F, et al. The UCSC Genome Browser Database: update 2006. Nucleic Acids Res. 2006;34:D590–D8.
Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote. Nucleic Acids Res. 2013;41:e108-e.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26:139–40.
Eden E, Navon R, Steinfeld I, Lipson D, Yakhini Z. GOrilla: a tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinformatics. 2009;10:48.
The authors are grateful to all the participants who contributed to the study. They also thank Dr. Barbara Cameron for her assistance with the Multiplex assay and Dr. Chris Brownly for his assistance with the flow cytometry and cell sorting experiments at the Biological Resources Imaging Laboratory (Flow Cytometry Facility, UNSW Australia, Australia). We would like to acknowledge The Danish National High-Throughput DNA Sequencing Centre for the sequencing.
This work was supported by a grant from the Cancer Institute New South Wales (09/RFG/2-21) to R.J. Cohn and D. Simar, an International Postgraduate Research Scholarship from UNSW Sydney to S. Daniel, and a research grant from the Danish Diabetes Academy supported by the Novo Nordisk Foundation to O. Fabre. The Novo Nordisk Foundation Center for Basic Metabolic Research is an independent Research Center at the University of Copenhagen partially funded by an unrestricted donation from the Novo Nordisk Foundation.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This project was approved by the Human Research Ethics Committee (UNSW Sydney, HREC-10/017) and performed in accordance with Helsinki Declaration procedures. The nature of the project was explained, and participants provided informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. CD8+ cells polarized activation in childhood cancer survivors (CCS). Figure S2. Intracellular signaling pathways involved in T cell polarization in response to direct and indirect irradiation. Figure S3. Intracellular signaling pathways and polarized activation in CD8+ cells exposed to conditioned media from irradiated adipocytes. Table S1. Differentially Methylated Genes. Table S2. Gene Ontology from differentially methylated genes. Table S3. a. Genes differentially expressed. b. Gene Ontology terms from differentially expressed genes. Table S4. Proteins identified in the supernatant of irradiated fibroblasts. (ZIP 4590 kb)
About this article
Cite this article
Daniel, S., Nylander, V., Ingerslev, L.R. et al. T cell epigenetic remodeling and accelerated epigenetic aging are linked to long-term immune alterations in childhood cancer survivors. Clin Epigenet 10, 138 (2018) doi:10.1186/s13148-018-0561-5
- Epigenetic aging
- DNA methylation
- cancer survivors
- T cell