Drug Targeting and Biomarkers in Head and Neck Cancers: Insights from Systems Biology Analyses

Potential Signature Therapeutic Biomarkers TOP2A, MAD2L1, and CDK1 in Colorectal Cancer: A Systems Biomedicine-Based Approach

Colorectal cancer is the third deadliest and fourth most diagnosed cancer. It is heterogeneously driven by varied mutations and mutagens, and thus, it is challenging for targeted therapy. The rapid advancement of high-throughput technology presents considerable opportunities for discovering new colon cancer biomarkers. In the present study, we have explored and identified the biomarkers based on molecular interactions. We curated cancer datasets that were not micro-dissected and performed gene expression analysis. The protein-protein interactions were curated, and a network was constructed for the up-regulated genes. The hub genes were analyzed using 12 different topological parameters. The correlation analysis selected TOP2A, CDK1, CCNB1, AURKA, and MAD2L1 as hub genes. Further, survival analysis was performed to determine the effectiveness of the hub gene on the patient's survival rate. Our findings explore various transcription factors such as E2F4, FOXM1, E2F6, MAX, and SIN3A, along with kinases CSNK2A1, MAPK14, CDK1, CDK4, and CDK2, as potential molecular signatures and aid researchers in understanding the pathophysiological mechanisms underlying CRC development and thus providing novel therapeutic and diagnostic recourse. Furthermore, investigating miRNAs, we focused on hsa-miR-215-5p, hsa-miR-192-5p, and hsa-miR-193b-3p due to their observed impact on a diverse set of colorectal cancer genes. Thereby, the current approach brings into light CRC- related genes at the RNA and protein levels that can potentially act as novel biomarkers opening doors to diagnostic and treatment purposes.

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

Identifying Biomarkers and Therapeutic Targets by Multiomic Analysis for HNSCC: Precision Medicine and Healthcare Management

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the major types of cancer, with 900,000 cases and over 400,000 deaths annually. It constitutes 3-4% of all cancers in Europe and western countries. As early diagnosis is the key to treating the disease, reliable biomarkers play an important role in the precision medicine of HNSCC. Despite treatments, the survival rate of cancer patients remains unchanged, and this is mainly due to the failure to detect the disease early. Thus, the objective of this study is to identify reliable biomarkers for head and neck cancers for better healthcare management. Methods: In this study, all available, curated human genes were screened for their expression against HNSCC TCGA patient samples using genomic and proteomic data by various bioinformatic approaches and datamining. Docking studies were performed using AutoDock or online virtual screening tools for identifying potential ligands. Results: Sixty genes were short-listed, and most of them show a consistently higher expression in head and neck patient samples at both the mRNA and the protein level. Irrespective of human papillomavirus (HPV) status, all of them show a higher expression in cancer samples. The higher expression of 30 genes shows adverse effects on patient survival. Out of the 60 genes, 12 genes have crystal structures and druggable potential. We show that genes such as GTF2H4, HAUS7, MSN, and MNDA could be targets of Pembrolizumab and Nivolumab, which are approved monoclonal antibodies for HNSCC. Conclusion: Sixty genes are identified as potential biomarkers for head and neck cancers based on their consistent and statistically significantly higher expression in patient samples. Four proteins have been identified as potential drug targets based on their crystal structure. However, the utility of these candidate genes has to be further tested using patient samples.

Identification of hub genes and drug candidates for NF2-related vestibular schwannoma by bioinformatics tools

Neurofibromatosis type 2 (NF2)-related vestibular schwannoma (NF2-VS) is a rare genetic disorder that results in bilateral acoustic neuromas. However, the exact pathogenesis of the disease is still unclear. This study aims to use bioinformatics analyses to identify potential hub genes and therapeutic. We retrieved the mRNA expression profiles (GSE108524 and GSE141801) of NF2-VS from the database, and selected the leading 25% genes with the most variance across samples for weighted correlation network analysis. Subsequently, we conducted gene ontology term and Kyoto Encyclopedia of Genes and Genomes signaling network enrichment analyses. The STRING database was employed for protein-protein interaction (PPI) axis construction. The mRNA-miRNA modulatory network was generated via the miRTarBase database. Differentially expressed genes (DEGs) were identified via the R package "limma" in both datasets, and hub genes were screened via intersection of common DEGs, candidate hub genes from the PPI axis, and candidate hub genes from the key module. Finally, common DEGs were uploaded onto the connectivity map database to determine drug candidates. Based on our observations, the blue module exhibited the most significant relation to NF2-VS, and it included the NF2 gene. Using enrichment analysis, we demonstrated that the blue modules were intricately linked to modulations of cell proliferation, migration, adhesion, junction, and actin skeleton. Overall, 356 common DEGs were screened in both datasets, and 33 genes carrying a degree > 15 were chosen as candidate hub genes in the PPI axis. Subsequently, 4 genes, namely, GLUL, CAV1, MYH11, and CCND1 were recognized as real hub genes. In addition, 10 drugs with enrichment scores < -0.7 were identified as drug candidates. Our conclusions offered a novel insight into the potential underlying mechanisms behind NF2-VS. These findings may facilitate the identification of novel therapeutic targets in the future.

Prognostic and immune correlation analysis of mitochondrial autophagy and aging-related genes in lung adenocarcinoma

PURPOSE

Mitophagy and aging (MiAg) are very important pathophysiological mechanisms contributing to tumorigenesis. MiAg-related genes have prognostic value in lung adenocarcinoma (LUAD). However, prognostic, and immune correlation studies of MiAg-related genes in LUAD are lacking.

METHODS

MiAg differentially expressed genes (DEGs) in LUAD were obtained from public sequencing datasets. A prognostic model including MiAg DEGs was constructed according to patients divided into low- and high-risk groups. Gene Ontology, gene set enrichment analysis, gene set variation analysis, CIBERSORT immune infiltration analysis, and clinical characteristic correlation analyses were performed for functional annotation and correlation of MiAgs with prognosis in patients with LUAD.

RESULTS

Seven MiAg DEGs of LUAD were identified: CAV1, DSG2, DSP, MYH11, NME1, PAICS, PLOD2, and the expression levels of these genes were significantly correlated (P < 0.05). The RiskScore of the MiAg DEG prognostic model demonstrated high predictive ability of overall survival of patients diagnosed with LUAD. Patients with high and low MiAg phenotypic scores exhibited significant differences in the infiltration levels of eight types of immune cells (P < 0.05). The multi-factor DEG regression model showed higher efficacy in predicting 5-year survival than 3- and 1-year survival of patients with LUAD.

CONCLUSIONS

Seven MiAg-related genes were identified to be significantly associated with the prognosis of patients diagnosed with LUAD. Moreover, the identified MiAg DEGs might affect the immunotherapy strategy of patients with LUAD.

