Decoding temporal heterogeneity in NSCLC through machine learning and prognostic model construction.

BACKGROUND: Non-small cell lung cancer (NSCLC) is a prevalent and heterogeneous disease with significant genomic variations between the early and advanced stages. The identification of key genes and pathways driving NSCLC tumor progression is critical for improving the diagnosis and treatment outcomes of this disease. METHODS: In this study, we conducted single-cell transcriptome analysis on 93,406 cells from 22 NSCLC patients to characterize malignant NSCLC cancer cells. Utilizing cNMF, we classified these cells into distinct modules, thus identifying the diverse molecular profiles within NSCLC. Through pseudotime analysis, we delineated temporal gene expression changes during NSCLC evolution, thus demonstrating genes associated with disease progression. Using the XGBoost model, we assessed the significance of these genes in the pseudotime trajectory. Our findings were validated by using transcriptome sequencing data from The Cancer Genome Atlas (TCGA), supplemented via LASSO regression to refine the selection of characteristic genes. Subsequently, we established a risk score model based on these genes, thus providing a potential tool for cancer risk assessment and personalized treatment strategies. RESULTS: We used cNMF to classify malignant NSCLC cells into three functional modules, including the metabolic reprogramming module, cell cycle module, and cell stemness module, which can be used for the functional classification of malignant tumor cells in NSCLC. These findings also indicate that metabolism, the cell cycle, and tumor stemness play important driving roles in the malignant evolution of NSCLC. We integrated cNMF and XGBoost to select marker genes that are indicative of both early and advanced NSCLC stages. The expression of genes such as CHCHD2, GAPDH, and CD24 was strongly correlated with the malignant evolution of NSCLC at the single-cell data level. These genes have been validated via histological data. The risk score model that we established (represented by eight genes) was ultimately validated with GEO data. CONCLUSION: In summary, our study contributes to the identification of temporal heterogeneous biomarkers in NSCLC, thus offering insights into disease progression mechanisms and potential therapeutic targets. The developed workflow demonstrates promise for future applications in clinical practice.

Crosstalk of different cell-death patterns predicts prognosis and drug sensitivity in glioma.

BACKGROUND: Glioma is a malignant brain tumor originating from glial cells, and there still a challenge to accurately predict the prognosis. Programmed cell death (PCD) plays a key role in tumorigenesis and immune response. However, the crosstalk and potential role of various PCDs in prognosis and tumor microenvironment remains unknown. Therefore, we comprehensively discussed the relationship between different models of PCD and the prognosis of glioma and provided new ideas for the optimal targeted therapy of glioma. MATERIALS AND METHODS: We compared and analyzed the role of 14 PCD patterns on the prognosis from different levels. We constructed the cell death risk score (CDRS) index and conducted a comprehensive analysis of CDRS and TME characteristics, clinical characteristics, and drug response. RESULTS: Effects of different PCDs at the genomic, functional, and immune microenvironment levels were discussed. CDRS index containing 6 gene signatures and a nomogram were established. High CDRS is associated with a worse prognosis. Through transcriptome and single-cell data, we found that patients with high CDRS showed stronger immunosuppressive characteristics. Moreover, the high-CDRS group was resistant to the traditional glioma chemotherapy drug Vincristine, but more sensitive to the Temozolomide and the clinical experimental drug Bortezomib. In addition, we identified 19 key potential therapeutic targets during malignant differentiation of tumor cells. CONCLUSION: Overall, we provide the first systematic description of the role of 14 PCDs in glioma. A new CDRS model was built to predict the prognosis and to provide a new idea for the targeted therapy of glioma.

Predicting drug synergy using a network propagation inspired machine learning framework.

