Rare variant associations with birth weight identify genes involved in adipose tissue regulation, placental function and insulin-like growth factor signalling.

Investigating the genetic factors influencing human birth weight may lead to biological insights into fetal growth and long-term health. Genome-wide association studies of birth weight have highlighted associated variants in more than 200 regions of the genome, but the causal genes are mostly unknown. Rare genetic variants with robust evidence of association are more likely to point to causal genes, but to date, only a few rare variants are known to influence birth weight. We aimed to identify genes that harbour rare variants that impact birth weight when carried by either the fetus or the mother, by analysing whole exome sequence data in UK Biobank participants. We annotated rare (minor allele frequency <0.1%) protein-truncating or high impact missense variants on whole exome sequence data in up to 234,675 participants with data on their own birth weight (fetal variants), and up to 181,883 mothers who reported the birth weight of their first child (maternal variants). Variants within each gene were collapsed to perform gene burden tests and for each associated gene, we compared the observed fetal and maternal effects. We identified 8 genes with evidence of rare fetal variant effects on birth weight, of which 2 also showed maternal effects. One additional gene showed evidence of maternal effects only. We observed 10/11 directionally concordant associations in an independent sample of up to 45,622 individuals (sign test P=0.01). Of the genes identified, IGF1R and PAPPA2 (fetal and maternal-acting) have known roles in insulin-like growth factor bioavailability and signalling. PPARG, INHBE and ACVR1C (all fetal-acting) have known roles in adipose tissue regulation and rare variants in the latter two also showed associations with favourable adiposity patterns in adults. We highlight the dual role of PPARG in both adipocyte differentiation and placental angiogenesis. NOS3, NRK, and ADAMTS8 (fetal and maternal-acting) have been implicated in both placental function and hypertension. Analysis of rare coding variants has identified regulators of fetal adipose tissue and fetoplacental angiogenesis as determinants of birth weight, as well as further evidence for the role of insulin-like growth factors.

Protein-truncating variants in BSN are associated with severe adult-onset obesity, type 2 diabetes and fatty liver disease.

Obesity is a major risk factor for many common diseases and has a substantial heritable component. To identify new genetic determinants, we performed exome-sequence analyses for adult body mass index (BMI) in up to 587,027 individuals. We identified rare loss-of-function variants in two genes (BSN and APBA1) with effects substantially larger than those of well-established obesity genes such as MC4R. In contrast to most other obesity-related genes, rare variants in BSN and APBA1 were not associated with normal variation in childhood adiposity. Furthermore, BSN protein-truncating variants (PTVs) magnified the influence of common genetic variants associated with BMI, with a common variant polygenic score exhibiting an effect twice as large in BSN PTV carriers than in noncarriers. Finally, we explored the plasma proteomic signatures of BSN PTV carriers as well as the functional consequences of BSN deletion in human induced pluripotent stem cell-derived hypothalamic neurons. Collectively, our findings implicate degenerative processes in synaptic function in the etiology of adult-onset obesity.

Mendelian randomization identifies circulating proteins as biomarkers for age at menarche and age at natural menopause.

Age at menarche (AAM) and age at natural menopause (ANM) are highly heritable traits and have been linked to various health outcomes. We aimed to identify circulating proteins associated with altered ANM and AAM using an unbiased two-sample Mendelian randomization (MR) and colocalization approach. By testing causal effects of 1,271 proteins on AAM, we identified 22 proteins causally associated with AAM in MR, among which 13 proteins (GCKR, FOXO3, SEMA3G, PATE4, AZGP1, NEGR1, LHB, DLK1, ANXA2, YWHAB, DNAJB12, RMDN1 and HPGDS) colocalized. Among 1,349 proteins tested for causal association with ANM using MR, we identified 19 causal proteins among which 7 proteins (CPNE1, TYMP, DNER, ADAMTS13, LCT, ARL and PLXNA1) colocalized. Follow-up pathway and gene enrichment analyses demonstrated links between AAM-related proteins and obesity and diabetes, and between AAM and ANM-related proteins and various types of cancer. In conclusion, we identified proteomic signatures of reproductive ageing in women, highlighting biological processes at both ends of the reproductive lifespan.

