Construction of the circRNA-miRNA-mRNA Regulatory Network of an Abdominal Aortic Aneurysm to Explore Its Potential Pathogenesis

Zhang, Hao; Bian, Ce; Tu, Simei; Yin, Fanxing; Guo, Panpan; Zhang, Jian; Wu, Yihao; Yin, Yuhan; Guo, Jiahui; Han, Yanshuo

doi:https://doi.org/10.1155/2021/9916881

Disease Markers

On this page

Abstract Introduction Materials and Methods Results Discussion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Circulating Non-Coding RNAs as Promising Candidates for Biomarker-Guided Therapeutics

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 9916881 | https://doi.org/10.1155/2021/9916881

Construction of the circRNA-miRNA-mRNA Regulatory Network of an Abdominal Aortic Aneurysm to Explore Its Potential Pathogenesis

Hao Zhang,¹Ce Bian,²Simei Tu,¹Fanxing Yin,¹Panpan Guo,¹Jian Zhang,³Yihao Wu,¹Yuhan Yin,¹Jiahui Guo,¹and Yanshuo Han¹

Academic Editor: Ayesha Obaid

Received19 May 2021

Revised29 Sept 2021

Accepted20 Oct 2021

Published05 Nov 2021

Abstract

Background. Abdominal aortic aneurysm (AAA) is a progressive cardiovascular disease, which is a permanent and localized dilatation of the abdominal aorta with potentially fatal consequence of aortic rupture. Dysregulation of circRNAs is correlated with the development of various pathological events in cardiovascular diseases. However, the function of circRNAs in abdominal aortic aneurysm (AAA) is unknown and remains to be explored. This study is aimed at determining the regulatory mechanisms of circRNAs in AAAs. This study was aimed at exploring the underlying molecular mechanisms of abdominal aortic aneurysms based on the competing endogenous RNA (ceRNA) regulatory hypothesis of circRNA, miRNA, and mRNA. Methods. The expression profiles of circRNAs (GSE144431), miRNAs (GSE62179), and mRNAs (GSE7084, GSE57691, and GSE47472) in human tissue sample from the aneurysm group and normal group were obtained from the Gene Expression Omnibus database, respectively. The circRNA-miRNA-mRNA network was constructed by using Cytoscape 3.7.2 software; then, the protein-protein interaction (PPI) network was constructed by using the STRING database, and the hub genes were identified by using the cytoHubba plug-in. The circRNA-miRNA-hub gene regulatory subnetwork was formed to understand the regulatory axis of hub genes in AAAs. Results. The present study identified 40 differentially expressed circRNAs (DECs) in the GSE144431, 90 differentially expressed miRNAs (DEmiRs) in the GSE62179, and 168 differentially expressed mRNAs (DEGs) with the same direction regulation (130 downregulated and 38 upregulated) in the GSE7084, GSE57691, and GSE47472 datasets identified regarding AAAs. The miRNA response elements (MREs) of three DECs were then predicted. Four overlapping miRNAs were obtained by intersecting the predicted miRNA and DEmiRs. Then, 17 overlapping mRNAs were obtained by intersecting the predicted target mRNAs of 4 miRNAs with 168 DEGs. Furthermore, the circRNA-miRNA-mRNA network was constructed through 3 circRNAs, 4 miRNAs, and 17 mRNAs, and three hub genes (SOD2, CCR7, and PGRMC1) were identified. Simultaneously, functional enrichment and pathway analysis were performed within genes in the circRNA-miRNA-mRNA network. Three of them (SOD2, CCR7, and PGRMC1) were suggested to be crucial based on functional enrichment, protein-protein interaction, and ceRNA network analysis. Furthermore, the expression of SOD2 and CCR7 may be regulated by hsa_circ_0011449/hsa_circ_0081968/hsa-let-7f-5p; the expression of PGRMC1 may be regulated by hsa_circ_0011449/hsa_circ_0081968-hsa-let-7f-5p/hsa-let-7e-5p. Conclusion. In conclusion, the ceRNA interaction axis we identified may be an important target for the treatment of abdominal aortic aneurysms. This study provided further understanding of the potential pathogenesis from the perspective of the circRNA-related competitive endogenous RNA network in AAAs.

1. Introduction

Abdominal aortic aneurysm (AAA) is one of the cardiovascular diseases mainly distributed in Western countries and Asia, which mostly occurs in elderly men [1–3]. Local permanent dilation and weakening of the abdominal aorta are important characteristics of AAA [4]. Typically, degeneration of the abdominal aorta is a manifestation of a systemic process characterized by inflammation, apoptosis of smooth muscle cells, and destruction of elastin and collagen in the media and adventitia [5, 6]. More specifically, it would arise three clinical symptoms, asymptomatic and ruptured signs, and the characteristics of permanent expansion [7]. Generally, the mortality rate of AAA is up to 85%, which is inseparable from the severe bleeding and unpredictability caused by its rupture, and there is no effective treatment at present [8]. Nevertheless, AAA mortality data based on 1990-2010 from the Global Burden of Disease Study showed that the mortality of AAA decreased by 12.7% [9], and 43 percent of AAAs patients who were not included in the AAA screening criteria died, of which 9% were men and 34% were women [10]. Most of AAA patients were associated with asymptomatic AAA; therefore, the majority of AAA are identified incidentally during an examination for another unrelated pathology [11]. Therefore, the method to prevent the formation of aneurysms and the slow progress of aneurysms has become a growing public concern [12, 13].

