Association of Pericardiac Adipose Tissue With Coronary Artery Disease

Background and Aim: Coronary artery disease (CAD) poses a worldwide health threat. Compelling evidence shows that pericardial adipose tissue (PAT), a brown-like adipose adjacent to the external surface of the pericardium, is associated with CAD. However, the specific molecular mechanisms of PAT in CAD are elusive. This study aims to characterize human PAT and explore its association with CAD.

Methods: We acquired samples of PAT from 31 elective cardiac surgery patients (17 CAD patients and 14 controls). The transcriptome characteristics were assessed in 5 CAD patients and 4 controls via RNA-sequencing. Cluster profile R package, String database, Cytoscape were applied to analyze the potential pathways and PPI-network key to DEGS, whereas the hubgenes were predicted via Metascape, Cytohubba, and MCODE. We use Cibersort, ENCORI, and DGIDB to predict immunoinfiltration, mRNA-miRNA target gene network, and search potential drugs targeting key DEGs. The predictable hubgenes and infiltrating inflammatory cells were validated in 22 patients (12 CAD samples and 10 control samples) through RT-qPCR and immunohistochemistry.

Results: A total of 147 different genes (104 up-regulated genes and 43 down-regulated genes) were identified in CAD patients. These different genes were associated with immunity and inflammatory dysfunction. Cibersort analysis showed monocytes and macrophages were the most common subsets in immune cells, whereas immunohistochemical results revealed there were more macrophages and higher proportion of M1 subtype cells in PAT of CAD patients. The PPI network and module analysis uncovered several crucial genes, defined as candidate genes, including Jun, ATF3, CXCR4, FOSB, CCl4, which were validated through RT-qPCR. The miRNA-mRNA network implicated hsa-miR-185-5p as diagnostic targets and drug-gene network showed colchicine, fenofibrate as potential therapeutic drugs, respectively.

Conclusion: This study demonstrates that PAT is mainly associated with the occurrence of CAD following the dysfunction of immune and inflammatory processes. The identified hubgenes, predicted drugs and miRNAs are promising biomarkers and therapeutic targets for CAD.

Introduction

Coronary artery disease (CAD) is a global health threat, particularly due to its high level of morbidity, which poses an enormous socioeconomic and medical burden (1). Obesity, a type of metabolic syndrome characterized by abnormal deposition of body fat with chronic inflammation of adipose tissue, is associated with multiple cardiovascular diseases, including CAD (2). According to the Framingham Heart Study, the risk factors for CAD are more associated with omental adipose tissue than subcutaneous adipose tissue (SAT) (3). These observations may be plausible owing to the differences in adipose tissue endocrine and metabolism. White adipose tissue (WAT) and brown adipose tissue (BAT) are the two major types of adipose tissues. Briefly, WAT comprises adipocytes with a large, single fat droplet and is presumed as the main depot for lipid storage, whereas BAT comprises several smaller fat droplets and numerous mitochondria and plays a role in heat production (4). In humans, BAT was thought to rapidly involute and essentially disappear within the first years after birth, only a small amount is found in the scapula, paraspinal, and around the heart and aorta in adulthood (5).

Pericardial adipose tissue (PAT), which refers to the fat surrounding the external surface of the pericardium, is supplied by the internal mammary artery (6). PAT covers the pericardium which is closely adjacent to the epicardial adipose tissue (EAT), coronary artery, and the heart. Scholars have suggested that PAT may play a vital role in cardiovascular disease (7–9). Previous reports in humans indicate that PAT may appear brown-like adipocyte in morphology, with distinct features different from WAT and BAT (10–12). Elsewhere, PAT, as a metabolically active endocrine local adipocyte depot, was found to be associated with coronary artery disease (CAD) through the production of free fatty acids and pro-and anti-inflammatory adipocytokines (13). Elevated PAT volume is known to be associated with coronary atherosclerosis, hypoadiponectinemia, inflammation and represents the highest risk factor for atherosclerosis (14). Other reports have further demonstrated the association of PAT with cardiovascular events and left ventricular remodeling (15, 16). However, whether the transcriptome of PAT changes during CAD, the molecular mechanism by which PAT mediates CAD progress, and the possibility to improve the function of PAT in CAD treatment remains elusive.

In the present study, we employed the RNA-sequencing (RNA-Seq) approach to explore the transcriptome characteristic shift in PAT from humans undergoing cardiac surgery with or without CAD. This was followed by the analysis of the functional enrichment, protein-protein interaction (PPI) network, hubgenes, and microRNA (miRNA)-mRNA regulatory network. We further predicted the potential drugs that target the key differentially expressed genes (DEGs). This study will deepen our understanding of PAT in the pathogenesis of CAD.

Study Patients

Subject Recruitment

This study complied with the Declaration of Helsinki and was approved by the ethics committee of Huadong Hospital Affiliated to Fudan University, Shanghai, China (2020K082). All patients signed written informed consent, and underwent preoperative coronary angiography. Control patients were referred for elective valve surgery and exhibited no significant CAD (a single lesion >50%) on preoperative coronary angiograms. Besides, CAD patients were referred for coronary artery bypass (CABG) surgery because of significant stenosis and surgical indications. Eventually, 31 patients, including 17 CAD patients undergoing coronary artery bypass grafting (CABG) and 14 control patients undergoing valve replacement or valve repair were enrolled for analysis.

Sample Collection and Preparation

We took a sample of PAT (1.0 g) adjacent to the pericardial surface during surgery (Figure 1A). The adipose sections collected from 5 CAD and 4 control patients were kept in 10% formalin for histological analysis. The other adipose tissues were immediately snap-frozen and stored at liquid nitrogen (–196°C). Five CAD and 4 control adipose tissue samples underwent RNA sequencing, while the rest were used in RT-PCR analyses.

FIGURE 1

Figure 1 (A) Localization of pericardial fat; (B) Heatmap results of DEGs; (C) Volcano plot results of DEGs; (D, E) The significantly enriched KEGG and GO terms that correspond to coding gene functions of upregulated and downregulated DEGs.

Histology and Immunohistochemistry

Formalin-fixed adipose tissues were embedded in paraffin and sectioned (8μm thick) and stained with hematoxylin and eosin (H&E). H&E-stained sections revealed the adipocyte size. To quantify M1 or M2 macrophages, cryosections were stained with rabbit anti-human CD11b (1:1000 dilution) (Abcam, UK) and rabbit anti-human CD68 (1:1000 dilution) (Abcam, UK) antibodies after which a goat anti-rabbit secondary antibody was conjugated to horseradish peroxidase. Immunohistochemical images were captured using a Zeiss LSM 5 Pascal microscope (M205FA, Zeiss, Oberkochen, Germany). Image-Pro Plus software (version6.0, MEDIA CYBERNETICS, USA) was employed to analyze the images.

