Identification of host transcriptome-guided repurposable drugs for SARS-CoV-1 infections and their validation with SARS-CoV-2 infections by using the integrated bioinformatics approaches

Discovery of common molecular signatures and drug repurposing for COVID-19/Asthma comorbidity: ACE2 and multi-partite networks

Angiotensin-converting enzyme 2 (ACE2) is identified as the functional receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of the ongoing global coronavirus disease-2019 (COVID-19) pandemic. This study aimed to elucidate potential therapeutic avenues by scrutinizing approved drugs through the identification of the genetic signature associated with SARS-CoV-2 infection in individuals with asthma. This exploration was conducted through an integrated analysis, encompassing interaction networks between the ACE2 receptor and common host (co-host) factors implicated in COVID-19/asthma comorbidity. The comprehensive analysis involved the identification of common differentially expressed genes (cDEGs) and hub-cDEGs, functional annotations, interaction networks, gene set variation analysis (GSVA), gene set enrichment analysis (GSEA), and module construction. Interaction networks were used to identify overlapping disease modules and potential drug targets. Computational biology and molecular docking analyzes were utilized to discern functional drug modules. Subsequently, the impact of the identified drugs on the expression of hub-cDEGs was experimentally validated using a mouse model. A total of 153 cDEGs or co-host factors associated with ACE2 were identified in the COVID-19 and asthma comorbidity. Among these, seven significant cDEGs and proteins - namely, HRAS, IFNG, JUN, CDH1, TLR4, ICAM1, and SCD-were recognized as pivotal host factors linked to ACE2. Regulatory network analysis of hub-cDEGs revealed eight top-ranked transcription factors (TFs) proteins and nine microRNAs as key regulatory factors operating at the transcriptional and post-transcriptional levels, respectively. Molecular docking simulations led to the proposal of 10 top-ranked repurposable drug molecules (Rapamycin, Ivermectin, Everolimus, Quercetin, Estradiol, Entrectinib, Nilotinib, Conivaptan, Radotinib, and Venetoclax) as potential treatment options for COVID-19 in individuals with comorbid asthma. Validation analysis demonstrated that Rapamycin effectively inhibited ICAM1 expression in the HDM-stimulated mice group (p < 0.01). This study unveils the common pathogenesis and genetic signature underlying asthma and SARS-CoV-2 infection, delineated by the interaction networks of ACE2-related host factors. These findings provide valuable insights for the design and discovery of drugs aimed at more effective therapeutics within the context of lung disease comorbidities.

Discovery of Bacterial Key Genes from 16S rRNA-Seq Profiles That Are Associated with the Complications of SARS-CoV-2 Infections and Provide Therapeutic Indications

SARS-CoV-2 infections, commonly referred to as COVID-19, remain a critical risk to both human life and global economies. Particularly, COVID-19 patients with weak immunity may suffer from different complications due to the bacterial co-infections/super-infections/secondary infections. Therefore, different variants of alternative antibacterial therapeutic agents are required to inhibit those infection-causing drug-resistant pathogenic bacteria. This study attempted to explore these bacterial pathogens and their inhibitors by using integrated statistical and bioinformatics approaches. By analyzing bacterial 16S rRNA sequence profiles, at first, we detected five bacterial genera and taxa (Bacteroides, Parabacteroides, Prevotella Clostridium, Atopobium, and Peptostreptococcus) based on differentially abundant bacteria between SARS-CoV-2 infection and control samples that are significantly enriched in 23 metabolic pathways. A total of 183 bacterial genes were found in the enriched pathways. Then, the top-ranked 10 bacterial genes (accB, ftsB, glyQ, hldD, lpxC, lptD, mlaA, ppsA, ppc, and tamB) were selected as the pathogenic bacterial key genes (bKGs) by their protein-protein interaction (PPI) network analysis. Then, we detected bKG-guided top-ranked eight drug molecules (Bemcentinib, Ledipasvir, Velpatasvir, Tirilazad, Acetyldigitoxin, Entreatinib, Digitoxin, and Elbasvir) by molecular docking. Finally, the binding stability of the top-ranked three drug molecules (Bemcentinib, Ledipasvir, and Velpatasvir) against three receptors (hldD, mlaA, and lptD) was investigated by computing their binding free energies with molecular dynamic (MD) simulation-based MM-PBSA techniques, respectively, and was found to be stable. Therefore, the findings of this study could be useful resources for developing a proper treatment plan against bacterial co-/super-/secondary-infection in SARS-CoV-2 infections.

