The temporal dynamics of vertebrate limb development, teratogenesis and evolution

Immobilization secondary to cell death of muscle precursors with a dual transcriptional signature contributes to the emu wing skeletal pattern

Limb reduction has occurred multiple times in tetrapod history. Among ratites, wing reductions range from mild vestigialization to complete loss, with emus (Dromaius novaehollandiae) serving as a model for studying the genetic mechanisms behind limb reduction. Here, we explore the developmental mechanisms underlying wing reduction in emu. Our analyses reveal that immobilization resulting from the absence of distal muscles contributes to skeletal shortening, fusion and left-right intraindividual variation. Expression analysis and single cell-RNA sequencing identify muscle progenitors displaying a dual lateral plate mesodermal and myogenic signature. These cells aggregate at the proximal region of wing buds and undergo cell death. We propose that this cell death, linked to the lack of distal muscle masses, underlines the morphological features and variability in skeletal elements due to reduced mechanical loading. Our results demonstrate that differential mobility during embryonic development may drive morphological diversification in vestigial structures.

Whole genome resequencing reveals candidate genes for postaxial polydactyly in Large White pigs

Polydactyly is a genetic abnormality that affects both pig welfare and industry profits. Despite efforts to explore the genetic basis of pig polydactyly, progress remains limited. In this study, we analyzed a group of Large White pigs with postaxial polydactyly, including 29 cases and 79 controls from 24 families. High-depth sequencing was performed on 20 pigs, while low-depth sequencing was improved through imputation for the remaining pigs. A genome-wide association study (GWAS) and genetic differentiation were conducted using the resequencing dataset, resulting in the identification of 48 significantly associated SNPs and 27 candidate regions. The genetic differentiation regions on chromosomes 5 and 18, which harbored GWAS-identified SNPs, were delineated as confidence regions. The confidence region at Chr18: 1.850-1.925 Mb covers the fifth intron of LMBR1, a gene that contains an important regulatory element for SHH, known as ZRS. Mutations in this ZRS have been found to cause polydactyly in animals and humans. Therefore, we propose LMBR1 as a prospective candidate gene for postaxial polydactyly. These findings emphasize the importance of exploring the role of ZRS within LMBR1 in the pathogenesis of polydactyly in pigs.

Developmental Genetic Basis of Hoxd9 Homeobox Domain Deletion in Pampus argenteus Pelvic Fin Deficiency

Pampus argenteus is important for commercial fishery catch species and is an emerging target for aquaculture production. Notably, P. argenteus has a bizarre morphology and lacks pelvic fins. However, the reason for the lack of pelvic fins remains unclear, ultimately leading to frequent upside-down floating of P. argenteus during breeding and marked consumption of physical energy. Some lineages, including whales, fugu, snakes, and seahorse, independently lost the pelvic appendages over evolutionary time. Do different taxa employ the same molecular genetic pathways when they independently evolve similar developmental morphologies? Through analysis of the gene responsible for appendage localization, Hoxd9, it was discovered that the Hox domain was absent in the Hoxd9 gene of P. argenteus, and the Hoxd9b gene lacked the Hox9 activation region, a feature not observed in the Hoxd9 gene of other fish species. Interestingly, those distinctive characteristics are not observed in the Hoxd9 gene of other fish species. To determine the association between the Hoxd9 gene characteristics and the pelvic fin deletion in P. argenteus, the full-length cDNA of the Hoxd9a gene was cloned, and morphological observations of the species' juveniles were performed using stereomicroscopy and scanning electron microscopy. Thereafter, the tissue localization of Hoxd9a in the species was analyzed at the gene and protein levels. Based on the results, deletion of the Hoxd9a structural domain possibly leads to disruptions in the protein translation and the pelvic fin localization in P. argenteus during its early ontogenetic developmental stage, resulting in the absence of pelvic fins.

Retinoids and developmental neurotoxicity: Utilizing toxicogenomics to enhance adverse outcome pathways and testing strategies

The use of genomic approaches in toxicological studies has greatly increased our ability to define the molecular profiles of environmental chemicals associated with developmental neurotoxicity (DNT). Integration of these approaches with adverse outcome pathways (AOPs), a framework that translates environmental exposures to adverse developmental phenotypes, can potentially inform DNT testing strategies. Here, using retinoic acid (RA) as a case example, we demonstrate that the integration of toxicogenomic profiles into the AOP framework can be used to establish a paradigm for chemical testing. RA is a critical regulatory signaling molecule involved in multiple aspects of mammalian central nervous system (CNS) development, including hindbrain formation/patterning and neuronal differentiation, and imbalances in RA signaling pathways are linked with DNT. While the mechanisms remain unresolved, environmental chemicals can cause DNT by disrupting the RA signaling pathway. First, we reviewed literature evidence of RA and other retinoid exposures and DNT to define a provisional AOP related to imbalances in RA embryonic bioavailability and hindbrain development. Next, by integrating toxicogenomic datasets, we defined a relevant transcriptomic signature associated with RA-induced developmental neurotoxicity (RA-DNT) in human and rodent models that was tested against zebrafish model data, demonstrating potential for integration into an AOP framework. Finally, we demonstrated how these approaches may be systematically utilized to identify chemical hazards by testing the RA-DNT signature against azoles, a proposed class of compounds that alters RA-signaling. The provisional AOP from this study can be expanded in the future to better define DNT biomarkers relevant to RA signaling and toxicity.

Whole embryo culture, transcriptomics and RNA interference identify TBX1 and FGF11 as novel regulators of limb development in the mouse

Identifying genes involved in vertebrate developmental processes and characterizing this involvement are daunting tasks, especially in the mouse where viviparity complicates investigations. Attempting to devise a streamlined approach for this type of study we focused on limb development. We cultured E10.5 and E12.5 embryos and performed transcriptional profiling to track molecular changes in the forelimb bud over a 6-hour time-window. The expression of certain genes was found to diverge rapidly from its normal path, possibly reflecting the activation of a stress-induced response. Others, however, maintained for up to 3 hours dynamic expression profiles similar to those seen in utero. Some of these resilient genes were known regulators of limb development. The implication of the others in this process was either unsuspected or unsubstantiated. The localized knockdown of two such genes, Fgf11 and Tbx1, hampered forelimb bud development, providing evidence of their implication. These results show that combining embryo culture, transcriptome analysis and RNA interference could speed up the identification of genes involved in a variety of developmental processes, and the validation of their implication.

