A machine learning platform to optimize the translation of personalized network models to the clinic

PURPOSE Dynamic network models predict clinical prognosis and inform therapeutic intervention by elucidating disease-driven aberrations at the systems level. However, the personalization of model predictions requires the profiling of multiple model inputs, which hampers clinical translation. PATIENTS AND METHODS We applied APOPTO-CELL, a prognostic model of apoptosis signaling, to showcase the establishment of computational platforms that require a reduced set of inputs. We designed two distinct and complementary pipelines: a probabilistic approach to exploit a consistent subpanel of inputs across the whole cohort (Ensemble) and a machine learning approach to identify a reduced protein set tailored for individual patients (Tree). Development was performed on a virtual cohort of 3,200,000 patients, with inputs estimated from clinically relevant protein profiles. Validation was carried out in an in-house stage III colorectal cancer cohort, with inputs profiled in surgical resections by reverse phase protein array (n = 120) and/or immunohistochemistry (n = 117). RESULTS Ensemble and Tree reproduced APOPTO-CELL predictions in the virtual patient cohort with 92% and 99% accuracy while decreasing the number of inputs to a consistent subset of three proteins (40% reduction) or a personalized subset of 2.7 proteins on average (46% reduction), respectively. Ensemble and Tree retained prognostic utility in the in-house colorectal cancer cohort. The association between the Ensemble accuracy and prognostic value (Spearman ρ = 0.43; P = .02) provided a rationale to optimize the input composition for specific clinical settings. Comparison between profiling by reverse phase protein array (gold standard) and immunohistochemistry (clinical routine) revealed that the latter is a suitable technology to quantify model inputs. CONCLUSION This study provides a generalizable framework to optimize the development of network-based prognostic assays and, ultimately, to facilitate their integration in the routine clinical workflow

Genomic Exploration of Distinct Molecular Phenotypes Steering Temozolomide Resistance Development in Patient-Derived Glioblastoma Cells

Chemotherapy using temozolomide is the standard treatment for patients with glioblastoma. Despite treatment, prognosis is still poor largely due to the emergence of temozolomide resistance. This resistance is closely linked to the widely recognized inter- and intra-tumoral heterogeneity in glioblastoma, although the underlying mechanisms are not yet fully understood. To induce temozolomide resistance, we subjected 21 patient-derived glioblastoma cell cultures to Temozolomide treatment for a period of up to 90 days. Prior to treatment, the cells’ molecular characteristics were analyzed using bulk RNA sequencing. Additionally, we performed single-cell RNA sequencing on four of the cell cultures to track the evolution of temozolomide resistance. The induced temozolomide resistance was associated with two distinct phenotypic behaviors, classified as “adaptive” (ADA) or “non-adaptive” (N-ADA) to temozolomide. The ADA phenotype displayed neurodevelopmental and metabolic gene signatures, whereas the N-ADA phenotype expressed genes related to cell cycle regulation, DNA repair, and protein synthesis. Single-cell RNA sequencing revealed that in ADA cell cultures, one or more subpopulations emerged as dominant in the resistant samples, whereas N-ADA cell cultures remained relatively stable. The adaptability and heterogeneity of glioblastoma cells play pivotal roles in temozolomide treatment and contribute to the tumor’s ability to survive. Depending on the tumor’s adaptability potential, subpopulations with acquired resistance mechanisms may arise.</p

Ex vivo drug sensitivity screening predicts response to temozolomide in glioblastoma patients and identifies candidate biomarkers

Background: Patient-derived glioma stem-like cells (GSCs) have become the gold-standard in neuro-oncological research; however, it remains to be established whether loss of in situ microenvironment affects the clinically-predictive value of this model. We implemented a GSC monolayer system to investigate in situ-in vitro molecular correspondence and the relationship between in vitro and patient response to temozolomide (TMZ). Methods: DNA/RNA-sequencing was performed on 56 glioblastoma tissues and 19 derived GSC cultures. Sensitivity to TMZ was screened across 66 GSC cultures. Viability readouts were related to clinical parameters of corresponding patients and whole-transcriptome data. Results: Tumour DNA and RNA sequences revealed strong similarity to corresponding GSCs despite loss of neuronal and immune interactions. In vitro TMZ screening yielded three response categories which significantly correlated with patient survival, therewith providing more specific prediction than the binary MGMT marker. Transcriptome analysis identified 121 genes related to TMZ sensitivity of which 21were validated in external datasets. Conclusion:GSCs retain patient-unique hallmark gene expressions despite loss of their natural environment. Drug screening using GSCs predicted patient response to TMZ more specifically than MGMT status, while transcriptome analysis identified potential biomarkers for this response. GSC drug screening therefore provides a tool to improve drug development and precision medicine for glioblastoma.</p

Search for dark matter produced in association with bottom or top quarks in √s = 13 TeV pp collisions with the ATLAS detector

A search for weakly interacting massive particle dark matter produced in association with bottom or top quarks is presented. Final states containing third-generation quarks and miss- ing transverse momentum are considered. The analysis uses 36.1 fb−1 of proton–proton collision data recorded by the ATLAS experiment at √s = 13 TeV in 2015 and 2016. No significant excess of events above the estimated backgrounds is observed. The results are in- terpreted in the framework of simplified models of spin-0 dark-matter mediators. For colour- neutral spin-0 mediators produced in association with top quarks and decaying into a pair of dark-matter particles, mediator masses below 50 GeV are excluded assuming a dark-matter candidate mass of 1 GeV and unitary couplings. For scalar and pseudoscalar mediators produced in association with bottom quarks, the search sets limits on the production cross- section of 300 times the predicted rate for mediators with masses between 10 and 50 GeV and assuming a dark-matter mass of 1 GeV and unitary coupling. Constraints on colour- charged scalar simplified models are also presented. Assuming a dark-matter particle mass of 35 GeV, mediator particles with mass below 1.1 TeV are excluded for couplings yielding a dark-matter relic density consistent with measurements