Metabolic control of puberty: 60 years in the footsteps of Kennedy and Mitra's seminal work.

An individual's nutritional status has a powerful effect on sexual maturation. Puberty onset is delayed in response to chronic energy insufficiency and is advanced under energy abundance. The consequences of altered pubertal timing for human health are profound. Late puberty increases the chances of cardiometabolic, musculoskeletal and neurocognitive disorders, whereas early puberty is associated with increased risks of adult obesity, type 2 diabetes mellitus, cardiovascular diseases and various cancers, such as breast, endometrial and prostate cancer. Kennedy and Mitra's trailblazing studies, published in 1963 and using experimental models, were the first to demonstrate that nutrition is a key factor in puberty onset. Building on this work, the field has advanced substantially in the past decade, which is largely due to the impressive development of molecular tools for experimentation and population genetics. In this Review, we discuss the latest advances in basic and translational sciences underlying the nutritional and metabolic control of pubertal development, with a focus on perspectives and future directions.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Associations of appetitive traits with growth velocities from infancy to childhood.

Several studies have reported associations between appetitive traits and weight gain during infancy or childhood, but none have directly compared these associations across both age periods. Here, we tested the associations between appetitive traits and growth velocities from birth to childhood. Appetitive trait data were collected using the Children's Eating Behaviour Questionnaire (CEBQ) in 149 children from the Cambridge Baby Growth Study at age 9-17 years. These participants also provided anthropometric measurements during infancy (birth, 3, 12, 18, and 24 months) and childhood (5 to 11 years). Standardized growth velocities (in weight, length/height, BMI, and body fat percentage) for 0-3 months, 3-24 months, and 24 months to childhood were estimated using individual linear-spline models. Associations between each of the eight CEBQ traits and each growth velocity were tested in separate multilevel linear regression models, adjusted for sex, age at CEBQ completion, and the corresponding birth measurement (weight, length, BMI, or body fat percentage). The three food-approach traits (food responsiveness, enjoyment of food and emotional overeating) were positively associated with infancy and childhood growth velocities in weight, BMI, and body fat percentage. By contrast, only one of the food-avoidant traits, satiety responsiveness, was negatively associated with all growth velocities. Significant associations were mostly of similar magnitude across all age periods. These findings reveal a broadly consistent relationship between appetitive traits with gains in weight and adiposity throughout infancy and childhood. Future interventions and strategies to prevent obesity may benefit from measuring appetitive traits in infants and children and targeting these as part of their programs.

Large-scale exome sequence analysis identifies sex- and age-specific determinants of obesity.

Obesity contributes substantially to the global burden of disease and has a significant heritable component. Recent large-scale exome sequencing studies identified several genes in which rare, protein-coding variants have large effects on adult body mass index (BMI). Here we extended such work by performing sex-stratified associations in the UK Biobank study (N∼420,000). We identified genes in which rare heterozygous loss-of-function increases adult BMI in women (DIDO1, PTPRG, and SLC12A5) and in men (SLTM), with effect sizes up to ∼8 kg/m2. This is complemented by analyses implicating rare variants in OBSCN and MADD for recalled childhood adiposity. The known functions of these genes, as well as findings of common variant genome-wide pathway enrichment analyses, suggest a role for neuron death, apoptosis, and DNA damage response mechanisms in the susceptibility to obesity across the life-course. These findings highlight the importance of considering sex-specific and life-course effects in the genetic regulation of obesity.

Genetic impacts on DNA methylation help elucidate regulatory genomic processes.

