Quantifying expression variability in single-cell RNA sequencing data

Transcriptional noise is an intrinsic feature of cell populations and plays a driving role in mammalian development, tissue homoeostasis and immune function. While expression heterogeneity, a phenotypic readout of transcriptional noise, has been broadly studied in prokaryotic model systems or by profiling individual genes, few whole-transcriptome studies in mammalian systems have been reported. The development of single-cell RNA sequencing technologies introduced powerful tools to investigate transcriptional differences between individual cells, therefore allowing the in-depth characterisation of expression variability. In this thesis, I computationally analysed single-cell RNA sequencing data to understand transcriptional variability and expanded a statistical model to avoid confounding effects when quantifying such variability. First, I profiled individual transcriptomes of CD4$^+$ T cells, identifying a global decrease in transcriptional variability upon immune activation. By extending this analysis across two sub-species of mice, I identified an evolutionarily conserved set of immune response genes for which transcriptional variability increases during ageing. I used a Bayesian modelling framework to quantify mean expression and transcriptional variability but due to a strong confounding effect between these two parameters, variability analysis was restricted to genes that are similarly expressed across the tested conditions. To address this problem, I extended the computational framework allowing the parallel assessment of changes in mean expression and variability. Within this Bayesian framework, I introduced a joint prior linking mean expression and variability parameters, which allowed a residual over-dispersion to be measured for each gene. This measure allowed me to statistically assess changes in variability even for genes with differences in mean expression between conditions. Finally, I applied the model to identify temporal changes in variability over the time-course of spermatogenesis. This unidirectional differentiation process involves several complex steps before mature sperm form from spermatogonial stem cells. When profiling changes in variability across this developmental time-course, peaks in variability are caused by rapid changes in gene expression along the differentiation trajectory. This thesis provides a deeper understanding of technical and biological factors that drive transcriptional variability and offers a basis for future research to characterise its role in health and disease.Funding was provided via the EMBL international PhD programm

Cytomapper: an R/bioconductor package for visualisation of highly multiplexed imaging data

SUMMARY: Highly multiplexed imaging technologies enable spatial profiling of dozens of biomarkers in situ. Here we describe cytomapper, a computational tool written in R, that enables visualisation of pixel- and cell-level information obtained by multiplexed imaging. To illustrate its utility, we analysed 100 images obtained by imaging mass cytometry from a cohort of type 1 diabetes patients. In addition, cytomapper includes a Shiny application that allows hierarchical gating of cells based on marker expression and visualisation of selected cells in corresponding images. AVAILABILITY AND IMPLEMENTATION: The cytomapper package can be installed via https://www.bioconductor.org/packages/release/bioc/html/cytomapper.html. Code for analysis and further instructions can be found at https://github.com/BodenmillerGroup/cytomapper_publication. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online

Correcting the Mean-Variance Dependency for Differential Variability Testing Using Single-Cell RNA Sequencing Data.

Cell-to-cell transcriptional variability in otherwise homogeneous cell populations plays an important role in tissue function and development. Single-cell RNA sequencing can characterize this variability in a transcriptome-wide manner. However, technical variation and the confounding between variability and mean expression estimates hinder meaningful comparison of expression variability between cell populations. To address this problem, we introduce an analysis approach that extends the BASiCS statistical framework to derive a residual measure of variability that is not confounded by mean expression. This includes a robust procedure for quantifying technical noise in experiments where technical spike-in molecules are not available. We illustrate how our method provides biological insight into the dynamics of cell-to-cell expression variability, highlighting a synchronization of biosynthetic machinery components in immune cells upon activation. In contrast to the uniform up-regulation of the biosynthetic machinery, CD4+ T cells show heterogeneous up-regulation of immune-related and lineage-defining genes during activation and differentiation.NE was funded by the European Molecular Biology Laboratory (EMBL) international PhD programme. ACR was funded by the MRC Skills Development Fellowship (MR/P014178/1). SR was funded by MRC grant MC_UP_0801/1. JCM was funded by core support of Cancer Research UK and EMBL. CAV was funded by The Alan Turing Institute, EPSRC grant EP/N510129/1

