Mutation and cell state compatibility is required and targetable in Ph+ acute lymphoblastic leukemia minimal residual disease.

Efforts to cure BCR::ABL1 B cell acute lymphoblastic leukemia (Ph+ ALL) solely through inhibition of ABL1 kinase activity have thus far been insufficient despite the availability of tyrosine kinase inhibitors (TKIs) with broad activity against resistance mutants. The mechanisms that drive persistence within minimal residual disease (MRD) remain poorly understood and therefore untargeted. Utilizing 13 patient-derived xenograft (PDX) models and clinical trial specimens of Ph+ ALL, we examined how genetic and transcriptional features co-evolve to drive progression during prolonged TKI response. Our work reveals a landscape of cooperative mutational and transcriptional escape mechanisms that differ from those causing resistance to first generation TKIs. By analyzing MRD during remission, we show that the same resistance mutation can either increase or decrease cellular fitness depending on transcriptional state. We further demonstrate that directly targeting transcriptional state-associated vulnerabilities at MRD can overcome BCR::ABL1 independence, suggesting a new paradigm for rationally eradicating MRD prior to relapse. Finally, we illustrate how cell mass measurements of leukemia cells can be used to rapidly monitor dominant transcriptional features of Ph+ ALL to help rationally guide therapeutic selection from low-input samples.

Cellular Mass Response to Therapy Correlates With Clinical Response for a Range of Malignancies.

PURPOSE: Cancer patients with advanced-stage disease have poor prognosis, typically having limited options for efficacious treatment, and genomics-based therapy guidance continues to benefit only a fraction of patients. Next-generation ex vivo approaches, such as cell mass-based response testing (MRT), offer an alternative precision medicine approach for a broader population of patients with cancer, but validation of clinical feasibility and potential impact remain necessary. MATERIALS AND METHODS: We evaluated the clinical feasibility and accuracy of using live-cell MRT to predict patient drug sensitivity. Using a unified measurement workflow with a 48-hour result turnaround time, samples were subjected to MRT after treatment with a panel of drugs in vitro. After completion of therapeutic course, clinical response data were correlated with MRT-based predictions of outcome. Specimens were collected from 104 patients with solid (n = 69) and hematologic (n = 35) malignancies, using tissue formats including needle biopsies, malignant fluids, bone marrow aspirates, and blood samples. Of the 81 (78%) specimens qualified for MRT, 41 (51%) patients receiving physician-selected therapies had treatments matched to MRT. RESULTS: MRT demonstrated high concordance with clinical responses with an odds ratio (OR) of 14.80 (P = .0003 [95% CI, 2.83 to 102.9]). This performance held for both solid and hematologic malignances with ORs of 20.67 (P = .0128 [95% CI, 1.45 to 1,375.57]) and 8.20 (P = .045 [95% CI, 0.77 to 133.56]), respectively. Overall, these results had a predictive accuracy of 80% (P = .0026 [95% CI, 65 to 91]). CONCLUSION: MRT showed highly significant correlation with clinical response to therapy. Routine clinical use is technically feasible and broadly applicable to a wide range of samples and malignancy types, supporting the need for future validation studies.

A pipeline for malignancy and therapy agnostic assessment of cancer drug response using cell mass measurements.

Functional precision medicine offers a promising complement to genomics-based cancer therapy guidance by testing drug efficacy directly on a patient's tumor cells. Here, we describe a workflow that utilizes single-cell mass measurements with inline brightfield imaging and machine-learning based image classification to broaden the clinical utility of such functional testing for cancer. Using these image-curated mass measurements, we characterize mass response signals for 60 different drugs with various mechanisms of action across twelve different cell types, demonstrating an improved ability to detect response for several slow acting drugs as compared with standard cell viability assays. Furthermore, we use this workflow to assess drug responses for various primary tumor specimen formats including blood, bone marrow, fine needle aspirates (FNA), and malignant fluids, all with reports generated within two days and with results consistent with patient clinical responses. The combination of high-resolution measurement, broad drug and malignancy applicability, and rapid return of results offered by this workflow suggests that it is well-suited to performing clinically relevant functional assessment of cancer drug response.

Mass measurements during lymphocytic leukemia cell polyploidization decouple cell cycle- and cell size-dependent growth.