BACKGROUND: Pinpointing genetic impacts on DNA methylation can improve our understanding of pathways that underlie gene regulation and disease risk. RESULTS: We report heritability and methylation quantitative trait locus (meQTL) analysis at 724,499 CpGs profiled with the Illumina Infinium MethylationEPIC array in 2358 blood samples from three UK cohorts. Methylation levels at 34.2% of CpGs are affected by SNPs, and 98% of effects are cis-acting or within 1 Mbp of the tested CpG. Our results are consistent with meQTL analyses based on the former Illumina Infinium HumanMethylation450 array. Both SNPs and CpGs with meQTLs are overrepresented in enhancers, which have improved coverage on this platform compared to previous approaches. Co-localisation analyses across genetic effects on DNA methylation and 56 human traits identify 1520 co-localisations across 1325 unique CpGs and 34 phenotypes, including in disease-relevant genes, such as USP1 and DOCK7 (total cholesterol levels), and ICOSLG (inflammatory bowel disease). Enrichment analysis of meQTLs and integration with expression QTLs give insights into mechanisms underlying cis-meQTLs (e.g. through disruption of transcription factor binding sites for CTCF and SMC3) and trans-meQTLs (e.g. through regulating the expression of ACD and SENP7 which can modulate DNA methylation at distal sites). CONCLUSIONS: Our findings improve the characterisation of the mechanisms underlying DNA methylation variability and are informative for prioritisation of GWAS variants for functional follow-ups. The MeQTL EPIC Database and viewer are available online at https://epicmeqtl.kcl.ac.uk .

Rare variants in the MECP2 gene in girls with central precocious puberty: a translational cohort study.

BACKGROUND: Identification of genetic causes of central precocious puberty have revealed epigenetic mechanisms as regulators of human pubertal timing. MECP2, an X-linked gene, encodes a chromatin-associated protein with a role in gene transcription. MECP2 loss-of-function mutations usually cause Rett syndrome, a severe neurodevelopmental disorder. Early pubertal development has been shown in several patients with Rett syndrome. The aim of this study was to explore whether MECP2 variants are associated with an idiopathic central precocious puberty phenotype. METHODS: In this translational cohort study, participants were recruited from seven tertiary centres from five countries (Brazil, Spain, France, the USA, and the UK). Patients with idiopathic central precocious puberty were investigated for rare potentially damaging variants in the MECP2 gene, to assess whether MECP2 might contribute to the cause of central precocious puberty. Inclusion criteria were the development of progressive pubertal signs (Tanner stage 2) before the age of 8 years in girls and 9 years in boys and basal or GnRH-stimulated LH pubertal concentrations. Exclusion criteria were the diagnosis of peripheral precocious puberty and the presence of any recognised cause of central precocious puberty (CNS lesions, known monogenic causes, genetic syndromes, or early exposure to sex steroids). All patients included were followed up at the outpatient clinics of participating academic centres. We used high-throughput sequencing in 133 patients and Sanger sequencing of MECP2 in an additional 271 patients. Hypothalamic expression of Mecp2 and colocalisation with GnRH neurons were determined in mice to show expression of Mecp2 in key nuclei related to pubertal timing regulation. FINDINGS: Between Jun 15, 2020, and Jun 15, 2022, 404 patients with idiopathic central precocious puberty (383 [95%] girls and 21 [5%] boys; 261 [65%] sporadic cases and 143 [35%] familial cases from 134 unrelated families) were enrolled and assessed. We identified three rare heterozygous likely damaging coding variants in MECP2 in five girls: a de novo missense variant (Arg97Cys) in two monozygotic twin sisters with central precocious puberty and microcephaly; a de novo missense variant (Ser176Arg) in one girl with sporadic central precocious puberty, obesity, and autism; and an insertion (Ala6_Ala8dup) in two unrelated girls with sporadic central precocious puberty. Additionally, we identified one rare heterozygous 3'UTR MECP2 insertion (36_37insT) in two unrelated girls with sporadic central precocious puberty. None of them manifested Rett syndrome. Mecp2 protein colocalised with GnRH expression in hypothalamic nuclei responsible for GnRH regulation in mice. INTERPRETATION: We identified rare MECP2 variants in girls with central precocious puberty, with or without mild neurodevelopmental abnormalities. MECP2 might have a role in the hypothalamic control of human pubertal timing, adding to the evidence of involvement of epigenetic and genetic mechanisms in this crucial biological process. FUNDING: Fundação de Amparo à Pesquisa do Estado de São Paulo, Conselho Nacional de Desenvolvimento Científico e Tecnológico, and the Wellcome Trust.

