Yield of genetic association signals from genomes, exomes and imputation in the UK Biobank.

Whole-genome sequencing (WGS), whole-exome sequencing (WES) and array genotyping with imputation (IMP) are common strategies for assessing genetic variation and its association with medically relevant phenotypes. To date, there has been no systematic empirical assessment of the yield of these approaches when applied to hundreds of thousands of samples to enable the discovery of complex trait genetic signals. Using data for 100 complex traits from 149,195 individuals in the UK Biobank, we systematically compare the relative yield of these strategies in genetic association studies. We find that WGS and WES combined with arrays and imputation (WES + IMP) have the largest association yield. Although WGS results in an approximately fivefold increase in the total number of assayed variants over WES + IMP, the number of detected signals differed by only 1% for both single-variant and gene-based association analyses. Given that WES + IMP typically results in savings of lab and computational time and resources expended per sample, we evaluate the potential benefits of applying WES + IMP to larger samples. When we extend our WES + IMP analyses to 468,169 UK Biobank individuals, we observe an approximately fourfold increase in association signals with the threefold increase in sample size. We conclude that prioritizing WES + IMP and large sample sizes rather than contemporary short-read WGS alternatives will maximize the number of discoveries in genetic association studies.

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.

Agonist antibody to guanylate cyclase receptor NPR1 regulates vascular tone.

Heart failure is a leading cause of morbidity and mortality1,2. Elevated intracardiac pressures and myocyte stretch in heart failure trigger the release of counter-regulatory natriuretic peptides, which act through their receptor (NPR1) to affect vasodilation, diuresis and natriuresis, lowering venous pressures and relieving venous congestion3-8. Recombinant natriuretic peptide infusions were developed to treat heart failure but have been limited by a short duration of effect9,10. Here we report that in a human genetic analysis of over 700,000 individuals, lifelong exposure to coding variants of the NPR1 gene is associated with changes in blood pressure and risk of heart failure. We describe the development of REGN5381, an investigational monoclonal agonist antibody that targets the membrane-bound guanylate cyclase receptor NPR1. REGN5381, an allosteric agonist of NPR1, induces an active-like receptor conformation that results in haemodynamic effects preferentially on venous vasculature, including reductions in systolic blood pressure and venous pressure in animal models. In healthy human volunteers, REGN5381 produced the expected haemodynamic effects, reflecting reductions in venous pressures, without obvious changes in diuresis and natriuresis. These data support the development of REGN5381 for long-lasting and selective lowering of venous pressures that drive symptomatology in patients with heart failure.

Deficiency of Peptidylglycine-alpha-amidating Monooxygenase, a Cause of Sarcopenic Diabetes Mellitus.

CONTEXT: Peptidylglycine-α-amidating monooxygenase (PAM) is a critical enzyme in the endocrine system responsible for activation, by amidation, of bioactive peptides. OBJECTIVE: To define the clinical phenotype of carriers of genetic mutations associated with impaired PAM-amidating activity (PAM-AMA). DESIGN: We used genetic and phenotypic data from cohort studies: the Malmö Diet and Cancer (MDC; 1991-1996; reexamination in 2002-2012), the Malmö Preventive Project (MPP; 2002-2006), and the UK Biobank (UKB; 2012). SETTING: Exome-wide association analysis was used to identify loss-of-function (LoF) variants associated with reduced PAM-AMA and subsequently used for association with the outcomes. PATIENTS OR OTHER PARTICIPANTS: This study included n∼4500 participants from a subcohort of the MDC (MDC-Cardiovascular cohort), n∼4500 from MPP, and n∼300,000 from UKB. MAIN OUTCOME MEASURES: Endocrine-metabolic traits suggested by prior literature, muscle mass, muscle function, and sarcopenia. RESULTS: Two LoF variants in the PAM gene, Ser539Trp (minor allele frequency: 0.7%) and Asp563Gly (5%), independently contributed to a decrease of 2.33 [95% confidence interval (CI): 2.52/2.15; P = 2.5E-140] and 0.98 (1.04/0.92; P = 1.12E-225) SD units of PAM-AMA, respectively. The cumulative effect of the LoF was associated with diabetes, reduced insulin secretion, and higher levels of GH and IGF-1. Moreover, carriers had reduced muscle mass and function, followed by a higher risk of sarcopenia. Indeed, the Ser539Trp mutation increased the risk of sarcopenia by 30% (odds ratio 1.31; 95% CI: 1.16/1.47; P = 9.8E-06), independently of age and diabetes. CONCLUSION: PAM-AMA genetic deficiency results in a prediabetic sarcopenic phenotype. Early identification of PAM LoF carriers would allow targeted exercise interventions and calls for novel therapies that restore enzymatic activity.

