Association of MicroRNA Biogenesis Genes Polymorphisms with Ischemic Stroke Susceptibility and Post-Stroke Mortality

Article information

J Stroke. 2018;20(1):110-121
Publication date (electronic) : 2018 January 31
doi : https://doi.org/10.5853/jos.2017.02586
aDepartment of Biomedical Science, CHA University College of Life Science, Seongnam, Korea
bDepartment of Emergency Medicine, CHA Bundang Medical Center, CHA University, Seongnam, Korea
cDepartment of Neurology, CHA Bundang Medical Center, CHA University, Seongnam, Korea
dDepartment of Neurosurgery, CHA Bundang Medical Center, CHA University, Seongnam, Korea
eDepartment of Internal Medicine, CHA University, Seongnam, Korea
Correspondence: Nam Keun Kim Department of Biomedical Science, CHA University College of Life Science, 335 Pangyo-ro, Bundang-gu, Seongnam 13488, Korea Tel: +82-31-881-7137 Fax: +82-31-881-7249 E-mail: nkkim@cha.ac.kr
Co-correspondence: Ok Joon Kim Department of Neurology, CHA Bundang Medical Center, CHA University, 59 Yatap-ro, Bundang-gu, Seongnam 13496, Korea Tel: +82-31-780-5481 Fax: +82-31-881-7249 E-mail: okjun77@cha.ac.kr
*These authors contributed equally to the manuscript as first author.
Received 2017 November 6; Revised 2017 December 30; Accepted 2018 January 19.

Abstract

Background and Purpose

MicroRNA (miRNA) expression has been examined in multiple conditions, including various cancers, neurological diseases, and cerebrovascular diseases, particularly stroke. Existing evidence indicates that miRNA biosynthesis and function play crucial roles in ischemic stroke physiology and pathology. In this study, we selected six known polymorphisms in miRNA-biogenesis genes; DICER rs13078A>T, rs3742330A>G; DROSHA rs10719T>C, rs6877842G>C; Ran GTPase (RAN) rs14035C>T; exportin 5 (XPO5) rs11077A>C.

Methods

We analyzed the associations between these polymorphisms and disease status and clinical factors in 585 ischemic stroke patients and 403 controls. Genotyping was performed with the polymerase chain reaction-restriction fragment length polymorphism method.

Results

The DICER rs3742330A>G (AA vs. AG+GG: adjusted odds ratio [AOR], 1.360; 95% confidence interval [CI], 1.024 to 1.807; P=0.034) and DROSHA rs10719T>C polymorphisms (TT vs. CC: AOR, 2.038; 95% CI, 1.113 to 3.730; P=0.021) were associated with ischemic stroke prevalence. During a mean follow-up of 4.80±2.11 years, 99 (5.91%) of the stroke patients died. In multivariate Cox proportional hazard regression models, a significant association was found between RAN rs14035 and survival of large artery disease patients with ischemic stroke (CC vs. TT: adjusted hazard ratio, 5.978; P=0.015).

Conclusions

An association was identified between the DICER and DROSHA polymorphisms and ischemic stroke. Specifically, polymorphisms (rs3742330 and rs10719) were more common in stroke patients, suggesting that they may be associated with an increased risk of ischemic stroke.

Introduction

Stroke is regarded as a complex, multifactorial, polygenic disease arising from a wide number of gene-gene and gene-environment interactions [1,2]. Multiple factors including hypertension, diabetes mellitus, smoking, hyperlipidemia, and hyperhomocysteinemia are associated with a higher risk of stroke [3,4]. Hyperhomocysteinemia, in particular, has been demonstrated to be an independent risk factor for ischemic stroke in several studies involving different ethnic groups [5,6].

MicroRNAs (miRNAs) are a class of endogenous, small, noncoding RNAs that pair with sites in 3’-untranslated regions (3’-UTRs) in mRNAs to downregulate their expression [7-9]. Previous studies have suggested that gene expression may be regulated by a small number of miRNAs [10-13]. To date, miRNA expression has been examined in patients with tumors [12,14,15], Alzheimer’s disease [16], Parkinson’s disease [17], schizophrenia [18], and stroke [19-21]. The evidence gathered to date indicates that miRNA biosynthesis plays crucial, physiological, and pathological roles [22-24]. Biosynthesis of miRNAs involves several miRNA biogenesis genes and occurs in multiple steps [7]. RNA polymerase II produces large primary miRNA transcripts (about 500 to 3,000 nucleotides) in the nucleus. The transcripts are processed by a multiprotein complex that includes DROSHA to form precursor miRNA (pre-miRNA) hairpins (about 60 to 100 nucleotides). After pre-miRNA has been exported to the cytoplasm by Ran GTPase (RAN) and exportin 5 (XPO5), it is further processed by DICER1, a polymerase II enzyme. Subsequently, the double-stranded miRNA duplex unwinds, forming an 18- to 24-nucleotide single-stranded, mature miRNA [7,25-27].

The present study tested the hypothesis that there is an association between miRNA biogenesis gene polymorphism and ischemic stroke risk. The objective was to investigate associations between six known miRNA biogenesis gene polymorphisms (DICER 3’-UTR rs13078A>T, DICER 3’-UTR rs3742330A>G, DROSHA 3’-UTR rs10719T>C, DROSHA 3’-UTR rs6877842G>C, RAN 3’-UTR rs14035C>T, and XPO5 3’-UTR rs11077A>C) and ischemic stroke and its risk factors.

Methods

Ethics statement

All study protocols were reviewed and approved by the Institutional Review Board of CHA Bundang Medical Center and followed the recommendations of the Declaration of Helsinki. Study subjects were recruited from the South Korean provinces of Seoul and Gyeonggi-do between 2000 and 2008. The Institutional Review Board of CHA Bundang Medical Center approved this genetic study in June 2000 (IRB No. 2013-09-073) and informed consent was obtained from study participants.

Study population

The Department of Neurology at CHA Bundang Medical Center, CHA University, referred 585 consecutive patients with ischemic stroke. Ischemic stroke was defined as a stroke (a clinical syndrome characterized by rapidly developing clinical symptoms and signs of focal or global loss of brain function) with evidence of cerebral infarction in clinically relevant areas of the brain according to magnetic resonance imaging (MRI) scan finding. Based on clinical manifestations and neuroimaging data, two neurologists classified all ischemic strokes into four causative subtypes using the Trial of Org 10172 in Acute Stroke Treatment (TOAST) criteria, as follows: (1) large artery disease (LAD), characterized by an infarction lesion ≥15 mm in diameter documented by MRI, and significant (>50%) stenosis of a major brain artery or a branch cortical artery documented by cerebral angiography with symptoms associated with that arterial territory; (2) small vessel disease (SVD), characterized by an infarction lesion <15 and ≥5 mm in diameter documented using MRI, and classic lacunar syndrome without evidence of cerebral cortical dysfunction or potentially detectable cardiac sources for embolism; (3) cardioembolism (CE) or arterial occlusions presumably due to an embolus arising in the heart, as detected by cardiac evaluation; and (4) undetermined pathogenesis, in which the cause of stroke could not be determined with any degree of confidence or involved >2 causes. Single and multiple (≥2 lesions) SVD cases were distinguished via brain MRI scans. The sizes and sites of cerebral infarctions were documented using MRI only. We selected 403 control subjects that were matched for sex ratio and age (within 5 years) in accordance with the patient group (Table 1). Controls were drawn from subjects visiting our hospitals during the same period for health examinations, including biochemical testing, electrocardiograms, and brain MRIs. Control subjects did not have a recent history of cerebrovascular disease or myocardial infarction. Exclusion criteria were the same as those used for the case group, as mentioned previously.

