Environment- and epigenome-wide association study of obesity in ‘Children of 1997’ birth cohort
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.sa0Introduction
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).
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).
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.
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.
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
-
Adipose tissue epigenetic profile in obesity-related dysglycemia-a systematic reviewFrontiers in Endocrinology 12:681649.https://doi.org/10.3389/fendo.2021.681649
-
Interventions to prevent global childhood overweight and obesity: a systematic reviewThe Lancet. Diabetes & Endocrinology 6:332–346.https://doi.org/10.1016/S2213-8587(17)30358-3
-
Multiple correlations and bonferroni’s correctionBiological Psychiatry 44:775–777.https://doi.org/10.1016/s0006-3223(98)00043-2
-
Genetics of sleep and insights into its relationship with obesityAnnual Review of Nutrition 41:223–252.https://doi.org/10.1146/annurev-nutr-082018-124258
-
Epigenomic elements analyses for promoters identify ESRRG as a new susceptibility gene for obesity-related traitsInternational Journal of Obesity 40:1170–1176.https://doi.org/10.1038/ijo.2016.44
-
Obesity and diabetes: the increased risk of cancer and cancer-related mortalityPhysiological Reviews 95:727–748.https://doi.org/10.1152/physrev.00030.2014
-
Age at menarche and adult body mass index: a Mendelian randomization studyInternational Journal of Obesity 42:1574–1581.https://doi.org/10.1038/s41366-018-0048-7
-
ConferenceEnvironment-wide association study (ewas) for type 2 diabetes in the marshfield personalized medicine research project biobankProceedings of the Pacific Symposium. pp. 200–211.https://doi.org/10.1142/9789814583220_0020
-
In utero exposure to gestational diabetes and adiposity: does breastfeeding make a difference?International Journal of Obesity 42:1317–1325.https://doi.org/10.1038/s41366-018-0077-2
-
Does cocoa/dark chocolate supplementation have favorable effect on body weight, body mass index and waist circumference? A systematic review, meta-analysis and dose-response of randomized clinical trialsCritical Reviews in Food Science and Nutrition 59:2349–2362.https://doi.org/10.1080/10408398.2018.1451820
-
Is a cutoff of 10 % appropriate for the change-in-estimate criterion of confounder identification?Journal of Epidemiology 24:161–167.https://doi.org/10.2188/jea.je20130062
-
Tobacco control policy in Hong KongHong Kong Medical Journal = Xianggang Yi Xue Za Zhi 22:96–97.https://doi.org/10.12809/hkmj164848
-
Cohort profile: family cohortInternational Journal of Epidemiology 46:e1.https://doi.org/10.1093/ije/dyu257
-
Informatics and data analytics to support exposome-based discovery for public healthAnnual Review of Public Health 38:279–294.https://doi.org/10.1146/annurev-publhealth-082516-012737
-
Nonnutritive sweetener consumption in humans: effects on appetite and food intake and their putative mechanismsThe American Journal of Clinical Nutrition 89:1–14.https://doi.org/10.3945/ajcn.2008.26792
-
Genetic variation in the raptor gene is associated with overweight but not hypertension in American men of Japanese ancestryAmerican Journal of Hypertension 28:508–517.https://doi.org/10.1093/ajh/hpu188
-
Coffee intake and risk of obesity, metabolic syndrome and type 2 diabetes: a Mendelian randomization studyInternational Journal of Epidemiology 44:551–565.https://doi.org/10.1093/ije/dyv083
-
Assessing the causal role of adiposity on disordered eating in childhood, adolescence, and adulthood: a Mendelian randomization analysisThe American Journal of Clinical Nutrition 106:764–772.https://doi.org/10.3945/ajcn.117.154104
-
Cohort profile: `` children of 1997'': a Hong Kong Chinese birth cohortInternational Journal of Epidemiology 41:611–620.https://doi.org/10.1093/ije/dyq243
-
Cohort profile: Hong Kong department of health elderly health service cohortInternational Journal of Epidemiology 45:64–72.https://doi.org/10.1093/ije/dyu227
-
Sleep and metabolism: an overviewInternational Journal of Endocrinology 2010:1–12.https://doi.org/10.1155/2010/270832
-
'Mendelian randomization': Can genetic epidemiology contribute to understanding environmental determinants of disease?International Journal of Epidemiology 32:1–22.https://doi.org/10.1093/ije/dyg070
-
Effects of poor and short sleep on glucose metabolism and obesity riskNature Reviews Endocrinology 5:253–261.https://doi.org/10.1038/nrendo.2009.23
-
Environment-wide association study on childhood obesity in the U.SEnvironmental Research 191:110109.https://doi.org/10.1016/j.envres.2020.110109
-
Early-Life environmental exposures and childhood obesity: an exposome-wide approachEnvironmental Health Perspectives 128:67009.https://doi.org/10.1289/EHP5975
-
Independent and joint effects of prenatal maternal smoking and maternal exposure to second-hand smoke on the development of adolescent obesity: a longitudinal studyJournal of Paediatrics and Child Health 50:908–915.https://doi.org/10.1111/jpc.12667
-
BookRegional Office for the Western P. The Asia-Pacific Perspective: Redefining Obesity and Its TreatmentHealth Communications Australia.
-
Investigating nutrition and lifestyle factors as determinants of abdominal obesity: an environment-wide studyInternational Journal of Obesity 41:340–347.https://doi.org/10.1038/ijo.2016.203
-
Genetically predicted milk consumption and bone health, ischemic heart disease and type 2 diabetes: A Mendelian randomization studyEuropean Journal of Clinical Nutrition 71:1008–1012.https://doi.org/10.1038/ejcn.2017.8
Article and author information
Author details
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
-
- 694
- views
-
- 89
- downloads
-
- 2
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Epidemiology and Global Health
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.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
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.