Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Evaluation of the Genetic Basis of Familial Aggregation of Pacemaker Implantation by a Large Next Generation Sequencing Panel

  • Patrícia B. S. Celestino-Soper ,

    Contributed equally to this work with: Patrícia B. S. Celestino-Soper, Anisiia Doytchinova

    Affiliation Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • Anisiia Doytchinova ,

    Contributed equally to this work with: Patrícia B. S. Celestino-Soper, Anisiia Doytchinova

    Affiliation Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • Hillel A. Steiner,

    Affiliations Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America, Baruch Padeh Medical Center, Poriya MP Lower Galilee, Israel, Faculty of Medicine in the Galilee, Bar-Ilan University, Safed, Israel

  • Andrea Uradu,

    Affiliation Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • Ty C. Lynnes,

    Affiliation Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • William J. Groh,

    Affiliation Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • John M. Miller,

    Affiliation Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • Hai Lin,

    Affiliations Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, United States of America, Center for Computational Biology and Bioinformatics, Indiana University Purdue University Indianapolis, Indianapolis, IN, United States of America

  • Hongyu Gao,

    Affiliation Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • Zhiping Wang,

    Affiliation Center for Computational Biology and Bioinformatics, Indiana University Purdue University Indianapolis, Indianapolis, IN, United States of America

  • Yunlong Liu,

    Affiliations Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, United States of America, Center for Computational Biology and Bioinformatics, Indiana University Purdue University Indianapolis, Indianapolis, IN, United States of America

  • Peng-Sheng Chen,

    Affiliation Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America

  • Matteo Vatta

    mvatta@iu.edu

    Affiliations Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, United States of America, Krannert Institute of Cardiology, Division of Cardiology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, United States of America

Correction

15 Jan 2016: Celestino-Soper PBS, Doytchinova A, Steiner HA, Uradu A, Lynnes TC, et al. (2016) Correction: Evaluation of the Genetic Basis of Familial Aggregation of Pacemaker Implantation by a Large Next Generation Sequencing Panel. PLOS ONE 11(1): e0147455. https://doi.org/10.1371/journal.pone.0147455 View correction

Abstract

Background

The etiology of conduction disturbances necessitating permanent pacemaker (PPM) implantation is often unknown, although familial aggregation of PPM (faPPM) suggests a possible genetic basis. We developed a pan-cardiovascular next generation sequencing (NGS) panel to genetically characterize a selected cohort of faPPM.

Materials and Methods

We designed and validated a custom NGS panel targeting the coding and splicing regions of 246 genes with involvement in cardiac pathogenicity. We enrolled 112 PPM patients and selected nine (8%) with faPPM to be analyzed by NGS.

Results

Our NGS panel covers 95% of the intended target with an average of 229x read depth at a minimum of 15-fold depth, reaching a SNP true positive rate of 98%. The faPPM patients presented with isolated cardiac conduction disease (ICCD) or sick sinus syndrome (SSS) without overt structural heart disease or identifiable secondary etiology. Three patients (33.3%) had heterozygous deleterious variants previously reported in autosomal dominant cardiac diseases including CCD: LDB3 (p.D117N) and TRPM4 (p.G844D) variants in patient 4; TRPM4 (p.G844D) and ABCC9 (p.V734I) variants in patient 6; and SCN5A (p.T220I) and APOB (p.R3527Q) variants in patient 7.

Conclusion

FaPPM occurred in 8% of our PPM clinic population. The employment of massive parallel sequencing for a large selected panel of cardiovascular genes identified a high percentage (33.3%) of the faPPM patients with deleterious variants previously reported in autosomal dominant cardiac diseases, suggesting that genetic variants may play a role in faPPM.

Introduction

Cardiac conduction disease (CCD) and sick sinus syndrome (SSS) are common indications for pacemaker implantation [14]. According to the current guidelines, genetic testing may be considered for patients with isolated CCD (ICCD) or CCD with concomitant congenital heart disease, especially with documented remarkable family history of cardiac conduction disease. Genetic testing for ICCD remains a Class IIb recommendation, as variants in SCN5A, which cause the majority of familial cardiac conduction disease cases, account for only about 5% of the disease burden [5]. While other genes including SCN1B, TRPM4, and HCN4 have been associated with idiopathic sinus node dysfunctions (SNDs) [5], it is unclear if more comprehensive testing would yield additional useful results.

The development of next generation sequencing (NGS) has brought an economical and efficient way to detect variants in thousands of target regions. Several targeting approaches have been developed for NGS including whole genome sequencing (WGS), whole exome sequencing (WES), and disease-targeted panels. WGS and WES methods are better suited for very rare disorders for which no other clinical genetic tests are available. In comparison, the gene panels can be designed with a specific clinical phenotype/s or disease/s in mind, allowing a more targeted approach for result interpretation and subsequent patient management recommendations. Unlike WGS and WES, gene panels also permit a better sequencing coverage of targeted regions, more accurate and reliable analyses, and overall lower cost [6, 7]. Here we describe the application of a custom NGS targeted enrichment method to implement a research panel for the study of selected patients with familial aggregation of permanent pacemaker implantation (faPPM). Currently, there is limited knowledge about the role of genetic predisposition to faPPM. Here we present the findings resulting from a comprehensive genetic approach using a large panel including genes associated with cardiomyopathies, arrhythmias, congenital heart defects, aortopathy, connective tissue disorders, Noonan spectrum disorders, pulmonary arterial hypertension, metabolic disorders with cardiac presentation, and lipid disorders.