Recently, several evidences have demonstrated that noncoding RNAs, such as long noncoding RNA (lncRNA), circular RNA (circRNA), and microRNA (miRNA), play an important role in the development of cardiovascular diseases [14–16]. CircRNA is spliced by exons at specific splicing sites to form a circular closed structure [17]. CircRNAs act as a competitive endogenous RNA (ceRNA) and compete with miRNA through its miRNA response element. miRNAs negatively regulate mRNA expression of protein-coding genes by binding to complementary sequences [18, 19]. Therefore, the circRNA-miRNA-mRNA interaction may be a significant mechanism in the occurrence and development of AAAs. According to our previous study, posttranscriptional modifications play important roles in modulating the functions of RNA involved in abdominal aorta aneurysm [20]. Yue et al. revealed the role of circCBFB/miR-28-5p/GRIA4/LYPD3 in the apoptosis of vascular smooth muscle cells (VSMCs), which conduced to the generation of new thoughts in AAA management [21]. Recently, Yang et al. investigated the role of circRNA CCDC66 in the pathogenesis of AAAs and determine its targeted miRNA and mRNA. The results show that CCDC66 will affect the proliferation of VSMCs through miR-342-3p/CCDC6 through circRNA CCDC66. They suggested that the overexpression of CCDC66 would induce the penetration of AAAs through the circCCDC66/miR-342-3p/CCDC66 pathway [22]. Meanwhile, Zhao et al. demonstrated that the pathogenesis of AAAs can be regulated by the CDR1as/miR-7/CKAP4 axis and affect the proliferation and apoptosis of primary vascular smooth muscle cells [23]. In summary, the circRNA-miRNA-mRNA network might play an important role in the pathogenesis of AAA, and circRNA has different targets and functions in different tissue cells.

Recently, Tian et al. performed microarray sequencing on human and rat abdominal aortic aneurysmal tissue and determined the differentially expressed lncRNA, miRNA, and mRNA by analyzing the public dataset and constructed the ceRNA network of lncRNA [24]. However, the specific targets and mechanisms of circRNA in AAAs have not been reported. Meanwhile, there are few reports on the ceRNA regulation mechanism of circRNA-miRNA-mRNA related to AAAs and on the interaction between ncRNA.

In this study, the microarray data of circRNA, miRNA, and mRNA were downloaded from the Gene Expression Comprehensive (GEO) database, and RStudio was used to perform differential expression analysis of the microarray data to obtain differentially expressed circRNAs (DECs), differentially expressed miRNAs (DEmiRs), and differentially expressed genes (DEGs). Cytoscape 3.7.2 constructs the circRNA-miRNA-mRNA network, and the STRING (search tool for interactive gene retrieval) online analysis tool constructs the PPI network. This study is aimed at exploring novel circRNAs and their mechanisms in human tissues of AAA patients. The purpose of this study is to further screen the ceRNA axis of the key circRNA-miRNA-mRNA in AAAs with microarray data collected from public databases and bioinformatics methods. Besides, the findings of this study promoted the understanding of the molecular mechanism of AAAs, which provided a new direction for targeted therapy of AAAs. The circRNA-miRNA-mRNA regulatory subnetwork was formed to understand the regulatory axis of hub genes in AAA. Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) are introduced to reveal the potential impact of key genes in the development of AAA. [25]. The flowchart is shown in Figure 1.

2. Materials and Methods

2.1. Data Extraction

