Bioinformatics and cancer target discovery

Comprehensive bioinformatics analysis of human cytomegalovirus pathway genes in pan-cancer

BACKGROUND

Human cytomegalovirus (HCMV) is a herpesvirus that can infect various cell types and modulate host gene expression and immune response. It has been associated with the pathogenesis of various cancers, but its molecular mechanisms remain elusive.

METHODS

We comprehensively analyzed the expression of HCMV pathway genes across 26 cancer types using the Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) databases. We also used bioinformatics tools to study immune invasion and tumor microenvironment in pan-cancer. Cox regression and machine learning were used to analyze prognostic genes and their relationship with drug sensitivity.

RESULTS

We found that HCMV pathway genes are widely expressed in various cancers. Immune infiltration and the tumor microenvironment revealed that HCMV is involved in complex immune processes. We obtained prognostic genes for 25 cancers and significantly found 23 key genes in the HCMV pathway, which are significantly enriched in cellular chemotaxis and synaptic function and may be involved in disease progression. Notably, CaM family genes were up-regulated and AC family genes were down-regulated in most tumors. These hub genes correlate with sensitivity or resistance to various drugs, suggesting their potential as therapeutic targets.

CONCLUSIONS

Our study has revealed the role of the HCMV pathway in various cancers and provided insights into its molecular mechanism and therapeutic significance. It is worth noting that the key genes of the HCMV pathway may open up new doors for cancer prevention and treatment.

Key genes and immune infiltration in chronic spontaneous urticaria: a study of bioinformatics and systems biology

Background

Chronic spontaneous urticaria (CSU) is defined by the spontaneous occurrence of wheals and/or angioedema for >6 weeks. The pathogenesis involves skin mast cells, but the complex causes of their activation remain to be characterized in detail.

Objectives

To explore disease-driving genes and biological pathways in CSU.

Methods

Two microarray data sets, e.g., GSE57178 and GSE72540, with mRNA information of skin from CSU patients, were downloaded from the Gene Expression Omnibus (GEO) database. An integrated bioinformatics pipeline including identification of differentially expressed genes (DEGs), functional enrichment analysis, protein-protein interaction (PPI) network analysis, co-expression and drug prediction analysis, and immune and stromal cells deconvolution analyses were applied to identify hub genes and key drivers of CSU pathogenesis.

Results

In total, we identified 92 up-regulated and 7 down-regulated genes in CSU lesions. These were significantly enriched in CSU-related pathways such as TNF, NF-κB, and JAK-STAT signaling. Based on PPI network modeling, four genes, i.e., IL-6, TLR-4, ICAM-1, and PTGS-2, were computationally identified as key pathogenic players in CSU. Immune infiltration analyses indicated that dendritic cells, Th2 cells, mast cells, megakaryocyte-erythroid progenitor, preadipocytes, and M1 macrophages were increased in lesional CSU skin.

Conclusion

Our results offer new insights on the pathogenesis of CSU and suggest that TNF, NF-κB, JAK-STAT, IL-6, TLR-4, ICAM-1, and PTGS-2 may be candidate targets for novel CSU treatments.

Novel tumor necrosis factor-related long non-coding RNAs signature for risk stratification and prognosis in glioblastoma

Background

Tumor necrosis factor (TNF) is an inflammatory cytokine that can coordinate tissue homeostasis by co-regulating the production of cytokines, cell survival, or death. It widely expresses in various tumor tissues and correlates with the malignant clinical features of patients. As an important inflammatory factor, the role of TNFα is involved in all steps of tumorigenesis and development, including cell transformation, survival, proliferation, invasion and metastasis. Recent research has showed that long non-coding RNAs (lncRNAs), defined as RNA transcripts >200 nucleotides that do not encode a protein, influence numerous cellular processes. However, little is known about the genomic profile of TNF pathway related-lncRNAs in GBM. This study investigated the molecular mechanism of TNF related-lncRNAs and their immune characteristics in glioblastoma multiforme (GBM) patients.

Methods

To identify TNF associations in GBM patients, we performed bioinformatics analysis of public databases - The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). The ConsensusClusterPlus, CIBERSORT, Estimate, GSVA and TIDE and first-order bias correlation and so on approaches were conducted to comprehensively characterize and compare differences among TNF-related subtypes.

Results

Based on the comprehensive analysis of TNF-related lncRNAs expression profiles, we constructed six TNF-related lncRNAs (C1RL-AS1, LINC00968, MIR155HG, CPB2-AS1, LINC00906, and WDR11-AS1) risk signature to determine the role of TNF-related lncRNAs in GBM. This signature could divide GBM patients into subtypes with distinct clinical and immune characteristics and prognoses. We identified three molecular subtypes (C1, C2, and C3), with C2 showing the best prognosis; otherwise, C3 showing the worst prognosis. Moreover, we assessed the prognostic value, immune infiltration, immune checkpoints, chemokines cytokines and enrichment analysis of this signature in GBM. The TNF-related lncRNA signature was tightly associated with the regulation of tumor immune therapy and could serve as an independent prognostic biomarker in GBM.

