Open Access

Identification of EGFLAM, SPATC1L and RNASE13 as novel susceptibility loci for aortic aneurysm in Japanese individuals by exome-wide association studies

  • 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: March 21, 2017     https://doi.org/10.3892/ijmm.2017.2927
  • Pages: 1091-1100
  • 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

We performed an exome-wide association study (EWAS) to identify genetic variants - in particular, low‑frequency or rare variants with a moderate to large effect size - that confer susceptibility to aortic aneurysm with 8,782 Japanese subjects (456 patients with aortic aneurysm, 8,326 control individuals) and with the use of Illumina HumanExome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays. The correlation of allele frequencies for 41,432 single nucleotide polymorphisms (SNPs) that passed quality control to aortic aneurysm was examined with Fisher's exact test. Based on Bonferroni's correction, a P-value of <1.21x10-6 was considered statistically significant. The EWAS revealed 59 SNPs that were significantly associated with aortic aneurysm. None of these SNPs was significantly (P<2.12x10-4) associated with aortic aneurysm by multivariable logistic regression analysis with adjustment for age, gender and hypertension, although 8 SNPs were related (P<0.05) to this condition. Examination of the correlation of these latter 8 SNPs to true or dissecting aortic aneurysm separately showed that rs1465567 [T/C (W229R)] of the EGF-like, fibronectin type III, and laminin G domains gene (EGFLAM) (dominant model; P=0.0014; odds ratio, 1.63) was significantly (P<0.0016) associated with true aortic aneurysm. We next performed EWASs for true or dissecting aortic aneurysm separately and found that 45 and 19 SNPs were significantly associated with these conditions, respectively. Multivariable logistic regression analysis with adjustment for covariates revealed that rs113710653 [C/T (E231K)] of the spermatogenesis- and centriole associated 1-like gene (SPATC1L) (dominant model; P=0.0002; odds ratio, 5.32) and rs143881017 [C/T (R140H)] of the ribonuclease A family member 13 gene (RNASE13) (dominant model; P=0.0006; odds ratio, 5.77) were significantly (P<2.78x10-4 or P<6.58x10-4, respectively) associated with true or dissecting aortic aneurysm, respectively. EGFLAM and SPATC1L may thus be susceptibility loci for true aortic aneurysm and RNASE13 may be such a locus for dissecting aneurysm in Japanese individuals.

Introduction

Aortic aneurysm is a serious condition that results from an atherosclerotic aorta and is a leading cause of mortality in humans (1). Studies on the genetic basis of familial aortic aneurysm have centered on the relationship between the condition to systemic connective tissue disorders such as Marfan syndrome (2) and Ehlers-Danlos syndrome (3). Although the molecular mechanism underlying nonsyndromic aortic aneurysm is complex and has not been determined definitively, several risk factors, including age, arteriosclerosis, hypertension and inflammatory or autoimmune diseases that affect the aorta, have been identified clinically (4,5). In addition to these conventional risk factors, recent studies have shown the importance of genetic factors in the development of sporadic aortic aneurysm by revealing a heritability of ~70% (6). Genome-wide association studies (GWASs) have uncovered several loci and genes that confer susceptibility to aortic aneurysm in European ancestry populations (712), but genetic variants that contribute to the development of this condition in Japanese individuals have not been identified definitively.

Genetic variants previously associated with aortic aneurysm typically have a minor allele frequency (MAF) of >10% and a small individual effect size (712). Given that these common variants explain only a small fraction of disease heritability, low-frequency (MAF of 0.5–5%) or rare (MAF of <0.5%) variants with a larger effect size may contribute to the genetic architecture of aortic aneurysm (13).

In the present study, we performed exome-wide association studies (EWASs) with the use of exome array-based genotyping methods to identify single nucleotide polymorphisms (SNPs) - in particular, low-frequency or rare coding variants with a moderate to large effect size - that confer susceptibility to aortic aneurysm in Japanese individuals. Given that most low-frequency or rare variants were not included in the arrays of previous GWASs, we used Illumina HumanExome-12 DNA Analysis BeadChip or Infinium Exome-24 BeadChip arrays, which provide coverage for functional SNPs including low-frequency or rare variants in entire exons.

Materials and methods

Study population

A total of 8,782 Japanese individuals (456 patients with aortic aneurysm, 8,326 controls) was examined. The subjects were recruited 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) either because they were experiencing 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 or in Tokyo or Kusatsu between 2011 and 2015; or from individuals who underwent autopsy at Tokyo Metropolitan Geriatric Hospital from 1995 to 2012.

