Heterogeneous combinatorial expression of Hoxd genes in single cells during limb development

BST2 promotes gastric cancer metastasis under the regulation of HOXD9 and PABPC1

Gastric cancer (GC) constitutes substantial cancer mortality worldwide. Several cancer types aberrantly express bone marrow stromal cell antigen 2 (BST2), yet its functional and underlying mechanisms in GC progression remain unknown. In our study, RNA sequencing data revealed that BST2 was transcriptionally activated by homeobox D9 (HOXD9). BST2 was significantly upregulated in GC tissues and promoted epithelial-mesenchymal transition and metastasis of GC. BST2 knockdown reversed HOXD9's oncogenic effect on GC metastasis. Moreover, BST2 messenger RNA stability could be enhanced by poly(A) binding protein cytoplasmic 1 (PABPC1) through the interaction between BST2 3'-UTR and PABPC1 in GC cells. PABPC1 promoted GC metastasis, which BST2 silencing attenuated in vitro and in vivo. In addition, positive correlations among HOXD9, BST2, and PABPC1 were established in clinical samples. Taken together, increased expression of BST2 induced by HOXD9 synergizing with PABPC1 promoted GC cell migration and invasion capacity.

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

A single-cell census of mouse limb development identifies complex spatiotemporal dynamics of skeleton formation

Limb development has long served as a model system for coordinated spatial patterning of progenitor cells. Here, we identify a population of naive limb progenitors and show that they differentiate progressively to form the skeleton in a complex, non-consecutive, three-dimensional pattern. Single-cell RNA sequencing of the developing mouse forelimb identified three progenitor states: naive, proximal, and autopodial, as well as Msx1 as a marker for the naive progenitors. In vivo lineage tracing confirmed this role and localized the naive progenitors to the outer margin of the limb, along the anterior-posterior axis. Sequential pulse-chase experiments showed that the progressive transition of Msx1⁺ naive progenitors into proximal and autopodial progenitors coincides with their differentiation to Sox9⁺ chondroprogenitors, which occurs along all the forming skeletal segments. Indeed, tracking the spatiotemporal sequence of differentiation showed that the skeleton forms progressively in a complex pattern. These findings suggest an alternative model for limb skeleton development.

Theories, laws, and models in evo-devo

Evolutionary developmental biology (evo‐devo) is the study of the evolution of developmental mechanisms. Here, I review some of the theories, models, and laws in evo‐devo, past and present. Nineteenth‐century evo‐devo was dominated by recapitulation theory and archetypes. It also gave us germ layer theory, the vertebral theory of the skull, floral organs as modified leaves, and the “inverted invertebrate” theory, among others. Newer theories and models include the frameshift theory, the genetic toolkit for development, the ABC model of flower development, the developmental hourglass, the zootype, Urbilateria, and the hox code. Some of these new theories show the influence of archetypes and recapitulation. Interestingly, recent studies support the old “primordial leaf,” “inverted invertebrate,” and “segmented head” theories. Furthermore, von Baer's first three laws may now need to be rehabilitated, and the hourglass model modified, in view of what Abzhanov has pointed out about the maternal‐zygotic transition. There are many supposed “laws” of evo‐devo but I argue that these are merely generalizations about trends in particular lineages. I argue that the “body plan” is an archetype, and is often used in such a way that it lacks any scientific meaning. Looking to the future, one challenge for evo‐devo will be to develop new theories and models to accommodate the wealth of new data from high‐throughput sequencing, including single‐cell sequencing. One step in this direction is the use of sophisticated in silico analyses, as in the “transcriptomic hourglass” models.

Laws and other universal concepts, past and present, are reviewed.

I show that many concepts focus on conserved aspects of development.

The puzzle remains as to how embryonic phenotype, natural selection, and developmental mechanisms can be aligned to give an integrated view of evolution and development.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Scalable control of developmental timetables by epigenetic switching networks

