Integration of multiple networks and pathways identifies cancer driver genes in pan-cancer analysis

A method to improve the prediction performance of cancer-gene association by screening negative training samples through gene network data

This work focuses on data sampling in cancer-gene association prediction. Currently, researchers are using machine learning methods to predict genes that are more likely to produce cancer-causing mutations. To improve the performance of machine learning models, methods have been proposed, one of which is to improve the quality of the training data. Existing methods focus mainly on positive data, i.e. cancer driver genes, for screening selection. This paper proposes a low-cancer-related gene screening method based on gene network and graph theory algorithms to improve the negative samples selection. Genetic data with low cancer correlation is used as negative training samples. After experimental verification, using the negative samples screened by this method to train the cancer gene classification model can improve prediction performance. The biggest advantage of this method is that it can be easily combined with other methods that focus on enhancing the quality of positive training samples. It has been demonstrated that significant improvement is achieved by combining this method with three state-of-the-arts cancer gene prediction methods.

Children's Oncology Group's 2023 blueprint for research: Germ cell tumors

Extracranial germ cell tumors (GCT) are a biologically diverse group of tumors occurring in children, adolescents, and young adults. The majority of patients have excellent outcomes, but treatment-related toxicities impact their quality of survivorship. A subset of patients succumbs to the disease. Current unmet needs include clarifying which patients can be safely observed after initial surgical resection, refinement of risk stratification to reduce chemotherapy burden in patients with standard-risk disease, and intensify therapy for patients with poor-risk disease. Furthermore, enhancing strategies for detection of minimal residual disease and early detection of relapse, particularly in serum tumor marker-negative histologies, is critical. Improving the understanding of the developmental and molecular origins of GCTs may facilitate discovery of novel targets. Future efforts should be directed toward assessing novel therapies in a biology-driven, biomarker-defined, histology-specific, risk-stratified patient population. Fragmentation of care between subspecialists restricts the unified study of these rare tumors. It is imperative that trials be conducted in collaboration with national and international cooperative groups, with harmonized data and biospecimen collection. Key priorities for the Children's Oncology Group (COG) GCT Committee include (a) better understanding the biology of GCTs, with a focus on molecular targets and mechanisms of treatment resistance; (b) strategic development of pediatric and young adult clinical trials; (c) understanding late effects of therapy and identifying individuals most at risk; and (d) prioritizing diversity, equity, and inclusion to reduce cancer health disparities and studying the impacts of social determinants of health on outcomes.

Individualized discovery of rare cancer drivers in global network context

Late advances in genome sequencing expanded the space of known cancer driver genes several-fold. However, most of this surge was based on computational analysis of somatic mutation frequencies and/or their impact on the protein function. On the contrary, experimental research necessarily accounted for functional context of mutations interacting with other genes and conferring cancer phenotypes. Eventually, just such results become ‘hard currency’ of cancer biology. The new method, NEAdriver employs knowledge accumulated thus far in the form of global interaction network and functionally annotated pathways in order to recover known and predict novel driver genes. The driver discovery was individualized by accounting for mutations’ co-occurrence in each tumour genome – as an alternative to summarizing information over the whole cancer patient cohorts. For each somatic genome change, probabilistic estimates from two lanes of network analysis were combined into joint likelihoods of being a driver. Thus, ability to detect previously unnoticed candidate driver events emerged from combining individual genomic context with network perspective. The procedure was applied to 10 largest cancer cohorts followed by evaluating error rates against previous cancer gene sets. The discovered driver combinations were shown to be informative on cancer outcome. This revealed driver genes with individually sparse mutation patterns that would not be detectable by other computational methods and related to cancer biology domains poorly covered by previous analyses. In particular, recurrent mutations of collagen, laminin, and integrin genes were observed in the adenocarcinoma and glioblastoma cancers. Considering constellation patterns of candidate drivers in individual cancer genomes opens a novel avenue for personalized cancer medicine.

Immunotherapeutic Value of MAP1LC3C and Its Candidate FDA-Approved Drugs Identified by Pan-Cancer Analysis, Virtual Screening and Sensitivity Analysis

