Identification of heart rate-associated loci and their effects on cardiac conduction and rhythm disorders

HCN4 Gene Polymorphisms Are Associated With Occurrence of Tachycardia-Induced Cardiomyopathy in Patients With Atrial Fibrillation

BACKGROUND

Tachycardia-induced cardiomyopathy (TIC) is a reversible cardiomyopathy induced by tachyarrhythmia, and the genetic background of the TIC is not well understood. The hyperpolarization-activated cyclic nucleotide-gated channel gene HCN4 is highly expressed in the conduction system where it is involved in heart rate control. We speculated that the HCN4 gene is associated with TIC.

METHODS

We enrolled 930 Japanese patients with atrial fibrillation (AF) for screening, 350 Japanese patients with AF for replication, and 1635 non-AF controls. In the screening AF set, we compared HCN4 single-nucleotide polymorphism genotypes between AF subjects with TIC (TIC, n=73) and without TIC (non-TIC, n=857). Of 17 HCN4 gene-tag single-nucleotide polymorphisms, rs7172796, rs2680344, rs7164883, rs11631816, and rs12905211 were significantly associated with TIC. Among them, only rs7164883 was independently associated with TIC after conditional analysis (TIC versus non-TIC: minor allele frequency, 26.0% versus 9.7%; P=1.62×10^-9; odds ratio=3.2).

RESULTS

We confirmed this association of HCN4 single-nucleotide polymorphism rs7164883 with TIC in the replication set (TIC=41 and non-TIC=309; minor allele frequency, 28% versus 9.9%; P=1.94×10^-6; odds ratio=3.6). The minor allele frequency of rs7164883 was similar in patients with AF and non-AF controls (11% versus 10.9%; P=0.908).

CONCLUSIONS

The HCN4 gene single-nucleotide polymorphism rs7164883 may be a new genetic marker for TIC in patients with AF.

Identifying blood-specific age-related DNA methylation markers on the Illumina MethylationEPIC® BeadChip