Identification of a two metastasis-related prognostic signature in the process of predicting the survival of laryngeal squamous cell carcinoma

Metastasis is a major cause of treatment failure and poor outcomes in cancer patients. The data used in the current study was downloaded from TCGA and GEO databases. Differentially expressed metastasis-related genes were identified and the biological functions were implemented. Kaplan-Meier analysis univariate, and, multivariate Cox regression analyses were performed to identify robust prognostic biomarkers, followed by construction of the risk model and nomogram. Gene set enrichment analysis was performed to identify pathways enriched in low- and high-risk groups. POLR2J3 and MYH11 were treated as prognostic biomarkers in LSCC and the risk model was constructed. Receiver operating characteristic curves revealed the good performance of the risk model. A nomogram with high accuracy was constructed, as evidenced by calibration and decision curves. Moreover, we found that the expressions of POLR2J3 and MYH11 was significantly higher in metastasis tissues compared with those in non-metastasis tissues by RT-qPCR and IHC. Our study identified novel metastasis-related prognostic biomarkers in LSCC and constructed a unique nomogram for predicting the prognosis of LSCC patients. Moreover, we explored the related mechanisms of metastasis-related genes in regulating LSCC.

HPV16 status predicts potential protein biomarkers and therapeutics in head and neck squamous cell carcinoma

Human papillomavirus (HPV) infection, especially HPV16, is one of the causative factors for the development of head and neck squamous cell (HNSC) carcinoma. HPV-positive and HPV-negative HNSC patients differ significantly in their molecular profiles and clinical features, so they should be evaluated differently depending on their HPV status. Given the tremendous variation in HNSC cancers depending on HPV, our goal in this study was to present biomarkers and treatment options tailored to the patient's HPV status. Gene expression levels of HPV16-positive and -negative patients were used as proxies, and the differential interactome algorithm was employed to identify the differential interacting proteins (DIPs). By assessing the prognostic capabilities and druggabilities of DIPs and their interacting partners (DIP-centered modules), we introduce eight modules as potential biomarkers specialized for either positive or negative phenotype. Finally, raloxifene was repositioned for the first time as a drug candidate for the treatment of HPV16-positive HNSC patients.

Bioinformatics screening of colorectal-cancer causing molecular signatures through gene expression profiles to discover therapeutic targets and candidate agents

BACKGROUND

Detection of appropriate receptor proteins and drug agents are equally important in the case of drug discovery and development for any disease. In this study, an attempt was made to explore colorectal cancer (CRC) causing molecular signatures as receptors and drug agents as inhibitors by using integrated statistics and bioinformatics approaches.

METHODS

To identify the important genes that are involved in the initiation and progression of CRC, four microarray datasets (GSE9348, GSE110224, GSE23878, and GSE35279) and an RNA_Seq profiles (GSE50760) were downloaded from the Gene Expression Omnibus database. The datasets were analyzed by a statistical r-package of LIMMA to identify common differentially expressed genes (cDEGs). The key genes (KGs) of cDEGs were detected by using the five topological measures in the protein-protein interaction network analysis. Then we performed in-silico validation for CRC-causing KGs by using different web-tools and independent databases. We also disclosed the transcriptional and post-transcriptional regulatory factors of KGs by interaction network analysis of KGs with transcription factors (TFs) and micro-RNAs. Finally, we suggested our proposed KGs-guided computationally more effective candidate drug molecules compared to other published drugs by cross-validation with the state-of-the-art alternatives of top-ranked independent receptor proteins.

RESULTS

We identified 50 common differentially expressed genes (cDEGs) from five gene expression profile datasets, where 31 cDEGs were downregulated, and the rest 19 were up-regulated. Then we identified 11 cDEGs (CXCL8, CEMIP, MMP7, CA4, ADH1C, GUCA2A, GUCA2B, ZG16, CLCA4, MS4A12 and CLDN1) as the KGs. Different pertinent bioinformatic analyses (box plot, survival probability curves, DNA methylation, correlation with immune infiltration levels, diseases-KGs interaction, GO and KEGG pathways) based on independent databases directly or indirectly showed that these KGs are significantly associated with CRC progression. We also detected four TFs proteins (FOXC1, YY1, GATA2 and NFKB) and eight microRNAs (hsa-mir-16-5p, hsa-mir-195-5p, hsa-mir-203a-3p, hsa-mir-34a-5p, hsa-mir-107, hsa-mir-27a-3p, hsa-mir-429, and hsa-mir-335-5p) as the key transcriptional and post-transcriptional regulators of KGs. Finally, our proposed 15 molecular signatures including 11 KGs and 4 key TFs-proteins guided 9 small molecules (Cyclosporin A, Manzamine A, Cardidigin, Staurosporine, Benzo[A]Pyrene, Sitosterol, Nocardiopsis Sp, Troglitazone, and Riccardin D) were recommended as the top-ranked candidate therapeutic agents for the treatment against CRC.

CONCLUSION

The findings of this study recommended that our proposed target proteins and agents might be considered as the potential diagnostic, prognostic and therapeutic signatures for CRC.

Integration of Mendelian randomisation and systems biology models to identify novel blood-based biomarkers for stroke

Stroke is the second largest cause of mortality in the world. Genome-wide association studies (GWAS) have identified some genetic variants associated with stroke risk, but their putative functional causal genes are unknown. Hence, we aimed to identify putative functional causal gene biomarkers of stroke risk. We used a summary-based Mendelian randomisation (SMR) approach to identify the pleiotropic associations of genetically regulated traits (i.e., gene expression and DNA methylation) with stroke risk. Using SMR approach, we integrated cis-expression quantitative loci (cis-eQTLs) and cis-methylation quantitative loci (cis-mQTLs) data with GWAS summary statistics of stroke. We also utilised heterogeneity in dependent instruments (HEIDI) test to distinguish pleiotropy from linkage from the observed associations identified through SMR analysis. Our integrative SMR analyses and HEIDI test revealed 45 candidate biomarker genes (FDR < 0.05; P_HEIDI>0.01) that were pleiotropically or potentially causally associated with stroke risk. Of those candidate biomarker genes, 10 genes (HTRA1, PMF1, FBN2, C9orf84, COL4A1, BAG4, NEK6, SH2B3, SH3PXD2A, ACAD10) were differentially expressed in genome-wide blood transcriptomics data from stroke and healthy individuals (FDR<0.05). Functional enrichment analysis of the identified candidate biomarker genes revealed gene ontologies and pathways involved in stroke, including "cell aging", "metal ion binding" and "oxidative damage". Based on the evidence of genetically regulated expression of genes through SMR and directly measured expression of genes in blood, our integrative analysis suggests ten genes as blood biomarkers of stroke risk. Furthermore, our study provides a better understanding of the influence of DNA methylation on the expression of genes linked to stroke risk.

