Peter Weeke

Common variants at SCN5A-SCN10A and HEY2 are associated with Brugada syndrome, a rare disease with high risk of sudden cardiac death

Brugada syndrome is a rare cardiac arrhythmia disorder, causally related to SCN5A mutations in around 20% of cases. Through a genome-wide association study of 312 individuals with Brugada syndrome and 1,115 controls, we detected 2 significant association signals at the SCN10A locus (rs10428132) and near the HEY2 gene (rs9388451). Independent replication confirmed both signals (meta-analyses: rs10428132, P = 1.0 × 10−68; rs9388451, P = 5.1 × 10−17) and identified one additional signal in SCN5A (at 3p21; rs11708996, P = 1.0 × 10−14). The cumulative effect of the three loci on disease susceptibility was unexpectedly large (Ptrend = 6.1 × 10−81). The association signals at SCN5A-SCN10A demonstrate that genetic polymorphisms modulating cardiac conduction can also influence susceptibility to cardiac arrhythmia. The implication of association with HEY2, supported by new evidence that Hey2 regulates cardiac electrical activity, shows that Brugada syndrome may originate from altered transcriptional programming during cardiac development. Altogether, our findings indicate that common genetic variation can have a strong impact on the predisposition to rare diseases.

Role of common and rare variants in SCN10A: results from the Brugada syndrome QRS locus gene discovery collaborative study

AIMS Brugada syndrome (BrS) remains genetically heterogeneous and is associated with slowed cardiac conduction. We aimed to identify genetic variation in BrS cases at loci associated with QRS duration. METHODS AND RESULTS A multi-centre study sequenced seven candidate genes (SCN10A, HAND1, PLN, CASQ2, TKT, TBX3, and TBX5) in 156 Caucasian SCN5A mutation-negative BrS patients (80% male; mean age 48) with symptoms (64%) and/or a family history of sudden death (47%) or BrS (18%). Forty-nine variants were identified: 18 were rare (MAF <1%) and non-synonymous; and 11/18 (61.1%), mostly in SCN10A, were predicted as pathogenic using multiple bioinformatics tools. Allele frequencies were compared with the Exome Sequencing and UK10K Projects. SKAT methods tested rare variation in SCN10A finding no statistically significant difference between cases and controls. Co-segregation analysis was possible for four of seven probands carrying a novel pathogenic variant. Only one pedigree (I671V/G1299A in SCN10A) showed co-segregation. The SCN10A SNP V1073 was, however, associated strongly with BrS [66.9 vs. 40.1% (UK10K) OR (95% CI) = 3.02 (2.35-3.87), P = 8.07 × 10-19]. Voltage-clamp experiments for NaV1.8 were performed for SCN10A common variants V1073, A1073, and rare variants of interest: A200V and I671V. V1073, A200V and I671V, demonstrated significant reductions in peak INa compared with ancestral allele A1073 (rs6795970). CONCLUSION Rare variants in the screened QRS-associated genes (including SCN10A) are not responsible for a significant proportion of SCN5A mutation negative BrS. The common SNP SCN10A V1073 was strongly associated with BrS and demonstrated loss of NaV1.8 function, as did rare variants in isolated patients.

Biobanks and Electronic Medical Records: Enabling Cost-Effective Research

Linking of data from electronic medical records to biological specimens enables cost-effective and rapid genomic analyses. The use of electronic medical record data linked to biological specimens in health care settings is expected to enable cost-effective and rapid genomic analyses. Here, we present a model that highlights potential advantages for genomic discovery and describe the operational infrastructure that facilitated multiple simultaneous discovery efforts.

Exome sequencing implicates an increased burden of rare potassium channel variants in the risk of drug-induced long QT interval syndrome.

OBJECTIVES The aim of this study was to test the hypothesis that rare variants are associated with drug-induced long QT interval syndrome (diLQTS) and torsades de pointes. BACKGROUND diLQTS is associated with the potentially fatal arrhythmia torsades de pointes. The contribution of rare genetic variants to the underlying genetic framework predisposing to diLQTS has not been systematically examined. METHODS We performed whole-exome sequencing on 65 diLQTS patients and 148 drug-exposed control subjects of European descent. We used rare variant analyses (variable threshold and sequence kernel association test) and gene-set analyses to identify genes enriched with rare amino acid coding (AAC) variants associated with diLQTS. Significant associations were reanalyzed by comparing diLQTS patients with 515 ethnically matched control subjects from the National Heart, Lung, and Blood Grand Opportunity Exome Sequencing Project. RESULTS Rare variants in 7 genes were enriched in the diLQTS patients according to the sequence kernel association test or variable threshold compared with drug-exposed controls (p < 0.001). Of these, we replicated the diLQTS associations for KCNE1 and ACN9 using 515 Exome Sequencing Project control subjects (p < 0.05). A total of 37% of the diLQTS patients also had 1 or more rare AAC variants compared with 21% of control subjects (p = 0.009), in a pre-defined set of 7 congenital long QT interval syndrome (cLQTS) genes encoding potassium channels or channel modulators (KCNE1, KCNE2, KCNH2, KCNJ2, KCNJ5, KCNQ1, AKAP9). CONCLUSIONS By combining whole-exome sequencing with aggregated rare variant analyses, we implicate rare variants in KCNE1 and ACN9 as risk factors for diLQTS. Moreover, diLQTS patients were more burdened by rare AAC variants in cLQTS genes encoding potassium channel modulators, supporting the idea that multiple rare variants, notably across cLQTS genes, predispose to diLQTS.

