bilayer-microfluidic-device-for-combinatorial-plug-production

Bilayer Microfluidic Device for Combinatorial Plug Production

Droplet microfluidics is a versatile tool that allows the execution of a large number of reactions in chemically distinct nanoliter compartments. Such ...

Eindhoven University of Technology

Screening drugs on patient biopsies from solid tumours has immense potential, but is challenging due to the small amount of available material. To address this, we present here a plug-based microfluidics platform for functional screening of drug combinations. Integrated Braille valves allow changing the plug composition on demand and enable collecting >1200 data points (56 different conditions with at least 20 replicates each) per biopsy. After deriving and validating efficient and specific drug combinations for two genetically different pancreatic cancer cell lines and xenograft mouse models, we additionally screen live cells from human solid tumours with no need for ex vivo culturing steps, and obtain highly specific sensitivity profiles. The entire workflow can be completed within 48 h at assay costs of less than US$ 150 per patient. We believe this can pave the way for rapid determination of optimal personalized cancer therapies.

A microfluidics platform for combinatorial drug screening on cancer biopsies.

The tumor microenvironment (TME) is a multifaceted ecosystem characterized by profound cellular heterogeneity, dynamicity, and complex intercellular cross-talk. The striking responses obtained with immune checkpoint blockers, i.e., antibodies targeting immune-cell regulators to boost antitumor immunity, have demonstrated the enormous potential of anticancer treatments that target TME components other than tumor cells. However, as checkpoint blockade is currently beneficial only to a limited fraction of patients, there is an urgent need to understand the mechanisms orchestrating the immune response in the TME to guide the rational design of more effective anticancer therapies. In this Mini Review, we give an overview of the methodologies that allow studying the heterogeneity of the TME from multi-omics data generated from bulk samples, single cells, or images of tumor-tissue slides. These include approaches for the characterization of the different cell phenotypes and for the reconstruction of their spatial organization and inter-cellular cross-talk. We discuss how this broader vision of the cellular heterogeneity and plasticity of tumors, which is emerging thanks to these methodologies, offers the opportunity to rationally design precision immuno-oncology treatments. These developments are fundamental to overcome the current limitations of targeted agents and checkpoint blockers and to bring long-term clinical benefits to a larger fraction of cancer patients.

Multi-Omics Profiling of the Tumor Microenvironment: Paving the Way to Precision Immuno-Oncology.

Mechanistic modeling of signaling pathways mediating patient-specific response to therapy can help to unveil resistance mechanisms and improve therapeutic strategies. Yet, creating such models for patients, in particular for solid malignancies, is challenging. A major hurdle to build these models is the limited material available that precludes the generation of large-scale perturbation data. Here, we present an approach that couples ex vivo high-throughput screenings of cancer biopsies using microfluidics with logic-based modeling to generate patient-specific dynamic models of extrinsic and intrinsic apoptosis signaling pathways. We used the resulting models to investigate heterogeneity in pancreatic cancer patients, showing dissimilarities especially in the PI3K-Akt pathway. Variation in model parameters reflected well the different tumor stages. Finally, we used our dynamic models to efficaciously predict new personalized combinatorial treatments. Our results suggest that our combination of microfluidic experiments and mathematical model can be a novel tool toward cancer precision medicine.

Patient-specific logic models of signaling pathways from screenings on cancer biopsies to prioritize personalized combination therapies.

Immunotherapy with checkpoint blockers (ICBs), aimed at unleashing the immune response toward tumor cells, has shown a great improvement in overall patient survival compared to standard therapy, but only in a subset of patients. While a number of recent studies have significantly improved our understanding of mechanisms playing an important role in the tumor microenvironment (TME), we still have an incomplete view of how the TME works as a whole. This hampers our ability to effectively predict the large heterogeneity of patients' response to ICBs. Systems approaches could overcome this limitation by adopting a holistic perspective to analyze the complexity of tumors. In this Mini Review, we focus on how an integrative view of the increasingly available multi-omics experimental data and computational approaches enables the definition of new systems-based predictive biomarkers. In particular, we will focus on three facets of the TME toward the definition of new systems biomarkers. First, we will review how different types of immune cells influence the efficacy of ICBs, not only in terms of their quantification, but also considering their localization and functional state. Second, we will focus on how different cells in the TME interact, analyzing how inter- and intra-cellular networks play an important role in shaping the immune response and are responsible for resistance to immunotherapy. Finally, we will describe the potential of looking at these networks as dynamic systems and how mathematical models can be used to study the rewiring of the complex interactions taking place in the TME.

Toward Systems Biomarkers of Response to Immune Checkpoint Blockers.