Identification of shared pathogenetic mechanisms between COVID-19 and IC through bioinformatics and system biology

COVID-19 increased global mortality in 2019. Cystitis became a contributing factor in SARS-CoV-2 and COVID-19 complications. The complex molecular links between cystitis and COVID-19 are unclear. This study investigates COVID-19-associated cystitis (CAC) molecular mechanisms and drug candidates using bioinformatics and systems biology. Obtain the gene expression profiles of IC (GSE11783) and COVID-19 (GSE147507) from the Gene Expression Omnibus (GEO) database. Identified the common differentially expressed genes (DEGs) in both IC and COVID-19, and extracted a number of key genes from this group. Subsequently, conduct Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis on the DEGs. Additionally, design a protein-protein interaction (PPI) network, a transcription factor gene regulatory network, a TF miRNA regulatory network, and a gene disease association network using the DEGs. Identify and extract hub genes from the PPI network. Then construct Nomogram diagnostic prediction models based on the hub genes. The DSigDB database was used to forecast many potential molecular medicines that are associated with common DEGs. Assess the precision of hub genes and Nomogram models in diagnosing IC and COVID-19 by employing Receiver Operating Characteristic (ROC) curves. The IC dataset (GSE57560) and the COVID-19 dataset (GSE171110) were selected to validate the models' diagnostic accuracy. A grand total of 198 DEGs that overlapped were found and chosen for further research. FCER1G, ITGAM, LCP2, LILRB2, MNDA, SPI1, and TYROBP were screened as the hub genes. The Nomogram model, built using the seven hub genes, demonstrates significant utility as a diagnostic prediction model for both IC and COVID-19. Multiple potential molecular medicines associated with common DEGs have been discovered. These pathways, hub genes, and models may provide new perspectives for future research into mechanisms and guide personalised and effective therapeutics for IC patients infected with COVID-19.

Biological analysis of the potential pathogenic mechanisms of Infectious COVID-19 and Guillain-Barré syndrome

Background

Guillain-Barré syndrome (GBS) is a medical condition characterized by the immune system of the body attacking the peripheral nerves, including those in the spinal nerve roots, peripheral nerves, and cranial nerves. It can cause limb weakness, abnormal sensations, and facial nerve paralysis. Some studies have reported clinical cases associated with the severe coronavirus disease 2019 (COVID-19) and GBS, but how COVID-19 affects GBS is unclear.

Methods

We utilized bioinformatics techniques to explore the potential genetic connection between COVID-19 and GBS. Differential expression of genes (DEGs) related to COVID-19 and GBS was collected from the Gene Expression Omnibus (GEO) database. By taking the intersection, we obtained shared DEGs for COVID-19 and GBS. Subsequently, we utilized bioinformatics analysis tools to analyze common DEGs, conducting functional enrichment analysis and constructing Protein-protein interaction networks (PPI), Transcription factors (TF) -gene networks, and TF-miRNA networks. Finally, we validated our findings by constructing the Receiver Operating Characteristic (ROC) curves.

Results

This study utilizes bioinformatics tools for the first time to investigate the close genetic relationship between COVID-19 and GBS. CAMP, LTF, DEFA1B, SAMD9, GBP1, DDX60, DEFA4, and OAS3 are identified as the most significant interacting genes between COVID-19 and GBS. In addition, the signaling pathway of NOD-like receptors is believed to be essential in the link between COVID-19 and GBS.

Genetic risk factors for severe and fatigue dominant long COVID and commonalities with ME/CFS identified by combinatorial analysis

BACKGROUND

Long COVID is a debilitating chronic condition that has affected over 100 million people globally. It is characterized by a diverse array of symptoms, including fatigue, cognitive dysfunction and respiratory problems. Studies have so far largely failed to identify genetic associations, the mechanisms behind the disease, or any common pathophysiology with other conditions such as myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) that present with similar symptoms.

METHODS