Combination therapy is a promising strategy for cancers, increasing therapeutic options and reducing drug resistance. Yet, systematic identification of efficacious drug combinations is limited by the combinatorial explosion caused by a large number of possible drug pairs and diseases. At present, machine learning techniques have been widely applied to predict drug combinations, but most studies rely on the response of drug combinations to specific cell lines and are not entirely satisfactory in terms of mechanism interpretability and model scalability. Here, we proposed a novel network propagation-based machine learning framework to predict synergistic drug combinations. Based on the topological information of a comprehensive drug-drug association network, we innovatively introduced an affinity score between drug pairs as one of the features to train machine learning models. We applied network-based strategy to evaluate their therapeutic potential to different cancer types. Finally, we identified 17 specific-, 21 general- and 40 broad-spectrum antitumor drug combinations, in which 69% drug combinations were validated by vitro cellular experiments, 83% drug combinations were validated by literature reports and 100% drug combinations were validated by biological function analyses. By quantifying the network relationships between drug targets and cancer-related driver genes in the human protein-protein interactome, we show the existence of four distinct patterns of drug-drug-disease relationships. We also revealed that 32 biological pathways were correlated with the synergistic mechanism of broad-spectrum antitumor drug combinations. Overall, our model offers a powerful scalable screening framework for cancer treatments.

Multi-level characteristics recognition of cancer core therapeutic targets and drug screening for a broader patient population.

Introduction: Target therapy for cancer cell mutation has brought attention to several challenges in clinical applications, including limited therapeutic targets, less patient benefits, and susceptibility to acquired due to their clear biological mechanisms and high specificity in targeting cancers with specific mutations. However, the identification of truly lethal synthetic lethal therapeutic targets for cancer cells remains uncommon, primarily due to compensatory mechanisms. Methods: In our pursuit of core therapeutic targets (CTTs) that exhibit extensive synthetic lethality in cancer and the corresponding potential drugs, we have developed a machine-learning model that utilizes multiple levels and dimensions of cancer characterization. This is achieved through the consideration of the transcriptional and post-transcriptional regulation of cancer-specific genes and the construction of a model that integrates statistics and machine learning. The model incorporates statistics such as Wilcoxon and Pearson, as well as random forest. Through WGCNA and network analysis, we identify hub genes in the SL network that serve as CTTs. Additionally, we establish regulatory networks for non-coding RNA (ncRNA) and drug-target interactions. Results: Our model has uncovered 7277 potential SL interactions, while WGCNA has identified 13 gene modules. Through network analysis, we have identified 30 CTTs with the highest degree in these modules. Based on these CTTs, we have constructed networks for ncRNA regulation and drug targets. Furthermore, by applying the same process to lung cancer and renal cell carcinoma, we have identified corresponding CTTs and potential therapeutic drugs. We have also analyzed common therapeutic targets among all three cancers. Discussion: The results of our study have broad applicability across various dimensions and histological data, as our model identifies potential therapeutic targets by learning multidimensional complex features from known synthetic lethal gene pairs. The incorporation of statistical screening and network analysis further enhances the confidence in these potential targets. Our approach provides novel theoretical insights and methodological support for the identification of CTTs and drugs in diverse types of cancer.

Integrating multiple machine learning algorithms for prognostic prediction of gastric cancer based on immune-related lncRNAs.

Background: Long non-coding RNAs (lncRNAs) play an important role in the immune regulation of gastric cancer (GC). However, the clinical application value of immune-related lncRNAs has not been fully developed. It is of great significance to overcome the challenges of prognostic prediction and classification of gastric cancer patients based on the current study. Methods: In this study, the R package ImmLnc was used to obtain immune-related lncRNAs of The Cancer Genome Atlas Stomach Adenocarcinoma (TCGA-STAD) project, and univariate Cox regression analysis was performed to find prognostic immune-related lncRNAs. A total of 117 combinations based on 10 algorithms were integrated to determine the immune-related lncRNA prognostic model (ILPM). According to the ILPM, the least absolute shrinkage and selection operator (LASSO) regression was employed to find the major lncRNAs and develop the risk model. ssGSEA, CIBERSORT algorithm, the R package maftools, pRRophetic, and clusterProfiler were employed for measuring the proportion of immune cells among risk groups, genomic mutation difference, drug sensitivity analysis, and pathway enrichment score. Results: A total of 321 immune-related lncRNAs were found, and there were 26 prognostic immune-related lncRNAs. According to the ILPM, 18 of 26 lncRNAs were selected and the risk score (RS) developed by the 18-lncRNA signature had good strength in the TCGA training set and Gene Expression Omnibus (GEO) validation datasets. Patients were divided into high- and low-risk groups according to the median RS, and the low-risk group had a better prognosis, tumor immune microenvironment, and tumor signature enrichment score and a higher metabolism, frequency of genomic mutations, proportion of immune cell infiltration, and antitumor drug resistance. Furthermore, 86 differentially expressed genes (DEGs) between high- and low-risk groups were mainly enriched in immune-related pathways. Conclusion: The ILPM developed based on 26 prognostic immune-related lncRNAs can help in predicting the prognosis of patients suffering from gastric cancer. Precision medicine can be effectively carried out by dividing patients into high- and low-risk groups according to the RS.

