Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 24 May 2019
Sec. Computational Genomics
This article is part of the Research Topic Machine Learning Techniques on Gene Function Prediction View all 48 articles

A Systems Analysis of the Relationships Between Anemia and Ischemic Stroke Rehabilitation Based on RNA-Seq Data

\r\nYingying Wang&#x;Yingying Wang1†Xingxian Huang&#x;Xingxian Huang2†Jianfeng LiuJianfeng Liu3Xuefei ZhaoXuefei Zhao4Haibo Yu*Haibo Yu2*Yunpeng Cai*Yunpeng Cai1*
  • 1Research Center for Biomedical Information Technology, Shenzhen Institutes of Advanced Technologies, Chinese Academy of Sciences, Shenzhen, China
  • 2Shenzhen Traditional Chinese Medicine Hospital, Shenzhen, China
  • 3Department of Neurology, The First Affiliated Hospital of Harbin Medical University, Harbin, China
  • 4Institute of Harbin Hematology & Oncology, The First Hospital of Harbin, Harbin, China

Ischemic stroke (IS) is one of the main causes of morbidity and disability worldwide due to its complex mechanism. Anemia was characterized as a risk factor of IS because the direct connection between central nervous system, blood supply, and tissue oxygen delivery. As the key oxygen-carrying molecule in the blood, hemoglobin (Hb) may be decisive in the destiny of penumbral area or influence the brain recovery and neurologic function, which could finally affect the outcome of IS. However, more detailed information on the expression levels of Hb related genes were still lacking possibly because the concentration of Hb was determined by the genes’ expression several hours ago, which may make the research more difficult to perform. This time gap between gene expressions and protein concentration could make these genes predictive bio-markers for IS outcome. In this study, we choose 28 IS patients, of which 12 were suffering from anemia. Statistical analysis results showed that the outcomes of the patients were different when dividing them into two groups characterized by Hb concentration. 2 sex and age matched patients were first chosen to perform RNA-seq analysis on, on two occasions at two different time points, after which the Hb counts were tested at least 24 h after sequencing. Results showed that the outcome of anemia patients was poor compared with non-anemia patients. Two other patients were then chosen for analysis which excluded the coincidence of other factors. The results showed that the low value of Hb under 13 g/dL in men were closely related to the poor outcome of IS patients. Differently expressed Hb related genes were tested and six genes were shown to be positively correlated with the recovery degree of IS patients: ELANE, FGF23, HBB, PIEZO1, RASA4, and PRTN3. Gene CPM was shown to be negatively correlated with clinical outcomes. All of the seven genes were validated to be related to strokes using real-time PCR or literature searches. Taken together, these genes could be considered as new predictors for the recovery of IS patients.

Introduction

Ischemic stroke (IS) is one of the main causes of morbidity and disability worldwide due to its complex mechanism. Many bio-markers and risk factors had been successfully identified using different methods (Jiang et al., 2018a,b, 2019). Similarly, alongside smoking, diabetes mellitus, hypertension, and hypercholesterolemia, anemia was characterized as ‘the fifth cardiovascular risk factor’ (Kaiafa et al., 2015). For example, the mortality rate was shown to be significantly higher in atherosclerosis-related IS patients suffering from anemia when admitted (Huang et al., 2009). The results of ASTRAL (Acute Stroke Registry and Analysis of Lausanne) showed that anemia on admission could predict both short-term and long-term outcomes in patients with IS and the risk of recurrent stroke was higher in anemia patients (Milionis et al., 2015).

The relationships between IS and anemia may be explained partly by the direct connection between central nervous system (CNS), blood supply, and tissue oxygen delivery (Dubyk et al., 2012). A sudden interruption of the oxygen supply to brain tissue was a crucial step in the pathophysiology of IS. The regaining of oxygen supply to ensure timely reperfusion or collateral perfusion would then determine the conditions of brain tissues (Kellert et al., 2011). As the key oxygen-carrying molecule in the blood, hemoglobin (Hb) may be decisive for the destiny of penumbral area or influence the brain recovery and neurologic function, which could finally affect the outcome of IS (Kimberly et al., 2013; Park et al., 2013). It had been shown that poor outcome and mortality after IS were strongly associated with low and further decreasing hemoglobin (Hb) and hematocrit (Hct) levels (Kellert et al., 2011). Another study found that the poor outcome of acute IS was related to the lower but not the higher end of the Hb, regardless of the time point and methods Hb concentrations been measured (Kimberly et al., 2013). It was widely accepted that women often had worse outcomes after suffering a stroke compared to men. A study focusing on the relationship between Hb and clinical outcome measured by mRS found that sex differences in stroke outcome were linked to lower Hb level, which was more prevalent in women (Park et al., 2013).

