Microfluidic Ecology Unravels the Genetic and Ecological Drivers of T4r Bacteriophage Resistance in E. coli: Insights into Biofilm-Mediated Evolution.

We use a microfluidic ecology which generates non-uniform phage concentration gradients and micro-ecological niches to reveal the importance of time, spatial population structure and collective population dynamics in the de novo evolution of T4r bacteriophage resistant motile E. coli. An insensitive bacterial population against T4r phage occurs within 20 hours in small interconnected population niches created by a gradient of phage virions, driven by evolution in transient biofilm patches. Sequencing of the resistant bacteria reveals mutations at the receptor site of bacteriophage T4r as expected but also in genes associated with biofilm formation and surface adhesion, supporting the hypothesis that evolution within transient biofilms drives de novo phage resistance.

Lung adenocarcinomas without driver genes converge to common adaptive strategies through diverse genetic, epigenetic, and niche construction evolutionary pathways.

Somatic evolution selects cancer cell phenotypes that maximize survival and proliferation in dynamic environments. Although cancer cells are molecularly heterogeneous, we hypothesized convergent adaptive strategies to common host selection forces can be inferred from patterns of epigenetic and genetic evolutionary selection in similar tumors. We systematically investigated gene mutations and expression changes in lung adenocarcinomas with no common driver genes (n = 313). Although 13,461 genes were mutated in at least one sample, only 376 non-synonymous mutations evidenced positive evolutionary selection with conservation of 224 genes, while 1736 and 2430 genes exhibited ≥ two-fold increased and ≥ 50% decreased expression, respectively. Mutations under positive selection are more frequent in genes with significantly altered expression suggesting they often "hardwire" pre-existing epigenetically driven adaptations. Conserved genes averaged 16-fold higher expression in normal lung tissue compared to those with selected mutations demonstrating pathways necessary for both normal cell function and optimal cancer cell fitness. The convergent LUAD phenotype exhibits loss of differentiated functions and cell-cell interactions governing tissue organization. Conservation with increased expression is found in genes associated with cell cycle, DNA repair, p53 pathway, epigenetic modifiers, and glucose metabolism. No canonical driver gene pathways exhibit strong positive selection, but extensive down-regulation of membrane ion channels suggests decreased transmembrane potential may generate persistent proliferative signals. NCD LUADs perform niche construction generating a stiff, immunosuppressive microenvironment through selection of specific collagens and proteases. NCD LUADs evolve to a convergent phenotype through a network of interconnected genetic, epigenetic, and ecological pathways.

Exploration of drug resistance mechanisms in triple negative breast cancer cells using a microfluidic device and patient tissues.

Chemoresistance is a major cause of treatment failure in many cancers. However, the life cycle of cancer cells as they respond to and survive environmental and therapeutic stress is understudied. In this study, we utilized a microfluidic device to induce the development of doxorubicin-resistant (DOXR) cells from triple negative breast cancer (TNBC) cells within 11 days by generating gradients of DOX and medium. In vivo chemoresistant xenograft models, an unbiased genome-wide transcriptome analysis, and a patient data/tissue analysis all showed that chemoresistance arose from failed epigenetic control of the nuclear protein-1 (NUPR1)/histone deacetylase 11 (HDAC11) axis, and high NUPR1 expression correlated with poor clinical outcomes. These results suggest that the chip can rapidly induce resistant cells that increase tumor heterogeneity and chemoresistance, highlighting the need for further studies on the epigenetic control of the NUPR1/HDAC11 axis in TNBC.

The contribution of evolvability to the eco-evolutionary dynamics of competing species.

Evolvability is the capacity of a population to generate heritable variation that can be acted upon by natural selection. This ability influences the adaptations and fitness of individual organisms. By viewing this capacity as a trait, evolvability is subject to natural selection and thus plays a critical role in eco-evolutionary dynamics. Understanding this role provides insight into how species respond to changes in their environment and how species coexistence can arise and be maintained. Here, we create a G-function model of competing species, each with a different evolvability. We analyze population and strategy (= heritable phenotype) dynamics of the two populations under clade initiation (when species are introduced into a population), evolutionary tracking (constant, small changes in the environment), adaptive radiation (availability of multiple ecological niches), and evolutionary rescue (extreme environmental disturbances). We find that when species are far from an eco-evolutionary equilibrium, faster-evolving species reach higher population sizes, and when species are close to an equilibrium, slower-evolving species are more successful. Frequent, minor environmental changes promote the extinction of species with small population sizes, regardless of their evolvability. When several niches are available for a species to occupy, coexistence is possible, though slower-evolving species perform slightly better than faster-evolving ones due to the well-recognized inherent cost of evolvability. Finally, disrupting the environment at intermediate frequencies can result in coexistence with cyclical population dynamics of species with different rates of evolution.

