Systems biology approaches to measure and model phenotypic heterogeneity in cancer

Sequential Chromogenic IHC: Spatial Analysis of Lymph Nodes Identifies Contact Interactions between Plasmacytoid Dendritic Cells and Plasmablasts

Recent clinical observations have emphasized the critical role that the spatial organization of immune cells in lymphoid structures plays in the success of cancer immunotherapy and patient survival. However, implementing sequential chromogenic IHC (scIHC) to analyze multiple biomarkers on a single tissue section has been limited because of a lack of a standardized, rigorous guide to the development of customized biomarker panels and a need for user-friendly analysis pipelines that can extract meaningful data. In this context, we provide a comprehensive guide for the development of novel biomarker panels for scIHC, using practical examples and illustrations to highlight the most common complications that can arise during the setup of a new biomarker panel, and provide detailed instructions on how to prevent and detect cross-reactivity between secondary reagents and carryover between detection antibodies. We also developed a novel analysis pipeline based on non-rigid tissue deformation correction, Cellpose-inspired automated cell segmentation, and computational network masking of low-quality data. We applied this biomarker panel and pipeline to study regional lymph nodes from patients with head and neck cancer, identifying novel contact interactions between plasmablasts and plasmacytoid dendritic cells in vivo. Given that Toll-like receptors, which are highly expressed in plasmacytoid dendritic cells, play a key role in vaccine efficacy, the significance of this cell-cell interaction decisively warrants further studies. In summary, this work provides a streamlined approach to the development of customized biomarker panels for scIHC that will ultimately improve our understanding of immune responses in cancer.

Significance

We present a comprehensive guide for developing customized biomarker panels to investigate cell-cell interactions in the context of immune responses in cancer. This approach revealed novel contact interactions between plasmablasts and plasmacytoid dendritic cells in lymph nodes from patients with head and neck cancer.

Analysis and modeling of cancer drug responses using cell cycle phase-specific rate effects

Identifying effective therapeutic treatment strategies is a major challenge to improving outcomes for patients with breast cancer. To gain a comprehensive understanding of how clinically relevant anti-cancer agents modulate cell cycle progression, here we use genetically engineered breast cancer cell lines to track drug-induced changes in cell number and cell cycle phase to reveal drug-specific cell cycle effects that vary across time. We use a linear chain trick (LCT) computational model, which faithfully captures drug-induced dynamic responses, correctly infers drug effects, and reproduces influences on specific cell cycle phases. We use the LCT model to predict the effects of unseen drug combinations and confirm these in independent validation experiments. Our integrated experimental and modeling approach opens avenues to assess drug responses, predict effective drug combinations, and identify optimal drug sequencing strategies.

A single-cell analysis of breast cancer cell lines to study tumour heterogeneity and drug response

Cancer cells within a tumour have heterogeneous phenotypes and exhibit dynamic plasticity. How to evaluate such heterogeneity and its impact on outcome and drug response is still unclear. Here, we transcriptionally profile 35,276 individual cells from 32 breast cancer cell lines to yield a single cell atlas. We find high degree of heterogeneity in the expression of biomarkers. We then train a deconvolution algorithm on the atlas to determine cell line composition from bulk gene expression profiles of tumour biopsies, thus enabling cell line-based patient stratification. Finally, we link results from large-scale in vitro drug screening in cell lines to the single cell data to computationally predict drug responses starting from single-cell profiles. We find that transcriptional heterogeneity enables cells with differential drug sensitivity to co-exist in the same population. Our work provides a framework to determine tumour heterogeneity in terms of cell line composition and drug response.

The Yin-Yang Dynamics in Cancer Pharmacogenomics and Personalized Medicine

The enormous heterogeneity of cancer systems has made it very challenging to overcome drug resistance and adverse reactions to achieve personalized therapies. Recent developments in systems biology, especially the perception of cancer as the complex adaptive system (CAS), may help meet the challenges by deciphering the interactions at various levels from the molecular, cellular, tissue-organ, to the whole organism. The ubiquitous Yin-Yang interactions among the coevolving components, including the genes and proteins, decide their spatiotemporal features at various stages from cancer initiation to metastasis. The Yin-Yang imbalances across different systems levels, from genetic mutations to tumor cells adaptation, have been related to the intra- and inter-tumoral heterogeneity in the micro- and macro-environments. At the molecular and cellular levels, dysfunctional Yin-Yang dynamics in the cytokine networks, mitochondrial activities, redox systems, apoptosis, and metabolism can contribute to tumor cell growth and escape of immune surveillance. Up to the organism and system levels, the Yin-Yang imbalances in the cancer microenvironments can lead to different phenotypes from breast cancer to leukemia. These factors may be considered the systems-based biomarkers and treatment targets. The features of adaptation and nonlinearity in Yin-Yang dynamical interactions should be addressed by individualized drug combinations, dosages, intensities, timing, and frequencies at different cancer stages. The comprehensive "Yin-Yang dynamics" framework would enable powerful approaches for personalized and systems medicine strategies.

