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

Cumulative effects of common genetic variants on risk of sudden cardiac death

BACKGROUND

Genome-wide association studies and candidate-gene based approaches have identified multiple common variants associated with increased risk of sudden cardiac death (SCD). However, the independent contribution of these individual loci to disease risk is modest.

OBJECTIVE

To investigate the cumulative effects of genetic variants previously associated with SCD risk.

METHODS

A total of 966 SCD cases from the Oregon-Sudden Unexpected Death Study and 1,926 coronary artery disease controls from the Wellcome Trust Case-Control Consortium were investigated. We generated genetic risk scores (GRS) for each trait composed of variants previously associated with SCD or with abnormalities in specific electrocardiographic traits such as QRS duration, QTc interval and heart rate. GRSs were calculated using a weighted approach based on the number of risk alleles weighted by the beta coefficients derived from the original studies. We also compared the highest and lowest quintiles for the GRS composed of SCD SNPs.

RESULTS

Increased cumulative risk was observed for a GRS composed of 14 SCD-SNPs (OR=1.17 [1.05-1.29], P = 0.002). The risk for SCD was 1.5 fold higher in the highest quintile when compared to the lowest quintile (OR = 1.46[1.11-1.92]). We did not observe significant associations with SCD for SNPs that determine electrocardiographic traits.

CONCLUSIONS

A modest but significant effect on SCD risk was identified for a GRS composed of 14 previously associated SCD SNPs. While next generation sequencing methodology will continue to identify additional novel variants, these findings represent proof of concept for the additive effects of gene variants on SCD risk.

Translational regulation shapes the molecular landscape of complex disease phenotypes

The extent of translational control of gene expression in mammalian tissues remains largely unknown. Here we perform genome-wide RNA sequencing and ribosome profiling in heart and liver tissues to investigate strain-specific translational regulation in the spontaneously hypertensive rat (SHR/Ola). For the most part, transcriptional variation is equally apparent at the translational level and there is limited evidence of translational buffering. Remarkably, we observe hundreds of strain-specific differences in translation, almost doubling the number of differentially expressed genes. The integration of genetic, transcriptional and translational data sets reveals distinct signatures in 3'UTR variation, RNA-binding protein motifs and miRNA expression associated with translational regulation of gene expression. We show that a large number of genes associated with heart and liver traits in human genome-wide association studies are primarily translationally regulated. Capturing interindividual differences in the translated genome will lead to new insights into the genes and regulatory pathways underlying disease phenotypes.

Genetic variation at the human connexin 43 locus but not at the connexin 40 locus is associated with left bundle branch block

Background

Bundle branch block (BBB) has been regarded as a disease of the conduction system, but occurs in mice lacking connexin 40 (expressed in atria, proximal conduction system) or connexin 43 (expressed in Purkinje cells, cardiomyocytes).

Objective

The aim of this paper is to explore whether BBB is heritable, and whether polymorphisms at connexin 40 and connexin 43 loci are associated with BBB.

Methods

To assess BBB heritability, we screened descendants of men with BBB in the population cohort ‘The Study of Men Born 1913’. DNA samples from 80-year-old men with extreme QRS-duration phenotypes were used to search for polymorphisms at connexin 40 and 43 loci. Associations between identified polymorphisms and BBB were evaluated in an independent cohort (INTERGENE).

Results

Seventy-seven men from ‘The Study of Men Born 1913’ with BBB had 116 descendants. Among the 76 participating descendants, 2 sons (6.4%) had BBB at 54 years of age. At the same age, 0.9% of men born in 1913 had BBB. We identified 6 single nucleotide polymorphisms (SNPs) in connexin 40 and 1 polymorphism in connexin 43. In the INTERGENE cohort, the connexin 43 polymorphism was associated with left BBB (LBBB) (4 of 35 LBBB vs 16 of 232 without BBB, χ²=7.4, p=0.03), but not with right BBB (RBBB) or overall BBB. None of the connexin 40 SNPs or haplotypes were associated with LBBB or RBBB.

Conclusions

These findings indicate that conduction by connexin 43 within the ventricular muscle distal to the specialised conduction system may be important for LBBB development.

Estrogen-related receptor α (ERRα) and ERRγ are essential coordinators of cardiac metabolism and function

Almost all cellular functions are powered by a continuous energy supply derived from cellular metabolism. However, it is little understood how cellular energy production is coordinated with diverse energy-consuming cellular functions. Here, using the cardiac muscle system, we demonstrate that nuclear receptors estrogen-related receptor α (ERRα) and ERRγ are essential transcriptional coordinators of cardiac energy production and consumption. On the one hand, ERRα and ERRγ together are vital for intact cardiomyocyte metabolism by directly controlling expression of genes important for mitochondrial functions and dynamics. On the other hand, ERRα and ERRγ influence major cardiomyocyte energy consumption functions through direct transcriptional regulation of key contraction, calcium homeostasis, and conduction genes. Mice lacking both ERRα and cardiac ERRγ develop severe bradycardia, lethal cardiomyopathy, and heart failure featuring metabolic, contractile, and conduction dysfunctions. These results illustrate that the ERR transcriptional pathway is essential to couple cellular energy metabolism with energy consumption processes in order to maintain normal cardiac function.

Drosophila and experimental neurology in the post-genomic era