A mathematical investigation of polyaneuploid cancer cell memory and cross-resistance in state-structured cancer populations.

The polyaneuploid cancer cell (PACC) state promotes cancer lethality by contributing to survival in extreme conditions and metastasis. Recent experimental evidence suggests that post-therapy PACC-derived recurrent populations display cross-resistance to classes of therapies with independent mechanisms of action. We hypothesize that this can occur through PACC memory, whereby cancer cells that have undergone a polyaneuploid transition (PAT) reenter the PACC state more quickly or have higher levels of innate resistance. In this paper, we build on our prior mathematical models of the eco-evolutionary dynamics of cells in the 2N+ and PACC states to investigate these two hypotheses. We show that although an increase in innate resistance is more effective at promoting cross-resistance, this trend can also be produced via PACC memory. We also find that resensitization of cells that acquire increased innate resistance through the PAT have a considerable impact on eco-evolutionary dynamics and extinction probabilities. This study, though theoretical in nature, can help inspire future experimentation to tease apart hypotheses surrounding how cross-resistance in structured cancer populations arises.

Escherichia coli survival in response to ciprofloxacin antibiotic stress correlates with increased nucleoid length and effective misfolded protein management.

The evolution of antibiotic resistance is a fundamental problem in disease management but is rarely quantified on a single-cell level owing to challenges associated with capturing the spatial and temporal variation across a population. To evaluate cell biological phenotypic responses, we tracked the single-cell dynamics of filamentous bacteria through time in response to ciprofloxacin antibiotic stress. We measured the degree of phenotypic variation in nucleoid length and the accumulation of protein damage under ciprofloxacin antibiotic and quantified the impact on bacterial survival. Increased survival was correlated with increased nucleoid length and the variation in this response was inversely correlated with antibiotic concentration. Survival time was also increased through clearance of misfolded proteins, an unexpected mechanism of stress relief deployed by the filamentous bacteria. Our results reveal a diverse range of survival tactics employed by bacteria in response to ciprofloxacin and suggest potential evolutionary routes to resistance.

Updating the Definition of Cancer.

Most definitions of cancer broadly conform to the current NCI definition: "Cancer is a disease in which some of the body's cells grow uncontrollably and spread to other parts of the body." These definitions tend to describe what cancer "looks like" or "does" but do not describe what cancer "is" or "has become." While reflecting past insights, current definitions have not kept pace with the understanding that the cancer cell is itself transformed and evolving. We propose a revised definition of cancer: Cancer is a disease of uncontrolled proliferation by transformed cells subject to evolution by natural selection. We believe this definition captures the essence of the majority of previous and current definitions. To the simplest definition of cancer as a disease of uncontrolled proliferation of cells, our definition adds in the adjective "transformed" to capture the many tumorigenic processes that cancer cells adopt to metastasize. To the concept of uncontrolled proliferation of transformed cells, our proposed definition then adds "subject to evolution by natural selection." The subject to evolution by natural selection modernizes the definition to include the genetic and epigenetic changes that accumulate within a population of cancer cells that lead to the lethal phenotype. Cancer is a disease of uncontrolled proliferation by transformed cells subject to evolution by natural selection.

Nuclear morphology predicts cell survival to cisplatin chemotherapy.

The emergence of chemotherapy resistance drives cancer lethality in cancer patients, with treatment initially reducing overall tumor burden followed by resistant recurrent disease. While molecular mechanisms underlying resistance phenotypes have been explored, less is known about the cell biological characteristics of cancer cells that survive to eventually seed the recurrence. To identify the unique phenotypic characteristics associated with survival upon chemotherapy exposure, we characterized nuclear morphology and function as prostate cancer cells recovered following cisplatin treatment. Cells that survived in the days and weeks after treatment and resisted therapy-induced cell death showed increasing cell size and nuclear size, enabled by continuous endocycling resulting in repeated whole genome doubling. We further found that cells that survive after therapy release were predominantly mononucleated and likely employ more efficient DNA damage repair. Finally, we show that surviving cancer cells exhibit a distinct nucleolar phenotype and increased rRNA levels. These data support a paradigm where soon after therapy release, the treated population mostly contains cells with a high level of widespread and catastrophic DNA damage that leads to apoptosis, while the minority of cells that have successful DDR are more likely to access a pro-survival state. These findings are consistent with accession of the polyaneuploid cancer cell (PACC) state, a recently described mechanism of therapy resistance and tumor recurrence. Our findings demonstrate the fate of cancer cells following cisplatin treatment and define key cell phenotypic characteristics of the PACC state. This work is essential for understanding and, ultimately, targeting cancer resistance and recurrence.