True aortic aneurysm was defined as a permanent localized dilation of the aorta with a ≥50% increase in diameter relative to the expected normal size of the artery or with a diameter of >5 cm (14). Dissecting aortic aneurysm was defined as separation of the aortic wall layers with resulting true and false lumens or as intramural hematoma (15). The subjects with aortic aneurysm (279 with true aneurysm and 181 with dissecting aneurysm (four had both conditions) were examined by chest and abdominal X-ray and echocardiography followed by contrast medium-enhanced computed tomography. Some subjects were also examined by aortic angiography. Individuals with Marfan syndrome, Ehlers-Danlos syndrome, bicuspid aortic valve disease, aortitis syndrome, connective tissue disorder, congenital malformations of the heart or vessels, pseudoaneurysm, or traumatic aneurysm were excluded from the study. The control individuals had no history of aortic, coronary, or peripheral arterial disease; ischemic or hemorrhagic stroke; intracranial aneurysm; or other atherosclerotic, thrombotic, embolic or hemorrhagic disorders. Autopsy cases without aortic aneurysm were excluded from 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, Hirosaki University Graduate School of Medicine, Tokyo Metropolitan Institute of Gerontology, and participating hospitals. Written informed consent was obtained from each participant or from families of the deceased subjects.

EWASs

Venous blood (5 or 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 these cells either with the use of a DNA extraction kit (Genomix supplied by Talent, Trieste, Italy, or SMITEST EX-R&D supplied by Medical and Biological Laboratories, Nagoya, Japan) or by standard protocols based on phenol-chloroform extraction and spin columns. In autopsy cases, genomic DNA was extracted from kidneys. EWASs were performed for the 456 subjects with aortic aneurysm (or for the 279 subjects with true aneurysm or 181 subjects with dissecting aneurysm) and the 8,326 control subjects with the use of a HumanExome-12 v1.1 or v1.2 DNA Analysis BeadChip or an Infinium Exome-24 v1.0 BeadChip (Illumina, 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 (16). SNPs contained in only one of the exome arrays (~3.6%) were excluded from analysis. We performed quality control (17) as follows. i) Genotyping data with a call rate of <97% were discarded, with the mean call rate for the remaining data being 99.9%. ii) Gender specification was checked for each sample, with samples for which the gender designation in the clinical records was inconsistent with genetic sex being discarded. iii) Duplicated samples and cryptic relatedness were checked by calculation of identity by descent; all pairs with DNA samples showing identity by descent of >0.1875 were inspected and one sample from each pair was excluded. iv) The frequency of heterozygosity of SNPs was calculated for all samples, with those found to have extremely low or high heterozygosity (>3 standard deviations from the mean) being 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 in control individuals deviated significantly (P<0.001) from Hardy-Weinberg equilibrium were excluded. vii) The genotype data for each EWAS were examined for population stratification by principal components analysis (18), and population outliers were excluded from the analysis. A total of 41,432 SNPs passed quality control and was subjected to analysis.

Statistical analysis

Quantitative data for characteristics of the study subjects were compared between patients with aortic aneurysm 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 patients with aortic aneurysm and control subjects with Fisher's exact test. Multivariable logistic regression analysis was performed with aortic aneurysm as a dependent variable and independent variables including age, gender (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) models, which were analyzed simultaneously with a single statistical model. To compensate for multiple comparisons of genotypes with aortic aneurysm, we applied Bonferroni's correction for statistical significance of association. Given that 41,432 SNPs were finally examined, a P-value of <1.21×10−6 (0.05/41,432) was considered statistically significant. A quantile-quantile plot for P-values of allele frequencies in the EWAS for aortic aneurysm is shown in Fig. 1. The inflation factor (λ) was 1.57. P-values for other comparisons were similarly adjusted by Bonferroni's correction. Statistical tests were performed with JMP Genomics version 6.0 software (SAS Institute, Cary, NC, USA).

Results

Characteristics of the subjects

The characteristics of the subjects enrolled in the study are shown in Table I. Age, the frequency of males, and the prevalence of hypertension, diabetes mellitus, dyslipidemia, chronic kidney disease and hyperuricemia were significantly greater in patients with aortic aneurysm than in control individuals.

Table I

Characteristics of the 8,782 study subjects.

Table I

Characteristics of the 8,782 study subjects.

CharacteristicsAortic aneurysmControlP-value
No. of subjects4568326
Age (years)74.7±13.557.3±13.5<0.0001
Gender (male/female, %)64.0/36.051.4/48.6<0.0001
Body mass index (kg/m2)23.2±3.523.1±3.50.7068
Current or former smoker (%)44.139.80.1777
Hypertension (%)94.540.7<0.0001
Diabetes mellitus (%)44.714.7<0.0001
Dyslipidemia (%)67.956.80.0006
Chronic kidney disease (%)41.917.5<0.0001
Hyperuricemia (%)25.716.10.0003

[i] Quantitative data are expressed as means ± SD and were compared between patients with aortic aneurysm and control individuals with the unpaired Student's t-test. Categorical data were compared with Fisher's exact test. Based on Bonferroni's correction, a P-value of <0.0056 (0.05/9) was considered statistically significant.

EWAS of aortic aneurysm

We examined the correlation of allele frequencies for 41,432 SNPs that passed quality control to aortic aneurysm using the Fisher's exact test. A Manhattan plot for the EWAS of aortic aneurysm is shown in Fig. 2. After Bonferroni's correction, 59 SNPs were found to be significantly (P<1.21×10−6) associated with aortic aneurysm (Table II). The genotype distributions of these SNPs were in Hardy-Weinberg equilibrium (P>0.001) both among patients with aortic aneurysm and among the control individuals (data not shown).

Table II

The 59 single nucleotide polymorphisms (SNPs) significantly (P<1.21×10−6) associated with aortic aneurysm in an exome-wide association study (EWAS).

Table II

The 59 single nucleotide polymorphisms (SNPs) significantly (P<1.21×10−6) associated with aortic aneurysm in an exome-wide association study (EWAS).

