In tetrapods, the scapular and pelvic girdles perform the important function of anchoring the limbs to the trunk of the body and facilitating the movement of each appendage. This shared function, however, is one of relatively few similarities between the scapula and pelvis, which have significantly different morphologies, evolutionary histories, embryonic origins, and underlying genetic pathways. The scapula evolved in jawless fish prior to the pelvis, and its embryonic development is unique among bones in that it is derived from multiple progenitor cell populations, including the dermomyotome, somatopleure, and neural crest. Conversely, the pelvis evolved several million years later in jawed fish, and it develops from an embryonic somatopleuric cell population. The genetic networks controlling the formation of the pelvis and scapula also share similarities and differences, with a number of genes shaping only one or the other, while other gene products such as PBX transcription factors act as hierarchical developmental regulators of both girdle structures. Here, we provide a detailed review of the cellular processes and genetic networks underlying pelvis and scapula formation in tetrapods, while also highlighting unanswered questions about girdle evolution and development.

ABSTRACT

During development cell fates are specified by tightly controlled gene expression programs. PBX TALE transcription factors control gene regulatory networks (GRN) that direct vertebrate tissue patterning and organ morphogenesis. How PBX1/2 proteins acquire context-specific functions, despite widespread embryonic expression of Pbx1/2 , remains elusive. In mouse limb buds, mesenchymal-specific loss of PBX1/2 or of the transcriptional regulator HAND2 results in similar phenotypes, suggesting that PBX1/2- and HAND2-dependent programs converge to control limb development. To investigate this scenario, we combined tissue-specific and temporally-controlled mutagenesis with multi-omics approaches using the murine hindlimb bud as a model. We reconstructed a GRN collaboratively directed by PBX1/2 and HAND2, demonstrating that Pbx1 - Hand2 genetically interact in vivo during hindlimb patterning, with PBX1 concomitantly acting as an upstream regulator of Hand2 . At organismal-level resolution the GRN is active within restricted subsets of posterior-proximal hindlimb mesenchymal cells, wherein Pbx1/2 and Hand2 are co-expressed with their target genes. Genome-wide profiling of PBX1 binding across multiple tissues further revealed that HAND2 selects a subset of PBX-bound regions to impart limb patterning functionality. This research elucidates mechanisms underlying limb bud-specific functions by PBX1/2, while informing general principles by which promiscuous transcription factors cooperate with select cofactors to instruct distinct developmental programs.

Despite the great diversity of vertebrate limb proportion and our deep understanding of the genetic mechanisms that drive skeletal elongation, little is known about how individual bones reach different lengths in any species. Here, we directly compare the transcriptomes of homologous growth cartilages of the mouse (Mus musculus) and bipedal jerboa (Jaculus jaculus), the latter of which has "mouse-like" arms but extremely long metatarsals of the feet. Intersecting gene-expression differences in metatarsals and forearms of the two species revealed that about 10% of orthologous genes are associated with the disproportionately rapid elongation of neonatal jerboa feet. These include genes and enriched pathways not previously associated with endochondral elongation as well as those that might diversify skeletal proportion in addition to their known requirements for bone growth throughout the skeleton. We also identified transcription regulators that might act as "nodes" for sweeping differences in genome expression between species. Among these, Shox2, which is necessary for proximal limb elongation, has gained expression in jerboa metatarsals where it has not been detected in other vertebrates. We show that Shox2 is sufficient to increase mouse distal limb length, and a nearby putative cis-regulatory region is preferentially accessible in jerboa metatarsals. In addition to mechanisms that might directly promote growth, we found evidence that jerboa foot elongation may occur in part by de-repressing latent growth potential. The genes and pathways that we identified here provide a framework to understand the modular genetic control of skeletal growth and the remarkable malleability of vertebrate limb proportion.

