The epigenetic clock and pubertal, neuroendocrine, psychiatric, and cognitive outcomes in adolescents
Clinical Epigenetics volume 10, Article number: 96 (2018)
Molecular aging biomarkers, such as epigenetic age predictors, predict risk factors of premature aging, and morbidity/mortality more accurately than chronological age in middle-aged and elderly populations. Yet, it remains elusive if such biomarkers are associated with aging-related outcomes earlier in life when individuals begin to diverge in aging trajectories. We tested if the Horvath epigenetic age predictor is associated with pubertal, neuroendocrine, psychiatric, and cognitive aging-related outcomes in a sample of 239 adolescents, 11.0–13.2 years-old.
Each year increase in epigenetic age acceleration (AA) was associated with 0.06 SD units higher weight-for-age, 0.08 SD units taller height-for-age, -0.09 SD units less missed from the expected adult height, 13 and 16% higher odds, respectively, for each stage increase in breast/genitals development on the Tanner Staging Questionnaire and pubertal stage on the Pubertal Development Scale, 4.2% higher salivary cortisol upon awakening, and 18 to 34% higher odds for internalizing and thought problems on the Child Behavior Checklist (p values < 0.045). AA was not significantly associated with cognition.
Our findings suggest that already in adolescence, AA is associated with physiological age acceleration, which may index risk of earlier aging. AA may identify individuals for preventive interventions decades before aging-related diseases become manifest.
There is an imperative to discover biomarkers of human aging, driven by demographic shifts and the striking inter-individual variation in trajectories in aging [1, 2]. Biomarkers of aging would allow more precise quantification of true biological age, permitting the identification of individuals at risk of aging-related diseases for preventive interventions decades before disorders become manifest.
Among the most promising molecular biomarkers of aging are those based on changes in DNA methylation (DNAm), the covalent addition of a methyl group primarily to cytosine linked to guanine by phosphate (CpG) sites . Variation in DNAm of 71 CpG sites in the whole blood  and of 353 CpG sites from multiple tissues and cell types  have been identified to predict chronological age with high accuracy (r > .91); note that the algorithms involve both increased and decreased methylation of CpG sites to predict aging. These, Hannum and Horvath epigenetic age predictors, demonstrate a median absolute difference between DNAm age and actual chronological age of up to 4.9  and 3.5 years  and have been validated in 19–101  and 0–100-year-old  individuals, respectively. In middle-aged and elderly individuals, the Horvath and Hannum-based measures of epigenetic age acceleration (AA) (higher DNAm age than actual chronological age) were associated with higher body-mass index , lower physical and cognitive fitness , Alzheimer’s disease , menopause , and increased risk of all-cause mortality [10, 11].
Studies in middle-aged to elderly populations are, however, confounded by the often decade-long processes of aging-related disease and aging in itself. Therefore, studies of aging might better focus earlier in life, when inter-individual differences in aging trajectories start to emerge, but before most age-related diseases become manifest . Such studies focusing on AA early in life are scarce. In one study, which tested associations between the Horvath epigenetic age predictor at birth, 7 and 17 years and physical growth and development among 400 to 1000 UK children, found that higher AA at birth predicted higher fat mass in childhood and adolescence, faster growth in weight and body mass index (BMI), slower growth in fat mass, and higher odds of increasing Tanner stage of testes development between childhood and adolescence . The same study also found that AA at age 7 was associated with increased height in childhood and adolescence, but slower growth in height between childhood and adolescence , suggesting earlier physiological maturation. In a study of 46 US adolescent girls using the Horvath epigenetic age predictor, AA at age 13 years was associated with higher salivary cortisol .
There are few data on associations between epigenetic age and aging-related biomarkers early in life. There is an absence of literature of other early life phenotypes well-known to be related to aging-related diseases and/or premature mortality, namely, psychiatric problems and cognitive functioning [15, 16]. We examined 11.0–13.2-year-old Finnish adolescents to determine whether the Horvath epigenetic age predictor is related to the tempo of markers of physical growth and development, hypothalamic-pituitary-adrenal (HPA) axis functioning, psychiatric problems, and cognition. Based on preclinical studies showing that fetal exposure to excess glucocorticoids ‘programs’ an offspring phenotype of accelerated risk factors for cardiometabolic and psychiatric disorders, likely mediated, in part, via altered DNAm profiles [17, 18], we hypothesized that AA would be associated with more advanced physical growth and development, higher diurnal and lesser suppressed salivary cortisol in response to dexamethasone, higher total and internalizing and externalizing psychiatric problems, and lower scores on neuropsychological tests of intelligence.
Characteristics of the sample are in Table 1. Pearson correlations between DNAm age and chronological age was 0.13 (p = 0.041) (Fig. 1). Pubertal (Pearson’s r |0.28–0.91|), neuroendocrine (|0.41–0.58| p < 0.001; salivary cortisol awakening response and nadir r = − 0.02,p = 0.74), psychiatric (0.14–0.71 p < 0.05; DSM-IV-oriented somatic and oppositional defiant problems r = 0.13,p = 0.06), and cognitive (0.48–0.85 p < 0.001) outcomes were correlated. Of the covariates, none was significantly associated with the adolescent AA (Additional file 1: Table S1).
AA and physical growth and development
Table 2 shows that in models adjusted for adolescent’s sex and the first three MDS components (model 1), each year increase in AA was associated with 0.06 SD unit higher weight-for-age, 0.08 SD unit taller height-for-age, and − 0.09 SD units less missed from the target adult height (p values < 0.02). Also each year increase in AA was associated with a more advanced Tanner stage of breast/genitals development and a more advanced pubertal stage on the PDS (13 and 16% odds to increase in stage per each year increase in AA, respectively; p values < 0.018). When adjusted for the other covariates (model 2) the associations remained significant (p values < 0.014) except for weight-for-age SD score, which became non-significant (p = 0.051) (Table 2). When we made adjustments for maternal self-reported age at menarche, all the significant associations remained significant (p < 0.042; data not shown). When corrected for multiple testing (tests across 7 outcomes), the association with height-for-age SD score (Bonferroni-corrected p = 0.035) and own current height-for-age SD score minus mid-parental target height SD score (Bonferroni-corrected p = 0.007) remained significant.
None of the physical growth associations varied by sex (p values > 0.052 for sex × AA interactions; data not shown).
AA and diurnal and dexamethasone suppressed salivary cortisol
In models adjusting for adolescent’s sex, the first three MDS components and time at awakening, for each year increase in AA, salivary cortisol at awakening increased by 4.2% (95% CI 0.6; 7.9, p = 0.021) (Fig. 2). This association survived covariate adjustments (model 2), and adjustments for the adolescents BMI-for-age SD score (Fig. 2; p values < 0.02), but not correction for multiple testing (tests across 4 cortisol outcomes; Bonferroni-corrected p = 0.08). Higher AA also associated with salivary cortisol awakening response, but this association was sex-specific (p = 0.02 for sex × AA interaction). Additional file 2: Figure S1 shows that in boys, for each year increase in AA, salivary cortisol awakening response decreased by − 4.9% (95% CI − 10.4; 0.09, p = 0.06), while in girls it increased by 3.1% (95% CI − 1.0; 7.3, p = 0.11); in neither of the groups was the change in salivary cortisol awakening response significant. There were no other significant associations with the salivary cortisol parameters (p values > 0.06) and no other sex interactions (p values> 0.06) (data not shown).
AA and psychiatric problems
Figure 3 shows that after the model 1 covariate adjustments, each year increase in AA was associated with 29% higher odds for internalizing problems, and 29 and 34% higher odds for anxious/depressed and withdrawn problems on the internalizing problems domain, respectively (p values< 0.008); each year, increase in AA was associated with 27 and 25% higher odds for DSM-IV oriented affective and anxiety problems (p values< 0.045); and each year, increase in AA was also associated with 18% higher odds for thought problems within the domain of other problems (p = 0.035). When adjusted further for the other covariates (model 2), the association of AA with anxiety problems became non-significant (p = 0.07), while the other significant associations remained unaffected (p values < 0.034) (Fig. 3). When corrected for multiple testing (tests across 4 internalizing, 6 DSM-IV, and 3 other problems domains), the associations with internalizing, anxious/depressed, withdrawn, and affective problems (Bonferroni-corrected p value< 0.018) remained significant.
There were no other significant associations with child psychiatric problems (p values> 0.09; data not shown), and no significant sex differences in these associations (p values> 0.06 for sex × AA interactions; data not shown).
AA and cognition
There were no significant associations between AA and estimated intelligence (Additional file 1: Table S2).
As we hypothesized, adolescents with higher AA, meaning higher DNAm age than chronological age, was associated with more advanced physical growth and development, higher salivary cortisol, and higher odds for psychiatric problems. Adolescents with higher AA were heavier- and taller-for-age, and closer to their expected target adult height, suggesting an earlier growth spurt and less remaining growth potential. Their pubertal stage of breast/genital development, according to the Tanner Staging Questionnaire and of secondary sex characteristics according to PDS, were also at a more advanced stage.
The advanced growth and maturation was accompanied by higher salivary cortisol upon awakening, and they had higher odds for displaying borderline clinically significant internalizing problems, in particular anxious/depressed and withdrawn problems, and affective and anxiety as well as thought problems. While these associations remained significant after adjustments for a number of important covariates, including maternal pregnancy and child perinatal characteristics and genetic population structure, only the associations with stature and with internalizing/affective problems remained significant after correction for multiple testing. Yet, the Bonferroni-correction may be an overly stringent method to account for the multiple testing problem here, as the outcomes, within each developmental domain were not independent of each other. These results thus suggest that adolescents whose biological DNAm age is more advanced than their chronological age display physiological age acceleration that may indicate risk of earlier aging.
Our findings are consistent with the life history theory that suggests that early development and early puberty are meaningful tradeoffs in conditions of environmental adversity [19, 20]. We suggest that AA and advanced growth and maturation are indicators of more advanced tempo of aging processes present from early life onwards. Indeed, more advanced physical growth and pubertal development have been shown to predict aging-related diseases, including cancers , cardio-metabolic disorders and their risk factors , and depression . While associations between cortisol, the effector hormone of the HPA axis, and health are complex, in general high cortisol concentrations, are associated with a number of physical and mental health adversities, such as obesity , sleep problems , anxiety, and depression  in studies in both children and adults. Moreover, longitudinal studies suggest that high cortisol levels in late adulthood increase the risk for cardiometabolic disorders and cardiovascular mortality [27, 28]. While childhood psychiatric problems are associated with a number of physical health adversities and psychosocial problems , there is evidence that they tend to track into adulthood . Even in those children whom psychiatric problems do not persist into adulthood or are subthreshold , and in those whom childhood psychiatric problems are parent-rated , are at increased risk for adverse adulthood outcomes, and hence, are likely more vulnerable for earlier aging.
Our findings with physical growth agree with one previous study  which demonstrated that higher AA at age 7 was associated with increased height in childhood and adolescence. Yet, the study found that higher AA at age 7 predicted slower, rather than faster height growth between childhood and adolescence, and that AA at ages 7 and 17 years was not associated with a number of markers of pubertal development, including peak height velocity and the Tanner Staging Questionnaire . The previous study also found that higher AA at birth predicted higher fat mass, faster growth in weight and BMI, slower growth in fat mass, and higher odds of increasing Tanner stage of testes development between childhood and adolescence . As DNAm undergoes age-related changes , a potential explanation for the somewhat discrepant study findings on physical growth and development between our and the previous study is the age-stage at which DNAm and physical growth were measured. The difference in findings does not, however, relate to tissue type as both studies measured DNAm in childhood/adolescence from venous blood. Yet, the p values in the previous study were at their best around 0.007, and none of the associations would have survived Bonferroni-correction for multiple testing.
Our findings for salivary cortisol are also in partial agreement with one previous study in adolescent girls, which demonstrated that higher AA, measured from salivary DNA, was associated with higher salivary cortisol measured across 2 days. This study did not, however, account for covariates, except for adolescent’s age at testing, and neither of the previous studies accounted for genetic population structure, which is strongly associated with methylation profiles . Hence, it remains unclear if the findings of the latter would have survived adjustments for important covariates, such as BMI, and if both of these previous studies would have survived adjustments for genetic population structure.
Our study also revealed novel findings related to psychiatric problems, namely, that higher AA was associated with internalizing/affective-type and thought problems, but not externalizing-type of problems. This finding may reflect statistical power, as internalizing problems in our sample were twice as prevalent as externalizing problems, or studying adolescents, a time when a marked rise in internalizing problems is observed . This pattern is, however, congruent with our other recent study. Interestingly, we have demonstrated that another epigenetic age biomarker, namely, lower DNAm gestational age than chronological gestational at birth predicts higher internalizing, but not externalizing problems in early childhood, though in a sex-specific manner such that boys are more vulnerable . AA was not associated with intelligence. This is somewhat surprising, as poorer childhood cognitive functioning is predictive of aging-related diseases, including dementia . Hence, because of the developmental changes in DNAm , we cannot rule out that AA at later developmental stages will change and become associated with cognitive function and perhaps change its links to behavioral problems.
The strengths of our study relate to a well-characterized cohort and availability of a number of aging-related phenotypes that we measured decades before the aging-related diseases become manifest. We were also able to account for a number of early life adversities and their proxies, such as maternal smoking and alcohol use during pregnancy, maternal age and BMI at delivery, mode of delivery, and the adolescent’s birth weight and gestational age. We also accounted for maternal glycyrrhizin in licorice use during pregnancy, which is a potent inhibitor of the placental glucocorticoid barrier enzyme (11-beta hydroxysteroid dehydrogenase type 2) and which may result in fetal overexposure to maternal circulating glucocorticoids; our study was originally designed to examine its associations with offspring developmental outcomes, and it is commonly consumed in young Finnish women (in our original cohort nearly 50%) . None of these early life factors were significantly associated with epigenetic AA in this sample, which is contrary to what we expected based on previous findings by us and others showing that these factors may exert adverse consequences on offspring neurodevelopment and HPA-axis functioning [37, 38]. With regard to maternal glycyrrhizin use, we have previously shown in a larger sample of this study cohort that high maternal intake of glycyrrhizin in licorice during pregnancy (> 500 mg/week) is associated with slightly shorter length of gestation , poorer performance in neurocognitive tests at ages 8 and 11–13 [40, 41], higher odds for having borderline clinically significant externalizing psychiatric problems at ages 8 and 11–13 [40, 41], higher diurnal and stress-induced salivary cortisol profiles at age 8 , and more advanced pubertal maturation in girls at age 11–13 . Whether maternal intake of glycyrrhizin in licorice during pregnancy is associated with other epigenetic biomarkers than the epigenetic age biomarker studied here, awaits further investigation.
The limitations of our study are the narrow age range of our sample, and hence, the small magnitude of the correlation between DNAm age and chronological age. The small magnitude of this correlation is, however, similar to the other two previous childhood epigenetic age studies [13, 14]. The narrow age range also limits generalizability from our findings to samples that differ in age from ours as DNAm changes with age. Our findings are also limited to DNAm in one tissue type and our study precludes generalizations beyond Finnish children. Further, our study design was cross-sectional which precludes testing developmental changes and causal inferences. We can neither address the possibility of selection bias resulting from sample attrition to genetic analyses. Also measuring genome-wide methylation with the most recent Illumina EPIC array that lacks 16 (4.5%) of the 353 CpG sites originally needed for Horvath DNAm age calculation which should be kept in mind when interpreting the study findings. Finally, while we measured psychiatric problems with a standardized, validated, and widely used mother-report , we cannot rule out potential information-bias embedded in the mother-report. Hence, future studies need to confirm whether our study findings on psychiatric problems also pertain to clinical diagnoses.
Our study shows that among 11.0–13.2-year-old adolescents, AA is associated with a number of markers that index risk for earlier aging, namely, more advanced physical growth and development, higher salivary cortisol upon awakening, and psychiatric problems. Our findings are consistent with the life history theory [19, 20] and lend credence to the proposition that AA may be used as a biomarker of aging already early in life [3, 13].
Glycyrrhizin in Licorice (GLAKU) is an urban community-based cohort comprising originally 1049 women and their healthy, singleton infants born in 1998 in Helsinki, Finland . Between 2009 and 2011, all initial cohort members who had given permission to be contacted, and whose addresses were traceable (N = 920, 87.70% of the original cohort in 1998) were invited to a follow-up. Of them, 692 (75.20%) could be contacted by phone, and 451 (65.20% of those contacted by phone) participated in a follow-up at the child’s mean age of 12.3 years (SD = 0.5, range 11.0–13.2 years). Of the participating adolescents, 243 donated blood for genetic analyses. After quality control procedures, 239 DNA samples remained for genetic analyses.
In comparison to those who participated in the follow-up but who did not donate blood or were excluded from genetic analyses based on quality control reasons (n = 212), the analytic sample with DNA (n = 239) had higher weight (mean difference (MD) = 3.27 kg, p < 0.01), height (MD = 2.72 cm, p < 0.01) and BMI (MD = 0.69 kg/m2, p = 0.03) at the adolescence follow-up, and their mothers had higher weight at delivery (MD = 1.76 kg, p = 0.05). The groups did not differ in chronological age at adolescent follow-up, sex, birth order, body size at birth, length of gestation, nor did their mothers differ in licorice or alcohol consumption or smoking during pregnancy, BMI at delivery or mode of delivery, age at menarche; their mothers and fathers did not differ in educational attainment or height (all p values> 0.07).
Ethics Committees of the City of Helsinki and the Uusimaa Hospital District approved the study protocol. Written informed consent was obtained from the mother at birth and from parent/guardian and adolescent at the follow-up.
Adolescent DNA methylation, epigenetic age, and blood cell count composition
Blood samples were collected, and DNA was extracted according to standard procedures. Methylation analyses were performed at the Max Planck Institute of Psychiatry in Munich, Germany. DNA was bisulphite-converted using the EZ-96 DNA Methylation kit (Zymo Research, Irvine, CA). Genome-wide methylation status of over 850,000 CpG sites was measured using the Illumina Infinium MethylationEPIC arrays (Illumina Inc., San Diego, CA) according to the manufacturer’s protocol. The arrays were scanned using the iScan System (Illumina Inc., San Diego, CA). The final dataset contained 812,943 CpGs.
We obtained DNAm-predicted age based on Horvath method  using the online epigenetic clock calculator (http://labs.genetics.ucla.edu/horvath/dnamage/). This calculator also incorporates information on blood cell counts for 6 cell types (granulocytes, monocytes, natural killer cells, B cells, CD4+ T cells, and CD8+ T cells) based on the Houseman method . Epigenetic age was calculated as the unstandardized residual from a linear regression of DNAm age on chronological age and 6 cell count types.
Adolescent genotyping and multi-dimensional scaling analysis
Genotyping was performed on Illumina Human OmniExpress Exome 1.2 bead chip (Illumina Inc., San Diego, CA) at the Tartu University, Estonia in September 2014 according to the standard protocols. Genomic coverage was extended by imputation using the 1000 Genomes Phase I integrated variant set (v3/April 2012; NCBI build 37/hg19) as the reference sample and IMPUTE2 software. Before imputing the following QC, filters were applied: SNP clustering probability for each genotype > 95%, Call rate > 95% individuals and markers (99% for markers with MAF < 5%), MAF > 1%, HWE p > 1*10–6. Moreover, heterozygosity, sex check, and relatedness checks were performed and any discrepancies were removed (N = 2). We performed multi-dimensional scaling (MDS) analysis on the identity by state matrix of quality-controlled genotypes. The first three components depicted the origin admixture and were included as covariates in the regression analyses . This information was available for 221 participants.
Physical growth and development
We used the following three measures of growth and development: (i) the difference between the child’s height-for-age standard-deviation (SD) score based on Finnish growth charts  (using the current measured height without shoes, measured with a Seca stadiometer) (model 213; Seca GmbH & Co KG, Hamburg, Germany) and midparental target height in SD units ; this is a measure of remaining growth potential and, consequently, the timing of the pubertal growth spurt. (ii) The Tanner Staging Questionnaire administered by a research nurse. The questionnaire uses schematic drawings of two secondary sex characteristics (pubic hair development in girls and boys and breast development in girls and development of genitalia in boys) and yields two 5-stage scores ranging from pre-pubertal (stage I) to post-pubertal (stage V) . (iii) The Pubertal Development Scale (PDS) is a self-report questionnaire on secondary sex characteristics (growth spurt, body hair that is not specifically pubic hair, and skin changes in girls and boys; menarche and breast development in girls; and facial hair and voice change in boys). The PDS yields one 4-stage score ranging from no development (I) to full completion of development (IV) .
We also measured weight in light clothing without shoes (model 8; Seca GmbH & Co KG) and calculated BMI (weight (kg)/height (m2)). We transformed the values into weight-for-age and BMI-for-age SD scores based on Finnish growth charts .
Diurnal and dexamethasone-suppressed salivary cortisol
Saliva samples were collected using cotton swabs (Salivette; Sarstedt, Nümbrecht, Germany). The adolescents were asked not to eat, drink, or brush their teeth 30 min before each sample, with an exception of drinking a small amount of water. We also specifically asked them to restrain from any caffeine including drinks (coffee, tea, cola). On the first of two consecutive days, samples were collected upon awakening and 15, 30, 45, and 60 min thereafter, at 12:00 noon, at 5:00 p.m., and at bedtime. Dexamethasone was administered after the bedtime saliva sample, and a sample was collected upon awakening the next day. We used a low dose of dexamethasone (3 μg/kg of total body weight) to detect individual variation in hypothalamic-pituitary-adrenal (HPA) axis suppression . Salivary cortisol concentrations were determined by solid-phase, time-resolved fluorescence immunoassay with fluorometric end-point detection (DELFIA; Wallac, Turku, Finland). The intraassay and interassay coefficients of variation varied between 4.0 and 7.7%, and the mean coefficient of variation between duplicate analyses was 5.9%.
Of the diurnal measures, we used/calculated the following parameters: cortisol at awakening, cortisol awakening response (peak value after awakening minus value upon awakening), nadir (minimum of diurnal values), and response to dexamethasone suppression test (value upon awakening on day two minus value upon awakening on day one).
Mothers completed the Child Behavior Checklist (CBCL/6-18), a standardized and validated rating scale screening for psychiatric problems . The scale is hierarchically structured, such that it yields a total problems score, which is first subdivided into internalizing and externalizing problems scores, and then internalizing problems score is further subdivided into anxious/depressed, withdrawn, and somatic complaints scores, and externalizing problems into rule-breaking and aggressive behavior scores; scores not included in internalizing and externalizing problems score, but which are embedded in the total problems are social, thought, and attention problems . The CBCL also yields six Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV)-oriented scores, namely, affective, anxiety, somatic, attention deficit hyperactivity, conduct, and oppositional-defiant problems scores . Following the CBCL manual, we used the 82nd percentile as the cutoff to identify adolescents with borderline clinically significant problems .
We used the short form of the Wechsler-Intelligence-Scale-for-Children-III  which included vocabulary, similarities, block design, and picture arrangement subtests. We used age-standardized scores to estimate age-standardized total intelligence, and verbal and performance intelligence quotients (IQs) .
Covariates and confounders
All analyses were adjusted for child’s sex and the first three MDS components to control for population structure (model 1). We made further adjustments for covariates previously associated with physical growth and development, salivary cortisol, psychiatric problems, and cognition in this cohort : the highest educational level of either parent (secondary or less/vocational/university) reported at adolescent follow-up, maternal age (years), and BMI (kg/m2) at delivery calculated from weight and height derived from medical records, maternal smoking (no/yes), weekly alcohol (no/yes) and glycyrrhizin in licorice (0–249 mg/week, 250–499 mg/week, ≥ 500 mg/week) consumption during pregnancy, delivery mode (vaginal/cesarean), parity (primiparous/multiparous), gestational length (weeks) as confirmed by ultrasound scans, and birth weight (grams) of the adolescent, derived from birth records (model 2). In addition, we conducted analyses of pubertal maturation adjusting for maternal self-reported age at menarche (years) as a crude proxy of the genetic component of pubertal development. Analyses of HPA-axis activity were additionally adjusted for time at awakening and time at dexamethasone intake as well as for child’s BMI-for-age SD score.
We use generalized linear models (GLM) to study associations between AA and outcomes, specifying Gaussian reference distribution for continuous (growth anthropometry, salivary cortisol, cognition), ordinal logistic for categorical (Tanner stages and PDS), and binary logistic reference distribution for dichotomous outcomes (psychiatric problems). All analyses were adjusted for covariates and confounders as described and we also tested if the associations between AA on outcomes varied by sex by including sex × AA interaction into the GLMs following main effects of these variables as the sexes may differ in epigenetic age, pubertal maturation, and the prevalence and etiology of psychiatric problems. We also report Bonferroni-corrected p values to account for multiple testing within each developmental domain. Statistical analyses were performed using IBM SPSS version 24.0.
- 95% CI:
95% confidence interval
Epigenetic age acceleration
Glycyrrhizin in Licorice Study
United Nations, Department of Economic and Social Affairs, Population Division. World Population Ageing 2015 (ST/ESA/SER.A/390).
Suzman R, et al. Health in an ageing world—what do we know? Lancet. 2015;385(9967):484–6.
Jones MJ, Goodman SJ, Kobor MS. DNA methylation and healthy human aging. Aging Cell. 2015;14(6):924–32.
Hannum G, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49(2):359–67.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14(10):R115.
Horvath S, et al. Obesity accelerates epigenetic aging of human liver. Proc Natl Acad Sci U S A. 2014;111(43):15538–43.
Marioni RE, et al. The epigenetic clock is correlated with physical and cognitive fitness in the Lothian birth cohort 1936. Int J Epidemiol. 2015;44(4):1388–96.
Levine ME, et al. Epigenetic age of the pre-frontal cortex is associated with neuritic plaques, amyloid load, and Alzheimer’s disease related cognitive functioning. Aging. 2015;7(12):1198–211.
Levine ME, et al. Menopause accelerates biological aging. Proc Natl Acad Sci. 2016;113(33):9327–32.
Marioni RE, et al. DNA methylation age of blood predicts all-cause mortality in later life. Genome Biol. 2015;16:25.
Chen BH, et al. DNA methylation-based measures of biological age: meta-analysis predicting time to death. Aging (Albany NY). 2016;8(9):1844–65.
Belsky DW, et al. Quantification of biological aging in young adults. Proc Natl Acad Sci USA. 2015;112(30):E4104–10.
Simpkin AJ, et al. The epigenetic clock and physical development during childhood and adolescence : longitudinal analysis from a UK birth cohort. Int J Epidemiol. 2017;46(2):549–58.
Davis EG, et al. Accelerated DNA methylation age in adolescent girls associations with elevated diurnal cortisol and reduced hippocampal volume. Transl Psych. 2017;7(March):10–4.
von Stumm S, et al. Childhood behavior problems and health at midlife: 35-year follow-up of a Scottish birth cohort. J Child Psychol Psychiatry. 2011;52(9):992–1001.
Calvin CM, et al. Childhood intelligence in relation to major causes of death in 68 year follow-up: prospective population study. BMJ. 2017;357:j2708.
Seckl JR, Meaney MJ. Glucocorticoid programming. Ann N Y Acad Sci. 2004;1032:63–84.
Reynolds RM, et al. Transmitting biological effects of stress in utero: implications for mother and offspring. Psychoneuroendocrinology. 2013;38(9):1843–9.
Stearns S. The evolution of life histories. New York: Oxford University Press; 1992.
Flatt T, Heyland A. Mechanisms of life history evolution: the genetics and physiology of life history traits and trade-offs. Oxford: OUP Oxford; 2011. p. 2011.
Day FR, et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat Genet. 2017;49(6):834–41.
Widen E, et al. Pubertal timing and growth influences cardiometabolic risk factors in adult males and females. Diabetes Care. 2012;35(4):850–6.
Graber JA. Pubertal timing and the development of psychopathology in adolescence and beyond. Horm Behav. 2013;64(2):262–9.
Incollingo Rodriguez AC, et al. Hypothalamic-pituitary-adrenal axis dysregulation and cortisol activity in obesity: a systematic review. Psychoneuroendocrinology. 2015;62:301–18.
Pesonen AK, et al. Sex-specific associations between sleep problems and hypothalamic-pituitary-adrenocortical axis activity in children. Psychoneuroendocrinology. 2012;37(2):238–48.
Zorn JV, et al. Cortisol stress reactivity across psychiatric disorders: a systematic review and meta-analysis. Psychoneuroendocrinology. 2017;77:25–36.
Vogelzangs N, et al. Urinary cortisol and six-year risk of all-cause and cardiovascular mortality. J Clin Endocrinol Metab. 2010;95(11):4959–64.
Schoorlemmer RM, et al. Relationships between cortisol level, mortality and chronic diseases in older persons. Clin Endocrinol. 2009;71(6):779–86.
National Research Council and Institute of Medicine. Preventing mental, emotional, and behavioral disorders among young people: progress and possibilities. In: O’Connell ME, Boat T, Warner KE, editors. Committee on the Prevention of Mental Disorders and Substance Abuse Among Children, Youth, and Young Adults: Research Advances and Promising Interventions. Board on Children, Youth, and Families, Division of Behavioral and Social Sciences and Education. Washington, DC: The National Academies Press; 2009.
Copeland WE, et al. Adult functional outcomes of common childhood psychiatric problems: a prospective, longitudinal study. JAMA Psychiatry. 2015;72(9):892–9.
Sourander A, et al. Recognition of psychiatric disorders, and self-perceived problems. A follow-up study from age 8 to age 18. J Child Psychol Psychiatry. 2005;46(10):1124–34.
Jung M, Pfeifer GP. Aging and DNA methylation. BMC Biol. 2015;13:7.
Barfield RT, et al. Accounting for population stratification in DNA methylation studies. Genet Epidemiol. 2014;38(3):231–41.
Wichstrom L. The emergence of gender difference in depressed mood during adolescence: the role of intensified gender socialization. Dev Psychol. 1999;35(1):232–45.
Suarez A, et al. The epigenetic clock at birth: associations with maternal antenatal depression and child psychiatric problems. J Am Acad Child Adolesc Psychiatry. 2018;57(5):321–328.e2.
Strandberg TE, et al. Birth outcome in relation to licorice consumption during pregnancy. Am J Epidemiol. 2001;153(11):1085–8.
Girchenko P, et al. Maternal early pregnancy obesity and related pregnancy and pre-pregnancy disorders: associations with child developmental milestones in the prospective PREDO study. Int J Obes. 2018;42:995–1007.
Kajantie E, et al. Body size at birth predicts hypothalamic-pituitary-adrenal axis response to psychosocial stress at age 60 to 70 years. J Clin Endocrinol Metab. 2007;92(11):4094–100.
Strandberg TE, et al. Preterm birth and licorice consumption during pregnancy. Am J Epidemiol. 2002;156(9):803–5.
Raikkonen K, et al. Maternal licorice consumption and detrimental cognitive and psychiatric outcomes in children. Am J Epidemiol. 2009;170(9):1137–46.
Raikkonen K, et al. Maternal licorice consumption during pregnancy and pubertal, cognitive, and psychiatric outcomes in children. Am J Epidemiol. 2017;185(5):317–28.
Raikkonen K, et al. Maternal prenatal licorice consumption alters hypothalamic-pituitary-adrenocortical axis function in children. Psychoneuroendocrinology. 2010;35(10):1587–93.
Achenbach, T.M. and L.A Rescorla, Manual for the ASEBA school-age forms and profiles. English, 2001: p. University of Vermont, Research Center for Childre-University of Vermont, Research Center for Childre.
Houseman EA, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012;13(1):86.
Saari A, et al. New Finnish growth references for children and adolescents aged 0 to 20 years: length/height-for-age, weight-for-length/height, and body mass index-for-age. Ann Med. 2011;43(3):235–48.
Sorva R, et al. Growth evaluation: parent and child specific height standards. Arch Dis Child. 1989;64(10):1483–7.
Morris NM, Udry JR. Validation of a self-administered instrument to assess stage of adolescent development. J Youth Adolescence. 1980;9(3):271–80.
Petersen AC, et al. A self-report measure of pubertal status: reliability, validity, and initial norms. J Youth Adolescence. 1988;17(2):117–33.
Martikainen S, et al. Physical activity and hypothalamic-pituitary-adrenocortical axis function in adolescents. Psychoneuroendocrinology. 2014;49(1):96–105.
Wechsler D. Wechsler intelligence scale for children–Third edition. San Antonio: The Psychological Corporation; 1991.
Kaufman AS, et al. Comparison of three WISC-III short forms: weighing psychometric, clinical, and practical factors. J Clin Child Psychol. 1996;25(1):97–105.
We thank the GLAKU study participants for their enthusiastic participation.
The GLAKU study is funded by the Academy of Finland, EraNetNeuron, EVO, University of Helsinki Research Funds, the Signe and Ane Gyllenberg, Emil Aaltonen, Novo Nordisk, Päivikki and Sakari Sohlberg, and Sigrid Juselius, the Finnish Medical and the British Heart Foundation, the European Commission Horizon 2020 Award SC1-2016-RTD-733280 RECAP and the Doctoral Programme of Psychology, Learning and Communication.
Availability of data and materials
Any interested researchers can obtain a de-identified dataset after having obtained an approval from the GLAKU Study Board. Data requests may be subject to further review by the national register authority and ethical committees. Any requests for data use should be addressed to GLAKU Study individual researchers.
Ethics approval and consent to participate
Ethics Committees of the City of Helsinki and the Uusimaa Hospital District approved the study protocol. Written informed consent was obtained from the mother at birth and from parent/guardian and adolescent at the follow-up.
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.
Table S1. Associations between epigenetic age acceleration and covariates in 11.0–13.2-year-old adolescents. Table S2. Associations between epigenetic age acceleration and cognition in 11.0–13.2-year-old adolescents. (DOCX 17 kb)
Figure S1. A scatterplot with a regression lines showing associations between epigenetic age acceleration and salivary cortisol awakening response in 11.0–13.2-year-old adolescent boys and girls. Epigenetic age acceleration is calculated as the residual from a linear regression where DNA methylation age is regressed on chronological age and adjusted for 6 cell types. Numbers showing percent increase in salivary cortisol upon awakening per 1 year increase in epigenetic age acceleration and 95% confidence intervals are derived from generalized linear models with Gaussian reference distribution and adjusted for three multidimensional scaling components from genome-wide data and time upon awakening. (PPTX 38 kb)
About this article
Cite this article
Suarez, A., Lahti, J., Czamara, D. et al. The epigenetic clock and pubertal, neuroendocrine, psychiatric, and cognitive outcomes in adolescents. Clin Epigenet 10, 96 (2018). https://doi.org/10.1186/s13148-018-0528-6
- Aging-related Diseases
- Pubertal Development Scale (PDS)
- Tanner Stage
- Salivary Cortisol Awakening Response
- Childhood Psychiatric Problems