Polygenic Prediction of Weight and Obesity Trajectories from Birth to Adulthood.

Severe obesity is a rapidly growing global health threat. Although often attributed to unhealthy lifestyle choices or environmental factors, obesity is known to be heritable and highly polygenic; the majority of inherited susceptibility is related to the cumulative effect of many common DNA variants. Here we derive and validate a new polygenic predictor comprised of 2.1 million common variants to quantify this susceptibility and test this predictor in more than 300,000 individuals ranging from middle age to birth. Among middle-aged adults, we observe a 13-kg gradient in weight and a 25-fold gradient in risk of severe obesity across polygenic score deciles. In a longitudinal birth cohort, we note minimal differences in birthweight across score deciles, but a significant gradient emerged in early childhood and reached 12 kg by 18 years of age. This new approach to quantify inherited susceptibility to obesity affords new opportunities for clinical prevention and mechanistic assessment.

Human Gain-of-Function MC4R Variants Show Signaling Bias and Protect against Obesity.

The melanocortin 4 receptor (MC4R) is a G protein-coupled receptor whose disruption causes obesity. We functionally characterized 61 MC4R variants identified in 0.5 million people from UK Biobank and examined their associations with body mass index (BMI) and obesity-related cardiometabolic diseases. We found that the maximal efficacy of ß-arrestin recruitment to MC4R, rather than canonical Gαs-mediated cyclic adenosine-monophosphate production, explained 88% of the variance in the association of MC4R variants with BMI. While most MC4R variants caused loss of function, a subset caused gain of function; these variants were associated with significantly lower BMI and lower odds of obesity, type 2 diabetes, and coronary artery disease. Protective associations were driven by MC4R variants exhibiting signaling bias toward ß-arrestin recruitment and increased mitogen-activated protein kinase pathway activation. Harnessing ß-arrestin-biased MC4R signaling may represent an effective strategy for weight loss and the treatment of obesity-related cardiometabolic diseases.

Trim28 Haploinsufficiency Triggers Bi-stable Epigenetic Obesity.

More than one-half billion people are obese, and despite progress in genetic research, much of the heritability of obesity remains enigmatic. Here, we identify a Trim28-dependent network capable of triggering obesity in a non-Mendelian, "on/off" manner. Trim28(+/D9) mutant mice exhibit a bi-modal body-weight distribution, with isogenic animals randomly emerging as either normal or obese and few intermediates. We find that the obese-"on" state is characterized by reduced expression of an imprinted gene network including Nnat, Peg3, Cdkn1c, and Plagl1 and that independent targeting of these alleles recapitulates the stochastic bi-stable disease phenotype. Adipose tissue transcriptome analyses in children indicate that humans too cluster into distinct sub-populations, stratifying according to Trim28 expression, transcriptome organization, and obesity-associated imprinted gene dysregulation. These data provide evidence of discrete polyphenism in mouse and man and thus carry important implications for complex trait genetics, evolution, and medicine.

Host and Environmental Factors Influencing Individual Human Cytokine Responses.

Differences in susceptibility to immune-mediated diseases are determined by variability in immune responses. In three studies within the Human Functional Genomics Project, we assessed the effect of environmental and non-genetic host factors of the genetic make-up of the host and of the intestinal microbiome on the cytokine responses in humans. We analyzed the association of these factors with circulating mediators and with six cytokines after stimulation with 19 bacterial, fungal, viral, and non-microbial metabolic stimuli in 534 healthy subjects. In this first study, we show a strong impact of non-genetic host factors (e.g., age and gender) on cytokine production and circulating mediators. Additionally, annual seasonality is found to be an important environmental factor influencing cytokine production. Alpha-1-antitrypsin concentrations partially mediate the seasonality of cytokine responses, whereas the effect of vitamin D levels is limited. The complete dataset has been made publicly available as a comprehensive resource for future studies. PAPERCLIP.

PTER is a N-acetyltaurine hydrolase that regulates feeding and obesity.