Cell size is believed to influence cell growth and metabolism. Consistently, several studies have revealed that large cells have lower mass accumulation rates per unit mass (i.e., growth efficiency) than intermediate-sized cells in the same population. Size-dependent growth is commonly attributed to transport limitations, such as increased diffusion timescales and decreased surface-to-volume ratio. However, separating cell size- and cell cycle-dependent growth is challenging. To address this, we monitored growth efficiency of pseudodiploid mouse lymphocytic leukemia cells during normal proliferation and polyploidization. This was enabled by the development of large-channel suspended microchannel resonators that allow us to monitor buoyant mass of single cells ranging from 40 pg (small pseudodiploid cell) to over 4,000 pg, with a resolution ranging from ∼1% to ∼0.05%. We find that cell growth efficiency increases, plateaus, and then decreases as cell cycle proceeds. This growth behavior repeats with every endomitotic cycle as cells grow into polyploidy. Overall, growth efficiency changes 33% throughout the cell cycle. In contrast, increasing cell mass by over 100-fold during polyploidization did not change growth efficiency, indicating exponential growth. Consistently, growth efficiency remained constant when cell cycle was arrested in G2 Thus, cell cycle is a primary determinant of growth efficiency. As growth remains exponential over large size scales, our work finds no evidence for transport limitations that would decrease growth efficiency.

Rapid and high-precision sizing of single particles using parallel suspended microchannel resonator arrays and deconvolution.

Measuring the size of micron-scale particles plays a central role in the biological sciences and in a wide range of industrial processes. A variety of size parameters, such as particle diameter, volume, and mass, can be measured using electrical and optical techniques. Suspended microchannel resonators (SMRs) are microfluidic devices that directly measure particle mass by detecting a shift in resonance frequency as particles flow through a resonating microcantilever beam. While these devices offer high precision for sizing particles by mass, throughput is fundamentally limited by the small dimensions of the resonator and the limited bandwidth with which changes in resonance frequency can be tracked. Here, we introduce two complementary technical advancements that vastly increase the throughput of SMRs. First, we describe a deconvolution-based approach for extracting mass measurements from resonance frequency data, which allows an SMR to accurately measure a particle's mass approximately 16-fold faster than previously possible, increasing throughput from 120 particles/min to 2000 particles/min for our devices. Second, we describe the design and operation of new devices containing up to 16 SMRs connected fluidically in parallel and operated simultaneously on the same chip, increasing throughput to approximately 6800 particles/min without significantly degrading precision. Finally, we estimate that future systems designed to combine both of these techniques could increase throughput by nearly 200-fold compared to previously described SMR devices, with throughput potentially as high as 24 000 particles/min. We envision that increasing the throughput of SMRs will broaden the range of applications for which mass-based particle sizing can be employed.

Quantification of somatic mutation flow across individual cell division events by lineage sequencing.

Mutation data reveal the dynamic equilibrium between DNA damage and repair processes in cells and are indispensable to the understanding of age-related diseases, tumor evolution, and the acquisition of drug resistance. However, available genome-wide methods have a limited ability to resolve rare somatic variants and the relationships between these variants. Here, we present lineage sequencing, a new genome sequencing approach that enables somatic event reconstruction by providing quality somatic mutation call sets with resolution as high as the single-cell level in subject lineages. Lineage sequencing entails sampling single cells from a population and sequencing subclonal sample sets derived from these cells such that knowledge of relationships among the cells can be used to jointly call variants across the sample set. This approach integrates data from multiple sequence libraries to support each variant and precisely assigns mutations to lineage segments. We applied lineage sequencing to a human colon cancer cell line with a DNA polymerase epsilon (POLE) proofreading deficiency (HT115) and a human retinal epithelial cell line immortalized by constitutive telomerase expression (RPE1). Cells were cultured under continuous observation to link observed single-cell phenotypes with single-cell mutation data. The high sensitivity, specificity, and resolution of the data provide a unique opportunity for quantitative analysis of variation in mutation rate, spectrum, and correlations among variants. Our data show that mutations arrive with nonuniform probability across sublineages and that DNA lesion dynamics may cause strong correlations between certain mutations.

Linking single-cell measurements of mass, growth rate, and gene expression.

Mass and growth rate are highly integrative measures of cell physiology not discernable via genomic measurements. Here, we introduce a microfluidic platform enabling direct measurement of single-cell mass and growth rate upstream of highly multiplexed single-cell profiling such as single-cell RNA sequencing. We resolve transcriptional signatures associated with single-cell mass and growth rate in L1210 and FL5.12 cell lines and activated CD8+ T cells. Further, we demonstrate a framework using these linked measurements to characterize biophysical heterogeneity in a patient-derived glioblastoma cell line with and without drug treatment. Our results highlight the value of coupled phenotypic metrics in guiding single-cell genomics.

Trisomy of a Down Syndrome Critical Region Globally Amplifies Transcription via HMGN1 Overexpression.

