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.

Loss of Cardiac Splicing Regulator RBM20 Is Associated With Early-Onset Atrial Fibrillation.

We showed an association between atrial fibrillation and rare loss-of-function (LOF) variants in the cardiac splicing regulator RBM20 in 2 independent cohorts. In a rat model with loss of RBM20, we demonstrated altered splicing of sarcomere genes (NEXN, TTN, TPM1, MYOM1, and LDB3), and differential expression in key cardiac genes. We identified altered sarcomere and mitochondrial structure on electron microscopy imaging and found compromised mitochondrial function. Finally, we demonstrated that 3 novel LOF variants in RBM20, identified in patients with atrial fibrillation, lead to significantly reduced splicing activity. Our results implicate alternative splicing as a novel proarrhythmic mechanism in the atria.

LiverZap: a chemoptogenetic tool for global and locally restricted hepatocyte ablation to study cellular behaviours in liver regeneration.

The liver restores its mass and architecture after injury. Yet, investigating morphogenetic cell behaviours and signals that repair tissue architecture at high spatiotemporal resolution remains challenging. We developed LiverZap, a tuneable chemoptogenetic liver injury model in zebrafish. LiverZap employs the formation of a binary FAP-TAP photosensitiser followed by brief near-infrared illumination inducing hepatocyte-specific death and recapitulating mammalian liver injury types. The tool enables local hepatocyte ablation and extended live imaging capturing regenerative cell behaviours, which is crucial for studying cellular interactions at the interface of healthy and damaged tissue. Applying LiverZap, we show that targeted hepatocyte ablation in a small region of interest is sufficient to trigger local liver progenitor-like cell (LPC)-mediated regeneration, challenging the current understanding of liver regeneration. Surprisingly, the LPC response is also elicited in adjacent uninjured tissue, at up to 100 µm distance to the injury. Moreover, dynamic biliary network rearrangement suggests active cell movements from uninjured tissue in response to substantial hepatocyte loss as an integral step of LPC-mediated liver regeneration. This precisely targetable liver cell ablation tool will enable the discovery of key molecular and morphogenetic regeneration paradigms.

Lineage tracing identifies heterogeneous hepatoblast contribution to cell lineages and postembryonic organ growth dynamics.

To meet the physiological demands of the body, organs need to establish a functional tissue architecture and adequate size as the embryo develops to adulthood. In the liver, uni- and bipotent progenitor differentiation into hepatocytes and biliary epithelial cells (BECs), and their relative proportions, comprise the functional architecture. Yet, the contribution of individual liver progenitors at the organ level to both fates, and their specific proportion, is unresolved. Combining mathematical modelling with organ-wide, multispectral FRaeppli-NLS lineage tracing in zebrafish, we demonstrate that a precise BEC-to-hepatocyte ratio is established (i) fast, (ii) solely by heterogeneous lineage decisions from uni- and bipotent progenitors, and (iii) independent of subsequent cell type-specific proliferation. Extending lineage tracing to adulthood determined that embryonic cells undergo spatially heterogeneous three-dimensional growth associated with distinct environments. Strikingly, giant clusters comprising almost half a ventral lobe suggest lobe-specific dominant-like growth behaviours. We show substantial hepatocyte polyploidy in juveniles representing another hallmark of postembryonic liver growth. Our findings uncover heterogeneous progenitor contributions to tissue architecture-defining cell type proportions and postembryonic organ growth as key mechanisms forming the adult liver.

Outlining cardiac ion channel protein interactors and their signature in the human electrocardiogram.

Protein-protein interactions are essential for normal cellular processes and signaling events. Defining these interaction networks is therefore crucial for understanding complex cellular functions and interpretation of disease-associated gene variants. We need to build a comprehensive picture of the interactions, their affinities and interdependencies in the specific organ to decipher hitherto poorly understood signaling mechanisms through ion channels. Here we report the experimental identification of the ensemble of protein interactors for 13 types of ion channels in murine cardiac tissue. Of these, we validated the functional importance of ten interactors on cardiac electrophysiology through genetic knockouts in zebrafish, gene silencing in mice, super-resolution microscopy and patch clamp experiments. Furthermore, we establish a computational framework to reconstruct human cardiomyocyte ion channel networks from deep proteome mapping of human heart tissue and human heart single-cell gene expression data. Finally, we integrate the ion channel interactome with human population genetics data to identify proteins that influence the electrocardiogram (ECG). We demonstrate that the combined channel network is enriched for proteins influencing the ECG, with 44% of the network proteins significantly associated with an ECG phenotype. Altogether, we define interactomes of 13 major cardiac ion channels, contextualize their relevance to human electrophysiology and validate functional roles of ten interactors, including two regulators of the sodium current (epsin-2 and gelsolin). Overall, our data provide a roadmap for our understanding of the molecular machinery that regulates cardiac electrophysiology.