Background: The autophagy pathway within the tumour microenvironment can be regulated to inhibit or promote tumour development. In the fight against tumour growth, immunotherapy induces an anti-tumour immune response, whereas autophagy modulates this immune response. A key protein in the autophagy pathway, microtubule-associated protein 1 light chain 3 (MAP1LC3), has recently become a hotspot for tumour research. As a relatively novel member, the function of MAP1LC3C in tumours still need to be investigated. Therefore, the goal of this study was to look into the possible link between MAP1LC3C and immunotherapy for 33 kinds of human malignancies by using pan-cancer analysis. Methods: High-throughput sequencing data from The Cancer Genome Atlas, Genotype-Tissue Expression Project and Cancer Cell Line Encyclopedia databases, combined with clinical data, were used to analyze the expression of MAP1LC3C in 33 types of cancer, as well as patient prognosis and neoplasm staging. Activity scores were calculated using ssGSEA to assess the MAP1LC3C activity in pan-cancer. Associations between MAP1LC3C and the tumour microenvironment, including immune cell infiltration and immunomodulators, were analyzed. Moreover, tumour tissue ImmuneScores and StromalScores were analyzed using the ESTIMATE algorithm. Additionally, associations between MAP1LC3C and tumour mutational burden/microsatellite instability, were investigated. Finally, based on the expression and structure of MAP1LC3C, the United States Food and Drug Administration (FDA)-approved drugs, were screened by virtual screening, molecular docking and NCI-60 drug sensitivity analysis. Results: Our study found that MAP1LC3C was differentially expressed in tumour and normal tissues in 23 of 33 human cancer types, among which MAP1LC3C had prognostic effects in 12 cancer types, and MAP1LC3C expression was significantly correlated with tumour stage in four cancer types. In addition, MAP1LC3C activity in 14 cancer types was consistent with changes in transcription levels. Moreover, MAP1LC3C strongly correlated with immune infiltration, immune modulators and immune markers. Finally, a number of FDA-approved drugs were identified via virtual screening and drug sensitivity analysis. Conclusion: Our study investigated the prognostic and immunotherapeutic value of MAP1LC3C in 33 types of cancer, and several FDA-approved drugs were identified to be highly related to MAP1LC3C and can be potential cancer therapeutic candidates.

Combining bioinformatics, network pharmacology and artificial intelligence to predict the mechanism of celastrol in the treatment of type 2 diabetes

BACKGROUND

Type 2 diabetes (T2D) is a common chronic disease with many serious complications. Celastrol can prevent and treat type 2 diabetes by reversing insulin resistance in a number of ways. However, the specific mechanisms by which celastrol prevents and treats T2D are not well understood. The aim of this study was to explore the key gene targets and potential signaling pathway mechanisms of celastrol for the treatment of T2D.

METHODS

GSE184050 was downloaded from the Gene Expression Omnibus online database. Blood samples from patients and healthy individuals with T2D were analyzed to identify differentially expressed genes (DEGs), and a protein-protein interaction network (PPI) was constructed. Key gene analysis of DEGs was performed using the MCODE plugin in Cystoscope as well as the Hubba plugin, and intersections were taken to obtain hub genes, which were displayed using a Venn diagram. Enrichment analysis was then performed via the ClueGo plugin in Cytoscape and validated using Gene Set Enrichment Analysis. The therapeutic targets of celastrol were then analyzed by pharmacophore network pharmacology, intersected to identify the therapeutic targets of celastrol, enriched for all targets, and intersected to obtain the signaling pathways for celastrol treatment. The protein structures of the therapeutic targets were predicted using the artificial intelligence AlphaFold2. Finally, molecular docking was used to verify whether celastrol could be successfully docked to the predicted targets.

RESULTS

618 DEGs were obtained, and 9 hub genes for T2D were identified by the MCODE and Hubba plug-ins, including ADAMTS15, ADAMTS7, ADAMTSL1, SEMA5B, ADAMTS8, THBS2, HBB, HBD and HBG2. The DEG-enriched signaling pathways mainly included the ferroptosis and TGF-beta signaling pathways. A total of 228 target genes were annotated by pharmacophore target analysis, and the therapeutic targets were identified, including S100A11, RBP3, HBB, BMP7 and IQUB, and 9 therapeutic signaling pathways were obtained by an intersectional set. The protein structures of the therapeutic targets were successfully predicted by AlphaFold2, and docking was validated using molecular docking.

CONCLUSION

Celastrol may prevent and treat T2D through key target genes, such as HBB, as well as signaling pathways, such as the TGF-beta signaling pathway and type II diabetes mellitus.

Application of Immune Infiltration Signature and Machine Learning Model in the Differential Diagnosis and Prognosis of Bone-Related Malignancies

Bone-related malignancies, such as osteosarcoma, Ewing's sarcoma, multiple myeloma, and cancer bone metastases have similar histological context, but they are distinct in origin and biological behavior. We hypothesize that a distinct immune infiltrative microenvironment exists in these four most common malignant bone-associated tumors and can be used for tumor diagnosis and patient prognosis. After sample cleaning, data integration, and batch effect removal, we used 22 publicly available datasets to draw out the tumor immune microenvironment using the ssGSEA algorithm. The diagnostic model was developed using the random forest. Further statistical analysis of the immune microenvironment and clinical data of patients with osteosarcoma and Ewing's sarcoma was carried out. The results suggested significant differences in the microenvironment of bone-related tumors, and the diagnostic accuracy of the model was higher than 97%. Also, high infiltration of multiple immune cells in Ewing's sarcoma was suggestive of poor patient prognosis. Meanwhile, increased infiltration of macrophages and B cells suggested a better prognosis for patients with osteosarcoma, and effector memory CD8 T cells and type 2 T helper cells correlated with patients' chemotherapy responsiveness and tumor metastasis. Our study revealed that the random forest diagnostic model based on immune infiltration can accurately perform the differential diagnosis of bone-related malignancies. The immune microenvironment of osteosarcoma and Ewing's sarcoma has an important impact on patient prognosis. Suppressing the highly inflammatory environment of Ewing's sarcoma and promoting macrophage and B cell infiltration may have good potential to be a novel adjuvant treatment option for osteosarcoma and Ewing's sarcoma.

