BioMed Research International
Volume 2018 (2018), Article ID 4862480, 12 pages
https://doi.org/10.1155/2018/4862480

Research Article

Atrial Structural Remodeling Gene Variants in Patients with Atrial Fibrillation

Rosa Doñate Puertas,¹ Gilles Millat,1,2,3 Isabelle Ernens,⁴ Vincent Gache,¹ Samuel Chauveau,⁵ Elodie Morel,⁵ Emilie Christin,¹ Nathalie Couturier,¹ Yvan Devaux,⁴ and Philippe Chevalier1,5

¹Institut NeuroMyoGène (INMG), UMR CNRS 5310-INSERM U1217/University of Lyon, Lyon, France
²Laboratory of Molecular Cardiogenetics, Hospices Civils de Lyon, Lyon, France
³NGS Sequencing Platform for Molecular Diagnosis, Hospices Civils de Lyon, Lyon, France
⁴Cardiovascular Research Unit, Luxembourg Institute of Health, Luxembourg
⁵Rhythmology Unit, Louis Pradel Cardiology Hospital, Hospices Civils de Lyon, Lyon, France

Correspondence should be addressed to Philippe Chevalier; philippe.chevalier@chu-lyon.fr

Received 3 May 2018; Revised 30 June 2018; Accepted 17 July 2018; Published 10 September 2018

Academic Editor: Diego Franco

Copyright © 2018 Rosa Doñate Puertas et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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) [12–14] 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 [16–18].

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/m² or a surface > 22 cm², 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].

Table 1: List of the genes included in our panel.

Gene
NM number
Gene
NM number
Gene
NM number
Gene
NM number
Gene
NM number

ABCC8
NM_001287174.1
FHOD3
NM_025135.4
IL10
NM_000572.2
MYL7
NM_021223.2
SMAD2
NM_005901.5

ABCC9
NM_005691.3
FKRP
NM_001039885.2
JPH2
NM_020433.4
MYOZ1
NM_021245.3
SMAD3
NM_005902.3

ACE
NM_000789.3
GATA4
NM_001308093.1
JUP
NM_002230.2
NODAL
NM_018055.4
SMAD4
NM_005359.5

AKAP9
NM_005751.4
GATA5
NM_080473.4
LEFTY1
NM_020997.3
NOS1AP
NM_014697.2
SMYD2
NM_020197.2

CER1
NM_005454.2
GATA6
NM_005257.5
LEFTY2
NM_003240.3
NOS3
NM_000603.4
TGFB1
NM_000660.5

CTGF
NM_001901.2
GJA1
NM_000165.4
LMNA
NM_170707.3
ORAI1
NM_032790.3
TGFB3
NM_003239.3

DES
NM_001927.3
GJA5
NM_005266.6
LTB2
NM_002341.1
PCSK6
NM_002570.3
TIMP-1
NM_003254.2

DSC2
NM_024422.4
GJC1
NM_001080383.1
MMP2
NM_004530.5
PITX2
NM_153427.2
TIMP-2
NM_003255.4

DSG2
NM_001943.4
GNB3
NM_002075.3
MMP9
NM_004994.2
PKP2
NM_001005242.2
TMEM43
NM_024334.2

DSP
NM_004415.3
HSP90AB1
NM_001271971.1
MYBPC2
NM_004533.3
PRRX1
NM_006902.4
TMPO
NM_003276.2

FGFR1
NM_023110.2
HSPE1
NM_002157.2
MYH7
NM_000257.3
SHOX2
NM_003030.4
TRPM4
NM_017636.3

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).

Figure 1: Filtering steps. Decision tree for exploration of genes related to atrial fibrillation using a next-generation sequencing approach to detect mutations, based on a custom AmpliSeq library and Ion Torrent PGM sequencing. Abbreviations: MAF = minor allele frequency; PGM = Personal Genome Machine.

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.

Table 2: Clinical parameters of the patients involved in the study.

All patients (n=94)

Ratio M/F
63/31

Age (years)
54.4 ( ± 12.0)

Age of AF onset (years)
48.01 ( ± 14.35)

BMI
27.7 ( ± 5.5)

        AF type

Paroxysmal
55

Persistent
21

Permanent
18

        Risk factors

Hypertension
23

CVA
16

LA surface (cm²)
26.8 ( ± 7.5)

Abbreviations. AF = atrial fibrillation; BMI= body mass index; CVA = cerebrovascular accident; LA = left atrium.