Exploring Precise Medication Strategies for OSCC Based on Single-Cell Transcriptome Analysis from a Dynamic Perspective.

At present, most patients with oral squamous cell carcinoma (OSCC) are in the middle or advanced stages at the time of diagnosis. Advanced OSCC patients have a poor prognosis after traditional therapy, and the complex heterogeneity of OSCC has been proven to be one of the main reasons. Single-cell sequencing technology provides a powerful tool for dissecting the heterogeneity of cancer. However, most of the current studies at the single-cell level are static, while the development of cancer is a dynamic process. Thus, understanding the development of cancer from a dynamic perspective and formulating corresponding therapeutic measures for achieving precise treatment are highly necessary, and this is also one of the main study directions in the field of oncology. In this study, we combined the static and dynamic analysis methods based on single-cell RNA-Seq data to comprehensively dissect the complex heterogeneity and evolutionary process of OSCC. Subsequently, for clinical practice, we revealed the association between cancer heterogeneity and the prognosis of patients. More importantly, we pioneered the concept of pseudo-time score of patients, and we quantified the levels of heterogeneity based on the dynamic development process to evaluate the relationship between the score and the survival status at the same stage, finding that it is closely related to the prognostic status. The pseudo-time score of patients could not only reflect the tumor status of patients but also be used as an indicator of the effects of drugs on the patients so that the medication strategy can be adjusted on time. Finally, we identified candidate drugs and proposed precision medication strategies to control the condition of OSCC in two respects: treatment and blocking.

Identification of Warning Transition Points from Hepatitis B to Hepatocellular Carcinoma Based on Mutation Accumulation for the Early Diagnosis and Potential Drug Treatment of HBV-HCC.

The accumulation of multiple genetic mutations is essential during the occurrence and development of hepatocellular carcinoma induced by hepatitis B (HBV-HCC), but understanding their cooperative effects and identifying the warning transition point from hepatitis B to HCC are challenges. In the genomic analysis of somatic mutations of the patient with HBV-HCC in a patient-specific protein-protein interaction (ps-PPI) network, we find mutation influence can propagate along the ps-PPI network. Therefore, in the article, we got the mutation cluster as a new research unit using the Random Walks with Restarts algorithm that is used to describe the efficient boundary of mutation influences. The connection of mutation cluster leads to dysregulation of signaling pathways corresponding to HCC, while dysregulated signaling pathways accumulate gradually and experience a process from quantitative to qualitative changes including a critical mutation cluster called transition point (TP) from hepatitis B to HCC. Moreover, two subtypes of HCC patients with different prognosis and their corresponding biological and clinical characteristics were identified according to TP. The poor prognosis HCC subtype was associated with significant metabolic pathway dysregulation and lower immune cell infiltration, while we also identified several preventive drugs to block the transformation of hepatitis B to hepatocellular carcinoma. The network-level study integrated multiomics data not only showed the sequence of multiple somatic mutations and their cooperative effect but also identified the warning transition point in HCC tumorigenesis for each patient. Our study provides new insight into exploring the cooperative molecular mechanism of chronic inflammatory malignancy in the liver and lays the foundation for the development of new approaches for early prediction and diagnosis of hepatocellular carcinoma and personalized targeted therapy.

Precision treatment exploration of breast cancer based on heterogeneity analysis of lncRNAs at the single-cell level.

