Genetic influences on resting electrocardiographic variables in older women: a twin study

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 combined linkage, microarray and exome analysis suggests MAP3K11 as a candidate gene for left ventricular hypertrophy

BACKGROUND

Electrocardiographic measures of left ventricular hypertrophy (LVH) are used as predictors of cardiovascular risk. We combined linkage and association analyses to discover novel rare genetic variants involved in three such measures and two principal components derived from them.

METHODS

The study was conducted among participants from the Erasmus Rucphen Family Study (ERF), a Dutch family-based sample from the southwestern Netherlands. Variance components linkage analyses were performed using Merlin. Regions of interest (LOD > 1.9) were fine-mapped using microarray and exome sequence data.

RESULTS

We observed one significant LOD score for the second principal component on chromosome 15 (LOD score = 3.01) and 12 suggestive LOD scores. Several loci contained variants identified in GWAS for these traits; however, these did not explain the linkage peaks, nor did other common variants. Exome sequence data identified two associated variants after multiple testing corrections were applied.

CONCLUSIONS

We did not find common SNPs explaining these linkage signals. Exome sequencing uncovered a relatively rare variant in MAPK3K11 on chromosome 11 (MAF = 0.01) that helped account for the suggestive linkage peak observed for the first principal component. Conditional analysis revealed a drop in LOD from 2.01 to 0.88 for MAP3K11, suggesting that this variant may partially explain the linkage signal at this chromosomal location. MAP3K11 is related to the JNK pathway and is a pro-apoptotic kinase that plays an important role in the induction of cardiomyocyte apoptosis in various pathologies, including LVH.

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.

Heritability of ECG Biomarkers in the Netherlands Twin Registry Measured from Holter ECGs

Introduction: The resting ECG is the most commonly used tool to assess cardiac electrophysiology. Previous studies have estimated heritability of ECG parameters based on these snapshots of the cardiac electrical activity. In this study we set out to determine whether analysis of heart rate specific data from Holter ECGs allows more complete assessment of the heritability of ECG parameters.

Methods and Results: Holter ECGs were recorded from 221 twin pairs and analyzed using a multi-parameter beat binning approach. Heart rate dependent estimates of heritability for QRS duration, QT interval, T_peak–T_end and T_height were calculated using structural equation modeling. QRS duration is largely determined by environmental factors whereas repolarization is primarily genetically determined. Heritability estimates of both QT interval and T_height were significantly higher when measured from Holter compared to resting ECGs and the heritability estimate of each was heart rate dependent. Analysis of the genetic contribution to correlation between repolarization parameters demonstrated that covariance of individual ECG parameters at different heart rates overlap but at each specific heart rate there was relatively little overlap in the genetic determinants of the different repolarization parameters.

Conclusions: Here we present the first study of heritability of repolarization parameters measured from Holter ECGs. Our data demonstrate that higher heritability can be estimated from the Holter than the resting ECG and reveals rate dependence in the genetic—environmental determinants of the ECG that has not previously been tractable. Future applications include deeper dissection of the ECG of participants with inherited cardiac electrical disease.

Heritabilities, proportions of heritabilities explained by GWAS findings, and implications of cross-phenotype effects on PR interval

Electrocardiogram (ECG) measurements are a powerful tool for evaluating cardiac function and are widely used for the diagnosis and prediction of a variety of conditions, including myocardial infarction, cardiac arrhythmias, and sudden cardiac death. Recently, genome-wide association studies (GWASs) identified a large number of genes related to ECG parameter variability, specifically for the QT, QRS, and PR intervals. The aims of this study were to establish the heritability of ECG traits, including indices of left ventricular hypertrophy, and to directly assess the proportion of those heritabilities explained by GWAS variants. These analyses were conducted in a large, Dutch family-based cohort study, the Erasmus Rucphen Family study using variance component methods implemented in the SOLAR (Sequential Oligogenic Linkage Analysis Routines) software package. Heritability estimates ranged from 34 % for QRS and Cornell voltage product to 49 % for 12-lead sum. Trait-specific GWAS findings for each trait explained a fraction of their heritability (17 % for QRS, 4 % for QT, 2 % for PR, 3 % for Sokolow–Lyon index, and 4 % for 12-lead sum). The inclusion of all ECG-associated single nucleotide polymorphisms explained an additional 6 % of the heritability of PR. In conclusion, this study shows that, although GWAS explain a portion of ECG trait variability, a large amount of heritability remains to be explained. In addition, larger GWAS for PR are likely to detect loci already identified, particularly those observed for QRS and 12-lead sum.

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.

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.

Heritability of Tpeak-Tend interval and T-wave amplitude: a twin study

BACKGROUND

Tpeak-Tend interval (TpTe) and T-wave amplitude (Tamp) carry diagnostic and prognostic information regarding cardiac morbidity and mortality. Heart rate and QT interval are known to be heritable traits. The heritability of T-wave morphology parameters such as TpTe and Tamp is unknown. TpTe and Tamp were evaluated in a large sample of twins.

METHODS AND RESULTS

Twins from the GEMINAKAR study (611 pairs, 246 monozygotic, 365 dizygotic; mean age, 38±11 years; 49% men) who had an ECG performed during 1997 to 2000 were included. Tamp was measured in leads V1 and V5. Duration variables (RR interval, QTpeak and QTend interval) were measured and averaged over 3 consecutive beats in lead V5. TpTe was calculated as the QTend- and QTpeak-interval difference. Heritability was assessed using structural equation models adjusting for age, sex, and body mass index. All models were reducible to a model of additive genetics and unique environment. All variables had considerable genetic components. Adjusted heritability estimates were as follows: TpTe, 46%; Tamp lead V1, 34%; Tamp lead V5, 47%; RR interval, 55%; QT interval, 67%; and Bazett-corrected QT interval, 42%.

