Pan-cancer analysis of pathway-based gene expression pattern at the individual level reveals biomarkers of clinical prognosis

Surgical Tumor Resection Deregulates Hallmarks of Cancer in Resected Tissue and the Surrounding Microenvironment

Surgery exposes tumor tissue to severe hypoxia and mechanical stress leading to rapid gene expression changes in the tumor and its microenvironment, which remain poorly characterized. We biopsied tumor and adjacent normal tissues from patients with breast (n = 81) and head/neck squamous cancers (HNSC; n = 10) at the beginning (A), during (B), and end of surgery (C). Tumor/normal RNA from 46/81 patients with breast cancer was subjected to mRNA-Seq using Illumina short-read technology, and from nine patients with HNSC to whole-transcriptome microarray with Illumina BeadArray. Pathways and genes involved in 7 of 10 known cancer hallmarks, namely, tumor-promoting inflammation (TNF-A, NFK-B, IL18 pathways), activation of invasion and migration (various extracellular matrix-related pathways, cell migration), sustained proliferative signaling (K-Ras Signaling), evasion of growth suppressors (P53 signaling, regulation of cell death), deregulating cellular energetics (response to lipid, secreted factors, and adipogenesis), inducing angiogenesis (hypoxia signaling, myogenesis), and avoiding immune destruction (CTLA4 and PDL1) were significantly deregulated during surgical resection (time points A vs. B vs. C). These findings were validated using NanoString assays in independent pre/intra/post-operative breast cancer samples from 48 patients. In a comparison of gene expression data from biopsy (analogous to time point A) with surgical resection samples (analogous to time point C) from The Cancer Genome Atlas study, the top deregulated genes were the same as identified in our analysis, in five of the seven studied cancer types. This study suggests that surgical extirpation deregulates the hallmarks of cancer in primary tumors and adjacent normal tissue across different cancers.

IMPLICATIONS

Surgery deregulates hallmarks of cancer in human tissue.

Targeting neddylation and sumoylation in chemoresistant triple negative breast cancer

Triple negative breast cancer (TNBC) accounts for 15-20% of breast cancer cases in the United States. Systemic neoadjuvant chemotherapy (NACT), with or without immunotherapy, is the current standard of care for patients with early-stage TNBC. However, up to 70% of TNBC patients have significant residual disease once NACT is completed, which is associated with a high risk of developing recurrence within two to three years of surgical resection. To identify targetable vulnerabilities in chemoresistant TNBC, we generated longitudinal patient-derived xenograft (PDX) models from TNBC tumors before and after patients received NACT. We then compiled transcriptomes and drug response profiles for all models. Transcriptomic analysis identified the enrichment of aberrant protein homeostasis pathways in models from post-NACT tumors relative to pre-NACT tumors. This observation correlated with increased sensitivity in vitro to inhibitors targeting the proteasome, heat shock proteins, and neddylation pathways. Pevonedistat, a drug annotated as a NEDD8-activating enzyme (NAE) inhibitor, was prioritized for validation in vivo and demonstrated efficacy as a single agent in multiple PDX models of TNBC. Pharmacotranscriptomic analysis identified a pathway-level correlation between pevonedistat activity and post-translational modification (PTM) machinery, particularly involving neddylation and sumoylation targets. Elevated levels of both NEDD8 and SUMO1 were observed in models exhibiting a favorable response to pevonedistat compared to those with a less favorable response in vivo. Moreover, a correlation emerged between the expression of neddylation-regulated pathways and tumor response to pevonedistat, indicating that targeting these PTM pathways may prove effective in combating chemoresistant TNBC.

Pathway Activation Analysis for Pan-Cancer Personalized Characterization Based on Riemannian Manifold

The pathogenesis of carcinoma is believed to come from the combined effect of polygenic variation, and the initiation and progression of malignant tumors are closely related to the dysregulation of biological pathways. Quantifying the alteration in pathway activation and identifying coordinated patterns of pathway dysfunction are the imperative part of understanding the malignancy process and distinguishing different tumor stages or clinical outcomes of individual patients. In this study, we have conducted in silico pathway activation analysis using Riemannian manifold (RiePath) toward pan-cancer personalized characterization, which is the first attempt to apply the Riemannian manifold theory to measure the extent of pathway dysregulation in individual patient on the tangent space of the Riemannian manifold. RiePath effectively integrates pathway and gene expression information, not only generating a relatively low-dimensional and biologically relevant representation, but also identifying a robust panel of biologically meaningful pathway signatures as biomarkers. The pan-cancer analysis across 16 cancer types reveals the capability of RiePath to evaluate pathway activation accurately and identify clinical outcome-related pathways. We believe that RiePath has the potential to provide new prospects in understanding the molecular mechanisms of complex diseases and may find broader applications in predicting biomarkers for other intricate diseases.

