Genome-wide screening of SARS-CoV-2 infection-related genes based on the blood leukocytes sequencing data set of patients with COVID-19

OmicPredict: a framework for omics data prediction using ANOVA-Firefly algorithm for feature selection

High-throughput technologies and machine learning (ML), when applied to a huge pool of medical data such as omics data, result in efficient analysis. Recent research aims to apply and develop ML models to predict a disease well in time using available omics datasets. The present work proposed a framework, 'OmicPredict', deploying a hybrid feature selection method and deep neural network (DNN) model to predict multiple diseases using omics data. The hybrid feature selection method is developed using the Analysis of Variance (ANOVA) technique and firefly algorithm. The OmicPredict framework is applied to three case studies, Alzheimer's disease, Breast cancer, and Coronavirus disease 2019 (COVID-19). In the case study of Alzheimer's disease, the framework predicts patients using GSE33000 and GSE44770 dataset. In the case study of Breast cancer, the framework predicts human epidermal growth factor receptor 2 (HER2) subtype status using Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) dataset. In the case study of COVID-19, the framework performs patients' classification using GSE157103 dataset. The experimental results show that DNN model achieved an Area Under Curve (AUC) score of 0.949 for the Alzheimer's (GSE33000 and GSE44770) dataset. Furthermore, it achieved an AUC score of 0.987 and 0.989 for breast cancer (METABRIC) and COVID-19 (GSE157103) datasets, respectively, outperforming Random Forest, Naïve Bayes models, and the existing research.

APOBEC3-related mutations in the spike protein-encoding region facilitate SARS-CoV-2 evolution

SARS-CoV-2 evolves gradually to cause COVID-19 epidemic. One of driving forces of SARS-CoV-2 evolution might be activation of apolipoprotein B mRNA editing catalytic subunit-like protein 3 (APOBEC3) by inflammatory factors. Here, we aimed to elucidate the effect of the APOBEC3-related viral mutations on the infectivity and immune evasion of SARS-CoV-2. The APOBEC3-related C > U mutations ranked as the second most common mutation types in the SARS-CoV-2 genome. mRNA expression of APOBEC3A (A3A), APOBEC3B (A3B), and APOBEC3G (A3G) in peripheral blood cells increased with disease severity. A3B, a critical member of the APOBEC3 family, was significantly upregulated in both severe and moderate COVID-19 patients and positively associated with neutrophil proportion and COVID-19 severity. We identified USP18 protein, a key molecule centralizing the protein-protein interaction network of key APOBEC3 proteins. Furthermore, mRNA expression of USP18 was significantly correlated to ACE2 and TMPRSS2 expression in the tissue of upper airways. Knockdown of USP18 mRNA significantly decreased A3B expression. Ectopic expression of A3B gene increased SARS-CoV-2 infectivity. C > U mutations at S371F, S373L, and S375F significantly conferred with the immune escape of SARS-CoV-2. Thus, APOBEC3, whose expression are upregulated by inflammatory factors, might promote SARS-CoV-2 evolution and spread via upregulating USP18 level and facilitating the immune escape. A3B and USP18 might be therapeutic targets for interfering with SARS-CoV-2 evolution.

Multi-omics analysis uncovered systemic lupus erythematosus and COVID-19 crosstalk

BACKGROUND

Studies have highlighted a possible crosstalk between the pathogeneses of COVID-19 and systemic lupus erythematosus (SLE); however, the interactive mechanisms remain unclear. We aimed to elucidate the impact of COVID-19 on SLE using clinical information and the underlying mechanisms of both diseases.

METHODS

RNA-seq datasets were used to identify shared hub gene signatures between COVID-19 and SLE, while genome-wide association study datasets were used to delineate the interaction mechanisms of the key signaling pathways. Finally, single-cell RNA-seq datasets were used to determine the primary target cells expressing the shared hub genes and key signaling pathways.

RESULTS