BACKGROUND: Breast cancer (BC) is a complex disease with high heterogeneity, which often leads to great differences in treatment results. Current common molecular typing method is PAM50, which shows positive results for precision medicine; however, room for improvement still remains because of the different prognoses of subtypes. Therefore, in this article, we used lncRNAs, which are more tissue-specific and developmental stage-specific than other RNAs, as typing markers and combined single-cell expression profiles to retype BC, to provide a new method for BC classification and explore new precise therapeutic strategies based on this method. METHODS: Based on lncRNA expression profiles of 317 single cells from 11 BC patients, SC3 was used to retype BC, and differential expression analysis and enrichment analysis were performed to identify biological characteristics of new subtypes. The results were validated for survival analysis using data from TCGA. Then, the downstream regulatory genes of lncRNA markers of each subtype were searched by expression correlation analysis, and these genes were used as targets to screen therapeutic drugs, thus proposing new precision treatment strategies according to the different subtype compositions of patients. RESULTS: Seven lncRNA subtypes and their specific biological characteristics are obtained. Then, 57 targets and 210 drugs of 7 subtypes were acquired. New precision medicine strategies were proposed according to the different compositions of patient subtypes. CONCLUSIONS: For patients with different subtype compositions, we propose a strategy to select different drugs for different patients, which means using drugs targeting multi subtype or combinations of drugs targeting a single subtype to simultaneously kill different cancer cells by personalized treatment, thus reducing the possibility of drug resistance and even recurrence.

New insights for precision treatment of glioblastoma from analysis of single-cell lncRNA expression.

INTRODUCTION: Glioblastoma (GBM) is a complex disease with high intratumoral heterogeneity, understanding the molecular characteristics of intratumoral heterogeneity accurately is the basis for precision treatment. Although the existing typing strategy based on tumor molecular characteristics has a positive effect, there is still room for improvement, which is mainly because the traditional typing is completed based on the sequencing data of tissue samples, that is, the obtained data are the average level of patient tumor tissues, masking the intratumoral heterogeneity of a single patient and cannot reflect the real level of patient tumor cells. At present, cancer molecular typing is mostly performed based on transcriptome (RNA-seq) without considering lncRNA molecules that are also tissue-specific and developmental stage-specific. Therefore, in this study, we used lncRNAs as typing markers and combined single-cell expression profiles to retype glioblastoma, providing new ideas for GBM molecular typing, and further analyzed the shortcomings of traditional therapies at the singlecell level based on typing results and proposed new precise therapeutic insights. METHODS: We downloaded GBM single-cell sequencing data from GSE84465 and performed a series of preprocessing. The intratumoral heterogeneity of patients at the single-cell level was revealed using t-SNE, and the room for improvement of the existing traditional histotyping method was revealed using heat map and density curve. Subsequently, to validate the feasibility of lncRNA typing, we compared the similarities and differences of expression patterns between lncRNAs and mRNAs in GBM cells. Then, we used the R package "Seurat" to perform unsupervised clustering of GBM cells for re-typing and performed a detailed analysis of the biological characteristics of each subtype, including differentially expressed lncRNAs and marker lncRNAs. For validation, we performed survival analysis on GBM tissue data from the TCGA database to reveal the impact of different subtypes on patient survival prognosis. Eventually, based on the results, we screened the therapeutic drugs of each subtype by targeting the downstream regulatory genes of lncRNAs and proposed a new precision therapeutic strategy. RESULTS: GBM has significant intratumoral heterogeneity at the single-cell level, with more than one traditional subtype highly expressed in each patient, which reflects the shortcomings of traditional histotyping. LncRNAs and mRNAs have similar expression patterns in GBM cells, and the expression coefficient of variation of lncRNAs is higher than that of mRNAs, meaning that lncRNAs will better reflect the intratumoral heterogeneity. GBM was reclassified into four subtypes by unsupervised clustering, with different subtypes having different biological characteristics. Survival analysis showed that patients with different subtype compositions had different prognostic outcomes, so different subtypes had different effects on patient prognosis. Based on this, 47 drugs were screened for treatment. There are both shared and unique drugs between different subtypes. A new precision treatment strategy was proposed: for patients with different subtypes, in addition to the combination of drugs targeting single subtype, drugs targeting multiple subtypes can also be selected. CONCLUSION: Intratumoral heterogeneity may lead to poor prognosis or recurrence after treatment, and more precise typing of GBM can be performed based on single-cell lncRNA expression profiles. The biological characteristics possessed by different subtypes will have different effects on patients, such as survival time. For different subtypes, there are both drugs targeting single subtype and drugs targeting multiple subtypes, and we prefer drugs targeting multiple subtypes because this strategy can maximize medication efficiency and reduce the types of medication to reduce risks and side effects.

