Tailor: Targeting heavy tails in flow cytometry data with fast, interpretable mixture modeling.

Automated clustering workflows are increasingly used for the analysis of high parameter flow cytometry data. This trend calls for algorithms which are able to quickly process tens of millions of data points, to compare results across subjects or time points, and to provide easily actionable interpretations of the results. To this end, we created Tailor, a model-based clustering algorithm specialized for flow cytometry data. Our approach leverages a phenotype-aware binning scheme to provide a coarse model of the data, which is then refined using a multivariate Gaussian mixture model. We benchmark Tailor using a simulation study and two flow cytometry data sets, and show that the results are robust to moderate departures from normality and inter-sample variation. Moreover, Tailor provides automated, non-overlapping annotations of its clusters, which facilitates interpretation of results and downstream analysis. Tailor is released as an R package, and the source code is publicly available at www.github.com/matei-ionita/Tailor.

An integrated enrichment system to facilitate isolation and molecular characterization of single cancer cells from whole blood.

Circulating tumor cells (CTCs) carry valuable biological information. While enumeration of CTCs in peripheral blood is an FDA-approved prognostic indicator of survival in metastatic prostate and other cancers, analysis of CTC phenotypic and genomic markers is needed to identify cancer origin and elucidate pathways that can guide therapeutic selection for personalized medicine. Given the emergence of single-cell mRNA sequencing technologies, a method is needed to isolate CTCs with high sensitivity and specificity as well as compatibility with downstream genomic analysis. Flow cytometry is a powerful tool to analyze and sort single cells, but pre-enrichment is required prior to flow sorting for efficient isolation of CTCs due to the extreme low frequency of CTCs in blood (one in billions of blood cells). While current enrichment technologies often require many steps and result in poor recovery, we demonstrate a magnetic separator and acoustic microfluidic focusing chip integrated system that enriches rare cells in-line with FACS™ (fluorescent activated cell sorting) and single-cell sequencing. This system analyzes, isolates, and index sorts single cells directly into 96-well plates containing reagents for Molecular Indexing (MI) and transcriptional profiling of single cells. With an optimized workflow using the integrated enrichment-FACS system, we performed a proof-of-concept experiment with spiked prostate cancer cells in peripheral blood and achieved: (i) a rapid one-step process to isolate rare cancer cells from lysed whole blood; (ii) an average of 92% post-enrichment cancer cell recovery (R2 = 0.9998) as compared with 55% recovery for a traditional benchtop workflow; and (iii) detection of differentially expressed genes at a single cell level that are consistent with reported cell-type dependent expression signatures for prostate cancer cells. These model system results lay the groundwork for applying our approach to human blood samples from prostate and other cancer patients, and support the enrichment-FACS system as a flexible solution for isolation and characterization of CTCs for cancer diagnosis. © 2018 International Society for Advancement of Cytometry.

Potent obatoclax cytotoxicity and activation of triple death mode killing across infant acute lymphoblastic leukemia.

Survival in infants younger than 1 year who have acute lymphoblastic leukemia (ALL) is inferior whether MLL is rearranged (R) or germline (G). MLL translocations confer chemotherapy resistance, and infants experience excess complications. We characterized in vitro sensitivity to the pan-antiapoptotic BCL-2 family inhibitor obatoclax mesylate in diagnostic leukemia cells from 54 infants with ALL/bilineal acute leukemia because of the role of prosurvival BCL-2 proteins in resistance, their imbalanced expression in infant ALL, and evidence of obatoclax activity with a favorable toxicity profile in early adult leukemia trials. Overall, half maximal effective concentrations (EC50s) were lower than 176 nM (the maximal plasma concentration [Cmax] with recommended adult dose) in 76% of samples, whether in MLL-AF4, MLL-ENL, or other MLL-R or MLL-G subsets, and regardless of patients' poor prognostic features. However, MLL status and partner genes correlated with EC50. Combined approaches including flow cytometry, Western blot, obatoclax treatment with death pathway inhibition, microarray analyses, and/or electron microscopy indicated a unique killing mechanism involving apoptosis, necroptosis, and autophagy in MLL-AF4 ALL cell lines and primary MLL-R and MLL-G infant ALL cells. This in vitro obatoclax activity and its multiple killing mechanisms across molecular cytogenetic subsets provide a rationale to incorporate a similarly acting compound into combination strategies to combat infant ALL.

