Bioscience Reports
Portland Press Ltd.
Effect of X-rays on transcript expression of rat brain microvascular endothelial cells: role of calcium signaling in X-ray-induced endothelium damage
Volume: 40, Issue: 4
DOI 10.1042/BSR20193760

Radiation-induced brain edema is a serious adverse effect of radiotherapy. Although there are many causes of radiation-induced brain edema, the pathogenesis is not clear and clinical treatment is not ideal. Therefore, knowing the differential expression of the brain microvascular endothelial cell (BMEC) transcriptome after brain radiotherapy may shed light on the pathogenesis of radiation-induced brain edema. The present study used RNA-Seq technique to identify 383 BMEC transcripts differentially expressed (many 2-fold or higher; P < 0.05) between control and X-ray–treated primary cultured rat BMECs. Compared with controls, X-ray–treated BMECs had 183 significantly up-regulated transcripts and 200 significantly down-regulated transcripts. The differentially expressed genes were associated with the biological processes of the cell cycle, apoptosis, vascular permeability, and extracellular junctions. The functional changes identified in the X-ray–treated BMECs included Ca2+ signaling, phosphoinositide 3-kinase–Akt signaling, and methionine degradation. These results indicated that transcript expression was substantially affected by radiation exposure and the proteins encoded by these differentially expressed genes may play a significant role in radiotherapy-induced brain edema. Our findings provide additional insight into the molecular mechanisms of radiation-induced brain edema and may be helpful in the development of clinical treatment of this adverse reaction to radiotherapy.

Wu, Fang, Zhang, Song, Wang, Chen, Du, Zheng, and Shen: Effect of X-rays on transcript expression of rat brain microvascular endothelial cells: role of calcium signaling in X-ray-induced endothelium damage


Radiation therapy is one of the most common treatments for malignant tumors of the head and neck. It greatly benefits patients that have nasopharyngeal carcinoma, primary brain tumor, or brain metastases [1], in particular, radiation therapy is currently the preferred treatment for patients with nasopharyngeal carcinoma [2]. However, recent clinical studies have found that brain microvascular endothelial cell (BMEC) damage and radiation-induced brain edema are common in patients with intracranial tumors after radiotherapy [3]. Radiation-induced cerebral edema may cause significant decrease in the utilization of glucose by brain tissue, leading to ischemia and hypoxia, which are the main causes of the exacerbation of clinical symptoms and neurological deficits [4,5]. The current treatment of radiation-induced brain injury in clinical practice is challenging. Despite many studies investigating radiation-induced brain injury, its pathogenesis is still unclear [6].

At present, it is believed that BMECs are highly sensitive to radiation and may be the most important effector cells of radiation-induced cerebral edema [7,8]. The apoptosis or death of BMECs can directly cause capillary structure abnormalities and blood–brain barrier (BBB) disintegration [9]. Intercellular junction types include tight junctions, desmosomes, adherens junctions, and gap junctions [10]. Tight junctions consist primarily of intrinsic transmembrane proteins, scaffold proteins, and regulatory molecules between cells [10]. Studies have shown that tight junctions between vascular endothelial cells play a critical role in the BBB [11]. Thus, most studies on the effects of radiation on blood vessels have focused on endothelial cells [7].

High-throughput technology enables comprehensive monitoring of biological systems and provides an entry point to clarify molecular mechanisms that are involved in cell responses to environment changes. Therefore, we investigated the effect of a dose of X-rays (20 Gy) on transcripts in rat BMECs [11,12] to provide a theoretical basis for understanding the molecular mechanisms underpinning radiotherapy-induced brain edema. The results of the present study may also offer potential targets for future drug development.

Materials and methods

Primary BMEC culture

Sprague Dawley rats (2–3 weeks old, 40–60 g body weight, male and female) were purchased from the Animal Center of Anhui Medical University, and were housed at 23 ± 2°C in separate cages with water, and fed ad libitum in a 12-h reverse light cycle. All animal experiments were performed in the department of physiology, and procedures used in the present study were conducted in accordance with the National Institutes of Health publication No. 8523 and were approved by the Animal Experimentation Ethics Committee of Anhui Medical University (No. LLSC20150048). Rats were anesthetized in a chamber containing 2% isoflurane mixed with 0.2 l/min 100% O2, and killed by an overdose of carbon dioxide (CO2), and then placed in 75% ethanol for 2 min to sterilize the body. The whole brain was removed and placed in a Petri dish containing cold phosphate-buffered saline. The cerebral hemispheres were then slowly rolled over dry filter paper to remove the pia mater and meningeal great vessels. The white matter, residual great vessels, and pia meninges were removed with fine dissection forceps, and the cerebral cortex was preserved. Next, the tissue was placed in 1 ml of Dulbecco’s modified Eagle’s medium (DMEM), and cut with ophthalmic scissors into 1-mm3 specimens. The tissue was then digested with 0.1% type II collagenase (containing 30 U/ml DNase I) for 1.5 h in a water bath at 37°C. After centrifugation at 1000 × g for 8 min at room temperature, the supernatant was removed. Bovine serum albumin (20%) was added, and the tissue was centrifuged again (1000 × g, 20 min, 4°C) to remove the upper nervous tissue and large vessels. Collagenase/dispase (2 ml 0.1%, containing 20 U/ml DNase I) was added to the remaining precipitate, which was allowed to digestfor 1 h at 37°C in a water bath. Following centrifugation (1000 × g, 8 min, room temperature), the tissue was suspended in 2 ml of DMEM and separated in a continuous gradient of 12 ml of 50% Percoll (25,000 g/l, 60 min, 4°C) with centrifugation The white and yellow layer observed near the bottom of the red blood cell layer contained the purified microvessels. After a wash, DMEM containing 20% fetal bovine serum (FBS) and 100 μg/ml heparin sodium was added to the purified microvessels, and the cells were inoculated onto a 35-mm Petri dish coated with rat tail collagen and placed in an incubator with 5% CO2 set at 37°C for 12–24 h. The next change of cell culture medium used fibroblast growth factor-basic (1 ng/ml) added to Medium M199.