Four Novel Prognostic Genes Related to Prostate Cancer Identified Using Co-expression Structure Network Analysis

Prostate cancer (PCa) is one of the most common malignancies for males, but very little is known about its pathogenesis. This study aimed to identify novel biomarkers associated with PCa prognosis and elucidate the underlying molecular mechanism. First, The Cancer Genome Atlas (TCGA) RNA-sequencing data were utilized to identify differentially expressed genes (DEGs) between tumor and normal samples. The DEGs were then applied to construct a co-expression and mined using structure network analysis. The magenta module that was highly related to the Gleason score (r = 0.46, p = 3e-26) and tumor stage (r = 0.38, p = 2e-17) was screened. Subsequently, all genes of the magenta module underwent function annotation. From the key module, CCNA2, CKAP2L, NCAPG, and NUSAP1 were chosen as the four candidate genes. Finally, internal (TCGA) and external data sets (GSE32571, GSE70770, and GSE141551) were combined to validate and predict the value of real hub genes. The results show that the above genes are up-regulated in PCa samples, and higher expression levels show significant association with higher Gleason scores and tumor T stage. Moreover, receiver operating characteristic curve and survival analysis validate the excellent value of hub genes in PCa progression and prognosis. In addition, the protein levels of these four genes also remain higher in tumor tissues when compared with normal tissues. Gene set enrichment analysis and gene set variation analysis for a single gene reveal the close relation with cell proliferation. Meanwhile, 11 small molecular drugs that have the potential to treat PCa were also screened. In conclusion, our research identified four potential prognostic genes and several candidate molecular drugs for treating PCa.

Analysis of the PD-1 Ligands Among Gastrointestinal Cancer Patients: Focus on Cancer Immunity

Many types of gastrointestinal cancer have shown promising outcomes after checkpoint blockade immunotherapy; however, it remains largely unclear about the expression profiles of programmed death 1 (PD-1) ligands (CD274 and PDCD1LG2) in the context of human pan-cancer. This work comprehensively analyzed the expression pattern of the PD-1 ligands and the clinical significance in the prognosis prediction among the seven types of gastrointestinal malignancies collected from The Cancer Genome Atlas (TCGA) and the Cancer Cell Line Encyclopedia (CCLE) database. Furthermore, the correlation of CD274/PDCD1LG2 with cancer immunity was also explored. The patients with liver hepatocellular carcinoma (LIHC) receiving cytokine-induced killer (CIK) cell immunotherapy at our cancer center were enrolled. CD274 and PDCD1LG2 displayed inconsistent gene expression levels among the diverse cancer cell lines. Typically, the abnormal expression level of CD274 and PDCD1LG2 was detected in both esophageal carcinoma (ESCA) and stomach adenocarcinoma (STAD), where PDCD1LG2 was related to the overall survival (OS) of the patients in ESCA (p = 0.015) and STAD (p = 0.025). High-serum CD274 and PDCD1LG2 levels predicted a worse survival in the patients with LIHC receiving CIK therapy. More importantly, the expression level of CD274 and PDCD1LG2 was significantly correlated with the degree of Estimation of STromal and Immune cells in MAlignant Tumor tissues using the Expression data (ESTIMATE). In addition, we found that CD274 and PDCD1LG2 were correlated with gene markers in tumor-infiltrating immune cells. Furthermore, the expression of CD274 and PDCD1LG2 was correlated with tumor mutation burden (TMB), microsatellite instability (MSI), mismatch repair (MMR), and DNA methyltransferase (DNMT) of different types of cancers. The present work comprehensively analyzed a RNA sequencing of the PD-1 ligands across the seven distinct types of gastrointestinal cancers, which provided clues for further studies in cancer immunity and development.

An integrated pan-cancer analysis of TFAP4 aberrations and the potential clinical implications for cancer immunity

Studies have shown that transcription factor activating enhancer binding protein 4 (TFAP4) plays a vital role in multiple types of cancer; however, the TFAP4 expression profile is still unknown, as is its value within the human pan‐cancer analysis. The present study comprehensively analysed TFAP4 expression patterns from 33 types of malignancies, along with the significance of TFAP4 for prognosis prediction and cancer immunity. TFAP4 displayed inconsistent levels of gene expression across the diverse cancer cell lines, and displayed abnormal expression within most malignant tumours, which closely corresponded to overall survival. More importantly, the TFAP4 level was also significantly related to the degree of tumour infiltration. TFAP4 was correlated using gene markers in tumour‐infiltrating immune cells and immune scores. TFAP4 expression was correlated with tumour mutation burden and microsatellite instability in different cancer types, and enrichment analyses identified TFAP4‐associated terms and pathways. The present study comprehensively analysed the expression of TFAP4 across 33 distinct types of cancers, which revealed that TFAP4 may possibly play a vital role during cancer formation and development. TFAP4 is related to differing degrees of immune infiltration within cancers, which suggests the potential of TFAP4 as an immunotherapy target in cancers. Our study demonstrated that TFAP4 plays an important role in tumorigenesis as a prognostic biomarker, which highlights the possibility of developing new targeted treatments.