COVID-19 may affect patients with SLE through hematologic involvement and exacerbated inflammatory responses. We identified 14 shared hub genes between COVID-19 and SLE that were significantly associated with interferon (IFN)-I/II. We also screened and obtained four core transcription factors related to these hub genes, confirming the regulatory role of the IFN-I/II-mediated Janus kinase/signal transducers and activators of transcription (JAK-STAT) signaling pathway on these hub genes. Further, SLE and COVID-19 can interact via IFN-I/II and IFN-I/II receptors, promoting the levels of monokines, including interleukin (IL)-6/10, tumor necrosis factor-α, and IFN-γ, and elevating the incidence rate and risk of cytokine release syndrome. Therefore, in SLE and COVID-19, both hub genes and core TFs are enriched within monocytes/macrophages.

CONCLUSIONS

The interaction between SLE and COVID-19 promotes the activation of the IFN-I/II-triggered JAK-STAT signaling pathway in monocytes/macrophages. These findings provide a new direction and rationale for diagnosing and treating patients with SLE-COVID-19 comorbidity.

Metacell-based differential expression analysis identifies cell type specific temporal gene response programs in COVID-19 patient PBMCs

By profiling gene expression in individual cells, single-cell RNA-sequencing (scRNA-seq) can resolve cellular heterogeneity and cell-type gene expression dynamics. Its application to time-series samples can identify temporal gene programs active in different cell types, for example, immune cells' responses to viral infection. However, current scRNA-seq analysis has limitations. One is the low number of genes detected per cell. The second is insufficient replicates (often 1-2) due to high experimental cost. The third lies in the data analysis-treating individual cells as independent measurements leads to inflated statistics. To address these, we explore a new computational framework, specifically whether "metacells" constructed to maintain cellular heterogeneity within individual cell types (or clusters) can be used as "replicates" for increasing statistical rigor. Toward this, we applied SEACells to a time-series scRNA-seq dataset from peripheral blood mononuclear cells (PBMCs) after SARS-CoV-2 infection to construct metacells, and used them in maSigPro for quadratic regression to find significantly differentially expressed genes (DEGs) over time, followed by clustering expression velocity trends. We showed that such metacells retained greater expression variances and produced more biologically meaningful DEGs compared to either metacells generated randomly or from simple pseudobulk methods. More specifically, this approach correctly identified the known ISG15 interferon response program in almost all PBMC cell types and many DEGs enriched in the previously defined SARS-CoV-2 infection response pathway. It also uncovered additional and more cell type-specific temporal gene expression programs. Overall, our results demonstrate that the metacell-pseudoreplicate strategy could potentially overcome the limitation of 1-2 replicates.

Nucleic acid biomarkers of immune response and cell and tissue damage in children with COVID-19 and MIS-C

Differential host responses in coronavirus disease 2019 (COVID-19) and multisystem inflammatory syndrome in children (MIS-C) remain poorly characterized. Here, we use next-generation sequencing to longitudinally analyze blood samples from pediatric patients with COVID-19 or MIS-C across three hospitals. Profiling of plasma cell-free nucleic acids uncovers distinct signatures of cell injury and death between COVID-19 and MIS-C, with increased multiorgan involvement in MIS-C encompassing diverse cell types, including endothelial and neuronal cells, and an enrichment of pyroptosis-related genes. Whole-blood RNA profiling reveals upregulation of similar pro-inflammatory pathways in COVID-19 and MIS-C but also MIS-C-specific downregulation of T cell-associated pathways. Profiling of plasma cell-free RNA and whole-blood RNA in paired samples yields different but complementary signatures for each disease state. Our work provides a systems-level view of immune responses and tissue damage in COVID-19 and MIS-C and informs future development of new disease biomarkers.

Role of E3 ubiquitin ligases and deubiquitinating enzymes in SARS-CoV-2 infection

