Multimodal spatiotemporal phenotyping of human retinal organoid development.

Organoids generated from human pluripotent stem cells provide experimental systems to study development and disease, but quantitative measurements across different spatial scales and molecular modalities are lacking. In this study, we generated multiplexed protein maps over a retinal organoid time course and primary adult human retinal tissue. We developed a toolkit to visualize progenitor and neuron location, the spatial arrangements of extracellular and subcellular components and global patterning in each organoid and primary tissue. In addition, we generated a single-cell transcriptome and chromatin accessibility timecourse dataset and inferred a gene regulatory network underlying organoid development. We integrated genomic data with spatially segmented nuclei into a multimodal atlas to explore organoid patterning and retinal ganglion cell (RGC) spatial neighborhoods, highlighting pathways involved in RGC cell death and showing that mosaic genetic perturbations in retinal organoids provide insight into cell fate regulation.

A gene expression atlas of embryonic neurogenesis in Drosophila reveals complex spatiotemporal regulation of lncRNAs.

Cell type specification during early nervous system development in Drosophila melanogaster requires precise regulation of gene expression in time and space. Resolving the programs driving neurogenesis has been a major challenge owing to the complexity and rapidity with which distinct cell populations arise. To resolve the cell type-specific gene expression dynamics in early nervous system development, we have sequenced the transcriptomes of purified neurogenic cell types across consecutive time points covering crucial events in neurogenesis. The resulting gene expression atlas comprises a detailed resource of global transcriptome dynamics that permits systematic analysis of how cells in the nervous system acquire distinct fates. We resolve known gene expression dynamics and uncover novel expression signatures for hundreds of genes among diverse neurogenic cell types, most of which remain unstudied. We also identified a set of conserved long noncoding RNAs (lncRNAs) that are regulated in a tissue-specific manner and exhibit spatiotemporal expression during neurogenesis with exquisite specificity. lncRNA expression is highly dynamic and demarcates specific subpopulations within neurogenic cell types. Our spatiotemporal transcriptome atlas provides a comprehensive resource for investigating the function of coding genes and noncoding RNAs during crucial stages of early neurogenesis.

The Drosophila embryo at single-cell transcriptome resolution.

By the onset of morphogenesis, Drosophila embryos consist of about 6000 cells that express distinct gene combinations. Here, we used single-cell sequencing of precisely staged embryos and devised DistMap, a computational mapping strategy to reconstruct the embryo and to predict spatial gene expression approaching single-cell resolution. We produced a virtual embryo with about 8000 expressed genes per cell. Our interactive Drosophila Virtual Expression eXplorer (DVEX) database generates three-dimensional virtual in situ hybridizations and computes gene expression gradients. We used DVEX to uncover patterned expression of transcription factors and long noncoding RNAs, as well as signaling pathway components. Spatial regulation of Hippo signaling during early embryogenesis suggests a mechanism for establishing asynchronous cell proliferation. Our approach is suitable to generate transcriptomic blueprints for other complex tissues.

Cell fixation and preservation for droplet-based single-cell transcriptomics.

BACKGROUND: Recent developments in droplet-based microfluidics allow the transcriptional profiling of thousands of individual cells in a quantitative, highly parallel and cost-effective way. A critical, often limiting step is the preparation of cells in an unperturbed state, not altered by stress or ageing. Other challenges are rare cells that need to be collected over several days or samples prepared at different times or locations. METHODS: Here, we used chemical fixation to address these problems. Methanol fixation allowed us to stabilise and preserve dissociated cells for weeks without compromising single-cell RNA sequencing data. RESULTS: By using mixtures of fixed, cultured human and mouse cells, we first showed that individual transcriptomes could be confidently assigned to one of the two species. Single-cell gene expression from live and fixed samples correlated well with bulk mRNA-seq data. We then applied methanol fixation to transcriptionally profile primary cells from dissociated, complex tissues. Low RNA content cells from Drosophila embryos, as well as mouse hindbrain and cerebellum cells prepared by fluorescence-activated cell sorting, were successfully analysed after fixation, storage and single-cell droplet RNA-seq. We were able to identify diverse cell populations, including neuronal subtypes. As an additional resource, we provide 'dropbead', an R package for exploratory data analysis, visualization and filtering of Drop-seq data. CONCLUSIONS: We expect that the availability of a simple cell fixation method will open up many new opportunities in diverse biological contexts to analyse transcriptional dynamics at single-cell resolution.

Evaluating eight newly identified susceptibility loci for nonsyndromic cleft lip with or without cleft palate in a Mesoamerican population.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is among the most frequently occurring congenital malformations worldwide. The number of genetic loci identified as being involved in NSCL/P etiology was recently increased by a large genome-wide meta-analysis of European and Asian samples. This meta-analysis confirmed all six previously recognized genetic susceptibility loci and identified six novel ones. METHODS: To investigate which of these 12 loci contribute to NSCL/P risk in an independent sample of distinct ethnicity, we performed a case-control association analysis in a sample of the Mesoamerican population. A total of 153 individuals with NSCL/P (cases) and 337 unaffected controls were included. Top single-nucleotide polymorphisms (SNPs) at 8 of the 12 loci (1p22.1, 1p36, 2p21, 3p11.1, 8q21.3, 13q31.1, 15q22, and 20q12) were analyzed using mass spectroscopy and restriction-length-fragment polymorphism analyses. In a previous study, we had analyzed the remaining four NSCL/P susceptibility regions (IRF6, 8q24, 10q25, and 17q22) in the same sample. RESULTS: Single-marker association analyses applying allelic, dominant, and recessive models revealed nominal significant associations for four of the eight loci, with two additional loci showing at least a trend of association in the hypothesized direction. CONCLUSION: In combination with results from our previous study using the same sample, our data suggest that the majority of the known NSCL/P susceptibility regions identified to date also confer risk for this malformation in the Mesoamerican population. Birth Defects Research (Part A) 100:43-47, 2014. © 2013 Wiley Periodicals, Inc.

Resequencing of VAX1 in patients with nonsyndromic cleft lip with or without cleft palate.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common of all congenital anomalies, and has a multifactorial etiology involving both environmental and genetic factors. Recent genome-wide association studies (GWAS) identified strong association between a locus on chromosome 10q25.3 and NSCL/P in European samples. One gene at 10q25.3, the ventral anterior homeobox 1 (VAX1) gene, is considered a strong candidate gene for craniofacial malformations. The purpose of the present study was to provide further evidence that VAX1 is the causal gene at the 10q25.3 locus through identification of an excess of rare mutations in patients with NSCL/P. METHODS: The 5'UTR, complete coding regions, and adjacent splice sites of the two known VAX1 isoforms were sequenced in 384 patients with NSCL/P and 384 controls of Central European descent. Observed variants were investigated with respect to familial cosegregation or de novo occurrence, and in silico analyses were performed to identify putative effects on the transcript or protein level. RESULTS: Eighteen single-base variants were found, 15 of them rare and previously unreported. In the long VAX1 isoform, predicted functionally relevant variants were observed more often in NSCL/P cases, although this difference was not significant (p = 0.17). Analysis of family members demonstrated incomplete cosegregation in most pedigrees. CONCLUSION: Our data do not support the hypothesis that highly penetrant rare variants in VAX1 are a cause of NSCL/P. To determine whether VAX1 is the causative gene at 10q25.3 further research, in particular into the biologic function of its long isoform, is warranted. Birth Defects Research (Part A), 2012.