For decades, the fruit fly, Drosophila melanogaster, has been among the premiere genetic model systems for probing fundamental neurobiology, including elucidation of mechanisms responsible for human neurologic disorders. Flies continue to offer virtually unparalleled versatility and speed for genetic manipulation, strong genomic conservation, and a nervous system that recapitulates a range of cellular and network properties relevant to human disease. I focus here on four critical challenges emerging from recent advances in our understanding of the genomic basis of human neurologic disorders where innovative experimental strategies are urgently needed: (1) pinpointing causal genes from associated genomic loci; (2) confirming the functional impact of allelic variants; (3) elucidating nervous system roles for novel or poorly studied genes; and (4) probing network interactions within implicated regulatory pathways. Drosophila genetic approaches are ideally suited to address each of these potential translational roadblocks, and will therefore contribute to mechanistic insights and potential breakthrough therapies for complex genetic disorders in the coming years. Strategic collaboration between neurologists, human geneticists, and the Drosophila research community holds great promise to accelerate progress in the post-genomic era.

Cardiac deficiency of single cytochrome oxidase assembly factor scox induces p53-dependent apoptosis in a Drosophila cardiomyopathy model

The heart is a muscle with high energy demands. Hence, most patients with mitochondrial disease produced by defects in the oxidative phosphorylation (OXPHOS) system are susceptible to cardiac involvement. The presentation of mitochondrial cardiomyopathy includes hypertrophic, dilated and left ventricular noncompaction, but the molecular mechanisms involved in cardiac impairment are unknown. One of the most frequent OXPHOS defects in humans frequently associated with cardiomyopathy is cytochrome c oxidase (COX) deficiency caused by mutations in COX assembly factors such as Sco1 and Sco2. To investigate the molecular mechanisms that underlie the cardiomyopathy associated with Sco deficiency, we have heart specifically interfered scox expression, the single Drosophila Sco orthologue. Cardiac-specific knockdown of scox reduces fly lifespan, and it severely compromises heart function and structure, producing dilated cardiomyopathy. Cardiomyocytes with low levels of scox have a significant reduction in COX activity and they undergo a metabolic switch from OXPHOS to glycolysis, mimicking the clinical features found in patients harbouring Sco mutations. The major cardiac defects observed are produced by a significant increase in apoptosis, which is dp53-dependent. Genetic and molecular evidence strongly suggest that dp53 is directly involved in the development of the cardiomyopathy induced by scox deficiency. Remarkably, apoptosis is enhanced in the muscle and liver of Sco2 knock-out mice, clearly suggesting that cell death is a key feature of the COX deficiencies produced by mutations in Sco genes in humans.

CASP3 genetic variants and susceptibility to atrial fibrillation in Chinese Han population

BACKGROUND

Caspase-3 plays an important role in the initiation and maintenance of atrial fibrillation (AF), but little is known about the role of CASP3 variants in the susceptibility to atrial fibrillation (AF). The purpose of this study was to comprehensively investigate the association between common genetic variants of CASP3 gene and AF in Chinese Han population.

METHODS AND RESULTS

We investigated the association of five variants in CASP3 and the risk of AF in 889 AF patients and 1015 controls. The genotype distribution of the rs4647602 was significantly different between patients with AF and controls (p<0.001), and the A allele frequency was significantly higher in AF patients than in controls (61.0% vs 53.2%; p<0.001). Compared with CC genotype carriers, subjects with AA genotype had significantly increased susceptibility to AF (OR=1.84, p<0.001). Multivariable logistic regression analysis showed that the rs4647602 polymorphism was significantly associated with risk of AF under dominant, recessive and additive genetic model (OR, 1.44-1.64; all p<0.001). There was no association between the other four SNPs (rs6948, rs2696056, rs4647602 and rs4647610) and risk of AF.

CONCLUSION

The rs4647602 polymorphism is independently associated with the risk of AF in Chinese Han population.

Genetic contribution to the variance of blood pressure and heart rate: a systematic review and meta-regression of twin studies

The genetic contribution of blood pressure and heart rate (HR) varied widely between studies. Demographic factors such as ethnicity, age and/or sex might explain some of the heterogeneity. We performed a systematic review focusing on four phenotypes: systolic blood pressure (SBP), diastolic blood pressure (DBP), HR and pulse pressure (PP). Meta-regression was conducted to analyze potential factors in relation to SBP and DBP heritability. A total of 10,613 independent twins that came from 17 studies were included in the analysis. The weighted mean value of heritability for SBP and DBP was 0.54 (95% CIs: 0.48-0.60) and 0.49 (95% CIs: 0.42-0.56). Comparatively, three studies of HR and four studies of PP heritability were limited for the heterogeneity test. Meta-regression showed that, on average, SBP heritability with additive genes/unique environment (AE) model tend to have a higher heritability than additive genes/shared environment/unique environment (ACE) model (coefficient = 0.0947, p = .0142). A similar result was found for DBP as well. No other factors such as sex, age, ethnicity, publication year were significantly associated with heritability variance. Our study shows heritability estimates based on twin studies of both SBP and DBP are around 50%, using an AE rather than an ACE model; the variance due to C ended up in A, suggesting that the AE model may overestimate heritability if a small contribution of shared environment exists.

New approaches to establish genetic causality

Cardiovascular medicine has evolved rapidly in the era of genomics with many diseases having primary genetic origins becoming the subject of intense investigation. The resulting avalanche of information on the molecular causes of these disorders has prompted a revolution in our understanding of disease mechanisms and provided new avenues for diagnoses. At the heart of this revolution is the need to correctly classify genetic variants discovered during the course of research or reported from clinical genetic testing. This review will address current concepts related to establishing the cause and effect relationship between genomic variants and heart diseases. A survey of general approaches used for functional annotation of variants will also be presented.

Fruit flies in biomedical research