Baseline characteristics of ischemic stroke cases and controls

Genotyping

DNA was extracted from leukocytes using a G-DEX II Genomic DNA Extraction kit (Intron Biotechnology, Seongnam, Korea) according to the manufacturer’s instructions. The six best-studied single nucleotide polymorphisms (SNPs) in the miRNA biogenesis genes were determined through a documentary search that included 3′-UTR SNPs (DICER rs13078A>T, rs3742330A>G, DROSHA rs10719T>C, rs6877842G>C, RAN rs14035C>T, and XPO5 rs11077A>C). The miRNA biogenesis gene polymorphisms were analyzed by the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method. The PCR conditions for miRNA biogenesis genes polymorphism analyses are presented in Supplementary Table 1. To validate RFLP findings, 30% of the PCR assays for each polymorphism were randomly selected and repeated, followed by DNA sequencing. Sequencing was performed using an ABI 3730×l DNA Analyzer (Applied Biosystems, Foster City, CA, USA). The concordance of these quality control samples with the RFLP results was 100%.

Post-stroke mortality

To evaluate the association between miRNA biogenesis gene polymorphisms and long-term prognosis after ischemic stroke, survival time from stroke onset to death was tracked. The dates of death for each stroke patient (n=585) were ascertained using death certificates from the Korean National Statistical Office. Patients who were alive on December 31, 2013 were excluded from the study.

Statistical analysis

Genotype and allele combination frequencies in ischemic stroke cases and controls were compared using multivariate logistic regression models and Fisher exact test, respectively. Allele frequencies were calculated to identify deviations from Hardy-Weinberg equilibrium using P=0.05 as a threshold. Odds ratios, adjusted odds ratios (AORs), and 95% confidence intervals (CIs) were used to measure the strength of association between various genotypes and ischemic stroke. The association between miRNA biogenesis gene SNPs and post-stroke mortality was evaluated using Cox proportional hazard regression. The proportional hazards assumption was tested using a log(–log[survival]) plot and interaction for follow-up time in a time-dependent Cox regression model, which was found to be satisfactory. For multivariate analyses, logistic regression analyses were used to adjust for possible confounders, including age, sex, hypertension, diabetes mellitus, hyperlipidemia, and smoking. Statistical significance was accepted at the P<0.05 level [28,29].

Results

Baseline characteristics

The demographic characteristics of the 585 stroke cases and 403 controls are presented in Table 1. Of the stroke and control samples, 41.5% and 41.7%, respectively, were men, and the mean ages of stroke cases and controls were 62.7±10.9 and 62.8±10.6 years, respectively.

Genotype frequencies of the miRNA biogenesis genes polymorphisms

Table 2 provides the genotype distributions of the six miRNA biogenesis gene polymorphisms in ischemic stroke cases and controls. The DICER rs3742330A>G polymorphism was associated with greater odds of ischemic stroke (AA vs. GG: AOR, 1.459; 95% CI, 1.000 to 2.126; P=0.050; and AA vs. AG+GG: AOR, 1.360; 95% CI, 1.024 to 1.807; P=0.034). The DROSHA rs10719T>C polymorphism was also associated with greater odds of ischemic stroke (TT vs. CC: AOR, 2.038; 95% CI, 1.113 to 3.730; P=0.021; and TT+TC vs. CC: AOR, 2.001; 95% CI, 1.106 to 3.621; P=0.022). By contrast, the XPO5 rs11077A>C polymorphism was associated with lower odds of stroke (AA vs. CC: AOR, 0.101; 95% CI, 0.011 to 0.951; P=0.045; and AA vs. AC+CC: AOR, 0.669; 95% CI, 0.473 to 0.945; P=0.023). The frequency of the DICER1 rs13078A>T, DROSHA ABI 3730×l DNA Analyzer C>G, and RAN rs14035C>T polymorphisms was not significantly different between stroke cases and controls. To examine whether the effect of each polymorphism was confined to a specific subtype, stroke patients were separated into three subgroups (LAD, SVD, and CE) according to TOAST classifications (Table 3). Comparisons were also performed with control subjects and single versus multiple SVD patients (Supplementary Table 2). LAD was not significantly associated with any of the polymorphisms examined. However, the DICER1 rs3742330A>G polymorphism was significantly associated with SVD (AA vs. AG: AOR, 1.705; 95% CI, 1.073 to 2.710; P=0.024; and AA vs. AG+GG: AOR, 1.616; 95% CI, 1.041 to 2.509; P=0.032). In addition, the DROSHA rs10719T>C polymorphism was significantly associated with CE (TT vs. CC: AOR, 3.451; 95% CI, 1.264 to 9.422; P=0.016; and TT+TC vs. CC: AOR, 3.499; 95% CI, 1.348 to 9.082; P=0.010).

Comparison of DICER, DROSHA, RAN, and XPO5 polymorphisms between ischemic stroke patients and controls subjects

Comparison of DICER, DROSHA, RAN, and XPO5 polymorphisms between ischemic stroke subtypes and controls subjects

Combined effects of miRNA biogenesis gene polymorphisms and clinical factors

Stratified analysis of each clinical factor was performed to confirm the influence of clinical factors on the occurrence of ischemic stroke. However, no significant clinical factors were found to affect ischemic stroke risk (Supplementary Table 3). Therefore, a combined effect analysis was conducted to ascertain the effect of stroke and genotype on the prevalence of ischemic stroke. A synergistic effect was found for ischemic stroke prevalence between clinical factors (hypertension and diabetes mellitus) and miRNA biogenesis gene polymorphisms (Figure 1). The DROSHA rs10719 CC genotype was associated with stroke in individuals with hypertension (AOR, 4.781; 95% CI, 1.981 to 11.54). Diabetes mellitus combined with the DROSHA rs10719 CC genotype yielded the most significant association with stroke (AOR, 12.05; 95% CI, 1.541 to 94.19). Other gene-clinical factor combinations were not significantly associated with ischemic stroke (Supplementary Tables 4-6). The effects of miRNA biogenesis gene genotypes on blood coagulation status were evaluated by measuring platelet proportion, prothrombin time, activated partial thromboplastin time (aPTT), fibrinogen, and antithrombin (Supplementary Table 7). It was found that the DROSHA rs10719 CC genotype was significantly associated with elevated aPTT (TT vs. CC: P=0.007; TC vs. CC: P=0.019) (Supplementary Figure 1A) and antithrombin (TC vs. CC: P=0.039) (Supplementary Figure 1B). The other coagulant factors did not exhibit any statistically significant associations with any of the tested genotypes.