Ever since its emergence in 2019, COVID-19 has rapidly disseminated worldwide, engendering a pervasive pandemic that has profoundly impacted healthcare systems and the socio-economic milieu. A plethora of studies has been conducted targeting its pathogenic virus, SARS-CoV-2, to find ways to combat COVID-19. The ubiquitin-proteasome system (UPS) is widely recognized as a crucial mechanism that regulates human biological activities by maintaining protein homeostasis. Within the UPS, the ubiquitination and deubiquitination, two reversible modifications, of substrate proteins have been extensively studied and implicated in the pathogenesis of SARS-CoV-2. The regulation of E3 ubiquitin ligases and DUBs(Deubiquitinating enzymes), which are key enzymes involved in the two modification processes, determines the fate of substrate proteins. Proteins associated with the pathogenesis of SARS-CoV-2 may be retained, degraded, or even activated, thus affecting the ultimate outcome of the confrontation between SARS-CoV-2 and the host. In other words, the clash between SARS-CoV-2 and the host can be viewed as a battle for dominance over E3 ubiquitin ligases and DUBs, from the standpoint of ubiquitin modification regulation. This review primarily aims to clarify the mechanisms by which the virus utilizes host E3 ubiquitin ligases and DUBs, along with its own viral proteins that have similar enzyme activities, to facilitate invasion, replication, escape, and inflammation. We believe that gaining a better understanding of the role of E3 ubiquitin ligases and DUBs in COVID-19 can offer novel and valuable insights for developing antiviral therapies.

A Survey of COVID-19 Diagnosis Using Routine Blood Tests with the Aid of Artificial Intelligence Techniques

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), causing a disease called COVID-19, is a class of acute respiratory syndrome that has considerably affected the global economy and healthcare system. This virus is diagnosed using a traditional technique known as the Reverse Transcription Polymerase Chain Reaction (RT-PCR) test. However, RT-PCR customarily outputs a lot of false-negative and incorrect results. Current works indicate that COVID-19 can also be diagnosed using imaging resolutions, including CT scans, X-rays, and blood tests. Nevertheless, X-rays and CT scans cannot always be used for patient screening because of high costs, radiation doses, and an insufficient number of devices. Therefore, there is a requirement for a less expensive and faster diagnostic model to recognize the positive and negative cases of COVID-19. Blood tests are easily performed and cost less than RT-PCR and imaging tests. Since biochemical parameters in routine blood tests vary during the COVID-19 infection, they may supply physicians with exact information about the diagnosis of COVID-19. This study reviewed some newly emerging artificial intelligence (AI)-based methods to diagnose COVID-19 using routine blood tests. We gathered information about research resources and inspected 92 articles that were carefully chosen from a variety of publishers, such as IEEE, Springer, Elsevier, and MDPI. Then, these 92 studies are classified into two tables which contain articles that use machine Learning and deep Learning models to diagnose COVID-19 while using routine blood test datasets. In these studies, for diagnosing COVID-19, Random Forest and logistic regression are the most widely used machine learning methods and the most widely used performance metrics are accuracy, sensitivity, specificity, and AUC. Finally, we conclude by discussing and analyzing these studies which use machine learning and deep learning models and routine blood test datasets for COVID-19 detection. This survey can be the starting point for a novice-/beginner-level researcher to perform on COVID-19 classification.

Immune responses of different COVID-19 vaccination strategies by analyzing single-cell RNA sequencing data from multiple tissues using machine learning methods

Multiple types of COVID-19 vaccines have been shown to be highly effective in preventing SARS-CoV-2 infection and in reducing post-infection symptoms. Almost all of these vaccines induce systemic immune responses, but differences in immune responses induced by different vaccination regimens are evident. This study aimed to reveal the differences in immune gene expression levels of different target cells under different vaccine strategies after SARS-CoV-2 infection in hamsters. A machine learning based process was designed to analyze single-cell transcriptomic data of different cell types from the blood, lung, and nasal mucosa of hamsters infected with SARS-CoV-2, including B and T cells from the blood and nasal cavity, macrophages from the lung and nasal cavity, alveolar epithelial and lung endothelial cells. The cohort was divided into five groups: non-vaccinated (control), 2*adenovirus (two doses of adenovirus vaccine), 2*attenuated (two doses of attenuated virus vaccine), 2*mRNA (two doses of mRNA vaccine), and mRNA/attenuated (primed by mRNA vaccine, boosted by attenuated vaccine). All genes were ranked using five signature ranking methods (LASSO, LightGBM, Monte Carlo feature selection, mRMR, and permutation feature importance). Some key genes that contributed to the analysis of immune changes, such as RPS23, DDX5, PFN1 in immune cells, and IRF9 and MX1 in tissue cells, were screened. Afterward, the five feature sorting lists were fed into the feature incremental selection framework, which contained two classification algorithms (decision tree [DT] and random forest [RF]), to construct optimal classifiers and generate quantitative rules. Results showed that random forest classifiers could provide relative higher performance than decision tree classifiers, whereas the DT classifiers provided quantitative rules that indicated special gene expression levels under different vaccine strategies. These findings may help us to develop better protective vaccination programs and new vaccines.