Immunologic signatures of response and resistance to nivolumab with ipilimumab in advanced metastatic cancer.

Identifying pan-tumor biomarkers that predict responses to immune checkpoint inhibitors (ICI) is critically needed. In the AMADEUS clinical trial (NCT03651271), patients with various advanced solid tumors were assessed for changes in intratumoral CD8 percentages and their response to ICI. Patients were grouped based on tumoral CD8 levels: those with CD8 <15% (CD8-low) received nivolumab (anti-PD-1) plus ipilimumab (anti-CTLA4) and those with CD8 ≥15% (CD8-high) received nivolumab monotherapy. 79 patients (72 CD8-low and 7 CD8-high) were treated. The disease control rate was 25.0% (18/72; 95% CI: 15.8-35.2) in CD8-low and 14.3% (1/7; 95% CI: 1.1-43.8) in CD8-high. Tumors from 35.9% (14/39; 95% CI: 21.8-51.4) of patients converted from CD8 <15% pretreatment to ≥15% after treatment. Multiomic analyses showed that CD8-low responders had an inflammatory tumor microenvironment pretreatment, enhanced by an influx of CD8 T cells, CD4 T cells, B cells, and macrophages upon treatment. These findings reveal crucial pan-cancer immunological features for ICI response in patients with metastatic disease.

Endothelial microparticles: sophisticated vesicles modulating vascular function.

Endothelial microparticles (EMPs) belong to a family of extracellular vesicles that are dynamic, mobile, biological effectors capable of mediating vascular physiology and function. The release of EMPs can impart autocrine and paracrine effects on target cells through surface interaction, cellular fusion, and, possibly, the delivery of intra-vesicular cargo. A greater understanding of the formation, composition, and function of EMPs will broaden our understanding of endothelial communication and may expose new pathways amenable for therapeutic manipulation.

Circulating endothelial progenitor cells are reduced in SLE in the absence of coronary artery calcification.

Circulating endothelial progenitor cells (EPCs) are reduced in patients with systemic lupus erythematosus (SLE). A reduced number of EPCs are associated with the presence of atherosclerosis in other populations. We sought to determine whether the reduction in EPC numbers in SLE is dependent on the presence of advanced coronary artery calcification (CAC). Patients with SLE had previous coronary calcium scores which placed them in either the >75th percentile or <25th percentile for their age. Seventeen patients with SLE and 13 healthy controls (HC) were included in the study. White blood cells were stained for EPC and progenitor cell markers including CD34, CD133, and VEGFR and analyzed by flow cytometry. SLE patients had repeated coronary imaging as well as carotid ultrasound. There was no difference in age between groups. SLE patients with advanced CAC were more likely to be hypertensive, to be smokers, and to have longer disease duration than SLE patients without CAC. SLE patients without evidence of CAC had a significantly lower number of EPCs (CD34+/CD133+/VEGFR+) compared to HC (median (IQR)) 0 (0, 6.7) vs. 10.2 (5.8, 12.3) (P = 0.02). Total numbers of PCs (CD133+/CD34+) were not significantly decreased in patients with SLE ((mean ± SEM) 1,007 ± 154 vs. 824 ± 170 (P = 0.20)). No significant difference was seen in EPC number between SLE patients without CAC and those with advanced CAC. Increased carotid intima-media thickness did not correlate with CAC or EPC number in SLE patients. Reduced numbers of EPCs in SLE patients may be observed compared to HC even in the absence of CAC. Differences in measured risk factor profiles and depletion of total circulating PCs do not fully explain this finding.

