Int J Med Sci 2023; 20(5):627-638. doi:10.7150/ijms.82503 This issue Cite

Research Paper

Transcriptomic Profiling of circRNAs in rat Hippocampus after Deep Hypothermic Circulatory Arrest

Tianlong Wang1#, Weidong Yan1#, Shengqiang Pei2, Mingru Zhang3, Qiaoni Zhang1, Yuan Teng1, Gang Liu1, Jian Wang1, Shujie Yan1, Bingyang Ji1 Corresponding address

1. Department of Cardiopulmonary Bypass, Fuwai Hospital, National Center for Cardiovascular Disease, State Key Laboratory of Cardiovascular Medicine, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing 100037, China
2. State Key Laboratory of Cardiovascular Disease, Beijing Key Laboratory for Molecular Diagnostics of Cardiovascular Diseases, Diagnostic Laboratory Service, Fuwai Hospital, National Center for Cardiovascular Diseases, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
3. Department of Anaesthesiology, Beijing Tongren Hospital, Capital Medical University, Beijing, China
# These authors contributed to this work equally.

Citation:
Wang T, Yan W, Pei S, Zhang M, Zhang Q, Teng Y, Liu G, Wang J, Yan S, Ji B. Transcriptomic Profiling of circRNAs in rat Hippocampus after Deep Hypothermic Circulatory Arrest. Int J Med Sci 2023; 20(5):627-638. doi:10.7150/ijms.82503. https://www.medsci.org/v20p0627.htm
Other styles

File import instruction

Abstract

Graphic abstract

Neurologic abnormalities occurring after deep hypothermic circulatory arrest (DHCA) remain a significant concern. However, molecular mechanisms leading to DHCA-related cerebral injury are still ill-defined. Circular RNAs (circRNAs) are a class of covalently closed non-coding RNAs and can play important roles in different types of cerebral injury. This study aimed to investigate circRNAs expression profiles in rat hippocampus after DHCA and explore the potential functions of circRNAs in DHCA-related cerebral injury. Hence, the DHCA procedure in rats was established and a transcriptomic profiling of circRNAs in rat hippocampus was done. As a result, a total of 35192 circRNAs were identified. Among them, 339 circRNAs were dysregulated, including 194 down-regulated and 145 up-regulated between DHCA and sham group. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed based on the host genes of all dysregulated circRNAs. Also, 4 circRNAs were validated by RT-qPCR (rno_circ_0028462, rno_circ_0037165, rno_circ_0045161 and rno_circ_0019047). Then a circRNA-microRNA (miRNA) interaction network involving 4 candidate circRNAs was constructed. Furthermore, functional enrichment analysis of the miRNA-targeting mRNAs of every candidate circRNA was conducted to gain insight into each of the 4 circRNAs. Our study provided a better understanding of circRNAs in the mechanisms of DHCA-related cerebral injury and some potential targets for neuroprotection.

Keywords: Circular RNAs, Deep hypothermic circulatory arrest, Cerebral injury, Rat, Hippocampus

Introduction

Deep hypothermic circulatory arrest (DHCA) is widely used in complex cardiothoracic surgeries, such as aortic dissection and congenital heart diseases [1]. Circulatory arrest can offer a bloodless surgical field to surgeons and deep hypothermia can decrease the metabolism and protect organs against global ischemia [2]. However, when blood contacts the artificial surface of DHCA circuits, the inflammation will be triggered. During the rewarming and reperfusion period of DHCA, it would induce ischemia/reperfusion (I/R) injury [3]. The rate of neurologic complication in adults who underwent aortic surgery with DHCA was 13.9% [4]. In one recent meta-analysis, the incidence of clinical and electroencephalogram seizures among pediatric patients with DHCA was 12.9% and 14.9%, respectively [5]. Neurologic abnormalities occurring after DHCA remain a significant concern.

In clinical practice, selective antegrade or retrograde cerebral perfusion are the common approaches for neuroprotection, but it is still controversial whether they can improve the nervous system outcomes [6, 7]. In basic medical research, although, some neuroprotective agents have been prone to improve DHCA-related cerebral injury in the animal models, the clinical effect and application still need to be evaluated [8, 9]. This study aimed to deeply evaluate cerebral injury after DHCA from the perspective of genetic changes.

Circular RNAs (circRNAs) are a class of non-coding RNA characterized by a covalently closed ring structure without 5′-3′ polarity [10]. A wide range of studies have demonstrated circRNAs can act as microRNA (miRNA) sponges, RNA-binding protein interacting agents and nuclear transcriptional regulator, thus controlling gene expression at transcriptional, post-transcriptional, and translational levels [11]. Currently, several studies have shown that there were differences in the expression of circRNAs in different types of cerebral injury [12-14]. Several key circRNAs can regulate the process of cerebral I/R injury through a variety of mechanisms via miRNA-messenger RNA (mRNA) axis [15].

Hence, we need to further explore dysregulated circRNAs in the central nervous system after DHCA. The hippocampus is associated with cognition, learning, and memory. Also, it is the most vulnerable to I/R injury [16]. Here, we established the DHCA procedure in rats and analyzed the differentially expressed circRNAs (DECs) in the hippocampus. Furthermore, the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and circRNA-miRNA network pattern analyses were conducted to identify the potential roles of DECs in DHCA-related cerebral injury. Taken together, our results may bring valuable information for researchers to further explore the functional roles of circRNAs and new targets in DHCA-related cerebral injury.

Materials and Methods

Animals

