<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Neurosci.</journal-id>
<journal-title>Frontiers in Neuroscience</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Neurosci.</abbrev-journal-title>
<issn pub-type="epub">1662-453X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fnins.2022.834683</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Neuroscience</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Transcriptomic Analysis of Rat Cerebral Cortex Reveals the Potential Mechanism of Electroacupuncture Opening Blood Brain Barrier</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Ma</surname> <given-names>Congcong</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1593174/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Gan</surname> <given-names>Lin</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1095150/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Hao</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1094716/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Ren</surname> <given-names>Li</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1673066/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Lin</surname> <given-names>Yubo</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1673069/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhao</surname> <given-names>Yibin</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1094714/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhang</surname> <given-names>Shanshan</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1017254/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Gong</surname> <given-names>Peng</given-names></name>
<uri xlink:href="http://loop.frontiersin.org/people/1094734/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Lin</surname> <given-names>Xianming</given-names></name>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1673075/overview"/>
</contrib>
</contrib-group>
<aff><institution>Key Laboratory of Acupuncture and Neurology of Zhejiang Province, The Third Affiliated Hospital of Zhejiang Chinese Medical University, Zhejiang Chinese Medical University</institution>, <addr-line>Hangzhou</addr-line>, <country>China</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Jee-Yeon Hwang, Creighton University, United States</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Yujie Chen, Army Medical University, China; Mootaz M. Salman, University of Oxford, United Kingdom</p></fn>
<corresp id="c001">&#x002A;Correspondence: Xianming Lin, <email>linxianming1966@163.com</email></corresp>
<fn fn-type="other" id="fn004"><p>This article was submitted to Neurogenomics, a section of the journal Frontiers in Neuroscience</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>24</day>
<month>02</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>16</volume>
<elocation-id>834683</elocation-id>
<history>
<date date-type="received">
<day>13</day>
<month>12</month>
<year>2021</year>
</date>
<date date-type="accepted">
<day>02</day>
<month>02</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2022 Ma, Gan, Wang, Ren, Lin, Zhao, Zhang, Gong and Lin.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Ma, Gan, Wang, Ren, Lin, Zhao, Zhang, Gong and Lin</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Therapeutic treatment options for central nervous system (CNS) diseases are greatly limited by the blood-brain barrier (BBB). Electroacupuncture (EA) can be used to induce an increase in BBB permeability on rats, providing a potential approach for the delivery of drugs from the systemic circulation into the brain. However, there remains a large gap in our knowledge regarding the impact of EA on brain gene expression. This work is focused on investigating the transcriptional changes of rat cerebral cortex following EA and expression changes in genes and bioinformatic analysis was performed. We found that the potential mechanism of EA opening BBB involves receptor-mediated/carrier-mediated endocytosis (RMT/CMT), and related genes include solute carrier (SLC) transporter genes and ATP-binding cassette (ABC) transporter genes. The results also suggested that EA may affect the expression of tight junction (TJ) proteins in endothelial cells by affecting integrin binding, autophagy pathway and calcium signaling pathway, thus further affecting the permeability of blood-brain barrier. Our results provide a valuable resource that will guide mechanism research of EA opening BBB and other ways to mediate drug delivery into the brain.</p>
</abstract>
<kwd-group>
<kwd>electroacupuncture</kwd>
<kwd>BBB</kwd>
<kwd>transcriptomic</kwd>
<kwd>mechanism</kwd>
<kwd>bioinformatics</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="5"/>
<equation-count count="0"/>
<ref-count count="44"/>
<page-count count="16"/>
<word-count count="8117"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>Introduction</title>
<p>Vascular endothelial cells in the brain are closely linked to each other through various junction proteins, and interact with pericytes and astrocytes to form a special blood-brain barrier (BBB) barrier system. BBB strictly limits the entry of neurotoxic substances, inflammatory factors and immune cells in the blood into central nervous system (CNS), and excludes the metabolites and neurotoxic substances in CNS from the brain. However, BBB also prevents certain drugs and large-molecule therapeutics from entering the brain (<xref ref-type="bibr" rid="B5">Banks, 2016</xref>). Clearly, higher levels of therapeutic drugs in the brain will be beneficial, reducing dose requirements and providing greater therapeutic indices. Reflecting this demand, many brain delivery technologies have been developed and tested in laboratory and clinical studies, comprising invasive technologies such as direct brain injection, intrathecal brain deliver, intracerebral grafts and deep brain stimulation and non-invasive technologies including nanoparticulate systems, focused ultrasound, biological mechanisms and intranasal brain delivery (<xref ref-type="bibr" rid="B32">Terstappen et al., 2021</xref>). A number of factors contribute to the physical barrier of BBB and the major are endothelial transport systems and cellular junction molecules.</p>
<p>Our previous study found that 2/100 Hz Electroacupuncture (EA) stimulation increased BBB permeability, reduced ZO-1 and occludin levels, and induced ultrastructural changes in TJ morphology (<xref ref-type="bibr" rid="B40">Zhang et al., 2020</xref>), but the molecular and signaling mechanisms of BBB opening mediated by EA are still lacking. Therefore, in this manuscript, transcriptome sequencing combined with bioinformatics analysis was used to explore the effect of specific stimulation mode electroacupuncture on the transcriptome of rat cerebral cortex, focusing on the molecules and signaling pathways related to the blood-brain barrier, in order to reveal the potential mechanism of electroacupuncture opening the blood-brain barrier, and provide potential targets and new research ideas for mediating central nervous system therapeutic drugs into the brain.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2.SS1">
<title>Experimental Animals</title>
<p>Adult male Sprague-Dawley rats (250&#x2013;300 g) were purchased from Shanghai Laboratory Animal Center, Chinese Academy of Sciences and raised in the Laboratory Animal Center of Zhejiang Chinese Medical University. Rats were housed five per cage with <italic>ad libitum</italic> access to food and water at constant temperature (25 &#x00B1; 2&#x00B0;C) with a standard 12:12 h light-dark cycle. This study was approved by the Animal Protection and Use Committee of Zhejiang University of Traditional Chinese Medicine. All procedures were carried out in accordance with the National Institutes of Health Guide for Care and Use of Laboratory Animals.</p>
</sec>
<sec id="S2.SS2">
<title>Electroacupuncture Stimulation</title>
<p>Rats in EA groups were treated with disposable sterile acupuncture needles (Beijing Zhongyan Taihe Medical Device Co., Ltd., China) at GV20 (Baihui) and GV26 (Shuigou), GV20 with a needle of 25 mm in length and 0.13 mm in diameter and GV26 of a 16 mm in length and 0.07 mm in diameter. Then the needle was stimulated with a self-made relay (mode 6 s, stop 6 s) connected to the acupuncture point stimulator (HANS-200, Nanjing Jinsheng, Ltd., China). The current intensity was 3 mA, the frequency was 2/100 Hz, and the stimulation lasted for 40 min. Control rats were treated with the same binding for 40 min.</p>
</sec>
<sec id="S2.SS3">
<title>Extraction of Cerebral Cortex Tissue</title>
<p>After EA, rats were anesthetized with pentobarbital (50 mg/kg). The animals were perfused with 0.9% normal saline through the left ventricle until the colorless liquid was obtained in the right atrium, and the liver became white. The brain was cut off, and the cortex tissue was quickly stripped, cut into small pieces, and put into a 1.5 mL EP tube filled with RNA later prepared in advance, so that the sample was completely immersed in the liquid, and then stored overnight at 4&#x00B0;C.</p>
</sec>
<sec id="S2.SS4">
<title>RNA-Seq Library Establishment and RNA-Seq</title>
<p>Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) following the manufacturer&#x2019;s procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, United States). The RNA integrity was assessed by Bioanalyzer 2100 (Agilent, CA, United States) with RIN number &#x003E; 7.0, and confirmed by electrophoresis with denaturing agarose gel. Poly (A) RNA is purified from 1 &#x03BC;g total RNA using Dynabeads Oligo (dT)25-61005 (Thermo Fisher Scientific, CA, United States) using two rounds of purification. Then the poly(A) RNA was fragmented into small pieces using Magnesium RNA Fragmentation Module (NEB, cat.e6150, United States) under 94&#x00B0;C 5&#x2013;7 min. Then the cleaved RNA fragments were reverse-transcribed to create the cDNA by SuperScript&#x2122; II Reverse Transcriptase (Invitrogen, cat. 1896649, United States), which were next used to synthesize U-labeled second-stranded DNAs with E. coli DNA polymerase I (NEB, cat.m0209, United States), RNase H (NEB, cat.m0297, United States), and dUTP Solution (Thermo Fisher Scientific, cat.R0133, United States). An A-base is then added to the blunt ends of each strand, preparing them for ligation to the indexed adapters. Each adapter contains a T-base overhang for ligating the adapter to the A-tailed fragmented DNA. Single- or dual-index adapters are ligated to the fragments, and size selection was performed with AMPureXP beads. After the heat-labile UDG enzyme (NEB, cat.m0280, United States) treatment of the U-labeled second-stranded DNAs, the ligated products are amplified with PCR by the following conditions: initial denaturation at 95&#x00B0;C for 3 min; 8 cycles of denaturation at 98&#x00B0;C for 15 s, annealing at 60&#x00B0;C for 15 s, and extension at 72&#x00B0;C for 30 s; and then final extension at 72&#x00B0;C for 5 min. The average insert size for the final cDNA library was 300 &#x00B1; 50 bp. At last, we performed the 2 &#x00D7; 150bp paired-end sequencing (PE150) on an Illumina Novaseq&#x2122; 6000 (LC-Bio Technology Co., Ltd., Hangzhou, China) following the vendor&#x2019;s recommended protocol.</p>
</sec>
<sec id="S2.SS5">
<title>Data Analysis</title>
<sec id="S2.SS5.SSS1">
<title>Sequence and Primary Analysis</title>
<p>Cutadapt software (<sup><xref ref-type="fn" rid="footnote1">1</xref></sup> version:cutadapt-1.9) was used to remove the reads that contained adaptor contamination. After removed the low quality bases and undetermined bases, we used HISAT2 software (<sup><xref ref-type="fn" rid="footnote2">2</xref></sup> version:hisat2-2.0.4) to map reads to the genome. The mapped reads of each sample were assembled using StringTie (<sup><xref ref-type="fn" rid="footnote3">3</xref></sup> version:stringtie-1.3.4d. Linux_ x86 _64). Then, all transcriptomes from all samples were merged to reconstruct a comprehensive transcriptome using gffcompare software (<sup><xref ref-type="fn" rid="footnote4">4</xref></sup> version:gffcompare-0.9.8. Linux_x86_64). After the final transcriptome was generated, StringTie and ballgown<sup><xref ref-type="fn" rid="footnote5">5</xref></sup> were used to estimate the expression levels of all transcripts and perform expression level for mRNAs by calculating FPKM.</p>
<p>Primary sequencing data produced by RNA-Seq (raw reads) were subjected to quality control (QC). Specific data preprocessing steps are as follows: remove reads with adaptor; reads containing N (N means undetermined base information) with a proportion of more than 5% were removed; remove low quality reads (the base number of <italic>Q</italic> &#x2264; 10 accounts for more than 20% of the whole read). Statistics of original sequencing, effective sequencing, Q20, Q30, GC content and comprehensive evaluation. The information of total reads and mapping ratio reads were shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap position="float" id="T1">
<label>TABLE 1</label>
<caption><p>The primers used in qPCR.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Primers</td>
<td valign="top" align="left">Forward</td>
<td valign="top" align="left">Reverse</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><italic>Cttn</italic></td>
<td valign="top" align="left">GCAGAAGGATCGGATGGACAAGAAC</td>
<td valign="top" align="left">TGTACTCAGGCTCAGGCTCACTAC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Slc9a3r2</italic></td>
<td valign="top" align="left">CATGCCGAAGTTGTTGCCAGAATC</td>
<td valign="top" align="left">CACTGCCATCCTCATTGTCCTTCTC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Slc38a3</italic></td>
<td valign="top" align="left">ATGGTGGTGGTGGAGGAGAAGTC</td>
<td valign="top" align="left">ATAGAAGGTGAGGTAGCCGAAGAGG</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Slc7a1</italic></td>
<td valign="top" align="left">CCGTATCCGCTGCTCCATTGAAC</td>
<td valign="top" align="left">GACCCTCTGTGAACCTTAAACCCATC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Pigg</italic></td>
<td valign="top" align="left">GCGTAGCGATTGAGGAGTCAGTTC</td>
<td valign="top" align="left">ATGACATCAACAAAGCCAGGGAGAC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Itpr3</italic></td>
<td valign="top" align="left">CATTGGTGCGTCTGGAGGAACTG</td>
<td valign="top" align="left">CTCGGTCTTGGTGATGTGCTTCTC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Ltbp3</italic></td>
<td valign="top" align="left">AGGGCTACACTCAAGACAACAACATC</td>
<td valign="top" align="left">CTCATCCAAGCACTCATCCACATCC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Dnaja1</italic></td>
<td valign="top" align="left">CAACCGAACCATAGTCATCACCTCTC</td>
<td valign="top" align="left">GCTCATCATCCTCATACGCTTCTCC</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Cblb</italic></td>
<td valign="top" align="left">ACAGACGCCACGATTTGCCTTC</td>
<td valign="top" align="left">GACCATTATCACAAGACCGAACAGGAG</td>
</tr>
<tr>
<td valign="top" align="left"><italic>Hspa1b</italic></td>
<td valign="top" align="left">AAGATCACCATCACCAACGACAAGG</td>
<td valign="top" align="left">CCTCTTTCTCAGCCAGCGTGTTAG</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="S2.SS5.SSS2">
<title>Cluster Analysis and Screening of Differentially Expressed Genes</title>
<p>Clustering is calculated by using hclust function in stats package in R language, z-score normalization is performed by using scale function, and visualization is performed by using pheatmap package (<xref ref-type="bibr" rid="B25">Murtagh and Legendre, 2014</xref>). The differentially expressed genes (DEGs) were selected with fold change = 1.2 and <italic>q</italic>-Value &#x003C; 0.05 by R package DESeq2<sup><xref ref-type="fn" rid="footnote6">6</xref></sup>.</p>
</sec>
<sec id="S2.SS5.SSS3">
<title>Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Analysis</title>
<p>Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment to the DEGs was performed using the OmicStudio tools at <ext-link ext-link-type="uri" xlink:href="https://www.omicstudio.cn/tool">https://www.omicstudio.cn/tool</ext-link>. GO enrichment analysis can predict the functional role of DEGs from three aspects: biological process (BP), cellular component (CC) and molecular function (MF), and KEGG analysis can determine the related pathways of these genes. Next, DEGs were imported into Metascape<sup><xref ref-type="fn" rid="footnote7">7</xref></sup> for a network of enrichment terms (<xref ref-type="bibr" rid="B42">Zhou et al., 2019</xref>).</p>
</sec>
<sec id="S2.SS5.SSS4">
<title>Gene Set Enrichment Analysis</title>
<p>Considering the interaction between genes, we performed GSEA enrichment analysis to connect EA intervention with changes in biological function by comparing genes with predefined gene sets in MSigDB, analyzing gene expression data, and obtains whether expression is enriched in some function (<xref ref-type="bibr" rid="B29">Subramanian et al., 2005</xref>). The selected gene set is c2. cp. kegg. v7.4 in MSigDB, screening with NES &#x003E; 1, <italic>P</italic> &#x003C; 0.05 and FDR &#x003C; 0.25. Clusterprofiler package in R language was used for analysis and enrichplot package for visualization (<xref ref-type="bibr" rid="B38">Yu et al., 2012</xref>).</p>
</sec>
<sec id="S2.SS5.SSS5">
<title>Protein&#x2013;Protein Interaction Network Analysis</title>
<p>Interactions between proteins are central to all biological functions, and protein-protein interaction networks (PPINs) provide further insights into the protein function. Protein&#x2013;protein interactions (PPIs) are experimentally detected at a large scale and used to determine which proteins in a network directly interact. The Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database<sup><xref ref-type="fn" rid="footnote8">8</xref></sup> is used to provide information about protein prediction and experimental interaction by setting the combination score &#x003E; 0.4 as the reliability threshold, and the cytoscape software were used to construct protein&#x2013;protein interaction networks and screen the hub genes using cytohubba (<xref ref-type="bibr" rid="B28">Shannon et al., 2003</xref>; <xref ref-type="bibr" rid="B9">Chin et al., 2014</xref>). Functional modules within the network were identified with the MCODE app (<xref ref-type="bibr" rid="B2">Bader and Hogue, 2003</xref>).</p>
</sec>
</sec>
<sec id="S2.SS6">
<title>Real-Time Quantitative PCR Analysis</title>
<p>According to the manufacturer&#x2019;s procedure, total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, United States). The primer sequence is listed in <xref ref-type="table" rid="T1">Table 1</xref>, the total RNA extracted was reversely transcribed into cDNA, and &#x03B2;-actin was used as an internal reference gene, using random hexamer primers (Takara bioInc., Shiga, Japan). In accordance with instructions, through cfx96 real-time system (bio rad laboratory Co., Ltd., Hercules, CA, United States), fast start universal SYBR Green master kit (Takara bioInc., China) and 25 &#x03BC;L reaction system for qPCR. Each reaction was performed three times and normalized to &#x03B2;- Actin gene expression. Cfx96 real-time system software is used to determine the CT value of each well and calculate the average value of three times. Relative quantification was determined by &#x0394;&#x0394;CT method (<xref ref-type="bibr" rid="B22">Livak and Schmittgen, 2001</xref>).</p>
</sec>
<sec id="S2.SS7">
<title>Statistical Analysis</title>
<p>Data in graphs are expressed as means &#x00B1; SEM. Student&#x2019;s <italic>t</italic>-test was used for comparisons between two groups. Comparison is considered significantly different if the <italic>P</italic>-value &#x003C; 0.05.</p>
</sec>
</sec>
<sec id="S3" sec-type="results">
<title>Results</title>
<sec id="S3.SS1">
<title>RNA-Seq of Transcriptome</title>
<p>In order to explore the possible genes involved in the mechanism of EA opening blood brain barrier, we collected the cerebral cortex of rats in EA group and control group, and analyzed the gene expression profile by RNA-Seq. The Q20% of each RNA-Seq sample reached 99.0% (<xref ref-type="table" rid="T2">Table 2</xref>), and more than 90% of the base reached mass fraction &#x2265; Q30.</p>
<table-wrap position="float" id="T2">
<label>TABLE 2</label>
<caption><p>The information of raw reads and valid mapping ratio for electroacupuncture and control groups in RNA-Seq.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Sample</td>
<td valign="top" align="center">Raw read</td>
<td valign="top" align="center">Valid read</td>
<td valign="top" align="center">Valid base</td>
<td valign="top" align="center">Valid ratio</td>
<td valign="top" align="center">Q20%</td>
<td valign="top" align="center">Q30%</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">control_1</td>
<td valign="top" align="center">50386002</td>
<td valign="top" align="center">49081690</td>
<td valign="top" align="center">7.36G</td>
<td valign="top" align="center">97.41</td>
<td valign="top" align="center">99.97</td>
<td valign="top" align="center">99.05</td>
</tr>
<tr>
<td valign="top" align="left">control_2</td>
<td valign="top" align="center">50532582</td>
<td valign="top" align="center">49165486</td>
<td valign="top" align="center">7.37G</td>
<td valign="top" align="center">97.29</td>
<td valign="top" align="center">99.96</td>
<td valign="top" align="center">99.04</td>
</tr>
<tr>
<td valign="top" align="left">control_3</td>
<td valign="top" align="center">46976536</td>
<td valign="top" align="center">45706082</td>
<td valign="top" align="center">6.86G</td>
<td valign="top" align="center">97.30</td>
<td valign="top" align="center">99.96</td>
<td valign="top" align="center">99.00</td>
</tr>
<tr>
<td valign="top" align="left">control_4</td>
<td valign="top" align="center">46213950</td>
<td valign="top" align="center">45019482</td>
<td valign="top" align="center">6.75G</td>
<td valign="top" align="center">97.42</td>
<td valign="top" align="center">99.96</td>
<td valign="top" align="center">98.98</td>
</tr>
<tr>
<td valign="top" align="left">control_5</td>
<td valign="top" align="center">51654782</td>
<td valign="top" align="center">50154864</td>
<td valign="top" align="center">7.52G</td>
<td valign="top" align="center">97.10</td>
<td valign="top" align="center">99.96</td>
<td valign="top" align="center">99.00</td>
</tr>
<tr>
<td valign="top" align="left">EA_1</td>
<td valign="top" align="center">47846688</td>
<td valign="top" align="center">46200880</td>
<td valign="top" align="center">6.93G</td>
<td valign="top" align="center">96.56</td>
<td valign="top" align="center">99.96</td>
<td valign="top" align="center">99.06</td>
</tr>
<tr>
<td valign="top" align="left">EA_2</td>
<td valign="top" align="center">53084710</td>
<td valign="top" align="center">51424978</td>
<td valign="top" align="center">7.71G</td>
<td valign="top" align="center">96.87</td>
<td valign="top" align="center">99.97</td>
<td valign="top" align="center">99.08</td>
</tr>
<tr>
<td valign="top" align="left">EA_3</td>
<td valign="top" align="center">40298052</td>
<td valign="top" align="center">39082794</td>
<td valign="top" align="center">5.86G</td>
<td valign="top" align="center">96.98</td>
<td valign="top" align="center">99.97</td>
<td valign="top" align="center">99.18</td>
</tr>
<tr>
<td valign="top" align="left">EA_4</td>
<td valign="top" align="center">47508576</td>
<td valign="top" align="center">46006426</td>
<td valign="top" align="center">6.90G</td>
<td valign="top" align="center">96.84</td>
<td valign="top" align="center">99.97</td>
<td valign="top" align="center">99.05</td>
</tr>
<tr>
<td valign="top" align="left">EA_5</td>
<td valign="top" align="center">50579690</td>
<td valign="top" align="center">49280184</td>
<td valign="top" align="center">7.39G</td>
<td valign="top" align="center">97.43</td>
<td valign="top" align="center">99.97</td>
<td valign="top" align="center">99.05</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="S3.SS2">
<title>Identification of Differentially Expressed Genes</title>
<p>The 22310 genes were successfully mapped and identified from RNA-Seq. Then, we began to identify differentially expressed genes (DEGs) by FC &#x003E; 1.2 and <italic>q</italic> &#x003C; 0.05 and cluster the differentially expressed genes. The heat map showed that the differentially expressed genes were significantly separated between the EA group and the control group (<xref ref-type="fig" rid="F1">Figure 1</xref>). According to this standard, we identified a total of 156 DEGs (including 70 up-regulated and 86 down regulated), which were displayed in the volcanic map, marked the most significantly up-regulated 10 genes and the most significantly down regulated 10 genes (<xref ref-type="fig" rid="F2">Figure 2A</xref>), and analyzed the correlation of these genes (<xref ref-type="fig" rid="F2">Figure 2B</xref>). In addition, we described the expression of significantly up-regulated and down-regulated genes in the two groups (<xref ref-type="fig" rid="F2">Figures 2C,D</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>Cluster heatmap of differentially expressed genes.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g001.tif"/>
</fig>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Volcano plot of DEGs <bold>(A)</bold>, the correlation of the first 10 significantly up-regulated genes and significantly down-regulated genes <bold>(B)</bold>, the expression of up-regulated genes in electroacupuncture group and control group <bold>(C)</bold>, and the expression of down-regulated genes in electroacupuncture group and control group <bold>(D)</bold>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g002.tif"/>
</fig>
<p>The blood-brain barrier is not a single entity, but a series of characteristics, which together allow the endothelial cells of the central nervous system to strictly regulate the movement of ions, molecules and cells between blood and nerve tissue. Central nervous system endothelial cells are connected together by tight junction (TJ) to form a high resistance paracellular barrier (<xref ref-type="bibr" rid="B15">Furuse, 2010</xref>). Compared with peripheral endothelial cells, central nervous system endothelial cells show a lower cell transport rate (<xref ref-type="bibr" rid="B16">Hallmann et al., 1995</xref>). The expression of leukocyte adhesion molecules (LAMS) in endothelial cells of the central nervous system is also low, which limits the immune monitoring of the central nervous system (<xref ref-type="bibr" rid="B14">Engelhardt, 2008</xref>). Central nervous system endothelial cells express a variety of transporters, including a large number of special SLC transporters and ABC transporters. Genes related to blood-brain barrier were divided into groups for different BBB properties: tight junction integrity, SLC transporters, ABC transporters, LAMs, other transporters, transcytosis, endocytosis, angiogenesis or other BBB-enriched (<xref ref-type="table" rid="T3">Table 3</xref>). Our results showed that among these genes, the genes significantly up-regulated were SLC transporters related gene <italic>Slc25a25</italic>, endocytosis related genes <italic>Ccdc117</italic>, <italic>Cdc42ep3</italic>, <italic>Cblb</italic>, and <italic>Hspa1b</italic>, and the genes significantly down regulated were Slc transporter related genes <italic>Slc7a1</italic>, <italic>Slc9a3r2</italic>, <italic>Slc38a3</italic>, <italic>Slco2b1</italic>, ABC transporter related gene <italic>Abcc6</italic>, transcytosis related gene <italic>Cavin2</italic>, and endocytosis related genes <italic>Cttn</italic>, <italic>Cdc45</italic>, <italic>Hspa12b.</italic></p>
<table-wrap position="float" id="T3">
<label>TABLE 3</label>
<caption><p>Genes related to blood-brain barrier.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">BBB properties</td>
<td valign="top" align="center"><italic>Gene symbol</italic></td>
<td valign="top" align="center">Log2 FC (EA/control)</td>
<td valign="top" align="center">Adjusted <italic>P</italic>-value</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Tight junction integrity</td>
<td valign="top" align="center"><italic>Cgn</italic></td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.982</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cgnl1</italic></td>
<td valign="top" align="center">&#x2013;0.09</td>
<td valign="top" align="center">0.954</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cldn5</italic></td>
<td valign="top" align="center">&#x2013;11.59</td>
<td valign="top" align="center">0.306</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Ctnna1</italic></td>
<td valign="top" align="center">&#x2013;0.09</td>
<td valign="top" align="center">0.869</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Ctnna2</italic></td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">0.925</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Ctnnb1</italic></td>
<td valign="top" align="center">&#x2013;0.12</td>
<td valign="top" align="center">0.660</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Jaml</italic></td>
<td valign="top" align="center">0.22</td>
<td valign="top" align="center">0.949</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Jam2</italic></td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.980</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Jam3</italic></td>
<td valign="top" align="center">&#x2013;0.04</td>
<td valign="top" align="center">0.978</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Ocln</italic></td>
<td valign="top" align="center">&#x2013;0.47</td>
<td valign="top" align="center">0.236</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Tjp1</italic></td>
<td valign="top" align="center">&#x2013;0.05</td>
<td valign="top" align="center">0.952</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Tjp2</italic></td>
<td valign="top" align="center">&#x2013;0.08</td>
<td valign="top" align="center">0.928</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Tjp3</italic></td>
<td valign="top" align="center">0.25</td>
<td valign="top" align="center">0.889</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Marveld1</italic></td>
<td valign="top" align="center">&#x2013;0.21</td>
<td valign="top" align="center">0.830</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Marveld2</italic></td>
<td valign="top" align="center">&#x2013;0.49</td>
<td valign="top" align="center">0.715</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Marveld3</italic></td>
<td valign="top" align="center">5.93</td>
<td valign="top" align="center">0.935</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Lsr</italic></td>
<td valign="top" align="center">&#x2013;0.34</td>
<td valign="top" align="center">0.268</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>F11r</italic></td>
<td valign="top" align="center">&#x2013;0.12</td>
<td valign="top" align="center">0.899</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Emp1</italic></td>
<td valign="top" align="center">&#x2013;0.01</td>
<td valign="top" align="center">0.996</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cgnl1</italic></td>
<td valign="top" align="center">&#x2013;0.09</td>
<td valign="top" align="center">0.954</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Pmp22</italic></td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.955</td>
</tr>
<tr>
<td valign="top" align="left">Slc transporters</td>
<td valign="top" align="center"><bold><italic>Slc7a1</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.80</bold></td>
<td valign="top" align="center"><bold>0.033</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Slc9a3r2</italic></bold></td>
<td valign="top" align="center">&#x2013;0.44</td>
<td valign="top" align="center"><bold>0.022</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Slc25a25</italic></bold></td>
<td valign="top" align="center"><bold>0.31</bold></td>
<td valign="top" align="center"><bold>0.013</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Slc38a3</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.29</bold></td>
<td valign="top" align="center"><bold>0.002</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Slco2b1</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.29</bold></td>
<td valign="top" align="center"><bold>0.010</bold></td>
</tr>
<tr>
<td valign="top" align="left">ABC transporters</td>
<td valign="top" align="center"><italic>Abca2</italic></td>
<td valign="top" align="center">&#x2013;0.23</td>
<td valign="top" align="center">0.022</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Abcb4</italic></td>
<td valign="top" align="center">&#x2013;0.24</td>
<td valign="top" align="center">0.116</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Abcc3</italic></td>
<td valign="top" align="center">&#x2013;0.57</td>
<td valign="top" align="center">0.148</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Abcc6</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.66</bold></td>
<td valign="top" align="center"><bold>0.015</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Abcd4</italic></td>
<td valign="top" align="center">&#x2013;0.14</td>
<td valign="top" align="center">0.924</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Abcg2</italic></td>
<td valign="top" align="center">&#x2013;0.25</td>
<td valign="top" align="center">0.658</td>
</tr>
<tr>
<td valign="top" align="left">LAMs</td>
<td valign="top" align="center"><italic>Sele</italic></td>
<td valign="top" align="center">&#x2013;1.26</td>
<td valign="top" align="center">0.836</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Selp</italic></td>
<td valign="top" align="center">0.32</td>
<td valign="top" align="center">0.979</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Icam1</italic></td>
<td valign="top" align="center">&#x2013;0.20</td>
<td valign="top" align="center">0.877</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Vcam1</italic></td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.980</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Alcam</italic></td>
<td valign="top" align="center">0.07</td>
<td valign="top" align="center">0.594</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Mcam</italic></td>
<td valign="top" align="center">&#x2013;0.07</td>
<td valign="top" align="center">0.962</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Ninj1</italic></td>
<td valign="top" align="center">&#x2013;0.05</td>
<td valign="top" align="center">0.980</td>
</tr>
<tr>
<td valign="top" align="left">Other transporters</td>
<td valign="top" align="center"><italic>Lrp1</italic></td>
<td valign="top" align="center">&#x2013;0.18</td>
<td valign="top" align="center">0.188</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Lrp10</italic></td>
<td valign="top" align="center">&#x2013;0.09</td>
<td valign="top" align="center">0.926</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Mfsd2a</italic></td>
<td valign="top" align="center">&#x2013;0.38</td>
<td valign="top" align="center">0.286</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Mfsd2b</italic></td>
<td valign="top" align="center">&#x2013;8.26</td>
<td valign="top" align="center">0.777</td>
</tr>
<tr>
<td valign="top" align="left">Transcytosis</td>
<td valign="top" align="center"><italic>Tfrc</italic></td>
<td valign="top" align="center">0.10</td>
<td valign="top" align="center">0.602</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cav1</italic></td>
<td valign="top" align="center">&#x2013;0.03</td>
<td valign="top" align="center">0.998</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cav2</italic></td>
<td valign="top" align="center">0.06</td>
<td valign="top" align="center">0.936</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cav3</italic></td>
<td valign="top" align="center">&#x2013;8.05</td>
<td valign="top" align="center">0.835</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cavin1</italic></td>
<td valign="top" align="center">&#x2013;0.42</td>
<td valign="top" align="center">0.659</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Cavin2</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.62</bold></td>
<td valign="top" align="center"><bold>0.035</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cavin3</italic></td>
<td valign="top" align="center">&#x2013;0.25</td>
<td valign="top" align="center">0.381</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Cavin4</italic></td>
<td valign="top" align="center">&#x2013;0.41</td>
<td valign="top" align="center">0.863</td>
</tr>
<tr>
<td valign="top" align="left">Endocytosis</td>
<td valign="top" align="center"><italic>Arf1</italic></td>
<td valign="top" align="center">&#x2013;0.02</td>
<td valign="top" align="center">0.989</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Arf2</italic></td>
<td valign="top" align="center">0.06</td>
<td valign="top" align="center">0.863</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Arf4</italic></td>
<td valign="top" align="center">0.13</td>
<td valign="top" align="center">0.349</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Arf5</italic></td>
<td valign="top" align="center">&#x2013;0.07</td>
<td valign="top" align="center">0.935</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Cttn</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>12.38</bold></td>
<td valign="top" align="center"><bold>0.000</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Ccdc117</italic></bold></td>
<td valign="top" align="center"><bold>0.42</bold></td>
<td valign="top" align="center"><bold>0.005</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Cdc42ep3</italic></bold></td>
<td valign="top" align="center"><bold>0.37</bold></td>
<td valign="top" align="center"><bold>0.008</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Cdc45</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.37</bold></td>
<td valign="top" align="center"><bold>0.033</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Cblb</italic></bold></td>
<td valign="top" align="center"><bold>0.74</bold></td>
<td valign="top" align="center"><bold>0.013</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Dnm1</italic></td>
<td valign="top" align="center">&#x2013;0.03</td>
<td valign="top" align="center">0.992</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Eps8</italic></td>
<td valign="top" align="center">&#x2013;0.17</td>
<td valign="top" align="center">0.374</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Eps15</italic></td>
<td valign="top" align="center">&#x2013;0.02</td>
<td valign="top" align="center">0.992</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Flot1</italic></td>
<td valign="top" align="center">&#x2013;0.12</td>
<td valign="top" align="center">0.879</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Flot2</italic></td>
<td valign="top" align="center">0.29</td>
<td valign="top" align="center">0.101</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Fyn</italic></td>
<td valign="top" align="center">0.12</td>
<td valign="top" align="center">0.352</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Hspa1b</italic></bold></td>
<td valign="top" align="center"><bold>2.67</bold></td>
<td valign="top" align="center"><bold>0.000</bold></td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><bold><italic>Hspa12b</italic></bold></td>
<td valign="top" align="center">&#x2013;<bold>0.64</bold></td>
<td valign="top" align="center"><bold>0.041</bold></td>
</tr>
<tr>
<td valign="top" align="left">Angiogenesis</td>
<td valign="top" align="center"><italic>Angpt4</italic></td>
<td valign="top" align="center">0.37</td>
<td valign="top" align="center">0.879</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Angpt1</italic></td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.980</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Angpt2</italic></td>
<td valign="top" align="center">&#x2013;0.10</td>
<td valign="top" align="center">0.970</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Egr1</italic></td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.961</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Egr2</italic></td>
<td valign="top" align="center">0.13</td>
<td valign="top" align="center">0.861</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Egr3</italic></td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.935</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Egr4</italic></td>
<td valign="top" align="center">&#x2013;0.16</td>
<td valign="top" align="center">0.886</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Fgfr1</italic></td>
<td valign="top" align="center">&#x2013;0.13</td>
<td valign="top" align="center">0.507</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Fgfr2</italic></td>
<td valign="top" align="center">&#x2013;0.04</td>
<td valign="top" align="center">0.995</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Itgb1</italic></td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.922</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Itgb2</italic></td>
<td valign="top" align="center">0.04</td>
<td valign="top" align="center">0.954</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Itgb3</italic></td>
<td valign="top" align="center">&#x2013;0.29</td>
<td valign="top" align="center">0.935</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Lgals1</italic></td>
<td valign="top" align="center">&#x2013;0.10</td>
<td valign="top" align="center">0.961</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Lgals2</italic></td>
<td valign="top" align="center">&#x2013;7.70</td>
<td valign="top" align="center">0.935</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Lgals3</italic></td>
<td valign="top" align="center">&#x2013;0.14</td>
<td valign="top" align="center">0.968</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Mmp9</italic></td>
<td valign="top" align="center">&#x2013;0.08</td>
<td valign="top" align="center">0.980</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Pdgfa</italic></td>
<td valign="top" align="center">&#x2013;0.07</td>
<td valign="top" align="center">0.953</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Pdgfb</italic></td>
<td valign="top" align="center">&#x2013;0.05</td>
<td valign="top" align="center">0.980</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Tgfb2</italic></td>
<td valign="top" align="center">0.14</td>
<td valign="top" align="center">0.692</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Vegfa</italic></td>
<td valign="top" align="center">&#x2013;0.13</td>
<td valign="top" align="center">0.949</td>
</tr>
<tr>
<td/>
<td valign="top" align="center"><italic>Vegfb</italic></td>
<td valign="top" align="center">&#x2013;0.31</td>
<td valign="top" align="center">0.605</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn><p><italic>Genes in bold refer to differentially expressed genes (DEGs) identified by FC &#x003E; 1.2 and q &#x003C; 0.05.</italic></p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="S3.SS3">
<title>Differential Expression Detected by Transcriptome Was Strongly Related to qPCR Results</title>
<p>Next, we set to examine the reliability of our RNA-Seq data using qPCR. We randomly selected three up-regulated genes (<italic>Hspa1b, Cblb, Dnaja1</italic>), two down-regulated genes (<italic>Itpr3, Ltbp3</italic>) and four genes related to BBB (<italic>Cttn, Abcc6, Slc9a3r2, Slc38a3</italic>) for qPCR verification. The results of qPCR showed that the expression of <italic>Hspa1b, Cblb</italic> and <italic>Dnaja1</italic> were all significantly up-regulated (<xref ref-type="fig" rid="F3">Figure 3B</xref>), whereas <italic>Itpr3</italic> and <italic>Ltbp3</italic> were significantly downregulated in EA rats vs. control rats (<xref ref-type="fig" rid="F3">Figure 3A</xref>), which are consistent with RNA-Seq data. Furthermore, the qPCR results of the four selected genes of interest were also consistent with the sequencing results (<xref ref-type="fig" rid="F3">Figure 3C</xref>). Thus, qPCR results provide evidence that RNA-Seq data for gene expression profiling are reliable.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>RNA-Seq results were validated by qPCR. Two randomly selected down-regulated DEG expression <bold>(A)</bold> from RNA-Seq and three up-regulated DEG <bold>(B)</bold> expression were detected by qPCR, and <bold>(C)</bold> four typical genes related to BBB. &#x002A;<italic>p</italic> &#x003C; 0.05, &#x002A;&#x002A;<italic>p</italic> &#x003C; 0.01, and &#x002A;&#x002A;&#x002A;<italic>p</italic> &#x003C; 0.00. Student&#x2019;s <italic>t</italic>-test was used for comparisons.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g003.tif"/>
</fig>
</sec>
<sec id="S3.SS4">
<title>Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Analysis of Differentially Expressed Genes</title>
<p>To analyze molecular mechanisms and signaling pathways further, we performed GO Term and KEGG pathway enrichment analysis on DEGs. As shown in <xref ref-type="fig" rid="F4">Figure 4</xref>, we show the GO items with top 10 counts and significant enrichment (<italic>P</italic> &#x003C; 0.05) in biological processes (BP), cell composition (CC) and molecular function (MF), respectively (<xref ref-type="fig" rid="F4">Figures 4A-C</xref>). Among them, the BP that may be related to the composition of BBB and material transport are ion transport (GO:0006811), transmembrane transport (GO:0055085), angiogenesis (GO:0001525), brain development (GO:0007420), potassium ion transport (GO:0006813); the CC are apical plasma membrane (GO:0016324), basolateral plasma membrane (GO:0016323), collagen-containing extracellular matrix (GO:0062023), focal adhesion (GO:0005925), caveola (GO:0005901), clathrin-coated pit (GO:0005905); the MF are ion channel activity (GO:0005216), potassium channel activity (GO:0005267). Then we list the statistics of the top 20 terms with the most significant enrichment (minimum <italic>p-</italic>Value) (<xref ref-type="fig" rid="F4">Figure 4D</xref>). Among these GO terms, what may be related to BBB are potassium channel activity (GO:0005267), angiogenesis (GO:0001525), calcium-release channel activity (GO:0015278), release of sequestered calcium ion into cytosol (GO:0051209), calcium ion transport into cytosol (GO:0060402), endothelial cell chemotaxis (GO:0035767), L-histidine transmembrane transporter activity (GO:0005290). All GO terms associated with the BBB are listed in <xref ref-type="table" rid="T4">Table 4</xref>.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Gene ontology items with top 10 counts and significant enrichment (<italic>P</italic> &#x003C; 0.05) in BP <bold>(A)</bold>, CC <bold>(B)</bold>, and MF <bold>(C)</bold>, and the statistics of the top 20 terms with the most significant enrichment <bold>(D)</bold>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g004.tif"/>
</fig>
<table-wrap position="float" id="T4">
<label>TABLE 4</label>
<caption><p>Details of all gene ontology terms related to blood-brain barrier.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">GO function</td>
<td valign="top" align="center">GO ID</td>
<td valign="top" align="left">GO term</td>
<td valign="top" align="center">DEGs</td>
<td valign="top" align="center"><italic>P</italic>-value</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">BP</td>
<td valign="top" align="center">GO:0006811</td>
<td valign="top" align="left">ion transport</td>
<td valign="top" align="center"><italic>Kcnt2, Ryr2, Kcng4, Kcnh4, Slco2b1, Kcnk4, Itpr3, Grin2c, Slc38a3</italic></td>
<td valign="top" align="center">0.0056</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0055085</td>
<td valign="top" align="left">transmembrane transport</td>
<td valign="top" align="center"><italic>Abcc6, Itpr3, Ryr2, Slco2b1, Kcnh4, Kcng4, Slc25a25, Slc7a1</italic></td>
<td valign="top" align="center">0.0248</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0001525</td>
<td valign="top" align="left">angiogenesis</td>
<td valign="top" align="center"><italic>Plcd3, Scg2, Sox18, Efna1, Ccn3, Tie1</italic></td>
<td valign="top" align="center">0.0014</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0007420</td>
<td valign="top" align="left">brain development</td>
<td valign="top" align="center"><italic>Rgs2, Cttn, Cxcl12, Slc38a3, Sphk2</italic></td>
<td valign="top" align="center">0.0251</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0006813</td>
<td valign="top" align="left">potassium ion transport</td>
<td valign="top" align="center"><italic>Kcnt2, Kcnh4, Kcng4, Kcnk4</italic></td>
<td valign="top" align="center">0.0087</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0071805</td>
<td valign="top" align="left">potassium ion transmembrane transport</td>
<td valign="top" align="center"><italic>Kcnt2, Kcnk4, Kcnh4, Kcng4</italic></td>
<td valign="top" align="center">0.0092</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0001568</td>
<td valign="top" align="left">blood vessel development</td>
<td valign="top" align="center"><italic>Tie1, Sox18, Sphk2</italic></td>
<td valign="top" align="center">0.0080</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0001938</td>
<td valign="top" align="left">positive regulation of endothelial cell proliferation</td>
<td valign="top" align="center"><italic>Cxcl12, Scg2, Cyba</italic></td>
<td valign="top" align="center">0.0090</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0070588</td>
<td valign="top" align="left">calcium ion transmembrane transport</td>
<td valign="top" align="center"><italic>Ryr2, Slc25a25, Itpr3</italic></td>
<td valign="top" align="center">0.0353</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0035767</td>
<td valign="top" align="left">endothelial cell chemotaxis</td>
<td valign="top" align="center"><italic>Coro1b, Ccn3</italic></td>
<td valign="top" align="center">0.0008</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0060402</td>
<td valign="top" align="left">calcium ion transport into cytosol</td>
<td valign="top" align="center"><italic>Ryr2, Itpr3</italic></td>
<td valign="top" align="center">0.0011</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0090280</td>
<td valign="top" align="left">positive regulation of calcium ion import</td>
<td valign="top" align="center"><italic>Cxcl12, Sphk2</italic></td>
<td valign="top" align="center">0.0044</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0006865</td>
<td valign="top" align="left">amino acid transport</td>
<td valign="top" align="center"><italic>Slc7a1, Slc38a3</italic></td>
<td valign="top" align="center">0.0339</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1903810</td>
<td valign="top" align="left">L-histidine import across plasma membrane</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0063</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:2000487</td>
<td valign="top" align="left">positive regulation of glutamine transport</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0063</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015817</td>
<td valign="top" align="left">histidine transport</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0063</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0051284</td>
<td valign="top" align="left">positive regulation of sequestering of calcium ion</td>
<td valign="top" align="center"><italic>Ryr2</italic></td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1901726</td>
<td valign="top" align="left">negative regulation of histone deacetylase activity</td>
<td valign="top" align="center"><italic>Sphk2</italic></td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0006867</td>
<td valign="top" align="left">asparagine transport</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1900924</td>
<td valign="top" align="left">negative regulation of glycine import across plasma membrane</td>
<td valign="top" align="center"><italic>Rgs2</italic></td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0060401</td>
<td valign="top" align="left">cytosolic calcium ion transport</td>
<td valign="top" align="center"><italic>Ryr2</italic></td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015819</td>
<td valign="top" align="left">lysine transport</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015822</td>
<td valign="top" align="left">ornithine transport</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0187</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0071603</td>
<td valign="top" align="left">endothelial cell-cell adhesion</td>
<td valign="top" align="center"><italic>Ccn3</italic></td>
<td valign="top" align="center">0.0187</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1903826</td>
<td valign="top" align="left">arginine transmembrane transport</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0187</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1901896</td>
<td valign="top" align="left">positive regulation of ATPase-coupled calcium transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Ryr2</italic></td>
<td valign="top" align="center">0.0249</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0030947</td>
<td valign="top" align="left">regulation of vascular endothelial growth factor receptor signaling pathway</td>
<td valign="top" align="center"><italic>Tmem204</italic></td>
<td valign="top" align="center">0.0249</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0089709</td>
<td valign="top" align="left">L-histidine transmembrane transport</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0249</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015809</td>
<td valign="top" align="left">arginine transport</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0043535</td>
<td valign="top" align="left">regulation of blood vessel endothelial cell migration</td>
<td valign="top" align="center"><italic>Efna1</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1990822</td>
<td valign="top" align="left">basic amino acid transmembrane transport</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0071896</td>
<td valign="top" align="left">protein localization to adherence junction</td>
<td valign="top" align="center"><italic>Scrib</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0097553</td>
<td valign="top" align="left">calcium ion transmembrane import into cytosol</td>
<td valign="top" align="center"><italic>Grin2c</italic></td>
<td valign="top" align="center">0.0371</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015807</td>
<td valign="top" align="left">L-amino acid transport</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0371</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0006868</td>
<td valign="top" align="left">glutamine transport</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0371</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015808</td>
<td valign="top" align="left">L-alanine transport</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0432</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:1903348</td>
<td valign="top" align="left">positive regulation of bicellular tight junction assembly</td>
<td valign="top" align="center"><italic>Nphp4</italic></td>
<td valign="top" align="center">0.0432</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0043615</td>
<td valign="top" align="left">astrocyte cell migration</td>
<td valign="top" align="center"><italic>Scrib</italic></td>
<td valign="top" align="center">0.0432</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0014808</td>
<td valign="top" align="left">release of sequestered calcium ion into cytosol by sarcoplasmic reticulum</td>
<td valign="top" align="center"><italic>Ryr2</italic></td>
<td valign="top" align="center">0.0492</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0038203</td>
<td valign="top" align="left">TORC2 signaling</td>
<td valign="top" align="center"><italic>Prr5l</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td valign="top" align="left">CC</td>
<td valign="top" align="center">GO:0016323</td>
<td valign="top" align="left">basolateral plasma membrane</td>
<td valign="top" align="center"><italic>Slc38a3, Slc7a1, Abcc6, Hspa1b, Scrib</italic></td>
<td valign="top" align="center">0.0087</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0062023</td>
<td valign="top" align="left">collagen-containing extracellular matrix</td>
<td valign="top" align="center"><italic>Mgp, Vwf, Ltbp3, Ccn3</italic></td>
<td valign="top" align="center">0.0072</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005901</td>
<td valign="top" align="left">caveola</td>
<td valign="top" align="center"><italic>Lipe, Cavin2, Htr2a</italic></td>
<td valign="top" align="center">0.0083</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005905</td>
<td valign="top" align="left">clathrin-coated pit</td>
<td valign="top" align="center"><italic>Cttn, Sh3bp4</italic></td>
<td valign="top" align="center">0.0424</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0062023</td>
<td valign="top" align="left">collagen-containing extracellular matrix</td>
<td valign="top" align="center"><italic>Mgp, Vwf, Ltbp3, Ccn3</italic></td>
<td valign="top" align="center">0.0071</td>
</tr>
<tr>
<td valign="top" align="left">MF</td>
<td valign="top" align="center">GO:0005216</td>
<td valign="top" align="left">ion channel activity</td>
<td valign="top" align="center"><italic>Itpr3, Kcnh4, Grin2c, Kcng4, Ryr2</italic></td>
<td valign="top" align="center">0.0058</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005267</td>
<td valign="top" align="left">potassium channel activity</td>
<td valign="top" align="center"><italic>Kcnk4, Kcnh4, Kcnt2, Kcng4</italic></td>
<td valign="top" align="center">0.0014</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005290</td>
<td valign="top" align="left">L-histidine transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc7a1, Slc38a3</italic></td>
<td valign="top" align="center">0.0002</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015278</td>
<td valign="top" align="left">calcium-release channel activity</td>
<td valign="top" align="center"><italic>Itpr3, Ryr2</italic></td>
<td valign="top" align="center">0.0014</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015171</td>
<td valign="top" align="left">amino acid transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc38a3, Slc7a1</italic></td>
<td valign="top" align="center">0.0133</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0098782</td>
<td valign="top" align="left">mechanosensitived potassium channel activity</td>
<td valign="top" align="center"><italic>Kcnk4</italic></td>
<td valign="top" align="center">0.0063</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005228</td>
<td valign="top" align="left">intracellular sodium activated potassium channel activity</td>
<td valign="top" align="center"><italic>Kcnt2</italic></td>
<td valign="top" align="center">0.0063</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0070089</td>
<td valign="top" align="left">chloride-activated potassium channel activity</td>
<td valign="top" align="center"><italic>Kcnt2</italic></td>
<td valign="top" align="center">0.0063</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015182</td>
<td valign="top" align="left">L-asparagine transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0187</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005219</td>
<td valign="top" align="left">ryanodine-sensitive calcium-release channel activity</td>
<td valign="top" align="center"><italic>Ryr2</italic></td>
<td valign="top" align="center">0.0187</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0048763</td>
<td valign="top" align="left">calcium-induced calcium release activity</td>
<td valign="top" align="center"><italic>Ryr2</italic></td>
<td valign="top" align="center">0.0249</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005220</td>
<td valign="top" align="left">inositol 1,4,5-trisphosphate-sensitive calcium-release channel activity</td>
<td valign="top" align="center"><italic>Itpr3</italic></td>
<td valign="top" align="center">0.0249</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0022849</td>
<td valign="top" align="left">glutamate-gated calcium ion channel activity</td>
<td valign="top" align="center"><italic>Grin2c</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015181</td>
<td valign="top" align="left">arginine transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0311</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015186</td>
<td valign="top" align="left">L-glutamine transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0371</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015174</td>
<td valign="top" align="left">basic amino acid transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc7a1</italic></td>
<td valign="top" align="center">0.0371</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0004972</td>
<td valign="top" align="left">NMDA glutamate receptor activity</td>
<td valign="top" align="center"><italic>Grin2c</italic></td>
<td valign="top" align="center">0.0371</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015180</td>
<td valign="top" align="left">L-alanine transmembrane transporter activity</td>
<td valign="top" align="center"><italic>Slc38a3</italic></td>
<td valign="top" align="center">0.0492</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0015271</td>
<td valign="top" align="left">outward rectifier potassium channel activity</td>
<td valign="top" align="center"><italic>Kcnt2</italic></td>
<td valign="top" align="center">0.0492</td>
</tr>
<tr>
<td/>
<td valign="top" align="center">GO:0005178</td>
<td valign="top" align="left">integrin binding</td>
<td valign="top" align="center"><italic>Vwf, Ccn3, Cxcl12</italic></td>
<td valign="top" align="center">0.0426</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Kyoto Encyclopedia of Genes and Genomes results predicted the signal pathways related to DEGs (<xref ref-type="fig" rid="F5">Figure 5A</xref>). There were 7 significantly enriched (<italic>P</italic> &#x003C; 0.05) pathways (<xref ref-type="fig" rid="F5">Figure 5B</xref>) and 15 pathways with more than or equal to 3 DEGs (<xref ref-type="table" rid="T5">Table 5</xref>). Calcium signaling pathway (4020) may be related to BBB under the intervention of EA and the DEGs involved in this signaling pathway are <italic>Plcd3</italic>, <italic>Ryr2</italic>, <italic>Itpr3</italic>, <italic>Htr2a</italic>, <italic>Grin2c</italic>, and <italic>Sphk2.</italic> In addition, apelin signaling pathway (4371) and platelet activation (4611) signaling pathway are also associated with EA intervention, involved DEGs are <italic>Itpr3, Lipe, Ryr2, Sphk2</italic> in the former and <italic>Gp1bb, Vwf, Itpr3</italic> in the later.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>Top 20 Kyoto Encyclopedia of Genes and Genomes enrichment pathways <bold>(A)</bold> and the significant enrichment pathways <bold>(B)</bold>, and the network of enrichment terms colored by cluster ID <bold>(C)</bold>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g005.tif"/>
</fig>
<table-wrap position="float" id="T5">
<label>TABLE 5</label>
<caption><p>Enrichment pathways with more than or equal to three differentially expressed genes.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<td valign="top" align="left">Pathway ID</td>
<td valign="top" align="left">Pathway name</td>
<td valign="top" align="left">Genes</td>
<td valign="top" align="center"><italic>P</italic>-value</td>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">4020</td>
<td valign="top" align="left">Calcium signaling pathway</td>
<td valign="top" align="left"><italic>Plcd3, Ryr2, Itpr3, Htr2a, Grin2c, Sphk2</italic></td>
<td valign="top" align="center">0.00</td>
</tr>
<tr>
<td valign="top" align="left">4141</td>
<td valign="top" align="left">Protein processing in endoplasmic reticulum</td>
<td valign="top" align="left"><italic>Hsph1, Hspa1b, Dnaja1, Dnajb1, Ppp1r15a</italic></td>
<td valign="top" align="center">0.00</td>
</tr>
<tr>
<td valign="top" align="left">4371</td>
<td valign="top" align="left">Apelin signaling pathway</td>
<td valign="top" align="left"><italic>Itpr3, Lipe, Ryr2, Sphk2</italic></td>
<td valign="top" align="center">0.01</td>
</tr>
<tr>
<td valign="top" align="left">4010</td>
<td valign="top" align="left">MAPK signaling pathway</td>
<td valign="top" align="left"><italic>LOC100912399, Mecom, Hspa1b, Efna1</italic></td>
<td valign="top" align="center">0.09</td>
</tr>
<tr>
<td valign="top" align="left">4713</td>
<td valign="top" align="left">Circadian entrainment</td>
<td valign="top" align="left"><italic>Ryr2, Itpr3, Grin2c</italic></td>
<td valign="top" align="center">0.02</td>
</tr>
<tr>
<td valign="top" align="left">4670</td>
<td valign="top" align="left">Leukocyte transendothelial migration</td>
<td valign="top" align="left"><italic>Cyba, Cd99, Cxcl12</italic></td>
<td valign="top" align="center">0.03</td>
</tr>
<tr>
<td valign="top" align="left">4724</td>
<td valign="top" align="left">Glutamatergic synapse</td>
<td valign="top" align="left"><italic>Grin2c, Itpr3, Slc38a3</italic></td>
<td valign="top" align="center">0.03</td>
</tr>
<tr>
<td valign="top" align="left">4611</td>
<td valign="top" align="left">Platelet activation</td>
<td valign="top" align="left"><italic>Gp1bb, Vwf, Itpr3</italic></td>
<td valign="top" align="center">0.04</td>
</tr>
<tr>
<td valign="top" align="left">4261</td>
<td valign="top" align="left">Adrenergic signaling in cardiomyocytes</td>
<td valign="top" align="left"><italic>Ryr2, Crem, Agtr2</italic></td>
<td valign="top" align="center">0.05</td>
</tr>
<tr>
<td valign="top" align="left">4921</td>
<td valign="top" align="left">Oxytocin signaling pathway</td>
<td valign="top" align="left"><italic>Ryr2, Rgs2, Itpr3</italic></td>
<td valign="top" align="center">0.06</td>
</tr>
<tr>
<td valign="top" align="left">4530</td>
<td valign="top" align="left">Tight junction</td>
<td valign="top" align="left"><italic>Scrib, LOC100912399, Cttn</italic></td>
<td valign="top" align="center">0.07</td>
</tr>
<tr>
<td valign="top" align="left">4360</td>
<td valign="top" align="left">Axon guidance</td>
<td valign="top" align="left"><italic>Efna1, Cxcl12, Robo3</italic></td>
<td valign="top" align="center">0.09</td>
</tr>
<tr>
<td valign="top" align="left">4024</td>
<td valign="top" align="left">cAMP signaling pathway</td>
<td valign="top" align="left"><italic>Grin2c, Lipe, Ryr2</italic></td>
<td valign="top" align="center">0.12</td>
</tr>
<tr>
<td valign="top" align="left">4060</td>
<td valign="top" align="left">Cytokine-cytokine receptor interaction</td>
<td valign="top" align="left"><italic>Cxcl12, Il17rc, Bmp3</italic></td>
<td valign="top" align="center">0.20</td>
</tr>
<tr>
<td valign="top" align="left">4080</td>
<td valign="top" align="left">Neuroactive ligand-receptor interaction</td>
<td valign="top" align="left"><italic>Htr2a, Agtr2, Grin2c</italic></td>
<td valign="top" align="center">0.33</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The term with the most statistical significance in the cluster is selected as the term representing the cluster, and then a rich subset of terms is selected and displayed as a network diagram to further determine the relationship between terms (<xref ref-type="fig" rid="F5">Figure 5C</xref>). Among them, calcium signaling pathway is related to the sensory perception of temperature stimulus.</p>
</sec>
<sec id="S3.SS5">
<title>Gene Set Enrichment Analysis</title>
<p>Gene Set Enrichment Analysis (GSEA) revealed the enrichment of EA responsive gene sets, identifying the regulation of autophagy as the top differentially modulated pathway, display NES of 1.564 (<italic>P</italic> = 0.00, FDR = 0.083). Consistently, GO analysis also revealed a change in TORC2 signaling (GO:0038203, <italic>P</italic> = 0.03) as shown in <xref ref-type="table" rid="T4">Table 4</xref>. MTOR (the mechanism target of rapamycin kinase) signal is the most well-known autophagy regulator and MTOR signal is composed of two multi-protein complexes with different structures and functions, mTORC1 and mTORC2.</p>
</sec>
<sec id="S3.SS6">
<title>Protein&#x2013;Protein Interactions Analysis</title>
<p>The PPI network was constructed by using DEGs (<xref ref-type="fig" rid="F6">Figure 6A</xref>). The three significant modules constructed from the PPI network of the DEGs including module 1 containing 5 nodes and 10 edges (MCODE score = 5), module 2 with 4 nodes and 6 edges (MCODE score = 4), and module 3 including 3 nodes and 3 edges (MCODE score = 3) (<xref ref-type="fig" rid="F6">Figure 6B</xref>). MCODE 1 combined is primarily associated with protein folding, MCODE 2 is mainly involved in rRNA pseudouridine synthesis while MCODE 3 is related to DNA replication.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption><p>Protein-protein interaction of the DEGs <bold>(A)</bold>, red represents up-regulation genes, blue represents down-regulation genes, and the thicker the line, the higher the connection score. Three functional modules within the PPI <bold>(B)</bold>. 10 hub genes <bold>(C)</bold> and their expression among samples <bold>(D)</bold>, correlation between 10 hub genes <bold>(E)</bold>.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g006.tif"/>
</fig>
<p>The 10 hub genes screened were <italic>Dnaja1</italic>, <italic>Dnajb1</italic>, <italic>Dnajb4</italic>, <italic>Hsph1</italic>, <italic>Chordc1, Nhp2</italic>, <italic>Utp11</italic>, <italic>Noc3l</italic>, <italic>Rsl1d1l1</italic> (<xref ref-type="fig" rid="F6">Figure 6C</xref>). According to the expression of hub genes in the samples (<xref ref-type="fig" rid="F6">Figure 6D</xref>), it is suggested that the expression of <italic>Hsph1</italic> and <italic>Dnaja1</italic> is the highest, followed by <italic>Dnajb1</italic>, <italic>Dnajb4</italic> and <italic>Utp11</italic>. However, <italic>Chordc1</italic>, <italic>Hspa12b</italic>, <italic>Nhp2</italic>, and <italic>Noc3l</italic> are less expressed in the samples, while <italic>Rsl1d1l1</italic> is very little. <xref ref-type="fig" rid="F6">Figure 6E</xref> shows the correlation of gene expression of hub genes.</p>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<title>Discussion</title>
<p>In the development of drugs for the treatment of CNS diseases, achieving adequate BBB transmission is a key challenge. The BBB is formed by brain microvascular endothelial cells that are sealed by tight junctions, which makes it an important obstacle to most brain therapies. Therefore, poor BBB permeability is the main reason for limiting the application of newly developed CNS therapeutic drugs in clinical practice.</p>
<p>Multiple ABC proteins are expressed in the BBB lumen and blood-oriented inner cortical membrane, which limits the permeability of a large number of toxins, including therapeutic agents (<xref ref-type="bibr" rid="B24">Miller, 2015</xref>). ABC transporters are ATP-driven efflux pumps of exogenous and endogenous metabolites, and their high expression in BBB leads to drug resistance in CNS. The decreased expression and/or functional activity of ABC-BBB transporters in Alzheimer&#x2019;s disease (AD) and Parkinson&#x2019;s disease (PD) patients (<xref ref-type="bibr" rid="B44">Zlokovic, 2011</xref>), and in the brain of AD animal models proved that the decreased expression and/or functional activity of ABC-BBB transporters led to the accumulation of amyloid beta peptide (A&#x03B2;) (<xref ref-type="bibr" rid="B10">Cirrito et al., 2005</xref>). However, the clinical potential of ABC transporters for disease management and drug delivery improvement remains unclear. EA can reduce the expression of <italic>Abcc6</italic> (<xref ref-type="table" rid="T3">Table 3</xref>), suggesting that further studies can explore whether electroacupuncture can mediate the clearance of harmful sediments such as A&#x03B2; and the delivery of CNS therapeutic drugs by affecting the expression of <italic>Abcc6</italic>. Furthermore, EA intervention can change the expression of multiple SLC, including <italic>Slc7a1</italic>, <italic>Slc9a3r2</italic>, <italic>Slc38a3</italic>, <italic>Slco2b1</italic>, and <italic>Slc25a25</italic> (<xref ref-type="table" rid="T3">Table 3</xref>). The non-invasive delivery of therapeutic drugs through BBB may use endogenous processes, such as adsorption-mediated cell transport, CMT and RMT (<xref ref-type="bibr" rid="B32">Terstappen et al., 2021</xref>). Among them, the CMT system is expressed by genes in the SLC transporter gene family, including more than 300 transporter genes encoding membrane-binding proteins, which are helpful for the trans-membrane transport of multiple substrates (<xref ref-type="bibr" rid="B21">Lin et al., 2015</xref>). In BBB, SLC proteins promote cross-cellular transport of multiple molecules, including carbohydrates, amino acids, monocarboxylic acids, hormones, fatty acids, nucleotides, organic anions, amines, choline and vitamins (<xref ref-type="bibr" rid="B43">Zlokovic, 2008</xref>; <xref ref-type="bibr" rid="B11">Daneman and Prat, 2015</xref>; <xref ref-type="bibr" rid="B26">Pardridge, 2015</xref>). It is worth noting that <italic>Cavin2</italic> is also a gene related to transcytosis, which can be significantly down regulated by electroacupuncture; Other genes related to endocytosis include <italic>Cttn</italic>, <italic>Cdc45</italic>, <italic>Hspa12b</italic>, <italic>Cdcc117</italic>, <italic>Cdc42ep3</italic>, Cblb and <italic>Hspa1b</italic> (<xref ref-type="table" rid="T3">Table 3</xref>). The role of these genes in maintaining the integrity of BBB should be further studied.</p>
<p>Basement membrane (BM) is one of the important components of BBB because endothelial cells are connected to the BM through integrins. In addition, integrins also interact with extracellular matrix proteins such as laminin, collagen and perlecan, and mediate signal transduction by activating ECM ligands, growth factors and growth factor receptors (<xref ref-type="bibr" rid="B4">Baeten and Akassoglou, 2011</xref>). Our GO analysis showed that the mechanism of EA intervention involved integrin binding (GO:0005178, <italic>P</italic> = 0.0426) and collagen-containing extracellular matrix (GO:0062023, <italic>P</italic> = 0.0071) (<xref ref-type="table" rid="T4">Table 4</xref>). Mice lacking &#x03B2;1- integrin in endothelial cells produce abnormal cadherin signaling, claudin-5 deficiency and immature BBB (<xref ref-type="bibr" rid="B35">Yamamoto et al., 2015</xref>). Similarly, mice lacking laminin secreted by astrocytes showed BBB breakdown (<xref ref-type="bibr" rid="B37">Yao et al., 2014</xref>).</p>
<p>It is well known that the TJ in the BBB is a dynamic structure. TJ proteins undergo changes in expression, subcellular localization, post-translational modifications and protein-protein interactions under physiological and pathophysiological conditions (<xref ref-type="bibr" rid="B17">Huber et al., 2001</xref>). Interactions between TJ proteins and intracellular signaling pathways have begun to be elucidated (<xref ref-type="bibr" rid="B6">Bazzoni and Dejana, 2004</xref>). The cultured epithelial cells in the medium lacking Ca<sup>2 +</sup> showed the loss of cell membrane ZO-1, ZO-2 and thromboxane, and the increase of related paracellular permeability (<xref ref-type="bibr" rid="B19">Klingler et al., 2000</xref>). Interestingly, although artificially increasing intracellular [Ca<sup>2 +</sup>] does not change BBB TJ (<xref ref-type="bibr" rid="B8">Brown et al., 2004</xref>), the treatment with Ca<sup>2 +</sup> channel blocker SKF 96365 can improve the permeability increase effect of hypoxia/hyperglycemia on BBB through a mechanism that seems to involve the recruitment of occludin into the cell membrane (<xref ref-type="bibr" rid="B7">Brown and Davis, 2005</xref>). The results of KEGG analysis showed that the most significant enrichment pathway was the calcium signaling pathway, and the differential genes associated with this pathway are <italic>Plcd3, Ryr2, Itpr3, Htr2a, Grin2c, Sphk2</italic> (<xref ref-type="fig" rid="F6">Figure 6</xref> and <xref ref-type="table" rid="T5">Table 5</xref>). <italic>Itpr3</italic> is an inositol 1,4,5-triphosphate receptor, type 3 (InsP3R3), encoding a receptor for inositol 1,4,5-trisphosphate, a second messenger that mediates the release of intracellular calcium. Although EA intervention did not significantly change the mRNA level of TJ proteins in this experiment, the previous study of our research group found that EA intervention could affect the expression level of ZO-1 and occluding (<xref ref-type="bibr" rid="B40">Zhang et al., 2020</xref>). Therefore, it is not excluded that EA regulates the subcellular localization and post-translational modification of TJ proteins through the calcium signaling pathway. The next step can study the effect of EA intervention on the protein expression level of TJ proteins and the possible role of calcium signaling in it.</p>
<p>Autophagy persists at the basic level of cells and is the mechanism of organelle or protein degradation and recycling (<xref ref-type="bibr" rid="B20">Lakhan et al., 2013</xref>). It is stimulated by a variety of pathological processes and participates in the elimination of organelle and protein decomposition, thereby avoiding excessive damage and cell dysfunction in various organs and cells (<xref ref-type="bibr" rid="B12">De Meyer et al., 2015</xref>). <italic>In vivo</italic> experiments showed that rapamycin (an autophagy inducer) led to the decrease of ZO-1 and the increase of BBB permeability, and confirmed the role of autophagy in regulating paracellular permeability (<xref ref-type="bibr" rid="B39">Zhang et al., 2018</xref>). However, studies have reported that CLDN5 (claudin 5) is abnormally accumulated in the cytoplasm of BMECs in stroke patients, accompanied by autophagy activation. Studies on zebrafish <italic>in vivo</italic> and <italic>in vitro</italic> cells showed that the BBB decomposition was caused by the redistribution of CLDN5 mediated by CAV1 (caveolin 1) into the cytoplasm under hypoxia. At the same time, autophagy is activated, which is mainly helpful degradation of CAV1 and CLDN5 in BMECs cytoplasm. These data suggest that autophagy plays a role in BBB integrity by regulating the redistribution of CLDN5 (<xref ref-type="bibr" rid="B36">Yang et al., 2021</xref>). Also, autophagy involves the dynamic disposition of intracellular and extracellular components, theoretically affecting BBB penetration. For example, autophagy-mediated lysosomal degradation pathway is involved in the degradation of PPR/<italic>p</italic> TRAIL in brain capillary endothelial cells and prevents them from penetrating BBB. The pre-inhibition of BBB autophagy by wortmannin loaded liposomes (Wtmn-Lip) can increase the accumulation of non-viral gene vector PPR in the brain without damaging the BBB tight junction (<xref ref-type="bibr" rid="B33">Wang et al., 2020</xref>). In the study of APP/PS1 mice, LSP (Lychee seed polyphenol) induced autophagy in bEnd.3 cells through AMPK/mTOR/ULK1 pathway and increased TJs expression (<xref ref-type="bibr" rid="B34">Xiong et al., 2021</xref>). The above data show that autophagy can regulate BBB permeability and material transport, but there is still no direct evidence to explain the exact relationship between the two, which may be related to diseases and the types of models used in the study. Our GSEA results showed that EA at specific stimulation mode could activate autophagy (<xref ref-type="fig" rid="F7">Figure 7</xref>) and inhibit TORC2 signal (GO: 0038203, <italic>P</italic> = 0.0311), and the genes involved were significantly down-regulated Prr5l (<xref ref-type="table" rid="T4">Table 4</xref>). Therefore, we believe that autophagy is closely related to the translocation and degradation of TJ proteins in BBB, which is also an idea to study the mechanism of EA opening BBB. However, the protein expression and mRNA expression are not completely consistent, because the sample of this experiment is the total RNA of rat cortical tissue, not a single cell, which is also one of the limitations of this experiment. Future studies should first confirm whether the autophagy pathway activated by EA is involved in the translocation and degradation of BBB TJ protein. Secondly, whether the degree of autophagy activation is enough to cause cell damage. Finally, is the autophagy activation after EA a direct effect or is it because EA changes the cellular distribution of TJ proteins (such as translocation from cell membrane to cytoplasm), thus promoting autophagy activation.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption><p>Gene set enrichment analysis for &#x201C;regulation of autophagy.&#x201D; Screening with NES &#x003E; 1, <italic>P</italic> &#x003C; 0.05 and FDR &#x003C; 0.25.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fnins-16-834683-g007.tif"/>
</fig>
<p>In addition, EA may have an effect on the clathrin-coated pit (GO:0005905, <italic>P</italic> = 0.0424), involving genes including <italic>Cttn</italic> and <italic>Sh3bp4</italic> (<xref ref-type="table" rid="T4">Table 4</xref>). Clathrin-coated pit is a part of the endomembrane system in the form of an invagination of a membrane upon which a clathrin coat forms, and that can be converted by vesicle budding into a clathrin-coated vesicle. Coated pits form on the plasma membrane, where they are involved in receptor-mediated selective transport of many proteins and other macromolecules across the cell membrane. Experimental studies on various animal models have shown that A&#x03B2; is mainly cleared through the BBB across blood vessels (70&#x2013;85%), while a small part is cleared through the ISF flow (<xref ref-type="bibr" rid="B3">Bading et al., 2002</xref>; <xref ref-type="bibr" rid="B13">Deane et al., 2004</xref>; <xref ref-type="bibr" rid="B31">Tarasoff-Conway et al., 2015</xref>) and the molecular mechanism of A&#x03B2; clearance through BBB has been clarified in more detail recently. In short, A&#x03B2; produced in the brain binds to LRP1 on the albumin side of BBB, resulting in its rapid internalization into endothelial cells and clearance through blood. Phosphatidylinositol-binding clathrin assembly protein (PICALM) is crucial for the internalization of LRP1 &#x2013; A&#x03B2; complex mediated by endothelial clathrin/PICALM, and guides the transport of endocytosis vesicles containing A&#x03B2; in endothelial cells by sequential fusion with Rab5 positive early endosomes and Rab11 positive sorting endosomes used for the endocytosis of BBB luminal (<xref ref-type="bibr" rid="B41">Zhao et al., 2015</xref>).</p>
<p>However, as we mentioned above, our sample is a mixture of all cells in the cortex, which can not directly prove that these genes and proteins are expressed in which kind of cells of the blood-brain barrier. Our research mainly focuses on the role of endothelial cells, and the blood-brain barrier is not a single entity. It is composed of astrocytes, pericytes and endothelial cells, which allow the endothelial cells of the central nervous system to strictly regulate the movement of ions, molecules and cells between blood and neural tissues. Astrocytes are the most abundant cells in the brain, which are very important in regulating brain and BBB function. A recent breakthrough study (<xref ref-type="bibr" rid="B18">Kitchen et al., 2020</xref>) proved that targeting astrocytes at the blood-brain barrier (BBB) and spinal cord barrier (BSCB) is a feasible option for the treatment of TBI and spinal cord injury. This study showed that pharmacological inhibition of these signal events could prevent the development of central nervous system edema and promote the functional recovery of injured rats, which was confirmed by <xref ref-type="bibr" rid="B30">Sylvain et al. (2021)</xref>. Recently, the molecular mechanism of clearing A&#x03B2; through BBB has been elucidated in more detail. Lymphatic pathway is a waste removal system, which plays an important role in dementia. It uses the unique system of perivascular channels to promote the effective removal of soluble proteins and metabolites from the central nervous system. AQP4 is highly expressed in astrocytes. The enrichment of perivascular AQP4 at BSCB/BBB shows that it plays a role in lymphatic function. The development of new drugs for this system will have great therapeutic potential (<xref ref-type="bibr" rid="B27">Salman et al., 2021</xref>). In addition, pericytes also play an important role in the maintenance of BBB function. In addition to communicating with astrocytes, pericytes also support the maintenance of BBB in postnatal brain. Pericyte deficient mutant mice showed increased permeability of BBB to low molecular weight and high molecular weight tracers, and pericytes can directly regulate the expression of transporters (<xref ref-type="bibr" rid="B1">Arvanitis et al., 2020</xref>).</p>
<p>Achieving adequate BBB transmission is a key challenge in the development of drugs for the treatment of central nervous system (CNS) diseases. However, most traditional BBB opening strategies are difficult to be applied in clinic because of their extensive and non-specific regulation of BBB, resulting in damage to normal brain tissue. Therefore, non-invasive methods that selectively and effectively cross the BBB, such as receptor-mediated cell transport and the use of neurophilic viruses, nanoparticles and exosomes, have become popular in recent years (<xref ref-type="bibr" rid="B32">Terstappen et al., 2021</xref>). For example, in the process of glioma treatment, the poor permeability of BBB is the main reason that restricts the application of newly developed treatment methods for glioma in clinic. Researchers have proposed a variety of strategies to target and effectively regulate the BBB of tumor, so as to obtain more effective glioma treatment (<xref ref-type="bibr" rid="B23">Luo and Shusta, 2020</xref>).</p>
</sec>
<sec id="S5" sec-type="conclusion">
<title>Conclusion</title>
<p>Our research will provide valuable resources for discovering and studying new targets and signaling pathways that mediate crosstalk between different cell types in NVU. This may lead to the development of new transgenic animal models &#x2013; BBB, NVU, and pluripotent stem cell models of different neurological diseases, as a valuable way to discover and test new drug delivery methods.</p>
</sec>
<sec id="S6" sec-type="data-availability">
<title>Data Availability Statement</title>
<p>The original contributions presented in the study are publicly available. This data can be found here: <ext-link ext-link-type="uri" xlink:href="https://www.ncbi.nlm.nih.gov/geo/">https://www.ncbi.nlm.nih.gov/geo/</ext-link>, GSE192885.</p>
</sec>
<sec id="S7">
<title>Ethics Statement</title>
<p>This study was approved by the Animal Protection and Use Committee of Zhejiang University of Traditional Chinese Medicine.</p>
</sec>
<sec id="S8">
<title>Author Contributions</title>
<p>CM was responsible for drafting and writing the manuscript. LG and HW analyzed and sorted out the data. LR and YL provided help for the experiment. SZ and PG contributed to make important modifications to the manuscript. XL approved the final manuscript to be published. All authors contributed to the article and approved the submitted version.</p>
</sec>
<sec id="conf1" sec-type="COI-statement">
<title>Conflict of Interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="pudiscl1" sec-type="disclaimer">
<title>Publisher&#x2019;s Note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<sec id="S9" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by the National Natural Science Foundation of China (No. 82174502) and the Key Laboratory of Acupuncture and Neurology of Zhejiang Province.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Arvanitis</surname> <given-names>C. D.</given-names></name> <name><surname>Ferraro</surname> <given-names>G. B.</given-names></name> <name><surname>Jain</surname> <given-names>R. K.</given-names></name></person-group> (<year>2020</year>). <article-title>The blood-brain barrier and blood-tumour barrier in brain tumours and metastases.</article-title> <source><italic>Nat. Rev. Cancer</italic></source> <volume>20</volume> <fpage>26</fpage>&#x2013;<lpage>41</lpage>. <pub-id pub-id-type="doi">10.1038/s41568-019-0205-x</pub-id> <pub-id pub-id-type="pmid">31601988</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bader</surname> <given-names>G. D.</given-names></name> <name><surname>Hogue</surname> <given-names>C. W.</given-names></name></person-group> (<year>2003</year>). <article-title>An automated method for finding molecular complexes in large protein interaction networks.</article-title> <source><italic>BMC Bioinform.</italic></source> <volume>4</volume>:<issue>2</issue>. <pub-id pub-id-type="doi">10.1186/1471-2105-4-2</pub-id> <pub-id pub-id-type="pmid">12525261</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bading</surname> <given-names>J. R.</given-names></name> <name><surname>Yamada</surname> <given-names>S.</given-names></name> <name><surname>Mackic</surname> <given-names>J. B.</given-names></name> <name><surname>Kirkman</surname> <given-names>L.</given-names></name> <name><surname>Miller</surname> <given-names>C.</given-names></name> <name><surname>Calero</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2002</year>). <article-title>Brain clearance of Alzheimer&#x2019;s amyloid-beta40 in the squirrel monkey: a SPECT study in a primate model of cerebral amyloid angiopathy.</article-title> <source><italic>J. Drug Target</italic></source> <volume>10</volume> <fpage>359</fpage>&#x2013;<lpage>368</lpage>. <pub-id pub-id-type="doi">10.1080/10611860290031831</pub-id> <pub-id pub-id-type="pmid">12164385</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Baeten</surname> <given-names>K. M.</given-names></name> <name><surname>Akassoglou</surname> <given-names>K.</given-names></name></person-group> (<year>2011</year>). <article-title>Extracellular matrix and matrix receptors in blood-brain barrier formation and stroke.</article-title> <source><italic>Dev. Neurobiol.</italic></source> <volume>71</volume> <fpage>1018</fpage>&#x2013;<lpage>1039</lpage>. <pub-id pub-id-type="doi">10.1002/dneu.20954</pub-id> <pub-id pub-id-type="pmid">21780303</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Banks</surname> <given-names>W. A.</given-names></name></person-group> (<year>2016</year>). <article-title>From blood-brain barrier to blood-brain interface: new opportunities for CNS drug delivery.</article-title> <source><italic>Nat. Rev. Drug Discov.</italic></source> <volume>15</volume> <fpage>275</fpage>&#x2013;<lpage>292</lpage>. <pub-id pub-id-type="doi">10.1038/nrd.2015.21</pub-id> <pub-id pub-id-type="pmid">26794270</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bazzoni</surname> <given-names>G.</given-names></name> <name><surname>Dejana</surname> <given-names>E.</given-names></name></person-group> (<year>2004</year>). <article-title>Endothelial cell-to-cell junctions: molecular organization and role in vascular homeostasis.</article-title> <source><italic>Physiol. Rev.</italic></source> <volume>84</volume> <fpage>869</fpage>&#x2013;<lpage>901</lpage>. <pub-id pub-id-type="doi">10.1152/physrev.00035.2003</pub-id> <pub-id pub-id-type="pmid">15269339</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname> <given-names>R. C.</given-names></name> <name><surname>Davis</surname> <given-names>T. P.</given-names></name></person-group> (<year>2005</year>). <article-title>Hypoxia/aglycemia alters expression of occludin and actin in brain endothelial cells.</article-title> <source><italic>Biochem. Biophys. Res. Commun.</italic></source> <volume>327</volume> <fpage>1114</fpage>&#x2013;<lpage>1123</lpage>. <pub-id pub-id-type="doi">10.1016/j.bbrc.2004.12.123</pub-id> <pub-id pub-id-type="pmid">15652512</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname> <given-names>R. C.</given-names></name> <name><surname>Mark</surname> <given-names>K. S.</given-names></name> <name><surname>Egleton</surname> <given-names>R. D.</given-names></name> <name><surname>Davis</surname> <given-names>T. P.</given-names></name></person-group> (<year>2004</year>). <article-title>Protection against hypoxia-induced blood-brain barrier disruption: changes in intracellular calcium.</article-title> <source><italic>Am. J. Physiol. Cell Physiol.</italic></source> <volume>286</volume> <fpage>C1045</fpage>&#x2013;<lpage>C1052</lpage>. <pub-id pub-id-type="doi">10.1152/ajpcell.00360.2003</pub-id> <pub-id pub-id-type="pmid">15075203</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chin</surname> <given-names>C. H.</given-names></name> <name><surname>Chen</surname> <given-names>S. H.</given-names></name> <name><surname>Wu</surname> <given-names>H. H.</given-names></name> <name><surname>Ho</surname> <given-names>C. W.</given-names></name> <name><surname>Ko</surname> <given-names>M. T.</given-names></name> <name><surname>Lin</surname> <given-names>C. Y.</given-names></name></person-group> (<year>2014</year>). <article-title>cytoHubba: identifying hub objects and sub-networks from complex interactome.</article-title> <source><italic>BMC Syst. Biol.</italic></source> <volume>8(Suppl. 4)</volume>:<issue>S11</issue>. <pub-id pub-id-type="doi">10.1186/1752-0509-8-S4-S11</pub-id> <pub-id pub-id-type="pmid">25521941</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cirrito</surname> <given-names>J. R.</given-names></name> <name><surname>Deane</surname> <given-names>R.</given-names></name> <name><surname>Fagan</surname> <given-names>A. M.</given-names></name> <name><surname>Spinner</surname> <given-names>M. L.</given-names></name> <name><surname>Parsadanian</surname> <given-names>M.</given-names></name> <name><surname>Finn</surname> <given-names>M. B.</given-names></name><etal/></person-group> (<year>2005</year>). <article-title>P-glycoprotein deficiency at the blood-brain barrier increases amyloid-beta deposition in an Alzheimer disease mouse model.</article-title> <source><italic>J. Clin. Invest.</italic></source> <volume>115</volume> <fpage>3285</fpage>&#x2013;<lpage>3290</lpage>. <pub-id pub-id-type="doi">10.1172/JCI25247</pub-id> <pub-id pub-id-type="pmid">16239972</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Daneman</surname> <given-names>R.</given-names></name> <name><surname>Prat</surname> <given-names>A.</given-names></name></person-group> (<year>2015</year>). <article-title>The blood-brain barrier.</article-title> <source><italic>Cold Spring Harb. Perspect. Biol.</italic></source> <volume>7</volume>:<issue>a020412</issue>.</citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Meyer</surname> <given-names>G. R.</given-names></name> <name><surname>Grootaert</surname> <given-names>M. O.</given-names></name> <name><surname>Michiels</surname> <given-names>C. F.</given-names></name> <name><surname>Kurdi</surname> <given-names>A.</given-names></name> <name><surname>Schrijvers</surname> <given-names>D. M.</given-names></name> <name><surname>Martinet</surname> <given-names>W.</given-names></name></person-group> (<year>2015</year>). <article-title>Autophagy in vascular disease.</article-title> <source><italic>Circ. Res.</italic></source> <volume>116</volume> <fpage>468</fpage>&#x2013;<lpage>479</lpage>.</citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deane</surname> <given-names>R.</given-names></name> <name><surname>Wu</surname> <given-names>Z.</given-names></name> <name><surname>Sagare</surname> <given-names>A.</given-names></name> <name><surname>Davis</surname> <given-names>J.</given-names></name> <name><surname>Du Yan</surname> <given-names>S.</given-names></name> <name><surname>Hamm</surname> <given-names>K.</given-names></name><etal/></person-group> (<year>2004</year>). <article-title>LRP/amyloid beta-peptide interaction mediates differential brain efflux of Abeta isoforms.</article-title> <source><italic>Neuron</italic></source> <volume>43</volume> <fpage>333</fpage>&#x2013;<lpage>344</lpage>. <pub-id pub-id-type="doi">10.1016/j.neuron.2004.07.017</pub-id> <pub-id pub-id-type="pmid">15294142</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Engelhardt</surname> <given-names>B.</given-names></name></person-group> (<year>2008</year>). <article-title>Immune cell entry into the central nervous system: involvement of adhesion molecules and chemokines.</article-title> <source><italic>J. Neurol. Sci.</italic></source> <volume>274</volume> <fpage>23</fpage>&#x2013;<lpage>26</lpage>. <pub-id pub-id-type="doi">10.1016/j.jns.2008.05.019</pub-id> <pub-id pub-id-type="pmid">18573502</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Furuse</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title>Molecular basis of the core structure of tight junctions.</article-title> <source><italic>Cold Spring Harb. Perspect. Biol.</italic></source> <volume>2</volume>:<issue>a002907</issue>. <pub-id pub-id-type="doi">10.1101/cshperspect.a002907</pub-id> <pub-id pub-id-type="pmid">20182608</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hallmann</surname> <given-names>R.</given-names></name> <name><surname>Mayer</surname> <given-names>D. N.</given-names></name> <name><surname>Berg</surname> <given-names>E. L.</given-names></name> <name><surname>Broermann</surname> <given-names>R.</given-names></name> <name><surname>Butcher</surname> <given-names>E. C.</given-names></name></person-group> (<year>1995</year>). <article-title>Novel mouse endothelial cell surface marker is suppressed during differentiation of the blood brain barrier.</article-title> <source><italic>Dev. Dyn.</italic></source> <volume>202</volume> <fpage>325</fpage>&#x2013;<lpage>332</lpage>. <pub-id pub-id-type="doi">10.1002/aja.1002020402</pub-id> <pub-id pub-id-type="pmid">7626790</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huber</surname> <given-names>J. D.</given-names></name> <name><surname>Witt</surname> <given-names>K. A.</given-names></name> <name><surname>Hom</surname> <given-names>S.</given-names></name> <name><surname>Egleton</surname> <given-names>R. D.</given-names></name> <name><surname>Mark</surname> <given-names>K. S.</given-names></name> <name><surname>Davis</surname> <given-names>T. P.</given-names></name></person-group> (<year>2001</year>). <article-title>Inflammatory pain alters blood-brain barrier permeability and tight junctional protein expression.</article-title> <source><italic>Am. J. Physiol. Heart Circ. Physiol.</italic></source> <volume>280</volume> <fpage>H1241</fpage>&#x2013;<lpage>H1248</lpage>. <pub-id pub-id-type="doi">10.1152/ajpheart.2001.280.3.H1241</pub-id> <pub-id pub-id-type="pmid">11179069</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kitchen</surname> <given-names>P.</given-names></name> <name><surname>Salman</surname> <given-names>M. M.</given-names></name> <name><surname>Halsey</surname> <given-names>A. M.</given-names></name> <name><surname>Clarke-Bland</surname> <given-names>C.</given-names></name> <name><surname>Macdonald</surname> <given-names>J. A.</given-names></name> <name><surname>Ishida</surname> <given-names>H.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Targeting Aquaporin-4 subcellular localization to treat central nervous system edema.</article-title> <source><italic>Cell</italic></source> <volume>181</volume>:<issue>e719</issue>. <pub-id pub-id-type="doi">10.1016/j.cell.2020.03.037</pub-id> <pub-id pub-id-type="pmid">32413299</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Klingler</surname> <given-names>C.</given-names></name> <name><surname>Kniesel</surname> <given-names>U.</given-names></name> <name><surname>Bamforth</surname> <given-names>S. D.</given-names></name> <name><surname>Wolburg</surname> <given-names>H.</given-names></name> <name><surname>Engelhardt</surname> <given-names>B.</given-names></name> <name><surname>Risau</surname> <given-names>W.</given-names></name></person-group> (<year>2000</year>). <article-title>Disruption of epithelial tight junctions is prevented by cyclic nucleotide-dependent protein kinase inhibitors.</article-title> <source><italic>Histochem. Cell Biol.</italic></source> <volume>113</volume> <fpage>349</fpage>&#x2013;<lpage>361</lpage>. <pub-id pub-id-type="doi">10.1007/s004180000143</pub-id> <pub-id pub-id-type="pmid">10883394</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lakhan</surname> <given-names>S. E.</given-names></name> <name><surname>Kirchgessner</surname> <given-names>A.</given-names></name> <name><surname>Tepper</surname> <given-names>D.</given-names></name> <name><surname>Leonard</surname> <given-names>A.</given-names></name></person-group> (<year>2013</year>). <article-title>Matrix metalloproteinases and blood-brain barrier disruption in acute ischemic stroke.</article-title> <source><italic>Front. Neurol.</italic></source> <volume>4</volume>:<issue>32</issue>. <pub-id pub-id-type="doi">10.3389/fneur.2013.00032</pub-id> <pub-id pub-id-type="pmid">23565108</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lin</surname> <given-names>L.</given-names></name> <name><surname>Yee</surname> <given-names>S. W.</given-names></name> <name><surname>Kim</surname> <given-names>R. B.</given-names></name> <name><surname>Giacomini</surname> <given-names>K. M.</given-names></name></person-group> (<year>2015</year>). <article-title>SLC transporters as therapeutic targets: emerging opportunities.</article-title> <source><italic>Nat. Rev. Drug Discov.</italic></source> <volume>14</volume> <fpage>543</fpage>&#x2013;<lpage>560</lpage>. <pub-id pub-id-type="doi">10.1038/nrd4626</pub-id> <pub-id pub-id-type="pmid">26111766</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Livak</surname> <given-names>K. J.</given-names></name> <name><surname>Schmittgen</surname> <given-names>T. D.</given-names></name></person-group> (<year>2001</year>). <article-title>Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method.</article-title> <source><italic>Methods</italic></source> <volume>25</volume> <fpage>402</fpage>&#x2013;<lpage>408</lpage>. <pub-id pub-id-type="doi">10.1006/meth.2001.1262</pub-id> <pub-id pub-id-type="pmid">11846609</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luo</surname> <given-names>H.</given-names></name> <name><surname>Shusta</surname> <given-names>E. V.</given-names></name></person-group> (<year>2020</year>). <article-title>Blood-brain barrier modulation to improve glioma drug delivery.</article-title> <source><italic>Pharmaceutics</italic></source> <volume>12</volume>:<issue>1085</issue>. <pub-id pub-id-type="doi">10.3390/pharmaceutics12111085</pub-id> <pub-id pub-id-type="pmid">33198244</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Miller</surname> <given-names>D. S.</given-names></name></person-group> (<year>2015</year>). <article-title>Regulation of ABC transporters blood-brain barrier: the good, the bad, and the ugly.</article-title> <source><italic>Adv. Cancer Res.</italic></source> <volume>125</volume> <fpage>43</fpage>&#x2013;<lpage>70</lpage>. <pub-id pub-id-type="doi">10.1016/bs.acr.2014.10.002</pub-id> <pub-id pub-id-type="pmid">25640266</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Murtagh</surname> <given-names>F.</given-names></name> <name><surname>Legendre</surname> <given-names>P.</given-names></name></person-group> (<year>2014</year>). <article-title>Ward&#x2019;s hierarchical agglomerative clustering method: which algorithms implement Ward&#x2019;s criterion?</article-title> <source><italic>J. Classification</italic></source> <volume>31</volume> <fpage>274</fpage>&#x2013;<lpage>295</lpage>.</citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pardridge</surname> <given-names>W. M.</given-names></name></person-group> (<year>2015</year>). <article-title>Blood-brain barrier endogenous transporters as therapeutic targets: a new model for small molecule CNS drug discovery.</article-title> <source><italic>Expert Opin. Ther. Targets</italic></source> <volume>19</volume> <fpage>1059</fpage>&#x2013;<lpage>1072</lpage>. <pub-id pub-id-type="doi">10.1517/14728222.2015.1042364</pub-id> <pub-id pub-id-type="pmid">25936389</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Salman</surname> <given-names>M. M.</given-names></name> <name><surname>Kitchen</surname> <given-names>P.</given-names></name> <name><surname>Halsey</surname> <given-names>A.</given-names></name> <name><surname>Wang</surname> <given-names>M. X.</given-names></name> <name><surname>Tornroth-Horsefield</surname> <given-names>S.</given-names></name> <name><surname>Conner</surname> <given-names>A. C.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>Emerging roles for dynamic aquaporin-4 subcellular relocalization in CNS water homeostasis.</article-title> <source><italic>Brain [Online ahead of print]</italic></source> awab311. <pub-id pub-id-type="doi">10.1093/brain/awab311</pub-id> <pub-id pub-id-type="pmid">34499128</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shannon</surname> <given-names>P.</given-names></name> <name><surname>Markiel</surname> <given-names>A.</given-names></name> <name><surname>Ozier</surname> <given-names>O.</given-names></name> <name><surname>Baliga</surname> <given-names>N. S.</given-names></name> <name><surname>Wang</surname> <given-names>J. T.</given-names></name> <name><surname>Ramage</surname> <given-names>D.</given-names></name><etal/></person-group> (<year>2003</year>). <article-title>Cytoscape: a software environment for integrated models of biomolecular interaction networks.</article-title> <source><italic>Genome Res.</italic></source> <volume>13</volume> <fpage>2498</fpage>&#x2013;<lpage>2504</lpage>. <pub-id pub-id-type="doi">10.1101/gr.1239303</pub-id> <pub-id pub-id-type="pmid">14597658</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Subramanian</surname> <given-names>A.</given-names></name> <name><surname>Tamayo</surname> <given-names>P.</given-names></name> <name><surname>Mootha</surname> <given-names>V. K.</given-names></name> <name><surname>Mukherjee</surname> <given-names>S.</given-names></name> <name><surname>Ebert</surname> <given-names>B. L.</given-names></name> <name><surname>Gillette</surname> <given-names>M. A.</given-names></name><etal/></person-group> (<year>2005</year>). <article-title>Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>102</volume> <fpage>15545</fpage>&#x2013;<lpage>15550</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0506580102</pub-id> <pub-id pub-id-type="pmid">16199517</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sylvain</surname> <given-names>N. J.</given-names></name> <name><surname>Salman</surname> <given-names>M. M.</given-names></name> <name><surname>Pushie</surname> <given-names>M. J.</given-names></name> <name><surname>Hou</surname> <given-names>H.</given-names></name> <name><surname>Meher</surname> <given-names>V.</given-names></name> <name><surname>Herlo</surname> <given-names>R.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>The effects of trifluoperazine on brain edema, aquaporin-4 expression and metabolic markers during the acute phase of stroke using photothrombotic mouse model.</article-title> <source><italic>Biochim. Biophys. Acta Biomembr.</italic></source> <volume>1863</volume>:<issue>183573</issue>. <pub-id pub-id-type="doi">10.1016/j.bbamem.2021.183573</pub-id> <pub-id pub-id-type="pmid">33561476</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tarasoff-Conway</surname> <given-names>J. M.</given-names></name> <name><surname>Carare</surname> <given-names>R. O.</given-names></name> <name><surname>Osorio</surname> <given-names>R. S.</given-names></name> <name><surname>Glodzik</surname> <given-names>L.</given-names></name> <name><surname>Butler</surname> <given-names>T.</given-names></name> <name><surname>Fieremans</surname> <given-names>E.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Clearance systems in the brain-implications for Alzheimer disease.</article-title> <source><italic>Nat. Rev. Neurol.</italic></source> <volume>11</volume> <fpage>457</fpage>&#x2013;<lpage>470</lpage>. <pub-id pub-id-type="doi">10.1038/nrneurol.2015.119</pub-id> <pub-id pub-id-type="pmid">26195256</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Terstappen</surname> <given-names>G. C.</given-names></name> <name><surname>Meyer</surname> <given-names>A. H.</given-names></name> <name><surname>Bell</surname> <given-names>R. D.</given-names></name> <name><surname>Zhang</surname> <given-names>W.</given-names></name></person-group> (<year>2021</year>). <article-title>Strategies for delivering therapeutics across the blood-brain barrier.</article-title> <source><italic>Nat. Rev. Drug Discov.</italic></source> <volume>20</volume> <fpage>362</fpage>&#x2013;<lpage>383</lpage>. <pub-id pub-id-type="doi">10.1038/s41573-021-00139-y</pub-id> <pub-id pub-id-type="pmid">33649582</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>X.</given-names></name> <name><surname>Yin</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>M.</given-names></name> <name><surname>Rao</surname> <given-names>J.</given-names></name> <name><surname>Wan</surname> <given-names>D.</given-names></name> <name><surname>Qiu</surname> <given-names>Y.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Autophagy inhibition changes the disposition of non-viral gene carriers during blood-brain barrier penetration and enhances TRAIL-induced apoptosis in brain metastatic tumor.</article-title> <source><italic>J. Control Release</italic></source> <volume>321</volume> <fpage>497</fpage>&#x2013;<lpage>508</lpage>. <pub-id pub-id-type="doi">10.1016/j.jconrel.2020.02.042</pub-id> <pub-id pub-id-type="pmid">32112851</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiong</surname> <given-names>R.</given-names></name> <name><surname>Zhou</surname> <given-names>X. G.</given-names></name> <name><surname>Tang</surname> <given-names>Y.</given-names></name> <name><surname>Wu</surname> <given-names>J. M.</given-names></name> <name><surname>Sun</surname> <given-names>Y. S.</given-names></name> <name><surname>Teng</surname> <given-names>J. F.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>Lychee seed polyphenol protects the blood-brain barrier through inhibiting Abeta(25-35)-induced NLRP3 inflammasome activation via the AMPK/mTOR/ULK1-mediated autophagy in bEnd.3 cells and APP/PS1 mice.</article-title> <source><italic>Phytother. Res.</italic></source> <volume>35</volume> <fpage>954</fpage>&#x2013;<lpage>973</lpage>. <pub-id pub-id-type="doi">10.1002/ptr.6849</pub-id> <pub-id pub-id-type="pmid">32893437</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yamamoto</surname> <given-names>H.</given-names></name> <name><surname>Ehling</surname> <given-names>M.</given-names></name> <name><surname>Kato</surname> <given-names>K.</given-names></name> <name><surname>Kanai</surname> <given-names>K.</given-names></name> <name><surname>Van Lessen</surname> <given-names>M.</given-names></name> <name><surname>Frye</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Integrin beta1 controls VE-cadherin localization and blood vessel stability.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>6</volume>:<issue>6429</issue>. <pub-id pub-id-type="doi">10.1038/ncomms7429</pub-id> <pub-id pub-id-type="pmid">25752958</pub-id></citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>Z.</given-names></name> <name><surname>Lin</surname> <given-names>P.</given-names></name> <name><surname>Chen</surname> <given-names>B.</given-names></name> <name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Xiao</surname> <given-names>W.</given-names></name> <name><surname>Wu</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2021</year>). <article-title>Autophagy alleviates hypoxia-induced blood-brain barrier injury via regulation of CLDN5 (claudin 5).</article-title> <source><italic>Autophagy</italic></source> <volume>17</volume> <fpage>3048</fpage>&#x2013;<lpage>3067</lpage>. <pub-id pub-id-type="doi">10.1080/15548627.2020.1851897</pub-id> <pub-id pub-id-type="pmid">33280500</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yao</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>Z. L.</given-names></name> <name><surname>Norris</surname> <given-names>E. H.</given-names></name> <name><surname>Strickland</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Astrocytic laminin regulates pericyte differentiation and maintains blood brain barrier integrity.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>5</volume>:<issue>3413</issue>. <pub-id pub-id-type="doi">10.1038/ncomms4413</pub-id> <pub-id pub-id-type="pmid">24583950</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>G.</given-names></name> <name><surname>Wang</surname> <given-names>L. G.</given-names></name> <name><surname>Han</surname> <given-names>Y.</given-names></name> <name><surname>He</surname> <given-names>Q. Y.</given-names></name></person-group> (<year>2012</year>). <article-title>clusterProfiler: an R package for comparing biological themes among gene clusters.</article-title> <source><italic>OMICS: J. Integrat. Biol.</italic></source> <volume>16</volume> <fpage>284</fpage>&#x2013;<lpage>287</lpage>. <pub-id pub-id-type="doi">10.1089/omi.2011.0118</pub-id> <pub-id pub-id-type="pmid">22455463</pub-id></citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>An</surname> <given-names>Q.</given-names></name> <name><surname>Wang</surname> <given-names>T.</given-names></name> <name><surname>Gao</surname> <given-names>S.</given-names></name> <name><surname>Zhou</surname> <given-names>G.</given-names></name></person-group> (<year>2018</year>). <article-title>Autophagy- and MMP-2/9-mediated Reduction and Redistribution of ZO-1 contribute to hyperglycemia-increased blood-brain barrier permeability during early reperfusion in stroke.</article-title> <source><italic>Neuroscience</italic></source> <volume>377</volume> <fpage>126</fpage>&#x2013;<lpage>137</lpage>.</citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Gong</surname> <given-names>P.</given-names></name> <name><surname>Zhang</surname> <given-names>J.</given-names></name> <name><surname>Mao</surname> <given-names>X.</given-names></name> <name><surname>Zhao</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>H.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Specific frequency electroacupuncture stimulation transiently enhances the permeability of the blood-brain barrier and induces tight junction changes.</article-title> <source><italic>Front. Neurosci.</italic></source> <volume>14</volume>:<issue>582324</issue>. <pub-id pub-id-type="doi">10.3389/fnins.2020.582324</pub-id> <pub-id pub-id-type="pmid">33122995</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>Z.</given-names></name> <name><surname>Sagare</surname> <given-names>A. P.</given-names></name> <name><surname>Ma</surname> <given-names>Q.</given-names></name> <name><surname>Halliday</surname> <given-names>M. R.</given-names></name> <name><surname>Kong</surname> <given-names>P.</given-names></name> <name><surname>Kisler</surname> <given-names>K.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Central role for PICALM in amyloid-beta blood-brain barrier transcytosis and clearance.</article-title> <source><italic>Nat. Neurosci.</italic></source> <volume>18</volume> <fpage>978</fpage>&#x2013;<lpage>987</lpage>.</citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>Y.</given-names></name> <name><surname>Zhou</surname> <given-names>B.</given-names></name> <name><surname>Pache</surname> <given-names>L.</given-names></name> <name><surname>Chang</surname> <given-names>M.</given-names></name> <name><surname>Khodabakhshi</surname> <given-names>A. H.</given-names></name> <name><surname>Tanaseichuk</surname> <given-names>O.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>Metascape provides a biologist-oriented resource for the analysis of systems-level datasets.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>10</volume>:<issue>1523</issue>.</citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zlokovic</surname> <given-names>B. V.</given-names></name></person-group> (<year>2008</year>). <article-title>The blood-brain barrier in health and chronic neurodegenerative disorders.</article-title> <source><italic>Neuron</italic></source> <volume>57</volume> <fpage>178</fpage>&#x2013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1016/j.neuron.2008.01.003</pub-id> <pub-id pub-id-type="pmid">18215617</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zlokovic</surname> <given-names>B. V.</given-names></name></person-group> (<year>2011</year>). <article-title>Neurovascular pathways to neurodegeneration in Alzheimer&#x2019;s disease and other disorders.</article-title> <source><italic>Nat. Rev. Neurosci.</italic></source> <volume>12</volume> <fpage>723</fpage>&#x2013;<lpage>738</lpage>. <pub-id pub-id-type="doi">10.1038/nrn3114</pub-id> <pub-id pub-id-type="pmid">22048062</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="footnote1">
<label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="https://cutadapt.readthedocs.io/en/stable/">https://cutadapt.readthedocs.io/en/stable/</ext-link></p></fn>
<fn id="footnote2">
<label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="https://daehwankimlab.github.io/hisat2/">https://daehwankimlab.github.io/hisat2/</ext-link></p></fn>
<fn id="footnote3">
<label>3</label>
<p><ext-link ext-link-type="uri" xlink:href="http://ccb.jhu.edu/software/stringtie/">http://ccb.jhu.edu/software/stringtie/</ext-link></p></fn>
<fn id="footnote4">
<label>4</label>
<p><ext-link ext-link-type="uri" xlink:href="http://ccb.jhu.edu/software/stringtie/gffcompare.shtml">http://ccb.jhu.edu/software/stringtie/gffcompare.shtml</ext-link></p></fn>
<fn id="footnote5">
<label>5</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.bioconductor.org/packages/release/bioc/html/ballgown.html">http://www.bioconductor.org/packages/release/bioc/html/ballgown.html</ext-link></p></fn>
<fn id="footnote6">
<label>6</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html">http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html</ext-link></p></fn>
<fn id="footnote7">
<label>7</label>
<p><ext-link ext-link-type="uri" xlink:href="http://metascape.org">http://metascape.org</ext-link></p></fn>
<fn id="footnote8">
<label>8</label>
<p><ext-link ext-link-type="uri" xlink:href="http://string-db.org/">http://string-db.org/</ext-link></p></fn>
</fn-group>
</back>
</article>
