An oncogenic phenoscape of colonic stem cell polarization.

Cancer cells are regulated by oncogenic mutations and microenvironmental signals, yet these processes are often studied separately. To functionally map how cell-intrinsic and cell-extrinsic cues co-regulate cell fate, we performed a systematic single-cell analysis of 1,107 colonic organoid cultures regulated by (1) colorectal cancer (CRC) oncogenic mutations, (2) microenvironmental fibroblasts and macrophages, (3) stromal ligands, and (4) signaling inhibitors. Multiplexed single-cell analysis revealed a stepwise epithelial differentiation phenoscape dictated by combinations of oncogenes and stromal ligands, spanning from fibroblast-induced Clusterin (CLU)+ revival colonic stem cells (revCSCs) to oncogene-driven LRIG1+ hyper-proliferative CSCs (proCSCs). The transition from revCSCs to proCSCs is regulated by decreasing WNT3A and TGF-ß-driven YAP signaling and increasing KRASG12D or stromal EGF/Epiregulin-activated MAPK/PI3K flux. We find that APC loss and KRASG12D collaboratively limit access to revCSCs and disrupt stromal-epithelial communication-trapping epithelia in the proCSC fate. These results reveal that oncogenic mutations dominate homeostatic differentiation by obstructing cell-extrinsic regulation of cell-fate plasticity.

Trellis tree-based analysis reveals stromal regulation of patient-derived organoid drug responses.

Patient-derived organoids (PDOs) can model personalized therapy responses; however, current screening technologies cannot reveal drug response mechanisms or how tumor microenvironment cells alter therapeutic performance. To address this, we developed a highly multiplexed mass cytometry platform to measure post-translational modification (PTM) signaling, DNA damage, cell-cycle activity, and apoptosis in >2,500 colorectal cancer (CRC) PDOs and cancer-associated fibroblasts (CAFs) in response to clinical therapies at single-cell resolution. To compare patient- and microenvironment-specific drug responses in thousands of single-cell datasets, we developed "Trellis"-a highly scalable, tree-based treatment effect analysis method. Trellis single-cell screening revealed that on-target cell-cycle blockage and DNA-damage drug effects are common, even in chemorefractory PDOs. However, drug-induced apoptosis is rarer, patient-specific, and aligns with cancer cell PTM signaling. We find that CAFs can regulate PDO plasticity-shifting proliferative colonic stem cells (proCSCs) to slow-cycling revival colonic stem cells (revCSCs) to protect cancer cells from chemotherapy.

Oncogenic KRAS Regulates Tumor Cell Signaling via Stromal Reciprocation.

Oncogenic mutations regulate signaling within both tumor cells and adjacent stromal cells. Here, we show that oncogenic KRAS (KRAS(G12D)) also regulates tumor cell signaling via stromal cells. By combining cell-specific proteome labeling with multivariate phosphoproteomics, we analyzed heterocellular KRAS(G12D) signaling in pancreatic ductal adenocarcinoma (PDA) cells. Tumor cell KRAS(G12D) engages heterotypic fibroblasts, which subsequently instigate reciprocal signaling in the tumor cells. Reciprocal signaling employs additional kinases and doubles the number of regulated signaling nodes from cell-autonomous KRAS(G12D). Consequently, reciprocal KRAS(G12D) produces a tumor cell phosphoproteome and total proteome that is distinct from cell-autonomous KRAS(G12D) alone. Reciprocal signaling regulates tumor cell proliferation and apoptosis and increases mitochondrial capacity via an IGF1R/AXL-AKT axis. These results demonstrate that oncogene signaling should be viewed as a heterocellular process and that our existing cell-autonomous perspective underrepresents the extent of oncogene signaling in cancer. VIDEO ABSTRACT.

Cell-type-specific signaling networks in heterocellular organoids.