Figure 1.

The effects of DROSHA rs10719 T>C variant on ischemic stroke development modulated by clinical factors. (A) The synergistic effect for ischemic stroke susceptibility in DROSHA rs10719CC with hypertension (adjusted odds ratio [AOR], 4.781), DROSHA rs10719CC with non-hypertension (AOR, 1.967), or DROSHA rs10719TT with hypertension (AOR, 2.415). (B) DROSHA rs10719 T>C was associated with elevated ischemic stroke prevalence, in the case of DROSHA rs10719CC, with diabetes mellitus (AOR, 12.046), DROSHA rs10719CC, with non-diabetes mellitus (AOR, 1.770), and DROSHA rs10719TT, with diabetes mellitus (AOR, 2.351).

Polymorphisms in miRNA biogenesis genes versus post-stroke mortality

To evaluate the association between miRNA biogenesis gene polymorphisms and post-stroke mortality, Cox regression analysis was performed on the 585 patients with total ischemic stroke according to TOAST subtype (Figure 2 and Supplementary Table 8). During a mean follow-up of 4.80±2.11 years, 99 of the stroke patients died. In the multivariate Cox proportional hazard regression models, a significant association was found between RAN rs14035 and survival of LAD patients with ischemic stroke (CC vs. TT: adjusted hazard ratio [HR], 5.978; P=0.015; and CC+CT vs. TT: adjusted HR, 3.946; P=0.034) (Figure 2). A significant association was also found between RAN rs14035 and SVD in our analysis of ischemic stroke subtypes (CC vs. TT: adjusted HR, 9.403; P=0.015; and CC+CT vs. TT: adjusted HR, 5.223; P=0.039) (Figure 2). However, survival analysis was performed by Cox proportional-hazards regression based on the stepwise method for confirming covariant effect. A stepwise Cox regression analysis of ischemic stroke-related survival is shown in Supplementary Table 9. Mortality in SVD subgroup of ischemic stroke cases was associated with age and RAN rs14039 polymorphism status.

Figure 2.

Survival plot from a Cox proportional hazards model with RAN rs14035C>T polymorphisms in ischemic stroke. Survival curve of patients grouped by large artery disease (LAD) subtype based on (A) RAN rs14035CC vs. RAN rs14035TT genotypes and (B) RAN rs14035CC+CT vs. RAN rs14035TT genotypes. In addition, survival curve of patients grouped by SVD subtype based on (C) RAN rs14035CC vs. RAN rs14035TT genotypes and (D) RAN rs14035CC+CT vs. RAN rs14035TT genotypes. HR, hazard ratio.

Supplemental data

Gene-gene interaction analyses were performed for miRNA biogenesis gene polymorphisms to identify combinations that have synergistic effects on stroke risk (Supplementary Table 10 and 11). Some variants and allele combinations exhibited significant associations. However, the meaning of these associations should be interpreted with caution because the sample size is rather small.

Discussion

A recent study indicated that miR-221 and miR-222 modulate the angiogenic properties of human umbilical vein endothelial cells [30]. However, the function and biosynthesis of miRNAs in endothelial cell biology remains unclear. Therefore, in this study of ischemic stroke, we focused on miRNA biogenesis genes, such as DICER1, DROSHA, XPO5, and RAN, which are related to endothelial miRNA expression and angiogenesis [31,32]. We evaluated six polymorphisms in DICER1, DROSHA, XPO5, and RAN genes essential for miRNA biosynthesis [33-35], in ischemic stroke cases and controls. We found that polymorphisms in DICER and DROSHA, both of which are involved in angiogenesis and coagulation mechanisms [36-38], were linked to ischemic stroke. Polymorphisms in DICER, a gene already known to play a role in vascular growth and genesis, displayed the strongest association with ischemic stroke. For example, the DICER rs3742330 GG genotype was significantly more frequent in both overall stroke cases and subtype SVD cases than in controls. Moreover, the interplay between the DICER rs3742330 GG genotype and hyperlipidemia status was elevated stroke prevalence. The roles of DICER in angiogenesis and vascular growth have previously been investigated, and there is a reported association between DICER expression and healthy endothelial cell growth [39]. Moreover, there is strong evidence that a functional DICER1-dependent pathway is essential for a healthy endothelial angiogenic response. All major steps of the angiogenic process, including adhesion, proliferation, migration, and capillary-like structure formation are compromised by disrupted DICER1 signaling in cerebromicrovascular endothelial cells [40,41], in addition to other cell types [32,38,42,43].

DICER and DROSHA play crucial roles in vertebrate development. DICER1-deficient mice die early in development, between embryonic days 12.5 and 14.5, displaying impaired blood vessel and yolk sac formation. Similarly, zebrafish DICER mutant embryos display abnormal morphogenesis during gastrulation, brain formation, somatogenesis, and heart development [38]. In addition, loss of DROSHA leads to vascular smooth muscle cells disorder followed by hypoplastic blood vessel walls, cardiomyopathy, and liver hemorrhage in mice between embryonic days 13.5 and 14.5, and causes embryonic mortality in affected mice [44]. A number of studies have reported that DICER and DROSHA polymorphisms, including rs3742330 and rs10719, affect disease development and patient survival in various cancers [34,45-49]. In addition, functional analysis of rs1057035, which resides in the 3’-UTR of DICER, has revealed that the polymorphism affects hsa-miR-574-3p targeting and DICER expression [50]. The DROSHA rs10719 polymorphism, which is located in the 3’-UTR of DROSHA, was associated with different DROSHA expression levels [33] and presented different binding efficiency for the target site of hsa-miR-27b [35,51]. Furthermore, previous studies reported that DROSHA rs6877842 and rs640831 polymorphisms affected miRNA expression levels [34,52].

There are several potential mechanisms linking DICER and DROSHA polymorphisms to ischemic disease. First, DICER and DROSHA polymorphisms may directly affect angiogenesis via endothelial cell growth or induce blood vessel defects in embryos, resulting in vascular abnormalities [39,44]. In a DICER and DROSHA knockout model, hemorrhaging during vascular smooth muscle cell development was observed [39,44]. Moreover, DICER silencing in endothelial cells modulated the expression of several genes involved in endothelial biology, including nitric oxide synthase 3, matrix metalloproteinase 2 (MMP-2), integrins-v and -1, fibronectin, endothelin receptor types A, endothelin 1, vascular endothelial cadherin, and caspase-3. Both integrins-v and -1 are implicated in angiogenesis and endothelial survival [53], and MMP-2 participates in autocrine processes that influence hypoxia-induced migration and apoptotic death in endothelial cells [54]. Additionally, DROSHA has a similar role to DICER in vascular smooth muscle cell survival through ERK1/2 and AKT regulation [44]. Furthermore, previous studies identified that DROSHA influenced the regulation of miRNA expression. Transcription of certain miRNAs does not require DICER, but does need DROSHA (e.g., miR-1225 and miR-228) [55]. In addition, other miRNAs such as miR-877, miR-1224, and miR-1226 are independent of the canonical miRNA biogenesis pathway but dependent on the splicing process by DROSHA [56-58].

