Transcriptional Trajectories in Mouse Limb Buds Reveal the Transition from Anterior-Posterior to Proximal-Distal Patterning at Early Limb Bud Stage

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Rapid and Quantitative Functional Interrogation of Human Enhancer Variant Activity in Live Mice

Functional analysis of non-coding variants associated with human congenital disorders remains challenging due to the lack of efficient in vivo models. Here we introduce dual-enSERT, a robust Cas9-based two-color fluorescent reporter system which enables rapid, quantitative comparison of enhancer allele activities in live mice of any genetic background. We use this new technology to examine and measure the gain- and loss-of-function effects of enhancer variants linked to limb polydactyly, autism, and craniofacial malformation. By combining dual-enSERT with single-cell transcriptomics, we characterize variant enhancer alleles at cellular resolution, thereby implicating candidate molecular pathways in pathogenic enhancer misregulation. We further show that independent, polydactyly-linked enhancer variants lead to ectopic expression in the same cell populations, indicating shared genetic mechanisms underlying non-coding variant pathogenesis. Finally, we streamline dual-enSERT for analysis in F0 animals by placing both reporters on the same transgene separated by a synthetic insulator. Dual-enSERT allows researchers to go from identifying candidate enhancer variants to analysis of comparative enhancer activity in live embryos in under two weeks.

Multi-species atlas resolves an axolotl limb development and regeneration paradox

Humans and other tetrapods are considered to require apical-ectodermal-ridge (AER) cells for limb development, and AER-like cells are suggested to be re-formed to initiate limb regeneration. Paradoxically, the presence of AER in the axolotl, a primary model organism for regeneration, remains controversial. Here, by leveraging a single-cell transcriptomics-based multi-species atlas, composed of axolotl, human, mouse, chicken, and frog cells, we first establish that axolotls contain cells with AER characteristics. Further analyses and spatial transcriptomics reveal that axolotl limbs do not fully re-form AER cells during regeneration. Moreover, the axolotl mesoderm displays part of the AER machinery, revealing a program for limb (re)growth. These results clarify the debate about the axolotl AER and the extent to which the limb developmental program is recapitulated during regeneration.

Sox, Fox, and Lmx1b binding sites differentially regulate a Gdf5-Associated regulatory region during elbow development

Introduction: The articulating ends of limb bones have precise morphology and asymmetry that ensures proper joint function. Growth differentiation factor 5 (Gdf5) is a secreted morphogen involved in cartilage and bone development that contributes to the architecture of developing joints. Dysregulation of Gdf5 results in joint dysmorphogenesis often leading to progressive joint degeneration or osteoarthritis (OA). The transcription factors and cis-regulatory modules (CRMs) that regulate Gdf5 expression are not well characterized. We previously identified a Gdf5-associated regulatory region (GARR) that contains predicted binding sites for Lmx1b, Osr2, Fox, and the Sox transcription factors. These transcription factors are recognized factors involved in joint morphogenesis and skeletal development. Methods: We used in situ hybridization to Gdf5, Col2A1, and the transcription factors of interest in developing chicken limbs to determine potential overlap in expression. We further analyzed scRNA-seq data derived from limbs and knees in published mouse and chicken datasets, identifying cells with coexpression of Gdf5 and the transcription factors of interest. We also performed site-directed mutatgenesis of the predicted transcription factor binding sites in a GARR-reporter construct and determined any change in activity using targeted regional electroporation (TREP) in micromass and embryonic chicken wing bioassays. Results: Gdf5 expression overlapped the expression of these transcription factors during joint development both by in situ hybridization (ISH) and scRNA-seq analyses. Within the GARR CRM, mutation of two binding sites common to Fox and Sox transcripstion factors reduced enhancer activity to background levels in micromass cultures and in ovo embryonic chicken wing bioassays, whereas mutation of two Sox-only binding sites caused a significant increase in activity. These results indicate that the Fox/Sox binding sites are required for activity, while the Sox-only sites are involved in repression of activity. Mutation of Lmx1b binding sites in GARR caused an overall reduction in enhancer activity in vitro and a dorsal reduction in ovo. Despite a recognized role for Osr2 in joint development, disruption of the predicted Osr2 site did not alter GARR activity. Conclusion: Taken together, our data indicates that GARR integrates positive, repressive, and asymmetrical inputs to fine-tune the expression of Gdf5 during elbow joint development.

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.

Single-cell response to Wnt signaling activation reveals uncoupling of Wnt target gene expression

