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

Mechanical deformation and death of circulating tumor cells in the bloodstream

The circulation of tumor cells through the bloodstream is a significant step in tumor metastasis. To better understand the metastatic process, circulating tumor cell (CTC) survival in the circulation must be explored. While immune interactions with CTCs in recent decades have been examined, research has yet to sufficiently explain some CTC behaviors in blood flow. Studies related to CTC mechanical responses in the bloodstream have recently been conducted to further study conditions under which CTCs might die. While experimental methods can assess the mechanical properties and death of CTCs, increasingly sophisticated computational models are being built to simulate the blood flow and CTC mechanical deformation under fluid shear stresses (FSS) in the bloodstream.Several factors contribute to the mechanical deformation and death of CTCs as they circulate. While FSS can damage CTC structure, diverse interactions between CTCs and blood components may either promote or hinder the next metastatic step-extravasation at a remote site. Overall understanding of how these factors influence the deformation and death of CTCs could serve as a basis for future experiments and simulations, enabling researchers to predict CTC death more accurately. Ultimately, these efforts can lead to improved metastasis-specific therapeutics and diagnostics specific in the future.

Insights into the mechanobiology of cancer metastasis via microfluidic technologies

During cancer metastasis, cancer cells will encounter various microenvironments with diverse physical characteristics. Changes in these physical characteristics such as tension, stiffness, viscosity, compression, and fluid shear can generate biomechanical cues that affect cancer cells, dynamically influencing numerous pathophysiological mechanisms. For example, a dense extracellular matrix drives cancer cells to reorganize their cytoskeleton structures, facilitating confined migration, while this dense and restricted space also acts as a physical barrier that potentially results in nuclear rupture. Identifying these pathophysiological processes and understanding their underlying mechanobiological mechanisms can aid in the development of more effective therapeutics targeted to cancer metastasis. In this review, we outline the advances of engineering microfluidic devices in vitro and their role in replicating tumor microenvironment to mimic in vivo settings. We highlight the potential cellular mechanisms that mediate their ability to adapt to different microenvironments. Meanwhile, we also discuss some important mechanical cues that still remain challenging to replicate in current microfluidic devices in future direction. While much remains to be explored about cancer mechanobiology, we believe the developments of microfluidic devices will reveal how these physical cues impact the behaviors of cancer cells. It will be crucial in the understanding of cancer metastasis, and potentially contributing to better drug development and cancer therapy.

Epithelial and mesenchymal phenotypes determine the dynamics of circulating breast tumor cells in microfluidic capillaries under chemotherapy-induced stress

Circulating tumor cells (CTCs) with different epithelial and mesenchymal phenotypes play distinct roles in the metastatic cascade. However, the influence of their phenotypic traits and chemotherapy on their transit and retention within capillaries remains unclear. To explore this, we developed a microfluidic device comprising 216 microchannels of different widths from 5 to 16 μm to mimic capillaries. This platform allowed us to study the behaviors of human breast cancer epithelial MCF-7 and mesenchymal MDA-MB-231 cells through microchannels under chemotherapy-induced stress. Our results revealed that when the cell diameter to microchannel width ratio exceeded 1.2, MCF-7 cells exhibited higher transit percentages than MDA-MB-231 cells under a flow rate of 0.13 mm/s. Tamoxifen (250 nM) reduced the transit percentage of MCF-7 cells, whereas 100 nM paclitaxel decreased transit percentages for both cell types. These differential responses were partially due to altered cell stiffness following drug treatments. When cells were entrapped at microchannel entrances, tamoxifen, paclitaxel, and high-flow stress (0.5 mm/s) induced a reduction in mitochondrial membrane potential (MMP) in MCF-7 cells. Tamoxifen treatment also elevated reactive oxygen species (ROS) levels in MCF-7 cells. Conversely, MMP and ROS levels in entrapped MDA-MB-231 cells remained unaffected. Consequently, the viability and proliferation of entrapped MCF-7 cells declined under these chemical and physical stress conditions. Our findings emphasize that phenotypically distinct CTCs may undergo selective filtration and exhibit varied responses to chemotherapy in capillaries, thereby impacting cancer metastasis outcomes. This highlights the importance of considering both cell phenotype and drug response to improve treatment strategies.