Alternately, DICER and DROSHA may indirectly affect miRNA regulation via RNA interference. As noted above, DICER and DROSHA are involved in miRNA biogenesis. miRNAs have multiple mRNA targets, are important regulators of gene expression, and play important roles in the initiation and progression of diverse diseases including leukemia, rheumatoid arthritis, and multiple sclerosis [59-62]. In particular, miRNAs are known to affect the immune system and vasculature in ischemic stroke [19,21,62]. At present, it is not known whether DICER polymorphisms affect stroke risk by affecting DICER enzyme function or via RNA interference. Further in vitro studies are needed to distinguish between these two hypotheses.

Interestingly, the results of the current study indicate a significant association between increased mortality after stroke and the RAN rs14035 C>T polymorphism, after adjusting for age, sex, hypertension, hyperlipidemia, and smoking status. As we did not have data on causes of death, we cannot be certain that the high mortality in patients with RAN rs14035 TT genotype was due to vascular events. However, analysis revealed that the RAN rs14035 TT genotype was significantly associated with the survival rate of ischemic stroke patients, supporting the likelihood of it also causing post-stroke mortality.

This was a case-control association study with 988 samples, but it has several limitations. First, although we found an association between XPO5 rs11077 polymorphisms and stroke risk, there is still no hypothesized mechanism for the role of this polymorphism in ischemic stroke prevalence. Second, the weak associations observed between the DICER rs13078, DROSHA rs6877842, and RAN rs14035 polymorphisms and ischemic stroke require replication. Third, some of the controls in our study were seeking medical attention; therefore, they were not completely healthy. However, recruitment of healthy participants with imaging and laboratory tests would markedly reduce the enrollment rate and including participants without imaging and laboratory tests may produce other vascular risk factor assessment biases. Finally, the study population was restricted to patients of Korean ethnicity.

Conclusions

We have identified an association between ischemic stroke susceptibility and polymorphisms in DICER rs3742330 and DROSHA rs10719, in addition to a significant association with the RAN rs14035 polymorphism in post-stroke mortality. These findings may encourage research efforts focusing on the role of DICER and DROSHA in vascular development. We postulate that the DICER rs3742330 and DROSHA rs10719 polymorphisms influence miRNA biosynthesis and therefore, miRNA post-transcriptional regulation during vascular endothelial cell growth, proliferation, and differentiation. However, the underlying mechanism remains to be elucidated in future research.

Supplementary materials

Supplementary materials related to this article can be found online at https://doi.org/10.5853/jos.2017.02586.

Supplementary Table 1.

PCR-RFLP condition for microRNA machinery genes polymorphism

jos-2017-02586-suppl1.pdf
Supplementary Table 2.

Comparison of genotype frequencies and AOR of DICER 6095 rs13078 A>T, DICER 9480 rs3742330 T>C, DROSHA -715 rs6877842 C>G, DROSHA 4576 rs10719 T>C, RAN 1857 rs14035 C>T, and XPO5 4485 rs11077 A>C polymorphisms between stroke subtype and controls

jos-2017-02586-suppl2.pdf
Supplementary Table 3.

AORs for ischemic stroke associated with microRNA machinery gene genotypes, stratified by clinical factors

jos-2017-02586-suppl3.pdf
Supplementary Table 4.

Adjusted odds ratios for ischemic stroke associated with DICER genotypes, combined by clinical factors

jos-2017-02586-suppl4.pdf
Supplementary Table 5.

Adjusted odds ratios for ischemic stroke associated with DROSHA genotypes, combined by clinical factors

jos-2017-02586-suppl5.pdf
Supplementary Table 6.

Adjusted odds ratios for ischemic stroke associated with RAN and XPO5 genotypes, combined by clinical factors

jos-2017-02586-suppl6.pdf
Supplementary Table 7.

Altered blood coagulation factors according to miRNA processing gene genotypes

jos-2017-02586-suppl7.pdf
Supplementary Table 8.

Genotype frequencies of microRNA processing genes polymorphism and ischemic stroke patients’ mortality according to TOAST subgroup

jos-2017-02586-suppl8.pdf
Supplementary Table 9.

Results of stepwise Cox regression analysis of ischemic stroke survival

jos-2017-02586-suppl9.pdf
Supplementary Table 10.

Frequency of DICERA and DROSHA genotype combinations predicted by multidimensional reduction in ischemic stroke cases and controls

jos-2017-02586-suppl10.pdf
Supplementary Table 11.

Allele combinations of DICER, DROSHA, RAN, and XPO5 polymorphisms between ischemic stroke patients and control subjects by multidimensional reduction method

jos-2017-02586-suppl11.pdf
Supplementary Figure 1.

Differences in activated partial thromboplastin time (aPTT) and antithrombin proportions based on DROSHA rs10719 T>C in ischemic stroke patients. Statistical analysis was performed using analysis of variance (ANOVA) test or Student t-test for each DROSHA rs10719 T>C genotype. (A) aPTT: the blood coagulation time was significantly different (P=0.005) between the DROSHA rs10719 TT (31.07±7.06), TC (31.20±7.30), and CC (36.23±35.28) genotypes. (B) Plasma antithrombin proportion: it was found that the DROSHA rs10719 T>C polymorphism affected the antithrombin proportion. The DROSHA rs10719CC genotype was associated with an elevated antithrombin percentage (97.32±27.29) compared with the DROSHA rs10719TT genotype (94.67±17.64), which had high antithrombin proportion relative to the DROSHA rs10719CC genotype (P=0.017). *P<0.05 calculated by ANOVA test; P<0.05 calculated by Student t-test.

jos-2017-02586-suppl12.pdf

Notes

Disclosure

The authors have no financial conflicts of interest.

Acknowledgements

This study was partially supported by National Research Foundation of Korea Grants funded by the Korean Government (NRF2015R1D1A1A09057432, NRF-2016R1D1A1B03930141, and NRF-2017R1D1A1B03029582) and by the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health and Welfare, Republic of Korea (grant number: HI16C1559).

References