Robust identification of common genomic biomarkers from multiple gene expression profiles for the prognosis, diagnosis, and therapies of pancreatic cancer

Pancreatic cancer (PC) is one of the leading causes of cancer-related death globally. So, identification of potential molecular signatures is required for diagnosis, prognosis, and therapies of PC. In this study, we detected 71 common differentially expressed genes (cDEGs) between PC and control samples from four microarray gene-expression datasets (GSE15471, GSE16515, GSE71989, and GSE22780) by using robust statistical and machine learning approaches, since microarray gene-expression datasets are often contaminated by outliers due to several steps involved in the data generating processes. Then we detected 8 cDEGs (ADAM10, COL1A2, FN1, P4HB, ITGB1, ITGB5, ANXA2, and MYOF) as the PC-causing key genes (KGs) by the protein-protein interaction (PPI) network analysis. We validated the expression patterns of KGs between case and control samples by box plot analysis with the TCGA and GTEx databases. The proposed KGs showed high prognostic power with the random forest (RF) based prediction model and Kaplan-Meier-based survival probability curve. The KGs regulatory network analysis detected few transcriptional and post-transcriptional regulators for KGs. The cDEGs-set enrichment analysis revealed some crucial PC-causing molecular functions, biological processes, cellular components, and pathways that are associated with KGs. Finally, we suggested KGs-guided five repurposable drug molecules (Linsitinib, CX5461, Irinotecan, Timosaponin AIII, and Olaparib) and a new molecule (NVP-BHG712) against PC by molecular docking. The stability of the top three protein-ligand complexes was confirmed by molecular dynamic (MD) simulation studies. The cross-validation and some literature reviews also supported our findings. Therefore, the finding of this study might be useful resources to the researchers and medical doctors for diagnosis, prognosis and therapies of PC by the wet-lab validation.

Immune Cells in Head-and-Neck Tumor Microenvironments

Head-and-neck cancers constitute a heterogeneous group of aggressive tumors with high incidence and low survival rates, collectively being the sixth most prevalent cancer type globally. About 90% of head-and-neck cancers are classified as squamous cell carcinomas (HNSCC). The innate and adaptive immune systems, indispensable for anti-cancer immune surveillance, largely define the rates of HNSCC emergence and progression. HNSCC microenvironments harbor multiple cell types that infiltrate the tumors and interact both with tumor cells and among themselves. Gradually, tumor cells learn to manipulate the immune system, either by adapting their own immunogenicity or through the release of immunosuppressive molecules. These interactions continuously evolve and shape the tumor microenvironment, both structurally and functionally, facilitating angiogenesis, proliferation and metastasis. Our understanding of this evolution is directly related to success in the development of advanced therapies. This review focuses on the key mechanisms that rule HNSCC infiltration, featuring particular immune cell types and their roles in the pathogenesis. A close focus on the tumor-immunity interactions will help identify new immunotherapeutic targets in patients with HNSCC.

Meta-Data Analysis to Explore the Hub of the Hub-Genes That Influence SARS-CoV-2 Infections Highlighting Their Pathogenetic Processes and Drugs Repurposing

The pandemic of SARS-CoV-2 infections is a severe threat to human life and the world economic condition. Although vaccination has reduced the outspread, but still the situation is not under control because of the instability of RNA sequence patterns of SARS-CoV-2, which requires effective drugs. Several studies have suggested that the SARS-CoV-2 infection causing hub differentially expressed genes (Hub-DEGs). However, we observed that there was not any common hub gene (Hub-DEGs) in our analyses. Therefore, it may be difficult to take a common treatment plan against SARS-CoV-2 infections globally. The goal of this study was to examine if more representative Hub-DEGs from published studies by means of hub of Hub-DEGs (hHub-DEGs) and associated potential candidate drugs. In this study, we reviewed 41 articles on transcriptomic data analysis of SARS-CoV-2 and found 370 unique hub genes or studied genes in total. Then, we selected 14 more representative Hub-DEGs (AKT1, APP, CXCL8, EGFR, IL6, INS, JUN, MAPK1, STAT3, TNF, TP53, UBA52, UBC, VEGFA) as hHub-DEGs by their protein-protein interaction analysis. Their associated biological functional processes, transcriptional, and post-transcriptional regulatory factors. Then we detected hHub-DEGs guided top-ranked nine candidate drug agents (Digoxin, Avermectin, Simeprevir, Nelfinavir Mesylate, Proscillaridin, Linifanib, Withaferin, Amuvatinib, Atazanavir) by molecular docking and cross-validation for treatment of SARS-CoV-2 infections. Therefore, the findings of this study could be useful in formulating a common treatment plan against SARS-CoV-2 infections globally.

A Bi(OTf)3-Promoted Hydrosulfonylation of Alkenes with Sulfonyl Hydrazides: An Approach to Branched Sulfones

The Bi(OTf)₃-promoted hydrosulfonylation of alkenes with sulfonyl hydrazides to produce branched sulfones is reported, in which various branched sulfones (>40 examples) have been prepared in moderate to good yields. The gram-scale reaction and synthesis of the experimental inhibitor precursor showed the potential application. A preliminary mechanistic study revealed that double-bond migration to form the α,β-conjugated alkene is crucial for this transformation.

Construction and Validation of a UPR-Associated Gene Prognostic Model for Head and Neck Squamous Cell Carcinoma

Our study is aimed at constructing and validating a UPR-associated gene signature to predict HNSCC prognosis. We obtained 544 samples of RNA sequencing data and clinical characteristics from TCGA database and randomly grouped the samples into training and testing cohorts (1 : 1 ratio). After identifying 14 UPR-associated genes with LASSO and univariate Cox regression analysis, HNSCC samples were categorized into low-risk (LR) and high-risk (HR) subgroups depending on the risk score. Our analyses indicated that low-risk patients had a much better prognosis in the training and testing cohorts. To predict the HNSCC prognosis with the 14 UPR-associated gene signatures, we incorporated the UPR gene risk score, N stage, M stage, and age into a nomogram model. We further explored the sensitivity to anticancer drugs by using the IC50 analysis in two subgroups from the Cancer Genome Project database. The outcomes showed that the AKT inhibitor III and sorafenib were sensitive anticancer drugs in HR and LR patients, respectively. The immune cell infiltration analysis and GSEA provided strong evidence for elucidating the molecular mechanisms of UPR-associated genes affecting HNSCC. In conclusion, the UPR-associated gene risk score, N stage, M stage, and age can serve as a robust model for predicting prognosis and can improve decision-making at the individual patient level.

Statistics and network-based approaches to identify molecular mechanisms that drive the progression of breast cancer

