key: cord-1046895-mymkul48
authors: Xia, Jing; Chen, Shuhan; Li, Yaping; Li, Hua; Gan, Minghong; Wu, Jiashuo; Prohaska, Clare Colette; Bai, Yang; Gao, Lu; Gu, Li; Zhang, Dongfang
title: Immune Response Is Key to Genetic Mechanisms of SARS-CoV-2 Infection With Psychiatric Disorders Based on Differential Gene Expression Pattern Analysis
date: 2022-02-04
journal: Front Immunol
DOI: 10.3389/fimmu.2022.798538
sha: 14413b14aaa234623d6c827e780e0defe7724fe2
doc_id: 1046895
cord_uid: mymkul48

Existing evidence demonstrates that coronavirus disease 2019 (COVID-19) leads to psychiatric illness, despite its main clinical manifestations affecting the respiratory system. People with mental disorders are more susceptible to COVID-19 than individuals without coexisting mental health disorders, with significantly higher rates of severe illness and mortality in this population. The incidence of new psychiatric diagnoses after infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is also remarkably high. SARS-CoV-2 has been reported to use angiotensin-converting enzyme-2 (ACE2) as a receptor for infecting susceptible cells and is expressed in various tissues, including brain tissue. Thus, there is an urgent need to investigate the mechanism linking psychiatric disorders to COVID-19. Using a data set of peripheral blood cells from patients with COVID-19, we compared this to data sets of whole blood collected from patients with psychiatric disorders and used bioinformatics and systems biology approaches to identify genetic links. We found a large number of overlapping immune-related genes between patients infected with SARS-CoV-2 and differentially expressed genes of bipolar disorder (BD), schizophrenia (SZ), and late-onset major depressive disorder (LOD). Many pathways closely related to inflammatory responses, such as MAPK, PPAR, and TGF-β signaling pathways, were observed by enrichment analysis of common differentially expressed genes (DEGs). We also performed a comprehensive analysis of protein–protein interaction network and gene regulation networks. Chemical–protein interaction networks and drug prediction were used to screen potential pharmacologic therapies. We hope that by elucidating the relationship between the pathogenetic processes and genetic mechanisms of infection with SARS-CoV-2 with psychiatric disorders, it will lead to innovative strategies for future research and treatment of psychiatric disorders linked to COVID-19.

The novel coronavirus disease 2019 (COVID-19) caused by the highly transmissible severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the biggest health problem facing the world today (1) . Its severe infectivity and fatality rate prompted the World Health Organization (WHO) to declare it a pandemic in March 2020. COVID-19 is a threat to human health and is wreaking havoc on other areas of daily life, such as social interactions and the global economy (2, 3) . Globally, by August 2021, there were about 210 million confirmed cases of COVID-19, including 4.4 million deaths (4) . There is some evidence that the determining factors of the severity of infection with SARS-CoV-2 are more possibly related to host factors but not viral genetic variation (5) , which has been confirmed by several clinical studies. Pre-existing conditions, such as lung disease, cardiovascular disease, kidney disease, type 2 diabetes, hypertension, and psychiatric illness, significantly increase the likelihood of severe symptoms as well as mortality in patients with COVID-19 (5) (6) (7) (8) (9) (10) . In a recent study, the prevalence of schizophrenia was 3.6% in 7,341 patients with COVID-19, much higher than the 0.66% prevalence in the general population (11) (12) (13) . Previous studies have shown that people with mental illness are two to seven times more likely to die after a respiratory infection, and severe mental health disorders have been shown to be a risk factor for increased mortality from COVID-19 (13) (14) (15) . Furthermore, extensive studies indicate that while COVID-19 is damaging to the respiratory system, it can also have negative mental health effects. Patients with COVID-19 initially show symptoms of headache, fatigue, unsteady gait, with severity depending on age, obesity, mental status, and other health-related conditions (16) . COVID-19 may be further complicated by neurological conditions, including polyneuritis, Guillain-Barre syndrome, and encephalomyelitis (17) . Due to the fact that preexisting psychiatric disorders elevate the risk of severe illness and death in COVID-19, it is necessary to explore how these disorders are related to COVID-19 infection. Hopefully, future treatments for SARS-CoV-2 infection can ameliorate these complications by modulating pathways between the disease and COVID-19.

COVID-19 is caused by SARS-CoV-2, which belongs to the coronavirus family (15, 18) . Its pathogenicity is effected via the angiotensin-converting enzyme-2 receptor (ACE2) (19) (20) (21) . Studies have shown that ACE2 greatly promotes the efficiency of SARS-CoV-2 replication. The ability of SARS-CoV-2 to gain intracellular entry depends on the presence of ACE2 protein (22) , leading to active infection (23) . The ACE2 receptor is present in the lungs, cardiovascular system, gut, kidney, liver, brain, and adipose tissue (24, 25) , as well as the medulla oblongata, glial cells, and central nervous system (17, 19, 26) . ACE2 has a variety of functions in the central nervous system, with significant regulatory roles in stress response, anxiety, cognition, brain damage, and neurogenesis (27, 28) . Thus, ACE2 may be the key to potentially shared mechanisms underlying COVID-19 and psychiatric disorders.

