Next Article in Journal
Length of Hospitalization-Related Differences and Associated Long-Term Prognosis of Patients with Cardiac Resynchronization Therapy: A Propensity Score-Matched Cohort
Previous Article in Journal
Minimally Invasive Surgical Repair of a Partial Atrioventricular Canal Defect in a 20-Year-Old Patient—A Case Report and Review of Literature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bioinformatics Analysis of Common Differential Genes of Viral Myocarditis and Dilated Cardiomyopathy: Screening for Potential Pharmacological Compounds

Department of Cardiology, Dushu Lake Hospital Affiliated to Soochow University, Suzhou 215000, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Cardiovasc. Dev. Dis. 2022, 9(10), 353; https://doi.org/10.3390/jcdd9100353
Submission received: 14 September 2022 / Revised: 9 October 2022 / Accepted: 11 October 2022 / Published: 15 October 2022

Abstract

:
(1) Background: The mechanism of viral myocarditis (VMC) progression to dilated cardiomyopathy (DCM) remains unclear. The aim of this study was to identify key genes in the progression of VMC to DCM, so as to find potential therapeutic drugs and provide insights for future research. (2) Methods: Differential expression analysis of GSE4172 and GSE17800 from the Gene Expression Omnibus (GEO) database was performed using GEO2R, which contained genome-wide analysis of myocardial biopsies from VMC and DCM, respectively. We used the Venn diagram analysis to screen the common differentially expressed genes (DEGs). GO functional enrichment analysis and KEGG pathway analysis were also performed. Then we conducted protein–protein interaction (PPI) networks using STRING and identified hub genes using Cytoscape. Finally, we used cMAP to screen out candidate compounds targeting these hub genes; (3) Results: In total, 2143 DEGs for VMC and 1365 DEGs for DCM were found. Then a total of 191 common DEGs were identified. Biological processes and pathway involved in these genes mainly include GABA−gated chloride ion channel activity and Rap1 signaling pathway. A total of 14 hub genes were identified. PPI network showed these hubs mainly enriched in regulation of WNT signaling pathway and GABA-gated chloride ion channel activity. Subgroup analysis of Severe VMC cohort revealed 10 hub genes which mainly clustered in GABA channel activity, extracellular matrix remodeling and sarcomere dysfunction. Using cMAP, we obtained top 10 potential medications, but only amlodipine is currently viable; (4) Conclusions: Our study finds the hub genes and reveals the important role of GABA-gated chloride ion channel, Rap1 signaling pathway, WNT signaling pathway, extracellular matrix remodeling and sarcomere dysfunction in the progression from VMC to DCM. Amlodipine is a potential viable drug in preventing the progression of VMC to DCM.

1. Introduction

Viral myocarditis (VMC) is a localized or diffuse inflammatory disease of the myocardium associated with viral infection, which may have an acute or chronic course [1]. Most patients have a good prognosis, but 10%–15% may progress to dilated cardiomyopathy (DCM) [2]. DCM is a myocardial disease characterized by unilateral or bilateral enlargement of heart cavities and decreased myocardial contractility. Its clinical manifestations are mainly progressive congestive heart failure, often accompanied by various types of arrhythmias, which has a serious impact on patients’ quality of life and life expectancy [3]. Currently, there is a lack of effective clinical treatment. Patients with dilated cardiomyopathy have only a 55% 5-year survival rate with current heart failure therapy, suggesting an urgent need for targeted therapy strategies [4].
The viruses known to cause myocarditis include coxsackie virus (A,B), echovirus, adenovirus, influenza virus (A,B), cytomegalovirus, and chiropoliovirus [5]. The hypothesis of the progression from VMC to DCM has been proposed for a long time. Kawai postmortem examined 44 DCM patients, 8 of whom had definite histological myocarditis [6]. It was found that neutralizing antibody titers of Coxsackie GROUP B (CB) virus were greater than 1:1024 in approximately 30% of DCM patients [7]. Of 165 serologically confirmed CB viral myocarditis patients with long-term follow-up, 13% developed DCM [8]. Due to the close relationship between VMC and DCM, it is often impossible to separate them completely in clinical practice, so some scholars believe that they are different stages of the same pathological process. However, the actual pathogenesis is not well understood. In recent years, with the great improvement of molecular biology techniques such as nucleic acid hybridization, polymerase chain reaction (PCR), endocardial biopsy, and modern immunohistochemistry, more evidence has been obtained on the mechanism of progression from VMC to DCM. Virus infection, immune response, and apoptosis are considered to be important factors [9,10,11]. However, the key genes in the pathogenesis from VMC to DCM have not been identified so far.
In this study, through bioinformatics, we searched for the key genes in the transition from viral myocarditis to non-familial dilated cardiomyopathy by screening the common differentially expressed genes (DEGs) between VMC and DCM, with the hope of finding potential drugs that could prevent the progression of VMC to DCM to some extent and providing clues for further research.

2. Materials and Methods

2.1. Acquisition of Data of Gene Expression Profiles