Materials, Methods, and Patient Selection Methodology

This study protocol was approved by the Institutional Review Board of the Indiana University School of Medicine. All participants have provided written informed consent to participate in this study. A total of 112 patients with permanent pacemaker (PPM) implantation were recruited from a pacemaker clinic. Among them, nine (8%) presented with ICCD or SSS without structural heart disease and had at least one first degree relative (FDR) with PPM, thus described to have faPPM and their DNAs were used for NGS analyses. Table 1 summarizes the characteristics of all recruited patients including those submitted for genetic analysis. Exclusion criteria for patients not submitted for genetic analysis are listed in Table A in S1 File. A description of each patient at the time of their study enrollment and a summary of their clinical course as available in the medical records is described in S1 File (patient description section).

Post-bioinformatics filtering of variants found in the 9 pacemaker individuals started with categorization of variants as either present or absent in the Human Gene Mutation Database (HGMD) [8]. Table H in S1 File summarizes the initial classification of HGMD variants, which was followed by re-classification, if needed, after careful literature analyses. The following literature evidence was used to re-classify the HGMD variants:

  1. Affects protein function
    1. Affects function in isolation: truncating variant in gene where loss of function is a known mechanism of disease; or variant found to be disease causing from linkage studies or co-segregation studies in large families. In all cases, there must be sufficient functional data to support a deleterious effect of the variant and the variant must have frequency of less than or equal to 5% in all combined populations in both the Exome Sequencing Project (ESP) and the 1000 Genomes project (1000G) (hg19 data from ANNOVAR [9]). These variants were confirmed by Sanger sequencing.
    2. May or may not affect function in isolation: de novo variant in a patient with the disease and unaffected parents; or co-segregation studies in small to mid-sized families. In all cases, there must be sufficient functional data to support a deleterious effect of the variant and the variant must have frequency of less than or equal to 5% in all combined populations in both the ESP and the 1000G.
    3. Affects function as a modifier: same requirements as the two categories above, but ESP or 1000G frequency must be above 5% and smaller than or equal to 50%.
  2. Variant of Unknown Significance (VUS): literature support for disease associated variant (without functional data) or literature support for deleterious effect of the variant (without disease association). In all cases, the variant must have frequency of less than or equal to 5% in all combined populations in both the ESP and the 1000G.
  3. Likely does not affect function/does not affect function: literature support for lack of pathogenicity; or frequency above 50% in ESP or 1000G for variant that would otherwise be classified in category 1.3; or frequency above 5% in ESP or 1000G for variant that would otherwise be classified in category 2.

Variants that were not found in the HGMD database were classified as VUSs. To filter out variants that would likely not affect protein function, we retained splicing (within 2 bp of splicing junction) and exonic variants (except synonymous) that had a frequency <5% in all combined populations in both the ESP and the 1000G (hg19 data from ANNOVAR [9]), and that had a damage score of 4 or above (see Table G in S1 File). The selected variants above were then subdivided into either being recurrent (found in 2 or more pacemaker patients) or non-recurrent. Recurrent variants that were found to be false positive (FP) in the Coriell samples used for validation were eliminated. Variants that were located within TNXB and ADAMTSL2 and within regions of known segmental duplications were excluded for potentially being FP variants. Frame-shift, premature stop, and splicing VUSs were confirmed by Sanger. For both HGMD and nonHGMD variants, a cut-off for the Minor Allele Frequency (MAF) of 5% was used because in addition to rare and highly penetrant variants, previous reports have implicated common functional genetic variants leading to arrhythmia susceptibility or variable ECG traits including SND. Moreover, it is well known that genetic arrhythmias, although mostly of monogenic and Mendelian basis, may be caused by multiple deleterious DNA variants as it happens in up to 8.5% of subjects with long QT syndrome (LQTS), usually leading to more severe phenotypes [10]. In summary, a 5% MAF cutoff was used to capture rare and highly penetrant deleterious variants, which alone could be responsible for the SND clinical presentation, along with well known or potential more common susceptibility variants, which could play a significant role in the development of SND. Moreover, the recently released “Standards and guidelines for the interpretation of sequence variants” issued by consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology also suggested that allele frequency >5% in Exome Sequencing Project, 1000 Genomes Project, or Exome Aggregation Consortium for a rare Mendelian disorder may be used as stand-alone evidence of benign impact [11]. Provided that PPM occurs at a relative high frequency in the population, that genetic variants may not be highly penetrant, and that multiple deleterious variants may be necessary to cause a clinically relevant presentation, our cut-off of <5% appears to be reasonable.

Information regarding panel design, library preparation, sequencing procedures, bioinformatics pipeline and data analysis, panel validation, additional post-bioinformatics variant filtering procedures, Sanger sequencing for confirmation of variants, and analytic performance assessment are found in S1 File (materials and methods section). Due to the nature of the study, we did not have the permission to contact the affected relatives of the index cases to obtain in-depth clinical information and to enroll them in the genetic study. Therefore, as also explained in the “Limitation Study” section, co-segregation analysis of the putative pathogenic variants was not possible.

Results