Mathematical modeling of chondrogenic pattern formation during limb development: Recent advances in continuous models

The phenomenon of chondrogenic pattern formation in the vertebrate limb is one of the best studied examples of organogenesis. Many different models, mathematical as well as conceptual, have been proposed for it in the last fifty years or so. In this review, we give a brief overview of the fundamental biological background, then describe in detail several models which aim to describe qualitatively and quantitatively the corresponding biological phenomena. We concentrate on several new models that have been proposed in recent years, taking into account recent experimental progress. The major mathematical tools in these approaches are ordinary and partial differential equations. Moreover, we discuss models with non-local flux terms used to account for cell-cell adhesion forces and a structured population model with diffusion. We also include a detailed list of gene products and potential morphogens which have been identified to play a role in the process of limb formation and its growth.

Robinow syndrome skeletal phenotypes caused by the WNT5AC83S variant are due to dominant interference with chondrogenesis

Heterozygous missense mutations in several genes in the WNT5A signaling pathway cause autosomal dominant Robinow syndrome 1 (DRS1). Our objective was to clarify the functional impact of a missense mutation in WNT5A on the skeleton, one of the main affected tissues in RS. We delivered avian replication competent retroviruses (RCAS) containing human wild-type WNT5A (wtWNT5A), WNT5AC83S variant or GFP/AlkPO4 control genes to the chicken embryo limb. Strikingly, WNT5AC83S consistently caused a delay in ossification and bones were more than 50% shorter and 200% wider than controls. In contrast, bone dimensions in wtWNT5A limbs were slightly affected (20% shorter, 25% wider) but ossification occurred on schedule. The dysmorphology of bones was established during cartilage differentiation. Instead of stereotypical stacking of chondrocytes, the WNT5AC83S-infected cartilage was composed of randomly oriented chondrocytes and that had diffuse, rather than concentrated Prickle staining, both signs of disrupted planar cell polarity (PCP) mechanisms. Biochemical assays revealed that C83S variant was able to activate the Jun N-terminal kinase-PCP pathway similar to wtWNT5A; however, the activity of the variant ligand was influenced by receptor availability. Unexpectedly, the C83S change caused a reduction in the amount of protein being synthesized and secreted, compared to wtWNT5A. Thus, in the chicken and human, RS phenotypes are produced from the C83S mutation, even though the variant protein is less abundant than wtWNT5A. We conclude the variant protein has dominant-negative effects on chondrogenesis leading to limb abnormalities.

Evolution of the avian digital pattern

Variation in digit number has occurred multiple times in the history of archosaur evolution. The five digits of dinosaur limbs were reduced to three in bird forelimbs, and were further reduced in the vestigial forelimbs of the emu. Regulation of digit number has been investigated previously by examining genes involved in anterior-posterior patterning in forelimb buds among emu (Dromaius novaehollandiae), chicken (Gallus gallus) and zebra finch (Taeniopygia guttata). It was described that the expression of posterior genes are conserved among these three birds, whereas expression of anterior genes Gli3 and Alx4 varied significantly. Here we re-examined the expression pattern of Gli3 and Alx4 in the forelimb of emu, chicken and zebra finch. We found that Gli3 is expressed in the anterior region, although its range varied among species, and that the expression pattern of Alx4 in forelimb buds is broadly conserved in a stage-specific manner. We also found that the dynamic expression pattern of the BMP antagonist Gremlin1 (Grem1) in limb buds, which is critical for autopodial expansion, was consistent with the digital pattern of emu, chicken and zebra finch. Furthermore, in emu, variation among individuals was observed in the width of Grem1 expression in forelimb buds, as well as in the adult skeletal pattern. Our results support the view that the signalling system that regulates the dynamic expression of Grem1 in the limb bud contributes substantially to variations in avian digital patterns.

'Building a perfect body': control of vertebrate organogenesis by PBX-dependent regulatory networks

Pbx genes encode transcription factors that belong to the TALE (three-amino-acid loop extension) superclass of homeodomain proteins. We have witnessed a surge in information about the roles of this gene family as leading actors in the transcriptional control of development. PBX proteins represent a clear example of how transcription factors can regulate developmental processes by combinatorial properties, acting within multimeric complexes to implement activation or repression of transcription depending on their interaction partners. Here, we revisit long-emphasized functions of PBX transcription factors as cofactors for HOX proteins, major architects of the body plan. We further discuss new knowledge on roles of PBX proteins in different developmental contexts as upstream regulators of Hox genes-as factors that interact with non-HOX proteins and can work independently of HOX-as well as potential pioneer factors. Committed to building a perfect body, PBX proteins govern regulatory networks that direct essential morphogenetic processes and organogenesis in vertebrate development. Perturbations of PBX-dependent networks can cause human congenital disease and cancer.

Fibroblast Growth Factor 10 and Vertebrate Limb Development

Early limb development requires fibroblast growth factor (Fgf)-mediated coordination between growth and patterning to ensure the proper formation of a functional organ. The apical ectodermal ridge (AER) is a domain of thickened epithelium located at the distal edge of the limb bud that coordinates outgrowth along the proximodistal axis. Considerable amount of work has been done to elucidate the cellular and molecular mechanisms underlying induction, maintenance and regression of the AER. Fgf10, a paracrine Fgf that elicits its biological responses by activating the fibroblast growth factor receptor 2b (Fgfr2b), is crucial for governing proximal distal outgrowth as well as patterning and acts upstream of the known AER marker Fgf8. A transgenic mouse line allowing doxycycline-based inducible and ubiquitous expression of a soluble form of Fgfr2b has been extensively used to identify the role of Fgfr2b ligands at different time points during development. Overexpression of soluble Fgfr2b (sFgfr2b) post-AER induction leads to irreversible loss of cellular β-catenin organization and decreased Fgf8 expression in the AER. A similar approach has been carried out pre-AER induction. The observed limb phenotype is similar to the severe proximal truncations observed in human babies exposed to thalidomide, which has been proposed to block the Fgf10-AER-Fgf8 feedback loop. Novel insights on the role of Fgf10 signaling in limb formation pre- and post-AER induction are summarized in this review and will be integrated with possible future investigations on the role of Fgf10 throughout limb development.