GenedbSNPNucleotide (amino acid) substitutionaChromosome: positionMAF (%)P-value (allele)Allele OR
CATSPER4rs11247866A/G (Q77R)1:261913030.4 9.82×10−1471.81
RNASE13rs143881017C/T (R140H)14:210338700.5 2.49×10−1442.77
RNASE10rs202109789G/A (G87S)14:205107300.2 4.32×10−1230.47
rs2582513A/G14:10494845339.9 1.72×10−1180.86
HEATR1rs193150310T/A (V1975D)1:2365547520.3 3.09×10−1180.76
KIAA1217rs10828663G/A (A807T)10:2452452510.4 2.68×10−1001.00
MTUS1rs3739407G/A (R148C)8:1775536638.4 2.32×10−870.98
OR5W2rs75634103G/A11:5591452310.4 2.74×10−871.11
ALPK1rs2074379A/G (I732M)4:11243174332.0 1.03×10−861.11
ATAD5rs11657270T/C (Y1419H)17:3088736918.1 3.24×10−641.13
ACAT2rs25683A/G (K211R)6:15977531118.8 1.90×10−610.86
ZNF474rs201335566G/A (R253Q)5:1221527480.5 4.22×10−441.29
ZNF804Brs6963781A/G (M1105V)7:893362955.1 1.85×10−420.90
LOC100506679rs5751416G/A22:4303682026.3 1.62×10−410.81
SSPOrs191064068G/A (R209H)7:1497777381.1 8.55×10−351.13
ARHGEF28rs536568A/C5:7393584145.8 3.61×10−331.03
TMEM2rs142154818G/A (T1062M)9:717006451.0 6.41×10−331.90
HLA-DMBrs151719A/G6:3293612325.7 1.11×10−301.03
CCDC66rs61747994T/C (L802S)3:566193999.8 3.77×10−300.92
rs3135365T/G6:3242147818.9 1.62×10−280.84
NAA25rs12231744C/T (R876K)12:11203925135.1 5.02×10−281.07
RALGAPA2rs142962992G/C (E1676D)20:205054350.9 1.77×10−261.10
NEU1rs13118T/A6:318595099.7 6.92×10−221.20
AXDND1rs41267592C/T (T627M)1:1794685240.3 1.42×10−200.63
PHYKPLrs146105181T/C (N88D)5:1782300160.2 2.15×10−201.46
PCDH8rs5030685A/G (V743A)13:528462090.3 9.92×10−202.72
SELErs5361T/G (S149R)1:11697319193.3 2.49×10−170.95
MOV10L1rs760749A/C (I454L)22:5011725727.8 2.97×10−171.02
HHLA1rs75623295C/G (T90R)8:1320988932.9 1.23×10−160.80
TUBB1rs6070697G/A (R307H)20:5902434712.1 2.17×10−160.92
ZNF708rs504280C/T (R66Q)19:212945777.4 2.16×10−150.96
TICRRrs79501973G/A (V1373I)15:8962442714.7 2.51×10−150.97
ADNPrs148496595C/G (D924E)20:508919420.3 2.71×10−140.68
FCARrs11666735G/A (D113N)19:548855013.2 3.10×10−141.03
rs2823962G/A21:1667391332.8 7.80×10−140.93
EGFLAMrs1465567T/C (W229R)5:3837043525.1 5.92×10−131.19
rs1480347G/A8:2048994617.3 3.49×10−121.07
UBE4Brs180983516G/A (R331H)1:101063790.8 4.50×10−120.60
rs448705A/G8:1783719312.4 5.42×10−121.01
rs11970286C/T6:11835921117.3 1.38×10−110.90
rs10047727T/C13:2174305142.7 6.67×10−100.90
rs507856C/T3:16173615838.3 7.71×10−90.93
SLC1A6rs7253812C/A19:1498269126.7 9.43×10−90.99
FGBrs1800789G/A4:15456159113.3 1.86×10−80.99
SLC9A4rs1014286A/G (S784G)2:10253264143.9 1.90×10−81.13
HECTD4rs2074356C/T12:11220759725.4 2.22×10−81.03
PKD1L1rs66755489G/A (P2021L)7:478350322.9 2.50×10−81.16
CAMSAP1rs201291561T/C (N1062S)9:1358214760.2 3.04×10−81.14
C7orf43rs3800952C/T (R353Q)7:1001603316.3 4.08×10−81.08
ZNF671rs3746207G/A (A149V)19:5772164012.6 5.90×10−80.99
RIN3rs7150931T/C14:9267169646.2 1.26×10−71.03
rs10805579G/A5:1912741810.5 1.48×10−71.03
rs12546220T/C8:6946149329.1 1.82×10−70.96
DRD2rs12363125C/T11:1134151946.2 1.89×10−70.88
MTUS2rs17571410G/A13:2900748141.4 2.35×10−70.91
GALNTL5rs11766982A/G7:15199641727.6 2.56×10−71.05
POLErs5745022C/T12:13263239320.6 2.82×10−70.96
CHATrs3810947A/G10:4961319743.0 3.19×10−70.97
LILRB5rs117421142A/G (I420T)19:542523831.0 1.16×10−61.29

{ label (or @symbol) needed for fn[@id='tfn2-ijmm-39-05-1091'] } Allele frequencies were analyzed with Fisher's exact test.

a Major allele/minor allele. MAF, minor allele frequency; OR, odds ratio.

Multivariable logistic regression analysis of the correlation of SNPs to aortic aneurysm

The relation of the 59 identified SNPs to aortic aneurysm was examined further by multivariable logistic regression analysis with adjustment for age, gender and the prevalence of hypertension. Although 8 SNPs were related (P<0.05) to aortic aneurysm, no SNP was significantly [P<2.12×10−4 (0.05/236)] associated with this condition (Table III). We then examined the correlation of the 8 identified SNPs to true or dissecting aortic aneurysm separately. Five SNPs were related (P<0.05) to true aortic aneurysm (Table IV), among which rs1465567 [T/C (W229R)] of EGF-like, fibronectin type III, and laminin G domains gene (EGFLAM) was significantly [P<0.0016 (0.05/32)] associated with this condition, with the minor C allele representing a risk factor. No SNP was found to be related to dissecting aortic aneurysm (data not shown).

