Association between plausible genetic factors and weight loss from GLP1-RA and bariatric surgery: a multi-ancestry study in 10 960 individuals from 9 biobanks ============================================================================================================================================================== * Jakob German * Mattia Cordioli * Veronica Tozzo * Sarah Urbut * Kadri Arumäe * Roelof A.J. Smit * Jiwoo Lee * Josephine H. Li * Adrian Janucik * Yi Ding * Akintunde Akinkuolie * Henrike Heyne * Andrea Eoli * Chadi Saad * Yasser Al-Sarraj * Rania Abdel-latif * Alexandra Barry * Zhe Wang * Estonian Biobank research team * Pradeep Natarajan * Samuli Ripatti * Anthony Philippakis * Lukasz Szczerbinski * Bogdan Pasaniuc * Hamdi Mbarek * Ruth J.F. Loos * Uku Vainik * Andrea Ganna ## Abstract Obesity is a significant public health concern. GLP-1 receptor agonists (GLP1-RA), predominantly in use as a type 2 diabetes treatment, are a promising pharmacological approach for weight loss, while bariatric surgery (BS) remains a durable, but invasive, intervention. Despite observed heterogeneity in weight loss effects, the genetic effects on weight loss from GLP1-RA and BS have not been extensively explored in large sample sizes, and most studies have focused on differences in race and ethnicity, rather than genetic ancestry. We studied whether genetic factors, previously shown to affect body weight, impact weight loss due to GLP1-RA therapy or BS in 10,960 individuals from 9 multi-ancestry biobank studies in 6 countries. The average weight change between 6 and 12 months from therapy initiation was -3.93% for GLP1-RA users, with marginal differences across genetic ancestries. For BS patients the weight change between 6 and 48 months from the operation was -21.17%. There were no significant associations between weight loss due to GLP1-RA and polygenic scores for BMI or type 2 diabetes or specific missense variants in the *GLP1R, PCSK1 and APOE* genes, after multiple-testing correction. However, a higher polygenic score for BMI was significantly linked to lower weight loss after BS (+0.7% for 1 standard deviation change in the polygenic score, P = 1.24×10-4). In contrast, higher weight at baseline was associated with greater weight loss. Our findings suggest that existing polygenic scores related to weight and type 2 diabetes and missense variants in the drug target gene do not have a large impact on GLP1-RA effectiveness. Our results also confirm the effectiveness of these treatments across all major continental ancestry groups considered. ## Introduction The obesity epidemic presents a significant public health burden, driving the need for effective treatment options.1,2 Despite 75% of adults with obesity having attempted to lose weight, most have not achieved lasting success.3 Amid these challenges, glucagon-like peptide-1 receptor agonists (GLP1-RA) have emerged as a promising solution. As first discovered in the late 1980s, glucagon-like peptide-1 (GLP-1) is linked to the incretin effect, which promotes insulin release from pancreatic beta cells, thus playing a crucial role in regulating blood sugar levels.4,5 This evidence led to the approval of the first GLP1-RA for treatment of type 2 diabetes in 2005.6–8 Nine years later, an association between GLP-1 and reduced food intake paved the way for the development and approval of the first GLP1-RA for obesity treatment.5,9 Clinical trials conducted in the following years demonstrated weight loss effects ranging from 6% to 16% for liraglutide and semaglutide, leading to a substantial increase in GLP1-RA usage worldwide10–12. Additional trials have further revealed cardiovascular13,14 and renal benefits15, with ongoing trials currently underway.16 Overall, anticipated trends suggest that by 2030, the US will witness approximately 30 million patients using GLP1-RA drugs, representing roughly 9% of its population, underscoring their potential in combating the obesity crisis and exploring new therapeutic applications.17–19 Meanwhile, bariatric surgery (BS) stands as a durable and effective treatment for individuals with severe obesity (BMI >= 40 kg/m2) and metabolic comorbidities, and consistent with its effectiveness, the threshold has been reduced to a BMI >=30 kg/m2.20,21 Meta-analyses combining results from randomized controlled trials and observational studies found that, on average, individuals undergoing BS experienced a weight loss of 26 kg.22,23 However, such an invasive procedure inherently carries increased risks of severe adverse events.24,25 Despite the robust weight loss effects observed with both GLP1-RA treatment and BS, there exists heterogeneity in their effectiveness. One out of every five patients undergoing BS may not achieve the desired weight loss within the first year or may experience weight regain within two years.26–29 The variability of GLP1-RA treatment effects on glycemic outcomes has been widely investigated and can be partially explained by clinical indicators like low β-cell function or C-peptide levels as well as genetic variants.30–35 However, with the exception of type 2 diabetes, sex, and baseline weight, the factors influencing the diversity in weight loss effects of these drugs are still poorly understood.36–38 Additionally, differences in patient adherence to treatment regimens may further contribute to this variability.39 Genetics has been considered a potential factor for heterogeneity in weight loss treatments, but its role remains largely unexplored.36 A previous study has examined the role of a common missense variant in the *GLP1R* genes in 57 women and found a weak link with GLP1-RA treatment response.32,40 A greater number and better-powered studies have been conducted in individuals with BS 41,42, however, a consistent genome-wide association for weight loss has yet to be identified. Finally, there has been limited exploration into the effectiveness of these obesity treatments across different genetic ancestries and across various healthcare settings.10,43 In our study we seek to i) characterize the observed GLP1-RA and BS body weight-lowering effects using real-world data from 10 960 individuals across 9 biobanks in 6 countries, with 6 major continental ancestry groups; ii) identify whether plausible genetic factors (i.e. polygenic scores for BMI and T2D and coding variants in *GLP1R*, *PCSK1* and *APOE* gene) associate with heterogeneity in GLP1-RA and BS weight loss effects, and thus identify patient groups that will benefit most from such treatments; iii) compare the effects of GLP1-RA treatment and BS on weight loss across major continental ancestry groups. ## Results ### GLP1-RA and bariatric surgery associated weight loss in 10 960 individuals from 9 cohorts and 6 major continental ancestry groups We investigated changes in body weight associated with GLP1-RA treatment and BS across 9 biobank studies: Helsinki University Hospital (HUS, Finland; N = 633), Estonian Biobank (ESTBB; N = 464), UK Biobank (UKBB, United Kingdom; N = 810), All of Us (AoU, USA; N = 559), BioMe Biobank (BioMe, USA; N = 2,170), Mass General Brigham Biobank (MGBB, USA; N = 2,141), Atlas Biobank (UCLA-ATLAS, USA; N = 1,445), Qatar Biobank (QBB, Qatar; N = 2,383), Bialystok Bariatric Surgery Study (BBSS, Poland, N = 355). We included individuals that were at least 18 years of age when starting a GLP1-RA treatment or undergoing BS. Individuals undergo GLP1-RA treatment for different medical reasons, not only for weight loss. We defined GLP1-RA usage based on prescription or purchase data (ATC codes starting with A10B*; including exenatide, liraglutide, lixisenatide, albiglutide, dulaglutide, semaglutide, and beinaglutide) covering at least 12 months, and excluding individuals with indication of treatment discontinuation within this time window. For BS, we considered individuals undergoing the following procedures: Roux-en-Y gastric bypass, sleeve gastrectomy, adjustable gastric band, vertical banded gastroplasty, biliopancreatic diversion with duodenal switch. More information can be found in the **Methods** section. We included only individuals who had an initial body weight measurement recorded at most 12 months prior to the baseline (defined as treatment initiation or surgery) and a follow-up measurement taken at least 6 but no more than 12 months after the baseline for GLP1-RA treatments, or within 48 months after baseline for BS. In case of multiple measurements within this time frame, we considered the median of such measurements (**Supplementary Figure 1** for a schematic of the approach used). We defined our main outcome of interest as the percentage change in body weight, calculated as the difference between the second body weight measurement (or median of multiple measurements) and the initial one, divided by the initial measurement. Our study population included a total of 6,750 GLP1-RA users and 4,210 individuals undergoing BS (baseline characteristics are reported in **Table 1**). Overall, GLP1-RA users 1) had a lower proportion of women (61.0% weighted average across studies) when compared to BS patients (73.1% weighted average across studies), 2) were on average older (57.30 *vs* 46.50), and 3) had a lower initial body weight (99.80 kg *vs* 116.96 kg). The proportion of individuals diagnosed with type 2 diabetes at time of GLP1-RA initiation was at least 48% across all studies and major continental ancestry groups. View this table: [Table 1](http://medrxiv.org/content/early/2024/09/12/2024.09.11.24313458/T1) Table 1 Baseline characteristics for 6 750 GLP1-RA users and 4 210 patients undergoing bariatric surgery included in our study population. For each cohort and major continental ancestry group (AFR: African, AMR: admixed American, EAS: East Asian, EUR: European, MID: Middle Eastern, SAS: South Asian), we report the total number of individuals eligible for analysis, the proportion of females, the mean (SD) body weight (in Kg) at baseline, the mean (SD) age (in years) at baseline, the prevalence of T2D, and (for GLP1-RA users only) the most commonly used type of medication and the proportion of semaglutide users. HUS, Helsinki University Hospital; ESTBB, Estonian Biobank; UKBB, UK Biobank; AoU, All of US; MGBB, Mass General Brigham Biobank; BioMe, BioMe; UCLA-ATLAS, Atlas Biobank; QBB, Qatar Biobank; BBSS, Bialystok Bariatric Surgery Study Among GLP1-RA users, we observed an average body weight change (**Figure 1**) across studies of -3.93% (ranging from -1.08% to -7.07%), in line with results from observational studies.44,45 Among BS patients, we observed an average change of -21.17% (-15.00% to -27.72%). ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/12/2024.09.11.24313458/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2024/09/12/2024.09.11.24313458/F1) Figure 1 Average percentage change in body weight for 6 750 GLP1-RA users and 4 210 bariatric surgery patients. Bars represent the mean percentage change in body weight within each study and ancestry, and the overall combined weighted mean change. Error bars represent standard errors. ### Demographic factors associated with GLP1-RA and bariatric surgery weight change We explored the effect of baseline characteristics on body weight changes associated with GLP1-RA treatment and BS. We fitted a multivariable linear model with percentage change in body weight as outcome and baseline body weight, sex, age at baseline, the first 20 principal genetic components and medication type (for GLP1-RA users only) as predictors. Each cohort was analyzed separately, stratifying by major continental ancestry group, and the effects were then meta-analyzed (**Table 2**). Baseline body weight was significantly associated with greater weight loss for both GLP1-RA and BS treatments (*βGLP1-RA* = -0.05 % weight change compared to baseline, *P* = 1.63 × 10-35; *βBS* = -0.14 % weight change compared to baseline, *P* = 5.49 × 10-69) and women had a significantly greater weight loss for both treatments (*βGLP1-RA* = -1.54 % weight change compared to baseline, *P* = 2.47 × 10-16; *βBS* = -2.66 % weight change compared to baseline, *P* = 6.83 × 10-10). Being older at treatment initiation was only associated with less weight loss for bariatric surgery (*βBS* = 0.15 % weight change compared to baseline, *P* = 2.54 × 10-15). The effects of baseline weight and sex on weight loss were highly heterogeneous across ancestries and studies (I2 for GLP1-RA= 0.99 and 0.83 for baseline weight and sex, respectively) (**Supplementary Figures 2-4 and Supplementary Tables 1-2**). View this table: [Table 2](http://medrxiv.org/content/early/2024/09/12/2024.09.11.24313458/T2) Table 2 Effect of baseline weight, sex and age at baseline of percentage weight changes associated with GLP1-RA and bariatric surgery. Multivariable model adjusting for the three variables described in the table and the first 20 principal genetic components and medication type (for GLP1-RA users only). Negative effects indicate larger weight loss. ### Marginal effect of ancestry on GLP1-RA and bariatric surgery associated weight loss We evaluate whether the effect of GLP1-RA and BS associated weight loss was significantly different across genetically-defined continental ancestry groups after accounting for differences in weight at baseline, age, sex and, only for GLP1-RA, medication type (**Table 3**). Only multi-ancestry biobanks (AoU, BioMe, MGBB, UCLA-ATLAS) were included in this analysis (N = 6 315). Compared to individuals of European ancestry, all other ancestry groups had lower change in body weight. However, this effect was statistically significant only among GLP1-RA users of African (*β* = 0.76 % weight change from baseline weight, *P* = 0.02) and admixed American (*β* = 0.87 % weight change from baseline weight, *P* = 3 × 10-3) ancestries. View this table: [Table 3](http://medrxiv.org/content/early/2024/09/12/2024.09.11.24313458/T3) Table 3 Associations between weight change and major continental ancestry groups. Meta-analysis effect sizes, standard errors and P values for the association between weight change and major continental ancestry groups (AFR: African, AMR: admixed American, EAS: East Asian, EUR: European, SAS: South Asian). Coefficients are derived from a multivariable model including percentage weight change as outcome and ancestry as categorical variable (with EUR as reference category), and adjusting for sex, age and weight at baseline and, for GLP1-RA only, medication type. Negative effects indicate larger weight loss ### Plausible genetic factors do not associate with weight change due to GLP1-RA We investigated the effect of genetic exposures with plausible effects on body weight changes. Genetic exposures were selected upon discussion between authors and before initiating analyses. First, we considered polygenic scores for BMI based on previous observation that a polygenic score was associated with BMI trajectories.46,47 We also consider the only replicated genome-wide variant associated with BMI change in the largest study of BMI trajectories (*rs429358* in *APOE*).46 Second, we considered a polygenic score for T2D based on the observation that biomarkers of type 2 diabetes were associated with effectiveness of GLP1-RA.31,33 We confirmed that both polygenic scores for BMI and T2Ds were significantly associated with their respective traits across all major continental ancestry groups (**Supplementary Table 3**) although the effects were larger in individuals of European ancestry. Third, we considered 13 SNPs in the *GLP1R* and *PCSK1* genes. We considered all missense variants in the *GLP1R* gene with at least 1% minor allele frequency in at least one of major continental ancestry group, as per gnomAD (version 4.1)48, three of which have been functional characterized as gain or loss of function (**Supplementary Table 4**).49 We also considered two non-synonymous variants in the *PCSK1* gene that have been associated with BMI variation in the population.50 Finally, we sought to replicate the effect of two variants associated with weight loss in one of the largest studies of BS42(**Methods**). After analyses, two missense variants in the *GLP1R* genes, *rs2295006* and *rs201672448* were too rare to be tested across studies and were not considered further. We estimated the effect of each genetic exposure on weight change in each cohort and major continental ancestry groups separately, fitting a linear model adjusted for baseline body weight, sex, age at baseline, first 20 principal genetic components, and medication type (for GLP1-RA users only), and we subsequently combined the effects through ancestry-specific and multi-ancestry meta-analyses (**Figure 2, Supplementary Tables 5-6**). As the approach of adjusting longitudinal change phenotypes for the baseline trait has been criticized for potentially inducing biases for genetic variants associated with the baseline trait51, we also considered a second model using post-treatment weight as outcome, including the same predictors as in the model above and an additional term modeling the interaction between baseline weight and each genetic exposure (**Supplementary Tables 7-8**). ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/12/2024.09.11.24313458/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2024/09/12/2024.09.11.24313458/F2) Figure 2 Effect of 15 genetic exposures on body weight changes associated with GLP1-RA treatment and bariatric surgery. a. Ancestry-specific and multi-ancestry meta-analysis effect sizes for association between percentage change in body weight and PGS for BMI and type 2 diabetes. Dots represent the percentage change in body weight per one standard deviation change in PGS, error bars represent the 95% confidence interval. Full dots represent statistical significance at P < 0.004. b. Multi-ancestry meta-analysis effect sizes for the association between percentage change in body weight and genotype at each locus. Dots represent the effect size, error bars represent the 95% confidence interval. PGS, polygenic score. Among GLP1-RA users, we did not observe any significant genetic exposure associated with weight loss after multiple testing correction (*P* < 0.004, Bonferroni correction for 13 exposures tested). Using the alternative model, we identified one missense variant in the *GLP1R* genes, *rs3765467*, for which we observed a statistically significant interaction with baseline weight (*βw0*rs3765467* = 0.09 kg weight change from baseline weight, *P* = 2.7 × 10-4) on the weight measured post GLP1-RA treatment after multiple testing correction. When considering nominal significance, another missense variant in the *GLP1R* genes*, rs6923761*, had a statistically significant interaction with baseline weight (*βw0*rs6923761* = 0.02 kg weight change from baseline weight, *P* = 0.01) on the weight measured post GLP1-RA treatment. We further tested whether the lack of associations observed was possibly due to high effect-size heterogeneity across ancestries, or due to lack of statistical power. We did not observe any significant heterogeneity across ancestries for any of the genetic exposure (average I2 across genetic exposures for GLP1-RA was 0.06 **- Supplementary Tables 9-10**). Moreover, given the current sample size for GLP1-RA users (N = 6 750) our power calculation showed we would have at least 80% power to detect a statistically significant effect (at *P* < 0.05) of at least a 0.3% change in body weight for 1 standard deviation (SD) change in the polygenic score, or for SNPs with a minor allele frequency of at least 1% (**Methods**). ### A polygenic score for BMI associates with weight change due to BS Among BS patients, a higher polygenic score for BMI was significantly associated with lower weight loss (*βBMI PGS* = 0.70 % weight change compared to baseline for 1 SD change in the polygenic score, *P* = 1.24 × 10-4) after multiple testing correction (*P* < 0.01, Bonferroni correction for 5 exposures tested), in the opposite direction than observed for weight at baseline (**Figure 2**). A higher polygenic score for T2D was nominally associated with lower weight loss (*βT2D PGS* = 0.45 % weight change compared to baseline for 1 SD change in the polygenic score, *P* = 0.03). We observed significant heterogeneity in the association between BMI PGS and weight loss across studies (I2 = 0.64, **Supplementary Table 10**). ## Discussion In this study we provide a comprehensive assessment of major demographic factors and plausible genetic factors on body weight changes associated with GLP1-RA use and BS. In line with previous work,36,52–55 our main outcome of interest was the percentage change in body weight between 6 and 12 months after GLP1-RA treatment initiation as well as between 6 and 48 months after BS. By considering studies from 6 different countries, we were able to highlight the heterogeneity in use of GLP1-RA across different healthcare systems. GLP1-RA users consisted mostly of individuals with type 2 diabetes, reflecting their original use for treatment of this condition. Biobank studies covering Boston (MGBB), New York (BioMe), Los Angeles (UCLA-ATLAS) and Helsinki (HUS) had overall lower rates of type 2 diabetes individuals compared to the other studies and lower weight at baseline, which, together with the availability of more updated electronic health records data, suggest a larger number of individuals being prescribed GLP1-RA for weight loss. However, care should be taken in interpreting variation in rate of type 2 diabetes across biobanks as this might reflect different approaches used for capturing this diagnosis from the electronic health records as well as different sampling strategies of biobanks. The comparison between GLP1-RA and BS is informative as it highlights the younger age, higher weight and lower percentage of men undergoing BS compared to GLP1-RA users. The weight reduction was also approximately 5 times higher among BS patients compared to GLP1-RA users. The weight reduction among GLP1-RA users was also more heterogeneous across biobanks compared to BS, but in line with what reported previously in the literature using real world data34 and lower than results from clinical trials.11,56 Our findings also highlight a consistent effect of baseline body weight and sex on weight loss associated with both GLP1-RA treatment and BS, suggesting that individuals with higher weight at baseline, as well as women benefit more from both interventions. These findings underscore the importance of considering individual characteristics when designing weight loss interventions. The main aim of this study was to evaluate the impact of genetic exposures on body weight changes induced by GLP1-RA or BS. At the current sample size, a genome-wide scan would have limited power to identify potential new signals. With 6,750 individuals, we would have at least 80% statistical power to identify, at a genome-wide significant threshold, a genetic variant with an allele frequency of at least 1% and an effect larger than 0.3% weight change per change in one SD, which is unlikely given previous effects of common variants on treatment response.51 Instead, we formulated a set of hypotheses based on the results from the largest genome-wide association studies of BMI and type 2 diabetes, and, explored biologically plausible mutations in several candidate genes. Polygenic scores represent the most powerful genetic predictors for BMI and type 2 diabetes and we confirm their significant association with the corresponding trait across all major continental ancestry groups. Thus, lack of associations with these polygenic scores would suggest that GLP1-RA and BS-related weight loss is independent of one’s genetic predisposition to BMI or type 2 diabetes, as explained by common genetic variants. Moreover, because the biological underpinning of weight change is partially distinct from weight at baseline, we also consider the few replicated genetic signals for BMI variability and BMI trajectories. For bariatric surgery we observed that genetically higher BMI was associated with weight gain, in the opposite direction for the phenotypic association with initial body weight. This might suggest that individuals with a higher genetic predisposition for elevated BMI might regain weight, after the operation, faster than those with lower genetic predisposition. To maintain weight loss, they may need additional lifestyle modifications, behavioral health therapy or adjuvant pharmacologic therapies as a GLP1-RA treatment.57,58 On the contrary, none of the genetic exposures were associated with weight loss among GLP1-RA users, after multiple-testing correction. These results were not attributable to high heterogeneity across studies or across major continental ancestry groups. Whether we directly tested the impact of genetic exposure on percentage body weight change or used post-treatment weight as outcome and tested the interaction between baseline weight and the genetic exposure, results did not significantly change. In light of the associations observed for BS and following the hypothesis that individuals with genetically higher BMI tend to re-gain weight faster, it would not be surprising to find, at larger sample size, a significant association between BMI polygenic score and weight loss due to GLP1-RA. However, our analyses exclude that this effect, should it exist, is unlikely to guide clinical treatment assignment. Our results showed that the considered genetic influence on weight loss in response to GLP1-RA treatment may be too small to impact clinical decision-making or patient outcomes. This statement should be revisited either by using more powerful polygenic scores or by considering individuals undergoing GLP1-RA treatment only for weight loss. Another aim of the study was to compare the effects of GLP1-RA treatment and BS on weight loss across major continental ancestry groups. Ancestry can provide more useful information about population health than racial and ethnic categories because it is not directly influenced by geographic, cultural, and sociopolitical forces.59 Nonetheless, ancestry is also partially correlated with self-reported race and ethnicity and socio-economic characteristics do differ across racial and ethnic groups. Overall, we did not find large differences in GLP1-RA and BS associated weight loss between major continental ancestry groups. We do observe that GLP1-RA users of African and admixed American ancestry had statistically significantly lower change in body weight compared to Europeans, after accounting for differences in sex, age, weight at baseline and medication type. The effects were consistent in BioMe and UCLA-ATLAS biobanks for admixed American ancestry, and additionally in AoU and MGBB for African ancestry. This might reflect different socio-economic characteristics between individuals of admixed American and European ancestry resulting, for example, in differences in access to treatment. The main strength of our study is its size and diversity across countries, healthcare systems and ancestries. Our study presents various limitations. First, we cannot be sure that all GLP1-RA users fulfilled their prescription. This is a potential issue only in studies with prescription information but not in studies with purchase information. Nonetheless weight loss was comparable across both study types. Moreover, we ensured that observed changes in body weight are due to the treatment by considering only measurements within one month from the last purchase or prescription. Second, we could not accurately estimate treatment adherence, but overcome this limitation by considering treatments of maximum 12 months, a period for which adherence has been observed to be around 65%.52 Third, we used ancestry labels based on major continental ancestry groups, while ancestry can be better characterized as a continuous quantity.59 This discretization reflects pragmatic considerations to enable comparable analyses across the different studies. In conclusion, our study suggests that GLP1-RA treatments work equally well in individuals carrying common non-synonymous mutations in the GLP1R gene and in individuals at high genetic risk for BMI and type 2 diabetes. It also suggests that BS is less effective among individuals with higher polygenic score for BMI, but the effect is not large enough to be clinically relevant. In fact, for each additional SD in the polygenic score for BMI, individuals undergoing BS have on average 0.7% lower weight loss from their initial body weight. Finally, these two major weight-loss interventions achieved sustained weight reduction regardless of an individual’s genetic ancestry. However, socio-economic characteristics correlating with ancestry should be further studied to better understand some of the remaining differences observed in this study. ## Methods ### Study population In the current study, we included samples from 10 960 individuals from the following 9 biobanks: Helsinki University Hospital (HUS, Finland; N = 633), Estonian Biobank (ESTBB, Estonia; N = 464), UK Biobank (UKBB, United Kingdom; N = 810), All of Us (AoU, USA; N = 559), BioMe (BioMe, USA; N = 2,170), Mass General Brigham Biobank (MGBB, USA; N = 2,141), Atlas Biobank (UCLA-ATLAS, USA; N = 1,445), Qatar Biobank (QBB, Qatar; N = 2,383), Bialystok Bariatric Surgery Study (BBSS, Poland, N = 355). The biobank studies include samples from (hospital) biobanks, prospective epidemiological and disease-based cohorts. Follow-up covers a total of 20 years with the earliest study starting follow-up in 2004 (AoU) and the latest study ending follow-up in 2024 (BioMe). In the **Supplementary Material** we provide a detailed description of the population selected in each study. ### Inclusion and exclusion criteria for GLP1-RA analysis For our GLP1-RA analysis we implemented the following inclusion and exclusion criteria. We required an initiation of a GLP1-RA treatment (ATC codes: A10BJ*, A10BX04, A10BX10, A10BX13, A10BX14) defined as the date of first medication purchase or medication prescription of GLP1-RA. Individuals were required to be on a GLP1-RA treatment for at least twelve months. This was defined through prescription length, regular drug prescription refills or drug purchases, depending on the dataset and healthcare setting. Individuals aged under 18 were excluded from the analysis. Furthermore, individuals were required to have at least one body weight measurement maximum one year prior to treatment initiation or within 14 days after initiation, if no measure before initiation was available. If we observed multiple weight measurements within this time frame, only the closest one to initiation was considered. This measurement defined our baseline weight variable. At least one body weight measurement between 26 weeks and 52 weeks was necessary for each individual to be included in the analysis. Individuals were excluded if they underwent BS prior to or during the first year of GLP1-RA treatment. ### Inclusion and exclusion criteria for BS analysis For our BS analysis we implemented the following inclusion and exclusion criteria. Individuals were included if they underwent any type of BS. The respected code definitions for BS varied across countries and healthcare systems and can be found in **Supplementary Table 11**. Individuals aged under 18 at baseline were excluded from the analysis. Analogously to the GLP1-RA analysis, individuals were required to have at least one body weight measurement maximum one year prior to surgery or within 14 days after it, if no measure before baseline was available. If we observed multiple weight measurements within this time frame, only the closest one to initiation was considered. This measurement defined our baseline weight variable. Additionally, at least one body weight measurement between 26 weeks and 208 weeks was necessary for each individual to be included in the analysis. ### Outcome definitions In our primary model, for both GLP1-RA and BS we defined the outcome to be the percentage change in body weight from baseline The baseline body weight measured closest to *T* and between -52 weeks and +2 weeks from *T* was defined as *W*. The second body weight (*mW1*) was defined as the median of body weight measurements between 26 weeks and 52 weeks from *T* for the GLP1-RA analysis and between 180 days and 1460 days for the BS analysis. This was to enhance robustness against outlier measurements occurring within these intervals. Therefore, the percentage change in body weight from baseline was calculated using ![Formula][1] In our secondary model, we used *mW1* as an outcome. ### Genetic exposures We identified 15 genetic variants (single nucleotide polymorphism, SNP) and two PGS of interest for our analysis. Among the included SNPs 3 are functionally characterized *GLP1R* variants (*rs10305492*, *rs146868158*, *rs6923761*) associated with random glucose levels that have shown decreased/increased response to different endogenous and exogenous GLP1-RA.49 Seven SNPs are *GLP1R* missense variants (*rs1042044*, *rs10305420*, *rs3765467*, *rs10305421*, *rs2295006*, *rs10305510*, *rs201672448*) with a frequency over 1% in at least one ancestry group in gnomAD (version 4.1).48 One SNP is the only robustly replicated lead genome-wide significant variant (*rs429358*, *ApoE*) in a GWAS of BMI change over time.46,60 Two SNPs are *PCSK1* missense variants (*rs6232*, *rs6235*) that have been found associated with BMI variations.50 And two SNPs (*rs728996*, *rs17702901*) were associated with an effect on excess body mass index loss amongst BS patients.42 A list of all genetic exposures and their frequencies across all major continental ancestry groups can be found in **Supplementary Table 4**. The included PGS are for body mass index (BMI) and type 2 diabetes (T2D).61 The PGS weights were taken from Weissbrod et al. 202261 and obtained from UK Biobank. The scores were computed in each study using PLINK 2.0.62 Because the PGS scores were derived from UK Biobank, we couldn’t use the same score in UK Biobank, instead we used Thompson et al. 2022.63 ### Statistical models We employed two linear regression models (A and B) to investigate the effect of the chosen genetic exposures on weight change after initiation of GLP1-RA treatment and BS, respectively. In both models we adjusted for the baseline weight, sex, age at initiation of treatment or surgery in years, the first 20 genetic principal components and study-specific covariates, such as the genotyping batch. In the GLP1-RA analysis we additionally adjusted for the medication type within the drug class. In model A % *weight change* was the outcome variable, while in model B the outcome was *mW1*. In model B we additionally included the interaction term between the genetic exposure variable and the baseline weight. #### Model A (primary model) ![Formula][2] #### Model B (secondary model) ![Formula][3] ### Meta-analyses We conducted ancestry-specific and multi-ancestry fixed-effect meta-analyses for each of the genetic exposures. Effect sizes were combined using fixed-effect inverse-variance-weighting (IVW), as implemented in the R package “meta”.64 For each exposure we reported the unadjusted *P* values, but considered as threshold for statistical significance *P* < 0.004 for GLP1-RA (i.e. P < 0.05/13, Bonferroni correction for 13 exposures tested) and *P* < 0.01 for BS (i.e. P < 0.05/5, Bonferroni correction for 5 exposures tested). We tested for heterogeneity in effect sizes between ancestries using Cochran’s Q-test65 and by inspecting the I2 statistics66, which is defined as: ![Formula][4] where *df* is equal to *number of studies* − 1. ### Power calculation To calculate the statistical power for the PGS effect sizes, we used a t-test for linear regression coefficient as implemented in the R package pwrs67 assuming a standard deviation for the outcome variable SDy = 7.8 (weighted combine SD across studies), a sample size of n = 6 750, a total number of predictors k = 25 and an adjusted r2 = 0.05, to calculate power at P < 0.05 significance level. For the coefficients of association with single variants, we estimated statistical power via the non-centrality parameter (NCP) of the chi-square distribution. We defined NCP = 2 f (1 – f) n β^2 where f is MAF, n is the effective sample size, and β is the expected effect size.68 We set a beta of 0.3 (% weight change from baseline) as the expected genetic effect on GLP1 response, based on findings from a recent large-scale study on genetic influences on drug response.51 ## Supporting information Supplementary Material [[supplements/313458_file02.pdf]](pending:yes) Supplementary Tables [[supplements/313458_file03.pdf]](pending:yes) ## Data Availability All data utilized in this study can be accessed by applying to each of the participating studies individually. Additional details are provided in the manuscript under the Data Availability section. ## Code Availability The study utilized previously published analysis tools as described in the **Methods** section. Additional code used for these analyses is available at [https://github.com/dsgelab/glp1-bs-genetics/](https://github.com/dsgelab/glp1-bs-genetics/) ## Data Availability The Helsinki biobank can provide access for research projects within the scope regulated by the Finnish Biobank Act, which is research utilizing the biobank samples or data for the purposes of promoting health, understanding the mechanisms of disease or developing products and treatment practices used in health and medical care. De-identified data of the MGBB that supports this study is available from the MGB Biobankportal at [https://www.massgeneralbrigham.org/en/research-and-innovation/participate-in-research/biobank/for-researchers](https://www.massgeneralbrigham.org/en/research-and-innovation/participate-in-research/biobank/for-researchers). Restrictions apply to the availability of these data, which are available to MGB-affiliated researchers via a formal application. UK Biobank data are available through a procedure described at [http://www.ukbiobank.ac.uk/using-the-resource/](http://www.ukbiobank.ac.uk/using-the-resource/). Registered researchers whose institutions have Data Use and Registration Agreements in place with All of Us that include the Controlled Tier can access genomic data. Clinical and genotype data from UCLA ATLAS Community Health Initiative patients, de-identified for research purposes, are accessible to UCLA-approved researchers through the Discovery Data Repository (DDR). BioMe data is available through a process described at [https://icahn.mssm.edu/research/ipm/programs/biome-biobank/researcher-faqs](https://icahn.mssm.edu/research/ipm/programs/biome-biobank/researcher-faqs). Pseudonymised data and/or biological samples can be accessed for research and development purposes in accordance with the Estonian Human Genome Research Act. To access data, the research proposal must be approved by the Scientific Advisory Committee of the Estonian Biobank as well as by the Estonian Committee on Bioethics and Human Research. For more details on data access and relevant documents, please see [https://genomics.ut.ee/en/content/estonian-biobank#dataaccess](https://genomics.ut.ee/en/content/estonian-biobank#dataaccess). For access to the Bialystok Bariatric Surgery Study please refer to the Medical University of Bialystok. The data and biosamples collected or generated by QBB will be made available to researchers employed within or otherwise contractually bound to public and private institutions that conduct scientific research and that meet the requirements detailed in the Qatar Biobank Research Access policy. The PGS weights used in this study are available in Weissbrod et al. 2022 and Thompson et al. 2022.61 ## Acknowledgements We gratefully acknowledge the resources provided by the Institute for Precision Health and participating patients from the UCLA ATLAS Community Health Initiative ([https://www.uclahealth.org/international-services/consulting-services/clinical-research-consultation/ucla-atlas-community-health-initiative](https://www.uclahealth.org/international-services/consulting-services/clinical-research-consultation/ucla-atlas-community-health-initiative)). The UCLA ATLAS Community Health Initiative in collaboration with UCLA ATLAS Precision Health Biobank is a program of the Institute for Precision Health, which directs and supports the biobanking and genotyping of biospecimen samples from participating patients from UCLA in collaboration with the David Geffen School of Medicine, UCLA Clinical and Translational Science Institute and UCLA Health. The ATLAS Community Health Initiative is supported by UCLA Health, the David Geffen School of Medicine and a grant from the UCLA Clinical and Translational Science Institute (UL1TR001881). This work was financially supported in part by National Institutes of Health awards U01HG011715, R01HG009120 and R01MH115676. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. We would furthermore like to thank all participants of HUS ([https://www.helsinginbiopankki.fi/en/Helsinki-Biobank](https://www.helsinginbiopankki.fi/en/Helsinki-Biobank); [https://www.hus.fi/en](https://www.hus.fi/en)), ESTBB ([https://genomics.ut.ee/en/content/estonian-biobank](https://genomics.ut.ee/en/content/estonian-biobank)), BBSS ([https://classic.clinicaltrials.gov/ct2/show/NCT04634591](https://classic.clinicaltrials.gov/ct2/show/NCT04634591)), QBB ([https://www.qphi.org.qa/](https://www.qphi.org.qa/)), BioMe ([https://icahn.mssm.edu/research/ipm/programs/biome-biobank](https://icahn.mssm.edu/research/ipm/programs/biome-biobank)), UKBB ([https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)), MGBB ([https://www.massgeneralbrigham.org/en/research-and-innovation/participate-in-research/biobank](https://www.massgeneralbrigham.org/en/research-and-innovation/participate-in-research/biobank)) and AOU ([https://allofus.nih.gov/](https://allofus.nih.gov/)) for their generous contribution. This work was further supported in part by funding from the Eric and Wendy Schmidt Center at the Broad Institute of MIT and Harvard. U.V. and K.A. have been funded by Estonian Research Council’s personal research funding start-up grant PSG759 and European Research Council’ starting grant under the grant agreement no 101117251 (OBECAUSE). J.H.L. is supported by NIDDK K23 DK131345. The activities of the EstBB are regulated by the Human Genes Research Act, which was adopted in 2000 specifically for the operations of the EstBB. Individual level data analysis in the EstBB was carried out under ethical approvals of 1.1-12/1409 and 1.1-12/2161 from the Estonian Committee on Bioethics and Human Research (Estonian Ministry of Social Affairs), using data according to release application 6-7/GI/18857 from the Estonian Biobank. ## Footnotes * + These authors jointly supervised the project * Received September 11, 2024. * Revision received September 11, 2024. * Accepted September 12, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Hruby, A. & Hu, F. B. The Epidemiology of Obesity: A Big Picture. Pharmacoeconomics 33, 673–689 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40273-014-0243-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25471927&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 2. 2.Boutari, C. & Mantzoros, C. S. A 2022 update on the epidemiology of obesity and a call to action: as its twin COVID-19 pandemic appears to be receding, the obesity and dysmetabolism pandemic continues to rage on. Metabolism 133, 155217 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.metabol.2022.155217&link_type=DOI) 3. 3.EASO Secretariat. European survey finds over three-quarters of adults with obesity have attempted to lose weight in the past year, but most have been unsuccessful. EASO [https://easo.org/european-survey-finds-over-three-quarters-of-adults-with-obesity-have-attempted-to-lose-weight-in-the-past-year-but-most-have-been-unsuccessful/](https://easo.org/european-survey-finds-over-three-quarters-of-adults-with-obesity-have-attempted-to-lose-weight-in-the-past-year-but-most-have-been-unsuccessful/) (2022). 4. 4.Creutzfeldt, W. The incretin concept today. Diabetologia 16, 75–85 (1979). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF01225454&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32119&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979GJ69400002&link_type=ISI) 5. 5.Kreymann, B., Williams, G., Ghatei, M. A. & Bloom, S. R. Glucagon-like peptide-1 7-36: a physiological incretin in man. Lancet 2, 1300–1304 (1987). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(87)91194-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2890903&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987L167000005&link_type=ISI) 6. 6.Nauck, M. A. et al. Normalization of fasting glycaemia by intravenous GLP-1 ([7-36 amide] or [7-37]) in type 2 diabetic patients. Diabet. Med. 15, 937–945 (1998). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/(SICI)1096-9136(1998110)15:11<937::AID-DIA701>3.0.CO;2-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9827848&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076834800007&link_type=ISI) 7. 7.Nauck, M. A. et al. Effects of glucagon-like peptide 1 on counterregulatory hormone responses, cognitive functions, and insulin secretion during hyperinsulinemic, stepped hypoglycemic clamp experiments in healthy volunteers. J. Clin. Endocrinol. Metab. 87, 1239–1246 (2002). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.87.3.1239&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11889194&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000174620600048&link_type=ISI) 8. 8.Amylin and Lilly Announce FDA Approval of BYETTA(TM) (Exenatide) Injection. Eli Lilly and Company [https://investor.lilly.com/news-releases/news-release-details/amylin-and-lilly-announce-fda-approval-byettatm-exenatide](https://investor.lilly.com/news-releases/news-release-details/amylin-and-lilly-announce-fda-approval-byettatm-exenatide). 9. 9.Secher, A. et al. The arcuate nucleus mediates GLP-1 receptor agonist liraglutide-dependent weight loss. J. Clin. Invest. 124, 4473–4488 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI75276&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25202980&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 10. 10.Ryan, D. H. et al. Long-term weight loss effects of semaglutide in obesity without diabetes in the SELECT trial. Nat. Med. (2024) doi:10.1038/s41591-024-02996-7. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-024-02996-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=38740993&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 11. 11.Wilding, J. P. H. et al. Once-Weekly Semaglutide in Adults with Overweight or Obesity. N. Engl. J. Med. 384, 989–1002 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2032183&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 12. 12.Rubino, D. M. et al. Effect of Weekly Subcutaneous Semaglutide vs Daily Liraglutide on Body Weight in Adults With Overweight or Obesity Without Diabetes: The STEP 8 Randomized Clinical Trial. JAMA 327, 138–150 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2021.23619&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35015037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 13. 13.Marso, S. P. et al. Liraglutide and Cardiovascular Outcomes in Type 2 Diabetes. N. Engl. J. Med. 375, 311–322 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1603827&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27295427&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 14. 14.Marso, S. P. et al. Semaglutide and Cardiovascular Outcomes in Patients with Type 2 Diabetes. N. Engl. J. Med. 375, 1834–1844 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1607141&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27633186&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 15. 15.Perkovic, V. et al. Effects of Semaglutide on Chronic Kidney Disease in Patients with Type 2 Diabetes. N. Engl. J. Med. 391, 109–121 (2024). 16. 16.Sattar, N. et al. Cardiovascular, mortality, and kidney outcomes with GLP-1 receptor agonists in patients with type 2 diabetes: a systematic review and meta-analysis of randomised trials. Lancet Diabetes Endocrinol 9, 653–662 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2213-8587(21)00203-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34425083&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 17. 17.Watanabe, J. H., Kwon, J., Nan, B. & Reikes, A. Trends in glucagon-like peptide 1 receptor agonist use, 2014 to 2022. J. Am. Pharm. Assoc. 64, 133–138 (2024). 18. 18.Melson, E., Ashraf, U., Papamargaritis, D. & Davies, M. J. What is the pipeline for future medications for obesity? Int. J. Obes. (2024) doi:10.1038/s41366-024-01473-y. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41366-024-01473-y&link_type=DOI) 19. 19.Morgan, J. P. The increase in appetite for obesity drugs. [https://www.jpmorgan.com/insights/global-research/current-events/obesity-drugs](https://www.jpmorgan.com/insights/global-research/current-events/obesity-drugs) (2023). 20. 20.Clapp, B. et al. American Society for Metabolic and Bariatric Surgery 2020 estimate of metabolic and bariatric procedures performed in the United States. Surg. Obes. Relat. Dis. 18, 1134–1140 (2022). 21. 21.After 30 Years — New Guidelines For Weight-Loss Surgery. American Society for Metabolic and Bariatric Surgery [https://asmbs.org/news\_releases/after-30-years-new-guidelines-for-weight-loss-surgery/](https://asmbs.org/news_releases/after-30-years-new-guidelines-for-weight-loss-surgery/) (2022). 22. 22.Gloy, V. L. et al. Bariatric surgery versus non-surgical treatment for obesity: a systematic review and meta-analysis of randomised controlled trials. BMJ 347, f5934 (2013). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNDcvb2N0MjJfMS9mNTkzNCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzEyLzIwMjQuMDkuMTEuMjQzMTM0NTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 23. 23.Salminen, P. et al. Effect of Laparoscopic Sleeve Gastrectomy vs Roux-en-Y Gastric Bypass on Weight Loss, Comorbidities, and Reflux at 10 Years in Adult Patients With Obesity: The SLEEVEPASS Randomized Clinical Trial. JAMA Surg. 157, 656–666 (2022). 24. 24.Lim, R., Beekley, A., Johnson, D. C. & Davis, K. A. Early and late complications of bariatric operation. Trauma Surg Acute Care Open 3, e000219 (2018). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToidHNhY28iO3M6NToicmVzaWQiO3M6MTE6IjMvMS9lMDAwMjE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTIvMjAyNC4wOS4xMS4yNDMxMzQ1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 25. 25.Ma, I. T. & Madura, J. A., 2nd. Gastrointestinal Complications After Bariatric Surgery. Gastroenterol. Hepatol. 11, 526–535 (2015). 26. 26.Azagury, D. et al. Heterogeneity of weight loss after gastric bypass, sleeve gastrectomy, and adjustable gastric banding. Surgery 165, 565–570 (2019). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 27. 27.Cadena-Obando, D. et al. Are there really any predictive factors for a successful weight loss after bariatric surgery? BMC Endocr. Disord. 20, 20 (2020). 28. 28.Çalık Başaran, N., Dotan, I. & Dicker, D. Post metabolic bariatric surgery weight regain: the importance of GLP-1 levels. Int. J. Obes. (2024) doi:10.1038/s41366-024-01461-2. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41366-024-01461-2&link_type=DOI) 29. 29.Lauti, M., Kularatna, M., Hill, A. G. & MacCormick, A. D. Weight Regain Following Sleeve Gastrectomy-a Systematic Review. Obes. Surg. 26, 1326–1334 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11695-016-2152-x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27048439&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 30. 30.Young, K. G. et al. Treatment effect heterogeneity following type 2 diabetes treatment with GLP1-receptor agonists and SGLT2-inhibitors: a systematic review. Commun. Med. 3, 131 (2023). 31. 31.Jones, A. G. et al. Markers of β-Cell Failure Predict Poor Glycemic Response to GLP-1 Receptor Agonist Therapy in Type 2 Diabetes. Diabetes Care 39, 250–257 (2016). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGlhY2FyZSI7czo1OiJyZXNpZCI7czo4OiIzOS8yLzI1MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzEyLzIwMjQuMDkuMTEuMjQzMTM0NTguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 32. 32.Dawed, A. Y. et al. Pharmacogenomics of GLP-1 receptor agonists: a genome-wide analysis of observational data and large randomised controlled trials. Lancet Diabetes Endocrinol 11, 33–41 (2023). 33. 33.Cardoso, P. et al. Phenotype-based targeted treatment of SGLT2 inhibitors and GLP-1 receptor agonists in type 2 diabetes. Diabetologia 67, 822–836 (2024). 34. 34.Brown, E., Wilding, J. P. H., Barber, T. M., Alam, U. & Cuthbertson, D. J. Weight loss variability with SGLT2 inhibitors and GLP-1 receptor agonists in type 2 diabetes mellitus and obesity: Mechanistic possibilities. Obes. Rev. 20, 816–828 (2019). 35. 35.Umapathysivam, M. M. et al. Type 2 Diabetes risk alleles in Peptidyl-glycine Alpha-amidating Monooxygenase influence GLP-1 levels and response to GLP-1 Receptor Agonists. medRxiv (2023) doi:10.1101/2023.04.07.23288197. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMy4wNC4wNy4yMzI4ODE5N3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTIvMjAyNC4wOS4xMS4yNDMxMzQ1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 36. 36.Jensterle, M., Rizzo, M., Haluzík, M. & Janež, A. Efficacy of GLP-1 RA Approved for Weight Management in Patients With or Without Diabetes: A Narrative Review. Adv. Ther. 39, 2452–2467 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12325-022-02153-x&link_type=DOI) 37. 37.Nathan, B. M. et al. Predictors of weight-loss response with glucagon-like peptide-1 receptor agonist treatment among adolescents with severe obesity. Clin. Obes. 6, 73–78 (2016). 38. 38.Vosoughi, K. et al. Association of Glucagon-like Peptide 1 Analogs and Agonists Administered for Obesity with Weight Loss and Adverse Events: A Systematic Review and Network Meta-analysis. EClinicalMedicine 42, 101213 (2021). 39. 39.Segal, J. B. et al. Assessing Heterogeneity of Treatment Effect in Real-World Data. Ann. Intern. Med. 176, 536–544 (2023). 40. 40.Jensterle, M., Pirš, B., Goričar, K., Dolžan, V. & Janež, A. Genetic variability in GLP-1 receptor is associated with inter-individual differences in weight lowering potential of liraglutide in obese women with PCOS: a pilot study. Eur. J. Clin. Pharmacol. 71, 817–824 (2015). 41. 41.van der Meer, R. et al. Genetic variants associated with weight loss and metabolic outcomes after bariatric surgery: A systematic review. Obes. Rev. 24, e13626 (2023). 42. 42.Aasbrenn, M. et al. Genetic Determinants of Weight Loss After Bariatric Surgery. Obes. Surg. 29, 2554–2561 (2019). 43. 43.NCD Risk Factor Collaboration (NCD-RisC). Worldwide trends in underweight and obesity from 1990 to 2022: a pooled analysis of 3663 population-representative studies with 222 million children, adolescents, and adults. Lancet 403, 1027–1050 (2024). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(23)02750-2&link_type=DOI) 44. 44.White, G. E. et al. Real-world weight-loss effectiveness of glucagon-like peptide-1 agonists among patients with type 2 diabetes: A retrospective cohort study. Obesity 31, 537–544 (2023). 45. 45.Inoue, K. et al. Long-term impact of liraglutide, a glucagon-like peptide-1 (GLP-1) analogue, on body weight and glycemic control in Japanese type 2 diabetes: an observational study. Diabetol. Metab. Syndr. 6, 95 (2014). 46. 46.Venkatesh, S. S. et al. The genetic architecture of changes in adiposity during adulthood. medRxiv (2023) doi:10.1101/2023.01.09.23284364. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMy4wMS4wOS4yMzI4NDM2NHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTIvMjAyNC4wOS4xMS4yNDMxMzQ1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 47. 47.Xu, J. et al. Exploring the clinical and genetic associations of adult weight trajectories using electronic health records in a racially diverse biobank: a phenome-wide and polygenic risk study. Lancet Digit Health 4, e604–e614 (2022). 48. 48.Karczewski, K. J. et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 49. 49.Lagou, V. et al. GWAS of random glucose in 476,326 individuals provide insights into diabetes pathophysiology, complications and treatment stratification. Nat. Genet. 55, 1448– 1461 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-023-01462-3&link_type=DOI) 50. 50.Nead, K. T. et al. Contribution of common non-synonymous variants in PCSK1 to body mass index variation and risk of obesity: a systematic review and meta-analysis with evidence from up to 331 175 individuals. Hum. Mol. Genet. 24, 3582–3594 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddv097&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25784503&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 51. 51.Sadler, M. C. et al. Leveraging large-scale biobank EHRs to enhance pharmacogenetics of cardiometabolic disease medications. medRxiv (2024) doi:10.1101/2024.04.06.24305415. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyNC4wNC4wNi4yNDMwNTQxNXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTIvMjAyNC4wOS4xMS4yNDMxMzQ1OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 52. 52.Weiss, T. et al. Real-world weight change, adherence, and discontinuation among patients with type 2 diabetes initiating glucagon-like peptide-1 receptor agonists in the UK. BMJ Open Diabetes Res Care 10, (2022). 53. 53.Ghusn, W. et al. Weight Loss Outcomes Associated With Semaglutide Treatment for Patients With Overweight or Obesity. JAMA Netw Open 5, e2231982 (2022). 54. 54.Xu, T. et al. Timing of Maximal Weight Reduction Following Bariatric Surgery: A Study in Chinese Patients. Front. Endocrinol. 11, 615 (2020). 55. 55.Monpellier, V. M., Janssen, I. M. C., Antoniou, E. E. & Jansen, A. T. M. Weight Change After Roux-en Y Gastric Bypass, Physical Activity and Eating Style: Is There a Relationship? Obes. Surg. 29, 526–533 (2019). 56. 56.Knop, F. K. et al. Oral semaglutide 50 mg taken once per day in adults with overweight or obesity (OASIS 1): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet 402, 705–719 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(23)01185-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37385278&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 57. 57.Hall, K. D. & Kahan, S. Maintenance of Lost Weight and Long-Term Management of Obesity. Med. Clin. North Am. 102, 183–197 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mcna.2017.08.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29156185&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 58. 58.Mok, J. et al. Safety and Efficacy of Liraglutide, 3.0 mg, Once Daily vs Placebo in Patients With Poor Weight Loss Following Metabolic Surgery: The BARI-OPTIMISE Randomized Clinical Trial. JAMA Surg. 158, 1003–1011 (2023). 59. 59.Borrell, L. N. et al. Race and Genetic Ancestry in Medicine - A Time for Reckoning with Racism. N. Engl. J. Med. 384, 474–480 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMms2029562&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33406325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 60. 60.Kemper, K. E. et al. Genetic influence on within-person longitudinal change in anthropometric traits in the UK Biobank. Nat. Commun. 15, 3776 (2024). 61. 61.Weissbrod, O. et al. Leveraging fine-mapping and multipopulation training data to improve cross-population polygenic risk scores. Nat. Genet. 54, 450–458 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01036-9&link_type=DOI) 62. 62.Chang, C. C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13742-015-0047-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25722852&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) 63. 63.Thompson, D. J., et al. UK Biobank release and systematic evaluation of optimised polygenic risk scores for 53 diseases and quantitative traits. bioRxiv (2022) doi:10.1101/2022.06.16.22276246. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2022.06.16.22276246&link_type=DOI) 64. 64.Schwarzer, G., Carpenter, J. R. & Rücker, G. Meta-Analysis with R. (Springer International Publishing, 2015). 65. 65.Cochran, W. G. The Combination of Estimates from Different Experiments. Biometrics 10, 101–129 (1954). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/3001666&link_type=DOI) 66. 66.Higgins, J. P. T. & Thompson, S. G. Quantifying heterogeneity in a meta-analysis. Stat. Med. 21, 1539–1558 (2002). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.1186&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12111919&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F12%2F2024.09.11.24313458.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176016900005&link_type=ISI) 67. 67.Bulus, M. A practical guide to statistical power and sample size calculations in R. [https://cran.r-project.org/web/packages/pwrss/vignettes/examples.html](https://cran.r-project.org/web/packages/pwrss/vignettes/examples.html) (2023). 68. 68.Kanai, M. et al. Insights from complex trait fine-mapping across diverse populations. bioRxiv (2021) doi:10.1101/2021.09.03.21262975. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2021.09.03.21262975&link_type=DOI) [1]: /embed/graphic-6.gif [2]: /embed/graphic-7.gif [3]: /embed/graphic-8.gif [4]: /embed/graphic-9.gif