Unique pelvic fin in a tetrapod-like fossil fish, and the evolution of limb patterning

The fossil fish Rhizodus hibberti, a member of the tetrapod stem group, shows a unique skeletal pattern in the pelvic fin. Rather than the highly conserved one-to-two pattern of a femur, tibia, and fibula (seen in all known tetrapods, including the extinct, fishlike members of the group), the fin of Rhizodus comprises a femur articulating distally with three bones, each with a distinct morphology. This reveals an early stage in the evolution of limb development, in which the processes patterning the proximal parts of the embryonic fin/limb (the stylopod and zeugopod) were not constrained in the way seen in living tetrapods and could produce more varied skeletal patterns in the adult.

All living tetrapods have a one-to-two branching pattern in the embryonic proximal limb skeleton, with a single element at the base of the limb (the humerus or femur) that articulates distally with two parallel radials (the ulna and radius or the tibia and fibula). This pattern is also seen in the fossilized remains of stem-tetrapods, including the fishlike members of the group, in which despite the absence of digits, the proximal parts of the fin skeleton clearly resemble those of later tetrapods. However, little is known about the developmental mechanisms that establish and canalize this highly conserved pattern. We describe the well-preserved pelvic fin skeleton of Rhizodus hibberti, a Carboniferous sarcopterygian (lobe-finned) fish, and member of the tetrapod stem group. In this specimen, three parallel radials, each robust with a distinct morphology, articulate with the femur. We review this unexpected morphology in a phylogenetic and developmental context. It implies that the developmental patterning mechanisms seen in living tetrapods, now highly constrained, evolved from mechanisms flexible enough to accommodate variation in the zeugopod (even between pectoral and pelvic fins), while also allowing each element to have a unique morphology.

Functional Characterization of the Osteoarthritis Genetic Risk Residing at ALDH1A2 Identifies rs12915901 as a Key Target Variant

OBJECTIVE

To identify the functional single-nucleotide polymorphisms (SNPs) and mechanisms conferring increased risk of hand osteoarthritis (OA) at the ALDH1A2 locus, which is a retinoic acid regulatory gene.

METHODS

Tissue samples from 247 patients with knee, hip, or hand OA who had undergone joint surgery were included. RNA-sequencing analysis was used to investigate differential expression of ALDH1A2 and other retinoic acid signaling pathway genes in cartilage. Expression of ALDH1A2 in joint tissues obtained from multiple sites was quantified using quantitative reverse transcription-polymerase chain reaction. Allelic expression imbalance (AEI) was measured by pyrosequencing. The consequences of ALDH1A2 depletion by RNA interference were assessed in primary human chondrocytes. In silico and in vitro analyses were used to pinpoint which, among 62 highly correlated SNPs, could account for the association at the locus.

RESULTS

ALDH1A2 expression was observed across multiple joint tissue samples, including osteochondral tissue from the hand. The expression of ALDH1A2 and of several retinoic acid signaling genes was different in diseased cartilage compared to non-diseased cartilage, with ALDH1A2 showing lower levels in OA cartilage. Experimental depletion of ALDH1A2 resulted in changes in the expression levels of a number of chondrogenic markers, including SOX9. In addition, reduced expression of the OA risk-conferring allele was witnessed in a number of joint tissues, with the strongest effect in cartilage. The intronic SNP rs12915901 recapitulated the AEI observed in patient tissues, while the Ets transcription factors were identified as potential mediators of this effect.

CONCLUSION

The ALDH1A2 locus seems to increase the risk of hand OA through decreased expression of ALDH1A2 in joint tissues, with the effect dependent on rs12915901. These findings indicate a mechanism that may now be targeted to modulate OA risk.

Expression of meis and hoxa11 in dipnoan and teleost fins provides new insights into the evolution of vertebrate appendages

Background

The concerted activity of Meis and Hoxa11 transcription factors is essential for the subdivision of tetrapod limbs into proximo-distal (PD) domains; however, little is know about the evolution of this patterning mechanism. Here, we aim to study the expression of meis and hoxa11 orthologues in the median and paired rayed fins of zebrafish and in the lobed fins of the Australian lungfish.

Results

First, a late phase of expression of meis1.1 and hoxa11b in zebrafish dorsal and anal fins relates with segmentation of endochondral elements in proximal and distal radials. Second, our zebrafish in situ hybridization results reveal spatial and temporal changes between pectoral and pelvic fins. Third, in situ analysis of meis1, meis3 and hoxa11 genes in Neoceratodus pectoral fins identifies decoupled domains of expression along the PD axis.

Conclusions

Our data raise the possibility that the origin of stylopod and zeugopod lies much deeper in gnathostome evolution and that variation in meis and hoxa11 expression has played a substantial role in the transformation of appendage anatomy. Moreover, these observations provide evidence that the Meis/Hoxa11 profile considered a hallmark of stylopod/zeugopod patterning is present in Neoceratodus.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0099-9) contains supplementary material, which is available to authorized users.

Glycol chitosan functionalized asenapine nanostructured lipid carriers for targeted brain delivery: Pharmacokinetic and teratogenic assessment

Blood brain barrier (BBB) is a complex, tight barrier between endothelial cells of cerebral blood vessels. It acts as a physical barrier and provides access to only those moieties which are necessary for proper brain functioning. However, this selective prudence also acts as a hindrance in therapeutic targeting of brain necessitating pharmaceutical intervention. Intranasal drug delivery is one such approach which we have exploited here for targeted brain delivery of asenapine by glycol chitosan coated nanostructured lipid carrier (GC-ANLC). The best formulation was characterized for particle size (184.2±5.59nm), zeta potential (18.83±1.18mV), entrapment efficiency (83.52±2.59%) and surface morphology (spherical and smooth). In-vitro drug-release study showed that Higuchi model (r²=0.9938, AIC=52.94) dictated asenapine release from GC-ANLC. Cell compatibility study suggested biocompatibility of GC-ANLC with A549 cell line as well as nasal epithelial cell membrane. After intranasal delivery, Charles-Foster rats demonstrated approximately 2.3 and 4 fold higher systemic and brain bioavailability of GC-ANLC compared to asenapine solution (ASM). Embryo fetal toxicity study was further conducted to investigate the teratogenic effect of GC-ANLC. In conclusion, prepared GC-ANLC could be used as a promising drug carrier for delivery of asenapine via intranasal route with better pharmacokinetic and safety profile.

