Brain Metastasis from EGFR-Mutated Non-Small Cell Lung Cancer: Secretion of IL11 from Astrocytes Up-Regulates PDL1 and Promotes Immune Escape.

Patients who have non-small cell lung cancer (NSCLC) with epidermal growth factor receptor (EGFR) mutations are more prone to brain metastasis (BM) and poor prognosis. Previous studies showed that the tumor microenvironment of BM in these patients is immunosuppressed, as indicated by reduced T-cell abundance and activity, although the mechanism of this immunosuppression requires further study. This study shows that reactive astrocytes play a critical role in promoting the immune escape of BM from EGFR-mutated NSCLC by increasing the apoptosis of CD8+ T lymphocytes. The increased secretion of interleukin 11(IL11) by astrocytes promotes the expression of PDL1 in BM, and this is responsible for the increased apoptosis of T lymphocytes. IL11 functions as a ligand of EGFR, and this binding activates EGFR and downstream signaling to increase the expression of PDL1, culminating in the immune escape of tumor cells. IL11 also promotes immune escape by binding to its intrinsic receptor (IL11Rα/glycoprotein 130 [gp130]). Additional in vivo studies show that the targeted inhibition of gp130 and EGFR suppresses the growth of BM and prolongs the survival time of mice. These results suggest a novel therapeutic strategy for treatment of NSCLC patients with EGFR mutations.

Starfysh integrates spatial transcriptomic and histologic data to reveal heterogeneous tumor-immune hubs.

Spatially resolved gene expression profiling provides insight into tissue organization and cell-cell crosstalk; however, sequencing-based spatial transcriptomics (ST) lacks single-cell resolution. Current ST analysis methods require single-cell RNA sequencing data as a reference for rigorous interpretation of cell states, mostly do not use associated histology images and are not capable of inferring shared neighborhoods across multiple tissues. Here we present Starfysh, a computational toolbox using a deep generative model that incorporates archetypal analysis and any known cell type markers to characterize known or new tissue-specific cell states without a single-cell reference. Starfysh improves the characterization of spatial dynamics in complex tissues using histology images and enables the comparison of niches as spatial hubs across tissues. Integrative analysis of primary estrogen receptor (ER)-positive breast cancer, triple-negative breast cancer (TNBC) and metaplastic breast cancer (MBC) tissues led to the identification of spatial hubs with patient- and disease-specific cell type compositions and revealed metabolic reprogramming shaping immunosuppressive hubs in aggressive MBC.

Coordinated Immune Cell Networks in the Bone Marrow Microenvironment Define the Graft versus Leukemia Response with Adoptive Cellular Therapy.

Understanding how intra-tumoral immune populations coordinate to generate anti-tumor responses following therapy can guide precise treatment prioritization. We performed systematic dissection of an established adoptive cellular therapy, donor lymphocyte infusion (DLI), by analyzing 348,905 single-cell transcriptomes from 74 longitudinal bone-marrow samples of 25 patients with relapsed myeloid leukemia; a subset was evaluated by protein-based spatial analysis. In acute myelogenous leukemia (AML) responders, diverse immune cell types within the bone-marrow microenvironment (BME) were predicted to interact with a clonally expanded population of ZNF683 + GZMB + CD8+ cytotoxic T lymphocytes (CTLs) which demonstrated in vitro specificity for autologous leukemia. This population, originating predominantly from the DLI product, expanded concurrently with NK and B cells. AML nonresponder BME revealed a paucity of crosstalk and elevated TIGIT expression in CD8+ CTLs. Our study highlights recipient BME differences as a key determinant of effective anti-leukemia response and opens new opportunities to modulate cell-based leukemia-directed therapy.

Cellstitch: 3D cellular anisotropic image segmentation via optimal transport.

BACKGROUND: Spatial mapping of transcriptional states provides valuable biological insights into cellular functions and interactions in the context of the tissue. Accurate 3D cell segmentation is a critical step in the analysis of this data towards understanding diseases and normal development in situ. Current approaches designed to automate 3D segmentation include stitching masks along one dimension, training a 3D neural network architecture from scratch, and reconstructing a 3D volume from 2D segmentations on all dimensions. However, the applicability of existing methods is hampered by inaccurate segmentations along the non-stitching dimensions, the lack of high-quality diverse 3D training data, and inhomogeneity of image resolution along orthogonal directions due to acquisition constraints; as a result, they have not been widely used in practice. METHODS: To address these challenges, we formulate the problem of finding cell correspondence across layers with a novel optimal transport (OT) approach. We propose CellStitch, a flexible pipeline that segments cells from 3D images without requiring large amounts of 3D training data. We further extend our method to interpolate internal slices from highly anisotropic cell images to recover isotropic cell morphology. RESULTS: We evaluated the performance of CellStitch through eight 3D plant microscopic datasets with diverse anisotropic levels and cell shapes. CellStitch substantially outperforms the state-of-the art methods on anisotropic images, and achieves comparable segmentation quality against competing methods in isotropic setting. We benchmarked and reported 3D segmentation results of all the methods with instance-level precision, recall and average precision (AP) metrics. CONCLUSIONS: The proposed OT-based 3D segmentation pipeline outperformed the existing state-of-the-art methods on different datasets with nonzero anisotropy, providing high fidelity recovery of 3D cell morphology from microscopic images.

Multimodal single-cell and whole-genome sequencing of small, frozen clinical specimens.

Single-cell genomics enables dissection of tumor heterogeneity and molecular underpinnings of drug response at an unprecedented resolution1-11. However, broad clinical application of these methods remains challenging, due to several practical and preanalytical challenges that are incompatible with typical clinical care workflows, namely the need for relatively large, fresh tissue inputs. In the present study, we show that multimodal, single-nucleus (sn)RNA/T cell receptor (TCR) sequencing, spatial transcriptomics and whole-genome sequencing (WGS) are feasible from small, frozen tissues that approximate routinely collected clinical specimens (for example, core needle biopsies). Compared with data from sample-matched fresh tissue, we find a similar quality in the biological outputs of snRNA/TCR-seq data, while reducing artifactual signals and compositional biases introduced by fresh tissue processing. Profiling sequentially collected melanoma samples from a patient treated in the KEYNOTE-001 trial12, we resolved cellular, genomic, spatial and clonotype dynamics that represent molecular patterns of heterogeneous intralesional evolution during anti-programmed cell death protein 1 therapy. To demonstrate applicability to banked biospecimens of rare diseases13, we generated a single-cell atlas of uveal melanoma liver metastasis with matched WGS data. These results show that single-cell genomics from archival, clinical specimens is feasible and provides a framework for translating these methods more broadly to the clinical arena.

CRISPRmap: Sequencing-free optical pooled screens mapping multi-omic phenotypes in cells and tissue.

Pooled genetic screens are powerful tools to study gene function in a high-throughput manner. Typically, sequencing-based screens require cell lysis, which limits the examination of critical phenotypes such as cell morphology, protein subcellular localization, and cell-cell/tissue interactions. In contrast, emerging optical pooled screening methods enable the investigation of these spatial phenotypes in response to targeted CRISPR perturbations. In this study, we report a multi-omic optical pooled CRISPR screening method, which we have named CRISPRmap. Our method combines a novel in situ CRISPR guide identifying barcode readout approach with concurrent multiplexed immunofluorescence and in situ RNA detection. CRISPRmap barcodes are detected and read out through combinatorial hybridization of DNA oligos, enhancing barcode detection efficiency, while reducing both dependency on third party proprietary sequencing reagents and assay cost. Notably, we conducted a multi-omic base-editing screen in a breast cancer cell line on core DNA damage repair genes involved in the homologous recombination and Fanconi anemia pathways investigating how nucleotide variants in those genes influence DNA damage signaling and cell cycle regulation following treatment with ionizing radiation or DNA damaging agents commonly used for cancer therapy. Approximately a million cells were profiled with our multi-omic approach, providing a comprehensive phenotypic assessment of the functional consequences of the studied variants. CRISPRmap enabled us to pinpoint likely-pathogenic patient-derived mutations that were previously classified as variants of unknown clinical significance. Furthermore, our approach effectively distinguished barcodes of a pooled library in tumor tissue, and we coupled it with cell-type and molecular phenotyping by cyclic immunofluorescence. Multi-omic spatial analysis of how CRISPR-perturbed cells respond to various environmental cues in the tissue context offers the potential to significantly expand our understanding of tissue biology in both health and disease.