Breast cancer (BC) is one of the most malignant tumors and the leading cause of cancer-related death in women worldwide. So, an in-depth investigation on the molecular mechanisms of BC progression is required for diagnosis, prognosis and therapies. In this study, we identified 127 common differentially expressed genes (cDEGs) between BC and control samples by analyzing five gene expression profiles with NCBI accession numbers GSE139038, GSE62931, GSE45827, GSE42568 and GSE54002, based-on two statistical methods LIMMA and SAM. Then we constructed protein-protein interaction (PPI) network of cDEGs through the STRING database and selected top-ranked 7 cDEGs (BUB1, ASPM, TTK, CCNA2, CENPF, RFC4, and CCNB1) as a set of key genes (KGs) by cytoHubba plugin in Cytoscape. Several BC-causing crucial biological processes, molecular functions, cellular components, and pathways were significantly enriched by the estimated cDEGs including at-least one KGs. The multivariate survival analysis showed that the proposed KGs have a strong prognosis power of BC. Moreover, we detected some transcriptional and post-transcriptional regulators of KGs by their regulatory network analysis. Finally, we suggested KGs-guided three repurposable candidate-drugs (Trametinib, selumetinib, and RDEA119) for BC treatment by using the GSCALite online web tool and validated them through molecular docking analysis, and found their strong binding affinities. Therefore, the findings of this study might be useful resources for BC diagnosis, prognosis and therapies.

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.

Disclosing Potential Key Genes, Therapeutic Targets and Agents for Non-Small Cell Lung Cancer: Evidence from Integrative Bioinformatics Analysis

Non-small-cell lung cancer (NSCLC) is considered as one of the malignant cancers that causes premature death. The present study aimed to identify a few potential novel genes highlighting their functions, pathways, and regulators for diagnosis, prognosis, and therapies of NSCLC by using the integrated bioinformatics approaches. At first, we picked out 1943 DEGs between NSCLC and control samples by using the statistical LIMMA approach. Then we selected 11 DEGs (CDK1, EGFR, FYN, UBC, MYC, CCNB1, FOS, RHOB, CDC6, CDC20, and CHEK1) as the hub-DEGs (potential key genes) by the protein–protein interaction network analysis of DEGs. The DEGs and hub-DEGs regulatory network analysis commonly revealed four transcription factors (FOXC1, GATA2, YY1, and NFIC) and five miRNAs (miR-335-5p, miR-26b-5p, miR-92a-3p, miR-155-5p, and miR-16-5p) as the key transcriptional and post-transcriptional regulators of DEGs as well as hub-DEGs. We also disclosed the pathogenetic processes of NSCLC by investigating the biological processes, molecular function, cellular components, and KEGG pathways of DEGs. The multivariate survival probability curves based on the expression of hub-DEGs in the SurvExpress web-tool and database showed the significant differences between the low- and high-risk groups, which indicates strong prognostic power of hub-DEGs. Then, we explored top-ranked 5-hub-DEGs-guided repurposable drugs based on the Connectivity Map (CMap) database. Out of the selected drugs, we validated six FDA-approved launched drugs (Dinaciclib, Afatinib, Icotinib, Bosutinib, Dasatinib, and TWS-119) by molecular docking interaction analysis with the respective target proteins for the treatment against NSCLC. The detected therapeutic targets and repurposable drugs require further attention by experimental studies to establish them as potential biomarkers for precision medicine in NSCLC treatment.

Multiomics Data Integration Identifies New Molecular Signatures for Abdominal Aortic Aneurysm and Aortic Occlusive Disease: Implications for Early Diagnosis, Prognosis, and Therapeutic Targets

Cardiovascular disease (CVD) is the leading cause of death among adults in developed countries. Among CVDs, abdominal aortic aneurysm (AAA) and aortic occlusive disease (AOD) are of great public health importance because of the high mortality rate in the elderly population. Despite significant molecular insights into AAA and AOD, the molecular mechanisms of these diseases remain unclear, and the current lack of robust diagnostic and prognostic biomarkers requires novel approaches to biomarker discovery and molecular targeting. In this study, we performed a comparative analysis of genome-wide expression data from patients with large AAA (n = 29), small AAA (n = 20), AOD (n = 9), and controls (n = 10). Specifically, we identified the differentially expressed genes and associated molecular pathways and biological processes (BPs) in each disease. Using a systems science approach, these data were linked to comprehensive human biological networks (i.e., protein-protein interaction, transcriptional regulatory, and metabolic networks) to identify molecular signatures of the salient mechanisms of AAA and AOD. Significant alterations in lipid metabolism and valine, leucine, and isoleucine metabolism, as well as neurodegenerative diseases and sex differences in the pathogenesis of AAA and AOD were identified. In the presence of aneurysm, size-dependent changes in lipid metabolism were observed. In addition, molecules and signaling pathways related to immunity, inflammation, infectious disease, and oxidative phosphorylation were identified in common. The results of the comparative and integrative analyzes revealed important clues to disease mechanisms and reporter molecules at various levels that warrant future development as potential prognostic biomarkers and putative therapeutic targets.

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

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is one of the most severe global pandemic due to its high pathogenicity and death rate starting from the end of 2019. Though there are some vaccines available against SAER-CoV-2 infections, we are worried about their effectiveness, due to its unstable sequence patterns. Therefore, beside vaccines, globally effective supporting drugs are also required for the treatment against SARS-CoV-2 infection. To explore commonly effective repurposable drugs for the treatment against different variants of coronavirus infections, in this article, an attempt was made to explore host genomic biomarkers guided repurposable drugs for SARS-CoV-1 infections and their validation with SARS-CoV-2 infections by using the integrated bioinformatics approaches. At first, we identified 138 differentially expressed genes (DEGs) between SARS-CoV-1 infected and control samples by analyzing high throughput gene-expression profiles to select drug target key receptors. Then we identified top-ranked 11 key DEGs (SMAD4, GSK3B, SIRT1, ATM, RIPK1, PRKACB, MED17, CCT2, BIRC3, ETS1 and TXN) as hub genes (HubGs) by protein-protein interaction (PPI) network analysis of DEGs highlighting their functions, pathways, regulators and linkage with other disease risks that may influence SARS-CoV-1 infections. The DEGs-set enrichment analysis significantly detected some crucial biological processes (immune response, regulation of angiogenesis, apoptotic process, cytokine production and programmed cell death, response to hypoxia and oxidative stress), molecular functions (transcription factor binding and oxidoreductase activity) and pathways (transcriptional mis-regulation in cancer, pathways in cancer, chemokine signaling pathway) that are associated with SARS-CoV-1 infections as well as SARS-CoV-2 infections by involving HubGs. The gene regulatory network (GRN) analysis detected some transcription factors (FOXC1, GATA2, YY1, FOXL1, TP53 and SRF) and micro-RNAs (hsa-mir-92a-3p, hsa-mir-155-5p, hsa-mir-106b-5p, hsa-mir-34a-5p and hsa-mir-19b-3p) as the key transcriptional and post- transcriptional regulators of HubGs, respectively. We also detected some chemicals (Valproic Acid, Cyclosporine, Copper Sulfate and arsenic trioxide) that may regulates HubGs. The disease-HubGs interaction analysis showed that our predicted HubGs are also associated with several other diseases including different types of lung diseases. Then we considered 11 HubGs mediated proteins and their regulatory 6 key TFs proteins as the drug target proteins (receptors) and performed their docking analysis with the SARS-CoV-2 3CL protease-guided top listed 90 anti-viral drugs out of 3410. We found Rapamycin, Tacrolimus, Torin-2, Radotinib, Danoprevir, Ivermectin and Daclatasvir as the top-ranked 7 candidate-drugs with respect to our proposed target proteins for the treatment against SARS-CoV-1 infections. Then, we validated these 7 candidate-drugs against the already published top-ranked 11 target proteins associated with SARS-CoV-2 infections by molecular docking simulation and found their significant binding affinity scores with our proposed candidate-drugs. Finally, we validated all of our findings by the literature review. Therefore, the proposed candidate-drugs might play a vital role for the treatment against different variants of SARS-CoV-2 infections with comorbidities, since the proposed HubGs are also associated with several comorbidities.

