# Associations of four biological age markers with child development: A multi-omic analysis in the European HELIX cohort

Oliver Robinson, Chungho E Lau, Sungyeon Joo, Sandra Andrusaityte, Eva Borras, Paula de Prado-Bert, Lida Chatzi, Hector C Keun, Regina Grazuleviciene, Kristine B Gutzkow, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Oliver Robinson, Chungho E Lau, Sungyeon Joo, Sandra Andrusaityte, Eva Borras, et al. Associa-
tions of four biological age markers with child development: A multi-omic analysis in the European
HELIX cohort. eLife, 2023, 12, pp.108036. 10.7554/elife.85104. hal-04302446

HAL Id: hal-04302446

https://hal.science/hal-04302446v1

Submitted on 23 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

# Associations of four biological age markers with child development: A multi-omic analysis in the European HELIX cohort

Oliver Robinson¹,²*, ChungHo E Lau¹, Sungyeon Joo¹, Sandra Andrusaityte³, Eva Borras⁴,⁵, Paula de Prado-Bert⁵,⁶,⁷, Lida Chatzi⁸, Hector C Keun⁹,¹⁰, Regina Grazuleviciene³, Kristine B Gutzkow¹¹, Lea Maitre⁵,⁶,⁷, Dries S Martens¹², Eduard Sabido⁴,⁵, Valérie Siroux¹³, Jose Urquiza⁵,⁶,⁷, Marina Vafeiadi¹⁴, John Wright¹⁵, Tim S Nawrot¹², Mariona Bustamante⁵,⁶,⁷, Martine Vrijheid⁵,⁶,⁷

¹Medical Research Council Centre for Environment and Health, Imperial College London, London, United Kingdom; ²Mohn Centre for Children's Health and Well-being, School of Public Health, Imperial College London, London, United Kingdom; ³Department of Environmental Science, Vytautas Magnus University, Kaunas, Lithuania; ⁴Center for Genomics Regulation, Barcelona Institute of Science and Technology, Barcelona, Spain; ⁵Universitat Pompeu Fabra (UPF), Barcelona, Spain; ⁶Institute for Global Health (ISGlobal), Barcelona, Spain; ⁷CIBER Epidemiologa y Salud Pública (CIBERESP), Madrid, Spain; ⁸Department of Preventive Medicine, Keck School of Medicine, University of Southern California, Los Angeles, United States; ⁹Division of Systems Medicine, Department of Metabolism, Digestion and Reproduction, Imperial College London, London, United Kingdom; ¹⁰Cancer Metabolism & Systems Toxicology Group, Division of Cancer, Department of Surgery & Cancer, Imperial College London, London, United Kingdom; ¹¹Division of Climate and Environmental Health, Norwegian Institute of Public Health, Oslo, Norway; ¹²Centre for Environmental Sciences, Hasselt University, Hasselt, Belgium; ¹³University Grenoble Alpes, Inserm U 1209, CNRS UMR 5309, Team of environmental epidemiology applied to the development and respiratory health, Institute for Advanced Biosciences, 38000 Grenoble, France; ¹⁴Department of Social Medicine, School of Medicine, University of Crete, Heraklion, Crete, Greece; ¹⁵Bradford Institute for Health Research, Bradford Teaching Hospitals NHS Foundation Trust, Bradford, United Kingdom

*For correspondence: o.robinson@imperial.ac.uk

**Competing interest:** The authors declare that no competing interests exist.

**Funding:** See page 24

**Received:** 22 November 2022

**Preprinted:** 23 January 2023

**Accepted:** 06 June 2023

**Published:** 06 June 2023

**Reviewing Editor:** Sara Hägg, Karolinska Institutet, Sweden

© Copyright Robinson et al. This article is distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use and redistribution provided that the original author and source are credited.

## Abstract

**Background:** While biological age in adults is often understood as representing general health and resilience, the conceptual interpretation of accelerated biological age in children and its relationship to development remains unclear. We aimed to clarify the relationship of accelerated biological age, assessed through two established biological age indicators, telomere length and DNA methylation age, and two novel candidate biological age indicators, to child developmental outcomes, including growth and adiposity, cognition, behavior, lung function and the onset of puberty, among European school-age children participating in the HELIX exposome cohort.

**Methods:** The study population included up to 1173 children, aged between 5 and 12 years, from study centres in the UK, France, Spain, Norway, Lithuania, and Greece. Telomere length
was measured through qPCR, blood DNA methylation, and gene expression was measured using microarray, and proteins and metabolites were measured by a range of targeted assays. DNA methylation age was assessed using Horvath's skin and blood clock, while novel blood transcriptome and 'immunometabolic' (based on plasma proteins and urinary and serum metabolites) clocks were derived and tested in a subset of children assessed six months after the main follow-up visit. Associations between biological age indicators with child developmental measures as well as health risk factors were estimated using linear regression, adjusted for chronological age, sex, ethnicity, and study centre. The clock derived markers were expressed as Δ age (i.e. predicted minus chronological age).

**Results:** Transcriptome and immunometabolic clocks predicted chronological age well in the test set (r=0.93 and *r*=0.84 respectively). Generally, weak correlations were observed, after adjustment for chronological age, between the biological age indicators.

Among associations with health risk factors, higher birthweight was associated with greater immunometabolic Δ age, smoke exposure with greater DNA methylation Δ age, and high family affluence with longer telomere length.

Among associations with child developmental measures, all biological age markers were associ-
ated with greater BMI and fat mass, and all markers except telomere length were associated with
greater height, at least at nominal significance (p<0.05). Immunometabolic Δ age was associated
with better working memory (p=4 e-3) and reduced inattentiveness (p=4 e-4), while DNA methyla-
tion Δ age was associated with greater inattentiveness (p=0.03) and poorer externalizing behaviors
(p=0.01). Shorter telomere length was also associated with poorer externalizing behaviors (p=0.03).

**Conclusions:** In children, as in adults, biological aging appears to be a multi-faceted process and adiposity is an important correlate of accelerated biological aging. Patterns of associations suggested that accelerated immunometabolic age may be beneficial for some aspects of child development while accelerated DNA methylation age and telomere attrition may reflect early detrimental aspects of biological aging, apparent even in children.

**Funding:** UK Research and Innovation (MR/S03532X/1); European Commission (grant agreement numbers: 308333; 874583).

Editor's evaluation

This is an important study that examined multiple biological age measures in children, which has
been lacking in literature. The findings of this study provided convincing evidence to interpret and
understand the aging and developmental processes in children.

**Introduction**

The field of geroscience proposes that biological aging, a set of interrelated molecular and cellular
changes associated with aging, drives the physiological deterioration that is the root of multiple age-
related health conditions (Sierra, 2016). Understanding the process of biological aging and devel-
oping markers to accurately assess biological age in individuals, holds great promise for public health
and biomedical research in general to develop interventions, even in childhood and early life, that
slow physiological decline and reduce the risk of chronic disease and disability in later life.

Telomere length, which shortens with age, is one of the most widely applied biological age markers primarily as it directly assesses a primary Hallmark of Aging (López-Otín et al., 2013; Jylhävä et al., 2017). More recently, high-throughput 'omic' methods, which provide simultaneous quantification of thousands of epigenetic marks, transcripts, proteins, and metabolites, have been used to develop 'biological clocks' that provide a global measure of changes with age at the molecular level (Rutledge et al., 2022). While biological clocks have been primarily trained on chronological age, 'age acceleration,' commonly defined as the difference between clock-predicted age and chronological age, has been associated with age-related phenotypes and mortality (Hertel et al., 2016; Perna et al., 2016; Marioni et al., 2015; Peters et al., 2015; Lehallier et al., 2019; van den Akker et al., 2020; Macdonald-Dunlop et al., 2022), indicating their utility as biological age markers. DNA methylation-based clocks, such as the clock of Horvath, 2013, have been extensively applied in large-scale studies and remain a research field under active development, with 'second generation' clocks further
**eLife digest** Although age is generally measured by the number of years since birth, many factors contribute to the rate at which a person physically ages. In adults, linking these measurements to age gives a measure of overall health and resilience. This 'biological age' offers a better prediction of remaining life and disease risk than the number of years lived.

Multiple factors can be used to calculate biological age, such as measuring the length of telomeres – protective caps on the end of chromosomes – which shorten as people age. The rate at which they shorten can give an indication of how quickly someone is ageing. Researchers can also study epigenetic factors: these mechanisms lead to certain genes being switched on or off, and they can be combined into a 'epigenetic clock' to assess biological age. However, compared with adults, the relationship between biological age and child health and developmental maturity is less well understood.

Robinson et al. studied 1,173 school-aged children from six European countries, measuring telomere length, epigenetic factors and other biological indicators related to metabolism and the immune system. The relationships between these factors and an array of child developmental measures such as height, weight, behaviour and the age of onset of puberty were established. The findings showed that biological age indicators are only weakly linked to each other in children. Despite this, biological age was related to greater amount of body fat across all tested indicators – which is also associated with biological age in adults and is an important determinant of lifespan.

Among several observed effects on development, analysis found that shorter telomere length and older epigenetic age were associated with greater behavioural problems, suggesting they may be detrimental to child development. On the other hand, a greater age due to metabolic and immune related changes was associated with greater cognitive and behavioural maturity. Environmental factors were also linked to biological ageing, with children exposed to smoking in their homes or while their mother was pregnant displaying an older epigenetic age.

Robinson et al. showed that biological ageing in children is multifaceted and can have both beneficial and harmful impacts on development. This knowledge is important for identifying early life risk factors that might influence healthy ageing in later life. Future work will help researchers to understand these complex interactions and the long-term consequences for health and well-being.

incorporating clinical biomarker and mortality information to improve their clinical utility (Lu et al., 2019; Belsky et al., 2020). Further clocks have been developed using transcriptome (Peters et al., 2015), metabolome (Robinson et al., 2020), and proteome (Lehallier et al., 2019) data, including those that specifically target immune-system-related proteins (Sayed et al., 2021). Generally, clocks have been found to be only weakly correlated with each other, suggesting that each clock captures different facets of biological aging (Belsky et al., 2018; Jansen et al., 2021).

While biological age in adults is intuitively understood as an overall indicator of general health and resilience, the conceptual interpretation of biological age acceleration in children is much less clear. Child development and aging may at first be considered opposing processes, representing growth and decay respectively. However, various related theoretical frameworks link the two processes: Under the developmental origin of health and disease hypothesis, the early life environment is a key determinant of ageing trajectories and disease risk in later life. Life-course models of ageing, supported by measures of physical and cognitive capability, view the childhood developmental phase as key to building up 'biological capital' and to determining how long capabilities and disease risk remain above critical thresholds in later years following the gradual decline phase of adult life (Kuh, 2007). Horvath's DNA methylation clock is currently the only clock trained to predict age throughout the lifespan, and many of the clock's CpGs are in genomic regions known to regulate development and differentiation (Horvath, 2013). However, unlike life-course models of physical function, the level of DNA methylation at the clock's CpGs changes in a predictable, unidirectional manner throughout the life-course, albeit at a much faster rate during childhood. This continuous molecular readout suggests that processes directing development are at least indirectly related to the detrimental processes in later-life and is consistent with quasi-programmed theories of aging such as antagonistic pleiotropy (Emery Thompson, 2022), whereby molecular functions that promote development, inadvertently lead to aging in later life (Horvath and Raj, 2018). Therefore, some authors have suggested that DNA
methylation-based age acceleration may be beneficial during childhood (Horvath and Raj, 2018; McEwen et al., 2020), reflecting greater physical maturity and build-up of biological capital.

Biological aging is conceived as a continuous balance of cellular damage, caused by both extrinsic environmental factors and by normal physiological processes, and resiliency mechanisms that protect against and compensate for this damage (Ferrucci et al., 2020). An alternative ‘wear and tear’ model would view cellular damage to occur continuously from birth and, since the epigenetic clock has been proposed to reflect the epigenomic maintenance system, a resiliency mechanism, DNA methylation age acceleration in children may, as in adults, represent a greater accumulation of epigenetic instability and, therefore, reduced biological capital. However, so far only a handful of studies have examined associations with developmental maturity in children (Binder et al., 2018; Simpkin et al., 2017; Suarez et al., 2018). Telomere length attrition is more rapid in early childhood during rapid somatic growth and more gradual in adulthood, with those with a shorter telomere length in childhood maintaining a lower telomere length into adulthood (Martens et al., 2021). While telomere length may serve as both a mitotic-clock and as a mediator of cellular stress (Coimbra et al., 2017), the associations reported between environmental stress in childhood and shorter telomere length suggest it reflects early-life cellular damage that may be carried into adulthood.

Little is known regarding the interpretation of biological age in children assessed at the transcriptome, proteomic and metabolomic levels, since few biological clocks are available for this age range using these data. To the best of our knowledge, only the study of Giallourou et al., 2020 has applied metabolomic data to provide a multivariate model of age in children, finding that growth-constrained infants lag in their metabolic maturity relative to their healthier peers. It is possible that biological clocks constructed using these data, particularly proteomics and metabolomics, support the life-course aging framework, where age acceleration in children represents a buildup of biological capital, since they are closer to the phenotype than the DNA-based epigenetic clocks and telomere length.

We hypothesized that age acceleration would be associated with child development. To test this and assess whether age acceleration is associated with beneficial or detrimental effects on child development, we have performed a comparative analysis of two established and two candidate

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Participant flowchart. See *Supplementary file 1* for details on quality control of molecular data at sample and feature levels.</figcaption></figure>
assessments of biological age within the pan-European Human Early Life Exposome (HELIX) cohort of children aged between 5 and 12 years. We systematically compared associations with developmental endpoints, including growth and adiposity, cognition, behavior, lung function and pubertal development, and common health risk factors, for telomere length, DNA methylation age, and two newly derived clocks: transcriptome age and immunometabolic age. Through this analysis, we aimed to clarify the interpretation of age acceleration in children, and more broadly develop new biological markers of overall developmental staging in children.

# Materials and methods

## Study population

This study population included children recruited into the European population-based HELIX exposome cohort (Maitre et al., 2018; Vrijheid et al., 2014), which was based on six on-going longitudinal population-based birth cohorts established in six countries across different parts of Europe (Born in Bradford [BiB; UK] Wright et al., 2013), Étude des Déterminants Pré et Postnatals du Développement et de la Santé de l'Enfant [EDEN; France] (Heude et al., 2016), Infancia y Medio Ambiente [INMA; Spain] (Guxens et al., 2012), Kaunas Cohort [KANC; Lithuania] (Grazuleviciene et al., 2009), Norwegian Mother, Father and Child Cohort Study [MoBa; Norway] (Magnus et al., 2016), and Mother-Child Cohort in Crete [RHEA; Greece] (Chatzi et al., 2017) covering singleton deliveries from 2003–2008. All children participated in a harmonized 'HELIX subcohort' clinical examination in their respective study centres during 2014–2015, where biological samples were collected. A subset of children (from all study centres apart from MoBa), attended a second clinical examination, as part of the 'HELIX panel study' approximately 6 months after the first 'HELIX subcohort' examination, where a similar suite of biological samples was collected. A full description of the HELIX follow-up methods and study population, including eligibility criteria and sample size calculations are available in Agier et al., 2019; Quanjer et al., 2012. In the current study, we included all children with available molecular data (Figure 1).