Taurine is a conditionally essential micronutrient and one of the most abundant amino acids in humans1-3. In endogenous taurine metabolism, dedicated enzymes are involved in the biosynthesis of taurine from cysteine and in the downstream metabolism of secondary taurine metabolites4,5. One taurine metabolite is N-acetyltaurine6. Levels of N-acetyltaurine are dynamically regulated by stimuli that alter taurine or acetate flux, including endurance exercise7, dietary taurine supplementation8 and alcohol consumption6,9. So far, the identities of the enzymes involved in N-acetyltaurine metabolism, and the potential functions of N-acetyltaurine itself, have remained unknown. Here we show that the body mass index associated orphan enzyme phosphotriesterase-related (PTER)10 is a physiological N-acetyltaurine hydrolase. In vitro, PTER catalyses the hydrolysis of N-acetyltaurine to taurine and acetate. In mice, PTER is expressed in the kidney, liver and brainstem. Genetic ablation of Pter in mice results in complete loss of tissue N-acetyltaurine hydrolysis activity and a systemic increase in N-acetyltaurine levels. After stimuli that increase taurine levels, Pter knockout mice exhibit reduced food intake, resistance to diet-induced obesity and improved glucose homeostasis. Administration of N-acetyltaurine to obese wild-type mice also reduces food intake and body weight in a GFRAL-dependent manner. These data place PTER into a central enzymatic node of secondary taurine metabolism and uncover a role for PTER and N-acetyltaurine in body weight control and energy balance.

Epigenetic inheritance of diet-induced and sperm-borne mitochondrial RNAs.

Spermatozoa harbour a complex and environment-sensitive pool of small non-coding RNAs (sncRNAs)1, which influences offspring development and adult phenotypes1-7. Whether spermatozoa in the epididymis are directly susceptible to environmental cues is not fully understood8. Here we used two distinct paradigms of preconception acute high-fat diet to dissect epididymal versus testicular contributions to the sperm sncRNA pool and offspring health. We show that epididymal spermatozoa, but not developing germ cells, are sensitive to the environment and identify mitochondrial tRNAs (mt-tRNAs) and their fragments (mt-tsRNAs) as sperm-borne factors. In humans, mt-tsRNAs in spermatozoa correlate with body mass index, and paternal overweight at conception doubles offspring obesity risk and compromises metabolic health. Sperm sncRNA sequencing of mice mutant for genes involved in mitochondrial function, and metabolic phenotyping of their wild-type offspring, suggest that the upregulation of mt-tsRNAs is downstream of mitochondrial dysfunction. Single-embryo transcriptomics of genetically hybrid two-cell embryos demonstrated sperm-to-oocyte transfer of mt-tRNAs at fertilization and suggested their involvement in the control of early-embryo transcription. Our study supports the importance of paternal health at conception for offspring metabolism, shows that mt-tRNAs are diet-induced and sperm-borne and demonstrates, in a physiological setting, father-to-offspring transfer of sperm mitochondrial RNAs at fertilization.

Smoking changes adaptive immunity with persistent effects.

Individuals differ widely in their immune responses, with age, sex and genetic factors having major roles in this inherent variability1-6. However, the variables that drive such differences in cytokine secretion-a crucial component of the host response to immune challenges-remain poorly defined. Here we investigated 136 variables and identified smoking, cytomegalovirus latent infection and body mass index as major contributors to variability in cytokine response, with effects of comparable magnitudes with age, sex and genetics. We find that smoking influences both innate and adaptive immune responses. Notably, its effect on innate responses is quickly lost after smoking cessation and is specifically associated with plasma levels of CEACAM6, whereas its effect on adaptive responses persists long after individuals quit smoking and is associated with epigenetic memory. This is supported by the association of the past smoking effect on cytokine responses with DNA methylation at specific signal trans-activators and regulators of metabolism. Our findings identify three novel variables associated with cytokine secretion variability and reveal roles for smoking in the short- and long-term regulation of immune responses. These results have potential clinical implications for the risk of developing infections, cancers or autoimmune diseases.

Human genetics shape the gut microbiome.

Host genetics and the gut microbiome can both influence metabolic phenotypes. However, whether host genetic variation shapes the gut microbiome and interacts with it to affect host phenotype is unclear. Here, we compared microbiotas across >1,000 fecal samples obtained from the TwinsUK population, including 416 twin pairs. We identified many microbial taxa whose abundances were influenced by host genetics. The most heritable taxon, the family Christensenellaceae, formed a co-occurrence network with other heritable Bacteria and with methanogenic Archaea. Furthermore, Christensenellaceae and its partners were enriched in individuals with low body mass index (BMI). An obese-associated microbiome was amended with Christensenella minuta, a cultured member of the Christensenellaceae, and transplanted to germ-free mice. C. minuta amendment reduced weight gain and altered the microbiome of recipient mice. Our findings indicate that host genetics influence the composition of the human gut microbiome and can do so in ways that impact host metabolism.

Mexican Biobank advances population and medical genomics of diverse ancestries.