The past decade has seen rapid development in DNA methylation (DNAm) microarrays, including the Illumina HumanMethylation27 and HumanMethylation450 (450K) chips, which have played an essential role in identifying and evaluating age-related (AR) DNAm markers in different tissues. Recently, a new array, the Illumina MethylationEPIC (EPIC) was introduced, with nearly double the number of probes as the 450K (∼850,000 probes). In this study, we test these newly added probes for age association using a large cohort of 754 DNAm profiles from blood samples assayed on the EPIC BeadChip, for individuals aged 0-88 years old. 52 AR CpG sites (Spearman's abs(rho) >0.6 and P-value <10^-83) were identified, 21 of which were novel sites and mapped to 18 genes, nine of which (LHFPL4, SLC12A8, EGFEM1P, GPR158, TAL1, KIAA1755, LOC730668, DUSP16, and FAM65C) have never previously been reported to be associated with age. The data were subsequently split into a 527-sample training set and a 227-sample testing set to build and validate two age prediction models using elastic net regression and multivariate regression. Elastic net regression selected 425 CpG markers with a mean absolute deviation (MAD) of 2.6 years based on the testing set. To build a multivariate linear regression model, AR CpG sites with R² > 0.5 at FDR < 0.05 were input into stepwise regression to select the best subset for age prediction. The resulting six CpG markers were linearly modelled with age and explained 81% of age-correlated variation in DNAm levels. Age estimation accuracy using bootstrap analysis was 4.5 years, with 95% confidence intervals of 4.56 to 4.57 years based on the testing set. These results suggest that EPIC BeadChip probes for age estimation fall within the range of probes found on the previous Illumina HumanMethylation platforms in terms of their age-prediction ability.

Transcriptional regulation of the cardiac conduction system

The rate and rhythm of heart muscle contractions are coordinated by the cardiac conduction system (CCS), a generic term for a collection of different specialized muscular tissues within the heart. The CCS components initiate the electrical impulse at the sinoatrial node, propagate it from atria to ventricles via the atrioventricular node and bundle branches, and distribute it to the ventricular muscle mass via the Purkinje fibre network. The CCS thereby controls the rate and rhythm of alternating contractions of the atria and ventricles. CCS function is well conserved across vertebrates from fish to mammals, although particular specialized aspects of CCS function are found only in endotherms (mammals and birds). The development and homeostasis of the CCS involves transcriptional and regulatory networks that act in an embryonic-stage-dependent, tissue-dependent, and dose-dependent manner. This Review describes emerging data from animal studies, stem cell models, and genome-wide association studies that have provided novel insights into the transcriptional networks underlying CCS formation and function. How these insights can be applied to develop disease models and therapies is also discussed.

Association between resting heart rate and incident diabetes risk: a Mendelian randomization study

AIMS

Observational studies indicated that resting heart rate (RHR) was associated with diabetes mellitus (DM) risk; however, it remains unclear whether the association between RHR and DM is causal. We aimed to examine whether there was causal association of RHR with DM risk.

METHODS

A prospective study including 16,201 middle-aged and older Chinese (7031 males and 9170 females) derived from the Dongfeng-Tongji cohort was performed. Cox proportional hazard regression models were conducted to estimate the associations between RHR and incident DM risk. In 7481 participants, 65 single nucleotide polymorphisms related to RHR were genotyped. A genetic risk score (GRS) of RHR was calculated based on the RHR-associated variants. The causal associations of RHR with DM risk were investigated by Mendelian randomization analysis.

RESULTS

During a mean (SD) follow-up of 4.5 (0.5) years, 1110 diabetes were identified. Compared with the referential RHR group (≤ 60 beats per minute [bpm]), individuals with RHR > 80 bpm have a higher incident diabetes risk, with a hazard ratio of 1.40 (95% confidence interval [CI], 1.05-1.88). With per SD increase in the weighted genetic risk score, the resting heart rate increased by 0.71 bpm (95% CI 0.49-0.93). By using the GRS to estimate the unconfounded effect, we found that higher resting heart rate did not have a causal effect on diabetes risk (OR 1.00 [95% CI 0.95-1.05]).

CONCLUSIONS

The present study supported a positive but not a causal association of RHR with incident diabetes risk. More studies are needed to verify our findings.

Using zebrafish to study skeletal genomics

While genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of skeletal diseases, animal models are required to identify causal mechanisms and to translate underlying biology into new therapies. Despite large-scale knockout mouse phenotyping efforts, the skeletal functions of most genes residing at GWAS-identified loci remain unknown, highlighting a need for complementary model systems to accelerate gene discovery. Over the past several decades, zebrafish (Danio rerio) has emerged as a powerful system for modeling the genetics of human diseases. In this review, our goal is to outline evidence supporting the utility of zebrafish for accelerating our understanding of human skeletal genomics, as well as gaps in knowledge that need to be filled for this purpose. We do this by providing a basic foundation of the zebrafish skeletal morphophysiology and phenotypes, and surveying evidence of skeletal gene homology and the use of zebrafish for post-GWAS analysis in other tissues and organs. We also outline challenges in translating zebrafish mutant phenotypes. Finally, we conclude with recommendations of future directions and how to leverage the large body of tools and knowledge of skeletal genetics in zebrafish for the needs of human skeletal genomic exploration. Due to their amenability to rapid genetic approaches, as well as the large number of conserved genetic and phenotypic features, there is a strong rationale supporting the use of zebrafish for human skeletal genomic studies.

Impact of Oxidative Stress on Long-Term Heart Rate Variability: Linear Versus Non-Linear Heart Rate Dynamics

BACKGROUND

Heart rate variability (HRV) is a widely used non-invasive and quantitative marker of cardiac autonomic control. Elevated oxidative stress (OS) and reduced HRV have been proven in specific disease subsets. However, the impact of OS on the long-term heart rate dynamics of both conventional linear and non-linear origin in the general population is not known.

METHODS

The 24-hour ambulatory electrocardiogram recordings and plasma 8-iso-prostaglandin F2α (8-iso-PGF2α) levels as an OS marker were acquired simultaneously in 71 consecutive patients. The conventional time and frequency domain HRV parameters and non-linear parameters were measured.

RESULTS

The 8-iso-PGF2α is a significant determinant of most long-term conventional time and frequency domain HRV parameters and standard deviation (SD1, perpendicular to the line of identity; SD2, along the line of identity) descriptors from Poincaré plot analysis, but not of non-linear complexity and fractal parameters. Patients with a high OS burden had lower absolute low-frequency and high-frequency powers during both the night and morning periods, with a significant decrease in high-frequency power in the morning.

CONCLUSIONS

Oxidative stress is one of the significant determinants of the HRV. The severity of OS is reflected in the conventional time and frequency domain HRV parameters, but not in the non-linear measurements.

Transcriptomic Profiling of the Developing Cardiac Conduction System at Single-Cell Resolution

RATIONALE

The cardiac conduction system (CCS) consists of distinct components including the sinoatrial node, atrioventricular node, His bundle, bundle branches, and Purkinje fibers. Despite an essential role for the CCS in heart development and function, the CCS has remained challenging to interrogate because of inherent obstacles including small cell numbers, large cell-type heterogeneity, complex anatomy, and difficulty in isolation. Single-cell RNA-sequencing allows for genome-wide analysis of gene expression at single-cell resolution.

OBJECTIVE

Assess the transcriptional landscape of the entire CCS at single-cell resolution by single-cell RNA-sequencing within the developing mouse heart.

METHODS AND RESULTS

Wild-type, embryonic day 16.5 mouse hearts (n=6 per zone) were harvested and 3 zones of microdissection were isolated, including: Zone I-sinoatrial node region; Zone II-atrioventricular node/His region; and Zone III-bundle branch/Purkinje fiber region. Tissue was digested into single-cell suspensions, cells isolated, mRNA reverse transcribed, and barcoded before high-throughput sequencing and bioinformatics analyses. Single-cell RNA-sequencing was performed on over 22 000 cells, and all major cell types of the murine heart were successfully captured including bona fide clusters of cells consistent with each major component of the CCS. Unsupervised weighted gene coexpression network analysis led to the discovery of a host of novel CCS genes, a subset of which were validated using fluorescent in situ hybridization as well as whole-mount immunolabeling with volume imaging (iDISCO+) in 3 dimensions on intact mouse hearts. Further, subcluster analysis unveiled isolation of distinct CCS cell subtypes, including the clinically relevant but poorly characterized transitional cells that bridge the CCS and surrounding myocardium.

CONCLUSIONS

Our study represents the first comprehensive assessment of the transcriptional profiles from the entire CCS at single-cell resolution and provides a characterization in the context of development and disease.

GWAS of QRS duration identifies new loci specific to Hispanic/Latino populations

BACKGROUND

The electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.

METHODS

We performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.

RESULTS

We identified six loci associated with QRS (P<5x10-8), including two novel loci: MYOCD, a nuclear protein expressed in the heart, and SYT1, an integral membrane protein. The top SNP in the MYOCD locus, intronic SNP rs16946539, was found in Hispanics/Latinos with a minor allele frequency (MAF) of 0.04, but is monomorphic in European and African descent populations. The most significant QRS duration association was with intronic SNP rs3922344 (P = 1.19x10-24) in SCN5A/SCN10A. Three other previously identified loci, CDKN1A, VTI1A, and HAND1, also exceeded the GWAS significance threshold among Hispanics/Latinos. A total of 27 of 32 previously identified QRS duration SNPs were shown to generalize in Hispanics/Latinos.

CONCLUSIONS

Our QRS duration GWAS, the first in Hispanic/Latino populations, identified two new loci, underscoring the utility of extending large scale genomic studies to currently under-examined populations.

Effects of Calcium, Magnesium, and Potassium Concentrations on Ventricular Repolarization in Unselected Individuals

BACKGROUND

Subclinical changes on the electrocardiogram are risk factors for cardiovascular mortality. Recognition and knowledge of electrolyte associations in cardiac electrophysiology are based on only in vitro models and observations in patients with severe medical conditions.

OBJECTIVES

This study sought to investigate associations between serum electrolyte concentrations and changes in cardiac electrophysiology in the general population.

METHODS

Summary results collected from 153,014 individuals (54.4% women; mean age 55.1 ± 12.1 years) from 33 studies (of 5 ancestries) were meta-analyzed. Linear regression analyses examining associations between electrolyte concentrations (mmol/l of calcium, potassium, sodium, and magnesium), and electrocardiographic intervals (RR, QT, QRS, JT, and PR intervals) were performed. The study adjusted for potential confounders and also stratified by ancestry, sex, and use of antihypertensive drugs.

RESULTS

Lower calcium was associated with longer QT intervals (-11.5 ms; 99.75% confidence interval [CI]: -13.7 to -9.3) and JT duration, with sex-specific effects. In contrast, higher magnesium was associated with longer QT intervals (7.2 ms; 99.75% CI: 1.3 to 13.1) and JT. Lower potassium was associated with longer QT intervals (-2.8 ms; 99.75% CI: -3.5 to -2.0), JT, QRS, and PR durations, but all potassium associations were driven by use of antihypertensive drugs. No physiologically relevant associations were observed for sodium or RR intervals.

CONCLUSIONS

The study identified physiologically relevant associations between electrolytes and electrocardiographic intervals in a large-scale analysis combining cohorts from different settings. The results provide insights for further cardiac electrophysiology research and could potentially influence clinical practice, especially the association between calcium and QT duration, by which calcium levels at the bottom 2% of the population distribution led to clinically relevant QT prolongation by >5 ms.

Genome-wide association meta-analysis of 30,000 samples identifies seven novel loci for quantitative ECG traits

Genome-wide association studies (GWAS) of quantitative electrocardiographic (ECG) traits in large consortia have identified more than 130 loci associated with QT interval, QRS duration, PR interval, and heart rate (RR interval). In the current study, we meta-analyzed genome-wide association results from 30,000 mostly Dutch samples on four ECG traits: PR interval, QRS duration, QT interval, and RR interval. SNP genotype data was imputed using the Genome of the Netherlands reference panel encompassing 19 million SNPs, including millions of rare SNPs (minor allele frequency < 5%). In addition to many known loci, we identified seven novel locus-trait associations: KCND3, NR3C1, and PLN for PR interval, KCNE1, SGIP1, and NFKB1 for QT interval, and ATP2A2 for QRS duration, of which six were successfully replicated. At these seven loci, we performed conditional analyses and annotated significant SNPs (in exons and regulatory regions), demonstrating involvement of cardiac-related pathways and regulation of nearby genes.

Biomarkers of Dietary Omega-6 Fatty Acids and Incident Cardiovascular Disease and Mortality

BACKGROUND

Global dietary recommendations for and cardiovascular effects of linoleic acid, the major dietary omega-6 fatty acid, and its major metabolite, arachidonic acid, remain controversial. To address this uncertainty and inform international recommendations, we evaluated how in vivo circulating and tissue levels of linoleic acid (LA) and arachidonic acid (AA) relate to incident cardiovascular disease (CVD) across multiple international studies.

METHODS

We performed harmonized, de novo, individual-level analyses in a global consortium of 30 prospective observational studies from 13 countries. Multivariable-adjusted associations of circulating and adipose tissue LA and AA biomarkers with incident total CVD and subtypes (coronary heart disease, ischemic stroke, cardiovascular mortality) were investigated according to a prespecified analytic plan. Levels of LA and AA, measured as the percentage of total fatty acids, were evaluated linearly according to their interquintile range (ie, the range between the midpoint of the first and fifth quintiles), and categorically by quintiles. Study-specific results were pooled using inverse-variance-weighted meta-analysis. Heterogeneity was explored by age, sex, race, diabetes mellitus, statin use, aspirin use, omega-3 levels, and fatty acid desaturase 1 genotype (when available).

RESULTS

In 30 prospective studies with medians of follow-up ranging 2.5 to 31.9 years, 15 198 incident cardiovascular events occurred among 68 659 participants. Higher levels of LA were significantly associated with lower risks of total CVD, cardiovascular mortality, and ischemic stroke, with hazard ratios per interquintile range of 0.93 (95% CI, 0.88-0.99), 0.78 (0.70-0.85), and 0.88 (0.79-0.98), respectively, and nonsignificantly with lower coronary heart disease risk (0.94; 0.88-1.00). Relationships were similar for LA evaluated across quintiles. AA levels were not associated with higher risk of cardiovascular outcomes; in a comparison of extreme quintiles, higher levels were associated with lower risk of total CVD (0.92; 0.86-0.99). No consistent heterogeneity by population subgroups was identified in the observed relationships.

CONCLUSIONS

In pooled global analyses, higher in vivo circulating and tissue levels of LA and possibly AA were associated with lower risk of major cardiovascular events. These results support a favorable role for LA in CVD prevention.

β2 -Adrenergic Receptor Gene Affects the Heart Rate Response of β-Blockers: Evidence From 3 Clinical Studies

β-Blockers' heart rate (HR)-lowering effect is an important determinant of the effectiveness for this class of drugs, yet it is variable among β-blocker-treated patients. To date, genetic studies have revealed several genetic signals associated with HR response to β-blockers. However, these genetic signals have not been consistently replicated across multiple independent cohorts. Here we sought to use data from 3 hypertension clinical trials to validate single-nucleotide polymorphisms (SNPs) previously associated with the HR response to β-blockers. Using linear regression analysis, we investigated the effects of 6 SNPs in 3 genes, including ADRB1, ADRB2, and GNB3, relative to the HR response following β-blocker used in the PEAR (n = 757), PEAR-2 (n = 368), and INVEST (n = 1401) trials, adjusting for baseline HR, age, sex, and ancestry. Atenolol was used in PEAR and INVEST, and metoprolol was used in PEAR-2. We found that rs1042714 and rs1042713 in ADRB2 were significantly associated with HR response to both β-blockers in whites (rs1042714 C-allele carriers, meta-analysis β = -0.95 beats per minute [bpm], meta-analysis P = 3×10^-4 ; rs1042713 A-allele carriers, meta-analysis β = -1.15 bpm, meta-analysis P = 2×10^-3 ). In conclusion, the results of our analyses provide strong evidence to support the hypothesis that rs1042714 and rs1042713 in the ADRB2 gene are important predictors of HR response to cardioselective β-blockade in hypertensive patient cohorts.

Transcriptome analysis of mouse and human sinoatrial node cells reveals a conserved genetic program

The rate of contraction of the heart relies on proper development and function of the sinoatrial node, which consists of a small heterogeneous cell population, including Tbx3⁺ pacemaker cells. Here, we have isolated and characterized the Tbx3⁺ cells from Tbx3 ^+/Venus knock-in mice. We studied electrophysiological parameters during development and found that Venus-labeled cells are genuine Tbx3⁺ pacemaker cells. We analyzed the transcriptomes of late fetal FACS-purified Tbx3⁺ sinoatrial nodal cells and Nppb-Katushka⁺ atrial and ventricular chamber cardiomyocytes, and identified a sinoatrial node-enriched gene program, including key nodal transcription factors, BMP signaling and Smoc2, the disruption of which in mice did not affect heart rhythm. We also obtained the transcriptomes of the sinoatrial node region, including pacemaker and other cell types, and right atrium of human fetuses, and found a gene program including TBX3, SHOX2, ISL1 and HOX family members, and BMP and NOTCH signaling components conserved between human and mouse. We conclude that a conserved gene program characterizes the sinoatrial node region and that the Tbx3 ^+/Venus allele provides a reliable tool for visualizing the sinoatrial node, and studying its development and function.

Phenotypes associated with genes encoding drug targets are predictive of clinical trial side effects

Only a small fraction of early drug programs progress to the market, due to safety and efficacy failures, despite extensive efforts to predict safety. Characterizing the effect of natural variation in the genes encoding drug targets should present a powerful approach to predict side effects arising from drugging particular proteins. In this retrospective analysis, we report a correlation between the organ systems affected by genetic variation in drug targets and the organ systems in which side effects are observed. Across 1819 drugs and 21 phenotype categories analyzed, drug side effects are more likely to occur in organ systems where there is genetic evidence of a link between the drug target and a phenotype involving that organ system, compared to when there is no such genetic evidence (30.0 vs 19.2%; OR = 1.80). This result suggests that human genetic data should be used to predict safety issues associated with drug targets.

Safety issues including side effects are one of the major factors causing failure of clinical trials in drug development. Here, the authors leverage information about phenotypes associated with variation in genes encoding drug targets to predict drug-treatment-related side effects.

Genome-Wide Association Studies of Hypertension and Several Other Cardiovascular Diseases

Genome-wide association studies (GWAS) have greatly expanded our understanding of the genetic architecture of cardiovascular diseases in the past decade. They have revealed hundreds of suggestive genetic loci that replicate known biological candidate genes and indicate the existence of a previously unsuspected new biology relevant to cardiovascular disorders. These data have been used successfully to create genetic risk scores that may improve risk prediction and the identification of susceptive individuals. Furthermore, these GWAS-identified novel pathways may herald a new era of novel drug development and stratification of patients. In this review, we will briefly summarize the literature on the candidate genes and signals discovered by GWAS on hypertension and coronary artery disease and discuss their implications on clinical medicine.

DNA damage-induced PARP1 activation confers cardiomyocyte dysfunction through NAD+ depletion in experimental atrial fibrillation

Atrial fibrillation (AF) is the most common clinical tachyarrhythmia with a strong tendency to progress in time. AF progression is driven by derailment of protein homeostasis, which ultimately causes contractile dysfunction of the atria. Here we report that tachypacing-induced functional loss of atrial cardiomyocytes is precipitated by excessive poly(ADP)-ribose polymerase 1 (PARP1) activation in response to oxidative DNA damage. PARP1-mediated synthesis of ADP-ribose chains in turn depletes nicotinamide adenine dinucleotide (NAD+), induces further DNA damage and contractile dysfunction. Accordingly, NAD+ replenishment or PARP1 depletion precludes functional loss. Moreover, inhibition of PARP1 protects against tachypacing-induced NAD+ depletion, oxidative stress, DNA damage and contractile dysfunction in atrial cardiomyocytes and Drosophila. Consistently, cardiomyocytes of persistent AF patients show significant DNA damage, which correlates with PARP1 activity. The findings uncover a mechanism by which tachypacing impairs cardiomyocyte function and implicates PARP1 as a possible therapeutic target that may preserve cardiomyocyte function in clinical AF.

Genomics of Cardiovascular Measures of Autonomic Tone

The autonomic nervous system exerts broad control over the involuntary functions of the human body through complex equilibrium between sympathetic and parasympathetic tone. Imbalance in this equilibrium is associated with a multitude of cardiovascular outcomes, including mortality. The cardiovascular static state of this equilibrium can be quantified using physiological parameters such as heart rate (HR), blood pressure, and by spectral analysis of HR variability. Here, we review the current state of knowledge of the genetic background of cardiovascular measurements of autonomic tone. For most parameters of autonomic tone, a large portion of variability is explained by genetic heritability. Many of the static parameters of autonomic tone have also been studied through candidate-gene approach, yielding some insight into how genotypes of adrenergic receptors affect variables such as HR. Genome-wide approaches in large cohorts similarly exist for static variables such as HR and blood pressure but less is known about the genetic background of the dynamic and more specific measurements, such as HR variability. Furthermore, because most autonomic measures are likely polygenic, pathway analyses and modeling of polygenic effects are critical. Future work will hopefully explain the control of autonomic tone and guide individualized therapeutic interventions.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

Resting heart rate and relation to disease and longevity: past, present and future

Assessment of heart rate has been used for millennia as a marker of health. Several studies have indicated that low resting heart rate (RHR) is associated with health and longevity, and conversely, a high resting heart to be associated with disease and adverse events. Longitudinal studies have shown a clear association between increase in heart rate over time and adverse events. RHR is a fundamental clinical characteristic and several trials have assessed the effectiveness of heart rate lowering medication, for instance beta-blockers and selective sinus node inhibition. Advances in technology have provided new insights into genetic factors related to RHR as well as insights into whether elevated RHR is a risk factor or risk marker. Recent animal research has suggested that heart rate lowering with sinus node inhibition is associated with increased lifespan. Furthermore, genome-wide association studies in the general population using Mendelian randomization have demonstrated a causal link between heart rate at rest and longevity. Furthermore, the development in personal digital devices such as mobile phones, fitness trackers and eHealth applications has made heart rate information and knowledge in this field as important as ever for the public as well as the clinicians. It should therefore be expected that clinicians and health care providers will be met by relevant questions and need of advice regarding heart rate information from patients and the public. The present review provides an overview of the current knowledge in the field of heart rate and health.

Detection of Shared Balancing Selection in the Absence of Trans-Species Polymorphism

Trans-species polymorphism has been widely used as a key sign of long-term balancing selection across multiple species. However, such sites are often rare in the genome and could result from mutational processes or technical artifacts. Few methods are yet available to specifically detect footprints of trans-species balancing selection without using trans-species polymorphic sites. In this study, we develop summary- and model-based approaches that are each specifically tailored to uncover regions of long-term balancing selection shared by a set of species by using genomic patterns of intraspecific polymorphism and interspecific fixed differences. We demonstrate that our trans-species statistics have substantially higher power than single-species approaches to detect footprints of trans-species balancing selection, and are robust to those that do not affect all tested species. We further apply our model-based methods to human and chimpanzee whole-genome sequencing data. In addition to the previously established major histocompatibility complex and malaria resistance-associated FREM3/GYPE regions, we also find outstanding genomic regions involved in barrier integrity and innate immunity, such as the GRIK1/CLDN17 intergenic region, and the SLC35F1 and ABCA13 genes. Our findings not only echo the significance of pathogen defense but also reveal novel candidates in maintaining balanced polymorphisms across human and chimpanzee lineages. Finally, we show that these trans-species statistics can be applied to and work well for an arbitrary number of species, and integrate them into open-source software packages for ease of use by the scientific community.

Genetics of Atrial Fibrilation: In Search of Novel Therapeutic Targets

Atrial fibrillation (AF) is the most frequent arrhythmogenic disease in humans, ranging from 2% in the general population and rising up to 10-12% in 80+ years. Genetic analyses of AF familiar cases have identified a series of point mutations in distinct ion channels, supporting a causative link. However, these genetic defects only explain a minority of AF patients. Genomewide association studies identified single nucleotide polymorphisms (SNPs), close to PITX2 on 4q25 chromosome, that are highly associated to AF. Subsequent GWAS studies have identified several new loci, involving additional transcription and growth factors. Furthermore, these risk 4q25 SNPs serve as surrogate biomarkers to identify AF recurrence in distinct surgical and pharmacological interventions. Experimental studies have demonstrated an intricate signalling pathway supporting a key role of the homeobox transcription factor PITX2 as a transcriptional regulator. Furthermore, cardiovascular risk factors such as hyperthyroidism, hypertension and redox homeostasis have been identified to modulate PITX2 driven gene regulatory networks. We provide herein a state-of-the-art review of the genetic bases of atrial fibrillation, our current understanding of the genetic regulatory networks involved in AF and its plausible usage for searching novel therapeutic targets.

Examining the current standards for genetic discovery and replication in the era of mega-biobanks

With the recent deluge of mega-biobank data, it is time to revisit what constitutes “replication” for genome-wide association studies. Many replication samples are unavailable or underpowered, therefore alternatives beyond strict statistical replication are needed until the required resources become available.

Genome-wide association studies (GWAS) have become a mainstay in genetics research to understand genotype-phenotype relationships. Following the second release of UK Biobank data and the flood of publications using these data, here the author revisits the standards for discovery, replication and follow-up in GWAS today.

Exploring Coronary Artery Disease GWAs Targets With Functional Links to Immunometabolism

Finding genetic variants that cause functional disruption or regulatory change among the many implicated GWAs variants remains a key challenge to translating the findings from GWAs to therapeutic treatments. Defining the causal mechanisms behind the variants require functional screening experiments that can be complex and costly. Prioritizing variants for functional characterization using techniques that capture important functional and regulatory elements can assist this. The genetic architecture of complex traits such as cardiovascular disease and type II diabetes comprise an enormously large number of variants of small effect contributing to heritability and spread throughout the genome. This makes it difficult to distinguish which variants or core genes are most relevant for prioritization and how they contribute to the regulatory networks that become dysregulated leading to disease. Despite these challenges, recent GWAs for CAD prioritized genes associated with lipid metabolism, coagulation and adhesion along with novel signals related to innate immunity, adipose tissue and, vascular function as important core drivers of risk. We focus on three examples of novel signals associated with CAD which affect risk through missense or UTR mutations indicating their potential for therapeutic modification. These variants play roles in adipose tissue function vascular function and innate immunity which form the cornerstones of immuno-metabolism. In addition we have explored the putative, but potentially important interactions between the environment, specifically food and nutrition, with respect to key processes.

Filamin-A as a Balance between Erk/Smad Activities During Cardiac Valve Development

Mitral valve prolapse (MVP) affects 2.4% of the population and has poorly understood etiology. Recent genetic studies have begun to unravel the complexities of MVP and through these efforts, mutations in the FLNA (Filamin-A) gene were identified as disease causing. Our in vivo and in vitro studies have validated these genetic findings and have revealed FLNA as a central regulator of valve morphogenesis. The mechanisms by which FLNA mutations result in myxomatous mitral valve disease are currently unknown, but may involve proteins previously associated with mutated regions of the FLNA protein, such as the small GTPase signaling protein, R-Ras. Herein, we report that Filamin-A is required for R-Ras expression and activation of the Ras-Mek-Erk pathway. Loss of the Ras/Erk pathway correlated with hyperactivation of pSmad2/3, increased extracellular matrix (ECM) production and enlarged mitral valves. Analyses of integrin receptors in the mitral valve revealed that Filamin-A was required for β1-integrin expression and provided a potential mechanism for impaired ECM compaction and valve enlargement. Our data support Filamin-A as a protein that regulates the balance between Erk and Smad activation and an inability of Filamin-A deficient valve interstitial cells to effectively remodel the increased ECM production through a β1-integrin mechanism. As a consequence, loss of Filamin-A function results in increased ECM production and generation of a myxomatous phenotype characterized by improperly compacted mitral valve tissue. Anat Rec, 302:117-124, 2019. © 2018 Wiley Periodicals, Inc.

Detecting past and ongoing natural selection among ethnically Tibetan women at high altitude in Nepal

Adaptive evolution in humans has rarely been characterized for its whole set of components, i.e. selective pressure, adaptive phenotype, beneficial alleles and realized fitness differential. We combined approaches for detecting polygenic adaptations and for mapping the genetic bases of physiological and fertility phenotypes in approximately 1000 indigenous ethnically Tibetan women from Nepal, adapted to high altitude. The results of genome-wide association analyses and tests for polygenic adaptations showed evidence of positive selection for alleles associated with more pregnancies and live births and evidence of negative selection for those associated with higher offspring mortality. Lower hemoglobin level did not show clear evidence for polygenic adaptation, despite its strong association with an EPAS1 haplotype carrying selective sweep signals.

Anderson on behalf of METASTROKE and the ISGC, Boehnke M, Laakso M, Atzmon G, Glaser B, Mercader JM, Gaulton K, Flannick J, Getz G, Florez JC. Type 2 diabetes genetic loci informed by multi-trait associations point to disease mechanisms and subtypes: A soft clustering analysis

BACKGROUND

Type 2 diabetes (T2D) is a heterogeneous disease for which (1) disease-causing pathways are incompletely understood and (2) subclassification may improve patient management. Unlike other biomarkers, germline genetic markers do not change with disease progression or treatment. In this paper, we test whether a germline genetic approach informed by physiology can be used to deconstruct T2D heterogeneity. First, we aimed to categorize genetic loci into groups representing likely disease mechanistic pathways. Second, we asked whether the novel clusters of genetic loci we identified have any broad clinical consequence, as assessed in four separate subsets of individuals with T2D.

METHODS AND FINDINGS

In an effort to identify mechanistic pathways driven by established T2D genetic loci, we applied Bayesian nonnegative matrix factorization (bNMF) clustering to genome-wide association study (GWAS) results for 94 independent T2D genetic variants and 47 diabetes-related traits. We identified five robust clusters of T2D loci and traits, each with distinct tissue-specific enhancer enrichment based on analysis of epigenomic data from 28 cell types. Two clusters contained variant-trait associations indicative of reduced beta cell function, differing from each other by high versus low proinsulin levels. The three other clusters displayed features of insulin resistance: obesity mediated (high body mass index [BMI] and waist circumference [WC]), "lipodystrophy-like" fat distribution (low BMI, adiponectin, and high-density lipoprotein [HDL] cholesterol, and high triglycerides), and disrupted liver lipid metabolism (low triglycerides). Increased cluster genetic risk scores were associated with distinct clinical outcomes, including increased blood pressure, coronary artery disease (CAD), and stroke. We evaluated the potential for clinical impact of these clusters in four studies containing individuals with T2D (Metabolic Syndrome in Men Study [METSIM], N = 487; Ashkenazi, N = 509; Partners Biobank, N = 2,065; UK Biobank [UKBB], N = 14,813). Individuals with T2D in the top genetic risk score decile for each cluster reproducibly exhibited the predicted cluster-associated phenotypes, with approximately 30% of all individuals assigned to just one cluster top decile. Limitations of this study include that the genetic variants used in the cluster analysis were restricted to those associated with T2D in populations of European ancestry.

CONCLUSION

Our approach identifies salient T2D genetically anchored and physiologically informed pathways, and supports the use of genetics to deconstruct T2D heterogeneity. Classification of patients by these genetic pathways may offer a step toward genetically informed T2D patient management.

PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

Exome-chip meta-analysis identifies novel loci associated with cardiac conduction, including ADAMTS6

BACKGROUND

Genome-wide association studies conducted on QRS duration, an electrocardiographic measurement associated with heart failure and sudden cardiac death, have led to novel biological insights into cardiac function. However, the variants identified fall predominantly in non-coding regions and their underlying mechanisms remain unclear.

RESULTS

Here, we identify putative functional coding variation associated with changes in the QRS interval duration by combining Illumina HumanExome BeadChip genotype data from 77,898 participants of European ancestry and 7695 of African descent in our discovery cohort, followed by replication in 111,874 individuals of European ancestry from the UK Biobank and deCODE cohorts. We identify ten novel loci, seven within coding regions, including ADAMTS6, significantly associated with QRS duration in gene-based analyses. ADAMTS6 encodes a secreted metalloprotease of currently unknown function. In vitro validation analysis shows that the QRS-associated variants lead to impaired ADAMTS6 secretion and loss-of function analysis in mice demonstrates a previously unappreciated role for ADAMTS6 in connexin 43 gap junction expression, which is essential for myocardial conduction.

CONCLUSIONS

Our approach identifies novel coding and non-coding variants underlying ventricular depolarization and provides a possible mechanism for the ADAMTS6-associated conduction changes.

miRNAS in cardiovascular diseases: potential biomarkers, therapeutic targets and challenges

Cardiovascular diseases (CVD) are the leading cause of morbidity and mortality in the world. Although considerable progress has been made in the diagnosis, treatment and prognosis of CVD, there is still a critical need for novel diagnostic biomarkers and new therapeutic interventions to decrease the incidence of this disease. Recently, there is increasing evidence that circulating miRNAs (miRNAs), i.e. endogenous, stable, single-stranded, short, non-coding RNAs, can be used as diagnostic biomarkers for CVD. Furthermore, miRNAs represent potential novel therapeutic targets for several cardiovascular disorders. In this review we provides an overview of the effects of several CVD; including heart failure, acute myocardial infarction, arrhythmias and pulmonary hypertension; on levels of circulating miRNAs. In addition, the use of miRNA as therapeutic targets is also discussed, as well as challenges and recommendations in their use in the diagnosis of CVD.

Promoter interactome of human embryonic stem cell-derived cardiomyocytes connects GWAS regions to cardiac gene networks

Long-range chromosomal interactions bring distal regulatory elements and promoters together to regulate gene expression in biological processes. By performing promoter capture Hi-C (PCHi-C) on human embryonic stem cell-derived cardiomyocytes (hESC-CMs), we show that such promoter interactions are a key mechanism by which enhancers contact their target genes after hESC-CM differentiation from hESCs. We also show that the promoter interactome of hESC-CMs is associated with expression quantitative trait loci (eQTLs) in cardiac left ventricular tissue; captures the dynamic process of genome reorganisation after hESC-CM differentiation; overlaps genome-wide association study (GWAS) regions associated with heart rate; and identifies new candidate genes in such regions. These findings indicate that regulatory elements in hESC-CMs identified by our approach control gene expression involved in ventricular conduction and rhythm of the heart. The study of promoter interactions in other hESC-derived cell types may be of utility in functional investigation of GWAS-associated regions.

Genomic approaches for the elucidation of genes and gene networks underlying cardiovascular traits

Genome-wide association studies have shed light on the association between natural genetic variation and cardiovascular traits. However, linking a cardiovascular trait associated locus to a candidate gene or set of candidate genes for prioritization for follow-up mechanistic studies is all but straightforward. Genomic technologies based on next-generation sequencing technology nowadays offer multiple opportunities to dissect gene regulatory networks underlying genetic cardiovascular trait associations, thereby aiding in the identification of candidate genes at unprecedented scale. RNA sequencing in particular becomes a powerful tool when combined with genotyping to identify loci that modulate transcript abundance, known as expression quantitative trait loci (eQTL), or loci modulating transcript splicing known as splicing quantitative trait loci (sQTL). Additionally, the allele-specific resolution of RNA-sequencing technology enables estimation of allelic imbalance, a state where the two alleles of a gene are expressed at a ratio differing from the expected 1:1 ratio. When multiple high-throughput approaches are combined with deep phenotyping in a single study, a comprehensive elucidation of the relationship between genotype and phenotype comes into view, an approach known as systems genetics. In this review, we cover key applications of systems genetics in the broad cardiovascular field.

Multi-ethnic genome-wide association study for atrial fibrillation

Atrial fibrillation (AF) affects more than 33 million individuals worldwide1 and has a complex heritability2. We conducted the largest meta-analysis of genome-wide association studies (GWAS) for AF to date, consisting of more than half a million individuals, including 65,446 with AF. In total, we identified 97 loci significantly associated with AF, including 67 that were novel in a combined-ancestry analysis, and 3 that were novel in a European-specific analysis. We sought to identify AF-associated genes at the GWAS loci by performing RNA-sequencing and expression quantitative trait locus analyses in 101 left atrial samples, the most relevant tissue for AF. We also performed transcriptome-wide analyses that identified 57 AF-associated genes, 42 of which overlap with GWAS loci. The identified loci implicate genes enriched within cardiac developmental, electrophysiological, contractile and structural pathways. These results extend our understanding of the biological pathways underlying AF and may facilitate the development of therapeutics for AF.

Functional Polymorphism in the MSI1 Gene Promoter Confers a Decreased Risk of Lung Cancer in Chinese by Reducing MSI1 Expression

Background:

Musashi1 (MSI1) is a characteristic stem cell marker that regulates the balance between cell self-renewal and differentiation. Evidence has identified MSI1 as a pivotal oncogenic regulator in diverse malignancies. However, little evidence uncovers the role of genetic variations of MSI1 gene in cancer etiology.

Objective:

The aim of this study was to investigate the association between genetic variants in the MSI1 gene and lung cancer risk.

Methods:

Based on a two-stage retrospective study with a total of 1559 patients with lung cancer and 1667 healthy controls, we evaluated the relevance between three putative functional SNPs in the MSI1 promoter (i.e., -2696T>C[rs7959801], -2297T>C[rs3742038] and -1081C>T[rs34570155]) and lung cancer risk.

Results:

We found that the SNP rs7959801T>C was significantly associated with lung cancer susceptibility. Compared to those with rs7959801TT wild-genotype, individuals with CT/CC variant genotypes exerted consistently beneficial roles in lung cancer risk in the discovery set (adjusted odd ratios [OR] = 0.67; 95% confidence interval [CI] = 0.57-0.80), and in the validation set (OR=0.69; 95%CI=0.54-0.88). Functional assays indicated that the allele transformation from T to C in rs7959801 of MSI1 gene arrestingly decreased its transcription activity in vitro. Furthermore, the expression levels of MSI1 were significantly lower in the patients with CT/CC variants than in those who were with TT genotype.

Conclusion:

Our findings suggested that the rs7959801T>C polymorphism in the MSI1 promoter conferred a decreased risk to lung cancer by reducing the expression of MSI1 and it may be a promising indicator for lung cancer predisposition.

Thirty loci identified for heart rate response to exercise and recovery implicate autonomic nervous system

Impaired capacity to increase heart rate (HR) during exercise (ΔHR^ex), and a reduced rate of recovery post-exercise (ΔHR^rec) are associated with higher cardiovascular mortality rates. Currently, the genetic basis of both phenotypes remains to be elucidated. We conduct genome-wide association studies (GWASs) for ΔHR^ex and ΔHR^rec in ~40,000 individuals, followed by replication in ~27,000 independent samples, all from UK Biobank. Six and seven single-nucleotide polymorphisms for ΔHR^ex and ΔHR^rec, respectively, formally replicate. In a full data set GWAS, eight further loci for ΔHR^ex and nine for ΔHR^rec are genome-wide significant (P ≤ 5 × 10⁻⁸). In total, 30 loci are discovered, 8 being common across traits. Processes of neural development and modulation of adrenergic activity by the autonomic nervous system are enriched in these results. Our findings reinforce current understanding of HR response to exercise and recovery and could guide future studies evaluating its contribution to cardiovascular risk prediction.

Genome-wide association studies have identified multiple loci for resting heart rate (HR) but the genetic factors associated with HR increase during and HR recovery after exercise are less well studied. Here, the authors examine both traits in a two-stage GWAS design in up to 67,257 individuals from UK Biobank.

Observational and Genetic Associations of Resting Heart Rate With Aortic Valve Calcium

It is unknown if lifelong exposure to increased hemodynamic stress from an elevated resting heart rate (HR) may contribute to aortic valve calcium (AVC). We performed multivariate regression analyses using data from 1,266 Framingham Heart Study (FHS) Offspring cohort participants and 6,764 Multi-Ethnic Study of Atherosclerosis (MESA) participants. We constructed a genetic risk score (GRS) for HR using summary-level data in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) AVC Consortium to investigate if there was evidence in favor of a causal relation. AVC was present in 39% of FHS Offspring cohort participants and in 13% of MESA cohort participants. In multivariate adjusted models, participants in the highest resting HR quartiles had significantly greater prevalence of AVC, with a prevalence ratio of 1.19 (95% confidence interval [CI] 0.99 to 1.44) for the FHS Offspring cohort and 1.32 (95% CI 1.12 to 1.63) for the MESA cohort, compared with those in the lowest quartile. There was a similar increase in the prevalence of AVC per standard deviation increase in resting HR in both FHS Offspring (prevalence ratio 1.08, 95% CI 1.01 to 1.15) and MESA (1.10, 95% CI 1.03 to 1.17). In contrast with these observational findings, a HR associated GRS was not significantly associated with AVC. Although our observational analysis indicates that a higher resting HR is associated with AVC, our genetic results do not support a causal relation. Unmeasured environmental and/or lifestyle factors associated with both increased resting HR and AVC that are not fully explained by covariates in our observational models may account for the association between resting HR and AVC.