Traditional chinese medicine syndromes classification associates with tumor cell and microenvironment heterogeneity in colorectal cancer: a single cell RNA sequencing analysis

BACKGROUND

Colorectal cancer (CRC) is one of the common gastrointestinal malignancies, tumor heterogeneity is the main cause of refractory CRC. Syndrome differentiation is the premise of individualized treatment of traditional Chinese medicine (TCM), but TCM syndrome lacks objective identification in CRC. This study is to investigate the correlation and significance of tumor heterogeneity and TCM syndromes classification in CRC.

METHODS

In this study, we using scRNA-seq technology, investigate the significance of tumor heterogeneity in TCM syndromes classification on CRC.

RESULTS

The results showed that 662 cells isolated from 11 primary CRC tumors are divided into 14 different cell clusters, and each cell subtype and its genes have different functions and signal transduction pathways, indicating significant heterogeneity. CRC tumor cell clusters have different proportions in Excess, Deficiency and Deficiency-Excess syndromes, and have their own characteristic genes, gene co-expression networks, gene functional interpretations as well as monocle functional evolution. Moreover, there were significant differences between the high expressions of MUC2, REG4, COL1A2, POSTN, SDPR, GPX1, ELF3, KRT8, KRT18, KRT19, FN1, SERPINE1, TCF4 and ZEB1 genes in Excess and Deficiency syndrome classification in CRC (P < 0.01).

CONCLUSIONS

The Excess and Deficiency syndromes classification may be related to tumor heterogeneity and its microenvironment in CRC.

Biophysical Informatics Approach For Quantifying Phenotypic Heterogeneity In Cancer Cell Migration In Confined Microenvironments

MOTIVATION

Cancer cell heterogeneity can manifest genetically and phenotypically. Bioinformatics methods have been used to analyze complex genomics and transcriptomics data, but have not been well-established for analyzing biophysical data of phenotypically heterogeneous tumor cells. Here, we take an informatics approach to analyze the biophysical data of MDA-MB-231 cells, a widely used breast cancer cell line, during their spontaneous migration through confined environments. Experimentally, we vary the constriction microchannel geometries (wide channel, short constriction, and long constriction) and apply drug treatments. We find that cells in the short constriction are similar in morphology to the cells in the wide channel. However, their fluorescence profiles are comparable to those in the long constriction. We demonstrate that the cell migratory phenotype is correlated more to mitochondria in a non-confined environment and more to actin in a confined environment. We demonstrate that the cells' migratory phenotypes are altered by ciliobrevin D, a dynein inhibitor, in both confined and non-confined environments. Overall, our approach elucidates phenotypic heterogeneity in cancer cells under confined microenvironments at single-cell resolution.

RESULTS

Here, we apply a bioinformatics approach to a single cell invasion assay. We demonstrate that this method can determine distinctions in morphology, cytoskeletal activities, and mitochondrial activities under various geometric constraints and for cells of different speeds. Our approach can be readily adapted to various heterogeneity studies for different types of input biophysical data. In addition, this approach can be applied to studies related to biophysical changes due to differences in external stimuli, such as treatment effects on cellular and subcellular activities, at single-cell resolution. Finally, as similar bioinformatics methods have been widely applied in studies of genetic heterogeneity, biophysical information extracted using this approach can be analyzed together with the genetic data to relate genetic and phenotypic heterogeneity.

AVAILABILITY

The data that support the findings of this study are available from the corresponding author upon reasonable request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Understanding breast cancer heterogeneity through non-genetic heterogeneity