Many scientists complain that the current funding situation is dire. Indeed, there has been an overall decline in support in funding for research from the National Institutes of Health and the National Science Foundation. Within the Drosophila field, some of us question how long this funding crunch will last as it demotivates principal investigators and perhaps more importantly affects the long-term career choice of many young scientists. Yet numerous very interesting biological processes and avenues remain to be investigated in Drosophila, and probing questions can be answered fast and efficiently in flies to reveal new biological phenomena. Moreover, Drosophila is an excellent model organism for studies that have translational impact for genetic disease and for other medical implications such as vector-borne illnesses. We would like to promote a better collaboration between Drosophila geneticists/biologists and human geneticists/bioinformaticians/clinicians, as it would benefit both fields and significantly impact the research on human diseases.

Novel mutation in the α-myosin heavy chain gene is associated with sick sinus syndrome

BACKGROUND

Recent genome-wide association studies have demonstrated an association between MYH6, the gene encoding α-myosin heavy chain (α-MHC), and sinus node function in the general population. Moreover, a rare MYH6 variant, R721W, predisposing susceptibility to sick sinus syndrome has been identified. However, the existence of disease-causing MYH6 mutations for familial sick sinus syndrome and their underlying mechanisms remain unknown.

METHODS AND RESULTS

We screened 9 genotype-negative probands with sick sinus syndrome families for mutations in MYH6 and identified an in-frame 3-bp deletion predicted to delete one residue (delE933) at the highly conserved coiled-coil structure within the binding motif to myosin-binding protein C in one patient. Co-immunoprecipitation analysis revealed enhanced binding of delE933 α-MHC to myosin-binding protein C. Irregular fluorescent speckles retained in the cytoplasm with substantially disrupted sarcomere striation were observed in neonatal rat cardiomyocytes transfected with α-MHC mutants carrying delE933 or R721W. In addition to the sarcomere impairments, delE933 α-MHC exhibited electrophysiological abnormalities both in vitro and in vivo. The atrial cardiomyocyte cell line HL-1 stably expressing delE933 α-MHC showed a significantly slower conduction velocity on multielectrode array than those of wild-type α-MHC or control plasmid transfected cells. Furthermore, targeted morpholino knockdown of MYH6 in zebrafish significantly reduced the heart rate, which was rescued by coexpressed wild-type human α-MHC but not by delE933 α-MHC.

CONCLUSIONS

The novel MYH6 mutation delE933 causes both structural damage of the sarcomere and functional impairments on atrial action propagation. This report reinforces the relevance of MYH6 for sinus node function and identifies a novel pathophysiology underlying familial sick sinus syndrome.

Loss of dihydrolipoyl succinyltransferase (DLST) leads to reduced resting heart rate in the zebrafish

The genetic underpinnings of heart rate regulation are only poorly understood. In search for genetic regulators of cardiac pacemaker activity, we isolated in a large-scale mutagenesis screen the embryonic lethal, recessive zebrafish mutant schneckentempo (ste). Homozygous ste mutants exhibit a severely reduced resting heart rate with normal atrio-ventricular conduction and contractile function. External electrical pacing reveals that defective excitation generation in cardiac pacemaker cells underlies bradycardia in ste^−/− mutants. By positional cloning and gene knock-down analysis we find that loss of dihydrolipoyl succinyltransferase (DLST) function causes the ste phenotype. The mitochondrial enzyme DLST is an essential player in the citric acid cycle that warrants proper adenosine-tri-phosphate (ATP) production. Accordingly, ATP levels are significantly diminished in ste^−/− mutant embryos, suggesting that limited energy supply accounts for reduced cardiac pacemaker activity in ste^−/− mutants. We demonstrate here for the first time that the mitochondrial enzyme DLST plays an essential role in the modulation of the vertebrate heart rate by controlling ATP production in the heart.

Pacemaker activity of the human sinoatrial node: an update on the effects of mutations in HCN4 on the hyperpolarization-activated current

Since 2003, several loss-of-function mutations in the HCN4 gene, which encodes the HCN4 protein, have been associated with sinus node dysfunction. In human sinoatrial node (SAN), HCN4 is the most abundant of the four isoforms of the HCN family. Tetramers of HCN subunits constitute the ion channels that conduct the hyperpolarization-activated "funny" current (If), which plays an important modulating role in SAN pacemaker activity. Voltage-clamp experiments on HCN4 channels expressed in COS-7, CHO and HEK-293 cells, as well as in Xenopus oocytes have revealed changes in the expression and kinetics of mutant channels, but the extent to which especially the kinetic changes would affect If flowing during a human SAN action potential often remains unresolved. In our contribution to the Topical Collection on Human Single Nucleotide Polymorphisms and Disease Diagnostics, we provide an updated review of the mutation-induced changes in the expression and kinetics of HCN4 channels and provide an overview of their effects on If during the time course of a human SAN action potential, as assessed in simulated action potential clamp experiments. Future research may solve apparent inconsistencies between data from clinical studies and data from in vitro and in silico experiments.

Animal and in silico models for the study of sarcomeric cardiomyopathies

Over the past decade, our understanding of cardiomyopathies has improved dramatically, due to improvements in screening and detection of gene defects in the human genome as well as a variety of novel animal models (mouse, zebrafish, and drosophila) and in silico computational models. These novel experimental tools have created a platform that is highly complementary to the naturally occurring cardiomyopathies in cats and dogs that had been available for some time. A fully integrative approach, which incorporates all these modalities, is likely required for significant steps forward in understanding the molecular underpinnings and pathogenesis of cardiomyopathies. Finally, novel technologies, including CRISPR/Cas9, which have already been proved to work in zebrafish, are currently being employed to engineer sarcomeric cardiomyopathy in larger animals, including pigs and non-human primates. In the mouse, the increased speed with which these techniques can be employed to engineer precise 'knock-in' models that previously took years to make via multiple rounds of homologous recombination-based gene targeting promises multiple and precise models of human cardiac disease for future study. Such novel genetically engineered animal models recapitulating human sarcomeric protein defects will help bridging the gap to translate therapeutic targets from small animal and in silico models to the human patient with sarcomeric cardiomyopathy.

