Sequencing of 640,000 exomes identifies GPR75 variants associated with protection from obesity

Introduction Obesity accounts for a substantial and growing burden of disease globally. Body adiposity is highly heritable, and human genetic studies can lead to biological and therapeutic insights. Rationale Whole-exome sequencing of hundreds of thousands of individuals is complementary to approaches used to date in obesity genetics and has the potential to identify rare protein-coding variants with large phenotypic impact. We sequenced the exomes of 645,626 individuals from the UK, the US, and Mexico and estimated associations of rare coding variants with body mass index (BMI), a measure of overall adiposity used to define obesity in clinical practice. We complemented exome sequencing with fine-mapping of common alleles, polygenic score analysis, and in vitro and in vivo modeling work. Results We identified 16 genes for which the burden of rare nonsynonymous variants was associated with BMI at exome-wide statistical significance (inverse-variance weighted meta-analysis P < 3.6 × 10−7), including associations at five brain-expressed G protein–coupled receptors (CALCR, MC4R, GIPR, GPR151, and GPR75). We observed an overrepresentation of genes highly expressed in the hypothalamus, a key center for the neuroendocrine regulation of energy balance. Protein-truncating variants in GPR75 were found in ~4/10,000 sequenced people and were associated with 1.8 kg/m2 lower BMI, 5.3 kg lower bodyweight, and 54% lower odds of obesity in heterozygous carriers. Knock out of Gpr75 in mice resulted in resistance to weight gain in a high-fat diet model, which was allele-dose dependent (25% and 44% lower weight gain, respectively, for heterozygous Gpr75−/+ mice and knockout Gpr75−/− mice compared with wild type) and accompanied by improved glycemic control and insulin sensitivity. Protein-truncating variants in CALCR were associated with higher BMI and obesity risk, whereas protein-truncating variants in GIPR and two missense alleles [Arg190→Gln (Arg190Gln), Glu288Gly], which we show result in loss of function in vitro, were associated with lower adiposity. Among monogenic obesity genes in the leptin-melanocortin pathway, heterozygous predicted loss-of-function variants in LEP, POMC, PCSK1, and MC4R (but not LEPR) were associated with higher BMI. Rare protein-truncating variants in UBR2, ANO4, and PCSK1 were associated with more than twofold higher odds of obesity in heterozygous carriers, similar to predicted-deleterious nonsynonymous variants in MC4R, which are considered the most common cause of monogenic obesity. Polygenic predisposition due to >2 million common genetic variants influenced the penetrance of obesity in rare variant carriers in an additive fashion. Conclusion These results suggest that inhibition of GPR75 may be a therapeutic strategy for obesity and illustrate the power of massive-scale exome sequencing for the identification of large-effect coding variant associations and drug targets for complex traits.