Computational identification of host genomic biomarkers highlighting their functions, pathways and regulators that influence SARS-CoV-2 infections and drug repurposing

The pandemic threat of COVID-19 has severely destroyed human life as well as the economy around the world. Although, the vaccination has reduced the outspread, but people are still suffering due to the unstable RNA sequence patterns of SARS-CoV-2 which demands supplementary drugs. To explore novel drug target proteins, in this study, a transcriptomics RNA-Seq data generated from SARS-CoV-2 infection and control samples were analyzed. We identified 109 differentially expressed genes (DEGs) that were utilized to identify 10 hub-genes/proteins (TLR2, USP53, GUCY1A2, SNRPD2, NEDD9, IGF2, CXCL2, KLF6, PAG1 and ZFP36) by the protein–protein interaction (PPI) network analysis. The GO functional and KEGG pathway enrichment analyses of hub-DEGs revealed some important functions and signaling pathways that are significantly associated with SARS-CoV-2 infections. The interaction network analysis identified 5 TFs proteins and 6 miRNAs as the key regulators of hub-DEGs. Considering 10 hub-proteins and 5 key TFs-proteins as drug target receptors, we performed their docking analysis with the SARS-CoV-2 3CL protease-guided top listed 90 FDA approved drugs. We found Torin-2, Rapamycin, Radotinib, Ivermectin, Thiostrepton, Tacrolimus and Daclatasvir as the top ranked seven candidate drugs. We investigated their resistance performance against the already published COVID-19 causing top-ranked 11 independent and 8 protonated receptor proteins by molecular docking analysis and found their strong binding affinities, which indicates that the proposed drugs are effective against the state-of-the-arts alternatives independent receptor proteins also. Finally, we investigated the stability of top three drugs (Torin-2, Rapamycin and Radotinib) by using 100 ns MD-based MM-PBSA simulations with the two top-ranked proposed receptors (TLR2, USP53) and independent receptors (IRF7, STAT1), and observed their stable performance. Therefore, the proposed drugs might play a vital role for the treatment against different variants of SARS-CoV-2 infections.

Integrative Multi-Omics Analysis Reveals Candidate Biomarkers for Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of cancer worldwide. Due to the lack of early detection and treatment, the survival rate of OSCC remains poor and the incidence of OSCC has not decreased during the past decades. To explore potential biomarkers and therapeutic targets for OSCC, we analyzed differentially expressed genes (DEGs) associated with OSCC using RNA sequencing technology. Methylation-regulated and differentially expressed genes (MeDEGs) of OSCC were further identified via an integrative approach by examining publicly available methylomic datasets together with our transcriptomic data. Protein-protein interaction (PPI) networks of MeDEGs were constructed and highly connected hub MeDEGs were identified from these PPI networks. Subsequently, expression and survival analyses of hub genes were performed using The Cancer Genome Atlas (TCGA) database and the Gene Expression Profiling Interactive Analysis (GEPIA) online tool. A total of 56 upregulated MeDEGs and 170 downregulated MeDEGs were identified in OSCC. Eleven hub genes with high degree of connectivity were picked out from the PPI networks constructed by those MeDEGs. Among them, the expression level of four hub genes (CTLA4, CDSN, ACTN2, and MYH11) were found to be significantly changed in the head and neck squamous carcinoma (HNSC) patients. Three hypomethylated hub genes (CTLA4, GPR29, and TNFSF11) and one hypermethylated hub gene (ISL1) were found to be significantly associated with overall survival (OS) of HNSC patients. Therefore, these hub genes may serve as potential DNA methylation biomarkers and therapeutic targets of OSCC.

Revisiting miRNA Association with Melanoma Recurrence and Metastasis from a Machine Learning Point of View

The diagnostic and prognostic value of miRNAs in cutaneous melanoma (CM) has been broadly studied and supported by advanced bioinformatics tools. From early studies using miRNA arrays with several limitations, to the recent NGS-derived miRNA expression profiles, an accurate diagnostic panel of a comprehensive pre-specified set of miRNAs that could aid timely identification of specific cancer stages is still elusive, mainly because of the heterogeneity of the approaches and the samples. Herein, we summarize the existing studies that report several miRNAs as important diagnostic and prognostic biomarkers in CM. Using publicly available NGS data, we analyzed the correlation of specific miRNA expression profiles with the expression signatures of known gene targets. Combining network analytics with machine learning, we developed specific non-linear classification models that could successfully predict CM recurrence and metastasis, based on two newly identified miRNA signatures. Subsequent unbiased analyses and independent test sets (i.e., a dataset not used for training, as a validation cohort) using our prediction models resulted in 73.85% and 82.09% accuracy in predicting CM recurrence and metastasis, respectively. Overall, our approach combines detailed analysis of miRNA profiles with heuristic optimization and machine learning, which facilitates dimensionality reduction and optimization of the prediction models. Our approach provides an improved prediction strategy that could serve as an auxiliary tool towards precision treatment.

Identification of a Novel Signature Predicting Overall Survival in Head and Neck Squamous Cell Carcinoma

