Identification of Blood Let-7e-5p as a Biomarker for Ischemic Stroke

Abstract

Circulating microRNAs (miRNAs) are emerging as novel disease biomarkers. Using a miRNA microarray, we previously showed that the whole blood level of let-7e-5p was significantly higher in ischemic stroke patients than in control subjects. However, the association between let-7e-5p expression and the occurrence of ischemic stroke remains unknown. In this study, we validated the expression levels of let-7e-5p in two case-control populations using miRNA TaqMan assays and further investigated the potential targets of let-7e-5p. The results suggest that the blood level of let-7e-5p was significantly higher in patients with ischemic stroke than in controls (p<0.05). Higher levels of let-7e-5p were associated with increased occurrence of ischemic stroke (adjusted OR, 1.89; 95% CI, 1.61~2.21, p<0.001) in the combined population. The addition of let-7e-5p to traditional risk factors led to an improvement in the area under the curve, which increased from 0.74 (95% CI, 0.70~0.78) to 0.82 (95% CI, 0.78~0.85), with a net reclassification improvement of 16.76% (p<0.0001) and an integrated discrimination improvement of 0.10 (p<0.0001) for patients with ischemic stroke. Bioinformatics prediction and cell experiments suggested that the expression levels of four genes enriched in the MAPK signaling pathway were down-regulated by let-7e-5p transfection. Specifically, the expression levels of the genes CASP3 and NLK were significantly lower in ischemic stroke patients than in controls and were negatively correlated with let-7e-5p expression. In summary, our study suggests the potential use of blood let-7e-5p as a biomarker for ischemic stroke and indicates its involvement in the related pathomechanism.

About the Authors:

Suli Huang

Contributed equally to this work with: Suli Huang, Ziquan Lv

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Ziquan Lv

Contributed equally to this work with: Suli Huang, Ziquan Lv

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Yi Guo

Affiliation: Department of Neurology, People’s Hospital of Shenzen, Guangdong, China

Lu Li

Affiliation: Key Laboratory of Environment and Health, Ministry of Education and Ministry of Environmental Protection, and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China

Yanwei Zhang

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Li Zhou

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Binyao Yang

Affiliation: State Key Lab of Respiratory Disease, Institute for Chemical Carcinogenesis, School of Public Health, Guangzhou Medical University, Guangzhou, Guangdong, China

Shuang Wu

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Ying Zhang

Affiliation: Department of Neurology, People’s Hospital of Shenzen, Guangdong, China

Changhui Xie

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Shanshan Li

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Jinquan Cheng

* E-mail: [email protected]

Affiliation: Shenzhen Center for Disease Control and Prevention, Shenzhen, Guangdong, China

Introduction

Globally, stroke is the second leading cause of death in people more than 60 years old, and it is the fifth leading cause of death in people between 15 and 59 years old [1]. In China, as in other countries, ischemic stroke is the most common type of stroke, accounting for 43% to 79% of the total number of strokes [2]. Ischemic stroke has several etiological origins, the most common one of which is the narrowing of the neck or head arteries through atherosclerosis. Atherosclerosis is a chronic disease process involving endothelial dysfunction, lipid disturbance, platelet activation, thrombosis, and inflammation. Multiple risk factors have also been identified for ischemic stroke, including age, male sex, smoking, history of hypertension and diabetes mellitus. However, a diagnosis biomarker for ischemic stroke has not been found. In recent years, advances in epigenomics have promoted the generation of many novel disease diagnosis biomarkers, among which microRNA (miRNA) have attracted extensive interest.