Interrogating the Human Diplome: Computational Methods, Emerging Applications, and Challenges

Human DNA sequencing protocols have revolutionized human biology, biomedical science, and clinical practice, but still have very important limitations. One limitation is that most protocols do not separate or assemble (i.e., "phase") the nucleotide content of each of the maternally and paternally derived chromosomal homologs making up the 22 autosomal pairs and the chromosomal pair making up the pseudo-autosomal region of the sex chromosomes. This has led to a dearth of studies and a consequent underappreciation of many phenomena of fundamental importance to basic and clinical genomic science. We discuss a few protocols for obtaining phase information as well as their limitations, including those that could be used in tumor phasing settings. We then describe a number of biological and clinical phenomena that require phase information. These include phenomena that require precise knowledge of the nucleotide sequence in a chromosomal segment from germline or somatic cells, such as DNA binding events, and insight into unique cis vs. trans-acting functionally impactful variant combinations-for example, variants implicated in a phenotype governed by compound heterozygosity. In addition, we also comment on the need for reliable and consensus-based diploid-context computational workflows for variant identification as well as the need for laboratory-based functional verification strategies for validating cis vs. trans effects of variant combinations. We also briefly describe available resources, example studies, as well as areas of further research, and ultimately argue that the science behind the study of human diploidy, referred to as "diplomics," which will be enabled by nucleotide-level resolution of phased genomes, is a logical next step in the analysis of human genome biology.

Single sample pathway analysis in metabolomics: performance evaluation and application

BACKGROUND

Single sample pathway analysis (ssPA) transforms molecular level omics data to the pathway level, enabling the discovery of patient-specific pathway signatures. Compared to conventional pathway analysis, ssPA overcomes the limitations by enabling multi-group comparisons, alongside facilitating numerous downstream analyses such as pathway-based machine learning. While in transcriptomics ssPA is a widely used technique, there is little literature evaluating its suitability for metabolomics. Here we provide a benchmark of established ssPA methods (ssGSEA, GSVA, SVD (PLAGE), and z-score) alongside the evaluation of two novel methods we propose: ssClustPA and kPCA, using semi-synthetic metabolomics data. We then demonstrate how ssPA can facilitate pathway-based interpretation of metabolomics data by performing a case-study on inflammatory bowel disease mass spectrometry data, using clustering to determine subtype-specific pathway signatures.

RESULTS

