# iTRAQ-Based Proteomics Analysis Reveals the Effect of Rhubarb in Rats with Ischemic Stroke

Authors: Xiangping Lin; Tao Liu; Pengfei Li; Zehui He; Yuanyuan Zhong; Hanjin Cui; Jiekun Luo; Yang Wang; Tao Tang
Journal: BioMed Research International
 (2018)

---

## Abstract

Background.Rhubarb, a traditional Chinese medicine, promotes viscera and remove blood stasis. Rhubarb is skilled in smoothening meridians, improving blood circulation which exhibits better efficacy on cerebral ischemic stroke. In this study, we aimed to analyze the underlying mechanisms of rhubarb which treated rats of middle cerebral artery occlusion (MCAO) model according to an iTRAQ-based proteomics and bioinformatics analysis. 30 rats were randomly allocated into three groups including sham group (SG), model group (MG), and rhubarb group (RG). Rhubarb group was given a gavage of rhubarb decoction at dose of 3 g/kg and the remaining groups were prepared with normal saline by gavage. Rats from MG and RG were induced into MCAO model. The effects of rhubarb were estimated by Modified Neurological Severity Score (mNSS) and cerebral infarct volume. The brain tissues were measured via the quantitative proteomic approach of iTRAQ coupled to liquid chromatography-tandem mass spectrometry (LC-MS/MS). Furthermore, the bioinformatics analysis of overlapping differentially expression proteins (DEPs) was conducted by DAVID, KEGG, and Cytoscape. Specific selective DEPs were validated by Western blotting. Rats treated with rhubarb after MCAO showed a significant reduction on mNSS and cerebral infarct volume compared with MG. In MG versus SG and RG versus MG, we identified a total of 4578 proteins, of which 287 were DEPs. There were 76 overlapping DEPs between MG versus SG and RG versus MG. Through bioinformatics analysis, 14 associated pathways were searched including cGMP-PKG signaling pathway, tuberculosis, synaptic vesicle cycle, amyotrophic lateral sclerosis, long-term potentiation, and so on. 76 overlapping DEPs mainly involved synaptic vesicle cycling biological processes based on GO annotation. Further, the selective overlapping DEPs were verified at the protein level by using Western blotting. Our present study reveals that rhubarb highlights promising neuroprotective effect. Rhubarb exerts novel therapeutic action via modulating multiple proteins, targets, and pathways.


---

## Body

## 1. Introduction

Stroke is the second leading cause of death and has remained one of the biggest killers globally in the last 15 years [1, 2]. Approximately 11.6 million ischemic stroke events are in low- and middle-income countries [3]. In 2013, there was an estimated 10.3 million new strokes, 67% of which were ischemic stroke in the United States [4]. The prevalence rate of ischemic stroke nearly doubled from 1990 to 2013 globally [5]. In view of the above facts, the government from each country tends to develop the prevention of ischemic stroke as a national priority.Ischemic stroke generally leads to a cascade of pathological mechanisms, contributing to the irreversible tissue damage including mitochondrial dysfunction, cerebral edema, cerebral hemorrhage, and neuronal death [6, 7]. Pathogenically, ischemic stroke accounts for 85% of all conditions caused by vessel occlusions [8]. These patients were provided antiplatelet or anticoagulant drugs at acute ischemic stroke onset [9]. Unfortunately, these drugs could increase the risk of fatal or disabling intracranial hemorrhage, which offset any benefits in current ischemic stroke [10]. Thus, the novel therapeutic agents with safety and efficacy have attracted much attention of neuroscientists and doctors.The clinical applications and experimental studies have reported that traditional Chinese medicines including Radix/rhizome notoginseng extract, Buyang Huanwu decoction, and Danhong injection are ascendant against convalescent ischemic stroke [11–13]. Nevertheless, few studies and clinical trials focus on the earlier phase of acute therapy for ischemic stroke. According to the traditional Chinese medicine theory, the core of ischemic stroke refers to binding phlegm-heat and blood stasis syndrome. In the light of this, rhubarb (Rheum palmatumL. orRheum tanguticumMaxim, dahuang in China) is skilled in smoothening meridians, improving blood circulation which exhibits better efficacy on cerebral ischemia [14]. Modern studies show that rhubarb protects the blood brain barrier and provides neuroprotective effect via an antioxidative molecular mechanism in brain injury [15, 16]. However, there is no explicit mechanism responsible by rhubarb treatment for ischemic stroke.Due to the effects of herbal medicine via multiple targets and multidirectional treatments, single pharmacological action cannot unfold the effects of rhubarb to treat ischemic stroke. Within recent years, omics-based methods are applied to uncover the multipathways of herbal medicine for disease treatment [17]. The proteomics is rapidly emerging as an essential tool to gain insights into complex systems at the level of protein expression. It has been widely used to detect novel targets, study interaction of proteins, and discover the associated mechanisms of action in the life sciences [18, 19]. The main method for proteomics study includes iTRAQ which combines with multidimensional liquid chromatography (LC) and mass spectroscopic (MS) analysis [20]. It has become a crucial tool in quantitative proteomics due to its high throughput, sensitivity, and accuracy [21, 22]. This approach enables comparison between normal and disease samples to uncover the dysregulation of proteins [23]. The iTRAQ has been frequently utilized in the studies of traditional Chinese herbs such as Gastrodia elata blume and Pseudostellaria heterophylla [24, 25]. Hence, it is feasible to reveal the multiple targeting effects of rhubarb in treating acute ischemic stroke-based on quantitative proteomics.In this study, the iTRAQ-based quantitative proteomic approach was applied to identify disease relevant proteins in rat model of ischemic stroke. We next used the bioinformatics analyses including KEGG, DAVID, and Cytoscape to exhibit the potential targets and molecular signal pathways of rhubarb to treat acute ischemic stroke. The further verifications were conducted by Western blotting. The present study could help to facilitate the rhubarb as a potential therapeutic agent for acute ischemic stroke.


## 2. Materials and Methods

### 2.1. Experimental Animals

Male Sprague-Dawley (SD) rats weighing 220–240g were purchased from Laboratory Animal Centre of Central South University (CSU) and allowed to acclimatize for at least 5 days. The environmental conditions were maintained in specific pathogen free grade. 30 rats were randomly allocated into three groups including sham group (SG), model group (MG), and rhubarb group (RG). RG was given a gavage of rhubarb decoction at dose of 3 g/kg and the rest groups were prepared with normal saline by gavage. The animal procedures were approved and verified by the Animal Ethics Committee of CSU.


### 2.2. Chemicals and Drugs

The antibodies of Synapsin-1 andβ-tubulin were from Proteintech Group (Chicago, USA). ERK1/2 antibody was from the American Abcam Company (Cambridge, UK). All the other reagents were of analytical grade or better. Rhubarbs (voucher specimen no. 2015101605, Gansu, China) were purchased from Pharmacy of Traditional Medicine of Xiangya Hospital Central South University (Changsha, Hunan). Dried rhubarb was crushed into powder and soaked in distilled water for 30 min. They were boiled for 10 min at 100°C. The extraction yield was around 1:12 w/v. The decoction was stored under 4°C. The method of identifying and quantifying the principal components on rhubarb decoction had been reported in our previous study [26].


### 2.3. Establishment of the Middle Cerebral Artery Occlusion (MCAO) Model

Rat was produced by transient intracranial MCAO, as formerly described [27]. Rats were placed in a supine position. The right common carotid artery (CCA), the right external carotid artery (ECA), and the right internal carotid artery (ICA) were exposed and isolated. A 4-0 nylon monofilament suture (Beijing Cinontech Co., Ltd., Beijing, China) coated with 1% poly-L-lysine was inserted from the CCA to the ICA and gently advanced to occlude the origin of the middle cerebral artery (MCA). The suture had reached the circle of Willis. After 2h of occlusion, the MCA-suture was carefully removed for reperfusion.The suture was withdrawn to allow reperfusion 24h after MCAO. Rats in SG received all the surgical procedures except that the arteries were not occluded. SG and MG were given physiological saline under the same conditions. After assessment of the neurological deficits, rats were sacrificed. All surgical operations were finished in a sterile environment with body temperatures maintained at 37± 0.5°C.