Computational drug repositioning based on the relationships between substructure-indication.

At present, computational methods for drug repositioning are mainly based on the whole structures of drugs, which limits the discovery of new functions due to the similarities between local structures of drugs. In this article, we, for the first time, integrated the features of chemical-genomics (substructure-domain) and pharmaco-genomics (domain-indication) based on the assumption that drug-target interactions are mediated by the substructures of drugs and the domains of proteins to identify the relationships between substructure-indication and establish a drug-substructure-indication network for predicting all therapeutic effects of tested drugs through only information on the substructures of drugs. In total, 83 205 drug-indication relationships with different correlation scores were obtained. We used three different verification methods to indicate the accuracy of the method and the reliability of the scoring system. We predicted all indications of olaparib using our method, including the known antitumor effect and unknown antiviral effect verified by literature, and we also discovered the inhibitory mechanism of olaparib toward DNA repair through its specific sub494 (o = C-C: C), as it participates in the low synthesis of the poly subfunction of the apoptosis pathway (hsa04210) by inhibiting the Inositol 1,4,5-trisphosphate receptor(s) (ITPRs) and hydrolyzing poly (ADP ribose) polymerases. ElectroCardioGrams of four drugs (quinidine, amiodarone, milrinone and fosinopril) demonstrated the effect of anti-arrhythmia. Unlike previous studies focusing on the overall structures of drugs, our research has great potential in the search for more therapeutic effects of drugs and in predicting all potential effects and mechanisms of a drug from the local structural similarity.

Identification of Early Warning Signals at the Critical Transition Point of Colorectal Cancer Based on Dynamic Network Analysis.

Colorectal cancer (CRC) is one of the leading causes of cancer-related death worldwide. Due to the lack of early diagnosis methods and warning signals of CRC and its strong heterogeneity, the determination of accurate treatments for CRC and the identification of specific early warning signals are still urgent problems for researchers. In this study, the expression profiles of cancer tissues and the expression profiles of tumor-adjacent tissues in 28 CRC patients were combined into a human protein-protein interaction (PPI) network to construct a specific network for each patient. A network propagation method was used to obtain a mutant giant cluster (GC) containing more than 90% of the mutation information of one patient. Next, mutation selection rules were applied to the GC to mine the mutation sequence of driver genes in each CRC patient. The mutation sequences from patients with the same type CRC were integrated to obtain the mutation sequences of driver genes of different types of CRC, which provide a reference for the diagnosis of clinical CRC disease progression. Finally, dynamic network analysis was used to mine dynamic network biomarkers (DNBs) in CRC patients. These DNBs were verified by clinical staging data to identify the critical transition point between the pre-disease state and the disease state in tumor progression. Twelve known drug targets were found in the DNBs, and 6 of them have been used as targets for anticancer drugs for clinical treatment. This study provides important information for the prognosis, diagnosis and treatment of CRC, especially for pre-emptive treatments. It is of great significance for reducing the incidence and mortality of CRC.

CHG: A Systematically Integrated Database of Cancer Hallmark Genes.