Metabolite patterns associated with individual response to supervised exercise therapy in patients with intermittent claudication.

Objective: Supervised exercise therapy (SET) is the first line treatment for intermittent claudication owing to peripheral arterial disease. Despite multiple randomized controlled trials proving the efficacy of SET, there are large differences in individual patient's responses. We used plasma metabolomics to identify potential metabolic influences on the individual response to SET. Methods: Primary metabolites, complex lipids, and lipid mediators were measured on plasma samples taken at before and after Gardner graded treadmill walking tests that were administered before and after 12 weeks of SET. We used an ensemble modeling approach to identify metabolites or changes in metabolites at specific time points that associated with interindividual variability in the functional response to SET. Specific time points analyzed included baseline metabolite levels before SET, dynamic metabolomics changes before SET, the difference in pre- and post-SET baseline metabolomics, and the difference (pre- and post-SET) of the dynamic (pre- and post-treadmill). Results: High levels of baseline anandamide levels pre- and post-SET were associated with a worse response to SET. Increased arachidonic acid (AA) and decreased levels of the AA precursor dihomo-γ-linolenic acid across SET were associated with a worse response to SET. Participants who were able to tolerate large increases in AA during acute exercise had longer, or better, walking times both before and after SET. Conclusions: We identified two pathways of relevance to individual response to SET that warrant further study: anandamide synthesis may activate endocannabinoid receptors, resulting in worse treadmill test performance. SET may train patients to withstand higher levels of AA, and inflammatory signaling, resulting in longer walking times. Clinical Relevance: This manuscript describes the use of metabolomic techniques to measure the interindividual effects of SET in patients with peripheral artery disease (PAD). We identified high levels of AEA are linked to CB1 signaling and activation of inflammatory pathways. This alters energy expenditure in myoblasts by decreasing glucose uptake and may induce an acquired skeletal muscle myopathy. SET may also help participants tolerate increased levels of AA and inflammation produced during exercise, resulting in longer walking times. This data will enhance understanding of the pathophysiology of PAD and the mechanism by which SET improves walking intolerance.

Circulating Donor Lung-specific Exosome Profiles Enable Noninvasive Monitoring of Acute Rejection in a Rodent Orthotopic Lung Transplantation Model.

BACKGROUND: There is a critical need for development of biomarkers to noninvasively monitor for lung transplant rejection. We investigated the potential of circulating donor lung-specific exosome profiles for time-sensitive diagnosis of acute rejection in a rat orthotopic lung transplant model. METHODS: Left lungs from Wistar transgenic rats expressing human CD63-GFP, an exosome marker, were transplanted into fully MHC-mismatched Lewis recipients or syngeneic controls. Recipient blood was collected between 4 h and 10 d after transplantation, and plasma was processed for exosome isolation by size exclusion column chromatography and ultracentrifugation. Circulating donor exosomes were profiled using antihuman CD63 antibody quantum dot on the nanoparticle detector and via GFP trigger on the nanoparticle flow cytometer. RESULTS: In syngeneic controls, steady-state levels of circulating donor exosomes were detected at all posttransplant time points. Allogeneic grafts lost perfusion by day 8, consistent with acute rejection. Levels of circulating donor exosomes peaked on day 1, decreased significantly by day 2, and then reached baseline levels by day 3. Notably, decrease in peripheral donor exosome levels occurred before grafts had histological evidence of acute rejection. CONCLUSIONS: Circulating donor lung-specific exosome profiles enable an early detection of acute rejection before histologic manifestation of injury to the pulmonary allograft. As acute rejection episodes are a major risk factor for the development of chronic lung allograft dysfunction, this biomarker may provide a novel noninvasive diagnostic platform that can translate into earlier therapeutic intervention for lung transplant patients.

Sotigalimab and/or nivolumab with chemotherapy in first-line metastatic pancreatic cancer: clinical and immunologic analyses from the randomized phase 2 PRINCE trial.