Mutation of the Na+/K+-ATPase Atp1a1a.1 causes QT interval prolongation and bradycardia in zebrafish

The genetic underpinnings that orchestrate the vertebrate heart rate are not fully understood yet, but of high clinical importance, since diseases of cardiac impulse formation and propagation are common and severe human arrhythmias. To identify novel regulators of the vertebrate heart rate, we deciphered the pathogenesis of the bradycardia in the homozygous zebrafish mutant hiphop (hip) and identified a missense-mutation (N851K) in Na⁺/K⁺-ATPase α1-subunit (atp1a1a.1). N851K affects zebrafish Na⁺/K⁺-ATPase ion transport capacity, as revealed by in vitro pump current measurements. Inhibition of the Na⁺/K⁺-ATPase in vivo indicates that hip rather acts as a hypomorph than being a null allele. Consequently, reduced Na⁺/K⁺-ATPase function leads to prolonged QT interval and refractoriness in the hip mutant heart, as shown by electrocardiogram and in vivo electrical stimulation experiments. We here demonstrate for the first time that Na⁺/K⁺-ATPase plays an essential role in heart rate regulation by prolonging myocardial repolarization.

Big data from electronic health records for early and late translational cardiovascular research: challenges and potential

Aims

Cohorts of millions of people's health records, whole genome sequencing, imaging, sensor, societal and publicly available data present a rapidly expanding digital trace of health. We aimed to critically review, for the first time, the challenges and potential of big data across early and late stages of translational cardiovascular disease research.

