Bioinformatics and machine learning methodologies to identify the effects of central nervous system disorders on glioblastoma progression

Neuroinflammation mediates the progression of neonate hypoxia-ischemia brain damage to Alzheimer's disease: a bioinformatics and experimental study

Background

Traumatic brain injury (TBI) can generally be divided into focal damage and diffuse damage, and neonate Hypoxia-Ischemia Brain Damage (nHIBD) is one of the causes of diffuse damage. Patients with nHIBD are at an increased risk of developing Alzheimer's disease (AD). However, the shared pathogenesis of patients affected with both neurological disorders has not been fully elucidated.

Purpose

We here aim to identify the shared molecular signatures between nHIBD and AD. We used an integrated analysis of the cortex gene expression data, targeting differential expression of genes related to the mechanisms of neurodegeneration and cognitive impairment following traumatic brain injury.

Methods

The gene expression profiles of Alzheimer's disease (GSE203206) and that of Neonate Hypoxia-Ischemia Brain Damage (GSE23317) were obtained from the Gene Expression Omnibus (GEO) database. After identifying the common differentially expressed genes (DEGs) of Alzheimer's disease and neonate Hypoxia-Ischemia Brain Damage by limma package analysis, five kinds of analyses were performed on them, namely Gene Ontology (GO) and pathway enrichment analysis, protein-protein interaction network, DEG-transcription factor interactions and DEG-microRNA interactions, protein-drug interactions and protein-disease association analysis, and gene-inflammation association analysis and protein-inflammation association analysis.

Results

In total, 12 common DEGs were identified including HSPB1, VIM, MVD, TUBB4A, AACS, ANXA6, DIRAS2, RPH3A, CEND1, KALM, THOP1, AREL1. We also identified 11 hub proteins, three central regulatory transcription factors, and three microRNAs encoded by the DEGs. Protein-drug interaction analysis showed that CYC1 and UQCRFS1 are associated with different drugs. Gene-disease association analysis shows Mammary Neoplasms, Neoplasm Metastasis, Schizophrenia, and Brain Ischemia diseases are the most relevant to the hub proteins we identified. Gene-inflammation association analysis shows that the hub gene AREL1 is related to inflammatory response, while the protein-inflammation association analysis shows that the hub proteins AKT1 and MAPK14 are related to inflammatory response.

Conclusion

This study provides new insights into the shared molecular mechanisms between AD and nHIBD. These common pathways and hub genes could potentially be used to design therapeutic interventions, reducing the likelihood of Alzheimer's disease development in survivors of neonatal Hypoxic-Ischemia brain injury.

Identification of RN7SK LncRNA as a novel biomarker in Alzheimer's disease using bioinformatics and expression analysis

Alzheimer's disease (AD) is a degenerative illness that accounts for the common type of dementia among adults over the age of 65. Despite extensive studies on the pathogenesis of the disease, early diagnosis of AD is still debatable. In this research, we performed bioinformatics approaches on the AD-related E-MTAB 6094 dataset to uncover new potential biomarkers for AD diagnosis. To achieve this, we performed in-depth in silico assays, including differentially expressed genes analysis, weighted gene co-expression network analyses, module-trait association analyses, gene ontology and pathway enrichment analyses, and hub genes network analyses. Finally, the expression of the identified candidate genes was evaluated in AD patients PBMC samples by qRT-PCR. Through computational analyses, we found that RN7SK LncRNA and its co-expressed genes of TNF, TNFAIP3, CCLT3, and FLT3 are from key genes in AD development that are associated with inflammatory responses. Our experimental validation revealed that RN7SK LncRNA and TNF were substantially up-regulated in AD samples (P = 0.006 and P = 0.023, respectively). Whereas, TNFAIP3 expression was significantly decreased (P = 0.016). However, the expression of CCL3 and FLT3 did not differ significantly between two groups (P = 0.396 and P = 0.521, respectively). In conclusion, in this study a novel LncRNA associated with AD pathogenesis were identified, which may provide new diagnostic biomarker for AD.

One path, two solutions: Network-based analysis identifies targetable pathways for the treatment of comorbid type II diabetes and neuropsychiatric disorders