Chemotherapy combined with immunotherapy has improved the treatment of certain solid tumors, but effective regimens remain elusive for pancreatic ductal adenocarcinoma (PDAC). We conducted a randomized phase 2 trial evaluating the efficacy of nivolumab (nivo; anti-PD-1) and/or sotigalimab (sotiga; CD40 agonistic antibody) with gemcitabine/nab-paclitaxel (chemotherapy) in patients with first-line metastatic PDAC ( NCT03214250 ). In 105 patients analyzed for efficacy, the primary endpoint of 1-year overall survival (OS) was met for nivo/chemo (57.7%, P = 0.006 compared to historical 1-year OS of 35%, n = 34) but was not met for sotiga/chemo (48.1%, P = 0.062, n = 36) or sotiga/nivo/chemo (41.3%, P = 0.223, n = 35). Secondary endpoints were progression-free survival, objective response rate, disease control rate, duration of response and safety. Treatment-related adverse event rates were similar across arms. Multi-omic circulating and tumor biomarker analyses identified distinct immune signatures associated with survival for nivo/chemo and sotiga/chemo. Survival after nivo/chemo correlated with a less suppressive tumor microenvironment and higher numbers of activated, antigen-experienced circulating T cells at baseline. Survival after sotiga/chemo correlated with greater intratumoral CD4 T cell infiltration and circulating differentiated CD4 T cells and antigen-presenting cells. A patient subset benefitting from sotiga/nivo/chemo was not identified. Collectively, these analyses suggest potential treatment-specific correlates of efficacy and may enable biomarker-selected patient populations in subsequent PDAC chemoimmunotherapy trials.

Effect of darbepoetin alfa on endothelial progenitor cells and vascular reactivity in chronic kidney disease.

Endothelial progenitor cells (EPCs) are thought to be important for maintaining normal vascular function. We conducted a prospective study evaluating the effect of the erythropoiesis-stimulating agent darbepoetin alfa on EPCs and vascular function in patients with chronic kidney disease (CKD), with or without diabetes. Thirty subjects with CKD (20 subjects with type II diabetes mellitus and 10 without diabetes mellitus) received weekly subcutaneous administration of darbepoetin alfa for 4 weeks. EPCs were measured at baseline and 2 and 4 weeks after drug administration. Vascular function was measured with brachial ultrasound and cell activity was measured with a cell proliferation assay. Cells expressing CD133, CD34, CD146 and CD146/31 were significantly elevated (all p < 0.05), flow-mediated vasodilatation increased 2.1%, 95% CI: (0.4%, 3.8%) and colony-forming units increased twofold, 95% CI: (1.7, 2.3) after 4 weeks of treatment with darbepoetin alfa. Subjects with diabetes exhibited an increase in a subset of EPCs (CD133( +) and 34(+), p < 0.01 and p = 0.06, respectively), vasodilatation and proliferation. In conclusion, the administration of darbepoetin alfa for 4 weeks increased a subset of EPCs, improved endothelial function and increased cell proliferation, including those with diabetes, which is consistent with a favorable improvement in vascular health.

Reconcile the debate over protective effects of BCG vaccine against COVID-19.

While awaiting the COVID-19 vaccines, researchers have been actively exploring the effectiveness of existing vaccines against the new virus, among which the BCG vaccine (Bacillus Calmette-Guérin) receives the most attention. While many reports suggest a potential role for BCG immunization in ameliorating SARS-CoV-2 infection, these findings remain controversial. With country-level COVID-19 outbreak data from Johns Hopkins University Coronavirus Resource Center, and BCG program data from World Atlas of BCG Policies and Practices and WHO/UNICE, we estimated a dynamic model to investigate the effect of BCG vaccination across time during the pandemic. Our results reconcile these varying reports regarding protection by BCG against COVID-19 in a variety of clinical scenarios and model specifications. We observe a notable protective effect of the BCG vaccine during the early stage of the pandemic. However, we do not see any strong evidence for protection during the later stages. We also see that a higher proportion of vaccinated young population may confer some level of communal protection against the virus in the early pandemic period, even when the proportion of vaccination in the older population is low. Our results highlight that while BCG may offer some protection against COVID-19, we should be cautious in interpreting the estimated effectiveness as it may vary over time and depend on the age structure of the vaccinated population.