Nine faPPM patients were studied using an Agilent HaloPlex custom pan-cardiovascular disease targeted panel for NGS. The NGS panel was designed to detect variants in the splicing and coding regions of 246 genes associated with cardiovascular disorders including cardiomyopathies, arrhythmias, congenital heart defects, aortopathy, connective tissue disorders, Noonan spectrum disorders, pulmonary arterial hypertension, metabolic disorders that afflict the heart, and lipid disorders. Commercially available DNA samples from seven genotype-known individuals from Coriell were sequenced to estimate target coverage and to optimize and validate the panel. Detailed results regarding coverage, sequence depth, inter and intra-run variability using Coriell and pacemaker samples, and the HaloPlex SNP performance for Coriell sequencing runs are found in S1 File (results section). Overall, we obtained consistent coverage within and between MiSeq sequence runs. The precision of the HaloPlex custom panel sequence runs and the variability in sequence depth among runs showed acceptable variability. Validation experiments performed with Coriell samples showed that our custom NGS platform and bioinformatics analysis pipeline can produce data with less than 8% false positive (FP) calls (range of 2.33 to 7.32) and less than 5% false negative (FN) calls (range of 0.63 to 4.23) for exonic and splicing regions of all 246 cardiovascular genes with a minimum of 15x sequencing depth (95% of targets, see Table B in S1 File) after removal of regions of segmental duplications. The overall pacemaker patient target region average sequence depth and the exonic and splicing only (focus regions of our analyses) sequence depth were calculated for our experiments and resulted in an average of 194x and 200x, respectively (Table B in S1 File). Finally, we observed that at a minimum of 15x sequence depth, an average of 95% of target regions were covered with a mean read depth of 205x for pacemaker patients (Table B in S1 File).

Distribution of Variants in ICCD or SSS

Table A in S2 File summarizes the distribution of variants found in each pacemaker patient. An average of 2,344 variants was obtained per patient, with 83.2% of variants mapping to non-exonic and non-splicing positions. Although not intentionally targeted, these intronic, untranslated region (UTR), and intergenic variants come from the library preparation process, which, as expected, does not limit the capture and enrichment process to the exact coordinates of the exons and splicing regions. Additionally, the large occurrence of non-exonic/non-splicing variants is expected given their tendency to be less conserved, supporting the observation that only 0.3% of non-exonic and/or non-splicing variants were documented in the Human Gene Mutation Database (HGMD). This supports the removal of non-HGMD (not present in the HGMD database) non-exonic/non-splicing variants, as described in S1 File (methods section). Likewise, an average of 392 exonic variants was obtained per patient, with nearly 60% of exonic variants being synonymous (silent) and almost 39% of exonic variants being non-synonymous (missense) (Table A in S2 File). Additionally, it was observed that only 2.9% of synonymous variants were found in HGMD, supporting that elimination of non-HGMD synonymous variants, as described in S1 File (methods section), is an acceptable step in the process of finding disease causing variants. A limited number of splicing, stop-gain, stop-loss, frameshift and non-frameshift variants were found in the pacemaker patients.

Classification of Variants in Genotype Unknown ICCD or SSS

Non-HGMD variants were filtered as described in S1 File (materials and methods section) to select a group of variants of unknown significance (VUSs) with potential involvement in the patients’ phenotype (Table C in S1 File). Although these variants have not been previously reported, they are rare in the general population and had a high prediction score of potentially affecting protein function. A minimum of 2 and maximum of 10 non-HGMD VUSs were found in all pacemaker patients (Table D in S1 File). It was hypothesized that genes associated with arrhythmias (Table B in S2 File) may be strong candidates in our pacemaker patients, given their ICCD or SSS presentation and a family history of pacemaker implantation. Eight of our 9 pacemaker patients had at least one VUS in an arrhythmia associated gene (Tables C and D in S1 File).

HGMD documented variants were categorized as described in S1 File (materials and methods section) to come up with a well annotated group of variants. After careful literature review, approximately 33% of HGMD variants were re-classified compared to their original classification (see Table D in S1 File). The final list of re-categorized HGMD VUSs is presented in Table C in S1 File, while the final list of re-categorized HGMD variants that affect protein function in isolation (category 1.1), variants that may or may not affect function in isolation (category 1.2), or variants that affect function as modifiers (category 1.3) is presented in Table E in S1 File. Three out of 9 patients had previously reported pathogenic variants that affect function in isolation (category 1.1); however, 3 of those variants in 2 patients are known to be pathogenic in a homozygous state, whereas the patients were heterozygous (Tables D and E in S1 File). Five out of 9 pacemaker patients had previously reported pathogenic variants that may or may not affect function in isolation (category 1.2), and two of our patients had two of these variants (Tables D and E in S1 File). Overall, patients that had previously reported deleterious variants that could explain their phenotype include those that had at least one variant in category 1.1 or at least two variants in category 1.2. This aggregation resulted in three out of nine of the pacemaker patients (33.3%) that had deleterious variants that could explain their phenotype (patients 4, 6, and 7). Patient 4 had a category 1.1 variant in the LDB3 gene—heterozygous c.349G>A (p.D117N), NM_001080116, rs121908338 (frequency of approximately 1% in the 1000G database and 0.7% in the ESP database), and a category 1.2 variant in the arrhythmia associated gene TRPM4 –heterozygous c.2531G>A (p.G844D), NM_017636, rs200038418 (frequency of approximately 0.05% in the 1000G database and 0.1% in the ESP database). Patient 6 also had the heterozygous TRPM4 above, in addition to a second category 1.2 variant in the arrhythmia associated gene ABCC9—heterozygous c.2200G>A (p.V734I), NM_005691, rs61688134 (frequency of approximately 1% in the 1000G database and 0.9% in the ESP database and 0.7% in the ExAC database). Patient 7 had a category 1.2 variant in the APOB gene—heterozygous c.10580G>A (p.R3527Q), NM_000384, rs5742904 (not found in the 1000G database and 0.07% frequency in the ESP database), and a category 1.2 variant in the arrhythmia associated gene SCN5A—heterozygous c.659C>T (p.T220I), NM_198056, rs45620037 (frequency of approximately 0.05% in the 1000G database and 0.03% in the ESP database).

