Open Access

Identification of TNFSF13, SPATC1L, SLC22A25 and SALL4 as novel susceptibility loci for atrial fibrillation by an exome‑wide association study

  • Authors:
    • Yoshiji Yamada
    • Jun Sakuma
    • Ichiro Takeuchi
    • Yoshiki Yasukochi
    • Kimihiko Kato
    • Mitsutoshi Oguri
    • Tetsuo Fujimaki
    • Hideki Horibe
    • Masaaki Muramatsu
    • Motoji Sawabe
    • Yoshinori Fujiwara
    • Yu Taniguchi
    • Shuichi Obuchi
    • Hisashi Kawai
    • Shoji Shinkai
    • Seijiro Mori
    • Tomio Arai
    • Masashi Tanaka
  • View Affiliations

  • Published online on: August 23, 2017     https://doi.org/10.3892/mmr.2017.7334
  • Pages: 5823-5832
  • Copyright: © Yamada et al. This is an open access article distributed under the terms of Creative Commons Attribution License.

Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

An exome‑wide association study (EWAS) was performed to identify genetic variants, particularly low‑frequency or rare coding variants with a moderate to large effect size, that confer susceptibility to atrial fibrillation in Japanese. The EWAS for atrial fibrillation was performed with 13,166 subjects (884 patients with atrial fibrillation and 12,282 controls) using an Illumina HumanExome‑12 DNA Analysis BeadChip or Infinium Exome‑24 BeadChip arrays. The association of atrial fibrillation with allele frequencies of 41,243 single nucleotide polymorphisms (SNPs) that passed quality control was examined with Fisher's exact test. Based on Bonferroni's correction, a P<1.21x10‑6 was considered statistically significant. The EWAS for atrial fibrillation revealed that 122 SNPs were significantly associated with this condition. The association of the identified SNPs to atrial fibrillation was further examined by multivariable logistic regression analysis with adjustment for age, sex and the prevalence of hypertension. Eight SNPs were related (P<0.01) to atrial fibrillation, among which three polymorphisms, rs11552708 [G/A (G67R)]of TNF superfamily member 13 (TNFSF13; dominant model; P=9.36x10‑9; odds ratio, 0.58), rs113710653 [C/T (E231 K)] of spermatogenesis and centriole associated 1 like (SPATC1L; dominant model; P=1.09x10‑5; odds ratio, 3.27), and rs11231397 [G/C (R300T)] of solute carrier family 22 member 25 (SLC22A25; additive model; P=3.71x10‑5; odds ratio, 1.77), were significantly (P<1.02x10‑4) associated with this condition. The minor T allele of rs113710653 and the minor C allele of rs11231397 were risk factors for atrial fibrillation, whereas the minor A allele of rs11552708 was protective against this condition. In addition, rs77538589 [C/T (G117R)] of SALL4 exhibited a tendency to be associated with atrial fibrillation (dominant model; P=0.0002; odds ratio, 1.88), with the minor T allele representing a risk factor for this condition. TNFSF13, SPATC1L, SLC22A25 and SALL4 may thus be novel susceptibility loci for atrial fibrillation in the Japanese population.

Introduction

Atrial fibrillation (AF) is the most common type of cardiac arrhythmia and a major public health problem. The estimated number of individuals with AF worldwide was 33.5 million in 2010 (1). The prevalence of AF is increasing and is estimated to double in the United States by 2050 (2). AF is associated with an increased risk of heart failure, thromboembolic diseases, such as cardioembolic stroke, and mortality (3). Although the molecular mechanism underlying AF is complex and has not been determined definitively, several risk factors, including age, smoking, obesity, hypertension, diabetes mellitus, valvular heart disease, coronary artery disease and heart failure, have been identified clinically (4). In addition to these conventional risk factors, recent studies have reported the importance of genetic factors in the development of AF (5). Genome-wide association studies (GWASs) and subsequent meta-analyses of such studies have identified several genes and loci that confer susceptibility to AF in populations of European ancestry (610). However, susceptibility genes for AF have not been determined definitively in Japanese populations.

The majority of genetic variants previously reported to be associated with AF have a minor allele frequency (MAF) of ≥5% and a small individual effect size. Given that these common variants explain only a fraction of the heritability of AF, it is likely that low-frequency (0.5%≤ MAF <5%) or rare (MAF <0.5%) variants with larger effect sizes contribute to the genetic architecture of AF (11).

An exome-wide association study (EWAS) with exome array-based genotyping was performed to identify single nucleotide polymorphisms (SNPs), especially low-frequency or rare coding variants with moderate to large effect sizes, that confer susceptibility to AF in Japanese individuals. Given that GWAS arrays of previous studies did not include most low-frequency or rare variants, Illumina human arrays that cover functional SNPs in entire exons including such variants were used in the current study.

Materials and methods

Study subjects

A total of 13,166 Japanese individuals (884 patients with AF and 12,282 controls) were examined. The subjects were recruited either from individuals who visited outpatient clinics of or were admitted to participating hospitals (Gifu Prefectural Tajimi Hospital, Tajimi; Gifu Prefectural General Medical Center, Gifu; Japanese Red Cross Nagoya First Hospital, Nagoya; Inabe General Hospital, Inabe; Hirosaki University Hospital and Hirosaki Stroke and Rehabilitation Center, Hirosaki, Japan) because of various symptoms or for an annual health checkup between 2002 and 2014; from community-dwelling individuals recruited to a population-based cohort study in Inabe between 2010 and 2014 (Health Care Center of Inabe General Hospital) or in Tokyo or Kusatsu between 2011 and 2015 (Tokyo Metropolitan Institute of Gerontology); or from cases of autopsy performed at the Tokyo Metropolitan Geriatric Hospital from 1995 to 2012.

Individuals with persistent or paroxysmal AF or with a history of AF were included in the AF patient group regardless of medical treatment. Patients with AF who had severe valvular heart disease, hypertrophic or dilated cardiomyopathy, or congenital heart disease were excluded from the study. The control individuals had no history of AF or other significant supraventricular or ventricular arrhythmias or of taking antiarrhythmic medication. Autopsy cases were not used as controls.

The study protocol complied with the Declaration of Helsinki and was approved by the Committees on the Ethics of Human Research of Mie University Graduate School of Medicine (Tsu, Japan), Hirosaki University Graduate School of Medicine (Hirosaki, Japan), Tokyo Metropolitan Institute of Gerontology (Tokyo, Japan) and participating hospitals. Written informed consent was obtained from each participant or families of the deceased subjects.