Methods and results

We sought exemplars based on literature reviews and expertise across the BigData@Heart Consortium. We identified formidable challenges including: data quality, knowing what data exist, the legal and ethical framework for their use, data sharing, building and maintaining public trust, developing standards for defining disease, developing tools for scalable, replicable science and equipping the clinical and scientific work force with new inter-disciplinary skills. Opportunities claimed for big health record data include: richer profiles of health and disease from birth to death and from the molecular to the societal scale; accelerated understanding of disease causation and progression, discovery of new mechanisms and treatment-relevant disease sub-phenotypes, understanding health and diseases in whole populations and whole health systems and returning actionable feedback loops to improve (and potentially disrupt) existing models of research and care, with greater efficiency. In early translational research we identified exemplars including: discovery of fundamental biological processes e.g. linking exome sequences to lifelong electronic health records (EHR) (e.g. human knockout experiments); drug development: genomic approaches to drug target validation; precision medicine: e.g. DNA integrated into hospital EHR for pre-emptive pharmacogenomics. In late translational research we identified exemplars including: learning health systems with outcome trials integrated into clinical care; citizen driven health with 24/7 multi-parameter patient monitoring to improve outcomes and population-based linkages of multiple EHR sources for higher resolution clinical epidemiology and public health.

Conclusion