A comparison of the characteristics of genotype positive and genotype negative patients is also included in Table 1. Notably, all genotype positive patients were Caucasian. While the average age at study recruitment was similar among all groups of patients studied, the genotype positive patients on average had their PPM implanted at an earlier age compared to all patients with pacemakers (49.3 ± 22.5 vs 55.8 ± 16.7 years) (Table 1). No statistically significant differences were found using Fisher’s test between the no family history and the positive family history subjects or between the positive family history and the positive family history + ICCD/SSS subjects for the categories analyzed (ethnicities, dependency, conduction disturbance, complete AVB, sick sinus syndrome, and atrial fibrillation). None of the genotype positive patients had co-morbidities at the time of PPM implantation. A list of major co-morbidities in patients submitted for genetic analysis is listed in Table F. The nine patients with ICCD or SSS without structural heart disease had clinical data available for 18.1 ± 9.7 years after device implantation. Six of them did not develop symptomatic congestive heart failure or coronary artery disease, although echocardiographic abnormalities such as left ventricular hypertrophy or diastolic dysfunction were common. There was one patient in the genotype positive group and two in the genotype negative group who developed coronary artery disease during follow-up. Detailed clinical follow-up data of all patients submitted for genetic analysis is listed in Table 2. We were able to locate ECGs from 5 of the 9 patients studied (see S1 Fig). Among them, the ECGs of patients 2 and 4 showed complete heart block. ECG of patient 6 showed ventricularly paced rhythm (ECG prior to pacemaker implant was not available). ECG of patient 9 showed 2:1 (most likely Mobitz II) AV block. The ECG of patient 8 showed profound sinus bradycardia, junctional rhythm with AV dissociation, consistent with sick sinus syndrome. The ECGs of the remaining 4 patients were not available. Out of the two patients with congenital heart block in the absence of congenital heart disease, one was genotype positive for the LDB3 and TRPM4 variants (patient 4) and the other one was negative. In the patient who was genotype negative (patient 2), while the etiology of the complete heart block was listed as unknown, suspected etiology included both congenital heart block and history of scarlet fever at age 1. Given that no deleterious variants were identified in her case, in retrospect, her conduction disorder could be due to her scarlet fever in childhood.

thumbnail
Table 2. Follow up after pacemaker insertion in patients who had ICCD or SSS without structural heart disease.

https://doi.org/10.1371/journal.pone.0143588.t002

Discussion

Conduction defects represent a major cause of morbidity and mortality worldwide. Pacemaker implantation or a dual device (pacemaker and ICD) implantation remains an appropriate option in subjects presenting with SND symptoms. However, in this study, we have observed that familial aggregation of PPM (faPPM) occurs in 9/112 (8%) patients evaluated in our clinic, thus representing a significant proportion of all individuals requiring a PPM. The occurrence of faPPM suggests a possible genetic predisposition to the development of SND, thus making those subjects candidates for genetic screening. Here, we show that a comprehensive genetic approach using a targeted pan-cardiovascular 246 genes panel could identify at least one previously reported pathogenic category 1.1 or at least two category 1.2 variants as the most likely explanation of their clinical phenotypes as described below.

The LDB3 gene mapping to chromosome 10q23.2 codes for the protein ZASP, which is highly expressed in the heart and skeletal muscle and is a crucial component of the sarcomeric Z-disks in binding critical sarcomeric proteins including alpha actinin-2, protein kinase C, and myozenin family proteins (OMIM 605906) [12, 13]. Variants in LDB3 have been associated with dilated cardiomyopathy (DCM), left ventricular noncompaction (LVNC), and myofibrillar myopathy. The heterozygous LDB3 variant seen in pacemaker patient 4, p.D117N, has been previously seen in heterozygosity in two unrelated individuals with DCM/LVNC and conduction defects [14]. Previous functional analyses in HEK293 and neonatal rat cardiomyocytes demonstrated that ZASP interacts with the SCN5A-coded sodium channel Nav1.5, and affects the channel voltage-dependent activation and inactivation, leading to loss of function of Nav1.5, consistent with the conduction defects seen in the variant harboring patients [15]. We propose that this variant could explain the cardiac phenotype seen in patient 4, alone or in conjunction with the TRPM4 variant described below.

