Key Genes of Immunity Associated with Pterygium and Primary Sjögren's Syndrome.

Pterygium and primary Sjögren's Syndrome (pSS) share many similarities in clinical symptoms and ocular pathophysiological changes, but their etiology is unclear. To identify the potential genes and pathways related to immunity, two published datasets, GSE2513 containing pterygium information and GSE176510 containing pSS information, were selected from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) of pterygium or pSS patients compared with healthy control conjunctiva, and the common DEGs between them were analyzed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted for common DEGs. The protein-protein interaction (PPI) network was constructed using the STRING database to find the hub genes, which were verified in clinical samples. There were 14 co-upregulated DEGs. The GO and KEGG analyses showed that these common DEGs were enriched in pathways correlated with virus infection, antigen processing and presentation, nuclear factor-kappa B (NF-κB) and Th17 cell differentiation. The hub genes (IL1R1, ICAM1, IRAK1, S100A9, and S100A8) were selected by PPI construction. In the era of the COVID-19 epidemic, the relationship between virus infection, vaccination, and the incidence of pSS and pterygium growth deserves more attention.

Peripheral Blood Genes Crosstalk between COVID-19 and Sepsis.

Severe coronavirus disease 2019 (COVID-19) has led to a rapid increase in death rates all over the world. Sepsis is a life-threatening disease associated with a dysregulated host immune response. It has been shown that COVID-19 shares many similarities with sepsis in many aspects. However, the molecular mechanisms underlying sepsis and COVID-19 are not well understood. The aim of this study was to identify common transcriptional signatures, regulators, and pathways between COVID-19 and sepsis, which may provide a new direction for the treatment of COVID-19 and sepsis. First, COVID-19 blood gene expression profile (GSE179850) data and sepsis blood expression profile (GSE134347) data were obtained from GEO. Then, we intersected the differentially expressed genes (DEG) from these two datasets to obtain common DEGs. Finally, the common DEGs were used for functional enrichment analysis, transcription factor and miRNA prediction, pathway analysis, and candidate drug analysis. A total of 307 common DEGs were identified between the sepsis and COVID-19 datasets. Protein-protein interactions (PPIs) were constructed using the STRING database. Subsequently, hub genes were identified based on PPI networks. In addition, we performed GO functional analysis and KEGG pathway analysis of common DEGs, and found a common association between sepsis and COVID-19. Finally, we identified transcription factor-gene interaction, DEGs-miRNA co-regulatory networks, and protein-drug interaction, respectively. Through ROC analysis, we identified 10 central hub genes as potential biomarkers. In this study, we identified SARS-CoV-2 infection as a high risk factor for sepsis. Our study may provide a potential therapeutic direction for the treatment of COVID-19 patients suffering from sepsis.

Recognition of Differentially Expressed Molecular Signatures and Pathways Associated with COVID-19 Poor Prognosis in Glioblastoma Patients.

Glioblastoma (GBM) is a type of brain cancer that is typically very aggressive and difficult to treat. Glioblastoma cases have been reported to have increased during COVID-19. The mechanisms underlying this comorbidity, including genomic interactions, tumor differentiation, immune responses, and host defense, are not completely explained. Therefore, we intended to investigate the differentially expressed shared genes and therapeutic agents which are significant for these conditions by using in silico approaches. Gene expression datasets of GSE68848, GSE169158, and GSE4290 studies were collected and analyzed to identify the DEGs between the diseased and the control samples. Then, the ontology of the genes and the metabolic pathway enrichment analysis were carried out for the classified samples based on expression values. Protein-protein interactions (PPI) map were performed by STRING and fine-tuned by Cytoscape to screen the enriched gene module. In addition, the connectivity map was used for the prediction of potential drugs. As a result, 154 overexpressed and 234 under-expressed genes were identified as common DEGs. These genes were found to be significantly enriched in the pathways involved in viral diseases, NOD-like receptor signaling pathway, the cGMP-PKG signaling pathway, growth hormone synthesis, secretion, and action, the immune system, interferon signaling, and the neuronal system. STAT1, CXCL10, and SAMDL were screened out as the top 03 out of the top 10 most critical genes among the DEGs from the PPI network. AZD-8055, methotrexate, and ruxolitinib were predicted to be the possible agents for the treatment. The current study identified significant key genes, common metabolic signaling networks, and therapeutic agents to improve our perception of the common mechanisms of GBM-COVID-19.

