Home Epithelial-mesenchymal transition-related genes in coronary artery disease
Article Open Access

Epithelial-mesenchymal transition-related genes in coronary artery disease

  • Xiang Xu , Renchao Zou , Xiaoyong Liu , Jia Liu and Qianqian Su EMAIL logo
Published/Copyright: April 22, 2022

Abstract

Epithelial-mesenchymal transition (EMT) is critical in the development of coronary artery disease (CAD). However, landscapes of EMT-related genes have not been fully established in CAD. We identified the differentially expressed mRNAs and lncRNAs (DElncRNAs) from the Gene Expression Omnibus database. Pearson’s correlation analysis, the least absolute shrinkage and selection operator regression, and support vector machine reverse feature elimination algorithms were used to screen EMT-related lncRNAs. The cis–trans regulatory networks were constructed based on EMT-related lncRNAs. Quantitative real-time polymerase chain reaction was performed to validate the expression of EMT-related genes in a cohort of six patients with CAD and six healthy controls. We further estimated the infiltration of the immune cells in CAD patients with five algorithms, and the correlation between EMT-related genes and infiltrating immune cells was analyzed. We identified eight EMT-related lncRNAs in CAD. The area under curve value was greater than 0.95. The immune analysis revealed significant CD8 T cells, monocytes, and NK cells in CAD and found that EMT-related lncRNAs were correlated with these immune cell subsets. Moreover, SNAI2, an EMT-TF gene, was found in the trans-regulatory network of EMT-related lncRNAs. Further, we found SNAI2 as a biomarker for the diagnosis of CAD but it also had a close correlation with immune cell subsets in CAD. Eight EMT-related lncRNAs and SNAI2 have important significance in the diagnosis of CAD patients.

1 Introduction

Coronary artery disease (CAD) is a common public health problem, mainly occurring in people over 45 years of age. CAD is now the leading cause of death in the United States, accounting for one in six deaths alone [1]. The American Heart Association has said that cardiovascular disease causes more than 17.3 million deaths a year; by 2030, the number of deaths will exceed 23.6 million [2]. CAD is the most common type of cardiovascular disease. Its pathogenesis is due to coronary artery atherosclerotic lesions resulting in the narrowing of vascular lumen causing obstruction and reduced blood supply to the myocardium, hypoxia or necrosis, and eventually leading to heart failure [3,4]. However, the pathogenesis of CAD is complex, and there are no apparent symptoms in the early stage. The results of the myocardial enzyme spectrum can be negative, which is only manifested by abnormal ST-T segment changes in the exercise plate electrocardiogram. Although coronary angiography is the gold standard for diagnosing CAD, its high cost and technical requirements, reliance on specific equipment, and potential risk of radionuclide radiation have limited its use. Meanwhile, it is neither practical nor ethical to perform invasive coronary angiography on low-risk patients [5]. The cost of blood biomarker detection is low and easy to promote [6]. Therefore, it is vital to search for more potential biomarkers for the diagnosis and treatment of CAD based on blood sequencing data.

Epithelial-mesenchymal transition (EMT) is a biological process in which epithelial cells are transformed into cells with a mesenchymal phenotype through a specific procedure. During this process, endothelial cells gradually lose their morphology and function and acquire the phenotypic characteristics of mesenchymal cells such as proliferation, migration, and collagen synthesis [7]. Recent studies have demonstrated that EMT plays a pivotal physiological and pathological role in the development and structural remodeling of the myocardium, blood vessels, and valves, suggesting that EMT may be a worthwhile target for preventing and treating cardiovascular diseases. For example, endocardial EMT generates valvular cells necessary for heart valve formation and complete septal formation [8]. Epicardial EMT also generates cardiac fibroblasts, vascular smooth muscle cells, and surrounding cardiomyocytes necessary for cardiac muscle growth and coronary angiogenesis [9]. In CAD, endothelial cells participate in the formation of fibroblasts through the mesenchymal transformation of epithelial cells and promote cardiac fibrosis [10]. In addition, the EMT of endothelial cells plays a crucial role in the process of atherogenesis [11]. These studies suggest that EMT genes have an important significance in the field of cardiovascular disease.

Long-coding RNAs (lncRNAs) are noncoding RNAs with a length of more than 200 nucleotides. They have a wide range of functional activities, including RNA decay, genetic regulation of gene expression, RNA splicing, microRNA (miRNA) regulation, and protein folding [12]. LncRNAs also play an essential role in forming atherosclerosis, CAD, and heart failure [13]. For example, several genome-wide association studies have found that some single nucleotide polymorphisms located in lncRNA-ANRIL are closely related to the susceptibility to atherosclerosis and are also the sites with the most potent genetic susceptibility in CAD [14,15]. However, among all the identified lncRNAs, only a few have been verified as being involved in the regulation of EMT. For example, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is a lncRNA competing with miRNAs, directly interacts with oncogenes and proteins, and is involved in the activation of Wnt/β-catenin, PI3K/Akt/mTOR – these are typical EMT-related signal pathways [16,17]. A novel study indicated that MALAT1/miRNA-203/Wnt5a axis was a potential regulate mechanism for CAD [18]. Currently, the number of EMT-related lncRNAs in CAD research is small. How EMT-related lncRNAs regulate the formation and progression of CAD remains unclear.

Therefore, based on the previous research, we constructed two machine learning algorithms: least absolute shrinkage and selection operator (LASSO) regression algorithm and support vector machine reverse feature elimination (SVM-RFE) algorithm to screen out EMT-related diagnostic lncRNAs in CAD patients. Meanwhile, we constructed cis–trans regulatory networks based on EMT-related lncRNAs and explored the potential EMT gene of related molecules in the cis–trans network and the target drugs and structures. We also investigated the correlation between EMT-related diagnostic signatures and immune cell subsets by immune analysis. Through bioinformatics methods, an in-depth excavation of the EMT genes as having a promoting role in coronary atherosclerosis and the potential signal pathways and molecular mechanisms for the prevention and treatment of CAD, can provide a new train of thought and targets.

2 Methods and materials

2.1 Data collection