mRNA Isolation and Real-Time PCR

Total RNA was isolated from adipose tissue biopsies using an RNA rapid extraction kit (TR205-200, Tianmo biotech, Beijing, China) following the manufacturer’s protocol. RT-PCR was performed using a cDNA reverse transcription kit (RR047A, TAKARA, Japan), Premix Ex Taq (RR420A, TAKARA, Japan) in 7500 Real-Time PCR system (Applied Biosystems, San Francisco, CA, USA). Standard and melting curves were generated in every plate for each gene to ensure that the reaction is efficient and specific. The cycle threshold value of β-actin acted as the internal control. The relative expression levels of different genes were analyzed via the 2−ΔΔCT method. Primer sequences were obtained from PrimerBank (http://pga.mgh.harvard.edu/ primerbank) as follows (Table 1): JUN forward TCCAAGTGCCGAAAAAGGAAG, JUN reverse CGAGTTCTGAGCTTTCAAGGT; ATF3 forward CGAGTTCTGAGCTTTCAAGGT, ATF3 reverse TTCTTTCTCGTCGCCTCTTTTT; CXCR4 forward ACTACACCGAGGAAATGGGCT, CXCR4 reverse CCCACAATGCCAGTTAAGAAGA; FOSB forward GCTGCAAGATCCCCTACGAAG, FOSB reverse ACGAAGAAGTGTACGAAGGGTT; CCL4 forward CTGTGCTGATCCCAGTGAATC, CCL4 reverse TCAGTTCAGTTCCAGGTCATACA; CXCL2 forward CTCAAGAATGGGCAGAAAGC, CXCL2 reverse CTCCTAAGTGATGCTCAAAC (17).

TABLE 1

Table 1 The primer sequences about hub genes.

RNA Sequencing

Paired-end libraries were synthesized using the TruSeq™ RNA Sample Preparation Kit (Illumina, USA) following the TruSeq™ RNA sample preparation guide. Briefly, the poly-A containing mRNA molecules were purified using poly-T oligo-attached magnetic beads. Thereafter, the mRNA was fragmented into small pieces using divalent cations under 94°C for 8 min. The cleaved RNA fragments were copied into first-strand cDNA using reverse transcriptase and random primers. Subsequently, the second strand cDNA was synthesized using DNA Polymerase I and RNase H. These cDNA fragments then underwent an end repair process, the addition of a single ‘A’ base, and ligation of the adapters. The products were purified and enriched with PCR to generate the final cDNA library. Clean libraries were quantified using a Qubit^® 2.0 Fluorometer (Life Technologies, USA), and validated by Agilent 2100 bioanalyzer (Agilent Technologies, USA) to confirm the insert size and evaluate the mole concentration. Cluster was generated using cBot with the library diluted to 10 pM and then was sequenced on the Illumina NovaSeq 6000 (Illumina, USA). Library construction and sequencing were performed by Sinotech Genomics Co., Ltd (Shanghai, China).

Data Processing of DEGs

Gene abundance was expressed as fragments per kilobase of exon per million reads mapped (FPKM). We employed the Stringtie software to count the fragments within each gene. TMM algorithm was applied for normalization. The DEGs between CAD and control patients were detected in the Illumina data collection software, whereby the P-value and |log₂FC| were calculated. Genes that met the cutoff criteria, P value<0.05 and |log₂FC|>1.0, considered as significantly modulated, were retained for subsequent analysis.

Enrichment Analysis

The Gene Ontology (GO) analysis and a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis were performed in the R package cluster Profiler (version 3.18.0) (18). GO enrichment analysis included biological processes (BP), cellular components (CC), and molecular functions (MF).

PPI Network Construction, Significant Module, and Hub Genes Analysis

The PPI network was first analyzed using Cytoscape (version 3.7.2, www.cytoscape.org) software after which the key genes in the PPI networks were identified using cytohubba (version 1.4.2), a plug-in of the Cytoscape software (19, 20). The Molecular Complex Detection tool (MCODE) (version 1.5.1) and Metascape (http://metascape.org/gp/) were employed to screen the significant module (21, 22).

Identification and Analysis of Significant Genes

A Venn diagram was delineated to identify significant common genes across “Metascape_MCODE”, “Cytoscape_MCODE”, and “Cytoscape_cytoHubba” by Veeny2.1 (https://bioinfogp.cnb.csic.es/tools/venny/). Summarized functions of significant genes were obtained via GeneCards (https://www.genecards.org/) (23).

Analyses of miRNA-mRNA Targets

We applied the miRNA-target tool ENCORI to predict the miRNA of DEGs (24). There were nine databases about miRNA-mRNA prediction. miRNAs predicted in at least two databases were selected as the potential target miRNAs of hubgenes. Then, the Cytoscape software was employed to assess the regulatory networks of the miRNA-mRNA pairs.

Immune Cell Infiltration

The immune cell components in adipose tissue were analyzed via CIBERSORT (25).

Prediction of Drugs Targeting DEGs

The Drug–Gene Interaction Database (DGIdb, www.dgidb.org) is a web resource that organizes and presents gene druggability information and drug-gene interactions from databases, articles, and web resources (26). Herein, we used DGIdb (version 3.0.2) to predict the potential drugs that target key DEGs confirmed via network module analysis. The following parameters were used: Preset filter, Food, and Drug Administration approved; advanced filters, source databases: all; gene categories, all; and interaction types, all. The interaction network was constructed using Cytoscape.

Statistics

Comparison of the clinical characteristics was achieved using the 2-tailed Student’s t-test for continuous variables or the χ2 test for dichotomous variables. A p-value is less than 0.05 denoted significance. Adipocyte size and positive cell numbers were compared via a 2-tailed Student’s t-test. For the qRT-PCR experiment, expression values relative to β-actin were compared using 2-way ANOVA and Tukey’s multiple comparisons test.

Results

Clinical Characteristics of Subjects

Table 2 outlines the clinical characteristics of 31 subjects. The mean age of the CAD group and control group were 66.06 ± 7.94 years and 60.93 ± 12.03 years, respectively (P = 0.165). Except for a trend towards more female patients in the control group (17.6% vs 42.9%, P = 0.124) and diabetes history in the CAD group (7/17 vs 2/14, P=0.101), we reported no difference in age, gender, BMI, comorbid conditions, and clinical biochemical characteristics. As such, our study groups were well matched. The clinical features of 9 cases who underwent RNA sequencing are outlined in Supplemental Table 1. Discriminant multivariate analysis demonstrated that the two groups were well matched.

TABLE 2

Table 2 Clinical characteristics of the subjects.

Depot-Specific Transcriptomic Profiles of PAT in CAD and Control Groups

RNA-seq was performed in 9 patients (CAD=5, Control=4), we first sought to identify depot-specific gene signatures in a pairwise manner. Using threshold criteria of fold change (FC) greater than or equal to 2 or less than or equal to –2 and p-value less than or equal to 0.05, 147 differentially expressed genes (DEGs) (104 upregulated and 43 downregulated) were identified in the CAD group relative to the control group (i.e., CAD vs. control groups) (Figures 1B, C). In addition, the top significantly changing genes were IL6, FOSB, CSF3, IGLC7, and PPP1R14C; these genes were associated with inflammation, nuclear transcription, cell differentiation, immunity, and neuronal activity.

Function and Pathway Enrichment Analyses

To explore the relative functions and pathways of DEGs in CAD vs. control groups, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis on those DEGs. The DEGs were significantly enriched in regulation of acute inflammatory response, regulation of complement activation, and complement activation classical pathway associated biological process (BP) terms, circulating immunoglobulin complex, immunoglobulin complex, and external side of plasma membrane associated component cell (CC) terms, antigen binding, RNA polymerase II-specific immunoglobulin receptor binding, and DNA-binding transcription activator activity associated molecule function (MF) terms (Figure 1E). KEGG pathway analysis demonstrated the enrichment of DEGs in cytokine-cytokine receptor interaction, IL-17 signaling pathway, and COVID-19 (Figure 1D).

PPI Network Construction and Module Analysis of DEGs

We systematically analyzed the biological functions of the obtained DEGs between the two groups using a PPI network of DEGs constructed via the STRING database and visualized using Cytoscape. The PPI network comprised 77 nodes (proteins) and 271 edges (interactions; Figure 2A). The 77 nodes included 71 up-regulated genes and 6 down-regulated genes. According to the topology score, Interleukin-6 (IL-6), Jun proto-oncogene (JUN), and prostaglandin-endoperoxide synthase 2 (PTGS2) activating transcription factor 3 (ATF3), chemokine (C-X-C motif) receptor 4 (CXCR4) were the top five genes (Figure 2B). Additionally, five MCODE modules were identified from the PPI network via cytoscape_CMODE; however, one submodule with a score >5 was extracted from the PPI network, which comprised 8 nodes and 26 edges (Figure 2C). In this module, IL-6, chemokine (C-X-C motif) ligand 2 (CXCL2), selectin E (SELE), C-C motif chemokine ligand 4 (CCL4), PTGS2 were the top five genes (Figure 2C). Three MCODE modules were extracted through Metascape analysis (Figure 2D).

FIGURE 2

Figure 2 (A) The protein-protein interaction network constructed based on the DEGs. The green heart represents the downregulated genes, the red parallelogram represents the upregulated genes. (B–D) The hubgenes predicted by cytohubba, Mcode, and metascape Mcode, respectively. (E) Overlap of hubgenes identified via the three methods; (F) The Pathway analysis results showing the six enriched hubgenes.

Identification and Analysis of Significant Genes

To reveal the most important hub gene, we filtered these genes using the VENN diagram. The VENN diagram revealed six significant common genes, including JUN, ATF3, CXCR4, FOSB, CCL4, and CXCL2 (Figure 2E). The functions of the six significant genes are summarized in Table 3. We shifted our focus to explore the potential function of these six key genes. Pathway analysis demonstrated that the six genes were mainly enriched in the chemokine-mediated signaling pathway, inflammatory response, transcription from RNA polymerase II promoter, G-protein coupled receptor signaling pathway, and positive regulation of transcription from RNA polymerase II promoter (Figure 2F). RT-qPCR analysis revealed that the relative expression levels of JUN (Figure 3A), FOSB (Figure 3B), ATF3 (Figure 3C), CCL4 (Figure 3D), and CXCR4 (Figure 3E) were significantly higher in the CAD group adipose tissue than in control group.

TABLE 3

Table 3 The enrichment results of GO and KEGG of hubgenes.

FIGURE 3

Figure 3 (A–E) RT-PCR verification of the expression levels of 5 genes (JUN, FOSB, ATF3, CCL4, CXCR4). The 2−ΔΔCT method was used to analyze the relative expression levels of various genes. * means 0.05 < p<0.1; ** means p < 0.05; *** means p < 0.01; **** means p < 0.001.

miRNA-mRNA Interaction Network

We explored the effect of PAT on CAD progression and its potential gene regulation mechanism. miRNA-target gene interaction pairs of reverse association were predicted via ENCORI according to the hub genes identified previously. Considering the identified miRNA-mRNA pairs, we compared the interaction network containing 70 miRNA-mRNA pairs and visualized them via the Cytoscape software. Through comparison of the targets of hub genes, we found CXCR4 to be a potential target of 23 miRNAs, including hsa-miR-185-5p. Also, FUN and FOSB were the potential targets of 15 miRNAs. The miRNA-gene regulation network is illustrated in Table 4 and Figure 4, respectively.

TABLE 4

Table 4 The predicted miRNA about 6 hubgenes.

FIGURE 4

Figure 4 The constructed mRNA-miRNA network. The green heart represents the downregulated genes, the red diamond represents the hubgenes.

Macrophages Are the Major Immune Infiltrating Cell Subset

GO and KEGG enrichment analysis has shown inflammatory response is the major procession in PAT of CAD patients. However, the component of immune cells in the pericardial adipose deposit is unknown. We employed CIBERSORT, a bioinformatics tool used to infer immune cell composition from RNA-seq datasets, to compute the relative frequency of 22 infiltrating immune cell subsets in the 9 cases. Results revealed that monocytes and macrophages were the most common immune cell subsets with mean fractions of 0.098 and 0.325, respectively (Figure 5). Further, we verified the inflammatory phenotype of human PAT from both CAD and control patients through immunohistochemistry analysis of markers of T cell (CD3), macrophage (CD68), and its M1 subtype (CD11b). We found that PAT from CAD patients had significantly more T cells (CD3+) and macrophages (CD68+) than that of control patients (Figure 6). Semi-quantification analysis revealed that the IOD/Area of macrophages and T-cells in PAT of CAD patients was significantly different (Table 5). These findings further proved that PAT exerts a potential inflammatory effect in the regulation of CAD.

FIGURE 5

Figure 5 Distinction of infiltrating immune cell subpopulations and levels between CHD/NCHD groups.

FIGURE 6

Figure 6 IHC staining of CD68, CD11b, and CD3. CD68, CD11b, and CD3 protein are stained in brown.

TABLE 5

Table 5 The semi-quantification analysis of infiltrating immune cell.

Drug Predictable Results

The search for potential drugs that can improve the function of PAT to intervene in the disease process of CAD propelled us to analyze the six key genes for potential drugs via DGIdb. Finally, five genes were considered as druggable genes, including JUN, CXCL2, and CXCR4; also we obtained 77 possible drugs. An interaction network including 82 nodes and 77 edges was constructed based on the five druggable genes (Figure 7).

FIGURE 7

Figure 7 The constructed mRNA-drug network. The purple oval represents the genes, and the light red heart represents the drugs.

Discussion

In this study, we have sequenced the transcriptome of PAT from five patients with CAD and four patients without CAD and analyzed the DEGs between two groups. A total of 147 DEGs (104 up-regulated genes and 43 down-regulated genes) were identified. We subjected the DEGs to GO functional and KEGG pathway enrichment analyses via Clusterprofile package, and obtained 40 significant GO terms and 55 significant KEGG pathways. GO term also demonstrated the interaction of these genes with immune and inflammatory cells, and we, therefore, employed the Cibersort online tool (https://cibersort.stanford.edu/) to analyze 22 types of immune cells infiltration. Based on the results, monocytes and macrophages were the most common immune cell subsets. Further immunohistochemical staining demonstrated that more macrophages were present in PAT from CAD patients, and the proportion of M1 subtype cells was higher. These findings imply that PAT inflammation potentially plays a crucial role in CAD progression.

Moreover, GO and KEGG pathway analysis results revealed the significant enrichment of DEGs in cytokine-cytokine receptor interaction, IL-17 signaling pathway, and COVID-19, all of which are inflammatory processes. Compelling evidence implicates the IL-17 signaling pathway in the pathogenesis of many inflammatory-driven conditions, including atherosclerosis (27–30). A study found that IL-17 maintained plaque stability through induction of proliferation of SMCs and collagen content in atherosclerotic plaques (29). Also, IL-17 had been revealed to potentially lower the expression of vascular cell adhesion molecule (VCAM)-1 in endothelial cells and prevent monocyte adherence and block T cell infiltration into plaques (31). Other reports show that IL-17 can induce the release of chemokines, including CXCL1, CXCL2, CXCL8, CXCL10. Subsequently, these chemokines can recruit neutrophils and monocytes to the atherosclerotic lesion (32, 33). Furthermore, IL-17 can stimulate macrophages to produce inflammatory cytokines, such as IL-6, TNF-a, and IL-1b (34–36), and induce apoptosis of vascular endothelial cells through activation of caspase-3 and caspase-9. These observations affirm that the IL-17 signaling pathway mediates atherosclerosis (37). We also found that other signaling pathways such as NF-kappa B signaling pathway, TNF signaling pathway, JAK-STAT signaling pathway, and Toll-like receptor signaling pathway were enriched in PAT of CAD patients. Of note, some proteins in these pathways can interplay with each other, thereby increasing the severity of inflammation and apoptosis (38). These results imply that the management of inflammation in CAD may require nonspecific inflammatory inhibitors.

The PPI network-integrated three modules in the analysis of six significant and reproducible genes (JUN, ATF3, CXCR4, FOSB, CCL4, CXCL2), and revealed differential expression between the CAD and control groups. Three core genes (CXCR4, CXCL2, and CCL4) were markedly enriched in the chemokine signaling pathway, whereas the others (Jun, FosB, and ATF3) were transcription factors. Chemokines are cytokines that mediate cell chemotaxis and stagnation as they bind to their respective cell surface receptors (39). Chemokines induce the aggregation of inflammatory cells to the inflammatory site, which is the main factor that trigger vascular inflammatory injury (39, 40). Previous reports indicate that CCL4 is highly expressed in atherosclerotic patients (41), and it potentially exerts a crucial role in the progression of atherosclerotic plaque (42–44). In addition, inhibition of CCL4 can stabilize atherosclerotic plaques by decreasing the expression of MMPs, inflammatory cell infiltration, circulation of inflammatory factors, and regulation of blood lipid levels (43). The C-X-C motif chemokine receptor 4 (CXCR4) is a cytokine receptor and mediates various biological processes (45). The ligand of CXCR4 is C-X-C motif chemokine ligand 12 (CXCL12) (46). CXCL12 and CXCR4 coordinatively play a pivotal role in atherosclerosis and arterial injury (47, 48). CXCR4 expression in leukocytes is closely associated with the vulnerability of atherosclerosis plaque (49). A recent study found that CXCR4-positive macrophages accumulated in tissue samples of human carotid plaques and, CXCR4 was expressed in both smooth muscle cell progenitors and endothelial progenitor cells in atherosclerotic plaque progression (46). The study by Puca AA showed that upregulation of CXCR4 decreases the development of atherogenic process. This is because it can skew macrophages to acquire an M2-resolving phenotype, maintains maintain arterial integrity, preserves preserve endothelial cell integrity, and restore the normal contractile SMC phenotype (50, 51). In contrast to CXCR4, CCL4 inhibition reduced the adhesiveness of coronary endothelial cells, which is an early sign of atherogenesis (43). Herein, we found that CCL4, CXCR4 derived from PAT was highly expressed at the mRNA level in CAD patients. More exploration of CCL4 and CXCR4 secreted by PAT would enhance the understanding of the pathophysiological mechanism of atherosclerosis and facilitate its utilization as a biomarker and intervention target for atherosclerosis.

Activator protein-1 (AP-1) is an important class of nuclear transcription factors in the body. It is a homodimer or heterodimer composed of JUN, FOS, ATF, and MAF protein families which exert biological effects (52). AP-1 is largely associated with cell proliferation, differentiation, apoptosis, and inflammation (53–55). Previously, the AP-1 cascade was found to potentially drive additional leukocyte recruitment to the developing atherosclerotic matrix (56). AP-1 could also promote the development of atherosclerosis by inducing endothelial cell death, the proliferation and migration of vascular smooth muscle cells (57, 58). In this study, we observed a significant increase in mRNA expression of AP-1 factors JUN, FOSB, and ATF3 from PAT in CAD patients. These results imply that PAT-derived AP-1 has potential effects on the progression of the coronary artery through cell proliferation, inflammation, and extracellular matrix remodeling.

Based on bioinformatics, we constructed the miRNA-mRNA network of core genes through the ENCORI database to further elucidate the regulatory mechanism of core genes. Among the predicted miRNAs, only miRNA-185-5p bound to three of the six core genes (FOSB, CCL4, and CXCR4). Previous studies had confirmed that miRNA-185-5p is widely involved in the proliferation, metastasis, and inflammation of various tumors (59–61). However, its regulatory role in inflammation and the progression of CAD remains elusive.

Through the prediction of the Drug-Gene Interaction database, we revealed some drugs that potentially interact with the core genes. Some drugs were used in the clinic but not to manage CAD. For instance, colchicine, an anti-inflammatory alkaloid, has been mainly used in gout diseases (62); it exerts inhibitory effects on leukocyte chemotaxis and adhesion, microtubule assembly, and reduces the production of inflammatory mediators (63, 64). Because of its anti-inflammatory properties, the feasibility of the value of colchicine in CAD had been investigated in several clinical trials. Most clinical trials suggest that colchicine is beneficial in patients with CAD or myocardial infarction (65–67). However, some studies revealed that colchicine could reduce major cardiovascular events (68), however, the overall mortality did not decrease. As such, colchicine may increase the mortality of non-cardiovascular events and limit its application in CAD (69), which warrants more large-scale clinical trials. Fenofibrate is a broad-spectrum lipid-lowering drug, which acts on PPAR to potentially reduce cholesterol and triglyceride (70, 71). Elsewhere, a study found that fibrates have potential anti-inflammatory effects, protect endothelial cell function, and improve insulin resistance (72). However, these drugs are not recommended for the secondary prevention of CAD. In clinical trials such as HHS, VA-HIT, and BIP, fibrates have proved to exhibit significant benefits in subgroups of patients, including those with insulin resistance or metabolic syndrome (73–75). Recent studies have further outlined that fibrates can reduce the risk of cardiovascular events and the level of uric acid in diabetic patients with dyslipidemia (76, 77). Collectively, fibrates may have great clinical prospects in some specific groups of patients with CAD.

In summary, this work has allowed for the identification of the transcriptome characteristic shift in PAT from humans with or without CAD and revealed the upregulation in inflammatory processes of PAT in CAD patients. The crucial genes, pathways, and drug target genes closely associated with CAD have been revealed through bioinformatics analyses. Particularly, the critical genes are potential biomarkers and therapeutic targets in CAD. Also, colchicine and fenofibrate are the most promising drugs for coronary artery disease.

Limitations of the Study

The primary limitation of the research is that we did not analyze samples from healthy subjects; this could have been more interesting for the PAT physiology study. Besides, the analyzed samples from older patients with some comorbidity who came in for various cardiac surgeries rather than no comorbidity. We also acknowledge that we had a rather small RNA-Seq sample size, which could not allow us to analyze samples grouped by obesity, age, and gender. This may have contributed to the CAD-associated transcriptomic shifts.

Data Availability Statement

The datasets presented in this study can be found in GEO (Gene Expression Ominbus), a public functional genomics data repository. The accession number(s) in the GEO can be found below: GSE179397.

Ethics Statement

The studies involving human participants were reviewed and approved by the ethics committee of Huadong Hospital Affiliated to Fudan University, Shanghai, China. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

ML and LQ have contributed equally to this work and share first authorship. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Science and Technology Commission of Shanghai Municipality (20140900600), National Natural Science Foundation of China (81770420, 91949126, 81970378), Shanghai Municipal Key Clinical Specialty (shslczdzk02801), National Key Research and Development Program of China (2020YFC2009001), and Center of geratic coronary artery disease.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We are grateful for the patients and their surgeons.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.724859/full#supplementary-material

References

1. Jemal A, Ward E, Hao Y, Thun M. Trends in the Leading Causes of Death in the United States, 1970-2002. JAMA (2005) 294(10):1255–9. doi: 10.1001/jama.294.10.1255

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Van Gaal LF, Mertens IL, De Block CE. Mechanisms Linking Obesity With Cardiovascular Disease. Nature (2006) 444(7121):875–80. doi: 10.1038/nature05487

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Fox CS, Massaro JM, Hoffmann U, Pou KM, Maurovich-Horvat P, Liu CY, et al. Abdominal Visceral and Subcutaneous Adipose Tissue Compartments: Association With Metabolic Risk Factors in the Framingham Heart Study. Circulation (2007) 116(1):39–48. doi: 10.1161/CIRCULATIONAHA.106.675355

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Fruhbeck G, Becerril S, Sainz N, Garrastachu P, Garcia-Velloso MJ. BAT: A New Target for Human Obesity? Trends Pharmacol Sci (2009) 30(8):387–96. doi: 10.1016/j.tips.2009.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Nedergaard J, Bengtsson T, Cannon B. Unexpected Evidence for Active Brown Adipose Tissue in Adult Humans. Am J Physiol Endocrinol Metab (2007) 293(2):E444–52. doi: 10.1152/ajpendo.00691.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Iacobellis G, Corradi D, Sharma AM. Epicardial Adipose Tissue: Anatomic, Biomolecular and Clinical Relationships With the Heart. Nat Clin Pract Cardiovasc Med (2005) 2(10):536–43. doi: 10.1038/ncpcardio0319

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Ding J, Hsu FC, Harris TB, Liu Y, Kritchevsky SB, Szklo M, et al. The Association of Pericardial Fat With Incident Coronary Heart Disease: The Multi-Ethnic Study of Atherosclerosis (MESA). Am J Clin Nutr (2009) 90(3):499–504. doi: 10.3945/ajcn.2008.27358

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kim TH, Yu SH, Choi SH, Yoon JW, Kang SM, Chun EJ, et al. Pericardial Fat Amount is an Independent Risk Factor of Coronary Artery Stenosis Assessed by Multidetector-Row Computed Tomography: The Korean Atherosclerosis Study 2. Obes (Silver Spring) (2011) 19(5):1028–34. doi: 10.1038/oby.2010.246

CrossRef Full Text | Google Scholar

9. Yang FS, Yun CH, Wu TH, Hsieh YC, Bezerra HG, Liu CC, et al. High Pericardial and Peri-Aortic Adipose Tissue Burden in Pre-Diabetic and Diabetic Subjects. BMC Cardiovasc Disord (2013) 13:98. doi: 10.1186/1471-2261-13-98

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Cheung L, Gertow J, Werngren O, Folkersen L, Petrovic N, Nedergaard J, et al. Human Mediastinal Adipose Tissue Displays Certain Characteristics of Brown Fat. Nutr Diabetes (2013) 3:e66. doi: 10.1038/nutd.2013.6

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sacks HS, Fain JN, Bahouth SW, Ojha S, Frontini A, Budge H, et al. Adult Epicardial Fat Exhibits Beige Features. J Clin Endocrinol Metab (2013) 98(9):E1448–55. doi: 10.1210/jc.2013-1265

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sacks HS, Fain JN, Holman B, Cheema P, Chary A, Parks F, et al. Uncoupling Protein-1 and Related Messenger Ribonucleic Acids in Human Epicardial and Other Adipose Tissues: Epicardial Fat Functioning as Brown Fat. J Clin Endocrinol Metab (2009) 94(9):3611–5. doi: 10.1210/jc.2009-0571

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Yamada H, Sata M. Role of Pericardial Fat: The Good, the Bad and the Ugly. J Cardiol (2015) 65(1):2–4. doi: 10.1016/j.jjcc.2014.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Greif M, Becker A, von Ziegler F, Lebherz C, Lehrke M, Broedl UC, et al. Pericardial Adipose Tissue Determined by Dual Source CT is a Risk Factor for Coronary Atherosclerosis. Arterioscler Thromb Vasc Biol (2009) 29(5):781–6. doi: 10.1161/ATVBAHA.108.180653

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Shah RV, Anderson A, Ding J, Budoff M, Rider O, Petersen SE, et al. Pericardial, But Not Hepatic, Fat by CT is Associated With CV Outcomes and Structure: The Multi-Ethnic Study of Atherosclerosis. JACC Cardiovasc Imaging (2017) 10(9):1016–27. doi: 10.1016/j.jcmg.2016.10.024

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Cai S, Wald R, Deva DP, Kiaii M, Ng MY, Karur GR, et al. Cardiac MRI Measurements of Pericardial Adipose Tissue Volumes in Patients on in-Centre Nocturnal Hemodialysis. J Nephrol (2020) 33(2):355–63. doi: 10.1007/s40620-019-00665-4

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wang X, Spandidos A, Wang H, Seed B. Primerbank: A PCR Primer Database for Quantitative Gene Expression Analysis, 2012 Update. Nucleic Acids Res (2012) 40(Database issue):D1144–9. doi: 10.1093/nar/gkr1013

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yu G, Wang LG, Han Y, He QY. Clusterprofiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res (2003) 13(11):2498–504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. Cytohubba: Identifying Hub Objects and Sub-Networks From Complex Interactome. BMC Syst Biol (2014) 8 Suppl 4:S11. doi: 10.1186/1752-0509-8-S4-S11

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bader GD, Hogue CW. An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinf (2003) 4:2. doi: 10.1186/1471-2105-4-2

CrossRef Full Text | Google Scholar

22. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D. Genecards: A Novel Functional Genomics Compendium With Automated Data Mining and Query Reformulation Support. Bioinformatics (1998) 14(8):656–64. doi: 10.1093/bioinformatics/14.8.656

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Li JH, Liu S, Zhou H, Qu LH, Yang JH. Starbase V2.0: Decoding MiRNA-CeRNA, MiRNA-NcRNA and Protein-RNA Interaction Networks From Large-Scale CLIP-Seq Data. Nucleic Acids Res (2014) 42(Database issue):D92–7. doi: 10.1093/nar/gkt1248

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust Enumeration of Cell Subsets From Tissue Expression Profiles. Nat Methods (2015) 12(5):453–7. doi: 10.1038/nmeth.3337

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Cotto KC, Wagner AH, Feng YY, Kiwala S, Coffman AC, Spies G, et al. Dgidb 3.0: A Redesign and Expansion of the Drug-Gene Interaction Database. Nucleic Acids Res (2018) 46(D1):D1068–D73. doi: 10.1093/nar/gkx1143

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Amara S, Lopez K, Banan B, Brown SK, Whalen M, Myles E, et al. Synergistic Effect of Pro-Inflammatory Tnfalpha and IL-17 in Periostin Mediated Collagen Deposition: Potential Role in Liver Fibrosis. Mol Immunol (2015) 64(1):26–35. doi: 10.1016/j.molimm.2014.10.021

PubMed Abstract | CrossRef Full Text | Google Scholar

28. de Boer OJ, van der Meer JJ, Teeling P, van der Loos CM, Idu MM, van Maldegem F, et al. Differential Expression of Interleukin-17 Family Cytokines in Intact and Complicated Human Atherosclerotic Plaques. J Pathol (2010) 220(4):499–508. doi: 10.1002/path.2667

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Gistera A, Robertson AK, Andersson J, Ketelhuth DF, Ovchinnikova O, Nilsson SK, et al. Transforming Growth Factor-Beta Signaling in T Cells Promotes Stabilization of Atherosclerotic Plaques Through an Interleukin-17-Dependent Pathway. Sci Transl Med (2013) 5(196):196ra00. doi: 10.1126/scitranslmed.3006133

CrossRef Full Text | Google Scholar

30. Yazdani MR, Khosropanah S, Doroudchi M. Interleukin-17 Production by CD4+CD45RO+Foxp3+ T Cells in Peripheral Blood of Patients With Atherosclerosis. Arch Med Sci Atheroscler Dis (2019) 4:e215–e24. doi: 10.5114/amsad.2019.87525

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Taleb S, Romain M, Ramkhelawon B, Uyttenhove C, Pasterkamp G, Herbin O, et al. Loss of SOCS3 Expression in T Cells Reveals a Regulatory Role for Interleukin-17 in Atherosclerosis. J Exp Med (2009) 206(10):2067–77. doi: 10.1084/jem.20090545

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Rao DA, Eid RE, Qin L, Yi T, Kirkiles-Smith NC, Tellides G, et al. Interleukin (IL)-1 Promotes Allogeneic T Cell Intimal Infiltration and IL-17 Production in a Model of Human Artery Rejection. J Exp Med (2008) 205(13):3145–58. doi: 10.1084/jem.20081661

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Weaver CT, Hatton RD. Interplay Between the TH17 and Treg Cell Lineages: A (Co-)Evolutionary Perspective. Nat Rev Immunol (2009) 9(12):883–9. doi: 10.1038/nri2660

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Erbel C, Chen L, Bea F, Wangler S, Celik S, Lasitschka F, et al. Inhibition of IL-17A Attenuates Atherosclerotic Lesion Development in Apoe-Deficient Mice. J Immunol (2009) 183(12):8167–75. doi: 10.4049/jimmunol.0901126

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Chen S, Crother TR, Arditi M. Emerging Role of IL-17 in Atherosclerosis. J Innate Immun (2010) 2(4):325–33. doi: 10.1159/000314626

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ng HP, Burris RL, Nagarajan S. Attenuated Atherosclerotic Lesions in Apoe-Fcgamma-Chain-Deficient Hyperlipidemic Mouse Model is Associated With Inhibition of Th17 Cells and Promotion of Regulatory T Cells. J Immunol (2011) 187(11):6082–93. doi: 10.4049/jimmunol.1004133

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhu F, Wang Q, Guo C, Wang X, Cao X, Shi Y, et al. IL-17 Induces Apoptosis of Vascular Endothelial Cells: A Potential Mechanism for Human Acute Coronary Syndrome. Clin Immunol (2011) 141(2):152–60. doi: 10.1016/j.clim.2011.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

38. McKellar GE, McCarey DW, Sattar N, McInnes IB. Role for TNF in Atherosclerosis? Lessons From Autoimmune Disease. Nat Rev Cardiol (2009) 6(6):410–7. doi: 10.1038/nrcardio.2009.57

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Blanchet X, Langer M, Weber C, Koenen RR, von Hundelshausen P. Touch of Chemokines. Front Immunol (2012) 3:175. doi: 10.3389/fimmu.2012.00175

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Mikolajczyk TP, Szczepaniak P, Vidler F, Maffia P, Graham GJ, Guzik TJ. Role of Inflammatory Chemokines in Hypertension. Pharmacol Ther (2020) 223:107799. doi: 10.1016/j.pharmthera.2020.107799

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Cagnin S, Biscuola M, Patuzzo C, Trabetti E, Pasquali A, Laveder P, et al. Reconstruction and Functional Analysis of Altered Molecular Pathways in Human Atherosclerotic Arteries. BMC Genomics (2009) 10:13. doi: 10.1186/1471-2164-10-13

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Osonoi Y, Mita T, Azuma K, Nakajima K, Masuyama A, Goto H, et al. Defective Autophagy in Vascular Smooth Muscle Cells Enhances Cell Death and Atherosclerosis. Autophagy (2018) 14(11):1991–2006. doi: 10.1080/15548627.2018.1501132

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Chang TT, Yang HY, Chen C, Chen JW. CCL4 Inhibition in Atherosclerosis: Effects on Plaque Stability, Endothelial Cell Adhesiveness, and Macrophages Activation. Int J Mol Sci (2020) 21(18):6567. doi: 10.3390/ijms21186567

CrossRef Full Text | Google Scholar

44. Munjal A, Khandia R. Atherosclerosis: Orchestrating Cells and Biomolecules Involved in its Activation and Inhibition. Adv Protein Chem Struct Biol (2020) 120:85–122. doi: 10.1016/bs.apcsb.2019.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Hattermann K, Mentlein R. An Infernal Trio: The Chemokine CXCL12 and its Receptors CXCR4 and CXCR7 in Tumor Biology. Ann Anat (2013) 195(2):103–10. doi: 10.1016/j.aanat.2012.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Merckelbach S, van der Vorst EPC, Kallmayer M, Rischpler C, Burgkart R, Doring Y, et al. Expression and Cellular Localization of CXCR4 and CXCL12 in Human Carotid Atherosclerotic Plaques. Thromb Haemost (2018) 118(1):195–206. doi: 10.1160/TH17-04-0271

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Sun K, Xiang X, Li N, Huang S, Qin X, Wu Y, et al. Gene-Diet Interaction Between SIRT6 and Soybean Intake for Different Levels of Pulse Wave Velocity. Int J Mol Sci (2015) 16(7):14338–52. doi: 10.3390/ijms160714338

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Guo Y, Garcia-Barrio MT, Wang L, Chen YE. Experimental Biology for the Identification of Causal Pathways in Atherosclerosis. Cardiovasc Drugs Ther (2016) 30(1):1–11. doi: 10.1007/s10557-016-6644-7

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Doring Y, Jansen Y, Cimen I, Aslani M, Gencer S, Peters LJF, et al. B-Cell-Specific CXCR4 Protects Against Atherosclerosis Development and Increases Plasma Igm Levels. Circ Res (2020) 126(6):787–8. doi: 10.1161/CIRCRESAHA.119.316142

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Puca AA, Carrizzo A, Spinelli C, Damato A, Ambrosio M, Villa F, et al. Single Systemic Transfer of a Human Gene Associated With Exceptional Longevity Halts the Progression of Atherosclerosis and Inflammation in Apoe Knockout Mice Through a CXCR4-Mediated Mechanism. Eur Heart J (2020) 41(26):2487–97. doi: 10.1093/eurheartj/ehz459

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Doring Y, Noels H, van der Vorst EPC, Neideck C, Egea V, Drechsler M, et al. Vascular CXCR4 Limits Atherosclerosis by Maintaining Arterial Integrity: Evidence From Mouse and Human Studies. Circulation (2017) 136(4):388–403. doi: 10.1161/CIRCULATIONAHA.117.027646

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Eferl R, Wagner EF. AP-1: A Double-Edged Sword in Tumorigenesis. Nat Rev Cancer (2003) 3(11):859–68. doi: 10.1038/nrc1209

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Lee JG, Heur M. Interleukin-1beta Enhances Cell Migration Through AP-1 and NF-KappaB Pathway-Dependent FGF2 Expression in Human Corneal Endothelial Cells. Biol Cell (2013) 105(4):175–89. doi: 10.1111/boc.201200077

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Jeon Y, Jung Y, Kim MC, Kwon HC, Kang KS, Kim YK, et al. Sargahydroquinoic Acid Inhibits Tnfalpha-Induced AP-1 and NF-Kappab Signaling in Hacat Cells Through Pparalpha Activation. Biochem Biophys Res Commun (2014) 450(4):1553–9. doi: 10.1016/j.bbrc.2014.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Hellweg CE, Spitta LF, Henschenmacher B, Diegeler S, Baumstark-Khan C. Transcription Factors in the Cellular Response to Charged Particle Exposure. Front Oncol (2016) 6:61. doi: 10.3389/fonc.2016.00061

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Adhikari N, Charles N, Lehmann U, Hall JL. Transcription Factor and Kinase-Mediated Signaling in Atherosclerosis and Vascular Injury. Curr Atheroscler Rep (2006) 8(3):252–60. doi: 10.1007/s11883-006-0081-1

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Ji R, Gu Y, Zhang J, Gao C, Gao W, Zang X, et al. TRIM7 Promotes Proliferation and Migration of Vascular Smooth Muscle Cells in Atherosclerosis Through Activating C-Jun/AP-1. IUBMB Life (2020) 72(2):247–58. doi: 10.1002/iub.2181

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Valente AJ, Irimpen AM, Siebenlist U, Chandrasekar B. Oxldl Induces Endothelial Dysfunction and Death via TRAF3IP2: Inhibition by HDL3 and AMPK Activators. Free Radic Biol Med (2014) 70:117–28. doi: 10.1016/j.freeradbiomed.2014.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Zhao J, Jiang Y, Zhang H, Zhou J, Chen L, Li H, et al. The SRSF1/Circatp5b/MiR-185-5p/HOXB5 Feedback Loop Regulates the Proliferation of Glioma Stem Cells via the IL6-Mediated JAK2/STAT3 Signaling Pathway. J Exp Clin Cancer Res (2021) 40(1):134. doi: 10.1186/s13046-021-01931-9

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Liu JQ, Deng M, Xue NN, Li TX, Guo YX, Gao L, et al. LncRNA KLF3-AS1 Suppresses Cell Migration and Invasion in ESCC by Impairing Mir-185-5p-Targeted KLF3 Inhibition. Mol Ther Nucleic Acids (2020) 20:231–41. doi: 10.1016/j.omtn.2020.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Sun CC, Zhang L, Li G, Li SJ, Chen ZL, Fu YF, et al. The Lncrna PDIA3P Interacts With MiR-185-5p to Modulate Oral Squamous Cell Carcinoma Progression by Targeting Cyclin D2. Mol Ther Nucleic Acids (2017) 9:100–10. doi: 10.1016/j.omtn.2017.08.015

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Dudkiewicz I, Brosh T, Perelman M, Salai M. Colchicine Inhibits Fracture Union and Reduces Bone Strength–In Vivo Study. J Orthop Res (2005) 23(4):877–81. doi: 10.1016/j.orthres.2004.11.014

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Deftereos S, Giannopoulos G, Papoutsidakis N, Panagopoulou V, Kossyvakis C, Raisakis K, et al. Colchicine and the Heart: Pushing the Envelope. J Am Coll Cardiol (2013) 62(20):1817–25. doi: 10.1016/j.jacc.2013.08.726

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Imazio M, Andreis A, Brucato A, Adler Y, De Ferrari GM. Colchicine for Acute and Chronic Coronary Syndromes. Heart (2020) 106(20):1555–60. doi: 10.1136/heartjnl-2020-317108

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Roubille F, Tardif JC. Colchicine for Secondary Cardiovascular Prevention in Coronary Disease. Circulation (2020) 142(20):1901–4. doi: 10.1161/CIRCULATIONAHA.120.051240

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Slomski A. Colchicine Reduces Cardiovascular Events in Chronic Coronary Disease. JAMA (2020) 324(16):1599. doi: 10.1001/jama.2020.20646

CrossRef Full Text | Google Scholar

67. Bates ER. In Patients With CAD, Adding Colchicine to Evidence-Based Therapies Reduces MACE Without Increasing Adverse Outcomes. Ann Intern Med (2021) 174(3):JC30. doi: 10.7326/ACPJ202103160-030

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Tong DC, Quinn S, Nasis A, Hiew C, Roberts-Thomson P, Adams H, et al. Colchicine in Patients With Acute Coronary Syndrome: The Australian COPS Randomized Clinical Trial. Circulation (2020) 142(20):1890–900. doi: 10.1161/CIRCULATIONAHA.120.050771

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Fiolet ATL, Opstal TSJ, Mosterd A, Eikelboom JW, Jolly SS, Keech AC, et al. Efficacy and Safety of Low-Dose Colchicine in Patients With Coronary Disease: A Systematic Review and Meta-Analysis of Randomized Trials. Eur Heart J (2021) 42(28):2765–75. doi: 10.1093/eurheartj/ehab115

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Maciejewski S, Hilleman D. Effectiveness of a Fenofibrate 145-Mg Nanoparticle Tablet Formulation Compared With the Standard 160-Mg Tablet in Patients With Coronary Heart Disease and Dyslipidemia. Pharmacotherapy (2008) 28(5):570–5. doi: 10.1592/phco.28.5.570

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Packard CJ. Overview of Fenofibrate. Eur Heart J (1998) 19 Suppl A:A62–5.

PubMed Abstract | Google Scholar

72. Balakumar P, Sambathkumar R, Mahadevan N, Muhsinah AB, Alsayari A, Venkateswaramurthy N, et al. Molecular Targets of Fenofibrate in the Cardiovascular-Renal Axis: A Unifying Perspective of its Pleiotropic Benefits. Pharmacol Res (2019) 144:132–41. doi: 10.1016/j.phrs.2019.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Rubins HB, Robins SJ, Collins D, Fye CL, Anderson JW, Elam MB, et al. Gemfibrozil for the Secondary Prevention of Coronary Heart Disease in Men With Low Levels of High-Density Lipoprotein Cholesterol. Veterans Affairs High-Density Lipoprotein Cholesterol Intervention Trial Study Group. N Engl J Med (1999) 341(6):410–8. doi: 10.1056/NEJM199908053410604

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Frick MH, Elo O, Haapa K, Heinonen OP, Heinsalmi P, Helo P, et al. Helsinki Heart Study: Primary-Prevention Trial With Gemfibrozil in Middle-Aged Men With Dyslipidemia. Safety of Treatment, Changes in Risk Factors, and Incidence of Coronary Heart Disease. N Engl J Med (1987) 317(20):1237–45. doi: 10.1056/NEJM198711123172001

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Bezafibrate Infarction Prevention (BIP) Study. Secondary Prevention by Raising HDL Cholesterol and Reducing Triglycerides in Patients With Coronary Artery Disease. Circulation (2000) 102(1):21–7. doi: 10.1161/01.cir.102.1.21

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Waldman B, Ansquer JC, Sullivan DR, Jenkins AJ, McGill N, Buizen L, et al. Effect of Fenofibrate on Uric Acid and Gout in Type 2 Diabetes: A Post-Hoc Analysis of the Randomised, Controlled FIELD Study. Lancet Diabetes Endocrinol (2018) 6(4):310–8. doi: 10.1016/S2213-8587(18)30029-9

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Kim NH, Han KH, Choi J, Lee J, Kim SG. Use of Fenofibrate on Cardiovascular Outcomes in Statin Users With Metabolic Syndrome: Propensity Matched Cohort Study. BMJ (2019) 366:l5125. doi: 10.1136/bmj.l5125

PubMed Abstract | CrossRef Full Text | Google Scholar