Latin America continues to be severely underrepresented in genomics research, and fine-scale genetic histories and complex trait architectures remain hidden owing to insufficient data1. To fill this gap, the Mexican Biobank project genotyped 6,057 individuals from 898 rural and urban localities across all 32 states in Mexico at a resolution of 1.8 million genome-wide markers with linked complex trait and disease information creating a valuable nationwide genotype-phenotype database. Here, using ancestry deconvolution and inference of identity-by-descent segments, we inferred ancestral population sizes across Mesoamerican regions over time, unravelling Indigenous, colonial and postcolonial demographic dynamics2-6. We observed variation in runs of homozygosity among genomic regions with different ancestries reflecting distinct demographic histories and, in turn, different distributions of rare deleterious variants. We conducted genome-wide association studies (GWAS) for 22 complex traits and found that several traits are better predicted using the Mexican Biobank GWAS compared to the UK Biobank GWAS7,8. We identified genetic and environmental factors associating with trait variation, such as the length of the genome in runs of homozygosity as a predictor for body mass index, triglycerides, glucose and height. This study provides insights into the genetic histories of individuals in Mexico and dissects their complex trait architectures, both crucial for making precision and preventive medicine initiatives accessible worldwide.

Diminishing benefits of urban living for children and adolescents' growth and development.

Optimal growth and development in childhood and adolescence is crucial for lifelong health and well-being1-6. Here we used data from 2,325 population-based studies, with measurements of height and weight from 71 million participants, to report the height and body-mass index (BMI) of children and adolescents aged 5-19 years on the basis of rural and urban place of residence in 200 countries and territories from 1990 to 2020. In 1990, children and adolescents residing in cities were taller than their rural counterparts in all but a few high-income countries. By 2020, the urban height advantage became smaller in most countries, and in many high-income western countries it reversed into a small urban-based disadvantage. The exception was for boys in most countries in sub-Saharan Africa and in some countries in Oceania, south Asia and the region of central Asia, Middle East and north Africa. In these countries, successive cohorts of boys from rural places either did not gain height or possibly became shorter, and hence fell further behind their urban peers. The difference between the age-standardized mean BMI of children in urban and rural areas was <1.1 kg m-2 in the vast majority of countries. Within this small range, BMI increased slightly more in cities than in rural areas, except in south Asia, sub-Saharan Africa and some countries in central and eastern Europe. Our results show that in much of the world, the growth and developmental advantages of living in cities have diminished in the twenty-first century, whereas in much of sub-Saharan Africa they have amplified.

An orexigenic subnetwork within the human hippocampus.

Only recently have more specific circuit-probing techniques become available to inform previous reports implicating the rodent hippocampus in orexigenic appetitive processing1-4. This function has been reported to be mediated at least in part by lateral hypothalamic inputs, including those involving orexigenic lateral hypothalamic neuropeptides, such as melanin-concentrating hormone5,6. This circuit, however, remains elusive in humans. Here we combine tractography, intracranial electrophysiology, cortico-subcortical evoked potentials, and brain-clearing 3D histology to identify an orexigenic circuit involving the lateral hypothalamus and converging in a hippocampal subregion. We found that low-frequency power is modulated by sweet-fat food cues, and this modulation was specific to the dorsolateral hippocampus. Structural and functional analyses of this circuit in a human cohort exhibiting dysregulated eating behaviour revealed connectivity that was inversely related to body mass index. Collectively, this multimodal approach describes an orexigenic subnetwork within the human hippocampus implicated in obesity and related eating disorders.

The ribosome's comeback: new role in body weight regulation.

The recent discovery of an association between ribosomal DNA (rDNA) copy number and body mass index (BMI) by Law et al. sheds light on a possible role of 45S rDNA in body-weight regulation. This finding opens new avenues for further investigations into the effect of rDNA on various human phenotypes.

An integrative framework to prioritize genes in more than 500 loci associated with body mass index.

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss.

Genetic modifiers of body mass index in individuals with cystic fibrosis.