High volumes of inherently diverse ('big') EHR data are beginning to disrupt the nature of cardiovascular research and care. Such big data have the potential to improve our understanding of disease causation and classification relevant for early translation and to contribute actionable analytics to improve health and healthcare.

Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

Analogy between non-alcoholic steatohepatitis (NASH) and hypertension: a stepwise patient-tailored approach for NASH treatment

Non-alcoholic steatohepatitis (NASH) is a common liver disorder worldwide. Although there has been improvement in our understanding of the natural history and pathogenesis of the disease, there is still no approved therapy for NASH. NASH shares many similarities with primary hypertension, in that both are extremely common disorders that can easily lead to serious complications if left untreated. Both conditions are viewed as "silent killers", because the disease can progress over a period of time prior to the occurrence of potentially deadly outcomes. While attempts to find the "miracle pill" for NASH are unrealistic, we can make an analogy with the "stepwise combination" approach developed over the last few decades for the treatment of hypertension. In the present review, we summarize some of the similarities in the concepts that underlie NASH and hypertension. The development of a stepwise patient-tailored method for the treatment of NASH is presented.

Heart rate as a predictor of cardiovascular risk

BACKGROUND

Heart rate (HR) is a predictor of cardiovascular, cerebrovascular and all-cause mortality in the general population, as well as in patients with cardio- and cerebrovascular diseases. We aimed to summarize current knowledge regarding the influence of HR on cardio- and cerebrovascular morbidity and mortality.