Caveolae are an essential component of the pathway for endothelial cell signaling associated with abrupt reduction of shear stress.

Abrupt cessation of flow representing the acute loss of shear stress (simulated ischemia) to flow-adapted pulmonary microvascular endothelial cells (PMVEC) leads to reactive oxygen species (ROS) generation that signals for EC proliferation. We evaluated the role of caveolin-1 on this cellular response with mouse PMVEC that were preconditioned for 72 h to laminar flow at 5 dyn/cm(2) followed by stop of flow ("ischemia"). Preconditioning resulted in a 2.7-fold increase in cellular expression of K(ATP) (K(IR) 6.2) channels but no change in expression level of caveolin-1, gp91(phox), or MAP kinases. The initial response to ischemia in wild type cells was cell membrane depolarization that was abolished by gene targeting of K(IR) 6.2. The subsequent response was increased ROS production associated with activation of NADPH oxidase (NOX2) and then phosphorylation of MAP kinases (Erk, JNK). After 24 h of ischemia in wild type cells, the cell proliferation index increased 2.5 fold and the % of cells in S+G(2)/M phases increased 6-fold. This signaling cascade (cell membrane depolarization, ROS production, MAP kinase activation and cell proliferation) was abrogated in caveolin-1 null PMVEC or by treatment of wild type cells with filipin. These studies indicate that caveolin-1 functions as a shear sensor in flow-adapted EC resulting in ROS-mediated cell signaling and endothelial cell proliferation following the abrupt reduction in flow.

Diabetes reduces bone marrow and circulating porcine endothelial progenitor cells, an effect ameliorated by atorvastatin and independent of cholesterol.

Bone marrow derived endothelial progenitor cells (EPCs) are early precursors of mature endothelial cells which replenish aging and damaged endothelial cells. The authors studied a diabetic swine model to determine if induction of DM adversely affects either bone marrow or circulating EPCs and whether a HMG-CoA reductase inhibitor (statin) improves development and recruitment of EPCs in the absence of cholesterol lowering. Streptozotocin was administered to Yorkshire pigs to induce DM. One month after induction, diabetic pigs were treated with atorvastatin (statin, n = 10), ezetimibe (n = 10) or untreated (n = 10) and evaluated for number of bone marrow and circulating EPCs and femoral artery endothelial function. There was no effect of either medication on cholesterol level. One month after induction of DM prior to administration of drugs, the number of bone marrow and circulating EPCs significantly decreased (P < 0.0001) compared to baseline. Three months after DM induction, the mean proportion of circulating EPCs significantly increased in the atorvastatin group, but not in the control or ezetimibe groups. The control group showed progressive reduction in percentage of flow mediated vasodilatation (no dilatation at 3 months) whereas the atorvastatin group and ezetimibe exhibited vasodilatation, 6% and 4% respectively. DM results in significant impairment of bone marrow and circulating EPCs as well as endothelial function. The effect is ameliorated, in part, by atorvastatin independent of its cholesterol lowering effect. These data suggest a model wherein accelerated atherosclerosis seen with DM may, in part, result from reduction in EPCs which may be ameliorated by treatment with a statin.

Hyperbaric oxygen stimulates vasculogenic stem cell growth and differentiation in vivo.

