Search

Abstract Recent innovations in single-cell RNA-sequencing (scRNA-seq) provide the technology to investigate biological questions at cellular resolution. Pooling cells from multiple individuals has become a common strategy, and droplets can subsequently be assigned to a specific individual by leveraging their inherent genetic differences. An implicit challenge with scRNA-seq is the occurrence of doublets—droplets containing two or more cells. We develop Demuxafy, a framework to enhance donor assignment and doublet removal through the consensus intersection of multiple demultiplexing and doublet detecting methods. Demuxafy significantly improves droplet assignment by separating singlets from doublets and classifying the correct individual.

Purpose: Genome-wide association studies have recently uncovered many loci associated with variation in intraocular pressure (IOP). Artificial intelligence (AI) can be used to interrogate the effect of specific genetic knockouts on the morphology of trabecular meshwork cells (TMCs) and thus, IOP regulation. Design: Experimental study. Subjects: Primary TMCs collected from human donors. Methods: Sixty-two genes at 55 loci associated with IOP variation were knocked out in primary TMC lines. All cells underwent high-throughput microscopy imaging after being stained with a 5-channel fluorescent cell staining protocol. A convolutional neural network was trained to distinguish between gene knockout and normal control cell images. The area under the receiver operator curve (AUC) metric was used to quantify morphological variation in gene knockouts to identify potential pathological perturbations. Main Outcome Measures: Degree of morphological variation as measured by deep learning algorithm accuracy of differentiation from normal controls. Results: Cells where LTBP2 or BCAS3 had been perturbed demonstrated the greatest morphological variation from normal TMCs (AUC 0.851, standard deviation [SD] 0.030; and AUC 0.845, SD 0.020, respectively). Of 7 multigene loci, 5 had statistically significant differences in AUC (P

Abstract The mechanisms by which DNA alleles contribute to disease risk, drug response, and other human phenotypes are highly context-specific, varying across cell types and different conditions. Human induced pluripotent stem cells are uniquely suited to study these context-dependent effects but cell lines from hundreds or thousands of individuals are required. Village cultures, where multiple induced pluripotent stem lines are cultured and differentiated in a single dish, provide an elegant solution for scaling induced pluripotent stem experiments to the necessary sample sizes required for population-scale studies. Here, we show the utility of village models, demonstrating how cells can be assigned to an induced pluripotent stem line using single-cell sequencing and illustrating that the genetic, epigenetic or induced pluripotent stem line-specific effects explain a large percentage of gene expression variation for many genes. We demonstrate that village methods can effectively detect induced pluripotent stem line-specific effects, including sensitive dynamics of cell states.

Abstract Complex diseases are heterogenous due to variation in their genetic and environmental underpinnings, leading to varied treatment responses. Genome‐wide association studies (GWAS) integrated with single‐cell expression quantitative trait loci analyses (eQTL) can pinpoint cell‐type specific candidate disease‐relevant genes and pathways. This knowledge can be applied to patient stratification and novel therapeutic target identification. Here, we describe the translational potential of cell‐type specific genetic regulation, using Crohn's disease as an example.

Summary: The hypothalamus is a region of the brain that plays an important role in regulating body functions and behaviors. There is a growing interest in human pluripotent stem cells (hPSCs) for modeling diseases that affect the hypothalamus. Here, we established an hPSC-derived hypothalamus organoid differentiation protocol to model the cellular diversity of this brain region. Using an hPSC line with a tyrosine hydroxylase (TH)-TdTomato reporter for dopaminergic neurons (DNs) and other TH-expressing cells, we interrogated DN-specific pathways and functions in electrophysiologically active hypothalamus organoids. Single-cell RNA sequencing (scRNA-seq) revealed diverse neuronal and non-neuronal cell types in mature hypothalamus organoids. We identified several molecularly distinct hypothalamic DN subtypes that demonstrated different developmental maturities. Our in vitro 3D hypothalamus differentiation protocol can be used to study the development of this critical brain structure and can be applied to disease modeling to generate novel therapeutic approaches for disorders centered around the hypothalamus.

Abstract Using latent variables in gene expression data can help correct unobserved confounders and increase statistical power for expression quantitative trait Loci (eQTL) detection. The probabilistic estimation of expression residuals (PEER) and principal component analysis (PCA) are widely used methods that can remove unwanted variation and improve eQTL discovery power in bulk RNA-seq analysis. However, their performance has not been evaluated extensively in single-cell eQTL analysis, especially for different cell types. Potential challenges arise due to the structure of single-cell RNA-seq data, including sparsity, skewness, and mean-variance relationship. Here, we show by a series of analyses that PEER and PCA require additional quality control and data transformation steps on the pseudo-bulk matrix to obtain valid latent variables; otherwise, it can result in highly correlated factors (Pearson's correlation r = 0.63 ~ 0.99). Incorporating valid PFs/PCs in the eQTL association model would identify 1.7 ~ 13.3% more eGenes. Sensitivity analysis showed that the pattern of change between the number of eGenes detected and fitted PFs/PCs varied significantly in different cell types. In addition, using highly variable genes to generate latent variables could achieve similar eGenes discovery power as using all genes but save considerable computational resources (~ 6.2-fold faster).

Age-related macular degeneration (AMD) is a leading cause of vision loss, and there is no approved treatment for AMD with geographic atrophy. Here, the authors used transcriptomic and proteomic analyses of patient induced pluripotent stem cell-derived retinal pigment epithelium to better understand disease mechanisms.

Abstract Background While single-cell transcriptional profiling has greatly increased our capacity to interrogate biology, accurate cell classification within and between datasets is a key challenge. This is particularly so in pluripotent stem cell-derived organoids which represent a model of a developmental system. Here, clustering algorithms and selected marker genes can fail to accurately classify cellular identity while variation in analyses makes it difficult to meaningfully compare datasets. Kidney organoids provide a valuable resource to understand kidney development and disease. However, direct comparison of relative cellular composition between protocols has proved challenging. Hence, an unbiased approach for classifying cell identity is required. Methods The R package, scPred, was trained on multiple single cell RNA-seq datasets of human fetal kidney. A hierarchical model classified cellular subtypes into nephron, stroma and ureteric epithelial elements. This model, provided in the R package DevKidCC ( github.com/KidneyRegeneration/DevKidCC ), was then used to predict relative cell identity within published kidney organoid datasets generated using distinct cell lines and differentiation protocols, interrogating the impact of such variations. The package contains custom functions for the display of differential gene expression within cellular subtypes. Results DevKidCC was used to directly compare between distinct kidney organoid protocols, identifying differences in relative proportions of cell types at all hierarchical levels of the model and highlighting variations in stromal and unassigned cell types, nephron progenitor prevalence and relative maturation of individual epithelial segments. Of note, DevKidCC was able to distinguish distal nephron from ureteric epithelium, cell types with overlapping profiles that have previously confounded analyses. When applied to a variation in protocol via the addition of retinoic acid, DevKidCC identified a consequential depletion of nephron progenitors. Conclusions The application of DevKidCC to kidney organoids reproducibly classifies component cellular identity within distinct single-cell datasets. The application of the tool is summarised in an interactive Shiny application, as are examples of the utility of in-built functions for data presentation. This tool will enable the consistent and rapid comparison of kidney organoid protocols, driving improvements in patterning to kidney endpoints and validating new approaches.

Abstract Background Advances in droplet-based single-cell RNA-sequencing (scRNA-seq) have dramatically increased throughput, allowing tens of thousands of cells to be routinely sequenced in a single experiment. In addition to cells, droplets capture cell-free “ambient” RNA predominantly caused by lysis of cells during sample preparation. Samples with high ambient RNA concentration can create challenges in accurately distinguishing cell-containing droplets and droplets containing ambient RNA. Current methods to separate these groups often retain a significant number of droplets that do not contain cells or empty droplets. Additionally, there are currently no methods available to detect droplets containing damaged cells, which comprise partially lysed cells, the original source of the ambient RNA. Results Here, we describe DropletQC, a new method that is able to detect empty droplets, damaged, and intact cells, and accurately distinguish them from one another. This approach is based on a novel quality control metric, the nuclear fraction, which quantifies for each droplet the fraction of RNA originating from unspliced, nuclear pre-mRNA. We demonstrate how DropletQC provides a powerful extension to existing computational methods for identifying empty droplets such as EmptyDrops. Conclusions We implement DropletQC as an R package, which can be easily integrated into existing single-cell analysis workflows.

Summary: To assess the transcriptomic profile of disease-specific cell populations, fibroblasts from patients with primary open-angle glaucoma (POAG) were reprogrammed into induced pluripotent stem cells (iPSCs) before being differentiated into retinal organoids and compared with those from healthy individuals. We performed single-cell RNA sequencing of a total of 247,520 cells and identified cluster-specific molecular signatures. Comparing the gene expression profile between cases and controls, we identified novel genetic associations for this blinding disease. Expression quantitative trait mapping identified a total of 4,443 significant loci across all cell types, 312 of which are specific to the retinal ganglion cell subpopulations, which ultimately degenerate in POAG. Transcriptome-wide association analysis identified genes at loci previously associated with POAG, and analysis, conditional on disease status, implicated 97 statistically significant retinal ganglion cell-specific expression quantitative trait loci. This work highlights the power of large-scale iPSC studies to uncover context-specific profiles for a genetically complex disease.

Abstract Background High throughput single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for exploring cellular heterogeneity among complex human cancers. scRNA-Seq studies using fresh human surgical tissue are logistically difficult, preclude histopathological triage of samples, and limit the ability to perform batch processing. This hindrance can often introduce technical biases when integrating patient datasets and increase experimental costs. Although tissue preservation methods have been previously explored to address such issues, it is yet to be examined on complex human tissues, such as solid cancers and on high throughput scRNA-Seq platforms. Methods Using the Chromium 10X platform, we sequenced a total of ~ 120,000 cells from fresh and cryopreserved replicates across three primary breast cancers, two primary prostate cancers and a cutaneous melanoma. We performed detailed analyses between cells from each condition to assess the effects of cryopreservation on cellular heterogeneity, cell quality, clustering and the identification of gene ontologies. In addition, we performed single-cell immunophenotyping using CITE-Seq on a single breast cancer sample cryopreserved as solid tissue fragments. Results Tumour heterogeneity identified from fresh tissues was largely conserved in cryopreserved replicates. We show that sequencing of single cells prepared from cryopreserved tissue fragments or from cryopreserved cell suspensions is comparable to sequenced cells prepared from fresh tissue, with cryopreserved cell suspensions displaying higher correlations with fresh tissue in gene expression. We showed that cryopreservation had minimal impacts on the results of downstream analyses such as biological pathway enrichment. For some tumours, cryopreservation modestly increased cell stress signatures compared to freshly analysed tissue. Further, we demonstrate the advantage of cryopreserving whole-cells for detecting cell-surface proteins using CITE-Seq, which is impossible using other preservation methods such as single nuclei-sequencing. Conclusions We show that the viable cryopreservation of human cancers provides high-quality single-cells for multi-omics analysis. Our study guides new experimental designs for tissue biobanking for future clinical single-cell RNA sequencing studies.

Human pluripotent stem cell (hPSC)-derived progenies are immature versions of cells, presenting a potential limitation to the accurate modelling of diseases associated with maturity or age. Hence, it is important to characterise how closely cells used in culture resemble their native counterparts. In order to select appropriate time points of retinal pigment epithelium (RPE) cultures that reflect native counterparts, we characterised the transcriptomic profiles of the hPSC-derived RPE cells from 1- and 12-month cultures. We differentiated the human embryonic stem cell line H9 into RPE cells, performed single-cell RNA-sequencing of a total of 16,576 cells to assess the molecular changes of the RPE cells across these two culture time points. Our results indicate the stability of the RPE transcriptomic signature, with no evidence of an epithelial–mesenchymal transition, and with the maturing populations of the RPE observed with time in culture. Assessment of Gene Ontology pathways revealed that as the cultures age, RPE cells upregulate expression of genes involved in metal binding and antioxidant functions. This might reflect an increased ability to handle oxidative stress as cells mature. Comparison with native human RPE data confirms a maturing transcriptional profile of RPE cells in culture. These results suggest that long-term in vitro culture of RPE cells allows the modelling of specific phenotypes observed in native mature tissues. Our work highlights the transcriptional landscape of hPSC-derived RPE cells as they age in culture, which provides a reference for native and patient samples to be benchmarked against.

Abstract Background The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) has provided a foundation for in vitro human disease modelling, drug development and population genetics studies. Gene expression plays a critical role in complex disease risk and therapeutic response. However, while the genetic background of reprogrammed cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells which would provide significant resolution. By integrating single cell RNA-sequencing (scRNA-seq) and population genetics, we apply a framework in which to evaluate cell type-specific effects of genetic variation on gene expression. Results Here, we perform scRNA-seq on 64,018 fibroblasts from 79 donors and map expression quantitative trait loci (eQTLs) at the level of individual cell types. We demonstrate that the majority of eQTLs detected in fibroblasts are specific to an individual cell subtype. To address if the allelic effects on gene expression are maintained following cell reprogramming, we generate scRNA-seq data in 19,967 iPSCs from 31 reprogramed donor lines. We again identify highly cell type-specific eQTLs in iPSCs and show that the eQTLs in fibroblasts almost entirely disappear during reprogramming. Conclusions This work provides an atlas of how genetic variation influences gene expression across cell subtypes and provides evidence for patterns of genetic architecture that lead to cell type-specific eQTL effects.

Inferring cell type proportions from transcriptomics data is affected by data transformation, normalization, choice of method and the markers used. Here, the authors use single-cell RNAseq datasets to evaluate the impact of these factors and propose guidelines to maximise deconvolution performance.

Healthspan and healthy aging are areas of research with potential socioeconomic impact. Here, the authors present the Medical Genome Reference Bank (MGRB) which consist of over 4,000 individuals aged 70 years and older without a history of the major age-related diseases and report on results from whole-genome sequencing and association analyses.

Abstract Single-cell RNA sequencing has enabled the characterization of highly specific cell types in many tissues, as well as both primary and stem cell-derived cell lines. An important facet of these studies is the ability to identify the transcriptional signatures that define a cell type or state. In theory, this information can be used to classify an individual cell based on its transcriptional profile. Here, we present scPred, a new generalizable method that is able to provide highly accurate classification of single cells, using a combination of unbiased feature selection from a reduced-dimension space, and machine-learning probability-based prediction method. We apply scPred to scRNA-seq data from pancreatic tissue, mononuclear cells, colorectal tumor biopsies, and circulating dendritic cells and show that scPred is able to classify individual cells with high accuracy. The generalized method is available at https://github.com/powellgenomicslab/scPred/.

Abstract A variety of methods have been developed to demultiplex pooled samples in a single cell RNA sequencing (scRNA-seq) experiment which either require hashtag barcodes or sample genotypes prior to pooling. We introduce scSplit which utilizes genetic differences inferred from scRNA-seq data alone to demultiplex pooled samples. scSplit also enables mapping clusters to original samples. Using simulated, merged, and pooled multi-individual datasets, we show that scSplit prediction is highly concordant with demuxlet predictions and is highly consistent with the known truth in cell-hashing dataset. scSplit is ideally suited to samples without external genotype information and is available at: https://github.com/jon-xu/scSplit

Summary: Langerhans cells (LC) are skin-resident antigen-presenting cells that regulate immune responses to epithelial microorganisms. Human papillomavirus (HPV) infection can promote malignant epithelial transformation. As LCs are considered important for controlling HPV infection, we compared the transcriptome of murine LCs from skin transformed by K14E7 oncoprotein and from healthy skin. We identified transcriptome heterogeneity at the single cell level amongst LCs in normal skin, associated with ontogeny, cell cycle, and maturation. We identified a balanced co-existence of immune-stimulatory and immune-inhibitory LC cell states in normal skin that was significantly disturbed in HPV16 E7-transformed skin. Hyperplastic skin was depleted of immune-stimulatory LCs and enriched for LCs with an immune-inhibitory gene signature, and LC-keratinocyte crosstalk was dysregulated. We identified reduced expression of interleukin (IL)-34, a critical molecule for LC homeostasis. Enrichment of an immune-inhibitory LC gene signature and reduced levels of epithelial IL-34 were also found in human HPV-associated cervical epithelial cancers.

Summary: DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Finding cell states and their transcriptional relatedness is a main outcome from analysing single-cell data. In developmental biology, determining whether cells are related in a differentiation lineage remains a major challenge. A seamless analysis pipeline from cell clustering to estimating the probability of transitions between cell clusters is lacking. Here, we present Single Cell Global fate Potential of Subpopulations (scGPS) to characterise transcriptional relationship between cell states. scGPS decomposes mixed cell populations in one or more samples into clusters (SCORE algorithm) and estimates pairwise transitioning potential (scGPS algorithm) of any pair of clusters. SCORE allows for the assessment and selection of stable clustering results, a major challenge in clustering analysis. scGPS implements a novel approach, with machine learning classification, to flexibly construct trajectory connections between clusters. scGPS also has a feature selection functionality by network and modelling approaches to find biological processes and driver genes that connect cell populations. We applied scGPS in diverse developmental contexts and show superior results compared to a range of clustering and trajectory analysis methods. scGPS is able to identify the dynamics of cellular plasticity in a user-friendly workflow, that is fast and memory efficient. scGPS is implemented in R with optimised functions using C++ and is publicly available in Bioconductor.

Summary: We assessed the pluripotency of human induced pluripotent stem cells (iPSCs) maintained on an automated platform using StemFlex and TeSR-E8 media. Analysis of transcriptome of single cells revealed similar expression of core pluripotency genes, as well as genes associated with naive and primed states of pluripotency. Analysis of individual cells from four samples consisting of two different iPSC lines each grown in the two culture media revealed a shared subpopulation structure with three main subpopulations different in pluripotency states. By implementing a machine learning approach, we estimated that most cells within each subpopulation are very similar between all four samples. The single-cell RNA sequencing analysis of iPSC lines grown in both media reports the molecular signature in StemFlex medium and how it compares to that observed in the TeSR-E8 medium. : Stem Cells Research; Transcriptomics; Automation Subject Areas: Stem Cells Research, Transcriptomics, Automation

Abstract Background Giant cell arteritis (GCA) is the most common form of vasculitis affecting elderly people. It is one of the few true ophthalmic emergencies but symptoms and signs are variable thereby making it a challenging disease to diagnose. A temporal artery biopsy is the gold standard to confirm GCA, but there are currently no specific biochemical markers to aid diagnosis. We aimed to identify a less invasive method to confirm the diagnosis of GCA, as well as to ascertain clinically relevant predictive biomarkers by studying the transcriptome of purified peripheral CD4+ and CD8+ T lymphocytes in patients with GCA. Methods We recruited 16 patients with histological evidence of GCA at the Royal Victorian Eye and Ear Hospital, Melbourne, Australia, and aimed to collect blood samples at six time points: acute phase, 2–3 weeks, 6–8 weeks, 3 months, 6 months and 12 months after clinical diagnosis. CD4+ and CD8+ T-cells were positively selected at each time point through magnetic-assisted cell sorting. RNA was extracted from all 195 collected samples for subsequent RNA sequencing. The expression profiles of patients were compared to those of 16 age-matched controls. Results Over the 12-month study period, polynomial modelling analyses identified 179 and 4 statistically significant transcripts with altered expression profiles (FDR

Covariance of gene expression pairs is due to a combination of shared genetic and environmental factors. Here the authors estimate the genetic correlation between highly heritable pairs and identify transcription factor control and chromatin interactions as possible mechanisms of correlation.

Expression quantitative trait locus (eQTL) detection has emerged as an important tool for unraveling of the relationship between genetic risk factors and disease or clinical phenotypes. Most studies use single marker linear regression to discover primary signals, followed by sequential conditional modeling to detect secondary genetic variants affecting gene expression. However, this approach assumes that functional variants are sparsely distributed and that close linkage between them has little impact on estimation of their precise location and the magnitude of effects. We describe a series of simulation studies designed to evaluate the impact of linkage disequilibrium (LD) on the fine mapping of causal variants with typical eQTL effect sizes. In the presence of multisite regulation, even though between 80 and 90% of modeled eSNPs associate with normally distributed traits, up to 10% of all secondary signals could be statistical artifacts, and at least 5% but up to one-quarter of credible intervals of SNPs within r2 > 0.8 of the peak may not even include a causal site. The Bayesian methods eCAVIAR and DAP (Deterministic Approximation of Posteriors) provide only modest improvement in resolution. Given the strong empirical evidence that gene expression is commonly regulated by more than one variant, we conclude that the fine mapping of causal variants needs to be adjusted for multisite influences, as conditional estimates can be highly biased by interference among linked sites, but ultimately experimental verification of individual effects is needed. Presumably similar conclusions apply not just to eQTL mapping, but to multisite influences on fine mapping of most types of quantitative trait.

Summary: The cellular and molecular profiles that govern the endothelial heterogeneity of the circulatory system have yet to be elucidated. Using a data-driven approach to study the endothelial compartment via single-cell RNA sequencing, we characterized cell subpopulations within and assigned them to a defined endothelial hierarchy. We show that two transcriptionally distinct endothelial populations exist within the aorta and, using two independent trajectory analysis methods, confirm that they represent transitioning cells rather than discrete cell types. Gene co-expression analysis revealed crucial regulatory networks underlying each population, including significant metabolic gene networks in progenitor cells. Using mitochondrial activity assays and phenotyping, we confirm that endovascular progenitors display higher mitochondrial content compared to differentiated endothelial cells. The identities of these populations were further validated against bulk RNA sequencing (RNA-seq) data obtained from normal and tumor-derived vasculature. Our findings validate the heterogeneity of the aortic endothelium and previously suggested hierarchy between progenitor and differentiated cells. : Lukowski et al. demonstrate the existence of two distinct endothelial cell populations in the aorta via single-cell sequencing. The data confirm that a progenitor population transitions to a mature endothelial cell, defining an endothelial hierarchy. Keywords: single-cell RNA sequencing, transcriptome, endothelial cell, aorta, hierarchy, endovascular progenitor

A Correction to this paper has been published: https://doi.org/10.1038/s41467-020-20288-9

Ataxia-telangiectasia (A-T) is a rare genetic disorder caused by loss of function of the ataxia-telangiectasia-mutated kinase and is characterized by a predisposition to cancer, pulmonary disease, immune deficiency and progressive degeneration of the cerebellum. As animal models do not faithfully recapitulate the neurological aspects, it remains unclear whether cerebellar degeneration is a neurodevelopmental or neurodegenerative phenotype. To address the necessity for a human model, we first assessed a previously published protocol for the ability to generate cerebellar neuronal cells, finding it gave rise to a population of precursors highly enriched for markers of the early hindbrain such as EN1 and GBX2, and later more mature cerebellar markers including PTF1α, MATH1, HOXB4, ZIC3, PAX6, and TUJ1. RNA sequencing was used to classify differentiated cerebellar neurons generated from integration-free A-T and control induced pluripotent stem cells. Comparison of RNA sequencing data with datasets from the Allen Brain Atlas reveals in vitro-derived cerebellar neurons are transcriptionally similar to discrete regions of the human cerebellum, and most closely resemble the cerebellum at 22 weeks post-conception. We show that patient-derived cerebellar neurons exhibit disrupted gene regulatory networks associated with synaptic vesicle dynamics and oxidative stress, offering the first molecular insights into early cerebellar pathogenesis of ataxia-telangiectasia.