Table III

Correlation of single nucleotide polymorphisms (SNPs) to aortic aneurysm as determined by multivariable logistic regression analysis.

Table III

Correlation of single nucleotide polymorphisms (SNPs) to aortic aneurysm as determined by multivariable logistic regression analysis.

SNPDominant
Recessive
Additive 1
Additive 2
P-valueOR
(95% CI)
P-valueOR
(95% CI)
P-valueOR
(95% CI)
P-value (95%OR
(95% CI)
rs143881017C/T (R140H)0.02083.00 (1.20–6.80)0.96650.02073.00 (1.20–6.79)0.9667
rs5751416G/A0.03510.79 (0.63–0.98)0.20340.06970.1055
rs142154818G/A (T1062M)0.04861.94 (1.00–3.50)0.70440.04521.97 (1.02–3.54)0.7078
rs13118T/A0.04151.32 (1.01–1.71)0.41820.02451.37 (1.04–1.78)0.4879
rs5030685A/G (V743A)0.02932.94 (1.13–6.76)0.60670.02273.12 (1.19–7.22)0.6090
rs1465567T/C (W229R)0.00041.49 (1.19–1.85)0.69320.00041.51 (1.20–1.90)0.2807
rs7253812C/A0.08340.05430.01891.31 (1.05–1.65)0.1650
rs7150931T/C0.50330.02241.36 (1.04–1.75)0.92250.0586

[i] Multivariable logistic regression analysis was performed with adjustment for age, gender and the prevalence of hypertension. Based on Bonferroni's correction, a P-value of <2.12×10−4 (0.05/236) was considered statistically significant. OR, odds ratio; CI, confidence interval.

Table IV

Correlation of single nucleotide polymorphisms (SNPs) to true aortic aneurysm as determined by multivariable logistic regression analysis.

Table IV

Correlation of single nucleotide polymorphisms (SNPs) to true aortic aneurysm as determined by multivariable logistic regression analysis.

SNPDominant
Recessive
Additive 1
Additive 2
P-valueOR
(95% CI)
P-valueOR
(95% CI)
P-valueOR
(95% CI)
P-valueOR
(95% CI)
rs143881017C/T (R140H)0.2899ND0.2899ND
rs5751416G/A0.04480.73 (0.54–0.99)0.12210.10300.0663
rs142154818G/A (T1062M)0.22140.82090.21540.8231
rs13118T/A0.02531.51 (1.05–2.14)0.60760.01671.57 (1.09–2.23)0.6967
rs5030685A/G (V743A)0.07790.73140.06670.7333
rs1465567T/C (W229R)0.00141.63 (1.21–2.21)0.66700.00141.66 (1.22–2.27)0.2911
rs7253812C/A0.01221.47 (1.09–1.99)0.22650.00321.60 (1.17–2.18)0.5864
rs7150931T/C0.73290.0431 1.45(1.01–2.04)0.27370.2350

[i] Multivariable logistic regression analysis was performed with adjustment for age, gender and the prevalence of hypertension. Based on Bonferroni's correction, P-values of <0.0016 (0.05/32) were considered statistically significant and are shown in bold. OR, odds ratio; CI, confidence interval; ND, not determined.

EWASs of true or dissecting aortic aneurysm

We next examined the relation of allele frequencies for the total of 41,432 SNPs to true or dissecting aortic aneurysm separately with the use of Fisher's exact test. After Bonferroni's correction, 45 or 19 SNPs were found to be significantly (P<1.21×10−6) associated with true (Table V) or dissecting (Table VI) aortic aneurysm, respectively. The genotype distributions of these SNPs were in Hardy-Weinberg equilibrium (P>0.001) both among patients with true or dissecting aortic aneurysm and among control individuals (data not shown).

Table V

The 45 single nucleotide polymorphisms (SNPs) significantly (P<1.21×10−6) associated with true aortic aneurysm in an exome-wide association study (EWAS).

Table V

The 45 single nucleotide polymorphisms (SNPs) significantly (P<1.21×10−6) associated with true aortic aneurysm in an exome-wide association study (EWAS).

GenedbSNPNucleotide (amino acid) substitutionaChromosome: positionMAF (%)P-value (allele)Allele OR
KIAA1217rs10828663G/A (A807T)10:2452452510.4 6.87×10−941.14
NRAPrs79461687G/T (H1246Q)10:1136062471.3 1.58×10−901.67
OR5W2rs75634103G/A11:5591452310.4 4.84×10−811.17
ATAD5rs11657270T/C (Y1419H)17:3088736918.1 9.29×10−700.98
rs9683944A/G4:13751200810.8 1.65×10−560.83
TMPRSS3rs928302C/T (V53I)21:4238997527.6 1.34×10−520.99
ZNF804Brs6963781A/G (M1105V)7:893362955.1 6.67×10−471.01
ZNF474rs201335566G/A (R253Q)5:1221527480.5 5.79×10−420.84
LOC100506679rs5751416G/A22:4303682026.3 2.46×10−260.71
ARHGEF28rs536568A/C5:7393584145.8 3.10×10−231.03
RALGAPA2rs142962992G/C (E1676D)20:205054350.9 8.85×10−230.72
LYSMD1rs79024247G/T (Q150K)1:1511609745.2 7.43×10−210.88
CNGA1rs192912733C/T (R493Q)4:479372230.7 8.13×10−211.17
MOV10L1rs760749A/C (I454L)22:5011725727.8 1.23×10−201.18
HLA-DMBrs151719A/G6:3293612325.7 1.96×10−191.17
TUBB1rs6070697G/A (R307H)20:5902434712.1 3.20×10−180.97
rs3135365T/G6:3242147818.9 3.94×10−170.83
CCDC33rs1484214A/C15:7428873249.3 6.60×10−171.09
ZNF708rs504280C/T (R66Q)19:212945777.4 3.03×10−150.98
rs2823962G/A21:1667391332.8 3.36×10−140.97
SGCZrs1037934G/A8:143990659.0 3.44×10−131.22
DUOX1rs199549867A/T (R569S)15:451419970.2 1.58×10−100.81
CTSCrs3888798T/C (I453V)11:8829404116.5 3.21×10−101.06
C10orf128rs118189413C/G (H67D)10:491669086.3 4.04×10−100.84
OR5V1rs9405124A/G6:2940103619.9 4.15×10−100.98
UBAP2Lrs143080179T/C (S641P)1:1542551631.3 4.20×10−100.66
CAMSAP1rs201291561T/C (N1062S)9:1358214760.2 5.67×10−101.87
TICRRrs79501973G/A (V1373I)15:8962442714.7 7.76×10−101.04
NCAM1rs7111410C/T11:11317856515.9 2.63×10−90.80
IRGQrs3817C/A19:4358604347.7 3.30×10−90.92
SSPOrs191064068G/A (R209H)7:1497777381.1 2.35×10−81.67
SLC9A4rs1014286A/G (S784G)2:10253264143.9 2.84×10−81.11
rs2138852A/G17:293763312.5 2.92×10−80.77
SPATC1Lrs113710653C/T (E231K)21:461619211.9 3.91×10−87.39
DRD2rs12363125C/T11:1134151946.2 4.99×10−80.87
CCT5rs201280643C/G (S373C)5:102625840.4 7.04×10−80.99
AXDND1rs41267592C/T (T627M)1:1794685240.3 7.85×10−81.03
rs507856C/T3:16173615838.3 1.10×10−71.02
rs962040A/G8:1545436930.1 1.21×10−70.96
rs11970286C/T6:11835921117.3 1.60×10−70.83
NINLrs199671123C/T (A796T)20:254769050.2 2.45×10−70.85
rs12531488C/T7:14519499316.0 7.91×10−71.14
MROH7rs143029488G/C (A1313P)1:547101520.8 1.05×10−61.34
AFAP1rs28406288G/C (C403S)4:78005000.1 1.10×10−6ND
GALMrs6741892A/T (N190Y)2:3868982820.1 1.12×10−60.93

{ label (or @symbol) needed for fn[@id='tfn6-ijmm-39-05-1091'] } Allele frequencies were analyzed with Fisher's exact test.

a Major allele/minor allele. MAF, minor allele frequency; OR, odds ratio; ND, not determined.

Table VI

The 19 single nucleotide polymorphisms (SNPs) significantly (P<1.21×10−6) associated with dissecting aortic aneurysm in an exome-wide association study (EWAS).

Table VI

The 19 single nucleotide polymorphisms (SNPs) significantly (P<1.21×10−6) associated with dissecting aortic aneurysm in an exome-wide association study (EWAS).

GenedbSNPNucleotide (amino acid) substitutionaChromosome: positionMAF (%)P-value (allele)OR
ATXN7rs3774729G/A (V862M)3:6399640646.8 3.81×10−350.91
KIAA1217rs10828663G/A (A807T)10:2452452510.4 1.34×10−270.78
RNASE13rs143881017C/T (R140H)14:210338700.5 5.12×10−264.48
INPP5Frs3736822A/G (I453V)10:1198063971.7 8.50×10−230.62
rs9683944A/G4:13751200810.8 4.61×10−171.17
rs9610342A/G22:3573453030.7 8.20×10−170.79
FAM98Crs3745962C/A (T240K)19:3840560416.7 8.99×10−150.96
ZNF474rs201335566G/A (R253Q)5:1221527480.5 1.02×10−121.95
rs3135365T/G6:3242147818.9 1.21×10−120.86
DEPDC7rs34161108G/A (A192T)11:330277956.2 2.29×10−120.74
ARHGEF28rs536568A/C5:7393584145.8 5.53×10−121.04
RALGPS1rs57728614G/T (G383C)9:1271965839.5 5.73×10−110.84
AIM1Lrs34370465C/T (R847H)1:2634411823.9 3.63×10−101.09
ANXA7rs3750575C/T (R419Q)10:733789335.7 4.20×10−101.06
AXDND1rs41267592C/T (T627M)1:1794685240.3 1.77×10−8ND
rs2138852A/G17:293763312.5 1.49×10−71.32
CHATrs78925077C/G (S119R)10:496221090.5 1.50×10−71.93
GPR156rs902790A/T (E512D)3:1201679294.8 2.24×10−70.71
SELErs5361T/G (S149R)1:1697319193.3 3.58×10−70.98

{ label (or @symbol) needed for fn[@id='tfn8-ijmm-39-05-1091'] } Allele frequencies were analyzed with Fisher's exact test.

a Major allele/minor allele. MAF, minor allele frequency; OR, odds ratio; ND, not determined.

Multivariable logistic regression analysis of the correlation of SNPs to true or dissecting aortic aneurysm

The corrrelation of the 45 identified SNPs to true aortic aneurysm was examined further by multivariable logistic regression analysis with adjustment for age, gender and the prevalence of hypertension. Among these SNPs, rs113710653 [C/T (E231K)] of the spermatogenesis and centriole-associated 1-like gene (SPATC1L) was significantly [P<2.78×10−4 (0.05/180)] associated with true aortic aneurysm, with the minor T allele representing a risk factor for this condition (Table VII). The correlation of the 19 identified SNPs to dissecting aortic aneurysm was also further examined by multivariable logistic regression analysis with adjustment for age, gender and the prevalence of hypertension. The SNP rs143881017 [C/T (R140H)] of the ribonuclease A family member 13 gene (RNASE13) was significantly [P<6.58×10−4 (0.05/76)] associated with dissecting aortic aneurysm, with the minor T allele representing a risk factor for this condition (Table VII).

Table VII

Relation of single nucleotide polymorphisms (SNPs) to true or dissecting aortic aneurysm as determined by multivariable logistic regression analysis.

Table VII

Relation of single nucleotide polymorphisms (SNPs) to true or dissecting aortic aneurysm as determined by multivariable logistic regression analysis.

SNPDominant
Recessive
Additive 1
Additive 2
P-valueOR (95% CI)P-valueOR (95% CI)P-valueOR (95% CI)P-valueOR (95% CI)
True aortic aneurysm
rs113710653 C/T (E231K)
0.00025.32 (2.33–11.14)0.92630.00025.34 (2.34–11.18)0.9299
Dissecting aortic aneurysm
rs143881017 C/T (R140H)
0.00065.77 (2.25–12.95)0.96150.00065.48 (2.26–12.97)0.9621

[i] Multivariable logistic regression analysis was performed with adjustment for age, gender and the prevalence of hypertension. Based on Bonferroni's correction, P-values of <2.78×10−4 (0.05/180) or <6.58×10−4 (0.05/76) were considered statistically significant for true or dissecting aortic aneurysm, respectively, and are shown in bold. OR, odds ratio; CI, confidence interval.

Correlation of SNPs to intermediate phenotypes of aortic aneurysm

Finally, we examined the correlation of three SNPs (rs1465567, rs113710653 and rs143881017) to intermediate phenotypes (hypertension, diabetes mellitus, hypertriglyceridemia, hypo-HDL-cholesterolemia, hyper-LDL-cholesterolemia, chronic kidney disease, obesity and hyperuricemia) of aortic aneurysm. No SNP was found to be significantly [P<0.0021 (0.05/24)] associated with intermediate phenotypes (data not shown).

Discussion

True and dissecting aneurysms of the aorta develop as a result of progressive weakening of the vessel wall. They are associated with characteristic histological features including medial degeneration, which involves degeneration and fragmentation of elastic fibers as well as loss of smooth muscle cells and an accumulation of basophilic ground substances (19). In the present study, we showed that rs1465567 [T/C (W229R)] of EGFLAM and rs113710653 [C/T (E231K)] of SPATC1L were significantly associated with true aortic aneurysm, whereas rs143881017 [C/T (R140H)] of RNASE13 was significantly associated with dissecting aortic aneurysm, in Japanese individuals. The minor alleles of these SNPs were all risk factors for these conditions.

The EGFLAM is located at chromosomal region 5p13.2-p13.1 (NCBI Gene, https://www.ncbi.nlm.nih.gov/gene) and is expressed in various tissues and organs including vascular smooth muscle (The Human Protein Atlas, http://www.protein-atlas.org). EGFLAM is an extracellular matrix-like protein that colocalizes with both dystrophin and dystroglycan to the synaptic cleft of the photoreceptor ribbon synapse in the retina and which directly interacts with dystroglycan. It plays an important role in interactions between the photoreceptor ribbon synapse and bipolar dendrites (20,21), and it is implicated in defective photoreceptor synaptic function associated with congenital muscular dystrophies such as muscle-eye-brain disease caused by defective glycosylation of α-dystroglycan (22). A genome-wide pharmacogenomics study identified EGFLAM as a potential susceptibility locus for citalopram-induced side effects (23). We have now shown that rs1465567 [T/C (W229R)] of EGFLAM was significantly associated with true aortic aneurysm, with the minor C allele representing a risk factor for this condition, although the molecular mechanism underlying this association remains unclear.

The SPATC1L is located at chromosomal region 21q22.3 (NCBI Gene) and is expressed in various tissues and organs including vascular smooth muscle (The Human Protein Atlas). SPATC1L is distributed in the cytoplasm, nucleus, and perinuclear region of cells, and it translocates to the sites of cell-cell junctions in response to stimulation of cells with the neuropeptide neurokinin A (24). Expression of SPATC1L was also found to modulate the response of cells to N-methyl-N′-nitro-N-nitrosoguanidine and may thereby protect cells from cell death induced by this DNA-damaging agent (25). We demonstrated that rs113710653 [C/T (E231K)] of SPATC1L was significantly associated with true aortic aneurysm, with the minor T allele representing a risk factor for this condition, although the functional relevance of this association remains to be elucidated.

RNASE13 is located at chromosomal region 14q11.2 (NCBI Gene) and is expressed at a high level in the epididymis (The Human Protein Atlas). A GWAS showed that an SNP (rs3748348) located in the vicinity of RNASE13 was associated with executive functioning resilience (26). Gene-based analyses also revealed a genome-wide significant association between RNASE13 and executive functioning resilience (27). We now showed that rs143881017 [C/T (R140H)] of RNASE13 was significantly associated with dissecting aortic aneurysm, with the minor T allele representing a risk factor for this condition, although the molecular mechanism underpinning this association remains unknown.

Previous GWASs identified the SNPs: rs10757278 of CDKN2BAS, rs7025486 of DAP2IP, rs1466535 of LRP1, rs2118181 of FBN1, rs6511720 of LDLR and rs599839 of SORT1 as susceptibility loci for aortic aneurysm (712). The MAFs of these SNPs were >10%, and the odds ratios were 0.8–1.8 (5.7–12.28). We now identified three novel loci that may confer susceptibility to true or dissecting aortic aneurysm, with the odds ratios (MAF, %) of rs1465567 of EGFLAM, rs113710653 of SPATC1L, and rs143881017 of RNASE13 being 1.63 (25.1%), 5.32 (1.9%),and 5.77 (0.5%), respectively. Although rs1465567 of EGFLAM was a common variant with a small effect size, rs113710653 of SPATC1L and rs143881017 of RNASE13 were low-frequency variants with moderate to large effect sizes.

There are some limitations to the present study: i) Given that the number of subjects with aortic aneurysm was relatively small and the results of the study were not replicated, our findings will require validation with other independent subject panels or in other ethnic groups. ii) It is possible that rs1465567 of EGFLAM, rs113710653 of SPATC1L, or rs143881017 of RNASE13 is in linkage disequilibrium with other polymorphisms in the same gene or in other nearby genes that are actually responsible for the development of true or dissecting aneurysm. iii) The functional relevance of these SNPs to the pathogenesis of true or dissecting aneurysm remains to be elucidated.