Table 3: List of putative pathogenic variations identified in a cohort 94 patients.

Case
Patient characteristics
Variation
Presence in databases
Nucleotide conservation prediction
Grantham Score
Protein domain

AF type
Dilated LA
Gene
Chr
Exon
Nucleotide change
Effect on protein
Pathogenicity
ClinVar
ExAC
PhyloP

2115
Paroxysmal
Yes
ABCC8
11
32
c.3941G>A
p.Arg1314His
Likely pathogenic

Yes
3.98
29
ABC transporter type 1, transmembrane domain

6211
Permanent
Yes
AKAP9
7
46
c.11229G>A
p.Met3743Ile
Likely pathogenic
RCV000310251.1 (Uncertain significance), RCV000171732.1 (Likely benign), RCV000362669.1 (Uncertain significance)
Yes
7.99
10
Pericentrin/AKAP-450 centrosomal targeting domain

6198
Permanent
Yes
DSG2
18
6
c.566C>T
p.Pro189Leu
Likely pathogenic

No
1.50
98
Cadherin; Cadherin-like

2235
Paroxysmal
No
DSP
6
23
c.3550C>T
p.Arg1184Trp
Likely pathogenic

Yes
3.98
101

4464
Permanent
Yes
DSP
6
24
c.7997G>A
p.Gly2666Asp
Likely pathogenic

Yes
10.00
94
Plectin repeat

2095
Paroxysmal
Yes
FHOD3
18
7
c.614T>C
p.Leu205Pro
Likely pathogenic

No
7.99
98
Rho GTPase-binding/formin homology 3 (GBD/FH3) domain; Armadillo-type fold

1885
Paroxysmal
Yes
FHOD3
18
8
c.776C>T
p.Thr259Met
Likely pathogenic

Yes
7.76
81
Rho GTPase-binding/formin homology 3 (GBD/FH3) domain; Armadillo-type fold

4162
Permanent
Yes
JPH2
20
2
c.764C>T
p.Ser255Leu
Likely pathogenic

No
9.18
145
Junctophilin

2186
Paroxysmal
Yes
MMP9
20
3
c.427C>T
p.Arg143Cys
Likely pathogenic

Yes
2.68
180
Peptidase M10, metallopeptidase; Peptidase, metallopeptidase; Peptidase M10A

1875
Paroxysmal
No
MYOZ1
10
3
c.167G>C
p.Gly56Ala
Likely pathogenic

Yes
7.76
60
Myozenin

2691
Paroxysmal
Yes
TMEM43
3
5
c.424G>A
p.Glu142Lys
Likely pathogenic
RCV000039386.3 (Uncertain significance), RCV000172593.3 (Likely benign), RCV000250239.1 (Uncertain significance)
Yes
8.11
56
Transmembrane protein 43 family

LA volume >32 ml/m² or a surface >22 cm²; specific standard terminologies—“pathogenic”, “likely pathogenic”, “uncertain significance”, “likely benign”, and “benign” were used to describe variants identified (Ref [28]); range of PhyloP score -20.0;10.0; range of Grantham score 0-215. Abbreviations: AF = atrial fibrillation; ExAC = Exome Aggregation Consortium; LA= left atrium.

Figure 2: High conservation across species. Multiple protein sequence alignments and the evolutionary conservation of each altered amino acid among species (H. sapiens, P. troglodytes, M. mulatta, C. lupus, and B. taurus).

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.

Table 4: List of variants identified in the prevalent arrhythmia-causing genes.

Case
Variation
Presence in databases
Prediction analysis of missense variants
Nucleotide conservation prediction
Grantham Score

Gene
NM_number
Chr
Exon/ intron
Nucleotide change
Effect on protein
Pathogenicity
ClinVar
ExAC
PolyPhen-2
SIFT
Mutation Taster
PhyloP

2095
ANK2
NM_001148.4
4
35
c.4315G>A
p.Gly1439Ser
Likely pathogenic

No
Possibly damaging
Tolerated
Disease causing
2.80
56

6198
KCNH2
NM_000238.3
7
11
c.2681G>A
p.Arg894His
Likely pathogenic

Yes (MAF ≤0.01%)
Probably damaging
Tolerated
Disease causing
2.68
29

1885
KCNH2
NM_000238.4
7
13
c.3052C>G
p.Pro1018Ala
Uncertain significance
RCV000181908.1 (Uncertain significance)
Yes (MAF ≤0.01%)
Benign
Tolerated
Disease causing
1.50
27

