Pitx1 directly modulates the core limb development program to implement hindlimb identity

Thalidomide-induced limb malformations: an update and reevaluation

Historically, thalidomide-induced congenital malformations have served as an important example of the enhanced susceptibility of developing embryos to chemical perturbation. The compound produced a wide variety of congenital malformations in humans, which were initially detected by an association with a relatively rare limb defect labeled phocomelia. Although true phocomelia in the most severe form is a transverse defect with intercalary absence of limb regions, it is proposed that thalidomide produces a longitudinal limb phenotype in humans under usual circumstances that can become transverse in severe cases with a preferential sensitivity of forelimb over hindlimb, preaxial over postaxial, and left more impacted than the corresponding non-autopod limb bones on the right. The thalidomide-induced limb phenotype in humans is described and followed by a hierarchical comparison with various laboratory animal species. Mechanistic studies have been hampered by the fact that only non-human primates and rabbits have malformations that are anatomically similar to humans. Included in this review are unpublished data on limb malformations produced by thalidomide in rhesus monkeys from experiments performed more than 50 years ago. The critical period in gestation for the induction of phocomelia may initiate prior to the development of the embryonic limb bud, which contrasts with other chemical and physical agents that are known to produce this phenotype. The importance of toxicokinetic parameters is reviewed including dose, enantiomers, absorption, distribution, and both non-enzymatic and enzymatic biotransformations. The limb embryopathy mechanism that provides a partial explanation of the limb phenotype is that cereblon binds to thalidomide creating a protein complex that ubiquitinates protein substrates (CRL4CRBN) that are not targets for the complex in the absence of the thalidomide. One of these neosubstrates is SALL4 which when mutated causes a syndrome that phenocopies aspects of thalidomide embryopathy. Other candidate neosubstrates for the complex that have been found in non-human species may contribute to an understanding of the limb defect including PLZF, p63, and various zinc finger transcription factors. It is proposed that it is important to consider the species-specificity of the compound when considering potential mechanistic pathways and that some of the more traditional mechanisms for explaining the embryopathy, such as anti-angiogenesis and redox perturbation, may contribute to a full understanding of this teratogen.

Construction of an atlas of transcription factor binding during mouse development identifies popular regulatory regions

Gene regulators physically associate with the genome, in a combinatorial fashion, to drive tissue-specific gene expression. Uncovering the genome-wide activity of all gene regulators across tissues is therefore needed to understand gene regulation during development. Here, we take a first step towards this goal. Using CUT&RUN, we systematically mapped genome-wide binding profiles of key transcription factors and co-factors that mediate ontogenetically relevant signaling pathways in select mouse tissues at two developmental stages. Computation of the datasets unveiled tissue- and time-specific activity for each gene regulator. We identified 'popular' regulatory regions that are bound by a multitude of regulators, which tend to be more evolutionarily conserved. Consistently, they lie near the transcription start site of genes for which dysregulation results in early embryonic lethality. Moreover, the human homologs of these regions are similarly bound by many gene regulators and are highly conserved, indicating a retained relevance for human development. This work constitutes a decisive step towards understanding how the genome is simultaneously read and used by gene regulators in a holistic fashion to drive embryonic development.

Functional genomics of human skeletal development and the patterning of height heritability

Underlying variation in height are regulatory changes to chondrocytes, cartilage cells comprising long-bone growth plates. Currently, we lack knowledge on epigenetic regulation and gene expression of chondrocytes sampled across the human skeleton, and therefore we cannot understand basic regulatory mechanisms controlling height biology. We first rectify this issue by generating extensive epigenetic and transcriptomic maps from chondrocytes sampled from different growth plates across developing human skeletons, discovering novel regulatory networks shaping human bone/joint development. Next, using these maps in tandem with height genome-wide association study (GWAS) signals, we disentangle the regulatory impacts that skeletal element-specific versus global-acting variants have on skeletal growth, revealing the prime importance of regulatory pleiotropy in controlling height variation. Finally, as height is highly heritable, and thus often the test case for complex-trait genetics, we leverage these datasets within a testable omnigenic model framework to discover novel chondrocyte developmental modules and peripheral-acting factors shaping height biology and skeletal growth.