Conclusion

This analysis provides a comprehensive understanding of the role of TNF-related characters, which may improve the clinical outcome of GBM patients.

Systemic analysis identifying PVT1 / DUSP13 axis for microvascular invasion in hepatocellular carcinoma

BACKGROUND

Microvascular invasion (MVI) is an independent detrimental risk factor for tumor recurrence and poor survival in hepatocellular carcinoma (HCC). Competitive endogenous RNA (ceRNA) networks play a pivotal role in the modulation of carcinogenesis and progression among diverse tumor types. However, whether the ceRNA mechanisms are engaged in promoting the MVI process in patients with HCC remains unknown.

METHODS

A ceRNA regulatory network was constructed based on RNA-seq data of patients with HCC from The Cancer Genome Atlas (TCGA) database. In total, 10 hub genes of the ceRNA network were identified using four algorithms: "MCC," "Degree," "Betweenness," and "Stress." Transcriptional expressions were verified by in situ hybridization using clinical samples. Interactions between ceRNA modules were validated by luciferase reporting assay. Logistic regression analysis, correlation analysis, enrichment analysis, promoter region analysis, methylation analysis, and immune infiltration analysis were performed to further investigate the molecular mechanisms and clinical transformation value.

RESULTS

The ceRNA regulatory network featuring a tumor invasion phenotype consisting of 3 long noncoding RNAs, 3 microRNAs, and 93 mRNAs was constructed using transcriptional data from the TCGA database. Systemic analysis and experimentally validation identified a ceRNA network (PVT1/miR-1258/DUSP13 axis) characterized by lipid regulatory potential, immune properties, and abnormal methylation states in patients with HCC and MVI. Meanwhile, 28 transcriptional factors were identified as potential promotors of PVT1 with 3 transcriptional factors MXD3, ZNF580, and KDM1A promising as therapeutic targets in patients with HCC and MVI. Furthermore, miR-1258 was an independent predictor for MVI in patients with HCC.

CONCLUSION

The PVT1/DUSP13 axis is significantly associated with MVI progression in HCC patients. This study provides new insight into mechanisms related to lipids, immune phenotypes, and abnormal epigenetics in oncology research.

Development and validation of a ferroptosis-related lncRNAs signature to predict prognosis and microenvironment for melanoma

Ferroptosis plays an important role in cancer. However, studies about ferroptosis-related lncRNAs (FRLs) in skin cutaneous melanoma (SKCM) are scarce. Moreover, the relationship between prognostic FRLs and tumor microenvironment (TME) in melanoma remains unclear. This study investigates the potential prognostic value of FRLs and their association with TME in SKCM. The RNA-sequencing data of SKCM were downloaded from The Cancer Genome Atlas (TCGA) database. Melanoma patients were randomly divided into training and testing groups in a 1:1 ratio. A signature composed of 19 FRLs was developed by the least absolute shrinkage and selection operator (LASSO) regression analysis to divide patients into a low-risk group with a better prognosis and a high-risk group with a poor prognosis. Multivariate Cox regression analysis suggested that the risk score was an independent prognostic factor. The Area Under Curve (AUC) value of the risk score reached 0.768 in the training group and 0.770 in the testing group. Subsequent analysis proved that immune-related signaling pathways were significantly enriched in the low-risk group. The tumor immune cell infiltration analysis demonstrated that melanoma in the high-risk group tended to be immunologically "cold". We identified a novel FRLs signature which could accurately predict the prognosis of patients with melanoma.

The cuproptosis-related signature associated with the tumor environment and prognosis of patients with glioma

Background

Copper ions are essential for cellular physiology. Cuproptosis is a novel method of copper-dependent cell death, and the cuproptosis-based signature for glioma remains less studied.

Methods

Several glioma datasets with clinicopathological information were collected from TCGA, GEO and CGGA. Robust Multichip Average (RMA) algorithm was used for background correction and normalization, cuproptosis-related genes (CRGs) were then collected. The TCGA-glioma cohort was clustered using ConsensusClusterPlus. Univariate Cox regression analysis and the Random Survival Forest model were performed on the differentially expressed genes to identify prognostic genes. The cuproptosis-signature was constructed by calculating CuproptosisScore using Multivariate Cox regression analysis. Differences in terms of genomic mutation, tumor microenvironment, and enrichment pathways were evaluated between high- or low-CuproptosisScore. Furthermore, drug response prediction was carried out utilizing pRRophetic.