Rare and Common Genetic Variation Underlying Atrial Fibrillation Risk.

Importance: Atrial fibrillation (AF) has a substantial genetic component. The importance of polygenic risk is well established, while the contribution of rare variants to disease risk warrants characterization in large cohorts. Objective: To identify rare predicted loss-of-function (pLOF) variants associated with AF and elucidate their role in risk of AF, cardiomyopathy (CM), and heart failure (HF) in combination with a polygenic risk score (PRS). Design, Setting, and Participants: This was a genetic association and nested case-control study. The impact of rare pLOF variants was evaluated on the risk of incident AF. HF and CM were assessed in cause-specific Cox regressions. End of follow-up was July 1, 2022. Data were analyzed from January to October 2023. The UK Biobank enrolled 502 480 individuals aged 40 to 69 years at inclusion in the United Kingdom between March 13, 2006, and October 1, 2010. UK residents of European ancestry were included. Individuals with prior diagnosis of AF were excluded from analyses of incident AF. Exposures: Rare pLOF variants and an AF PRS. Main Outcomes and Measures: Risk of AF and incident HF or CM prior to and subsequent to AF diagnosis. Results: A total of 403 990 individuals (218 489 [54.1%] female) with a median (IQR) age of 58 (51-63) years were included; 24 447 were diagnosed with incident AF over a median (IQR) follow-up period of 13.3 (12.4-14.0) years. Rare pLOF variants in 6 genes (TTN, RPL3L, PKP2, CTNNA3, KDM5B, and C10orf71) were associated with AF. Of these, TTN, RPL3L, PKP2, CTNNA3, and KDM5B replicated in an external cohort. Combined with high PRS, rare pLOF variants conferred an odds ratio of 7.08 (95% CI, 6.03-8.28) for AF. Carriers with high PRS also had a substantial 10-year risk of AF (16% in female individuals and 24% in male individuals older than 60 years). Rare pLOF variants were associated with increased risk of CM both prior to AF (hazard ratio [HR], 3.13; 95% CI, 2.24-4.36) and subsequent to AF (HR, 2.98; 95% CI, 1.89-4.69). Conclusions and Relevance: Rare and common genetic variation were associated with an increased risk of AF. The findings provide insights into the genetic underpinnings of AF and may aid in future genetic risk stratification.

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Divergent role of Mitochondrial Amidoxime Reducing Component 1 (MARC1) in human and mouse.

Recent human genome-wide association studies have identified common missense variants in MARC1, p.Ala165Thr and p.Met187Lys, associated with lower hepatic fat, reduction in liver enzymes and protection from most causes of cirrhosis. Using an exome-wide association study we recapitulated earlier MARC1 p.Ala165Thr and p.Met187Lys findings in 540,000 individuals from five ancestry groups. We also discovered novel rare putative loss of function variants in MARC1 with a phenotype similar to MARC1 p.Ala165Thr/p.Met187Lys variants. In vitro studies of recombinant human MARC1 protein revealed Ala165Thr substitution causes protein instability and aberrant localization in hepatic cells, suggesting MARC1 inhibition or deletion may lead to hepatoprotection. Following this hypothesis, we generated Marc1 knockout mice and evaluated the effect of Marc1 deletion on liver phenotype. Unexpectedly, our study found that whole-body Marc1 deficiency in mouse is not protective against hepatic triglyceride accumulation, liver inflammation or fibrosis. In attempts to explain the lack of the observed phenotype, we discovered that Marc1 plays only a minor role in mouse liver while its paralogue Marc2 is the main Marc family enzyme in mice. Our findings highlight the major difference in MARC1 physiological function between human and mouse.

