- 1Department of Cardiology, The Affiliated Hospital of Qingdao University, Qingdao, China
- 2Department of Genetics and Cell Biology, Basic Medical College, Qingdao University, Qingdao, China
- 3Department of Dermatology, The Affiliated Haici Hospital of Qingdao University, Qingdao, China
- 4Department of Cardiology, The Affiliated Cardiovascular Hospital of Qingdao University, Qingdao, China
Background: Increasing evidence indicating that coronavirus disease 2019 (COVID-19) increased the incidence and related risks of pericarditis and whether COVID-19 vaccine is related to pericarditis has triggered research and discussion. However, mechanisms behind the link between COVID-19 and pericarditis are still unknown. The objective of this study was to further elucidate the molecular mechanisms of COVID-19 with pericarditis at the gene level using bioinformatics analysis.
Methods: Genes associated with COVID-19 and pericarditis were collected from databases using limited screening criteria and intersected to identify the common genes of COVID-19 and pericarditis. Subsequently, gene ontology, pathway enrichment, protein–protein interaction, and immune infiltration analyses were conducted. Finally, TF–gene, gene–miRNA, gene–disease, protein–chemical, and protein–drug interaction networks were constructed based on hub gene identification.
Results: A total of 313 common genes were selected, and enrichment analyses were performed to determine their biological functions and signaling pathways. Eight hub genes (IL-1β, CD8A, IL-10, CD4, IL-6, TLR4, CCL2, and PTPRC) were identified using the protein–protein interaction network, and immune infiltration analysis was then carried out to examine the functional relationship between the eight hub genes and immune cells as well as changes in immune cells in disease. Transcription factors, miRNAs, diseases, chemicals, and drugs with high correlation with hub genes were predicted using bioinformatics analysis.
Conclusions: This study revealed a common gene interaction network between COVID-19 and pericarditis. The screened functional pathways, hub genes, potential compounds, and drugs provided new insights for further research on COVID-19 associated with pericarditis.
1 Introduction
COVID-19 is an atypical respiratory disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has triggered a global pandemic and caused significant loss of life and property (1). Globally, as of 21 January, there were 774,395,593 confirmed cases of COVID-19, including 7,023,271 deaths, reported by the World Health Organization (WHO) (https://covid19.who.int/). Although more than 80% of patients with COVID-19 present with asymptomatic infection or mild to moderate self-resolving symptoms, more than 15% of patients still develop into severe cases, manifested as severe pneumonia or acute respiratory distress syndrome, and even multiple organ failure (2, 3). In addition to respiratory diseases, cardiovascular complications have gradually become a major threat for patients with COVID-19. Pericarditis is the most common pericardial disease worldwide; the pericardium provides fixation and physical protection for the heart, such as slowing down the impact of heart contraction on the surrounding blood vessels and preventing the spread of pulmonary and thoracic infections (4). The etiology of pericarditis may be infectious (bacterial or viral) or noninfectious (systemic inflammatory disease or post-cardiac injury syndrome) (5). Viral infection is an important cause of pericarditis and studies have shown that it can be an early complication of COVID-19. Notably, the incidence rate of pericarditis has increased by at least 15 times after SARS-CoV-2 infection than before COVID-19, and estimates of excess cases associated with vaccination also indicate a burden associated with pericarditis (6, 7). Over 13 billion doses of COVID-19 vaccines have been administered; furthermore, several passive surveillance systems have indicated that the risk of pericarditis increased after COVID-19 vaccination, especially in young men, but authoritative research claimed that the incidence was rare (8). Considering benefits and risks, vaccination should be firmly supported, but strengthening the surveillance of adverse events following vaccination and continuing to study the mechanistic relationship between COVID-19 and pericarditis are still essential.
Increasing evidence suggests that immune responses and potential immune markers may be associated with COVID-19 severity. Differences in innate immune system components lead to heterogeneity in the COVID-19 disease spectrum (9). An imbalanced immune response during viral invasion is an important immunopathological mechanism in severe diseases (10). After SARS-CoV-2 infection, immune effector cells release a large number of pro-inflammatory cytokines, triggering a cytokine storm that causes important immunopathological events, such as ARDS and multiple organ failure (11). In recent years, immune checkpoints have led to breakthroughs and progress in cancer treatment; however, their application is still limited due to immune-related adverse events, such as cardiotoxicity. The onset and progression of pericarditis in the cardiotoxicity brought on by immunotherapy are intimately tied to the unrestricted regulation of the immune system (12). The predisposing factors and pathogenesis of pericarditis remain unclear, and may be related to viral infections or autoimmune-inflammatory diseases. Under the influence of exogenous triggers, infections may lead to an autoimmune response in susceptible hosts by activating innate immunity (13, 14). Therefore, gaining a comprehensive and in-depth understanding of interactions between viruses and the human immune system is necessary. Furthermore, studying the impact and mechanism of the clinical outcomes of COVID-19 and pericarditis is crucial to promote the research and development of vaccines with reduced side effects.
In recent years, with the rapid development of high-throughput biotechnology, the use of gene interaction networks in bioinformatics research has become increasingly convenient. The construction of gene interaction network not only helps to further understand various biological processes from the perspective of systems, but also can be widely applied to explore the pathogenesis of diseases. However, massive data cannot be verified one by one to explain the mechanism, so gene enrichment studies are needed to classify differential genes, so as to filter redundant data and screen out more valuable functional information (15). Based on data sources and algorithms, methods for gene functional enrichment analysis can be roughly divided into four categories: over-representative analysis (ORA), function set scoring (FCS), pathway topology (PT), and network topology (NT). Enrichment analyses commonly used include GO enrichment analysis, KEGG enrichment analysis and gene set enrichment analysis (GSEA). GSEA consists of three key elements: calculating enrichment scores, evaluating significance, and adjusting for multiple hypothesis tests. Weighted Kolmogorov Smirnov (WKS) test was used in GSEA to obtain the statistical value of the functional set of the tested gene, and there are other statistical algorithms, such as χ2-test, Mean test, Median test, Wilcoxon rank sum test, etc. (16, 17). The schematic overview for GSEA can be found in Supplementary Figure 1 (https://www.gsea-msigdb.org/gsea/) (16).
In this study, we employed a wide range of bioinformatics techniques to identify the common genes of COVID-19 and pericarditis and analyzed their enrichment pathways and functions. Protein–protein interaction (PPI) networks were constructed to identify hub genes and further analyze the interaction networks of transcription factors (TFs), microRNAs (miRNAs), chemicals, and drugs. The immune response can serve as a resection point for studying the common pathogenesis of comorbidities. To uncover molecular regulatory networks and investigate the relationship between hub genes and immune cells, immune infiltration analysis was employed. This study provides new insights for exploring the pathophysiological connections and immune mechanisms and excavating the potential biomarkers and therapeutic targets for COVID-19 and pericarditis. The overall flowchart of the study is shown in Figure 1.
2 Materials and methods
2.1 Dataset preparation
By searching the DisGeNET (https://www.disgenet.org/) (18), comparative toxicogenomics database (CTD) (http://ctdbase.org/) (19) and GeneCards (https://www.genecards.org/) (20) databases, we identified genes related to pericarditis and COVID-19. We selected supplementary datasets from Gene Expression Omnibus (GEO) of the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/geo/) (21). GSE164805, platform number GPL26963, the whole genome transcriptome of peripheral blood mononuclear cells was analyzed on five healthy controls and ten COVID-19 patients (22).
2.2 Identification of the common genes of COVID-19 and pericarditis
Based on the scoring standards of the different databases, we collected the top 500 genes from the DisGeNET, CTD, and GeneCards databases when the number was greater than 500. Online GEO analysis tool GEO2R (www.ncbi.nlm.nih.gov/geo/geo2r/) was used to analyze sample data for differential gene expression (21). We utilized GEO2R to identify the differentially expressed genes (DEGs) with a false discovery rate (FDR) < 0.00001 and |log fold-change| > 1 for GSE164805. Subsequently, we integrated these two parts of genes related to COVID-19 and then took the intersection of COVID-19 and pericarditis to obtain common genes using an online Venn tool (http://jvenn.toulouse.inra.fr/app/example.html) (23).
2.3 GO and KEGG pathway enrichment analyses
In order to investigate the probable biological connection between COVID-19 and pericarditis, the clusterProfiler software (version 3.14.3) was used to conduct Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis for the common genes (24). Biological processes (BP), cellular elements (CC), and molecular functions (MF) were all included in the GO analysis. The top 10 GO and top 20 in KEGG items with the lowest p-values were shown as bubble diagrams using an online platform (http://www.bioinformatics.com.cn) for data processing and visualization.
2.4 PPI network analysis and hub genes screening
Identifying unknown protein functional modules from PPI networks is crucial for understanding protein function and interpreting key data in cell biology. PPI network analysis is a promising strategy that can provide a deeper and more comprehensive insight into the relationships between various diseases from the standpoint of protein interactions (25, 26). The online analytical tool STRING (https://string-db.org/) was to study protein interactions, systematically collect and integrate physical regulatory interactions and functional relationships between proteins (27). We constructed a PPI network based on the score greater than 0.4 and analyzed and visualized the results using Cytoscape 3.9.1, which is an open-source project designed to integrate high-throughput data and molecular interaction networks into a unitive framework (28). CytoHubba (https://apps.cytoscape.org/apps/cytohubba), an important Cytoscape plugin for network topology analysis, uses 11 methods for studying key genes from different perspectives.
2.5 Immune infiltration analysis
Immune cells exhibit specific patterns of infiltration and residence. Studying the infiltration status can provide a better understanding of their role and mechanism in disease pathogenesis and can thus be applied to the discovery of new treatment strategies for many diseases (29). The CIBERSORT tool, based on the linear support vector regression, decomposes the expression matrix of subtypes of human immune cells for immune-immersion analysis (30). The proportion of immune cells in GSE164805 was calculated, along with the relevance between immune cells and hub genes, as well as each immune cell.
2.6 Identification of TFs and miRNAs
TFs are proteins to recognize special DNA sequences and are key cellular components forming complex regulatory systems to control gene expression (31). NetworkAnalyst (http://www.networkanalyst.ca) is to conduct complex meta-analyses for gene expression and is suitable for data processing and analysis in the context of PPI networks (32). The construction of the TF–genes was based on the JASPAR database (http://jaspar.genereg.net), which includes TF-binding profiles of multiple species from six taxonomic groups (33). MiRNAs regulate protein expression by binding to TF; research on the interaction network of TF-miRNAs was conducted using the RegNetwork database (http://www.regnetworkweb.org/) (34). In addition to studying the role of target genes and miRNAs with TF, we carried out topology analysis and construction of gene–miRNA networks based on miRTarBase v8.0 (https://miRTarBase.cuhk.edu.cn/) (35).
2.7 Analysis of gene–disease interaction networks
DisGeNET integrates and standardizes disease-related genes and variant data, covering the whole spectrum of human diseases as well as normal and abnormal features (36). Gene–disease network was established was to study diseases related to COVID-19 and pericarditis using the NetworkAnalyst platform.
2.8 Analysis of protein–chemical and protein–drug interaction networks
Constructing protein–chemical and protein–drug networks is conducive to predicting the target information of drugs and chemicals relevant to COVID-19 and pericarditis. In the NetworkAnalyst platform, the corresponding compounds and drugs were identified and obtained using the CTD and DrugBank database.
3 Results
3.1 Identification of common genes of COVID-19 and pericarditis
By searching the DisGeNET, CTD, and GeneCards databases, we identified genes related to COVID-19 and pericarditis. To improve the integration and standardization of the data, we summarized the top 500 genes in each database according to their scoring standards. If the original data were less than 500, all retrieved data were included. Using this rule, we obtained 51, 500, and 500 pericarditis-related genes from the DisGeNET, CTD, and GeneCards databases, respectively. Subsequently, 930 pericarditis-related genes were identified by merging and de-duplicating the results from the three databases.
Using the same method, we obtained 1236 COVID-19-related genes from DisGeNET, CTD, and GeneCards. In addition, we gained 494 COVID-related genes under the settings: FDR < 0.00001 and log fold-change > 1. By merging and deduplicating the results from the three databases and the GEO dataset, we obtained 1711 COVID-related genes. Finally, 313 common genes between COVID-19 and pericarditis were identified through intersections (Figure 2) (Table 1; Supplementary Tables 1–3).
3.2 GO and KEGG pathway enrichment analyses
Based on common genes, GO and KEGG enrichment were carried out to examine their biological roles and signal pathways. Figure 3 displayed the top 10 terms in the BP, CC, and MF categories with the lowest p-values. BP terms mainly involved the regulation of cytokine production and inflammatory response and were associated with the proliferation of immune cells, involving in lymphocytes, leukocytes, and mononuclear cells (Figure 3A). CC terms mainly revealed the external side of the plasma membrane and some lumens, such as secretory granules, cytoplasmic vesicles, and the endoplasmic reticulum (Figure 3B). MF terms mainly demonstrated the activity of signaling receptors and immune receptors which are crucial for the binding of various factors, including proteases, cytokines, and chemokines (Figure 3C). Furthermore, KEGG analysis showed that most pathways were involved in immune- and infection-related diseases, such as COVID-19, influenza, measles, tuberculosis, rheumatoid arthritis, hepatitis and inflammatory bowel disease. Notably, multiple immune-related pathways were also enriched, including cytokine receptor interaction, T helper 17 (Th17) cell differentiation, interleukin 17 (IL-17), tumor necrosis factor (TNF), and Toll-like receptor signaling pathways (Figures 3D, E). All results were visualized using bubble plots, which manifested that common genes might be involved in immune-related functions and pathways, thereby affecting the progression of COVID-19 and pericarditis (Tables 2, 3).
Figure 3 GO and KEGG enrichment analysis of the common genes. (A) Biological processes in bubble chart. (B) Cellular component in bubble chart. (C) Molecular function in bubble chart. (D) The top 20 signaling pathways of KEGG in bubble chart. (E) The top 20 signaling pathways of KEGG in bar graph. The color of the dots reflects the size of the p-values, and the size of the dots reflects the number of annotated genes in the bubble charts. Different colors represent different pathway classifications, and the horizontal axis represents the number of genes per pathway in a bar graph.
3.3 PPI network and hub gene analyses
The common genes of COVID-19 and pericarditis were imported to STRING to create a PPI network, and then were uploaded into Cytoscape for comprehensive analysis to forecast gene interactions and associated pathways. CytoHubba is a plugin that identifies hub nodes and provides 11 analysis algorithms to calculate and sort nodes in the network. We use seven algorithms to calculate the top 20 and then take intersections to screen hub genes: Maximal Clique Centrality (MCC), Maximum Neighborhood Component (MNC), Degree, Closeness, Radiality, Stress and Edge Percolated Component (EPC) (Figure 4A) (Table 4). Interleukin 1 beta (IL-1β), cluster of differentiation 8 antigen (CD8A), interleukin 10 (IL-10), cluster of differentiation 4 (CD4), interleukin 6 (IL-6), Toll-like receptor 4 (TLR4), chemokine ligand 2 (CCL2), and Protein Tyrosine Phosphatase Receptor Type C (PTPRC) were among the top 20 genes identified from the seven algorithm scores. The area under the curve (AUC) results were to assess the specificity and sensitivity of the eight hub genes to COVID-19 using receiver operating characteristic (ROC) curve analysis. The AUC values of the hub genes were greater than 0.75, except for IL-10, indicating that these genes may be potential biomarkers and have a high diagnostic value for disease (Figures 4B–I).
Figure 4 Screening and validation of hub genes. (A) The Venn diagram shows eight overlapping hub genes screened via the six algorithms. (B–I) Results of the ROC curve analysis and AUC values of hub genes in the COVID-19 dataset.
3.4 Immune infiltration analysis
Investigating immune cell infiltration patterns in COVID-19 patients was using the CIBERSORT algorithm. The proportions of 22 immune cells with COVID-19 are shown in Figure 5A. Figure 5B illustrates the distribution of 22 immune cells in COVID-19, and the infiltration of plasma cells, memory resting CD4 T cells, monocytes, M0 macrophages, resting mast cells and neutrophils in tissues from patients with COVID-19 is considerably higher than that in normal tissue (p < 0.05). Compared to the normal group, the proportion of activated CD8 T cells and natural killer (NK) cells in patients with COVID-19 is lower. Moreover, Figure 6 depicts the relationships between the 22 immune cells. Neutrophils, naive CD4 T cells, memory B cells, M0 macrophages, and monocytes all showed negative correlations with CD8 T cells. NK cell activation was negatively correlated with M0 macrophage, monocyte, and dendritic cell activation. The relevance between neutrophils, memory B cells, and M0 macrophages was positive. Naïve CD4 T cells were positively related to gamma delta T cells, memory B cells, and M0 macrophages. Resting memory CD4 T cells and plasma cells showed the positive correlation with resting NK cells.
Figure 5 Immune infiltration analysis. (A) The ratio of 22 immune cells in COVID-19 and control. (B) The proportion and comparison of immune cells in COVID-19 and control.
The relationship between immune cells and hub genes is depicted in Figures 7, 8. CCL2 exhibited a positive correlation with naive B cells and resting mast cells but a negative correlation with eosinophils. CD4 was statistically positively correlated with CD8 T cells, activated NK cells and resting dendritic cells, but negatively connected with naïve CD4 T cells, monocytes and M0 macrophages. CD8A was statistically positively relevant to CD8 T cells, activated NK cells and resting dendritic cells, but negatively connected with gamma delta T cells, naïve CD4 T cells and activated dendritic cells. Activated NK cells and eosinophils statistically linked positively with IL-1β, but resting NK and resting mast cells statistically related negatively with IL-1β. Statistically, there was a positive correlation between IL-6 and M0 macrophages but a negative association with activated NK cells. IL-10 was statistically positively relevant to activated memory CD4 T cells and plasma cells. PTPRC was statistically positively associated with activated NK cells but, negatively correlated to plasma and activated dendritic cells. While CD8 T cells and active NK cells were inversely connected with TLR4, resting memory CD4 T cells, monocytes, resting mast cells, and M0 macrophages were positively related to TLR4.
Figure 7 The correlation between hub genes and immune cells. The correlations of hub genes (CCL2, CD4, CD8A and IL-1β) with 22 immune cells were determined using p < 0.05 as the screening criterion.
Figure 8 The correlation between hub genes and immune cells. The correlations of hub genes (IL-6, IL-10, PTPRC and TLR4) with 22 immune cells were determined using p < 0.05 as the screening criterion.
3.5 Construction of gene regulatory networks
To identify the main variations at the transcriptional level and further study key protein regulatory molecules, we employed a network-based approach to decipher the regulatory TFs and miRNAs. 38 TFs were connected with the eight hub genes; these TFs included MEF2A, POU2F2, CREB1, PPARG, YY1, NR2F1, JUN, FOXC1, NR3C1, and RELA (Figure 9A). Six hub genes corresponded to 35 miRNAs; the miRNAs binding to multiple hub genes were hsa-mir-21-5p, hsa-mir-26b-5p, hsa-mir-24-3p, hsa-mir-335-5p, hsa-mir-1-3p, hsa-mir-146a-5p, hsa-mir-146b-5p, hsa-mir-124-3p, hsa-mir-106a-5p, hsa-mir-155-5p, hsa-mir-98-5p, and hsa-let-7c-5p (Figure 9B).
Figure 9 (A) TF–gene interaction network analysis. Dots represent hub genes; square dots represent transcription factors. Darker colors indicate stronger associations. (B) Gene–miRNA interaction network analysis. Dots represent hub genes; square dots represent miRNAs. Darker colors indicate stronger associations.
3.6 Gene–disease interaction network
The development of technology and solutions for disease treatment begins with studying the links between diseases and genes; the interrelationships between different diseases usually require one or more similar genes (37). Based on DisGeNET, the results showed that the gene–disease network was linked to at least three hub genes. The following diseases had the strongest coordination with the hub genes studied: rheumatoid arthritis, glomerulonephritis, hyperalgesia, inflammation, liver cirrhosis, reperfusion injury, schizophrenia, and major depressive disorder (Figure 10). Notably, these diseases are mostly related to inflammation or immune responses, which have implications for the development of mechanisms and treatment methods for COVID-19 and pericarditis.
Figure 10 Gene–disease interaction network analysis. Dots represent hub genes; square dots represent diseases related to hub genes.
3.7 Protein–chemical and protein–drug Interaction networks
Constructing protein–chemical and protein–drug interaction networks contributes to the exploration of the biological functions of proteins in cells and the research of potential drugs. Only the chemicals linked to at least four hub genes are displayed in Figure 9. The top eight chemicals were methotrexate, antirheumatic agents, nickel, tretinoin, arsenic, benzo(a)pyrene, cadmium, and dexamethasone, demonstrating their tight association with COVID-19 and pericarditis (Figure 11A). Protein–drug network indicates that drugs related to IL-1β and IL-10 may have broader scope for study, with AV411 having potential associations with two genes (Figure 11B).
Figure 11 Protein–chemical and protein–drug interaction network analyses. (A) The interaction between hub genes and potential chemicals. (B) The interaction between hub genes and potential drugs. Dots represent hub genes; square dots represent chemicals or drugs.
4 Discussion
Increasing evidence has linked cardiovascular disease to increased morbidity and mortality from COVID-19, and the burden is evident even among patients who are not hospitalized (38, 39). COVID-19 is a real-time global pandemic, and this virus infection is also a pathogenic factor of pericarditis. COVID-19 has significantly increased the risk of pericarditis, and because multiple monitoring reports suggest that the COVID-19 vaccine may also increase the likelihood of people suffering from pericarditis, many studies on vaccines and pericarditis have been conducted in different countries and regions (40, 41). Currently, most reports on COVID-19 and pericarditis focus on the epidemiology and vaccines; however, research on the potential mechanisms of comorbidity remains lacking. Therefore, our research aimed to reveal the etiology and mechanism of COVID-19 and pericarditis from the perspective of molecular regulation, based on network data mining and bioinformatics analysis.
GO is a bioinformatics resource that provides gene product functions and uses ontology to represent biological knowledge and it can identify the biological process of shared genes in this study (42). The enrichment results in the BP and CC were associated with immune cell proliferation, including that of white blood cells, lymphocytes, and monocytes, suggesting their involvement in regulating cytokine production and inflammatory responses. Lymphocyte count and cytokine levels are closely correlated with disease severity, which is of great significance for the early diagnosis, treatment, and prognosis of COVID-19 (43). Gerd et al. examined the immune cell spectrum in the cerebrospinal fluid with COVID-19 and the results indicated an expansion of dedifferentiated monocytes and interferon signature of leukocytes (44). Anti-inflammatory treatment can reduce the thickness of pericardial late gadolinium enhancement, alleviate pericardial and systemic inflammation, and improve the physiological status and symptoms of pericarditis (45). Studies have focused on the prognostic value of inflammatory markers in active pericarditis, suggesting that high-sensitivity C-reactive protein may affect the intensity and duration of pericarditis (46). MF results mainly involve the binding of proteases and chemokines as well as the activity of signaling and immune receptors. Assessing the chemokine status after SARS-CoV-2 infection and detecting the “immune signature” is crucial for individual risk stratification (47). Studying the extensive cytokine releasing syndrome in COVID-19 may be helpful for targeting chemokines and growth factors as therapeutic drugs; furthermore, autoantibodies targeting chemokines may inhibit the potentially harmful immune response observed in patients with COVID-19 (48, 49). KEGG connects genomic and higher-order functional information from the perspective of genes and molecular networks and annotates up-to-date gene catalogs and functions (50). The top 20 pathways in KEGG involved immune-related pathways such as COVID-19, influenza, hepatitis, rheumatoid arthritis, and inflammatory bowel disease, and involved in immune pathways such as Th17 cell differentiation, IL-17, TNF, and Toll-like receptors. TNF-α blockers have made important progress in the treatment of idiopathic recurrent pericarditis (51). Christian et al. found that after virus clearance, clonally expanded Th17 cells remained in the lungs, which was associated with the potentially pathogenic cytokine expression profile of IL-17, and interacted with cytotoxic CD8 T cells and macrophages (52). As key regulatory factors of the innate immune system, Toll-like receptors recognize viral particles and induce the secretion of pro-inflammatory cytokines, which may also be potential targets for vaccine production (53). In this study, we constructed the gene network to obtain the common differential genes of COVID-19 and pericarditis, so as to obtain the hub gene and explore the potential mechanism correlation between them. GO term functional enrichment and KEGG pathway enrichment were applied to identify functional changes caused by differential genes and their effects in the pathway. However, due to the need for a clear threshold for GO/KEGG enrichment of the common genes, it is possible to miss genes with significant biological significance. GSEA, which may not require the clear threshold and is based on overall trend analysis, can be implemented in the future. Meanwhile, Weighted correlation network analysis (WGCNA) can serve as a supplement to gene network studies and can analyze gene modules that coordinate expression.
Immune cell infiltration showed that patients with COVID-19 had significantly higher levels of plasma cells, resting memory CD4 rest cells, monocytes, M0 macrophages, resting mast cells, and neutrophils than the healthy population. Persistent antibody protection, produced by memory B cells and long-lived plasma cells, is the pillar of the “arms race” between vaccines immunity and the constantly mutating SARS-CoV-2 virus infection, and long-lived plasma cells in bone marrow tissue are the source of these persistent “memory” antibodies during acute infection (54). Some studies have evaluated whether pre-existing cross-reactive memory T cells affect vaccine immunity, and the results displayed that subjects with memory CD4 T cells have stronger antibody responses to vaccines (55). Notably, long-lasting memory T cells responded to SARS-CoV-2 and exhibited substantial cross-reactivity with the N protein of SARS-CoV-2 in patients recovering from SARS-CoV-2 in 2003 (56). Excessive infiltration of macrophages and monocytes into organs is a critical driver of severe COVID-19, and the activation of pulmonary macrophages from infiltrating monocytes results in the recruitment of cytotoxic effector cells and the release of pro-inflammatory cytokines (57). Macrophages activate inflammasomes, which oppose host infection and promote tissue repair by releasing interleukin and inducing pyroptosis; however, it should be noted that macrophage activation syndrome induced by macrophage dysfunction may cause damage to the host (58, 59). The high density of mast cells is related to the activation and release of proteases, which are affected by soluble factors released by T cells with the help of stem cell factors. The proliferation and activation of mast cells are manifestations of inflammatory cell changes in severe and lethal SARS-CoV-2 infection; thus, regulating mast cells and their pro-inflammatory mediators may be a potentially effective treatment for COVID-19 (60, 61). The characteristics of neutrophils in severe COVID-19 include the formation and degradation of neutrophil extracellular traps, expansion and infiltration of neutrophils into the lungs, and activation and immune suppression of neutrophil subsets in the circulatory system (62). The abnormal response of neutrophils after infection with SARS-CoV-2 may be related to uncontrolled viral replication and exacerbated inflammation. Assessing the number, function, and status of neutrophils are crucial for distinguishing the disease severity and identifying the clinical deterioration risk (63). Many studies have been conducted on the progression of immune cells in pericarditis. Neutrophils and macrophages produce a large number of cytokines through the activation of inflammasomes, which contribute to the immune pathogenesis of recurrent pericarditis (64). Interference with neutrophil chemotaxis and adhesion, reduction of recruitment to damaged tissues, and superoxide production are effective strategies and mechanisms for the anti-inflammatory treatment of pericarditis (65). The pericardial interstitial cells of patients with pericarditis exhibit senescent features that induce structural remodeling of the pericardium, such as increased collagen matrix secretion and calcium deposition, promotion of chemotaxis of monocytes/lymphocytes, and recruitment of inflammatory factors (66). Moreover, many clinical samples have been collected from patients with pericarditis. Pleural biopsy revealed pleural pericarditis accompanied by lymphoplasmacytic inflammation, such as IgG4-positive plasma cells, and detection of pericardial fluid indicated that the histamine receptor depended on mast cells to infiltrate the pericardial tissue and was involved in the inflammatory reaction (67, 68). Notably, this study analyzed the difference of immune infiltration in COVID-19 and identified the types of immune cells associated with hub gene. With the accumulation of research data related to pericarditis, immune infiltration analysis of pericarditis can be carried out in this study.
Based on the PPI network and topological analysis, IL-1β, CD8A, IL-10, CD4, IL-6, TLR4, CCL2, and PTPRC were identified as hub genes. Since the COVID-19 pandemic, many trials have found that the COVID-19 group had higher levels of IL-1, IL-6, and IL-10 than the control group. Logistic regression and ROC analyses have revealed that these cytokines have a predictive effect on disease severity (69, 70). A genome-wide association study showed that patients with critical COVID-19 had significantly greater blood IL-6 expression levels than patients without symptoms and that an allele change at the rs2069837 site can reduce IL-6 levels to prevent critical conditions (71). Anti-interleukin (IL)-1 drugs have been developed and used to treat autoimmune and rheumatic immune diseases. As one of the family members with the strongest pro-inflammatory effects, IL-1β is considered the therapeutic target for recurrent idiopathic pericarditis (72). Research has been conducted on the etiology, immune mechanisms, and treatment of tuberculous pericarditis, and showed that IL-10 levels are elevated in the pericardium and blood (73). Myocardial fibrosis is considered a non-negligible feature of constrictive pericarditis. IL-6 mediates abnormal Ca2+ handling and induces atrial fibrosis in sterile pericarditis rats (74). Patients with severe COVID-19 exhibit less pronounced increases in TLR4 expression on CD14 monocytes than those with mild COVID-19, which is related to activation of TLR4/NF-κB axis after lipopolysaccharide stimulation (75). TLR4 signaling pathway also regulates myocardial fibrosis by inhibiting its target genes (76). The detection of specific chemokines in the plasma at the mRNA and protein levels suggests that higher concentrations of CCL2 are associated with the severity of COVID-19, which has potential as a prognostic factor (77). The immunopathological changes in the spleen of patients with COVID-19 are also worthy of attention as they involve the functions of plasma cells and monocytes/macrophages and a decrease in CD8A abundance (78). Through the analysis of transcriptome data, PTPRC was shown to be an important inflammatory and immunomodulatory signature in COVID-19, and that it has high binding efficiency with related drugs in clinical transformation research (79). There have been some studies on CCL2, CD8A, and PTPRC in cardiovascular diseases; however, their roles in pericarditis require further research. Overall, these hub genes may be potential immune regulatory pivots in COVID-19 and pericarditis. In addition, there is a very interesting issue worth discussing. The expression of pro-inflammatory cytokines such as IL-1β and IL-6 may change over time, and their dynamics may be potential predictors of disease (80). Currently, database-based research considers gene expression at different time points as a whole, and more in-depth studies in the future will focus on dynamic changes in gene expression to achieve precise intervention for diseases.
We constructed TF-gene and gene-miRNA interaction networks to better understand the molecular regulation between COVID-19 and pericarditis. CREB1, YY1, FOXC1, and NR3C1 were the TFs having the strongest correlation to the hub genes. We analyzed the transcriptome RNA-seq data related to COVID-19 and used bioinformatics to decode the molecular tags and pathways of the host cell response to SARS CoV-2. The genes YY1 and CREB1 may co-regulate autophagy to affect severe conditions. FOXC1 and YY1 may have good binding affinities to candidate drugs. The NR3C1-CXCL8-neutrophil axis may determine the severity of COVID-19 (81–83). The miRNAs that strongly interacted with the hub genes were hsa-mir-335-5p, hsa-mir-1-3p, hsa-mir-106a-5p, and hsa-mir-98-5p. In different studies, peripheral blood mononuclear cells, serum samples, and bronchial aspirates from patents with COVID-19 and healthy individuals were collected. Sequencing analysis has shown that miR-1-3p is involved in the regulation of autophagy and has high specificity and sensitivity for predicting mortality (84, 85). MiR-335-5p is regulated by angiotensin-converting enzyme and histone deacetylase and is involved in drug development to interfere with host-virus interactions (86). TMPRSS2 is a potential therapeutic target for COVID-19, and miR-98-5p is a regulatory factor of TMPRSS2 that originates from two types of endothelial cells in the lungs and umbilical vein (87). Due to the lack of transcriptome sequencing and network information analysis of pericarditis, the regulatory roles of these TFs and miRNAs in pericarditis need to be further elucidated.
According to the gene-disease interaction network, COVID-19 combined with pericarditis can damage the heart, liver, kidney, and other organs; trigger inflammation and rheumatism; and cause neurological and psychiatric diseases. Several cases have reported that simultaneous onset of glomerulonephritis and pericarditis in patients with rheumatic immunity and viral infection (88, 89), and many studies have focused on COVID-19 vaccine that may increase the risk of glomerulonephritis similar to pericarditis (90, 91). The humoral immune response of patients with liver cirrhosis after COVID-19 vaccination is being explored, and case suggest that constrictive pericarditis may appear as a comorbidity in patients with liver cirrhosis (92, 93). The increased incidence of rheumatoid arthritis in patients with COVID-19 may be due to the impaired function of the autoimmune system and the iatrogenic effect of immunosuppressants, and genome-wide cross-trait analysis shows that higher genetic susceptibility to rheumatoid arthritis also increases the risk of COVID-19 (94, 95). A similar mechanism reveals the possibility of antirheumatic drugs as a potential treatment for COVID-19. The effects of drugs on pericardial contractions secondary to rheumatoid arthritis have also been studied in patients with pericarditis (96, 97). Globally, major depressive disorder and anxiety disorder cases have increased by 27.6% and 25.6%, respectively, as a result of the COVID-19 pandemic, resulting in 49.4 million and 44.5 million DALYs (disability-adjusted life years) (98). Immune dysfunction caused by infection can aggravate mental sequelae, and studies on the influence of COVID-19 on mental health have found that the levels of inflammatory markers are directly proportional to depression severity of depression (99). Tryptophan metabolism may be correlated with the potential susceptibility to depression, and tryptophan supplementation may improve depressive symptoms in patients with COVID-19 treated with drugs that can affect tryptophan metabolism (100). Numerous investigations have shown that COVID-19 increases the risk of schizophrenia, and willingness to be vaccinated is related to the severity of psychiatric symptoms (101, 102). A case report of clozapine-related pericarditis in a patient with refractory schizophrenia during the drug titration phase suggested that great attention should be paid to the side effects of antipsychotics and antidepressants in patients with COVID-19 and pericarditis (103).
Protein–chemical interaction networks indicated that methotrexate, antirheumatic agents, nickel, tretinoin, arsenic, benzo(a)pyrene, cadmium, and dexamethasone have a high correlation with hub genes. As a specific immunosuppressive drug, methotrexate impairs immunogenicity and raises the risk of infection and poor prognosis (104). Interruption of methotrexate for two weeks enhances antibody responses in patients with immune-mediated inflammatory diseases after vaccination (105). IL-6 and IL-1β are pivotal targets of antirheumatic agents, and there is evidence that blocking the IL-6 receptor can reduce lung involvement and acute cardiovascular complications in patients with COVID-19 by inhibiting the systemic inflammatory response (106). Several clinical trials have been conducted to evaluate the long-term prognosis of COVID-19 with different doses of dexamethasone and whether the clinical benefits are related to different respiratory support modes (107, 108). Methotrexate has previously been used for the treatment of purulent pericarditis in rheumatoid arthritis; however, methotrexate-induced pericarditis and pericardial effusion should be considered (109, 110). In addition, clinical guidelines indicate that non-steroidal anti-inflammatory drugs (NSAIDs) such as aspirin are recommended as effective drugs for the first-line treatment of pericarditis (111). Notably, when the screening scope of the interaction network is expanded, NSAIDs such as aspirin and ibuprofen can also be searched. Protein–drug interaction networks revealed that some drugs, such as AV411, minocycline, rilonacept, canakinumab, XOMA 052, and VX-765, exert therapeutic effects by targeting hub genes. AV411 reduces opioid withdrawal by inhibiting glial pro-inflammatory responses, whereas minocycline prevents potentially fatal arrhythmias by inhibiting pro-inflammatory cytokines and poly (ADP-ribose) polymerase-1 associated with SARS-CoV-2 replication (112, 113). Rilonacept is a trap for IL-1β and has been shown in clinical trials to inhibit recurrent pericarditis episodes and prevent the recurrence of pericarditis (114). Canakinumab, a human monoclonal antibody targeting IL-1β, is associated with the reduction of serum C-reactive protein level and the improvement of overall mortality in COVID-19; case reports showed that canakinumab can reduce the risk of recurrence of systemic disease-related pericarditis (115, 116). As the neutralizing antibody to IL-1β, XOMA 052 has a rapid onset and sustained control of intraocular inflammation, and VX-765 ameliorates myocardial reperfusion injury by inhibiting caspase-1 activity and reducing lactate dehydrogenase release (117, 118). In summary, some chemicals and drugs predicted based on hub genes have been proven in clinical trials and experimental studies of COVID-19 and pericarditis, whereas others deserve further exploration.
This study had several limitations. There are currently no suitable microarray or RNA sequencing data for pericarditis, resulting in a lack of available datasets that may prevent the acquisition and identification of sufficient DEGs. In addition, our study was purely based on bioinformatics analysis and requires subsequent in vivo and in vitro to confirm the validity of the results, as well as to fully evaluate the biological function of the hub gene and the clinical value of the drug.
5 Conclusion
In recent years, the topic of increased risk of pericarditis caused by COVID-19 has triggered a large number of studies and heated discussions, but there is still a lack of exploration and research on the mechanism of COVID-19 and pericarditis. The immunological mechanisms and common genes linked to COVID-19 and pericarditis were identified in this investigation. The eight hub genes (IL-1β, CD8A, IL-10, CD4, IL-6, TLR4, CCL2, and PTPRC) are relatively mature and have been extensively studied in immune regulation, and some also have the potential to affect immune functions. Thus, COVID-19 and pericarditis exhibit complex interactions. The enrichment analysis and various interaction networks constructed and analyzed in this study revealed the molecular mechanisms of COVID-19 and pericarditis from multiple perspectives. Based on the analysis, some potential compounds and drugs were predicted. However, further research on their functions and mechanisms is required to provide new ideas for identify potential biomarkers and explore appropriate treatment methods.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.
Author contributions
DL: Formal Analysis, Investigation, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing. RC: Visualization, Writing – review & editing. CH: Writing – review & editing. GZ: Writing – review & editing. ZL: Writing – review & editing. XX: Writing – review & editing. BW: Writing – review & editing. BL: Writing – review & editing, Writing – original draft. XC: Writing – review & editing, Writing – original draft.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by The National Natural Science Foundation of China (no. 82172574, no. 81871231), The Natural Science Foundation of Shandong Province (no. ZR2020MH016, no. ZR202209280042), Research Planning Project of Shandong Higher Medical Education Research Center (no. YJKT202171), Shandong Taishan Scholars Young Experts Program (No. tsqn202103056), Project of China Association of Chinese Medicine (no. 2021HH-006)
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2024.1264856/full#supplementary-material
Supplementary Figure 1 | Schematic overview for gene set enrichment analysis. GSEA is a computational method that determines whether an a priori defined set of genes shows statistically significant, concordant differences between two biological states. An expression data set sorted by phenotype, heat map and the gene tags, and plot of the running sum for dataset, including the location of the maximum enrichment score (ES) and the leading-edge subset (https://www.gsea-msigdb.org/gsea/) (16).
Abbreviations
COVID-19, coronavirus disease 2019; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2; WHO, World Health Organization; ARDS, acute respiratory distress syndrome; miRNA, micro-RNA; TFs, transcription factors; PPI, protein–protein interaction; ORA, over-representative analysis; FCS, function set scoring; PT, pathway topology; NT, network topology; WKS, Weighted Kolmogorov Smirnov; CTD, comparative toxicogenomics database; GEO, Gene Expression Omnibus; FDR, false discovery rate; CC, cellular components; BP, biological processes; MF, molecular functions; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes; ROC, receiver operating characteristic; AUC, area under the curve; Th17, T helper 17; IL-17, interleukin 17; TNF, tumor necrosis factor; MCC, Maximal Clique Centrality; MNC, Maximum Neighborhood Component; IL-1β, interleukin 1 beta; CD8A, cluster of differentiation 8 antigen; IL-10, interleukin 10; CD4, cluster of differentiation 4; IL-6, interleukin 6; TLR4, Toll-like receptor 4; CCL2, chemokine ligand 2; PTPRC, Protein Tyrosine Phosphatase Receptor Type C; NK, natural killer; WGCNA, Weighted correlation network analysis; NSAIDs, non-steroidal anti-inflammatory drugs.
References
1. Wiersinga WJ, Rhodes A, Cheng AC, Peacock SJ, Prescott HC. Pathophysiology, transmission, diagnosis, and treatment of coronavirus disease 2019 (COVID-19): A review. Jama. (2020) 324:782–93. doi: 10.1001/jama.2020.12839
2. Guan WJ, Zhong NS. Clinical characteristics of COVID-19 in China. Reply. N Engl J Med. (2020) 382:1861–2. doi: 10.1056/NEJMc2005203
3. Xu Z, Shi L, Wang Y, Zhang J, Huang L, Zhang C, et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir Med. (2020) 8:420–2. doi: 10.1016/s2213-2600(20)30076-x
4. Troughton RW, Asher CR, Klein AL. Pericarditis. Lancet. (2004) 363:717–27. doi: 10.1016/s0140-6736(04)15648-1
5. Imazio M, Gaita F, LeWinter M. Evaluation and treatment of pericarditis: A systematic review. Jama. (2015) 314:1498–506. doi: 10.1001/jama.2015.12763
6. Fairweather D, Beetler DJ, Di Florio DN, Musigk N, Heidecker B, Cooper LT Jr. COVID-19, myocarditis and pericarditis. Circ Res. (2023) 132:1302–19. doi: 10.1161/circresaha.123.321878
7. Le Vu S, Bertrand M, Jabagi MJ, Botton J, Drouin J, Baricault B, et al. Age and sex-specific risks of myocarditis and pericarditis following COVID-19 messenger RNA vaccines. Nat Commun. (2022) 13:3633. doi: 10.1038/s41467-022-31401-5
8. Wong HL, Hu M, Zhou CK, Lloyd PC, Amend KL, Beachler DC, et al. Risk of myocarditis and pericarditis after the COVID-19 mRNA vaccination in the USA: a cohort study in claims databases. Lancet. (2022) 399:2191–9. doi: 10.1016/s0140-6736(22)00791-7
9. Schultze JL, Aschenbrenner AC. COVID-19 and the human innate immune system. Cell. (2021) 184:1671–92. doi: 10.1016/j.cell.2021.02.029
10. Brodin P. Immune determinants of COVID-19 disease presentation and severity. Nat Med. (2021) 27:28–33. doi: 10.1038/s41591-020-01202-8
11. Li X, Geng M, Peng Y, Meng L, Lu S. Molecular immune pathogenesis and diagnosis of COVID-19. J Pharm Anal. (2020) 10:102–8. doi: 10.1016/j.jpha.2020.03.001
12. Altan M, Toki MI, Gettinger SN, Carvajal-Hausdorf DE, Zugazagoitia J, Sinard JH, et al. Immune checkpoint inhibitor-associated pericarditis. J Thorac Oncol. (2019) 14:1102–8. doi: 10.1016/j.jtho.2019.02.026
13. Cantarini L, Imazio M, Brucato A, Lucherini OM, Galeazzi M. Innate versus acquired immune response in the pathogenesis of recurrent idiopathic pericarditis. Autoimmun Rev. (2010) 9:436–40. doi: 10.1016/j.autrev.2009.11.021
14. Davies AJ. Immunological tolerance and the autoimmune response. Autoimmun Rev. (2008) 7:538–43. doi: 10.1016/j.autrev.2008.04.007
15. Bu D, Luo H, Huo P, Wang Z, Zhang S, He Z, et al. KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. (2021) 49:W317–w25. doi: 10.1093/nar/gkab447
16. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
17. Hung JH, Yang TH, Hu Z, Weng Z, DeLisi C. Gene set enrichment analysis: performance evaluation and usage guidelines. Brief Bioinform. (2012) 13:281–91. doi: 10.1093/bib/bbr049
18. Piñero J, Saüch J, Sanz F, Furlong LI. The DisGeNET cytoscape app: Exploring and visualizing disease genomics data. Comput Struct Biotechnol J. (2021) 19:2960–7. doi: 10.1016/j.csbj.2021.05.015
19. Davis AP, Wiegers TC, Johnson RJ, Sciaky D, Wiegers J, Mattingly CJ. Comparative toxicogenomics database (CTD): update 2023. Nucleic Acids Res. (2023) 51:D1257–d62. doi: 10.1093/nar/gkac833
20. Safran M, Rosen N, Twik M, BarShir R, Stein TI, Dahary D, et al. The geneCards suite. In: Abugessaisa I, Kasukawa T, editors. Practical guide to life science databases. Springer Nature Singapore, Singapore (2021). p. 27–56.
21. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. (2013) 41:D991–5. doi: 10.1093/nar/gks1193
22. Zhang Q, Meng Y, Wang K, Zhang X, Chen W, Sheng J, et al. Inflammation and antiviral immune response associated with severe progression of COVID-19. Front Immunol. (2021) 12:631226. doi: 10.3389/fimmu.2021.631226
23. Bardou P, Mariette J, Escudié F, Djemiel C, Klopp C. jvenn: an interactive Venn diagram viewer. BMC Bioinf. (2014) 15:293. doi: 10.1186/1471-2105-15-293
24. 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
25. Li T, Wernersson R, Hansen RB, Horn H, Mercer J, Slodkowicz G, et al. A scored human protein-protein interaction network to catalyze genomic interpretation. Nat Methods. (2017) 14:61–4. doi: 10.1038/nmeth.4083
26. Chen H, Cai Y, Ji C, Selvaraj G, Wei D, Wu H. AdaPPI: identification of novel protein functional modules via adaptive graph convolution networks in a protein-protein interaction network. Brief Bioinform. (2023) 24(1):bbac523. doi: 10.1093/bib/bbac523
27. Szklarczyk D, Kirsch R, Koutrouli M, Nastou K, Mehryary F, Hachilif R, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. (2023) 51:D638–d46. doi: 10.1093/nar/gkac1000
28. 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
29. Makowski L, Chaib M, Rathmell JC. Immunometabolism: From basic mechanisms to translation. Immunol Rev. (2020) 295:5–14. doi: 10.1111/imr.12858
30. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
31. Lambert SA, Jolma A, Campitelli LF, Das PK, Yin Y, Albu M, et al. The human transcription factors. Cell. (2018) 172:650–65. doi: 10.1016/j.cell.2018.01.029
32. Xia J, Gill EE, Hancock RE. NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data. Nat Protoc. (2015) 10:823–44. doi: 10.1038/nprot.2015.052
33. Fornes O, Castro-Mondragon JA, Khan A, van der Lee R, Zhang X, Richmond PA, et al. JASPAR 2020: update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. (2020) 48:D87–d92. doi: 10.1093/nar/gkz1001
34. Liu ZP, Wu C, Miao H, Wu H. RegNetwork: an integrated database of transcriptional and post-transcriptional regulatory networks in human and mouse. Database (Oxford). (2015) 2015:bav095. doi: 10.1093/database/bav095
35. Huang HY, Lin YC, Cui S, Huang Y, Tang Y, Xu J, et al. miRTarBase update 2022: an informative resource for experimentally validated miRNA-target interactions. Nucleic Acids Res. (2022) 50:D222–d30. doi: 10.1093/nar/gkab1079
36. Piñero J, Ramírez-Anguita JM, Saüch-Pitarch J, Ronzano F, Centeno E, Sanz F, et al. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. (2020) 48:D845–d55. doi: 10.1093/nar/gkz1021
37. Mahmud SMH, Al-Mustanjid M, Akter F, Rahman MS, Ahmed K, Rahman MH, et al. Bioinformatics and system biology approach to identify the influences of SARS-CoV-2 infections to idiopathic pulmonary fibrosis and chronic obstructive pulmonary disease patients. Brief Bioinform. (2021) 22(5):bbab115. doi: 10.1093/bib/bbab115
38. Clerkin KJ, Fried JA, Raikhelkar J, Sayer G, Griffin JM, Masoumi A, et al. COVID-19 and cardiovascular disease. Circulation. (2020) 141:1648–55. doi: 10.1161/circulationaha.120.046941
39. 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
40. Pillay J, Gaudet L, Wingert A, Bialy L, Mackie AS, Paterson DI, et al. Incidence, risk factors, natural history, and hypothesised mechanisms of myocarditis and pericarditis following covid-19 vaccination: living evidence syntheses and review. Bmj. (2022) 378:e069445. doi: 10.1136/bmj-2021-069445
41. Diaz GA, Parsons GT, Gering SK, Meier AR, Hutchinson IV, Robicsek A. Myocarditis and pericarditis after vaccination for COVID-19. Jama. (2021) 326:1210–2. doi: 10.1001/jama.2021.13443
42. Drabkin H, Hill D, Carbon S, Dietze H, Lomax J, Osumi-Sutherland D, et al. Gene Ontology Consortium: going forward. Nucleic Acids Res. (2015) 43:D1049–56. doi: 10.1093/nar/gku1179
43. Idiz UO, Yurttas TT, Degirmencioglu S, Orhan B, Erdogan E, Sevik H, et al. Immunophenotyping of lymphocytes and monocytes and the status of cytokines in the clinical course of Covid-19 patients. J Med Virol. (2022) 94:4744–53. doi: 10.1002/jmv.27917
44. Heming M, Li X, Räuber S, Mausberg AK, Börsch AL, Hartlehnert M, et al. Neurological manifestations of COVID-19 feature T cell exhaustion and dedifferentiated monocytes in cerebrospinal fluid. Immunity. (2021) 54:164–75.e6. doi: 10.1016/j.immuni.2020.12.011
45. Feng D, Glockner J, Kim K, Martinez M, Syed IS, Araoz P, et al. Cardiac magnetic resonance imaging pericardial late gadolinium enhancement and elevated inflammatory markers can predict the reversibility of constrictive pericarditis after antiinflammatory medical therapy: a pilot study. Circulation. (2011) 124:1830–7. doi: 10.1161/circulationaha.111.026070
46. Kumar AK, Yesilyaprak A, Furqan MM, Jain V, Montane B, Imazio M, et al. Prognostic value of inflammatory markers in idiopathic recurrent pericarditis. J Am Coll Cardiol. (2022) 79:1644–5. doi: 10.1016/j.jacc.2022.02.016
47. Coperchini F, Chiovato L, Ricci G, Croce L, Magri F, Rotondi M. The cytokine storm in COVID-19: Further advances in our understanding the role of specific chemokines involved. Cytokine Growth Factor Rev. (2021) 58:82–91. doi: 10.1016/j.cytogfr.2020.12.005
48. Xu ZS, Shu T, Kang L, Wu D, Zhou X, Liao BW, et al. Temporal profiling of plasma cytokines, chemokines and growth factors from mild, severe and fatal COVID-19 patients. Signal Transduct Target Ther. (2020) 5:100. doi: 10.1038/s41392-020-0211-1
49. Bordon Y. Autoantibodies against chemokines linked with better disease outcomes in COVID-19. Nat Rev Immunol. (2023) 23:203. doi: 10.1038/s41577-023-00860-2
50. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. doi: 10.1093/nar/28.1.27
51. Nieto González JC, Monteagudo Saez I, López-Longo FJ, Serrano B, Mata Martínez C, Carreño Pérez L. Idiopathic recurrent pericarditis treated successfully with tumour necrosis factor alpha blocking agents (anti-TNF-α). Clin Exp Rheumatol. (2013) 31:776–8.
52. Zhao Y, Kilian C, Turner JE, Bosurgi L, Roedl K, Bartsch P, et al. Clonal expansion and activation of tissue-resident memory-like Th17 cells expressing GM-CSF in the lungs of severe COVID-19 patients. Sci Immunol. (2021) 6:eabf6692. doi: 10.1126/sciimmunol.abf6692
53. Khanmohammadi S, Rezaei N. Role of Toll-like receptors in the pathogenesis of COVID-19. J Med Virol. (2021) 93:2735–9. doi: 10.1002/jmv.26826
54. Nguyen DC, Lamothe PA, Woodruff MC, Saini AS, Faliti CE, Sanz I, et al. COVID-19 and plasma cells: Is there long-lived protection? Immunol Rev. (2022) 309:40–63. doi: 10.1111/imr.13115
55. Mateus J, Dan JM, Zhang Z, Rydyznski Moderbacher C, Lammers M, Goodwin B, et al. Low-dose mRNA-1273 COVID-19 vaccine generates durable memory enhanced by cross-reactive T cells. Science. (2021) 374:eabj9853. doi: 10.1126/science.abj9853
56. Le Bert N, Tan AT, Kunasegaran K, Tham CYL, Hafezi M, Chia A, et al. SARS-CoV-2-specific T cell immunity in cases of COVID-19 and SARS, and uninfected controls. Nature. (2020) 584:457–62. doi: 10.1038/s41586-020-2550-z
57. Knoll R, Schultze JL, Schulte-Schrepping J. Monocytes and macrophages in COVID-19. Front Immunol. (2021) 12:720109. doi: 10.3389/fimmu.2021.720109
58. Sefik E, Qu R, Junqueira C, Kaffe E, Mirza H, Zhao J, et al. Inflammasome activation in infected macrophages drives COVID-19 pathology. Nature. (2022) 606:585–93. doi: 10.1038/s41586-022-04802-1
59. Merad M, Martin JC. Pathological inflammation in patients with COVID-19: a key role for monocytes and macrophages. Nat Rev Immunol. (2020) 20:355–62. doi: 10.1038/s41577-020-0331-4
60. Schaller T, Märkl B, Claus R, Sholl L, Hornick JL, Giannetti MP, et al. Mast cells in lung damage of COVID-19 autopsies: A descriptive study. Allergy. (2022) 77:2237–9. doi: 10.1111/all.15293
61. Theoharides TC, Guerra L, Patel K. Successful treatment of a patient with severe COVID-19 using an integrated approach addressing mast cells and their mediators. Int J Infect Dis. (2022) 118:164–6. doi: 10.1016/j.ijid.2022.02.049
62. Hazeldine J, Lord JM. Neutrophils and COVID-19: active participants and rational therapeutic targets. Front Immunol. (2021) 12:680134. doi: 10.3389/fimmu.2021.680134
63. Dennison D, Al Khabori M, Al Mamari S, Aurelio A, Al Hinai H, Al Maamari K, et al. Circulating activated neutrophils in COVID-19: An independent predictor for mechanical ventilation and death. Int J Infect Dis. (2021) 106:155–9. doi: 10.1016/j.ijid.2021.03.066
64. Lopalco G, Rigante D, Cantarini L, Imazio M, Lopalco A, Emmi G, et al. The autoinflammatory side of recurrent pericarditis: Enlightening the pathogenesis for a more rational treatment. Trends Cardiovasc Med. (2021) 31:265–74. doi: 10.1016/j.tcm.2020.04.006
65. Bonaventura A, Montecucco F. Inflammation and pericarditis: Are neutrophils actors behind the scenes? J Cell Physiol. (2019) 234:5390–8. doi: 10.1002/jcp.27436
66. Han L, Li X, Zhang G, Xu Z, Gong D, Lu F, et al. Pericardial interstitial cell senescence responsible for pericardial structural remodeling in idiopathic and postsurgical constrictive pericarditis. J Thorac Cardiovasc Surg. (2017) 154:966–75.e4. doi: 10.1016/j.jtcvs.2017.03.115
67. Sekiguchi H, Horie R, Utz JP, Ryu JH. IgG4-related systemic disease presenting with lung entrapment and constrictive pericarditis. Chest. (2012) 142:781–3. doi: 10.1378/chest.11-2608
68. Maslińska D, Laure-Kamionowska M, Rzodkiewicz P, Lukasik-Wojtecka E, Szukiewicz D, Karolczak M, et al. Histamine in pericarditis of children with congenital heart malformations. Inflammation Res. (2010) 59 Suppl 2:S259–61. doi: 10.1007/s00011-009-0121-8
69. Schultheiß C, Willscher E, Paschold L, Gottschick C, Klee B, Henkes SS, et al. The IL-1β, IL-6, and TNF cytokine triad is associated with post-acute sequelae of COVID-19. Cell Rep Med. (2022) 3:100663. doi: 10.1016/j.xcrm.2022.100663
70. Han H, Ma Q, Li C, Liu R, Zhao L, Wang W, et al. Profiling serum cytokines in COVID-19 patients reveals IL-6 and IL-10 are disease severity predictors. Emerg Microbes Infect. (2020) 9:1123–30. doi: 10.1080/22221751.2020.1770129
71. Gong B, Huang L, He Y, Xie W, Yin Y, Shi Y, et al. A genetic variant in IL-6 lowering its expression is protective for critical patients with COVID-19. Signal Transduct Target Ther. (2022) 7:112. doi: 10.1038/s41392-022-00923-1
72. Imazio M, Lazaros G, Gattorno M, LeWinter M, Abbate A, Brucato A, et al. Anti-interleukin-1 agents for pericarditis: a primer for cardiologists. Eur Heart J. (2022) 43:2946–57. doi: 10.1093/eurheartj/ehab452
73. Isiguzo G, Du Bruyn E, Howlett P, Ntsekhe M. Diagnosis and management of tuberculous pericarditis: what is new? Curr Cardiol Rep. (2020) 22:2. doi: 10.1007/s11886-020-1254-1
74. Liao J, Zhang S, Yang S, Lu Y, Lu K, Wu Y, et al. Interleukin-6-mediated-ca(2+) handling abnormalities contributes to atrial fibrillation in sterile pericarditis rats. Front Immunol. (2021) 12:758157. doi: 10.3389/fimmu.2021.758157
75. Dorneles GP, Teixeira PC, Peres A, Rodrigues Júnior LC, da Fonseca SG, Monteiro MC, et al. Endotoxin tolerance and low activation of TLR-4/NF-κB axis in monocytes of COVID-19 patients. J Mol Med (Berl). (2023) 101:183–95. doi: 10.1007/s00109-023-02283-x
76. Xiao Y, Qiao W, Wang X, Sun L, Ren W. MiR-146a mediates TLR-4 signaling pathway to affect myocardial fibrosis in rat constrictive pericarditis model. J Thorac Dis. (2021) 13:935–45. doi: 10.21037/jtd-20-2716
77. Pius-Sadowska E, Niedźwiedź A, Kulig P, Baumert B, Sobuś A, Rogińska D, et al. CXCL8, CCL2, and CMV seropositivity as new prognostic factors for a severe COVID-19 course. Int J Mol Sci. (2022) 23(19):11338. doi: 10.3390/ijms231911338
78. Russell CD, Valanciute A, Gachanja NN, Stephen J, Penrice-Randal R, Armstrong SD, et al. Tissue proteomic analysis identifies mechanisms and stages of immunopathology in fatal COVID-19. Am J Respir Cell Mol Biol. (2022) 66:196–205. doi: 10.1165/rcmb.2021-0358OC
79. Lawal B, Kuo YC, Rachmawati Sumitra M, Wu ATH, Huang HS. Identification of a novel immune-inflammatory signature of COVID-19 infections, and evaluation of pharmacokinetics and therapeutic potential of RXn-02, a novel small-molecule derivative of quinolone. Comput Biol Med. (2022) 148:105814. doi: 10.1016/j.compbiomed.2022.105814
80. Liu Z, Li J, Chen D, Gao R, Zeng W, Chen S, et al. Dynamic interleukin-6 level changes as a prognostic indicator in patients with COVID-19. Front Pharmacol. (2020) 11:1093. doi: 10.3389/fphar.2020.01093
81. Huang J, Wang Y, Zha Y, Zeng X, Li W, Zhou M. Transcriptome analysis reveals hub genes regulating autophagy in patients with severe COVID-19. Front Genet. (2022) 13:908826. doi: 10.3389/fgene.2022.908826
82. Islam T, Rahman MR, Aydin B, Beklen H, Arga KY, Shahjaman M. Integrative transcriptomics analysis of lung epithelial cells and identification of repurposable drug candidates for COVID-19. Eur J Pharmacol. (2020) 887:173594. doi: 10.1016/j.ejphar.2020.173594
83. Park JH, Lee HK. Re-analysis of single cell transcriptome reveals that the NR3C1-CXCL8-neutrophil axis determines the severity of COVID-19. Front Immunol. (2020) 11:2145. doi: 10.3389/fimmu.2020.02145
84. Giannella A, Riccetti S, Sinigaglia A, Piubelli C, Razzaboni E, Di Battista P, et al. Circulating microRNA signatures associated with disease severity and outcome in COVID-19 patients. Front Immunol. (2022) 13:968991. doi: 10.3389/fimmu.2022.968991
85. Molinero M, Benítez ID, González J, Gort-Paniello C, Moncusí-Moix A, Rodríguez-Jara F, et al. Bronchial aspirate-based profiling identifies microRNA signatures associated with COVID-19 and fatal disease in critically ill patients. Front Med (Lausanne). (2021) 8:756517. doi: 10.3389/fmed.2021.756517
86. Teodori L, Sestili P, Madiai V, Coppari S, Fraternale D, Rocchi MBL, et al. MicroRNAs bioinformatics analyses identifying HDAC pathway as a putative target for existing anti-COVID-19 therapeutics. Front Pharmacol. (2020) 11:582003. doi: 10.3389/fphar.2020.582003
87. Matarese A, Gambardella J, Sardu C, Santulli G. miR-98 regulates TMPRSS2 expression in human endothelial cells: key implications for COVID-19. Biomedicines. (2020) 8(11):462. doi: 10.3390/biomedicines8110462
88. Luo Y, Wen Y, Arevalo Molina AB, Dahal P, Leuprecht L, Bsrat M. Macrophage activation syndrome, glomerulonephritis, pericarditis, and retinal vasculitis as initial presentation of systemic lupus erythematosus. Case Rep Med. (2018) 2018:5979386. doi: 10.1155/2018/5979386
89. Tono T, Takayama Y, Ogawa E, Hoshiyama T, Tono N, Aihara S, et al. [An adult case of acute pericarditis and glomerulonephritis associated with human parvovirus B19 infection]. Kansenshogaku Zasshi. (2013) 87:752–5. doi: 10.11150/kansenshogakuzasshi.87.752
90. Sekar A, Campbell R, Tabbara J, Rastogi P. ANCA glomerulonephritis after the Moderna COVID-19 vaccination. Kidney Int. (2021) 100:473–4. doi: 10.1016/j.kint.2021.05.017
91. Park HJ, An WS, Rha SH, Kim SE, Lee SM. Minimal change glomerulonephritis following the second dose of the Moderna COVID-19 vaccine. Qjm. (2022) 115:490–1. doi: 10.1093/qjmed/hcac094
92. Biliotti E, Caioli A, Sorace C, Lionetti R, Milozzi E, Taibi C, et al. Humoral Immune Response after COVID-19 mRNA Vaccination in Patients with Liver Cirrhosis: A Prospective Real-Life Single Center Study. Biomedicines. (2023) 11(5):1320. doi: 10.3390/biomedicines11051320
93. Shteinshnaider M, Gorelik O, Almoznino-Sarafian D, Alon I, Modai D, Golik A, et al. Constrictive pericarditis in a patient with Turner's syndrome and liver cirrhosis: coincidence or part of the syndrome? J Clin Gastroenterol. (2006) 40:451–2. doi: 10.1097/00004836-200605000-00019
94. Yao M, Huang X, Guo Y, Zhao JV, Liu Z. Disentangling the common genetic architecture and causality of rheumatoid arthritis and systemic lupus erythematosus with COVID-19 outcomes: Genome-wide cross trait analysis and bidirectional Mendelian randomization study. J Med Virol. (2023) 95:e28570. doi: 10.1002/jmv.28570
95. Favalli EG, Ingegnoli F, De Lucia O, Cincinelli G, Cimaz R, Caporali R. COVID-19 infection and rheumatoid arthritis: Faraway, so close! Autoimmun Rev. (2020) 19:102523. doi: 10.1016/j.autrev.2020.102523
96. Thould AK. Constrictive pericarditis in rheumatoid arthritis. Ann Rheum Dis. (1986) 45:89–94. doi: 10.1136/ard.45.2.89
97. Schatz A, Trankle C, Yassen A, Chipko C, Rajab M, Abouzaki N, et al. Resolution of pericardial constriction with Anakinra in a patient with effusive-constrictive pericarditis secondary to rheumatoid arthritis. Int J Cardiol. (2016) 223:215–6. doi: 10.1016/j.ijcard.2016.08.131
98. Damian F, Ana M, Jamileh S, Zheng P, Charlie A, David M, et al. Global prevalence and burden of depressive and anxiety disorders in 204 countries and territories in 2020 due to the COVID-19 pandemic. Lancet. (2021) 398:1700–12. doi: 10.1016/s0140-6736(21)02143-7
99. Mazza MG, De Lorenzo R, Conte C, Poletti S, Vai B, Bollettini I, et al. Anxiety and depression in COVID-19 survivors: Role of inflammatory and clinical predictors. Brain Behav Immun. (2020) 89:594–600. doi: 10.1016/j.bbi.2020.07.037
100. Shader RI. COVID-19, interferons, and depression: A commentary. Psychiatry Res. (2020) 291:113198. doi: 10.1016/j.psychres.2020.113198
101. Baranova A, Cao H, Zhang F. Severe COVID-19 increases the risk of schizophrenia. Psychiatry Res. (2022) 317:114809. doi: 10.1016/j.psychres.2022.114809
102. Jia F, Dong C, Guo H, Liu X, Zheng X, Wang L, et al. Psychiatric symptoms before and after COVID-19 vaccination: A cohort study of hospitalized schizophrenia patients. Asian J Psychiatr. (2022) 78:103319. doi: 10.1016/j.ajp.2022.103319
103. De Berardis D, Campanella D, Serroni N, Rapini G, Olivieri L, Fornaro M, et al. Clozapine-related pericarditis during titration phase in a patient with resistant schizophrenia and concomitant valproate treatment: a case report. J Clin Psychopharmacol. (2014) 34:649–51. doi: 10.1097/jcp.0000000000000179
104. Sparks JA, Tedeschi SK. Improving COVID-19 vaccine immunogenicity by interrupting methotrexate treatment. Lancet Respir Med. (2022) 10:813–5. doi: 10.1016/s2213-2600(22)00224-7
105. Abhishek A, Boyton RJ, Peckham N, McKnight Á, Coates LC, Bluett J, et al. Effect of a 2-week interruption in methotrexate treatment versus continued treatment on COVID-19 booster vaccine immunity in adults with inflammatory conditions (VROOM study): a randomised, open label, superiority trial. Lancet Respir Med. (2022) 10:840–50. doi: 10.1016/s2213-2600(22)00186-2
106. Capecchi PL, Lazzerini PE, Volterrani L, Mazzei MA, Rossetti B, Zanelli G, et al. Antirheumatic agents in covid-19: is IL-6 the right target? Ann Rheum Dis. (2021) 80:e2. doi: 10.1136/annrheumdis-2020-217523
107. Horby P, Lim WS, Emberson JR, Mafham M, Bell JL, Linsell L, et al. Dexamethasone in hospitalized patients with covid-19. N Engl J Med. (2021) 384:693–704. doi: 10.1056/NEJMoa2021436
108. Bouadma L, Mekontso-Dessap A, Burdet C, Merdji H, Poissy J, Dupuis C, et al. High-dose dexamethasone and oxygen support strategies in intensive care unit patients with severe COVID-19 acute hypoxemic respiratory failure: the COVIDICUS randomized clinical trial. JAMA Intern Med. (2022) 182:906–16. doi: 10.1001/jamainternmed.2022.2168
109. Lutman C, Finocchiaro G, Abate E, Milo M, Morassi P, Sinagra G. Purulent pericarditis in rheumatoid arthritis treated with rituximab and methotrexate. J Cardiovasc Med (Hagerstown). (2014) 15:880–2. doi: 10.2459/jcm.0000000000000160
110. Forbat LN, Hancock BW, Gershlick AH. Methotrexate-induced pericarditis and pericardial effusion; first reported case. Postgrad Med J. (1995) 71:244–5. doi: 10.1136/pgmj.71.834.244
111. Lilly LS. Treatment of acute and recurrent idiopathic pericarditis. Circulation. (2013) 127:1723–6. doi: 10.1161/circulationaha.111.066365
112. Hutchinson MR, Lewis SS, Coats BD, Skyba DA, Crysdale NY, Berkelhammer DL, et al. Reduction of opioid withdrawal and potentiation of acute opioid analgesia by systemic AV411 (ibudilast). Brain Behav Immun. (2009) 23:240–50. doi: 10.1016/j.bbi.2008.09.012
113. Diana G, Strollo R, Diana D, Strollo M, Galassi AR, Crea F. Cardiac safety and potential efficacy: two reasons for considering minocycline in place of azithromycin in COVID-19 management. Eur Heart J Cardiovasc Pharmacother. (2021) 7:e53–e4. doi: 10.1093/ehjcvp/pvaa049
114. Klein AL, Imazio M, Cremer P, Brucato A, Abbate A, Fang F, et al. Phase 3 trial of interleukin-1 trap rilonacept in recurrent pericarditis. N Engl J Med. (2021) 384:31–41. doi: 10.1056/NEJMoa2027892
115. Ao G, Wang Y, Li A, Tran C, Yang Q. The effect of canakinumab on clinical outcomes in patients with COVID-19: A meta-analysis. J Infect. (2022) 84:834–72. doi: 10.1016/j.jinf.2022.03.011
116. Kougkas N, Fanouriakis A, Papalopoulos I, Bertsias G, Avgoustidis N, Repa A, et al. Canakinumab for recurrent rheumatic disease associated-pericarditis: a case series with long-term follow-up. Rheumatol (Oxford). (2018) 57:1494–5. doi: 10.1093/rheumatology/key077
117. Gül A, Tugal-Tutkun I, Dinarello CA, Reznikov L, Esen BA, Mirza A, et al. Interleukin-1β-regulating antibody XOMA 052 (gevokizumab) in the treatment of acute exacerbations of resistant uveitis of Behcet's disease: an open-label pilot study. Ann Rheum Dis. (2012) 71:563–6. doi: 10.1136/annrheumdis-2011-155143
118. Audia JP, Yang XM, Crockett ES, Housley N, Haq EU, O'Donnell K, et al. Caspase-1 inhibition by VX-765 administered at reperfusion in P2Y(12) receptor antagonist-treated rats provides long-term reduction in myocardial infarct size and preservation of ventricular function. Basic Res Cardiol. (2018) 113:32. doi: 10.1007/s00395-018-0692-z
Keywords: COVID-19, pericarditis, common genes, bioinformatics analysis, immunology
Citation: Li D, Chen R, Huang C, Zhang G, Li Z, Xu X, Wang B, Li B and Chu X-M (2024) Comprehensive bioinformatics analysis and systems biology approaches to identify the interplay between COVID-19 and pericarditis. Front. Immunol. 15:1264856. doi: 10.3389/fimmu.2024.1264856
Received: 21 July 2023; Accepted: 08 February 2024;
Published: 22 February 2024.
Edited by:
Pietro Ghezzi, University of Urbino Carlo Bo, ItalyReviewed by:
Diego Ulisse Pizzagalli, University of Italian Switzerland, SwitzerlandSabeena Mustafa, King Abdullah International Medical Research Center (KAIMRC), Saudi Arabia
Copyright © 2024 Li, Chen, Huang, Zhang, Li, Xu, Wang, Li and Chu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xian-Ming Chu, chuxianming@qdu.edu.cn; Bing Li, libing_516619@163.com
†ORCID: Xianming Chu, orcid.org/0000-0003-4755-162X