Comorbid diseases complicate patient outcomes and escalate healthcare costs, necessitating the need for a deeper mechanistic understanding. Neuropsychiatric disorders (NPDs) such as Neurotic Disorder, Major Depression, Bipolar Disorder, Anxiety Disorder, and Schizophrenia significantly exacerbate Type 2 Diabetes Mellitus (DM2), often leading to suboptimal treatment outcomes. The neurobiological mechanisms underlying this comorbidity remain poorly understood. To address this gap, we developed a novel pathway-based network computational framework to identify critical shared disease mechanisms between DM2 and these five prevalent comorbid NPDs. Our approach involves reconstructing an integrated DM2 ∩ NPDs KEGG pathway-pathway network and employs two complementary analytical methods, including the "minimum path to comorbidity" method to identify the shortest path fostering comorbid development. This analysis uncovered shared pathways like the PI3K-Akt signaling pathway and highlighted key nodes such as calcium signaling, MAPK, estrogen signaling, and apoptosis pathways. Dysregulation of these pathways likely contributes to the development of DM2-NPDs comorbidity. These findings have significant clinical implications, as they identify promising therapeutic targets that could lead to more effective treatments addressing both DM2 and NPDs simultaneously. Our model not only elucidates the intricate molecular interactions driving this comorbidity but also identifies promising therapeutic targets, paving the way for innovative treatment strategies. Additionally, the framework developed in this study can be adapted to study other complex comorbid conditions, advancing personalized medicine for comorbidities and improving patient care.

Integrative Multi-Omics Analysis Reveals Critical Molecular Networks Linking Intestinal-System Diseases to Colorectal Cancer Progression

Background/Objectives: Colorectal cancer (CRC) frequently co-occurs with intestinal system diseases (ISDs), yet their molecular interplay remains poorly understood. We employed a comprehensive bioinformatics approach to elucidate shared genetic signatures and pathways between CRC and ISDs. Methods: We systematically analyzed 12 microarray and RNA-seq datasets encompassing 989 samples across seven ISDs and CRC. Differentially expressed genes (DEGs) were identified using Limma and DESeq2. Functional enrichment analysis was performed using clusterProfiler. Protein-protein interaction networks were constructed via STRING and visualized with Cytoscape to identify hub genes. Clinical significance of shared genes was further assessed through survival analysis and validated by immunohistochemistry staining of 30 paired CRC-normal tissue samples. Results: Integrating bioinformatics and machine learning approaches, we uncovered 160 shared DEGs (87 upregulated, 73 downregulated), which predominantly enriched cell metabolism, immune homeostasis, gut-brain communication, and inflammation pathways. Network analysis revealed nine key hub proteins linking CRC and ISDs, with seven upregulated (CD44, MYC, IL17A, CXCL1, FCGR3A, SPP1, and IL1A) and two downregulated (CXCL12 and CCL5). Survival analysis demonstrated the prognostic potential of these shared genes, while immunohistochemistry confirmed their differential expression in CRC tissues. Conclusions: Our findings unveil potential biomarkers and therapeutic targets, providing insights into ISD-influenced CRC progression and offering a robust foundation for improved diagnostic and treatment strategies in ISD-associated CRC.

GeneAI 3.0: powerful, novel, generalized hybrid and ensemble deep learning frameworks for miRNA species classification of stationary patterns from nucleotides

Due to the intricate relationship between the small non-coding ribonucleic acid (miRNA) sequences, the classification of miRNA species, namely Human, Gorilla, Rat, and Mouse is challenging. Previous methods are not robust and accurate. In this study, we present AtheroPoint's GeneAI 3.0, a powerful, novel, and generalized method for extracting features from the fixed patterns of purines and pyrimidines in each miRNA sequence in ensemble paradigms in machine learning (EML) and convolutional neural network (CNN)-based deep learning (EDL) frameworks. GeneAI 3.0 utilized five conventional (Entropy, Dissimilarity, Energy, Homogeneity, and Contrast), and three contemporary (Shannon entropy, Hurst exponent, Fractal dimension) features, to generate a composite feature set from given miRNA sequences which were then passed into our ML and DL classification framework. A set of 11 new classifiers was designed consisting of 5 EML and 6 EDL for binary/multiclass classification. It was benchmarked against 9 solo ML (SML), 6 solo DL (SDL), 12 hybrid DL (HDL) models, resulting in a total of 11 + 27 = 38 models were designed. Four hypotheses were formulated and validated using explainable AI (XAI) as well as reliability/statistical tests. The order of the mean performance using accuracy (ACC)/area-under-the-curve (AUC) of the 24 DL classifiers was: EDL > HDL > SDL. The mean performance of EDL models with CNN layers was superior to that without CNN layers by 0.73%/0.92%. Mean performance of EML models was superior to SML models with improvements of ACC/AUC by 6.24%/6.46%. EDL models performed significantly better than EML models, with a mean increase in ACC/AUC of 7.09%/6.96%. The GeneAI 3.0 tool produced expected XAI feature plots, and the statistical tests showed significant p-values. Ensemble models with composite features are highly effective and generalized models for effectively classifying miRNA sequences.