MATERIALS AND METHODS

PubMed, MEDLINE, Ovid and EMBASE databases were searched for large follow-up studies or meta-analysis published between January 1990 and September 2017 in the English language using the following keyword "heart rate," "resting heart rate," "mortality," "outcome," "hypertension," "heart failure," "ischaemic heart disease," "coronary heart disease" and "stroke."

RESULTS

The relationship between increased HR and cardio- and cerebrovascular morbidity and mortality has been reported in a large number of studies, and the results regarding this association are concurrent. This connection is generally stronger in men than in women. The increase in HR usually occurs in parallel with elevation of blood pressure and metabolic disturbances (insulin resistance, dyslipidaemia). However, even after adjustment for the most important cardiovascular risk factors, HR remained an independent predictor of adverse events in global population or in patients with cardio- and cerebrovascular diseases.

CONCLUSION

HR has an important negative effect on cardio- and cerebrovascular morbidity and mortality. Future longitudinal investigations should clarify HR significance and optimal HR reduction for primary and secondary prevention in cardio- and cerebrovascular events.

Genetic study links components of the autonomous nervous system to heart-rate profile during exercise

Heart rate (HR) responds to exercise by increasing during exercise and recovering after exercise. As such, HR is an important predictor of mortality that researchers believe is modulated by the autonomic nervous system. However, the mechanistic basis underlying inter-individual differences has yet to be explained. Here, we perform a large-scale genome-wide analysis of HR increase and HR recovery in 58,818 UK Biobank individuals. Twenty-five independent SNPs in 23 loci are identified to be associated (p < 8.3 × 10⁻⁹) with HR increase or HR recovery. A total of 36 candidate causal genes are prioritized that are enriched for pathways related to neuron biology. No evidence is found of a causal relationship with mortality or cardiovascular diseases. However, a nominal association with parental lifespan requires further study. In conclusion, the findings provide new biological and clinical insight into the mechanistic underpinnings of HR response to exercise. The results also underscore the role of the autonomous nervous system in HR recovery.

