Macrophage WDFY3, a protector against autoimmunity.

Efficient efferocytosis is essential for maintaining homeostasis. Excessive apoptotic cell (AC) death and impaired macrophage efferocytosis lead to autoantigen release and autoantibody production, immune activation, and organ damage. It remains unclear whether these immunogenic autoantigens are the sole cause of increased autoimmunity or if efferocytosis of ACs directly influences macrophage function, impacting their ability to activate T cells and potentially amplifying autoimmune responses. Additionally, it has not been established if enhancing macrophage efferocytosis or modulating macrophage responses to AC engulfment can be protective in autoimmune-like disorders. Our previous work showed WDFY3 is crucial for efficient macrophage efferocytosis. This study reveals that myeloid knockout of Wdfy3 exacerbates autoimmunity in young mice with increased AC burden by systemic injections of ACs and in middle-aged mice developing spontaneous autoimmunity, whereas ectopic overexpression of WDFY3 suppresses autoimmunity in these models. Macrophages, as efferocytes, can activate T cells and the inflammasome upon engulfing ACs, which are suppressed by overexpressing WDFY3. This work uncovered the role of WDFY3 as a protector against autoimmunity by promoting macrophage efferocytosis thus limiting autoantigen production, as well as mitigating T cell activation and inflammasome activation.

NOTCH3 p.Arg1231Cys is markedly enriched in South Asians and associated with stroke.

The genetic factors of stroke in South Asians are largely unexplored. Exome-wide sequencing and association analysis (ExWAS) in 75 K Pakistanis identified NM_000435.3(NOTCH3):c.3691 C > T, encoding the missense amino acid substitution p.Arg1231Cys, enriched in South Asians (alternate allele frequency = 0.58% compared to 0.019% in Western Europeans), and associated with subcortical hemorrhagic stroke [odds ratio (OR) = 3.39, 95% confidence interval (CI) = [2.26, 5.10], p = 3.87 × 10-9), and all strokes (OR [CI] = 2.30 [1.77, 3.01], p = 7.79 × 10-10). NOTCH3 p.Arg231Cys was strongly associated with white matter hyperintensity on MRI in United Kingdom Biobank (UKB) participants (effect [95% CI] in SD units = 1.1 [0.61, 1.5], p = 3.0 × 10-6). The variant is attributable for approximately 2.0% of hemorrhagic strokes and 1.1% of all strokes in South Asians. These findings highlight the value of diversity in genetic studies and have major implications for genomic medicine and therapeutic development in South Asian populations.

Identification and characterization of human GDF15 knockouts.

Growth differentiation factor 15 (GDF15) is a secreted protein that regulates food intake, body weight and stress responses in pre-clinical models1. The physiological function of GDF15 in humans remains unclear. Pharmacologically, GDF15 agonism in humans causes nausea without accompanying weight loss2, and GDF15 antagonism is being tested in clinical trials to treat cachexia and anorexia. Human genetics point to a role for GDF15 in hyperemesis gravidarum, but the safety or impact of complete GDF15 loss, particularly during pregnancy, is unknown3-7. Here we show the absence of an overt phenotype in human GDF15 loss-of-function carriers, including stop gains, frameshifts and the fully inactivating missense variant C211G3. These individuals were identified from 75,018 whole-exome/genome-sequenced participants in the Pakistan Genomic Resource8,9 and recall-by-genotype studies with family-based recruitment of variant carrier probands. We describe 8 homozygous ('knockouts') and 227 heterozygous carriers of loss-of-function alleles, including C211G. GDF15 knockouts range in age from 31 to 75 years, are fertile, have multiple children and show no consistent overt phenotypes, including metabolic dysfunction. Our data support the hypothesis that GDF15 is not required for fertility, healthy pregnancy, foetal development or survival into adulthood. These observations support the safety of therapeutics that block GDF15.

Low MBOAT7 expression, a genetic risk for MASH, promotes a profibrotic pathway involving hepatocyte TAZ upregulation.

BACKGROUND AND AIMS: The common genetic variant rs641738 C>T is a risk factor for metabolic dysfunction-associated steatotic liver disease and metabolic dysfunction-associated steatohepatitis (MASH), including liver fibrosis, and is associated with decreased expression of the phospholipid-remodeling enzyme MBOAT7 (LPIAT1). However, whether restoring MBOAT7 expression in established metabolic dysfunction-associated steatotic liver disease dampens the progression to liver fibrosis and, importantly, the mechanism through which decreased MBOAT7 expression exacerbates MASH fibrosis remain unclear. APPROACH AND RESULTS: We first showed that hepatocyte MBOAT7 restoration in mice with diet-induced steatohepatitis slows the progression to liver fibrosis. Conversely, when hepatocyte-MBOAT7 was silenced in mice with established hepatosteatosis, liver fibrosis but not hepatosteatosis was exacerbated. Mechanistic studies revealed that hepatocyte-MBOAT7 restoration in MASH mice lowered hepatocyte-TAZ (WWTR1), which is known to promote MASH fibrosis. Conversely, hepatocyte-MBOAT7 silencing enhanced TAZ upregulation in MASH. Finally, we discovered that changes in hepatocyte phospholipids due to MBOAT7 loss-of-function promote a cholesterol trafficking pathway that upregulates TAZ and the TAZ-induced profibrotic factor Indian hedgehog (IHH). As evidence for relevance in humans, we found that the livers of individuals with MASH carrying the rs641738-T allele had higher hepatocyte nuclear TAZ, indicating higher TAZ activity and increased IHH mRNA. CONCLUSIONS: This study provides evidence for a novel mechanism linking MBOAT7-LoF to MASH fibrosis, adds new insight into an established genetic locus for MASH, and, given the druggability of hepatocyte TAZ for MASH fibrosis, suggests a personalized medicine approach for subjects at increased risk for MASH fibrosis due to inheritance of variants that lower MBOAT7.

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.

Whole-genome sequencing uncovers two loci for coronary artery calcification and identifies ARSE as a regulator of vascular calcification.

Coronary artery calcification (CAC) is a measure of atherosclerosis and a well-established predictor of coronary artery disease (CAD) events. Here we describe a genome-wide association study (GWAS) of CAC in 22,400 participants from multiple ancestral groups. We confirmed associations with four known loci and identified two additional loci associated with CAC (ARSE and MMP16), with evidence of significant associations in replication analyses for both novel loci. Functional assays of ARSE and MMP16 in human vascular smooth muscle cells (VSMCs) demonstrate that ARSE is a promoter of VSMC calcification and VSMC phenotype switching from a contractile to a calcifying or osteogenic phenotype. Furthermore, we show that the association of variants near ARSE with reduced CAC is likely explained by reduced ARSE expression with the G allele of enhancer variant rs5982944. Our study highlights ARSE as an important contributor to atherosclerotic vascular calcification, and a potential drug target for vascular calcific disease.

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.