Expression profile analysis of circulating microRNAs and their effects on ion channels in Chinese atrial fibrillation patients

OBJECTIVE

To investigate the changes in expression profile of circulating microRNAs (miRNAs) and the regulatory effect of atrial fibrilation (AF)-related miRNAs on ion channels.

METHODS

112 patients with AF were assigned into observation group, and another 112 non-AF people were assigned into control group. Total plasma RNAs were extracted from patients' blood samples. Differentially expressed miRNA-1s were transfected into primary-cultured neonatal rat cardiac myocytes.

RESULTS

Compared with control group, significant differences were observed in 15 kinds of miRNAs in observation group. Down-regulation of the expression of miRNAs included hsa-miR-328, hsa-miR-145, hsa-miR-222, hsa-miR-1, hsa-miR-162, hsa-miR-432, and hsa-miR-493b; Up-regulation of the expression included hsa-miR634, hsa-miR-664, hsa-miR-9, hsa-miR-152, hsa-miR-19, hsa-miR-454, hsa-miR-146, and hsa-miR-374a. The expression level of CACNB2 protein in miRNA-1 group was significantly lower than that in blank control group, negative control group, MTmiRNA-1 group, AMO-1 group and miRNA-1+AMO-1 cotransfection group (P < 0.05), while in AMO-1 group, the expression level of CACNB2 protein was significantly higher than that in other groups (P < 0.05). These results indicated that transfected miRNA-1 could significantly inhibit the expression of CACNB2 protein.

CONCLUSIONS

Circulating miRNAs can be used in studies concerning on the regulation mechanism of the occurrence and development of AF. MiRNA-1 can decrease the intracellular Ca(2+) concentration and prevent the AF.

Gain-of-function mutations in potassium channel subunit KCNE2 associated with early-onset lone atrial fibrillation

AIMS

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Disturbances in cardiac potassium conductance are considered as one of the disease mechanisms in AF. We aimed to investigate if mutations in potassium-channel β-subunits KCNE2 and KCNE3 are associated with early-onset lone AF.

METHODS & RESULTS

The coding regions of KCNE2 and KCNE3 were bidirectionally sequenced in 192 unrelated patients diagnosed with early-onset lone AF (<40 years). Two nonsynonymous missense mutations were identified in KCNE2 (M23L and I57T). Both mutations were absent in a healthy control group (n=1500 alleles). Electrophysiological investigations were performed for both mutations in combination with candidate pore-forming α-subunits KV7.1, KV11.1, KV4.3 and KV1.5. A significant gain-of-function effect was observed upon coexpression with KV7.1 and KV7.1+KCNE1. Confocal imaging found no differences in subcellular localization. No disease-suspected mutations were identified in KCNE3.

CONCLUSION

We identified two KCNE2 gain-of-function missense mutations that seem to increase the susceptibility of early-onset lone AF. These results confirm previous findings indicating that gain-of-function in the slow delayed rectifier potassium current might be involved in the pathogenesis of AF.

Prioritizing causal disease genes using unbiased genomic features

BACKGROUND

Cardiovascular disease (CVD) is the leading cause of death in the developed world. Human genetic studies, including genome-wide sequencing and SNP-array approaches, promise to reveal disease genes and mechanisms representing new therapeutic targets. In practice, however, identification of the actual genes contributing to disease pathogenesis has lagged behind identification of associated loci, thus limiting the clinical benefits.

RESULTS

To aid in localizing causal genes, we develop a machine learning approach, Objective Prioritization for Enhanced Novelty (OPEN), which quantitatively prioritizes gene-disease associations based on a diverse group of genomic features. This approach uses only unbiased predictive features and thus is not hampered by a preference towards previously well-characterized genes. We demonstrate success in identifying genetic determinants for CVD-related traits, including cholesterol levels, blood pressure, and conduction system and cardiomyopathy phenotypes. Using OPEN, we prioritize genes, including FLNC, for association with increased left ventricular diameter, which is a defining feature of a prevalent cardiovascular disorder, dilated cardiomyopathy or DCM. Using a zebrafish model, we experimentally validate FLNC and identify a novel FLNC splice-site mutation in a patient with severe DCM.

CONCLUSION

Our approach stands to assist interpretation of large-scale genetic studies without compromising their fundamentally unbiased nature.

Pacemaker activity of the human sinoatrial node: effects of HCN4 mutations on the hyperpolarization-activated current

The hyperpolarization-activated 'funny' current, If, plays an important modulating role in the pacemaker activity of the human sinoatrial node (SAN). If is carried by hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, which are tetramers built of four HCN subunits. In human SAN, HCN4 is the most abundant of the four isoforms of the HCN family. Since 2003, several loss-of-function mutations in the HCN4 gene, which encodes the HCN4 protein, or in the KCNE2 gene, which encodes the MiRP1 accessory β-subunit, have been associated with sinus node dysfunction. Voltage-clamp experiments on HCN4 channels expressed in COS-7 cells, Xenopus oocytes, or HEK-293 cells have revealed changes in the expression and kinetics of mutant channels, but the extent to which these changes would affect If flowing during a human SAN action potential is unresolved. Here, we review the changes in expression and kinetics of HCN4 mutant channels and provide an overview of their effects on If during the time course of a human SAN action potential, both under resting conditions and upon adrenergic stimulation. These effects are assessed in simulated action potential clamp experiments, with action potentials recorded from isolated human SAN pacemaker cells as command potential and kinetics of If based on voltage-clamp data from these cells. Results from in vitro and in silico experiments show several inconsistencies with clinical observations, pointing to challenges for future research.