MiRNAs are short, endogenous, and non-coding RNAs that regulate gene expression at the posttranscriptional level by binding to 3-untranslated regions of their target mRNAs [3,4]. miRNAs have emerged as important regulators and fine-tuners of a range of pathophysiological cellular effects and molecular signaling pathways involved in atherosclerosis [5]. Moreover, miRNAs are also involved in the regulation of stroke-related cellular and molecular networks, and the identification of these miRNAs may provide new insights into the underlying mechanisms of stroke. For example, miRNAs are involved in the modulation of cerebral edema [6,7], a common consequence of stroke that contributes to tissue and cellular damage following ischemia. Microarray studies in animal models have also shown that numerous miRNAs are rapidly altered following ischemic preconditioning, which leads to the modulation of several key signaling pathways and an increase in cytoprotection and cellular regeneration in the affected brain tissue [8,9]. Several studies have begun to analyze miRNAs in plasma or serum to identify biomarkers for disease risk prediction or diagnosis [10–12]. However, identifying biomarkers using plasma or serum miRNAs has limitations because of the extremely low concentration of circulating miRNAs. Seeking to circumvent this problem, we hypothesized that miRNAs from peripheral whole blood could also be used. Peripheral blood contains immune cells, which have close contact with the vascular wall, a key component of the inflammatory response and are involved in the development and progression of atherosclerosis. Hoekstra et al. reported the dysregulation of several miRNAs in peripheral blood mononuclear cells from patients with coronary artery disease [13]. Another study revealed that 157 miRNAs were differentially expressed in the whole blood of stroke patients and controls [14]. Ontological analysis predicted that the targets of deregulated miRNAs were involved in angiogenesis, hypoxia, endothelial cell regulation, and immune response. Limited by a small sample size, however, the study discussed above did not determine the association between miRNA expression and stroke. In addition, the functional role of circulating miRNAs in ischemic stroke remains unknown.

In this study, we aimed to determine whether the whole blood level of let-7e-5p is differentially expressed in patients with ischemic stroke. We previously explored the miRNA expression profiles of 25 ischemic stroke patients and 25 control subjects using miRNA microarray and found that 455 miRNAs were differentially expressed (unpublished data, p<0.05). Of the differentially expressed miRNAs, let-7e-5p became our focus after reviewing the literature. For example, Li et al. found that the plasma level of let-7e-5p was significantly higher in essential hypertension patients than in control subjects [15], and let-7e-5p was also involved in the response of macrophages to lipopolysaccharide [16]. Because of the important role of let-7e-5p in the cardiovascular system, we selected let-7e-5p for further investigation and analyzed its potential as a biomarker for ischemic stroke.

Materials and Methods

Study design and populations

We determined the blood expression level of let-7e-5p in two independent case-control populations. Ischemic stroke patients from the stage I population were recruited from Shenzhen No. 2 Hospital between April 2013 and December 2013. Patients from the stage II population were consecutively recruited from Shenzhen People’s Hospital between June 2012 and December 2014. All subjects were inpatients diagnosed with ischemic stroke for the first time and were recruited based on the appearance of a new and abrupt focal neurological deficit, with neurological symptoms and signs persisting for more than 24 h, as described previously [17]. Ischemic stroke was confirmed by positive findings on head CT or MRI according to the International Classification of Disease (9th Revision, codes 430 to 438). Patients with a history of stroke, peripheral arterial occlusive disease or cancer were excluded from this study. Subjects without medical history of cerebrovascular diseases were selected as controls during a physical examination at the hospital and were matched with patients by age and sex. Approximately 5 mL of vein blood samples were collected from each participant in EDTA-anticoagulant tubes and stored at −80°C until use. Blood samples from patients with ischemic stroke were collected within 12 h after hospital admission.

For all participants, structured questionnaires were used by trained interviewers to collect information about demographic characteristics and clinical biochemistry. The ethics committee of the Shenzhen Center for Disease Control and Prevention approved this study, and written informed consent was obtained from each participant.

Total RNA isolation

Approximately 200 μL of whole blood was used from each sample to isolate total RNAs using the mirVana PARIS miRNA Isolation Kit (Ambion 1556, Austin, TX) according to the manufacturer’s protocol but with slight modifications, namely, each sample was treated twice with acid-phenol chloroform. The total RNAs of cultured cells were isolated using Trizol reagent (Invitrogen, AL, USA). To detect the expression levels of genes, the total RNA of each sample was isolated from 500 μL of whole blood using Trizol LS reagent (Invitrogen, AL, USA) according to the manufacturer’s protocol.

qRT-PCR assays

For miRNA detection, the input RNAs were reverse transcribed (RT) in a small-scale reaction using the TaqMan miRNA Reverse Transcription Kit (Applied BioSystems, Foster City, CA) following the manufacturer’s protocol (5 μL total volume with 1 μL of input RNA; components other than the input RNA were prepared as a large-volume master mix). RT products were diluted 1:5 and subjected to qPCR in triplicate using the TaqMan miRNA Assay Kit (Applied BioSystems, Foster City, CA) according to the manufacturer’s protocol in a small-scale reaction (10 μL total volume with 4.5 μL of diluted RT products; components other than input RT products were prepared as a large-volume master mix). MiRNA expression levels were normalized to U6 and calculated using the equation 2^-ΔCt, where ΔCt = cycle threshold (Ct) (miRNA)−Ct _(U6).

