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.

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

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.

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.

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.

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.

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.

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.

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.

A Protein-Truncating HSD17B13 Variant and Protection from Chronic Liver Disease.

BACKGROUND: Elucidation of the genetic factors underlying chronic liver disease may reveal new therapeutic targets. METHODS: We used exome sequence data and electronic health records from 46,544 participants in the DiscovEHR human genetics study to identify genetic variants associated with serum levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST). Variants that were replicated in three additional cohorts (12,527 persons) were evaluated for association with clinical diagnoses of chronic liver disease in DiscovEHR study participants and two independent cohorts (total of 37,173 persons) and with histopathological severity of liver disease in 2391 human liver samples. RESULTS: A splice variant (rs72613567:TA) in HSD17B13, encoding the hepatic lipid droplet protein hydroxysteroid 17-beta dehydrogenase 13, was associated with reduced levels of ALT (P=4.2×10-12) and AST (P=6.2×10-10). Among DiscovEHR study participants, this variant was associated with a reduced risk of alcoholic liver disease (by 42% [95% confidence interval {CI}, 20 to 58] among heterozygotes and by 53% [95% CI, 3 to 77] among homozygotes), nonalcoholic liver disease (by 17% [95% CI, 8 to 25] among heterozygotes and by 30% [95% CI, 13 to 43] among homozygotes), alcoholic cirrhosis (by 42% [95% CI, 14 to 61] among heterozygotes and by 73% [95% CI, 15 to 91] among homozygotes), and nonalcoholic cirrhosis (by 26% [95% CI, 7 to 40] among heterozygotes and by 49% [95% CI, 15 to 69] among homozygotes). Associations were confirmed in two independent cohorts. The rs72613567:TA variant was associated with a reduced risk of nonalcoholic steatohepatitis, but not steatosis, in human liver samples. The rs72613567:TA variant mitigated liver injury associated with the risk-increasing PNPLA3 p.I148M allele and resulted in an unstable and truncated protein with reduced enzymatic activity. CONCLUSIONS: A loss-of-function variant in HSD17B13 was associated with a reduced risk of chronic liver disease and of progression from steatosis to steatohepatitis. (Funded by Regeneron Pharmaceuticals and others.).

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.

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.

ALG9 Mutation Carriers Develop Kidney and Liver Cysts.

BACKGROUND: Mutations in PKD1 or PKD2 cause typical autosomal dominant polycystic kidney disease (ADPKD), the most common monogenic kidney disease. Dominantly inherited polycystic kidney and liver diseases on the ADPKD spectrum are also caused by mutations in at least six other genes required for protein biogenesis in the endoplasmic reticulum, the loss of which results in defective production of the PKD1 gene product, the membrane protein polycystin-1 (PC1). METHODS: We used whole-exome sequencing in a cohort of 122 patients with genetically unresolved clinical diagnosis of ADPKD or polycystic liver disease to identify a candidate gene, ALG9, and in vitro cell-based assays of PC1 protein maturation to functionally validate it. For further validation, we identified carriers of ALG9 loss-of-function mutations and noncarrier matched controls in a large exome-sequenced population-based cohort and evaluated the occurrence of polycystic phenotypes in both groups. RESULTS: Two patients in the clinically defined cohort had rare loss-of-function variants in ALG9, which encodes a protein required for addition of specific mannose molecules to the assembling N-glycan precursors in the endoplasmic reticulum lumen. In vitro assays showed that inactivation of Alg9 results in impaired maturation and defective glycosylation of PC1. Seven of the eight (88%) cases selected from the population-based cohort based on ALG9 mutation carrier state who had abdominal imaging after age 50; seven (88%) had at least four kidney cysts, compared with none in matched controls without ALG9 mutations. CONCLUSIONS: ALG9 is a novel disease gene in the genetically heterogeneous ADPKD spectrum. This study supports the utility of phenotype characterization in genetically-defined cohorts to validate novel disease genes, and provide much-needed genotype-phenotype correlations.