EWAS

Venous blood (5–7 ml) was collected into tubes containing 50 mmol/l ethylenediaminetetraacetic acid (disodium salt), peripheral blood leukocytes were isolated, and genomic DNA was extracted from the cells with the use of a DNA extraction kit (Genomix supplied by Talent Srl, Trieste, Italy; or SMITEST EX-R&D supplied by Medical & Biological Laboratories, Co., Ltd., Nagoya, Japan) or by standard protocols based on phenol-chloroform extraction and spin columns. Lysed blood cells were mixed with an equal volume of a phenol-chloroform mixture. Following centrifugation at 16,000 × g for 5 min at room temperature, two distinct phases were formed. The upper aqueous phase was collected, and genomic DNA was further purified by a spin column-based DNA extraction method that is based on the fact that DNA binds to the solid phase of silica under certain conditions. A buffer solution (Buffer B3; Takara Bio, Inc., Otsu, Japan) was added to the DNA sample along with 100% ethanol. This solution was transferred to a spin column, and the column was centrifuged at 11,000 × g for 1 min at room temperature. A wash buffer [Buffer BW (first wash) or Buffer B5 (second wash); Takara Bio, Inc.] was added to the column, and the column was centrifuged at 11,000 × g for 1 min at room temperature again. An elution buffer (Buffer BE; Takara Bio, Inc.) was added to the column, and the column was centrifuged at 11,000 × g for 1 min at room temperature again. Genomic DNA was collected from the bottom of the column. In autopsy cases, genomic DNA was extracted from kidneys. The EWAS was performed for 884 subjects with AF and 12,282 control individuals with the use of a HumanExome-12 v1.1 or v1.2 DNA Analysis BeadChip or Infinium Exome-24 v1.0 BeadChip (Illumina, Inc., San Diego, CA, USA). These exome arrays include putative functional exonic variants selected from >12,000 individual exome or whole-genome sequences. The exonic content consists of ~244,000 SNPs representing diverse populations, including European, African, Chinese and Hispanic individuals (12). SNPs contained in only one of the exome arrays (~3.6%) were excluded from analysis. Quality control (13) was performed as follows: i) Genotyping data with a call rate of <97% were discarded, with the mean call rate for the remaining data at 99.9%; ii) gender specification was checked for all samples, and those for which gender in the clinical records was inconsistent with genetic sex were discarded; iii) duplicate samples and cryptic relatedness were checked by calculation of identity by descent, all pairs of DNA samples with identity by descent of >0.1875 were inspected and one sample from each pair was excluded; iv) the frequency of SNP heterozygosity was calculated for all samples, and those with extremely low or high heterozygosity (>3 standard deviations from the mean) were discarded; v) SNPs in sex chromosomes or mitochondrial DNA were excluded from the analysis, as were nonpolymorphic SNPs or SNPs with a MAF of <0.1%; vi) SNPs whose genotype distributions deviated significantly (P<0.001) from Hardy-Weinberg equilibrium in control individuals were discarded; vii) the genotype data for the EWAS were examined for population stratification by principal components analysis (14), and population outliers were excluded from the analysis. A total of 41,243 SNPs passed quality control and was subjected to analysis.

Statistical analysis

For analysis of characteristics of the study subjects, quantitative data were compared between patients with AF and control individuals with the unpaired Student's t test. Categorical data were compared between the two groups with Fisher's exact test. Allele frequencies were estimated by the gene counting method, and Fisher's exact test was applied to identify departure from Hardy-Weinberg equilibrium. Allele frequencies of SNPs were compared between subjects with AF and controls with Fisher's exact test in the EWAS. Multivariable logistic regression analysis was performed with AF as a dependent variable and independent variables including age, sex (0, woman; 1, man), the prevalence of hypertension (0, no history of this condition; 1, positive history), and genotype of each SNP. Genotypes of SNPs were assessed according to dominant [0, AA; 1, AB+BB (A, major allele; B, minor allele)], recessive (0, AA+AB; 1, BB), and additive genetic models, and the P-value, odds ratio, and 95% confidence interval were calculated. Additive models comprised additive 1 (0, AA; 1, AB; 0, BB) and additive 2 (0, AA; 0, AB; 1, BB) scenarios, which were analyzed simultaneously with a single statistical model. The relation of SNPs to the prevalence of intermediate phenotypes of AF was compared among genotypes with Fisher's exact test (2×2) or Pearson's chi-square test (2×3). To compensate for multiple comparisons of genotypes with AF, we applied Bonferroni's correction for statistical significance of association. Given that 41,243 SNPs were analyzed, the significance level was set at P<1.21×10−6 (0.05/41,243) for the EWAS. A quantile-quantile plot for P-values of allele frequencies in the EWAS for AF is shown in Fig. 1. The inflation factor (λ) was 1.52. Bonferroni's correction was also applied to other statistical comparisons as indicated. Statistical tests were performed with JMP Genomics version 6.0 software (SAS Institute, Inc., Cary, NC, USA).

Results

Characteristics of the study subjects

The characteristics of the subjects enrolled in the study are presented in Table I. Age, the number of men, and the prevalence of hypertension, diabetes mellitus, chronic kidney disease and hyperuricemia were significantly greater in patients with AF than in control subjects.

Table I.

Characteristics of the 13,166 study subjects.

Table I.

Characteristics of the 13,166 study subjects.

CharacteristicAtrial fibrillationControlP-value
No. of subjects88412,282
Age (years)75.8±12.359.4±13.2<0.0001
Sex (male/female, %)67.8/32.256.8/43.2<0.0001
Body mass index (kg/m2)23.4±3.723.3±3.50.7412
Current or former smoker (%)33.636.30.2576
Hypertension (%)83.951.1<0.0001
Diabetes mellitus (%)5123.2<0.0001
Dyslipidemia (%)65.3620.188
Chronic kidney disease (%)40.323<0.0001
Hyperuricemia (%)31.917<0.0001

[i] Quantitative data are presented as the mean ± standard deviation and were compared between patients with atrial fibrillation and control individuals using the unpaired Student's t test. Categorical data were compared between the two groups using Fisher's exact test. Based on Bonferroni's correction, P<0.0056 (0.05/9) was considered to indicate a statistically significant difference.

