By the onset of morphogenesis, 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 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.

Drosophila is a premier model system for understanding the molecular mechanisms of development. By the onset of morphogenesis, ~6000 cells express distinct gene combinations according to embryonic position. Despite extensive mRNA in situ screens, combinatorial gene expression within individual cells is largely unknown. Therefore, it is difficult to comprehensively identify the coding and non-coding transcripts that drive patterning and to decipher the molecular basis of cellular identity. Here, we single-cell sequence precisely staged embryos, measuring >3100 genes per cell. We produce a 'transcriptomic blueprint' of development -a virtual embryo where 3D locations of sequenced cells are confidently identified.Our "Drosophila-Virtual-Expression-eXplorer" performs virtual in situ hybridizations and computes expression gradients. Using DVEX, we predict spatial expression and discover patterned lncRNAs. DEVX is sensitive enough to detect subtle evolutionary changes in expression patterns between Drosophila species. We believe DVEX is a prototype for powerful single cell studies in complex tissues.

BackgroundRecent 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.MethodsHere, 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.ResultsBy 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.ConclusionsWe 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12915-017-0383-5) contains supplementary material, which is available to authorized users.

Background: Recent developments in droplet-based microfluidics allow the transcriptional profiling of

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.