Front. Med. Frontiers in Medicine Front. Med. 2296-858X Frontiers Media S.A. 10.3389/fmed.2023.1191354 Medicine Original Research Exploring potential biomarkers and therapeutic targets of long COVID-associated inflammatory cardiomyopathy Qi Peng1† Huang Mengjie2† Zhu Haiyan1* 1Department of Emergency, First Medical Center of Chinese PLA General Hospital, Beijing, China 2Department of Nephrology, First Medical Center of Chinese PLA General Hospital, Beijing, China

Edited by: Arthur Sá Ferreira, University Center Augusto Motta, Brazil

Reviewed by: Tejasav S. Sehrawat, Mayo Clinic, United States; Mohadeseh Poudineh, Zanjan University of Medical Sciences, Iran

 *Correspondence: Haiyan Zhu, zhuhaiyan301@sina.com

†These authors have contributed equally to this work and share first authorship

 29 06 2023 2023 10 1191354 22 03 2023 15 06 2023 Copyright © 2023 Qi, Huang and Zhu. 2023 Qi, Huang and Zhu

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Background

The negative impact of long COVID on social life and human health is increasingly prominent, and the elevated risk of cardiovascular disease in patients recovering from COVID-19 has also been fully confirmed. However, the pathogenesis of long COVID-related inflammatory cardiomyopathy is still unclear. Here, we explore potential biomarkers and therapeutic targets of long COVID-associated inflammatory cardiomyopathy.

 Methods

Datasets that met the study requirements were identified in Gene Expression Omnibus (GEO), and differentially expressed genes (DEGs) were obtained by the algorithm. Then, functional enrichment analysis was performed to explore the basic molecular mechanisms and biological processes associated with DEGs. A protein–protein interaction (PPI) network was constructed and analyzed to identify hub genes among the common DEGs. Finally, a third dataset was introduced for validation.

 Results

Ultimately, 3,098 upregulated DEGs and 1965 downregulated DEGs were extracted from the inflammatory cardiomyopathy dataset. A total of 89 upregulated DEGs and 217 downregulated DEGs were extracted from the dataset of convalescent COVID patients. Enrichment analysis and construction of the PPI network confirmed VEGFA, FOXO1, CXCR4, and SMAD4 as upregulated hub genes and KRAS and TXN as downregulated hub genes. The separate dataset of patients with COVID-19 infection used for verification led to speculation that long COVID-associated inflammatory cardiomyopathy is mainly attributable to the immune-mediated response and inflammation rather than to direct infection of cells by the virus.

 Conclusion

Screening of potential biomarkers and therapeutic targets sheds new light on the pathogenesis of long COVID-associated inflammatory cardiomyopathy as well as potential therapeutic approaches. Further clinical studies are needed to explore these possibilities in light of the increasingly severe negative impacts of long COVID.

 COVID-19 long COVID inflammatory cardiomyopathy bioinformatic analysis differentially expressed genes Ministry of Science and Technology10.13039/100007225 National Natural Science Foundation of China10.13039/501100001809 Beijing Natural Science Foundation10.13039/501100004826 section-at-acceptance Infectious Diseases: Pathogenesis and Therapy Background

COVID-19 is the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2); this disease has a significant impact on global public health and the health of individuals (1, 2). As of February 2023, the World Health Organization reported 750 million confirmed cases of COVID-19, including more than 6.8 million fatal cases. According to the World Health Organization, 80–85% of COVID-19 infections are mild or asymptomatic, with recovery similar to that of seasonal influenza. However, an increasing number of reports indicate that many patients experience persistent negative effects after recovery from the initial stage of the disease, resulting in long-term health problems (3, 4). This phenomenon was first described on Twitter as ‘long COVID’ in May 2020; other reports have described it as ‘postacute sequelae of COVID-19’ (5) or ‘chronic COVID syndrome’ (6). Long COVID is a multisystem disease that occurs in at least 10% of COVID-19 infections. More than 200 symptoms affecting multiple organ systems have been identified (5). Fatigue, dyspnea, cough, loss of sense of smell, brain fog and dysgeusia are the most common symptoms of long COVID, but there have also been reports of organ-specific injuries involving the cardiovascular and neuropsychiatric systems, as well as diabetes symptoms (7–10). Recent studies have advanced many hypotheses regarding the pathogenesis of long COVID (11–18), but the exact etiology of long COVID is still unclear. The importance of developing and verifying biomarkers that can be used to diagnose long COVID cannot be overemphasized (5). Such biomarkers would not only aid in diagnosis but also potentially provide targets for treatment and prevention. Inflammatory cardiomyopathy is one of the more common components of long COVID (19, 20). Typical symptoms and signs of inflammatory cardiomyopathy include chest pain, dyspnea, fatigue, palpitations, syncope, and cardiogenic shock. It can also manifest as sudden cardiac death (21), accounting for approximately 10% of sudden cardiac deaths in young individuals (age < 35 years) (22). Inflammatory cardiomyopathy is mainly mediated by viral infection but can also be induced by systemic immune-mediated disease (23). In a prospective study, 60% of convalescent COVID-19 patients had myocardial inflammation, and 78% had abnormal cardiac magnetic resonance imaging results (24). In another study involving 153,760 patients with COVID-19, the patients exhibited an increased cardiovascular disease event burden for up to 12 months, even those patients without any cardiovascular disease basis (10). These data suggest that a deeper understanding of the pathogenesis of cardiovascular disease in convalescent COVID-19 patients is needed, and the negative challenges posed by long COVID urgently require treatment approaches that are still lacking (25, 26). To explore the potential biomarkers and therapeutic targets of long COVID-related inflammatory cardiomyopathy, we screened Gene Expression Omnibus to find appropriate datasets and screened out DEGs. To explore the basic molecular mechanism of the biological processes of the DEGs, functional enrichment analysis was performed. Then, a protein–protein interaction network was constructed to obtain hub genes. Finally, a third dataset was used to verify and make reasonable inferences. The workflow of this study is shown in Figure 1. We used bioinformatics techniques to mine more valuable information from the raw data in public databases to inform the subsequent experimental research design. This approach can also provide a better understanding of the mechanism of long COVID-related inflammatory cardiomyopathy and more effectively promote rapid diagnosis and drug treatment in the future.

The workflow used in this study to explore potential biomarkers and therapeutic targets of long COVID-associated inflammatory cardiomyopathy.

 Methods Microarray data

