Next Article in Journal
R-α-Lipoic Acid and 4-Phenylbutyric Acid Have Distinct Hypolipidemic Mechanisms in Hepatic Cells
Next Article in Special Issue
Weighted Gene Co-Expression Network Analysis Reveals Key Genes and Potential Drugs in Abdominal Aortic Aneurysm
Previous Article in Journal
Relationship between Diet, Microbiota, and Healthy Aging
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomedicines
Volume 8
Issue 8
10.3390/biomedicines8080288
biomedicines-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Regulatory Network Analysis of Perivascular Adipose Tissue of Abdominal Aortic Aneurysm Identifies Master Regulators of Key Pathogenetic Pathways

by
Luca Piacentini
*,
Mattia Chiesa
and
Gualtiero Ivanoe Colombo
Immunology and Functional Genomics Unit, Centro Cardiologico Monzino, IRCCS, 20138 Milan, Italy
*
Author to whom correspondence should be addressed.
Biomedicines 2020, 8(8), 288; https://0-doi-org.brum.beds.ac.uk/10.3390/biomedicines8080288
Submission received: 1 July 2020 / Revised: 30 July 2020 / Accepted: 12 August 2020 / Published: 14 August 2020
(This article belongs to the Special Issue Aortic Aneurysms: Vascular Remodeling and Repair)
Browse Figures
Review Reports Versions Notes

Abstract

:
The lack of medical therapy to treat abdominal aortic aneurysm (AAA) stems from our inadequate understanding of the mechanisms underlying AAA pathogenesis. To date, the only available treatment option relies on surgical intervention, which aims to prevent AAA rupture. Identifying specific regulators of pivotal pathogenetic mechanisms would allow the development of novel treatments. With this work, we sought to identify regulatory factors associated with co-expressed genes characterizing the diseased perivascular adipose tissue (PVAT) of AAA patients, which is crucially involved in AAA pathogenesis. We applied a reverse engineering approach to identify cis-regulatory elements of diseased PVAT genes, the associated transcription factors, and upstream regulators. Finally, by analyzing the topological properties of the reconstructed regulatory disease network, we prioritized putative targets for AAA interference treatment options. Overall, we identified NFKB1, SPIB, and TBP as the most relevant transcription factors, as well as MAPK1 and GSKB3 protein kinases and RXRA nuclear receptor as key upstream regulators. We showed that these factors could regulate different co-expressed gene subsets in AAA PVAT, specifically associated with both innate and antigen-driven immune response pathways. Inhibition of these factors may represent a novel option for the development of efficient immunomodulatory strategies to treat AAA.

1. Introduction

Despite intensive efforts over the last decades, abdominal aortic aneurysm (AAA) remains an elusive disease for which no effective treatment aiming to hinder or reduce its growth is yet available [1,2]. This reflects our incomplete understanding of the etiology and pathogenetic mechanisms leading to the development and evolution of AAA. Mechanistic studies mainly rely on animal models, which despite the large number of different methods available to induce an “acute” form of AAA, have not been individually able to thoroughly elucidate the pathogenesis of the disease, which instead presents a natural history characterized by a multifactorial, slow, and chronic process [3]. Thus, integration of in vivo models with knowledge from preclinical research on human patients is needed for a greater understanding of the processes underlying AAA [2].
A successful approach to dissecting complex phenotypes is molecular profiling, which allows a large-scale exploration of pathological processes in diseased tissues without the need for an a priori selection of the factors to be tested; therefore, it is potentially capable of uncovering new and unrecognized causes of disease onset or progression [4,5]. We have recently used this approach to explore the transcriptome of the perivascular adipose tissue (PVAT) in patients with abdominal aortic diseases [6,7]. PVAT has received growing interest in the study of large artery diseases, either because it has a fundamental role in the regulation of vascular physiology [8,9,10] or because its dysfunction is recognized to affect the development of both dilated and atherosclerotic aortic diseases [11]. Indeed, by comparing the adipose layer of a dilated abdominal aorta with that of a non-dilated aortic neck in each patient, we revealed locally restricted gene expression patterns characterizing the dilated PVAT of AAA. Overall, these genes were functionally associated with inflammatory and innate or adaptive immune responses, which along with other relevant pathways, including cell-death and extracellular matrix degradation, led us to hypothesize that AAA is an immunological disease with a possible underlying autoimmune component [6]. However, the regulatory elements underlying the pivotal pathogenetic processes associated with PVAT in AAA patients remained to be defined.
Reverse engineering from RNA expression data is a valuable and grounded approach that allows reconstruction of gene regulatory networks by identifying cis-regulatory elements, which are the targets of sequence-specific transcription factors (TFs) [12,13]. TFs are often defined as the master regulators of cellular processes because they can control the simultaneous expression of many genes at once [14]. TFs regulate the transcription of their target genes by interacting with other TFs and co-factors to constitute transcriptional complexes, but can also be directly targeted by other regulators acting upstream of the TFs, e.g., protein kinases. Identifying cis-regulatory elements and their respective TFs from co-expressed genes can, thus, provide important insights into the regulatory mechanism affecting a specific biological process.
In the present work, we aimed to identify the regulators that may control the most prominent pathogenetic processes we found associated with PVAT in AAA patients. For this purpose, we first applied a reverse engineering approach on co-expressed genes, distinguishing the dilated PVAT (diseased) from the non-dilated (healthy) aortic neck to find cis-regulatory elements and associated TFs. Then, we sought upstream regulators that could directly affect (e.g., activate) the TFs identified above or help form transcriptional complexes. Finally, exploring the topological properties of the resulting regulatory network, we outlined the most likely putative targets for AAA interference treatment options.

2. Materials and Methods

2.1. Study Population and Gene Expression Data

This study relies on the cohort of 30 AAA patients characterized in the previous study by Piacentini et al. [6]. Patient features, including demographic data, risk factors, medications, and exclusion criteria, have been extensively described in the aforementioned work. AAA patients enrolled in the study underwent elective surgery at the Centro Cardiologico Monzino IRCCS, Milan, Italy, between 2010 and 2014. Elective repair of AAA was done in compliance with the international and national guidelines for the care and treatment of AAA [15,16]. From each patient, adipose tissue specimens were collected at the time of surgery, which included the periaortic adipose tissue obtained from the aortic neck proximal to the aneurismal sac (non-dilated PVAT) and periaortic adipose tissue surrounding the aneurysmal sac (dilated PVAT).
The expression data produced in the study have been made publicly available at the NCBI’s GEO repository and can be accessed at https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE119717.

2.2. Selection of Differentially Expressed (DE) Genes in Dilated PVAT of AAA

Based on our previous results, we selected 172 unique genes positively associated with dilated AAA, corresponding to the most significant set of DE probes (transcripts), which distinguished dilated (diseased) from non-dilated (healthy) PVAT. Supplementary Data 1a shows the annotation of these 172 DE genes. We used these DE genes as the seed gene list for the identification of cis-regulatory elements and candidate TFs.

2.3. Identification of Cis-Regulatory Elements and TFs