BACKGROUND: The analysis of cancer diversity based on a logical framework of hallmarks has greatly improved our understanding of the occurrence, development and metastasis of various cancers. METHODS: We designed Cancer Hallmark Genes (CHG) database which focuses on integrating hallmark genes in a systematic, standard way and annotates the potential roles of the hallmark genes in cancer processes. Following the conceptual criteria description of hallmark function the keywords for each hallmark were manually selected from the literature. Candidate hallmark genes collected were derived from 301 pathways of KEGG database by Lucene and manually corrected. RESULTS: Based on the variation data, we finally identified the hallmark genes of various types of cancer and constructed CHG. And we also analyzed the relationships among hallmarks and potential characteristics and relationships of hallmark genes based on the topological structures of their networks. We manually confirm the hallmark gene identified by CHG based on literature and database. We also predicted the prognosis of breast cancer, glioblastoma multiforme and kidney papillary cell carcinoma patients based on CHG data. CONCLUSIONS: In summary, CHG, which was constructed based on a hallmark feature set, provides a new perspective for analyzing the diversity and development of cancers.

Bioinformatics profiling utilized a nine immune-related long noncoding RNA signature as a prognostic target for pancreatic cancer.

PURPOSE: To identify an immune-related long noncoding RNA (lncRNA) signature with potential prognostic value for patients with pancreatic cancer. METHODS: Pancreatic cancer samples with available clinical information and whole genomic mRNA expression data obtained from The Cancer Genome Atlas (TCGA) were enrolled in our research. The immune score of each sample was calculated according to the expression level of immune-related genes and used to identify the most promising immune-related lncRNAs. According to the risk score developed from screened immune-related lncRNAs, the high- and low-risk groups were separated on the basis of the median risk score. The prediction reliability was further evaluated in the validation set and combination set. Both gene set enrichment analysis (GSEA) and principal component analysis (PCA) were performed for functional annotation, and the microenvironment cell population record was applied to evaluate the immune composition and purity of the tumor. RESULTS: A cohort of 176 samples was included in this study. A total of 163 immune-related lncRNAs were collected according to Pearson correlation analyses between immune score and lncRNA expression |R| > 0.5, P < 0.01). Nine immune-related lncRNAs (AL138966.2, AL133520.1, AC142472.1, AC127024.5, AC116913.1, AC083880.1, AC124016.1, AC008443.5, and AC092171.5) with the most significant prognostic values (P < 0.01) were identified. In the training set, it was observed that patients in the low-risk group showed longer overall survival (OS) than those in the high-risk group (P < 0.001); meanwhile, similar results were found in the validation set, combination set and various stratified sets (P < 0.05, P < 0.001, P < 0.05, respectively). Moreover, the signature was identified as an independent prognostic factor and significantly associated with the OS of pancreatic cancer. The area under curve (AUC) of the receiver operating characteristic curve (ROC curve) for the nine lncRNA signature in predicting the 2-year survival rate was 0.703. In addition, the low-risk and high-risk groups displayed different distributed patterns in PCA and different immune statuses in the GSEA. The signature indicated decreased purity of the tumor by implying a lower proportion of cancer cells along with an increasing enrichment of fibroblasts, myeloid dendritic cells, and monocytic lineage cells. CONCLUSIONS: Our research suggests that the immune-related lncRNA signature possesses latent prognostic value for patients with pancreatic cancer and may provide new information for immunological research and treatment in pancreatic cancer.

Discovering Proangiogenic Drugs in Ischemic Stroke Based on the Relationship between Protein Domain and Drug Substructure.

As an important protective mechanism against cerebral ischemia, angiogenesis has become a topic of interest in the treatment of ischemic stroke with the challenge that few drugs promote angiogenesis. Previous studies of the identification of drug-target interactions mainly focused on the overall structures of drugs and proteins, which limited the discovery of novel structure drugs. In this article, we proposed a new strategy for discovering proangiogenic drugs based on the assumption that drug-protein interaction is mediated by substructure and domain. First, we identified substructure-domain relationships according to the known drug-protein interactions and established the drug-substructure-domain-protein relationships of genes that are proangiogenic in brain tissue and expressed significantly during ischemic stroke. Then we quantified the intensity of interaction between each drug and each protein. Finally, we obtained 540 interactive relationships between 238 drugs and 54 genes, establishing a drug-gene network with two patterns of independent and complex drug-gene interactions. Both of the patterns facilitated finding not only drugs with the same overall structure but also drugs with a different overall structure based on the same or a similar protein spectrum. In addition, we analyzed the mechanism of action of each predicted drug and extracted drugs with similar mechanisms. In vitro, our results showed that azelnidipine, azilsartan, lercanidipine, nafcillin, and vortioxetine enhanced bEnd.3 cell proliferation, migration, tube formation, and the expression of angiogenic marker PCNA. Azelnidipine, azilsartan, lercanidipine, and nafcillin increased the level of expression of proangiogenic factor VEGF. Unlike previous studies focusing on the overall structures of drugs, our research highlighted local structural similarity, which has great potential in the search for more proangiogenic drugs with novel structures.