Polygenic risk score for ACE-inhibitor-associated cough based on the discovery of new genetic loci.

AIMS: To search for sequence variants associated with ACEi discontinuation and to test their association with ACEi-associated adverse drug reactions (ADRs). METHODS AND RESULTS: A genome-wide association study (GWAS) on ACEi discontinuation was conducted, including 33 959 ACEi-discontinuers and 44 041 controls. Cases were defined as persons who switched from an ACEi treatment to an angiotensin receptor blocker. Controls were defined as persons who continued ACEi treatment for at least 1 year. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were computed for ACEi discontinuation risk by mixed model regression analysis. Summary statistics from the individual cohorts were meta-analyzed with a fixed-effects model. To test for association with specific ACEi-associated ADRs, any genome-wide significant (P < 5 × 10-8) ACEi discontinuation variants was tested for association with ACEi-associated cough and angioedema. A polygenetic risk score (PRS) based on ACEi discontinuation GWAS data was constructed and tested for association with ACEi-associated cough and angioedema in two population-based samples. In total, seven genetic genome-wide loci were identified, of which six were previously unreported. The strongest association with ACEi discontinuation was at 20q13.3 (NTSR1; OR: 1.21; 95% CI: 1.17-1.24; P = 2.1 × 10-34). Five of seven lead variants were associated with ACEi-associated cough, whereas none were associated with ACEi-associated angioedema. The ACEi discontinuation PRS was associated with ACEi-associated cough in a dose-response manner but not with ACEi-associated angioedema. ACEi discontinuation was genetically correlated with important causes for cough, including gastro-esophageal reflux disease, allergic rhinitis, hay fever, and asthma, which indicates partly shared genetic underpinning between these traits. CONCLUSION: This study showed the advantage of using prescription patterns to discover genetic links with ADRs. In total, seven genetic loci that associated with ACEi discontinuation were identified. There was evidence of a strong association between our ADR phenotype and ACEi-associated cough. Taken together, these findings increase insight into the pathophysiological processes that underlie ACEi-associated ADRs.

A Multi-Site Assessment of Anesthetic Overdose, Hypothermic Shock, and Electrical Stunning as Methods of Euthanasia for Zebrafish (Danio rerio) Embryos and Larvae.

Euthanasia in zebrafish (Danio rerio) younger than 5 days post fertilization (dpf) is poorly described in the literature, and standardized protocols are lacking, most likely because larvae not capable of independent feeding are often not protected under national legislations. We assessed the euthanasia efficacy in laboratories in different countries of a one hour anesthetic overdose immersion with buffered lidocaine hydrochloride (1 g/L, with or without 50 mL/L of ethanol), buffered tricaine (1 g/L), clove oil (0.1%), benzocaine (1 g/L), or 2-phenoxyethanol (3 mL/L), as well as the efficacy of hypothermic shock (one hour immersion) and electrical stunning (for one minute), on zebrafish at <12 h post fertilization (hpf), 24 hpf, and 4 dpf. Based on the survival/recovery rates 24 h after treatment, the most effective methods were clove oil, lidocaine with ethanol, and electrical stunning. For 4 dpf larvae, signs of aversion during treatment demonstrated that all anesthetics, except lidocaine, induced aversive behavior. Therefore, the most suited euthanasic treatment was lidocaine hydrochloride 1 g/L, buffered with 2 g/L of sodium bicarbonate and mixed with 50 mL/L of ethanol, which euthanized both embryos and larvae in an efficient and stress-free manner. Electrical stunning also euthanized embryos and larvae efficiently and without signs of aversion; this method needs further assessment in other laboratories to draw firm conclusions.

Integration of multiple imaging platforms to uncover cardiovascular defects in adult zebrafish.