Results

Two subclusters based on CRGs were identified. Patients in cluster2 had better clinical outcomes. The cuproptosis-signature was constructed based on CuproptosisScore. Patients with higher CuproptosisScore had higher WHO grades and worse prognosis, while patients with lower grades were more likely to develop IDH mutations or MGMT methylation. Univariate and Multivariate Cox regression analysis demonstrated CuproptosisScore was an independent prognostic factor. The accuracy of the signature in prognostic prediction was further confirmed in 11 external validation datasets. In groups with high-CuproptosisScore, PIK3CA, MUC16, NF1, TTN, TP53, PTEN, and EGFR showed high mutation frequency. IDH1, TP53, ATRX, CIC, and FUBP1 demonstrated high mutation frequency in low-CuproptosisScore group. The level of immune infiltration increased as CuproptosisScore increased. SubMap analysis revealed patients with high-CuproptosisScore may respond to anti-PD-1 therapy. The IC50 values of Bexarotene, Bicalutamide, Bortezomib, and Cytarabine were lower in the high-CuproptosisScore group than those in the low-CuproptosisScore group. Finally, the importance of IGFBP2 in TCGA-glioma cohort was confirmed.

Conclusion

The current study revealed the novel cuproptosis-based signature might help predict the prognosis, biological features, and appropriate treatment for patients with glioma.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

Expression and clinical implications of basic leucine zipper ATF-like transcription factor 2 in breast cancer

Background

Basic leucine zipper ATF-like transcription factor 2 (BATF2) has been reported to participate in the occurrence and development of some malignancies. Herein, we aimed to explore the expression pattern and clinical implications of BATF2 in breast cancer (BC).

Methods

We assessed the differences in BATF2 mRNA expression between cancerous and noncancerous tissues in BC using GEPIA and UALCAN data and in BATF2 protein expression pattern using Human Protein Atlas (HPA) data. BATF2 co-expression networks were analyzed in Coexpedia. The association between the differentially expressed BATF2 mRNA and BC prognosis was assessed using UALCAN, OSbrca, and GEPIA databases. In external validations, BATF2 protein expression in BC tissues was quantitated using a tissue microarray and immunohistochemistry (IHC) analysis, and BATF2 mRNA expression in serum and serum-derived exosomes of BC patients using real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR).

Results

No difference in the BATF2 mRNA expression level was found between cancerous and noncancerous tissues in BC based on databases. There were low-to-moderate levels of increases in BATF2 protein expressions in BC cases from the HPA cohort. BATF2 mRNA expression was negatively correlated with androgen receptor (AR) and positively correlated with BRCA2 DNA repair associated (BRCA2), marker of proliferation Ki-67 (Mki67), and tumor protein p53 (TP53) expressions. Generally, BATF2 mRNA exhibited a non-significant association with BC prognosis; yet the subgroup analyses showed that triple-negative breast cancer (TNBC) patients with high BATF2 mRNA expressions had a longer overall survival (OS). Our IHC analysis revealed a positive rate of BATF2 protein expression of 46.90%, mainly located in the nucleus of cancer cells in BC, and the OS of BC patients with high BATF2 protein expressions was prolonged. The positive rates of BATF2 mRNA expressions in the serum and exosomes were 45.00 and 41.67%, respectively. Besides, the AUCs of serum and exosomal BATF2 mRNA for BC diagnosis were 0.8929 and 0.8869, respectively.

Conclusions

BC patients exhibit low-to-moderate expressions in BATF2 mRNA expression levels in cancerous tissues. The high BATF2 protein expression can be a potential indicator of a better BC prognosis. Serum and exosomal BATF2 mRNA levels also serve as promising noninvasive biomarkers for BC diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08785-6.

Development and Validation of an Mesenchymal-Related Long Non-Coding RNA Prognostic Model in Glioma