The datasets of gene expression profiles GSE4172 and GSE17800 were acquired by logging in the Gene Expression Omnibus (GEO) database (accessed on 5 June 2022, https://www.ncbi.nlm.nih.gov/gds/) [12,13]. In GSE4172, there were 8 samples from VMC patients infected with parvovirus B19 and 4 samples from controls. In GSE17800, the number of samples from DCM patients and controls was 40 and 8, respectively. Both of the two chips were from endomyocardial biopsies (EMBs), the detailed information of which is shown in Table 1.

2.2. Processing of Raw Data and Screening of Differentially Expressed Genes (DEGs)

We used fold change and corrected p-values to draw Volcano plots. Boxplots were drawn using the R package ggplot2; Expression heatmaps were drawn using the R package heatmap. DEGs were analyzed using the GEO2R tool. We defined p < 0.05 and |log2 fold change (FC)| > 0.58 (namely fold change >1.5 times) as statistically significant DEGs.

2.3. Screening of Common Differentially Expressed Genes

The differential genes selected from the two chips underwent Venn diagram analysis and a minority of DEGs shared by both were defined as common DEGs, which will be used for subsequent analysis.

2.4. Gene Ontology (GO) Enrichment Analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Analysis

GO functional enrichment analysis which included molecular function (MF), biological process (BP), and cellular component [CC] analysis and KEGG pathway analysis of the DEGs were performed using the R package clusterProfiler.

2.5. Construction of Protein–Protein Interaction (PPI) Networks and Screening of Hub Genes

PPI network analysis was performed for the common DEGs of VMC and DCM by using String (accessed on 6 June 2022, https://www.string-db.org/) to predict the interaction between proteins encoded by DEGs. Confidence interaction score was set at 0.15 for the significant criterion. According to the result of PPI network, we used Cytoscape and the MCODE module to screen the hub genes.

2.6. Screening of Potential Pharmacological Targets

The ConnectivityMap (cMAP) (accessed on 6 June 2022, https://clue.io) database contains data on changes in gene expression profiles caused by more than 30,000 small-molecule compounds acting on multiple cell lines and can be used to compare similarities between drug-induced gene profiles and gene expressions, with connectivity scores ranging from −100 to 100. A positive score suggests that compounds can cause changes similar to those in uploaded genes; a negative score suggests that the compound causes changes opposite to those in the uploaded genes. Therefore, for those up-regulated genes, a negative score suggests that the compound may have a therapeutic effect on the disease. On the contrary, for those down-regulated genes, a positive score suggests that the agent may have a therapeutic effect on the disease. Agents with absolute value of the connectivity scores >80 were considered promising predictors.

2.7. Subgroup Analysis

We additionally selected patients with severe myocarditis (sVMC) from GSE4172 for subgroup analysis. We defined sVMC as an EF < 40%. Four patients met the definition of sVMC.

3. Results

3.1. Screening of Differentially Expressed Genes

3.1.1. Screening of Differentially Expressed Genes in GSE417

A total of 2143 DEGs were selected from GSE4172, including 1217 up-regulated genes and 926 downregulated genes. A part of the major DEGs is presented in Table 2. Clustering analysis of these DEGs was performed, as shown in the volcano plot (Figure 1A). Data normalization and cross-comparability were also evaluated. As shown in Figure 1B, the selected samples were centered and the numerical distribution met the standard, indicating that the microarray data were of high quality and had cross-comparability. Based on the p-value, the former 50 DEGs with the lowest p-values were analyzed in the heatmap (Figure 1C).

3.1.2. Screening of Differentially Expressed Genes in GSE17800

Totally, 1365 DEGs were selected from GSE17800, of which 869 were upregulated genes and 496 were downregulated genes. A part of major DEGs is presented in Table 3. Similarly, clustering analysis, data normalization and cross comparability and heatmap of these DEGs were also performed. As a result, the volcano plot is shown in Figure 2A. The box plot showing a concentrated sample distribution is presented in Figure 2B. The heatmap of the former 50 DEGs with the lowest p-values was shown in Figure 2C.

3.1.3. Screening of Common Differentially Expressed Genes

Using Venn analysis, a total of 191 common DEGs were screened from GSE4172 and GSE17800 (Figure 3A). Among these genes, 66 were common up-regulated DEGs (Figure 3B) and 70 were common down-regulated DEGs (Figure 3C). The former ten common DEGs with the highest average |log2FC| are shown in Table 4.

3.2. GO Enrichment and KEGG Pathway Analysis

The gene ontology enrichment analysis of the common up-regulated DEGs of VMC and DCM were mainly enriched in biological process (BP), cellular localization (CL), and molecular function (MF), especially in “negative regulation of anion transport” (GO-ID:1903792, p = 4.99 × 10−5), “GABA receptor complex” (GO-ID:1902710, p = 0.001), “GABA−gated chloride ion channel activity” (GO-ID:0022851, p = 5.51 × 10−4) and “xenobiotic transmembrane transporter activity” (GO-ID:0042910, p = 5.51 × 10−4). KEGG showed that these genes were mainly involved in GABAergic synapse (Table 5, Figure 4A). These results suggest that GABA ion channel transport function and xenobiologic transmembrane transport function are enhanced in patients with VMC and DCM compared with normal controls.
The gene ontology enrichment analysis of the common down-regulated DEGs of VMC and DCM were mainly enriched in biological process (BP) and molecular function (MF), especially in “ventricular cardiac muscle tissue morphogenesis” (GO-ID:0055010, p = 1.96 × 10−4) and “ventricular cardiac muscle tissue development” (GO-ID:0003229, p = 2.94 × 10−4). KEGG enrichment results showed that these genes were mainly involved in Rap1 signaling pathway (Table 6, Figure 4B). These results suggest that ventricular myocardial tissue development and Rap1 signaling pathway are inhibited in patients with VMC and DCM.

3.3. PPI Network Analysis and Screening of Hub Genes

To further explore the roles of these common DEGs, we used protein–protein interaction data from STRING to construct a network which contained 109 nodes and 306 edges (Figure 5). By importing the results of STRING into Cytoscape, we identified 14 hub genes: FGFR2, AQP4, LTBP2, BMP6, WISP1, SFRP4, FRZB, ABCG2, HBEGF, UCHL1, VSNL1, GABRB1, GABRA4, and RGS4 (Figure 6), of which FGFR2 and AQP4 were down-regulated DEGs and others were up-regulated DEGs. The proteins in the network were mainly enriched in regulation of WNT signaling pathway (GO-ID:0030111, p = 0.0413), cellular response to stimulus (GO-ID:0051716, p = 0.0262), and GABA-gated chloride ion channel activity (GO-ID:0022851, p = 0.0467).

3.4. Screening for Potential Pharmacological Targets

Potential targeted drugs of these hub genes were searched in the cMAP database, and the drugs were sorted and screened according to their connectivity score. The top 10 compounds suggested for these common DEGs were HG-6-64-01, selamectin, tetrindole, PD-198306, elvitegravir, amsacrine, GSK-461364, gefitinib, amlodipine, and AZ-628 (Table 7).

3.5. Subgroup Analysis

From the GSE4172 dataset, four VMC patients with EF < 40% formed the sVMC subgroup. Their differential gene analysis with four normal controls revealed a total of 2121 DEGs. There were 180 common DEGs for sVMC and DCM (Figure 7A). GO/KEGG analysis showed that these genes were mainly enriched in the collagen-containing extracellular matrix (Figure 7B). A PPI network containing 143 nodes and 447 edges were constructed. By importing these data into Cytoscape, we identified 10 hub genes (4 upregulated and 6 downregulated): VSNL1, GABRA4, GABRB1, ACTG2, SCN1A, NFASC, AQP4, MYH6, COL1A1 and CSRP2, respectively (Figure 7C). The proteins encoded by these genes were mainly clustered in the GABA-A receptor complex (GO:1902711) and myosin filaments (GO:0032982).

4. Discussion

Dilated cardiomyopathy is a heterogeneous group of diseases that can be divided into two categories: primary (familial inheritance) and secondary. Viral myocarditis has traditionally been considered the most common cause of secondary dilated cardiomyopathy [14,15]. However, the pathogenesis of transition from viral myocarditis to dilated cardiomyopathy is poorly understood, and the key genes and pathways among this pathophysiological process await further discovery. In this study, we first obtained the common DEGs between VMC and DCM based on gene chips from myocardial biopsy. A total of 191 DEGs of VMC and DCM were unearthed.
When compared with the studies that originally generated the reevaluated datasets, the original study that generated the GSE4172 dataset [12] found that cysteine-rich angiogenic inducer 61 (CYR61) and adiponectin (APN) played an important role in inflammatory cardiomyopathy caused by parvovirus B19-infected viral myocarditis. However, in our study, the APN gene was not among the DEGs. Although CYR61 was among the DEGs, it was not found in the final common DEGs, indicating that although CYR61 and APN may play an important role in viral myocarditis, they are not the key genes in the progression of viral myocarditis to dilated cardiomyopathy. The original study that generated the GSE4172 dataset [13] compared the differential genes of responders and non-responders to immunoadsorption with subsequent immunoglobulin G substitution (IA/IgG) treatment in 40 DCM patients with normal controls. Because our study overall compared genetic differences between DCM patients and normal controls, the results of the DEGs differed from this study.
GO enrichment analysis and KEGG pathway analysis showed that the common upregulated GEGs are mainly involved in GABA receptor complex, GABA−gated chloride ion channel activity, xenobiotic transmembrane transporter activity, and GABAergic synapse, indicating that GABA-related pathway functions are enhanced in VMC and DCM patients. This is a novel finding. Currently, it is known that γ-aminobutyric acid (GABA) is widely distributed in the cardiovascular system activity regulatory region of the central nervous system (CNS) and directly regulates cardiac sympathetic nerve activity (C-SNA) [16,17]. On the other hand, there is strong evidence that cardiac sympathetic hyperactivation is associated with DCM [18,19,20]. Therefore, whether GABA dysregulation in VMC leads to the development of DCM through this pathophysiological mechanism is worth exploring.
For the common downregulated DEGs, GO enrichment analysis found these genes were mainly enriched in ventricular cardiac muscle tissue morphogenesis and ventricular cardiac muscle tissue development, which means that ventricular myocardial tissue development are suppressed in patients with VMC and DCM. The KEGG pathway analysis found that these downregulated DEGs were mainly enriched in Rap1 signaling pathway. Rap1 (Ras-proximate-1) is a small GTPase protein which plays a crucial role in mediating cAMP signaling in isolated cardiac tissues and cell lines [21]. Similar to our study, a recent study in found a significant enrichment in Rap1 signaling pathway in VMC mice [22]. Further, basic studies show that suppression of Rap1 impairs cardiac myofibrils and conduction system in zebrafish [23]. Our study suggests that inhibition of the Rap1 signaling pathway in VMC may contribute to the progression of DCM.
Fourteen genes were identified as hub genes in the progression from VMC to DCM. Of them, FGFR2 (Fibroblast growth factor receptor 2) is a tyrosine-protein kinase that plays an essential role in the regulation of cell proliferation, differentiation, migration and apoptosis [24]. AQP4 (Aquaporin-4) forms a water-specific channel and plays an important role in brain water homeostasis and in glymphatic solute transport [25]. LTBP2 (Latent-transforming growth factor beta-binding protein 2) may play an integral structural role in elastic-fiber architectural organization and/or assembly [26]. BMP6 (Bone morphogenetic protein 6) plays an important role in bone formation and iron metabolism [27]. WISP1 (WNT1 inducible signaling pathway protein 1) is a downstream regulator in the Wnt/Frizzled-signaling pathway and is associated with cell survival [28]. SFRP4 (Secreted frizzled-related protein 4) and FRZB (frizzled related protein) are both soluble frizzled-related proteins (sFRPS) that function as modulators of Wnt signaling through direct interaction with Wnts. [29,30]. ABCG2 (ATP-binding cassette subfamily G member 2) is a broad substrate-specificATP-binding cassette transporter that plays an important role in porphyrin homeostasis and the cellular export of heme [31]. HBEGF (heparin binding EGF-like growth factor) is a Growth factor that may be involved in macrophage-mediated cellular proliferation and fibroblast mitosis [32]. UCHL1 (Ubiquitin carboxyl-terminal hydrolase isozyme L1) involved both in the processing of ubiquitin precursors and of ubiquitinated proteins. [33]. VSNL1 (Visinin-like protein 1) regulates the inhibition of rhodopsin phosphorylation in a calcium-dependent manner [34]. GABRB1 (Gamma-aminobutyric acid receptor subunit beta-1) and GABRA4 (Gamma-aminobutyric acid receptor subunit alpha-4) are both components of receptor for GABA, which functions as a ligand-gated chloride channel, and play an important role in ion transport [35,36]. RGS4 (Regulator of G-protein signaling 4) inhibits signal transduction by increasing the GTPase activity of G protein alpha subunits, which is associated with negative regulation of cardiac muscle cell development [37]. PPI network showed that these hub genes were mainly clustered in regulation of WNT signaling pathway and GABA-gated chloride ion channel activity. Consistent with our study, several studies have reported the downregulation of the Wnt pathway in animal models of DCM. In contrast, activation of the Wnt pathway in different ways helps prevent DCM development [38,39,40]. Our findings suggest that the suppression of the Wnt signaling pathway may play an important role in the progression of VMC to DCM. These hub genes may be potential therapeutic targets for the protection against progression of VMC to DCM.
As shown in Table 7, the top 10 agents predicted by cMAP for these hub genes were more concentrated in the RAF/MEK/MAPK signaling axis. To evaluate the feasibility of these drugs, HG-6-64-01 has not yet entered clinical trials. Selamectin is a potent nematocide but is highly toxic (neurotoxic). Tetrindole, as a selective inhibitor of monoamine oxidase A, is a new type of antidepressant but has a marked negative effect on detoxification. PD198306 is a selective inhibitor of MEK-1/2, which has not been used in clinical practice. Elvitegravir is a human immunodeficiency virus (HIV) integrase inhibitor which is used largely in a four-drug combination with cobicistat, emtricitabine, and tenofovir as therapy of HIV infection, but no cardiovascular applications have been reported. Az-628 is an RAF inhibitor that has not yet entered clinical trials. Amsacrine is an acridine derivative that is active in the treatment of acute leukemias and lymphomas, but no cardiovascular applications have been reported. GSK461364 is an ATP-competitive inhibitor of polo-like kinase 1 (Plk1) that currently has only one completed phase I clinical trial for non-Hodgkin’s lymphoma. Gefitinib, an EGFR inhibitor, is widely used in the treatment of non-small cell lung cancer but not cardiovascular disease. Amlodipine, as a calcium channel blocker, is one of the most widely used antihypertensive drugs. Although calcium channel blockers have not been shown to be beneficial for the treatment of patients with heart failure, several studies have shown a favorable effect of amlodipine on nonischemic dilated cardiomyopathy [41,42]. In addition, amlodipine has been reported to have a protective effect on myocardial injury in mice with congestive heart failure caused by viral myocarditis, and this therapeutic effect may be in part resulting from the inhibition of NO overproduction [42]. In conclusion, only amlodipine is currently viable among these potential therapeutic agents. Clinical trials are needed to verify its effectiveness in preventing the progression of viral myocarditis to dilated cardiomyopathy.
Because of evidence of a strong association between severe myocarditis and dilated cardiomyopathy, we created a subgroup of patients with severe myocarditis. Excitingly, the subgroup analysis yielded some novel findings. Firstly, the common DEGs of cVMC and DCM were enriched in the collagen-containing extracellular matrix, which was not found before. Currently, there is evidence that extracellular matrix remodeling plays an important role in dilated cardiomyopathy [43]. Our study suggests that viral myocarditis may develop dilated cardiomyopathy through changes in extracellular matrix related genes. Second, AQP4, VSNL1, GABRA4, and GABRB1 were still among the hub genes obtained by subgroup analysis, indicating that these genes were strongly associated with dilated cardiomyopathy caused by myocarditis. Among the newly discovered hub genes, ACTG2 was up-regulated and the rest were down-regulated. ACTG2 is involved in coding for actin, which is a major component of the sarcomere and plays an important role in cardiac contraction [44]. SCN1A encodes Sodium channel protein type 1 subunit alpha, which mediates the membrane permeability of voltage-dependent sodium channels and plays an important role in the water–sodium balance [45]. NFASC encodes Neurofascin, an ankyrin-binding protein that promotes cell adhesion [46]. Myosin-6 encoded by MYH6 is an actin-based motor protein that plays an important role in myocardial contraction [47]. Col1A1 encodes the collagen alpha-1 (1) chain, which plays an important role in extracellular matrix and fiber formation [48]. CSRP2 encodes cysteine and glycine-rich protein 2, which promotes smooth muscle cell proliferation and dedifferentiation [49]. Thirdly, the PPI network showed that these genes were clustered in the GABA-A receptor complex and myosin filaments. This is consistent with our previous findings of global hub genes clustering at GABA receptors, whereas myosin filament clustering is a novel finding. Combined with the functions of these hub genes, we can infer that severe myocarditis develops dilated cardiomyopathy through genetic changes related to sarcomere dysfunction and extracellular matrix remodeling. However, these need to be verified experimentally.
For dilated cardiomyopathy, the known and recognized mutated genes need to be taken seriously. Currently the major mutations identified in familial dilated cardiomyopathy include TIN, LMNA, MYH7, MYH6, TNNT2, ACTC1, BAG3, DSP, MYBPC3, RBM20, SCN5A, and TPM1. Our differential gene analysis did detect MYH7 and MYH6, and MYH6 was also among the hub genes of the severe myocarditis cohort. However, no other known genes were detected, probably because most of the included patients in the original study producing GSE17800 had sporadic dilated cardiomyopathy.
There are several limitations in this study. Firstly, the data set for myocarditis in this study used parvovirus B19 myocarditis samples, and whether it is representative is debatable. However, this is the largest genome-wide sample of viral myocarditis currently available, and the Venn method of screening for common DEGs has itself largely eliminated confounding factors such as different viral infections. Secondly, there were no RNAseq data, which limits the detection of moderate/weak expression changes. Unfortunately, there are currently no high-throughput datasets from myocardial biopsy samples from human viral myocarditis. The lack of RNAseq data does reduce the sensitivity of this study, but in turn increases the specificity of the detection of differential expression. Thirdly, our study lacks protein expression data and thus the most direct representation of the biological process of VMC to DCM. We will add the content of proteomics in a subsequent study. Fourthly, although we screened out potential therapeutic compounds targeting hub genes, most of them seemed infeasible, and thus further verification by animal experiments is required.

5. Conclusions

In conclusion, our study is the first to identify hub genes in the progression of VMC to DCM, and reveals the important role of GABA-gated chloride ion channel activity, RAP1 signaling pathway, and WNT signaling pathway. Hub genes of severe VMC and DCM include VSNL1, GABRA4, GABRB1, ACTG2, SCN1A, NFASC, AQP4, MYH6, COL1A1, and CSRP2, which mainly cluster in extracellular matrix remodeling and sarcomere dysfunction. These genes can be considered as genetic markers of early high-risk myocarditis. Amlodipine is a potential viable drug in preventing the progression of VMC to DCM. Further basic studies are needed to verify the role of these genes and pathways in the progression of VMC to DCM.

Author Contributions

Conceptualization, J.Z.; methodology, M.X.; software, M.X.; validation, Y.Z.; formal analysis, M.X.; investigation, M.X.; resources, M.X.; data curation, M.X.; writing—original draft preparation, J.Z.; writing—review and editing, T.C.; visualization, M.X.; supervision, Y.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from National Natural Science Foundation of China (81873486), the Science and Technology Development Program of Jiangsu Province-Clinical Frontier Technology (BE2022754), Clinical Medicine Expert Team (Class A) of Jinji Lake Health Talents Program of Suzhou Industrial Park(SZYQTD202102), Suzhou Key Discipline for Medicine(SZXK202129), Demonstration of Scientific and Technological Innovation Project (SKY2021002), Suzhou Dedicated Project on Diagnosis and Treatment Technology of Major Diseases(LCZX202132), Research on Collaborative Innovation of medical engineering combination(SZM2021014). The funders had no roles in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

GSE4172: accessed on 5 June 2022, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE4172. GSE17800: accessed on 5 June 2022, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi.

Acknowledgments

The providers of the datasets GSE4172 and GSE17800.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pollack, A.; Kontorovich, A.R.; Fuster, V.; Dec, G.W. Viral myocarditis--diagnosis, treatment options, and current controversies. Nat. Rev. Cardiol. 2015, 12, 670–680. [Google Scholar] [CrossRef] [PubMed]
  2. Sole, M.J.; Liu, P. Viral myocarditis: A paradigm for understanding the pathogenesis and treatment of dilated cardiomyopathy. J. Am. Coll. Cardiol. 1993, 22, 99A–105A. [Google Scholar] [CrossRef]
  3. Schultheiss, H.-P.; Fairweather, D.; Caforio, A.L.P.; Escher, F.; Hershberger, R.; Lipshultz, S.E.; Liu, P.P.; Matsumori, A.; Mazzanti, A.; McMurray, J.; et al. Dilated cardiomyopathy. Nat. Rev. Dis. Prim. 2019, 5, 32. [Google Scholar] [CrossRef] [PubMed]
  4. Merlo, M.; Cannatà, A.; Gobbo, M.; Stolfo, D.; Elliott, P.M.; Sinagra, G. Evolving concepts in dilated cardiomyopathy. Eur. J. Heart Fail. 2017, 20, 228–239. [Google Scholar] [CrossRef] [Green Version]
  5. Yajima, T.; Knowlton, K.U. Viral myocarditis: From the perspective of the virus. Circulation 2009, 119, 2615–2624. [Google Scholar] [CrossRef] [Green Version]
  6. Kawai, C.; Matsumori, A.; Fujiwara, H. Myocarditis and Dilated Cardiomyopathy. Annu. Rev. Med. 1987, 38, 221–239. [Google Scholar] [CrossRef]
  7. Reyes, M.P.; Lerner, A.M. Coxsackievirus myocarditis—With special reference to acute and chronic effects. Prog. Cardiovasc. Dis. 1985, 27, 373–394. [Google Scholar] [CrossRef]
  8. Kereiakes, D.J.; Parmley, W.W. Myocarditis and cardiomyopathy. Am. Heart J. 1984, 108, 1318–1326. [Google Scholar] [CrossRef]
  9. Shi, Y.; Chen, C.; Lisewski, U.; Wrackmeyer, U.; Radke, M.; Westermann, D.; Sauter, M.; Tschöpe, C.; Poller, W.; Klingel, K.; et al. Cardiac Deletion of the Coxsackievirus-Adenovirus Receptor Abolishes Coxsackievirus B3 Infection and Prevents Myocarditis In Vivo. J. Am. Coll. Cardiol. 2009, 53, 1219–1226. [Google Scholar] [CrossRef] [Green Version]
  10. Kaya, Z.; Katus, H.A.; Rose, N.R. Cardiac troponins and autoimmunity: Their role in the pathogenesis of myocarditis and of heart failure. Clin. Immunol. 2010, 134, 80–88. [Google Scholar] [CrossRef]
  11. Wu, R.Z.; Wu, H.P.; Wang, Z.Q.; Rong, X.; Jiang, J.B.; Qiu, H.X.; Chen, Q.; Xiang, R.L. Effects of carvedilol on cardiomyocyte apoptosis in autoimmune myocarditis in mice. Zhonghua Er Ke Za Zhi 2011, 49, 376–380. [Google Scholar] [PubMed]
  12. Wittchen, F.; Suckau, L.; Witt, H.; Skurk, C.; Lassner, D.; Fechner, H.; Sipo, I.; Ungethüm, U.; Ruiz, P.; Pauschinger, M.; et al. Genomic expression profiling of human inflammatory cardiomyopathy (DCMi) suggests novel therapeutic targets. Klin. Wochenschr. 2006, 85, 257–271. [Google Scholar] [CrossRef] [Green Version]
  13. Ameling, S.; Herda, L.R.; Hammer, E.; Steil, L.; Teumer, A.; Trimpert, C.; Dörr, M.; Kroemer, H.K.; Klingel, K.; Kandolf, R.; et al. Myocardial gene expression profiles and cardiodepressant autoantibodies predict response of patients with dilated cardiomyopathy to immunoadsorption therapy. Eur. Heart J. 2012, 34, 666–675. [Google Scholar] [CrossRef] [PubMed]
  14. Pauschinger, M.; Doerner, A.; Kuehl, U.; Schwimmbeck, P.L.; Poller, W.; Kandolf, R.; Schultheiss, H.-P. Enteroviral RNA Replication in the Myocardium of Patients With Left Ventricular Dysfunction and Clinically Suspected Myocarditis. Circulation 1999, 99, 889–895, Erratum in: Circulation 2010, 122, e388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Cooper, L.T. Myocarditis. N. Engl. J. Med. 2009, 360, 1526–1538. [Google Scholar] [CrossRef] [Green Version]
  16. Samuels, M.A. The Brain–Heart Connection. Circulation 2007, 116, 77–84. [Google Scholar] [CrossRef]
  17. Kimura, K.; Ieda, M.; Fukuda, K. Development, Maturation, and Transdifferentiation of Cardiac Sympathetic Nerves. Circ. Res. 2012, 110, 325–336. [Google Scholar] [CrossRef]
  18. Ohshima, S.; Isobe, S.; Izawa, H.; Nanasato, M.; Ando, A.; Yamada, A.; Yamada, K.; Kato, T.S.; Obata, K.; Noda, A.; et al. Cardiac Sympathetic Dysfunction Correlates With Abnormal Myocardial Contractile Reserve in Dilated Cardiomyopathy Patients. J. Am. Coll. Cardiol. 2005, 46, 2061–2068. [Google Scholar] [CrossRef] [Green Version]
  19. Kasama, S.; Toyama, T.; Hatori, T.; Sumino, H.; Kumakura, H.; Takayama, Y.; Ichikawa, S.; Suzuki, T.; Kurabayashi, M. Evaluation of cardiac sympathetic nerve activity and left ventricular remodelling in patients with dilated cardiomyopathy on the treatment containing carvedilol. Eur. Heart J. 2007, 28, 989–995. [Google Scholar] [CrossRef]
  20. Parthenakis, F.I.; Patrianakos, A.; Prassopoulos, V.; Papadimitriou, E.; Nikitovic, D.; Karkavitsas, N.S.; Vardas, P.E. Relation of cardiac sympathetic innervation to proinflammatory cytokine levels in patients with heart failure secondary to idiopathic dilated cardiomyopathy. Am. J. Cardiol. 2003, 91, 1190–1194. [Google Scholar] [CrossRef]
  21. Jeyaraj, S.C.; Unger, N.T.; Chotani, M.A. Rap1 GTPases: An emerging role in the cardiovasculature. Life Sci. 2011, 88, 645–652. [Google Scholar] [CrossRef] [Green Version]
  22. Wang, Y.; Wei, C.-X.; Shao, L.-Q.; Zhao, M. MiRNA Signaling in Viral Myocarditis Novel and Unique Pathological Features. Acta Cardiol. Sin. 2018, 34, 77–86. [Google Scholar] [CrossRef] [PubMed]
  23. Dong, W.; Yang, Z.; Yang, F.; Wang, J.; Zhuang, Y.; Xu, C.; Zhang, B.; Tian, X.-L.; Liu, D. Suppression of Rap1 Impairs Cardiac Myofibrils and Conduction System in Zebrafish. PLoS ONE 2012, 7, e50960. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Katoh, Y.; Katoh, M. FGFR2-related pathogenesis and FGFR2-targeted therapeutics (Review). Int. J. Mol. Med. 2009, 23, 307–311. [Google Scholar] [CrossRef] [Green Version]
  25. Vandebroek, A.; Yasui, M. Regulation of AQP4 in the Central Nervous System. Int. J. Mol. Sci. 2020, 21, 1603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Oklü, R.; Hesketh, R. The latent transforming growth factor beta binding protein (LTBP) family. Biochem. J. 2000, 352, 601–610. [Google Scholar] [CrossRef] [PubMed]
  27. Andriopoulos, B., Jr.; Corradini, E.; Xia, Y.; Faasse, S.A.; Chen, S.; Grgurevic, L.; Knutson, M.D.; Pietrangelo, A.; Vukicevic, S.; Lin, H.Y.; et al. BMP6 is a key endogenous regulator of hepcidin expression and iron metabolism. Nat. Genet. 2009, 41, 482–487. [Google Scholar] [CrossRef] [Green Version]
  28. Su, F.; Overholtzer, M.; Besser, D.; Levine, A.J. WISP-1 attenuates p53-mediated apoptosis in response to DNA damage through activation of the Akt kinase. Genes Dev. 2002, 16, 46–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Berndt, T.; Craig, T.A.; Bowe, A.E.; Vassiliadis, J.; Reczek, D.; Finnegan, R.; De Beur, S.M.J.; Schiavi, S.C.; Kumar, R. Secreted frizzled-related protein 4 is a potent tumor-derived phosphaturic agent. J. Clin. Investig. 2003, 112, 785–794. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, S.; Krinks, M.; Lin, K.; Luyten, F.P.; Moos, M. Frzb, a Secreted Protein Expressed in the Spemann Organizer, Binds and Inhibits Wnt-8. Cell 1997, 88, 757–766. [Google Scholar] [CrossRef]
  31. Suzuki, M.; Suzuki, H.; Sugimoto, Y.; Sugiyama, Y. ABCG2 Transports Sulfated Conjugates of Steroids and Xenobiotics. J. Biol. Chem. 2003, 278, 22644–22649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Smith, E.R.; Hewitson, T.D. HBEGF: An EGF-like growth factor with FGF23-like activity? Kidney Int. 2021, 99, 539–542. [Google Scholar] [CrossRef] [PubMed]
  33. Bilguvar, K.; Tyagi, N.K.; Ozkara, C.; Tuysuz, B.; Bakircioglu, M.; Choi, M.; Delil, S.; Caglayan, A.O.; Baranoski, J.F.; Erturk, O.; et al. Recessive loss of function of the neuronal ubiquitin hydrolase UCHL1 leads to early-onset progressive neurodegeneration. Proc. Natl. Acad. Sci. USA 2013, 110, 3489–3494. [Google Scholar] [CrossRef] [Green Version]
  34. Polymeropoulos, M.H.; Ide, S.; Soares, M.B.; Lennon, G.G. Sequence Characterization and Genetic Mapping of the Human VSNL1 Gene, a Homologue of the Rat Visinin-like Peptide RNVP1. Genomics 1995, 29, 273–275. [Google Scholar] [CrossRef]
  35. Polymeropoulos, M.; Rath, D.; Xiao, H.; Merril, C. Dinucleotide repeat polymorphism at the human β1 subunit of the GABAAreceptor gene (GABRB1). Nucleic Acids Res. 1991, 19, 6345. [Google Scholar] [CrossRef] [Green Version]
  36. Mclean, P.J.; Farb, D.H.; Russek, S.J. Mapping of the α4 subunit gene (GABRA4) to human chromosome 4 defines an α2—α4—β1—γ1 gene cluster: Further evidence that modern GABAA receptor gene clusters are derived from an ancestral cluster. Genomics 1995, 26, 580–586. [Google Scholar] [CrossRef]
  37. Guo, J.; Hang, P.; Yu, J.; Li, W.; Zhao, X.; Sun, Y.; Fan, Z.; Du, Z. The association between RGS4 and choline in cardiac fibrosis. Cell Commun. Signal. 2021, 19, 46. [Google Scholar] [CrossRef] [PubMed]
  38. Qiu, Z.; Chen, W.; Liu, Y.; Jiang, B.; Yin, L.; Chen, X. LncRNA AC061961.2 overexpression inhibited endoplasmic reticulum stress induced apoptosis in dilated cardiomyopathy rats and cardiomyocytes via activating wnt/β-catenin pathway. J. Recept. Signal Transduct. 2020, 41, 494–503. [Google Scholar] [CrossRef]
  39. Lu, D.; Bao, D.; Dong, W.; Liu, N.; Zhang, X.; Gao, S.; Ge, W.; Gao, X.; Zhang, L. Dkk3 prevents familial dilated cardiomyopathy development through Wnt pathway. Lab. Investig. 2015, 96, 239–248. [Google Scholar] [CrossRef] [Green Version]
  40. Le Dour, C.; Macquart, C.; Sera, F.; Homma, S.; Bonne, G.; Morrow, J.P.; Worman, H.J.; Muchir, A. Decreased WNT/β-catenin signalling contributes to the pathogenesis of dilated cardiomyopathy caused by mutations in the lamin a/C gene. Hum. Mol. Genet. 2016, 26, 333–343. [Google Scholar] [CrossRef]
  41. Zhang, X.; Hintze, T.H. Amlodipine releases nitric oxide from canine coronary microvessels: An unexpected mechanism of action of a calcium channel-blocking agent. Circulation 1998, 97, 576–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Wang, W.Z.; Matsumori, A.; Yamada, T.; Shioi, T.; Okada, I.; Matsui, S.; Sato, Y.; Suzuki, H.; Shiota, K.; Sasayama, S. Beneficial effects of amlodipine in a murine model of congestive heart failure induced by viral myocarditis. A possible mechanism through inhibition of nitric oxide production. Circulation 1997, 95, 245–251. [Google Scholar] [CrossRef] [PubMed]
  43. Parikh, M.; Seidman, M.; Oudit, G.Y. Critical Role of Extracellular Matrix Remodelling in Patients With Dilated Cardiomyopathy: Lessons From Connective Tissue Disorders. Can. J. Cardiol. 2022, 38, 309–310. [Google Scholar] [CrossRef]
  44. Wang, W.; Wang, Y.; Piao, H.; Li, B.; Zhu, Z.; Li, D.; Wang, T.; Liu, K. Bioinformatics Analysis Reveals MicroRNA-193a-3p Regulates ACTG2 to Control Phenotype Switch in Human Vascular Smooth Muscle Cells. Front. Genet. 2021, 11, 572707. [Google Scholar] [CrossRef] [PubMed]
  45. Ozturk, S.; Yetkin, E. SCN1A Mutation or Cross Talk? The Connection Between the Heart and Brain. Pediatr. Neurol. 2016, 63, e3. [Google Scholar] [CrossRef]
  46. Monfrini, E.; Straniero, L.; Bonato, S.; Compagnoni, G.M.; Bordoni, A.; Dilena, R.; Rinchetti, P.; Silipigni, R.; Ronchi, D.; Corti, S.; et al. Neurofascin (NFASC) gene mutation causes autosomal recessive ataxia with demyelinating neuropathy. Park. Relat. Disord. 2019, 63, 66–72. [Google Scholar] [CrossRef]
  47. Heissler, S.M.; Selvadurai, J.; Bond, L.M.; Fedorov, R.; Kendrick-Jones, J.; Buss, F.; Manstein, D.J. Kinetic properties and small-molecule inhibition of human myosin-6. FEBS Lett. 2012, 586, 3208–3214. [Google Scholar] [CrossRef] [Green Version]
  48. Pan, X.; Chen, X.; Ren, Q.; Yue, L.; Niu, S.; Li, Z.; Zhu, R.; Chen, X.; Jia, Z.; Zhen, R.; et al. Single-cell transcriptomics identifies Col1a1 and Col1a2 as hub genes in obesity-induced cardiac fibrosis. Biochem. Biophys. Res. Commun. 2022, 618, 30–37. [Google Scholar] [CrossRef]
  49. Sagave, J.F.; Moser, M.; Ehler, E.; Weiskirchen, S.; Stoll, D.; Günther, K.; Büttner, R.; Weiskirchen, R. Targeted disruption of the mouse Csrp2gene encoding the cysteine- and glycine-rich LIM domain protein CRP2 result in subtle alteration of cardiac ultrastructure. BMC Dev. Biol. 2008, 8, 80. [Google Scholar] [CrossRef]
Figure 1. Differentially expressed genes (DEG) analysis on viral myocarditis (VMC). (A) Volcanic plots of gene expression of VMC in GSE4172 using ggplot2 package of R 3.6.3. Red represents upregulated DEGs, blue represents downregulated DEGs, grey represents genes which are not differentially expressed. (B) Cross comparability evaluation of microarray data using ggplot2 package of R 3.6.3. (C) Heat map of the former 50 DEGs using ComplexHeatmap package of R 3.6.3.
Figure 1. Differentially expressed genes (DEG) analysis on viral myocarditis (VMC). (A) Volcanic plots of gene expression of VMC in GSE4172 using ggplot2 package of R 3.6.3. Red represents upregulated DEGs, blue represents downregulated DEGs, grey represents genes which are not differentially expressed. (B) Cross comparability evaluation of microarray data using ggplot2 package of R 3.6.3. (C) Heat map of the former 50 DEGs using ComplexHeatmap package of R 3.6.3.
Jcdd 09 00353 g001
Figure 2. Differentially expressed genes (DEG) analysis on dilated cardiomyopathy (DCM). (A) Volcanic plots of gene expression of DCM in GSE17800 using ggplot2 package of R 3.6.3. Red represents upregulated DEGs, blue represents downregulated DEGs, grey represents genes which are not differentially expressed. (B) Cross comparability evaluation of microarray data using ggplot2 package of R 3.6.3. (C) Heat map of the former 50 DEGs using ComplexHeatmap package of R 3.6.3.
Figure 2. Differentially expressed genes (DEG) analysis on dilated cardiomyopathy (DCM). (A) Volcanic plots of gene expression of DCM in GSE17800 using ggplot2 package of R 3.6.3. Red represents upregulated DEGs, blue represents downregulated DEGs, grey represents genes which are not differentially expressed. (B) Cross comparability evaluation of microarray data using ggplot2 package of R 3.6.3. (C) Heat map of the former 50 DEGs using ComplexHeatmap package of R 3.6.3.
Jcdd 09 00353 g002
Figure 3. Venn analysis of common differentially expressed genes (DEG) of VMC and DCM using ggplot2 package of R 3.6.3. (A) Venn diagram of common DEGs of VMC and DCM. (B) Venn diagram of common up-regulated DEGs of VMC and DCM. (C) Venn diagram of common down-regulated DEGs of VMC and DCM. (Filter criteria: |log2FC| > 0.58; p-Value < 0.05. For upregulated DEGs, the log2FC >0; For downregulated DEGs, the log2FC <0).
Figure 3. Venn analysis of common differentially expressed genes (DEG) of VMC and DCM using ggplot2 package of R 3.6.3. (A) Venn diagram of common DEGs of VMC and DCM. (B) Venn diagram of common up-regulated DEGs of VMC and DCM. (C) Venn diagram of common down-regulated DEGs of VMC and DCM. (Filter criteria: |log2FC| > 0.58; p-Value < 0.05. For upregulated DEGs, the log2FC >0; For downregulated DEGs, the log2FC <0).
Jcdd 09 00353 g003
Figure 4. Results of GO enrichment analysis and KEGG pathway analysis for common differentially expressed genes (DEG) of VMC and DCM using clusterProfiler package of R 3.6.3. (A) Results of GO enrichment analysis and KEGG pathway analysis for common up-regulated DEGs of VMC and DCM. (B) Results of GO enrichment analysis and KEGG pathway analysis for common down-regulated DEGs of VMC and DCM.
Figure 4. Results of GO enrichment analysis and KEGG pathway analysis for common differentially expressed genes (DEG) of VMC and DCM using clusterProfiler package of R 3.6.3. (A) Results of GO enrichment analysis and KEGG pathway analysis for common up-regulated DEGs of VMC and DCM. (B) Results of GO enrichment analysis and KEGG pathway analysis for common down-regulated DEGs of VMC and DCM.
Jcdd 09 00353 g004
Figure 5. Results of the protein–protein interaction (PPI) network analysis of common differentially expressed genes (DEG) of VMC and DCM constructed using STRING.
Figure 5. Results of the protein–protein interaction (PPI) network analysis of common differentially expressed genes (DEG) of VMC and DCM constructed using STRING.
Jcdd 09 00353 g005
Figure 6. Identification of common differentially expressed hub genes (DEG) of VMC and DCM using Cytoscape and MCODE plugin.
Figure 6. Identification of common differentially expressed hub genes (DEG) of VMC and DCM using Cytoscape and MCODE plugin.
Jcdd 09 00353 g006
Figure 7. Subgroup analysis for severe myocarditis (sVMC). (A) Venn plot of common differentially expressed genes (DEGs) of sVMC and DCM. (B) Results of GO/KEGG enrichment analysis of the common DEGs of sVMC and DCM. (C) Identification of hub common DEGs of sVMC and DCM by using Cytoscape and MCODE plugin.
Figure 7. Subgroup analysis for severe myocarditis (sVMC). (A) Venn plot of common differentially expressed genes (DEGs) of sVMC and DCM. (B) Results of GO/KEGG enrichment analysis of the common DEGs of sVMC and DCM. (C) Identification of hub common DEGs of sVMC and DCM by using Cytoscape and MCODE plugin.
Jcdd 09 00353 g007
Table 1. Detailed data of GSE4172 and GSE17800.
Table 1. Detailed data of GSE4172 and GSE17800.
Sequence Number of ChipGSE4172GSE17800
PlatformGPL570GPL570
Disease VMCDCM
Chip providerPatricia Ruiz Lab of Max Planck Institute for Molecular Genetics, GermanyFunktionelle Genomforschung LAb of Universitätsmedizin Greifswald, Germany
AddressIhnestrasse 65, BerlinJahnstraße 15a, Greifswald
Research objectHumanHuman
Experiment typeExpression profiling by arrayExpression profiling by array
Sample typeMyocardial biopsyMyocardial biopsy
Biopsy methodEMBs from the right ventricular septum by standard procedureEMBs in accordance with the Dallas criteria
Number of chip samples (used/total)12/1248/48
Number of cases/controls8/440/8
Clinical diagnosisPVB19 infected viral myocarditisDilated cardiomyopathy in chronic heart failure
Age (years)52.5 ± 1749 ± 10
Genetic DCM-None available
NYHA classification (n)
II
III
None available21
19
LVEF (%)42 ± 9.933 ± 6
Time of uploading chipPublic on 1 December 2012Public on 31 May 2006
Table 2. A part of the major differentially expressed genes in GSE4172.
Table 2. A part of the major differentially expressed genes in GSE4172.
GeneLog (Fold Change (FC))p-ValueAdjusted p-Value
(Adj. p-Value)
Ribosomal protein S4 Y-linked 1
(RPS4Y1)
4.641601440.0456743640.329086230
PRKAG2 antisense RNA 1
(PRKAG2-AS1)
3.301305935.04 ×10−60.036385392
Long intergenic non-protein coding RNA 619 (LINC00619)3.177810810.0013899130.111583803
Lysine demethylase 5D
(KDM5D)
3.078365730.0470575840.331496365
Cytochrome P450 family 7 subfamily A member 1 (CYP7A1)−2.96020190.0095826480.190923364
Adiponectin, C1Q and collagen domain containing (ADIPOQ)−2.90036380.0059910970.163809624
Selectin E (SELE)2.88204850.0016437190.116984950
Carboxyl ester lipase (CEL)−2.84750790.0017689020.116984950
Odd-skipped related transcription factor 1 (OSR1)−2.77891560.0012939610.110794568
Death domain containing 1 (DTHD1)−2.71526710.0044489660.150065986
Table 3. A part of the major differentially expressed genes in GSE17800.
Table 3. A part of the major differentially expressed genes in GSE17800.
GeneLog [Fold Change (FC)]p-ValueAdjusted p-Value
(Adj. p-Value)
Natriuretic peptide B
(NPPB)
4.050156830.0001151400.020903873
Natriuretic peptide A
(NPPA)
2.756228230.0005004000.033619529
Secretoglobin family 1D member 2
(SCGB1D2)
−2.579022886.78 × 10−50.019568598
BRAF-activated non-protein coding RNA (BANCR)2.418820110.0101474900.144758837
SH3 domain containing GRB2 like 2
(SH3GL2)
2.353135333.26 ×10−50.015148092
Dickkopf WNT signaling pathway inhibitor 2 (DKK2)2.193909350.0015087770.057371767
Tryptase alpha/beta 1 (TPSAB1)2.168416002.59 × 10−70.001867599
Tubulin tyrosine ligase like 2
(TTLL2)
−2.081520920.0002915710.027683044
Transmembrane protein with EGF like and two follistatin like domains 2
(TMEFF2)
−2.062657300.0034093970.08621005
Ankyrin repeat domain 33B
(ANKRD33B)
2.047543021.72 × 10−50.010465682
Table 4. A part of the common differentially expressed genes of VMC and DCM.
Table 4. A part of the common differentially expressed genes of VMC and DCM.
GeneGSE4172GSE17800
LogFCp-ValueAdj. p-ValueLogFCp-ValueAdj. p-Value
Secretoglobin family 1D member 2S
(SCGB1D2)
−1.754986170.042043050.31896103−2.579022886.78 × 10−50.01956859
Tubulin tyrosine ligase like 2 (TTLL2)−1.826732520.038905540.30862780−2.081520920.000291570.02768304
Carbohydrate sulfotransferase 9
(CHST9)
−2.195660550.016265510.22897841−1.635102550.002833990.07848970
Ankyrin repeat domain 33B (ANKRD33B)1.689000710.009592450.190923362.047543021.72 × 10−50.01046568
Attractin like 1
(ATRNL1)
1.910973250.001973450.118759411.728847009.96 × 10−70.00359621
leucine rich repeat transmembrane neuronal 4 (LRRTM4)−2.398783600.022268520.25514298−1.220081700.001916670.06546316
LOC101927256−2.291638920.001929500.11875941−1.300366080.007370640.12564681
C3orf67−2.016703340.016704820.23158961−1.560474720.005735560.11059997
Sphingosine-1-phosphate phosphatase 2 (SGPP2)−2.457056590.033096880.29296064−0.978538120.049711220.29650900
AC124997.1−1.498920260.046626490.33039811−1.881923870.000511510.03366817
VMC, viral myocarditis; DCM, dilated cardiomyopathy.
Table 5. GO enrichment and KEGG pathway analysis for the common up-regulated differentially expressed genes of VMC and DCM.
Table 5. GO enrichment and KEGG pathway analysis for the common up-regulated differentially expressed genes of VMC and DCM.
CategoryIDTermCountp-Value
BPGO:1903792negative regulation of anion transport3/484.99 × 10−5
GO:0042220response to cocaine3/483.81 × 10−4
GO:0060359response to ammonium ion4/483.98 × 10−4
GO:0051956negative regulation of amino acid transport2/484.20 × 10−4
GO:0060078regulation of postsynaptic membrane potential4/484.56 × 10−4
CCGO:1902711GABA-A receptor complex2/510.001
GO:1902710GABA receptor complex2/510.001
GO:0097730non-motile cilium3/510.007
GO:0034707chloride channel complex2/510.007
GO:0034702ion channel complex4/510.008
MFGO:0022851GABA-gated chloride ion channel activity2/485.51 × 10−4
GO:0042910xenobiotic transmembrane transporter activity2/485.51 × 10−4
GO:0004890GABA-A receptor activity2/480.001
GO:0099095ligand-gated anion channel activity2/480.001
GO:0016917GABA receptor activity2/480.002
KEGGhsa04727GABAergic synapse3/250.002
hsa05032Morphine addiction3/250.003
hsa05033Nicotine addiction2/250.007
hsa05219Bladder cancer2/250.007
Table 6. GO enrichment and KEGG pathway analysis for the common down-regulated differentially expressed genes of VMC and DCM.
Table 6. GO enrichment and KEGG pathway analysis for the common down-regulated differentially expressed genes of VMC and DCM.
CategoryIDTermCountp-Value
BPGO:0055010ventricular cardiac muscle tissue morphogenesis3/441.96 × 10−4
GO:0003229ventricular cardiac muscle tissue development3/442.94 × 10−4
GO:0055008cardiac muscle tissue morphogenesis3/445.74 × 10−4
GO:0003208cardiac ventricle morphogenesis3/446.77 × 10−4
GO:0048639positive regulation of developmental growth4/449.11 × 10−4
MFGO:0004714transmembrane receptor protein tyrosine kinase activity2/440.010
GO:0016829lyase activity3/440.011
GO:0005518collagen binding2/440.012
GO:0019199transmembrane receptor protein kinase activity2/440.016
GO:0038024cargo receptor activity2/440.019
KEGGhsa04015Rap1 signaling pathway4/230.003
Table 7. Top 10 prediction results from cMAP for the prevention against progression of VMC to DCM.
Table 7. Top 10 prediction results from cMAP for the prevention against progression of VMC to DCM.
ScoreNameDescriptionComment
95.81HG-6-64-01RAF inhibitor, Abl kinase inhibitor, ephrin receptor inhibitor, KIT inhibitor, MAP kinase inhibitor, MEK inhibitor, p38 MAPK inhibitor, src inhibitoruncharacterized
94.20selamectinnematocideneurotoxicity
92.97tetrindolemonoamine oxidase inhibitorantidepressant
92.63PD-198306MAP kinase inhibitor, MEK inhibitoruncharacterized
91.78elvitegravirHIV integrase inhibitor, HIV inhibitorclassic cocktail drug ingredients
91.60amsacrinetopoisomerase inhibitor, DNA intercalating drugin clinical trial for Acute Myeloid Leukemia
91.53GSK-461364PLK inhibitorin clinical trial for Non-Hodgkin Lymphoma
91.45gefitinibEGFR inhibitorfirst-generation targeted drug for non-small cell lung cancer
91.13amlodipinebreast cancer resistance protein inhibitor, calcium channel blocker, calcium channel inhibitor, L-type calcium channel blockerclassic CCB class antihypertensive drug
90.41AZ-628RAF inhibitoruncharacterized
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Xu, M.; Chen, T.; Zhou, Y. Bioinformatics Analysis of Common Differential Genes of Viral Myocarditis and Dilated Cardiomyopathy: Screening for Potential Pharmacological Compounds. J. Cardiovasc. Dev. Dis. 2022, 9, 353. https://doi.org/10.3390/jcdd9100353

AMA Style

Zhang J, Xu M, Chen T, Zhou Y. Bioinformatics Analysis of Common Differential Genes of Viral Myocarditis and Dilated Cardiomyopathy: Screening for Potential Pharmacological Compounds. Journal of Cardiovascular Development and Disease. 2022; 9(10):353. https://doi.org/10.3390/jcdd9100353

Chicago/Turabian Style

Zhang, Junyi, Mingzhu Xu, Tan Chen, and Yafeng Zhou. 2022. "Bioinformatics Analysis of Common Differential Genes of Viral Myocarditis and Dilated Cardiomyopathy: Screening for Potential Pharmacological Compounds" Journal of Cardiovascular Development and Disease 9, no. 10: 353. https://doi.org/10.3390/jcdd9100353

APA Style

Zhang, J., Xu, M., Chen, T., & Zhou, Y. (2022). Bioinformatics Analysis of Common Differential Genes of Viral Myocarditis and Dilated Cardiomyopathy: Screening for Potential Pharmacological Compounds. Journal of Cardiovascular Development and Disease, 9(10), 353. https://doi.org/10.3390/jcdd9100353

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

Article Metrics

Back to TopTop