1. Kluijtmans LA, Young IS, Boreham CA, Murray L, McMaster D, McNulty H, et al. Genetic and nutritional factors contributing to hyperhomocysteinemia in young adults. Blood 2003;101:2483–2488.
2. Gellekink H, den Heijer M, Heil SG, Blom HJ. Genetic determinants of plasma total homocysteine. Semin Vasc Med 2005;5:98–109.
3. Kitamura A, Iso H, Naito Y, Iida M, Konishi M, Folsom AR, et al. High-density lipoprotein cholesterol and premature coronary heart disease in urban Japanese men. Circulation 1994;89:2533–2539.
4. Iso H, Naito Y, Sato S, Kitamura A, Okamura T, Sankai T, et al. Serum triglycerides and risk of coronary heart disease among Japanese men and women. Am J Epidemiol 2001;153:490–499.
5. Kosch A, Koch HG, Heinecke A, Kurnik K, Heller C, Nowak-Göttl U, et al. Increased fasting total homocysteine plasma levels as a risk factor for thromboembolism in children. Thromb Haemost 2004;91:308–314.
6. Lee R, Frenkel EP. Hyperhomocysteinemia and thrombosis. Hematol Oncol Clin North Am 2003;17:85–102.
7. Ambros V. The functions of animal microRNAs. Nature 2004;431:350–355.
8. Kim JH, Li LH, Cai H, Nguyen VH, Min JJ, Shin BA, et al. miRNA-105 and -128 function as rheostats modulating MMP-2 activities by downregulation of TIMP-2 and upregulation of MT1-MMP. Genes Genom 2016;38:217–223.
9. Lee G, Choi YC, Byun Y, Yoon S, Jeong Y, Yoon J, et al. The seed sequence is necessary but insufficient for downregulation of target genes by miR-608. Genes Genom 2016;38:567–572.
10. Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet 2007;8:93–103.
11. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell 2004;116:281–297.
12. Guarnieri DJ, DiLeone RJ. MicroRNAs: a new class of gene regulators. Ann Med 2008;40:197–208.
13. Carthew RW, Sontheimer EJ. Origins and mechanisms of miRNAs and siRNAs. Cell 2009;136:642–655.
14. Nicoloso MS, Calin GA. MicroRNA involvement in brain tumors: from bench to bedside. Brain Pathol 2008;18:122–129.
15. Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, et al. MicroRNA expression profiles classify human cancers. Nature 2005;435:834–838.
16. Hill JM, Pogue AI, Lukiw WJ. Pathogenic microRNAs common to brain and retinal degeneration: recent observations in Alzheimer’s disease and age-related macular degeneration. Front Neurol 2015;6:232.
17. Kim J, Inoue K, Ishii J, Vanti WB, Voronov SV, Murchison E, et al. A microRNA feedback circuit in midbrain dopamine neurons. Science 2007;317:1220–1224.
18. Beveridge NJ, Tooney PA, Carroll AP, Gardiner E, Bowden N, Scott RJ, et al. Dysregulation of miRNA 181b in the temporal cortex in schizophrenia. Hum Mol Genet 2008;17:1156–1168.
19. Khoshnam SE, Winlow W, Farbood Y, Moghaddam HF, Farzaneh M. Emerging roles of microRNAs in ischemic stroke: as possible therapeutic agents. J Stroke 2017;19:166–187.
20. Jeyaseelan K, Lim KY, Armugam A. MicroRNA expression in the blood and brain of rats subjected to transient focal ischemia by middle cerebral artery occlusion. Stroke 2008;39:959–966.
21. Dharap A, Bowen K, Place R, Li LC, Vemuganti R. Transient focal ischemia induces extensive temporal changes in rat cerebral microRNAome. J Cereb Blood Flow Metab 2009;29:675–687.
22. Bazan HA, Hatfield SA, O’Malley CB, Brooks AJ, Lightell D Jr, Woods TC. Acute loss of miR-221 and miR-222 in the atherosclerotic plaque shoulder accompanies plaque rupture. Stroke 2015;46:3285–2387.
23. McManus DD, Freedman JE. MicroRNAs in platelet function and cardiovascular disease. Nat Rev Cardiol 2015;12:711–717.
24. Ksiazek-Winiarek DJ, Kacperska MJ, Glabinski A. MicroRNAs as novel regulators of neuroinflammation. Mediators Inflamm 2013;2013:172351.
25. Zamore PD, Haley B. Ribo-gnome: the big world of small RNAs. Science 2005;309:1519–1524.
26. Lund E, Güttinger S, Calado A, Dahlberg JE, Kutay U. Nuclear export of microRNA precursors. Science 2004;303:95–98.
27. Bartel DP, Chen CZ. Micromanagers of gene expression: the potentially widespread influence of metazoan microRNAs. Nat Rev Genet 2004;5:396–400.
28. Kim OJ, Hong SH, Oh SH, Kim TG, Min KT, Oh D, et al. Association between VEGF polymorphisms and homocysteine levels in patients with ischemic stroke and silent brain infarction. Stroke 2011;42:2393–2402.
29. Jeon YJ, Kim OJ, Kim SY, Oh SH, Oh D, Kim OJ, et al. Association of the miR-146a, miR-149, miR-196a2, and miR-499 polymorphisms with ischemic stroke and silent brain infarction risk. Arterioscler Thromb Vasc Biol 2013;33:420–430.
30. Poliseno L, Tuccoli A, Mariani L, Evangelista M, Citti L, Woods K, et al. MicroRNAs modulate the angiogenic properties of HUVECs. Blood 2006;108:3068–3071.
31. Kuehbacher A, Urbich C, Zeiher AM, Dimmeler S. Role of Dicer and Drosha for endothelial microRNA expression and angiogenesis. Circ Res 2007;101:59–68.
32. Saxena A, Tabin CJ. miRNA-processing enzyme Dicer is necessary for cardiac outflow tract alignment and chamber septation. Proc Natl Acad Sci U S A 2010;107:87–91.
33. Saeki M, Watanabe M, Inoue N, Tokiyoshi E, Takuse Y, Arakawa Y, et al. DICER and DROSHA gene expression and polymorphisms in autoimmune thyroid diseases. Autoimmunity 2016;49:514–522.
34. Rotunno M, Zhao Y, Bergen AW, Koshiol J, Burdette L, Rubagotti M, et al. Inherited polymorphisms in the RNA-mediated interference machinery affect microRNA expression and lung cancer survival. Br J Cancer 2010;103:1870–1874.
35. Zhang Y, Cao AL, Dong C. rs10719 polymorphism located within DROSHA 3’-untranslated region is responsible for development of primary hypertension by disrupting binding with microRNA-27b. Med Sci Monit 2017;23:911–918.
36. Sung H, Lee KM, Choi JY, Han S, Lee JY, Li L, et al. Common genetic polymorphisms of microRNA biogenesis pathway genes and risk of breast cancer: a case-control study in Korea. Breast Cancer Res Treat 2011;130:939–951.
37. Albinsson S, Suarez Y, Skoura A, Offermanns S, Miano JM, Sessa WC. MicroRNAs are necessary for vascular smooth muscle growth, differentiation, and function. Arterioscler Thromb Vasc Biol 2010;30:1118–1126.
38. Yang WJ, Yang DD, Na S, Sandusky GE, Zhang Q, Zhao G. Dicer is required for embryonic angiogenesis during mouse development. J Biol Chem 2005;280:9330–9335.
39. Ungvari Z, Tucsek Z, Sosnowska D, Toth P, Gautam T, Podlutsky A, et al. Aging-induced dysregulation of dicer1-dependent microRNA expression impairs angiogenic capacity of rat cerebromicrovascular endothelial cells. J Gerontol A Biol Sci Med Sci 2013;68:877–891.
40. Asada S, Takahashi T, Isodono K, Adachi A, Imoto H, Ogata T, et al. Downregulation of Dicer expression by serum withdrawal sensitizes human endothelial cells to apoptosis. Am J Physiol Heart Circ Physiol 2008;295:H2512–H2521.
41. Suárez Y, Fernández-Hernando C, Yu J, Gerber SA, Harrison KD, Pober JS, et al. Dicer-dependent endothelial microRNAs are necessary for postnatal angiogenesis. Proc Natl Acad Sci U S A 2008;105:14082–14087.
42. Suárez Y, Fernández-Hernando C, Pober JS, Sessa WC. Dicer dependent microRNAs regulate gene expression and functions in human endothelial cells. Circ Res 2007;100:1164–1173.
43. Lin J, Horikawa Y, Tamboli P, Clague J, Wood CG, Wu X. Genetic variations in microRNA-related genes are associated with survival and recurrence in patients with renal cell carcinoma. Carcinogenesis 2010;31:1805–1812.
44. Fan P, Chen Z, Tian P, Liu W, Jiao Y, Xue Y, et al. miRNA biogenesis enzyme Drosha is required for vascular smooth muscle cell survival. PLoS One 2013;8e60888.
45. Li X, Tian X, Zhang B, Zhang Y, Chen J. Variation in dicer gene is associated with increased survival in T-cell lymphoma. PLoS One 2012;7e51640.
46. Clague J, Lippman SM, Yang H, Hildebrandt MA, Ye Y, Lee JJ, et al. Genetic variation in MicroRNA genes and risk of oral premalignant lesions. Mol Carcinog 2010;49:183–189.
47. Merritt WM, Bar-Eli M, Sood AK. The dicey role of Dicer: implications for RNAi therapy. Cancer Res 2010;70:2571–2574.
48. Ma H, Yuan H, Yuan Z, Yu C, Wang R, Jiang Y, et al. Genetic variations in key microRNA processing genes and risk of head and neck cancer: a case-control study in Chinese population. PLoS One 2012;7e47544.
49. Khan S, Greco D, Michailidou K, Milne RL, Muranen TA, Heikkinen T, et al. MicroRNA related polymorphisms and breast cancer risk. PLoS One 2014;9e109973.
50. Stupack DG, Cheresh DA. Apoptotic cues from the extracellular matrix: regulators of angiogenesis. Oncogene 2003;22:9022–9029.
51. Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. The role of site accessibility in microRNA target recognition. Nat Genet 2007;39:1278–1284.
52. Martin-Guerrero I, Gutierrez-Camino A, Lopez-Lopez E, Bilbao-Aldaiturriaga N, Pombar-Gomez M, Ardanaz M, et al. Genetic variants in miRNA processing genes and pre-miRNAs are associated with the risk of chronic lymphocytic leukemia. PLoS One 2015;10e0118905.
53. Ben-Yosef Y, Miller A, Shapiro S, Lahat N. Hypoxia of endothelial cells leads to MMP-2-dependent survival and death. Am J Physiol Cell Physiol 2005;289:C1321–C1331.
54. Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro-RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A 2002;99:15524–15529.
55. Havens MA, Reich AA, Duelli DM, Hastings ML. Biogenesis of mammalian microRNAs by a non-canonical processing pathway. Nucleic Acids Res 2012;40:4626–4640.
56. Sibley CR, Seow Y, Saayman S, Dijkstra KK, El Andaloussi S, Weinberg MS, et al. The biogenesis and characterization of mammalian microRNAs of mirtron origin. Nucleic Acids Res 2012;40:438–448.
57. Dai L, Chen K, Youngren B, Kulina J, Yang A, Guo Z, et al. Cytoplasmic Drosha activity generated by alternative splicing. Nucleic Acids Res 2016;44:10454–10466.
58. Liu YP, Vink MA, Westerink JT, Ramirez de Arellano E, Konstantinova P, Ter Brake O, et al. Titers of lentiviral vectors encoding shRNAs and miRNAs are reduced by different mechanisms that require distinct repair strategies. RNA 2010;16:1328–1339.
59. Pauley KM, Satoh M, Chan AL, Bubb MR, Reeves WH, Chan EK. Upregulated miR-146a expression in peripheral blood mononuclear cells from rheumatoid arthritis patients. Arthritis Res Ther 2008;10:R101.
60. Du C, Liu C, Kang J, Zhao G, Ye Z, Huang S, et al. MicroRNA miR-326 regulates TH-17 differentiation and is associated with the pathogenesis of multiple sclerosis. Nat Immunol 2009;10:1252–1259.
61. O’Connell RM, Rao DS, Chaudhuri AA, Baltimore D. Physiological and pathological roles for microRNAs in the immune system. Nat Rev Immunol 2010;10:111–122.
62. Liu DZ, Tian Y, Ander BP, Xu H, Stamova BS, Zhan X, et al. Brain and blood microRNA expression profiling of ischemic stroke, intracerebral hemorrhage, and kainate seizures. J Cereb Blood Flow Metab 2010;30:92–101.

