Skip to main content


DNA methylation profiles of elderly individuals subjected to indentured childhood labor and trauma

Article metrics



Childhood trauma is associated with increased vulnerability to mental and somatic disorders later in life. Epigenetic modifications such as DNA methylation are one potential mechanism through which such long-lasting impairments/consequences can be explained. The aim of the present study was to investigate whether childhood trauma is associated with long-term DNA methylation alterations in old age.


We assessed genome-wide DNA methylation profiles in a cohort of former indentured child laborers (“Verdingkinder”) who suffered severe childhood adversities (N = 30; M age = 75.9 years), and compared them to control group with similar demographic characteristics (N = 15, M age = 72.8 years). DNA was isolated from epithelial buccal cells and hybridized to the Illumina Infinium 450 k DNA methylation array, which provides coverage of 485,000 methylation sites.


After accounting for batch effects, age, gender and multiple testing, 71 differentially methylated CpG positions were identified between the two groups. They were annotated among others to genes involved in neuronal projections and neuronal development. Some of the identified genes with differential methylation (DLG associated protein 2, mechanistic target of rapamycin) have previously been associated with traumatic stress.


The results indicate specific epigenetic alterations in elderly individuals who were subjected to childhood adversities. Psychiatric and somatic comorbidities as well as differences in buccal epithelial cells proportion may contribute to the observed epigenetic differences.


Childhood adversities are highly prevalent [1, 2]. In the World Mental Health surveys, prevalence of childhood maltreatment - one type of childhood adversity - has been estimated at 8% for physical abuse, 4.4% for neglect, and 1.6% for sexual abuse [2]. Exposure to childhood adversities is associated with increased vulnerability to a number of mental disorders in adulthood, including major depression, anxiety, suicidal ideation, and posttraumatic stress disorder (PTSD) [35]. In addition, childhood trauma has been linked to higher incidence of general medical conditions and increased disability [3]. While it has been consistently demonstrated that childhood adversities are associated with long-term health consequences, the molecular mechanisms through which this occurs remain unknown. It has been hypothesized that a combination of genetic predisposition, epigenetic modifications, and changes in immune and hormonal parameters may be the mechanisms through which childhood trauma can lead to an increased risk of morbidity later in life [6].

Epigenetic modifications include, among others, DNA methylation [7]. So far only few studies have explored DNA methylation differences between victims of childhood trauma and controls on the whole genome level [813]. In an attempt to do so, Yang et al. identified differential methylation in 2868 CpG sites in DNA from saliva of maltreated children (n = 96) and demographically matched controls (n = 96) [11]. Genes relevant to biological processes such as neurogenesis and axonal guidance, as well as a number of disease biomarkers differed in their DNA methylation. In another study, Suderman et al. assessed whole blood DNA methylation profiles of 45 years old individuals subjected to childhood abuse (n = 12) and compared these to the methylation profiles of healthy controls (n = 28) [10]. Here, the authors found DNA methylation alterations in 997 promoters, mostly related to developmental and regulational functions, and enrichment in genes associated with the WNT signaling pathway. Yet another study looking at children who were institutionalized at birth (n = 14) and comparing them to non-institutionalized controls (n = 14), found differential methylation in 914 CpG sites located in genes involved in cell signaling, inflammation, but also in neuronal development and communication [12]. Some studies, however, failed to detect significant DNA methylation differences, such as the one conducted by Smith and colleagues who explored peripheral blood DNA methylation profiles in 110 African American individuals, stratified according to childhood trauma and PTSD diagnosis. In that investigation, no significant epigenome-wide differences could be found as a result of childhood abuse. However, CpG sites in five genes, mainly related to immune dysregulation, were differentially methylated in subjects suffering from PTSD in comparison to subjects without PTSD [13].

Overall, previous research has provided evidence for differentially methylated regions associated with previous traumatization and potentially explaining later morbidity. However, most of these studies have been carried out in children or young/middle aged adults, although there is evidence, that long-term epigenetic alterations can also be observed in older individuals. In a study conducted on elderly Dutch individuals subjected to another type of early-life stress (prenatal famine) during World War II and comparing them to their siblings, Tobi et al. found evidence for differential DNA methylation extending to regulatory regions, as well pathways related to growth and metabolism [14]. In contrast, in an investigation of elderly Finnish individuals who experienced separation from their families during the Second World War (n = 83) and demographically matched controls (n = 83), experiment-wide differences in blood DNA methylation profiles were found only in association with depressive symptoms and not with early life separation [15].

