Familial aggregation and heritability of electrocardiographic intervals and heart rate in a rural Chinese population

Differential genome-wide associated variants and enriched pathways of ECG parameters among people with versus without HIV

OBJECTIVES

People with HIV (PWH) are more likely to develop ECG abnormalities. Substantial evidence exists for genetic contribution to ECG parameters among general population. However, whether and how would host genome associate with ECG parameters among PWH is unclear. Our research aims to analyze and compare genetic variants, mapped genes, and enriched pathways of ECG parameters among PWH and HIV-negative controls.

DESIGN

A cross-sectional study.

METHOD

We performed a large original genome-wide association study (GWAS) of ECG parameters among PWH ( n  = 1730) and HIV-negative controls ( n =  3746). Genome-wide interaction analyses were also conducted.

RESULTS

A total of 18 novel variants were detected among PWH, six for PR interval including rs76345397 at ATL2 , 11 for QRS duration including rs10483994 at KCNK10 and rs2478830 at JCAD , and one for QTc interval (rs9815364). Among HIV-negative controls, we identified variants located at previously reported ECG-related genes ( SCN5A , CNOT1 ). Genetic variants had a significant interaction with HIV infection ( P  < 5 × 10 -8 ), implying that HIV infection and host genome might jointly influence ECG parameters. Mapped genes for PR interval and QRS duration among PWH were enriched in the biological process of viral genome replication and host response to virus, respectively, whereas enriched pathways for PR interval among HIV-negative controls were in the cellular component of voltage-gated sodium channel complex.

CONCLUSION

The present GWAS indicated a distinctive impact of host genome on quantitative ECG parameters among PWH. Different from HIV-negative controls, host genome might influence the cardiac electrical activity by interfering with HIV viral infection, production, and latency among PWH.

Heritability of a resting heart rate in a 20-year follow-up family cohort with GWAS data: Insights from the STANISLAS cohort

BACKGROUND

The association between resting heart rate (HR) and cardiovascular outcomes, especially heart failure, is now well established. However, whether HR is mainly an integrated marker of risk associated with other features, or rather a genetic origin risk marker, is still a matter for debate. Previous studies reported a heritability ranging from 14% to 65%.

DESIGN

We assessed HR heritability in the STANISLAS family-study, based on the data of four visits performed over a 20-year period, and adjusted for most known confounding effects.

METHODS

These analyses were conducted using a linear mixed model, adjusted on age, sex, tea or coffee consumption, beta-blocker use, physical activity, tobacco use, and alcohol consumption to estimate the variance captured by additive genetic effects, via average information restricted maximum likelihood analysis, with both self-reported pedigree and genetic relatedness matrix (GRM) calculated from genome-wide association study data.

RESULTS

Based on the data of all visits, the HR heritability (h2) estimate was 23.2% with GRM and 24.5% with pedigree. However, we found a large heterogeneity of HR heritability estimations when restricting the analysis to each of the four visits (h2 from 19% to 39% using pedigree, and from 14% to 32% using GRM). Moreover, only a little part of variance was explained by the common household effect (<5%), and half of the variance remained unexplained.

CONCLUSION

Using a comprehensive analysis based on a family cohort, including the data of multiple visits and GRM, we found that HR variability is about 25% from genetic origin, 25% from repeated measures and 50% remains unexplained.

A genome-wide association and polygenic risk score study on abnormal electrocardiogram in a Chinese population

Electrocardiography is a common and widely-performed medical examination based on the measurement and evaluation of electrocardiogram (ECG) to assess the up-to-date cardiac rhythms and thus suggest the health conditions of cardiovascular system and on a larger level the individual’s wellness. Abnormal ECG assessment from the detection of abnormal heart rhythms may have clinical implications including blood clots in formation, ongoing heart attack, coronary artery blockage, etc. Past genetic-phenotypic research focused primarily on the physical parameters of ECG but not the medical evaluation. To unbiasedly uncover the underlying links of genetic variants with normal vs. abnormal ECG assessment, a genome-wide association study (GWAS) is carried out in a 1006-participant cohort of Chinese population effectively genotyped for 243487 single nucleotide polymorphisms (SNPs). Both age and sex are influential factors, and six novel SNPs are identified for potential association with abnormal ECG. With the selected SNPs, a polygenic risk score (PRS) differentiates the case–control subgroups, and correlates well with increased risk of abnormal ECG. The findings are reproduced in an independent validation cohort. The derived PRS may function as a potential biomarker for prospectively screening the high-risk subgroup of heart issues in the Chinese population.

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.