The association between the timing, intensity and magnitude of adolescent growth and body composition in early adulthood.

OBJECTIVES: There's paucity of longitudinal studies assessing the role of adolescent growth on adult body composition in developing countries. The aims of this study were to assess the association between adolescent change in height, weight and BMI and early adult height, weight, body fat and lean mass. METHODS: Magnitude, timing and intensity of height, weight and BMI growth were modelled for participants from the Birth to Thirty (Bt30) cohort (7-23 years). Early adult height, weight, BMI and DXA-derived body composition were obtained 1881 black participants (21-24 years). Linear regression analyses were used to assess associations. RESULTS: Adolescents with an earlier onset of puberty were heavier in childhood and had an earlier timing and faster weight gain velocity in late adolescence. The intensity of adolescent weight gain was positively associated with adult BMI and fat mass index (FMI) in females. Early timing of adolescent BMI gain was associated with increased weight and BMI in adult females and FMI in adult males. Achieving peak weight velocity around age at peak height velocity was associated with lower BMI and fat mass in both sexes. CONCLUSION: This study confirms the adverse consequences of excessive weight gain prior to puberty, which is associated with an earlier and faster resurgence in weight gain velocity in early adulthood. Factors that contribute to an asynchronous timing of ages of peak weight and peak height velocities may accentuate the risk of adult obesity.

Penetrance of pathogenic genetic variants associated with premature ovarian insufficiency.

Premature ovarian insufficiency (POI) affects 1% of women and is a leading cause of infertility. It is often considered to be a monogenic disorder, with pathogenic variants in ~100 genes described in the literature. We sought to systematically evaluate the penetrance of variants in these genes using exome sequence data in 104,733 women from the UK Biobank, 2,231 (1.14%) of whom reported at natural menopause under the age of 40 years. We found limited evidence to support any previously reported autosomal dominant effect. For nearly all heterozygous effects on previously reported POI genes, we ruled out even modest penetrance, with 99.9% (13,699 out of 13,708) of all protein-truncating variants found in reproductively healthy women. We found evidence of haploinsufficiency effects in several genes, including TWNK (1.54 years earlier menopause, P = 1.59 × 10-6) and SOHLH2 (3.48 years earlier menopause, P = 1.03 × 10-4). Collectively, our results suggest that, for the vast majority of women, POI is not caused by autosomal dominant variants either in genes previously reported or currently evaluated in clinical diagnostic panels. Our findings, plus previous studies, suggest that most POI cases are likely oligogenic or polygenic in nature, which has important implications for future clinical genetic studies, and genetic counseling for families affected by POI.

Prevalence of Deleterious Variants in MC3R in Patients With Constitutional Delay of Growth and Puberty.