The TRPM4 gene, located on chromosome 19q13.33, encodes for a calcium-activated ion channel highly expressed in the heart, prostate and colon (OMIM 606936) [16]. Variants in the TRPM4 gene have been associated with ICCD, Brugada syndrome, and heart block. The heterozygous TRPM4 variant seen in pacemaker patients 4 and 6, p.G844D, has been previously reported in heterozygosity in multiple families with individuals affected by CCD, right bundle branch block (RBBB) or left anterior hemiblock, with variable penetrance. Functional analysis showed that p.G844D significantly increased current amplitudes in transfected HEK293 cells [17, 18].

The ABCC9 gene maps to chromosome 12p12.1 and encodes the ATP-dependent potassium channel SUR2, which belongs to the family of ABC transporters transporting molecules across membranes (OMIM 601439). Variants in the ABCC9 gene have been associated with DCM, Cantú syndrome, early repolarization syndrome, Brugada syndrome, myocardial infarction (MI), and atrial fibrillation. The heterozygous ABCC9 variant (p.V734I) in pacemaker patient 6 has been previously seen in heterozygosity in individuals with early acute MI [19, 20]. In addition, Dr. Charles Antzelevitch’s group has reported the ABCC9 variant (p.V734I) in four subjects with early repolarization (ERS) and bradycardia. Functional analysis of ABCC9-V734I co-expressed with KCNJ11-WT in human embryonic kidney cell line TSA201 revealed a gain of function in IK-ATP due to a reduced sensitivity of the ATP-sensitive potassium channel (KATP) to ATP, leading to a five-fold increase in KCNJ11-WT function [21]. Moreover, analysis using synchrotron radiation X-ray scattering for the ABCC9 p.V734I variant also demonstrated that it alters protein-protein interaction, which is critical for the structural integrity of the KATP channel complex [21]. Interestingly, one index case harboring the ABCC9 p.V734I variant also harbors a frameshift variant in SCN5A (c.3890_3891insA). The patient, a 40 year old male, presented with sinus bradycardia, first degree atrioventricular block (AVB), ventricular bigeminy and a global ER pattern [21]. Taken together, the functional studies, the variant frequency in major databases and the fact that it was identified along with a putative SCN5A deleterious variant, we have categorize the ABCC9 variant p.V734I as 1.2 because it may not be sufficient in isolation to produce a dramatic clinical effect, but in combination with the TRPM4 variant, it may explain the phenotype of our pacemaker patient 6.

The APOB gene maps to chromosome 2p24.1 and encodes the Apo B-100 protein, which regulates the binding and internalization of low density lipoproteins (LDLs) particles (OMIM 107730). Variants in the APOB gene have been associated with hypercholesterolemia, hypobetalipoproteinemia, apolipoprotein B deficiency, hypertriglyceridemia, and ischemic stroke. The heterozygous APOB variant seen in pacemaker patient 7, p.R3527Q, has been previously seen in heterozygosity in individuals with hypercholesterolemia, MI, and ischemic heart disease (summarized in OMIM 107730). It is known that in rabbit hearts hypercholesterolaemia leads to degeneration of the AV node tissue and the cardiac neural tissue [22]. In addition, aging and hypercholesterolemia in rabbit hearts were associated with the development of atrial tachyarrhythmias [23], while in humans, aging is associated with the increased lipid and collagen deposition around the sinoatrial (SA) node, leading to conduction abnormalities and bradyarrhythmias [24]. However, the role of this variant in the development of conduction disease is currently unknown.

The SCN5A gene mapping to chromosome 3p22.2, encodes the alpha subunit Nav1.5 of the cardiac voltage-gated sodium channel governing the phase 0 of the cardiac action potential (OMIM 600163). Variants in the SCN5A gene have been associated with CCD, SSS, DCM, Brugada syndrome (BrS), long QT syndrome type 3 (LQT3), atrial and ventricular fibrillation, sudden death (adults and infants), arrhythmogenic right ventricular dysplasia, drug-induced arrhythmia, heart block, nodal rhythm defect, and early repolarization syndrome, among others. Previous studies have shown that loss of function variants in SCN5A, are usually associated with BrS and ICCD; whereas gain of function variants are usually associated with the LQT3 [2528]. However, several variants in SCN5A have also been associated with more than one clinical phenotype even within the same family [29]. The heterozygous SCN5A variant seen in pacemaker patient 7, p.T220I, has been previously seen in an individual with congenital SSS that was compound heterozygote for the SCN5A p.T220I and p.R1623X variants (in trans). Functional studies on SCN5A p.T220I in mammalian cell lines showed a small degree of inactivation delay compared to wild-type cells (the p.R1623X variant encoded non-functional sodium channels, almost inexistent in the plasma membrane) [30, 31]. The p.T220I variant was also found in two relatives with DCM [32], an individual with early onset atrial fibrillation [33], as well as an individual with progressive SND who was a compound heterozygote for the SCN5A p.T220I and p.1048SfsX97 variants [34]. Altogether, although the role of the known pathogenic APOB variant in the development of SND is currently unknown, it may be possible that the altered lipid metabolism caused by the APOB variant could have exacerbated the effect of the SCN5A variant. Thus, the combination of the APOB variant and the SCN5A variant could have been sufficient to exert the effect leading to the cardiac phenotype seen in patient 7. However, further studies are warranted to confirm this hypothesis.