1885
SCN1B
NM_001037.4
19
intron 1
c.40+2T>G

Likely pathogenic

No

Specific standard terminologies—“pathogenic”, “likely pathogenic”, “uncertain significance”, “likely benign” and “benign” were used to describe variants identified [28]; range of PhyloP score -20.0;10.0; range of Grantham score 0-215. Abbreviations: ExAC = Exome Aggregation Consortium; MAF = minor allele frequency.

Table 5: Genes associated with cardiac diseases.

Gene
Protein
Cardiac diseases classification
Protein localisation

LQTS
SQTS
BrS
SIDS
CDM
CHD
AF

ABCC8
Sulfonylurea receptor 1
X
–
–
–
X
X
X
Sarcolemma

AKAP9
A-kinase anchor protein 9
X
–
X
X
X
–
–
Centrosome

DSG2
Desmoglein 2
–
–
X
–
X
–
X
Desmosome

DSP
Desmoplakin
X
–
X
X
X
–
X
Desmosome

FHOD3
FH1/FH2 domain-containing protein 3
–
–
–
–
X
–
–
Z-disk

JPH2
Junctophilin 2
–
–
X
–
X
–
X
Sarcoplasmic reticulum

MMP9
Matrix metalloproteinase-9
X
–
–
–
X
X
X
Extracellular matrix

MYOZ1
Myozenin-1
–
–
–
–
X
–
X
Z-disk

TMEM43
Transmembrane protein 43
–
–
X
–
X
–
–
Transmembrane

The presence of an (X) indicates involvement of the gene in each specific cardiac disease classification. The localisation of the protein encoded by each disease gene is also described. Abbreviations: BrS = Brugada syndrome; CDM = cardiomyopathies; CHD = congenital heart diseases; LQTS = long QT syndrome; SIDS = sudden infant death syndrome; SQTS = short QT syndrome.

Figure 3: Atrial fibrillation disease genes. A schematic of proteins encoded by genes related to atrial fibrillation and their subcellular localization. Proteins participate in many diverse biological processes of cardiomyocytes/fibroblasts.

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].

(a)

(b)

(c)

Figure 4: AKAP9 is required for myonuclear positioning in C2C12 myotubes. (a) Representative immunofluorescence images of control (no transfection and Scramble siRNA treated cells) or AKAP9-depleted C2C12 myotubes (using 3 individual siRNA, 30 nM each) differentiated for 5 days and immunostained for myosin heavy chain (green) and 49,6-diamidino-2-phenylindole (red). Scale bar, 160 um. (b) 7× magnifications of rectangles shown in images (a). Scale bar, 160 um. (c) Histogram of percentage of total C2C12 myotubes with aggregated nuclei control (no transfection and Scramble siRNA treated cells) or AKAP9-depleted C2C12 myotubes (using 3 individual siRNA, 30 nM each) differentiated for 5 days. Center lines show the medians; box limits indicate the 25th and 75th percentiles as determined by R software; whiskers extend to 5th and 95th percentiles, outliers are represented by dots; width of the boxes is proportional to the square root of the sample size; data points are plotted as open circles. n = 6, 12, 9, 9, 7 sample points. Student’s t-tests were performed between scrambled siRNA and experimental condition. Asterisk, P, 0.05; two asterisks, P, 0.01; ns: nonsignificant.

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)

References

S. S. Chugh, R. Havmoeller, K. Narayanan et al., “Worldwide epidemiology of atrial fibrillation: a global burden of disease 2010 study,” Circulation, vol. 129, no. 8, pp. 837–847, 2014.

P. Kirchhof, G. Y. H. Lip, I. C. Van Gelder et al., “Comprehensive risk reduction in patients with atrial fibrillation: Emerging diagnostic and therapeutic options—a report from the 3rd Atrial Fibrillation Competence NETwork/European Heart Rhythm Association consensus conference,” Europace, vol. 14, no. 1, pp. 8–27, 2012.

S. Nattel, B. Burstein, and D. Dobrev, “Atrial remodeling and atrial fibrillation: mechanisms and implications.,” Circulation: Arrhythmia and Electrophysiology, vol. 1, no. 1, pp. 62–73, 2008.

P.-S. Chen, L. S. Chen, M. C. Fishbein, S.-F. Lin, and S. Nattel, “Role of the autonomic nervous system in atrial fibrillation: pathophysiology and therapy,” Circulation Research, vol. 114, no. 9, pp. 1500–1515, 2014.