Article information Continued

Figure 1.

The effects of DROSHA rs10719 T>C variant on ischemic stroke development modulated by clinical factors. (A) The synergistic effect for ischemic stroke susceptibility in DROSHA rs10719CC with hypertension (adjusted odds ratio [AOR], 4.781), DROSHA rs10719CC with non-hypertension (AOR, 1.967), or DROSHA rs10719TT with hypertension (AOR, 2.415). (B) DROSHA rs10719 T>C was associated with elevated ischemic stroke prevalence, in the case of DROSHA rs10719CC, with diabetes mellitus (AOR, 12.046), DROSHA rs10719CC, with non-diabetes mellitus (AOR, 1.770), and DROSHA rs10719TT, with diabetes mellitus (AOR, 2.351).

Figure 2.

Survival plot from a Cox proportional hazards model with RAN rs14035C>T polymorphisms in ischemic stroke. Survival curve of patients grouped by large artery disease (LAD) subtype based on (A) RAN rs14035CC vs. RAN rs14035TT genotypes and (B) RAN rs14035CC+CT vs. RAN rs14035TT genotypes. In addition, survival curve of patients grouped by SVD subtype based on (C) RAN rs14035CC vs. RAN rs14035TT genotypes and (D) RAN rs14035CC+CT vs. RAN rs14035TT genotypes. HR, hazard ratio.

Table 1.

Baseline characteristics of ischemic stroke cases and controls