To identify regulatory elements, we relied on a reverse engineering approach to infer transcriptional regulatory network underlying the 172 DE genes by cis-regulatory sequence analysis. This analysis was performed in the Cytoscape environment v3.7.1 [17] through the iRegulon software v1.3 [18]. Briefly, iRegulon performs a “rank-and-recovery” procedure. The ranking step allows genes (RefSeq annotation) to be ranked for a library of positional weight matrices (PWMs), which represent matrices of regulatory motifs. Then, for each gene a regulatory search space around the transcription start site (TSS) is scanned for homotypic cis-regulatory modules (CRM) using a hidden Markov model across multiple vertebrate species. A ranked list of genes is generated, with the most likely genomic target of a specific motif at the top of the ranking. In the recovery step, the enrichment of genes from the input gene list (i.e., DE genes) is tested against the gene rankings by calculating the area under the cumulative recovery curve (AUC) in the top of the ranking (3%), which is then normalized into a normalized enrichment score (NES). The method exploits a wide collection of 18 databases of 9713 non-redundant TF motifs and 3 databases of 1120 chromatin immunoprecipitation sequencing (ChIP-seq) signals along the genome (i.e., tracks).
Parameters used for the ranking step were: motive collection: “10K (9713 PWMs)” (i.e., the most extensive motif collections); track collection: “1120 ChIP-seq tracks (ENCODE raw signals)” (i.e., full collection of ChIP-seq data against TFs); putative regulatory region: “20 kb centered around TSS” (i.e., it may return promoter-based and/or distal regulators); motif rankings database: “20 kb centered around TSS (7 species)” (including the mammalian species Bos taurus, Canis familiaris, Mus musculus, Monodelphis domestica, Pan troglodytes, Macaca mulatta, and Rattus norvegicus, as well as considering the conservation among them).
Parameters for the recovery step were: enrichment score threshold = “3.0” (raising the post-hoc threshold to 3.5); receiver operating characteristic (ROC) threshold for AUC calculation = “0.03”; rank threshold = “5000”; minimum identity between orthologous genes = “0”; maximum false discovery rate (FDR) on motif similarity = “0.001”.
Resulting motifs and tracks were ranked according to NES and labeled with the motif or track identifier (ID) of their original database. Motifs that shared a higher level of similarity were grouped into clusters using the default method (see [19]). The motif-to-TF association procedure for the algorithm, based on motif similarity and orthology, returns the TFs that more likely bind to the enriched motif. The TFs with the highest levels of confidence (i.e., TFs recorded as having a direct annotation, meaning that the PWM was determined for a certain TF in that species) were selected as the most reliable candidate TFs for each cluster of motifs or tracks.

2.4. Inferring the TFs’ Upstream Regulatory Factors

To infer regulatory modules upstream of the “hub” candidate TFs, additional proteins that directly connect to TFs were firstly identified (i) by using experimentally reported protein–protein interactions (PPI) or protein complexes; and then (ii) by identifying protein kinases regulating the above extended transcriptional complexes. The Expression2kinase (X2K) software v1.6.1207 [20] was used with default parameters for both protein network expansion and kinase retrieval, except for the “allow a maximum of 10,000 node links from a node”, “allow a maximum of 100,000 interactions from an article”, and “allow a minimum of 1 article reporting a specific interaction” options, which were enabled to ensure a higher quality of interactions in outputted additional proteins. The method for drawing PPI exploits experimentally validated mammalian interactions from 18 databases containing more than 24,000 proteins and almost 390,000 interactions. The kinase–substrate interactions are instead from other sources for a consolidated dataset of 14,374 interactions from more than 3400 publications on 436 kinases (see [20] for details).

2.5. Topological Analysis

Topological analysis was performed on two different networks through Network Analyzer v3.3.2 [21] and CentiScape 2.2 software [22] into the Cytoscape environment v3.7.1. The first analyzed network was the whole regulatory network (cf. Figure 2) drawn to connect candidate TFs (regulators; source nodes) with their target DE genes (regulated; target nodes). This facilitated visual inspection of clusters of shared and unique target regulated genes and their relationships. To disentangle the complexity of the network and to find the most relevant “hub” TFs, three centrality measures, i.e., connectivity degree (aka degree), betweenness, and radiality, were evaluated. The higher the values of these topological indexes among nodes, the higher the relevance of a specific node in the network in relation to the others [23]. TFs ranked in the upper tertile for the three centrality measures were selected as the “hub” master regulators.
The second topological analysis was performed on the regulatory transcriptional complex network made of “hub” TFs, additional or intermediate proteins, and upstream kinases (cf. Figure 3). The network was treated as a directed network because the regulator–regulated target relationships were highlighted by directed edges. Measures of in-going and out-going direct connections (edges) of a node were defined as the in-degree and out-degree, respectively. Additionally, radiality, stress, betweenness, bridging, centroid, closeness, eccentricity, eigenvector centralities were also assessed. The relevance of the topological index for each node was weighted through the z-score calculation, which assesses how a specific observation (node) moves away from the mean of the total observations, measured in terms of standard deviations from the mean. A node with a z-score ≥2 was deemed as very relevant for that topological index.
For an easier interpretation of topological indexes [22,23], it should be considered that nodes with a large degree of connectivity are defined as “hubs”, since they are likely crucial factors that might play key (e.g., causative) roles in the biological context (e.g., disease) of interest [24,25]. Nodes with high betweenness (and with similar but not equal stress, bridging, closeness, and eigenvector values) are defined as “bottlenecks”, which are likely important factors that can hold together communicating proteins or genes and which might be relevant as organizing (or very central) regulatory molecules [26]. Nodes with high eccentricity show how easily a protein or gene can be functionally influenced by all other proteins or genes in the network. High centroid values instead show those nodes that are functionally able to organize subclusters of proteins or genes, thus possibly coordinating the activity of nodes with high connectivity. A node with high radiality, high eccentricity, and high closeness provides a consistent indication that it plays a central position in the network.

2.6. Linking Transcriptional Complex Clusters with AAA Pathogenetic Biological Functions

Functional association of transcriptional complex clusters (each including a subset of the seed list of DE genes and additional or intermediate proteins and kinases) with AAA pathogenetic biological processes or pathways was performed by measuring the overlap of the transcriptional complex cluster gene sets with the results of the original gene set enrichment analysis (GSEA; see Supplementary Materials of [6]). The collection of transcriptional complex cluster gene sets was generated in the Gene Matrix Transposed file format (*.gmt) and imported into the Enrichment Map software v3.2.1 [27] to visualize and measure the overlap of these gene sets with the enrichment network drawn from GSEA, which displayed all of the significant Gene Ontology–biological process (GO-BP) orpathway gene sets that were suggested to be pathogenically associated with dilated PVAT in AAA. The consistency of the overlap was tested through a hypergeometric test. Associations were deemed to be significant for adjusted p-values < 0.01 (Benjamini–Hochberg method for multiple testing correction) and with a gene set overlap corresponding to ≥5% of the genes for each specific transcriptional complex cluster. To facilitate visualization of the most significant relationships of each transcriptional complex cluster with their associated GO-BP or pathways, an enrichment subnetwork was subsequently generated.

3. Results

For this work, we selected the 172 unique genes we found overexpressed in dilated PVAT compared with non-dilated PVAT in AAA patients (see Supplementary Materials of [6]). Gene annotation is reported in Data Supplement 1a. These genes were positively associated with the layer of adipose tissue surrounding the dilated (diseased) abdominal aorta and included several factors of the inflammatory and immune responses, which were claimed as being the most relevant for the pathogenesis of AAA [28,29]. We, thus, sought the regulatory sequence motifs (i.e., target gene nucleotide sequences used to control its expression) and tracks (i.e., ChIP-seq signals along the genome) for these 172 genes and reconstructed a gene regulatory network.

3.1. Identification of Motifs and Tracks of Genes Over-Expressed in Dilated PVAT of AAA Patients

By cis-regulatory sequence analysis, which leverages the use of combined multiple collections of motifs and tracks, we identified 30 motifs and 1 track with NES values >3.5 and AUC values >0.06. The selected NES and AUC thresholds were chosen to reduce the probability of recovering false-positive associations, which would, therefore, ensure more robust and accurate results [18]. Exploring the motif enrichment results, we observed that subgroups of the 30 motifs showed a high sequence overlap, which were, thus, grouped into 8 different clusters based on similarity (M1 to 8; Data Supplement 1b). Notably, each cluster included one TF with a direct annotation, which represents the highest level of confidence for a motif-to-TF association. Thus, we ranked the TF with the direct annotation as the most reliable candidate TF for each related motif or track cluster (Table 1, Figure 1, and Data Supplement 1b).

3.1.1. Regulatory Network of TFs and Over-Expressed Genes in Dilated PVAT

To visually explore the relationships between candidate TFs and their relative target genes, we drew a gene regulatory network consisting of TFs as source nodes, presenting direct evidence for significant motifs with the highest NES values, as well as the DE genes as target nodes (Figure 2). The network shows the complex relationships between TFs and their targets, with TFs displaying both shared and unique target genes.

3.1.2. Selection of “hub” TFs through Topological Analysis

Although all of the identified TFs presented a significant association with DE genes, we tried to figure out the regulatory elements that more likely associate with the most relevant pathogenetic processes in dilated PVAT of AAA. For this purpose, we applied a network topology analysis to identify TFs representing the most relevant “hub” in the whole regulatory network. With this approach, we found that three TFs, i.e., TBP, NFKB1, and SPIB, displayed the highest values for all the assessed topological centrality measures (Table 2). The advantage of this approach is that we could extract key TFs by taking into account both the statistics of motif and track identification and TF “centrality” ranking in a complex gene regulatory network. In particular, the degree index indicates the relevance that SPIB, NFKB1, and TBP have within the network—they displayed the highest connectivity by interacting with large numbers of the 172 DE genes (104, 91, and 74, respectively), suggesting a central regulatory role.

