Quantitative phosphoproteomics to unravel the cellular response to chemical stressors with different modes of action.

Damage to cellular macromolecules and organelles by chemical exposure evokes activation of various stress response pathways. To what extent different chemical stressors activate common and stressor-specific pathways is largely unknown. Here, we used quantitative phosphoproteomics to compare the signaling events induced by four stressors with different modes of action: the DNA damaging agent: cisplatin (CDDP), the topoisomerase II inhibitor: etoposide (ETO), the pro-oxidant: diethyl maleate (DEM) and the immunosuppressant: cyclosporine A (CsA) administered at an equitoxic dose to mouse embryonic stem cells. We observed major differences between the stressors in the number and identity of responsive phosphosites and the amplitude of phosphorylation. Kinase motif and pathway analyses indicated that the DNA damage response (DDR) activation by CDDP occurs predominantly through the replication-stress-related Atr kinase, whereas ETO triggers the DDR through Atr as well as the DNA double-strand-break-associated Atm kinase. CsA shares with ETO activation of CK2 kinase. Congruent with their known modes of action, CsA-mediated signaling is related to down-regulation of pathways that control hematopoietic differentiation and immunity, whereas oxidative stress is the most prominent initiator of DEM-modulated stress signaling. This study shows that even at equitoxic doses, different stressors induce distinctive and complex phosphorylation signaling cascades.

Vesicle-bound EBV-BART13-3p miRNA in circulation distinguishes nasopharyngeal from other head and neck cancer and asymptomatic EBV-infections.

Cell-free microRNA (miRNA) in biofluids released by tumors in either protein or vesicle-bound form, represent promising minimally-invasive cancer biomarkers. However, a highly abundant non-tumor background in human plasma and serum complicates the discovery and detection of tumor-selective circulating miRNAs. We performed small RNA sequencing on serum and plasma RNA from Nasopharyngeal Carcinoma (NPC) patients. Collectively, Epstein Barr virus-encoded miRNAs, more so than endogenous miRNAs, signify presence of NPC. However, RNAseq-based EBV miRNA profiles differ between NPC patients, suggesting inter-tumor heterogeneity or divergent secretory characteristics. We determined with sensitive qRT-PCR assays that EBV miRNAs BART7-3p, BART9-3p and BART13-3p are actively secreted by C666.1 NPC cells bound to extracellular vesicles (EVs) and soluble ribonucleoprotein complexes. Importantly, these miRNAs are expressed in all primary NPC tumor biopsies and readily detected in nasopharyngeal brushings from both early and late-stage NPC patients. Increased levels of BART7-3p, BART9-3p and particularly BART13-3p, distinguish NPC patient sera from healthy controls. Receiver operating characteristic curve analysis using sera from endemic NPC patients, other head and neck cancers and individuals with asymptomatic EBV-infections reveals a superior diagnostic performance of EBV miRNAs over anti-EBNA1 IgA serology and EBV-DNA load (AUC 0.87-0.96 vs 0.86 and 0.66 respectively). The high specificity of circulating EBV-BART13-3p (97%) for NPC detection is in agreement with active secretion from NPC tumor cells. We conclude EV-bound BART13-3p in circulation is a promising, NPC-selective, biomarker that should be considered as part of a screening strategy to identify NPC in endemic regions.

The functional impact of variants of uncertain significance in BRCA2.

PURPOSE: Genetic testing has uncovered large numbers of variants in the BRCA2 gene for which the clinical significance is unclear. Cancer risk prediction of these variants of uncertain significance (VUS) can be improved by reliable assessment of the extent of impairment of the tumor suppressor function(s) of BRCA2. METHODS: Here, we evaluated the performance of the mouse embryonic stem cell (mESC)-based functional assay on an extensive set of BRCA2 missense variants. RESULTS: Whereas all 20 nonpathogenic (class 1/2) variants were able to complement the cell lethal phenotype induced by loss of endogenous mouse Brca2, only 1 out of 15 pathogenic (class 4/5) variants (p.Gly2609Asp) was able to do so. However, in this variant the major tumor suppressive activity of BRCA2, i.e., homology directed repair (HDR), was severely abrogated. Among 43 evaluated VUS (class 3), 7 were unable to complement the lethal phenotype of mouse Brca2 loss while 7 other variants displayed a more severe reduction of HDR activity than observed for class 1/ 2 variants. CONCLUSION: The mESC-based BRCA2 functional assay can reliably determine the functional impact of VUS, distinguish between pathogenic and nonpathogenic variants, and may contribute to improved cancer risk estimation for BRCA2 VUS carriers.