Prior to the start of HELIX, all six cohorts had undergone the required evaluation by national ethics committees and obtained all the required permissions for their cohort recruitment and follow-up visits. Each cohort also confirmed that relevant informed consent and approval were in place for the secondary use of data from pre-existing data. The work in HELIX was covered by new ethical approvals in each country and at enrolment in the new follow-up, participants were asked to sign a new informed consent form. Additionally, the current study was approved by the Imperial College Research Ethics Committee (Reference: 19IC5567).

## Biological sample collection and processing

Blood was collected at the end of the clinical examination of the child to ensure an approximate 3 hr (median=3.5 hr, SD=1.1 hr) fasting time since the last meal. Blood samples were collected using a 'butterfly' vacuum clip and local anesthetic and processed into a variety of sample matrices, including plasma, whole blood for RNA extraction (Tempus tubes - Life Technologies, USA), red cells, and buffy coat for DNA extraction. These samples were frozen at -80 °C under optimized and standardized procedures until analysis.

DNA was obtained from children's peripheral blood (buffy coat) collected in EDTA tubes. DNA was extracted using the Chemagen kit (Perkin Elmer, USA) in batches of 12 samples within each cohort. DNA concentration was determined in a Nanodrop 1000 UV-Vis Spectrophotometer (Thermo Fisher Scientific, USA) and also with Quant-iTTM PicoGreen dsDNA Assay Kit (Life Technologies, USA). DNA extraction was repeated in around 8% of the blood samples as the DNA quantity or quality of the first extraction was low. Less than 1.5% of the samples were finally excluded due to low quality.

RNA was extracted from whole blood samples collected in Tempus tubes (Thermo Fisher Scientific, USA) using MagMAX for Stabilized Blood Tubes RNA Isolation Kit. The quality of RNA was evaluated with a 2100 Bioanalyzer (Agilent Technologies, USA), and the concentration with a NanoDrop 1000 UV-Vis Spectrophotometer. Samples classified as good RNA quality (78.67%) had a similar RNA pattern at visual inspection in the Bioanalyzer, an RNA Integrity Number (RIN)>5, and a concentration >10 ng/ul. Mean values (standard deviation, SD) for the RIN, concentration (ng/ul), Nanodrop 260/280 ratio and Nanodrop 260/230 ratio were: 7.05 (0.72), 109.07 (57.63), 2.15 (0.16), and 0.61 (0.41).
During the clinical examination, two spot urine samples (one before bedtime and one first morning void) were brought by the participants to the research centre in cool packs and stored at 4 °C until processing. Urine samples of the night before the visit and the first-morning void on the day of the visit were combined to provide two urine samples, representing the last night-time and first morning voids, which were collected on the evening and morning before the clinical examination and were subsequently pooled to generate a more representative sample of the last 24 hr for metabolomic analysis (N=1107). Either the night-time void (N=37) or morning void (N=48) sample was analyzed in cases where a pooled sample was missing (Maitre et al., 2018).

## Telomere length measurement

Blood average relative telomere length was measured by a modified qPCR protocol as described previously (Cawthon, 2009). Telomere and single copy-gene reaction mixture and PCR cycles used can be found in Martens et al., 2016. All measurements were performed in triplicate on a 7900HT Fast Real-Time PCR System (Applied Biosystems) in a 384-well format. On each run, a six-point serial dilution of pooled DNA was run to assess PCR efficiency as well as eight inter-run calibrators to account for the inter-run variability. Relative telomere lengths were calculated using qBase software (Biogazelle, Zwijnaarde, Belgium) and were expressed as the ratio of telomere copy number to single-copy gene number (T/S) relative to the average T/S ratio of the entire sample set. We achieved CV's within triplicates of the telomere runs, single-copy gene runs, and T/S ratios of 0.84%, 0.43%, and 6.4%, respectively.

## DNA methylation

Blood DNA methylation was assessed with the Infinium HumanMethylatio450 beadchip (Illumina, USA) at the University of Santiago de Compostela - Spanish National Genotyping Center (CeGen-USC) (Spain). 700 ng of DNA were bisulfite-converted using the EZ 96-DNA kit (Zymo Research, USA) following the manufacturer's standard protocol. All samples of the study were randomized considering sex and cohort. In addition, each plate contained a HapMap control sample and 24 HELIX inter-plate duplicates were included.

After an initial inspection of the quality of the methylation data with the MethylAid package (van Iterson et al., 2014), probes with a call rate <95% based on a detection p-value of 1e-16, and samples with a call rate <98% were removed (Lehne et al., 2015). Samples with discordant sex were eliminated from the study as well as duplicates with inconsistent genotypes and samples with inconsistent genotypes with respect to existing genome-wide genotyping array data. Methylation data were normalized using the functional normalization method with prior background correction with Noob (Fortin et al., 2014). Then, some probes were filtered out: control probes, probes to detect single nucleotide polymorphisms (SNPs), probes to detect methylation in non-CpG sites, probes located in sexual chromosomes, cross hybridizing probes (Chen et al., 2013), probes containing an SNP at any position of the sequence with a minor allele frequency (MAF) >5% and probes with an SNP at the CpG site or at the single base extension (SBE) at any MAF in the combined population from 1000 Genomes Project (Supplementary file 1). Batch effect (slide) was corrected using the ComBat R package (Johnson et al., 2007). CpGs were annotated with the IlluminaHumanMethylation450kanno. ilmn12.hg19 R package (Hansen KD, 2012).

## Transcriptome analysis

Gene expression was assessed using the GeneChip Human Transcriptome Array 2.0 (HTA 2.0) from Affymetrix (USA) at the University of Santiago de Compostela (USC) (Spain). Briefly, RNA samples were concentrated or evaporated in order to reach the required RNA input concentration (200 ng of total RNA). Amplified and biotinylated sense-strand DNA targets were generated from total RNA. Microarrays were hybridized according to the Affymetrix recommendations using the Affymetrix labeling and hybridization kits. All samples were randomized within each batch considering sex and cohort. Two different types of control RNA samples (HeLa and FirstChoice Human Brain Reference RNA (Thermo Fisher Scientific, USA)) were included in each batch, but they were hybridized only in the first batches.