An integrated in-silico Pharmaco-BioInformatics approaches to identify synergistic effects of COVID-19 to HIV patients

BACKGROUND

With high inflammatory states from both COVID-19 and HIV conditions further result in complications. The ongoing confrontation between these two viral infections can be avoided by adopting suitable management measures.

PURPOSE

The aim of this study was to figure out the pharmacological mechanism behind apigenin's role in the synergetic effects of COVID-19 to the progression of HIV patients.

METHOD

We employed computer-aided methods to uncover similar biological targets and signaling pathways associated with COVID-19 and HIV, along with bioinformatics and network pharmacology techniques to assess the synergetic effects of apigenin on COVID-19 to the progression of HIV, as well as pharmacokinetics analysis to examine apigenin's safety in the human body.

RESULT

Stress-responsive, membrane receptor, and induction pathways were mostly involved in gene ontology (GO) pathways, whereas apoptosis and inflammatory pathways were significantly associated in the Kyoto encyclopedia of genes and genomes (KEGG). The top 20 hub genes were detected utilizing the shortest path ranked by degree method and protein-protein interaction (PPI), as well as molecular docking and molecular dynamics simulation were performed, revealing apigenin's strong interaction with hub proteins (MAPK3, RELA, MAPK1, EP300, and AKT1). Moreover, the pharmacokinetic features of apigenin revealed that it is an effective therapeutic agent with minimal adverse effects, for instance, hepatoxicity.

CONCLUSION

Synergetic effects of COVID-19 on the progression of HIV may still be a danger to global public health. Consequently, advanced solutions are required to give valid information regarding apigenin as a suitable therapeutic agent for the management of COVID-19 and HIV synergetic effects. However, the findings have yet to be confirmed in patients, suggesting more in vitro and in vivo studies.

Implication of lncRNA ZBED3-AS1 downregulation in acquired resistance to Temozolomide and glycolysis in glioblastoma

Temozolomide (TMZ) is the recommended drug for glioblastoma (GBM) treatment, but its clinical effect is restricted due to drug resistance. This research studies the effects of long non-coding RNA (lncRNA) ZBED3-AS1 and its related molecules on acquired TMZ resistance in glioblastoma (GBM). ZBED3-AS1 was identified to be downregulated in TMZ-resistant GBM cells by analyzing GSE113510 and GSE100736 datasets. ZBED3-AS1 downregulation was detected in TMZ-resistant GBM tissues and cell lines (U251/TMZ and U87/TMZ). ZBED3-AS1 knockdown promoted, whereas its overexpression suppressed TMZ resistance, viability and mobility, and glycolytic activity of TMZ-resistant cells. ZBED3-AS1 bound to Spi-1 proto-oncogene (SPI1) but did not affect its expression. Instead, it blocked SPI1-mediated transcriptional activation of thrombomodulin (THBD). SPI1 and THBD increased TMZ resistance and glycolysis in TMZ-resistant cells. Either ZBED3-AS1 overexpression or SPI1 knockdown in U87/TMZ cells blocked the growth of orthotopic and subcutaneous xenograft tumors in nude mice. In conclusion, this study demonstrates that ZBED3-AS1 downregulation and THBD activation is linked to increased TMZ resistance and glycolysis in GBM cells.

Bioinformatics and System Biological Approaches for the Identification of Genetic Risk Factors in the Progression of Cardiovascular Disease