Despite the widespread adoption of organoids as biomimetic tissue models, methods to comprehensively analyze cell-type-specific post-translational modification (PTM) signaling networks in organoids are absent. Here, we report multivariate single-cell analysis of such networks in organoids and organoid cocultures. Simultaneous analysis by mass cytometry of 28 PTMs in >1 million single cells derived from small intestinal organoids reveals cell-type- and cell-state-specific signaling networks in stem, Paneth, enteroendocrine, tuft and goblet cells, as well as enterocytes. Integrating single-cell PTM analysis with thiol-reactive organoid barcoding in situ (TOBis) enables high-throughput comparison of signaling networks between organoid cultures. Cell-type-specific PTM analysis of colorectal cancer organoid cocultures reveals that shApc, KrasG12D and Trp53R172H cell-autonomously mimic signaling states normally induced by stromal fibroblasts and macrophages. These results demonstrate how standard mass cytometry workflows can be modified to perform high-throughput multivariate cell-type-specific signaling analysis of healthy and cancerous organoids.

Microenvironmental Heterogeneity Parallels Breast Cancer Progression: A Histology-Genomic Integration Analysis.

BACKGROUND: The intra-tumor diversity of cancer cells is under intense investigation; however, little is known about the heterogeneity of the tumor microenvironment that is key to cancer progression and evolution. We aimed to assess the degree of microenvironmental heterogeneity in breast cancer and correlate this with genomic and clinical parameters. METHODS AND FINDINGS: We developed a quantitative measure of microenvironmental heterogeneity along three spatial dimensions (3-D) in solid tumors, termed the tumor ecosystem diversity index (EDI), using fully automated histology image analysis coupled with statistical measures commonly used in ecology. This measure was compared with disease-specific survival, key mutations, genome-wide copy number, and expression profiling data in a retrospective study of 510 breast cancer patients as a test set and 516 breast cancer patients as an independent validation set. In high-grade (grade 3) breast cancers, we uncovered a striking link between high microenvironmental heterogeneity measured by EDI and a poor prognosis that cannot be explained by tumor size, genomics, or any other data types. However, this association was not observed in low-grade (grade 1 and 2) breast cancers. The prognostic value of EDI was superior to known prognostic factors and was enhanced with the addition of TP53 mutation status (multivariate analysis test set, p = 9 × 10-4, hazard ratio = 1.47, 95% CI 1.17-1.84; validation set, p = 0.0011, hazard ratio = 1.78, 95% CI 1.26-2.52). Integration with genome-wide profiling data identified losses of specific genes on 4p14 and 5q13 that were enriched in grade 3 tumors with high microenvironmental diversity that also substratified patients into poor prognostic groups. Limitations of this study include the number of cell types included in the model, that EDI has prognostic value only in grade 3 tumors, and that our spatial heterogeneity measure was dependent on spatial scale and tumor size. CONCLUSIONS: To our knowledge, this is the first study to couple unbiased measures of microenvironmental heterogeneity with genomic alterations to predict breast cancer clinical outcome. We propose a clinically relevant role of microenvironmental heterogeneity for advanced breast tumors, and highlight that ecological statistics can be translated into medical advances for identifying a new type of biomarker and, furthermore, for understanding the synergistic interplay of microenvironmental heterogeneity with genomic alterations in cancer cells.

Cell-specific labeling enzymes for analysis of cell-cell communication in continuous co-culture.

We report the orthologous screening, engineering, and optimization of amino acid conversion enzymes for cell-specific proteomic labeling. Intracellular endoplasmic-reticulum-anchored Mycobacterium tuberculosis diaminopimelate decarboxylase (DDC(M.tub-KDEL)) confers cell-specific meso-2,6-diaminopimelate-dependent proliferation to multiple eukaryotic cell types. Optimized lysine racemase (Lyr(M37-KDEL)) supports D-lysine specific proliferation and efficient cell-specific isotopic labeling. When ectopically expressed in discrete cell types, these enzymes confer 90% cell-specific isotopic labeling efficiency after 10 days of co-culture. Moreover, DDC(M.tub-KDEL) and Lyr(M37-KDEL) facilitate equally high cell-specific labeling fidelity without daily media exchange. Consequently, the reported novel enzyme pairing can be used to study cell-specific signaling in uninterrupted, continuous co-cultures. Demonstrating the importance of increased labeling stability for addressing novel biological questions, we compare the cell-specific phosphoproteome of fibroblasts in direct co-culture with epithelial tumor cells in both interrupted (daily media exchange) and continuous (no media exchange) co-cultures. This analysis identified multiple cell-specific phosphorylation sites specifically regulated in the continuous co-culture. Given their applicability to multiple cell types, continuous co-culture labeling fidelity, and suitability for long-term cell-cell phospho-signaling experiments, we propose DDC(M.tub-KDEL) and Lyr(M37-KDEL) as excellent enzymes for cell-specific labeling with amino acid precursors.