Weighted gene co-expression network-based identification of genetic effect of mRNA vaccination and previous infection on SARS-CoV-2 infection.

To investigate the effect conferred by vaccination and previous infection against SARS-CoV-2 infection in molecular level, weighted gene co-expression network analysis was applied to screen vaccination, prior infection and Omicron infection-related gene modules in 46 Omicron outpatients and 8 controls, and CIBERSORT algorithm was used to infer the proportions of 22 subsets of immune cells. 15 modules were identified, where the brown module showed positive correlations with Omicron infection (r = 0.35, P = 0.01) and vaccination (r = 0.62, P = 1 × 10-6). Enrichment analysis revealed that LILRB2 was the unique gene shared by both phosphatase binding and MHC class I protein binding. Pathways including "B cell receptor signaling pathway" and "FcγR-mediated phagocytosis" were enriched in the vaccinated samples of the highly correlated LILRB2. LILRB2 was also identified as the second hub gene through PPI network, after LCP2. In conclusion, attenuated LILRB2 transcription in PBMC might highlight a novel target in overcoming immune evasion and improving vaccination strategies.

Bioinformatic analysis and preliminary validation of potential therapeutic targets for COVID-19 infection in asthma patients.

BACKGROUND: Severe acute respiratory syndrome coronavirus 2 causes coronavirus disease 19 (COVID-19). The number of confirmed cases of COVID-19 is also rapidly increasing worldwide, posing a significant challenge to human safety. Asthma is a risk factor for COVID-19, but the underlying molecular mechanisms of the asthma-COVID-19 interaction remain unclear. METHODS: We used transcriptome analysis to discover molecular biomarkers common to asthma and COVID-19. Gene Expression Omnibus database RNA-seq datasets (GSE195599 and GSE196822) were used to identify differentially expressed genes (DEGs) in asthma and COVID-19 patients. After intersecting the differentially expressed mRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to identify the common pathogenic molecular mechanism. Bioinformatic methods were used to construct protein-protein interaction (PPI) networks and identify key genes from the networks. An online database was used to predict interactions between transcription factors and key genes. The differentially expressed long noncoding RNAs (lncRNAs) in the GSE195599 and GSE196822 datasets were intersected to construct a competing endogenous RNA (ceRNA) regulatory network. Interaction networks were constructed for key genes with RNA-binding proteins (RBPs) and oxidative stress-related proteins. The diagnostic efficacy of key genes in COVID-19 was verified with the GSE171110 dataset. The differential expression of key genes in asthma was verified with the GSE69683 dataset. An asthma cell model was established with interleukins (IL-4, IL-13 and IL-17A) and transfected with siRNA-CXCR1. The role of CXCR1 in asthma development was preliminarily confirmed. RESULTS: By intersecting the differentially expressed genes for COVID-19 and asthma, 393 common DEGs were obtained. GO and KEGG enrichment analyses of the DEGs showed that they mainly affected inflammation-, cytokine- and immune-related functions and inflammation-related signaling pathways. By analyzing the PPI network, we obtained 10 key genes: TLR4, TLR2, MMP9, EGF, HCK, FCGR2A, SELP, NFKBIA, CXCR1, and SELL. By intersecting the differentially expressed lncRNAs for COVID-19 and asthma, 13 common differentially expressed lncRNAs were obtained. LncRNAs that regulated microRNAs (miRNAs) were mainly concentrated in intercellular signal transduction, apoptosis, immunity and other related functional pathways. The ceRNA network suggested that there were a variety of regulatory miRNAs and lncRNAs upstream of the key genes. The key genes could also bind a variety of RBPs and oxidative stress-related genes. The key genes also had good diagnostic value in the verification set. In the validation set, the expression of key genes was statistically significant in both the COVID-19 group and the asthma group compared with the healthy control group. CXCR1 expression was upregulated in asthma cell models, and interference with CXCR1 expression significantly reduced cell viability. CONCLUSIONS: Key genes may become diagnostic and predictive biomarkers of outcomes in COVID-19 and asthma. Video Abstract.

Molecular signatures in the progression of COVID-19 severity.

