Molecular Mechanisms of Cardiac Injury Associated With Myocardial SARS-CoV-2 Infection

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread around the world. The development of cardiac injury is a common condition in patients with COVID-19, but the pathogenesis remains unclear. The RNA-Seq dataset (GSE150392) comparing expression profiling of mock human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) and SARS-CoV-2-infected hiPSC-CMs was obtained from Gene Expression Omnibus (GEO). We identified 1,554 differentially expressed genes (DEGs) based on GSE150392. Gene set enrichment analysis (GSEA), Gene ontology (GO) analysis, and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that immune-inflammatory responses were activated by SARS-CoV-2, while muscle contraction, cellular respiration, and cell cycle of hiPSC-CMs were inhibited. A total of 15 hub genes were identified according to protein–protein interaction (PPI), among which 11 upregulated genes were mainly involved in cytokine activation related to the excessive inflammatory response. Moreover, we identified potential drugs based on these hub genes. In conclusion, SARS-CoV-2 infection of cardiomyocytes caused a strong defensive response, leading to excessive immune inflammation, cell hypoxia, functional contractility reduction, and apoptosis, ultimately resulting in myocardial injury.

Introduction

Coronavirus disease 2019 (COVID-19) is a viral pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (1). In December 2019, COVID-19 was first reported in Wuhan of China as a severe unknown form of pneumonia. Subsequently, a global pandemic was declared by the WHO in March 2020. At the time of preparing this manuscript, there were more than 5,00,000 fatalities caused by COVID-19.

Acute myocardial damage is the most common described cardiovascular (CV) complication in patients with COVID-19 (2). In multivariable-adjusted models, cardiac injury has been identified as a significant and independent risk factor [hazard ratios (HR) = 4.26] associated with mortality (3). Cytokine storms (caused by acute systemic inflammation) (1, 4), hypoxemia (5), and pathogen-mediated damage (6, 7) were considered as the potential mechanisms responsible for CV complications in COVID-19. However, the exact mechanism of SARS-CoV-2 infection-related myocardial injury remains unclear.

The study of Arun Sharma et al. reported that SARS-CoV-2 directly infected human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) in vitro, inducing contractility depletion, and apoptosis (8). They obtained the expression profiling of three SARS-CoV-2 infected hiPSC-CMs samples and three mock hiPSC-CMs samples by high throughput sequencing and deposit them on the National Center for Biotechnology Information (NCBI, USA) Gene Expression Omnibus (GEO) database (GSE150392). In this study, we performed a detailed bioinformatics analysis of the GSE150392 RNA-seq data to further examine the specific mechanisms of myocardial damage caused by SARS-CoV-2.

Materials and Methods

Data Sources

We obtained the RNA-seq data (GSE150392) from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus GEO database. The GSE150392 dataset consists of six samples divided into two groups. The Cov group included three hiPSC-CMs (human-induced pluripotent stem cell-derived cardiomyocytes) samples pretreated with SARS-COV-2 for 72 h (Cov-1-3: GSM4548303-5), and the Mock group included three hiPSC-CMs samples with a mock treatment without virus (Mock1-3: GSM4548306-8).

Identification of Differentially Expressed Genes