CONTEXT: The melanocortin 3 receptor (MC3R) has recently emerged as a critical regulator of pubertal timing, linear growth, and the acquisition of lean mass in humans and mice. In population-based studies, heterozygous carriers of deleterious variants in MC3R report a later onset of puberty than noncarriers. However, the frequency of such variants in patients who present with clinical disorders of pubertal development is currently unknown. OBJECTIVE: This work aimed to determine whether deleterious MC3R variants are more frequently found in patients clinically presenting with constitutional delay of growth and puberty (CDGP) or normosmic idiopathic hypogonadotropic hypogonadism (nIHH). METHODS: We examined the sequence of MC3R in 362 adolescents with a clinical diagnosis of CDGP and 657 patients with nIHH, experimentally characterized the signaling properties of all nonsynonymous variants found and compared their frequency to that in 5774 controls from a population-based cohort. Additionally, we established the relative frequency of predicted deleterious variants in individuals with self-reported delayed vs normally timed menarche/voice-breaking in the UK Biobank cohort. RESULTS: MC3R loss-of-function variants were infrequent but overrepresented in patients with CDGP (8/362 [2.2%]; OR = 4.17; P = .001). There was no strong evidence of overrepresentation in patients with nIHH (4/657 [0.6%]; OR = 1.15; P = .779). In 246 328 women from the UK Biobank, predicted deleterious variants were more frequently found in those self-reporting delayed (aged ≥16 years) vs normal age at menarche (OR = 1.66; P = 3.90E-07). CONCLUSION: We have found evidence that functionally damaging variants in MC3R are overrepresented in individuals with CDGP but are not a common cause of this phenotype.

Leptin and IGF-1 in Infancy Are Associated With Variants in DHCR7 and CYP2R1 That Relate With Type 1 Diabetes and 25OHD.

CONTEXT: Vitamin D has been variably implicated in risk of developing type 1 diabetes based on cohorts of at-risk individuals. Emergent type 1 diabetes in childhood is putatively preceded by altered growth. OBJECTIVE: We explored whether polymorphisms in vitamin D metabolism genes modify risk of type 1 diabetes via effects on growth in a prospective, population-based cohort of infants. METHODS: The Cambridge Baby Growth Study enrolled newborns from Cambridgeshire, UK, for follow-up in infancy. In 612 infants, we genotyped single nucleotide polymorphisms in vitamin D metabolism genes that relate with type 1 diabetes: rs10741657 and rs12794714 in CYP2R1, rs12785878 in DHCR7, and rs10877012 in CYP27B1. Multivariate linear regression analyses tested associations between genotypes and anthropometric indices (weight, length, and skinfold thickness) or growth-related hormones (C-peptide, IGF-1, and leptin) in infancy. RESULTS: Birth weight showed borderline associations with the diabetes risk-increasing alleles in CYP2R1, rs10741657 (ß = -.11, P = .02) and rs12794714 (ß = -.09, P = .04). The risk-increasing allele rs12794714 was also associated with higher IGF-1 levels at age 24 months (ß = .30, P = .01). At age 3 months, the risk-increasing allele rs12785878 in DHCR7, known to negatively associate with 25-hydroxyvitamin D levels, showed a positive association with leptin levels (ß = .23, P = .009), which was pronounced in girls (P = .004) vs boys (P = .7). CONCLUSION: The vitamin D metabolism genes DHCR7 and CYP2R1 might influence infancy leptin and IGF-1 levels respectively. These findings open the possibility for a developmental role of vitamin D that is mediated by growth-related hormones with implications for the onset of type 1 diabetes autoimmunity.

Pre-pubertal accelerometer-assessed physical activity and timing of puberty in British boys and girls: the Millennium Cohort Study.

BACKGROUND: Early puberty timing is associated with adverse health outcomes. We aimed to examine prospective associations between objectively measured physical activity and puberty timing in boys and girls. METHODS: In the UK Millennium Cohort Study, physical activity volume and intensities at 7 years were measured using accelerometers. Status of several pubertal traits and age at menarche were reported at 11, 14 and 17 years. Age at menarche in girls was categorized into tertiles. Other puberty traits were categorized into earlier or later than the median ages calculated from probit models, separately in boys and girls. Multivariable regression models, with adjustment for maternal and child characteristics including body mass index (BMI) at age 7 years as potential confounders, were performed to test the associations of total daily activity counts and fractions of activity counts across intensities (in compositional models) with puberty timing, separately in boys (n = 2531) and girls (n = 3079). RESULTS: Higher total daily activity counts were associated with lower risks for earlier (vs later) growth spurt, body hair growth, skin changes and menarche in girls, and more weakly with lower risks for earlier skin changes and voice breaking in boys (odds ratios = 0.80-0.87 per 100 000 counts/day). These associations persisted on additional adjustment for BMI at 11 years as a potential mediator. No association with puberty timing was seen for any physical activity intensity (light, moderate or vigorous). CONCLUSIONS: More physical activity regardless of intensity may contribute to the avoidance of earlier puberty timing, independently of BMI, particularly in girls.