However, more detailed information on the expression levels of Hb related genes was still lacking which might be caused by the fact that Hb was in mature red blood cells with no nucleus. In other words, the concentration of Hb was determined by the genes’ expression several hours ago which may make the research more difficult to perform. However, this time gap between gene expressions and protein concentration could make these genes as predictive bio-markers for IS outcome. Compared with commonly used clinical scales such as NIH Stroke Scale (NIHSS) score, modified Rankin Scale (mRS), and biochemical tests such as routine blood tests that could describe the current recovery degree of patients, the bio-markers on gene levels could act as an early predictor for the IS patients which may help reduce the chance or severity of disability and chance of recurrence.

In this study, we chose 28 IS patients, of which 12 were suffering from anemia. Statistical analysis results showed that the outcome of the patients was different when dividing them into two groups characterized by Hb concentration. 2 sex and age matched patients were first chosen to perform RNA-seq analyses twice on different time points, after which the Hb counts were tested at least 24 h after sequencing. Results showed that the outcome of anemia patients was poor compared with the other patients. Then two other patients were chosen to exclude other factors besides Hbs. The results showed that the low value of Hb under 13 g/dL in men was closely related to the poor outcome of IS patients. Different expressed Hb related genes were tested and six genes were shown to be positively correlated with the recovery degree of IS patients: ELANE, FGF23, HBB, PIEZO1, RASA4, and PRTN3. Gene CPM was shown to be correlated with clinical outcome negatively. All of the seven genes were validated to be related to stroke using real-time PCR or literature search. Taken these together, these genes could be considered as new predictors for the recovery of IS patients.

Materials and Methods

The framework of this work could be divided into several steps shown in Figure 1, which will be elaborated on in the remainder of this section.

FIGURE 1
www.frontiersin.org

Figure 1. Framework of this work.

Clinical Samples

Whole blood samples were collected from 46 IS patients under different stages to perform analyses on three levels as follows:

(1) Clinical information analyses group (C-group): all the 28 IS patients in this group were in non-acute period (at least 10 days after the first onset occurred). They were divided into two sub-groups based on whether they were diagnosis as anemia according to the World Health Organization (WHO) criteria (Hb concentration < 13 g/dL in men and < 12 g/dl in women). The anemia ischemic stroke group (ISA) was composed of 12 patients while the non-anemia ischemic stroke group (non-ISA) contained 16 patients.

(2) RNA-Seq data analyses group (R-group): four patients in C-group were chosen including 3 non-ISA patients (marked as P1, P3, and P4) and 1 ISA patient (marked as P2). The basic clinical characters (including age, time intervals, and medical history) of P1 and P2 were similar. P3 and P4 were chosen to exclude the influence of initial condition (mRs, NIHSS, and myodynamia) when admitted.

(3) PCR analyses group (P-group): nine IS patients whose first onset occurred in May, 2018 and 9 control were chosen to perform real-time PCR for the validation of selected bio-markers. The patients were divided into two groups: 3 in ISA, 6 in non-ISA. 5 TIA patients and 4 healthy people were divided into control group.

Data Generation

Besides Hb concentration, different types of data sets were generated from the three clinical sample groups as follows:

(1) C-group: the NIH Stroke Scale (NIHSS) score, pre-stroke disability, and clinical outcome defined using the modified Rankin Scale (mRS), BI (Barthel Index), and Myodynamia (including four levels) of all the 28 patients were collected (see Supplementary Table S1 for details). The higher scores of NIHSS, mRS or Myodynamia indicated the worse outcome of IS patients while a lower score of BI indicated a worse outcome.

(2) R-group: Illumina Hiseq2500 V4 was used to perform the RNAseq analyses on the whole blood samples for each patients on two different time points. The raw sequence data reported in this paper has been deposited in the Genome Sequence Archive (Wang et al., 2017b) in Big Data Center Members (2018), Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, with accession number CRA001420 and is publicly accessible at http://bigd.big.ac.cn/gsa, as well as in the NCBI SRA repository with accession number PRJNA533989. The HB counts were tested and recorded at least 24 h after sequencing.

