Abstract Abdominal aortic aneurysm (AAA) is a prevalent aortic disease that causes high mortality due to asymptomatic gradual expansion and sudden rupture. The underlying molecular mechanisms and effective pharmaceutical therapy for preventing AAA progression have not been fully identified. In this study, we identified the key modules and hub genes involved in AAA growth from the [34]GSE17901 dataset in the Gene Expression Omnibus (GEO) database through the weighted gene co-expression network analysis (WGCNA). Key genes were further selected and validated in the mouse dataset ([35]GSE12591) and human datasets ([36]GSE7084, [37]GSE47472, and [38]GSE57691). Finally, we predicted drug candidates targeting key genes using the Drug–Gene Interaction database. Overall, we identified key modules enriched in the mitotic cell cycle, GTPase activity, and several metabolic processes. Seven key genes (CCR5, ADCY5, ADCY3, ACACB, LPIN1, ACSL1, UCP3) related to AAA progression were identified. A total of 35 drugs/compounds targeting the key genes were predicted, which may have the potential to prevent AAA progression. Keywords: abdominal aortic aneurysm, weighted gene co-expression network, key module, hub gene, functional enrichment, drug–gene prediction 1. Introduction Abdominal aortic aneurysm (AAA) is a localized dilation or bulging of the abdominal aorta, commonly occurring in the infrarenal region [[39]1]. Most patients with AAA remain asymptomatic for years or even decades. It is estimated that around 200,000 AAA rupture cases are diagnosed worldwide annually, and the mortality after rupture remains around 80% [[40]2,[41]3,[42]4]. Currently, AAA requiring intervention, e.g., large aneurysms with a diameter more than 5.5 cm, aneurysms that expand rapidly in a short period, or aneurysms that compromise the perfusion to distant organs are indicated for open surgical or endovascular aortic repair. However, the outcomes from these measures are not so satisfactory [[43]5,[44]6]. For patients with small AAAs or those who are not eligible for AAA repair, close aneurysm surveillance and adjuvant therapy are recommended [[45]5]. So far, no effective pharmacological treatments have been developed to prevent AAA growth or rupture [[46]7,[47]8]. Hence, there is a need to elucidate the possible mechanisms of AAA progression and explore corresponding pharmaceutical treatments. A number of preclinical mouse AAA models have been developed to understand the pathogenesis of AAA [[48]9,[49]10]. Among these models, angiotensin II-infused ApoE^−/− mice are the commonly used [[50]11,[51]12,[52]13,[53]14,[54]15]. Although the inherent pathology of aneurysm is different between mice and humans, it shares some of the important properties of human AAA, like pronounced inflammatory responses and aortic rupture [[55]11,[56]12,[57]13,[58]14,[59]15]. Based on the findings from mouse models and human samples, AAA is currently accepted as an inflammation-driven disease, as many related processes (such as infiltration of macrophages, neutrophils, B cells and T cells, and activation of inflammatory pathways) were found both in humans and mice [[60]16,[61]17,[62]18,[63]19]. Overactivation of the inflammatory response leads to the destruction of aortic media through the release of proteolytic enzymes and the death of vascular smooth muscle cells, which further promote AAA development [[64]20]. Several studies based on the high-throughput microarray profiling further confirmed the involvement of the above biological processes in AAA, including the immune response, chronic inflammation, and reactive oxygen species [[65]21,[66]22,[67]23]. Dozens of genes related to AAA development were identified through gene expression profiles [[68]24,[69]25,[70]26]. However, these studies exclusively focused on the differentially expressed genes (DEGs) between AAA and control groups, which ignored some key genes that are highly correlated to specific sample traits of AAA. Weighted gene co-expression network analysis (WGCNA) is a bioinformatics algorithm developed by Horvath et al. [[71]27]. By constructing a scale-free weighted network, WGCNA can investigate biologically meaningful gene sets connected to sample features and explore inner module hub genes that are highly associated inside the co-expression module. WGCNA has been successfully used to identify key modules and hub genes related to cardiovascular diseases, such as atherosclerosis, heart failure, and acute myocardial infarction [[72]28,[73]29,[74]30]. So far, data collected at different time points of AAA progression have not been subjected to WGCNA analysis to identify the critical modules and hub genes. In this study, WGCNA analysis was performed using the explore dataset [75]GSE17901 in the Gene Expression Omnibus (GEO) database. Key modules of AAA development and hub genes in each module were identified. Gene functional enrichment analysis of key modules was applied to show their potential biological activities. Hub genes were screened in the STRING database and further selected in the Cytoscape software (San Diego, CA, USA). Key genes from hub genes were validated using mouse AAA model [76]GSE12591 dataset and human AAA sample [77]GSE7084, [78]GSE47472, and [79]GSE57691 datasets. Candidate drugs for AAA treatment were screened in the Drug Gene Interaction Database (DGIdb) based on the above-identified key genes. 2. Materials and Methods 2.1. Data Sources and Preprocessing The workflow of this study is shown in [80]Figure 1. Datasets related to AAA—[81]GSE17901, [82]GSE12591, [83]GSE7084, [84]GSE47472 and [85]GSE57691 ([86]Table 1) were downloaded from the GEO database (accessed on 1 April 2020 from [87]https://www.ncbi.nlm.nih.gov/geo/). In the explore dataset [88]GSE17901 [[89]26], aortic samples were taken on day 7, day 14, and day 28 from ApoE^−/− mice treated by angiotensin II or saline. The diameters of the treated aortas increased throughout the 28-day course, which we defined as the progression of AAA, so samples with AAA (n = 18) were selected for weighted gene co-expression network (WGCNA) analysis. Mouse dataset ([90]GSE12591) and human datasets ([91]GSE7084, [92]GSE47472, and [93]GSE57691) were used to validate the hub genes. The [94]GSE12591 dataset included 18 mouse aortas exposed to saline (n = 6) or angiotensin II (n = 12) infusion [[95]25]. The [96]GSE7084 included control samples (n = 10) and AAA samples from patients (n = 9) [[97]24]. The [98]GSE47472 contained AAA neck specimen (n = 14) and normal aortic tissue from organ donors (n = 8). The [99]GSE57691 included AAA samples (n = 49) and normal aortic specimens of organ donors (n = 10) [[100]31]. Each dataset was processed by background correction, including removal of batch effect using the sva R package (version 3.12) and quantile normalization with the limma R package (version 3.38.3) [[101]32] for further analysis. Figure 1. [102]Figure 1 [103]Open in a new tab Flowchart of analysis in the study. [104]GSE17901 was a mouse dataset containing AAA samples collected on day 7, day 14 and day 28, which was used for exploring the key modules and hub genes related to AAA progression. Hub genes were identified through the STRING database and Cytoscape software (San Diego, CA, USA). Key genes were further selected from the hub genes and validated in the mouse ([105]GSE12591) and human ([106]GSE7084, [107]GSE47472 and [108]GSE57691) AAA datasets. Finally, potential drugs or compounds targeting these key genes were screened in the DGIdb database. AAA: abdominal aortic aneurysm. The flowchart was created with BioRender.com (accessed on 11 April 2021). Table 1. GSE datasets included in the study. Catalog. GSE Dataset Organism Sample Number * PMID Explore dataset [109]GSE17901 Mouse AAA day7: 7, AAA day14: 5, AAA day28: 6 21712436 Validate dataset [110]GSE12591 Mouse Control: 6, AAA: 5 19580648 [111]GSE7084 Human Donor: 10, AAA: 9 17634102 [112]GSE47472 Human Donor: 8, AAA: 14 NA [113]GSE57691 Human Donor: 10, AAA: 49 NA [114]Open in a new tab *: Number of samples (control or AAA) used in this study; NA: not applicable. 2.2. Construction of WGCNA The WGCNA R package (version 1.69) was used to perform the weighted co-expression network analysis. Genes with the top 25% variance from the explore dataset [115]GSE17901 were selected for the following analysis step. Using the pick Soft Threshold function, the soft-thresholding power was determined and used to construct a scale-free network. Thereafter, gene co-expression modules were identified using the one-step network construction method and labeled with different colors. The reassign threshold was set at 0.25, and the minimum number of genes in each module was 30. 2.3. Selection of Key Modules Corresponding to Sample Traits To explore the key modules that are significantly associated with sample traits of AAA, we calculated the relevancy between module eigengene (ME), which summarizes each module’s expression profiles. The correlation results were shown using the ggcorrplot R package (version 0.1.3) [[116]33]. Furthermore, Gene Significance (GS) was quantified by the absolute value of the association between the gene expression and sample trait. In every module, measurement of module membership (MM) was defined as the correlation of the ME and gene expression profile. Modules with high significance (p-value < 0.05) and relationships (correlation >0.6 or <−0.6) were defined as key modules of AAA and used for hub gene selection. 2.4. Functional Enrichment Analysis of the Key Modules To understand the biological activities of genes in key modules, we conducted Gene Ontology (GO) function enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis with the clusterProfiler R package (version 3.10) [[117]34]. Adjusted p-value < 0.05 was considered a statistically significant difference in enrichment analysis, and the top 10 of each analysis were extracted for visualization. 2.5. Identification of Hub Genes in the Key Modules Hub genes are those that have a high degree of intramodular connectivity. In this study, hub genes were defined as the top 10% of genes from key modules with the highest connectivity. We uploaded them into the search tool for the retrieval of the interacting genes (STRING) website (accessed on 1 May 2020 from [118]www.string-db.org) for protein–protein interaction analysis, choosing the confidence >0.4 [[119]35]. Cytoscape software (San Diego, CA, USA) was used for network visualization and hub gene selection [[120]36]. The top 10 hub genes in each module were selected with the maximal clique centrality (MCC) method using cytoHubba plugin software in Cytoscape (San Diego, CA, USA) [[121]37]. 2.6. Hub Genes Validation and Key Genes Selection The validation of hub genes was performed by comparing the normalized gene expression value between control and AAA groups. The validated datasets [122]GSE12591, [123]GSE7084, [124]GSE47472, and [125]GSE57691 were downloaded from the GEO database, and data were preprocessed as mentioned before. In the [126]GSE12591 mouse dataset, the gene expression of the selected hub genes in AAA and controls were compared, and genes with p < 0.05 were confirmed as the key genes. In the [127]GSE7084, [128]GSE47472, and [129]GSE57691 human datasets, genes were extracted as described for dataset [130]GSE12591. Genes with p < 0.05 were confirmed as the key genes. Common genes in both the mouse dataset and human datasets were defined as the final key genes. 2.7. Predication of Drug–Gene Interaction The Drug–Gene Interaction Database (DGIdb) (accessed on 8 June 2020 from [131]http://www.dgidb.org/) is an online database of drug–gene interaction data aggregated from various sources, including several drug databases (DrugBank, PharmGKB, ChEMBL), clinical trial databases, and literature from PubMed [[132]38]. The selected key genes that were considered the potential pharmaceutical targets for AAA treatment were imported into DGIdb to explore existing drugs or small organic compounds. Results were displayed using the R packages ggplot2 (version 3.2.1) [[133]39] and ggalluvial (version 0.11.1) [[134]40]. 2.8. Statistical Analysis To define the statistical significance of differences between the two groups, we performed analysis using a non-parametric test or t-test based on data distribution characteristics. All analyses were conducted with R software (version 3.5.5). p-value < 0.05 was assigned statistical significance. 3. Results 3.1. Construction of Weighted Gene Co-Expression Network After cleaning the data in the explore dataset [135]GSE17901 by WGCNA, 5408 genes from 17 samples were analyzed for co-expression network construction. A scale-free network was constructed with a soft-threshold at nine, and a correlation coefficient threshold set at 0.85 ([136]Figure 2A), and 15 related co-expression modules were obtained ([137]Figure 2B). Four main clusters were observed. The turquoise module (1394 genes) was the biggest cluster, followed by the blue module (897 genes), brown module (793 genes), and yellow module (586 genes). All the ungrouped genes (199 genes) were included in the grey module. Figure 2. [138]Figure 2 [139]Open in a new tab Construction of gene co-expression network by WGCNA. (A) Determination of soft-thresholding power for scale-free network construction. Here, we set the coefficient threshold at 0.85, and the soft-threshold was 9; (B) cluster analysis of the dendrogram and identification of co-expressed modules. In this study, we got 15 related co-expression modules. 3.2. Construction of Module-Trait Relationships and Detection of Key Modules The related sample traits (time—day 7, day 14, day 28; dissection of abdominal aorta) were obtained from the sample information in the [140]GSE17901 dataset ([141]Figure S1A). The relationships between these traits and each module were defined by the correlation between ME and sample traits ([142]Figure 3, [143]Figure S1B). These results indicated that three modules (blue, green, and brown) were strongly related to the time trait, representing the progression of AAA ([144]Figure 3, [145]Figure S2A–C). Blue and green modules also significantly correlated with the dissection sample trait ([146]Figure 3, [147]Figure S2D–E). Thus, the blue (897 genes), green (436 genes), and brown (793 genes) modules were defined as the key modules that were highly correlated with AAA. Figure 3. [148]Figure 3 [149]Open in a new tab Identification of the key modules associated with AAA progression. Green, blue and brown modules were highly correlated (correlation > 0.6 or −0.6 and p-value < 0.01) to the time of sample collecting which stands for AAA progression. Besides, green and blue modules were also related to the dissection happening in the AAA sample (correlation > 0.6 or −0.6 and p-value < 0.01). AAA: abdominal aortic aneurysm. 3.3. Functional Enrichment Analysis of Genes in the Module To investigate the biological functions of key modules related to sample traits, we conducted GO and KEGG enrichment analysis for genes in every key module. The GO analysis showed that genes in the blue modules were mainly involved in the organelle fission, regulation of mitotic cell cycle, and nuclear division related to cell development or differentiation ([150]Figure 4A). The green module was involved in GTPase activity ([151]Figure 4B), and the brown module was clustered in cellular metabolic processes, especially cofactor metabolism, purine-containing compound metabolism, and purine nucleotide metabolism ([152]Figure 4C). The results of the KEGG analysis revealed that the blue module was enriched in fluid shear stress and atherosclerosis pathway, highly related to the progression of AAA ([153]Figure 5A). Genes in the green module were enriched in the regulation of lipolysis in the adipocyte pathway and the pancreatic secretion pathway ([154]Figure 5B). The brown module was enriched in the citrate cycle (TCA cycle) pathway ([155]Figure 5C). Figure 4. [156]Figure 4 [157]Open in a new tab Gene ontology enrichment analysis of key modules of AAA progression. (A) blue module; (B) green module; (C) brown module. Count—the number of genes in the given GO term. Gene ration—the percentage of total genes in the given GO term. Figure 5. [158]Figure 5 [159]Open in a new tab KEGG pathway enrichment analysis of key modules. (A) blue module; (B) green module; (C) brown module. Count—the number of genes in the given KEGG pathway. Rich factor—the ratio of the number of genes annotated in a pathway to the number of all genes annotated in this pathway. 3.4. Identification of Hub Genes in the Key Modules To explore the hub genes that regulate AAA development, we imported the top 10% of genes with the highest connectivity into the String online database for protein–protein interaction detection, and networks were formed in Cytoscape (San Diego, CA, USA) (the PPI networks were stored in the NDEx: accessed on 11 December 2020 from [160]https://bit.ly/37XZZWh; [161]https://bit.ly/3a7Q2sc; [162]https://bit.ly/38fyckz). With the cytoHubba plugin using the MCC method, the top 10 hub genes were identified in the key modules, namely, in the blue module (Ccr5, Fpr2, Ccr2, Fpr1, P2ry12, Hcar1, Ppbp, Aif1, Sirpb1b, Clec4n), green module (Gnai1, Adcy5, Adcy3, Rnase2a, Cxcl13, Clca1, Ear10, Ear1, Npr1, Ccl11), and brown module (Lpl, Dgat2, Fasn, Acacb, Lpin1, Acsl1, Mogat1, Lep, Ucp3, Pdk4) ([163]Table 2). Table 2. Top 10 ranked genes in key modules with the MCC method in cytoHubba. Catalog Key Modules Blue Green Brown Top 10 Gene Ccr5 Gnai1 Lpl Fpr2 Adcy5 Dgat2 Ccr2 Adcy3 Fasn Fpr1 Rnase2a Acacb P2ry12 Cxcl13 Lpin1 Hcar1 Clca1 Acsl1 Ppbp Ear10 Mogat1 Aif1 Ear1 Lep Sirpb1b Npr1 Ucp3 Clec4n Ccl11 Pdk4 [164]Open in a new tab 3.5. Hub Genes Validation and Key Genes Selection To further validate and evaluate the hub genes identified through the above analysis, the mouse dataset [165]GSE12591 was checked using the same mouse angiotensin II-induced AAA model as [166]GSE17901. In the blue module, Ccr5 and P2ry12 were significantly upregulated in the AAA group ([167]Figure 6A), and Hcar1 was significantly down-regulated in the AAA group ([168]Figure 6A). In the green module, Adcy5 and Adcy3 were the two significantly expressed genes ([169]Figure 6B). All significantly expressed genes (Dgat2, Fasn, Acacb, Lpin1, Acsl1, Mogat1, Ucp3, Pdk4) in the brown module were down-regulated in the AAA group ([170]Figure 6C). In the human AAA datasets [171]GSE7084, [172]GSE47472, and [173]GSE57691, all of the significantly expressed genes were identified by comparing organ donors and AAA patients ([174]Table 3). Considering the individual differences within each sample, genes expressed significantly in every human dataset were defined as human key genes. Finally, CCR5, ADCY5, ADCY3, ACACB, LPIN1, ACSL1, and UCP3 were the common genes that showed up both in the mouse AAA dataset and human AAA datasets and these were selected as the key genes in AAA progression. Figure 6. [175]Figure 6 [176]Open in a new tab Validation of gene expression from hub genes in mouse dataset [177]GSE12591. (A) Ccr5, P2ry12 and Hcar1 were differentially expressed in the blue module; (B) Adcy5 and Adcy3 were differentially expressed in the green module; (C) Dgat2, Fasn, Acacb, Lpin1, Acsl1, Mogat1, Ucp3 and Pdk4 were differentially expressed in the brown module. *: p < 0.05, **: p < 0.01, ***: p < 0.001 (Wilcoxon rank-sum test). Table 3. Significantly expressed hub genes in human AAA datasets. Datasets Key Modules Blue Green Brown [178]GSE7084 CCR5, CCR2, FPR2, FPR1, AIF1 GNAI1, RNASE2, NPR1 NA [179]GSE47472 CCR2, FPR2, PPBP GNAI1, RNASE2, CLCA1, LYVE1 LPIN1, UCP3 [180]GSE57691 CCR2, FPR2, PPBP, CLEC6A, SIRPB1 ADCY5, ADCY3, CXCL13, CLCA1, CCL11 ACACB, LPIN1, ACSL1, LEP Human CCR5, CCR2, FPR2, PPBP, AIF1, CLEC6A, SIRPB1, FPR1 GNAI1, RNASE2, NPR1, CLCA1, LYVE1, ADCY5, ADCY3, CXCL13, CCL11 ACACB, LPIN1, ACSL1, LEP, UCP3 [181]Open in a new tab 3.6. Predication of Drug-Gene Interaction The seven key genes CCR5, ADCY5, ADCY3, ACACB, LPIN1, ACSL1, and UCP3 were used as the potential druggable targets for AAA treatment. The drug–gene interaction results from the DGIdb database revealed 35 potential target drugs/compounds for AAA treatment. Of these, 23 drugs targeted CCR5, among which maraviroc had the highest score of prediction; seven drugs targeted ACACB, two drugs each targeted ACSL1 and ADCY5, and one drug targeted LPIN1 ([182]Figure 7, [183]Table S1). No potential drugs could be identified for ADCY3 and UCP3 Figure 7. [184]Figure 7 [185]Open in a new tab Drug–gene interaction prediction of key genes. Five key genes—ACACB, ACSL1, ADCY5, CCR5 and LPIN1 were targeted in the DGIdb database. A total of 35 potential target drugs/compounds were predicted from the database. AMP: Adenosine monophosphate; ATP: Adenosine triphosphate. 4. Discussion In the present study, we used WGCNA analysis to identify the key genes involved in AAA progression and the drugs that target these genes, which could be potentially effective for the repression of AAA growth. WGCNA was performed on the available mouse dataset ([186]GSE17901), where AAA samples were obtained at day 7, day 14, and day 28 from ApoE^−/− mice treated by angiotensin II or saline. We identified three modules (blue, green, and brown) as key modules that correlated closely with AAA growth. In these three modules, we further identified hub genes using Cytoscape software (San Diego, CA, USA) and validated the model in mouse and human datasets. Seven genes—CCR5, ADCY5, ADCY3, ACACB, LPIN1, ACSL1, and UCP3 were identified as the key genes in AAA progression. Finally, using the DGIdb database, we identified 35 drugs as potential candidates/compounds that could target the key genes and yield beneficial effects in treating AAA. WGCNA is a systematic biological method that describes the gene co-expression pattern between different samples. It identifies gene sets with highly coordinated variations. The candidate biomarkers or targets of the disease are based on the connectivity between gene modules and sample traits. Compared to the traditional differential gene expression analyses, which focus solely on genes characterizing the difference between groups, WGCNA groups co-expressed genes in an unbiased manner into modules that can be connected to sample traits. Among the 15 co-expression modules obtained by WGCNA, the blue, green, and brown modules were mostly related to the AAA progression. The enrichment analysis of these key modules’ biological functions and pathways revealed that genes in the blue module were mainly enriched in the cellular process, particularly the regulation of the mitotic cell cycle. This has also been reported in several studies. For instance, Butt et al. performed peripheral blood transcriptome profiling of individuals with AAA and healthy donors. They described that significantly expressed genes were enriched in this GO term [[187]41]. Another study showed that the mitotic cell cycle was also significantly associated with dilated aortic perivascular adipose tissue [[188]16]. The most enriched pathway of the blue module in KEGG was fluid shear stress and atherosclerosis. Several studies have shown the association of atherosclerosis disease with AAA [[189]1,[190]42]. Shear stress induced by abnormal blood flow was also previously reported to contribute to the growth or rupture of AAA [[191]43,[192]44]. The GO analysis of the green module showed that the biological process of GTPase activity was involved in AAA development. Dysregulation of GTPase activity would influence normal functions of endothelial cells and vascular smooth cells, including re-endothelialization, cell migration, and proliferation [[193]45,[194]46]. KEGG pathway enrichment of genes in the green module demonstrated that the regulation of lipolysis in the adipocyte pathway is also engaged in AAA growth. Adventitia of the aorta which contains the mass of adipocytes is a new direction of AAA research. One recent study revealed the key regulatory factors in perivascular adipose tissue of AAA [[195]19]. Another study further proved that the increase in AAA diameter was correlated with lipid-related processes in the adventitia [[196]18]. Results from functional enrichment analysis of the brown module indicated that some metabolic processes or pathways are also involved in AAA progression. In our study, cofactor metabolism was the most enriched process. This is in agreement with previously published studies that have shown that cofactors like cobalamin (vitamin B12) and glutathione could slow down the progression of AAA to some extent [[197]47,[198]48]. These findings confirm the involvement of the mitotic cell cycle, GTPase activity, and metabolic process in the pathogenesis of AAA. The hub genes in the present study were selected by a combined analysis of gene intramodular connectivity and protein–protein interaction in the STRING database and Cytoscape software (San Diego, CA, USA). These selected hub genes were further confirmed in mouse and human datasets with gene differential expression analysis. Seven key genes were eventually identified—CCR5 from the blue module, ADCY5 and ADCY3 from the green module, ACACB, LPIN1, ACSL1, and UCP3 from the brown module. The vital role of CCR5, C-C motif chemokine receptor 5, in HIV-1 infection has been accepted since the discovery of this receptor [[199]49]. It is expressed in many immune cells, including macrophages, T cells, and natural killer cells. CCR5 and its ligands regulate the inflammatory response by affecting the biological activities of the above-mentioned immune cells [[200]50]. The results from [201]GSE12591 identifying Ccr5 as a differential gene upregulated in the mouse aortas with aneurysms [[202]25]. CCR5 signaling in the macrophage pathway was enriched by functional analyses of differential genes in [203]GSE7084 [[204]24]. Furthermore, patients with AAA frequently have CCR5 Delta 32 deletion mutations and are vulnerable to rupture of aneurysms [[205]51]. Thus, CCR5 may be a potential biomarker for AAA progression and an indication of rupture. The hub gene ADCY5 (mouse—Adcy5) in the green module was related to mouse AAA progression and dissection. This was consistent with the findings by Phillips et al. which showed Adcy5 was one of the differentially expressed genes in the murine dissecting AAA [[206]52]. ADCY3 is an enzyme that regulates the cyclic adenosine monophosphate (cAMP). Besides its role in AAA progression, loss of ADCY3 increases the risk of obesity and type 2 diabetes [[207]53], and the single nucleotide polymorphisms of this gene are related to hypertension [[208]54], which are the risk factors leading to the initiation of AAA [[209]55]. LPIN1, ACSL1, and UCP3 were related to adipocyte differentiation and muscle growth [[210]56,[211]57,[212]58,[213]59,[214]60], so dysregulation of these three genes may lead to AAA initiation, growth or rupture, as adipocytes residing in the perivascular tissue, and vascular smooth muscle cells play an important role in the development of AAA [[215]19,[216]60]. According to the reviewed literature, the remaining key gene ACACB had no apparent connection with AAA. This, however, requires further investigation to clarify its function in AAA progression. So far, there is no effective drug therapy for the prevention of AAA progression or rupture. In this study, seven key genes were identified and used for predicting drug-gene interactions. A total of potential 35 drugs or compounds were presented in the DGIdb database. Most of these targeted the CCR5 gene. We checked these 35 candidates from the literature and ClinicalTrials.gov (accessed on 18 July 2020 from [217]https://clinicaltrials.gov/), the largest clinical trials database containing over 329,000 trials worldwide. Five targetable drugs (PF-05175157, firsocostat, and metformin targeting ACACB; maraviroc targeting CCR5; rosiglitazone targeting LPIN1) were found to be used for AAA treatment. PF-05175157 and firsocostat are two novel acetyl-CoA carboxylase (ACC) inhibitors for lipid disorders [[218]61,[219]62], which could potentially rebalance dysregulated lipid metabolism in AAA to limit the development of the disease. Metformin is the first-line oral antidiabetic drug [[220]63]. It also has proven effects on cardiovascular diseases through the reduction of inflammation and oxidative stress [[221]64,[222]65,[223]66]. Several epidemiological studies have indicated that the use of metformin use could decrease yearly AAA growth [[224]67,[225]68]. Though maraviroc is a CCR5 antagonist prescribed for HIV-1 treatment, it could also be applied for AAA treatment since it was reported that maraviroc could reduce cardiovascular risk by modulation of atherosclerotic progression in vivo and in vitro [[226]69,[227]70]. Rosiglitazone (RGZ) is a potent peroxisome proliferator-activated receptor-γ (PPAR-γ) agonist that can protect against ischemia/reperfusion injury due to its anti-inflammatory effects [[228]71]. It has been reported that RGZ reduces stent-induced neointimal formation by decreasing the inflammatory responses and vascular smooth muscle hyperplasia [[229]72]. Through the same anti-inflammatory effect, RGZ could also inhibit the growth and rupture of mouse aortic aneurysms induced by angiotensin II and high cholesterol [[230]73]. No drugs could be predicted for the ADCY3 and UCP3 genes. These two gene candidates will have to be evaluated as potential targets in AAA treatment in further studies. Though our study is the first that performed WGCNA analysis with samples collected at different points of time in AAA growth, this study still has some limitations. Firstly, upon screening of the public database mouse dataset, [231]GSE17901 was the only dataset available that allowed us to follow gene function over time and was used as an exploration dataset for WGCNA analysis. As a result, the sample size used for WGCNA analysis (n = 17) just passed the minimum official criteria (n = 15), therefore there may be noise for the biological network construction. The angiotensin II-induced AAA in mice may share similar features with human AAA, but the inherent pathology is different and thus, our results should be interpreted with caution. This study has indeed predicted interesting key genes involved in the progression of AAA and potentially useful drugs, however these findings should be validated further with in vitro and in vivo models of AAA. In summary, this study identified key co-expression modules, key genes, and several critical biological processes related to AAA progression. With drug–gene interaction prediction, target drugs or compounds may provide the possibility of developing a medical treatment for AAA. 5. Conclusions Our study using WGCNA analyses revealed seven key genes (CCR5, ADCY5, ADCY3, ACACB, LPIN1, ACSL1, UCP3) in three modules correlated to AAA progression. Mitotic cell cycle, GTPase activity, and metabolic process were involved in the pathogenesis of AAA. The therapeutic potential of several predicted drugs for the treatment of AAA could be further explored. Acknowledgments