Neuron-specific enolase (NSE) is a well-known biomarker for the diagnosis, prognosis and treatment monitoring of small-cell lung cancer (SCLC). Nevertheless, its clinical applicability is limited since serum NSE levels are influenced by hemolysis, leading to falsely elevated results. Therefore, this study aimed to develop a hemolysis correction equation and evaluate its role in SCLC diagnostics. Two serum pools were spiked with increasing amounts of hemolysate obtained from multiple individuals. A hemolysis correction equation was obtained by analyzing the relationship between the measured NSE concentration and the degree of hemolysis. The equation was validated using intentionally hemolyzed serum samples, which showed that the correction was accurate for samples with an H-index up to 30 μmol/L. Correction of the measured NSE concentration in patients suspected of lung cancer caused an increase in AUC and a significantly lower cut-off value for SCLC detection when compared to uncorrected results. Therefore, a hemolysis correction equation should be used to correct falsely elevated NSE concentrations. Results of samples with an H-index above 30 μmol/L should not be reported to clinicians. Application of the equation illustrates the importance of hemolysis correction in SCLC diagnostics and questions the correctness of the currently used diagnostic cut-off value.

Correction of the NSE concentration in hemolyzed serum samples improves its diagnostic accuracy in small-cell lung cancer.

Cancer cells can leverage several cell-intrinsic and -extrinsic mechanisms to escape immune system recognition. The inherent complexity of the tumor microenvironment, with its multicellular and dynamic nature, poses great challenges for the extraction of biomarkers of immune response and immunotherapy efficacy. Here, we use RNA-sequencing (RNA-seq) data combined with different sources of prior knowledge to derive system-based signatures of the tumor microenvironment, quantifying immune-cell composition and intra- and intercellular communications. We applied multi-task learning to these signatures to predict different hallmarks of immune responses and derive cancer-type-specific models based on interpretable systems biomarkers. By applying our models to independent RNA-seq data from cancer patients treated with PD-1/PD-L1 inhibitors, we demonstrated that our method to Estimate Systems Immune Response (EaSIeR) accurately predicts therapeutic outcome. We anticipate that EaSIeR will be a valuable tool to provide a holistic description of immune responses in complex and dynamic systems such as tumors using available RNA-seq data.

Interpretable systems biomarkers predict response to immune-checkpoint inhibitors.

Engineering and computational advances have opened many new avenues in cancer research, particularly when being exploited in interdisciplinary approaches. For example, the combination of microfluidics, novel sequencing technologies, and computational analyses has been crucial to enable single-cell assays, giving a detailed picture of tumor heterogeneity for the very first time. In a similar way, these 'tech' disciplines have been elementary for generating large data sets in multidimensional cancer 'omics' approaches, cell-cell interaction screens, 3D tumor models, and tissue level analyses. In this review we summarize the most important technology and computational developments that have been or will be instrumental for transitioning classical cancer research to a large data-driven, high-throughput, high-content discipline across all biological scales.

Technological and computational advances driving high-throughput oncology.

Identification of actionable mutations in advanced stage non-squamous non-small-cell lung cancer (NSCLC) patients is recommended by guidelines as it enables treatment with targeted therapies. In current practice, mutations are identified by next-generation sequencing of tumor DNA (tDNA-NGS), which requires tissue biopsies of sufficient quality. Alternatively, circulating tumor DNA (ctDNA) could be used for mutation analysis. This prospective, multicenter study establishes the diagnostic value of ctDNA analysis by droplet digital PCR (ctDNA-ddPCR) in patients with primary lung cancer. CtDNA from 458 primary lung cancer patients was analyzed using a panel of multiplex ddPCRs for EGFR (Ex19Del, G719S, L858R, L861Q and S768I), KRAS G12/G13 and BRAF V600 mutations. For 142 of 175 advanced stage non-squamous NSCLC patients tDNA-NGS results were available to compare to ctDNA-ddPCR. tDNA-NGS identified 98 mutations, of which ctDNA-ddPCR found 53 mutations (54%), including 32 of 45 (71%) targetable driver mutations. In 2 of these 142 patients, a mutation was found by ctDNA-ddPCR only. In 33 advanced stage patients lacking tDNA-NGS results, ctDNA-ddPCR detected 15 additional mutations, of which 7 targetable. Overall, ctDNA-ddPCR detected 70 mutations and tDNA-NGS 98 mutations in 175 advanced NSCLC patients. Using an up-front ctDNA-ddPCR strategy, followed by tDNA-NGS only if ctDNA-ddPCR analysis is negative, increases the number of mutations found from 98 to 115 (17%). At the same time, up-front ctDNA-ddPCR reduces tDNA-NGS analyses by 40%, decreasing the need to perform (additional) biopsies.

Up-front mutation detection in circulating tumor DNA by droplet digital PCR has added diagnostic value in lung cancer.