For the detection of gene expression, first strand cDNA was synthesized from equal amounts of RNA using the Super Script ^™ First-Strand Synthesis System (Life Technologies, USA). Quantitative PCR was performed in triplicate with the 7500 Real-time PCR System (Applied Biosystems, USA) using SYBR Premix Ex Taq^™ (Takara, Dalian, China). Gene expression values were calculated using the comparative quantitative method (the -ΔΔCT method) and normalized to values obtained from the amplification of beta-actin. PCR amplification was performed using the sets of primers designed by Primer 5.0, and the abbreviation list for the gene names were listed in S1 Table.

Target gene prediction and pathway analysis

The miRNA target genes were predicted using three bioinformatics tools: TargetScan (http://www.targetscan.org/), miRanda (http://www.microrna.org/microrna/home.do) and Pictar (http://pictar.mdc-berlin.de/). The pathway enrichment of predicted genes was analyzed using the DAVID tool (http://david.abcc.ncifcrf.gov/).

Cell culture and transfection

The U937 cell line, purchased from the Procell company (Wuhan, China), was cultured in RPMI-1640 medium (GIBCO) supplemented with 10% fetal calf serum (GIBCO) and penicillin/streptomycin in a humidified atmosphere containing 5% CO₂ at 37°C. Approximately 50 nM let-7e-5p mimic or control mimic purchased from RIBObio (Guangzhou RIBObio company) was transfected into U937 cells for 48 h using the riboFECT CP transfection kit (Guangzhou RIBObio company) according to the manufacturer’s protocol. All experiments were conducted with cells at a logarithmic stage of growth.

Statistical analysis

The normal distribution of data was tested using the 1-sample Kolmogorov–Smirnov test. Continuous variables were expressed as the mean ± SD or median (25th–75th quartile), and categorical variables were expressed as frequency. Differences in clinical characteristics between cases and controls were examined by the χ² test for categorical variables, Student’s t test for normally distributed data, or the Mann–Whitney U test for skewed data. The difference in circulating let-7e-5p expression levels or mRNA levels of predicted genes between cases and controls was examined using Student’s t test. Differences in mRNA levels of genes among the transfected cells were analyzed by one-way ANOVA. The odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to assess the association between let-7e-5p and ischemic stroke using the multivariate logistic regression model after adjusting for conventional risk factors, including age, sex, smoking, drinking, hypertension and diabetes mellitus. The ability of let-7e-5p to discriminate patients with ischemic stroke from controls or to reclassify patients with ischemic stroke was determined by the receiver operating characteristic curve and reclassification analysis using 2 models: a baseline model including conventional risk factors, and a model that incorporated the factors from the baseline model and let-7e-5p expression. Correlations between miRNA and gene expression levels in the blood were analyzed using the Pearson correlation test. Correlations between miRNA expression and platelet parameters were analyzed using the Spearman correlation test. All statistical analyses were performed using SPSS 11.0 software (Statistical Package for the Social Sciences, Chicago, IL). A value of p<0.05 was considered significant (2-tailed).

Results

General characteristics of the study population

The general characteristics of the two case-control populations are shown in Table 1. For the stage I population, the levels of high-density lipoprotein cholesterol (HDL-c), low-density lipoprotein cholesterol (LDL-c), and platelet distribution width (PDW) were significantly lower than in the controls (p<0.05). The levels of mean platelet volume (MPV), smoking frequency, drinking, and history of hypertension or diabetes mellitus were significantly higher in patients (p<0.05). The other variables showed no significant difference.

[Figure omitted. See PDF.]

Table 1. General characteristics of the study population.

https://doi.org/10.1371/journal.pone.0163951.t001

For the second population, the levels of fasting glucose, MPV, and history of hypertension and diabetes mellitus were significantly higher in patients than in control subjects (p<0.05), whereas the level of HDL-c was significantly lower in patients (p<0.05). The other variables showed no differences between the two groups.

Association of the whole blood level of let-7e-5p with the occurrence of ischemic stroke

We detected the whole blood level of let-7e-5p in the first case-control population. As shown in Fig 1, the results suggested that the expression level of let-7e-5p was significantly higher in ischemic stroke patients than in control subjects (Fig 1A, log-transformed expression levels relative to U6, 3.92± 0.74 vs. 3.14±0.98, p<0.001). In the second case-control population, we found that the expression level of let-7e-5p remained significantly higher in ischemic stroke patients than in control subjects (Fig 1B, log-transformed expression levels relative to U6, 4.99±1.04 vs. 4.07±1.37, p<0.001).

[Figure omitted. See PDF.]

Fig 1. The expression levels of let-7e-5p in stage I and II populations.

Let-7e-5p expression was detected in (A) stage I (n = 44 patients and controls) and (B) stage II (n = 302 patients and controls) population. The relative expression levels were normalized to U6 and then log-transformed. The whiskers of the plots represent the 2.5 to 97.5 percentiles.

https://doi.org/10.1371/journal.pone.0163951.g001

Pooled analysis was performed by combining the two populations. The results showed that higher levels of let-7e-5p were associated with an increased risk of ischemic stroke (adjusted OR, 1.89; 95% CI, 1.61~2.21; P<0.001; Table 2) after adjusting for conventional risk factors, including age, sex, smoking, drinking, history of hypertension and diabetes mellitus.

[Figure omitted. See PDF.]

Table 2. Association of the expression level of let-7e-5p with ischemic stroke.

https://doi.org/10.1371/journal.pone.0163951.t002

Diagnostic value of let-7e-5p by receiver operating characteristic curve and reclassification analysis

The diagnostic value of let-7e-5p was evaluated in the combined population. As shown in Table 3, the area under the receiver operating characteristic curve was 0.74 (95% CI, 0.70~0.78) for the baseline model and significantly increased to 0.82 (95% CI, 0.78~0.85) when let-7e-5p was added to the prediction model (p<0.05). The net reclassification improvement and integrated discrimination improvement were computed to determine the ability of miRNAs to reclassify patients misclassified by the baseline model. let-7e-5p was able to reclassify a significant portion of patients, with a net reclassification improvement (NRI) of 16.76% (95% CI, 11.16%~22.37%; p<0.0001). The integrated discrimination improvement (IDI) for let-7e-5p was 0.10 (95% CI, 0.08~0.13; p<0.0001).

[Figure omitted. See PDF.]

Table 3. AUC, NRI, and IDI calculations for let-7e-5p in the combined population.

https://doi.org/10.1371/journal.pone.0163951.t003

Correlation between let-7e-5p expression and platelet parameters

Platelet parameters play an important role in the pathogenesis of vascular diseases [18]. We further evaluated the correlation between let-7e-5p expression and platelet parameters in patients from the combined population, assessing platelet count value (PLT), MPV, PDW and plateletcrit (PCT). The results suggested that let-7e-5p was negatively correlated with PLT (r = -0.20, p<0.001), MPV (r = -0.192, p = 0.001) and PCT (r = -0.278, p<0.001). Conversely, let-7e-5p expression was positively correlated with PDW (r = 0.235, p<0.001), as shown in Fig 2.

[Figure omitted. See PDF.]

Fig 2. Correlations between let-7e-5p expression and platelet parameters.

The correlation was evaluated by the Spearman test in the combined patient population.

https://doi.org/10.1371/journal.pone.0163951.g002

Identification of let-7e-5p target genes

Using bioinformatics tools, we found that a total of 1800 genes were predicted to be targets of let-7e-5p. These genes were input into DAVID (http://david.abcc.ncifcrf.gov/) to obtain information about the signaling pathways in which the target genes may enriched. As shown in Table 4, the results suggested that 23 genes were enriched in the MAPK signaling pathway (p<0.001), 13 genes were enriched in the p53 signaling pathway (p = 0.017), 13 genes were enriched in the adherens junction (p = 0.042), and 14 genes were enriched in the TGF-beta signaling pathway (p = 0.047). Because most of the genes were enriched in the MAPK signaling pathway and because the relationship between let-7e-5p and the MAPK signaling pathway remains largely unknown, we selected these 23 genes for further investigation.

[Figure omitted. See PDF.]

Table 4. The enriched signal transduction pathways of predicted genes by let-7e-5p.

https://doi.org/10.1371/journal.pone.0163951.t004

To investigate the regulation of genes by let-7e-5p, we transfected the miRNA mimics into U937 cells and detected the changes in gene expression. As shown in Fig 3, the mRNA levels of the ETS-domain protein (ELK4), caspase 3 (CASP3), nemo-like kinase (NLK), and tumor protein p53 (TP53) were all significantly lower in the miRNA-transfected cells than in the control group. No difference was found between the mimic control and control groups for these four genes (p>0.05). The expression levels of the other 19 genes showed no significant differences among the three groups of treated cells (data not shown) (p>0.05).

[Figure omitted. See PDF.]

Fig 3. The mRNA levels of the target genes in different cell groups.

U937 cells were transfected with 50 nM control mimics (negative control) or 50 nM let-7e-5p mimics. The normal control is normal cultured cells. *p<0.05 compared to the negative control; **p<0.01 compared to the negative control.

https://doi.org/10.1371/journal.pone.0163951.g003

Correlations between let-7e-5p and target gene expressions

To further investigate the regulation of genes by let-7e-5p, we detected the expression levels of let-7e-5p and the 23 genes in the whole blood of 20 ischemic stroke patients and 20 control subjects from the second study population. The expression level of let-7e-5p was significantly higher in patients than in control subjects (p = 0.042, shown in S1 Fig), whereas the mRNA levels of 12 genes were found to be significantly lower in the patient group than in the control group (Fig 4, p<0.05). These genes were activating transcription factor 2 (ATF2), CASP3, cell division cycle 25B (CDC25B), ELK4, fibroblast growth factor receptor 2 (FGFR2), mitogen-activated protein kinase kinase kinase kinase 4 (MAP4K4), nerve growth factor (NGF), NLK, protein tyrosine phosphatase, receptor type R (PTPRR), RAS guanyl releasing protein 1 (RASGRP1), transforming growth factor beta receptor I (TGFBR1) and TP53. The correlation analysis conducted with the entire population further suggested that let-7e-5p expression was negatively correlated with ATF2 (r = -0.322, p = 0.049), CASP3 (r = -0.417, p = 0.020), FGFR2 (r = -0.427, p = 0.015), NLK (r = -0.527, p = 0.003), PTPN7 (r = -0.425, p = 0.009), RASGRP1 (r = -0.515, p = 0.004) and TGFBR1 (r = -0.496, p = 0.003). The other genes showed no significant correlation with let-7e-5p expression, as shown in Table 5.

[Figure omitted. See PDF.]

Table 5. Correlation analysis of let-7e-5p and gene expression.

https://doi.org/10.1371/journal.pone.0163951.t005

[Figure omitted. See PDF.]

Fig 4. The whole blood mRNA levels of the target genes in control subjects (n = 20) and ischemic stroke patients (n = 20).

The relative expression levels were normalized to beta-actin. *p<0.05; **p<0.01.

https://doi.org/10.1371/journal.pone.0163951.g004

Discussion

In this study, we found that the whole blood level of let-7e-5p was significantly higher in patients with ischemic stroke than in control patients in a Chinese population, and a higher level of let-7e-5p was associated with increased incidence of ischemic stroke. ROC curve and reclassification analyses suggested that the addition of let-7e-5p to traditional risk factors improved the predictive value of ischemic stroke. Bioinformatics prediction and target gene expression analysis showed that let-7e-5p may be involved in the pathogenesis of ischemic stroke by regulating CASP3 and NLK expression, two genes enriched in the MAPK signaling pathway. Collectively, our study demonstrated that let-7e-5p may act as a novel biomarker for ischemic stroke, the biological function of which warrants further investigation.

A biomarker should be non-invasive, easy to identify, and have good sensitivity and specificity. Plasma/serum miRNAs have been extensively studied and have provided a series of excellent biomarkers for the clinical diagnosis or prognosis of disease. However, these biomarkers also face technical limitations because of the extremely low concentrations of circulating miRNAs in the plasma or serum, and some miRNAs with low expression are not detectable. In addition, the absence of ribosomal RNAs in these samples makes assessment of RNA quality difficult. This is also an obstacle to the standardization and translation of the use of plasma/serum miRNAs in a clinical context. Seeking to circumvent this problem, we hypothesized that miRNAs from peripheral blood could be used. Differentially expressed miRNAs in the whole blood of patients with coronary heart disease [19,20] or peripheral arterial disease [21] have been documented. In this study, we demonstrated that the whole blood level of let-7e-5p was significantly higher in ischemic stroke patients and was independently associated with an increased risk of the disease after adjusting for conventional risk factors. The diagnosis potential of let-7e-5p was then evaluated by ROC and reclassification analyses. The addition of let-7e-5p to the traditional risk factor model increased the AUC from 0.74 to 0.82 (p<0.05). Furthermore, let-7e-5p was able to reclassify a significant portion of patients, with an NRI of 16.76% and IDI of 0.10; these results indicate good diagnostic values. Moreover, the correlation between let-7e-5p expression and the platelet parameters indicates the involvement of let-7e-5p in dysfunctional platelets. Our study was the first to demonstrate that the expression level of let-7e-5p in whole blood may serve as a useful noninvasive circulating biomarker for ischemic stroke. While we were preparing this manuscript, a study was published reporting that the expression levels of let-7e-5p in the serum and cerebral spinal fluid of acute ischemic stroke patients were also significantly higher than in the control subjects and showed a specificity as high as 73.4% and a sensitivity of 82.8% in ischemic stroke patients at the acute stage [22]. In addition, it has been suggested that the peripheral blood level of let-7e-5p exhibits good diagnostic value for peripheral arterial disease [21]. In addition, the plasma level of let-7e-5p was higher in hypertension patients than in control subjects [15,23]. In summary, the limited number of clinical studies about let-7e-5p have provided evidence that in addition to being a biomarker for ischemic stroke, let-7e-5p may also act as a molecular biomarker for other cardiovascular diseases. These findings suggest the involvement of let-7e-5p in the cardiovascular disease-related pathomechanism.

Few studies have reported on the biological function of miRNA let-7e-5p. One study demonstrated that let-7e may protect PC12 cells, a nerve cell line, against apoptosis following anoxia/reoxygenation injury by negatively regulating the expression ofCASP3[24]. The involvement of let-7e-5p in immunity has been documented by Quinn et al., who reported that let-7e-5p could regulate the TLR4 signaling pathway during the development of endotoxin tolerance at the receptor, signalling pathway, and gene transcription and translational levels [25]. In addition, Androulidaki et al. also demonstrated that let-7e-5p could repress TLR4, which is a critical protein for lipopolysaccharide-driven TLR signaling, and could regulate endotoxin sensitivity and tolerance [16]. Collectively, these studies only provide limited information about the functional role of let-7e-5p in cardiovascular diseases. In this study, we predicted target genes of let-7e-5p using bioinformatics tools and then tested the miRNA-gene regulation using cell transfection experiments and gene expression analysis of human whole blood samples. We found that two genes, CASP3 and NLK, were likely regulated by let-7e-5p. In cells transfected with let-7e-5p mimics, the mRNA levels of CASP3 and NLK were both significantly lower than in negative control. In addition, in human whole blood samples, the mRNA levels of both CASP3 and NLK were lower in ischemic stroke patients than in control subjects, and they were both negatively correlated with let-7e-5p expression. Our results were consistent with those of a previous study showing the regulation relationship between let-7e-5p and CASP3 [24]; however, this is the first report of the regulation of the gene NLK by let-7e-5p. NLK is an evolutionarily conserved MAP kinase-related kinase and plays a role in multiple processes due to its capacity to regulate a diverse array of signaling pathways, including the Wnt/β-catenin, Activin, IL-6, and Notch signaling pathways [26]. A previous study also demonstrated that NLK was involved in neuronal apoptosis after traumatic brain injury [27]. Therefore, we propose that the differentially expressed let-7e-5p may be functionally relevant to the deregulated apoptosis state by targeting genes CASP3 and NLK. Further investigation is required to confirm the target binding sites by luciferase reporter gene experiments and to evaluate the biological consequences of the artificial disturbance of let-7e-5p expression or the expression of CASP3 and NLK genes.

This study has several strengths. First, the expression levels of let-7e-5p were replicated in two independent case-control studies, and the large sample size provides sufficient statistical power for multivariate analysis. Second, we tested miRNA-gene regulation by conducting cell transfection experiments and miRNA-gene expression correlation analysis in human blood samples, the results of which provide a plausible biological explanation for the epidemiological findings in our study. However, the limitations should also be addressed. First, because this is a case-control study, potential selection bias may influence the interpretation of the results. Second, because circulating miRNAs are often correlated with each other and the combined miRNA panels may provide better sensitivity and specificity for diagnosis, as mentioned by previous studies [28,29], more miRNAs need to be detected to identify miRNA profiles with better prediction efficiency. Third, further investigations are needed to locate the binding sites between let-7e-5p and the target genes identified in this study. The biological consequences for the cerebrovascular diseases following changes in miRNA or target gene expression also merit further investigation. Finally, the causal relationship between the change of let-7e-5p expression and the occurrence of stroke also merits further investigation in a prospective cohort study.

In conclusion, this study suggests that the blood level of let-7e-5p is significantly higher in ischemic stroke patients and is associated with the occurrence of ischemic stroke. Moreover, let-7e-5p regulates the expression of CASP3 and NLK and may be involved in the pathogenesis of ischemic stroke. Further investigations are needed to explore the biological relevance of our findings.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. The expression levels of let-7e-5p in ischemic stroke patients (n = 20) and control subjects (n = 20).

The relative expression levels were normalized to U6 and then log-transformed. The data are expressed as mean±SD.

https://doi.org/10.1371/journal.pone.0163951.s001

(DOC)

S1 Table. Real-time PCR primer sequences for target genes.

https://doi.org/10.1371/journal.pone.0163951.s002

(DOC)

Acknowledgments

We are particularly grateful to all the ischemic stroke patients and volunteers who participated in this study and to the medical personnel of People’s Hospital and No. 2 Hospital of Shenzhen for their kind assistance in collecting questionnaires and blood samples.

Author Contributions

1. Conceptualization: SLH JQC.

2. Formal analysis: SLH ZQL.

3. Funding acquisition: SLH ZQL JQC.

4. Investigation: SLH ZQL CHX SSL.

5. Project administration: JQC.

6. Resources: LL YWZ BYY SW YZ JQC LZ YG.

7. Supervision: JQC.

8. Validation: SLH ZQL.

9. Writing – original draft: SLH ZQL.

10. Writing – review & editing: LL BYY LZ YG.

Citation: Huang S, Lv Z, Guo Y, Li L, Zhang Y, Zhou L, et al. (2016) Identification of Blood Let-7e-5p as a Biomarker for Ischemic Stroke. PLoS ONE 11(10): e0163951. https://doi.org/10.1371/journal.pone.0163951

References

1. Johnston SC, Mendis S, Mathers CD. Global variation in stroke burden and mortality: estimates from monitoring, surveillance, and modelling. Lancet Neurol. 2009; 8: 345–354. pmid:19233730

2. Liu L, Wang D, Wong KS, Wang Y. Stroke and stroke care in China: huge burden, significant workload, and a national priority. Stroke. 2011; 42: 3651–3654. pmid:22052510

3. Ambros V. The functions of animal microRNAs. Nature. 2004; 431: 350–355. pmid:15372042

4. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004; 116: 281–297. pmid:14744438

5. Feinberg MW, Moore KJ. MicroRNA Regulation of Atherosclerosis. Circ Res. 2016; 118: 703–720. pmid:26892968

6. Manley GT, Fujimura M, Ma T, Noshita N, Filiz F, Bollen AW, et al. Aquaporin-4 deletion in mice reduces brain edema after acute water intoxication and ischemic stroke. Nat Med. 2000; 6: 159–163. pmid:10655103

7. 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. pmid:18258830

8. Dharap A, Vemuganti R. Ischemic pre-conditioning alters cerebral microRNAs that are upstream to neuroprotective signaling pathways. J Neurochem. 2010; 113: 1685–1691. pmid:20402966

9. Lee ST, Chu K, Jung KH, Yoon HJ, Jeon D, Kang KM, et al. MicroRNAs induced during ischemic preconditioning. Stroke. 2010; 41: 1646–1651. pmid:20576953

10. Fichtlscherer S, De Rosa S, Fox H, Schwietz T, Fischer A, Liebetrau C, et al. Circulating microRNAs in patients with coronary artery disease. Circ Res. 2010; 107: 677–684. pmid:20595655

11. Huang S, Chen M, Li L, He M, Hu D, Zhang X, et al. Circulating MicroRNAs and the occurrence of acute myocardial infarction in Chinese populations. Circ Cardiovasc Genet. 2014; 7: 189–198. pmid:24627568

12. Zampetaki A, Kiechl S, Drozdov I, Willeit P, Mayr U, Prokopi M, et al. Plasma microRNA profiling reveals loss of endothelial miR-126 and other microRNAs in type 2 diabetes. Circ Res. 2010; 107: 810–817. pmid:20651284

13. Hoekstra M, van der Lans CA, Halvorsen B, Gullestad L, Kuiper J, Aukrust P, et al. The peripheral blood mononuclear cell microRNA signature of coronary artery disease. Biochem Biophys Res Commun. 2010; 394: 792–797. pmid:20230787

14. Tan KS, Armugam A, Sepramaniam S, Lim KY, Setyowati KD, Wang CW, et al. Expression profile of MicroRNAs in young stroke patients. PLoS One. 2009; 4: e7689. pmid:19888324

15. Li S, Zhu J, Zhang W, Chen Y, Zhang K, Popescu LM, et al. Signature microRNA expression profile of essential hypertension and its novel link to human cytomegalovirus infection. Circulation. 2011; 124: 175–184. pmid:21690488

16. Androulidaki A, Iliopoulos D, Arranz A, Doxaki C, Schworer S, Zacharioudaki V, et al. The kinase Akt1 controls macrophage response to lipopolysaccharide by regulating microRNAs. Immunity. 2009; 31: 220–231. pmid:19699171

17. Zhang Y, Tong Y, Ding H, Zhang H, Geng Y, Zhang R, et al. Two novel susceptibility SNPs for ischemic stroke using exome sequencing in Chinese Han population. Mol Neurobiol. 2014; 49: 852–862. pmid:24122314

18. Ross R. The pathogenesis of atherosclerosis: a perspective for the 1990s. Nature. 1993; 362: 801–809. pmid:8479518

19. Sinnaeve PR, Donahue MP, Grass P, Seo D, Vonderscher J, Chibout SD, et al. Gene expression patterns in peripheral blood correlate with the extent of coronary artery disease. PLoS One. 2009; 4: e7037. pmid:19750006

20. Grayson BL, Wang L, Aune TM. Peripheral blood gene expression profiles in metabolic syndrome, coronary artery disease and type 2 diabetes. Genes Immun. 2011; 12: 341–351. pmid:21368773

21. Stather PW, Sylvius N, Wild JB, Choke E, Sayers RD, Bown MJ. Differential microRNA expression profiles in peripheral arterial disease. Circ Cardiovasc Genet. 2013; 6: 490–497. pmid:24129592

22. Peng G, Yuan Y, Wu S, He F, Hu Y, Luo B. MicroRNA let-7e Is a Potential Circulating Biomarker of Acute Stage Ischemic Stroke. Transl Stroke Res. 2015; 6: 437–445. pmid:26415639

23. Cengiz M, Karatas OF, Koparir E, Yavuzer S, Ali C, Yavuzer H, et al. Differential expression of hypertension-associated microRNAs in the plasma of patients with white coat hypertension. Medicine (Baltimore). 2015; 94: e693. pmid:25837765

24. Peng G, Yuan Y, He Q, Wu W, Luo BY. MicroRNA let-7e regulates the expression of caspase-3 during apoptosis of PC12 cells following anoxia/reoxygenation injury. Brain Res Bull. 2011; 86: 272–276. pmid:21827835

25. Quinn EM, Wang J, Redmond HP. The emerging role of microRNA in regulation of endotoxin tolerance. J Leukoc Biol. 2012; 91: 721–727. pmid:22389313

26. Ishitani T, Ishitani S. Nemo-like kinase, a multifaceted cell signaling regulator. Cell Signal. 2013; 25: 190–197. pmid:23000342

27. Li Z, Cui G, Wang J, Yu Z, Zhao L, Lv Z. Nemo-like kinase (NLK) involves in neuronal apoptosis after traumatic brain injury. Cell Mol Neurobiol. 2012; 32: 381–389. pmid:22127415

28. Hu Z, Chen X, Zhao Y, Tian T, Jin G, Shu Y, et al. Serum microRNA signatures identified in a genome-wide serum microRNA expression profiling predict survival of non-small-cell lung cancer. J Clin Oncol. 2010; 28: 1721–1726. pmid:20194856

29. Zampetaki A, Willeit P, Tilling L, Drozdov I, Prokopi M, Renard JM, et al. Prospective Study on Circulating MicroRNAs and Risk of Myocardial Infarction. J Am Coll Cardiol. 2012; 60: 290–299. pmid:22813605