3.2. Identification of Transcriptional Complexes and Upstream Regulators

We refined our search for regulatory molecules by including those that physically connect to and act on the identified “hub” TFs. To this end, we first extended the regulatory network by retrieving intermediate direct interactors of SPIB, NFKB1, and TBP, and by drawing putative transcriptional complexes. Then, we inferred upstream protein kinases that could regulate those transcriptional complexes.

3.2.1. Connecting Additional Proteins to TFs through Protein–Protein Interactions (PPI)

By leveraging experimentally validated PPI, we identified 28 intermediate proteins (annotated in Data Supplement 1d) directly linked to SPIB, NFKB1, and TBP. Then, we drew a directional network to show the relationships (i.e., from regulator-to-regulated target) among the “hub” TFs and all the retrieved intermediate proteins (Figure 3). NKFB1 and TBP displayed the highest direct connections, with 8 out-degree (outcoming connections) and 19 in-degree (incoming connections) and with 12 out-degree and 15 in-degree connections, respectively. SPIB presented 4 out- and 3 in-degree connections but displayed a direct edge towards TBP, suggesting that their possible interaction constitutes a transcriptional complex together with other proteins. Notably, we found 18 other TFs among the additional proteins (i.e., AR, BCL3, CEBPB, CREBBP, E2F1, ELF3, ESR1, FOS, JUN, KLF5, NCOA1, NCOA6, NR3C1, REL, RELA, RXRA, SP1, and SPI1), which may participate in forming such transcriptional complexes (Figure 3; green ovals). Consistently, some of them were also predicted to associate with motif clusters M2 (i.e., BCL3, E2F1, REL, RELA), M3 (i.e., SPI1), and M8 (i.e., CEBPB; cf. Data Supplement 1b), strengthening the idea that these TFs may cooperate with “hub” TFs to target common sequence motifs and coordinate transcription of DE genes. Furthermore, relying on the Transcription co-Factor DataBase (TcoF-DB) [30], we could annotate 7 proteins that act as transcriptional co-factors (i.e., that may regulate transcription by interacting with TFs), but in contrast to TFs these do not bind directly to regulatory DNA regions (Figure 3; orange ovals). Six of them (i.e., CTNNB1, HDAC1, HMGB1, RUVBL2, SIN3A, and TRIP4) have experimental evidence for both involvement in transcriptional regulation and for presence in the cell nucleus (classified as class high-confidence); RUVBL1 has evidence: inferred from electronic annotation” (classified as class 2) for involvement in transcriptional regulation.

3.2.2. Identification of Protein Kinases Upstream Transcriptional Complexes.

The expansion of the network to protein–TF interactions allowed an increase of the number of possible regulatory molecules that can influence specific functions in the dilated PVAT of AAA. The search for further kinase–substrate interactions can add an upstream level of control of gene expression and helps to identify key target regulators. Indeed, by kinase enrichment analysis, we found 42 associations with protein kinases (nominal p-values < 0.01), of which 28 stood correction for multiple testing (adjusted p-values < 0.01; Data Supplement 1c). To focus on the most significant kinases, we selected those that were ranked in the upper tertile of the distribution according to the “combined score”, i.e., CHUK, CSNK2A1, CSNK2A2, GSK3B, MAPK1, MAPK14, MAPK3, MAPK8, PRKDC, and TAF1 (Figure 3, light blue ovals; annotation in Data Supplement 1d).
Network analysis showed that CHUK, CSNK2A1, GSK3B, and PRKDC were directly associated with NFKB1; TAF1 targeted TBP; and CSNK2A1, CSNK2A2, GSK3B, MAPK3, MAPK8, and MAPK14 interacted with SPIB. Specifically, GSK3B, MAPK1, and CSNK2A1 displayed the highest out-degree values (17, 15, and 13, respectively), which suggests that they can be functionally relevant for most of the protein–TF interactions in the regulatory network.

3.2.3. Topological Analysis of Regulatory Transcriptional-Complex

Topological analysis of the regulatory transcriptional complex network (cf. Figure 3), inferred the relative importance that specific proteins and genes may have as regulators in the network (Data Supplement 1e). First, we observed that centrality measures with similar meanings were highly correlated, indicating their consistency (Figure 4). Indeed, out-degree, radiality, and closeness, which indicate the possibility that a node is functionally relevant for several others, showed high positive correlations, with Pearson’s coefficient r values ranging from 0.89 to 0.99 (p-values < 0.001). Similarly, betweenness, bridging, and stress indexes, whose values if elevated suggest that a gene or protein likely connects (i.e., holds together) pivotal regulatory molecules, were positively correlated, with r values ranging from 0.65 to 0.97 (p-values < 0.001).
Specifically, we showed that NFKB1 showed the highest degree (total connectivity), in-degree, and eccentricity, suggesting that it is likely affected by the activity of many other proteins. Interestingly, REL, which is part of the well-known NK-kB complex together with NFKB1, has high eccentricity as well, although with lower global connectivity. NFKB1 also presents high betweenness and stress, suggesting the capability to hold together communicating proteins and to play a role as an organizing regulatory molecule. Similarly, RELA, which is also a subunit of the NF-kB complex, has the highest betweenness, stress, and bridging centrality values, leading to the hypothesis that it can serve in an organizing module along with NFKB1, “bridging” another possible set of regulatory proteins. For example, RELA targets TBP, which also displays very high betweenness and stress indexes and may serve as a central regulatory protein. Another remarkable finding concerns the role of MAPK1. Indeed, we can postulate that MAPK1 could be functionally essential for several other proteins, as indicated by its eigenvector index, out-degree (n = 15), and high closeness and radiality values compared to the mean of the network, which overall suggests a prominent regulatory role for this kinase in the network. Similarly, GSK3B showed a high out-degree (n = 17), radiality, and closeness over the network mean. Together these results provide a consistent suggestion that both MAPK1 and GSK3B kinases might have central regulatory roles in the network. CREBBP, as with NFKB1, RELA, and TBP, showed a high-stress index, indicating that it may be functionally relevant in connecting other regulatory molecules. CREBBP presents a high total degree of connectivity (15 in-degree and 13 out-degree connections), with a direct connection towards NFKB1 and incoming connections from RELA and SPIB. TBP displayed high stress and betweenness, which underlie its important roles as a “hub” TF and an organizing regulatory protein. Finally, RXRA showed the most consistent bridging index together with RELA, indicating its possible role as a connector of other central regulatory proteins. Indeed, RXRA is directly linked to very central proteins, such as the TFs NFKB1, RELA, and TBP, and the kinases MAPK1 and GSK3B. In line with this, we observed that the connection between RELA and RXRA showed the highest “edge-betweenness” value (33; with network mean ± SD = 8 ± 5, over 311 total edges), which indicates that this connection separates highly interconnected subgraph clusters.

3.3. Association of Regulatory Subnetworks with AAA Pathogenetic Pathways