Glioma is well known as the most aggressive and prevalent primary malignant tumor in the central nervous system. Molecular subtypes and prognosis biomarkers remain a promising research area of gliomas. Notably, the aberrant expression of mesenchymal (MES) subtype related long non-coding RNAs (lncRNAs) is significantly associated with the prognosis of glioma patients. In this study, MES-related genes were obtained from The Cancer Genome Atlas (TCGA) and the Ivy Glioblastoma Atlas Project (Ivy GAP) data sets of glioma, and MES-related lncRNAs were acquired by performing co-expression analysis of these genes. Next, Cox regression analysis was used to establish a prognostic model, that integrated ten MES-related lncRNAs. Glioma patients in TCGA were divided into high-risk and low-risk groups based on the median risk score; compared with the low-risk groups, patients in the high-risk group had shorter survival times. Additionally, we measured the specificity and sensitivity of our model with the ROC curve. Univariate and multivariate Cox analyses showed that the prognostic model was an independent prognostic factor for glioma. To verify the predictive power of these candidate lncRNAs, the corresponding RNA-seq data were downloaded from the Chinese Glioma Genome Atlas (CGGA), and similar results were obtained. Next, we performed the immune cell infiltration profile of patients between two risk groups, and gene set enrichment analysis (GSEA) was performed to detect functional annotation. Finally, the protective factors DGCR10 and HAR1B, and risk factor SNHG18 were selected for functional verification. Knockdown of DGCR10 and HAR1B promoted, whereas knockdown of SNHG18 inhibited the migration and invasion of gliomas. Collectively, we successfully constructed a prognostic model based on a ten MES-related lncRNAs signature, which provides a novel target for predicting the prognosis for glioma patients.

Autophagy-related circRNA evaluation reveals hsa_circ_0001747 as a potential favorable prognostic factor for biochemical recurrence in patients with prostate cancer

Prostate cancer (PCa) is a common high-incidence malignancy in men, some of whom develop biochemical recurrence (BCR) in the advanced stage. However, there are currently no accurate prognostic indicators of BCR in PCa. The aim of our study was to identify an autophagy-related circular RNA prognostic factor of BCR for patients with PCa. In this study, immunochemistry revealed that the classic autophagy marker MAP1LC3B was positively correlated with Gleason score. Least absolute shrinkage and selector operator regression were conducted to develop a novel prognostic model with tenfold cross-validation and an L1 penalty. Five autophagy-related circRNA signatures were included in the prognostic model. Patients with PCa were ultimately divided into high- and low-risk groups, based on the median risk score. Patients with PCa, who had a high risk score, were more likely to develop BCR in a shorter period of time. Univariate and multivariate Cox regression analyses demonstrated that the risk score was an independent variable for predicting BCR in PCa. In addition, a prognostic nomogram integrated with the risk score and numerous clinicopathological parameters was developed to accurately predict 3- and 5-year BCR of patients with PCa. Finally, the hsa_circ_0001747 signature was selected for further experimental verification in vitro and in vivo, which showed that downregulated hsa_circ_0001747 might facilitate PCa via augmenting autophagy. Our findings indicate that the autophagy-related circRNA signature hsa_circ_0001747 may serve as a promising indicator for BCR prediction in patients with PCa.

Anthelmintic activity of a nanoformulation based on thiophenes identified in Tagetes patula L. (Asteraceae) against the small ruminant nematode Haemonchus contortus

The synthesis of thiophenic compounds, previously identified in Tagetes patula, revealed that 4-(5'-(hydroxymethyl)-[2,2'-bithiophene]-5-yl)but-3-yn-1-ol), or simply Thio1, has a pronounced in vitro anthelmintic effect against Haemonchus contortus, showing 100% efficacy in the egg hatch and larval development tests presenting EC₅₀ = 0.1731 mg.mL^-1 and EC₅₀ = 0.3243 mg.mL^-1, respectively. So, this compound was selected to preparation of a nanostructured formulation to be orally administered to Santa Inês sheep. In general, from the fecal egg count reduction test (FECRT), it was observed that the product kept the parasitic load in the digestive tract of the hosts stable, with eggs per gram of faeces (EPG) counts having a mean value < 3,000 (EPG_mean = 2167.1, efficacy = 36,45%), thus protecting the animals from health risks caused by a massive nematode infestation. To better understand the mode of action of this thiophene derivative, in silico molecular modeling studies were carried out with the glutamate-activated chloride channel (GluCl), a well-known molecular target of anthelmintic compounds. Based on the affinity score (GlideScore = -5.7 kcal.mol^-1) and the proposed binding mode, Thio1 could be classified as a potential GluCl ligand, justifying the promising results observed in the anthelmintic assays.

A Ferroptosis-Related lncRNAs Signature Predicts Prognosis and Immune Microenvironment for Breast Cancer

Background: Ferroptosis, a regulated cell death which is driven by the iron-dependent peroxidation of lipids, plays an important role in cancer. However, studies about ferroptosis-related Long non-coding RNAs (lncRNAs) in breast cancer (BC) are limited. Besides, the prognostic role of ferroptosis-related lncRNAs and their relationship to immune microenvironment in breast cancer remain unclear. This study aimed to explore the potential prognostic value of ferroptosis-related lncRNAs and their relationship to immune microenvironment in breast cancer.