In conclusion, rs1465567 of EGFLAM and rs113710653 of SPATC1L may be susceptibility loci for true aortic aneurysm and rs143881017 of RNASE13 may be such a locus for dissecting aortic aneurysm in Japanese individuals. Determination of genotypes for these SNPs may prove informative for assessment of the genetic risk for these conditions in Japanese individuals.

Acknowledgments

The present study was supported by CREST (H25-H30), Japan Science and Technology Agency (to Y.Y., J.S. and I.T.) and by Japan Society for the Promotion of Science KAKENHI grants JP15H04772 (to Y.Y.), JP25242062 (to M.T.) and JP16H01872 (to M.T.).

References

1 

Chiesa R, Melissano G, Civilini E, de Moura ML, Carozzo A and Zangrillo A: Ten years experience of thoracic and thoracoabdominal aortic aneurysm surgical repair: Lessons learned. Ann Vasc Surg. 18:514–520. 2004. View Article : Google Scholar : PubMed/NCBI

2 

Mizuguchi T, Collod-Beroud G, Akiyama T, Abifadel M, Harada N, Morisaki T, Allard D, Varret M, Claustres M, Morisaki H, et al: Heterozygous TGFBR2 mutations in Marfan syndrome. Nat Genet. 36:855–860. 2004. View Article : Google Scholar : PubMed/NCBI