Herein, we explored the cross talk between COVID-19 and bipolar disorder (BD), schizophrenia (SZ), and late-onset major depressive disorder (LOD). BD is characterized by alternating hypomanic and manic episodes and is among the most serious mental health disorders globally, typically associated with an early age of onset (29, 30) . SZ is a multifactorial neuropsychiatric disorder attributed to a developmental disorder of the brain, with symptoms usually occurring during adolescence through young adulthood (31, 32) . LOD refers to major depressive disorder that manifests typically after the age of 50 with depressive symptoms and characteristics different from those of early-onset patients (33, 34) . We used data sets for COVID-19 and each psychiatric disorder to obtain common differentially expressed genes (DEGs) and common immune DEGs. Through the analysis of these DEGs, we explored the pathway and gene expression ontologies to study their biological significance. Proteinprotein, transcription, and post-transcription expression interaction networks were built with common DEGs. Chemical-protein interaction and drug prediction for comorbidities were screened.

Bioinformatics along with system biology methods is commonly used to study comorbidity complexities (35) (36) (37) . We collected data sets from a public data resource. To determine the DEGs in each data set, gene expression analysis was performed. Common DEGs between COVID-19 and BD, SZ, and LOD were identified. These common DEGs were then used to establish proteinprotein interaction (PPI) networks and identify cell signaling pathways and functional gene ontology (GO) to investigate the biological significance. We used common DEGs to create gene regulatory networks (GRN): DEG-microRNA (miRNA) network and DEG-transcription factor (TF) network. Lastly, we identified protein-chemical interactions and screened suggested drugs. The overall workflow is depicted in Table 1 .

We collected selected human gene expression raw data sets for BD, SZ, and LOD from the Gene Expression Omnibus (GEO) of the National Center for Biotechnology Information (NCBI). The GEO accession number of the COVID-19 data set is GSE163151, and the samples of this data set were obtained from whole peripheral blood cell samples of individuals with COVID-19 along with healthy control subjects. Samples of selected psychiatric disorders data sets having association numbers GSE46416 (BD patients and healthy controls), GSE54913 (SZ patients and healthy controls), and GSE76826 (LOD patients and healthy controls) were collected by analyzing the peripheral blood of the patients and healthy controls.