Multi-Approach Bioinformatics Analysis of Curated Omics Data Provides a Gene Expression Panorama for Multiple Cancer Types

Studies describing the expression patterns and biomarkers for the tumoral process increase in number every year. The availability of new datasets, although essential, also creates a confusing landscape where common or critical mechanisms are obscured amidst the divergent and heterogeneous nature of such results. In this work, we manually curated the Gene Expression Omnibus using rigorous filtering criteria to select the most homogeneous and highest quality microarray and RNA-seq datasets from multiple types of cancer. By applying systems biology approaches, combined with machine learning analysis, we investigated possible frequently deregulated molecular mechanisms underlying the tumoral process. Our multi-approach analysis of 99 curated datasets, composed of 5,406 samples, revealed 47 differentially expressed genes in all analyzed cancer types, which were all in agreement with the validation using TCGA data. Results suggest that the tumoral process is more related to the overexpression of core deregulated machinery than the underexpression of a given gene set. Additionally, we identified gene expression similarities between different cancer types not described before and performed an overall survival analysis using 20 cancer types. Finally, we were able to suggest a core regulatory mechanism that could be frequently deregulated.

DORGE: Discovery of Oncogenes and tumoR suppressor genes using Genetic and Epigenetic features

Data-driven discovery of cancer driver genes, including tumor suppressor genes (TSGs) and oncogenes (OGs), is imperative for cancer prevention, diagnosis, and treatment. Although epigenetic alterations are important for tumor initiation and progression, most known driver genes were identified based on genetic alterations alone. Here, we developed an algorithm, DORGE (Discovery of Oncogenes and tumor suppressoR genes using Genetic and Epigenetic features), to identify TSGs and OGs by integrating comprehensive genetic and epigenetic data. DORGE identified histone modifications as strong predictors for TSGs, and it found missense mutations, super enhancers, and methylation differences as strong predictors for OGs. We extensively validated DORGE-predicted cancer driver genes using independent functional genomics data. We also found that DORGE-predicted dual-functional genes (both TSGs and OGs) are enriched at hubs in protein-protein interaction and drug-gene networks. Overall, our study has deepened the understanding of epigenetic mechanisms in tumorigenesis and revealed previously undetected cancer driver genes.

Clinical Implications of Aberrant PD-1 and CTLA4 Expression for Cancer Immunity and Prognosis: A Pan-Cancer Study

Combination therapy with inhibitors of cytotoxic T lymphocyte-associated protein (CTLA)4 and programmed death (PD)-1 has demonstrated efficacy in cancer patients. However, there is little information on CTLA4 and PD-1 expression levels and their clinical significance across diverse cancers. In this study, we addressed this question by analyzing PD-1 and CTLA4 levels in 33 different types of cancer along with their prognostic significance using The Cancer Genome Atlas (TCGA) and Cancer Cell Line Encyclopedia datasets. Liver hepatocellular carcinoma (LIHC) patients receiving cytokine-induced killer cell (CIK) immunotherapy at Sun Yat-sen University cancer center were enrolled for survival analysis. The correlation between PD-1/CTLA4 expression and cancer immunity was also analyzed. The results showed that PD-1 and CTLA4 transcript levels varied across cancer cell lines, with aberrant expression detected in certain cancer types; Kaplan–Meier analysis with the Cox proportional hazards model showed that this was closely related to overall survival in breast invasive carcinoma, glioblastoma multiforme, head and neck squamous cell carcinoma, acute myeloid leukemialymphoma, uterine corpus endometrial carcinoma, and uveal melanoma in TCGA. High serum PD-1 and CTLA4 levels predicted better survival in LIHC patients receiving CIK therapy. PD-1 and CTLA4 levels were found to be significantly correlated with the degree of tumor cell infiltration using Tumor Immune Estimation Resource, Estimating Relative Subsets of RNA Transcripts, and Estimation of Stromal and immune Cells in Malignant Tumor Tissues Using Expression Data as well as with tumor-infiltrating lymphocyte marker expression; they were also related to tumor mutation burden, microsatellite instability, mismatch repair, and the expression of DNA methyltransferases in some cancer types. Gene set enrichment analysis of 33 cancer types provided further evidence for associations between PD-1/CTLA4 levels and cancer development and immunocyte infiltration. Thus, PD-1 and CTLA4 play important roles in tumorigenesis and tumor immunity and can serve as prognostic biomarkers in different cancer types.

When information is the treatment? Precision medicine in healthcare