Gourraud and colleagues previously studied familial aggregation of PM implantation in the background of families with inherited form of progressive CCD (PCCD), which suggested a possible genetic involvement in this clinical phenotype [35]. In Gourraud’s study, sequencing of 3 genes (SCN5A, SCN1B, and TRPM4) was performed and one previously reported SCN5A variant was detected. In the present report, we have expanded the study of PM implantation irrespective of clinical etiology in the 112 initial patients of our cohort. We found familial aggregation in ICCD and SSS and expanded the genetic analysis to 246 genes using NGS. Overall, of the nine ICCD/SSS pacemaker patients sequenced for 246 genes, three had deleterious variants (category 1.2) in genes (SCN5A and TRPM4) previously associated with CCD or idiopathic SND [5]. Notably, a category 1.2 variant in ABCC9 was found in one of the patients harboring the TRPM4 variant. Besides category 1.1 and 1.2 variants, 44 variants in HGMD were re-categorized as likely genetic modifiers, 59 were re-categorized as VUSs, and 259 were re-categorized to likely not affect function. Although the aforementioned variants, along with the non-HGMD filtered variants (Table C in S1 File), are unlikely to play a significant role in our patients’ phenotype in isolation, they may still cumulatively contribute to the overall clinical presentation or may help reach a threshold for disease manifestation.

Conclusion

NGS is increasingly being established as a powerful molecular genetics tool for the discovery and detection of variants associated with human disease. We have designed and tested an effective NGS panel that has high sequence depth and low FP and FN rates (average of 4.7% and 1.5%, respectively in splicing and exonic regions with at least 15x sequence depth) for exonic and splicing regions of 246 selected cardiovascular genes. The use of this panel in a small cohort of nine ICCD or SSS individuals with faPPM, along with the detection of a molecular cause in 33% of subjects, suggests that genetic testing in the faPPM population may aid in identifying their genetic etiology and in advising other family members at risk.

Study Limitations

Our study has several limitations. A major limitation of the study is that we did not have the permission to contact the affected relatives of the family probands to collect clinical information and specimens for the genetic study. In addition to including only retrospective medical record review, the patient population of our pacemaker clinic is primarily of Caucasian ethnicity; thus validating studies will need to be done in other ethnic groups as well. For ethical reasons we could not contact the relative with the implanted pacemaker to verify the implantation, however the index patient did report the age of implantation of his/her relatives during a follow up phone call, increasing the reliability of the information provided. In the validation of the panel reported, we have not performed sequencing runs with the same individual that would allow for analysis of inter-run variability. We are in the process of implementing such analysis for our clinically available panels.

Supporting Information

S1 Fig. ECG of the study patients.

We were able to locate ECGs from 5 of the 9 patients studied. Among them, the ECGs of patients 2, 4 and 6 showed complete heart block. ECG of patient 9 showed 2:1 AV block and complete heart block. The ECG of patient 8 showed junctional rhythm with AV dissociation. The ECG of the remaining 4 patients were not available.

https://doi.org/10.1371/journal.pone.0143588.s001

(TIF)

S2 Fig. Pair-wise comparison of the coverage of 246 target genes.

Pearson correlation coefficient of the coverage of 246 target genes was calculated for each pair of the experiment. The red dots represent Pearson correlation coefficients between the expected target region coverage and the observed target region coverage of each experiment/sample; the blue dots show correlation coefficients between pairs of samples. The darker blue color indicates overlapped points, the darker the color, the more number of overlapped points.

https://doi.org/10.1371/journal.pone.0143588.s002

(TIF)

S3 Fig. Percentage of target bases covered at the indicated coverage thresholds.

Each color represents results from the sequence run of a Coriell or pacemaker implanted patient.

https://doi.org/10.1371/journal.pone.0143588.s003

(TIF)

S4 Fig. Correlation of reads passing filter (PF) with obtained sequence depth.

https://doi.org/10.1371/journal.pone.0143588.s004

(TIF)

S1 File. Includes supplemental patient data; supplemental materials, methods, and patient selection methodology; supplemental results; supplement references, and 12 tables.

Table A. Exclusion criteria for genetic testing in patients with a first degree relative with a pacemaker. Table B. HaloPlex NGS depth of coverage for Coriell and pacemaker samples. Table C. Selected HGMD and nonHGMD VUSs in pacemaker patients. Table D. Pacemaker variants filtering analysis and classification (number of variants). Table E. HGMD variants with disease association in pacemaker patients. Table F. Major co-morbidities at the time of pacemaker implantation in patients with ICCD or SSS without structural heart disease. Table G. Variant annotation file description. Table H. HGMD initial variant classification. Table I. Primers for Sanger/Big Dye variant confirmation. Table J. HaloPlex intra-run performance. Table K. Overall Performance of SNP variant calling in genotype known Coriell samples. Table L. Selective analyses of SNP performance in Coriell samples.

https://doi.org/10.1371/journal.pone.0143588.s005

(DOCX)

S2 File. Includes 3 tables.

Table A. Pacemaker samples distribution of NGS variants. Table B. Pan-cardiovascular panel genes. Table C. HaloPlex custom target regions.

https://doi.org/10.1371/journal.pone.0143588.s006

(XLSX)

Acknowledgments