Transcriptomic insights into the genetic basis of mammalian limb diversity

Background

From bat wings to whale flippers, limb diversification has been crucial to the evolutionary success of mammals. We performed the first transcriptome-wide study of limb development in multiple species to explore the hypothesis that mammalian limb diversification has proceeded through the differential expression of conserved shared genes, rather than by major changes to limb patterning. Specifically, we investigated the manner in which the expression of shared genes has evolved within and among mammalian species.

Results

We assembled and compared transcriptomes of bat, mouse, opossum, and pig fore- and hind limbs at the ridge, bud, and paddle stages of development. Results suggest that gene expression patterns exhibit larger variation among species during later than earlier stages of limb development, while within species results are more mixed. Consistent with the former, results also suggest that genes expressed at later developmental stages tend to have a younger evolutionary age than genes expressed at earlier stages. A suite of key limb-patterning genes was identified as being differentially expressed among the homologous limbs of all species. However, only a small subset of shared genes is differentially expressed in the fore- and hind limbs of all examined species. Similarly, a small subset of shared genes is differentially expressed within the fore- and hind limb of a single species and among the forelimbs of different species.

Conclusions

Taken together, results of this study do not support the existence of a phylotypic period of limb development ending at chondrogenesis, but do support the hypothesis that the hierarchical nature of development translates into increasing variation among species as development progresses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0902-6) contains supplementary material, which is available to authorized users.

Cellular and molecular drivers of differential organ growth: insights from the limbs of Monodelphis domestica

A fundamental question in biology is "how is growth differentially regulated during development to produce organs of particular sizes?" We used a new model system for the study of differential organ growth, the limbs of the opossum (Monodelphis domestica), to investigate the cellular and molecular basis of differential organ growth in mammals. Opossum forelimbs grow much faster than hindlimbs, making opossum limbs an exceptional system with which to study differential growth. We first used the great differences in opossum forelimb and hindlimb growth to identify cellular processes and molecular signals that underlie differential limb growth. We then used organ culture and pharmacological addition of FGF ligands and inhibitors to test the role of the Fgf/Mitogen-activated protein kinases (MAPK) signaling pathway in driving these cellular processes. We found that molecular signals from within the limb drive differences in cell proliferation that contribute to the differential growth of the forelimb and hindlimbs of opossums. We also found that alterations in the Fgf/MAPK pathway can generate differences in cell proliferation that mirror those observed between wild-type forelimb and hindlimbs of opossums and that manipulation of Fgf/MAPK signaling affects downstream focal adhesion-extracellular matrix (FA-ECM) and Wnt signaling in opossum limbs. Taken together, these findings suggest that evolutionary changes in the Fgf/MAPK pathway could help drive the observed differences in cell behaviors and growth in opossum forelimb and hindlimbs.

Functional and Comparative Genomics of Hoxa2 Gene cis-Regulatory Elements: Evidence for Evolutionary Modification of Ancestral Core Element Activity

Hoxa2 is an evolutionarily conserved developmental regulatory gene that functions to specify rhombomere (r) and pharyngeal arch (PA) identities throughout the Osteichthyes. Japanese medaka (Oryzias latipes) hoxa2a, like orthologous Hoxa2 genes from other osteichthyans, is expressed during embryogenesis in r2–7 and PA2-7, whereas the paralogous medaka pseudogene, ψhoxa2b, is expressed in noncanonical Hoxa2 domains, including the pectoral fin buds. To understand the evolution of cis-regulatory element (CRE) control of gene expression, we conducted eGFP reporter gene expression studies with extensive functional mapping of several conserved CREs upstream of medaka hoxa2a and ψhoxa2b in transient and stable-line transgenic medaka embryos. The CREs tested were previously shown to contribute to directing mouse Hoxa2 gene expression in r3, r5, and PA2-4. Our results reveal the presence of sequence elements embedded in the medaka hoxa2a and ψhoxa2b upstream enhancer regions (UERs) that mediate expression in r4 and the PAs (hoxa2a r4/CNCC element) or in r3–7 and the PAs ψhoxa2b r3–7/CNCC element), respectively. Further, these elements were shown to be highly conserved among osteichthyans, which suggests that the r4 specifying element embedded in the UER of Hoxa2 is a deeply rooted rhombomere specifying element and the activity of this element has been modified by the evolution of flanking sequences that redirect its activity to alternative developmental compartments.

Geometric Morphometrics on Gene Expression Patterns Within Phenotypes: A Case Example on Limb Development

How the genotype translates into the phenotype through development is critical to fully understand the evolution of phenotypes. We propose a novel approach to directly assess how changes in gene expression patterns are associated with changes in morphology using the limb as a case example. Our method combines molecular biology techniques, such as whole-mount in situ hybridization, with image and shape analysis, extending the use of Geometric Morphometrics to the analysis of nonanatomical shapes, such as gene expression domains. Elliptical Fourier and Procrustes-based semilandmark analyses were used to analyze the variation and covariation patterns of the limb bud shape with the expression patterns of two relevant genes for limb morphogenesis, Hoxa11 and Hoxa13. We devised a multiple thresholding method to semiautomatically segment gene domains at several expression levels in large samples of limb buds from C57Bl6 mouse embryos between 10 and 12 postfertilization days. Besides providing an accurate phenotyping tool to quantify the spatiotemporal dynamics of gene expression patterns within developing structures, our morphometric analyses revealed high, non-random, and gene-specific variation undergoing canalization during limb development. Our results demonstrate that Hoxa11 and Hoxa13, despite being paralogs with analogous functions in limb patterning, show clearly distinct dynamic patterns, both in shape and size, and are associated differently with the limb bud shape. The correspondence between our results and already well-established molecular processes underlying limb development confirms that this morphometric approach is a powerful tool to extract features of development regulating morphogenesis. Such multilevel analyses are promising in systems where not so much molecular information is available and will advance our understanding of the genotype–phenotype map. In systematics, this knowledge will increase our ability to infer how evolution modified a common developmental pattern to generate a wide diversity of morphologies, as in the vertebrate limb.