Raw data were extracted with the Affymetrix AGCC software and normalized with the GCCN (SST-RMA) algorithm at the gene level (http://tools.thermofisher.com/content/sfs/brochures/sst_gccn_whitepaper.pdf). Annotation of transcripts clusters (TCs) to genes was done with the Affymetrix
Expression Console software using the HTA-2.0 Transcript Cluster Annotations Release na36 (hg19). A transcript cluster is defined as a group of one or more probes covering a region of the genome reflecting all the exonic transcription evidence known for the region and corresponding to a known or putative gene. Four samples with discordant sex were detected with the MassiR R package (*Buck-berry et al., 2014*) and excluded. Control probes, and TCs in sexual chromosomes and without chromosome information were filtered out. Batch effect (slide) was corrected using the ComBat R package (*Johnson et al., 2007*). To determine the TC call rate, 10 constitutive or best probes based on probe scoring and cross-hybridization potential were selected per TC. Probe Detection Above Background (DABG) p-values were computed based on the rank order against the background probe set intensities. Probe level p-values were combined into a TC level p-value using the Fisher equation. TCs with a DABG p-value <0.05 were defined as detected. Three samples with low call rate (<40%) as well as TCs with a call rate <1% were excluded from the dataset (*Supplementary file 1*). Gene expression values were log2 transformed.

## Proteome analysis

Plasma protein levels were assessed using the antibody-based multiplexed platform from Luminex. Three kits targeting 43 unique candidate proteins were selected (Thermo Fisher Scientific, USA): Cytokines 30-plex (Catalog Number (CN): LHC6003M), Apolipoprotein 5-plex (CN: LHP0001M) and Adipokine 15-plex (CN: LHC0017M).

All samples were randomized and blocked by the cohort prior to measurement. For quantification, an eight-point calibration curve per plate was performed with protein standards provided in the Luminex kit and following procedures described by the vendor. Commercial heat-inactivated, sterile-filtered plasma from human male AB plasma (Sigma-Aldrich, USA) was used as constant samples to control for intra- and inter-plate variability. Four control samples were added per plate. All samples, including controls, were diluted 1/2 for the 30-plex kit, 1/4 for the 15-plex kit, and 1/2500 for the 5-plex kit.

Raw intensities obtained with the xMAP and Luminex system for each plasma sample were converted to pg/ml using the calculated standard curves of each plate and accounting for the dilutions made prior to measurement. The percentages of coefficients of variation (CV%) for each protein by plate ranged from 3-6%. The limit of detection (LOD) and the lower and upper limit of quantification (LOQ1 and LOQ2, respectively) were estimated by plate, and then averaged. Only proteins with >30% of measurements in the linear range of quantification were kept in the database and the others were removed. Seven proteins were measured twice (in two different multiplex kits). We kept the measure with higher quality. The 36 proteins that passed the quality control criteria mentioned above were log2 transformed (*Vives-Usano et al., 2020*). Then, the plate batch effect was corrected by subtracting the plate-specific average for each protein minus the overall average of all plates for that protein. After that, values below the LOQ1 and above the LOQ2 were imputed using a truncated normal distribution implemented in the truncdist R package (*Nadarajah and Kotz, 2006*). Twenty samples were excluded due to having ten or more proteins out of the linear range of quantification (*Supplementary file 1*).

## Metabolomic analysis

The AbsoluteIDQTM p180 kit was chosen for serum analysis as it is a standardized, targeted LC-MS/MS assay, widely used for large-scale epidemiology studies and its inter-laboratory reproducibility has been demonstrated by several independent laboratories (*Siskos et al., 2017*). Serum samples were quantified using the AbsoluteIDQTM p180 kit following the manufacturer's protocol (User Manual UM_p180_AB_SCIEX_9, Biocrates Life Sciences AG) using LC-MS/MS; an Agilent HPLC 1100 liquid chromatography coupled to a SCIEX QTRAP 6500 triple quadrupole mass spectrometer. A full description of the HELIX metabolomics methods and data can be found elsewhere (*Lau et al., 2018*).

Briefly, the kit allows for the targeted analysis of 188 metabolites in the classes of amino acids, biogenic amines, acylcarnitines, glycerophospholipids, sphingolipids, and the sum of hexoses, covering a wide range of analytes and metabolic pathways in one targeted assay. The kit consists of a single sample processing procedure, with two separate analytical runs, a combination of liquid chromatography (LC) and flow injection analysis (FIA) coupled to tandem mass spectrometry (MS/MS). Isotopically labeled and chemically homologous internal standards were used for quantification.
The AbsoluteIDQ p180 data of serum samples were acquired in 18 batches. Every analytical batch, in a 96-well plate format, included up to 76 randomized cohort samples. Also in every analytical batch, three sets of quality control samples were included, the NIST SRM 1950 plasma reference material (in four replicates), a commercially available serum QC material (CQC in two replicates, SeraLab, S-123-M-27485), and the QCs provided by the manufacturer in three concentration levels. The NIST SRM 1950 reference was used as the main quality control sample for the LC-MS/MS analysis. Coefficients of variation (CVs) for each metabolite were calculated based on the NIST SRM 1950 and also the limits of detection (LODs) were also used to assess the analytical performance of individual metabolites. Metabolite exclusion was based on a metabolite variable meeting two conditions: (1) CV of over 30% and (2) over 30% of the data are below LOD. Eleven out of the 188 serum metabolites detected were excluded as a result, leaving 177 serum metabolites to be used for further statistical analysis (Supplementary file 1). The mean coefficient of variation across the 177 LC-MS/MS detected serum metabolites was 16%. We also excluded one HELIX sample, which was hemolyzed.

Urinary metabolic profiles were acquired using ¹ H NMR spectroscopy according to Lau et al., 2018. In brief one-dimensional 600 MHz ¹ H NMR spectra of urine samples from each cohort were acquired on the same Bruker Avance III spectrometer operating at 14.1 Tesla within a period of 1 month. The spectrometer was equipped with a Bruker SampleJet system, and a 5 mm broad-band inverse configuration probe maintained at 300 K. Prior to analysis, cohort samples were randomized. Deuterated 3-(trimethylsilyl)-[2,2,3,3-d4]-propionic acid sodium salt (TSP) was used as an internal reference. Aliquots of the study pooled quality control (QC) sample were used to monitor analytical performance throughout the run and were analyzed at an interval of every 23 samples (i.e. four QC samples per well plate). The ¹ H NMR spectra were acquired using a standard one-dimensional solvent suppression pulse sequence. 44 metabolites were identified and quantified as described (Lau et al., 2018). The urinary NMR showed excellent analytical performance, the mean coefficient of variation across the 44 NMR-detected urinary metabolites was 11%. Data were normalized using the median fold change normalization method (Dieterle et al., 2006), which takes into account the distribution of relative levels of all 44 metabolites compared to the reference sample in determining the most probable dilution factor. An offset of ½ of the minimal value was applied and then concentration levels were expressed as log2.

## Building biological clocks

Child epigenetic age was calculated based on Horvath's Skin and Blood clock (Horvath et al., 2018) using the methylclock R package (Pelegí-Sisó et al., 2021).

New transcriptome and immunometabolic clocks were trained against chronological age on transcriptome data and concatenated proteomic and metabolomic data respectively, from the HELIX subcohort children through elastic net regression, using the glmnet R package (Zou and Hastie, 2005). All 'omic data was first mean-centred and univariate scaled. To tune hyperparameters alpha and lambda, we performed a line search for alpha between 0 and 1, in 0.1 increments, and each time found the optimal value of lambda based on minimization of cross-validated mean squared error, using the cvfits function and 10-fold cross-validation. The best-performing combination of alpha and lambda was reserved for fitting the final model.

Transcriptome data and concatenated proteomic and metabolomic data from the HELIX panel study children were reserved for testing performance (Pearson's r and mean absolute error with chronological age) of the derived clocks. Paired, one-tailed t-tests were used to test if biological age measures increased between the HELIX subcohort and subsequent HELIX panel clinical examinations.

## Developmental measurements

During the HELIX subcohort examination, height and weight were measured using regularly calibrated instruments and converted to BMI and height age-and-sex-standardized z-scores (zBMI and zHeight) using the international World Health Organization (WHO) reference curves (de Onis et al., 2007). Bioelectric impedance analyses were performed with the Bodystat 1500 (Bodystat Ltd.) equipment after 5 min of lying down. The proportion of fat mass was calculated using published age- and race-specific equations validated for use in children (Clasey et al., 2011).

Trained fieldwork technicians measured three cognitive domains in children using a battery of computer-based tests: fluid intelligence (Raven Coloured Progressive Matrices Test [CPM]), attention
function (Attention Network Test [ANT]), and working memory (N-Back task). Complete outcome descriptions are provided in Julvez et al., 2021. The CPM comprised a total of 36 items and we used the total number of correct responses as the outcome. A higher CPM scoring indicates better fluid intelligence. Fluid intelligence is the ability to solve novel reasoning problems and depends only minimally on prior learning. For ANT, we used the outcome of hit reaction time standard error (HRT-SE), a measure of response speed consistency throughout the test. A high HRT-SE indicates highly variable reaction time during the attention task and is considered a measure of inattentiveness (Forns et al., 2014). As the main parameter of N-Back, we used d prime (d') from the three-back colors test, a measure derived from signal detection theory calculated by subtracting the z-score of the false alarm rate from the z-score of the hit rate. A higher d' indicates more accurate test performance, i.e., better working memory (Forns et al., 2014). All examiners were previously trained following a standardized assessment protocol by the study expert psychologist. Furthermore, during the pilot phase, a coordinator visited each cohort site and checked for any potential error committed by the previously trained examiners.

Parents completed questionnaires related to child's behavior, including the Conner rating scale (N=1287) and child behavior checklist (CBCL, N=1298), within a week before the follow-up visit at 6-11 years of age. The 99-item CBCL/6-18 version for school children was used to obtain standardized parent reports of children's problem behaviors, translated and validated in each native language of the participating six cohort populations (Achenbach, 2001). The parents responded along a three-point scale with the code of 0 if the item is not true of the child, one for sometimes true, and two for often true. The internalizing score includes the subscales of emotionally reactive and anxious/depressed symptoms, as well as somatic complaints and symptoms of being withdrawn. The externalizing score includes attention problems and aggressive behaviors.

Lung function was measured by a spirometry test (EasyOne spirometer; NDD [New Diagnostic Design], Zurich, Switzerland), by trained research technicians using a standardized protocol. The child, sitting straight and equipped with a nose clip, was asked to perform at least six maneuvers (if possible). Details of the exclusion of unacceptable maneuvers and validation of acceptable spirometer curves is fully described in Agier et al., 2019. FEV1 percent predicted values were computed using the reference equations estimated by the Global Lung Initiative (Quanjer et al., 2012), standardized by age, height, sex, and ethnicity.

Parents of children aged 8 years or older completed an additional questionnaire based on the pubertal development scale (PDS) (Petersen et al., 1988). Boys were asked whether growth has not begun, barely begun, is definitely underway, or has finished on five dimensions: body hair, facial hair, voice change, skin change, and a growth spurt. Girls were asked the same questions about body hair, skin change, breast development, and growth spurt. Responses were coded on four-point scales (1=no development and 4=completed development). For girls, a yes-no question about the onset of menarche is weighted more heavily (1=no and 4=yes). For both genders, ratings are then averaged to create an overall score for physical maturation. Due to the young age of participants, we took the average scores and created a binary variable, to define whether puberty had started (PDS >1) or not (PDS=1).

## Covariates

During pregnancy and in the childhood HELIX subcohort examination information on the following key covariates was collected: cohort study centre (BiB, EDEN, INMA, MoBa, KANC, and RHEA), self-reported maternal education (primary school, secondary school, and a university degree or higher), self-reported ancestry (White European, Asian and Pakistani, or other), birth weight (continuous, kg), gestational age at delivery (continuous in weeks).

Information about the children's habitual diet was collected via a semi-quantitative food-frequency questionnaire (FFQ) covering the child's habitual diet, which was filled in by the parent attending the examination appointment. The FFQ, covering the past year, was developed by the HELIX research group, translated, and applied to all cohorts. For the Mediterranean Diet Quality Index (KIDMED index) (Serra-Majem et al., 2004), items positively associated with the Mediterranean diet pattern (11 items) were assigned a value of +1, while those negatively associated with the Mediterranean diet pattern (four items) were assigned a value of -1. The scores for all 15 items were summed, resulting in a total KIDMED score ranging from -4-11, with higher scores reflecting greater adherence to a Mediterranean diet.
The smoking status of the mother at any point during pregnancy was categorized into ‘non-active smoker, or ‘active smoker.’ Global exposure of the child to environmental tobacco smoke was defined based on the questionnaires completed by the parents into: ‘no exposure,’ no exposure at home neither in other places; ‘exposure:’ exposure in at least one place, at home or outside. Moderate-to-vigorous physical activity variable was created based on the physical activity questionnaire developed by the HELIX research group. It was defined as the amount of time children spent doing physical activities with intensity above three metabolic equivalent tasks (METs) and is expressed in units of min/day.

Family Affluence Score (FAS) (Boyce et al., 2006) was included based on questions from the subcohort questionnaire. A composite FAS score was calculated based on the responses to the next four items: (1) Does your family own a car, van, or truck? (2) Do you have your own bedroom for yourself? (3) During the past 12 months, how many times did you travel away on holiday with your family? (4) How many computers does your family own? A three-point ordinal scale was used, where FAS low (score 0, 1, 2) indicates low affluence, FAS medium (score 3, 4, 5) indicates middle affluence, and FAS high (score 6, 7, 8, 9) indicates high affluence FAS.

Family social capital-related questions were included in the HELIX questionnaire to capture different aspects of social capital, relating both to the cognitive (feelings about relationships) and structural (number of friends, number of organizations) dimensions and to bonding capital (close friends and family), bridging capital (neighborhood connections, looser ties) and linking capital (ties across power levels; for example, political membership). Family social capital was categorized into low, medium, and high based on terciles.

## Statistical analysis

All statistical analyses described here were performed among the HELIX subcohort children only. Since there were few missing covariate data (Supplementary file 2), a complete-case analysis was performed. Correlations between biological age measures and chronological age were calculated using Pearson's correlations. Partial correlations, adjusted for chronological age and cohort study centre, were applied to assess correlations between biological age measures.

In analysis with health risk factors and developmental outcomes, relative telomere length was multiplied by -1 to provide directions of effect consistent with the biological age clocks and univariate scaled to express effects in terms of SD change in telomere length. The markers derived from omic-based biological clocks were expressed as Δ age (clock-predicted age – chronological age). Associations between the biological age markers and developmental measures were estimated using linear regression, or logistic regression for the onset of puberty, with the developmental measure as the dependent variable. CBCL scores were log-transformed to achieve an approximately normal distribution. Continuous outcomes, apart from the BMI and height z-scores, were mean-centered and univariate scaled for the purposes of graphical representation. Associations between health risk factors and biological age markers were estimated using linear regression with the biological age marker as the dependent variable. All regression analyses were adjusted for chronological age, sex, ethnicity, and study centre.

We preformed four sensitivity analyses: Firstly, we repeated the analysis with health outcomes stratified by child sex, since the relationship between biological age and development may differ between boys and girls. Secondly, we further adjusted regression models for estimated cell counts (CD4T, CD8T, monocytes, B cells, NK cells, neutrophils, and eosinophils), since it has been proposed for epigenetic clocks that cell proportion adjustments allow estimation of effects on the intrinsic cellular aging rate, rather than the extrinsic rate outputted by blood-based biological clocks, which may be partly determined by age-related changes in cell composition (Chen et al., 2016). Blood cell type proportion was estimated from DNA methylation data using the He et al., 2014 reference panel as implemented in meffil package (Costa et al., 2015). Third, we assessed the effects of further adjustment for health risk factors identified as associated with any of the biological age markers (family affluence and social capital, birthweight, maternal active smoking, and child passive smoking). In our main analysis, we have not adjusted for these factors as our assumption is that the effects of health risk factors on child development are mediated through biological age. However, an alternative assumption is that health risk factors exert independent effects on both biological age and developmental outcomes, which would require adjustment for these factors to estimate the direct effects of biological age on developmental outcomes. Finally, we stratified by study centre to
check the consistency of effects among observed associations with developmental outcomes across study centres.

Due to the exploratory nature of the analysis, we report and discuss associations significant at both the nominal significance threshold (unadjusted p<0.05) and after correction for 5% false discovery rate using the *Benjamini and Hochberg, 1995* method, calculated across all computed associations (i.e. multiple testing corrected p-value <0.05).

We performed overrepresentation analyses (ORA) among KEGG and REACTOME pathways and gene ontology (GO) sets of all transcripts contributing to the transcriptome clock using the Consensus-pathDB online tool (http://consensuspathdb.org/). A pathway or GO set was considered significantly enriched if FDR-corrected p-values were smaller than 0.05 and included at least 3 genes. Additionally, to assess concordance with gene expression changes with age in adults, we tested the enrichment of all transcripts contributing to the transcriptome clock among age-associated transcripts reported by *Peters et al., 2015*, using a hypergeometric test using the R 'phyper' function.

All analyses were performed in R version 4.1.2.

# Results

## Sample characteristics

We used blood or urine-derived measurements from the pan-European HELIX cohort. This included blood telomere length (N=1162), blood DNA methylation (N=1173, 450 K CpGs), blood gene expression (N=1007, 50 K genes), and proteins and metabolites (N=1152, 36 plasma proteins, 177 serum metabolites, and 44 urinary metabolites), with 869 children overlapping across all measurements. Each subsample included around 55% boys, 89% children of white European ancestry, and a mean age of around 8 years (range 5–12 years). Around 51% of mothers of the HELIX children in each subsample had a high education level. The HELIX cohort included children from six study centres based in the UK, Spain, Greece, Lithuania, France, and Norway, with each centre contributing between 11 and 24% to each subsample (see Table 1 for sample characteristics).

## Biological age marker performance

We included two established markers of biological age, telomere length, DNA methylation age, and developed two new candidate biological age markers, transcriptome age and 'immunometabolic' age (Figure 2). DNA methylation age was calculated using the published Skin and blood Horvath clock (Horvath et al., 2018) to allow greater comparison to the wider literature, including in adults. We previously reported this epigenetic clock to show the best performance in chronological age prediction within the HELIX cohort (de Prado-Bert et al., 2021). Since no published applicable transcriptome, proteome or metabolome clocks were available for the age range of our sample, we trained two new biological clocks using these data in the HELIX cohort, through elastic net regression and cross-validation. We combined the proteome and metabolome data into a single immunometabolic age clock, since the available proteomic data included biomarkers targeting both metabolic and inflammatory functions, both omic types represent final products of gene regulation, and since the metabolic and immune systems are closely linked (Hotamisligil, 2017).

The correlation between telomere length and chronological age was weak but statistically significant (r=-0.07, p=0.02). Correlations with chronological age were r=0.85 for DNA methylation age, r=0.94 for transcriptome age, and r=0.86 for immunometabolic age (Figure 2).

We validated the transcriptome and immunometabolic clocks using cross-validation within the HELIX subcohort (cross-validated r of 0.87 and 0.82, respectively) and further tested in a subset of children who attended a second clinic visit approximately 0.5 years after the main follow-up visit (standard deviation (SD)=0.18 years) as part of the HELIX panel study. Correlations in this test set were r=0.93 for transcriptome age (N=128) and r=0.84 for immunometabolic age (N=151) (Figure 2). Predicted biological age increased by mean of 0.33 years (SD=0.58) for transcriptome age (t-test, p=3 e-5) and a mean of 0.22 years (SD: 0.59 years) for immunometabolic age (t-test, p=2 e-5) between the first and second visits (Figure 2—figure supplement 1). Correlations were significant (p<0.05) within each study centre for both clocks, except for immunometabolic age for children from the BiB (UK) cohort (Figure 2—figure supplement 2 and Figure 2—figure supplement 3).
<table><caption>**Table 1.** Summary Statistics for the study population.</caption><thead><tr><th></th><th>Telomere Length</th><th>DNA methylation age</th><th>Transcript-ome age</th><th>Immuno-metabolic age</th></tr><tr><th></th><th>N (%) or Mean (SD)</th><th>N (%) or Mean (SD)</th><th>N (%) or Mean (SD)</th><th>N (%) or Mean (SD)</th></tr></thead><tbody><tr><td>N</td><td>1162</td><td>1173</td><td>1007</td><td>1152</td></tr><tr><td colspan="5">Demographic factors</td></tr><tr><td>Age (years)</td><td>7.84 (1.54)</td><td>7.84 (1.54)</td><td>7.90 (1.50)</td><td>7.86 (1.55)</td></tr><tr><td>Sex-Male</td><td>639 (55)</td><td>644 (54.9)</td><td>547 (54.3)</td><td>628 (54.5)</td></tr><tr><td>Sex-Female</td><td>523 (45)</td><td>529 (45.1)</td><td>460 (45.7)</td><td>524 (45.5)</td></tr><tr><td>Ethnicity-White</td><td>1039 (89.4)</td><td>1048 (89.3)</td><td>905 (89.9)</td><td>1032 (89.6)</td></tr><tr><td>Ethnicity-Pakistani/ Asian</td><td>96 (8.3)</td><td>98 (8.4)</td><td>76 (7.5)</td><td>93 (8.1)</td></tr><tr><td>Ethnicity -Other</td><td>27 (2.3)</td><td>27 (2.3)</td><td>26 (2.6)</td><td>27 (2.3)</td></tr><tr><td>Cohort-BIB</td><td>200 (17.2)</td><td>203 (17.3)</td><td>162 (16.1)</td><td>191 (16.6)</td></tr><tr><td>Cohort-EDEN</td><td>145 (12.5)</td><td>146 (12.4)</td><td>109 (10.8)</td><td>149 (12.9)</td></tr><tr><td>Cohort-INMA</td><td>212 (18.2)</td><td>215 (18.3)</td><td>184 (18.3)</td><td>201 (17.4)</td></tr><tr><td>Cohort-KANC</td><td>196 (16.9)</td><td>198 (16.9)</td><td>151 (15)</td><td>197 (17.1)</td></tr><tr><td>Cohort-MOBA</td><td>211 (18.2)</td><td>212 (18.1)</td><td>245 (24.3)</td><td>222 (19.3)</td></tr><tr><td>Cohort-RHEA</td><td>198 (17)</td><td>199 (17)</td><td>156 (15.5)</td><td>192 (16.7)</td></tr><tr><td colspan="5">Prenatal factors</td></tr><tr><td>maternal non-active smoker during pregnancy</td><td>988 (85)</td><td>998 (85.1)</td><td>859 (85.3)</td><td>981 (85.2)</td></tr><tr><td>Maternal active smoker during pregnancy</td><td>174 (15)</td><td>175 (14.9)</td><td>148 (14.7)</td><td>171 (14.8)</td></tr><tr><td>Birthweight (kg)</td><td>3.37 (0.5)</td><td>3.37 (0.5)</td><td>3.38 (0.52)</td><td>3.38 (0.5)</td></tr><tr><td>Gestational age (weeks)</td><td>39.57 (1.67)</td><td>39.58 (1.67)</td><td>39.59 (1.75)</td><td>39.59 (1.66)</td></tr><tr><td colspan="5">Family Capital</td></tr><tr><td>Maternal Education (low)</td><td>165 (14.7)</td><td>166 (14.7)</td><td>140 (14.4)</td><td>157 (14.1)</td></tr><tr><td>Maternal Education (medium)</td><td>391 (34.8)</td><td>394 (34.8)</td><td>328 (33.8)</td><td>391 (35.1)</td></tr><tr><td>Maternal Education (high)</td><td>568 (50.5)</td><td>573 (50.6)</td><td>503 (51.8)</td><td>565 (50.8)</td></tr><tr><td>Family Affluence (low)</td><td>133 (11.5)</td><td>135 (11.5)</td><td>112 (11.1)</td><td>128 (11.1)</td></tr><tr><td>Family Affluence (medium)</td><td>462 (39.8)</td><td>466 (39.8)</td><td>394 (39.2)</td><td>450 (39.1)</td></tr><tr><td>Family Affluence (high)</td><td>565 (48.7)</td><td>570 (48.7)</td><td>499 (49.7)</td><td>572 (49.7)</td></tr><tr><td>Family Social Capital (low)</td><td>513 (47.7)</td><td>516 (47.5)</td><td>422 (45.8)</td><td>496 (46.7)</td></tr><tr><td>Family Social Capital (medium)</td><td>264 (24.6)</td><td>269 (24.8)</td><td>228 (24.7)</td><td>259 (24.4)</td></tr><tr><td>Family Social Capital (high)</td><td>298 (27.7)</td><td>301 (27.7)</td><td>272 (29.5)</td><td>307 (28.9)</td></tr><tr><td colspan="5">Child factors</td></tr></tbody></table>

Table 1 continued on next page
<table><caption>Table 1 continued</caption><thead><tr><th></th><th>Telomere Length</th><th>DNA methylation age</th><th>Transcript-ome age</th><th>Immuno-metabolic age</th></tr></thead><tbody><tr><td>No passive smoke exposure</td><td>723 (63.8)</td><td>732 (63.9)</td><td>639 (64.5)</td><td>718 (63.8)</td></tr><tr><td>Passive smoke exposure</td><td>411 (36.2)</td><td>413 (36.1)</td><td>351 (35.5)</td><td>407 (36.2)</td></tr><tr><td>Physical Activity-Low</td><td>418 (36.9)</td><td>420 (36.8)</td><td>349 (35.3)</td><td>416 (37.1)</td></tr><tr><td>Physical Activity-Medium</td><td>336 (29.7)</td><td>341 (29.9)</td><td>295 (29.9)</td><td>330 (29.4)</td></tr><tr><td>Physical Activity-High</td><td>378 (33.4)</td><td>381 (33.4)</td><td>344 (34.8)</td><td>375 (33.5)</td></tr><tr><td>KIDMED diet score</td><td>2.81 (1.77)</td><td>2.82 (1.78)</td><td>2.88 (1.77)</td><td>2.84 (1.76)</td></tr><tr><td>Developmental measures</td><td></td><td></td><td></td><td></td></tr><tr><td>Height z-score</td><td>0.4 (0.97)</td><td>0.39 (0.98)</td><td>0.39 (0.96)</td><td>0.4 (0.98)</td></tr><tr><td>BMI z-score</td><td>0.43 (1.2)</td><td>0.43 (1.2)</td><td>0.4 (1.15)</td><td>0.42 (1.18)</td></tr><tr><td>Adiposity (BIA fat-mass %)</td><td>6.76 (4.01)</td><td>6.77 (4.01)</td><td>6.52 (3.9)</td><td>6.72 (3.95)</td></tr><tr><td>Working memory (3-back d')</td><td>1.1 (1.01)</td><td>1.1 (1.01)</td><td>1.13 (1)</td><td>1.1 (1.01)</td></tr><tr><td>Inattentiveness (ANT-HRT)</td><td>301.97 (90.38)</td><td>301.93 (90.46)</td><td>297.69 (89.36)</td><td>301.35 (89.84)</td></tr><tr><td>Fluid Intelligence (CPM)</td><td>25.87 (6.33)</td><td>25.86 (6.32)</td><td>26.12 (6.26)</td><td>25.95 (6.3)</td></tr><tr><td>Internalizing behaviors (CBCL)</td><td>6.49 (5.9)</td><td>6.48 (5.9)</td><td>6.36 (5.89)</td><td>6.52 (5.87)</td></tr><tr><td>Externalizing behaviors (CBCL)</td><td>6.81 (6.5)</td><td>6.82 (6.51)</td><td>6.67 (6.49)</td><td>6.74 (6.42)</td></tr><tr><td>Lung Function (FEV1)</td><td>99.26 (13.46)</td><td>99.25 (13.47)</td><td>99.16 (13.02)</td><td>99.17 (13.47)</td></tr><tr><td>Puberty not started</td><td>250 (46.6)</td><td>252 (46.5)</td><td>254 (49.7)</td><td>260 (48)</td></tr><tr><td>Puberty started (PDS >1)</td><td>287 (53.4)</td><td>290 (53.5)</td><td>257 (50.3)</td><td>282 (52)</td></tr></tbody></table>

The immunometabolic age clock was composed of 135 predictors including 20 proteins, 79 serum metabolites, and 36 urinary metabolites (Supplementary file 3). The transcriptome clock was composed of 1445 genes, 652 of which were annotated to Gene Symbols ( Supplementary file 4). The transcriptome clock genes were enriched (false discovery rate (FDR)-corrected p<0.05) in 'ribosome' and 'ribosome biogenesis' KEGG pathways (Supplementary file 5) and the following level 2 Gene Ontology biological process terms: 'leukocyte activation,' 'movement of cell or subcellular component,' 'leukocyte migration,' 'cell activation,' and 'secretion by cell' (Supplementary file 6). We also tested the enrichment of transcriptome clock predictors among genes reported by a large meta-analysis of age in adults (Peters et al., 2015): among the 1406 reported age-associated genes that could be matched to our measured genes, 43 were included in our transcriptome clock (hypergeometric enrichment test, p=0.052). We note that since a common definition of markers of biological age is that they should be associated with age-related disease and mortality (Moskalev, 2019) these new clocks may only currently be considered 'candidate' biological age markers. However, we have referred to both the established and candidate markers as biological age markers throughout to simplify the presentation.

**Figure 3** shows partial correlations, adjusted for chronological age and study centre, between the biological age markers. Only null to weak correlations were observed, with significant correlations between telomere length and DNA methylation age ($r=-0.06$, p=0.04) and between transcriptome age and immunometabolic age ($r=0.08$, p=0.01).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Study design schematic. Source data for reproducing correlation plots are provided in **Figure 2—source data 1**. The online version of this article includes the following source data and figure supplement(s) for figure 2:</figcaption></figure>

**Source data 1.** Source data for reproducing correlation plots in **Figure 2**.

Figure 2 continued on next page
Figure 2 continued

**Figure supplement 1.** Comparison between immunometabolic and transcriptome age between first and second study visits.

**Figure supplement 2.** Age Prediction by study centre of transcriptome age.

**Figure supplement 3.** Age Prediction by study centre of immunometabolic age.

## Biological clock associations with health risk factors

*Table 2* shows associations, adjusted for chronological age, sex, study centre and ethnicity, between health risk factors and the biological age markers. The markers derived from omic-based biological clocks are expressed as Δ age (clock-predicted age – chronological age) and since the adjustment set included chronological age, effects can be interpreted as years of age acceleration as often defined (Jansen *et al.*, 2021).

Nominally significant associations were observed for the following health risk factors: Telomere length was longer among girls compared to boys (p=3 e-06) and among children of high affluence families (p=0.008). DNA methylation Δ age was higher among children of mothers who actively smoked during pregnancy (p=0.018) and children exposed to passive smoke (p=0.023), while

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Correlations between biological age indicators. Heatmap shows partial Pearson's correlations, adjusted for chronological age and study centre. \* indicates p<0.05. Source data for reproducing plots is provided in **Figure 3—source data 1**.</figcaption></figure>

The online version of this article includes the following source data for figure 3:

**Source data 1.** Source data for reproducing Figure 3.
<caption>**Table 2.** Associations between health risk factors and biological age measures.</caption>

Estimates were calculated using linear regression, adjusted for chronological age, sex, ethnicity, and study centre. Bold indicates p<0.05 and *indicates FDR <5%. Telomere length is expressed as a standard deviation (SD) decrease in length (multiplied by -1) to provide estimates indicative of accelerated biological age, as the other biological age indicators. Telomere Length N=1162, DNA methylation age N=1173, Transcriptome age N=1007, Immunometabolic age N=1152.

<table><thead><tr><th rowspan="2"></th><th colspan="2">Telomere Length</th><th colspan="2">DNA methylation age</th><th colspan="2">Transcriptome age</th><th colspan="2">Immunometabolic age</th></tr><tr><th>SD Decrease (95% CI)</th><th>p-value</th><th>Increase in years Δ Age (95% CI)</th><th>p-value</th><th>Increase in years Δ Age (95% CI)</th><th>p-value</th><th>Increase in years Δ Age (95% CI)</th><th>p-value</th></tr></thead><tbody><tr><td>Sex-Male</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Sex-Female</td><td>-0.27 (-0.39, -0.16)</td><td>3.30E-06*</td><td>0.07 (-0.01, 0.16)</td><td>0.1</td><td>0 (-0.01, 0.02)</td><td>0.73</td><td>0.06 (-0.01, 0.13)</td><td>0.086</td></tr><tr><td colspan="9"><strong>Prenatal factors</strong></td></tr><tr><td>maternal non-active smoker during pregnancy</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="9"><strong>Maternal active smoker during pregnancy</strong></td></tr><tr><td>Birthweight (kg)</td><td>-0.098 (-0.218, 0.023)</td><td>0.11</td><td>-0.021 (-0.114, 0.072)</td><td>0.66</td><td>0.005 (-0.01, 0.02)</td><td>0.51</td><td>0.102 (0.027, 0.177)</td><td>0.0075</td></tr><tr><td>Gestational age (weeks)</td><td>-0.012 (-0.048, 0.024)</td><td>0.52</td><td>0.013 (-0.015, 0.041)</td><td>0.35</td><td>0 (-0.005, 0.004)</td><td>0.89</td><td>0.018 (-0.005, 0.04)</td><td>0.12</td></tr><tr><td colspan="9"><strong>Family Capital</strong></td></tr><tr><td>Maternal Education (low)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Maternal Education (medium)</td><td>-0.06 (-0.26, 0.13)</td><td>0.53</td><td>0.02 (-0.14, 0.17)</td><td>0.84</td><td>0.01 (-0.02, 0.03)</td><td>0.61</td><td>0.08 (-0.04, 0.2)</td><td>0.21</td></tr><tr><td>Maternal Education (high)</td><td>-0.1 (-0.29, 0.1)</td><td>0.32</td><td>-0.07 (-0.22, 0.08)</td><td>0.37</td><td>0 (-0.02, 0.03)</td><td>0.85</td><td>0.12 (0, 0.24)</td><td>0.051</td></tr><tr><td>Family Affluence (low)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Family Affluence (medium)</td><td>-0.15 (-0.34, 0.05)</td><td>0.13</td><td>-0.11 (-0.26, 0.03)</td><td>0.13</td><td>0 (-0.03, 0.02)</td><td>0.85</td><td>0.02 (-0.1, 0.14)</td><td>0.8</td></tr></tbody></table>

Table 2 continued on next page
<table><caption>Table 2 continued</caption><thead><tr><th></th><th>Telomere Length</th><th>DNA methylation age</th><th>TranScriptome age</th><th>Immunometabolic age</th></tr></thead><tbody><tr><td>Family Affluence (high)</td><td>-0.27 (-0.47, -0.07)</td><td><strong>0.0081</strong></td><td>-0.14 (-0.29, 0.02)</td><td>0.083</td><td>0.01 (-0.01, 0.04)</td><td>0.35</td><td>0.09 (-0.04, 0.21)</td><td>0.17</td></tr><tr><td>Family Social Capital (low)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Family Social Capital (medium)</td><td>-0.06 (-0.21, 0.09)</td><td>0.45</td><td>-0.03 (-0.14, 0.09)</td><td>0.62</td><td>0.02 (0.01, 0.04)</td><td><strong>0.012</strong></td><td>-0.04 (-0.14, 0.05)</td><td>0.36</td></tr><tr><td>Family Social Capital (high)</td><td>-0.15 (-0.3, 0)</td><td>0.054</td><td>-0.12 (-0.23, 0)</td><td><strong>0.048</strong></td><td>0.02 (0.01, 0.04)</td><td><strong>0.011</strong></td><td>-0.06 (-0.15, 0.04)</td><td>0.25</td></tr><tr><td>Child factors</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No passive smoke exposure</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Passive smoke exposure</td><td>0.05 (-0.08, 0.18)</td><td>0.42</td><td>0.11 (0.02, 0.21)</td><td><strong>0.023</strong></td><td>0.01 (0, 0.03)</td><td>0.16</td><td>-0.01 (-0.09, 0.07)</td><td>0.76</td></tr><tr><td>Physical Activity-Low</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Physical Activity-Medium</td><td>0.09 (-0.06, 0.23)</td><td>0.25</td><td>-0.08 (-0.2, 0.03)</td><td>0.15</td><td>-0.01 (-0.03, 0.01)</td><td>0.17</td><td>0.03 (-0.06, 0.12)</td><td>0.56</td></tr><tr><td>Physical Activity-High</td><td>0.14 (-0.01, 0.29)</td><td>0.067</td><td>-0.1 (-0.22, 0.01)</td><td>0.08</td><td>0 (-0.02, 0.01)</td><td>0.69</td><td>-0.06 (-0.15, 0.04)</td><td>0.24</td></tr><tr><td>KIDMED diet score</td><td>-0.03 (-0.064, 0.005)</td><td>0.092</td><td>0.005 (-0.022, 0.031)</td><td>0.74</td><td>0.004 (-0.001, 0.008)</td><td>0.10</td><td>-0.005 (-0.027, 0.016)</td><td>0.64</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Associations between biological age measures and developmental measures. Estimates were calculated using linear regression, adjusted for chronological age, sex, ethnicity, and study centre. \*indicates FDR <5%. Telomere length is expressed as a standard deviation (SD) decrease in length (multiplied by -1) to provide estimates indicative of accelerated biological age, as the other biological age indicators. Error bars show 95% confidence intervals. See **Table 3** for numbers included in each analysis and exact point estimates and confidence intervals.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 4:

**Figure supplement 1.** Associations between biological age measures and developmental measures, stratified by sex.

**Figure supplement 2.** Associations between telomere length and developmental measures adjusted for (A) chronological age, sex, ethnicity, and study centre; (B) as for A plus estimated cell counts; (C) as for A plus family affluence and social capital, birthweight, maternal active smoking, and child passive smoking; (D) as for C plus estimated cell counts.

**Figure supplement 3.** Associations between DNA methylation Δ age and developmental measures adjusted for (A) chronological age, sex, ethnicity, and study centre; (B) as for A plus estimated cell counts; (C) as for A plus family affluence and social capital, birthweight, maternal active smoking, and child passive smoking; (D) as for C plus estimated cell counts.

**Figure supplement 4.** Associations between transcriptome Δ age and developmental measures adjusted for (A) chronological age, sex, ethnicity, and study centre; (B) as for A plus estimated cell counts; (C) as for A plus family affluence and social capital, birthweight, maternal active smoking, and child passive smoking; (D) as for C plus estimated cell counts.

**Figure supplement 5.** Associations between immunometabolic Δ age and developmental measures adjusted for (A) chronological age, sex, ethnicity, and study centre; (B) as for A plus estimated cell counts; (C) as for A plus family affluence and social capital, birthweight, maternal active smoking, and child passive smoking; (D) as for C plus estimated cell counts.

**Figure supplement 6.** Associations between biological age measures and developmental measures, stratified by study centre (adjusted for chronological age, sex, and ethnicity).

DNA methylation Δ age was lower among children from families with high social capital (p=0.048). Conversely, transcriptome Δ age was positively associated with medium and high (p=0.011) family social capital. Immunometabolic Δ age was associated with higher birthweight (p=0.0075). Only the association between longer telomere length and female sex passed the FDR correction.

## Biological age associations with development

Figure 4 and Table 3 shows associations, adjusted for chronological age, sex, study centre and ethnicity, between the biological age markers and developmental outcomes related to growth and adiposity, cognition, behavior, lung function, and the onset of puberty.

Several developmental outcomes were associated with biological age markers after FDR correction: DNA methylation and immunometabolic Δ age were associated with greater height z-score (p=6 e-6 and p=4 e-11, respectively) and greater fat mass % (p=0.0004 and p=5 e-6, respectively). All biological age markers were associated with greater BMI z-score (telomere length p=8 e-4, DNA methylation Δ age p=8 e-5, transcriptome Δ age p=0.005, immunometabolic Δ age p=4 e-19). Furthermore,
<table><caption>**Table 3.** Associations between biological age measures and developmental measures. Estimates were calculated using linear regression, adjusted for chronological age, sex, ethnicity, and study centre.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Telomere Length</th><th colspan="3">DNA methylation age</th><th colspan="3">Transcriptome age</th><th colspan="3">Immunometabolic age</th></tr><tr><th>SD increase / odds ratio * per SD shortening (95% CI)</th><th>p-value</th><th>N</th><th>SD increase / odds ratio per year increase in Δ age (95% CI)</th><th>p-value</th><th>N</th><th>SD increase / odds ratio per year increase in Δ age (95% CI)</th><th>p-value</th><th>N</th><th>SD increase/ odds ratio per year increase in Δ age (95% CI)</th><th>p-value</th></tr></thead><tbody><tr><td>Height z-score</td><td>1162</td><td>0 (-0.05, 0.06)</td><td>0.89</td><td>1173</td><td>0.17 (0.09, 0.24)</td><td>6.20E-06*</td><td>1007</td><td>0.66 (0.13, 1.18)</td><td>0.014</td><td>1152</td><td>0.31 (0.22, 0.4)</td><td>4.30E-11*</td></tr><tr><td>BMI z-score</td><td>1162</td><td>0.12 (0.05, 0.19)</td><td>0.00082*</td><td>1173</td><td>0.18 (0.09, 0.27)</td><td>7.70E-05*</td><td>1007</td><td>0.9 (0.27, 1.53)</td><td>0.005*</td><td>1152</td><td>0.5 (0.4, 0.61)</td><td>3.80E-19*</td></tr><tr><td>Adiposity (BIA fat-mass %)</td><td>1153</td><td>0.07 (0.02, 0.12)</td><td>0.0093*</td><td>1164</td><td>0.12 (0.05, 0.19)</td><td>0.0004*</td><td>999</td><td>0.45 (-0.05, 0.94)</td><td>0.079</td><td>1144</td><td>0.2 (0.11, 0.28)</td><td>5.00E-06*</td></tr><tr><td>Working memory (3-back d')<sup>†</sup></td><td>882</td><td>-0.03 (-0.09, 0.03)</td><td>0.37</td><td>890</td><td>0.05 (-0.03, 0.13)</td><td>0.19</td><td>784</td><td>0.06 (-0.51, 0.63)</td><td>0.84</td><td>876</td><td>0.15 (0.05, 0.26)</td><td>0.0036*</td></tr><tr><td>Inattentiveness (ANT-HRT)</td><td>1142</td><td>-0.01 (-0.05, 0.04)</td><td>0.78</td><td>1153</td><td>0.07 (0.01, 0.13)</td><td>0.03</td><td>997</td><td>-0.39 (-0.85, 0.08)</td><td>0.11</td><td>1135</td><td>-0.14 (-0.22, -0.06)</td><td>5.00E-04*</td></tr><tr><td>Fluid Intelligence (CPM)</td><td>1156</td><td>-0.03 (-0.07, 0.01)</td><td>0.21</td><td>1167</td><td>-0.03 (-0.08, 0.02)</td><td>0.21</td><td>1001</td><td>0.08 (-0.3, 0.47)</td><td>0.67</td><td>1147</td><td>0.06 (-0.01, 0.12)</td><td>0.08</td></tr><tr><td>Internalizing Behaviors (CBCL)</td><td>1156</td><td>0.03 (-0.03, 0.08)</td><td>0.33</td><td>1166</td><td>0.06 (-0.01, 0.13)</td><td>0.094</td><td>1002</td><td>0.4 (-0.12, 0.93)</td><td>0.13</td><td>1146</td><td>-0.09 (-0.18, 0)</td><td>0.053</td></tr><tr><td>Externalizing behaviors (CBCL)</td><td>1156</td><td>0.06 (0, 0.11)</td><td>0.032</td><td>1166</td><td>0.09 (0.02, 0.16)</td><td>0.01</td><td>1002</td><td>0.45 (-0.06, 0.97)</td><td>0.083</td><td>1146</td><td>-0.01 (-0.1, 0.08)</td><td>0.82</td></tr><tr><td>Lung Function (FEV1)</td><td>911</td><td>-0.01 (-0.07, 0.05)</td><td>0.75</td><td>921</td><td>-0.07 (-0.15, 0.01)</td><td>0.085</td><td>795</td><td>0.47 (-0.06, 1.01)</td><td>0.08</td><td>907</td><td>0.07 (-0.03, 0.17)</td><td>0.16</td></tr><tr><td>Puberty onset<sup>†</sup></td><td>537</td><td>0.92 (0.76, 1.11)</td><td>0.36</td><td>542</td><td>1.25 (0.99, 1.57)</td><td>0.058</td><td>511</td><td>1.84 (0.41, 8.44)</td><td>0.42</td><td>542</td><td>1.41 (1.01, 1.97)</td><td>0.046</td></tr></tbody></table>

Bold indicates p<0.05 and *indicates FDR <5%.

*Odds ratio provided for puberty onset only.

†Not available in the Lithuanian KANC cohort.

‡Only assessed in children over 8 years old.
immunometabolic Δ age was associated after FDR correction with better working memory (p=0.0036) and reduced inattentiveness (p=5 e-4).

Associations at the nominal significance (p<0.05) level was observed for increases in height z-score with transcriptome Δ age (p=0.014), shorter telomere length with increased fat mass % (p=0.009), and DNA methylation Δ age with greater inattentiveness (p=0.03). Both shorter telomere length and DNA methylation Δ age were associated with greater externalizing behaviors (p=0.032 and p=0.01, respectively). Among a smaller subset of children (Table 1) aged over 8 years, we observed a nominally significant association between immunometabolic Δ age and odds of onset of puberty (Odds Ratio: 1.41, 95% CI: 1.01, 1.97, p=0.046).

No significant associations with lung function were observed, but like the patterns of associations observed with cognitive and behavioral outcomes, there was a trend for a negative association with DNA methylation Δ age (p=0.085) and a positive association with immunometabolic Δ age (p=0.16).

## Sensitivity analysis

In sensitivity analysis. We first stratified by sex and generally observed similar associations among boys and girls, apart from the following differences (Figure 4—figure supplement 1): Associations between shorter telomere length and BMI z-score and adiposity were stronger among boys. For DNA methylation Δ age, associations with poorer externalizing and internalizing behaviours were only apparent among boys. For transcriptome Δ age, stronger associations among boys were observed with BMI z-score, adiposity, and poorer externalizing and internalizing behaviors. Conversely, we observed an association between transcriptome Δ age and reduced inattentiveness among girls only. Immunometabolic Δ age was more strongly associated with reduced inattentiveness among girls and also associated with greater odds of puberty onset among girls only.

Second, we additionally adjusted our models by estimating cell counts to determine the influence of cell composition on associations with developmental outcomes (Figure 4—figure supplement 2B, Figure 4—figure supplement 3B, Figure 4—figure supplement 4B and Figure 4—figure supplement 5B; Supplementary file 7). Associations were generally little changed: For DNA methylation Δ age, associations were attenuated with adiposity and growth outcomes although all remained FDR significant and the association with externalizing behavior was slightly attenuated. For transcriptome, Δ age associations with adiposity and growth outcomes and lung function increased slightly and the association with greater lung function became nominally significant.

Third, we assessed the effects of further adjustment for health risk factors (family affluence and social capital, birthweight, maternal active smoking, and child passive smoking) since health risk factors could be independently associated with both biological age and developmental outcomes (Figure 4—figure supplement 2C, Figure 4—figure supplement 3C, Figure 4—figure supplement 4C and Figure 4—figure supplement 5C; Supplementary file 7). Associations were generally little changed, expect for an attenuation of the association between telomere length and externalizing behavior, while conversely the association between DNA methylation Δ age and externalizing behavior was slightly strengthened.

Finally, we stratified by cohort study centre, observing generally consistent directions of effect among developmental outcomes of at least at nominal significance (p<0.05) in the main pooled analysis (Figure 4—figure supplement 6).

## Discussion

In a large sample of European children, we have analyzed two established and two candidate measures of biological age, derived from molecular features at different levels of biological organization, in relation to developmental outcomes and health risk factors. We assessed two established biological age markers, telomere length and DNA methylation age, and derived two new measures, transcriptome age, and immunometabolic age. Despite finding only null to weak correlations between the measures, we found all measures to be positively associated with greater BMI and adiposity, and both DNA methylation Δ age and immunometabolic Δ age were associated with taller height. While immunometabolic Δ age was associated with greater cognitive maturity including greater working memory and attentiveness, conversely DNA methylation Δ age was nominally associated with greater
inattentiveness and both DNA methylation Δ age and shorter telomere length were nominally associ-
ated with poorer externalizing behaviors.

BMI has consistently been associated with accelerated aging in adults across a diverse range of biological age markers (Peters et al., 2015; Jansen et al., 2021; Nevalainen et al., 2017; Fiorito et al., 2019) underlining the integral role of metabolism in aging. Indeed, a recent large study of Dutch adults found BMI to be the only health risk factor tested associated with accelerated aging across five biological age clocks, including telomere length, DNA methylation, transcriptome, proteome, and metabolomic age markers (Jansen et al., 2021). Here, we show that the link between BMI and multiple dimensions of accelerated aging is also apparent in children. Energy and nutrient intake influence all Ageing Hallmarks and multiple lines of evidence link increased adiposity to shorter lifespan (López-Otín et al., 2016). These effects appear to be partially mediated through evolutionarily conserved nutrient-sensing systems such as the mTOR signaling pathway, which promote anti-aging cellular repair mechanisms, at the expense of growth and metabolism, in response to lower nutrient availability (López-Otín et al., 2016). Furthermore, excess adiposity increases generalized inflammation and oxidative stress (Suzuki et al., 2003; Minamino et al., 2009), which may have direct effects on age markers, particularly telomere length and DNA methylation age acceleration.

The observed associations between greater height with biological age may indicate developmental maturity. Height is generally considered reflective of a beneficial early-life environment (Davey Smith et al., 2000), however, evidence for an association with lifespan is mixed (Davey Smith et al., 2000; Tanisawa et al., 2018), with a recent meta-analysis suggesting a U-shaped relationship with all-cause mortality (Li et al., 2021). Greater comparative height at 10 years was also inversely associated with longevity in a recent large-scale Medelian randomization study (Huang et al., 2021). Furthermore, there is some evidence that the link between height and longevity may be mediated through the insulin-like growth factor-1 signaling pathway (Tanisawa et al., 2018; He et al., 2014). The associations may also be interpreted as greater rates of growth and anabolism exerting greater 'wear and tear' on cellular structures. Two other studies have also observed an association between height and DNA methylation age acceleration in children (Simpkin et al., 2017; Suarez et al., 2018).

Despite similarities in associations with growth and adiposity measures, patterns of association across cognitive and behavioral domains varied across biological age markers, underlying the view of biological aging as a multi-faceted process. Immunometabolic Δ age was associated with greater cognitive maturity, fitting the life-course model of greater accumulation of biological capital during the build-up phase of development. Immunometabolic Δ age may be considered a phenotypic summary measure of metabolic and immune system maturity, and these cognitive developmental associations suggest that it may also be generalizable to the overall developmental stage. On the other hand, DNA methylation Δ age was related to relative immaturity in attentiveness and externalizing behavior. A previous Finish study of children aged between 11 and 13 years also reported associations between DNA methylation age acceleration and behavioral problems (Suarez et al., 2018). Similarly, shorter telomere length was associated here with greater externalizing behaviors, although not with any cognitive domains. Four other studies have examined the link between shorter telomere length and externalizing behaviors (Costa et al., 2015; Daoust et al., 2023; Wojcicki et al., 2015; Kroenke et al., 2011), with all, except one (Kroenke et al., 2011), also reporting an association.

Overall patterns of associations between risk factors and biological age measures also suggest the detrimental nature of accelerated aging in children assessed through telomere length and DNA methylation Δ age, and the potentially beneficial nature of advanced immunometabolic Δ age. Both prenatal maternal active smoking and child passive smoking were associated with DNA methylation Δ age, while greater birthweight was associated with immunometabolic Δ age. We examined maternal education level, family affluence, and social capital which broadly represent the three forms of interconvertible capital (cultural, economic, and social) proposed by Bourdieu, 1986. It has been theorized that biological capital represents a fourth type of human capital, and that the conversion across these forms of capital underlies inequalities in aging trajectories (Vineis and Kelly-Irving, 2019). Nominally significant associations between higher family affluence with longer telomere length and high social capital with a younger DNA methylation age indicate that age acceleration assessed through these measures does not represent an accumulation of biological capital. Generally, directions of effect for immunometabolic Δ age were in the opposite direction which may suggest it represents greater biological capital.
Girls were found to have longer telomere lengths than boys. Women have been consistently found to have longer telomere lengths (Gardner et al., 2014) although the few generally smaller studies in children have been inconsistent (Wojcicki et al., 2015; Buxton et al., 2011; Ly et al., 2019; Okuda et al., 2002). No other biological age markers were associated with sex, which contrasts with the study of Jansen et al., 2021 in adults which reported accelerated biological age in men across all measures tested except for proteomic age. Indeed, the phenomenon of accelerated DNA methylation age in men is well-established (Horvath and Raj, 2018), consistent with lower life-expectancies for men. Although it is not known if these biological age differences are due to biological mechanisms or greater prevalence of disease risk factors among men, our data in children before the divergence of risk factor prevalence could indicate a biological mechanism for telomere sex differences and a risk factor-mediated mechanism for other biological age markers. Interestingly, we observed differences in associations between biological age and development between boys and girls, with some consistency across markers: Both shorter telomere length and transcriptome Δ age and were more strongly associated with adiposity in boys, DNA methylation and transcriptome Δ age showed stronger associations among boys with poorer behavior, while in girls both transcriptome and immunometabolic Δ age showed stronger associations with improved attentiveness. Given observed sexual dimorphism in both developmental rates (León et al., 2014) and biological age measures through a variety of proposed mechanisms (Hägg and Jylhävä, 2021), it may be unsurprising that relationship between biological age and development also differs between the sexes.

Furthermore, we observed that immunometabolic Δ age was associated with greater odds of puberty onset, driven by effects observed among girls only. We did not observe any further significant associations with the onset of puberty, however, the sample size in the subset of children was small compared to the other developmental measures. There was also suggestive evidence for associations between DNA methylation Δ age with the onset of puberty with associations close to the nominal significance threshold. Three previous studies have reported associations between DNA methylation age acceleration and puberty onset and stage (Binder et al., 2018; Simpkin et al., 2017; Suarez et al., 2018), and one study has reported associations between shorter telomere length and puberty onset (Koss et al., 2020). However, the directions of effect for telomere length in our study were in the opposite direction. While earlier age at puberty is representative of more advanced physical maturation, it has been associated with metabolic diseases in later life, including cancers (Cancer", 2012) and all-cause mortality (Charalampopoulos et al., 2014).

We found transcriptome data to be highly accurate in predicting chronological age, including in a test set of children assessed six months later, demonstrating that gene expression tracks closely with age in children, even over this relatively short period. We analyzed biological pathways and processes enriched among transcript clusters contributing to the transcriptome clock, observing the integral role of ribosome and ribosome biogenesis pathways, central to protein synthesis, and biological processes including the immunity-related processes leukocyte migration and activation, and cell movement, activation, and secretion. Strikingly, gene expression in adults is similarly characterized by downregulation of ribosomal genes and enrichment of expression in immune-related genes (Peters et al., 2015). This indicates that similar to DNA methylation changes (Horvath and Raj, 2018), there is some overlap in gene expression related to both development in children and aging in adults. Although formal testing of enrichment of genes contributing to the transcriptome clock presented here among age-associated genes in adults showed enrichment at only borderline statistical significance, the transcriptome clock predictors are an underrepresentation of the full profile of gene expression associated with age in children, due to the sparsity enforced during the variable selection training process.

Despite associations with growth and adiposity measures, transcriptome age generally showed weaker associations with other developmental outcomes than for the other biological age markers. While this in part can be attributed to the slightly smaller sample size for children with transcriptome data, it is also likely due to the high accuracy in predicting the chronological age of the transcriptome clock, resulting in lower variation in the portion of transcriptome age that is not explained by chronological age, further reducing statistical power. This makes it challenging to judge the relevance of transcriptome age, if any, to developmental endpoints, which may be mixed since a non-significant direction of effects were observed with both maturity in attention and lung function, yet relative immaturity in behavior. In fact, training clocks using chronological age, which while providing an accessible route to understanding molecular changes associated with age, does pose limitations
generally for inference regarding biological aging. Particularly for high-dimensional data such as DNA methylation, it has been shown that it is possible to predict chronological age near-perfectly (Zhang et al., 2019), thereby limiting information on biological age and its variation. For this reason, newer epigenetic clocks have included clinical and mortality data, to improve clinical relevance and sensitivity to risk factors (Lu et al., 2019; Belsky et al., 2020), which should be considered in future studies developing clocks in children.

Other limitations include the cross-sectional design of the main analysis, which limits inference regarding the directionality of associations and allows the possibility of age-associated environmental factors to influence the clock development. Furthermore, there were differences in age by study centre. Children from the EDEN cohort study centre were generally older, which likely introduced a degree of cohort bias into the age modeling. For this reason, we adjusted all associations by study centre and additionally assessed age correlation within each study centre. Although cohorts were recruited from the general population, certain ethnicities or socio-economically disadvantaged groups may have been under-represented, limiting generalizability somewhat. A bias towards the over-representation of White ethnic groups is an issue generally with the development of biological clocks, which means associations observed with ethnicity should be interpreted cautiously. While the DNA methylation and transcriptome data were representative of the full genome, our coverage of the metabolome and proteome was limited to targeted assays. For analysis of age at the gene expression and protein/metabolite levels we developed new clocks, and these remain to be validated as true biological age indicators, through testing association with age-related disease and mortality in adults. While biological age markers exist for these data types, they are not appropriate for applying in our dataset, since the same model predictors (i.e. metabolites and proteins) were not included in the assays used here and/or the markers were not trained in pediatric populations, which will drastically reduce their accuracy as molecular changes within in childhood cannot be assumed to follow the same relation with age as in adulthood: For instance, DNA methylation shows a logarithmic dependence with age during childhood and a linear dependence in adulthood (Horvath, 2013).

However, the strengths of this study include the large population sample, drawing from six countries from around Europe, increasing generalizability, and the integration of rich molecular data and a broad range of developmental outcomes into a single systematic analysis. Although our age range was somewhat limited, missing the infancy and adolescent periods, the age range covered a key childhood period, where energy expenditure (an indicator of the level of overall physiology) has entered a period of steady increase following more rapid increases during infancy and before stabilization during adolescence (Pontzer et al., 2021).

In conclusion, in this large Pan-European study we have found that four indicators of biological age, representing complementary molecular processes, were all associated with BMI after controlling for chronological age, indicating that adiposity is an important correlate of accelerated biological aging in children. We developed a highly accurate 'transcriptome age' clock although it was found to be relatively insensitive to other development phenotypes. We found that immunometabolic Δ age was associated with cognitive maturity fitting a buildup of the biological capital model of aging in children, while shorter telomere length and DNA methylation Δ age was associated with greater behavioral problems suggesting a 'wear and tear' model of aging in children. Our findings contribute to the interpretation and understanding of biological age measures in children, crucial for clinical and epidemiological research into early life risk factors for adverse aging trajectories. Future long-term studies should investigate associations between age acceleration in children and adults to further test the antagonistic pleiotropy hypothesis.

## Acknowledgements

The authors are grateful to all the participating families in the six countries who took part in this study and to all the field workers for their dedication and efficiency. ISGlobal acknowledges support from the Spanish Ministry of Science and Innovation through the 'Centro de Excelencia Severo Ochoa 2019-2023' Program (CEX 2018-000806 S), and support from the Generalitat de Catalunya through the CERCA Programme. The CRG/UPF Proteomics Unit is part of the Spanish Infrastructure for Omics Technologies (ICTS OmicsTech) and it is supported by 'Secretaria d'Universitats i Recerca del Departament d'Economia i Coneixement de la Generalitat de Catalunya' (2017SGR595). We also acknowledge the support of the Spanish Ministry of Science and Innovation to the EMBL partnership, the
Centro de Excelencia Severo Ochoa, and the CERCA Programme/Generalitat de Catalunya. Some figures were created with BioRender.com.

## Additional information

### Funding

<table><thead><tr><th>Funder</th><th>Grant reference number</th><th>Author</th></tr></thead><tbody><tr><td>UK Research and Innovation</td><td>MR/S03532X/1</td><td>Oliver Robinson</td></tr><tr><td>European Commission</td><td>308333 874583</td><td>Oliver Robinson<br>ChungHo E Lau<br>Sandra Andrusaityte<br>Eva Borras<br>Paula de Prado-Bert<br>Lida Chatzi<br>Hector C Keun<br>Regina Grazuleviciene<br>Kristine B Gutzkow<br>Lea Maitre<br>Dries S Martens<br>Eduard Sabido<br>Valérie Siroux<br>Jose Urquiza<br>Marina Vafeiadi<br>John Wright<br>Tim S Nawrot<br>Mariona Bustamante</td></tr><tr><td>Secretaria d'Universitats i Recerca del Departament d'Economia i Coneixement de la Generalitat de Catalunya</td><td>2017SGR595</td><td>Jose Urquiza<br>Lea Maitre<br>Mariona Bustamante<br>Martine Vrijheid</td></tr><tr><td>Departament de Salut de la Generalitat de Catalunya</td><td>SLT017/20/000119</td><td>Jose Urquiza<br>Lea Maitre<br>Mariona Bustamante<br>Martine Vrijheid</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

Oliver Robinson, Conceptualization, Formal analysis, Funding acquisition, Writing - original draft, Project administration, Writing - review and editing; ChungHo E Lau, Formal analysis, Methodology, Writing - review and editing; Sungyeon Joo, Formal analysis, Writing - review and editing; Sandra Andrusaityte, Lea Maitre, Resources, Data curation, Investigation, Writing - review and editing; Eva Borras, Resources, Formal analysis, Writing - review and editing; Paula de Prado-Bert, Valérie Siroux, Methodology, Writing - review and editing; Lida Chatzi, Resources, Funding acquisition, Investigation, Writing - review and editing; Hector C Keun, Eduard Sabido, Resources, Data curation, Writing - review and editing; Regina Grazuleviciene, Resources, Data curation, Funding acquisition, Investigation; Kristine B Gutzkow, John Wright, Resources, Data curation, Funding acquisition, Investigation, Writing - review and editing; Dries S Martens, Resources, Investigation, Methodology, Writing - review and editing; Jose Urquiza, Data curation, Writing - review and editing; Marina Vafeiadi, Resources, Investigation, Writing - review and editing; Tim S Nawrot, Resources, Methodology, Writing - review and editing; Mariona Bustamante, Data curation, Formal analysis, Methodology, Writing - review and editing; Martine Vrijheid, Conceptualization, Resources, Data curation, Funding acquisition, Investigation, Methodology, Writing - review and editing

### Author ORCIDs

Oliver Robinson http://orcid.org/0000-0002-4735-0468
Kristine B Gutzkow http://orcid.org/0000-0002-6716-5921
Eduard Sabido http://orcid.org/0000-0001-6506-7714

## Ethics

Human subjects: Prior to the start of HELIX, all six cohorts had undergone the required evaluation by national ethics committees and obtained all the required permissions for their cohort recruitment and follow-up visits. Each cohort also confirmed that relevant informed consent and approval were in place for secondary use of data from pre-existing data. The work in HELIX was covered by new ethical approvals in each country and at enrolment in the new follow-up, participants were asked to sign a new informed consent form. Additionally, the current study was approved by the Imperial College Research Ethics Committee (Reference: 19IC5567).

## Decision letter and Author response

Decision letter https://doi.org/10.7554/eLife.85104.sa1

Author response https://doi.org/10.7554/eLife.85104.sa2

## Additional files

### Supplementary files

* MDAR checklist

* Source code 1. R script for all data analyses.

* Supplementary file 1. Number of samples and features before and after the quality control process.

* Supplementary file 2. Proportion of covariates missing for each biological age marker.

* Supplementary file 3. Immunometabolic age clock coefficients.

* Supplementary file 4. Transcriptome age clock coefficients.

* Supplementary file 5. Overrepresentation analysis in ConsesuspathDB against KEGG and REACTOME pathways, of all transcripts contributing to the transcriptome clock.

* Supplementary file 6. Overrepresentation analysis in ConsesuspathDB against Gene Ontology (GO) biological process terms, of all transcripts contributing to the transcriptome clock.

* Supplementary file 7. Associations between biological age measures and developmental measures, in main analysis (model 1) and sensitivity analyses (models 2-4).

### Data availability

Due to data protection regulations in each participating country and participant data use agreements, human subject data used in this project cannot be freely shared. The raw data supporting the current study are available on request subject to ethical and legislative review. The "HELIX Data External Data Request Procedures" are available with the data inventory in this website: http://www.projecthelix.eu/data-inventory. The document describes who can apply to the data and how, the timings for approval and the conditions to data access and publication. Researchers who have an interest in using data from this project for reproducibility or in using data held in general in the HELIX data warehouse for research purposes can apply for access to data. Interested researchers should fill in the application protocol found in ANNEX I at https://www.projecthelix.eu/files/helix_external_data_request_procedures_final.pdf and send this protocol to helixdata@isglobal.org. The applications are received by the HELIX Coordinator, and are processed and approved by the HELIX Project Executive Committee. All code used for data analysis has been provided as supplementary material. Deidentified dataset for generation of figures 1 and 2 has been provided as a supplementary dataset.

## References

* Achenbach TM. 2001. Manual for ASEBA School-Age Forms & Profiles University of Vermont, Research Center for Children, Youth & Families.

* Agier L, Basagaña X, Maitre L, Granum B, Bird PK, Casas M, Oftedal B, Wright J, Andrusaityte S, de Castro M, Cequier E, Chatzi L, Donaire-Gonzalez D, Grazuleviciene R, Haug LS, Sakhi AK, Leventakou V, McEachan R, Nieuwenhuijsen M, Petraviciene I, et al. 2019. Early-life exposome and lung function in children in Europe: an analysis of data from the longitudinal, population-based HELIX cohort. *The Lancet. Planetary Health* 3:e81–e92. DOI: https://doi.org/10.1016/S2542-5196(19)30010-5, PMID: 30737192
Belsky DW, Moffitt TE, Cohen AA, Corcoran DL, Levine ME, Prinz JA, Schaefer J, Sugden K, Williams B, Poulton R, Caspi A. 2018. Eleven telomere, epigenetic clock, and biomarker-composite quantifications of biological aging: Do they measure the same thing. *American Journal of Epidemiology* **187**:1220–1230. DOI: https://doi.org/10.1093/aje/kwx346, PMID: 29149257

Belsky DW, Caspi A, Arseneault L, Baccarelli A, Corcoran DL, Gao X, Hannon E, Harrington HL, Rasmussen LJ, Houts R, Huffman K, Kraus WE, Kwon D, Mill J, Pieper CF, Prinz JA, Poulton R, Schwartz J, Sugden K, Vokonas P, et al. 2020. Quantification of the pace of biological aging in humans through a blood test, the DunedinPoAm DNA methylation algorithm. *eLife* **9**:e54870. DOI: https://doi.org/10.7554/eLife.54870, PMID: 32367804

Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: a practical and powerful approach to multiple testing. *Journal of the Royal Statistical Society* **57**:289–300. DOI: https://doi.org/10.1111/j.2517-6161.1995.tb02031.x

Binder AM, Corvalan C, Mericq V, Pereira A, Santos JL, Horvath S, Shepherd J, Michels KB. 2018. Faster ticking rate of the epigenetic clock is associated with faster pubertal development in girls. *Epigenetics* **13**:85–94. DOI: https://doi.org/10.1080/15592294.2017.1414127, PMID: 29235933

Bourdieu P. 1986. The forms of capital. *Handbook of Theory and Research for the Sociology of Education* Westport, CT: Greenwood. p. 241–258.

Boyce W, Torsheim T, Currie C, Zambon A. 2006. The family affluence scale as a measure of national wealth: validation of an adolescent self-report measure. *Social Indicators Research* **78**:473–487. DOI: https://doi.org/10.1007/s11205-005-1607-6

Buckberry S, Bent SJ, Bianco-Miotto T, Roberts CT. 2014. massiR: a method for predicting the sex of samples in gene expression Microarray Datasets. *Bioinformatics* **30**:2084–2085. DOI: https://doi.org/10.1093/bioinformatics/btu161, PMID: 24659105

Buxton JL, Walters RG, Visvikis-Siest S, Meyre D, Froguel P, Blakemore AIF. 2011. Childhood obesity is associated with shorter leukocyte telomere length. *The Journal of Clinical Endocrinology and Metabolism* **96**:1500–1505. DOI: https://doi.org/10.1210/jc.2010-2924, PMID: 21349907

Cancer* C. 2012. Menarche, menopause, and breast cancer risk: individual participant meta-analysis, including 118 964 women with breast cancer from 117 epidemiological studies. *The Lancet Oncology* **13**:1141–1151. DOI: https://doi.org/10.1016/S1470-2045(12)70425-4

Cawthon RM. 2009. Telomere length measurement by a novel monochrome multiplex quantitative PCR method. *Nucleic Acids Research* **37**:e21. DOI: https://doi.org/10.1093/nar/gkn1027, PMID: 19129229

Charalampopoulos D, McLoughlin A, Elks CE, Ong KK. 2014. Age at menarche and risks of all-cause and cardiovascular death: a systematic review and meta-analysis. *American Journal of Epidemiology* **180**:29–40. DOI: https://doi.org/10.1093/aje/kwu113, PMID: 24920784

Chatzi L, Leventakou V, Vafeiadi M, Koutra K, Roumeliotaki T, Chalkiadaki G, Karachaliou M, Daraki V, Kyriklaki A, Kampouri M, Fthenou E, Sarri K, Vassilaki M, Fasoulaki M, Bitsios P, Koutis A, Stephanou EG, Kogevinas M. 2017. Cohort profile: the mother-child cohort in Crete, Greece (Rhea study). *International Journal of Epidemiology* **46**:1392–1393k. DOI: https://doi.org/10.1093/ije/dyx084, PMID: 29040580

Chen Y, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, Gallinger S, Hudson TJ, Weksberg R. 2013. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 Microarray. *Epigenetics* **8**:203–209. DOI: https://doi.org/10.4161/epi.23470, PMID: 23314698

Chen BH, Marioni RE, Colicino E, Peters MJ, Ward-Caviness CK, Tsai PC, Roetker NS, Just AC, Demerath EW, Guan W, Bressler J, Fornage M, Studenski S, Vandiver AR, Moore AZ, Tanaka T, Kiel DP, Liang L, Vokonas P, Schwartz J, et al. 2016. DNA methylation-based measures of biological age: meta-analysis predicting time to death. *Aging* **8**:1844–1865. DOI: https://doi.org/10.18632/aging.101020, PMID: 27690265

Clasey JL, Bradley KD, Bradley JW, Long DE, Griffith JR. 2011. A new BIA equation estimating the body composition of young children. *Obesity* **19**:1813–1817. DOI: https://doi.org/10.1038/oby.2011.158, PMID: 21681223

Coimbra BM, Carvalho CM, Moretti PN, Mello MF, Belangero SI. 2017. Stress-related telomere length in children: A systematic review. *Journal of Psychiatric Research* **92**:47–54. DOI: https://doi.org/10.1016/j.jpsychires.2017.03.023, PMID: 28407508

Costa DdeS, Rosa DVF, Barros AGA, Romano-Silva MA, Malloy-Diniz LF, Mattos P, de Miranda DM. 2015. Telomere length is highly inherited and associated with hyperactivity-Impulsivity in children with attention deficit/hyperactivity disorder. *Frontiers in Molecular Neuroscience* **8**:28. DOI: https://doi.org/10.3389/fnmol.2015.00028, PMID: 26217174

Daoust AR, Thakur A, Kotelnikova Y, Kleiber ML, Singh SM, Hayden EP. 2023. Associations between children's telomere length, parental Intrusiveness, and the development of early externalizing behaviors. *Child Psychiatry and Human Development* **54**:672–682. DOI: https://doi.org/10.1007/s10578-021-01279-3, PMID: 34727279

Davey Smith G, Hart C, Upton M, Hole D, Gillis C, Watt G, Hawthorne V. 2000. Height and risk of death among men and women: Aetiological implications of associations with cardiorespiratory disease and cancer mortality. *Journal of Epidemiology and Community Health* **54**:97–103. DOI: https://doi.org/10.1136/jech.54.2.97, PMID: 10715741

de Onis M, Onyango AW, Borghi E, Siyam A, Nishida C, Siekmann J. 2007. Development of a WHO growth reference for school-aged children and adolescents. *Bulletin of the World Health Organization* **85**:660–667. DOI: https://doi.org/10.2471/blt.07.043497, PMID: 18026621

de Prado-Bert P, Ruiz-Arenas C, Vives-Usano M, Andrusaityte S, Cadiou S, Carracedo Á, Casas M, Chatzi L, Dadvand P, González JR, Grazuleviciene R, Gutzkow KB, Haug LS, Hernandez-Ferrer C, Keun HC, Lepeule J,
Maitre L, McEachan R, Nieuwenhuijsen MJ, Pelegrí D, et al. 2021. The early-life exposome and epigenetic age acceleration in children. *Environment International* **155**:106683. DOI: https://doi.org/10.1016/j.envint.2021.106683, PMID: 34144479

Dieterle F, Ross A, Schlotterbeck G, Senn H. 2006. Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. application in 1H NMR Metabonomics. *Analytical Chemistry* **78**:4281–4290. DOI: https://doi.org/10.1021/ac051632c, PMID: 16808434

Emery Thompson M. 2022. Evolutionary approaches in aging research. *Cold Spring Harbor Perspectives in Medicine* **12**:a041195. DOI: https://doi.org/10.1101/cshperspect.a041195, PMID: 36041879

Ferrucci L, Gonzalez-Freire M, Fabbri E, Simonsick E, Tanaka T, Moore Z, Salimi S, Sierra F, de Cabo R. 2020. Measuring biological aging in humans: A quest. *Aging Cell* **19**:e13080. DOI: https://doi.org/10.1111/acel.13080, PMID: 31833194

Fiorito G, McCrory C, Robinson O, Carmeli C, Ochoa-Rosales C, Zhang Y, Colicino E, Dugué PA, Artaud F, McKay GJ, Jeong A, Mishra PP, Nøst TH, Krogh V, Panico S, Sacerdote C, Tumino R, Palli D, Matullo G, Guarrera S, et al. 2019. Socioeconomic position, lifestyle habits and biomarkers of epigenetic aging: a multi-cohort analysis. *Aging* **11**:2045–2070. DOI: https://doi.org/10.18632/aging.101900, PMID: 31009935

Forns J, Esnaola M, López-Vicente M, Suades-González E, Alvarez-Pedrerol M, Julvez J, Grellier J, Sebastián-Gallés N, Sunyer J. 2014. The N-back test and the attentional network task as measures of child neuropsychological development in Epidemiological studies. *Neuropsychology* **28**:519–529. DOI: https://doi.org/10.1037/neu0000085, PMID: 24819069

Fortin J-P, Labbe A, Lemire M, Zanke BW, Hudson TJ, Fertig EJ, Greenwood CM, Hansen KD. 2014. Functional normalization of 450K methylation array data improves replication in large cancer studies. *Genome Biology* **15**:503. DOI: https://doi.org/10.1186/s13059-014-0503-2, PMID: 25599564

Gardner M, Bann D, Wiley L, Cooper R, Hardy R, Nitsch D, Martin-Ruiz C, Shiels P, Sayer AA, Barbieri M, Bekaert S, Bischoff C, Brooks-Wilson A, Chen W, Cooper C, Christensen K, De Meyer T, Deary I, Der G, Diez Roux A, et al. 2014. Gender and telomere length: systematic review and meta-analysis. *Experimental Gerontology* **51**:15–27. DOI: https://doi.org/10.1016/j.exger.2013.12.004, PMID: 24365661

Giallourou N, Fardus-Reid F, Panic G, Veselkov K, McCormick BJJ, Olortegui MP, Ahmed T, Mduma E, Yori PP, Mahfuz M, Svensen E, Ahmed MMM, Colston JM, Kosek MN, Swann JR. 2020. Metabolic maturation in the first 2 years of life in resource-constrained settings and its association with postnatal growths. *Science Advances* **6**:eaay5969. DOI: https://doi.org/10.1126/sciadv.aay5969, PMID: 32284996

Grazuleviciene R, Danileviciute A, Nadisauskiene R, Vencloviene J. 2009. Maternal smoking, GSTM1 and GSTT1 polymorphism and susceptibility to adverse pregnancy outcomes. *International Journal of Environmental Research and Public Health* **6**:1282–1297. DOI: https://doi.org/10.3390/ijerph6031282, PMID: 19440446

Guxens M, Ballester F, Espada M, Fernández MF, Grimalt JO, Ibarluzea J, Olea N, Rebagliato M, Tardón A, Torrent M, Vioque J, Vrijheid M, Sunyer J, INMA Project. 2012. Cohort profile: the INMA--INFancia Y Medio Ambiente--(Environment and childhood) project. *International Journal of Epidemiology* **41**:930–940. DOI: https://doi.org/10.1093/ije/dyr054, PMID: 21471022

Hägg S, Jylhävä J. 2021. Sex differences in biological aging with a focus on human studies. *eLife* **10**:e63425. DOI: https://doi.org/10.7554/eLife.63425, PMID: 33982659

Hansen KD AM. 2012. Illuminahumanmethylation450Kmanifest: Annotation for Illumina's 450K methylation Arrays. 3.17. Bioconductor. https://bioconductor.org/packages/release/data/annotation/html/IlluminaHuma
nMethylation450kmanifest.html

He Q, Morris BJ, Grove JS, Petrovitch H, Ross W, Masaki KH, Rodriguez B, Chen R, Donlon TA, Willcox DC, Willcox BJ, Forloni G. 2014. Shorter men live longer: Association of height with longevity and Foxo3 genotype in American men of Japanese ancestry. *PLOS ONE* **9**:e94385. DOI: https://doi.org/10.1371/journal.pone.0094385

Hertel J, Friedrich N, Wittfeld K, Pietzner M, Budde K, Van der Auwera S, Lohmann T, Teumer A, Völzke H, Nauck M, Grabe HJ. 2016. Measuring biological age via metabonomics: the metabolic age score. *Journal of Proteome Research* **15**:400–410. DOI: https://doi.org/10.1021/acs.jproteome.5b00561, PMID: 26652958

Heude B, Forhan A, Slama R, Douhaud L, Bedel S, Saurel-Cubizolles MJ, Hankard R, Thiebaugeorges O, De Agostini M, Annesi-Maesano I, Kaminski M, Charles MA, EDEN mother-child cohort study group. 2016. Cohort profile: the EDEN mother-child cohort on the prenatal and early postnatal determinants of child health and development. *International Journal of Epidemiology* **45**:353–363. DOI: https://doi.org/10.1093/ije/dyv151, PMID: 26283636

Horvath S. 2013. DNA methylation age of human tissues and cell types. *Genome Biology* **14**:R115. DOI: https://doi.org/10.1186/gb-2013-14-10-r115, PMID: 24138928

Horvath S, Oshima J, Martin GM, Lu AT, Quach A, Cohen H, Felton S, Matsuyama M, Lowe D, Kabacik S, Wilson JG, Reiner AP, Maierhofer A, Flunkert J, Aviv A, Hou L, Baccarelli AA, Li Y, Stewart JD, Whitsel EA, et al. 2018. Epigenetic clock for skin and blood cells applied to Hutchinson Gilford Progeria syndrome and ex Vivo studies. *Aging* **10**:1758–1775. DOI: https://doi.org/10.18632/aging.101508, PMID: 30048243

Horvath S, Raj K. 2018. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. *Nature Reviews. Genetics* **19**:371–384. DOI: https://doi.org/10.1038/s41576-018-0004-3, PMID: 29643443

Hotamisligil GS. 2017. Inflammation, metaflammation and immunometabolic disorders. *Nature* **542**:177–185. DOI: https://doi.org/10.1038/nature21363, PMID: 28179656

Huang S-Y, Yang Y-X, Chen S-D, Li H-Q, Zhang X-Q, Kuo K, Tan L, Feng L, Dong Q, Zhang C, Yu J-T. 2021. Investigating causal relationships between exposome and human longevity: a Mendelian randomization analysis. *BMC Medicine* **19**:150. DOI: https://doi.org/10.1186/s12916-021-02030-4, PMID: 34281550
Jansen R, Han LK, Verhoeven JE, Aberg KA, van den Oord EC, Milaneschi Y, Penninx BW. 2021. An integrative study of five biological clocks in somatic and mental health. *eLife* **10**:10. DOI: https://doi.org/10.7554/eLife.59479

Johnson WE, Li C, Rabinovic A. 2007. Adjusting batch effects in microarray expression data using empirical Bayes methods. *Biostatistics* **8**:118–127. DOI: https://doi.org/10.1093/biostatistics/kxj037, PMID: 16632515

Julvez J, López-Vicente M, Warembourg C, Maitre L, Philippat C, Gütkow KB, Guxens M, Evandt J, Andrusaityte S, Burgaleta M, Casas M, Chatzi L, de Castro M, Donaire-González D, Gražulevičienė R, Hernandez-Ferrer C, Heude B, Mceachan R, Mon-Williams M, Nieuwenhuijsen M, et al. 2021. Early life multiple exposures and child cognitive function: A multi-centric birth cohort study in six European countries. *Environmental Pollution* **284**:117404. DOI: https://doi.org/10.1016/j.envpol.2021.117404, PMID: 34077897

Jylhävä J, Pedersen NL, Hägg S. 2017. Biological age predictors. *EBioMedicine* **21**:29–36. DOI: https://doi.org/10.1016/j.ebiom.2017.03.046, PMID: 28396265

Koss KJ, Schneper LM, Brooks-Gunn J, McLanahan S, Mitchell C, Notterman DA. 2020. Early puberty and telomere length in preadolescent girls and mothers. *The Journal of Pediatrics* **222**:193–199.. DOI: https://doi.org/10.1016/j.jpeds.2020.02.075, PMID: 32586523

Kroenke CH, Epel E, Adler N, Bush NR, Obradovic J, Lin J, Blackburn E, Stamperdahl JL, Boyce WT. 2011. Autonomic and adrenocortical reactivity and buccal cell telomere length in kindergarten children. *Psychosomatic Medicine* **73**:533–540. DOI: https://doi.org/10.1097/PSY.0b013e318229acfc, PMID: 21873585

Kuh D. 2007. A life course approach to healthy aging, frailty, and capability. *The Journals of Gerontology. Series A, Biological Sciences and Medical Sciences* **62**:717–721. DOI: https://doi.org/10.1093/gerona/62.7.717, PMID: 17634317

Lau C-HE, Siskos AP, Maitre L, Robinson O, Athersuch TJ, Want EJ, Urquiza J, Casas M, Vafeiadi M, Roumeliotaki T, McEachan RRC, Azad R, Haug LS, Meltzer HM, Andrusaityte S, Petraviciene I, Grazuleviciene R, Thomsen C, Wright J, Slama R, et al. 2018. Determinants of the urinary and serum metabolome in children from six European populations. *BMC Medicine* **16**:202. DOI: https://doi.org/10.1186/s12916-018-1190-8, PMID: 30404627

Lehallier B, Gate D, Schaum N, Nanasi T, Lee SE, Yousef H, Moran Losada P, Berdnik D, Keller A, Verghese J, Sathyan S, Franceschi C, Milman S, Barzilai N, Wyss-Coray T. 2019. Undulating changes in human plasma proteome profiles across the lifespan. *Nature Medicine* **25**:1843–1850. DOI: https://doi.org/10.1038/s41591-019-0673-2

Lehne B, Drong AW, Loh M, Zhang W, Scott WR, Tan S-T, Afzal U, Scott J, Jarvelin M-R, Elliott P, McCarthy MI, Kooner JS, Chambers JC. 2015. A coherent approach for analysis of the Illumina Humanmethylation450 Beadchip improves data quality and performance in epigenome-wide association studies. *Genome Biology* **16**:37. DOI: https://doi.org/10.1186/s13059-015-0600-x, PMID: 25853392

León I, Cimadevilla JM, Tascón L. 2014. Developmental gender differences in children in a virtual spatial memory task. *Neuropsychology* **28**:485–495. DOI: https://doi.org/10.1037/neu0000054, PMID: 24588700

Li Q, Liu Y, Sun X, Li H, Cheng C, Liu L, Liu F, Zhou Q, Guo C, Tian G, Qie R, Han M, Huang S, Li L, Wang B, Zhao Y, Ren Y, Zhang M, Hu D, Wu J, et al. 2021. Dose-response association between adult height and all-cause mortality: a systematic review and meta-analysis of cohort studies. *European Journal of Public Health* **31**:652–658. DOI: https://doi.org/10.1093/eurpub/ckaa213

López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. 2013. The hallmarks of aging. *Cell* **153**:1194–1217. DOI: https://doi.org/10.1016/j.cell.2013.05.039, PMID: 23746838

López-Otín C, Galluzzi L, Freije JMP, Madeo F, Kroemer G. 2016. Metabolic control of longevity. *Cell* **166**:802–821. DOI: https://doi.org/10.1016/j.cell.2016.07.031, PMID: 27518560

Lu AT, Quach A, Wilson JG, Reiner AP, Aviv A, Raj K, Hou L, Baccarelli AA, Li Y, Stewart JD, Whitsel EA, Assimes TL, Ferrucci L, Horvath S. 2019. DNA methylation GrimAge strongly predicts lifespan and healthspan. *Aging* **11**:303–327. DOI: https://doi.org/10.18632/aging.101684, PMID: 30669119

Ly K, Walker C, Berry S, Snell R, Marks E, Thayer Z, Atatao-Carr P, Morton S. 2019. Telomere length in early childhood is associated with sex and ethnicity. *Scientific Reports* **9**:10359. DOI: https://doi.org/10.1038/s41598-019-46338-x, PMID: 31316106

Macdonald-Dunlop E, Taba N, Klarić L, Frkatović A, Walker R, Hayward C, Esko T, Haley C, Fischer K, Wilson JF, Joshi PK. 2022. A catalogue of omics biological ageing clocks reveals substantial commonality and associations with disease risk. *Aging* **14**:623–659. DOI: https://doi.org/10.18632/aging.203847, PMID: 35073279

Magnus P, Birke C, Vejrup K, Haugan A, Alsaker E, Daltveit AK, Handal M, Haugen M, Høiseth G, Knudsen GP, Paltiel L, Schreuder P, Tambs K, Vold L, Stoltenberg C. 2016. Cohort profile update: the Norwegian mother and child cohort study (MoBa). *International Journal of Epidemiology* **45**:382–388. DOI: https://doi.org/10.1093/ije/dyw029, PMID: 27063603

Maitre L, de Bont J, Casas M, Robinson O, Aasvang GM, Agier L, Andrušaitytė S, Ballester F, Basagaña X, Borràs E, Brochot C, Bustamante M, Carracedo A, de Castro M, Dedele A, Donaire-Gonzalez D, Estivill X, Evandt J, Fossati S, Giorgis-Allemand L, et al. 2018. Human Early Life Exposome (HELIX) study: a European population-based exposome cohort. *BMJ Open* **8**:e021311. DOI: https://doi.org/10.1136/bmjopen-2017-021311, PMID: 30206078

Marioni RE, Shah S, McRae AF, Ritchie SJ, Muniz-Terrera G, Harris SE, Gibson J, Redmond P, Cox SR, Pattie A, Corley J, Taylor A, Murphy L, Starr JM, Horvath S, Visscher PM, Wray NR, Deary IJ. 2015. The epigenetic clock is correlated with physical and cognitive fitness in the Lothian Birth Cohort 1936. *International Journal of Epidemiology* **44**:1388–1396. DOI: https://doi.org/10.1093/ije/dyu277, PMID: 25617346
Martens DS, Plusquin M, Gyselaers W, De Vivo I, Nawrot TS. 2016. Maternal pre-pregnancy body mass index and newborn telomere length. BMC Medicine **14**:148. DOI: https://doi.org/10.1186/s12916-016-0689-0, PMID: 27751173

Martens DS, Van Der Stukken C, Derom C, Thiery E, Bijnens EM, Nawrot TS. 2021. Newborn Telomere length predicts later life telomere length: tracking telomere length from birth to child- and adulthood. EBioMedicine **63**:103164. DOI: https://doi.org/10.1016/j.ebiom.2020.103164, PMID: 33422989

McEwen LM, O'Donnell KJ, McGill MG, Edgar RD, Jones MJ, Maclsaac JL, Lin DTS, Ramadori K, Morin A, Gladish N, Garg E, Unternaehrer E, Pokhvisneva I, Karnani N, Kee MZL, Klengel T, Adler NE, Barr RG, Letourneau N, Giesbrecht GF, et al. 2020. The PedBE clock accurately estimates DNA methylation age in pediatric buccal cells. PNAS **117**:23329-23335. DOI: https://doi.org/10.1073/pnas.1820843116, PMID: 31611402

Minamino T, Orimo M, Shimizu I, Kunieda T, Yokoyama M, Ito T, Nojima A, Nabetani A, Oike Y, Matsubara H, Ishikawa F, Komuro I. 2009. A crucial role for adipose tissue P53 in the regulation of insulin resistance. Nature Medicine **15**:1082-1087. DOI: https://doi.org/10.1038/nm.2014, PMID: 19718037

Moskalev A. 2019. Biomarkers of human aging. Moskalev A (Ed). *Introduction, in Biomarkers of Human Aging* Cham: Springer Cham. DOI: https://doi.org/10.1007/978-3-030-24970-0

Nadarajah S, Kotz S. 2006. The Exponentiated type distributions. Acta Applicandae Mathematicae **92**:97-111. DOI: https://doi.org/10.1007/s10440-006-9055-0

Nevalainen T, Kananen L, Marttila S, Jylhävä J, Mononen N, Kähönen M, Raitakari OT, Hervonen A, Jylhä M, Lehtimäki T, Hurme M. 2017. Obesity accelerates epigenetic aging in middle-aged but not in elderly individuals. *Clinical Epigenetics* **9**:20. DOI: https://doi.org/10.1186/s13148-016-0301-7, PMID: 28289477

Okuda K, Bardeguez A, Gardner JP, Rodriguez P, Ganesh V, Kimura M, Skurnick J, Awad G, Aviv A. 2002. Telomere length in the newborn. *Pediatric Research* **52**:377-381. DOI: https://doi.org/10.1203/00006450-200209000-00012, PMID: 12193671

Pelegí-Sisó D, de Prado P, Ronkainen J, Bustamante M, González JR. 2021. Methylclock: a bioconductor package to estimate DNA methylation age. *Bioinformatics* **37**:1759-1760. DOI: https://doi.org/10.1093/bioinformatics/btaa825, PMID: 32960939

Perna L, Zhang Y, Mons U, Holleczek B, Saum KU, Brenner H. 2016. Epigenetic age acceleration predicts cancer, cardiovascular, and all-cause mortality in a German case cohort. *Clinical Epigenetics* **8**:64. DOI: https://doi.org/10.1186/s13148-016-0228-z, PMID: 27274774

Peters MJ, Joehanes R, Pilling LC, Schurmann C, Conneely KN, Powell J, Reinmaa E, Sutphin GL, Zhernakova A, Schramm K, Wilson YA, Kobes S, Tukiainen T, Ramos YF, Göring HHH, Fornage M, Liu Y, Gharib SA, Stranger BE, De Jager PL, et al. 2015. The transcriptional landscape of age in human peripheral blood. *Nature Communications* **6**:8570. DOI: https://doi.org/10.1038/ncomms9570, PMID: 26490707

Petersen AC, Crockett L, Richards M, Boxer A. 1988. A self-report measure of pubertal status: reliability, validity, and initial norms. *Journal of Youth and Adolescence* **17**:117-133. DOI: https://doi.org/10.1007/BF01537962, PMID: 24277579

Pontzer H, Yamada Y, Sagayama H, Ainslie PN, Andersen LF, Anderson LJ, Arab L, Baddou I, Bedu-Addo K, Blaak EE, Blanc S, Bonomi AG, Bouten CVC, Bovet P, Buchowski MS, Butte NF, Camps SG, Close GL, Cooper JA, Cooper R, et al. 2021. Daily energy expenditure through the human life course. *Science* **373**:808-812. DOI: https://doi.org/10.1126/science.abe5017, PMID: 34385400

Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, Culver BH, Enright PL, Hankinson JL, Ip MSM, Zheng J, Stocks J, ERS Global Lung Function Initiative. 2012. Multi-ethnic reference values for spirometry for the 3-95-yr age range: the global lung function 2012 equations. *The European Respiratory Journal* **40**:1324-1343. DOI: https://doi.org/10.1183/09031936.00080312, PMID: 22743675

Robinson O, Chadeau Hyam M, Karaman I, Climaco Pinto R, Ala-Korpela M, Handakas E, Fiorito G, Gao H, Heard A, Jarvelin MR, Lewis M, Pazoki R, Polidoro S, Tzoulaki I, Wielscher M, Elliott P, Vineis P. 2020. Determinants of accelerated metabolomic and epigenetic aging in a UK cohort. *Aging Cell* **19**:e13149. DOI: https://doi.org/10.1111/acel.13149, PMID: 32363781

Rutledge J, Oh H, Wyss-Coray T. 2022. Measuring biological age using omics data. *Nature Reviews. Genetics* **23**:715-727. DOI: https://doi.org/10.1038/s41576-022-00511-7, PMID: 35715611

Sayed N, Huang Y, Nguyen K, Krejciova-Rajaniemi Z, Grawe AP, Gao T, Tibshirani R, Hastie T, Alpert A, Cui L, Kuznetsova T, Rosenberg-Hasson Y, Ostan R, Monti D, Lehallier B, Shen-Orr SS, Maecker HT, Dekker CL, Wyss-Coray T, Franceschi C, et al. 2021. An inflammatory aging clock (iAge) based on deep learning tracks multimorbidity, immunosenescence, frailty and cardiovascular aging. *Nature Aging* **1**:598-615. DOI: https://doi.org/10.1038/s43587-021-00082-y

Serra-Majem L, Ribas L, Ngo J, Ortega RM, García A, Pérez-Rodrigo C, Aranceta J. 2004. Food, youth and the Mediterranean diet in Spain. Development of KIDMED, Mediterranean Diet Quality Index in children and adolescents. *Public Health Nutrition* **7**:931-935. DOI: https://doi.org/10.1079/PHN2004556

Sierra F. 2016. The emergence of geroscience as an interdisciplinary approach to the enhancement of health span and life span. *Cold Spring Harbor Perspectives in Medicine* **6**:a025163. DOI: https://doi.org/10.1101/cshperspect.a025163, PMID: 26931460

Simpkin AJ, Howe LD, Tilling K, Gaunt TR, Lyttleton O, McArdle WL, Ring SM, Horvath S, Smith GD, Relton CL. 2017. The epigenetic clock and physical development during childhood and adolescence: longitudinal analysis from a UK birth cohort. *International Journal of Epidemiology* **46**:549-558. DOI: https://doi.org/10.1093/ije/dyw307, PMID: 28089957
Siskos AP, Jain P, Römisch-Margl W, Bennett M, Achaintre D, Asad Y, Marney L, Richardson L, Koulman A, Griffin JL, Raynaud F, Scalbert A, Adamski J, Prehn C, Keun HC. 2017. Interlaboratory reproducibility of a targeted metabolomics platform for analysis of human serum and plasma. *Analytical Chemistry* **89**:656–665. DOI: https://doi.org/10.1021/acs.analchem.6b02930, PMID: 27959516

Suarez A, Lahti J, Czamara D, Lahti-Pulkkinen M, Girchenko P, Andersson S, Strandberg TE, Reynolds RM, Kajantie E, Binder EB, Raikkonen K. 2018. The epigenetic clock and pubertal, Neuroendocrine, psychiatric, and cognitive outcomes in adolescents. *Clinical Epigenetics* **10**:96. DOI: https://doi.org/10.1186/s13148-018-0528-6, PMID: 30021623

Suzuki K, Ito Y, Ochiai J, Kusuhara Y, Hashimoto S, Tokudome S, Kojima M, Wakai K, Toyoshima H, Tamakoshi K, Watanabe Y, Hayakawa N, Maruta M, Watanabe M, Kato K, Ohta Y, Tamakoshi A, JACC Study Group. 2003. Relationship between obesity and serum markers of oxidative stress and inflammation in Japanese. *Asian Pacific Journal of Cancer Prevention* **4**:259–266 PMID: 14507248.

Tanisawa K, Hirose N, Arai Y, Shimokata H, Yamada Y, Kawai H, Kojima M, Obuchi S, Hirano H, Suzuki H, Fujiwara Y, Taniguchi Y, Shinkai S, Ihara K, Sugaya M, Higuchi M, Arai T, Mori S, Sawabe M, Sato N, et al. 2018. Inverse association between height-increasing alleles and extreme longevity in Japanese women. *The Journals of Gerontology* **73**:588–595. DOI: https://doi.org/10.1093/gerona/glx155

van den Akker EB, Trompet S, Barkey Wolf JJH, Beekman M, Suchiman HED, Deelen J, Asselbergs FW, Boersma E, Cats D, Elders PM, Geleijnse JM, Ikram MA, Kloppenburg M, Mei H, Meulenbelt I, Mooijaart SP, Nelissen RGHH, Netea MG, Penninx BWJH, Slofstra M, et al. 2020. Metabolic age based on the BBMRI-NL (1) H-NMR Metabolomics repository as biomarker of age-related disease. *Circulation. Genomic and Precision Medicine* **13**:541–547. DOI: https://doi.org/10.1161/CIRCGEN.119.002610, PMID: 33079603

van Iterson M, Tobi EW, Slieker RC, den Hollander W, Luijk R, Slagboom PE, Heijmans BT. 2014. Methylaid: visual and interactive quality control of large Illumina 450K Datasets. *Bioinformatics* **30**:3435–3437. DOI: https://doi.org/10.1093/bioinformatics/btu566, PMID: 25147358

Vineis P, Kelly-Irving M. 2019. Biography and biological capital. *European Journal of Epidemiology* **34**:979–982. DOI: https://doi.org/10.1007/s10654-019-00539-w, PMID: 31342230

Vives-Usano M, Hernandez-Ferrer C, Maitre L, Ruiz-Arenas C, Andrusaityte S, Borràs E, Carracedo Á, Casas M, Chatzi L, Coen M, Estivill X, González JR, Grazuleviciene R, Gutzkow KB, Keun HC, Lau C-HE, Cadiou S, Lepeule J, Mason D, Quintela I, et al. 2020. In utero and childhood exposure to tobacco smoke and multi-layer molecular signatures in children. *BMC Medicine* **18**:243. DOI: https://doi.org/10.1186/s12916-020-01686-8, PMID: 32811491

Vrijheid M, Slama R, Robinson O, Chatzi L, Coen M, van den Hazel P, Thomsen C, Wright J, Athersuch TJ, Avellana N, Basagaña X, Brochot C, Bucchini L, Bustamante M, Carracedo A, Casas M, Estivill X, Fairley L, van Gent D, Gonzalez JR, et al. 2014. The Human Early-Life Exposome (HELIX): project rationale and design. *Environmental Health Perspectives* **122**:535–544. DOI: https://doi.org/10.1289/ehp.1307204, PMID: 24610234

Wojcicki JM, Heyman MB, Elwan D, Shiboski S, Lin J, Blackburn E, Epel E. 2015. Telomere length is associated with oppositional defiant behavior and maternal clinical depression in Latino preschool children. *Translational Psychiatry* **5**:e581. DOI: https://doi.org/10.1038/tp.2015.71, PMID: 26080316

Wright J, Small N, Raynor P, Tuffnell D, Bhopal R, Cameron N, Fairley L, Lawlor DA, Parslow R, Petherick ES, Pickett KE, Waiblinger D, West J, Born in Bradford Scientific Collaborators Group. 2013. Cohort profile: The Born in Bradford multi-ethnic family cohort study. *International Journal of Epidemiology* **42**:978–991. DOI: https://doi.org/10.1093/ije/dys112, PMID: 23064411

Zhang Q, Vallerga CL, Walker RM, Lin T, Henders AK, Montgomery GW, He J, Fan D, Fowdar J, Kennedy M, Pitcher T, Pearson J, Halliday G, Kwok JB, Hickie I, Lewis S, Anderson T, Silburn PA, Mellick GD, Harris SE, et al. 2019. Improved precision of epigenetic clock estimates across tissues and its implication for biological ageing. *Genome Medicine* **11**:54. DOI: https://doi.org/10.1186/s13073-019-0667-1, PMID: 31443728

Zou H, Hastie T. 2005. Regularization and variable selection via the elastic net. *Journal of the Royal Statistical Society Series B* **67**:301–320. DOI: https://doi.org/10.1111/j.1467-9868.2005.00503.x