Background

Cardiovascular disease (CVD) is the combination of coronary heart disease, myocardial infarction, rheumatic heart disease, and peripheral vascular disease of the heart and blood vessels. It is one of the leading deadly diseases that causes one-third of the deaths yearly in the globe. Additionally, the risk factors associated with it make the situation more complex for cardiovascular patients, which lead them towards mortality, but the genetic association between CVD and its risk factors is not clearly explored in the global literature. We addressed this issue and explored the linkage between CVD and its risk factors.

Methods

We developed an analytical approach to reveal the risk factors and their linkages with CVD. We used GEO microarray datasets for the CVD and other risk factors in this study. We performed several analyses including gene expression analysis, diseasome analysis, protein-protein interaction (PPI) analysis, and pathway analysis for discovering the relationship between CVD and its risk factors. We also examined the validation of our study using gold benchmark databases OMIM, dbGAP, and DisGeNET.

Results

We observed that the number of 32, 17, 53, 70, and 89 differentially expressed genes (DEGs) is overlapped between CVD and its risk factors of hypertension (HTN), type 2 diabetes (T2D), hypercholesterolemia (HCL), obesity, and aging, respectively. We identified 10 major hub proteins (FPR2, TNF, CXCL8, CXCL1, IL1B, VEGFA, CYBB, PTGS2, ITGAX, and CCR5), 12 significant functional pathways, and 11 gene ontological pathways that are associated with CVD. We also found the connection of CVD with its risk factors in the gold benchmark databases. Our experimental outcomes indicate a strong association of CVD with its risk factors of HTN, T2D, HCL, obesity, and aging.

Conclusions

Our computational approach explored the genetic association of CVD with its risk factors by identifying the significant DEGs, hub proteins, and signaling and ontological pathways. The outcomes of this study may be further used in the lab-based analysis for developing the effective treatment strategies of CVD.

Integrated Analysis of Gene Co-Expression Network and Prediction Model Indicates Immune-Related Roles of the Identified Biomarkers in Sepsis and Sepsis-Induced Acute Respiratory Distress Syndrome

Sepsis is a series of clinical syndromes caused by immunological response to severe infection. As the most important and common complication of sepsis, acute respiratory distress syndrome (ARDS) is associated with poor outcomes and high medical expenses. However, well-described studies of analysis-based researches, especially related bioinformatics analysis on revealing specific targets and underlying molecular mechanisms of sepsis and sepsis-induced ARDS (sepsis/se-ARDS), still remain limited and delayed despite the era of data-driven medicine. In this report, weight gene co-expression network based on data from a public database was constructed to identify the key modules and screen the hub genes. Functional annotation by enrichment analysis of the modular genes also demonstrated the key biological processes and signaling pathway; among which, extensive immune-involved enrichment was remarkably associated with sepsis/se-ARDS. Based on the differential expression analysis, least absolute shrink and selection operator, and multivariable logistic regression analysis of the screened hub genes, SIGLEC9, TSPO, CKS1B and PTTG3P were identified as the candidate biomarkers for the further analysis. Accordingly, a four-gene-based model for diagnostic prediction assessment was established and then developed by sepsis/se-ARDS risk nomogram, whose efficiency was verified by calibration curves and decision curve analyses. In addition, various machine learning algorithms were also applied to develop extra models based on the four genes. Receiver operating characteristic curve analysis proved the great diagnostic and predictive performance of these models, and the multivariable logistic regression of the model was still found to be the best as further verified again by the internal test, training, and external validation cohorts. During the development of sepsis/se-ARDS, the expressions of the identified biomarkers including SIGLEC9, TSPO, CKS1B and PTTG3P were all regulated remarkably and generally exhibited notable correlations with the stages of sepsis/se-ARDS. Moreover, the expression levels of these four genes were substantially correlated during sepsis/se-ARDS. Analysis of immune infiltration showed that multiple immune cells, neutrophils and monocytes in particular, might be closely involved in the process of sepsis/se-ARDS. Besides, SIGLEC9, TSPO, CKS1B and PTTG3P were considerably correlated with the infiltration of various immune cells including neutrophils and monocytes during sepsis/se-ARDS. The discovery of relevant gene co-expression network and immune signatures might provide novel insights into the pathophysiology of sepsis/se-ARDS.