X-ray treatment

Rat BMECs were cultured in Medium M199 containing 10% FBS, 100 U/ml penicillin, and 100 μg/ml streptomycin in a humidified atmosphere of 5% CO2 at 37°C. When cell confluence reached 60%, the BMECs were irradiated with X-rays at room temperature for 8 min (2.5 Gy/min; total 20 Gy). After being cultured for an additional 24 h, the cells were used in following experiments [12,13].

RNA extraction

The collected samples were incubated with lysis reagent for 5 min at room temperature. For every milliliter of lysis reagent, 0.2 ml of chloroform was added, and the samples were mixed using an up and down motion for 15 s and then incubated for 2–3 min at room temperature. The samples were centrifuged at 12,000 × g for 10 min at 4°C. The separated aqueous phase was transferred to a fresh tube, and 0.5 ml of isopropanol was added for a 10-min incubation at room temperature. The sample was then centrifuged at 12,000 × g for 10 min at 4°C. After being washed with 75% ethanol and centrifuged at 7500 × g for 5 min at 4°C, the supernatant was removed and air-dried at room temperature. The remaining precipitate was resuspended in 15–30 μl of diethylpyrocarbonate (DEPC)-treated water. The RNA suspension was mixed with a 1–10 part by volume of 3 M sodium acetate (pH 5.5) and 1 part by volume of isopropanol. After being incubated for 20 min at room temperature and centrifuged at 12,500 × g for 10 min at 4°C, the pellet was washed with ice-cold 70% ethanol in DEPC-treated water and centrifuged again at 10,000 × g for 5 min at 4°C. The supernatant was removed and dried at room temperature. Finally, the RNA was resuspended in DEPC-treated water. The concentration (ng/μl) and purity (absorption ratio of 260 to 230 nm) of the extracted RNA were determined using a spectrophotometer to ensure a 260/230 absorption ratio of 2.0. The RNA samples were then sequenced.

Next-generation sequencing and data analysis