For our study, we downloaded the microarray gene expression profiling data of CAD from the Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) database with accession number GSE113079 [19]. The platform for GSE113079 was GPL20115, Agilent-067406 Human CBC lncRNA + mRNA microarray V4.0 (Probe name version), which contained peripheral blood mononuclear cells of 93 CAD patients and 48 healthy controls. Two hundred EMT-related genes were obtained from the Molecular Signatures Database (MsigDB, http://www.broad.mit.edu/gsea/msigdb/). Besides, 1,639 genes related to TFs were acquired from the database of The Human Transcription Factors (TFBS, http://tfbsdb.systemsbiology.net/).

2.2 Differentially expressed analysis

The limma package in R was used to identify the differentially expressed lncRNAs (DElncRNAs) and DE genes (DEGs) between CAD and normal samples (Tables S1 and S2). The lncRNAs/mRNAs met the selection standards of |(Fold change (FC))| >1.5 and false discovery rate (FDR) <0.01 and were considered as DElncRNAs/DEGs for further study.

2.3 Correlation analysis

By mating the listed 200 EMT-related genes in the MsigDB database, DE EMT genes for CAD were identified. Then, Pearson’s correlation analysis was operated between the harvested DE-EMTs and DElncRNAs expression data in samples to identify the EMT-related lncRNAs according to the correlation coefficient and P values (|Cor| >0.8 and P < 0.05) (Table S3).

2.4 Diagnostic value of EMT-lncRNAs and SNAI2

To explore the diagnostic ability of EMT-lncRNAs mentioned above, receiver operating characteristic (ROC) analysis was first performed using the R package pROC, and the EMT-related lncRNAs with area under curve (AUC) >0.95 were screened for further study. After filtration of EMT-related lncRNAs, candidate diagnostic lncRNAs for CAD were selected via an integrated analysis of two algorithms consisting of LASSO and SVM-RFE. Logistic regression was performed on diagnostic lncRNAs and the SNAI2 gene, respectively, to construct a logistic regression diagnostic model, and the bias residual diagram was drawn (Figure S1a and b). Five-fold cross validation was used to evaluate the performance of the diagnostic signature. Moreover, the diagnostic value of EMT-related lncRNAs and SNAI2 was assessed by ROC curve analysis using the pROC package in the R language.

2.5 EMT-related lncRNAs categorization

Based on modifications of the previous classification [20], we classified lncRNAs according to their gene positions related to the most proximal protein-coding genes. First, based on whether they intersect a protein-coding gene, the lncRNA genes were regarded as intergenic and genic. Furthermore, intergenic lncRNAs were categorized into two groups depending on whether they were their transcribed from the same or opposite strands: convergent (IC) and divergent (ID). Genic lncRNAs were separated into genic exonic (genic exonic same strand [GES] and genic exonic antisense [GEAS]), genic intronic (genic intronic same strand [GIS] and genic intronic antisense [GIAS]), or overlapping (genic overlapping same strand [GOS] and genic overlapping antisense [GOAS]) based on whether they overlapped with the exons or introns of a protein-coding gene.

2.6 The regulatory mechanisms of diagnostic EMT-lncRNAs

It is reported that lncRNAs regulated transcription of the genes nearby by acting in cis- and tans-manners. For the cis-regulation manner, we first selected the genes located on the same chromosome within a 300 kb region upstream or downstream of the lncRNAs. Subsequently, the Pearson analysis method was performed to analyze the correlation between the harvested lncRNAs and their corresponding genes under the selection criteria of |Cor| >0.3 and P < 0.05.

For trans prediction, we focused on the fact that lncRNAs might regulate the expression levels of TFs by the trans manner. After selecting the genes correlated with lncRNAs by the Pearson method (|Cor| >0.8 and P < 0.05), we further overlapped these genes with identified DEGs and TFs to obtain trans-regulated genes. A lncRNA–mRNA network that included EMT-lncRNAs, cis- and trans-regulated genes was constructed and visualized by the Cytoscape software.

2.7 Functional enrichment analysis

To explore the latent biological functions and pathways, gene ontology (GO) annotation and kyoto encyclopedia of genes and genomes (KEGG) pathway analyses were employed on the DE-EMTs, cis-regulated genes, trans-regulated genes, and SNAI2-regulated genes of CAD, respectively. The optional pathways related to CAD were predicted by the Comparative Toxicogenomics Database (CTD, http://ctdbase.org). Genes related to CAD were predicted by the DisGeNET database (https://www.disgenet.org/home/). The KEGG pathways both in the CAD database and KEGG analysis were introduced into the lncRNA–mRNA network to establish a lncRNA–mRNA-pathway network for CAD.

2.8 Immunity analysis and its correlation with key genes

We used cell type identification by estimating relative subsets of RNA transcripts algorithm (CIBERSORT) [21] for immune infiltration. R script downloaded from CIBERSORT website (https://cibersort.stanford.edu/). After obtaining the immune cell expression matrix according to the instructions on the CIBERSORT website, we used the “ggplot2” software package to create a cumulative histogram that visually showed the proportion of 22 immune cell infiltrates in CAD patients. We also used the “vioplot” package to draw violin plots showing differences in the expression of 22 infiltrating immune cells. We used the “corrplot” software package in R software to calculate Pearson correlation coefficients among immune cells, and displayed the results by using a correlation heat map. Pearson correlation coefficients and p value between identified key genes and infiltrating immune cells were calculated by “cor” and “Hmisc” software packages and then visualized by the “ggcorrplot” software package. In addition, single-sample gene set enrichment analysis (ssGSEA) [22], microenvironment cell population (MCP)-counter algorithm [23], EPIC [24], and QuanTIseq [25] algorithms were also used to compare and assess cellular components between the high SNAI2 and low SNAI2 gene expression groups. The differences in the immune response under different algorithms were uncovered using a Heatmap.

2.9 The drug–gene prediction

The genes cis- and trans-regulated by diagnostic lncRNAs were supposed to be the promising drug targets for searching drugs through the Drug–Gene Interaction database (DGIdb, https://dgidb.genome.wustl.edu/) that contained the DGI information of several databases [26]. The drug–gene network was visualized by the Cytoscape tool.

2.10 Study population

A total of six CAD patients and six healthy controls were recruited as a validation cohort in the Second Affiliated Hospital of Kunming Medical University, Kunming, China. This study is in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of the Second Affiliated Hospital of Kunming Medical University (No. PJ-2022-14). All subjects were aged >18 years. The diagnosis of CAD is based on the European Society of Cardiology criteria [27], which specifies at least one vascular lesion with a narrowing of lumen diameter greater than 50% using coronary angiography. Healthy subjects had no history of chronic disease and CAD. All participants signed informed consent after receiving a complete study explanation.

2.11 RNA extraction and quantitative real-time polymerase chain reaction (qPCR)

The expression of target genes was detected by qPCR to validate our microarray data. Briefly, EDTA vacuum anticoagulant vessels were used to collect 4–8 mL of peripheral blood from patients with CAD and healthy controls in the morning. The total RNA was extracted from peripheral blood cells of patients with CAD using TRIzol reagent (Invitrogen) and reverse-transcribed into cDNA (SweScript RT I Frist Strand cDNA Synthesis Kit, Servicebio, Wuhan, China). qPCR was performed using 2× Universal Blue SYBR Green qPCR Master Mix (Servicebio, Wuhan, China) following the manufacturer’s instructions. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was measured as an internal control. The reaction conditions were 95°C for 1 min, 40 cycles at 95°C for 20 s, 55°C for 20 s, and 72°C for 30 s. The melting curve was analyzed for each sample. The average value in each duplicate was used to calculate the relative amount of target genes using the 2−ΔΔCt method. All primers were synthesized by TSINGKE (TSINGKE Biotechnology CO., Ltd, Beijing, China). The primer sequences of the tested genes are shown in Table 1. The mRNA levels of lncRNAs between CAD and normal groups were compared using t test (P < 0.05) using R.

Table 1

Primers used for qPCR analysis of eight EMT-related lncRNAs and SNAI2

Gene names Position Primer sequence
CTD-2089N3.3 Forward GACGACCACCACGCAGGAG
Reverse GGAGAGGGGGAACAAGGCT
AC113167.2 Forward TTGGCTTCCTTACCTTGAGT
Reverse CCGATTATCCTTTTTCTTCC
LINC02747 Forward TGTGTGCGTCCTCCCTAAA
Reverse AGCAGAGACAGAGCCGGTT
RP11-1152H15.1 Forward GGTCCACACTGCTTTTATG
Reverse CTCTAGTCCTTCGGCTCAA
LINC02833 Forward CAGACAAAATCAAAATGGAG
Reverse TAAGGTGGGAGTAAGGAAAC
AC109460.4 Forward GACTGTCCCTAATTTCCCCT
Reverse TGTCTGACTCGTTCTCCTCA
LINC01775 Forward AGACAGGAGGGCTGGGGT
Reverse GGAGGCTGAGGCAGGAGA
RP11-103H7.3 Forward CAGTTTCTGCCCACAATGA
Reverse GTGTAACCTCTCCAGCCCT
SNAI2 Forward AAACTACAGCGAACTGGACA
Reverse ATAGAGATACGGGGAAATAA
H-GAPDH Forward CCCATCACCATCTTCCAGG
Reverse CATCACGCCACAGTTTCCC

2.12 Statistical analysis

The subcellular localization of diagnostic EMT-related lncRNAs was predicted by the LncLocator online tool [28]. The clusterProfiler package in R was utilized to perform GO and KEGG analyses. The statistical analyses were performed using R (R Foundation for Statistical Computing, Vienna, Austria, version 4.0.3) (http://www.R-project.org/) and GraphPad Prism (Version 8 GraphPad Software, La Jolla, CA, USA). P value <0.05 was considered as statistically significant.

  1. Ethics statement: The study was carried out in accordance with the recommendations of the Ethics Committee of the Second Affiliated Hospital of Kunming Medical University (No. PJ-2022-14) and was approved by said committee. All participants signed informed consent after receiving a complete study explanation.

  2. Informed consent: A preprint has previously been published and PREPRINT (Version 2) available at Research Square, and DOI is https://doi.org/10.21203/rs.3.rs-940366/v2.

3 Results

3.1 Identification of EMT-related genes

We performed the differentially expressed analysis on the GSE113079 data set. As shown in Figure 1a, 5,955 DElncRNAs were identified between CAD and normal samples under |(FC)| >1.5 and FDR <0.01 with 3,067 were upregulated and 2,888 downregulated. Meanwhile, we screened 2,868 DEGs between the two groups, including 1,540 upregulated and 1,328 downregulated DEGs (Figure 1c). The expressed levels of DElncRNAs and DEGs were shown in the heatmap plot and displayed in Figure 1b and d, respectively. The 32 DEGs related to EMT were generated by overlapping 200 EMT genes in the MsigDB database and preselected DE-EMTs, in which 21 were upregulated and 11 were downregulated (Figure 1e and f).

Figure 1 
                  Identification of DE-EMTs in CAD. (a) Volcano plot of lncRNAs expression between CAD and normal groups. (b) Volcano plot of mRNAs expression between CAD and normal groups. (c) Heatmap of DElncRNAs between CAD and normal groups. (d) Heatmap of DEGs between CAD and normal groups. (e) Venn diagram was used for the intersection of DEGs and EMT genes. (f) Heatmap of DE-EMTs between CAD and normal groups. (g) GO enrichment analysis of DE-EMTs: the larger the bubble and longer columns represent the more genes enriched in this function, the deeper the color of the bubble and bars, the smaller the P value. (h) KEGG enrichment analysis of DE-EMTs.
Figure 1

Identification of DE-EMTs in CAD. (a) Volcano plot of lncRNAs expression between CAD and normal groups. (b) Volcano plot of mRNAs expression between CAD and normal groups. (c) Heatmap of DElncRNAs between CAD and normal groups. (d) Heatmap of DEGs between CAD and normal groups. (e) Venn diagram was used for the intersection of DEGs and EMT genes. (f) Heatmap of DE-EMTs between CAD and normal groups. (g) GO enrichment analysis of DE-EMTs: the larger the bubble and longer columns represent the more genes enriched in this function, the deeper the color of the bubble and bars, the smaller the P value. (h) KEGG enrichment analysis of DE-EMTs.

To further reveal the potential mechanism of the role of the above 32 DE-EMT in CAD, we performed functional enrichment analysis by Metascape software. A total of 262 biological process (BP) terms, 26 molecular function (MF) terms, and 26 cellular component terms were enriched in the GO system (Figure 1g). We mainly focused on the enrichment results of GO-BP categories. The “response to wounding” was the most significantly enriched term, speculating that these genes may be involved in the regulation of the response after coronary artery damage. Next, “extracellular matrix organization” and “extracellular structure organization” were significantly enriched, suggesting that these genes may be associated with the acquisition of a mesenchymal phenotype by epithelial cells after EMT. The enrichment of terms related to cell adhesion (“regulation of cell adhesion,” “cell-matrix adhesion,” “positive regulation of cell adhesion,” “cell-substrate adhesion,” “regulation of cell-substrate adhesion,” etc.) and cytoskeletal (“actin cytoskeleton organization,” “regulation of actin cytoskeleton organization,” “regulation of cytoskeleton organization,” “positive regulation of cytoskeleton organization,” etc.) regulation and extracellular matrix disassembly could also be confirming the possibility of this speculation. Surprisingly, we found that terms related to the dynamic developmental processes of the vasculature (“blood vessel morphogenesis,” “blood vessel development,” “vasculature development,” “angiogenesis,” etc.) and cardio (“heart development,” “semi-lunar valve development,” “heart morphogenesis,” “heart valve morphogenesis,” “heart valve development,” etc.) were also significantly enriched, implying that DE-EMT may be involved in the process of CAD onset and development. Besides, terms related to the regulation of biological processes such as activation and differentiation of B cells, lymphocytes, myeloid leukocytes, and terms related to immune response (“negative regulation of immune system process,” “cell activation involved in immune response,” “humoral immune response,” etc.) were also closely related to these genes. The enrichment results on the GO system are available in Table S4. KEGG pathway enrichment showed that 32 DE-EMT were involved in a total of 24 pathways (Figure 1h; Table S5). Among them, “PI3K-Akt signaling pathway” and “focal adhesion” were the two most significantly enriched pathways. Moreover, these genes were also linked to multiple CAD-related cardiac disease pathways, such as “IL-17 signaling pathway,” “PI3K-Akt signaling pathway,” and “Fluid shear stress and atherosclerosis.” This evidence further suggested that the 32 CAD-related DE-EMT genes may play a role in the CAD process through a certain EMT mechanism.

3.2 Identification of EMT-related lncRNAs

Next, we screened the lncRNAs related to EMT in CAD with correlation analysis. According to the screening criteria, 1,141 EMT-related lncRNAs were identified. The Top50 EMT-lncRNAs are illustrated in Figure 2a. Genic lncRNAs stand for the largest category (57.2%) of EMT-related lncRNAs (GEAS = 7.8% and GES = 7%, GIAS = 19% and GIS = 18.4%, GOAS = 2.8% and GOS = 2.3%), following intergenic lncRNAs were 42.7% (IC = 19.9% and ID = 22.8%) (Figure 2b).

Figure 2 
                  Identification of EMT-related lncRNAs in CAD. (a) Correlation analyses for EMT-related lncRNAs. (b) Classification bar diagram of EMT-related lncRNAs.
Figure 2

Identification of EMT-related lncRNAs in CAD. (a) Correlation analyses for EMT-related lncRNAs. (b) Classification bar diagram of EMT-related lncRNAs.

3.3 Construction of an EMT-related lncRNAs diagnostic signature for CAD

To further detect the diagnostic ability of these EMT-related lncRNAs, the AUC value of each EMT-related lncRNAs was analyzed. Two hundred twenty-three EMT-related lncRNAs were screened with an AUC value above 0.95 (Table S6). LASSO regression analysis and SVM-RFE algorithm were used to identify the optimal diagnostic lncRNAs in the GSE113079 data set and establish the risk signature for CAD. Sixteen EMT-related lncRNAs were screened via the LASSO analysis, which intersected with 34 EMT-related lncRNAs obtained from the SVM-RFE algorithm to identify 11 diagnostic lncRNAs for CAD (Table S7, Figure 3a–e). After annotating the diagnostic lncRNAs using the Rsubread package in R, we obtained eight lncRNAs that were used to construct a diagnostic signature for CAD (Table S8), showing accuracy and specificity for the diagnosis of CAD (AUC = 1) (Figure 3f). Besides, the AUC value of each diagnostic lncRNAs was greater than 0.95, which exhibited a better ability to distinguish CAD patients from normal individuals (Figure 3g). Subcellular localization of each lncRNA determines the regulatory models. To investigate the subcellular localization of the diagnostic lncRNAs, we assessed LncLocator online platforms. We uncovered that these diagnostic lncRNAs were mainly located in the cytosol and cytoplasm (Figure S2)

Figure 3 
                  EMT-related lncRNAs diagnostic signature for CAD. (a) LASSO coefficient profiles of the 16 EMT-related lncRNAs selected by the optimal lambda. (b) The diagnostic signature selection of optimal parameter (lambda) in LASSO model. (c and d) Results of SVM-RFE algorithms: The point highlighted indicates the lowest error rate, and the corresponding EMT-related lncRNAs at this point are the best signature selected by SVM-RFE. (e) Venn diagram of overlap EMT-related lncRNAs selected by LASSO and SVM-RFE algorithms. (f) The ROC curve of predicted outcomes of eight EMT-related lncRNAs diagnostic signature by a logistic regression model. (g) ROC analysis results for eight EMT-related lncRNAs in CAD.
Figure 3

EMT-related lncRNAs diagnostic signature for CAD. (a) LASSO coefficient profiles of the 16 EMT-related lncRNAs selected by the optimal lambda. (b) The diagnostic signature selection of optimal parameter (lambda) in LASSO model. (c and d) Results of SVM-RFE algorithms: The point highlighted indicates the lowest error rate, and the corresponding EMT-related lncRNAs at this point are the best signature selected by SVM-RFE. (e) Venn diagram of overlap EMT-related lncRNAs selected by LASSO and SVM-RFE algorithms. (f) The ROC curve of predicted outcomes of eight EMT-related lncRNAs diagnostic signature by a logistic regression model. (g) ROC analysis results for eight EMT-related lncRNAs in CAD.

3.4 Establishment of cis- and trans-regulatory network

Previous studies indicated that lncRNAs regulated gene expression via local (cis) and long-distance (trans) mechanisms [29]. In this study, we identified seven diagnostic lncRNAs regulated their nearby genes via the cis-regulatory manner, except RP11-103H7.3 (Table 2). Among them, only CTD-2089N3.3 were significantly correlated with their corresponding gene EMB via Pearson analysis under |Cor| >0.3 and P value <0.05 (Figure 4a). Based on the median expression level of EMB, we divided the CAD patients in the GSE113079 data set into the high-expressed and low-expressed EMB groups. GSEA result suggested that several metabolic- and tumor-related pathways were associated with high-expressed EMB groups, including “fatty acid metabolism,” “pyrimidine metabolism,” “mTOR signaling pathway,” and “TGF BETA signaling pathway.” In contrast, “calcium signaling pathway,” “complement, and coagulation cascades,” “neuroactive ligand–receptor interaction,” and “olfactory transduction” were involved in the low-expressed EMB group (Figure S3).

Table 2

Seven EMT-related diagnostic lncRNAs regulated their nearby genes via the cis-regulatory manner

Gene lncRNA Symbol Cor P value Chromsome Distance
EMB p22710 CTD-2089N3.3 −0.46194075 8.12 × 10−9 chr5 130,172
MEF2C p40724_v4 AC113167.2 0.156713247 0.063484128 chr5 −294,003
MYEOV p2383 LINC02747 0.23717925 0.004628671 chr11 0
KCNK10 p4630 RP11-1152H15.1 −0.251517089 0.002625124 chr14 0
UNC79 p4940 LINC02833 −0.073789859 0.384523053 chr14 0
LAT p5961 AC109460.4 −0.244589931 0.003466511 chr16 −9,677
LMNB2 p8126 LINC01775 0.05339002 0.529495266 chr19 −1,974
Figure 4 
                  The cis and trans network of the EMT-related lncRNAs. (a) Cis-regulation gene of lncRNA CTD-2089N3.3 in the chromosome. The X-axis represents lncRNA position in chromosome and the Y-axis represents correlation coefficient of lncRNA and the potential “cis” gene. The red line represents the genome width of lncRNA and blue point represents the position of potential “cis” gene. (b) Heatmap of the correlations between eight EMT-related lncRNAs and trans-regulated genes. (c) Venn diagram of 33 “trans” genes regulated by eight EMT-related lncRNAs. (d) The cis and trans network of the eight EMT-related lncRNAs (red) in CAD and their target mRNAs (green). (e and f) KEGG enrichment analysis of 34 cis and trans genes regulated by EMT-related lncRNAs.
Figure 4

The cis and trans network of the EMT-related lncRNAs. (a) Cis-regulation gene of lncRNA CTD-2089N3.3 in the chromosome. The X-axis represents lncRNA position in chromosome and the Y-axis represents correlation coefficient of lncRNA and the potential “cis” gene. The red line represents the genome width of lncRNA and blue point represents the position of potential “cis” gene. (b) Heatmap of the correlations between eight EMT-related lncRNAs and trans-regulated genes. (c) Venn diagram of 33 “trans” genes regulated by eight EMT-related lncRNAs. (d) The cis and trans network of the eight EMT-related lncRNAs (red) in CAD and their target mRNAs (green). (e and f) KEGG enrichment analysis of 34 cis and trans genes regulated by EMT-related lncRNAs.

By combining 685 genes correlated with lncRNAs with DEGs and identified TFs, 33 genes were identified to be regulated by diagnostic lncRNAs via trans manner, in which SNAI2 was founded to be a DE-EMT (Figure 4b and c). Then, a lncRNA–mRNA regulatory network was constructed that contained diagnostic lncRNAs, cis- and trans-regulated genes, which consisted of 42 nodes and 93 edges (Figure 4d). These genes in the regulatory network were mainly involved in nervous development and vitamin metabolism by GO analysis (Figure S4a). Combining the pathways related to these genes identified with the KEGG analysis and related to CAD development in the CTD database, “maturity onset diabetes of the young” and “transcriptional misregulation in cancer” were discovered (Figure 4e and f, Table 3). Hence, these two KEGG pathways were introduced into the lncRNA–mRNA regulatory network to establish a lncRNA–mRNA pathway network for CAD that included 46 nodes and 98 edges (Figure S4b).

Table 3

Combining the pathways related to potential “cis” and “trans” genes identified with the KEGG analysis and related to CAD development in the CTD database

# Input DiseaseName PathwayName PathwayID Inference GeneSymbol
Coronary heart disease Coronary disease Maturity onset diabetes of the young KEGG:hsa0495 0 HNF1A
Coronary heart disease Coronary disease Transcriptional misregulation in cancer KEGG:hsa0520 2 MMP3
Coronary heart disease Coronary disease Transcriptional misregulation in cancer KEGG:hsa0520 2 PLAU

3.5 Prediction of regulatory genes of SNAI2

Based on the above results, we uncovered that SNAI2 was found to be a DE-EMT and TFs among all genes regulated by diagnostic lncRNAs. In our study, SNAI2 is obviously higher expressed in CAD groups than normal groups (p = 5.8 × 10–15; Figure 5a). Considering the importance of SNAI2, we also detect the diagnostic ability of SNAI2 in CAD patients. The AUC value of SNAI2 was 0.902 (Figure 5b). In addition, we used five-fold cross validation to evaluate the reliability of the SNAI2 gene. First, we randomly divided the samples into five parts, of which four parts were used as training sets to build the logistic regression model, and the rest were used to verify the model. This process was then repeated five times to reduce errors and improve the sensitivity of the model. The AUC values of the five models were 0.9479, 0.9144, 0.9391, 0.7692, and 0.8766, respectively, indicating that the models had good explanatory power (Figure 5c). Besides, GSEA was performed to investigate the latent biological functions. “Calcium signaling pathway,” “linoleic acid metabolism,” “neuroactive ligand receptor interaction,” and “olfactory transduction” were mainly associated with the high-expressed SNAI2 group. “RNA degradation,” “splicesome,” “fatty acid metabolism,” and “histone metabolism” were involved in the low-expressed SNAI2 group (Figure S5). Moreover, we overlapped 234 genes regulated by SNAI2 acquired from the TFBS database and 1,576 genes related to CAD acquired from the DisGeNET database to obtain 21 genes regulated by SNAI2 for CAD (Figure 5d). Functional enrichment analysis determined that the harvested 21 genes were concerted on the “insulin secretion,” “peptide hormone secretion,” and “long-chain fatty acid biosynthetic process” (Figure 5e). No pathways were detected by the KEGG analysis.

Figure 5 
                  Clinical performance of SNAI2 in CAD. (a) SNAI2 is highly expressed in CAD groups and low in normal groups. (b) ROC analysis results for SNAI2 in CAD. (c) The ROC curve of predicted outcomes of SNAI2 diagnostic signature by a logistic regression model. (d) Venn diagram of 21 CAD-related genes regulated by SNAI2. (e) GO enrichment analysis of 21 CAD-related genes regulated by SNAI2.
Figure 5

Clinical performance of SNAI2 in CAD. (a) SNAI2 is highly expressed in CAD groups and low in normal groups. (b) ROC analysis results for SNAI2 in CAD. (c) The ROC curve of predicted outcomes of SNAI2 diagnostic signature by a logistic regression model. (d) Venn diagram of 21 CAD-related genes regulated by SNAI2. (e) GO enrichment analysis of 21 CAD-related genes regulated by SNAI2.

3.6 Prediction of the target drugs of genes in the cis- and trans-regulatory network

Next, the target drugs of genes regulated by diagnostic lncRNAs were predicted by the DGIdb database. Through DGIdb prediction, a total of 483 drug–gene pairs were identified, and a target-drug network for CAD was constructed, including five genes and 476 drugs. Four hundred fifty-nine drugs interacted with VDR, which might be promising to treat patients with CAD (Figure S6a). The structures of these drugs are illustrated in Figure S6b–i.

3.7 Immune analysis of EMT-lncRNAs and SNAI2

Enrichment analysis showed that DE-EMT gene was enriched in inflammatory response-related pathways. Therefore, we evaluated the type and fraction of immune cell infiltration between CAD patients and normal samples in the data set using the CIBERSORT algorithm. The relative proportion of immune cell subtypes is shown in the cumulative histogram (Figure 6a). Our results found an apparent proportion of CD8 T cells, NK cells activated, and monocytes. Moreover, the infiltration of CD8 T cells and NK cell activated were decreased, and the infiltration of monocytes was increased in CAD patients (Figure 6b). By principal component analysis (PCA), immune cell fractions in CAD patients and normal controls showed intergroup bias and individual differences (Figure 6c). In the correlation heatmap (Figure 6d), we found that CD8 T cells were negatively correlated with monocytes and macrophages M0, and positively correlated with NK cells activated. This is consistent with the correlation between seven EMT-related lncRNAs and the immune cells we found, except lncRNA AC109460.4 (Figure 7a). In addition, we also conducted an immune analysis of SNAI2. Then, we divided the samples into high and low groups according to the expression level of SNAI2. It was found that in the high expression level group, the infiltration of monocytes was decreased. In contrast, the infiltration of NK cells activated and CD8 T cells were increased, which was similar to immune cell infiltration in CAD patients (Figure 7b). The heatmap of immune cell compositions based on CIBERSORT, quanTIseq, ssGSEA, MCP-counter, and EPIC algorithms are shown in Figure 8. It was found that CD8 T cells, monocytes, and NK cells activated had similar immune cell infiltration trends in the CAD and SNAI2 gene high expression groups.

Figure 6 
                  Evaluation and visualization of immune cell infiltration by CIBERSORT algorithm. (a) Immune cell types and ratios of CAD groups. (b) Violin plot comparing immune cell compositions between CAD groups and normal groups. (c) PCA for immune cell compositions between CAD groups and normal groups. (d) Pearson correlation heatmap between infiltrating immune cell subpopulations.
Figure 6

Evaluation and visualization of immune cell infiltration by CIBERSORT algorithm. (a) Immune cell types and ratios of CAD groups. (b) Violin plot comparing immune cell compositions between CAD groups and normal groups. (c) PCA for immune cell compositions between CAD groups and normal groups. (d) Pearson correlation heatmap between infiltrating immune cell subpopulations.

Figure 7 
                  Immune analysis of SNAI2. (a) Pearson correlation heatmap between eight EMT-related diagnostic signatures and infiltrating immune cells. (b) Boxplot comparing immune cell compositions between high SNAI2 expression level groups and low SNAI2 expression level groups.
Figure 7

Immune analysis of SNAI2. (a) Pearson correlation heatmap between eight EMT-related diagnostic signatures and infiltrating immune cells. (b) Boxplot comparing immune cell compositions between high SNAI2 expression level groups and low SNAI2 expression level groups.

Figure 8 
                  Heatmap for immune cell compositions based on ssGSEA, CIBERSORT, quanTIseq, MCP-counter, and EPIC algorithms among different groups.
Figure 8

Heatmap for immune cell compositions based on ssGSEA, CIBERSORT, quanTIseq, MCP-counter, and EPIC algorithms among different groups.

3.8 Validation of EMT-lncRNAs and SNAI2 expression

To assess the accuracy of our results, qPCR was used to detect the expression of eight EMT-related lncRNAs and SNAI2 in peripheral blood samples from six CAD patients and six normal patients (Figure 9). The results showed that the expression patterns of eight EMT-related lncRNAs and SNAI2 were consistent with microarray data in peripheral blood of CAD patients and normal controls. The results show that CTD-2089N3.3 (P = 0.0154 < 0.05), AC113167.2 (P = 0.0170 < 0.05), LINC02747 (P = 0.0166 < 0.05), RP11-1152H15.1 (P = 0.0151 < 0.05), LINC02833 (P = 0.0164 < 0.05), LINC01775 (P = 0.0159 < 0.05), RP11-103H7.3 (P = 0.0164 < 0.05), and SNAI2 (P = 0.0212 < 0.05) were significantly higher in the CAD groups than in the normal control groups. The expression level of AC109460.4 was increased in patients with normal control groups compared with CAD patients (P = 0.0133 < 0.05).

Figure 9 
                  The relative expression of eight EMT-related lncRNAs and SNAI2 in the validation cohort: (a) CTD-2089N3.3, (b) AC113167.2, (c) LINC02747, (d) RP11-1152H15.1, (e) LINC02833, (f) LINC01775, (g) RP11-103H7.3, (h) SNAI2, and (i) AC109460.4.
Figure 9

The relative expression of eight EMT-related lncRNAs and SNAI2 in the validation cohort: (a) CTD-2089N3.3, (b) AC113167.2, (c) LINC02747, (d) RP11-1152H15.1, (e) LINC02833, (f) LINC01775, (g) RP11-103H7.3, (h) SNAI2, and (i) AC109460.4.

4 Discussion

EMT plays a critical physiological and pathological role in developing the cardiovascular system, vascular tissue remodeling, and heart valve disease during the embryonic period [30]. However, more research has been focused on the impact of the EMT in tumor development and treatment. In contrast, few studies have explored the diagnostic value of EMT-related genes or lncRNAs in CAD. Hence, exploring diagnostic biomarkers of EMT-related genes or lncRNAs in CAD is urgent.

Our analyses uncovered 32 EMT-related DEGs in CAD. KEGG pathway analysis of these DE-EMTs was mainly enriched in the PI3K-Akt signaling pathway. Several reports have shown that PI3K-Akt pathway is a canonical EMT signaling pathway [31,32]. Meanwhile, we found this signaling pathway plays an essential role in CAD. A recent study indicated that miRNA-26a-5p activated the PI3K-Akt pathway by targeting phosphatase and tensin homolog (PTEN) and affected the proliferation and apoptosis of endothelial cells isolated from CAD mice [33]. A comparative study also reported that miR-26a-5p could activate the PI3K-Akt signaling pathway through the inhibition of PTEN, thereby protecting against myocardial defect/reperfusion injury [34]. These studies have confirmed that activating the PI3K-Akt signaling pathway can prevent myocardial ischemia-reperfusion in animal models. Other studies have also suggested that regulation of the PI3K-Akt signaling pathway plays a vital role in inhibiting myocardial fibrosis, apoptosis, and the inflammatory response [35,36].

In this study, we performed a coexpression analysis between EMT genes and DElncRNAs through paired lncRNA and mRNA expression data in CAD patients from GEO. Eight differently expressed EMT-related lncRNAs were found to be diagnosis factor for CAD patients. After a literature review, we found no research had been conducted about the mechanisms of the eight lncRNAs except LINC02747. Previous studies have reported that LINC02747 can upregulate the expression of TFE3 by absorbing miR-608 and ultimately promote the proliferation of clear cell renal cell carcinoma (ccRCC) [37]. Gu et al. indicated that miR-608 exerts anti-inflammatory effects by targeting ELANE in monocytes [38]. Our results showed that monocytes were more expressed in the CAD group, so whether the regulation of LINC02747-mir608-ELANE might achieve the reversal of inflammatory response in CAD patients is not clear. Other seven EMT-related lncRNAs have not been reported in relevant studies, and reports on how lncRNAs interact with EMT genes have been rarer. However, many “cis” and “trans” genes are involved in the formation and development of CAD in the cis-trans regulatory network. For example, EMB, as a “cis” gene, was enriched in the mTOR signaling pathway in our GSEA analysis. This pathway is closely associated with atherosclerosis, and the pro-inflammatory response of monocytes in CAD requires activation of mTOR [39]. Among “trans” genes, some studies have reported that VDR gene polymorphisms leads to the development and formation of CAD by affecting changes in serum levels of 25(OH) vitamin D. [40,41]. A previous study reported VDR in regulating inflammation through inhibiting the NF-ĸB pathway and activating autophagy [42]. EBF4 gene promotes the elevation of Cu and leads to the progression of CAD by affecting copper-related DNA methylation sites [43]. CTCF gene is essential for cardiogenesis and to inhibit cardiomyocytes apoptosis, and can be applied as a therapeutic target for the treatment of heart failure in future [44,45]. FLI1 gene is also reported to be closely related to immune dysfunction and platelet disorders [46]. Despite the lack of direct support from literature, we speculated that these cis–trans genes, under the regulation of lncRNA, affect the formation and development of CAD through immune microenvironment, cell apoptosis, platelet dysfunction, and other ways. So far, there has been no study on the role of EMT-related lncRNA in CAD diagnosis. In our study, EMT-related genes with high specificity were identified by bioinformatics methods, and these genes in CAD groups and normal groups were validated by qPCR method. These findings may provide valuable insights into the future diagnosis and treatment of CAD.

The presence of immune cells in the infarct area is vital for initiating the repair process of the injured heart tissue. Temporal and spatial regulation of inflammation after infarction is crucial [47,48]. We evaluated the type and fraction of immune cell infiltration between the CAD patients and normal samples in the data set using the CIBERSORT algorithm. Our results found CD8 T cells and NK cells share a decreased infiltration, and the infiltration of monocytes was increased in CAD patients, which was similar to the previous results [49,50,51]. In this GEO data set, CD8 T cells and NK cells are favorable factors for preventing CAD, and it is likely that monocytes promote the occurrence of CAD. Previous studies have suggested that the imbalance of immune regulation is an essential factor for promoting atherosclerosis, heart failure, and chronic kidney disease by monocytes cells [48]. CD8 T cells play a dual role in atherosclerosis. On the one hand, CD8 T cells can secret many inflammatory cytokines to accelerate the inflammatory response and increase the instability of atherosclerotic plaques. On the other hand, cytotoxic activity against antigen-presenting cells and the presence of regulatory CD8 T cell subsets could suppress immunity and limit atherosclerosis [52]. Ong et al. suggested that NK cells appear to protect the development of cardiac fibrosis by preventing the accumulation of specific inflammatory groups in the heart and directly restricting collagen formation in cardiac fibroblasts [53]. Although the results of our study are similar to these researches, the mechanism of the immune system is still very complex and some results in the immunotherapy of CAD are not ideal. We need a lot of clinical studies to demonstrate the underlying mechanism. Besides, we also found that except AC109460.4, the other seven lncRNAs related to EMT were significantly negatively correlated with CD8 T cell and NK cell and positively correlated with Treg and monocytes. The results of AC109460.4 were just the opposite. The association between these lncRNAs and the innate immune system is still unclear. More in vivo and in vitro studies are needed to explain the interaction mechanism between these lncRNAs and immune cells in CAD.

It is generally believed that lncRNAs can act in “trans” to regulate TFs mediated chromatin remodeling and transcription [54]. These lncRNAs recruit protein factors to enhancer and regulate enhancer activity [55]. We constructed cis- and trans-regulatory networks based on these eight signatures. In the trans-regulatory network, we obtained 33 differentially expressed TF genes. The most surprising discovery was the screening of SNAI2, an EMT-TF gene (the gene coding product was the transcription factor Snai2). Our results indicated that SNAI2 was not only significantly highly expressed in CAD patients but also strongly positively correlated with LINC01775 and CTD-2089N3.3. The ROC curve showed that the SNAI2 could be a potential biomarker for diagnosing CAD. Additionally, we also validated the high expression of SNAI2 gene in the CAD group by using the qPCR method. As a classic EMT-TF gene, SNAI2 has recently been shown to be involved in a broader range of biological processes, including tumor metastasis, heart development, cell differentiation, vascular remodeling, and DNA damage repair [56,57,58]. Previous studies have reported that the deletion of protein arginine methyltransferase 1 leads to the accumulation of p53, and enhancing the degradation of SNAI2 can limit the formation of cardiac fibroblasts, coronary smooth muscle cells, and pericytes [59]. Meanwhile, Cooley et al. reported that, by grafting mouse veins to the femoral artery in mice to simulate human coronary artery bypass grafting, the results showed that TGF-β/Smad2/3-Snai2-mediated EMT plays a crucial role in venous graft vessel remodeling [60]. These studies have indicated that high expression of SNAI2 can promote the formation of vascular endothelium to EMT and vascular remodeling, which is one of the vital factors in the formation of CAD. At present, the role of SNAI2 in CAD has not been reported. Several studies have proven that the vascular endothelial EMT process is involved in atherosclerosis, post-stent stenosis, pulmonary hypertension, and coronary artery remodeling [61,62,63]. Additionally, the role of EMT can be seen in a range of cell types involved in immunity, such as lymphocytes, NK cells, and myeloid cells, which contribute to inflammatory responses in diverse pathophysiological processes. Ricciardi et al. have reported a decreased viability and proliferation of NK cells and T cells after coculture with cancer cell lines in which EMT had been induced [64]. In our study, SNAI2 correlated with infiltration of monocytes, CD8 T cells, and NK cells was activated. Previous studies have suggested that SNAI2 deletion in mice leads to impaired development of the T-lymphatic system [65]. Subsequent studies also confirmed that Snai2 plays a vital influence in regulating CD8 T cells and targets genes with functions for T cells [66]. Furthermore, our results indicated that the difference in these immune cell infiltrations in the SNAI2 high expression group was similar to the results of CAD patients. These immune cells have been researched to play a role in the formation, erosion, and rupture of coronary plaques [67,68]. In summary, we inferred that SNAI2 might play a significant role in the occurrence of CAD by regulating innate and adaptive immunity through these immune cells. To confirm our conclusion, more experimental mechanistic research should be carried out in future studies.

Our study should acknowledge some limitations. First, these EMT-related lncRNAs were investigated in data sets with no access to individual patients’ characteristics; thus, we cannot adjust the ROC curve for traditional cardiovascular risk factors. A prospective cohort recruiting CAD patients is needed to confirm the predictive value of EMT-related lncRNAs. Second, the MF details of SNAI2 and EMT-related lncRNAs in the progression of CAD have not been further studied. Therefore, molecular biological experiments and flow cytometry analysis are required to validate these findings, and another external validation based on a larger sample is needed.

5 Conclusion

In conclusion, this comprehensive bioinformatic analysis revealed that SNAI2 and EMT-related lncRNAs could be reliable biomarkers for diagnosing CAD and may be used for decision-making in the treatment of CAD patients. At the same time, based on the eight EMT-related lncRNAs, we constructed the cis- and trans-regulatory networks of CAD. Furthermore, the immune analysis suggested that these biomarkers were closely related to immune cells and CAD. These findings provide references for clinicians to understand the molecular mechanism of interaction between CAD and EMT and develop individualized treatment for CAD patients.

Abbreviations

CAD

coronary artery disease

EMT

epithelial-mesenchymal transition

lncRNAs

long-coding RNAs

MALAT1

metastasis associated lung adenocarcinoma transcript 1

miRNAs

microRNAs

SNPs

single nucleotide polymorphisms

SVM-RFE

support vector machine reverse feature elimination

GSEA

gene set enrichment analysis

ssGSEA
BP

biological processes

MF

molecular function

CC

cellular components

KEGG

kyoto encyclopedia of genes and genomes

DEGs

differentially expressed genes

DE-EMTs

differentially expressed EMT genes

DElncRNAs

differentially expressed lncRNAs

FC

fold-change

CIBERSORT

cell type identification by estimating relative subsets of RNA transcripts algorithm

PCA

principal component analysis

ccRCC

clear cell renal cell carcinoma

GEO

gene expression omnibus

LASSO

least absolute shrinkage and selection operator

FDR

false discovery rate

AUC

the area under the curve

ROC

receiver operating characteristic

PRMT1

protein arginine methyltransferase 1

CABG

coronary artery bypass grafting.


# These authors contributed equally to this work.


Acknowledgments

Not applicable.

  1. Funding information: This study was supported by the Scientific Research Project of The Second Affiliated Hospital of Kunming Medical University (No. 2021yk013).

  2. Author contributions: Xiang Xu designed the experiments, performed the experiments, analyzed the data, prepared figures, reviewed drafts of the paper, approved the final draft. Qianqian Su designed the experiments, performed the experiments, reviewed drafts of the paper, approved the final draft. Renchao Zou analyzed the data, contributed analysis tools, prepared figures, approved the final draft. Xiaoyong Liu analyzed the data, prepared figures, approved the final draft. Jia Liu performed the experiments, reviewed drafts of the paper, approved the final draft.

  3. Conflict of interest: The authors declare that they have no competing interests.

  4. Data availability statement: The data used to support the findings of this study are available from the corresponding author upon request. We thank the contributors of the GEO databases for the availability of the data. The gene expression profiles in the GSE113079 dataset were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE113079).

References

[1] Adapala RK, Kanugula AK, Paruchuri S, Chilian WM, Thodeti CK. TRPV4 deletion protects heart from myocardial infarction-induced adverse remodeling via modulation of cardiac fibroblast differentiation. Basic Res Cardiol. 2020;115(2):14. 10.1007/s00395-020-0775-5.Search in Google Scholar PubMed PubMed Central

[2] Writing Group M, Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, et al. Heart disease and stroke statistics-2016 update: a report from the American Heart Association. Circulation. 2016;133(4):e38-360. 10.1161/CIR.0000000000000350.Search in Google Scholar PubMed

[3] Li GM, Zhang CL, Rui RP, Sun B, Guo W. Bioinformatics analysis of common differential genes of coronary artery disease and ischemic cardiomyopathy. Eur Rev Med Pharmacol Sci. 2018;22(11):3553–69. 10.26355/eurrev_201806_15182.Search in Google Scholar PubMed

[4] Chen CH, Wang SS, Wei EI, Chu TY, Hsieh PC. Hyaluronan enhances bone marrow cell therapy for myocardial repair after infarction. Mol Ther. 2013;21(3):670–9. 10.1038/mt.2012.268.Search in Google Scholar PubMed PubMed Central

[5] De Marco E, Vacchiano G, Frati P, La Russa R, Santurro A, Scopetti M, et al. Evolution of post-mortem coronary imaging: from selective coronary arteriography to post-mortem CT-angiography and beyond. Radiol Med. 2018;123(5):351–8. 10.1007/s11547-018-0855-x.Search in Google Scholar PubMed

[6] Rusnak J, Fastner C, Behnes M, Mashayekhi K, Borggrefe M, Akin I. Biomarkers in stable coronary artery disease. Curr Pharm Biotechnol. 2017;18(6):456–71. 10.2174/1389201018666170630120805.Search in Google Scholar PubMed

[7] Georgakopoulos-Soares I, Chartoumpekis DV, Kyriazopoulou V, Zaravinos A. EMT factors and metabolic pathways in cancer. Front Oncol. 2020;10:499. 10.3389/fonc.2020.00499.Search in Google Scholar PubMed PubMed Central

[8] von Gise A, Pu WT. Endocardial and epicardial epithelial to mesenchymal transitions in heart development and disease. Circ Res. 2012;110(12):1628–45. 10.1161/CIRCRESAHA.111.259960.Search in Google Scholar PubMed PubMed Central

[9] Zhou B, von Gise A, Ma Q, Hu YW, Pu WT. Genetic fate mapping demonstrates contribution of epicardium-derived cells to the annulus fibrosis of the mammalian heart. Dev Biol. 2010;338(2):251–61. 10.1016/j.ydbio.2009.12.007.Search in Google Scholar PubMed PubMed Central

[10] Wilhelmi T, Xu X, Tan X, Hulshoff MS, Maamari S, Sossalla S, et al. Serelaxin alleviates cardiac fibrosis through inhibiting endothelial-to-mesenchymal transition via RXFP1. Theranostics. 2020;10(9):3905–24. 10.7150/thno.38640.Search in Google Scholar PubMed PubMed Central

[11] Souilhol C, Harmsen MC, Evans PC, Krenning G. Endothelial-mesenchymal transition in atherosclerosis. Cardiovasc Res. 2018;114(4):565–77. 10.1093/cvr/cvx253.Search in Google Scholar PubMed

[12] Khan S, Masood M, Gaur H, Ahmad S, Syed MA. Long non-coding RNA: an immune cells perspective. Life Sci. 2021;271:119152. 10.1016/j.lfs.2021.119152.Search in Google Scholar PubMed

[13] Uchida S, Dimmeler S. Long noncoding RNAs in cardiovascular diseases. Circ Res. 2015;116(4):737–50. 10.1161/CIRCRESAHA.116.302521.Search in Google Scholar PubMed

[14] Cunnington MS, Santibanez Koref M, Mayosi BM, Burn J, Keavney B. Chromosome 9p21 SNPs associated with multiple disease phenotypes correlate with ANRIL expression. PLoS Genet. 2010;6(4):e1000899. 10.1371/journal.pgen.1000899.Search in Google Scholar PubMed PubMed Central

[15] Broadbent HM, Peden JF, Lorkowski S, Goel A, Ongen H, Green F, et al. Susceptibility to coronary artery disease and diabetes is encoded by distinct, tightly linked SNPs in the ANRIL locus on chromosome 9p. Hum Mol Genet. 2008;17(6):806–14. 10.1093/hmg/ddm352.Search in Google Scholar PubMed

[16] Uthman YA, Ibrahim KG, Abubakar B, Bello MB, Malami I, Imam MU, et al. MALAT1: A promising therapeutic target in metastatic colorectal cancer. Biochem Pharmacol. 2021;190:114657. 10.1016/j.bcp.2021.114657.Search in Google Scholar PubMed

[17] Xu W, Ding M, Wang B, Cai Y, Guo C, Yuan C. Molecular mechanism of the canonical oncogenic lncRNA MALAT1 in gastric cancer. Curr Med Chem. 2021;28:8800–9. 10.2174/0929867328666210521213352.Search in Google Scholar PubMed

[18] Cai Q, Gao ML, Huang LS, Chen HS, Pan LH. MALAT1/miRNA-203/Wnt5a: a potential mechanism for regulating coronary artery disease. Int J Cardiol. 2021;329:48. 10.1016/j.ijcard.2020.12.046.Search in Google Scholar PubMed

[19] Li L, Wang L, Li H, Han X, Chen S, Yang B, et al. Characterization of LncRNA expression profile and identification of novel LncRNA biomarkers to diagnose coronary artery disease. Atherosclerosis. 2018;275:359–67. 10.1016/j.atherosclerosis.2018.06.866.Search in Google Scholar PubMed

[20] Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22(9):1775–89. 10.1101/gr.132159.111.Search in Google Scholar PubMed PubMed Central

[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. 10.1038/nmeth.3337.Search in Google Scholar PubMed PubMed Central

[22] Yi M, Nissley DV, McCormick F, Stephens RM. ssGSEA score-based Ras dependency indexes derived from gene expression data reveal potential Ras addiction mechanisms with possible clinical implications. Sci Rep. 2020;10(1):10258. 10.1038/s41598-020-66986-8.Search in Google Scholar PubMed PubMed Central

[23] Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218. 10.1186/s13059-016-1070-5.Search in Google Scholar PubMed PubMed Central

[24] Racle J, de Jonge K, Baumgaertner P, Speiser DE, Gfeller D. Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. Elife. 2017;6. 10.7554/eLife.26476.Search in Google Scholar PubMed PubMed Central

[25] Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019;11(1):34. 10.1186/s13073-019-0638-6.Search in Google Scholar PubMed PubMed Central

[26] Freshour SL, Kiwala S, Cotto KC, Coffman AC, McMichael JF, Song JJ, et al. Integration of the drug-gene interaction database (DGIdb 4.0) with open crowdsource efforts. Nucleic Acids Res. 2021;49(D1):D1144-D51. 10.1093/nar/gkaa1084.Search in Google Scholar PubMed PubMed Central

[27] Task Force M, Montalescot G, Sechtem U, Achenbach S, Andreotti F, Arden C, et al. 2013 ESC guidelines on the management of stable coronary artery disease: the task force on the management of stable coronary artery disease of the European Society of Cardiology. Eur Heart J. 2013;34(38):2949–3003. 10.1093/eurheartj/eht296.Search in Google Scholar PubMed

[28] Cao Z, Pan X, Yang Y, Huang Y, Shen HB. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics. 2018;34(13):2185–94. 10.1093/bioinformatics/bty085.Search in Google Scholar PubMed

[29] Chen LL. Linking long noncoding RNA localization and function. Trends Biochem Sci. 2016;41(9):761–72. 10.1016/j.tibs.2016.07.003.Search in Google Scholar PubMed

[30] Liu X, Wang Y, Liu F, Zhang M, Song H, Zhou B, et al. Wdpcp promotes epicardial EMT and epicardium-derived cell migration to facilitate coronary artery remodeling. Sci Signal. 2018;11(519). 10.1126/scisignal.aah5770.Search in Google Scholar PubMed

[31] Dong WQ, Chao M, Lu QH, Chai WL, Zhang W, Chen XY, et al. Prohibitin overexpression improves myocardial function in diabetic cardiomyopathy. Oncotarget. 2016;7(1):66–80. 10.18632/oncotarget.6384.Search in Google Scholar PubMed PubMed Central

[32] Xu Z, Jia K, Wang H, Gao F, Zhao S, Li F, et al. METTL14-regulated PI3K/Akt signaling pathway via PTEN affects HDAC5-mediated epithelial-mesenchymal transition of renal tubular cells in diabetic kidney disease. Cell Death Dis. 2021;12(1):32. 10.1038/s41419-020-03312-0.Search in Google Scholar PubMed PubMed Central

[33] Jing R, Zhong QQ, Long TY, Pan W, Qian ZX. Downregulated miRNA-26a-5p induces the apoptosis of endothelial cells in coronary heart disease by inhibiting PI3K/AKT pathway. Eur Rev Med Pharmacol Sci. 2019;23(11):4940–7. 10.26355/eurrev_201906_18084.Search in Google Scholar PubMed

[34] Xing X, Guo S, Zhang G, Liu Y, Bi S, Wang X, et al. miR-26a-5p protects against myocardial ischemia/reperfusion injury by regulating the PTEN/PI3K/AKT signaling pathway. Braz J Med Biol Res. 2020;53(2):e9106. 10.1590/1414-431X20199106.Search in Google Scholar PubMed PubMed Central

[35] Li X, Sun S, Chen D, Yuan T, Chen Y, Wang D, et al. Puerarin attenuates the endothelial-mesenchymal transition induced by oxidative stress in human coronary artery endothelial cells through PI3K/AKT pathway. Eur J Pharmacol. 2020;886:173472. 10.1016/j.ejphar.2020.173472.Search in Google Scholar PubMed

[36] Guan BF, Dai XF, Huang QB, Zhao D, Shi JL, Chen C, et al. Icariside II ameliorates myocardial ischemia and reperfusion injury by attenuating inflammation and apoptosis through the regulation of the PI3K/AKT signaling pathway. Mol Med Rep. 2020;22(4):3151–60. 10.3892/mmr.2020.11396.Search in Google Scholar PubMed PubMed Central

[37] Ju X, Sun Y, Zhang F, Wei X, Wang Z, He X. Long non-coding RNA LINC02747 promotes the proliferation of clear cell renal cell carcinoma by inhibiting miR-608 and activating TFE3. Front Oncol. 2020;10:573789. 10.3389/fonc.2020.573789.Search in Google Scholar PubMed PubMed Central

[38] Gu W, Wen D, Lu H, Zhang A, Wang H, Du J, et al. MiR-608 Exerts anti-inflammatory effects by targeting ELANE in monocytes. J Clin Immunol. 2020;40(1):147–57. 10.1007/s10875-019-00702-8.Search in Google Scholar PubMed

[39] Gao S, Liu W, Zhuo X, Wang L, Wang G, Sun T, et al. The activation of mTOR is required for monocyte pro-inflammatory response in patients with coronary artery disease. Clin Sci (Lond). 2015;128(8):517–26. 10.1042/CS20140427.Search in Google Scholar PubMed

[40] Tabaei S, Motallebnezhad M, Tabaee SS. Vitamin D receptor (VDR) gene polymorphisms and risk of coronary artery disease (CAD): systematic review and meta-analysis. Biochem Genet. 2021;59(4):813–36. 10.1007/s10528-021-10038-x.Search in Google Scholar PubMed

[41] Moradi N, Fadaei R, Ahmadi R, Mohammad MH, Shahmohamadnejad S, Tavakoli-Yaraki M, et al. Role of serum MMP-9 levels and vitamin D receptor polymorphisms in the susceptibility to coronary artery disease: an association study in Iranian population. Gene. 2017;628:295–300. 10.1016/j.gene.2017.07.060.Search in Google Scholar PubMed

[42] Bakke D, Sun J. Ancient nuclear receptor VDR with new functions: microbiome and inflammation. Inflamm Bowel Dis. 2018;24(6):1149–54. 10.1093/ibd/izy092.Search in Google Scholar PubMed PubMed Central

[43] Long P, Wang Q, Zhang Y, Zhu X, Yu K, Jiang H, et al. Profile of copper-associated DNA methylation and its association with incident acute coronary syndrome. Clin Epigenetics. 2021;13(1):19. 10.1186/s13148-021-01004-w.Search in Google Scholar PubMed PubMed Central

[44] Gomez-Velazquez M, Badia-Careaga C, Lechuga-Vieco AV, Nieto-Arellano R, Tena JJ, Rollan I, et al. CTCF counter-regulates cardiomyocyte development and maturation programs in the embryonic heart. PLoS Genet. 2017;13(8):e1006985. 10.1371/journal.pgen.1006985.Search in Google Scholar PubMed PubMed Central

[45] Zeng Z, Huang N, Zhang Y, Wang Y, Su Y, Zhang H, et al. CTCF inhibits endoplasmic reticulum stress and apoptosis in cardiomyocytes by upregulating RYR2 via inhibiting S100A1. Life Sci. 2020;242:117158. 10.1016/j.lfs.2019.117158.Search in Google Scholar PubMed

[46] Antony-Debre I, Bluteau D, Itzykson R, Baccini V, Renneville A, Boehlen F, et al. MYH10 protein expression in platelets as a biomarker of RUNX1 and FLI1 alterations. Blood. 2012;120(13):2719–22. 10.1182/blood-2012-04-422352.Search in Google Scholar PubMed

[47] Jiang YX, Yang SW, Li PA, Luo X, Li ZY, Hao YX, et al. The promotion of the transformation of quiescent gastric cancer stem cells by IL-17 and the underlying mechanisms. Oncogene. 2017;36(9):1256–64. 10.1038/onc.2016.291.Search in Google Scholar PubMed PubMed Central

[48] Dounousi E, Duni A, Naka KK, Vartholomatos G, Zoccali C. The innate immune system and cardiovascular disease in ESKD: monocytes and natural killer cells. Curr Vasc Pharmacol. 2021;19(1):63–76. 10.2174/1570161118666200628024027.Search in Google Scholar PubMed

[49] Yang Y, Xu X. Identification of key genes in coronary artery disease: an integrative approach based on weighted gene co-expression network analysis and their correlation with immune infiltration. Aging (Albany NY). 2021;13:8306–19. 10.18632/aging.202638.Search in Google Scholar PubMed PubMed Central

[50] Jabir NR, Firoz CK, Ahmed F, Kamal MA, Hindawi S, Damanhouri GA, et al. Reduction in CD16/CD56 and CD16/CD3/CD56 natural killer cells in coronary artery disease. Immunol Invest. 2017;46(5):526–35. 10.1080/08820139.2017.1306866.Search in Google Scholar PubMed

[51] Yan W, Zhou L, Wen S, Duan Q, Huang F, Tang Y, et al. Differential loss of natural killer cell activity in patients with acute myocardial infarction and stable angina pectoris. Int J Clin Exp Pathol. 2015;8(11):14667–75.Search in Google Scholar

[52] van Duijn J, Kuiper J, Slutter B. The many faces of CD8+ T cells in atherosclerosis. Curr Opin Lipidol. 2018;29(5):411–6. 10.1097/MOL.0000000000000541.Search in Google Scholar PubMed

[53] Ong S, Rose NR, Cihakova D. Natural killer cells in inflammatory heart disease. Clin Immunol. 2017;175:26–33. 10.1016/j.clim.2016.11.010.Search in Google Scholar PubMed PubMed Central

[54] Guttman M, Rinn JL. Modular regulatory principles of large non-coding RNAs. Nature. 2012;482(7385):339–46. 10.1038/nature10887.Search in Google Scholar PubMed PubMed Central

[55] Peng Z, Zhang C, Duan C. Functions and mechanisms of long noncoding RNAs in lung cancer. Onco Targets Ther. 2016;9:4411–24. 10.2147/OTT.S109549.Search in Google Scholar PubMed PubMed Central

[56] Hussen BM, Shoorei H, Mohaqiq M, Dinger ME, Hidayat HJ, Taheri M, et al. The impact of non-coding RNAs in the epithelial to mesenchymal transition. Front Mol Biosci. 2021;8:665199. 10.3389/fmolb.2021.665199.Search in Google Scholar PubMed PubMed Central

[57] Heallen T, Zhang M, Wang J, Bonilla-Claudio M, Klysik E, Johnson RL, et al. Hippo pathway inhibits Wnt signaling to restrain cardiomyocyte proliferation and heart size. Science. 2011;332(6028):458–61. 10.1126/science.1199010.Search in Google Scholar PubMed PubMed Central

[58] Osoegawa K, Iovannisci DM, Lin B, Parodi C, Schultz K, Shaw GM, et al. Identification of novel candidate gene loci and increased sex chromosome aneuploidy among infants with conotruncal heart defects. Am J Med Genet A. 2014;164A(2):397–406. 10.1002/ajmg.a.36291.Search in Google Scholar PubMed PubMed Central

[59] Jackson-Weaver O, Ungvijanpunya N, Yuan Y, Qian J, Gou Y, Wu J, et al. PRMT1-p53 Pathway controls epicardial EMT and invasion. Cell Rep. 2020;31(10):107739. 10.1016/j.celrep.2020.107739.Search in Google Scholar PubMed PubMed Central

[60] Cooley BC, Nevado J, Mellad J, Yang D, St Hilaire C, Negro A, et al. TGF-beta signaling mediates endothelial-to-mesenchymal transition (EndMT) during vein graft remodeling. Sci Transl Med. 2014;6(227):227ra34. 10.1126/scitranslmed.3006927.Search in Google Scholar PubMed PubMed Central

[61] Cheng SL, Shao JS, Behrmann A, Krchma K, Towler DA. Dkk1 and MSX2-Wnt7b signaling reciprocally regulate the endothelial-mesenchymal transition in aortic endothelial cells. Arterioscler Thromb Vasc Biol. 2013;33(7):1679–89. 10.1161/ATVBAHA.113.300647.Search in Google Scholar PubMed PubMed Central

[62] Kato H, Fu YY, Zhu J, Wang L, Aafaqi S, Rahkonen O, et al. Pulmonary vein stenosis and the pathophysiology of “upstream” pulmonary veins. J Thorac Cardiovasc Surg. 2014;148(1):245–53. 10.1016/j.jtcvs.2013.08.046.Search in Google Scholar PubMed

[63] Wu X, Du X, Yang Y, Liu X, Liu X, Zhang N, et al. Inhibition of miR-122 reduced atherosclerotic lesion formation by regulating NPAS3-mediated endothelial to mesenchymal transition. Life Sci. 2021;265:118816. 10.1016/j.lfs.2020.118816.Search in Google Scholar PubMed

[64] Ricciardi M, Zanotto M, Malpeli G, Bassi G, Perbellini O, Chilosi M, et al. Epithelial-to-mesenchymal transition (EMT) induced by inflammatory priming elicits mesenchymal stromal cell-like immune-modulatory properties in cancer cells. Br J Cancer. 2015;112(6):1067–75. 10.1038/bjc.2015.29.Search in Google Scholar PubMed PubMed Central

[65] Pioli PD, Dahlem TJ, Weis JJ, Weis JH. Deletion of Snai2 and Snai3 results in impaired physical development compounded by lymphocyte deficiency. PLoS One. 2013;8(7):e69216. 10.1371/journal.pone.0069216.Search in Google Scholar PubMed PubMed Central

[66] Pioli PD, Whiteside SK, Weis JJ, Weis JH. Snai2 and Snai3 transcriptionally regulate cellular fitness and functionality of T cell lineages through distinct gene programs. Immunobiology. 2016;221(5):618–33. 10.1016/j.imbio.2016.01.007.Search in Google Scholar PubMed PubMed Central

[67] Manduteanu I, Simionescu M. Inflammation in atherosclerosis: a cause or a result of vascular disorders? J Cell Mol Med. 2012;16(9):1978–90. 10.1111/j.1582-4934.2012.01552.x.Search in Google Scholar PubMed PubMed Central

[68] Jabir NR, Tabrez S. Cardiovascular disease management through restrained inflammatory responses. Curr Pharm Des. 2016;22(7):940–6. 10.2174/1381612822666151209153823.Search in Google Scholar PubMed

Received: 2021-09-03
Revised: 2022-02-26
Accepted: 2022-03-22
Published Online: 2022-04-22

© 2022 Xiang Xu et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. AMBRA1 attenuates the proliferation of uveal melanoma cells
  3. A ceRNA network mediated by LINC00475 in papillary thyroid carcinoma
  4. Differences in complications between hepatitis B-related cirrhosis and alcohol-related cirrhosis
  5. Effect of gestational diabetes mellitus on lipid profile: A systematic review and meta-analysis
  6. Long noncoding RNA NR2F1-AS1 stimulates the tumorigenic behavior of non-small cell lung cancer cells by sponging miR-363-3p to increase SOX4
  7. Promising novel biomarkers and candidate small-molecule drugs for lung adenocarcinoma: Evidence from bioinformatics analysis of high-throughput data
  8. Plasmapheresis: Is it a potential alternative treatment for chronic urticaria?
  9. The biomarkers of key miRNAs and gene targets associated with extranodal NK/T-cell lymphoma
  10. Gene signature to predict prognostic survival of hepatocellular carcinoma
  11. Effects of miRNA-199a-5p on cell proliferation and apoptosis of uterine leiomyoma by targeting MED12
  12. Does diabetes affect paraneoplastic thrombocytosis in colorectal cancer?
  13. Is there any effect on imprinted genes H19, PEG3, and SNRPN during AOA?
  14. Leptin and PCSK9 concentrations are associated with vascular endothelial cytokines in patients with stable coronary heart disease
  15. Pericentric inversion of chromosome 6 and male fertility problems
  16. Staple line reinforcement with nebulized cyanoacrylate glue in laparoscopic sleeve gastrectomy: A propensity score-matched study
  17. Retrospective analysis of crescent score in clinical prognosis of IgA nephropathy
  18. Expression of DNM3 is associated with good outcome in colorectal cancer
  19. Activation of SphK2 contributes to adipocyte-induced EOC cell proliferation
  20. CRRT influences PICCO measurements in febrile critically ill patients
  21. SLCO4A1-AS1 mediates pancreatic cancer development via miR-4673/KIF21B axis
  22. lncRNA ACTA2-AS1 inhibits malignant phenotypes of gastric cancer cells
  23. circ_AKT3 knockdown suppresses cisplatin resistance in gastric cancer
  24. Prognostic value of nicotinamide N-methyltransferase in human cancers: Evidence from a meta-analysis and database validation
  25. GPC2 deficiency inhibits cell growth and metastasis in colon adenocarcinoma
  26. A pan-cancer analysis of the oncogenic role of Holliday junction recognition protein in human tumors
  27. Radiation increases COL1A1, COL3A1, and COL1A2 expression in breast cancer
  28. Association between preventable risk factors and metabolic syndrome
  29. miR-29c-5p knockdown reduces inflammation and blood–brain barrier disruption by upregulating LRP6
  30. Cardiac contractility modulation ameliorates myocardial metabolic remodeling in a rabbit model of chronic heart failure through activation of AMPK and PPAR-α pathway
  31. Quercitrin protects human bronchial epithelial cells from oxidative damage
  32. Smurf2 suppresses the metastasis of hepatocellular carcinoma via ubiquitin degradation of Smad2
  33. circRNA_0001679/miR-338-3p/DUSP16 axis aggravates acute lung injury
  34. Sonoclot’s usefulness in prediction of cardiopulmonary arrest prognosis: A proof of concept study
  35. Four drug metabolism-related subgroups of pancreatic adenocarcinoma in prognosis, immune infiltration, and gene mutation
  36. Decreased expression of miR-195 mediated by hypermethylation promotes osteosarcoma
  37. LMO3 promotes proliferation and metastasis of papillary thyroid carcinoma cells by regulating LIMK1-mediated cofilin and the β-catenin pathway
  38. Cx43 upregulation in HUVECs under stretch via TGF-β1 and cytoskeletal network
  39. Evaluation of menstrual irregularities after COVID-19 vaccination: Results of the MECOVAC survey
  40. Histopathologic findings on removed stomach after sleeve gastrectomy. Do they influence the outcome?
  41. Analysis of the expression and prognostic value of MT1-MMP, β1-integrin and YAP1 in glioma
  42. Optimal diagnosis of the skin cancer using a hybrid deep neural network and grasshopper optimization algorithm
  43. miR-223-3p alleviates TGF-β-induced epithelial-mesenchymal transition and extracellular matrix deposition by targeting SP3 in endometrial epithelial cells
  44. Clinical value of SIRT1 as a prognostic biomarker in esophageal squamous cell carcinoma, a systematic meta-analysis
  45. circ_0020123 promotes cell proliferation and migration in lung adenocarcinoma via PDZD8
  46. miR-22-5p regulates the self-renewal of spermatogonial stem cells by targeting EZH2
  47. hsa-miR-340-5p inhibits epithelial–mesenchymal transition in endometriosis by targeting MAP3K2 and inactivating MAPK/ERK signaling
  48. circ_0085296 inhibits the biological functions of trophoblast cells to promote the progression of preeclampsia via the miR-942-5p/THBS2 network
  49. TCD hemodynamics findings in the subacute phase of anterior circulation stroke patients treated with mechanical thrombectomy
  50. Development of a risk-stratification scoring system for predicting risk of breast cancer based on non-alcoholic fatty liver disease, non-alcoholic fatty pancreas disease, and uric acid
  51. Tollip promotes hepatocellular carcinoma progression via PI3K/AKT pathway
  52. circ_0062491 alleviates periodontitis via the miR-142-5p/IGF1 axis
  53. Human amniotic fluid as a source of stem cells
  54. lncRNA NONRATT013819.2 promotes transforming growth factor-β1-induced myofibroblastic transition of hepatic stellate cells by miR24-3p/lox
  55. NORAD modulates miR-30c-5p-LDHA to protect lung endothelial cells damage
  56. Idiopathic pulmonary fibrosis telemedicine management during COVID-19 outbreak
  57. Risk factors for adverse drug reactions associated with clopidogrel therapy
  58. Serum zinc associated with immunity and inflammatory markers in Covid-19
  59. The relationship between night shift work and breast cancer incidence: A systematic review and meta-analysis of observational studies
  60. LncRNA expression in idiopathic achalasia: New insight and preliminary exploration into pathogenesis
  61. Notoginsenoside R1 alleviates spinal cord injury through the miR-301a/KLF7 axis to activate Wnt/β-catenin pathway
  62. Moscatilin suppresses the inflammation from macrophages and T cells
  63. Zoledronate promotes ECM degradation and apoptosis via Wnt/β-catenin
  64. Epithelial-mesenchymal transition-related genes in coronary artery disease
  65. The effect evaluation of traditional vaginal surgery and transvaginal mesh surgery for severe pelvic organ prolapse: 5 years follow-up
  66. Repeated partial splenic artery embolization for hypersplenism improves platelet count
  67. Low expression of miR-27b in serum exosomes of non-small cell lung cancer facilitates its progression by affecting EGFR
  68. Exosomal hsa_circ_0000519 modulates the NSCLC cell growth and metastasis via miR-1258/RHOV axis
  69. miR-455-5p enhances 5-fluorouracil sensitivity in colorectal cancer cells by targeting PIK3R1 and DEPDC1
  70. The effect of tranexamic acid on the reduction of intraoperative and postoperative blood loss and thromboembolic risk in patients with hip fracture
  71. Isocitrate dehydrogenase 1 mutation in cholangiocarcinoma impairs tumor progression by sensitizing cells to ferroptosis
  72. Artemisinin protects against cerebral ischemia and reperfusion injury via inhibiting the NF-κB pathway
  73. A 16-gene signature associated with homologous recombination deficiency for prognosis prediction in patients with triple-negative breast cancer
  74. Lidocaine ameliorates chronic constriction injury-induced neuropathic pain through regulating M1/M2 microglia polarization
  75. MicroRNA 322-5p reduced neuronal inflammation via the TLR4/TRAF6/NF-κB axis in a rat epilepsy model
  76. miR-1273h-5p suppresses CXCL12 expression and inhibits gastric cancer cell invasion and metastasis
  77. Clinical characteristics of pneumonia patients of long course of illness infected with SARS-CoV-2
  78. circRNF20 aggravates the malignancy of retinoblastoma depending on the regulation of miR-132-3p/PAX6 axis
  79. Linezolid for resistant Gram-positive bacterial infections in children under 12 years: A meta-analysis
  80. Rack1 regulates pro-inflammatory cytokines by NF-κB in diabetic nephropathy
  81. Comprehensive analysis of molecular mechanism and a novel prognostic signature based on small nuclear RNA biomarkers in gastric cancer patients
  82. Smog and risk of maternal and fetal birth outcomes: A retrospective study in Baoding, China
  83. Let-7i-3p inhibits the cell cycle, proliferation, invasion, and migration of colorectal cancer cells via downregulating CCND1
  84. β2-Adrenergic receptor expression in subchondral bone of patients with varus knee osteoarthritis
  85. Possible impact of COVID-19 pandemic and lockdown on suicide behavior among patients in Southeast Serbia
  86. In vitro antimicrobial activity of ozonated oil in liposome eyedrop against multidrug-resistant bacteria
  87. Potential biomarkers for inflammatory response in acute lung injury
  88. A low serum uric acid concentration predicts a poor prognosis in adult patients with candidemia
  89. Antitumor activity of recombinant oncolytic vaccinia virus with human IL2
  90. ALKBH5 inhibits TNF-α-induced apoptosis of HUVECs through Bcl-2 pathway
  91. Risk prediction of cardiovascular disease using machine learning classifiers
  92. Value of ultrasonography parameters in diagnosing polycystic ovary syndrome
  93. Bioinformatics analysis reveals three key genes and four survival genes associated with youth-onset NSCLC
  94. Identification of autophagy-related biomarkers in patients with pulmonary arterial hypertension based on bioinformatics analysis
  95. Protective effects of glaucocalyxin A on the airway of asthmatic mice
  96. Overexpression of miR-100-5p inhibits papillary thyroid cancer progression via targeting FZD8
  97. Bioinformatics-based analysis of SUMOylation-related genes in hepatocellular carcinoma reveals a role of upregulated SAE1 in promoting cell proliferation
  98. Effectiveness and clinical benefits of new anti-diabetic drugs: A real life experience
  99. Identification of osteoporosis based on gene biomarkers using support vector machine
  100. Tanshinone IIA reverses oxaliplatin resistance in colorectal cancer through microRNA-30b-5p/AVEN axis
  101. miR-212-5p inhibits nasopharyngeal carcinoma progression by targeting METTL3
  102. Association of ST-T changes with all-cause mortality among patients with peripheral T-cell lymphomas
  103. LINC00665/miRNAs axis-mediated collagen type XI alpha 1 correlates with immune infiltration and malignant phenotypes in lung adenocarcinoma
  104. The perinatal factors that influence the excretion of fecal calprotectin in premature-born children
  105. Effect of femoral head necrosis cystic area on femoral head collapse and stress distribution in femoral head: A clinical and finite element study
  106. Does the use of 3D-printed cones give a chance to postpone the use of megaprostheses in patients with large bone defects in the knee joint?
  107. lncRNA HAGLR modulates myocardial ischemia–reperfusion injury in mice through regulating miR-133a-3p/MAPK1 axis
  108. Protective effect of ghrelin on intestinal I/R injury in rats
  109. In vivo knee kinematics of an innovative prosthesis design
  110. Relationship between the height of fibular head and the incidence and severity of knee osteoarthritis
  111. lncRNA WT1-AS attenuates hypoxia/ischemia-induced neuronal injury during cerebral ischemic stroke via miR-186-5p/XIAP axis
  112. Correlation of cardiac troponin T and APACHE III score with all-cause in-hospital mortality in critically ill patients with acute pulmonary embolism
  113. LncRNA LINC01857 reduces metastasis and angiogenesis in breast cancer cells via regulating miR-2052/CENPQ axis
  114. Endothelial cell-specific molecule 1 (ESM1) promoted by transcription factor SPI1 acts as an oncogene to modulate the malignant phenotype of endometrial cancer
  115. SELENBP1 inhibits progression of colorectal cancer by suppressing epithelial–mesenchymal transition
  116. Visfatin is negatively associated with coronary artery lesions in subjects with impaired fasting glucose
  117. Treatment and outcomes of mechanical complications of acute myocardial infarction during the Covid-19 era: A comparison with the pre-Covid-19 period. A systematic review and meta-analysis
  118. Neonatal stroke surveillance study protocol in the United Kingdom and Republic of Ireland
  119. Oncogenic role of TWF2 in human tumors: A pan-cancer analysis
  120. Mean corpuscular hemoglobin predicts the length of hospital stay independent of severity classification in patients with acute pancreatitis
  121. Association of gallstone and polymorphisms of UGT1A1*27 and UGT1A1*28 in patients with hepatitis B virus-related liver failure
  122. TGF-β1 upregulates Sar1a expression and induces procollagen-I secretion in hypertrophic scarring fibroblasts
  123. Antisense lncRNA PCNA-AS1 promotes esophageal squamous cell carcinoma progression through the miR-2467-3p/PCNA axis
  124. NK-cell dysfunction of acute myeloid leukemia in relation to the renin–angiotensin system and neurotransmitter genes
  125. The effect of dilution with glucose and prolonged injection time on dexamethasone-induced perineal irritation – A randomized controlled trial
  126. miR-146-5p restrains calcification of vascular smooth muscle cells by suppressing TRAF6
  127. Role of lncRNA MIAT/miR-361-3p/CCAR2 in prostate cancer cells
  128. lncRNA NORAD promotes lung cancer progression by competitively binding to miR-28-3p with E2F2
  129. Noninvasive diagnosis of AIH/PBC overlap syndrome based on prediction models
  130. lncRNA FAM230B is highly expressed in colorectal cancer and suppresses the maturation of miR-1182 to increase cell proliferation
  131. circ-LIMK1 regulates cisplatin resistance in lung adenocarcinoma by targeting miR-512-5p/HMGA1 axis
  132. LncRNA SNHG3 promoted cell proliferation, migration, and metastasis of esophageal squamous cell carcinoma via regulating miR-151a-3p/PFN2 axis
  133. Risk perception and affective state on work exhaustion in obstetrics during the COVID-19 pandemic
  134. lncRNA-AC130710/miR-129-5p/mGluR1 axis promote migration and invasion by activating PKCα-MAPK signal pathway in melanoma
  135. SNRPB promotes cell cycle progression in thyroid carcinoma via inhibiting p53
  136. Xylooligosaccharides and aerobic training regulate metabolism and behavior in rats with streptozotocin-induced type 1 diabetes
  137. Serpin family A member 1 is an oncogene in glioma and its translation is enhanced by NAD(P)H quinone dehydrogenase 1 through RNA-binding activity
  138. Silencing of CPSF7 inhibits the proliferation, migration, and invasion of lung adenocarcinoma cells by blocking the AKT/mTOR signaling pathway
  139. Ultrasound-guided lumbar plexus block versus transversus abdominis plane block for analgesia in children with hip dislocation: A double-blind, randomized trial
  140. Relationship of plasma MBP and 8-oxo-dG with brain damage in preterm
  141. Identification of a novel necroptosis-associated miRNA signature for predicting the prognosis in head and neck squamous cell carcinoma
  142. Delayed femoral vein ligation reduces operative time and blood loss during hip disarticulation in patients with extremity tumors
  143. The expression of ASAP3 and NOTCH3 and the clinicopathological characteristics of adult glioma patients
  144. Longitudinal analysis of factors related to Helicobacter pylori infection in Chinese adults
  145. HOXA10 enhances cell proliferation and suppresses apoptosis in esophageal cancer via activating p38/ERK signaling pathway
  146. Meta-analysis of early-life antibiotic use and allergic rhinitis
  147. Marital status and its correlation with age, race, and gender in prognosis of tonsil squamous cell carcinomas
  148. HPV16 E6E7 up-regulates KIF2A expression by activating JNK/c-Jun signal, is beneficial to migration and invasion of cervical cancer cells
  149. Amino acid profiles in the tissue and serum of patients with liver cancer
  150. Pain in critically ill COVID-19 patients: An Italian retrospective study
  151. Immunohistochemical distribution of Bcl-2 and p53 apoptotic markers in acetamiprid-induced nephrotoxicity
  152. Estradiol pretreatment in GnRH antagonist protocol for IVF/ICSI treatment
  153. Long non-coding RNAs LINC00689 inhibits the apoptosis of human nucleus pulposus cells via miR-3127-5p/ATG7 axis-mediated autophagy
  154. The relationship between oxygen therapy, drug therapy, and COVID-19 mortality
  155. Monitoring hypertensive disorders in pregnancy to prevent preeclampsia in pregnant women of advanced maternal age: Trial mimicking with retrospective data
  156. SETD1A promotes the proliferation and glycolysis of nasopharyngeal carcinoma cells by activating the PI3K/Akt pathway
  157. The role of Shunaoxin pills in the treatment of chronic cerebral hypoperfusion and its main pharmacodynamic components
  158. TET3 governs malignant behaviors and unfavorable prognosis of esophageal squamous cell carcinoma by activating the PI3K/AKT/GSK3β/β-catenin pathway
  159. Associations between morphokinetic parameters of temporary-arrest embryos and the clinical prognosis in FET cycles
  160. Long noncoding RNA WT1-AS regulates trophoblast proliferation, migration, and invasion via the microRNA-186-5p/CADM2 axis
  161. The incidence of bronchiectasis in chronic obstructive pulmonary disease
  162. Integrated bioinformatics analysis shows integrin alpha 3 is a prognostic biomarker for pancreatic cancer
  163. Inhibition of miR-21 improves pulmonary vascular responses in bronchopulmonary dysplasia by targeting the DDAH1/ADMA/NO pathway
  164. Comparison of hospitalized patients with severe pneumonia caused by COVID-19 and influenza A (H7N9 and H1N1): A retrospective study from a designated hospital
  165. lncRNA ZFAS1 promotes intervertebral disc degeneration by upregulating AAK1
  166. Pathological characteristics of liver injury induced by N,N-dimethylformamide: From humans to animal models
  167. lncRNA ELFN1-AS1 enhances the progression of colon cancer by targeting miR-4270 to upregulate AURKB
  168. DARS-AS1 modulates cell proliferation and migration of gastric cancer cells by regulating miR-330-3p/NAT10 axis
  169. Dezocine inhibits cell proliferation, migration, and invasion by targeting CRABP2 in ovarian cancer
  170. MGST1 alleviates the oxidative stress of trophoblast cells induced by hypoxia/reoxygenation and promotes cell proliferation, migration, and invasion by activating the PI3K/AKT/mTOR pathway
  171. Bifidobacterium lactis Probio-M8 ameliorated the symptoms of type 2 diabetes mellitus mice by changing ileum FXR-CYP7A1
  172. circRNA DENND1B inhibits tumorigenicity of clear cell renal cell carcinoma via miR-122-5p/TIMP2 axis
  173. EphA3 targeted by miR-3666 contributes to melanoma malignancy via activating ERK1/2 and p38 MAPK pathways
  174. Pacemakers and methylprednisolone pulse therapy in immune-related myocarditis concomitant with complete heart block
  175. miRNA-130a-3p targets sphingosine-1-phosphate receptor 1 to activate the microglial and astrocytes and to promote neural injury under the high glucose condition
  176. Review Articles
  177. Current management of cancer pain in Italy: Expert opinion paper
  178. Hearing loss and brain disorders: A review of multiple pathologies
  179. The rationale for using low-molecular weight heparin in the therapy of symptomatic COVID-19 patients
  180. Amyotrophic lateral sclerosis and delayed onset muscle soreness in light of the impaired blink and stretch reflexes – watch out for Piezo2
  181. Interleukin-35 in autoimmune dermatoses: Current concepts
  182. Recent discoveries in microbiota dysbiosis, cholangiocytic factors, and models for studying the pathogenesis of primary sclerosing cholangitis
  183. Advantages of ketamine in pediatric anesthesia
  184. Congenital adrenal hyperplasia. Role of dentist in early diagnosis
  185. Migraine management: Non-pharmacological points for patients and health care professionals
  186. Atherogenic index of plasma and coronary artery disease: A systematic review
  187. Physiological and modulatory role of thioredoxins in the cellular function
  188. Case Reports
  189. Intrauterine Bakri balloon tamponade plus cervical cerclage for the prevention and treatment of postpartum haemorrhage in late pregnancy complicated with acute aortic dissection: Case series
  190. A case of successful pembrolizumab monotherapy in a patient with advanced lung adenocarcinoma: Use of multiple biomarkers in combination for clinical practice
  191. Unusual neurological manifestations of bilateral medial medullary infarction: A case report
  192. Atypical symptoms of malignant hyperthermia: A rare causative mutation in the RYR1 gene
  193. A case report of dermatomyositis with the missed diagnosis of non-small cell lung cancer and concurrence of pulmonary tuberculosis
  194. A rare case of endometrial polyp complicated with uterine inversion: A case report and clinical management
  195. Spontaneous rupturing of splenic artery aneurysm: Another reason for fatal syncope and shock (Case report and literature review)
  196. Fungal infection mimicking COVID-19 infection – A case report
  197. Concurrent aspergillosis and cystic pulmonary metastases in a patient with tongue squamous cell carcinoma
  198. Paraganglioma-induced inverted takotsubo-like cardiomyopathy leading to cardiogenic shock successfully treated with extracorporeal membrane oxygenation
  199. Lineage switch from lymphoma to myeloid neoplasms: First case series from a single institution
  200. Trismus during tracheal extubation as a complication of general anaesthesia – A case report
  201. Simultaneous treatment of a pubovesical fistula and lymph node metastasis secondary to multimodal treatment for prostate cancer: Case report and review of the literature
  202. Two case reports of skin vasculitis following the COVID-19 immunization
  203. Ureteroiliac fistula after oncological surgery: Case report and review of the literature
  204. Synchronous triple primary malignant tumours in the bladder, prostate, and lung harbouring TP53 and MEK1 mutations accompanied with severe cardiovascular diseases: A case report
  205. Huge mucinous cystic neoplasms with adhesion to the left colon: A case report and literature review
  206. Commentary
  207. Commentary on “Clinicopathological features of programmed cell death-ligand 1 expression in patients with oral squamous cell carcinoma”
  208. Rapid Communication
  209. COVID-19 fear, post-traumatic stress, growth, and the role of resilience
  210. Erratum
  211. Erratum to “Tollip promotes hepatocellular carcinoma progression via PI3K/AKT pathway”
  212. Erratum to “Effect of femoral head necrosis cystic area on femoral head collapse and stress distribution in femoral head: A clinical and finite element study”
  213. Erratum to “lncRNA NORAD promotes lung cancer progression by competitively binding to miR-28-3p with E2F2”
  214. Retraction
  215. Expression and role of ABIN1 in sepsis: In vitro and in vivo studies
  216. Retraction to “miR-519d downregulates LEP expression to inhibit preeclampsia development”
  217. Special Issue Computational Intelligence Methodologies Meets Recurrent Cancers - Part II
  218. Usefulness of close surveillance for rectal cancer patients after neoadjuvant chemoradiotherapy
Downloaded on 17.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/med-2022-0476/html
Scroll to top button