Integrated Bioinformatics Analysis and Verification of Gene Targets for Myocardial Ischemia-Reperfusion Injury

Background

Myocardial ischemia-reperfusion injury (MIRI) has become a thorny and unsolved clinical problem. The pathological mechanisms of MIRI are intricate and unclear, so it is of great significance to explore potential hub genes and search for some natural products that exhibit potential therapeutic efficacy on MIRI via targeting the hub genes.

Methods

First, the differential expression genes (DEGs) from GSE58486, GSE108940, and GSE115568 were screened and integrated via a robust rank aggregation algorithm. Then, the hub genes were identified and verified by the functional experiment of the MIRI mice. Finally, natural products with protective effects against MIRI were retrieved, and molecular docking simulations between hub genes and natural products were performed.

Results

230 integrated DEGs and 9 hub genes were identified. After verification, Emr1, Tyrobp, Itgb2, Fcgr2b, Cybb, and Fcer1g might be the most significant genes during MIRI. A total of 75 natural products were discovered. Most of them (especially araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin) showed good ability to bind the hub genes.

Conclusions

Emr1, Tyrobp, Itgb2, Fcgr2b, Cybb, and Fcer1g might be critical in the pathological process of MIRI, and the natural products (araloside C, glycyrrhizic acid, ophiopogonin D, polyphyllin I, and punicalagin) targeting these hub genes exhibited potential therapeutic efficacy on MIRI. Our findings provided new insights to explore the mechanism and treatments for MIRI and revealed new therapeutic targets for natural products with protective properties against MIRI.

Bioinformatics and Network-based Approaches for Determining Pathways, Signature Molecules, and Drug Substances connected to Genetic Basis of Schizophrenia etiology

Knowledge of heterogeneous etiology and pathophysiology of schizophrenia (SZP) is reasonably inadequate and non-deterministic due to its inherent complexity and underlying vast dynamics related to genetic mechanisms. The evolution of large-scale transcriptome-wide datasets and subsequent development of relevant, robust technologies for their analyses show promises toward elucidating the genetic basis of disease pathogenesis, its early risk prediction, and predicting drug molecule targets for therapeutic intervention. In this research, we have scrutinized the genetic basis of SZP through functional annotation and network-based system biology approaches. We have determined 96 overlapping differentially expressed genes (DEGs) from 2 microarray datasets and subsequently identified their interconnecting networks to reveal transcriptome signatures like hub proteins (FYN, RAD51, SOCS3, XIAP, AKAP13, PIK3C2A, CBX5, GATA3, EIF3K, and CDKN2B), transcription factors and miRNAs. In addition, we have employed gene set enrichment to highlight significant gene ontology (e.g., positive regulation of microglial cell activation) and relevant pathways (such as axon guidance and focal adhesion) interconnected to the genes associated with SZP. Finally, we have suggested candidate drug substances like Luteolin HL60 UP as a possible therapeutic target based on these key molecular signatures.

Network based systems biology approach to identify diseasome and comorbidity associations of Systemic Sclerosis with cancers

Systemic Sclerosis (SSc) is an autoimmune disease associated with changes in the skin's structure in which the immune system attacks the body. A recent meta-analysis has reported a high incidence of cancer prognosis including lung cancer (LC), leukemia (LK), and lymphoma (LP) in patients with SSc as comorbidity but its underlying mechanistic details are yet to be revealed. To address this research gap, bioinformatics methodologies were developed to explore the comorbidity interactions between a pair of diseases. Firstly, appropriate gene expression datasets from different repositories on SSc and its comorbidities were collected. Then the interconnection between SSc and its cancer comorbidities was identified by applying the developed pipelines. The pipeline was designed as a generic workflow to demonstrate a premise comorbid condition that integrate regarding gene expression data, tissue/organ meta-data, Gene Ontology (GO), Molecular pathways, and other online resources, and analyze them with Gene Set Enrichment Analysis (GSEA), Pathway enrichment and Semantic Similarity (SS). The pipeline was implemented in R and can be accessed through our Github repository: https://github.com/hiddenntreasure/comorbidity. Our result suggests that SSc and its cancer comorbidities share differentially expressed genes, functional terms (gene ontology), and pathways. The findings have led to a better understanding of disease pathways and our developed methodologies may be applied to any set of diseases for finding any association between them. This research may be used by physicians, researchers, biologists, and others.

