Gui-Jiang Wei1,2
Ming-Qing Yuan1Li-He Jiang1Yu-Lan Lu3Chun-Hong Liu3Hong-Cheng Luo3Hua-Tuo Huang3Zong-Quan Qi1*Ye-Sheng Wei1,2*- 1Department of Cell Biology, Medical College of Guangxi University, Nanning, China
- 2Department of Medical Laboratory, Affiliated Hospital of Guilin Medical University, Guilin, China
- 3Department of Medical Laboratory, Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
miRNAs are small non-coding RNAs modulating gene expression, and variants in miRNA genes are involved in the pathogenesis of ischemic stroke (IS). However, the effect of miR-34a polymorphisms on IS susceptibility has rarely been reported. In the present study, we investigated the association between rs12128240, rs2666433, and rs6577555 of the miR-34a gene and IS susceptibility. Snapshot assay was used to detect miR-34a polymorphisms in 548 IS patients and 560 controls. Relative expression of miR-34a was measured by quantitative real-time PCR. We found that rs2666433 was associated with a significantly increased risk of IS (AA vs. GG: OR = 1.61, 95% CI = 1.05–2.52, P = 0.031; AA vs. GG+GA: OR = 1.58, 95% CI = 1.05–2.45, P = 0.026). For the IS subtypes, rs2666433 was associated with large artery atherosclerosis (AA vs. GG: OR = 2.09, 95% CI = 1.16–3.51, P = 0.007; AA vs. GG+GA: OR = 2.02, 95% CI = 1.15–3.33, P = 0.007; A vs. G: OR = 1.36, 95% CI = 1.07–1.81, P = 0.021). Additionally, the level of miR-34a was significantly up-regulated in IS patients compared to the controls (P < 0.001), and patients with rs2666433 AA genotype had a higher level of miR-34a than those with GG+GA genotypes (P < 0.001). Furthermore, increased level of homocysteine was observed in IS patients compared to the controls (P < 0.001), especially in patients carrying the rs2666433AA genotype compared to those carrying the rs2666433 GG+GA genotypes (P < 0.001). However, no significant association between rs12128240 or rs6577555 and IS was found. Collectively, our study found the association between miR-34a polymorphisms and the risk of IS among the Chinese population. The results may provide an explanation for etiology of IS and a potential biomarker or therapeutic target for IS.
HIGHLIGHTS
- MiR-34a rs2666433 polymorphism was associated with an increased risk of ischemic stroke.
-The level of miR-34a was significantly up-regulated in ischemic stroke patients compared with controls, and patients with rs2666433 AA genotype had a higher level miR-34a than those with GG+GA genotypes.
-Furthermore, increased level of homocysteine was showed in IS patients compared to controls, and in patients carrying the rs2666433AA compared to those carrying the rs2666433 GG+GA.
Introduction
Stroke is a leading cause of mortality among the elderly population in the world (Adeloye, 2014; Kim, 2014; Robert and Zamzami, 2014). It can be divided into hemorrhagic and ischemic strokes (IS), with the latter being the most common form of stroke, accounting for 42–79% of all strokes (Zhang et al., 2003). In China, because of its high incidence, IS has exceeded the incidence of heart diseases and cancers, bringing great economic and health burdens to families and society (Liu et al., 2011). Atherosclerosis is an inflammatory reaction that involves aggregation of leukocytes, vascular smooth muscle cell and endothelial senescence (Libby, 2002). The etiology of atherosclerosis is complex but generally involves the formation of atherosclerotic plaques. Evidence showed cardiovascular obstruction causes coronary artery disease and cerebrovascular obstruction leads to IS (Weber and Noels, 2011). Multiple factors, including diabetes mellitus, hypertension, smoking, hyperhomocysteinemia and hyperlipidemia, are associated with a higher risk of IS (Tuhrim, 2000). Nevertheless, these risk factors just clarify some parts of IS pathogenesis. Clinical and epidemiological studies suggest genetic factors play an important role in IS pathogenesis (Jerrard et al., 2003).
MiRNAs are small non-coding RNAs which bind to the 3′-untranslated region (3′-UTR) of target messenger RNAs (mRNAs), and then promote mRNAs degradation or inhibit translation (Bhaskaran and Mohan, 2014). Many studies have shown that they play important roles in various diseases such as neurodegenerative disorders (Karnati et al., 2015), cancer (Wang et al., 2017), cardiovascular diseases (Srinivasan and Das, 2015) and stroke (Khoshnam et al., 2017). The miR-34a gene, located at the chromosome 1p36 locus, is a modulator of the vascular smooth muscle cell and endothelial senescence (Badi et al., 2015; Ito et al., 2010). The up-regulated expression of miR-34a was observed in both atherosclerosis and IS (Raitoharju et al., 2011; Liang and Lou, 2016), and inhibition of miR-34a prevents endothelial cell apoptosis in the setting of atherosclerosis (Li et al., 2018). These outcomes suggest that miR-34a is a potential IS biomarker or a new target for IS therapy.
Single nucleotide polymorphisms (SNPs) are the most common variations at specific base positions of the human genome. The expression of miRNAs can be affected by SNPs locating in the coding genes of miRNAs, and subsequently change a variety of biological functions (Kim et al., 2016). Previous works have shown that the susceptibility of individuals to multiple human diseases may be modulated by SNPs of the miR-34a gene such as osteosarcoma (Lv et al., 2014), colon cancer (Jiang et al., 2017), type 2 diabetes (Sun et al., 2018) and IS (Choi et al., 2016). However, the association between miR-34a SNPs and IS susceptibility has not been reported in the Chinese population. To examine the association between miR-34a SNPs and IS susceptibility, we performed a case-control study. Three SNPs of the miR-34a gene, namely rs12128240, rs2666433 and rs6577555, were selected based on location, minor allele frequency (MAF) and previous reports. Meanwhile, we also assessed the impact of miR-34a polymorphisms on the expression of miR-34a and other IS-related factors.
Materials and Methods
Study Population
Five hundred and forty-eight patients in this study were recruited from the Department of Neurology, Affiliated Hospital of Youjiang Medical University for Nationalities, Guangxi, China in the period from January 2017 to December 2017. Based on clinical symptoms, neurological examination, magnetic resonance imaging or/and computed tomography, an IS patient was diagnosed by two experienced clinical neurologists. According to Trial of ORG 10172 in Acute Stroke Treatment (TOAST), only patients with subtypes of small artery occlusion (SAO) and large artery atherosclerosis (LAA) were included in this study (Adams et al., 1993). Exclusion criteria were listed below: recurrent IS, hemorrhagic stroke, tumorous, traumatic, infectious and cardiogenic cerebrovascular diseases, cerebrovascular malformation and severe mental disorder. To check the neurological status of the patients, National Institutes of Health Stroke Scale (NIHSS) on admission was used to scale IS severity (Tilley et al., 1996). To determine the relationship between the miR-34a gene polymorphisms and IS severity, IS patients were divided into two subgroups according to the NIHSS score (<6 and ≥6) (O’Donnell et al., 2010).
Five hundred and sixty gender and age-matched controls were recruited from the Health Examination Department of Affiliated Hospital of Youjiang Medical University for Nationalities, Guangxi, China in the period from January 2017 to December 2017. Exclusion criteria were listed below: stroke history, cancer, informatory, genetic, autoimmune, neurological or cardiogenic diseases. All the patients and controls came from the Guangxi Zhuang Autonomous Region of China.
Clinical information for patients and controls comprising age, gender, hypertension, diabetes mellitus, smoking and stroke history were recorded on admission or receiving health examination, respectively. Hypertension was defined as receiving anti-hypertensive treatment currently or having a mean diastolic blood pressure ≥ 90 mm Hg or/and a systolic blood pressure ≥ 140 mm Hg. Smokers referred to those smoking daily. Blood samples of controls were collected during health examination, and those of patients were collected within 24 h after onset of symptoms on admission. Total cholesterol (TCH), homocysteine, triglyceride, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), very low-density lipoprotein cholesterol (VLDL-C) and high sensitivity C-reactive protein (hsCRP) were determined.
The study was approved by the Ethics Committee of Youjiang Medical University for Nationalities and conformed to the guidelines set forth by the Declaration of Helsinki. Meanwhile, written informed consent was obtained from all patients and controls prior to entering the study.
SNPs Selection and Genotyping
The criteria to select the SNPs in this study include: (I) Minor allele frequency is greater than 10% in the southern Han Chinese population; (II) The SNPs should locate in the regulatory regions of the miR-34a gene; (III) The SNPs have been reported by previous references. According to previous reports (Kim et al., 2012), 555 SNPs were found by searching upstream 2 kb and downstream 500 b of the miR-34a gene in the National Center for Biotechnology Information SNP database1. Three SNPs (rs2666433/NC_000001.11:g.9153118A>G and rs6577555/NC_000001.11:g.9152228A>C in the promoter; rs12128240/NC_000001.11:g.9151244C>T in the 3′-flanking region) with MAF greater than 10% in southern Han Chinese population were selected.
DNA extraction kits (Yaneng, China) were used to extract genomic DNA from whole blood samples and the genomic DNA was stored at –80°C. Primers for genotyping SNPs of the miR-34a gene were designed by Primer Premier 6 and synthesized by Shanghai Sangon Corporation. Primers for genotyping SNPs of the miR-34a gene are listed in Supplementary Table 1. After amplifying the miR-34a gene by polymerase chain reaction, the rs12128240, rs2666433 and rs6577555 polymorphisms were genotyped by Snapshot assay. Finally, Sanger DNA sequencing was used for verifying genotyping results of Snapshot assay.
Quantitative Analysis of miR-34a
The analysis was conducted on a subgroup of 79 IS patients and 79 controls, who were chosen according to the results of genotyping. Among the 79 IS patients, 39 patients had AA genotype and 40 patients had GG+GA genotypes. There was a similar balance in the controls. Total RNA was extracted from 200 μL plasma of controls and patients following the instruction of RNA extraction kit (Qiagen, Germany). Synthetic cel-miR-39 was used for normalization of sample variation based on the previous studies (Hauser et al., 2012; Wang et al., 2018), and validated in our previous study. The concentration of RNA was determined by Nanodrop 2000 Spectrophotometer (Thermo Fisher Scientific, United States). Reverse transcription kit (Ribobio, China) was used to transcribe 100 ng of total RNA into complementary DNA. The products of reverse transcription were stored at –80°C. Quantitative real-time PCR was performed by Mir-X miRNA qRT-PCR SYBR Kit (Takara, China), on Applied Biosystems 7900HT system (ABI, United States) and followed these steps: 95°C for 10 s, 40 cycles of 95°C for 5 s and 60°C for 20 s. According to the result of quantitative real-time PCR for miR-148b conducted previously in our laboratory, miR-148b was chosen as a negative control in this study. The primers for miR-34a/39/148b were designed and synthesized by Ribobio Corporation. Reactions were performed in triplicate. Relative expression levels of miR-34a and miR-148b were calculated by 2Ct(cel-miR-39)-Ct(miR-34a) and 2Ct(cel-miR-39)-Ct(miR-148b), respectively.
Statistical Analysis
Statistical analysis was performed by Statistical Package for Social Science (version 17.0). Hardy-Weinberg equilibrium testing and Categorical variables comparison were proceeded by two-sided chi-square test. The normality of continuous variables was assessed by the Shapiro–Wilk test. Continuous variables were analyzed by the Mann–Whitney U test or the Student’s t-test. Logistic regression was used for evaluating the association between miR-34a polymorphisms and IS risk, and odds ratio (OR), 95% confidence interval (CI) and P-values were adjusted by gender, age, smoker, diabetes and hypertension. The Benjamini–Hochberg(B-H) method was applied to control the false discovery rate. Genotypes and alleles distribution of rs2666433 in different populations were obtained from the website below2. Online bioinformatics software was used to predict the potential targets of miR-34a3. An online tool SHEsis50 was used for Haplotype analysis (SNPs order: rs12128240-rs2666433-rs6577555). P < 0.05 indicated statistical significance. The power calculation was performed by the GAS Power Calculator. Because the MAF of rs6577555 was lower than rs12128240 and rs2666433, the power calculation was performed according to rs6577555. Under the parameters of disease prevalence = 0.016 (Jiang et al., 2017), MAF = 0.20, significance level = 0.05, OR = 1.5 and sample size (cases = 548, controls = 560), we had 95% statistical power to identify a convincing association between SNPs and IS risk.
Results
Clinical Features
The clinical features of controls and IS patients are shown in Table 1. Smoking, hypertension, diabetes and high levels of TCH, homocysteine, hsCRP, LDL-C were more common in the IS patients than the controls (P < 0.001). In contrast, the level of HDL-C in the IS patients was significantly lower than that in the controls (P < 0.001). Furthermore, as for gender, age, triglyceride and VLDL-C, there were no significant differences between IS patients and the controls (P = 0.739, P = 0.284, P = 0.432 and P = 0.498, respectively). As shown in Supplementary Table 3, as for gender, age and smoker, there were no significant differences between the NHISS < 6 patients and the NHISS ≥ 6 patients (P = 0.541, P = 0.192 and P = 0.130, respectively). However, as for hypertension and diabetes mellitus, there were significant differences between the NHISS < 6 patients and the NHISS ≥ 6 patients (P = 0.011 and P = 0.008, respectively).
Table 1. Clinical characteristics of the study population.
Association Between the miR-34a Polymorphisms and the Risk of IS
Sequencing diagrams of rs12128240, rs2666433, and rs6577555 are shown in Supplementary Figure 1. The association between the rs12128240, rs2666433 or rs6577555 polymorphism and IS risk is shown in Table 2. The genotypes of rs12128240, rs2666433, and rs6577555 in both controls (P = 0.434, P = 0.992, and P = 0.975, respectively) and cases (P = 0.491, P = 0.318 and P = 0.866, respectively) were in Hardy-Weinberg equilibrium. We found that the rs2666433 AA genotype was associated with a significantly increased IS risk compared with the GG genotype (AA vs. GG: OR = 1.61, 95% CI = 1.05–2.52, P = 0.031). Similarly, we also observed a significantly increased risk of IS in AA vs. GG+GA model analysis (AA vs. GG+GA: OR = 1.58, 95% CI = 1.05–2.45, P = 0.026). However, the statistical significance was lost after correcting the P-values by B-H method. Additionally, no significant association between the rs12128240 or rs6577555 and the risk of IS was observed.
Table 2. Association between the miR-34a polymorphisms and risk of IS.
Association Between the miR-34a Polymorphisms and LAA or SAO Subtypes of IS
The association between the miR-34a polymorphisms and LAA or SAO subtypes of IS was investigated according to the TOAST classification. As shown in Table 3, AA genotype of rs2666433 was associated with a significantly increased risk of IS LAA subtype compared with GG genotype, and this association remained significant in AA vs. GG+GA model analysis (AA vs. GG+GA: OR = 2.02, 95% CI = 1.15–3.33, P = 0.007). Notably, the statistical significance remained after correcting the P-values by B-H method. Additionally, A allele of rs2666433 was also associated with a significantly increased risk of IS LAA subtype compared with G allele (AA vs. GG: OR = 2.09, 95% CI = 1.16–3.51, P = 0.007; A vs. G: OR = 1.36, 95% CI = 1.07–1.81, P = 0.021). In contrast, the rs2666433 was not associated with SAO. Meanwhile, the rs12128240 and rs6577555 were not associated with LAA and SAO. The clinical features of LAA and SAO subgroups are shown in Supplementary Table 2.
Table 3. Association of miR-34a polymorphisms with LAA and SAO subtypes of IS.
Association Between the miR-34a Polymorphisms and the NIHSS Score
To analyze the association between the miR-34a polymorphisms and IS severity, we further investigated the association between the rs12128240, rs2666433 or rs6577555 and the NIHSS score on admission. As shown in Supplementary Table 3, of the 548 patients, 312 individuals were classified as having a mild stroke (NHISS < 6) and 236 as having a severe stroke (NHISS ≥ 6). However, as shown in Table 4, we did not find a significant association between the miR-34a polymorphisms and IS severity (P > 0.05).
Table 4. Association between the miR-34a polymorphisms and the NIHSS score.
Linkage Disequilibrium and Haplotype Analysis
Linkage disequilibrium and haplotype analysis were also performed. As shown in Table 5, seven possible haplotypes were listed, the rs12128240 was in strong linkage disequilibrium with the rs2666433 (D′ = 0.904) and the rs6577555 (D′ = 0.846), and the rs2666433 was also in strong linkage disequilibrium with the rs6577555 (D′ = 0.801). Moreover, CGC and CAC were identified as the two main haplotypes in both controls (21.5 and 24.3%, respectively) and cases (37.6 and 36.6%, respectively). Nevertheless, we did not find a significant association between miR-34a haplotypes and the risk of IS.
Table 5. Haplotype analysis of the miR-34a polymorphisms with risk of IS.
Association Between the rs2666433 Polymorphism and the Expression Level of miR-34a
The plasma levels of miR-34a and miR-148b were examined among the IS patients and the controls. The average RNA integrity number among samples was 8.6. As shown in Figure 1, the plasma level of miR-34a was significantly up-regulated in IS patients compared with the controls (P < 0.001). Particularly, we found that patients with AA genotype of rs2666433 had a higher level of miR-34a than those with GG+GA genotypes (P < 0.001). Controls with AA genotype of rs2666433 also had a higher level of miR-34a than those with GG+GA genotypes, but the difference was not significant (P = 0.067). As a negative control, the plasma level of miR-148b showed no significant differences among IS patients, controls and different genotype groups (P > 0.05). The clinical characteristics and genotypes distribution of qPCR subgroups are shown in Supplementary Table 4.
Figure 1. Association between rs2666433 polymorphism and expression level of miR-34a. The effect of rs2666433 polymorphism on the expression of miR-34a in controls and IS patients. (A) Increased level of miR-34a in IS patients compared to controls (P < 0.001), especially in IS patients carrying the rs2666433 AA compared with those carrying the rs2666433 GG+GA (P < 0.001). (B) As a negative control, plasma levels of miR-148b showed no significant differences among IS patients, controls and different genotype groups (P > 0.05). Plasma levels of miRNAs were normalized to cel-miR-39. The upper and lower boundaries of the boxes represent the 25th and 75th percentiles of the expression value. The upper and lower whiskers represent the maximum and minimum values. The horizontal line within the box represents the median value. All P-values were adjusted by gender, age, smoker, diabetes and hypertension.
Association Between the rs2666433 Polymorphism and the Serum Level of Homocysteine
We further analyzed the association between rs2666433 polymorphism and serum level of homocysteine. As shown in Figure 2B, the serum level of homocysteine was significantly up-regulated in IS patients compared with the controls (P < 0.001). Additionally, after correlation analysis of the rs2666433 polymorphism and the serum level of homocysteine, we found that patients with the rs2666433 AA genotype had a higher level of homocysteine than those with GG+GA genotypes (P < 0.001). Controls with AA genotype of rs2666433 also had a higher level of homocysteine than those with GG+GA genotypes, but the difference was not significant (P = 0.098). As shown in Figure 2A, after predicting the potential targets of miR-34a by software Targetscan, Methylenetetrahydrofolate reductase (MTHFR) was found to be a potential target of miR-34a. The relevant information of bioinformatics analysis is shown in Supplementary Table 5.
Figure 2. Association between rs2666433 polymorphism and serum level of homocysteine. The effect of rs2666433 polymorphism on the serum level of homocysteine in controls and IS patients. (A) Bioinformatics software predicted that MTHFR gene was a target of miR-34a. (B) Increased level of homocysteine in IS patients compared to controls (P < 0.001), especially in IS patients carrying the rs2666433AA compared to those carrying the rs2666433 GG+GA (P < 0.001). The upper and lower boundaries of the boxes represent the 25th and 75th percentiles of the homocysteine value. The upper and lower whiskers represent the maximum and minimum values. The horizontal line within the box represents the median value. MTHFR: methylenetetrahydrofolate reductase; UTR: untranslated region.
Genotypes and Alleles Distribution of rs2666433 in Different Populations
Considering the importance of rs2666433 polymorphism in the pathogenesis of IS, we further compared the genotypes and alleles distribution of rs2666433 in 11 different populations. As shown in Table 6, the genotypes and alleles distribution of rs2666433 in the present study were significantly different from 1KGP-ESN, 1KGP-ASW, 1KGP-CLM, 1KGP-PUR, 1KGP-TSI, 1KGP-CEU, and 1KGP-PJL populations (P < 0.05). Nevertheless, we did not find significant difference after comparing with 1KGP-CHS, 1KGP-CHB, and 1KGP- JPT populations (P > 0.05).
Table 6. Genotypes and alleles distribution of rs2666433 in different populations.
Discussion
To the best of our knowledge, this is the first study which investigated the association between rs12128240, rs2666433, and rs6577555 of the miR-34a gene and IS susceptibility. The main findings of the present study are listed as below: (I) MiR-34a rs2666433 polymorphism was associated with an increased risk of IS. (II) The level of miR-34a was significantly up-regulated in IS patients compared with the controls, and patients with rs2666433 AA genotype had a higher level of miR-34a than those with GG+GA genotypes. (III) Increased level of homocysteine was observed in IS patients compared to the controls, especially in patients carrying the rs2666433AA genotype compared to those carrying the rs2666433 GG+GA genotypes. Taken together, these findings indicate that the rs2666433 polymorphism may contribute to the susceptibility of IS.
In China, there are millions of people suffering from IS, which has brought great economic and health burdens to many families and society (Liu et al., 2011). Clinical and epidemiological studies suggested genetic factors played an important role in IS pathogenesis (Jerrard et al., 2003). As previously reported, miR-34a, a tumor suppressor miRNA in many species, suppresses proliferation and metastasis of tumor cells by silencing target mRNAs (Li R. et al., 2015; Li T. et al., 2015; Liang et al., 2015). In 2011, analysis of miRNAs expression profile performed by Raitoharju et al. (2011) revealed that miR-34a is upregulated in human atherosclerotic plaques. These results suggested the possibility of using miR-34a as a biomarker in atherosclerosis-related diseases such as IS. Subsequently, research conducted in 2015 showed overexpression of miRNA-34a in blood samples of acute IS patients and a rat model of middle cerebral occlusion compared to the controls (Liang and Lou, 2016). In contrast, inhibition of miR-34a prevented endothelial cell apoptosis in the setting of atherosclerosis (Li et al., 2018). In this study, overexpression of miR-34a in IS patients was also observed, validating the results of previous studies. Taking these results together, miR-34a may be related to the etiology of IS, but not every individual with up-regulation of miR-34a eventually suffered from IS, indicating that genetic variants like SNPs may also be involved in the etiology of IS.
Single nucleotide polymorphisms located in miRNA have been widely reported as the cause of abnormal functions of the miRNAs. Previously, some studies have reported that SNPs located in the miR-34a gene were associated with susceptibility to multiple human diseases such as osteosarcoma (Lv et al., 2014), colon cancer (Jiang et al., 2017), type 2 diabetes (Sun et al., 2018) and IS (Choi et al., 2016). To date, few studies focused on the association between miR-34a gene polymorphisms and IS susceptibility. In this study, we failed to identify a positive association between the rs6577555 polymorphism and IS or IS subtypes. Consistent with our data, Choi et al. (2016) also could not find any rs6577555 polymorphism frequency differences between IS and control groups in a Korean population. It may be ascribed to a similar distribution of SNPs among different East Asian populations. Table 6 supported our assumption, since the genotypes and alleles distribution of rs2666433 in the Guangxi population was found to be similar to those in 1KGP-CHS, 1KGP-CHB and 1KGP- JPT populations. So as far as we know, this is the first report of the association between the rs12128240 or rs2666433 and IS risk. In the present study, it was found that the rs12128240 polymorphism is not related to IS or IS subtypes. There is one possible explanation for this result. Rs12128240 is a sequence variant locating in 3′ flanking region of the miR-34a gene, which lacks important gene regulation related elements. Promisingly, we found a positive association between the rs2666433 polymorphism and IS or IS subtypes. AA genotype of rs2666433 was associated with a significantly increased risk of IS and LAA subtype compared with GG genotype. Moreover, A allele was also associated with a significantly increased risk of LAA subtype compared with G allele. These results indicated that rs2666433 may be involved in the etiology of IS.
More and more studies revealed that SNPs of the gene promoter can modulate the promoter transcription activity and the gene expression level by affecting binding efficiency for transcriptional factors (Hamaoka et al., 2013; Ma et al., 2016). Additionally, using bioinformatics analysis, the rs2666433 polymorphism in the promoter may affect the binding of transcription factors (Sun et al., 2018). In consideration of the evidence, we deduced that rs2666433 may influence the expression level of miR-34a. Subsequently, we detected the expression level of miR-34a and assessed miR-34a expression with rs2666433 polymorphism. In the current study, the plasma level of miR-34a was significantly up-regulated in IS patients compared to controls. Additionally, patients with rs2666433 AA genotype had a higher level of miR-34a than those with GG+GA genotypes. These outcomes support our assumption.
Based on a number of previously published clinical evidence, a significantly increased homocysteine level has been regarded as a crucial role in the formation of atherosclerosis and IS (Sacco et al., 2004; Fang et al., 2014; Chen et al., 2015). MTHFR is a key enzyme which plays an important role in the metabolism of homocysteine. A lower expression level of MTHFR causes accumulation of homocysteine, and then results in a higher level of homocysteine (Moll and Varga, 2015; He W. et al., 2017; Liu et al., 2018). According to research conducted by Zhang et al. (2015), as a member of the miR-34 family, miR-34b may influence the serum level of homocysteine by binding the 3′-UTR of MTHFR mRNA. After predicting the potential targets of miR-34a by software Targetscan, we found that MTHFR was also a potential target of miR-34a. Inspired by the observations above, we investigated the association between miR-34a rs2666433 polymorphism and serum level of homocysteine. In the present study, the serum level of homocysteine was significantly up-regulated in IS patients, and patients with the rs2666433 AA genotype had a significantly higher level of homocysteine than those with GG+GA genotypes. All the results above suggested that miR-34a rs2666433 AA genotype might upregulate the level of miR-34a, which in turn downregulated the level of MTHFR, resulted in the accumulation of homocysteine and IS occurrence.
In this study, controls with AA genotype of rs2666433 also had higher levels of homocysteine and miR-34a than those with GG+GA genotypes, but the difference was not significant. There was an explanation for this phenomenon. As we know, the miRNA expression can be regulated by the polymorphism in the promoter of a gene. Conversely, the polymorphism in the promoter of a gene is not the only factor that could affect the miRNA expression. Many studies reported that the long non-coding RNAs including Lnc-MEG3 (Huang et al., 2017a) and Lnc-OC1 (Tao et al., 2018), circular RNAs including circ-GFRA1 (He R. et al., 2017) and circ-NFIX (Xu et al., 2018) and coding RNA including PDL1 (Huang et al., 2017b) and LDHA (Xiao et al., 2016), can function as competing endogenous RNAs to regulate miR-34a expression by sponging miR-34a. Under normal conditions, the normal expression of these genes does not affect the expression of miR-34a. However, under pathological conditions, down-regulation of these genes can further up-regulate the expression of miR-34a based on the up-regulatory effect of the miR-34a polymorphism. Therefore, a gene-gene interaction analysis in the future will better reveal the roles of miR-34a and its regulatory genes in the etiology of ischemic stroke.
In the present study, we also investigated the association between the miR-34a polymorphisms and IS severity. However, no significant association was found between the miR-34a polymorphisms and IS severity. Additionally, 7 possible haplotypes were found in our study, and CGC and CAC were the two main haplotypes in both controls and cases. Nevertheless, we did not find a significant association between miR-34a haplotypes and the risk of IS.
There are some limitations in the current study. First, our results are just representative for IS patients of Guangxi Chinese in southwest China, data from other patient cohorts is also needed to verify our results. Second, because of selection bias, the results may be influenced by control and case ascertainment. Third, it should be acknowledged that the differences in characteristics among qPCR subgroups may have a potential effect on the expression profile of miR-34a.
In summary, the present study found the association between miR-34a polymorphisms and the risk of IS in the Guangxi Chinese in southwest China. The results may provide some evidence for etiology of IS and a potential biomarker or therapeutic target for IS. In order to validate the results, further studies in ethnically disparate populations are going to be carried out.
Ethics Statement
The study got the approval by the Ethics Committee of Youjiang Medical College, and written informed consent was obtained from all patients and controls prior to entering the study.
Author Contributions
G-JW and M-QY wrote the manuscript. G-JW, M-QY, L-HJ, Y-LL, C-HL, H-CL, and H-TH collected samples and performed the experiments. Z-QQ and Y-SW conceived and designed the experiments.
Funding
This study was supported by National Natural Science Foundation of China (Nos. 81560552, 81260234), Natural Science Foundation of Guangxi Zhuang Autonomous Region (CN) (2017JJA180826), Innovation Project of Guangxi Graduate Education (CN) (201601009) and Key Laboratory Open Project Fund of Guangxi Zhuang Autonomous Region (CN) (kfkt20160064).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2019.00432/full#supplementary-material
Footnotes
- ^http://www.ncbi.nlm.nih.gov/SNP
- ^http://www.internationalgenome.org
- ^http://www.targetscan.org/mamm_31/
References
Adams, H. P., Bendixen, B. H., Kappelle, L. J., Biller, J., Love, B. B., Gordon, D. L., et al. (1993). Classification of subtype of acute ischemic stroke. Stroke 24, 35–41. doi: 10.1161/01.STR.24.1.35
Adeloye, D. (2014). An estimate of the incidence and prevalence of stroke in Africa: a systematic review and meta-analysis. PLoS One 9:e100724. doi: 10.1371/journal.pone.0100724
Badi, I., Burba, I., Ruggeri, C., Zeni, F., Bertolotti, M., Scopece, A., et al. (2015). MicroRNA-34a induces vascular smooth muscle cells senescence by sirt1 downregulation and promotes the expression of age-associated pro-inflammatory secretory factors. J. Gerontol. A Biol. Sci. Med. Sci. 70, 1304–1311. doi: 10.1093/gerona/glu180
Bhaskaran, M., and Mohan, M. (2014). MicroRNAs: history, biogenesis, and their evolving role in animal development and disease. Vet. Pathol. 51, 759–774. doi: 10.1177/0300985813502820
Chen, Y., Liu, R., Zhang, G., Yu, Q., Jia, M., Zheng, C., et al. (2015). Hypercysteinemia promotes atherosclerosis by reducing protein S-nitrosylation. Biomed. Pharmacother. 70, 253–259. doi: 10.1016/j.biopha.2015.01.030
Choi, G. H., Ko, K. H., Kim, J. O., Kim, J., Oh, S. H., Han, I. B., et al. (2016). Association of miR-34a, miR-130a, miR-150 and miR-155 polymorphisms with the risk of ischemic stroke. Int. J. Mol. Med. 38, 345–356. doi: 10.3892/ijmm.2016.2609
Fang, P., Zhang, D., Cheng, Z., Yan, C., Jiang, X., Kruger, W. D., et al. (2014). Hyperhomocysteinemia potentiates hyperglycemia-induced inflammatory monocyte differentiation and atherosclerosis. Diabetes Metab. Res. Rev. 63, 4275–4290. doi: 10.2337/db14-0809
Hamaoka, K., Nagoshi, S., Sugawara, K., Naiki, K., Uchida, Y., Inao, M., et al. (2013). SNPs in the promoter region of the osteopontin gene as a possible host factor for sex difference in hepatocellular carcinoma development in patients with HCV. Hepatol. Int. 7, 683–692. doi: 10.1007/s12072-012-9404-1
Hauser, S., Wulfken, L. M., Holdenrieder, S., Moritz, R., Ohlmann, C. H., Jung, V., et al. (2012). Analysis of serum microRNAs (miR-26a-2∗, miR-191, miR-337-3p and miR-378) as potential biomarkers in renal cell carcinoma. Cancer Epidemiol. 36, 391–394. doi: 10.1016/j.canep.2012.04.001
He, R., Liu, P., Xie, X., Zhou, Y., Liao, Q., Xiong, W., et al. (2017). circGFRA1 and GFRA1 act as ceRNAs in triple negative breast cancer by regulating miR-34a. J. Exp. Clin. Cancer Res. 36:145. doi: 10.1186/s13046-017-0614-1
He, W., Lu, M., Li, G., Sun, Z., Liu, D., and Gu, L. (2017). Methylene tetrahydrofolate reductase (MTHFR) rs868014 polymorphism regulated by mir-1203 associates with risk and short term outcome of ischemic stroke. Cell. Physiol. Biochem. 41, 701–710. doi: 10.1159/000458429
Huang, X., Gao, Y., Qin, J., and Lu, S. (2017a). The mechanism of long non-coding RNA MEG3 for hepatic ischemia-reperfusion: mediated by miR-34a/Nrf2 signaling pathway. J. Cell. Biochem. 119, 1163–1172. doi: 10.1002/jcb.26286
Huang, X., Xie, X., Hua, W., Xiao, X., Lu, Y., Zhi, T., et al. (2017b). PDL1 And LDHA act as ceRNAs in triple negative breast cancer by regulating miR-34a. J. Exp. Clin. Cancer Res. 36:129. doi: 10.1186/s13046-017-0593-2
Ito, T., Yagi, S., and Yamakuchi, M. (2010). MicroRNA-34a regulation of endothelial senescence. Biochem. Biophys. Res. Commun. 398, 735–740. doi: 10.1016/j.bbrc.2010.07.012
Jerrard, P., Cloud, G., Hassan, A., and Markus, H. S. (2003). Evaluating the genetic component of ischemic stroke subtypes: a family history study. Stroke 34, 1364–1369. doi: 10.1161/01.str.0000069723.17984.fd
Jiang, H., Ge, F., Hu, B., Wu, L., Yang, H., and Wang, H. (2017). rs35301225 polymorphism in miR-34a promotes development of human colon cancer by deregulation of 3′UTR in E2F1 in chinese population. Cancer Cell Int. 17:39. doi: 10.1186/s12935-017-0402-1
Karnati, H. K., Panigrahi, M. K., Gutti, R. K., Greig, N. H., and Tamargo, I. A. (2015). miRNAs: key players in neurodegenerative disorders and epilepsy. J. Alzheimers Dis. 48, 563–580. doi: 10.3233/jad-150395
Khoshnam, S. E., Winlow, W., Farbood, Y., Moghaddam, H. F., and Farzaneh, M. (2017). Emerging roles of micrornas in ischemic stroke: as possible therapeutic agents. J. Stroke 19, 166–187. doi: 10.5853/jos.2016.01368
Kim, J., Choi, G. H., Ko, K. H., Kim, J. O., Oh, S. H., Park, Y. S., et al. (2016). Association of the single nucleotide polymorphisms in micrornas 130b, 200b, and 495 with ischemic stroke susceptibility and post-stroke mortality. PLoS One 11:e0162519. doi: 10.1371/journal.pone.0162519
Kim, J. S. (2014). Stroke in Asia: a global disaster. Int. J. Stroke 9, 856–857. doi: 10.1111/ijs.12317
Kim, Y. O., Kim, S. Y., Yun, D. H., and Lee, S. W. (2012). Association between ABCB1 polymorphisms and ischemic stroke in korean population. Exp. Neurobiol. 21, 164–171. doi: 10.5607/en.2012.21.4.164
Li, R., Shi, X., Ling, F., Wang, C., Liu, J., Wang, W., et al. (2015). MiR-34a suppresses ovarian cancer proliferation and motility by targeting AXL. Tumour Biol. 36, 7277–7283. doi: 10.1007/s13277-015-3445-8
Li, T., Li, L., Li, D., Wang, S., and Sun, J. (2015). MiR-34a inhibits oral cancer progression partially by repression of interleukin-6-receptor. Int. J. Clin. Exp. Pathol. 8, 1364–1373.
Li, Y., Zhang, K., and Mao, W. (2018). Inhibition of miR34a prevents endothelial cell apoptosis by directly targeting HDAC1 in the setting of atherosclerosis. Mol. Med. Rep. 17, 4645–4650. doi: 10.3892/mmr.2018.8411
Liang, J., Li, Y., Daniels, G., Sfanos, K., De Marzo, A., Wei, J., et al. (2015). LEF1 targeting EMT in prostate cancer invasion is regulated by miR-34a. Mol. Cancer Res. 13, 681–688. doi: 10.1158/1541-7786.mcr-14-0503
Liang, T. Y., and Lou, J. Y. (2016). Increased expression of mir-34a-5p and clinical association in acute ischemic stroke patients and in a rat model. Med. Sci. Monit. 22, 2950–2955. doi: 10.12659/MSM.900237
Liu, L., Wang, D., Wong, K. S., and Wang, Y. (2011). Stroke and stroke care in China: huge burden, significant workload, and a national priority. Stroke 42, 3651–3654. doi: 10.1161/strokeaha.111.635755
Liu, X., Wang, L., Chi, H., Wang, J., Zheng, Q., Li, J., et al. (2018). The SNP Rs915014 in MTHFR regulated by MiRNA associates with atherosclerosis. Cell. Physiol. Biochem. 45, 1149–1155. doi: 10.1159/000487355
Lv, H., Pei, J., Liu, H., Wang, H., and Liu, J. (2014). A polymorphism site in the premiR34a coding region reduces miR34a expression and promotes osteosarcoma cell proliferation and migration. Mol. Med. Rep. 10, 2912–2916. doi: 10.3892/mmr.2014.2582
Ma, X., Jiao, X., Ma, Z., and Hejtmancik, J. F. (2016). Polymorphism rs7278468 is associated with age-related cataract through decreasing transcriptional activity of the CRYAA promoter. Sci. Rep. 6:23206. doi: 10.1038/srep23206
Moll, S., and Varga, E. A. (2015). Homocysteine and MTHFR mutations. Circulation 132, e6–e9. doi: 10.1161/circulationaha.114.013311
O’Donnell, M. J., Xavier, D., Liu, L., Zhang, H., Chin, S. L., Rao-Melacini, P., et al. (2010). Risk factors for ischaemic and intracerebral haemorrhagic stroke in 22 countries (the INTERSTROKE study): a case-control study. Lancet 376, 112–123. doi: 10.1016/s0140-6736(10)60834-3
Raitoharju, E., Lyytikainen, L. P., Levula, M., Oksala, N., Mennander, A., Tarkka, M., et al. (2011). miR-21, miR-210, miR-34a, and miR-146a/b are up-regulated in human atherosclerotic plaques in the tampere vascular study. Atherosclerosis 219, 211–217. doi: 10.1016/j.atherosclerosis.2011.07.020
Robert, A. A., and Zamzami, M. M. (2014). Stroke in Saudi Arabia: a review of the recent literature. Pan Afr. Med. J. 17:14. doi: 10.11604/pamj.2014.17.14.3015
Sacco, R. L., Anand, K., Lee, H. S., Boden, B., Stabler, S., Allen, R., et al. (2004). Homocysteine and the risk of ischemic stroke in a triethnic cohort: the northern manhattan study. Stroke 35, 2263–2269. doi: 10.1161/01.str.0000142374.33919.92
Srinivasan, H., and Das, S. (2015). Mitochondrial miRNA (MitomiR): a new player in cardiovascular health. Can. J. Physiol. Pharmacol. 93, 855–861. doi: 10.1139/cjpp-2014-0500
Sun, Y., Peng, R., Li, A., Zhang, L., Liu, H., Peng, H., et al. (2018). Sequence variation in microRNA-34a is associated with diabetes mellitus susceptibility in a southwest chinese han population. Int. J. Exp. Pathol. 11, 1637–1644.
Tao, F., Tian, X., Lu, M., and Zhang, Z. (2018). A novel lncRNA, Lnc-OC1, promotes ovarian cancer cell proliferation and migration by sponging miR-34a and miR-34c. J. Genet. Genomic 45, 137–145. doi: 10.1016/j.jgg.2018.03.001
Tilley, B. C., Marler, J., Geller, N. L., Lu, M., Legler, J., Brott, T., et al. (1996). Use of a global test for multiple outcomes in stroke trials with application to the national institute of neurological disorders and stroke t-pa stroke trial. Stroke 27, 2136–2142. doi: 10.1161/01.STR.27.11.2136
Wang, H., Zhi, H., Ma, D., and Li, T. (2017). MiR-217 promoted the proliferation and invasion of glioblastoma by repressing YWHAG. Cytokine 92, 93–102. doi: 10.1016/j.cyto.2016.12.013
Wang, R., Wang, C. F., Qin, H. M., Lu, Y. L., Wei, G. J., Huang, H. T., et al. (2018). Association between polymorphisms in the promoter region of miR-17-92 cluster and systemic lupus erythematosus in a chinese population. J. Cell. Mol. Med. 22, 4016–4020. doi: 10.1111/jcmm.13672
Weber, C., and Noels, H. (2011). Atherosclerosis: current pathogenesis and therapeutic options. Nat. Med. 17, 1410–1422. doi: 10.1038/nm.2538
Xiao, X., Huang, X., Ye, F., Chen, B., Song, C., Wen, J., et al. (2016). The miR-34a-LDHA axis regulates glucose metabolism and tumor growth in breast cancer. Sci. Rep. 6:21735. doi: 10.1038/srep21735
Xu, H., Zhang, Y., Qi, L., Ding, L., Jiang, H., and Yu, H. (2018). NFIX circular RNA promotes glioma progression by regulating miR-34a-5p via notch signaling pathway. Front. Mol. Neurosci. 11:225. doi: 10.3389/fnmol.2018.00225
Zhang, L. F., Yang, J., Hong, Z., Yuan, G. G., Zhou, B. F., Zhao, L. C., et al. (2003). Proportion of different subtypes of stroke in China. Stroke 34, 2091–2096. doi: 10.1161/01.str.0000087149.42294.8c
Zhang, W., Lin, W. Q., Cao, H. F., Li, C. Y., and Li, F. (2015). Association of a miR-34b binding site single nucleotide polymorphism in the 3′-untranslated region of the methylenetetrahydrofolate reductase gene with susceptibility to male infertility. Genet. Mol. Res. 14, 12196–12204. doi: 10.4238/2015.October.9.8
Keywords: ischemic stroke, miR-34a, single nucleotide polymorphism, genotype, risk
Citation: Wei G-J, Yuan M-Q, Jiang L-H, Lu Y-L, Liu C-H, Luo H-C, Huang H-T, Qi Z-Q and Wei Y-S (2019) A Genetic Variant of miR-34a Contributes to Susceptibility of Ischemic Stroke Among Chinese Population. Front. Physiol. 10:432. doi: 10.3389/fphys.2019.00432
Received: 30 November 2018; Accepted: 28 March 2019;
Published: 24 April 2019.
Edited by:
Bernhard Schaller, University of Zurich, SwitzerlandReviewed by:
Michela Masè, Fondazione Bruno Kessler, ItalyKatarzyna Bogunia-Kubik, Ludwik Hirszfeld Institute of Immunology and Experimental Therapy (PAN), Poland
Copyright © 2019 Wei, Yuan, Jiang, Lu, Liu, Luo, Huang, Qi and Wei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zong-Quan Qi, 13627769496@163.com Ye-Sheng Wei, yeshengwei66@yeah.net