Methods: RNA-sequencing data of female breast cancer patients were downloaded from TCGA database. 937 patients were randomly separated into training or validation cohort in 2:1 ratio. Ferroptosis-related lncRNAs were screened by Pearson correlation analysis with 239 reported ferroptosis-related genes. A ferroptosis-related lncRNAs signature was constructed with univariate and multivariate Cox regression analyses in the training cohort, and its prognostic value was further tested in the validation cohort.

Results: An 8-ferroptosis-related-lncRNAs signature was developed by multivariate Cox regression analysis to divide patients into two risk groups. Patients in the high-risk group had worse prognosis than patients in the low-risk group. Multivariate Cox regression analysis showed the risk score was an independent prognostic indicator. Receiver operating characteristic curve (ROC) analysis proved the predictive accuracy of the signature. The area under time-dependent ROC curve (AUC) reached 0.853 at 1 year, 0.802 at 2 years, 0.740 at 5 years in the training cohort and 0.791 at 1 year, 0.778 at 2 years, 0.722 at 5 years in the validation cohort. Further analysis demonstrated that immune-related pathways were significantly enriched in the high-risk group. Analysis of the immune cell infiltration landscape showed that breast cancer in the high-risk group tended be immunologically “cold”.

Conclusion: We identified a novel ferroptosis-related lncRNA signature which could precisely predict the prognosis of breast cancer patients. Ferroptosis-related lncRNAs may have a potential role in the process of anti-tumor immunity and serve as therapeutic targets for breast cancer.

An autophagy-related long non-coding RNA prognostic signature accurately predicts survival outcomes in bladder urothelial carcinoma patients

In this study, we analyzed the prediction accuracy of an autophagy-related long non-coding RNA (lncRNA) prognostic signature using bladder urothelial carcinoma (BLCA) patient data from The Cancer Genome Atlas (TCGA) database. Univariate and multivariate Cox regression analyses showed significant correlations between five autophagy-related lncRNAs, LINC02178, AC108449.2, Z83843.1, FAM13A-AS1 and USP30-AS1, and overall survival (OS) among BCLA patients. The risk scores based on the autophagy-related lncRNA prognostic signature accurately distinguished high- and low-risk BCLA patients that were stratified according to age; gender; grade; and AJCC, T, and N stages. The autophagy-related lncRNA signature was an independent prognostic predictor with an AUC value of 0.710. The clinical nomogram with the autophagy-related lncRNA prognostic signature showed a high concordance index of 0.73 and accurately predicted 1-, 3-, and 5-year survival times among BCLA patients in the high- and low-risk groups. The lncRNA-mRNA co-expression network contained 77 lncRNA-mRNA links among 5 lncRNAs and 49 related mRNAs. Gene set enrichment analysis showed that cancer- and autophagy-related pathways were significantly enriched in the high-risk group, and immunoregulatory pathways were enriched in the low-risk group. These findings demonstrate that an autophagy-related lncRNA signature accurately predicts the prognosis of BCLA patients.

Large-scale in-silico identification of a tumor-specific antigen pool for targeted immunotherapy in triple-negative breast cancer

Since the advent of cetuximab, clinical cancer treatment has evolved from the standard, relatively nonspecific chemo- and radiotherapy with significant cytotoxic side effects towards immunotherapeutic approaches with selective, target-mechanism-based effects. Antibody therapies as the most successful form of cancer immunotherapy led to approved treatments for specific cancer types with increased patient survival. Thus, the identification of tumor antigens with high immunogenicity is in central focus now. In this study, we applied computational methods to comprehensively discover overexpressed molecular targets with high therapeutic relevance for clinical, immunotherapeutic cancer treatment in triple-negative breast cancer (TNBC). By actively modeling potential negative side effects utilizing expression data of 29 different, normal human tissues, we were able to develop a highly-specific coverage of TNBC patients with RNA targets. We identified here more than 400 potential tumor-specific antigens suitable for targeted therapy, including several already identified as potential targets for TNBC and other solid tumors. A specific cocktail of MAGEB4, CT83, TLX3, ACTL8, PRDM13 achieved almost 94% patient coverage in TNBC. Overall, these results show that our approach can identify and prioritize TNBC targets suitable for targeted therapy. Therefore, our method has the potential to lead to new and more effective immunotherapeutic cancer treatment.

In silico Leishmania proteome mining applied to identify drug target potential to be used to treat against visceral and tegumentary leishmaniasis