3 

Schwarze U, Schievink WI, Petty E, Jaff MR, Babovic-Vuksanovic D, Cherry KJ, Pepin M and Byers PH: Haploinsufficiency for one COL3A1 allele of type III procollagen results in a phenotype similar to the vascular form of Ehlers-Danlos syndrome, Ehlers-Danlos syndrome type IV. Am J Hum Genet. 69:989–1001. 2001. View Article : Google Scholar : PubMed/NCBI

4 

Nordon IM, Hinchliffe RJ, Loftus IM and Thompson MM: Pathophysiology and epidemiology of abdominal aortic aneurysms. Nat Rev Cardiol. 8:92–102. 2011. View Article : Google Scholar

5 

Ye Z, Austin E, Schaid DJ and Kullo IJ: A multi-locus genetic risk score for abdominal aortic aneurysm. Atherosclerosis. 246:274–279. 2016. View Article : Google Scholar : PubMed/NCBI

6 

Harrison SC, Holmes MV, Agu O and Humphries SE: Genome wide association studies of abdominal aortic aneurysms - biological insights and potential translation applications. Atherosclerosis. 217:47–56. 2011. View Article : Google Scholar : PubMed/NCBI

7 

Helgadottir A, Thorleifsson G, Magnusson KP, Grétarsdottir S, Steinthorsdottir V, Manolescu A, Jones GT, Rinkel GJ, Blankensteijn JD, Ronkainen A, et al: The same sequence variant on 9p21 associates with myocardial infarction, abdominal aortic aneurysm and intracranial aneurysm. Nat Genet. 40:217–224. 2008. View Article : Google Scholar : PubMed/NCBI