Deciphering the RNA landscape by RNAome sequencing.

Current RNA expression profiling methods rely on enrichment steps for specific RNA classes, thereby not detecting all RNA species in an unperturbed manner. We report strand-specific RNAome sequencing that determines expression of small and large RNAs from rRNA-depleted total RNA in a single sequence run. Since current analysis pipelines cannot reliably analyze small and large RNAs simultaneously, we developed TRAP, Total Rna Analysis Pipeline, a robust interface that is also compatible with existing RNA sequencing protocols. RNAome sequencing quantitatively preserved all RNA classes, allowing cross-class comparisons that facilitates the identification of relationships between different RNA classes. We demonstrate the strength of RNAome sequencing in mouse embryonic stem cells treated with cisplatin. MicroRNA and mRNA expression in RNAome sequencing significantly correlated between replicates and was in concordance with both existing RNA sequencing methods and gene expression arrays generated from the same samples. Moreover, RNAome sequencing also detected additional RNA classes such as enhancer RNAs, anti-sense RNAs, novel RNA species and numerous differentially expressed RNAs undetectable by other methods. At the level of complete RNA classes, RNAome sequencing also identified a specific global repression of the microRNA and microRNA isoform classes after cisplatin treatment whereas all other classes such as mRNAs were unchanged. These characteristics of RNAome sequencing will significantly improve expression analysis as well as studies on RNA biology not covered by existing methods.

Divergent Molecular and Cellular Responses to Low and High-Dose Ionizing Radiation.

Cancer risk after ionizing radiation (IR) is assumed to be linear with the dose; however, for low doses, definite evidence is lacking. Here, using temporal multi-omic systems analyses after a low (LD; 0.1 Gy) or a high (HD; 1 Gy) dose of X-rays, we show that, although the DNA damage response (DDR) displayed dose proportionality, many other molecular and cellular responses did not. Phosphoproteomics uncovered a novel mode of phospho-signaling via S12-PPP1R7, and large-scale dephosphorylation events that regulate mitotic exit control in undamaged cells and the G2/M checkpoint upon IR in a dose-dependent manner. The phosphoproteomics of irradiated DNA double-strand breaks (DSBs) repair-deficient cells unveiled extended phospho-signaling duration in either a dose-dependent (DDR signaling) or independent (mTOR-ERK-MAPK signaling) manner without affecting signal magnitude. Nascent transcriptomics revealed the transcriptional activation of genes involved in NRF2-regulated antioxidant defense, redox-sensitive ERK-MAPK signaling, glycolysis and mitochondrial function after LD, suggesting a prominent role for reactive oxygen species (ROS) in molecular and cellular responses to LD exposure, whereas DDR genes were prominently activated after HD. However, how and to what extent the observed dose-dependent differences in molecular and cellular responses may impact cancer development remain unclear, as the induction of chromosomal damage was found to be dose-proportional (10-200 mGy).

The Extended ToxTracker Assay Discriminates Between Induction of DNA Damage, Oxidative Stress, and Protein Misfolding.

Chemical exposure of cells may damage biomolecules, cellular structures, and organelles thereby jeopardizing cellular homeostasis. A multitude of defense mechanisms have evolved that can recognize specific types of damaged molecules and will initiate distinct cellular programs aiming to remove the damage inflicted and prevent cellular havoc. As a consequence, quantitative assessment of the activity of the cellular stress responses may serve as a sensitive reporter for the induction of specific types of damage. We have previously developed the ToxTracker assay, a mammalian stem cell-based genotoxicity assay employing two green fluorescent protein reporters specific for DNA damage and oxidative stress. We have now expanded the ToxTracker assay with an additional four reporter cell lines to include monitoring of additional stress signaling pathways. This panel of six green fluorescent protein reporters is able to discriminate between different primary reactivity of chemicals being their ability to react with DNA and block DNA replication, induce oxidative stress, activate the unfolded protein response, or cause a general P53-dependent cellular stress response. Extensive validation using the compound library suggested by the European Centre for the Validation of Alternative Methods (ECVAM) and a large panel of reference chemicals shows that the ToxTracker assay has an outstanding sensitivity and specificity. In addition, we developed Toxplot, a dedicated software tool for automated data analysis and graphical representation of the test results. Rapid and reliable identification by the ToxTracker assay of specific biological reactivity can significantly improve in vitro human hazard assessment of chemicals.

