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

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

Genomics and phenomics of body mass index reveals a complex disease network.

Elevated body mass index (BMI) is heritable and associated with many health conditions that impact morbidity and mortality. The study of the genetic association of BMI across a broad range of common disease conditions offers the opportunity to extend current knowledge regarding the breadth and depth of adiposity-related diseases. We identify 906 (364 novel) and 41 (6 novel) genome-wide significant loci for BMI among participants of European (N~1.1 million) and African (N~100,000) ancestry, respectively. Using a BMI genetic risk score including 2446 variants, 316 diagnoses are associated in the Million Veteran Program, with 96.5% showing increased risk. A co-morbidity network analysis reveals seven disease communities containing multiple interconnected diseases associated with BMI as well as extensive connections across communities. Mendelian randomization analysis confirms numerous phenotypes across a breadth of organ systems, including conditions of the circulatory (heart failure, ischemic heart disease, atrial fibrillation), genitourinary (chronic renal failure), respiratory (respiratory failure, asthma), musculoskeletal and dermatologic systems that are deeply interconnected within and across the disease communities. This work shows that the complex genetic architecture of BMI associates with a broad range of major health conditions, supporting the need for comprehensive approaches to prevent and treat obesity.

Gene Sequencing Identifies Perturbation in Nitric Oxide Signaling as a Nonlipid Molecular Subtype of Coronary Artery Disease.

BACKGROUND: A key goal of precision medicine is to disaggregate common, complex diseases into discrete molecular subtypes. Rare coding variants in the low-density lipoprotein receptor gene (LDLR) are identified in 1% to 2% of coronary artery disease (CAD) patients, defining a molecular subtype with risk driven by hypercholesterolemia. METHODS: To search for additional subtypes, we compared the frequency of rare, predicted loss-of-function and damaging missense variants aggregated within a given gene in 41 081 CAD cases versus 217 115 controls. RESULTS: Rare variants in LDLR were most strongly associated with CAD, present in 1% of cases and associated with 4.4-fold increased CAD risk. A second subtype was characterized by variants in endothelial nitric oxide synthase gene (NOS3), a key enzyme regulating vascular tone, endothelial function, and platelet aggregation. A rare predicted loss-of-function or damaging missense variants in NOS3 was present in 0.6% of cases and associated with 2.42-fold increased risk of CAD (95% CI, 1.80-3.26; P=5.50×10-9). These variants were associated with higher systolic blood pressure (+3.25 mm Hg; [95% CI, 1.86-4.65]; P=5.00×10-6) and increased risk of hypertension (adjusted odds ratio 1.31; [95% CI, 1.14-1.51]; P=2.00×10-4) but not circulating cholesterol concentrations, suggesting that, beyond lipid pathways, nitric oxide synthesis is a key nonlipid driver of CAD risk. CONCLUSIONS: Beyond LDLR, we identified an additional nonlipid molecular subtype of CAD characterized by rare variants in the NOS3 gene.

Analyzing human knockouts to validate GPR151 as a therapeutic target for reduction of body mass index.

Novel drug targets for sustained reduction in body mass index (BMI) are needed to curb the epidemic of obesity, which affects 650 million individuals worldwide and is a causal driver of cardiovascular and metabolic disease and mortality. Previous studies reported that the Arg95Ter nonsense variant of GPR151, an orphan G protein-coupled receptor, is associated with reduced BMI and reduced risk of Type 2 Diabetes (T2D). Here, we further investigate GPR151 with the Pakistan Genome Resource (PGR), which is one of the largest exome biobanks of human homozygous loss-of-function carriers (knockouts) in the world. Among PGR participants, we identify eleven GPR151 putative loss-of-function (plof) variants, three of which are present at homozygosity (Arg95Ter, Tyr99Ter, and Phe175LeufsTer7), with a cumulative allele frequency of 2.2%. We confirm these alleles in vitro as loss-of-function. We test if GPR151 plof is associated with BMI, T2D, or other metabolic traits and find that GPR151 deficiency in complete human knockouts is not associated with clinically significant differences in these traits. Relative to Gpr151+/+ mice, Gpr151-/- animals exhibit no difference in body weight on normal chow and higher body weight on a high-fat diet. Together, our findings indicate that GPR151 antagonism is not a compelling therapeutic approach to treatment of obesity.

Identification of genetic effects underlying type 2 diabetes in South Asian and European populations.

