- 1Department of Neurology, Sichuan Provincial People’s Hospital, University of Electronic Science and Technology of China, Chengdu, China
- 2Department of Neurology, Clinical Medical College, The First Affiliated Hospital of Chengdu Medical College, Chengdu, China
- 3School of Bioscience and Technology, Chengdu Medical College, Chengdu, China
- 4International Clinical Research Center, Chengdu Medical College, Chengdu, China
- 5School of Public Health, Chengdu Medical College, Chengdu, China
Stroke is one of the major causes of death and long-term disability, of which acute ischemic stroke (AIS) is the most common type. Although circular RNA (circRNA) expression profiles of AIS patients have been reported to be significantly altered in blood and peripheral blood mononuclear cells, the role of exosome-containing circRNAs after AIS is still unknown. Plasma exosomes from 10 AIS patients and 10 controls were isolated, and through microarray and bioinformatics analysis, the profile and putative function of circRNAs in the plasma exosomes were studied. A total of 198 circRNAs were differentially quantified (|log2 fold change| ≥ 1.00, p < 0.05) between AIS patients and controls. The levels of 12 candidate circRNAs were verified by qRT-PCR, and the quantities of 10 of these circRNAs were consistent with the data of microarray. The functions of host genes of differentially quantified circRNAs, including RNA and protein process, focal adhesion, and leukocyte transendothelial migration, were associated with the development of AIS. As a miRNA sponge, differentially quantified circRNAs had the potential to regulate pathways related to AIS, like PI3K-Akt, AMPK, and chemokine pathways. Of 198 differentially quantified circRNAs, 96 circRNAs possessing a strong translational ability could affect cellular structure and activity, like focal adhesion, tight junction, and endocytosis. Most differentially quantified circRNAs were predicted to bind to EIF4A3 and AGO2—two RNA-binding proteins (RBPs)—and to play a role in AIS. Moreover, four of ten circRNAs with verified levels by qRT-PCR (hsa_circ_0112036, hsa_circ_0066867, hsa_circ_0093708, and hsa_circ_0041685) were predicted to participate in processes of AIS, including PI3K-Akt, AMPK, and chemokine pathways as well as endocytosis, and to be potentially useful as diagnostic biomarkers for AIS. In conclusion, plasma exosome-derived circRNAs were significantly differentially quantified between AIS patients and controls and participated in the occurrence and progression of AIS by sponging miRNA/RBPs or translating into proteins, indicating that circRNAs from plasma exosomes could be crucial molecules in the pathogenesis of AIS and promising candidates as diagnostic biomarkers and therapeutic targets for the condition.
Introduction
As a common vascular disease, stroke is one of the major causes of death and long-term disability, challenging the current healthcare system worldwide (Hu et al., 2017; Szegedi et al., 2017). Of all stroke cases, ischemic stroke (IS) accounts for more than 80% of them (Benjamin et al., 2018). In humans, IS results from various factors or diseases, including hypertension, large-vessel atherosclerosis and rupture of an atherosclerotic plaque, cardioembolism, and lacunar infarcts induced by small-vessel disease. However, the underlying mechanism and relationship of these factors or diseases with stroke remain incompletely known, which restricts safer and more efficient treatment to stroke. To treat patients with IS, early reperfusion of blocked arteries is most frequently performed, even though it is constrained by a narrow time window and risk of bleeding (Snow, 2016; Alberts, 2017; Powers et al., 2018). Therefore, novel targets for diagnostics or therapeutics are urgently needed for patients with IS.
Exosomes, nanosized vesicles (∼30–150 nm) with lipid bilayer membranes, are secreted by various cells and naturally present in the blood (Caby et al., 2005; Couzin, 2005; Raposo and Stoorvogel, 2013). When incorporating into target recipient cells, exosomes release cargos, like proteins, lipids, DNA, mRNAs, and noncoding RNAs (Hon et al., 2017; Skotland et al., 2017; Mousavi et al., 2019), into the cellular microenvironment of recipient cells (Cocucci et al., 2009; Simons and Raposo, 2009). Therefore, exosomes play a crucial role in cell-to-cell communication and intercellular signal transduction (Zhang et al., 2017) and are emerging as promising targets for the pathogenesis of diseases and potent tools for diagnosis and therapy. Among the various cargos in exosomes, circular RNA (circRNA) is attracting more attention due to its unique structure. CircRNAs are a class of noncoding RNAs with covalent bonds between 3ʹ head and 5ʹ tail ends to produce a circular pattern (Memczak et al., 2013), a more stable one than that of linear RNAs (Suzuki et al., 2006). Moreover, circRNAs are highly conserved and abundantly expressed in human cells (Rybak-Wolf et al., 2015) and characterized by distinct tissue-specific expression (Memczak et al., 2013). The high stability, evolutionary conservation, and abundance of circRNAs in various species endow them with numerous different potential functions in the regulation of gene expression, such as miRNA sponges, interactions with RNA-binding proteins (RBPs), or translation into proteins (Han et al., 2018). Consequently, circRNAs have been considered important biological regulators for understanding the molecular mechanisms of disease and identifying effective diagnostic biomarkers or therapeutic targets (Liu W. et al., 2017).
CircRNAs have been reported to be involved in the progression of many diseases, including Alzheimer’s disease (Lukiw, 2013), cardiovascular diseases (Du et al., 2017a; Jiang et al., 2019), atherosclerosis (Burd et al., 2010), stroke (Peng et al., 2019; Li et al., 2021), and various cancers (Li et al., 2015; Yao et al., 2017). Recently, circRNAs have been considered to be potential biomarkers or targets for diagnosis or treatment of acute IS (AIS) (Lu et al., 2020; Zuo et al., 2020). CircUCK2, as a miRNA sponge, markedly reduced IS infarct volume and improved neurological impairment in the mouse model of focal cerebral ischemia and reperfusion (Chen W. et al., 2020), while circHECTD1 aggravated cerebral infarction volume and neuronal apoptosis in the mouse model of IS (Dai et al., 2021). Moreover, circRNAs exhibit dynamic expression patterns in a series of physiological and pathological conditions and possess important regulatory effects. Expression profiles of circRNAs from peripheral blood mononuclear cells and blood of patients with AIS are significantly different from those of controls (Li S. et al., 2020; Dong et al., 2020; Ostolaza et al., 2020).
However, the profile and regulatory network of circRNAs derived from plasma exosomes of patients with AIS remain unclear. In this prospective investigation, a case-controlled study was conducted, and the circRNA profile of plasma exosomes was compared between AIS patients and controls through microarray. By bioinformatics approaches, the functions and signaling pathways of differentially quantified circRNAs were analyzed. Based on the results of microarray, 12 circRNAs with high fold changes among differentially quantified circRNAs were chosen for validation by qRT-PCR, and their roles in the pathogenesis of AIS and diagnostic values for AIS were also explored. The data of the present study shed light on the possible roles of exosome-derived circRNAs in regulating the pathogenesis of AIS and provide new potential targets for diagnosis and therapy of AIS.
Materials and Methods
Study Subjects
Ten adult patients with AIS (five males and five females) were selected from the First Affiliated Hospital of Chengdu Medical College in 2020. AIS was diagnosed according to the Guidelines for the Prevention of Stroke in Patients with Stroke and Transient Ischemic Attack (Kernan et al., 2014). Ten controls approximately age-matched were enrolled from the health management center of the First Affiliated Hospital of Chengdu Medical College in 2020. The demographic description of AIS patients and controls is shown in Table 1. The investigation was approved by the Ethical Committee of First Affiliated Hospital of Chengdu Medical College. All study participants or their authorized representatives signed informed consent forms.
Inclusion criteria are as follows: (1) aged 50–90 years old, (2) within 24 h after onset, and (3) have agreed to sign the informed consent. Exclusion criteria are as follows: (1) have received thrombolysis and/or thrombectomy treatment and (2) have other concomitant diseases, including acute infection, immune diseases, neurodegenerative diseases, and tumors.
Isolation and Characterization of Plasma Exosomes
Blood samples from 10 AIS patients and 10 controls included in the present study were collected into EDTA tubes by venipuncture. After centrifuging at 2,000 rcf for 15 min at room temperature, a clear top layer was harvested to a labeled tube (BS-50-M, Biosharp, China) and stored at −80°C. Plasma exosomes were isolated using the exoRNeasy Serum/Plasma Maxi kit (77064, Qiagen Sciences Inc., Germany) following the manufacturer’s instruction. The resulting exosome pellet was suspended in PBS for further RNA or protein extraction. The exosome fraction isolated from plasma was characterized by western blot, nanoparticle tracking analysis (NTA), and transmission electron microscopy (TEM).
Western Blot
Plasma exosomes suspended in PBS were mixed with the same volume of RIPA lysis buffer, and exosome proteins were purified from each sample and separated by SDS-PAGE gel and then transferred to a polyvinylidene fluoride (PVDF) membrane. The membrane was blocked using 5% BSA and then treated with antibodies against CD63 (ab134045, Abcam, United Kingdom) or TSG101 (ab133586, Abcam, United Kingdom) at 4°C for 12 h. Horseradish peroxidase-conjugated secondary antibody (SBI, United States) was then applied, and the blots were developed with enhanced chemiluminescence reagents.
NTA
Isolated exosomes were diluted in PBS and analyzed using the ZetaView (S/N 17-310, Particle Metrix, Germany) with NTA software (ZetaView 8.04.02). Triplicate measurements were recorded for each sample. Size distribution and concentration profiles were averaged to derive the representative size distribution profiles.
TEM
Plasma exosomes suspended in PBS were fixed in 2% paraformaldehyde. The fixed sample was absorbed onto formvar-coated copper grids for 5 min at room temperature. After being rinsed in distilled water, samples were stained with methyl cellulose uranyl acetate for 1 min at room temperature. Excess liquid was wicked off the grid using filter paper, and grids were stored at room temperature until imaging. Imaging was performed using Tecnai G2 Spirit BioTwin (FEI, United States).
Total RNA Extraction From Plasma Exosomes
By exoRNeasy Serum/Plasma Maxi Kit (77064, Qiagen Sciences Inc., Germany), total RNA was extracted from plasma exosomes isolated from the 10 AIS patients and 10 controls following the manufacturer’s recommendation. RNA integrity of total RNA was inspected by an Agilent Bioanalyzer 4200 (Agilent Technologies, Santa Clara, CA, United States) (data not shown).
CircRNA Microarray Profiling
SurePrint G3 Human ceRNA microarray (4 × 180K, design ID: 085499) was made by Agilent Technologies and contained probes interrogating 84,569 circRNAs. Total RNAs were incubated with RNase R to eliminate linear RNAs. The enriched circRNAs were transcribed into fluorescent complementary RNA (cRNA), which then were hybridized onto the circRNA microarray. After hybridization, the signal of the circRNA microarray was scanned using the Agilent Microarray Scanner (Agilent Technologies Inc., Santa Clara, CA, United States).
Microarray Data Analysis
The scanned information was extracted using the Agilent Feature Extraction software 10.7 (Agilent Technologies Inc., Santa Clara, CA, United States). Raw data were normalized by quantile algorithm, R package “limma”. Differentially quantified circRNAs were identified through volcano plot filtering. Hierarchical clustering analysis was conducted to display the circRNA profile pattern between patients with AIS and controls.
Validation of Selected Differentially Quantified circRNAs by qRT-PCR
The microarray results were verified by qRT-PCR using samples from the same subjects (10 AIS patients and 10 controls). Total RNA from plasma exosomes was incubated with RNase R to eliminate linear RNAs and reverse transcribed by RevertAid First-Strand cDNA Synthesis Kit (K1622, Thermo Scientific, United States). CircRNA levels were determined by qRT-PCR via TB Green® Premix Ex Taq™ II (RR820Q, TaKaRa, Japan) with the BIO-RAD system (MiniOpticon, United States) according to the instruction of the manufacturer. The reaction system was 25 μl in volume and consisted of 12.5 μl of TB Green Premix Ex Taq II, 2 μl of cDNA, 1 μl of forward and reverse primers (10 μM), and 8.5 μl of RNase-free water. The optimum reaction conditions were as follows: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. Twelve circRNAs with high fold changes from the microarray results were chosen for validation by qRT-PCR with GAPDH as a normalization reference. The specific primers for target circRNAs are described in Table 2 and synthesized by Sangon Biotech (Shanghai, China). The primers for GAPDH were purchased from Sangon Biotech (B662104, Sangon Biotech, China). All reactions were performed in triplicate. The threshold cycle (Ct) was determined using the default threshold settings, and the average Ct value was used to calculate the relative fold changes of the 12 circRNAs by the 2−ΔΔCt method (Livak and Schmittgen, 2001).
Bioinformatics Analysis and Target Prediction
Host/target genes of differentially quantified circRNAs were illuminated by Gene Ontology (GO, http://www.geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG, http://kobas.cbi.pku.edu.cn/) analyses. Target miRNAs of circRNAs were predicted using miRanda (http://www.microrna.org/microrna/home.do) and TargetScan (http://www.targetscan.org/vert_72/). MiRNA response elements as regulators of mRNAs were predicted through TargetScan (http://www.targetscan.org/vert_72/) and microT-CDS (http://www.microrna.gr/microT-CDS). All circRNA–miRNA–mRNA interaction networks of differentially quantified circRNAs were visualized by Cytoscape 3.8.2. The tissue specificity of differentially quantified circRNAs was annotated via TSCD (http://gb.whu.edu.cn/TSCD/). The types and sites of RBPs binding to differentially quantified circRNAs were annotated by circinteractome (https://circinteractome.nia.nih.gov/). The probability of differentially quantified circRNAs coding for proteins was estimated by circBank (http://www.circbank.cn/), in which circRNAs with a coding probability >0.99 were analyzed.
Receiver Operating Characteristic Curve Analysis
The levels of differentially quantified circRNAs verified by qRT-PCR were applied to generate ROC curves for the 10 AIS patients and 10 controls. The area under the curve (AUC) was calculated to assess the predictive value of the selected circRNAs for AIS diagnosis.
Statistical Analysis
The differentially quantified circRNAs were defined as |log2 fold change| ≥ 1.00 and p < 0.05. The p-value denotes the significance of GO and KEGG pathway analyses (p < 0.05). All statistical analyses were performed by SPSS (26.0, IBM, United States). All data were represented as the mean ± standard deviation (
Results
Characterization of Plasma Exosome
The exosomes isolated from plasma of patients with AIS and controls had obvious detectable levels of CD63 and TSG101 by western blot (Figure 1A). The typical bilayer of exosomes was observed (Figures 1B,C) through TEM, and the median size of exosomes isolated from plasma was ∼137 nm for controls and ∼141 nm for patients with AIS (Figures 1D,E, respectively). These data confirmed that we successfully isolated plasma exosomes from blood.
FIGURE 1. Identification and characterization of plasma exosomes. (A) Detection of CD63 and TSG101 as exosome markers by western blot. (B,C) Representative images of transmission electron microscopy. (D,E) Representative size of exosomes from controls and AIS patients, respectively.
Profile of CircRNAs Derived From Plasma Exosomes
The profile of 84,569 human circRNAs was generated by SurePrint G3 Human ceRNA microarray using plasma exosomes from patients with AIS and controls. After fold-change filtering, 198 circRNAs were differentially quantified between patients with AIS and controls, including 72 upregulated (log2 fold change ≥ 1.00, p < 0.05) and 126 downregulated (log2 fold change ≤ −1.00, p < 0.05) circRNAs. Hierarchical clustering of the levels between AIS patients and controls is illustrated in Figure 2A. The volcano plot showed that significant variation was observed in circRNA levels between AIS patients and controls (Figure 2B).
FIGURE 2. CircRNA profiles of plasma exosomes from AIS patients (AIS) and controls (CON). (A) Heat map of circRNAs differentially quantified between AIS patients and controls. (B) Volcano plot of circRNA profiles between AIS patients and controls. Red dots represent significantly upregulated circRNAs (log2 fold change ≥ 1.00 and p < 0.05, n = 10, independent two-sample t-test), and green dots represent significantly downregulated circRNAs (log2 fold change ≤ −1.00 and p < 0.05, n = 10, independent two-sample t-test).
The distribution patterns of differentially quantified circRNAs in the chromosomes are shown in Figure 3A, indicating that the profile of circRNAs in patients with AIS was significantly different from that in controls. Moreover, the tissue-specific analysis showed that the differentially quantified circRNAs were mainly derived from the digestive system (esophagus, esophagogastric junction, intestine, and liver), the cardiovascular system (blood vessel and heart), and lungs (Figure 3B).
FIGURE 3. Characterization of circRNA profiles between AIS patients and controls. (A) Bar diagram of the distribution of differentially quantified circRNAs in different chromosomes. (B) Top 15 tissue-specific distributions of differentially quantified circRNAs.
GO and KEGG Pathway Analyses of Host Genes of Differentially Quantified circRNAs
The different profiles of circRNAs between patients with AIS and controls could be the result of disordered expression of host genes generating circRNAs; thus, GO and KEGG pathway analyses were performed to evaluate the putative functions of host genes. According to the GO analysis shown in Figure 4A, the host genes of differentially quantified circRNAs, in the cellular component, were mainly involved in ribonucleoprotein (ribonucleoprotein complex, cytoplasmic ribonucleoprotein granule, and ribonucleoprotein granule) and cell junction (focal adhesion and cellsubstrate junction). The molecular functions related to host genes mainly included translation regulator activity and binding (ribosome, ribonucleoprotein, and cadherin), and the major biological processes of host genes were associated with RNA regulation. The KEGG analysis showed that intracellular physiological processes (like RNA transport and protein processing in endoplasmic reticulum) and intercellular interaction (like focal adhesion and leukocyte transendothelial migration) could be affected, indicating cellular functions disrupted in AIS (Figure 4B).
FIGURE 4. GO annotations and KEGG pathway analyses for host genes of the differentially quantified circRNAs. (A) GO annotations. The GO analysis categorized the mRNAs into different groups under the theme of biological process, cellular component, and molecular function. The gene count for each GO category is listed on the left. (B) Top 10 KEGG pathways. The vertical axis indicates the pathway category, and the horizontal axis indicates the −log10 p-value of a pathway (Fisher’s exact test).
Validation of Levels of Differentially Quantified circRNAs
To validate the results from microarray, the levels of 12 differentially quantified circRNAs with high fold changes were determined using qRT-PCR. Among 12 selected differentially quantified circRNAs, the levels of 10 circRNAs were significantly different (p < 0.05) between AIS patients and controls, including five upregulated (hsa_circ_0066867, hsa_circ_0093708, hsa_circ_0032222, hsa_circ_0059662, and hsa_circ_0041685) and five downregulated (hsa_circ_0131433, hsa_circ_0123103, hsa_circ_0112036, hsa_circ_0113001, and hsa_circ_0050840) circRNAs (Figures 5A,B). The data showed that the levels of differentially quantified circRNAs by qRT-PCR were consistent with the results from microarray, supporting its reliability.
FIGURE 5. Validation and network diagram of differentially quantified circRNAs. (A,B) Relative fold changes of 12 differentially quantified circRNAs using qRT-PCR (n = 10). *, p < 0.05. (C) Network of circRNA-miRNA of 198 differentially quantified circRNAs. (D) Network of circRNA-miRNA-mRNA of the 10 differentially quantified circRNAs verified by qRT-PCR.
Predicted circRNA-miRNA-mRNA Networks of Differentially Quantified circRNAs
As one of the mechanisms regulating gene expression, circRNA sponges miRNA and represses the effects of miRNA on mRNA. The miRNAs predicted to be sponged by 198 differentially quantified circRNAs formed an intricate network (Figure 5C). To explore the potential pathways involved with the differentially quantified circRNAs, verified differentially quantified circRNAs, including five upregulated (hsa_circ_0066867, hsa_circ_0093708, hsa_circ_0032222, hsa_circ_0059662, and hsa_circ_0041685) and five downregulated (hsa_circ_0131433, hsa_circ_0123103, hsa_circ_0112036, hsa_circ_0113001, and hsa_circ_0050840) circRNAs, were analyzed. Figure 5D showed the predicted network of circRNA–miRNA–mRNA, suggesting that verified differentially quantified circRNAs could regulate the expression of various target genes by sponging multiple miRNAs.
GO and KEGG Pathway Analyses of Target Genes of Differentially Quantified circRNAs
Then, GO and KEGG pathway analyses of target genes of 10 differentially quantified circRNAs verified by qRT-PCR were performed to study the potential roles of these circRNAs in AIS. As to the cellular component of GO annotation, the target genes were enriched in the nucleus (nucleus, nucleoplasm, and nuclear lumen), organelle (intracellular/intracellular membrane-bounded/membrane-bounded organelle), and membrane (intrinsic/integral component of membrane) (Figure 6A). The molecular functions were mainly related to DNA binding and cis-regulatory region, and the biological processes were associated with the regulation of macromolecules and nucleic acids (RNA and DNA). The KEGG analysis showed that metabolic (AMPK signaling pathway and PI3K-Akt signaling pathway), endocrine-related (vasopressin-regulated water reabsorption, adrenergic signaling in cardiomyocytes, and relaxin signaling pathway), and inflammatory (chemokine signaling pathway) pathways were enriched (Figure 6B).
FIGURE 6. GO annotations and KEGG pathway analyses for target genes of the 10 differentially quantified circRNAs verified by qRT-PCR. (A) GO annotations. The GO analysis categorized the mRNAs into different groups under the theme of biological process, cellular component, and molecular function. The gene count for each GO category is listed on the left. (B) Top 10 KEGG pathways. The vertical axis indicates the pathway category, and the horizontal axis indicates the −log10 p-value of a pathway (Fisher’s exact test).
According to circRNA–miRNA–mRNA networks, by GO and KEGG analyses, several differentially quantified circRNAs verified by qRT-PCR were involved in processes and pathways related to AIS. Hsa_circ_0112036 sponging miR-24-3p and regulating the levels of STRADB/BCL2L11/CDKN1B and hsa_circ_0066867 impacting miR-182-5p-FOXO3 were predicted in the AMPK and PI3K-Akt pathways. Hsa_circ_0093708, via miR-4269-CREB5 or miR-4533-AQP4, and hsa_circ_0041685, via miR-4267-DCTN2, played roles in the vasopressin-regulated water reabsorption. Moreover, hsa_circ_0066867 or hsa_circ_0041685 potentially regulated chemokine signaling through miR-6737-5p-CCL2 or miR-3192-5p-CXCL12, respectively.
Coding Probability of Differentially Quantified circRNAs
Certain circRNAs could be translated into proteins and participate in many processes of diseases. Thus, the coding probability of the differentially quantified circRNAs was analyzed. Of all 198 differentially quantified circRNAs, 96 circRNAs had strong coding potential (coding probability > 0.99), indicating that these molecules, through translation, could contribute to the process of AIS (Figure 7A).
FIGURE 7. GO annotations and KEGG pathway analyses for proteins potentially coded by differentially quantified circRNAs. (A) Distribution of the coding probability of differentially quantified circRNAs. (B) GO annotations. The GO analysis categorized the genes into different groups under the theme of biological process, cellular component, and molecular function. The gene count for each GO category is listed on the left. (C) Top 10 KEGG pathways. The vertical axis indicates the pathway category, and the horizontal axis indicates the −log10 p-value of a pathway (Fisher’s exact test).
To further investigate the role of proteins translated from differentially quantified circRNAs with strong coding potential in AIS, GO, and KEGG pathway analyses for these proteins were performed. The results of GO analysis showed that cellular components were mainly located in the intracellular (cytoplasm and organelle) and extracellular (vesicle and exosome) compartments. The molecular functions of these proteins were related to RNA and protein binding (enzyme/cadherin/cell adhesion molecule binding), and the biological processes were mainly related to localization (protein and cell) and cellular homeostasis (cellular component organization, organelle assembly, microtubule-based process, and cytoskeleton organization) (Figure 7B). The KEGG pathway analysis showed that intracellular activities (endocytosis, RNA transport, phagosome, and protein processing in endoplasmic reticulum) and interactions between cells (focal adhesion and tight junction) could be affected by proteins predicted to be translated from circRNAs (Figure 7C). Of the 10 verified differentially quantified circRNAs, hsa_circ_0041685 was predicted to be translated into RABEP1 and to impact the process of endocytosis.
Prediction of RBPs Binding to Differentially Quantified circRNAs
Besides sponging miRNAs, circRNAs could bind to RBPs to regulate the expression of genes and processes of diseases. Among 38 RBPs, eukaryotic initiation factor 4A-3 (EIF4A3) had the highest count of binding sites for the differentially quantified circRNAs, while Argonaute 2 (AGO2) and fragile X mental retardation protein (FMRP) ranked second and third places, respectively. Moreover, of 198 differentially quantified circRNAs, 170 circRNAs were predicted to bind to EIF4A3, 89 circRNAs to bind to AGO2, and 65 circRNAs to bind to fused in sarcoma (FUS) (Figure 8).
Evaluation of Diagnostic Value of circRNAs in AIS With ROC Curve Analysis
ROC curve analysis was applied to evaluate the potential diagnostic value of the differentially quantified circRNAs derived from plasma exosomes. Based on the results of GO and KEGG analyses above, hsa_circ_0112036, hsa_circ_0066867, hsa_circ_0093708, and hsa_circ_0041685 were chosen for this analysis. The AUC of these four circRNAs ranged from 0.760 to 0.810 (p < 0.05), with the highest AUC found for hsa_circ_0112036 (AUC = 0.810) (Figure 9).
FIGURE 9. ROC curves of four differentially quantified circRNAs verified by qRT-PCR, including (A) hsa_circ_0112036, (B) hsa_circ_0041685, (C) hsa_circ_0093708, and (D) hsa_circ_0066867 (n = 10). The red diagonal lines in the graphs represent the diagnostic threshold value of circRNAs for AIS (AUC = 0.500).
Discussion
Stroke causes high mortality and disability rates globally (Hu et al., 2017). As the most common type of stroke, IS could be the result of various factors or diseases, including advanced age, smoking, hypertension, hyperlipidemia, and diabetes mellitus (Guzik and Bushnell, 2017). It has been reported that circRNA profiles of peripheral blood and peripheral blood mononuclear cells from patients with AIS were markedly different from those of controls (Li S. et al., 2020; Dong et al., 2020), indicating that circRNAs could be regulatory factors in AIS. However, it remains unknown whether circRNAs carried by exosomes mediate the signaling transduction and participate in the occurrence and progression of AIS. Hence, the present study aimed to elucidate the profile and potential function of circRNAs derived from plasma exosomes of patients with AIS.
By microarray, 198 differentially quantified circRNAs were found between patients with AIS and controls, evidencing that the profile of circRNAs derived from plasma exosomes was affected by AIS. CircRNAs are highly conserved and characterized by a stable structure and distinct tissue-specific expression (Memczak et al., 2013; Beermann et al., 2016). A tissue-specific circRNA database showed that the major origins of differentially quantified circRNAs with tissue specificity included the digestive system, the cardiovascular system, and lungs. The diseases contributing to AIS, like hypertension, atherosclerosis, and diabetes mellitus, cause marked impairment of the cardiovascular system, disrupting the expression of circRNAs (Wu et al., 2017; Zhao et al., 2020; Fu and Sun, 2021). The aberrant profile of circRNAs specific to the digestive system could be in part the result of the eating habits of patients with AIS. According to the latest literature, a high-fat diet could impair the colonic epithelium and alter the metabolic capacity of the microbiota, increasing the level of circulating trimethylamine N-oxide and the risk of atherosclerosis (Yoo and Zieba, 2021), which could partially explain the changes in the profile of circRNAs specific to the digestive system in AIS. Moreover, changes in the expression of lung-specific circRNAs could be contributing to the altered function of platelets in AIS. Platelet activation and platelet–leukocyte aggregation are independent determinants for AIS (Schmalbach et al., 2015). The lungs have been identified as a primary site of terminal platelet production and an organ with considerable hematopoietic potential (Lefrançais et al., 2017). This suggests that dysfunctional platelet from the lung could partially account for the association between aberrant expression of pulmonary circRNAs and AIS. These results indicate that, except for the well-known cardiovascular system, the digestive system and lungs could be non-negligible aspects in the pathogenesis of AIS.
Eukaryotic circRNAs are most commonly generated by back-splicing from pre-mRNAs of host genes—a process in which a downstream 5ʹ splice site joins with an upstream 3ʹ splice site—and regulate their expression via feedback (Memczak et al., 2013; Ashwal-Fluss et al., 2014; Lasda and Parker, 2014; Ragan et al., 2019). Thus, the profile of exosomal circRNAs could be an indicator of the expression profile of host genes. According to GO annotations of host genes, the ribosome and ribonucleoprotein could be impacted in AIS, disordering RNA translation/regulation and protein processing in the endoplasmic reticulum, which was highlighted by KEGG analysis. It has been reported that endoplasmic reticulum stress induced by abnormal protein processing takes place in atherosclerosis (Tabas, 2010; Huang et al., 2018) and hypertension (Naiel et al., 2019; Liu et al., 2020), for which the underlying mechanism could be the dysfunction of vascular endothelial cells by endoplasmic reticulum stress, like excessive apoptosis (Sun et al., 2015; Carlisle et al., 2016). Moreover, by KEGG pathway analysis, focal adhesion and leukocyte transendothelial migration were significantly enriched. Focal adhesions—contact points for the cell with the extracellular matrix—regulate communication of the cell with the surrounding extracellular environment and signaling of diverse cellular processes, including proliferation, migration, apoptosis, spreading, and differentiation (Carragher and Frame, 2004), also found in patients with large-artery atherosclerotic stroke by Xiao et al. (2021). There are a variety of molecules associated with focal adhesion, such as focal adhesion kinase (FAK), integrin, talin, and vinculin, and their activities are related to the development of hypertension (Sugimura et al., 2010; Sen et al., 2011; Jia et al., 2017), atherosclerosis (von Essen et al., 2016; Murphy et al., 2019), and thrombosis (Hitchcock et al., 2008), contributing to the progression of AIS. Moreover, leukocytes, especially monocytes, migrate from blood to the subendothelial space and form foam cells with lipoproteins, a key point to aggravate inflammation of the atherosclerotic plaque (Hansson and Libby, 2006; Gisterå and Hansson, 2017). Based on the data mentioned above, the diseases or factors causing AIS, including hypertension, atherosclerosis, and thrombosis, were significantly accompanied by altered profiles of host genes and circRNAs, demonstrating that circRNAs derived from plasma exosomes could be useful biomarkers to evaluate AIS or diseases/factors related to it. In addition, the data on GO and KEGG pathway analyses for host genes of differentially quantified circRNAs were different from those of previous research on blood of patients with AIS (Li S. et al., 2020), indicating that the profile of circRNAs from diverse origins was markedly distinct.
Though the function of most circRNAs remains unknown, a few circRNAs described in physiological and pathological processes generally act as miRNA/protein sponges or translate into proteins (Hansen et al., 2013; Zheng et al., 2016; Pamudurti et al., 2017). To evaluate the function of differentially quantified circRNAs found in the present study, by in silico analysis, the putative roles of exosome-derived circRNAs were predicted in AIS.
As miRNA sponges, circRNAs bind to miRNAs and consequently inhibit their function, which is a new mechanism in the regulation of miRNA activity and gene expression (Hansen et al., 2013; Memczak et al., 2013). In the present study, the interaction between differentially quantified circRNAs and miRNAs formed a complex network, suggesting that circRNAs could regulate various stages of AIS via sponging miRNAs. MiR-382-5p, containing binding sites for hsa_circ_0069594 and hsa_circ_0128924, is involved in the NFIA signaling regulating cholesterol homeostasis and inflammation (Hu et al., 2015). MiR-92a (potentially sponged by hsa_circ_0024722, hsa_circ_0059965, and hsa_circ_0062949) and miR-221/222 (potentially sponged by hsa_circ_0139214 and hsa_circ_0108959) participated in the course of atherosclerosis by targeting KLF2 mRNA, which modulates shear stress genes (Wu et al., 2011), and sustaining plaque stability (Bazan et al., 2015), respectively. Furthermore, miR-298 (potentially sponged by hsa_circ_0097102 and hsa_circ_0014293) exacerbated cerebral ischemia/reperfusion injury by targeting ACT1 mRNA (Sun et al., 2018), while miR-210 (potentially sponged by hsa_circ_0140265 and hsa_circ_0040760) targets CASP8AP2 mRNA and ISCU1/2 mRNA to suppress apoptosis and improve cerebral injury following AIS (Chan et al., 2009; Kim et al., 2009).
To explore the function of genes regulated by the circRNA–miRNA–mRNA interaction network, we performed GO and KEGG pathway analyses for the target genes of 10 differentially quantified circRNAs verified by qRT-PCR, including five upregulated (hsa_circ_0066867, hsa_circ_0093708, hsa_circ_0032222, hsa_circ_0059662, and hsa_circ_0041685) and five downregulated (hsa_circ_0131433, hsa_circ_0123103, hsa_circ_0112036, hsa_circ_0113001, and hsa_circ_0050840) circRNAs. The results showed that, as competing endogenous RNAs, differentially quantified circRNAs, through misregulating RNA and DNA, would impact cellular components, including the nucleus, organelle, and membrane, and modify the function and fate of cells/tissues. AIS induces neuron apoptosis (Jiang et al., 2018), while autophagy could rescue neurons from apoptosis and confer neuroprotection in AIS (Wang et al., 2012; Wang et al., 2014). These processes involved organelles, like the mitochondria and endoplasmic reticulum (Chen et al., 2002; Yorimitsu et al., 2006; Heath-Engel et al., 2008; Guan et al., 2018), and signaling pathways enriched by KEGG analysis, including PI3K-Akt and AMPK pathways (Sheng et al., 2010; Jiang et al., 2015; Hou et al., 2018; Wang et al., 2020). Similarly, Li et al. 2020b also predicted that circRNAs from the blood could regulate apoptosis in AIS. Based on the circRNA–miRNA–mRNA network, hsa_circ_0112036 and hsa_circ_0066867 have been predicted to regulate mRNAs related to AMPK and PI3K-Akt signaling, like STRADB, BCL2L11, CDKN1B, and FOXO3, and be involved with the process of AIS. Moreover, KEGG analysis highlighted pathways related to vasopressin-regulated water reabsorption and chemokine signaling. It has been reported that vasopressin is associated with stroke-related edema (Ameli et al., 2014) and that aquaporin 4 (AQP4), the most ubiquitous water channel in the central nervous system and abundantly expressed in astrocytes, participates in the vasopressin-regulated water reabsorption. Zheng et al. 2017 reported that astrocyte cell injury would be ameliorated by downregulation of AQP4 expression in cerebral IS. In the present study, the hsa_circ_0093708-miR-4533-AQP4 axis was predicted, and the increased level of hsa_circ_0093708 in plasma exosomes of AIS patients suggested that AIS patients may have upregulated the expression of AQP4 and, therefore, be prone to brain injury. Chemokine signaling was indispensable for leukocyte migration, which is a crucial event in the development of inflammation and atherosclerosis (Ruytinx et al., 2018; Yan et al., 2021). In this context, we found that hsa_circ_0066867 and hsa_circ_0041685, significantly increased in AIS patients, potentially regulate chemokine signaling through miR-6737-5p-CCL2 and miR-3192-5p-CXCL12, respectively. Previous studies showed that CCL2 promoted IS via chemokine signaling (Li L. et al., 2020) and that the level of serum CXCL12 was positively correlated with stroke severity (Liu et al., 2015). Therefore, circRNAs from plasma exosomes could potentially prevent and/or assist in the treatment of AIS. Besides, the results on target genes of differentially quantified circRNAs were different from those of Dong et al. (2020), in which circRNAs of peripheral blood mononuclear cells from patients with AIS were involved in inflammation and immunity. We speculated that, compared with circRNAs from a single source, circRNAs from plasma exosomes originating from various cells/tissues could provide more comprehensive information as to AIS.
Since researchers found that certain proteins translated by circRNAs participate in the processes of human diseases (Yang et al., 2018; Zhang et al., 2018), a variety of circRNAs have been described as coding RNAs, and the resulting proteins may play biological roles in the emergence and progression of human diseases. The translational potential of differentially quantified circRNAs were assessed, and the functions of circRNA-translated proteins were also analyzed by GO and KEGG. We found that 96 of 198 differentially quantified circRNAs had strong translational potential. On one hand, the proteins that might be translated from differentially quantified circRNAs with translational potential may contribute to intracellular metabolic processes and functions of cells, including RNA binding, organelle assembly, microtubule-based process, and protein processing in the endoplasmic reticulum. On the other hand, these proteins were involved in the communication between cells via several ways, like endocytosis, phagosome, focal adhesion, and tight junction. Endocytosis—a ubiquitous physiological process mediating nutrient uptake, receptor internalization, and signaling, essential events for cell growth and survival—has been reported to cause neuronal death and exacerbate brain damage during AIS (Troulinaki and Tavernarakis, 2012; Tejeda et al., 2019; Diaz-Guerra, 2021). RABEP1 is an essential and rate-limiting component of the endosome fusion regulating early endosomal transport in endocytosis. Upregulated hsa_circ_0041685 might result in increased level of RABEP1 and accelerate endocytosis, playing a role in the development of AIS. Focal adhesion, predicted by GO and KEGG analyses of host genes as well, played an important role in various inducers of AIS, including hypertension (Jia et al., 2017), atherosclerosis (von Essen et al., 2016), and thrombosis (Hitchcock et al., 2008). Chen et al. reported that, in a mouse model, targeting of FAK could be a potential treatment for early IS (Chen et al., 2018). Tight junction is a selectively permeable barrier that generally represents the rate-limiting step of paracellular transport, and, as one of the pathophysiological features of AIS, loss of blood–brain barrier tight junction integrity results in vasogenic edema, hemorrhagic transformation, and increased mortality (Abdullahi et al., 2018; Yang et al., 2019). Thus, we suspected that exosome-derived circRNAs could be a new treatment approach to AIS. Furthermore, there were several events highlighted by GO and KEGG analyses of both host genes and proteins encoded by circRNAs, suggesting that disordered expression of host genes could affect the profile of circRNAs in AIS, and vice versa.
RBPs assemble the ribonucleoprotein complexes to bind RNA sequences by interacting with specific cis-regulatory elements and affect the expression and function of their target RNAs (Janga and Mittal, 2011). Increased evidence indicates that many circRNAs interact with RBPs (Hansen et al., 2013; Ashwal-Fluss et al., 2014; Abdelmohsen et al., 2017) and that the interactions between circRNAs and RBPs are also deemed to be an essential element involved in gene transcriptional regulation (Hansen et al., 2013; Memczak et al., 2013), circRNA translation (Yang et al., 2017), and RBP sponging (Du et al., 2017b). Therefore, the distribution of putative RBPs interacting with differentially quantified circRNAs was investigated in the present study. The data showed that the main RBPs binding to differentially quantified circRNAs were EIF4A3 and AGO2. EIF4A3, a core component of the exon junction complex, regulated neuronal cell injury by targeting cGLIS3 in AIS (Jiang et al., 2021). AGO2 is a critical component of the RNA-induced silencing complex and, as a consequence, a master regulator of miRNA-dependent gene silencing. AGO2-associated miRNA profiles were modified in the brain of a rat stroke model (Liu X. S. et al., 2017). Besides, RBPs can be critical factors for the promotion of circRNA transmission from cells (Janas et al., 2015) and can serve as intracellular inducers of circRNA loading in exosomes (O'Leary et al., 2017), which could assist circRNAs in mediating communication between cells via exosomes.
CircRNAs are stable and abundant and function as effective diagnostic biomarkers of stroke, like hsa_circ_0001599 (Li et al., 2021) and hsa_circ_0141720 (Chen Y. et al., 2020). In order to explore new circRNAs as potential biomarkers of AIS, ROC curves of four circRNAs (hsa_circ_0112036, hsa_circ_0066867, hsa_circ_0093708, and hsa_circ_0041685) predicted to play a role in the pathogenesis of AIS were analyzed. The data showed that the levels of the four circRNAs could significantly differentiate AIS patients from controls, with hsa_circ_0112036 possessing the highest AUC. These results indicated that the four circRNAs could be further explored as promising biomarkers for AIS diagnosis.
Taken together, the data of this study demonstrate that circRNAs derived from plasma exosomes were differentially quantified between AIS patients and controls. Importantly, these results revealed that exosomal circRNAs from plasma, especially hsa_circ_0112036, hsa_circ_0066867, hsa_circ_0093708, and hsa_circ_0041685, potentially participate in the progression of AIS via sponging miRNA/RBPs or encoding proteins, may be explored as biomarkers for the diagnosis of AIS, and may also be potential targets for therapeutic interventions. However, these findings are based on only 10 AIS patients and 10 controls, and more samples are necessary to assess and validate the data in the future.
Data Availability Statement
The microarray datasets generated during the current study have been deposited in the NCBI Gene Expression Omnibus (GEO) repository and are accessible via GEO series accession number GSE195442 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE195442).
Ethics Statement
The investigation was approved by the Ethical Committee of the First Affiliated Hospital of Chengdu Medical College. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JY, JH, and KC conceived the study; YG, KL, MY (sixth author), HC, and MY (eighth author) collected the samples; JH and YL performed isolation and characterization of plasma exosomes; JY, JH, and KC performed validation of circRNA levels by qRT-PCR; JY and KC analyzed the data and wrote the manuscript. All authors have read and agreed on the contents of the manuscript.
Funding
This study was supported by grants from the National Natural Science Foundation of China (81870940, 82171295), the Health Commission of Sichuan Province (18ZD008) and the Natural Science Foundation of Chengdu Medical College (CYZ17-06).
Conflict of Interest
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.
Publisher’s Note
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.
References
Abdelmohsen, K., Panda, A. C., Munk, R., Grammatikakis, I., Dudekula, D. B., De, S., et al. (2017). Identification of HuR Target Circular RNAs Uncovers Suppression of PABPN1 Translation by CircPABPN1. RNA Biol. 14 (3), 361–369. doi:10.1080/15476286.2017.1279788
Abdullahi, W., Tripathi, D., and Ronaldson, P. T. (2018). Blood-brain Barrier Dysfunction in Ischemic Stroke: Targeting Tight Junctions and Transporters for Vascular protection. Am. J. Physiology-Cell Physiol. 315 (3), C343–c356. doi:10.1152/ajpcell.00095.2018
Alberts, M. J. (2017). Stroke Treatment with Intravenous Tissue-type Plasminogen Activator. Circulation 135 (2), 140–142. doi:10.1161/circulationaha.116.025724
Ameli, P. A., Ameli, N. J., Gubernick, D. M., Ansari, S., Mohan, S., Satriotomo, I., et al. (2014). Role of Vasopressin and its Antagonism in Stroke Related Edema. J. Neurosci. Res. 92 (9), 1091–1099. doi:10.1002/jnr.23407
Ashwal-Fluss, R., Meyer, M., Pamudurti, N. R., Ivanov, A., Bartok, O., Hanan, M., et al. (2014). circRNA Biogenesis Competes with Pre-mRNA Splicing. Mol. Cel 56 (1), 55–66. doi:10.1016/j.molcel.2014.08.019
Bazan, H. A., Hatfield, S. A., O’Malley, C. B., Brooks, A. J., Lightell, D., and Woods, T. C. (2015). Acute Loss of miR-221 and miR-222 in the Atherosclerotic Plaque Shoulder Accompanies Plaque Rupture. Stroke 46 (11), 3285–3287. doi:10.1161/strokeaha.115.010567
Beermann, J., Piccoli, M.-T., Viereck, J., and Thum, T. (2016). Non-coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol. Rev. 96 (4), 1297–1325. doi:10.1152/physrev.00041.2015
Benjamin, E. J., Virani, S. S., Callaway, C. W., Chamberlain, A. M., Chang, A. R., Cheng, S., et al. (2018). Heart Disease and Stroke Statistics-2018 Update: A Report from the American Heart Association. Circulation 137 (12), e67–e492. doi:10.1161/cir.0000000000000558
Burd, C. E., Jeck, W. R., Liu, Y., Sanoff, H. K., Wang, Z., and Sharpless, N. E. (2010). Expression of Linear and Novel Circular Forms of an INK4/ARF-Associated Non-coding RNA Correlates with Atherosclerosis Risk. Plos Genet. 6 (12), e1001233. doi:10.1371/journal.pgen.1001233
Caby, M.-P., Lankar, D., Vincendeau-Scherrer, C., Raposo, G., and Bonnerot, C. (2005). Exosomal-like Vesicles Are Present in Human Blood Plasma. Int. Immunol. 17 (7), 879–887. doi:10.1093/intimm/dxh267
Carlisle, R. E., Werner, K. E., Yum, V., Lu, C., Tat, V., Memon, M., et al. (2016). Endoplasmic Reticulum Stress Inhibition Reduces Hypertension through the Preservation of Resistance Blood Vessel Structure and Function. J. Hypertens. 34 (8), 1556–1569. doi:10.1097/hjh.0000000000000943
Carragher, N. O., and Frame, M. C. (2004). Focal Adhesion and Actin Dynamics: a Place where Kinases and Proteases Meet to Promote Invasion. Trends Cel Biol. 14 (5), 241–249. doi:10.1016/j.tcb.2004.03.011
Chan, S. Y., Zhang, Y.-Y., Hemann, C., Mahoney, C. E., Zweier, J. L., and Loscalzo, J. (2009). MicroRNA-210 Controls Mitochondrial Metabolism during Hypoxia by Repressing the Iron-Sulfur Cluster Assembly Proteins ISCU1/2. Cel Metab. 10 (4), 273–284. doi:10.1016/j.cmet.2009.08.015
Chen, D., Wei, L., Liu, Z.-R., Yang, J. J., Gu, X., Wei, Z. Z., et al. (2018). Pyruvate Kinase M2 Increases Angiogenesis, Neurogenesis, and Functional Recovery Mediated by Upregulation of STAT3 and Focal Adhesion Kinase Activities after Ischemic Stroke in Adult Mice. Neurotherapeutics 15 (3), 770–784. doi:10.1007/s13311-018-0635-2
Chen, S.-D., Lee, J.-M., Yang, D.-I., Nassief, A., and Hsu, C. Y. (2002). Combination Therapy for Ischemic Stroke. Am. J. Cardiovasc. Drugs 2 (5), 303–313. doi:10.2165/00129784-200202050-00003
Chen, W., Wang, H., Feng, J., and Chen, L. (2020a). Overexpression of circRNA circUCK2 Attenuates Cell Apoptosis in Cerebral Ischemia-Reperfusion Injury via miR-125b-5p/GDF11 Signaling. Mol. Ther. - Nucleic Acids 22, 673–683. doi:10.1016/j.omtn.2020.09.032
Chen, Y., Wang, B., Liu, W., Xu, P., and Song, L. (2020b). Diagnostic Value of Serum Hsa_circ_0141720 in Patients with Acute Ischemic Stroke. Clin. Lab. 66 (8). doi:10.7754/Clin.Lab.2020.191266
Cocucci, E., Racchetti, G., and Meldolesi, J. (2009). Shedding Microvesicles: Artefacts No More. Trends Cel Biol. 19 (2), 43–51. doi:10.1016/j.tcb.2008.11.003
Couzin, J. (2005). The Ins and Outs of Exosomes. Science 308 (5730), 1862–1863. doi:10.1126/science.308.5730.1862
Dai, Q., Ma, Y., Xu, Z., Zhang, L., Yang, H., Liu, Q., et al. (2021). Downregulation of Circular RNA HECTD1 Induces Neuroprotection against Ischemic Stroke through the microRNA-133b/TRAF3 Pathway. Life Sci. 264, 118626. doi:10.1016/j.lfs.2020.118626
Diaz-Guerra, M. (2021). Excitotoxicity-induced Endocytosis as a Potential Target for Stroke Neuroprotection. Neural Regen. Res. 16 (2), 300–301. doi:10.4103/1673-5374.290892
Dong, Z., Deng, L., Peng, Q., Pan, J., and Wang, Y. (2020). CircRNA Expression Profiles and Function Prediction in Peripheral Blood Mononuclear Cells of Patients with Acute Ischemic Stroke. J. Cel Physiol 235 (3), 2609–2618. doi:10.1002/jcp.29165
Du, W. W., Yang, W., Chen, Y., Wu, Z.-K., Foster, F. S., Yang, Z., et al. (2017a). Foxo3 Circular RNA Promotes Cardiac Senescence by Modulating Multiple Factors Associated with Stress and Senescence Responses. Eur. Heart J. 38 (18), ehw001–1412. doi:10.1093/eurheartj/ehw001
Du, W. W., Zhang, C., Yang, W., Yong, T., Awan, F. M., and Yang, B. B. (2017b). Identifying and Characterizing circRNA-Protein Interaction. Theranostics 7 (17), 4183–4191. doi:10.7150/thno.21299
Fu, Y., Sun, C., Li, Q., Qian, F., Li, C., Xi, X., et al. (2021). Differential RNA Expression Profiles and Competing Endogenous RNA-Associated Regulatory Networks during the Progression of Atherosclerosis. Epigenomics 13 (2), 99–112. doi:10.2217/epi-2020-0252
Gisterå, A., and Hansson, G. K. (2017). The Immunology of Atherosclerosis. Nat. Rev. Nephrol. 13 (6), 368–380. doi:10.1038/nrneph.2017.51
Guan, R., Zou, W., Dai, X., Yu, X., Liu, H., Chen, Q., et al. (2018). Mitophagy, a Potential Therapeutic Target for Stroke. J. Biomed. Sci. 25 (1), 87. doi:10.1186/s12929-018-0487-4
Guzik, A., and Bushnell, C. (2017). Stroke Epidemiology and Risk Factor Management. CONTINUUM: Lifelong Learn. NeurologyCerebrovascular Disease) 23 (1), 15–39. doi:10.1212/con.0000000000000416
Han, B., Chao, J., and Yao, H. (2018). Circular RNA and its Mechanisms in Disease: From the Bench to the Clinic. Pharmacol. Ther. 187, 31–44. doi:10.1016/j.pharmthera.2018.01.010
Hansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., et al. (2013). Natural RNA Circles Function as Efficient microRNA Sponges. Nature 495 (7441), 384–388. doi:10.1038/nature11993
Hansson, G. K., and Libby, P. (2006). The Immune Response in Atherosclerosis: a Double-Edged Sword. Nat. Rev. Immunol. 6 (7), 508–519. doi:10.1038/nri1882
Heath-Engel, H. M., Chang, N. C., and Shore, G. C. (2008). The Endoplasmic Reticulum in Apoptosis and Autophagy: Role of the BCL-2 Protein Family. Oncogene 27 (50), 6419–6433. doi:10.1038/onc.2008.309
Hitchcock, I. S., Fox, N. E., Prévost, N., Sear, K., Shattil, S. J., and Kaushansky, K. (2008). Roles of Focal Adhesion Kinase (FAK) in Megakaryopoiesis and Platelet Function: Studies Using a Megakaryocyte Lineage-specific FAK Knockout. Blood 111 (2), 596–604. doi:10.1182/blood-2007-05-089680
Hon, K. W., Abu, N., Ab Mutalib, N.-S., and Jamal, R. (2017). Exosomes as Potential Biomarkers and Targeted Therapy in Colorectal Cancer: A Mini-Review. Front. Pharmacol. 8, 583. doi:10.3389/fphar.2017.00583
Hou, Y., Wang, K., Wan, W., Cheng, Y., Pu, X., and Ye, X. (2018). Resveratrol Provides Neuroprotection by Regulating the JAK2/STAT3/PI3K/AKT/mTOR Pathway after Stroke in Rats. Genes Dis. 5 (3), 245–255. doi:10.1016/j.gendis.2018.06.001
Hu, X., De Silva, T. M., Chen, J., and Faraci, F. M. (2017). Cerebral Vascular Disease and Neurovascular Injury in Ischemic Stroke. Circ. Res. 120 (3), 449–471. doi:10.1161/circresaha.116.308427
Hu, Y.-W., Zhao, J.-Y., Li, S.-F., Huang, J.-L., Qiu, Y.-R., Ma, X., et al. (2015). RP5-833A20.1/miR-382-5p/NFIA-dependent Signal Transduction Pathway Contributes to the Regulation of Cholesterol Homeostasis and Inflammatory Reaction. Atvb 35 (1), 87–101. doi:10.1161/atvbaha.114.304296
Huang, A., Patel, S., McAlpine, C., and Werstuck, G. (2018). The Role of Endoplasmic Reticulum Stress-Glycogen Synthase Kinase-3 Signaling in Atherogenesis. Ijms 19 (6), 1607. doi:10.3390/ijms19061607
Janas, T., Janas, M. M., Sapoń, K., and Janas, T. (2015). Mechanisms of RNA Loading into Exosomes. FEBS Lett. 589 (13), 1391–1398. doi:10.1016/j.febslet.2015.04.036
Janga, S. C., and Mittal, N. (2011). Construction, Structure and Dynamics of post-transcriptional Regulatory Network Directed by RNA-Binding Proteins. Adv. Exp. Med. Biol. 722, 103–117. doi:10.1007/978-1-4614-0332-6_7
Jia, D., Zhu, Q., Liu, H., Zuo, C., He, Y., Chen, G., et al. (2017). Osteoprotegerin Disruption Attenuates HySu-Induced Pulmonary Hypertension through Integrin α V β 3/FAK/AKT Pathway Suppression. Circ. Cardiovasc. Genet. 10 (1), 1591. doi:10.1161/circgenetics.116.001591
Jiang, Q., Su, D.-Y., Wang, Z.-Z., Liu, C., Sun, Y.-N., Cheng, H., et al. (2021). Retina as a Window to Cerebral Dysfunction Following Studies with circRNA Signature during Neurodegeneration. Theranostics 11 (4), 1814–1827. doi:10.7150/thno.51550
Jiang, S., Guo, C., Zhang, W., Che, W., Zhang, J., Zhuang, S., et al. (2019). The Integrative Regulatory Network of circRNA, microRNA, and mRNA in Atrial Fibrillation. Front. Genet. 10, 526. doi:10.3389/fgene.2019.00526
Jiang, S., Li, T., Ji, T., Yi, W., Yang, Z., Wang, S., et al. (2018). AMPK: Potential Therapeutic Target for Ischemic Stroke. Theranostics 8 (16), 4535–4551. doi:10.7150/thno.25674
Jiang, T., Yu, J.-T., Zhu, X.-C., Zhang, Q.-Q., Tan, M.-S., Cao, L., et al. (2015). Ischemic Preconditioning Provides Neuroprotection by Induction of AMP-Activated Protein Kinase-dependent Autophagy in a Rat Model of Ischemic Stroke. Mol. Neurobiol. 51 (1), 220–229. doi:10.1007/s12035-014-8725-6
Kernan, W. N., Ovbiagele, B., Black, H. R., Bravata, D. M., Chimowitz, M. I., Ezekowitz, M. D., et al. (2014). Guidelines for the Prevention of Stroke in Patients with Stroke and Transient Ischemic Attack. Stroke 45 (7), 2160–2236. doi:10.1161/str.0000000000000024
Lasda, E., and Parker, R. (2014). Circular RNAs: Diversity of Form and Function. Rna 20 (12), 1829–1842. doi:10.1261/rna.047126.114
Lefrançais, E., Ortiz-Muñoz, G., Caudrillier, A., Mallavia, B., Liu, F., Sayah, D. M., et al. (2017). The Lung Is a Site of Platelet Biogenesis and a Reservoir for Haematopoietic Progenitors. Nature 544 (7648), 105–109. doi:10.1038/nature21706
Li, L., Lou, W., Li, H., Zhu, Y., and Huang, X. a. (2020a). Upregulated C-C Motif Chemokine Ligand 2 Promotes Ischemic Stroke via Chemokine Signaling Pathway. Ann. Vasc. Surg. 68, 476–486. doi:10.1016/j.avsg.2020.04.047
Li, P., Chen, S., Chen, H., Mo, X., Li, T., Shao, Y., et al. (2015). Using Circular RNA as a Novel Type of Biomarker in the Screening of Gastric Cancer. Clinica Chim. Acta 444, 132–136. doi:10.1016/j.cca.2015.02.018
Li, S., Chen, L., Xu, C., Qu, X., Qin, Z., Gao, J., et al. (2020b). Expression Profile and Bioinformatics Analysis of Circular RNAs in Acute Ischemic Stroke in a South Chinese Han Population. Sci. Rep. 10 (1), 10138. doi:10.1038/s41598-020-66990-y
Li, S., Hu, W., Deng, F., Chen, S., Zhu, P., Wang, M., et al. (2021). Identification of Circular RNA Hsa_circ_0001599 as a Novel Biomarker for Large-Artery Atherosclerotic Stroke. DNA Cel Biol. 40 (3), 457–468. doi:10.1089/dna.2020.5662
Liu, G., Wu, F., Jiang, X., Que, Y., Qin, Z., Hu, P., et al. (2020). Inactivation of Cys 674 in SERCA2 Increases BP by Inducing Endoplasmic Reticulum Stress and Soluble Epoxide Hydrolase. Br. J. Pharmacol. 177 (8), 1793–1805. doi:10.1111/bph.14937
Liu, P., Xiang, J.-W., and Jin, S.-X. (2015). Serum CXCL12 Levels Are Associated with Stroke Severity and Lesion Volumes in Stroke Patients. Neurol. Res. 37 (10), 853–858. doi:10.1179/1743132815y.0000000063
Liu, W., Zhang, J., Zou, C., Xie, X., Wang, Y., Wang, B., et al. (2017a). Microarray Expression Profile and Functional Analysis of Circular RNAs in Osteosarcoma. Cell Physiol Biochem 43 (3), 969–985. doi:10.1159/000481650
Liu, X. S., Fan, B. Y., Pan, W. L., Li, C., Levin, A. M., Wang, X., et al. (2017b). Identification of miRNomes Associated with Adult Neurogenesis after Stroke Using Argonaute 2-based RNA Sequencing. RNA Biol. 14 (5), 488–499. doi:10.1080/15476286.2016.1196320
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262
Lu, D., Ho, E. S., Mai, H., Zang, J., Liu, Y., Li, Y., et al. (2020). Identification of Blood Circular RNAs as Potential Biomarkers for Acute Ischemic Stroke. Front. Neurosci. 14, 81. doi:10.3389/fnins.2020.00081
Lukiw, W. J. (2013). Circular RNA (circRNA) in Alzheimer's Disease (AD). Front. Genet. 4, 307. doi:10.3389/fgene.2013.00307
Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., et al. (2013). Circular RNAs Are a Large Class of Animal RNAs with Regulatory Potency. Nature 495 (7441), 333–338. doi:10.1038/nature11928
Mousavi, S., Moallem, R., Hassanian, S. M., Sadeghzade, M., Mardani, R., Ferns, G. A., et al. (2019). Tumor‐derived Exosomes: Potential Biomarkers and Therapeutic Target in the Treatment of Colorectal Cancer. J. Cel Physiol 234 (8), 12422–12432. doi:10.1002/jcp.28080
Murphy, J. M., Jeong, K., Rodriguez, Y. A. R., Kim, J.-H., Ahn, E.-Y. E., and Lim, S.-T. S. (2019). FAK and Pyk2 Activity Promote TNF-α and IL-1β-mediated Pro-inflammatory Gene Expression and Vascular Inflammation. Sci. Rep. 9 (1), 7617. doi:10.1038/s41598-019-44098-2
Naiel, S., Carlisle, R. E., Lu, C., Tat, V., and Dickhout, J. G. (2019). Endoplasmic Reticulum Stress Inhibition Blunts the Development of Essential Hypertension in the Spontaneously Hypertensive Rat. Am. J. Physiology-Heart Circulatory Physiol. 316 (5), H1214–h1223. doi:10.1152/ajpheart.00523.2018
O'Leary, V. B., Smida, J., Matjanovski, M., Brockhaus, C., Winkler, K., Moertl, S., et al. (2017). The circRNA Interactome-Innovative Hallmarks of the Intra- and Extracellular Radiation Response. Oncotarget 8 (45), 78397–78409. doi:10.18632/oncotarget.19228
Ostolaza, A., Blanco-Luquin, I., Urdánoz-Casado, A., Rubio, I., Labarga, A., Zandio, B., et al. (2020). Circular RNA Expression Profile in Blood According to Ischemic Stroke Etiology. Cell Biosci 10, 34. doi:10.1186/s13578-020-00394-3
Pamudurti, N. R., Bartok, O., Jens, M., Ashwal-Fluss, R., Stottmeister, C., Ruhe, L., et al. (2017). Translation of CircRNAs. Mol. Cel 66 (1), 9–21. doi:10.1016/j.molcel.2017.02.021
Peng, X., Jing, P., Chen, J., and Xu, L. (2019). The Role of Circular RNA HECTD1 Expression in Disease Risk, Disease Severity, Inflammation, and Recurrence of Acute Ischemic Stroke. J. Clin. Lab. Anal. 33 (7), e22954. doi:10.1002/jcla.22954
Powers, W. J., Rabinstein, A. A., Ackerson, T., Adeoye, O. M., Bambakidis, N. C., Becker, K., et al. (2018). 2018 Guidelines for the Early Management of Patients with Acute Ischemic Stroke: A Guideline for Healthcare Professionals from the American Heart Association/American Stroke Association. Stroke 49 (3), e46–e110. doi:10.1161/str.0000000000000158
Ragan, C., Goodall, G. J., Shirokikh, N. E., and Preiss, T. (2019). Insights into the Biogenesis and Potential Functions of Exonic Circular RNA. Sci. Rep. 9 (1), 2048. doi:10.1038/s41598-018-37037-0
Raposo, G., and Stoorvogel, W. (2013). Extracellular Vesicles: Exosomes, Microvesicles, and Friends. J. Cel Biol 200 (4), 373–383. doi:10.1083/jcb.201211138
Ruytinx, P., Proost, P., Van Damme, J., and Struyf, S. (2018). Chemokine-Induced Macrophage Polarization in Inflammatory Conditions. Front. Immunol. 9, 1930. doi:10.3389/fimmu.2018.01930
Rybak-Wolf, A., Stottmeister, C., Glažar, P., Jens, M., Pino, N., Giusti, S., et al. (2015). Circular RNAs in the Mammalian Brain Are Highly Abundant, Conserved, and Dynamically Expressed. Mol. Cel 58 (5), 870–885. doi:10.1016/j.molcel.2015.03.027
Schmalbach, B., Stepanow, O., Jochens, A., Riedel, C., Deuschl, G., and Kuhlenbäumer, G. (2015). Determinants of Platelet-Leukocyte Aggregation and Platelet Activation in Stroke. Cerebrovasc. Dis. 39 (3-4), 176–180. doi:10.1159/000375396
Sen, S., Tewari, M., Zajac, A., Barton, E., Sweeney, H. L., and Discher, D. E. (2011). Upregulation of Paxillin and Focal Adhesion Signaling Follows Dystroglycan Complex Deletions and Promotes a Hypertensive State of Differentiation. Eur. J. Cel Biol. 90 (2-3), 249–260. doi:10.1016/j.ejcb.2010.06.005
Sheng, R., Zhang, L.-S., Han, R., Liu, X.-Q., Gao, B., and Qin, Z.-H. (2010). Autophagy Activation Is Associated with Neuroprotection in a Rat Model of Focal Cerebral Ischemic Preconditioning. Autophagy 6 (4), 482–494. doi:10.4161/auto.6.4.11737
Simons, M., and Raposo, G. (2009). Exosomes - Vesicular Carriers for Intercellular Communication. Curr. Opin. Cel Biol. 21 (4), 575–581. doi:10.1016/j.ceb.2009.03.007
Skotland, T., Sandvig, K., and Llorente, A. (2017). Lipids in Exosomes: Current Knowledge and the Way Forward. Prog. Lipid Res. 66, 30–41. doi:10.1016/j.plipres.2017.03.001
Snow, S. J. (2016). Stroke and T-PA - Triggering New Paradigms of Care. N. Engl. J. Med. 374 (9), 809–811. doi:10.1056/NEJMp1514696
Sugimura, K., Fukumoto, Y., Nawata, J., Wang, H., Onoue, N., Tada, T., et al. (2010). Hypertension Promotes Phosphorylation of Focal Adhesion Kinase and Proline-Rich Tyrosine Kinase 2 in Rats: Implication for the Pathogenesis of Hypertensive Vascular Disease. Tohoku J. Exp. Med. 222 (3), 201–210. doi:10.1620/tjem.222.201
Sun, H., Zhong, D., Wang, C., Sun, Y., Zhao, J., and Li, G. (2018). MiR-298 Exacerbates Ischemia/Reperfusion Injury Following Ischemic Stroke by Targeting Act1. Cel Physiol Biochem 48 (2), 528–539. doi:10.1159/000491810
Sun, Y., Zhang, T., Li, L., and Wang, J. (2015). Induction of Apoptosis by Hypertension via Endoplasmic Reticulum Stress. Kidney Blood Press. Res. 40 (1), 41–51. doi:10.1159/000368481
Suzuki, H., Zuo, Y., Wang, J., Zhang, M. Q., Malhotra, A., and Mayeda, A. (2006). Characterization of RNase R-Digested Cellular RNA Source that Consists of Lariat and Circular RNAs from Pre-mRNA Splicing. Nucleic Acids Res. 34 (8), e63. doi:10.1093/nar/gkl151
Szegedi, I., Szapáry, L., Csécsei, P., Csanádi, Z., and Csiba, L. (2017). Potential Biological Markers of Atrial Fibrillation: A Chance to Prevent Cryptogenic Stroke. Biomed. Res. Int. 2017, 1–10. doi:10.1155/2017/8153024
Tabas, I. (2010). The Role of Endoplasmic Reticulum Stress in the Progression of Atherosclerosis. Circ. Res. 107 (7), 839–850. doi:10.1161/circresaha.110.224766
Tejeda, G. S., Esteban‐Ortega, G. M., San Antonio, E., Vidaurre, Ó. G., and Díaz‐Guerra, M. (2019). Prevention of Excitotoxicity‐induced Processing of BDNF Receptor TrkB‐FL Leads to Stroke Neuroprotection. EMBO Mol. Med. 11 (7), e9950. doi:10.15252/emmm.201809950
Troulinaki, K., and Tavernarakis, N. (2012). Endocytosis and Intracellular Trafficking Contribute to Necrotic Neurodegeneration inC. Elegans. EMBO J. 31 (3), 654–666. doi:10.1038/emboj.2011.447
von Essen, M., Rahikainen, R., Oksala, N., Raitoharju, E., Seppälä, I., Mennander, A., et al. (2016). Talin and Vinculin Are Downregulated in Atherosclerotic Plaque; Tampere Vascular Study. Atherosclerosis 255, 43–53. doi:10.1016/j.atherosclerosis.2016.10.031
Wang, M.-M., Zhang, M., Feng, Y.-S., Xing, Y., Tan, Z.-X., Li, W.-B., et al. (2020). Electroacupuncture Inhibits Neuronal Autophagy and Apoptosis via the PI3K/AKT Pathway Following Ischemic Stroke. Front. Cel. Neurosci. 14, 134. doi:10.3389/fncel.2020.00134
Wang, P., Guan, Y.-F., Du, H., Zhai, Q.-W., Su, D.-F., and Miao, C.-Y. (2012). Induction of Autophagy Contributes to the Neuroprotection of Nicotinamide Phosphoribosyltransferase in Cerebral Ischemia. Autophagy 8 (1), 77–87. doi:10.4161/auto.8.1.18274
Wang, P., Xu, T.-Y., Wei, K., Guan, Y.-F., Wang, X., Xu, H., et al. (2014). ARRB1/β-arrestin-1 Mediates Neuroprotection through Coordination of BECN1-dependent Autophagy in Cerebral Ischemia. Autophagy 10 (9), 1535–1548. doi:10.4161/auto.29203
Won Kim, H., Haider, H. K., Jiang, S., and Ashraf, M. (2009). Ischemic Preconditioning Augments Survival of Stem Cells via miR-210 Expression by Targeting Caspase-8-Associated Protein 2. J. Biol. Chem. 284 (48), 33161–33168. doi:10.1074/jbc.M109.020925
Wu, N., Jin, L., and Cai, J. (2017). Profiling and Bioinformatics Analyses Reveal Differential Circular RNA Expression in Hypertensive Patients. Clin. Exp. Hypertens. 39 (5), 454–459. doi:10.1080/10641963.2016.1273944
Wu, W., Xiao, H., Laguna-Fernandez, A., Villarreal, G., Wang, K.-C., Geary, G. G., et al. (2011). Flow-Dependent Regulation of Krüppel-like Factor 2 Is Mediated by MicroRNA-92a. Circulation 124 (5), 633–641. doi:10.1161/circulationaha.110.005108
Xiao, Q., Yin, R., Wang, Y., Yang, S., Ma, A., Pan, X., et al. (2021). Comprehensive Analysis of Peripheral Exosomal circRNAs in Large Artery Atherosclerotic Stroke. Front. Cel Dev. Biol. 9, 685741. doi:10.3389/fcell.2021.685741
Yan, Y., Thakur, M., van der Vorst, E. P. C., Weber, C., and Döring, Y. (2021). Targeting the Chemokine Network in Atherosclerosis. Atherosclerosis 330, 95–106. doi:10.1016/j.atherosclerosis.2021.06.912
Yang, C., Hawkins, K. E., Doré, S., and Candelario-Jalil, E. (2019). Neuroinflammatory Mechanisms of Blood-Brain Barrier Damage in Ischemic Stroke. Am. J. Physiology-Cell Physiol. 316 (2), C135–c153. doi:10.1152/ajpcell.00136.2018
Yang, Y., Fan, X., Mao, M., Song, X., Wu, P., Zhang, Y., et al. (2017). Extensive Translation of Circular RNAs Driven by N6-Methyladenosine. Cell Res 27 (5), 626–641. doi:10.1038/cr.2017.31
Yang, Y., Gao, X., Zhang, M., Yan, S., Sun, C., Xiao, F., et al. (2018). Novel Role of FBXW7 Circular RNA in Repressing Glioma Tumorigenesis. J. Natl. Cancer Inst. 110 (3), 304–315. doi:10.1093/jnci/djx166
Yao, J.-T., Zhao, S.-H., Liu, Q.-P., Lv, M.-Q., Zhou, D.-X., Liao, Z.-J., et al. (2017). Over-expression of CircRNA_100876 in Non-small Cell Lung Cancer and its Prognostic Value. Pathol. - Res. Pract. 213 (5), 453–456. doi:10.1016/j.prp.2017.02.011
Yoo, W., Zieba, J. K., Foegeding, N. J., Torres, T. P., Shelton, C. D., Shealy, N. G., et al. (2021). High-fat Diet-Induced Colonocyte Dysfunction Escalates Microbiota-Derived Trimethylamine N -oxide. Science 373 (6556), 813–818. doi:10.1126/science.aba3683
Yorimitsu, T., Nair, U., Yang, Z., and Klionsky, D. J. (2006). Endoplasmic Reticulum Stress Triggers Autophagy. J. Biol. Chem. 281 (40), 30299–30304. doi:10.1074/jbc.M607007200
Zhang, H., Deng, T., Liu, R., Bai, M., Zhou, L., Wang, X., et al. (2017). Exosome-delivered EGFR Regulates Liver Microenvironment to Promote Gastric Cancer Liver Metastasis. Nat. Commun. 8, 15016. doi:10.1038/ncomms15016
Zhang, M., Zhao, K., Xu, X., Yang, Y., Yan, S., Wei, P., et al. (2018). A Peptide Encoded by Circular Form of LINC-PINT Suppresses Oncogenic Transcriptional Elongation in Glioblastoma. Nat. Commun. 9 (1), 4475. doi:10.1038/s41467-018-06862-2
Zhao, W., Liang, J., Chen, Z., Diao, Y., and Miao, G. (2020). Combined Analysis of circRNA and mRNA Profiles and Interactions in Patients with Diabetic Foot and Diabetes Mellitus. Int. Wound J. 17 (2), 1183–1193. doi:10.1111/iwj.13420
Zheng, L., Cheng, W., Wang, X., Yang, Z., Zhou, X., and Pan, C. (2017). Overexpression of MicroRNA-145 Ameliorates Astrocyte Injury by Targeting Aquaporin 4 in Cerebral Ischemic Stroke. Biomed. Res. Int. 2017, 1–9. doi:10.1155/2017/9530951
Zheng, Q., Bao, C., Guo, W., Li, S., Chen, J., Chen, B., et al. (2016). Circular RNA Profiling Reveals an Abundant circHIPK3 that Regulates Cell Growth by Sponging Multiple miRNAs. Nat. Commun. 7, 11215. doi:10.1038/ncomms11215
Keywords: circRNA, plasma exosome, profile, acute ischemic stroke, ceRNA
Citation: Yang J, Hao J, Lin Y, Guo Y, Liao K, Yang M, Cheng H, Yang M and Chen K (2022) Profile and Functional Prediction of Plasma Exosome-Derived CircRNAs From Acute Ischemic Stroke Patients. Front. Genet. 13:810974. doi: 10.3389/fgene.2022.810974
Received: 08 November 2021; Accepted: 31 January 2022;
Published: 14 March 2022.
Edited by:
Gabriel Adelman Cipolla, Federal University of Paraná, BrazilReviewed by:
Jiankun Zang, Jinan University, ChinaVenkata Garikipati, The Ohio State University, United States
Copyright © 2022 Yang, Hao, Lin, Guo, Liao, Yang, Cheng, Yang and Chen. 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.
*Correspondence: Kejie Chen, Y2tqOTMwQDEyNi5jb20=