Wnt signaling drives nuclear translocation of β-catenin and its subsequent association with the DNA-bound TCF/LEF transcription factors, which dictate target gene specificity by recognizing Wnt responsive elements across the genome. β-Catenin target genes are therefore thought to be collectively activated upon Wnt pathway stimulation. However, this appears in contrast with the non-overlapping patterns of Wnt target gene expression in several contexts, including early mammalian embryogenesis. Here we followed Wnt target gene expression in human embryonic stem cells after Wnt pathway stimulation at a single-cell resolution. Cells changed gene expression program over time consistent with three key developmental events: i) loss of pluripotency, ii) induction of Wnt target genes, and iii) mesoderm specification. Contrary to our expectation, not all cells displayed equal amplitude of Wnt target gene activation; rather, they distributed in a continuum from strong to weak responders when ranked based on the expression of the target AXIN2. Moreover, high AXIN2 did not always correspond to elevated expression of other Wnt targets, which were activated in different proportions in individual cells. The uncoupling of Wnt target gene expression was also identified in single cell transcriptomics profiling of other Wnt-responding cell types, including HEK293T, murine developing forelimbs, and human colorectal cancer. Our finding underlines the necessity to identify additional mechanisms that explain the heterogeneity of the Wnt/β-catenin-mediated transcriptional outputs in single cells.

Digits in a dish: An in vitro system to assess the molecular genetics of hand/foot development at single-cell resolution

In vitro models allow for the study of developmental processes outside of the embryo. To gain access to the cells mediating digit and joint development, we identified a unique property of undifferentiated mesenchyme isolated from the distal early autopod to autonomously re-assemble forming multiple autopod structures including: digits, interdigital tissues, joints, muscles and tendons. Single-cell transcriptomic analysis of these developing structures revealed distinct cell clusters that express canonical markers of distal limb development including: Col2a1, Col10a1, and Sp7 (phalanx formation), Thbs2 and Col1a1 (perichondrium), Gdf5, Wnt5a, and Jun (joint interzone), Aldh1a2 and Msx1 (interdigital tissues), Myod1 (muscle progenitors), Prg4 (articular perichondrium/articular cartilage), and Scx and Tnmd (tenocytes/tendons). Analysis of the gene expression patterns for these signature genes indicates that developmental timing and tissue-specific localization were also recapitulated in a manner similar to the initiation and maturation of the developing murine autopod. Finally, the in vitro digit system also recapitulates congenital malformations associated with genetic mutations as in vitro cultures of Hoxa13 mutant mesenchyme produced defects present in Hoxa13 mutant autopods including digit fusions, reduced phalangeal segment numbers, and poor mesenchymal condensation. These findings demonstrate the robustness of the in vitro digit system to recapitulate digit and joint development. As an in vitro model of murine digit and joint development, this innovative system will provide access to the developing limb tissues facilitating studies to discern how digit and articular joint formation is initiated and how undifferentiated mesenchyme is patterned to establish individual digit morphologies. The in vitro digit system also provides a platform to rapidly evaluate treatments aimed at stimulating the repair or regeneration of mammalian digits impacted by congenital malformation, injury, or disease.

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Identification of Transcription Factor Networks during Mouse Hindlimb Development

Mammalian hindlimb development involves a variety of cells and the regulation of spatiotemporal molecular events, but regulatory networks of transcription factors contributing to hindlimb morphogenesis are not well understood. Here, we identified transcription factor networks during mouse hindlimb morphology establishment through transcriptome analysis. We used four stages of embryonic hindlimb transcription profiles acquired from the Gene Expression Omnibus database (GSE30138), including E10.5, E11.5, E12.5 and E13.5, to construct a gene network using Weighted Gene Co-expression Network Analysis (WGCNA), and defined seven stage-associated modules. After filtering 7625 hub genes, we further prioritized 555 transcription factors with AnimalTFDB3.0. Gene ontology enrichment showed that transcription factors of different modules were enriched in muscle tissue development, connective tissue development, embryonic organ development, skeletal system morphogenesis, pattern specification process and urogenital system development separately. Six regulatory networks were constructed with key transcription factors, which contribute to the development of different tissues. Knockdown of four transcription factors from regulatory networks, including Sox9, Twist1, Snai2 and Klf4, showed that the expression of limb-development-related genes was also inhibited, which indicated the crucial role of transcription factor networks in hindlimb development.

Transgenic force sensors and software to measure force transmission across the mammalian nuclear envelope in vivo

Nuclear mechanotransduction is a growing field with exciting implications for the regulation of gene expression and cellular function. Mechanical signals may be transduced to the nuclear interior biochemically or physically through connections between the cell surface and chromatin. To define mechanical stresses upon the nucleus in physiological settings, we generated transgenic mouse strains that harbour FRET-based tension sensors or control constructs in the outer and inner aspects of the nuclear envelope. We knocked-in a published esprin-2G sensor to measure tensions across the LINC complex and generated a new sensor that links the inner nuclear membrane to chromatin. To mitigate challenges inherent to fluorescence lifetime analysis in vivo, we developed software (FLIMvivo) that markedly improves the fitting of fluorescence decay curves. In the mouse embryo, the sensors responded to cytoskeletal relaxation and stretch applied by micro-aspiration. They reported organ-specific differences and a spatiotemporal tension gradient along the proximodistal axis of the limb bud, raising the possibility that mechanical mechanisms coregulate pattern formation. These mouse strains and software are potentially valuable tools for testing and refining mechanotransduction hypotheses in vivo.