South Asians are at high risk of developing type 2 diabetes (T2D). We carried out a genome-wide association meta-analysis with South Asian T2D cases (n = 16,677) and controls (n = 33,856), followed by combined analyses with Europeans (neff = 231,420). We identify 21 novel genetic loci for significant association with T2D (P = 4.7 × 10-8 to 5.2 × 10-12), to the best of our knowledge at the point of analysis. The loci are enriched for regulatory features, including DNA methylation and gene expression in relevant tissues, and highlight CHMP4B, PDHB, LRIG1 and other genes linked to adiposity and glucose metabolism. A polygenic risk score based on South Asian-derived summary statistics shows ~4-fold higher risk for T2D between the top and bottom quartile. Our results provide further insights into the genetic mechanisms underlying T2D, and highlight the opportunities for discovery from joint analysis of data from across ancestral populations.

Genome-wide analysis of blood lipid metabolites in over 5000 South Asians reveals biological insights at cardiometabolic disease loci.

BACKGROUND: Genetic, lifestyle, and environmental factors can lead to perturbations in circulating lipid levels and increase the risk of cardiovascular and metabolic diseases. However, how changes in individual lipid species contribute to disease risk is often unclear. Moreover, little is known about the role of lipids on cardiovascular disease in Pakistan, a population historically underrepresented in cardiovascular studies. METHODS: We characterised the genetic architecture of the human blood lipidome in 5662 hospital controls from the Pakistan Risk of Myocardial Infarction Study (PROMIS) and 13,814 healthy British blood donors from the INTERVAL study. We applied a candidate causal gene prioritisation tool to link the genetic variants associated with each lipid to the most likely causal genes, and Gaussian Graphical Modelling network analysis to identify and illustrate relationships between lipids and genetic loci. RESULTS: We identified 253 genetic associations with 181 lipids measured using direct infusion high-resolution mass spectrometry in PROMIS, and 502 genetic associations with 244 lipids in INTERVAL. Our analyses revealed new biological insights at genetic loci associated with cardiometabolic diseases, including novel lipid associations at the LPL, MBOAT7, LIPC, APOE-C1-C2-C4, SGPP1, and SPTLC3 loci. CONCLUSIONS: Our findings, generated using a distinctive lipidomics platform in an understudied South Asian population, strengthen and expand the knowledge base of the genetic determinants of lipids and their association with cardiometabolic disease-related loci.

Genetic analysis in European ancestry individuals identifies 517 loci associated with liver enzymes.

Serum concentration of hepatic enzymes are linked to liver dysfunction, metabolic and cardiovascular diseases. We perform genetic analysis on serum levels of alanine transaminase (ALT), alkaline phosphatase (ALP) and gamma-glutamyl transferase (GGT) using data on 437,438 UK Biobank participants. Replication in 315,572 individuals from European descent from the Million Veteran Program, Rotterdam Study and Lifeline study confirms 517 liver enzyme SNPs. Genetic risk score analysis using the identified SNPs is strongly associated with serum activity of liver enzymes in two independent European descent studies (The Airwave Health Monitoring study and the Northern Finland Birth Cohort 1966). Gene-set enrichment analysis using the identified SNPs highlights involvement in liver development and function, lipid metabolism, insulin resistance, and vascular formation. Mendelian randomization analysis shows association of liver enzyme variants with coronary heart disease and ischemic stroke. Genetic risk score for elevated serum activity of liver enzymes is associated with higher fat percentage of body, trunk, and liver and body mass index. Our study highlights the role of molecular pathways regulated by the liver in metabolic disorders and cardiovascular disease.

Correction: A missense variant in Mitochondrial Amidoxime Reducing Component 1 gene and protection against liver disease.

[This corrects the article DOI: 10.1371/journal.pgen.1008629.].

Deficiency of macrophage PHACTR1 impairs efferocytosis and promotes atherosclerotic plaque necrosis.