C. S. Fox, H. Parise, R. B. D'Agostino et al., “Parental atrial fibrillation as a risk factor for atrial fibrillation in offspring,” Journal of the American Medical Association, vol. 291, no. 23, pp. 2851–2855, 2004.

Y.-H. Chen, S.-J. Xu, S. Bendahhou et al., “KCNQ1 gain-of-function mutation in familial atrial fibrillation,” Science, vol. 299, no. 5604, pp. 251–254, 2003.

D. M. Hodgson-Zingman, M. L. Karst, L. V. Zingman et al., “Atrial natriuretic peptide frameshift mutation in familial atrial fibrillation,” The New England Journal of Medicine, vol. 359, no. 2, pp. 158–165, 2008.

P. G. A. Volders, Q. Zhu, C. Timmermans et al., “Mapping a novel locus for familial atrial fibrillation on chromosome 10p11-q21,” Heart Rhythm, vol. 4, no. 4, pp. 469–475, 2007.

P. T. Ellinor, J. T. Shin, R. K. Moore, D. M. Yoerger, and C. A. MacRae, “Locus for atrial fibrillation maps to chromosome 6q14-16,” Circulation, vol. 107, no. 23, pp. 2880–2883, 2003.

I. L. Thibodeau, J. Xu, Q. Li et al., “Paradigm of genetic mosaicism and lone atrial fibrillation: Physiological characterization of a connexin 43-deletion mutant identified from atrial tissue,” Circulation, vol. 122, no. 3, pp. 236–244, 2010.

J. Jabbari, M. S. Olesen, L. Yuan et al., “Common and rare variants in SCN10A modulate the risk of atrial fibrillation,” Circulation: Cardiovascular Genetics, vol. 8, no. 1, pp. 64–73, 2015.

D. F. Gudbjartsson, D. O. Arnar, A. Helgadottir et al., “Variants conferring risk of atrial fibrillation on chromosome 4q25,” Nature, vol. 448, no. 7151, pp. 353–357, 2007.

P. T. Ellinor, K. L. Lunetta, N. L. Glazer et al., “Common variants in KCNN3 are associated with lone atrial fibrillation,” Nature Genetics, vol. 42, no. 3, pp. 240–244, 2010.

P. T. Ellinor, K. L. Lunetta, C. M. Albert et al., “Meta-analysis identifies six new susceptibility loci for atrial fibrillation,” Nature Genetics, vol. 44, no. 6, pp. 670–675, 2012.

B. Sikkema-Raddatz, L. F. Johansson, E. N. de Boer et al., “Targeted next-generation sequencing can replace sanger sequencing in clinical diagnostics,” Human Mutation, vol. 34, no. 7, pp. 1035–1042, 2013.

C.-T. Tsai, C.-S. Hsieh, S.-N. Chang et al., “Next-generation sequencing of nine atrial fibrillation candidate genes identified novel de novo mutations in patients with extreme trait of atrial fibrillation,” Journal of Medical Genetics, vol. 52, no. 1, pp. 28–36, 2015.

H. Lin, M. F. Sinner, J. A. Brody et al., “Targeted sequencing in candidate genes for atrial fibrillation: The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Targeted Sequencing Study,” Heart Rhythm, vol. 11, no. 3, pp. 452–457, 2014.

P. Weeke, B. Parvez, M. Blair et al., “Candidate gene approach to identifying rare genetic variants associated with lone atrial fibrillation,” Heart Rhythm, vol. 11, no. 1, pp. 46–52, 2014.

A. Boldt, U. Wetzel, J. Lauschke et al., “Fibrosis in left atrial tissue of patients with atrial fibrillation with and without underlying mitral valve disease,” Heart, vol. 90, no. 4, pp. 400–405, 2004.

A. Frustaci, C. Chimenti, F. Bellocci, E. Morgante, M. A. Russo, and A. Maseri, “Histological substrate of atrial biopsies in patients with lone atrial fibrillation,” Circulation, vol. 96, no. 4, pp. 1180–1184, 1997.

Q.-Q. Cui, W. Zhang, H. Wang et al., “Assessment of atrial electromechanical coupling and influential factors in nonrheumatic paroxysmal atrial fibrillation,” Clinical Cardiology, vol. 31, no. 2, pp. 74–78, 2008.

D. F. Gudbjartsson, H. Holm, P. Sulem et al., “A frameshift deletion in the sarcomere gene MYL4 causes early-onset familial atrial fibrillation,” European Heart Journal, vol. 38, no. 1, pp. 27–34, 2017.