Reproducible automated phosphopeptide enrichment using magnetic TiO2 and Ti-IMAC.

Reproducible, comprehensive phosphopeptide enrichment is essential for studying phosphorylation-regulated processes. Here, we describe the application of hyper-porous magnetic TiO2 and Ti-IMAC microspheres for uniform automated phosphopeptide enrichment. Combining magnetic microspheres with a magnetic particle-handling robot enables rapid (45 min), reproducible (r2 ≥ 0.80) and high-fidelity (>90% purity) phosphopeptide purification in a 96-well format. Automated phosphopeptide enrichment demonstrates reproducible synthetic phosphopeptide recovery across 2 orders of magnitude, "well-to-well" quantitative reproducibility indistinguishable to internal SILAC standards, and robust "plate-to-plate" reproducibility across 5 days of independent enrichments. As a result, automated phosphopeptide enrichment enables statistical analysis of label-free phosphoproteomic samples in a high-throughput manner. This technique uses commercially available, off-the-shelf components and can be easily adopted by any laboratory interested in phosphoproteomic analysis. We provide a free downloadable automated phosphopeptide enrichment program to facilitate uniform interlaboratory collaboration and exchange of phosphoproteomic data sets.

Cross-domain inhibition of TACE ectodomain.

Proteolytic release from the cell surface is an essential activation event for many growth factors and cytokines. TNF-α converting enzyme (TACE) is a membrane-bound metalloprotease responsible for solubilizing many pathologically significant membrane substrates and is an attractive therapeutic target for the treatment of cancer and arthritis. Prior attempts to antagonize cell-surface TACE activity have focused on small-molecule inhibition of the metalloprotease active site. Given the highly conserved nature of metalloprotease active sites, this paradigm has failed to produce a truly specific TACE inhibitor and continues to obstruct the clinical investigation of TACE activity. We report the bespoke development of a specific TACE inhibitory human antibody using "two-step" phage display. This approach combines calculated selection conditions with antibody variable-domain exchange to direct individual antibody variable domains to desired epitopes. The resulting "cross-domain" human antibody is a previously undescribed selective TACE antagonist and provides a unique alternative to small-molecule metalloprotease inhibition.

Plastic persisters: revival stem cells in colorectal cancer.

Colorectal cancer (CRC) is traditionally considered to be a genetically driven disease. However, nongenetic plasticity has recently emerged as a major driver of tumour initiation, metastasis, and therapy response in CRC. Central to these processes is a recently discovered cell type, the revival colonic stem cell (revCSC). In contrast to traditional proliferative CSCs (proCSCs), revCSCs prioritise survival over propagation. revCSCs play an essential role in primary tumour formation, metastatic dissemination, and nongenetic chemoresistance. Current evidence suggests that CRC tumours leverage intestinal stem cell plasticity to both proliferate (via proCSCs) when unchallenged and survive (via revCSCs) in response to cell-extrinsic pressures. Although revCSCs likely represent a major source of therapeutic failure in CRC, our increasing knowledge of this important stem cell fate provides novel opportunities for therapeutic intervention.

Functional analysis of cell plasticity using single-cell technologies.

Metazoan organisms are heterocellular systems composed of hundreds of different cell types, which arise from an isogenic genome through differentiation. Cellular 'plasticity' further enables cells to alter their fate in response to exogenous cues and is involved in a variety of processes, such as wound healing, infection, and cancer. Recent advances in cellular model systems, high-dimensional single-cell technologies, and lineage tracing have sparked a renaissance in plasticity research. Here, we discuss the definition of cell plasticity, evaluate state-of-the-art model systems and techniques to study cell-fate dynamics, and explore the application of single-cell technologies to obtain functional insights into cell plasticity in healthy and diseased tissues. The integration of advanced biomimetic model systems, single-cell technologies, and high-throughput perturbation studies is enabling a new era of research into non-genetic plasticity in metazoan systems.

The MondoA-dependent TXNIP/GDF15 axis predicts oxaliplatin response in colorectal adenocarcinomas.