Efferocytosis, the process through which apoptotic cells (ACs) are cleared through actin-mediated engulfment by macrophages, prevents secondary necrosis, suppresses inflammation, and promotes resolution. Impaired efferocytosis drives the formation of clinically dangerous necrotic atherosclerotic plaques, the underlying etiology of coronary artery disease (CAD). An intron of the gene encoding PHACTR1 contains rs9349379 (A>G), a common variant associated with CAD. As PHACTR1 is an actin-binding protein, we reasoned that if the rs9349379 risk allele G causes lower PHACTR1 expression in macrophages, it might link the risk allele to CAD via impaired efferocytosis. We show here that rs9349379-G/G was associated with lower levels of PHACTR1 and impaired efferocytosis in human monocyte-derived macrophages and human atherosclerotic lesional macrophages compared with rs9349379-A/A. Silencing PHACTR1 in human and mouse macrophages compromised AC engulfment, and Western diet-fed Ldlr-/- mice in which hematopoietic Phactr1 was genetically targeted showed impaired lesional efferocytosis, increased plaque necrosis, and thinner fibrous caps - all signs of vulnerable plaques in humans. Mechanistically, PHACTR1 prevented dephosphorylation of myosin light chain (MLC), which was necessary for AC engulfment. In summary, rs9349379-G lowered PHACTR1, which, by lowering phospho-MLC, compromised efferocytosis. Thus, rs9349379-G may contribute to CAD risk, at least in part, by impairing atherosclerotic lesional macrophage efferocytosis.

Heterozygous ABCG5 Gene Deficiency and Risk of Coronary Artery Disease.

BACKGROUND: Familial sitosterolemia is a rare Mendelian disorder characterized by hyperabsorption and decreased biliary excretion of dietary sterols. Affected individuals typically have complete genetic deficiency-homozygous loss-of-function (LoF) variants-in the ABCG5 or ABCG8 genes and have substantially elevated plasma sitosterol and LDL (low-density lipoprotein) cholesterol (LDL-C) levels. The impact of partial genetic deficiency of ABCG5 or ABCG8-as occurs in heterozygous carriers of LoF variants-on LDL-C and risk of coronary artery disease (CAD) has remained uncertain. METHODS: We first recruited 9 sitosterolemia families, identified causative LoF variants in ABCG5 or ABCG8, and evaluated the associations of these ABCG5 or ABCG8 LoF variants with plasma phytosterols and lipid levels. We next assessed for LoF variants in ABCG5 or ABCG8 in CAD cases (n=29 321) versus controls (n=357 326). We tested the association of rare LoF variants in ABCG5 or ABCG8 with blood lipids and risk for CAD. Rare LoF variants were defined as protein-truncating variants with minor allele frequency <0.1% in ABCG5 or ABCG8. RESULTS: In sitosterolemia families, 7 pedigrees harbored causative LoF variants in ABCG5 and 2 pedigrees in ABCG8. Homozygous LoF variants in either ABCG5 or ABCG8 led to marked elevations in sitosterol and LDL-C. Of those sitosterolemia families, heterozygous carriers of ABCG5 LoF variants exhibited increased sitosterol and LDL-C levels compared with noncarriers. Within large-scale CAD case-control cohorts, prevalence of rare LoF variants in ABCG5 and in ABCG8 was ≈0.1% each. ABCG5 heterozygous LoF variant carriers had significantly elevated LDL-C levels (25 mg/dL [95% CI, 14-35]; P=1.1×10-6) and were at 2-fold increased risk of CAD (odds ratio, 2.06 [95% CI, 1.27-3.35]; P=0.004). By contrast, ABCG8 heterozygous LoF carrier status was not associated with increased LDL-C or risk of CAD. CONCLUSIONS: Although familial sitosterolemia is traditionally considered as a recessive disorder, we observed that heterozygous carriers of an LoF variant in ABCG5 had significantly increased sitosterol and LDL-C levels and a 2-fold increase in risk of CAD.

Validating a non-invasive, ALT-based non-alcoholic fatty liver phenotype in the million veteran program.