Baf250a orchestrates an epigenetic pathway to repress the Nkx2.5-directed contractile cardiomyocyte program in the sinoatrial node

The sinoatrial node (SAN) is essential for rhythmic beating of the heart; however, our understanding of what controls proper functioning of the SAN remains primitive. To explore molecular control of SAN function, we specifically deleted Baf250a, a key regulatory component of the ATP-dependent chromatin remodeling complex SWI/SNF, in the SAN. Deletion of Baf250a in the SAN led to sinus bradycardia. Time series analysis of dysregulated genes after deletion of Baf250a reveals a transcriptional hierarchy maintaining pacemaker cell identity, i.e., Baf250a activates the expression of Tbx3, and Baf250a, Tbx3 and histone deacetylase 3 coordinately repress the expression of Nkx2.5. Disruption of this repressive pathway switches on expression of Nkx2.5, which stimulates expression of Gata4 and Tbx5. These three cardiac transcription factors further turn on a contractile cardiomyocyte program in the SAN, which eventually leads to sick sinus disease (SSD). Our study suggests that disruption of key genetic pathways regulating cardiac lineage segregation may cause SSD and cardiac arrhythmias in general.

A genome-wide association study to identify genomic modulators of rate control therapy in patients with atrial fibrillation

For many patients with atrial fibrillation, ventricular rate control with atrioventricular (AV) nodal blockers is considered first-line therapy, although response to treatment is highly variable. Using an extreme phenotype of failure of rate control necessitating AV nodal ablation and pacemaker implantation, we conducted a genome-wide association study (GWAS) to identify genomic modulators of rate control therapy. Cases included 95 patients who failed rate control therapy. Controls (n = 190) achieved adequate rate control therapy with ≤2 AV nodal blockers using a conventional clinical definition. Genotyping was performed on the Illumina 610-Quad platform, and results were imputed to the 1000 Genomes reference haplotypes. A total of 554,041 single-nucleotide polymorphisms (SNPs) met criteria for minor allele frequency (>0.01), call rate (>95%), and quality control, and 6,055,224 SNPs were available after imputation. No SNP reached the canonical threshold for significance for GWAS of p <5 × 10(-8). Sixty-three SNPs with p <10(-5) at 6 genomic loci were genotyped in a validation cohort of 130 cases and 157 controls. These included 6q24.3 (near SAMD5/SASH1, p = 9.36 × 10(-8)), 4q12 (IGFBP7, p = 1.75 × 10(-7)), 6q22.33 (C6orf174, p = 4.86 × 10(-7)), 3p21.31 (CDCP1, p = 1.18 × 10(-6)), 12p12.1 (SOX5, p = 1.62 × 10(-6)), and 7p11 (LANCL2, p = 6.51 × 10(-6)). However, none of these were significant in the replication cohort or in a meta-analysis of both cohorts. In conclusion, we identified several potentially important genomic modulators of rate control therapy in atrial fibrillation, particularly SOX5, which was previously associated with heart rate at rest and PR interval. However, these failed to reach genome-wide significance.

Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization

The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal Mendelian Long QT Syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals we identified 35 common variant QT interval loci, that collectively explain ∼8-10% of QT variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 novel QT loci in 298 unrelated LQTS probands identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode for proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies novel candidate genes for ventricular arrhythmias, LQTS,and SCD.

Annotation of loci from genome-wide association studies using tissue-specific quantitative interaction proteomics

Genome-wide association studies (GWAS) have identified thousands of loci associated wtih complex traits, but it is challenging to pinpoint causal genes in these loci and to exploit subtle association signals. We used tissue-specific quantitative interaction proteomics to map a network of five genes involved in the Mendelian disorder long QT syndrome (LQTS). We integrated the LQTS network with GWAS loci from the corresponding common complex trait, QT interval variation, to identify candidate genes that were subsequently confirmed in Xenopus laevis oocytes and zebrafish. We used the LQTS protein network to filter weak GWAS signals by identifying single nucleotide polymorphisms (SNPs) in proximity to genes in the network supported by strong proteomic evidence. Three SNPs passing this filter reached genome-wide significance after replication genotyping. Overall, we present a general strategy to propose candidates in GWAS loci for functional studies and to systematically filter subtle association signals using tissue-specific quantitative interaction proteomics.

Atrial fibrillation associated with ivabradine treatment: meta-analysis of randomised controlled trials

Objective

To quantify any risk of atrial fibrillation (AF) associated with ivabradine treatment by meta-analysis of clinical trial data.

Methods

Medline, Embase, Web of Knowledge and the Cochrane central register of controlled trials were searched for double-blinded randomised controlled trials of ivabradine with a minimum follow-up period of 4 weeks. For studies where AF data were unpublished, safety data were obtained from the European Medicines Agency (EMeA) website and personal communications. Studies were appraised for risk of bias using components recommended by the Cochrane Collaboration. Meta-analyses were performed of relative risk of AF and absolute risk difference of AF per year of treatment. The main outcome measure was incident AF during the follow-up period.

Results

AF data were available from 11 studies: one from the published report, six from the EMeA and four from personal communications. Ivabradine treatment was associated with a relative risk of AF of 1.15 (95% CI 1.07 to 1.24, p=0.0027) among 21 571 patients in the meta-analysis. From this we estimated that the number needed to harm for ivabradine would be 208 (95% CI 122 to 667) per year of treatment.

