- 1Department of Cardiology, the Affiliated Cardiovascular Hospital of Kunming Medical University (Fuwai Yunnan Cardiovascular Hospital), Kunming, China
- 2Department of Orthopedic and Trauma Surgery, the First Affiliated Hospital, Guangxi Medical University, Nanning, China
Purpose: To identify biomarkers associated with CD8+ T cells in coronary artery disease (CAD) and initially explore their potential role in the tumor immune microenvironment.
Materials and Methods: CAD-related datasets GSE12288, GSE34198, and GSE66360, were downloaded from the GEO database. First, GSVA was performed based on the GSE12288 dataset. Then WGCNA analysis was performed to identify the most relevant module and candidate hub gene for CD8+ T cells, followed by GO and KEGG analysis of this module. Secondly, the relationship between candidate hub genes and CD8+ T cells was verified using GSE34198 and GSE66360, which led to the identification of hub genes. The relationship of hub genes with CD8+ T cells in cancer was analyzed using the TIMER database. Methylation analysis of hub genes was performed using the DiseaseMeth database. CAD, pan-cancer, pan-cell lines, and pan-normal tissues, correlations between hub genes. In addition, potential drugs and TFs associated with hub genes were predicted, and the ceRNA network was constructed. Finally, GSEA was performed separately for hub genes.
Results: CAD was shown to be associated with immune response by GSVA analysis. WGCNA identified the blue module as most related to CD8+ T cells and identified nine candidate hub genes. The relevance of CAD to immunity was further confirmed by GO and KEGG analysis of the module. Two additional datasets validated and identified three hub genes (FBXO7, RAD23A, and MKRN1) that significantly correlated with CD8+ T cells. In addition, we found that hub genes were positively associated with CD8+ T cells in TGCT, THCA, and KICH cancers by our analysis. Moreover, the hub gene was differentially methylated. We also analyzed the correlation between hub genes in CAD, different cancers, different cell lines, and different normal tissues. The results of all the analyses showed a positive correlation between them. Finally, we successfully constructed hub gene-associated TF-gene and ceRNA networks and predicted 11 drugs associated with hub genes. GSEA suggests that hub genes are related to multiple immune response processes.
Conclusion: FBXO7, RAD23A, and MKRN1 are significantly associated with CD8+ T cells in CAD and multiple cancers and may act through immune responses in CAD and cancer.
Introduction
Coronary atherosclerotic heart disease, or coronary artery disease (CAD), is the most common form of cardiovascular disease. CAD is the most common cause of human death, both domestically and abroad (1). To date, the gold standard for the diagnosis of CAD is the performance of coronary angiography, which has the highest accuracy. Still, this means of confirming the diagnosis is invasive and can only be done by cardiologists who have specialized studies in the cardiac intervention (2). Moreover, pharmacological treatment is essential throughout, in addition to medical intervention and surgical bypass surgery (3). In addition, cancer has a very high morbidity and mortality rate, and recent studies have found that cancer patients often have a combination of heart disease, especially CAD (4). However, little has been reported about the relationship between CAD and cancer. Therefore, the exploration of potential biomarkers associated with CAD and the exploration of possible mechanisms and common therapeutic targets between CAD and cancer deserve further investigation.
Atherosclerosis is the most basic pathological mechanism of CAD (5, 6). The immune response plays an integral role in the formation as well as in the development of atherosclerosis (7, 8). Among them, T lymphocytes are the primary immune cell type in atherosclerotic plaques, with a predominance of cytotoxic CD8+ T cells (9, 10). Animal experiments have shown that cytotoxic CD8+ T cells promote the formation of necrotic cores and exacerbate the progression of unfavorable plaques (11). It has been found that antibody-mediated depletion of CD8+ T cells ameliorates atherosclerosis (11). Also, CD8+ T cells can induce apoptosis and inflammatory responses by secreting cytokines such as tumor necrosis factor-α and IFN-γ (12, 13). In addition, CD8+ T cell function is tightly controlled by immune checkpoint proteins. Regulators of immune checkpoint proteins, such as the E3-ligase Casitas B-cell lymphoma B (CBL-B), control the progression of human atheromatous plaques mainly by limiting the activation CD8+ T cells (14). In contrast, co-suppressors that inhibit immune checkpoints, such as cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), programmed cell death protein 1 (PD-1), and programmed cell death one ligand 1 (PD-L1), can increase cytotoxic CD8+ T cells and thus aggravate atherosclerotic plaques (15). Tumor immunotherapy most notably enhances the activity of CD8+ T cells to obtain a sustained and effective anticancer response (16). In the tumor microenvironment (TME), CTLA-4 inhibits the initiation of initial T cells, while PD-1 and PD-L1 accelerate effector T cell depletion (17). Blocking CTLA-4 or PD-1 enhances CD8+ T cell-mediated immune responses (18). In addition, immunomodulatory antibodies or cancer vaccines can restart CD8+ T-cell activity (19). In conclusion, immunotherapy has been widely used in the clinical treatment of cancer but is limited to laboratory studies on CAD.
Bioinformatics analysis methods have been utilized to conduct relevant explorations to gain insight into the underlying pathophysiological mechanisms of diseases. Among them, weighted gene co-expression network analysis (WGCNA) is a widely used bioinformatics analysis tool that can mine the gene expression modules most relevant to the immune response to disease, from which hub genes can be further screened (20). On the other hand, CIBERSORT is another bioinformatics analysis method that can identify the type of immune cell composition (21). In recent years, the combined application of WGCNA and CIBERSORT to explore potential biomarkers has been commonly used in oncology. However, few studies have used it to investigate potential biomarkers associated with CAD immune responses, especially to explore biomarkers related to CD8+ T cells.
We used the Gene Expression Omnibus (GEO) database to download the gene expression dataset GSE12288 from patients with CAD and first analyzed the enrichment of all sample genes in this dataset using GSVA. This dataset was then analyzed using WGCNA to identify the modules most associated with CD8+ T cell infiltration levels and candidate hub genes. Then, GO, and KEGG analyses were performed on this module to confirm the association of CAD with CD8+ T cells. In addition, we further explored the relevance of these candidate hub genes to CD8+ T cells using other datasets, GSE34198 and GSE66360, thus identifying three hub genes in CAD with CD8+ T cells. Finally, we also explored their role in different cancer types. This is the first time we have used WGCNA combined with immune infiltration analysis to identify potential molecular markers associated with CD8+ T cells in CAD while further exploring their possible immune mechanisms in CAD and different cancers. This discovery may provide a valid theoretical basis for the future development of novel immunotherapies for CAD and cancer and provide new insights into the possible link between CAD and cancer.
Materials and Methods
Gene Set Variation Analysis
Gene Set Variation Analysis (GSVA) is a non-parametric, unsupervised analysis method (22). It evaluated the enrichment of different signaling pathways of genes in the GSE12288 dataset in different samples. Gene sets C2 (containing GO information) and C5 (containing KEGG information) were first downloaded using the msigdbr package. Finally, gene set variation analysis was performed using the limma and GSVA packages.
Acquisition and Processing of Gene Expression Data
The GSE12288 dataset was retrieved and downloaded from Gene Expression Omnibus (GEO, http://www.ncbi.), containing 110 coronary heart disease samples and 112 standard human samples (23). The dataset was first normalized using the normalized between arrays function of the limma package. Then the conversion between probe IDs and gene symbols was performed according to the GPL96 platform information. Probes without gene symbols were removed, and for multiple probes corresponding to the same gene symbol, the average expression value of each search was taken. We use the coefficient of variation values to select the genes with the most significant variation, select genes with gene variation coefficients greater than 0.1, and then use these genes to construct the network.
Evaluation of Infiltration of CAD Immune Cells
The CIBERSORT algorithm is often used to calculate the proportion of 22 immune cell types (21). To understand the infiltration of immune cells in the samples in the dataset GSE12288, we used the R package “CIBERSORT” to calculate the raw data of GSE12288.
Construction of Co-Expression Network
Based on the expression of genes with coefficients of variation greater than 0.1 in the GSE12288 dataset, we used the R package “WGCNA” to construct a gene co-expression network (20). The steps are as follows: we first convert the expression levels of individual transcripts into a similarity matrix based on the Pearson correlation values after gene pairing. Then, the similarity matrix is transformed into an adjacency matrix according to the formula as amn=|cMn|β (cMn=Pearson correlation between paired genes; amn=adjacency relationship between paired genes). When the power of β = 3, the adjacency matrix can be transformed into a topological overlap matrix. Finally, we choose the module minimum cutoff value of 30 and use the bottom-up algorithm to divide the genes with similar expression patterns into different modules.
Building Module Feature Relationships
Component analysis was performed for each module according to its characteristics. We used Pearson’s test to analyze the correlation between module features and T-cell subtypes. When p<0.05, we considered the module to have a significant correlation with T-cell subtypes. We selected the module with the highest correlation coefficient with CD8+ T cells as the hub module.
Functional and Pathway Enrichment Analysis
To further elucidate the biological functions and signaling pathways of the identified genes in the hub module, we performed gene ontology (GO) analysis (24) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis (25) using the “cluster profile” R package (26). Among them, GO analysis includes molecular functions (MF), cellular components (CC), and biological processes (BP).
Identification and Validation of Hub Genes
We selected candidate hub genes based on module connectivity and clinical trait relationships for each gene in the hub module. Module connectivity was defined as the absolute value of Pearson correlation between genes (module membership). Clinical trait relationship is defined as the total value of Pearson correlation between each gene and trait (gene significance). For candidate hub genes, we set Module Membership>0.8 and gene significance>0.5. To verify whether there exists a significant positive correlation between candidate hub genes and CD8+ T, and thus identify reliable hub genes. We downloaded the GSE34198 (27) and GSE66360 (28) datasets and calculated the CD8+ T cells in each specimen using the cybersport algorithm. Spearman correlation between CD8+ T cell infiltration level and hub gene expression was also calculated. We defined hub genes as candidate hub genes significantly correlated with CD8+ T cells in the GSE12288, GSE34198, and GSE66360 datasets. In addition, we further validated the Spearman correlation between hub gene expression levels and CD8+ T cell infiltration levels using the Tumor Immune Evaluation Resource (TIMER) database. We visualized the results with “ggplots2”.
Hub Gene Methylation Levels
Methylation gives each CD8+ T cell subpopulation its unique epigenetic fingerprint. Thus these cells can rapidly produce corresponding effects when they reencounter cognate antigens. DiseaseMeth 2.0 database stores abnormal methylation data for various human diseases (29). Therefore, we used it to explore the relationship between FBXO7, RAD23A, and MKRN1, three hub genes expression, and methylation levels in CAD.
Hub Gene Correlation
In addition to the GSE12288 dataset, we also used two other datasets, GSE34198 and GSE66360, to synthesize the correlation between hub genes and then used “ggplots2” to visualize the analysis results. To further demonstrate the correlation between hub genes, we also downloaded high-throughput expression data of cancer, cell lines, and normal tissues from three databases, including TCGA, CCLE, and GETx, respectively. By calculating the correlation between hub genes in different cancers, different cell lines, and different normal tissues. Finally, “ggplots2” was used for visualization.
The Linkage Between Hub Genes
mRNA controls many critical biological processes, and the functions it performs are closely related to its localization in the cell. mRNA located database is a comprehensive and advanced model for accurate eukaryotic mRNA subcellular localization expression. We used it to analyze the location of these three hub genes in the cell. Transcription factors regulate transcriptional processes by recognizing specific DNA sequences, thus participating in various complex biological processes. The ChEA3 database integrates RNA-seq data from ENCODE, ReMap, GTEx, TCGA, and ARCHS4, using which common transcription factors for multiple genes can be predicted. Since the mechanism of the role of ceRNA in CAD is still unclear, we used the ENCORE database to predict the microRNA of HUB genes in reverse and the LncRNA of the common microRNA of HUB genes, and finally constructed the ceRNA network. In addition, we completed the GO annotation data required for three different aspects of GOSemSim from BP, CC, and MF levels of the corresponding species. The functional similarity among hub genes was searched for by comprehensively assessing the gene’s function in the three groups.
Small Molecule Drug
The CMAP (Connectivity Map) database has the most comprehensive data on the transcriptome of drug interference therapies (30). We placed FBXO7, RAD23A, and MKRN1 into CMAP, defining a negative connectivity score as a potential therapeutic agent. Therefore, we used enrichment <-0.7 and P<0.0 1 for screening.
Gene Set Enrichment Analysis
Gene set enrichment analysis (GSEA) can be used to clarify whether a set of gene sets is significantly different between two biological states (31). To further explore the potential mechanisms by which the three identified hub genes affect CAD, we performed a single-gene GSEA analysis for this purpose. In the GSE12288 dataset, according to the expression levels of hub genes, their correlation coefficients with other genes were calculated separately using spearman analysis, and the correlation coefficients were sorted as the input gene list of the clusterProfilerR package. The h.all.v7.4.sytmbols.gmt from the Molecular Signature Database (MSigDB) was selected as the reference gene set, and a P-regulation value <0.05 was used as the screening criterion.
Results
Gene Set Variation Analysis of Gene Expression Data
To investigate the potential role of genes in the GSE12288 dataset associated with CAD, we performed a GSVA analysis (Figures 1, 2). These genes were significantly enriched in some immune-related pathways, such as il6 jak stat3 signaling, interferon-gamma response, and il2 stat5 signaling.
Figure 1 Biological processes and pathways associated with all sample genes in the GSE12288 dataset related to CAD. Hallmark section of GSVA.
Gene Expression Data for CAD
We obtained microarray expression data for coronary artery disease from the GEO database. We selected 664 genes with coefficients of variation greater than 0.1 for the subsequent analysis step.
Evaluation of CAD-Infiltrating Immune Cells
We analyzed the corresponding expression data in the dataset GSE12288 using the R program “CIBERSORT” to determine the proportion of different immune cell isoforms in each sample in the dataset. Afterward, the proportions of seven T-cell subtypes in each piece were selected as trait data for WGCNA.
Gene Co-Expression Network of CAD
Based on the expression levels of 664 genes in the dataset GSE12288, we constructed a gene co-expression network using the R package “WGCNA.” Then, we calculated the average linkage and Pearson correlation coefficients and performed cluster analysis on all samples in this dataset. We successfully constructed a scale-first network with β=3 as the soft threshold power (Figure 3). We also used a dynamic hybrid cut method to build a hierarchical clustering tree. Each leaf represents a gene, and each branch represents a module, which aggregates all genes with similar expression levels. We then summed the functionally equivalent modules into one large module and obtained four modules (Figure 3).
Figure 3 WGCNA. (A) Scale-free fit index analysis of 1-20 soft threshold power (left) as well as mean connectivity analysis (right); (B) Gene classification into different modules by hierarchical clustering, with different colors representing different modules; (C) Heat map showing module features associated with T-cell infiltration; (D) Scatter plot of genes in blue modules.
Hub Module Identification and Functional Enrichment Analysis
The blue module of the four modules had the most significant relationship with T-cell CD8 (CD8+ T cells) (R2 = 0.27, P=4E-05; Figure 3). Therefore, we defined the blue module as the hub module. In this module, the genes with higher linkage were identified as potentially essential genes most associated with CD8+ T cells. Based on the truncation criteria (module affiliation > 0.8 and gene significance > 0.5), nine genes were selected (PAX8, ATP6V0C, MMP14, PSMF1, RAD23A, PRDX6, MKRN1, YBX3, FBXO7) as candidate hub genes (Figure 3). Next, the genes in the blue module were analyzed by GO and KEGG analysis using the “cluster profile” R package. In BP, module genes were significantly enriched in regulating protein catabolic process, cellular amino acid metabolic process, and positive regulation of intracellular transport (Figure 4). In CC, modular genes were significantly enriched in proteasome complex, endopeptidase complex, and melanosome (Figure 4). In MF, modular genes were increased dramatically in ubiquitin-protein ligase binding, ubiquitin-like protein ligase binding, and proteasome binding (Figure 4). In KEGG, modular genes were significantly enriched in Huntington disease, Proteasome, and Parkinson’s disease (Figure 4).
Figure 4 GO and KEGG. To clarify the functions and signaling pathways of genes in the blue module, we performed not only GO analysis containing biological processes (A), cellular components (B), and molecular functions (C), but also KEGG pathway (D) analysis.
Identification and Validation of Hub Genes
We used three datasets, GSE12288, GSE34198, and GSE66360, to identify and validate the correlation between PAX8, ATP6V0C, MMP14, PSMF1, RAD23A, PRDX6, MKRN1, YBX3, and CD8+ T cell infiltration levels, and finally determined that FBXO7, RAD23A, and MKRN1 could serve as reliable Hub genes. The analysis results showed that all three hub genes, FBXO7, RAD23A, and MKRN1, were significantly and positively correlated with the degree of infiltration of CD8+ T cells in these three datasets (Figures 5A). To further explore the expression values of FBXO7, RAD23A, and MKRN1 genes concerning CD8+ T cells in different cancers, the correlations of FBXO7, RAD23A, and MKRN1 expression values with CD8+ T cells in the TIMER database were next obtained and downloaded and visualized using the “ggplot “ R package for visualization. It was found that the expression values of FBXO7 were also significantly and positively correlated with the degree of CD8+ T cell infiltration in several cancers (Figure 6), followed by RAD23A (Figure 6) and MKRN1 (Figure 6). These analyses confirm that the identified HUB genes are closely associated with the level of CD8+ T-cell infiltration and play an essential role in the tumor immune microenvironment.
Figure 5 Scatter plots of hub gene expression versus CD8+ T cell infiltration levels in three different datasets: (A-C) GSE12288. (D-F) GSE34198. (G-I) GSE66360.
Figure 6 The relationship between the three hub gene expression levels and the degree of CD8+ T cell infiltration in different cancers was statistically significant at P<0.05: (A) FBXO7. (B) RAD23A. (C) MKRN1.
Methylation Levels of Hub Genes
Data consistent with methylation levels of the 3 HUB genes in coronary artery disease were obtained from DiseaseMeth 2.0. analysis showed that FBXO7 showed significantly lower methylation levels (Figure 7), and RAD23A had substantially higher methylation levels (Figure 7), whereas the methylation results for MKRN1 were not statistically significant (Figure 7).
Figure 7 Methylation analysis of three hub genes (A) FBXO70, (B) RAD23A, and (C) MKRN1 in patients with coronary artery disease (left) versus normal subjects (right).
Hub Gene Correlation
The analysis showed a significant positive correlation between FBXO7, RAD23A, and MKRN1 in the three datasets GSE12288, GSE34198, and GSE66360 (Figures 8A). We next analyzed the correlation between FBXO7, RAD23A, and MKRN1 in different cancers. Except for USC, MESO, and CHOL, FBXO7 was significantly positively correlated with MKRN1 in the remaining cancers (Figure 9). MKRN1 was significantly positively correlated with RAD23A in the remaining cancers except for MESO, LGG, UCS, CHOL, READ, and BRCA (Figure 9). Except for UVM, CHOL, UCS, SKCM, MKRN1 was significantly positively correlated with RAD23A in the remaining cancers (Figure 9). The correlation between FBXO7 and MKRN1 in cell lines is shown in Figure 9, the correlation between FBXO7 and RAD23A in cell lines is shown in Figure 9, and the correlation between MKRN1 and RAD23A in cell lines is shown in Figure 9. FBXO7 was significantly positively correlated with MKRN1 in cardiac tissues (Figure 9), FBXO7 was significantly positively correlated with RAD23A in cardiac tissues (Figure 9), and MKRN1 was significantly positively correlated with RAD23A in cardiac tissues (Figure 9).
Figure 8 Relationship between the expression levels of three hub genes in three different gene sets (A-C) GSE12288. (D-F) GSE34198. (G-I) GSE66360.
Figure 9 Relationship between the expression levels of three hub genes in different cancer types (A-C), different cell lines (D-F), and normal tissues (G-I).
Linkage Between Hub Genes
The mRNA subcellular localization results showed that FBXO7 was mainly localized in Nucleus and cytoplasm (Figure 10), RAD23A was mainly localized in cytoplasm and Nucleus (Figure 10), and MKRN1 was mainly localized in Nucleus and cytoplasm (Figure 10). We identified 19 common transcription factors associated with FBXO7, RAD23A and MKRN1 from the CHEA3 database (Figure 10). Three common microRNAs associated with FBXO7, RAD23A and MKRN1 (hsa-miR-335-5p, hsa-miR-340-5p and hsa-miR-514a-3p) were identified from the ENCORI database, and common LncRNAs for common microRNAs (AC021078.1, NEAT1, LINC00943, TUG1) (Figure 10). The results of functional similarity analysis showed that there were functional similarities among FBXO7, RAD23A and MKRN1, with FBXO7 occupying the most important position (Figure 10).
Figure 10 Subcellular localization of 3 hub genes (A) FBXO70, (B) RAD23A, and (C) MKRN1. (D) A network composed of target gene transferrin (Gen-TF). The blue diamond is the gene, and the pink triangle is the TF; (E) CircRNA-miRNA-Hubgene network. The network consists of four CircRNAs, three miRNAs, and three Hub genes. The yellow V-shape represents mRNAs, the pink diamond represents miRNAs, and the purple circle represents differentially expressed CircRNAs; (F) Functional similarity analysis of the three hub genes.
Small Molecule Drug Analysis
We entered FBXO7, RAD23A, and MKRN1 into the camp database and eventually identified 11 drugs that could be used to treat coronary artery disease, including doxorubicin, gabapentin, suxibuzone, prasterone, kawain, 5255229, decline, nalbuphine, MK-801, trimethoprim, and GW-8510 (Table 1).
Gene Enrichment of GSEA
Based on “h.all.v7.4.entrez.gmt” as the reference genome, we found the complete list of gene sets enriched in samples related to FBXO7 (Figure 11), RAD23A (Figure 11), and MKRN1 (Figure 11). Based on “c2.cp.kegg.v7.4.entrez.gmt” as the reference genome, we found the complete list of gene sets enriched in samples related to FBXO7 (Figure 11), RAD23A (Figure 11), and MKRN1 (Figure 11). Then, we selected the gene sets associated with immunity and CAD from the complete list for further analysis. When we used “h.all.v7.4.entrez.gmt” as the reference genome, the analysis revealed that the samples of the highly correlated FBXO7 (Figure 12) were enriched for allograft rejection, interferon-gamma response, and tnfa signaling via NF-kB in 3 sets of gene sets; In contrast, three of the five gene sets enriched in samples from the highly correlated RAD23A (Figure 12) were the same as those from FBXO7, in addition to the inflammatory response and the interferon-alpha response; Six gene sets were enriched in the highly correlated MKRN1 (Figure 12) samples, except for five of them, the same as the RAD23A enrichment results, and another set of pi3k Akt mtor signaling. And when we used “c2.cp.kegg.v7.4.entrez.gmt” as the reference genome, three sets of gene sets were enriched in the samples of the highly associated FBXO7 (Figure 12), including dilated cardiomyopathy, hypertrophic cardiomyopathy hcm, and nod like receptor signaling pathway; Six gene sets were enriched in samples with highly correlated RAD23A (Figure 12), three of which were improved in the same way as FBXO7 and the other three in the B cell receptor signaling pathway, graft versus host disease, and arrhythmogenic right ventricular cardiomyopathy arvc; Four gene sets were enriched in the highly correlated MKRN1 (Figure 12) samples, except for three of them, the same as the FBXO7 enrichment results and one set of B cell receptor signaling pathways.
Figure 11 Gene set enrichment analysis (GSEA) was performed to analyze hallmark and KEGG. The entire gene set enriched in samples containing FBXO70 (A, D), RAD23A (B, E), and MKRN1 (C, F) was highly expressed.
Figure 12 Gene sets are associated with immunity or heart disease. The gene sets related to immunity or heart disease were highly correlated in samples enriched for FBXO70 (A, D), RAD23A (B, E), and MKRN1 (C, F).
Discussion
Despite the current high level of medical care available, CAD and cancer remain important causes of human mortality worldwide. A growing number of studies show that immune responses are essential in developing both CAD and cancer. It is well known that immunotherapy is an integral part of anti-cancer treatment. Several common types of immunotherapies, such as immune checkpoint blockade (ICB), chimeric antigen receptor (CAR)-T cell therapy, cytokine therapy, immunomodulatory antibodies, monoclonal antibodies, and cancer vaccines, have achieved better anti-cancer responses (32, 33). In contrast, no immunotherapies have been identified that can be applied to the clinical treatment of CAD. In recent years, findings exploring immunotherapeutic approaches for CAD, such as the Canakinumab Anti-Inflammatory Thrombosis Outcome Study (CANTOS) and the Cardiovascular Inflammation Reduction Trial (CIRT), have shown that only blocking specific signaling pathways associated with atherogenesis can be genuinely therapeutic (34, 35). This finding has further stimulated interest in exploring potential target biomarkers for immune correlates of CAD. In addition, patients with cancer combined with CAD may have accelerated CAD progression due to the pro-inflammatory response and hypercoagulability brought about by cancer itself (36). Moreover, anticancer therapy in cancer patients may also increase the risk of developing CAD due to its cardiotoxic effects (37). However, there is little evidence regarding the exact mechanism of interaction between CAD and cancer to date.
This study used WGCNA in the GSE12288 dataset to find the most relevant modules to CD8+ T cells for the following analysis. FBXO7, RAD23A, and MKRN1 were verified and identified as hub genes by using GSE34198 and GSE66360 datasets. Their expression levels were positively correlated with CD8+ T-cell infiltration in CAD and multiple cancers. Among them, FBXO7, as a member of the F-box protein family, can be involved in the pathological process of Parkinson’s disease (PD) development through various mechanisms such as indirect ubiquitination and disruption of autophagic function (38, 39). FBXO7 is a component of the SCF (SKP1/cullin-1/F-box protein)-E3 type ubiquitin ligase complex, which leads to the production of the corresponding substrate ubiquitination (40). Previous correlations have reported that FBXO7 is expressed at high levels in cancer and is a potential oncogene (41). It is associated with immune escape and immunotherapy resistance in cancer patients (42). RAD23A, also known as RAD23 and HHR23A, acts as a junction protein during ubiquitin-mediated proteasomal degradation, increasing the degradation of target proteins and enhancing the response of CD8+ T cells (43). Human tumor virus E6-associated protein (E6AP) acts as an E3 ligase in ubiquitination. HHR23A participates in the ubiquitination process by binding to E6AP (44). rAD23A increases cancer chemoresistance and the chance of recurrence by upregulating autophagy (45). Makorin Ringinger Protein 1 (MKRN1) is an E3-type ubiquitin ligase. It has been shown that MKRN1 promotes p21 protein ubiquitination and proteasome pathway degradation, thereby preventing intermittent hypoxia (IH)-induced myocardial apoptosis (46). MKRN1 can disrupt p53 and thus apoptosis through ubiquitination (47). P53 is known to regulate autophagy, according to previous studies (48). Of course, further experiments are needed to verify whether there is a link between MKRN1, p53, and autophagy. MKRN1 positively regulates biological processes mediated by Wnt/β-linked proteins and plays an essential role in the proliferation, migration, and invasion of cancer cells (49). In addition, it has been reported that MKRN1 may be involved in regulating the development of gastric carcinogenesis (50). depletion of MKRN1 may inhibit tumor growth (51). It has been previously found that ubiquitination is a physiological process by which ubiquitin binds to target proteins to cause a series of reactions in substrates. It plays a crucial role in regulating immune responses (52). Moreover, ubiquitination plays an essential role in the pathogenesis of CAD, and its expression level increases with the severity of CAD (52, 53). On the other hand, autophagy is a conserved evolutionary process that provides energy and biological components by digesting its dysfunctional organelles (54, 55). Early correlations reported that autophagy could be stimulated in surviving cardiomyocytes to produce autophagy for self-preservation in the context of long-term chronic myocardial ischemia (56, 57). In addition, ubiquitination and autophagy, two major pathways of intracellular protein degradation, can be presented to the surface of MCH class I molecules, thereby activating CD8+ T cells and thus inducing an immune response in the organism (58). The ubiquitination of three hub genes, all associated with E3 ligases. Targeting the E3 ligase inhibits tumor progression (59). In the early stages of cancer, autophagy inhibits carcinogenesis development, while it promotes tumor survival and metastasis in later stages (60). Thus, our study provides a potential direction for the involvement of FBXO7, RAD23A, and MKRN1 in the immune response to CAD through ubiquitination or autophagy. Meanwhile, they may play an oncogenic role in cancer and are closely related to CD8+ T cells.
In addition, we found that FBXO7, RAD23A, and MKRN1 were significantly positively correlated in cardiac tissues and most cancers. No correlation studies have been found on the role of these three core genes in the association between CAD and cancer. Our study provides a new direction for the joint research between CAD and cancer while possibly providing possible options for optimizing the treatment of potential CAD in cancer patients. Of course, these need further experimental validation and clinical studies. It has now been confirmed that the efficacy of immunotherapy such as ICB and CAR T-cell therapy for tumors is mainly influenced by epigenetic regulation. And epigenetic modifications associated with CD8+ T cells can enhance immunotherapy (61). As we all know, epigenetics primarily focuses on methylation modifications (62). Therefore, we performed methylation analysis of the three core genes and showed that except for MKRN1, which was not statistically significant, FBXO7 showed hypomethylation, while RAD23A was the exact opposite. This provides a new theory to study the methylation process of FBXO7 and RAD23A and improve potential CAD immune efficacy in cancer patients.
To further understand the linkage of hub genes, we identified 19 common transcription factors associated with FBXO7, RAD23A, and MKRN1 by constructing a TF-gene network. FOXO3 has been reported to be a lifespan-associated gene that increases mortality in elderly patients with CAD (63). At first, FOXO3 was mainly considered a tumor suppressor. However, recent studies have found that it is also involved in maintaining starter cells in leukemia and colon cancer metastasis (64, 65). It was found that blocking STAT3 activity may be an effective strategy for treating CAD (66). STAT3 is often over-activated in various cancers and is associated with tumor immunosuppression and poor prognosis (67). Targeting STAT3 is emerging as a potentially promising therapeutic modality for many cancers (68). These two transcription factors may serve as targets for CAD-cancer interactions. Another 17 transcription factors are not studied for their association with CAD. In addition, there is growing evidence that competitive endogenous RNAs (ceRNAs) play a key role in CAD and cancer (69). However, there is still a lack of ceRNA networks that can predict the diagnosis and treatment of CAD. In the present study, we identified three common microRNAs (hsa-miR-335-5p, hsa-miR-340-5p, and hsa-miR-514a-3p) associated with FBXO7, RAD23A, and MKRN1. It has been suggested that miR-335-5p may regulate the process of cardiomyocyte differentiation by activating the WNT and TGFβ signaling pathways (70). It can be used as a potential biomarker for various cancers, for example, osteosarcoma, breast cancer, and non-Hodgkin’s lymphoma (71–73). In lung cancer, reduced expression of has-miR-340-5p may mediate cisplatin resistance (74). miR-514a-3p promotes melanoma growth (75). All three microRNAs have been studied in cancer in a limited number of ways. However, none of them is associated with CAD so far. At the same time, we also reverse inferred four common LncRNAs (AC021078.1, NEAT1, LINC00943, TUG1) of the common microRNAs. Among them, it was found that NEAT1 may inhibit CAD cell apoptosis by activating the miR-140-3p/MAPK1 pathway (76). neat1 promotes the progression and metastasis of several cancers (77, 78). linc00943 is associated with gastric cancer and clear cell renal cell carcinoma (79, 80). tug1 regulates endothelial cell proliferation and migration through the Wnt pathway to stimulate diabetic atherosclerosis (81). tUG1 is involved in tumorigenesis and tumor progression (82). The ceRNA regulatory network constructed in this study may play a role in CAD and cancer development and progression. The TF-gene and ceRNA network suggests that three hub genes may act as a bridge between CAD and cancer, which of course requires further experimental studies.
Meanwhile, to predict potentially effective therapeutic agents for CAD, we applied the CMAP database to identify 11 therapeutic agents that may reverse abnormally high expression of hub genes associated with CAD. Among them, doxorubicin, although improving prognosis and mortality in cancer patients, increased cardiotoxic effects (83). Gabapentin alleviates the painful reaction caused after coronary artery bypass grafting (CABG) in patients with CAD (84). Delsoline plays an essential role in cardiac remodeling and heart failure progression after myocardial infarction (85). For the 11 drugs identified in this study, no relevant studies have been found so far in cardiovascular disease, except for the above three drugs for which there are few relevant reports. Our findings may provide an essential basis for future exploration of effective drugs for CAD treatment.
In addition to the above studies, to further understand the potential biological pathways between CAD and cancer, we performed GSVA analysis on all samples in the dataset GSE12288. Genes in this dataset were significantly enriched in biological pathways such as il6 jak stat3 signaling, interferon-gamma response (IFN-γ), and il2 stat5 signaling. Single gene GSEA analysis of the screened three hub genes was also performed. The results showed that they were significantly enriched in multiple pathways such as Inflammatory response, interferon-alpha response, interferon-gamma response, TNF signaling via NF-kB, allograft rejection, and pi3k Akt mTOR signaling, but also in various cardiomyopathies such as dilated cardiomyopathy and hypertrophic cardiomyopathy. A literature review revealed that flavanones could inhibit the IL6/JAK/STAT3/SOCS3 signaling pathway and improve atherosclerosis (86). The IL6/JAK/STAT3 signaling pathway plays an essential role in the proliferation, survival, and invasion of tumor cells (87, 88). Meanwhile, targeting the IL-6/JAK/STAT3 signaling axis downregulates the expression levels of the CD8+ T cell surface receptors PD-1 or PD-L1 (89). Thus, IL-6/JAK/STAT3 may directly inhibit tumor cell growth or indirectly enhance the antitumor immune response effects of immune checkpoint inhibitors. CD8+ T cells may promote atherosclerosis by increasing the expression level of IFN-γ (11). CD8+ T cells can also induce cytotoxic effects in cancer cells by releasing IFN-γ (90). In the presence of IFN-γ, cancer cells can release several immunosuppressive mediators, including PD-L1 and transcriptional activator 3 (STAT3). These mediators negatively feedback to inhibit the release of IFN-γ, thereby suppressing the activity of CD8+ T cells (91, 92). IL-2 activates STAT5 in CD8+ T cells (93). H9T, an agonist of IL-2, leads to increased expression of CD8+ T cells. In animal models of melanoma and acute lymphoblastic leukemia, CAR-modified CD8+ T cells expanded with H9T have strong antitumor activity (94). The pro-inflammatory cytokine IFN-α can lead to atherosclerosis flare-up and progression (95). the antitumor effect of IFN-α is mediated mainly by CD8+ T cells (96). the inflammatory response has been considered a cancer marker and is closely associated with cancer progression (97). In addition, it is an essential pathogenic mechanism of atherosclerosis (5). MicroRNA-499 was found to cause endothelial cell injury by activating tnfa signaling via NF-kB (98). And NF-kB is also involved in tumorigenesis and development (99). Triclosan (TCS) can damage endothelial cells by inhibiting the PI3K/Akt/mTOR signaling pathway (100). The PI3K/Akt/mTOR pathway is usually activated in cancer, and inhibition may have an anti-cancer effect (101). These findings reaffirm that the development of CAD is associated with multiple immune processes mediated by CD8+ T cells. Moreover, these immune processes play a crucial role in cancer.
In conclusion, we identified FBXO7, RAD23A, and MKRN1 as biomarkers in CAD associated with immune responses mediated by CD8+ T cells. This finding helps to understand the part of the immune response in the development of CAD. At the same time, they play the role of potential therapeutic targets in various cancers. FBXO7, RAD23A and MKRN1 may act as hub genes linking CAD and cancer in terms of immunity. It may provide some theoretical basis for identifying new drug targets and developing new therapeutic approaches to reduce the incidence of CAD in cancer patients. Of course, there are some limitations to this study. On the one hand, FBXO7, RAD23A, and MKRN1 have not been reported in CAD. The mechanism of their roles in CAD needs to be studied in depth. On the other hand, the functions of these three hub genes in CAD and cancer are only based on our bioinformatics analysis, and future in vitro or in vivo experiments are needed to validate these findings further.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary materials.
Author Contributions
JZ, YW, XX, YanW, and SZ substantially contributed to conception and design. SZ and YW conducted the literature search. SZ and YW contributed to the acquisition of data or analysis and interpretation of data. SZ wrote the article. YW performed data analysis and drafted. JZ revised the article. All authors gave the final approval of the version to be submitted.
Funding
This work was supported by scientific research funding from the 2018 Science and Technology Plan Project of the Yunnan Provincial Science and Technology Department. Grant Program [2018FE001 (–285)].
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 wish to thank all the patients and the staff members of the medical staff for their cooperation.
References
1. Wang H, Naghavi M, Allen C, Barber R M, Bhutta Z A, Carter A, et al Global, Regional, and National Life Expectancy, All-Cause Mortality, and Cause-Specific Mortality for 249 Causes of Death, 1980-2015: A Systematic Analysis for the Global Burden of Disease Study 2015. Lancet (2016) 388(10053):1459–544. doi: 10.1016/S0140-6736(16)31012-1
2. Özcan C, Deleskog A, Schjerning Olsen AM, Nordahl Christensen H, Lock Hansen M, Hilmar Gislason G. Coronary Artery Disease Severity and Long-Term Cardiovascular Risk in Patients With Myocardial Infarction: A Danish Nationwide Register-Based Cohort Study. Eur Heart J Cardiovasc Pharmacother (2018) 4(1):25–35. doi: 10.1093/ehjcvp/pvx009
3. Virani SS, Alonso A, Benjamin EJ, Bittencourt MS, Callaway CW, Carson AP, et al. Heart Disease and Stroke Statistics-2020 Update: A Report From the American Heart Association. Circulation (2020) 141(9):e139–596. doi: 10.1161/CIR.0000000000000757
4. Koene RJ, Prizment AE, Blaes A, Konety SH. Shared Risk Factors in Cardiovascular Disease and Cancer. Circulation (2016) 133(11):1104–14. doi: 10.1161/CIRCULATIONAHA.115.020406
5. Weber C, Noels H. Atherosclerosis: Current Pathogenesis and Therapeutic Options. Nat Med (2011) 17(11):1410–22. doi: 10.1038/nm.2538
6. Gisterå A, Hansson GK. The Immunology of Atherosclerosis. Nat Rev Nephrol (2017) 13(6):368–80. doi: 10.1038/nrneph.2017.51
7. Hansson GK, Hermansson A. The Immune System in Atherosclerosis. Nat Immunol (2011) 12(3):204–12. doi: 10.1038/ni.2001
8. Libby P, Lichtman AH, Hansson GK. Immune Effector Mechanisms Implicated in Atherosclerosis: From Mice to Humans. Immunity (2013) 38(6):1092–104. doi: 10.1016/j.immuni.2013.06.009
9. Bergström I, Backteman K, Lundberg A, Ernerudh J, Jonasson L. Persistent Accumulation of Interferon-γ-Producing CD8+CD56+ T Cells in Blood From Patients With Coronary Artery Disease. Atherosclerosis (2012) 224(2):515–20. doi: 10.1016/j.atherosclerosis.2012.07.033
10. Fernandez DM, Rahman AH, Fernandez NF, Chudnovskiy A, Amir ED, Amadori L, et al. Single-Cell Immune Landscape of Human Atherosclerotic Plaques. Nat Med (2019) 25(10):1576–88. doi: 10.1038/s41591-019-0590-4
11. Kyaw T, Winship A, Tay C, Kanellakis P, Hosseini H, Cao A, et al. Cytotoxic and Proinflammatory CD8+ T Lymphocytes Promote Development of Vulnerable Atherosclerotic Plaques in apoE-Deficient Mice. Circulation (2013) 127(9):1028–39. doi: 10.1161/CIRCULATIONAHA.112.001347
12. Nakajima C, Uekusa Y, Iwasaki M, Yamaguchi N, Mukai T, Gao P, et al. A Role of Interferon-Gamma (IFN-Gamma) in Tumor Immunity: T Cells With the Capacity to Reject Tumor Cells are Generated But Fail to Migrate to Tumor Sites in IFN-Gamma-Deficient Mice. Cancer Res (2001) 61(8):3399–405.
13. van Horssen R, Ten Hagen TL, Eggermont AM. TNF-Alpha in Cancer Treatment: Molecular Insights, Antitumor Effects, and Clinical Utility. Oncologist (2006) 11(4):397–408. doi: 10.1634/theoncologist.11-4-397
14. Seijkens TTP, Poels K, Meiler S, van Tiel CM, Kusters PJH, Reiche M, et al. Deficiency of the T Cell Regulator Casitas B-Cell Lymphoma-B Aggravates Atherosclerosis by Inducing CD8+ T Cell-Mediated Macrophage Death. Eur Heart J (2019) 40(4):372–82. doi: 10.1093/eurheartj/ehy714
15. Poels K, van Leent MMT, Boutros C, Tissot H, Roy S, Meerwaldt AE, et al. Immune Checkpoint Inhibitor Therapy Aggravates T Cell-Driven Plaque Inflammation in Atherosclerosis. JACC CardioOncol (2020) 2(4):599–610. doi: 10.1016/j.jaccao.2020.08.007
16. Schillaci R, Salatino M, Cassataro J, Proietti CJ, Giambartolomei GH, Rivas MA, et al. Immunization With Murine Breast Cancer Cells Treated With Antisense Oligodeoxynucleotides to Type I Insulin-Like Growth Factor Receptor Induced an Antitumoral Effect Mediated by a CD8+ Response Involving Fas/Fas Ligand Cytotoxic Pathway. J Immunol (2006) 176(6):3426–37. doi: 10.4049/jimmunol.176.6.3426
17. Du X, Tang F, Liu M, Su J, Zhang Y, Wu W, et al. A Reappraisal of CTLA-4 Checkpoint Blockade in Cancer Immunotherapy. Cell Res (2018) 28(4):416–32. doi: 10.1038/s41422-018-0011-0
18. Topalian SL, Taube JM, Anders RA, Pardoll DM. Mechanism-Driven Biomarkers to Guide Immune Checkpoint Blockade in Cancer Therapy. Nat Rev Cancer (2016) 16(5):275–87. doi: 10.1038/nrc.2016.36
19. Borst J, Ahrends T, Bąbała N, Melief CJM, Kastenmüller W. CD4(+) T Cell Help in Cancer Immunology and Immunotherapy. Nat Rev Immunol (2018) 18(10):635–47. doi: 10.1038/s41577-018-0044-0
20. Langfelder P, Horvath S. WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559
21. 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
22. Hänzelmann S, Castelo R, Guinney J. GSVA: Gene Set Variation Analysis for Microarray and RNA-Seq Data. BMC Bioinf (2013) 14:7. doi: 10.1186/1471-2105-14-7
23. Sinnaeve PR, Donahue MP, Grass P, Seo D, Vonderscher J, Chibout SD, et al. Gene Expression Patterns in Peripheral Blood Correlate With the Extent of Coronary Artery Disease. PloS One (2009) 4(9):e7037. doi: 10.1371/journal.pone.0007037
24. Harris M A, Clark J I, Ireland A, Lomax J, Ashburner M, Collins R, et al The Gene Ontology (GO) Project in 2006. Nucleic Acids Res (2006) 34(Database issue):D322–6. doi: 10.1093/nar/gkj021
25. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res (2000) 28(1):27–30. doi: 10.1093/nar/28.1.27
26. 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
27. Zbuk K, Xie C, Young R, Heydarpour M, Pare G, Davis AD, et al. BRCA2 Variants and Cardiovascular Disease in a Multi-Ethnic Study. BMC Med Genet (2012) 13:56. doi: 10.1186/1471-2350-13-56
28. Muse ED, Kramer ER, Wang H, Barrett P, Parviz F, Novotny MA, et al. A Whole Blood Molecular Signature for Acute Myocardial Infarction. Sci Rep (2017) 7(1):12268. doi: 10.1038/s41598-017-12166-0
29. Xiong Y, Wei Y, Gu Y, Zhang S, Lyu J, Zhang B, et al. DiseaseMeth Version 2.0: A Major Expansion and Update of the Human Disease Methylation Database. Nucleic Acids Res (2017) 45(D1):D888–d95. doi: 10.1093/nar/gkw1123
30. Subramanian A, Narayan R, Corsello SM, Peck DD, Natoli TE, Lu X, et al. A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles. Cell (2017) 171(6):1437–52.e17. doi: 10.1016/j.cell.2017.10.049
31. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc Natl Acad Sci U.S.A. (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
32. Yang Y. Cancer Immunotherapy: Harnessing the Immune System to Battle Cancer. J Clin Invest (2015) 125(9):3335–7. doi: 10.1172/JCI83871
33. Lohmueller J, Finn OJ. Current Modalities in Cancer Immunotherapy: Immunomodulatory Antibodies, CARs and Vaccines. Pharmacol Ther (2017) 178:31–47. doi: 10.1016/j.pharmthera.2017.03.008
34. Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory Therapy With Canakinumab for Atherosclerotic Disease. N Engl J Med (2017) 377(12):1119–31. doi: 10.1056/NEJMoa1707914
35. Ridker PM, Everett BM, Pradhan A, MacFadyen JG, Solomon DH, Zaharris E, et al. Low-Dose Methotrexate for the Prevention of Atherosclerotic Events. N Engl J Med (2019) 380(8):752–62. doi: 10.1056/NEJMoa1809798
36. Lee JM, Yoon CH. Acute Coronary Stent Thrombosis in Cancer Patients: A Case Series Report. Korean Circ J (2012) 42(7):487–91. doi: 10.4070/kcj.2012.42.7.487
37. Hess CN, Roe MT. Treatment of Coronary Artery Disease in Cancer Survivors: An Emerging Challenge. Coron Artery Dis (2017) 28(1):1–2. doi: 10.1097/MCA.0000000000000444
38. Randle SJ, Laman H. Structure and Function of Fbxo7/PARK15 in Parkinson's Disease. Curr Protein Pept Sci (2017) 18(7):715–24. doi: 10.2174/1389203717666160311121433
39. Lu J, Wu M, Yue Z. Autophagy and Parkinson's Disease. Adv Exp Med Biol (2020) 1207:21–51. doi: 10.1007/978-981-15-4272-5_2
40. Zhou ZD, Lee JCT, Tan EK. Pathophysiological Mechanisms Linking F-Box Only Protein 7 (FBXO7) and Parkinson's Disease (PD). Mutat Res Rev Mutat Res (2018) 778:72–8. doi: 10.1016/j.mrrev.2018.10.001
41. Laman H. Fbxo7 Gets Proactive With Cyclin D/Cdk6. Cell Cycle (2006) 5(3):279–82. doi: 10.4161/cc.5.3.2403
42. Shen JZ, Qiu Z, Wu Q, Zhang G, Harris R, Sun D, et al. A FBXO7/EYA2-SCF(FBXW7) Axis Promotes AXL-Mediated Maintenance of Mesenchymal and Immune Evasion Phenotypes of Cancer Cells. Mol Cell (2022) 82(6):1123–39.e8. doi: 10.1016/j.molcel.2022.01.022
43. Andersson HA, Passeri MF, Barry MA. Rad23 as a Reciprocal Agent for Stimulating or Repressing Immune Responses. Hum Gene Ther (2005) 16(5):634–41. doi: 10.1089/hum.2005.16.634
44. Kumar S, Talis AL, Howley PM. Identification of HHR23A as a Substrate for E6-Associated Protein-Mediated Ubiquitination. J Biol Chem (1999) 274(26):18785–92. doi: 10.1074/jbc.274.26.18785
45. Tan X, Wang HC, Liang RY, Chuang SM. Human Rad23A Plays a Regulatory Role in Autophagy. Biochem Biophys Res Commun (2016) 478(4):1772–9. doi: 10.1016/j.bbrc.2016.09.025
46. Bai X, Yang H, Pu J, Zhao Y, Jin Y, Yu Q. MKRN1 Ubiquitylates P21 to Protect Against Intermittent Hypoxia-Induced Myocardial Apoptosis. Oxid Med Cell Longev (2021) 2021:9360339. doi: 10.1155/2021/9360339
47. Lee EW, Lee MS, Camus S, Ghim J, Yang MR, Oh W, et al. Differential Regulation of P53 and P21 by MKRN1 E3 Ligase Controls Cell Cycle Arrest and Apoptosis. EMBO J (2009) 28(14):2100–13. doi: 10.1038/emboj.2009.164
48. Song YM, Lee WK, Lee YH, Kang ES, Cha BS, Lee BW. Metformin Restores Parkin-Mediated Mitophagy, Suppressed by Cytosolic P53. Int J Mol Sci (2016) 17(1):122. doi: 10.3390/ijms17010122
49. Lee HK, Lee EW, Seo J, Jeong M, Lee SH, Kim SY, et al. Ubiquitylation and Degradation of Adenomatous Polyposis Coli by MKRN1 Enhances Wnt/β-Catenin Signaling. Oncogene (2018) 37(31):4273–86. doi: 10.1038/s41388-018-0267-3
50. Ko A, Shin JY, Seo J, Lee KD, Lee EW, Lee MS, et al. Acceleration of Gastric Tumorigenesis Through MKRN1-Mediated Posttranslational Regulation of P14arf. J Natl Cancer Inst (2012) 104(21):1660–72. doi: 10.1093/jnci/djs424
51. Lee EW, Kim JH, Ahn YH, Seo J, Ko A, Jeong M, et al. Ubiquitination and Degradation of the FADD Adaptor Protein Regulate Death Receptor-Mediated Apoptosis and Necroptosis. Nat Commun (2012) 3:978. doi: 10.1038/ncomms1981
52. Herrmann J, Edwards WD, Holmes DR Jr., Shogren KL, Lerman LO, Ciechanover A, et al. Increased Ubiquitin Immunoreactivity in Unstable Atherosclerotic Plaques Associated With Acute Coronary Syndromes. J Am Coll Cardiol (2002) 40(11):1919–27. doi: 10.1016/S0735-1097(02)02564-0
53. Chen SM, Li YG, Wang DM, Zhang GH, Tan CJ. Expression of Heme Oxygenase-1, Hypoxia Inducible Factor-1alpha, and Ubiquitin in Peripheral Inflammatory Cells From Patients With Coronary Heart Disease. Clin Chem Lab Med (2009) 47(3):327–33. doi: 10.1515/CCLM.2009.073
54. Levine B, Kroemer G. Autophagy in the Pathogenesis of Disease. Cell (2008) 132(1):27–42. doi: 10.1016/j.cell.2007.12.018
55. Mizushima N, Levine B, Cuervo AM, Klionsky DJ. Autophagy Fights Disease Through Cellular Self-Digestion. Nature (2008) 451(7182):1069–75. doi: 10.1038/nature06639
56. Yan L, Vatner DE, Kim SJ, Ge H, Masurekar M, Massover WH, et al. Autophagy in Chronically Ischemic Myocardium. Proc Natl Acad Sci U.S.A. (2005) 102(39):13807–12. doi: 10.1073/pnas.0506843102
57. Kanamori H, Takemura G, Goto K, Maruyama R, Ono K, Nagao K, et al. Autophagy Limits Acute Myocardial Infarction Induced by Permanent Coronary Artery Occlusion. Am J Physiol Heart Circ Physiol (2011) 300(6):H2261–71. doi: 10.1152/ajpheart.01056.2010
58. Warnatsch A, Bergann T, Krüger E. Oxidation Matters: The Ubiquitin Proteasome System Connects Innate Immune Mechanisms With MHC Class I Antigen Presentation. Mol Immunol (2013) 55(2):106–9. doi: 10.1016/j.molimm.2012.10.007
59. Liu J, Peng Y, Zhang J, Long J, Liu J, Wei W. Targeting SCF E3 Ligases for Cancer Therapies. Adv Exp Med Biol (2020) 1217:123–46. doi: 10.1007/978-981-15-1025-0_9
60. Li X, He S, Ma B. Autophagy and Autophagy-Related Proteins in Cancer. Mol Cancer (2020) 19(1):12. doi: 10.1186/s12943-020-1138-4
61. Zebley CC, Gottschalk S, Youngblood B. Rewriting History: Epigenetic Reprogramming of CD8(+) T Cell Differentiation to Enhance Immunotherapy. Trends Immunol (2020) 41(8):665–75. doi: 10.1016/j.it.2020.06.008
62. Dawson MA, Kouzarides T. Cancer Epigenetics: From Mechanism to Therapy. Cell (2012) 150(1):12–27. doi: 10.1016/j.cell.2012.06.013
63. Willcox BJ, Morris BJ, Tranah GJ, Chen R, Masaki KH, He Q, et al. Longevity-Associated FOXO3 Genotype and its Impact on Coronary Artery Disease Mortality in Japanese, Whites, and Blacks: A Prospective Study of Three American Populations. J Gerontol A Biol Sci Med Sci (2017) 72(5):724–8. doi: 10.1093/gerona/glw196
64. Coomans de Brachène A, Demoulin JB. FOXO Transcription Factors in Cancer Development and Therapy. Cell Mol Life Sci (2016) 73(6):1159–72. doi: 10.1007/s00018-015-2112-y
65. Hornsveld M, Dansen TB, Derksen PW, Burgering BMT. Re-Evaluating the Role of FOXOs in Cancer. Semin Cancer Biol (2018) 50:90–100. doi: 10.1016/j.semcancer.2017.11.017
66. Plens-Galaska M, Szelag M, Collado A, Marques P, Vallejo S, Ramos-González M, et al. Genome-Wide Inhibition of Pro-Atherogenic Gene Expression by Multi-STAT Targeting Compounds as a Novel Treatment Strategy of CVDs. Front Immunol (2018) 9:2141. doi: 10.3389/fimmu.2018.02141
67. Kortylewski M, Kujawski M, Wang T, Wei S, Zhang S, Pilon-Thomas S, et al. Inhibiting Stat3 Signaling in the Hematopoietic System Elicits Multicomponent Antitumor Immunity. Nat Med (2005) 11(12):1314–21. doi: 10.1038/nm1325
68. Zou S, Tong Q, Liu B, Huang W, Tian Y, Fu X. Targeting STAT3 in Cancer Immunotherapy. Mol Cancer (2020) 19(1):145. doi: 10.1186/s12943-020-01258-7
69. Ji Y, Yan T, Zhu S, Wu R, Zhu M, Zhang Y, et al. The Integrative Analysis of Competitive Endogenous RNA Regulatory Networks in Coronary Artery Disease. Front Cardiovasc Med (2021) 8:647953. doi: 10.3389/fcvm.2021.647953
70. Kay M, Soltani BM, Aghdaei FH, Ansari H, Baharvand H. Hsa-miR-335 Regulates Cardiac Mesoderm and Progenitor Cell Differentiation. Stem Cell Res Ther (2019) 10(1):191. doi: 10.1186/s13287-019-1249-2
71. Chen Y, Chen Q, Zou J, Zhang Y, Bi Z. Construction and Analysis of a ceRNA−ceRNA Network Reveals Two Potential Prognostic Modules Regulated by Hsa−Mir−335−5p in Osteosarcoma. Int J Mol Med (2018) 42(3):1237–46. doi: 10.3892/ijmm.2018.3709
72. Fernandes M, Marques H, Teixeira AL, Medeiros R. ceRNA Network of lncRNA/miRNA as Circulating Prognostic Biomarkers in Non-Hodgkin Lymphomas: Bioinformatic Analysis and Assessment of Their Prognostic Value in an NHL Cohort. Int J Mol Sci (2021) 23(1):201. doi: 10.3390/ijms23010201
73. Mahmoudian M, Razmara E, Mahmud Hussen B, Simiyari M, Lotfizadeh N, Motaghed H, et al. Identification of a six-microRNA Signature as a Potential Diagnostic Biomarker in Breast Cancer Tissues. J Clin Lab Anal (2021) 35(11):e24010. doi: 10.1002/jcla.24010
74. Cui F, Hao ZX, Li J, Zhang YL, Li XK, He JX. SOX2 Mediates Cisplatin Resistance in Small-Cell Lung Cancer With Downregulated Expression of hsa-miR-340-5p. Mol Genet Genomic Med (2020) 8(5):e1195. doi: 10.1002/mgg3.1195
75. Stark MS, Bonazzi VF, Boyle GM, Palmer JM, Symmons J, Lanagan CM, et al. miR-514a Regulates the Tumour Suppressor NF1 and Modulates BRAFi Sensitivity in Melanoma. Oncotarget (2015) 6(19):17753–63. doi: 10.18632/oncotarget.3924
76. Zhang H, Ji N, Gong X, Ni S, Wang Y. NEAT1/miR-140-3p/MAPK1 Mediates the Viability and Survival of Coronary Endothelial Cells and Affects Coronary Atherosclerotic Heart Disease. Acta Biochim Biophys Sin (Shanghai) (2020) 52(9):967–74. doi: 10.1093/abbs/gmaa087
77. Zhang M, Weng W, Zhang Q, Wu Y, Ni S, Tan C, et al. The lncRNA NEAT1 Activates Wnt/β-Catenin Signaling and Promotes Colorectal Cancer Progression via Interacting With DDX5. J Hematol Oncol (2018) 11(1):113. doi: 10.1186/s13045-018-0656-7
78. Park MK, Zhang L, Min KW, Cho JH, Yeh CC, Moon H, et al. NEAT1 is Essential for Metabolic Changes That Promote Breast Cancer Growth and Metastasis. Cell Metab (2021) 33(12):2380–97.e9. doi: 10.1016/j.cmet.2021.11.011
79. Xu Y, Ji T, An N, Wang X, Zhang H, Xu F. LINC00943 is Correlated With Gastric Cancer and Regulates Cancer Cell Proliferation and Chemosensitivity via hsa-miR-101-3p. Int J Clin Oncol (2021) 26(9):1650–60. doi: 10.1007/s10147-021-01945-5
80. Zhang Y, Dai J, Huang W, Chen Q, Chen W, He Q, et al. Identification of a Competing Endogenous RNA Network Related to Immune Signature in Clear Cell Renal Cell Carcinoma. Aging (Albany NY) (2021) 13(24):25980–6002. doi: 10.18632/aging.203784
81. Yan HY, Bu SZ, Zhou WB, Mai YF. TUG1 Promotes Diabetic Atherosclerosis by Regulating Proliferation of Endothelial Cells via Wnt Pathway. Eur Rev Med Pharmacol Sci (2018) 22(20):6922–9. doi: 10.26355/eurrev_201810_16162
82. Zhou H, Sun L, Wan F. Molecular Mechanisms of TUG1 in the Proliferation, Apoptosis, Migration and Invasion of Cancer Cells. Oncol Lett (2019) 18(5):4393–402. doi: 10.3892/ol.2019.10848
83. Hershman DL, McBride RB, Eisenberger A, Tsai WY, Grann VR, Jacobson JS. Doxorubicin, Cardiac Risk Factors, and Cardiac Toxicity in Elderly Patients With Diffuse B-Cell non-Hodgkin's Lymphoma. J Clin Oncol (2008) 26(19):3159–65. doi: 10.1200/JCO.2007.14.1242
84. Ucak A, Onan B, Sen H, Selcuk I, Turan A, Yilmaz AT. The Effects of Gabapentin on Acute and Chronic Postoperative Pain After Coronary Artery Bypass Graft Surgery. J Cardiothorac Vasc Anesth (2011) 25(5):824–9. doi: 10.1053/j.jvca.2010.11.017
85. Li GH, Shi Y, Chen Y, Sun M, Sader S, Maekawa Y, et al. Gelsolin Regulates Cardiac Remodeling After Myocardial Infarction Through DNase I-Mediated Apoptosis. Circ Res (2009) 104(7):896–904. doi: 10.1161/CIRCRESAHA.108.172882
86. Wiejak J, Dunlop J, Mackay SP, Yarwood SJ. Flavanoids Induce Expression of the Suppressor of Cytokine Signalling 3 (SOCS3) Gene and Suppress IL-6-Activated Signal Transducer and Activator of Transcription 3 (STAT3) Activation in Vascular Endothelial Cells. Biochem J (2013) 454(2):283–93. doi: 10.1042/BJ20130481
87. Zhang N, Zeng Y, Du W, Zhu J, Shen D, Liu Z, et al. The EGFR Pathway is Involved in the Regulation of PD-L1 Expression via the IL-6/JAK/STAT3 Signaling Pathway in EGFR-Mutated non-Small Cell Lung Cancer. Int J Oncol (2016) 49(4):1360–8. doi: 10.3892/ijo.2016.3632
88. Bu LL, Yu GT, Wu L, Mao L, Deng WW, Liu JF, et al. STAT3 Induces Immunosuppression by Upregulating PD-1/PD-L1 in HNSCC. J Dent Res (2017) 96(9):1027–34. doi: 10.1177/0022034517712435
89. Johnson DE, O'Keefe RA, Grandis JR. Targeting the IL-6/JAK/STAT3 Signalling Axis in Cancer. Nat Rev Clin Oncol (2018) 15(4):234–48. doi: 10.1038/nrclinonc.2018.8
90. Thomas DA, Massagué J. TGF-Beta Directly Targets Cytotoxic T Cell Functions During Tumor Evasion of Immune Surveillance. Cancer Cell (2005) 8(5):369–80. doi: 10.1016/j.ccr.2005.10.012
91. Attili I, Karachaliou N, Bonanno L, Berenguer J, Bracht J, Codony-Servat J, et al. STAT3 as a Potential Immunotherapy Biomarker in Oncogene-Addicted non-Small Cell Lung Cancer. Ther Adv Med Oncol (2018) 10:1758835918763744. doi: 10.1177/1758835918763744
92. Huang Q, Xia J, Wang L, Wang X, Ma X, Deng Q, et al. miR-153 Suppresses IDO1 Expression and Enhances CAR T Cell Immunotherapy. J Hematol Oncol (2018) 11(1):58. doi: 10.1186/s13045-018-0600-x
93. Liu Y, Zhou N, Zhou L, Wang J, Zhou Y, Zhang T, et al. IL-2 Regulates Tumor-Reactive CD8(+) T Cell Exhaustion by Activating the Aryl Hydrocarbon Receptor. Nat Immunol (2021) 22(3):358–69. doi: 10.1038/s41590-020-00850-9
94. Mo F, Yu Z, Li P, Oh J, Spolski R, Zhao L, et al. An Engineered IL-2 Partial Agonist Promotes CD8(+) T Cell Stemness. Nature (2021) 597(7877):544–8. doi: 10.1038/s41586-021-03861-0
95. Lai JH, Hung LF, Huang CY, Wu DW, Wu CH, Ho LJ. Mitochondrial Protein CMPK2 Regulates IFN Alpha-Enhanced Foam Cell Formation, Potentially Contributing to Premature Atherosclerosis in SLE. Arthritis Res Ther (2021) 23(1):120. doi: 10.1186/s13075-021-02470-6
96. Horton HM, Anderson D, Hernandez P, Barnhart KM, Norman JA, Parker SE. A Gene Therapy for Cancer Using Intramuscular Injection of Plasmid DNA Encoding Interferon Alpha. Proc Natl Acad Sci U.S.A. (1999) 96(4):1553–8. doi: 10.1073/pnas.96.4.1553
97. Mantovani A, Allavena P, Sica A, Balkwill F. Cancer-Related Inflammation. Nature (2008) 454(7203):436–44. doi: 10.1038/nature07205
98. Zhang YH, He K, Shi G. Effects of MicroRNA-499 On the Inflammatory Damage of Endothelial Cells During Coronary Artery Disease Via the Targeting of PDCD4 Through the NF-κβ/TNF-α Signaling Pathway. Cell Physiol Biochem (2017) 44(1):110–24. doi: 10.1159/000484588
99. Pikarsky E, Porat RM, Stein I, Abramovitch R, Amit S, Kasem S, et al. NF-kappaB Functions as a Tumour Promoter in Inflammation-Associated Cancer. Nature (2004) 431(7007):461–6. doi: 10.1038/nature02924
100. Zhang M, Zhu R, Zhang L. Triclosan Stimulates Human Vascular Endothelial Cell Injury via Repression of the PI3K/Akt/mTOR Axis. Chemosphere (2020) 241:125077. doi: 10.1016/j.chemosphere.2019.125077
Keywords: coronary artery disease (CAD), CD8+ T cells, hub genes, weighted gene co-expression network analysis (WGCNA), CIBERSORT, cancer
Citation: Zhao S, Wu Y, Wei Y, Xu X and Zheng J (2022) Identification of Biomarkers Associated With CD8+ T Cells in Coronary Artery Disease and Their Pan-Cancer Analysis. Front. Immunol. 13:876616. doi: 10.3389/fimmu.2022.876616
Received: 15 February 2022; Accepted: 24 May 2022;
Published: 21 June 2022.
Edited by:
Mohamed S. Abdel-Hakeem, Emory University, United StatesReviewed by:
Ahmed F El-Yazbi, Alexandria University, EgyptIshtiaque Ahammad, National Institute of Biotechnology, Bangladesh
Copyright © 2022 Zhao, Wu, Wei, Xu and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jialin Zheng, emhlbmdqbDE2M19jb21AMTYzLmNvbQ==
†These authors have contributed equally to this work and share first authorship