The Relationship between Gene Network Structure and Expression Variation among Individuals and Species

Variation among individuals is a prerequisite of evolution by natural selection. As such, identifying the origins of variation is a fundamental goal of biology. We investigated the link between gene interactions and variation in gene expression among individuals and species using the mammalian limb as a model system. We first built interaction networks for key genes regulating early (outgrowth; E9.5-11) and late (expansion and elongation; E11-13) limb development in mouse. This resulted in an Early (ESN) and Late (LSN) Stage Network. Computational perturbations of these networks suggest that the ESN is more robust. We then quantified levels of the same key genes among mouse individuals and found that they vary less at earlier limb stages and that variation in gene expression is heritable. Finally, we quantified variation in gene expression levels among four mammals with divergent limbs (bat, opossum, mouse and pig) and found that levels vary less among species at earlier limb stages. We also found that variation in gene expression levels among individuals and species are correlated for earlier and later limb development. In conclusion, results are consistent with the robustness of the ESN buffering among-individual variation in gene expression levels early in mammalian limb development, and constraining the evolution of early limb development among mammalian species.

Thalidomide reduces mechanical hyperalgesia and depressive-like behavior induced by peripheral nerve crush in mice

BACKGROUND

It has been shown that chronic pain is able to induce depressive disorders in humans, in part, due to peripheral inflammation that reaches the central nervous system. However, the mechanisms involved remain to be established. The purpose of this study was to investigate whether sciatic nerve crush could produce depression-like behaviors, in addition to pain-related behaviors, in mice. Once confirmed, this model was used to investigate tumor necrosis factor-α (TNF-α) as a key mediator involved in the pathophysiology of both pain and depression.

EXPERIMENTAL APPROACH

Male Swiss mice were divided into three groups, naïve, sham and operated. In the operated group, the sciatic nerve was crushed. Following surgery, animals from the operated group were treated daily by oral gavage (p.o.) with saline (10 ml/kg), fluoxetine (20 mg/kg) or thalidomide (10 mg/kg) for 15 days. Mechanical hyperalgesia was evaluated every 3 days by von Frey filaments and depressive-like behavior was assessed at the end of day 15, using the tail suspension test (TST) and the forced swimming test (FST). Then, samples from the prefrontal cortex, hippocampus and sciatic nerve were processed to measure TNF-α levels by enzyme-linked immunosorbent assay (ELISA).

RESULTS

Crush caused significant mechanical hyperalgesia and depressive-like behaviors and increased TNF-α levels in the sciatic nerve, prefrontal cortex and hippocampus of operated animals. Treatment with fluoxetine or thalidomide reversed crush-induced mechanical hyperalgesia, depressive-like behaviors and the increased TNF-α levels in the sciatic nerve, prefrontal cortex and hippocampus.

CONCLUSIONS

The sciatic nerve crush model represents a good model to study to mechanisms underlying both pain and depressive-like behaviors. Furthermore, inhibitors of TNF-α synthesis, like thalidomide, have a potential to treat depressive disorders associated with neuropathic pain.

T396I mutation of mouse Sufu reduces the stability and activity of Gli3 repressor

Hedgehog signaling is primarily transduced by two transcription factors: Gli2, which mainly acts as a full-length activator, and Gli3, which tends to be proteolytically processed from a full-length form (Gli3^FL) to an N-terminal repressor (Gli3^REP). Recent studies using a Sufu knockout mouse have indicated that Sufu is involved in regulating Gli2 and Gli3 activator and repressor activity at multiple steps of the signaling cascade; however, the mechanism of specific Gli2 and Gli3 regulation remains to be elucidated. In this study, we established an allelic series of ENU-induced mouse strains. Analysis of one of the missense alleles, Sufu^T396I, showed that Thr³⁹⁶ residue of Sufu played a key role in regulation of Gli3 activity. Sufu^T396I/T396I embryos exhibited severe polydactyly, which is indicative of compromised Gli3 activity. Concomitantly, significant quantitative reductions of unprocessed Gli3 (Gli3^FL) and processed Gli3 (Gli3^REP) were observed in vivo as well as in vitro. Genetic experiments showed that patterning defects in the limb buds of Sufu^T396I/T396I were rescued by a constitutive Gli3^REP allele (Gli3^∆699), strongly suggesting that Sufu^T396I reduced the truncated Gli3 repressor. In contrast, Sufu^T396I qualitatively exhibited no mutational effects on Gli2 regulation. Taken together, the results of this study show that the Thr³⁹⁶ residue of Sufu is specifically required for regulation of Gli3 but not Gli2. This implies a novel Sufu-mediated mechanism in which Gli2 activator and Gli3 repressor are differentially regulated.

ESCRT-II/Vps25 constrains digit number by endosome-mediated selective modulation of FGF-SHH signaling

Sorting and degradation of receptors and associated signaling molecules maintain homeostasis of conserved signaling pathways during cell specification and tissue development. Yet, whether machineries that sort signaling proteins act preferentially on different receptors and ligands in different contexts remains mysterious. Here, we show that Vacuolar protein sorting 25, Vps25, a component of ESCRT-II (Endosomal Sorting Complex Required for Transport II), directs preferential endosome-mediated modulation of FGF signaling in limbs. By ENU-induced mutagenesis, we isolated a polydactylous mouse line carrying a hypomorphic mutation of Vps25 (Vps25(ENU)). Unlike Vps25-null embryos we generated, Vps25(ENU/ENU) mutants survive until late gestation. Their limbs display FGF signaling enhancement and consequent hyperactivation of the FGF-SHH feedback loop causing polydactyly, whereas WNT and BMP signaling remain unperturbed. Notably, Vps25(ENU/ENU) Mouse Embryonic Fibroblasts exhibit aberrant FGFR trafficking and degradation; however, SHH signaling is unperturbed. These studies establish that the ESCRT-II machinery selectively limits FGF signaling in vertebrate skeletal patterning.