(3) P-group: the total RNA was extracted from blood samples using TRIzol LS according to the manufacturer’s protocol. The purified RNA was conveyed to cDNA with Golden 1st cDNA Synthesis Kit (with gDNA Remover). A SYBR Green fluorescent quantitative PCR Kit was used for real-time PCR analysis. The amplification program used in reactions were: 95°C for 15 min, 95°C for 10 s, 40 cycles of 30 s 60°C annealing.

The study was approved by SIAT Institutional Review Committee with IRB number SIAT-IRB-16515-H0107 and all the procedures were in accordance with the SIAT-IRB guidelines and the Declaration of Helsinki.

Statistical Analyses for C-Group

The baseline characteristics of the patients were compared between the ISA and non-ISA groups using the Student’s t- test for continuous variables and the Pearson χ2 tests for categorical variables. The categorical variables including sex and medical history for each group of patients were reported as percentages, and continuous variables were reported as mean ± SD.

Regression analyses were performed using linear models to explore the relationships between the Hb count and NIHSS, mRS, BI, and Myodynamia. A small Pr( > | t| ) indicated a significant correlation between Hb counts and the corresponding factors mentioned above.

Bioinformatics Analyses for R-Group

(1) Differential expression analyses: the sequencing raw datasets generated in R-group were analyzed using Tophat2 and Cufflinks. The differential expressed genes (DEG) between the two different time points for a same patient were identified using Cuffdiff with q-value not over 0.05. The biological DEG changes were further analyzed using Fold Change (FC) methods, which was the ratio between two time points based on the calculation of the average expression values. The common DEGs were selected as initial candidate bio-markers.

(2) Hemoglobin pattern analyses: the expression of 141 hemoglobin related genes retrieved by using ‘hemoglobin’ as key word in NCBI Gene with at least one expression over 0 in the 8 samples we collected were selected. A co-expression network was constructed by calculating the Pearson correlation coefficient for each pair of the 141 genes. All the pairs with p-value not over 0.05 were kept to construct the HB-network (hemoglobin network).

(3) Candidate HB-network (CHB-network) construction: the common DEGs were mapped into the HB-network. All the genes connecting to the common DEGs in the HB-network were kept to form the CHB-network.

(4) Extended CHB-network (ECHB-network) construction and analyses: all the genes in the CHB-network were mapped into the HB-network. Any genes connecting to any of them were kept to construct the ECHB-network. The topological analyses were performed using a R package “iGraph.” Each gene in the ECHB-network were considered as candidate bio-markers, an expression changed score (ECS) were calculated as follows:

ECS=(E2E1)/E1

Of which, E2 presented the expression value of the gene on time-point 2 while E1 presented the expression value of the gene on time-point 1. This score reflected the changes each genes occurs between the two time points that could help estimate the biological roles of different genes. DAVID was used to perform the functional analyses for these bio-markers.

Real-Time PCR Analysis for P-Group

The relative differences in expression between the four different groups were measured using ΔΔcycle time (CT) values: the CT values of the interested genes were semi-quantitative corrected using the internal reference gene β-actin. The expression of each gene was reported as mean ± SD. The Wilcoxon rank test was used to perform the statistical analysis.

Results and Discussion

Baseline Characteristics of Patients in C-Group

The baseline characteristics of the 28 patients in C-group were listed in Table 1. There were no significant differences (with p-value over 0.05) between the two groups on the characteristics of sex, lateral upper limb-distal, hypertension, diabetes mellitus, hyperlipemia, hyperuricemia, cerebral arteriosclerosis, atrial fibrillation, coronary disease, and carotid atherosclerosis. The regression analysis between Hb count and these baseline characters showed similar results: there were no significant associations between Hb counts and sex, Lateral upper limb-proximal, Lateral upper limb-distal, and all the Medical history besides pulmonary infection. It was interesting that the relationships between stroke, Hb counts, and pulmonary infection were different from other related diseases which might be explained by a ‘double hit model’ showing reasons of stroke-induced respiratory distress syndrome (Mascia, 2009).

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics.

Compared with this, the patients in ISA group were older than those in non-ISA group. Importantly, the higher NIHSS, mRs, and lower BI of ISA group compared with non-ISA showing the worse outcome of ISA patients. Besides, the myodynamia scores of the ISA group were lower than the non-ISA group. These results exhibited the important role of anemia played in the recovery of IS. The patients with lower HB counts were likely to have a worse outcome. The regression analysis results listed in Table 1 show the close relationship between Hb counts and the outcome of IS patients with smaller Pr( > | t| ) values.