We searched the GEO database with ‘SARS-CoV-2’ and ‘COVID-19’ as keywords, with ‘Experiment type’ limited to ‘Expression profiling by array’, and ‘Organism’ limited to ‘Homo sapiens’. From the original literature, we selected datasets that included both convalescent COVID-19 patients and healthy controls. Then, the datasets that met the principles of randomized controls and ethical requirements were retained, and finally, the GSE211378 dataset was selected. This dataset used the GPL32571 nCounter Host Response Panel v1.0, which included 160 convalescent COVID-19 patients (age:49.3 ± 6.2) with molecular or serological laboratory evidence of COVID-19 in the past and who had fully recovered from COVID-19. Except for loss of taste or smell, there is no other COVID-19 symptom that lasts for ≥28 days, or molecular detection that was negative for ≥14 days after recovery, and 40 healthy controls (age:50.9 ± 14.5) were selected. Then, the same search strategy was used to search for ‘inflammatory cardiomyopathy’ as the keyword, and the ‘GSE4172’ dataset, which included 8 patients with inflammatory cardiomyopathy. Inflammatory cardiomyopathy was confirmed through endocardial myocardial biopsy based on clinical, morphological and functional diagnostic criteria. The functional classification of the cases according to New York Heart Association was II. Four healthy controls were selected. The platform used was GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. The two datasets were downloaded through R software (version 4.2.1), 1 and the original data from the two datasets were read using the ‘affy’ package (27). The RMA algorithm was used for background correction and data normalization.

 Identification of DEGs

The ‘limma’ package (28) in R was used to group and label the disease groups and the healthy control groups in the two datasets. The DEGs between the groups were obtained. Assuming that the differences between convalescent COVID-19 patients and normal controls might be relatively small, the criteria were set as (1) | log2 (fold-change)| > 0.1 and (2) p < 0.05 to screen out more meaningful differentially expressed genes. The genes that met these conditions were identified as having significant differences (28, 29). Among the selected DEGs, genes with logFC >0.1 were considered to be upregulated, while genes with logFC < −0.1 were considered to be downregulated.

 Go and KEGG pathway enrichment analysis of DEGs

GO functional analysis is a systematic method for gene annotation in the context of RNA and protein expression and is a powerful bioinformatics tool for classifying gene expression and its characteristics (30). KEGG is an online database of cellular, enzymatic and biochemical pathways that is used to determine which pathways may be impacted by changes in gene expression (31). The clusterProfiler package (32) in R was used for DEG enrichment analysis. The ‘org.Hs.eg.db’ package in R was used for ID conversion (33), and the ‘ggplot2’ package was used to visualize the enrichment analysis results.

 Construction of the PPI network

The DEGs common to the two datasets were obtained by Venn diagram analysis, and then the co-upregulated DEGs and co-downregulated DEGs were uploaded to the STRING database (34) 2 to construct the PPI network, which is helpful for mining the core regulatory genes. Finally, the STRING results were imported into Cytoscape software (35) to test the potential correlation between these DEGs, and the results were outputted as a visual molecular interaction network.

 Screening of hub genes

Hub genes are defined as those proteins or DEGs that have higher degrees of association with others. The cytoHubba plug-in (36) in Cytoscape was used to identify hub genes among the common DEGs. Different plug-in algorithms may have different selection bias. Therefore, the DMNC, MCC, EPC, and MNC algorithms were used to calculate the hub genes, and the genes identified by all of the algorithms were defined as the final hub genes. The results are visualized in a Venn diagram.

 Analysis of hub genes

The Human Protein Atlas (37) 3 is a public database of gene expression profiles in different human organs from which the basic RNA and protein expression levels of specific genes can be identified. The PANTHER classification system (38) 4 has a comprehensive, annotated gene family system library that can classify proteins (and their genes). The selected hub genes were analyzed using the above databases.

 Verification of hub genes

Some reports have noted that the persistence of the virus may be a driving factor for long COVID, and it has even been shown that active virus or viral components can still be found in convalescent COVID-19 patients 12 months after the initial infection. In addition, inflammatory cardiomyopathy is mainly mediated by viral infection. Regardless of whether long COVID-associated inflammatory cardiomyopathy is mediated only by COVID-19 infection or also by the immune-mediated response and inflammation, the exact mechanism is still unclear. To further clarify this mechanism, we introduce data from a third set of COVID-19-infected patients for verification.

 Results Identification of DEGs

The two datasets were processed by the ‘limma’ package in R. In the GSE4172 dataset of inflammatory cardiomyopathy, the number of genes conforming to (1) log2 (fold-change) > 0.1 and (2) p <0.05 was 3,098, which was defined as upregulated DEGs. The number of genes conforming to (1) log2 (fold-change) < −0.1 and (2) p < 0.05 was 1965; these genes were defined as downregulated DEGs. Among the convalescent COVID-19 patients in the GSE211378 dataset, the number of genes conforming to (1) log2 (fold-change) > 0.1 and (2) p < 0.05 was 89; these were defined as upregulated DEGs. The number of genes conforming to (1) log2 (fold-change) < −0.1 and (2) p < 0.05 was 217; these genes were defined as downregulated DEGs. The results of differential analysis were visualized by the ‘ggplot2’ package (Figure 2).

DEGs in two datasets. (A) Volcano plot of GSE4172; red represents upregulated genes, blue represents downregulated genes, and gray represents nonsignificantly expressed genes. (B) Volcano plot of GSE211378. The criteria for statistically significant differences for DEGs were adjusted |log2(FC)| > 0.1 and p < 0.05 for expression.

 Go and KEGG pathway enrichment analyses of DEGs

Enrichment analysis of DEGs between the two datasets was performed with the clusterProfiler package in R. The biological process (BP) terms enriched among the DEGs in inflammatory cardiomyopathy (GSE4172) were mainly mRNA processing, generation of precursor metabolites and energy, RNA splicing, proteasome-mediated ubiquitin-dependent protein catabolic process and regulation of mRNA metabolic process. The enriched cell component (CC) terms were mainly mitochondrial matrix, nuclear speck, focal adhesion, cell-substrate junction and spliceosomal complex. The enriched molecular function (MF) terms were mainly ubiquitin protein ligase binding, ubiquitin-like protein ligase binding, transcription coregulator activity, oxidoreduction-driven active transmembrane, transporter activity and cytochrome-c oxidase activity. The enriched KEGG pathways were mainly nonalcoholic fatty liver disease, amyotrophic lateral sclerosis, protein processing in endoplasmic reticulum, Alzheimer’s disease and Huntington’s disease. The BP terms enriched among the DEGs of convalescent COVID-19 patients in GSE211378 were mainly cytokine-mediated signaling pathway, positive regulation of cytokine production, leukocyte proliferation, regulation of leukocyte proliferation and positive regulation of leukocyte proliferation. The enriched CC terms were mainly external side of plasma membrane, secretory granule membrane, membrane microdomain, membrane raft and plasma membrane signaling receptor complex. The enriched MF terms were mainly cytokine receptor binding, cytokine activity, cytokine receptor activity, immune receptor activity and receptor ligand activity. The enriched KEGG pathways were mainly cytokine–cytokine receptor interaction, lipid and atherosclerosis, JAK–STAT signaling pathway, viral protein interaction with cytokine and cytokine receptor and Th17 cell differentiation. Thus, many of the DEGs of the convalescent COVID-19 patients are related to cytokines, and research on the mechanism of long COVID can focus on this aspect. The enrichment results for the two datasets are visualized in the form of bubbles (Figure 3) and listed in a three-line table format (Tables 1, 2).

