Engineered matrices reveal stiffness-mediated chemoresistance in patient-derived pancreatic cancer organoids.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by its fibrotic and stiff extracellular matrix. However, how the altered cell/extracellular-matrix signalling contributes to the PDAC tumour phenotype has been difficult to dissect. Here we design and engineer matrices that recapitulate the key hallmarks of the PDAC tumour extracellular matrix to address this knowledge gap. We show that patient-derived PDAC organoids from three patients develop resistance to several clinically relevant chemotherapies when cultured within high-stiffness matrices mechanically matched to in vivo tumours. Using genetic barcoding, we find that while matrix-specific clonal selection occurs, cellular heterogeneity is not the main driver of chemoresistance. Instead, matrix-induced chemoresistance occurs within a stiff environment due to the increased expression of drug efflux transporters mediated by CD44 receptor interactions with hyaluronan. Moreover, PDAC chemoresistance is reversible following transfer from high- to low-stiffness matrices, suggesting that targeting the fibrotic extracellular matrix may sensitize chemoresistant tumours. Overall, our findings support the potential of engineered matrices and patient-derived organoids for elucidating extracellular matrix contributions to human disease pathophysiology.

Germline-mediated immunoediting sculpts breast cancer subtypes and metastatic proclivity.

Tumors with the same diagnosis can have different molecular profiles and response to treatment. It remains unclear when and why these differences arise. Somatic genomic aberrations occur within the context of a highly variable germline genome. Interrogating 5870 breast cancer lesions, we demonstrated that germline-derived epitopes in recurrently amplified genes influence somatic evolution by mediating immunoediting. Individuals with a high germline-epitope burden in human epidermal growth factor receptor 2 (HER2/ERBB2) are less likely to develop HER2-positive breast cancer compared with other subtypes. The same holds true for recurrent amplicons defining three aggressive estrogen receptor (ER)-positive subgroups. Tumors that overcome such immune-mediated negative selection are more aggressive and demonstrate an "immune cold" phenotype. These data show that the germline genome plays a role in dictating somatic evolution.

Osteosarcoma PDX-Derived Cell Line Models for Preclinical Drug Evaluation Demonstrate Metastasis Inhibition by Dinaciclib through a Genome-Targeted Approach.

PURPOSE: Models to study metastatic disease in rare cancers are needed to advance preclinical therapeutics and to gain insight into disease biology. Osteosarcoma is a rare cancer with a complex genomic landscape in which outcomes for patients with metastatic disease are poor. As osteosarcoma genomes are highly heterogeneous, multiple models are needed to fully elucidate key aspects of disease biology and to recapitulate clinically relevant phenotypes. EXPERIMENTAL DESIGN: Matched patient samples, patient-derived xenografts (PDX), and PDX-derived cell lines were comprehensively evaluated using whole-genome sequencing and RNA sequencing. The in vivo metastatic phenotype of the PDX-derived cell lines was characterized in both an intravenous and an orthotopic murine model. As a proof-of-concept study, we tested the preclinical effectiveness of a cyclin-dependent kinase inhibitor on the growth of metastatic tumors in an orthotopic amputation model. RESULTS: PDXs and PDX-derived cell lines largely maintained the expression profiles of the patient from which they were derived despite the emergence of whole-genome duplication in a subset of cell lines. The cell lines were heterogeneous in their metastatic capacity, and heterogeneous tissue tropism was observed in both intravenous and orthotopic models. Single-agent dinaciclib was effective at dramatically reducing the metastatic burden. CONCLUSIONS: The variation in metastasis predilection sites between osteosarcoma PDX-derived cell lines demonstrates their ability to recapitulate the spectrum of the disease observed in patients. We describe here a panel of new osteosarcoma PDX-derived cell lines that we believe will be of wide use to the osteosarcoma research community.

The G1/S transition is promoted by Rb degradation via the E3 ligase UBR5.

Mammalian cells make the decision to divide at the G1/S transition in response to diverse signals impinging on the retinoblastoma protein Rb, a cell cycle inhibitor and tumor suppressor. Rb is inhibited by two parallel pathways. In the canonical pathway, Cyclin D-Cdk4/6 kinase complexes phosphorylate and inactivate Rb. In the second, recently discovered pathway, Rb concentration decreases during G1 to promote cells progressing through the G1/S transition. However, the mechanisms underlying this second pathway are unknown. Here, we found that the Rb concentration drop in G1 and recovery in S/G2 is controlled by phosphorylation-dependent protein degradation. In early G1 phase, un- and hypo-phosphorylated Rb is targeted by the E3 ligase UBR5. UBR5 knockout cells have higher Rb concentrations in early G1, exhibit a lower G1/S transition rate, and are more sensitive to Cdk4/6 inhibition. This last observation suggests that UBR5 inhibition can strengthen the efficacy of Cdk4/6 inhibitor-based cancer therapies.