Differential Expression Analyses in R-Group

Differential expression analyses were performed to find the biomarkers of IS based on the two whole blood samples extracted from same patients on different time points in R-group. Patients P1 in non-ISA and P2 in ISA were chosen for this purpose since their basic clinical characters were similar as follows:

(1) Age: both under 65 years old. Besides age-matched, this basic feature could exclude the results caused by aging which was shown to be an important factor for the IS outcome.

(2) Time intervals: the two blood extracting time-points of P1 were 32 and 56 days after the first onset of IS while for P2 were 34 and 54 days, correspondingly. The similar time intervals of the two patients could exclude the influence caused by different stages of IS which had been shown to be of great importance in the IS outcome (Wang et al., 2017a).

(3) Medical history: both P1 and P2 suffered from hypertension which was considered as an independent risk factor for the worse outcome of IS. However, the outcomes of the two patients were different: P1 was shown to have better outcome with smaller mRs, higher BI and myodynamia compared with P2 (Table 2). The minimum HB counts of P1 and P2 were different on both of the two time points which indicated us that this could a possible reason for the different outcomes.

TABLE 2
www.frontiersin.org

Table 2. Patients information.

Despite the above correlation, we still cannot make the conclusion that different outcomes were related with HB counts since the initial conditions of P1 and P2 were different: the NIHSS of P1 was 2 while P2 was 5 showing that the condition of P1 was much better than P2 when admitted. To exclude this influence, we chose two other patients (marked as P3 and P4) from the non-ISA group with higher NIHSS scores (9 and 10) respectively. Interestingly, we found similar results as follows:

(1) P3 was shown to have better outcome than P2 with lower mRs, higher BI and myodynamia when discharged. The minimum HB count of P3 was 142, higher than P2.

(2) P4 was similar to P2 on myodynamia evaluations. However, the smaller mRs and higher BI also showed better outcome of P4 than P2. Similarly as P1, P3, the minimum HB count of P4 was 136, higher than P2. The outcome of P4 may be a result of the patient’s age, as this patient was 86 years old, commonly considered by most to be very old.

Under the threshold of q-value ≤ 0.05, the numbers of DEGs were small: 0 for P1, 8 for P2, 48 for P3, and 30 for P4 between the two time points. The common DEGs between P2 and P3 were PRKCSH, HBD, and HBB. The common DEGs between P2 and P4 were DHPS, CFD, IRF9, FTH1P10, and HBB. HBB was the only common DEGs among P2, P3, and P4. The fold change calculating for HBB showed that for better outcome of P1 and P3, the fold change value was positive (2.42268 for P1 and 4.59503 for P3). Compared with this, the fold change value for worse outcomes was negative (−3.43035 for P2 and −3.35511 for P4). This indicated us that the decrease on HBB expression levels were closely related to the worse outcome of IS, especially in anemia or older IS patients.

Correlations Between Candidate Bio-Markers and IS Outcome

The HB-network was composed of 141 nodes (hemoglobin related genes) and 977 edges between them. HBB was found to be connected to nine genes as follows: PIEZO1, FGF23, ELANE, PRTN3, FN3KRP, TERT, HFE, TNF, and REN (see Figure 2 for details). These 10 genes formed the CHB-network. 6 of these genes were shown to be related to stroke as follows: (1) PIEZO1 was involved in hypertension-dependent arterial re-modeling (Retailleau et al., 2015). (2) FGF23 was a risk factor for overall stroke (Wright et al., 2014). (3) TERT was a determinant of risk of IS in the Atherosclerosis Risk in Communities (ARIC) study (Bressler et al., 2015). (4) HFE may play a role in modifying the relationship between smoking and stroke (Njajou et al., 2002). (5) TNF was differently expressed between stroke and controls (Bokhari et al., 2014). (6) REN was an important part in the renin angiotensin aldosterone system (RAAS), which played important roles in acute IS. Besides, the change of RAAS was shown to be related to the increased blood pressure (Back et al., 2015).

FIGURE 2
www.frontiersin.org

Figure 2. Network composed of 23 Hb related genes. Nine genes directly connected with HBB were shown by the black border lines. Different shapes represent the relationships between these genes and IS as follows: HBB-diamond. To be validated-rectangle; stroke biomarker-ellipse, stroke-related risk factor-octagon.