BACKGROUND & AIMS: Given ongoing challenges in non-invasive non-alcoholic liver disease (NAFLD) diagnosis, we sought to validate an ALT-based NAFLD phenotype using measures readily available in electronic health records (EHRs) and population-based studies by leveraging the clinical and genetic data in the Million Veteran Program (MVP), a multi-ethnic mega-biobank of US Veterans. METHODS: MVP participants with alanine aminotransferases (ALT) >40 units/L for men and >30 units/L for women without other causes of liver disease were compared to controls with normal ALT. Genetic variants spanning eight NAFLD risk or ALT-associated loci (LYPLAL1, GCKR, HSD17B13, TRIB1, PPP1R3B, ERLIN1, TM6SF2, PNPLA3) were tested for NAFLD associations with sensitivity analyses adjusting for metabolic risk factors and alcohol consumption. A manual EHR review assessed performance characteristics of the NAFLD phenotype with imaging and biopsy data as gold standards. Genetic associations with advanced fibrosis were explored using FIB4, NAFLD Fibrosis Score and platelet counts. RESULTS: Among 322,259 MVP participants, 19% met non-invasive criteria for NAFLD. Trans-ethnic meta-analysis replicated associations with previously reported genetic variants in all but LYPLAL1 and GCKR loci (P<6x10-3), without attenuation when adjusted for metabolic risk factors and alcohol consumption. At the previously reported LYPLAL1 locus, the established genetic variant did not appear to be associated with NAFLD, however the regional association plot showed a significant association with NAFLD 279kb downstream. In the EHR validation, the ALT-based NAFLD phenotype yielded a positive predictive value 0.89 and 0.84 for liver biopsy and abdominal imaging, respectively (inter-rater reliability (Cohen's kappa = 0.98)). HSD17B13 and PNPLA3 loci were associated with advanced fibrosis. CONCLUSIONS: We validate a simple, non-invasive ALT-based NAFLD phenotype using EHR data by leveraging previously established NAFLD risk-associated genetic polymorphisms.

Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis.

We investigated type 2 diabetes (T2D) genetic susceptibility via multi-ancestry meta-analysis of 228,499 cases and 1,178,783 controls in the Million Veteran Program (MVP), DIAMANTE, Biobank Japan and other studies. We report 568 associations, including 286 autosomal, 7 X-chromosomal and 25 identified in ancestry-specific analyses that were previously unreported. Transcriptome-wide association analysis detected 3,568 T2D associations with genetically predicted gene expression in 687 novel genes; of these, 54 are known to interact with FDA-approved drugs. A polygenic risk score (PRS) was strongly associated with increased risk of T2D-related retinopathy and modestly associated with chronic kidney disease (CKD), peripheral artery disease (PAD) and neuropathy. We investigated the genetic etiology of T2D-related vascular outcomes in the MVP and observed statistical SNP-T2D interactions at 13 variants, including coronary heart disease (CHD), CKD, PAD and neuropathy. These findings may help to identify potential therapeutic targets for T2D and genomic pathways that link T2D to vascular outcomes.

The effect of LRRK2 loss-of-function variants in humans.

Human genetic variants predicted to cause loss-of-function of protein-coding genes (pLoF variants) provide natural in vivo models of human gene inactivation and can be valuable indicators of gene function and the potential toxicity of therapeutic inhibitors targeting these genes1,2. Gain-of-kinase-function variants in LRRK2 are known to significantly increase the risk of Parkinson's disease3,4, suggesting that inhibition of LRRK2 kinase activity is a promising therapeutic strategy. While preclinical studies in model organisms have raised some on-target toxicity concerns5-8, the biological consequences of LRRK2 inhibition have not been well characterized in humans. Here, we systematically analyze pLoF variants in LRRK2 observed across 141,456 individuals sequenced in the Genome Aggregation Database (gnomAD)9, 49,960 exome-sequenced individuals from the UK Biobank and over 4 million participants in the 23andMe genotyped dataset. After stringent variant curation, we identify 1,455 individuals with high-confidence pLoF variants in LRRK2. Experimental validation of three variants, combined with previous work10, confirmed reduced protein levels in 82.5% of our cohort. We show that heterozygous pLoF variants in LRRK2 reduce LRRK2 protein levels but that these are not strongly associated with any specific phenotype or disease state. Our results demonstrate the value of large-scale genomic databases and phenotyping of human loss-of-function carriers for target validation in drug discovery.

A missense variant in Mitochondrial Amidoxime Reducing Component 1 gene and protection against liver disease.

Analyzing 12,361 all-cause cirrhosis cases and 790,095 controls from eight cohorts, we identify a common missense variant in the Mitochondrial Amidoxime Reducing Component 1 gene (MARC1 p.A165T) that associates with protection from all-cause cirrhosis (OR 0.91, p = 2.3*10-11). This same variant also associates with lower levels of hepatic fat on computed tomographic imaging and lower odds of physician-diagnosed fatty liver as well as lower blood levels of alanine transaminase (-0.025 SD, 3.7*10-43), alkaline phosphatase (-0.025 SD, 1.2*10-37), total cholesterol (-0.030 SD, p = 1.9*10-36) and LDL cholesterol (-0.027 SD, p = 5.1*10-30) levels. We identified a series of additional MARC1 alleles (low-frequency missense p.M187K and rare protein-truncating p.R200Ter) that also associated with lower cholesterol levels, liver enzyme levels and reduced risk of cirrhosis (0 cirrhosis cases for 238 R200Ter carriers versus 17,046 cases of cirrhosis among 759,027 non-carriers, p = 0.04) suggesting that deficiency of the MARC1 enzyme may lower blood cholesterol levels and protect against cirrhosis.

The promise and reality of therapeutic discovery from large cohorts.

Technological advances in rapid data acquisition have transformed medical biology into a data mining field, where new data sets are routinely dissected and analyzed by statistical models of ever-increasing complexity. Many hypotheses can be generated and tested within a single large data set, and even small effects can be statistically discriminated from a sea of noise. On the other hand, the development of therapeutic interventions moves at a much slower pace. They are determined from carefully randomized and well-controlled experiments with explicitly stated outcomes as the principal mechanism by which a single hypothesis is tested. In this paradigm, only a small fraction of interventions can be tested, and an even smaller fraction are ultimately deemed therapeutically successful. In this Review, we propose strategies to leverage large-cohort data to inform the selection of targets and the design of randomized trials of novel therapeutics. Ultimately, the incorporation of big data and experimental medicine approaches should aim to reduce the failure rate of clinical trials as well as expedite and lower the cost of drug development.

Genetic Interleukin 6 Signaling Deficiency Attenuates Cardiovascular Risk in Clonal Hematopoiesis.

BACKGROUND: Clonal hematopoiesis of indeterminate potential (CHIP) refers to clonal expansion of hematopoietic stem cells attributable to acquired leukemic mutations in genes such as DNMT3A or TET2. In humans, CHIP associates with prevalent myocardial infarction. In mice, CHIP accelerates atherosclerosis and increases IL-6/IL-1ß expression, raising the hypothesis that IL-6 pathway antagonism in CHIP carriers would decrease cardiovascular disease (CVD) risk. METHODS: We analyzed exome sequences from 35 416 individuals in the UK Biobank without prevalent CVD, to identify participants with DNMT3A or TET2 CHIP. We used the IL6R p.Asp358Ala coding mutation as a genetic proxy for IL-6 inhibition. We tested the association of CHIP status with incident CVD events (myocardial infarction, coronary revascularization, stroke, or death), and whether it was modified by IL6R p.Asp358Ala. RESULTS: We identified 1079 (3.0%) individuals with CHIP, including 432 (1.2%) with large clones (allele fraction >10%). During 6.9-year median follow-up, CHIP associated with increased incident CVD event risk (hazard ratio, 1.27 [95% CI, 1.04-1.56], P=0.019), with greater risk from large CHIP clones (hazard ratio, 1.59 [95% CI, 1.21-2.09], P<0.001). IL6R p.Asp358Ala attenuated CVD event risk among participants with large CHIP clones (hazard ratio, 0.46 [95% CI, 0.29-0.73], P<0.001) but not in individuals without CHIP (hazard ratio, 0.95 [95% CI, 0.89-1.01], P=0.08; Pinteraction=0.003). In 9951 independent participants, the association of CHIP status with myocardial infarction similarly varied by IL6R p.Asp358Ala (Pinteraction=0.036). CONCLUSIONS: CHIP is associated with increased risk of incident CVD. Among carriers of large CHIP clones, genetically reduced IL-6 signaling abrogated this risk.

Genome-wide association analysis of venous thromboembolism identifies new risk loci and genetic overlap with arterial vascular disease.

Venous thromboembolism is a significant cause of mortality1, yet its genetic determinants are incompletely defined. We performed a discovery genome-wide association study in the Million Veteran Program and UK Biobank, with testing of approximately 13 million DNA sequence variants for association with venous thromboembolism (26,066 cases and 624,053 controls) and meta-analyzed both studies, followed by independent replication with up to 17,672 venous thromboembolism cases and 167,295 controls. We identified 22 previously unknown loci, bringing the total number of venous thromboembolism-associated loci to 33, and subsequently fine-mapped these associations. We developed a genome-wide polygenic risk score for venous thromboembolism that identifies 5% of the population at an equivalent incident venous thromboembolism risk to carriers of the established factor V Leiden p.R506Q and prothrombin G20210A mutations. Our data provide mechanistic insights into the genetic epidemiology of venous thromboembolism and suggest a greater overlap among venous and arterial cardiovascular disease than previously thought.