Discovering Common Pathogenic Mechanisms of COVID-19 and Parkinson Disease: An Integrated Bioinformatics Analysis

Coronavirus disease 2019 (COVID-19) has emerged since December 2019 and was later characterized as a pandemic by WHO, imposing a major public health threat globally. Our study aimed to identify common signatures from different biological levels to enlighten the current unclear association between COVID-19 and Parkinson's disease (PD) as a number of possible links, and hypotheses were reported in the literature. We have analyzed transcriptome data from peripheral blood mononuclear cells (PBMCs) of both COVID-19 and PD patients, resulting in a total of 81 common differentially expressed genes (DEGs). The functional enrichment analysis of common DEGs are mostly involved in the complement system, type II interferon gamma (IFNG) signaling pathway, oxidative damage, microglia pathogen phagocytosis pathway, and GABAergic synapse. The protein-protein interaction network (PPIN) construction was carried out followed by hub detection, revealing 10 hub genes (MX1, IFI27, C1QC, C1QA, IFI6, NFIX, C1S, XAF1, IFI35, and ELANE). Some of the hub genes were associated with molecular mechanisms such as Lewy bodies-induced inflammation, microglia activation, and cytokine storm. We investigated regulatory elements of hub genes at transcription factor and miRNA levels. The major transcription factors regulating hub genes are SOX2, XAF1, RUNX1, MITF, and SPI1. We propose that these events may have important roles in the onset or progression of PD. To sum up, our analysis describes possible mechanisms linking COVID-19 and PD, elucidating some unknown clues in between.

A machine learning approach utilizing DNA methylation as an accurate classifier of COVID-19 disease severity

Since the onset of the COVID-19 pandemic, increasing cases with variable outcomes continue globally because of variants and despite vaccines and therapies. There is a need to identify at-risk individuals early that would benefit from timely medical interventions. DNA methylation provides an opportunity to identify an epigenetic signature of individuals at increased risk. We utilized machine learning to identify DNA methylation signatures of COVID-19 disease from data available through NCBI Gene Expression Omnibus. A training cohort of 460 individuals (164 COVID-19-infected and 296 non-infected) and an external validation dataset of 128 individuals (102 COVID-19-infected and 26 non-COVID-associated pneumonia) were reanalyzed. Data was processed using ChAMP and beta values were logit transformed. The JADBio AutoML platform was leveraged to identify a methylation signature associated with severe COVID-19 disease. We identified a random forest classification model from 4 unique methylation sites with the power to discern individuals with severe COVID-19 disease. The average area under the curve of receiver operator characteristic (AUC-ROC) of the model was 0.933 and the average area under the precision-recall curve (AUC-PRC) was 0.965. When applied to our external validation, this model produced an AUC-ROC of 0.898 and an AUC-PRC of 0.864. These results further our understanding of the utility of DNA methylation in COVID-19 disease pathology and serve as a platform to inform future COVID-19 related studies.

XAF1 Protects Host against Emerging RNA Viruses by Stabilizing IRF1-Dependent Antiviral Immunity