Single-nucleus transcriptomics reveal the differentiation trajectories of periosteal skeletal/stem progenitor cells in bone regeneration

Bone regeneration is mediated by skeletal stem/progenitor cells (SSPCs) that are mainly recruited from the periosteum after bone injury. The composition of the periosteum and the steps of SSPC activation and differentiation remain poorly understood. Here, we generated a single-nucleus atlas of the periosteum at steady state and of the fracture site during the early stages of bone repair (https://fracture-repair-atlas.cells.ucsc.edu). We identified periosteal SSPCs expressing stemness markers (Pi16 and Ly6a/SCA1) and responding to fracture by adopting an injury-induced fibrogenic cell (IIFC) fate, prior to undergoing osteogenesis or chondrogenesis. We identified distinct gene cores associated with IIFCs and their engagement into osteogenesis and chondrogenesis involving Notch, Wnt, and the circadian clock signaling, respectively. Finally, we show that IIFCs are the main source of paracrine signals in the fracture environment, suggesting a crucial paracrine role of this transient IIFC population during fracture healing. Overall, our study provides a complete temporal topography of the early stages of fracture healing and the dynamic response of periosteal SSPCs to injury, redefining our knowledge of bone regeneration.

Development of a ZRS Reporter System for the Newt (Cynops pyrrhogaster) During Terrestrial Limb Regeneration

BACKGROUND

Newts, a type of urodele amphibian, offer remarkable insights into regenerative medicine due to their extraordinary tissue regeneration capabilities-a challenging feat in humans. During limb regeneration of adult newts, fascinating cellular and molecular processes are revealed, including scarless healing, de-differentiation of mature cells, and regeneration of limbs and digits. Sonic hedgehog (Shh), crucial for vertebrate limb development, is regulated by the zone of polarizing activity regulatory sequence (ZRS) in the limb bud zone of polarizing activity (ZPA). The metamorphosed (terrestrial) newt can reactivate Shh during regeneration, facilitating proper limb patterning. Cell types capable of regulating the ZRS in metamorphosed newts remain unknown. The identification of such cell types provides invaluable insight into novel regenerative mechanisms.

OBJECTIVE

In this study, we developed the first newt ZRS reporter.

METHODS

We isolated and characterized the newt ZRS enhancer (nZRS), identifying conserved DNA binding sites. Several binding sites with medical relevance were conserved in the newt ZRS. In functional analysis, we developed a system composed of a transgenic nZRS reporter newt and a new newt anti-Shh antibody, which allowed Shh monitoring during limb regeneration.

RESULTS

We identified a group of Schwann cells capable of ZRS reporter and Shh protein expression during terrestrial limb regeneration.

CONCLUSIONS

This system provides a valuable in vivo approach for future genetic studies of patterning during limb regeneration.

Single-cell sequencing of the vermiform appendix during development identifies transcriptional relationships with appendicitis in preschool children

BACKGROUND

The development of the human vermiform appendix at the cellular level, as well as its function, is not well understood. Appendicitis in preschool children, although uncommon, is associated with a high perforation rate and increased morbidity.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) on the human appendix during fetal and pediatric stages as well as preschool-age inflammatory appendices. Transcriptional features of each cell compartment were discussed in the developing appendix. Cellular interactions and differentiation trajectories were also investigated. We compared scRNA-seq profiles from preschool appendicitis to those of matched healthy controls to reveal disease-associated changes. Bulk transcriptomic data, immunohistochemistry, and real-time quantitative PCR were used to validate the findings.

RESULTS