To further investigate the long-term epigenetic changes after childhood traumatization that might explain the higher risk for morbidity in later life, we assessed genome-wide DNA methylation profiles in a group of elderly, previously indentured Swiss child laborers (i.e., “Verdingkinder”) and compared them to a group of controls, consisting of individuals with similar age and gender ratio. The “Verdingkinder” were children who up until the 1950s were removed from their biological families due to different reasons (for example poverty, being born out of wedlock or being orphaned). They were sent to live and work with farmers and many of these children experienced physical, emotional and/or sexual abuse during the period of indentured labor [16, 17]. Therefore, they represent an ideal cohort to be studied in terms of long-term consequences of early childhood trauma.


Study group

The “case” sample of former indentured child laborers belonged to a larger study group from a longitudinal project, aiming at evaluating a range of physical and mental health parameters in traumatized individuals [1820]. For DNA methylation comparison, a control group with similar demographic characteristics was additionally recruited. The study was approved by the Cantonal Ethic Commission of Zurich – study number KEK-ZH-Nr. 2012–0245. Informed consent was obtained by all participants in the study. Inclusion criteria for both the former indentured child labor group and the control group were: at least 65 years of age, voluntary participation, (Swiss) German speaking, upbringing in rural areas of Switzerland. Specific inclusion criteria for the former indentured child labor group further included: at least one period of indentured child labor. Specific inclusion criteria for the control group were: no indentured child labor or adoptive care, no major childhood trauma.

Assessment of childhood trauma

Former child laborers were asked about their age at the beginning and end of the indentured child labor experience, as well as the number of different families in which they were placed. To evaluate the presence and extent of childhood trauma, the childhood trauma questionnaire – short form (CTQ-SF) was used in both child laborers and the controls [21, 22]. The CTQ-SF is a retrospective self-report 28-item questionnaire (consisting of 25 clinical items and 3 validity items). Subscales scores for emotional abuse, physical abuse, sexual abuse, emotional neglect and physical neglect are included, each consisting of 5 items with severity scores ranging from 5 to 25. In addition, a summary score (ranging from 25 to 125) for total severity of childhood abuse and neglect can be computed by summing up all 5 subscales. The CTQ-SF has proven to be a reliable measure of childhood abuse and neglect in adult samples in a number of validation studies. A recent validation study of a German version of the CTQ-SF conducted in a sample of Swiss patients and non-clinical participants [23] confirmed the original five-factor model of the original English version and further found high internal consistencies for all scales (physical abuse =0.83, emotional abuse =0.83, sexual abuse =0.9, emotional neglect =0.91), apart from physical neglect (α = 0.53).

DNA methylation array

Epithelial buccal cells were sampled using Isohelix DNA Buccal Swabs, which have been shown to yield a relatively homogenous population of squamous epithelial cells, with a smaller proportion of blood cells [24]. DNA on the swab heads was stabilized with Isohelix Dri-Capsules (Cell Projects Ltd, Kent, United Kingdom). None of the individuals from whom buccal swabs were collected showed any signs of oral infection, which could lead to higher cellular heterogeneity. Buccal swabs were sent to the “Barts and The London Genome Centre” in the United Kingdom for further processing. DNA was isolated from the buccal swabs using the Isohelix buccal DNA isolation kit (Cell Projects Ltd) and purified with the Zymo ZR-96 DNA clean-up kit (Zymo Research Corporation, Irvine, United States). Before and after purification, DNA concentration was measured with the Qubit 2.0 Fluorometer (Life Technologies). DNA integrity was estimated with agarose gel electrophoresis. After DNA bisulfite conversion using the EZ DNA Methylation kit (Zymo Research Corporation), DNA concentration was again controlled for using Nanodrop (Thermo Scientific). Bisulfite converted DNA was hybridized to the Illumina Infinium 450 k DNA methylation array providing coverage of 485,000 methylation sites for each analyzed sample [25]. The array’s coverage extends to 99% of RefSeq genes, including multiple gene probes, 96% of CpG islands, CpG island shores, and additional content [26].

DNA methylation data analysis