Early development of infant gut microbiota in relation to breastfeeding and human milk oligosaccharides.

Background: Infant gut microbiota composition is influenced by various factors early in life. Here, we investigate associations between infant gut microbiome development, infant age, breastfeeding duration, and human milk oligosaccharides (HMO) composition in breastmilk. Methods: A total of 94 mother-infant pairs were recruited as part of the Cambridge Baby Growth and Breastfeeding Study (CBGS-BF) (Cambridge, UK). Infant stool samples (n = 337) were collected at 2 week, 6 week, 3 month, and 6 month of age. The 16S rRNA V3-V4 rRNA region was sequenced using MiSeq Illumina to determine microbiota composition and diversity. Mother's hindmilk samples were collected at birth, 2 week, 6 week, 3 month, and 6 month postpartum. Concentrations of five neutral [2'FL, 3'FL, lacto-N-fucopentaose 1 (LNFP1), LNnT, LNT] and two acidic (3'SL, and 6'SL) HMOs were measured in all milk samples using High-Performance Anion-Exchange Chromatography with Pulsed Amperometric Detection (HPAEC-PAD). We explored the associations between infant gut microbiome parameters and age, duration of exclusive breastfeeding (EBF), and levels of individual HMOs. Results: Bifidobacterium was the most abundant genus in infant stool at all-time points, irrespective of breastfeeding duration, with an overall mean relative abundance of 70%. The relative abundance of B. bifidum in stool from infants who were breastfed for longer than 6 months was significantly higher compared to the infant breastfed up to 3 months (p = 0.0285). Alpha-diversity (both Shannon and ASV-level Richness) of infant gut microbiota showed a biphasic change with infant age, decreasing from 2 weeks until 3 months and then increasing until 6 months of age. Bifidobacterium relative abundance was associated with higher concentrations of 2'FL and LNFP1 in breastmilk across all time-points (p = 0.049 and 0.017, respectively), with trends toward a higher abundance of B. longum species. No significant association with Bifidobacterium was found for breastmilk LNnT, 3'SL, and 6'SL levels. Conclusion: Our study is in line with previous data demonstrating that EBF duration in the first months of life impacts infant gut microbiota composition. The observed links between specific HMOs in breastmilk and bacteria in infant stool provide evidence of how mother's milk affects infant microbiome development.

Butyrate in Human Milk: Associations with Milk Microbiota, Milk Intake Volume, and Infant Growth.

Butyrate in human milk (HM) has been suggested to reduce excessive weight and adipo-sity gains during infancy. However, HM butyrate's origins, determinants, and its influencing mechanism on weight gain are not completely understood. These were studied in the prospective longitudinal Cambridge Baby Growth and Breastfeeding Study (CBGS-BF), in which infants (n = 59) were exclusively breastfed for at least 6 weeks. Infant growth (birth, 2 weeks, 6 weeks, 3 months, 6 months, and 12 months) and HM butyrate concentrations (2 weeks, 6 weeks, 3 months, and 6 months) were measured. At age 6 weeks, HM intake volume was measured by deuterium-labelled water technique and HM microbiota by 16S sequencing. Cross-sectionally at 6 weeks, HM butyrate was associated with HM microbiota composition (p = 0.036) although no association with the abundance of typical butyrate producers was detected. In longitudinal analyses across all time points, HM butyrate concentrations were overall negatively associated with infant weight and adiposity, and associations were stronger at younger infant ages. HM butyrate concentration was also inversely correlated with HM intake volume, supporting a possible mechanism whereby butyrate might reduce infant growth via appetite regulation and modulation of HM intake.