We eventually built the entire regulatory network obtained by reverse engineering, including DE genes, candidate TFs, intermediate protein interactors, and upstream kinase regulators. By selecting the first neighbors of each “hub” TF throughout the network, we generated three “transcriptional clusters” to infer specific functional association with AAA-related pathways. This approach allowed matching each transcriptional cluster with those GO-BP or pathways associated with dysfunctional PVAT that could have a pathogenetic role in the development or progression of AAA (see Supplementary Materials in [6]). The significance of the overlapping genes, which measure the consistency of the association between each transcriptional cluster and specific GO-BP pathways, was estimated through a hypergeometric test. We found 69, 54, and 32 GO-BP pathways with a significant number of genes overlapping those of the transcriptional clusters of NFKB1, SPIB, and TBP, respectively (adjusted I-values < 0.01 and overlap size threshold ≥ 5% for the genes of the transcriptional cluster; Data Supplement 2). Interestingly, the 3 transcriptional clusters presented both shared unique associations with pathogenetic GO-BP pathways. To summarize results and reduce redundancy of the overlapping GO-BP pathways, we drew an enrichment network (Figure 5). We found that the NFKB1 transcriptional cluster was uniquely associated with “regulation of lymphocytes proliferation” (specifically T-cells), “regulation of protein secretion”, and “vasculature development”. The SPIB transcriptional cluster was univocally associated with “regulation of phagocytosis”, “granulocyte or neutrophil chemotaxis”, “humoral immune response”, and “Fc receptor-mediated signaling”. TBP displayed a distinctive association with “regulation of cytokine’s biosynthetic process”. NFKB1 and TBP transcriptional clusters shared GO-BP pathways involved in TLR signaling (in particular, TRAF6-mediated induction of NFkB and MAP kinases upon endosomal TLR activation) and in the regulation of the JAK/STAT cascade. NFKB1 and SPIB transcriptional clusters shared GO-BP pathways involved in the regulation of lymphocyte activation and differentiation and inflammatory response, as well as in leukocyte cell–cell adhesion. Finally, we found common GO-BP pathways between the three transcriptional clusters, e.g., “response to molecule of bacterial origin”, “leukocyte chemotaxis”, “regulation of ERK1 and ERK2 cascade”, and “IL-4 and IL-13 signaling”.

4. Discussion

The pathogenic mechanisms responsible for AAA formation and expansion are largely unknown. PVAT is a well-known vessel homeostasis regulator that can have a major role in vascular disease pathogenesis [11,31]. Mounting evidence suggests that dysfunctional PVAT may play a role in vascular diseases, including AAA development [32]. To date, large-scale gene expression profiling of AAA has been carried out only on full-thickness aortic walls, excluding the adipose layer of the vessel [33]. Through genome-wide expression studies, we have recently gained insights into the important role that PVAT can play in both atherosclerotic and non-atherosclerotic abdominal aortic diseases in humans [6,7]. In AAA patients, our data suggested that an altered immune response in PVAT, alongside other concurring mechanisms, is a crucial pathogenetic element that leads to the progression of AAA, and ultimately to its rupture, by amplifying inflammation and degenerative mechanisms (e.g., loss of aorta structural integrity) [6].
With this work, we inferred the regulatory molecules that govern the most prominent pathogenetic processes in PVAT of AAA patients. Our study consistently shows that subsets of co-expressed genes characterizing the diseased PVAT in AAA present common cis-regulatory elements, which are the targets of specific TFs.
The reconstruction of regulatory networks proves to be particularly enlightening, as the disease phenotype is only occasionally the result of a single gene or protein effector, but is more commonly the combination of multiple pathobiological effectors acting together in complex relationships [24]. Analyzing molecular networks through their topological properties can help to identify both crucial disease genes (hubs) and potential drug targets [24,34]. However, it is worth noting that pharmacological targets are not necessarily the network “hub nodes”, because these are often essential factors with a strong influence on the phenotype of the cell or organism. Hub genes or proteins should be targeted only in particular subsets of cells with a strong and direct pathogenetic role. Otherwise, interference with the activity of essential genes may influence physiological processes that are crucial to the normal functioning of the organism. On the contrary, proteins that act as connectors (or presenting lower degrees of connectivity) have a lower impact on the overall structure of the network, thus resulting in likely safer candidate disease targets [35,36].
We identified SPIB, NFKB1, and TBP (i.e., “hub” TFs) as the master regulators of the resulting gene regulatory network, as they have the largest connectivity with the co-expressed genes associated with diseased PVAT. Therefore, we hypothesized that these TFs could be responsible for regulating the most impactful pathogenetic mechanisms of dysfunctional PVAT in AAA patients. Additionally, we identified additional proteins that can directly interact with the “hub” TFs to form transcriptional complexes, as well as protein kinases that can regulate the state of activation of “hub” TFs and transcriptional complexes.
We showed that the transcriptional clusters of SPIB, NFKB1, and TBP have strong functional associations with mechanisms of the innate and adaptive immune response, which we assumed to be the main drivers of the (auto)immune mechanisms in AAA PVAT. Specifically, we associated the NFKB1 transcriptional cluster with the positive regulation of lymphocyte proliferation, and together with TBP with the expression of genes involved in the innate immune response, i.e., the toll-like receptor (TLR) signaling. This is consistent with the well-known role of NF-kB complex activation in guiding both lymphocyte function, including cell proliferation and differentiation [37], and TLR signaling [38], which ultimately ends in NF-kB-triggered downstream transcription of genes that allow for an effective response to primary stimuli, such as antigens (or self-antigens) and danger- or pathogen-associated molecular patterns [39]. Additionally, NFKB1, REL, and RELA may differently combine to form heterodimers and carry out both unique and overlapping roles in T-cell proliferation during different stages of the cell cycle [40]. The relevance of the NF-kB complex in AAA pathogenesis has already been shown in in vivo AAA-induced models. Activation of the NF-kB complex has been proposed as a key factor inducing macrophage infiltration and osteoclastogenic differentiation [41,42] and affecting the inflammatory response in other cell types (e.g., vascular smooth muscle cells and mesenchymal cells) [43,44], ultimately leading to aortic inflammation and vessel wall degeneration. Nonetheless, our data suggest that lymphocyte activation and proliferation mechanisms could be crucial for amplifying the local antigen-driven immune response in the PVAT of AAA patients. Interfering with NK-kB signaling [45] in immune cells may, thus, have a significant impact on the evolution of AAA. Although immune-suppressive treatment should be carefully evaluated [46], selectively inhibiting NF-kB signaling in activated and proliferating lymphocytes in dysfunctional AAA PVAT may be a targeted therapeutic intervention without compromising the healthy effect mediated by NF-kB in other immune cells.
Mechanisms of the innate and adaptive immune response can also be mediated by SPIB and its cognate SPI1, which we found to be directly connected to SPIB in the regulatory network. Both SPIB and SPI1 are required for all of the signaling pathways by which B-lymphocytes sense and respond to local environmental stimuli, including antigens and molecules acting through TLRs [47]. Accordingly, we found that the SPIB transcriptional cluster was associated with response to TLR engagement, humoral immunity, and lymphocyte activation in PVAT of AAA. Interestingly, we also found a significant association with Fcγ receptor signaling, which underlies a localized antigen-specific humoral immune response. This signaling mechanism can be activated by B cells or neutrophils in the presence of immunoglobulin–antigen aggregates [48,49], and has been involved in the pathogenesis of AAA [50]. The downregulation of Fc receptor actions has been proposed as a therapeutic approach for other inflammatory and autoimmune diseases [51]. Targeting Fc receptor signaling regulators by manipulating the SPIB/SPI1 complex could be tested to block or minimize the local self-sustaining chronic inflammation observed in AAA PVAT.
Evaluating topological indexes of the extended regulatory network, we found that the protein kinases MAPK1 and GSK3B and the nuclear receptor RXRA (a type of retinoid X receptor) appear to play key roles in AAA pathogenetic mechanisms, and thus may be intended for interference therapy, as they are “non-hub” connecting proteins that could regulate signaling through the “hub” TFs. Both MAPKs and GSK3s have been recognized as having considerable roles in the immune response. MAPKs are strongly involved in innate immunity, including signaling related to TLRs [52]. Additionally, the MAPK/ERK pathway is an important modulator of matrix metalloproteinases during AAA formation in in vivo models, and its inhibition could represent a possible therapeutic approach to prevent AAA formation [53]. Consistently, we found that NFKB1, SPIB, and TBP transcriptional clusters were associated with the regulation of ERK1 and ERK2 cascade in the PVAT, reinforcing the concept that specific protein kinases play central roles in the regulation of inflammation and immune response in AAA. A direct role of GSK3s in AAA pathogenesis has not been established, although GSK3s were shown to regulate the production of pro- and anti-inflammatory cytokines and to influence the proliferation, differentiation, and survival of T-cells. GSK3s mainly work through the regulation of critical transcription factors, including NF-kB, the inhibition of which has been proposed for the treatment of several pathological conditions with underlying altered immune responses in animal models [54]. Evidence concerning the roles of retinoid X receptors in AAA is also limited [55]. It has been proposed that RXR activation may have a beneficial role in AAA by inhibiting the angiotensin type 1 receptor in vascular smooth muscle cells, a factor which is known to affect AAA development in angiotensin II–induced AAA models [56]. However, RXRA is also known to be involved in the regulation of innate immunity [57], and thus its inhibition could suppress immune-mediated pathogenetic mechanisms in AAA PVAT.
Since SPIB, TBP, MAPK1, GSK3B, and RXRA have not yet been associated with AAA in humans, they may be considered as newly identified targets for disease treatment.
Other proteins or kinases with less connectivity in the regulatory network may be potential therapeutic targets for AAA [35]. An example is histone deacetylase 1 (HDCA1), which despite its low centrality in our regulatory network, directly interacts with both NFKB1 and TBP, and thus might be considered as a possible target for interference treatments. Consistently, the use of HDACs inhibitors has been demonstrated to be effective in in vivo AAA models and has been proposed for treatment in humans [58]. Notably, HDACs have many TFs as natural substrates, including NF-kB, and are well-known regulators of the T-cell immune response. HDAC inhibition enhances T-regulatory cell survival and immune-modulating functions, and affects the development of T-cell response, including T-cell proliferation in response to antigen stimulation [59].
In summary, innate and adaptive immune responses are characterized by complex crosstalking mechanisms [60] and are both involved in AAA pathogenesis [28]. Additionally, communication between perivascular adipocytes and immune cells may also play a role in these complex relationships. For example, adipocytes secrete soluble factors (e.g., adipokines, chemokine, or pro-inflammatory interleukins), which may trigger signaling pathways in target cells involving the activation of the NF-kB complex, among others. Immune cells, in turn, produce inflammatory molecules and recognize antigens presented by the class I or II human leukocyte antigens (HLAs) expressed on adipocytes [8,61,62].
Innovative strategies for treating AAA may include both immune modulation to stimulate local anti-inflammatory mechanisms and targeting of specific pathogenetic lymphocytes and their interactions with other cells, such as antigen-presenting cells [63,64,65]. In our previous research, while showing that the activation of several pathways can contribute to the evolution of the disease, we assumed that immune response processes were essential to sustain a chronic inflammatory cycle that characterizes AAA, since they were involved in both the early and later stages of AAA [6]. It is tempting to speculate that interventions on these specific pathogenetic pathways could be therapeutically beneficial, because they could limit the main harmful mechanism that appears to be necessary for all phases of AAA. Notably, an immune phenotyping analysis in human AAA samples recently found that T-lymphocytes are the primary cell leukocyte population in AAA, with the largest concentration in PVAT, and that these PVAT-associated T-lymphocytes correlated with the severity of the disease [66]. Assessing the presence of a locally restricted, antigen-driven clonal expansion in PVAT of AAA could, thus, be the goal of future studies, which will allow for the precise detection of pathogenic lymphocytes for direct interference treatment [45], possibly through the above-suggested disease gene targets.
To our knowledge, this is the first study that describes the regulatory elements that may contribute to controlling the major pathogenetic processes observed in PVAT of AAA patients. Our study takes advantage of a reverse engineering approach that emerged as a valuable tool for elucidating cellular function and dysregulation in pathological contexts, and which promises to increase our ability to identify potential therapeutic targets and disease biomarkers as well [13].
Our work also has limitations. It relies on inferential analysis and suggests putative sensitive AAA target genes, for which formal evidence should still be provided showing that they can be used for effective interference treatments. Furthermore, our work has focused on the potential regulators of the most significant overexpressed genes restricted to the diseased PVAT of our AAA patient cohort, without ruling out the possibility that other regulatory mechanisms may still be important in the development and progression of the disease. Finally, we did not evaluate any possible interactions with genetic variants known to be associated with AAA [67]. However, given their putative relevance in aneurysm diseases, it may be important in future research to explore the existence of specific polymorphisms that could affect gene expression and offer insight into the molecular pathogenesis of AAA.

