Methods for predicting weight gain in patients
A machine learning approach using ML-GRS integrates genetic and non-genetic data to predict obesity phenotypes and tailor treatments, addressing inefficiencies in identifying responders and non-responders to anti-obesity medications, enhancing treatment efficacy and reducing adverse effects.
Patent Information
- Authority / Receiving Office
- WO · WO
- Patent Type
- Applications
- Current Assignee / Owner
- Filing Date
- 2025-09-19
- Publication Date
- 2026-03-26
AI Technical Summary
Existing methods struggle to accurately identify individuals who will respond to weight-loss interventions and are susceptible to adverse effects of anti-obesity medications, leading to inefficiencies in treatment and potential adverse effects.
A computer-implemented method using a machine learning assisted gene risk score (ML-GRS) integrates genetic and non-genetic data to predict obesity phenotypes, such as abnormal satiation, and stratify subjects into responders and non-responders to specific pharmacological interventions, including or excluding GLP-1 agonists based on genetic and metabolic profiles.
The method achieves sensitivity and specificity of at least 65% in predicting obesity phenotypes and tailoring treatments, reducing adverse effects and improving treatment efficacy.
Smart Images

Figure US2025047225_26032026_PF_FP_ABST
Abstract
Description
IN THE UNITED STATES PATENT & TRADEMARK RECEIVING OFFICEINTERNATIONAL PCT PATENT APPLICATIONMETHODS FOR PREDICTING WEIGHT GAIN IN PATIENTSCROSS-REFERENCE TO RELATED APPLICATIONS
[0001] This application claims the benefit of priority to U.S. Provisional Application No. 63 / 696,765 filed September 19, 2024, U.S. Provisional Application No. 63 / 696,716 filed September 19, 2024, U.S. Provisional Application No. 63 / 696,761 filed September 19, 2024, U.S. Provisional Application No. 63 / 696,767 filed September 19, 2024, U.S. Provisional Application No. 63 / 715,240 filed November 1 2024, U.S. Provisional Application No. 63 / 744,706 filed January 13, 2025, U.S. Provisional Application No. 63 / 798,977 filed May 2, 2025, U.S. Provisional Application No. 63 / 798,916 filed May 2, 2025, U.S. Provisional Application No. 63 / 798,788 filed May 2, 2025, U.S. Provisional Application No. 63 / 799,026 filed May 2, 2025 and U.S. Provisional Application No. 63 / 799,127 filed May 2, 2025, each of which is hereby incorporated by reference in its entirety for all purposes.STATEMENT REGARDING FEDERAL FUNDING
[0002] This invention was made with government support under DK067071 and DK114460 awarded by the National Institutes of Health. The government has certain rights in the invention.FIELD
[0003] The present disclosure relates to methods for ascertaining a subject’s obesity phenotype using a multi-omics approach or assessing individuals and populations for genetic variants of interest that relate to a particular phenotypic trait (e.g., obesity) and using polygenic risk scores and machine learning to determine the presence of said phenotypic trait based on said genetic variants of interest. The present disclosure also provides methods for ascertaining a subject’s obesity phenotype using a multi-omics approach as well as methods for stratifying subjects with obesity as responders and non-responders to specific pharmacological interventions for treating obesity.BACKGROUND
[0004] Human genetics drive our propensity to have certain characteristics. These characteristics are often complex and result from differences in several parts of our genome.Now that technologies are available to identify genetic differences between individuals (variants), the challenge becomes identifying which variants are relevant to a particular characteristic, such as susceptibility to a disease. Being able to identify relevant variants and how they contribute to a characteristic enables prediction of, among other things, disease risk. Methods that do this identification and predict the likelihood of a characteristic are Polygenic Risk Scores (PRSs).
[0005] Obesity is a chronic and complex disease affecting almost 600 million adults worldwide (see Collaborators GO, N. Engl. J. Med., 377: 13-27 (2017). Obesity complexity is related to a multi-factorial imbalance between energy intake and energy expenditure, leading to excess energy storage (Hill, Wyatt, and Peters 2012). Most interventions aim to reduce caloric intake; however, there is considerable heterogeneity in all responses. Furthermore, several processes that favor weight gain affect this energy balance equilibrium after weight loss (Aronne et al. 2021; Christoff ersen et al. 2022). Consequently, understanding how to control this energy balance in humans is essential for developing successful weight loss interventions.
[0006] The stages of appetite have been identified as hunger, satiation, and postprandial satiety (Blundell and Halford 1994). Satiation is the process that controls the meal size and leads to meal termination (Cannon and Washburn 1993). Satiation is the result of volume-dependent signals from the gastrointestinal tract, as well as enteroendocrine hormones released by the interaction of the gut wall with nutrients (Cifuentes and Acosta 2022). It has been defined as a sense of fullness after eating or even symptoms such as nausea and bloating, and it has been quantified by the number of calories ingested to reach fullness (Camilleri 2015; Halawi et al. 2017). Further, abnormal satiation has been identified as one of the identified phenotypes of obesity (Acosta, Camilleri, Shin, Vazquez -Roque, et al. 2015) and patients with obesity and abnormal satiation are more prone to gain weight compared to patients with normal satiation. Finally, these patients' responses to all obesity therapies may differ. For example, patients with abnormal satiation were able to lose more weight with phentermine-topiramate than patients with normal satiation (Acosta, Camilleri, Shin, Vazquez -Roque, et al. 2015).
[0007] Moreover, anti-obesity medications can cause adverse effects in patients when administered alone or in combination. For example, while glucagon-like peptide-1 receptor agonists (GLP-lRAs) are effective therapies for obesity and type 2 diabetes mellitus, gastrointestinal adverse effects (AEs) are common in patients administered GLP-lRAs alone or in combination. These AEs, in turn, limit AOM (e.g., GLP-1RA) tolerability.
[0008] As a result, it has become essential to generate means for swiftly identifying individuals who would respond to weight-loss interventions such as anti-obesity medications, yet alsodetermine who among such individuals may be susceptible to the adverse effects of anti-obesity medications. The present disclosure aims to solve the aforementioned needs.SUMMARY
[0009] In one aspect, provided herein is a computer-implemented method for generating a machine learning assisted gene risk score (ML-GRS) for predicting an obesity phenotype of interest in a subject suffering from obesity, the method comprising: (a) receiving in a computer system, a genetic dataset comprising a plurality of single nucleotide polymorphisms (SNPs) located in, around or near a set of genes obtained from samples obtained from each subject from a population of subjects, wherein each subject in the population of subjects is obese, wherein each gene in the set of genes is known or suspected to play a role in obesity; (b) generating by the computer system, a contribution factor for each SNP located in, around or near a gene from the set of genes from the genetic dataset for a subject from the population of subjects, wherein the contribution factor is a numerical value, vector, matrix, or function of various biological factors that represents the contribution or effect that each SNP has on a specific trait known to be associated with the obesity phenotype of interest; (c) calculating by the computer system, a gene risk score (GRS) for each gene in the set of genes by integrating the contribution factors generated in step (b) for each SNP for each gene from the set of genes and SNP genotype information located in, around, or near each gene in the set of genes, wherein the GRS for each gene represents the cumulative effect of the variation in SNPs located in, around, or near the gene in that particular subject has on the function or regulation of the gene so as to increase or decrease the risk of the subject possessing a specific trait known to be associated with the obesity phenotype of interest; (d) iteratively performing by the computer system, steps (a)-(c) for each additional subject from the population of subjects to generate a set of GRSs for each gene in the set of genes from the population of subjects, wherein each gene has a GRS for each sex; (e) integrating by the computer system, the GRSs from step (d) with non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects using a machine learning model to calculate a machine learning assisted gene risk score (ML-GRS), wherein the ML-GRS predicts if any one subject possesses the specific trait known to be associated with the obesity phenotype of interest by combining gene risk scores and additional biological, psychological, or environmental measurements of the any one subject; (f) training by the computer system, the ML-GRS from step (e) using a classifier model to predict subjects from a training set of samples from a training population of obese subjects as being positive for the specific trait known to be associated with the obesityphenotype of interest, thereby generating a trained ML-GRS; and (g) determining by the computer system, if a test subject is positive for the specific trait known to be associated with the obesity phenotype of interest by: (i) calculating by the computer system the ML-GRS for the test subject using steps (a)-(e) on a sample obtained from the test subject; and (ii) comparing by the computer system, the ML-GRS for the test subject to the trained ML-GRS, wherein the test subject is positive for the specific trait known to be associated with the obesity phenotype ofinterest if the ML-GRS for the subjectis above the sex-specific 75thpercentile for the specific trait known to be associated with the obesity phenotype of interest or negative for the specific trait known to be associated with the obesity phenotype of interest if the ML-GRS for the subject is below the sex-specific 75thpercentile for the specific trait known to be associated with the obesity phenotype ofinterest. In some cases, the SNP genotype information comprises risk allele dosage, major allele presence, occurrence of a de novo variant, insertion, deletion, or genome rearrangement. In some cases, each gene in the set of genes is known or suspected to play a role in appetite regulation, energy expenditure, lipid metabolism or adipogenesis. In some cases, the subject is obese if the subject possesses a body mass index (BMI) > 30 kg / m2with or without type 2 diabetes. In some cases, the ML-GRS for the obesity phenotype of interest predicts the specific trait known to be associated with the obesity phenotype of interest with a sensitivity and / or specificity is at least 65%. In some cases, the ML-GRS predicts the specific trait known to be associated with the obesity phenotype of interest with an area under the curve (AUC) of at least 0.65, 0.7, 0.75, 0.8, 0.85, 0.9 0.95 or 1. In some cases, the ML-GRS predicts the specific trait known to be associated with the obesity phenotype of interest with positive predictive value of at least 50%. In some cases, the contribution factor is a function or algorithm that utilizes output of a regression analysis or a statistical analysis on the genetic dataset as well as one or more additional parameters, covering regression against one or more traits using one or more types of genetic datasets selected from the group consisting of a genome-wide data microarray chip (GWAS), targeted sequencing (exome or targeted genetic panel), variants detected by qPCR through targeted amplification, whole genome sequencing (WGS), targeted or untargeted sequencing of genomic re-arrangement, deletions, duplications, repeat extensions, and complex genotyping of HLA and CYP genes. In some cases, the one or more additional parameters relate to biological factors regarding genetics, epigenetics, and gene regulation for each SNP. In some cases, the one or more additional parameters are selected from the group consisting of proximity of the SNP to the specific gene, known or theoretical mechanistic role of the SNP metadata regarding the SNP in relationship to the gene, and any combination thereof. In some cases, the known or theoretical mechanistic role of the SNP isselected from the group consisting of synonymous mutation, nonsynonymous mutation, frameshift mutation and nonsense mutation. In some cases, the known or theoretical mechanistic role of the SNP affects the three-dimensional structure of the translated protein binding sites, protein-protein interaction sites and modification sites, wherein the modification sites are selected from the group consisting of phosphorylation, glycosylation and proteolytic site. In some cases, the mechanistic role comprises whether or not the SNP is within an intron, exon, enhancer region or regulatory region of the gene, a cis-regulatory element, a promoter region, a non-coding exonic region, a coding exonic region, intronic region, splice site, transcription factor binding site, epigenetic modification site, at a remote genomic location involved in a three-dimensional chromatin contact with the gene, or a cell-type-specific topological-associated domain that contains the gene. In some cases, the cis-regulatory element is an enhancer or insulator. In some cases, the epigenetic modification site comprises DNA methylation modifications, histone modifications, and / or chromatin accessibility. In some cases, the non-genetic and / or multi-omic data is selected from the group consisting of metabolomic data, proteomic data, peptidomic data, epigenetic data, microbiome data, results from one or more questionnaires and any combination thereof. In some cases, the epigenetic data comprises a presence or absence of DNA modifications selected from the group consisting of DNA methylation modifications, histone modifications and chromatin accessibility. In some cases, the microbiome data is from any site on the subj ect or the subj ecf s environment. In some cases, the machine learning model in step (e) comprises a forward feature selection, backward feature selection, or random feature sampling algorithm and a random forest, GBM, or linear predictor that iteratively selects a feature with a desired training accuracy using the random forest predictor in multiple rounds until no further training accuracy improvement is detected. In some cases, the GRSs of step (e) are normalized. In some cases, each GSR provides an indication of a role or effect that the SNPs located around or near a specific gene may play on the specific trait known to be associated with the obesity phenotype of interest. In some cases, the SNP is considered to be located around or near a specific gene if it is located within at least 500,000 kilobases (kb) of the specific gene. In some cases, the method further comprises performing a genetic analysis on a sample obtained from the subject suffering from or suspected of suffering from obesity prior to step (a). In some cases, the genetic analysis comprises obtaining sequence reads from the whole or portions of the whole genome of the subject. In some cases, the genetic analysis comprises performing a genotyping method selected from the group consisting of restriction fragment length polymorphism identification (RFLPI), random amplified polymorphic detection (RAPD), amplified fragment lengthpolymorphism detection (AFLPD), polymerase chain reaction (PCR), DNA sequencing, RNA sequencing, allele specific oligonucleotide (ASO) probes, and hybridization to microarrays or beads. In some cases, the method further comprises performing by the computer system, a regression analysis and / or statistical analysis on the genetic dataset from step (a) prior to step (b). In some cases, the regression analysis is a ridge regression or least absolute shrinkage and selection operator (LASSO) regression. In some cases, beta values obtained from the regression analyses performed on the genetic datasets are used to generate the contribution factor in step (b). In some cases, the statistical analysis is a genome-wide association study (GWAS). In some cases, p-values and / or regression weights obtained from the statistical analyses performed on the genetic datasets are used to generate the contribution factors in step (b). In some cases, the obesity phenotype of interest is selected from the group consisting of hungry brain (abnormal satiation), hungry gut (abnormal satiety), slow bum (slow metabolism) or hedonic / emotional eating). In some cases, the obesity phenotype of interest is a hungry brain (abnormal satiation) or hungry gut (abnormal satiety). In some cases, the specific trait in step (b) is total kcal consumed to satiation (CTS) at an ad libitum meal. In some cases, following step (e), the subject is CTS positive if the subject has a ML-GRS for CTS above the sex-specific seventyfifth percentile or the 4thquartile, while the subject is CTS negative if the subject has a ML- GRS for CTS below the sex-specific seventy-fifth percentile. In some cases, the plurality of genes in step (a) are two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MB0AT4 (GOAT), ADRA2A, GNB3, GPBAR1, LEP, LEPR, SH2B1, SIM1, NC0A1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LING01, FANCL, CADM2, PTBP2, NUDT3, 0LFM4, H0XB5, NPY1R, NPY2R, NPY4R, NPY5R, NR1H4, SLC6A4, NTS, UCP3, ADIPOQ, CCK, CCKR, DPP4, FGFR4, GHRL, MC4R, MC3R, POMC, AGRI5, GIP, GIPR, CNR1, I- A AH, GCG, CEL A2A, PPARG, TFAP2B, APOE, TNFRSF 11 A or DYRK1B. In some cases, the obesity phenotype of interest is a hungry gut (abnormal satiety). In some cases, the specific trait is increased or accelerated baseline gastric emptying for the subject as compared to a control subject, wherein the control subject is not obese. In some cases, the set of genes in step (a) consists of genes known to play role in conferring the hungry gut phenotype. In some cases, the obesity phenotype of interest is a slow burn. In some cases, the specific trait is decreased or low resting energy expenditure (REE) for the subject as compared to a control, wherein the control subject is not obese. In some cases, the set of genes in step (a) consists of genes known to play role in conferring the slow burn phenotype. In some cases, the obesity phenotype of interest is an emotional eating phenotype. In some cases, the specific trait is a finding or result from a behavioral questionnaireindicative of anxiety or emotional eating for the subject. In some cases, the questionnaire is the Hospital Anxiety and Depression Scale (HADS) questionnaire. In some cases, the plurality of genes in step (a) consists of genes known to a play role in conferring the emotional eating phenotype.
[0010] In another aspect, provided herein is a computer-implemented method for predicting or diagnosing a subject suffering from or suspected of suffering from obesity as possessing abnormal satiation, the method comprising: (a) receiving in a computer system, a genetic dataset comprising a plurality of single nucleotide polymorphisms (SNPs) located in, around or near a set of genes obtained from samples obtained from each subject from a population of subjects, wherein each subject in the population of subjects is obese, wherein the set of genes comprises GLP1R, UCP2, FTO, TCF7L2, MBOAT4 (GOAT), ADRA2A, GNB3, GPBAR1, LEP, NPY1R, NPY2R, NPY4R, NPY5R, NR1H4, SLC6A4, NTS, UCP3, ADIPOQ, CCK,CCKR,DPP4, FGFR4, GHRL, MC4R, MC3R, POMC, AGRP, G1P, GIPR, CNR1, FA AH, GCG, CI ’L A2A, PPARG, TFAP2B, APOE, LEPR, SH2B1, SIM1, NCOA1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO1, FANCL, CADM2, PTBP2, NUDT3, OLFM4, HOXB5, TNFRSF11 A and DYRK1B; (b) generating by the computer system, a contribution factor for each SNP located in or near a gene from the set of genes from the genetic dataset for a subject from the population of subjects, wherein the contribution factor is a numerical value, vector, matrix, or function of various biological factors that represents the contribution or effect that each SNP has on total kcal consumed to satiation (CTS) at an ad libitum meal; (c) calculating by the computer system, a gene risk score (GRS) for each gene in the set of genes by integrating the contribution factors generated in step (b) for each SNP for each gene from the set of genes and SNP genotype information located in, around or near each gene in the set of genes, wherein the GRS for each gene represents the cumulative effect of the variation in SNPs located in, around, or near the gene in that particular subject has on the function or regulation of the gene so as to increase or decrease the risk of being above a sex-specific threshold total kcal consumed to satiation (CTS) at an ad libitum meal; (d) iteratively performing by the computer system, steps (a)-(c) for each additional subject from the population of subjects to generate a set of GRSs for each gene in the set of genes, wherein each gene has a GRS for each sex; and(e) integrating by the computer system, the GRSs from step (d) with non-genetic and / or multi- omic data obtained from the samples obtained from each subject from the population of subjects using a machine learning model to calculate a machine learning assisted gene risk score (ML-GRS), wherein the ML-GRS predicts the CTS by combining gene risk scores andother biological, psychological, or environmental measurements of the subject; (f) training by the computer system, the ML-GRS from step (e) using a classifier model to predict subjects from a training set of samples from a training population of obese subjects as being positive for being above a sex-specific threshold total kcal consumed to satiation (CTS) at an ad libitum meal, thereby generating a trained ML-GRS; and (g) determining by the computer system, if a test subject is positive for the specific trait known to be associated with the obesity phenotype of interest by: (i) calculating by the computer system the ML-GRS for the test subject using steps (a)-(e) on a sample obtained from the test subject; and (ii) comparing by the computer system, the ML-GRS for the test subject to the trained ML-GRS, wherein the test subject is positive for CTS if the ML-GRS for the subject is above the sex-specific 75thpercentile for CTS or negative for CTS if the ML-GRS for the subject is below the sex-specific 75thpercentile for CTS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample.
[0011] In another aspect, provided herein is a method for treating obesity in a subject in need thereof, the method comprising: (a) determining if the subject possesses abnormal satiation using the method of embodiment 46 on a genetic dataset obtained from a sample obtained the subject to determine a ML-GRS for abnormal satiation; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the ML-GRS determined for the subject indicates that the subject is positive for CTS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the ML-GRS determined for the subject indicates that the subject is negative for CTS.
[0012] In another aspect, provided herein is a method for treating obesity in a subject in need thereof, the method comprising: (a) determining the presence, absence or level of a plurality of gastrointestinal (GI) peptides, a plurality of metabolites, and / or a plurality of genetic variants in a sample obtained from a subject suffering from obesity as well as satiety, satiation, resting energy expenditure and results on a behavioral questionnaire for the subject, thereby generating an obesity analyte signature for the sample; (b) populating a predictive machine learning model with the obesity analyte signature for the subject; and (c) utilizing the predictive machine learning model to predict an obesity phenotype of the subject suffering from obesity, wherein the obesity phenotype is selected from the group consisting of hungry brain (abnormal satiation), hungry gut (abnormal satiety), slow bum (slow metabolism) or hedonic / emotional eating); and (d) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the predictive machine learning model to predicts that the subject possesses an abnormal satiation phenotype or administering a pharmacotherapy that may include a GLP-1agonist to the subject if the predictive machine learning model to predicts that the subject does not possess an abnormal satiation phenotype.
[0013] In some cases, the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate. In some cases, the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatidey maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006, BMF-650 and any combination thereof.
[0014] In yet another aspect, provided herein is a method of treating obesity in a subject in need thereof, the method comprising: (a) detecting an obesity analyte signature in a sample obtained from a subject; wherein the obesity analyte signature is indicative of an obesity phenotype of the subject; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the obesity analyte signature detected in the sample obtained from the subject indicates that the subject possesses an abnormal satiation phenotype or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the obesity analyte signature detected in the sample obtained from the subject indicates that the subject does not possess an abnormal satiation phenotype. In some cases, the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS- 4571, ASC-30, TTP273, MLX-7006, BMF-650 and any combination thereof.
[0015] In another aspect, provided herein is a method of identifying a subject suffering from obesity as a non-responder to treatment with a GLP-1 agonist, the method comprising detecting an obesity phenotype of the subject; and identifying the subject as a non-responder to GLP-1 agonist treatment if the subject is determined to possess an abnormal satiation obesity phenotype.
[0016] In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the obesity analyte signature comprises the presence of serotonin, glutamine, isocaproic acid, alloisoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic acid, tyrosine, and PYY, an absence of 1-methylhi stine, gamma-amino-n-butyric-acid, phenylalanine, ghrelin, and a HADS questionnaire result that does not indicate an anxiety subscale (HADS-A; e.g., includes a HADS-A questionnaire result). In some cases, the subject consumes above the seventy-fifth percentile amount of calories at an ad libitum meal. In some cases, the obesity analyte signature comprises the presence of 1-methylhi stine, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, and phenylalanine, an absence of serotonin, glutamine, gamma-amino-n-butyric-acid, isocaproic acid, hexanoic acid, tyrosine, ghrelin, PYY, and does not include a HADS questionnaire result that indicates an anxiety subscale. In some cases, the subject has increased or accelerated baseline gastric emptying as compared to a control subject, wherein the control subject is not obese. In some cases, the obesity analyte signature comprises the presence of 1-methylhistine, serotonin, glutamine, gamma-amino-n-butyric-acid, isocaproic acid, allo-isoleucine, alanine, tyrosine, ghrelin, PYY, an absence of hydroxyproline, beta-aminoisobutyric-acid, hexanoic acid, and phenylalanine, and does includes a HADS questionnaire result that indicates an anxiety subscale. In some cases, the subject has decreased or low resting energy expenditure (REE) for the subject as compared to a control, wherein the control subject is not obese. In some cases, the obesity analyte signature comprises the presence of serotonin, an absence of 1-methylhistine, glutamine, gamma-amino-n-butyric-acid, isocaproic acid, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic acid, tyrosine, phenylalanine, ghrelin, and PYY, and does includes a HADS questionnaire result that indicates an anxiety subscale.
[0017] In one aspect, provided herein is a computer-implemented method for predicting a subject’s calories to satiation (CTS), the method comprising: (a) receiving in a computer system, genetic data obtained from a sample obtained from a subject as well as anthropometric data from the subject; (b) calculating by the computer system, a gene risk score (GRS) for each respective gene in a set of genes from the genetic data obtained from the sample obtained from the subject, wherein the GRS for each respective gene in the set of genes represents a cumulative effect of single nucleotide polymorphisms (SNPs) located in, around, or near the respective gene has on the function or regulation of the respective gene in the subject so as to increase or decrease the risk of the subject being above a sex-specific threshold of total kcal consumed to satiation (CTS) at an ad libitum meal, wherein the set of genes comprises, consists essentially of or consists of the set of genes listed in Table 5; (c) normalizing by the computer system, the anthropometric data and the GRS for each respective gene from the set of genes to training samples from a CTSGRS classifier to generate normalized GRSs and anthropometricdata for the subject; and (d) inputting by the computer system, the normalized GRSs and anthropometric data for the subject into the CTSGRS classifier to generate a value ranging from 0 to 1, wherein a value of greater than or equal to 0.5 classifies the subject as having a high CTSGRS and a value of less than 0.5 classifies the subject as having a low CTSGRS, wherein a high CTSGRS indicates a probability of being above the sex-specific 75thpercentile for CTS and a low CTSGRS indicates a probability of being below the sex-specific 75th percentile for CTS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the sexspecific threshold of total kcal consumed to satiation (CTS) at an ad libitum meal in (b) is the sex-specific 75thpercentile. In some cases, the sex-specific 75thpercentile for males is about 1370 kCal and for females is about 975 kCal. In some cases, the anthropometric data comprises, consists essentially of or consists of height, weight, age and sex. In some cases, the set of genes comprises, consists essentially of or consists of GLP1R, UCP2, FTO, TCF7L2, ADRA2A, LEPR, SH2B1, SIM1, PCSK1 and TNFRSFUA. In some cases, the CTSGRS classifier was generated by: (a) receiving in a computer system, a genetic dataset comprising a plurality of single nucleotide polymorphisms (SNPs) located in, around or near a set of genes obtained from samples obtained from each subject from a population of subjects, wherein each subject in the population of subjects is obese, wherein the set of genes comprises, consists essentially of or consists of the genes listed in Table 5; (b) generating by the computer system, a contribution factor for each SNP located in, around or near a gene from the set of genes from the genetic dataset for a subject from the population of subjects, wherein the contribution factor is a numerical value, vector, matrix, or function of various biological factors that represents the contribution or effect that each SNP has on total kcal consumed to satiation (CTS) at an ad libitum meal; (c) calculating by the computer system a gene risk score (GRS) for each gene in the set of genes by integrating the contribution factors generated in step (b) for each SNP for each gene from the set of genes and SNP genotype information located in, around or near each gene in the set of genes, wherein the GRS for each gene represents the cumulative effect of the variation in SNPs located in, around, or near the gene in that particular subject has on the function or regulation of the gene so as to increase or decrease the risk of being above a sexspecific threshold total kcal consumed to satiation (CTS) at an ad libitum meal; (d) iterativelyperforming by the computer system, steps (a)-(c) for each additional subject from the population of subjects to generate a set of GRSs for each gene in the set of genes, wherein each gene has a GRS for each sex; and (e) integrating by the computer system, the GRSs from step (d) with non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects using a machine learning model to calculate a machine learning assisted gene risk score (ML-GRS) for CTS (CTSGRS), wherein the CTSGRS predicts the CTS at an ad libitum meal for each subject from the population of subjects by combining gene risk scores and the non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects; (f) training by the computer system, the CTSGRS from step (e) using a classifier model to generate the CTSGRS classifier, wherein the CTSGRS classifier generates a value ranging from 0 to 1 for each subject, wherein a value of greater than or equal to 0.5 classifies the subject as having a high CTSGRS and a value of less than 0.5 classifies the subject as having a low CTSGRS, wherein a high CTSGRS indicates a probability of the subject being above the sex-specific 75thpercentile for CTS and a low CTSGRS indicates a probability of the subject being below the sex-specific 75th percentile for CTS. In some cases, the SNP genotype information comprises risk allele dosage, major allele presence, occurrence of a de novo variant, insertion, deletion, or genome rearrangement. In some cases, each gene in the set of genes is known or suspected to play a role in appetite regulation, energy expenditure, lipid metabolism or adipogenesis. In some cases, the subject is obese if the subject possesses a body mass index (BMI) > 30 kg / m2with or without type 2 diabetes. In some cases, the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with a sensitivity and / or specificity is at least 65%. In some cases, the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with an area under the curve (AUC) of at least 0.65, 0.7, 0.75, 0.8, 0.85, 0.9 0.95 or 1. In some cases, the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with positive predictive value of at least 50%. In some cases, the contribution factor is a function or algorithm that utilizes output of a regression analysis or a statistical analysis on the genetic dataset on CTS. In some cases, the genetic dataset is selected from the group consisting of a genome-wide data microarray chip (GWAS), targeted sequencing (exome or targeted genetic panel), variants detected by qPCR through targeted amplification, whole genome sequencing (WGS), targeted or untargeted sequencing of genomic re-arrangement, deletions, duplications, repeat extensions, and complex genotyping of HLA and CYP genes. In some cases, the non-genetic data is anthropometric data selected from the group consisting of age, weight, height, sex and a combination thereof. In some cases,the non-genetic and / or multi-omic data is selected from the group consisting of metabolomic data, proteomic data, peptidomic data, epigenetic data, microbiome data, results from one or more questionnaires and any combination thereof. In some cases, the epigenetic data comprises a presence or absence of DNA modifications selected from the group consisting of DNA methylation modifications, histone modifications and chromatin accessibility. In some cases, the microbiome data is from any site on the subject or the subject’s environment. In some cases, the machine learning model in step (e) comprises a forward feature selection, backward feature selection, or random feature sampling algorithm and a random forest, GBM, or linear predictor that iteratively selects CTS with a desired training accuracy using the random forest predictor in multiple rounds until no further training accuracy improvement is detected. In some cases, the GRSs of step (e) are normalized. In some cases, each GSR provides an indication of a role or effect that the SNPs located around or near a specific gene may play on the specific trait known to be associated with CTS at an ad libitum meal. In some cases, the SNP is considered to be located around or near a specific gene if it is located within at least 500,000 kilobases (kb) of the specific gene. In some cases, the method further comprises performing a genetic analysis on a sample obtained from the subject suffering from or suspected of suffering from obesity prior to step (a). In some cases, the genetic analysis comprises obtaining sequence reads from the whole or portions of the whole genome of the subject. In some cases, the genetic analysis comprises performing a genotyping method selected from the group consisting of restriction fragment length polymorphism identification (RFLPI), random amplified polymorphic detection (RAPD), amplified fragment length polymorphism detection (AFLPD), polymerase chain reaction (PCR), DNA sequencing, RNA sequencing, allele specific oligonucleotide (ASO) probes, and hybridization to microarrays or beads. In some cases, the method further comprises performing by the computer system, a regression analysis and / or statistical analysis on the genetic dataset from step (a) prior to step (b). In some cases, the regression analysis is a ridge regression or least absolute shrinkage and selection operator (LASSO) regression. In some cases, beta values obtained from the regression analyses performed on the genetic datasets are used to generate the contribution factor in step (b). In some cases, the statistical analysis is a genome-wide association study (GWAS). In some cases, p-values and / or regression weights obtained from the statistical analyses performed on the genetic datasets are used to generate the contribution factors in step (b).
[0018] In another aspect provided herein is a method of predicting weight gain in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein ahigh CTSGRS indicates that the subject is above the sex-specific 75thpercentilefor CTS and is predicted to gain more weight over time as compared to a subject with a low CTSGRS. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0019] In yet another aspect, provided herein is a method for treating weight gain in a subject in need thereof, the method comprising: (a) determining if a subject possesses a high calorie to satiation (CTS) using any of the methods provided herein on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate. In some cases, the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006, BMF-650 and any combination thereof.
[0020] In another aspect provided herein is a method of predicting obesity in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein ahigh CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to develop obesity over time as compared to a subject with a low CTSGRS.
[0021] In yet another aspect, provided herein is a method for treating obesity in a subject in need thereof, the method comprising: (a) determining if a subject possesses a high calorie to satiation (CTS) using any of the methods provided herein on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate. In some cases, the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650.
[0022] In one embodiment, provided herein is a method for treating weight gain in a subject in need thereof that comprises administering a pharmacotherapy that may include a GLP-1 agonist to a subject if the subject possesses a low calories to satiation or administering a pharmacotherapy that does not include a GLP-1 agonist if the subject possesses a high calories to satiation. The method can comprise: (a) determining if a subject possesses a high calorie to satiation (CTS) using any of the methods provided herein on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity.In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate. In some cases, the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650.
[0023] In another aspect, provided herein is a method of predicting response to a bariatric surgical intervention in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a high CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to respond to the bariatric surgical intervention as compared to a subject with a low CTSGRS. In some cases, the response to the bariatric surgical intervention is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4 or 5 years following the bariatric surgical intervention. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0024] In another aspect, provided herein is a method of predicting response to endoscopic sleeve gastroplasty in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a low CTSGRS indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to endoscopic sleeve gastroplasty as compared to a subject with a high CTSGRS. In some cases, the response to endoscopic sleeve gastroplasty is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following endoscopic sleeve gastroplasty, endoscopy or surgical intervention that changes stomach size or shape. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subjecthas a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0025] In another aspect, provided herein is a method of treating diabetes or prediabetes in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a low CTSGRS (i.e., hungry gut positive (HG+ or ML-GRS<0.5)) indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to a GLP-1 receptor agonist as compared to a subject with a high CTSGRS ((i.e., hungry gut negative (HG- or ML-GRS>0.5)); and administering the GLP-1 receptor agonist. In some cases, the response to the GLP-1 receptor agonist is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4 or 5 years following administration of the GLP-1 receptor agonist. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes. In some cases, the administration of the GLP-1 agonist to the subject results in a reduction in the subject’s HbAlc. In some cases, the reduction in the subject’s HbAlc is at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4 or 5 years following administration of the GLP-1 receptor agonist. In some cases, the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650.
[0026] In another aspect, provided herein is a method of selecting a treatment for a subject suffering from diabetes or prediabetes, the method comprising determining the CTSGRS of a subject using any of the methods of provided herein, wherein a low CTSGRS (i.e., hungry gut positive (HG+ or ML-GRS<0.5)) indicates that the subject is a responder to a GLP-1 receptor agonist as compared to a subject with a high CTSGRS ((i.e., hungry gut negative (HG- or ML- GRS>0.5)). In some cases, the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650.
[0027] In another aspect provided herein is a method of predicting obesity in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein ahigh CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to develop obesity over time as compared to a subject with a low CTSGRS.
[0028] In yet another aspect, provided herein is a method for treating obesity in a subject in need thereof, the method comprising: (a) determining if a subject possesses a high calorie to satiation (CTS) using any of the methods provided herein on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate. In some cases, the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650.
[0029] In yet another aspect, provided herein is a method of predicting one or more adverse effects (AEs) in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a high CTSGRS (i.e., hungry gut positive (HG- or ML-GRS>0.5)) indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to experience adverse effects when administered a GLP-1 receptor agonist aloneor in combination with one or more additional therapeutic agents as compared to a subject with a low CTSGRS ((i.e., hungry gut positive (HG+ or ML-GRS<0.5)). In some cases, the one or more AE is selected from the group consisting of diarrhea, constipation, abdominal discomfort, nausea, vomiting, headache and any combination thereof. In some cases, the AE is nausea. In some cases, the AE is nausea and headache. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes. In some cases, the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650. In some cases, the one or more additional therapeutic agents is an antihyperglycemic medications (OAMs). In some cases, the 0AM is insulin. In some cases, the one or more additional therapeutic agents is a sodium-glucose cotransporter-2 (SGLT-2) inhibitor. In some cases, the SGLT-2 inhibitor is selected from the group consisting of empagliflozin, canagliflozin, dapagliflozin, ertugliflozin, bexagliflozin and sotagliflozin. In some cases, the one or more additional therapeutic agents is a dipeptidyl peptidase-4 (DPP-4) inhibitor. In some cases, the DPP-4 inhibitor is selected from the group consisting of sitagliptin, linagliptin, saxagliptin and alogliptin.BRIEF DESCRIPTION OF THE FIGURES
[0030] FIG. 1 depicts the study protocol for in-person physiological testing. After an 8-hour fast, all participants in the study provided a plasma blood sample, had physiologic in-person testing to measure food intake, which included a standardized 320 kcal breakfast, gastric emptying for solids for 4 hours, and an ad libitum meal.
[0031] FIGs 2A-2H depicts the inter-individual variability postprandial satiety and gender differences. FIG. 2A illustrates gastric emptying curve differences between males and females. Half time for gastric emptying for females was 1 hour and 51 min (right vertical dotted line) and for males was 1 hour and 25 minutes (left vertical dotted line). FIG. 2B illustrates visual analogue score (VAS) for hunger in males and females before and after breakfast. FIG. 2C illustrates VAS-fullness in males and females before and after breakfast. FIG. 2D illustrates VAS-satisfaction in males and females before and after breakfast. FIG. 2E illustrates caloriesconsumed by macronutrient composition between females and males. FIG. 2F illustrates total calories to fullness in ad libitum meal in all participants. FIG. 2G illustrates frequency histograms of the density of calories to fullness in males. FIG. 2H illustrates frequency histograms of the density of calories to fullness in females.
[0032] FIG. 3 illustrates determinants postprandial satiety variability. Trait variations for each input variable are derived from separate (non-hierarchical) regression models. Values represent the adjusted proportion of variance explained (R2), and error bars show the 95% CI. R2 values were derived using multivariable linear regression.
[0033] FIGs 4A-4H shows the characterization of abnormal satiation. FIG. 4A shows half gastric emptying time of participants with normal and abnormal satiation. FIGs 4B-4D shows plasma hormone levels (FIG. 4B-Ghrelin; FIG. 4C-PYY; and FIG. 4D-GLP-1) at baseline and 15, 45 and 90 minutes after breakfast (320 kcal). FIGs 4E-4H shows appetite sensations including satisfaction (FIG. 4E), hunger (FIG. 4F), fullness (FIG. 4G), and desire to eat (FIG. 4H) measured by visual analogue scales before and after breakfast. Data is presented as mean and standard error of the mean.
[0034] FIGs 5A-5B show ROC curves of the predictive model for abnormal satiation (i.e., the machine learning (ML)-driven full phenotype test (FPT)) in the training cohort (FIG. 5A) and in a validation cohort (FIG. 5B).
[0035] FIGs 6A-6B shows total body weight loss % by prediction group with liraglutide (FIG. 6A) and placebo (FIG. 6B). True is responder to GLP-1.
[0036] FIG. 7 depicts a diagram of the polygenic risk score (PRS) assessment of the present invention that integrates expert guidance, machine learning, and non-genetic / multi-omic data in order to determine the obesity phenotype of a subject suffering from obesity.
[0037] FIGs 8A-8B show ROC curves of the ML-driven polygenic risk score (PRS) model for appetite in the training cohort (FIG. 8A) and in a validation cohort (FIG. 8B). FIGs. 8C-8D show ROC curves of the PRS model for GLP-1 responder in the training cohort (FIG. 8C) and validation cohort (FIG. 8D)
[0038] FIGs 9A-9B show total body weight loss % in groups predicted as responder vs. nonresponder using the PRS model for appetite for patients treated with liraglutide (3mg; FIG. 9A) or placebo (FIG. 9B). True is responder to GLP-1.
[0039] FIGs 10A-10D shows an experimental design that includes data collection by deep phenotype testing in 717 participants (FIG. 10A), feature extraction for satiation characterization (FIG. 10B), model training and validation for gene risk score development (FIG. 10C) and model extrapolation for prediction of response to anti-obesity medicine (AOM)interventions (FIG. 10D). FIG. 10A illustrates the rigorous testing day starting at 7 am, including resting energy expenditure assessment, blood sampling, radio-labeled breakfast, gastric emptying scans, body composition by DEXA imaging, ad libitum meal test, and behavioral questionnaires. FIG. 10B shows male and (n=179) and female (n=538) satiation distribution by calories to satiation, highlighting factors like genetics, medical records, metabolites, and hormones that contribute to individual satiation responses. FIG. IOC demonstrates the creation and validation of a polygenic risk score (PRS) and machine-learning optimized gene risk score (GRS) for satiation. FIG. 10D depicts the model's effectiveness in predicting response to obesity interventions using outcomes in a 52-week randomized, placebo- controlled trial of phentermine-topiramate ER and a 16-week randomized, placebo-controlled trial of liraglutide by categorizing patients according to CTS by ad libitum meal or by genotype groups.
[0040] FIGs 11A-11F shows a comprehensive analysis of satiation variability and influencing factors that includes inter-individual variability in satiation (FIG. 11 A), variability of satiation explained by input variables (FIG. 11B), sex -related effects on satiation parameters (FIG. 11C), anthropometric influence on satiation (FIG. 11D), body composition influence on satiation (FIG. HE) and questionnaire-based satiation patterns (FIG. HF). FIG. HA illustrates the distribution of satiation responses among 717 participants, highlighting the wide range of caloric intake to satiation. FIG. 11B depicts the proportional variance explained (R2) for satiation variability attributed to individual input variables, obtained from separate non- hierarchical regression models. Error bars represent the 95% confidence interval. R2 values were calculated through multivariable linear regression. FIG. 11C is a violin plot that showcases the distinct impact of sex on satiation parameters, revealing higher caloric intake required to reach fullness in males (p<0.001). FIG. HD shows a univariate linear regression analyses display the weak correlation between calories to satiation and height in 717 participants. FIG. HE shows univariate linear regression analyses display the correlation between calories to satiation and total fat percent in 395 participants that completed a DEXA scan. FIG. HF demonstrates the limited contribution of behavioral questionnaires to explain the variability of objective measurements of satiation as demonstrated by the weak correlation between calories to satiation and scores in the components of the three-factor eating questionnaire: uncontrolled eating, cognitive restrain, and emotional eating completed in 215 participants.
[0041] FIG. 12A shows comparative analysis of high and low satiation groups by sex and illustrates the comprehensive distinction between high and low calories to satiation, furthercategorized as high calories to satiation (n=179) with >977 kcal for females (n=134) and >1374 kcal for males (n=45), and low calories to satiation (n=161) with <650 kcal for females (n=134) and <927 kcal for males (n=44). FIG. 12B shows gastric emptying time by Gender in High and Low Calories to Satiation with significant differences in solid gastric emptying half-time in females (p<0.001) but not in males (p=0.15). FIG. 12C shows fullness by visual analog scale (VAS) after consuming a 320-kcal breakfast - High Calories to Satiation (n=70) vs. Low Calories to Satiation (n=79). FIG. 12D shows hunger levels after a 320-kcal Breakfast - High Calories to Satiation (n=70) vs. Low Calories to Satiation (n=79). FIG. 12E shows profile of Gastrointestinal Hormones Fasting and Postprandial in High and Low Calories to Satiation: Profile of Peptide YY (PYY) Levels in High (n=83) and Low Satiation (n=59) showing not significant difference between groups at any time point. FIG. 12F shows profile of GLP-1 Levels in High (n=123) and Low Calories to Satiation (n=101) showing only significant differences between groups at 90 minutes. FIG. 12G shows profile of Ghrelin Levels in High (n=85) and Low Calories to Satiation (n=58) showing not significant difference between groups at any time point. FIG. 12H shows the distribution of eating behavior scores in the three-factor eating questionnaire (n=215) (cognitive restraint, uncontrolled eating, emotional eating) between High and Low Calories to Satiation.
[0042] FIG. 13A shows a Genome-Wide Association Study (GWAS) for Calories to Satiation. This Manhattan plot displays the results of a GWAS investigating genetic variants associated with calories to satiation in the ad libitum meal. The absence of significant peaks across the genome highlights that no genetic variants reached the threshold for genome-wide significance (p < 5 x 10A-8).
[0043] FIG. 13B shows a volcano plot of association with calories to satiation traits. Each point represents the association of a variability of a SNP with CTS among male (blue) and female (red) subjects. The x-axis shows the magnitude and direction of the association and the y-axis the significance of the association.
[0044] FIG. 13C shows a beeswarm plot of Shapley Additive Explanations (SHAP) values for the CTSGRS model representing the contribution of specific features to model output. Each point for a given feature represents a sample with its horizontal position representing the magnitude and direction that feature effects the probability of high calories to satiation. The color scale indicates the relative value of that feature for that sample as ranked among all other samples. Features are ranked by the mean absolute SHAP value, i.e., magnitude of effect on predictions.
[0045] FIGs 14A-14G illustrates the development and validation of a machine learning assisted gene risk score for high calories to satiation in people with obesity. FIG. 14A showsrisk allele distribution for satiation. The number of risk alleles associated with higher calories to satiation was normally distributed in our population. Each circle represents the mean number of calories to satiation during the ad libitum meal test for each risk allele score. The solid line represents the linear regression for the risk allele score and calories to satiation (R2 = 0.55; p coefficient, 6.78; 95% CI, 5.46 to 8.11). FIG. 14B shows receiver operating characteristics curve for the Machine Learning Assisted Gene Risk Score in the training (n=483), validation (n=50), and independent validation cohort (n=l 10). FIG. 14C shows caloric intake in males and females according to their Machine Learning Assisted Gene Risk Score (CTSGRS) groups. Participants predicted as low CTSGRS had a lower caloric intake than participants predicted as high CTSGRS. FIGS 14D-14G shows that there were no differences in BMI, or other parameters associated with postprandial satiety when dividing participants according to their CTSGRS group.
[0046] FIG. 15 shows a flowchart diagram for the training, validation, and independent testing cohort for development of the machine-learning assisted gene risk score for calories to satiation.
[0047] FIG. 16 shows a flowchart diagram for a 12-month, placebo-controlled, randomized clinical trial with phentermine-topiramate extended-release.
[0048] FIG. 17 shows a flowchart diagram for a 16-week, placebo-controlled, randomized clinical trial with liraglutide.
[0049] FIGs 18A-18F shows weight loss response by baseline ad libitum buffet meal energy intake and machine-learning assisted gene risk score prediction in response to placebo or phentermine-topiramate ER and placebo vs. liraglutide. FIG. 18A shows weight loss in all patients that participated in the 52-week randomized clinical trial with placebo and phentermine-topiramate ER and underwent satiation testing at baseline. FIG. 18B shows weight loss by sex-stratified satiation by ad libitum buffet meal in patients assigned to placebo and phentermine-topiramate ER, where patients assigned to phentermine and high satiation by ad libitum buffet meal had the greatest total body weight lost percentage. FIG. 18C shows weight loss by sex-stratified satiation by CTSGRS classified as low or high in patients assigned to placebo and phentermine-topiramate ER, where patients assigned to phentermine-topiramate ER and high CTSGRS had the greatest total body weight lost percentage. FIG. 18D shows weight loss in all patients that participated in the 16-week randomized clinical trial of liraglutide vs placebo and underwent satiation testing at baseline. FIG. 18E shows weight loss by sex-stratified satiation by ad libitum buffet meal in patients assigned to liraglutide vs placebo, where patients assigned to liraglutide and low satiation by ad libitum buffet meal hadthe greatest total body weight lost percentage. FIG. 18F shows weight loss by sex-stratified satiation by CTSGRS in patients assigned to liraglutide vs placebo, where patients assigned to liraglutide and low satiation by CTSGRS had the greatest total body weight lost percentage.
[0050] FIGs 19A-19C shows performance of the CTSGRS model in predicting weight loss outcomes of a 12-month, placebo-controlled, randomized clinical trial with phenterminetopiramate extended-release. FIG. 19A shows the CTSGRS model had an AUC of 0.71 for predicting responders achieving at least 15% of total body weight loss at 12 months in participants that were assigned to phentermine-topiramate extended-release. FIG. 19B shows violin plots showing individual weight loss outcomes for participants that completed the 12- month trial stratified by treatment allocation and CTSGRS score. FIG. 19C shows the proportion of responders achieving a total body weight loss superior to 5, 10, and 15% by treatment allocation and CTSGRS group at 12-months.
[0051] FIGs 20A-20C shows performance of the CTSGRS model in predicting weight loss outcomes of a 16-week, placebo-controlled, randomized clinical trial with liraglutide. FIG. 20A shows the CTSGRS model had an AUC of 0.75 for predicting responders achieving at least 4% of total body weight loss at 16 weeks in participants assigned to phentermine-topiramate extended-release. FIG. 20B shows violin plots showing individual weight loss outcomes for participants that completed the 16-week trial stratified by treatment allocation and CTSGRS score. FIG. 20C shows the proportion of responders achieving a total body weight loss superior to 4, 5, and 6% by treatment allocation and CTSGRS group at 16-weeks.
[0052] FIG. 21A shows total body weight loss percentage (TBWL%) at 3, 6, 9 and 12 months, and FIG. 21B shows individual TBWL% outcomes at 12 Months. Abbreviations: TBWL%: Total body weight loss percentage; HG+: Hungry Gut Positive; HG-: Hungry Gut Negative.
[0053] FIG. 22A shows a representational overview of the genotype and CTS data collected; samples from females marked in red and from males in blue.
[0054] FIG. 22B shows a GWAS performed between genotypes in samples collected from females and CTS (w=359), while FIG. 22C shows the same as FIG. 22B, but with only male samples (n=124).
[0055] FIG. 22D shows selection of parameters using the SNPs from 41 gene (Table 5) regions and beta values from either Bf or Bmand Gene Risk Score function, where the function GRSgcomputes the gene risk score for gene g on sample S using GWAS weights B E [Bf, Bm}, kgis the number of genotyped SNPs around gene g, Contpcalculates the contribution factor using the genomic context of SNP i and the regression weight / ?£for the SNP in B, and 6s iis thedosage of risk alleles of SNP i in sample S. Here, Cont^ = if |Pos(i) — Pos(^) | < 1 x 106, otherwise 0, where Pos is the coordinate of the SNP or the coordinate of the end of gene g nearest SNP i .
[0056] FIG. 22E shows feature matrix F (w=658) containing all GRSs calculated using both Bf and Bmalong with anthropometric features A. The target matrix T contains a 1 for HS+ or 0 for HS- for each sample when CTS is greater than or less than the sex-specific CTS threshold, respectively. The FSeifeatures are chosen based on these data using backward feature selection, which iteratively removes features that cumulatively improve the AUC until no improvement is seen.
[0057] FIG. 22F shows the final model training and testing, where STe= {s |s in phentop trial}, Sv= {s|s in lira trial}, STr= {s|s in neither trial}. The CTSGRS is an SVC model trained samples STr(w=483) using the selected features FSeiagainst target variable T, oversampled using SMOTE to create an additional 225 synthetic samples (total training / / =708). Features are normalized using a Yeo power transform fit on the training features: Norm(») = Fit(FTr, A). The performance of the CTSGRS, including ROC-AUC, SN, and SP, was measured using the outputs of CTSGRS(FS6;) against T using data from samples 5Teor Sv.
[0058] FIG. 23 shows a flow chart of how participants were chosen for prospective study detailed in Example 6.
[0059] FIG. 24A shows weight change after five years follow-up based on satiation status (altered vs normal). FIG. 24B-24D shows weight change after five years follow-up based on satiation status after 5 years of follow-up by gender (FIG. 24B), and BMI classification (FIG. 24C) and age groups (FIG. 24D). Data is shown as mean and standard error of the mean.
[0060] FIG. 25 shows weight change at each follow-up year based on satiation status (altered vs normal). Data is shown as mean and standard error of the mean.
[0061] FIG. 26 shows weight change at each follow-up year based on CTSGRS (High CTSGRS vs High CTSGRS). Data is shown as mean and standard error of the mean.
[0062] FIG. 27 shows the odds ratio of upper gastrointestinal symptoms in patients with obesity and abnormal satiation.
[0063] FIGs 28A-28C shows machine learning analysis to explain the variability of caloric intake to satiation (CTS). (FIG. 28A) SHAP values from the non-genetic model highlight sex as the most influential variable, with moderate contributions from fat percentage and height, and minimal contributions from hormonal measures (postprandial PYY and GLP-1, AUC).(FIG. 28B) SHAP analysis of the genetic risk score (GRS) model identifies key genetic contributors. (FIG. 28C) Combining genetic and non-genetic features significantly improves model performance, with SHAP values showing that sex remains the dominant factor, while GRS enhances predictive accuracy alongside anthropometric features. In all panels, the color of the data points represents the direction of influence (e.g., higher or lower feature values), and the size of the points indicates the strength of each feature's impact on CTS variability. All variables were standardized, and LASSO regularization was applied to ensure robust feature selection.
[0064] FIG. 29 shows prediction of response to bariatric surgery and weight loss after one year.
[0065] FIG. 30 shows response to bariatric surgery prediction accuracy.
[0066] FIG. 31 depicts weight changes after bariatric surgery by CTSGRS.
[0067] FIG. 32 depicts TBWL% by timepoint after RYGB according to carrier state and CTSGRS model predictions. Total body weight loss percentage (TBWL%) change during 15 years after Roux-en-Y gastric bypass (RYGB). Graph represents Least Sq Mean ± SEM at each timepoint.
[0068] FIG. 33 shows weight loss outcomes for patients with high or low ML-GRS (CTSGRS) scores.
[0069] FIG. 34 shows calories to satiation (CTS) across the four quartiles for ML-GRS (restricted CTSGRS) scores.
[0070] FIG. 35 shows calories to satiation (CTS) across the four quartiles for ML-GRS (restricted CTSGRS) scores.
[0071] FIG. 36 shows linear correlation between the ML-GRS (restricted CTSGRS) score and calories to satiation (CTS) during an ad libitum buffet meal.
[0072] FIG. 37 shows sex-stratified linear correlation between the ML-GRS (restricted CTSGRS) score and calories to satiation (CTS) during an ad libitum buffet meal.
[0073] FIG. 38 shows linear correlation between the ML-GRS (restricted CTSGRS) score and gastric emptying half-life.
[0074] FIG. 39 shows sex-stratified linear correlation between the ML-GRS (restricted CTSGRS) score and gastric emptying half-life.
[0075] FIG. 40 depicts regression analysis displaying the correlation between time to 50% gastric emptying and genetic risk score for calories consumed to reach satiation (CTSGRS) in 50 participants
[0076] FIG. 41 shows total body weight loss percentage (TBWL%) of participants (HG+ and HG-) at 12 months after semaglutide initiation. Data shown as median and IQR Abbreviations: HG+= Hungry Gut Positive, HG-= Hungry Gut Negative, TBWL%= Total Body Weight Loss percentage.
[0077] FIG. 42 shows participants (HG+ and HG-) achieving a Total Body Weight Loss Percentage (TBWL%) > 5, 10, 15 and 20, at 12 months after semaglutide initiation.
[0078] FIG. 43 shows A HbAlc% of participants (HG+ and HG-) at 12 months after semaglutide initiation.DETAILED DESCRIPTIONDefinitions
[0079] While the following terms are believed to be well understood by one of ordinary skill in the art, the following definitions are set forth to facilitate explanation of the presently disclosed subject matter.
[0080] As used herein, the term “a” or “an” can refer to one or more of that entity, i.e., can refer to a plural referents. As such, the terms “a” or “an”, “one or more” and “at least one” can be used interchangeably herein. In addition, reference to “an element” by the indefinite article “a” or “an” does not exclude the possibility that more than one of the elements is present, unless the context clearly requires that there is one and only one of the elements.
[0081] Unless the context requires otherwise, throughout the present specification and claims, the word “comprise” and variations thereof, such as, “comprises” and “comprising” are to be construed in an open, inclusive sense that is as “including, but not limited to”.
[0082] The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
[0083] Reference throughout this specification to “one embodiment” or “an embodiment” means that a particular feature, structure or characteristic described in connection with the embodiment may be included in at least one embodiment of the present disclosure. Thus, the appearances of the phrases “in one embodiment” or “in an embodiment” in various places throughout this specification may not necessarily all referring to the same embodiment. It is appreciated that certain features of the disclosure, which are, for clarity, described in the context of separate embodiments, may also be provided in combination in a single embodiment. Conversely, various features of the disclosure, which are, for brevity, described in the context of a single embodiment, may also be provided separately or in any suitable sub-combination.
[0084] The terms "ad libitum diet" as used herein refer to a diet where the amount of daily calories intake of a subject is not restricted to a particular value. A subject following an ad libitum diet is free to eat till fullness.
[0085] The terms “Glucagon-like peptide- 1 receptor agonist” or “GLP-1 receptor agonist” as used herein can be used interchangeably with the terms “GLP-1 agonist” or “GLP-I analog”. Said terms can also be referred to as incretin mimetics. All of the aforementioned terms can refer to agents that act as agonists of the GLP-1 receptor and can work by activating the GLP- 1 receptor.
[0086] The term “postprandial satiety” as used herein can be interchangeably with “satiety” and refers to the sensation of fullness after a meal termination that perdures through time until hunger returns. Postprandial satiety may overlap with hunger or desire to eat.
[0087] As used herein, the term “genotype” can refer to the complete set of genetic material in an organism. Ther term “genotype” can also be used to refer to the alleles or variants an individual carries in a particular gene or genetic location.
[0088] As used herein, the term “allele(s)” can mean any of one or more alternative forms of a gene, all of which alleles relate to at least one trait or characteristic. In a diploid cell, the two alleles of a given gene occupy corresponding loci on a pair of homologous chromosomes. Since the present disclosure, in embodiments, relates to quantitative trait loci (QTLs), i.e., genomic regions that may comprise one or more genes or regulatory sequences, it is in some instances more accurate to refer to “haplotype” (i.e., an allele of a chromosomal segment) instead of “allele”, however, in those instances, the term “allele” should be understood to comprise the term “haplotype”. An “effect allele” as used herein can refer to the less common allele(s) at a particular location in the genome. As used herein, “dosage” can refer to the number of effect alleles in a genotype.
[0089] As used herein, the term “nucleotide change” can refers to, e.g., nucleotide substitution, deletion, and / or insertion, as is well understood in the art. For example, mutations contain alterations that produce silent substitutions, additions, or deletions, but do not alter the properties or activities of the encoded protein or how the proteins are made.
[0090] As used herein, the term “protein modification” refers to, e.g., amino acid substitution, amino acid modification, deletion, and / or insertion, as is well understood in the art.
[0091] As used herein, the term “phenotype” can refer to the observable characteristics of an individual cell, cell culture, organism, or group of organisms, which can result from the interaction between that individual’s genetic makeup (i.e., genotype) and the environment.
[0092] As used herein, the term “gene” can refer to any segment of DNA associated with a biological function. Thus, genes can include, but are not limited to, coding sequences and / or the regulatory sequences required for their expression. Genes can also include non-expressed DNA segments that, for example, form recognition sequences for other proteins. Genes can be obtained from a variety of sources, including cloning from a source of interest or synthesizing from known or predicted sequence information, and may include sequences designed to have desired parameters.
[0093] As used herein, the term “SNP” refers to “single nucleotide polymorphisms”, which can refer to substitution of a single nucleotide at a specific position in the genome. The specific location in the genome can possess multiple single-nucleotide alleles. The use of this term in this document should not be construed to exclude the use of any other type of polymorphism in addition to or in place of a SNP such as, for example, sequence insertions, deletions, inversions, and other sequence replacements. Single-nucleotide polymorphisms may fall within coding sequences of genes, non-coding regions of genes, or in the intergenic regions (regions between genes). SNPs within a coding sequence may not necessarily change the amino acid sequence of the protein that is produced, due to degeneracy of the genetic code. SNPs in the coding region can be of two types: synonymous SNPs and nonsynonymous SNPs. Synonymous SNPs do not affect the protein sequence, while nonsynonymous SNPs change the amino acid sequence of protein.
[0094] As used herein, the term “homologous” or “homologue” or “orthologue” is known in the art and can refer to related sequences that share a common ancestor or family member and are determined based on the degree of sequence identity. The terms “homology,” “homologous,” “substantially similar” and “corresponding substantially” are used interchangeably herein. They refer to nucleic acid fragments wherein changes in one or more nucleotide bases do not affect the ability of the nucleic acid fragment to mediate gene expression or produce a certain phenotype. These terms also refer to modifications of the nucleic acid fragments of the instant disclosure such as deletion or insertion of one or more nucleotides that do not substantially alter the functional properties of the resulting nucleic acid fragment relative to the initial, unmodified fragment. It is therefore understood, as those skilled in the art will appreciate, that the disclosure encompasses more than the specific exemplary sequences. These terms describe the relationship between a gene found in one species, subspecies, variety, cultivar or strain and the corresponding or equivalent gene in another species, subspecies, variety, cultivar or strain. For purposes of this disclosure homologous sequences are compared. “Homologous sequences” or “homologues” or “orthologues” arethought, believed, or known to be functionally related. A functional relationship may be indicated in any one of a number of ways, including, but not limited to: (a) degree of sequence identity and / or (b) the same or similar biological function. Preferably, both (a) and (b) are indicated. Homology can be determined using software programs readily available in the art, such as those discussed in Current Protocols in Molecular Biology (F.M. Ausubel et al., eds., 1987) Supplement 30, section 7.718, Table 7.71. Some alignment programs are MacVector (Oxford Molecular Ltd, Oxford, U.K.), ALIGN Plus (Scientific and Educational Software, Pennsylvania) and AlignX (Vector NTI, Invitrogen, Carlsbad, CA). Another alignment program is Sequencher (Gene Codes, Ann Arbor, Michigan), using default parameters.
[0095] The terms “substantially reduced” and “substantially less” are used interchangeably herein and, when referring to an expression level or amount or an activity level of a protein or enzyme, can refer to a lowering of said amount or activity by a percentage or range of percentages as compared to or versus a control or reference level or activity of said protein or enzyme. The terms “substantially reduced” and “substantially less” can refer to a lowering of an amount or level of a protein or enzyme or an activity of an enzyme by at least, at most, exactly or about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 21%, 22%, 23%, 24%, 25%, 26%, 27%, 28%, 29%, 30%, 30%,31%, 32%, 33%, 34%, 35%, 36%, 37%, 38%, 39%, 40%, 41%, 42%, 43%, 44%, 45%, 46%,47%, 48%, 49%, 50%, 51%, 52%, 53%, 54%, 55%, 56%, 57%, 58%, 59%, 60%, 61%, 62%,63%, 64%, 65%, 66%, 67%, 68%, 69%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%,79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%,95%, 96%, 97%, 98%, 99% or 100% as compared to or versus a control or reference (e.g., a control or reference level or activity of said protein or enzyme). The terms “substantially reduced” and “substantially less” can refer to a lowering of an amount or level of a protein or enzyme or activity of an enzyme (e.g., enzymatic activity) by l%-5%, 10%-l 5%, 15%-20%, 20%-25%, 25%-30%, 30%-35%, 35%-40%, 40%-45%, 45%-50%, 50%-55%, 55%-60%, 60%-65%, 65%-70%, 70%-75%, 75%-80%, 80%-85%, 85%-90%, 90%-95% or 95%-100%, inclusive of the endpoints, as compared to or versus a control or reference (e.g., a control or reference level or activity of said protein or enzyme). The terms “substantially reduced” and “substantially less” can also mean that the amount of a protein or enzyme or the activity of an enzyme can be at least, at most, exactly or about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%,11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 21%, 22%, 23%, 24%, 25%, 26%,27%, 28%, 29%, 30%, 30%, 31%, 32%, 33%, 34%, 35%, 36%, 37%, 38%, 39%, 40%, 41%,42%, 43%, 44%, 45%, 46%, 47%, 48%, 49%, 50%, 51%, 52%, 53%, 54%, 55%, 56%, 57%,58%, 59%, 60%, 61%, 62%, 63%, 64%, 65%, 66%, 67%, 68%, 69%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99% or 100% of the amount of control or reference version of said protein or enzyme or the activity of said enzyme. The terms “substantially reduced” and “substantially less” can also mean that the amount of a protein or enzyme or the activity of an enzyme is l%-5%, 10%-l 5%, 15%-20%, 20%-25%, 25%-30%, 30%-35%, 35%-40%, 40%-45%, 45%-50%, 50%-55%, 55%-60%, 60%-65%, 65%-70%, 70%-75%, 75%-80%, 80%-85%, 85%-90%, 90%-95% or 95%-100%, inclusive of the endpoints, of the amount of a control or reference version of said protein or enzyme or the activity of an enzyme.. With regards to a level or amount of a protein or enzyme, the control or reference can be a level or amount of said protein or enzyme in a control or reference cell. In one embodiment, the tested protein or enzyme in a control of reference cell does not have a heterologous modification. With regards to activity of an enzyme, the control or reference can be the activity of said protein or enzyme in a control or reference cell. In one embodiment, the tested protein or enzyme in a control of reference cell does not have a heterologous modification.
[0096] The level or activity of a protein or enzyme provided herein can be measured within a cell or after extraction and / or isolation from a cell (e.g., in vitro). In some cases, the level or amount of a gene encoding a protein of interest is measured or determined. The level or amount of a gene provided herein can be measured within a cell or after extraction from a cell (e.g., in vitro). In some cases, the activity of an enzyme encoded by a gene provided herein is measured or determined. The activity (e.g., specific activity) of an enzyme encoded by a gene provided herein can be measured within a cell or after extraction from a cell (e.g., in vitro). The assay utilized to measure the level or amount of expression of a gene or protein provided herein can be high-throughput in nature. The assay utilized to measure the activity of an enzyme encoded by a gene provided herein can be high-throughput in nature.
[0097] The level or amount of a gene provided herein can be measured using any assay known in the art for measuring a level or amount of a gene at the nucleic acid level. Examples of suitable assays for determining or measuring the levels of nucleic acid (e.g., a gene provided herein) can be selected from microarray analysis, RT-PCR such as quantitative RT-PCR (qRT- PCR), serial analysis of gene expression (SAGE), RNA-seq, Northern Blot, digital molecular barcoding technology, for example, Nanostring Counter Analysis, and TaqMan quantitative PCR assays. Other methods of mRNA detection and quantification can be applied, such as mRNA in situ hybridization. mRNA in situ hybridization can be measured using QuantiGeneViewRNA (Affymetrix), which uses probe sets for each mRNA that bind specifically to an amplification system to amplify the hybridization signals; these amplified signals can be visualized using a standard fluorescence microscope or imaging system. This system for example can detect and measure transcript levels in heterogeneous samples. In one embodiment, the level, presence or amount of gene or genetic variant (e.g., SNP) provided herein can be determined using an array. The array can be a “Genotyping SNP-Chip”.
[0098] The level or amount of a protein encoded by a gene provided herein can be measured using any assay known in the art for measuring a level or amount at the protein level. Examples of suitable assays for determining or measuring the levels of protein (e.g., encoded by a gene provided herein) can be selected from quantitative mass spectrometry or immunoassays including, for example, immunohistochemistry, ELISA, Western blot, immunoprecipitation, Luminex® assay, and the like, where a biomarker detection agent such as an antibody, for example, a labeled antibody, specifically binds a protein encoded by a gene provided herein and permits, for example, relative or absolute ascertaining of the amount of a protein in a sample or a cell. The level or amount of an enzyme encoded by a gene provided herein or of the gene itself that has been heterologously modified as provided herein can be compared to the level or amount of the same enzyme or gene that has not been heterologously modified as described herein and the percentage of the level or amount of the modified enzyme or gene vs. the non-modified enzyme or gene can be determined.
[0099] The activity of an enzyme encoded by a gene provided herein can be measured using any assay known in the art for measuring enzyme activity. Examples of suitable assays for determining enzyme activity can be any kinase assay known in the art such as, for example, biochemical kinase assays commercially available from EMD Millipore (e.g., FRET-based HTRF assays), eBioscience (e.g., Instant One cell signaling assays), Life Technologies (LanthaScreen or Omnia kinase assays), Symansis (e.g., Multikinase assay array), Abeam or Promega (e.g., the ADP-Glo Kinase Assay).
[0100] As used herein, the term “quantitative trait” can refer to a numerical measurement of a physical, chemical, or other aspect, trait or characteristic of an individual.
[0101] As used herein, the term “contribution factor” can refer to the apparent amount that the dosage at a particular SNP contributes to a quantitative trait.
[0102] As used herein, the terms “polygenic risk score (PRS)”; “polygenic score (PGS)” or “polygenic index (PGI)” can be used simultaneously and can refer to a numerical value indicating the magnitude of risk of a particular phenotype based upon the combined effect of SNPs from multiple locations in the genome. In some cases, the polygenicscore (PGS), polygenic risk score (PRS) or polygenic index (PGI) can be a number that summarizes the estimated effect of many genetic variants on an individual's phenotype, typically calculated as a weighted sum of trait-associated alleles.
[0103] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. Although methods and materials similar or equivalent to those described herein can be used to practice the invention, suitable methods and materials are described below. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.Overview
[0104] Provided herein are methods and systems for predicting whether or not a subject possesses a specific trait. To accomplish this, the methods and systems provided herein rely on data to up-or-down weight a genetic variant, which is located around, within or near a gene known or suspected to play in role in generating the specific trait (e.g., single nucleotide polymorphism (SNP)), contribution to the specific trait, uses systems biology to construct those weights, and uses an open number of gene risk scores (GRSs) for a set of genes known or suspected to play in role in generating the specific trait that are determined by machine learning (ML) feature selection techniques. As such, the methods and system provided herein for determining whether or not a subject possesses a specific trait do not select discrete genetic variants (e.g., SNPs). As alluded to, the genetic variants can be SNPs, but they could also be some other type of genetic variant such as, for example, an epigenetic feature. The epigenetic feature could be selected from methylation site(s), histone modification(s) or chromatin accessibility site(s). The specific trait could be a trait associated with an obesity phenotype as described herein. The specific trait associated with an obesity phenotype as provided herein could be selected from the group consisting of total kcal consumed to satiation (CTS) at an ad libitum meal, gastric emptying rate, resting energy expenditure (REE), results from a behavioral questionnaire indicative of anxiety or emotional eating and any combination thereof. The specific trait could also be, for example, weight loss or comorbidity change in response to weight loss, or response to an obesity intervention. The obesity intervention could be an anti-obesity medication (AOM), a surgical intervention (e.g., bariatric surgical intervention), a bariatric non-surgical and / or minimally invasive intervention, or a behavioral and / or lifestyle intervention. The AOM could be phentermine-topiramate or a GLP-1 agonist. The data used to up or down weight a genetic variant can be genetic and / or non-genetic data. The genetic data can include, but is not limited to, cell-type- specific epigenetic activity, gene regulatory activity, or chromatin structure and / or genomic positioning. In any method or system provided herein, a GRS is a multidimensional integration function with dimensions of genotype, trait contribution, genomic contribution as input. In any method or system provided herein, the method or system can further comprise iteratively removing genes from the set of genes whose gene risk scores are deemed not informative for the specific trait and subsequently combine the remaining GRSs with other non-genetic and / or multi-omic data for a subject with machine learning models to generate a machine-learning assisted GRS for the specific trait. The ML-GRS can then be used to predict or diagnose a de novo subject as possessing the specific trait for which the ML-GRS was constructed. The non-genetic and / or multi-omic data can be selected from the group consisting of metabolomic data, proteomic data, peptidomic data, epigenetic data, microbiome data, results from one or more questionnaires and any combination thereof. The epigenetic data can be the presence or absence of DNA modifications such as methylation, histone modifications, chromatin accessibility, etc.
[0105] In one aspect, provided herein are methods and systems for predicting whether or not a subject possesses a specific trait associated with an obesity phenotype of an individual or subject with obesity from one or more samples obtained from said individual or subject with obesity. In one embodiment, the methods and systems utilize machine learning in order to predict the obesity phenotype of the subject with obesity. In one embodiment, the methods and system utilizes contribution factors and gene risk scores (GRS) to determine the obesity phenotype of the subject with obesity. The gene risk scores can integrate expert guidance, machine learning, and non-genetic / multi-omic data. In one embodiment, the obesity phenotype of an individual with obesity can be ascertained using an adaptive, contribution factor based GRS as provided herein. The adaptive, contribution factor based GRS can be tailored to and / or used for determining or diagnosing a specific obesity phenotype such as, for example, hungry brain, hungry gut, slow bum or hedonic eating as described in US20210072259A1, which is herein incorporated by reference in its entirety. The adaptive, contribution factor based GRS can utilize expert guidance, machine learning and non-genetic / multi-omic data. The non-genetic / multi-omic data can comprise data obtained from determining the metabolome, the proteome, and / or the peptidome of an individual (e.g., an individual with obesity) of a sample obtained from the obese subject. Insome cases, the non-genetic / multi-omic data can be identified in a sample obtained from individual with obesity.
[0106] In one embodiment, the methods and systems for ascertaining an obesity phenotype of an individual or subject utilizes a regression, machine learning driven model for appetite. The regression machine learning model can be a multi-variate logistic regression, machine learning driven model. In one embodiment, the multi-variate logistic regression, machine learning driven model can be for appetite and can utilize a myriad of genetic and non- genetic multi-omic data. The myriad of genetic and non-genetic multi-omic data used in a multi-variate logistic regression machine learning-driven model for appetite provided herein can be referred to as the full phenotype tests or FPTs.
[0107] Also provided herein are methods and systems for predicting or stratifying individuals with obesity into sub-populations of responders and non-responders for one or more interventions. In one embodiment, a machine learning based system or method provided herein for use in determining the obesity phenotype of a subject or individual with obesity can be used to select said subject or individual with obesity for treatment with one or more interventions for treating obesity. In this manner, the machine learning based system or method provided herein for use in determining the obesity phenotype of a subject or individual with obesity can serve as a companion diagnostic to an intervention for treating obesity. The one or more interventions can be selected from the group consisting of a pharmacological intervention, a surgical intervention (e.g., bariatric surgical intervention), a bariatric non-surgical and / or minimally invasive intervention, a weight loss device, a diet intervention, a lifestyle intervention, a behavior intervention, and a microbiome intervention).
[0108] A lifestyle intervention for treating an obesity phenotype or weight gain as determined using any of the methods provided herein can be tailored to the obesity phenotype or a feature of said obesity phenotype (e.g., calories to satiation or accelerated gastric emptying) as described in WO2023220618A1, which is hereby incorporated by reference in its entirety for all purposes. The lifestyle intervention can be any lifestyle intervention described in WO2023220618A1, which is hereby incorporated by reference in its entirety for all purposes. For example, a lifestyle intervention useful for a subject determined to possess an abnormal satiation obesity phenotype or a high CTSGRS or high CTS using any method described herein can entail administration of a predetermined diet that is a volumetries diet comprising increased dietary fiber. In some cases, the increased dietary fiber is represented by a dietary fiber intake of greater than 3.5 g / kcal, wherein g is grams of food. In some cases, the increased dietary fiber is represented by a dietary fiber intake of greater than 4 g / kcal, wherein g is grams of food. Insome cases, the increased dietary fiber is represented by a dietary fiber content above 20 to 45 g / 10 MJ or 25 to 45 g / 10 MJ, wherein g is grams of food. In some cases, the energy density of the volumetries diet is less than 115 kcal / 100 g, less than 112 kcal / 100 g, less than 110 kcal / 100 g, less than 105 kcal / 100 g, less than 100 kcal / 100 g, less than 95 kcal / 100 g or less than 90 kcal / 100 g, wherein g is grams of food. In some cases, the energy density of the volumetries diet is between 90 and 115 kcal / 100 g, between 90 and 112 kcal / 100 g, between 90 and 110 kcal / 100 g, between 90 and 105 kcal / 100 g or between 90 and 100 kcal / 100 g, wherein g is grams of food. In some cases, the increased dietary fiber is as compared to a recommended dietary fiber intake for the subject's demographic group. In some cases, the subject is female, and the recommended dietary fiber intake is 21-25 grams of dietary fiber per day. In some cases, the subject is a male, and the recommended dietary fiber intake is 30-38 grams of dietary fiber per day. In some cases, the Caloric content of the volumetric diet comprising increased fiber content is set to 500 Calories less than the basal resting energy expenditure (REE) of the subject. In some cases, the REE of the subject is measured by indirect calorimetry. In some cases, the subject consumes 1-2 meals per day. In some cases, the subject is permitted a single serving of a fruit or vegetable in addition to each meal. In some cases, the lifestyle intervention further comprises administrating an exercise regimen to the subject. In some cases, the exercise regimen comprises a combination of resistance training and cardiovascular training for at least 150 minutes per week. In some cases, the method further comprises administering one or more pharmacological interventions based on the subject's obesity phenotype. In some cases, the one or more pharmacological interventions is phentermine-topiramate pharmacotherapy.
[0109] In one embodiment, provided herein is a method of treating obesity in a subject in need thereof, the method comprising: (a) determining whether or not the subject possesses an abnormal satiation or hungry brain phenotype; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the obesity analyte signature detected in the sample obtained from the subject indicates that the subject possesses an abnormal satiation phenotype or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the obesity analyte signature detected in the sample obtained from the subject indicates that the subject does not possess an abnormal satiation phenotype. In some cases, step (b) comprises administering phentermine-topiramate pharmacotherapy if the obesity analyte signature detected in the sample obtained from the subject indicates that the subject possesses an abnormal satiation phenotype. In some cases, step (a) is performed by detecting an obesity analyte signature in a sample obtained from a subject; wherein the obesity analyte signature is indicative of an obesity phenotype of the subject. In some cases, step (a) is performed using aML-GRS developed for predicting the abnormal satiation or hungry brain phenotype as provided herein. In some cases, step (a) is performed using a machine learning, full phenotype test for abnormal satiation or hungry brain phenotype as provided herein. The GLP-1 agonist can be any agent known in the art whose mechanism of action comprises activating the GLP- 1 receptor. GLP-1 agonists for use in a method and / or system provided herein can include, but are not limited to, exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SALO 112, GS-4571, ASC- 30, TTP273, MLX-7006, BMF-650 or any combination thereof.
[0110] In one embodiment, provided herein is a method of identifying a subject suffering from obesity as a non-responder to treatment with a GLP-1 agonist, the method comprising detecting an obesity phenotype of the subject; and identifying the subject as a non- responder to GLP-1 agonist treatment if the subject is determined to possess an abnormal satiation obesity phenotype. In some cases, a non-responder to GLP-1 agonist treatment may be a responder to phentermine-topiramate pharmacotherapy if the subject is determined to possess an abnormal satiation obesity phenotype. In one embodiment, the detection of the obesity phenotype of the subject can comprise utilizing a machine learning based system or method provided herein for use in determining the obesity phenotype of a subject or individual with obesity. In one embodiment, the detection of the obesity phenotype of the subject is performed using a ML-GRS developed for predicting a specific trait associated with the abnormal satiation or hungry brain phenotype as provided herein. The specific trait can be total kcal consumed to satiation (CTS) at an ad libitum meal. In one embodiment, the detection of the obesity phenotype of the subject is performed using by determining using a machine learning, full phenotype test for abnormal satiation or hungry brain phenotype as provided herein. The GLP-1 agonist can be any agent known in the art whose mechanism of action comprises activating the GLP-1 receptor. GLP-1 agonists for use in a method and / or system provided herein can include, but are not limited to, exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601,ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006, BMF-650 or any combination thereof.
[0111] In another embodiment, provided herein is a method of predicting one or more adverse effects (AEs) in a subject following administration of an anti-obesity medicine (AOM). In some cases, the AOM is administered alone. In some cases, the AOM was administered with one or more additional therapeutic agents. In some cases, the method of predicting one or more adverse effects (AEs) in a subject following administration of an anti-obesity medicine (AOM) alone or in combination with one or more additional therapeutic agents comprises determining the CTSGRS of a subject using any of the method provided herein. If the subject has a high CTSGRS (i.e., hungry gut positive (HG- or ML-GRS>0.5)) then that subject is predicted to experience adverse effects when administered a GLP-1 receptor agonist as compared to a subject with a low CTSGRS (i.e., hungry gut positive (HG+ or ML-GRS<0.5)). The one or more AEs is selected from the group consisting of diarrhea, constipation, abdominal discomfort, nausea, vomiting, headache and any combination thereof. In some cases, the AE is nausea. In some cases, the AE is nausea and headache. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes. In some cases, subjects being administered a GLP-1RA who are found to a possess a high CTSGRS alter their use of said GLP- 1RA. The alteration can be permanent or temporary discontinuation. The alteration can be to alter the frequency of administration. For example, subjects being administered a GLP-1RA who are found to a possess a high CTSGRS may alter their use from daily administration to lower frequency. The lower frequency may be weekly or bi-weekly administration. In some cases, the AOM is a GLP-1 receptor agonist. The GLP-1 receptor agonist can be selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SALO 112, GS-4571, ASC- 30, TTP273, MLX-7006 and BMF-650. In some cases, the one or more additional therapeutic agents can be an antihyperglycemic medications (OAMs). The 0AM can be insulin. In some cases, the one or more additional therapeutic agents can be a sodium-glucose co-transporter-2 (SGLT-2) inhibitor. The SGLT-2 inhibitor can be selected from the group consisting of empagliflozin, canagliflozin, dapagliflozin, ertugliflozin, bexagliflozin and sotagliflozin. In some cases, the one or more additional therapeutic agents is a dipeptidyl peptidase-4 (DPP -4)inhibitor. The DPP -4 inhibitor can be selected from the group consisting of sitagliptin, linagliptin, saxagliptin and alogliptin.
[0112] In one embodiment, provided herein is a method for predicting a response to a bariatric procedure or any surgical intervention that changes stomach size or shape that comprises determining the CTSGRS of a subject and selecting the subject as a responder or nonresponder to the bariatric procedure or the surgical intervention that changes stomach size or shape based on the determined CTSGRS of the subject. The CTSGRS of the subject can be determined using any method provided herein for determining a CTSGRS. The bariatric intervention can be any invasive or non-invasive bariatric intervention known in the art and / or provided herein. In some cases, the bariatric procedure is a bariatric surgical intervention. The bariatric surgical intervention can be selected from the group consisting of sleeve gastrectomy, Roux-en-Y gastric bypass, adjustable gastric banding (AGB), biliopancreatic diversion with duodenal switch (BPD-DS), single-anastomosis duodeno-ileostomy with sleeve (SADI-S), one-anastomosis gastric bypass (OAGB), sleeve gastrectomy with transit bipartition (SG-TB) and single-anastomosis sleeve ileal (SASI) bypass. In some cases, the bariatric procedure is a bariatric non-surgical and / or minimally invasive intervention. The bariatric non-surgical and / or minimally invasive intervention can be selected from the group consisting of endoscopic sleeve gastroplasty (ESG), gastric plication, intragastric balloons, AspireAssist and bypass liners. In some cases, the bariatric non-surgical and / or minimally invasive intervention is ESG. In some cases, the response to the bariatric procedure or the surgical intervention that changes stomach size or shape is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following the bariatric procedure or the surgical intervention that changes stomach size or shape. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0113] In one embodiment, provided herein is a method for treating obesity, the method comprising: determining the CTSGRS of a subject; selecting the subject as a responder or nonresponder to a bariatric procedure or a surgical intervention that changes stomach size or shape based on the determined CTSGRS of the subject; and performing or administering the bariatric procedure or a surgical intervention that changes stomach size or shape based on the determined CTSGRS of the subject. The CTSGRS of the subject can be determined using any methodprovided herein for determining a CTSGRS. The bariatric intervention can be any invasive or non-invasive bariatric intervention known in the art and / or provided herein. In some cases, the bariatric procedure is a bariatric surgical intervention. In one embodiment, if the subject is determined to have a high CTSGRS as described herein, the subject is selected for a bariatric surgical intervention and the bariatric surgical intervention is performed on the subject. The bariatric surgical intervention can be selected from the group consisting of sleeve gastrectomy, Roux-en-Y gastric bypass, adjustable gastric banding (AGB), biliopancreatic diversion with duodenal switch (BPD-DS), single-anastomosis duodeno-ileostomy with sleeve (SADI-S), one-anastomosis gastric bypass (OAGB), sleeve gastrectomy with transit bipartition (SG-TB) and single-anastomosis sleeve ileal (SASI) bypass. In some cases, the bariatric procedure is a bariatric non-surgical and / or minimally invasive intervention. The bariatric non-surgical and / or minimally invasive intervention can be selected from the group consisting of endoscopic sleeve gastroplasty (ESG), gastric plication, intragastric balloons, AspireAssist and bypass liners. In some cases, the bariatric non-surgical and / or minimally invasive intervention is ESG. In one embodiment, if the subject is determined to have a low CTSGRS as described herein, the subject is selected for ESG and the ESG is performed on the subject. In some cases, the response to the bariatric procedure or the surgical intervention that changes stomach size or shape is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following the bariatric procedure or the surgical intervention that changes stomach size or shape. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0114] In another aspect, provided herein is a method of predicting response to a bariatric surgical intervention in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a high CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to respond to the bariatric surgical intervention as compared to a subject with a low CTSGRS. In some cases, In some cases, the method further comprises determining the presence or absence of genetic variants in components of the leptin-melanocortin pathway (LMP). In some cases, a subject is selected for treatment with RYGB if said subject is determined to lack genetic variants in the LMP pathway and possess a low CTSGRS. the response to the bariatric surgical intervention isevidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4 or 5 years following the bariatric surgical intervention. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0115] In another aspect, provided herein is a method of treating obesity in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein ahigh CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to respond to a bariatric surgical intervention as compared to a subject with a low CTSGRS;; performing the bariatric surgical intervention on the subject if the subject is determined to possess a high CTSGRS. The bariatric surgical intervention can be selected from the group consisting of sleeve gastrectomy, Roux-en-Y gastric bypass (RYGB), adjustable gastric banding (AGB), biliopancreatic diversion with duodenal switch (BPD-DS), singleanastomosis duodeno-ileostomy with sleeve (SADI-S), one-anastomosis gastric bypass (OAGB), sleeve gastrectomy with transit bipartition (SG-TB) and single-anastomosis sleeve ileal (SASI) bypass. In some cases, the response to the bariatric surgical intervention is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4 or 5 years following the bariatric surgical intervention. In some cases, the method further comprises determining the presence or absence of genetic variants in components of the leptin-melanocortin pathway (LMP). In some cases, a subject is selected for treatment with RYGB if said subject is determined to lack genetic variants in the LMP pathway and possess a low CTSGRS. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0116] In another aspect, provided herein is a method of predicting response to endoscopic sleeve gastroplasty in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a low CTSGRS indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to endoscopic sleeve gastroplasty as compared to a subject with a high CTSGRS. In some cases,the response to endoscopic sleeve gastroplasty is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following endoscopic sleeve gastroplasty, endoscopy or surgical intervention that changes stomach size or shape. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0117] In another aspect, provided herein is a method of treating obesity in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a low CTSGRS indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to ESG as compared to a subject with a high CTSGRS;; performing the ESG on the subject if the subject is determined to possess a low CTSGRS. In some cases, the response to endoscopic sleeve gastroplasty is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following endoscopic sleeve gastroplasty, endoscopy or surgical intervention that changes stomach size or shape. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subj ect has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0118] A subject or individual can be a mammal. Any type of mammal can be assessed and / or treated as described herein. Examples of mammals that can be assessed and / or treated as described herein include, without limitation, primates (e.g., humans and monkeys), dogs, cats, horses, cows, pigs, sheep, rabbits, mice, and rats. In some cases, the mammal can be a human. In some cases, a mammal can be a mammal with obesity. For example, humans with obesity can be assessed for intervention (e.g., lifestyle intervention and / or a pharmacological intervention) responsiveness and treated with one or more interventions as described herein. In cases where the mammal is a human, the human can be of any race. For example, a human can be Caucasian or Asian.
[0119] Any appropriate method can be used to identify a mammal or individual as being overweight (e.g., as having obesity). In some cases, calculating body mass index (BMI),measuring waist and / or hip circumference, health history (e.g., weight history, weight-loss efforts, exercise habits, eating patterns, other medical conditions, medications, stress levels, and / or family health history), physical examination (e.g., measuring your height and examining your abdomen), percentage of body fat and distribution, percentage of visceral and organs fat, metabolic syndrome, and / or obesity related comorbidities can be used to identify mammals (e.g., humans) as having obesity. For example, a BMI of greater than about 30 kg / m2can be used to identify mammals (e.g., Caucasian humans) as having obesity with or without a comorbidity such as type 2 diabetes. For example, a BMI of greater than about 27 kg / m2with a co-morbidity can be used to identify mammals (e.g., Asian humans) as having obesity.
[0120] The individual can also have one or more obesity-related (e.g., weight-related) co-morbidities. Examples of weight-related co-morbidities include, without limitation, hypertension, type 2 diabetes, dyslipidemia, obstructive sleep apnea, gastroesophageal reflux disease, weight baring joint arthritis, cancer, non-alcoholic fatty liver disease, nonalcoholic steatohepatitis, depression, anxiety, and atherosclerosis (coronary artery disease and / or cerebrovascular disease). In some cases, the methods and materials described herein can be used to treat one or more obesity-related co-morbidities.
[0121] A sample can be any type of sample that can be obtained from a subject. In some cases, a sample can be a biological sample. In some cases, a sample can contain obesity analytes (e.g., DNA, RNA, proteins, peptides, metabolites, hormones, and / or exogenous compounds (e.g., medications). Examples of samples that can be assessed as described herein include, without limitation, fluid samples (e.g., blood, serum, plasma, urine, saliva, sweat, or tears), breath samples, cellular samples (e.g., buccal samples), tissue samples (e.g., adipose samples), stool samples, gastrointestinal mucosa samples. In some cases, a sample (e.g., a blood sample) can be collected while the mammal is fasting (e.g., a fasting sample such as a fasting blood sample). In some cases, a sample can be processed (e.g., to extract and / or isolate obesity analytes).Gene Risk Score Model
[0122] In one embodiment, the obesity phenotype of an individual can be ascertained using an adaptive, contribution factor based gene risk score (GRS) as provided herein. In some cases, the individual has been diagnosed with obesity. In some cases, the individual has not been diagnosed with obesity. In some cases, the individual has a BMI of greater than about 30 kg / m2with or without a co-morbidity such as type 2 diabetes. In some cases, the individual has a BMI of greater than about 27 kg / m2with a co-morbidity. In some cases, the individual has aBMI of less than or equal to about 30 kg / m2with or without a co-morbidity such as type 2 diabetes. A GRS as provided herein can utilize a basic regression model or a statistical model as known in the art. In some cases, an adaptive, contribution factor based GRS can be generated for a specific obesity phenotype. The specific obesity phenotype can be any obesity phenotype known in the art, such as, for example hungry brain, hungry gut, slow burn or hedonic / emotional eating as described in US20210072259A1, which is herein incorporated by reference in its entirety. In some cases, provided herein is a suite of adaptive, contribution factor based GRSs that can integrate data obtained from a sample obtained from a subject with obesity to ascertain said subject’s obesity phenotype. The suite of adaptive, contribution factor based GRSs can comprise adaptive, contribution factor based GRSs developed for each obesity phenotype (e.g., hungry brain, hungry gut, slow bum or hedonic eating). The suite of adaptive, contribution factor based GRSs can be part of a system provided herein for determining a subject’s obesity phenotype. Any obesity phenotype specific, adaptive, contribution factor based GRS can utilize expert guidance, machine learning and non-genetic / multi-omic data. The non-genetic / multi-omic data can comprise measuring or assessing the metabolome, the proteome, and / or the peptidome of an individual (e.g., an individual with obesity). In some cases, the non-genetic / multi-omic data can be identified in a sample obtained from an individual with obesity. The multi-omic data obtained from a sample obtained from an individual with obesity can represent an obesity analyte signature for that individual.
[0123] In one aspect, provided herein is a method of predicting or diagnosing an obesity phenotype of subject suffering from or suspected of suffering from obesity that utilizes a machine learning assisted gene risk score (ML-GRS). In some cases, the method of predicting or diagnosing an obesity phenotype of a subject suffering from or suspected of suffering from obesity that utilizes a ML-GRS is a computer implemented method. In some cases, the method of predicting or diagnosing an obesity phenotype of a subject suffering from or suspected of suffering from obesity that utilizes a ML-GRS is an in silico method. The GRS for use in this method can be an adaptive, contribution factor based GRS that can utilize expert guidance for selecting specific genes of interest for calculating gene risk scores, machine learning and non- genetic / multi-omic data as provided herein (e.g., Example 4). A subject can be determined to be suffering from obesity if they are determined to possess a BMI of greater than about 30 kg / m2 with or without a co-morbidity or a BMI of greater than about 27 kg / m2 with a comorbidity. Generation of a ML-GRS for use in predicting an obesity phenotype of interest (e.g., hungry brain (abnormal satiation), hungry gut (abnormal satiety), slow burn (slow metabolism) or hedonic / emotional eating) in a subject suffering from obesity can comprise: (a) receiving ina computer system, a genetic dataset comprising a plurality of single nucleotide polymorphisms (SNPs) located in, around or near a set of genes obtained from samples obtained from each subject from a population of subjects, wherein each subject in the population of subjects is obese, and wherein each gene in the set of genes is known or suspected to play a role in obesity; (b) generating by the computer system, a contribution factor for each SNP located in, around or near a gene from the set of genes from the genetic dataset for a subject from the population of subjects; (c) calculating by the computer system, a gene risk score (GRS) for each gene in the set of genes by integrating the contribution factors generated in step (b) for each SNP for each gene from the set of genes and SNP genotype information located in, around, or near each gene in the set of genes, wherein the GRS for each gene represents the cumulative effect of the variation in SNPs located in, around, or near the gene in that particular subject has on the function or regulation of the gene so as to increase or decrease the risk of the subject possessing a specific trait known to be associated with the obesity phenotype of interest; (d) iteratively performing by the computer system, steps (a)-(c) for each additional subject from the population of subjects to generate a set of GRSs for each gene in the set of genes from the population of subjects, wherein each gene has a GRS for each sex; (e) integrating by the computer system, the GRSs from step (d) with non-genetic and / or multi-omic data, which can be any such data provided herein and / or known in the art, obtained from the samples obtained from each subject from the population of subjects using a machine learning model to calculate a machine learning assisted gene risk score (ML-GRS) for the specific trait known to be associated with the genotype of interest; (f) training by the computer system, the ML-GRS from step (e) using a classifier model to predict subjects from a training set of samples from a training population of obese subjects as being positive for the specific trait known to be associated with the obesity phenotype of interest, thereby generating a trained ML-GRS. In some cases, the method further comprises: (g) determining by the computer system, if a test subject is positive for the specific trait known to be associated with the obesity phenotype of interest by: (i) calculating by the computer system the ML-GRS for the test subject using steps (a)-(e) on a sample obtained from the test subject; and (ii) comparing by the computer system, the ML-GRS for the test subject to the trained ML-GRS. The test subject can be determined to be negative for the specific trait known to be associated with the obesity phenotype of interest if the ML-GRS for the test subject is below a sex-specific percentile. In some cases, the sexspecific percentile is at least the 60th, 65th, 70th, 75th, 80th, 85th, 90th, 95thor 99thpercentile for the specific trait known to be associated with the obesity phenotype of interest. In some cases, the sex-specific percentile is between the 60th-65th, 65th-70th, 70th-75th, 75th-80th, 80th-85th, 85th-90th, 90th-95thor 95th-99thpercentile for the specific trait known to be associated with the obesity phenotype of interest. The test subject can be determined to be negative for the specific trait known to be associated with the obesity phenotype of interest if the ML-GRS for the test subject is below a sex-specific percentile. In some cases, the sex-specific percentile is at least the 10th, 15th, 20th, 25th, 30th, 35th, 40th, 45th, 50th, 55thor 60thpercentile for the specific trait known to be associated with the obesity phenotype of interest. In some cases, the sex-specific percentile is between the 10th- 15th, 15th-20th, 20th-25th, 25th-30th, 30th-35th, 35th-40th, 40th-45thor 45th-50thpercentile for the specific trait known to be associated with the obesity phenotype of interest. In some cases, the test subject can be determined to be positive for the specific trait known to be associated with the obesity phenotype of interest if the ML-GRS for the test subj ect is above the sex-specific 75thpercentile for the specific trait known to be associated with the obesity phenotype of interest or negative for the specific trait known to be associated with the obesity phenotype of interest if the ML-GRS for the subject is below the sex-specific 75thpercentile for the specific trait known to be associated with the obesity phenotype of interest. In this way, the final ML-GRS model can provide a binary decision on whether or not the test subject possesses the specific trait or not. The final ML-GRS model or classifier can generate a value ranging from 0 to 1, wherein a value of greater than or equal to 0.5 classifies the subject as having a high GRS for the specific trait and a value of less than 0.5 classifies the subject as having a low GRS for the specific trait, wherein a high GRS indicates a probability of being above the sex-specific 75thpercentile for the specific trait and a low GRS indicates a probability of being below the sex-specific 75th percentile for the specific trait. In some cases, the SNP genotype information can comprise risk allele dosage, major allele presence, occurrence of a de novo variant, insertion, deletion, or genome rearrangement. In some cases, each gene in the set of genes is known or suspected to play a role in appetite regulation, energy expenditure, lipid metabolism or adipogenesis.
[0124] In some cases, the classifier model in step (e) can be selected from the group consisting of a support vector classifier (SVC), gradient boosted machines (GBM), Neural networks (NN), convolutional neural networks (CNN), Logistic regression (e.g., LASSO and Ridge regression), classification and regression trees (CRT), Bayesian inference, and clustering methods.
[0125] The contribution factor can be a numerical value, vector, matrix, or function of various biological factors that represents the contribution or effect that each SNP has on a specific trait known to be associated with the obesity phenotype of interest. The various biological factors can be non-genetic and / or multi-omic data. The non-genetic and / or multiomic data can be selected from the group consisting of metabolomic data, proteomic data, peptidomic data, epigenetic data, microbiome data, results from one or more questionnaires and any combination thereof. The epigenetic data can comprise a presence or absence of DNA modifications selected from the group consisting of methylation modifications, histone modifications and chromatin accessibility. The microbiome data can be from any site on a subject or a subject’s environment. For example, but in no way limiting, the microbiome data can be from a stool sample, a skin sample, or a sample from the subject’s environment.
[0126] In some cases, the contribution factor is a function or algorithm that utilizes output of a regression analysis or a statistical analysis on the genetic dataset as well as one or more additional co-factors, covering regression against one or more traits using one or more types of genetic datasets selected from the group consisting of a genome-wide data microarray chip (GWAS), targeted sequencing (exome or targeted genetic panel), variants detected by qPCR through targeted amplification, whole genome sequencing (WGS), targeted or untargeted sequencing of genomic re-arrangement, deletions, duplications, repeat extensions, and complex genotyping of HLA and CYP genes. In some cases, the results from multiple regressions can be combined into a single contribution factor. In some cases, the one or more additional co-factors relate to biological factors regarding genetics, epigenetics, and gene regulation for each SNP. In some cases, the one or more additional co-factors are selected from the group consisting of proximity of the SNP to the specific gene, known or theoretical mechanistic role of the SNP metadata regarding the SNP in relationship to the gene, and any combination thereof. The known or theoretical mechanistic role of the SNP can be selected from the group consisting of synonymous mutation, nonsynonymous mutation, frameshift mutation and nonsense mutation. The known or theoretical mechanistic role of the SNP can affect the three-dimensional structure of the translated protein binding sites, protein-protein interaction sites and modification sites. The modification sites can be selected from the group consisting of phosphorylation, glycosylation and proteolytic site. In some cases, the mechanistic role comprises whether or not the SNP is within an intron, exon, enhancer region or regulatory region of the gene, a cis-regulatory element, a promoter region, a non-coding exonic region, a coding exonic region, intronic region, splice site, transcription factor binding site, epigenetic modification site, at a remote genomic location involved in a three-dimensional chromatin contact with the gene, or a cell-type-specific topological-associated domain that contains the gene. The cis-regulatory element can be an enhancer or insulator. In some cases, the epigenetic modification site comprises methylation modifications, histone modifications, and / or chromatin accessibility.
[0127] In some cases, any method or system provided herein for generating a ML-GRS step (e) comprises a forward feature selection, backward feature selection, or random feature sampling algorithm and a random forest, GBM, or linear predictor that iteratively selects a feature with a desired training accuracy using the random forest predictor in multiple rounds until no further training accuracy improvement is detected.
[0128] In some cases, each GRS generated in a method or system provided herein for generating a ML-GRS provides an indication of a role or effect that the SNPs located around or near a specific gene may play on the specific trait known to be associated with the obesity phenotype of interest.
[0129] In some cases, any method or system provided herein for generating a ML-GRS further comprises performing a genetic analysis on a sample obtained from the subject suffering from or suspected of suffering from obesity prior to step (a). In some cases, the genetic analysis comprises obtaining sequence reads from the whole or portions of the whole genome of the subject. In some cases, the genetic analysis comprises performing a genotyping method selected from the group consisting of restriction fragment length polymorphism identification (RFLPI), random amplified polymorphic detection (RAPD), amplified fragment length polymorphism detection (AFLPD), polymerase chain reaction (PCR), DNA sequencing, RNA sequencing, allele specific oligonucleotide (ASO) probes, and hybridization to microarrays or beads.
[0130] In some cases, any method or system provided herein for generating a ML-GRS further comprises performing by the computer system, a regression analysis and / or statistical analysis on the genetic dataset from step (a) prior to step (b). In some cases, the regression analysis is a ridge regression or least absolute shrinkage and selection operator (LASSO) regression. In some cases, the beta values obtained from the regression analyses performed on the genetic datasets are used to generate the contribution factor in step (b). In some cases, the statistical analysis is a genome-wide association study (GWAS). In some cases, p-values and / or regression weights obtained from the statistical analyses performed on the genetic datasets are used to generate the contribution factors in step (b). In some cases, performing the regression analysis or statistical analysis on the genetic dataset occurs following the genetic analysis on the sample as described herein such that the genetic dataset obtained from the genetic analysis (e.g., genotyping method) is used for the regression or statistical analysis. In some cases, the method for generating the ML-GRS further comprises performing a regression analysis on the genetic dataset obtained from the sample obtained from the subject prior to generating the contribution factor. The regression analysis can be any regression analysis known in the art andappropriate for the variables (i.e., predictor variable(s) and response variable) utilized in the analysis. In one embodiment, the regression analysis is a ridge regression or least absolute shrinkage and selection operator (LASSO) regression. In one embodiment, beta values obtained from the regression analyses (e.g., LASSO regression) performed on the genetic datasets are used to generate the contribution factors as described herein (e.g., Example 4). In some cases, the method for generating the ML-GRS further comprises performing a statistical analysis on the genetic dataset obtained from the sample obtained from the subject prior to generating the contribution factor. The statistical analysis can be any statistical analysis known in the art and appropriate for determining the correlation between SNPs and a particular obesity phenotype (e.g., hungry brain, hungry gut, emotional eating and slow burn). In one embodiment, the statistical analysis is a genome-wide association study (GWAS). In one embodiment, p-values obtained from the statistical analyses (e.g., GWAS) performed on the genetic datasets are used to generate the contribution factors as described herein.
[0131] In one embodiment, the contribution factor and / or gene risk score can be determined or calculated using the methods / formulas depicted in FIG. 22A-22E. In some cases, the prior to determining a contribution factor and / or gene risk score, a first GWAS can be performed between genotypes in samples collected from females and a specific trait (e.g., CTS) and a second GWAS can be performed between genotypes in samples collected from males and a specific trait (e.g., CTS). Parameters can then be selected using the SNPs from expert guided or candidate genes (e.g., 41 gene regions for CTSGRS) and beta values from either Bf or Bmand Gene Risk Score function. In some cases, the Gene Risk Score function GRSgcomputes the gene risk score for gene g on sample S using GWAS weights BG{Bf, Bm}, kgis the number of genotyped SNPs around gene g. Contgcan calculate the contribution factor for a gene in any method provided herein using the genomic context of SNP i and the regression weight βifor the SNP in B, and 5s, i which can be the dosage of risk alleles of SNP i in sample S. Further to this, Contg=Pi if |Pos(i)-Pos(g)|< 1 x 106, otherwise 0, where Pos is the coordinate of the SNP or the coordinate of the end of gene g nearest SNP i. Any method provided herein for generating a ML GRS can further comprise a Feature Matrix. The feature matrix F can contain all GRSs calculated using both Brand Bmalong with anthropometric features A (e.g., age, weight, height and / or sex). The target matrix T can contain a 1 for being positive for a specific trait (e.g., HS+ for abnormal satiation) or 0 for being negative for a specific trait (e.g., HS- for abnormal satiation) for each sample when the specific trait (e.g., CTS) is greater than or less than the sex-specific specific trait (e.g., CTS) threshold, respectively. The Fsei features can be chosen based on these data using backward feature selection, which begins using allfeatures and iteratively removes features whose removal maximally improves the AUC until no features remain. The selected features Fsei can be those that produced the maximum AUC.
[0132] It should be noted that an alternative for using SNPs in a method for generating an ML-GRS for a specific trait or phenotype of a subject (e.g., a specific trait associated with an obesity phenotype of interest) as provided herein can entail using epigenetic data form the subject or population of subjects. Examples of epigenetic data that can be used include, but is not limited to, epigenetic data on histone modification sites. Histone modifications that can used to generate a ML-GRS as provided herein can be selected from the group consisting of acetylation, methylation, phosphorylation, ubiquitylation, GlcNAcylation, citrullination, crotonylation, sumoylation, isomerization and combinations thereof.
[0133] The ML-GRS in any method provided herein can predict if any one subject possesses the specific trait known to be associated with the obesity phenotype of interest by combining gene risk scores and additional biological, psychological, or environmental measurements of the any one subject. The additional biological, psychological or environmental measurements can include any non-genetic and / or multi-omic data provided herein.
[0134] In some cases, the gene risk scores can be normalized prior to training the SVC. The normalization can be performed using a Yeo power transform and can entail creating synthetic samples for training using, for example, SMOTE sampling. In one embodiment, the GSRs are normalized with a Yeo power transformation and the height of the individual using a support vector classifier (SVC). The SVC was trained using a linear kernel and a C value of 0.001 using SMOTE oversampling.
[0135] In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a sensitivity and / or specificity is at least 65%. In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a sensitivity and / or specificity is at least 72%. In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with an area under the curve (AUC) of at least 0.60, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9 0.95 or 1. In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with positive predict value of at least 50%. The prediction can be as compared to a control. In one embodiment, the control is a patient or population of patients previously determined to possess the specific obesity phenotype the ML-GRS is developed to predict using any method known in the art such as those taught in US20210072259A1 and / or WO2022246203A1, both which are herein incorporated by reference in their entireties. The obesity of interest can be selected from thegroup consisting of hungry brain (abnormal satiation), hungry gut (abnormal satiety), slow burn (slow metabolism) and hedonic / emotional eating.
[0136] As noted previously herein, the specific trait on which an ML-GRS can be built as described herein can vary. In some cases, the ML-GRS can be built around response to an obesity treatment. The obesity treatment could be a pharmacological or non-pharmacological treatment. The pharmacological treatment can be response to phen-top treatment or response to GLP-1 agonist treatment.
[0137] In one embodiment, the obesity phenotype of interest is the hungry brain or abnormal satiation obesity phenotype. Further to this embodiment, the specific trait in step (b) is consuming above the seventy-fifth percentile amount of calories (e.g., kcal) at an ad libitum meal (CTS). The set of genes in step (a) comprises or consists of the genes listed in Table 5. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MBOAT4 (GOAT), ADRA2A, GNB3, GPBAR1, LEP, LEPR, SH2B1, SIM1, NCOA1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO 1, FANCL, CADM2, PTBP2, NUDT3, OLFM4, HOXB5, NPY1R, NPY2R, NPY4R, NPY5R, NR1H4, SLC6A4, NTS, UCP3, ADIPOQ, CCK, CCKR, DPP4, FGFR4, GHRL, MC4R, MC3R, POMC, AGRP, GIP, GIPR, CNR1, FAAH, GCG, CELA2A, PPARG, TFAP2B, A PO E, TNFRSF11 A or DYRK1B. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MBOAT4 (GOAT), ADRA2A, GNB3, LEPR, SH2B1, SIM1, NCOA1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO1, FANCL, CADM2, PTBP2, NUDT3, OLFM4, HOXB5, NPY1R, NPY5R, NR1H4, UCP3, FGFR4, GHRL, MC4R, POMC, AGRP, GIPR or TFAP2B. In some cases, steps (a)-(d) produce a defined set of GRSs across the set of genes. In one embodiment, steps (a)-(d) produce GRSs across the set of genes that optimally maximize the true positive rate while minimizing the false positive rate for predicting the hungry brain phenotype. In one embodiment, steps (a)-(d) produce GRSs across a plurality of genes that are predictive of the hungry brain phenotype with the desired sensitivity, specificity, AUC and / or PPV as provided herein. In some cases, the plurality of genes for which steps (a)- (d) produce a defined set of GRSs consists essentially of, consists of or comprises SIM1, PCSK1, SH2B1, LEPR, UCP2, FTO, TCF7L2, GLP1R, TNFRSF11A, MC4R, CCK, CCKR, MC3R, GIP, GIPR and ADRA2A or any combination thereof. A ML-GRS generated for determining if a test subject possesses abnormal satiation (i.e., high CTS as described herein such as in Example 4) can be referred to as CTSGRS. In one embodiment, the obesity phenotypeof interest is the hungry gut or abnormal satiety obesity phenotype. Further to this embodiment, the specific trait in step (b) is an increased or accelerated baseline gastric emptying for the subject as compared to a control (e.g., individual without obesity. The gastric emptying can be measured using any method known in the art such as, for example, scintigraphy and can be represented as GE Tl / 4 or GE Tl / 2. The set of genes in step (a) comprises or consists of genes known to play role in conferring the hungry gut phenotype. The set of genes in step (a) comprises or consists of the genes listed in Table 5. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MB0AT4 (GOAT), ADRA2A, GNB3, GPBAR1, LEP, LEPR, SH2B1, SIM1, NC0A1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LING01, FANCL, CADM2, PTBP2, NUDT3, 0LFM4, H0XB5, NPY1R, NPY2R, NPY4R, NPY5R, NR1H4, SLC6A4, NTS, UCP3, ADIPOQ, CCK, CCKR, DPP4, FGFR4, GHRL, MC4R, MC3R, POMC, AGRP, GIP, GIPR, CNR1, FAAH, GCG, CELA2A, PPARG, TF AP2B, APOE, TNFR.SF11 A or DYRK 1B. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MB0AT4 (GOAT), ADRA2A, GNB3, LEPR, SH2B1, SIM1, NC0A1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO 1, FANCL, CADM2, PTBP2, NUDT3, 0LFM4, H0XB5, NPY1R, NPY5R, NR1H4, UCP3, FGFR4, GHRL, MC4R, POMC, AGRP, GIPR or TFAP2B. In some cases, steps (a)-(d) produce a defined set of GRSs across a plurality of genes. In one embodiment, steps (a)-(d) produce GRSs across a plurality of genes that optimally maximize the true positive rate while minimizing the false positive rate for predicting the hungry gut phenotype. In one embodiment, steps (a)-(d) produce GRSs across a plurality of genes that are predictive of the hungry gut phenotype with the desired sensitivity, specificity, AUC and / or PPV as provided herein.
[0138] In one embodiment, provided herein is a computer-implemented method for predicting a subject’s calories to satiation (CTS), the method comprising: (a) receiving in a computer system, genetic data obtained from a sample obtained from a subject as well as anthropometric data from the subject; (b) calculating by the computer system, a gene risk score (GRS) for each respective gene in a set of genes from the genetic data obtained from the sample obtained from the subject, wherein the GRS for each respective gene in the set of genes represents a cumulative effect of single nucleotide polymorphisms (SNPs) located in, around, or near the respective gene has on the function or regulation of the respective gene in the subj ect so as to increase or decrease the risk of the subject being above a sex-specific threshold of total kcal consumed to satiation (CTS) at an ad libitum meal, wherein the set of genes comprises,consists essentially of or consists of the set of genes listed in Table 5; (c) normalizing by the computer system, the anthropometric data and the GRS for each respective gene from the set of genes to training samples from a CTSGRS classifier to generate normalized GRSs and anthropometric data for the subject; and (d) inputting by the computer system, the normalized GRSs and anthropometric data for the subject into the CTSGRS classifier to generate a value ranging from 0 to 1, wherein a value of greater than or equal to 0.5 classifies the subject as having a high CTSGRS and a value of less than 0.5 classifies the subject as having a low CTSGRS, wherein a high CTSGRS indicates a probability of being above the sex-specific 75thpercentile for CTS and a low CTSGRS indicates a probability of being below the sex-specific 75th percentile for CTS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the sex-specific threshold of total kcal consumed to satiation (CTS) at an ad libitum meal in (b) is the sex-specific 75thpercentile. In some cases, the sex-specific 75thpercentile for males is about 1370 kCal and for females is about 975 kCal. In some cases, the anthropometric data comprises, consists essentially of or consists of height, weight, age and sex. In some cases, the set of genes comprises, consists essentially of or consists of GLP1R, UCP2, FTO, TCF7L2, ADRA2A, LEPR, SH2B1, SIM1, PCSK1, CCK, CCKR, MC4R, MC3R, GIP, GIPR and TNFRSF11A.
[0139] In some cases, the CTSGRS classifier was generated by: (a) receiving in a computer system, a genetic dataset comprising a plurality of single nucleotide polymorphisms (SNPs) located in, around or near a set of genes obtained from samples obtained from each subject from a population of subjects, wherein each subject in the population of subjects is obese, wherein the set of genes comprises, consists essentially of or consists of the genes listed in Table 5; (b) generating by the computer system, a contribution factor for each SNP located in, around or near a gene from the set of genes from the genetic dataset for a subject from the population of subjects, wherein the contribution factor is a numerical value, vector, matrix, or function of various biological factors that represents the contribution or effect that each SNP has on total kcal consumed to satiation (CTS) at an ad libitum meal; (c) calculating by the computer system a gene risk score (GRS) for each gene in the set of genes by integrating the contribution factors generated in step (b) for each SNP for each gene from the set of genes andSNP genotype information located in, around or near each gene in the set of genes, wherein the GRS for each gene represents the cumulative effect of the variation in SNPs located in, around, or near the gene in that particular subject has on the function or regulation of the gene so as to increase or decrease the risk of being above a sex-specific threshold total kcal consumed to satiation (CTS) at an ad libitum meal; (d) iteratively performing by the computer system, steps (a)-(c) for each additional subject from the population of subjects to generate a set of GRSs for each gene in the set of genes, wherein each gene has a GRS for each sex; and (e) integrating by the computer system, the GRSs from step (d) with non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects using a machine learning model to calculate a machine learning assisted gene risk score (ML-GRS) for CTS (CTSGRS), wherein the CTSGRS predicts the CTS at an ad libitum meal for each subject from the population of subjects by combining gene risk scores and the non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects; (f) training by the computer system, the CTSGRS from step (e) using a classifier model to generate the CTSGRS classifier, wherein the CTSGRS classifier generates a value ranging from 0 to 1 for each subject, wherein a value of greater than or equal to 0.5 classifies the subject as having a high CTSGRS and a value of less than 0.5 classifies the subject as having a low CTSGRS, wherein a high CTSGRS indicates a probability of the subject being above the sex-specific 75thpercentile for CTS and a low CTSGRS indicates a probability of the subject being below the sexspecific 75th percentile for CTS. In some cases, the SNP genotype information comprises risk allele dosage, major allele presence, occurrence of a de novo variant, insertion, deletion, or genome rearrangement. In some cases, the subject is obese if the subject possesses a body mass index (BMI) > 30 kg / m2with or without type 2 diabetes. In some cases, the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with a sensitivity and / or specificity is at least 65%. In some cases, the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with an area under the curve (AUC) of at least 0.65, 0.7, 0.75, 0.8, 0.85, 0.9 0.95 or 1. In some cases, the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with positive predictive value of at least 50%. In some cases, the contribution factor is a function or algorithm that utilizes output of a regression analysis or a statistical analysis on the genetic dataset on CTS. In some cases, the genetic dataset is selected from the group consisting of a genome-wide data microarray chip (GWAS), targeted sequencing (exome or targeted genetic panel), variants detected by qPCR through targeted amplification, whole genome sequencing (WGS), targeted or untargeted sequencing of genomic re-arrangement, deletions,duplications, repeat extensions, and complex genotyping of HLA and CYP genes. In some cases, the non-genetic data is anthropometric data selected from the group consisting of age, weight, height, sex and a combination thereof. In some cases, the machine learning model in step (e) comprises a forward feature selection, backward feature selection, or random feature sampling algorithm and a random forest, GBM, or linear predictor that iteratively selects CTS with a desired training accuracy using the random forest predictor in multiple rounds until no further training accuracy improvement is detected. In some cases, the GRSs of step (e) are normalized. In some cases, each GSR provides an indication of a role or effect that the SNPs located around or near a specific gene may play on the specific trait known to be associated with CTS at an ad libitum meal. In some cases, the SNP is considered to be located around or near a specific gene if it is located within at least 500,000 kilobases (kb) of the specific gene. In some cases, the method further comprises performing a genetic analysis on a sample obtained from the subject suffering from or suspected of suffering from obesity prior to step (a). In some cases, the genetic analysis comprises obtaining sequence reads from the whole or portions of the whole genome of the subject. In some cases, the genetic analysis comprises performing a genotyping method selected from the group consisting of restriction fragment length polymorphism identification (RFLPI), random amplified polymorphic detection (RAPD), amplified fragment length polymorphism detection (AFLPD), polymerase chain reaction (PCR), DNA sequencing, RNA sequencing, allele specific oligonucleotide (ASO) probes, and hybridization to microarrays or beads. In some cases, the method further comprises performing by the computer system, a regression analysis and / or statistical analysis on the genetic dataset from step (a) prior to step (b). In some cases, the regression analysis is a ridge regression or least absolute shrinkage and selection operator (LASSO) regression. In some cases, beta values obtained from the regression analyses performed on the genetic datasets are used to generate the contribution factor in step (b). In some cases, the statistical analysis is a genome-wide association study (GWAS). In some cases, p-values and / or regression weights obtained from the statistical analyses performed on the genetic datasets are used to generate the contribution factors in step (b).
[0140] In one embodiment provided herein is a method of predicting weight gain or obesity in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a high CTSGRS indicates that the subject is above the sexspecific 75thpercentile for CTS and is predicted to gain more weight or develop obesity over time as compared to a subject with a low CTSGRS. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) ofless than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0141] In another embodiment, provided herein is a method for treating weight gain or obesity in a subject in need thereof, the method comprising: (a) determining if a subject possesses a high calorie to satiation (CTS) using any of the methods provided herein on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and (b) administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS. In some cases, the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes. In some cases, the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate. In some cases, the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide and tirzepatide.
[0142] In another embodiment, provided herein is a method of predicting upper gastrointestinal (GI) symptoms in a subject. In some cases, the method of predicting upper GI symptoms comprises determining the CTSGRS of a subject using any of the method provided herein. If the subject has a high CTSGRS (i.e., hungry gut positive (HG- or ML-GRS>0.5)) then that subject is predicted to experience upper GI symptoms as compared to a subject with a low CTSGRS (i.e., hungry gut positive (HG+ or ML-GRS<0.5)). In some cases, the upper GI symptom is nausea. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes. In some cases, the subject has not been previously diagnosed with obesity.
[0143] In one embodiment, provided herein is a method of predicting response to bariatric surgery in a subject, the method comprising determining the CTSGRS of a subject using any of the methods provided herein, wherein a high CTSGRS indicates that the subject is abovethe sex-specific 75thpercentile for CTS and is predicted to respond to bariatric surgery as compared to a subject with a low CTSGRS. In some cases, the response to bariatric surgery is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4 or 5 years following bariatric surgery. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0144] In one embodiment, provided herein is a method of predicting response to endoscopic sleeve gastroplasty in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a low CTSGRS indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to endoscopic sleeve gastroplasty as compared to a subject with a high CTSGRS. In some cases, the response to endoscopic sleeve gastroplasty is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time. In some cases, the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following endoscopic sleeve gastroplasty. In some cases, the subject has not been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
[0145] In another embodiment, provided herein is a method of predicting one or more adverse effects (AEs) in a subject following administration of an anti-obesity medicine (AOM). In some cases, the AOM is administered alone. In some cases, the AOM was administered with one or more additional therapeutic agents. In some cases, the method of predicting one or more adverse effects (AEs) in a subject following administration of an anti-obesity medicine (AOM) alone or in combination with one or more additional therapeutic agents comprises determining the CTSGRS of a subject using any of the method provided herein. If the subject has a high CTSGRS (i.e., hungry gut positive (HG- or ML-GRS>0.5)) then that subject is predicted to experience adverse effects when administered a GLP-1 receptor agonist as compared to a subject with a low CTSGRS (i.e., hungry gut positive (HG+ or ML-GRS<0.5)). The one or more AEs is selected from the group consisting of diarrhea, constipation, abdominal discomfort, nausea, vomiting, headache and any combination thereof. In some cases, the AE is nausea. Insome cases, the AE is nausea and headache. In some cases, the subject has been previously diagnosed with obesity. In some cases, the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes. In some cases, subjects being administered a GLP-1RA who are found to a possess a high CTSGRS alter their use of said GLP- 1RA. The alteration can be permanent or temporary discontinuation. The alteration can be to alter the frequency of administration. For example, subjects being administered a GLP-1RA who are found to possess a high CTSGRS may alter their use from daily administration to lower frequency. The lower frequency may be weekly or bi-weekly administration. In some cases, the AOM is a GLP-1 receptor agonist. The GLP-1 receptor agonist can be selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SALO 112, GS-4571, ASC- 30, TTP273, MLX-7006 and BMF-650. In some cases, the one or more additional therapeutic agents can be an antihyperglycemic medications (OAMs). The 0AM can be insulin. In some cases, the one or more additional therapeutic agents can be a sodium-glucose cotransporter-2 (SGLT-2) inhibitor. The SGLT-2 inhibitor can be selected from the group consisting of empagliflozin, canagliflozin, dapagliflozin, ertugliflozin, bexagliflozin and sotagliflozin. In some cases, the one or more additional therapeutic agents is a dipeptidyl peptidase-4 (DPP -4) inhibitor. The DPP -4 inhibitor can be selected from the group consisting of sitagliptin, linagliptin, saxagliptin and alogliptin.
[0146] In one embodiment, the obesity phenotype of interest is the slow bum obesity phenotype. Further to this embodiment, the specific trait in step (b) is a decreased or low resting energy expenditure (REE) for the subject as compared to a control (e.g., individual without obesity. The REE can be measured using any method known in the art such as, for example, indirect calorimetry and can be represented as kcal / 24hrs. The set of genes in step (a) comprises or consists of genes known to play role in conferring the slow burn phenotype. The set of genes in step (a) comprises or consists of the genes listed in Table 5. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MB0AT4 (GOAT), ADRA2A, GNB3, GPBAR1, LEP, LEPR, SH2B1, SIM1, NCOA1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO1, FANCL, CADM2, PTBP2, NUDT3, OLFM4, H0XB5, NPY1R, NPY2R, NPY4R, NPY5R, NR1H4, SLC6A4, NTS, UCP3, ADIPOQ, CCK, CCKR,DPP4, FGFR4, GHRL, MC4R, MC3R, POMC, AGRP, GIP, GIP , CNR1, FA AH, GCG, CELA2A, PPARG, TFAP2B, APOE, TNFRSF11A or DYRK1B. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MB0AT4 (GOAT), ADRA2A, GNB3, LEPR, SH2B1, SIM1, NC0A1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LING01, FANCL, CADM2, PTBP2, NUDT3, 0LFM4, H0XB5, NPY1R, NPY5R, NR1H4, UCP3, FGFR4, GHRL, MC4R, POMC, AG P, GIPR or TFAP2B. In some cases, steps (a)-(d) produce a defined set of GRSs across a plurality of genes. In one embodiment, steps (a)-(d) produce GRSs across a plurality of genes that optimally maximize the true positive rate while minimizing the false positive rate for predicting the slow burn phenotype. In one embodiment, steps (a)-(d) produce GRSs across a plurality of genes that are predictive of the slow burn phenotype with the desired sensitivity, specificity, AUC and / or PPV as provided herein.
[0147] In one embodiment, the obesity phenotype of interest is the emotional eating obesity phenotype. Further to this embodiment, the specific trait in step (b) is a finding or result from a behavioral questionnaire indicative of anxiety or emotional eating for the subject. The questionnaire can be any questionnaire known in the art for assessing anxiety and / or emotional eating, for example, the Hospital Anxiety and Depression Scale (HADS) questionnaire. The set of genes in step (a) comprises or consists of genes known to play a role in conferring the emotional eating phenotype. The set of genes in step (a) comprises or consists of the genes listed in Table 5. In some cases, the plurality of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MBOAT4 (GOAT), ADRA2A, GNB3, GPBAR1, LEP, LEPR, SH2B1, SIM1, NCOA1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO 1, FANCL, CADM2, PTBP2, NUDT3, OLFM4, HOXB5, NPY1R, NPY2R, NPY4R, NPY5R, NR1H4, SLC6A4, NTS, UCP3, ADIPOQ, CCK, CCKR, DPP4, FGFR4, GHRL, MC4R, MC3R, POMC, AGRP, GIP, GIPR, CNR1, FAAH, GCG, CELA2A, PPARG, TFAP2B, APOE, TNFRSF1 1 A or DYRK1B. In some cases, the set of genes can comprise or consist of two or more genes selected from GLP1R, UCP2, FTO, TCF7L2, MBOAT4 (GOAT), ADRA2A, GNB3, LEPR, SH2B1, SIM1, NCOA1 (SRC1), PCSK1, TMEM18, NEGRI, BDNF, GPRC5B, GNPDA2, MTCH2, KCTD15, SEC16B, FAIM2, TNNI3K, LINGO 1, FANCL, CADM2, PTBP2, NUDT3, OLFM4, HOXB5, NPY1R, NPY5R, NR1H4, UCP3, FGFR4, GHRL, MC4R, POMC, AGRP, GIPR or TFAP2B. In some cases, steps (a)-(d) produce a defined set of GRSs across a plurality of genes. In one embodiment, steps (a)-(d) produce GRSsacross a plurality of genes that optimally maximize the true positive rate while minimizing the false positive rate for predicting the emotional eating phenotype. In one embodiment, steps (a)- (d) produce GRSs across a plurality of genes that are predictive of the emotional eating phenotype with the desired sensitivity, specificity, AUC and / or PPV as provided herein.
[0148] It is to be noted that the phenotype, Hungry Brain, as used herein can be defined based on a threshold applied to separate a range of CTS values into categories above the 75%ile. However, other thresholds and logics can define additional phenotypes. One example is the Slow Bum phenotype, which can be defined as the subjects with a resting energy expenditure < the 25%ile. The threshold and logic can be further altered to express phenotypes, such as < 10%ile for an ‘extreme’ phenotype, or between a range, such as 40%ile to 60%ile of CTS being a ‘borderline hungry gut’ . A %ile may be optimized to separate distinct phenotypes, which could result in a %ile threshold from 1-99%. The threshold may or may not be sexspecific, but could be further specified, such as having a different CTS threshold for subjects of different genetic ancestries or having the presence or absence of another designation, such as having uncontrolled type II diabetes. The phenotype can also be defined by a categorical designation of a subject in the absence of any ranged values to which a threshold is applied, such as the presence of coronary artery disease or congestive heart failure. Finally, the phenotype can be defined by a combination of one or more ranged values and designations, such as an ‘emotional hunger-slow burn’ using %ile thresholds being applied to both HADS score and resting energy expenditure, or an ‘early onset hungry brain’ wherein a subject has > 50%ile CTS and history of obesity prior to age 13. The ML-GRS can be tailored to any of these approaches to specifying a phenotype at every stage of development: gene selection, trait & contribution factor selection, feature selection, and ML model training.
[0149] In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a sensitivity of about, at most, at least 65%, 66%, 67%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95% 96%, 97%, 98% or 99%. In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a sensitivity of at least about 65%, at least about 66%, at least about 67%, at least about 68%, at least about 69%, at least about 70%, at least about 71%, at least about 72%, at least about 73%, at least about 74%, at least about 75%, at least about 76%, at least about 77%, at least about 78%, at least about 79%, at least about 80%, at least about 81%, at least about 82%, at least about 83%, at least about 84%, at least about 85%, at least about 86%, at least about 87%, at least about 88%, at least about 89%, at least about 90%, at least about 91%, atleast about 92%, at least about 93%, at least about 94%, at least about 95%, at least about 96%, at least about 97%, at least about 98%, at least about 99%, up to 100%, and all values in between.
[0150] In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a specificity of about, at most, at least 65%, 66%, 67%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95% 96%, 97%, 98% or 99%. In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a specificity of at least about 65%, at least about 66%, at least about 67%, at least about 68%, at least about 69%, at least about 70%, at least about 71%, at least about 72%, at least about 73%, at least about 74%, at least about 75%, at least about 76%, at least about 77%, at least about 78%, at least about 79%, at least about 80%, at least about 81%, at least about 82%, at least about 83%, at least about 84%, at least about 85%, at least about 86%, at least about 87%, at least about 88%, at least about 89%, at least about 90%, at least about 91%, at least about 92%, at least about 93%, at least about 94%, at least about 95%, at least about 96%, at least about 97%, at least about 98%, at least about 99%, up to 100%, and all values in between.
[0151] In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with positive predictive value of about, at most, at least 50%, 51%, 52%, 53%, 54%, 55%, 56%, 57%, 58%, 59%, 60%, 61%, 62%, 63%, 64%, 65%, 66%, 67%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95% 96%, 97%, 98% or 99%. In some cases, the ML-GRS for the obesity phenotype of interest predicts the obesity phenotype of interest with a PPV of at least about 50%, at least about 51%, at least about 52%, at least about 53%, at least about 54%, at least about 55%, at least about 56%, at least about 57%, at least about 58%, at least about 59%, at least about 60%, at least about 61%, at least about 62%, at least about 63%, at least about 64%, at least about 65%, at least about 66%, at least about 67%, at least about 68%, at least about 69%, at least about 70%, at least about 71%, at least about 72%, at least about 73%, at least about 74%, at least about 75%, at least about 76%, at least about 77%, at least about 78%, at least about 79%, at least about 80%, at least about 81%, at least about 82%, at least about 83%, at least about 84%, at least about 85%, at least about 86%, at least about 87%, at least about 88%, at least about 89%, at least about 90%, at least about 91%, at least about 92%, at least about 93%, at least about 94%, at least about95%, at least about 96%, at least about 97%, at least about 98%, at least about 99%, up to 100%, and all values in between.
[0152] A SNP can be considered to be located around or near a specific gene if it is located within at least, at most, about or exactly 1, 10, 50, 100, 150, 200, 250, 300, 350, 400, 450, 500, 600, 700, 800, 900, 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 15000, 20000, 30000, 40000, 50000, 60000, 70000, 80000, 90000, 100000, 200000, 300000, 400000, 500000, 600000, 700000, 800000, 900000 or 1000000 kilobases (KBs) of the specific gene. A SNP is considered to be located around or near a specific gene if it is located within about Ikb to about lOOkb, about lOOkb to about 250kb, about 250 kb to about 500kb, about 500kb to about lOOOkb, about lOOOkb to about 2000kb, about 2000kb to about 5000kb, about 5000kb to about lOOOOkb, about lOOOOkb to about 20000kb, about 20000kb to about 50000kb, about 50000kb to about lOOOOOkb, about lOOOOOkb to about 500000kb or about 500000kb to about 1000000 kbs of the specific gene.Machine Learning, Full Phenotype Test Model
[0153] In yet another aspect, provided herein is a method, which can be computer implemented, for determining an obesity phenotype of a subject suffering from obesity, the method comprising: (a) determining the presence, absence or level of a plurality of gastrointestinal (GI) peptides, a plurality of metabolites, and / or a plurality of genetic variants in a sample obtained from a subject suffering from obesity, thereby generating an obesity analyte signature for the sample; (b) determining satiety, satiation, resting energy expenditure and results on a behavioral questionnaire for the subject; (c) populating a predictive machine learning model with the obesity analyte signature of step (a) and results from step (b) for the subject; and (d) utilizing the predictive machine learning model to predict an obesity phenotype of the subject suffering from obesity. In some cases, the predictive machine learning model is a multivariate logistic regression machine learning model. The multivariate logistic regression can be a least absolute shrinkage and selection operator (LASSO) regression. In some cases, the obesity phenotype is selected from the group consisting of abnormal satiation (hungry brain), abnormal satiety (hungry gut); hedonic eating (emotional hunger) and slow metabolism (slow bum). In some cases, utilization of the predictive machine learning model predicts the obesity phenotype of the mammal suffering from obesity with a sensitivity of about, at most, at least 65%, 66%, 67%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95% 96%, 97%, 98% or 99%. In some cases, utilization of the predictive machine learning model predicts theobesity phenotype of the mammal suffering from obesity with a specificity of about, at most, at least 65%, 66%, 67%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95% 96%, 97%, 98% or 99%. In some cases, utilization of the predictive machine learning model predicts the obesity phenotype of the mammal suffering from obesity with a positive predictive value of about, at most, at least 50%, 51%, 52%, 53%, 54%, 55%, 56%, 57%, 58%, 59%, 60%, 61%, 62%, 63%, 64%, 65%, 66%, 67%, 70%, 71%, 72%, 73%, 74%, 75%, 76%, 77%, 78%, 79%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95% 96%, 97%, 98% or 99%. Satiety can be determined by measuring gastric emptying. Gastric emptying can be measured using any method known in the art such as, for example, scintigraphy and can be represented as GE Tl / 4 or GE Tl / 2. Satiation can be determined by measuring calories to fullness or the amount of calories (e.g., kcal) consumed at an ad libitum meal. Resting energy expenditure can be measured using any method known in the art such as, for example, indirect calorimetry and can be represented as kcal / 24hrs. The behavioral questionnaire can be any questionnaire known in the art for assessing anxiety and / or emotional eating, for example, the Hospital Anxiety and Depression Scale (HADS) questionnaire.Multi-Omic Data
[0154] As described previously herein, any of the machine learning based methods provided herein for determining an individual’s obesity phenotype can utilize or integrate multi-omic data. The multi-omic data utilized or integrated into a machine learning obesity phenotyping method provided herein can comprise genetic or non-genetic multi-omic data obtained from an individual or a sample obtained from an individual as required by said machine learning obesity phenotyping method provided herein. The multi-omic data can be selected from the group consisting of metabolomic data, genomic data, proteomic data, peptidomic data and any combination thereof. In some cases, the multi-omic data obtained for an individual or a sample obtained from an individual represents an obesity analyte signature for that individual.
[0155] An obesity analyte signature can include any appropriate analyte. Examples of analytes that can be included in an obesity analyte signature described herein include, without limitation, DNA, RNA, proteins, peptides, metabolites, hormones, and exogenous compounds (e.g., medications). In some cases, the obesity analyte signature can be obtained by detecting the presence, absence, or level of one or more metabolites, detecting the presence, or absence, or level one or more peptides (e.g., gastrointestinal peptides), and / or detecting the presence orabsence of one or more single nucleotide polymorphisms (SNPs). An obesity analyte signature can be evaluated using any appropriate methods. For example, metabolomics, genomics, microbiome, proteomic, peptidomics, and behavioral questionnaires can be used to evaluate and / or identify an obesity analyte signature described herein.
[0156] A metabolite can be any metabolite that is associated with obesity. In some cases, a metabolite can be an amino-compound. In some cases, a metabolite can be a neurotransmitter. In some cases, a metabolite can be a fatty acid (e.g., a short chain fatty acid). In some cases, a metabolite can be an amino compound. In some cases, a metabolite can be a bile acid. Examples of metabolites that can be used to determine the obesity analyte signature in a sample (e.g., in a sample obtained from an obese mammal) include, without limitation, 1- methylhi stine, serotonin, glutamine, gamma-amino-n-butyric-acid, isocaproic, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, phenylalanine .gamma.-aminobutyric acid, acetic, histidine, LCA, ghrelin, ADRA2A, cholesterol, glucose, acetylcholine, propionic, CDCA, PYY, ADRA2C, insulin, adenosine, isobutyric, 1- methylhistidine, DCA, CCK, GNB3, glucagon, aspartate, butyric, 3 -methylhistidine, UDCA, GLP-1, FTO, leptin, dopamine, valeric, asparagine, HDCA, GLP-2, MC4R, adiponectin, D- serine, isovaleric, phosphoethanolamine, CA, glucagon, TCF7L2, glutamate, hexanoic, arginine, GLCA, oxyntomodulin, 5-HTTLPR, glycine, octanoic, carnosine, GCDCA, neurotensin, HTR2C, myristic, taurine, GDCA, FGF, UCP2, norepinephrine, palmitic, anserine, GUDCA, GIP, UCP3, serotonin, palmitoleic, serine, GHDCA, OXM, GPBAR1, taurine, palmitelaidic, glutamine, GCA, FGF19, NR1H4, stearic, ethanolamine, TLCA, FGF21, FGFR4, oleic, glycine, TCDCA, LDL, elaidic, aspartic acid, TDCA, insulin, GLP-1, linoleic, sarcosine, TUDCA, glucagon, CCK, a-linolenic, proline, THDCA, amylin, arachidonic, alpha-aminoadipic-acid, TCA, pancreatic polypeptide, eicosapentaenoic, DHCA, neurotensin, docosahexaenoic, alpha-amino-N-butyric-acid, THCA, ornithine, GLP-1 receptor, triglycerides, cystathionine 1, GOAT, cystine, DPP4, lysine, methionine, valine, isoleucine, leucine, homocystine, tryptophan, citrulline, glutamic acid, beta-alanine, threonine, hydroxylysine 1, acetone, and acetoacetic acid. In some cases, an obesity analyte signature can include 1-methylhi stine, serotonin, glutamine, gamma-amino-n-butyric-acid, isocaproic, allo- isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, and phenylalanine.
[0157] A gastrointestinal peptide can be any gastrointestinal peptide that is associated with obesity. In some cases, a gastrointestinal peptide can be a peptide hormone. In some cases, a gastrointestinal peptide can be released from gastrointestinal cells in response to feeding.Examples of gastrointestinal peptides that can be used to determine the obesity analyte signature in a sample (e.g., in a sample obtained from an obese mammal) include, without limitation, ghrelin, peptide tyrosine tyrosine (PYY), cholecystokinin (CCK), glucagon-like peptide- 1 (GLP-1), GLP-2, glucagon, oxyntomodulin, neurotensin, fibroblast growth factor (FGF), GIP, OXM, FGF19, and pancreatic polypeptide.
[0158] A SNP can be any SNP that is associated with obesity. A SNP can be in a coding sequence (e.g., in a gene) or a non-coding sequence. For example, in cases where a SNP is in a coding sequence, the coding sequence can be any appropriate coding sequence. Examples of coding sequences that a SNP associated with obesity can be in or near include, without limitation, ADRA2A, ADRA2C, BDNF, CADM2, GNB3, FTO, MC4R, TCF7L2, 5-HTTLPR, HTR2C, UCP2, UCP3, GPBAR1, NR1H4, FAIM2, FANCL, FGFR4, PYY, GLP-1, CCK, leptin, adiponectin, neurotensin, ghrelin, GLP-1 receptor, GOAT, GPRC5B, GNPDA2, DPP4, POMC, NPY, AGRP, SERT, BDNF, SLC6A4, DRD2, LEP, LEPR, UCP1, KLF14, NPC1, LYPLAL1, ADRB2, ADRB3, BBS1, ACSL6, ADARB2, ADCY8, ADH1B, AJAP1, ATP2C2, ATP6V0D2, C21orf7, CAMKMT, CAP2, CASC4, CD48, CDC42SE2, CDYL, CES5AP1, CLMN, CNPY4, COL19A1, COL27A1, COL4A3, COROIC, CPZ, CTIF, DAAM2, DCHS2, DOCKS, EGFLAM, FAM125B, FAM71E2, FRMD3, GALNTL4, GLT1D1, HHAT, HOXB5, KCTD15, KRT23, LEP, LEPR,LHPP, LINC00578, LINC00620, LINGO1, LIPC, LOC100128714, LOC100287160, LOC100289473,LOC100293612|LINC00620, LOC100506869, LOC100507053, LOC100507053|ADH1A, LOCI 005070531 ADH, LOC100507443, LOCI 00996571|CYYR1, LOC152225, LOC255130, LPAR1, LUZP2, MCM7, MICAL3, MMS19, MTCH2, MYBPC1, NEGRI, NR2F2-AS1, NSMCE2, NTN1, NUDT3, O3FAR1, OAZ2, OLFM4, OSBP2, P4HA2, PADI1, PARD3B, PARK2, PCDH15, PCSK1, PIEZO2, PKIB, PRH1-PRR4, PTBP2, PTPRD, RALGPS1|ANGPTL2, RPS24P10, RTN4RL1, RYR2, SCN2A, SEC16B, SEMA3C, SEMA5A, SFMBT2, SGCG, SH2B1, SIM1, SLC22A15, SLC2A2, SLCO1B1, SMOC2 ,SNCAIP, SNX18, NCOA1 (SRC1), SRRM4, SUSD1, TBC1D16, TCERG1L, TENM3, TJP3, TLL1, TMEM9B, TMEM18, TNNI3K, TPM1, VTI1A, VWF, WWOX, WWTR1, ZFYVE28, ZNF3, ZNF609, and ZSCAN21. Examples of SNPS that can be used to determine the obesity analyte signature in a sample (e.g., in a sample obtained from an obese mammal) include, without limitation, rs657452, rsl 1583200, rs2820292, rsl 1126666, rsl 1688816, rsl528435, rs7599312, rs6804842, rs2365389, rs3849570, rsl6851483, rsl7001654, rsl 1727676, rs2033529, rs9400239, rsl3191362, rsl 167827, rs2245368, rs2033732, rs4740619, rs6477694, rsl928295, rsl0733682, rs7899106, rsl7094222, rsl 1191560, rs7903146, rs2176598,rsl2286929, rsl 1057405, rsl0132280, rsl2885454, rs3736485, rs758747, rs2650492, rs9925964, rsl000940, rsl808579, rs7243357, rsl7724992, rs977747, rsl460676, rsl7203016, rsl3201877, rsl441264, rs7164727, rs2080454, rs9914578, rs2836754, rs492400, rsl6907751, rs9374842, rs9641123, rs9540493, rs4787491, rs6465468, rs7239883, rs3101336, rsl2566985, rsl2401738, rsl 1165643, rsl7024393, rs543874, rsl3021737, rsl0182181, rsl016287, rs2121279, rsl3078960, rsl516725, rsl0938397, rsl3107325, rs2112347, rs205262, rs2207139, rsl7405819, rsl0968576, rs4256980, rsl 1030104, rs3817334, rs7138803, rsl2016871, rsl2429545, rsl 1847697, rs7141420, rsl6951275, rsl2446632, rs3888190, rsl558902, rsl2940622, rs6567160, rs29941, rs2075650, rs2287019, rs3810291, rs7715256, rs2176040, rs6091540, rsl800544, Ins-Del-322 , rs5443, rsl 129649, rsl047776, rs9939609, rsl7782313, rs7903146, rs4795541, rs3813929, rs518147, rsl414334, rs659366 , -3474, rs2075577, rsl5763, rsl626521, rsl 1554825, rs4764980, rs434434, rs351855, and rs2234888. Examples of SNPS that can be used to determine the obesity analyte signature in a sample (e.g., in a sample obtained from an obese mammal) include, without limitation, can include any or all of the SNPs associated with an obesity analyte signature as described in WO2022246203 Al, which is herein incorporated by reference in its entirety. For example, SNPS that can be used to determine the obesity analyte signature in a sample (e.g., in a sample obtained from an obese mammal) include, without limitation, rsl664232, rsl 1118997, rs9342434, rs2335852, rsl 1020655, rsl885034, rs7277175, rs6923761, rs7903146, rsl7782313, rs3813929, rsl047776 and any combination thereof.
[0159] An obesity analyte signature described herein can include any appropriate combination of analytes. For example, when an obesity analyte signature includes 14 analytes, the analytes can include 1-methylhi stine, serotonin, glutamine, gamma-amino-n-butyric-acid, isocaproic, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, phenylalanine, ghrelin, and PYY. For example, when an obesity analyte signature includes 9 analytes, the analytes can include HTR2C, GNB3, FTO, isocaproic, beta- aminoisobutyric-acid, butyric, allo-isoleucine, tryptophan, and glutamine.
[0160] Any appropriate method can be used to detect the presence, absence, or level of an obesity analyte within a sample. For example, mass spectrometry (e.g., triple-stage quadrupole mass spectrometry coupled with ultra-performance liquid chromatography (UPLC)), radioimmunoassays, enzyme-linked immunosorbent assays, sequencing techniques (e.g., PCR-based sequencing techniques), and / or restriction fragment length polymorphism (RFLP) can be used to determine the presence, absence, or level of one or more analytes in a sample.
[0161] In some cases, identifying the obesity phenotype can include obtaining results from one or more questionnaires. A questionnaire can be associated with obesity. In some cases, a questionnaire can be answered at the time of the assessment. In some cases, a questionnaire can be answered prior to the time of assessment. For example, when a questionnaire is answered prior to the time of the assessment, the questionnaire results can be obtained by reviewing patient history (e.g., a medical chart). A questionnaire can be a behavioral questionnaire (e.g., psychological welfare questionnaires, questionnaire for assessing anxiety and / or depression, alcohol use questionnaires, eating behavior questionnaires, body image questionnaires, physical activity level questionnaires, and weight management questionnaires). Examples of questionnaires that can be used to determine the obesity phenotype of a mammal (e.g., an obese mammal) include, without limitation, The Hospital Anxiety and Depression Scale (HADS) questionnaire, The Hospital Anxiety and Depression Inventory questionnaire, The Questionnaire on Eating and Weight Patterns, The Weight Efficacy Life-Style (WEL) Questionnaire, Three-Factor Eating Questionnaire (TFEQ), and The Multidimensional Body-Self Relations Questionnaire. For example, a questionnaire can be a HADS questionnaire.
[0162] In some cases, an obesity analyte signature can include the presence of serotonin, glutamine, isocaproic, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, and PYY. For example, a hungry brain (i.e., abnormal satiation) obesity phenotype can have an obesity analyte signature that includes the presence of serotonin, glutamine, isocaproic, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, and PYY. For example, a hungry brain (i.e., abnormal satiation) obesity phenotype can have an obesity analyte signature that has an absence of (e.g., lacks the presence of) 1-methylhi stine, gamma-amino-n-butyric-acid, phenylalanine, ghrelin, and includes a HADS questionnaire result that does not indicate an anxiety subscale (HADS-A; e.g., includes a HADS-A questionnaire result).
[0163] In some cases, an obesity analyte signature can include the presence of 1- methylhi stine, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, and phenylalanine. For example, a hungry gut (i.e., abnormal satiety) obesity phenotype can have an obesity analyte signature that includes the presence of 1-methylhi stine, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, and phenylalanine. For example, a hungry gut (i.e., abnormalsatiety) obesity phenotype can have an obesity analyte signature that has an absence of (e.g., lacks the presence of) serotonin, glutamine, gamma-amino-n-butyric-acid,isocaproic, hexanoic, tyrosine, ghrelin, PYY, and does not include a HAD S questionnaire result that indicates an anxiety subscale (e.g., does not include a HADS-A questionnaire result).
[0164] In some cases, an obesity analyte signature can include the presence of serotonin and can include a HADS-A questionnaire. For example, an emotional or hedonic eating (e.g., behavioral eating) obesity phenotype can have an obesity analyte signature that includes serotonin and includes a HADS-A questionnaire result. For example, an emotional or hedonic eating (e.g., behavioral eating) obesity phenotype can have an obesity analyte signature that has an absence of (e.g., lacks the presence of) 1-methylhi stine, glutamine, gamma-amino-n- butyric-acid, isocaproic, allo-isoleucine, hydroxyproline, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, phenylalanine, ghrelin, and PYY.
[0165] In some cases, an obesity analyte signature can include the presence of 1- methylhi stine, glutamine, gamma-amino-n-butyric-acid, isocaproic, allo-isoleucine, beta- aminoisobutyric-acid, alanine, hexanoic, tyrosine, phenylalanine, PYY, and includes a HADS- A questionnaire result. For example, a large fasting gastric volume obesity phenotype can have an obesity analyte signature that includes 1-methylhi stine, glutamine, gamma-amino-n-butyric- acid, isocaproic, allo-isoleucine, beta-aminoisobutyric-acid, alanine, hexanoic, tyrosine, phenylalanine, PYY, and includes a HADS-A questionnaire result. For example, a large fasting gastric volume obesity phenotype can have an obesity analyte signature that has an absence of (e.g., lacks the presence of) serotonin, hydroxyproline, and ghrelin.
[0166] In some cases, an obesity analyte signature can include the presence of serotonin, beta-aminoisobutyric-acid, alanine, hexanoic, phenylalanine, and includes a HADS- A questionnaire. For example, an obesity phenotype that is mixed can have an obesity analyte signature that includes the presence of serotonin, beta-aminoisobutyric-acid, alanine, hexanoic, phenylalanine, and includes a HADS-A questionnaire result. For example, an obesity phenotype that is mixed can have an obesity analyte signature that has an absence of (e.g., lacks the presence of) 1-methylhi stine, glutamine, gamma-amino-n-butyric-acid, isocaproic, allo- isoleucine, and hydroxyproline.
[0167] In some cases, an obesity analyte signature can include the presence of 1- methylhi stine, serotonin, glutamine, gamma-amino-n-butyric-acid, isocaproic, allo-isoleucine, alanine, tyrosine, ghrelin, PYY, and includes a HADS-A questionnaire result. For example, a low resting energy expenditure (e.g., slow burn) obesity phenotype can have an obesity analyte signature that includes the presence of 1-methylhi stine, serotonin, glutamine, gamma-amino- n-butyric-acid, isocaproic, allo-isoleucine, alanine, tyrosine, ghrelin, PYY, and includes a HADS-A questionnaire result. For example, a low resting energy expenditure (e.g., slow burn)obesity phenotype can have an obesity analyte signature that has an absence of (e.g., lacks the presence of) hydroxyproline, beta-aminoisobutyric-acid, hexanoic, and phenylalanine.
[0168] In some cases, identifying the obesity phenotype also can include identifying one or more additional variables and / or one or more additional assessments. For example, identifying the obesity phenotype also can include assessing the microbiome of a mammal (e.g., an obese mammal). For example, identifying the obesity phenotype also can include assessing leptin levels. For example, identifying the obesity phenotype also can include assessing the metabolome of a mammal (e.g., an obese mammal). For example, identifying the obesity phenotype also can include assessing the genome of a mammal (e.g., an obese mammal). For example, identifying the obesity phenotype also can include assessing the proteome of a mammal (e.g., an obese mammal). For example, identifying the obesity phenotype also can include assessing the peptidome of a mammal (e.g., an obese mammal).Clinical Uses
[0169] Once the obesity phenotype of the mammal has been identified using any of the methods provided herein, the obesity phenotype can be used to select a treatment option for the mammal. In some cases, determining the obesity phenotype of a subject using the methods provided herein can be used to identify the subject as being responsive or non-responsive to a particular intervention (e.g., pharmacological intervention, bariatric procedure, weight loss device, diet intervention, lifestyle intervention, behavior intervention, and / or microbiome intervention) . For example, use of the M-GRS and FPT methods provided herein can be used to identify that the subject possesses a hungry brain phenotype that can serve to identify or select the subject as a non-responder to treatment with a GLP-1 agonist or a responder to phentermine-topiramate pharmacotherapy. For example, once a mammal is identified as being responsive to one or more interventions (e.g., pharmacological intervention, surgical intervention, bariatric procedure, weight loss device, diet intervention, lifestyle intervention, behavior intervention, and / or microbiome intervention) based, at least in part, on an obesity phenotype, which is based, at least in part, on an obesity analyte signature in the sample or CTSGRS score of the mammal, the mammal can be administered or instructed to self-administer one or more interventions. The GLP-1 agonist can be any agent known in the art whose mechanism of action comprises activating the GLP-1 receptor. GLP-1 agonists for use in a method and / or system provided herein can include, but are not limited to, exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735,ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI- 7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650 or any combination thereof.
[0170] When treating obesity in a subject as described herein, the intervention can be effective to reduce the weight, reduce the waist circumference and / or slow or prevent weight gain of the subject. For example, treatment described herein can be effective to reduce the weight (e.g., the total body weight) of a subject with obesity by at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 21%, 22%, 23%, 24%, 25%, 26%, 27%, 28%, 29%, 30%, 31%, 32%, 33%, 34%, 35%, 36%, 37%, 38%, 39%, 40%, 41%, 42%, 43%, 44% or 45%. Treatment described herein can be effective to reduce the weight (e.g., the total body weight) of a with obesity subject by about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 21%, 22%, 23%, 24%, 25%, 26%, 27%, 28%, 29%, 30%, 31%, 32%, 33%, 34%, 35%, 36%, 37%, 38%, 39%, 40%, 41%, 42%, 43%, 44% or 45%. Treatment described herein can be effective to reduce the weight (e.g., the total body weight) of a subject with obesity by at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 21%, 22%, 23%, 24%, 25%, 26%, 27%, 28%, 29%, 30%, 31%, 32%, 33%, 34%, 35%, 36%, 37%, 38%, 39%, 40%, 41%, 42%, 43%, 44% or 45%. In some cases, the intervention described herein can be effective to reduce the weight (e.g., the total body weight) of a subject with obesity by at least 3%, at least 5%, at least 8%, at least 10%, at least 12%, at least 15%, at least 18%, at least 20%, at least 22%, at least 25%, at least 28%, at least 30%, at least 33%, at least 36%, at least 39%, or at least 40%). For example, the intervention described herein can be effective to reduce the weight (e.g., the total body weight) of a subject with obesity by from about 3% to about 40% (e.g., from about 3% to about 35%, from about 3% to about 30%, from about 3% to about 25%, from about 3% to about 20%, from about 3% to about 15%, from about 3% to about 10%, from about 3% to about 5%, from about 5% to about 40%, from about 10% to about 40%, from about 15% to about 40%, from about 20% to about 40%, from about 25% to about 40%, from about 35% to about 40%, from about 5% to about 35%, from about 10% to about 30%, from about 15% to about 25%, or from about 18% to about 22%). For example, the intervention described herein can be effective to reduce the weight (e.g., the total body weight) of a subject with obesity by from about 3 kg to about 100 kg (e.g., about 5 kg to about 100 kg, about 8 kg to about 100 kg, about 10 kg to about 100 kg, about 15kg to about 100 kg, about 20 kg to about 100 kg, about 30 kg to about 100 kg, about 40 kg to about 100 kg, about 50 kg to about 100 kg, about 60 kg to about 100 kg, about 70 kg to about 100 kg, about 80 kg to about 100 kg, about 90 kg to about 100 kg, about 3 kg to about 90 kg, about 3 kg to about 80 kg, about 3 kg to about 70 kg, about 3 kg to about 60 kg, about 3 kg to about 50 kg, about 3 kg to about 40 kg, about 3 kg to about 30 kg, about 3 kg to about 20 kg, about 3 kg to about 10 kg, about 5 kg to about 90 kg, about 10 kg to about 75 kg, about 15 kg to about 50 kg, about 20 kg to about 40 kg, or about 25 kg to about 30 kg). For example, the intervention described herein can be effective to reduce the waist circumference of a subject with obesity by from about 1 inches to about 10 inches (e.g., about 1 inches to about 9 inches, about 1 inches to about 8 inches, about 1 inches to about 7 inches, about 1 inches to about 6 inches, about 1 inches to about 5 inches, about 1 inches to about 4 inches, about 1 inches to about 3 inches, about 1 inches to about 2 inches, about 2 inches to about 10 inches, about 3 inches to about 10 inches, about 4 inches to about 10 inches, about 5 inches to about 10 inches, about 6 inches to about 10 inches, about 7 inches to about 10 inches, about 8 inches to about 10 inches, about 9 inches to about 10 inches, about 2 inches to about 9 inches, about 3 inches to about 8 inches, about 4 inches to about 7 inches, or about 5 inches to about 7 inches. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a reduction of the waist circumference of the subject with obesity of at least 1, 2, 3, 4, 5, 6, 7, 8, 9 or 10 inches. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a reduction of the waist circumference of the subject with obesity of at most 1, 2, 3, 4, 5, 6, 7, 8, 9 or 10 inches. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a reduction of the waist circumference of the subject with obesity of about 1, 2, 3, 4, 5, 6, 7, 8, 9 or 10 inches.
[0171] For example, any treatment described herein can be effective to result in a total body weight lost (TBWL) of the subject with obesity of at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Any treatment described herein can be effective to result in a total body weight lost (TBWL) of the subject with obesity of at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Any treatment described herein can be effective to result in a total body weight lost (TBWL) of the subject with obesity of about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%,14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Any treatment described herein can be effective to result in a total body weight lost (TBWL) of the subject with obesity of 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Any treatment described herein can be effective to result in a total body weight lost (TBWL) of the subject with obesity of from 1% to 5%, from 5% to 10%, from 10% to 15%, from 15% to 20%, from 20% to 25%, from 25% to 30%, from 30% to 35%, from 35% to 40%, from 40% to 45%, from 45% to 50%, from 50% to 55%, from 55% to 60%, from 60% to 65%, from 65% to 70%, from 70% to 75% or from 75% to 80%. In one embodiment, administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of from 2% to 8%. In one embodiment, administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of at least 3%. In one embodiment, administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of at least 6%.
[0172] Administration of any intervention provided herein to a subject determined to possess an obesity phenotype or CTSGRS score using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of at least 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95,99 or 100 kg. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of at most 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 99 or 100 kg. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 99 or 100 kg. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype using a method described herein can result in a total body weight lost (TBWL) of the subject with obesity of from about 3 kg to about 100 kg, about 5 kg to about100 kg, about 8 kg to about 100 kg, about 10 kg to about 100 kg, about 15 kg to about 100 kg, about 20 kg to about 100 kg, about 30 kg to about 100 kg, about 40 kg to about 100 kg, about50 kg to about 100 kg, about 60 kg to about 100 kg, about 70 kg to about 100 kg, about 80 kg to about 100 kg, about 90 kg to about 100 kg, about 3 kg to about 90 kg, about 3 kg to about 80 kg, about 3 kg to about 70 kg, about 3 kg to about 60 kg, about 3 kg to about 50 kg, about 3 kg to about 40 kg, about 3 kg to about 30 kg, about 3 kg to about 20 kg, about 3 kg to about 10 kg, about 5 kg to about 90 kg, about 10 kg to about 75 kg, about 15 kg to about 50 kg, about 20 kg to about 40 kg, or about 25 kg to about 30 kg.
[0173] Administration of any intervention provided herein to a subject determined to possess an obesity phenotype or CTSGRS score using a method described herein can result in a decrease in the waist circumference of the subject with obesity of at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype or CTSGRS score using a method described herein can result in a decrease in the waist circumference of the subject with obesity of at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype or CTSGRS score using a method described herein can result in a decrease in the waist circumference of the subject with obesity of about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype or CTSGRS score using a method described herein can result in a decrease in the waist circumference of the subject with obesity of 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80%. Administration of any intervention provided herein to a subject determined to possess an obesity phenotype or CTSGRS score using a method described herein can result in a decrease in the waist circumference of the subject with obesity of from 1% to 5%, from 5% to 10%, from 10% to 15%, from 15% to 20%, from 20% to 25%, from 25% to 30%, from 30% to 35%, from 35% to 40%, from 40% to 45%, from 45% to 50%, from 50% to 55%, from 55% to 60%, from 60% to 65%, from 65% to 70%, from 70% to 75% or from 75% to 80%.
[0174] In some cases, the change or alteration in the body weight (e.g., total body weight lost (TBWL)) of a subject with an obesity phenotype or CTSGRS score as described herein, the waist circumference of the subject, the characteristics that define the obesity phenotype of the subject (e.g., CTS, gastric emptying, resting energy expenditure (REE),indication of anxiety on HADS, etc.) or any combination thereof that results from administration of any intervention provided herein to the subject results persists for a defined or specified amount of time. The defined or specified amount of time can be at least 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63, 66, 67, 70, 73, 76, 79, 82, 85, 88, 91, 94, 97, 99, 102, 105, 108, 111, 114, 117 or 120 months. The defined or specified amount of time can be at most 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63, 66, 67, 70, 73, 76, 79, 82, 85, 88, 91, 94, 97, 99, 102, 105, 108, 111, 114, 117 or 120 months. The defined or specified amount of time can be about 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63, 66, 67, 70, 73, 76, 79, 82, 85, 88, 91, 94, 97, 99, 102, 105, 108, 111, 114, 117 or 120 months. In some cases, the defined or specified amount of time can be at least 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65 or 70 years.
[0175] Individualized pharmacological interventions for the treatment of obesity (e.g., based on the obesity phenotypes assessed using a method described herein) can include any one or more (e.g., 1, 2, 3, 4, 5, 6, or more) pharmacotherapies (e.g., individualized pharmacotherapies). A pharmacotherapy can include any appropriate pharmacotherapy. In some cases, a pharmacotherapy can be an obesity pharmacotherapy. In some cases, a pharmacotherapy can be an appetite suppressant. In some cases, a pharmacotherapy can be an anticonvulsant. In some cases, a pharmacotherapy can be a GLP-1 agonist. The GLP-1 agonist can be any agent known in the art whose mechanism of action comprises activating the GLP- 1 receptor. GLP-1 agonists for use in a method and / or system provided herein can include, but are not limited to, exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SALO 112, GS-4571, ASC- 30, TTP273, MLX-7006 and BMF-650 or any combination thereof. In some cases, a pharmacotherapy can be an antidepressant. In some cases, a pharmacotherapy can be an opioid antagonist. In some cases, a pharmacotherapy can be a controlled release pharmacotherapy. For example, a controlled release pharmacotherapy can be an extended release (ER) and / or a slow release (SR) pharmacotherapy. In some cases, a pharmacotherapy can be a lipase inhibitor. In some cases, a pharmacotherapy can be a DPP4 inhibitor. In some cases, a pharmacotherapy can be a SGLT2 inhibitor. In some cases, a pharmacotherapy can be a dietary supplement. Examples of pharmacotherapies that can be used in an individualizedpharmacological intervention as described herein include, without limitation, orlistat, phentermine, topiramate, lorcaserin, naltrexone, bupropion, liraglutide, exenatide, metformin, pramlitide, Januvia, canagliflozin, dexamphetamines, prebiotics, probiotics, postbiotics, Ginkgo biloba, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT- 388, VK-2735, ecnoglutide, utresglutide, mazdutide, amy cretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT- 075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650 and combinations thereof. For example, combination pharmacological interventions for the treatment of obesity (e.g., based on the obesity phenotypes as described herein) can include phentermine-topiramate ER, naltrexone-bupropion SR, phentermine-lorcaserin, lorcaserin-liraglutide, and lorcarserin-januvia. In some cases, a pharmacotherapy can be administered as described elsewhere (see, e.g., Sjostrom et al., 1998 Lancet 352: 167-72; Hollander et al., 1998 Diabetes Care 21 : 1288-94; Davidson et al., 1999 IAMA 281 :235-42; Gadde et al., 2011 Lancet 377: 1341-52; Smith et al., 2010 New Engl. I. Med. 363:245-256; Apovian et al., 2013 Obesity 21 :935-43; Pi-Sunyer et al., 2015 New Engl. I. Med. 373: 11-22; and Acosta et al., 2015 Clin Gastroenterol Hepatol. 13:2312-9).
[0176] When a mammal is identified as having an obesity phenotype or CTSGRS score using a machine learning method provided herein (e.g., ML-driven, full phenotype method or adaptive ML-GRS method) that is responsive to treatment with one or more pharmacotherapies, the mammal can be administered or instructed to self-administer one or more pharmacotherapies. In some cases, when a mammal is identified as having a hungry gut (e.g., abnormal satiety) phenotype, based, at least in part, on an obesity analyte signature, the mammal can be administered or instructed to self-administer one or more GLP-1 agonists (e.g., liraglutide or semaglutide) to treat the obesity. In some cases, when a mammal is identified as having a hungry brain (e.g., abnormal satiation) phenotype, based, at least in part, on a hungry brain ML-GRS provided herein, the mammal can be administered or instructed to self- administer phentermine-topiramate pharmacotherapy and / or lorcaserin pharmacotherapy to treat the obesity. In some cases, when a mammal is identified as having a hungry brain (e.g., abnormal satiation) phenotype, based, at least in part, on a hungry brain ML-GRS provided herein, the mammal can be administered or instructed to self-administer a pharmacotherapy that is not a GLP-1 agonist. In other words, a subject who is diagnosed or determined to possess a hungry brain (i.e., abnormal satiation) phenotype can be classified or selected as a nonresponder to GLP-1 agonist treatment. In some cases, when a mammal is identified as havinga slow burn (e.g., energy expenditure) phenotype, based, at least in part, on a slow bum ML- GRS provided herein, the mammal can be administered or instructed to self-administer phentermine to treat the obesity. In some cases, when a mammal is identified as having a hedonic eating (e.g., behavioral / psychological eating) phenotype, based, at least in part, on a hedonic eating ML-GRS provided herein, the mammal can be administered or instructed to self-administer naltrexone / bupropion to treat the obesity. The GLP-1 agonist can be any agent known in the art whose mechanism of action comprises activating the GLP-1 receptor. GLP-1 agonists for use in a method and / or system provided herein can include, but are not limited to, exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide, maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT- 075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650 or any combination thereof.
[0177] In some cases, one or more pharmacotherapies described herein can be administered to an obese mammal as a combination therapy with one or more additional agents / therapies used to treat obesity. For example, a combination therapy used to treat an obese mammal (e.g., a human with obesity) can include administering to the mammal one or more pharmacotherapies described herein and one or more obesity treatments such as weightloss surgeries (e.g., gastric bypass surgery, laparoscopic adjustable gastric banding (LAGB), biliopancreatic diversion with duodenal switch, and a gastric sleeve), vagal nerve blockade, endoscopic devices (e.g., intragastric balloons or endoliners, magnets), endoscopic sleeve gastroplasty, and / or gastric or duodenal ablations. For example, a combination therapy used to treat an obese mammal (e.g., a human with obesity) can include administering to the mammal one or more pharmacotherapies described herein and one or more obesity therapies such as exercise modifications (e.g., increased physical activity, change in type of physical activity, timing and frequency of physical activity), dietary modifications (e.g., reduced-calorie diet), behavioral modifications, lifestyle intervention, commercial weight loss programs, wellness programs, and / or wellness devices (e.g. dietary tracking devices and / or physical activity tracking devices). In cases where one or more pharmacotherapies described herein are used in combination with one or more additional agents / therapies used to treat obesity, the one or more additional agents / therapies used to treat obesity can be administered / performed at the same time or independently. For example, the one or more pharmacotherapies described herein canbe administered first, and the one or more additional agents / therapies used to treat obesity can be administered / performed second, or vice versa.
[0178] In some cases, efficacy of any intervention provided herein can be ascertained by a change or alteration in the body weight (e.g., total body weight lost (TBWL)) of the subject as provided herein, waist circumference, the characteristics that define the hungry brain phenotype or any combination thereof. For example, a subject determined to possess a hungry brain phenotype displayed a decrease in the number of calories consumed to reach fullness following administration of any intervention provided herein as compared to the number of calories consumed to reach fullness for the same subject prior to administration of the intervention. The number of calories consumed to reach fullness following administration of any intervention described herein can be at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% of the number of calories consumed to reach fullness for the same subject prior to administration of the intervention. The number of calories consumed to reach fullness following administration of any intervention described herein can be at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% of the number of calories consumed to reach fullness for the same subject prior to administration of the intervention. The number of calories consumed to reach fullness following administration of any intervention described herein can be about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% of the number of calories consumed to reach fullness for the same subject prior to administration of the intervention. The number of calories consumed to reach fullness following administration of any intervention described herein can be from 1% to 5%, from 5% to 10%, from 10% to 15%, from 15% to 20%, from 20% to 25%, from 25% to 30%, from 30% to 35%, from 35% to 40%, from 40% to 45%, from 45% to 50%, from 50% to 55%, from 55% to 60%, from 60% to 65%, from 65% to 70%, from 70% to 75% or from 75% to 80% of the number of calories consumed to reach fullness for the same subject prior to administration of the intervention. The number of kcal consumed to reach fullness following administration of any intervention described herein can be reduced by at least 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900 or 100 kcal consumed to reach fullness for the same subject prior to administration of the intervention. The number of kcal consumed to reach fullness following administration of any intervention described herein can be reducedby at most 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900 or 100 kcal consumed to reach fullness for the same subject prior to administration of the intervention. The number of kcal consumed to reach fullness following administration of the intervention described herein can be reduced by about 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900 or 100 kcal consumed to reach fullness for the same subject prior to administration of the intervention. The fullness both pre- and -post intervention administration, the fullness can be a maximal fullness (MTV) or normal or usual fullness (VTF). The 'maximal' fullness (MTV) or 'usual' fullness (VTF) can be as measured in a nutrient drink test or to mixed meal (solids) in an ad libitum buffet meal.
[0179] In some cases, efficacy of any intervention provided herein can be ascertained by a change or alteration in the body weight (e.g., total body weight lost (TBWL)) of the subject as provided herein, waist circumference, the characteristics that define the hungry gut phenotype or any combination thereof. For example, a subject determined to possess a hungry gut phenotype displayed a decrease in gastric emptying (e.g., gastric half-emptying time (GE ti / 2)) following administration of any intervention provided herein as compared to the gastric emptying (e.g., gastric half-emptying time (GE ti / 2)) for the same subj ect prior to administration of the intervention. Gastric emptying in a subject determined to possess a hungry gut phenotype following administration of any intervention described herein can be decreased by at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the gastric emptying for the same subject prior to administration of the intervention. Gastric emptying in a subject determined to possess a hungry gut phenotype following administration of any intervention described herein can be decreased by at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the gastric emptying for the same subject prior to administration of the intervention. Gastric emptying in a subject determined to possess a hungry gut phenotype following administration of any intervention described herein can be decreased by about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the gastric emptying for the same subj ect prior to administration of the intervention. Gastric emptying in a subject determined to possess a hungry gut phenotype following administration of any intervention described herein can be decreased from 1% to 5%, from 5%to 10%, from 10% to 15%, from 15% to 20%, from 20% to 25%, from 25% to 30%, from 30% to 35%, from 35% to 40%, from 40% to 45%, from 45% to 50%, from 50% to 55%, from 55% to 60%, from 60% to 65%, from 65% to 70%, from 70% to 75% or from 75% to 80% versus the gastric emptying for the same subject prior to administration of the intervention.
[0180] In some cases, efficacy of any intervention provided herein can be ascertained by a change or alteration in the body weight (e.g., total body weight lost (TBWL)) of the subject as provided herein, waist circumference, the characteristics that define the emotional hunger phenotype or any combination thereof. For example, a subject determined to possess an emotional hunger phenotype displayed a change in the hospital anxiety and depression score (HADS) questionnaire following administration of any intervention provided herein as compared to the same subject prior to administration of the intervention such that the HADS questionnaire result does not indicate a positive an anxiety component following administration of the intervention. In some cases, the change in the hospital anxiety and depression score (HADS) questionnaire following administration of the intervention can be a decrease in the score on the HADS questionnaire by about, at most or at least 1, 2, 3 or 4 as compared to the same subject prior to administration of the intervention such that the HADS questionnaire result.
[0181] In some cases, efficacy of any intervention provided herein can be ascertained by a change or alteration in the body weight (e.g., total body weight lost (TBWL)) of the subject as provided herein, waist circumference, the characteristics that define the slow bum phenotype (e.g., lean muscle mass and / or REE) or any combination thereof. For example, a subject determined to possess a slow burn phenotype displayed an increase in resting energy expenditure (REE); indirect calorimetry; either % or kcal / day) following administration of any intervention provided herein as compared to the REE (either % or kcal / day) for the same subject prior to administration of the intervention. REE in a subject determined to possess a slow burn phenotype following administration of any intervention described herein can be increased by at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the REE for the same subject prior to administration of the intervention. REE in a subject determined to possess a slow burn phenotype following administration of any intervention described herein can be increased by at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the REE for the same subject prior to administration of the intervention. REE in a subject determined to possess a slow bumphenotype following administration of any intervention described herein can be increased by about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the REE for the same subject prior to administration of the intervention. REE in a subject determined to possess a slow bum phenotype following administration of any intervention described herein can be increased from 1% to 5%, from 5% to 10%, from 10% to 15%, from 15% to 20%, from 20% to 25%, from 25% to 30%, from 30% to 35%, from 35% to 40%, from 40% to 45%, from 45% to 50%, from 50% to 55%, from 55% to 60%, from 60% to 65%, from 65% to 70%, from 70% to 75% or from 75% to 80% versus the REE for the same subject prior to administration of the intervention. In some cases, a subject determined to possess a slow burn phenotype displayed an increase in lean muscle mass following administration of any intervention provided herein as compared to the lean muscle mass for the same subject prior to administration of the intervention. Lean muscle mass in a subject determined to possess a slow burn phenotype following administration of any intervention described herein can be increased by at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the lean muscle mass for the same subject prior to administration of the intervention. Lean muscle mass in a subject determined to possess a slow burn phenotype following administration of any intervention described herein can be increased by at most 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the lean muscle mass for the same subject prior to administration of the intervention. Lean muscle mass in a subject determined to possess a slow bum phenotype following administration of any intervention described herein can be increased by about 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 11%, 12%, 13%, 14%, 15%, 16%, 17%, 18%, 19%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75% or 80% versus the lean muscle mass for the same subject prior to administration of the intervention. Lean muscle mass in a subject determined to possess a slow bum phenotype following administration of any intervention described herein can be increased from 1% to 5%, from 5% to 10%, from 10% to 15%, from 15% to 20%, from 20% to 25%, from 25% to 30%, from 30% to 35%, from 35% to 40%, from 40% to 45%, from 45% to 50%, from 50% to 55%, from 55% to 60%, from 60% to 65%, from 65% to 70%, from 70% to 75% or from 75% to 80% versus the lean muscle mass for the same subject prior to administration of the intervention.Systems
[0182] A computer implemented method described herein may be performed on a single computing machine, a virtual machine, a distributed computing system that includes multiples nodes of computing machines, or any other suitable arrangement of computing devices.
[0183] In one embodiment, any of the methods provided herein for determining the obesity phenotype of a subject suffering from obesity may be perform by a system comprising (a) one or more processors; (b) one or more memories or computer-readable medium operatively coupled to at least one of the one or more processors and having instructions stored thereon that, when executed by at least one of the one or more processors, cause the system to performing a method provided herein for determining the obesity phenotype (i.e., ML, FPT method or adaptive, contribution factor based ML-GRS method); and c) one or more instruments in communication with at least one of the one or more processors, wherein the instruments, upon receipt of instructions sent by the at least one of the one or more processors, perform the method provided herein for determining the obesity phenotype (i.e., ML, FPT method or adaptive, contribution factor based ML-GRS method). The one or more processors can be one or more computers and coupled together can form a. computer system for implementing any method provided herein.
[0184] In some embodiments, the computing system operates as a standalone device or may be connected (e.g., networked) to other machines. In a networked deployment, the machine may operate in the capacity of a server machine or a client machine in a serv er-cl lent network environment, or as a peer machine in a peer-to-peer (or distributed) network environment.
[0185] 1 he structure of the computer system may correspond to any software, hardware, or combined components including but not limited to, a client device, a computing server, and various engines, interfaces, terminals, and machines shown.
[0186] By way of example, the computer system may comprise a personal computer (PC), a tablet PC, a set-top box (STB), a personal digital assistant (PDA), a cellular telephone, a smartphone, a web appliance, a network router, an internet of things (loT) device, a switch or bridge, or any machine capable of executing instructions that specify actions to be taken by that machine. Further, the term “machine” and “computer” may also be taken to include any collection of machines that individually or jointly execute instructions to perform any one or more of the methodologies discussed herein.
[0187] Instructions can be any directions, commands, or orders that may be stored in different forms, such as equipment-readable instructions, programming instructions including source code, and other communication signals and orders. Instructions may be used in a general sense and are not limited to machine-readable codes.
[0188] One and more methods described herein improve the operation speed of the processors used in a system provided herein and reduces the space required for the memory. For example, the machine learning methods described herein can reduce the complexity of the computation of the processors by applying one or more novel techniques that simplify the steps in training, reaching convergence, and generating results of the processors. The algorithms described herein also reduces the size of the models and datasets to reduce the storage space requirement for memory'.
[0189] The performance of certain of the operations may be distributed among the more than one processors, not only residing within a single machine, but deployed across a number of machines. In some example embodiments, the one or more processors or processor- implemented engines may be located in a single geographic location (e.g., within a home environment, an office environment, or a server farm). In other example embodiments, the one or more processors or processor-implemented engines may be distributed across a number of geographic locations. Even though in the specification or the claims may refer some processes to be performed by a processor, this should be construed to include a joint operation of multiple distributed processors.
[0190] The term “computer-readable medium” should be taken to include a single medium or multiple media (e.g., a centralized or distributed database, or associated caches and servers) able to store instructions. The computer-readable medium may include any medium that is capable of storing instructions for execution by the processors and that cause the processors to perform any one or more of the methodologies disclosed herein. The computer- readable medium may include, but not be limited to, data repositories in the form of solid-state memories, optical media, and magnetic media. The computer-readable medium does not include a transitory medium such as a propagating signal or a carrier wave and is, thus non- transitory.EXAMPLES
[0191] The present disclosure is further illustrated by reference to the following Examples. However, it should be noted that these Examples, like the embodiments describedabove, are illustrative and are not to be construed as restricting the scope of the invention in any way.Example 1-Abnormal Satiation as a predictor for GLP-1 receptor agonist treatmentBackground
[0192] Obesity is a disease of energy imbalance, where food intake regulation takes a central role. The food intake cycle can be divided in stages that are hunger, satiation and satiety. Multi-omics biomarkers in obesity can provide insight into its etiology and pathophysiology links with chronic illnesses, which may eventually help approach obesity more efficiently and individually. (Aleksandrova et al. 2020). These omics-based biomarkers also allow for the analysis of multiple factors for prediction at the individual level. These high-throughput profiling technologies have explored the connections between a high body mass index or adiposity and the human gut microbiota (Turnbaugh et al. 2009; Dethlefsen and Reiman 2011; Santacruz et al. 2009), plasma metabolome (Cirulli et al. 2019), and host epigenome (Wahl et al. 2017) and genome (Farooqi and O'Rahilly 2007). Genetic and epigenetic modifications that result in transcriptome, proteome, and metabolome changes may have a role in the pathophysiology of energy imbalance seen in obesity.Objective
[0193] The purpose of this study was to uncover the distinctive markers of abnormal satiation using a multi-omic profile of satiation in patients with obesity.Materials and MethodsStudy population
[0194] The study included baseline characteristics from participants recruited between 2010 and 2021. Participants were recruited at Mayo Clinic in Rochester, MN. Ethical approval for the study was obtained from the Mayo Clinic Institutional Review Board (IRB#: 10- 007083, 17-003449, and 19-007485) and registered on ClinicalTrials.gov (NCT00433641, NCT003374956, and NCT04073394). Characteristics measured at baseline from all trials were recorded. Eligible participants were patients with obesity (BMI > 30 kg / m2); otherwise, healthy with no unstable psychiatric disease, aged 18-75 years. Exclusion criteria included abdominal surgery (other than an appendectomy, laparoscopic cholecystectomy, caesarian section, or tubal ligation); positive history of chronic gastrointestinal diseases or systemic disease that could affect gastrointestinal motility; use of medications that may alter gastrointestinal motility, appetite, or absorption; significant psychiatric dysfunction based upon HospitalAnxiety and Depression Scale (HADS), positive alcoholism screening test (AUDIT-C) and the Questionnaire on Eating and Weight Patterns (binge eating disorders and bulimia). All participants provided written informed consent for study participation and further use of their data, including genetic information.
[0195] The validation cohort included participants from a double-blind, placebo- controlled, parallel-group trial testing liraglutide subcutaneous (SQ) 3mg or placebo for a total treatment period of 16 weeks. Protocol and results have been published elsewhere (Maselli et al. 2022, in press).In-person physiologic testing visit
[0196] All participants were screened to complete inclusion and exclusion criteria and were e-mailed a pre-visit study pack with questionnaires (i.e., HADS, AUDIT-C, and the Questionnaire on Eating and Weight Patterns). Questionnaires were completed at home and returned to the study staff at the phenotype visit. Patients were asked to report with prior 8- hour fasting and abstain from caffeine the day prior. Participants arrived between 6:00 - and 8:00 am at the Clinical Research Trials Unit (CRTU). Heart rate and blood pressure were measured using an automated blood-pressure monitor while participants fasted. Participant weight, height, and hip and waist circumference were measured using standard clinical techniques. Following the anthropometric measurements, participants had a baseline fasting blood draw, consumed a standardized radiolabel breakfast for gastric emptying, and 4 hours later had an ad libitum meal for lunch (FIG. 1). Table 1 shows the baseline discovery cohort characteristics measured in the in-person physiological testing.
[0197] Table 1 Baseline discovery cohort characteristics measured in the in-person physiological testing day.Outcome variables and sample collection, handling, and analysisGastric emptying for solids
[0198] A standard breakfast of 320 kcal and 30% fat consisting of two 99mTc- radiolabeled eggs, toast, and 80 mL of skim milk was given to participants. Images were obtained immediately after radiolabeled meal ingestion and every 30 minutes for the first 2 hours, then every hour for the next 2 hours (total 4 hours after the radiolabeled meal) (Camilleri et al. 2012). Gastric emptying (GE) is summarized by the half-emptying time (GE T1 / 2) in minutes.Appetite sensations
[0199] The visual analog scales (VAS) was measured using a validated, standard, 100mm appetite VAS for hunger, fullness, desire to eat, and satisfaction score (Flint et al. 2000; Calderon et al. 2020). Appetite VAS was assessed 15 minutes before breakfast and, then every 30 minutes for the first 2 hours.Satiation test: Ad libitum meal
[0200] The ad libitum meal included: vegetable lasagna [Stouffers®, Nestle USA, Inc, Solon, OH; nutritional analysis of each 326g box: 420kcal, 17g protein (16% of energy), 38g carbohydrate (37% of energy), and 22g fat (47% of energy)]; vanilla pudding [Hunts®, Kraft Foods North America, Tarrytown, NY; nutritional analysis of each 99g carton: 130kcal, 1gprotein (3% of energy), 21g carbohydrate (65% of energy), and 4.5g fat (32% of energy)]; and skim milk [nutritional analysis of each 236mL carton: 90kcal, 8g protein (36% of energy), 13g carbohydrate (64% of energy), and Og fat]. The total amount (g and kcal) of food consumed and the kcal of each macronutrient at the ad libitum meal were analyzed by a registered dietitian using validated software (ProNutra 3.0; Viocare Technologies Inc, Princeton, NJ).
[0201] Participants were served the prepared meal. Once the subject had partially consumed the meal, they were inquired if additional lasagna, pudding, or milk as desired. Subjects were asked to use any food products until they achieved maximum fullness. As the subject completed each meal component, the empty plate, empty pudding container, and / or liquid volume of milk were weighed back. Once the subject reached maximum fullness, the total weight of food products consumed was summarized. The primary measurement was total kcal consumed to fullness (CTF).
[0202] Table 2. Baseline characteristics of participants with normal and abnormal satiation.Sample collection and DNA extraction
[0203] Baseline blood samples following an overnight fast were collected from each subject into 6ml EDTA blood collection tubes. Serum was kept at room temperature and immediately prepared for subsequent analysis of routine laboratory analyses or aliquoted andstored at - 80 °C until further analyses. Plasma was obtained from EDTA tubes, immediately placed on ice, and centrifuged within 10 min (1500 g, 4 °C, 10 min). Plasma samples were aliquoted and stored at - 80 °C until further analyses.Metabolomics
[0204] The Mayo Clinic Metabolomics Core conducted a quantitative analysis of amino acids and amino acid metabolites in plasma samples using a liquid chromatographytandem mass spectrometry (LC-MS / MS) method with Waters MassTrack Amino Acid Analysis Solution (Waters, Milford, MA). (Lanza et al. 2010) Plasma samples and amino acid reference standards were prepared using Waters' MassTrak Amino Acid (AA) Analysis Solution kit, with minor modifications for detection on a mass spectrometer. A 10-point standard concentration curve was made from the calibration standard solution to calculate amino acid concentrations in plasma samples. A total of 46 amino acids and metabolites were able to be measured. Table 3 lists metabolites measured in this study.
[0205] Table 3 List of metabolites measured in this study.Genotyping
[0206] According to the manufacturer's instructions, DNA was extracted from whole blood as previously described (Acosta, Camilleri, Shin, Vazquez-Roque, et al. 2015) byTaqMan® SNP Genotyping Assays (Applied Biosystems, Foster City, CA) using 10-20ng DNA. Following PCR amplification, end reactions were read on the ABI 7500 Fast Real-Time PCR system using Sequence Detection Software version 1.3.1 (Applied Biosystems). SNPs previously reported to be associated with variation of BMI, postprandial satiation, and with pharmacological response to sibutramine or involved in the control of gastrointestinal functions in a preliminary study were extracted from the full set of genome-wide genotypes using PLINK and were tested for association with gastric emptying using linear regression methods. Selected SNPs are described in Table 4.
[0207] Table 4 Minor allele frequencies of selected genes and reference for selection in this study.Predictive Model Creation and Validation
[0208] The abnormal satiation phenotype was then determined for patients by creating a machine-learning (ML) model for appetite using the full phenotype tests (FPTs). To generatethe ML predictive model, two independent cohorts of patients with obesity were included who each completed the following validated tests at baseline: a) satiation, studied by ad libitum buffet meal (kcal consumed to maximal fullness), and b) gastric emptying of solids studied by scintigraphy. A fasting blood sample was collected for satiety (PYY, GLP1) hormones (ELISA), metabolomics (Mass Spectometry), and DNA (SNP Array). In the first cohort of a total of 799 patients (mean BMI of 38.5 ± 4.38 kg / m2), a multivariate logistic regression ML- driven model for appetite was developed by using the FPTs with all baseline data collected. The linear predictor was calculated with the model's parameter estimates to fit y=a+bx. The ML, FPT model was then validated in an independent cohort of 141 patients (BMI: 39.09 ± 7.40 kg / m2).Statistical Analysis
[0209] Categorical variables are presented as percentages, and all continuous data are summarized as mean and standard deviation (SD). Descriptive statistics were used to show demographic, anthropometric, and physiological parameters. Statistical testing of distributions (Abnormal postprandial satiety vs. normal postprandial satiety) were assessed by a two sampled unpaired t-test if data followed a Gaussian distribution or Mann Whitney test if nonGaussian distribution. All tests were two-tailed, and a p-value <0.05 was considered statistically significant. Multivariable regressions were applied to investigate the association between satiation as continues variable and our exposures (demographics, anthropometries, vital signs, questionnaires, metabolites, and hormones), and R2values were reported. Statistical analysis was performed using JMP®, Version 14.3.0 (SAS Institute Inc., Cary, NC, 1989-2019) and visualized with GraphPad Prism (Version 8).ResultsCohort Measurements and Heterogeneity of Satiation
[0210] To evaluate food intake, and particularly satiation, adults were recruited with obesity (body mass index >30 mg / kg2) aged between 21 to 70 years to complete a full day of in-person physiologic testing (FIG. 1). 717 subjects were studied during a full day of in-person physiological testing. Participants were mostly middle-aged females (75%) and White Americans, with a mean BMI of 36.2 kg / m2(Table 1). The findings were validated with baseline physiological characteristics of 163 adults with obesity recruited for a weight loss trial, methods and results of the trial have been previously published (Maselli et al. 2022, in press).
[0211] All participants reported after an 8-hour fasting and had a standardized breakfast of 320kcal with 30% to measure gastric emptying for 4 hours and complete appetite ratings in this period. The mean GE T ’A was 127.0 ± 32.7 min, females have a slower gastric emptyingwhen compared to males (132.7 ± 31.5 min vs. 109.6 ± 30.3 min; p<0.001) (FIG. 2A). With the same standardized meal, females reported lower levels of hunger and higher levels of fullness and satisfaction than males (FIGs 2B-D). Following their last scan, participants had their satiation test which consisted of an ad libitum meal. Participants were served the prepared meal and were asked to eat until reaching fullness. The mean calories to fullness were 917.4 ± 315.5 kcal, and females consume fewer calories to reach fullness (835.2 ± 259.1 kcal vs. 1164.2 ± 340.5 kcal; p<0.001) (FIG. 2E-H). These highlights the huge heterogeneity in food intake in patients with obesity and the gender differences which required to assess these components separately between females and males.Predicting Individual Satiation within Patients with Obesity.
[0212] Participants were asked to complete behavioral questionnaires and were assessed by a physician to confirm their medical history. A targeted metabolite panel (Table 3) was measured along with enteroendocrine hormones (i.e., glucagon like peptide 1 [GLP-1], peptide YY [PYY], and ghrelin). All data at different -omics level was used to explain the variability of satiation. The data showed that among the input variables baseline characteristics (age, gender, and race) can explain the most significant proportion of postprandial satiety variance (r2= 0.24) (FIG. 3).Characterization of Abnormal Satiation
[0213] Next, the cohort was divided into quartiles by gender and participants with normal and abnormal satiation were compared. In the normal satiation group, females consumed less than 650 kcal to reach fullness and males less than 927 kcal. While in the abnormal satiation group, females consumed more than 977 kcal to reach fullness and males more than 1374 kcal (Table 2). First, gastric emptying between groups was evaluated, participants with abnormal satiation have a faster gastric emptying with 111.7 ± 29.8 min half emptying time compared with 124.5 ± 32.7 min in the normal satiation group (FIG. 4A). The hormonal curves between groups was then evaluated. There was no difference in ghrelin, but participants with abnormal satiation were observed have higher levels of GLP-1 (11.98 ± 8.5 ug / ml vs. 9.36 ± 6.9 ug / ml) and PYY (11.98 ± 8.5 ug / ml vs. 9.36 ± 6.9 ug / ml) after 90 minutes of a standardized meal of 320 kcal (30% fat) (FIGs. 4B-D). Finally, appetite sensations between groups were compared, and participants with abnormal satiation had more hunger and less satisfaction at all time points after their breakfast, lower fullness at baseline and at 60, 90, and 120 minutes (FIGs. 4E-G). And in general, participants with abnormal satiation wanted toeat more than their control group at all timepoints. There are large interpersonal differences in satiation, but this classification allowed characterization of this group of participants who present a different subjective response measured with appetite scales and objective response measured with gastric emptying and postprandial hormonal levels to the same meal.Prediction and Validation of Normal and Abnormal Satiation
[0214] Next, clinical and -omics factors were integrated into an algorithm that differentiate normal vs. abnormal satiation (i.e., the machine-learning (ML) model for appetite using the full phenotype tests (FPTs). To this end, the algorithm was developed on the main cohort of 799 participants, and performance was evaluated in the validation cohort of 141 subjects that completed the same assessments. Integrating demographics, anthropometries, and fasting blood parameters, the ML-driven, full phenotype biomarker had an AUC of >0.80, while the validation cohort had an AUC of 0.84 (FIG. 5).Validation of Normal and Abnormal Satiation discrimination for prediction of response to liraglutide
[0215] To assess the clinical relevance of the predictor model, we tested the predictions with short term response to liraglutide were tested using the ML-driven, full phenotype biomarker. Previously, baseline characteristics have not been able to discriminate responders from non-responders. Energy intake has been used as a clinical measure to assess the effect of liraglutide on weight loss(Van Can et al. 2014; Horowitz et al. 2012; Flint, Kapitza, and Zdravkovic 2013). And weight loss was found to be associated with a decrease in kcal consumed during ad libitum meal in two studies(Van Can et al. 2014; Flint, Kapitza, and Zdravkovic 2013). Responders were classified based on a clinically significant response of greater than 5% total body weight lost and the ML-driven, full phenotype biomarker model was used to discriminate responders and non-responders. In particular, the performance of the FPT (n=60) was tested in predicting response to treatment with liraglutide (LIRA) 3mg SQ / day from a randomized, placebo-controlled trial of 16 weeks. A response was defined as >5% total body weight loss [TBWL%], Analysis was done using student’s T-test to compare among the predicted groups. In participants with liraglutide, the ML-driven, full phenotype biomarker model had an AUC of 0.70 to identify responders (FIG. 6A). To confirm its utility to detect response to liraglutide rather than just weight loss, the model was tested in patients with placebo and the model had an AUC of 0.63 (CI 0.53 to 0.73) (FIG. 6B). Finally, the actual total body weight loss percentage was evaluated in participants predicted as responders (i.e., patients with a non-hungry brain or non-abnormal satiation phenotype) and non-responders (i.e., patients with a hungry brain or abnormal satiation phenotype) and the TBWL% in groups predicted asresponder vs. non-responder using FPT was -7.6+3.4 vs. -3.6+3.3; p<0.001 (FIG. 6C). The same testing was performed in the placebo group, and here, the difference was non-significant between groups (FIG. 6D).Discussion
[0216] Individual aspects of satiation were explored at different omic levels and a prediction model for abnormal satiation was developed. Basic variables, such as demographics or anthropometries, cannot explain observed variability and are insufficient in predicting abnormal postprandial satiety, highlighting the importance of deep phenotyping incorporating many variables in characterizing physiological components of food intake. The ML, FPT prediction model developed in this Example had an AU00.80 and was validated in an independent cohort of patients with obesity. And this allowed discrimination of responders to GLP-1 analogues in a short-term placebo-controlled trial.
[0217] Individual variability in predisposition to weight gain and resistance to weight reduction may partly be explained by inter-individual variability in all stages of food intake, hunger, satiation, and satiety. (Gibbons et al. 2019) Abnormal satiation is thought to be a risk factor for overconsumption (Drapeau et al. 2007) and is associated with obesity. (Acosta, Camilleri, Shin, Vazquez -Roque, et al. 2015) .Example 2- Development of Obesity Phenotype Polygenic Risk Score / Machine Learning PredictorBackground
[0218] Human genetics drive our propensity to have certain characteristics. These characteristics can often be complex and result from differences in several parts of our genome. Now that technologies are available to identify genetic differences between individuals (variants), the challenge becomes identifying which variants are relevant to a particular characteristic, such as susceptibility to a disease. Being able to identify relevant variants and how they contribute to a characteristic enables prediction of, among other things, disease risk. One method known in that art that can do this identification and can predict the likelihood of a characteristic are Polygenic Risk Scores (PRSs).
[0219] A PRS typically starts with a basic regression or statistical model to identify the relevance of a SNP. Using a regression model, the relevance is the effect size of a SNP (i.e., beta value of the regression). Using a statistical model, the relevance is whether the p-value of association between the SNP and the particular trait is less that a certain threshold. Whichever the approach, the information about the variants in a particular sample are combined. Forregression models, it may be the sum of the beta values for SNPs with an effect allele, while for statistical models, the count of the number of SNPs below the p-value threshold with an effect allele. Additional transformations can be involved, such as various forms of normalization, as well as numerous quality control steps, such as accounting for linkage disequilibrium. However, these approaches have several shortcomings.
[0220] One basic shortcoming is that identifying which variants are relevant can be ambiguous due to the structure of genetic inheritance. Further, the low statistical power of the basic genetic analyses underlying PRSs can create PRSs that are irrelevant to the characteristic. Because the number of samples required to have sufficient statistical power to overcome these sources of ambiguity is infeasible, a different approach was warranted.Objective
[0221] The purpose of this study was to develop a new type of PRS that integrates expert guidance, machine learning, and non-genetic / multi-omic data that addresses the shortcomings of standard genetic analysis and can serve as a simpler, DNA only, surrogate to the ML, FTP obesity phenotype predictor developed in Example 1. The approach taken in this Example side stepped resolving ambiguity by allowing it to be present in initial steps of generating a PRS and in later steps used expert and data-derived guidance to minimize ambiguity. Further, while the approach used herein uses some standard PRS components, such as quality control steps to control from sample relatedness and regression between genetics of a set of samples and a particular trait, the approach used herein avoided the irrelevance of naively scoring a particular characteristic directly and instead sought to first estimate the contribution of an individual’s genetics to a set of traits associated with the particular characteristic.Materials and Methods
[0222] Conceptually, the adaptive contribution-factor-based PRS of this Example was developed as follows:
[0223] Rather than using regression beta values across the genome naively, the beta in the regression model was transformed using an algorithm that encoded biological factors about genetics, epigenetics, gene regulation, and similar, for each variant. The output numerical value representing the contribution of a variant to a trait after this algorithm runs was termed the Contribution Factor. This was step (3) in the diagram shown in FIG. 7.
[0224] Rather than using SNPs independently, the effect variants at SNPs may have on a particular gene were summarized, creating an intermediary Gene Risk Score (GRS). To generate a GRS, the dosage of variants around a gene were integrated with the contributionfactor of each SNP. The resulting value was a risk score that represented the likelihood that the variants around that gene change the function or regulation of that gene, which in turn plays a role in the characteristic.
[0225] Because the gene risk scores reflected a particular trait, and not the characteristic to be predicted, extra steps were taken to identify gene risk scores that best explained the target characteristic. Expert guidance and machine learning techniques were combined to iteratively remove irrelevant genes (rather than variants individually). This iterative process ended when the gene list was whittled down to the shortest list of genes whose gene risk scores for a set of samples best explained the characteristic in those samples. This was the Adaptive steps (4), (5), and (6) in the diagram shown in FIG. 7.
[0226] Finally, the PRS combined these intermediary gene risk scores for all selected traits as well as non-genetic and / or multi-omic data (step (7) in FIG. 7) to predict the target characteristic. This prediction was the adaptive contribution-factor-based polygenic risk score (step (8) in FIG. 7).
[0227] In this Example, an adaptive contribution-factor-based PRS was developed to predict the Hungry Brain (i.e., abnormal satiation) obesity phenotype. Hungry Brain is a type of obesity driven by a defect in satiation: the brain does not properly signal that the stomach is full, leading to overconsumption of calories. Hungry Brain is defined as consuming above the 75%ile amount of calories at an ad libitum meal. SNP variant information from quality- controlled samples from 497 of the 799-patient cohort (mean BMI of 38.5 ± 4.38 kg / m2) referenced in Example, was obtained using a SNP chip. Quantitative traits were also collected for these 497 patients. A regression between the variants in these samples and a trait called ‘buffet meal’ was performed. This ‘buffet meal’ trait was the number of calories consumed at an ad libitum meal. Because this trait was substantially different in males vs females, separate regressions for male samples and female samples were performed rather than including them as a cofactor. In particular, a LASSO regression through the PLINK 2.0 tool using default settings was used to perform the regression analyses.
[0228] Subsequently, contribution factors and gene risk scores were created, and the adaptive step of the process outlined in FIG. 7 was performed. The contribution factor was derived from the beta value but was the output of a function that also takes other parameters into account. The other parameters included the proximity of the variant to the gene, apparent mechanistic role of the variant (e.g., as within an intron vs exon vs enhancer region), or metadata about the SNP in relationship to the gene, such as the likelihood that it alters expression of the gene. The contribution factor and dosage information of SNPs within eachgene and wi thing 500KB of that gene were then integrated to create a gene risk score for each gene. Importantly, the GRS was not calculated for all known genes, but only for an expert- guided set of 22 genes. In fact, for each of these genes, two GRSs were created, one using contribution factors based on the male regression model and the other based on the female regression model. Using forward feature selection and random forest predictors, a plurality of GRSs across a plurality of the genes were identified that are predictive of the Hungry Brain phenotype in the training cohort, with the remaining GRSs not providing any further benefit.
[0229] Finally, a PRS using the selected GRSs was created and validated in an independent cohort of 141 patients (BMI: 39.09 4^7.40 kg / m2). The final PRS was built using the 13 selected GRSs normalized with a Yeo power transformation and the height of the individual using a support vector classifier (SVC). The SVC was trained using a linear kernel and a C value of 0.001 using SMOTE oversampling of the 497 samples. The final PRS was then the trained classifier’s predicted probability of Hungry Brain (e.g., abnormal satiation).
[0230] Subsequently, the performance of the PRS (n=40) model in predicting response to treatment with liraglutide (LIRA) 3mg SQ / day from a randomized, placebo-controlled trial of 16 weeks was tested. Like in Example 1, a response was defined as >5% total body weight loss [TBWL%], Analysis was done using student’s T-test to compare among the predicted groups.Results and Conclusions
[0231] In the training cohort the PRS model had an AUC>0.80. In the independent validation cohort, the PRS had an AUC=0.70 with a sensitivity of 0.71, specificity of 0.65, and PPV of 0.53 was measured (see FIGs 8A-8B). Note that because the phenotype was defined as being above the 75%ile of calories consumed, a priori a PPV of 0.25 for random guessing. When used to predict response to liraglutide, the PRS an AUC=0.69. In patients treated with liraglutide, the TBWL% in groups predicted as responder vs. non-responder using the PRS was -8.0+4.4 vs. -5.6+3.3; p=0.056 (FIGs 9A-9B). There was no difference among the predicted groups in patients who received placebo.
[0232] The ML-driven, full phenotype biomarker and its PRS surrogate for appetite can predict response to treatment with a GLP-1 analogue. The findings support a novel approach for precision medicine for obesity. The PRS model might be a simpler test to deploy in the clinical practice.References
[0233] Acosta, Andres, Michael Camilleri, Andrea Shin, Maria I. Vazquez -Roque, Johanna Iturrino, Duane Burton, Jessica O'Neill, Deborah Eckert, and Alan R. Zinsmeister. 2015. 'Quantitative gastrointestinal and psychological traits associated with obesity and response to weight-loss therapy', Gastroenterology, 148: 537-46. e4.
[0234] Acosta, Andres, Michael Camilleri, Andrea Shin, Maria I Vazquez-Roque, Johanna Iturrino, Ian R Lanza, K Sreekumaran Nair, Duane Burton, Jessica O'Neill, and Deborah Eckert. 2015. 'Association of UCP-3 rs 1626521 with obesity and stomach functions in humans', Obesity, 23: 898-906.
[0235] Aleksandrova, K., C. Egea Rodrigues, A. Floegel, and W. Ahrens. 2020. 'Omics Biomarkers in Obesity: Novel Etiological Insights and Targets for Precision Prevention', Curr Obes Rep, 9: 219-30.
[0236] Aronne, Louis J, Kevin D Hall, John M. Jakicic, Rudolph L Leibel, Michael R Lowe, Michael Rosenbaum, and Samuel Klein. 2021. 'Describing the weight-reduced state: Physiology, behavior, and interventions', Obesity, 29: S9-S24.
[0237] Blundell, JE, and JCG Halford. 1994. 'Regulation of nutrient supply: the brain and appetite control', Proceedings of the Nutrition Society, 53: 407-18.
[0238] Calderon, G., A. McRae, J. Rievaj, J. Davis, I. Zandvakili, S. Linker-Nord, D. Burton, G. Roberts, F. Reimann, B. Gedulin, A. Vella, N. F. LaRusso, M. Camilleri, F. M. Gribble, and A. Acosta. 2020. 'Ileo-colonic delivery of conjugated bile acids improves glucose homeostasis via colonic GLP-1 -producing enteroendocrine cells in human obesity and diabetes', EBioMedicine, 55: 102759.
[0239] Camilleri, M., J. Iturrino, A. Bharucha, D. Burton, A. Shin, I. D. Jeong, and A. Zinsmeister. 2012. 'Performance characteristics of scintigraphic measurement of gastric emptying of solids in healthy participants', Neurogastroenterology and motility: the official journal of the European Gastrointestinal Motility Society, 24: 1076.
[0240] Camilleri, Michael. 2015. 'Peripheral mechanisms in appetite regulation', Gastroenterology, 148: 1219-33.
[0241] Cannon, Walter B, and Anton L Washbum. 1993. 'An Explanation of Hunger 1', Obesity research, 1 : 494-500.
[0242] Christoff ersen, Berit Ostergaard, Guillermo Sanchez-Delgado, Linu Mary John,Donna H Ryan, Kirsten Raun, and Eric Ravussin. 2022. 'Beyond appetite regulation: Targeting energy expenditure, fat oxidation, and lean mass preservation for sustainable weight loss', Obesity, 30: 841-57.
[0243] Choi SW, Mak TS, O'Reilly PF. Tutorial: a guide to performing polygenic risk score analyses. Nat Protoc. 2020 Sep;15(9):2759-2772. doi: 10.1038 / s41596-020-0353-l. Epub 2020 Jul 24. PMID: 32709988; PMCID: PMC7612115.
[0244] Cifuentes, Lizeth, and Andres Acosta. 2022. 'Homeostatic regulation of food intake', Clinics and Research in Hepatology and Gastroenterology, 46: 101794.
[0245] Cirulli, E. T., L. Guo, C. Leon Swisher, N. Shah, L. Huang, L. A. Napier, E. F. Kirkness, T. D. Spector, C. T. Caskey, B. Thorens, J. C. Venter, and A. Telenti. 2019. 'Profound Perturbation of the Metabolome in Obesity Is Associated with Health Risk', Cell Metab, 29: 488-500. e2.
[0246] Collaborators, G. B. D. Obesity, A. Afshin, M. H. Forouzanfar, M. B. Reitsma, P. Sur, K. Estep, A. Lee, L. Marczak, A. H. Mokdad, M. Moradi-Lakeh, M. Naghavi, J. S. Salama, T. Vos, K. H. Abate, C. Abbafati, M. B. Ahmed, Z. Al-Aly, A. Alkerwi, R. Al- Raddadi, A. T. Amare, A. Amberbir, A. K. Amegah, E. Amini, S. M. Amrock, R. M. Anjana,J. Amlov, H. Asayesh, A. Banerjee, A. Barac, E. Baye, D. A. Bennett, A. S. Beyene, S. Biadgilign, S. Biryukov, E. Bjertness, D. J. Boneya, I. Campos-Nonato, J. J. Carrero, P. Cecilio,K. Cercy, L. G. Ciobanu, L. Comaby, S. A. Damtew, L. Dandona, R. Dandona, S. D. Dharmaratne, B. B. Duncan, B. Eshrati, A. Esteghamati, V. L. Feigin, J. C. Fernandes, T. Furst, T. T. Gebrehiwot, A. Gold, P. N. Gona, A. Goto, T. D. Habtewold, K. T. Hadush, N. Hafezi- Nejad, S. I. Hay, M. Horino, F. Islami, R. Kamal, A. Kasaeian, S. V. Katikireddi, A. P. Kengne, C. N. Kesavachandran, Y. S. Khader, Y. H. Khang, J. Khubchandani, D. Kim, Y. J. Kim, Y. Kinfu, S. Kosen, T. Ku, B. K. Defo, G. A. Kumar, H. J. Larson, M. Leinsalu, X. Liang, S. S. Lim, P. Liu, A. D. Lopez, R. Lozano, A. Majeed, R. Malekzadeh, D. C. Malta, M. Mazidi, C. McAlinden, S. T. McGarvey, D. T. Mengistu, G. A. Mensah, G. B. M. Mensink, H. B. Mezgebe, E. M. Mirrakhimov, U. O. Mueller, J. J. Noubiap, C. M. Obermeyer, F. A. Ogbo, M. O. Owolabi, G. C. Patton, F. Pourmalek, M. Qorbani, A. Rafay, R. K. Rai, C. L. Ranabhat, N. Reinig, S. Safiri, J. A. Salomon, J. R. Sanabria, I. S. Santos, B. Sartorius, M. Sawhney, J. Schmidhuber, A. E. Schutte, M. I. Schmidt, S. G. Sepanlou, M. Shamsizadeh, S. Sheikhbahaei, M. J. Shin, R. Shiri, I. Shiue, H. S. Roba, D. A. S. Silva, J. I. Silverberg, J. A. Singh, S. Stranges, S. Swaminathan, R. Tabares-Seisdedos, F. Tadese, B. A. Tedla, B. S. Tegegne, A. S. Terkawi, J. S. Thakur, M. Tonelli, R. Topor-Madry, S. Tyrovolas, K. N. Ukwaja, O. A. Uthman, M. Vaezghasemi, T. Vasankari, V. V. Vlassov, S. E. Vollset, E. Weiderpass, A. Werdecker, J. Wesana, R. Westerman, Y. Yano, N. Yonemoto, G. Yonga, Z. Zaidi, Z. M. Zenebe, B. Zipkin, and C. J. L. Murray. 2017. Health Effects of Overweight and Obesity in 195 Countries over 25 Years', N Engl J Med, 377: 13-27.
[0247] Dethlefsen, Les, and David A Reiman. 2011. 'Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation', Proceedings of the National Academy of Sciences, 108: 4554-61.
[0248] Drapeau, Vicky, Neil King, Marion Hetherington, Eric Doucet, John Blundell, and Angelo Tremblay. 2007. 'Appetite sensations and satiety quotient: predictors of energy intake and weight loss', Appetite, 48: 159-66.
[0249] Farooqi, I. S., and S. O'Rahilly. 2007. 'Genetic factors in human obesity', Obes Rev, 8 Suppl 1 : 37-40.
[0250] Flint, A, C Kapitza, and M Zdravkovic. 2013. 'The once-daily human GLP-1 analogue liraglutide impacts appetite and energy intake in patients with type 2 diabetes after short-term treatment', Diabetes, Obesity and Metabolism, 15: 958-62.
[0251] Flint, A., A. Raben, J. E. Blundell, and A. Astrup. 2000. 'Reproducibility, power and validity of visual analogue scales in assessment of appetite sensations in single test meal studies', Int J Obes Relat Metab Disord, 24: 38-48.
[0252] Galmes, S., A. Palou, and F. Serra. 2020. 'Increased Risk of High Body Fat and Altered Lipid Metabolism Associated to Suboptimal Consumption of Vitamin A Is Modulated by Genetic Variants rs5888 (SCARB1), rsl800629 (UCP1) and rs659366 (UCP2)', Nutrients, 12.
[0253] Gibbons, Catherine, Mark Hopkins, Kristine Beaulieu, Pauline Oustric, and John E Blundell. 2019. 'Issues in measuring and interpreting human appetite (satiety / satiation) and its contribution to obesity', Current Obesity Reports, 8: 77-87.
[0254] Halawi, H., M. Camilleri, A. Acosta, M. Vazquez -Roque, I. Oduyebo, D. Burton, I. Busciglio, and A. R. Zinsmeister. 2017. 'Relationship of gastric emptying or accommodation with satiation, satiety, and postprandial symptoms in health', Am J Physiol Gastrointest Liver Physiol, 313: G442-g47.
[0255] Hill, James O, Holly R Wyatt, and John C Peters. 2012. 'Energy balance and obesity', Circulation, 126: 126-32.
[0256] Horowitz, Michael, Anne Flint, Karen L Jones, Charlotte Hindsberger, Mads F Rasmussen, Christoph Kapitza, Selena Doran, Thomas Jax, Milan Zdravkovic, and Ian M Chapman. 2012. 'Effect of the once-daily human GLP-1 analogue liraglutide on appetite, energy intake, energy expenditure and gastric emptying in type 2 diabetes', Diabetes Research and Clinical Practice, 97: 258-66.
[0257] Lanza, Ian R, Shucha Zhang, Lawrence E Ward, Helen Karakelides, Daniel Raftery, and K Sreekumaran Nair. 2010. 'Quantitative metabolomics by 1 H-NMR and LC- MS / MS confirms altered metabolic pathways in diabetes', PLoS One, 5: el0538.
[0258] Marginean, C., C. O. Marginean, C. Banescu, L. E. Meli], F. Tripon, and M. lancu. 2019. 'The relationship among GNB3 rs5443, PNPLA3 rs738409, GCKRrs780094 gene polymorphisms, type of maternal gestational weight gain and neonatal outcomes (STROBE- compliant article)', Medicine (Baltimore), 98: el6414.
[0259] Maselli, Daniel, Jessica Atieh, Matthew M. Clark, Deborah Eckert, Ann Taylor,Paula Carlson, Duane D. Burton, Irene Busciglio, W. Scott Harmsen, Adrian Vella, Andres Acosta, and Michael Camilleri. 2022, in press. 'Effects of liraglutide on gastrointestinal functions and weight in obesity: a randomized clinical and pharmacogenomic trial', Obesity.
[0260] Miranda, R. C., S. B. Vetter, J. P. Genro, P. D. Campagnolo, V. S. Mattevi, M. R. Vitolo, and S. Almeida. 2015. 'SLC6A14 and 5-HTR2C polymorphisms are associated with food intake and nutritional status in children', Clin Biochem, 48: 1277-82.
[0261] Papathanasopoulos, A., M. Camilleri, P. J. Carlson, A. Vella, S. J. Nord, D. D. Burton, S. T. Odunsi, and A. R. Zinsmeister. 2010. 'A preliminary candidate genotypeintermediate phenotype study of satiation and gastric motor function in obesity', Obesity (Silver Spring), 18: 1201-11.
[0262] Salopuro, T., L. Pulkkinen, J. Lindstrom, M. Kolehmainen, A. M. Tolppanen, J. G. Eriksson, T. T. Valle, S. Aunola, P. Ilanne-Parikka, S. Keinanen-Kiukaanniemi, J. Tuomilehto, M. Laakso, and M. Uusitupa. 2009. 'Variation in the UCP2 and UCP3 genes associates with abdominal obesity and serum lipids: the Finnish Diabetes Prevention Study', BMC Med Genet, 10: 94.
[0263] Santacruz, Arlette, Ascension Marcos, Julia Wamberg, Amelia Marti, Miguel Martin-Matillas, Cristina Campoy, Luis A Moreno, Oscar Veiga, Carlos Redondo-Figuero, and Jesus M Garagorri. 2009. 'Interplay between weight loss and gut microbiota composition ioverweight adolescents', Obesity, 17: 1906-15.
[0264] Tumbaugh, Peter J, Micah Hamady, Tanya Yatsunenko, Brandi L Cantarel, Alexis Duncan, Ruth E Ley, Mitchell L Sogin, William J Jones, Bruce A Roe, and Jason P Affourtit. 2009. 'A core gut microbiome in obese and lean twins', Nature, 457: 480-84.
[0265] van Abeelen, A. F., M. de Krom, J. Hendriks, D. E. Grobbee, R. A. Adan, and Y. T. van der Schouw. 2008. 'Variations in the uncoupling protein-3 gene are associated with specific obesity phenotypes', Eur J Endocrinol, 158: 669-76.
[0266] Van Can, J, B Sloth, CB Jensen, A Flint, EE Blaak, and WHM Saris. 2014. 'Effects of the once-daily GLP-1 analog liraglutide on gastric emptying, glycemic parameters, appetite and energy metabolism in obese, non-diabetic adults', International Journal of Obesity, 38: 784-93.
[0267] Wahl, Simone, Alexander Drong, Benjamin Lehne, Marie Loh, William R Scott, Sonja Kunze, Pei-Chien Tsai, Janina S Ried, Weihua Zhang, and Youwen Yang. 2017. 'Epigenome-wide association study of body mass index, and the adverse outcomes of adiposity', Nature, 541 : 81-86.
[0268] Wrzosek, M., A. Sawicka, M. Wrzosek, P. Pietkiewicz, M. Talalaj, and G. Nowicka. 2019. 'Age at onset of obesity, transcription factor 7-like 2 (TCF7L2) rs7903146 polymorphism, adiponectin levels and the risk of type 2 diabetes in obese patients', Arch Med Sci, 15: 321-29.
[0269] Yu, K., L. Li, L. Zhang, L. Guo, and C. Wang. 2020. 'Association between MC4R rs 17782313 genotype and obesity: A meta-analysis', Gene, 733: 144372.Example 3- A polygenic risk score to identify abnormal satiation predicts response to phentermine-topiramateBackground
[0270] Anti-obesity medications have highly variable outcomes. The inventors proposed an obesity classification based on phenotypes to address this heterogeneity in weight loss response. In a 2-week randomized clinical trial (RCT), the inventors showed that one of the phenotypes, abnormal satiation, measured with calories to fullness (CTF) during an ad libitum buffet meal (>75th percentile) had superior weight loss outcomes with phenterminetopiramate ER (phen-top). However, current methods used to identify phenotypes are timeconsuming, invasive, expensive, and limited to a few academic centers. Here, a biomarker for abnormal satiation has been developed and validated and its efficacy was assessed to predict weight loss outcomes in a RCT of phen-top vs placebo.Methods
[0271] The biomarker for abnormal satiation was developed in an independent cohort of 497 patients using a polygenic risk score (PRS) derived from an expert-guided set of SNPs in 22 genes, and demographic characteristics. This model was validated in a new independent cohort of 141 patients. A post-hoc analysis of a 12-month, double-blinded, single-center RCT was conducted that randomized 80 patients with obesity to either placebo (n=38) or phen-top (n=42). These patients completed phenotype testing at baseline. Wilcoxon tests were used tocompare the total body weight loss (TBWL) among patients with or without abnormal satiation, defined by either CTF in an ad libitum buffet meal or with the biomarker, in patients assigned to phen-top or placebo.Results
[0272] In the independent validation cohort, the biomarker for abnormal satiation had a ROC=0.72. In the trial, at 12-months, the phen-top group had a TBWL of 16.4 kg (10.5), as compared to 6 kg (11.8) (p=0.005) in placebo. At 12 months, within patients assigned to phen- top, those with abnormal satiation defined by CTF (40%) had a higher TBWL than the normal satiation group [23.1 kg (11.9) vs. 12.8 kg (8.0), p=0.02]. The biomarker performance for predicting a TBWL >15% at 12 months in patients assigned to phen-top had an ROC=0.88, and those with a positive test (20%) had a higher TBWL than those with a negative test [30 kg (3.9) vs. 12.8 kg (1.9), p=0.006]. There was no difference in the placebo group with either approach.
[0273] Conclusion: Abnormal satiation phenotype, measured by CTF in ad libitum buffet meal or by a polygenic risk score, predicts response to phentermine-topiramate ER.Example 4- A machine-learning-assisted gene risk score to identify abnormal satiation and it use in predicting response to anti-obesity interventionsObjective
[0274] This example aims to explore the variability of human satiation, a key component of energy intake, and its potential implications for precision obesity management. Through this, satiation was revealed to be highly variable among individuals, and while baseline characteristics, anthropometries, body composition and hormones, contribute to this variability, they alone cannot account for the totality of the variance in satiation. To address this gap, the role of a polygenic risk score was explored, which demonstrated a robust association with satiation. Furthermore, a machine-learning gene risk score was introduced to predict satiation and to assess this score in the prediction of responses to anti-obesity medications.MethodsStudy population
[0275] Baseline characteristics were collected from participants with obesity recruited between 2010 and 2021 at Mayo Clinic, Rochester, MN. Eligible participants were adults aged 18 to 75 years, with obesity (BMI > 30 kg / m2), with or without type 2 diabetes, and otherwise healthy. Exclusion criteria included a history of abdominal surgery (except for appendectomy,laparoscopic cholecystectomy, cesarean section, or tubal ligation), chronic gastrointestinal diseases, systemic diseases affecting gastrointestinal motility, and the use of medications that may alter gastrointestinal motility, appetite, or absorption. Ethical approval for the use of data and samples from prior studies for the characterization and development of a biomarker was obtained from the Mayo Clinic Institutional Review Board (IRB#: 17-000838). All participants provided written informed consent for study participation, including the use of their genetic information in future studies.
[0276] Study Cohorts:
[0277] This study included two distinct cohorts of participants. The initial study cohort comprised a total of 717 participants who completed a full day of in-person deep-phenotyping. This initial cohort was further divided into sub-cohort 1, utilized for training the machinelearning models, and a sub-cohort 2 used for validation of the models. The validation cohort consisted of 62 participants that also were part of the 52-week, double-blind, placebo- controlled, randomized clinical trial comparing phentermine-topiramate extended-release versus placebo (NCT04408586).
[0278] The second cohort consisted of 110 participants who underwent deep- phenotyping and were enrolled in a 16-week, double-blind, placebo-controlled, randomized clinical trial comparing liraglutide versus placebo6(NCT03523273). This cohort served as the independent testing cohort. Upon excluding related samples and samples from outlier genetic backgrounds, the validation cohort for the genetic models was refined to a total of 50 samples.Study Protocol for Deep Phenotyping and sample collection, and handling,In-person physiologic testing visit
[0279] Participants underwent screening to meet inclusion and exclusion criteria and were sent a pre-visit study pack containing questionnaires (Hospital Anxiety and Depression Scale [HADS] questionnaire, Questionnaire on Eating and Weight Patterns, and the Three- Factor Eating Questionnaire [TFEQ-R21]) to complete at home. Upon arriving at the Clinical Research Trials Unit (CRTU) at Mayo Clinic, Rochester, MN, between 6:00 am and 8:00 am, participants were measured for heart rate and blood pressure while fasting (i.e., overnight fast). Standard clinical techniques were used to measure weight, height, hip, and waist circumference. Following the anthropometric measurements, participants underwent a baseline fasting blood draw, consumed a standardized radiolabeled breakfast to measure gastric emptying, and had an ad libitum meal for lunch (FIG. 10A).Outcome variables and sample collection, handling, and analysis
[0280] Gastric emptying for solids
[0281] Participants received a standard breakfast comprising 320 kcal and 30% fat, including two 99mTc-radiolabeled eggs, toast, and 80 mL of skim milk. Images were taken immediately after ingestion and at regular intervals for a total of 4 hours. Gastric emptying (GE) was summarized by the half-emptying time (GE Tl / 2) in minutes.38(Camilleri et al. 2012).
[0282] Fasting and postprandial gastric volumes
[0283] Participants gastric volumes were measured by single photon emission computed tomography imaging (SPECT) of the stomach after intravenous injection of 99mTc- pertechnetate during fasting and after 300mL Ensure drink39.
[0284] Satiation test: Ad libitum meal
[0285] The ad libitum meal included: vegetable lasagna [Stouffers®, Nestle USA, Inc, Solon, OH; nutritional analysis of each 326g box: 420kcal, 17g protein (16% of energy), 38g carbohydrate (37% of energy), and 22g fat (47% of energy)]; vanilla pudding [Hunts®, Kraft Foods North America, Tarrytown, NY; nutritional analysis of each 99g carton: 130kcal, 1g protein (3% of energy), 21g carbohydrate (65% of energy), and 4.5g fat (32% of energy)]; and skim milk [nutritional analysis of each 236mL carton: 90kcal, 8g protein (36% of energy), 13g carbohydrate (64% of energy), and 0g fat]. The total amount (g and kcal) consumed and the kcal of each macronutrient at the ad libitum meal were analyzed by a registered dietitian using validated software (ProNutra 3.0; Viocare Technologies Inc, Princeton, NJ). Participants were served the prepared meal. Subjects were asked to consume any food products until they achieved maximum satiation. A visual analog scale was used to assess maximal satiation (rate 0 to 5). As the subject completed each meal component, the empty plate, empty pudding container, and / or liquid volume of milk were weighed back. Once the subject reached maximum satiation, the total weight of food products consumed was summarized. The primary measurement was total kcal consumed to satiation (CTS).
[0286] Eating behaviors
[0287] The Three-Factor Eating Questionnaire (TFEQ-R21) is a validated questionnaire designed to assess eating and weight control behaviors40. Participants completed the TFEQ-R21 on their baseline assessment visit. The TFEQ-R21 is a 21 -item instrument that measures three domains of eating behavior: cognitive restraint (CR), uncontrolled eating (UE), and emotional eating (EE). The first twenty items are rated on a 4-point Likert scale, and item 21 is answered through an 8-point Likert scale. Before calculating domain scores, items 1-16 were reverse coded, and item 21 was recorded as follows: 1-2 scores as 1; 3-4 as 2; 5-6 as 3; 7-8 as 4. Domain scores were then calculated as a means of all items within each domain;hence, domain scores ranged from 1 to 4 (CR [six items], UE [nine items], and EE [six items]), with higher scores being indicative of greater CR, UE, and EE.
[0288] Sample collection
[0289] Informed consent was obtained from all participants. Baseline blood samples following an overnight fast were collected from each subject into 6ml EDTA blood collection tubes. Serum was kept at room temperature and immediately prepared for subsequent analysis of routine laboratory analyses or aliquoted and stored at - 80 °C until further analyses.
[0290] DNA extraction
[0291] The samples were thawed at room temperature and a lysis buffer was added. The samples were then incubated at 56°C for 10 minutes. Proteinase K was added, and the samples were incubated at 56°C for an additional 30 minutes. Phenol / chloroform was added, and the samples were vortexed to mix. The samples were then centrifuged at 10,000 rpm for 10 minutes. The aqueous layer was transferred to a new tube and ethanol was added. The samples were vortexed to mix and then centrifuged at 10,000 rpm for 10 minutes. The supernatant was removed, and the pellet was air-dried. The pellet was then resuspended in the TE buffer. The DNA was purified using a commercial DNA purification kit. The kit included a lysis buffer, proteinase K, phenol / chloroform, and ethanol. The DNA was purified according to the manufacturer's instructions. The DNA was quantified using a NanoDrop spectrophotometer. The DNA concentration was determined by measuring the absorbance of the DNA at 260 nm.
[0292] Metabolomics
[0293] Blood samples were collected from participants and stored at -80°C until ready for analysis. The samples were thawed at room temperature and vortexed to mix. The samples were analyzed using liquid chromatography-mass spectrometry (LC-MS / MS). The samples were first separated on an LC column using a gradient elution method. The eluted compounds were then detected by mass spectrometry and the peak areas were quantified. The metabolites were identified by their mass-to-charge (m / z) ratios and retention times.
[0294] Hormones
[0295] In this example, the levels of appetite-regulating hormones were measured, Peptide YY (PYY), glucagon -like peptide- 1 (GLP-1), and ghrelin, in human blood samples using enzyme-linked immunosorbent assay (ELISA). ELISA is a highly sensitive and specific technique that allows for the quantification of target hormones in plasma. Monoclonal antibodies specific to PYY, GLP-1, and ghrelin were immobilized onto microtiter plates, and the blood samples were added to allow hormone binding. Subsequently, a second enzyme-linked antibody was added, initiating a color-producing reaction proportional to the hormone concentration.
[0296] Genome-Wide Association Study (GWAS)
[0297] Genotyping was carried out using the Illumina Infinium BeadChip technology. Quality control measures were implemented to ensure data accuracy and reliability, and samples with low call rates or ambiguous results were excluded from further analysis. Prior to analysis, further metrics were implemented to ensure only accurate genotype calls were used. SNPs with a call rate < 0.95 and samples with call rate<0.95 were removed. The relatedness between each sample was assessed using the kinship coefficient. One subject from the following related pairs was retained: twins, parent / offspring, and full siblings. The sample with the higher call rate was retained. Genomic ancestry was assessed using STRUCTURE. NonEuropean samples were removed (e.g., >20% non-European). Calls for SNPs with heterozygosity < lxlOA-6, MAF of 1%, and an information score of < O.X. SNP flips between samples were handled automatically by the analysis software.
[0298] P-values for each SNP were derived using a set of 483 samples. These samples have both genotyping data for 2.7M SNPs typed using Illumina’s Omni Exome array and CTS and exclude related samples and samples from outlier genetic backgrounds. Using PLINK 211, linear regression was performed using age as a covariate. Regression was performed separately for males and females as they have different CTS endpoints for high calories to satiation (FIG.22B)
[0299] Candidate Gene Selection & Polygenic Risk Score
[0300] A set of candidate genes to build a polygenic risk score were selected from a literature review based on their biological relevance to obesity and previous evidence from genetic studies (Table 5). Genes involved in appetite regulation, energy expenditure, lipid metabolism, and adipogenesis were prioritized for inclusion in the study. This resulted in a set of 41 genes and 19,161 unique annotated variants, of which 2,637 were SNPs contained in the Omni Exome array. The weighted risk alleles were correlated with CTS for each SNP that had genetic variability across all samples and the top 10% most correlated SNPS (n = 227) were selected as those most likely involved in appetite regulation. Subsequently, the PRS for each individual was calculated by summing the weighted risk alleles of the chosen SNPs.
[0301] Table 5. List of selected candidate genes for the development of the polygenic risk score for satiation.
[0302] Gene Risk Score Development
[0303] In preparation for training a machine-learning gene risk score model to predict high calories to satiation (CTSGRS), a method for the scoring of genes was created (FIG. 22A- 22F). This score numerically represents the overall contribution of variants in and around a candidate gene to CTS and was done separately for males and females. To calculate this score, contributing factors at variants near the 41 candidate genes were investigated. These factors (deemed “contribution factors”) included the resulting beta value from linear regression analysis regressing SNP dosage on CTS (using sex-specific endpoints) (FIG. 22B) and weighting each SNP based on the proximity to the gene, and its presence within the gene (FIG. 22D). SNPs within one megabase pairs flanking each candidate gene were included (Table 5).This process resulted in two numerical representations for each gene, which were separately derived for males and females based on the beta values from the endpoints of the GWAS. Each numerical representation was termed a Gene Risk Score (GRS). There were two GRSs for each gene, one for males and one for females.
[0304] Table 6. Demographic characteristics and parameters of energy balance regulation by low and high calories consumed to satiation (CTS).Abbreviations used: AOM, anti-obesity medications; BMI, body mass index; CTS, calories to satiation; Hx, history; VAS, visual analog scale. Categorical variables were compared with
[0305] To develop a CTSGRS to predict high CTS, the GRSs were considered in addition to demographic and anthropometric sources (e.g., age, sex, weight) (all of which were considered features). A random forest model with backward feature selection was used to identify the most informative features as measured by training AUC (see FIG. 22E). In brief, starting with a single feature, the feature with the best training accuracy using a random forest model was selected. Then the remaining feature with the best training accuracy when combined with the previously selected feature(s) were selected. This was continued until no accuracy improvement was observed. This method selected height, weight, age, sex, and 16 genetic risk scores associated with 10 genes: SIMl*, PCSK1, SH2B1, LEPR, UCP2*, FTO*, TCF7L2*, GLP1R, TNFRSF11A, and ADRA2A. Genes marked with an asterisk had only the GRS derived from the female endpoint selected (see Table 5 for the number of SNPs included in each GRS). Note that because of the algorithm’s greedy approach, features were selected when they were both predictive and contributed novel information. Thus, some expected genes, such as MC4R, were not selected because functionally related genes already contained similar predictive information, not because such genes were not predictive of CTS.
[0306] Scores for selected features were generated for all samples with both genetic information and CTS: 50 samples in the Phentermine / Topiramate trial; 110 in the Liraglutide trial; and 483 that participated in neither trial.
[0307] Machine-Learning Model Training To Predict High Calories to Satiation
[0308] The CTSGRS was trained in 483 samples using a Support Vector Classifier (SVC) with the 16 selected GRSs, height, weight, age, and sex to predict high calories to satiation positive individuals (HS+, which can also be referred to as high CTSGRS). HS+ was defined as those individuals whose CTS was above the sex-specific 75th %ile (Male = 1372.4275, Female = 977.3325) and the others as HS- (can also be referred to as low CTSGRS). A binary indicator of sex was then included along with the genetic risk scores. Before training, the scores were normalized using a Yeo power transform41and conducted SMOTE42sampling to create 225 synthetic samples for training. An SVC model 10-fold cross-validation on these oversampled data was trained, with each fold’s test data normalized based on the training data from that fold (see FIG. 22F). A mean AUC = 0.82 was measured across the 10 folds. The final SVC model used all the training samples. All SVC models were trained using the Python sklearn.svm package with a linear kernel, C = 3, gamma = 0.01, and a balanced class weight.The final CTSGRS model used oversampling and Yeo power transform normalization for all 483 training samples and 225 synthetic (708 total) samples and yielded both an HS+ / - prediction and a probability of HS+ (p(HS+)).
[0309] Machine-Learning Model Validation
[0310] To evaluate the model, the performance of the final SVC satiation model’s prediction accuracy was tested on two placebo-controlled trials. These samples included 50 from a 52-week phentermine / topiramate trial (Test) which also underwent calories to satiation testing and had genotyping data. These samples’ data were normalized the same as the training (i.e., the Yeo power transformer fit to the 708 training samples) and the HS+ (high CTSGRS) class was set using the same CTS thresholds (75th%ile CTS thresholds). The outcomes of the samples of the treatment arm of the study was also evaluated for those predicted to HS+ (high CTSGRS) and HS- (low CTSGRS) using a threshold of p(HS+) > 0.5 to define HS+. The Total Body Weight Loss of individuals was measured at 12-, 26-, 36-, and 52 weeks in each group under the hypothesis that predicted HS+ individuals are better responders (>15%TBWL at 52 weeks). Comparison of predicted and actual HS+ / - is shown in Table 12A-12B.
[0311] Machine-Learning Model Independent Testing and Validation
[0312] The performance of the final SVC model’s predictions on the Liraglutide trial samples was then me...
Claims
CLAIMSWhat is claimed:
1. A computer-implemented method for predicting a subj ect’ s calories to satiation (CTS), the method comprising:(a) receiving in a computer system, genetic data obtained from a sample obtained from a subject as well as anthropometric data from the subject;(b) calculating by the computer system, a gene risk score (GRS) for each respective gene in a set of genes from the genetic data obtained from the sample obtained from the subject, wherein the GRS for each respective gene in the set of genes represents a cumulative effect of single nucleotide polymorphisms (SNPs) located in, around, or near the respective gene has on the function or regulation of the respective gene in the subject so as to increase or decrease the risk of the subject being above a sex-specific threshold of total kcal consumed to satiation (CTS) at an ad libitum meal, wherein the set of genes comprises, consists essentially of or consists of the set of genes listed in Table 5;(c) normalizing by the computer system, the anthropometric data and the GRS for each respective gene from the set of genes to training samples from a CTSGRS classifier to generate normalized GRSs and anthropometric data for the subject; and(d) inputting by the computer system, the normalized GRSs and anthropometric data for the subject into the CTSGRS classifier to generate a value ranging from 0 to 1, wherein a value of greater than or equal to 0.5 classifies the subject as having a high CTSGRS and a value of less than 0.5 classifies the subject as having a low CTSGRS, wherein a high CTSGRS indicates a probability of being above the sex-specific 75thpercentile for CTS and a low CTSGRS indicates a probability of being below the sexspecific 75th percentile for CTS.
2. The computer-implemented method of claim 1, wherein the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample.
3. The computer-implemented method of claim 1 or 2, wherein the subject has not been previously diagnosed with obesity.
4. The computer-implemented method of claim 1 or 2, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
5. The computer-implemented method of claim 1 or 2, wherein the subject has been previously diagnosed with obesity.
6. The computer-implemented method of claim 1 or 2, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
7. The computer-implemented method of any one of the above claims, wherein the sexspecific threshold of total kcal consumed to satiation (CTS) at an ad libitum meal in (b) is the sex-specific 75thpercentile.
8. The computer-implemented method of any one of the above claims, wherein the sexspecific 75thpercentile for males is about 1370 kCal and for females is about 975 kCal.
9. The computer-implemented method of any one of the above claims, wherein the anthropometric data comprises, consists essentially of or consists of height, weight, age and sex.
10. The computer-implemented method of any one of the above claims, wherein the set of genes comprises, consists essentially of or consists of GLP1R, UCP2, FTO, TCF7L2, ADRA2A, LEPR, SH2B1, SIM1, PCSK1 and TNFRSF11A.
11. The computer-implemented method of any one of the above claims, wherein the CTSGRS classifier was generated by:(a) receiving in a computer system, a genetic dataset comprising a plurality of single nucleotide polymorphisms (SNPs) located in, around or near a set of genes obtained from samples obtained from each subject from a population of subjects, wherein each subject in the population of subjects is obese, wherein the set of genes comprises, consists essentially of or consists of the genes listed in Table 5;(b) generating by the computer system, a contribution factor for each SNP located in, around or near a gene from the set of genes from the genetic dataset for a subject from the population of subjects, wherein the contribution factor is a numerical value, vector, matrix, or function of various biological factors that represents the contribution or effect that each SNP has on total kcal consumed to satiation (CTS) at an ad libitum meal;(c) calculating by the computer system a gene risk score (GRS) for each gene in the set of genes by integrating the contribution factors generated in step (b) for each SNP for each gene from the set of genes and SNP genotype information located in, around or near each gene in the set of genes, wherein the GRS for each gene represents the cumulative effect of the variation in SNPs located in, around, or near the gene in that particular subject has on the function or regulation of the gene so as to increase ordecrease the risk of being above a sex-specific threshold total kcal consumed to satiation (CTS) at an ad libitum meal;(d) iteratively performing by the computer system, steps (a)-(c) for each additional subject from the population of subjects to generate a set of GRSs for each gene in the set of genes, wherein each gene has a GRS for each sex; and(e) integrating by the computer system, the GRSs from step (d) with non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects using a machine learning model to calculate a machine learning assisted gene risk score (ML-GRS) for CTS (CTSGRS), wherein the CTSGRS predicts the CTS at an ad libitum meal for each subject from the population of subjects by combining gene risk scores and the non-genetic and / or multi-omic data obtained from the samples obtained from each subject from the population of subjects;(f) training by the computer system, the CTSGRS from step (e) using a classifier model to generate the CTSGRS classifier, wherein the CTSGRS classifier generates a value ranging from 0 to 1 for each subject, wherein a value of greater than or equal to 0.5 classifies the subject as having a high CTSGRS and a value of less than 0.5 classifies the subject as having a low CTSGRS, wherein a high CTSGRS indicates a probability of the subject being above the sex-specific 75thpercentile for CTS and a low CTSGRS indicates a probability of the subject being below the sex-specific 75th percentile for CTS.
12. The computer-implemented method of claim 11, wherein the SNP genotype information comprises risk allele dosage, major allele presence, occurrence of a de novo variant, insertion, deletion, or genome rearrangement.
13. The computer-implemented method of claim 11 or 12, wherein each gene in the set of genes is known or suspected to play a role in appetite regulation, energy expenditure, lipid metabolism or adipogenesis.
14. The computer-implemented method of claim 11 or 12, wherein the subject is obese if the subject possesses a body mass index (BMI) > 30 kg / m2with or without type 2 diabetes.
15. The computer-implemented method of any one of the above claims, wherein the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with a sensitivity and / or specificity is at least 65%.
16. The computer-implemented method of any one of the above claims, wherein the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with an area under the curve (AUC) of at least 0.65, 0.7, 0.75, 0.8, 0.85, 0.9 0.95 or 1.
17. The computer-implemented method of any one of claims 11-16, wherein the CTSGRS classifier predicts the subject as having a CTS above the sex-specific 75thpercentile for CTS with positive predictive value of at least 50%.
18. The computer-implemented method of any one of claims 11-17, wherein the contribution factor is a function or algorithm that utilizes output of a regression analysis or a statistical analysis on the genetic dataset on CTS,19. The computer-implemented method of claim 18, wherein the genetic dataset is selected from the group consisting of a genome-wide data microarray chip (GWAS), targeted sequencing (exome or targeted genetic panel), variants detected by qPCR through targeted amplification, whole genome sequencing (WGS), targeted or untargeted sequencing of genomic re-arrangement, deletions, duplications, repeat extensions, and complex genotyping of HLA and CYP genes.
20. The computer-implemented method of any one of claims 11-19, wherein the non- genetic data is anthropometric data selected from the group consisting of age, weight, height, sex and a combination thereof;21. The computer-implemented method of any one of claims 11-20, wherein the non- genetic and / or multi-omic data is selected from the group consisting of metabolomic data, proteomic data, peptidomic data, epigenetic data, microbiome data, results from one or more questionnaires and any combination thereof.
22. The computer-implemented method of claim 21, wherein the epigenetic data comprises a presence or absence of DNA modifications selected from the group consisting of DNA methylation modifications, histone modifications and chromatin accessibility.
23. The computer-implemented method of claim 21, wherein the microbiome data is from any site on the subject or the subject’s environment.
24. The computer-implemented method of any one of claims 11-23, wherein the machine learning model in step (e) comprises a forward feature selection, backward feature selection, or random feature sampling algorithm and a random forest, GBM, or linear predictor that iteratively selects CTS with a desired training accuracy using the random forest predictor in multiple rounds until no further training accuracy improvement is detected.
25. The computer-implemented method of any one of claims 11-24, wherein the GRSs of step (e) are normalized.
26. The computer-implemented method of any one of claims 11-25, wherein each GSR provides an indication of a role or effect that the SNPs located around or near a specific gene may play on the specific trait known to be associated with CTS at an ad libitum meal.
27. The computer-implemented method of any one of claims 11-26, wherein the SNP is considered to be located around or near a specific gene if it is located within at least 500,000 kilobases (kb) of the specific gene.
28. The computer-implemented method of any one of claims 11-27, further comprising performing a genetic analysis on a sample obtained from the subject suffering from or suspected of suffering from obesity prior to step (a).
29. The computer-implemented method of claim 28, wherein the genetic analysis comprises obtaining sequence reads from the whole or portions of the whole genome of the subject.
30. The computer-implemented method of claim 28, wherein the genetic analysis comprises performing a genotyping method selected from the group consisting of restriction fragment length polymorphism identification (RFLPI), random amplified polymorphic detection (RAPD), amplified fragment length polymorphism detection (AFLPD), polymerase chain reaction (PCR), DNA sequencing, RNA sequencing, allele specific oligonucleotide (ASO) probes, and hybridization to microarrays or beads.
31. The computer-implemented method of any one of claims 11-30, further comprising performing by the computer system, a regression analysis and / or statistical analysis on the genetic dataset from step (a) prior to step (b).
32. The computer-implemented method of claim 31, wherein the regression analysis is a ridge regression or least absolute shrinkage and selection operator (LASSO) regression.
33. The computer-implemented method of claim 31, wherein beta values obtained from the regression analyses performed on the genetic datasets are used to generate the contribution factor in step (b).
34. The computer-implemented method of claim 31, wherein the statistical analysis is a genome-wide association study (GWAS).
35. The computer-implemented method of claim 34, wherein p-values and / or regression weights obtained from the statistical analyses performed on the genetic datasets are used to generate the contribution factors in step (b).
36. A method of predicting weight gain in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a high CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to gain more weight over time as compared to a subject with a low CTSGRS.
37. The method of claim 36, wherein the subject has not been previously diagnosed with obesity.
38. The method of claim 36, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
39. The method of claim 36, wherein the subject has been previously diagnosed with obesity.
40. The method of claim 36, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
41. A method for treating weight gain in a subject in need thereof, the method comprising: a. determining if a subject possesses a high calorie to satiation (CTS) using the method of any one of claims 1-35 on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and b. administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS.
42. The method of claim 41, wherein the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample.
43. The method of claim 41 or 42, wherein the subject has not been previously diagnosed with obesity.
44. The method of claim 41 or 42, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
45. The method of claim 41 or 42, wherein the subject has been previously diagnosed with obesity.
46. The method of claim 41 or 42, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
47. The method of any one of claims 41-46, wherein the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate.
48. The method of any one of claims 41-46, wherein the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK- 2735, ecnoglutide, utresglutide, mazdutide, amy cretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS- 4571, ASC-30, TTP273, MLX-7006 and BMF-650.
49. A method of predicting obesity in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a high CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to develop obesity over time as compared to a subject with a low CTSGRS.
50. The method of claim 49, wherein the subject has not been previously diagnosed with obesity.
51. The method of claim 49, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
52. A method for treating obesity in a subject in need thereof, the method comprising: a. determining if a subject possesses a high calorie to satiation (CTS) using the method of any one of claims 1-35 on a genetic dataset obtained from a sample obtained from the subject to determine a CTSGRS for the subject; and b. administering a pharmacotherapy to the subject that does not include a GLP-1 agonist if the CTSGRS determined for the subject indicates that the subject has a high CTSGRS or administering a pharmacotherapy that may include a GLP-1 agonist to the subject if the CTSGRS determined for the subject indicates that the subject has a low CTSGRS.
53. The method of claim 52, wherein the sample is selected from the group consisting of a blood sample, a saliva sample, a urine sample, a breath sample, and a stool sample.
54. The method of claim 52 or 53, wherein the subject has not been previously diagnosed with obesity.
55. The method of claim 52 or 53, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
56. The method of claim 52 or 53, wherein the subject has been previously diagnosed with obesity.
57. The method of claim 52 or 53, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
58. The method of any one of claims 52-57, wherein the pharmacotherapy that does not include a GLP-1 agonist is phentermine / topiramate.
59. The method of any one of claims 52-57, wherein the pharmacotherapy that does include a GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK- 2735, ecnoglutide, utresglutide, mazdutide, amy cretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS- 4571, ASC-30, TTP273, MLX-7006 and BMF-650.
60. A method of predicting response to a bariatric surgical intervention in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a high CTSGRS indicates that the subject is above the sexspecific 75thpercentile for CTS and is predicted to respond to the bariatric surgical intervention as compared to a subject with a low CTSGRS.
61. The method of claim 60, wherein the response to the bariatric surgical intervention is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time.
62. The method of claim 61, wherein the defined period of time is at least 1, 2, 3, 4 or 5 years following the bariatric surgical intervention.
63. The method of claim 60, wherein the subject has not been previously diagnosed with obesity.
64. The method of claim 60, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
65. The method of claim 60, wherein the subject has been previously diagnosed with obesity.
66. The method of claim 60, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
67. A method of predicting response to endoscopic sleeve gastroplasty in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a low CTSGRS indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to endoscopic sleeve gastroplasty as compared to a subject with a high CTSGRS.
68. The method of claim 67, wherein the response to endoscopic sleeve gastroplasty is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time.
69. The method of claim 68, wherein the defined period of time is at least 1, 2, 3, 4, 5 or 6 months following endoscopic sleeve gastroplasty.
70. The method of claim 67, wherein the subject has not been previously diagnosed with obesity.
71. The method of claim 67, wherein the subject has a body mass index (BMI) of less than or equal to 30 kg / m2with or without type 2 diabetes.
72. The method of claim 67, wherein the subject has been previously diagnosed with obesity.
73. The method of claim 72, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes.
74. A method of treating diabetes or prediabetes in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a low CTSGRS indicates that the subject is below the sex-specific 75thpercentile for CTS and is predicted to respond to a GLP-1 receptor agonist as compared to a subject with a high CTSGRS; and administering the GLP-1 receptor agonist.
75. The method of claim 74, wherein the response to the GLP-1 receptor agonist is evidenced by a percentage of total body weight loss (TBWL) of at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time.
76. The method of claim 75, wherein the defined period of time is at least 1, 2, 3, 4 or 5 years following administration of the GLP-1 receptor agonist.
77. The method of claim 74, wherein the subject has been previously diagnosed with obesity.
78. The method of claim 74, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes.
79. The method of any one of claims 74-78, wherein the administration of the GLP-1 agonist to the subject results in a reduction in the subject’s HbAlc.
80. The method of claim 79, wherein the reduction in the subject’s HbAlc is at least 5%, 10%, 15%, 20%, 25% or 30% after a defined period of time.
81. The method of claim 80, wherein the defined period of time is at least 1, 2, 3, 4 or 5 years following administration of the GLP-1 receptor agonist.
82. The method of any one of claims 74-81, wherein the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF- 650.
83. A method of selecting a treatment for a subject suffering from diabetes or prediabetes, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a low CTSGRS indicates that the subject is a responder to a GLP-1 receptor agonist as compared to a subject with a high CTSGRS.
84. The method of claim 83, wherein the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF-650.
85. A method of predicting one or more adverse effects (AEs) in a subject, the method comprising determining the CTSGRS of a subject using the method of any one of claims 1-35, wherein a high CTSGRS indicates that the subject is above the sex-specific 75thpercentile for CTS and is predicted to experience adverse effects when administered a GLP-1 receptor agonist alone or in combination with one or more additional therapeutic agents as compared to a subject with a low CTSGRS.
86. The method of claim 85, wherein the one or more AE is selected from the group consisting of diarrhea, constipation, abdominal discomfort, nausea, vomiting, headache and any combination thereof.
87. The method of claim 85 or 86, wherein the AE is nausea.
88. The method of claim 85, wherein the AE is nausea and headache.
89. The method of claim 85, wherein the subject has been previously diagnosed with obesity.
90. The method of claim 85, wherein the subject has a body mass index (BMI) of greater than 30 kg / m2with or without type 2 diabetes or prediabetes.
91. The method of any one of claims 85-90, wherein the GLP-1 agonist is selected from the group consisting of exenatide, liraglutide, albiglutide, dulaglutide, lixisenatide, semaglutide, tirzepatide maritide, CT-388, VK-2735, ecnoglutide, utresglutide, mazdutide, amycretin, dapiglutide, retatrutide, cagrisema, cotadutide, orforglipron, pemyidutide, danuglipron, aleniglipron, AZN5004, RGT-075, KAI-7535, suryodutide, GZR18, KAI-9531, MDR-001, HDM1002, VCT220, CT-996, TERN-601, ID110521156, XW014, SAL0112, GS-4571, ASC-30, TTP273, MLX-7006 and BMF- 650.
92. The method of any one of claims 85-91, wherein the one or more additional therapeutic agents is an antihyperglycemic medications (OAMs).
93. The method of claim 92, wherein the 0 AM is insulin.
94. The method of any one of claims 85-91, wherein the one or more additional therapeutic agents is a sodium-glucose cotransporter-2 (SGLT-2) inhibitor.
95. The method of claim 94, wherein the SGLT-2 inhibitor is selected from the group consisting of empagliflozin, canagliflozin, dapagliflozin, ertugliflozin, bexagliflozin and sotagliflozin.
96. The method of any one of claims 85-91, wherein the one or more additional therapeutic agents is a dipeptidyl peptidase-4 (DPP-4) inhibitor.
97. The method of claim 96, wherein the DPP-4 inhibitor is selected from the group consisting of sitagliptin, linagliptin, saxagliptin and alogliptin.
Citation Information
Patent Citations
Polygenic risk scores for predicting disease complications and / or response to therapy
US20210254164A1
Machine learning platform for generating risk models
US20210375392A1
Polygenic risk score for in vitro fertilization
US20220367063A1