CRUX, a platform for visualising, exploring and analysing cancer genome cohort data

To better understand how tumours develop, identify prognostic biomarkers and find new treatments, researchers have generated vast catalogues of cancer genome data. However, these datasets are complex, so interpreting their important features requires specialized computational skills and analytical tools, which presents a significant technical challenge. To address this, we developed CRUX, a platform for exploring genomic data from cancer cohorts. CRUX enables researchers to perform common analyses including cohort comparisons, biomarker discovery, survival analysis, and to create visualisations including oncoplots and lollipop charts. CRUX simplifies cancer genome analysis in several ways: (i) it has an easy-to-use graphical interface; (ii) it enables users to create custom cohorts, as well as analyse precompiled public and private user-created datasets; (iii) it allows analyses to be run locally to address data privacy concerns (though an online version is also available) and (iv) it makes it easy to use additional specialized tools by exporting data in the correct formats. We showcase CRUX's capabilities with case studies employing different types of cancer genome analysis, demonstrating how it can be used flexibly to generate valuable insights into cancer biology. CRUX is freely available at https://github.com/CCICB/CRUX and https://ccicb.shinyapps.io/crux (DOI: 10.5281/zenodo.8015714).

Nursing Informatics and Epigenetics: Methodological Considerations for Big Data Analysis

Nursing informatics requires an understanding of patient-centered data and clinical workflow, and epigenetic research requires an understanding of data analysis. The purpose of this article is to document the methodology that nursing informatics specialists can use to conduct epigenetic research and subsequently strengthen patient-centered care. A pilot study of a secondary methylation data analysis using The Cancer Genome Atlas data from individuals with colon cancer is utilized to illustrate the methodology. The steps for conducting the study using public and free resources are discussed. These steps include finding a data source; downloading and analyzing differentially methylated regions; annotating differentially methylated region, gene ontology and function analysis; and reporting results. A model of epigenetic testing workflow is provided, as is a list of publicly available data and analysis sources that can be used to conduct epigenetic research.

Obesity-Associated Differentially Methylated Regions in Colon Cancer

Obesity with adiposity is a common disorder in modern days, influenced by environmental factors such as eating and lifestyle habits and affecting the epigenetics of adipose-based gene regulations and metabolic pathways in colorectal cancer (CRC). We compared epigenetic changes of differentially methylated regions (DMR) of genes in colon tissues of 225 colon cancer cases (154 non-obese and 71 obese) and 15 healthy non-obese controls by accessing The Cancer Genome Atlas (TCGA) data. We applied machine-learning-based analytics including generalized regression (GR) as a confirmatory validation model to identify the factors that could contribute to DMRs impacting colon cancer to enhance prediction accuracy. We found that age was a significant predictor in obese cancer patients, both alone (p = 0.003) and interacting with hypomethylated DMRs of ZBTB46, a tumor suppressor gene (p = 0.008). DMRs of three additional genes: HIST1H3I (p = 0.001), an oncogene with a hypomethylated DMR in the promoter region; SRGAP2C (p = 0.006), a tumor suppressor gene with a hypermethylated DMR in the promoter region; and NFATC4 (p = 0.006), an adipocyte differentiating oncogene with a hypermethylated DMR in an intron region, are also significant predictors of cancer in obese patients, independent of age. The genes affected by these DMR could be potential novel biomarkers of colon cancer in obese patients for cancer prevention and progression.

Multiomics Topic Modeling for Breast Cancer Classification

The integration of transcriptional data with other layers of information, such as the post-transcriptional regulation mediated by microRNAs, can be crucial to identify the driver genes and the subtypes of complex and heterogeneous diseases such as cancer. This paper presents an approach based on topic modeling to accomplish this integration task. More specifically, we show how an algorithm based on a hierarchical version of stochastic block modeling can be naturally extended to integrate any combination of 'omics data. We test this approach on breast cancer samples from the TCGA database, integrating data on messenger RNA, microRNAs, and copy number variations. We show that the inclusion of the microRNA layer significantly improves the accuracy of subtype classification. Moreover, some of the hidden structures or "topics" that the algorithm extracts actually correspond to genes and microRNAs involved in breast cancer development and are associated to the survival probability.

High mutation burden in the checkpoint and micro-RNA processing genes in myelodysplastic syndrome

A number of sequencing studies identified the prognostic impact of somatic mutations in myelodysplastic syndrome (MDS). However the majority of them focused on methylation regulation, apoptosis and proliferation genes. Despite the number of experimental studies published on the role of micro-RNA processing and checkpoint genes in the development of MDS, the clinical data about mutational landscape in these genes is limited. We performed a pilot study which evaluated mutational burden in these genes and their association with common MDS mutations. High prevalence of mutations was observed in the genes studied: 54% had mutations in DICER1, 46% had mutations in LAG3, 20% in CTLA4, 23% in B7-H3, 17% in DROSHA, 14% in PD-1 and 3% in PD-1L. Cluster analysis that included these mutations along with mutations in ASXL1, DNMT3A, EZH2, IDH1, RUNX1, SF3B1, SRSF2, TET2 and TP53 effectively predicted overall survival in the study group (HR 4.2, 95%CI 1.3-13.6, p = 0.016). The study results create the rational for incorporating micro-RNA processing and checkpoint genes in the sequencing panels for MDS and evaluate their role in the multicenter studies.