5. Conclusions

After a long-lasting period of research, our understanding of AAA is still in its “adolescence”. Despite a clear knowledge of the pathological hallmarks, the enigmatic “trait” characterizing AAA is still unknown, and is a relentless obstacle to its overall comprehension and the potential for effective care.
The translational gap between preclinical disease models and successful clinical studies prompts the research on AAA to explore and find novel candidates for treatments. Due to its critical role in AAA pathogenesis, PVAT promises to be a reliable target for testing innovative treatment options.
With this work, we identified master regulators of prominent pathogenetic processes associated with PVAT of AAA patients, i.e., altered immune response, including antigen-specific lymphocytes activation or proliferation and TLR signaling. Through the reconstruction of a gene regulatory network and the associated upstream regulators, we also suggested novel possible targets that may be considered for locally restricted interference treatments of AAA.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/2227-9059/8/8/288/s1.

Author Contributions

Conceptualization, L.P.; methodology, L.P., and M.C.; formal analysis, L.P.; writing—original draft preparation, L.P. and G.I.C.; visualization, L.P. and M.C.; supervision, G.I.C.; project administration, L.P. and G.I.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Italian Ministry of Health (Research Projects RC Nos. 2600658, 2627621, and 2631196).

Acknowledgments

We thank Rita Spirito, Claudio Saccu, and José Pablo Werba for providing original specimens.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AAAAbdominal aortic aneurysm
AUCArea under the cumulative recovery curve
ChIPChromatin immunoprecipitation
DEDifferentially expressed
GO-BPGene Ontology–biological process
GSEAGene set enrichment analysis
NESNormalized enrichment score
PWMPositional weight matrix
PPIProtein–protein interaction
PVATPerivascular adipose tissue
TcoFTranscription co-factor
TFTranscription factor
TLRToll-like receptors
TSSTranscription start site

