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.

The Phenotypic Spectrum of COL4A3 Heterozygotes.

Introduction: The penetrance and phenotypic spectrum of autosomal dominant Alport Syndrome (ADAS), affecting 1 in 106, remains understudied. Methods: Using data from 174,418 participants in the Geisinger MyCode/DiscovEHR study, an unselected health system-based cohort with whole exome sequencing, we identified 403 participants who were heterozygous for likely pathogenic COL4A3 variants. Phenotypic data was evaluated using International Classification of Diseases (ICD) codes, laboratory data, and chart review. To evaluate the phenotypic spectrum of genetically-determined ADAS, we matched COL4A3 heterozygotes 1:5 to nonheterozygotes using propensity scores by demographics, hypertension, diabetes, and nephrolithiasis. Results: COL4A3 heterozygotes were at significantly increased risks of hematuria, decreased estimated glomerular filtration rate (eGFR), albuminuria, and kidney failure (P < 0.05 for all comparisons) but not bilateral sensorineural hearing loss (P = 0.9). Phenotypic severity was more severe for collagenous domain glycine missense variants than protein truncating variants (PTVs). For example, patients with Gly695Arg (n = 161) had markedly increased risk of dipstick hematuria (odds ratio [OR] 9.50; 95% confidence interval [CI]: 6.32, 14.28) and kidney failure (OR 7.02; 95% CI: 3.48, 14.16) whereas those with PTVs (n = 119) had moderately increased risks of dipstick hematuria (OR 1.64; 95% CI: 1.03, 2.59) and kidney failure (OR 3.44; 95% CI: 1.28, 9.22). Less than a third of patients had albuminuria screening completed, and fewer than 1 of 3 were taking inhibitors of the renin-angiotensin-aldosterone system. Conclusion: This study demonstrates a wide spectrum of phenotypic severity in ADAS due to COL4A3 with phenotypic variability by genotype. Future studies are needed to evaluate the impact of earlier diagnosis, appropriate evaluation, and treatment of ADAS.

GWAS for the composite traits of hematuria and albuminuria.

