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°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 singlenucleus 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.

Single-cell and single-nucleus RNA sequencing have been widely adopted in studies of heterogeneous tissues to estimate their cellular composition and obtain transcriptional profiles of individual cells. However, the current fragmentary understanding of artefacts introduced by sample preparation protocols impedes the selection of optimal workflows and compromises data interpretation. To bridge this gap, we compared performance of several workflows applied to adult mouse kidneys. Our study encompasses two tissue dissociation protocols, two cell preservation methods, bulk tissue RNA sequencing, single-cell and three single-nucleus RNA sequencing workflows for the 10x Genomics Chromium platform. These experiments enable a systematic comparison of recovered cell types and their transcriptional profiles across the workflows and highlight protocol-specific biases important for the experimental design and data interpretation.

Current approaches to stage chronic liver diseases have limited utility to directly predict liver cancer risk. Here, we employed single nucleus RNA sequencing (snRNA-seq) to characterize the cellular microenvironment of healthy and chronically injured pre-malignant livers using two distinct mouse models. Analysis of 40,748 hepatic nuclei unraveled a previously uncharacterized disease-associated hepatocyte transcriptional state (daHep). These cells were absent in healthy livers, but were increasingly prevalent as chronic liver disease progressed towards hepatocarcinogenesis. Gene expression deconvolution of 1,439 human liver transcriptomes from publicly available datasets revealed that daHep frequencies highly correlate with current histopathological liver disease staging systems. Importantly, we show that high daHep levels precede carcinogenesis in mice and humans and predict a higher risk of hepatocellular carcinoma (HCC) development. This novel transcriptional signature with diagnostic and, more importantly, prognostic significance has the potential to change the way chronic liver disease patients are staged, surveilled and risk-stratified.

The alpaca classic grey phenotype is of particular interest to the industry. Until now, there were only indirect data suggesting that the KIT gene was involved in the classic grey phenotype. All exons of KIT in three black and three classic silvergrey alpacas were sequenced. Five non-synonymous SNPs were observed. There was only one SNP found that was present only in the silvergrey alpacas, and this was also the only SNP predicted to be damaging. This variant results in a change of a glycine (Gly) to an arginine (Arg) at amino acid position 126 (c.376G>A), occurring in the second Ig-like domain of the extracellular domain of KIT. Basic protein modelling predicted that this variant is likely destabilising. Therefore, an additional 488 alpacas were genotyped for this SNP using the tetra-primer amplification refractory mutation system PCR (Tetra-primer ARMS-PCR). All classic grey alpacas were observed to be heterozygous, and 99.3% of non-grey dark base colour alpacas were found to be homozygous for the wildtype allele in this position. These results confirm that the classic grey phenotype in alpacas is the result of a c.376G>A (p.Gly126Arg) SNP in exon 3 of KIT. These data also support the hypothesis that the grey phenotype is autosomal dominant and that the mutation is most likely homozygous lethal.

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.

Tumour heterogeneity is a major obstacle to the success of cancer treatment. An accurate understanding and recognition of tumour heterogeneity is critical in the clinical management of cancer patients. Here, we utilised single-cell RNA sequencing (scRNA-seq) to uncover the intra- and inter-tumoural heterogeneity of liver metastases from a patient with metastatic uveal melanoma. The two metastases analysed were largely infiltrated by non-cancerous cells with significant variability in the proportion of different cell types. Analysis of copy number variations (CNVs) showed gain of 8q and loss of 6q in both tumours, but loss of chromosome 3 was only detected in one of the tumours. SNP array revealed a uniparental isodisomy 3 in the tumour with two copies of chromosome 3, indicating a re-gain of chromosome 3 during the development of the metastatic disease. In addition, both tumours harboured subclones with additional CNVs. Pathway enrichment analysis of differentially expressed genes revealed that cancer cells in the metastasis with isodisomy 3 showed up-regulation in epithelial-mesenchymal transition and myogenesis related genes. In contrast, upregulation in interferon signalling was observed in the metastasis with monosomy 3 and increased T-cell infiltrate. This study highlights the complexity and heterogeneity of different metastases within an individual case of uveal melanoma.

No abstract

No abstract