Pathologic subtyping of tissue biopsies is the gold standard for the diagnosis of lung cancer (LC), which could be complicated in cases of e.g. inconclusive tissue biopsies or unreachable tumors. The diagnosis of LC could be supported in a minimally invasive manner using protein tumor markers (TMs) and circulating tumor DNA (ctDNA) measured in liquid biopsies (LBx). This study evaluates the performance of LBx-based decision-support algorithms for the diagnosis of LC and subtyping into small- and non-small-cell lung cancer (SCLC and NSCLC) aiming to directly impact clinical practice.

Liquid biopsy-based decision support algorithms for diagnosis and subtyping of lung cancer.

Prostate cancer (PCa) is a highly prevalent cancer type with a heterogeneous prognosis. An accurate assessment of tumor aggressiveness can pave the way for tailored treatment strategies, potentially leading to better outcomes. While tumor aggressiveness is typically assessed based on invasive methods (e.g., biopsy), radiogenomics, combining diagnostic imaging with genomic information can help uncover aggressive (imaging) phenotypes, which in turn can provide non-invasive advice on individualized treatment regimens. In this study, we carried out a parallel analysis on both imaging and transcriptomics data in order to identify features associated with clinically significant PCa (defined as an ISUP grade ≥ 3), subsequently evaluating the correlation between them. Textural imaging features were extracted from multi-parametric MRI sequences (T2W, DWI, and DCE) and combined with DCE-derived parametric pharmacokinetic maps obtained using magnetic resonance dispersion imaging (MRDI). A transcriptomic analysis was performed to derive functional features on transcription factors (TFs), and pathway activity from RNA sequencing data, here referred to as transcriptomic features. For both the imaging and transcriptomic features, different machine learning models were separately trained and optimized to classify tumors in either clinically insignificant or significant PCa. These models were validated in an independent cohort and model performance was used to isolate a subset of relevant imaging and transcriptomic features to be further investigated. A final set of 31 imaging features was correlated to 33 transcriptomic features obtained on the same tumors. Five significant correlations ( < 0.05) were found, of which, three had moderate strength (|r| ≥ 0.5). The strongest significant correlations were seen between a perfusion-based imaging feature-MRDI  median-and the activities of the TFs STAT6 (-0.64) and TFAP2A (-0.50). A higher-order T2W textural feature was also significantly correlated to the activity of the TF STAT6 (-0.58). STAT6 plays an important role in controlling cell proliferation and migration. Loss of the AP2alpha protein expression, quantified by TFAP2A, has been strongly associated with aggressiveness and progression in PCa. According to our findings, a combination of texture features extracted from T2W and DCE, as well as perfusion-based pharmacokinetic features, can be considered for the prediction of clinically significant PCa, with the pharmacokinetic MRDI  feature being the most correlated with the underlying transcriptomic information. These results highlight a link between quantitative imaging features and the underlying transcriptomic landscape of prostate tumors.

Radiogenomics Analysis Linking Multiparametric MRI and Transcriptomics in Prostate Cancer.

Over 10% of men will be diagnosed with prostate cancer during their lifetime. Arising from luminal cells of the prostatic acinus, prostate cancer is influenced by multiple cells in its microenvironment. To expand our knowledge and explore means to prevent and treat the disease, it is important to understand what drives the onset and early stages of prostate cancer. In this study, we developed an agent-based model of a prostatic acinus including its microenvironment, to allow for  studying of prostate cancer development. The model was based on prior reports and in-house data of tumor cells cocultured with cancer-associated fibroblasts (CAF) and protumor and/or antitumor macrophages. Growth patterns depicted by the model were pathologically validated on hematoxylin and eosin slide images of human prostate cancer specimens. We identified that stochasticity of interactions between macrophages and tumor cells at early stages strongly affect tumor development. In addition, we discovered that more systematic deviations in tumor development result from a combinatorial effect of the probability of acquiring mutations and the tumor-promoting abilities of CAFs and macrophages.  modeled tumors were then compared with 494 patients with cancer with matching characteristics, showing strong association between predicted tumor load and patients' clinical outcome. Our findings suggest that the likelihood of tumor formation depends on a combination of stochastic events and systematic characteristics. While stochasticity cannot be controlled, information on systematic effects may aid the development of prevention strategies tailored to the molecular characteristics of an individual patient.

Federica Eduati

Department of Biomedical Engineering,
Institute for Complex Molecular Systems,
Department of Biomedical Engineering, Institute for Complex Molecular Systems

Eindhoven University of Technology

Bilayer Microfluidic Device for Combinatorial Plug Production

Federica Eduati

Department of Biomedical Engineering,Institute for Complex Molecular Systems,Department of Biomedical Engineering, Institute for Complex Molecular Systems

Eindhoven University of Technology

Bilayer Microfluidic Device for Combinatorial Plug Production

Department of Biomedical Engineering,
Institute for Complex Molecular Systems,
Department of Biomedical Engineering, Institute for Complex Molecular Systems