Microfluidic single-cell migration chip reveals insights into the impact of extracellular matrices on cell movement

Cell migration is a complex process that plays a crucial role in normal physiology and pathologies such as cancer, autoimmune diseases, and mental disorders. Conventional cell migration assays face limitations in tracking a large number of individual migrating cells. To address this challenge, we have developed a high-throughput microfluidic cell migration chip, which seamlessly integrates robotic liquid handling and computer vision to swiftly monitor the movement of 3200 individual cells, providing unparalleled single-cell resolution for discerning distinct behaviors of the fast-moving cell population. This study focuses on the ECM's role in regulating cellular migration, utilizing this cutting-edge microfluidic technology to investigate the impact of ten different ECMs on triple-negative breast cancer cell lines. We found that collagen IV, collagen III, and collagen I coatings were the top enhancers of cell movement. Combining these ECMs increased cell motility, but the effect was sub-additive. Furthermore, we examined 87 compounds and found that while some compounds inhibited migration on all substrates, significantly distinct effects on differently coated substrates were observed, underscoring the importance of considering ECM coating. We also utilized cells expressing a fluorescent actin reporter and observed distinct actin structures in ECM-interacting cells. ScRNA-Seq analysis revealed that ECM coatings induced EMT and enhanced cell migration. Finally, we identified genes that were particularly up-regulated by collagen IV and the selective inhibitors successfully blocked cell migration on collagen IV. Overall, the study provides insights into the impact of various ECMs on cell migration and dynamics of cell movement with implications for developing therapeutic strategies to combat diseases related to cell motility.

Compression drives diverse transcriptomic and phenotypic adaptations in melanoma

Physical forces are prominent during tumor progression. However, it is still unclear how they impact and drive the diverse phenotypes found in cancer. Here, we apply an integrative approach to investigate the impact of compression on melanoma cells. We apply bioinformatics to screen for the most significant compression-induced transcriptomic changes and investigate phenotypic responses. We show that compression-induced transcriptomic changes are associated with both improvement and worsening of patient prognoses. Phenotypically, volumetric compression inhibits cell proliferation and cell migration. It also induces organelle stress and intracellular oxidative stress and increases pigmentation in malignant melanoma cells and normal human melanocytes. Finally, cells that have undergone compression become more resistant to cisplatin treatment. Our findings indicate that volumetric compression is a double-edged sword for melanoma progression and drives tumor evolution.

Biophysical informatics reveals distinctive phenotypic signatures and functional diversity of single-cell lineages

MOTIVATION

In this work, we present an analytical method for quantifying both single-cell morphologies and cell network topologies of tumor cell populations and use it to predict 3D cell behavior.

RESULTS

We utilized a supervised deep learning approach to perform instance segmentation on label-free live cell images across a wide range of cell densities. We measured cell shape properties and characterized network topologies for 136 single-cell clones derived from the YUMM1.7 and YUMMER1.7 mouse melanoma cell lines. Using an unsupervised clustering algorithm, we identified six distinct morphological subclasses. We further observed differences in tumor growth and invasion dynamics across subclasses in an in vitro 3D spheroid model. Compared to existing methods for quantifying 2D or 3D phenotype, our analytical method requires less time, needs no specialized equipment and is capable of much higher throughput, making it ideal for applications such as high-throughput drug screening and clinical diagnosis.

AVAILABILITY AND IMPLEMENTATION

https://github.com/trevor-chan/Melanoma_NetworkMorphology.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Extracellular matrix mechanobiology in cancer cell migration

