Anne Senabouth

ascend: R package for analysis of single cell RNA-seq data

Summary ascend is an R package comprised of fast, streamlined analysis functions optimized to address the statistical challenges of single cell RNA-seq. The package incorporates novel and established methods to provide a flexible framework to perform filtering, quality control, normalization, dimension reduction, clustering, differential expression and a wide-range of plotting. ascend is designed to work with scRNA-seq data generated by any high-throughput platform, and includes functions to convert data objects between software packages. Availability The R package and associated vignettes are freely available at https://github.com/IMB-Computational-Genomics-Lab/ascend. Contact joseph.powell@uq.edu.au Supplementary information An example dataset is available at ArrayExpress, accession number E-MTAB-6108

Single cell RNA sequencing of stem cell-derived retinal ganglion cells

We used single cell sequencing technology to characterize the transcriptomes of 1,174 human embryonic stem cell-derived retinal ganglion cells (RGCs) at the single cell level. The human embryonic stem cell line BRN3B-mCherry (A81-H7), was differentiated to RGCs using a guided differentiation approach. Cells were harvested at day 36 and prepared for single cell RNA sequencing. Our data indicates the presence of three distinct subpopulations of cells, with various degrees of maturity. One cluster of 288 cells showed increased expression of genes involved in axon guidance together with semaphorin interactions, cell-extracellular matrix interactions and ECM proteoglycans, suggestive of a more mature RGC phenotype.

Single-cell RNA-seq of human induced pluripotent stem cells reveals cellular heterogeneity and cell state transitions between subpopulations

Heterogeneity of cell states represented in pluripotent cultures has not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing can be used to identify how individual pluripotent cells function. Here, we present results from the analysis of single-cell RNA sequencing data from 18,787 individual WTC-CRISPRi human induced pluripotent stem cells. We developed an unsupervised clustering method and, through this, identified four subpopulations distinguishable on the basis of their pluripotent state, including a core pluripotent population (48.3%), proliferative (47.8%), early primed for differentiation (2.8%), and late primed for differentiation (1.1%). For each subpopulation, we were able to identify the genes and pathways that define differences in pluripotent cell states. Our method identified four transcriptionally distinct predictor gene sets composed of 165 unique genes that denote the specific pluripotency states; using these sets, we developed a multigenic machine learning prediction method to accurately classify single cells into each of the subpopulations. Compared against a set of established pluripotency markers, our method increases prediction accuracy by 10%, specificity by 20%, and explains a substantially larger proportion of deviance (up to threefold) from the prediction model. Finally, we developed an innovative method to predict cells transitioning between subpopulations and support our conclusions with results from two orthogonal pseudotime trajectory methods.

Determining cell fate specification and genetic contribution to cardiac disease risk in hiPSC-derived cardiomyocytes at single cell resolution

The majority of genetic loci underlying common disease risk act through changing genome regulation, and are routinely linked to expression quantitative trait loci, where gene expression is measured using bulk populations of mature cells. A crucial step that is missing is evidence of variation in the expression of these genes as cells progress from a pluripotent to mature state. This is especially important for cardiovascular disease, as the majority of cardiac cells have limited properties for renewal postneonatal. To investigate the dynamic changes in gene expression across the cardiac lineage, we generated RNA-sequencing data captured from 43,168 single cells progressing through in vitro cardiac-directed differentiation from pluripotency. We developed a novel and generalized unsupervised cell clustering approach and a machine learning method for prediction of cell transition. Using these methods, we were able to reconstruct the cell fate choices as cells transition from a pluripotent state to mature cardiomyocytes, uncovering intermediate cell populations that do not progress to maturity, and distinct cell trajectories that terminate in cardiomyocytes that differ in their contractile forces. Second, we identify new gene markers that denote lineage specification and demonstrate a substantial increase in their utility for cell identification over current pluripotent and cardiogenic markers. By integrating results from analysis of the single cell lineage RNA-sequence data with population-based GWAS of cardiovascular disease and cardiac tissue eQTLs, we show that the pathogenicity of disease-associated genes is highly dynamic as cells transition across their developmental lineage, and exhibit variation between cell fate trajectories. Through the integration of single cell RNA-sequence data with population-scale genetic data we have identified genes significantly altered at cell specification events providing insights into a context-dependent role in cardiovascular disease risk. This study provides a valuable data resource focused on in vitro cardiomyocyte differentiation to understand cardiac disease coupled with new analytical methods with broad applications to single-cell data.

Single-Cell Profiling Identifies Key Pathways Expressed by iPSCs Cultured in Different Commercial Media

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.

Generation of human neural retina transcriptome atlas by single cell RNA sequencing

The retina is a highly specialized neural tissue that senses light and initiates image processing. Although the functional organisation of specific cells within the retina has been well-studied, the molecular profile of many cell types remains unclear in humans. To comprehensively profile cell types in the human retina, we performed single cell RNA-sequencing on 20,009 cells obtained post-mortem from three donors and compiled a reference transcriptome atlas. Using unsupervised clustering analysis, we identified 18 transcriptionally distinct clusters representing all known retinal cells: rod photoreceptors, cone photoreceptors, Müller glia cells, bipolar cells, amacrine cells, retinal ganglion cells, horizontal cells, retinal astrocytes and microglia. Notably, our data captured molecular profiles for healthy and early degenerating rod photoreceptors, and revealed a novel role of MALAT1 in putative rod degeneration. We also demonstrated the use of this retina transcriptome atlas to benchmark pluripotent stem cell-derived cone photoreceptors and an adult Müller glia cell line. This work provides an important reference with unprecedented insights into the transcriptional landscape of human retinal cells, which is fundamental to our understanding of retinal biology and disease.