This study used DNA samples from the NINDS Repository (catalog.coriell.org/ninds). NINDS Repository sample numbers corresponding to the samples used are: NA12878, NA19240, NA12003, NA19449, NA19982, NA19704 and NA11931. Samples from the Indiana Biobank, which receives government support under a cooperative agreement grant (UL1TR000006) awarded by the National Center for Advancing Translational Research (NCATS) and the Lilly Endowment, were used in this study. Study data were collected and managed using REDCap electronic data capture tools hosted at Indiana University [36]. REDCap (Research Electronic Data Capture) is a secure, web-based application designed to support data capture for research studies, providing 1) an intuitive interface for validated data entry; 2) audit trails for tracking data manipulation and export procedures; 3) automated export procedures for seamless data downloads to common statistical packages; and 4) procedures for importing data from external sources. We would like to thank Lynne Foreman, RN, David Adams, BSEE, David Wagner, BSEE, and Susan Straka, RN for their assistance in obtaining records for these patients.

Author Contributions

Conceived and designed the experiments: MV HAS PSC PBSC. Performed the experiments: PBSC TCL. Analyzed the data: HL HG ZW YL MV PBSC AD AU. Contributed reagents/materials/analysis tools: HL HG ZW YL MV PBSC. Wrote the paper: PBSC AD HAS WJG JMM HG YL PSC MV.