The extracellular matrix (ECM) is pivotal in modulating tumor progression. Besides chemically stimulating tumor cells, it also offers physical support that orchestrates the sequence of events in the metastatic cascade upon dynamically modulating cell mechanosensation. Understanding this translation between matrix biophysical cues and intracellular signaling has led to rapid growth in the interdisciplinary field of cancer mechanobiology in the last decade. Substantial efforts have been made to develop novel in vitro tumor mimicking platforms to visualize and quantify the mechanical forces within the tissue that dictate tumor cell invasion and metastatic growth. This review highlights recent findings on tumor matrix biophysical cues such as fibrillar arrangement, crosslinking density, confinement, rigidity, topography, and non-linear mechanics and their implications on tumor cell behavior. We also emphasize how perturbations in these cues alter cellular mechanisms of mechanotransduction, consequently enhancing malignancy. Finally, we elucidate engineering techniques to individually emulate the mechanical properties of tumors that could help serve as toolkits for developing and testing ECM-targeted therapeutics on novel bioengineered tumor platforms. STATEMENT OF SIGNIFICANCE: Disrupted ECM mechanics is a driving force for transitioning incipient cells to life-threatening malignant variants. Understanding these ECM changes can be crucial as they may aid in developing several efficacious drugs that not only focus on inducing cytotoxic effects but also target specific matrix mechanical cues that support and enhance tumor invasiveness. Designing and implementing an optimal tumor mimic can allow us to predictively map biophysical cue-modulated cell behaviors and facilitate the design of improved lab-grown tumor models with accurately controlled structural features. This review focuses on the abnormal changes within the ECM during tumorigenesis and its implications on tumor cell-matrix mechanoreciprocity. Additionally, it accentuates engineering approaches to produce ECM features of varying levels of complexity which is critical for improving the efficiency of current engineered tumor tissue models.

A microfluidic-informatics assay for quantitative physical occlusion measurement in sickle cell disease

Sickle cell disease (SCD) is a genetic condition that causes abnormalities in hemoglobin mechanics. Those affected are at high risk of vaso-occlusive crisis (VOC), which can induce life-threatening symptoms. The development of measurements related to vaso-occlusion facilitates the diagnosis of the patient's disease state. To complement existing readouts, we design a microfluidic-informatics analytical system with varied confined geometries for the quantification of sickle cell disease occlusion. We detect an increase in physical occlusion events in the most severe hemoglobin SS group. We use bioinformatics and modeling to quantify the in vitro disease severity score (DSS) of individual patients. We also show the potential effect of hydration, clinically recommended for crisis management, on reducing the disease severity of high-risk patients. Overall, we demonstrate the device as an easy-to-use assay for quick occlusion information extraction with a simple setup and minimal additional instruments. We show the device can provide physical readouts distinct from clinical data. We also show the device sensitivity in separate samples from patients with different disease severity. Finally, we demonstrate the system as a potential platform for testing the effectiveness of therapeutic strategies (e.g. hydration) on reducing sickle cell disease severity.

Biglycan Promotes Cancer Stem Cell Properties, NFκB Signaling and Metastatic Potential in Breast Cancer Cells

Simple Summary

Breast cancer stem cells (BCSCs) are a small sub-population of cells within tumors with high metastatic potential. We identified biglycan (BGN) as a prospective molecular target in BCSCs that regulates the aggressive phenotypes of these cells. These findings establish a foundation for the development of therapeutics against BGN to eliminate BCSCs and prevent metastatic breast cancer.

Abstract

It is a major challenge to treat metastasis due to the presence of heterogenous BCSCs. Therefore, it is important to identify new molecular targets and their underlying molecular mechanisms in various BCSCs to improve treatment of breast cancer metastasis. Here, we performed RNA sequencing on two distinct co-existing BCSC populations, ALDH⁺ and CD29^hi CD61⁺ from PyMT mammary tumor cells and detected upregulation of biglycan (BGN) in these BCSCs. Genetic depletion of BGN reduced BCSC proportions and tumorsphere formation. Furthermore, BCSC associated aggressive traits such as migration and invasion were significantly reduced by depletion of BGN. Glycolytic and mitochondrial metabolic assays also revealed that BCSCs exhibited decreased metabolism upon loss of BGN. BCSCs showed decreased activation of the NFκB transcription factor, p65, and phospho-IκB levels upon BGN ablation, indicating regulation of NFκB pathway by BGN. To further support our data, we also characterized CD24⁻/CD44⁺ BCSCs from human luminal MCF-7 breast cancer cells. These CD24⁻/CD44⁺ BCSCs similarly exhibited reduced tumorigenic phenotypes, metabolism and attenuation of NFκB pathway after knockdown of BGN. Finally, loss of BGN in ALDH⁺ and CD29^hi CD61⁺ BCSCs showed decreased metastatic potential, suggesting BGN serves as an important therapeutic target in BCSCs for treating metastasis of breast cancer.