Our GWAS of hematuria in the UK Biobank identified 6 loci, some of which overlap with loci for albuminuria suggesting pleiotropy. Since clinical syndromes are often defined by combinations of traits, generating a combined phenotype can improve power to detect loci influencing multiple characteristics. Thus the composite trait of hematuria and albuminuria was chosen to enrich for glomerular pathologies. Cases had both hematuria defined by ICD codes and albuminuria defined as uACR > 3 mg/mmol. Controls had neither an ICD code for hematuria nor an uACR > 3 mg/mmol. 2429 cases and 343,509 controls from the UK Biobank were included. eGFR was lower in cases compared to controls, with the exception of the comparison in females using CKD-EPI after age adjustment. Variants at 4 loci met genome-wide significance with the following nearest genes: COL4A4, TRIM27, ETV1 and CUBN. TRIM27 is part of the extended MHC locus. All loci with the exception of ETV1 were replicated in the Geisinger MyCode cohort. The previous GWAS of hematuria reported COL4A3-COL4A4 variants and HLA-B*0801 within MHC, which is in linkage disequilibrium with the TRIM27 variant (D' = 0.59). TRIM27 is highly expressed in the tubules. Additional loci included a coding sequence variant in CUBN (p.Ala2914Val, MAF = 0.014 (A), p = 3.29E-8, OR = 2.09, 95% CI = 1.61-2.72). Overall, GWAS for the composite trait of hematuria and albuminuria identified 4 loci, 2 of which were not previously identified in a GWAS of hematuria.

The Phenotypic Spectrum of COL4A3 Heterozygotes.

Most data on Alport Syndrome (AS) due to COL4A3 are limited to families with autosomal recessive AS or severe manifestations such as focal segmental glomerulosclerosis (FSGS). Using data from 174,418 participants in the Geisinger MyCode/DiscovEHR study, an unselected health system-based cohort with whole exome sequencing, we identified 403 participants (0.2%) who were heterozygous for likely pathogenic COL4A3 variants. Phenotypic data was evaluated using International Classification of Diseases (ICD) codes, laboratory data, and chart review. To evaluate the phenotypic spectrum of genetically-determined autosomal dominant AS, we matched COL4A3 heterozygotes 1:5 to non-heterozygotes using propensity scores by demographics, hypertension, diabetes, and nephrolithiasis. COL4A3 heterozygotes were at significantly increased risks of hematuria, decreased estimated glomerular filtration rate (eGFR), albuminuria, and end-stage kidney disease (ESKD) (p<0.05 for all comparisons) but not bilateral sensorineural hearing loss (p=0.9). Phenotypic severity tended to be more severe among patients with glycine missense variants located within the collagenous domain. For example, patients with Gly695Arg (n=161) had markedly increased risk of dipstick hematuria (OR 9.47, 95% CI: 6.30, 14.22) and ESKD diagnosis (OR 7.01, 95% CI: 3.48, 14.12) whereas those with PTVs (n=119) had moderately increased risks of dipstick hematuria (OR 1.63, 95% CI: 1.03, 2.58) and ESKD diagnosis (OR 3.43, 95% CI: 1.28, 9.19). Less than a third of patients had albuminuria screening completed, and fewer than 1/3 were taking inhibitors of the renin-angiotensin-aldosterone system (RAASi). Future studies are needed to evaluate the impact of earlier diagnosis, appropriate evaluation, and treatment of ADAS.

Prediction of Kv11.1 potassium channel PAS-domain variants trafficking via machine learning.

Congenital long QT syndrome (LQTS) is characterized by a prolonged QT-interval on an electrocardiogram (ECG). An abnormal prolongation in the QT-interval increases the risk for fatal arrhythmias. Genetic variants in several different cardiac ion channel genes, including KCNH2, are known to cause LQTS. Here, we evaluated whether structure-based molecular dynamics (MD) simulations and machine learning (ML) could improve the identification of missense variants in LQTS-linked genes. To do this, we investigated KCNH2 missense variants in the Kv11.1 channel protein shown to have wild type (WT) like or class II (trafficking-deficient) phenotypes in vitro. We focused on KCNH2 missense variants that disrupt normal Kv11.1 channel protein trafficking, as it is the most common phenotype for LQTS-associated variants. Specifically, we used computational techniques to correlate structural and dynamic changes in the Kv11.1 channel protein PAS domain (PASD) with Kv11.1 channel protein trafficking phenotypes. These simulations unveiled several molecular features, including the numbers of hydrating waters and hydrogen bonding pairs, as well as folding free energy scores, that are predictive of trafficking. We then used statistical and machine learning (ML) (Decision tree (DT), Random forest (RF), and Support vector machine (SVM)) techniques to classify variants using these simulation-derived features. Together with bioinformatics data, such as sequence conservation and folding energies, we were able to predict with reasonable accuracy (≈75%) which KCNH2 variants do not traffic normally. We conclude that structure-based simulations of KCNH2 variants localized to the Kv11.1 channel PASD led to an improvement in classification accuracy. Therefore, this approach should be considered to complement the classification of variant of unknown significance (VUS) in the Kv11.1 channel PASD.

Individuals heterozygous for ALG8 protein-truncating variants are at increased risk of a mild cystic kidney disease.

ALG8 protein-truncating variants (PTVs) have previously been described in patients with polycystic liver disease and in some cases cystic kidney disease. Given a lack of well-controlled studies, we determined whether individuals heterozygous for ALG8 PTVs are at increased risk of cystic kidney disease in a large, unselected health system-based observational cohort linked to electronic health records in Pennsylvania (Geisinger-Regeneron DiscovEHR MyCode study). Out of 174,172 patients, 236 were identified with ALG8 PTVs. Using ICD-based outcomes, patients with these variants were significantly at increased risk of having any kidney/liver cyst diagnosis (Odds Ratio 2.42, 95% confidence interval: 1.53-3.85), cystic kidney disease (3.03, 1.26-7.31), and nephrolithiasis (1.89, 1.96-2.97). To confirm this finding, blinded radiology review of computed tomography and magnetic resonance imaging studies was completed in a matched cohort of 52 thirty-plus year old ALG8 PTV heterozygotes and related non-heterozygotes. ALG8 PTV heterozygotes were significantly more likely to have cystic kidney disease, defined as four or more kidney cysts (57.7% vs. 7.7%), or bilateral kidney cysts (69.2% vs. 15.4%), but not one or more liver cyst (11.5% vs. 7.7%). In publicly available UK Biobank data, ALG8 PTV heterozygotes were at significantly increased risk of ICD code N28 (other disorders of kidney/ureter) (3.85% vs. 1.33%). ALG8 PTVs were not associated with chronic kidney disease or kidney failure in the MyCode study or the UK Biobank data. Thus, PTVs in ALG8 result in increased risk of a mild cystic kidney disease phenotype.

Exome Sequencing of a Clinical Population for Autosomal Dominant Polycystic Kidney Disease.

Importance: Most studies of autosomal dominant polycystic kidney disease (ADPKD) genetics have used kidney specialty cohorts, focusing on PKD1 and PKD2. These can lead to biased estimates of population prevalence of ADPKD-associated gene variants and their phenotypic expression. Objective: To determine the prevalence of ADPKD and contributions of PKD1, PKD2, and other genes related to cystic kidney disease in a large, unselected cohort. Design, Setting, and Participants: This retrospective observational study used an unselected health system-based cohort in central and northeast Pennsylvania with exome sequencing (enrolled from 2004 to 2020) and electronic health record data (up to October 2021). The genotype-first approach included the entire cohort and the phenotype-first approach focused on patients with ADPKD diagnosis codes, confirmed by chart and imaging review. Exposures: Loss-of-function (LOF) variants in PKD1, PKD2, and other genes associated with cystic kidney disease (ie, ALG8, ALG9, DNAJB11, GANAB, HNF1B, IFT140, SEC61B, PKHD1, PRKCSH, SEC63); likely pathogenic missense variants in PKD1 and PKD2. Main Outcomes and Measures: Genotype-first analysis: ADPKD diagnosis code (Q61.2, Q61.3, 753.13, 753.12); phenotype-first analysis: presence of a rare variant in PKD1, PKD2, or other genes associated with cystic kidney disease. Results: Of 174 172 patients (median age, 60 years; 60.6% female; 93% of European ancestry), 303 patients had ADPKD diagnosis codes, including 235 with sufficient chart review data for confirmation. In addition to PKD1 and PKD2, LOF variants in IFT140, GANAB, and HNF1B were associated with ADPKD diagnosis after correction for multiple comparisons. Among patients with LOF variants in PKD1, 66 of 68 (97%) had ADPKD; 43 of 43 patients (100%) with LOF variants in PKD2 had ADPKD. In contrast, only 24 of 77 patients (31.2%) with a PKD1 missense variant previously classified as "likely pathogenic" had ADPKD, suggesting misclassification or variable penetrance. Among patients with ADPKD diagnosis confirmed by chart review, 180 of 235 (76.6%) had a potential genetic cause, with the majority being rare variants in PKD1 (127 patients) or PKD2 (34 patients); 19 of 235 (8.1%) had variants in other genes associated with cystic kidney disease. Of these 235 patients with confirmed ADPKD, 150 (63.8%) had a family history of ADPKD. The yield for a genetic determinant of ADPKD was higher for those with a family history of ADPKD compared with those without family history (91.3% [137/150] vs 50.6% [43/85]; difference, 40.7% [95% CI, 29.2%-52.3%]; P < .001). Previously unreported PKD1, PKD2, and GANAB variants were identified with pedigree data suggesting pathogenicity, and several PKD1 missense variants previously reported as likely pathogenic appeared to be benign. Conclusions and Relevance: This study demonstrates substantial genetic and phenotypic variability in ADPKD among patients within a regional health system in the US.

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.).