Pharmacologic and Genetic Downregulation of Proprotein Convertase Subtilisin/Kexin Type 9 and Survival From Sepsis.

OBJECTIVES: Treatments that prevent sepsis complications are needed. Circulating lipid and protein assemblies-lipoproteins play critical roles in clearing pathogens from the bloodstream. We investigated whether early inhibition of proprotein convertase subtilisin/kexin type 9 (PCSK9) may accelerate bloodstream clearance of immunogenic bacterial lipids and improve sepsis outcomes. DESIGN: Genetic and clinical epidemiology, and experimental models. SETTING: Human genetics cohorts, secondary analysis of a phase 3 randomized clinical trial enrolling patients with cardiovascular disease (Evaluation of Cardiovascular Outcomes After an Acute Coronary Syndrome During Treatment With Alirocumab [ODYSSEY OUTCOMES]; NCT01663402), and experimental murine models of sepsis. PATIENTS OR SUBJECTS: Nine human cohorts with sepsis (total n = 12,514) were assessed for an association between sepsis mortality and PCSK9 loss-of-function (LOF) variants. Incident or fatal sepsis rates were evaluated among 18,884 participants in a post hoc analysis of ODYSSEY OUTCOMES. C57BI/6J mice were used in Pseudomonas aeruginosa and Staphylococcus aureus bacteremia sepsis models, and in lipopolysaccharide-induced animal models. INTERVENTIONS: Observational human cohort studies used genetic PCSK9 LOF variants as instrumental variables. ODYSSEY OUTCOMES participants were randomized to alirocumab or placebo. Mice were administered alirocumab, a PCSK9 inhibitor, at 5 mg/kg or 25 mg/kg subcutaneously, or isotype-matched control, 48 hours prior to the induction of bacterial sepsis. Mice did not receive other treatments for sepsis. MEASUREMENTS AND MAIN RESULTS: Across human cohort studies, the effect estimate for 28-day mortality after sepsis diagnosis associated with genetic PCSK9 LOF was odds ratio = 0.86 (95% CI, 0.67-1.10; p = 0.24). A significant association was present in antibiotic-treated patients. In ODYSSEY OUTCOMES, sepsis frequency and mortality were infrequent and did not significantly differ by group, although both were numerically lower with alirocumab vs. placebo (relative risk of death from sepsis for alirocumab vs. placebo, 0.62; 95% CI, 0.32-1.20; p = 0.15). Mice treated with alirocumab had lower endotoxin levels and improved survival. CONCLUSIONS: PCSK9 inhibition may improve clinical outcomes in sepsis in preventive, pretreatment settings.

Converging evidence from exome sequencing and common variants implicates target genes for osteoporosis.

In this study, we leveraged the combined evidence of rare coding variants and common alleles to identify therapeutic targets for osteoporosis. We undertook a large-scale multiancestry exome-wide association study for estimated bone mineral density, which showed that the burden of rare coding alleles in 19 genes was associated with estimated bone mineral density (P < 3.6 × 10-7). These genes were highly enriched for a set of known causal genes for osteoporosis (65-fold; P = 2.5 × 10-5). Exome-wide significant genes had 96-fold increased odds of being the top ranked effector gene at a given GWAS locus (P = 1.8 × 10-10). By integrating proteomics Mendelian randomization evidence, we prioritized CD109 (cluster of differentiation 109) as a gene for which heterozygous loss of function is associated with higher bone density. CRISPR-Cas9 editing of CD109 in SaOS-2 osteoblast-like cell lines showed that partial CD109 knockdown led to increased mineralization. This study demonstrates that the convergence of common and rare variants, proteomics and CRISPR can highlight new bone biology to guide therapeutic development.

Activin E-ACVR1C cross talk controls energy storage via suppression of adipose lipolysis in mice.