We used a combinatorial analysis approach to identify combinations of genetic variants significantly associated with the development of long COVID and to examine the biological mechanisms underpinning its various symptoms. We compared two subpopulations of long COVID patients from Sano Genetics' Long COVID GOLD study cohort, focusing on patients with severe or fatigue dominant phenotypes. We evaluated the genetic signatures previously identified in an ME/CFS population against this long COVID population to understand similarities with other fatigue disorders that may be triggered by a prior viral infection. Finally, we also compared the output of this long COVID analysis against known genetic associations in other chronic diseases, including a range of metabolic and neurological disorders, to understand the overlap of pathophysiological mechanisms.

RESULTS

Combinatorial analysis identified 73 genes that were highly associated with at least one of the long COVID populations included in this analysis. Of these, 9 genes have prior associations with acute COVID-19, and 14 were differentially expressed in a transcriptomic analysis of long COVID patients. A pathway enrichment analysis revealed that the biological pathways most significantly associated with the 73 long COVID genes were mainly aligned with neurological and cardiometabolic diseases. Expanded genotype analysis suggests that specific SNX9 genotypes are a significant contributor to the risk of or protection against severe long COVID infection, but that the gene-disease relationship is context dependent and mediated by interactions with KLF15 and RYR3. Comparison of the genes uniquely associated with the Severe and Fatigue Dominant long COVID patients revealed significant differences between the pathways enriched in each subgroup. The genes unique to Severe long COVID patients were associated with immune pathways such as myeloid differentiation and macrophage foam cells. Genes unique to the Fatigue Dominant subgroup were enriched in metabolic pathways such as MAPK/JNK signaling. We also identified overlap in the genes associated with Fatigue Dominant long COVID and ME/CFS, including several involved in circadian rhythm regulation and insulin regulation. Overall, 39 SNPs associated in this study with long COVID can be linked to 9 genes identified in a recent combinatorial analysis of ME/CFS patient from UK Biobank. Among the 73 genes associated with long COVID, 42 are potentially tractable for novel drug discovery approaches, with 13 of these already targeted by drugs in clinical development pipelines. From this analysis for example, we identified TLR4 antagonists as repurposing candidates with potential to protect against long term cognitive impairment pathology caused by SARS-CoV-2. We are currently evaluating the repurposing potential of these drug targets for use in treating long COVID and/or ME/CFS.

CONCLUSION

This study demonstrates the power of combinatorial analytics for stratifying heterogeneous populations in complex diseases that do not have simple monogenic etiologies. These results build upon the genetic findings from combinatorial analyses of severe acute COVID-19 patients and an ME/CFS population and we expect that access to additional independent, larger patient datasets will further improve the disease insights and validate potential treatment options in long COVID.

Pharmacogenomic Landscape of Ivermectin and Selective Antioxidants: Exploring Gene Interplay in the Context of Long COVID

COVID-19 pandemic has caused widespread panic and fear among the global population. As such, repurposing drugs are being used as viable therapeutic options due to the limited effective treatments for Long COVID symptoms. Ivermectin is one of the emerging repurposed drugs that has been shown effective to have antiviral effects in clinical trials. In addition, antioxidant compounds are also gaining attention due to their capabilities of reducing inflammation and severity of symptoms. Due to the absence of knowledge in pharmacogenomics and modes of actions in the human body for these compounds, this study aims to provide a pharmacogenomic profile for the combination of ivermectin and six selected antioxidants (epigallocatechin gallate (EGCG), curcumin, sesamin, anthocyanins, quercetin, and N-acetylcysteine (NAC)) as potentially effective regimens for long COVID symptoms. Results showed that there were 12 interacting genes found among the ivermectin, 6 antioxidants, and COVID-19. For network pharmacology, the 12 common interacting genes/proteins had the highest associations with Pertussis pathway, AGE-RAGE signaling pathway in diabetic complications, and colorectal cancer in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Disease analyses also revealed that the top three relevant diseases with COVID-19 infections were diabetes mellitus, ischemia, reperfusion injury. We also identified 6 potential target microRNAs (miRNAs) of the 12 commonly curated genes used as molecular biomarkers for COVID-19 treatments. The established pharmacogenomic network, disease analyses, and identified miRNAs could facilitate developments of effective regimens for chronic sequelae of COVID-19 especially in this post-pandemic era. However, further studies and clinical trials are needed to substantiate the effectiveness and dosages for COVID-19 treatments.

Bioinformatics and system biology approach to identify the influences of SARS-CoV-2 on metabolic unhealthy obese patients

