The rapid improvements in single cell sequencing technologies and analyses methods afford greater scope for dissecting organoid cultures composed of multiple cell types and create an opportunity to interrogate these models to understand tissue biology, cellular behaviour and interactions. To this end, retinal organoids generated from human embryonic stem cells (hESCs) were analysed by single cell RNA-Sequencing at three time points of differentiation. Combinatorial data from all time points revealed the presence of nine clusters, five of which corresponded to key retinal cell types, namely retinal pigment epithelium (RPE), retinal ganglion cells (RGCs), cone and rod photoreceptors and Müller glia cells. The remaining four clusters expressed genes typical of mitotic cells, extracellular matrix (ECM) components and those involved in retinal homeostasis. The cell clustering analysis revealed the decreasing presence of mitotic cells and RGCs, formation of a distinct RPE cluster, the emergence of cone and rod photoreceptors from photoreceptor precursors and an increasing number of Müller Glia cells over time. The pseudotime analysis resembled the order of cell birth during retinal development, with the mitotic cluster commencing the trajectory and the large majority of Müller glia being the latest. Together, these data demonstrate the feasibility and potential of single cell RNA-Seq to dissect the inherent complexity of the organoids and the orderly birth of key retinal cell types. Overall design: A hESC (H9) cell line harbouring a CRX-GFP reporter was differentiated to retinal organoids 25. Samples were collected at 60, 90 and 200 days, dissociated, partitioned into single cells using the Fluidigm C1 Single-Cell mRNA-Seq HT IFC and processed for scRNA-Seq.
Deconstructing Retinal Organoids: Single Cell RNA-Seq Reveals the Cellular Components of Human Pluripotent Stem Cell-Derived Retina.
Cell line, Subject, Time
View SamplesDeath of photoreceptors and/or Retinal Pigment Epithelium (RPE) cells is a common cause of age related and inherited retinal dystrophies, thus their replenishment from renewable stem cell sources is a well sought therapeutic goal. Human pluripotent stem cells provide a useful cell source in view of their limitless self-renewal capacity and potential to differentiate into all key retinal cell types either in isolation or as part of three dimensional retinal organoids. Photoreceptor precursors have been isolated from differentiating human pluripotent stem cells either through application of cell surface markers or fluorescent reporter approaches and shown to share a transcriptional profile akin to foetal photoreceptors. In this study we investigated the transcriptional profile of CRX+ photoreceptor precursors derived from human embryonic stem cells (hESC) using single cell RNA sequencing and their engraftment capacity in an animal model of retinitis pigmentosa (C3H/rd1). Single cell RNA seq analysis revealed the presence of dominant cell cluster which displayed the hallmarks of early cone photoreceptor expression. When transplanted subretinally into the C3H/rd1 mice, the Crx positive cells settled next to the inner nuclear layer of host retina, matured into cone photoreceptors and made connections with the inner neurones of the host retina. Cellular transfer between the host retina and donor photoreceptors was investigated and shown to be minimal. Together our data provide valuable molecular insights into the transcriptional profile of human pluripotent stem cells derived CRX+ photoreceptor precursors and indicate their usefulness as a source of transplantable cone photoreceptors. Overall design: CRX-GFP human ESC line was differentiated to retinal organoids. At day 90 CRX+ and CRX- cells were purified by flow activated cell sorting and subjected to single cell RNA-seq. RNA-seq of bulk CRX+ and CRX- from the same experiment was carried out in parallel.
CRX Expression in Pluripotent Stem Cell-Derived Photoreceptors Marks a Transplantable Subpopulation of Early Cones.
Specimen part, Subject
View SamplesMitochondrial DNA (mtDNA) mutations cause inherited diseases and are implicated in the pathogenesis of common late-onset disorders, but it is not clear how they arise and propagate in the humans. Here we show that mtDNA mutations are present in primordial germ cells (PGCs) within healthy female human embryos. Close scrutiny revealed the signature of selection against non-synonymous variants in the protein-coding region, tRNA gene variants, and variants in specific regions of the non-coding D-loop. In isolated single PGCs we saw a profound reduction in the cellular mtDNA content, with discrete mitochondria containing ~5 mtDNA molecules during early germline development. Single cell deep mtDNA sequencing showed rare variants reaching higher heteroplasmy levels in later PGCs, consistent with the observed genetic bottleneck, and predicting >80% levels within isolated organelles. Genome-wide RNA-seq showed a progressive upregulation of genes involving mtDNA replication and transcription, linked to a transition from glycolytic to oxidative metabolism. The metabolic shift exposes deleterious mutations to selection at the organellar level during early germ cell development. In this way, the genetic bottleneck prevents the relentless accumulation of mtDNA mutations in the human population predicted by Muller's ratchet. Mutations escaping this mechanism will, however, show massive shifts in heteroplasmy levels within one human generation, explaining the extreme phenotypic variation seen in human pedigrees with inherited mtDNA disorders. Overall design: RNA-Seq and NGS analysis to investigate transcriptomes and mtDNA sequences of fetal hPGCs
Segregation of mitochondrial DNA heteroplasmy through a developmental genetic bottleneck in human embryos.
No sample metadata fields
View SamplesUsing a combination of cell sorting and microarray analysis, we identified almost 200 genes as having a high level of expression in the notochord.
Integrated microarray and ChIP analysis identifies multiple Foxa2 dependent target genes in the notochord.
Sex
View SamplesHigh-fat diets are associated with increased obesity and metabolic disease in mice and humans. Here we used analysis of variance (ANOVA) to scrutinize a microarray data set consisting of 10 inbred strains of mice from both sexes fed atherogenic high-fat and control chow diets. An overall F-test was applied to the 40 unique groups of strain-diet-sex to identify 15,288 genes with altered transcription. Bootstrapping k-means clustering separated these changes into four strain-dependent expression patterns, including two sex-related profiles and two diet-related profiles. Sex-induced effects correspond to secretion (males) or fat and energy metabolism (females), whereas diet-induced changes relate to neurological processes (chow) or immune response (high-fat). The full set of pairwise contrasts for differences between strains within sex (90 different statistical tests) uncovered 32,379 total changes. These differences were unevenly distributed across strains and between sexes, indicating that strain-specific responses to high-fat diet differ between sexes. Correlations between expression levels and 8 obesity-related traits identified 5,274 associations between transcript abundance and measured phenotypic endpoints. From this number, 2,678 genes are positively correlated with total cholesterol levels and associate with immune-related categories while 2,596 genes are negatively correlated with cholesterol and connect to cholesterol synthesis.
Practical applications of the bioinformatics toolbox for narrowing quantitative trait loci.
Sex
View SamplesDamage to and/or loss of sensory neurons can result in debilitating neuropathies that often have a dramatic impact on quality of life. The cellular mechanisms involved in the response of neurons and glia to such pathological insults are poorly understood. Investigation has shown that peripheral glia play critical roles in both the degenerative and regenerative processes that are involved in the responses to peripheral nerve damage. The vast majority of studies have focused primarily on myelinating Schwann cells], with the result that very little is known regarding how the non-myelinating glia that ensheath axons and neuronal somas respond to nerve damage. This is a significant knowledge gap, given that over 80% of cutaneous fibers are unmyelinated, that they transduce such important modalities as itch, pain, temperature, touch and pressure, and that they are affected in many prevalent peripheral neuropathies. It is the goal of this study to shed light on the genetic programs involved in the responses of non-myelinating glia roles to nerve degeneration. We utilized RNA-seq to identify genes that were differentially expressed in the larval head during the process of sensory neuron ablation and axon degeneration in both wild-type larvae and in larvae that do not have peripheral glia (foxd3 mutants) using a selective, conditional approach. Overall, the information regarding differential gene expression in these conditions will provide a basis for further investigation into the cellular processes that underlie pathophysiological responses of neurons and glia to sensory nerve damage. Overall design: mRNA levels were determined using biological triplicate samples from five sets of samples. Three sets from wild-type: control, 2 hrs of metronidazole treatment and 5 hrs of metronidazole treatment. And two sets from foxd3 mutants: control and 5hrs of metronidazole treatment.
Transcriptome Analysis of Chemically-Induced Sensory Neuron Ablation in Zebrafish.
No sample metadata fields
View SamplesBivalent histone domains have been proposed to contribute to pluripotency in embryonic stem cells, suggesting an epigenetic mechanism may regulate stem cell behavior in general. Here we compare histone modifications in two other stem cells derived from the blastocyst. We show that extraembryonic stem cells have little repressive lysine 27 trimethylation and few bivalent domains. Thus, bivalent domains are not a common mechanism for maintaining the undifferentiated state in blastocyst-derived stem cells and alternative mechanisms must mediate transcriptional repression in extraembryonic cells. We show that lysine 9 trimethylation, but not DNA methylation, is likely to fulfill this role. Intriguingly, although we do detect bivalent domains in pluripotent cells in the early mouse embryo, the epigenetic status of extraembryonic cells does not entirely reflect their in vitro stem cell counterparts. Therefore, differences in epigenetic regulation between lineage progenitors in vivo and in vitro may arise during selection for self-renewal in vitro.
Distinct histone modifications in stem cell lines and tissue lineages from the early mouse embryo.
Cell line
View SamplesVery little is known about how intervertebral disc (IVD) is formed or maintained. Members of the TGF- superfamily are secreted signaling proteins that regulate many aspects of development including cellular differentiation. We recently showed that deletion of Tgfbr2 in Col2a expressing tissue results in alterations in development of IVD annulus fibrosus. The results suggested TGF- has an important role in regulating development of the axial skeleton, however, the mechanistic basis of TGF- action in these specialized joints is not known. One of the hurdles to understanding development of IVD is a lack of known markers. To identify genes that are enriched in the developing IVD and to begin to understand the mechanism of TGF- action in IVD development, we undertook a global analysis of gene expression comparing gene expression profiles in developing vertebrae and IVD. We also compared expression profiles in tissues from wild type and Tgfbr2 mutant mice. Lists of IVD and vertebrae enriched genes were generated. Expression patterns for several genes were verified either through in situ hybridization or literature/ database searches resulting in a list of genes that can be used as markers of IVD. Cluster analysis using genes listed under the Gene Ontology terms multicellular organism development and pattern specification indicated that mutant IVD more closely resembled vertebrae than wild type IVD. We propose TGF- has two functions in IVD development: 1) to prevent chondrocyte differentiation in the presumptive IVD and 2) to promote differentiation of annulus fibrosus from sclerotome. We have identified genes that are enriched in the IVD and regulated by TGF- that warrant further investigation as regulators of IVD development.
Molecular profiling of the developing mouse axial skeleton: a role for Tgfbr2 in the development of the intervertebral disc.
Specimen part
View SamplesVery little is known about how intervertebral disc (IVD) is formed or maintained. Members of the TGF- superfamily are secreted signaling proteins that regulate many aspects of development including cellular differentiation. We recently showed that deletion of Tgfbr2 in Col2a expressing tissue results in alterations in development of IVD annulus fibrosus. The results suggested TGF- has an important role in regulating development of the axial skeleton, however, the mechanistic basis of TGF- action in these specialized joints is not known. To understand the mechanism of TGF- action in IVD development, we undertook a global analysis of gene expression comparing gene expression profiles in sclerotome cultures treated with TGF- or BMP4. As expected, treatment with BMP4 resulted in up-regulation of cartilage marker genes including Acan, Sox 5, Sox6, and Sox9. In contrast, treatment with TGF-1 did not regulate expression of cartilage markers but instead resulted in up-regulation of many IVD markers including Fmod and Adamtsl2. We propose TGF- has two functions in IVD development: 1) to prevent chondrocyte differentiation in the presumptive IVD and 2) to promote differentiation of annulus fibrosus from sclerotome. We have identified genes that are enriched in the IVD and regulated by TGF- that warrant further investigation as regulators of IVD development.
Molecular profiling of the developing mouse axial skeleton: a role for Tgfbr2 in the development of the intervertebral disc.
No sample metadata fields
View SamplesThe role of microRNAs (miRNA) in first cell fate choice of the preimplantation mouse embryo remains unresolved, as gene expression and knockout data are conflicting. This cell fate choice generates the extraembryonic lineage of the trophoblast and the embryonic lineage of the epiblast (inner cell mass). The trophoblast differentiates into polar and mural cells, where polar cells contribute to placental development and mural cells to the implantation process and Reicherts membrane. The inner cell mass further differentiates into the epiblast and primitive endoderm. We used stem cell lines that can be derived from the trophoblast and epiblast lineages to address the role of miRNAs in early lineage cell fate specification and determination. Using embryonic stem cells (ESC) and trophoblast stem cells (TSC) as starting and ending states of cell development we identified a network of TSC expressed miRNAs that are enriched in ESC targets mRNA. We used constitutive and inducible expression of TSC enriched miRNAs in ESC and show that they can drive cell morphology and gene expression profiles similar to trophoblast. Additionally we show that this process required HDAC2 inhibition and is miRNA specific, as cardiac specific miR-1 could not induce trophoblast under these conditions. In contrast to embryo derived and Cdx2 induced trophoblast cells, miRNAs promote a mural TE like cell phenotype. Transplantation into preimplantation mouse embryos showed that miRNA-induced trophoblast preferentially contributes to the mural trophoblast in both the blastocyst and the Reicherts membrane. Our data support a role for miRNAs and HDACs in the specification of the trophoblast and potentially the polar and mural cell types.
Overexpression of Trophoblast Stem Cell-Enriched MicroRNAs Promotes Trophoblast Fate in Embryonic Stem Cells.
No sample metadata fields
View Samples