Functional screening of amplification outlier oncogenes in organoid models of early tumorigenesis.

Somatic copy number gains are pervasive across cancer types, yet their roles in oncogenesis are insufficiently evaluated. This inadequacy is partly due to copy gains spanning large chromosomal regions, obscuring causal loci. Here, we employed organoid modeling to evaluate candidate oncogenic loci identified via integrative computational analysis of extreme copy gains overlapping with extreme expression dysregulation in The Cancer Genome Atlas. Subsets of "outlier" candidates were contextually screened as tissue-specific cDNA lentiviral libraries within cognate esophagus, oral cavity, colon, stomach, pancreas, and lung organoids bearing initial oncogenic mutations. Iterative analysis nominated the kinase DYRK2 at 12q15 as an amplified head and neck squamous carcinoma oncogene in p53-/- oral mucosal organoids. Similarly, FGF3, amplified at 11q13 in 41% of esophageal squamous carcinomas, promoted p53-/- esophageal organoid growth reversible by small molecule and soluble receptor antagonism of FGFRs. Our studies establish organoid-based contextual screening of candidate genomic drivers, enabling functional evaluation during early tumorigenesis.

Serine starvation silences estrogen receptor signaling through histone hypoacetylation.

Loss of estrogen receptor (ER) pathway activity promotes breast cancer progression, yet how this occurs remains poorly understood. Here, we show that serine starvation, a metabolic stress often found in breast cancer, represses estrogen receptor alpha (ERα) signaling by reprogramming glucose metabolism and epigenetics. Using isotope tracing and time-resolved metabolomic analyses, we demonstrate that serine is required to maintain glucose flux through glycolysis and the TCA cycle to support acetyl-CoA generation for histone acetylation. Consequently, limiting serine depletes histone H3 lysine 27 acetylation (H3K27ac), particularly at the promoter region of ER pathway genes including the gene encoding ERα, ESR1. Mechanistically, serine starvation impairs acetyl-CoA-dependent gene expression by inhibiting the entry of glycolytic carbon into the TCA cycle and down-regulating the mitochondrial citrate exporter SLC25A1, a critical enzyme in the production of nucleocytosolic acetyl-CoA from glucose. Consistent with this model, total H3K27ac and ERα expression are suppressed by SLC25A1 inhibition and restored by acetate, an alternate source of acetyl-CoA, in serine-free conditions. We thus uncover an unexpected role for serine in sustaining ER signaling through the regulation of acetyl-CoA metabolism.

A microwell platform for high-throughput longitudinal phenotyping and selective retrieval of organoids.

Organoids are powerful experimental models for studying the ontogeny and progression of various diseases including cancer. Organoids are conventionally cultured in bulk using an extracellular matrix mimic. However, bulk-cultured organoids physically overlap, making it impossible to track the growth of individual organoids over time in high throughput. Moreover, local spatial variations in bulk matrix properties make it difficult to assess whether observed phenotypic heterogeneity between organoids results from intrinsic cell differences or differences in the microenvironment. Here, we developed a microwell-based method that enables high-throughput quantification of image-based parameters for organoids grown from single cells, which can further be retrieved from their microwells for molecular profiling. Coupled with a deep learning image-processing pipeline, we characterized phenotypic traits including growth rates, cellular movement, and apical-basal polarity in two CRISPR-engineered human gastric organoid models, identifying genomic changes associated with increased growth rate and changes in accessibility and expression correlated with apical-basal polarity. A record of this paper's transparent peer review process is included in the supplemental information.

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes.

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for studying cellular differentiation, but accurately tracking cell fate transitions can be challenging, especially in disease conditions. Here we introduce PhyloVelo, a computational framework that estimates the velocity of transcriptomic dynamics by using monotonically expressed genes (MEGs) or genes with expression patterns that either increase or decrease, but do not cycle, through phylogenetic time. Through integration of scRNA-seq data with lineage information, PhyloVelo identifies MEGs and reconstructs a transcriptomic velocity field. We validate PhyloVelo using simulated data and Caenorhabditis elegans ground truth data, successfully recovering linear, bifurcated and convergent differentiations. Applying PhyloVelo to seven lineage-traced scRNA-seq datasets, generated using CRISPR-Cas9 editing, lentiviral barcoding or immune repertoire profiling, demonstrates its high accuracy and robustness in inferring complex lineage trajectories while outperforming RNA velocity. Additionally, we discovered that MEGs across tissues and organisms share similar functions in translation and ribosome biogenesis.

Deterministic evolution and stringent selection during preneoplasia.

