Exploring the role of genetic variability and exposure to bisphenols and parabens on excess body weight in Spanish children

Gene-environment interaction studies are emerging as a promising tool to shed light on the reasons for the rapid increase in excess body weight (overweight and obesity). We aimed to investigate the influence of several polymorphisms on excess weight in Spanish children according to a short- and long-term exposure to bisphenols and parabens, combining individual approach with the joint effect of them. This case-control study included 144 controls and 98 cases children aged 3–12 years. Thirty SNPs in genes involved in obesity-related pathways, xenobiotic metabolism and hormone systems were genotyped using the GSA microchip technology and qPCRs with Taqman® probes. Levels of bisphenols and parabens in urine and hair were used to assess short- and long-term exposure, respectively, via UHPLC-MS/MS system. LEPR rs9436303 was identified as a relevant risk variant for excess weight (ORDom:AAvsAG+GG=2.65, p<0.001), and this effect persisted across exposure-stratified models. For long-term exposure, GPX1 rs1050450 was associated with increased excess weight at low single paraben exposure (ORGvsA=2.00, p=0.028, p-interaction=0.016), whereas LEPR rs1137101 exhibited a protective function at high co-exposure (ORDom:AAvsAG+GG=0.17, p=0.007, p-interaction=0.043). ESR2 rs3020450 (ORDom:GGvsAG+AA=5.17, p=0.020, p-interaction=0.028) and CYP2C19 rs4244285 (ORDom:GGvsAG+AA=3.54, p=0.039, p-interaction=0.285) were identified as predisposing variants at low and high co-exposure, respectively. In short-term exposure, higher odds were observed for INSIG2 rs7566605 at high bisphenol exposure (ORCvsG=2.97, p=0.035, p-interaction=0.017) and for GSTP1 rs1695 at low levels (ORDom:AAvsAG+GG=5.38, p=0.016, p-interaction=0.016). At low and medium co-exposure, SH2B1 rs7498665 (ORAvsG=0.17, p=0.015, p-interaction=0.085) and MC4R rs17782313 (ORAvsG=0.10, p=0.023, p-interaction=0.045) displayed a protective effect, whereas ESR2 rs3020450 maintained its contributing role (ORGvsA=3.12, p=0.030, p-interaction=0.010). Our findings demonstrate for the first time that understanding the genetic variation in excess weight and how the level of exposure to bisphenols and parabens might interact with it, is crucial for a more in-depth comprehension of the complex polygenic and multifactorial aetiology of overweight and obesity.