New therapeutic strategies against leishmaniasis are desirable, since the treatment against disease presents problems, such as the toxicity, high cost and/or parasite resistance. As consequence, new antileishmanial compounds are necessary to be identified, as presenting high activity against Leishmania, but low toxicity in mammalian hosts. In the present study, a Leishmania proteome mining strategy was developed, in order to select new drug targets with low homology to human proteins, but that are considered relevant for the parasite' survival. Results showed a hypothetical protein, which was functionally annotated as a glucosidase-like protein, as presenting such characteristics. This protein was associated with the metabolic network of the N-Glycan biosynthesis pathway in Leishmania, and two specific inhibitors - acarbose and miglitol - were predicted to be potential targets against it. In this context, miglitol [1-(2-Hydroxyethyl)-2-(hydroxymethyl)piperidine-3,4,5-triol] was tested against stationary promastigotes and axenic amastigotes of the Leishmania amazonensis and L. infantum species, and results showed high values of antileishmanial inhibition against both parasite species. Miglitol showed also efficacy in the treatment of Leishmania-infected macrophages; thus denoting its potential use as an antileishmanial candidate. In conclusion, this work presents a new drug target identified by a proteome mining strategy associated with bioinformatics tools, and suggested its use as a possible candidate to be applied in the treatment against disease.

GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses

Tremendous amount of RNA sequencing data have been produced by large consortium projects such as TCGA and GTEx, creating new opportunities for data mining and deeper understanding of gene functions. While certain existing web servers are valuable and widely used, many expression analysis functions needed by experimental biologists are still not adequately addressed by these tools. We introduce GEPIA (Gene Expression Profiling Interactive Analysis), a web-based tool to deliver fast and customizable functionalities based on TCGA and GTEx data. GEPIA provides key interactive and customizable functions including differential expression analysis, profiling plotting, correlation analysis, patient survival analysis, similar gene detection and dimensionality reduction analysis. The comprehensive expression analyses with simple clicking through GEPIA greatly facilitate data mining in wide research areas, scientific discussion and the therapeutic discovery process. GEPIA fills in the gap between cancer genomics big data and the delivery of integrated information to end users, thus helping unleash the value of the current data resources. GEPIA is available at http://gepia.cancer-pku.cn/.

Computational selection of antibody-drug conjugate targets for breast cancer

The selection of therapeutic targets is a critical aspect of antibody-drug conjugate research and development. In this study, we applied computational methods to select candidate targets overexpressed in three major breast cancer subtypes as compared with a range of vital organs and tissues. Microarray data corresponding to over 8,000 tissue samples were collected from the public domain. Breast cancer samples were classified into molecular subtypes using an iterative ensemble approach combining six classification algorithms and three feature selection techniques, including a novel kernel density-based method. This feature selection method was used in conjunction with differential expression and subcellular localization information to assemble a primary list of targets. A total of 50 cell membrane targets were identified, including one target for which an antibody-drug conjugate is in clinical use, and six targets for which antibody-drug conjugates are in clinical trials for the treatment of breast cancer and other solid tumors. In addition, 50 extracellular proteins were identified as potential targets for non-internalizing strategies and alternative modalities. Candidate targets linked with the epithelial-to-mesenchymal transition were identified by analyzing differential gene expression in epithelial and mesenchymal tumor-derived cell lines. Overall, these results show that mining human gene expression data has the power to select and prioritize breast cancer antibody-drug conjugate targets, and the potential to lead to new and more effective cancer therapeutics.

Identifying new targets in leukemogenesis using computational approaches

There is a need to identify novel targets in Acute Lymphoblastic Leukemia (ALL), a hematopoietic cancer affecting children, to improve our understanding of disease biology and that can be used for developing new therapeutics. Hence, the aim of our study was to find new genes as targets using in silico studies; for this we retrieved the top 10% overexpressed genes from Oncomine public domain microarray expression database; 530 overexpressed genes were short-listed from Oncomine database. Then, using prioritization tools such as ENDEAVOUR, DIR and TOPPGene online tools, we found fifty-four genes common to the three prioritization tools which formed our candidate leukemogenic genes for this study. As per the protocol we selected thirty training genes from PubMed. The prioritized and training genes were then used to construct STRING functional association network, which was further analyzed using cytoHubba hub analysis tool to investigate new genes which could form drug targets in leukemia. Analysis of the STRING protein network built from these prioritized and training genes led to identification of two hub genes, SMAD2 and CDK9, which were not implicated in leukemogenesis earlier. Filtering out from several hundred genes in the network we also found MEN1, HDAC1 and LCK genes, which re-emphasized the important role of these genes in leukemogenesis. This is the first report on these five additional signature genes in leukemogenesis. We propose these as new targets for developing novel therapeutics and also as biomarkers in leukemogenesis, which could be important for prognosis and diagnosis.

In-silico prediction of drug targets, biological activities, signal pathways and regulating networks of dioscin based on bioinformatics

Background

Inverse docking technology has been a trend of drug discovery, and bioinformatics approaches have been used to predict target proteins, biological activities, signal pathways and molecular regulating networks affected by drugs for further pharmacodynamic and mechanism studies.

Methods