Chemotherapy, the standard of care treatment for cancer patients with advanced disease, has been increasingly recognized to activate host immune responses to produce durable outcomes. Here, in colorectal adenocarcinoma (CRC) we identify oxaliplatin-induced Thioredoxin-Interacting Protein (TXNIP), a MondoA-dependent tumor suppressor gene, as a negative regulator of Growth/Differentiation Factor 15 (GDF15). GDF15 is a negative prognostic factor in CRC and promotes the differentiation of regulatory T cells (Tregs), which inhibit CD8 T-cell activation. Intriguingly, multiple models including patient-derived tumor organoids demonstrate that the loss of TXNIP and GDF15 responsiveness to oxaliplatin is associated with advanced disease or chemotherapeutic resistance, with transcriptomic or proteomic GDF15/TXNIP ratios showing potential as a prognostic biomarker. These findings illustrate a potentially common pathway where chemotherapy-induced epithelial oxidative stress drives local immune remodeling for patient benefit, with disruption of this pathway seen in refractory or advanced cases.

GEODESIC SINKHORN FOR FAST AND ACCURATE OPTIMAL TRANSPORT ON MANIFOLDS.

Efficient computation of optimal transport distance between distributions is of growing importance in data science. Sinkhorn-based methods are currently the state-of-the-art for such computations, but require On2 computations. In addition, Sinkhorn-based methods commonly use an Euclidean ground distance between datapoints. However, with the prevalence of manifold structured scientific data, it is often desirable to consider geodesic ground distance. Here, we tackle both issues by proposing Geodesic Sinkhorn-based on diffusing a heat kernel on a manifold graph. Notably, Geodesic Sinkhorn requires only O(nlog⁡n) computation, as we approximate the heat kernel with Chebyshev polynomials based on the sparse graph Laplacian. We apply our method to the computation of barycenters of several distributions of high dimensional single cell data from patient samples undergoing chemotherapy. In particular, we define the barycentric distance as the distance between two such barycenters. Using this definition, we identify an optimal transport distance and path associated with the effect of treatment on cellular data.

Highlights from the 1st European cancer dependency map symposium and workshop.

The systematic identification of tumour vulnerabilities through perturbational experiments on cancer models, including genome editing and drug screens, is playing a crucial role in combating cancer. This collective effort is known as the Cancer Dependency Map (DepMap). The 1st European Cancer Dependency Map Symposium (EuroDepMap), held in Milan last May, featured talks, a roundtable discussion, and a poster session, showcasing the latest discoveries and future challenges related to the DepMap. The symposium aimed to facilitate interactions among participants across Europe, encourage idea exchange with leading experts, and present their work and future projects. Importantly, it sparked discussions on future endeavours, such as screening more complex cancer models and accounting for tumour evolution.

Activation priming and cytokine polyfunctionality modulate the enhanced functionality of low-affinity CD19 CAR T cells.

We recently described a low-affinity second-generation CD19 chimeric antigen receptor (CAR) CAT that showed enhanced expansion, cytotoxicity, and antitumor efficacy compared with the high-affinity (FMC63-based) CAR used in tisagenlecleucel, in preclinical models. Furthermore, CAT demonstrated an excellent toxicity profile, enhanced in vivo expansion, and long-term persistence in a phase 1 clinical study. To understand the molecular mechanisms behind these properties of CAT CAR T cells, we performed a systematic in vitro characterization of the transcriptomic (RNA sequencing) and protein (cytometry by time of flight) changes occurring in T cells expressing low-affinity vs high-affinity CD19 CARs following stimulation with CD19-expressing cells. Our results show that CAT CAR T cells exhibit enhanced activation to CD19 stimulation and a distinct transcriptomic and protein profile, with increased activation and cytokine polyfunctionality compared with FMC63 CAR T cells. We demonstrate that the enhanced functionality of low-affinity CAT CAR T cells is a consequence of an antigen-dependent priming induced by residual CD19-expressing B cells present in the manufacture.

High-dimensional functional phenotyping of preclinical human CAR T cells using mass cytometry.