AIMS: Mammalian models have been instrumental in investigating adult heart function and human disease. However, electrophysiological differences with human hearts and high costs motivate the need for non-mammalian models. The zebrafish is a well-established genetic model to study cardiovascular development and function; however, analysis of cardiovascular phenotypes in adult specimens is particularly challenging as they are opaque. METHODS AND RESULTS: Here, we optimized and combined multiple imaging techniques including echocardiography, magnetic resonance imaging, and micro-computed tomography to identify and analyse cardiovascular phenotypes in adult zebrafish. Using alk5a/tgfbr1a mutants as a case study, we observed morphological and functional cardiovascular defects that were undetected with conventional approaches. Correlation analysis of multiple parameters revealed an association between haemodynamic defects and structural alterations of the heart, as observed clinically. CONCLUSION: We report a new, comprehensive, and sensitive platform to identify otherwise indiscernible cardiovascular phenotypes in adult zebrafish.

Quantitative proteome comparison of human hearts with those of model organisms.

Delineating human cardiac pathologies and their basic molecular mechanisms relies on research conducted in model organisms. Yet translating findings from preclinical models to humans present a significant challenge, in part due to differences in cardiac protein expression between humans and model organisms. Proteins immediately determine cellular function, yet their large-scale investigation in hearts has lagged behind those of genes and transcripts. Here, we set out to bridge this knowledge gap: By analyzing protein profiles in humans and commonly used model organisms across cardiac chambers, we determine their commonalities and regional differences. We analyzed cardiac tissue from each chamber of human, pig, horse, rat, mouse, and zebrafish in biological replicates. Using mass spectrometry-based proteomics workflows, we measured and evaluated the abundance of approximately 7,000 proteins in each species. The resulting knowledgebase of cardiac protein signatures is accessible through an online database: atlas.cardiacproteomics.com. Our combined analysis allows for quantitative evaluation of protein abundances across cardiac chambers, as well as comparisons of cardiac protein profiles across model organisms. Up to a quarter of proteins with differential abundances between atria and ventricles showed opposite chamber-specific enrichment between species; these included numerous proteins implicated in cardiac disease. The generated proteomics resource facilitates translational prospects of cardiac studies from model organisms to humans by comparisons of disease-linked protein networks across species.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease.

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Early-onset atrial fibrillation patients show reduced left ventricular ejection fraction and increased atrial fibrosis.

Atrial fibrillation (AF) has traditionally been considered an electrical heart disease. However, genetic studies have revealed that the structural architecture of the heart also play a significant role. We evaluated the functional and structural consequences of harboring a titin-truncating variant (TTNtv) in AF patients, using cardiac magnetic resonance (CMR). Seventeen early-onset AF cases carrying a TTNtv, were matched 1:1 with non-AF controls and a replication cohort of early-onset AF cases without TTNtv, and underwent CMR. Cardiac volumes and left atrial late gadolinium enhancement (LA LGE), as a fibrosis proxy, were measured by a blinded operator. Results: AF cases with TTNtv had significantly reduced left ventricular ejection fraction (LVEF) compared with controls (57 ± 4 vs 64 ± 5%, P < 0.001). We obtained similar findings in early-onset AF patients without TTNtv compared with controls (61 ± 4 vs 64 ± 5%, P = 0.02). We furthermore found a statistically significant increase in LA LGE when comparing early-onset AF TTNtv cases with controls. Using state-of-the-art CMR, we found that early-onset AF patients, irrespective of TTNtv carrier status, had reduced LVEF, indicating that early-onset AF might not be as benign as previously thought.

Early sarcomere and metabolic defects in a zebrafish pitx2c cardiac arrhythmia model.

Atrial fibrillation (AF) is the most common type of cardiac arrhythmia. The major AF susceptibility locus 4q25 establishes long-range interactions with the promoter of PITX2, a transcription factor gene with critical functions during cardiac development. While many AF-linked loci have been identified in genome-wide association studies, mechanistic understanding into how genetic variants, including those at the 4q25 locus, increase vulnerability to AF is mostly lacking. Here, we show that loss of pitx2c in zebrafish leads to adult cardiac phenotypes with substantial similarities to pathologies observed in AF patients, including arrhythmia, atrial conduction defects, sarcomere disassembly, and altered cardiac metabolism. These phenotypes are also observed in a subset of pitx2c+/- fish, mimicking the situation in humans. Most notably, the onset of these phenotypes occurs at an early developmental stage. Detailed analyses of pitx2c loss- and gain-of-function embryonic hearts reveal changes in sarcomeric and metabolic gene expression and function that precede the onset of cardiac arrhythmia first observed at larval stages. We further find that antioxidant treatment of pitx2c-/- larvae significantly reduces the incidence and severity of cardiac arrhythmia, suggesting that metabolic dysfunction is an important driver of conduction defects. We propose that these early sarcomere and metabolic defects alter cardiac function and contribute to the electrical instability and structural remodeling observed in adult fish. Overall, these data provide insight into the mechanisms underlying the development and pathophysiology of some cardiac arrhythmias and importantly, increase our understanding of how developmental perturbations can predispose to functional defects in the adult heart.