Modeling cancer's ecological and evolutionary dynamics.

In this didactic paper, we present a theoretical modeling framework, called the G-function, that integrates both the ecology and evolution of cancer to understand oncogenesis. The G-function has been used in evolutionary ecology, but has not been widely applied to problems in cancer. Here, we build the G-function framework from fundamental Darwinian principles and discuss how cancer can be seen through the lens of ecology, evolution, and game theory. We begin with a simple model of cancer growth and add on components of cancer cell competition and drug resistance. To aid in exploration of eco-evolutionary modeling with this approach, we also present a user-friendly software tool. By the end of this paper, we hope that readers will be able to construct basic G function models and grasp the usefulness of the framework to understand the games cancer plays in a biologically mechanistic fashion.

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter.

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

A life history model of the ecological and evolutionary dynamics of polyaneuploid cancer cells.

Therapeutic resistance is one of the main reasons for treatment failure in cancer patients. The polyaneuploid cancer cell (PACC) state has been shown to promote resistance by providing a refuge for cancer cells from the effects of therapy and by helping them adapt to a variety of environmental stressors. This state is the result of aneuploid cancer cells undergoing whole genome doubling and skipping mitosis, cytokinesis, or both. In this paper, we create a novel mathematical framework for modeling the eco-evolutionary dynamics of state-structured populations and use this framework to construct a model of cancer populations with an aneuploid and a PACC state. Using in silico simulations, we explore how the PACC state allows cancer cells to (1) survive extreme environmental conditions by exiting the cell cycle after S phase and protecting genomic material and (2) aid in adaptation to environmental stressors by increasing the cancer cell's ability to generate heritable variation (evolvability) through the increase in genomic content that accompanies polyploidization. In doing so, we demonstrate the ability of the PACC state to allow cancer cells to persist under therapy and evolve therapeutic resistance. By eliminating cells in the PACC state through appropriately-timed PACC-targeted therapies, we show how we can prevent the emergence of resistance and promote cancer eradication.

Stochastic models of Mendelian and reverse transcriptional inheritance in state-structured cancer populations.

Recent evidence suggests that a polyaneuploid cancer cell (PACC) state may play a key role in the adaptation of cancer cells to stressful environments and in promoting therapeutic resistance. The PACC state allows cancer cells to pause cell division and to avoid DNA damage and programmed cell death. Transition to the PACC state may also lead to an increase in the cancer cell's ability to generate heritable variation (evolvability). One way this can occur is through evolutionary triage. Under this framework, cells gradually gain resistance by scaling hills on a fitness landscape through a process of mutation and selection. Another way this can happen is through self-genetic modification whereby cells in the PACC state find a viable solution to the stressor and then undergo depolyploidization, passing it on to their heritably resistant progeny. Here, we develop a stochastic model to simulate both of these evolutionary frameworks. We examine the impact of treatment dosage and extent of self-genetic modification on eco-evolutionary dynamics of cancer cells with aneuploid and PACC states. We find that under low doses of therapy, evolutionary triage performs better whereas under high doses of therapy, self-genetic modification is favored. This study generates predictions for teasing apart these biological hypotheses, examines the implications of each in the context of cancer, and provides a modeling framework to compare Mendelian and non-traditional forms of inheritance.

Robots as models of evolving systems.

Experimental robobiological physics can bring insights into biological evolution. We present a development of hybrid analog/digital autonomous robots with mutable diploid dominant/recessive 6-byte genomes. The robots are capable of death, rebirth, and breeding. We map the quasi-steady-state surviving local density of the robots onto a multidimensional abstract "survival landscape." We show that robot death in complex, self-adaptive stress landscapes proceeds by a general lowering of the robotic genetic diversity, and that stochastically changing landscapes are the most difficult to survive.

It doesn't always pay to be fit: success landscapes.

Landscapes play an important role in many areas of biology, in which biological lives are deeply entangled. Here we discuss a form of landscape in evolutionary biology which takes into account (1) initial growth rates, (2) mutation rates, (3) resource consumption by organisms, and (4) cyclic changes in the resources with time. The long-term equilibrium number of surviving organisms as a function of these four parameters forms what we call a success landscape, a landscape we would claim is qualitatively different from fitness landscapes which commonly do not include mutations or resource consumption/changes in mapping genomes to the final number of survivors. Although our analysis is purely theoretical, we believe the results have possibly strong connections to how we might treat diseases such as cancer in the future with a deeper understanding of the interplay between resource degradation, mutation, and uncontrolled cell growth.