References

  1. Golledge, J. Abdominal aortic aneurysm: Update on pathogenesis and medical treatments. Nat. Rev. Cardiol. 2019, 16, 225–242. [Google Scholar] [CrossRef] [PubMed]
  2. Lindeman, J.H.; Matsumura, J.S. Pharmacologic Management of Aneurysms. Circ. Res. 2019, 124, 631–646. [Google Scholar] [CrossRef] [PubMed]
  3. Poulsen, J.L.; Stubbe, J.; Lindholt, J.S. Animal Models Used to Explore Abdominal Aortic Aneurysms: A Systematic Review. Eur. J. Vasc. Endovasc. Surg. 2016, 52, 487–499. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tromp, G.; Kuivaniemi, H. Developments in genomics to improve understanding, diagnosis and management of aneurysms and peripheral artery disease. Eur. J. Vasc. Endovasc. Surg. 2009, 38, 676–682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Casamassimi, A.; Federico, A.; Rienzo, M.; Esposito, S.; Ciccodicola, A. Transcriptome Profiling in Human Diseases: New Advances and Perspectives. Int. J. Mol. Sci. 2017, 18, 1652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Piacentini, L.; Werba, J.P.; Bono, E.; Saccu, C.; Tremoli, E.; Spirito, R.; Colombo, G.I. Genome-Wide Expression Profiling Unveils Autoimmune Response Signatures in the Perivascular Adipose Tissue of Abdominal Aortic Aneurysm. Arter. Thromb. Vasc. Biol. 2019, 39, 237–249. [Google Scholar] [CrossRef]
  7. Piacentini, L.; Saccu, C.; Bono, E.; Tremoli, E.; Spirito, R.; Colombo, G.I.; Werba, J.P. Gene-expression profiles of abdominal perivascular adipose tissue distinguish aortic occlusive from stenotic atherosclerotic lesions and denote different pathogenetic pathways. Sci. Rep. 2020, 10, 6245. [Google Scholar] [CrossRef]
  8. Rajsheker, S.; Manka, D.; Blomkalns, A.L.; Chatterjee, T.K.; Stoll, L.L.; Weintraub, N.L. Crosstalk between perivascular adipose tissue and blood vessels. Curr. Opin. Pharmacol. 2010, 10, 191–196. [Google Scholar] [CrossRef] [Green Version]
  9. Brown, N.K.; Zhou, Z.; Zhang, J.; Zeng, R.; Wu, J.; Eitzman, D.T.; Chen, Y.E.; Chang, L. Perivascular adipose tissue in vascular function and disease: A review of current research and animal models. Arter. Thromb. Vasc. Biol. 2014, 34, 1621–1630. [Google Scholar] [CrossRef] [Green Version]
  10. Gil-Ortega, M.; Somoza, B.; Huang, Y.; Gollasch, M.; Fernández-Alfonso, M.S. Regional differences in perivascular adipose tissue impacting vascular homeostasis. Trends Endocrinol. Metab. 2015, 26, 367–375. [Google Scholar] [CrossRef]
  11. Queiroz, M.; Sena, C.M. Perivascular adipose tissue in age-related vascular disease. Ageing Res. Rev. 2020, 59, 101040. [Google Scholar] [CrossRef] [PubMed]
  12. Emmert-Streib, F.; Dehmer, M.; Haibe-Kains, B. Gene regulatory networks and their applications: Understanding biological and medical problems in terms of networks. Front. Cell Dev. Biol. 2014, 2, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Lefebvre, C.; Rieckhof, G.; Califano, A. Reverse-engineering human regulatory networks. Wiley Interdiscip. Rev. Syst. Biol. Med. 2012, 4, 311–325. [Google Scholar] [CrossRef] [PubMed]
  14. Lambert, S.A.; Jolma, A.; Campitelli, L.F.; Das, P.K.; Yin, Y.; Albu, M.; Chen, X.; Taipale, J.; Hughes, T.R.; Weirauch, M.T. The Human Transcription Factors. Cell 2018, 172, 650–665. [Google Scholar] [CrossRef] [PubMed]
  15. Moll, F.L.; Powell, J.T.; Fraedrich, G.; Verzini, F.; Haulon, S.; Waltham, M.; van Herwaarden, J.A.; Holt, P.J.E.; van Keulen, J.W.; Rantner, B.; et al. Management of abdominal aortic aneurysms clinical practice guidelines of the European society for vascular surgery. Eur. J. Vasc. Endovasc. Surg. 2011, 41, S1–S58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Italian Journal of Vascular and Endovascular Surgery 23 March 2016 (1 Suppl 1)—Minerva Medica Journals. Available online: https://www.minervamedica.it/en/journals/vascular-endovascular-surgery/issue.php?cod=R46Y2016S01 (accessed on 30 July 2020).
  17. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  18. Janky, R.; Verfaillie, A.; Imrichová, H.; Van de Sande, B.V.; Standaert, L.; Christiaens, V.; Hulselmans, G.; Herten, K.; Sanchez, M.N.; Potier, D.; et al. iRegulon: From a Gene List to a Gene Regulatory Network Using Large Motif and Track Collections. PLoS Comput. Biol. 2014, 10, e1003731. [Google Scholar] [CrossRef] [Green Version]
  19. Verfaillie, A.; Imrichova, H.; Janky, R.; Aerts, S. iRegulon and i-cisTarget: Reconstructing Regulatory Networks Using Motif and Track Enrichment. Curr. Protoc. Bioinform. 2015, 52, 2.16.1–2.16.39. [Google Scholar] [CrossRef]
  20. Chen, E.Y.; Xu, H.; Gordonov, S.; Lim, M.P.; Perkins, M.H.; Ma’ayan, A. Expression2Kinases: mRNA profiling linked to multiple upstream regulatory layers. Bioinformatics 2012, 28, 105–111. [Google Scholar] [CrossRef] [Green Version]
  21. Assenov, Y.; Ramírez, F.; Schelhorn, S.-E.; Lengauer, T.; Albrecht, M. Computing topological parameters of biological networks. Bioinformatics 2008, 24, 282–284. [Google Scholar] [CrossRef] [Green Version]
  22. Scardoni, G.; Petterlini, M.; Laudanna, C. Analyzing biological network parameters with CentiScaPe. Bioinformatics 2009, 25, 2857–2859. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Doncheva, N.T.; Assenov, Y.; Domingues, F.S.; Albrecht, M. Topological analysis and interactive visualization of biological networks and protein structures. Nat. Protoc. 2012, 7, 670–685. [Google Scholar] [CrossRef] [PubMed]
  24. Barabási, A.-L.; Gulbahce, N.; Loscalzo, J. Network medicine: A network-based approach to human disease. Nat. Rev. Genet. 2011, 12, 56–68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Stelzl, U.; Worm, U.; Lalowski, M.; Haenig, C.; Brembeck, F.H.; Goehler, H.; Stroedicke, M.; Zenkner, M.; Schoenherr, A.; Koeppen, S.; et al. A human protein-protein interaction network: A resource for annotating the proteome. Cell 2005, 122, 957–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yu, H.; Kim, P.M.; Sprecher, E.; Trifonov, V.; Gerstein, M. The Importance of Bottlenecks in Protein Networks: Correlation with Gene Essentiality and Expression Dynamics. PLoS Comput. Biol. 2007, 3, e59. [Google Scholar] [CrossRef] [PubMed]
  27. Merico, D.; Isserlin, R.; Stueker, O.; Emili, A.; Bader, G.D. Enrichment map: A network-based method for gene-set enrichment visualization and interpretation. PLoS ONE 2010, 5, e13984. [Google Scholar] [CrossRef]
  28. Dale, M.A.; Ruhlman, M.K.; Baxter, B.T. Inflammatory cell phenotypes in AAAs: Their role and potential as targets for therapy. Arter. Thromb. Vasc. Biol. 2015, 35, 1746–1755. [Google Scholar] [CrossRef] [Green Version]
  29. Folkesson, M.; Vorkapic, E.; Gulbins, E.; Japtok, L.; Kleuser, B.; Welander, M.; Länne, T.; Wågsäter, D. Inflammatory cells, ceramides, and expression of proteases in perivascular adipose tissue adjacent to human abdominal aortic aneurysms. J. Vasc. Surg. 2017, 65, 1171–1179.e1. [Google Scholar] [CrossRef] [Green Version]
  30. Schmeier, S.; Alam, T.; Essack, M.; Bajic, V.B. TcoF-DB v2: Update of the database of human and mouse transcription co-factors and transcription factor interactions. Nucleic Acids Res. 2017, 45, D145–D150. [Google Scholar] [CrossRef] [Green Version]
  31. Cao, Z.F.H.; Stoffel, E.; Cohen, P. Role of Perivascular Adipose Tissue in Vascular Physiology and Pathology. Hypertension 2017, 69, 770–777. [Google Scholar] [CrossRef]
  32. Horimatsu, T.; Kim, H.W.; Weintraub, N.L. The Role of Perivascular Adipose Tissue in Non-atherosclerotic Vascular Disease. Front. Physiol. 2017, 8, 969. [Google Scholar] [CrossRef] [PubMed]
  33. Estrelinha, M.; Hinterseher, I.; Kuivaniemi, H. Gene expression studies in human abdominal aortic aneurysm. Rev. Vasc. Med. 2014, 2, 77–82. [Google Scholar] [CrossRef]
  34. Hao, T.; Wang, Q.; Zhao, L.; Wu, D.; Wang, E.; Sun, J. Analyzing of Molecular Networks for Human Diseases and Drug Discovery. Curr. Top. Med. Chem. 2018, 18, 1007–1014. [Google Scholar] [CrossRef] [PubMed]
  35. Safari-Alighiarloo, N.; Taghizadeh, M.; Rezaei-Tavirani, M.; Goliaei, B.; Peyvandi, A.A. Protein-protein interaction networks (PPI) and complex diseases. Gastroenterol. Hepatol. Bed Bench 2014, 7, 17–31. [Google Scholar]
  36. Rivas, J.D.L.; Prieto, C. Protein Interactions: Mapping Interactome Networks to Support Drug Target Discovery and Selection. In Bioinformatics and Drug Discovery; Larson, R.S., Ed.; Methods in Molecular Biology; Humana Press: Totowa, NJ, USA, 2012; pp. 279–296. ISBN 978-1-61779-965-5. [Google Scholar]
  37. Gerondakis, S.; Siebenlist, U. Roles of the NF-κB Pathway in Lymphocyte Development and Function. Cold Spring Harb. Perspect. Biol. 2010, 2, a000182. [Google Scholar] [CrossRef] [Green Version]
  38. Kawai, T.; Akira, S. Signaling to NF-kappaB by Toll-like receptors. Trends Mol. Med. 2007, 13, 460–469. [Google Scholar] [CrossRef]
  39. Hayden, M.S.; West, A.P.; Ghosh, S. NF-κB and the immune response. Oncogene 2006, 25, 6758–6780. [Google Scholar] [CrossRef] [Green Version]
  40. Gerondakis, S.; Grumont, R.; Gugasyan, R.; Wong, L.; Isomura, I.; Wong, H.; Banerjee, A. Unravelling the complexities of the NF-kappaB signalling pathway using mouse knockout and transgenic models. Oncogene 2006, 25, 6781–6799. [Google Scholar] [CrossRef] [Green Version]
  41. Shiraya, S.; Miwa, K.; Aoki, M.; Miyake, T.; Oishi, M.; Kataoka, K.; Ohgi, S.; Ogihara, T.; Kaneda, Y.; Morishita, R. Hypertension Accelerated Experimental Abdominal Aortic Aneurysm Through Upregulation of Nuclear Factor κB and Ets. Hypertension 2006, 48, 628–636. [Google Scholar] [CrossRef] [Green Version]
  42. Takei, Y.; Tanaka, T.; Kent, K.C.; Yamanouchi, D. Osteoclastogenic Differentiation of Macrophages in the Development of Abdominal Aortic Aneurysms. Arterioscler. Thromb. Vasc. Biol. 2016, 36, 1962–1971. [Google Scholar] [CrossRef] [Green Version]
  43. Cui, R.; Tieu, B.; Recinos, A.; Tilton, R.G.; Brasier, A.R. RhoA mediates angiotensin II-induced phospho-Ser536 nuclear factor kappaB/RelA subunit exchange on the interleukin-6 promoter in VSMCs. Circ. Res. 2006, 99, 723–730. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Ijaz, T.; Sun, H.; Pinchuk, I.V.; Milewicz, D.M.; Tilton, R.G.; Brasier, A.R. Deletion of NF-κB/RelA in Angiotensin II-Sensitive Mesenchymal Cells Blocks Aortic Vascular Inflammation and Abdominal Aortic Aneurysm Formation. Arter. Thromb. Vasc. Biol. 2017, 37, 1881–1890. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Miyake, T.; Miyake, T.; Kurashiki, T.; Morishita, R. Molecular Pharmacological Approaches for Treating Abdominal Aortic Aneurysm. Ann. Vasc. Dis. 2019, 12, 137–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Lindeman, J.H.N.; Rabelink, T.J.; van Bockel, J.H. Immunosuppression and the abdominal aortic aneurysm: Doctor Jekyll or Mister Hyde? Circulation 2011, 124, e463–e465. [Google Scholar] [CrossRef] [Green Version]
  47. Willis, S.N.; Tellier, J.; Liao, Y.; Trezise, S.; Light, A.; O’Donnell, K.; Garrett-Sinha, L.A.; Shi, W.; Tarlinton, D.M.; Nutt, S.L. Environmental sensing by mature B cells is controlled by the transcription factors PU.1 and SpiB. Nat. Commun. 2017, 8, 1426. [Google Scholar] [CrossRef]
  48. Rosales, C. Fcγ Receptor Heterogeneity in Leukocyte Functional Responses. Front. Immunol. 2017, 8, 280. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, Y.; Jönsson, F. Expression, Role, and Regulation of Neutrophil Fcγ Receptors. Front. Immunol. 2019, 10, 1958. [Google Scholar] [CrossRef]
  50. Shi, Y.; Yang, C.-Q.; Wang, S.-W.; Li, W.; Li, J.; Wang, S.-M. Characterization of Fc gamma receptor IIb expression within abdominal aortic aneurysm. Biochem. Biophys. Res. Commun. 2017, 485, 295–300. [Google Scholar] [CrossRef]
  51. Ben Mkaddem, S.; Benhamou, M.; Monteiro, R.C. Understanding Fc Receptor Involvement in Inflammatory Diseases: From Mechanisms to New Therapeutic Tools. Front. Immunol. 2019, 10, 811. [Google Scholar] [CrossRef] [Green Version]
  52. Arthur, J.S.C.; Ley, S.C. Mitogen-activated protein kinases in innate immunity. Nat. Rev. Immunol. 2013, 13, 679–692. [Google Scholar] [CrossRef]
  53. Ghosh, A.; DiMusto, P.D.; Ehrlichman, L.K.; Sadiq, O.; McEvoy, B.; Futchko, J.S.; Henke, P.K.; Eliason, J.L.; Upchurch, G.R. The role of extracellular signal-related kinase during abdominal aortic aneurysm formation. J. Am. Coll. Surg. 2012, 215, 668–680.e1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Beurel, E.; Michalek, S.M.; Jope, R.S. Innate and adaptive immune responses regulated by glycogen synthase kinase-3 (GSK3). Trends Immunol. 2010, 31, 24–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Neels, J.G.; Hassen-Khodja, R.; Chinetti, G. Nuclear receptors in abdominal aortic aneurysms. Atherosclerosis 2020, 297, 87–95. [Google Scholar] [CrossRef] [PubMed]
  56. Takeda, K.; Ichiki, T.; Funakoshi, Y.; Ito, K.; Takeshita, A. Downregulation of angiotensin II type 1 receptor by all-trans retinoic acid in vascular smooth muscle cells. Hypertension 2000, 35, 297–302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Núñez, V.; Alameda, D.; Rico, D.; Mota, R.; Gonzalo, P.; Cedenilla, M.; Fischer, T.; Boscá, L.; Glass, C.K.; Arroyo, A.G.; et al. Retinoid X receptor α controls innate inflammatory responses through the up-regulation of chemokine expression. Proc. Natl. Acad. Sci. USA 2010, 107, 10626–10631. [Google Scholar] [CrossRef] [Green Version]
  58. Galán, M.; Varona, S.; Orriols, M.; Rodríguez, J.A.; Aguiló, S.; Dilmé, J.; Camacho, M.; Martínez-González, J.; Rodriguez, C. Induction of histone deacetylases (HDACs) in human abdominal aortic aneurysm: Therapeutic potential of HDAC inhibitors. Dis. Model. Mech. 2016, 9, 541–552. [Google Scholar] [CrossRef] [Green Version]
  59. Akimova, T.; Beier, U.H.; Liu, Y.; Wang, L.; Hancock, W.W. Histone/protein deacetylases and T-cell immune responses. Blood 2012, 119, 2443–2451. [Google Scholar] [CrossRef] [Green Version]
  60. Iwasaki, A.; Medzhitov, R. Control of adaptive immunity by the innate immune system. Nat. Immunol. 2015, 16, 343–353. [Google Scholar] [CrossRef]
  61. Kim, H.W.; de Chantemèle, E.J.B.; Weintraub, N.L. Perivascular Adipocytes in Vascular Disease. Arterioscler. Thromb. Vasc. Biol. 2019, 39, 2220–2227. [Google Scholar] [CrossRef]
  62. Huh, J.Y.; Park, Y.J.; Ham, M.; Kim, J.B. Crosstalk between Adipocytes and Immune Cells in Adipose Tissue Inflammation and Metabolic Dysregulation in Obesity. Mol. Cells 2014, 37, 365–371. [Google Scholar] [CrossRef]
  63. Wang, S.K.; Murphy, M.P. Immune Modulation as a Treatment for Abdominal Aortic Aneurysms. Circ. Res. 2018, 122, 925–927. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Chang, T.W.; Gracon, A.S.A.; Murphy, M.P.; Wilkes, D.S. Exploring autoimmunity in the pathogenesis of abdominal aortic aneurysms. Am. J. Physiol. Heart Circ. Physiol. 2015, 309, H719–H727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bluestone, J.A.; Bour-Jordan, H. Current and future immunomodulation strategies to restore tolerance in autoimmune diseases. Cold Spring Harb. Perspect. Biol. 2012, 4, a007542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Sagan, A.; Mikolajczyk, T.P.; Mrowiecki, W.; MacRitchie, N.; Daly, K.; Meldrum, A.; Migliarino, S.; Delles, C.; Urbanski, K.; Filip, G.; et al. T Cells Are Dominant Population in Human Abdominal Aortic Aneurysms and Their Infiltration in the Perivascular Tissue Correlates With Disease Severity. Front. Immunol. 2019, 10, 1979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Pinard, A.; Jones, G.T.; Milewicz, D.M. Genetics of Thoracic and Abdominal Aortic Diseases. Circ. Res. 2019, 124, 588–606. [Google Scholar] [CrossRef]