8 

Gretarsdottir S, Baas AF, Thorleifsson G, Holm H, den Heijer M, de Vries JP, Kranendonk SE, Zeebregts CJ, van Sterkenburg SM, Geelkerken RH, et al: Genome-wide association study identifies a sequence variant within the DAB2IP gene conferring susceptibility to abdominal aortic aneurysm. Nat Genet. 42:692–697. 2010. View Article : Google Scholar : PubMed/NCBI

9 

Bown MJ, Jones GT, Harrison SC, Wright BJ, Bumpstead S, Baas AF, Gretarsdottir S, Badger SA, Bradley DT, Burnand K, et al CARDIoGRAM Consortium; Global BPgen Consortium; DIAGRAM Consortium; VRCNZ Consortium: Abdominal aortic aneurysm is associated with a variant in low-density lipoprotein receptor-related protein 1. Am J Hum Genet. 89:619–627. 2011. View Article : Google Scholar : PubMed/NCBI

10 

LeMaire SA, McDonald ML, Guo DC, Russell L, Miller CC III, Johnson RJ, Bekheirnia MR, Franco LM, Nguyen M, Pyeritz RE, et al: Genome-wide association study identifies a susceptibility locus for thoracic aortic aneurysms and aortic dissections spanning FBN1 at 15q21.1. Nat Genet. 43:996–1000. 2011. View Article : Google Scholar : PubMed/NCBI

11 