Profoundly more data-intensive than conventional medicine, precision medicine's distinctive informational needs present new challenges for healthcare management. Data protection and privacy law are key determinants in precision medicine's future. This article examines legal and regulatory barriers to the incorporation of precision medicine into healthcare. Specific attention is paid to analyzing recent health privacy laws, court cases, and medical device regulations. Considering the challenges identified, recommendations and guidance are crafted for health leaders with reference to domestic and international initiatives.

NCAPG Is a Promising Therapeutic Target Across Different Tumor Types

Background

With the advent of CRISPR-Cas9 genome editing tool in gene therapy, identification of aberrantly expressed genes is of great value across various cancer types. Since a large number of patients may benefit from molecular targeted gene therapy. The purpose of this study was to identify aberrantly expressed genes across various cancer types, analyze prospective mechanisms and their correlation with survival outcomes.

Results

NCAPG was highly expressed in The Cancer Genome Atlas (TCGA) database, which includes the transcriptomes of 6,647 cancer and 647 normal tissue samples from 16 cancer types. Furthermore, a predicted NCAPG overexpression rate was also observed at the protein level in 16 tumor types. Importantly, high NCAPG level was significantly associated with unfavorable survival in various cancer types such as hepatocellular carcinoma (HCC), breast, lung or ovarian cancer. The multivariate analyses demonstrated that NCAPG, TNM, and Barcelona Clinic Liver Cancer (BCLC) staging were independent risk factors for mortality of patients with HCC. Moreover, functional and pathway enrichment analysis suggested that NCAPG was closely correlated with the pathways of cell cycle, cellular senescence, and mismatch repair. By weighted gene co-expression network analysis (WGCNA), we identified NCAPG as a hub gene in the turquoise module mostly related to the survival time of HCC samples.

Conclusion

To our knowledge, this study represents a comprehensive RNA-Seq analysis of several tumor types, revealing NCAPG as a promising molecular target. NCAPG overexpression may play important roles in carcinogenesis and progression of tumors via regulating tumor-related pathways, thereby broadening the understanding of the pathogenic mechanisms and highlighting the possibility of developing novel targeted therapeutics.

Pan-cancer mapping of differential protein-protein interactions

Deciphering the variations in the protein interactome is required to reach a systems-level understanding of tumorigenesis. To accomplish this task, we have considered the clinical and transcriptome data on >6000 samples from The Cancer Genome Atlas for 12 different cancers. Utilizing the gene expression levels as a proxy, we have identified the differential protein-protein interactions in each cancer type and presented a differential view of human protein interactome among the cancers. We clearly demonstrate that a certain fraction of proteins differentially interacts in the cancers, but there was no general protein interactome profile that applied to all cancers. The analysis also provided the characterization of differentially interacting proteins (DIPs) representing significant changes in their interaction patterns during tumorigenesis. In addition, DIP-centered protein modules with high diagnostic and prognostic performances were generated, which might potentially be valuable in not only understanding tumorigenesis, but also developing effective diagnosis, prognosis, and treatment strategies.

Theranostic application of miR-429 in HER2+ breast cancer

Human epidermal growth factor receptor 2 (HER2) is overexpressed/amplified in one third of breast cancers (BCs), and is associated with the poorer prognosis and the higher metastatic potential in BC. Emerging evidences highlight the role of microRNAs (miRNAs) in the regulation of several cellular processes, including BC.

Methods: Here we identified, by in silico approach, a group of three miRNAs with central biological role (high degree centrality) in HER2+ BC. We validated their dysregulation in HER2+ BC and we analysed their functional role by in vitro approaches on selected cell lines and by in vivo experiments in an animal model.

Results: We found that their expression is dysregulated in both HER2+ BC cell lines and human samples. Focusing our study on the only upregulated miRNA, miR-429, we discovered that it acts as an oncogene and its upregulation is required for HER2+ cell proliferation. It controls the metastatic potential of HER2+ BC subtype by regulating migration and invasion of the cell.

Conclusions: In HER2+ BC oncogenic miR-429 is able to regulate HIF1α pathway by directly targeting VHL mRNA, a molecule important for the degradation of HIF1α. The overexpression of miR-429, observed in HER2+ BC, causes increased proliferation and migration of the BC cells. More important, silencing miR-429 succeeds in delaying tumor growth, thus miR-429 could be proposed as a therapeutic probe in HER2+ BC tumors.

Horizontal and vertical integrative analysis methods for mental disorders omics data

In recent biomedical studies, omics profiling has been extensively conducted on various types of mental disorders. In most of the existing analyses, a single type of mental disorder and a single type of omics measurement are analyzed. In the study of other complex diseases, integrative analysis, both vertical and horizontal integration, has been conducted and shown to bring significantly new insights into disease etiology, progression, biomarkers, and treatment. In this article, we showcase the applicability of integrative analysis to mental disorders. In particular, the horizontal integration of bipolar disorder and schizophrenia and the vertical integration of gene expression and copy number variation data are conducted. The analysis is based on the sparse principal component analysis, penalization, and other advanced statistical techniques. In data analysis, integration leads to biologically sensible findings, including the disease-related gene expressions, copy number variations, and their associations, which differ from the “benchmark” analysis. Overall, this study suggests the potential of integrative analysis in mental disorder research.