Direct functional consequences of ZRS enhancer mutation combine with secondary long range SHH signalling effects to cause preaxial polydactyly

Sonic hedgehog (SHH) plays a central role in patterning numerous embryonic tissues including, classically, the developing limb bud where it controls digit number and identity. This study utilises the polydactylous Silkie (Slk) chicken breed, which carries a mutation in the long range limb-specific regulatory element of SHH, the ZRS. Using allele specific SHH expression analysis combined with quantitative protein analysis, we measure allele specific changes in SHH mRNA and concentration of SHH protein over time. This confirms that the Slk ZRS enhancer mutation causes increased SHH expression in the posterior leg mesenchyme. Secondary consequences of this increased SHH signalling include increased FGF pathway signalling and growth as predicted by the SHH/GREM1/FGF feedback loop and the Growth/Morphogen models. Manipulation of Hedgehog, FGF signalling and growth demonstrate that anterior-ectopic expression of SHH and induction of preaxial polydactyly is induced secondary to increased SHH signalling and Hedgehog-dependent growth directed from the posterior limb. We predict that increased long range SHH signalling acts in combination with changes in activation of SHH transcription from the Slk ZRS allele. Through analysis of the temporal dynamics of anterior SHH induction we predict a gene regulatory network which may contribute to activation of anterior SHH expression from the Slk ZRS.

•

Overexpression of posterior SHH in the limb bud can cause preaxial polydactyly.

•

Increased activation of SHH/GREM/FGF feedback and growth induces Slk preaxial polydactyly.

•

Autoregulated expression of SHH can occur within 1.5–2 h in the limb bud.

DLX5, FGF8 and the Pin1 isomerase control ΔNp63α protein stability during limb development: a regulatory loop at the basis of the SHFM and EEC congenital malformations

Ectrodactyly, or Split-Hand/Foot Malformation (SHFM), is a congenital condition characterized by the loss of central rays of hands and feet. The p63 and the DLX5;DLX6 transcription factors, expressed in the embryonic limb buds and ectoderm, are disease genes for these conditions. Mutations of p63 also cause the ectodermal dysplasia–ectrodactyly–cleft lip/palate (EEC) syndrome, comprising SHFM. Ectrodactyly is linked to defects of the apical ectodermal ridge (AER) of the developing limb buds. FGF8 is the key signaling molecule in this process, able to direct proximo-distal growth and patterning of the skeletal primordial of the limbs. In the limb buds of both p63 and Dlx5;Dlx6 murine models of SHFM, the AER is poorly stratified and FGF8 expression is severely reduced. We show here that the FGF8 locus is a downstream target of DLX5 and that FGF8 counteracts Pin1–ΔNp63α interaction. In vivo, lack of Pin1 leads to accumulation of the p63 protein in the embryonic limbs and ectoderm. We show also that ΔNp63α protein stability is negatively regulated by the interaction with the prolyl-isomerase Pin1, via proteasome-mediated degradation; p63 mutant proteins associated with SHFM or EEC syndromes are resistant to Pin1 action. Thus, DLX5, p63, Pin1 and FGF8 participate to the same time- and location-restricted regulatory loop essential for AER stratification, hence for normal patterning and skeletal morphogenesis of the limb buds. These results shed new light on the molecular mechanisms at the basis of the SHFM and EEC limb malformations.

Molecular analysis of regulative events in the developing chick limb

The developing chick limb has the remarkable ability to regulate for the loss of large amounts of mesenchyme and maintain a normal limb pattern in early (Hamburger and Hamilton Stage 19; E3) limbs. How the limb can regulate for tissue loss and why this ability is lost as development proceeds (after Hamburger and Hamilton Stage 21; E3.5) is unclear. We have investigated the origins of cells involved in regulative processes and, for the first time, the molecular changes occurring, and find striking differences between developmental time points just 0.5 days apart. We demonstrate that subtle changes in cell dispersal and cell proliferation occur in HH St21 limbs but not in HH St19 limbs and also demonstrate that there is no net replacement of removed tissue at either HH St21 or St19. We further show that changes in the Fgf8/Shh/Bmp4/Gremlin signaling pathway together with the appearance of distal Hox gene activation coincide with the limbs' ability to regulate for large tissue loss. We also demonstrate that following small tissue loss, limbs can regulate for missing tissue to produce normal pattern with no net replacement of missing tissue, as seen in limbs following large tissue loss. Our results indicate the regulative ability of the limb is not due to changes in cell proliferation, cell lineage nor replacement of the missing tissue - regulative ability is reliant upon the signaling environment remaining.

The origin of the tetrapod limb: from expeditions to enhancers

More than three centuries ago natural philosophers, and later anatomists, recognized a fundamental organization to the skeleton of tetrapod limbs. Composed of three segments, stylopod, zeugopod, and autopod, this pattern has served as the basis for a remarkably broad adaptive radiation from wings and flippers to hands and digging organs. A central area of inquiry has been tracing the origins of the elements of this Bauplan in the fins of diverse fish. Can equivalents of the three segments, and the developmental processes that pattern them, be seen in fish fins? In addition, if so, how do these data inform theories of the transformation of fins into limbs? Answers to these questions come from linking discoveries in paleontology with those of developmental biology and genetics. Burgeoning discoveries in the regulatory biology of developmental genes and in the genomics of diverse species offer novel data to investigate these classical questions.

BMP-mediated functional cooperation between Dlx5;Dlx6 and Msx1;Msx2 during mammalian limb development