Raw intensity data files (IDAT files) served as an input and were analyzed in the R environment ( [27]. Data import, quality control and normalization steps of the minfi Bioconductor package were used [28, 29]. Strict quality criteria were introduced in order to decrease the variability of the data set and to allow more reliable detection of true intergroup differences. Samples with fraction of failed probes per sample > 0.05 (one sample) or samples, which in spot-checks of mean fluorescence showed values above a threshold of 0.01 (7 samples), were excluded from the analysis. Probes exhibiting detection p-values > 0.01 in one or more samples were excluded (n = 451). Background correction and normalization were carried out with the Illumina method implemented in minfi (bg.correct = TRUE, normalize = “controls”). The probes on the Illumina Infinium 450 k array measuring single nucleotide polymorphisms (SNPs) (n = 65) or non-CpG methylation (n = 3026) were filtered out. Probes containing SNPs at the target CpG site or at single base extension position were excluded and SNPs with any minor allele frequency were considered (n = 17534). The SNP annotation integrated in the minfi package was used [30]. Finally, probes on the sex chromosomes were also filtered out (n = 11361). The minfi preprocessed data were analyzed with CpGassoc package to assess differences in DNA methylation profiles between the study groups [31]. A linear fixed effects model was applied for each individual CpG site with age, gender and batch included as covariates in the analysis. The false discovery rate (FDR) calculated according to the Benjamini-Hochberg method was set to 0.05. In a secondary set of analyses major depressive disorder (MDD), positive PTSD screen, coronary heart disease, diabetes, length of education and estimated proportion of buccal epithelial cells were also included as covariates. The proportion of buccal epithelial cells versus leukocytes in the buccal swabs was estimated according to a previously published method [32]. Finally, the data were checked for the possible effect of methylation outliers through visually inspecting CpG plots and including the minfi function fixMethOutliers (object, K = −3, verbose = FALSE) in the analysis. 2 of the top reported 71 associations from the original analysis were excluded by controlling for outliers (cg16669619 and cg14895646) and they were not mapped to genes included in the downstream gene enrichment analysis, suggesting that methylation outliers did not affect the observed results. For statistical analysis, M-values, which represent log2-ratios of methylated and unmethylated probes intensities and have been shown to have higher statistical validity, were used [33]. The M-value method has been shown to have superior performance in detection rate and true positive rate of both unmethylated and highly methylated CpG sites [33]. For graphic representation of the data, beta values were employed due to their more intuitive visual interpretation.

Gene enrichment analysis

Genes in the former indentured child labor group showing significant CpG site differential methylation were analyzed for annotation to biological pathways and their overlap using the Database for Annotation, Visualization and Integration Discovery (DAVID) [34], ToppGene Suite [35] and Pathway Studio [36]. The aim of the analysis was to explore potential biological meaning, even in cases of limited statistical power. DAVID database functional annotation clustering was performed for genes with differential methylation between the two groups. Clusters with enrichment score (ES) > 1.5 and p < 0.05 were reported, which is more conservative relative to the generally suggested threshold enrichment score of > 1.3 [34]. Genes from the human genome were used as background. In the ToppGene Suite, the ToppFun module was used to analyze functional enrichment in biological processes of differentially methylated genes. Enrichment for genes from input ≥ 3 and FDR ≤ 0.05 after Benjamini-Hochberg correction for multiple testing is reported. In Pathway Studio gene set enrichment analysis was carried out with a cut-off of for genes from input ≥ 3 and p < 0.05. No correction was applied for the number of CpG sites genes have measured on the microarray.


Former child laborers and controls did not differ significantly based on age (mean ± SD 75.9 ± 5.8 years for the former child labor group compared to 72.8 ± 5.8 years for the control group) and gender (14 F/16 M for the former child labor group compared to 8 F/7 M for the control group). There were no significant differences between the two groups in regards to household income, alcohol or substance dependence/abuse, blood pressure, and/or cholesterol levels. PTSD symptoms, coronary heart disease and diabetes were more frequent in the former child labor group, in accordance with the expected higher prevalence of these disorders after childhood trauma. The control group reported more years of education. Positive depressive symptoms screen was more commonly observed in the former child labor group, but the two groups did not differ significantly in the incidence of MDD diagnosis. Smoking status was known for only a subset of the individuals involved in the study (over 90% of both former indentured child laborers and controls for whom the smoking status was known were non-smokers). We assessed the methylation levels of cg05575921 (annotated to the aryl hydrocarbon receptor repressor), whose methylation was recently found to be a sensitive and specific marker for smoking status in adults [37]. Mean DNA methylation levels of cg05575921 did not differ between the two groups (p-value 0.61), implying similarity of the smoking status between the two groups as a whole. Medication status was also known only for a subset of the individuals included in the study, but due to the advanced age of the study population over 80% of participants with known status in both groups were taking a medication. Childhood trauma severity scores were significantly higher in the former child labor group in comparison to the control group. Demographic and childhood trauma characteristics of the participants are presented in Table 1 and Additional file 1: Table S1.

Table 1 Summary of the demographic and childhood trauma characteristics of the subjects included in the study

Differential methylation between former indentured child laborers and controls

We found 71 differentially methylated CpG positions between the former indentured child laborer and control group, after accounting for age, gender, batch effects, and multiple comparisons. The list of differentially methylated positions meeting experiment-wide significance set at FDR < 0.05 is presented in Table 2, while β-values in the two groups for these positions are presented in Additional file 1: Table S3. One gene – the src kinase associated phosphoprotein 2 (SKAP2) - included 5 differentially methylated CpG positions; 41 genes included one differentially methylated position; and 25 differentially methylated positions were located in intergenic regions. Box-plots presenting beta levels between the two groups in the four most significant differentially methylated positions are presented in Fig. 1.

Table 2 List of differentially methylated CpG positions between the two groups
Fig. 1

Box plots representing the top four differentially methylated positions between the former indentured child labor and control groups

Enrichment analysis of differentially methylated genes between the former indentured child labor and control groups

Functional annotation clustering of genes with differentially methylated CpG sites using the DAVID database revealed two clusters (Table 3). One of them included genes related to the regulation of cell projection and cellular component organization (ES 1.74, genes: MTOR; PLXNB1; ROBO1), while the second involved genes related to neuronal projection/dendrites (ES 1.65, genes: SH3PXD2A; CACNA1H; DLGAP2; DNAH1; MTOR; ROBO1). Using the ToppGene Suite, we found enrichment of genes with differentially methylated CpG sites implicated among others in telencephalon development (ROBO1; DLX1; SKI; GNG12; OGDH) and olfactory bulb development (ROBO1; SKI; OGDH) (Additional file 1: Table S2). Overlap of Pathway Studio analysis data with the other gene enrichment tools included enrichment in genes related to multicellular organismal development (DLX1; DLX5; ESRRG; PLXNB1; ROBO1), dendrites (CACNA1H; DLGAP2; CNTNAP4; MTOR) and axon guidance (CACNA1H; PLXNB1; ROBO1; DLX5).

Table 3 Functional annotation cluster analysis of genes differentially methylated in the former indentured child labor group versus the control group according to DAVID with p < 0.05

Effect of psychiatric and somatic comorbidities, education duration and proportion of buccal epithelial cells on differentially methylated genes between the groups

Finally, we assessed the role of some major psychiatric (MDD, positive PTSD screen) and somatic (coronary heart disease, diabetes) comorbidities, education and estimated proportion of buccal epithelial cells on the most significant DNA methylation differences between the former child labor and the control group. DNA methylation differences were evaluated between the four most significantly differentially methylated CpG sites in the primary analysis, as well as CpG sites in genes showing consistent enrichment with the applied gene enrichment tools. We included MDD, positive PTSD screen, coronary heart disease, diabetes, education duration and estimated proportion of buccal epithelial cells as covariates. P-values for the assessed CpG sites varied between 9.47 × 10−6 and 0.00193 (Tables 4 and 5).

Table 4 Effect of adjusting for psychiatric and somatic comorbidities, length of education and estimated proportion of buccal epithelial cells on the top four differentially methylated CpG positions between the two groups
Table 5 Effect of adjusting for psychiatric and somatic comorbidities, length of education and estimated proportion of buccal epithelial cells on differentially methylated CpG positions between the two groups implicated consistently by the applied gene enrichment tools


In the present study we assessed differences in DNA methylation patterns between former indentured child laborers and controls with similar demographic characteristics. Overall, we found DNA methylation differences in 71 CpG positions meeting the experiment-wide significance criteria of FDR < 0.05 between the two groups. Psychiatric and somatic comorbidities and estimated proportion of epithelial cells contribiute partially to the detected differences.

The gene that showed the strongest difference in methylation patterns was SKAP2, which included 5 differentially methylated positions. SKAP2 is an adaptor protein that plays an important role in src signaling and is involved in a wide range of intracellular processes, such as suppression of cell migration and tumor invasion by inhibition of actin polymerization [3840]. The gene has also been shown to negatively affect the phosphorylation of alpha-synuclein - a protein critically involved in neurodegenerative disorders such as Parkinson’s disease - after cellular stress [41]. Interestingly, the src-kinase for which SKAP2 is an adaptor protein, converges with the WNT pathway recently implicated in DNA methylation alterations in 45 years old individuals subjected to childhood maltreatment [10, 42]. In our study, all of the 5 detected SKAP2 differentially methylated positions were located in the gene body and all of them showed hypermethylation in the former indentured child labor group. DNA methylation in gene promoters has long been shown to interfere with gene transcription [7]. A role of DNA methylation in gene bodies for gene expression has been suggested, but the exact mechanism and effect are still not conclusively clarified [43, 44]. Since we found methylation changes in CpG sites located both in promotor and gene body regions (Table 1), the effect on gene and protein expression cannot be directly predicted and warrants further investigation.

We further found enrichment of genes with differentially methylated CpG sites in cell/neuronal projection (CACNA1H; DLGAP2; MTOR; ROBO1; SH3PXD2A), dendrites (CACNA1H; DLGAP2; CNTNAP4; MTOR), axon guidance (CACNA1H; PLXNB1; ROBO1; DLX5), brain development (ROBO1; DLX1; SKI; GNG12; OGDH) and multicellular organismal development (DLX1; DLX5; ESRRG; PLXNB1; ROBO1). Our findings are in line with previous studies also reporting links of childhood trauma with genes related to cell signaling and neuronal development [811]. In particular, DLG associated protein 2 (DLGAP2) and mechanistic target of rapamycin (MTOR) have previously been associated with traumatic stress. In a rat model of PTSD increased DLGAP2 DNA methylation levels and decreased mRNA expression levels have been observed [45]. MTOR has been shown to regulate fear memory reconsolidation [46] and its activation in a rat PTSD model has been demonstrated [47]. Our study focused on elderly individuals subjected to childhood trauma, while most of previous investigations have focused on young to middle aged adults [48]. The aging process itself is associated with DNA methylation alterations and traumatic stress can accelerate epigenetic aging [49]. However, it has been shown that some early-life induced DNA methylation changes are protected from erasure associated with age [48]. The extent to which age-associated DNA methylation alterations obliterate the ones induced by early-life trauma in our study population is not possible to predict.

About 30% of the differentially methylated CpG positions detected were located in intergenic regions. A previous study by Yang et al. investigating abused children found about 20% of the differentially methylated positions to be in intergenic regions [11]. Changes in DNA methylation in intergenic regions may play an important role in the genome organization and gene expression regulation through the binding sites for transcription factors located in them [50].

Our study has several limitations. These include the relatively small sample size and the cross-sectional design of the study. Due to the old age of the participants, comorbidities and medication intake were very common, making controlling for potential confounding factors difficult. Even though we assessed the role of certain psychiatric and somatic comorbidities, additional controls and replication in larger study cohorts are needed. In addition, we investigated a unique group of individuals subjected to prolonged and severe traumatic experience in their childhood, and it may not be possible to extrapolate detected differences directly to other study populations. Finally, we used a peripheral tissue – buccal epithelial cells for detection of DNA methylation patterns. While epigenetic patterns are tissue-specific, DNA methylation alterations in response to early life adversity have been shown to be a system-wide phenomenon, and studying peripheral tissues may yield important clues to the biological pathways alterations induced [51].


Altogether, we found differences in DNA methylation profiles of elderly individuals subjected to prolonged and severe childhood trauma in comparison to controls with similar demographic characteristics. They encompass an adaptor protein for the src-kinase, genes related to development, cellular and neuronal projection. Other psychiatric and somatic comorbidities and estimated proportion of buccal epithelial cells appear to partially contribute to the here observed DNA methylation differences.



3’ untranslated region


5’ untranslated region


Childhood trauma questionnaire – short form


Database for annotation, visualization and integration discovery


DLG associated protein 2


Enrichment score




False discovery rate

IDAT file:

Raw intensity data file




Major depressive disorder


Mechanistic target of rapamycin


Posttraumatic stress disorder


Standard deviation


Src kinase associated phosphoprotein 2


Single nucleotide polymorphism


Within 1500 bps of a transcription start site


Within 200 bps of a transcription start site


  1. 1.

    Boyce P, Harris A. Childhood adversity, trauma and abuse: context and consequences. Aust N Z J of Psychiatry. 2011;45:608–10.

  2. 2.

    Kessler RC, McLaughlin KA, Green JG, Gruber MJ, Sampson NA, Zaslavsky AM, Aguilar-Gaxiola S, Alhamzawi AO, Alonso J, Angermeyer M, Benjet C, Bromet E, Chatterji S, de Girolamo G, Demyttenaere K, Fayyad J, Florescu S, Gal G, Gureje O, Haro JM, Hu CY, Karam EG, Kawakami N, Lee S, Lépine JP, Ormel J, Posada-Villa J, Sagar R, Tsang A, Ustün TB, Vassilev S, Viana MC, Williams DR. Childhood adversities and adult psychopathology in the WHO World Mental Health Surveys. Br J Psychiatry. 2010;197:378–85.

  3. 3.

    Cuijpers P, Smit F, Unger F, Stikkelbroek Y, ten Have M, de Graaf R. The disease burden of childhood adversities in adults: a population-based study. Child Abuse Negl. 2011;35:937–45.

  4. 4.

    Green JG, McLaughlin KA, Berglund PA, Gruber MJ, Sampson NA, Zaslavsky AM, Kessler RC. Childhood adversities and adult psychiatric disorders in the national comorbidity survey replication I: associations with first onset of DSM-IV disorders. Arch Gen Psychiatry. 2010;67:113–23.

  5. 5.

    McLaughlin KA, Green JG, Gruber MJ, Sampson NA, Zaslavsky AM, Kessler RC. Childhood adversities and adult psychiatric disorders in the national comorbidity survey replication II: associations with persistence of DSM-IV disorders. Arch Gen Psychiatry. 2010;67:124–32.

  6. 6.

    Ehlert U. Enduring psychobiological effects of childhood adversity. Psychoneuroendocrinology. 2013;38:1850–7.

  7. 7.

    Klose RJ, Bird AP. Genomic DNA methylation: the mark and its mediators. Trends Biochem Sci. 2006;31:89–97.

  8. 8.

    Labonté B, Suderman M, Maussion G, Navaro L, Yerko V, Mahar I, Bureau A, Mechawar N, Szyf M, Meaney MJ, Turecki G. Genome-wide epigenetic regulation by early-life trauma. Arch Gen Psychiatry. 2012;69:722–31.

  9. 9.

    Mehta D, Klengel T, Conneely KN, Smith AK, Altmann A, Pace TW, Rex-Haffner M, Loeschner A, Gonik M, Mercer KB, Bradley B, Müller-Myhsok B, Ressler KJ, Binder EB. Childhood maltreatment is associated with distinct genomic and epigenetic profiles in posttraumatic stress disorder. Proc Natl Acad Sci U S A. 2013;110:8302–7.

  10. 10.

    Suderman M, Borghol N, Pappas JJ, Pinto Pereira SM, Pembrey M, Hertzman C, Power C, Szyf M. Childhood abuse is associated with methylation of multiple loci in adult DNA. BMC Med Genomics. 2014;7:13.

  11. 11.

    Yang BZ, Zhang H, Ge W, Weder N, Douglas-Palumberi H, Perepletchikova F, Gelernter J, Kaufman J. Child abuse and epigenetic mechanisms of disease risk. Am J Prev Med. 2013;44:101–7.

  12. 12.

    Naumova OY, Lee M, Koposov R, Szyf M, Dozier M, Grigorenko EL. Differential patterns of whole-genome DNA methylation in institutionalized children and children raised by their biological parents. Dev Psychopathol. 2012;24:143–55.

  13. 13.

    Smith AK, Conneely KN, Kilaru V, Mercer KB, Weiss TE, Bradley B, Tang Y, Gillespie CF, Cubells JF, Ressler KJ. Differential immune system DNA methylation and cytokine regulation in post-traumatic stress disorder. Am J Med Genet B. 2011;156:700–8.

  14. 14.

    Tobi EW, Goeman JJ, Monajemi R, Gu H, Putter H, Zhang Y, Slieker RC, Stok AP, Thijssen PE, Müller F, van Zwet EW, Bock C, Meissner A, Lumey LH, Eline Slagboom P, Heijmans BT. DNA methylation signatures link prenatal famine exposure to growth and metabolism. Nat Commun. 2014;5:5592.

  15. 15.

    Khulan B, Manning JR, Dunbar DR, Seckl JR, Raikkonen K, Eriksson JG, Drake AJ. Epigenomic profiling of men exposed to early-life stress reveals DNA methylation differences in association with current mental state. Transl Psychiatry. 2014;4:e448.

  16. 16.

    Leuenberger M. Versorgt und vergessen: Ehemalige Verdingkinder erzählen. Zurich: Rotpunktverlag. 2008.

  17. 17.

    Wohlwend L, Honeggerm A. Gestohlene Seelen: Verdingkinder in der Schweiz. Frauenfeld: Huber; 2004.

  18. 18.

    Burri A, Maercker A, Krammer S, Simmen-Janevska K. Childhood trauma and PTSD symptoms increase the risk of cognitive impairment in a sample of former indentured child laborers in old age. PLoS One. 2013;8:e57826.

  19. 19.

    Kuhlman KR, Maercker A, Bachem R, Simmen K, Burri A. Developmental and contextual factors in the role of severe childhood trauma in geriatric depression: the sample case of former indentured child laborers. Child Abuse Negl. 2013;37:969–78.

  20. 20.

    Maercker A, Hilpert P, Burri A. Childhood trauma and resilience in old age: applying a context model of resilience to a sample of former indentured child laborers. Aging Ment Health. 2016;20:616–26.

  21. 21.

    Bernstein DP, Fink L. Childhood trauma questionnaire: a retrospective self-report manual. San Antonio: Harcourt Brace & Company; 1998.

  22. 22.

    Bernstein DP, Stein JA, Newcomb MD, Walker E, Pogge D, Ahluvalia T, Stokes J, Handelsman L, Medrano M, Desmond D, Zule W. Development and validation of a brief screening version of the childhood trauma questionnaire. Child Abuse Negl. 2003;27:169–90.

  23. 23.

    Karos K, Niederstrasser N, Abidi L, Bernstein DP, Bader K. Factor structure, reliability, and known groups validity of the German version of the Childhood Trauma Questionnaire (Short-form) in Swiss patients and nonpatients. J Child Sex Abus. 2014;23:418–30.

  24. 24.

    Martino D, Loke YJ, Gordon L, Ollikainen M, Cruickshank MN, Saffery R, Craig JM. Longitudinal, genome-scale analysis of DNA methylation in twins from birth to 18 months of age reveals rapid epigenetic change in early life and pair-specific effects of discordance. Genome Biol. 2013;14:R42.

  25. 25.

    Jones MJ, Farré P, McEwen LM, Macisaac JL, Watt K, Neumann SM, Emberly E, Cynader MS, Virji-Babul N, Kobor MS. Distinct DNA methylation patterns of cognitive impairment and trisomy 21 in Down syndrome. BMC Med Genet. 2013;6:58.

  26. 26.

    Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, Fan JB, Shen R. High density DNA methylation array with single CpG site resolution. Genomics. 2011;98:288–95.

  27. 27.

    R_Development_Core_Team. R. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.

  28. 28.

    Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.

  29. 29.

    Hansen KD, Ayree M, Irizary RA, Jaffe AE, Maksimovic J, Houseman A, Fortin JP, Triche T. Package “minfi”: Analyze Illumina’s 450 k methylation arrays. R package version 1.10.2. 2014.

  30. 30.

    Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.

  31. 31.

    Barfield RT, Kilaru V, Smith AK, Conneely KN. CpGassoc: an R function for analysis of DNA methylation microarray data. Bioinformatics. 2012;28:1280–1.

  32. 32.

    Eipel M, Mayer F, Arent T, Ferreira MR, Birkhofer C, Gerstenmaier U, Costa IG, Ritz-Timme S, Wagner W. Epigenetic age predictions based on buccal swabs are more precise in combination with cell type-specific DNA methylation signatures. Aging (Albany NY). 2016;8:1034–48.

  33. 33.

    Du P, Zhang X, Huang CC, Jafari N, Kibbe WA, Hou L, Lin SM. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinf. 2010;11:587.

  34. 34.

    da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

  35. 35.

    Chen J, Bardes EE, Aronow BJ, Jegga AG. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009;37:W305–11.

  36. 36.

    Nikitin A, Egorov S, Daraselia N, Mazo I. Pathway studio-the analysis and navigation of molecular networks. Bioinformatics. 2003;19:2155–7.

  37. 37.

    Philibert R, Hollenbeck N, Andersen E, Osborn T, Gerrard M, Gibbons FX, Wang K. A quantitative epigenetic approach for the assessment of cigarette consumption. Front Psychol. 2015;6:656.

  38. 38.

    Kouroku Y, Soyama A, Fujita E, Urase K, Tsukahara T, Momoi T. RA70 is a src kinase-associated protein expressed ubiquitously. Biochem Biophys Res Commun. 1998;252:738–42.

  39. 39.

    Marie-Cardine A, Verhagen AM, Eckerskorn C, Schraven B. SKAP-HOM, a novel adaptor protein homologous to the FYN-associated protein SKAP55. FEBS Lett. 1998;435:55–60.

  40. 40.

    Shimamura S, Sasaki K, Tanaka M. The Src substrate SKAP2 regulates actin assembly by interacting with WAVE2 and cortactin proteins. J Biol Chem. 2013;288:1171–83.

  41. 41.

    Takahashi T, Yamashita H, Nagano Y, et al. Identification and characterization of a novel Pyk2/related adhesion focal tyrosine kinase-associated protein that inhibits alpha-synuclein phosphorylation. J Biol Chem. 2003;278:42225–33.

  42. 42.

    Taylor S, Bagrodia S. Src and Wnt converge to seal cell’s fate. Mol Cell. 2002;10:10–1.

  43. 43.

    Ball MP, Li JB, Gao Y, Lee JH, LeProust EM, Park IH, Xie B, Daley GQ, GM C. Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat Biotechnol. 2009;27:361–8.

  44. 44.

    Lou S, Lee HM, Qin H, Li JW, Gao Z, Liu X, Chan LL, Kl Lam V, So WY, Wang Y, Lok S, Wang J, Ma RC, Tsui SK, Chan JC, Chan TF, Yip KY. Whole-genome bisulfite sequencing of multiple individuals reveals complementary roles of promoter and gene body methylation in transcriptional regulation. Genome Biol. 2014;15:408.

  45. 45.

    Chertkow-Deutsher Y, Cohen H, Klein E, Ben-Shachar D. DNA methylation in vulnerability to post-traumatic stress in rats: evidence for the role of the post-synaptic density protein Dlgap2. Int J Neuropsychopharmacol. 2010;13:347–59.

  46. 46.

    Blundell J, Kouser M, Powell CM. Systemic inhibition of mammalian target of rapamycin inhibits fear memory reconsolidation. Neurobiol Learn Mem. 2008;90:28–35.

  47. 47.

    Sui ZX, Liu H, Wang HT, Yang ZL, Liu JG, Xu AJ, Lu HX. Alteration of apoptosis and Akt/mTOR signal pathway in hippocampal neurons of rat with post-traumatic stress. Sichuan Da Xue Xue Bao Yi Xue Ban. 2014;45:221–4.

  48. 48.

    Murgatroyd C, Wu Y, Bockmühl Y, Spengler D. The Janus face of DNA methylation in aging. Aging (Albany NY). 2010;2:107–10.

  49. 49.

    Boks MP, van Mierlo HC, Rutten BP, Radstake TR, De Witte L, Geuze E, Horvath S, Schalkwyk LC, Vinkers CH, Broen JC, Vermetten E. Longitudinal changes of telomere length and epigenetic age related to traumatic stress and post-traumatic stress disorder. Psychoneuroendocrinology. 2015;51:506–12.

  50. 50.

    Hodges E, Molaro A, Dos Santos CO, Thekkat P, Song Q, Uren PJ, Park J, Butler J, Rafii S, McCombie WR, Smith AD, Hannon GJ. Directional DNA methylation changes and complex intermediate states accompany lineage specificity in the adult hematopoietic compartment. Mol Cell. 2011;44:17–28.

  51. 51.

    Szyf M, Bick J. DNA methylation: a mechanism for embedding early life experiences in the genome. Child Dev. 2013;84:49–57.

Download references


We thank all individuals who participated in the current study. We thank the “Barts and The London Genome Centre” in the United Kingdom for performing the Illumina Infinium 450 k DNA methylation arrays experiments. Financial support for the study was provided by the “Hedwig Widmer Foundation”. AB reports an Ambizione personal career fellowship from the Swiss National Science Foundation, as well as a project grant from the Swisslife Jubiläumsstiftung.

Availability of data and materials

The ethical permission for the study does not include a permission to make underlying patient-specific data publicly available.

Authors’ contributions

AM and AB conceived and designed the study with input from SW, EG and ZM. AM, AB, AK, ZM, EG and SW were involved in the process of patients’ recruitment and biological samples collection. ZM, MR, TW, EG, AB were involved in data analysis and interpretation with input from AM, SW and AK. ZM drafted the manuscript, which was revised and approved by all authors.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

This study was approved by Cantonal Ethic Commission of Zurich. Informed consent was obtained from all participants.

Author information

Correspondence to Zoya Marinova.

Additional file

Additional file 1:

Supplementary data. (DOCX 29 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Marinova, Z., Maercker, A., Küffer, A. et al. DNA methylation profiles of elderly individuals subjected to indentured childhood labor and trauma. BMC Med Genet 18, 21 (2017) doi:10.1186/s12881-017-0370-2

Download citation


  • Childhood Trauma
  • Buccal Epithelial Cell
  • Somatic Comorbidities
  • ToppGene Suite
  • Positive Ptsd Screen