We further extended the CHB-network by finding out the neighbor of the 9 nodes and got other 14 nodes: AGT, CD59, CPM, ICAM1, MMP9, POU5F1, PRKCA, SMG5, UGT1A1, CSNK2A1, RASA4, GP1BA, TUBGCP3, and PKM (see Figure 2 for details). These 24 genes formed the ECHB-network. Their functional analysis results based on KEGG pathway using DAVID were shown in Supplementary Table S2. seven genes were shown to be related to stroke as follows: (1) AGT: its polymorphisms were proven to be related to the risk of IS in a meta-analysis in the Chinese population (Gao et al., 2015). (2) CD59: 50% patients with three additional mutations in CD59 were shown to have recurrent strokes (Tabib et al., 2017). (3) ICAM1: a neuro-inflammatory biomarker in post-stroke (Deddens et al., 2013). (4) MMP9: the combination therapies of MMP-9 inhibitor along with tPA was proven to be beneficial in IS (Chaturvedi and Kaczmarek, 2014). (5) PRKCA was related to the blood pressure (Shahin and Johnson, 2016). (6) RASA4 was involved in Ras signaling pathway which contributed to neuro-protective signaling cascades in stroke (Shi et al., 2011). (7) GP1BA was considered as one of the candidate ‘stroke risk’ genes affecting hemostasis (Stankovic and Majkic-Singh, 2010).

As shown in Table 2, the patients P1 and P3 had better outcomes compared with P2 and P4. We than compared the ranking of these genes based on their ECSs to find the correlations between them and the clinical outcome. As shown in Table 3, six genes were shown to be positively correlated with the recovery degree of IS patients: ELANE, FGF23, HBB, PIEZO1, RASA4, and PRTN3. Gene CPM was shown to be correlated with clinical outcome negatively. Of which, the relationships between stroke and PIEZO1, FGF23, RASA4 had been validated in former researches as mentioned above (Shi et al., 2011; Wright et al., 2014; Retailleau et al., 2015).

TABLE 3
www.frontiersin.org

Table 3. Expression changed score and correlations with IS outcome for seven hemoglobin related genes.

Biological Validation of Candidate Bio-Markers

Real-time PCR was performed for the following genes, which were identified as candidate bio-markers without supporting literature: ELANE, HBB, PRTN3, CPM, FN3KRP, POU5F1, SMG5, UGT1A1, CSNK2A1, TUBGCP3, and PKM. Of which, ELANE, HBB, PRTN3, POU5F1, and PKM were shown to be up-regulated in IS patients in both ISA and non-ISA groups compared with the control group. It was interesting to find out that the rank of the genes’ mean expression values were: ISA, non-ISA, TIA, control (see Figures 3A–E for details). Compared with this, CPM was shown to have trend to the contray with the mean expression rank as: control, TIA, non-ISA, ISA (see Figure 3F), which was in accordance with our results listed in Table 3. Taken together, we believe that these genes could be considered as new predictors for the recovery of IS patients.

FIGURE 3
www.frontiersin.org

Figure 3. RT-PCR results for six genes with significant differences between different sample groups; (A) ELANE; (B) HBB; (C) PRTN3; (D) POU5F1; (E) PKM; (F) CPM.

Ethics Statement

The study was approved by SIAT Institutional Review Committee with IRB number SIAT-IRB-16515-H0107 and all the procedures were in accordance with the SIAT-IRB guidelines and the Declaration of Helsinki.

Author Contributions

YW performed the bio-informatics analyses and wrote the manuscript. XH collected all the blood samples for sequencing. JL analyzed the clinical information. XZ performed the PCR experiment. HY directed the clinical analyses and revised the manuscript. YC directed the bio-informatics analyses and revised the manuscript. All authors agreed to be accountable for the content of the work.

Funding

This study was partly funded by National Natural Science Foundation of China (61702496 and 81601575), Shenzhen Technology R&D Program (JSGG20170413152936281), and Shenzhen Science and Technology Research Funding (20170502171625936).

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/fgene.2019.00456/full#supplementary-material

TABLE S1 | Clinical information of patients.

TABLE S2 | Common pathways of different groups.

References

Back, C., Thiesen, K. L., Skovgaard, K., Edvinsson, L., Jensen, L. T., Larsen, V. A., et al. (2015). RAAS and stress markers in acute ischemic stroke: preliminary findings. Acta Neurol. Scand. 131, 132–139. doi: 10.1111/ane.12298