Familial Clustering of Cardiac Conduction Defects and Pacemaker Insertion

BACKGROUND

The etiopathogenesis of electrocardiographic bundle branch and atrioventricular blocks is not fully understood. We investigated familial clustering of cardiac conduction defects and pacemaker insertion in the FHS (Framingham Heart Study). Additionally, we assessed familial clustering of pacemaker insertion in the Danish general population.

METHODS

In FHS, we used multivariable-adjusted logistic regression models to investigate the association of parental atrioventricular block (PR interval, ≥0.2 s), complete bundle branch block (QRS, ≥0.12 s), or pacemaker insertion with the occurrence of cardiac conduction abnormalities in their offspring. The Danish nationwide administrative registries were interrogated to assess the relations of parental pacemaker insertion with offspring pacemaker insertion.

RESULTS

In FHS (n=371 cases with first-degree atrioventricular block, complete bundle branch block, or pacemaker insertion, and 1471 age- and sex-matched controls), individuals with at least 1 affected parent with a conduction defect had a 1.65-fold odds (odds ratio, 95% CI, 1.32-2.07) for manifesting an atrioventricular block and a 1.62-fold odds (95% CI, 1.08-2.42) for developing a complete bundle branch block. If at least 1 parent had any electrocardiographic conduction defect or pacemaker insertion, the offspring had a 1.62-fold odds (95% CI, 1.31-2.00) for experiencing any of these conditions. In Denmark (n=2 824 199 individuals; 5397 incident pacemaker implantations), individuals with at least 1 first-degree relative with history of pacemaker insertion had a multivariable-adjusted 1.68-fold (incidence rate ratio, 95% CI, 1.49-1.89) risk of undergoing a pacemaker insertion. If the affected relative was ≤45 years of age, the incidence rate ratio was markedly increased to 51.0 (95% CI, 32.7-79.9).

CONCLUSIONS

Cardiac conduction blocks and risk for pacemaker insertion cluster within families. A family history of conduction system disturbance or pacemaker insertion should trigger increased awareness of a similar propensity in other family members, especially so when the conduction system disease occurs at a younger age.

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.

A Comparison of Heritability Estimates by Classical Twin Modeling and Based on Genome-Wide Genetic Relatedness for Cardiac Conduction Traits

Twin studies have found that ~50% of variance in electrocardiogram (ECG) traits can be explained by genetic factors. However, genetic variants identified through genome-wide association studies explain less than 10% of the total trait variability. Some have argued that the equal environment assumption for the classical twin model might be invalid, resulting in inflated narrow-sense heritability (h 2) estimates, thus explaining part of the 'missing h 2'. Genomic relatedness restricted maximum likelihood (GREML) estimation overcomes this issue. This method uses both family data and genome-wide coverage of common SNPs to determine the degree of relatedness between individuals to estimate both h 2 explained by common SNPs and total h 2. The aim of the current study is to characterize more reliably than previously possible ECG trait h 2 using GREML estimation, and to compare these outcomes to those of the classical twin model. We analyzed ECG traits (heart rate, PR interval, QRS duration, RV5+SV1, QTc interval, Sokolow-Lyon product, and Cornell product) in up to 3,133 twins from the TwinsUK cohort and derived h 2 estimates by both methods. GREML yielded h 2 estimates between 47% and 68%. Classical twin modeling provided similar h 2 estimates, except for the Cornell product, for which the best fit included no genetic factors. We found no evidence that the classical twin model leads to inflated h 2 estimates. Therefore, our study confirms the validity of the equal environment assumption for monozygotic and dizygotic twins and supports the robust basis for future studies exploring genetic variants responsible for the variance of ECG traits.

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.

Identification of three novel genetic variations associated with electrocardiographic traits (QRS duration and PR interval) in East Asians