Staged developmental mapping and X chromosome transcriptional dynamics during mouse spermatogenesis.

Male gametes are generated through a specialised differentiation pathway involving a series of developmental transitions that are poorly characterised at the molecular level. Here, we use droplet-based single-cell RNA-Sequencing to profile spermatogenesis in adult animals and at multiple stages during juvenile development. By exploiting the first wave of spermatogenesis, we both precisely stage germ cell development and enrich for rare somatic cell-types and spermatogonia. To capture the full complexity of spermatogenesis including cells that have low transcriptional activity, we apply a statistical tool that identifies previously uncharacterised populations of leptotene and zygotene spermatocytes. Focusing on post-meiotic events, we characterise the temporal dynamics of X chromosome re-activation and profile the associated chromatin state using CUT&RUN. This identifies a set of genes strongly repressed by H3K9me3 in spermatocytes, which then undergo extensive chromatin remodelling post-meiosis, thus acquiring an active chromatin state and spermatid-specific expression

Whole-body single-cell sequencing reveals transcriptional domains in the annelid larval body.

Animal bodies comprise diverse arrays of cells. To characterise cellular identities across an entire body, we have compared the transcriptomes of single cells randomly picked from dissociated whole larvae of the marine annelid Platynereis dumerilii. We identify five transcriptionally distinct groups of differentiated cells, each expressing a unique set of transcription factors and effector genes that implement cellular phenotypes. Spatial mapping of cells into a cellular expression atlas, and wholemount in situ hybridisation of group-specific genes reveals spatially coherent transcriptional domains in the larval body, comprising e.g. apical sensory-neurosecretory cells vs. neural/epidermal surface cells. These domains represent new, basic subdivisions of the annelid body based entirely on differential gene expression, and are composed of multiple, transcriptionally similar cell types. They do not represent clonal domains, as revealed by developmental lineage analysis. We propose that the transcriptional domains that subdivide the annelid larval body represent families of related cell types that have arisen by evolutionary diversification. Their possible evolutionary conservation makes them a promising tool for evo-devo research. (167/250).KA and JM were supported by the Marie Curie COFUND programme from the European Commission and by EMBL core funding. NE, PC, VB, and DA were supported by core funding from EMBL. KA, HMV, PYB, PV were supported by the Advanced grant “Brain Evo-Devo” from the European Research Council. JCM was supported by core funding from EMBL and Cancer Research UK

An end-to-end workflow for multiplexed image processing and analysis

Simultaneous profiling of the spatial distributions of multiple biological molecules at single-cell resolution has recently been enabled by the development of highly multiplexed imaging technologies. Extracting and analyzing biologically relevant information contained in complex imaging data requires the use of a diverse set of computational tools and algorithms. Here, we report the development of a user-friendly, customizable, and interoperable workflow for processing and analyzing data generated by highly multiplexed imaging technologies. The steinbock framework supports image pre-processing, segmentation, feature extraction, and standardized data export. Each step is performed in a reproducible fashion. The imcRtools R/Bioconductor package forms the bridge between image processing and single-cell analysis by directly importing data generated by steinbock. The package further supports spatial data analysis and integrates with tools developed within the Bio-conductor project. Together, the tools described in this workflow facilitate analyses of multiplexed imaging raw data at the single-cell and spatial level

cytoviewer: an R/Bioconductor package for interactive visualization and exploration of highly multiplexed imaging data

Background: Highly multiplexed imaging enables single-cell-resolved detection of numerous biological molecules in their spatial tissue context. Interactive visualization of multiplexed imaging data is crucial at any step of data analysis to facilitate quality control and the spatial exploration of single cell features. However, tools for interactive visualization of multiplexed imaging data are not available in the statistical programming language R. Results: Here, we describe cytoviewer, an R/Bioconductor package for interactive visualization and exploration of multi-channel images and segmentation masks. The cytoviewer package supports flexible generation of image composites, allows side-by-side visualization of single channels, and facilitates the spatial visualization of single-cell data in the form of segmentation masks. As such, cytoviewer improves image and segmentation quality control, the visualization of cell phenotyping results and qualitative validation of hypothesis at any step of data analysis. The package operates on standard data classes of the Bioconductor project and therefore integrates with an extensive framework for single-cell and image analysis. The graphical user interface allows intuitive navigation and little coding experience is required to use the package. We showcase the functionality and biological application of cytoviewer by analysis of an imaging mass cytometry dataset acquired from cancer samples. Conclusions: The cytoviewer package offers a rich set of features for highly multiplexed imaging data visualization in R that seamlessly integrates with the workflow for image and single-cell data analysis. It can be installed from Bioconductor viaISSN:1471-210