In the present paper, inverse docking technology was applied to screen potential targets from potential drug target database (PDTD). Then, the corresponding gene information of the obtained drug-targets was applied to predict the related biological activities, signal pathways and processes networks of the compound by using MetaCore platform. After that, some most relevant regulating networks were considered, which included the nodes and relevant pathways of dioscin.

Results

71 potential targets of dioscin from humans, 7 from rats and 8 from mice were screened, and the prediction results showed that the most likely targets of dioscin were cyclin A2, calmodulin, hemoglobin subunit beta, DNA topoisomerase I, DNA polymerase lambda, nitric oxide synthase and UDP-N-acetylhexosamine pyrophosphorylase, etc. Many diseases including experimental autoimmune encephalomyelitis of human, temporal lobe epilepsy of rat and ankylosing spondylitis of mouse, may be inhibited by dioscin through regulating immune response alternative complement pathway, G-protein signaling RhoB regulation pathway and immune response antiviral actions of interferons, etc. The most relevant networks (5 from human, 3 from rat and 5 from mouse) indicated that dioscin may be a TOP1 inhibitor, which can treat cancer though the cell cycle– transition and termination of DNA replication pathway. Dioscin can down regulate EGFR and EGF to inhibit cancer, and also has anti-inflammation activity by regulating JNK signaling pathway.

Conclusions

The predictions of the possible targets, biological activities, signal pathways and relevant regulating networks of dioscin provide valuable information to guide further investigation of dioscin on pharmacodynamics and molecular mechanisms, which also suggests a practical and effective method for studies on the mechanism of other chemicals.

Biology and bioinformatics of myeloma cell

Multiple myeloma (MM) is a plasma cell disorder that occurs in about 10% of all hematologic cancers. The majority of patients (99%) are over 50 years of age when diagnosed. In the bone marrow (BM), stromal and hematopoietic stem cells (HSCs) are responsible for the production of blood cells. Therefore any destruction or/and changes within the BM undesirably impacts a wide range of hematopoiesis, causing diseases and influencing patient survival. In order to establish an effective therapeutic strategy, recognition of the biology and evaluation of bioinformatics models for myeloma cells are necessary to assist in determining suitable methods to cure or prevent disease complications in patients. This review presents the evaluation of molecular and cellular aspects of MM such as genetic translocation, genetic analysis, cell surface marker, transcription factors, and chemokine signaling pathways. It also briefly reviews some of the mechanisms involved in MM in order to develop a better understanding for use in future studies.

Target discovery from data mining approaches

Data mining of available biomedical data and information has greatly boosted target discovery in the 'omics' era. Target discovery is the key step in the biomarker and drug discovery pipeline to diagnose and fight human diseases. In biomedical science, the 'target' is a broad concept ranging from molecular entities (such as genes, proteins and miRNAs) to biological phenomena (such as molecular functions, pathways and phenotypes). Within the context of biomedical science, data mining refers to a bioinformatics approach that combines biological concepts with computer tools or statistical methods that are mainly used to discover, select and prioritize targets. In response to the huge demand of data mining for target discovery in the 'omics' era, this review explicates various data mining approaches and their applications to target discovery with emphasis on text and microarray data analysis. Two emerging data mining approaches, chemogenomic data mining and proteomic data mining, are briefly introduced. Also discussed are the limitations of various data mining approaches found in the level of database integration, the quality of data annotation, sample heterogeneity and the performance of analytical and mining tools. Tentative strategies of integrating different data sources for target discovery, such as integrated text mining with high-throughput data analysis and integrated mining with pathway databases, are introduced.

European Organisation for Research and Treatment of Cancer (EORTC) Pathobiology Group standard operating procedure for the preparation of human tumour tissue extracts suited for the quantitative analysis of tissue-associated biomarkers

With the new concept of 'individualized treatment and targeted therapies', tumour tissue-associated biomarkers have been given a new role in selection of cancer patients for treatment and in cancer patient management. Tumour biomarkers can give support to cancer patient stratification and risk assessment, treatment response identification, or to identifying those patients who are expected to respond to certain anticancer drugs. As the field of tumour-associated biomarkers has expanded rapidly over the last years, it has become increasingly apparent that a strong need exists to establish guidelines on how to easily disintegrate the tumour tissue for assessment of the presence of tumour tissue-associated biomarkers. Several mechanical tissue (cell) disruption techniques exist, ranging from bead mill homogenisation and freeze-fracturing through to blade or pestle-type homogenisation, to grinding and ultrasonics. Still, only a few directives have been given on how fresh-frozen tumour tissues should be processed for the extraction and determination of tumour biomarkers. The PathoBiology Group of the European Organisation for Research and Treatment of Cancer therefore has devised a standard operating procedure for the standardised preparation of human tumour tissue extracts which is designed for the quantitative analysis of tumour tissue-associated biomarkers. The easy to follow technical steps involved require 50-300 mg of deep-frozen cancer tissue placed into small size (1.2 ml) cryogenic tubes. These are placed into the shaking flask of a Mikro-Dismembrator S machine (bead mill) to pulverise the tumour tissue in the capped tubes in the deep-frozen state by use of a stainless steel ball, all within 30 s of exposure. RNA is isolated from the pulverised tissue following standard procedures. Proteins are extracted from the still frozen pulverised tissue by addition of Tris-buffered saline to obtain the cytosol fraction of the tumour or by the Tris buffer supplemented with the non-ionic detergent Triton X-100, and, after high-speed centrifugation, are found in the tissue supernatant. The resulting tissue cell debris sediment is a rich source of genomic DNA.