A Topic Modeling Analysis of TCGA Breast and Lung Cancer Transcriptomic Data

Topic modeling is a widely used technique to extract relevant information from large arrays of data. The problem of finding a topic structure in a dataset was recently recognized to be analogous to the community detection problem in network theory. Leveraging on this analogy, a new class of topic modeling strategies has been introduced to overcome some of the limitations of classical methods. This paper applies these recent ideas to TCGA transcriptomic data on breast and lung cancer. The established cancer subtype organization is well reconstructed in the inferred latent topic structure. Moreover, we identify specific topics that are enriched in genes known to play a role in the corresponding disease and are strongly related to the survival probability of patients. Finally, we show that a simple neural network classifier operating in the low dimensional topic space is able to predict with high accuracy the cancer subtype of a test expression sample.

ELMER v.2: an R/Bioconductor package to reconstruct gene regulatory networks from DNA methylation and transcriptome profiles

Motivation

DNA methylation has been used to identify functional changes at transcriptional enhancers and other cis-regulatory modules (CRMs) in tumors and other disease tissues. Our R/Bioconductor package ELMER (Enhancer Linking by Methylation/Expression Relationships) provides a systematic approach that reconstructs altered gene regulatory networks (GRNs) by combining enhancer methylation and gene expression data derived from the same sample set.

Results

We present a completely revised version 2 of ELMER that provides numerous new features including an optional web-based interface and a new Supervised Analysis mode to use pre-defined sample groupings. We show that Supervised mode significantly increases statistical power and identifies additional GRNs and associated Master Regulators, such as SOX11 and KLF5 in Basal-like breast cancer.

Availability and implementation

ELMER v.2 is available as an R/Bioconductor package at http://bioconductor.org/packages/ELMER/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Schwann cell reprogramming and lung cancer progression: a meta-analysis of transcriptome data

Schwann cells were identified in the tumor surrounding area prior to initiate the invasion process underlying connective tissue. These cells promote cancer invasion through direct contact, while paracrine signaling and matrix remodeling are not sufficient to proceed. Considering the intertwined structure of signaling, regulatory, and metabolic processes within a cell, we employed a genome-scale biomolecular network. Accordingly, a meta-analysis of Schwann cells associated transcriptomic datasets was performed, and the core information on differentially expressed genes (DEGs) was obtained by statistical analyses. Gene set over-representation analyses was performed on core DEGs to identify significantly functional and pathway enrichment analysis between Schwann cells and, lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). DEGs were further integrated with genome-scale human biomolecular networks. miRNAs were proposed by the reconstruction of a transcriptional and post-transcriptional regulatory network. Moreover, microarray-based transcriptome profiling was performed, and the prognostic power of selected dedifferentiated Schwann cell biomolecules was predicted. We observed that pathways associated with Schwann cells dedifferentiation was overexpressed in lung cancer samples. However, genes associated with Schwann cells migration inhibition system were downregulated. Besides, miRNA targeting those pathways were also deregulated. In this study, we report valuable data for further experimental and clinical analysis, because the proposed biomolecules have significant potential as systems biomarkers for screening or for therapeutic purposes in perineural invasion of lung cancer.

Correlation between renin-angiotensin system (RAS) related genes, type 2 diabetes, and cancer: Insights from metanalysis of transcriptomics data

Although studies have provided significant evidence about the role of RAS in mediating cancer risk in type 2 diabetes mellitus (DM), conclusions about the central molecular mechanisms underlying this disease remain to be reached, because this type of information requires an integrative multi-omics approach. In the current study, meta-analysis was performed on type 2 diabetes and breast, bladder, liver, pancreas, colon and rectum cancer-associated transcriptome data, and reporter biomolecules were identified at RNA, protein, and metabolite levels using the integration of gene expression profiles with genome-scale biomolecular networks in diabetes samples. This approach revealed that RAS biomarkers could be associated with cancer initiation and progression, which include metabolites (particularly, aminoacyl-tRNA biosynthesis and ABC transporters) as novel biomarker candidates and potential therapeutic targets. We detected downregulation and upregulation of differentially expressed genes (DEGs) in blood, pancreatic islets, liver and skeletal muscle from normal and diabetic patients. DEGs were combined with 211 renin-angiotensin-system related genes. Upregulated genes were enriched using Pathway analysis of cancer in pancreatic islets, blood and skeletal muscle samples. It seems that the changes in mRNA are contributing to the phenotypic changes in carcinogenesis, or that they are as a result of the phenotypic changes associated with the malignant transformation. Our analyses showed that Ctsg and Ednrb are downregulated in cancer samples. However, by immunohistochemistry experiments we observed that EDNRB protein showed increased expression in tumor samples. It is true that alterations in mRNA expression do not always reflect alterations in protein expression, since post-translational changes can occur in proteins. In this study, we report valuable data for further experimental and clinical analysis, because the proposed biomolecules have significant potential as systems biomarkers for screening or for therapeutic purposes in type 2 diabetes and cancer-associated pathways.