Abstract

Atrial fibrillation (AF) is a common arrhythmia for which the genetic studies mainly focused on the genes involved in electrical remodeling, rather than left atrial muscle remodeling. To identify rare variants involved in atrial myopathy using mutational screening, a high-throughput next-generation sequencing (NGS) workflow was developed based on a custom AmpliSeq™ panel of 55 genes potentially involved in atrial myopathy. This workflow was applied to a cohort of 94 patients with AF, 76 with atrial dilatation and 18 without. Bioinformatic analyses used NextGENe® software and in silico tools for variant interpretation. The AmpliSeq custom-made panel efficiently explored 96.58% of the targeted sequences. Based on in silico analysis, 11 potentially pathogenic missense variants were identified that were not previously associated with AF. These variants were located in genes involved in atrial tissue structural remodeling. Three patients were also carriers of potential variants in prevalent arrhythmia-causing genes, usually associated with AF. Most of the variants were found in patients with atrial dilatation (n=9, 82%). This NGS approach was a sensitive and specific method that identified 11 potentially pathogenic variants, which are likely to play roles in the predisposition to left atrial myopathy. Functional studies are needed to confirm their pathogenicity.

1. Introduction

Atrial fibrillation (AF) is the most frequent arrhythmia, affecting 30 million individuals worldwide [1]. Advanced age and hypertension, which can damage the left atrium (LA), are the main predisposing risk factors for AF [2]. A plethora of evidence suggests that the onset of most AF types is facilitated by LA remodeling, i.e., atrial myopathy [3]. Ion-channel, neural, and structural remodeling of the LA muscle has been widely documented [4] and numerous studies have found a genetic predisposition and a highly heritable component associated with AF risk [5].

In the past 20 years, the genetic basis for AF was established through studies evaluating familial AF [6, 7], linkage [8, 9], candidate genes [10, 11], and genome-wide association studies (GWAS) [1214] that reported common and rare variants in genes encoding ion-channels, gap junction proteins, and signaling molecules. Recently, next-generation sequencing (NGS) technologies have advanced in terms of sensibility, specificity, practicability, and the cost to rapidly screen large numbers of genes. Massively parallel NGS approaches, including gene panels, whole exome sequencing, or whole genome sequencing, are beginning to supplant Sanger sequencing [15]. Thus, sequencing candidate genes might be the best approach to reveal variations in AF-associated genes [1618].

The available molecular data only account for a limited percentage of the genes involved in AF, mainly those involved in ion-channel remodeling. Atrial myocardial damage is characterized by atrial fibrosis [19], inflammatory infiltrates [20], altered cell-to-cell adhesion and mechanical coupling [21], and abnormal contractions [22]. To identify variants in the genes coding for proteins potentially involved in atrial tissue rather than ion-channel remodeling, we designed a fast protocol utilizing a custom AmpliSeq panel and Ion Personal Genome Machine (PGM) Sequencer to sequence 55 atrial myopathy candidate genes in a prospective cohort of 94 patients, 76 with and 18 without atrial dilatation. Patients carrying pathogenic or likely pathogenic variants were also screened against a homemade panel of prevalent arrhythmia-causing genes, mainly involved in electrical remodeling.

2. Materials and Methods

2.1. Patients

The cohort included 94 nonvalvular patients with AF prospectively recruited from the Louis Pradel Cardiology Hospital (Hospices Civils de Lyon, Lyon, France). The ethics committee of Lyon approved the study and informed consent was obtained from each patient prior to enrollment (DC2015-2566). Individuals older than 18 years with a confirmed diagnosis of paroxysmal/persistent/permanent AF but without significant underlying heart disease, left ventricular dysfunction (left ventricular ejection fraction <50%), valvular heart diseases, or other systemic/metabolic diseases were included in the cohort. The presence of AF was determined by ECG or Holter recordings. Paroxysmal AF was defined as self-terminating, in most cases within 48 hours. Persistent AF lasted longer than 7 days and was terminated by either pharmacologic intervention or electrical cardioversion. For permanent AF, rhythm control interventions were not pursued [23]. Left atrial dilatation was defined as a volume superior to 32 ml/m2 or a surface > 22 cm2, measured by transthoracic echocardiography.

2.2. NGS Strategy

