Publications
Some recent publications are shown below. For a full list of publications, please visit Google Scholar.
PREPRINT; Spatial transcriptomics reveals ovarian cancer subclones with distinct tumour microenvironments
High-grade serous ovarian carcinoma (HGSOC) is characterised by recurrence, chemotherapy resistance and overall poor prognosis. Genetic heterogeneity of tumour cells and the microenvironment of the tumour have been hypothesised as key determinants of treatment resistance and relapse. Here, using a combination of spatial and single cell transcriptomics (10x Visium and Chromium platforms), we examine tumour genetic heterogeneity and infiltrating populations of HGSOC samples from eight patients with variable response to neoadjuvant chemotherapy. By inferring gross copy number alterations (CNAs), we identified distinct tumour subclones co-existing within individual tumour sections. These tumour subclones have unique CNA profiles and spatial locations within each tumour section, which were further validated by ultra-low-pass whole genome sequencing. Differential expression analysis between subclones within the same section identified both tumour cell intrinsic expression differences and markers indicative of different infiltrating cell populations. The gene sets differentially expressed between subclones were significantly enriched for genes encoding plasma membrane and secreted proteins, indicative of subclone-specific microenvironments. Furthermore, we identified tumour derived ligands with variable expression levels between subclones that correlated or anticorrelated with various non-malignant cell infiltration patterns. We highlight several of these that are potentially direct tumour-stroma/immune cell relationships as the non-malignant cell type expresses a cognate receptor for the tumour derived ligand. These include predictions of CXCL10-CXCR3 mediated recruitment of T and B cells to associate with the subclones of one patient and CD47-SIRPA mediated exclusion of macrophages from association with subclones of another. Finally, we show that published HGSOC molecular subtype signatures associated with prognosis are heterogeneously expressed across tumour sections and that areas containing different tumour subclones with different infiltration patterns can match different subtypes. Our study highlights the high degree of intratumoural subclonal and infiltrative heterogeneity in HGSOC which will be critical to better understand resistance and relapse.
Systematic assessment of tissue dissociation and storage biases in single-cell and single-nucleus RNA-seq workflows.
Background Single-cell RNA sequencing has been widely adopted to estimate the cellular composition of heterogeneous tissues and obtain transcriptional profiles of individual cells. Multiple approaches for optimal sample dissociation and storage of single cells have been proposed as have single-nuclei profiling methods. What has been lacking is a systematic comparison of their relative biases and benefits. Results Here, we compare gene expression and cellular composition of single-cell suspensions prepared from adult mouse kidney using two tissue dissociation protocols. For each sample, we also compare fresh cells to cryopreserved and methanol-fixed cells. Lastly, we compare this single-cell data to that generated using three single-nucleus RNA sequencing workflows. Our data confirms prior reports that digestion on ice avoids the stress response observed with 37 degrees C dissociation. It also reveals cell types more abundant either in the cold or warm dissociations that may represent populations that require gentler or harsher conditions to be released intact. For cell storage, cryopreservation of dissociated cells results in a major loss of epithelial cell types; in contrast, methanol fixation maintains the cellular composition but suffers from ambient RNA leakage. Finally, cell type composition differences are observed between single-cell and single-nucleus RNA sequencing libraries. In particular, we note an underrepresentation of T, B, and NK lymphocytes in the single-nucleus libraries. Conclusions Systematic comparison of recovered cell types and their transcriptional profiles across the workflows has highlighted protocol-specific biases and thus enables researchers starting single-cell experiments to make an informed choice.
Predicting cell-to-cell communication networks using NATMI
Development of high throughput single-cell sequencing technologies has made it cost-effective to profile thousands of cells from diverse samples containing multiple cell types. To study how these different cell types work together, here we develop NATMI (Network Analysis Toolkit for Multicellular Interactions). NATMI uses connectomeDB2020 (a database of 2293 manually curated ligand-receptor pairs with literature support) to predict and visualise cell-to-cell communication networks from single-cell (or bulk) expression data. Using multiple published single-cell datasets we demonstrate how NATMI can be used to identify (i) the cell-type pairs that are communicating the most (or most specifically) within a network, (ii) the most active (or specific) ligand-receptor pairs active within a network, (iii) putative highly-communicating cellular communities and (iv) differences in intercellular communication when profiling given cell types under different conditions. Furthermore, analysis of the Tabula Muris (organism-wide) atlas confirms our previous prediction that autocrine signalling is a major feature of cell-to-cell communication networks, while also revealing that hundreds of ligands and their cognate receptors are co-expressed in individual cells suggesting a substantial potential for self-signalling. Single cell expression data allows for inferring cell-cell communication between cells expressing ligands and those expressing their cognate receptors. Here the authors present an updated and curated database of ligand-receptor pairs and a Python-based toolkit to construct and analyse communication networks from single cell and bulk expression data.
An atlas of human long non-coding RNAs with accurate 5′ ends.
Long non-coding RNAs (lncRNAs) are largely heterogeneous and functionally uncharacterized. Here, using FANTOM5 cap analysis of gene expression (CAGE) data, we integrate multiple transcript collections to generate a comprehensive atlas of 27,919 human lncRNA genes with high-confidence 5' ends and expression profiles across 1,829 samples from the major human primary cell types and tissues. Genomic and epigenomic classification of these lncRNAs reveals that most intergenic lncRNAs originate from enhancers rather than from promoters. Incorporating genetic and expression data, we show that lncRNAs overlapping trait-associated single nucleotide polymorphisms are specifically expressed in cell types relevant to the traits, implicating these lncRNAs in multiple diseases. We further demonstrate that lncRNAs overlapping expression quantitative trait loci (eQTL)-associated single nucleotide polymorphisms of messenger RNAs are co-expressed with the corresponding messenger RNAs, suggesting their potential roles in transcriptional regulation. Combining these findings with conservation data, we identify 19,175 potentially functional lncRNAs in the human genome.