The earliest events during human tumour initiation, although poorly characterized, may hold clues to malignancy detection and prevention1. Here we model occult preneoplasia by biallelic inactivation of TP53, a common early event in gastric cancer, in human gastric organoids. Causal relationships between this initiating genetic lesion and resulting phenotypes were established using experimental evolution in multiple clonally derived cultures over 2 years. TP53 loss elicited progressive aneuploidy, including copy number alterations and structural variants prevalent in gastric cancers, with evident preferred orders. Longitudinal single-cell sequencing of TP53-deficient gastric organoids similarly indicates progression towards malignant transcriptional programmes. Moreover, high-throughput lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programmes repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a marked degree of phenotypic convergence in premalignant epithelial organoids. These data imply predictability in the earliest stages of tumorigenesis and show evolutionary constraints and barriers to malignant transformation, with implications for earlier detection and interception of aggressive, genome-instable tumours.

Germline-mediated immunoediting sculpts breast cancer subtypes and metastatic proclivity.

Cancer represents a broad spectrum of molecularly and morphologically diverse diseases. Individuals with the same clinical diagnosis can have tumors with drastically different molecular profiles and clinical response to treatment. It remains unclear when these differences arise during disease course and why some tumors are addicted to one oncogenic pathway over another. Somatic genomic aberrations occur within the context of an individual's germline genome, which can vary across millions of polymorphic sites. An open question is whether germline differences influence somatic tumor evolution. Interrogating 3,855 breast cancer lesions, spanning pre-invasive to metastatic disease, we demonstrate that germline variants in highly expressed and amplified genes influence somatic evolution by modulating immunoediting at early stages of tumor development. Specifically, we show that the burden of germline-derived epitopes in recurrently amplified genes selects against somatic gene amplification in breast cancer. For example, individuals with a high burden of germline-derived epitopes in ERBB2, encoding human epidermal growth factor receptor 2 (HER2), are significantly less likely to develop HER2-positive breast cancer compared to other subtypes. The same holds true for recurrent amplicons that define four subgroups of ER-positive breast cancers at high risk of distant relapse. High epitope burden in these recurrently amplified regions is associated with decreased likelihood of developing high risk ER-positive cancer. Tumors that overcome such immune-mediated negative selection are more aggressive and demonstrate an "immune cold" phenotype. These data show the germline genome plays a previously unappreciated role in dictating somatic evolution. Exploiting germline-mediated immunoediting may inform the development of biomarkers that refine risk stratification within breast cancer subtypes.

Development and characterization of new patient-derived xenograft (PDX) models of osteosarcoma with distinct metastatic capacities.

Models to study metastatic disease in rare cancers are needed to advance preclinical therapeutics and to gain insight into disease biology, especially for highly aggressive cancers with a propensity for metastatic spread. Osteosarcoma is a rare cancer with a complex genomic landscape in which outcomes for patients with metastatic disease are poor. As osteosarcoma genomes are highly heterogeneous, a large panel of models is needed to fully elucidate key aspects of disease biology and to recapitulate clinically-relevant phenotypes. We describe the development and characterization of osteosarcoma patient-derived xenografts (PDXs) and a panel of PDX-derived cell lines. Matched patient samples, PDXs, and PDX-derived cell lines were comprehensively evaluated using whole genome sequencing and RNA sequencing. PDXs and PDX-derived cell lines largely maintained the expression profiles of the patient from which they were derived despite the emergence of whole-genome duplication (WGD) in a subset of cell lines. These cell line models were heterogeneous in their metastatic capacity and their tissue tropism as observed in both intravenous and orthotopic models. As proof-of-concept study, we used one of these models to test the preclinical effectiveness of a CDK inhibitor on the growth of metastatic tumors in an orthotopic amputation model. Single-agent dinaciclib was effective at dramatically reducing the metastatic burden in this model.

Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts.

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

Targeted profiling of human extrachromosomal DNA by CRISPR-CATCH.

Extrachromosomal DNA (ecDNA) is a common mode of oncogene amplification but is challenging to analyze. Here, we adapt CRISPR-CATCH, in vitro CRISPR-Cas9 treatment and pulsed field gel electrophoresis of agarose-entrapped genomic DNA, previously developed for bacterial chromosome segments, to isolate megabase-sized human ecDNAs. We demonstrate strong enrichment of ecDNA molecules containing EGFR, FGFR2 and MYC from human cancer cells and NRAS ecDNA from human metastatic melanoma with acquired therapeutic resistance. Targeted enrichment of ecDNA versus chromosomal DNA enabled phasing of genetic variants, identified the presence of an EGFRvIII mutation exclusively on ecDNAs and supported an excision model of ecDNA genesis in a glioblastoma model. CRISPR-CATCH followed by nanopore sequencing enabled single-molecule ecDNA methylation profiling and revealed hypomethylation of the EGFR promoter on ecDNAs. We distinguished heterogeneous ecDNA species within the same sample by size and sequence with base-pair resolution and discovered functionally specialized ecDNAs that amplify select enhancers or oncogene-coding sequences.