GO terms and KEGG pathway enrichment. (A) Details of GO terms and KEGG pathway enrichment for GSE4172. The colored dots represent the p value for that term, with red representing greater significance. The size of the dot represents the number of genes involved. (B) Details of GO terms and KEGG pathway enrichment for GSE211378.

Details of GO terms and KEGG pathway enrichment in GSE4172.

Ontology ID Description GeneRatio BgRatio p value p. Adjust
BP GO:0006397 mRNA processing 164/3742 500/18800 4.08e-12 2.57e-08
BP GO:0006091 Generation of precursor metabolites and energy 156/3742 494/18800 3.31e-10 1.04e-06
BP GO:0008380 RNA splicing 141/3742 438/18800 5.6e-10 1.17e-06
BP GO:0043161 Proteasome-mediated ubiquitin-dependent protein catabolic process 133/3742 414/18800 1.99e-09 3.14e-06
BP GO:1903311 Regulation of mRNA metabolic process 101/3742 294/18800 3.65e-09 4.6e-06
CC GO:0005759 Mitochondrial matrix 153/3883 473/19594 5.08e-11 1.4e-08
CC GO:0016607 Nuclear speck 137/3883 411/19594 5.2e-11 1.4e-08
CC GO:0005925 Focal adhesion 139/3883 419/19594 5.5e-11 1.4e-08
CC GO:0030055 Cell-substrate junction 140/3883 428/19594 1.42e-10 2.71e-08
CC GO:0005681 Spliceosomal complex 76/3883 194/19594 3.46e-10 5.28e-08
MF GO:0031625 Ubiquitin protein ligase binding 108/3824 298/18410 4.02e-10 2.98e-07
MF GO:0044389 Ubiquitin-like protein ligase binding 113/3824 317/18410 4.97e-10 2.98e-07
MF GO:0003712 Transcription coregulator activity 150/3824 497/18410 3.4e-07 0.0001
MF GO:0015453 Oxidoreduction-driven active transmembrane transporter activity 32/3824 72/18410 5.13e-06 0.0015
MF GO:0004129 Cytochrome-c oxidase activity 13/3824 19/18410 9.97e-06 0.0022
KEGG hsa04932 Non-alcoholic fatty liver disease 64/1779 155/8164 2.8e-08 6.86e-06
KEGG hsa05014 Amyotrophic lateral sclerosis 123/1779 364/8164 4.58e-08 6.86e-06
KEGG hsa04141 Protein processing in endoplasmic reticulum 68/1779 171/8164 6.16e-08 6.86e-06
KEGG hsa05010 Alzheimer disease 124/1779 384/8164 6.82e-07 4.19e-05
KEGG hsa05016 Huntington disease 103/1779 306/8164 7.24e-07 4.19e-05

Details of GO terms and KEGG pathway enrichment in GSE211378.

Ontology ID Description GeneRatio BgRatio p value p. Adjust
BP GO:0019221 Cytokine-mediated signaling pathway 107/298 486/18800 1.85e-94 7.49e-91
BP GO:0001819 Positive regulation of cytokine production 73/298 475/18800 3.89e-51 7.86e-48
BP GO:0070661 Leukocyte proliferation 59/298 330/18800 5.31e-45 7.16e-42
BP GO:0070663 Regulation of leukocyte proliferation 50/298 254/18800 3.08e-40 3.12e-37
BP GO:0070665 Positive regulation of leukocyte proliferation 42/298 158/18800 1.04e-39 8.4e-37
CC GO:0009897 External side of plasma membrane 47/299 455/19594 1.46e-25 4.82e-23
CC GO:0030667 Secretory granule membrane 27/299 312/19594 3.65e-13 6.01e-11
CC GO:0098857 Membrane microdomain 27/299 327/19594 1.12e-12 1.23e-10
CC GO:0045121 Membrane raft 26/299 326/19594 6.56e-12 5.4e-10
CC GO:0098802 Plasma membrane signaling receptor complex 25/299 313/19594 1.64e-11 1.08e-09
MF GO:0005126 Cytokine receptor binding 60/298 272/18410 5.35e-51 2.8e-48
MF GO:0005125 Cytokine activity 48/298 235/18410 5.11e-39 1.34e-36
MF GO:0004896 Cytokine receptor activity 33/298 97/18410 4.62e-35 6.82e-33
MF GO:0140375 Immune receptor activity 38/298 148/18410 5.22e-35 6.82e-33
MF GO:0048018 Receptor ligand activity 50/298 489/18410 7.6e-26 7.95e-24
KEGG hsa04060 Cytokine-cytokine receptor interaction 80/260 295/8164 9.7e-55 2.45e-52
KEGG hsa05417 Lipid and atherosclerosis 48/260 215/8164 4.95e-28 6.26e-26
KEGG hsa04630 JAK–STAT signaling pathway 42/260 162/8164 2.38e-27 2.01e-25
KEGG hsa04061 Viral protein interaction with cytokine and cytokine receptor 34/260 100/8164 1.35e-26 8.55e-25
KEGG hsa04659 Th17 cell differentiation 34/260 108/8164 2.53e-25 1.28e-23
Construction of PPI network

The upregulated and downregulated DEGs between the two datasets were identified and then analyzed using the ‘VennDiagram’ package in R to obtain 12 upregulated common DEGs and 16 downregulated common DEGs (Figure 4A). Then, the obtained common DEGs were imported into the STRING database to generate a PPI network, and the generated information was imported into Cytoscape software to visualize the PPI network (Figures 4B,C).

Venn diagram and PPI network. (A) The Venn diagram shows the co-DEGs of the two datasets. (B) A PPI network was constructed using co-upregulated genes. (C) A PPI network was constructed using co-downregulated genes.

 Screening of hub genes