Introduction: The severe acute respiratory syndrome coronavirus 2 (SARS-COV-2) has posed a significant challenge to individuals' health. Increasing evidence shows that patients with metabolic unhealthy obesity (MUO) and COVID-19 have severer complications and higher mortality rate. However, the molecular mechanisms underlying the association between MUO and COVID-19 are poorly understood. Methods: We sought to reveal the relationship between MUO and COVID-19 using bioinformatics and systems biology analysis approaches. Here, two datasets (GSE196822 and GSE152991) were employed to extract differentially expressed genes (DEGs) to identify common hub genes, shared pathways, transcriptional regulatory networks, gene-disease relationship and candidate drugs. Results: Based on the identified 65 common DEGs, the complement-related pathways and neutrophil degranulation-related functions are found to be mainly affected. The hub genes, which included SPI1, CD163, C1QB, SIGLEC1, C1QA, ITGAM, CD14, FCGR1A, VSIG4 and C1QC, were identified. From the interaction network analysis, 65 transcription factors (TFs) were found to be the regulatory signals. Some infections, inflammation and liver diseases were found to be most coordinated with the hub genes. Importantly, Paricalcitol, 3,3',4,4',5-Pentachlorobiphenyl, PD 98059, Medroxyprogesterone acetate, Dexamethasone and Tretinoin HL60 UP have shown possibility as therapeutic agents against COVID-19 and MUO. Conclusion: This study provides new clues and references to treat both COVID-19 and MUO.

In-silico study for the identification of potential destabilizers between the spike protein of SARS-CoV-2 and human ACE-2

The emergence of the new SARS-CoV-2 virus, which causes the disease known as COVID-19, has generated a pandemic that has plunged the world into a health crisis. The infection process is triggered by the direct binding of the receptor-binding domain (RBD) of the spike (S) protein of SARS-CoV-2 to the angiotensin-converting enzyme 2 (ACE2) of the host cell. In the present study, virtual screening techniques such as molecular docking, molecular dynamics, calculation of free energy using the GBSA method, prediction of drug similarity, pharmacokinetic, and toxicological properties of various ligands interacting with the RBD-ACE2 complex were applied. The ligands radotinib, hinokiflavone, and ginkgetin were identified as potential destabilizers of the RBD-ACE2 interaction, which could produce their pharmacological effect by interacting at an allosteric site of ACE2, with affinity energy values of -10.2 ± 0.1, -9.8 ± 0.0, and -9.4 ± 0.0 kcal/mol, indicating strong receptor affinity. The complex with hinokiflavone showed the highest conformational stability and rigidity of the dynamic simulation and also obtained the best binding free energy of the three molecules, with an energy of -215.86 kcal/mol.

Bioinformatics and system biology approach to identify the influences among COVID-19, ARDS and sepsis

Background Severe coronavirus disease 2019 (COVID -19) has led to severe pneumonia or acute respiratory distress syndrome (ARDS) worldwide. we have noted that many critically ill patients with COVID-19 present with typical sepsis-related clinical manifestations, including multiple organ dysfunction syndrome, coagulopathy, and septic shock. The molecular mechanisms that underlie COVID-19, ARDS and sepsis are not well understood. The objectives of this study were to analyze potential molecular mechanisms and identify potential drugs for the treatment of COVID-19, ARDS and sepsis using bioinformatics and a systems biology approach. Methods Three RNA-seq datasets (GSE171110, GSE76293 and GSE137342) from Gene Expression Omnibus (GEO) were employed to detect mutual differentially expressed genes (DEGs) for the patients with the COVID-19, ARDS and sepsis for functional enrichment, pathway analysis, and candidate drugs analysis. Results We obtained 110 common DEGs among COVID-19, ARDS and sepsis. ARG1, FCGR1A, MPO, and TLR5 are the most influential hub genes. The infection and immune-related pathways and functions are the main pathways and molecular functions of these three diseases. FOXC1, YY1, GATA2, FOXL, STAT1 and STAT3 are important TFs for COVID-19. mir-335-5p, miR-335-5p and hsa-mir-26a-5p were associated with COVID-19. Finally, the hub genes retrieved from the DSigDB database indicate multiple drug molecules and drug-targets interaction. Conclusion We performed a functional analysis under ontology terms and pathway analysis and found some common associations among COVID-19, ARDS and sepsis. Transcription factors-genes interaction, protein-drug interactions, and DEGs-miRNAs coregulatory network with common DEGs were also identified on the datasets. We believe that the candidate drugs obtained in this study may contribute to the effective treatment of COVID-19.