Genomic DNA samples underwent NGS using a custom AmpliSeq design (Life Technologies, Carlsbad, CA, USA) created using Ion AmpliSeq designer software. In the first step, the criteria for gene selection were based on the previously reported transcriptome of atrial tissue in patients with AF [24]. We found that 1,627 genes had altered basal expression levels in the LA tissue of patients with AF compared with the control group. The significantly enriched Gene Ontology biological process “anatomical structure morphogenesis” contained the highest number of genes, and this was in line with changes in structure that occur when the human heart remodels following AF development (i.e., left atrial dilatation and interstitial fibrosis). We then selected the most dysregulated genes to build a homemade gene panel. In the second step, genes were selected, using PubMed, based on their documented or potential involvement in structural remodeling. A candidate ID gene list was generated using the search terms: “structural remodeling”, “AF fibrosis”, “AF conduction”, and “AF inflammation”. Articles concerning the structural remodeling of AF were included predominantly in the list. The genetic panel was made of 55 genes potentially involved in structural heart disease (Table 1). The design allowed analysis of all coding exons of selected genes (padding ±30 bp). Library preparation and Ion Torrent PGM sequencing were performed as previously reported [25, 26]. Selected patients carrying pathogenic or likely pathogenic variants in this panel were further tested by NGS using a second custom panel designed to identify disease-causing variants in 38 known arrhythmia-causing genes [27].

2.3. Bioinformatic Analyses

Bioinformatic analyses were performed using a homemade pipeline based on NextGENe v.2.3.4.2 (SoftGenetics, State College, PA, USA) and Alamut® 2.7.1 (Interactive Biosoftware, Rouen, France) software, as previously reported [25, 26]. Identified gene variants (i.e., missense, nonsynonymous, splice site, insertions, and deletions) were further analysed using the filtering steps shown in Figure 1. According to reported guidelines, specific standard terms [“pathogenic”, “likely pathogenic”, “uncertain significance”, “likely benign”, and “benign”] were used to evaluate the pathogenicity of variants identified in the studied genes [28]. Likely pathogenic variants (single nucleotide variants) were verified by conventional dideoxy sequencing using the BigDye® Terminator v.3.1 Cycle Sequencing Kit (Life Technologies, Carlsbad, CA, USA) and an ABI 3730 automatic sequencer (Life Technologies).