Integrative Analysis of Cancer Omics Data for Prognosis Modeling

Prognosis modeling plays an important role in cancer studies. With the development of omics profiling, extensive research has been conducted to search for prognostic markers for various cancer types. However, many of the existing studies share a common limitation by only focusing on a single cancer type and suffering from a lack of sufficient information. With potential molecular similarity across cancer types, one cancer type may contain information useful for the analysis of other types. The integration of multiple cancer types may facilitate information borrowing so as to more comprehensively and more accurately describe prognosis. In this study, we conduct marginal and joint integrative analysis of multiple cancer types, effectively introducing integration in the discovery process. For accommodating high dimensionality and identifying relevant markers, we adopt the advanced penalization technique which has a solid statistical ground. Gene expression data on nine cancer types from The Cancer Genome Atlas (TCGA) are analyzed, leading to biologically sensible findings that are different from the alternatives. Overall, this study provides a novel venue for cancer prognosis modeling by integrating multiple cancer types.

Extracting predictors for lung adenocarcinoma based on Granger causality test and stepwise character selection

Background

Lung adenocarcinoma is the most common type of lung cancer, with high mortality worldwide. Its occurrence and development were thoroughly studied by high-throughput expression microarray, which produced abundant data on gene expression, DNA methylation, and miRNA quantification. However, the hub genes, which can be served as bio-markers for discriminating cancer and healthy individuals, are not well screened.

Result

Here we present a new method for extracting gene predictors, aiming to obtain the least predictors without losing the efficiency. We firstly analyzed three different expression microarrays and constructed multi-interaction network, since the individual expression dataset is not enough for describing biological behaviors dynamically and systematically. Then, we transformed the undirected interaction network to directed network by employing Granger causality test, followed by the predictors screened with the use of the stepwise character selection algorithm. Six predictors, including TOP2A, GRK5, SIRT7, MCM7, EGFR, and COL1A2, were ultimately identified. All the predictors are the cancer-related, and the number is very small fascinating diagnosis. Finally, the validation of this approach was verified by robustness analyses applied to six independent datasets; the precision is up to 95.3% ∼ 100%.

Conclusion

Although there are complicated differences between cancer and normal cells in gene functions, cancer cells could be differentiated in case that a group of special genes expresses abnormally. Here we presented a new, robust, and effective method for extracting gene predictors. We identified as low as 6 genes which can be taken as predictors for diagnosing lung adenocarcinoma.

An Integrated Approach to Unravel Hidradenitis Suppurativa Etiopathogenesis

Hidradenitis suppurativa/acne inversa (HS) is a chronic inflammatory disease involving hair follicles that presents with painful nodules, abscesses, fistulae, and hypertrophic scars, typically occurring in apocrine gland bearing skin. Establishing a diagnosis of HS may take up to 7 years after disease onset. HS severely impairs the quality of life of patients and its high frequency causes significant costs for health care system. HS patients have an increased risk of developing associated diseases, such as inflammatory bowel diseases and spondyloarthropathies, thereby suggesting a common pathophysiological mechanism. Familial cases, which are around 35% of HS patients, have allowed the identification of susceptibility genes. HS is perceived as a complex disease where environmental factors trigger chronic inflammation in the skin of genetically predisposed individuals. Despite the efforts made to understand HS etiopathogenesis, the exact mechanisms at the basis of the disease need to be still unraveled. In this review, we considered all OMICs studies performed on HS and observed that OMICs contribution in the context of HS appeared as not clear enough and/or rich of useful clinical information. Indeed, most studies focused only on one aspect—genome, transcriptome, or proteome—of the disease, enrolling small numbers of patients. This is quite limiting for the genetic studies, from different geographical areas and looking at a few aspects of HS pathogenesis without any integration of the findings obtained or a comparison among different studies. A strong need for an integrated approach using OMICs tools is required to discover novel actors involved in HS etiopathogenesis. Moreover, we suggest the constitution of consortia to enroll a higher number of patients to be analyzed following common and consensus OMICs strategies. Comparison and integration with the findings present in the OMICs repositories are mandatory. In a theoretic pipeline, the Skin-OMICs profile obtained from each HS patient should be compared and integrated with repositories and literature data by using appropriate InterOMICs approach. The final goal is not only to improve the knowledge of HS etiopathogenesis but also to provide novel tools to the clinicians with the eventual aim of offering a tailored treatment for HS patients.

Pan-Cancer analysis of the expression and regulation of matrisome genes across 32 tumor types

The microenvironment plays a central role in cancer, and neoplastic cells actively shape it to their needs by complex arrays of extracellular matrix (ECM) proteins, enzymes, cytokines and growth factors collectively referred to as the matrisome. Studies on the cancer matrisome have been performed for single or few neoplasms, but a more systematic analysis is still missing. Here we present a Pan-Cancer study of matrisome gene expression in 10,487 patients across 32 tumor types, supplemented with transcription factors (TFs) and driver genes/pathways regulating each tumor's matrisome. We report on 919 TF-target pairs, either used specifically or shared across tumor types, and their prognostic significance, 40 master regulators, 31 overarching regulatory pathways and the potential for druggability with FDA-approved cancer drugs. These results provide a comprehensive transcriptional architecture of the cancer matrisome and suggest the need for development of specific matrisome-targeting approaches for future therapies.