Our analysis identified 76 cell types in total and described the cellular atlas of the developing appendix. We discovered the potential role of the BMP signaling pathway in appendiceal epithelium development and identified HOXC8 and PITX2 as the specific regulons of appendix goblet cells. Higher pericyte coverage, endothelial angiogenesis, and goblet mucus scores together with lower epithelial and endothelial tight junction scores were found in the preschool appendix, which possibly contribute to the clinical features of preschool appendicitis. Preschool appendicitis scRNA-seq profiles revealed that the interleukin-17 signaling pathway may participate in the inflammation process.

CONCLUSIONS

Our study provides new insights into the development of the appendix and deepens the understanding of appendicitis in preschool children.

Single nuclei transcriptomics reveal the differentiation trajectories of periosteal skeletal/stem progenitor cells in bone regeneration

Bone regeneration is mediated by skeletal stem/progenitor cells (SSPCs) that are mainly recruited from the periosteum after bone injury. The composition of the periosteum and the steps of SSPC activation and differentiation remain poorly understood. Here, we generated a single-nuclei atlas of the periosteum at steady-state and of the fracture site during early stages of bone repair ( https://fracture-repair-atlas.cells.ucsc.edu ). We identified periosteal SSPCs expressing stemness markers ( Pi16 and Ly6a /SCA1) and responding to fracture by adopting an injury-induced fibrogenic cell (IIFC) fate, prior to undergoing osteogenesis or chondrogenesis. We identified distinct gene cores associated with IIFCs and their engagement into osteogenesis and chondrogenesis involving Notch, Wnt and the circadian clock signaling respectively. Finally, we show that IIFCs are the main source of paracrine signals in the fracture environment, suggesting a crucial paracrine role of this transient IIFC population during fracture healing. Overall, our study provides a complete temporal topography of the early stages of fracture healing and the dynamic response of periosteal SSPCs to injury, redefining our knowledge of bone regeneration.

Chromatin profiling identifies chondrocyte-specific Sox9 enhancers important for skeletal development

The transcription factor SRY-related HMG box 9 (Sox9) is essential for chondrogenesis. Mutations in and around SOX9 cause campomelic dysplasia (CD) characterized by skeletal malformations. Although the function of Sox9 in this context is well studied, the mechanisms that regulate Sox9 expression in chondrocytes remain to be elucidated. Here, we have used genome-wide profiling to identify 2 Sox9 enhancers located in a proximal breakpoint cluster responsible for CD. Enhancer activity of E308 (located 308 kb 5' upstream) and E160 (located 160 kb 5' upstream) correlated with Sox9 expression levels, and both enhancers showed a synergistic effect in vitro. While single deletions in mice had no apparent effect, simultaneous deletion of both E308 and E160 caused a dwarf phenotype, concomitant with a reduction of Sox9 expression in chondrocytes. Moreover, bone morphogenetic protein 2-dependent chondrocyte differentiation of limb bud mesenchymal cells was severely attenuated in E308/E160 deletion mice. Finally, we found that an open chromatin region upstream of the Sox9 gene was reorganized in the E308/E160 deletion mice to partially compensate for the loss of E308 and E160. In conclusion, our findings reveal a mechanism of Sox9 gene regulation in chondrocytes that might aid in our understanding of the pathophysiology of skeletal disorders.

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Sall genes regulate hindlimb initiation in mouse embryos

Vertebrate limbs start to develop as paired protrusions from the lateral plate mesoderm at specific locations of the body with forelimb buds developing anteriorly and hindlimb buds posteriorly. During the initiation process, limb progenitor cells maintain active proliferation to form protrusions and start to express Fgf10, which triggers molecular processes for outgrowth and patterning. Although both processes occur in both types of limbs, forelimbs (Tbx5), and hindlimbs (Isl1) utilize distinct transcriptional systems to trigger their development. Here, we report that Sall1 and Sall4, zinc finger transcription factor genes, regulate hindlimb initiation in mouse embryos. Compared to the 100% frequency loss of hindlimb buds in TCre; Isl1 conditional knockouts, Hoxb6Cre; Isl1 conditional knockout causes a hypomorphic phenotype with only approximately 5% of mutants lacking the hindlimb. Our previous study of SALL4 ChIP-seq showed SALL4 enrichment in an Isl1 enhancer, suggesting that SALL4 acts upstream of Isl1. Removing 1 allele of Sall4 from the hypomorphic Hoxb6Cre; Isl1 mutant background caused loss of hindlimbs, but removing both alleles caused an even higher frequency of loss of hindlimbs, suggesting a genetic interaction between Sall4 and Isl1. Furthermore, TCre-mediated conditional double knockouts of Sall1 and Sall4 displayed a loss of expression of hindlimb progenitor markers (Isl1, Pitx1, Tbx4) and failed to develop hindlimbs, demonstrating functional redundancy between Sall1 and Sall4. Our data provides genetic evidence that Sall1 and Sall4 act as master regulators of hindlimb initiation.

New findings on the genetic basis of feathered legs in chickens: association of CUBN gene mutations with feathered-leg phenotype

Chickens are the most thoroughly domesticated vertebrate species, and after long-continued natural and artificial selection, they now show rich phenotypic diversity. In particular, feathered legs present in domestic chickens are a characteristic that is carefully selected by advanced breeders. Previous studies have identified the key mutations responsible for feathered legs on chromosomes 13 and 15; however, not all chickens can be easily distinguished based on these two markers. In this study, whole-genome resequencing of 29 Bamaxiaogu chickens (BXCs) yielded 12,201,978 valid single-nucleotide polymorphisms (SNPs) and 2,792,426 valid insertions and deletions (InDels). Population structure analysis based on SNPs revealed that the test samples came from the same natural population. Based on these findings, we used an SNP- and InDel-based genome-wide association study (GWAS) to investigate the genetic basis of feathered legs in chickens. GWAS results revealed that 2 SNPs located in the introns of cubilin (CUBN; SNP1, chr2:19885382T>A) and recombinant Ras suppressor protein 1 (RSU1) genes (SNP2, chr2:20002551G>A), as well as an InDel (InDel1, chr2:19884383TG>T) on CUBN, were all significantly associated with the presence of feathered legs. Diagnostic testing demonstrated that SNP1 effectively differentiated between chickens with feathered legs and those with clean legs (leg without feathers) within the BXC population and may thus be considered an effective marker of feathered legs in BXC. In contrast, other loci did not show the same discriminatory power. This study not only presents a new variant of feathered legs but also provides valuable novel insights into the underlying mechanisms of variation in the feathered-legs trait among chickens.

Identification of potential key variants in mandibular premolar hypodontia through whole-exome sequencing

Determining genotype-phenotype correlations in patients with hypodontia is important for understanding disease pathogenesis, although only a few studies have elucidated it. We aimed to identify genetic variants linked to non-syndromic bilateral mandibular second premolar hypodontia in a Korean population for the first time by specifying the phenotype of hypodontia. Twenty unrelated individuals with non-syndromic bilateral mandibular second premolar hypodontia were enrolled for whole-exome sequencing. Using a tooth agenesis gene set panel consisting of 112 genes based on literature, potential candidate variants were screened through variant filtering and prioritization. We identified 13 candidate variants in 12 genes, including a stop-gain variant (c.4750C>T) in LAMA3. Through the functional enrichment analysis of the prioritized genes, several terms related to tooth development were enriched in a protein-protein interaction network of candidate genes for mandibular premolar hypodontia. The hypodontia group also had approximately 2-fold as many mutated variants in all four genes related to these key terms, which are CDH1, ITGB4, LAMA3, LAMB3, as those in the 100 healthy control group individuals. The relationship between enriched terms and pathways and mandibular premolar hypodontia was also investigated. In addition, we identified some known oligodontia variants in patients with hypodontia, strengthening the possibility of synergistic effects in other genes. This genetic investigation may be a worthwhile preliminary attempt to reveal the pathogenesis of tooth agenesis and sets a background for future studies.

Cell-specific alterations in Pitx1 regulatory landscape activation caused by the loss of a single enhancer

Developmental genes are frequently controlled by multiple enhancers sharing similar specificities. As a result, deletions of such regulatory elements have often failed to reveal their full function. Here, we use the Pitx1 testbed locus to characterize in detail the regulatory and cellular identity alterations following the deletion of one of its enhancers (Pen). By combining single cell transcriptomics and an in-embryo cell tracing approach, we observe an increased fraction of Pitx1 non/low-expressing cells and a decreased fraction of Pitx1 high-expressing cells. We find that the over-representation of Pitx1 non/low-expressing cells originates from a failure of the Pitx1 locus to coordinate enhancer activities and 3D chromatin changes. This locus mis-activation induces a localized heterochrony and a concurrent loss of irregular connective tissue, eventually leading to a clubfoot phenotype. This data suggests that, in some cases, redundant enhancers may be used to locally enforce a robust activation of their host regulatory landscapes.

Joint disease-specificity at the regulatory base-pair level

Given the pleiotropic nature of coding sequences and that many loci exhibit multiple disease associations, it is within non-coding sequence that disease-specificity likely exists. Here, we focus on joint disorders, finding among replicated loci, that GDF5 exhibits over twenty distinct associations, and we identify causal variants for two of its strongest associations, hip dysplasia and knee osteoarthritis. By mapping regulatory regions in joint chondrocytes, we pinpoint two variants (rs4911178; rs6060369), on the same risk haplotype, which reside in anatomical site-specific enhancers. We show that both variants have clinical relevance, impacting disease by altering morphology. By modeling each variant in humanized mice, we observe joint-specific response, correlating with GDF5 expression. Thus, we uncouple separate regulatory variants on a common risk haplotype that cause joint-specific disease. By broadening our perspective, we finally find that patterns of modularity at GDF5 are also found at over three-quarters of loci with multiple GWAS disease associations.

Pitx genes in development and disease

Homeobox genes encode sequence-specific transcription factors (SSTFs) that recognize specific DNA sequences and regulate organogenesis in all eukaryotes. They are essential in specifying spatial and temporal cell identity and as a result, their mutations often cause severe developmental defects. Pitx genes belong to the PRD class of the highly evolutionary conserved homeobox genes in all animals. Vertebrates possess three Pitx paralogs, Pitx1, Pitx2, and Pitx3 while non-vertebrates have only one Pitx gene. The ancient role of regulating left-right (LR) asymmetry is conserved while new functions emerge to afford more complex body plan and functionalities. In mouse, Pitx1 regulates hindlimb tissue patterning and pituitary development. Pitx2 is essential for the development of the oral cavity and abdominal wall while regulates the formation and symmetry of other organs including pituitary, heart, gut, lung among others by controlling growth control genes upon activation of the Wnt/ß-catenin signaling pathway. Pitx3 is essential for lens development and migration and survival of the dopaminergic neurons of the substantia nigra. Pitx gene mutations are linked to various congenital defects and cancers in humans. Pitx gene family has the potential to offer a new approach in regenerative medicine and aid in identifying new drug targets.

Control of mouse limb initiation and antero-posterior patterning by Meis transcription factors

Meis1 and Meis2 are homeodomain transcription factors that regulate organogenesis through cooperation with Hox proteins. Elimination of Meis genes after limb induction has shown their role in limb proximo-distal patterning; however, limb development in the complete absence of Meis function has not been studied. Here, we report that Meis1/2 inactivation in the lateral plate mesoderm of mouse embryos leads to limb agenesis. Meis and Tbx factors converge in this function, extensively co-binding with Tbx to genomic sites and co-regulating enhancers of Fgf10, a critical factor in limb initiation. Limbs with three deleted Meis alleles show proximal-specific skeletal hypoplasia and agenesis of posterior skeletal elements. This failure in posterior specification results from an early role of Meis factors in establishing the limb antero-posterior prepattern required for Shh activation. Our results demonstrate roles for Meis transcription factors in early limb development and identify their involvement in previously undescribed interaction networks that regulate organogenesis.

Establishing the pattern of the vertebrate limb

The vertebrate limb continues to serve as an influential model of growth, morphogenesis and pattern formation. With this Review, we aim to give an up-to-date picture of how a population of undifferentiated cells develops into the complex pattern of the limb. Focussing largely on mouse and chick studies, we concentrate on the positioning of the limbs, the formation of the limb bud, the establishment of the principal limb axes, the specification of pattern, the integration of pattern formation with growth and the determination of digit number. We also discuss the important, but little understood, topic of how gene expression is interpreted into morphology.

A single cell transcriptional atlas of early synovial joint development

Synovial joint development begins with the formation of the interzone, a region of condensed mesenchymal cells at the site of the prospective joint. Recently, lineage-tracing strategies have revealed that Gdf5-lineage cells native to and from outside the interzone contribute to most, if not all, of the major joint components. However, there is limited knowledge of the specific transcriptional and signaling programs that regulate interzone formation and fate diversification of synovial joint constituents. To address this, we have performed single cell RNA-Seq analysis of 7329 synovial joint progenitor cells from the developing murine knee joint from E12.5 to E15.5. By using a combination of computational analytics, in situ hybridization and in vitro characterization of prospectively isolated populations, we have identified the transcriptional profiles of the major developmental paths for joint progenitors. Our freely available single cell transcriptional atlas will serve as a resource for the community to uncover transcriptional programs and cell interactions that regulate synovial joint development.

Evolutionary Selection and Constraint on Human Knee Chondrocyte Regulation Impacts Osteoarthritis Risk

During human evolution, the knee adapted to the biomechanical demands of bipedalism by altering chondrocyte developmental programs. This adaptive process was likely not without deleterious consequences to health. Today, osteoarthritis occurs in 250 million people, with risk variants enriched in non-coding sequences near chondrocyte genes, loci that likely became optimized during knee evolution. We explore this relationship by epigenetically profiling joint chondrocytes, revealing ancient selection and recent constraint and drift on knee regulatory elements, which also overlap osteoarthritis variants that contribute to disease heritability by tending to modify constrained functional sequence. We propose a model whereby genetic violations to regulatory constraint, tolerated during knee development, lead to adult pathology. In support, we discover a causal enhancer variant (rs6060369) present in billions of people at a risk locus (GDF5-UQCC1), showing how it impacts mouse knee-shape and osteoarthritis. Overall, our methods link an evolutionarily novel aspect of human anatomy to its pathogenesis.

Developmental and Evolutionary Allometry of the Mammalian Limb Skeleton

The variety of limb skeletal proportions enables a remarkable diversity of behaviors that include powered flight in bats and flipper-propelled swimming in whales using extremes of a range of homologous limb architectures. Even within human limbs, bone lengths span more than an order of magnitude from the short finger and toe bones to the long arm and leg bones. Yet all of this diversity arises from embryonic skeletal elements that are each a very similar size at formation. In this review article, I survey what is and is not yet known of the development and evolution of skeletal proportion at multiple hierarchical levels of biological organization. These include the cellular parameters of skeletal elongation in the cartilage growth plate, genes associated with differential growth, and putative gene regulatory mechanisms that would allow both covariant and independent evolution of the forelimbs and hindlimbs and of individual limb segments. Although the genetic mechanisms that shape skeletal proportion are still largely unknown, and most of what is known is limited to mammals, it is becoming increasingly apparent that the diversity of bone lengths is an emergent property of a complex system that controls elongation of individual skeletal elements using a genetic toolkit shared by all.

Retinoic Acid Promotes Retinoic Acid Signaling by Suppression of Pitx1 In Tendon Cells: A Possible Mechanism of a Clubfoot-Like Phenotype Induced by Retinoic Acid

Background

The pathogenesis of idiopathic congenital clubfoot (CCF) is unknown. Although some familial patients have Pitx1 mutations, and the Pitx1+/− genotype causes a clubfoot-like phenotype in mice, the mechanism of Pitx1-induced CCF is unknown.

Material/Methods

We used tibialis anterior tendon samples to detect the expression of Pitx1 in idiopathic and neurogenic clubfoot patients. After obtaining Sprague-Dawley (SD) rat Achilles tendon cells, the expression of Pitx1 was knocked down by SiRNA. After 48 h of culture, mass spectrometry was used to quantitatively analyze proteins. Then, Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used to assess the downstream pathway of PITX1. The relationship between Pitx1 and the promoter region of deacetylase 1 (Sirtuin-1 and Sirt1) was examined by luciferase and ChIP assays.

Results

We found that Pitx1 expression in the tendon samples of idiopathic CCF patients was downregulated. Mass spectrometry analysis revealed that the inhibition of Pitx1 induced the downregulation of Sirt1 expression in tendon cells. Luciferase and ChIP assays confirmed that Pitx1 binds to the promoter region of SIRT1 and promotes Sirt1 gene transcription. Further results showed that, after the inhibition of Pitx1 in tendon cells, CRABP2 acetylation increased, the nuclear import of CRABP2 was enhanced, and the expression of RARβ2 increased. After the inhibition of Pitx1, RARβ2 expression was further increased by RA treatment in tendon cells. In the presence of retinoic acid, the expression of Pitx1 was inhibited in tendon cells.

Conclusions

Pitx1 binds to the promoter region of SIRT1 and promotes the transcription of SIRT1. Positive feedback occurs between RA signaling and Pitx1.

Pigeon foot feathering reveals conserved limb identity networks

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

Similarities and differences in the regulation of HoxD genes during chick and mouse limb development

In all tetrapods examined thus far, the development and patterning of limbs require the activation of gene members of the HoxD cluster. In mammals, they are regulated by a complex bimodal process that controls first the proximal patterning and then the distal structure. During the shift from the former to the latter regulation, this bimodal regulatory mechanism allows the production of a domain with low Hoxd gene expression, at which both telomeric (T-DOM) and centromeric regulatory domains (C-DOM) are silent. These cells generate the future wrist and ankle articulations. We analyzed the implementation of this regulatory mechanism in chicken, i.e., in an animal for which large morphological differences exist between fore- and hindlimbs. We report that although this bimodal regulation is globally conserved between the mouse and the chick, some important modifications evolved at least between these two model systems, in particular regarding the activity of specific enhancers, the width of the TAD boundary separating the two regulations, and the comparison between the forelimb versus hindlimb regulatory controls. At least one aspect of these regulations seems to be more conserved between chick and bats than with mouse, which may relate to the extent to which forelimbs and hindlimbs of these various animals differ in their morphologies.

A comparison of Hox gene regulation during the development of limbs in birds and mammals reveals that whereas the characteristic bimodal regulatory system, based on large chromatin domains, is largely conserved between these morphologically distinct structures, some differences are revealed in the way this is implemented in various vertebrates.

The shapes of limbs vary greatly among tetrapod species, even between the forelimbs and hindlimbs of the same animal. Hox genes regulate the proper growth and patterning of tetrapod limbs. In order to evaluate whether variations in the complex regulation of a cluster of Hox genes—the Hoxd genes—during limb development contribute to the differences in limb shape, we compared their transcriptional control during limb bud development in the forelimbs and hindlimbs of mouse and chicken embryos. We found that the regulatory mechanism underlying Hoxd gene expression is highly conserved, but some clear differences exist. For instance, we observed a variation in the topologically associating domain (TAD; a self-interacting genomic region) boundary interval between the mouse and the chick, as well as differences in the activity of a conserved enhancer element situated within the telomeric regulatory domain. In contrast to the mouse, the chicken enhancer has a stronger activity in the forelimb buds than in the hindlimb buds, which is correlated with the striking differences in the mRNA levels of the genes. We conclude that differences in both the timing and duration of TAD activities and in the width of their boundary may parallel the important decrease in Hoxd gene transcription in chick hindlimb buds versus forelimb buds. These differences may also account for the slightly distinct regulatory strategies implemented by mammals and birds at this locus.

Dynamic 3D chromatin architecture contributes to enhancer specificity and limb morphogenesis

The regulatory specificity of enhancers and their interaction with gene promoters is thought to be controlled by their sequence and the binding of transcription factors. By studying Pitx1, a regulator of hindlimb development, we show that dynamic changes in chromatin conformation can restrict the activity of enhancers. Inconsistent with its hindlimb-restricted expression, Pitx1 is controlled by an enhancer (Pen) that shows activity in forelimbs and hindlimbs. By Capture Hi-C and three-dimensional modeling of the locus, we demonstrate that forelimbs and hindlimbs have fundamentally different chromatin configurations, whereby Pen and Pitx1 interact in hindlimbs and are physically separated in forelimbs. Structural variants can convert the inactive into the active conformation, thereby inducing Pitx1 misexpression in forelimbs, causing partial arm-to-leg transformation in mice and humans. Thus, tissue-specific three-dimensional chromatin conformation can contribute to enhancer activity and specificity in vivo and its disturbance can result in gene misexpression and disease.

Shh signaling influences the phenotype of Pitx1-/- hindlimbs

Forelimbs (FLs) and hindlimbs (HLs) develop under the instructive and integrated guidance of signaling centers and transcription factor (TF) action. The development of structures specific to each limb type depends on the limb-specific modulation of these integrated components. Pitx1 is a transcription factor gene expressed in HL, absent in FL, and required for HL-specific patterning and development, in particular for formation of anterior HL skeletal elements. Pitx1 achieves this function by direct TF action on the core limb program, which is largely shared between FL and HL. Shh signaling plays a crucial role in anterior-posterior (AP) patterning in both FL and HL. The present work assessed the relationship between Shh signaling and Pitx1 action for AP patterning. We found that reducing the gene dosage of Shh in the context of the Pitx1-/- HL decreases the severity of the Pitx1-/- phenotype, in particular, the loss of anterior limb structures and the shortening of femur length. However, this did not rescue HL-specific patterning features. Thus, Pitx1 action integrates Shh signaling but not for limb-type-specific patterning.

Regulatory integration of Hox factor activity with T-box factors in limb development

In tetrapods, Tbx4, Tbx5 and Hox cluster genes are crucial for forelimb and hindlimb development and mutations in these genes are responsible for congenital limb defects. The molecular basis of their integrated mechanisms of action in the context of limb development remains poorly understood. We studied Tbx4 and Hoxc10 owing to their overlapping loss-of-function phenotypes and colocalized expression in mouse hindlimb buds. We report an extensive overlap between Tbx4 and Hoxc10 genome occupancy and their putative target genes. Tbx4 and Hoxc10 interact directly with each other, have the ability to bind to a previously unrecognized T-box-Hox composite DNA motif and show synergistic activity when acting on reporter genes. Pitx1, the master regulator for hindlimb specification, also shows extensive genomic colocalization with Tbx4 and Hoxc10. Genome occupancy by Tbx4 in hindlimb buds is similar to Tbx5 occupancy in forelimbs. By contrast, another Hox factor, Hoxd13, also interacts with Tbx4/Tbx5 but antagonizes Tbx4/Tbx5-dependent transcriptional activity. Collectively, the modulation of Tbx-dependent activity by Hox factors acting on common DNA targets may integrate different developmental processes for the balanced formation of proportionate limbs.