Systems Biology and Bioinformatics approach to Identify blood based signatures molecules and drug targets of patient with COVID-19

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection results in the development of a highly contagious respiratory ailment known as new coronavirus disease (COVID-19). Despite the fact that the prevalence of COVID-19 continues to rise, it is still unclear how people become infected with SARS-CoV-2 and how patients with COVID-19 become so unwell. Detecting biomarkers for COVID-19 using peripheral blood mononuclear cells (PBMCs) may aid in drug development and treatment. This research aimed to find blood cell transcripts that represent levels of gene expression associated with COVID-19 progression. Through the development of a bioinformatics pipeline, two RNA-Seq transcriptomic datasets and one microarray dataset were studied and discovered 102 significant differentially expressed genes (DEGs) that were shared by three datasets derived from PBMCs. To identify the roles of these DEGs, we discovered disease-gene association networks and signaling pathways, as well as we performed gene ontology (GO) studies and identified hub protein. Identified significant gene ontology and molecular pathways improved our understanding of the pathophysiology of COVID-19, and our identified blood-based hub proteins TPX2, DLGAP5, NCAPG, CCNB1, KIF11, HJURP, AURKB, BUB1B, TTK, and TOP2A could be used for the development of therapeutic intervention. In COVID-19 subjects, we discovered effective putative connections between pathological processes in the transcripts blood cells, suggesting that blood cells could be used to diagnose and monitor the disease’s initiation and progression as well as developing drug therapeutics.

Exploration of the core protein network under endometriosis symptomatology using a computational approach

BACKGROUND

Endometriosis is defined by implantation and invasive growth of endometrial tissue in extra-uterine locations causing heterogeneous symptoms, and a unique clinical picture for each patient. Understanding the complex biological mechanisms underlying these symptoms and the protein networks involved may be useful for early diagnosis and identification of pharmacological targets.

METHODS

In the present study, we combined three approaches (i) a text-mining analysis to perform a systematic search of proteins over existing literature, (ii) a functional enrichment analysis to identify the biological pathways in which proteins are most involved, and (iii) a protein-protein interaction (PPI) network to identify which proteins modulate the most strongly the symptomatology of endometriosis.

RESULTS

Two hundred seventy-eight proteins associated with endometriosis symptomatology in the scientific literature were extracted. Thirty-five proteins were selected according to degree and betweenness scores criteria. The most enriched biological pathways associated with these symptoms were (i) Interleukin-4 and Interleukin-13 signaling (p = 1.11 x 10^-16), (ii) Signaling by Interleukins (p = 1.11 x 10^-16), (iii) Cytokine signaling in Immune system (p = 1.11 x 10^-16), and (iv) Interleukin-10 signaling (p = 5.66 x 10^-15).

CONCLUSION

Our study identified some key proteins with the ability to modulate endometriosis symptomatology. Our findings indicate that both pro- and anti-inflammatory biological pathways may play important roles in the symptomatology of endometriosis. This approach represents a genuine systemic method that may complement traditional experimental studies. The current data can be used to identify promising biomarkers for early diagnosis and potential therapeutic targets.

Comprehensive Analysis of Alteration Landscape and Its Clinical Significance of Mitochondrial Energy Metabolism Pathway-Related Genes in Lung Cancers

Background

Mitochondria are the energy factories of cells. The abnormality of mitochondrial energy metabolism pathways is closely related to the occurrence and development of lung cancer. The abnormal genes in mitochondrial energy metabolism pathways might be the novel targets and biomarkers to diagnose and treat lung cancers.

Method

Genes in major mitochondrial energy metabolism pathways were obtained from the KEGG database. The transcriptomic, mutation, and clinical data of lung cancers were obtained from The Cancer Genome Atlas (TCGA) database. Genes and clinical biomarkers were mined that affected lung cancer survival. Gene enrichment analysis was performed with ClusterProfiler and the gene set enrichment analysis (GSEA). STRING database and Cytoscape were used for protein-protein interaction (PPI) analysis. The diagnostic biomarker pattern of lung cancer was optimized, and its accuracy was verified with 10-fold cross-validation. The four genes screened by logistic regression model were verified by western blot in 5 pairs of lung cancer specimens collected in hospital.