EWAS of AF

The association of AF with allele frequencies of 41,243 SNPs that passed quality control was determined using Fisher's exact test. A Manhattan plot for the EWAS of AF is presented in Fig. 2. Following Bonferroni's correction, 122 SNPs were significantly [P<1.21×10−6 (0.05/41,243)] associated with AF (Table II). The genotype distributions of these SNPs were in Hardy-Weinberg equilibrium (P>0.001) both among subjects with AF and among control individuals (data not shown).

Table II.

The 122 SNPs significantly associated with atrial fibrillation by exome-wide association analysis.

Table II.

The 122 SNPs significantly associated with atrial fibrillation by exome-wide association analysis.

GeneSNPNucleotide substitution (amino acid)aChromosome: PositionMAF (%)P-value (allele)Allele OR
N4BP2rs61748749T/G (S1353R)4: 401221702.3 <1.0×10−2550.74
DNAH17rs690844A/C (I1742M)17: 785018382.5 8.94×10−2551.31
HELZrs184499441C/T (G1288R)17: 671143801.6 4.40×10−2401.08
rs7828656A/C8: 1164542544.8 1.87×10−2201.1
SLA2rs221308T/C20: 3664799546.7 4.11×10−2180.99
SSPOrs55976638G/T7: 1498328814.3 2.99×10−1821.31
TCEB3Brs2010834A/C (F254C)18: 4703450424.5 2.60×10−1411.1
FANCLrs149731356T/C (T224A)2: 581657450.7 7.35×10−1300.75
PIEZO1rs143004911G/A (R333C)16: 887379572.7 5.07×10−1220.87
TTYH2rs9899862C/A (D423E)17: 742530906.3 1.31×10−960.86
TUBB3rs2302898A/G16: 8993238623.5 1.13×10−870.87
SLCO6A1rs17150488T/C (K381R)5: 1024387510.6 3.30×10−811.05
GMDSrs9378305C/T6: 170305641.4 1.20×10−540.99
RP1L1rs79329877T/C8: 106227019.2 2.05×10−531.05
FCMRrs150080259T/G (S61R)1: 2069029761.5 1.07×10−491.02
RTKN2rs7090884A/G10: 6220226733.2 4.97×10−461.05
UTP4rs193164904A/G (I534V)16: 691631310.2 1.90×10−451.58
SNAPC1rs74810099T/G (M36R)14: 617625672.8 1.99×10−431.02
ALPK2rs3809977G/T (P1174H)18: 5853666617.2 3.05×10−420.93
CSPG4rs137981794T/C (D1936G)15: 756767120.4 4.05×10−42ND
MDN1rs9294445A/G (Y3423H)6: 896927637.5 1.15×10−410.99
SETD7rs6814310C/A4: 13953437448.4 1.64×10−410.97
PLA2G4Ers4924595T/C (N400S)15: 4199540818.2 1.87×10−351.01
KIF7rs117123311C/G (S788R)15: 896422333.1 1.03×10−321
CTBP2rs3781411C/T (R298Q)10: 1250268677 1.46×10−300.9
GATA2rs78245253G/C (A250P)3: 1284858504.6 1.37×10−270.79
DLGAP1rs3745051C/T18: 372917534.4 1.64×10−260.86
rs1711393T/C4: 17807408831 1.70×10−261.14
SLAMF7rs117009784A/C (R96S)1: 1607500535.8 3.86×10−261.07
USP32rs8079220C/T17: 6023937220.5 3.56×10−250.91
rs8011192T/G14: 2831941233.1 8.46×10−250.99
IMPDH1rs201001000G/A (T369M)7: 1283950030.1 3.85×10−240.79
ADRA1Ars151273238G/A (T391M)8: 267703780.2 2.12×10−230.71
TNFSF13rs11552708G/A (G67R)17: 755923836.6 1.32×10−220.54
SLC18A3rs118107581A/G (I426V)10: 496120162 1.62×10−221.11
NFATC2rs12479626T/C (H446R)20: 514756564.7 2.79×10−221.06
TENM4rs3812723C/T (V396I)11: 788630311.7 1.03×10−210.89
EPN1rs200478642C/T (P203L)19: 556893010.1 1.52×10−210.35
HNRNPCrs17197037A/G14: 2125749529.1 3.26×10−210.99
TMX4rs2076015T/C (R303G)20: 798239441.1 4.49×10−211.16
FOXN4rs140167217G/A (S308F)12: 1092817780.4 4.51×10−191.12
CEP152rs145138194G/A (S894F)15: 487601481.2 6.92×10−180.73
FREM2rs114864077C/T (P128L)13: 386877270.1 1.04×10−171.81
CPA6rs4737845T/C8: 6754270743.3 1.87×10−171.04
KIF15rs146292440G/A (R1199H)3: 448431350.6 3.33×10−170.57
MFSD1rs3765083A/G (I230V)3: 15881965430.3 4.19×10−171.05
BAHD1rs3743143A/G (E26G)15: 4045854110.5 5.85×10−170.99
rs1395821A/G4: 14712639839.6 6.48×10−171.05
BRWD1rs2183573G/A (P1511S)21: 3920237946 1.11×10−161.06
CD69rs199676648G/A (R32C)12: 97563900.3 2.93×10−160.9
HRrs12675375C/T (G337D)8: 2212743243.4 8.47×10−161.02
SOAT1rs143616084G/A (R292Q)1: 1793450080.1 6.19×10−15ND
JMJD1Crs149833441T/C (K878E)10: 632084911.2 9.49×10−151.44
VWDErs848016A/G (F142S)7: 1238917736.7 1.41×10−141.05
VPS13Drs143833298G/A (R830Q)1: 122760770.8 1.59×10−140.86
SPATC1Lrs113710653C/T (E231K)21: 461619211.9 1.60×10−145.36
SNX19rs117834100C/A (G416C)11: 1309146946.9 1.83×10−140.89
rs9854207A/C3: 2757282540.3 3.41×10−141.04
ARHGAP8rs5766113A/G22: 4485554336.8 7.18×10−141.03
rs4407763G/A7: 6755370535.8 1.14×10−131.08
SLC22A25rs11231397G/C (R300T)11: 6318374943.7 1.53×10−131.55
XIRP2rs77219745G/A (G1839D)2: 1672469087 2.89×10−130.81
MCM10rs7905784A/T (T541S)10: 131923562.2 5.31×10−131.05
HIST1H2ACrs198823G/T6: 2612270522.9 1.60×10−120.86
rs10102598G/A8: 4717812837.7 1.61×10−120.96
VPS13Crs77555508G/A (S1798F)15: 619407260.6 2.47×10−121.44
ADCY3rs7586879C/T2: 2489410844.7 1.18×10−110.99
CTC1rs183966301G/A (A1025V)17: 82293840.2 1.84×10−110.56
SALL4rs77538589C/T (G117R)20: 517921344.3 2.24×10−112.61
ADCY7rs201661947G/A (A475T)16: 503044140.2 3.37×10−11ND
TP53INP1rs896854G/A8: 9494828330.7 3.76×10−110.96
TMEM245rs2271877C/T (A314T)9: 10909113211.5 8.13×10−110.98
FCRL1rs149740001A/T (K103I)1: 1578038561 8.38×10−111.22
SCYL2rs200554353T/C (M256T)12: 1003125681.6 8.73×10−110.36
TMCO3rs185071949C/T (P14L)13: 1134956220.4 1.97×10−101.02
WDR27rs3734905C/T6: 16955888619.5 3.00×10−100.99
NGBrs117207261C/G (Q60E)14: 772692380.5 3.59×10−101.24
rs6695567A/G1: 5316391343.6 4.89×10−101.05
FAPrs151314911C/T2: 1621737862.9 5.26×10−101.02
rs13277113A/G8: 1149167731.9 5.36×10−101.01
ACER1rs72981971T/C (M74V)19: 631227936.5 5.41×10−101
FREM2rs2496425T/C (F1070S)13: 3869055338.3 7.29×10−101.01
ASB13rs138695721A/C (V139G)10: 56490710.6 7.67×10−101.09
rs10943716T/C6: 8054335640.2 1.05×10−91.05
CCDC168rs1449707A/G (I3015T)13: 10274165344.4 1.82×10−90.98
ADGRV1rs2366928A/G (K3471E)5: 9072891817.9 2.17×10−90.97
MDN1rs115931523G/A (T3130M)6: 896959870.2 2.44×10−91.58
CD96rs140727933A/G (Y11C)3: 1115422800.2 2.63×10−91.16
rs4965121G/C15: 979755628.5 2.69×10−90.97
KNL1rs11858113T/C (M598T)15: 4062197927.4 3.29×10−90.99
OR4X2rs7120775C/G (Y27*)11: 4824518416.9 3.65×10−91.05
TRPM2rs144412484A/G (E450G)21: 443909340.4 3.68×10−90.82
MGAT5rs66523341C/T2: 13436480437 4.17×10−90.93
GCOM1rs4774980G/A15: 576916776 5.40×10−90.98
CSMD2rs1874045T/C (K2096R)1: 3360592545.1 7.56×10−91.08
ADAT1rs200524721G/C (Q167H)16: 756127850.3 8.75×10−91.12
rs4420065T/C1: 6569577813.1 1.52×10−81.02
NLRX1rs149129258C/A (P262Q)11: 1191740340.5 2.31×10−80.99
DNAAF3rs890871A/G (L280P)19: 5516068725.2 3.01×10−80.98
ZNF25rs150582814T/C (Y202C)10: 379528930.1 3.01×10−8ND
CMYA5rs62621915C/T (L1038F)5: 7973187749.3 5.62×10−80.93
SYDE2rs141587551C/A (D173Y)1: 852004806.6 6.20×10−81.03
SLC15A5rs3915247C/T12: 1624767947.6 7.01×10−80.97
CDC42BPGrs3741395T/C (Q1135R)11: 648300345.6 8.31×10−81.11
rs8030485G/A15: 7911659035.2 9.03×10−80.98
rs2564486G/T18: 5957906112.9 9.24×10−81.05
SLC4A4rs1062677A/C (I1074L)4: 715678286.1 1.20×10−70.91
STEAP1Brs17364464A/G7: 2247443412.2 2.03×10−70.99
KLF17rs11210969T/A (I35N)1: 4412937522.5 2.83×10−71.11
ADAMTS13rs78977446C/T (S903L)9: 1334457964.5 3.10×10−71.24
ZNF879rs17078988A/G (T112A)5: 17903228233.8 4.42×10−71.01
rs1464833T/C7: 12547560810.2 4.58×10−71.01
PKD1L1rs10951936A/T7: 4788208421 7.76×10−71.01
SNX32rs200684568G/A (G179R)11: 658507870.5 9.00×10−70.56
NTF3rs6332G/A12: 549446641.8 9.12×10−71.09
EFHD1rs4072149T/C2: 23264396736.7 9.18×10−71.01
URB2rs3811473G/T (G778V)1: 22963694615.8 9.48×10−71.05
CCDC71rs4955419A/T (Q317L)3: 491632594.6 9.91×10−71.13
rs543588T/G13: 3009684831.4 1.01×10−60.9
TRIM40rs757259G/A (E244K)6: 3014776516 1.06×10−61.12
rs3129264T/C6: 331338258.3 1.07×10−60.8
SEMA6Ars12516652G/T (D567E)5: 1164755521.9 1.14×10−61.2

