Environment- and epigenome-wide association study of obesity in ‘Children of 1997’ birth cohort

  1. Jie Zhao  Is a corresponding author
  2. Bohan Fan
  3. Jian Huang
  4. Benjamin John Cowling
  5. Shiu Lun Ryan Au Yeung
  6. Andrea Baccarelli
  7. Gabriel M Leung
  8. C Mary Schooling
  1. School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong
  2. Singapore Institute for Clinical Sciences (SICS), Agency for Science, Technology and Research (A*STAR), Singapore
  3. Mailman School of Public Health, Columbia University, United States
  4. City University of New York, School of Public Health and Health policy, United States

Abstract

Background:

Increasing childhood obesity is a global issue requiring potentially local solutions to ensure it does not continue into adulthood. We systematically identified potentially modifiable targets of obesity at the onset and end of puberty in Hong Kong, the most economically developed major Chinese city.

Methods:

We conducted an environment-wide association study (EWAS) and an epigenome-wide association study of obesity to systematically assess associations with body mass index (BMI) and waist–hip ratio (WHR) in Hong Kong’s population-representative ‘Children of 1997’ birth cohort. Univariable linear regression was used to select exposures related to obesity at ~11.5 years (BMI and obesity risk n ≤ 7119, WHR n = 5691) and ~17.6 years (n = 3618) at Bonferroni-corrected significance, and multivariable regression to adjust for potential confounders followed by replicated multivariable regression (n = 308) and CpG by CpG analysis (n = 286) at ~23 years. Findings were compared with evidence from published randomized controlled trials (RCTs) and Mendelian randomization (MR) studies.

Results:

At ~11.5 and~17.6 years the EWAS identified 14 and 37 exposures associated with BMI, as well as 7 and 12 associated with WHR, respectively. Most exposures had directionally consistent associations at ~23 years. Maternal second-hand smoking, maternal weight, and birth weight were consistently associated with obesity. Diet (including dairy intake and artificially sweetened beverages), physical activity, snoring, binge eating, and earlier puberty were positively associated with BMI at ~17.6 years, while eating before sleep was inversely associated with BMI at ~17.6 years. Findings for birth weight, dairy intake, and binge eating are consistent with available evidence from RCTs or MR studies. We found 17 CpGs related to BMI and 17 to WHR.

Conclusions:

These novel insights into potentially modifiable factors associated with obesity at the outset and the end of puberty could, if causal, inform future interventions to improve population health in Hong Kong and similar Chinese settings.

Funding:

This study including the follow-up survey and epigenetics testing was supported by the Health and Medical Research Fund Research Fellowship, Food and Health Bureau, Hong Kong SAR Government (#04180097). The DNA extraction of the samples used for epigenetic testing was supported by CFS-HKU1.

Editor's evaluation

This study presents a valuable finding on the association of environmental and epigenomic factors with obesity in the adolescent population, using a specialized cohort. The evidence supporting the claims of the authors is solid, although this study is inevitably open to residual confounding because of the limitation of observational study. The work will be of interest to both clinicians and researchers working on obesity and related metabolic disorders.

https://doi.org/10.7554/eLife.82377.sa0

Introduction

With improving living standards and socioeconomic development, non-communicable chronic diseases pose a heavy burden on society in both developed and developing countries (Lozano et al., 2012). Obesity is a well-established risk factor for multiple chronic diseases, including cardiovascular disease, diabetes, and cancer (Gallagher and LeRoith, 2015). According to the World Health Organization (WHO), obesity is defined as ‘abnormal or excessive fat accumulation that presents a risk to health’ (World Health Organization, 2000). Body mass index (BMI) is the most appropriate measure for overweight and obesity because the cut-offs account for age, sex, and ethnicity (World Health Organization, 2000). Obesity has increased substantially in many settings, including in Hong Kong. Obesity is complicated with multifactorial risk factors, such as socioeconomic position (SEP), mood disturbance, and genetic factors (Cardel et al., 2020), so there is a need to evaluate their roles in obesity comprehensively and systematically (Manrai et al., 2017). Moreover, given most studies on targets of obesity are conducted in a western setting (Cardel et al., 2020), a comprehensive assessment of modifiable factors of early life obesity in a non-western setting, with a different social structure, provides a valuable opportunity to identify novel exposures.

Environment-wide association studies (EWAS) enable us to assess a variety of exposures across the human environmental exposome in a high-throughput manner (Hall et al., 2013), similar to genome-wide association studies for genetic associations. Previous EWAS have been performed on outcomes, such as type 2 diabetes (Patel et al., 2010), cardiovascular disease (Zhuang et al., 2018), and childhood obesity in western settings (Vrijheid et al., 2020; Uche et al., 2020) but not in a Chinese setting. In the previous EWAS of childhood obesity in the US (6–17 years old), UK, and Europe (6–11 years old), second-hand smoking was related to higher childhood BMI, whilst some other exposures, such as vitamins, were not consistently associated with obesity in these settings (Vrijheid et al., 2020; Uche et al., 2020). Observational studies are open to confounding by SEP, thus assessing associations in a different social context can triangulate the evidence concerning early life obesity.

In addition to the environmental factors, it is increasingly realized that epigenetic factors, which are also modifiable, may play an important role in obesity (Huang et al., 2018a). DNA methylation, which refers to the addition of a methyl group to the 5′ position of a cytosine residue of the DNA, is the most frequently examined epigenetic modification (Fall et al., 2017). DNA methylation may modulate gene expression and thereby influence susceptibility to obesity or obesity-related chronic disease (Fall et al., 2017). Epigenome-wide association study provides an approach to identify the related epigenetic loci in a comprehensive way. For example, DNA methylation at cg06500161 was previously identified as related to obesity in the US (Huang et al., 2018a). Unlike genetic variants, DNA methylation is modifiable and may change in response to environmental factors or disease and therefore might be open to confounding by these factors (Fall et al., 2017; Fraga et al., 2005). As such, findings from western settings may not be generalizable to Chinese populations.

In this situation, Hong Kong, a non-western developed setting, can provide unique insights into health determinants. Most Chinese people in Hong Kong are first-, second-, or third-generation migrants from the neighbouring province of Guangdong in southern China. Dietary habits of people in Hong Kong are similar to those in southern China, although also influenced by western culture (Leung et al., 2003). Lifestyle in Hong Kong also differs from more commonly studied western populations on some important attributes, for example, active smoking among Chinese mothers was rare while maternal exposure to second-hand smoking during pregnancy was common before the smoking ban in public and workplaces was implemented in 2007 (Lee, 2016). Moreover, most current theories concerning the aetiology of and disparities in chronic diseases originate from observations in long-term developed populations of European descent. However, in Hong Kong the economic transition from pre- to post-industrial living conditions has occurred within one lifetime of the older people (Leung et al., 2017), whereas children today in Hong Kong represent the first generation of Chinese to grow up in a post-industrial Chinese setting, which is unrivalled anywhere in the world (Schooling et al., 2012). As such, a study in young Chinese people in Hong Kong, a different setting provides ‘a sentinel for populations currently experiencing very rapid economic development’ (Schooling et al., 2016), which may help identify whether these associations reflect SEP within a specific context or are biologically based as well as having the potential to identify any attributes relevant to the majority of the global population but not necessarily evident in more commonly studied western populations, such as maternal birthplace (Schooling et al., 2010). For example, in the ‘Children of 1997’ birth cohort, a large cohort in Hong Kong, the associations of sugar-sweetened beverages (Zhang et al., 2020), breastfeeding (Hui et al., 2018), milk consumption frequency (Lin et al., 2012), sleep duration (Wang et al., 2019), and parental smoking (Kwok et al., 2010) with childhood and adolescent obesity have been examined, with some important differences detected. The associations for breastfeeding in Hong Kong are much more similar to those seen in randomized controlled trials (RCTs) than those typically seen in western settings (Hui et al., 2018). To take advantage of this unique setting, we conducted an environment- and epigenome-wide association study, to identify further potential drivers of obesity. We focused on puberty because it is an important stage involving a re-orientation from childhood priorities to adulthood (Karlberg, 1989). Exposures associated with obesity at puberty may be important for health in later life (Richardson et al., 2020). Considering that exposures related to obesity at the outset and at the end of puberty may be different, and the associations of the same exposure with obesity may vary by age, we conducted the EWAS at the onset and the end of puberty.

Methods

Participants

The study takes advantage of the ‘Children of 1997’ birth cohort, a large (n = 8327) population-representative Chinese cohort in Hong Kong (Schooling et al., 2012). The participants were originally recruited shortly after birth in April and May 1997 at all of the 49 Maternal and Child Health Centers (MCHCs) in Hong Kong, which provide free check-ups and immunizations. The study included 88% of births in the relevant period. A self-administered questionnaire in Chinese was used at baseline to collect information on family, education, birth characteristics, infant feeding, and second-hand smoke exposure. The initial study was designed to provide a short-term assessment of the effects of second-hand smoking and included follow-up via the MCHCs until 18 months. In 2005, funded by the Health and Health Services Research Fund (HHSRF) and Health and Medical Research Fund (HMRF) we extended the information on this cohort via record linkage to include infant characteristics, serious morbidity, childhood obesity, pubertal development, history of migration, and SEP; with regular updates on subsequent growth obtained from the Student Health Service, including annual height and weight measurements from age ~6 years, and in this study, we used height and weight at age ~11.5 years. In 2007, with support from The University of Hong Kong University Research Committee Strategic Research Theme of Public Health, we instituted a program to re-establish and maintain direct contact with the cohort through direct mailing (newsletters, birthday cards, and seasonal cards) and the mass media (press conference and a full-length television documentary). We have since conducted three questionnaires/telephone surveys and an in-person Biobank Clinical follow-up in one visit at age ~17.6 years (Phase 1 in 2013–2016 included 3460 people with mean age 17.5 years and Phase 2 in the second half of 2017 included 158 people at mean age 19.5 years) with their blood samples stored. In 2020 (at age ~23 years), we conducted a follow-up survey to obtain updated information on anthropometric measurements. The number of participants in each age is shown in Figure 1—figure supplement 1.

Outcomes

At ~11.5 years, BMI was calculated from height and weight measurements records provided by the Student Health Service, the Department of Health. Waist–hip ratio (WHR) was calculated based on waist and hip circumference collected in Survey I conducted in 2008–2009. At ~17.6 years, in both phases of the Biobank Clinical follow-up, BMI was assessed by bio-electrical impedance analysis (BIA) with a Tanita segmental body composition monitor (Tanita BC-545, Tanita Co, Tokyo, Japan). Waist and hip circumference measurements were made using a tape twice following a standard protocol by trained technicians and nurses. Given BMI has a more accepted cut-off value than WHR for children, we classified obesity as BMI ≥20.89 kg/m2 for boys and BMI ≥21.20 kg/m2 for girls at age 11.5; and BMI ≥24.73 kg/m2 for men and BMI ≥24.85 kg/m2 for women at age 17.6 (Cole et al., 2000). In the follow-up survey at ~23 years, questionnaires were sent to 700 participants randomly selected from those with blood samples available and with BMI below the 25th centile or above the 75th centile. The questionnaires were accompanied by clear instructions on anthropometric measurement and a tape measure (the same as used in the Biobank Clinical follow-up). In total, 308 participants replied and provided their waist, hip, height, and body weight.

Assessment of DNA methylation

DNA methylation was conducted in 286 participants randomly selected from the 308 participants in the follow-up survey. DNA were extracted from buffy coat samples previously stored at −80° using EZI DNA blood kit (QIAGEN) with magnetic particle technology. DNA methylation was assessed using the Illumina Methylation EPIC Beadchip, which interrogated the methylation status of over 850,000 CpG sites. We conducted quality control using the ‘ewastools’ package (Heiss and Just, 2018), which included an evaluation of control metrics monitoring the various experimental steps, such as bisulfite conversion or staining and a sex check comparing actual sex to the records. After sample-level quality control, we excluded two samples that had a sex mismatch, so 286 samples (168 women and 118 men) were included in the analysis. We corrected for dye bias using RELIC (Xu et al., 2017), without normalization. At the probe level, we excluded non-CpG probes and probes located on the sex chromosomes; a total of 843,393 probes remained for analyses.

Exposures and categorization

Supplementary file 1 shows exposure categorization and data sources. After excluding exposures with missing values ≥50%, we included 123 exposures for BMI and 115 exposures for WHR at ~11.5 years, and 441 exposures for BMI and WHR at ~17.6 years using information from the original study, record linkage, the three surveys, that is Survey I (2008–2009), Survey II (2010–2012), and Survey III (2011–2012), as well as the Biobank Clinical follow-up. The exposures considered for obesity at ~11.5 years were classified into 12 categories, including baseline characteristics, SEP, family history, paternal information, maternal information, infant feeding and caring, diet (measured at Survey I at ~11.5 years old), health status (referring to physical health condition; details of the questions can be found in Supplementary file 1), parents’ health status, physical activity, lifestyle, and home facilities and pets. The exposures at 17.6 years were classified into 16 categories: baseline characteristics, SEP, family history, paternal information, maternal information, infant feeding and caring, diet (measured at the Biobank Clinical follow-up at ~17.6 years), children’s use of medications, children’s health status, parent’s health status, physical activity, home facilities and pets, moods and feelings, academic performance, sleep, and pubertal timing.

Statistical analysis

In the EWAS, similar to genome-wide association studies (Barrera-Gómez et al., 2017), first we used univariable linear regression to assess associations of each of the exposures with the measures of obesity at ages ~11.5 and ~17.6 years. We conducted the analysis in people with both exposure and outcome available, specifically, in up to 7119 participants for BMI at ~11.5 years, 5691 participants for WHR at ~11.5 years, and 3618 participants for BMI and WHR at ~17.6 years; their baseline characteristics at different ages are shown in Table 1. We only considered exposures reaching Bonferroni-corrected significance (e.g. p < 0.05/441 = 1.2 × 10−4 for obesity at ~17.6 years) to account for multiple testing (Curtin and Schulz, 1998). Second, we used multivariable linear regression controlling for potential confounders (sex, housing type at birth, household income at birth, maternal second-hand smoking during pregnancy, maternal age at birth, maternal education, maternal birth place, and the interaction of maternal education with maternal birthplace [Schooling et al., 2010]) at age ~11.5 and ~17.6 years, and excluded exposures that had over 50% of change-in-estimates ratios (Lee, 2014). Third, to assess whether the associations differed by age, we checked the associations for the selected exposures from the earlier age groups (~11.5 and ~17.6 years) in the follow-up survey (n = 308) at age ~23 years and compared the direction of associations with those at earlier age groups (~11.5 and ~17.6 years). Associations with consistent directions of associations in earlier age groups (~11.5 or ~17.6 years) with those at ~23 years suggest a consistent association by age. Additionally, to account for the time lag between the age at which physical measurements were taken and age at exposure collection, we included a time difference variable in the adjusted models for BMI. Finally, exposures that remained after controlling for confounders, were compared with the evidence from existing RCTs and Mendelian randomization (MR) studies, a study design which uses genetic variants as instrument and provides less confounded associations (Smith and Ebrahim, 2003).

Table 1
Baseline characteristics of Hong Kong’s ‘Children of 1997’ birth cohort included in the environment-wide association study (EWAS) of obesity at 11.5, 17.6, and 23 years.
CharacteristicsAt ~11.5 years(n ≤ 7119)At ~17.6 years(n = 3618)At ~23 years(n = 308)
Socio-demographics
SexBoys37321825123
Girls33871793185
Maternal birthplaceNot Hong Kong (migrant)27071509128
Hong Kong41852089179
Highest maternal educationGrade 9 or below28731421124
Grade 10–1130601575127
Grade 12 or above106460057
Household income per head in quintiles at recruitment1st quintile (HK$ 1746 ± 419)126059945
2nd quintile (HK$ 2853 ± 325)128864749
3rd quintile (HK$ 4365 ± 557)126964658
4th quintile (HK$ 6826 ± 883)124365465
5th quintile (HK$ 14,945 ± 15,620)124067467
Maternal age at delivery (mean [SD])30.2 (4.7)30.5 (4.6)30.8 (4.4)
Anthropometrics
Mean height in cm (SD)Boys166.4 (7.5)171.8 (6.0)173.1 (6.3)
Girls158.3 (5.9)159.5 (5.4)160.6 (5.2)
Mean weight in kg (SD)Boys54.6 (11.3)62.7 (11.8)67.3 (15.0)
Girls48.5 (8.0)52.6 (8.9)56.0 (30.5)
Mean BMI (SD)Boys19.2 (3.7)21.1 (3.7)22.4 (4.8)
Girls18.3 (3.1)20.6 (3.3)21.7 (11.5)
Mean WHR (SD)Boys0.84 (0.07)0.79 (0.07)0.87 (0.07)
Girls0.79 (0.06)0.75 (0.05)0.78 (0.07)

In the epigenome-wide association study, considering the heteroscedasticity in methylation beta-values (Du et al., 2010), we used robust linear regression models to assess the epigenome-wide association of each CpG with BMI and WHR at age ~23 years. We adjusted for age at blood draw for DNA methylation, age at follow-up survey, sex, cell type proportion, methylation assay batches, maternal second-hand smoking during pregnancy, maternal education, maternal birthplace, and household income at birth. The significance was considered as p < 1 × 10−6, genome-wide significance (5 × 10−8) was not used given the relatively small sample size in the epigenome-wide association study. To estimate genomic inflation, we used a Bayesian method that estimates inflation more accurately in epigenome-wide association studies based on the empirical null distributions (van Iterson et al., 2017), implemented using the R package ‘bacon’.

Ethical approval

This study complies with the Declaration of Helsinki. Since our participants are children, informed (non-written) consent for the original survey and subsequent record linkage was obtained from the parents, next of kin, caretakers, or guardians (informants) on behalf of the participants by the informant agreeing and subsequently completing the questionnaire at enrollment, this manner of obtaining consent was approved by The University of Hong Kong Medical Faculty Ethics Committee over 20 years ago. Informed written consent for subsequent Surveys and in-person follow-up was obtained from a parent or guardian, or at ages 18+ years from the participant. Ethical approval for this study, including the follow-up survey at ~23 years and comprehensive health-related analyses, was obtained from the University of Hong Kong-Hospital Authority Hong Kong West Cluster, Joint Institutional Review Board, Hong Kong Special Administrative Region, China (reference numbers: UW13-367; UW19-367).

Results

Figures 1 and 2 show the association of each exposure with BMI and WHR at ~11.5 and ~17.6 years. At ~11.5 years, 18 associations with BMI and 19 associations with WHR remained after Bonferroni correction (Figure 1). Of these 18 associations with BMI, 14 associations with BMI remained after controlling for confounders, 13 showed significant association with obesity risk at 11.5 years, and 11 exposures had concordant direction of associations with BMI at ~23 years (Table 2). Of the 19 associations with WHR at ~11.5 years, 7 exposures remained after controlling for confounders, and 6 had the same direction of association with WHR at ~23 years (Table 3). At ~17.6 years, 37 associations with BMI and 19 with WHR remained after Bonferroni correction (Figure 2). Of these 37 associations with BMI, all remained after controlling for confounders, 27 exposures were associated with obesity risk, and 32 exposures had the same direction of association with BMI at ~23 years (Table 4). Of the 19 associations with WHR at ~17.6 years, 12 remained after controlling for confounders, and 7 had the same directions of association with WHR at ~23 years (Table 5).

Figure 1 with 1 supplement see all
Associations of all exposures with body mass index (BMI) and waist–hip ratio (WHR) at age ~11.5 years in the univariable regression in Hong Kong’s ‘Children of 1997’ birth cohort.

(a) Associations of all exposures with body mass index (BMI) at age ~11.5 years in the univariable regression in Hong Kong’s ‘Children of 1997’ birth cohort. (b) Associations of all exposures with waist–hip ratio (WHR) at age ~11.5 years in the univariable regression in Hong Kong’s ‘Children of 1997’ birth cohort. DM, diabetes; edu, education; ASB, artificially sweetened beverages; PA, physical activity; SHS, second-hand smoking; freq, frequency. In total, we included 123 exposures for BMI at 11.5 years and 115 exposures for WHR at 11.5 years. The cut-off lines indicate Bonferroni-corrected p thresholds (p < 0.05/123 = 4.07 × 10−4 for BMI, p < 0.05/115 = 4.35 × 10−4 for WHR).

Associations of all exposures with body mass index (BMI) and waist–hip ratio (WHR) at age ~17.6 years in the univariable regression in 3618 participants of Hong Kong’s ‘Children of 1997’ birth cohort in the Biobank Clinical follow-up.

(a) Associations of all exposures with body mass index (BMI) at age ~17.6 years in the univariable regression in 3618 participants of Hong Kong’s ‘Children of 1997’ birth cohort in the Biobank Clinical follow-up. (b) Associations of all exposures with waist–hip ratio (WHR) at age ~17.6 years in the univariable regression in 3618 participants of Hong Kong’s ‘Children of 1997’ birth cohort in the Biobank Clinical follow-up. DM, diabetes; ASB, artificially sweetened beverages; PA, physical activity; SHS, second-hand smoking; freq, frequency. The cut-off lines indicate Bonferroni-corrected p thresholds (p < 0.05/441 exposures = 1.2 × 10−4 for BMI and WHR).

Table 2
Associations of selected exposures with body mass index (BMI) and obesity risk at ~11.5 years, and associations of those exposures with BMI at ~17.6 years and with BMI at ~23 years in participants of Hong Kong’s ‘Children of 1997’ birth cohort.
Group namesVariable descriptionWith BMI at ~11.5 years*With obesity risk at ~11.5 years*With BMI at ~17.6 yearsWith BMI at ~23 years
Betap valueOR (95% CI)Betap valueBetap value
Baseline characteristicsSex0.855.65E−222.06 (1.82, 2.32)0.231.99E−030.258.44E−01
Baseline characteristicsBirth weight0.405.40E−151.23 (1.15, 1.32)0.567.84E−060.613.78E−01
Baseline characteristicsSmall for gestational age: birth weight <10% by sex and gestational week distribution in singletons−0.689.93E−060.71 (0.56, 0.89)−0.425.65E−02−0.976.46E−01
SEPNo. of smokers at home0.152.83E−021.15 (1.06, 1.25)0.204.26E−021.073.06E−01
Paternal informationPaternal weight0.051.58E−181.03 (1.02, 1.04)0.041.73E−090.045.68E−01
Maternal informationMaternal weight0.091.84E−491.06 (1.05, 1.07)0.091.30E−280.185.23E−02
Maternal informationPostnatal smoking hygiene (i.e. second-hand smoking by timing [pre- and/or postnatal])0.107.79E−051.06 (1.02, 1.10)0.142.51E−040.363.67E−01
Maternal informationSecond-hand smoke by sources and timing0.082.69E−051.05 (1.02, 1.08)0.103.31E−040.234.54E−01
Maternal informationMother exposed to second-hand smoke during pregnancy0.154.07E−031.08 (1.01, 1.16)0.239.70E−050.325.85E−01
Family historyPaternal diabetes0.695.82E−031.31 (0.93, 1.86)0.312.10E−018.742.08E−04
Family historyMaternal diabetes1.452.47E−052.21 (1.42, 3.44)0.362.23E−01−0.897.67E−01
Family historyGestational diabetes0.682.28E−041.51 (1.18, 1.93)0.571.84E−01−0.839.37E−01
DietWater: frequency of consumption in the last week0.329.45E−061.20 (1.06, 1.36)0.371.89E−041.272.38E−01
LifestyleHaving meals: hours spent yesterday−0.402.06E−030.80 (0.65, 0.99)−0.421.65E−021.713.59E−01
  1. *

    Association of selected exposures with BMI and obesity risk at age ~11.5 years in up to 7119 participants after controlling for confounders.

  2. Association of selected exposures for BMI at age ~11.5 years with BMI at age ~17.6 years (n = 3618) and ~23 years (n = 308) after controlling for confounders.

Table 3
Associations of selected exposures with waist–hip ratio (WHR) at ~11.5 years, and associations of those exposures with WHR at ~17.6 years and with WHR at ~23 years in participants of Hong Kong’s ‘Children of 1997’ birth cohort.
Group namesVariable descriptionWith WHR at ~11.5 years*With WHR at ~17.6 yearsWith WHR at ~23 years
Betap valueBetap valueBetap value
Baseline characteristicsSex0.056.58E−1710.041.28E−880.0824.30E−18
SEPHousing type at birth (home ownership)−0.012.70E−02−0.0018.51E−01−0.014.20E−01
SEPMaternal education level−0.014.90E−05−0.0042.03E−01−0.0019.02E−01
SEPMother’s birthplace−0.011.18E−01−0.012.76E−01−0.0067.97E−01
Maternal informationMaternal weight0.007.98E−100.0011.13E−060.0016.58E−02
Family historyMaternal diabetes0.032.30E−050.011.27E−01−0.0177.82E−01
LifestyleHaving meals: hours spent yesterday−0.018.15E−09−4.73E−031.30E−01−0.0028.11E−01
  1. *

    Association of selected exposures with WHR at age ~11.5 years in up to 7119 participants after controlling for confounders.

  2. Association of selected exposures for WHR at age ~11.5 years with BMI at age ~17.6 years (n = 3618) and ~23 years (n = 308) after controlling for confounders.

Table 4
Associations of selected exposures with body mass index (BMI) and obesity risk at ~17.6 years, and associations of those exposures with BMI at ~23 years in participants from Hong Kong’s ‘Children of 1997’ birth cohort.
Group namesVariable descriptionWith BMI at ~17.6 years*With obesity risk at ~17.6 years*With BMI at ~23 years
Betap valueOR (95% CI)Betap value
Baseline characteristicsSex0.567.84E−061.72 (1.38, 2.14)0.258.44E−01
Baseline characteristicsBirth weight0.231.99E−031.1 (0.97, 1.25)0.613.78E−01
Baseline characteristicsTwin−1.771.05E−030.17 (0.02, 1.27)−3.735.28E−01
Paternal informationPaternal weight0.041.73E−091.03 (1.02, 1.04)0.045.68E−01
Maternal informationMaternal weight0.091.30E−281.05 (1.03, 1.06)0.185.23E−02
Maternal informationMother exposed to second-hand smoke during pregnancy0.239.70E−051.12 (0.99, 1.27)0.029.78E−01
Maternal informationPostnatal smoking hygiene (i.e. second-hand smoking by timing [pre- and/or postnatal])0.096.88E−021.72 (1.38, 2.14)0.781.57E−01
DietChocolate: frequency of consumption in the past week−0.349.02E−060.73 (0.63, 0.84)0.296.92E−01
DietChocolate: no. of servings/time−0.191.09E−040.83 (0.73, 0.95)0.088.69E−01
DietSweets: frequency of consumption in the past week−0.254.18E−040.83 (0.73, 0.95)−0.148.25E−01
DietSweets: no. of servings/time−0.175.42E−050.85 (0.77, 0.94)−0.176.85E−01
DietChinese tea: frequency of consumption in the past week0.271.69E−051.25 (1.13, 1.38)0.296.30E−01
DietChinese tea: no. of cup/time0.172.08E−071.14 (1.08, 1.2)0.185.74E−01
DietChinese herbal tea: frequency of consumption in the past week0.323.16E−041.18 (1.02, 1.36)1.111.82E−01
DietChinese herbal tea: no. of cup/time0.129.17E−041.13 (1.05, 1.22)0.265.01E−01
DietCoffee: frequency of consumption in the past week0.302.90E−051.17 (1.04, 1.32)0.386.05E−01
DietWater: no. of cup/time0.251.07E−091.13 (1.05, 1.22)0.511.97E−01
DietEnergy drinks: no. of cup/time0.141.21E−041.13 (0.99, 1.3)0.664.84E−01
DietWhen consume milk, how often reduced-fat/skim milk0.602.56E−240.7 (0.63, 0.77)1.645.35E−03
DietWhen consume soy milk/flavored drinks, how often lower sugar0.485.14E−110.73 (0.64, 0.82)0.782.55E−01
DietWhen consume soft drinks, how often diet/artificially sweetened drinks0.501.75E−100.7 (0.62, 0.8)0.623.61E−01
DietEating: 1 hr before sleeping−0.762.85E−060.65 (0.47, 0.9)−1.045.04E−01
Children’s HealthParents’ questionnaire: do you consider your child now to be: very thin/underweight; a little bit thin/underweight; about the right weight; a little bit fat/overweight; very fat/overweight; don’t know2.510.00E+007.45 (6.09, 9.12)2.134.12E−04
Children’s HealthSnoring in the past 4 weeks (child)0.769.72E−080.51 (0.4, 0.65)1.303.18E−01
Children’s HealthGrowth problem (child)0.367.41E−081.5 (1.37, 1.64)−0.029.76E−01
Children’s HealthDiabetes (child)0.324.62E−041.29 (1.15, 1.45)0.019.94E−01
Physical activityBall games: no. of times you done past week0.224.04E−031.03 (0.92, 1.17)1.299.89E−02
Physical activityGymnasium: no. of times you done past week0.443.10E−051.14 (0.97, 1.34)0.258.18E−01
Physical activityWhat did you normally do at lunch in the past 7 days0.345.50E−051.22 (1.06, 1.39)−0.228.12E−01
Physical activityHow many nights did you do sports in the past 7 days0.161.13E−020.98 (0.88, 1.08)0.653.43E−01
Physical activityTuesday: how often you did physical activity0.182.15E−031.01 (0.92, 1.11)−0.049.48E−01
Physical activityThursday: how often you did physical activity0.212.86E−041.08 (0.98, 1.2)1.509.65E−03
Moods and feelingsHas he/she gone on eating binges3.225.23E−155.19 (3.18, 8.46)3.373.06E−01
SleepCough or snore loudly0.671.62E−131.43 (1.25, 1.63)1.152.15E−01
SleepHave bad dreams−0.271.87E−030.77 (0.64, 0.91)−0.992.31E−01
Pubertal timingAge of onset of pubic hair−0.562.41E−100.78 (0.67, 0.92)−0.851.24E−03
Pubertal timingAge of menarche (girls)−0.812.46E−270.54 (0.45, 0.64)−1.112.18E−01
  1. *

    Association of selected exposures with BMI and obesity risk at age ~17.6 years in 3618 participants after controlling for confounders.

  2. Association of selected exposures for BMI at age ~17.6 years with BMI at ~23 years (n = 308) after controlling for confounders.

Table 5
Associations of selected exposures with waist–hip ratio (WHR) at ~17.6 years, and associations of those exposures with WHR at ~23 years in participants from Hong Kong’s ‘Children of 1997’ birth cohort.
Group namesVariable descriptionWith WHR at ~17.6 years*With WHR at ~23 years
Betap valueBetap value
Baseline characteristicsSex0.0451.28E−880.0835.71E−18
Maternal informationMaternal weight0.0011.13E−060.0019.38E−02
SEPHousing type at birth (home ownership)−0.0062.70E−02−0.0066.10E−01
DietBiscuits: frequency of consumption in the past week−0.0032.60E−020.0009.36E−01
DietChocolate: frequency of consumption in the past week−0.0041.86E−03−0.0035.20E−01
DietChocolate: no. of servings/time−0.0035.03E−050.0009.11E−01
DietSweets: frequency of consumption in the past week−0.0031.07E−020.0009.80E−01
DietWhen consume soft drinks, how often diet/artificially sweetened drinks0.0051.04E−04–0.0089.36E−02
Children’s HealthParents’ questionnaire: do you consider your child now to be: very thin/underweight; a little bit thin/underweight; about the right weight; a little bit fat/overweight; very fat/overweight; don’t know0.0213.55E−590.0153.02E−04
Children’s HealthWhether experienced musculoskeletal pain in past 4 weeks0.0061.57E−02−0.0018.94E−01
SleepCough or snore loudly0.0072.30E−060.0134.76E−02
SleepSleep disturbance (child)0.0051.93E−030.0036.28E−01
  1. *

    Association of selected exposures with WHR at age ~17.6 years in 3618 participants after controlling for confounders.

  2. Association of selected exposures for WHR at age ~17.6 years with WHR at ~23 years (n = 308) after controlling for confounders.

Specifically, for obesity at ~11.5 years, we found sex (being male), higher birth weight, maternal second-hand smoking, higher parental weight, family history of diabetes, gestational diabetes, and more water consumption were associated with higher BMI, while being small for gestational age and spending more time having meals were associated with lower BMI at ~11.5 years (Table 2). Except for paternal diabetes which was not significant, the rest of exposures all showed consistent associations with obesity risk (Table 2). However, the associations for family history of diabetes and time spent on meals showed inconsistent directions of associations for BMI ~23 years (Table 2). Regarding WHR, the associations were generally consistent with those seen for BMI and showed consistent directions of associations with those at ~17.6 and ~23 years, except for maternal diabetes (Table 3).

For obesity at ~17.6 years, in addition to some shared factors, including sex, birth weight, parental weight, and maternal second-hand smoking, we found some aspects of diet (i.e. more artificially sweetened beverage [ASB], lower-sugar soy milk, reduced-fat/skim milk, Chinese herbal tea, Chinese tea, energy drinks, coffee, and fish consumptions), physical activity, health status (i.e. diabetes, growth problem and snoring), earlier puberty and binge eating associated with higher BMI at ~17.6 years (Table 4). Being a twin, sweets consumption, chocolate consumption, eating before sleep, and having bad dreams were associated with lower BMI at ~17.6 years (Table 4). Regarding the association with obesity risk at 17.6 years, birth weight, being a twin, maternal second-hand smoking, physical activity, and energy drinks intake were not significant but the direction of association was consistent (Table 4). In addition, most exposures had the same direction of association at ~23 years. However, the association of chocolate consumption with BMI at ~23 years was in the other direction. Regarding WHR, the associations were generally consistent with those seen for BMI. Sex, drinking ASB, children’s health status, and coughing or snoring during sleep were also related to WHR at ~17.6 and ~23 years (Table 5). The associations of selected exposures with BMI at ~11.5 and ~17.6 years were similar after adjusting for the time difference between age at anthropometric measurements and age of exposure collection.

Regarding the comparison with RCTs and MR studies, we found RCTs on dark chocolate consumption (Kord-Varkaneh et al., 2019), water consumption promotion (Muckelbauer et al., 2009), and physical activity (Bleich et al., 2018), and MR studies related to drinking coffee (Nordestgaard et al., 2015), dairy intake (Huang et al., 2018b), binge eating (Reed et al., 2017), physical activity (Carrasquilla et al., 2022), snoring (Campos et al., 2020), puberty (Gill et al., 2018; Bell et al., 2018), birth weight (Zanetti et al., 2018), and maternal obesity (Bond et al., 2022; Richmond et al., 2017; Table 6). The available evidence from both RCTs and MR studies show that more physical activity lowers BMI (Bleich et al., 2018), and MR studies also suggest that dairy intake (Yang et al., 2017; Huang et al., 2018b) and binge eating (Reed et al., 2017) are associated with higher BMI.

Table 6
Evidence from published systematic reviews, randomized controlled trials (RCTs) and Mendelian randomization (MR) studies regarding the role of exposures selected in our environment-wide association study (EWAS) in obesity.
ExposurePublished RCTs studiesPublished MR studies
Water consumptionIntervention on promoting water consumption showed no effect on body mass index (BMI; Muckelbauer et al., 2009).NA
Reduced-fat/skim milk consumptionNAMR study on different types of milk consumption (i.e. reduced-fat, skimmed, reduced-sugar milk) among children is lacking. Nevertheless, two MR studies suggested higher dairy intake was associated with higher adult’s BMI (Huang et al., 2018b; Yang et al., 2017).
Coffee consumptionNAGenetically predicted more coffee intake was not associated with obesity, BMI, or waist circumference in two large adult population cohorts (Nordestgaard et al., 2015). Also, most MR studies do not support a role of caffeine consumption on BMI or waist circumstance (Nordestgaard et al., 2015; Cornelis and Munafo, 2018).
Chocolate consumptionMeta-analysis of RCTs did not support a significant effect of cocoa/dark chocolate supplementation on body weight or BMI (Kord-Varkaneh et al., 2019).NA
DiabetesNAMR study supported genetic predisposition to higher childhood BMI was associated with risk of type 2 diabetes (Geng et al., 2018).
Binge eatingNAThe MR study suggests a bi-directional association, that is, more binge eating and overeating are associated with higher BMI in later life, and higher children’s BMI is associated with more binge eating (Reed et al., 2017).
Physical exercisesA systematic review shows school based RCTs targeting physical activity, or physical activity combined with diet interventions, were effective in reducing BMI among children (Bleich et al., 2018).Genetically predicted more physical activity was associated with lower BMI (Carrasquilla et al., 2022).
SnoringNAMR study suggests genetically predicted BMI is related to snoring (Campos et al., 2020).
Pubertal timingNAGenetically predicted earlier age at puberty was related to higher BMI (Gill et al., 2018), but pleiotropy might exist (Gill et al., 2018). The association attenuated towards the null after controlling for prior BMI (Bell et al., 2018). It is also possible that genetically predicted BMI was associated with earlier age at puberty (Chen et al., 2019).
Birth weightNAMR analyses indicated positive causal associations of birthweight with BMI in the UKB (Zanetti et al., 2018).
Maternal adiposity during pregnancyNAUsing BMI polygenic risk scores calculated from maternal non-transmitted alleles, previous MR studies using mother–offspring pairs from two large UK cohorts did not support causal associations between maternal pre/early pregnancy BMI and offspring adolescent adiposity (Bond et al., 2022; Richmond et al., 2017).

In the epigenome-wide association study of 286 participants we identified 17 CpGs for BMI at ~23 years in the genes RBM16, SCN2B, SLC24A4, TECPR2, KSR1, RPTOR, GTF3C3, ZNF827, TXNDC15, C2, and RPS6KA2 and 17 for WHR in the genes LANCL2, C6orf195, MIR4535, CTRL, LYRM9, DCDC2, DIRC3, RPS6KA2, LPP, NFIC, MIR7641-2, ZNF141, RNF213, and OPA3 (Tables 7 and 8; Figures 3 and 4). cg14630200 in RPS6KA2 was a shared CpG for both BMI and WHR. The genetic inflation factors (lambda) were 1.006 and 1.005, respectively.

Epigenome-wide association with body mass index (BMI) at age ~23 years in 286 participants from Hong Kong’s ‘Children of 1997’ birth cohort in the follow-up survey.
Figure 3—source data 1

The associations for depicting Figure 3 in the study.

https://cdn.elifesciences.org/articles/82377/elife-82377-fig3-data1-v1.csv
Table 7
Associations of CpGs with body mass index (BMI) at age ~23 years in 286 participants from Hong Kong’s ‘Children of 1997’ birth cohort in the follow-up survey.
CpGGeneBetap value
cg00379266RBM16−69.79.8E−07
cg01525498RPTOR−40.24.6E−10
cg02200725-−18.94.2E−07
cg04563671SCN2B14.81.6E−07
cg07243141-24.69.0E−07
cg08205320-−21.92.6E−08
cg09639369SLC24A4−15.75.7E−07
cg11937362-1099.5E−07
cg12484266C2835.7E−09
cg14630200RPS6KA2−45.17.1E−12
cg14850190KSR1103.24.3E−07
cg18681028-−48.57.6E−16
cg18893311-30.87.5E−08
cg22197830TXNDC15−68.84.9E−07
cg25356423ZNF827−141.6E−07
cg25903177-−31.45.8E−07
cg26347007GTF3C3−45.88.0E−07
  1. CpGs reaching significance of p < 1 × 10−6 were shown in the table.

Table 8
Associations of CpGs with waist–hip ratio (WHR) at age ~23 years in 286 participants from Hong Kong’s ‘Children of 1997’ birth cohort in the follow-up survey.
CpGGeneBetap value
cg00952960LANCL2−0.543.0E−11
cg04405211C6orf1950.359.7E−07
cg05059349-0.263.7E−07
cg07416310DCDC2−0.458.8E−08
cg09422806NFIC0.53.4E−08
cg09684856DIRC3−0.791.4E−07
cg09907395RNF2131.092.2E−07
cg11344771LPP−0.791.0E−08
cg14630200RPS6KA2−0.465.4E−07
cg17007012MIR4535−0.356.6E−11
cg17043713ZNF1410.558.1E−08
cg17769811CTRL−0.761.3E−09
cg19481727MIR7641-2−0.365.4E−07
cg19579176-−0.359.9E−07
cg22364817-−1.466.7E−07
cg25097095LYRM9−0.418.8E−09
cg27384074OPA30.551.7E−07
  1. CpGs reaching significance of p < 1 × 10−6 were shown in the table.

Epigenome-wide association with waist–hip ratio (WHR) at age ~23 years in 286 participants from Hong Kong’s ‘Children of 1997’ birth cohort in the follow-up survey.

Source code file: source code for the environment- and epigenome-wide association study analyses.

Figure 4—source data 1

The associations for depicting Figure 4 in the study.

https://cdn.elifesciences.org/articles/82377/elife-82377-fig4-data1-v1.csv

Discussion

In this environment- and epigenome-wide association study, we systematically examined associations of over 400 exposures with obesity in a unique Chinese birth cohort, as well as the association of DNA methylation with obesity. Building on the previous studies in this birth cohort (Zhang et al., 2020; Hui et al., 2018; Lin et al., 2012; Kwok et al., 2010), we not only confirmed established risk factors, such as maternal second-hand smoking (Wang et al., 2014), but also added by identifying novel exposures not reported in previous EWAS in western settings (Vrijheid et al., 2020; Uche et al., 2020), such as consumption of ASB and soymilk. The comparison with RCTs or MR studies support a role of higher birth weight, dairy intake, binge eating, and possibly earlier puberty in obesity. We also identified several CpGs related to BMI and WHR in young Chinese, as reported in other populations (Kvaløy et al., 2018).

Our study found that maternal second-hand smoking was consistently associated with obesity at different ages, which is consistent with the concerns repeatedly raised in previous studies (Kwok et al., 2010; Wang et al., 2014), and adds support to the policy of banning smoking cigarettes and alternative smoking products in all indoor areas including workplaces and public places, as well as certain outdoor areas, such as open areas of schools, leisure facilities, bathing beaches, and public transport facilities in Hong Kong (WSC, 2022). Maternal weight is another maternal factor related to higher BMI and WHR consistently at different ages before adulthood. However, recent MR studies do not support a role of maternal overweight in offspring obesity (Bond et al., 2022; Richmond et al., 2017). Gestational diabetes was also identified to be associated with obesity at ~11.5 years, and the positive association remained for obesity at ~17.6 and ~23 years. It would be worthwhile to test its role in MR studies.

Regarding dietary factors, as the dietary assessments were more comprehensively conducted in the Biobank Clinical follow-up, the identified dietary factors were mainly for obesity at ~17.6 years. Interestingly, we found that children who consume more ASB have higher BMI, which is consistent with meta-analyses of cohort studies (Qin et al., 2020; Rousham et al., 2022). Consistent with a previous study in this birth cohort (Zhang et al., 2020), we did not find an association of sugar-sweetened beverages with obesity. The different associations for ASB and sugar-sweetened beverages might be because few consumed sugar-sweetened beverages regularly (6.8% consumed daily) (Zhang et al., 2020), while many consumed ASB (43% participants reported consumption). Consistent with our EWAS, an EWAS in the US also found that consumption of aspartame, a synthetic non-nutritive sweetener, was positively associated with abdominal obesity (Wulaningsih et al., 2017). ASB intake may induce appetite for similar sweet foods, leading to excess energy intake (Mattes and Popkin, 2009). The consistency across settings suggests this association is less likely to be confounded. However, whether it can be used as a target of intervention needs to be tested in a randomized, placebo-controlled trial (Rousham et al., 2022).

Another interesting finding is that milk consumption was not related to obesity at ~11.5 years, while reduced-fat/skim milk consumption was associated with higher BMI at ~17.6 years, with a consistent direction of associations for BMI at ~23 years. Our findings are consistent with a previous study in this cohort, which showed milk consumption frequency was not associated with BMI at 13 years (Lin et al., 2012), however, the previous study did not assess the specific type of milk. Our findings are different from a cross-sectional study in Portugal which shows more skimmed or semi-skimmed milk consumption was associated with lower abdominal obesity (Abreu et al., 2014). An explanation is that the observed associations of reduced-fat or skim milk with higher BMI could be due to increased muscle mass rather than body fat mass, or residual confounding by SEP. Alternatively, it might be due to reverse causality as young people with higher BMI might be more motivated to consume a specific diet. Interestingly, our findings are more consistent with an MR study suggesting genetically predicted higher dairy intake was associated with higher BMI (Huang et al., 2018b).

We also found tea or coffee consumption were associated with higher BMI at ~17.6 years. RCTs of coffee or tea consumption are scarce among children and adolescents because they require long-term adherence. MR studies do not suggest that coffee consumption affects obesity (Nordestgaard et al., 2015; Cornelis and Munafo, 2018), so the observed association might be due to confounding or chance. Similarly, the associations of chocolate and sweets intake with lower BMI at ~17.6 years, as well as physical activity with higher BMI at ~17.6 years are not consistent with RCTs of dark chocolate consumption (Kord-Varkaneh et al., 2019) and physical activity (Bleich et al., 2018), or MR studies of physical activity (Carrasquilla et al., 2022), and might be due to confounding or reverse causality.

Echoing the increasing attention to the role of mood and emotion in obesity control (Cardel et al., 2020), we found that binge eating was associated with higher BMI at ~17.6 years, with a consistent direction of association in the follow-up, consistent with an MR study (Reed et al., 2017). Our findings are also in line with the National Institute for Health and Care Excellence (NICE) guidance which also included binge eating in the consideration of children’s weight management (NICE, 2013). The underlying mechanism has not been clarified, but in general mental wellbeing may be linked to obesity via the neurohormonal weight control network concerning the hypothalamus (Sharma and Kavuru, 2010; Spiegel et al., 2009) as well as via psychosocial factors, lifestyle and behaviour.

As regards lifestyle, consistent with previous observational studies (Wang et al., 2019), we found sleep might play a role in childhood obesity at ~17.6 years. Despite a lack of RCTs, MR findings suggest sleep deprivation may be a causal factor for obesity (Wang et al., 2019; Dashti and Ordovás, 2021). Our study also shows coughing or snoring at night had a positive association with childhood BMI and WHR, which has not been identified in previous EWAS (Vrijheid et al., 2020; Uche et al., 2020). MR studies suggest genetically predicted BMI is positively associated with snoring (Campos et al., 2020), while the association of snoring with BMI is less clear. As such, we cannot exclude the possibility of reverse causality in our observation.

Consistent with previous studies (Lai et al., 2021), we found that earlier pubertal age for girls was related to higher BMI at ~17.6 years. Consistently, a previous study in this birth cohort suggested that maternal age at puberty was associated with offspring BMI at puberty (Lai et al., 2016). However, the findings from MR studies are controversial, with evidence showing genetically predicted earlier age at puberty related to higher BMI (Gill et al., 2018) while another MR study showing puberty timing has a small influence on BMI (Bell et al., 2018). Meanwhile, we cannot exclude a relation in the other direction (Chen et al., 2019); clarifying the bi-directional association would be worthwhile in future studies.

In the epigenome-wide association study, we found DNA methylation at RPS6KA2 was associated with both BMI and WHR, consistent with the previous epigenome-wide association studies of obesity in different populations (Kvaløy et al., 2018). Our study also identified several other genes, such as ZNF827, MIR7641-2, RAPTOR, KSR1, GTF3C3, and NFIC, whose role in obesity or obesity-related disorders has been consistently shown in previous studies. For example, ZNF827, MIR7641-2, and RAPTOR have been reported to be related to obesity (Huang et al., 2015; Dong et al., 2016) and/or overweight (Morris et al., 2015). KSR1 has been reported to be related to the regulation of glucose homeostasis (Klutho et al., 2011), and GTF3C3- to obesity-related dysglycaemia (Andrade et al., 2021). NFIC, which encodes nuclear factor I-C, regulates adipocyte differentiation (Zhou et al., 2017). Opa3, a novel regulator of mitochondrial function, controls thermogenesis and abdominal fat mass (Wells et al., 2012). The consistency of our study with other studies in different settings with different confounding structures suggests these association are less likely to be a product of confounding.

Strengths and limitations

To our knowledge, this study is the first study comprehensively assessing environmental factors related to obesity at the outset and at end of puberty in Asians, including some exposures specifically relevant in Asians, such as soymilk intake. We also replicated the associations in a follow-up survey and compared our findings with those from studies of different designs. Nevertheless, several limitations exist. First, the sample size for the EWAS and epigenome-wide association study is relatively small. Replication in a larger study is needed. The evidence from RCTs and MR is mainly in adults, which restricted the comparison with our findings. For better comparison, evidence from RCTs and MR conducted at similar ages are needed. Second, misclassification is possible for the exposures, which typically biases towards the null (Rothman, 2008). The use of questionnaires to ascertain exposures is prone to recall bias and social desirability bias, however, some previous studies within this cohort have suggested accurate reporting (Kwok et al., 2010). To minimize these possible biases, exposures measurements were collected using standard protocols and equipment with clear instructions. After accounting for the time difference between age at anthropometric measurements and age of exposure collection, we also obtained similar results. Third, we obtained similar results for BMI and obesity risk, however, given the lack of high-quality evidence about the cut-off values for waist circumference and waist-to-hip ratio in Asian children and adolescents, we did not perform logistic regression on central obesity risk. Fourth, the inconsistency between some of our findings and previous studies, such as chocolate, sweets, tea, and coffee consumption, should be interpreted cautiously. It may not only reflect differences between the West and China (Hong Kong), but also may be due to changes in structural socioeconomic and environmental factors, as well as changes in living environment, family relationships, social and community networks, housing and the care environment. Fifth, we only collected blood samples at the Biobank Clinical follow-up (age ~17.6 years), so we only conducted the epigenome-wide association study for DNA methylation at ~17.6 years. It would be worthwhile to examine the association of DNA methylation at different ages with obesity. Finally, although we controlled for several confounders in the multivariable analysis, residual confounding may still exist given our study is observational. Comparing our study with evidence from MR studies which are less likely to be confounded (Reed et al., 2017; Gill et al., 2018; Bell et al., 2018; Chen et al., 2019), we found a consistent direction of association for dairy intake (Huang et al., 2018b) and binge eating (Reed et al., 2017) being associated with higher BMI. Evidence for some exposures, such as tea and chocolate consumption, is still lacking; further MR studies are needed to assess causality.

Implications

In this study, we not only confirmed established risk factors, such as maternal second-hand smoking, but also identified several factors not reported or not examined in previous EWAS in western countries (Vrijheid et al., 2020; Uche et al., 2020), such as ASB consumption and soymilk intake. The comparison with RCTs or MR studies to some extent supports a role of dairy intake and binge eating, suggesting these factors or their drivers (e.g., sex hormones as a potential driver of binge eating) might be considered as potential targets for intervention. Other factors, such as soymilk intake, need to be tested in RCTs. We also identified several methylation loci related to obesity. Our study based on the unique setting in Hong Kong, provides potential drivers of obesity applicable to Hong Kong Chinese, with relevance to health policy interventions and future research.

Conclusions

This study takes advantage of the unique setting of Hong Kong and provides more insight about the role of environmental exposures and epigenetics in early life obesity. If these associations are found to be causal, they may provide novel intervention targets to improve population health.

Reporting

The study conforms with the STROBE checklist, which was attached as a supplementary file.

Data availability

The data used in this study are based on the 'Children of 1997' Birth cohort, maintained by the School of Public Health, The University of Hong Kong. With the approved ethics for this study, the individual participant data cannot be made freely available online. Interested parties can access the data used in this study upon reasonable request, with approval by the birth cohort team. As part of this process, researchers will be required to submit a project proposal for approval, to ensure the data is being used responsibly, ethically, and for scientifically sound projects. Requesters should be employees of a recognized academic institution, health service organization, or charitable research organization with experience in medical research. Requestors should be able to demonstrate, through their peer-reviewed publications in the area of interest, their ability to carry out the proposed study. Source data files have been uploaded for each of the results figures (Figures 1–4) showing the model summary data for plotting the Manhattan plots in environment-wide and epigenome-wide associations with BMI and WHR. Source code for the analyses has been uploaded as Source Code.

References

    1. Lee JPM
    (2016) Tobacco control policy in Hong Kong
    Hong Kong Medical Journal = Xianggang Yi Xue Za Zhi 22:96–97.
    https://doi.org/10.12809/hkmj164848
    1. Lozano R
    2. Naghavi M
    3. Foreman K
    4. Lim S
    5. Shibuya K
    6. Aboyans V
    7. Abraham J
    8. Adair T
    9. Aggarwal R
    10. Ahn SY
    11. Alvarado M
    12. Anderson HR
    13. Anderson LM
    14. Andrews KG
    15. Atkinson C
    16. Baddour LM
    17. Barker-Collo S
    18. Bartels DH
    19. Bell ML
    20. Benjamin EJ
    21. Bennett D
    22. Bhalla K
    23. Bikbov B
    24. Bin Abdulhak A
    25. Birbeck G
    26. Blyth F
    27. Bolliger I
    28. Boufous S
    29. Bucello C
    30. Burch M
    31. Burney P
    32. Carapetis J
    33. Chen H
    34. Chou D
    35. Chugh SS
    36. Coffeng LE
    37. Colan SD
    38. Colquhoun S
    39. Colson KE
    40. Condon J
    41. Connor MD
    42. Cooper LT
    43. Corriere M
    44. Cortinovis M
    45. de Vaccaro KC
    46. Couser W
    47. Cowie BC
    48. Criqui MH
    49. Cross M
    50. Dabhadkar KC
    51. Dahodwala N
    52. De Leo D
    53. Degenhardt L
    54. Delossantos A
    55. Denenberg J
    56. Des Jarlais DC
    57. Dharmaratne SD
    58. Dorsey ER
    59. Driscoll T
    60. Duber H
    61. Ebel B
    62. Erwin PJ
    63. Espindola P
    64. Ezzati M
    65. Feigin V
    66. Flaxman AD
    67. Forouzanfar MH
    68. Fowkes FGR
    69. Franklin R
    70. Fransen M
    71. Freeman MK
    72. Gabriel SE
    73. Gakidou E
    74. Gaspari F
    75. Gillum RF
    76. Gonzalez-Medina D
    77. Halasa YA
    78. Haring D
    79. Harrison JE
    80. Havmoeller R
    81. Hay RJ
    82. Hoen B
    83. Hotez PJ
    84. Hoy D
    85. Jacobsen KH
    86. James SL
    87. Jasrasaria R
    88. Jayaraman S
    89. Johns N
    90. Karthikeyan G
    91. Kassebaum N
    92. Keren A
    93. Khoo J-P
    94. Knowlton LM
    95. Kobusingye O
    96. Koranteng A
    97. Krishnamurthi R
    98. Lipnick M
    99. Lipshultz SE
    100. Ohno SL
    101. Mabweijano J
    102. MacIntyre MF
    103. Mallinger L
    104. March L
    105. Marks GB
    106. Marks R
    107. Matsumori A
    108. Matzopoulos R
    109. Mayosi BM
    110. McAnulty JH
    111. McDermott MM
    112. McGrath J
    113. Mensah GA
    114. Merriman TR
    115. Michaud C
    116. Miller M
    117. Miller TR
    118. Mock C
    119. Mocumbi AO
    120. Mokdad AA
    121. Moran A
    122. Mulholland K
    123. Nair MN
    124. Naldi L
    125. Narayan KMV
    126. Nasseri K
    127. Norman P
    128. O’Donnell M
    129. Omer SB
    130. Ortblad K
    131. Osborne R
    132. Ozgediz D
    133. Pahari B
    134. Pandian JD
    135. Rivero AP
    136. Padilla RP
    137. Perez-Ruiz F
    138. Perico N
    139. Phillips D
    140. Pierce K
    141. Pope CA
    142. Porrini E
    143. Pourmalek F
    144. Raju M
    145. Ranganathan D
    146. Rehm JT
    147. Rein DB
    148. Remuzzi G
    149. Rivara FP
    150. Roberts T
    151. De León FR
    152. Rosenfeld LC
    153. Rushton L
    154. Sacco RL
    155. Salomon JA
    156. Sampson U
    157. Sanman E
    158. Schwebel DC
    159. Segui-Gomez M
    160. Shepard DS
    161. Singh D
    162. Singleton J
    163. Sliwa K
    164. Smith E
    165. Steer A
    166. Taylor JA
    167. Thomas B
    168. Tleyjeh IM
    169. Towbin JA
    170. Truelsen T
    171. Undurraga EA
    172. Venketasubramanian N
    173. Vijayakumar L
    174. Vos T
    175. Wagner GR
    176. Wang M
    177. Wang W
    178. Watt K
    179. Weinstock MA
    180. Weintraub R
    181. Wilkinson JD
    182. Woolf AD
    183. Wulf S
    184. Yeh P-H
    185. Yip P
    186. Zabetian A
    187. Zheng Z-J
    188. Lopez AD
    189. Murray CJL
    190. AlMazroa MA
    191. Memish ZA
    (2012) Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the global burden of disease study 2010
    Lancet 380:2095–2128.
    https://doi.org/10.1016/S0140-6736(12)61728-0
  1. Book
    1. Rothman KJ
    (2008)
    Modern Epidemiology
    Wolters Kluwer.
  2. Book
    1. World Health Organization
    (2000)
    Regional Office for the Western P. The Asia-Pacific Perspective: Redefining Obesity and Its Treatment
    Health Communications Australia.

Article and author information

Author details

  1. Jie Zhao

    School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, Hong Kong
    Contribution
    Conceptualization, Resources, Data curation, Supervision, Investigation, Methodology, Writing – original draft, Writing – review and editing
    Contributed equally with
    Bohan Fan
    For correspondence
    janezhao@hku.hk
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1564-0057
  2. Bohan Fan

    School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, Hong Kong
    Contribution
    Formal analysis, Visualization, Writing – original draft
    Contributed equally with
    Jie Zhao
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5349-9461
  3. Jian Huang

    Singapore Institute for Clinical Sciences (SICS), Agency for Science, Technology and Research (A*STAR), Singapore, Singapore
    Contribution
    Formal analysis, Validation
    Competing interests
    No competing interests declared
  4. Benjamin John Cowling

    School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, Hong Kong
    Contribution
    Investigation, Writing – review and editing
    Competing interests
    Benjamin Cowling received research funding from Fosun Pharma, and received honoraria from AstraZeneca, Fosun Pharma, GlaxoSmithKline, Moderna, Pfizer, Roche and Sanofi Pasteur. The author has no other competing interests to declare
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6297-7154
  5. Shiu Lun Ryan Au Yeung

    School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, Hong Kong
    Contribution
    Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6136-1836
  6. Andrea Baccarelli

    Mailman School of Public Health, Columbia University, New York, United States
    Contribution
    Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
  7. Gabriel M Leung

    School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, Hong Kong
    Contribution
    Resources, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  8. C Mary Schooling

    1. School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong, Hong Kong
    2. City University of New York, School of Public Health and Health policy, New York, United States
    Contribution
    Resources, Investigation, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9933-5887

Funding

Health and Medical Research Fund (04180097)

  • Jie Zhao

Health and Medical Research Fund (CFS-HKU1)

  • Shiu Lun Ryan Au Yeung

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

Acknowledgements

This study was supported by the Health and Medical Research Fund Research Fellowship, Food and Health Bureau, Hong Kong SAR Government (#04180097). The DNA extraction was supported by CFS-HKU1. We would like to sincerely thank Food and Health Bureau for funding this project. We would also like thank all the participants and research staff in this project.

Ethics

This study complies with the Declaration of Helsinki. Since our participants are children, informed (non-written) consent for the original survey and subsequent record linkage was obtained from the parents, next of kin, caretakers, or guardians (informants) on behalf of the participants by the informant agreeing and subsequently completing the questionnaire at enrollment, this manner of obtaining consent was approved by The University of Hong Kong Medical Faculty Ethics Committee over 20 years ago. Informed written consent for subsequent surveys and in-person follow-up was obtained from a parent or guardian, or at ages 18+ years from the participant. Ethical approval for this study, including the follow-up survey at ~23 years and comprehensive health-related analyses, was obtained from the University of Hong Kong-Hospital Authority Hong Kong West Cluster, Joint Institutional Review Board, Hong Kong Special Administrative Region, China (reference numbers: UW13-367; UW19-367).

Copyright

© 2023, Zhao, Fan et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 707
    views
  • 91
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Jie Zhao
  2. Bohan Fan
  3. Jian Huang
  4. Benjamin John Cowling
  5. Shiu Lun Ryan Au Yeung
  6. Andrea Baccarelli
  7. Gabriel M Leung
  8. C Mary Schooling
(2023)
Environment- and epigenome-wide association study of obesity in ‘Children of 1997’ birth cohort
eLife 12:e82377.
https://doi.org/10.7554/eLife.82377

Share this article

https://doi.org/10.7554/eLife.82377

Further reading

    1. Epidemiology and Global Health
    Xiaoning Wang, Jinxiang Zhao ... Dong Liu
    Research Article

    Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.

    1. Epidemiology and Global Health
    2. Microbiology and Infectious Disease
    Amanda C Perofsky, John Huddleston ... Cécile Viboud
    Research Article

    Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.