•

In-depth characterization of matrisome gene expression and regulation in 10,487 patients across 32 human tumor types.

•

Identification of transcription factor (TF) and “master regulators” governing each cancer’s matrisome.

•

Analysis unveils therapeutic possibilities and suggests new treatments by repurposing of FDA-approved cancer drugs.

In silico perturbation of drug targets in pan-cancer analysis combining multiple networks and pathways

The knowledge of cancer cell response to conventional therapies is crucial in order to choose the correct therapy of patients affected by cancer. The major problem is generally attributed to the lack of specific biological processes able to predict the therapy efficacy. Here, we optimized a computational method for the analysis of gene networks able to detect and quantify the effects of a drug in a pan-cancer study. Overall, our method, using several network topological measures has identified a cancer gene network with a key role in biological processes. The gene network, able to classify with a good performance cancer vs normal samples, was modulated in silico to evaluate the effects of new or approved drugs. This computational model could offer an interesting hint to decipher molecular mechanisms contributing to resistance or inefficacy of drugs.

Translating cancer genomics into precision medicine with artificial intelligence: applications, challenges and future perspectives

In the field of cancer genomics, the broad availability of genetic information offered by next-generation sequencing technologies and rapid growth in biomedical publication has led to the advent of the big-data era. Integration of artificial intelligence (AI) approaches such as machine learning, deep learning, and natural language processing (NLP) to tackle the challenges of scalability and high dimensionality of data and to transform big data into clinically actionable knowledge is expanding and becoming the foundation of precision medicine. In this paper, we review the current status and future directions of AI application in cancer genomics within the context of workflows to integrate genomic analysis for precision cancer care. The existing solutions of AI and their limitations in cancer genetic testing and diagnostics such as variant calling and interpretation are critically analyzed. Publicly available tools or algorithms for key NLP technologies in the literature mining for evidence-based clinical recommendations are reviewed and compared. In addition, the present paper highlights the challenges to AI adoption in digital healthcare with regard to data requirements, algorithmic transparency, reproducibility, and real-world assessment, and discusses the importance of preparing patients and physicians for modern digitized healthcare. We believe that AI will remain the main driver to healthcare transformation toward precision medicine, yet the unprecedented challenges posed should be addressed to ensure safety and beneficial impact to healthcare.

An Introduction to the Mathematical Modeling in the Study of Cancer Systems Biology

BACKGROUND

Frequently occurring in cancer are the aberrant alterations of regulatory onco-metabolites, various oncogenes/epigenetic stochasticity, and suppressor genes, as well as the deficient mismatch repair mechanism, chronic inflammation, or those deviations belonging to the other cancer characteristics. How these aberrations that evolve overtime determine the global phenotype of malignant tumors remains to be completely understood. Dynamic analysis may have potential to reveal the mechanism of carcinogenesis and can offer new therapeutic intervention.

AIMS

We introduce simplified mathematical tools to model serial quantitative data of cancer biomarkers. We also highlight an introductory overview of mathematical tools and models as they apply from the viewpoint of known cancer features.

METHODS

Mathematical modeling of potentially actionable genomic products and how they proceed overtime during tumorigenesis are explored. This report is intended to be instinctive without being overly technical.

RESULTS

To date, many mathematical models of the common features of cancer have been developed. However, the dynamic of integrated heterogeneous processes and their cross talks related to carcinogenesis remains to be resolved.

CONCLUSIONS

In cancer research, outlining mathematical modeling of experimentally obtained data snapshots of molecular species may provide insights into a better understanding of the multiple biochemical circuits. Recent discoveries have provided support for the existence of complex cancer progression in dynamics that span from a simple 1-dimensional deterministic system to a stochastic (ie, probabilistic) or to an oscillatory and multistable networks. Further research in mathematical modeling of cancer progression, based on the evolving molecular kinetics (time series), could inform a specific and a predictive behavior about the global systems biology of vulnerable tumor cells in their earlier stages of oncogenesis. On this footing, new preventive measures and anticancer therapy could then be constructed.

Synstable Fusion: A Network-Based Algorithm for Estimating Driver Genes in Fusion Structures