Down syndrome (DS, trisomy 21) is associated with developmental abnormalities and increased leukemia risk. To reconcile chromatin alterations with transcriptome changes, we performed paired exogenous spike-in normalized RNA and chromatin immunoprecipitation sequencing in DS models. Absolute normalization unmasks global amplification of gene expression associated with trisomy 21. Overexpression of the nucleosome binding protein HMGN1 (encoded on chr21q22) recapitulates transcriptional changes seen with triplication of a Down syndrome critical region on distal chromosome 21, and HMGN1 is necessary for B cell phenotypes in DS models. Absolute exogenous-normalized chromatin immunoprecipitation sequencing (ChIP-Rx) also reveals a global increase in histone H3K27 acetylation caused by HMGN1. Transcriptional amplification downstream of HMGN1 is enriched for stage-specific programs of B cells and B cell acute lymphoblastic leukemia, dependent on the developmental cellular context. These data offer a mechanistic explanation for DS transcriptional patterns and suggest that further study of HMGN1 and RNA amplification in diverse DS phenotypes is warranted.

Microfluidic active loading of single cells enables analysis of complex clinical specimens.

A fundamental trade-off between flow rate and measurement precision limits performance of many single-cell detection strategies, especially for applications that require biophysical measurements from living cells within complex and low-input samples. To address this, we introduce 'active loading', an automated, optically-triggered fluidic system that improves measurement throughput and robustness by controlling entry of individual cells into a measurement channel. We apply active loading to samples over a range of concentrations (1-1000 particles µL-1), demonstrate that measurement time can be decreased by up to 20-fold, and show theoretically that performance of some types of existing single-cell microfluidic devices can be improved by implementing active loading. Finally, we demonstrate how active loading improves clinical feasibility for acute, single-cell drug sensitivity measurements by deploying it to a preclinical setting where we assess patient samples from normal brain, primary and metastatic brain cancers containing a complex, difficult-to-measure mixture of confounding biological debris.

Microfluidic platform for characterizing TCR-pMHC interactions.

The physical characteristics of the T cell receptor (TCR)-peptide-major histocompatibility complex (pMHC) interaction are known to play a central role in determining T cell function in the initial stages of the adaptive immune response. State-of-the-art assays can probe the kinetics of this interaction with single-molecular-bond resolution, but this precision typically comes at the cost of low throughput, since the complexity of these measurements largely precludes "scaling up." Here, we explore the feasibility of detecting specific TCR-pMHC interactions by flowing T cells past immobilized pMHC and measuring the reduction in cell speed due to the mechanical force of the receptor-ligand interaction. To test this new fluidic measurement modality, we fabricated a microfluidic device in which pMHC-coated beads are immobilized in hydrodynamic traps along the length of a serpentine channel. As T cells flow past the immobilized beads, their change in speed is tracked via microscopy. We validated this approach using two model systems: primary CD8+ T cells from an OT-1 TCR transgenic mouse with beads conjugated with H-2Kb:SIINFEKL, and Jurkat T cells with beads conjugated with anti-CD3 and anti-CD28 antibodies.

Determining therapeutic susceptibility in multiple myeloma by single-cell mass accumulation.

Multiple myeloma (MM) has benefited from significant advancements in treatment that have improved outcomes and reduced morbidity. However, the disease remains incurable and is characterized by high rates of drug resistance and relapse. Consequently, methods to select the most efficacious therapy are of great interest. Here we utilize a functional assay to assess the ex vivo drug sensitivity of single multiple myeloma cells based on measuring their mass accumulation rate (MAR). We show that MAR accurately and rapidly defines therapeutic susceptibility across human multiple myeloma cell lines to a gamut of standard-of-care therapies. Finally, we demonstrate that our MAR assay, without the need for extended culture ex vivo, correctly defines the response of nine patients to standard-of-care drugs according to their clinical diagnoses. This data highlights the MAR assay in both research and clinical applications as a promising tool for predicting therapeutic response using clinical samples.

Drug sensitivity of single cancer cells is predicted by changes in mass accumulation rate.

Assays that can determine the response of tumor cells to cancer therapeutics could greatly aid the selection of drug regimens for individual patients. However, the utility of current functional assays is limited, and predictive genetic biomarkers are available for only a small fraction of cancer therapies. We found that the single-cell mass accumulation rate (MAR), profiled over many hours with a suspended microchannel resonator, accurately defined the drug sensitivity or resistance of glioblastoma and B-cell acute lymphocytic leukemia cells. MAR revealed heterogeneity in drug sensitivity not only between different tumors, but also within individual tumors and tumor-derived cell lines. MAR measurement predicted drug response using samples as small as 25 µl of peripheral blood while maintaining cell viability and compatibility with downstream characterization. MAR measurement is a promising approach for directly assaying single-cell therapeutic responses and for identifying cellular subpopulations with phenotypic resistance in heterogeneous tumors.