SARS-CoV-2 is the causative agent of COVID-19 that has infected over 642 million and killed over 6.6 million people around the globe. Underlying a wide range of clinical manifestations of this disease, from moderate to extremely severe systemic conditions, could be genes or pathways differentially expressing in the hosts. It is therefore important to gain insights into pathways involved in COVID-19 pathogenesis and host defense and thus understand the host response to this pathogen at the physiological and molecular level. To uncover genes and pathways involved in the differential clinical manifestations of this disease, we developed a novel gene co-expression network based pipeline that uses gene expression obtained from different SARS-CoV-2 infected human tissues. We leveraged the network to identify novel genes or pathways that likely differentially express and could be physiologically significant in the COVID-19 pathogenesis and progression but were deemed statistically non-significant and therefore not further investigated in the original studies. Our network-based approach aided in the identification of co-expression modules enriched in differentially expressing genes (DEGs) during different stages of COVID-19 and enabled discovery of novel genes involved in the COVID-19 pathogenesis, by virtue of their transcript abundance and association with genes expressing differentially in modules enriched in DEGs. We further prioritized by considering only those enriched gene modules that have most of their genes differentially expressed, inferred by the original studies or this study, and document here 7 novel genes potentially involved in moderate, 2 in severe, 48 in extremely severe COVID-19, and 96 novel genes involved in the progression of COVID-19 from severe to extremely severe conditions. Our study shines a new light on genes and their networks (modules) that drive the progression of COVID-19 from moderate to extremely severe condition. These findings could aid development of new therapeutics to combat COVID-19.

Potential links between COVID-19 and periodontitis: a bioinformatic analysis based on GEO datasets.

BACKGROUND: 2019 Coronavirus disease (COVID-19) is an infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The COVID-19 pandemic has already had a serious influence on human existence, causing a huge public health concern for countries all around the world. Because SARS-CoV-2 infection can be spread by contact with the oral cavity, the link between oral illness and COVID-19 is gaining traction. Through bioinformatics approaches, we explored the possible molecular mechanisms linking the COVID-19 and periodontitis to provide the basis and direction for future research. METHODS: Transcriptomic data from blood samples of patients with COVID-19 and periodontitis was downloaded from the Gene Expression Omnibus database. The shared differentially expressed genes were identified. The analysis of Gene Ontology, Kyoto Encyclopedia of Genesand Genomes pathway, and protein-protein interaction network was conducted for the shared differentially expressed genes. Top 5 hub genes were selected through Maximal Clique Centrality algorithm. Then mRNA-miRNA network of the hub genes was established based on miRDB database, miRTarbase database and Targetscan database. The Least absolute shrinkage and selection operator regression analysis was used to discover possible biomarkers, which were then investigated in relation to immune-related genes. RESULTS: Fifty-six shared genes were identified through differential expression analysis in COVID-19 and periodontitis. The function of these genes was enriched in regulation of hormone secretion, regulation of secretion by cell. Myozenin 2 was identified through Least absolute shrinkage and selection operator regression Analysis, which was down-regulated in both COVID-19 and periodontitis. There was a positive correlation between Myozenin 2 and the biomarker of activated B cell, memory B cell, effector memory CD4 T cell, Type 17 helper cell, T follicular helper cell and Type 2 helper cell. CONCLUSION: By bioinformatics analysis, Myozenin 2 is predicted to correlate to the pathogenesis and immune infiltrating of COVID-19 and periodontitis. However, more clinical and experimental researches are needed to validate the function of Myozenin 2.

Plasma metabolomics and gene regulatory networks analysis reveal the role of nonstructural SARS-CoV-2 viral proteins in metabolic dysregulation in COVID-19 patients.

Metabolomic analysis of blood plasma samples from COVID-19 patients is a promising approach allowing for the evaluation of disease progression. We performed the metabolomic analysis of plasma samples of 30 COVID-19 patients and the 19 controls using the high-performance liquid chromatography (HPLC) coupled with tandem mass spectrometric detection (LC-MS/MS). In our analysis, we identified 103 metabolites enriched in KEGG metabolic pathways such as amino acid metabolism and the biosynthesis of aminoacyl-tRNAs, which differed significantly between the COVID-19 patients and the controls. Using ANDSystem software, we performed the reconstruction of gene networks describing the potential genetic regulation of metabolic pathways perturbed in COVID-19 patients by SARS-CoV-2 proteins. The nonstructural proteins of SARS-CoV-2 (orf8 and nsp5) and structural protein E were involved in the greater number of regulatory pathways. The reconstructed gene networks suggest the hypotheses on the molecular mechanisms of virus-host interactions in COVID-19 pathology and provide a basis for the further experimental and computer studies of the regulation of metabolic pathways by SARS-CoV-2 proteins. Our metabolomic analysis suggests the need for nonstructural protein-based vaccines and the control strategy to reduce the disease progression of COVID-19.