Genome-wide analysis provides genetic evidence that ACE2 influences COVID-19 risk and yields risk scores associated with severe disease.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) enters human host cells via angiotensin-converting enzyme 2 (ACE2) and causes coronavirus disease 2019 (COVID-19). Here, through a genome-wide association study, we identify a variant (rs190509934, minor allele frequency 0.2-2%) that downregulates ACE2 expression by 37% (P = 2.7 × 10-8) and reduces the risk of SARS-CoV-2 infection by 40% (odds ratio = 0.60, P = 4.5 × 10-13), providing human genetic evidence that ACE2 expression levels influence COVID-19 risk. We also replicate the associations of six previously reported risk variants, of which four were further associated with worse outcomes in individuals infected with the virus (in/near LZTFL1, MHC, DPP9 and IFNAR2). Lastly, we show that common variants define a risk score that is strongly associated with severe disease among cases and modestly improves the prediction of disease severity relative to demographic and clinical factors alone.

Whole-genome sequencing reveals host factors underlying critical COVID-19.

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2-4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes-including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)-in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease.

Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study.

BACKGROUND: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. METHODS: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record-derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. RESULTS: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. CONCLUSIONS: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier; NCT03394859.

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

Large-scale human exome sequencing can identify rare protein-coding variants with a large impact on complex traits such as body adiposity. We sequenced the exomes of 645,626 individuals from the United Kingdom, the United States, and Mexico and estimated associations of rare coding variants with body mass index (BMI). We identified 16 genes with an exome-wide significant association with BMI, including those encoding five brain-expressed G protein-coupled receptors (CALCR, MC4R, GIPR, GPR151, and GPR75). Protein-truncating variants in GPR75 were observed in ~4/10,000 sequenced individuals and were associated with 1.8 kilograms per square meter lower BMI and 54% lower odds of obesity in the heterozygous state. Knock out of Gpr75 in mice resulted in resistance to weight gain and improved glycemic control in a high-fat diet model. Inhibition of GPR75 may provide a therapeutic strategy for obesity.