For the purpose of determining the DEGs in data sets, we used the limma R package, to identify DEGs (42) . The Benjamini-Hochberg false discovery rate method can effectively discover genes which are statistically significant and limit false positives. Genes exhibiting an adjusted P of <0.05 along with a log 2 foldchange of ≥1 are determined as DEGs. DEGs of all data sets were identified in this study, including one for SARS-CoV-2 and three for psychiatric disorders. The shared DEGs in psychiatric illnesses and SARS-CoV-2 were also screened. Additionally, DEGs which are shared in SARS-CoV-2 immune and psychiatric illness were determined by ImmPort (http://www. immport.org/).

Gene set enrichment analysis, which can clarify potential biological mechanisms, consisted of pathway analysis and gene ontology analysis (9, 39) . Pathways with the adjusted P-value <0.05 were eligible. We selected pathway and ontology terms by KOBAS-i (http://kobas.cbi.pku.edu.cn/) (39) . Pathways were obtained from four data resources including Reactome, KEGG, Panther, and BioCyc in KOBAS-i.

The PPI (38) is commonly used to explain the association between diseases, so we utilized DEGs to create PPI networks. We used version 11.5 of STRING (https://www.string-db.org/) to insert common DEGs to generate PPI networks. The confidence score was used for the PPI network using the STRING platform with a medium confidence score of 0.400. We performed the Markov cluster algorithm (MCL) function in STRING to identify gene clusters. For a superior visual representation of the network and for the purpose of identifying hub genes, the PPIs are analyzed through Cytoscape (https://cytoscape.org/). 

We used GRN to explore the underlying mechanisms of disease. DEG-miRNA cross talk networks and DEG-TF cross talk networks from common DEGs were constructed by Network Analyst (40) . TarBase (43) and miRTarBase (44) data resources were employed for DEG-miRNA cross talk networks, and the JASPAR (45) data resource was used in TF-DEG cross talk network analysis. DEG-miRNA and DEG-TF networks were filtered with a betweenness value of 150 and 4,000, respectively, and were illustrated in Cytoscape. Regulation at the posttranscriptional level was identified by these interaction networks.

We used common DEGs between COVID-19 and psychiatric illnesses to construct the chemical-protein interaction networks and a list of suggested therapeutic drugs. The chemical-protein interaction networks from common DEGs were built by the Comparative Toxicogenomics Database (CTD) in Network Analyst (40) . The networks were filtered with a degree value of 10 and were illustrated in Cytoscape. The list of drugs derived by common DEGs was identified by the web-utility L1000FDW (https://maayanlab. cloud/L1000FWD/) (41) , which can determine potential repurposable pharmacological agents. An adjusted P-value (q-value) <0.05 was used as a threshold for significance. Among them, the top 5 drugs which have been marketed are shown in our results.

DEGs were identified for each data set by using adjusted P-values (adj P) <0.05 along with an absolute value of logFC (log foldchange) of >1. Table 1 ). Figure 1 demonstrates notable genes for SZ, BD, and LOD. The red dots in the volcano plots designate genes remarkably upregulated and the blue dots designate remarkably downregulated genes. The top 10 gene names with the lowest adj P-value are represented expressly in the figure. We conducted a comparative assessment to determine DEGs that are shared in COVID-19 with other diseases. Figure 1D illustrates the number DEGs shared in these conditions. Our results reveal that shared DEGs with COVID-19 are more common in SZ and BD but not with LOD. Of particular note is CELSR1, a common differential gene of COVID-19, SZ, and LOD.

Common DEGs for Immune Responses of COVID-19, BD, SZ, and LOD

We identified DEGs from immune response gene sets from patients with COVID-19, BD, SZ, and LOD and identified DEGs common among them ( Figure 2A) , with heatmaps demonstrating the relationship between these common DEGs.

To reflect the relationship between unaffected and disease states, we constructed heatmaps from two different perspectives: adj Pvalue ( Figure 2B ) and the values of log fold-change ( Figure 2C ). 

We identified DEG-related pathways shared by COVID-19 and BD, SZ, and LOD, analyzing the possible roles of these genes in the involved pathways using enrichment analysis combining all DEGs from peripheral blood from patients with COVID-19. After integrating the pathways in these data resources, we plotted the bubble map of the top 25 most significant pathways according to adjusted P, as illustrated in Figure 3 . For example, Pathways Eukaryotic Translation Elongation, Signal Transduction, and Signaling by Activin are the most significantly enriched pathways in BD, SZ, and LOD, respectively ( Figure 3 ). As exhibited in Figure 3A and Supplementary Table 2 , there are many pathways associated with the pathogenic mechanisms of various viruses, i.e., viral mRNA translation, influenza viral RNA transcription along with replication, human papillomavirus infection, influenza life cycle, and influenza infection in COVID and BD. In addition, many pathways closely related to inflammatory responses were observed. For instance, the PI3K-Akt signaling cascade, PPAR signaling cascade, and MAPK signaling cascade were found to be closely related in both COVID and SZ ( Figure 3B and Supplementary Table 3 ). We also found many pathways related to TGF-b and FGF/FGFR signaling between COVID and LOD ( Figure 3C and Supplementary Table 4 ).

Next, we performed GO pathway analysis using the same common DEGs to explore biological process categories. Figure 4 shows the top 25 GO pathways of COVID-19 relative to BD, SZ, and LOD. The significantly enriched GO terms are identified by high logarithmic value of the adj P-value. Based on Figure 4 , protein binding is a prominent ontology pathway between COVID-19 and BD. This GO pathway was also significantly enriched in SZ and LOD (Supplementary Tables 5-7) .

We built a PPI network with the shared DEGs between COVID-19 and BD, SZ, and LOD. As illustrated in Figure 5A , the PPI network is created by shared DEGs in whole blood along with immune response cells from individuals infected with COVID-19, as well as patients with BD, SZ, and LOD. The purple circles indicate proteins engaged in each of the psychiatric disorders, as well as COVID-19. The cross talk between proteins is designated by the edges. The size and color depth of the circles designate the degree of protein interrelatedness (number of connecting edges), with stronger relationships indicated by more connections with hub proteins (46) . Figure 5B is a network of hub proteins and 

We obtained a clear network of interactions between DEG and miRNA by using TarBase and miRTarBase data resources, as shown in Figure 6 . Circles in the figure designate DEGs, and diamonds indicate miRNAs. One circle is connected with more than one diamond, with the node degree referring to connection numbers between the node and other nodes in the network. Circles and diamonds in Figure 6 are sized and colored to indicate the degree of the nodes, with the larger the size and the darker the color of a node, the more nodes connected to it, and such nodes are considered to have a critical role in the network. We found that in the BD-COVID-19 interaction ( Figure 6A ), genes RRM2, IFNAR2, and SLC7A11 and miRNAs hsa-mir-30a-5p, hsa-mir-93-5p, and hsa-mir-192-5p have more significant effects. As displayed in Figures 6D-F , we also constructed an interaction network between DEGs and transcription factors (TF) through the JASPER data resource, using the same DEGs as the DEG-miRNA network. Circles and triangles designate DEGs and TFs, respectively. Likewise, the size of the circles and triangles, as well as their shades of colors, depends on the degree of the nodes. Figure 6E illustrates the DEG-TF network of BD-COVID-19 demonstrating that BDNF, RUNX1T1, CDK1, TOP2B, and TEDDM1 constitute the overexpressed genes, and FOXC1 and GATA2 are among the significant TFs.

A chemical-protein interaction network is an important research tool for understanding the function of proteins, which is helpful for advancing drug discovery. We created a chemical-protein cross talk network using shared DEGs for COVID-19 with BD, SZ, and LOD using peripheral blood and immune cells, shown in Figure 7 . The squares in Figure 7 designate chemicals that influence the expression of the gene, and round nodes represent protein-coding genes. The degree of the node reflects the number of connections the node harbors with other nodes in the network. The color from yellow to purple and node size indicate the degree of node from low to high. Figure 7A illustrates the chemical-protein cross talk network between BD and COVID-19. Significant proteins in this network include SLC7A11, TRIB3, RRM2, EGF, and NUSAP1. Figure 7B illustrates the chemical-protein network between COVID-19 and SZ. The overexpressed proteins in this cross talk include CDK1, PPARG, EDN1, CXCL1, and CCNE1. Figure 7C shows the overexpressed proteins between LOD and COVID-19, including SRXN1, SPRY1, EPAS1, and CYP4F3. Valproic acid is the most highly enriched chemical observed in BD, SZ, and LOD networks with COVID-19. Table 2 shows the predicted potential drugs that could treat COVID-19 patients with BD, SZ, or LOD, arranged based on similarity scores. The top drug with the most anti-similar signature for BD was capsaicin. The top drug for SZ was racecadotril. Erythromycin was the most significant drug for LOD.

The focus of this paper is to explore the pathogenetic processes and genetic mechanisms between COVID-19 and psychiatric disorders from different biological perspectives. Recently, transcriptome data have often been used in studies related to comorbidities of COVID-19 (9, 10) . Through a comprehensive analysis of multiple transcriptomic data, we found that the most shared DEGs were found between COVID-19 with BD and SZ, while only 28 shared DEGs were found between COVID-19 and LOD. Only one gene, Cadherin EGF LAG Seven-Pass G-Type Receptor 1 (CELSR1), has been found to be shared between COVID-19 and SZ and LOD. CELSR1 is related to B-cell proliferation, cell-cell adhesion, and central nervous system development (48) . Many clinical manifestations are associated with these processes. Immune cells in COVID-19 patients produce inflammatory cytokines such as IL-2 and IL-6 to stimulate B-cell proliferation (49) . Infection with SARS-CoV-2 also causes the increased expression of endothelial cell adhesion molecules, inducing dysfunction of the coagulation cascade (50) . We also analyzed immune responses of individuals with COVID-19 with BD, SZ, and LOD to understand how COVID-19 affects the immune system response process in patients with psychiatric disorders. One hundred and fifteen DEGs from our COVID-19 immune response data set determined by ImmPort were identified. COVID-19 immune response and SZ share one DEG (activin A receptor type II-like 1, or ACVRL1) that is upregulated and two DEGs (CXCR2, CXCL1) that are both downregulated, all of which are strongly associated with the pathogenesis of COVID-19. ACVRL1 belongs to the family of TGF-b I receptors, and its abnormal changes significantly affect the regulation of the TGF-b pathway, contributing to the development of pulmonary hypertension (51, 52) . The imbalance of TGF-b biological activity predisposes individuals to infection with COVID-19 (53) . C-X-C chemokine receptor type 2 protein (CXCR2) and its ligand C-X-C motif chemokine 1 (CXCL1) are involved in various inflammatory diseases (54) (55) (56) (57) (58) . In COVID-19 patients, CXCL1 binds to CXCR2 to promote immune cell chemotaxis, causing a cytokine storm that worsens the symptoms of patients (59, 60) . One upregulated DEG (EGF) and one downregulated DEG (IFNAR2) have been uncovered in COVID-19 responses in BD. The Interferon receptor 2 (IFNAR2) gene is involved in the synthesis of interferon which has antiviral effects, with low gene expression levels in patients with severe COVID-19 (61) . In order to illustrate the relationship of mechanisms between COVID-19 and BD, SZ, and LOD, we performed enrichment analyses by using DEGs to discover the biological processes of pathways and gene ontology. The common DEGs used for enrichment analysis were derived from the patients with BD, SZ, and LOD and peripheral blood and immune responses from patients infected with COVID-19. Pathways enriched by genes in COVID-19 and BD are found to be related to viral infections, including the influenza virus and human papillomavirus. During the life cycle of the virus, viral RNA transcription, replication, and viral mRNA translation are highly enriched. This is likely to be the main mechanism by which BD patients are susceptible to SARS-CoV-2 (62). Alternatively, we observed that metabolism of proteins and selenoamino acid metabolism are enriched to varying degrees. We identified multiple cell signaling pathways including PI3K-Akt, MAPK, PPAR, and p53 signaling important between COVID-19 and SZ. The PI3K-Akt signaling pathway responds to extracellular signals and promotes metabolism, proliferation, cell survival, growth, and angiogenesis (63) (64) (65) . The MAPK signaling pathway is also highly enriched, which is related to the regulation of cell growth, differentiation, stress adaptation to the environment, inflammation, and other remarkable cellular physiological/pathological processes (66) (67) (68) (69) . PPAR signaling modulates the expression of genes that regulate lipid metabolism, fat formation, and maintenance of metabolic homeostasis along with inflammation, and induces anticancer effects in numerous human tumors (70) (71) (72) (73) (74) . The p53 signaling pathway is also identified, which plays an important role in restoring the proliferation of cells that have DNA damage (75, 76) . Moreover, some potential pathways have been noted such as cytokine-cytokine receptor interaction and signaling by receptor tyrosine kinases. TGF-b and fibroblast growth factors (FGF) and fibroblast growth factor receptor (FGFR) signaling pathways are very significant for LOD and COVID-19. Members of the TGF-b family are expressed in tissue-specific forms at different time points and, therefore, play an indispensable role in the development, homeostasis, and repair of most tissues in the body (77) (78) (79) . All immune cells, entailing B cells, T cells, dendritic cells, and macrophages, secrete TGF-b, which in turn negatively modulates the proliferation, differentiation, and activation of immune cells through other cytokines (80) (81) (82) (83) (84) . In addition, its signal disturbance has been linked to autoimmune diseases, inflammation, and cancer (85) . FGFs bind to FGFRs, activating downstream signaling pathways, which play an important role in biological processes promoting or inhibiting division (86) (87) (88) . TGF-b and FGF/FGFR signaling pathways influence biological processes through their interactions (89) (90) (91) .

The gene ontology analysis reveals protein binding, exocytosis, positive regulation of neurotransmitter secretion, and protein targeting-related pathways for BD. Neutrophil chemotaxis and negative regulation of cell adhesion pathways are the key ontological biological processes in COVID-19 and SZ. Neutrophils play an important role in both sterile and infective inflammatory responses. As the first-line cells of inflammatory responses, neutrophil chemotaxis is a significant manifestation of the anti-inflammatory ability of the body (92) . Cell adhesion plays a pivotal role in maintaining the stability of bronchial epithelium mediating airway inflammation in acute and chronic lung inflammation (93) (94) (95) (96) . Increased chemotaxis of inflammatory cells and upregulation of cellular adhesion molecule expression are also present in the response of the body to COVID-19 (97, 98) . In addition to lung development, cell differentiation, and regulation of timing of cell differentiation, there are also activin receptor signaling pathway, protein binding, ATP binding, and apoptotic nuclear change ontology pathways common to LOD and COVID-19.

With the purpose of systematically analyzing the interaction between a large number of proteins and understanding the physiological mechanisms of disease states, we created a PPI network. We found that epidermal growth factor (EGF), which regulates intracellular signaling molecules, was the highest expressed gene in the interactive network. EGF can promote the survival, maturation, and differentiation of neurons, especially midbrain dopaminergic neurons, and its dysfunction may affect the physiological and pathological conditions of BD (99, 100) . A previous study on COVID-19 suggested that EGF has a potential role in lung injury in infection with SARS-CoV-2 (101). Furthermore, cyclin-dependent kinase 1 (Cdk1) has the highest connectivity among hub genes. In response to SARS-CoV-2 infection, the activity of CDK1 and other cycle-related kinases decreased significantly in host cells (102, 103) . It plays an indispensable role in regulating mitosis, cell cycle progression, apoptosis, and cell pluripotency (104, 105) . The cell cycle-related genes have been used as biomarkers in research investigating the diagnosis and treatment of schizophrenia (106) (107) (108) . Many overexpressed genes, such as UBE2C, ATM, RRM2, UHRF1, and CHEK1, are related to immune system function and provide potential mechanisms in SARS-CoV-2 infection (109) (110) (111) (112) (113) . A further in-depth study of these genes may contribute to the clinical diagnosis and treatment of COVID-19 and its comorbidities with psychiatric disorders.

We have found relationships between DEG-miRNA and DEG-TF interactions in COVID-19 and psychiatric disorders. We identified miRNAs with key roles including hsa-mir-30a-5p, hsamir-93-5p, hsa-mir-155-5p, hsa-mir-92a-3p, hsa-mir-335-5p, hsamir-124-3p, and hsa-mir-6849-3p. These miRNAs are important in the pathogenesis of mental illnesses and also play a role in other major diseases. For example, hsa-mir-124-3p plays a remarkable role in the development of colon cancer, nasopharyngeal cancer, and gliomas (114) . Furthermore, our research indicates that transcription factors FOXC1, GATA2, YY1, NR2C2, and E2F4 were overexpressed in common DEGs of psychiatric disorders and COVID-19. Previous studies have suggested that the functions of the FOX and GATA families affect the pathophysiological mechanisms of psychiatric illnesses (115, 116) . Using our chemical-protein interaction networks, we were able to predict drugs by the web-utility L1000FDW to screen candidate drug molecules. Through chemical-protein networks, we found valproic acid, cyclosporine, estradiol, and tretinoin, all of which have potential effects in fighting COVID-19. Valproic acid had shown good clinical value in the treatment of COVID-19 and Alzheimer's disease (117, 118) . Drug prediction results showed that capsaicin, racecadotril, and erythromycin had the strongest effect according to similarity score. In one recent study, capsaicin showed strong activity in treating mental illness (119) . In addition, capsaicin and cefixime both have a strong ability to penetrate the blood-brain barrier (BBB) (120, 121) , and vincristine can also penetrate the BBB to a certain extent (122, 123) . This will provide implications and new ideas for the clinical prevention along with treatment of COVID-19 and its co-occurrence with BD, SZ, and LOD.

To explore the pathogenetic processes and genetic mechanisms of comorbidities of COVID-19 with psychiatric disorders, we identified common DEGs in BD, SZ, and LOD with COVID-19 and its immune response. Common DEGs in psychiatric disorders with COVID-19 and its immune response have been integrated to identify the pathways and genetic ontological terms. We found that TGF-b, MAPK, PPAR, and other immune-related pathways were highly enriched. We created PPI interaction networks to understand the relationship between important proteins in COVID-19 patients and BD, SZ, and LOD. Hub proteins were identified, with CDK1 being the most highly expressed protein. Transcriptional and posttranscriptional analyses were also performed to uncover the involvement of potential miRNAs and transcription factors and to identify important factors for infection with SARS-CoV-2. Finally, we created chemical-protein networks and drug prediction to identify a number of potentially drug molecules. The genes, pathways, and networks identified in this study are in accordance with previous studies, implying that there are multiple mechanisms of comorbidity between COVID-19 and BD, SZ, and LOD. We hope our work will lead to further investigation and provide new ideas and treatment strategies in the clinical fight against COVID-19.

The data sets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

The Architecture of the SARS-CoV-2 RNA Genome Inside Virion

Neurological Insights of COVID-19 Pandemic

No Participant Left Behind: Conducting Science During COVID-19

COVID-19 in People With Schizophrenia: Potential Mechanisms Linking Schizophrenia to Poor Prognosis

COVID-19 and Angiotensin-Converting Enzyme Inhibitors and Angiotensin Receptor Blockers: What Is the Evidence?

Network-Based Identification Genetic Effect of SARS-CoV-2 Infections to Idiopathic Pulmonary Fibrosis (IPF) Patients

Updates on What ACS Reported: Emerging Evidences of COVID-19 With Nervous System Involvement

Predicting COVID-19-Comorbidity Pathway Crosstalk-Based Targets and Drugs: Towards Personalized COVID-19 Management

A Novel Multi-Omics-Based Highly Accurate Prediction of Symptoms, Comorbid Conditions, and Possible Long-Term Complications of COVID-19

Effect of Underlying Comorbidities on the Infection and Severity of COVID-19 in Korea: A Nationwide Case-Control Study

Annual Prevalence and Incidence of Schizophrenia and Similar Psychotic Disorders in the Republic of Korea: A National Health Insurance Data-Based Study

Increased Risk of COVID-19 Infection and Mortality in People With Mental Disorders: Analysis From Electronic Health Records in the United States

Association Between Mental Health Disorders and Mortality Among Patients With COVID-19 in 7 Countries: A Systematic Review and Meta-Analysis

Severe Mental Illness is Associated With Increased Mortality and Severe Course of COVID-19

Potential Neurological Symptoms of COVID-19

Evidence of the COVID-19 Virus Targeting the CNS: Tissue Distribution, Host-Virus Interaction, and Proposed Neurotropic Mechanisms

Characteristics and Prognosis of COVID-19 in Patients With COPD

Angiotensin-Converting Enzyme 2: SARS-CoV-2 Receptor and Regulator of the Renin-Angiotensin System: Celebrating the 20th Anniversary of the Discovery of ACE2

ACE2 Inhibits Breast Cancer Angiogenesis via Suppressing the VEGFa/VEGFR2/ERK Pathway

A Pneumonia Outbreak Associated With a New Coronavirus of Probable Bat Origin

Angiotensin-Converting Enzyme 2: A Functional Receptor for SARS Coronavirus

SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor

Tissue Distribution of ACE2 Protein, the Functional Receptor for SARS Coronavirus. A First Step in Understanding SARS Pathogenesis

Expression of the SARS-CoV-2 Cell Receptor Gene ACE2 in a Wide Variety of Human Tissues

Angiotensin-Converting Enzyme 2 Overexpression Improves Central Nitric Oxide-Mediated Sympathetic Outflow in Chronic Heart Failure

ACE2 in Brain Physiology and Pathophysiology: Evidence From Transgenic Animal Models

ACE2 Activation Protects Against Cognitive Decline and Reduces Amyloid Pathology in the Tg2576 Mouse Model of Alzheimer's Disease

Global Burden of Disease in Young People Aged 10-24 Years: A Systematic Analysis

Grand Challenges in Global Mental Health

The Neurodevelopmental Hypothesis of Schizophrenia

Non-Genetic Factors in Schizophrenia

The Neurobiology of Depression in Later-Life: Clinical, Neuropsychological, Neuroimaging and Pathophysiological Features

Early and Late Onset Depression in Young and Middle Aged Adults: Differential Symptomatology, Characteristics and Risk Factors?

Bioinformatics and System Biology Approach to Identify the Influences of COVID-19 on Cardiovascular and Hypertensive Comorbidities

Discovering Common Pathogenetic Processes Between COVID-19 and Diabetes Mellitus by Differential Gene Expression Pattern Analysis

Genetic Mechanisms of COVID-19 and Its Association With Smoking and Alcohol Consumption

Protein-Protein Association Networks With Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets

KOBAS-I: Intelligent Prioritization and Exploratory Visualization of Biological Functions for Gene Enrichment Analysis

NetworkAnalyst for Statistical, Visual and Network-Based Meta-Analysis of Gene Expression Data

L1000FWD: Fireworks Visualization of Drug-Induced Transcriptomic Signatures

Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies

TarBase 6.0: Capturing the Exponential Growth of miRNA Targets With Experimental Support

Updates to the Experimentally Validated microRNA-Target Interaction Database

Update of the Open-Access Database of Transcription Factor Binding Profiles

What Properties Characterize the Hub Proteins of the Protein-Protein Interaction Network of Saccharomyces Cerevisiae?

Cytohubba: Identifying Hub Objects and Sub-Networks From Complex Interactome

Neurosurgical Practice During Coronavirus Disease 2019 (COVID-19) Pandemic

Expansion of Plasmablasts and Loss of Memory B Cells in Peripheral Blood From COVID-19 Patients With Pneumonia

Molecular and Functional Analysis Identifies ALK-1 as the Predominant Cause of Pulmonary Hypertension Related to Hereditary Haemorrhagic Telangiectasia

Transforming Growth Factor-Beta Receptor Mutations and Pulmonary Arterial Hypertension in Childhood

Features of Cytokine Storm Identified by Distinguishing Clinical Manifestations in COVID-19

Potential Roles and Targeted Therapy of the CXCLs/CXCR2 Axis in Cancer and Inflammatory Diseases

Targeted Deletion of CXCR2 in Myeloid Cells Alters the Tumor Immune Environment to Improve Antitumor Immunity

Macrophage/epithelial Cell CCL2 Contributes to Rhinovirus-Induced Hyperresponsiveness and Inflammation in a Mouse Model of Allergic Airways Disease

Smoking Cessation and Cardiovascular Disease Risk Factors: Results From the Third National Health and Nutrition Examination Survey

Inhibition of CXCR2 Profoundly Suppresses Inflammation-Driven and Spontaneous Tumorigenesis

Imbalanced Host Response to SARS-CoV-2 Drives Development of COVID-19

Chemokines and Eicosanoids Fuel the Hyperinflammation Within the Lungs of Patients With Severe COVID-19

Genetic Mechanisms of Critical Illness in COVID-19

Viral and Cellular mRNA Translation in Coronavirus-Infected Cells

PI3K/Akt/mTOR and Ras/Raf/MEK/ERK Signaling Pathways Inhibitors as Anticancer Agents: Structural and Pharmacological Perspectives

The Fibroblast Growth Factor Signaling Pathway

Distinct Functions of AKT Isoforms in Breast Cancer: A Comprehensive Review

Signaling Pathway of MAPK/ERK in Cell Proliferation, Differentiation, Migration, Senescence and Apoptosis

Interactions Between Canonical Wnt Signaling Pathway and MAPK Pathway Regulate Differentiation, Maturation and Function of Dendritic Cells

Cell Signaling and Stress Responses

The P38 MAPK Inhibitors for the Treatment of Inflammatory Diseases and Cancer

PPARs as Nuclear Receptors for Nutrient and Energy Metabolism

The CD36-Pparg Pathway in Metabolic Disorders

Feedback Regulation of Arid5a and Ppar-g2 Maintains Adipose Tissue Homeostasis

Oxidative Stress and Immune Dysfunction in Down Syndrome / Maziana Mahamood

Expression and Function of PPARs in Cancer Stem Cells

The P53 Pathway: What Questions Remain to be Explored?

The Role of P53 in Metabolic Regulation

TGF-b Signaling is Dynamically Regulated During the Alveolarization of Rodent and Human Lungs

Relative Roles of TGF-Beta1 and Wnt in the Systemic Regulation and Aging of Satellite Cell Responses

Psychiatric Disorder Comorbidities of SARS-CoV-2

Dynamic Regulation of Tgf-B Signaling by Tif1g: A Computational Approach

Interplay Between Craniofacial Stem Cells and Immune Stimulus

Molecular Therapy for Disk Degeneration and Pain

The CREB/CRTC2 Pathway Modulates Autoimmune Disease by Promoting Th17 Differentiation

TGF-br Inhibitor SB431542 Restores Immune Suppression Induced by Regulatory B-T Cell Axis and Decreases Tumour Burden in Murine Fibrosarcoma

The Influences of TGF-b1 Upon the Human Adenocarcinoma Cell of Lung A549 and Cellular Immunity

Mechanisms of TGF-b Signaling in Regulation of Cell Growth and Differentiation

FGF Controls Epithelial-Mesenchymal Transitions During Gastrulation by Regulating Cell Division and Apicobasal Polarity

The Ubiquitin Ligase Mlin41 Temporally Promotes Neural Progenitor Cell Maintenance Through FGF Signaling

In Vitro and in Vivo Analyses Reveal Profound Effects of Fibroblast Growth Factor 16 as a Metabolic Regulator

Toward Pathogenesis of Apert Cleft Palate: FGF, FGFR, and Tgfb Genes Are Differentially Expressed in Sequential Stages of Human Palatal Shelf Fusion

FGFR1-WNT-TGF-b Signaling in Prostate Cancer Mouse Models Recapitulates Human Reactive Stroma

SIG-02. Rational Targeting of Pro-Invasive FGFR Signaling in Medulloblastoma

Vitamin C and Immune Function

Expression of the Cell Adhesion Molecule CD146/MCAM in Non-Small Cell Lung Cancer

Dipeptidase-1 Is an Adhesion Receptor for Neutrophil Recruitment in Lungs and Liver

Cell Adhesion Molecules and Potential for Pharmacological Intervention in Lung Inflammation

Differential Contribution of Adhesion Molecules to Th1 and Th2 Cell-Mediated Lung and Bowel Inflammation

Inflammatory Profiles and Clinical Features of COVID-19 Survivors Three Months After Discharge in Wuhan

Elevated Serum Endothelial Cell Adhesion Molecules Expression in COVID-19 Patients

Association Between EGF A61G Polymorphism and Bipolar Disorder in Korean Population

The P107/E2F Pathway Regulates Fibroblast Growth Factor 2 Responsiveness in Neural Precursor Cells

COVID-19: The Potential Treatment of Pulmonary Fibrosis Associated With SARS-CoV-2 Infection

The Global Phosphorylation Landscape of SARS-CoV-2 Infection

Actionable Cytopathogenic Host Responses of Human Alveolar Type 2 Cells to SARS-CoV-2 -ScienceDirect

Cdk1-Dependent Phosphorylation of Cdc13 Coordinates Telomere Elongation During Cell-Cycle Progression

Cell Cycle Regulation of DNA Double-Strand Break End Resection by Cdk1-Dependent Dna2 Phosphorylation

The Cell Cycle-Related Genes as Biomarkers for Schizophrenia

Expression Profiling of Fibroblasts Identifies Cell Cycle Abnormalities in Schizophrenia

High EGFR_1 Inside-Out Activated Inflammation-Induced Motility Through SLC2A1-CCNB2-HMMR-KIF11-NUSAP1-PRC1-UBE2C

miR-34a in Extracellular Vesicles From Bone Marrow Mesenchymal Stem Cells Reduces Rheumatoid Arthritis Inflammation via the Cyclin I/ATM/ATR/p53 Axis

Independent Prognostic Implications of RRM2 in Lung Adenocarcinoma

Uhrf1 and Dnmt1 Loss Induces an Immune Response in Zebrafish Livers Due to Viral Mimicry by Transposable Elements

Combination of CHEK1/2 Inhibition and Ionizing Radiation Results in Abscopal Tumor Response Through Increased Micronuclei Formation

Sulforaphane Promotes Apoptosis, and Inhibits Proliferation and Self-Renewal of Nasopharyngeal Cancer Cells by Targeting STAT Signal Through miRNA-124-3p

Loss of Foxg1 Impairs the Development of Cortical SST-Interneurons Leading to Abnormal Emotional and Social Behaviors

The Integrated Landscape of Causal Genes and Pathways in Schizophrenia

Valproic Acid in Prevention and Treatment of COVID-19

Valproic Acid as a Promising Agent to Combat Alzheimer's Disease

Gut Microbiota Mediates the Preventive Effects of Dietary Capsaicin Against Depression-Like Behavior Induced by Lipopolysaccharide in Mice

Transient Activation of Specific Neurons in Mice by Selective Expression of the Capsaicin Receptor

Pore-Forming Proteins as Mediators of Novel Epigenetic Mechanism of Epilepsy

P-Glycoprotein in the Blood-Brain Barrier of Mice Influences the Brain Penetration and Pharmacological Activity of Many Drugs

Disruption of the Mouse Mdr1a P-Glycoprotein Gene Leads to a Deficiency in the Blood-Brain Barrier and to Increased Sensitivity to Drugs

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