Here, we present a comprehensive protocol for the generation and functional characterization of chimeric antigen receptor (CAR) T cells and their products by mass cytometry in a reproducible and scalable manner. We describe the production of CAR T cells from human peripheral blood mononuclear cells. We then detail a three-step staining protocol with metal-labeled antibodies and the subsequent mass cytometry analysis. This protocol allows simultaneous characterization of CAR T cell intracellular signaling, activation, proliferation, cytokine production, and phenotype in a single assay.

Thiol isomerases negatively regulate the cellular shedding activity of ADAM17.

ADAM17 (where ADAM is 'a disintegrin and metalloproteinase') can rapidly modulate cell-surface signalling events by the proteolytic release of soluble forms of proligands for cellular receptors. Many regulatory pathways affect the ADAM17 sheddase activity, but the mechanisms for the activation are still not clear. We have utilized a cell-based ADAM17 assay to show that thiol isomerases, specifically PDI (protein disulfide isomerase), could be responsible for maintaining ADAM17 in an inactive form. Down-regulation of thiol isomerases, by changes in the redox environment (for instance as elicited by phorbol ester modulation of mitochondrial reactive oxygen species) markedly enhanced ADAM17 activation. On the basis of ELISA binding studies with novel fragment antibodies against ADAM17 we propose that isomerization of the disulfide bonds in ADAM17, and the subsequent conformational changes, form the basis for the modulation of ADAM17 activity. The shuffling of disulfide bond patterns in ADAMs has been suggested by a number of recent adamalysin crystal structures, with distinct disulfide bond patterns altering the relative orientations of the domains. Such a mechanism is rapid and reversible, and the role of thiol isomerases should be investigated further as a potential factor in the redox regulation of ADAM17.

Deciphering Organoids: High-Dimensional Analysis of Biomimetic Cultures.

Organoids are self-organising stem cell-derived ex vivo cultures widely adopted as biomimetic models of healthy and diseased tissues. Traditional low-dimensional experimental methods such as microscopy and bulk molecular analysis have generated remarkable biological insights from organoids. However, as complex heterocellular systems, organoids are especially well-positioned to take advantage of emerging high-dimensional technologies. In particular, single-cell methods offer considerable opportunities to analyse organoids at unprecedented scale and depth, enabling comprehensive characterisation of cellular processes and spatial organisation underpinning organoid heterogeneity. This review evaluates state-of-the-art analytical methods applied to organoids, discusses the latest advances in single-cell technologies, and speculates on the integration of these two rapidly developing fields.

Intra- and intercellular signaling pathways associated with drug-induced cardiac pathophysiology.

Cardiac physiology and homeostasis are maintained by the interaction of multiple cell types, via both intra- and intercellular signaling pathways. Perturbations in these signaling pathways induced by oncology therapies can reduce cardiac function, ultimately leading to heart failure. As cancer survival increases, related cardiovascular complications are becoming increasingly prevalent, thus identifying the perturbations and cell signaling drivers of cardiotoxicity is increasingly important. Here, we discuss the homotypic and heterotypic cellular interactions that form the basis of intra- and intercellular cardiac signaling pathways, and how oncological agents disrupt these pathways, leading to heart failure. We also highlight the emerging systems biology techniques that can be applied, enabling a deeper understanding of the intra- and intercellular signaling pathways across multiple cell types associated with cardiovascular toxicity.

Multiplexed single-cell analysis of organoid signaling networks.

Organoids are biomimetic tissue models comprising multiple cell types and cell states. Post-translational modification (PTM) signaling networks control cellular phenotypes and are frequently dysregulated in diseases such as cancer. Although signaling networks vary across cell types, there are limited techniques to study cell type-specific PTMs in heterocellular organoids. Here, we present a multiplexed mass cytometry (MC) protocol for single-cell analysis of PTM signaling and cell states in organoids and organoids co-cultured with fibroblasts and leukocytes. We describe how thiol-reactive organoid barcoding in situ (TOBis) enables 35-plex and 126-plex single-cell comparison of organoid cultures and provide a cytometry by time of flight (CyTOF) signaling analysis pipeline (CyGNAL) for computing cell type-specific PTM signaling networks. The TOBis MC protocol takes ~3 d from organoid fixation to data acquisition and can generate single-cell data for >40 antibodies from millions of cells across 126 organoid cultures in a single MC run.