References

  1. 1. Dobrzynski H, Boyett MR, Anderson RH. New insights into pacemaker activity: promoting understanding of sick sinus syndrome. Circulation. 2007; 115(14):1921–1932. pmid:17420362
  2. 2. Michaëlsson M, Jonzon A, Riesenfeld T. Isolated congenital complete atrioventricular block in adult life. A prospective study. Circulation. 1995; 92(3):442–449. pmid:7634461
  3. 3. Smits JP, Veldkamp MW, Wilde AA. Mechanisms of inherited cardiac conduction disease. Europace. 2005; 7(2):122–137. pmid:15763526
  4. 4. Balmer C, Fasnacht M, Rahn M, Molinari L, Bauersfeld U. Long-term follow up of children with congenital complete atrioventricular block and the impact of pacemaker therapy. Europace. 2002; 4(4):345–349. pmid:12408251
  5. 5. Ackerman MJ, Priori SG, Willems S, Berul C, Brugada R, Calkins H, et al. HRS/EHRA expert consensus statement on the state of genetic testing for the channelopathies and cardiomyopathies. Heart Rhythm. 2011; 8:1308–1339. pmid:21787999
  6. 6. Rehm HL. Disease-targeted sequencing: a cornerstone in the clinic. Nat Rev Genet. 2013; 14(4):295–300. pmid:23478348
  7. 7. Rehm HL, Bale SJ, Bayrak-Toydemir P, Berg JS, Brown KK, Deignan JL, et al. ACMG clinical laboratory standards for next-generation sequencing. Genet Med. 2013; 15(9):733–747. pmid:23887774
  8. 8. Human Gene Mutation Database. HGMD Professional 2014.2, BIOBASE, Beverly, MA. Available: http://www.hgmd.org/.
  9. 9. ANNOVAR. Available: http://www.openbioinformatics.org/annovar/annovar_download.html.
  10. 10. Itoh H, Shimizu W, Hayashi K, Yamagata K, Sakaguchi T, Ohno S, et al. Long QT syndrome with compound mutations is associated with a more severe phenotype: a Japanese multicenter study. Heart Rhythm. 2010;7(10):1411–1418. pmid:20541041
  11. 11. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–423. pmid:25741868
  12. 12. Lin X, Ruiz J, Bajraktari I, Ohman R, Banerjee S, Gribble K, et al. Z-disc-associated, alternatively spliced, PDZ motif-containing protein (ZASP) mutations in the actin-binding domain cause disruption of skeletal muscle actin filaments in myofibrillar myopathy. J Biol Chem. 2014; 289(19):13615–13626. pmid:24668811
  13. 13. Frey N, Olson EN. Calsarcin-3, a novel skeletal muscle-specific member of the calsarcin family, interacts with multiple Z-disc proteins. J Biol Chem. 2002; 277(16):13998–14004. pmid:11842093
  14. 14. Vatta M, Mohapatra B, Jimenez S, Sanchez X, Faulkner G, Perles Z, et al. Mutations in Cypher/ZASP in patients with dilated cardiomyopathy and left ventricular non-compaction. J Am Coll Cardiol. 2003; 42(11):2014–2027. pmid:14662268
  15. 15. Xi Y, Ai T, De Lange E, Liz Z, Wu G, Brunelli L, et al. Loss of function of hNav1.5 by a ZASP1 mutation associated with intraventricular conduction disturbances in left ventricular noncompaction. Circ Arrhythm Electrophysiol. 2012; 5(5):1017–1026. pmid:22929165
  16. 16. Xu XZ, Moebius F, Gill DL, Montell C. Regulation of melastatin, a TRP-related protein, through interaction with a cytoplasmic isoform. Proc Natl Acad Sci U S A. 2001; 98(19):10692–10697. pmid:11535825
  17. 17. Liu H, El Zein L, Kruse M, Guinamard R, Beckmann A, Bozio A, et al. Gain-of-function mutations in TRPM4 cause autosomal dominant isolated cardiac conduction disease. Circ Cardiovasc Genet. 2010; 3(4):374–385. pmid:20562447
  18. 18. Stallmeyer B, Zumhagen S, Denjoy I, Duthoit G, Hébert JL, Ferrer X, et al. Mutational spectrum in the Ca(2+)—activated cation channel gene TRPM4 in patients with cardiac conductance disturbances. Hum Mutat. 2012; 33(1):109–117. pmid:21887725
  19. 19. Minoretti P, Falcone C, Aldeghi A, Olivieri V, Mori F, Emanuele E, et al. A novel Val734Ile variant in the ABCC9 gene associated with myocardial infarction. Clin Chim Acta. 2006; 370(1–2):124–128. pmid:16563363
  20. 20. Smith KJ, Chadburn AJ, Adomaviciene A, Minoretti P, Vignali L, Emanuele E, et al. Coronary spasm and acute myocardial infarction due to a mutation (V734I) in the nucleotide binding domain 1 of ABCC9. Int J Cardiol. 2013; 168(4):3506–3513. pmid:23739550
  21. 21. Hu D, Barajas-Martínez H, Terzic A, Park S, Pfeiffer R, Burashnikov E, et al. ABCC9 is a novel Brugada and early repolarization syndrome susceptibility gene. Int J Cardiol. 2014;171(3):431–442. pmid:24439875
  22. 22. Senges J, Weihe E, Brachmann J, Pelzer D, Nimmrich H, Kübler W. Effects of hypercholesterolaemia without ischemia on some electrophysiological and ultrastructural properties of the rabbit heart. J Mol Cell Cardiol. 1981; 13(3):253–264. pmid:7265247
  23. 23. Lee KT, Hsieh CC, Tsai WC, Tang PW, Liu IH, Chai CY, et al. Characteristics of atrial substrates for atrial tachyarrhythmias induced in aged and hypercholesterolemic rabbits. Pacing Clin Electrophysiol. 2012; 35(5):544–552. pmid:22443495
  24. 24. Chow GV, Marine JE, Fleg JL. Epidemiology of arrhythmias and conduction disorders in older adults. Clin Geriatr Med. 2012; 28(4):539–553. pmid:23101570
  25. 25. Wang Q, Shen J, Splawski I, Atkinson D, Li Z, Robinson JL, et al. SCN5A mutations associated with an inherited cardiac arrhythmia, long QT syndrome. Cell. 1995; 80(5):805–811. pmid:7889574
  26. 26. Chen Q, Kirsch GE, Zhang D, Brugada R, Brugada J, Brugada P, et al. Genetic basis and molecular mechanism for idiopathic ventricular fibrillation. Nature. 1998; 392(6673):293–296. pmid:9521325
  27. 27. Schott JJ, Alshinawi C, Kyndt F, Probst V, Hoorntje TM, Hulsbeek M, et al. Cardiac conduction defects associate with mutations in SCN5A. Nat Genet. 1999; 23(1):20–21. pmid:10471492
  28. 28. Tan HL, Bink-Boelkens MT, Bezzina CR, Viswanathan PC, Beaufort-Krol GC, van Tintelen PJ, et al. A sodium-channel mutation causes isolated cardiac conduction disease. Nature. 2001; 409(6823):1043–1047. pmid:11234013
  29. 29. Tan BH, Iturralde-Torres P, Medeiros-Domingo A, Nava S, Tester DJ, Valdivia CR, et al. A novel C-terminal truncation SCN5A mutation from a patient with sick sinus syndrome, conduction disorder and ventricular tachycardia. Cardiovasc Res. 2007; 76(3):409–417. pmid:17897635
  30. 30. Benson DW, Wang DW, Dyment M, Knilans TK, Fish FA, Strieper MJ, et al. Congenital sick sinus syndrome caused by recessive mutations in the cardiac sodium channel gene (SCN5A). J Clin Invest. 2003; 112(7):1019–1028. pmid:14523039
  31. 31. Gui J, Wang T, Jones RP, Trump D, Zimmer T, Lei M. Multiple loss-of-function mechanisms contribute to SCN5A-related familial sick sinus syndrome. PLoS One. 2010; 5(6):e10985. pmid:20539757
  32. 32. Olson TM, Michels VV, Ballew JD, Reyna SP, Karst ML, Herron KJ, et al. Sodium channel mutations and susceptibility to heart failure and atrial fibrillation. JAMA. 2005; 293(4):447–454. pmid:15671429
  33. 33. Olesen MS, Yuan L, Liang B, Holst AG, Nielsen N, Nielsen JB, et al. High prevalence of long QT syndrome-associated SCN5A variants in patients with early-onset lone atrial fibrillation. Circ Cardiovasc Genet. 2012; 5(4):450–459. pmid:22685113
  34. 34. Baskar S, Ackerman MJ, Clements D, Mayuga KA, Aziz PF. Compound heterozygous mutations in the SCN5A-encoded Nav1.5 cardiac sodium channel resulting in atrial standstill and His-Purkinje system disease. J Pediatr. 2014; 165(5):1050–1052. pmid:25171853
  35. 35. Gourraud JB, Kyndt F, Fouchard S, Rendu E, Jaafar P, Gully C, et al. Identification of a strong genetic background for progressive cardiac conduction defect by epidemiological approach. Heart. 2012;98(17):1305–1310. pmid:22717692
  36. 36. Harris PA, Taylor R, Thielke R, Payne J, Gonzalez N, Conde JG. Research electronic data capture (REDCap)—a metadata-driven methodology and workflow process for providing translational research informatics support. J Biomed Inform. 2009; 42(2):377–381. pmid:18929686