By using the cytoHubba plug-in of Cytoscape to extract hub genes, to avoid bias caused by different algorithms, four algorithms, DMNC, MCC, EPC, and MNC, were used to process the information, and the four groups of results were merged through the ‘VennDiagram’ package to confirm the final set of hub genes (Figure 5). VEGFA, FOXO1, CXCR4 and SMAD4 were identified as upregulated hub genes, and KRAS and TXN were identified as downregulated hub genes.

Hub genes were merged by the MCC, DMNC, MNC, and EPC algorithms. (A) VEGFA, FOXO1, CXCR4 and SMAD4 were upregulated hub genes. (B) KRAS and TXN were downregulated hub genes.

 Hub gene analysis

To explore the cell type specificity of hub genes in the human myocardium, the expression specificity of hub genes was predicted by using the Human Protein Atlas. 5 The database prediction was generated by integrated network analysis of a large amount of publicly available RNAseq data; the results are shown in Figure 6. The PANTHER classification system 6 can be used for one-stop annotation analysis of hub genes. The results are listed in Table 3.

Human Protein Atlas was used to analyze hub genes. The plot shows the enrichment prediction score for each cell type profiled in heart muscle. This score is the mean correlation between the selected gene and the 3 reference transcripts selected to represent each cell type profiled within the tissue. (A) VEGFA, (B) FOXO1, (C) CXCR4, (D) SMAD4, (E) KRAS, and (F) TXN.

The PANTHER classification system of hub genes.

ID Gene name Gene ID PANTHER family/subfamily PANTHER protein class
VEGFA Vascular endothelial growth factor A HUMAN|HGNC = 12,680|UniProtKB=P15692 VASCULAR ENDOTHELIAL GROWTH FACTOR A (PTHR12025:SF5) Growth factor
FOXO1 Forkhead box protein O1 HUMAN|HGNC = 3,819|UniProtKB = Q12778 FORKHEAD BOX PROTEIN O1 (PTHR45767:SF1) Winged helix/forkhead transcription factor
CXCR4 C-X-C chemokine receptor type 4 HUMAN|HGNC = 2,561|UniProtKB=P61073 C-X-C CHEMOKINE RECEPTOR TYPE 4 (PTHR10489:SF594) Cell adhesion molecule
SMAD4 Mothers against decapentaplegic homolog 4 HUMAN|HGNC = 6,770|UniProtKB = Q13485 MOTHERS AGAINST DECAPENTAPLEGIC HOMOLOG 4 (PTHR13703:SF63) DNA-binding transcription factor
KRAS GTPase KRas HUMAN|HGNC = 6,407|UniProtKB=P01116 GTPASE KRAS (PTHR24070:SF388) Small GTPase
TXN Thioredoxin HUMAN|HGNC = 12,435|UniProtKB=P10599 THIOREDOXIN (PTHR10438:SF18) Oxidoreductase
Verifying hub genes by ROC curve analysis

To verify the diagnostic performance of the hub genes, we performed ROC curve analysis based on two GEO datasets. The results showed that all hub genes had good diagnostic performance in the inflammatory cardiomyopathy dataset, and the area under the ROC curve (AUC) was more than 85%. The relatively small difference between the patients who recovered from COVID-19 and the control group may be small may affect the diagnostic performance of the hub genes, but the results still showed that the area under the ROC curve of all hub genes in the COVID-19 recovery patient dataset was more than 60%. Therefore, we conclude that VEGFA, FOXO1, CXCR4, SMAD4, KRAS, and TXN are key genes in the common pathogenesis of long COVID and inflammatory cardiomyopathy, with great diagnostic value (Figure 7).

The GEO dataset was used to verify the diagnostic performance of hub genes. (A) ROC curve of hub genes in the GSE4172 dataset. (B) ROC curve of hub genes in the GSE211378 dataset.

 Discussion