{ label (or @symbol) needed for fn[@id='tfn2-mmr-16-05-5823'] } Allele frequencies were analyzed with Fisher's exact test. P<1.21×10−6 was considered to indicate a significant association.

a Major allele/minor allele. SNP, single nucleotide polymorphism; MAF, minor allele frequency; OR, odds ratio; ND, not determined. SNPs without gene names are not located within the genes.

Multivariable logistic regression analysis of the relation of SNPs to AF

The association of AF with 122 SNPs identified by the EWAS was further examined by multivariable logistic regression analysis with adjustment for age, sex and the prevalence of hypertension. Eight SNPs were identified to be associated with AF (P<0.01 in at least one genetic model; Table III). Among these eight SNPs, rs11552708 [G/A (G67R)] of TNF superfamily member 13 (TNFSF13), rs113710653 [C/T (E231 K)] of spermatogenesis and centriole associated 1 like (SPATC1L), and rs11231397 [G/C (R300T)] of solute carrier family 22 member 25 (SLC22A25) were significantly [P<1.02×10−4 (0.05/488)] associated with AF. The minor T allele of rs113710653 and the minor C allele of rs11231397 were risk factors for AF, whereas the minor A allele of rs11552708 was protective against this condition (Table III). In addition, the association of rs77538589 [C/T (G117R)] of spalt like transcription factor 4 (SALL4) with AF almost achieved significance, with the minor T allele representing a risk factor for this condition.