Response of the heart rate (HR) to exercise is associated with cardiac fitness and risk of cardiac death. Here, in a genome-wide association study, Verweij et al. identify 23 loci for HR increase during exercise or HR recovery, and highlight pleiotropy with blood pressure by polygenic risk score analysis.

Genome-Wide Association Approach Identified Novel Genetic Predictors of Heart Rate Response to β-Blockers

Background

For many indications, the negative chronotropic effect of β‐blockers is important to their efficacy, yet the heart rate (HR) response to β‐blockers varies. Herein, we sought to use a genome‐wide association approach to identify novel single nucleotide polymorphisms (SNPs) associated with HR response to β‐blockers.

Methods and Results

We first performed 4 genome‐wide association analyses for HR response to atenolol (a β1‐adrenergic receptor blocker) as: (1) monotherapy or (2) add‐on therapy, in 426 whites and 273 blacks separately from the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) study. A meta‐analysis was then performed between the genome‐wide association analysis performed in PEAR atenolol monotherapy and add‐on therapy, in each race separately, using the inverse variance method assuming fixed effects. From this analysis, SNPs associated with HR response to atenolol at a P<1E‐05 were tested for replication in whites (n=200) and blacks (n=168) treated with metoprolol (a β1‐adrenergic receptor blocker). From the genome‐wide association meta‐analyses, SNP rs17117817 near olfactory receptor family10 subfamily‐p‐member1 (OR10P1), and SNP rs2364349 in sorting nexin‐9 (SNX9) replicated in blacks. The combined studies meta‐analysis P values for the rs17117817 and rs2364349 reached genome‐wide significance (rs17117817G‐allele; Meta‐β=5.53 beats per minute, Meta‐P=2E‐09 and rs2364349 A‐allele; Meta‐β=3.5 beats per minute, Meta‐P=1E‐08). Additionally, SNPs in the OR10P1 and SNX9 gene regions were also associated with HR response in whites.