XIAP-associated factor 1 (XAF1) is an interferon (IFN)-stimulated gene (ISG) that enhances IFN-induced apoptosis. However, it is unexplored whether XAF1 is essential for the host fighting against invaded viruses. Here, we find that XAF1 is significantly upregulated in the host cells infected with emerging RNA viruses, including influenza, Zika virus (ZIKV), and SARS-CoV-2. IFN regulatory factor 1 (IRF1), a key transcription factor in immune cells, determines the induction of XAF1 during antiviral immunity. Ectopic expression of XAF1 protects host cells against various RNA viruses independent of apoptosis. Knockout of XAF1 attenuates host antiviral innate immunity in vitro and in vivo, which leads to more severe lung injuries and higher mortality in the influenza infection mouse model. XAF1 stabilizes IRF1 protein by antagonizing the CHIP-mediated degradation of IRF1, thus inducing more antiviral IRF1 target genes, including DDX58, DDX60, MX1, and OAS2. Our study has described a protective role of XAF1 in the host antiviral innate immunity against RNA viruses. We have also elucidated the molecular mechanism that IRF1 and XAF1 form a positive feedback loop to induce rapid and robust antiviral immunity. IMPORTANCE Rapid and robust induction of antiviral genes is essential for the host to clear the invaded viruses. In addition to the IRF3/7-IFN-I-STAT1 signaling axis, the XAF1-IRF1 positive feedback loop synergistically or independently drives the transcription of antiviral genes. Moreover, XAF1 is a sensitive and reliable gene that positively correlates with the viral infection, suggesting that XAF1 is a potential diagnostic marker for viral infectious diseases. In addition to the antitumor role, our study has shown that XAF1 is essential for antiviral immunity. XAF1 is not only a proapoptotic ISG, but it also stabilizes the master transcription factor IRF1 to induce antiviral genes. IRF1 directly binds to the IRF-Es of its target gene promoters and drives their transcriptions, which suggests a unique role of the XAF1-IRF1 loop in antiviral innate immunity, particularly in the host defect of IFN-I signaling such as invertebrates.

Identifying novel host-based diagnostic biomarker panels for COVID-19: a whole-blood/nasopharyngeal transcriptome meta-analysis

BACKGROUND

Regardless of improvements in controlling the COVID-19 pandemic, the lack of comprehensive insight into SARS-COV-2 pathogenesis is still a sophisticated challenge. In order to deal with this challenge, we utilized advanced bioinformatics and machine learning algorithms to reveal more characteristics of SARS-COV-2 pathogenesis and introduce novel host response-based diagnostic biomarker panels.

METHODS

In the present study, eight published RNA-Seq datasets related to whole-blood (WB) and nasopharyngeal (NP) swab samples of patients with COVID-19, other viral and non-viral acute respiratory illnesses (ARIs), and healthy controls (HCs) were integrated. To define COVID-19 meta-signatures, Gene Ontology and pathway enrichment analyses were applied to compare COVID-19 with other similar diseases. Additionally, CIBERSORTx was executed in WB samples to detect the immune cell landscape. Furthermore, the optimum WB- and NP-based diagnostic biomarkers were identified via all the combinations of 3 to 9 selected features and the 2-phases machine learning (ML) method which implemented k-fold cross validation and independent test set validation.

RESULTS

The host gene meta-signatures obtained for SARS-COV-2 infection were different in the WB and NP samples. The gene ontology and enrichment results of the WB dataset represented the enhancement in inflammatory host response, cell cycle, and interferon signature in COVID-19 patients. Furthermore, NP samples of COVID-19 in comparison with HC and non-viral ARIs showed the significant upregulation of genes associated with cytokine production and defense response to the virus. In contrast, these pathways in COVID-19 compared to other viral ARIs were strikingly attenuated. Notably, immune cell proportions of WB samples altered in COVID-19 versus HC. Moreover, the optimum WB- and NP-based diagnostic panels after two phases of ML-based validation included 6 and 8 markers with an accuracy of 97% and 88%, respectively.

CONCLUSIONS

Based on the distinct gene expression profiles of WB and NP, our results indicated that SARS-COV-2 function is body-site-specific, although according to the common signature in WB and NP COVID-19 samples versus controls, this virus also induces a global and systematic host response to some extent. We also introduced and validated WB- and NP-based diagnostic biomarkers using ML methods which can be applied as a complementary tool to diagnose the COVID-19 infection from non-COVID cases.