The connection between different tissues is vital for the development and function of any organs and systems. In the musculoskeletal system, the attachment of elastic tendons to stiff bones poses a mechanical challenge that is solved by the formation of a transitional tissue, which allows the transfer of muscle forces to the skeleton without tearing. Here, we show that tendon-to-bone attachment cells are bi-fated, activating a mixture of chondrocyte and tenocyte transcriptomes, which is regulated by sharing regulatory elements with these cells and by Krüppel-like factors transcription factors (KLF). To uncover the molecular identity of attachment cells, we first applied high-throughput RNA sequencing to murine humeral attachment cells. The results, which were validated by in situ hybridization and single-molecule in situ hybridization, reveal that attachment cells express hundreds of chondrogenic and tenogenic genes. In search for the underlying mechanism allowing these cells to express these genes, we performed ATAC sequencing and found that attachment cells share a significant fraction of accessible intergenic chromatin areas with either tenocytes or chondrocytes. Epigenomic analysis further revealed transcriptional enhancer signatures for the majority of these regions. We then examined a subset of these regions using transgenic mouse enhancer reporter. Results verified the shared activity of some of these enhancers, supporting the possibility that the transcriptome of attachment cells is regulated by enhancers with shared activities in tenocytes or chondrocytes. Finally, integrative chromatin and motif analyses, as well as the transcriptome data, indicated that KLFs are regulators of attachment cells. Indeed, blocking the expression of Klf2 and Klf4 in the developing limb mesenchyme led to abnormal differentiation of attachment cells, establishing these factors as key regulators of the fate of these cells. In summary, our findings show how the molecular identity of bi-fated attachment cells enables the formation of the unique transitional tissue that connect tendon to bone. More broadly, we show how mixing the transcriptomes of two cell types through shared enhancers and a dedicated set of transcription factors can lead to the formation of a new cell fate that connects them.

The mechanical challenge of attaching elastic tendons to stiff bones is solved by the formation of a unique transitional tissue. Here, we show that murine tendon-to-bone attachment cells are bi-fated, activating a mixture of chondrocyte and tenocyte transcriptomes, under regulation of shared regulatory elements and Krüppel-like factors (KLFs) transcription factors. High-throughput bulk and single-cell RNA sequencing of humeral attachment cells revealed expression of hundreds of chondrogenic and tenogenic genes, which was validated by in situ hybridization and single-molecule ISH. ATAC sequencing showed that attachment cells share accessible intergenic chromatin areas with either tenocytes or chondrocytes. Epigenomic analysis revealed enhancer signatures for most of these regions. Transgenic mouse enhancer reporter assays verified the shared activity of some of these enhancers. Finally, integrative chromatin and motif analyses and transcriptomic data implicated KLFs as regulators of attachment cells. Indeed, blocking expression of both Klf2 and Klf4 in developing limb mesenchyme impaired their differentiation.

Recent studies have reported evidence suggesting that portions of contemporary human genomes introgressed from archaic hominin populations went to high frequencies due to positive selection. However, no study to date has specifically addressed the postintrogression population dynamics of these putative cases of adaptive introgression. Here, for the first time, we specifically define cases of immediate adaptive introgression (iAI) in which archaic haplotypes rose to high frequencies in humans as a result of a selective sweep that occurred shortly after the introgression event. We define these cases as distinct from instances of selection on standing introgressed variation (SI), in which an introgressed haplotype initially segregated neutrally and subsequently underwent positive selection. Using a geographically diverse data set, we report novel cases of selection on introgressed variation in living humans and shortlist among these cases those whose selective sweeps are more consistent with having been the product of iAI rather than SI. Many of these novel inferred iAI haplotypes have potential biological relevance, including three that contain immune-related genes in West Siberians, South Asians, and West Eurasians. Overall, our results suggest that iAI may not represent the full picture of positive selection on archaically introgressed haplotypes in humans and that more work needs to be done to analyze the role of SI in the archaic introgression landscape of living humans.

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.

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10^-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r²  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Conserved genomic sequences disrupted in humans may underlie uniquely human phenotypic traits. We identified and characterized 10,032 human-specific conserved deletions (hCONDELs). These short (average 2.56 base pairs) deletions are enriched for human brain functions across genetic, epigenomic, and transcriptomic datasets. Using massively parallel reporter assays in six cell types, we discovered 800 hCONDELs conferring significant differences in regulatory activity, half of which enhance rather than disrupt regulatory function. We highlight several hCONDELs with putative human-specific effects on brain development, including HDAC5, CPEB4, and PPP2CA. Reverting an hCONDEL to the ancestral sequence alters the expression of LOXL2 and developmental genes involved in myelination and synaptic function. Our data provide a rich resource to investigate the evolutionary mechanisms driving new traits in humans and other species.