Systems biology approach reveals a common molecular basis for COVID-19 and non-alcoholic fatty liver disease (NAFLD).

BACKGROUND: Patients with non-alcoholic fatty liver disease (NAFLD) may be more susceptible to coronavirus disease 2019 (COVID-19) and even more likely to suffer from severe COVID-19. Whether there is a common molecular pathological basis for COVID-19 and NAFLD remains to be identified. The present study aimed to elucidate the transcriptional alterations shared by COVID-19 and NAFLD and to identify potential compounds targeting both diseases. METHODS: Differentially expressed genes (DEGs) for COVID-19 and NAFLD were extracted from the GSE147507 and GSE89632 datasets, and common DEGs were identified using the Venn diagram. Subsequently, we constructed a protein-protein interaction (PPI) network based on the common DEGs and extracted hub genes. Then, we performed gene ontology (GO) and pathway analysis of common DEGs. In addition, transcription factors (TFs) and miRNAs regulatory networks were constructed, and drug candidates were identified. RESULTS: We identified a total of 62 common DEGs for COVID-19 and NAFLD. The 10 hub genes extracted based on the PPI network were IL6, IL1B, PTGS2, JUN, FOS, ATF3, SOCS3, CSF3, NFKB2, and HBEGF. In addition, we also constructed TFs-DEGs, miRNAs-DEGs, and protein-drug interaction networks, demonstrating the complex regulatory relationships of common DEGs. CONCLUSION: We successfully extracted 10 hub genes that could be used as novel therapeutic targets for COVID-19 and NAFLD. In addition, based on common DEGs, we propose some potential drugs that may benefit patients with COVID-19 and NAFLD.

Interaction-based transcriptome analysis via differential network inference.

Gene-based transcriptome analysis, such as differential expression analysis, can identify the key factors causing disease production, cell differentiation and other biological processes. However, this is not enough because basic life activities are mainly driven by the interactions between genes. Although there have been already many differential network inference methods for identifying the differential gene interactions, currently, most studies still only use the information of nodes in the network for downstream analyses. To investigate the insight into differential gene interactions, we should perform interaction-based transcriptome analysis (IBTA) instead of gene-based analysis after obtaining the differential networks. In this paper, we illustrated a workflow of IBTA by developing a Co-hub Differential Network inference (CDN) algorithm, and a novel interaction-based metric, pivot APC2. We confirmed the superior performance of CDN through simulation experiments compared with other popular differential network inference algorithms. Furthermore, three case studies are given using colorectal cancer, COVID-19 and triple-negative breast cancer datasets to demonstrate the ability of our interaction-based analytical process to uncover causative mechanisms.

Bioinformatics and systems-biology analysis to determine the effects of Coronavirus disease 2019 on patients with allergic asthma.

Background: The coronavirus disease (COVID-19) pandemic has posed a significant challenge for global health systems. Increasing evidence shows that asthma phenotypes and comorbidities are major risk factors for COVID-19 symptom severity. However, the molecular mechanisms underlying the association between COVID-19 and asthma are poorly understood. Therefore, we conducted bioinformatics and systems biology analysis to identify common pathways and molecular biomarkers in patients with COVID-19 and asthma, as well as potential molecular mechanisms and candidate drugs for treating patients with both COVID-19 and asthma. Methods: Two sets of differentially expressed genes (DEGs) from the GSE171110 and GSE143192 datasets were intersected to identify common hub genes, shared pathways, and candidate drugs. In addition, murine models were utilized to explore the expression levels and associations of the hub genes in asthma and lung inflammation/injury. Results: We discovered 157 common DEGs between the asthma and COVID-19 datasets. A protein-protein-interaction network was built using various combinatorial statistical approaches and bioinformatics tools, which revealed several hub genes and critical modules. Six of the hub genes were markedly elevated in murine asthmatic lungs and were positively associated with IL-5, IL-13 and MUC5AC, which are the key mediators of allergic asthma. Gene Ontology and pathway analysis revealed common associations between asthma and COVID-19 progression. Finally, we identified transcription factor-gene interactions, DEG-microRNA coregulatory networks, and potential drug and chemical-compound interactions using the hub genes. Conclusion: We identified the top 15 hub genes that can be used as novel biomarkers of COVID-19 and asthma and discovered several promising candidate drugs that might be helpful for treating patients with COVID-19 and asthma.