ZFP281 drives a mesenchymal-like dormancy program in early disseminated breast cancer cells that prevents metastatic outgrowth in the lung.

Increasing evidence shows that cancer cells can disseminate from early evolved primary lesions much earlier than the classical metastasis models predicted. Here, we reveal at a single-cell resolution that mesenchymal-like (M-like) and pluripotency-like programs coordinate dissemination and a long-lived dormancy program of early disseminated cancer cells (DCCs). The transcription factor ZFP281 induces a permissive state for heterogeneous M-like transcriptional programs, which associate with a dormancy signature and phenotype in vivo. Downregulation of ZFP281 leads to a loss of an invasive, M-like dormancy phenotype and a switch to lung metastatic outgrowth. We also show that FGF2 and TWIST1 induce ZFP281 expression to induce the M-like state, which is linked to CDH1 downregulation and upregulation of CDH11. We found that ZFP281 not only controls the early dissemination of cancer cells but also locks early DCCs in a dormant state by preventing the acquisition of an epithelial-like proliferative program and consequent metastases outgrowth.

Most cancers carry a substantial deleterious load due to Hill-Robertson interference.

Cancer genomes exhibit surprisingly weak signatures of negative selection (Martincorena et al., 2017; Weghorn, 2017). This may be because selective pressures are relaxed or because genome-wide linkage prevents deleterious mutations from being removed (Hill-Robertson interference; Hill and Robertson, 1966). By stratifying tumors by their genome-wide mutational burden, we observe negative selection (dN/dS ~ 0.56) in low mutational burden tumors, while remaining cancers exhibit dN/dS ratios ~1. This suggests that most tumors do not remove deleterious passengers. To buffer against deleterious passengers, tumors upregulate heat shock pathways as their mutational burden increases. Finally, evolutionary modeling finds that Hill-Robertson interference alone can reproduce patterns of attenuated selection and estimates the total fitness cost of passengers to be 46% per cell on average. Collectively, our findings suggest that the lack of observed negative selection in most tumors is not due to relaxed selective pressures, but rather the inability of selection to remove deleterious mutations in the presence of genome-wide linkage.

"Fateful" encounter: Lineage tracing meets phylogeny to unravel mysteries of cancer progression.

In a recent issue of Cell, Yang et al. utilize evolvable cellular barcodes to investigate the evolutionary trajectories of murine lung adenocarcinoma. Reconstructing the life histories of these tumors based on cellular barcodes reveals stringent selection and phenotypic differences across subclonal lineages.

Combinatorial immunotherapies overcome MYC-driven immune evasion in triple negative breast cancer.

Few patients with triple negative breast cancer (TNBC) benefit from immune checkpoint inhibitors with complete and durable remissions being quite rare. Oncogenes can regulate tumor immune infiltration, however whether oncogenes dictate diminished response to immunotherapy and whether these effects are reversible remains poorly understood. Here, we report that TNBCs with elevated MYC expression are resistant to immune checkpoint inhibitor therapy. Using mouse models and patient data, we show that MYC signaling is associated with low tumor cell PD-L1, low overall immune cell infiltration, and low tumor cell MHC-I expression. Restoring interferon signaling in the tumor increases MHC-I expression. By combining a TLR9 agonist and an agonistic antibody against OX40 with anti-PD-L1, mice experience tumor regression and are protected from new TNBC tumor outgrowth. Our findings demonstrate that MYC-dependent immune evasion is reversible and druggable, and when strategically targeted, may improve outcomes for patients treated with immune checkpoint inhibitors.

Single-cell analyses define a continuum of cell state and composition changes in the malignant transformation of polyps to colorectal cancer.

To chart cell composition and cell state changes that occur during the transformation of healthy colon to precancerous adenomas to colorectal cancer (CRC), we generated single-cell chromatin accessibility profiles and single-cell transcriptomes from 1,000 to 10,000 cells per sample for 48 polyps, 27 normal tissues and 6 CRCs collected from patients with or without germline APC mutations. A large fraction of polyp and CRC cells exhibit a stem-like phenotype, and we define a continuum of epigenetic and transcriptional changes occurring in these stem-like cells as they progress from homeostasis to CRC. Advanced polyps contain increasing numbers of stem-like cells, regulatory T cells and a subtype of pre-cancer-associated fibroblasts. In the cancerous state, we observe T cell exhaustion, RUNX1-regulated cancer-associated fibroblasts and increasing accessibility associated with HNF4A motifs in epithelia. DNA methylation changes in sporadic CRC are strongly anti-correlated with accessibility changes along this continuum, further identifying regulatory markers for molecular staging of polyps.