The electrocardiogram has several advantages in detecting cardiac arrhythmia-it is readily available, noninvasive and cost-efficient. Recent genome-wide association studies have identified single-nucleotide polymorphisms that are associated with electrocardiogram measures. We performed a genome-wide association study using Korea Association Resource data for the discovery phase (Phase 1, n = 6805) and two consecutive replication studies in Japanese populations (Phase 2, n = 2285; Phase 3, n = 5010) for QRS duration and PR interval. Three novel loci were identified: rs2483280 (PRDM16 locus) and rs335206 (PRDM6 locus) were associated with QRS duration, and rs17026156 (SLC8A1 locus) correlated with PR interval. PRDM16 was recently identified as a causative gene of left ventricular non-compaction and dilated cardiomyopathy in 1p36 deletion syndrome, which is characterized by heart failure, arrhythmia and sudden cardiac death. Thus, our finding that a PRDM16 SNP is linked to QRS duration strongly implicates PRDM16 in cardiac function. In addition, C allele of rs17026156 increases PR interval (beta ± SE, 2.39 ± 0.40 ms) and exists far more frequently in East Asians (0.46) than in Europeans and Africans (0.05 and 0.08, respectively).

Common SCN10A variants modulate PR interval and heart rate response during atrial fibrillation

AIMS

SCN10A encodes the sodium channel Nav1.8 implicated by genome-wide association studies as a modulator of atrioventricular conduction (PR interval). In a cohort of patients with atrial fibrillation (AF), we examined whether there was an association between common variants in SCN10A and both the PR interval during normal sinus rhythm and the heart rate response during AF.

METHODS AND RESULTS

Patients prospectively enrolled in the Vanderbilt AF registry with electrocardiograms in normal sinus rhythm and/or AF within 1 year of enrollment were genotyped for two common SCN10A variants rs6795970 and rs12632942. Both variants were associated with the PR interval duration in a gene-dose effect on unadjusted analysis; after adjustment for the covariates age, gender, body mass index, hypertension, congestive heart failure, and medication usage, the association remained for rs6795970 only (P = 0.012, partial R(2) = 0.0139). On unadjusted analysis, heart rate response during AF was associated with rs6795970 (P = 0.035, partial R(2) = 0.015), but not with rs12632942 (P = 0.89), and neither association was significant after adjustment for covariates.

CONCLUSION

The common variant rs6795970 in SCN10A is associated with the PR interval duration among healthy patients and those with AF. In addition, this single nucleotide polymorphism trended towards an association with heart rate response during AF indicating the importance of this common SCN10A polymorphism as a marker of atrioventricular conduction.

Common genetic variation modulating cardiac ECG parameters and susceptibility to sudden cardiac death

Sudden cardiac death (SCD) is a prevalent cause of death in Western societies. Genome-wide association studies (GWAS) conducted over the last few years have uncovered common genetic variants modulating risk of SCD. Furthermore, GWAS studies uncovered several loci impacting on heart rate and ECG indices of conduction and repolarization, as measures of cardiac electrophysiological function and likely intermediate phenotypes of SCD risk. We here review these recent developments and their implications for the identification of novel molecular pathways underlying normal electrophysiological function and susceptibility to SCD.

The emerging genetic landscape underlying cardiac conduction system function

Proper function of an organized Cardiac Conduction System (CCS) is vital to the survival of metazoans ranging from fly to man. The routine use of non-invasive electrocardiogram measures in the diagnosis and monitoring of cardiovascular health has established a trove of reliable CCS functional data in both normal and diseased cardiac states. Recent combination of echocardiogram (ECG) data with genome-wide association studies has identified genomic regions implicated in ECG variability which impact CCS function. In this study, we review the substantial recent progress in this area, highlighting the identification of novel loci, confirming the importance of previously implicated loci in CCS function, and exploring potential links between genes with important roles in developmental processes and variation in function of the CCS.

Several common variants modulate heart rate, PR interval and QRS duration

Electrocardiographic measures are indicative of the function of the cardiac conduction system. To search for sequence variants that modulate heart rate, PR interval and QRS duration in individuals of European descent, we performed a genome-wide association study in approximately 10,000 individuals and followed up the top signals in an additional approximately 10,000 individuals. We identified several genome-wide significant associations (with P < 1.6 x 10(-7)). We identified one locus for heart rate (MYH6), four for PR interval (TBX5, SCN10A, CAV1 and ARHGAP24) and four for QRS duration (TBX5, SCN10A, 6p21 and 10q21). We tested for association between these loci and subjects with selected arrhythmias in Icelandic and Norwegian case-control sample sets. We observed correlations between TBX5 and CAV1 and atrial fibrillation (P = 4.0 x 10(-5) and P = 0.00032, respectively), between TBX5 and advanced atrioventricular block (P = 0.0067), and between SCN10A and pacemaker implantation (P = 0.0029). We also replicated previously described associations with the QT interval.