Biomedicines 08 00288 g001 550
Figure 1. Target nucleotide sequences. Logos of cis-regulatory elements corresponding to the highest enriched motif of each cluster are shown. Motif IDs are reported on top of each sequence. The logo for tracks, i.e., CHD1, is not available.
Figure 1. Target nucleotide sequences. Logos of cis-regulatory elements corresponding to the highest enriched motif of each cluster are shown. Motif IDs are reported on top of each sequence. The logo for tracks, i.e., CHD1, is not available.
Biomedicines 08 00288 g001
Biomedicines 08 00288 g002 550
Figure 2. Gene regulatory network. The network shows the connectivity between candidate transcription factors (source nodes; green hexagons) and differentially expressed (DE) genes (target nodes; pink ovals). Edge color refers to clusters of DE genes based on the connections to a specific TF.
Figure 2. Gene regulatory network. The network shows the connectivity between candidate transcription factors (source nodes; green hexagons) and differentially expressed (DE) genes (target nodes; pink ovals). Edge color refers to clusters of DE genes based on the connections to a specific TF.
Biomedicines 08 00288 g002
Biomedicines 08 00288 g003 550
Figure 3. Extended regulatory network. The network shows the relationships (directed dashed edges) among the hub transcriptional factors, the additional proteins inferred by protein–protein interactions, and the most significant kinases identified by kinase enrichment analysis. Node shapes distinguish hub TFs (hexagon) from all the other inferred proteins (ovals). Node colors reflect node types or functions: hub TFs = red; inferred TFs = green; kinases = light blue; TcoFs = orange; inferred molecules with other functions = light orange. Also see Data Supplement 1d for detailed annotation.
Figure 3. Extended regulatory network. The network shows the relationships (directed dashed edges) among the hub transcriptional factors, the additional proteins inferred by protein–protein interactions, and the most significant kinases identified by kinase enrichment analysis. Node shapes distinguish hub TFs (hexagon) from all the other inferred proteins (ovals). Node colors reflect node types or functions: hub TFs = red; inferred TFs = green; kinases = light blue; TcoFs = orange; inferred molecules with other functions = light orange. Also see Data Supplement 1d for detailed annotation.
Biomedicines 08 00288 g003
Biomedicines 08 00288 g004a 550Biomedicines 08 00288 g004b 550
Figure 4. (a) Topological analysis of the extended regulatory network. Each panel represents the scatterplot of the index-over-index topological measures calculated for the extended regulatory network. To draw the representative scatterplots, only one index was chosen among those that were highly correlated and with very similar meanings. The presence of a non-uniform distribution of nodes, with most nodes having a low index and a few having a high index, identifies those nodes that clearly differ from the average measure of the index throughout the network (red and orange dots or labels). Red, orange, and blue colors, respectively, refer to nodes with z-scores > 2.0, between 1.5 and 2.0, and <1.5. Detailed results are reported in Data Supplement 1e. (b) Correlation plot of the topological indexes. Indexes with similar meanings show a highly significant Pearson’s correlation (p-values < 0.001). Numbers and circle sizes (from smaller to bigger) refer to correlation coefficient r values. Positive and negative correlations are displayed by a gradient color from white (low) to blue (high) and from white (low) to red (high), respectively. To simplify visualization, crosses mark non-significant correlations.
Figure 4. (a) Topological analysis of the extended regulatory network. Each panel represents the scatterplot of the index-over-index topological measures calculated for the extended regulatory network. To draw the representative scatterplots, only one index was chosen among those that were highly correlated and with very similar meanings. The presence of a non-uniform distribution of nodes, with most nodes having a low index and a few having a high index, identifies those nodes that clearly differ from the average measure of the index throughout the network (red and orange dots or labels). Red, orange, and blue colors, respectively, refer to nodes with z-scores > 2.0, between 1.5 and 2.0, and <1.5. Detailed results are reported in Data Supplement 1e. (b) Correlation plot of the topological indexes. Indexes with similar meanings show a highly significant Pearson’s correlation (p-values < 0.001). Numbers and circle sizes (from smaller to bigger) refer to correlation coefficient r values. Positive and negative correlations are displayed by a gradient color from white (low) to blue (high) and from white (low) to red (high), respectively. To simplify visualization, crosses mark non-significant correlations.
Biomedicines 08 00288 g004aBiomedicines 08 00288 g004b
Biomedicines 08 00288 g005 550
Figure 5. Enrichment network of transcriptional clusters. NFKB1, SPIB, and TBP transcriptional clusters (squared green nodes) and pathogenetic Gene Ontology–biological process (GO-BP) pathways associated with abdominal aortic aneurysm perivascular adipose tissue (red nodes) are connected by dashed orange edges, which indicate significant overlaps between two gene sets. Edge thickness (from thinner to thicker) is proportional to the number of overlapping genes. GO-BP pathway node colors (from lighter to darker) and sizes (from smallest to bigger) are proportional to the original normalized enrichment score calculated by gene set enrichment analysis. Groups of redundant gene sets were manually circled and labeled by relevant overview GO-BP pathways terms.
Figure 5. Enrichment network of transcriptional clusters. NFKB1, SPIB, and TBP transcriptional clusters (squared green nodes) and pathogenetic Gene Ontology–biological process (GO-BP) pathways associated with abdominal aortic aneurysm perivascular adipose tissue (red nodes) are connected by dashed orange edges, which indicate significant overlaps between two gene sets. Edge thickness (from thinner to thicker) is proportional to the number of overlapping genes. GO-BP pathway node colors (from lighter to darker) and sizes (from smallest to bigger) are proportional to the original normalized enrichment score calculated by gene set enrichment analysis. Groups of redundant gene sets were manually circled and labeled by relevant overview GO-BP pathways terms.
Biomedicines 08 00288 g005
Table
Table 1. Summary of enriched motifs and tracks aggregated in clusters.
Table 1. Summary of enriched motifs and tracks aggregated in clusters.
ClusterTFNESAUC# Targets# Motifs/Tracks
M1TBP6.200.098745
M2NFKB15.910.095917
T1CHD15.230.130291
M3SPIB5.120.0851044
M4SRF4.660.080536
M5BCL64.140.074553
M6PAX34.070.073482
M7ATF23.580.068461
M8CEBPE3.530.067502
Transcription factors (TFs) with direct annotation, normalized enrichment score (NES), and area under the cumulative recovery curve (AUC) values refer to the highest enriched motif or track of each cluster. Motif (M) and track (T) clusters are ordered by NES. #, number of.
Table
Table 2. Summary of topological indexes used to evaluate the centrality of each candidate TF in the reconstructed gene regulatory network.
Table 2. Summary of topological indexes used to evaluate the centrality of each candidate TF in the reconstructed gene regulatory network.
TF (Cluster)DegreeBetweenness CentralityRadiality
SPIB (M3)1040.363.41
NFKB1 (M2)910.223.24
TBP (M1)740.143.02
BCL6 (M5)550.082.77
SRF (M4)530.082.75
CEBPE (M8)500.082.71
PAX3 (M6)480.052.68
ATF2 (M7)460.062.65
CHD1 (T1)290.042.43
Candidate transcription factors (TFs) of each cluster are ranked by degree. The selected “hub” TFs are highlighted with bold characters.

Share and Cite

MDPI and ACS Style

Piacentini, L.; Chiesa, M.; Colombo, G.I. Gene Regulatory Network Analysis of Perivascular Adipose Tissue of Abdominal Aortic Aneurysm Identifies Master Regulators of Key Pathogenetic Pathways. Biomedicines 2020, 8, 288. https://0-doi-org.brum.beds.ac.uk/10.3390/biomedicines8080288

AMA Style

Piacentini L, Chiesa M, Colombo GI. Gene Regulatory Network Analysis of Perivascular Adipose Tissue of Abdominal Aortic Aneurysm Identifies Master Regulators of Key Pathogenetic Pathways. Biomedicines. 2020; 8(8):288. https://0-doi-org.brum.beds.ac.uk/10.3390/biomedicines8080288

Chicago/Turabian Style

Piacentini, Luca, Mattia Chiesa, and Gualtiero Ivanoe Colombo. 2020. "Gene Regulatory Network Analysis of Perivascular Adipose Tissue of Abdominal Aortic Aneurysm Identifies Master Regulators of Key Pathogenetic Pathways" Biomedicines 8, no. 8: 288. https://0-doi-org.brum.beds.ac.uk/10.3390/biomedicines8080288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop