Single-cell mRNA profiling reveals heterogeneous combinatorial expression of Hoxd genes during limb development

A global analysis of gene expression during development reveals specific transcription patterns associated with the emergence of various cell types, tissues and organs. These heterogeneous patterns are instrumental to ensure the proper formation of the different parts of our body, as shown by the phenotypic effects generated by functional genetic approaches. However, variations at the cellular level can be observed within each structure or organ. In the developing mammalian limbs, expression of Hoxd genes is differentially controlled in space and time in cells that will pattern the digits and the arms. Here we analyze single-cell transcriptomes of limb bud cells and show that Hox genes are expressed in specific combinations that match particular cell types. In the presumptive digits, we find that the expression of Hoxd gene is unbalanced, despite their common genomic proximity to known global enhancers, often expressing only a subset of the five genes transcribed in these cells. We also report that combinatorial expression follows a pseudo-time sequence, suggesting that a progression in combinatorial expression may be associated with cellular diversity in developing digits.