P. Kirchhof, S. Benussi, D. Kotecha et al., “2016 ESC Guidelines for the management of atrial fibrillation developed in collaboration with EACTS,” Europace, vol. 18, no. 11, pp. 1609–1678, 2016.

R. Doñate Puertas, E. Meugnier, C. Romestaing et al., “Atrial fibrillation is associated with hypermethylation in human left atrium, and treatment with decitabine reduces atrial tachyarrhythmias in spontaneously hypertensive rats,” Translational Research, vol. 184, pp. 57–67.e5, 2017.

G. Millat, V. Chanavat, and R. Rousson, “Evaluation of a new NGS method based on a custom AmpliSeq library and Ion Torrent PGM sequencing for the fast detection of genetic variations in cardiomyopathies,” Clinica Chimica Acta, vol. 433, pp. 266–271, 2014.

G. Millat, V. Chanavat, and R. Rousson, “Evaluation of a New High-Throughput Next-Generation Sequencing Method Based on a Custom AmpliSeq™ Library and Ion Torrent PGM™ Sequencing for the Rapid Detection of Genetic Variations in Long QT Syndrome,” Molecular Diagnosis & Therapy, vol. 18, no. 5, pp. 533–539, 2014.

V. Chanavat, A. Janin, and G. Millat, “A fast and cost-effective molecular diagnostic tool for genetic diseases involved in sudden cardiac death,” Clinica Chimica Acta, vol. 453, pp. 80–85, 2016.

S. Richards, N. Aziz, and S. Bale, “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,” Genetics in Medicine, vol. 17, no. 5, pp. 405–423, 2015.

I. A. Adzhubei, S. Schmidt, L. Peshkin et al., “A method and server for predicting damaging missense mutations,” Nature Methods, vol. 7, no. 4, pp. 248-249, 2010.

P. Kumar, S. Henikoff, and P. C. Ng, “Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm,” Nature Protocols, vol. 4, no. 7, pp. 1073–1082, 2009.

J. M. Schwarz, C. Rödelsperger, M. Schuelke, and D. Seelow, “MutationTaster evaluates disease-causing potential of sequence alterations,” Nature Methods, vol. 7, no. 8, pp. 575-576, 2010.

R. Grantham, “Amino acid difference formula to help explain protein evolution,” Science, vol. 185, no. 4154, pp. 862–864, 1974.

R. C. Edgar, “MUSCLE: multiple sequence alignment with high accuracy and high throughput,” Nucleic Acids Research, vol. 32, no. 5, pp. 1792–1797, 2004.

B. Cadot, V. Gache, E. Vasyutina, S. Falcone, C. Birchmeier, and E. R. Gomes, “Nuclear movement during myotube formation is microtubule and dynein dependent and is regulated by Cdc42, Par6 and Par3,” EMBO Reports, vol. 13, no. 8, pp. 741–749, 2012.

S. Rivero, J. Cardenas, M. Bornens, and R. M. Rios, “Microtubule nucleation at the cis-side of the golgi apparatus requires AKAP450 and GM130,” EMBO Journal, vol. 28, no. 8, pp. 1016–1028, 2009.

P. Gimpel, Y. L. Lee, R. M. Sobota et al., “Nesprin-1α-Dependent Microtubule Nucleation from the Nuclear Envelope via Akap450 Is Necessary for Nuclear Positioning in Muscle Cells,” Current Biology, vol. 27, no. 19, pp. 2999–3009.e9, 2017.

N. R. Tucker, S. Clauss, and P. T. Ellinor, “Common variation in atrial fibrillation: Navigating the path from genetic association to mechanism,” Cardiovascular Research, vol. 109, no. 4, pp. 493–501, 2016.

S. R. Cunha, T. J. Hund, S. Hashemi et al., “Defects in ankyrin-based membrane protein targeting pathways underlie atrial fibrillation,” Circulation, vol. 124, no. 11, pp. 1212–1222, 2011.

K. Hayashi, T. Konno, H. Tada et al., “Functional Characterization of Rare Variants Implicated in Susceptibility to Lone Atrial Fibrillation,” Circulation: Arrhythmia and Electrophysiology, vol. 8, no. 5, pp. 1095–1104, 2015.

W. L. Henry, J. Morganroth, A. S. Pearlman et al., “Relation between echocardiographically determined left atrial size and atrial fibrillation,” Circulation, vol. 53, no. 2, pp. 273–279, 1976.

