Faculty Bibliography

The effectiveness of most cancer targeted therapies is short-lived. Tumors often develop resistance that might be overcome with drug combinations. However, the number of possible combinations is vast, necessitating data-driven approaches to find optimal patient-specific treatments. Here we report AstraZeneca's large drug combination dataset, consisting of 11,576 experiments from 910 combinations across 85 molecularly characterized cancer cell lines, and results of a DREAM Challenge to evaluate computational strategies for predicting synergistic drug pairs and biomarkers. 160 teams participated to provide a comprehensive methodological development and benchmarking. Winning methods incorporate prior knowledge of drug-target interactions. Synergy is predicted with an accuracy matching biological replicates for >60% of combinations. However, 20% of drug combinations are poorly predicted by all methods. Genomic rationale for synergy predictions are identified, including ADAM17 inhibitor antagonism when combined with PIK3CB/D inhibition contrasting to synergy when combined with other PI3K-pathway inhibitors in PIK3CA mutant cells.

We studied the trajectories of polymers being advected while diffusing in a pressure driven flow along a periodic pillar nanostructure known as nanoscale deterministic lateral displacement (nanoDLD) array. We found that polymers follow different trajectories depending on their length, flow velocity and pillar array geometry, demonstrating that nanoDLD devices can be used as a continuous polymer fractionation tool. As a model system, we used double-stranded DNA (dsDNA) with various contour lengths and demonstrated that dsDNA in the range of 100-10 000 base pairs (bp) can be separated with a size-selective resolution of 200 bp. In contrast to spherical colloids, a polymer elongates by shear flow and the angle of polymer trajectories with respect to the mean flow direction decreases as the mean flow velocity increases. We developed a phenomenological model that explains the qualitative dependence of the polymer trajectories on the gap size and on the flow velocity. Using this model, we found the optimal separation conditions for dsDNA of different sizes and demonstrated the separation and extraction of dsDNA fragments with over 75% recovery and 3-fold concentration. Importantly, this velocity dependence provides a means of fine-tuning the separation efficiency and resolution, independent of the nanoDLD pillar geometry.