To identify modifier loci underlying variation in body mass index (BMI) in persons with cystic fibrosis (pwCF), we performed a genome-wide association study (GWAS). Utilizing longitudinal height and weight data, along with demographic information and covariates from 4,393 pwCF, we calculated AvgBMIz representing the average of per-quarter BMI Z scores. The GWAS incorporated 9.8M single nucleotide polymorphisms (SNPs) with a minor allele frequency (MAF) > 0.005 extracted from whole-genome sequencing (WGS) of each study subject. We observed genome-wide significant association with a variant in FTO (FaT mass and Obesity-associated gene; rs28567725; p value = 1.21e-08; MAF = 0.41, ß = 0.106; n = 4,393 individuals) and a variant within ADAMTS5 (A Disintegrin And Metalloproteinase with ThromboSpondin motifs 5; rs162500; p value = 2.11e-10; MAF = 0.005, ß = -0.768; n = 4,085 pancreatic-insufficient individuals). Notably, BMI-associated variants in ADAMTS5 occur on a haplotype that is much more common in African (AFR, MAF = 0.183) than European (EUR, MAF = 0.006) populations (1000 Genomes project). A polygenic risk score (PRS) calculated using 924 SNPs (excluding 17 in FTO) showed significant association with AvgBMIz (p value = 2.2e-16; r2 = 0.03). Association between variants in FTO and the PRS correlation reveals similarities in the genetic architecture of BMI in CF and the general population. Inclusion of Black individuals in whom the single-gene disorder CF is much less common but genomic diversity is greater facilitated detection of association with variants that are in LD with functional SNPs in ADAMTS5. Our results illustrate the importance of population diversity, particularly when attempting to identify variants that manifest only under certain physiologic conditions.

Once-Weekly Semaglutide in Persons with Obesity and Knee Osteoarthritis.

BACKGROUND: Weight reduction has been shown to alleviate symptoms of osteoarthritis of the knee, including pain. The effect of glucagon-like peptide-1 receptor agonists on outcomes in knee osteoarthritis among persons with obesity has not been well studied. METHODS: We conducted a 68-week, double-blind, randomized, placebo-controlled trial at 61 sites in 11 countries. Participants with obesity (a body-mass index [BMI; the weight in kilograms divided by the square of the height in meters] of ≥30) and a clinical and radiologic diagnosis of moderate knee osteoarthritis with at least moderate pain were randomly assigned, in a 2:1 ratio, to receive once-weekly subcutaneous semaglutide (2.4 mg) or placebo, in addition to counseling on physical activity and a reduced-calorie diet. The primary end points were the percentage change in body weight and the change in the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) pain score (on a scale of 0 to 100, with higher scores reflecting worse outcomes) from baseline to week 68. A key confirmatory secondary end point was the physical-function score on the 36-Item Short Form Health Survey (SF-36), version 2 (on a scale of 0 to 100, with higher scores indicating greater well-being). RESULTS: A total of 407 participants were enrolled. The mean age was 56 years, the mean BMI 40.3, and the mean WOMAC pain score 70.9. A total of 81.6% of the participants were women. The mean change in body weight from baseline to week 68 was -13.7% with semaglutide and -3.2% with placebo (P<0.001). The mean change in the WOMAC pain score at week 68 was -41.7 points with semaglutide and -27.5 points with placebo (P<0.001). Participants in the semaglutide group had a greater improvement in SF-36 physical-function score than those in the placebo group (mean change, 12.0 points vs. 6.5 points; P<0.001). The incidence of serious adverse events was similar in the two groups. Adverse events that led to permanent discontinuation of the trial regimen occurred in 6.7% of the participants in the semaglutide group and in 3.0% in the placebo group, with gastrointestinal disorders being the most common reason for discontinuation. CONCLUSIONS: Among participants with obesity and knee osteoarthritis with moderate-to-severe pain, treatment with once-weekly injectable semaglutide resulted in significantly greater reductions in body weight and pain related to knee osteoarthritis than placebo. (Funded by Novo Nordisk; STEP 9 ClinicalTrials.gov number, NCT05064735.).

Polygenic scores for complex traits are associated with changes in concentration of circulating lipid species.

Understanding perturbations in circulating lipid levels that often occur years or decades before clinical symptoms may enhance our understanding of disease mechanisms and provide novel intervention opportunities. Here, we assessed if polygenic scores (PGSs) for complex traits could detect lipid dysfunctions related to the traits and provide new biological insights. We constructed genome-wide PGSs (approximately 1 million genetic variants) for 50 complex traits in 7,169 Finnish individuals with routine clinical lipid profiles and lipidomics measurements (179 lipid species). We identified 678 associations (P < 9.0 × 10-5) involving 26 traits and 142 lipids. Most of these associations were also validated with the actual phenotype measurements where available (89.5% of 181 associations where the trait was available), suggesting that these associations represent early signs of physiological changes of the traits. We detected many known relationships (e.g., PGS for body mass index (BMI) and lysophospholipids, PGS for type 2 diabetes and triacyglycerols) and those that suggested potential target for prevention strategies (e.g., PGS for venous thromboembolism and arachidonic acid). We also found association of PGS for favorable adiposity with increased sphingomyelins levels, suggesting a probable role of sphingomyelins in increased risk for certain disease, e.g., venous thromboembolism as reported previously, in favorable adiposity despite its favorable metabolic effect. Altogether, our study provides a comprehensive characterization of lipidomic alterations in genetic predisposition for a wide range of complex traits. The study also demonstrates potential of PGSs for complex traits to capture early, presymptomatic lipid alterations, highlighting its utility in understanding disease mechanisms and early disease detection.