Identification of COVID-19-Specific Immune Markers Using a Machine Learning Method

Notably, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a tight relationship with the immune system. Human resistance to COVID-19 infection comprises two stages. The first stage is immune defense, while the second stage is extensive inflammation. This process is further divided into innate and adaptive immunity during the immune defense phase. These two stages involve various immune cells, including CD4⁺ T cells, CD8⁺ T cells, monocytes, dendritic cells, B cells, and natural killer cells. Various immune cells are involved and make up the complex and unique immune system response to COVID-19, providing characteristics that set it apart from other respiratory infectious diseases. In the present study, we identified cell markers for differentiating COVID-19 from common inflammatory responses, non-COVID-19 severe respiratory diseases, and healthy populations based on single-cell profiling of the gene expression of six immune cell types by using Boruta and mRMR feature selection methods. Some features such as IFI44L in B cells, S100A8 in monocytes, and NCR2 in natural killer cells are involved in the innate immune response of COVID-19. Other features such as ZFP36L2 in CD4⁺ T cells can regulate the inflammatory process of COVID-19. Subsequently, the IFS method was used to determine the best feature subsets and classifiers in the six immune cell types for two classification algorithms. Furthermore, we established the quantitative rules used to distinguish the disease status. The results of this study can provide theoretical support for a more in-depth investigation of COVID-19 pathogenesis and intervention strategies.

SARS-CoV-2 potential drugs, drug targets, and biomarkers: a viral-host interaction network-based analysis

COVID-19 is a global pandemic impacting the daily living of millions. As variants of the virus evolve, a complete comprehension of the disease and drug targets becomes a decisive duty. The Omicron variant, for example, has a notably high transmission rate verified in 155 countries. We performed integrative transcriptomic and network analyses to identify drug targets and diagnostic biomarkers and repurpose FDA-approved drugs for SARS-CoV-2. Upon the enrichment of 464 differentially expressed genes, pathways regulating the host cell cycle were significant. Regulatory and interaction networks featured hsa-mir-93-5p and hsa-mir-17-5p as blood biomarkers while hsa-mir-15b-5p as an antiviral agent. MYB, RRM2, ERG, CENPF, CIT, and TOP2A are potential drug targets for treatment. HMOX1 is suggested as a prognostic biomarker. Enhancing HMOX1 expression by neem plant extract might be a therapeutic alternative. We constructed a drug-gene network for FDA-approved drugs to be repurposed against the infection. The key drugs retrieved were members of anthracyclines, mitotic inhibitors, anti-tumor antibiotics, and CDK1 inhibitors. Additionally, hydroxyquinone and digitoxin are potent TOP2A inhibitors. Hydroxyurea, cytarabine, gemcitabine, sotalol, and amiodarone can also be redirected against COVID-19. The analysis enforced the repositioning of fluorouracil and doxorubicin, especially that they have multiple drug targets, hence less probability of resistance.

XGBoost-Based Feature Learning Method for Mining COVID-19 Novel Diagnostic Markers

In December 2019, an outbreak of novel coronavirus pneumonia spread over Wuhan, Hubei Province, China, which then developed into a significant global health public event, giving rise to substantial economic losses. We downloaded throat swab expression profiling data of COVID-19 positive and negative patients from the Gene Expression Omnibus (GEO) database to mine novel diagnostic biomarkers. XGBoost was used to construct the model and select feature genes. Subsequently, we constructed COVID-19 classifiers such as MARS, KNN, SVM, MIL, and RF using machine learning methods. We selected the KNN classifier with the optimal MCC value from these classifiers using the IFS method to identify 24 feature genes. Finally, we used principal component analysis to classify the samples and found that the 24 feature genes could effectively be used to classify COVID-19-positive and negative patients. Additionally, we analyzed the possible biological functions and signaling pathways in which the 24 feature genes were involved by GO and KEGG enrichment analyses. The results demonstrated that these feature genes were primarily enriched in biological functions such as viral transcription and viral gene expression and pathways such as Coronavirus disease-COVID-19. In summary, the 24 feature genes we identified were highly effective in classifying COVID-19 positive and negative patients, which could serve as novel markers for COVID-19.