Integrated network-based multiple computational analyses for identification of co-expressed candidate genes associated with neurological manifestations of COVID-19.

'Tripartite network' (TN) and 'combined gene network' (CGN) were constructed and their hub-bottleneck and driver nodes (44 genes) were evaluated as 'target genes' (TG) to identify 21 'candidate genes' (CG) and their relationship with neurological manifestations of COVID-19. TN was developed using neurological symptoms of COVID-19 found in literature. Under query genes (TG of TN), co-expressed genes were identified using pair-wise mutual information to genes available in RNA-Seq autopsy data of frontal cortex of COVID-19 victims. CGN was constructed with genes selected from TN and co-expressed in COVID-19. TG and their connecting genes of respective networks underwent functional analyses through findings of their enrichment terms and pair-wise 'semantic similarity scores' (SSS). A new integrated 'weighted harmonic mean score' was formulated assimilating values of SSS and STRING-based 'combined score' of the selected TG-pairs, which provided CG-pairs with properties of CGs as co-expressed and 'indispensable nodes' in CGN. Finally, six pairs sharing seven 'prevalent CGs' (ADAM10, ADAM17, AKT1, CTNNB1, ESR1, PIK3CA, FGFR1) showed linkages with the phenotypes (a) directly under neurodegeneration, neurodevelopmental diseases, tumour/cancer and cellular signalling, and (b) indirectly through other CGs under behavioural/cognitive and motor dysfunctions. The pathophysiology of 'prevalent CGs' has been discussed to interpret neurological phenotypes of COVID-19.

CORN-Condition Orientated Regulatory Networks: bridging conditions to gene networks.

