We use a microfabricated ecology with a doxorubicin gradient and population fragmentation to produce a strong Darwinian selective pressure that drives forward the rapid emergence of doxorubicin resistance in multiple myeloma (MM) cancer cells. RNA sequencing of the resistant cells was used to examine (i) emergence of genes with high de novo substitution densities (i.e., hot genes) and (ii) genes never substituted (i.e., cold genes). The set of cold genes, which were 21% of the genes sequenced, were further winnowed down by examining excess expression levels. Both the most highly substituted genes and the most highly expressed never-substituted genes were biased in age toward the most ancient of genes. This would support the model that cancer represents a revision back to ancient forms of life adapted to high fitness under extreme stress, and suggests that these ancient genes may be targets for cancer therapy.
Bacteria can rapidly evolve resistance to antibiotics via the SOS response, a state of high-activity DNA repair and mutagenesis. We explore here the first steps of this evolution in the bacterium Escherichia coli. Induction of the SOS response by the genotoxic antibiotic ciprofloxacin changes the E. coli rod shape into multichromosome-containing filaments. We show that at subminimal inhibitory concentrations of ciprofloxacin the bacterial filament divides asymmetrically repeatedly at the tip. Chromosome-containing buds are made that, if resistant, propagate nonfilamenting progeny with enhanced resistance to ciprofloxacin as the parent filament dies. We propose that the multinucleated filament creates an environmental niche where evolution can proceed via generation of improved mutant chromosomes due to the mutagenic SOS response and possible recombination of the new alleles between chromosomes. Our data provide a better understanding of the processes underlying the origin of resistance at the single-cell level and suggest an analogous role to the eukaryotic aneuploidy condition in cancer.
We demonstrate that a microfabricated bump array can concentrate genomic-length DNA molecules efficiently at continuous, high flow velocities, up to 40 μm/s, if the single-molecule DNA globule has a sufficiently large shear modulus. Increase in the shear modulus is accomplished by compacting the DNA molecules to minimal coil size using polyethylene glycol (PEG) derived depletion forces. We map out the sweet spot, where concentration occurs, as a function of PEG concentration and flow speed using a combination of theoretical analysis and experiment. Purification of DNA from enzymatic reactions for next-generation DNA-sequencing libraries will be an important application of this development.
Microfluidic deterministic lateral displacement (DLD) arrays have been applied for fractionation and analysis of cells in quantities of ~100 μL of blood, with processing of larger quantities limited by clogging in the chip. In this paper, we (i) demonstrate that this clogging phenomenon is due to conventional platelet-driven clot formation, (ii) identify and inhibit the two dominant biological mechanisms driving this process, and (iii) characterize how further reductions in clot formation can be achieved through higher flow rates and blood dilution. Following from these three advances, we demonstrate processing of 14 mL equivalent volume of undiluted whole blood through a single DLD array in 38 minutes to harvest PC3 cancer cells with ~86% yield. It is possible to fit more than 10 such DLD arrays on a single chip, which would then provide the capability to process well over 100 mL of undiluted whole blood on a single chip in less than one hour.
We describe a microfluidic device for on-chip chemical processing, such as staining, and subsequent washing of cells. The paper introduces "separator walls" to increase the on-chip incubation time and to improve the quality of washing. Cells of interest are concentrated into a treatment stream of chemical reagents at the first separator wall for extended on-chip incubation without causing excess contamination at the output due to diffusion of the unreacted treatment chemicals, and then are directed to the washing stream before final collections. The second separator wall further reduces the output contamination from diffusion to the washing stream. With this approach, we demonstrate on-chip leukocyte staining with Rhodamine 6G and washing. The results suggest that other conventional biological and analytical processes could be replaced by the proposed device.
One of the predictions of game theory is that cooperative behaviours are vulnerable to exploitation by selfish individuals, but this result seemingly contradicts the survival of cooperation observed in nature. In this review, we will introduce game theoretical concepts that lead to this conclusion and show how the spatial competition dynamics between microorganisms can be used to model the survival and maintenance of cooperation. In particular, we focus on how Escherichia coli bacteria with a growth advantage in stationary phase (GASP) phenotype maintain a proliferative phenotype when faced with overcrowding to gain a fitness advantage over wild-type populations. We review recent experimental approaches studying the growth dynamics of competing GASP and wild-type strains of E. coli inside interconnected microfabricated habitats and use a game theoretical approach to analyse the observed inter-species interactions. We describe how the use of evolutionary game theory and the ideal free distribution accurately models the spatial distribution of cooperative and selfish individuals in spatially heterogeneous environments. Using bacteria as a model system of cooperative and selfish behaviours may lead to a better understanding of the competition dynamics of other organisms-including tumour-host interactions during cancer development and metastasis.
Preventing relapse is the major challenge to effective therapy in cancer. Within the tumour, stromal (ST) cells play an important role in cancer progression and the emergence of drug resistance. During cancer treatment, the fitness of cancer cells can be enhanced by ST cells because their molecular signalling interaction delays the drug-induced apoptosis of cancer cells. On the other hand, competition among cancer and ST cells for space or resources should not be ignored. We explore the population dynamics of multiple myeloma (MM) versus bone marrow ST cells by using an experimental microecology that we call the death galaxy, with a stable drug gradient and connected microhabitats. Evolutionary game theory is a quantitative way to capture the frequency-dependent nature of interactive populations. Therefore, we use evolutionary game theory to model the populations in the death galaxy with the gradients of pay-offs and successfully predict the future densities of MM and ST cells. We discuss the possible clinical use of such analysis for predicting cancer progression.
Drug development faces its nemesis in the form of drug resistance. The rate of bacterial resistance to antibiotics, or tumor resistance to chemotherapy decisively depends on the surrounding heterogeneous tissue. However, in vitro drug testing is almost exclusively done in well stirred, homogeneous environments. Recent advancements in microfluidics and microfabrication introduce opportunities to develop in vitro culture models that mimic the complex in vivo tissue environment. In this review, we will first discuss the design principles underlying such models. Then we will demonstrate two types of microfluidic devices that combine stressor gradients, cell motility, large population of competing/cooperative cells and time varying dosage of drugs. By incorporating ideas from how natural selection and evolution move drug resistance forward, we show that drug resistance can occur at much greater rates than in well-stirred environments. Finally, we will discuss the future direction of in vitro microbial culture models and how to extend the lessons learned from microbial systems to eukaryotic cells.
The transport of objects in microfluidic arrays of obstacles is a surprisingly rich area of physics and statistical mechanics. Tom Duke's mastery of these areas had a major impact in the development of biotechnology which uses these ideas at an increasing scale. We first review how biological objects are transported in fluids at low Reynolds numbers, including a discussion of electrophoresis, then concentrate on the separation of objects in asymmetric arrays, sometimes called Brownian ratchets when diffusional symmetry is broken by the structures. We move beyond this to what are called deterministic arrays where non-hydrodynamic forces in asymmetric arrays allow for extraordinary separation, and we look to the future of using these unusual arrays at the nanoscale and at the hundreds of micrometre scale. The emphasis is on how the original ideas of Tom Duke drove this work forward.
Political correctness urges us to state how wonderful it is to work with biologists and how, just as the lion will someday lie down with the lamb, so will interdisciplinary work, where biologists and physicists are mixed together in light, airy buildings designed to force socialization, give rise to wonderful new science. But it has been said that the only drive in human nature stronger than the sex drive is the drive to censor and suppress, and so I claim that it is OK for physicists and biologists to maintain a wary distance from each other, so that neither one censors or suppresses the wild ideas of the other.
The issue of resistance to targeted drug therapy is of pressing concern, as it constitutes a major barrier to progress in managing cancer. One important aspect is the role of stochasticity in determining the nature of the patient response. We examine two particular experiments. The first measured the maximal response of melanoma to targeted therapy before the resistance causes the tumor to progress. We analyze the data in the context of a Delbruck-Luria type scheme, wherein the continued growth of preexistent resistant cells are responsible for progression. We show that, aside from a finite fraction of resistant cell-free patients, the maximal response in such a scenario would be quite uniform. To achieve the measured variability, one is necessarily led to assume a wide variation from patient to patient of the sensitive cells' response to the therapy. The second experiment is an in vitro system of multiple myeloma cells. When subject to a spatial gradient of a chemotherapeutic agent, the cells in the middle of the system acquire resistance on a rapid (two-week) timescale. This finding points to the potential important role of cell-to-cell differences, due to differing local environments, in addition to the patient-to-patient differences encountered in the first part. See all articles in this Cancer Research section, "Physics in Cancer Research."
Do genetically closely related organisms under identical, but strong selection pressure converge to a common resistant genotype or will they diverge to different genomic solutions? This question gets at the heart of how rough is the fitness landscape in the local vicinity of two closely related strains under stress. We chose a Growth Advantage in Stationary Phase (GASP) E scherichia coli strain to address this question because the GASP strain has very similar fitness to the wild-type (WT) strain in the absence of metabolic stress but in the presence of metabolic stress continues to divide and does not enter into stationary phase. We find that under strong antibiotic selection pressure by the fluoroquinolone antibiotic ciprofloxacin in a complex ecology that the GASP strain rapidly evolves in under 20 h missense mutation in gyrA only 2 amino acids removed from the WT strain indicating a convergent solution, yet does not evolve the other 3 mutations of the WT strain. Further the GASP strain evolves a prophage e14 excision which completely inhibits biofilm formation in the mutant strain, revealing the hidden complexity of E. coli evolution to antibiotics as a function of selection pressure. We conclude that there is a cryptic roughness to fitness landscapes in the absence of stress.
The emergence of resistance to chemotherapy by cancer cells, when combined with metastasis, is the primary driver of mortality in cancer and has proven to be refractory to many efforts. Theory and computer modeling suggest that the rate of emergence of resistance is driven by the strong selective pressure of mutagenic chemotherapy and enhanced by the motility of mutant cells in a chemotherapy gradient to areas of higher drug concentration and lower population competition. To test these models, we constructed a synthetic microecology which superposed a mutagenic doxorubicin gradient across a population of motile, metastatic breast cancer cells (MDA-MB-231). We observed the emergence of MDA-MB-231 cancer cells capable of proliferation at 200 nM doxorubicin in this complex microecology. Individual cell tracking showed both movement of the MDA-MB-231 cancer cells toward higher drug concentrations and proliferation of the cells at the highest doxorubicin concentrations within 72 h, showing the importance of both motility and drug gradients in the emergence of resistance.
Metastasis, the truly lethal aspect of cancer, occurs when metastatic cancer cells in a tumor break through the basement membrane and penetrate the extracellular matrix. We show that MDA-MB-231 metastatic breast cancer cells cooperatively invade a 3D collagen matrix while following a glucose gradient. The invasion front of the cells is a dynamic one, with different cells assuming the lead on a time scale of 70 h. The front cell leadership is dynamic presumably because of metabolic costs associated with a long-range strain field that precedes the invading cell front, which we have imaged using confocal imaging and marker beads imbedded in the collagen matrix. We suggest this could be a quantitative assay for an invasive phenotype tracking a glucose gradient and show that the invading cells act in a cooperative manner by exchanging leaders in the invading front.
We report an integrated nanochannel/nanoelectrode sensor for the detection of DNA using alternating currents. We find that DNA can be detected using platinum as the metal for the detecting electrodes, with a signal to noise ratio exceeding 10. We argue that the signal is at least in part electrochemical in nature, thus holds the promise to yield a sequence-dependent signal. However, we also find that for large voltages, DNA attaches irreversibly to the driving electrodes.
To investigate the transition from non-cancerous to metastatic from a physical sciences perspective, the Physical Sciences-Oncology Centers (PS-OC) Network performed molecular and biophysical comparative studies of the non-tumorigenic MCF-10A and metastatic MDA-MB-231 breast epithelial cell lines, commonly used as models of cancer metastasis. Experiments were performed in 20 laboratories from 12 PS-OCs. Each laboratory was supplied with identical aliquots and common reagents and culture protocols. Analyses of these measurements revealed dramatic differences in their mechanics, migration, adhesion, oxygen response, and proteomic profiles. Model-based multi-omics approaches identified key differences between these cells' regulatory networks involved in morphology and survival. These results provide a multifaceted description of cellular parameters of two widely used cell lines and demonstrate the value of the PS-OC Network approach for integration of diverse experimental observations to elucidate the phenotypes associated with cancer metastasis.
Circulating tumor cells (CTCs) and circulating clusters of cancer and stromal cells have been identified in the blood of patients with malignant cancer and can be used as a diagnostic for disease severity, assess the efficacy of different treatment strategies and possibly determine the eventual location of metastatic invasions for possible treatment. There is thus a critical need to isolate, propagate and characterize viable CTCs and clusters of cancer cells with their associated stroma cells. Here, we present a microfluidic device for mL/min flow rate, continuous-flow capture of viable CTCs from blood using deterministic lateral displacement (DLD) arrays. We show here that a DLD array device can isolate CTCs from blood with capture efficiency greater than 85% CTCs at volumetric flow rates of up to 10 mL/min with no effect on cell viability.
DNA is the central storage molecule of genetic information in the cell, and reading that information is a central problem in biology. While sequencing technology has made enormous advances over the past decade, there is growing interest in platforms that can readout genetic information directly from long single DNA molecules, with the ultimate goal of single-cell, single-genome analysis. Such a capability would obviate the need for ensemble averaging over heterogeneous cellular populations and eliminate uncertainties introduced by cloning and molecular amplification steps (thus enabling direct assessment of the genome in its native state). In this review, we will discuss how the information contained in genomic-length single DNA molecules can be accessed via physical confinement in nanochannels. Due to self-avoidance interactions, DNA molecules will stretch out when confined in nanochannels, creating a linear unscrolling of the genome along the channel for analysis. We will first review the fundamental physics of DNA nanochannel confinement--including the effect of varying ionic strength--and then discuss recent applications of these systems to genomic mapping. Apart from the intense biological interest in extracting linear sequence information from elongated DNA molecules, from a physics view these systems are fascinating as they enable probing of single-molecule conformation in environments with dimensions that intersect key physical length-scales in the 1 nm to 100 µm range.