Supporting Clinical COVID-19 Diagnosis with Routine Blood Tests Using Tree-Based Entropy Structured Self-Organizing Maps

Data classification is an automatic or semi-automatic process that, utilizing artificial intelligence algorithms, learns the variable and class relationships of a dataset for use a posteriori in situations where the class result is unknown. For many years, work on this topic has been aimed at increasing the hit rates of algorithms. However, when the problem is restricted to applications in healthcare, besides the concern with performance, it is also necessary to design algorithms whose results are understandable by the specialists responsible for making the decisions. Among the problems in the field of medicine, a current focus is related to COVID-19: AI algorithms may contribute to early diagnosis. Among the available COVID-19 data, the blood test is a typical procedure performed when the patient seeks the hospital, and its use in the diagnosis allows reducing the need for other diagnostic tests that can impact the detection time and add to costs. In this work, we propose using self-organizing map (SOM) to discover attributes in blood test examinations that are relevant for COVID-19 diagnosis. We applied SOM and an entropy calculation in the definition of a hierarchical, semi-supervised and explainable model named TESSOM (tree-based entropy-structured self-organizing maps), in which the main feature is enhancing the investigation of groups of cases with high levels of class overlap, as far as the diagnostic outcome is concerned. Framing the TESSOM algorithm in the context of explainable artificial intelligence (XAI) makes it possible to explain the results to an expert in a simplified way. It is demonstrated in the paper that the use of the TESSOM algorithm to identify attributes of blood tests can help with the identification of COVID-19 cases. It providing a performance increase in 1.489% in multiple scenarios when analyzing 2207 cases from three hospitals in the state of São Paulo, Brazil. This work is a starting point for researchers to identify relevant attributes of blood tests for COVID-19 and to support the diagnosis of other diseases.

Functional Roles and Targets of COVID-19 in Blood Cells Determined Using Bioinformatics Analyses

Background

Coronavirus disease 2019 (COVID-19), caused by novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a global epidemic with a high mortality rate. In this study, our goal was to identify the function and associated targets of SARS-CoV-2 from circulating monocytes in the blood and peripheral blood mononuclear cell (PBMC) dataset of patients with COVID-19.

Methods

The Gene Expression Omnibus database (GSE164805 and GSE180594) was used to identify differentially expressed genes (DEGs). Gene ontology function analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses of the DEGs were performed using the DAVID database.

Results

Gene ontology analysis of DEG revealed that GSE164805 and GSE180594 were involved in the regulation of cell migration, upregulation of cell proliferation, and in the activation of the mitogen-activated protein kinase signaling pathway. Kyoto Encyclopedia of Genes and Genomes analysis of GSE164805 revealed that the DEGs were enriched in peroxisome, melanogenesis, and actin regulation. Peroxisome genes were highly expressed in patients with mild and severe COVID-19. Bioinformatics analysis to compare GSE180594 and public data for the single-cell atlas of the peripheral immune response in patients with COVID-19 showed that interferon-associated genes were highly increased in acute COVID-19 PBMC and in CD14+ and CD16+ monocytes from patients with COVID-19.

Conclusions

We comprehensively analyzed the blood cell gene expression profile data of patients with COVID-19 using bioinformatics methods to preliminary understand the functions and associated targets of DEGs in the blood cells of these patients. Thus, our data provide targets for potential therapies against COVID-19.

Differential Co-Expression Network Analysis Reveals Key Hub-High Traffic Genes as Potential Therapeutic Targets for COVID-19 Pandemic

Background

The recent emergence of COVID-19, rapid worldwide spread, and incomplete knowledge of molecular mechanisms underlying SARS-CoV-2 infection have limited development of therapeutic strategies. Our objective was to systematically investigate molecular regulatory mechanisms of COVID-19, using a combination of high throughput RNA-sequencing-based transcriptomics and systems biology approaches.

