Single-Cell Map of Diverse Immune Phenotypes in the Breast Tumor Microenvironment.

Knowledge of immune cell phenotypes in the tumor microenvironment is essential for understanding mechanisms of cancer progression and immunotherapy response. We profiled 45,000 immune cells from eight breast carcinomas, as well as matched normal breast tissue, blood, and lymph nodes, using single-cell RNA-seq. We developed a preprocessing pipeline, SEQC, and a Bayesian clustering and normalization method, Biscuit, to address computational challenges inherent to single-cell data. Despite significant similarity between normal and tumor tissue-resident immune cells, we observed continuous phenotypic expansions specific to the tumor microenvironment. Analysis of paired single-cell RNA and T cell receptor (TCR) sequencing data from 27,000 additional T cells revealed the combinatorial impact of TCR utilization on phenotypic diversity. Our results support a model of continuous activation in T cells and do not comport with the macrophage polarization model in cancer. Our results have important implications for characterizing tumor-infiltrating immune cells.

Full-length haplotype reconstruction to infer the structure of heterogeneous virus populations.

Next-generation sequencing (NGS) technologies enable new insights into the diversity of virus populations within their hosts. Diversity estimation is currently restricted to single-nucleotide variants or to local fragments of no more than a few hundred nucleotides defined by the length of sequence reads. To study complex heterogeneous virus populations comprehensively, novel methods are required that allow for complete reconstruction of the individual viral haplotypes. Here, we show that assembly of whole viral genomes of ∼8600 nucleotides length is feasible from mixtures of heterogeneous HIV-1 strains derived from defined combinations of cloned virus strains and from clinical samples of an HIV-1 superinfected individual. Haplotype reconstruction was achieved using optimized experimental protocols and computational methods for amplification, sequencing and assembly. We comparatively assessed the performance of the three NGS platforms 454 Life Sciences/Roche, Illumina and Pacific Biosciences for this task. Our results prove and delineate the feasibility of NGS-based full-length viral haplotype reconstruction and provide new tools for studying evolution and pathogenesis of viruses.

Spatial analysis of recurrent glioblastoma reveals perivascular niche organization.

Tumor evolution is driven by genetic variation; however, it is the tumor microenvironment (TME) that provides the selective pressure contributing to evolution in cancer. Despite high histopathological heterogeneity within glioblastoma (GBM), the most aggressive brain tumor, the interactions between the genetically distinct GBM cells and the surrounding TME are not fully understood. To address this, we analyzed matched primary and recurrent GBM archival tumor tissues with imaging-based techniques aimed to simultaneously evaluate tumor tissues for the presence of hypoxic, angiogenic, and inflammatory niches, extracellular matrix (ECM) organization, TERT promoter mutational status, and several oncogenic amplifications on the same slide and location. We found that the relationships between genetic and TME diversity are different in primary and matched recurrent tumors. Interestingly, the texture of the ECM, identified by label-free reflectance imaging, was predictive of single-cell genetic traits present in the tissue. Moreover, reflectance of ECM revealed structured organization of the perivascular niche in recurrent GBM, enriched in immunosuppressive macrophages. Single-cell spatial transcriptomics further confirmed the presence of the niche-specific macrophage populations and identified interactions between endothelial cells, perivascular fibroblasts, and immunosuppressive macrophages. Our results underscore the importance of GBM tissue organization in tumor evolution and highlight genetic and spatial dependencies.

Cell identity revealed by precise cell cycle state mapping links data modalities.

Several methods for cell cycle inference from sequencing data exist and are widely adopted. In contrast, methods for classification of cell cycle state from imaging data are scarce. We have for the first time integrated sequencing and imaging derived cell cycle pseudo-times for assigning 449 imaged cells to 693 sequenced cells at an average resolution of 3.4 and 2.4 cells for sequencing and imaging data respectively. Data integration revealed thousands of pathways and organelle features that are correlated with each other, including several previously known interactions and novel associations. The ability to assign the transcriptome state of a profiled cell to its closest living relative, which is still actively growing and expanding opens the door for genotype-phenotype mapping at single cell resolution forward in time.

Peristromal niches protect lung cancers from targeted therapies through a combined effect of multiple molecular mediators.

Targeted therapies directed against oncogenic signaling addictions, such as inhibitors of ALK in ALK+ NSCLC often induce strong and durable clinical responses. However, they are not curative in metastatic cancers, as some tumor cells persist through therapy, eventually developing resistance. Therapy sensitivity can reflect not only cell-intrinsic mechanisms but also inputs from stromal microenvironment. Yet, the contribution of tumor stroma to therapeutic responses in vivo remains poorly defined. To address this gap of knowledge, we assessed the contribution of stroma-mediated resistance to therapeutic responses to the frontline ALK inhibitor alectinib in xenograft models of ALK+ NSCLC. We found that stroma-proximal tumor cells are partially protected against cytostatic effects of alectinib. This effect is observed not only in remission, but also during relapse, indicating the strong contribution of stroma-mediated resistance to both persistence and resistance. This therapy-protective effect of the stromal niche reflects a combined action of multiple mechanisms, including growth factors and extracellular matrix components. Consequently, despite improving alectinib responses, suppression of any individual resistance mechanism was insufficient to fully overcome the protective effect of stroma. Focusing on shared collateral sensitivity of persisters offered a superior therapeutic benefit, especially when using an antibody-drug conjugate with bystander effect to limit therapeutic escape. These findings indicate that stroma-mediated resistance might be the major contributor to both residual and progressing disease and highlight the limitation of focusing on suppressing a single resistance mechanism at a time.

PhysiCOOL: A generalized framework for model Calibration and Optimization Of modeLing projects.

In silico models of biological systems are usually very complex and rely on a large number of parameters describing physical and biological properties that require validation. As such, parameter space exploration is an essential component of computational model development to fully characterize and validate simulation results. Experimental data may also be used to constrain parameter space (or enable model calibration) to enhance the biological relevance of model parameters. One widely used computational platform in the mathematical biology community is PhysiCell, which provides a standardized approach to agent-based models of biological phenomena at different time and spatial scales. Nonetheless, one limitation of PhysiCell is the lack of a generalized approach for parameter space exploration and calibration that can be run without high-performance computing access. Here, we present PhysiCOOL, an open-source Python library tailored to create standardized calibration and optimization routines for PhysiCell models.

Virtual alignment of pathology image series for multi-gigapixel whole slide images.

Interest in spatial omics is on the rise, but generation of highly multiplexed images remains challenging, due to cost, expertise, methodical constraints, and access to technology. An alternative approach is to register collections of whole slide images (WSI), generating spatially aligned datasets. WSI registration is a two-part problem, the first being the alignment itself and the second the application of transformations to huge multi-gigapixel images. To address both challenges, we developed Virtual Alignment of pathoLogy Image Series (VALIS), software which enables generation of highly multiplexed images by aligning any number of brightfield and/or immunofluorescent WSI, the results of which can be saved in the ome.tiff format. Benchmarking using publicly available datasets indicates VALIS provides state-of-the-art accuracy in WSI registration and 3D reconstruction. Leveraging existing open-source software tools, VALIS is written in Python, providing a free, fast, scalable, robust, and easy-to-use pipeline for registering multi-gigapixel WSI, facilitating downstream spatial analyses.

Addressing persistent challenges in digital image analysis of cancerous tissues.

The National Cancer Institute (NCI) supports many research programs and consortia, many of which use imaging as a major modality for characterizing cancerous tissue. A trans-consortia Image Analysis Working Group (IAWG) was established in 2019 with a mission to disseminate imaging-related work and foster collaborations. In 2022, the IAWG held a virtual hackathon focused on addressing challenges of analyzing high dimensional datasets from fixed cancerous tissues. Standard image processing techniques have automated feature extraction, but the next generation of imaging data requires more advanced methods to fully utilize the available information. In this perspective, we discuss current limitations of the automated analysis of multiplexed tissue images, the first steps toward deeper understanding of these limitations, what possible solutions have been developed, any new or refined approaches that were developed during the Image Analysis Hackathon 2022, and where further effort is required. The outstanding problems addressed in the hackathon fell into three main themes: 1) challenges to cell type classification and assessment, 2) translation and visual representation of spatial aspects of high dimensional data, and 3) scaling digital image analyses to large (multi-TB) datasets. We describe the rationale for each specific challenge and the progress made toward addressing it during the hackathon. We also suggest areas that would benefit from more focus and offer insight into broader challenges that the community will need to address as new technologies are developed and integrated into the broad range of image-based modalities and analytical resources already in use within the cancer research community.

Interactions between ploidy and resource availability shape clonal interference at initiation and recurrence of glioblastoma.

Glioblastoma (GBM) is the most aggressive form of primary brain tumor. Complete surgical resection of GBM is almost impossible due to the infiltrative nature of the cancer. While no evidence for recent selection events have been found after diagnosis, the selective forces that govern gliomagenesis are strong, shaping the tumor's cell composition during the initial progression to malignancy with late consequences for invasiveness and therapy response. We present a mathematical model that simulates the growth and invasion of a glioma, given its ploidy level and the nature of its brain tissue micro-environment (TME), and use it to make inferences about GBM initiation and response to standard-of-care treatment. We approximate the spatial distribution of resource access in the TME through integration of in-silico modelling, multi-omics data and image analysis of primary and recurrent GBM. In the pre-malignant setting, our in-silico results suggest that low ploidy cancer cells are more resistant to starvation-induced cell death. In the malignant setting, between first and second surgery, simulated tumors with different ploidy compositions progressed at different rates. Whether higher ploidy predicted fast recurrence, however, depended on the TME. Historical data supports this dependence on TME resources, as shown by a significant correlation between the median glucose uptake rates in human tissues and the median ploidy of cancer types that arise in the respective tissues (Spearman r = -0.70; P = 0.026). Taken together our findings suggest that availability of metabolic substrates in the TME drives different cell fate decisions for cancer cells with different ploidy and shapes GBM disease initiation and relapse characteristics.

Sparcle: assigning transcripts to cells in multiplexed images.

Motivation: Imaging-based spatial transcriptomics has the power to reveal patterns of single-cell gene expression by detecting mRNA transcripts as individually resolved spots in multiplexed images. However, molecular quantification has been severely limited by the computational challenges of segmenting poorly outlined, overlapping cells and of overcoming technical noise; the majority of transcripts are routinely discarded because they fall outside the segmentation boundaries. This lost information leads to less accurate gene count matrices and weakens downstream analyses, such as cell type or gene program identification. Results: Here, we present Sparcle, a probabilistic model that reassigns transcripts to cells based on gene covariation patterns and incorporates spatial features such as distance to nucleus. We demonstrate its utility on both multiplexed error-robust fluorescence in situ hybridization, single-molecule FISH data, probabilistic cell typing in situ sequencing, spatially resolved transcript amplicon readout mapping and MERFISH from Vizgen. Sparcle improves transcript assignment, providing more realistic per-cell quantification of each gene, better delineation of cell boundaries and improved cluster assignments. Critically, our approach does not require an accurate segmentation and is agnostic to technological platform. Availability and implementation: The code is available at: https://github.com/sandhya212/Sparcle_for_spot_reassignments. Contact: sandhya.prabhakaran@moffitt.org. Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Developing tools for analyzing and viewing multiplexed images.

Dr. Prabhakaran and Dr Gatenbee are research scientists in Anderson's lab and have developed Mistic, a publicly available tool that simultaneously views multiplexed images and assists in gaining biological and clinical insights into patients' data. They discuss the role of mathematical modeling in translational cancer research and clinical decision making and describe how mathematical modeling fits into the data science definition.

Mistic: An open-source multiplexed image t-SNE viewer.

Understanding the complex ecology of a tumor tissue and the spatiotemporal relationships between its cellular and microenvironment components is becoming a key component of translational research, especially in immuno-oncology. The generation and analysis of multiplexed images from patient samples is of paramount importance to facilitate this understanding. Here, we present Mistic, an open-source multiplexed image t-SNE viewer that enables the simultaneous viewing of multiple 2D images rendered using multiple layout options to provide an overall visual preview of the entire dataset. In particular, the positions of the images can be t-SNE or UMAP coordinates. This grouped view of all images allows an exploratory understanding of the specific expression pattern of a given biomarker or collection of biomarkers across all images, helps to identify images expressing a particular phenotype, and can help select images for subsequent downstream analysis. Currently, there is no freely available tool to generate such image t-SNEs.

IgA-Dominated Humoral Immune Responses Govern Patients' Outcome in Endometrial Cancer.

Recent studies suggest that B cells could play an important role in the tumor microenvironment. However, the role of humoral responses in endometrial cancer remains insufficiently investigated. Using a cohort of 107 patients with different histological subtypes of endometrial carcinoma, we evaluated the role of coordinated humoral and cellular adaptive immune responses in endometrial cancer. Concomitant accumulation of T, B, and plasma cells at tumor beds predicted better survival. However, only B-cell markers corresponded with prolonged survival specifically in high-grade endometrioid type and serous tumors. Immune protection was associated with class-switched IgA and, to a lesser extent, IgG. Expressions of polymeric immunoglobulin receptor (pIgR) by tumor cells and its occupancy by IgA were superior predictors of outcome and correlated with defects in methyl-directed DNA mismatch repair. Mechanistically, pIgR-dependent, antigen-independent IgA occupancy drove activation of inflammatory pathways associated with IFN and TNF signaling in tumor cells, along with apoptotic and endoplasmic reticulum stress pathways, while thwarting DNA repair mechanisms. Together, these findings suggest that coordinated humoral and cellular immune responses, characterized by IgA:pIgR interactions in tumor cells, determine the progression of human endometrial cancer as well as the potential for effective immunotherapies. SIGNIFICANCE: This study provides new insights into the crucial role of humoral immunity in human endometrial cancer, providing a rationale for designing novel immunotherapies against this prevalent malignancy. See related commentary by Osorio and Zamarin, p. 766.

Ovarian cancer immunogenicity is governed by a narrow subset of progenitor tissue-resident memory T cells.

Despite repeated associations between T cell infiltration and outcome, human ovarian cancer remains poorly responsive to immunotherapy. We report that the hallmarks of tumor recognition in ovarian cancer-infiltrating T cells are primarily restricted to tissue-resident memory (TRM) cells. Single-cell RNA/TCR/ATAC sequencing of 83,454 CD3+CD8+CD103+CD69+ TRM cells and immunohistochemistry of 122 high-grade serous ovarian cancers shows that only progenitor (TCF1low) tissue-resident T cells (TRMstem cells), but not recirculating TCF1+ T cells, predict ovarian cancer outcome. TRMstem cells arise from transitional recirculating T cells, which depends on antigen affinity/persistence, resulting in oligoclonal, trogocytic, effector lymphocytes that eventually become exhausted. Therefore, ovarian cancer is indeed an immunogenic disease, but that depends on ∼13% of CD8+ tumor-infiltrating T cells (∼3% of CD8+ clonotypes), which are primed against high-affinity antigens and maintain waves of effector TRM-like cells. Our results define the signature of relevant tumor-reactive T cells in human ovarian cancer, which could be applicable to other tumors with unideal mutational burden.

Dirichlet Process Mixture Model for Correcting Technical Variation in Single-Cell Gene Expression Data.

We introduce an iterative normalization and clustering method for single-cell gene expression data. The emerging technology of single-cell RNA-seq gives access to gene expression measurements for thousands of cells, allowing discovery and characterization of cell types. However, the data is confounded by technical variation emanating from experimental errors and cell type-specific biases. Current approaches perform a global normalization prior to analyzing biological signals, which does not resolve missing data or variation dependent on latent cell types. Our model is formulated as a hierarchical Bayesian mixture model with cell-specific scalings that aid the iterative normalization and clustering of cells, teasing apart technical variation from biological signals. We demonstrate that this approach is superior to global normalization followed by clustering. We show identifiability and weak convergence guarantees of our method and present a scalable Gibbs inference algorithm. This method improves cluster inference in both synthetic and real single-cell data compared with previous methods, and allows easy interpretation and recovery of the underlying structure and cell types.

HIV Haplotype Inference Using a Propagating Dirichlet Process Mixture Model.

This paper presents a new computational technique for the identification of HIV haplotypes. HIV tends to generate many potentially drug-resistant mutants within the HIV-infected patient and being able to identify these different mutants is important for efficient drug administration. With the view of identifying the mutants, we aim at analyzing short deep sequencing data called reads. From a statistical perspective, the analysis of such data can be regarded as a nonstandard clustering problem due to missing pairwise similarity measures between non-overlapping reads. To overcome this problem we propagate a Dirichlet Process Mixture Model by sequentially updating the prior information from successive local analyses. The model is verified using both simulated and real sequencing data.

Probabilistic inference of viral quasispecies subject to recombination.

RNA viruses exist in their hosts as populations of different but related strains. The virus population, often called quasispecies, is shaped by a combination of genetic change and natural selection. Genetic change is due to both point mutations and recombination events. We present a jumping hidden Markov model that describes the generation of viral quasispecies and a method to infer its parameters from next-generation sequencing data. The model introduces position-specific probability tables over the sequence alphabet to explain the diversity that can be found in the population at each site. Recombination events are indicated by a change of state, allowing a single observed read to originate from multiple sequences. We present a specific implementation of the expectation maximization (EM) algorithm to find maximum a posteriori estimates of the model parameters and a method to estimate the distribution of viral strains in the quasispecies. The model is validated on simulated data, showing the advantage of explicitly taking the recombination process into account, and applied to reads obtained from a clinical HIV sample.