17Background. Through precise implementation of distinct cell type specification programs, 18 differentially regulated in both space and time, complex patterns emerge during 19 organogenesis. Thanks to its easy experimental accessibility, the developing chicken limb has 20 long served as a paradigm to study vertebrate pattern formation. Through decades' worth of 21 research, we now have a firm grasp on the molecular mechanisms driving limb formation at 22 the tissue-level. However, to elucidate the dynamic interplay between transcriptional cell type 23 specification programs and pattern formation at its relevant cellular scale, we lack 24 appropriately resolved molecular data at the genome-wide level. Here, making use of droplet-25 based single-cell RNA-sequencing, we catalogue the developmental emergence of distinct 26 tissue types and their transcriptome dynamics in the distal chicken limb, the so-called 27 autopod, at cellular resolution.
28Results. Using single-cell RNA-sequencing technology, we sequenced a total of 17,628 cells 29 coming from three key developmental stages of chicken autopod patterning. Overall, we 30 identified 23 cell populations with distinct transcriptional profiles. Amongst them were small, 31 albeit essential populations like the apical ectodermal ridge, demonstrating the ability to 32 detect even rare cell types. Moreover, we uncovered the existence of molecularly distinct sub-33 populations within previously defined compartments of the developing limb, some of which 34 have important signaling functions during autopod pattern formation. Finally, we inferred 35 gene co-expression modules that coincide with distinct tissue types across developmental 36 time, and used them to track patterning-relevant cell populations of the forming digits.37Conclusions. We provide a comprehensive functional genomics resource to study the 38 molecular effectors of chicken limb patterning at cellular resolution. Our single-cell 39 transcriptomic atlas captures all major cell populations of the developing autopod, and 40 highlights the transcriptional complexity in many of its components. Finally, integrating our 41 data-set with other single-cell transcriptomics resources will enable researchers to assess 42 molecular similarities in orthologous cell types across the major tetrapod clades, and provide 43 an extensive candidate gene list to functionally test cell-type-specific drivers of limb 44 morphological diversification. 45 46 Background 50 Embryonic pattern formation relies on the tight coordination of numerous developmental 51 processes, across multiple scales of complexity. From seemingly homogenous progenitor 52 populations, different cell types get specified and arranged in intricate patterns, to give rise to 53 functional tissues and organs. As progenitors mostly share a common genome, this 54 phenotypic specialization relies on the precise execution of distinct gene regulatory networks, 55 to enable cell type specification and ensuing pattern formation [1-3]. Slight deviations in 56 these p...