No abstract

Highlights d Proteogenomic characterization reveals the functional impact of genomic alterations d Phosphoproteomics uncovers putative therapeutic targets downstream of KRAS d Multiomics links endothelial cell remodeling and glycolysis to immune exclusion d Proteomics and glycoproteomics reveal candidates for early detection or intervention

Purpose: Non-small cell lung cancer (NSCLC) is the most common cause of cancer-related deaths worldwide. There is an unmet need to develop novel clinically relevant models of NSCLC to accelerate identification of drug targets and our understanding of the disease.Experimental Design: Thirty surgically resected NSCLC primary patient tissue and 35 previously established patient-derived xenograft (PDX) models were processed for organoid culture establishment. Organoids were histologically and molecularly characterized by cytology and histology, exome sequencing, and RNA-sequencing analysis. Tumorigenicity was assessed through subcutaneous injection of organoids in NOD/SCID mice. Organoids were subjected to drug testing using EGFR, FGFR, and MEK-targeted therapies.Results: We have identified cell culture conditions favoring the establishment of short-term and long-term expansion of NSCLC organoids derived from primary lung patient and PDX tumor tissue. The NSCLC organoids recapitulated the histology of the patient and PDX tumor. They also retained tumorigenicity, as evidenced by cytologic features of malignancy, xenograft formation, preservation of mutations, copy number aberrations, and gene expression profiles between the organoid and matched parental tumor tissue by whole-exome and RNA sequencing. NSCLC organoid models also preserved the sensitivity of the matched parental tumor to targeted therapeutics, and could be used to validate or discover biomarkerdrug combinations.Conclusions: Our panel of NSCLC organoids closely recapitulates the genomics and biology of patient tumors, and is a potential platform for drug testing and biomarker validation.

TTK protein kinase (TTK), also known as Monopolar spindle 1 (MPS1), is a key regulator of the spindle assembly checkpoint (SAC), which functions to maintain genomic integrity. TTK has emerged as a promising therapeutic target in human cancers, including triple-negative breast cancer (TNBC). Several TTK inhibitors (TTKis) are being evaluated in clinical trials, and an understanding of the mechanisms mediating TTKi sensitivity and resistance could inform the successful development of this class of agents. We evaluated the cellular effects of the potent clinical TTKi CFI-402257 in TNBC models. CFI-402257 induced apoptosis and potentiated aneuploidy in TNBC lines by accelerating progression through mitosis and inducing mitotic segregation errors. We used genome-wide CRISPR/Cas9 screens in multiple TNBC cell lines to identify mechanisms of resistance to CFI-402257. Our functional genomic screens identified members of the anaphase-promoting complex/cyclosome (APC/C) complex, which promotes mitotic progression following inactivation of the SAC. Several screen candidates were validated to confer resistance to CFI-402257 and other TTKis using CRISPR/Cas9 and siRNA methods. These findings extend the observation that impairment of the APC/C enables cells to tolerate genomic instability caused by SAC inactivation, and support the notion that a measure of APC/C function could predict the response to TTK inhibition. Indeed, an APC/C gene expression signature is significantly associated with CFI-402257 response in breast and lung adenocarcinoma cell line panels. This expression signature, along with somatic alterations in genes involved in mitotic progression, represent potential biomarkers that could be evaluated in ongoing clinical trials of CFI-402257 or other TTKis.

Multi-omics experiments are increasingly commonplace in biomedical research, and add layers of complexity to experimental design, data integration, and analysis. R and Bioconductor provide a generic framework for statistical analysis and visualization, as well as specialized data classes for a variety of high-throughput data types, but methods are lacking for integrative analysis of multi-omics experiments. The MultiAssayExperiment software package, implemented in R and leveraging Bioconductor software and design principles, provides for the coordinated representation of, storage of, and operation on multiple diverse genomics data. We provide the unrestricted multiple ‘omics data for each cancer tissue in The Cancer Genome Atlas (TCGA) as ready-to-analyze MultiAssayExperiment objects, and demonstrate in these and other datasets how the software simplifies data representation, statistical analysis, and visualization. The MultiAssayExperiment Bioconductor package reduces major obstacles to efficient, scalable and reproducible statistical analysis of multi-omics data and enhances data science applications of multiple omics datasets.