CONCLUSIONS

RR interval, QT interval, Tamp, and TpTe interval are heritable ECG parameters.

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.

Novel genes for QTc interval. How much heritability is explained, and how much is left to find?

The corrected QT (QTc) interval is a complex quantitative trait, believed to be influenced by several genetic and environmental factors. It is a strong prognostic indicator of cardiovascular mortality in patients with and without cardiac disease. More than 700 mutations have been described in 12 genes (LQT1-LQT12) involved in congenital long QT syndrome. However, the heritability (genetic contribution) of QTc interval in the general population cannot be adequately explained by these long QT syndrome genes. In order to further investigate the genetic architecture underlying QTc interval in the general population, genome-wide association studies, in which up to one million single nucleotide polymorphisms are assayed in thousands of individuals, are now being employed and have already led to the discovery of variants in seven novel loci and five loci that are known to cause congenital long or short QT syndrome. Here we show that a combined risk score using 11 of these loci explains about 10% of the heritability of QTc. Additional discovery of both common and rare variants will yield further etiological insight and accelerate clinical applications.

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.

Effects of long-term physical activity on cardiac structure and function: a twin study

Previous studies have shown that athletic training or other physical activity causes structural and functional adaptations in the heart, but less is known how long-term physical activity affects heart when genetic liability and childhood environment are taken into account. The aim of this study was to investigate the effects of long-term physical activity vs. inactivity on cardiac structure and function in twin pairs discordant for physical activity for 32 years. Twelve same-sex twin pairs (five monozygotic and seven dizygotic, 50-67 years) were studied as a part of the TWINACTIVE study. Discordance in physical activity was initially determined in 1975 and it remained significant throughout the follow-up. At the end of the follow-up in 2007, resting echocardiographic and electrocardiographic measurements were performed. During the follow-up period, the active co-twins were on average 8.2 (SD 4.0) MET hours/day more active than their inactive co-twins (p < 0.001). At the end of the follow-up, resting heart rate was lower in the active than inactive co-twins [59 (SD 5) vs. 68 (SD 10) bpm, p=0.03]. The heart rate-corrected QT interval was similar between the co-twins. Also, there was a tendency for left ventricular mass per body weight to be greater and T wave amplitude in lead II to be higher in the active co-twins (18% and 15%, respectively, p=0.08 for both). Similar trends were found for both monozygotic and dizygotic twin pairs. In conclusion, the main adaptation to long- term physical activity is lowered resting heart rate, even after partially or fully controlling for genetic liability and childhood environment. Key pointsThe main adaptation to long-term physical activity is lowering of resting heart rate, even after controlling for genetic liability.VO2peak is increased in the active co-twins compared with their inactive co-twins and accordingly, also submaximal heart rates during the clinical exercise test are lower in physically active co-twins.There is a tendency for increased LVM per body weight and heightened T wave amplitude in the active co-twins.

Electrocardiographic and other clinical correlates of walking ability in older women

The purpose of this study was to examine how resting electrocardiographic (ECG) and other clinical variables, which can be included in a routine clinical examination, predict walking ability in older women. Three hundred and twenty women (63-75 years) without overt cardiac diseases and apparent mobility limitations were studied. Measurements performed were clinical examination (standard 12-lead resting ECG, assessment of physical activity level, presence of chronic diseases, use of beta-blockers, body mass index (BMI), ability to squat, resting blood pressure) and six-minute walking test. Participants walked 533+/-75 m in the six-minute walking test. The best electrocardiographic predictors of long walking distance were high TV(5) and TII, but their explanation rates were small (4.5% and 3.8%, respectively). In hypertensive participants (systolic blood pressure=SBP> or =160 mmHg), the respective values were 9.3% and 5.8%. The best predictors of long walking distance were ability to squat without limitations and low BMI (15.5% and 13.6%, respectively). Altogether the studied variables explained 36% of the variation in walking distance. The data gathered in clinical examination give useful information for the assessment of walking ability in relatively healthy older women. Resting ECG does not give clinically significant additional information for the assessment in subjects without overt cardiac disease.

Electrocardiographic indices of left ventricular hypertrophy and repolarization phase share the same genetic influences: a twin study

BACKGROUND

Both left ventricular hypertrophy (LVH) and repolarization phase (RP) are known to be attributable to genetic influences, but less is known whether they share same genetic influences. The aim of this study was to investigate to what extent individual differences in electrocardiographic (ECG) LVH and RP are explained by genetic and environmental influences and whether these influences are shared between these two traits.

METHODS

Resting ECG recordings were obtained from 186 monozygotic and 203 dizygotic female twin individuals, aged 63 to 76 years. Latent factors, called LVH and RP, were formed to condense the information obtained from LVH indices (Cornell voltage and Cornell product) and T-wave amplitudes (leads V(5) and II), respectively. Multivariate quantitative genetic modeling was used both to decompose the phenotypic variances into additive genetic, common environmental, and unique environmental influences, and for the calculation of genetic and environmental correlations between LVH and RP.

RESULTS

Additive genetic influences explained 16% of individual differences in LVH and 74% in RP. The remaining individual differences were explained by both common and unique environmental influences. The genetic correlation and unique environmental correlation between LVH and RP were -0.93 and -0.05, respectively.

CONCLUSIONS

In older women without overt cardiac diseases, RP is under stronger genetic control than LVH. The majority of genetic influences are shared between LVH and RP whereas environmental influences are mainly specific to each.