A transcriptional regulatory network (TRN) is a collection of transcription regulators with their associated downstream genes, which is highly condition-specific. Understanding how cell states can be programmed through small molecules/drugs or conditions by modulating the whole gene expression system granted us the potential to amend abnormal cells and cure diseases. Condition Orientated Regulatory Networks (CORN, https://qinlab.sysu.edu.cn/home) is a library of condition (small molecule/drug treatments and gene knockdowns)-based transcriptional regulatory sub-networks (TRSNs) that come with an online TRSN matching tool. It allows users to browse condition-associated TRSNs or match those TRSNs by inputting transcriptomic changes of interest. CORN utilizes transcriptomic changes data after specific conditional treatment in cells, and in vivo transcription factor (TF) binding data in cells, by combining TF binding information and calculations of significant expression alterations of TFs and genes after the conditional treatments, TRNs under the effect of different conditions were constructed. In short, CORN associated 1805 different types of specific conditions (small molecule/drug treatments and gene knockdowns) to 9553 TRSNs in 25 human cell lines, involving 204TFs. By linking and curating specific conditions to responsive TRNs, the scientific community can now perceive how TRNs are altered and controlled by conditions alone in an organized manner for the first time. This study demonstrated with examples that CORN can aid the understanding of molecular pathology, pharmacology and drug repositioning, and screened drugs with high potential for cancer and coronavirus disease 2019 (COVID-19) treatments.

Transcriptome Analysis of Peripheral Blood Mononuclear Cells Response in Patients with Severe COVID-19 Reveals Crucial Genes Regulating Protein Ubiquitination.

BACKGROUND We sought to further our understanding of the biological characteristics underlying severe COVID-19. MATERIAL AND METHODS RNA sequencing (RNA-Seq) analysis was used to evaluate peripheral blood mononuclear cells from 4 patients with severe COVID-19 and 4 healthy controls. We performed gene expression analyses to detect differentially expressed genes (DEGs). Enrichment analyses were performed to identify their molecular processes and signaling pathways, and the protein-protein interaction network was constructed to extract the core gene cluster. The investigation of protein-chemical interactions and regulatory signatures for specific regulatory checkpoints and powerful chemical agents was then conducted for these essential genes. Finally, we used single-cell RNA-Seq analysis from an online platform to show how these genes were expressed differently, depending on the kind of cell. RESULTS A total of 268 DEGs were found. The biological process of protein ubiquitination was later discovered to be highly influenced by the core gene cluster (ITCH, TRIM21, RNF130, FBXO11, UBE2J1, and ASB16) at the transcriptome level. Six transcription factors, FNIC, FOXA1, YY1, GATA2, MET2A, and FOXC1, as well as miRNAs hsa-miR-1-3p and hsa-miR-27a-3p were identified. We found a potent chemical agent, copper sulfate, may regulate protein ubiquitination genes cooperatively, and the genes regulating protein ubiquitination could be expressed highly on the macrophages. CONCLUSIONS Taken together, we suggest that protein ubiquitination is a crucial functional process in patients with severe COVID-19. This study will give a deeper insight into biological characteristics and progression of COVID-19 and facilitate development of novel therapeutics, leading to significant advancements in the COVID-19 pandemic.

Differential expression analysis of mRNAs, lncRNAs, and miRNAs expression profiles and construction of ceRNA networks in PEDV infection.

BACKGROUND: Porcine Epidemic Diarrhea Virus (PEDV) is a coronavirus that seriously affects the swine industry. MicroRNAs and long noncoding RNAs are two relevant non-coding RNAs (ncRNAs) class and play crucial roles in a variety of physiological processes. Increased evidence indicates a complex interaction between mRNA and ncRNA. However, our understanding of the function of ncRNA involved in host-PEDV interaction is limited. RESULTS: A total of 1,197 mRNA transcripts, 539 lncRNA transcripts, and 208 miRNA transcripts were differentially regulated at 24 h and 48 h post-infection. Gene ontology (GO) and KEGG pathway enrichment analysis showed that DE mRNAs and DE lncRNAs were mainly involved in biosynthesis, innate immunity, and lipid metabolism. Moreover, we constructed a miRNA-mRNA-pathway network using bioinformatics, including 12 DE mRNAs, 120 DE miRNAs, and 11 pathways. Finally, the target genes of DE miRNAs were screened by bioinformatics, and we constructed immune-related lncRNA-miRNA-mRNA ceRNA networks. Then, the selected DE genes were validated by qRT-PCR, which were consistent with the results from RNA-Seq data. CONCLUSIONS: This study provides the comprehensive analysis of the expression profiles of mRNAs, lncRNAs, and miRNAs during PEDV infection. We characterize the ceRNA networks which can provide new insights into the pathogenesis of PEDV.

Weighted Gene Co-Expression Network Analysis to Identify Potential Biological Processes and Key Genes in COVID-19-Related Stroke.

The purpose of this research was to explore the underlying biological processes causing coronavirus disease 2019- (COVID-19-) related stroke. The Gene Expression Omnibus (GEO) database was utilized to obtain four COVID-19 datasets and two stroke datasets. Thereafter, we identified key modules via weighted gene co-expression network analysis, following which COVID-19- and stroke-related crucial modules were crossed to identify the common genes of COVID-19-related stroke. The common genes were intersected with the stroke-related hub genes screened via Cytoscape software to discover the critical genes associated with COVID-19-related stroke. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis for common genes associated with COVID-19-related stroke, and the Reactome database was used to annotate and visualize the pathways involved in the key genes. Two COVID-19-related crucial modules and one stroke-related crucial module were identified. Subsequently, the top five genes were screened as hub genes after visualizing the genes of stroke-related critical module using Cytoscape. By intersecting the COVID-19- and stroke-related crucial modules, 28 common genes for COVID-19-related stroke were identified. ITGA2B and ITGB3 have been further identified as crucial genes of COVID-19-related stroke. Functional enrichment analysis indicated that both ITGA2B and ITGB3 were involved in integrin signaling and the response to elevated platelet cytosolic Ca2+, thus regulating platelet activation, extracellular matrix- (ECM-) receptor interaction, the PI3K-Akt signaling pathway, and hematopoietic cell lineage. Therefore, platelet activation, ECM-receptor interaction, PI3K-Akt signaling pathway, and hematopoietic cell lineage may represent the potential biological processes associated with COVID-19-related stroke, and ITGA2B and ITGB3 may be potential intervention targets for COVID-19-related stroke.

Screening the hub genes and analyzing the mechanisms in discharged COVID-19 patients retesting positive through bioinformatics analysis.

BACKGROUND: After encountering COVID-19 patients who test positive again after discharge, our study analyzed the pathogenesis to further assess the risk and possibility of virus reactivation. METHODS: A separate microarray was acquired from the Gene Expression Omnibus (GEO), and its samples were divided into two groups: a "convalescent-RTP" group consisting of convalescent and "retesting positive" (RTP) patients (group CR) and a "healthy-RTP" group consisting of healthy control and RTP patients (group HR). The enrichment analysis was performed with R software, obtaining the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Subsequently, the protein-protein interaction (PPI) networks of each group were established, and the hub genes were discovered using the cytoHubba plugin. RESULTS: In this study, 6622 differentially expressed genes were identified in the group CR, among which RAB11B-AS1, DISP1, MICAL3, PSMG1, and DOCK4 were up-regulated genes, and ANAPC1, IGLV1-40, SORT1, PLPPR2, and ATP1A1-AS1 were down-regulated. 7335 genes were screened in the group HR, including the top 5 up-regulated genes ALKBH6, AMBRA1, MIR1249, TRAV18, and LRRC69, and the top 5 down-regulated genes FAM241B, AC018529.3, AL031963.3, AC006946.1, and FAM149B1. The GO and KEGG analysis of the two groups revealed a significant enrichment in immune response and apoptosis. In the PPI network constructed, group CR and group HR identified 10 genes, respectively, and TP53BP1, SNRPD1, and SNRPD2 were selected as hub genes. CONCLUSIONS: Using the messenger ribonucleic acid (mRNA) expression data from GSE166253, we found TP53BP1, SNRPD1, and SNRPD2 as hub genes in RTP patients, which is vital to the management and prognostic prediction of RTP patients.

Deregulation of ceRNA Networks in Frontal Cortex and Choroid Plexus of Brain during SARS-CoV-2 Infection Aggravates Neurological Manifestations: An Insight from Bulk and Single-Cell Transcriptomic Analyses.

Although transcriptomic studies of SARS-CoV-2-infected brains have depicted variability in gene expression, the landscape of deregulated cell-specific regulatory circuits has not been elucidated yet. Hence, bulk and single-cell RNA-seq data are analyzed to gain detailed insights. Initially, two ceRNA networks with 19 and 3 differentially expressed (DE) hub lncRNAs are reconstructed in SARS-CoV-2 infected Frontal Cortex (FC) and Choroid Plexus (CP), respectively. Functional and pathway enrichment analyses of downstream mRNAs of deregulated ceRNA axes demonstrate impairment of neurological processes. Mapping of hub lncRNA-mRNA pairs from bulk RNA-seq with snRNA-seq data has indicated that NORAD, NEAT1, and STXBP5-AS1 are downregulated across 4, 4, and 2 FC cell types, respectively. At the same time, MIRLET7BHG and MALAT1 are upregulated in excitatory neurons of FC and neurons of CP, respectively. Here, it is hypothesized that downregulation of NORAD, NEAT1, and STXBP5-AS1, and upregulation of MIRLET7BHG and MALAT1 might deregulate respectively 51, 6, and 37, and 31 and 19 mRNAs in cell types of FC and CP. Afterward, 13 therapeutic miRNAs are traced that might safeguard against deregulated lncRNA-mRNA pairs of NORAD, NEAT1, and MIRLET7BHG in FC. This study helps to explain the plausible mechanism of post-COVID neurological manifestation and also to devise therapeutics against it.

Exploring attractor bifurcations in Boolean networks.

BACKGROUND: Boolean networks (BNs) provide an effective modelling formalism for various complex biochemical phenomena. Their long term behaviour is represented by attractors-subsets of the state space towards which the BN eventually converges. These are then typically linked to different biological phenotypes. Depending on various logical parameters, the structure and quality of attractors can undergo a significant change, known as a bifurcation. We present a methodology for analysing bifurcations in asynchronous parametrised Boolean networks. RESULTS: In this paper, we propose a computational framework employing advanced symbolic graph algorithms that enable the analysis of large networks with hundreds of Boolean variables. To visualise the results of this analysis, we developed a novel interactive presentation technique based on decision trees, allowing us to quickly uncover parameters crucial to the changes in the attractor landscape. As a whole, the methodology is implemented in our tool AEON. We evaluate the method's applicability on a complex human cell signalling network describing the activity of type-1 interferons and related molecules interacting with SARS-COV-2 virion. In particular, the analysis focuses on explaining the potential suppressive role of the recently proposed drug molecule GRL0617 on replication of the virus. CONCLUSIONS: The proposed method creates a working analogy to the concept of bifurcation analysis widely used in kinetic modelling to reveal the impact of parameters on the system's stability. The important feature of our tool is its unique capability to work fast with large-scale networks with a relatively large extent of unknown information. The results obtained in the case study are in agreement with the recent biological findings.