Synthetic Lethality-based Identification of Targets for Anticancer Drugs in the Human Signaling Network.

Chemotherapy agents can cause serious adverse effects by attacking both cancer tissues and normal tissues. Therefore, we proposed a synthetic lethality (SL) concept-based computational method to identify specific anticancer drug targets. First, a 3-step screening strategy (network-based, frequency-based and function-based screening) was proposed to identify the SL gene pairs by mining 697 cancer genes and the human signaling network, which had 6306 proteins and 62937 protein-protein interactions. The network-based screening was composed of a stability score constructed using a network information centrality measure (the average shortest path length) and the distance-based screening between the cancer gene and the non-cancer gene. Then, the non-cancer genes were extracted and annotated using drug-target interaction and drug description information to obtain potential anticancer drug targets. Finally, the human SL data in SynLethDB, the existing drug sensitivity data and text-mining were utilized for target validation. We successfully identified 2555 SL gene pairs and 57 potential anticancer drug targets. Among them, CDK1, CDK2, PLK1 and WEE1 were verified by all three aspects and could be preferentially used in specific targeted therapy in the future.

Predicting hepatocellular carcinoma through cross-talk genes identified by risk pathways.

Hepatocellular carcinoma (HCC) is the most frequent type of liver cancer with poor survival rate and high mortality. Despite efforts on the mechanism of HCC, new molecular markers are needed for exact diagnosis, evaluation and treatment. Here, we combined transcriptome of HCC with networks and pathways to identify reliable molecular markers. Through integrating 249 differentially expressed genes with syncretic protein interaction networks, we constructed a HCC-specific network, from which we further extracted 480 pivotal genes. Based on the cross-talk between the enriched pathways of the pivotal genes, we finally identified a HCC signature of 45 genes, which could accurately distinguish HCC patients with normal individuals and reveal the prognosis of HCC patients. Among these 45 genes, 15 showed dysregulated expression patterns and a part have been reported to be associated with HCC and/or other cancers. These findings suggested that our identified 45 gene signature could be potential and valuable molecular markers for diagnosis and evaluation of HCC.

Biological evaluation and energetic analyses of novel GSK-3ß inhibitors.

Glycogen synthase kinase-3 beta (GSK-3ß) is involved in multiple signaling pathways. Consistent with its critical roles in normal cells, abnormalities in GSK-3ß activity have been implicated in diabetes, heart disease, Parkinson disease, and Alzheimer's disease. In this study, a series of new scaffolds of small molecule inhibitors of GSK-3ß were identified by virtual screening and bioassay. Candidates that adhere to drug-like criteria from a virtual library of compounds were tested using computational docking studies. Twenty selected compounds were tested, which led to the discovery of two hits. Compound 14 (IC50 = 8.48 µM) and compound 19 (IC50 = 2.19 µM) were identified with high affinity. Molecular dynamics (MD) simulations, in conjunction with molecular mechanics/Poisson-Boltzmann surface area binding free-energy analysis, were employed to gain insight into the binding modes and energetics of GSK-3ß inhibitors. The detailed analysis of molecular dynamics results shows that Ile62, Val70, Tyr134, and Leu188 in GSK-3ß are key residues responsible to the binding of compound 14 and compound 19. Importantly, our results also validated this combined virtual screening and biophysical technique approach to discovery kinase inhibitors, which may be applied for future inhibitor discovery work for GSK-3ß.

Complex integrated analysis of lncRNAs-miRNAs-mRNAs in oral squamous cell carcinoma.