Characteristic Cases (n=585) Controls (n=403) P*
Male sex 243 (41.5) 168 (41.7) 1.000
Age (yr) 62.72±10.91 62.79±10.61 0.249
Smokers 206 (35.2) 133 (33.0) 0.654
Hypertension 367 (62.7) 164 (40.7) 0.0002
Diabetes mellitus 156 (26.7) 52 (12.9) <0.0001
Hyperlipidemia 178 (30.4) 94 (23.3) 0.069
tHcy (μmol/L) 11.18±6.81 10.06±4.20 0.004
Folate (nmol/L) 6.99±5.14 8.55±5.96 <0.0001
Vitamin B12 (pg/mL) 747.10±616.70 744.70±669.67 0.905
Total cholesterol (mg/dL) 191.09±40.39 193.28±37.74 0.353
Triglycerides (mg/dL) 154.39±114.25 147.24±90.61 0.592
PLT (103/μL) 249.22±87.57 243.10±67.35 0.240
PT (sec) 11.78±0.98 11.77±0.80 0.875
aPTT (sec) 30.47±4.43 33.42±18.58 <0.001
Fibrinogen (mg/dL) 424.74±130.82 94.45±44.03 <0.001
Antithrombin (%) 94.09±18.82 400.17±120.45 <0.001
BUN (mg/dL) 15.86±6.10 15.78±4.93 0.828
Uric acid (mg/dL) 4.67±1.55 4.64±1.44 0.733

Values are presented as number (%) or mean±SD.

tHcy. total homocysteine; PLT, platelet count; PT, prothrombin time; aPTT, activated partial thromboplastin time; BUN, blood urea nitrogen.

*

P-values were calculated using the two-sided t-test for continuous variables and the chi-square test for categorical variables.

Table 2.

Comparison of DICER, DROSHA, RAN, and XPO5 polymorphisms between ischemic stroke patients and controls subjects

Genotype Controls (n=403) Cases (n=585) COR (95% CI) P* P AOR (95% CI) P§ P
DICER rs13078 A>T
 AA 360 (89.3) 527 (90.1) 1.000 (reference) 1.000 (reference)
 AT 43 (10.7) 55 (9.4) 0.874 (0.574–1.331) 0.530 0.530 0.926 (0.596–1.439) 0.733 0.733
 TT 0 3 (0.5) NA NA 0.994 0.994
 Dominant (AA vs. AT+TT) 0.921 (0.608–1.398) 0.700 0.700 0.978 (0.633–1.511) 0.920 0.920
 Recessive (AA+AT vs. TT) NA NA 0.994 0.994
 HWE-P 0.258 0.238
DICER rs3742330 A>G
 AA 148 (36.7) 169 (28.9) 1.000 (reference) 1.000 (reference)
 AG 180 (44.7) 280 (47.9) 1.362 (1.020–1.820) 0.036 0.129 1.313 (0.969–1.779) 0.079 0.237
 GG 75 (18.6) 136 (23.2) 1.588 (1.110–2.272) 0.011 0.043 1.459 (1.000–2.126) 0.050 0.100
 Dominant (AA vs. AG+GG) 1.429 (1.090–1.872) 0.010 0.057 1.360 (1.024–1.807) 0.034 0.102
 Recessive (AA+AG vs. GG) 1.325 (0.966–1.817) 0.081 0.135 1.254 (0.902–1.745) 0.178 0.356
 HWE-P 0.125 0.337
DROSHA rs6877842 C>G
 CC 371 (92.1) 548 (93.7) 1.000 (reference) 1.000 (reference)
 CG 31 (7.7) 36 (6.2) 0.786 (0.478–1.294) 0.344 0.503 0.785 (0.467–1.320) 0.361 0.542
 GG 1 (0.2) 1 (0.2) 0.677 (0.042–10.858) 0.783 0.783 0.769 (0.046–12.813) 0.855 0.994
 Dominant (CC vs. CG+GG) 0.783 (0.479–1.279) 0.328 0.394 0.784 (0.470–1.309) 0.352 0.422
 Recessive (CC+CG vs. GG) 0.688 (0.043–11.038) 0.792 0.792 0.766(0.046–12.743) 0.852 0.994
 HWE-P 0.680 0.614
DROSHA rs10719 T>C
 TT 228 (56.6) 304 (52.0) 1.000 (reference) 1.000 (reference)
 TC 158 (39.2) 235 (40.2) 1.116 (0.856–1.454) 0.419 0.503 1.102 (0.835–1.455) 0.492 0.590
 CC 17 (4.2) 46 (7.9) 2.029 (1.134–3.633) 0.017 0.043 2.038 (1.113–3.730) 0.021 0.994
 Dominant (TT vs. TC+CC) 1.204 (0.933–1.554) 0.153 0.306 1.193 (0.913–1.558) 0.196 0.294
 Recessive (TT+TC vs. CC) 1.938 (1.094–3.432) 0.023 0.115 2.001 (1.106–3.621) 0.022 0.132
 HWE-P 0.107 0.950
RAN rs14035 C>T
 CC 240 (59.6) 369 (63.1) 1.000 (reference) 1.000 (reference)
 CT 149 (37.0) 192 (32.8) 0.838 (0.641–1.097) 0.198 0.396 0.803 (0.606–1.064) 0.127 0.254
 TT 14 (3.5) 24 (4.1) 1.115 (0.566–2.198) 0.753 0.783 1.106 (0.545–2.244) 0.780 0.994
 Dominant (CC vs. CT+TT) 0.862 (0.664–1.118) 0.263 0.394 0.830 (0.632–1.091) 0.181 0.294
 Recessive (CC+CT vs. TT) 1.189 (0.607–2.327) 0.614 0.768 1.198 (0.597–2.403) 0.611 0.917
 HWE-P 0.114 0.876
XPO5 rs11077 A>C
 AA 319 (79.2) 497 (85.0) 1.000 (reference) 1.000 (reference)
 AC 79 (19.6) 87 (14.9) 0.707 (0.505–0.989) 0.043 0.129 0.707 (0.497–1.005) 0.053 0.237
 CC 5 (1.2) 1 (0.2) 0.128 (0.015–1.104) 0.062 0.103 0.101 (0.011–0.951) 0.045 0.100
 Dominant (AA vs. AC+CC) 0.672 (0.483–0.936) 0.019 0.057 0.669 (0.473–0.945) 0.023 0.102
 Recessive (AA+AC vs. CC) 0.136 (0.016–1.171) 0.069 0.135 0.116 (0.013–1.078) 0.058 0.174
 HWE-P 0.965 0.161

Values are presented as number (%).

RAN, Ran GTPase; XPO5, exportin 5; COR, crude odds ratio; CI, confidence interval; AOR, adjusted odds ratio; NA, not available; HWE, Hardy-Weinberg equilibrium.

*

Calculated by chi-square test according to genotype frequencies;

P-value calculated by false discovery rate test;

Odds ratios adjusted for age, sex, hypertension, diabetes mellitus, hyperlipidemia, and smoking status;