During development, progenitor cells follow timetables for differentiation that span many cell generations. These developmental timetables are robustly encoded by the embryo, yet scalably adjustable by evolution, facilitating variation in organism size and form. Epigenetic switches, involving rate-limiting activation steps at regulatory gene loci, control gene activation timing in diverse contexts, and could profoundly impact the dynamics of gene regulatory networks controlling developmental lineage specification. Here, we develop a mathematical framework to model regulatory networks with genes controlled by epigenetic switches. Using this framework, we show that such epigenetic switching networks uphold developmental timetables that robustly span many cell generations, and enable the generation of differentiated cells in precisely defined numbers and fractions. Changes to epigenetic switching networks can readily alter the timing of developmental events within a timetable, or alter the overall speed at which timetables unfold, enabling scalable control over differentiated population sizes. With their robust, yet flexibly adjustable nature, epigenetic switching networks could represent central targets on which evolution acts to manufacture diversity in organism size and form.

Single-cell transcriptomics of the Drosophila wing disc reveals instructive epithelium-to-myoblast interactions

In both vertebrates and invertebrates, generating a functional appendage requires interactions between ectoderm-derived epithelia and mesoderm-derived cells. To investigate such interactions, we used single-cell transcriptomics to generate a temporal cell atlas of the Drosophila wing disc from two developmental time points. Using these data, we visualized gene expression using a multilayered model of the wing disc and cataloged ligand–receptor pairs that could mediate signaling between epithelial cells and adult muscle precursors (AMPs). We found that localized expression of the fibroblast growth factor ligands, Thisbe and Pyramus, in the disc epithelium regulates the number and location of the AMPs. In addition, Hedgehog ligand from the epithelium activates a specific transcriptional program within adjacent AMP cells, defined by AMP-specific targets Neurotactin and midline, that is critical for proper formation of direct flight muscles. More generally, our annotated temporal cell atlas provides an organ-wide view of potential cell–cell interactions between epithelial and myogenic cells.

Single-cell transcriptome profiling of an adult human cell atlas of 15 major organs

BACKGROUND

As core units of organ tissues, cells of various types play their harmonious rhythms to maintain the homeostasis of the human body. It is essential to identify the characteristics of cells in human organs and their regulatory networks for understanding the biological mechanisms related to health and disease. However, a systematic and comprehensive single-cell transcriptional profile across multiple organs of a normal human adult is missing.

RESULTS

We perform single-cell transcriptomes of 84,363 cells derived from 15 tissue organs of one adult donor and generate an adult human cell atlas. The adult human cell atlas depicts 252 subtypes of cells, including major cell types such as T, B, myeloid, epithelial, and stromal cells, as well as novel COCH+ fibroblasts and FibSmo cells, each of which is distinguished by multiple marker genes and transcriptional profiles. These collectively contribute to the heterogeneity of major human organs. Moreover, T cell and B cell receptor repertoire comparisons and trajectory analyses reveal direct clonal sharing of T and B cells with various developmental states among different tissues. Furthermore, novel cell markers, transcription factors, and ligand-receptor pairs are identified with potential functional regulations in maintaining the homeostasis of human cells among tissues.

CONCLUSIONS

The adult human cell atlas reveals the inter- and intra-organ heterogeneity of cell characteristics and provides a useful resource in uncovering key events during the development of human diseases in the context of the heterogeneity of cells and organs.

Atlas of Musculoskeletal Stem Cells with the Soft and Hard Tissue Differentiation Architecture

Although being of utmost importance for human health and mobility, stem cell identity and hierarchical organization of musculoskeletal progenitors remain largely unexplored. Here, cells from E10.5, E12.5, and E15.5 murine limbs are analyzed by high throughput single-cell RNA sequencing to illustrate the cellular architecture during limb development. Single-cell transcriptional profiling demonstrates the identity and differentiation architecture of musculoskeletal stem cells (MSSC), soft and hard tissue progenitors through expression pattern of musculoskeletal markers (scleraxis [Scx], Hoxd13, Sox9, and Col1a1). This is confirmed by genetic in vivo lineage tracing. Moreover, single-cell analyses of Scx knockout mice tissues illustrates that Scx regulates MSSC self-renewal and proliferation potential. A high-throughput and low-cost multi-tissues RNA sequencing strategy further provides evidence that musculoskeletal system tissues, including muscle, bone, meniscus, and cartilage, are all abnormally developed in Scx knockout mice. These results establish the presence of an indispensable limb Scx+Hoxd13+ MSSC population and their differentiation into soft tissue progenitors (Scx+Col1a1+) and hard tissue progenitors (Scx+Sox9+). Collectively, this study paves the way for systematically decoding the complex molecular mechanisms and cellular programs of musculoskeletal tissues morphogenesis in limb development and regeneration.

Geminin is required for Hox gene regulation to pattern the developing limb

Development of the complex structure of the vertebrate limb requires carefully orchestrated interactions between multiple regulatory pathways and proteins. Among these, precise regulation of 5' Hox transcription factor expression is essential for proper limb bud patterning and elaboration of distinct limb skeletal elements. Here, we identified Geminin (Gmnn) as a novel regulator of this process. A conditional model of Gmnn deficiency resulted in loss or severe reduction of forelimb skeletal elements, while both the forelimb autopod and hindlimb were unaffected. 5' Hox gene expression expanded into more proximal and anterior regions of the embryonic forelimb buds in this Gmnn-deficient model. A second conditional model of Gmnn deficiency instead caused a similar but less severe reduction of hindlimb skeletal elements and hindlimb polydactyly, while not affecting the forelimb. An ectopic posterior SHH signaling center was evident in the anterior hindlimb bud of Gmnn-deficient embryos in this model. This center ectopically expressed Hoxd13, the HOXD13 target Shh, and the SHH target Ptch1, while these mutant hindlimb buds also had reduced levels of the cleaved, repressor form of GLI3, a SHH pathway antagonist. Together, this work delineates a new role for Gmnn in modulating Hox expression to pattern the vertebrate limb.

Modeling Hypoxia-Induced Neuropathies Using a Fast and Scalable Human Motor Neuron Differentiation System

Human motor neuron (MN) diseases encompass a spectrum of disorders. A critical barrier to dissecting disease mechanisms is the lack of appropriate human MN models. Here, we describe a scalable, suspension-based differentiation system to generate functional human MN diseases in 3 weeks. Using this model, we translated recent findings that mRNA mis-localization plays a role in disease development to the human context by establishing a membrane-based system that allows efficient fractionation of MN cell soma and neurites. In response to hypoxia, used to mimic diabetic neuropathies, MNs upregulated mitochondrial transcripts in neurites; however, mitochondria were decreased. These data suggest that hypoxia may disrupt translation of mitochondrial mRNA, potentially leading to neurite damage and development of neuropathies. We report the development of a novel human MN model system to investigate mechanisms of disease affecting soma and/or neurites that facilitates the rapid generation and testing of patient-specific MN diseases.

Accelerated Evolution of Limb-Related Gene Hoxd11 in the Common Ancestor of Cetaceans and Ruminants (Cetruminantia)

Reduced numbers of carpal and tarsal bones (wrist and ankle joints) are extensively observed in the clade of Cetacea and Ruminantia (Cetruminantia). Homebox D11 (Hoxd11) is one of the important genes required for limb development in mammals. Mutations in Hoxd11 can lead to defects in particular bones of limbs, including carpus and tarsus. To test whether evolutionary changes in Hoxd11 underlie the loss of these bones in Cetruminantia, we sequenced and analyzed Hoxd11 coding sequences and compared them with other 5′ HoxA and HoxD genes in a taxonomic coverage of Cetacea, Ruminantia and other mammalian relatives. Statistical tests on the Hoxd11 sequences found an accelerated evolution in the common ancestor of cetaceans and ruminants, which coincided with the reduction of carpal and tarsal bones in this clade. Five amino acid substitutions (G222S, G227A, G229S, A240T and G261V) and one amino acid deletion (G254Del) occurred in this lineage. In contrast, other 5′ HoxA and HoxD genes do not show this same evolutionary pattern, but instead display a highly conserved pattern of evolution in this lineage. Accelerated evolution of Hoxd11, but not other 5′ HoxA and HoxD genes, is probably related to the reduction of the carpal and tarsal bones in Cetruminantia. Moreover, we found two amino acid substitutions (G110S and D223N) in Hoxd11 that are unique to the lineage of Cetacea, which coincided with hindlimb loss in the common ancestor of cetaceans. Our results give molecular evidence of Hoxd11 adaptive evolution in cetaceans and ruminants, which could be correlated with limb morphological adaptation.