BackgroundLong non-coding RNA (lncRNA) expression has been implicated in a range of molecular mechanisms that are central in cancer. However, lncRNA expression has not yet been comprehensively characterized in acute myeloid leukemia (AML). Here, we assess to what extent lncRNA expression is prognostic of AML patient overall survival (OS) and determine if there are indications of lncRNA-based molecular subtypes of AML.MethodsWe performed RNA sequencing of 274 intensively treated AML patients in a Swedish cohort and quantified lncRNA expression. Univariate and multivariate time-to-event analysis was applied to determine association between individual lncRNAs with OS. Unsupervised statistical learning was applied to ascertain if lncRNA-based molecular subtypes exist and are prognostic.ResultsThirty-three individual lncRNAs were found to be associated with OS (adjusted p value < 0.05). We established four distinct molecular subtypes based on lncRNA expression using a consensus clustering approach. LncRNA-based subtypes were found to stratify patients into groups with prognostic information (p value < 0.05). Subsequently, lncRNA expression-based subtypes were validated in an independent patient cohort (TCGA-AML). LncRNA subtypes could not be directly explained by any of the recurrent cytogenetic or mutational aberrations, although associations with some of the established genetic and clinical factors were found, including mutations in NPM1, TP53, and FLT3.ConclusionLncRNA expression-based four subtypes, discovered in this study, are reproducible and can effectively stratify AML patients. LncRNA expression profiling can provide valuable information for improved risk stratification of AML patients.Electronic supplementary materialThe online version of this article (10.1186/s13045-018-0596-2) contains supplementary material, which is available to authorized users.

Risk stratification of acute myeloid leukemia (AML) patients needs improvement. Several AML risk classification models based on somatic mutations or gene-expression profiling have been proposed. However, systematic and independent validation of these models is required for future clinical implementation. We performed whole-transcriptome RNA-sequencing and panel-based deep DNA sequencing of 23 genes in 274 intensively treated AML patients (Clinseq-AML). We also utilized the The Cancer Genome Atlas (TCGA)-AML study (N=142) as a second validation cohort. We evaluated six previously proposed molecular-based models for AML risk stratification and two revised risk classification systems combining molecular- and clinical data. Risk groups stratified by five out of six models showed different overall survival in cytogenetic normal-AML patients in the Clinseq-AML cohort (P-value<0.05; concordance index >0.5). Risk classification systems integrating mutational or gene-expression data were found to add prognostic value to the current European Leukemia Net (ELN) risk classification. The prognostic value varied between models and across cohorts, highlighting the importance of independent validation to establish evidence of efficacy and general applicability. All but one model replicated in the Clinseq-AML cohort, indicating the potential for molecular-based AML risk models. Risk classification based on a combination of molecular and clinical data holds promise for improved AML patient stratification in the future.

Resistance to targeted therapies is an important clinical problem in HER2-positive (HER2+) breast cancer. “Drug-tolerant persisters” (DTP), a subpopulation of cancer cells that survive via reversible, nongenetic mechanisms, are implicated in resistance to tyrosine kinase inhibitors (TKI) in other malignancies, but DTPs following HER2 TKI exposure have not been well characterized. We found that HER2 TKIs evoke DTPs with a luminal-like or a mesenchymal-like transcriptome. Lentiviral barcoding/single-cell RNA sequencing reveals that HER2+ breast cancer cells cycle stochastically through a “pre-DTP” state, characterized by a G0-like expression signature and enriched for diapause and/or senescence genes. Trajectory analysis/cell sorting shows that pre-DTPs preferentially yield DTPs upon HER2 TKI exposure. Cells with similar transcriptomes are present in HER2+ breast tumors and are associated with poor TKI response. Finally, biochemical experiments indicate that luminal-like DTPs survive via estrogen receptor–dependent induction of SGK3, leading to rewiring of the PI3K/AKT/mTORC1 pathway to enable AKT-independent mTORC1 activation. Significance: DTPs are implicated in resistance to anticancer therapies, but their ontogeny and vulnerabilities remain unclear. We find that HER2 TKI-DTPs emerge from stochastically arising primed cells (“pre-DTPs”) that engage either of two distinct transcriptional programs upon TKI exposure. Our results provide new insights into DTP ontogeny and potential therapeutic vulnerabilities. This article is highlighted in the In This Issue feature, p. 873