Background: Head and neck squamous cell carcinoma (HNSCC) is a highly heterogeneous tumor with a high incidence and poor prognosis. Therefore, effective predictive models are needed to evaluate patient outcomes and optimize treatment. Methods: Robust Rank Aggregation (RRA) method was used to identify highly robust differentially-expressed genes (DEGs) between HNSCC and normal tissue in 9 GEO and TCGA datasets. Univariate Cox regression analysis and Lasso Cox regression analysis were performed to identify DEGs related to the Overall survival (OS) and to construct a prognostic gene signature (HNSCCSig). External validation was performed using GSE65858 dataset. Moreover, comprehensive bioinformatics analyses were used to identify the association between HNSCCSig and tumor immune environment. Results: A total of 257 reliable DEGs were identified by differentially analysis result of TCGA and GSE65858 datasets. The HNSCCSig including 7 mRNAs (SLURP1, SCARA5, CLDN10, MYH11, CXCL13, HLF, and ITGA3) were developed and validated to identify high-risk group who had a worse OS than low-risk group in TCGA and GSE65858 datasets. Cox regression analysis showed that the HNSCCSig could independently predict OS in both the TCGA and the GSE65858 datasets. Further research demonstrated that the infiltration bundance of CD8 + T cells, B cells, neutrophils, and NK cells were significantly lower in the high-risk group. A nomogram was also constructed by combining the HNSCCSig and clinical characters. Conclusion: We established and validated the HNSCCSig consisting of SLURP1, SCARA5, CLDN10, MYH11, CXCL13, HLF, and ITGA3. A nomogram combining HNSCCSig and some clinical parameters was constructed to identify high-risk HNSCC-patients with poor prognosis.

Peptidomics and proteogenomics: background, challenges and future needs

INTRODUCTION

With available genomic data and related information, it is becoming possible to better highlight mutations or genomic alterations associated with a particular disease or disorder. The advent of high-throughput sequencing technologies has greatly advanced diagnostics, prognostics, and drug development.

AREAS COVERED

Peptidomics and proteogenomics are the two post-genomic technologies that enable the simultaneous study of peptides and proteins/transcripts/genes. Both technologies add a remarkably large amount of data to the pool of information on various peptides associated with gene mutations or genome remodeling. Literature search was performed in the PubMed database and is up to date.

EXPERT OPINION

This article lists various techniques used for peptidomic and proteogenomic analyses. It also explains various bioinformatics workflows developed to understand differentially expressed peptides/proteins and their role in disease pathogenesis. Their role in deciphering disease pathways, cancer research, and biomarker discovery using biofluids is highlighted. Finally, the challenges and future requirements to overcome the current limitations for their effective clinical use are also discussed.

Interaction between DNMT3B and MYH11 via hypermethylation regulates gastric cancer progression

BACKGROUND

Gastric cancer (GC) has an unwelcoming prognosis when diagnosed at an advanced stage. The purpose of this study was to examine the expression of myosin heavy chain 11 (MYH11) in GC and mechanisms related.

METHODS

The MYH11 expression in GC was investigated via the SangerBox platform. MYH11 expression in GC tissues and cell lines was examined by immunohistochemistry, RT-qPCR, and western blot. The relationship between MYH11 expression and patients' prognosis was analyzed. The effects of MYH11 on the biological behaviors of GC cells were investigated by gain-of-function experiments. Bioinformatics analysis was used to find genes with relevance to MYH11 expression in GC. The relationship was verified by luciferase and ChIP-qPCR assays, followed by rescue assay validation. The causes of MYH11 downregulation in GC were verified by quantitative methylation-specific PCR. Finally, the effect of MYH11 on tumor growth was examined.

RESULTS

MYH11 was downregulated in GC and predicted poor prognoses. MYH11 reverted the malignant phenotype of GC cells. MYH11 repressed the TNFRSF14 expression by binding to the TNFRSF14 promoter. TNFRSF14 reversed the inhibitory effect of MYH11 on the malignant phenotype of GC cells. The methylation of the MYH11 promoter was elevated in GC, which was correlated with the elevated DNMT3B in GC. Overexpression of DNMT3B repressed transcription of MYH11 by promoting its methylation. Also, MYH11 upregulation inhibited tumor growth.

CONCLUSION

DNMT3B inhibits MYH11 expression by promoting its DNA methylation, thereby attenuating the repressive effect of MYH11 on the transcriptional of TNFRSF14 and promoting the progression of GC.

MicroRNA-1252-5p, regulated by Myb, inhibits invasion and epithelial-mesenchymal transition of pancreatic cancer cells by targeting NEDD9

MicroRNAs (miRNAs) are known to be involved in the development and progression of pancreatic cancer (PAC). The expression levels and roles of miR-1252-5p in PAC remain unclear. Quantitative real-time PCR and in situ hybridization were used to detect miR-1252-5p expression in PAC cells and human tissues. We studied the gain and loss of function of miR-1252-5p in the PAC cell lines in vitro and in vivo. The direct targets of miR-1252-5p were analyzed using public databases and a dual-luciferase reporter assay. Expression levels of miR-1252-5p are downregulated in PAC cell lines and tissue samples, and its expression is negatively associated with adverse clinical features and poor prognosis. In vitro and in vivo experiments show that miR-1252-5p overexpression inhibits the proliferation, migration, invasion, and epithelial-mesenchymal transition of PAC cells, and miR-1252-5p knockdown enhances these biological behaviors. MiR-1252-5p negatively regulates neural precursor cell expressed, developmentally downregulated 9 (NEDD9) by directly binding its 3'- untranslated region. Further mechanism research revealed that the SRC/STAT3 pathway is involved in miR-1252-5p/NEDD9 mediation of PAC's biological behaviors. We also verified that Myb inhibited miR-1252-5p by directly binding at its promoter. MiR-1252-5p may act as a tumor-suppressing miRNA in PAC and may be a potential therapeutic target for PAC patients.

Cancer Stem Cell Transcriptome Profiling Reveals Seed Genes of Tumorigenesis: New Avenues for Cancer Precision Medicine

Cancer stem-like cells (CSCs) possess the ability to self-renew and differentiate, and they are among the major factors driving tumorigenesis, metastasis, and resistance to chemotherapy. Therefore, it is critical to understand the molecular substrates of CSC biology so as to discover novel molecular biosignatures that distinguish CSCs and tumor cells. Here, we report new findings and insights by employing four transcriptome datasets associated with CSCs, with CSC and tumor samples from breast, lung, oral, and ovarian tissues. The CSC samples were analyzed to identify differentially expressed genes between CSC and tumor phenotypes. Through comparative profiling of expression levels in different cancer types, we identified 17 "seed genes" that showed a mutual differential expression pattern. We showed that these seed genes were strongly associated with cancer-associated signaling pathways and biological processes, the immune system, and the key cancer hallmarks. Further, the seed genes presented significant changes in their expression profiles in different cancer types and diverse mutation rates, and they also demonstrated high potential as diagnostic and prognostic biomarkers in various cancers. We report a number of seed genes that represent significant potential as "systems biomarkers" for understanding the pathobiology of tumorigenesis. Seed genes offer a new innovation avenue for potential applications toward cancer precision medicine in a broad range of cancers in oncology in the future.