Non­invasive prostate cancer detection by measuring miRNA variants (isomiRs) in urine extracellular vesicles.

In many cancer types, the expression and function of ~22 nucleotide-long microRNAs (miRNA) is deregulated. Mature miRNAs can be stably detected in extracellular vesicles (EVs) in biofluids, therefore they are considered to have great potential as biomarkers. In the present study, we investigated whether miRNAs have a distinct expression pattern in urine-EVs of prostate cancer (PCa) patients compared to control males. By next generation sequencing, we determined the miRNA expression in a discovery cohort of 4 control men and 9 PCa patients. miRNAs were validated by using a stemloop RT-PCR in an independent cohort of 74 patients (26 control and 48 PCa-patients). Whereas standard mapping protocols identified > 10 PCa associated miRNAs in urinary EVs, miR-21, miR-375 and miR-204 failed to robustly discriminate for disease in a validation study with RT-PCR-detection of mature miRNA sequences. In contrast, we observed that miRNA isoforms (isomiRs) with 3' end modifications were highly discriminatory between samples from control men and PCa patients. Highly differentially expressed isomiRs of miR-21, miR-204 and miR-375 were subsequently validated in an independent group of 74 patients. Receiver-operating characteristic analysis was performed to evaluate the diagnostic performance of three isomiRs, resulting in a 72.9% sensitivity with a high (88%) specificity and an area under the curve (AUC) of 0.866. In comparison, prostate specific antigen had an AUC of 0.707 and measuring the mature form of these miRNAs yielded a lower 70.8% sensitivity and 72% specificity (AUC 0.766). We propose that isomiRs may carry discriminatory information which is useful to generate stronger biomarkers.

Global phosphoproteome profiling reveals unanticipated networks responsive to cisplatin treatment of embryonic stem cells.

Cellular responses to DNA-damaging agents involve the activation of various DNA damage signaling and transduction pathways. Using quantitative and high-resolution tandem mass spectrometry, we determined global changes in protein level and phosphorylation site profiles following treatment of SILAC (stable isotope labeling by amino acids in cell culture)-labeled murine embryonic stem cells with the anticancer drug cisplatin. Network and pathway analyses indicated that processes related to the DNA damage response and cytoskeleton organization were significantly affected. Although the ATM (ataxia telangiectasia mutated) and ATR (ATM and Rad3-related) consensus sequence (S/T-Q motif) was significantly overrepresented among hyperphosphorylated peptides, about half of the >2-fold-upregulated phosphorylation sites based on the consensus sequence were not direct substrates of ATM and ATR. Eleven protein kinases mainly belonging to the mitogen-activated protein kinase (MAPK) family were identified as being regulated in their kinase domain activation loop. The biological importance of three of these kinases (cyclin-dependent kinase 7 [CDK7], Plk1, and KPCD1) in the protection against cisplatin-induced cytotoxicity was demonstrated by small interfering RNA (siRNA)-mediated knockdown. Our results indicate that the cellular response to cisplatin involves a variety of kinases and phosphatases not only acting in the nucleus but also regulating cytoplasmic targets, resulting in extensive cytoskeletal rearrangements. Integration of transcriptomic and proteomic data revealed a poor correlation between changes in the relative levels of transcripts and their corresponding proteins, but a large overlap in affected pathways at the levels of mRNA, protein, and phosphoprotein. This study provides an integrated view of pathways activated by genotoxic stress and deciphers kinases that play a pivotal role in regulating cellular processes other than the DNA damage response.