Biological Insights from Geometry and Structure of Single-Cell Data

Understanding the behavior of a cell requires that its molecular constituents, such as mRNA or protein levels, be profiled quantitatively. Typically, these measurements are performed in bulk and represent values aggregated from thousands of cells. Insights from such data can be very useful, but the loss of single-cell resolution can prove misleading for heterogeneous tissues and in diseases like cancer. Recently, technological advances have allowed us to profile multiple cellular parameters simultaneously at single-cell resolution, for thousands to millions of cells. While this provides an unprecedented opportunity to learn new biology, analyzing such massive and high-dimensional data requires efficient and accurate computational tools to extract the underlying biological phenomena. Such methods must take into account biological properties such as non-linear dependencies between measured parameters. In this dissertation, I contribute to the development of tools from harmonic analysis and computational geometry to study the shape and geometry of single-cell data collected using mass cytometry and single-cell RNA sequencing (scRNA-seq). In particular, I focus on diffusion maps, which can learn the underlying structure of the data by modeling cells as lying on a low-dimensional phenotype manifold embedded in high dimensions. Diffusion maps allow non-linear transformation of the data into a low-dimensional Euclidean space, in which pairwise distances robustly represent distances in the high-dimensional space. In addition to the underlying geometry, this work also attempts to study the shape of the data using archetype analysis. Archetype analysis characterizes extreme states in the data and complements traditional approaches such as clustering. It facilitates analysis at the boundary of the data enabling potentially novel insights about the system. I use these tools to study how the negative costimulatory molecules Ctla4 and Pdcd1 affect T-cell differentiation. Negative costimulatory molecules play a vital role in attenuating T-cell activation, in order to maintain activity within a desired physiological range and prevent autoimmunity. However, their potential role in T cell differentiation remains unknown. In this work, I analyze mass cytometry data profiling T cells in control and Ctla4- or Pdcd1-deficient mice and analyze differences using the tools above. I find that genetic loss of Ctla4 constrains CD4+ T-cell differentiation states, whereas loss of Pdcd1 subtly constrains CD8+ T-cell differentiation states. I propose that negative costimulatory molecules place limits on maximal protein expression levels to restrain differentiation states. I use similar approaches to study breast cancer cells, which are profiled using scRNA-seq as they undergo the pathological epithelial-to-mesenchymal transition (EMT). For this work, I introduce Markov Affinity based Graph Imputation of Cells (MAGIC), a novel algorithm designed in our lab to denoise and impute sparse single-cell data. The mRNA content of each cell is currently massively undersampled by scRNA-seq, resulting in 'zero' expression values for the majority of genes in a large fraction of cells. MAGIC circumvents this problem by using a diffusion process along the data to share information between similar cells and thereby denoise and impute expression values. In addition to MAGIC, I apply archetype analysis to study various cellular stages during EMT, and I find novel biological processes in the previously unstudied intermediate states. The work presented here introduces a mathematical modeling framework and advanced geometric tools to analyze single-cell data. These ideas can be generally applied to various biological systems. Here, I apply them to answer important biological questions in T cell differentiation and EMT. The obtained knowledge has applications in our basic understanding of the process of EMT, T cell biology and in cancer treatment

Proliferation tracing with single-cell mass cytometry optimizes generation of stem cell memory-like T cells.

Selective differentiation of naive T cells into multipotent T cells is of great interest clinically for the generation of cell-based cancer immunotherapies. Cellular differentiation depends crucially on division state and time. Here we adapt a dye dilution assay for tracking cell proliferative history through mass cytometry and uncouple division, time and regulatory protein expression in single naive human T cells during their activation and expansion in a complex ex vivo milieu. Using 23 markers, we defined groups of proteins controlled predominantly by division state or time and found that undivided cells account for the majority of phenotypic diversity. We next built a map of cell state changes during naive T-cell expansion. By examining cell signaling on this map, we rationally selected ibrutinib, a BTK and ITK inhibitor, and administered it before T cell activation to direct differentiation toward a T stem cell memory (TSCM)-like phenotype. This method for tracing cell fate across division states and time can be broadly applied for directing cellular differentiation

Why one-size-fits-all vaso-modulatory interventions fail to control glioma invasion: in silico insights

There is an ongoing debate on the therapeutic potential of vaso-modulatory interventions against glioma invasion. Prominent vasculature-targeting therapies involve functional tumour-associated blood vessel deterioration and normalisation. The former aims at tumour infarction and nutrient deprivation medi- ated by vascular targeting agents that induce occlusion/collapse of tumour blood vessels. In contrast, the therapeutic intention of normalising the abnormal structure and function of tumour vascular net- works, e.g. via alleviating stress-induced vaso-occlusion, is to improve chemo-, immuno- and radiation therapy efficacy. Although both strategies have shown therapeutic potential, it remains unclear why they often fail to control glioma invasion into the surrounding healthy brain tissue. To shed light on this issue, we propose a mathematical model of glioma invasion focusing on the interplay between the mi- gration/proliferation dichotomy (Go-or-Grow) of glioma cells and modulations of the functional tumour vasculature. Vaso-modulatory interventions are modelled by varying the degree of vaso-occlusion. We discovered the existence of a critical cell proliferation/diffusion ratio that separates glioma invasion re- sponses to vaso-modulatory interventions into two distinct regimes. While for tumours, belonging to one regime, vascular modulations reduce the tumour front speed and increase the infiltration width, for those in the other regime the invasion speed increases and infiltration width decreases. We show how these in silico findings can be used to guide individualised approaches of vaso-modulatory treatment strategies and thereby improve success rates

State Differentiation by Transient Truncation in Coupled Threshold Dynamics

Dynamics with a threshold input--output relation commonly exist in gene, signal-transduction, and neural networks. Coupled dynamical systems of such threshold elements are investigated, in an effort to find differentiation of elements induced by the interaction. Through global diffusive coupling, novel states are found to be generated that are not the original attractor of single-element threshold dynamics, but are sustained through the interaction with the elements located at the original attractor. This stabilization of the novel state(s) is not related to symmetry breaking, but is explained as the truncation of transient trajectories to the original attractor due to the coupling. Single-element dynamics with winding transient trajectories located at a low-dimensional manifold and having turning points are shown to be essential to the generation of such novel state(s) in a coupled system. Universality of this mechanism for the novel state generation and its relevance to biological cell differentiation are briefly discussed.Comment: 8 pages. Phys. Rev. E. in pres

Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics.

BackgroundSingle-cell transcriptomics allows researchers to investigate complex communities of heterogeneous cells. It can be applied to stem cells and their descendants in order to chart the progression from multipotent progenitors to fully differentiated cells. While a variety of statistical and computational methods have been proposed for inferring cell lineages, the problem of accurately characterizing multiple branching lineages remains difficult to solve.ResultsWe introduce Slingshot, a novel method for inferring cell lineages and pseudotimes from single-cell gene expression data. In previously published datasets, Slingshot correctly identifies the biological signal for one to three branching trajectories. Additionally, our simulation study shows that Slingshot infers more accurate pseudotimes than other leading methods.ConclusionsSlingshot is a uniquely robust and flexible tool which combines the highly stable techniques necessary for noisy single-cell data with the ability to identify multiple trajectories. Accurate lineage inference is a critical step in the identification of dynamic temporal gene expression