Intricacy in treatment and diagnosis of breast cancer has been an obstacle due to genotype and phenotype heterogeneity. Understanding of non-genetic heterogeneity mechanisms along with considering role of genetic heterogeneity may fill the gaps in landscape painting of heterogeneity. The main factors contribute to non-genetic heterogeneity including: transcriptional pulsing/bursting or discontinuous transcriptions, stochastic partitioning of components at cell division and various signal transduction from tumor ecosystem. Throughout this review, we desired to provide a conceptual framework focused on non-genetic heterogeneity, which has been intended to offer insight into prediction, diagnosis and treatment of breast cancer.

Coupled Feedback Loops Involving PAGE4, EMT and Notch Signaling Can Give Rise to Non-genetic Heterogeneity in Prostate Cancer Cells

Non-genetic heterogeneity is emerging as a crucial factor underlying therapy resistance in multiple cancers. However, the design principles of regulatory networks underlying non-genetic heterogeneity in cancer remain poorly understood. Here, we investigate the coupled dynamics of feedback loops involving (a) oscillations in androgen receptor (AR) signaling mediated through an intrinsically disordered protein PAGE4, (b) multistability in epithelial-mesenchymal transition (EMT), and c) Notch-Delta-Jagged signaling mediated cell-cell communication, each of which can generate non-genetic heterogeneity through multistability and/or oscillations. Our results show how different coupling strengths between AR and EMT signaling can lead to monostability, bistability, or oscillations in the levels of AR, as well as propagation of oscillations to EMT dynamics. These results reveal the emergent dynamics of coupled oscillatory and multi-stable systems and unravel mechanisms by which non-genetic heterogeneity in AR levels can be generated, which can act as a barrier to most existing therapies for prostate cancer patients.

Gene expression profiles of inflammatory breast cancer reveal high heterogeneity across the epithelial-hybrid-mesenchymal spectrum

•

No unique genome signature or molecular therapy exists for inflammatory breast cancer (IBC), a highly aggressive breast cancer with a 5-year survival rate of less than 30%.

•

We show that various gene lists proposed as molecular footprints of IBC have no overlap and thus very limited predictive accuracy in identifying IBC samples.

•

We observed that single-sample gene set enrichment analysis (ssGSEA) of IBC samples along the epithelial-hybrid-mesenchymal spectrum can help IBC identification.

•

IBC samples robustly displayed a higher coefficient of variation in terms of EMT scores, as compared to non-IBC samples.

•

Higher heterogeneity along the epithelial-hybrid-mesenchymal spectrum can be regarded to be a hallmark of IBC and a possibly useful biomarker.

Inflammatory breast cancer (IBC) is a highly aggressive breast cancer that metastasizes largely via tumor emboli, and has a 5-year survival rate of less than 30%. No unique genomic signature has yet been identified for IBC nor has any specific molecular therapeutic been developed to manage the disease. Thus, identifying gene expression signatures specific to IBC remains crucial. Here, we compare various gene lists that have been proposed as molecular footprints of IBC using different clinical samples as training and validation sets and using independent training algorithms, and determine their accuracy in identifying IBC samples in three independent datasets. We show that these gene lists have little to no mutual overlap, and have limited predictive accuracy in identifying IBC samples. Despite this inconsistency, single-sample gene set enrichment analysis (ssGSEA) of IBC samples correlate with their position on the epithelial-hybrid-mesenchymal spectrum. This positioning, together with ssGSEA scores, improves the accuracy of IBC identification across the three independent datasets. Finally, we observed that IBC samples robustly displayed a higher coefficient of variation in terms of EMT scores, as compared to non-IBC samples. Pending verification that this patient-to-patient variability extends to intratumor heterogeneity within a single patient, these results suggest that higher heterogeneity along the epithelial-hybrid-mesenchymal spectrum can be regarded to be a hallmark of IBC and a possibly useful biomarker.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

Measurement and models accounting for cell death capture hidden variation in compound response

Cancer cell sensitivity or resistance is almost universally quantified through a direct or surrogate measure of cell number. However, compound responses can occur through many distinct phenotypic outcomes, including changes in cell growth, apoptosis, and non-apoptotic cell death. These outcomes have divergent effects on the tumor microenvironment, immune response, and resistance mechanisms. Here, we show that quantifying cell viability alone is insufficient to distinguish between these compound responses. Using an alternative assay and drug-response analysis amenable to high-throughput measurement, we find that compounds with identical viability outcomes can have very different effects on cell growth and death. Moreover, additive compound pairs with distinct growth/death effects can appear synergistic when only assessed by viability. Overall, these results demonstrate an approach to incorporating measurements of cell death when characterizing a pharmacologic response.