Large scale clinical exome sequencing uncovers the scope and severity of skin disorders associated with MC1R genetic variants.

PURPOSE: Genetic variation in MC1R is a main determinant of red hair color (RHC) phenotype and confers susceptibility to skin disorders. METHODS: We assessed the effects and function of MC1R variants identified in our clinical cohort of 135,947 participants with available exome sequencing using phenome-wide association scan (PheWAS). Expression and function of several variants were evaluated. RESULTS: We found 24 nonsense and 215 missense variants in MC1R. Many common missense MC1R variants are strongly associated with skin disorders including skin cancer; however, each variant shows different penetrance and expressivity. Severity of skin phenotype was well correlated with the magnitude of functional defect measured as receptor expression and α-MSH stimulated cAMP production. Remarkably, MC1R deletions and nonsense variants are only weakly associated with milder skin phenotypes. CONCLUSION: Our comprehensive assessment of all MC1R variants in a large cohort clearly establish that individuals with some missense variants are more susceptible to severe skin disorders than those with MC1R deletions or nonsense variants.

Rare Coding Variants Associated With Electrocardiographic Intervals Identify Monogenic Arrhythmia Susceptibility Genes: A Multi-Ancestry Analysis.

BACKGROUND: Alterations in electrocardiographic (ECG) intervals are well-known markers for arrhythmia and sudden cardiac death (SCD) risk. While the genetics of arrhythmia syndromes have been studied, relations between electrocardiographic intervals and rare genetic variation at a population level are poorly understood. METHODS: Using a discovery sample of 29 000 individuals with whole-genome sequencing from Trans-Omics in Precision Medicine and replication in nearly 100 000 with whole-exome sequencing from the UK Biobank and MyCode, we examined associations between low-frequency and rare coding variants with 5 routinely measured electrocardiographic traits (RR, P-wave, PR, and QRS intervals and corrected QT interval). RESULTS: We found that rare variants associated with population-based electrocardiographic intervals identify established monogenic SCD genes (KCNQ1, KCNH2, and SCN5A), a controversial monogenic SCD gene (KCNE1), and novel genes (PAM and MFGE8) involved in cardiac conduction. Loss-of-function and pathogenic SCN5A variants, carried by 0.1% of individuals, were associated with a nearly 6-fold increased odds of the first-degree atrioventricular block (P=8.4×10-5). Similar variants in KCNQ1 and KCNH2 (0.2% of individuals) were associated with a 23-fold increased odds of marked corrected QT interval prolongation (P=4×10-25), a marker of SCD risk. Incomplete penetrance of such deleterious variation was common as over 70% of carriers had normal electrocardiographic intervals. CONCLUSIONS: Our findings indicate that large-scale high-depth sequence data and electrocardiographic analysis identifies monogenic arrhythmia susceptibility genes and rare variants with large effects. Known pathogenic variation in conventional arrhythmia and SCD genes exhibited incomplete penetrance and accounted for only a small fraction of marked electrocardiographic interval prolongation.