Body fat distribution is a heritable risk factor for cardiovascular and metabolic disease. In humans, rare Inhibin beta E (INHBE, activin E) loss-of-function variants are associated with a lower waist-to-hip ratio and protection from type 2 diabetes. Hepatic fatty acid sensing promotes INHBE expression during fasting and in obese individuals, yet it is unclear how the hepatokine activin E governs body shape and energy metabolism. Here, we uncover activin E as a regulator of adipose energy storage. By suppressing ß-agonist-induced lipolysis, activin E promotes fat accumulation and adipocyte hypertrophy and contributes to adipose dysfunction in mice. Mechanistically, we demonstrate that activin E elicits its effect on adipose tissue through ACVR1C, activating SMAD2/3 signaling and suppressing PPARG target genes. Conversely, loss of activin E or ACVR1C in mice increases fat utilization, lowers adiposity, and drives PPARG-regulated gene signatures indicative of healthy adipose function. Our studies identify activin E-ACVR1C as a metabolic rheostat promoting liver-adipose cross talk to restrain excessive fat breakdown and preserve fat mass during prolonged fasting, a mechanism that is maladaptive in obese individuals.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Genetic variation of the blood coagulation regulator tissue factor pathway inhibitor and venous thromboembolism among middle-aged and older adults: A population-based cohort study.

Background: Tissue factor is the main initiator of blood coagulation, and tissue factor pathway inhibitor (TFPI) is the primary inhibitor of the initiation of blood coagulation.The genetic variation of TFPI and the relation to venous thromboembolism (VTE), that is, venous thrombosis and pulmonary embolism, remains to be clarified. This exome sequencing study aimed to determine the molecular epidemiology of the TFPI gene and the relation to VTE in a large population-based cohort of middle-aged and older adults. Methods: The exomes of TFPI were analyzed for variants in 28,794 subjects without previous VTE (born 1923-1950, 60% women), who participated in the Malmö Diet and Cancer Study (1991-1996). Patients were followed until the first event of VTE, death, or 2018. Qualifying variants were defined as loss-of-function or nonbenign (PolyPhen-2) missense variants with minor allele frequency less than 0.1%. Results: No common variant was associated with VTE. Nine rare variants (two loss-of-function and seven nonbenign missense) were classified as qualifying and included in collapsing analysis. Prevalence of qualifying variants was 0.09%. Five individuals with VTE compared to 17 individuals without VTE carried one qualifying variant. Cox multivariate regression analysis adjusted for age, sex, body mass index, systolic blood pressure, smoking and alcohol consumption, rs6025, rs1799963, and ancestry showed a hazard ratio of 2.9 (95% CI, 1.2-7.1) for rare qualifying variants. Conclusion: Rare qualifying TFPI variants were associated with VTE, suggesting that rare variants in TFPI contribute to the development of VTE. The qualifying TFPI gene variants were very rare, suggesting a constrained gene.

Rare and common genetic determinants of metabolic individuality and their effects on human health.

Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.

Increased soluble urokinase plasminogen activator levels modulate monocyte function to promote atherosclerosis.

People with kidney disease are disproportionately affected by atherosclerosis for unclear reasons. Soluble urokinase plasminogen activator receptor (suPAR) is an immune-derived mediator of kidney disease, levels of which are strongly associated with cardiovascular outcomes. We assessed suPAR's pathogenic involvement in atherosclerosis using epidemiologic, genetic, and experimental approaches. We found serum suPAR levels to be predictive of coronary artery calcification and cardiovascular events in 5,406 participants without known coronary disease. In a genome-wide association meta-analysis including over 25,000 individuals, we identified a missense variant in the plasminogen activator, urokinase receptor (PLAUR) gene (rs4760), confirmed experimentally to lead to higher suPAR levels. Mendelian randomization analysis in the UK Biobank using rs4760 indicated a causal association between genetically predicted suPAR levels and atherosclerotic phenotypes. In an experimental model of atherosclerosis, proprotein convertase subtilisin/kexin-9 (Pcsk9) transfection in mice overexpressing suPAR (suPARTg) led to substantially increased atherosclerotic plaques with necrotic cores and macrophage infiltration compared with those in WT mice, despite similar cholesterol levels. Prior to induction of atherosclerosis, aortas of suPARTg mice excreted higher levels of CCL2 and had higher monocyte counts compared with WT aortas. Aortic and circulating suPARTg monocytes exhibited a proinflammatory profile and enhanced chemotaxis. These findings characterize suPAR as a pathogenic factor for atherosclerosis acting at least partially through modulation of monocyte function.