The protocols received institutional review and got approval from the Institutional Animal Care and Use Committee, Fuwai Hospital, Chinese Academy of Medical Sciences (FW-2021-0005). All experimental procedures complied with the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health. Sprague-Dawley rats were kept under standard laboratory conditions, within free access to food and water (provided by the HFK Bioscience, China). Rats (age, 12-14 weeks; weight, 450-550 g) were randomly allocated into 2 groups: sham group, DHCA group (n = 5, each group).

Deep hypothermic circulatory arrest

DHCA procedures were established as previously described [8, 9]. In DHCA group, rats were anesthetized with 2% sevoflurane and incubated with a 16-G endotracheal tube. Then rats were mechanically ventilated at 80 breaths per minute with 10 mL/kg tidal volume. Mean arterial blood pressure (MAP) was monitored through the left femoral artery. The tail artery and the right external jugular vein were cannulated and connected to a DHCA circuit, which contained a reservoir, a membrane oxygenator, and a heat exchanger. The DHCA circuit was primed with 12 mL of 6% hydroxyethyl starch and 2 mL saline with 150 IU heparin. After heparinization (500 IU/kg), cardiopulmonary bypass was initiated at a flow rate of 160-180 mL/kg/min for 10 min. Blood flow was directed from the jugular vein through silicon tubes to the membrane oxygenator and back to rat via the tail artery. The target deep temperature was 18 °C. After 30 minutes of systemic cooling, DHCA was induced by draining the blood to the reservoir and lasted nearly 45 minutes, which was confirmed by MAP = 0. In the rewarming phase, the temperature increased gradually by rewarming the blood in the DHCA circuits. The rewarming phase lasted more than 60 minutes. After that the rats were subjected to 40 minutes for reperfusion to recovery. Then the DHCA circuit was weaned off within another 20 minutes. Finally, rats were ventilated for 30 minutes without cardiopulmonary bypass support. During the whole procedure, MAP was more than 50 mm Hg. Rats in the sham group were only anesthetized, cannulated, and heparinized. At the end of the experiment, euthanasia was performed by cervical dislocation under deep anesthesia. Hippocampus was harvested in liquid nitrogen and stored at -80 °C immediately for further analysis.

RNA extraction

RNA was extracted from frozen tissue using Trizol (Invitrogen, Beijing, China) following the manufacturer's protocol. RNA purity was measured by spectrophotometry, with all samples having A260/280 ≥ 1.9. The RNA samples were stored at -80°C for sequencing and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analyses.

Library construction and sequencing

Library construction and sequencing was carried out by a professional company (Novogene, Beijing, China). Briefly, ribosomal RNAs and linear RNAs were removed by using Ribo-zero rRNA removal kit (Epicentre, USA) and RNase R (Epicentre, USA) respectively to enrich circRNAs. The enriched circRNAs were randomly fragmented into 250~300 bp by RNA Fragmentation Reagents (Ambion, USA). Subsequently, first-strand cDNA and second-strand cDNA were synthesized on the fragments, following end-repair, A-tailing, ligation of adaptors and PCR amplification by using NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, USA). At last, the PCR products were purified using the AMPure XP beads (Beckman Coulter, USA), and libraries were evaluated in the Agilent 2100 bioanalyzer (Agilent, USA) and quantified using qPCR. After library construction, Illumina PE150 (Illumina, USA) was used for further sequencing.

Pre‑processing of data and genomic alignment