PubMed Abstract | CrossRef Full Text | Google Scholar

Big Data Center Members. (2018). Database resources of the BIG data center in 2018. Nucleic Acids Res. 46, D14–D20. doi: 10.1093/nar/gkx897

PubMed Abstract | CrossRef Full Text | Google Scholar

Bokhari, F. A., Shakoori, T. A., Butt, A., and Ghafoor, F. (2014). TNF-alpha: a risk factor for ischemic stroke. J. Ayub Med. Coll. Abbottabad 26, 111–114.

Google Scholar

Bressler, J., Franceschini, N., Demerath, E. W., Mosley, T. H., Folsom, A. R., and Boerwinkle, E. (2015). Sequence variation in telomerase reverse transcriptase (TERT) as a determinant of risk of cardiovascular disease: the atherosclerosis risk in communities (ARIC) study. BMC Med. Genet. 16:52. doi: 10.1186/s12881-015-0194-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaturvedi, M., and Kaczmarek, L. (2014). Mmp-9 inhibition: a therapeutic strategy in ischemic stroke. Mol. Neurobiol. 49, 563–573. doi: 10.1007/s12035-013-8538-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Deddens, L. H., van Tilborg, G. A., van der Toorn, A., van der Marel, K., Paulis, L. E., van Bloois, L., et al. (2013). MRI of ICAM-1 upregulation after stroke: the importance of choosing the appropriate target-specific particulate contrast agent. Mol. Imaging Biol. 15, 411–422. doi: 10.1007/s11307-013-0617-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Dubyk, M. D., Card, R. T., Whiting, S. J., Boyle, C. A., Zlotkin, S. H., and Paterson, P. G. (2012). Iron deficiency anemia prevalence at first stroke or transient ischemic attack. Can. J. Neurol. Sci. 39, 189–195.

Google Scholar

Gao, T., Huang, L., Fu, Q., and Bai, Y. (2015). Association of polymorphisms in the AGT gene(M235T, T174M) with ischemic stroke in the Chinese population. JRAAS 16, 681–686. doi: 10.1177/1470320315583600

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, W. Y., Chen, I. C., Meng, L., Weng, W. C., and Peng, T. I. (2009). The influence of anemia on clinical presentation and outcome of patients with first-ever atherosclerosis-related ischemic stroke. J. Clin. Neurosci. 16, 645–649. doi: 10.1016/j.jocn.2008.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Ding, Y., Tang, J., and Guo, F. (2018a). MDA-SKF: similarity kernel fusion for accurately discovering mirna-disease association. Front. Genet. 9:618. doi: 10.3389/fgene.2018.00618

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Xiao, Y., Ding, Y., Tang, J., and Guo, F. (2018b). FKL-Spa-LapRLS: an accurate method for identifying human microRNA-disease association. BMC Genomics 19(Suppl. 10):911. doi: 10.1186/s12864-018-5273-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Xiao, Y., Ding, Y., Tang, J., and Guo, F. (2019). Discovering cancer subtypes via an accurate fusion strategy on multiple profile data. Front. Genet. 10:20. doi: 10.3389/fgene.2019.00020

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaiafa, G., Kanellos, I., Savopoulos, C., Kakaletsis, N., Giannakoulas, G., and Hatzitolios, A. I. (2015). Is anemia a new cardiovascular risk factor? Int. J. Cardiol. 186, 117–124. doi: 10.1016/j.ijcard.2015.03.159

PubMed Abstract | CrossRef Full Text | Google Scholar

Kellert, L., Martin, E., Sykora, M., Bauer, H., Gussmann, P., Diedler, J., et al. (2011). Cerebral oxygen transport failure?: decreasing hemoglobin and hematocrit levels after ischemic stroke predict poor outcome and mortality: stroke: relevant impact of hemoglobin, hematocrit and transfusion (STRAIGHT)–an observational study. Stroke 42, 2832–2837. doi: 10.1161/STROKEAHA.110.606665

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimberly, W. T., Lima, F. O., O’Connor, S., and Furie, K. L. (2013). Sex differences and hemoglobin levels in relation to stroke outcomes. Neurology 80, 719–724. doi: 10.1212/WNL.0b013e31828250ff

PubMed Abstract | CrossRef Full Text | Google Scholar