Some Questions and Answers About the Role of Hox Temporal Collinearity in Vertebrate Axial Patterning

The vertebrate anterior-posterior (A-P = craniocaudal) axis is evidently made by a timing mechanism. Evidence has accumulated that tentatively identifies the A-P timer as being or involving Hox temporal collinearity (TC). Here, I focus on the two current competing models based on this premise. Common features and points of dissent are examined and a common model is distilled from what remains. This is an attempt to make sense of the literature.

Single-cell transcriptional logic of cell-fate specification and axon guidance in early-born retinal neurons

Retinal ganglion cells (RGCs), cone photoreceptors (cones), horizontal cells and amacrine cells are the first classes of neurons produced in the retina. However, an important question is how this diversity of cell states is transcriptionally produced. Here, we profiled 6067 single retinal cells to provide a comprehensive transcriptomic atlas showing the diversity of the early developing mouse retina. RNA velocities unveiled the dynamics of cell cycle coordination of early retinogenesis and define the transcriptional sequences at work during the hierarchical production of early cell-fate specification. We show that RGC maturation follows six waves of gene expression, with older-generated RGCs transcribing increasing amounts of guidance cues for young peripheral RGC axons that express the matching receptors. Spatial transcriptionally deduced features in subpopulations of RGCs allowed us to define novel molecular markers that are spatially restricted. Finally, the isolation of such a spatially restricted population, ipsilateral RGCs, allowed us to identify their molecular identity at the time they execute axon guidance decisions. Together, these data represent a valuable resource shedding light on transcription factor sequences and guidance cue dynamics during mouse retinal development.

Anatomic Origin of Osteochondrogenic Progenitors Impacts Sensitivity to EWS-FLI1-Induced Transformation

Ewing sarcomas predominantly arise in pelvic and stylopod bones (i.e., femur and humerus), likely as a consequence of EWS-FLI1 oncogene-induced transformation of mesenchymal stem/progenitor cells (MSCs). MSCs located in the embryonic superficial zone cells (eSZ) of limbs express anatomically distinct posterior Hox genes. Significantly, high expression of posterior HOXD genes, especially HOXD13, is a hallmark of Ewing sarcoma. These data drove our hypothesis that Hox genes in posterior skeleton MSCs contribute to Ewing sarcoma tumorigenesis. We isolated eSZ cells from stylopod and zeugopod (i.e., tibia/fibula, radius/ulna) bones, from wild-type and Hoxd13 mutant embryos, and tested the impact of EWS-FLI1 transduction on cell proliferation, gene expression, and tumorigenicity. Our data demonstrate that both stylopod and zeugopod eSZ cells tolerate EWS-FLI1 but that stylopod eSZ cells are relatively more susceptible, demonstrating changes in proliferation and gene expression consistent with initiation of malignant transformation. Significantly, loss of Hoxd13 had no impact, showing that it is dispensable for the initiation of EWS-FLI1-induced transformation in mouse MSCs. These findings show that MSCs from anatomically distinct sites are differentially susceptible to EWS-FLI1-induced transformation, supporting the premise that the dominant presentation of Ewing sarcoma in pelvic and stylopod bones is attributable to anatomically-defined differences in MSCs.