Identification of host genomic biomarkers from multiple transcriptomics datasets for diagnosis and therapies of SARS-CoV-2 infections

The pandemic of COVID-19 is a severe threat to human life and the global economy. Despite the success of vaccination efforts in reducing the spread of the virus, the situation remains largely uncontrolled due to the random mutation in the RNA sequence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which demands different variants of effective drugs. Disease-causing gene-mediated proteins are usually used as receptors to explore effective drug molecules. In this study, we analyzed two different RNA-Seq and one microarray gene expression profile datasets by integrating EdgeR, LIMMA, weighted gene co-expression network and robust rank aggregation approaches, which revealed SARS-CoV-2 infection causing eight hub-genes (HubGs) including HubGs; REL, AURKA, AURKB, FBXL3, OAS1, STAT4, MMP2 and IL6 as the host genomic biomarkers. Gene Ontology and pathway enrichment analyses of HubGs significantly enriched some crucial biological processes, molecular functions, cellular components and signaling pathways that are associated with the mechanisms of SARS-CoV-2 infections. Regulatory network analysis identified top-ranked 5 TFs (SRF, PBX1, MEIS1, ESR1 and MYC) and 5 miRNAs (hsa-miR-106b-5p, hsa-miR-20b-5p, hsa-miR-93-5p, hsa-miR-106a-5p and hsa-miR-20a-5p) as the key transcriptional and post-transcriptional regulators of HubGs. Then, we conducted a molecular docking analysis to determine potential drug candidates that could interact with HubGs-mediated receptors. This analysis resulted in the identification of top-ranked ten drug agents, including Nilotinib, Tegobuvir, Digoxin, Proscillaridin, Olysio, Simeprevir, Hesperidin, Oleanolic Acid, Naltrindole and Danoprevir. Finally, we investigated the binding stability of the top-ranked three drug molecules Nilotinib, Tegobuvir and Proscillaridin with the three top-ranked proposed receptors (AURKA, AURKB, OAS1) by using 100 ns MD-based MM-PBSA simulations and observed their stable performance. Therefore, the findings of this study might be useful resources for diagnosis and therapies of SARS-CoV-2 infections.

MicroRNA-155 is a main part of proinflammatory puzzle during severe coronavirus disease 2019 (COVID-19)

Genetic and epigenetic parameters play critical roles in determining the outcomes of the severe acute respiratory syndrome coronavirus type 19 (SARS-CoV-2) infection. MicroRNAs (miRNAs) are an important part of the epigenetic factors that regulate several functions of the immune cells and also viruses. Accordingly, the molecules can regulate expression of the immune cell proteins and virus in the host cells. Among the miRNAs, miRNA-155 (miR-155) is well-studied in patients suffering from severe coronavirus disease 2019 (COVID-19). It has been reported that the SARS-CoV-2 infected patients may be directed to induce a cytokine storm or severe proinflammatory responses. This review article discusses the pathological roles of miR-155 during COVID-19 infection.

Discovering common pathogenetic processes between COVID-19 and sepsis by bioinformatics and system biology approach

Corona Virus Disease 2019 (COVID-19), an acute respiratory infectious disease caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has spread rapidly worldwide, resulting in a pandemic with a high mortality rate. In clinical practice, we have noted that many critically ill or critically ill patients with COVID-19 present with typical sepsis-related clinical manifestations, including multiple organ dysfunction syndrome, coagulopathy, and septic shock. In addition, it has been demonstrated that severe COVID-19 has some pathological similarities with sepsis, such as cytokine storm, hypercoagulable state after blood balance is disrupted and neutrophil dysfunction. Considering the parallels between COVID-19 and non-SARS-CoV-2 induced sepsis (hereafter referred to as sepsis), the aim of this study was to analyze the underlying molecular mechanisms between these two diseases by bioinformatics and a systems biology approach, providing new insights into the pathogenesis of COVID-19 and the development of new treatments. Specifically, the gene expression profiles of COVID-19 and sepsis patients were obtained from the Gene Expression Omnibus (GEO) database and compared to extract common differentially expressed genes (DEGs). Subsequently, common DEGs were used to investigate the genetic links between COVID-19 and sepsis. Based on enrichment analysis of common DEGs, many pathways closely related to inflammatory response were observed, such as Cytokine-cytokine receptor interaction pathway and NF-kappa B signaling pathway. In addition, protein-protein interaction networks and gene regulatory networks of common DEGs were constructed, and the analysis results showed that ITGAM may be a potential key biomarker base on regulatory analysis. Furthermore, a disease diagnostic model and risk prediction nomogram for COVID-19 were constructed using machine learning methods. Finally, potential therapeutic agents, including progesterone and emetine, were screened through drug-protein interaction networks and molecular docking simulations. We hope to provide new strategies for future research and treatment related to COVID-19 by elucidating the pathogenesis and genetic mechanisms between COVID-19 and sepsis.