Trajectory of exonic variant discovery in a large clinical population: implications for variant curation.

PURPOSE: Precision health initiatives and reduced sequencing costs are driving large-scale human genome analyses. Genetic variant curation is a bottleneck in clinical applications. The burden of variant curation can be high for newly discovered variants because they are less likely to have undergone previous clinical annotation; the rate of discovery of genetic variants in large clinical populations has not been empirically determined. METHODS: We determined the rate of accrual of unique sequence variants in 90,000 exome sequences. Separate analyses were done for 17,267 autosomal genes and a subset of 74 actionable genes; the effect of relatedness in the cohort was also determined. RESULTS: Variant discovery showed a nonlinear growth pattern. The rate of unique variant accrual decreased as the database size increased; by 90,000 exomes 97% of all projected coding and splicing variants had been observed. Variants in 74 actionable genes showed a similar pattern. Family relatedness slightly reduced the rate of discovery of unique variants. CONCLUSION: The heaviest burden of interpretation for genetic variants occurs early and diminishes as the database size increases. Our data provide a framework for scaling pathogenic genetic variant discovery and curation, a critical element of patient care in the era of precision health.

Cilia have high cAMP levels that are inhibited by Sonic Hedgehog-regulated calcium dynamics.

Protein kinase A (PKA) phosphorylates Gli proteins, acting as a negative regulator of the Hedgehog pathway. PKA was recently detected within the cilium, and PKA activity specifically in cilia regulates Gli processing. Using a cilia-targeted genetically encoded sensor, we found significant basal PKA activity. Using another targeted sensor, we measured basal ciliary cAMP that is fivefold higher than whole-cell cAMP. The elevated basal ciliary cAMP level is a result of adenylyl cyclase 5 and 6 activity that depends on ciliary phosphatidylinositol (3,4,5)-trisphosphate (PIP3), not stimulatory G protein (Gαs), signaling. Sonic Hedgehog (SHH) reduces ciliary cAMP levels, inhibits ciliary PKA activity, and increases Gli1. Remarkably, SHH regulation of ciliary cAMP and downstream signals is not dependent on inhibitory G protein (Gαi/o) signaling but rather Ca2+ entry through a Gd3+-sensitive channel. Therefore, PIP3 sustains high basal cAMP that maintains PKA activity in cilia and Gli repression. SHH activates Gli by inhibiting cAMP through a G protein-independent mechanism that requires extracellular Ca2+ entry.

CCL20 is up-regulated in non-alcoholic fatty liver disease fibrosis and is produced by hepatic stellate cells in response to fatty acid loading.

BACKGROUND: Nonalcoholic fatty liver disease (NAFLD) is a prevalent complication of extreme obesity. Loading of the liver with fat can progress to inflammation and fibrosis including cirrhosis. The molecular factors involved in the progression from simple steatosis to fibrosis remain poorly understood. METHODS: Gene expression profiling using microarray, PCR array, and RNA sequencing was performed on RNA from liver biopsy tissue from patients with extreme obesity. Patients were grouped based on histological findings including normal liver histology with no steatosis, lobular inflammation, or fibrosis, and grades 1, 2, 3, and 4 fibrosis with coexistent steatosis and lobular inflammation. Validation of expression was conducted using quantitative PCR. Serum analysis was performed using ELISA. Expression analysis of hepatocytes and hepatic stellate cells in response to lipid loading were conducted in vitro using quantitative PCR and ELISA. RESULTS: Three orthogonal methods to profile human liver biopsy RNA each identified the chemokine CCL20 (CC chemokine ligand 20 or MIP-3 alpha) gene as one of the most up-regulated transcripts in NAFLD fibrosis relative to normal histology, validated in a replication group. CCL20 protein levels in serum measured in 224 NAFLD patients were increased in severe fibrosis (p < 0.001), with moderate correlation of hepatic transcript levels and serum levels. Expression of CCL20, but not its cognate receptor CC chemokine receptor 6, was significantly (p < 0.001) increased in response to fatty acid loading in LX-2 hepatic stellate cells, with relative increases greater than those in HepG2 hepatocyte cells. CONCLUSIONS: These results suggest that expression of CCL20, an important inflammatory mediator, is increased in NAFLD fibrosis. CCL20 serves as a chemoattractant molecule for immature dendritic cells, which have been shown to produce many of the inflammatory molecules that mediate liver fibrosis. These data also point to hepatic stellate cells as a key cell type that may respond to lipid loading of the liver.