OBJECTIVES: This study aims to reveal regulatory network of lncRNAs-miRNAs-mRNAs in oral squamous cell carcinoma (OSCC) through gene expression data. MATERIAL AND METHODS: Differentially expressed lncRNAs, miRNAs and mRNAs (cut-off: False discovery rate (FDR)<0.05 and |fold change|>1.5) were unveiled by package edgeR of R. Cox regression analysis was performed to screen prognostic factors in OSCC related with overall survival (OS) and relapse-free survival (RFS). Protein-protein interaction (PPI) network was constructed for differentially expressed mRNAs using BioGRID, HPRD and DIP. Key hub genes were identified from top 100 differentially expressed mRNAs ranked by betweenness centrality using recursive feature elimination. LncRNA-miRNA and miRNA-mRNA regulatory network were constructed and combined into ceRNAs regulatory network. Gene ontology biological terms and Kyoto Encyclopedia of Genes and Genomes pathways were identified using Fisher's exact test. RESULTS: A total of 929 differentially expressed mRNAs, 23 differentially expressed lncRNAs and 29 differentially expressed miRNAs were identified. 59 mRNAs, 6 miRNAs (hsa-mir-133a-1, hsa-mir-1-2, hsa-mir-486, hsa-mir-135b, hsa-mir-196b, hsa-mir-193b) and 6 lncRNAs (C10orf91, C2orf48, SFTA1P, FLJ41941,PART1,TTTY14) were related with OS; and 52 mRNAs, 4 miRNAs (hsa-mir-133a-1, hsa-mir-135b, hsa-mir-196b, hsa-mir-193b) and 2 lncRNAs (PART1, TTTY14) were associated with RFS. A support vector machine (SVM) classifier containing 37 key hub genes was obtained. A ceRNA regulatory network containing 417 nodes and 696 edges was constructed. ECM-receptor interaction, cytokine-cytokine receptor interaction, focal adhesion, arachidonic acid metabolism, and p53 signaling pathway were significantly enriched in the network. CONCLUSION: These findings uncover the pathogenesis of OSCC and might provide potential therapeutic targets.

Identification of breast cancer recurrence risk factors based on functional pathways in tumor and normal tissues.

The recurrence of breast cancer (BC) is a serious therapeutic problem, and the risk factors for recurrence urgently need to be identified. In this study, we examined the functional pathways in tumor and normal tissues to more comprehensively identify biomarkers for the risk of BC recurrence. We collected tumor and normal tissue gene expression profiles of recurrent BC patients and non-recurrent BC patients from the TCGA database.We derived an expression interval (mean ± 1.96SD) based on non-recurrent patients rather than a single value, such as a mean or median. If the expression of a gene was significantly different from its normal expression interval, it was considered a differentially expressed gene. Eight pathways that significantly distinguished recurrent and non-recurrent BC patients were obtained based on 65% accuracy, and these pathways were all associated with the immune response and sensitivity to drugs. The genes in these eight pathways were also used to analyze survival, and the significance level reached 0.003 in an independent dataset (p = 0.02 in tumor and p = 0.03 in normal tissue). Our results reveal that the integration of tumor and normal tissue functional analyses can comprehensively enhance the understanding of BC prognosis.

IGSA: Individual Gene Sets Analysis, including Enrichment and Clustering.

Analysis of gene sets has been widely applied in various high-throughput biological studies. One weakness in the traditional methods is that they neglect the heterogeneity of genes expressions in samples which may lead to the omission of some specific and important gene sets. It is also difficult for them to reflect the severities of disease and provide expression profiles of gene sets for individuals. We developed an application software called IGSA that leverages a powerful analytical capacity in gene sets enrichment and samples clustering. IGSA calculates gene sets expression scores for each sample and takes an accumulating clustering strategy to let the samples gather into the set according to the progress of disease from mild to severe. We focus on gastric, pancreatic and ovarian cancer data sets for the performance of IGSA. We also compared the results of IGSA in KEGG pathways enrichment with David, GSEA, SPIA, ssGSEA and analyzed the results of IGSA clustering and different similarity measurement methods. Notably, IGSA is proved to be more sensitive and specific in finding significant pathways, and can indicate related changes in pathways with the severity of disease. In addition, IGSA provides with significant gene sets profile for each sample.