Associations between breast milk intake volume, macronutrient intake and infant growth in a longitudinal birth cohort: the Cambridge Baby Growth and Breastfeeding Study (CBGS-BF).

Growth patterns of breastfed infants show substantial inter-individual differences, partly influenced by breast milk (BM) nutritional composition. However, BM nutritional composition does not accurately indicate BM nutrient intakes. This study aimed to examine the associations between both BM intake volumes and macronutrient intakes with infant growth. Mother-infant dyads (n 94) were recruited into the Cambridge Baby Growth and Breastfeeding Study (CBGS-BF) from a single maternity hospital at birth; all infants received exclusive breast-feeding (EBF) for at least 6 weeks. Infant weight, length and skinfolds thicknesses (adiposity) were repeatedly measured from birth to 12 months. Post-feed BM samples were collected at 6 weeks to measure TAG (fat), lactose (carbohydrate) (both by 1H-NMR) and protein concentrations (Dumas method). BM intake volume was estimated from seventy infants between 4 and 6 weeks using dose-to-the-mother deuterium oxide (2H2O) turnover. In the full cohort and among sixty infants who received EBF for 3+ months, higher BM intake at 6 weeks was associated with initial faster growth between 0 and 6 weeks (ß + se 3·58 + 0·47 for weight and 4·53 + 0·6 for adiposity gains, both P < 0·0001) but subsequent slower growth between 3 and 12 months (ß + se - 2·27 + 0·7 for weight and -2·65 + 0·69 for adiposity gains, both P < 0·005). BM carbohydrate and protein intakes at 4-6 weeks were positively associated with early (0-6 weeks) but tended to be negatively related with later (3-12 months) adiposity gains, while BM fat intake showed no association, suggesting that carbohydrate and protein intakes may have more functional relevance to later infant growth and adiposity.

Damaging missense variants in IGF1R implicate a role for IGF-1 resistance in the etiology of type 2 diabetes.

Type 2 diabetes (T2D) is a heritable metabolic disorder. While population studies have identified hundreds of common genetic variants associated with T2D, the role of rare (frequency < 0.1%) protein-coding variation is less clear. We performed exome sequence analysis in 418,436 (n = 32,374 T2D cases) individuals in the UK Biobank. We identified previously reported genes (GCK, GIGYF1, HNF1A) in addition to missense variants in ZEB2 (n = 31 carriers; odds ratio [OR] = 5.5 [95% confidence interval = 2.5-12.0]; p = 6.4 × 10-7), MLXIPL (n = 245; OR = 2.3 [1.6-3.2]; p = 3.2 × 10-7), and IGF1R (n = 394; OR = 2.4 [1.8-3.2]; p = 1.3 × 10-10). Carriers of damaging missense variants within IGF1R were also shorter (-2.2 cm [-1.8 to -2.7]; p = 1.2 × 10-19) and had higher circulating insulin-like growth factor-1 (IGF-1) protein levels (2.3 nmol/L [1.7-2.9]; p = 2.8 × 10-14), indicating relative IGF-1 resistance. A likely causal role of IGF-1 resistance was supported by Mendelian randomization analyses using common variants. These results increase understanding of the genetic architecture of T2D and highlight the growth hormone/IGF-1 axis as a potential therapeutic target.

Lipid profiling analyses from mouse models and human infants.

This protocol outlines a translational lipidomic approach to discover lipid biomarkers that could predict morphometric body and histological organ measurements (e.g., weight and adiposity gains) during specific stages of life (e.g., early life). We describe procedures ranging from animal experimentation and histological analyses to downstream analytical steps through lipid profiling, both in mice and humans. This protocol represents a reliable and versatile approach to translate and validate candidate lipid biomarkers from animal models to a human cohort. For complete details on the use and execution of this protocol, please refer to Olga et al. (2021).