Mapping the human genetic architecture of COVID-19.

The genetic make-up of an individual contributes to the susceptibility and response to viral infection. Although environmental, clinical and social factors have a role in the chance of exposure to SARS-CoV-2 and the severity of COVID-191,2, host genetics may also be important. Identifying host-specific genetic factors may reveal biological mechanisms of therapeutic relevance and clarify causal relationships of modifiable environmental risk factors for SARS-CoV-2 infection and outcomes. We formed a global network of researchers to investigate the role of human genetics in SARS-CoV-2 infection and COVID-19 severity. Here we describe the results of three genome-wide association meta-analyses that consist of up to 49,562 patients with COVID-19 from 46 studies across 19 countries. We report 13 genome-wide significant loci that are associated with SARS-CoV-2 infection or severe manifestations of COVID-19. Several of these loci correspond to previously documented associations to lung or autoimmune and inflammatory diseases3-7. They also represent potentially actionable mechanisms in response to infection. Mendelian randomization analyses support a causal role for smoking and body-mass index for severe COVID-19 although not for type II diabetes. The identification of novel host genetic factors associated with COVID-19 was made possible by the community of human genetics researchers coming together to prioritize the sharing of data, results, resources and analytical frameworks. This working model of international collaboration underscores what is possible for future genetic discoveries in emerging pandemics, or indeed for any complex human disease.

Collider bias correction for multiple covariates in GWAS using robust multivariable Mendelian randomization.

Genome-wide association studies (GWAS) have identified many genetic loci associated with complex traits and diseases in the past 20 years. Multiple heritable covariates may be added into GWAS regression models to estimate direct effects of genetic variants on a focal trait, or to improve the power by accounting for environmental effects and other sources of trait variations. When one or more covariates are causally affected by both genetic variants and hidden confounders, adjusting for them in GWAS will produce biased estimation of SNP effects, known as collider bias. Several approaches have been developed to correct collider bias through estimating the bias by Mendelian randomization (MR). However, these methods work for only one covariate, some of which utilize MR methods with relatively strong assumptions, both of which may not hold in practice. In this paper, we extend the bias-correction approaches in two aspects: first we derive an analytical expression for the collider bias in the presence of multiple covariates, then we propose estimating the bias using a robust multivariable MR (MVMR) method based on constrained maximum likelihood (called MVMR-cML), allowing the presence of invalid instrumental variables (IVs) and correlated pleiotropy. We also established the estimation consistency and asymptotic normality of the new bias-corrected estimator. We conducted simulations to show that all methods mitigated collider bias under various scenarios. In real data analyses, we applied the methods to two GWAS examples, the first a GWAS of waist-hip ratio with adjustment for only one covariate, body-mass index (BMI), and the second a GWAS of BMI adjusting metabolomic principle components as multiple covariates, illustrating the effectiveness of bias correction.

Body mass index-dependent shifts along large-scale gradients in human cortical organization explain dietary regulatory success.

Making healthy dietary choices is essential for keeping weight within a normal range. Yet many people struggle with dietary self-control despite good intentions. What distinguishes neural processing in those who succeed or fail to implement healthy eating goals? Does this vary by weight status? To examine these questions, we utilized an analytical framework of gradients that characterize systematic spatial patterns of large-scale neural activity, which have the advantage of considering the entire suite of processes subserving self-control and potential regulatory tactics at the whole-brain level. Using an established laboratory food task capturing brain responses in natural and regulatory conditions (N = 123), we demonstrate that regulatory changes of dietary brain states in the gradient space predict individual differences in dietary success. Better regulators required smaller shifts in brain states to achieve larger goal-consistent changes in dietary behaviors, pointing toward efficient network organization. This pattern was most pronounced in individuals with lower weight status (low-BMI, body mass index) but absent in high-BMI individuals. Consistent with prior work, regulatory goals increased activity in frontoparietal brain circuits. However, this shift in brain states alone did not predict variance in dietary success. Instead, regulatory success emerged from combined changes along multiple gradients, showcasing the interplay of different large-scale brain networks subserving dietary control and possible regulatory strategies. Our results provide insights into how the brain might solve the problem of dietary control: Dietary success may be easier for people who adopt modes of large-scale brain activation that do not require significant reconfigurations across contexts and goals.