Long COVID may be caused by a variety of factors. A series of studies have proposed a variety of hypotheses about the pathogenesis of long COVID, including the persistence of viruses in tissues (12, 39), immune dysfunction (13, 14), destruction of the microbiota (40, 41), autoimmune response (15–18), coagulation and endothelial dysfunction (17, 18), and neurological signal dysfunction (42). With the increasing number of patients recovering from acute COVID-19 infection, the negative impact of long COVID on social life and individual health has become increasingly prominent. It has been reported that at 7 months after COVID-19 infection, patients’ quality of life remains diminished; 45.2% of convalescent COVID-19 patients need a reduced workload, and 22.3% of patients cannot return to work at all (7). This requires not only clinical trials on potential biomarkers and therapeutic targets for long COVID as soon as possible but also extensive public campaigns to make the public aware of the risks and consequences of long COVID (5). Although there are currently some diagnostic tools for long COVID-related diseases (43, 44), most of them are still being developed. There has been only some preliminary progress in the study of biomarkers of long COVID (45, 46). At present, there is an extreme lack of effective means to diagnose long COVID and evaluate its response to treatment. In this study, we performed bioinformatics analysis on long COVID-associated inflammatory cardiomyopathy. Through raw data mining and analysis, we obtained potential effective biomarkers and therapeutic targets. To verify and further explore the pathogenesis of long COVID-related inflammatory cardiomyopathy, we first introduced a third dataset for parallel verification. The GEO database was searched using ‘SARS-CoV-2’ and ‘COVID-19’ as keywords. From the original literature, we selected the datasets including both patients with COVID-19 infection and healthy controls, identified the datasets that met the principles of randomized control and ethical requirements, and finally selected the dataset GSE164805. The GSE164805 dataset includes 10 patients with COVID-19 infection and 5 healthy controls. By comparing the COVID-19 infection dataset (GSE164805) with the inflammatory cardiomyopathy dataset (GSE4172), we identified RRP12, FTSJ3, DCAF13, WDR46, DHX32, NLE1, POLR3B, NSUN4, HNRNPK, and HNRNPM as upregulated hub genes and RPS27, FAU, RPL35, RPL35A, RPL38, RPS11, RPS28, RPS19, RPL19, and RPL18 as downregulated hub genes. Through the comparison of hub genes between the three sets of datasets, it can be seen that the inflammatory cardiomyopathy-related targets are distinct between patients with COVID-19 infection and convalescent COVID-19 patients. Therefore, it is reasonable to speculate that long COVID-associated inflammatory cardiomyopathy is mainly attributable to the immune-mediated response and inflammation rather than direct viral infection of cells. Several other relevant reports support this theoretical inference (5, 47, 48). Second. we identified four upregulated hub genes, VEGFA, FOXO1, CXCR4, and SMAD4, and two downregulated hub genes, KRAS and TXN. Through Human Protein Atlas and the PANTHER classification system, we gained a preliminary and comprehensive understanding of the hub genes. We also conducted a detailed literature survey to further explore the hub genes. According to the literature, vascular endothelial growth factor A (VEGFA) is a homodimeric vasoactive glycoprotein and a key mediator of angiogenesis (49). Cardiomyocytes are the main source of VEGFA and the target of VEGFA (50–54). Studies have found that VEGFA plays a key role in triggering cardiac angiogenesis after acute myocardial infarction (55), but these studies are still limited to animal experiments. At present, there are no data on the expression and release of VEGFA in human cardiomyocytes. Further studies are needed to determine the mechanism of VEGFA in humans and whether this effect is valuable for patients recovering from COVID-19. Forkhead box O1 (FOXO1) is a representative member of the FOXO family and has key transcriptional regulatory activity (56). FOXO1 can regulate a variety of targets, such as genes involved in apoptosis and autophagy, antioxidant enzymes, cell cycle arrest genes, and metabolic and immunomodulators (57, 58). Due to its extensive endogenous expression, it may be an important target for the treatment and prevention of various diseases (59). C-X-C motif chemokine receptor 4 (CXCR4) is a cell chemokine receptor that plays an important role in development, hematopoiesis and immune surveillance through signal transduction induced by its ligand CXCL12. The pleiotropic effects of CXCR4 are well illustrated by the complexity of its biological functions and are mediated by mechanisms including receptor crosstalk, receptor and ligand isoforms, and unconventional ligands. Although the CXCL12–CXCR4 signaling pathway has been extensively studied, its atypical pathways remain to be fully elucidated (60). In addition, CXCR4 has become an interesting target in autoimmunity and inflammation based on its involvement in the chemotaxis of leukocytes to inflammatory sites (61, 62), the inappropriate retention (63) of activated innate inflammatory cells at inflammatory sites and its regulatory roles in autoimmune diseases (64, 65). Due to the great value of CXCR4 as a therapeutic target, the CXCR4 antagonist AMD3100 has also been approved by the US FDA (66, 67). Through ongoing research, CXCR4 has become a target for drug research and in the development of various diseases. SMAD family member 4 (SMAD4) is a member of the Smad family of signal transduction proteins, which are phosphorylated and activated by transmembrane serine–threonine receptor kinases through several pathways in response to transforming growth factor beta (TGF-β) signaling. At present, research on SMAD4 is mainly focused on its role as a target of tumor therapy (68), but some studies have confirmed that this protein also has activity in the late stage of COVID-19 infection. The mechanism of the association between long COVID and the development of cardiovascular disease may be through the activation of TGF-β signaling through the Smad pathway to induce subsequent myocardial fibrosis and scar formation (10). Since SMAD4 is abundantly expressed in the heart and vascular system during embryogenesis and plays a key role in normal cardiovascular development, the absence of SMAD4 leads to various cardiovascular developmental defects (69, 70). Thus, whether SMAD4 plays a specific role in inflammatory cardiomyopathy and cardiovascular disease during recovery from COVID-19 deserves further experimental research. The proto-oncogene KRAS is considered to be the most common gene driver of oncogenesis in human cancers (71). Overactivation of KRAS signaling has been shown to enhance the secretion of interleukin-6, which is necessary for tumorigenesis and tumor development (72). However, whether KRAS is involved in the development of inflammatory cardiomyopathy in patients recovering from COVID-19 and what role it plays in this process have not been published. Thioredoxin (TXN) is a hydrogen carrier protein that is widely expressed in organisms. TXN deficiency is associated with myocardial cell injury (73). This mechanism of this association may be the induction of autophagy in TXN-deficient cardiomyocytes, which leads to severe inflammation and myocardial cell damage (74). In summary, previous studies have, to a certain extent, confirmed the effectiveness and potential value of the hub genes identified in this study.

One limitation of our study is that the evidence for the identified potential biomarkers and therapeutic targets is limited to the theoretical level, and further experiments are needed to verify their activities. Another limitation is that there is no publicly available gene expression profile dataset of COVID-19-induced myocarditis. We did examine a dataset of parvovirus B19-related inflammatory cardiomyopathy. However, considering that the pathogenesis of inflammatory cardiomyopathy induced by COVID-19 infection and inflammatory cardiomyopathy in patients recovering from COVID-19 is not exactly the same, we used bioinformatics methods to analyze and extract the information hidden in the intergroup data, which will help to obtain effective potential biomarkers and therapeutic targets for long COVID-related inflammatory cardiomyopathy. Due to the particularity of SARS-CoV-2 infection, the development of many clinical studies has been limited to varying degrees. This makes the role of bioinformatics analysis more prominent in the study of COVID-19. Some recent studies (75) have also predicted potential therapeutic drugs for the ischemic cardiomyopathy associated with COVID-19 through bioinformatics methods. This research provides meaningful exploration methods for the study of long COVID-related diseases.

 Conclusion

Based on our bioinformatics analysis and previous research, this study identified potential targets (VEGFA, FOXO1, CXCR4, SMAD4, KRAS, TXN) that may be helpful for the diagnosis and treatment of long COVID-related inflammatory cardiomyopathy. The screening of these targets provides new directions for the study of the pathogenesis and treatment of long COVID-related inflammatory cardiomyopathy. Further clinical research is needed to explore these possibilities to address the increasingly serious negative impacts of long COVID.

 Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

 Author contributions

PQ was the major contributor to the writing of the manuscript. MJH evaluated and revised the manuscript. HZ guided the writing of the manuscript and critically reviewed it. All authors contributed to the article and approved the submitted version.

 Funding