The frequency of each variant in the general population was examined using the disease database ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/) and the population database Exome Aggregation Consortium (ExAC) (http://exac.broadinstitute.org/). In silico tools used for missense variant interpretation included PolyPhen-2 [29], SIFT [30], and MutationTaster [31]. The grade of evolutionary nucleotide conservation was determined by PhyloP scores (http://compgen.cshl.edu/phast/). The protein evolution was predicted with the Grantham score [32]. The protein domains affected by the single nucleotide changes were also described. Multiple protein sequences across species were aligned using the program MUSCLE [33] version 3.6.

2.4. Quantification Methods

Nuclear positioning was quantified in mammalian myotubes containing at least five nuclei, and myotubes were classified aggregated when more than 70% of the nuclei did not align along the same axis.

2.5. Transfections

Myoblasts were transfected with siRNA using Lipofectamine RNAiMAX (Invitrogen): siRNA sequences (Ambion): 5′-GCUCUAAACAUGAUUCAAGTT-3′ (AKPA9-#1); 5′-CGAUGGUAGAAUUCUUAGATT-3′ (AKPA9-#2); 5′-GCCAAGCUUGUCCAUUGAU-TT-3′ (AKPA9-#3).

2.6. Cell Culture

C2C12 myoblasts were grown and differentiated for 5 days as described before [34].

2.7. Statistical Analysis

Student’s t-tests were performed. Differences were considered statistically significant when P< 0.01.

3. Results

Clinical features of the 63 men and 31 women included in the cohort are listed in Table 2. The median age at the time of inclusion for AF probands was 54.4 years (range: 42–66 years). Paroxysmal AF was the most common type and 80.8% of patients with AF presented with left atrial dilatation. Particularly, patients developing permanent AF presented left atrial dilatation. Our AmpliSeq custom-made panel explored 96.58% of targeted sequences. Six runs, containing 16 DNA samples each, were performed and the coverage statistics were comparable between each run. The strategy for filtering (Figure 1) led to the identification of 11 putative pathogenic variants not previously reported in patients with AF (Table 3). Each variant was present in a single patient. Nine variants were found in patients with AF and left atrial dilatation and two in patients without atrial myopathy. Three variants were not reported in the ExAC consortium. All putative pathogenic missense variants were predicted to disrupt protein function by PolyPhen-2 (score ranges: 0 to 1), SIFT, and MutationTaster as “probably damaging” (0.85 = the threshold), “deleterious”, and “disease causing”, respectively. The PhyloP highlighted that mutated nucleotides detected in cases 6211, 4464, 2095, 1885, 4162, 1875, and 2691 were highly conserved. Comparisons based on the physical or chemical properties of amino acids showed the candidate JPH2 variant (p.Ser255Leu) and MMP9 variant (p.Arg143Cys) had high Grantham differences, suggesting that these missense variants could be pathogenic. A multialignment of proteins (Figure 2) showed that all altered amino acids had high evolutionary conservation across species, suggesting that they could be functionally important.

The 11 patients identified with variants involved in structural remodeling were further screened using an arrhythmia panel with genes known to be associated with AF [27]. Three of these patients were also carriers of likely pathogenic variants in AF-associated genes (Table 4). Left atrial dilatation was also a characteristic of these patients. Only eight patients were carriers of likely pathogenic variants in atrial myopathy genes. An overview of AF-associated genes is displayed in Table 5. The majority of these genes are linked with other cardiac diseases. The cellular localization of proteins encoded by candidate genes is shown in Figure 3.

AKAP9 encodes a scaffolding protein involved in Golgi apparatus integrity and Golgi-related microtubules nucleation [35]. It has been recently shown that AKAP9 can contribute to recruit microtubule-organizing center factors at the membrane of myonuclei [36]. We validated AKAP9-dependent myonuclei positioning in a muscle cells context using C2C12 myoblast and quantify myonuclei aggregation in AKAP9-depleted myotubes using 3 different siRNA (Figure 4). AKAP9-depleted myotubes significantly increase myonuclei aggregation phenotype (up to 30%) within myotubes (Figure 4(c)) without affecting myoblast fusion or myotubes differentiation (Figures 4(a) and 4(b)), confirming a microtubule integrity regulation by an AKAP9-dependant mechanism in a muscle cells context [36].

4. Discussion

This study identified 11 potentially pathogenic variants in patients with AF, using a simple and fast NGS mutation detection approach. In contrast with previous studies, our method focused on the identification of candidate gene variants not previously linked to AF-structural remodeling genes. The role of genetic factors in the development of AF, a complex and multifactorial arrhythmia, is increasingly recognized. At least 14 genetic loci revealed by GWAS are known to increase the risk of AF in populations [37], but these variants only explain a small fraction of the interindividual risk for AF. Most identified genetic loci are associated with genes of electrical remodeling, such as KCNN3 [13], or developmental genes, such as PITX2 [12]. However, a meta-analysis of GWAS suggested additional candidate AF loci, such as genes involved in structural components (SYNE2, MYOZ1, and SYNPO2L) [14]. The NGS represents a high-throughput, rapid, and low-cost strategy for the systematic detection of genomic variants involved in AF. Our NGS approach was based on a custom AmpliSeq design to detect variants in structural remodeling genes. The filtering strategy allowed us to identify 11 rare variants. For all variants, in silico tools were used to predict the possible pathogenic impact of an amino acid substitution on the structure and function of the human proteins. This predicted deleterious impact of these variants was strengthened by the evolutionary conservation of the altered amino acids.

Our initial hypothesis was that structural genes could be involved in atrial remodeling as much as ion-channel ones. Three likely pathogenic variants were in ion-channel genes previously associated with AF. Defects were found in ANK2, which encodes a multifunctional cytoskeletal adaptor [38], KCNH2, which encodes a potassium voltage-gated channel, and SCN1B, which encodes the β-subunit of the sodium channel [39]. Evaluation of the missense variants using both segregation data and in vitro systems may help better understand the pathogenicity. The substitution at the splice donor site of the SCN1B intron 1, which was not reported in the ExAC consortium, is expected to yield a nonsense-mediated decay mechanism, resulting in a reduction of protein and haploinsufficiency. Several studies have shown that atrial dilatation is an independent risk factor for the development of AF [40]. In a recent study of eight patients with AF and a frameshift deletion in MYL4, six subjects developed LA dilatation during the follow-up [22]. In the present study, 82% of the novel variants were found in patients with LA dilatation, reinforcing the suggestion that these variants could be involved in LA structural damage. Most of the identified genes were previously linked to other cardiac diseases (Table 5). AKAP9, FHOD3, and TMEM43 were not previously associated with AF in the literature, but they were linked with other cardiac diseases.

The majority of the new variants found in the present study are located in genes encoding a broad category of proteins. These proteins are involved in many diverse biological processes related to structural remodeling of the extracellular matrix, the sarcolemma, the cytoskeleton, desmosome, sarcomere, the sarcoplasmic reticulum, and nucleus. Upregulation of MMP9, a profibrotic and proinflammatory molecule, contributes to atrial extracellular matrix remodeling [41], which is associated with the development of AF [42]. In the sarcolemmal ATP-sensitive potassium channels of the cardiomyocytes, ABCC8 encodes the regulatory sulfonylurea receptor 1. Proteins involved in the desmosome structure include that encoded by DSG2 and DSP. DSG2 is more expressed in LA of patients with AF than control subjects as previously described [24]. Transcriptional network of cardiac rhythm driven by TBX5 and modulated by PITX2 regulates Scn5a, Gja1, Ryr2, Dsp, and Atp2a2 genes [43]. Some of the proteins associated with the selected variants contribute to the structure or function of the sarcomere, with FHOD3 playing a role in regulation of the actin filament assembly [44]. The cell structure gene MYOZ1 encodes myozenin-1, which is a skeletal muscle Z line protein involved in stabilizing the sarcomere [45]. In addition, JPH2 encodes a cardiac structural protein contributing to the formation of the junctional membrane complex architecture that links the sarcoplasmic reticulum with the plasma membrane in cardiomyocytes [46]. The JPH2 mutation is thought to cause AF because of impaired stabilization of ryanodine receptor Ca2+ channels [47]. The inner nuclear membrane contains associated proteins, including that encoded by TMEM43, which is associated with lamin A/C and emerin [48]. AKAP9, a scaffolding protein involved in Golgi apparatus integrity and Golgi-related microtubules nucleation [35], is known to be the long QT syndrome-causative gene [49]. Our results confirmed an altered microtubule network in absence of AKAP9 as inhibition of AKAP9 results in increased aggregation phenotype in myotubes [36]. Consequences of AKAP9 knockdown on remaining pool of microtubule-associate-partners remain to be determines. One can speculate that forces exerted by muscle molecular motors could be remodel in absence (or mutated forms) of AKAP9 and could contribute to alteration of microtubule network dynamic [50, 51]. Microtubules networks are mechanically involved in cardiomyocyte contraction [52]. It will be of interest to analyse resulting network depending on different AKAP9 variant and skeletal muscle cells could be used as a «simplified muscle model» to screen for the effect on microtubule dynamics of different variant of AKAP9 found in cardiac muscles.

Each of these variants is involved in different pathways. The link between these variants and the effect on gene expression is unclear. A recent study has found that the SNP rs2595104 associated with AF regulates PITX2c expression via interaction with TFAP2a [53]. MiRNAs are part of the molecular alterations in LA occurring in patients with atrial remodeling [54]. One might consider that a variant could regulate miRNA in AF patients [55]. Cumulative evidence suggests that response to therapy may be genotype dependent. For example, SNP on chromosome 4q25 associated with AF modulates response to antiarrhythmic therapy [56]. This work opens research directions to establish personalised therapies according to individual genomic data as in cancer patients [57].

5. Conclusions

Eleven rare or novel potentially pathogenic variants were identified using the NGS method in patients with nonvalvular AF, mainly in those with atrial dilatation. Validation studies are needed to confirm the involvement of these variants in atrial structural remodeling. This approach (Figure S1), based on genes involved in atrial structural remodeling, may help uncover new mechanisms underlying AF. In addition, candidate gene approaches based on disease physiopathology should be encouraged.

List of Abbreviations

AF:Atrial fibrillation
ExAC:Exome aggregation consortium
GWAS:Genome-wide association studies
LA:Left atrium
NGS:Next-generation sequencing
PGM:Personal Genome Machine.

Data Availability

The sequencing data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

An earlier version of this work was presented at Printemps de la Cardiologie 2018, 13th European Cardiac Arrhythmia Society Congress, CNIC Conference “Atrial fibrillation: from mechanisms to population science,” and 18th Annual Cardiologists Conference.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Authors’ Contributions

All the authors provided substantial contributions to the conception, design, and acquisition and analysis of the data. Pr. Philippe Chevalier is responsible for the integrity of the work.

Acknowledgments

The authors would like to thank the Laboratory of Molecular Cardiogenetics (Hospices Civils de Lyon, Lyon, France) for their technical assistance.

Supplementary Materials

Supplementary Figure 1: schematic overview of our approach in atrial fibrillation (AF) patients. According to our previous transcriptomic data on left atrium (LA) in AF patients that highlighted the structural genes [24], a next-generation sequencing (NGS) approach has been applied to a cohort of 94 AF patients. Eleven potentially pathogenic variants were identified in AF patients, mainly in those with atrial dilatation, through a custom-made panel of 55 genes potentially involved in atrial myopathy. These variants were located in genes involved in atrial tissue structural remodeling. (Supplementary Materials)