Targeting the tumor stroma: integrative analysis reveal GATA2 and TORYAIP1 as novel prognostic targets in breast and ovarian cancer

Tumor stroma interaction is known to take a crucial role in cancer growth and progression. In the present study, it was performed gene expression analysis of stroma samples with ovarian and breast cancer through an integrative analysis framework to identify common critical biomolecules at multiomics levels. Gene expression datasets were statistically analyzed to identify common differentially expressed genes (DEGs) by comparing tumor stroma and normal stroma samples. The integrative analyses of DEGs indicated that there were 59 common core genes, which might be feasible to be potential marks for cancer stroma targeted strategies. Reporter molecules (i.e. receptor, transcription factors and miRNAs) were determined through a statistical test employing the hypergeometric probability density function. Afterward, the tumor microenvironment protein-protein interaction and the generic network were reconstructed by using identified reporter molecules and common core DEGs. Through a systems medicine approach, it was determined that hub biomolecules, AR, GATA2, miR-124, TOR1AIP1, ESR1, EGFR, STAT1, miR-192, GATA3, COL1A1, in tumor microenvironment generic network. These molecules were also identified as prognostic signatures in breast and ovarian tumor samples via survival analysis. According to literature searching, GATA2 and TORYAIP1 might represent potential biomarkers and candidate drug targets for the stroma targeted cancer therapy applications.

Identification and validation of a prognostic signature and combination drug therapy for immunotherapy of head and neck squamous cell carcinoma

Immunotherapy has become a promising therapeutic option for Head and neck squamous cell carcinoma (HNSC). However, only a small percentage of patients could benefit from it, and the overall prognosis was far from satisfactory. In this study, by comprehensively computational analyses of hundreds of HNSC samples, a prognostic signature composed of 13 immune-related genes (IRGs) was constructed. The results of the analyses in multiple datasets indicated that our signature had high predictive accuracy and could serve as an independent prognostic predictor. Based on this signature and multiple clinical variables, we also established a prognostic nomogram to quantitatively predict the survival risk of individual patients. Moreover, this signature could accurately predict survival, reflect the immune microenvironment, and predict immunotherapy efficacy among HNSC patients. Two potential drugs (doxorubicin and daunorubicin) were also identified via Connectivity Map and molecular docking, which could be used for HNSC combination therapy. Taken together, we developed and validated a robust IRG-based prognostic signature to monitor the prognosis of HNSC, which could provide a solid foundation for individualized cancer immunotherapy.

Identification of hub genes and small molecule therapeutic drugs related to breast cancer with comprehensive bioinformatics analysis

Breast cancer is one of the most common malignant tumors among women worldwide and has a high morbidity and mortality. This research aimed to identify hub genes and small molecule drugs for breast cancer by integrated bioinformatics analysis. After downloading multiple gene expression datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, 283 overlapping differentially expressed genes (DEGs) significantly enriched in different cancer-related functions and pathways were obtained using LIMMA, VennDiagram and ClusterProfiler packages of R. We then analyzed the topology of protein–protein interaction (PPI) network with overlapping DEGs and further obtained six hub genes (RRM2, CDC20, CCNB2, BUB1B, CDK1, and CCNA2) from the network via STRING and Cytoscape. Subsequently, we conducted genes expression verification, genetic alterations evaluation, immune infiltration prediction, clinicopathological parameters analysis, identification of transcriptional and post-transcriptional regulatory molecules, and survival analysis for these hub genes. Meanwhile, 29 possible drug candidates (e.g., Cladribine, Gallium nitrate, Alvocidib, 1β-hydroxyalantolactone, Berberine hydrochloride, Nitidine chloride) were identified from the DGIdb database and the GSE85871 dataset. In addition, some transcription factors and miRNAs (e.g., E2F1, PTTG1, TP53, ZBTB16, hsa-miR-130a-3p, hsa-miR-204-5p) targeting hub genes were identified as key regulators in the progression of breast cancer. In conclusion, our study identified six hub genes and 29 potential drug candidates for breast cancer. These findings may advance understanding regarding the diagnosis, prognosis and treatment of breast cancer.

Integrated Bioinformatics Analysis of mRNAs and miRNAs Identified Potential Biomarkers of Oral Squamous Cell Carcinoma

Background: Oral cancer is a frequently encountered neoplasm of the head and neck region, being the eighth most common type of human malignancy worldwide. Despite improvement in its control, morbidity and mortality, rates have improved little in the past decades. The present investigations about gene interaction and pathways still could not clear the appearance and development of oral squamous cell carcinoma (OSCC), completely. The aim of this study is to investigate the key genes and microRNAs interaction in OSCC. Materials and Methods: The microarray datasets GSE13601 and GSE98463, including mRNA and miRNA profiles, were extracted from the GEO database and were analyzed using GEO2R. Functional and pathway enrichment analyses were performed by using the DAVID database. The protein-protein interaction (PPI) network was constructed and analyzed using STRING database and Cytoscape software, respectively. Finally, miRDB was applied to predict the targets of the differentially expressed miRNAs (DEMs). Results: Totally, 97 differentially expressed genes (DEGs) were found in OSCC, including 66 up-regulated and 31 down-regulated genes. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that up-regulated genes were significantly enriched in movement of cell or subcellular component, cell adhesion, biological adhesion, cellular localization, apoptotic signaling pathway, while the downregulated genes were enriched in muscle system process and oxidation-reduction process. From the PPI network, the top 10 nodes with the highest degree were detected as hub genes. In addition, 18 DEMs were screened, which included 7 up-regulated and 11 down-regulated miRNAs. STAT1 was potentially targeted by three miRNAs, including has-miR- 6825-5P, has-miR-4495, and has-miR-5580-3P. Conclusion: The roles of DEMs such as hsa-mir-5580-3p in OSCC through interactions with DEGs CD44, ACLY, ACTR3, STAT1, LAMC2 and YWHAZ may offer a suitable candidate biomarker pattern for diagnosis, prognosis and treatment processes in OSCC.

Novel molecular signatures and potential therapeutics in renal cell carcinomas: Insights from a comparative analysis of subtypes

Renal cell carcinomas (RCCs) are among the highest causes of cancer mortality. Although transcriptome profiling studies in the last decade have made significant molecular findings on RCCs, effective diagnosis and treatment strategies have yet to be achieved due to lack of adequate screening and comparative profiling of RCC subtypes. In this study, a comparative analysis was performed on RNA-seq based transcriptome data from each RCC subtype, namely clear cell RCC (KIRC), papillary RCC (KIRP) and kidney chromophobe (KICH), and mutual or subtype-specific reporter biomolecules were identified at RNA, protein, and metabolite levels by the integration of expression profiles with genome-scale biomolecular networks. This approach revealed already-known biomarkers in RCCs as well as novel biomarker candidates and potential therapeutic targets. Our findings also pointed out the incorporation of the molecular mechanisms of KIRC and KIRP, whereas KICH was shown to have distinct molecular signatures. Furthermore, considering the Dipeptidyl Peptidase 4 (DPP4) receptor as a potential therapeutic target specific to KICH, several drug candidates such as ZINC6745464 were identified through virtual screening of ZINC molecules. In this study, we reported valuable data for further experimental and clinical efforts, since the proposed molecules have significant potential for screening and therapeutic purposes in RCCs.