This research was supported by the National Key R&D Program of the Ministry of Science and Technology (2019YFC0119600), National Natural Science Foundation of China (82000631), Beijing Natural Science Foundation (7222169), Young Elite Scientist Sponsorship Program by CAST 2020QNRC001 (to MJH), and the Military Medical Youth Special Project of PLA General Hospital (QNF19035).

 Conflict of interest

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

 Publisher’s note

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

 References Global Burden of Disease 2020 Health Financing Collaborator Network . Tracking development assistance for health and for COVID-19: a review of development assistance, government, out-of-pocket, and other private spending on health for 204 countries and territories, 1990-2050. Lancet. (2021) 398:1317–43. doi: 10.1016/S0140-6736(21)01258-7, PMID: 34562388 Lal A Erondu NA Heymann DL Gitahi G Yates R. Fragmented health systems in COVID-19: rectifying the misalignment between global health security and universal health coverage. Lancet. (2021) 397:61–7. doi: 10.1016/S0140-6736(20)32228-5, PMID: 33275906 Al-Aly Z Bowe B Xie Y. Long COVID after breakthrough SARS-CoV-2 infection. Nat Med. (2022) 28:1461–7. doi: 10.1038/s41591-022-01840-0, PMID: 35614233 Ballering AV van Zon SKR Olde Hartman TC Rosmalen JGM. Persistence of somatic symptoms after COVID-19 in the Netherlands: an observational cohort study. Lancet. (2022) 400:452–61. doi: 10.1016/S0140-6736(22)01214-4, PMID: 35934007 Davis HE McCorkell L Vogel JM Topol EJ. Long COVID: major findings, mechanisms and recommendations. Nature reviews. Microbiology. (2023) 21:133–46. doi: 10.1038/s41579-022-00846-2 Baig AM . Chronic COVID syndrome: need for an appropriate medical terminology for long-COVID and COVID long-haulers. J Med Virol. (2021) 93:2555–6. doi: 10.1002/jmv.26624, PMID: 33095459 Davis HE Assaf GS McCorkell L Wei H Low RJ Re'em Y . Characterizing long COVID in an international cohort: 7 months of symptoms and their impact. EClinicalMedicine. (2021) 38:101019. doi: 10.1016/j.eclinm.2021.101019, PMID: 34308300 Mancini DM Brunjes DL Lala A Trivieri MG Contreras JP Natelson BH. Use of cardiopulmonary stress testing for patients with unexplained dyspnea post-coronavirus disease. JACC Heart Fail. (2021) 9:927–37. doi: 10.1016/j.jchf.2021.10.002, PMID: 34857177 Xie Y Al-Aly Z. Risks and burdens of incident diabetes in long COVID: a cohort study. Lancet Diabetes Endocrinol. (2022) 10:311–21. doi: 10.1016/S2213-8587(22)00044-4, PMID: 35325624 Xie Y Xu E Bowe B Al-Aly Z. Long-term cardiovascular outcomes of COVID-19. Nat Med. (2022) 28:583–90. doi: 10.1038/s41591-022-01689-3, PMID: 35132265 Swank Z Senussi Y Manickas-Hill Z Yu XG Li JZ Alter G . Persistent circulating SARS-CoV-2 spike is associated with post-acute COVID-19 sequelae. Clin Infect Dis. (2023) 76:e487–90. doi: 10.1093/cid/ciac722 Proal AD VanElzakker MB. Long COVID or post-acute sequelae of COVID-19 (PASC): an overview of biological factors that May contribute to persistent symptoms. Front Microbiol. (2021) 12:698169. doi: 10.3389/fmicb.2021.698169, PMID: 34248921 Glynne P Tahmasebi N Gant V Gupta R. Long COVID following mild SARS-CoV-2 infection: characteristic T cell alterations and response to antihistamines. J Invest Med. (2022) 70:61–7. doi: 10.1136/jim-2021-002051, PMID: 34611034 Phetsouphanh C Darley DR Wilson DB Howe A Munier CML Patel SK . Immunological dysfunction persists for 8 months following initial mild-to-moderate SARS-CoV-2 infection. Nat Immunol. (2022) 23:210–6. doi: 10.1038/s41590-021-01113-x, PMID: 35027728 Su Y Yuan D Chen DG Ng RH Wang K Choi J . Multiple early factors anticipate post-acute COVID-19 sequelae. Cells. (2022) 185:881–95.e20. doi: 10.1016/j.cell.2022.01.014, PMID: 35216672 Haffke M Freitag H Rudolf G Seifert M Doehner W Scherbakov N . Endothelial dysfunction and altered endothelial biomarkers in patients with post-COVID-19 syndrome and chronic fatigue syndrome (ME/CFS). J Transl Med. (2022) 20:138. doi: 10.1186/s12967-022-03346-2, PMID: 35317812 Charfeddine S Ibn Hadj Amor H Jdidi J Torjmen S Kraiem S Hammami R . Long COVID 19 syndrome: is it related to microcirculation and endothelial dysfunction? Insights from TUN-EndCOV study. Front Cardiovasc Med. (2021) 8:745758. doi: 10.3389/fcvm.2021.745758, PMID: 34917659 Pretorius E Venter C Laubscher GJ Kotze MJ Oladejo SO Watson LR . Prevalence of symptoms, comorbidities, fibrin amyloid microclots and platelet pathology in individuals with long COVID/post-acute sequelae of COVID-19 (PASC). Cardiovasc Diabetol. (2022) 21:148. doi: 10.1186/s12933-022-01579-5, PMID: 35933347 Małek ŁA Marczak M Miłosz-Wieczorek B Konopka M Braksator W Drygas W . Cardiac involvement in consecutive elite athletes recovered from Covid-19: a magnetic resonance study. J Magn Reson Imaging. (2021) 53:1723–9. doi: 10.1002/jmri.27513, PMID: 33474768 Rajpal S Tong MS Borchers J Zareba KM Obarski TP Simonetti OP . Cardiovascular magnetic resonance findings in competitive athletes recovering from COVID-19 infection. JAMA Cardiol. (2021) 6:116–8. doi: 10.1001/jamacardio.2020.4916, PMID: 32915194 Kociol RD Cooper LT Fang JC Moslehi JJ Pang PS Sabe MA . Recognition and initial Management of Fulminant Myocarditis: a scientific statement from the American Heart Association. Circulation. (2020) 141:e69–92. doi: 10.1161/CIR.0000000000000745, PMID: 31902242 Corrado D Basso C Thiene G. Sudden cardiac death in young people with apparently normal heart. Cardiovasc Res. (2001) 50:399–408. doi: 10.1016/S0008-6363(01)00254-1, PMID: 11334844 Tschöpe C Ammirati E Bozkurt B Caforio ALP Cooper LT Felix SB . Myocarditis and inflammatory cardiomyopathy: current evidence and future directions. Nat Rev Cardiol. (2021) 18:169–93. doi: 10.1038/s41569-020-00435-x, PMID: 33046850 Puntmann VO Carerj ML Wieters I Fahim M Arendt C Hoffmann J . Outcomes of cardiovascular magnetic resonance imaging in patients recently recovered from coronavirus disease 2019 (COVID-19). JAMA Cardiol. (2020) 5:1265–73. doi: 10.1001/jamacardio.2020.3557, PMID: 32730619 Briggs A Vassall A. Count the cost of disability caused by COVID-19. Nature. (2021) 593:502–5. doi: 10.1038/d41586-021-01392-2, PMID: 34040204 Alwan NA . The road to addressing long Covid. Science (New York, NY). (2021) 373:491–3. doi: 10.1126/science.abg7113 Gautier L Cope L Bolstad BM Irizarry RA. Affy--analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. (2004) 20:307–15. doi: 10.1093/bioinformatics/btg405, PMID: 14960456 Ritchie ME Phipson B Wu D Hu Y Law CW Shi W . Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007, PMID: 25605792 Shankar R Leimanis ML Newbury PA Liu K Xing J Nedveck D . Gene expression signatures identify paediatric patients with multiple organ dysfunction who require advanced life support in the intensive care unit. EBioMedicine. (2020) 62:103122. doi: 10.1016/j.ebiom.2020.103122, PMID: 33248372 Ashburner M Ball CA Blake JA Botstein D Butler H Cherry JM . Gene ontology: tool for the unification of biology. Gene Ontol Consort Nat Genet. (2000) 25:25–9. doi: 10.1038/75556, PMID: 10802651 Altermann E Klaenhammer TR. PathwayVoyager: pathway mapping using the Kyoto encyclopedia of genes and genomes (KEGG) database. BMC Genomics. (2005) 6:60. doi: 10.1186/1471-2164-6-60, PMID: 15869710 Yu G Wang LG Han Y He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. (2012) 16:284–7. doi: 10.1089/omi.2011.0118, PMID: 22455463 Walter W Sánchez-Cabo F Ricote M. GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics. (2015) 31:2912–4. doi: 10.1093/bioinformatics/btv300, PMID: 25964631 Szklarczyk D Gable AL Nastou KC Lyon D Kirsch R Pyysalo S . The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. (2021) 49:D605–12. doi: 10.1093/nar/gkaa1074, PMID: 33237311 Shannon P Markiel A Ozier O Baliga NS Wang JT Ramage D . Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303, PMID: 14597658 Chin CH Chen SH Wu HH Ho CW Ko MT Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. (2014) 8:S11. doi: 10.1186/1752-0509-8-S4-S11 Uhlen M Karlsson MJ Zhong W Tebani A Pou C Mikes J . A genome-wide transcriptomic analysis of protein-coding genes in human blood cells. Science (New York, NY). (2019) 366:eaax9198. doi: 10.1126/science.aax9198 Mi H Muruganujan A Huang X Ebert D Mills C Guo X . Protocol update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat Protoc. (2019) 14:703–21. doi: 10.1038/s41596-019-0128-8, PMID: 30804569 Swank Z Senussi Y Manickas-Hill Z Yu XG Li JZ Alter G . Persistent circulating severe acute respiratory syndrome coronavirus 2 spike is associated with post-acute coronavirus disease 2019 sequelae. Clin Infect Dis. (2023) 76:e487–90. doi: 10.1093/cid/ciac722, PMID: 36052466 Yeoh YK Zuo T Lui GC Zhang F Liu Q Li AY . Gut microbiota composition reflects disease severity and dysfunctional immune responses in patients with COVID-19. Gut. (2021) 70:698–706. doi: 10.1136/gutjnl-2020-323020, PMID: 33431578 Liu Q Mak JWY Su Q Yeoh YK Lui GC Ng SSS . Gut microbiota dynamics in a prospective cohort of patients with post-acute COVID-19 syndrome. Gut. (2022) 71:544–52. doi: 10.1136/gutjnl-2021-325989, PMID: 35082169 Spudich S Nath A. Nervous system consequences of COVID-19. Science (New York, NY). (2022) 375:267–9. doi: 10.1126/science.abm2052 Stavileci B Özdemir E Özdemir B Ereren E Cengiz M. De-novo development of fragmented QRS during a six-month follow-up period in patients with COVID-19 disease and its cardiac effects. J Electrocardiol. (2022) 72:44–8. doi: 10.1016/j.jelectrocard.2022.02.012, PMID: 35306293 Grist JT Collier GJ Walters H Kim M Chen M Abu Eid G . Lung abnormalities detected with hyperpolarized (129)Xe MRI in patients with long COVID. Radiology. (2022) 305:709–17. doi: 10.1148/radiol.220069, PMID: 35608443 Galán M Vigón L Fuertes D Murciano-Antón MA Casado-Fernández G Domínguez-Mateos S . Persistent overactive cytotoxic immune response in a Spanish cohort of individuals with long-COVID: identification of diagnostic biomarkers. Front Immunol. (2022) 13:848886. doi: 10.3389/fimmu.2022.848886, PMID: 35401523 Peluso MJ Deeks SG Mustapic M Kapogiannis D Henrich TJ Lu S . SARS-CoV-2 and mitochondrial proteins in neural-derived exosomes of COVID-19. Ann Neurol. (2022) 91:772–81. doi: 10.1002/ana.26350, PMID: 35285072 Afrin LB Weinstock LB Molderings GJ. Covid-19 hyperinflammation and post-Covid-19 illness may be rooted in mast cell activation syndrome. Int J Infect Dis. (2020) 100:327–32. doi: 10.1016/j.ijid.2020.09.016, PMID: 32920235 López CJ . Post-COVID-19 syndrome (PC19S): chronic reactive Endotheliitis and disseminated vascular disease. Acta Med Port. (2020) 33:859. doi: 10.20344/amp.14612, PMID: 33496258 Braile M Marcella S Cristinziano L Galdiero MR Modestino L Ferrara AL . VEGF-A in cardiomyocytes and heart diseases. Int J Mol Sci. (2020) 21:5294. doi: 10.3390/ijms21155294, PMID: 32722551 Taimeh Z Loughran J Birks EJ Bolli R. Vascular endothelial growth factor in heart failure. Nat Rev Cardiol. (2013) 10:519–30. doi: 10.1038/nrcardio.2013.94 Tang J Wang J Kong X Yang J Guo L Zheng F . Vascular endothelial growth factor promotes cardiac stem cell migration via the PI3K/Akt pathway. Exp Cell Res. (2009) 315:3521–31. doi: 10.1016/j.yexcr.2009.09.026, PMID: 19800880 Tao Z Chen B Tan X Zhao Y Wang L Zhu T . Coexpression of VEGF and angiopoietin-1 promotes angiogenesis and cardiomyocyte proliferation reduces apoptosis in porcine myocardial infarction (MI) heart. Proc Natl Acad Sci U S A. (2011) 108:2064–9. doi: 10.1073/pnas.1018925108, PMID: 21245320 Conway EM Collen D Carmeliet P. Molecular mechanisms of blood vessel growth. Cardiovasc Res. (2001) 49:507–21. doi: 10.1016/S0008-6363(00)00281-9 Borgers M Voipio-Pulkki L Izumo S. Apoptosis. Cardiovasc Res. (2000) 45:525–7. doi: 10.1016/S0008-6363(99)00404-6 Zhao T Zhao W Chen Y Liu L Ahokas RA Sun Y. Differential expression of vascular endothelial growth factor isoforms and receptor subtypes in the infarcted heart. Int J Cardiol. (2013) 167:2638–45. doi: 10.1016/j.ijcard.2012.06.127, PMID: 22818386 Kandula V Kosuru R Li H Yan D Zhu Q Lian Q . Forkhead box transcription factor 1: role in the pathogenesis of diabetic cardiomyopathy. Cardiovasc Diabetol. (2016) 15:44. doi: 10.1186/s12933-016-0361-1, PMID: 26956801 Murtaza G Khan AK Rashid R Muneer S Hasan SMF Chen J. FOXO transcriptional factors and long-term living. Oxidative Med Cell Longev. (2017) 2017:1–8. doi: 10.1155/2017/3494289 Wang S Xia P Huang G Zhu P Liu J Ye B . FoxO1-mediated autophagy is required for NK cell development and innate immunity. Nat Commun. (2016) 7:11023. doi: 10.1038/ncomms11023, PMID: 27010363 Xing YQ Li A Yang Y Li XX Zhang LN Guo HC. The regulation of FOXO1 and its role in disease progression. Life Sci. (2018) 193:124–31. doi: 10.1016/j.lfs.2017.11.030 Britton C Poznansky MC Reeves P. Polyfunctionality of the CXCR4/CXCL12 axis in health and disease: implications for therapeutic interventions in cancer and immune-mediated diseases. FASEB J. (2021) 35:e21260. doi: 10.1096/fj.202001273R, PMID: 33715207 Döring Y Pawig L Weber C Noels H. The CXCL12/CXCR4 chemokine ligand/receptor axis in cardiovascular disease. Front Physiol. (2014) 5:212. doi: 10.3389/fphys.2014.00212 Debnath B Xu S Grande F Garofalo A Neamati N. Small molecule inhibitors of CXCR4. Theranostics. (2013) 3:47–75. doi: 10.7150/thno.5376, PMID: 23382786 Isles HM Herman KD Robertson AL Loynes CA Prince LR Elks PM . The CXCL12/CXCR4 signaling Axis retains neutrophils at inflammatory sites in zebrafish. Front Immunol. (2019) 10:1784. doi: 10.3389/fimmu.2019.01784, PMID: 31417560 De Ceuninck F Duguet F Aussy A Laigle L Moingeon P. IFN-α: a key therapeutic target for multiple autoimmune rheumatic diseases. Drug Discov Today. (2021) 26:2465–73. doi: 10.1016/j.drudis.2021.06.010, PMID: 34224903 McHugh J . Hitting the right spot on CXCR4. Nat Rev Rheumatol. (2019) 15:514. doi: 10.1038/s41584-019-0287-6, PMID: 31371800 Adlere I Caspar B Arimont M Dekkers S Visser K Stuijt J . Modulators of CXCR4 and CXCR7/ACKR3 function. Mol Pharmacol. (2019) 96:737–52. doi: 10.1124/mol.119.117663, PMID: 31548340 De Clercq E . Potential clinical applications of the CXCR4 antagonist bicyclam AMD3100. Mini Rev Med Chem. (2005) 5:805–24. doi: 10.2174/1389557054867075, PMID: 16178723 Zhao M Mishra L Deng CX. The role of TGF-β/SMAD4 signaling in cancer. Int J Biol Sci. (2018) 14:111–23. doi: 10.7150/ijbs.23230, PMID: 29483830 Zhang P Hou S Chen J Zhang J Lin F Ju R . Smad4 deficiency in smooth muscle cells initiates the formation of aortic aneurysm. Circ Res. (2016) 118:388–99. doi: 10.1161/CIRCRESAHA.115.308040, PMID: 26699655 Yang J Wang J Zeng Z Qiao L Zhuang L Jiang L . Smad4 is required for the development of cardiac and skeletal muscle in zebrafish. Differentiation. (2016) 92:161–8. doi: 10.1016/j.diff.2016.06.005, PMID: 27477184 Huang L Guo Z Wang F Fu L. KRAS mutation: from undruggable to druggable in cancer. Signal Transduct Target Ther. (2021) 6:386. doi: 10.1038/s41392-021-00780-4, PMID: 34776511 Taniguchi K Karin M. IL-6 and related cytokines as the critical lynchpins between inflammation and cancer. Semin Immunol. (2014) 26:54–74. doi: 10.1016/j.smim.2014.01.001, PMID: 24552665 Yang J Hamid S Cai J Liu Q Xu S Zhang Z. Selenium deficiency-induced thioredoxin suppression and thioredoxin knock down disbalanced insulin responsiveness in chicken cardiomyocytes through PI3K/Akt pathway inhibition. Cell Signal. (2017) 38:192–200. doi: 10.1016/j.cellsig.2017.07.012, PMID: 28734787 Yang J Gong Y Cai J Liu Q Zhang Y Zheng Y . Dysfunction of thioredoxin triggers inflammation through activation of autophagy in chicken cardiomyocytes. Biofactors. (2020) 46:579–90. doi: 10.1002/biof.1625, PMID: 32031748 Luo S Zhang X Xiao X Luo W Yang Z Tang S . Exploring potential biomarkers and molecular mechanisms of ischemic cardiomyopathy and COVID-19 comorbidity based on bioinformatics and systems biology. Int J Mol Sci. (2023) 24:6511. doi: 10.3390/ijms24076511, PMID: 37047484 Abbreviations COVID

coronavirus disease

 SARS-CoV-2

severe acute respiratory syndrome coronavirus 2

 GEO

gene expression omnibus

 PPI

Protein–protein interaction

 DEG

differentially expressed gene

 BP

biological process

 CC

cell component

 MF

molecular function

 VEGFA

vascular endothelial growth factor A

 FOXO1

forkhead box O1

 CXCR4

C-X-C motif chemokine receptor 4

 SMAD4

SMAD family member 4

 KRAS

KRAS proto-oncogene

 TXN

thioredoxin

1 http://r-project.org/

2 http://string-db.org/

3 https://www.proteinatlas.org/

4 http://pantherdb.org/

5 https://www.proteinatlas.org/

6 http://pantherdb.org/