Reappraisal of variants previously linked with sudden infant death syndrome: results from three population-based cohorts.

We aimed to investigate the pathogenicity of cardiac ion channel variants previously associated with SIDS. We reviewed SIDS-associated variants previously reported in databases and the literature in three large population-based cohorts; The ExAC database, the Inter99 study, and the UK Biobank (UKBB). Variants were classified according to the American College of Medical Genetics and Genomics (ACMG) guidelines. Of the 92 SIDS-associated variants, 59 (64%) were present in ExAC, 18 (20%) in Inter99, and 24 (26%) in UKBB. Using the Inter99 cohort, we found no difference in J-point amplitude and QTc-interval between carriers and non-carriers for 14/18 variants. There was no difference in the risk of syncope (P = 0.32), malignant ventricular arrhythmia (P = 0.96), and all-cause mortality (P = 0.59) between carriers and non-carriers. The ACMG guidelines reclassified 75% of all variants as variant-of-uncertain significance, likely benign, and benign. We identified ~2/3 of variants previously associated with SIDS and found no significant associations with electrocardiographic traits, syncope, malignant ventricular arrhythmia, or all-cause mortality. These data indicate that many of these variants are not highly penetrant, monogenic causes of SIDS and underline the importance of frequent reappraisal of genetic variants to avoid future misdiagnosis.

Correction: Loss of the Mia40a oxidoreductase leads to hepato-pancreatic insufficiency in zebrafish.

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

Loss of the Mia40a oxidoreductase leads to hepato-pancreatic insufficiency in zebrafish.

Development and function of tissues and organs are powered by the activity of mitochondria. In humans, inherited genetic mutations that lead to progressive mitochondrial pathology often manifest during infancy and can lead to death, reflecting the indispensable nature of mitochondrial biogenesis and function. Here, we describe a zebrafish mutant for the gene mia40a (chchd4a), the life-essential homologue of the evolutionarily conserved Mia40 oxidoreductase which drives the biogenesis of cysteine-rich mitochondrial proteins. We report that mia40a mutant animals undergo progressive cellular respiration defects and develop enlarged mitochondria in skeletal muscles before their ultimate death at the larval stage. We generated a deep transcriptomic and proteomic resource that allowed us to identify abnormalities in the development and physiology of endodermal organs, in particular the liver and pancreas. We identify the acinar cells of the exocrine pancreas to be severely affected by mutations in the MIA pathway. Our data contribute to a better understanding of the molecular, cellular and organismal effects of mitochondrial deficiency, important for the accurate diagnosis and future treatment strategies of mitochondrial diseases.

Rare truncating variants in the sarcomeric protein titin associate with familial and early-onset atrial fibrillation.

A family history of atrial fibrillation constitutes a substantial risk of developing the disease, however, the pathogenesis of this complex disease is poorly understood. We perform whole-exome sequencing on 24 families with at least three family members diagnosed with atrial fibrillation (AF) and find that titin-truncating variants (TTNtv) are significantly enriched in these patients (P = 1.76 × 10-6). This finding is replicated in an independent cohort of early-onset lone AF patients (n = 399; odds ratio = 36.8; P = 4.13 × 10-6). A CRISPR/Cas9 modified zebrafish carrying a truncating variant of titin is used to investigate TTNtv effect in atrial development. We observe compromised assembly of the sarcomere in both atria and ventricle, longer PR interval, and heterozygous adult zebrafish have a higher degree of fibrosis in the atria, indicating that TTNtv are important risk factors for AF. This aligns with the early onset of the disease and adds an important dimension to the understanding of the molecular predisposition for AF.

PKD Phosphorylation as Novel Pathway of KV11.1 Regulation.