Rare genetic variants previously associated with congenital forms of long QT syndrome have little or no effect on the QT interval

AIMS We studied whether variants previously associated with congenital long QT syndrome (cLQTS) have an effect on the QTc interval in a Danish population sample. Furthermore, we assessed whether carriers of variants in cLQTS-associated genes are more prone to experience syncope compared with non-carriers and whether carriers have an increased mortality compared with non-carriers. METHODS AND RESULTS All genetic variants previously associated with cLQTS were surveyed using the Human Gene Mutation Database. We screened a Danish population-based sample with available whole-exome sequencing data (n = 870) and genotype array data (n = 6161) for putative cLQTS genetic variants. In total, 33 of 1358 variants previously reported to associate with cLQTS were identified. Of these, 10 variants were found in 8 or more individuals. Electrocardiogram results showed normal mean QTc intervals in carriers compared with non-carriers. Syncope data analysis between variant and non-variant carriers showed that 4 of 227 (1.8%) and 95 of 5861 (1.6%) individuals, respectively, had experienced syncope during follow-up (P = 0.80). There was no significant difference in overall mortality rates between carriers [7/217 (3.2%)] and non-carriers [301/6453 (4.7%)] (P = 0.24). CONCLUSION We present QTc data and register data, indicating that 26 cLQTS-associated variants neither had any effect on the QTc intervals nor on syncope propensity or overall mortality. Based on the frequency of individual gene variants, we suggest that the 10 variants frequently identified, assumed to relate to cLQTS, are less likely to associate with a dominant monogenic form of the disease.

SCN10A/Nav1.8 modulation of peak and late sodium currents in patients with early onset atrial fibrillation

AIMS To test the hypothesis that vulnerability to atrial fibrillation (AF) is associated with rare coding sequence variation in the SCN10A gene, which encodes the voltage-gated sodium channel isoform NaV1.8 found primarily in peripheral nerves and to identify potentially disease-related mechanisms in high-priority rare variants using in-vitro electrophysiology. METHODS AND RESULTS We re-sequenced SCN10A in 274 patients with early onset AF from the Vanderbilt AF Registry to identify rare coding variants. Engineered variants were transiently expressed in ND7/23 cells and whole-cell voltage clamp experiments were conducted to elucidate their functional properties. Resequencing SCN10A identified 18 heterozygous rare coding variants (minor allele frequency ≤1%) in 18 (6.6%) AF probands. Four probands were carriers of two rare variants each and 14 were carriers of one coding variant. Based on evidence of co-segregation, initial assessment of functional importance, and presence in ≥1 AF proband, three variants (417delK, A1886V, and the compound variant Y158D-R814H) were selected for functional studies. The 417delK variant displayed near absent current while A1886V and Y158D-R814H exhibited enhanced peak and late (INa-L) sodium currents; both Y158D and R818H individually contributed to this phenotype. CONCLUSION Rare SCN10A variants encoding Nav1.8 were identified in 6.6% of patients with early onset AF. In-vitro electrophysiological studies demonstrated profoundly altered function in 3/3 high-priority variants. Collectively, these data strongly support the hypothesis that rare SCN10A variants may contribute to AF susceptibility.

Whole exome sequencing identifies a causal RBM20 mutation in a large pedigree with familial dilated cardiomyopathy.

Background—Whole exome sequencing is a powerful technique for Mendelian disease gene discovery. However, variant prioritization remains a challenge. We applied whole exome sequencing to identify the causal variant in a large family with familial dilated cardiomyopathy of unknown pathogenesis. Methods and Results—A large family with autosomal dominant, familial dilated cardiomyopathy was identified. Exome capture and sequencing were performed in 3 remotely related, affected subjects predicted to share <0.1% of their genomes by descent. Shared variants were filtered for rarity, evolutionary conservation, and predicted functional significance, and remaining variants were filtered against 71 locally generated exomes. Variants were also prioritized using the Variant Annotation Analysis and Search Tool. Final candidates were validated by Sanger sequencing and tested for segregation. There were 664 shared heterozygous nonsense, missense, or splice site variants, of which 26 were rare (minor allele frequency ⩽0.001 or not reported) in 2 public databases. Filtering against internal exomes reduced the number of candidates to 2, and of these, a single variant (c.1907 G>A) in RBM20, segregated with disease status and was absent in unaffected internal reference exomes. Bioinformatic prioritization with Variant Annotation Analysis and Search Tool supported this result. Conclusions—Whole exome sequencing of remotely related dilated cardiomyopathy subjects from a large, multiplex family, followed by systematic filtering, identified a causal RBM20 mutation without the need for linkage analysis.

