Article Text
Abstract
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, χ2=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.
- GENETICS
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
Key MESSAGES
What is already known about this subject?
Bundle branch block is associated with an increased risk of sudden cardiac death. Reduced levels of connexin 40 are associated with bundle branch block and reduced levels of connexin 43 are associated with increased risk of ventricular arrhythmias.
What does this study add?
The present study indicates that bundle branch block might be due to a diffuse myocardial process within the myocardium and not only due to ischaemia of the specialised conduction system. It also indicated that there might be a genetic susceptibility for bundle branch block development.
How might this impact on clinical practice?
This might lead to an increased awareness that bundle branch block could be due to a diffuse myocardial process and thereby stimulate further research in this field. This further research might identify responders to cardiac resynchronisation therapy (CRT) more easily.
Introduction
The clinical implications of bundle branch block (BBB) have long been a subject of dispute.1 The ECG alone does not allow us to determine whether the conduction disturbance is a function of a focal lesion in main bundle branches or the result of a more diffuse myocardial disease producing a type of parietal block.1 ,2
In a previous study of a Swedish general population sample of men followed from 50 to 80 years of age, The Study of Men Born 1913, we found 86% of the BBB to be acquired and not correlated to ischaemic heart disease or its risk factors.3 Similar findings have been observed in a large Finnish data set of more than 10 000 participants,4 whereas in the US Framingham data set, the majority of BBB found was associated with cardiovascular disease or its risk factors.5 ,6 These results support the theory that BBB is a marker of a progressive disease that also affects conduction within the myocardium.3 ,4
The development of cardiac resynchronisation therapy for patients with BBB and heart failure has stimulated research on BBB,7 ,8 but the underlying nature of this progressive disorder is still mostly unknown. It is also unsolved whether BBB should be regarded as a single phenotype or if right BBB (RBBB) and left BBB (LBBB) are different entities. Therefore, we analysed all three, BBB, RBBB and LBBB, in the present study.
Molecular studies on animal and on human data have shown that variation in connexin 40 and connexin 43 can cause BBB and is associated with cardiomyopathy.9–11 Connexin 40 is expressed in the atria and in the proximal parts of the conduction system, coexpressed with connexin 43 in Purkinje fibres, and mice lacking connexin 40 have prolonged PQ interval, AV interval and QRS duration due to reduced conduction within the Purkinje fibres.12 It has also been shown to produce a phenotype of RBBB.13 In contrast, connexin 43 is expressed in the Purkinje fibres and working cardiomyocytes. In heterozygous knock-out mice with reduced connexin 43 levels, ventricular conduction was markedly reduced and had increased risk of ventricular arrhythmias.14 Moreover, human studies on heart failure have shown that connexin 43 expression is reduced in end-stage heart failure,10 and this altered expression is associated with an increased risk of ventricular arrhythmias.11
The hypothesis for the present study is that there is a hereditary predisposition to acquire BBB and that predisposition is associated with polymorphisms in genes coding for connexin 40 and connexin 43. This hypothesis is tested in three successive steps:
Heredity of the ECG phenotype BBB or prolonged QRS duration is tested in the descendants of men with BBB from the population sample of The Study of Men Born 1913.
Genetic variation at the connexin 40 and connexin 43 loci is explored by genotyping samples from The Study of Men Born 1913 with and without BBB.
The association of identified polymorphisms with BBB is subsequently tested in the independent INTERGENE sample.
Methods
Three samples for testing the three different steps of our hypothesis were used. The relation between populations for the analysis of heredity and the resequencing analysis is depicted in figure 1. The association analysis was performed in an independent sample as described below.
Patients for step 1: heredity—The Study of Men Born 1913 and their descendants
The Study of Men Born 1913 is a randomly selected general population of men followed since 1963.15 For this paper, we use data from 1963 until 1993, when 361 men (then aged 80 years) were still alive; 212 (59%) participated in the follow-up examination. Among the participating men, a diagnosis of BBB was made in 82 patients during 30 years of follow-up: 60 with RBBB and 22 with LBBB, as previously reported.3 The majority (86%) had acquired BBB unrelated to ischaemic heart disease. Of these 82 men, we were able to construct pedigrees in 77 of the men's families; the families of the remaining 5 men could not be traced. The 77 men had 116 sons and daughters, all of whom were invited to an investigation including physical examination, standardised history forms, spirometry and 12-lead ECG recording. Of the 116 descendants, 76 (66%) were willing to participate and they formed the basis of the report on the ECG phenotype (table 1).
Patients for step 2: resequencing
In order to separate the ECG phenotype optimally for the genetic analysis, we chose two groups from The Study of Men Born 1913.3 The first group consisted of 10 men with BBB (8 had RBBB and 2 had LBBB) who had the longest QRS duration from the youngest age and the second group was a contrast group of 10 men at age 80 without BBB and with the shortest possible QRS duration. The 10 men with BBB had a mean QRS duration of 136 ms that was present at the mean age of 63.5 years. The 10 men in the contrast group had a mean QRS duration of 78 ms and a mean age of 80 years. In addition to these 20 participants, 2 descendants identified as having BBB, 1 with LBBB and 1 with RBBB, were included in the screening panel, resulting in 10 participants without BBB and 12 with BBB.
Patients for step 3: association
In total, 3600 participants aged 25–70 years were randomly sampled from individuals living in Västra Götalands region as controls into the INTERGENE study cohort. Clinical examinations, ECG and collection of blood samples for DNA analysis were carried out as described.16 Out of the INTERGENE study cohort, 121 participants had an ECG indicating a possible BBB (QRS over 120 ms) and an adequate blood sample for DNA analysis. ECG classification by two cardiologists (see ECG below) resulted in 108 participants with BBB for the present study, 35 with LBBB and 73 with RBBB. For each case with possible BBB, 2 controls from the population cohort without clinically overt cardiovascular disease and without BBB, matched for age and sex, were selected (table 2).
Genetic analysis—resequencing
DNA was extracted from whole blood (if available) or serum. The screening panel was studied by resequencing by di-deoxy sequencing, as described elsewhere.17 A screening panel of 44 chromosomes resulted in 99% power to find any biallelic polymorphism with an allelic frequency of 10%, and 90% power for polymorphisms with an allelic frequency of 5%. For the connexin 40 gene, all described exons were analysed, as well as 1200 bp of the promoter; for the connexin 43 gene, all exons as well as 1000 bp of the promoter were sequenced.
Genetic analysis—genotyping
Genotyping was performed by TaqMan assays on the 7900 HT system (Applied Biosystems, Life Technologies, New York, USA). Primer and probe data are available from the authors on request.
ECG
Standard 12-lead ECGs were recorded with participants at rest in the supine position. Paper speed was 50 mm/s, and calibration was 1 mV:10 mm. ECGs with a QRS duration above 120 ms were considered for the study. Two independent cardiologists (PE and MD), classified the ECGs as RBBB, LBBB or not BBB (including hypertrophy and Wolff-Parkinson-White syndrome, which were excluded from further analyses). LBBB was defined as QRS duration above 120 ms, PQ interval above 120 ms, upright and slurred R-waves in leads I, V5, V6 and QS, or rS pattern in lead V1.18 The only deviation from the interpretation proposed by AHA19 is that we also included LBBB with a small q in leads V5 and V6. RBBB20 was defined as QRS duration above 120 ms, PQ interval above 120 ms, rSR′ in lead V1 or V2, and S-waves in lead I and either lead V5 or V6. RBBB was further classified as having concomitant left anterior hemiblock (LAH) if the QRS axis was below −30° and concomitant left posterior hemiblock (LPH) if the QRS axis was above 90°.20 Discrepancies were resolved by consensus. Bifascicular block was defined as LBBB or RBBB with concomitant LAH or LPH.
Statistical analyses and bioinformatics
IBM SPSS V.20.0 (IBM Corp, Armonk, New York, USA) was used for statistical analyses. Continuous variables are expressed as mean±SE of the mean. Proportions are expressed as percentages. Allele frequencies were estimated from genotype data. The χ2 test was used to assess association per genotype and to evaluate deviations from the Hardy-Weinberg equilibrium. Testing for linkage disequilibrium (LD) was performed by the maximum likelihood method using the linkage utility program EH for estimation of haplotype frequencies.21 Haplotype effects on BBB and QRS were studied with the THESIAS program 3.1.22 Computer search for activation of potential transcription factor binding sites by the single nucleotide polymorphisms (SNPs) was performed with SNPInspector, searching the MatInspector database.23 Computer search for potential influence on messenger RNA (mRNA) stability was performed with UTRdb.24 LD between the identified polymorphism at the connexin 43 locus and the SNP rs1015451 identified by previous genome-wide association studies was performed by PLINK.25
Missing data
In the association analysis, one patient had missing data for the connexin 40 679 genotype and one patient for the connexin 43 1400 A/insertion genotype.
Ethics statement
The protocol was approved by the ethical board of the University of Gothenburg, and patients provided written informed consent.
Results
Step 1: heredity
Two sons with complete BBB were identified (6.4% of all participating men). The only complete RBBB was in a 56-year-old son with a RBBB with QRS duration of 132 ms. His father had had an RBBB with QRS duration of 135 ms when measured at age 67. No BBB or disturbances in ECG measurements were found in the daughters (table 1). None of the two sons with complete BBB had any symptoms or ECG findings indicative of ischaemic heart disease or heart failure.
Step 2: resequencing
By resequencing connexin 40, six SNPs were identified at positions −992, −679, −55, +197, +312 and +2469 relative to transcription start and the resulting mRNA in GenBank (BC013313). All SNPs were in Hardy-Weinberg equilibrium. Frequencies and corresponding rs numbers are shown in table 3 and linkage equilibrium in table 4. Searching for putative transcription factor binding sites using SNPInspector revealed that although the connexin 40 SNPs altered transcription factor binding sites for sites located in the connexin 40 promoter (−992, −679, −55) and in the first intron (+197, +312), none of the listed transcription factors has been shown to be present in the heart (SNPInspector website search, 4 June 2014). The connexin 40 2469G>A SNP was located in the 3′ flanking region, 56 bp downstream of last amino acid, and did not affect mRNA stability according to UTRdb.
Resequencing of connexin 43 resulted in the identification of only one polymorphism, which was a single nucleotide insertion/deletion polymorphism in the non-coding part of the second exon at position 1400 relative to transcription start and resulting mRNA in GenBank (NM_000165.3). The connexin 43 1400A/ins was in Hardy-Weinberg equilibrium, had a minor allele frequency of 0.04 and is annotated rs67678923. The inserted A is inserted immediately after the coding sequence in the first nucleotide position in the 3′-untranslated region (UTR).
Step 3: Association
LBBB was more evenly distributed among men and women, whereas RBBB was more common in men (table 2). Hypertension, angina pectoris and hypercholesterolaemia were equally distributed among LBBB and RBBB and, as expected, they were more prevalent in these subjects than in the control group without cardiovascular disease. Subjects with LBBB more often had a history of myocardial infarction, diabetes and irregular pulse compared to the control group but also compared to subjects with RBBB (table 2).
Genotype distribution of the connexin 43 1400A/ins polymorphism was associated with LBBB (χ2=7.3, p=0.03) but not with RBBB, and thus not with total BBB (table 5). Owing to the frequency of the connexin 43 variant, a recessive model cannot be tested in the present material (only one subject homozygous for the variant). Testing with allelic OR is also hampered by low numbers. In a dominant model OR is 1.87 (0.59–5.97), β=0.626, p=0.29. If adjusted for age and sex the OR is 1.95 (0.60–6.32), β=0.67, p=0.26.When analysing the connexin 43 with QRS duration as a continuous variable, no significant difference was observed, although the single homozygous individual for connexin 43 1400A/ins polymorphism had a very long QRS duration.
None of the individually identified SNPs in connexin 40 were associated with BBB, either LBBB, RBBB or total BBB (table 5).
Based on connexin 40 genotype data and haplotypes of the screening panel, these haplotypes were tagged by four SNPs in the INTERGENE cohort. Four haplotypes occurred with a frequency of over 5% and these haplotypes corresponded to 95% of alleles; there were seven haplotypes with a frequency over 1% (table 6).
To further explore any relation of conduction and genetic variation at the connexin 40 locus, we performed a regression analysis of QRS duration as a continuous variable and haplotypes of connexin 40. Haplotype analysis of connexin 40 showed no relation to QRS duration. Individually assigned haplotypes of connexin 40 and QRS duration are displayed in figure 2.
Discussion
The association of LBBB and connexin 43 1400A/ins polymorphism together with this and previous studies on heritability, indicate that BBB development can be modulated by genetic factors and could be the result of a progressive disorder within the myocardium and not confined to the specialised conduction system.
In our study on heritability, we report the first study focused on BBB in fathers and their descendants in their 50s. From our previous study of the fathers, The Study of Men Born 1913,3 we know that at age 50 the prevalence of BBB among men is 1.2%. The sons had a longer QRS duration compared with the daughters. This, and the finding that two (6.4%) of the sons also had BBB, but none of the daughters did, might indicate a male predominance of BBB. Previous studies have shown a male predominance of RBBB but not for LBBB.26 The male predominance of BBB may be a paraphenomenon with cardiovascular disease, in keeping with the fact that women generally develop cardiovascular disease later than men; alternatively, it may represent a sex-linked inheritance of a conduction disorder unrelated to cardiovascular disease. Unfortunately, the present study cannot give us the answer to this question. Although not focused on BBB, previous twin and family studies have reported a heritability of QRS duration of 0.36–0.43.27 ,28
We chose a candidate gene approach to explore the association between heritable factors and QRS duration. On the basis of knock-out mice and functional studies,9 ,11 we chose connexin 40 and connexin 43. With our extreme phenotypes approach, we identified 6 SNPs in connexin 40 and 1 insertion/deletion polymorphism in connexin 43. This is lower than expected, despite the fact that the present study had reasonable power (over 99% for biallelic variants with a frequency above 10%). Our screening approach with extreme phenotypes might influence the number of identified variants, as we only intended to identify plausible functional SNPs associated with QRS width. It is likely that a random sample might identify more genetic variants.
The finding of connexin 43 and the association to LBBB is supported by genome-wide association studies (GWAS) as well as by functional data. In a large GWAS, genetic markers close to the gene for connexin 43, GJA1, showed a strong association with QRS duration.29 The allele frequency for our identified connexin 43 polymorphism, rs67678923, and the strongest marker of the GJA1 locus in the GWAS, rs1015451, differs, and therefore the linkage expressed as R2 and D′ differs; using PLINK in 1000 Genomes project (1092 individuals in phase 1), R2 was 0.02 and D′ was 0.89.25
From a functional point of view, there is a growing body of evidence that connexin 43 is a major player in cardiac conduction.14 ,30 Connexin 43 can be regulated at the post-transcriptional level, thereby affecting connexin 43 turnover and function.30 Ischaemia induces dephosphorylation of connexin 43 and its translocation from the gap junction to intracellular stores, with reduced conduction as a result.31 The connexin 43 insertion/deletion polymorphism (rs144617620) has an inserted A immediately after the coding sequence in the first nucleotide position in the 3′-UTR. The association of connexin 43 polymorphism and LBBB might be explained by post-translational regulation of connexin 43 with reduced mRNA stability and connexin 43 levels, and thereby conduction. Functional characterisation of the identified connexin 43 polymorphism would of course strengthen our finding and minimise the risk that it only represents a chance finding or is in LD with a functional variant. However, such an analysis is beyond the scope of the present paper, and should be carried out in separate studies.
Genetic variation at the connexin 40 locus has previously been shown to be associated with atrial fibrillation, but not with BBB or QRS duration.32 As these SNPs were not associated with our phenotype of interest, their further characterisation is beyond the scope of this paper.
Whether to study BBB as a group or RBBB and LBBB as separate entities is a matter of dispute. Most studies looking at LBBB and RBBB as separate entities are based on hospitalised patients; the majority have ischaemic heart disease where RBBB is more frequent than LBBB, due to myocardial scarring after myocardial infarction.33 ,34 Data from our epidemiological study3 suggest that BBB develops gradually. We hypothesised that this could be due to a progressive conduction disorder and decided to use a genetic approach without any assumptions about whether this could affect LBBB, RBBB or all types of BBB.
The present paper has some strengths and weaknesses. As shown by our previous study, BBB is a late-onset phenotype.3 The access to 80-year-old men without BBB represents true controls for people without BBB, which is hard to get. For the same reason, it is possible that the adult descendants in the heritability study were too young to develop BBB. Reassessment at an older age would probably increase the probability of finding more BBBs among them. In this study we only had ECG information for the fathers, but data for both parents would improve heritability assessment, especially on X linked traits. Regarding the association study, the INTERGENE study was a case–control study on ischaemic heart disease. Controls in the present study were derived from the population sample without prior cardiovascular disease and were thus less likely to have cardiovascular risk factors. The high prevalence of ischaemic heart disease and its risk factors in the present study does not necessarily reflect the true prevalence in people with or without BBB in the general population. Moreover, although a specific hypothesis was tested in this small study, the finding of connexin 43 polymorphism associated with LBBB must be replicated in larger studies with better power.
In conclusion, BBB is the result of a prolonged conduction within the heart. This can be caused by acute ischaemic lesions to the conduction system, that is, myocardial infarction. The heritability of longer QRS duration, and the association of LBBB and connexin 43 polymorphism, indicate that it can also be caused by a progressive conduction disorder within the myocardium distal to the specialised conduction system.
Acknowledgments
The authors thank research nurse Helena Dellborg for assistance in the descendants study. We thank Cecilia Lundholm for technical assistance in the sequencing and genotyping study. They also thank the Genomics Core Facility at the Sahlgrenska Academy, University of Gothenburg, for providing genotyping facilities.
References
Footnotes
Contributors PL drafted the article, designed the study including genetic analysis and performed statistical analyses. BA designed and performed genetic analyses and bioinformatics. P-OH and HE were responsible for the 1913 study. DT was responsible for the INTERGENE study. MD interpreted ECGs. PE designed the study, was responsible for the descendants study and interpreted ECGs. All authors were involved in the interpretation of results and reviewing the manuscript before submission.
Funding The study was funded by grants from the Swedish state under the LUA/ALF agreement.
Competing interests None.
Ethics approval Ethical board of Gothenburg University.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement Primer and probe data are available to readers on request by email to the corresponding author.