We hypothesized that oxidative stress from hyperbaric oxygen (HBO(2), 2.8 ATA for 90 min daily) exerts a trophic effect on vasculogenic stem cells. In a mouse model, circulating stem/progenitor cell (SPC) recruitment and differentiation in subcutaneous Matrigel were stimulated by HBO(2) and by a physiological oxidative stressor, lactate. In combination, HBO(2) and lactate had additive effects. Vascular channels lined by CD34(+) SPCs were identified. HBO(2) and lactate accelerated channel development, cell differentiation based on surface marker expression, and cell cycle entry. CD34(+) SPCs exhibited increases in thioredoxin-1 (Trx1), Trx reductase, hypoxia-inducible factors (HIF)-1, -2, and -3, phosphorylated mitogen-activated protein kinases, vascular endothelial growth factor, and stromal cell-derived factor-1. Cell recruitment to Matrigel and protein synthesis responses were abrogated by N-acetyl cysteine, dithioerythritol, oxamate, apocynin, U-0126, neutralizing anti-vascular endothelial growth factor, or anti-stromal cell-derived factor-1 antibodies, and small inhibitory RNA to Trx reductase, lactate dehydrogenase, gp91(phox), HIF-1 or -2, and in mice conditionally null for HIF-1 in myeloid cells. By causing an oxidative stress, HBO(2) activates a physiological redox-active autocrine loop in SPCs that stimulates vasculogenesis. Thioredoxin system activation leads to elevations in HIF-1 and -2, followed by synthesis of HIF-dependent growth factors. HIF-3 has a negative impact on SPCs.

Targeting EIF4E signaling with ribavirin in infant acute lymphoblastic leukemia.

The poor outcomes in infant acute lymphoblastic leukemia (ALL) necessitate new treatments. Here we discover that EIF4E protein is elevated in most cases of infant ALL and test EIF4E targeting by the repurposed antiviral agent ribavirin, which has anticancer properties through EIF4E inhibition, as a potential treatment. We find that ribavirin treatment of actively dividing infant ALL cells on bone marrow stromal cells (BMSCs) at clinically achievable concentrations causes robust proliferation inhibition in proportion with EIF4E expression. Further, we find that ribavirin treatment of KMT2A-rearranged (KMT2A-R) infant ALL cells and the KMT2A-AFF1 cell line RS4:11 inhibits EIF4E, leading to decreases in oncogenic EIF4E-regulated cell growth and survival proteins. In ribavirin-sensitive KMT2A-R infant ALL cells and RS4:11 cells, EIF4E-regulated proteins with reduced levels of expression following ribavirin treatment include MYC, MCL1, NBN, BCL2 and BIRC5. Ribavirin-treated RS4:11 cells exhibit impaired EIF4E-dependent nuclear to cytoplasmic export and/or translation of the corresponding mRNAs, as well as reduced phosphorylation of the p-AKT1, p-EIF4EBP1, p-RPS6 and p-EIF4E signaling proteins. This leads to an S-phase cell cycle arrest in RS4:11 cells corresponding to the decreased proliferation. Ribavirin causes nuclear EIF4E to re-localize to the cytoplasm in KMT2A-AFF1 infant ALL and RS4:11 cells, providing further evidence for EIF4E inhibition. Ribavirin slows increases in peripheral blasts in KMT2A-R infant ALL xenograft-bearing mice. Ribavirin cooperates with chemotherapy, particularly L-asparaginase, in reducing live KMT2A-AFF1 infant ALL cells in BMSC co-cultures. This work establishes that EIF4E is broadly elevated across infant ALL and that clinically relevant ribavirin exposures have preclinical activity and effectively inhibit EIF4E in KMT2A-R cases, suggesting promise in EIF4E targeting using ribavirin as a means of treatment.

Abundant anti-apoptotic BCL-2 is a molecular target in leukaemias with t(4;11) translocation.