BACKGROUND/AIMS: The voltage-gated potassium channel KV11.1 has been originally cloned from the brain and is expressed in a variety of tissues. The role of phosphorylation for channel function is a matter of debate. In this study, we aimed to elucidate the extent and role of protein kinase D mediated phosphorylation. METHODS: We employed mass spectrometry, whole-cell patch clamp electrophysiology, confocal microscopy, site-directed mutagenesis, and western blotting. RESULTS: Using brain tissue from rat and mouse, we mapped several phosphorylated KV11.1 residues by LC-MS mass spectrometry and identified protein kinase D (PKD1) as possible regulatory kinase. Co-expression of KV11.1 with PKD1 reduced current amplitudes without altering protein levels or surface expression of the channel. Based on LC-MS results from in vivo and HEK293 cell experiments we chose four KV11.1 mutant candidates for further functional analysis. Ablation of the putative PKD phosphorylation site in the mutant S284A increased the maximal current indicating S284 as a main PKD target in KV11.1. CONCLUSIONS: Our data might help mitigating a long-standing controversy in the field regarding PKC regulation of KV11.1. We propose that PKD1 mediates the PKC effects on KV11.1 and we found that PKD targets S284 in the N-terminus of the channel.

Microtubule Regulation of Kv7 Channels Orchestrates cAMP-Mediated Vasorelaxations in Rat Arterial Smooth Muscle.

Microtubules can regulate GPCR (G protein-coupled receptor) signaling in various cell types. In vascular smooth muscle, activation of the ß-adrenoceptor leads to production of cAMP to mediate a vasorelaxation. Little is known about the role of microtubules in smooth muscle, and given the importance of this pathway in vascular smooth muscle cells, we investigated the role of microtubule stability on ß-adrenoceptor signaling in rat renal and mesenteric arteries. In isometric tension experiments, incubation with the microtubule inhibitors colchicine and nocodazole enhanced isoprenaline-mediated relaxations of renal and mesenteric arteries that the microtubule stabilizer, paclitaxel, prevented. Sharp microelectrode experiments showed that colchicine treatment caused increased hyperpolarization of mesenteric artery segments in response to isoprenaline. Application of the Kv7 channel blocker, XE991, attenuated the effect of colchicine on isoprenaline relaxations, whereas iberiotoxin-a BKCa channel blocker-had no effect. In addition, colchicine improved the relaxations to the Kv7.2 to 7.5 activator, S-1, in both renal and mesenteric artery segments compared with dimethyl sulfoxide incubation. We determined that increased mesenteric artery myocytes treated with colchicine showed increased Kv7.4 membrane expression, but Western blot analysis showed no change in total Kv7.4 protein. This study is the first to show microtubule disruption improves the ß-adrenoceptor-mediated relaxations of mesenteric and renal arteries and determine this enhancement to be because of increased membrane expression of the Kv7 voltage-gated potassium channels.

Termination of Vernakalant-Resistant Atrial Fibrillation by Inhibition of Small-Conductance Ca2+-Activated K+ Channels in Pigs.

BACKGROUND: Evidence has emerged that small-conductance Ca2+-activated K+ (SK) channels constitute a new target for treatment of atrial fibrillation (AF). SK channels are predominantly expressed in the atria as compared with the ventricles. Various marketed antiarrhythmic drugs are limited by ventricular adverse effects and efficacy loss as AF progresses. METHODS AND RESULTS: A total of 43 pigs were used for the studies. AF reversion in conscious long-term tachypaced pigs: Pigs were subjected to atrial tachypacing (7 Hz) until they developed sustained AF that could not be reverted by vernakalant 4 mg/kg (18.8±3.3 days of atrial tachypacing). When the SK channel inhibitor AP14145 was tested in these animals, vernakalant-resistant AF was reverted to sinus rhythm, and reinduction of AF by burst pacing (50 Hz) was prevented in 8 of 8 pigs. Effects on refractory period and AF duration in open chest pigs: The effects of AP14145 and vernakalant on the effective refractory periods and acute burst pacing-induced AF were examined in anaesthetized open chest pigs. Both vernakalant and AP14145 significantly prolonged atrial refractoriness and reduced AF duration without affecting the ventricular refractoriness or blood pressure in pigs subjected to 7 days atrial tachypacing, as well as in sham-operated control pigs. CONCLUSIONS: SK currents play a role in porcine atrial repolarization, and pharmacological inhibition of these with AP14145 demonstrates antiarrhythmic effects in a vernakalant-resistant porcine model of AF. These results suggest SK channel blockers as potentially interesting anti-AF drugs.