Transcriptome profiling revealed early vascular smooth muscle cell gene activation following focal ischemic stroke in female rats – comparisons with males
BMC Genomics volume 21, Article number: 883 (2020)
Women account for 60% of all stroke deaths and are more often permanently disabled than men, despite their higher observed stroke incidence. Considering the clinical population affected by stroke, an obvious drawback is that many pre-clinical and clinical studies only investigate young males. To improve therapeutic translation from bench to bedside, we believe that it is advantageous to include both sexes in experimental models of stroke. The aims of this study were to identify early cerebral vascular responses to ischemic stroke in females, compare the differential gene expression patterns with those seen in males, and identify potential new therapeutic targets.
Transient middle cerebral artery occlusion (tMCAO) was used to induce stroke in both female and male rats, the middle cerebral arteries (MCAs) were isolated 3 h post reperfusion and RNA was extracted. Affymetrix whole transcriptome expression profiling was performed on female (n = 12) MCAs to reveal differentially expressed genes. In total, 1076 genes had an increased expression and 879 genes a decreased expression in the occluded MCAs as compared with the control MCAs from female rats. An enrichment of genes related to apoptosis, regulation of transcription, protein autophosphorylation, inflammation, oxidative stress, and tissue repair and recovery were seen in the occluded MCA. The high expression genes chosen for qPCR verification (Adamts4, Olr1, JunB, Fosl1, Serpine1, S1pr3, Ccl2 and Socs3) were all shown to be upregulated in the same manner in both females and males after tMCAO (p < 0.05; n = 23). When comparing the differentially expressed genes in female MCAs (occluded and non-occluded) with our previous findings in males after tMCAO, a total of 297 genes overlapped (all groups had 32 genes in common).
The cascades of processes initiated in the vasculature following reperfusion are complex. Dynamic gene expression alterations were observed in the occluded MCAs, and to a less pronounced degree in the non-occluded MCAs. Dysregulation of inflammation and blood-brain barrier breakdown are possible pharmacological targets. The sample of genes (< 1% of the differentially expressed genes) validated for this microarray did not reveal any sex differences. However, sex differences might be observed for other gene targets.
Ischemic stroke is one of the leading causes of death and disability in the world . Although the incidence is higher in men, women account for 60% of all stroke deaths and are more often permanently disabled than men . Thrombolysis, the only available non-invasive treatment for stroke, has in some studies been shown to have a better effect in women than in men [3, 4]. Although reperfusion by thrombolysis or thrombectomy has been shown effective in salvaging neurological function, restoration of blood flow and reduction of damages to the blood-brain barrier (BBB) increases the risk of hemorrhagic transformation and edema, which may be potentially fatal complications . In the case of thrombolysis, the risk of these adverse effects does not outweigh the benefits past 4.5 h post stroke, which limits the use in clinical practice .
Despite intense research efforts during several decades with more than 1000 compounds tested and numerous interventions that have shown promise in pre-clinical studies, all failed in the clinical studies . Some of the main reasons proposed for this “translational roadblock” in stroke treatment are related to the fact that the majority of pre-clinical studies have been performed in young healthy male rodents, a clear drawback when considering the clinical population affected by stroke. By including both sexes in experimental models of stroke, it may be possible to more accurately represent the clinical scenario and thus improve therapeutic translation from bench to bedside. Despite increasing awareness of the importance of sex differences, the majority of pre-clinical and clinical studies are still performed on males .
Developing effective treatment strategies for both men and women requires a deeper understanding of sex differences in the underlying mechanisms of ischemic injury. In experimental stroke models, female animals have smaller ischemic areas and better functional outcomes, and this difference is nullified by ovariectomy, suggesting that female sex hormones (estrogen and progesterone) are responsible . The protective effect of estrogen has been shown to be multifactorial, acting on both the vasculature and neurons . Sex hormones do not fully account for all sex differences, it has also been demonstrated that neuronal apoptosis pathways differ between males and females  and that male neurons are also more sensitive towards nitrosative stress . We have demonstrated that there are differences in the cerebrovascular receptor expression in males versus females both in human brain vessels and in rats after transient middle cerebral artery occlusion (tMCAO, a standard method for this type of experimental stroke) [12, 13].
Ischemic stroke is primarily a vascular disease and we hypothesize that reperfusion and subsequent protection of the brain against hemorrhage, inflammation and edema by targeting the cerebral arteries is the first step towards successful stroke treatment . The BBB consists of endothelial cells with continuous tight junctions, which offer protection against the pathogens, toxins and reduce the influence of the peripheral immune system in the brain. The endothelial cells are supported by the vascular smooth muscle cells (VSMCs), astrocytes, pericytes and the extracellular matrix (ECM). This system is disrupted after stroke due to formation of reactive oxygen species (ROS) and subsequent inflammatory processes . This allows peripheral inflammatory cells to migrate across the BBB and cause further destruction to the brain tissue. In addition to endothelial damage, reperfusion puts a strain on the VSMCs, causing enhanced vasocontractile responses which reduce perfusion . The VSMCs also express inflammatory cytokines in response to ischemia-reperfusion, including metalloproteinases, which contribute to recruitment of inflammatory cells and further BBB breakdown . This has been verified ex vivo in both rodents following experimental stroke [16, 17] and humans .
The present study was conducted to examine the early cerebrovascular processes of vascular damage after stroke in females and subsequently examine if sex differences and similarities in these responses occurring in the cerebral vessel wall exist. After a stroke, there is enhanced activation of phosphorylated extracellular signal-regulated kinase 1 and 2 (pERK1/2) in the cerebral vasculature already after a few minutes which reach even higher levels at 3 h . In order to examine which genes were activated in the early stroke stages (at 3 h), we performed whole-transcriptome expression profiling on middle cerebral arteries (MCAs) of female rats after tMCAO-induced ischemia. This was also done to identify activated biological processes and pathways locally in the MCAs which potentially could be targeted for vascular protection after stroke. In addition to the microarray, we validated selected high-expression genes that potentially are involved in reperfusion injury in male and female rats, and compared the differentially expressed genes in MCAs from females (current study) with MCAs from males (previous study, Grell et al. ) after tMCAO to contribute to basic knowledge of vascular wall processes in both sexes after stroke.
Physiological parameters for both sexes
The body weight was significant lower in females (243 ± 9 g) than males (323 ± 16 g) (p < 0. 05), although they were of the same age. During the occlusion and reperfusion, the blood flow over the MCA distribution area was measured with a laser Doppler flowmetry probe [21, 22]. Insertion of the intraluminal filament resulted in a mean reduction of blood flow by 73 ± 10% in females and by 76 ± 9% in males (p > 0.05). Withdrawal of the filament after two hours of occlusion resulted in a mean increase of blood flow by 65 ± 13% in females and by 49 ± 17% in males (p > 0.05). This resulted in a blood flow approximate to the level observed prior to the occlusion. Prior to the occlusion, body temperature, blood pressure, blood glucose, pH, paCO2 and paO2 were measured. These parameters were within the physiological range and did not differ between the sexes (data not shown).
Whole transcriptome expression profiling in females
In the microarray analysis, bilateral MCA segments (occluded (nsample = 5) and non-occluded (nsample = 6)) from 6 female rats subjected to tMCAO were analyzed. The arteries were removed 3 h post-reperfusion. MCAs from 6 healthy female rats were used as controls. In total, 1076 genes showed an increased expression and 879 genes had a decreased expression in the occluded MCAs as compared with the control MCAs. In the contralateral non-occluded MCAs, 111 genes had an increased expression and 92 genes had a decreased expression. The expressions of 80 of the differentially expressed genes were increased in both the occluded and non-occluded MCAs while the expressions of 67 of the differentially expressed genes were decreased in both the occluded and non-occluded MCAs (Fig. 1, Fig. 2).
Gene and protein ontology enrichment analysis in females
Occluded MCAs compared with control arteries
With the significantly differentially expressed genes in the microarray, a gene ontology (GO) enrichment analysis was performed to identify activated biological processes. When comparing the occluded MCAs to control MCAs, 91 GO terms within the biological process domain were overrepresented (the top findings are presented in Table 1). Examples of other overrepresented GO terms of interest than those presented in Table 1 include: ‘Regulation of epithelial cell migration’ (GO:0010632; annotation of 47 differentially expressed genes), ‘Regulation of neuron projection development’ (GO:0010975; annotation of 91 differentially expressed genes) and ‘Regulation of cytoskeleton organization’ (GO:0051493; annotation of 83 differentially expressed genes). Protein ANalysis THrough Evolutionary Relationships (PANTHER) and Reactome enrichment analyses were also carried out. The genes differentially expressed in the occluded MCAs in relation to control MCAs were overrepresented for 17 PANTHER protein classes, 9 PANTHER pathways and 9 Reactome pathways (Fig. 3a).
Across analyses, an enrichment of genes related to apoptosis, regulation of transcription, protein autophosphorylation, inflammation, oxidative stress, and tissue repair and recovery could be seen (Table 1 and Figs. 3, 4 and 5). Within the overrepresented PANTHER protein class ‘Intercellular signal molecule’ (PC00207), chemokine (C-C motif) ligand 2 (Ccl2) had the highest fold change (Fig. 3b). Ccl2 is also annotated to the overrepresented GO terms ‘Cellular response to interleukin-1’ (GO:0071347) and ‘Cellular response to tumor necrosis factor’ (GO:0071356) amongst other cytokines (Fig. 5). Within the overrepresented Reactome pathway ‘Cytokine signaling in immune system’ (R-RNO-1280215), suppressor of cytokine signaling 3 (Socs3) had the highest fold change (Fig. 3c). This was also the case amongst the kinase modulators (Table 2).
Occluded and non-occluded MCAs compared with control arteries
The number of differentially expressed genes after experimental stroke was considerable higher in the occluded MCAs (1955 differentially expressed genes) than in the non-occluded MCAs where only 203 differentially expressed genes were identified (Fig. 1). To reveal if the biological processes activated in the occluded MCAs are similar to those activated in the non-occluded MCAs, we looked for overlapping overrepresented GO terms. When using the predefined selection criteria, no overlap was observed. When no predefined selection criteria were applied, overlap was observed between the two experimental groups. A selection of overlapping overrepresented GO terms can be found in Table 3 (relaxed criteria). In addition to the GO enrichment analysis, 2 PANTHER pathways (‘Angiogenesis’ (P00005) and ‘Toll receptor signaling pathway’ (P00054)) and 3 PANTHER protein classes (‘Protein-binding activity modulator’ (PC00095), ‘Metabolite interconversion enzyme’ (PC00262) and ‘Gene-specific transcriptional regulator’ (PC00264)) were overrepresented in the list of differentially expressed genes for both the occluded and non-occluded MCAs when compared with control MCAs (Fig. 3a).
Focusing on the gene-specific transcriptional regulators, 108 and 15 differentially expressed genes in the occluded and non-occluded MCAs were annotated to this GO term, respectively. Of these, twelve differentially expressed genes were found for both experimental groups (e.g. PR domain zinc finger protein 4 (Prdm4), Runt-related transcription factor 1 (Runx1) and signal transducer and activator of transcription 3 (Stat3)). Focusing on the occluded MCAs, 94 of the gene-specific transcriptional regulators were DNA-binding transcription factors (PC00218). Ten of those can more specifically be categorized as basic leucine zipper transcription factors (PC00056, Table 2). For the 108 gene-specific transcriptional regulators, we identified 3 protein-protein interaction clusters formed by 32 of these regulators (e.g. Runx1, Fos-like antigen 1 (Fosl1) and Jun B proto-oncogene (JunB); Fig. 4a). For each cluster, we identified overrepresented GO terms within the biological process domain. The regulators forming cluster 2 (Fosl1 and JunB is a part of this cluster) are involved in transcription and cytokine response (Fig. 4b).
qPCR for validation of target genes for both sexes
To validate the results from the microarray, eight high-expression genes of interest were chosen for quantitative real-time polymerase chain reaction (qPCR) analysis in MCAs from both female and male rats (a new set of animals were operated on in both sexes, 5 stroke females and 6 stroke males were included in the analysis as well 6 controls of each sex). The expression of the target genes in the occluded MCAs, non-occluded MCAs and control MCAs from the microarray analysis can be found in Fig. 2 and Fig. 6a. To get an increased understanding of what biological processes the target genes are involved in, we categorized them into preselected GO terms. The target genes are all involved in response to stress. Other GO terms of interest include cell communication, defense response and response to cytokine (Fig. 6b).
For qPCR, glyceraldehyde 3-phosphate dehydrogenase (Gapdh) and actin B (ActB) were used as reference genes; an equal stable high expression was confirmed throughout the groups (data not shown). All of the 8 analyzed genes (Ccl2, oxidized low-density lipoprotein receptor 1 (Olr1), a disintegrin and metalloproteinase with thrombospondin type 1 motif, 4 (Adamts4), serine protease inhibitor, clade E, member 1 (Serpine1), sphingosine 1 phosphate receptor 3 (S1pr3), Socs3, JunB and Fosl1) were significantly upregulated in the occluded MCAs compared with control MCAs (Figs. 7 and 8). In addition, Ccl2, Socs3, Fosl1, JunB and Serpine1 were also upregulated in the non-occluded MCAs as compared to control MCAs. Sex did not have a significant effect on the expression of any of the 8 analyzed genes (p = 0.11–0.87, Figs. 7 and 8).
Cross-analysis to reveal sex similarities
The comparison with findings presented by Grell et al.  revealed overlap between the differential expressed genes identified in the occluded versus non-occluded MCAs from Wistar Kyoto (WKY) male rats, the occluded MCAs versus control MCAs from Wistar female rats and the non-occluded MCAs versus control MCAs from Wistar female rats. Thirty-two genes were differentially expressed in each of the experimental groups (Fig. 9). An overview of those genes is presented in Table 4.
This is to our knowledge the first study to explore possible mechanisms of reperfusion injury after stroke in females by performing a whole transcriptome microarray on isolated MCAs 3 h after the end of a 2-h tMCAO. This unique approach has allowed us to identify potential targets for treatment of vascular injury associated with stroke. The tMCAO model was used to mimic stroke conditions and the cerebral arteries were examined at an early time point to explore the events that occur in VSMCs and endothelial cells that lead to BBB breakdown and endothelial dysfunction. At a similar early time point, we have demonstrated that there is increased expression of phosphorylated protein kinase C-δ (pPKCδ) and pERK1/2 but not of phosphorylated p38, c-jun N-terminal kinase (pJNK) and pPKCα . This has been verified by protein studies  and by specific blockade, using inhibitors of the MEK-ERK pathway which resulted in reduced brain damage . Within this study, whole transcriptome expression profiling revealed a considerable number of differentially expressed genes (more than 2000) in the MCAs from female rats after tMCAO. Interestingly, almost 300 of these genes were also differentially expressed in MCAs from male rats after tMCAO.
The present study was also designed to examine if sex differences affect early gene expression patterns of stroke-associated genes. From the microarray analysis, we found a markedly increased expression of Adamts4, Socs3, Fosl1, JunB, Serpine1, Olr1, Ccl2 and S1pr3 in the ischemic artery compared to non-ischemic artery at 3-h post reperfusion in female rats. In a further validation of the genes we compared both sexes and found a similar degree of increase in expression. Thus, the results of this part of the study confirmed an increased expression of genes in the vasculature that contributes to inflammation, endothelial dysfunction and breakdown of the BBB.
Adamts4 is of interest since it for the first time is shown to be upregulated in the cerebral arteries in response to stroke. Adamts4 is a metalloproteinase that has been shown to have a regulatory function in inflammation and an ability to inhibit angiogenesis [24, 25]. Other metalloproteinases were also found to be overexpressed in the microarray analysis of the occluded MCA (matrix metallopeptidase 19 (Mmp19), Adamts1, Adamts5, Adamts9, ADAM metallopeptidase domain 8 (Adam8) and Adam9), but not to the same degree as Adamts4. Adamts4 has been shown to promote neurite-outgrowth through the MEK pathway, independent of its proteolytic capacities . It is highly expressed in VSMCs in macrophage-rich areas in the thoracic aorta of patients with thoracic aortic dissection and increased serum levels have been associated with carotid atherosclerotic plaque vulnerability [27, 28]. Based on knowledge about the proteolytic capacities of Adamts4 from earlier studies, one could hypothesize that Adamts4 and other metalloproteinases contribute to breakdown of the supporting ECM of the arteries, negatively affecting the integrity of the cerebral arteries, increasing the risk of hemorrhagic transformation and edema . Several of the metalloproteinases and ECM genes were upregulated at 6 and 24 h after subarachnoid hemorrhage (SAH) , indicating a process that is ongoing. The time course of the Mmp8, Mmp9 and Mmp13 were carefully monitored in the MCAs immediately after the SAH up until 48 h in male rats .
Olr1 and Ccl2 expression have earlier been shown to be upregulated in cerebral arteries following tMCAO in both normotensive and hypertensive male rats . Here, we report that this occurs equally well in male and female MCAs after experimental stroke. The Olr1 gene codes for a scavenger receptor located on endothelial cells, VSMCs and macrophages and are one of the key molecules in atherogenesis. Oxidized low-density lipoprotein activates the receptor, leading to recruitment of inflammatory cells by increasing chemokine and adhesion molecule levels. Although this process has been well characterized in the settings of developing atherosclerosis, it has not yet been established what role Olr1 plays in the acute phase after ischemia-reperfusion and in particular in the affected MCAs.
Ccl2 is a chemokine which recruits monocytes and basophils. Monocytes in turn, secrete further inflammatory cytokines, aggravating the damage further. Clot formation, especially fibrin, triggers activation of monocytes  and endothelial cells . Previous studies have shown that the expression of Ccl2 is highly increased in the brain after stroke, and that silencing the Ccl2 gene is protective in stroke models . Here we show that the expression is highly elevated in both sexes early after reperfusion which formerly has been seen in hypertensive males 6 h post reperfusion but never shown to be equally increased at this early time point after a tMCAO or in both sexes. We have previously reported on the increased presence of the chemokine (C-X-C motif) ligand 2 (Cxcl2) and interleukin-6 (Il-6) at gene level in experimental SAH at 24 h , and followed the time course from 0 to 48 h of several other cytokines such as tumor necrosis factor α (TNFα) and Il-1β, and chemokine ligands Cxcl1 and Ccl20, in addition to Il-6 and Cxcl2 . Again, this clearly shows that the process is ongoing in the cerebral vessel walls at least 48 h after experimental SAH.
Another gene of interest is S1pr3 which has an increased expression after ischemia-reperfusion in cerebral arteries. Sphingosine 1-phosphate (S1P) signaling via S1pr3 has been shown to have both protective effects on endothelial function in ischemia/reperfusion in the heart  as well as vasocontractile properties when located on the VSMCs in cerebral arteries . Activation of S1pr3 has also been shown to negatively influence endothelial barrier function .
Transcription factors and intracellular signaling pathways are most likely more effective treatment targets when it comes to inflammation as opposed to blocking individual specific cytokines or chemokines. Inflammation is crucial for repair and recovery after stroke, but in the initial phase, reduction of inflammation has been shown to be beneficial. By targeting the arteries, we believe that we can reduce the initial recruitment and migration of inflammatory cells. Our data showed that the selected inflammatory regulatory genes Socs3, JunB and Fosl1 were overexpressed in both the occluded and non-occluded artery compared to healthy controls, indicating that inflammation is not exclusive to the occluded artery. JunB and Fosl1 are part of the activator protein complex, which make up the transcription factor activator protein 1 (AP-1). AP-1 has been suggested to be responsible for activating inflammation and atherogenesis in response to changes in sheer stress, which could explain why these genes here are shown to be upregulated early after reperfusion . The upstream activator of the AP-1 complex is the stress-activated protein kinases/Jun amino-terminal kinases (SAPK/JNK), which has earlier been shown to be activated in cerebral VSMCs 24 h after SAH, tMCAO and ex vivo organ culture of cerebral arteries . In experimental SAH, we have additionally observed the dynamic activation of two other transcription factors, activating transcription factor–2 (ATF-2) and ETS transcription factor (ELK-1), during the first two days after the induction of the hemorrhage in male rats .
The current study focuses solely on the early cerebrovascular processes of vascular damage after a tMCAO-induced stroke. Investigation of long-term effects after a tMCAO-induced stroke, as well as stroke induced by other experimental procedures, would be beneficial on the path to discover new therapeutic targets and to improve our understanding of the stroke pathophysiology.
This study may pave the way towards novel ways to modify cerebrovascular gene regulation and provide novel therapeutic targets to improve cerebrovascular disorders like stroke. Our findings revealed that following tMCAO there are dynamic alterations in gene expression in the immediate MCA, but also on the contralateral side, albeit to a less pronounced degree. Focusing on the high expression gene targets, we observed no significant sex-dependent difference (pattern of changes in males was similar to those seen in females). The genes activated at this early 3 h time point suggested gene activation predominantly involved in transcription, BBB and endothelial function, and in inflammation . Our study was exploratory and designed to examine if there are sex similarities and differences in gene expression. The study on females provides a better knowledge base for further research as well as increasing the chance of novel findings. Examining early gene expression is just a minor step in understanding the complex cascade of events that are initiated by ischemia-reperfusion , but we believe that it can provide a base for further investigation of differential gene and protein expression at later time points.
We aimed at identifying early cerebral vascular responses to ischemic stroke in female rats and identifying potential new therapeutic targets. We also aimed at comparing the differential gene expression patterns observed in MCAs from female rats with those seen in males.
Overview of study design and inclusion and exclusion criteria
The study was designed as a case-control study with the occluded and non-occluded MCAs (after tMCAO) as case groups and the MCAs from healthy rats (non-operated) as control group, and the primary study outcome was differential gene expression (gene expression in cases relative to controls). The sample size was set based on experience from earlier studies conducted in the lab. The experimental unit was individual animal, only animals surviving the operation until the 3 h’ time point were included and the confounding factors age, sex, estrous cycle and environmental factors were controlled for. The confounding factor sex was investigated as an explanatory variable. To reduce pain, suffering and distress, animals were operated under full anesthesia and sacrificed under sedation. The following exclusion criteria were applied during animal experiments: (i) Not satisfactory occlusion or reperfusion, (ii) no reduction in neuroscore, or (iii) death during or after the procedure; and during microarray or qPCR: (i) Not satisfactory purity or concentration of RNA. All inclusion and exclusion criteria were specified prior to the experiments, except for the satisfactory occlusion threshold which was specified during the experiments. Ten female rats and 6 male rats were excluded due to death during the procedure or failing to meet the exclusion/inclusion criteria. For the microarray, three samples were excluded due to poor RNA integrity and 1 sample was excluded due to poor amplification. For the qPCR, seven samples were excluded at this stage due to low RNA-concentration and 3 samples due to failed amplification (Fig. S1).
12-week old female (n = 42) and male (n = 18) Wistar rats were obtained from Charles River (Charles River laboratories, Sulzfeld, Germany). After applying the inclusion and exclusion criteria, 23 female and 12 male rats remained (Fig. S1). Animals were housed in environmental enriched cages with a 12 h light/dark cycle (lights on 7 am-7 pm), climate- and humidity- controlled environment and with free access to food and water in the University of Lunds animal facility. Prior to initiating the experimental procedures, animals were subjected to an acclimatization period. All experimental procedures were performed in accordance with the ARRIVE guidelines, the ethical guidelines of the International Association for the Study of Pain regulations on animal welfare and the National Institutes of Health guidelines for the care and use of laboratory animals. The experimental procedures have been previously approved by the Institutional Animal Care and Use Committee of the University of Lund. The rats were randomly divided into control and tMCAO experiments. The female rats were monitored daily with vaginal smears for a minimum of two consecutive cycles prior to sacrifice to determine the estrous phase. Samples were collected with a saline soaked cotton swab, transferred onto microscope slides, air-dried and stained with hematoxylin-eosin. The samples were collected at the same time-point each day. Female rats that were in the high estrogen phase, proestrus, were not included to reduce variation due to hormone fluctuations. The methodology has been further described by Goldman et al. .
Transient middle cerebral artery occlusion
MCA occlusion was performed by using the intraluminal filament technique . Anesthesia was induced by 3.5% isoflurane in N2O:O2 (70:30) and maintained by continuous inhalation of 2.5% isoflurane in N2O:O2 (70:30). An arterial tail catheter was inserted to monitor blood gases and mean arterial blood pressure during the occlusion. A rectal thermometer connected to a homoeothermic blanket was used to maintain body temperature at 37 °C during the procedure. A laser Doppler probe (Perimed, Järfälla, Sweden) was fixed to the skull 6 mm laterally of the midline and 1 mm posterior to bregma. An incision was made in midline of the neck, exposing the carotid artery. The external carotid and common carotid were permanently ligated with sutures. A rubber-coated monofilament (Doccol Corporation, Redlands, CA, USA) was inserted through an incision in the common carotid artery through the internal carotid artery until a sudden drop in cerebral blood flow was observed in the area supplied by the MCA (as measured with Laser Doppler flowmetry). Animals with a minimum flow reduction of 50% were included in the study. The filament was secured by sutures and the surgical areas closed while the anesthesia was discontinued, and the animal allowed recovering. The method has been described in details before [21, 22].
Prior to reperfusion, the rats were evaluated with the 6-point neuroscore [39, 40]. Only animals showing a score of 3–4 in the test were re-anesthetized after 120 min and the filament removed, which resulted in an increase in cerebral blood flow (as measured with Laser Doppler flowmetry). Animals with a minimum flow increase of 30% were included in the study. The animals were then allowed to recover with free access to food and water for 3 h.
Middle cerebral artery isolation
Three hours post reperfusion; the animals were euthanized with carbon dioxide and decapitation. The brain was immediately removed and both the right and left MCAs (distal part; length = 5 mm, diameter 0.2 mm) were carefully dissected out. The MCAs were carefully cleaned from surrounding connective tissue and blood, frozen on dry ice and stored in − 80 °C. Inclusion criteria were successful tMCAO operation as judged from the Laser Doppler analysis of regional cerebral blood flow (rCBF) and neuroscore evaluation. All animals fulfilling this survived the 3 h' until sacrifice. The subsequent handling and analysis of the removed MCA segments were semi-blinded meaning that the tissue and group were unknown to the analyst from this time point until the time where the data were analyzed.
The extraction of RNA was performed with the same method for both the microarray and PCR. The RNA was isolated with the Nucleospin miRNA isolation kit (Machery-Nagel, Düren, Germany), following the manufactures instructions for extraction of total RNA. The artery-samples were first homogenized in in Lysing matrix D tubes containing 1.4 mm ceramic spheres (MP Biomedicals, CA, USA) and lysis buffer (ML buffer) from the NucleoSpin kit on dry ice in a FastPrep-24™ 5G instrument (MP Biomedicals, USA) with 3x20sec cycles.
After RNA extraction, the amount of RNA was quantified using a NanoDrop 2000 UV-Vis spectrophotometer (ThermoFisher Scientific, MA, USA). A ratio of sample absorbance at 260 nm and 280 nm in the range of 1.7 to 2.1 was accepted. The RNA which was to be used for microarray analysis was concentrated with a Scan Speed 32 speed vacuum concentrator (Labogene, Denmark). The concentration and quality of the concentrated RNA was determined with a NanoDrop ND1000 spectrophotometer (ThermoFisher Scientific, MA, USA).
Whole genome microarray
Affymetrix whole-transcriptome expression profiling was processed by Swegene centre for integrative biology (SCIBLU) genomics, Affymetrix unit at Lund University, Sweden.
MCAs from 6 stroke females and 6 female controls were analyzed. The integrity of the RNA was measured with the Agilent Bioanalyzer (Agilent Technologies, CA, USA). From a total of 100 ng RNA, single stranded complimentary DNA (cDNA) was synthesized using primers containing a T7 promoter sequence. The single stranded cDNA was converted to double stranded DNA and used as a template for in vitro transcription, producing complimentary RNA (cRNA) . After purification, sense-strand cDNA was synthesized and purified. The sense strand cDNA was fragmented and labeled and loaded onto Affymetrix GeneChip rat gene 2.0 ST arrays (ThermoFisher Scientific, MA, USA). This was followed by hybridization for 16 h at 45 °C in an Affymetrix Gene Chip Hybridization 645 oven. The array was scanned using the Affymetrix GeneChip scanner 3000 7G.
Complimentary DNA synthesis and quantitative PCR
MCAs from 5 female stroke rats, 6 male stroke rats and 6 controls of both sexes were analyzed. After RNA extraction, the amount of RNA was quantified using a NanoDrop 2000 UV-Vis spectrophotometer (ThermoFisher Scientific, MA, USA). cDNA was synthesized using the RT2 First Strand Kit (Qiagen, Germany) according to the manufacturer’s protocol.
The QuantStudio 12 K Flex real-time PCR system (ThermoFisher Scientific, MA, USA) was used for the qPCR. Taqman gene expression assays for Ccl2 (Rn00580555_m1), Olr1 (Rn00591116_m1), Adamts4 (Rn02103282_s1), Serpine1 (Rn01481341_m1), S1pr3 (Rn01757498_m1), Socs3 (Rn01470502_g1), JunB (Rn00572994_s1) and Fosl1 (Rn00564119_m1) were purchased from ThermoFisher Scientific, MA, USA.
The qPCR was performed in a 10-μl reaction volume containing TaqMan 2× universal PCR master mix (ThermoFisher Scientific, MA, USA), 20× TaqMan gene expression assay, RNase-free water and 2 μl cDNA using the QuantStudio 12 K Flex real-time PCR system (ThermoFisher Scientific, MA, USA) with ROX as a passive reference. A no-template control with RNase-free water instead of cDNA was used as negative control for all TaqMan gene expression assays. An inter-plate control for all TaqMan gene expression assays was used to control the thermal cycling between plates. ActB (Rn00667869_m1) and Gapdh (Rn01775763_g1) acted as housekeeping genes. All TaqMan gene expression assays were pipetted in triplicates for each sample.
Analysis and statistics
Analysis of microarray data
Basic Affymetrix GeneChip analysis and experimental quality control were performed using the Expression Console Software (v1.1.2), and the Robust Multi-array analysis method was used for probe summarization and data normalization (quantile normalization and log transformation). Data filtration was done for probe sets having a value less than the median values of the negative control in 80% of total samples.
Significance analysis of microarray was performed using the TMEV software (v4.0). Differentially expressed genes with a false detection rate (q-value) of zero were selected and single probe sets with more than one annotation were excluded. For downstream analyses, we performed enrichment analyses using the PANTHER Classification System [42,43,44]. Enriched PANTHER pathways and protein classes (v15), Reactome pathways (v65) and GO biological process terms (released 2020-02-21) were identified using Fisher’s Exact Test and the p-values were adjusted for multiple testing using the conservative Bonferroni correction. P < 0.05 was considered statistically significant.
For the biological process GO terms, only the most specialized term within each hierarchical group was included and GO terms with less than 30 annotated genes from the list of differentially expressed genes were excluded. For the PANTHER protein classes and Reactome pathways, only the least specialized protein class/pathway within each hierarchical group was included (cut-off: at least 10 annotated genes). On the contrary, only the most specialized pathway within each hierarchical group was included for the PANTHER pathways (cut-off: at least 5 annotated genes). Subsequently, overlap between overrepresented GO terms for the differentially expressed genes in the occluded MCAs and non-occluded MCAs both compared with control MCAs were identified.
To explore interactions between selected gene products, we utilized the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database v11.0 . Only experimentally determined and database curated Rattus norvegicus protein-protein interactions were used to create the networks. We defined a cluster as a network formed by at least 5 interacting proteins.
Gene selection and GO term categorization
Based on the results from the microarray, eight genes were selected for validation with qPCR. Ccl2, S1pr3, Socs3, Serpine1, JunB and Fosl1 were selected due to occurrence in biological processes that were significantly enriched in the GO analysis in combination with analysis of fold change and review of literature. Adamts4 was selected due to a very high relative expression, and Olr1 was selected due to the availability of previous data in the same model in hypertensive and normotensive males. In addition to qPCR, we categorized the selected genes into preselected GO biological process terms using the PANTHER Classification System [42,43,44] and the statistical software R (v4.0.2) to examine their involvement in specific biological processes.
In the PCR-validation, the mean computed cycle threshold (CT) for the triplicates of each sample was calculated. The results were normalized against the inter-plate control. The mean CT-value from the housekeeping genes ActB and Gapdh were subtracted from the mean CT value of the target gene to be able to compare samples that contained a differing amount of total-RNA. This relative value is hereafter referred to as delta-CT (dCT). For each gene, dCT values more than 3 standard deviations away from the mean were considered outliers and excluded (only one observation was considered an outlier across the datasets). Subsequently, the effects of sex (males and females) and experimental group (control MCAs, non-occluded MCAs and occluded MCAs) on dCT were investigated by means of linear mixed-effect models with restricted maximum likelihood using the statistical software R (v4.0.2) and the nlme R package (v3.1.148) . Sex and experimental group were fitted as fixed effects (fixed effects were omitted from the model if not reaching significance) and animal ID as a random effect to consider individual animal differences. To assess normality and equality of variance, we created and examined Q-Q plots and residuals versus fitted values plots, respectively. These data are presented using boxplots and statistical significance was set at p < 0.05. Since a higher dCT value corresponds with a lower expression, the axes in the figures were reversed to offer a more intuitive reading of the figures.
Cross-analysis with findings from Grell et al. 
To reveal sex similarities between males and females at a larger scale, we compared the findings from the current study with findings obtained by Grell et al. . Here, they investigated differential expression of genes in the occluded MCAs compared with the non-occluded MCAs in WKY male rats. They used the same stroke model (tMCAO), the rats were the same age (12-weeks) and the samples were processed by the same center using the same method (SCIBLU). Even though the groups are similar in many ways, some of the variations observed might be due to strain differences rather than sex differences. Therefore, we decided to focus solely on the overlapping differentially expressed genes.
Availability of data and materials
The datasets generated or analysed during the current study are available in the Gene Expression Omnibus (GEO) repository, GSE162072, and in the published article [and its supplementary information file ‘Rehnstrom_et_al_supplementary_files.xlsx’].
- ActB :
- Adam :
ADAM metallopeptidase domains
- Adam8 :
ADAM metallopeptidase domain 8
- Adamts :
A disintegrin and metalloproteinase with thrombospondin motifs
- Adamts4 :
A disintegrin and metalloproteinase with thrombospondin type 1 motif, 4
Activator protein 1
Activating transcription factor–2
- Ccl :
Chemokine (C-C motif) ligand
- Ccl2 :
Chemokine (C-C motif) ligand 2
- Chi3l1 :
Chitinase-3-like protein 1
- Cxcl :
Chemokine (C-X-C motif) ligand
Delta cycle threshold
ETS transcription factor
- Fkbp3 :
FK506 binding protein 3
- Fosl1 :
Fos-like antigen 1
- Gapdh :
Glyceraldehyde 3-phosphate dehydrogenase
- Icam1 :
Intercellular adhesion molecule 1
- Il :
- Il-6 :
Jun amino-terminal kinases
- JunB :
Jun B proto-oncogene
- Nfkb1 :
Nuclear factor of kappa light polypeptide gene enhancer in B-cells 1
- Rela :
Nuclear factor kappaB subunit p65
Middle cerebral artery
- Mmp :
- Mmp19 :
Matrix metallopeptidase 19
- Olr1 :
Oxidized low-density lipoprotein receptor 1
Protein ANalysis THrough Evolutionary Relationships
Phosphorylated extracellular signal-regulated kinase 1 and 2
Phosphorylated c-jun N-terminal kinase
Phosphorylated protein kinase C
Phosphorylated protein kinase C-δ
- Prdm4 :
PR domain zinc finger protein 4
Quantitative real-time polymerase chain reaction
Regional cerebral blood flow
Reactive oxygen species
- Runx1 :
Runt-related transcription factor 1
- S1pr3 :
Sphingosine 1 phosphate receptor 3
Stress-activated protein kinases
- Serpine1 :
Serine protease inhibitor, clade E, member 1
- Socs3 :
Suppressor of cytokine signaling 3
- Stat3 :
Signal transducer and activator of transcription 3
Search Tool for the Retrieval of Interacting Genes/Proteins
- Ttc39b :
Tetratricopeptide repeat domain 39B
Transient middle cerebral artery occlusion
- TNFα :
Tumor necrosis factor α
Vascular smooth muscle cell
- Zc3h12a :
Zinc finger CCCH type containing 12A
Benjamin EJ, Virani SS, Callaway CW, Chamberlain AM, Chang AR, Cheng S, et al. Heart disease and stroke Statistics-2018 update: a report from the American Heart Association. Circulation. 2018;137:e67–492.
Appelros P, Stegmayr B, Terént A. Sex differences in stroke epidemiology: a systematic review. Stroke. 2009;40:1082–90.
Ahnstedt H, McCullough LD, Cipolla MJ. The importance of considering sex differences in translational stroke research. Transl Stroke Res. 2016;7:261–73. https://doi.org/10.1007/s12975-016-0450-1.
Kent DM, Price LL, Ringleb P, Hill MD, Selker HP. Sex-based differences in response to recombinant tissue plasminogen activator in acute ischemic stroke: a pooled analysis of randomized clinical trials. Stroke. 2005;36:62–5.
Sandoval KE, Witt KA. Blood-brain barrier tight junction permeability and ischemic stroke. Neurobiol Dis. 2008;32:200–19.
Wardlaw JM, Murray V, Berge E, del Zoppo GJ. Thrombolysis for acute ischaemic stroke. Cochrane Database Syst Rev. 2014;2014:CD000213.
O’Collins VE, Macleod MR, Donnan GA, Horky LL, van der Worp BH, Howells DW. 1,026 experimental treatments in acute stroke. Ann Neurol. 2006;59:467–77.
Beery AK, Zucker I. Sex bias in neuroscience and biomedical research. Neurosci Biobehav Rev. 2011;35:565–72.
Simpkins JW, Rajakumar G, Zhang YQ, Simpkins CE, Greenwald D, Yu CJ, et al. Estrogens may reduce mortality and ischemic damage caused by middle cerebral artery occlusion in the female rat. J Neurosurg. 1997;87:724–30.
Zhang L, Li P-P, Feng X, Barker JL, Smith SV, Rubinow DR. Sex-related differences in neuronal cell survival and signaling in rats. Neurosci Lett. 2003;337:65–8.
Du L, Bayir H, Lai Y, Zhang X, Kochanek PM, Watkins SC, et al. Innate gender-based proclivity in response to cytotoxicity and programmed cell death pathway. J Biol Chem. 2004;279:38563–70.
Ahnstedt H, Cao L, Krause DN, Warfvinge K, Säveland H, Nilsson OG, et al. Male-female differences in upregulation of vasoconstrictor responses in human cerebral arteries. PLoS One. 2013;8:e62698.
Ahnstedt H, Mostajeran M, Blixt FW, Warfvinge K, Ansar S, Krause DN, et al. U0126 attenuates cerebral vasoconstriction and improves long-term neurologic outcome after stroke in female rats. J Cereb blood flow Metab Off J Int Soc Cereb Blood Flow Metab. 2015;35:454–60.
Edvinsson LIH, Povlsen GK. Vascular plasticity in cerebrovascular disorders. J Cereb blood flow Metab Off J Int Soc Cereb Blood Flow Metab. 2011;31:1554–71.
Yin K-J, Deng Z, Hamblin M, Zhang J, Chen YE. Vascular PPARδ protects against stroke-induced brain injury. Arterioscler Thromb Vasc Biol. 2011;31:574–81.
Vikman P, Ansar S, Edvinsson L. Transcriptional regulation of inflammatory and extracellular matrix-regulating genes in cerebral arteries following experimental subarachnoid hemorrhage in rats. Laboratory investigation J Neurosurg. 2007;107:1015–22.
Vikman P, Ansar S, Henriksson M, Stenman E, Edvinsson L. Cerebral ischemia induces transcription of inflammatory and extracellular-matrix-related genes in rat cerebral arteries. Exp Brain Res. 2007;183:499–510.
Vikman P, Edvinsson L. Gene expression profiling in the human middle cerebral artery after cerebral ischemia. Eur J Neurol. 2006;13:1324–32.
Ansar S, Edvinsson L. Subtype activation and interaction of protein kinase C and mitogen-activated protein kinase controlling receptor expression in cerebral arteries and microvessels after subarachnoid hemorrhage. Stroke. 2008;39:185–90.
Grell A-S, Mostajeran M, Frederiksen SD, Edvinsson L, Ansar S. Cerebrovascular gene expression in spontaneously hypertensive rats after transient middle cerebral artery occlusion. Neuroscience. 2017;367:219–32.
Mostajeran M, Wetterling F, Blixt FW, Edvinsson L, Ansar S. Acute mitogen-activated protein kinase 1/2 inhibition improves functional recovery and vascular changes after ischaemic stroke in rat-monitored by 9.4 T magnetic resonance imaging. Acta Physiol (Oxf). 2018;223:e12985.
Mostajeran M, Edvinsson L, Warfvinge K, Singh R, Ansar S. Inhibition of mitogen-activated protein kinase 1/2 in the acute phase of stroke improves long-term neurological outcome and promotes recovery processes in rats. Acta Physiol (Oxf). 2017;219:814–24.
Parker BL, Larsen MR, Edvinsson LIH, Povlsen GK. Signal transduction in cerebral arteries after subarachnoid hemorrhage-a phosphoproteomic approach. J Cereb blood flow Metab Off J Int Soc Cereb Blood Flow Metab. 2013;33:1259–69.
Wågsäter D, Björk H, Zhu C, Björkegren J, Valen G, Hamsten A, et al. ADAMTS-4 and -8 are inflammatory regulated enzymes expressed in macrophage-rich areas of human atherosclerotic plaques. Atherosclerosis. 2008;196:514–22.
Hsu Y-P, Staton CA, Cross N, Buttle DJ. Anti-angiogenic properties of ADAMTS-4 in vitro. Int J Exp Pathol. 2012;93:70–7.
Hamel MG, Ajmo JM, Leonardo CC, Zuo F, Sandy JD, Gottschall PE. Multimodal signaling by the ADAMTSs (a disintegrin and metalloproteinase with thrombospondin motifs) promotes neurite extension. Exp Neurol. 2008;210:428–40.
Ren P, Zhang L, Xu G, Palmero LC, Albini PT, Coselli JS, et al. ADAMTS-1 and ADAMTS-4 levels are elevated in thoracic aortic aneurysms and dissections. Ann Thorac Surg. 2013;95:570–7.
Dong H, Du T, Premaratne S, Zhao CX, Tian Q, Li Y, et al. Relationship between ADAMTS4 and carotid atherosclerotic plaque vulnerability in humans. J Vasc Surg. 2018;67:1120–6.
Kenagy RD, Min S-K, Clowes AW, Sandy JD. Cell death-associated ADAMTS4 and versican degradation in vascular tissue. J Histochem Cytochem Off J Histochem Soc. 2009;57:889–97.
Campbell RA, Vieira-de-Abreu A, Rowley JW, Franks ZG, Manne BK, Rondina MT, et al. Clots are potent triggers of inflammatory cell gene expression: indications for timely fibrinolysis. Arterioscler Thromb Vasc Biol. 2017;37:1819–27.
Qi J, Kreutzer DL. Fibrin activation of vascular endothelial cells. Induction of IL-8 expression. J Immunol. 1995;155:867–76.
Hughes PM, Allegrini PR, Rudin M, Perry VH, Mir AK, Wiessner C. Monocyte chemoattractant protein-1 deficiency is protective in a murine stroke model. J Cereb blood flow Metab Off J Int Soc Cereb Blood Flow Metab. 2002;22:308–17.
Theilmeier G, Schmidt C, Herrmann J, Keul P, Schäfers M, Herrgott I, et al. High-density lipoproteins and their constituent, sphingosine-1-phosphate, directly protect the heart against ischemia/reperfusion injury in vivo via the S1P3 lysophospholipid receptor. Circulation. 2006;114:1403–9.
Salomone S, Yoshimura S, Reuter U, Foley M, Thomas SS, Moskowitz MA, et al. S1P3 receptors mediate the potent constriction of cerebral arteries by sphingosine-1-phosphate. Eur J Pharmacol. 2003;469:125–34.
Singleton PA, Moreno-Vinasco L, Sammani S, Wanderling SL, Moss J, Garcia JGN. Attenuation of vascular permeability by methylnaltrexone: role of mOP-R and S1P3 transactivation. Am J Respir Cell Mol Biol. 2007;37:222–31.
Demicheva E, Hecker M, Korff T. Stretch-induced activation of the transcription factor activator protein-1 controls monocyte chemoattractant protein-1 expression during arteriogenesis. Circ Res. 2008;103:477–84.
Goldman JM, Murr AS, Cooper RL. The rodent estrous cycle: characterization of vaginal cytology and its utility in toxicological studies. Birth Defects Res B Dev Reprod Toxicol. 2007;80:84–97.
Memezawa H, Minamisawa H, Smith ML, Siesjö BK. Ischemic penumbra in a model of reversible middle cerebral artery occlusion in the rat. Exp Brain Res. 1992;89:67–78.
Bederson JB, Pitts LH, Tsuji M, Nishimura MC, Davis RL, Bartkowski H. Rat middle cerebral artery occlusion: evaluation of the model and development of a neurologic examination. Stroke. 1986;17:472–6.
Engelhorn T, Goerike S, Doerfler A, Okorn C, Forsting M, Heusch G, et al. The angiotensin II type 1-receptor blocker candesartan increases cerebral blood flow, reduces infarct size, and improves neurologic outcome after transient cerebral ischemia in rats. J Cereb blood flow Metab Off J Int Soc Cereb Blood Flow Metab. 2004;24:467–74.
Van Gelder RN, von Zastrow ME, Yool A, Dement WC, Barchas JD, Eberwine JH. Amplified RNA synthesized from limited quantities of heterogeneous cDNA. Proc Natl Acad Sci U S A. 1990;87:1663–7.
Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13:2129–41.
Mi H, Muruganujan A, Ebert D, Huang X, Thomas PD. PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019;47:D419–26.
Mi H, Muruganujan A, Huang X, Ebert D, Mills C, Guo X, et al. Protocol update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat Protoc. 2019;14:703–21.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. _nlme: Linear and Nonlinear Mixed Effects Models_. R package version 3.1–148. 2020. https://cran.r-project.org/package=nlme.
We would like to acknowledge SCIBLU genomics, the Affymetrix unit at Lund University in Sweden for processing of the microarray technology. Anne-Sofie Jagd Grell is thanked for training MR in wet lab techniques.
The study was supported by the Knut and Alice Wallenberg Foundation (awards grants to basic research), the Swedish Heart-Lung Foundation (awards grants to vascular research) and the Swedish Research Council (grant no 5958, support all scientific fields), as well as by a grant from Lund University Hospital, Sweden. The funding bodies provided salary for some of the researchers and/or monetary support to cover the animal costs (e.g. housing, surgery, anesthesia products and gases), and the microarray and qPCR analysis. The funding bodies were not involved in the design, analysis and reporting of the study. Open Access funding provided by Lund University.
Ethics approval and consent to participate
All experiments and procedures were performed in accordance with guidelines for the European Community Council Directive (2010/63/EU) for Protection of Vertebrate Animals Used for Experimental and other Scientific Purposes. The experimental procedures were approved by Lund-Malmö Institutional Ethics Committee under the auspices of the Swedish National Department of Agriculture (ethical permit number M153–15). The study complies with the ARRIVE guidelines (Animal Research: Reporting in Vivo Experiments).
Consent for publication
LE is the CSO and founder of the Edvince research company which aims at developing and commercializing new methods for acute treatment of stroke. However, this did not affect the outcomes of the study. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Rehnström, M., Frederiksen, S.D., Ansar, S. et al. Transcriptome profiling revealed early vascular smooth muscle cell gene activation following focal ischemic stroke in female rats – comparisons with males. BMC Genomics 21, 883 (2020). https://doi.org/10.1186/s12864-020-07295-2
- Focal cerebral ischemia
- Gene regulation
- Female rats
- Sex differences
- Transcription factors
- Endothelial function
- Pathway analysis