While GSEA-based and z-score methods outperformed the others in terms of recall, clustering/dimensionality reduction-based methods provided higher precision at moderate-to-high effect sizes. A case study applying ssPA to inflammatory bowel disease data demonstrates how these methods yield a much richer depth of interpretation than conventional approaches, for example by clustering pathway scores to visualise a pathway-based patient subtype-specific correlation network. We also developed the sspa python package (freely available at https://pypi.org/project/sspa/ ), providing implementations of all the methods benchmarked in this study.

CONCLUSION

This work underscores the value ssPA methods can add to metabolomic studies and provides a useful reference for those wishing to apply ssPA methods to metabolomics data.

Determination of a DNA repair-related gene signature with potential implications for prognosis and therapeutic response in pancreatic adenocarcinoma

Background

Pancreatic adenocarcinoma (PAAD) is one of the leading causes of cancer death worldwide. Alterations in DNA repair-related genes (DRGs) are observed in a variety of cancers and have been shown to affect the development and treatment of cancers. The aim of this study was to develop a DRG-related signature for predicting prognosis and therapeutic response in PAAD.

Methods

We constructed a DRG signature using least absolute shrinkage and selection operator (LASSO) Cox regression analysis in the TCGA training set. GEO datasets were used as the validation set. A predictive nomogram was constructed based on multivariate Cox regression. Calibration curve and decision curve analysis (DCA) were applied to validate the performance of the nomogram. The CIBERSORT and ssGSEA algorithms were utilized to explore the relationship between the prognostic signature and immune cell infiltration. The pRRophetic algorithm was used to estimate sensitivity to chemotherapeutic agents. The CellMiner database and PAAD cell lines were used to investigate the relationship between DRG expression and therapeutic response.

Results

We developed a DRG signature consisting of three DRGs (RECQL, POLQ, and RAD17) that can predict prognosis in PAAD patients. A prognostic nomogram combining the risk score and clinical factors was developed for prognostic prediction. The DCA curve and the calibration curve demonstrated that the nomogram has a higher net benefit than the risk score and TNM staging system. Immune infiltration analysis demonstrated that the risk score was positively correlated with the proportions of activated NK cells and monocytes. Drug sensitivity analysis indicated that the signature has potential predictive value for chemotherapy. Analyses utilizing the CellMiner database showed that RAD17 expression is correlated with oxaliplatin. The dynamic changes in three DRGs in response to oxaliplatin were examined by RT-qPCR, and the results show that RAD17 is upregulated in response to oxaliplatin in PAAD cell lines.

Conclusion

We constructed and validated a novel DRG signature for prediction of the prognosis and drug sensitivity of patients with PAAD. Our study provides a theoretical basis for further unraveling the molecular pathogenesis of PAAD and helps clinicians tailor systemic therapies within the framework of individualized treatment.

Molecular Biomarkers in Cancer

Molecular cancer biomarkers are any measurable molecular indicator of risk of cancer, occurrence of cancer, or patient outcome. They may include germline or somatic genetic variants, epigenetic signatures, transcriptional changes, and proteomic signatures. These indicators are based on biomolecules, such as nucleic acids and proteins, that can be detected in samples obtained from tissues through tumor biopsy or, more easily and non-invasively, from blood (or serum or plasma), saliva, buccal swabs, stool, urine, etc. Detection technologies have advanced tremendously over the last decades, including techniques such as next-generation sequencing, nanotechnology, or methods to study circulating tumor DNA/RNA or exosomes. Clinical applications of biomarkers are extensive. They can be used as tools for cancer risk assessment, screening and early detection of cancer, accurate diagnosis, patient prognosis, prediction of response to therapy, and cancer surveillance and monitoring response. Therefore, they can help to optimize making decisions in clinical practice. Moreover, precision oncology is needed for newly developed targeted therapies, as they are functional only in patients with specific cancer genetic mutations, and biomarkers are the tools used for the identification of these subsets of patients. Improvement in the field of cancer biomarkers is, however, needed to overcome the scientific challenge of developing new biomarkers with greater sensitivity, specificity, and positive predictive value.

SWAN pathway-network identification of common aneuploidy-based oncogenic drivers

Haploinsufficiency drives Darwinian evolution. Siblings, while alike in many aspects, differ due to monoallelic differences inherited from each parent. In cancer, solid tumors exhibit aneuploid genetics resulting in hundreds to thousands of monoallelic gene-level copy-number alterations (CNAs) in each tumor. Aneuploidy patterns are heterogeneous, posing a challenge to identify drivers in this high-noise genetic environment. Here, we developed Shifted Weighted Annotation Network (SWAN) analysis to assess biology impacted by cumulative monoallelic changes. SWAN enables an integrated pathway-network analysis of CNAs, RNA expression, and mutations via a simple web platform. SWAN is optimized to best prioritize known and novel tumor suppressors and oncogenes, thereby identifying drivers and potential druggable vulnerabilities within cancer CNAs. Protein homeostasis, phospholipid dephosphorylation, and ion transport pathways are commonly suppressed. An atlas of CNA pathways altered in each cancer type is released. These CNA network shifts highlight new, attractive targets to exploit in solid tumors.

Molecular Subclassification Based on Crosstalk Analysis Improves Prediction of Prognosis in Colorectal Cancer

The poor performance of single-gene lists for prognostic predictions in independent cohorts has limited their clinical use. Here, we employed a pathway-based approach using embedded biological features to identify reproducible prognostic markers as an alternative. We used pathway activity score, sure independence screening, and K-means clustering analyses to identify and cluster colorectal cancer patients into two distinct subgroups, G2 (aggressive) and G1 (moderate). The differences between these two groups with respect to survival, somatic mutation, pathway activity, and tumor-infiltration by immunocytes were compared. These comparisons revealed that the survival rates in the G2 subgroup were significantly reduced compared to that in the G1 subgroup; further, the mutational burden rates in several oncogenes, including KRAS, DCLK1, and EPHA5, were significantly higher in the G2 subgroup than in the G1 subgroup. The enhanced activity of the critical pathways such as MYC and epithelial-mesenchymal transition may also lead to the progression of colorectal cancer. Taken together, we established a novel prognostic classification system that offers meritorious insights into the hallmarks of colorectal cancer.