Pan-ancestry exome-wide association analyses of COVID-19 outcomes in 586,157 individuals.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19), a respiratory illness that can result in hospitalization or death. We used exome sequence data to investigate associations between rare genetic variants and seven COVID-19 outcomes in 586,157 individuals, including 20,952 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome wide or when specifically focusing on (1) 13 interferon pathway genes in which rare deleterious variants have been reported in individuals with severe COVID-19, (2) 281 genes located in susceptibility loci identified by the COVID-19 Host Genetics Initiative, or (3) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, and results are publicly available through the Regeneron Genetics Center COVID-19 Results Browser.

Genome-wide association analysis of serum alanine and aspartate aminotransferase, and the modifying effects of BMI in 388k European individuals.

Serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) are biomarkers for liver health. Here we report the largest genome-wide association analysis to date of serum ALT and AST levels in over 388k people of European ancestry from UK biobank and DiscovEHR. Eleven million imputed markers with a minor allele frequency (MAF) ≥ 0.5% were analyzed. Overall, 300 ALT and 336 AST independent genome-wide significant associations were identified. Among them, 81 ALT and 61 AST associations are reported for the first time. Genome-wide interaction study identified 9 ALT and 12 AST independent associations significantly modified by body mass index (BMI), including several previously reported potential liver disease therapeutic targets, for example, PNPLA3, HSD17B13, and MARC1. While further work is necessary to understand the effect of ALT and AST-associated variants on liver disease, the weighted burden of significant BMI-modified signals is significantly associated with liver disease outcomes. In summary, this study identifies genetic associations which offer an important step forward in understanding the genetic architecture of serum ALT and AST levels. Significant interactions between BMI and genetic loci not only highlight the important role of adiposity in liver damage but also shed light on the genetic etiology of liver disease in obese individuals.

Evaluation of the MC4R gene across eMERGE network identifies many unreported obesity-associated variants.