The Dlx and Msx homeodomain transcription factors play important roles in the control of limb development. The combined disruption of Msx1 and Msx2, as well as that of Dlx5 and Dlx6, lead to limb patterning defects with anomalies in digit number and shape. Msx1;Msx2 double mutants are characterized by the loss of derivatives of the anterior limb mesoderm which is not observed in either of the simple mutants. Dlx5;Dlx6 double mutants exhibit hindlimb ectrodactyly. While the morphogenetic action of Msx genes seems to involve the BMP molecules, the mode of action of Dlx genes still remains elusive. Here, examining the limb phenotypes of combined Dlx and Msx mutants we reveal a new Dlx-Msx regulatory loop directly involving BMPs. In Msx1;Dlx5;Dlx6 triple mutant mice (TKO), beside the expected ectrodactyly, we also observe the hallmark morphological anomalies of Msx1;Msx2 double mutants suggesting an epistatic role of Dlx5 and Dlx6 over Msx2. In Msx2;Dlx5;Dlx6 TKO mice we only observe an aggravation of the ectrodactyly defect without changes in the number of the individual components of the limb. Using a combination of qPCR, ChIP and bioinformatic analyses, we identify two Dlx/Msx regulatory pathways: 1) in the anterior limb mesoderm a non-cell autonomous Msx-Dlx regulatory loop involves BMP molecules through the AER and 2) in AER cells and, at later stages, in the limb mesoderm the regulation of Msx2 by Dlx5 and Dlx6 occurs also cell autonomously. These data bring new elements to decipher the complex AER-mesoderm dialogue that takes place during limb development and provide clues to understanding the etiology of congenital limb malformations.

Digit patterning during limb development as a result of the BMP-receptor interaction

Turing models have been proposed to explain the emergence of digits during limb development. However, so far the molecular components that would give rise to Turing patterns are elusive. We have recently shown that a particular type of receptor-ligand interaction can give rise to Schnakenberg-type Turing patterns, which reproduce patterning during lung and kidney branching morphogenesis. Recent knockout experiments have identified Smad4 as a key protein in digit patterning. We show here that the BMP-receptor interaction meets the conditions for a Schnakenberg-type Turing pattern, and that the resulting model reproduces available wildtype and mutant data on the expression patterns of BMP, its receptor, and Fgfs in the apical ectodermal ridge (AER) when solved on a realistic 2D domain that we extracted from limb bud images of E11.5 mouse embryos. We propose that receptor-ligand-based mechanisms serve as a molecular basis for the emergence of Turing patterns in many developing tissues.

Conserved cis-regulatory regions in a large genomic landscape control SHH and BMP-regulated Gremlin1 expression in mouse limb buds

Background

Mouse limb bud is a prime model to study the regulatory interactions that control vertebrate organogenesis. Major aspects of limb bud development are controlled by feedback loops that define a self-regulatory signalling system. The SHH/GREM1/AER-FGF feedback loop forms the core of this signalling system that operates between the posterior mesenchymal organiser and the ectodermal signalling centre. The BMP antagonist Gremlin1 (GREM1) is a critical node in this system, whose dynamic expression is controlled by BMP, SHH, and FGF signalling and key to normal progression of limb bud development. Previous analysis identified a distant cis-regulatory landscape within the neighbouring Formin1 (Fmn1) locus that is required for Grem1 expression, reminiscent of the genomic landscapes controlling HoxD and Shh expression in limb buds.

Results

Three highly conserved regions (HMCO1-3) were identified within the previously defined critical genomic region and tested for their ability to regulate Grem1 expression in mouse limb buds. Using a combination of BAC and conventional transgenic approaches, a 9 kb region located ~70 kb downstream of the Grem1 transcription unit was identified. This region, termed Grem1 Regulatory Sequence 1 (GRS1), is able to recapitulate major aspects of Grem1 expression, as it drives expression of a LacZ reporter into the posterior and, to a lesser extent, in the distal-anterior mesenchyme. Crossing the GRS1 transgene into embryos with alterations in the SHH and BMP pathways established that GRS1 depends on SHH and is modulated by BMP signalling, i.e. integrates inputs from these pathways. Chromatin immunoprecipitation revealed interaction of endogenous GLI3 proteins with the core cis-regulatory elements in the GRS1 region. As GLI3 is a mediator of SHH signal transduction, these results indicated that SHH directly controls Grem1 expression through the GRS1 region. Finally, all cis-regulatory regions within the Grem1 genomic landscape locate to the DNAse I hypersensitive sites identified in this genomic region by the ENCODE consortium.

Conclusions

This study establishes that distant cis-regulatory regions scattered through a larger genomic landscape control the highly dynamic expression of Grem1, which is key to normal progression of mouse limb bud development.

The molecular basis of human congenital limb malformations

This review focuses predominantly on the human congenital malformations caused by alterations affecting the morphoregulatory gene networks that control early limb bud patterning and outgrowth. Limb defects are among the most frequent congenital malformations in humans that are caused by genetic mutations or teratogenic effects resulting either in abnormal, loss of, or additional skeletal elements. Spontaneous and engineered mouse models have been used to identify and study the molecular alterations and disrupted gene networks that underlie human congenital limb malformations. More recently, mouse genetics has begun to reveal the alterations that affect the often-large cis-regulatory landscapes that control gene expression in limb buds and cause devastating effects on limb bud development. These findings have paved the way to identifying mutations in cis-regulatory regions as causal to an increasing number of congenital limb malformations in humans. In these cases, no mutations in the coding region of a presumed candidate were previously detected. This review highlights how the current understanding of the molecular gene networks and interactions that control mouse limb bud development provides insight into the etiology of human congenital limb malformations.

Hand asymmetry does not relate to key life history traits in post-menopausal contemporary Finnish women

Associations between fluctuating asymmetry (FA, a putative marker of developmental instability, DI) and life history traits have received a great deal of attention in the non-human literature. However, the patterns found are very heterogeneous and generalizations are difficult to make. In humans, only a few studies have related FA to life histories and fitness. In this paper we study such relationships using hand FA and several key life history traits in 209 post-menopausal Finnish women born between 1946 and 1958. Asymmetry measurements were based on scans of the hands and the life histories of these women were collected using questionnaires. No significant associations were detected and trends were opposite to expectations. We did find evidence for directional asymmetry, as traits in the right hand were larger on average. This may be due to handedness, questioning the usefulness of hand FA as a measure of DI. We conclude that future studies in humans should carefully examine the usefulness of traits as measures of DI.

The making of differences between fins and limbs