Conclusions

AF is a substantially more common side effect of ivabradine treatment than one patient in 10 000, the risk presently reported in the product literature. The incidence of AF has not routinely been reported in clinical trials of ivabradine.

Genome-wide identification of expression quantitative trait loci (eQTLs) in human heart

In recent years genome-wide association studies (GWAS) have uncovered numerous chromosomal loci associated with various electrocardiographic traits and cardiac arrhythmia predisposition. A considerable fraction of these loci lie within inter-genic regions. The underlying trait-associated variants likely reside in regulatory regions and exert their effect by modulating gene expression. Hence, the key to unraveling the molecular mechanisms underlying these cardiac traits is to interrogate variants for association with differential transcript abundance by expression quantitative trait locus (eQTL) analysis. In this study we conducted an eQTL analysis of human heart. For a total of 129 left ventricular samples that were collected from non-diseased human donor hearts, genome-wide transcript abundance and genotyping was determined using microarrays. Each of the 18,402 transcripts and 897,683 SNP genotypes that remained after pre-processing and stringent quality control were tested for eQTL effects. We identified 771 eQTLs, regulating 429 unique transcripts. Overlaying these eQTLs with cardiac GWAS loci identified novel candidates for studies aimed at elucidating the functional and transcriptional impact of these loci. Thus, this work provides for the first time a comprehensive eQTL map of human heart: a powerful and unique resource that enables systems genetics approaches for the study of cardiac traits.

Systems genetics: From GWAS to disease pathways

Most common diseases are complex, involving multiple genetic and environmental factors and their interactions. In the past decade, genome-wide association studies (GWAS) have successfully identified thousands of genetic variants underlying susceptibility to complex diseases. However, the results from these studies often do not provide evidence on how the variants affect downstream pathways and lead to the disease. Therefore, in the post-GWAS era the greatest challenge lies in combining GWAS findings with additional molecular data to functionally characterize the associations. The advances in various ~omics techniques have made it possible to investigate the effect of risk variants on intermediate molecular levels, such as gene expression, methylation, protein abundance or metabolite levels. As disease aetiology is complex, no single molecular analysis is expected to fully unravel the disease mechanism. Multiple molecular levels can interact and also show plasticity in different physiological conditions, cell types and disease stages. There is therefore a great need for new integrative approaches that can combine data from different molecular levels and can help construct the causal inference from genotype to phenotype. Systems genetics is such an approach; it is used to study genetic effects within the larger scope of systems biology by integrating genotype information with various ~omics datasets as well as with environmental and physiological variables. In this review, we describe this approach and discuss how it can help us unravel the molecular mechanisms through which genetic variation causes disease. This article is part of a Special Issue entitled: From Genome to Function.

Insight into genetic determinants of resting heart rate

BACKGROUND

Recent studies suggested that resting heart rate (RHR) might be an independent predictor of cardiovascular mortality and morbidity. Nonetheless, the interrelation between RHR and cardiovascular diseases is not clear. In order to resolve this puzzle, the importance of genetic determinants of RHR has been recently suggested, but it needs to be further investigated.

OBJECTIVE

The aim of this study was to estimate the contribution of common genetic variations on RHR using Genome Wide Association Study.

METHODS

We performed a Genome Wide Association Study in an isolated population cohort of 1737 individuals, the Italian Network on Genetic Isolates - Friuli Venezia Giulia (INGI-FVG). Moreover, a haplotype analysis was performed. A regression tree analysis was run to highlight the effect of each haplotype combination on the phenotype.

RESULTS

A significant level of association (p<5 × 10(-8)) was detected for Single Nucleotide Polymorphisms (SNPs) in two genes expressed in the heart: MAML1 and CANX. Founding that the three different variants of the haplotype, which encompass both genes, yielded a phenotypic correlation. Indeed, a haplotype in homozygosity is significantly associated with the lower quartile of RHR (RHR ≤ 58 bpm). Moreover no significant association was found between cardiovascular risk factors and the different haplotype combinations.

CONCLUSION

Mastermind-like 1 and Calnexin were found to be associated with RHR. We demonstrated a relation between a haplotype and the lower quartile of RHR in our populations. Our findings highlight that genetic determinants of RHR may be implicated in determining cardiovascular diseases and could allow a better risk stratification.

Cohort Profile: Estonian Biobank of the Estonian Genome Center, University of Tartu