Table III.

Association of SNPs to atrial fibrillation as determined by multivariable logistic regression analysis.

Table III.

Association of SNPs to atrial fibrillation as determined by multivariable logistic regression analysis.

DominantRecessiveAdditive 1Additive 2




SNPNucleotide substitution (amino acid)P-valueOR (95% CI)P-valueOR (95% CI)P-valueOR (95% CI)P-valueOR (95% CI)
rs11552708G/A (G67R) 9.36×10−90.58 (0.49–0.70) 3.44×10−50.54 (0.39–0.73) 3.31×10−60.63 (0.52–0.77) 3.66×10−80.43 (0.30–0.59)
rs2076015T/C (R303G)0.1339 0.00281.38 (1.12–1.69)0.6253 0.00381.41 (1.12–1.78)
rs113710653C/T (E231K) 1.09×10−53.27 (2.00–5.14)0.6731 9.01×10−63.32 (2.02–5.21)0.6802
rs11231397G/C (R300T)0.00161.40 (1.13–1.73)0.00041.51 (1.20–1.88)0.031.28 (1.02–1.60) 3.71×10−51.77 (1.35–2.31)
rs77219745G/A (G1839D)0.2458 0.0076<0.01 (ND)0.4205 0.0072<0.01 (ND)
rs77538589C/T (G117R)0.00021.88 (1.36–2.56)0.4697 0.00021.90 (1.37–2.58)0.4825
rs141587551C/A (D173Y)0.7965 0.00463.39 (1.50–7.07)0.6564 0.00483.37 (1.49–7.02)
rs3129264T/C0.04220.78 (0.61–0.99)0.0016<0.01 (0.00–0.36)0.1138 0.0014<0.01 (ND)

[i] Multivariable logistic regression analysis was performed with adjustment for age, sex, and the prevalence of hypertension. Based on Bonferroni's correction, P<1.02×10−4 (0.05/488) were considered statistically significant and are shown in bold. SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval; ND, not determined.

Association of SNPs with intermediate phenotypes of AF

The association of the four SNPs, rs11552708, rs113710653, rs11231397 and rs77538589, with intermediate phenotypes of AF was examined. The rs11552708 SNP was significantly associated with hyperuricemia [P<0.0016 (0.05/32); Table IV].

Table IV.

P-values of the association between SNPs and intermediate phenotypes of atrial fibrillation.

Table IV.

P-values of the association between SNPs and intermediate phenotypes of atrial fibrillation.

SNPNucleotide substitution (amino acid)HypertensionDMHyper-TGHypo-HDLHyper-LDLCKDObesityHU
rs11552708G/A (G67R)0.99430.74390.82960.06220.89020.07540.4720.0006
rs113710653C/T (E231K)0.12780.52730.14650.68620.27320.42150.99860.1295
rs11231397G/C (R300T)0.24710.94450.120.16830.94150.95710.19930.9519
rs77538589C/T (G117R)0.11790.81130.46750.66560.58440.88150.87410.6006

[i] The prevalence of each phenotype was compared among genotypes with Fisher's exact test (2×2) or Pearson's chi-square test (2×3). Based on Bonferroni's correction, P<0.0016 (0.05/32) were considered statistically significant and are shown in bold. SNP, single nucleotide polymorphism; DM, diabetes mellitus; hyper-TG, hypertriglyceridemia; hypo-HDL, hypo-high density lipoprotein-cholesterolemia; hyper-LDL, hyper-low density lipoprotein-cholesterolemia; CKD, chronic kidney disease; HU, hyperuricemia.

Association of genes and SNPs identified in the present study with phenotypes previously examined in GWASs