Volumetric compression develops noise-driven single-cell heterogeneity

Recent studies have revealed that extensive heterogeneity of biological systems arises through various routes ranging from intracellular chromosome segregation to spatiotemporally varying biochemical stimulations. However, the contribution of physical microenvironments to single-cell heterogeneity remains largely unexplored. Here, we show that a homogeneous population of non-small-cell lung carcinoma develops into heterogeneous subpopulations upon application of a homogeneous physical compression, as shown by single-cell transcriptome profiling. The generated subpopulations stochastically gain the signature genes associated with epithelial-mesenchymal transition (EMT; VIM, CDH1, EPCAM, ZEB1, and ZEB2) and cancer stem cells (MKI67, BIRC5, and KLF4), respectively. Trajectory analysis revealed two bifurcated paths as cells evolving upon the physical compression, along each path the corresponding signature genes (epithelial or mesenchymal) gradually increase. Furthermore, we show that compression increases gene expression noise, which interplays with regulatory network architecture and thus generates differential cell-fate outcomes. The experimental observations of both single-cell sequencing and single-molecule fluorescent in situ hybridization agrees well with our computational modeling of regulatory network in the EMT process. These results demonstrate a paradigm of how mechanical stimulations impact cell-fate determination by altering transcription dynamics; moreover, we show a distinct path that the ecology and evolution of cancer interplay with their physical microenvironments from the view of mechanobiology and systems biology, with insight into the origin of single-cell heterogeneity.

A mesoscale mechanical model of cellular interactions

Computational models of cell mechanics allow the precise interrogation of cell shape change. These morphological changes are required for cells to survive in diverse tissue environments. Here, we present a mesoscale mechanical model of cell-substrate interactions using the level set method based on experimentally measured parameters. By implementing a viscoelastic mechanical equivalent circuit, we accurately model whole-cell deformations that are important for a variety of cellular processes. To effectively model shape changes as a cell interacts with a substrate, we have included receptor-mediated adhesion, which is governed by catch-slip bond behavior. The effect of adhesion was explored by subjecting cells to a variety of different substrates including flat, curved, and deformable surfaces. Finally, we increased the accuracy of our simulations by including a deformable nucleus in our cells. This model sets the foundation for further exploration into computational analyses of multicellular interactions.

Morphodynamic signatures of MDA-MB-231 single cells and cell doublets undergoing invasion in confined microenvironments

Cancer cell metastasis is a major factor in cancer-related mortality. During the process of metastasis, cancer cells exhibit migratory phenotypes and invade through pores in the dense extracellular matrix. However, the characterization of morphological and subcellular features of cells in similar migratory phenotypes and the effects of geometric confinement on cell morphodynamics are not well understood. Here, we investigate the phenotypes of highly aggressive MDA-MB-231 cells in single cell and cell doublet (an initial and simplified collective state) forms in confined microenvironments. We group phenotypically similar single cells and cell doublets and characterize related morphological and subcellular features. We further detect two distinct migratory phenotypes, fluctuating and non-fluctuating, within the fast migrating single cell group. In addition, we demonstrate an increase in the number of protrusions formed at the leading edge of cells after invasion through geometric confinement. Finally, we track the short and long term effects of varied degrees of confinement on protrusion formation. Overall, our findings elucidate the underlying morphological and subcellular features associated with different single cell and cell doublet phenotypes and the impact of invasion through confined geometry on cell behavior.