Framework From a Multidisciplinary Approach for Transitioning Variants of Unknown Significance From Clinical Genetic Testing in Kidney Disease to a Definitive Classification.

Introduction: Monogenic causes in over 300 kidney-associated genes account for approximately 12% of end stage kidney disease (ESKD) cases. Advances in sequencing and large customized panels enable the noninvasive diagnosis of monogenic kidney disease at relatively low cost, thereby allowing for more precise management for patients and their families. A major challenge is interpreting rare variants, many of which are classified as variants of unknown significance (VUS). We present a framework in which we thoroughly evaluated and provided evidence of pathogenicity for HNF1B-p.Arg303His, a VUS returned from clinical diagnostic testing for a kidney transplant candidate. Methods: A blueprint was designed by a multidisciplinary team of clinicians, molecular biologists, and diagnostic geneticists. The blueprint included using a health system-based cohort with genetic and clinical information to perform deep phenotyping of VUS heterozygotes to identify the candidate VUS and rule out other VUS, examination of existing genetic databases, as well as functional testing. Results: Our approach demonstrated evidence for pathogenicity for HNF1B-p.Arg303His by showing similar burden of kidney manifestations in this variant to known HNF1B pathogenic variants, and greater burden compared to noncarriers. Conclusion: Determination of a molecular diagnosis for the example family allows for proper surveillance and management of HNF1B-related manifestations such as kidney disease, diabetes, and hypomagnesemia with important implications for safe living-related kidney donation. The candidate gene-variant pair also allows for clinical biomarker testing for aberrations of linked pathways. This working model may be applicable to other diseases of genetic etiology.

Calcium, Its Regulatory Hormones, and Their Causal Role on Blood Pressure: A Two-Sample Mendelian Randomization Study.

CONTEXT: Vitamin D (Vit-D), parathyroid hormone (PTH), and fibroblast growth factor 23 (FGF23) are the major calciotropic hormones involved in the regulation of blood calcium levels from the intestine, kidney, and bone through a tight endocrine feedback loop system. Altered levels of calcium itself or through the effect of its regulatory hormones could affect blood pressure (BP), but the exact mechanisms remain unclear. OBJECTIVE: To evaluate whether a causal relationship exists between serum calcium level and/or the regulatory hormones involved in its homeostasis with BP, we performed a two-sample Mendelian randomization (MR) study. METHODS: From 4 large genome-wide association studies (GWAS) we obtained independent (r2 < 0.001) single nucleotide polymorphisms (SNPs) associated with serum calcium (119 SNPs), Vit-D (78 SNPs), PTH (5 SNPs), and FGF23 (5 SNPs), to investigate through MR their association with systolic BP (SBP) and diastolic BP (DBP) in a Swedish urban-based study, the Malmö Diet and Cancer study (n = 29 298). Causality was evaluated by the inverse variance weighted method (IVW) and weighted median, while MR Egger and MR-PRESSO were used as sensitivity analyses. RESULTS: Genetically predicted serum calcium level was found to be associated with DBP (IVW: beta = 0.10, SE = 0.04, P = 0.007) and SBP (IVW: beta = 0.07, SE = 0.04, P = 0.04). Genetically predicted Vit-D and PTH showed no association with the traits, while FGF23 was inversely associated with SBP (IVW: beta = -0.11, SE = 0.04, P = 0.01), although this association lost statistical significance in sensitivity analysis. CONCLUSION: Our study shows a direct association between genetically predicted calcium level and DBP, and a weaker association with SBP. No such clear association was found for genetically predicted calciotropic hormone levels. It is of interest to detect which target genes involved in calcium homeostasis mediate the effect of calcium on BP, particularly for improving personalized intervention strategies.