Challenges in measuring and understanding biological noise.

Biochemical reactions are intrinsically stochastic, leading to variation in the production of mRNAs and proteins within cells. In the scientific literature, this source of variation is typically referred to as 'noise'. The observed variability in molecular phenotypes arises from a combination of processes that amplify and attenuate noise. Our ability to quantify cell-to-cell variability in numerous biological contexts has been revolutionized by recent advances in single-cell technology, from imaging approaches through to 'omics' strategies. However, defining, accurately measuring and disentangling the stochastic and deterministic components of cell-to-cell variability is challenging. In this Review, we discuss the sources, impact and function of molecular phenotypic variability and highlight future directions to understand its role.N.E. was supported by the European Molecular Biology Laboratory (EMBL) International PhD Programme. M.D.M. was supported by Wellcome Trust grant 105045/Z/14/Z to J.C.M. J.C.M. was supported by core funding from EMBL and Cancer Research UK (award number 17197)

Cytomapper: an R/bioconductor package for visualisation of highly multiplexed imaging data

Highly multiplexed imaging technologies enable spatial profiling of dozens of biomarkersin situ. Standard data processing pipelines quantify cell-specific features and generate object segmentation masks as well as multi-channel images. Therefore, multiplexed imaging data can be visualised across two layers of information: pixel-intensities represent the spatial expression of biomarkers across an image while segmented objects visualise cellular morphology, interactions and cell phenotypes in their microenvironment.Here we describecytomapper, a computational tool that enables visualisation of pixel- and cell-level information obtained by multiplexed imaging. The package is written in the statistical programming language R, integrates with the image and single-cell analysis infrastructure of the Bioconductor project, and allows visualisation of single to hundreds of images in parallel. Usingcytomapper, expression of multiple markers is displayed as composite images, segmentation masks are coloured based on cellular features, and selected cells can be outlined in images based on their cell type, among other functions. We illustrate the utility ofcytomapperby analysing 100 images obtained by imaging mass cytometry from a cohort of type 1 diabetes patients and healthy individuals. In addition,cytomapperincludes a Shiny application that allows hierarchical gating of cells based on marker expression and visualisation of selected cells in corresponding images. Together,cytomapperoffers tools for diverse image and single-cell visualisation approaches and supports robust cell phenotyping via gating

Multiplexed imaging mass cytometry of the chemokine milieus in melanoma characterizes features of the response to immunotherapy

Intratumoral immune cells are crucial for tumor control and anti-tumor responses during immunotherapy. Immune cell trafficking into tumors is mediated by chemokines, which are expressed and secreted upon various stimuli and interact with specific receptors. To broadly characterize chemokine expression and function in tumors, we have used multiplex mass cytometry-based imaging of protein markers and RNA transcripts to analyze the chemokine landscape and immune infiltration in metastatic melanoma samples. Tumors that lacked immune infiltration were devoid of most chemokines and exhibited particularly low levels of antigen presentation and inflammation. Infiltrated tumors were characterized by expression of multiple chemokines. CXCL9 and CXCL10 were often localized in patches associated with dysfunctional T cells expressing CXCL13 which was strongly associated with B cell patches and follicles. TCF7+ naïve-like T cells, which predict response to immunotherapy, were enriched in the vicinity of B cell patches and follicles. Our data highlight the strength of RNA and protein co-detection which was critical to deconvolve specialized immune microenvironments in inflamed tumors based on chemokine expression. Our findings further suggest that the formation of tertiary lymphoid structures is accompanied by naïve and naive- like T cell recruitment, which ultimately boosts anti-tumor activity