Clinical and prognostic significance of MYH11 in lung cancer

Myosin heavy chain 11 (MYH11), encoded by the MYH11 gene, is a protein that participates in muscle contraction through the hydrolysis of adenosine triphosphate. Although previous studies have demonstrated that MYH11 gene expression levels are downregulated in several types of cancer, its expression levels have rarely been investigated in lung cancer. The present study aimed to explore the clinical significance and prognostic value of MYH11 expression levels in lung cancer and to further study the underlying molecular mechanisms of the function of this gene. The Oncomine database showed that the MYH11 expression levels were decreased in lung cancer compared with those noted in the normal lung tissue (P<0.05). Kaplan-Meier plotter results revealed that the decreased MYH11 expression levels were correlated with poor prognosis in lung cancer patients. Among the lung cancer cases with gene alteration of MYH11, mutation was the most common of all alteration types. Coexpedia and Metascape analyses revealed that the target genes were primarily enriched in ‘muscle contraction’, ‘contractile fiber part’, ‘actin cytoskeleton’ and the ‘adherens junction’. These results indicated that MYH11 is a potential novel drug target and prognostic indicator of lung cancer.

Discovering Biomarkers and Pathways Shared by Alzheimer's Disease and Ischemic Stroke to Identify Novel Therapeutic Targets

Background and objectives: Alzheimer's disease (AD) is a progressive neurodegenerative disease that results in severe dementia. Having ischemic strokes (IS) is one of the risk factors of the AD, but the molecular mechanisms that underlie IS and AD are not well understood. We thus aimed to identify common molecular biomarkers and pathways in IS and AD that can help predict the progression of these diseases and provide clues to important pathological mechanisms. Materials and Methods: We have analyzed the microarray gene expression datasets of IS and AD. To obtain robust results, combinatorial statistical methods were used to analyze the datasets and 26 transcripts (22 unique genes) were identified that were abnormally expressed in both IS and AD. Results: Gene Ontology (GO) and KEGG pathway analyses indicated that these 26 common dysregulated genes identified several altered molecular pathways: Alcoholism, MAPK signaling, glycine metabolism, serine metabolism, and threonine metabolism. Further protein-protein interactions (PPI) analysis revealed pathway hub proteins PDE9A, GNAO1, DUSP16, NTRK2, PGAM2, MAG, and TXLNA. Transcriptional and post-transcriptional components were then identified, and significant transcription factors (SPIB, SMAD3, and SOX2) found. Conclusions: Protein-drug interaction analysis revealed PDE9A has interaction with drugs caffeine, γ-glutamyl glycine, and 3-isobutyl-1-methyl-7H-xanthine. Thus, we identified novel putative links between pathological processes in IS and AD at transcripts levels, and identified possible mechanistic and gene expression links between IS and AD.

Identification of Prognostic Biomarker Signatures and Candidate Drugs in Colorectal Cancer: Insights from Systems Biology Analysis

Background and objectives: Colorectal cancer (CRC) is the second most common cause of cancer-related death in the world, but early diagnosis ameliorates the survival of CRC. This report aimed to identify molecular biomarker signatures in CRC. Materials and Methods: We analyzed two microarray datasets (GSE35279 and GSE21815) from the Gene Expression Omnibus (GEO) to identify mutual differentially expressed genes (DEGs). We integrated DEGs with protein–protein interaction and transcriptional/post-transcriptional regulatory networks to identify reporter signaling and regulatory molecules; utilized functional overrepresentation and pathway enrichment analyses to elucidate their roles in biological processes and molecular pathways; performed survival analyses to evaluate their prognostic performance; and applied drug repositioning analyses through Connectivity Map (CMap) and geneXpharma tools to hypothesize possible drug candidates targeting reporter molecules. Results: A total of 727 upregulated and 99 downregulated DEGs were detected. The PI3K/Akt signaling, Wnt signaling, extracellular matrix (ECM) interaction, and cell cycle were identified as significantly enriched pathways. Ten hub proteins (ADNP, CCND1, CD44, CDK4, CEBPB, CENPA, CENPH, CENPN, MYC, and RFC2), 10 transcription factors (ETS1, ESR1, GATA1, GATA2, GATA3, AR, YBX1, FOXP3, E2F4, and PRDM14) and two microRNAs (miRNAs) (miR-193b-3p and miR-615-3p) were detected as reporter molecules. The survival analyses through Kaplan–Meier curves indicated remarkable performance of reporter molecules in the estimation of survival probability in CRC patients. In addition, several drug candidates including anti-neoplastic and immunomodulating agents were repositioned. Conclusions: This study presents biomarker signatures at protein and RNA levels with prognostic capability in CRC. We think that the molecular signatures and candidate drugs presented in this study might be useful in future studies indenting the development of accurate diagnostic and/or prognostic biomarker screens and efficient therapeutic strategies in CRC.

Network-based approach to identify molecular signatures and therapeutic agents in Alzheimer's disease

Alzheimer's disease (AD) is a dynamic degeneration of the brain with progressive dementia. Considering the uncertainties in its molecular mechanism, in the present study, we employed network-based integrative analyses, and aimed to explore the key molecules and their associations with small drugs to identify potential biomarkers and therapeutic agents for the AD. First of all, we studied a transcriptome dataset and identified 1521 differentially expressed genes (DEGs). Integration of transcriptome data with protein-protein and transcriptional regulatory interactions resulted with central (hub) proteins (UBA52, RAC1, CREBBP, AR, RPS11, SMAD3, RPS6, RPL12, RPL15, and UBC), regulatory transcription factors (FOXC1, GATA2, YY1, FOXL1, NFIC, E2F1, USF2, SRF, PPARG, and JUN) and microRNAs (mir-335-5p, mir-26b-5p, mir-93-5p, mir-124-3p, mir-17-5p, mir-16-5p, mir-20a-5p, mir-92a-3p, mir-106b-5p, and mir-192-5p) as key signaling and regulatory molecules associated with transcriptional changes for the AD. Considering these key molecules as potential therapeutic targets and Connectivity Map (CMap) architecture, candidate small molecular agents (such as STOCK1N-35696) were identified as novel potential therapeutics for the AD. This study presents molecular signatures at RNA and protein levels which might be useful in increasing discernment of the molecular mechanisms, and potential drug targets and therapeutics to design effective treatment strategies for the AD.