BACKGROUND/OBJECTIVES: Melanocortin-4 receptor (MC4R) plays an essential role in food intake and energy homeostasis. More than 170 MC4R variants have been described over the past two decades, with conflicting reports regarding the prevalence and phenotypic effects of these variants in diverse cohorts. To determine the frequency of MC4R variants in large cohort of different ancestries, we evaluated the MC4R coding region for 20,537 eMERGE participants with sequencing data plus additional 77,454 independent individuals with genome-wide genotyping data at this locus. SUBJECTS/METHODS: The sequencing data were obtained from the eMERGE phase III study, in which multisample variant call format calls have been generated, curated, and annotated. In addition to penetrance estimation using body mass index (BMI) as a binary outcome, GWAS and PheWAS were performed using median BMI in linear regression analyses. All results were adjusted for principal components, age, sex, and sites of genotyping. RESULTS: Targeted sequencing data of MC4R revealed 125 coding variants in 1839 eMERGE participants including 30 unreported coding variants that were predicted to be functionally damaging. Highly penetrant unreported variants included (L325I, E308K, D298N, S270F, F261L, T248A, D111V, and Y80F) in which seven participants had obesity class III defined as BMI ≥ 40 kg/m2. In GWAS analysis, in addition to known risk haplotype upstream of MC4R (best variant rs6567160 (P = 5.36 × 10-25, Beta = 0.37), a novel rare haplotype was detected which was protective against obesity and encompassed the V103I variant with known gain-of-function properties (P = 6.23 × 10-08, Beta = -0.62). PheWAS analyses extended this protective effect of V103I to type 2 diabetes, diabetic nephropathy, and chronic renal failure independent of BMI. CONCLUSIONS: MC4R screening in a large eMERGE cohort confirmed many previous findings, extend the MC4R pleotropic effects, and discovered additional MC4R rare alleles that probably contribute to obesity.

Electronic health record analysis identifies kidney disease as the leading risk factor for hospitalization in confirmed COVID-19 patients.

BACKGROUND: Empirical data on conditions that increase risk of coronavirus disease 2019 (COVID-19) progression are needed to identify high risk individuals. We performed a comprehensive quantitative assessment of pre-existing clinical phenotypes associated with COVID-19-related hospitalization. METHODS: Phenome-wide association study (PheWAS) of SARS-CoV-2-positive patients from an integrated health system (Geisinger) with system-level outpatient/inpatient COVID-19 testing capacity and retrospective electronic health record (EHR) data to assess pre-COVID-19 pandemic clinical phenotypes associated with hospital admission (hospitalization). RESULTS: Of 12,971 individuals tested for SARS-CoV-2 with sufficient pre-COVID-19 pandemic EHR data at Geisinger, 1604 were SARS-CoV-2 positive and 354 required hospitalization. We identified 21 clinical phenotypes in 5 disease categories meeting phenome-wide significance (P<1.60x10-4), including: six kidney phenotypes, e.g. end stage renal disease or stage 5 CKD (OR = 11.07, p = 1.96x10-8), six cardiovascular phenotypes, e.g. congestive heart failure (OR = 3.8, p = 3.24x10-5), five respiratory phenotypes, e.g. chronic airway obstruction (OR = 2.54, p = 3.71x10-5), and three metabolic phenotypes, e.g. type 2 diabetes (OR = 1.80, p = 7.51x10-5). Additional analyses defining CKD based on estimated glomerular filtration rate, confirmed high risk of hospitalization associated with pre-existing stage 4 CKD (OR 2.90, 95% CI: 1.47, 5.74), stage 5 CKD/dialysis (OR 8.83, 95% CI: 2.76, 28.27), and kidney transplant (OR 14.98, 95% CI: 2.77, 80.8) but not stage 3 CKD (OR 1.03, 95% CI: 0.71, 1.48). CONCLUSIONS: This study provides quantitative estimates of the contribution of pre-existing clinical phenotypes to COVID-19 hospitalization and highlights kidney disorders as the strongest factors associated with hospitalization in an integrated US healthcare system.

Gi-coupled receptor activation potentiates Piezo2 currents via Gßγ.

Mechanically activated Piezo2 channels are key players in somatosensory touch, but their regulation by cellular signaling pathways is poorly understood. Dorsal root ganglion (DRG) neurons express a variety of G-protein-coupled receptors that modulate the function of sensory ion channels. Gi-coupled receptors are generally considered inhibitory, as they usually decrease excitability. Paradoxically, activation of Gi-coupled receptors in DRG neurons sometimes induces mechanical hypersensitivity, the mechanism of which is not well understood. Here, we find that activation of Gi-coupled receptors potentiates mechanically activated currents in DRG neurons and heterologously expressed Piezo2 channels, but inhibits Piezo1 currents in heterologous systems in a Gßγ-dependent manner. Pharmacological inhibition of kinases downstream of Gßγ, phosphoinositide 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) also abolishes the potentiation of Piezo2 currents. Local injection of sumatriptan, an agonist of the Gi-coupled serotonin 1B/1D receptors, increases mechanical sensitivity in mice, and the effect is abolished by inhibiting PI3K and MAPK. Hence, our studies illustrate an indirect mechanism of action of Gßγ to sensitize Piezo2 currents and alter mechanosensitivity after activation of Gi-coupled receptors.