Methods

RNA-Seq data from peripheral blood mononuclear cells (PBMCs) of healthy persons, mild and severe 17 COVID-19 patients were analyzed to generate a gene expression matrix. Weighted gene co-expression network analysis (WGCNA) was used to identify co-expression modules in healthy samples as a reference set. For differential co-expression network analysis, module preservation and module-trait relationships approaches were used to identify key modules. Then, protein-protein interaction (PPI) networks, based on co-expressed hub genes, were constructed to identify hub genes/TFs with the highest information transfer (hub-high traffic genes) within candidate modules.

Results

Based on differential co-expression network analysis, connectivity patterns and network density, 72% (15 of 21) of modules identified in healthy samples were altered by SARS-CoV-2 infection. Therefore, SARS-CoV-2 caused systemic perturbations in host biological gene networks. In functional enrichment analysis, among 15 non-preserved modules and two significant highly-correlated modules (identified by MTRs), 9 modules were directly related to the host immune response and COVID-19 immunopathogenesis. Intriguingly, systemic investigation of SARS-CoV-2 infection identified signaling pathways and key genes/proteins associated with COVID-19's main hallmarks, e.g., cytokine storm, respiratory distress syndrome (ARDS), acute lung injury (ALI), lymphopenia, coagulation disorders, thrombosis, and pregnancy complications, as well as comorbidities associated with COVID-19, e.g., asthma, diabetic complications, cardiovascular diseases (CVDs), liver disorders and acute kidney injury (AKI). Topological analysis with betweenness centrality (BC) identified 290 hub-high traffic genes, central in both co-expression and PPI networks. We also identified several transcriptional regulatory factors, including NFKB1, HIF1A, AHR, and TP53, with important immunoregulatory roles in SARS-CoV-2 infection. Moreover, several hub-high traffic genes, including IL6, IL1B, IL10, TNF, SOCS1, SOCS3, ICAM1, PTEN, RHOA, GDI2, SUMO1, CASP1, IRAK3, HSPA5, ADRB2, PRF1, GZMB, OASL, CCL5, HSP90AA1, HSPD1, IFNG, MAPK1, RAB5A, and TNFRSF1A had the highest rates of information transfer in 9 candidate modules and central roles in COVID-19 immunopathogenesis.

Conclusion

This study provides comprehensive information on molecular mechanisms of SARS-CoV-2-host interactions and identifies several hub-high traffic genes as promising therapeutic targets for the COVID-19 pandemic.

Modifying effects of TNF-α, IL-6 and VDR genes on the development risk and the course of COVID-19. Pilot study

OBJECTIVES

COVID-19 continues to range around the world and set morbidity and mortality antirecords. Determining the role of genetic factors in the development of COVID-19 may contribute to the understanding of the pathogenetic mechanisms that lead to the development of complications and fatalities in this disease. The aim of our study was to analyze the effect of TNF-α (rs1800629), IL-6 (rs1800795) and VDR (rs731236 and rs1544410) genes variants on the development risk and the course of COVID-19 in intensive care patients.

METHODS

The study group included 31 patients with diagnosis "viral COVID-19 pneumonia". All patients underwent standard daily repeated clinical, instrumental and laboratory examinations. Determination of IL-6, TNF-α, and VDR genes variants was performed using the PCR-RFLP method.

RESULTS

It was found a significant increase in the rate of the CC genotype and C allele (38.7 vs. 12.0% and 0.6 vs. 0.4%, respectively) of the IL-6 gene in all patients of the study in comparison with population frequencies. There was a significantly higher rate of heterozygous genotypes TC and GA of the VDR gene in group of died patients. The rs1800629 variant of the TNF-α gene is associated with the need for respiratory support and its longer duration in patients with COVID-19.

CONCLUSIONS

The obtained results support a hypothesis about the influence of variants of IL-6, TNF-α and VDR genes on severity of COVID-19. However, in order to draw definite conclusions, further multifaceted research in this area are need.