Background

Applying directed acyclic graph (DAG) models to proteogenomic data has been shown effective for detecting causal biomarkers of complex diseases. However, there remain unsolved challenges in DAG learning to jointly model binary clinical outcome variables and continuous biomarker measurements.

Results

In this paper, we propose a new tool, DAGBagM, to learn DAGs with both continuous and binary nodes. By using appropriate models, DAGBagM allows for either continuous or binary nodes to be parent or child nodes. It employs a bootstrap aggregating strategy to reduce false positives in edge inference. At the same time, the aggregation procedure provides a flexible framework to robustly incorporate prior information on edges.

Conclusions

Through extensive simulation experiments, we demonstrate that DAGBagM has superior performance compared to alternative strategies for modeling mixed types of nodes. In addition, DAGBagM is computationally more efficient than two competing methods. When applying DAGBagM to proteogenomic datasets from ovarian cancer studies, we identify potential protein biomarkers for platinum refractory/resistant response in ovarian cancer. DAGBagM is made available as a github repository at https://github.com/jie108/dagbagM .

Among the pleotropic roles of transforming growth factor-β (TGFβ) signaling in cancer, its impact on genomic stability is least understood. Inhibition of TGFβ signaling increases use of alternative end joining (alt-EJ), an error-prone DNA repair process that typically functions as a "backup" pathway if double-strand break repair by homologous recombination or nonhomologous end joining is compromised. However, the consequences of this functional relationship on therapeutic vulnerability in human cancer remain unknown. Here, we show that TGFβ broadly controls the DNA damage response and suppresses alt-EJ genes that are associated with genomic instability. Mechanistically based TGFβ and alt-EJ gene expression signatures were anticorrelated in glioblastoma, squamous cell lung cancer, and serous ovarian cancer. Consistent with error-prone repair, more of the genome was altered in tumors classified as low TGFβ and high alt-EJ, and the corresponding patients had better outcomes. Pan-cancer analysis of solid neoplasms revealed that alt-EJ genes were coordinately expressed and anticorrelated with TGFβ competency in 16 of 17 cancer types tested. Moreover, regardless of cancer type, tumors classified as low TGFβ and high alt-EJ were characterized by an insertion-deletion mutation signature containing short microhomologies and were more sensitive to genotoxic therapy. Collectively, experimental studies revealed that loss or inhibition of TGFβ signaling compromises the DNA damage response, resulting in ineffective repair by alt-EJ. Translation of this mechanistic relationship into gene expression signatures identified a robust anticorrelation that predicts response to genotoxic therapies, thereby expanding the potential therapeutic scope of TGFβ biology.

Summary

A primary goal of the US National Cancer Institute's Ras initiative at the Frederick National Laboratory for Cancer Research is to develop methods to quantify RAS signaling to facilitate development of novel cancer therapeutics. We use targeted proteomics technologies to develop a community resource consisting of 256 validated multiple reaction monitoring (MRM)-based, multiplexed assays for quantifying protein expression and phosphorylation through the receptor tyrosine kinase, MAPK, and AKT signaling networks. As proof of concept, we quantify the response of melanoma (A375 and SK-MEL-2) and colorectal cancer (HCT-116 and HT-29) cell lines to BRAF inhibition by PLX-4720. These assays replace over 60 Western blots with quantitative mass spectrometry-based assays of high molecular specificity and quantitative precision, showing the value of these methods for pharmacodynamic measurements and mechanism of action studies. Methods, fit-for-purpose validation, and results are publicly available as a resource for the community at assays.cancer.gov.

Motivation

A lack of quantitative, multiplexable assays for phosphosignaling limits comprehensive investigation of aberrant signaling in cancer and evaluation of novel treatments. To alleviate this limitation, we sought to develop assays using targeted mass spectrometry for quantifying protein expression and phosphorylation through the receptor tyrosine kinase, MAPK, and AKT signaling networks. The resulting assays provide a resource for replacing over 60 Western blots in examining cancer signaling and tumor biology with high molecular specificity and quantitative rigor.

Adoption of targeted mass spectrometry (MS) approaches such as multiple reaction monitoring (MRM) to study biological and biomedical questions is well underway in the proteomics community. Successful application depends on the ability to generate reliable assays that uniquely and confidently identify target peptides in a sample. Unfortunately, there is a wide range of criteria being applied to say that an assay has been successfully developed. There is no consensus on what criteria are acceptable and little understanding of the impact of variable criteria on the quality of the results generated. Publications describing targeted MS assays for peptides frequently do not contain sufficient information for readers to establish confidence that the tests work as intended or to be able to apply the tests described in their own labs. Guidance must be developed so that targeted MS assays with established performance can be made widely distributed and applied by many labs worldwide. To begin to address the problems and their solutions, a workshop was held at the National Institutes of Health with representatives from the multiple communities developing and employing targeted MS assays. Participants discussed the analytical goals of their experiments and the experimental evidence needed to establish that the assays they develop work as intended and are achieving the required levels of performance. Using this "fit-for-purpose" approach, the group defined three tiers of assays distinguished by their performance and extent of analytical characterization. Computational and statistical tools useful for the analysis of targeted MS results were described. Participants also detailed the information that authors need to provide in their manuscripts to enable reviewers and readers to clearly understand what procedures were performed and to evaluate the reliability of the peptide or protein quantification measurements reported. This paper presents a summary of the meeting and recommendations.

To provide a detailed analysis of the molecular components and underlying mechanisms associated with ovarian cancer, we performed a comprehensive mass-spectrometry-based proteomic characterization of 174 ovarian tumors previously analyzed by The Cancer Genome Atlas (TCGA), of which 169 were high-grade serous carcinomas (HGSCs). Integrating our proteomic measurements with the genomic data yielded a number of insights into disease, such as how different copy-number alternations influence the proteome, the proteins associated with chromosomal instability, the sets of signaling pathways that diverse genome rearrangements converge on, and the ones most associated with short overall survival. Specific protein acetylations associated with homologous recombination deficiency suggest a potential means for stratifying patients for therapy. In addition to providing a valuable resource, these findings provide a view of how the somatic genome drives the cancer proteome and associations between protein and post-translational modification levels and clinical outcomes in HGSC. VIDEO ABSTRACT.