The association of the four genes and four SNPs identified in the present study with phenotypes previously examined in GWASs available in public databases [GWAS Catalog (http://www.ebi.ac.uk/gwas) and GWAS Central (http://www.gwascentral.org/browser)] was determined. None of the genes or SNPs were identified to be associated with AF in previous studies (data not shown).

Discussion

The development of AF has been reported to be associated with endothelial dysfunction, including the upregulation of adhesion molecules, and increased inflammation and oxidative stress (15). Such endothelial dysfunction also promotes the electrophysiological remodeling observed in AF (16) and accelerates atrial ectopy through electrical excitation of cells near the pulmonary vein (17). In addition, coronary atherosclerosis may injure atrial tissue as result of myocardial ischemia and give rise to histological changes, including myocyte hypertrophy, necrosis and apoptosis (18,19). These processes may promote atrial remodeling with structural, functional, electrical, and metabolic consequences that underlie the development of AF (18,19).

The current study demonstrated that rs11552708 [G/A (G67R)] of TNFSF13, rs113710653 [C/T (E231K)] of SPATC1L and rs11231397 [G/C (R300T)] of SLC22A25 were significantly associated with AF in Japanese individuals. The minor T allele of rs113710653 and the minor C allele of rs11231397 were risk factors for AF, whereas the minor A allele of rs11552708 was protective against this condition. The rs77538589 [C/T (G117R)] SNP of SALL4 also showed a tendency to be associated to AF, although the association did not achieve statistical significance, with the minor T allele representing a risk factor for this condition.

The TNFSF13 gene is located at chromosomal region 17p13.1 (http://www.ncbi.nlm.nih.gov/gene/8741) and is expressed in various tissues and organs including heart muscle (The Human Protein Atlas, www.proteinatlas.org/ENSG00000161955-TNFSF13/tissue). TNFSF13 is a regulator of the development and function of B cells (20). It has been suggested that TNFSF13 has important roles in normal immune responses, and also in the development of autoimmune and inflammatory diseases. Plasma levels of TNFSF13 were reported to be increased in patients with multiple sclerosis (21), Sjögren's syndrome (22), systemic lupus erythematosus (23) and rheumatoid arthritis (24). The levels of TNFSF13 were also demonstrated to be significantly higher in the synovial fluid of patients with rheumatoid arthritis than in the serum, suggesting that this protein may be produced locally at the site of inflammation (25). Previous GWASs identified TNFSF13 as a susceptibility locus for immunoglobulin A nephropathy (26). Chronic vascular inflammation is an important feature of atherosclerosis (27), and atherosclerosis was previously reported to be associated with increased expression of TNFSF13. The plasma levels of TNFSF13 were demonstrated to be greater in patients with coronary artery disease than in healthy controls, and TNFSF13 immunoreactivity was detected in aggregated platelets within the ruptured plaque of patients with myocardial infarction (28). Pronounced TNFSF13 immunostaining was also observed in macrophages from carotid plaques (28). The current study demonstrated that rs11552708 [G/A (G67R)] of TNFSF13 was significantly associated with AF, with the minor A allele being protective against this condition. Chronic vascular inflammation, atherosclerosis, endothelial dysfunction and subsequent atrial remodeling are important in the pathogenesis of AF. Given that TNFSF13 is implicated in vascular inflammation and atherosclerosis (28), the association of TNFSF13 rs11552708 [G/A (G67R)] with AF may be attributable to the effect of this gene on atrial remodeling.

The SPATC1L gene is located at chromosomal region 21q22.3 (http://www.ncbi.nlm.nih.gov/gene/84221) and is expressed in various tissues and organs including heart muscle (The Human Protein Atlas, www.proteinatlas.org/ENSG00000160284-SPATC1L/tissue). SPATC1L is distributed in the cytoplasm, nucleus and perinuclear region of cells, and it translocates to sites of cell-cell junctions in response to cell stimulation with neurokinin A (29). Expression of SPATC1L modulates the response of human cells to N-methyl-N'-nitro-N-nitrosoguanidine, suggesting that the encoded protein may have a role in protecting cells from cell death induced by this DNA-damaging agent (30). The current study demonstrated that rs113710653 [C/T (E231K)] of SPATC1L was significantly associated with AF, with the minor T allele representing a risk factor for this condition, although the underlying molecular mechanism of this association remains unclear.

The SLC22A25 gene is located at chromosomal region 11q12.3 (http://www.ncbi.nlm.nih.gov/gene/387601) and is abundantly expressed in the liver (The Human Protein Atlas, www.proteinatlas.org/ENSG00000196600-SLC22A25/tissue). SLC22 is a large family of transmembrane organic cation and anion transporter proteins expressed predominantly in kidney and liver that mediate the uptake and excretion of environmental toxins, endogenous substances and drugs in the body (31). The results of the present study indicated that rs11231397 [G/C (R300T)] of SLC22A25 was significantly associated with AF, with the minor C allele representing a risk factor for this condition, although the functional relevance of this association remains unknown.

The SALL4 gene is located at chromosomal region 20q13.2 (http://www.ncbi.nlm.nih.gov/gene/57167) and is expressed in various tissues and organs including heart muscle (The Human Protein Atlas, www.proteinatlas.org/ENSG00000101115-SALL4/tissue). Mutations in SALL4 have been identified as a cause of Okihiro syndrome, which is characterized by developmental defects in limbs and multiple organs including the heart (32,33). Expression of SALL4 was detected in the left ventricular myocardium and interventricular septum of the mouse heart at embryonic day 10.5 (34). The current study demonstrated that rs77538589 [C/T (G117R)] of SALL4 tended to be associated with AF, with the minor T allele representing a risk factor for this condition. Given that SALL4 may have an important role in heart development, the association of rs77538589 of SALL4 with AF may reflect an effect of this gene on atrial structure or electrical function.

A meta-analysis of GWASs for AF (10) identified 10 SNPs (rs6666258 of potassium calcium-activated channel subfamily N member 3, rs3903239 of paired related homeobox 1, rs6817105 of paired like homeodomain 2, rs2040862 of Wnt family member 8A, rs3807989 of caveolin 1, rs10821415 of chromosome 9 open reading frame 3, rs10824026 of synaptopodin 2 like, rs1152591 of spectrin repeat containing nuclear envelope protein 2, rs7164883 of hyperpolarization activated cyclic nucleotide gated potassium channel 4 and rs2106261 of zinc finger homeobox 3) as susceptibility loci for this condition. The MAF of these SNPs ranged from 13.1–44.7% and the odds ratio from 0.87–1.64 (10). In the current study, four novel loci that may confer susceptibility to AF were identified, with the allele odds ratio (MAF, %) of rs11552708 of TNFSF13, rs113710653 of SPATC1L, rs11231397 of SLC22A25 and rs77538589 of SALL4 being 0.54 (36.6%), 5.36 (1.9%), 1.55 (43.7%) and 2.61 (4.3%), respectively. Both rs11552708 of TNFSF13 and rs11231397 of SLC22A25 were thus common variants with a low effect size, whereas rs113710653 of SPATC1L and rs77538589 of SALL4 were low-frequency variants with a moderate to high effect size.

There are several limitations to the present study: i) Given that the results of the study were not replicated, the findings will require validation in other independent subject panels or other ethnic groups; ii) it is possible that rs11552708 of TNFSF13, rs113710653 of SPATC1L, rs11231397 of SLC22A25 or rs77538589 of SALL4 is in linkage disequilibrium with other polymorphisms in the same gene or other nearby genes that are actually responsible for the development of AF; and iii) the functional relevance of these SNPs to the pathogenesis of AF remains to be elucidated.

In conclusion, rs11552708 of TNFSF13, rs113710653 of SPATC1L, rs11231397 of SLC22A25, and rs77538589 of SALL4 may be susceptibility loci for AF in Japanese individuals. Determination of genotypes of these SNPs may prove informative for assessment of the genetic risk for AF in Japanese individuals.

Acknowledgements

This work was supported by CREST (grant no. JPMJCR1302), Japan Science and Technology Agency (to Y.Y, J.S. and I.T.) and by Japan Society for the Promotion of Science KAKENHI grants no. JP15H04772 (to Y.Y.), no. JP25242062 (to M.T.) and no. JP16H01872 (to M.T.).

References

1 

Chugh SS, Havmoeller R, Narayanan K, Singh D, Rienstra M, Benjamin EJ, Gillum RF, Kim YH, McAnulty JH Jr, Zheng ZJ, et al: Worldwide epidemiology of atrial fibrillation: A global burden of disease 2010 study. Circulation. 129:837–847. 2014. View Article : Google Scholar : PubMed/NCBI

2 

Go AS, Hylek EM, Phillips KA, Chang Y, Henault LE, Selby JV and Singer DE: Prevalence of diagnosed atrial fibrillation in adults: National implications for rhythm management and stroke prevention: The An ticoagulation and risk factors in atrial fibrillation (ATRIA) study. JAMA. 285:2370–2375. 2001. View Article : Google Scholar : PubMed/NCBI

3 

Lip GY, Tse HF and Lane DA: Atrial fibrillation. Lancet. 379:648–661. 2012. View Article : Google Scholar : PubMed/NCBI

4 

Schnabel RB, Sullivan LM, Levy D, Pencina MJ, Massaro JM, D'Agostion RB Sr, Newton-Cheh C, Yamamoto JF, Magnani JW, Tadros TM, et al: Development of a risk score for atrial fibrillation (Framingham Heart Study): A community-based cohort study. Lancet. 373:739–745. 2009. View Article : Google Scholar : PubMed/NCBI

5 

Lubitz SA, Ozcan C, Magnani JW, Kääb S, Benjamin EJ and Ellinor PT: Genetics of atrial fibrillation: Implications for future research directions and personalized medicine. Circ Arrhythm Electrophysiol. 3:291–299. 2010. View Article : Google Scholar : PubMed/NCBI

6 

Gudbjartsson DF, Arnar DO, Helgadottir A, Gretarsdottir S, Holm H, Sigurdsson A, Jonasdottir A, Baker A, Thorleifsson G, Kristjansson K, et al: Variants conferring risk of atrial fibrillation on chromosome 4q25. Nature. 448:353–357. 2007. View Article : Google Scholar : PubMed/NCBI

7 

Gudbjartsson DF, Holm H, Gretarsdottir S, Thorleifsson G, Walters GB, Thorgeirsson G, Gulcher J, Mathiesen EB, Njølstad I, Nyrnes A, et al: A sequence variant in ZFHX3 on 16q22 associates with atrial fibrillation and ischemic stroke. Nat Genet. 41:876–878. 2009. View Article : Google Scholar : PubMed/NCBI

8 

Benjamin EJ, Rice KM, Arking DE, Pfeufer A, van Noord C, Smith AV, Schnabel RB, Bis JC, Boerwinkle E, Sinner MF, et al: Variants in ZFHX3 are associated with atrial fibrillation in individuals of European ancestry. Nat Genet. 41:879–881. 2009. View Article : Google Scholar : PubMed/NCBI

9 

Ellinor PT, Lunetta KL, Glazer NL, Pfeufer A, Alonso A, Chung MK, Sinner MF, de Bakker PI, Mueller M, Lubitz SA, et al: Common variants in KCNN3 are associated with lone atrial fibrillation. Nat Genet. 42:240–244. 2010. View Article : Google Scholar : PubMed/NCBI

10 

Ellinor PT, Lunetta KL, Albert CM, Glazer NL, Ritchie MD, Smith AV, Arking DE, Müller-Nurasyid M, Krijthe BP, Lubitz SA, et al: Meta-analysis identifies six new susceptibility loci for atrial fibrillation. Nat Genet. 44:670–675. 2012. View Article : Google Scholar : PubMed/NCBI

11 

Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, et al: Finding the missing heritability of complex diseases. Nature. 461:747–753. 2009. View Article : Google Scholar : PubMed/NCBI

12 

Grove ML, Yu B, Cochran BJ, Haritunians T, Bis JC, Taylor KD, Hansen M, Borecki IB, Cupples LA, Fornage M, et al: Best practices and joint calling of the HumanExome BeadChip: The CHARGE consortium. PLoS One. 8:e680952013. View Article : Google Scholar : PubMed/NCBI

13 

Anderson CA, Pettersson FH, Clarke GM, Cardon LR, Morris AP and Zondervan KT: Data quality control in genetic case-control association studies. Nat Protoc. 5:1564–1573. 2010. View Article : Google Scholar : PubMed/NCBI

14 

Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA and Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 38:904–909. 2006. View Article : Google Scholar : PubMed/NCBI

15 

O'Neal WT, Efird JT, Yeboah J, Nazarian S, Alonso A, Heckbert SR and Soliman EZ: Brachial flow-mediated dilation and incident atrial fibrillation: The multi-ethnic study of atherosclerosis. Arterioscler Thromb Vasc Biol. 34:2717–2720. 2014. View Article : Google Scholar : PubMed/NCBI

16 

Kim YM, Guzik TJ, Zhang YH, Zhang MH, Kattach H, Ratnatunga C, Pillai R, Channon KM and Casadei B: A myocardial Nox2 containing NAD (P)H oxidase contributes to oxidative stress in human atrial fibrillation. Circ Res. 97:629–636. 2005. View Article : Google Scholar : PubMed/NCBI

17 

Haïssaguerre M, Jaïs P, Shah DC, Takahashi A, Hocini M, Quiniou G, Garrigue S, Le Mouroux A, Le Métayer P and Clémenty J: Spontaneous initiation of atrial fibrillation by ectopic beats origination in the pulmonary veins. N Engl J Med. 339:659–666. 1998. View Article : Google Scholar : PubMed/NCBI

18 

Herringa J, Van Der Kuip DA, Hofman A, Kors JA, van Rooij FJ, Lip GY and Witteman JC: Subclinical atherosclerosis and risk of atrial fibrillation: The rotterdam study. Arch Intern Med. 167:382–387. 2007. View Article : Google Scholar : PubMed/NCBI

19 

Casaclang-Verzosa G, Gresh BJ and Tsang TS: Structural and functional remodeling of the left atrium: Clinical and therapeutic implications for atrial fibrillation. J Am Coll Cardiol. 51:1–11. 2008. View Article : Google Scholar : PubMed/NCBI

20 

Dillon SR, Gross JA, Ansell SM and Novak AJ: An APRIL to remember: Novel TNF ligands as therapeutic targets. Nat Rev Drug Discov. 5:235–246. 2006. View Article : Google Scholar : PubMed/NCBI

21 

Thangarajh M, Masterman T, Rot U, Duvefelt K, Brynedal B, Karrenbauer VD and Hillert J: Increased levels of APRIL (a proliferation-inducing ligand) mRNA in multiple sclerosis. J Neuroimmunol. 167:210–214. 2005. View Article : Google Scholar : PubMed/NCBI

22 

Jonsson MV, Szodoray P, Jellestad S, Jonsson R and Skarstein K: Association between circulating levels of the novel TNF family members APRIL and BAFF and lymphoid organization in primary Sjögren's syndrome. J Clin Immunol. 25:189–201. 2005. View Article : Google Scholar : PubMed/NCBI

23 

Koyama T, Tsukamoto H, Miyagi Y, Himeji D, Otsuka J, Miyagawa H, Harada M and Horiuchi T: Raised serum APRIL levels in patients with systemic lupus erythematosus. Ann Rheum Dis. 64:1065–1067. 2005. View Article : Google Scholar : PubMed/NCBI

24 

Cheema GS, Roschke V, Hilbert DM and Stohl W: Elevated serum B lymphocyte stimulator levels in patients with systemic immune-based rheumatic diseases. Arthritis Rheum. 44:1313–1319. 2001. View Article : Google Scholar : PubMed/NCBI

25 

Tan SM, Xu D, Roschke V, Perry JW, Arkfeld DG, Ehresmann GR, Migone TS, Hilbert DM and Stohl W: Local production of B lymphocyte stimulator protein and APRIL in arthritic joints of patients with inflammatory arthritis. Arthritis Rheum. 48:982–992. 2003. View Article : Google Scholar : PubMed/NCBI

26 

Kiryluk K, Li Y, Scolari F, Sanna-Cherchi S, Choi M, Verbitsky M, Fasel D, Lata S, Prakash S, Shapiro S, et al: Discovery of new risk loci for IgA nephropathy implicates genes involved in immunity against intestinal pathogens. Nat Genet. 46:1187–1196. 2014. View Article : Google Scholar : PubMed/NCBI

27 

Libby P: Inflammation in atherosclerosis. Nature. 420:868–874. 2002. View Article : Google Scholar : PubMed/NCBI

28 

Sandberg WJ, Otterdal K, Gullestad L, Halvorsen B, Ragnarsson A, Frøland SS, Damås JK, Oie E, Aukrust P, Hansson GK and Yndestad A: The tumour necrosis factor superfamily ligand APRIL (TNFSF13) is released upon platelet activation and expressed in atherosclerosis. Thromb Haemost. 102:704–710. 2009.PubMed/NCBI

29 

Lecat S, Matthes HW, Pepperkok R, Simpson JC and Galzi JL: A fluorescent live imaging screening assay based on translocation criteria identifies novel cytoplasmic proteins implicated in G protein-coupled receptor signaling pathways. Mol Cell Proteomics. 14:1385–1399. 2015. View Article : Google Scholar : PubMed/NCBI

30 

Fry RC, Svensson JP, Valiathan C, Wang E, Hogan BJ, Bhattacharya S, Bugni JM, Whittaker CA and Samson LD: Genomic predictors of interindividual differences in response to DNA damaging agents. Genes Dev. 22:2621–2626. 2008. View Article : Google Scholar : PubMed/NCBI

31 

Jacobsson JA, Haitina T, Lindblom J and Fredriksson R: Identification of six putative human transporters with structural similarity to the drug transporter SLC22 family. Genomics. 90:595–609. 2007. View Article : Google Scholar : PubMed/NCBI

32 

Kohlhase J, Heinrich M, Schubert L, Liebers M, Kispert A, Laccone F, Turnpenny P, Winter RM and Reardon W: Okihiro syndrome is caused by SALL4 mutations. Hum Mol Genet. 11:2979–2987. 2002. View Article : Google Scholar : PubMed/NCBI

33 

Sakaki-Yumoto M, Kobayashi C, Sato A, Fujimura S, Matsumoto Y, Takasato M, Kodama T, Aburatani H, Asashima M, Yoshida N and Nishinakamura R: The murine homolog of SALL4, a causative gene in Okihiro syndrome, is essential for embryonic stem cell proliferation, and cooperates with Sall1 in anorectal, heart, brain and kidney development. Development. 133:3005–3013. 2006. View Article : Google Scholar : PubMed/NCBI

34 

Koshiba-Takeuchi K, Takeuchi JK, Arruda EP, Kathiriya IS, Mo R, Hui CC, Srivastava D and Bruneau BG: Cooperative and antagonistic interactions between Sall4 and Tbx5 pattern the mouse limb and heart. Nat Genet. 38:175–183. 2006. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

November-2017
Volume 16 Issue 5

Print ISSN: 1791-2997
Online ISSN:1791-3004

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Yamada Y, Sakuma J, Takeuchi I, Yasukochi Y, Kato K, Oguri M, Fujimaki T, Horibe H, Muramatsu M, Sawabe M, Sawabe M, et al: Identification of TNFSF13, SPATC1L, SLC22A25 and SALL4 as novel susceptibility loci for atrial fibrillation by an exome‑wide association study. Mol Med Rep 16: 5823-5832, 2017
APA
Yamada, Y., Sakuma, J., Takeuchi, I., Yasukochi, Y., Kato, K., Oguri, M. ... Tanaka, M. (2017). Identification of TNFSF13, SPATC1L, SLC22A25 and SALL4 as novel susceptibility loci for atrial fibrillation by an exome‑wide association study. Molecular Medicine Reports, 16, 5823-5832. https://doi.org/10.3892/mmr.2017.7334
MLA
Yamada, Y., Sakuma, J., Takeuchi, I., Yasukochi, Y., Kato, K., Oguri, M., Fujimaki, T., Horibe, H., Muramatsu, M., Sawabe, M., Fujiwara, Y., Taniguchi, Y., Obuchi, S., Kawai, H., Shinkai, S., Mori, S., Arai, T., Tanaka, M."Identification of TNFSF13, SPATC1L, SLC22A25 and SALL4 as novel susceptibility loci for atrial fibrillation by an exome‑wide association study". Molecular Medicine Reports 16.5 (2017): 5823-5832.
Chicago
Yamada, Y., Sakuma, J., Takeuchi, I., Yasukochi, Y., Kato, K., Oguri, M., Fujimaki, T., Horibe, H., Muramatsu, M., Sawabe, M., Fujiwara, Y., Taniguchi, Y., Obuchi, S., Kawai, H., Shinkai, S., Mori, S., Arai, T., Tanaka, M."Identification of TNFSF13, SPATC1L, SLC22A25 and SALL4 as novel susceptibility loci for atrial fibrillation by an exome‑wide association study". Molecular Medicine Reports 16, no. 5 (2017): 5823-5832. https://doi.org/10.3892/mmr.2017.7334