Raw reads were filtered by discarding reads containing adapter, ploy-N or low-quality bases. Then the clean reads were mapped to the reference genome sequence (ftp://ftp.ensembl.org/pub/release-97/gtf/rattus_norvegicus/) using HISAT2 (v2.0.5) [17], and the standard SAM files were generated.

CircRNA identification and differential circRNA expression analysis

A total of two tools, find_circ and CIRI were used to identify circRNAs [18, 19]. Since the circRNAs could not be mapped to the reference genome directly, find_circ extracted unmapped reads and selected both ends of these reads (default 20 bp) as anchor reads. Each pair of anchor reads was aligned to the reference genome sequence again. If the alignment positions of the two anchors were reversed in the linear direction, then the reads were extended to find the junction position of circRNAs. At last, the sequences on both sides containing GT/AG splicing sites would be considered circRNAs candidates. CIRI analyzed CIGAR scores and scanned paired chiastic clipping signals in SAM files. Then the candidate circRNAs could be identified based on junction reads and GT-AG splicing signals, following systematic filtering steps to eliminate potential false positives. Normalized read counts of the identified circRNAs were reported as transcripts per million (TPM). DESseq2 algorithm was conducted to obtain DECs between sham and DHCA groups. The criteria were defined as |log2 fold-change (log2FC) | ≥ 1.5 and p < 0.05.

GO annotation and KEGG pathway analyses

To detect the potential functions of DECs, we conducted GO and KEGG enrichment analyses on the host genes of all DEGs. DAVID online tool was chosen for generating analysis results [20]. GO and KEGG terms with p < 0.05 were considered significant enrichment.

RT-qPCR

Total RNA was isolated from stored RNA samples and was then converted into cDNA using PrimeScript™ RT reagent kit with gDNA Eraser (Catalog. #R323-01, Vazyme, Nanjing, China). RT-qPCR was performed using the SYBR Premix Ex Taq™ II (Catalog. #Q511-02, Vazyme, Nanjing, China) according to the manufacturer's protocol. RT-qPCR conditions included: pre-denaturation (95˚C, 5 min, 1 cycle) and PCR reaction (95˚C, 10 sec, 60˚C, 30 sec, 40 cycles in total) followed by a dissolution curve. β-actin was designed as an internal control. The 2-∆∆CT method was used to calculate the relative quantification of DEGs expression levels. Primer sequences used in this study were listed in Supplemental Table S1.

Agarose gel electrophoresis for the validated circRNAs

Agarose gel electrophoresis procedure was conducted as previously described [21, 22]. Firstly, we performed PCR enrichment of the divergent primers of 4 circular RNAs in both cDNA and gDNA by using 1.1 × T3 Super PCR Mix (TsingKe, Nanjing, China). PCR conditions included: pre-denaturation (98˚C, 30 sec, 1 cycle) and PCR reaction (98˚C, 15 sec, 60˚C, 10 sec, 72˚C, 10 sec, 35 cycles in total) followed by 72˚C for 7 min in 1 cycle. Then, the PCR products were checked on 1% agarose gel electrophoresis. The results are presented in Supplemental Figure S1.

Construction of circRNA-associated ceRNA network

Because circRNAs function as miRNA sponges, we further predicted the target miRNAs of the selected circRNAs as well as the miRNA-targeting mRNAs. CircRNA-miRNA interactions were predicted with miRanda (v3.3a), RNAhybrid (v2.1.1) and TargetScan (v8.0). Then, Cytoscape software (v3.9.1) was conducted to construct the whole circRNA-miRNA interaction network. And TargetScan (v8.0) and miRDB (v6.0), were used to identify target mRNAs of miRNAs in the circRNA-miRNA network. Finally, these target genes were subjected to GO and KEGG analyses.

Statistical analysis

Data processing and statistical analysis were performed using R (v4.2.1; https://www.R-project.org/ (accessed on 10 July 2022)). The experiment and RT-qPCR results were presented as mean ± SEM. Student's t-tests were used for comparisons be-tween the two groups. P < 0.05 were considered to be statistically significant.

Results

Transcriptomic analysis of circRNAs in rat hippocampus after DHCA

The whole workflow of our study was shown in Figure 1. In this study, rats were randomly divided into sham and DHCA groups. The physiologic parameters of the rats throughout the experimental period were illustrated in Figure 2. After experiencing 45 minutes DHCA, hippocampus tissue was harvested for circRNAs screening. Finally, a total of 35192 circRNAs were identified from all samples. The length, genomic location, and chromosome distribution of the identified circRNAs were consistent (Figure 3A-C).TPM was used as the final expression value of circRNAs. The TPM distribution of each sample was similar following data pre-treatment (Figure 3D). Then, differential expression analyses were conducted to figure out dysregulated circRNAs between the two groups. Under the criterion as |log2FC| ≥ 1.5 and p < 0.05, we found 339 significant DECs, of which 194 were down-regulated and 145 were up-regulated (Figure 3E). A heat map was generated to visualize the differential expression distribution for all DECs (Figure 3F).

Functional enrichment analysis on the differentially expressed circRNAs

To explore the potential functions of DECs, thus GO and KEGG enrichment analyses were conducted. GO analysis enriched 87 terms (p < 0.05) within the biological processes, molecular functions, and cellular components, including neuron projection developpment, negative regulation of oxidative stress-induced neuron death, positive regulation of dendrite development, and brain development. These results suggested that these DECs may play an essential role in brain and neuron development. The top 10 GO terms in each category were selected for visualization (Figure 4A). 26 KEGG pathways (p < 0.05) were enriched indicating DEGs were associated with T cell receptor signaling pathway, autophagy, cyclic adenosine monophosphate (cAMP) signaling pathway, cyclic guano-sine monophosphate (cGMP)-cGMP-dependent protein kinase (PKG) signaling pathway, and mammalian target of rapamycin (mTOR) signaling pathway. The top 15 KEGG pathways were demonstrated in Figure 4B.

 Table 1 

The characterizations of DECs selected for RT-qPCR verification.

IDlog2FCPValuesource_gene_idsource_geneChrstartendstrandgenomic_
length
circRNA_
length
Type
rno_circ_00284625.823.26E-06ENSRNOG00000016351Frrs12220435619220441299+5680563exon
rno_circ_00213405.691.91E-05ENSRNOG00000020514Sox61185722676185727014+4338395intron
rno_circ_00352205.645.23E-05n/a--4132094677132113888-19211344intergenic_region
rno_circ_00173585.549.05E-05n/a--184885317248865729+12557493intergenic_region
rno_circ_00497475.426.14E-05ENSRNOG00000008709Arhgap3283327781233296545+18733306exon
rno_circ_00280475.331.42E-03ENSRNOG00000019671Rsbn12206406385206407059+674674exon
rno_circ_00235355.247.08E-04ENSRNOG00000014779Pdcd41274628954274642449+13495320exon
rno_circ_00240625.202.31E-03ENSRNOG00000018191Oprm114350713043508768+1638874exon
rno_circ_00371655.131.45E-03ENSRNOG00000013072Plxna445965088159651296-415415intron
rno_circ_00387395.133.35E-03ENSRNOG00000060950--51267855412682441-3887292intron
rno_circ_00045515.001.76E-03ENSRNOG00000001782Osbpl11117081041670811760-1344256exon
rno_circ_00372475.005.38E-03ENSRNOG00000061080Kmt2c2462199146222632+2718246exon
rno_circ_00274634.986.43E-03ENSRNOG00000010695Pdgfc2180036669180046454+9785255intron
rno_circ_00025854.942.38E-03ENSRNOG00000010257Cuedc1107550883375513546+4713380exon
rno_circ_00300034.847.73E-03ENSRNOG00000045639RGD1560883244349434444835+9892319exon
rno_circ_0024647-4.459.57E-03ENSRNOG00000031391Ceacam1618077563980775966-327327exon
rno_circ_0054079-4.483.81E-03ENSRNOG00000018931Dis3l299402756794076328+48761188exon
rno_circ_0033398-4.487.98E-03n/a--34527734645288132+10786334intergenic_region
rno_circ_0043052-4.552.41E-03ENSRNOG00000032391Esyt26144211591144219597+8006354exon
rno_circ_0031837-4.566.18E-03ENSRNOG00000022162
ENSRNOG00000017583
Pbx3;
Mapkap1
31349685213510737+13885486intron
rno_circ_0045161-4.632.10E-03ENSRNOG00000005063Atl169233169692349697+18001274exon
rno_circ_0019047-4.635.37E-03ENSRNOG00000017724Sf3b3194339422643401075+6849241exon
rno_circ_0035410-4.691.75E-03ENSRNOG00000006509Srgap34144656178144656714-536334exon
rno_circ_0039004-4.714.72E-03ENSRNOG00000000145Pik3r35135105373135121657+16284235exon
rno_circ_0033973-4.801.02E-03ENSRNOG00000012580Ccdc14136393044263961037-30595290exon
rno_circ_0007116-5.049.41E-04ENSRNOG00000003743Dars134507454545086416-11871393exon
rno_circ_0031737-5.591.93E-05n/a--3131470577131475654+5077403intergenic_region
rno_circ_0030500-5.598.71E-05n/a--27375014973750789+640640intergenic_region
rno_circ_0004035-7.891.75E-05ENSRNOG00000030840Cadm21143243014397333-73032441exon
rno_circ_0016962-9.721.04E-06ENSRNOG00000012899Rbbp81831754823176344+862862intron

DECs, differentially expressed circRNAs; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

 Figure 1 

Workflow for the whole study. DHCA, Deep hypothermic circulatory arrest; RT-qPCR, reverse transcription-quantitative polymerase chain reaction.

Int J Med Sci Image
 Figure 2 

Physiological parameters for the rats throughout the experimental period. Temperature (A), heart rate (HR, B) and mean arterial pressure (MAP, C) were measured over time in the sham and DHCA groups (values: mean ± SEM). bpm, beat per minute. P > 0.05 was not shown.

Int J Med Sci Image
 Figure 3 

Expression profiling data of circRNAs in the rat hippocampal after DHCA. (A) Length distribution of the identified circRNAs in all samples. (B) Genomic location distribution of the identified circRNAs in all samples. (C) Chromosome distribution of the identified circRNA in all samples. (D) The boxplot of normalized expression values in all samples. (E) Volcano plot analysis for DECs (P < 0.05 and |log2FC| ≥ 1.5 are set as the cut-off criteria). Red, blue, and green points represent circRNAs that are down-regulated, up-regulated, and not significantly different. (F) Heatmap presents the expression of all DECs. Red and blue indicates increased and decreased expression, respectively. DHCA, deep hypothermic circulatory arrest; TPM, transcripts per million; DECs, differentially expressed circRNAs; log2FC, log2 fold-change.

Int J Med Sci Image

Validation of circRNA by RT-qPCR

To validate the expression profiling data, we chose 15 circRNAs in the top 30 up-regulated DECs and 15 circRNAs in the top 30 down-regulated DECs for further RT-qPCR validation by using previous hippocampus samples (n = 4, each group). The characterizations of these 30 DECs were listed in Table 1. Finally, RT-qPCR analysis revealed 4 DECs were differentially expressed in the DHCA samples compared to the sham group. Rno_circ_0028462 and rno_circ_0037165 were up-regulated, while rno_circ_0045161 and rno_circ_0019047 were down-regulated significantly (Figure 5A and B), which were also in accordance with our expression profiling results (Figure 5C). The genomic location and sequence information of the 4 validated circRNAs were visualized by using Integrated Genome Viewer and presented in Supplemental Table 2.

CircRNA-miRNA interaction network

It is well known that circRNAs can bind target miRNA and then influence the miRNA-targeted mRNA. Thus, 4 candidate circRNAs verified by RT‑qPCR were selected to construct a representative circRNA-miRNA network (Figure 6).

In addition, KEGG and GO analyses for miRNA-targeted mRNAs of every DEC were conducted to gain insight into each of the 4 DECs. GO analysis suggested that the miRNA-targeted mRNAs of the 4 DECs were mainly involved in transcription-associated terms. KEGG analysis enriched many key regulatory signaling pathways, like tumor necrosis factor (TNF), forkhead box O (FOXO), Hippo, cGMP-PKG, mitogen-activated protein kinase (MAPK) and wingless / integrated (Wnt) signaling pathways. The top 10 enriched GO terms in the three categories and the top 15 KEGG pathways were shown in Figure 7 and Figure 8 respectively.

 Figure 4 

GO and KEGG enrichment analyses on DECs in the hippocampal after DHCA. (A) Top 10 GO terms in the molecular function, cellular component and biological process categories of DECs in the rat hippocampal after DHCA. (B) Top 15 KEGG pathways of DECs in the rat hippocampal after DHCA. GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; DECs, differentially expressed circRNAs; DHCA, deep hypothermic circulatory arrest.

Int J Med Sci Image
 Figure 5 

Validation of differentially expressed circRNAs. (A-B) The high expression of rno_circRNA_0028462 and rno_circRNA_0037165 and the low expression of rno_circRNA_0045161 and rno_circRNA_0019047 were confirmed with RT-qPCR method (n = 4, each group). (C) The RT-qPCR data of 4 verified circRNAs displayed the tendency in accordance with our expression profiling results. *p < 0.05, **p < 0.001.

Int J Med Sci Image
 Figure 6 

CircRNA-miRNA interaction network. The 4 DECs (rno_circRNA_0028462, rno_circRNA_0037165, rno_circRNA_0045161 and rno_circRNA_0019047) verified by RT‑qPCR were selected to construct a representative circRNA-miRNA network. Red represents up-regulated DECs. Green represents down-regulated DECs. Blue represents the targeted miRNAs. DECs, differentially expressed circRNAs.

Int J Med Sci Image
 Figure 7 

GO and KEGG enrichment analyses of the miRNA-targeted mRNAs of 2 confirmed up-regulated DECs. (A, C) Top 10 GO terms in the molecular function, cellular component and biological process categories of the miRNA-targeted mRNAs of rno_circRNA_0028462 and rno_circRNA_0037165, respectively. (B, D) Top 15 KEGG pathways of the miRNA-targeted mRNAs of rno_circRNA_0028462 and rno_circRNA_0037165, respectively. GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; DECs, differentially expressed circRNAs.

Int J Med Sci Image
 Figure 8 

GO and KEGG enrichment analyses of the miRNA-targeted mRNAs of 2 confirmed down-regulated DECs. (A, C) Top 10 GO terms in the molecular function, cellular component and biological process categories of the miRNA-targeted mRNAs of rno_circRNA_0045161 and rno_circRNA_0019047, respectively. (B, D) Top 15 KEGG pathways of the miRNA-targeted mRNAs of rno_circRNA_0045161 and rno_circRNA_0019047, respectively. GO: Gene Ontology; KEGG: Kyoto Encyclopedia of Genes and Genomes; DECs, differentially expressed circRNAs.

Int J Med Sci Image

Discussion

To provide insights into the understanding of the DHCA pathophysiology, we analyzed the expression profiles of circRNAs in rat hippocampus via a DCHA rat model. A total of 35192 circRNAs were identified by using high-throughput sequencing technology. Among them, 339 circRNAs were dysregulated, including 194 down-regulated and 145 up-regulated between DHCA and sham groups. Furthermore, we verified 4 circRNAs by RT-qPCR (rno_circ_0028462, rno_circ_0037165, rno_circ_0045161 and rno_circ_0019047). Our results will provide new targets for interventions to attenuate cerebral injury induced by DHCA.

Various underlying mechanisms take part in cerebral injury after DHCA, including systemic inflammatory response syndrome, oxidative stress, and blood brain barrier (BBB) damage [23-25]. During the circulatory arrest period, the cerebral blood supply is suspended and restored during the reperfusion period, which increases the risk of I/R injury to brain. At present, a large number of studies have shown that circRNAs regulated ischemia-induced brain injury [15, 26]. Our results confirmed that the expression profiles of circRNAs in rat hippocampus changed after the DHCA procedure, which was consistent with other studies [27]. However, the number of total circRNAs and dysregulated circRNAs in our results was far more than in the previous study [27]. Many reasons attribute to this difference, such as the differences in sequencing methods and experimental environments. But we think that one of the reasons was most likely that our sample size in each group (n = 5) is larger than that (n = 3) in the previous study.

Via splicing machinery competition between linear and circular transcripts, circRNAs could interact with their host genes and affect their expression [28]. Additionally, some circRNAs remain in the nuclei and regulate the expression of their host genes by interacting with chromatin and/or their promoters [29]. Therefore, we examined the potential functions of our DECs by running GO and KEGG enrichment analyses on their host genes. GO results showed that DECs were associated with neuronal development. Prior studies have demonstrated that circRNAs are enriched in the mammalian brain and closely linked to neuronal development such as proliferation and differentiation of neural stem cells [30, 31], thus our result was reasonable. KEGG analysis included various signaling pathways (cAMP, cGMP-PKG, mTOR signaling pathways), which may involve in the regulation of DHCA-related cerebral injury. For example, activation of either the cAMP pathway or the cGMP pathway can phosphorylate the cAMP response element binding protein (CREB) [32]. In DHCA, increased phosphorylated CREB may attenuate the severity of brain injury [33]. The mTOR can be activated by the phosphatidylinositol 3' -kinase (PI3K)-protein kinase B (PKB/AKT) pathway, decreasing apoptosis and protecting the brain against I/R injury [34]. The mTOR also can be regulated by AMP-activated protein kinase (AMPK), resulting in the promoting of autophagy and ameliorating cerebral I/R injury [35]. In addition, our results indicated that several DECs were enriched in the autophagy pathway. So far, few studies investigated the autophagy of the central nervous system under DHCA condition. But researchers presented that the role of autophagy in DHCA-mediated lung injury varies with time recently. Autophagy showed protective inhibition at 3 hours after DHCA, but overactivated at 6 hours after DHCA, aggravating lung injury [36]. Together with the DECs in the autophagy pathway of the hippocampus, it is necessary to further evaluate the autophagy in cerebral injury induced by DHCA in the near future.

It is widely accepted that circRNAs can act as miRNA sponge or competing endogenous RNA (ceRNA) [11]. However, research on the ceRNA function of circRNAs in DHCA is limited. Therefore, we predicted the miRNAs of 4 candidate circRNAs verified by RT-qPCR (rno_circ_0028462, rno_circ_0037165, rno_circ_0045161 and rno_circ_0019047). Wang et al. analyzed the expression of miRNA in piglet hippocampus after DHCA and a number of dysregulated miRNAs were identified, such as miR-10b, miR-23a and miR-27a [37]. These miRNAs were also included in our prediction of the 4 candidate circRNAs. For example, rno_circ_0045161 was the sponge of miR-10b-3p and miR-23a-5p and rno_circ_0028462 was the sponge of miR-27a-3p. MiR-10b-3p was downregulated during focal cerebral ischemia and could alleviate I/R injury by targeting programmed cell death 5 (PDCD5) [38]. MiR-23a-5p is involved in white matter remodeling when facing stroke [39]. Besides, miR-23a-5p was up-regulated and exacerbated intestinal I/R injury by promoting oxidative stress via targeting peroxisome proliferators activated receptor alpha (PPARα) [40]. MiR-27a-3p takes part in the process against apoptosis and oxidative stress to alleviate brain I/R injury symptoms [41]. Additionally, miR-27-3p could induce expression of claudin-5 and occludin by downregulating glycogen synthase kinase 3 beta (GSK3ß) and activating the Wnt/ß-catenin signaling pathway to enhance the BBB maintenance [42]. It should be noticed that both rno_circ_0045161 and rno_circ_0037165 could act as a let-7 family sponge. Let-7 family members, which are highly conserved across species in functions and sequence, are closely associated with various brain diseases such as Parkinson's disease and traumatic brain injury [43, 44]. They also serve as important modulators to regulate several immune processes, such as promoting the polarization of macrophages into the anti-inflammatory M2 phenotype and directly targeting the cytokines IL-6 and IL-10 [45]. Besides, let-7 can reduce cell apoptosis, for example by regulating Fas expression and the sensitivity of Fas-mediated apoptosis [46].

To further explore the functions of 4 candidate circRNAs, functional enrichment analysis of the miRNA-targeting mRNAs of every candidate circRNAs were conducted. The KEGG analysis identified many signaling pathways such as TNF, FOXO, Hippo, cGMP-PKG, MAPK, and Wnt signaling pathways, which all have complex interactions. Most of these pathways are involved in the cerebral I/R injury process [47-49]. Interestingly, some of these pathways were also obtained in the functional enrichment of all DECs, like cGMP-PKG pathway. It was hypothesized that circRNAs can play the same role in DHCA both in regulating its host gene and acting as ceRNA.

There were still some limitations in our study. Although we selected a total of 30 DECs for RT-qPCR validation, only 4 candidate DECs were verified successfully. Both the designed primers and the low expression level of some circRNAs may contribute to this phenomenon. Besides, due to the limited hippocampus tissue, the downstream miRNAs and miRNA-targeting mRNAs of 4 candidate circRNAs were not selected for further verification and we also failed to explore whether the 4 candidate circRNAs would affect the expression of linear RNAs. It is needed to elucidate the underlying molecular mechanisms between these 4 candidate circRNAs and DHCA in future investigation.

Conclusion

Our study analyzed the expression profiles of circRNAs in rat hippocampus after DCHA. A total of 35192 circRNAs were identified. Among them, 339 circRNAs were dysregulated, including 194 down-regulated and 145 up-regulated between DHCA and sham groups. After RT-qPCR validation, we confirmed the high expression of rno_circ_0028462 and rno_circ_0037165 and the low expression of rno_circ_0045161 and rno_circ_0019047. Taken together, our results provided a better understanding of circRNAs in DHCA-related cerebral injury.

Supplementary Material

Supplementary figure and tables.

Attachment

Acknowledgements

This study was supported by the National Natural Science Foundation of China (Grant number: 82070479).

Author Contributions

Conceptualization, B.J. and T.W.; methodology, T.W., W.Y., S.P., M.Z., Y.T., J.W., G.L. and S.Y.; software, T.W.; validation, T.W., W.Y. and M.Z.; formal analysis, T.W and W.Y.; investigation, Q.Z. and J.B.; resources, W.Y. and S.P.; data curation, T.W and S.P.; writing—original draft preparation, T.W.; writing—review and editing, W.Y. and Q.Z.; visualization, Y.T., J.W. and S.Y.; supervision, G.L. and B.J.; project administration, B.J.; funding acquisition, W.Y. and B.J. All authors have read and agreed to the published version of the manuscript.

Competing Interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

1. Liu M, Zeng Q, Li Y, Liu G, Ji B. Neurologic recovery after deep hypothermic circulatory arrest in rats: A description of a long-term survival model without blood priming. Artif Organs. 2019;43:551-60

2. Tan SZ, Singh S, Austin NJ, Alfonso Palanca J, Jubouri M, Girardi LN. et al. Duration of deep hypothermic circulatory arrest for aortic arch surgery: is it a myth, fiction, or scientific leap?. J Cardiovasc Surg (Torino). 2022;63:243-53

3. Su XW, Undar A. Brain protection during pediatric cardiopulmonary bypass. Artif Organs. 2010;34:E91-102

4. Itagaki S, Chikwe J, Sun E, Chu D, Toyoda N, Egorova N. Impact of Cerebral Perfusion on Outcomes of Aortic Surgery: The Society of Thoracic Surgeons Adult Cardiac Surgery Database Analysis. Ann Thorac Surg. 2020;109:428-35

5. Alkhatip A, Kamel MG, Farag EM, Elayashy M, Farag A, Yassin HM. et al. Deep Hypothermic Circulatory Arrest in the Pediatric Population Undergoing Cardiac Surgery With Electroencephalography Monitoring: A Systematic Review and Meta-Analysis. J Cardiothorac Vasc Anesth. 2021;35:2875-88

6. Gatti G, Benussi B, Curro P, Forti G, Rauber E, Minati A. et al. The Risk of Neurological Dysfunctions after Deep Hypothermic Circulatory Arrest with Retrograde Cerebral Perfusion. J Stroke Cerebrovasc Dis. 2017;26:3009-19

7. Tian DH, Wan B, Bannon PG, Misfeld M, LeMaire SA, Kazui T. et al. A meta-analysis of deep hypothermic circulatory arrest alone versus with adjunctive selective antegrade cerebral perfusion. Ann Cardiothorac Surg. 2013;2:261-70

8. Liu M, Li Y, Gao S, Yan S, Zhang Q, Liu G. et al. A novel target to reduce microglial inflammation and neuronal damage after deep hypothermic circulatory arrest. J Thorac Cardiovasc Surg. 2020;159:2431-44 e7

9. Yan W, Gao S, Zhang Q, Qi J, Liu G, Teng Y. et al. AdipoRon Inhibits Neuroinflammation Induced by Deep Hypothermic Circulatory Arrest Involving the AMPK/NF-kappaB Pathway in Rats. Pharmaceutics. 2022 14

10. Patop IL, Wust S, Kadener S. Past, present, and future of circRNAs. EMBO J. 2019;38:e100836

11. Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20:675-91

12. Chen Z, Wang H, Zhong J, Yang J, Darwazeh R, Tian X. et al. Significant changes in circular RNA in the mouse cerebral cortex around an injury site after traumatic brain injury. Exp Neurol. 2019;313:37-48

13. Mehta SL, Pandi G, Vemuganti R. Circular RNA Expression Profiles Alter Significantly in Mouse Brain After Transient Focal Ischemia. Stroke. 2017;48:2541-8

14. Wang XZ, Li S, Liu Y, Cui GY, Yan FL. Construction of circRNA-Mediated Immune-Related ceRNA Network and Identification of Circulating circRNAs as Diagnostic Biomarkers in Acute Ischemic Stroke. J Inflamm Res. 2022;15:4087-104

15. Yuan L, Chen W, Xiang J, Deng Q, Hu Y, Li J. Advances of circRNA-miRNA-mRNA regulatory network in cerebral ischemia/reperfusion injury. Exp Cell Res. 2022;419:113302

16. Shim JK, Ma Q, Zhang Z, Podgoreanu MV, Mackensen GB. Effect of pregabalin on cerebral outcome after cardiopulmonary bypass with deep hypothermic circulatory arrest in rats. J Thorac Cardiovasc Surg. 2014;148:298-303

17. Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357-60

18. Gao Y, Wang J, Zhao F. CIRI: an efficient and unbiased algorithm for de novo circular RNA identification. Genome Biol. 2015;16:4

19. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A. et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333-8

20. Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC. et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022

21. Huang Q, Guo H, Wang S, Ma Y, Chen H, Li H. et al. A novel circular RNA, circXPO1, promotes lung adenocarcinoma progression by interacting with IGF2BP1. Cell Death Dis. 2020;11:1031

22. Zhang Y, Liu M, Shuidong O, Hu QL, Guo DC, Chen HY. et al. Detection and identification of avian, duck, and goose reoviruses by RT-PCR: goose and duck reoviruses are part of the same genogroup in the genus Orthoreovirus. Arch Virol. 2006;151:1525-38

23. Arvola O, Haapanen H, Herajarvi J, Anttila T, Puistola U, Karihtala P. et al. Remote Ischemic Preconditioning Reduces Cerebral Oxidative Stress Following Hypothermic Circulatory Arrest in a Porcine Model. Semin Thorac Cardiovasc Surg. 2016;28:92-102

24. Bartels K, Ma Q, Venkatraman TN, Campos CR, Smith L, Cannon RE. et al. Effects of deep hypothermic circulatory arrest on the blood brain barrier in a cardiopulmonary bypass model-a pilot study. Heart Lung Circ. 2014;23:981-4

25. Tu LN, Timms AE, Kibiryeva N, Bittel D, Pastuszko A, Nigam V. et al. Transcriptome profiling reveals activation of inflammation and apoptosis in the neonatal striatum after deep hypothermic circulatory arrest. J Thorac Cardiovasc Surg. 2019;158:882-90 e4

26. Yang K, Zeng L, Ge A, Wang S, Zeng J, Yuan X. et al. A systematic review of the research progress of non-coding RNA in neuroinflammation and immune regulation in cerebral infarction/ischemia-reperfusion injury. Front Immunol. 2022;13:930171

27. Li YA, Liu ZG, Zhang YP, Hou HT, He GW, Xue LG. et al. Differential expression profiles of circular RNAs in the rat hippocampus after deep hypothermic circulatory arrest. Artif Organs. 2021;45:866-80

28. Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014;159:134-47

29. Hsiao KY, Sun HS, Tsai SJ. Circular RNA - New member of noncoding RNA with novel functions. Exp Biol Med (Maywood). 2017;242:1136-41

30. Rybak-Wolf A, Stottmeister C, Glazar P, Jens M, Pino N, Giusti S. et al. Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol Cell. 2015;58:870-85

31. Yang Q, Wu J, Zhao J, Xu T, Zhao Z, Song X. et al. Correction to: circular RNA expression profiles during the differentiation of mouse neural stem cells. BMC Syst Biol. 2019;13:14

32. Kelly MP. Cyclic nucleotide signaling changes associated with normal aging and age-related diseases of the brain. Cell Signal. 2018;42:281-91

33. Markowitz SD, Mendoza-Paredes A, Liu H, Pastuszko P, Schultz SP, Schears GJ. et al. Response of brain oxygenation and metabolism to deep hypothermic circulatory arrest in newborn piglets: comparison of pH-stat and alpha-stat strategies. Ann Thorac Surg. 2007;84:170-6

34. Hou Y, Wang K, Wan W, Cheng Y, Pu X, Ye X. Resveratrol provides neuroprotection by regulating the JAK2/STAT3/PI3K/AKT/mTOR pathway after stroke in rats. Genes Dis. 2018;5:245-55

35. You JY, Liu XW, Bao YX, Shen ZN, Wang Q, He GY. et al. A novel phosphodiesterase 9A inhibitor LW33 protects against ischemic stroke through the cGMP/PKG/CREB pathway. Eur J Pharmacol. 2022;925:174987

36. Kong M, Wei D, Li X, Zhu X, Hong Z, Ni M. et al. The dynamic changes in autophagy activity and its role in lung injury after deep hypothermic circulatory arrest. J Cell Mol Med. 2022;26:1113-27

37. Wang XH, Yao DX, Luan XS, Wang Y, Liu HX, Liu B. et al. MicroRNA expression in the hippocampal CA1 region under deep hypothermic circulatory arrest. Neural Regen Res. 2019;14:2003-10

38. Ye X, Fang H, Feng Y, Qian H, Wu Z, Chen S. MiR-10b-3p Protects Cerebral I/R Injury through Targeting Programmed Cell Death 5 (PDCD5). Crit Rev Eukaryot Gene Expr. 2021;31:85-98

39. Li Y, Liu Z, Song Y, Pan JJ, Jiang Y, Shi X. et al. M2 microglia-derived extracellular vesicles promote white matter repair and functional recovery via miR-23a-5p after cerebral ischemia in mice. Theranostics. 2022;12:3553-73

40. Li LX, Yin LH, Gao M, Xu LN, Qi Y, Peng JY. MiR-23a-5p exacerbates intestinal ischemia-reperfusion injury by promoting oxidative stress via targeting PPAR alpha. Biochem Pharmacol. 2020;180:114194

41. Li W, Zhu Q, Xu X, Hu X. MiR-27a-3p suppresses cerebral ischemia-reperfusion injury by targeting FOXO1. Aging (Albany NY). 2021;13:11727-37

42. Harati R, Hammad S, Tlili A, Mahfood M, Mabondzo A, Hamoudi R. miR-27a-3p regulates expression of intercellular junctions at the brain endothelium and controls the endothelial barrier permeability. PLoS One. 2022;17:e0262152

43. Goh SY, Chao YX, Dheen ST, Tan EK, Tay SS. Role of MicroRNAs in Parkinson's Disease. Int J Mol Sci. 2019 20

44. Sajja V, Jablonska A, Haughey N, Bulte JWM, Stevens RD, Long JB. et al. Sphingolipids and microRNA Changes in Blood following Blast Traumatic Brain Injury: An Exploratory Study. J Neurotrauma. 2018;35:353-61

45. Slota JA, Booth SA. MicroRNAs in Neuroinflammation: Implications in Disease Pathogenesis, Biomarker Discovery and Therapeutic Applications. Noncoding RNA. 2019 5

46. Wang S, Tang Y, Cui H, Zhao X, Luo X, Pan W. et al. Let-7/miR-98 regulate Fas and Fas-mediated apoptosis. Genes Immun. 2011;12:149-54

47. Gong P, Zhang Z, Zou C, Tian Q, Chen X, Hong M. et al. Hippo/YAP signaling pathway mitigates blood-brain barrier disruption after cerebral ischemia/reperfusion injury. Behav Brain Res. 2019;356:8-17

48. Xu D, Li F, Xue G, Hou K, Fang W, Li Y. Effect of Wnt signaling pathway on neurogenesis after cerebral ischemia and its therapeutic potential. Brain Res Bull. 2020;164:1-13

49. Zhan L, Li D, Liang D, Wu B, Zhu P, Wang Y. et al. Activation of Akt/FoxO and inactivation of MEK/ERK pathways contribute to induction of neuroprotection against transient global cerebral ischemia by delayed hypoxic postconditioning in adult rats. Neuropharmacology. 2012;63:873-82

Author contact

Corresponding address Corresponding author: Bingyang Ji, MD, PhD, Department of Cardiopulmonary Bypass, National Center for Cardiovascular Disease and Fuwai Hospital, No. 167 Beilishi Road, Xicheng District, 10010, Beijing, China. Phone: +86-10-88398285. Email: jibingyangcom


Received 2023-1-9
Accepted 2023-3-15
Published 2023-4-1


Citation styles

APA
Wang, T., Yan, W., Pei, S., Zhang, M., Zhang, Q., Teng, Y., Liu, G., Wang, J., Yan, S., Ji, B. (2023). Transcriptomic Profiling of circRNAs in rat Hippocampus after Deep Hypothermic Circulatory Arrest. International Journal of Medical Sciences, 20(5), 627-638. https://doi.org/10.7150/ijms.82503.

ACS
Wang, T.; Yan, W.; Pei, S.; Zhang, M.; Zhang, Q.; Teng, Y.; Liu, G.; Wang, J.; Yan, S.; Ji, B. Transcriptomic Profiling of circRNAs in rat Hippocampus after Deep Hypothermic Circulatory Arrest. Int. J. Med. Sci. 2023, 20 (5), 627-638. DOI: 10.7150/ijms.82503.

NLM
Wang T, Yan W, Pei S, Zhang M, Zhang Q, Teng Y, Liu G, Wang J, Yan S, Ji B. Transcriptomic Profiling of circRNAs in rat Hippocampus after Deep Hypothermic Circulatory Arrest. Int J Med Sci 2023; 20(5):627-638. doi:10.7150/ijms.82503. https://www.medsci.org/v20p0627.htm

CSE
Wang T, Yan W, Pei S, Zhang M, Zhang Q, Teng Y, Liu G, Wang J, Yan S, Ji B. 2023. Transcriptomic Profiling of circRNAs in rat Hippocampus after Deep Hypothermic Circulatory Arrest. Int J Med Sci. 20(5):627-638.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Popup Image