The RNA samples were sequenced using the high-throughput Illumina HiSeq 2500 sequencing platform. The sequencing length was 150 base pairs. The RNA-Seq FastQ raw data were pruned using Trimmomatic to remove adaptors and lower mass readings [14]. The quality of clean data was assessed using FastQC software [15]. The quality-approved data were mapped to the reference genome of the rat (National Center for Biotechnology Information [NCBI] genome assembly version Rnor_6.0) using HISAT2 (v2.0.13) and annotated with the annotation file (.gtf) NCBI Rnor_6.0 [16,17]. The transcript expression level was then calculated using the most common method, that is, Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) [18]. The results were used to compare the differences in transcript expression among the samples. The RNA-Seq quantification software Kallisto was used to obtain a count of known mRNAs [19,20]. The differences in the number of reads were analyzed using edgR package software ( [21]. The gene significance levels between the two groups were calculated using a negative binomial model. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on differentially expressed transcriptomes, with a 2-fold change cutoff [22].

GO and KEGG pathway enrichment analyses

The functional pathway enrichment of the proteins encoded by candidate genes were analyzed, and these genes were annotated using the Metascape database [9]. GO annotations were performed using the Metascape online tool on the screened differentially expressed transcripts. The KEGG pathway analysis of differentially expressed genes (DEGs) was also performed using the Metascape online analysis database. The DEGs that were significantly up- or down-regulated as determined using RNA-Seq data were further analyzed, and a two-sided P value < 0.05 was considered statistically significant.

Protein–protein interaction (PPI) network integration

PPI analysis can provide new insights into protein function. It may also help to reveal the general organization principles of functional cellular networks. The online search tool STRING ( was used to find interacting genes/proteins and to reveal functional relationships between proteins [23]. We constructed PPI networks of high-expression genes and of low-expression genes to predict interactions between the selected genes. Cytoscape ( is widely used to make biomolecular interaction networks into models to construct PPI networks of differentially expressed mRNAs [24]. Most previously acquired biological networks have been found to be subject to scale-free attribution [25]. Therefore, we used the molecular complex detection (MCODE) in Cytoscape software (version 3.6.0) to screen modules with MCODE scores ≥ 3 and nodes ≥ 3 in the PPI network [26]. We used topology analysis to analyze the connectivity of nodes in the PPI network so as to obtain a higher degree of important nodes (central proteins) [27]. The top 10 hub genes were selected for association. The functional enrichment analysis of individual modules was performed using Metascape, with a threshold considered statistically significant of P<0.05.

Quantitative real-time polymerase chain reaction (qPCR)

We used qPCR experiment to test mRNA expression change. Briefly, cDNA was synthesized by adding total RNA into a kit of HiScript III RT SuperMix for qPCR (Vazyme, China). We mixed total RNA (4 μl), 4 × gDNA wiper Mix (4 μl) and RNase-free ddH2O (8 μl) together. The mixture was heated to 42°C for 2 min in a water bath. Then we added 4 μl HiScript III qRT SuperMix containing buffer, dNTP, HiScript III reverse transcriptase, RNase inhibitor and Random primers/Oligo (dT) 20VN primer mix into the heated mixture. Reaction was heated to 50°C for 15 min and then 85°C for 5 s. The Roche LightCycler 480 II PCR equipment (Roche, Switzerland) was used for qPCR reaction. Thermocycler steps were 5 min at 95°C, followed by 40 cycles of 15 s at 95°C, 30 s at 60°C. The primer sequences are listed in Table 1. The 2−ΔΔCT method was used to analyze relative gene expression. Gapdh was used as a control.

Table 1
The primer sequences used in qPCR
GeneForward primerReverse primer


All data were expressed as mean ± SEM. Two-tailed Mann–Whitney U-test was used to determine statistical significance. Differences were considered significant with a value of P<0.05.


Quality of RNA sequencing

We sequenced and obtained raw data from six primary cultured rat BMEC samples through the Illumina platform. After trimming and filtering, the quality of data was good. The median score of the reads in all samples was above 30, and all samples passed the quality check (Figure 1A).

Sequencing quality
Figure 1
(A) Box plot showing the sequencing quality of the six sample reads. The X-axis is the base position for each reading (bp indicates base pair). The Y-axis is the sequence quality score. A median score of less than 20 in the sequencing quality score means poor sequencing quality, and above 30 means good sequencing quality. (B) Mapping ratio of each sample in control (C1-C3) and X-ray–treated (S1-S3) groups. Gray bars represent samples that were sequenced. The abscissa represents the sample number and the ordinate represents the percentage of the alignment.Sequencing quality

The trimmed high-quality reads were used for further data analysis. Mapping and quantification were conducted using the standard pipelines of HISAT2 and StringTie, which are based on the rat reference genome annotation file of Ensembl and its release 92 of Rnor_6.0 [17]. On average, the mapping ratios for all six samples were high (>90%), indicating that the reference genome was appropriate and non-contaminated (Figure 1B) [28]. We examined the rationality and experimental reliability of the sample selection by analyzing the correlation of mRNA expression levels using correlation coefficient maps among the sequenced samples (Figure 2A) and using principal component analysis (Figure 2B). The results showed that the correlations among the groups were relatively high and that the distribution of the samples was consistent with the treated groups. There was a significant difference between the treated group and the control group. The results indicated that the experimental data were of both high quality and high reliability.

Transcript expression correlation and principal component analysis
Figure 2
(A) Transcript-based sample correlation analysis showing the transcript expression level correlation between control (C1-C3) and X-ray–treated (S1-S3) groups. The abscissa indicates the sample name, the ordinate indicates the corresponding sample name, and the color indicates the correlation coefficient size. (B) Principal component analysis based on log2 (FPKM+1) in control (C) and X-ray–treated (S) groups. Each point in the coordinate system represents a sample, with red dots representing control (C1-C3) groups and blue triangles representing X-ray–treated (S1-S3) groups. PC1 is the main component 1, and PC2 is the main component 2. The percentage refers to the contribution of this principal component to the overall variance. The distance between the points reflects the difference between the samples they correspond to.Transcript expression correlation and principal component analysis

Identification of DEGs

We investigated the overall transcript expression levels of BMECs in the control and X-ray–treated groups by using log10 (FPKM+1) values and found that the overall changes were not significantly different from the control groups (Figure 3A). Cluster analysis of the differentially expressed mRNAs showed significant differences among the groups (Figure 3B). A total of 383 genes were significantly differentially expressed, of which 183 were up- and 200 were down-regulated. The results of the differentially expressed mRNAs are displayed as a volcano map in Figure 4A. In primary cultured rat BMECs irradiated with X-rays, the top 10 up-regulated DEGs were Ass1, Prkar2b, Agtpbp1, Dnajc27, Ap3b1, Tbc1d31, Slc1a3, Osbpl5, Acsbg1 and Yif1b. The top 10 down-regulated DEGs were Cttn, Ppp1r18, Znf106, Ltbp2, Anln, Adgrb2, Prom1, Aggf1, Klhl15 and LOC108349244 (Figure 4B).

Overview of expression profiles of transcripts
Figure 3
(A) Box plot showing an overview of the expression profiles of transcripts in control (C1-C3) and X-ray–treated (S1-S3) groups. X-axis shows the sample name; Y-axis is the log2 (FPKM+1) value. (B) Transcript-based hierarchical clustering analysis combining tissues with similar natures. Scaled to log2 (FPKM+1) expression values by different color intensities, the Z scores of each transcript are normalized, and green and red represent low and high expression levels, respectively.Overview of expression profiles of transcripts
Differential expression of data between two sets of samples
Figure 4
(A) Volcano map. Red points represent upregulated transcripts screened on the basis of absolute fold change ≥ 2.0 and a corrected P value of < 0.05. Green points represent the expression of transcripts that were down-regulated, screened on the basis of absolute fold change ≥ 2.0 and a corrected P value of < 0.05. The black points represent transcripts with no significant change. (B) The most substantially (P<0.05) up-regulated and down-regulated transcripts based on fold changes following irradiation of endothelial cells (n=3).Differential expression of data between two sets of samples

DEGs reflect changed biological functions

The analysis of DEGs by GO functional annotation for the category biological process revealed that the up-regulated transcripts were mainly involved in regulation of lymphocyte activation, positive regulation of apoptotic cell clearance, positive regulation of cell migration, and response to nutrient levels. The down-regulated transcripts were primarily involved in positive regulation of cell cycle, reproductive structure development, microtubule-based process and learning or memory. Both of the up- and down-regulated transcripts were involved in response to glucocotticoid, regulation of developmental growth and cellular response to interferon-gamma (Figure 5A).

GO functional category enrichment following irradiation of endothelial cells
Figure 5
Heat maps of the top 20 enriched (A) biological processes, (B) cell composition, and (C) molecular functions following irradiation of endothelial cells. The result on the left is from the down-regulated gene list, and the right is from the up-regulated gene list. The ordinate represents the corresponding Go entry. Colors represent P values (-log10 (P)).GO functional category enrichment following irradiation of endothelial cells

For the cellular component category, the up-regulated transcripts were mainly involved in the composition of cells in some components of the biofilm and plasma membrane protein complexes. The downregulated transcripts were mainly involved in some cell cycle–related components, such as chromosomes, telomeres, and central granules (Figure 5B).

For the category of molecular function, the up-regulated transcripts were mainly involved in receptor modulator activity, chemokine activity, drug transmembrane transporter activity, protein tyrosine kinase activator activity, cationic transmembrane transporter activity, and cytokine receptor activity. The down-regulated transcripts were primarily involved in cell adhesion molecule binding, histone acetyltransferase binding and protein domain-specific junctions, catalytic activity acting on DNA, nucleotide transferase activity, nuclear hormone receptor binding, chromatin binding, and unfolded protein binding (Figure 5C).

We found that the largest number of DEGs were related to the metabolic pathway. Interestingly, several important pathways, including the Ca2+ signaling pathway, phosphoinositide 3-kinase (PI3K)–Akt signaling pathway, and methionine degradation, were identified to be in the top 20 enriched pathway terms (Figure 6). Because Ca2+ plays an important role in endothelial cell function and vascular diseases, we examined the up-regulated transcripts that were enriched in the Ca2+ signaling pathway, and the results are shown in Figure 7.

KEGG enrichment results
Figure 6
Scatterplot of the top 20 KEGG enrichment results of differentially expressed transcripts in each pairwise comparison that were annotated in the particular pathway term. The X-axis indicates the rich factor, and the Y-axis indicates pathway. RichFactor is the ratio of differentially expressed gene numbers to all gene numbers annoted in the pathway term. A high RichFactor represents greater intensiveness. A lower P value represents greater intensity.KEGG enrichment results
Ca2+ signaling pathway enrichment based on Kyoto Encyclopedia of Genes and Genomes results
Figure 7
Red color represents a gene that is both up-regulated and involved in Ca2+ signaling pathway.Ca2+ signaling pathway enrichment based on Kyoto Encyclopedia of Genes and Genomes results

PPI network integration

We constructed a PPI network encoded by the DEGs (Figure 8A). After we analyzed the PPI network of significant differentially expressed transcripts, four modules were obtained. These genes participate in the main pathways of positively regulated cell cycle processes: DNA damage response, signal transduction by p53 class mediator resulting in cell cycle arrest, positive regulation of cell adhesion, mitotic prometaphase, positive regulation of leukocyte chemotaxis, and chromatin organization (Figure 8B) [29]. Furthermore, as shown in Figure 9, the 10 most notable genes showing significant interactions were Cdkn1a, Ccl2, Ngf, Prkar2b, Cdc42, Hif1a, Rhog, Bub3, Vegfa, and Plk4.

Protein–protein interaction
Figure 8
(A) Interaction network and (B) modules. Circles represent genes, lines represent protein interactions between genes, and line colors represent evidence of interactions between proteins.Protein–protein interaction
Top 10 hub genes
Figure 9
Edges represent the protein–protein associations. Red color indicates genes with high scores. Yellow and orange color indicate genes with middle scores, increasing respectively.Top 10 hub genes

Examination of top gene expression by qPCR

To confirm the change of mRNA expression analyzed by RNA-Seq, we used qPCR to examine the mRNA expression levels of top five genes, including top five up-regulated genes (Ass1, Prkar2b, Agtpbp1, Dnajc27, Ap3b1) and top five down-regulated genes (Cttn, Ppp1r18, Znf106, Ltbp2, Anln). Our data showed that the expression levels of top five up-regulated genes (Ass1, Prkar2b, Agtpbp1, Dnajc27, Ap3b1) were significantly increased (Figure 10A), but top five down-regulated genes (Cttn, Ppp1r18, Znf106, Ltbp2, Anln) were significantly decreased (Figure 10B) in X-rays treatment group compared with control group. Our result of qPCR suggested that the findings from RNA-Seq were trustable.

Expression changes of top five differential expression genes resolved by RNA sequence
Figure 10
The mRNA relative expression levels of top five up-regulated (A) and down-regulated (B) genes in control (Ctrl) and X-ray–treated (X-ray, 20 Gy) brain microvascular endothelial cell of rats. Values are shown as the mean ± SEM (n=3–6). *P<0.05 for X-ray vs. Ctrl group.Expression changes of top five differential expression genes resolved by RNA sequence


The cause of radiation-induced brain edema varies and the related pathogenesis is complex, which makes it difficult to treat in the clinic [8]. Two main cell types are sensitive to radiation damage: oligodendrocytes and vascular endothelial cells [7]. Injury to these cell types will directly lead to the demyelination of nerve fibers and the enhancement of capillary permeability and will eventually trigger a partial collapse of the BBB system. Therefore, it is important to understand the differential expression of the genes related to the BMEC transcriptome after brain radiotherapy and the effects of this differential expression on radiation-induced brain edema.

In the present study, we found that the DEGs were associated with endothelial cell function, and the expression changes of DEGs could be confirmed by qPCR experiments. These findings may be the first evidence showing a linkage between changes in BMEC transcriptome and function after radiation exposure. Among the top 10 down-regulated transcripts, Aggf1 encodes an angiogenic factor that promotes angiogenesis and the proliferation of endothelial cells [28,30]. The Cttn-encoded protein is located in the cytoplasm where the cell matrix is in contact. This protein regulates the interaction between components of the adhesive knot and is associated with cytoskeleton and cell adhesion structures essential for sphingosine-1-phosphate–mediated endothelial cell barrier enhancement [31,32]. Prom1 encodes a pentaspan transmembrane glycoprotein. This protein related to membrane processes is commonly expressed on adult stem cells and may influence cell differentiation, proliferation, and apoptosis [33]. Anln encodes an actin-binding protein that plays a role in cell growth, migration, and cytokinesis. The protein is thought to regulate actin cytoskeletal dynamics [34]. The protein encoded by Aggf1 has an effect of promoting endothelial cell proliferation, and may decrease the proliferation of BMECs when it is down-regulated after radiation irradiation [30]. Because they are involved in cytoskeletal composition and cell adhesion, the proteins encoded by Ltbp2, Cttn, Klhl15, Prom1, and Anln may change cell morphology after irradiation, resulting in decreased connectivity between cells [10,32,35–37]. Few studies have previously indicated that these genes are involved in the development and progression of radiation-induced brain edema. We believed that these genes may potentially play an important role in this process.

The results of our GO analysis suggested that radiation exposure increased apoptosis and vascular permeability, affected cell morphology, impeded the normal cell cycle, weakened cell proliferation, and impaired cell connections. All of these may affect vascular endothelial cell function and be involve in radiation-induced endothelial cell dysfunction.

The enriched KEGG pathways of the DEGs included the Ca2+ signaling pathway, PI3K–Akt signaling pathway, and methionine degradation. Under various pathological conditions, dysregulation of intracellular Ca2+ may induce sharp rises in the intracellular free Ca2+ ion concentration, leading to Ca2+ overload. Studies have shown that radiation (X-rays, gamma rays, etc.) can induce a large Ca2+ influx, intracellular Ca2+ overload, and cell dysfunction [38]. We found that irradiated BMECs showed changes in the Ca2+ signaling pathway, suggesting that radiation-induced intracellular Ca2+ dysregulation may disrupt the BBB integrity to cause or exacerbate brain edema.

We constructed and analyzed a PPI network encoded by differentially expressed transcripts. The following 10 closely related genes were identified: Cdkn1a, Ccl2, Ngf, Prkar2b, Cdc42, Hif1a, Rhog, Bub3, Vegfa, and Plk4. The proteins encoded by these genes were key nodes in the PPI network and may be the main proteins associated with brain edema. Our PPI network analysis also showed that Prkar2b, Cdkn1a, Vegfa, Ccl2, Bub3, and Ngf were up-regulated, while Hif1a, Cdc42, Plk4, and Rhog were down-regulated by radiation treatment.

There are very few reports showed the directly evidence that the Prkar2b, Plk4, and Bub3 involved in BMECs damage. However, Cdkn1a-encoded protein plays a key role in proliferation, migration and tube formation in response to hypoxia in endothelial cells [39]. In addition, others demonstrated that CDKN1A protein can inhibit cultured BMECs apoptosis [35]. Our study showed that Cdkn1a mRNA was significantly increased in the radiation treatment group, which suggests that increased Cdkn1a expression may be protective against X-ray-induced cell injury. The Vegfa gene is a member of the PDGF/VEGF growth factor family [40]. The Vegfa-encoded protein induces proliferation and migration of vascular endothelial cells and is essential for both physiological and pathological angiogenesis. It also has anti-angiogenic and vascular permeability-inducing functions [41,42]. We believe that up-regulated VEGFA protein may lead to increased microvascular permeability in the brain. During neuroinflammation, VEGF is elevated in the central nervous system, resulting in a change in the junction complex of the BBB to cause destruction of the BBB [43]. In addition, VEGF can down-regulate the expression of claudin-5 and induce the reorganization of F-actin microfilaments, which increases the centripetal tension within the cells, leading to cell retraction [44]. Ccl2 gene is also up-regulated by radiation treatment in BMECs. Block of CCL2 protein decreased the permeability of BBB in cultured endothelial cells during ischemia/reperfusion [36], while another study reported that CCL2 treatment alone did not disrupt the barrier sufficiently to change BBB permeability in BMECs [37]. Our data displayed the increased mRNA expression of CCL2 in the X-rays treatment group. But further study is needed to clarify the effect of CCL2 on X-ray-induced BMECs injury. In human dermal microvascular endothelial cells, one report showed that Ngf-encoded protein treatment enhanced cell proliferation [40], and other study showed that in hypoxic retinal endothelium, NGF attenuated cell apoptosis, and may be as a potential target for proliferative retinopathies [41]. Our data showed that Ngf mRNA expression was significantly up-regulated in the X-rays treatment group. Therefore, based on these evidences, we hypothesize that increased expression of NGF in the radiation treatment maybe has a protective effect on the BMECs. Moreover, one study reported that inhibition of Hif1a-encoded protein attenuated BBB damage in acute cerebral ischemia, and we also found that Hif1a mRNA was down-regulated by radiation treatment in the present study [42], while deficiency of Cdc42-encoded protein also displayed increased vascular permeability in vivo and activation of CDC42 might protect endothelial barrier [43,45]. Moreover, it is interesting that one study showed that Rhog-encoded protein played an essential role of angiogenesis in vascular endothelial cells via mediating the CDC42 [46]. Whereas we found that both Cdc42 and Rhog genes were down-regulated by the radiation treatment in the present study. Thus, the reduction of Cdc42 and Rhog genes expression may be one reason for X-ray-induced endothelial cell injury.


In summary, we used next-generation sequencing technology to determine the expression profile of the transcriptome in rat BMECs after radiation exposure to explore the mechanism underlying radiotherapy-induced BMEC damage. Our bioinformatics analysis revealed DEGs and changed signaling pathways in BMECs after radiation exposure and their significance in radiation-induced brain edema. These results point to several specific genes as being potential candidates of the molecular mechanisms underpinning radiation-induced brain edema and provide new potential targets for clinical treatment and prevention of radiation-induced cerebral edema.

Competing Interests

The authors declare that there are no competing interests associated with the manuscript.


This work was supported by Natural Science Foundation of China [grant numbers 81371284, 8197102295, 81960662 and U1732157]; Natural Science Foundation of Anhui Province [grant number 1808085MH246]; Anhui Province Science and Technology Innovation Project Demonstration Project [grant number 201707d08050003]; Yunnan Provincial Science and Technology Department [grant number 2019FE001(-133)]; the Open Research Foundation of Yunnan Key Laboratory of Bioactive Peptides in Yunnan Province [grant number AMHD-2017-3]; the Open Research Foundation of Yunnan Key Laboratory of Pharmacology for Nature Products [grant number 2019G001]; the Yunnan Provincial Education Department [grant number 2018JS163].

Author Contribution

Q.W., Y.F., X.Z., F.S., Y.W., H.C., J.D. and B.S. designed experiments and analyzed data. Q.W., Y.F., F.S., C.Z. and B.S. wrote the manuscript. J.D. and B.S. supervised the project. All of the authors read and approved the final manuscript.



blood–brain barrier


brain microvascular endothelial cell


differentially expressed genes




molecular complex detection


phosphoinositide 3-kinase


protein–protein interaction



Asai A., and Kawamoto K. (2008) . Radiation-induced brain injury. Brain Nerve60, , pp.123–129


Yeh S.A., , Tang Y., , Lui C.C., , Huang Y.J., and Huang E.Y. (2005) . Treatment outcomes and late complications of 849 patients with nasopharyngeal carcinoma treated with radiotherapy alone. Int. J. Radiat. Oncol. Biol. Phys.62, , pp.672–679, doi: 10.1016/j.ijrobp.2004.11.002


Balentova S., and Adamkov M. (2015) . Molecular, Cellular and Functional Effects of Radiation-Induced Brain Injury: A Review. Int. J. Mol. Sci.16, , pp.27796–27815, doi: 10.3390/ijms161126068


Jha S.K. (2003) . Cerebral Edema and its Management. Med. J. Armed. Forces India59, , pp.326–331, doi: 10.1016/S0377-1237(03)80147-8


Greene-Schloesser D., , Robbins M.E., , Peiffer A.M., , Shaw E.G., , Wheeler K.T., and Chan M.D. (2012) . Radiation-induced brain injury: A review. Front. Oncol.2, , pp.73, doi: 10.3389/fonc.2012.00073


Jiang J.D., , Ma Z.M., and Liu Y.S. (2008) . Expression of aquaporin-4 protein in rats with acute radiation-induced cerebral edema. Zhong Nan Da Xue Xue Bao Yi Xue Ban33, , pp.252–256


Greene-Schloesser D., , Moore E., and Robbins M.E. (2013) . Molecular pathways: radiation-induced cognitive impairment. Clin. Cancer Res.19, , pp.2294–2300, doi: 10.1158/1078-0432.CCR-11-2903


Pena L.A., , Fuks Z., and Kolesnick R.N. (2000) . Radiation-induced apoptosis of endothelial cells in the murine central nervous system: protection by fibroblast growth factor and sphingomyelinase deficiency. Cancer Res.60, , pp.321–327


Zhang Y., , Gao Z., , Wang D., , Zhang T., , Sun B., , Mu L., et al. (2014) . Accumulation of natural killer cells in ischemic brain tissues and the chemotactic effect of IP-10. J. Neuroinflammation11, , pp.79, doi: 10.1186/1742-2094-11-79


Franke W.W. (2009) . Discovering the molecular components of intercellular junctions–a historical view. Cold Spring Harb. Perspect. Biol.1, , pp.a003061, doi: 10.1101/cshperspect.a003061


Choi S.H., , Kim A.R., , Nam J.K., , Kim J.M., , Kim J.Y., , Seo H.R., et al. (2018) . Tumour-vasculature development via endothelial-to-mesenchymal transition after radiotherapy controls CD44v6(+) cancer cell and macrophage polarization. Nat. Commun.9, , pp.5108, doi: 10.1038/s41467-018-07470-w


Hu S.Y., , Duan H.F., , Li Q.F., , Yang Y.F., , Chen J.L., , Wang L.S., et al. (2009) . Hepatocyte growth factor protects endothelial cells against gamma ray irradiation-induced damage. Acta Pharmacol. Sin.30, , pp.1415–1420, doi: 10.1038/aps.2009.133


Hu S., , Gao Y., , Zhou H., , Kong F., , Xiao F., , Zhou P., et al. (2017) . New insight into mitochondrial changes in vascular endothelial cells irradiated by gamma ray. Int. J. Radiat. Biol.93, , pp.470–476, doi: 10.1080/09553002.2017.1286048


Bolger A.M., , Lohse M., and Usadel B. (2014) . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics30, , pp.2114–2120, doi: 10.1093/bioinformatics/btu170


Hwang S.G., , Kim K.H., , Lee B.M., and Moon J.C. (2018) . Transcriptome analysis for identifying possible gene regulations during maize root emergence and formation at the initial growth stage. Genes Genomics40, , pp.755–766, doi: 10.1007/s13258-018-0687-z


Peterson C.S., , Huang S., , Lee S.A., , Ferguson A.V., and Fry W.M. (2018) . The transcriptome of the rat subfornical organ is altered in response to early postnatal overnutrition. IBRO Rep.5, , pp.17–23, doi: 10.1016/j.ibror.2018.06.001


Pertea M., , Kim D., , Pertea G.M., , Leek J.T., and Salzberg S.L. (2016) . Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc.11, , pp.1650–1667, doi: 10.1038/nprot.2016.095


Yu Y., , Fuscoe J.C., , Zhao C., , Guo C., , Jia M., , Qing T., et al. (2014) . A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages. Nat. Commun.5, , pp.3230, doi: 10.1038/ncomms4230


Yang J., , Xu Z., , Gao J., , Liao C., , Wang P., , Liu Y., et al. (2018) . Evaluation of early acute radiation-induced brain injury: Hybrid multifunctional MRI-based study. Magn. Reson. Imaging54, , pp.101–108, doi: 10.1016/j.mri.2018.08.012


Bray N.L., , Pimentel H., , Melsted P., and Pachter L. (2016) . Erratum: Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol.34, , pp.888, doi: 10.1038/nbt0816-888d


Robinson M.D., , McCarthy D.J., and Smyth G.K. (2010) . edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics26, , pp.139–140, doi: 10.1093/bioinformatics/btp616


Park B.J., , Lord D., and Hart J.D. (2010) . Bias properties of Bayesian statistics in finite mixture of negative binomial regression models in crash data analysis. Accid. Anal. Prev.42, , pp.741–749, doi: 10.1016/j.aap.2009.11.002


von Mering C., , Huynen M., , Jaeggi D., , Schmidt S., , Bork P., and Snel B. (2003) . STRING: a database of predicted functional associations between proteins. Nucleic Acids Res.31, , pp.258–261, doi: 10.1093/nar/gkg034


Smoot M.E., , Ono K., , Ruscheinski J., , Wang P.L., and Ideker T. (2011) . Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics27, , pp.431–432, doi: 10.1093/bioinformatics/btq675


Lim J., , Hao T., , Shaw C., , Patel A.J., , Szabo G., , Rual J.F., et al. (2006) . A protein-protein interaction network for human inherited ataxias and disorders of Purkinje cell degeneration. Cell125, , pp.801–814, doi: 10.1016/j.cell.2006.03.032


Cao L., , Chen Y., , Zhang M., , Xu D.Q., , Liu Y., , Liu T., et al. (2018) . Identification of hub genes and potential molecular mechanisms in gastric cancer by integrated bioinformatics analysis. PeerJ6, , pp.e5180, doi: 10.7717/peerj.5180


He X., and Zhang J. (2006) . Why do hubs tend to be essential in protein networks?PLos Genet.2, , pp.e88, doi: 10.1371/journal.pgen.0020088


Zhu J., , Chen G., , Zhu S., , Li S., , Wen Z., , Bin L., et al. (2016) . Identification of Tissue-Specific Protein-Coding and Noncoding Transcripts across 14 Human Tissues Using RNA-seq. Sci. Rep.6, , pp.28400, doi: 10.1038/srep28400


Stelzl U., , Worm U., , Lalowski M., , Haenig C., , Brembeck F.H., , Goehler H., et al. (2005) . A human protein-protein interaction network: a resource for annotating the proteome. Cell122, , pp.957–968, doi: 10.1016/j.cell.2005.08.029


Tian X.L., , Kadaba R., , You S.A., , Liu M., , Timur A.A., , Yang L., et al. (2004) . Identification of an angiogenic factor that when mutated causes susceptibility to Klippel-Trenaunay syndrome. Nature427, , pp.640–645, doi: 10.1038/nature02320


von Holleben M., , Gohla A., , Janssen K.P., , Iritani B.M., and Beer-Hammer S. (2011) . Immunoinhibitory adapter protein Src homology domain 3 lymphocyte protein 2 (SLy2) regulates actin dynamics and B cell spreading. J. Biol. Chem.286, , pp.13489–13501, doi: 10.1074/jbc.M110.155184


Dudek S.M., , Chiang E.T., , Camp S.M., , Guo Y., , Zhao J., , Brown M.E., et al. (2010) . Abl tyrosine kinase phosphorylates nonmuscle Myosin light chain kinase to regulate endothelial barrier function. Mol. Biol. Cell21, , pp.4042–4056, doi: 10.1091/mbc.e09-10-0876


Macaire P., , Nadhari M., , Greiss H., , Godwin A., , Elhanfi O., , Sainudeen S., et al. (2014) . Internet remote control of pump settings for postoperative continuous peripheral nerve blocks: a feasibility study in 59 patients. Ann. Fr. Anesth. Reanim.33, , pp.e1–e7, doi: 10.1016/j.annfar.2013.11.017


Gbadegesin R.A., , Hall G., , Adeyemo A., , Hanke N., , Tossidou I., , Burchette J., et al. (2014) . Mutations in the gene that encodes the F-actin binding protein anillin cause FSGS. J. Am. Soc. Nephrol.25, , pp.1991–2002, doi: 10.1681/ASN.2013090976


Nayak G., , Odaka Y., , Prasad V., , Solano A.F., , Yeo E.J., , Vemaraju S., et al. (2018) . Developmental vascular regression is regulated by a Wnt/beta-catenin, MYC and CDKN1A pathway that controls cell proliferation and cell death. Development145, , pp.1–12, doi: 10.1242/dev.154898


Dimitrijevic O.B., , Stamatovic S.M., , Keep R.F., and Andjelkovic A.V. (2006) . Effects of the chemokine CCL2 on blood-brain barrier permeability during ischemia-reperfusion injury. J. Cerebral Blood Flow Metab.: Off. J. Int. Soc. Cerebral Blood Flow Metab.26, , pp.797–810, doi: 10.1038/sj.jcbfm.9600229


Roberts T.K., , Eugenin E.A., , Lopez L., , Romero I.A., , Weksler B.B., , Couraud P.O., et al. (2012) . CCL2 disrupts the adherens junction: implications for neuroinflammation. Lab. Invest.92, , pp.1213–1233, doi: 10.1038/labinvest.2012.80


Spielberg H., , June C.H., , Blair O.C., , Nystrom-Rosander C., , Cereb N., and Deeg H.J. (1991) . UV irradiation of lymphocytes triggers an increase in intracellular Ca2+ and prevents lectin-stimulated Ca2+ mobilization: evidence for UV- and nifedipine-sensitive Ca2+ channels. Exp. Hematol.19, , pp.742–748


Shu X., , Mao Y., , Li Z., , Wang W., , Chang Y., , Liu S., et al. (2019) . MicroRNA93 regulates angiogenesis in peripheral arterial disease by targeting CDKN1A. Mol. Med. Rep.19, , pp.5195–5202


Raychaudhuri S.K., , Raychaudhuri S.P., , Weltman H., and Farber E.M. (2001) . Effect of nerve growth factor on endothelial cell biology: proliferation and adherence molecule expression on human dermal microvascular endothelial cells. Arch. Dermatol. Res.293, , pp.291–295, doi: 10.1007/s004030100224


Troullinaki M., , Alexaki V.I., , Mitroulis I., , Witt A., , Klotzsche-von Ameln A., , Chung K.J., et al. (2019) . Nerve growth factor regulates endothelial cell survival and pathological retinal angiogenesis. J. Cell. Mol. Med.23, , pp.2362–2371, doi: 10.1111/jcmm.14002


Shen Y., , Gu J., , Liu Z., , Xu C., , Qian S., , Zhang X., et al. (2018) . Inhibition of HIF-1alpha Reduced Blood Brain Barrier Damage by Regulating MMP-2 and VEGF During Acute Cerebral Ischemia. Front. Cell. Neurosci.12, , pp.288, doi: 10.3389/fncel.2018.00288


Ramchandran R., , Mehta D., , Vogel S.M., , Mirza M.K., , Kouklis P., and Malik A.B. (2008) . Critical role of Cdc42 in mediating endothelial barrier protection in vivo. Am. J. Physiol. Lung Cell. Mol. Physiol.295, , pp.L363–L369, doi: 10.1152/ajplung.90241.2008


Awata T., , Inoue K., , Kurihara S., , Ohkubo T., , Watanabe M., , Inukai K., et al. (2002) . A common polymorphism in the 5′-untranslated region of the VEGF gene is associated with diabetic retinopathy in type 2 diabetes. Diabetes51, , pp.1635–1639, doi: 10.2337/diabetes.51.5.1635


Waschke J., , Burger S., , Curry F.R., , Drenckhahn D., and Adamson R.H. (2006) . Activation of Rac-1 and Cdc42 stabilizes the microvascular endothelial barrier. Histochem. Cell Biol.125, , pp.397–406, doi: 10.1007/s00418-005-0080-2


Abraham S., , Scarcia M., , Bagshaw R.D., , McMahon K., , Grant G., , Harvey T., et al. (2015) . A Rac/Cdc42 exchange factor complex promotes formation of lateral filopodia and blood vessel lumen morphogenesis. Nat. Commun.6, , pp.7286, doi: 10.1038/ncomms8286 of X-rays on transcript expression of rat brain microvascular endothelial cells: role of calcium signaling in X-ray-induced endothelium damage&author=Qibing Wu,Yang Fang,Xinchen Zhang,Fei Song,Yang Wang,Hongbo Chen,Juan Du,Chang-Bo Zheng,Bing Shen,&keyword=brain microvascular endothelial cell,Ca2+ signal,radioactive brain edema,RNA-Seq,&subject=Cardiovascular System & Vascular Biology,Pharmacology & Toxicology,Research Articles,