Therapeutic targets: progress of their exploration and investigation of their characteristics

Modern drug discovery is primarily based on the search and subsequent testing of drug candidates acting on a preselected therapeutic target. Progress in genomics, protein structure, proteomics, and disease mechanisms has led to a growing interest in and effort for finding new targets and more effective exploration of existing targets. The number of reported targets of marketed and investigational drugs has significantly increased in the past 8 years. There are 1535 targets collected in the therapeutic target database compared with approximately 500 targets reported in a 1996 review. Knowledge of these targets is helpful for molecular dissection of the mechanism of action of drugs and for predicting features that guide new drug design and the search for new targets. This article summarizes the progress of target exploration and investigates the characteristics of the currently explored targets to analyze their sequence, structure, family representation, pathway association, tissue distribution, and genome location features for finding clues useful for searching for new targets. Possible "rules" to guide the search for druggable proteins and the feasibility of using a statistical learning method for predicting druggable proteins directly from their sequences are discussed.

Serological identification and bioinformatics analysis of immunogenic antigens in multiple myeloma

Identifying appropriate tumor antigens is critical to the development of successful specific cancer immunotherapy. Serological analysis of tumor antigens by a recombinant cDNA expression library (SEREX) allows the systematic cloning of tumor antigens recognized by the spontaneous autoantibody repertoire of cancer patients. We applied SEREX to the cDNA expression library of cell line HMy2, which led to the isolation of six known characterized genes and 12 novel genes. Known genes, including ring finger protein 167, KLF10, TPT1, p02 protein, cDNA FLJ46859 fis, and DNMT1, were related to the development of different tumors. Bioinformatics was performed to predict 12 novel MMSA (multiple myeloma special antigen) genes. The prediction of tumor antigens provides potential targets for the immunotherapy of patients with multiple myeloma (MM) and help in the understanding of carcinogenesis. Crude lysate ELISA methodology indicated that the optical density value of MMSA-3 and MMSA-7 were significantly higher in MM patients than in healthy donors. Furthermore, SYBR Green real-time PCR showed that MMSA-1 presented with a high number of copy messages in MM. In summary, the antigens identified in this study may be potential candidates for diagnosis and targets for immunotherapy in MM.

Progress and problems in the exploration of therapeutic targets

Drugs exert their therapeutic effect by binding and regulating the activity of a particular protein or nucleic acid target. A large number of targets have been explored for drug discovery. Continuous effort has been directed at the search for new targets and more-extensive exploration of existing targets. Knowledge of these targets facilitates the understanding of molecular mechanisms of drugs and the effort required for drug discovery and target searches. Areas of progress, current focuses of research and development and the difficulties in target exploration are reviewed. The characteristics of the currently explored targets and their correlation to the level of difficulty for target exploration are analyzed. From these characteristics, simple rules can be derived for estimating the difficulty level of target exploration. The feasibility of predicting druggable proteins by using simple rules and sequence-derived physicochemical properties is also discussed.

Oncology exploration: charting cancer medicinal chemistry space

Approaches for the experimental determination of protein-ligand molecular interactions are reliant on the quality of the compounds being tested. The application of large, randomly designed combinatorial libraries has given way to the creation of more-focused 'drug-like' libraries. Prior to synthesis, we wish to screen the potential compounds to remove undesired chemical moieties and to be within a required range of physiochemical properties. We have used a principal-component analysis (PCA) computational approach to analyze the 3D descriptor space of active and non-active (hit-like) cancer medicinal chemistry compounds. We define hit-like those molecules passing the unmodified OpenEye FILTER program. Our analysis indicates that these compounds occupy quite different regions in space. Cancer-active compounds exist in a much greater volume of space than generic hit-like space and most of them fail the commonly applied filters for orally bioavailable drugs. This is of great significance when designing orally bioavailable cancer target drugs.