C.-L. Chen, S. K. S. Huang, J.-L. Lin et al., “Upregulation of matrix metalloproteinase-9 and tissue inhibitors of metalloproteinases in rapid atrial pacing-induced atrial fibrillation,” Journal of Molecular and Cellular Cardiology, vol. 45, no. 6, pp. 742–753, 2008.

M. Li, G. Yang, B. Xie, K. Babu, and C. Huang, “Changes in matrix metalloproteinase-9 levels during progression of atrial fibrillation,” Journal of International Medical Research, vol. 42, no. 1, pp. 224–230, 2014.

R. D. Nadadur, M. T. Broman, B. Boukens et al., “Pitx2 modulates a Tbx5-dependent gene regulatory network to maintain atrial rhythm,” Science Translational Medicine, vol. 8, no. 354, 2016.

M. Kan-o, R. Takeya, K. Taniguchi, Y. Tanoue, R. Tominaga, and H. Sumimoto, “Expression and subcellular localization of mammalian formin Fhod3 in the embryonic and adult heart,” PLoS ONE, vol. 7, no. 4, 2012.

N. Frey, J. A. Richardson, and E. N. Olson, “Calsarcins, a novel family of sarcomeric calcineurin-binding proteins,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 97, no. 26, pp. 14632–14637, 2000.

A. Garbino and X. H. Wehrens, “Emerging role of junctophilin-2 as a regulator of calcium handling in the heart,” Acta Pharmacologica Sinica, vol. 31, no. 9, pp. 1019–1021, 2010.

D. L. Beavers, W. Wang, S. Ather et al., “Mutation E169K in junctophilin-2 causes atrial fibrillation due to impaired RyR2 stabilization,” Journal of the American College of Cardiology, vol. 62, no. 21, pp. 2010–2019, 2013.

L. Bengtsson and H. Otto, “LUMA interacts with emerin and influences its distribution at the inner nuclear membrane,” Journal of Cell Science, vol. 121, no. 4, pp. 536–548, 2008.

L. Chen, M. L. Marquardt, D. J. Tester, K. J. Sampson, M. J. Ackerman, and R. S. Kass, “Mutation of an A-kinase-anchoring protein causes long-QT syndrome,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 104, no. 52, pp. 20990–20995, 2007.

V. Gache, E. R. Gomes, and B. Cadot, “Microtubule motors involved in nuclear movement during skeletal muscle differentiation,” Molecular Biology of the Cell (MBoC), vol. 28, no. 7, pp. 865–874, 2017.

T. Metzger, V. Gache, M. Xu et al., “MAP and kinesin-dependent nuclear positioning is required for skeletal muscle function,” Nature, vol. 484, no. 7392, pp. 120–124, 2012.

P. Robison, M. A. Caporizzo, H. Ahmadzadeh et al., “Detyrosinated microtubules buckle and bear load in contracting cardiomyocytes,” Science, vol. 352, no. 6284, Article ID aaf0659, 2016.

J. Ye, N. R. Tucker, L.-C. Weng, S. Clauss, S. A. Lubitz, and P. T. Ellinor, “A Functional Variant Associated with Atrial Fibrillation Regulates PITX2c Expression through TFAP2a,” American Journal of Human Genetics, vol. 99, no. 6, pp. 1281–1291, 2016.

R. Doñate Puertas, A. Jalabert, E. Meugnier et al., “Analysis of the microRNA signature in left atrium from patients with valvular heart disease reveals their implications in atrial fibrillation,” PLoS ONE, vol. 13, no. 5, p. e0196666, 2018.

S. Hoffmann, S. Clauss, I. M. Berger et al., “Coding and non-coding variants in the SHOX2 gene in patients with early-onset atrial fibrillation,” Basic Research in Cardiology, vol. 111, no. 3, article 36, pp. 1–15, 2016.

B. Parvez, J. Vaglio, S. Rowan et al., “Symptomatic response to antiarrhythmic drug therapy is modulated by a common single nucleotide polymorphism in atrial fibrillation,” Journal of the American College of Cardiology, vol. 60, no. 6, pp. 539–545, 2012.

E. Piñeiro-Yáñez, M. Reboiro-Jato, G. Gómez-López et al., “PanDrugs: a novel method to prioritize anticancer drug treatments according to individual genomic data,” Genome Medicine, vol. 10, no. 1, 2018.