Genetic investigations of sudden unexpected deaths in infancy using next-generation sequencing of 100 genes associated with cardiac diseases

Sudden infant death syndrome (SIDS) is the most frequent manner of post-perinatal death among infants. One of the suggested causes of the syndrome is inherited cardiac diseases, mainly channelopathies, that can trigger arrhythmias and sudden death. The purpose of this study was to investigate cases of sudden unexpected death in infancy (SUDI) for potential causative variants in 100 cardiac-associated genes. We investigated 47 SUDI cases of which 38 had previously been screened for variants in RYR2, KCNQ1, KCNH2 and SCN5A. Using the Haloplex Target Enrichment System (Agilent) and next-generation sequencing (NGS), the coding regions of 100 genes associated with inherited channelopathies and cardiomyopathies were captured and sequenced on the Illumina MiSeq platform. Sixteen (34%) of the SUDI cases had variants with likely functional effects, based on conservation, computational prediction and allele frequency, in one or more of the genes screened. The possible effects of the variants were not verified with family or functional studies. Eight (17%) of the SUDI cases had variants in genes affecting ion channel functions. The remaining eight cases had variants in genes associated with cardiomyopathies. In total, one third of the SUDI victims in a forensic setting had variants with likely functional effect that presumably contributed to the cause of death. The results support the assumption that channelopathies are important causes of SUDI. Thus, analysis of genes associated with cardiac diseases in SUDI victims is important in the forensic setting and a valuable supplement to the clinical investigation in all cases of sudden death.

A genome-wide association study of heparin-induced thrombocytopenia using an electronic medical record

Heparin-induced thrombocytopenia (HIT) is an unpredictable, potentially catastrophic adverse effect of heparin treatment resulting from an immune response to platelet factor 4 (PF4)/heparin complexes. No genome-wide evaluations have been performed to identify potential genetic influences on HIT. Here, we performed a genome-wide association study (GWAS) and candidate gene study using HIT cases and controls identified using electronic medical records (EMRs) coupled to a DNA biobank and attempted to replicate GWAS associations in an independent cohort. We subsequently investigated influences of GWAS-associated single nucleotide polymorphisms (SNPs) on PF4/heparin antibodies in non-heparin treated individuals. In a recessive model, we observed significant SNP associations (odds ratio [OR] 18.52; 95% confidence interval [CI] 6.33-54.23; p=3.18×10(-9)) with HIT near the T-Cell Death-Associated Gene 8 (TDAG8). These SNPs are in linkage disequilibrium with a missense TDAG8 SNP. TDAG8 SNPs trended toward an association with HIT in replication analysis (OR 5.71; 0.47-69.22; p=0.17), and the missense SNP was associated with PF4/heparin antibody levels and positive PF4/heparin antibodies in non-heparin treated patients (OR 3.09; 1.14-8.13; p=0.02). In the candidate gene study, SNPs at HLA-DRA were nominally associated with HIT (OR 0.25; 0.15-0.44; p=2.06×10(-6)). Further study of TDAG8 and HLA-DRA SNPs is warranted to assess their influence on the risk of developing HIT.

Whole-exome sequencing in familial atrial fibrillation

AIMS Positional cloning and candidate gene approaches have shown that atrial fibrillation (AF) is a complex disease with familial aggregation. Here, we employed whole-exome sequencing (WES) in AF kindreds to identify variants associated with familial AF. METHODS AND RESULTS WES was performed on 18 individuals in six modestly sized familial AF kindreds. After filtering very rare variants by multiple metrics, we identified 39 very rare and potentially pathogenic variants [minor allele frequency (MAF) ≤0.04%] in genes not previously associated with AF. Despite stringent filtering >1 very rare variants in the 5/6 of the kindreds were identified, whereas no plausible variants contributing to familial AF were found in 1/6 of the kindreds. Two candidate AF variants in the calcium channel subunit genes (CACNB2 and CACNA2D4) were identified in two separate families using expression data and predicted function. CONCLUSION By coupling family data with exome sequencing, we identified multiple very rare potentially pathogenic variants in five of six families, suggestive of a complex disease mechanism, whereas none were identified in the remaining AF pedigree. This study highlights some important limitations and challenges associated with performing WES in AF including the importance of having large well-curated multi-generational pedigrees, the issue of potential AF misclassification, and limitations of WES technology when applied to a complex disease.