Metadata analysis to explore hub of the hub-genes highlighting their functions, pathways and regulators for cervical cancer diagnosis and therapies

Cervical cancer (CC) is considered as the fourth most common women cancer globally.that shows malignant features of local infiltration and invasion into adjacent organs and tissues. There are several individual studies in the literature that explored CC-causing hub-genes (HubGs), however, we observed that their results are not so consistent. Therefore, the main objective of this study was to explore hub of the HubGs (hHubGs) that might be more representative CC-causing HubGs compare to the single study based HubGs. We reviewed 52 published articles and found 255 HubGs/studied-genes in total. Among them, we selected 10 HubGs (CDK1, CDK2, CHEK1, MKI67, TOP2A, BRCA1, PLK1, CCNA2, CCNB1, TYMS) as the hHubGs by the protein-protein interaction (PPI) network analysis. Then, we validated their differential expression patterns between CC and control samples through the GPEA database. The enrichment analysis of HubGs revealed some crucial CC-causing biological processes (BPs), molecular functions (MFs) and cellular components (CCs) by involving hHubGs. The gene regulatory network (GRN) analysis identified four TFs proteins and three miRNAs as the key transcriptional and post-transcriptional regulators of hHubGs. Then, we identified hHubGs-guided top-ranked FDA-approved 10 candidate drugs and validated them against the state-of-the-arts independent receptors by molecular docking analysis. Finally, we investigated the binding stability of the top-ranked three candidate drugs (Docetaxel, Temsirolimus, Paclitaxel) by using 100 ns MD-based MM-PBSA simulations and observed their stable performance. Therefore the finding of this study might be the useful resources for CC diagnosis and therapies.

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.

Integrated bioinformatics and statistical approaches to explore molecular biomarkers for breast cancer diagnosis, prognosis and therapies

Integrated bioinformatics and statistical approaches are now playing the vital role in identifying potential molecular biomarkers more accurately in presence of huge number of alternatives for disease diagnosis, prognosis and therapies by reducing time and cost compared to the wet-lab based experimental procedures. Breast cancer (BC) is one of the leading causes of cancer related deaths for women worldwide. Several dry-lab and wet-lab based studies have identified different sets of molecular biomarkers for BC. But they did not compare their results to each other so much either computationally or experimentally. In this study, an attempt was made to propose a set of molecular biomarkers that might be more effective for BC diagnosis, prognosis and therapies, by using the integrated bioinformatics and statistical approaches. At first, we identified 190 differentially expressed genes (DEGs) between BC and control samples by using the statistical LIMMA approach. Then we identified 13 DEGs (AKR1C1, IRF9, OAS1, OAS3, SLCO2A1, NT5E, NQO1, ANGPT1, FN1, ATF6B, HPGD, BCL11A, and TP53INP1) as the key genes (KGs) by protein-protein interaction (PPI) network analysis. Then we investigated the pathogenetic processes of DEGs highlighting KGs by GO terms and KEGG pathway enrichment analysis. Moreover, we disclosed the transcriptional and post-transcriptional regulatory factors of KGs by their interaction network analysis with the transcription factors (TFs) and micro-RNAs. Both supervised and unsupervised learning's including multivariate survival analysis results confirmed the strong prognostic power of the proposed KGs. Finally, we suggested KGs-guided computationally more effective seven candidate drugs (NVP-BHG712, Nilotinib, GSK2126458, YM201636, TG-02, CX-5461, AP-24534) compared to other published drugs by cross-validation with the state-of-the-art alternatives top-ranked independent receptor proteins. Thus, our findings might be played a vital role in breast cancer diagnosis, prognosis and therapies.