'Evo-devo', an interdisciplinary field based on developmental biology, includes studies on the evolutionary processes leading to organ morphologies and functions. One fascinating theme in evo-devo is how fish fins evolved into tetrapod limbs. Studies by many scientists, including geneticists, mathematical biologists, and paleontologists, have led to the idea that fins and limbs are homologous organs; now it is the job of developmental biologists to integrate these data into a reliable scenario for the mechanism of fin-to-limb evolution. Here, we describe the fin-to-limb transition based on key recent developmental studies from various research fields that describe mechanisms that may underlie the development of fins, limb-like fins, and limbs.

A function for all posterior Hoxd genes during digit development?

BACKGROUND

Four posterior Hoxd genes, from Hoxd13 to Hoxd10, are collectively regulated during the development of tetrapod digits. Besides the well-documented role of Hoxd13, the function of the neighboring genes has been difficult to evaluate due to the close genetic linkage and potential regulatory interferences. We used a combination of five small nested deletions in cis, involving from two to four consecutive genes of the Hoxd13 to Hoxd9 loci, in mice, to evaluate their combined functional importance.

RESULTS

We show that deletions leading to a gain of function of Hoxd13, via regulatory re-allocation, generate abnormal phenotypes, in agreement with the dominant negative role of this gene. We also show that Hoxd10, Hoxd11, and Hoxd12 all seem to play a genuine role in digit development, though less compelling than that of Hoxd13. In contrast, the nearby Hoxd9 contributed no measurable function in digits.

CONCLUSIONS

We conclude that a slight and transient deregulation of Hoxd13 expression can readily affect the relative lengths of limb segments and that all posterior Hoxd genes likely contribute to the final limb morphology. We discuss the difficulty to clearly assess the functional share of individual genes within such a gene family, where closely located neighbors, coding for homologous proteins, are regulated by a unique circuitry and all contribute to shape the distal parts of our appendages.

An N-terminal G11A mutation in HOXD13 causes synpolydactyly and interferes with Gli3R function during limb pre-patterning

Synpolydactyly (SPD) is a distal limb anomaly characterized by incomplete digit separation and the presence of supernumerary digits in the syndactylous web. This phenotype has been associated with mutations in the homeodomain or polyalanine tract of the HOXD13 gene. We identified a novel mutation (G11A) in HOXD13 that is located outside the previously known domains and affects the intracellular half life of the protein. Misexpression of HOXD13(G11A) in the developing chick limb phenocopied the human SPD phenotype. Finally, we demonstrated through in vitro studies that this mutation has a destabilizing effect on GLI3R uncovering an unappreciated mechanism by which HOXD13 determines the patterning of the limb.

Limb regeneration: a new development?

Salamander limb regeneration is a classical model of tissue morphogenesis and patterning. Through recent advances in cell labeling and molecular analysis, a more precise, mechanistic understanding of this process has started to emerge. Long-standing questions include to what extent limb regeneration recapitulates the events observed in mammalian limb development and to what extent are adult- or salamander- specific aspects deployed. Historically, researchers studying limb development and limb regeneration have proposed different models of pattern formation. Here we discuss recent data on limb regeneration and limb development to argue that although patterning mechanisms are likely to be similar, cell plasticity and signaling from nerves play regeneration-specific roles.

A symphony of regulations centered on p63 to control development of ectoderm-derived structures

The p53-related transcription factor p63 is critically important for basic cellular functions during development of the ectoderm and derived structure and tissues, including skin, limb, palate, and hair. On the one side, p63 is required to sustain the proliferation of keratinocyte progenitors, while on the other side it is required for cell stratification, commitment to differentiate, cell adhesion, and epithelial-mesenchymal signaling. Molecules that are components or regulators of the p63 pathway(s) are rapidly being identified, and it comes with no surprise that alterations in the p63 pathway lead to congenital conditions in which the skin and other ectoderm-derived structures are affected. In this paper, we summarize the current knowledge of the molecular and cellular regulations centered on p63, derived from the comprehension of p63-linked human diseases and the corresponding animal models, as well as from cellular models and high-throughput molecular approaches. We point out common themes and features, that allow to speculate on the possible role of p63 downstream events and their potential exploitation in future attempts to correct the congenital defect in preclinical studies.

Gene set assembly for quantitative prediction of developmental toxicity in the embryonic stem cell test

The embryonic stem cell test (EST) is an in vitro method for predicting developmental toxicity based on compound-induced inhibition of embryonic stem cell (ESC) differentiation. We previously described how gene expression analysis in the EST can be used to describe normal ESC differentiation as well as identify compound developmental toxicity, by means of our differentiation track algorithm. In this study, we combined raw data from our three previous studies in a new integrated analysis, to identify a gene set that allows for improved prediction. By evaluating predictions of 100,000 randomly selected gene sets, we identified which genes contribute significantly to the prediction reliability. By additional cross-validation, we identified a set of 52 genes that allows for improved prediction of toxicity. The correlation between the predictions using this gene set and the magnitude of the EST endpoint was 0.85, and the gene set predicted developmental toxicity with 83% accuracy (area under the curve 89%). If compounds with ineffective data because of a too low tested concentration or too much variation between samples were excluded, even 100% accuracy could be reached based on 15 compounds. This novel gene set consists mainly of genes involved in the stem cell differentiation or other developmental processes. We expect that this set can be of use in future studies aimed at improving the EST for risk assessment, thus making a next step towards regulatory implementation of this method.

Pbx homeodomain proteins: TALEnted regulators of limb patterning and outgrowth

Limb development has long provided an excellent model for understanding the genetic principles driving embryogenesis. Studies utilizing chick and mouse have led to new insights into limb patterning and morphogenesis. Recent research has centered on the regulatory networks underlying limb development. Here, we discuss the hierarchical, overlapping, and iterative roles of Pbx family members in appendicular development that have emerged from genetic analyses in the mouse. Pbx genes are essential in determining limb bud positioning, early bud formation, limb axes establishment and coordination, and patterning and morphogenesis of most elements of the limb and girdle. Pbx proteins directly regulate critical effectors of limb and girdle development, including morphogen-encoding genes like Shh in limb posterior mesoderm, and transcription factor-encoding genes like Alx1 in pre-scapular domains. Interestingly, at least in limb buds, Pbx appear to act not only as Hox cofactors, but also in the upstream control of 5' HoxA/D gene expression.