Bradley DT, Hughes AE, Badger SA, Jones GT, Harrison SC, Wright BJ, Bumpstead S, Baas AF, Grétarsdóttir S, Burnand K, et al: A variant in LDLR is associated with abdominal aortic aneurysm. Circ Cardiovasc Genet. 6:498–504. 2013. View Article : Google Scholar : PubMed/NCBI

12 

Jones GT, Bown MJ, Gretarsdottir S, Romaine SP, Helgadottir A, Yu G, Tromp G, Norman PE, Jin C, Baas AF, et al: A sequence variant associated with sortilin-1 (SORT1) on 1p13.3 is independently associated with abdominal aortic aneurysm. Hum Mol Genet. 22:2941–2947. 2013. View Article : Google Scholar : PubMed/NCBI

13 

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

14 

Johnston KW, Rutherford RB, Tilson MD, Shah DM, Hollier L and Stanley JC: Suggested standards for reporting on arterial aneurysms. Subcommittee on Reporting Standards for Arterial Aneurysms, Ad Hoc Committee on Reporting Standards, Society for Vascular Surgery and North American Chapter, International Society for Cardiovascular Surgery. J Vasc Surg. 13:452–458. 1991. View Article : Google Scholar : PubMed/NCBI

15 

Olsson C, Thelin S, Ståhle E, Ekbom A and Granath F: Thoracic aortic aneurysm and dissection: Increasing prevalence and improved outcomes reported in a nationwide population-based study of more than 14,000 cases from 1987 to 2002. Circulation. 114:2611–2618. 2006. View Article : Google Scholar : PubMed/NCBI

16 

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

17 

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

18 

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

19 

Barbour JR, Spinale FG and Ikonomidis JS: Proteinase systems and thoracic aortic aneurysm progression. J Surg Res. 139:292–307. 2007. View Article : Google Scholar : PubMed/NCBI

20 

Sato S, Omori Y, Katoh K, Kondo M, Kanagawa M, Miyata K, Funabiki K, Koyasu T, Kajimura N, Miyoshi T, et al: Pikachurin, a dystroglycan ligand, is essential for photoreceptor ribbon synapse formation. Nat Neurosci. 11:923–931. 2008. View Article : Google Scholar : PubMed/NCBI

21 

Kanagawa M, Omori Y, Sato S, Kobayashi K, Miyagoe-Suzuki Y, Takeda S, Endo T, Furukawa T and Toda T: Post-translational maturation of dystroglycan is necessary for pikachurin binding and ribbon synaptic localization. J Biol Chem. 285:31208–31216. 2010. View Article : Google Scholar : PubMed/NCBI

22 

Hu H, Li J, Zhang Z and Yu M: Pikachurin interaction with dystroglycan is diminished by defective O-mannosyl glycosylation in congenital muscular dystrophy models and rescued by LARGE overexpression. Neurosci Lett. 489:10–15. 2011. View Article : Google Scholar :

23 

Adkins DE, Clark SL, Åberg K, Hettema JM, Bukszár J, McClay JL, Souza RP and van den Oord EJ: Genome-wide pharmacogenomic study of citalopram-induced side effects in STAR*D. Transl Psychiatry. 2:e1292012. View Article : Google Scholar

24 

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

25 

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

26 

Mukherjee S, Kim S, Gibbons LE, Nho K, Risacher SL, Glymour MM, Habeck C, Lee GJ, Mormino E, Ertekin-Taner N, et al: Alzheimer's Disease Neuroimaging Initiative: Genetic architecture of resilience of executive functioning. Brain Imaging Behav. 6:621–633. 2012. View Article : Google Scholar : PubMed/NCBI

27 

Mukherjee S, Kim S, Ramanan VK, Gibbons LE, Nho K, Glymour MM, Ertekin-Taner N, Montine TJ, Saykin AJ and Crane PK: Alzheimer's Disease Neuroimaging Initiative: Gene-based GWAS and biological pathway analysis of the resilience of executive functioning. Brain Imaging Behav. 8:110–118. 2014. View Article : Google Scholar :

28 

van't Hof FN, Ruigrok YM, Lee CH, Ripke S, Anderson G, de Andrade M, Baas AF, Blankensteijn JD, Böttinger EP, Bown MJ, et al Aneurysm Consortium; Vascular Research Consortium of New Zealand: Shared genetic risk factors of intracranial, abdominal, and thoracic aneurysms. J Am Heart Assoc. 5:e0026032016. View Article : Google Scholar

Related Articles

Journal Cover

May-2017
Volume 39 Issue 5

Print ISSN: 1107-3756
Online ISSN:1791-244X

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 EGFLAM, SPATC1L and RNASE13 as novel susceptibility loci for aortic aneurysm in Japanese individuals by exome-wide association studies. Int J Mol Med 39: 1091-1100, 2017
APA
Yamada, Y., Sakuma, J., Takeuchi, I., Yasukochi, Y., Kato, K., Oguri, M. ... Tanaka, M. (2017). Identification of EGFLAM, SPATC1L and RNASE13 as novel susceptibility loci for aortic aneurysm in Japanese individuals by exome-wide association studies. International Journal of Molecular Medicine, 39, 1091-1100. https://doi.org/10.3892/ijmm.2017.2927
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 EGFLAM, SPATC1L and RNASE13 as novel susceptibility loci for aortic aneurysm in Japanese individuals by exome-wide association studies". International Journal of Molecular Medicine 39.5 (2017): 1091-1100.
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 EGFLAM, SPATC1L and RNASE13 as novel susceptibility loci for aortic aneurysm in Japanese individuals by exome-wide association studies". International Journal of Molecular Medicine 39, no. 5 (2017): 1091-1100. https://doi.org/10.3892/ijmm.2017.2927