Germline Mutations in CIDEB and Protection against Liver Disease.

BACKGROUND: Exome sequencing in hundreds of thousands of persons may enable the identification of rare protein-coding genetic variants associated with protection from human diseases like liver cirrhosis, providing a strategy for the discovery of new therapeutic targets. METHODS: We performed a multistage exome sequencing and genetic association analysis to identify genes in which rare protein-coding variants were associated with liver phenotypes. We conducted in vitro experiments to further characterize associations. RESULTS: The multistage analysis involved 542,904 persons with available data on liver aminotransferase levels, 24,944 patients with various types of liver disease, and 490,636 controls without liver disease. We found that rare coding variants in APOB, ABCB4, SLC30A10, and TM6SF2 were associated with increased aminotransferase levels and an increased risk of liver disease. We also found that variants in CIDEB, which encodes a structural protein found in hepatic lipid droplets, had a protective effect. The burden of rare predicted loss-of-function variants plus missense variants in CIDEB (combined carrier frequency, 0.7%) was associated with decreased alanine aminotransferase levels (beta per allele, -1.24 U per liter; 95% confidence interval [CI], -1.66 to -0.83; P = 4.8×10-9) and with 33% lower odds of liver disease of any cause (odds ratio per allele, 0.67; 95% CI, 0.57 to 0.79; P = 9.9×10-7). Rare coding variants in CIDEB were associated with a decreased risk of liver disease across different underlying causes and different degrees of severity, including cirrhosis of any cause (odds ratio per allele, 0.50; 95% CI, 0.36 to 0.70). Among 3599 patients who had undergone bariatric surgery, rare coding variants in CIDEB were associated with a decreased nonalcoholic fatty liver disease activity score (beta per allele in score units, -0.98; 95% CI, -1.54 to -0.41 [scores range from 0 to 8, with higher scores indicating more severe disease]). In human hepatoma cell lines challenged with oleate, CIDEB small interfering RNA knockdown prevented the buildup of large lipid droplets. CONCLUSIONS: Rare germline mutations in CIDEB conferred substantial protection from liver disease. (Funded by Regeneron Pharmaceuticals.).

Thrombotic risk determined by rare and common SERPINA1 variants in a population-based cohort study.

BACKGROUND: Severe alpha-1-antitrypsin deficiency (AATD), phenotype PiZZ, was associated with venous thromboembolism (VTE) in a case-control study. OBJECTIVES: This study aimed to determine the genetic variation in the SERPINA1 gene and a possible thrombotic risk of these variants in a population-based cohort study. PATIENTS/METHODS: The coding sequence of SERPINA1 was analyzed for the Z (rs28929474), S (rs17580), and other qualifying variants in 28,794 subjects without previous VTE (born 1923-1950, 60% women), who participated in the Malmö Diet and Cancer study (1991-1996). Individuals were followed from baseline until the first event of VTE, death, or 2018. RESULTS: Resequencing the coding sequence of SERPINA1 identified 84 variants in the total study population, 21 synonymous, 62 missense, and 1 loss-of-function variant. Kaplan-Meier analysis showed that homozygosity for the Z allele increased the risk of VTE whereas heterozygosity showed no effect. The S (rs17580) variant was not associated with VTE. Thirty-one rare variants were qualifying and included in collapsing analysis using the following selection criteria, loss of function, in frame deletion or non-benign (PolyPhen-2) missense variants with minor allele frequency (MAF) <0.1%. Combining the rare qualifying variants with the Z variant showed that carrying two alleles (ZZ or compound heterozygotes) showed increased risk. Cox regression analysis revealed an adjusted hazard ratio of 4.5 (95% confidence interval 2.0-10.0) for combinations of the Z variant and rare qualifying variants. One other variant (rs141620200; MAF = 0.002) showed an increased risk of VTE. CONCLUSIONS: The SERPINA1 ZZ genotype and compound heterozygotes for severe AATD are rare but associated with VTE in a population-based Swedish study.