Microarray data of AAAs were downloaded from the National Center for Biotechnology Information (NCBI) GEO database (https://www.ncbi.nlm.nih.gov/geo/), with keywords searched such as “Abdominal aortic aneurysm” (AAA) and “Noncoding RNA” (circRNA, miRNA, and mRNA). During the retrieving process, “Homo sapiens” was set as the filtering condition; as a result, a total of five datasets were included in this study, namely, GSE144431, GSE62179, GSE7084, GSE47472, and GSE57691. As for circRNA expression data (GSE144431 dataset), the type of data was ncRNA profiling by array and the data of the microarray platform was 074301 Arraystar Human ncRNA microarray V2 (platform: GPL21825). In terms of the GSE62179 dataset, the type of data was miRNA profiling by array and the data of the microarray platform was Agilent-021827 Human miRNA Microarray G4470C (Feature Number version) (platform: GPL14767). Similarly, as for the three datasets (GSE7084, GSE47472, and GSE57691), the type of data was expression profiling by array and the data of the microarray platform was HG-U133_Plus_2 Affymetrix Human Genome U133 Plus 2.0 Array (platform: GPL570), Illumina HumanHT-12 V4.0 expression beadchip (platform: GPL10558), and Illumina HumanHT-12 V4.0 expression beadchip (platform: GPL10558), respectively. Each dataset was composed of AAA samples and healthy aorta samples. Series matrix files for differential expression analysis and other function analysis were downloaded from the Gene Expression Omnibus repository database. The box plot of datasets is shown in Supplementary Figure 1.

2.2. Differential Expression Analysis

GEO2R is a GEO database online analysis tool that can directly perform differential expression analysis. We directly exported the analysis results to make differential analysis. GEO2R and R package were employed to identify the differential expression of genes, and the corrected data matrix of each dataset was downloaded directly from the GEO database, with the differentially expressed gene (DEG) analysis method using the “Linear Models for Microarray Data (limma)” R package for the datasets. The annotation information of DEGs, differentially expressed circRNAs (DECs), and differentially expressed miRNAs (DEmiRs) of the three datasets were obtained with the GEO2R online software. The selection criteria for DECs included , and value < 0.1 was considered to indicate a statistically significant difference. The selection criteria for DEGs included , and value < 0.05 was considered to indicate a statistically significant difference. There was statistical significance for the selection of this threshold, and those genes that were up- and downregulated can also be selected for the subsequent analysis. The R software package was also adopted to draw a heat map of differentially expressed circRNAs and differentially expressed mRNAs.

In this study, these differentially expressed genes in the three mRNA databases were divided into the upregulated and downregulated groups. Subsequently, the Venn diagram was employed to obtain mRNAs that were regulated in the same regulated pattern, and Venny (version 2.1.0) (https://bioinfogp.cnb.csic.es/tools/venny/), the online mapping software, was utilized to create the Venn diagram. The R software package was used for functional enrichment and pathway analysis of differential genes. Due to the different types of miRNA data and different miRNAs in AAA groups were extracted compared with normal groups with the aim of exploring DEmiRs in different aneurysm cells, and the network was constructed at the end.

2.3. Protein-Protein Interaction Network Analysis of DEGs

STRING protein interaction prediction (version110.b) (STRING: functional protein association networks (string-db.org)) was used to search for known protein interactions and predict protein interactions. The intersections of the differentially expressed mRNA target genes were taken in a consistent manner from three mRNA datasets and were imported into the database. After the points without interaction were hidden, the data were imported into Cytoscape (version 3.7.2) to visualize the protein-protein interaction (PPI) network.

2.4. Prediction of circRNA-miRNA Pairs

The interaction relationships between circRNAs and miRNA were predicted with circBank (http://www.circbank.cn/help.html) [26]. circBank was employed to predict circRNA. Overlapping miRNAs were obtained by intersecting the predicted circRNAs and DEmiRs.

2.5. Prediction of Target Genes of miRNAs

The miRNA-mRNA interactions were based on the following three miRNA target prediction databases: miRTarBase (https://maayanlab.cloud/Harmonizome/resource/MiRTarBase), TargetScan (http://www.targetscan.org/vert_72/), and miRDB (http://mirdb.org/). Perl (version 5.32.0) was adopted to analyze DEmiRs and the targeted mRNAs. Overlapping mRNAs were obtained by intersecting the predicted mRNAs and DEmiRs.

2.6. CircRNA-miRNA-mRNA ceRNA Network Construction

In this study, the intersection of the three databases was selected as the results. The ceRNA network was constructed based on the commonly interactive differentially expressed mRNAs with differentially expressed miRNAs and differentially expressed circRNAs. The circRNA-miRNA-mRNA ceRNA network was visualized with Cytoscape (version 3.7.2) (https://cytoscape.org/).

2.7. Functional Enrichment and Pathway Analysis

To explore the main functions and pathways of DEGs derived from differentially expressed circRNA target genes, differentially expressed mRNA target genes in a consistent manner, and target genes of the circRNA network, the Database for Annotation, Visualization and Integrated Discovery online tool (version 6.8) (DAVID: https://david.ncifcrf.gov/) was introduced for the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. The gene lists were downloaded for the biological function information.

2.8. Differential Expression Analysis of Key Genes

The plug-in cytoHubba in the Cytoscape software is visualization software that obtains the dense relationship through the degree, betweenness centrality, and closeness centrality algorithms. To obtain the key genes in the PPI network, we identified the hub genes in the PPI network through cytoHubba. Key genes were obtained through the intersection of Venn diagram PPI network hub genes and ceRNA network genes. To accurately explore the differential expression of key genes in the three datasets, we used the R language software package to draw box plots of key genes.

3. Results

3.1. Differential Expression Analysis

According to the threshold criteria ( and ), the GSE144431 datasets contained circRNA data, and a total of 382 differentially expressed circRNAs were visualized by hierarchical clustering heat map analysis (Figure 2(a)) and volcano plot analysis (Figure 3(a)), with a view to endowing the selection results with much significance. A total of 40 samples ranking for the top 40 with the most significant differences were selected, including 20 upregulated DECs and 20 downregulated DECs (Table 1).

(a)

(b)

(c)

(d)

Figure 2

Hierarchical clustering and heat map analysis of differentially expressed circRNAs (a) and mRNAs (b–d). (a) GSE144431, circRNA; (b) GSE7084, mRNA; (c) GSE57691, mRNA; (d) GSE47472, mRNA. The green color indicates low expression, the red indicates high expression. Column: an aortic sample. Differentially expressed circRNA molecules were screened under the cut-off criteria and the value (). Differentially expressed mRNA molecules were screened under the cut-off criteria and the value ().

(a)

(b)

(c)

(d)

The mRNA datasets of GSE7084, GSE57691, and GSE47472 contained mRNA microarray data. All differentially expressed mRNAs in GSE7084, GSE57691, and GSE47472 were selected for hierarchical clustering heat map analysis and visualization (Figures 2(b)–2(d)); volcano plots are shown in Figures 3(b)–3(d). Simultaneously, the common mRNAs were extracted from three differentially expressed datasets for the first time, and a total of 215 common DEGs were obtained at the end. Among all results, the top 20 items with the most significant values were selected and included in Table 1. As for 215 common DEGs, it could be inferred from that Protein Phosphatase 1 Catalytic Subunit Beta (PPP1CB) and Transmembrane Protein 47 (TMEM47) were upregulated in the GSE47472 and downregulated in the GSE7084 and GSE57691. Similarly, G Protein-Coupled Receptor Class C Group 5 Member C (GPRC5C) was upregulated in the GSE57691 and downregulated in the GSE7084 and GSE47472 (Table 1). Therefore, the adverse regulated DEG results could not be regarded as the content of further analysis. In order to make the research results obtained more accurate, those DEGs with a different regulated trend in the three microarray databases were excluded. Subsequently, the DEGs with the same regulated trend in the three databases were intersected and 168 DEGs were selected. Among them, 130 DEGs were downregulated (Figure 4(a)) and 38 DEGs upregulated (Figure 4(b)). The GO analysis of 168 DEGs involves the biological process and cellular component, and biological process included vasoconstriction, negative regulation of blood vessel diameter, and cell-matrix adhesion (Figure 4(c)). The KEGG analysis of 168 differentially expressed genes included the calcium signaling pathway, vascular smooth muscle contraction, and focal adhesion (Figure 4(d)); 47 common DEGs of different regulation methods in the three datasets are listed in Supplement Table 1.

(a)

(b)

(c)

(d)

The GSE62179 dataset contained miRNA data, and there were 90 differentially expressed miRNAs in the control group compared with the AAA groups of M1 macrophages, M2 macrophages, and smooth muscle cells (SMCs) (Table 2).

3.2. PPI Network for DEGs

When the consistent differentially expressed mRNAs in the GSE7084, GSE47472, and GSE57691 were imported to the STRING website, it could be found that there were 131 interaction relationships for 168 DEGs, which were adopted to construct the PPI network. After the points that are not connected and interacted were hidden, the PPI network file was downloaded. In the present network, AHR could interact with MYOC, PTGS2, and CXCL8; CCR7 could interact with PRKCDBP, IL2RB, ITK, IL21R, PTGS2, and CXCL8; SOD2 could interact with CXCL8, DNM1L, ESD, MRPL33, PTGS2, and TAGLN (Figure 5). Furthermore, cytoHubba was used to identify hub genes with high accuracy according to the network topology. After the filtering condition was set according to the maximal clique centrality (Mcc), the determined top 30 hub genes could be obtained, including IL2RB, ITK, IL21R, MYLK, PTPN22, CXCL8, PTGS2, ITGB5, CD3D, AGTR1, CLEC12A, SORBS1, CCR7, STOM, CALD1, PDGFA, PGRMC1, RAB37, MGAM, SOD2 and etc, which may indicate potential regulatory networks and regulatory pathways.

3.3. Prediction of circRNA-miRNA Pairs

The circRNA-miRNA pairs corresponding to 40 DECs were predicted by using the Circular RNA Interactome online software (circBank), which is used for online prediction of miRNA bound to circRNA by combining with circBank. The 2230 miRNA interaction pairs combined with circRNA were predicted by using circBank; meanwhile, only the last five results of miRNA bound to circRNA were shown here (Supplementary Table 2). The predicted results were intersected with the differentially expressed miRNAs of the dataset, and 22 miRNAs could be obtained.

3.4. Prediction of Target Genes of miRNAs

Moreover, there were 957 mRNAs predicted to bind to 22 miRNAs by compilations of the miRTarBase database, TargetScan database, and miRDB database. The results obtained were intersected with 168 DEGs, and as a result, there were 17 potential mRNAs included.

3.5. Construction of the circRNA-miRNA-mRNA Network

After the points where circRNA-miRNA-mRNA did not connect and interact were removed, circRNA-miRNA-mRNA regulatory networks were constructed. There were 3 DECs, 4 DEmiRs, and 17 DEGs in the ceRNA network (Figure 6), and KIAA0930 and PGRMC1 can be targeted by two miRNAs (hsa-let-7f-5p and hsa-let-7e-5p). These two DEmiRs could also interact with hsa_circ_001149 and hsa_circ_0081968, which would form the following ceRNA relationships, namely, hsa_circ_0081968/hsa_circ_0011449-hsa-let-7f-5p/hsa-let-7e-5p-KIAA0930/PGRMC1. Therefore, hsa_circ_005073/hsa_circ_0081968-hsa-miR-107-FGFRL1/PLEKHF2 and hsa_circ_0081968/hsa_circ_0011449-hsa-let-7f-5p-SOD2 ceRNA networks also could be obtained.

Eventually, the circRNA-miRNA-hub gene subnetwork with six regulatory modules, including the hsa_circ_0011449/hsa-let-7f-5p/SOD2 regulatory axis, hsa_circ_0081968/hsa-let-7f-5p/SOD2 regulatory axis, hsa_circ_0011449/hsa-let-7f-5p/CCR7 regulatory axis, hsa_circ_0081968/hsa-let-7f-5p/CCR7 regulatory axis, hsa_circ_0081968/hsa-let-7e-5p/PGRMC1 regulatory axis, and hsa_circ_0081968/hsa-let-7f-5p/PGRMC1 regulatory axis, was constructed to depict the relationship among circRNAs, miRNAs, and hub genes.

3.6. Functional Enrichment and Pathway Analysis

DECs included 20 differentially expressed circRNAs, while DEGs included 168 consistent differentially expressed mRNAs. GO and KEGG analyses were applied on DECs and DEGs, respectively. These DECs were analyzed through the three processes of Gene Ontology (GO) analysis, including biological processes (BP), cellular components (CC), and molecular functions (MF). KEGG analysis was applied to find the functional pathway of those genes, with the findings showing that DEGs in 20 differentially expressed circRNAs were enriched into eight GO biological process terms, e.g., nucleoplasm and actin cytoskeleton and 10 KEGG pathways such as the activation of Rac and Beta2 integrin cell surface interactions (Table 3).

Due to the fact that consistent overlapping DEGs were collected from three mRNA microarray datasets, they were divided into two parts for GO and KEGG analysis (the downregulation group and the upregulation group). In the downregulation group, there were 23 GO biological process terms and 5 KEGG pathways, e.g., positive regulation of cell migration (MYOC, NTRK3, PDGFA, ILK, LGR6, and MYLK), caveola assembly (CAV2 and PACSIN2), cGMP-PKG signaling pathway (EDNRA, PLN, AGTR1, ATP1A2, and MYLK), and vascular smooth muscle contraction (EDNRA, ARHGEF12, CALD1, AGTR1, PPP1R12B, and MYLK) (Table 4). In the upregulation group, there were 14 GO biological process terms and 3 KEGG pathways, e.g., T cell receptor signaling pathway (ITK, PDE4B, PTPN22, CD3, and PAG1), negative regulation of cell proliferation (CXCL8, RASSF5, SOD2, and PTGS2), cytokine-cytokine receptor interaction (CXCL8, IL2RB, IL21R, CCR7, and IL12RB1), and chemokine signaling pathway (ITK, CXCL8, and CCR7) (Supplementary Table 3).

To reveal the biological mechanism or functional pathway of the ceRNA network, DEGs in the ceRNA network were also subjected to GO/KEGG analysis, with the findings of functional enrichment showing that there were six GO biological process terms and two KEGG pathways, e.g., liver development (AK4 and SOD2), COP9 signalosome (NCKIPSD and COPS8), and transport vesicle (PLEKHF2 and FGFRL1) for GO analysis and peroxisome (SOD2) and FoxO signaling pathway (SOD2) for the KEGG pathway (Table 5).

3.7. Differential Expression Analysis of Key Genes

The top 20 genes obtained by the degree, betweenness centrality, and closeness centrality algorithms in cytoHubba are listed in Figure 7(a); overlapping genes of the PPI network and ceRNA network were selected as hub genes SOD2, CCR7, and PGRMC1 (Figure 7(b)); the differential expression analysis of the hub genes in the three datasets is shown in Figures 7(c)–7(e). Compared with the control group, SOD2 was upregulated in the three datasets, CCR7 was upregulated in the three datasets, and PGRMC1 was downregulated in the three datasets.

(a)

(b)

(c)

(d)

(e)

Figure 7

Differential expression analysis of hub genes from PPI network and ceRNA network genes. (a) Hub gene relationship networks from the PPI network. (b) Venn diagram analysis of hub genes from PPI network and ceRNA network genes. (c) Box plot of the differentially expressed gene CCR7 in the three datasets. (d) Box plot of the differentially expressed gene PGRMC1 in the three datasets. (e) Box plot of the differentially expressed gene SOD2 in the three datasets. Green: normal samples; purple: aneurysm samples; : number of samples.

4. Discussion

In this study on AAAs, it has been confirmed that noncoding RNAs exert significant impacts on the pathogenesis of the disease [27, 28]. CircRNA is a stable noncoding RNA that has long been neglected by transcriptomics due to the lack of a 5 cap and a 3 polyadenylated tail. As recent studies, there are many endogenous circRNAs in mammalian cells, some of which show high abundances and evolutionary conservation. It has been initially indicated that circRNAs could mediate miRNA functions (e.g., via sponging) and control important events in transcription (e.g., RNA folding and endonuclease protection) [19]. In addition, it has been demonstrated in current evidence that circRNAs contain multiple microRNA response elements (MREs) that can bind to miRNAs, commonly called “miRNA sponges,” which could relieve the targeted inhibition of downstream mRNAs by miRNAs [29–31]. In this study, the underlying molecular mechanisms of AAAs have been investigated based on ceRNA regulatory hypothesis of circRNAs, miRNA, and mRNA. The expression profiles of circRNAs, miRNAs, and mRNAs in human tissue samples from the aneurysm group and the normal group have been identified.

In the three RNA (circRNA, miRNA, and mRNA) microarray datasets, 17 DEGs have been identified in the circRNA-miRNA-mRNA network, and they form a ceRNA regulatory network. Hub genes, as key genes, could exert decisive impacts on the biological processes. The regulation of other genes in the relevant pathway would be frequently affected by hub genes. Therefore, it is of great significance to find the hub genes for potential targets or research hotspots. CytoHubba is effective software of the Cytoscape plug-in, which is used to identify hub genes with high accuracy according to the network topology. Within bioinformatics analysis, three hub genes (SOD2, CCR7, and PGRMC1) were identified.

Strauss et al. suggested that SOD2 may be related to the formation of abdominal aortic aneurysm through oxidation reaction [32]. However, there are few reports on the effect of SOD2 genes in AAAs. Superoxide dismutase (SOD) is a ubiquitous antioxidant enzyme that catalytically converts the superoxide radical to hydrogen peroxide (H₂O₂) [33, 34]. The copper/zinc- (SOD1) and manganese- (MnSOD, SOD2) requiring superoxide dismutases are important and major antioxidant enzymes located in the mitochondria and exert decisive impacts on the progression of tumors [35, 36]. For instance, epigenetic silencing of SOD2 in KAS 6/1 human multiple myeloma cells can increase cell proliferation [37]. The blocking of SOD2 expression would significantly inhibit TNF-α-induced cell proliferation in A549 and H1299 cells in vitro. Therefore, TNF-α-mediated lung inflammation can upregulate SOD2 expression in lung adenocarcinoma, and macrophages contribute to SOD2 upregulation by secreting TNF-α [38]. In cardiovascular disease, SOD2 stands at the forefront against mitochondrial ROS in vascular smooth muscle cells (VSMCs) via its preferential localization to the mitochondria, which conduces to the potential modification of the probability of vascular calcification [39] initiation or progression. Moreover, SOD2 augmentation may be a promising therapeutic strategy for the prevention of lesion formation in proliferative vascular diseases such as restenosis [40]. Madamanchi et al. suggested that SOD2 could regulate SMC quiescence by suppressing divergent mitogenic signaling pathways, and the dysregulation of these enzymes under pathophysiological conditions may induce SMC hyperplasia and hypertrophy [41]. The expression level of SOD2 will also have an influence on the migration of VSMCs and the formation of intima. Mechanically, Qu et al. suggested that peroxisome proliferator-activated receptor-gamma coactivator-1alpha (PGC-1alpha) might upregulate the expression of the mitochondrial antioxidant enzyme SOD2, which would inhibit VSMC migration and neointimal formation after vascular injury in rats [42]. On the contrary, it has been reported in other studies that the loss of heterozygous SOD2 would aggravate chronic intermittent hypoxia-induced lung inflammation and vascular remodeling through the mtROS-NLRP3 signaling pathway [43]. Finally, the increase in SOD2 expression in genotyped primary human umbilical vein endothelial cells (ECs) induced by fluid shear stress (FSS) would effectively stabilize their antiatherosclerotic phenotype [44]. Besides, it has also been found in this study that SOD2 may exert its impacts through hsa_circ_0011449/hsa_circ_0081968/hsa-let-7f-5p. To sum up, it has been indicated in those studies that the dysregulation of SOD2 can cause related arterial diseases or changes in the endothelial cell or smooth muscle cell function and phenotype. Therefore, the research on SOD2 and its functions may become a popular trend for the treatment of AAAs. Considering the important role that SOD2 may play in AAAs, some papers have been found to explain the effects of SOD2 on arteries.

At the same time, it is also the most significant among the differentially expressed genes in AAAs. hsa-let-7f-5p can form a regulatory network with many differentially expressed genes, for example, interaction with AHR (Aryl Hydrocarbon Receptor). Besides, Fibroblast Growth Factor Receptor-Like 1 (FGFRL1) and Pleckstrin Homology and FYVE Domain Containing 2 (PLEKHF2) may affect the transportation of intracellular substances to promote aneurysm formation, through the hsa_circ_0005073/hsa_circ_0081968/hsa-miR-107 regulated axis. There are relatively few studies on PLEKHF2 and FGFRL1 to indicate that they may exert significant impacts on the occurrence and development of AAAs. E. Shamsara and J. Shamsara have confirmed that the amplification of the PLEKHF2 gene is related to the decreased survival of patients with prostate cancer [45]. FGFRL1 would be significantly upregulated in patients with ovarian cancer, and high FGFRL1 expression is associated with a poor prognosis. The loss of FGFRL1 function significantly affects the proliferation, apoptosis, and migration of ovarian cancer cells in vitro and tumor growth in vivo [46]. For hsa-miR-107, Wang et al. have found that Integral Membrane Protein 2C (ITM2C) would be underexpressed and miR-107-5p overexpressed in acute aortic dissection tissues [47]. Meanwhile, the overexpression of miR-107-5p would promote cell proliferation and inhibit cell apoptosis in rat aortic smooth muscle cells (RASMCs).

CCR7 is a proinflammatory cytokine and is found in human atherosclerotic plaques [48]. It has been found that the expression of CCR7 is dramatically downregulated in human carotid atherosclerotic plaques [49]. Wan et al. have assessed the protein expression of CCR7 in AAAs through immunohistochemistry, with the findings showing that CCR7 is obviously upregulated in AAA compared with the control healthy aorta [50]. Recently, the homeostatic chemokines, CCL19 and CCL21 and their receptor CCR7, have been linked to atherogenesis. The current findings suggest that CCR7 is immunohistochemically localized to macrophages and vascular SMCs via carotid atherosclerosis and highlight that although CLL19 and CCL21 are signaling through CCR7, they may have different effects on macrophages and SMC [51]. In this study, it has also been found that the hsa_circ_0081968/hsa-let-7f-5p/CCR7 regulatory axis may be an important target for the treatment of AAAs. Cai et al. have also revealed that a lncRNA/miRNA/mRNA interaction network may be involved in the pathogenesis of pulmonary arterial hypertension (PAH) and suggested hsa-let-7e-5p and CCR7 can be regarded as potential biomarkers for PAH [52].

However, there are some limitations in this study. Firstly, although it has been shown that SOD2 can cause arterial-related diseases, there is no research to prove the relationship between abdominal aortic aneurysm and SOD2. Secondly, the way of action has been predicted based on the measured RNA network, which, however, has not been confirmed (dual-luciferase reporter gene analysis, gene overexpression, or gene knockout). Although several related genes have been screened out in this study for the first time, further in vitro clinical research and in vivo experiments shall be carried out to confirm its expression and functional mechanism in terms of AAAs.

At present, there are relatively few articles on the circRNA-miRNA-mRNA network regulating AAA, and most of them are ceRNA regulatory networks based on lncRNA. Based on Gene Expression Omnibus, we detected six dysregulations of circRNA-miRNA-mRNA axes in the AAA condition compared with healthy aortic condition. The novelty of this study is that the circRNA-miRNA-mRNA network has been constructed for the first time through the GEO database. However, given that these findings are only based on bioinformatics models, it is necessary to conduct a thorough study to verify the possible effects of these 6 axes in AAAs. Further studies are needed to elucidate the biological function of these circRNAs in the pathogenesis of AAA.

Data Availability

The GEO dataset data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities (grant number DUT19RC(3)076), the National Natural Science Foundation of China (grant number 81600370), and the China Postdoctoral Science Foundation (grant number 2018M640270).

Supplementary Materials

Supplementary Figure 1: box plot of each sample in the four datasets. The blue line in the box represents the median, and circles represent outliers: (A) GSE144431; (B) GSE47472; (C) GSE7084; (D) GSE57691. Supplementary Table 1: common DEGs of different regulation methods in the three datasets. Supplementary Table 2: prediction of partial miRNAs binding to circRNAs. Supplementary Table 3: Gene Ontology term enrichment and genome pathway enrichment analysis for differentially expressed mRNAs of the upregulation group. (Supplementary Materials)

References

G. R. Upchurch Jr. and T. A. Schaub, “Abdominal aortic aneurysm,” American Family Physician, vol. 73, pp. 1198–1204, 2006.
View at: Google Scholar
S. A. Han, J. H. Joh, and H. C. Park, “Risk factors for abdominal aortic aneurysm in the Korean population,” Annals of Vascular Surgery, vol. 41, pp. 135–140, 2017.
View at: Publisher Site | Google Scholar
M. Stuntz, “Modeling the burden of abdominal aortic aneurysm in the USA in 2013,” Cardiology, vol. 135, no. 2, pp. 127–131, 2016.
View at: Publisher Site | Google Scholar
R. Umebayashi, H. A. Uchida, and J. Wada, “Abdominal aortic aneurysm in aged population,” Aging, vol. 10, no. 12, pp. 3650-3651, 2018.
View at: Publisher Site | Google Scholar
Lin Lin and C. H. Yi, “Infrarenal abdominal aortic aneurysm causing abdominal pain,” The American Journal of the Medical Sciences, vol. 360, 2020.
View at: Publisher Site | Google Scholar
M. Wassef, G. R. Upchurch Jr., H. Kuivaniemi, R. W. Thompson, and M. D. Tilson 3rd., “Challenges and opportunities in abdominal aortic aneurysm research,” Journal of Vascular Surgery, vol. 45, no. 1, pp. 192–198, 2007.
View at: Publisher Site | Google Scholar
L. B. Schwartz, Z. K. Baldwin, and M. A. Curi, “The changing face of abdominal aortic aneurysm management,” Annals of Surgery, vol. 238, Supplement 6, pp. S56–S66, 2003.
View at: Publisher Site | Google Scholar
M. Gawenda and J. Brunkwall, “Ruptured abdominal aortic Aneurysm,” Deutsches Ärzteblatt International, vol. 109, pp. 727–732, 2012.
View at: Publisher Site | Google Scholar
R. Lozano, M. Naghavi, K. Foreman et al., “Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010,” The Lancet, vol. 380, no. 9859, pp. 2095–2128, 2012.
View at: Publisher Site | Google Scholar
H. Abdulameer, H. Al Taii, S. G. Al-Kindi, and R. Milner, “Epidemiology of fatal ruptured aortic aneurysms in the United States (1999-2016),” Journal of Vascular Surgery, vol. 69, no. 2, pp. 378–384.e2, 2019, e372.
View at: Publisher Site | Google Scholar
M. Sprynger, M. Willems, H. Van Damme, B. Drieghe, J. C. Wautrecht, and M. Moonen, “Screening program of abdominal aortic aneurysm,” Angiology, vol. 70, no. 5, pp. 407–413, 2019.
View at: Publisher Site | Google Scholar
M. U. Ali, D. Fitzpatrick-Lewis, J. Miller et al., “Screening for abdominal aortic aneurysm in asymptomatic adults,” Journal of Vascular Surgery, vol. 64, no. 6, pp. 1855–1868, 2016.
View at: Publisher Site | Google Scholar
US Preventive Services Task Force, D. K. Owens, K. W. Davidson et al., “Screening for abdominal aortic Aneurysm,” JAMA, vol. 322, no. 22, pp. 2211–2218, 2019.
View at: Publisher Site | Google Scholar
F. Zhang, R. Zhang, X. Zhang et al., “Comprehensive analysis of circRNA expression pattern and circRNA-miRNA-mRNA network in the pathogenesis of atherosclerosis in rabbits,” Aging, vol. 10, no. 9, pp. 2266–2283, 2018.
View at: Publisher Site | Google Scholar
S. Huang, X. Li, H. Zheng et al., “Loss of super-enhancer-regulated circRNA Nfix induces cardiac regeneration after myocardial infarction in adult mice,” Circulation, vol. 139, no. 25, pp. 2857–2876, 2019.
View at: Publisher Site | Google Scholar
T. Barwari, A. Joshi, and M. Mayr, “MicroRNAs in cardiovascular disease,” Journal of the American College of Cardiology, vol. 68, no. 23, pp. 2577–2584, 2016.
View at: Publisher Site | Google Scholar
X. Zhai, Y. Zhang, S. Xin, P. Cao, and J. Lu, “Insights into the involvement of circular RNAs in autoimmune diseases,” Frontiers in Immunology, vol. 12, article 622316, 2021.
View at: Publisher Site | Google Scholar
Y. Shi, X. Jia, and J. Xu, “The new function of circRNA: translation,” Clinical & Translational Oncology, vol. 22, no. 12, pp. 2162–2169, 2020.
View at: Publisher Site | Google Scholar
T. B. Hansen, T. I. Jensen, B. H. Clausen et al., “Natural RNA circles function as efficient microRNA sponges,” Nature, vol. 495, no. 7441, pp. 384–388, 2013.
View at: Publisher Site | Google Scholar
Y. Han, H. Zhang, C. Bian et al., “Circular RNA expression: its potential regulation and function in abdominal aortic aneurysms,” Oxidative Medicine and Cellular Longevity, vol. 2021, Article ID 9934951, 21 pages, 2021.
View at: Publisher Site | Google Scholar
J. Yue, T. Zhu, J. Yang et al., “CircCBFB-mediated miR-28-5p facilitates abdominal aortic aneurysm via LYPD3 and GRIA4,” Life Sciences, vol. 253, article 117533, 2020.
View at: Publisher Site | Google Scholar
R. Yang, Z. Wang, G. Meng, and L. Hua, “Circular RNA CCDC66 facilitates abdominal aortic aneurysm through the overexpression of CCDC66,” Cell Biochemistry and Function, vol. 38, no. 7, pp. 830–838, 2020.
View at: Publisher Site | Google Scholar
F. Zhao, T. Chen, and N. Jiang, “CDR1as/miR-7/CKAP4 axis contributes to the pathogenesis of abdominal aortic aneurysm by regulating the proliferation and apoptosis of primary vascular smooth muscle cells,” Experimental and Therapeutic Medicine, vol. 19, pp. 3760–3766, 2020.
View at: Publisher Site | Google Scholar
L. Tian, X. Hu, Y. He, Z. Wu, D. Li, and H. Zhang, “Construction of lncRNA-miRNA-mRNA networks reveals functional lncRNAs in abdominal aortic aneurysm,” Experimental and Therapeutic Medicine, vol. 16, pp. 3978–3986, 2018.
View at: Publisher Site | Google Scholar
M. Kanehisa and S. Goto, “KEGG: Kyoto Encyclopedia of Genes and Genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000.
View at: Publisher Site | Google Scholar
M. Liu, Q. Wang, J. Shen, B. B. Yang, and X. Ding, “Circbank: a comprehensive database for circRNA with standard nomenclature,” RNA Biology, vol. 16, no. 7, pp. 899–905, 2019.
View at: Publisher Site | Google Scholar
D. Y. Li, A. Busch, H. Jin et al., “H19 induces abdominal aortic aneurysm development and progression,” Circulation, vol. 138, no. 15, pp. 1551–1568, 2018.
View at: Publisher Site | Google Scholar
S. Kumar, R. A. Boon, L. Maegdefessel, S. Dimmeler, and H. Jo, “Role of noncoding RNAs in the pathogenesis of abdominal aortic aneurysm,” Circulation Research, vol. 124, no. 4, pp. 619–630, 2019.
View at: Publisher Site | Google Scholar
L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. P. Pandolfi, “A _ceRNA_ Hypothesis: The Rosetta Stone of a Hidden RNA Language?” Cell, vol. 146, no. 3, pp. 353–358, 2011.
View at: Publisher Site | Google Scholar
Z. Wang, C. Ren, L. Yang et al., “Silencing of circular RNA_0000326 inhibits cervical cancer cell proliferation, migration and invasion by boosting microRNA-338-3p-dependent down-regulation of CDK4,” Aging, vol. 13, no. 6, pp. 9119–9134, 2021.
View at: Publisher Site | Google Scholar
C. Yu, D. Li, Q. Yan et al., “Circ_0005927 inhibits the progression of colorectal cancer by regulating miR-942-5p/BATF2 axis,” Cancer Management and Research, vol. 13, pp. 2295–2306, 2021.
View at: Publisher Site | Google Scholar
E. Strauss, J. Tomczak, R. Staniszewski, and G. Oszkinis, “Associations and interactions between variants in selenoprotein genes, selenoprotein levels and the development of abdominal aortic aneurysm, peripheral arterial disease, and heart failure,” PLoS One, vol. 13, no. 9, article e0203350, 2018.
View at: Publisher Site | Google Scholar
L. Miao and D. K. St. Clair, “Regulation of superoxide dismutase genes: implications in disease,” Free Radical Biology & Medicine, vol. 47, no. 4, pp. 344–356, 2009.
View at: Publisher Site | Google Scholar
Z. Yu, K. Morimoto, J. Yu, W. Bao, Y. Okita, and K. Okada, “Endogenous superoxide dismutase activation by oral administration of riboflavin reduces abdominal aortic aneurysm formation in rats,” Journal of Vascular Surgery, vol. 64, no. 3, pp. 737–745, 2016.
View at: Publisher Site | Google Scholar
A. Ashtekar, D. Huk, A. Magner, K. M. D. La Perle, L. Boucai, and L. S. Kirschner, “Alterations in Sod2-induced oxidative stress affect endocrine cancer progression,” The Journal of Clinical Endocrinology and Metabolism, vol. 103, no. 11, pp. 4135–4145, 2018.
View at: Publisher Site | Google Scholar
F. Ding, K. Sun, N. Sun, Q. Jiang, M. Cao, and Z. Wu, “iTRAQ-based proteomics reveals SOD2 as a potential salivary biomarker in liver cancer,” The International Journal of Biological Markers, vol. 34, pp. 221–231, 2019.
View at: Publisher Site | Google Scholar
D. R. Hodge, B. Peng, C. Pompeia et al., “Epigenetic silencing of manganese superoxide dismutase (SOD-2) in KAS 6/1 human multiple myeloma cells increases cell proliferation,” Cancer Biology & Therapy, vol. 4, no. 5, pp. 585–592, 2005.
View at: Publisher Site | Google Scholar
X. Han, X. Liu, X. Wang et al., “TNF‐α‐dependent lung inflammation upregulates superoxide dismutase-2 to promote tumor cell proliferation in lung adenocarcinoma,” Molecular Carcinogenesis, vol. 59, no. 9, pp. 1088–1099, 2020.
View at: Publisher Site | Google Scholar
Y. T. Tsai, H. Y. Yeh, C. T. Chao, and C. K. Chiang, “Superoxide dismutase 2 (SOD2) in vascular calcification: a focus on vascular smooth muscle cells, calcification pathogenesis, and therapeutic strategies,” Oxidative Medicine and Cellular Longevity, vol. 2021, Article ID 6675548, 9 pages, 2021.
View at: Publisher Site | Google Scholar
J. N. Wang, N. Shi, and S. Y. Chen, “Manganese superoxide dismutase inhibits neointima formation through attenuation of migration and proliferation of vascular smooth muscle cells,” Free Radical Biology & Medicine, vol. 52, no. 1, pp. 173–181, 2012.
View at: Publisher Site | Google Scholar
N. R. Madamanchi, S. K. Moon, Z. S. Hakim et al., “Differential activation of mitogenic signaling pathways in aortic smooth muscle cells deficient in superoxide dismutase isoforms,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 25, no. 5, pp. 950–956, 2005.
View at: Publisher Site | Google Scholar
A. Qu, C. Jiang, M. Xu et al., “PGC-1α attenuates neointimal formation via inhibition of vascular smooth muscle cell migration in the injured rat carotid artery,” American Journal of Physiology. Cell Physiology, vol. 297, no. 3, pp. C645–C653, 2009.
View at: Publisher Site | Google Scholar
J. Q. Song, L. Y. Jiang, C. P. Fu et al., “Heterozygous SOD2 deletion deteriorated chronic intermittent hypoxia-induced lung inflammation and vascular remodeling through mtROS-NLRP3 signaling pathway,” Acta Pharmacologica Sinica, vol. 41, no. 9, pp. 1197–1207, 2020.
View at: Publisher Site | Google Scholar
A. R. Asif, M. Hecker, and M. Cattaruzza, “Disinhibition of SOD-2 expression to compensate for a genetically determined NO deficit in endothelial cells--brief report,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 2009, p. 29, 2009.
View at: Publisher Site | Google Scholar
E. Shamsara and J. Shamsara, “Bioinformatics analysis of the genes involved in the extension of prostate cancer to adjacent lymph nodes by supervised and unsupervised machine learning methods: The role of SPAG1 and PLEKHF2,” Genomics, vol. 112, no. 6, pp. 3871–3882, 2020.
View at: Publisher Site | Google Scholar
H. Tai, Z. Wu, S. Sun, Z. Zhang, and C. Xu, “FGFRL1 promotes ovarian cancer progression by crosstalk with hedgehog signaling,” Journal of Immunology Research, vol. 2018, Article ID 7438608, 11 pages, 2018.
View at: Publisher Site | Google Scholar
Z. Wang, X. Zhuang, B. Chen, D. Feng, G. Li, and M. Wei, “The role of miR-107 as a potential biomarker and cellular factor for acute aortic dissection,” DNA and Cell Biology, vol. 39, no. 10, pp. 1895–1906, 2020.
View at: Publisher Site | Google Scholar
J. K. Damås, C. Smith, E. Øie et al., “Enhanced expression of the homeostatic chemokines CCL19 and CCL21 in clinical and experimental Atherosclerosis,” Arteriosclerosis, Thrombosis, and Vascular Biology, vol. 27, no. 3, pp. 614–620, 2007.
View at: Publisher Site | Google Scholar
T. Nickel, S. Pfeiler, C. Summo et al., “oxLDL downregulates the dendritic cell homing factors CCR7 and CCL21,” Mediators of Inflammation, vol. 2012, Article ID 320953, 10 pages, 2012.
View at: Publisher Site | Google Scholar
L. Wan, J. Huang, H. Ni, and G. Yu, “Screening key genes for abdominal aortic aneurysm based on gene expression omnibus dataset,” BMC Cardiovascular Disorders, vol. 18, no. 1, p. 34, 2018.
View at: Publisher Site | Google Scholar
B. Halvorsen, T. B. Dahl, L. M. Smedbakken et al., “Increased levels of CCR7 ligands in carotid atherosclerosis: different effects in macrophages and smooth muscle cells,” Cardiovascular Research, vol. 102, no. 1, pp. 148–156, 2014.
View at: Publisher Site | Google Scholar
M. Cai, X. Li, H. Dong, Y. Wang, and X. Huang, “CCR7and its related molecules may be potential biomarkers of pulmonary arterial hypertension,” FEBS Open Bio, vol. 11, no. 6, pp. 1565–1578, 2021.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Hao Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

451

Downloads

676

Citations