Cardiac differentiation at single cell resolution reveals a requirement of hypertrophic signaling for HOPX transcription

Differentiation into diverse cell lineages requires the orchestration of gene regulatory networks guiding diverse cell fate choices. Utilizing human pluripotent stem cells, we measured expression dynamics of 17,718 genes from 43,168 cells across five-time points over a thirty-day time-course of in vitro cardiac- directed differentiation. Unsupervised clustering and lineage prediction algorithms were used to map fate choices and transcriptional networks underlying cardiac differentiation. We leveraged this resource to identify strategies for controlling in vitro differentiation as it occurs in vivo. HOPX, a non-DNA binding homeodomain protein essential for heart development in vivo was identified as dysregulated in vitro derived cardiomyocytes. Utilizing genetic gain and loss of function approaches, we dissect the transcriptional complexity of the HOPX locus and identify the requirement of hypertrophic signaling for HOPX transcription in hPSC-derived cardiomyocytes. This work provides a single cell dissection of the transcriptional landscape of cardiac differentiation for broad applications of stem cells in cardiovascular biology.Differentiation into diverse cell lineages requires orchestration of gene regulatory networks guiding cell fate choices. Here, we present the dissection of cellular composition and gene networks from transcriptomic data of 43,168 cells across five discrete time points during cardiac-directed differentiation. We utilize unsupervised clustering and implement a lineage trajectory prediction algorithm that integrates transcription factor networks to predict cell fate progression of 15 subpopulations that correlate with germ layer and cardiovascular differentiation in vivo. These data reveal transcriptional networks underlying lineage derivation of mesoderm, definitive endoderm, vascular endothelium, cardiac precursors, and definitive cell types that comprise cardiomyocytes and a previously unrecognized cardiac outflow tract population. Single cell analysis of genetic regulators governing cardiac fate diversification identified the non-DNA binding homeodomain protein, HOPX, as functionally necessary for endothelial specification. Our findings further implicate dysregulation of HOPX during in vitro cardiac-directed differentiation underlying the molecular and physiological immaturity of stem cell-derived cardiomyocytes.

Detection of HPV E7 Transcription at Single-Cell Resolution in Epidermis

Persistent human papillomavirus (HPV) infection is responsible for at least 5% of human malignancies. Most HPV-associated cancers are initiated by the HPV16 genotype, as confirmed by detection of integrated HPV DNA in cells of oral and anogenital epithelial cancers. However, single-cell RNA sequencing may enable prediction of HPV involvement in carcinogenesis at other sites. We conducted single-cell RNA sequencing on keratinocytes from a mouse transgenic for the E7 gene of HPV16 and showed sensitive and specific detection of HPV16-E7 mRNA, predominantly in basal keratinocytes. We showed that increased E7 mRNA copy number per cell was associated with increased expression of E7 induced genes. This technique enhances detection of active viral transcription in solid tissue and may clarify possible linkage of HPV infection to development of squamous cell carcinoma.

Single-Cell Transcriptome Sequencing Of 18,787 Human Induced Pluripotent Stem Cells Identifies Differentially Primed Subpopulations

For pluripotent stem cells, transcriptional profiling is central to discovering the key genes and gene networks governing the undifferentiated state. However, the heterogeneity of cell states represented in pluripotent cultures have not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing (scRNA-seq) can be used to increase our understanding of how individual pluripotent cells function. Here, we present the scRNA-seq results of 18,787 individual WTC CRISPRi human induced pluripotent stem cells. Four subpopulations were distinguishable on the basis of their pluripotent state including: quiescent (48.3%), proliferative (47.8%), early-primed for differentiation (2.8%) and late-primed for differentiation (1.1%). We identified novel genes and pathways defining each of the subpopulations and developed a multigenic prediction model to accurately classify single cells into subpopulations. This study provides a benchmark single cell dataset that expands our understanding of the cellular complexity of pluripotency.Heterogeneity of cell states represented in pluripotent cultures have not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing can be used to identify how individual pluripotent cells function. Here, we present results from the analysis of single-cell RNA sequencing data from 18,787 individual WTC CRISPRi human induced pluripotent stem cells. We developed an unsupervised clustering method, and through this identified four subpopulations distinguishable on the basis of their pluripotent state including: a core pluripotent population (48.3%), proliferative (47.8%), early-primed for differentiation (2.8%) and late-primed for differentiation (1.1%). For each subpopulation we were able to identify the genes and pathways that define differences in pluripotent cell states. Our method identified four transcriptionally distinct predictor gene sets comprised of 165 unique genes that denote the specific pluripotency states; and using these sets, we developed a multigenic machine learning prediction method to accurately classify single cells into each of the subpopulations. Compared against a set of established pluripotency markers, our method increases prediction accuracy by 10%, specificity by 20%, and explains a substantially larger proportion of deviance (up to 3-fold) from the prediction model. Finally, we developed an innovative method to predict cells transitioning between subpopulations, and support our conclusions with results from two orthogonal pseudotime trajectory methods.