Chemotherapy resistance from imbalanced apoptosis regulation may contribute to poor outcome in leukaemias with t(4;11). Anti-apoptotic BCL-2 expression and target modulation were characterized in cell lines with t(4;11) and BCL-2 expression was examined in MLL and non-MLL infant/paediatric leukaemia cases by Western blot analysis and/or real-time polymerase chain reaction. Cytotoxicity of Genasensetrade mark (Oblimersen Sodium, G3139) alone or combined with cytotoxic drugs was assessed by MTT [(3-4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide] assays of the cell lines, applying pharmacostatistical response surface modelling of drug interactions. Apoptosis and cell cycle were evaluated by flow cytometry in RS4:11 cells. Primary leukaemias and cell lines with t(4;11) expressed abundant BCL2 mRNA and protein. Variable, sometimes substantial BCL2 mRNA was detected in other leukaemia subtypes. G3139 reduced BCL2 mRNA and protein in RS4:11 cells. The most sensitive cell line to single-agent G3139 was RS4:11. Low G3139 concentrations sensitized RS4:11 and MV4-11 cells to select anti-leukaemia cytotoxic drugs. In RS4:11 cells, combining G3139 with doxorubicin (ADR) increased active caspase 3 and TUNEL staining compared to ADR alone, indicating greater apoptosis, and G3139 increased S-phase progression. The abundant BCL-2 affords a molecular target in leukaemias with t(4;11). G3139 exhibits preclinical activity and synergy with select cytotoxic agents in RS4:11 and MV4-11 cells, and these effects occur through apoptosis.

Cytometric fingerprinting: quantitative characterization of multivariate distributions.

Recent technological advances in flow cytometry instrumentation provide the basis for high-dimensionality and high-throughput biological experimentation in a heterogeneous cellular context. Concomitant advances in scalable computational algorithms are necessary to better utilize the information that is contained in these high-complexity experiments. The development of such tools has the potential to expand the utility of flow cytometric analysis from a predominantly hypothesis-driven mode to one of discovery, or hypothesis-generating research. A new method of analysis of flow cytometric data called Cytometric Fingerprinting (CF) has been developed. CF captures the set of multivariate probability distribution functions corresponding to list-mode data and then "flattens" them into a computationally efficient fingerprint representation that facilitates quantitative comparisons of samples. An experimental and synthetic data were generated to act as reference sets for evaluating CF. Without the introduction of prior knowledge, CF was able to "discover" the location and concentration of spiked cells in ungated analyses over a concentration range covering four orders of magnitude, to a lower limit on the order of 10 spiked events in a background of 100,000 events. We describe a new method for quantitative analysis of list-mode cytometric data. CF includes a novel algorithm for space subdivision that improves estimation of the probability density function by dividing space into nonrectangular polytopes. Additionally it renders a multidimensional distribution in the form of a one-dimensional multiresolution hierarchical fingerprint that creates a computationally efficient representation of high dimensionality distribution functions. CF supports both the generation and testing of hypotheses, eliminates sources of operator bias, and provides an increased level of automation of data analysis.

An integrated flow cytometry-based platform for isolation and molecular characterization of circulating tumor single cells and clusters.

Comprehensive molecular analysis of rare circulating tumor cells (CTCs) and cell clusters is often hampered by low throughput and purity, as well as cell loss. To address this, we developed a fully integrated platform for flow cytometry-based isolation of CTCs and clusters from blood that can be combined with whole transcriptome analysis or targeted RNA transcript quantification. Downstream molecular signature can be linked to cell phenotype through index sorting. This newly developed platform utilizes in-line magnetic particle-based leukocyte depletion, and acoustic cell focusing and washing to achieve >98% reduction of blood cells and non-cellular debris, along with >1.5 log-fold enrichment of spiked tumor cells. We could also detect 1 spiked-in tumor cell in 1 million WBCs in 4/7 replicates. Importantly, the use of a large 200µm nozzle and low sheath pressure (3.5 psi) minimized shear forces, thereby maintaining cell viability and integrity while allowing for simultaneous recovery of single cells and clusters from blood. As proof of principle, we isolated and transcriptionally characterized 63 single CTCs from a genetically engineered pancreatic cancer mouse model (n = 12 mice) and, using index sorting, were able to identify distinct epithelial and mesenchymal sub-populations based on linked single cell protein and gene expression.