### 2.4. Neurological Evaluation and TTC Staining

Neurological functional assessments were performed at 24 h after MCAO using Modified Neurological Severity Score (mNSS) [28, 29]. Neurological function was graded on a scale of 0 to 18 (normal score, 0; maximal deficit score, 18). The more severe the injury, the higher the scores of the rats. The dead rats would be excluded. The brain tissues were rapidly removed and frozen for 20 min at -20°C. Frozen cerebra were sliced into 2 mm thick coronal section and incubated in 2,3,5-triphenyltetrazolium chloride (TTC) at 37°C for 30 min. Then the infarct volume of each section was measured with an image analysis system (Image-ProPlus 6.0, Media Cybernetics, MD, USA). % Infarct volume = (white infarct area ×thickness) / (whole slice area × thickness) ×100% [30].


### 2.5. Proteomics

#### 2.5.1. Protein Extraction

Briefly, the total protein from each cerebral tissue sample was ground to powder with liquid nitrogen and dissolved in lysis solution. The cell lysate was centrifuged at 12,000 rpm for 20 min and the supernatant was collected. All extraction steps were carried out at 4°C.


#### 2.5.2. iTRAQ Labeling

From each sample, 200 mg of protein was reduced, alkylated, and trypsin-digested overnight at 37°C. The labeling was performed according to the manufacturer's instructions for the iTRAQ Reagents 8-plex kit (AB Sciex, Foster City, CA, USA). Samples from MG were labeled with tag 113. The SG was labeled with 114. RG was labeled with tag 116.


#### 2.5.3. Peptide Fractionation with Strong Cation Exchange (SCX)

For SCX chromatography using the Shimadzu LC-20AD HPLC Pump system (Shimadzu Corp, Kyoto, Japan), the peptides from digestion were reconstituted with buffer A (20mM HCOONH4, PH 10) and loaded onto a 2 × 150 mm Gemini-NX C18 SCX column containing 3-μm particles (Phenomenex) with buffer B (20 mM HCOONH4 80% acetone PH 10). Elution was monitored by measuring absorbance at 214nm/280 nm. The flow rate was 200μl/min. Eluted peptides were collected as 24 fractions and vacuum-dried. Each final fraction was resuspended in 50% trifluoroacetic acid.


#### 2.5.4. LC-MS/MS Analysis

The peptides were analyzed on a Q Exzctive mass spectrometer (Thermo Scientific) and identified using the Thermo Dionex Ultimate 3000 RSLC nanosystem. The binary buffer system was consisting of 0.1% formic acid (buffer A) and 84% CAN in 0.1% formic acid (buffer B). Each fraction was dissolved in 0.1% formic acid and 2% ACN. The supernatant was loaded onto a PepMap C18 RP column (2μm, 75 μm × 150 mm, 100 A). The peptides separated under a 65min linear gradient from 4% buffer B to 90% buffer B at a constant flow rate of 300 nl/min.The eluted peptides were directly detected by Q Exactive online. MS data were acquired in a data-dependent mode using a top 20 high-energy collisional dissociation (HCD) method, which dynamically chose the most abundant precursor ions from the scan range (350–1800 m/z). The dynamic exclusion was set to 40 s, and the automatic obtaining control (AGC) target value was 3E6. The underfill ratio was defined as 0.1%.


#### 2.5.5. Proteomics Data Analysis

The obtained raw data were searched with Protein Pilot Software (AB Sciex, Foster, CA, USA; version 5.0.1) against the UniProt rat database. The search parameters were the following: iTRAQ 8 plex peptide labeled was selected for sample type and trypsin was selected as the digestion enzyme. Biological modifications were selected as ID focus. Only proteins with unused value more than 1.3 were considered for further analysis (equal to protein confidence was set at 95%). The data were based on a false discovery rate (FDR) ≤1% confidence for protein identification.


### 2.6. Bioinformatics Analysis