Conclusions

This study highlights OR10P1 and SNX9 as novel genes associated with changes in HR in response to β‐blockers.

Clinical Trial Registration

URL: http://www.clinicaltrials.gov. Unique identifier: NCT00246519.

Differential effects of inhibitory G protein isoforms on G protein-gated inwardly rectifying K+ currents in adult murine atria

G protein-gated inwardly rectifying K⁺ (GIRK) channels are the major inwardly rectifying K⁺ currents in cardiac atrial myocytes and an important determinant of atrial electrophysiology. Inhibitory G protein α-subunits can both mediate activation via acetylcholine but can also suppress basal currents in the absence of agonist. We studied this phenomenon using whole cell patch clamping in murine atria from mice with global genetic deletion of Gα_i2, combined deletion of Gα_i1/Gα_i3, and littermate controls. We found that mice with deletion of Gα_i2 had increased basal and agonist-activated currents, particularly in the right atria while in contrast those with Gα_i1/Gα_i3 deletion had reduced currents. Mice with global genetic deletion of Gα_i2 had decreased action potential duration. Tissue preparations of the left atria studied with a multielectrode array from Gα_i2 knockout mice showed a shorter effective refractory period, with no change in conduction velocity, than littermate controls. Transcriptional studies revealed increased expression of GIRK channel subunit genes in Gα_i2 knockout mice. Thus different G protein isoforms have differential effects on GIRK channel behavior and paradoxically Gα_i2 act to increase basal and agonist-activated GIRK currents. Deletion of Gα_i2 is potentially proarrhythmic in the atria.

Fifteen Genetic Loci Associated With the Electrocardiographic P Wave

BACKGROUND

The P wave on an ECG is a measure of atrial electric function, and its characteristics may serve as predictors for atrial arrhythmias. Increased mean P-wave duration and P-wave terminal force traditionally have been used as markers for left atrial enlargement, and both have been associated with increased risk of atrial fibrillation. Here, we explore the genetic basis of P-wave morphology through meta-analysis of genome-wide association study results for P-wave duration and P-wave terminal force from 12 cohort studies.

METHODS AND RESULTS

We included 44 456 individuals, of which 6778 (16%) were of African ancestry. Genotyping, imputation, and genome-wide association study were performed at each study site. Summary-level results were meta-analyzed centrally using inverse-variance weighting. In meta-analyses of P-wave duration, we identified 6 significant (P<5×10^-8) novel loci and replicated a prior association with SCN10A. We identified 3 loci at SCN5A, TBX5, and CAV1/CAV2 that were jointly associated with the PR interval, PR segment, and P-wave duration. We identified 6 novel loci in meta-analysis of P-wave terminal force. Four of the identified genetic loci were significantly associated with gene expression in 329 left atrial samples. Finally, we observed that some of the loci associated with the P wave were linked to overall atrial conduction, whereas others identified distinct phases of atrial conduction.

CONCLUSIONS

We have identified 6 novel genetic loci associated with P-wave duration and 6 novel loci associated with P-wave terminal force. Future studies of these loci may aid in identifying new targets for drugs that may modify atrial conduction or treat atrial arrhythmias.

If Channel as an Emerging Therapeutic Target for Cardiovascular Diseases: A Review of Current Evidence and Controversies

In 2015, non-communicable diseases accounted for 39.5 million (70%) of the total 56.4 million deaths that occurred globally, of which 17.7 million (45%) were due to cardiovascular diseases. An elevated heart rate is considered to be one of the independent predictors and markers of future cardiovascular diseases. A variety of experimental and epidemiological studies have found that atherosclerosis, heart failure, coronary artery disease, stroke, and arrhythmia are linked to elevated heart rate. Although there are established drugs to reduce the heart rate, these drugs have undesirable side effects. Hence, the development of new drugs that selectively inhibit the heart rate is considered necessary. In the search for such drugs, almost four decades ago the I_f channel, also known as the “funny channel,” emerged as a novel site for the selective inhibition of heart rate. These I_f channels, with a mixed sodium and potassium inward current, have been identified in the sinoatrial node of the heart, which mediates the slow diastolic depolarization of the pacemaker of the spontaneous rhythmic cells. The hyperpolarization-activated cyclic nucleotide-gated (HCN) subfamily is primarily articulated in the heart and neurons that are encoded by a family of four genes (HCN1-4) and they identify the funny channel. Of these, HCN-4 is the principal protein in the sinoatrial node. Currently, funny channel inhibition is being targeted for the treatment and prevention of cardiovascular diseases such as atherosclerosis and stroke. A selective I_f channel inhibitor named ivabradine was discovered for clinical use in treating heart failure and coronary artery disease. However, inconsistencies regarding the clinical effects of ivabradine have been reported in the literature, suggesting the need for a rigorous analysis of the available evidence. The objective of this review is therefore to assess the current advances in targeting the I_f channel associated with ivabradine and related challenges.

Genome-wide association study of heart rate and its variability in Hispanic/Latino cohorts

BACKGROUND

Although time-domain measures of heart rate variability (HRV) are used to estimate cardiac autonomic tone and disease risk in multiethnic populations, the genetic epidemiology of HRV in Hispanics/Latinos has not been characterized.

OBJECTIVE

The purpose of this study was to conduct a genome-wide association study of heart rate (HR) and its variability in the Hispanic Community Health Study/Study of Latinos, Multi-Ethnic Study of Atherosclerosis, and Women's Health Initiative Hispanic SNP-Health Association Resource project (n = 13,767).

METHODS

We estimated HR (bpm), standard deviation of normal-to-normal interbeat intervals (SDNN, ms), and root mean squared difference in successive, normal-to-normal interbeat intervals (RMSSD, ms) from resting, standard 12-lead ECGs. We estimated associations between each phenotype and 17 million genotyped or imputed single nucleotide polymorphisms (SNPs), accounting for relatedness and adjusting for age, sex, study site, and ancestry. Cohort-specific estimates were combined using fixed-effects, inverse-variance meta-analysis. We investigated replication for select SNPs exceeding genome-wide (P <5 × 10-8) or suggestive (P <10-6) significance thresholds.

RESULTS

Two genome-wide significant SNPs replicated in a European ancestry cohort, 1 one for RMSSD (rs4963772; chromosome 12) and another for SDNN (rs12982903; chromosome 19). A suggestive SNP for HR (rs236352; chromosome 6) replicated in an African-American cohort. Functional annotation of replicated SNPs in cardiac and neuronal tissues identified potentially causal variants and mechanisms.

CONCLUSION

This first genome-wide association study of HRV and HR in Hispanics/Latinos underscores the potential for even modestly sized samples of non-European ancestry to inform the genetic epidemiology of complex traits.