We used R (version 4.0.1, https://www.R-project.org/) package “DESeq2” (version 1.28.1) to determine differentially expressed genes (DEGs) (9) (Padj <0.01), | (log2FoldChange | > 2) between SARS-CoV-2 infected hiPSC-CMs (Cov group) and Mock hiPSC-CMs (Mock group). The R package “ggplot2” was used to graph the Gene expression values boxplot and minus-vs.-add (MA) plot. Padj is the p-value using the Benjamini-Hochberg procedure correction, i.e., false discovery rate (FDR).

Gene Set Enrichment Analysis of All Detected Genes

Gene set enrichment analysis (GSEA, a joint project of UC San Diego and Broad Institute, USA) was performed using the GSEA software (version 4.0.3) implementation in our study for identifying potential hallmarks of SARS-CoV-2 infected hiPSC-CMs (10). The annotated gene sets of “h.all.v7.1.symbols.gmt” were adopted from The Molecular Signatures Database (MSigDB, a joint project of UC San Diego and Broad Institute, USA; http://www.broad.mit.edu/gsea/msigdb/index.jsp). The collapse dataset to gene symbols was “False.” The permutation type was “gene set.” Other parameters were set to default values. Gene sets were considered statistically significant when nominal p-value <0.05, FDR was <0.05, and normalized enrichment score (NES) > 1.5.

Enrichment Analyses of DEGs

Gene ontology (GO) analysis (biological processes, molecular function, and cellular component) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment were analyzed using R package clusterProfiler (version 3.16) (11). Gene ontology analysis (Immune System Process), KEGG pathways functionally grouped networks, and ClinVar human diseases analysis using the CluGO (version 2.5.7) (12) and CluePedia (version 1.5.7) (13) apps of Cytoscape Software (version 3.8) (14). After inputting DEGs in CluGO, select an analysis method in “Ontologies/Pathways” with the following settings: Statistical Test Used was “Enrichment/Depletion (Two-sided hypergeometric test),” Correction Method Used was “Benjamini-Hochberg,” GO Fusion was “false,” and GO Group was “true.” The one with the smallest Padj value was retained among the results of the same GO group. Padj <0.05 was considered as the cut-off criteria.

Protein–Protein Interaction Network

The online tool of a search tool for the retrieval of interacting genes (STRING, Swiss Institute of Bioinformatics, Switzerland; https://string-db.org) (15) was applied to establish a Protein–Protein Interactions (PPIs) of DEGs with the score (median confidence) >0.4. Cytoscape and CytoHubba (version 0.1) (16) were used to visualize the PPI network and identify hub genes (Top 15). GeneMANIA online database (https://genemania.org/) was used to analyze the hub genes (17).

Identifying Drug Candidates

We used the Drug Signatures database DSigDB (University of Colorado at Denver and Health Sciences Center, USA) tool of Enrichr (https://maayanlab.cloud/Enrichr/) to identify drug candidates based on hub genes (18). Drug candidates are ranked by a combined score from highest to lowest.

Results

Severe Acute Respiratory Syndrome Coronavirus 2 Infection Caused a Large Number of Gene Expression Changes in hiPSC-CMs

The boxplot of expression values of all transcripts indicated similar whole transcriptome expression in each sample (Figure 1A) and sample to sample distance heatmap revealed a relatively clear distinction between samples in the two groups (Figure 1B). Based on the criteria of Padj <0.01, and |logFC(fold change)| > 2 (Figure 1C), a total of 1,554 DEGs was identified from GSE150392, including 726 upregulated genes and 828 downregulated genes. Sharma et al. (8) performed a similar DEGs analysis in their previous study, which was not significantly different from the results of this study.

FIGURE 1

Figure 1. Differential expression of data between two sets of samples. (A) Gene expression values boxplot of each sample after normalization using DESeq method; (B) Sample to sample distance heatmap of all samples in Cov and Mock group; (C) MA plot of genes (The blue points represent genes screened based on |fold change| > 2 and a corrected Padj of < 0.01).

Gene Set Enrichment Analysis Showed the SARS-CoV-2 Infection Was Mainly Associated With Immune-Inflammatory Response Activity

To obtain insight into the effect of SARS-CoV-2 on the heart, GSEA was used to map into hallmarks between two groups (Cov vs Mock), 14 significant gene sets were shown (Table 1; Figures 2a–o). Gene set enrichment analysis demonstrated that many inflammation-related gene sets, such as TNFα signaling via NF-κB, interferon-γ response, interferon-α response, inflammatory response, IL6-JAK-STAT3 signaling, IL2-STAT5 signaling, hypoxia, P53 pathway, and apoptosis gene sets were enriched by SARS-CoV-2 infection in hiPSC-CMs. Oxidative phosphorylation, E2F targets, G2M checkpoint, myogenesis, and MYC targets v1 gene sets associated with cell cycle, aerobic respiration, and myogenesis were inhibited in hiPSC-CMs with SARS-CoV-2 infection.

TABLE 1

Table 1. The most significant gene sets of Cov vs. Mock in gene set enrichment analysis (GSEA).

FIGURE 2

Figure 2. Gene set enrichment analysis (GSEA) of SARS-CoV-2 infection in hiPSC-CMs. (a–i) Nine representative gene sets enriched in hiPSC-CMs with SARS-CoV-2 infection; (j–n) Five representative gene sets inhibited in hiPSC-CMs with SARS-CoV-2 infection; (o) Ranked gene list correlation profile. The red color indicates the “COV” group, the blue color indicates the “MOCK” group, the black line indicates the genes in the GENE SET, and the green line consists of the enrichment fraction of each gene.

Gene Ontology Term Enrichment Analyses and ClinVar Human Diseases Analysis of DEGs

Gene ontology analysis of DEGs was divided into four functional groups, including biological processes (BP), cell composition (CC), molecular function (MF), and Immune System Process (ISP). The top five results of BP, CC, and MF are shown in Figure 3.

FIGURE 3

Figure 3. Gene ontology (GO) analysis (BP, CC, and MF) of differentially expressed genes (DEGs). (A) GO terms of upregulated DEGs. (B) GO terms of downregulated DEGs.

In the BP group, the upregulated genes were mainly enriched in response to the virus, response to lipopolysaccharide, response to molecule of bacterial origin, leukocyte cell–cell adhesion, and defense response to the virus. The downregulated genes were mainly concentrated in the muscle system process, muscle contraction, oxidative phosphorylation, striated muscle contraction, and respiratory electron transport chain.

In the CC group, the upregulated genes were mainly enriched in receptor complex. The downregulated genes were mainly concentrated in myofibril, contractile fiber, sarcomere, I band, and the inner mitochondrial membrane protein complex.

In the MF group, the upregulated genes were mainly involved in terms of cytokine receptor binding, cytokine activity, DNA-binding transcription activator activity (RNA polymerase II-specific), DNA-binding transcription activator activity, and chemokine activity. The downregulated genes were mainly enriched in actin binding, NADH dehydrogenase activity, NADH dehydrogenase (ubiquinone) activity, NADH dehydrogenase (quinone) activity, and structural constituent of muscle.

The top five ISP results were shown in Table 2. Type I interferon signaling pathway, regulation of adaptive immune response, neutrophil migration, regulation of type 2 immune response, and CD4-positive, α-β T cell cytokine production were enriched by upregulated genes. Positive regulation of megakaryocyte differentiation, regulation of megakaryocyte differentiation, and thymus development were enriched by downregulated genes.

TABLE 2

Table 2. Gene ontology (GO) (immune system process) analysis of differentially expressed genes (DEGs).

ClinVar human diseases analysis indicates that DEGs caused by SARS-CoV-2 infection are significantly involved in terms of myocardial infarction 1, non-syndromic hearing loss and deafness, cardiomyopathy, limb-girdle muscular dystrophy, mitochondrial complex IV deficiency, myofibrillar myopathy, and primary ciliary dyskinesia (Table 3).

TABLE 3

Table 3. ClinVar human diseases analysis of DEGs.

KEGG Enrichment Analysis of DEGs

The signaling pathways of DEGs were shown in Figure 4. The data were imported into Cytoscape to calculate the topological characteristics of the network and determine each node. Cytokine–cytokine receptor interaction, TNF signaling pathway, influenza A, NF-κB signaling pathway, osteoclast differentiation, measles, viral–protein interaction with cytokine and cytokine receptor, IL-17 signaling pathway, Toll-like receptor signaling pathway, and rheumatoid arthritis were significantly enriched by upregulated DEGs (Figure 4A). Thermogenesis, Parkinson's disease, oxidative phosphorylation, Huntington disease, cardiac muscle contraction, non-alcoholic fatty liver disease (NAFLD), retrograde endocannabinoid signaling, adrenergic signaling in cardiomyocytes, dilated cardiomyopathy (DCM), and hypertrophic cardiomyopathy (HCM) were significantly enriched by downregulated DEGs (Figure 4B). The top seven p-value KEGG pathways with their target genes were shown (Figure 4C).

FIGURE 4

Figure 4. The Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis of DEGs. (A) KEGG pathways of upregulated DEGs. (B) KEGG pathways of downregulated DEGs. (C) Representative pathways to genes network of upregulated DEGs validated genes (in red) targeted by KEGG pathways.

PPI Network Analysis and Hub Genes Recognition

The top 15 genes with the highest interaction degrees were identified, including 11 upregulated genes and four downregulated genes (Table 4). Then those genes were used for constructing the PPI network (Figure 5A). A total of 11 upregulated hub genes showed the complex PPI network with the co-expression of 66.42%, co-localization of 9.39%, genetic interactions of 8.72%, predicted of 8.6%, shared protein domains of 3.99%, pathway of 1.8%, and physical interactions of 1.09% (Figure 5C). Cytokine activity, chemokine activity, inflammatory response, leukocyte chemotaxis, and lipopolysaccharide-mediated signaling pathway were identified as the main function of those genes. Furthermore, four downregulated hub genes showed the complex PPI network with the physical interactions of 67.64%, co-expression of 13.50%, predicted of 6.35%, co-localization of 6.17%, the pathway of 4.35%, genetic interactions of 1.4%, and shared protein domains of 0.59% (Figure 5B). Those genes were related to mitosis, nuclear division, metaphase/anaphase transition of the cell cycle, and regulation of ubiquitin–protein ligase activity.

TABLE 4

Table 4. Top 15 hub genes with the highest interaction degrees in protein–protein interaction (PPI) network analysis.

FIGURE 5

Figure 5. Protein–protein interaction (PPI) network of hub genes. (A) The PPI network of the top 15 hub genes created by STRING; (B) PP1 networks and function analyses of the four downregulated hub genes; (C) PPI networks and function analyses of the 11 upregulated hub genes. Inner circles represent the input genes and outer circles correspond to GeneMANIA proposed hub genes, the size of the circles indicates the correlation with the input genes.

Identification of Candidate Drugs

The study further involved 15 hub genes that were used to screen potential drugs for the treatment of SARS-CoV-2 infection-related myocardial injury. According to the combined score, the top 10 candidate drugs were generated (Table 5). Acetovanillone, antimycin A, proscillaridin, chromium (III) chloride, and hyperforin were identified by 11 upregulated hub genes. Irinotecan hydrochloride, troglitazone, 67526-95-8 (thapsigargin), genistein, and testosterone were identified by four downregulated hub genes.

TABLE 5

Table 5. Drug candidates combined with hub genes.

Discussion

Severe acute respiratory syndrome coronavirus 2, SARS-CoV, and MERS-CoV (Middle East respiratory syndrome coronavirus) are enveloped positive RNA viruses, belong to the Coronaviridae family β genera (19, 20). Severe acute respiratory syndrome coronavirus 2 and SARS-CoV utilize the same receptor angiotensin-converting enzyme 2 (ACE2) for invading human bodies. Those coronaviruses spread widely among people and lead to potentially fatal diseases (21, 22). Coronavirus disease 2019 caused by SARS-CoV-2 has been spreading in 219 countries/areas/territories, with over 50,000,000 confirmed cases (WHO statistics as of November 9, 2020) (23).

The cardiac injury in COVID-19 was generally defined as the elevation of cardiac-specific biomarkers (e.g., troponin concentration above the 99th percentile upper reference limit) (1, 2). It occurs in roughly 8–13% (24, 25) in confirmed cases and as high as 23–44% in severe patients (26–30). Studies have shown that mortality was markedly higher in patients with cardiac injury than in patients without cardiac injury (51.2–59.6 vs. 4.5–8.9%) (3, 31). Besides, patients with underlying cardiovascular disease (CVD) are more likely to develop acute myocardial injury (29, 31). Concurrent occurrence of underlying CVD and myocardial injury causes significantly higher mortality than patients with underlying CVD but without cardiac injury (69.4 vs. 13.3%) (31), suggesting that myocardial injury played a greater role in the fatal outcome of COVID-19 than the presence of underlying CVD itself.

Angiotensin-converting enzyme 2 is an important target for SARS-CoV and SARS-CoV-2 infection, and the heart is one of the organs which express ACE2 (32). During the SARS epidemic, the SARS-CoV viral RNA was detected on autopsied human heart samples (33). Similarly, Guido Tavazzi et al. reported the first case of myocardial localization of SARS-CoV-2 (34). Furthermore, hiPSC-CMs are susceptible to SARS-CoV-2 infection in vitro (8). Those studies all suggested that SARS-CoV-2 may directly infect the heart in the human body through ACE2. However, the pathogenesis of SARS-CoV-2 infection-related acute myocardial injury is still unknown. Previous studies suggested that pathogen-mediated direct myocardial injury, acute systemic inflammatory response (1, 4), and low blood oxygen levels (7) may be the most critical causes of myocardial damage in patients with COVID-19. In this study, we described several specific mechanisms of cardiac injury caused by SARS-CoV-2 direct infection.

We used R software and bioinformatics to deeply analyze the RNA-Seq dataset GSE150293, which compared the gene expression between hiPSC-CMs and SARS-CoV-2-infected hiPSC-CMs samples. The results identified 1554 DEGs, including 726 upregulated genes and 828 downregulated genes. In order to reduce the bias, we first performed a GSEA analysis of all detected genes. Gene set enrichment analysis revealed that SARS-CoV-2 infection activated the inflammatory response, and inhibited oxidative phosphorylation and myogenesis, resulting in hypoxia and apoptosis in cardiomyocytes. Subsequently, GO analysis of DEGs revealed specific mechanisms of SARS-CoV-2 impairing myocardial function. Gene ontology analysis of upregulated DEGs suggested several specific mechanisms of inflammatory response activation, including promotion of receptor complexes expression, enhancement of cytokine and chemokine activity, promotion of cytokine binding to receptors, and boosted leukocyte adhesion. Gene ontology analysis of downregulated DEGs showed muscle system process, muscle contraction, oxidative phosphorylation, striated muscle contraction, and respiratory electron transport chain were significantly inhibited by SARS-CoV-2 infection. The expression of myofibril, contractile fiber, sarcomere, I band, and inner mitochondrial membrane protein complex was also suppressed. In the MF group, SARS-CoV-2 inhibits actin binding, NADH dehydrogenase (such as ubiquinone and quinone) activity, and structural constituent of muscle. Kyoto encyclopedia of genes and genomes pathway analysis showed thermogenesis, oxidative phosphorylation, cardiac muscle contraction, retrograde endocannabinoid signaling, and adrenergic signaling in cardiomyocytes were enriched with downregulated DEGs. Those results reveal specific mechanisms by which SARS-CoV-2 infection causes respiration dysfunction and muscle contraction disorders in hiPSC-CMs.

Previous studies have shown that SARS-CoV-2 infection activates multiple immune responses (35, 36), and to clarify the specific immune response processes involved, we performed ISP enrichment analysis. Similar to previous studies, ISP analysis of upregulated DEGs indicated that adaptive and type-2 immune responses were activated by SARS-CoV-2. In addition, neutrophil migration and CD4-positive α-β T cell cytokine production were promoted. Similarly, GSEA and KEGG pathway analysis showed that immune response correlative signal pathways were activated by SARS-CoV-2, including TNFα, IL6-JAK-STAT3, IL2-STAT5, NF-κB, IL17, and Toll-like receptor signaling pathway. Consistent with the inflammatory cytokines detected in the blood of patients with COVID-19 (1, 37), hiPSC-CMs produce and activate various cytokines, including TNFα, chemokines, interleukins (IL), and interferons after SARS-CoV-2 infection. The massive cytokine release reflects an excessive immune defense, which may cause harmful heart damage. These findings complement our understanding of the human immune response mechanism triggered by SARS-CoV-2 (35) and contribute to the treatment of patients with COVID-19. Similar to our study, Omar Pacha et al. suggested IL17 is immunologically plausible as a target to prevent ARDS in COVID-19 (38).

Furthermore, ISP analysis of downregulated DEGs showed positive regulation of megakaryocyte differentiation, and thymus development were repressed by SARS-CoV-2 infection. These mechanisms may be responsible for the significant reduction of platelet count (39) and peripheral blood T cells (37) in patients with severe COVID-19.

After removing confounding results of genetic disorders, ClinVar human diseases analysis showed SARS-CoV-2 infection may be associated with myocardial infarction (40, 41) and cardiomyopathy (42), which are similar to the clinically reported cardiac complications in patients with COVID-19. Cardiac injury often results in chronic heart failure due to the loss of cardiomyocytes (43). Therefore, we speculate that some patients with COVID-19 combined with myocardial damage may have long-term cardiac insufficiency after recovery.

The top genes with the highest degree of interaction in the PPI network are considered to be hub genes. The identification of hub genes may be critical for the diagnosis and treatment of myocardial injury in patients with COVID-19. In this study, functional analysis of 11 upregulated hub genes showed association with cytokine activity, chemokine activity, inflammatory response, leukocyte chemotaxis, and lipopolysaccharide-mediated signaling pathway. These results were highly similar to the GO analysis of all upregulated DEGs, suggesting that these 11 upregulated hub genes play an irreplaceable role in the SARS-CoV-2-induced myocardial injury process.

It is generally accepted that the regeneration of cardiomyocytes was highest in human infants (44, 45) and declines to relative stagnation in adults (46, 47). Restoring the mitotic activity of cardiomyocytes is considered to be the ultimate solution for the treatment of irreversible heart failure (48). The up-regulation of positive cell cycle regulators is one of the endogenous mechanisms for cardiomyocyte proliferation and regeneration (49). Adenoviral transfection of E2F2 expression in mouse cardiomyocytes leads to the proliferation of adult cardiomyocytes (50). Overexpression of G2/M transition promoter cyclin A2 in ischemic porcine hearts promotes recovery after ischemic injury and induced cytokinesis (51). In this study, GSEA showed that cell cycle-related gene sets such as E2F targets, G2/M checkpoint, myogenesis, and MYC targets v1 were inhibited by SARS-CoV-2 infection. The analysis of the PPI network showed that SARS-CoV-2 inhibits the regenerative potency of hiPSC-CMs by inhibiting the metaphase/anaphase transition of the mitotic cycle, suggesting that the downregulated four hub genes (CDK1/UBE2C/CDC20/AURKB) may be important genes for myogenesis. More research is needed to explain the role of SARS-CoV-2 infection-related mitotic pathways modification in cardiomyocytes, which may contribute to cardiac function recovery in patients with COVID-19.

This study identified 10 drug molecules by hub genes, including acetovanillone, antimycin A, proscillaridin chromium (III) chloride, hyperforin, irinotecan hydrochloride, troglitazone, thapsigargin, genistein, and testosterone. Among these 10 drugs, the cytotoxic drugs chromium (III) chloride and irinotecan hydrochloride should not be considered as candidates. Antimycin A (52), proscillaridin (53, 54), troglitazone (55, 56), thapsigargin (57, 58), genistein (59, 60), and testosterone (61) were considered as potential drugs for COVID-19 in previous studies. Acetovanillone and hyperforin have not been studied for the treatment of COVID-19. The therapeutic effect of these eight candidate drugs on myocardial injury in COVID-19 needs further study.

Overall, our study suggested that SARS-CoV-2 infection can induce a strong immune-inflammatory response, reduce contractility, increase hypoxia, and induce apoptosis in cardiomyocytes. Excessive inflammatory response, respiration dysfunction, and contraction disorders work together to cause cardiac injury. This study supports that direct infection of cardiomyocytes by SARS-CoV-2 is the primary cause of myocardial injury. It should be noted that the occurrence of myocardial injury may be a sign of SARS-CoV-2 entering the circulation, which also explains the acute kidney injury is more common among patients with cardiac injury (3).

The current study has several limitations. The study included only one available data packet of six hiPSC-CMs samples in vitro, which could not fully reflect the effect of SARS-CoV-2 on the heart in vivo. Only bioinformatics analysis was conducted in this study, further clinical reports or in vitro studies are needed to support our conclusions.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150392.

Author Contributions

XL, LL, and LZ designed analyses and drafted the manuscript. XL and LL performed the bioinformatics analysis. LZ reviewed the manuscript and provided funding support. All authors contributed to the final manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (Grant/Award Numbers: 81970723).

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

1. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. (2020) 395:497–506. doi: 10.1016/S0140-6736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Bansal M. Cardiovascular disease and COVID-19. Diabetes Metab Syndr. (2020) 14:247–50. doi: 10.1016/j.dsx.2020.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Shi S., Qin M., Shen B., Cai Y., Liu T., Yang F., et al. (2020). Association of cardiac injury with mortality in hospitalized patients with COVID-19 in Wuhan, China. JAMA Cardiol. 5:802–10. doi: 10.1001/jamacardio.2020.0950

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. (2020) 395:1054–62. doi: 10.1016/S0140-6736(20)30566-3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Cameron MJ, Bermejo-Martin JF, Danesh A, Muller MP, Kelvin DJ. Human immunopathogenesis of severe acute respiratory syndrome (SARS). Virus Res. (2008) 133:13–9. doi: 10.1016/j.virusres.2007.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Li B, Yang J, Zhao F, Zhi L, Wang X, Liu L, et al. Prevalence and impact of cardiovascular metabolic diseases on COVID-19 in China. Clin Res Cardiol. (2020) 109:531–8. doi: 10.1007/s00392-020-01626-9

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Xiong TY, Redwood S, Prendergast B, Chen M. Coronaviruses and the cardiovascular system: acute and long-term implications. Eur Heart J. (2020) 41:1798–800. doi: 10.1093/eurheartj/ehaa231

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sharma A, Garcia G Jr, Wang Y, Plummer JT, Morizono K, Arumugaswami V, et al. Human iPSC-derived cardiomyocytes are susceptible to SARS-CoV-2 infection. Cell Rep Med. (2020) 1:100052. doi: 10.1016/j.xcrm.2020.100052

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text

10. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

11. 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

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Bindea G, Galon J, Mlecnik B. CluePedia Cytoscape plugin: pathway insights using integrated experimental and in silico data. Bioinformatics. (2013) 29:661–3. doi: 10.1093/bioinformatics/btt019

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. (2019) 47:D607–13. doi: 10.1093/nar/gky1131

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Chin CH, Chen SH, Wu HH, Ho CW, Ko MT, Lin CY. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. (2014) 8(Suppl 4):S11. doi: 10.1186/1752-0509-8-S4-S11

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. (2010) 38:W214–20. doi: 10.1093/nar/gkq537

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Xie Z, Bailey A, Kuleshov MV, Clarke DJB, Evangelista JE, Jenkins SL, et al. Gene set knowledge discovery with enrichr. Curr Protoc. (2021) 1:e90. doi: 10.1002/cpz1.90

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Li F. Structure, function, and evolution of coronavirus spike proteins. Annu Rev Virol. (2016) 3:237–61. doi: 10.1146/annurev-virology-110615-042301

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Alsaadi EAJ, Jones IM. Membrane binding proteins of coronaviruses. Future Virol. (2019) 14:275–86. doi: 10.2217/fvl-2018-0144

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Al-Tawfiq JA, Auwaerter PG. Healthcare-associated infections: the hallmark of Middle East respiratory syndrome coronavirus with review of the literature. J Hosp Infect. (2019) 101:20–9. doi: 10.1016/j.jhin.2018.05.021

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Palacios Cruz M, Santos E, Velazquez Cervantes MA, Leon Juarez M. COVID-19, a worldwide public health emergency. Rev Clin Esp. (2020). 221:55–61. doi: 10.1016/j.rce.2020.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

23. WHO. WHO Coronavirus Disease (COVID-19) Dashboard. Geneva: World Health Organization (2020). Available online at: https://covid19.who.int/

24. Lippi G, Plebani M. Laboratory abnormalities in patients with COVID-2019 infection. Clin Chem Lab Med. (2020) 58:1131–4. doi: 10.1515/cclm-2020-0198

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Yang X, Yu Y, Xu J, Shu H, Xia J, Liu H, et al. Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: a single-centered, retrospective, observational study. Lancet Respir Med. (2020) 8:475–81. doi: 10.1016/S2213-2600(20)30079-5

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chen T, Wu D, Chen H, Yan W, Yang D, Chen G, et al. Clinical characteristics of 113 deceased patients with coronavirus disease 2019: retrospective study. BMJ. (2020) 368:m1091. doi: 10.1136/bmj.m1091

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Deng Q, Hu B, Zhang Y, Wang H, Zhou X, Hu W, et al. Suspected myocardial injury in patients with COVID-19: evidence from front-line clinical observation in Wuhan, China. Int J Cardiol. (2020) 311:116–21. doi: 10.1016/j.ijcard

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. (2020) 382:1708–20. doi: 10.1056/NEJMoa2002032

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Lala A, Johnson KW, Januzzi JL, Russak AJ, Paranjpe I, Richter F, et al. Prevalence and impact of myocardial injury in patients hospitalized with COVID-19 infection. J Am Coll Cardiol. (2020). 76:533–46. doi: 10.1016/j.jacc.2020.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Wang D, Hu B, Hu C, Zhu F, Liu X, Zhang J, et al. Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan, China. JAMA. (2020) 323:1061–9. doi: 10.1001/jama.2020.1585

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Guo T, Fan Y, Chen M, Wu X, Zhang L, He T, et al. Cardiovascular implications of fatal outcomes of patients with coronavirus disease 2019 (COVID-19). JAMA Cardiol. (2020) 5:811–8. doi: 10.1001/jamacardio.2020.1017

PubMed Abstract | CrossRef Full Text | Google Scholar

32. South AM, Diz DI, Chappell MC. COVID-19, ACE2, and the cardiovascular consequences. Am J Physiol Heart Circ Physiol. (2020) 318:H1084–90. doi: 10.1152/ajpheart.00217.2020

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Oudit GY, Kassiri Z, Jiang C, Liu PP, Poutanen SM, Penninger JM, et al. SARS-coronavirus modulation of myocardial ACE2 expression and inflammation in patients with SARS. Eur J Clin Invest. (2009) 39:618–25. doi: 10.1111/j.1365-2362.2009.02153.x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Tavazzi G, Pellegrini C, Maurelli M, Belliato M, Sciutti F, Bottazzi A, et al. Myocardial localization of coronavirus in COVID-19 cardiogenic shock. Eur J Heart Fail. (2020) 22:911–5. doi: 10.1002/ejhf.1828

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Catanzaro M, Fagiani F, Racchi M, Corsini E, Govoni S, Lanni C. Immune response in COVID-19: addressing a pharmacological challenge by targeting pathways triggered by SARS-CoV-2. Signal Transduct Target Ther. (2020) 5:84. doi: 10.1038/s41392-020-0191-1

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Wu F, Zhao S, Yu B, Chen YM, Wang W, Song ZG, et al. A new coronavirus associated with human respiratory disease in China. Nature. (2020) 579:265–9. doi: 10.1038/s41586-020-2008-3

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Qin C, Zhou L, Hu Z, Zhang S, Yang S, Tao Y, et al. Dysregulation of immune response in patients with COVID-19 in Wuhan, China. Clin Infect Dis. (2020) 71:762–8. doi: 10.1093/cid/ciaa248

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Pacha O, Sallman MA, Evans SE. COVID-19: a case for inhibiting IL-17? Nat Rev Immunol. (2020) 20:345–6. doi: 10.1038/s41577-020-0328-z

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Lippi G, Plebani M, Henry BM. Thrombocytopenia is associated with severe coronavirus disease 2019 (COVID-19) infections: a meta-analysis. Clin Chim Acta. (2020) 506:145–8. doi: 10.1016/j.cca.2020.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

40. De Rosa S, Spaccarotella C, Basso C, Calabro MP, Curcio A, Filardi PP, et al. Reduction of hospitalizations for myocardial infarction in Italy in the COVID-19 era. Eur Heart J. (2020) 41:2083–8. doi: 10.1093/eurheartj/ehaa409

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Mahmud E, Dauerman HL, Welt FGP, Messenger JC, Rao SV, Grines C, et al. Management of acute myocardial infarction during the COVID-19 pandemic. J Am Coll Cardiol. (2020). 96:336–45. doi: 10.1016/j.jacc.2020.04.039

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Arentz M, Yim E, Klaff L, Lokhandwala S, Riedo FX, Chong M, et al. Characteristics and outcomes of 21 critically ill patients with COVID-19 in Washington State. JAMA. (2020) 323:1612–4. doi: 10.1001/jama.2020.4326

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Bergmann O, Bhardwaj RD, Bernard S, Zdunek S, Barnabe-Heider F, Walsh S, et al. Evidence for cardiomyocyte renewal in humans. Science. (2009) 324:98–102. doi: 10.1126/science.1164680

PubMed Abstract | CrossRef Full Text | Google Scholar

44. MacMahon HE. Hyperplasia and regeneration of the myocardium in infants and in children. Am J Pathol. (1937) 13:845–54.

PubMed Abstract | Google Scholar

45. Zerbini C, Weinberg DS, Perez-Atayde AR. DNA ploidy analysis of myocardial hyperplasia. Hum Pathol. (1992) 23:1427–30. doi: 10.1016/0046-8177(92)90064-a

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Mollova M, Bersell K, Walsh S, Savla J, Das LT, Park SY, et al. Cardiomyocyte proliferation contributes to heart growth in young humans. Proc Natl Acad Sci USA. (2013) 110:1446–51. doi: 10.1073/pnas.1214608110

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Bergmann O, Zdunek S, Felker A, Salehpour M, Alkass K, Bernard S, et al. Dynamics of cell generation and turnover in the human heart. Cell. (2015) 161:1566–75. doi: 10.1016/j.cell.2015.05.026

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Laflamme MA, Murry CE. Heart regeneration. Nature. (2011) 473:326–35. doi: 10.1038/nature10147

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Yester JW, Kuhn B. Mechanisms of cardiomyocyte proliferation and differentiation in development and regeneration. Curr Cardiol Rep. (2017) 19:13. doi: 10.1007/s11886-017-0826-1

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Ebelt H, Zhang Y, Kampke A, Xu J, Schlitt A, Buerke M, et al. E2F2 expression induces proliferation of terminally differentiated cardiomyocytes in vivo. Cardiovasc Res. (2008) 80:219–26. doi: 10.1093/cvr/cvn194

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Shapiro SD, Ranjan AK, Kawase Y, Cheng RK, Kara RJ, Bhattacharya R, et al. Cyclin A2 induces cardiac regeneration after myocardial infarction through cytokinesis of adult cardiomyocytes. Sci Transl Med. (2014) 6, 224ra227. doi: 10.1126/scitranslmed.3007668

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Yang CW, Peng TT, Hsu HY, Lee YZ, Wu SH, Lin WH, et al. Repurposing old drugs as antiviral agents for coronaviruses. Biomed J. (2020) 43:368–74. doi: 10.1016/j.bj.2020.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Aishwarya S, Gunasekaran K, Margret AA. Computational gene expression profiling in the exploration of biomarkers, non-coding functional RNAs and drug perturbagens for COVID-19. J Biomol Struct Dyn. (2020) 2020:1–16. doi: 10.1080/07391102.2020.1850360

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Wang Z, Yang L. Chinese herbal medicine: fighting SARS-CoV-2 infection on all fronts. J Ethnopharmacol. (2021) 270:113869. doi: 10.1016/j.jep.2021.113869

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Chandra A, Gurjar V, Ahmed MZ, Alqahtani AS, Qamar I, Singh N. Exploring potential inhibitor of SARS-CoV2 replicase from FDA approved drugs using insilico drug discovery methods. J Biomol Struct Dyn. (2021) 1–8. doi: 10.1080/07391102.2020.1871416

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Shahjaman M, Rezanur Rahman M, Rabiul Auwul M. A network-based systems biology approach for identification of shared Gene signatures between male and female in COVID-19 datasets. Inform Med Unlocked. (2021) 25:100702. doi: 10.1016/j.imu.2021.100702

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Al-Beltagi S, Preda CA, Goulding LV, James J, Pu J, Skinner P, et al. Thapsigargin is a broad-spectrum inhibitor of major human respiratory viruses: coronavirus, respiratory syncytial virus and influenza A virus. Viruses. (2021) 13:234. doi: 10.3390/v13020234

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Shaban MS, Muller C, Mayr-Buro C, Weiser H, Meier-Soelch J, Albert BV, et al. Multi-level inhibition of coronavirus replication by chemical ER stress. Nat Commun. (2021) 12:5536. doi: 10.1038/s41467-021-25551-1

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Elfiky AA. Natural products may interfere with SARS-CoV-2 attachment to the host cell. J Biomol Struct Dyn. (2021) 39:3194–203. doi: 10.1080/07391102.2020.1761881

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Gour A, Manhas D, Bag S, Gorain B, Nandi U. Flavonoids as potential phytotherapeutics to combat cytokine storm in SARS-CoV-2. Phytother Res. (2021) 35:4258–83. doi: 10.1002/ptr.7092

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Al-Lami RA, Urban RJ, Volpi E, Algburi AMA, Baillargeon J. Sex hormones and novel corona virus infectious disease (COVID-19). Mayo Clin Proc. (2020) 95:1710–4. doi: 10.1016/j.mayocp.2020.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar