FunSpace: A functional and spatial analytic approach to cell imaging data using entropy measures.

Spatial heterogeneity in the tumor microenvironment (TME) plays a critical role in gaining insights into tumor development and progression. Conventional metrics typically capture the spatial differential between TME cellular patterns by either exploring the cell distributions in a pairwise fashion or aggregating the heterogeneity across multiple cell distributions without considering the spatial contribution. As such, none of the existing approaches has fully accounted for the simultaneous heterogeneity caused by both cellular diversity and spatial configurations of multiple cell categories. In this article, we propose an approach to leverage spatial entropy measures at multiple distance ranges to account for the spatial heterogeneity across different cellular organizations. Functional principal component analysis (FPCA) is applied to estimate FPC scores which are then served as predictors in a Cox regression model to investigate the impact of spatial heterogeneity in the TME on survival outcome, potentially adjusting for other confounders. Using a non-small cell lung cancer dataset (n = 153) as a case study, we found that the spatial heterogeneity in the TME cellular composition of CD14+ cells, CD19+ B cells, CD4+ and CD8+ T cells, and CK+ tumor cells, had a significant non-zero effect on the overall survival (p = 0.027). Furthermore, using a publicly available multiplexed ion beam imaging (MIBI) triple-negative breast cancer dataset (n = 33), our proposed method identified a significant impact of cellular interactions between tumor and immune cells on the overall survival (p = 0.046). In simulation studies under different spatial configurations, the proposed method demonstrated a high predictive power by accounting for both clinical effect and the impact of spatial heterogeneity.

A platform-independent framework for phenotyping of multiplex tissue imaging data.

Multiplex imaging is a powerful tool to analyze the structural and functional states of cells in their morphological and pathological contexts. However, hypothesis testing with multiplex imaging data is a challenging task due to the extent and complexity of the information obtained. Various computational pipelines have been developed and validated to extract knowledge from specific imaging platforms. A common problem with customized pipelines is their reduced applicability across different imaging platforms: Every multiplex imaging technique exhibits platform-specific characteristics in terms of signal-to-noise ratio and acquisition artifacts that need to be accounted for to yield reliable and reproducible results. We propose a pixel classifier-based image preprocessing step that aims to minimize platform-dependency for all multiplex image analysis pipelines. Signal detection and noise reduction as well as artifact removal can be posed as a pixel classification problem in which all pixels in multiplex images can be assigned to two general classes of either I) signal of interest or II) artifacts and noise. The resulting feature representation maps contain pixel-scale representations of the input data, but exhibit significantly increased signal-to-noise ratios with normalized pixel values as output data. We demonstrate the validity of our proposed image preprocessing approach by comparing the results of two well-accepted and widely-used image analysis pipelines.

Dynamical analysis of the fission yeast cell cycle via Markov chain.

The cell cycle is a complex network involved in the regulation of cell growth and proliferation. Intrinsic molecular noise in gene expression in the cell cycle network can generate fluctuations in protein concentration. How the cell cycle network maintains its robust transitions between cell cycle phases in the presence of these fluctuations remains unclear. To understand the complex and robust behavior of the cell cycle system in the presence of intrinsic noise, we developed a Markov model for the fission yeast cell cycle system. We quantified the effect of noise on gene and protein activity and on the probability of transition between different phases of the cell cycle. Our analysis shows how network perturbations decide the fate of the cell. Our model predicts that the cell cycle pathway (subsequent transitions from [Formula: see text]) is the most robust and probable pathway among all possible trajectories in the cell cycle network. We performed a sensitivity analysis to find correlations between protein interaction weights and transition probabilities between cell cycle phases. The sensitivity analysis predicts how network perturbations affect the transition probability between different cell cycle phases and, consequently, affect different cell fates, thus, forming testable in vitro/in vivo hypotheses. Our simulation results agree with published experimental findings and reveal how noise in the cell cycle regulatory network can affect cell cycle progression.

A stochastic model of size control in the budding yeast cell cycle.

BACKGROUND: Cell size is a key characteristic that significantly affects many aspects of cellular physiology. There are specific control mechanisms during cell cycle that maintain the cell size within a range from generation to generation. Such control mechanisms introduce substantial variabilities to important properties of the cell cycle such as growth and division. To quantitatively study the effect of such variability in progression through cell cycle, detailed stochastic models are required. RESULTS: In this paper, a new hybrid stochastic model is proposed to study the effect of molecular noise and size control mechanism on the variabilities in cell cycle of the budding yeast Saccharomyces cerevisiae. The proposed model provides an accurate, yet computationally efficient approach for simulation of an intricate system by integrating the deterministic and stochastic simulation schemes. The developed hybrid stochastic model can successfully capture several key features of the cell cycle observed in experimental data. In particular, the proposed model: 1) confirms that the majority of noise in size control stems from low copy numbers of transcripts in the G1 phase, 2) identifies the size and time regulation modules in the size control mechanism, and 3) conforms with phenotypes of early G1 mutants in exquisite detail. CONCLUSIONS: Hybrid stochastic modeling approach can be used to provide quantitative descriptions for stochastic properties of the cell cycle within a computationally efficient framework.

The spatial structure of the tumor immune microenvironment can explain and predict patient response in high-grade serous carcinoma.

Ovarian cancer is the deadliest gynecological malignancy, and therapeutic options and mortality rates over the last three decades have largely not changed. Recent studies indicate that the composition of the tumor immune microenvironment (TIME) influences patient outcomes. To improve spatial understanding of the TIME, we performed multiplexed ion beam imaging on 83 human high-grade serous carcinoma tumor samples, identifying about 160,000 cells across 23 cell types. For 77 of these samples meeting inclusion criteria, we generated composition features based on cell type proportions, spatial features based on the distances between cell types, and spatial network features representing cell interactions and cell clustering patterns, which we linked to traditional clinical and immunohistochemical variables and patient overall survival (OS) and progression-free survival (PFS) outcomes. Among these features, we found several significant univariate correlations, including B-cell contact with M1 macrophages (OS hazard ratio HR=0.696, p=0.011, PFS HR=0.734, p=0.039). We then used high-dimensional random forest models to evaluate out-of-sample predictive performance for OS and PFS outcomes and to derive relative feature importance scores for each feature. The top model for predicting low or high PFS used TIME composition and spatial features and achieved an average AUC (area under the receiver-operating characteristic curve) score of 0.71. The results demonstrate the importance of spatial structure in understanding how the TIME contributes to treatment outcomes. Furthermore, the present study provides a generalizable roadmap for spatial analyses of the TIME in ovarian cancer research.

The spatial structure of the tumor immune microenvironment can explain and predict patient response in high-grade serous carcinoma.

Despite ovarian cancer being the deadliest gynecological malignancy, there has been little change to therapeutic options and mortality rates over the last three decades. Recent studies indicate that the composition of the tumor immune microenvironment (TIME) influences patient outcomes but are limited by a lack of spatial understanding. We performed multiplexed ion beam imaging (MIBI) on 83 human high-grade serous carcinoma tumors - one of the largest protein-based, spatially-intact, single-cell resolution tumor datasets assembled - and used statistical and machine learning approaches to connect features of the TIME spatial organization to patient outcomes. Along with traditional clinical/immunohistochemical attributes and indicators of TIME composition, we found that several features of TIME spatial organization had significant univariate correlations and/or high relative importance in high-dimensional predictive models. The top performing predictive model for patient progression-free survival (PFS) used a combination of TIME composition and spatial features. Results demonstrate the importance of spatial structure in understanding how the TIME contributes to treatment outcomes. Furthermore, the present study provides a generalizable roadmap for spatial analyses of the TIME in ovarian cancer research.

A hybrid stochastic model of the budding yeast cell cycle.

The growth and division of eukaryotic cells are regulated by complex, multi-scale networks. In this process, the mechanism of controlling cell-cycle progression has to be robust against inherent noise in the system. In this paper, a hybrid stochastic model is developed to study the effects of noise on the control mechanism of the budding yeast cell cycle. The modeling approach leverages, in a single multi-scale model, the advantages of two regimes: (1) the computational efficiency of a deterministic approach, and (2) the accuracy of stochastic simulations. Our results show that this hybrid stochastic model achieves high computational efficiency while generating simulation results that match very well with published experimental measurements.