Gene fusion structure is a class of common somatic mutational events in cancer genomes, which are often formed by chromosomal mutations. Identifying the driver gene(s) in a fusion structure is important for many downstream analyses and it contributes to clinical practices. Existing computational approaches have prioritized the importance of oncogenes by incorporating prior knowledge from gene networks. However, different methods sometimes suffer different weaknesses when handling gene fusion data due to multiple issues such as fusion gene representation, network integration, and the effectiveness of the evaluation algorithms. In this paper, Synstable Fusion (SYN), an algorithm for computationally evaluating the fusion genes, is proposed. This algorithm uses network-based strategy by incorporating gene networks as prior information, but estimates the driver genes according to the destructiveness hypothesis. This hypothesis balances the two popular evaluation strategies in the existing studies, thereby providing more comprehensive results. A machine learning framework is introduced to integrate multiple networks and further solve the conflicting results from different networks. In addition, a synchronous stability model is established to reduce the computational complexity of the evaluation algorithm. To evaluate the proposed algorithm, we conduct a series of experiments on both artificial and real datasets. The results demonstrate that the proposed algorithm performs well on different configurations and is robust when altering the internal parameter settings.

In silico identification of drug target pathways in breast cancer subtypes using pathway cross-talk inhibition

Background

Despite great development in genome and proteome high-throughput methods, treatment failure is a critical point in the management of most solid cancers, including breast cancer (BC). Multiple alternative mechanisms upon drug treatment are involved to offset therapeutic effects, eventually causing drug resistance or treatment failure.

Methods

Here, we optimized a computational method to discover novel drug target pathways in cancer subtypes using pathway cross-talk inhibition (PCI). The in silico method is based on the detection and quantification of the pathway cross-talk for distinct cancer subtypes. From a BC data set of The Cancer Genome Atlas, we have identified different networks of cross-talking pathways for different BC subtypes, validated using an independent BC dataset from Gene Expression Omnibus. Then, we predicted in silico the effects of new or approved drugs on different BC subtypes by silencing individual or combined subtype-derived pathways with the aim to find new potential drugs or more effective synergistic combinations of drugs.

Results

Overall, we identified a set of new potential drug target pathways for distinct BC subtypes on which therapeutic agents could synergically act showing antitumour effects and impacting on cross-talk inhibition.

Conclusions

We believe that in silico methods based on PCI could offer valuable approaches to identifying more tailored and effective treatments in particular in heterogeneous cancer diseases.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1535-2) contains supplementary material, which is available to authorized users.

How do cancer cells replenish their fuel supply?

BACKGROUND

Multiple genetic changes, availability of cellular nutrients and metabolic alterations play a pivotal role in oncogenesis AIMS: We focus on cancer cell's metabolic properties, and we outline the cross talks between cellular oncogenic growth pathways in cancer metabolism. The review also provides a synopsis of the relevant cancer drugs targeting metabolic activities that are at various stages of clinical development.

METHODS

We review literature published within the last decade to include select articles that have highlighted energy metabolism crucial to the development of cancer phenotypes.

RESULTS

Cancer cells maintain their potent metabolism and keep a balanced redox status by enhancing glycolysis and autophagy and rerouting Krebs cycle intermediates and products of β-oxydation.

CONCLUSIONS

The processes underlying cancer pathogenesis are extremely complex and remain elusive. The new field of systems biology provides a mathematical framework in which these homeostatic dysregulation principles may be examined for better understanding of cancer phenotypes. Knowledge of key players in cancer-related metabolic reprograming may pave the way for new therapeutic metabolism-targeted drugs and ultimately improve patient care.

Radiation-Induced Gene Expression Changes in High and Low Grade Breast Cancer Cell Types

BACKGROUND

There is extensive scientific evidence that radiation therapy (RT) is a crucial treatment, either alone or in combination with other treatment modalities, for many types of cancer, including breast cancer (BC). BC is a heterogeneous disease at both clinical and molecular levels, presenting distinct subtypes linked to the hormone receptor (HR) status and associated with different clinical outcomes. The aim of this study was to assess the molecular changes induced by high doses of ionizing radiation (IR) on immortalized and primary BC cell lines grouped according to Human epidermal growth factor receptor (HER2), estrogen, and progesterone receptors, to study how HR status influences the radiation response. Our genomic approach using in vitro and ex-vivo models (e.g., primary cells) is a necessary first step for a translational study to describe the common driven radio-resistance features associated with HR status. This information will eventually allow clinicians to prescribe more personalized total doses or associated targeted therapies for specific tumor subtypes, thus enhancing cancer radio-sensitivity.

METHODS

Nontumorigenic (MCF10A) and BC (MCF7 and MDA-MB-231) immortalized cell lines, as well as healthy (HMEC) and BC (BCpc7 and BCpcEMT) primary cultures, were divided into low grade, high grade, and healthy groups according to their HR status. At 24 h post-treatment, the gene expression profiles induced by two doses of IR treatment with 9 and 23 Gy were analyzed by cDNA microarray technology to select and compare the differential gene and pathway expressions among the experimental groups.

RESULTS

We present a descriptive report of the substantial alterations in gene expression levels and pathways after IR treatment in both immortalized and primary cell cultures. Overall, the IR-induced gene expression profiles and pathways appear to be cell-line dependent. The data suggest that some specific gene and pathway signatures seem to be linked to HR status.

CONCLUSIONS

Genomic biomarkers and gene-signatures of specific tumor subtypes, selected according to their HR status and molecular features, could facilitate personalized biological-driven RT treatment planning alone and in combination with targeted therapies.