Results

In total, 188 mitochondrial energy metabolism pathway-related genes (MMRGs) were included in this study. GSEA analysis found that MMRGs in the lung cancer group were mainly enriched in the metabolic pathway of oxidative phosphorylation and electron respiratory transport chain compared to the control group. Age did not affect the mutation frequency of MMRGs. Comparative analysis of these 188 MMRGs identified 43 differentially expressed MMRGs (24 upregulated and 19 downregulated) in the lung cancer group compared to the control group. The survival analysis of these 43 differentially expressed MMRGs found that the survival time was better in the low-expressed GAPDHS group than that in the high-expressed GAPDHS group of lung cancers. The advanced age, high expression of GAPDHS, low expressions of ACSBG1 and CYP4A11, and ACOX3 mutation were biomarkers of poor prognosis in lung cancers. PPI analysis showed that proteins such as GAPDH and GAPDHS interacted with many proteins in mitochondrial metabolic pathways. A four-MMRG-signature model (y = 0.0069∗ACADL - 0.001∗ALDH18A1 - 0.0405∗CPT1B + 0.0008∗PPARG - 1.625) was established to diagnose lung cancer with the accuracy up to 98.74%, AUC value up to 0.992, and a missed diagnosis rate of only 0.6%. Western blotting showed that ALDH18A1 and CPT1B proteins were significantly overexpressed in the lung cancer group (p < 0.05), and ACADL and PPARG proteins were slightly underexpressed in the lung cancer group (p < 0.05), which were consistent with the results of their corresponding mRNA expressions.

Conclusion

Mitochondrial energy metabolism pathway alterations are the important hallmarks of lung cancer. Age did not increase the risk of MMRG mutation. High expression of GAPDHS, low expression of ACSBG1, low expression of CYP4A11, mutated ACOX3, and old age predict a poor prognosis of lung cancer. Four differentially expressed MMRGs (ACADL, ALDH18A1, CPT1B, and PPARG) established a logistic regression model, which could effectively diagnose lung cancer. At the protein level, ALDH18A1 and CPT1B were significantly upregulated, and ACADL and PPARG were slightly underexpressed, in the lung cancer group compared to the control group, which were consistent with the results of their corresponding mRNA expressions.

Bioinformatics and machine learning approach identifies potential drug targets and pathways in COVID-19

Current coronavirus disease-2019 (COVID-19) pandemic has caused massive loss of lives. Clinical trials of vaccines and drugs are currently being conducted around the world; however, till now no effective drug is available for COVID-19. Identification of key genes and perturbed pathways in COVID-19 may uncover potential drug targets and biomarkers. We aimed to identify key gene modules and hub targets involved in COVID-19. We have analyzed SARS-CoV-2 infected peripheral blood mononuclear cell (PBMC) transcriptomic data through gene coexpression analysis. We identified 1520 and 1733 differentially expressed genes (DEGs) from the GSE152418 and CRA002390 PBMC datasets, respectively (FDR < 0.05). We found four key gene modules and hub gene signature based on module membership (MMhub) statistics and protein-protein interaction (PPI) networks (PPIhub). Functional annotation by enrichment analysis of the genes of these modules demonstrated immune and inflammatory response biological processes enriched by the DEGs. The pathway analysis revealed the hub genes were enriched with the IL-17 signaling pathway, cytokine-cytokine receptor interaction pathways. Then, we demonstrated the classification performance of hub genes (PLK1, AURKB, AURKA, CDK1, CDC20, KIF11, CCNB1, KIF2C, DTL and CDC6) with accuracy >0.90 suggesting the biomarker potential of the hub genes. The regulatory network analysis showed transcription factors and microRNAs that target these hub genes. Finally, drug-gene interactions analysis suggests amsacrine, BRD-K68548958, naproxol, palbociclib and teniposide as the top-scored repurposed drugs. The identified biomarkers and pathways might be therapeutic targets to the COVID-19.