High-throughput measurement of single-cell growth rates using serial microfluidic mass sensor arrays.

Methods to rapidly assess cell growth would be useful for many applications, including drug susceptibility testing, but current technologies have limited sensitivity or throughput. Here we present an approach to precisely and rapidly measure growth rates of many individual cells simultaneously. We flow cells in suspension through a microfluidic channel with 10-12 resonant mass sensors distributed along its length, weighing each cell repeatedly over the 4-20 min it spends in the channel. Because multiple cells traverse the channel at the same time, we obtain growth rates for >60 cells/h with a resolution of 0.2 pg/h for mammalian cells and 0.02 pg/h for bacteria. We measure the growth of single lymphocytic cells, mouse and human T cells, primary human leukemia cells, yeast, Escherichia coli and Enterococcus faecalis. Our system reveals subpopulations of cells with divergent growth kinetics and enables assessment of cellular responses to antibiotics and antimicrobial peptides within minutes.

Biophysical changes reduce energetic demand in growth factor-deprived lymphocytes.

Cytokine regulation of lymphocyte growth and proliferation is essential for matching nutrient consumption with cell state. Here, we examine how cellular biophysical changes that occur immediately after growth factor depletion promote adaptation to reduced nutrient uptake. After growth factor withdrawal, nutrient uptake decreases, leading to apoptosis. Bcl-xL expression prevents cell death, with autophagy facilitating long-term cell survival. However, autophagy induction is slow relative to the reduction of nutrient uptake, suggesting that cells must engage additional adaptive mechanisms to respond initially to growth factor depletion. We describe an acute biophysical response to growth factor withdrawal, characterized by a simultaneous decrease in cell volume and increase in cell density, which occurs before autophagy initiation and is observed in both FL5.12 Bcl-xL cells depleted of IL-3 and primary CD8(+) T cells depleted of IL-2 that are differentiating toward memory cells. The response reduces cell surface area to minimize energy expenditure while conserving biomass, suggesting that the biophysical properties of cells can be regulated to promote survival under conditions of nutrient stress.

A microfluidic platform enabling single-cell RNA-seq of multigenerational lineages.

We introduce a microfluidic platform that enables off-chip single-cell RNA-seq after multi-generational lineage tracking under controlled culture conditions. We use this platform to generate whole-transcriptome profiles of primary, activated murine CD8+ T-cell and lymphocytic leukemia cell line lineages. Here we report that both cell types have greater intra- than inter-lineage transcriptional similarity. For CD8+ T-cells, genes with functional annotation relating to lymphocyte differentiation and function--including Granzyme B--are enriched among the genes that demonstrate greater intra-lineage expression level similarity. Analysis of gene expression covariance with matched measurements of time since division reveals cell type-specific transcriptional signatures that correspond with cell cycle progression. We believe that the ability to directly measure the effects of lineage and cell cycle-dependent transcriptional profiles of single cells will be broadly useful to fields where heterogeneous populations of cells display distinct clonal trajectories, including immunology, cancer, and developmental biology.

Rapid and continuous magnetic separation in droplet microfluidic devices.

We present a droplet microfluidic method to extract molecules of interest from a droplet in a rapid and continuous fashion. We accomplish this by first marginalizing functionalized super-paramagnetic beads within the droplet using a magnetic field, and then splitting the droplet into one droplet containing the majority of magnetic beads and one droplet containing the minority fraction. We quantitatively analysed the factors which affect the efficiency of marginalization and droplet splitting to optimize the enrichment of magnetic beads. We first characterized the interplay between the droplet velocity and the strength of the magnetic field and its effect on marginalization. We found that marginalization is optimal at the midline of the magnet and that marginalization is a good predictor of bead enrichment through splitting at low to moderate droplet velocities. Finally, we focused our efforts on manipulating the splitting profile to improve the enrichment provided by asymmetric splitting. We designed asymmetric splitting forks that employ capillary effects to preferentially extract the bead-rich regions of the droplets. Our strategy represents a framework to optimize magnetic bead enrichment methods tailored to the requirements of specific droplet-based applications. We anticipate that our separation technology is well suited for applications in single-cell genomics and proteomics. In particular, our method could be used to separate mRNA bound to poly-dT functionalized magnetic microparticles from single cell lysates to prepare single-cell cDNA libraries.