§

P-value calculated by logistics regression analysis.

Table 3.

Comparison of DICER, DROSHA, RAN, and XPO5 polymorphisms between ischemic stroke subtypes and controls subjects

Genotype Controls (n=403) Stroke: LAD (n=200) AOR (95% CI)* P P Stroke: SVD (n=149) AOR (95% CI)* P P Stroke: CE (n=54) AOR (95% CI)* P P
DICER rs13078 A>T
 AA 360 (89.3) 186 (93.0) 1.000 (reference) 130 (87.2) 1.000 (reference) 49 (90.7) 1.000 (reference)
 AT 43 (10.7) 14 (7.0) 0.597 (0.309–1.156) 0.126 0.378 19 (12.8) 1.368 (0.747–2.506) 0.310 0.465 4 (7.4) 0.621 (0.212–1.823) 0.386 0.701
 TT 0 0 NA 0 NA 1 (1.9) NA 0.995 0.996
 Dominant (AA vs. AT+TT) 0.597 (0.309–1.156) 0.126 0.378 1.368 (0.747–2.506) 0.310 0.465 0.778 (0.291–2.078) 0.617 0.779
 Recessive (AA+AT vs. TT) NA NA NA 0.995 0.998
DICER rs3742330 A>G
 AA 148 (36.7) 62 (31.0) 1.000 (reference) 40 (26.8) 1.000 (reference) 18 (33.3) 1.000 (reference)
 AG 180 (44.7) 89 (44.5) 1.084 (0.715–1.643) 0.703 0.844 78 (52.3) 1.705 (1.073–2.710) 0.024 0.144 20 (37.0) 0.885 (0.448–1.750) 0.726 0.726
 GG 75 (18.6) 49 (24.5) 1.448 (0.882–2.375) 0.143 0.600 31 (20.8) 1.412 (0.785–2.539) 0.249 0.623 16 (29.6) 1.535 (0.723–3.257) 0.264 0.792
 Dominant (AA vs. AG+GG) 1.187 (0.809–1.741) 0.380 0.535 1.616 (1.041–2.509) 0.032 0.192 1.099 (0.597–2.021) 0.763 0.779
 Recessive (AA+AG vs. GG) 1.385 (0.899–2.133) 0.140 0.655 1.041 (0.635–1.708) 0.873 0.993 1.681 (0.879–3.215) 0.116 0.348
DROSHA rs6877842 C>G
 CC 371 (92.1) 188 (94.0) 1.000 (reference) 143 (96.0) 1.000 (reference) 48 (88.9) 1.000 (reference)
 CG 31 (7.7) 12 (6.0) 0.777 (0.379–1.589) 0.489 0.786 5 (3.4) 0.428 (0.159–1.153) 0.093 0.279 6 (11.1) 1.416 (0.555–3.610) 0.467 0.701
 GG 1 (0.2) 0 NA 0.995 0.995 1 (0.7) 3.405 (0.201–57.78) 0.396 0.660 0 NA 0.996 0.996
 Dominant (CC vs. CG+GG) 0.758 (0.371–1.546) 0.446 0.535 0.505 (0.201–1.267) 0.145 0.394 1.389 (0.546–3.535) 0.491 0.779
 Recessive (CC+CG vs. GG) NA 0.995 0.995 3.438 (0.203–58.176) 0.392 0.960 NA 0.996 0.998
DROSHA rs10719 T>C
 TT 228 (56.6) 106 (53.0) 1.000 (reference) 78 (52.3) 1.000 (reference) 29 (53.7) 1.000 (reference)
 TC 158 (39.2) 82 (41.0) 1.129 (0.777–1.640) 0.524 0.786 60 (40.3) 1.089 (0.720–1.647) 0.686 0.686 18 (33.3) 0.865 (0.462–1.620) 0.650 0.726
 CC 17 (4.2) 12 (6.0) 1.646 (0.717–3.778) 0.240 0.600 11 (7.4) 1.786 (0.759–4.202) 0.184 0.623 7 (13.0) 3.451 (1.264–9.422) 0.016 0.096
 Dominant (TT vs. TC+CC) 1.170 (0.816–1.678) 0.392 0.535 1.156 (0.777–1.718) 0.475 0.535 1.086 (0.610–1.935) 0.779 0.779
 Recessive (TT+TC vs. CC) 1.588 (0.708–3.562) 0.262 0.655 1.674 (0.730–3.838) 0.224 0.960 3.499 (1.348–9.082) 0.010 0.060
RAN rs14035 C>T
 CC 240 (59.9) 119 (59.5) 1.000 (reference) 89 (59.7) 1.000 (reference) 39 (72.2) 1.000 (reference)
 CT 149 (37.0) 74 (37.0) 0.990 (0.680–1.442) 0.958 0.958 54 (36.2) 0.852 (0.560–1.296) 0.454 0.545 13 (24.1) 0.556 (0.283–1.092) 0.088 0.360
 TT 14 (3.5) 7 (3.5) 0.886 (0.327–2.397) 0.811 0.995 6 (4.0) 1.189 (0.417–3.391) 0.747 0.934 2 (3.7) 0.848 (0.180–3.985) 0.834 0.996
 Dominant (CC vs. CT+TT) 0.986 (0.684–1.421) 0.938 0.938 0.879 (0.585–1.320) 0.535 0.535 0.572 (0.301–1.088) 0.088 0.264
 Recessive (CC+CT vs. TT) 0.881 (0.334–2.321) 0.798 0.995 1.341 (0.480–3.742) 0.576 0.960 0.999 (0.217–4.585) 0.998 0.998
XPO5 rs11077 A>C
 AA 319 (79.2) 171 (85.5) 1.000 (reference) 125 (83.9) 1.000 (reference) 48 (88.9) 1.000 (reference)
 AC 79 (19.6) 29 (14.5) 0.665 (0.407–1.087) 0.103 0.378 24 (16.1) 0.756 (0.446–1.280) 0.297 0.465 6 (11.1) 0.490 (0.199–1.204) 0.120 0.360
 CC 5 (1.2) 0 NA 0.995 0.995 0 NA 0.993 0.934 0 NA 0.995 0.996
 Dominant (AA vs. AC+CC) 0.619 (0.380–1.008) 0.054 0.324 0.708 (0.420–1.196) 0.197 0.394 0.454 (0.185–1.112) 0.084 0.264
 Recessive (AA+AC vs. CC) NA 0.995 0.995 NA 0.993 0.993 NA 0.996 0.998

Values are presented as number (%).

RAN, Ran GTPase; XPO5, exportin 5; LAD, large artery disease; AOR, adjusted odds ratio; CI, confidence interval; SVD, small vessel disease; CE, cardioembolism; NA, not available.

*

Odds ratios adjusted for age, sex, hypertension, diabetes mellitus, hyperlipidemia, and smoking status;

P-value calculated by logistics regression analysis;

P-value calculated by false discovery rate test.