The Estonian Biobank cohort is a volunteer-based sample of the Estonian resident adult population (aged ≥18 years). The current number of participants-close to 52000--represents a large proportion, 5%, of the Estonian adult population, making it ideally suited to population-based studies. General practitioners (GPs) and medical personnel in the special recruitment offices have recruited participants throughout the country. At baseline, the GPs performed a standardized health examination of the participants, who also donated blood samples for DNA, white blood cells and plasma tests and filled out a 16-module questionnaire on health-related topics such as lifestyle, diet and clinical diagnoses described in WHO ICD-10. A significant part of the cohort has whole genome sequencing (100), genome-wide single nucleotide polymorphism (SNP) array data (20 000) and/or NMR metabolome data (11 000) available (http://www.geenivaramu.ee/for-scientists/data-release/). The data are continuously updated through periodical linking to national electronic databases and registries. A part of the cohort has been re-contacted for follow-up purposes and resampling, and targeted invitations are possible for specific purposes, for example people with a specific diagnosis. The Estonian Genome Center of the University of Tartu is actively collaborating with many universities, research institutes and consortia and encourages fellow scientists worldwide to co-initiate new academic or industrial joint projects with us.

Atrial fibrillation and microRNAs

Atrial fibrillation (AF) is the most common sustained arrhythmia, especially in the elderly, and has a significant genetic component. Recently, several independent investigators have demonstrated a functional role for small non-coding RNAs (microRNAs) in the pathophysiology of this cardiac arrhythmia. This report represents a systematic and updated appraisal of the main studies that established a mechanistic association between specific microRNAs and AF, focusing both on the regulation of electrical and structural remodeling of cardiac tissue.

From GWAS to function: genetic variation in sodium channel gene enhancer influences electrical patterning

The electrical activity of the heart depends on the correct interplay between key transcription factors and cis-regulatory elements, which together regulate the proper heterogeneous expression of genes encoding for ion channels and other proteins. Genome-wide association studies of ECG parameters implicated genetic variants in the genes for these factors and ion channels modulating conduction and depolarization. Here, we review recent insights into the regulation of localized expression of ion channel genes and the mechanism by which a single-nucleotide polymorphism (SNP) associated with alterations in cardiac conduction patterns in humans affects the transcriptional regulation of the sodium channel genes, SCN5A and SCN10A. The identification of regulatory elements of electrical activity genes helps to explain the impact of genetic variants in non-coding regulatory DNA sequences on regulation of cardiac conduction and the predisposition for cardiac arrhythmias.

Genetics of sudden cardiac death caused by ventricular arrhythmias

Sudden cardiac death (SCD) resulting from ventricular tachyarrhythmia is a major contributor to mortality. Clinical management of SCD, currently based on clinical markers of SCD risk, can be improved by integrating genetic information. The identification of multiple disease-causing gene variants has already improved patient management and increased our understanding of the rare Mendelian diseases associated with SCD risk in the young, but marked variability in disease severity suggests that additional genetic modifiers exist. Next-generation DNA sequencing could be crucial to the discovery of SCD-associated genes, but large data sets can be difficult to interpret. SCD usually occurs in patients with an average age of 65 years who have complex cardiac disease stemming from multiple, common, acquired disorders. Heritable factors are largely unknown, but are likely to have a role in determining the risk of SCD in these patients. Numerous genetic loci have been identified that affect electrocardiogram indices, which are regarded as intermediate phenotypes for tachyarrhythmia. These loci could help to identify new molecules and pathways affecting cardiac electrical function. These loci are often located in intergenic regions, so our evolving understanding of the noncoding regulatory regions of the genome are likely to aid in the identification of novel genes that are important for cardiac electrical function and possibly SCD.

The Rotterdam Study: 2014 objectives and design update

The Rotterdam Study is a prospective cohort study ongoing since 1990 in the city of Rotterdam in The Netherlands. The study targets cardiovascular, endocrine, hepatic, neurological, ophthalmic, psychiatric, dermatological, oncological, and respiratory diseases. As of 2008, 14,926 subjects aged 45 years or over comprise the Rotterdam Study cohort. The findings of the Rotterdam Study have been presented in over a 1,000 research articles and reports (see www.erasmus-epidemiology.nl/rotterdamstudy ). This article gives the rationale of the study and its design. It also presents a summary of the major findings and an update of the objectives and methods.

Beyond GWASs: illuminating the dark road from association to function

Genome-wide association studies (GWASs) have enabled the discovery of common genetic variation contributing to normal and pathological traits and clinical drug responses, but recognizing the precise targets of these associations is now the major challenge. Here, we review recent approaches to the functional follow-up of GWAS loci, including fine mapping of GWAS signal(s), prioritization of putative functional SNPs by the integration of genetic epidemiological and bioinformatic methods, and in vitro and in vivo experimental verification of predicted molecular mechanisms for identifying the targeted genes. The majority of GWAS-identified variants fall in noncoding regions of the genome. Therefore, this review focuses on strategies for assessing likely mechanisms affected by noncoding variants; such mechanisms include transcriptional regulation, noncoding RNA function, and epigenetic regulation. These approaches have already accelerated progress from genetic studies to biological knowledge and might ultimately guide the development of prognostic, preventive, and therapeutic measures.

Activation of histone deacetylase-6 induces contractile dysfunction through derailment of α-tubulin proteostasis in experimental and human atrial fibrillation

BACKGROUND

Atrial fibrillation (AF) is characterized by structural remodeling, contractile dysfunction, and AF progression. Histone deacetylases (HDACs) influence acetylation of both histones and cytosolic proteins, thereby mediating epigenetic regulation and influencing cell proteostasis. Because the exact function of HDACs in AF is unknown, we investigated their role in experimental and clinical AF models.

METHODS AND RESULTS

Tachypacing of HL-1 atrial cardiomyocytes and Drosophila pupae hearts significantly impaired contractile function (amplitude of Ca(2+) transients and heart wall contractions). This dysfunction was prevented by inhibition of HDAC6 (tubacin) and sirtuins (nicotinamide). Tachypacing induced specific activation of HDAC6, resulting in α-tubulin deacetylation, depolymerization, and degradation by calpain. Tachypacing-induced contractile dysfunction was completely rescued by dominant-negative HDAC6 mutants with loss of deacetylase activity in the second catalytic domain, which bears α-tubulin deacetylase activity. Furthermore, in vivo treatment with the HDAC6 inhibitor tubastatin A protected atrial tachypaced dogs from electric remodeling (action potential duration shortening, L-type Ca(2+) current reduction, AF promotion) and cellular Ca(2+)-handling/contractile dysfunction (loss of Ca(2+) transient amplitude, sarcomere contractility). Finally, atrial tissue from patients with AF also showed a significant increase in HDAC6 activity and reduction in the expression of both acetylated and total α-tubulin.

CONCLUSIONS

AF induces remodeling and loss of contractile function, at least in part through HDAC6 activation and subsequent derailment of α-tubulin proteostasis and disruption of the cardiomyocyte microtubule structure. In vivo inhibition of HDAC6 protects against AF-related atrial remodeling, disclosing the potential of HDAC6 as a therapeutic target in clinical AF.

Human heart rate: heritability of resting and stress values in twin pairs, and influence of genetic variation in the adrenergic pathway at a microribonucleic acid (microrna) motif in the 3'-UTR of cytochrome b561 [corrected]

OBJECTIVES

The goal of this study was to understand the role of genetic variation in the catecholamine biosynthetic pathway for control of human heart rate (HR).

BACKGROUND

Human HR is an integrated cardiovascular trait predictive of morbidity and survival. Because the autonomic pathway exerts rapid control over the heart, we probed the role of heredity in the control of HR, focusing on a component of the autonomic sympathetic pathway already predictive of outflow responses: cytochrome b561 (CYB561), the electron shuttle in catecholamine vesicle membranes for transmitter biosynthesis.

METHODS

We studied hereditary control of HR with the twin pair design, at rest and during environmental (cold) stress. Single nucleotide polymorphism disruption of a microribonucleic acid (microRNA) recognition motif in the human CYB561 3'-UTR was identified computationally, and its differential effect on gene expression was demonstrated in a transfected luciferase reporter/3'-UTR variant. We exposed stem cell-derived human embryoid bodies to the microRNA mimic or antagomir oligonucleotides, and we observed the effects on contraction rate in proto-hearts.

RESULTS

Substantial heritability (h(2)) was demonstrated by using twin pair variance components for both basal/resting HR (h(2) 50.9 ± 6.4% of trait variation, p = 2.47 × 10(-10)) and stress-augmented HR (h(2) 55.1 ± 5.9%, p = 8.79 × 10(-13)), and the 2 HR traits shared genetic determination (genetic covariance ρG 0.747 ± 0.058, p = 2.85 × 10(-9)). CYB561 displayed 1 common genetic variant in the transcript region: A+1485G (rs3087776), in the 3'-UTR, 1485 bp downstream of the termination codon, in a conserved region, with the A-allele ancestral in primates. In a twin/sibling sample (n = 576), A+1485G influenced HR, both at rest (p = 0.010) and after environmental stress (p = 0.002), with the minor (A) allele displaying a recessive effect with lower HR. The effect of A+1485G on HR was extended by meta-analysis into 2 additional population samples (total n = 2,579), and the influence remained directionally consistent and significant (p = 0.007). A+1485G disrupted a microRNA (human microribonucleic acid-1294 [hsa-miR-1294]) recognition motif in the 3'-UTR, as demonstrated by a transfected luciferase reporter/human 3'-UTR variant system in 2 different neuronal/neuroendocrine cell types. The microRNA effect was further documented by cotransfection of an hsa-miR-1294 mimic, yielding an exaggerated decline in expression of the A-allele (better match) reporter (p = 4.3 × 10(-5)). Similar findings of differential 3'-UTR allelic susceptibility to hsa-miR-1294 were noted during expression of the full-length human CYB561 messenger ribonucleic acid with its cognate 3'-UTR. Finally, exposure of stem cell-derived human embryoid bodies to hsa-miR-1294 mimic or antagomir oligonucleotides yielded directionally opposite effects on contraction rate in proto-hearts.

CONCLUSIONS

HR is a substantially heritable trait, with genetic influence by variation in the adrenergic pathway, here shown for messenger ribonucleic acid translational control at the CYB561 step of transmitter formation. The results have implications for potentially modifiable autonomic pathways that influence this risk trait in the population.

Genomics of cardiac electrical function

Proper generation and conduction of the cardiac electrical impulse is essential for the continuous coordinated contraction of the heart. Dysregulation of cardiac electrical function may lead to cardiac arrhythmias, which constitute a huge medical and social burden. Identifying the genetic factors underlying cardiac electrical activity serves the double purpose of allowing the early identification of individuals at risk for arrhythmia and discovering new potential therapeutic targets for prevention. The aim of this review is to provide an overview of the genes and genetic loci linked thus far to cardiac electrical function and arrhythmia. These genes and loci have been primarily uncovered through studies on the familial rhythm disorders and through genome-wide association studies on electrocardiographic parameters in large sets of the general population. An overview of all genes and loci with their respective effect is given.

Large-scale gene function analysis with the PANTHER classification system

The PANTHER (protein annotation through evolutionary relationship) classification system (http://www.pantherdb.org/) is a comprehensive system that combines gene function, ontology, pathways and statistical analysis tools that enable biologists to analyze large-scale, genome-wide data from sequencing, proteomics or gene expression experiments. The system is built with 82 complete genomes organized into gene families and subfamilies, and their evolutionary relationships are captured in phylogenetic trees, multiple sequence alignments and statistical models (hidden Markov models or HMMs). Genes are classified according to their function in several different ways: families and subfamilies are annotated with ontology terms (Gene Ontology (GO) and PANTHER protein class), and sequences are assigned to PANTHER pathways. The PANTHER website includes a suite of tools that enable users to browse and query gene functions, and to analyze large-scale experimental data with a number of statistical tests. It is widely used by bench scientists, bioinformaticians, computer scientists and systems biologists. In the 2013 release of PANTHER (v.8.0), in addition to an update of the data content, we redesigned the website interface to improve both user experience and the system's analytical capability. This protocol provides a detailed description of how to analyze genome-wide experimental data with the PANTHER classification system.