The networks were generated and analyzed by using Cytoscape 3.5.1. We used ClueGO, a plugin of Cytoscape, to constitute the network for Gene Ontology (GO) mapping and annotation. Kyoto encyclopedia of genes and genomes (KEGG) pathways analysis was performed using DAVID Bioinformatics Resources v6.8 (https://david.ncifcrf.gov/).Pvalue ≤0.05 was considered to indicate significant pathways.


### 2.7. Western Blotting

The frozen cerebral tissues (n=4 per group) were collected in ice-cold lysis buffer. The lysates were centrifuged at 12000 g for 10 min at 4°C. The protein concentration was measured by the BCA protein assay kit (Well Biotechnology Company, Changsha, China). The proteins were loaded and separated by 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to polyvinylidene fluoride (PVDF) membranes. The membrane was blocked at room temperature with a 5% milk-Tris-buffered saline and Tween 20 (TBST) solution for 1 h and then incubated with the anti-Syn1 (1:2000), anti-ERK1/2 (1:10000), and anti-β-tubulin (1:1000) overnight at 4°C. After being washed thrice with TBST (15 min each), the membranes were incubated with horseradish peroxidase-labeled secondary antibody (1:6000, Proteintech Group, Chicago, USA). After incubation finished, the membranes were washed three times (15 min each) and visualized using enhanced chemiluminescence detection system (ECL; Thermo Fisher Scientific, Shanghai, China) for 5 min before exposure to X-ray film. Western blotting bands were scanned and analyzed with UN-Scan-It 6.1 software (Silk Scientific, Inc., Orem, UT) [31].


### 2.8. Statistics Analysis

All data are expressed as mean ± standard deviation (SD) and evaluated by one-way ANOVA with SPSS version 18.0. The results in each group were compared by t-test. Ap value<0.05 was considered to be statistically significant.


## 2.1. Experimental Animals

Male Sprague-Dawley (SD) rats weighing 220–240g were purchased from Laboratory Animal Centre of Central South University (CSU) and allowed to acclimatize for at least 5 days. The environmental conditions were maintained in specific pathogen free grade. 30 rats were randomly allocated into three groups including sham group (SG), model group (MG), and rhubarb group (RG). RG was given a gavage of rhubarb decoction at dose of 3 g/kg and the rest groups were prepared with normal saline by gavage. The animal procedures were approved and verified by the Animal Ethics Committee of CSU.


## 2.2. Chemicals and Drugs

The antibodies of Synapsin-1 andβ-tubulin were from Proteintech Group (Chicago, USA). ERK1/2 antibody was from the American Abcam Company (Cambridge, UK). All the other reagents were of analytical grade or better. Rhubarbs (voucher specimen no. 2015101605, Gansu, China) were purchased from Pharmacy of Traditional Medicine of Xiangya Hospital Central South University (Changsha, Hunan). Dried rhubarb was crushed into powder and soaked in distilled water for 30 min. They were boiled for 10 min at 100°C. The extraction yield was around 1:12 w/v. The decoction was stored under 4°C. The method of identifying and quantifying the principal components on rhubarb decoction had been reported in our previous study [26].


## 2.3. Establishment of the Middle Cerebral Artery Occlusion (MCAO) Model

Rat was produced by transient intracranial MCAO, as formerly described [27]. Rats were placed in a supine position. The right common carotid artery (CCA), the right external carotid artery (ECA), and the right internal carotid artery (ICA) were exposed and isolated. A 4-0 nylon monofilament suture (Beijing Cinontech Co., Ltd., Beijing, China) coated with 1% poly-L-lysine was inserted from the CCA to the ICA and gently advanced to occlude the origin of the middle cerebral artery (MCA). The suture had reached the circle of Willis. After 2h of occlusion, the MCA-suture was carefully removed for reperfusion.The suture was withdrawn to allow reperfusion 24h after MCAO. Rats in SG received all the surgical procedures except that the arteries were not occluded. SG and MG were given physiological saline under the same conditions. After assessment of the neurological deficits, rats were sacrificed. All surgical operations were finished in a sterile environment with body temperatures maintained at 37± 0.5°C.


## 2.4. Neurological Evaluation and TTC Staining

Neurological functional assessments were performed at 24 h after MCAO using Modified Neurological Severity Score (mNSS) [28, 29]. Neurological function was graded on a scale of 0 to 18 (normal score, 0; maximal deficit score, 18). The more severe the injury, the higher the scores of the rats. The dead rats would be excluded. The brain tissues were rapidly removed and frozen for 20 min at -20°C. Frozen cerebra were sliced into 2 mm thick coronal section and incubated in 2,3,5-triphenyltetrazolium chloride (TTC) at 37°C for 30 min. Then the infarct volume of each section was measured with an image analysis system (Image-ProPlus 6.0, Media Cybernetics, MD, USA). % Infarct volume = (white infarct area ×thickness) / (whole slice area × thickness) ×100% [30].


## 2.5. Proteomics

### 2.5.1. Protein Extraction

Briefly, the total protein from each cerebral tissue sample was ground to powder with liquid nitrogen and dissolved in lysis solution. The cell lysate was centrifuged at 12,000 rpm for 20 min and the supernatant was collected. All extraction steps were carried out at 4°C.


### 2.5.2. iTRAQ Labeling

From each sample, 200 mg of protein was reduced, alkylated, and trypsin-digested overnight at 37°C. The labeling was performed according to the manufacturer's instructions for the iTRAQ Reagents 8-plex kit (AB Sciex, Foster City, CA, USA). Samples from MG were labeled with tag 113. The SG was labeled with 114. RG was labeled with tag 116.


### 2.5.3. Peptide Fractionation with Strong Cation Exchange (SCX)

For SCX chromatography using the Shimadzu LC-20AD HPLC Pump system (Shimadzu Corp, Kyoto, Japan), the peptides from digestion were reconstituted with buffer A (20mM HCOONH4, PH 10) and loaded onto a 2 × 150 mm Gemini-NX C18 SCX column containing 3-μm particles (Phenomenex) with buffer B (20 mM HCOONH4 80% acetone PH 10). Elution was monitored by measuring absorbance at 214nm/280 nm. The flow rate was 200μl/min. Eluted peptides were collected as 24 fractions and vacuum-dried. Each final fraction was resuspended in 50% trifluoroacetic acid.


### 2.5.4. LC-MS/MS Analysis

The peptides were analyzed on a Q Exzctive mass spectrometer (Thermo Scientific) and identified using the Thermo Dionex Ultimate 3000 RSLC nanosystem. The binary buffer system was consisting of 0.1% formic acid (buffer A) and 84% CAN in 0.1% formic acid (buffer B). Each fraction was dissolved in 0.1% formic acid and 2% ACN. The supernatant was loaded onto a PepMap C18 RP column (2μm, 75 μm × 150 mm, 100 A). The peptides separated under a 65min linear gradient from 4% buffer B to 90% buffer B at a constant flow rate of 300 nl/min.The eluted peptides were directly detected by Q Exactive online. MS data were acquired in a data-dependent mode using a top 20 high-energy collisional dissociation (HCD) method, which dynamically chose the most abundant precursor ions from the scan range (350–1800 m/z). The dynamic exclusion was set to 40 s, and the automatic obtaining control (AGC) target value was 3E6. The underfill ratio was defined as 0.1%.


### 2.5.5. Proteomics Data Analysis

The obtained raw data were searched with Protein Pilot Software (AB Sciex, Foster, CA, USA; version 5.0.1) against the UniProt rat database. The search parameters were the following: iTRAQ 8 plex peptide labeled was selected for sample type and trypsin was selected as the digestion enzyme. Biological modifications were selected as ID focus. Only proteins with unused value more than 1.3 were considered for further analysis (equal to protein confidence was set at 95%). The data were based on a false discovery rate (FDR) ≤1% confidence for protein identification.


## 2.5.1. Protein Extraction

Briefly, the total protein from each cerebral tissue sample was ground to powder with liquid nitrogen and dissolved in lysis solution. The cell lysate was centrifuged at 12,000 rpm for 20 min and the supernatant was collected. All extraction steps were carried out at 4°C.


## 2.5.2. iTRAQ Labeling

From each sample, 200 mg of protein was reduced, alkylated, and trypsin-digested overnight at 37°C. The labeling was performed according to the manufacturer's instructions for the iTRAQ Reagents 8-plex kit (AB Sciex, Foster City, CA, USA). Samples from MG were labeled with tag 113. The SG was labeled with 114. RG was labeled with tag 116.


## 2.5.3. Peptide Fractionation with Strong Cation Exchange (SCX)

For SCX chromatography using the Shimadzu LC-20AD HPLC Pump system (Shimadzu Corp, Kyoto, Japan), the peptides from digestion were reconstituted with buffer A (20mM HCOONH4, PH 10) and loaded onto a 2 × 150 mm Gemini-NX C18 SCX column containing 3-μm particles (Phenomenex) with buffer B (20 mM HCOONH4 80% acetone PH 10). Elution was monitored by measuring absorbance at 214nm/280 nm. The flow rate was 200μl/min. Eluted peptides were collected as 24 fractions and vacuum-dried. Each final fraction was resuspended in 50% trifluoroacetic acid.


## 2.5.4. LC-MS/MS Analysis

The peptides were analyzed on a Q Exzctive mass spectrometer (Thermo Scientific) and identified using the Thermo Dionex Ultimate 3000 RSLC nanosystem. The binary buffer system was consisting of 0.1% formic acid (buffer A) and 84% CAN in 0.1% formic acid (buffer B). Each fraction was dissolved in 0.1% formic acid and 2% ACN. The supernatant was loaded onto a PepMap C18 RP column (2μm, 75 μm × 150 mm, 100 A). The peptides separated under a 65min linear gradient from 4% buffer B to 90% buffer B at a constant flow rate of 300 nl/min.The eluted peptides were directly detected by Q Exactive online. MS data were acquired in a data-dependent mode using a top 20 high-energy collisional dissociation (HCD) method, which dynamically chose the most abundant precursor ions from the scan range (350–1800 m/z). The dynamic exclusion was set to 40 s, and the automatic obtaining control (AGC) target value was 3E6. The underfill ratio was defined as 0.1%.


## 2.5.5. Proteomics Data Analysis

The obtained raw data were searched with Protein Pilot Software (AB Sciex, Foster, CA, USA; version 5.0.1) against the UniProt rat database. The search parameters were the following: iTRAQ 8 plex peptide labeled was selected for sample type and trypsin was selected as the digestion enzyme. Biological modifications were selected as ID focus. Only proteins with unused value more than 1.3 were considered for further analysis (equal to protein confidence was set at 95%). The data were based on a false discovery rate (FDR) ≤1% confidence for protein identification.


## 2.6. Bioinformatics Analysis

The networks were generated and analyzed by using Cytoscape 3.5.1. We used ClueGO, a plugin of Cytoscape, to constitute the network for Gene Ontology (GO) mapping and annotation. Kyoto encyclopedia of genes and genomes (KEGG) pathways analysis was performed using DAVID Bioinformatics Resources v6.8 (https://david.ncifcrf.gov/).Pvalue ≤0.05 was considered to indicate significant pathways.


## 2.7. Western Blotting

The frozen cerebral tissues (n=4 per group) were collected in ice-cold lysis buffer. The lysates were centrifuged at 12000 g for 10 min at 4°C. The protein concentration was measured by the BCA protein assay kit (Well Biotechnology Company, Changsha, China). The proteins were loaded and separated by 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to polyvinylidene fluoride (PVDF) membranes. The membrane was blocked at room temperature with a 5% milk-Tris-buffered saline and Tween 20 (TBST) solution for 1 h and then incubated with the anti-Syn1 (1:2000), anti-ERK1/2 (1:10000), and anti-β-tubulin (1:1000) overnight at 4°C. After being washed thrice with TBST (15 min each), the membranes were incubated with horseradish peroxidase-labeled secondary antibody (1:6000, Proteintech Group, Chicago, USA). After incubation finished, the membranes were washed three times (15 min each) and visualized using enhanced chemiluminescence detection system (ECL; Thermo Fisher Scientific, Shanghai, China) for 5 min before exposure to X-ray film. Western blotting bands were scanned and analyzed with UN-Scan-It 6.1 software (Silk Scientific, Inc., Orem, UT) [31].


## 2.8. Statistics Analysis

All data are expressed as mean ± standard deviation (SD) and evaluated by one-way ANOVA with SPSS version 18.0. The results in each group were compared by t-test. Ap value<0.05 was considered to be statistically significant.


## 3. Results

### 3.1. Rats Treated with Rhubarb Show Attenuated Infarction Volume Rate and Improved Neurological Deficits

After 24 h reperfusion, 3 rats were dead and neurological assessments were conducted for the remaining rats. TTC staining visually showed the histological changes of cerebral tissues (Figure1(a)). SG appeared uniform red in color, indicating no infarction. Infraction volume of percentage showed a dramatic increase in MG (18.82% ± 5.23%) compared with SG (Figure 1(b)). Rhubarb significantly reduced acute infarction volume ratio (11.38% ± 4.04%) after 24 h reperfusion. As can be seen in Figure 1(c), obvious differences in regard to MCAO abnormality were exhibited between SG and MG (p<0.01). Rats treated with rhubarb showed a significant reduction of neurological deficits on mNSS compared with MG (p<0.05).Figure 1
(a) Representative TTC stained brain sections of sham group (n=3), model group (n=5), and rhubarb group (n=5) treated ischemic animals. (b) Infarct volume was determined 24 h after MCAO. (c) There were significant differences in neurological deficits at 24 h after MCAO between sham group (n=9), model group (n=10), and rhubarb group (n=8). Data are presented as mean ± SD. *p<0.05, **p<0.01, and ***p<0.001.
(a) (b) (c)


### 3.2. Quantitative Proteomic Analysis of Brains with MCAO Rats

In order to comprehensively understand the mechanism of rhubarb in ischemic stroke, iTRAQ analysis was performed on the protein extracts from rat brain with or without rhubarb treatment. After merging data from the respective replicates, 4578 proteins were identified and showed different expression (fold change ≥ 1.2, FDR ≤ 1%). Only the upregulated proteins with 1.5 ≤ an average ratio change ≤ 5 or downregulated proteins with an average ratio change ≤ 0.667. Peptides (95%) ≥ 2 and unused value ≥ 2 in both replicates were considered to be significantly expressed and used for the following analysis.In MG versus SG and RG versus MG, 287 proteins were differentially expression proteins (DEPs). There were 147 DEPs in RG versus MG and 140 DEPs in MG versus SG, respectively. Figure2(a) illustrated the proteins between the RG versus MG and MG versus SG by a Venn diagram. Only 64 DEPs were changed in MG versus SG uniquely, whereas 71 DEPs were regulated by rhubarb treatment alone. 76 overlapping DEPs showed altered protein expression level in both experimental sets and rhubarb group (Figure 2(b)). Table 1 summarized detailed information of 76 overlapping DEPs, including the UniProt accession number, gene name, peptides, fold change, and unused value. 53 overlapping DEPs were specifically upregulated, while 23 overlapping DEPs were downregulated after rhubarb treatment. Further, overlapping DEPs were subjected to functional analysis with reference to KEGG database.Table 1
Quantitative information of the overlapping DEPs in RG versus MG.
No. Uniprot Accession Gene name Unused % Cov Protein identity Peptides AVG. (95) (95%) 1 G3V7C6 Tubb4b 2 76.64 Tubulin beta-4B chain 274 2.43 2 Q9WTY2 Rem2 2.04 4.98 GTP-binding protein REM 2 2 1.88 3 Q04940 Nrgn 4.01 46.15 Neurogranin 7 3.81 4 P10818 Cox6a1 4.08 27.93 Cytochrome c oxidase subunit 6A1, mitochondrial 3 2.21 5 P04094 Penk 4.1 11.15 Proenkephalin-A 3 1.70 6 P22057 Ptgds 4.23 29.63 Prostaglandin-H2 D-isomerase 6 1.76 7 P35332 Hpcal4 6.4 49.21 Hippocalcin-like protein 4 12 2.15 8 P60841 Ensa 7.12 47.11 Alpha-endosulfine 7 2.99 9 Q3ZAV2 Ybx1 7.86 27.30 Nuclease-sensitive element-binding protein 1 4 1.60 10 G3V997 Dcx 8.55 20.82 Neuronal migration protein doublecortin 9 1.80 11 Q6AXS4 Atp6ap2 9.02 16.00 Renin receptor 8 1.73 12 P37377 Snca 9.85 80.71 Alpha-synuclein 22 2.37 13 O35796 C1qbp 10.07 28.32 Complement component 1 Q subcomponent-binding protein, mitochondrial 8 1.63 14 B7X6I3 Cend1 10.49 56.38 C38 protein 19 2.00 15 F8WFH6 Fam131b 11.76 20.28 Protein FAM131B 7 1.89 16 Q5U2U8 Bag3 12 17.42 Bcl2-associated athanogene 3 9 1.60 17 P35213 Ywhab 12.04 78.46 14-3-3 protein beta/alpha 65 1.92 18 O70257 Stx7 12.21 39.85 Syntaxin-7 13 3.47 19 A0A0G2K0M8 Ncam1 12.94 50.30 Neural cell adhesion molecule 1 77 2.44 20 P07825 Syp 13.95 40.07 Synaptophysin 19 1.55 21 P07171 Calb1 16.2 41.38 Calbindin 16 1.51 22 D3ZUX5 Chchd3 16.88 35.68 MICOS complex subunit 13 3.34 23 P62982 Rps27a 17.28 54.49 Ubiquitin-40S ribosomal protein S27a 23 1.62 24 D4A0T0 Ndufb10 18.05 61.93 Protein Ndufb10 15 2.43 25 Q09073 Slc25a5 18.08 61.07 ADP/ATP translocase 2 46 2.03 26 A0A0H2UHV6 Ppp3r1 18.44 61.25 Calcineurin subunit B type 1 17 1.55 27 D4A5L9 LOC679794 18.53 71.43 Protein LOC679794 24 3.31 28 P61589 Rhoa 19.6 64.77 Transforming protein RhoA 24 2.52 29 F1M378 Unc13a 20.94 6.01 Protein unc-13 homolog A 11 4.26 30 O88767 Park7 21.87 91.53 Protein deglycase DJ-1 22 1.76 31 D4AA63 Ubqln2 24.33 27.12 Protein Ubqln2 22 1.53 32 P14604 Echs1 24.5 53.45 Enoyl-CoA hydratase, mitochondrial 20 2.09 33 P08461 Dlat 27.5 35.60 Dihydrolipoyllysine-residue acetyltransferase component of pyruvate dehydrogenase complex, mitochondrial 25 3.38 34 P63086 Mapk1 31.61 49.16 Mitogen-activated protein kinase 1 26 1.63 35 Q6P6R2 Dld 33.08 54.03 Dihydrolipoyl dehydrogenase, mitochondrial 31 1.98 36 Q05175 Basp1 34.26 86.82 Brain acid soluble protein 1 82 4.18 37 P62161 Calm1 34.4 89.93 Calmodulin 80 2.64 38 P11275 Camk2a 35.04 64.44 Calcium/calmodulin-dependent protein kinase type II subunit alpha 76 1.73 39 P35565 Canx 35.34 37.39 Calnexin 31 1.50 40 P07936 Gap43 37.7 84.51 Neuromodulin 48 3.07 41 Q9Z0W5 Pacsin1 38.67 56.46 Protein kinase C and casein kinase substrate in neurons protein 1 50 1.61 42 P21707 Syt1 40.02 45.61 Synaptotagmin-1 37 1.79 43 Q5RKI0 Wdr1 42.43 46.37 WD repeat-containing protein 1 33 1.53 44 P63102 Ywhaz 43.94 85.31 14-3-3 protein zeta/delta 77 1.90 45 P06761 Hspa5 57.41 55.20 78 kDa glucose-regulated protein 57 1.55 46 P05708 Hk1 76.09 47.17 Hexokinase-1 80 1.50 47 Q63198 Cntn1 77.28 48.29 Contactin-1 77 1.74 48 D3ZMI4 Epb41l1 82.1 34.24 Band 4.1-like protein 1 52 1.77 49 P09951 Syn1 91.64 72.30 Synapsin-1 152 2.34 50 F1LQ81 Nsf 104.81 69.22 N-ethylmaleimide sensitive fusion protein, isoform CRA_b 86 1.84 51 D3ZQQ5 Dnm1 112.88 68.06 Dynamin-1 107 1.52 52 G3V984 Bsn 155.27 32.42 Protein bassoon 111 1.53 53 F1M779 Cltc 171.28 58.69 Clathrin heavy chain 213 1.76 54 A0A0G2K839 Clcn5 2 4.41 Chloride channel protein 3 0.02 55 G3V8L6 Kcna6 2.13 6.04 Potassium voltage-gated channel subfamily A member 6 2 0.65 56 D3ZV91 L3hypdh 4 7.34 Protein L3hypdh 2 0.48 57 F1LSV4 Sptlc2 4 5.18 Protein Sptlc2 3 0.52 58 A0A0G2K1A2 Mpo 4.01 3.76 Protein Mpo 3 0.65 59 P28572 Slc6a9 4.94 8.62 Sodium- and chloride-dependent glycine transporter 1 5 0.60 60 P55009 Aif1 5.28 23.13 Allograft inflammatory factor 1 5 0.60 61 Q6AYR8 Scrn2 5.8 13.00 Secernin-2 6 0.43 62 A0A0G2JWX4 Krt2 6.06 9.36 Keratin, type II cytoskeletal 2 epidermal 4 0.45 63 G3V6M3 Syt2 6.27 29.38 Synaptotagmin II 17 0.48 64 Q9R0I8 Pip4k2a 7.18 22.66 Phosphatidylinositol 5-phosphate 4-kinase type-2 alpha 14 0.46 65 Q7TP91 Surf4 10 9.55 Ab1-205 9 0.55 66 P0DMW1 Hspa1b 16.76 41.34 Heat shock 70 kDa protein 1B 73 0.51 67 P60203 Plp1 17.43 31.41 Myelin proteolipid protein 50 0.20 68 Q63345 Mog 18.16 36.73 Myelin-oligodendrocyte glycoprotein 14 0.39 69 A0A0G2JWM2 Sirt2 31.35 52.32 NAD-dependent protein deacetylase sirtuin-2 34 0.61 70 Q63357 Myo1d 31.36 21.27 Unconventional myosin-Id 27 0.48 71 E9PTX9 Slc12a2 34.28 21.28 Protein Slc12a2 26 0.44 72 P02688 Mbp 35.6 66.67 Myelin basic protein 88 0.23 73 F1LRZ7 Nefh 54.4 41.26 Neurofilament heavy polypeptide 55 0.47 74 P19527 Nefl 55.56 58.49 Neurofilament light polypeptide 87 0.32 75 P13233 Cnp 62.44 74.52 2′,3′-cyclic-nucleotide 3′-phosphodiesterase 83 0.47 76 P12839 Nefm 79.88 49.88 Neurofilament medium polypeptide 86 0.49Figure 2
The number of differentially expressed proteins (DEPs) and their overlap was exhibited in Venn diagram (a). 76 proteins with upregulated and downregulated expressions between MG versus SG and RG versus MG was showed in our study (b). Significant changes in the protein ratios were changed at least 1.5-fold.
(a) (b)


### 3.3. KEGG Database Analysis

The analysis was performed by matching the overlapping DEPs to the proteins annotated with KEGG pathway database via DAVID Bioinformatics Resources v6.8. The identified pathways were potentially affected by the modification of the abundance of proteins in the rat brain. Then 75 overlapping DEPs were identified and 1 overlapping DEP was not by DAVID. Among these overlapping DEPs, only 46 (56%) overlapping DEPs had a KEGG Orthology ID and were involved in 26 pathways. 14 pathways were statistically significant withp value < 0.05. The KEGG analysis results were shown in Figure 3 and Table 2. There are almost all closely related to Alzheimer's disease, oocyte meiosis, Parkinson's disease, oxytocin signaling pathway, cGMP-PKG signaling pathway, tuberculosis, synaptic vesicle cycle, amyotrophic lateral sclerosis (ALS), long-term potentiation (LTP), carbon metabolism, neurotrophin signaling pathway, prion diseases, amphetamine addiction, and glioma.Table 2
KEGG pathways associated with the 76 overlapping DEPs by DAVID.
Term Count PValue Genes rno04114: Oocyte meiosis 6 2.47E-04 P63102, P62161, P35213, A0A0H2UHV6, P63086, P11275 rno05010: Alzheimer’s disease 6 0.00251 P62161, P10818, D4A0T0, A0A0H2UHV6, P63086, P37377 rno05012: Parkinson’s disease 5 0.008091 O88767, Q09073, P10818, D4A0T0, P37377 rno04921: Oxytocin signaling pathway 5 0.009657 P61589, P62161, A0A0H2UHV6, P63086, P11275 rno04022: cGMP-PKG signaling pathway 5 0.011878 P61589, P62161, Q09073, A0A0H2UHV6, P63086 rno05152: Tuberculosis 5 0.014949 P61589, P62161, A0A0H2UHV6, P63086, P11275 rno05014: Amyotrophic lateral sclerosis (ALS) 4 0.002947 P19527, P12839, A0A0H2UHV6, F1LRZ7 rno04721: Synaptic vesicle cycle 4 0.004142 F1LQ81, P21707, F1M779, F1M378 rno04720: Long-term potentiation 4 0.004732 P62161, A0A0H2UHV6, P63086, P11275 rno01200: Carbon metabolism 4 0.025472 P05708, Q6P6R2, P14604, P08461 rno04722: Neurotrophin signaling pathway 4 0.028277 P61589, P62161, P63086, P11275 rno05020: Prion diseases 3 0.012162 P06761, A0A0G2K0M8, P63086 rno05031: Amphetamine addiction 3 0.044482 P62161, A0A0H2UHV6, P11275 rno05214: Glioma 3 0.045744 P62161, P63086, P11275Figure 3
KEGG pathway analysis of overlapping DEPs identified by iTRAQ analysis.


### 3.4. Protein-Protein Interaction (PPI) Network Analysis

Further PPI network analysis was conduct by Cytoscape. PPI network (Figure4) showed that several significant proteins might function in response to ischemic stroke. Of these 75 identified overlapping DEPs, a total of 52 were linked into a protein interaction network and 23 did not show any link. The network analysis revealed potential relationship between Mapk1, Syn1, Calm1, Syp, Gap43, Hspa5, Mbp, and so on. The network view summarized the network of predicted associations for particular proteins. It provided us a crucial modeling scaffolds and data reduction way which made us got insight into the potential mechanisms of rhubarb in brain.Figure 4
Protein network connectivity of identified proteins in MCAO rat with or without rhubarb.


### 3.5. GO Analysis

To detect the potential role of rhubarb in MCAO rats, we allocated 76 overlapping DEPs to biological process, cell component, and molecular function based on GO annotation. Biological processes that overlapping DEPs participated were representatively involved in synaptic vesicle cycle (28%), neurotransmitter transport (22%), neurofilament bundle assembly (16%), presynaptic process involved in chemical synaptic transmission (11%), oxidative phosphorylation (10%), mitochondrial membrane organization (5%), substantia nigra development (3%), negative regulation of transcription from RNA polymerase II promoter in response to stress (3%), and response to monoamine (2%) (Figures5(a) and 5(b)).Figure 5
GO annotation analysis of the overlapping DEPs. The color had a one-to-one correspondence between (a), (c), (e) and (b), (d), (f), respectively. ((a) and (b)) GO annotation enriched by ClueGO according to biological process. ((c) and (d)) GO annotation enriched by ClueGO according to cellular components. ((e) and (f)) GO annotation enriched by ClueGO molecular function. Functionally grouped network with terms linked by nodes, and functionally related groups were partially overlapped; the node size represented the term enrichment significance.
(a) (b) (c) (d) (e) (f)Cell component suggested that proteins were located-representatively in the axon (48%), growth cone (8%) neurofilament (8%), melanosome (8%), cell projection cytoplasm (8%), Ruffle membrane (4%), respiratory chain (4%), main axon (4%), myelin sheath (4%), and smooth endoplasmic reticulum (4%) (Figures5(c) and 5(d)).Most overlapping DEPs for the molecular functions classification were associated with syntaxin-1 binding (43%), calcium-dependent protein binding (15%), glutamate receptor binding (14%), clathrin binding (14%), and calmodulin binding (14%) (Figures5(e) and 5(f)).


### 3.6. Rhubarb Increased the Expression Level of ERK1/2 and Syn1

To confirm the results from the iTRAQ analysis, protein expression levels were evaluated for the 2 identified proteins Mapk1 (also means ERK2) and Synapse1 (Syn1) using Western blotting (n=4, each group) (Figure6). The results showed that the ratio of Syn1/β-tubulin and ERK1/2/β-tubulin were decreased in both brain tissues with CIR. Meanwhile, the protein expressions of Syn1 and ERK were increased significantly compared with MG. In comparison of SG, the increased protein relative expression of rhubarb group had not statistic difference in both Syn1 and ERK. In addition, quantitative protein levels by iTRAQ demonstrated the same trend for Syn1 and ERK in Western blotting.Figure 6
Western blotting analyses of Syn1 and ERK1/2. Proteins were examined for four times and normalized toβ-tubulin levels (loading control) for quantitative analysis. The levels of protein expression are presented as means ± SD. ∗p<0.05 compared with model.


## 3.1. Rats Treated with Rhubarb Show Attenuated Infarction Volume Rate and Improved Neurological Deficits

After 24 h reperfusion, 3 rats were dead and neurological assessments were conducted for the remaining rats. TTC staining visually showed the histological changes of cerebral tissues (Figure1(a)). SG appeared uniform red in color, indicating no infarction. Infraction volume of percentage showed a dramatic increase in MG (18.82% ± 5.23%) compared with SG (Figure 1(b)). Rhubarb significantly reduced acute infarction volume ratio (11.38% ± 4.04%) after 24 h reperfusion. As can be seen in Figure 1(c), obvious differences in regard to MCAO abnormality were exhibited between SG and MG (p<0.01). Rats treated with rhubarb showed a significant reduction of neurological deficits on mNSS compared with MG (p<0.05).Figure 1
(a) Representative TTC stained brain sections of sham group (n=3), model group (n=5), and rhubarb group (n=5) treated ischemic animals. (b) Infarct volume was determined 24 h after MCAO. (c) There were significant differences in neurological deficits at 24 h after MCAO between sham group (n=9), model group (n=10), and rhubarb group (n=8). Data are presented as mean ± SD. *p<0.05, **p<0.01, and ***p<0.001.
(a) (b) (c)


## 3.2. Quantitative Proteomic Analysis of Brains with MCAO Rats

In order to comprehensively understand the mechanism of rhubarb in ischemic stroke, iTRAQ analysis was performed on the protein extracts from rat brain with or without rhubarb treatment. After merging data from the respective replicates, 4578 proteins were identified and showed different expression (fold change ≥ 1.2, FDR ≤ 1%). Only the upregulated proteins with 1.5 ≤ an average ratio change ≤ 5 or downregulated proteins with an average ratio change ≤ 0.667. Peptides (95%) ≥ 2 and unused value ≥ 2 in both replicates were considered to be significantly expressed and used for the following analysis.In MG versus SG and RG versus MG, 287 proteins were differentially expression proteins (DEPs). There were 147 DEPs in RG versus MG and 140 DEPs in MG versus SG, respectively. Figure2(a) illustrated the proteins between the RG versus MG and MG versus SG by a Venn diagram. Only 64 DEPs were changed in MG versus SG uniquely, whereas 71 DEPs were regulated by rhubarb treatment alone. 76 overlapping DEPs showed altered protein expression level in both experimental sets and rhubarb group (Figure 2(b)). Table 1 summarized detailed information of 76 overlapping DEPs, including the UniProt accession number, gene name, peptides, fold change, and unused value. 53 overlapping DEPs were specifically upregulated, while 23 overlapping DEPs were downregulated after rhubarb treatment. Further, overlapping DEPs were subjected to functional analysis with reference to KEGG database.Table 1
Quantitative information of the overlapping DEPs in RG versus MG.
No. Uniprot Accession Gene name Unused % Cov Protein identity Peptides AVG. (95) (95%) 1 G3V7C6 Tubb4b 2 76.64 Tubulin beta-4B chain 274 2.43 2 Q9WTY2 Rem2 2.04 4.98 GTP-binding protein REM 2 2 1.88 3 Q04940 Nrgn 4.01 46.15 Neurogranin 7 3.81 4 P10818 Cox6a1 4.08 27.93 Cytochrome c oxidase subunit 6A1, mitochondrial 3 2.21 5 P04094 Penk 4.1 11.15 Proenkephalin-A 3 1.70 6 P22057 Ptgds 4.23 29.63 Prostaglandin-H2 D-isomerase 6 1.76 7 P35332 Hpcal4 6.4 49.21 Hippocalcin-like protein 4 12 2.15 8 P60841 Ensa 7.12 47.11 Alpha-endosulfine 7 2.99 9 Q3ZAV2 Ybx1 7.86 27.30 Nuclease-sensitive element-binding protein 1 4 1.60 10 G3V997 Dcx 8.55 20.82 Neuronal migration protein doublecortin 9 1.80 11 Q6AXS4 Atp6ap2 9.02 16.00 Renin receptor 8 1.73 12 P37377 Snca 9.85 80.71 Alpha-synuclein 22 2.37 13 O35796 C1qbp 10.07 28.32 Complement component 1 Q subcomponent-binding protein, mitochondrial 8 1.63 14 B7X6I3 Cend1 10.49 56.38 C38 protein 19 2.00 15 F8WFH6 Fam131b 11.76 20.28 Protein FAM131B 7 1.89 16 Q5U2U8 Bag3 12 17.42 Bcl2-associated athanogene 3 9 1.60 17 P35213 Ywhab 12.04 78.46 14-3-3 protein beta/alpha 65 1.92 18 O70257 Stx7 12.21 39.85 Syntaxin-7 13 3.47 19 A0A0G2K0M8 Ncam1 12.94 50.30 Neural cell adhesion molecule 1 77 2.44 20 P07825 Syp 13.95 40.07 Synaptophysin 19 1.55 21 P07171 Calb1 16.2 41.38 Calbindin 16 1.51 22 D3ZUX5 Chchd3 16.88 35.68 MICOS complex subunit 13 3.34 23 P62982 Rps27a 17.28 54.49 Ubiquitin-40S ribosomal protein S27a 23 1.62 24 D4A0T0 Ndufb10 18.05 61.93 Protein Ndufb10 15 2.43 25 Q09073 Slc25a5 18.08 61.07 ADP/ATP translocase 2 46 2.03 26 A0A0H2UHV6 Ppp3r1 18.44 61.25 Calcineurin subunit B type 1 17 1.55 27 D4A5L9 LOC679794 18.53 71.43 Protein LOC679794 24 3.31 28 P61589 Rhoa 19.6 64.77 Transforming protein RhoA 24 2.52 29 F1M378 Unc13a 20.94 6.01 Protein unc-13 homolog A 11 4.26 30 O88767 Park7 21.87 91.53 Protein deglycase DJ-1 22 1.76 31 D4AA63 Ubqln2 24.33 27.12 Protein Ubqln2 22 1.53 32 P14604 Echs1 24.5 53.45 Enoyl-CoA hydratase, mitochondrial 20 2.09 33 P08461 Dlat 27.5 35.60 Dihydrolipoyllysine-residue acetyltransferase component of pyruvate dehydrogenase complex, mitochondrial 25 3.38 34 P63086 Mapk1 31.61 49.16 Mitogen-activated protein kinase 1 26 1.63 35 Q6P6R2 Dld 33.08 54.03 Dihydrolipoyl dehydrogenase, mitochondrial 31 1.98 36 Q05175 Basp1 34.26 86.82 Brain acid soluble protein 1 82 4.18 37 P62161 Calm1 34.4 89.93 Calmodulin 80 2.64 38 P11275 Camk2a 35.04 64.44 Calcium/calmodulin-dependent protein kinase type II subunit alpha 76 1.73 39 P35565 Canx 35.34 37.39 Calnexin 31 1.50 40 P07936 Gap43 37.7 84.51 Neuromodulin 48 3.07 41 Q9Z0W5 Pacsin1 38.67 56.46 Protein kinase C and casein kinase substrate in neurons protein 1 50 1.61 42 P21707 Syt1 40.02 45.61 Synaptotagmin-1 37 1.79 43 Q5RKI0 Wdr1 42.43 46.37 WD repeat-containing protein 1 33 1.53 44 P63102 Ywhaz 43.94 85.31 14-3-3 protein zeta/delta 77 1.90 45 P06761 Hspa5 57.41 55.20 78 kDa glucose-regulated protein 57 1.55 46 P05708 Hk1 76.09 47.17 Hexokinase-1 80 1.50 47 Q63198 Cntn1 77.28 48.29 Contactin-1 77 1.74 48 D3ZMI4 Epb41l1 82.1 34.24 Band 4.1-like protein 1 52 1.77 49 P09951 Syn1 91.64 72.30 Synapsin-1 152 2.34 50 F1LQ81 Nsf 104.81 69.22 N-ethylmaleimide sensitive fusion protein, isoform CRA_b 86 1.84 51 D3ZQQ5 Dnm1 112.88 68.06 Dynamin-1 107 1.52 52 G3V984 Bsn 155.27 32.42 Protein bassoon 111 1.53 53 F1M779 Cltc 171.28 58.69 Clathrin heavy chain 213 1.76 54 A0A0G2K839 Clcn5 2 4.41 Chloride channel protein 3 0.02 55 G3V8L6 Kcna6 2.13 6.04 Potassium voltage-gated channel subfamily A member 6 2 0.65 56 D3ZV91 L3hypdh 4 7.34 Protein L3hypdh 2 0.48 57 F1LSV4 Sptlc2 4 5.18 Protein Sptlc2 3 0.52 58 A0A0G2K1A2 Mpo 4.01 3.76 Protein Mpo 3 0.65 59 P28572 Slc6a9 4.94 8.62 Sodium- and chloride-dependent glycine transporter 1 5 0.60 60 P55009 Aif1 5.28 23.13 Allograft inflammatory factor 1 5 0.60 61 Q6AYR8 Scrn2 5.8 13.00 Secernin-2 6 0.43 62 A0A0G2JWX4 Krt2 6.06 9.36 Keratin, type II cytoskeletal 2 epidermal 4 0.45 63 G3V6M3 Syt2 6.27 29.38 Synaptotagmin II 17 0.48 64 Q9R0I8 Pip4k2a 7.18 22.66 Phosphatidylinositol 5-phosphate 4-kinase type-2 alpha 14 0.46 65 Q7TP91 Surf4 10 9.55 Ab1-205 9 0.55 66 P0DMW1 Hspa1b 16.76 41.34 Heat shock 70 kDa protein 1B 73 0.51 67 P60203 Plp1 17.43 31.41 Myelin proteolipid protein 50 0.20 68 Q63345 Mog 18.16 36.73 Myelin-oligodendrocyte glycoprotein 14 0.39 69 A0A0G2JWM2 Sirt2 31.35 52.32 NAD-dependent protein deacetylase sirtuin-2 34 0.61 70 Q63357 Myo1d 31.36 21.27 Unconventional myosin-Id 27 0.48 71 E9PTX9 Slc12a2 34.28 21.28 Protein Slc12a2 26 0.44 72 P02688 Mbp 35.6 66.67 Myelin basic protein 88 0.23 73 F1LRZ7 Nefh 54.4 41.26 Neurofilament heavy polypeptide 55 0.47 74 P19527 Nefl 55.56 58.49 Neurofilament light polypeptide 87 0.32 75 P13233 Cnp 62.44 74.52 2′,3′-cyclic-nucleotide 3′-phosphodiesterase 83 0.47 76 P12839 Nefm 79.88 49.88 Neurofilament medium polypeptide 86 0.49Figure 2
The number of differentially expressed proteins (DEPs) and their overlap was exhibited in Venn diagram (a). 76 proteins with upregulated and downregulated expressions between MG versus SG and RG versus MG was showed in our study (b). Significant changes in the protein ratios were changed at least 1.5-fold.
(a) (b)


## 3.3. KEGG Database Analysis

The analysis was performed by matching the overlapping DEPs to the proteins annotated with KEGG pathway database via DAVID Bioinformatics Resources v6.8. The identified pathways were potentially affected by the modification of the abundance of proteins in the rat brain. Then 75 overlapping DEPs were identified and 1 overlapping DEP was not by DAVID. Among these overlapping DEPs, only 46 (56%) overlapping DEPs had a KEGG Orthology ID and were involved in 26 pathways. 14 pathways were statistically significant withp value < 0.05. The KEGG analysis results were shown in Figure 3 and Table 2. There are almost all closely related to Alzheimer's disease, oocyte meiosis, Parkinson's disease, oxytocin signaling pathway, cGMP-PKG signaling pathway, tuberculosis, synaptic vesicle cycle, amyotrophic lateral sclerosis (ALS), long-term potentiation (LTP), carbon metabolism, neurotrophin signaling pathway, prion diseases, amphetamine addiction, and glioma.Table 2
KEGG pathways associated with the 76 overlapping DEPs by DAVID.
Term Count PValue Genes rno04114: Oocyte meiosis 6 2.47E-04 P63102, P62161, P35213, A0A0H2UHV6, P63086, P11275 rno05010: Alzheimer’s disease 6 0.00251 P62161, P10818, D4A0T0, A0A0H2UHV6, P63086, P37377 rno05012: Parkinson’s disease 5 0.008091 O88767, Q09073, P10818, D4A0T0, P37377 rno04921: Oxytocin signaling pathway 5 0.009657 P61589, P62161, A0A0H2UHV6, P63086, P11275 rno04022: cGMP-PKG signaling pathway 5 0.011878 P61589, P62161, Q09073, A0A0H2UHV6, P63086 rno05152: Tuberculosis 5 0.014949 P61589, P62161, A0A0H2UHV6, P63086, P11275 rno05014: Amyotrophic lateral sclerosis (ALS) 4 0.002947 P19527, P12839, A0A0H2UHV6, F1LRZ7 rno04721: Synaptic vesicle cycle 4 0.004142 F1LQ81, P21707, F1M779, F1M378 rno04720: Long-term potentiation 4 0.004732 P62161, A0A0H2UHV6, P63086, P11275 rno01200: Carbon metabolism 4 0.025472 P05708, Q6P6R2, P14604, P08461 rno04722: Neurotrophin signaling pathway 4 0.028277 P61589, P62161, P63086, P11275 rno05020: Prion diseases 3 0.012162 P06761, A0A0G2K0M8, P63086 rno05031: Amphetamine addiction 3 0.044482 P62161, A0A0H2UHV6, P11275 rno05214: Glioma 3 0.045744 P62161, P63086, P11275Figure 3
KEGG pathway analysis of overlapping DEPs identified by iTRAQ analysis.


## 3.4. Protein-Protein Interaction (PPI) Network Analysis

Further PPI network analysis was conduct by Cytoscape. PPI network (Figure4) showed that several significant proteins might function in response to ischemic stroke. Of these 75 identified overlapping DEPs, a total of 52 were linked into a protein interaction network and 23 did not show any link. The network analysis revealed potential relationship between Mapk1, Syn1, Calm1, Syp, Gap43, Hspa5, Mbp, and so on. The network view summarized the network of predicted associations for particular proteins. It provided us a crucial modeling scaffolds and data reduction way which made us got insight into the potential mechanisms of rhubarb in brain.Figure 4
Protein network connectivity of identified proteins in MCAO rat with or without rhubarb.


## 3.5. GO Analysis

To detect the potential role of rhubarb in MCAO rats, we allocated 76 overlapping DEPs to biological process, cell component, and molecular function based on GO annotation. Biological processes that overlapping DEPs participated were representatively involved in synaptic vesicle cycle (28%), neurotransmitter transport (22%), neurofilament bundle assembly (16%), presynaptic process involved in chemical synaptic transmission (11%), oxidative phosphorylation (10%), mitochondrial membrane organization (5%), substantia nigra development (3%), negative regulation of transcription from RNA polymerase II promoter in response to stress (3%), and response to monoamine (2%) (Figures5(a) and 5(b)).Figure 5
GO annotation analysis of the overlapping DEPs. The color had a one-to-one correspondence between (a), (c), (e) and (b), (d), (f), respectively. ((a) and (b)) GO annotation enriched by ClueGO according to biological process. ((c) and (d)) GO annotation enriched by ClueGO according to cellular components. ((e) and (f)) GO annotation enriched by ClueGO molecular function. Functionally grouped network with terms linked by nodes, and functionally related groups were partially overlapped; the node size represented the term enrichment significance.
(a) (b) (c) (d) (e) (f)Cell component suggested that proteins were located-representatively in the axon (48%), growth cone (8%) neurofilament (8%), melanosome (8%), cell projection cytoplasm (8%), Ruffle membrane (4%), respiratory chain (4%), main axon (4%), myelin sheath (4%), and smooth endoplasmic reticulum (4%) (Figures5(c) and 5(d)).Most overlapping DEPs for the molecular functions classification were associated with syntaxin-1 binding (43%), calcium-dependent protein binding (15%), glutamate receptor binding (14%), clathrin binding (14%), and calmodulin binding (14%) (Figures5(e) and 5(f)).


## 3.6. Rhubarb Increased the Expression Level of ERK1/2 and Syn1

To confirm the results from the iTRAQ analysis, protein expression levels were evaluated for the 2 identified proteins Mapk1 (also means ERK2) and Synapse1 (Syn1) using Western blotting (n=4, each group) (Figure6). The results showed that the ratio of Syn1/β-tubulin and ERK1/2/β-tubulin were decreased in both brain tissues with CIR. Meanwhile, the protein expressions of Syn1 and ERK were increased significantly compared with MG. In comparison of SG, the increased protein relative expression of rhubarb group had not statistic difference in both Syn1 and ERK. In addition, quantitative protein levels by iTRAQ demonstrated the same trend for Syn1 and ERK in Western blotting.Figure 6
Western blotting analyses of Syn1 and ERK1/2. Proteins were examined for four times and normalized toβ-tubulin levels (loading control) for quantitative analysis. The levels of protein expression are presented as means ± SD. ∗p<0.05 compared with model.


## 4. Discussion

In this study, we first confirmed the effective treatment of rhubarb which significantly reduced infarction volume ratio and neurological assessments on MCAO model rats. Next, we performed a quantitative proteomics and bioinformatics analysis to explore the disease-related protein abundance. Further comparing RG with MG, Western blotting identified the changes of candidate proteins which were in accord with the results of proteomics. Our study may suggest that rhubarb is expected to be a novel drug for ischemic stroke.To the best of our knowledge, few studies reveal the mechanism of rhubarb on acute ischemic stroke. To address this concern, the iTRAQ-based proteomics analysis was utilized to exploit multiple targets and pathways following ischemic stroke in rats model of MCAO. As a result, we identified a total of 4578 proteins, of which 287 were DEPs. Among the three groups, there were 76 overlapping DEPs. 76 overlapping DEPs were part of altered protein after cerebral ischemia and rhubarb treatment interfered with these 76 proteins at the same time. Furthermore, there were 53 upregulated proteins and 23 downregulated proteins after rhubarb treatment. Then for these 76 overlapping DEPs a series of bioinformatics analysis were conducted. By GO annotation, it involves a variety of biological processes, including synaptic vesicle cycle and neurotransmitter transport, among others both in MG and RG.In KEGG pathways analysis, it exhibited that cerebral ischemia and rhubarb treatment involve oxytocin signaling pathway, cGMP-PKG signaling pathway, synaptic vesicle cycle, and LTP in top 14 pathways. Oxytocin (OT) which is a typical hormone responds to many pathological conditions following ischemia-reperfusion. According to altered expression patterns of GABAA receptor subunit and the kinetics of GABA-induced chloride ion inﬂux, OT conferred neuroprotection against ischemic stroke [32]. OT also could inhibit the calpain to improve morphology damage in neurons after ischemic stroke [33]. Through the activating cGMP-dependent protein kinas (PKG) and resultant phosphorylation of cellular proteins, cGMP play a role of a second messenger which exerts most of its action in central nervous system [34]. The activation of the cGMP-PKG signaling pathway enhanced LTP at thalamic inputs to the lateral amygdala (LA). It was suggested that synaptic plasticity and fear memory consolidation might be promoted via NO-cGMP-PKG in the LA and in part by activating the ERK/MAPK signaling cascade [35]. Rhubarb may exert a multitarget and multipathway regulation mechanism through these overlapping DEPs in ischemic stroke.Therefore, we chose Mapk1 and Syn1 for confirming the reproducibility of our studies and the excellent quality of our data. In iTRAQ analysis and the results of Western blotting, Mapk1 and Syn1 were obviously upregulated after rhubarb treatment. Mapk1 (also known as MAPK1, p42MAPK) and Syn1 have important roles in synaptic transmission and plasticity [36, 37].The synaptic damage and plasticity could affect the development of treatments to facilitate rehabilitation of patients with ischemic cerebral damage [38]. Animal evidence indicates that it opens a time-limited window of neuroplasticity following a stroke which is the greatest gains in recovery [39].Mapk signaling cascade is critical for both the memory consolidation and the long-term neuronal plasticity. ERK1 and ERK2 are activated in neurons in response to excitatory glutamatergic signaling, which controls many forms of synaptic plasticity [40–42]. Studies also demonstrate that ERK signaling plays an important role in the induction of LTP which is an activity-dependent strengthening of synaptic efficacy [43, 44]. ERK1/2 has also been shown to phosphorylate Syn1 which is a major substrate in nerve terminals [45]. Syn1, one of the synapsins family, is synaptic vesicle phosphoproteins regulating neural development, plasticity, and synaptic transmission [37, 46]. The ERK1/2-dependent phosphorylation of Syn1 is stimulated hippocampus-dependent learning. The increase of neurotransmitter release selectively in LTP facilitated enhancements in hippocampus-dependent learning [47]. Furthermore, Syn1 is required by ERK1/2-dependent synaptic plasticity [48]. Knocking out Syn1 will inhibit the enhancements of learning, presynaptic plasticity, and LTP [49]. Study also shows that Syn1 is the main presynaptic effector for certain forms of LTP triggered by the activation of the H-ras-Mapk pathway [43]. By activating the ERK1/2-Syn1 signaling pathway, the collagen nanofibrous scaffolds promoted to the presynaptic maturation [50]. Our Western blotting results suggested that ERK1/2 and Syn1 were significantly increased by the treatment of rhubarb compared with the model group. There are robustly associated with the synaptic vesicle cycle and presynaptic process involved in chemical synaptic transmission regulated by rhubarb following ischemic stroke. Thus, it is concluded that rhubarb may be essential for the development of synaptic plasticity via ERK1/2-Syn1 signal pathway.Rhubarb is one of the important Chinese herbal medicines. Its major active components including rhein and chrysophanol are the anthraquinones which exist in combined (glycosides) and free (aglycones) forms. The rhubarb and its extracts are involved in decreasing the content of NO in brain tissue and antianoxia on ischemic stroke [51, 52]. The study reports that anthraquinones present protective effect against focal cerebral ischemia via resisting the aggregation and adhesion of platelet or decreasing plasminogen and fibrinogen [53]. Rhubarb aglycone presents protective effect according to improving cerebral infraction and neuronal apoptosis in cerebral ischemia/reperfusion [54]. Chrysophanol, an active component of rhubarb aglycones, protects the damage of brain tissues following ischemic stroke via inhibiting the activation of inflammasome [55]. These studies and our research support that rhubarb may represent a potential neuroprotection modulating multiple proteins and pathways for ischemic stroke. In a further study, we should investigate more upstream or downstream molecular mechanisms from the DEPs in detail.


## 5. Conclusions

The present study reveals that rhubarb highlights a promising neuroprotective effect on acute MCAO rats. Rhubarb exerts the treatment via modulating multiple proteins, targets, and pathways in a rat model of ischemic stroke. By the verification of molecular biology technique, the rhubarb exhibits essential impact on synaptic plasticity. Our study might provide a novel therapeutic agent of cerebral ischemic stroke.



---

*Source: 6920213-2018-07-16.xml*