To develop a map of cell-cell communication mediated by extracellular RNA (exRNA), the NIH Extracellular RNA Communication Consortium created the exRNA Atlas resource (https://exrna-atlas.org). The Atlas version 4P1 hosts 5,309 exRNA-seq and exRNA qPCR profiles from 19 studies and a suite of analysis and visualization tools. To analyze variation between profiles, we apply computational deconvolution. The analysis leads to a model with six exRNA cargo types (CT1, CT2, CT3A, CT3B, CT3C, CT4), each detectable in multiple biofluids (serum, plasma, CSF, saliva, urine). Five of the cargo types associate with known vesicular and non-vesicular (lipoprotein and ribonucleoprotein) exRNA carriers. To validate utility of this model, we re-analyze an exercise response study by deconvolution to identify physiologically relevant response pathways that were not detected previously. To enable wide application of this model, as part of the exRNA Atlas resource, we provide tools for deconvolution and analysis of user-provided case-control studies.

Individual cells in clonal populations often respond differently to environmental changes; for binary phenotypes, such as cell death, this can be measured as a fractional response. These types of responses have been attributed to cell-intrinsic stochastic processes and variable abundances of biochemical constituents, such as proteins, but the influence of organelles is still under investigation. We use the response to TNF-related apoptosis inducing ligand (TRAIL) and a new statistical framework for determining parameter influence on cell-to-cell variability through the inference of variance explained, DEPICTIVE, to demonstrate that variable mitochondria abundance correlates with cell survival and determines the fractional cell death response. By quantitative data analysis and modeling we attribute this effect to variable effective concentrations at the mitochondria surface of the pro-apoptotic proteins Bax/Bak. Further, our study suggests that inhibitors of anti-apoptotic Bcl-2 family proteins, used in cancer treatment, may increase the diversity of cellular responses, enhancing resistance to treatment.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease where substantial heterogeneity in clinical presentation urgently requires a better stratification of patients for the development of drug trials and clinical care. In this study we explored stratification through a crowdsourcing approach, the DREAM Prize4Life ALS Stratification Challenge. Using data from >10,000 patients from ALS clinical trials and 1479 patients from community-based patient registers, more than 30 teams developed new approaches for machine learning and clustering, outperforming the best current predictions of disease outcome. We propose a new method to integrate and analyze patient clusters across methods, showing a clear pattern of consistent and clinically relevant sub-groups of patients that also enabled the reliable classification of new patients. Our analyses reveal novel insights in ALS and describe for the first time the potential of a crowdsourcing to uncover hidden patient sub-populations, and to accelerate disease understanding and therapeutic development.

Extracellular vesicles (EVs) offer many opportunities in early-stage disease diagnosis, treatment monitoring, and precision therapy owing to their high abundance in bodily fluids, accessibility from liquid biopsy, and presence of nucleic acid and protein cargo from their cell of origin. Despite their growing promise, isolation of EVs for analysis remains a labor-intensive and time-consuming challenge given their nanoscale dimensions (30-200 nm) and low buoyant density. Here, we report a simple, size-based EV separation technology that integrates 1024 nanoscale deterministic lateral displacement (nanoDLD) arrays on a single chip capable of parallel processing sample fluids at rates of up to 900 μL h-1. Benchmarking the nanoDLD chip against commonly used EV isolation technologies, including ultracentrifugation (UC), UC plus density gradient, qEV size-exclusion chromatography (Izon Science), and the exoEasy Maxi Kit (QIAGEN), we demonstrate a superior yield of ∼50% for both serum and urine samples, representing the ability to use smaller input volumes to achieve the same number of isolated EVs, and a concentration factor enhancement of up to ∼3× for both sample types, adjustable to ∼60× for urine through judicious design. Further, RNA sequencing was carried out on nanoDLD- and UC-isolated EVs from prostate cancer (PCa) patient serum samples, resulting in a higher gene expression correlation between replicates for nanoDLD-isolated EVs with enriched miRNA, decreased rRNA, and the ability to detect previously reported RNA indicators of aggressive PCa. Taken together, these results suggest nanoDLD as a promising alternative technology for fast, reproducible, and automatable EV-isolation.

BACKGROUND:The phenotypes of cancer cells are driven in part by somatic structural variants. Structural variants can initiate tumors, enhance their aggressiveness, and provide unique therapeutic opportunities. Whole-genome sequencing of tumors can allow exhaustive identification of the specific structural variants present in an individual cancer, facilitating both clinical diagnostics and the discovery of novel mutagenic mechanisms. A plethora of somatic structural variant detection algorithms have been created to enable these discoveries; however, there are no systematic benchmarks of them. Rigorous performance evaluation of somatic structural variant detection methods has been challenged by the lack of gold standards, extensive resource requirements, and difficulties arising from the need to share personal genomic information. RESULTS:To facilitate structural variant detection algorithm evaluations, we create a robust simulation framework for somatic structural variants by extending the BAMSurgeon algorithm. We then organize and enable a crowdsourced benchmarking within the ICGC-TCGA DREAM Somatic Mutation Calling Challenge (SMC-DNA). We report here the results of structural variant benchmarking on three different tumors, comprising 204 submissions from 15 teams. In addition to ranking methods, we identify characteristic error profiles of individual algorithms and general trends across them. Surprisingly, we find that ensembles of analysis pipelines do not always outperform the best individual method, indicating a need for new ways to aggregate somatic structural variant detection approaches. CONCLUSIONS:The synthetic tumors and somatic structural variant detection leaderboards remain available as a community benchmarking resource, and BAMSurgeon is available at https://github.com/adamewing/bamsurgeon .

We characterize the transcriptional splicing landscape of a prostate cancer cell line treated with a previously identified synergistic drug combination. We use a combination of third generation long-read RNA sequencing technology and short-read RNAseq to create a high-fidelity map of expressed isoforms and fusions to quantify splicing events triggered by treatment. We find strong evidence for drug-induced, coherent splicing changes which disrupt the function of oncogenic proteins, and detect novel transcripts arising from previously unreported fusion events.

Whole cell responses arise from coordinated interactions between diverse human gene products functioning within various pathways underlying sub-cellular processes (SCP). Lower level SCPs interact to form higher level SCPs, often in a context specific manner to give rise to whole cell function. We sought to determine if capturing such relationships enables us to describe the emergence of whole cell functions from interacting SCPs. We developed the Molecular Biology of the Cell Ontology based on standard cell biology and biochemistry textbooks and review articles. Currently, our ontology contains 5,384 genes, 753 SCPs and 19,180 expertly curated gene-SCP associations. Our algorithm to populate the SCPs with genes enables extension of the ontology on demand and the adaption of the ontology to the continuously growing cell biological knowledge. Since whole cell responses most often arise from the coordinated activity of multiple SCPs, we developed a dynamic enrichment algorithm that flexibly predicts SCP-SCP relationships beyond the current taxonomy. This algorithm enables us to identify interactions between SCPs as a basis for higher order function in a context dependent manner, allowing us to provide a detailed description of how SCPs together can give rise to whole cell functions. We conclude that this ontology can, from omics data sets, enable the development of detailed SCP networks for predictive modeling of emergent whole cell functions.