A conserved mechanism for gating in an ionotropic glutamate receptor.

Ionotropic glutamate receptor (iGluR) channels control synaptic activity. The crystallographic structure of GluA2, the prototypical iGluR, reveals a clamshell-like ligand-binding domain (LBD) that closes in the presence of glutamate to open a gate on the pore lining α-helix. How LBD closure leads to gate opening remains unclear. Here, we show that bending the pore helix at a highly conserved alanine residue (Ala-621) below the gate is responsible for channel opening. Substituting Ala-621 with the smaller more flexible glycine resulted in a basally active, nondesensitizing channel with ∼39-fold increase in glutamate potency without affecting surface expression or binding. On GluA2(A621G), the partial agonist kainate showed efficacy similar to a full agonist, and competitive antagonists CNQX and DNQX acted as a partial agonists. Met-629 in GluA2 sits above the gate and is critical in transmitting LBD closure to the gate. Substituting Met-629 with the flexible glycine resulted in reduced channel activity and glutamate potency. The pore regions in potassium channels are structurally similar to iGluRs. Whereas potassium channels typically use glycines as a hinge for gating, iGluRs use the less flexible alanine as a hinge at a similar position to maintain low basal activity allowing for ligand-mediated gating.

Weight-independent effects of roux-en-Y gastric bypass on glucose homeostasis via melanocortin-4 receptors in mice and humans.

BACKGROUND & AIMS: Roux-en-Y gastric bypass (RYGB) improves glucose homeostasis independently of changes in body weight by unknown mechanisms. Melanocortin-4 receptors (MC4R) have weight-independent effects on glucose homeostasis, via autonomic neurons, and also might contribute to weight loss after RYGB. We investigated whether MC4Rs mediate effects of RYGB, such as its weight-independent effects on glucose homeostasis, in mice and humans. METHODS: We studied C57BL/6 mice with diet-induced obesity, MC4R-deficient mice, and mice that re-express MC4R specifically in autonomic neurons after RYGB or sham surgeries. We also sequenced the MC4R locus in patients undergoing RYGB to investigate diabetes resolution in carriers of rare MC4R variants. RESULTS: MC4Rs in autonomic brainstem neurons (including the parasympathetic dorsal motor vagus) mediated improved glucose homeostasis independent of changes in body weight. In contrast, MC4Rs in cholinergic preganglionic motor neurons (sympathetic and parasympathetic) mediated RYGB-induced increased energy expenditure and weight loss. Increased energy expenditure after RYGB is the predominant mechanism of weight loss and confers resistance to weight gain from a high-fat diet, the effects of which are MC4R-dependent. MC4R-dependent effects of RYGB still occurred in mice with Mc4r haplosufficiency, and early stage diabetes resolved at a similar rate in patients with rare variants of MC4R and noncarriers. However, carriers of MC4R (I251L), a rare variant associated with increased weight loss after RYGB and increased basal activity in vitro, were more likely to have early and weight-independent resolution of diabetes than noncarriers, indicating a role for MC4Rs in the effects of RYGB. CONCLUSIONS: MC4Rs in autonomic neurons mediate beneficial effects of RYGB, including weight-independent improved glucose homeostasis, in mice and humans.