Mascia, L. (2009). Acute lung injury in patients with severe brain injury: a double hit model. Neurocrit. Care 11, 417–426. doi: 10.1007/s12028-009-9242-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Milionis, H., Papavasileiou, V., Eskandari, A., D’Ambrogio-Remillard, S., Ntaios, G., and Michel, P. (2015). Anemia on admission predicts short- and long-term outcomes in patients with acute ischemic stroke. Int. J. Stroke 10, 224–230. doi: 10.1111/ijs.12397

PubMed Abstract | CrossRef Full Text | Google Scholar

Njajou, O. T., Hollander, M., Koudstaal, P. J., Hofman, A., Witteman, J. C., Breteler, M. M., et al. (2002). Mutations in the hemochromatosis gene (HFE) and stroke. Stroke 33, 2363–2366.

Google Scholar

Park, Y. H., Kim, B. J., Kim, J. S., Yang, M. H., Jang, M. S., Kim, N., et al. (2013). Impact of both ends of the hemoglobin range on clinical outcomes in acute ischemic stroke. Stroke 44, 3220–3222. doi: 10.1161/STROKEAHA.113.002672

PubMed Abstract | CrossRef Full Text | Google Scholar

Retailleau, K., Duprat, F., Arhatte, M., Ranade, S. S., Peyronnet, R., Martins, J. R., et al. (2015). Piezo1 in smooth muscle cells is involved in hypertension-dependent arterial remodeling. Cell Rep. 13, 1161–1171. doi: 10.1016/j.celrep.2015.09.072

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahin, M. H., and Johnson, J. A. (2016). Mechanisms and pharmacogenetic signals underlying thiazide diuretics blood pressure response. Curr. Opin. Pharmacol. 27, 31–37. doi: 10.1016/j.coph.2016.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, G. X., Andres, D. A., and Cai, W. (2011). Ras family small GTPase-mediated neuroprotective signaling in stroke. Cent. Nerv. Syst. Agents Med. Chem. 11, 114–137.

PubMed Abstract | Google Scholar

Stankovic, S., and Majkic-Singh, N. (2010). Genetic aspects of ischemic stroke: coagulation, homocysteine, and lipoprotein metabolism as potential risk factors. Crit. Rev. Clin. Lab. Sci. 47, 72–123. doi: 10.3109/10408361003791520

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabib, A., Karbian, N., and Mevorach, D. (2017). Demyelination, strokes, and eculizumab: lessons from the congenital CD59 gene mutations. Mol. Immunol. 89, 69–72. doi: 10.1016/j.molimm.2017.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Liu, J., Yu, H., and Cai, Y. (2017a). “Analyses of the Dual Immune Roles Cytokines Play in Ischemic Stroke,” in Health Information Science. HIS 2017. Lecture Notes in Computer Science, Vol. 10594, eds S. Siuly, et al. (Cham: Springer).

Google Scholar

Wang, Y., Song, F., Zhu, J., Zhang, S., Yang, Y., Chen, T., et al. (2017b). GSA: genome sequence archivesup. Genomics Proteo. Bioinform. 15, 14–18. doi: 10.1016/j.gpb.2017.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, C. B., Dong, C., Stark, M., Silverberg, S., Rundek, T., Elkind, M. S., et al. (2014). Plasma FGF23 and the risk of stroke: the Northern Manhattan Study (NOMAS). Neurology 82, 1700–1706. doi: 10.1212/WNL.0000000000000410

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: hemoglobin, ischemic stroke, bioinformatics, anemia, predictive

Citation: Wang Y, Huang X, Liu J, Zhao X, Yu H and Cai Y (2019) A Systems Analysis of the Relationships Between Anemia and Ischemic Stroke Rehabilitation Based on RNA-Seq Data. Front. Genet. 10:456. doi: 10.3389/fgene.2019.00456

Received: 02 February 2019; Accepted: 30 April 2019;
Published: 24 May 2019.

Edited by:

Dariusz Mrozek, Silesian University of Technology, Poland

Reviewed by:

Fei Guo, Tianjin University, China
Haijing Wang, Microsoft (United States), United States

Copyright © 2019 Wang, Huang, Liu, Zhao, Yu and Cai. 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: Haibo Yu, MTM2MDMwNjYwOThAMTYzLmNvbQ==; Yunpeng Cai, eXAuY2FpQHNpYXQuYWMuY24=

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.