Protein dynamics implications of the low- and high-temperature denaturation of myoglobin.

We reinvestigate a simple model used in the literature concerning the thermodynamic analysis of protein cold denaturation. We derive an exact thermodynamic expression for cold denaturation and give a better approximation than exists in the literature for predicting cold denaturation temperatures in the two-state model. We discuss the "dark-side" implications of this work for previous temperature-dependent protein dynamics experiments and discuss microfluidic experimental technologies, which could explore the thermal stability range of proteins below the bulk freezing point of water.

Emergent Field-Driven Robot Swarm States.

We present an ecology-inspired form of active matter consisting of a robot swarm. Each robot moves over a planar dynamic resource environment represented by a large light-emitting diode array in search of maximum light intensity; the robots deplete (dim) locally by their presence the local light intensity and seek maximum light intensity. Their movement is directed along the steepest local light intensity gradient; we call this emergent symmetry breaking motion "field drive." We show there emerge dynamic and spatial transitions similar to gas, crystalline, liquid, glass, and jammed states as a function of robot density, resource consumption rates, and resource recovery rates. Paradoxically the nongas states emerge from smooth, flat resource landscapes, not rough ones, and each state can directly move to a glassy state if the resource recovery rate is slow enough, at any robot density.

Scaling of deterministic lateral displacement devices to a single column of bumping obstacles.

We describe a deterministic lateral displacement (DLD) for particle separation with only a single column of bumping features. The bifurcation of fluid streams at obstacles is not set by the "tilt" of columns with respect to macroscopic current flow, but rather by the fluidic resistances for lateral flow at each obstacle. With one column of 14 bumping features and corresponding inlet/outlet channels, the single-column DLD can separate particles with diameters of 4.8 µm and 9.9 µm at 30 µL min-1, with an area of only 0.37 mm × 1.5 mm (0.55 mm2). The large-cell output contains over 99% of the 9.9 µm particles and only 0.2% of the 4.8 m particles. The throughput per area of 54 µL min-1 per mm2 represents a 10× increase over previous selective harvesting reports for microfluidic devices in a similar particle size range.

Deterministic Lateral Displacement: Challenges and Perspectives.

The advent of microfluidics in the 1990s promised a revolution in multiple industries from healthcare to chemical processing. Deterministic lateral displacement (DLD) is a continuous-flow microfluidic particle separation method discovered in 2004 that has been applied successfully and widely to the separation of blood cells, yeast, spores, bacteria, viruses, DNA, droplets, and more. Deterministic lateral displacement is conceptually simple and can deliver consistent performance over a wide range of flow rates and particle concentrations. Despite wide use and in-depth study, DLD has not yet been fully elucidated or optimized, with different approaches to the same problem yielding varying results. We endeavor here to provide up-to-date expert opinion on the state-of-art and current fundamental, practical, and commercial challenges with DLD as well as describe experimental and modeling opportunities. Because these challenges and opportunities arise from constraints on hydrodynamics, fabrication, and operation at the micro- and nanoscale, we expect this Perspective to serve as a guide for the broader micro- and nanofluidic community to identify and to address open questions in the field.

An in vitro tumor swamp model of heterogeneous cellular and chemotherapeutic landscapes.

The heterogenous, highly metabolic stressed, poorly irrigated, solid tumor microenvironment - the tumor swamp - is widely recognized to play an important role in cancer progression as well as the development of therapeutic resistance. It is thus important to create realistic in vitro models within the therapeutic pipeline that can recapitulate the fundamental stress features of the tumor swamp. Here we describe a microfluidic system which generates a chemical gradient within connected microenvironments achieved through a static diffusion mechanism rather than active pumping. We show that the gradient can be stably maintained for over a week. Due to the accessibility and simplicity of the experimental platform, the system allows for not only well-controlled continuous studies of the interactions among various cell types at single-cell resolution, but also parallel experimentation for time-resolved downstream cellular assays on the time scale of weeks. This approach enables simple, compact implementation and is compatible with existing 6-well imaging technology for simultaneous experiments. As a proof-of-concept, we report the co-culture of a human bone marrow stromal cell line and a bone-metastatic prostate cancer cell line using the presented device, revealing on the same chip a transition in cancer cell survival as a function of drug concentration on the population level while exhibiting an enrichment of poly-aneuploid cancer cells (PACCs) as an evolutionary consequence of high stress. The device allows for the quantitative study of cancer cell dynamics on a stress landscape by real-time monitoring of various cell types with considerable experimental throughput.