Single cell RNA-sequencing and RNA-tomography of the avian embryo extending body axis.

Introduction: Vertebrate body axis formation initiates during gastrulation and continues within the tail bud at the posterior end of the embryo. Major structures in the trunk are paired somites, which generate the musculoskeletal system, the spinal cord-forming part of the central nervous system, and the notochord, with important patterning functions. The specification of these different cell lineages by key signalling pathways and transcription factors is essential, however, a global map of cell types and expressed genes in the avian trunk is missing. Methods: Here we use high-throughput sequencing approaches to generate a molecular map of the emerging trunk and tailbud in the chick embryo. Results and Discussion: Single cell RNA-sequencing (scRNA-seq) identifies discrete cell lineages including somites, neural tube, neural crest, lateral plate mesoderm, ectoderm, endothelial and blood progenitors. In addition, RNA-seq of sequential tissue sections (RNA-tomography) provides a spatially resolved, genome-wide expression dataset for the avian tailbud and emerging body, comparable to other model systems. Combining the single cell and RNA-tomography datasets, we identify spatially restricted genes, focusing on somites and early myoblasts. Thus, this high-resolution transcriptome map incorporating cell types in the embryonic trunk can expose molecular pathways involved in body axis development.

When Bigger Is Better: 3D RNA Profiling of the Developing Head in the Catshark Scyliorhinus canicula.

We report the adaptation of RNA tomography, a technique allowing spatially resolved, genome-wide expression profiling, to a species occupying a key phylogenetic position in gnathostomes, the catshark Scyliorhinus canicula. We focused analysis on head explants at an embryonic stage, shortly following neural tube closure and of interest for a number of developmental processes, including early brain patterning, placode specification or the establishment of epithalamic asymmetry. As described in the zebrafish, we have sequenced RNAs extracted from serial sections along transverse, horizontal and sagittal planes, mapped the data onto a gene reference taking advantage of the high continuity genome recently released in the catshark, and projected read counts onto a digital model of the head obtained by confocal microscopy. This results in the generation of a genome-wide 3D atlas, containing expression data for most protein-coding genes in a digital model of the embryonic head. The digital profiles obtained for candidate forebrain regional markers along antero-posterior, dorso-ventral and left-right axes reproduce those obtained by in situ hybridization (ISH), with expected relative organizations. We also use spatial autocorrelation and correlation as measures to analyze these data and show that they provide adequate statistical tools to extract novel expression information from the model. These data and tools allow exhaustive searches of genes exhibiting any predefined expression characteristic, such a restriction to a territory of interest, thus providing a reference for comparative analyses across gnathostomes. This methodology appears best suited to species endowed with large embryo or organ sizes and opens novel perspectives to a wide range of evo-devo model organisms, traditionally counter-selected on size criterion.