Alx1 Deficient Mice Recapitulate Craniofacial Phenotype and Reveal Developmental Basis of ALX1-Related Frontonasal Dysplasia

Identification and characterization of intermediate states in mammalian neural crest cell epithelial to mesenchymal transition and delamination

Epithelial to mesenchymal transition (EMT) is a cellular process that converts epithelial cells to mesenchymal cells with migratory potential in developmental and pathological processes. Although originally considered a binary event, EMT in cancer progression involves intermediate states between a fully epithelial and a fully mesenchymal phenotype, which are characterized by distinct combinations of epithelial and mesenchymal markers. This phenomenon has been termed epithelial to mesenchymal plasticity (EMP), however, the intermediate states remain poorly described and it's unclear whether they exist during developmental EMT. Neural crest cells (NCC) are an embryonic progenitor cell population that gives rise to numerous cell types and tissues in vertebrates, and their formation and delamination is a classic example of developmental EMT. However, whether intermediate states also exist during NCC EMT and delamination remains unknown. Through single-cell RNA sequencing of mouse embryos, we identified intermediate NCC states based on their transcriptional signature and then spatially defined their locations in situ in the dorsolateral neuroepithelium. Our results illustrate the importance of cell cycle regulation and functional role for the intermediate stage marker Dlc1 in facilitating mammalian cranial NCC delamination and may provide new insights into mechanisms regulating pathological EMP.

Identification and characterization of intermediate states in mammalian neural crest cell epithelial to mesenchymal transition and delamination

Epithelial to mesenchymal transition (EMT) is a cellular process that converts epithelial cells to mesenchymal cells with migratory potential in both developmental and pathological processes. Although originally considered a binary event, EMT in cancer progression involves intermediate states between a fully epithelial and a fully mesenchymal phenotype, which are characterized by distinct combinations of epithelial and mesenchymal markers. This phenomenon has been termed epithelial to mesenchymal plasticity (EMP), however, the intermediate states remain poorly described and it's unclear whether they exist during developmental EMT. Neural crest cells (NCC) are an embryonic progenitor cell population that gives rise to numerous cell types and tissues in vertebrates, and their formation is a classic example of developmental EMT. An important feature of NCC development is their delamination from the neuroepithelium via EMT, following which NCC migrate throughout the embryo and undergo differentiation. NCC delamination shares similar changes in cellular state and structure with cancer cell invasion. However, whether intermediate states also exist during NCC EMT and delamination remains unknown. Through single cell RNA sequencing, we identified intermediate NCC states based on their transcriptional signature and then spatially defined their locations in situ in the dorsolateral neuroepithelium. Our results illustrate the progressive transcriptional and spatial transitions from premigratory to migratory cranial NCC during EMT and delamination. Of note gene expression and trajectory analysis indicate that distinct intermediate populations of NCC delaminate in either S phase or G2/M phase of the cell cycle, and the importance of cell cycle regulation in facilitating mammalian cranial NCC delamination was confirmed through cell cycle inhibition studies. Additionally, transcriptional knockdown revealed a functional role for the intermediate stage marker Dlc1 in regulating NCC delamination and migration. Overall, our work identifying and characterizing the intermediate cellular states, processes, and molecular signals that regulate mammalian NCC EMT and delamination furthers our understanding of developmental EMP and may provide new insights into mechanisms regulating pathological EMP.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Developmental origin of the mammalian premaxilla

The evolution of jaws has played a major role in the success of vertebrate expansion into a wide variety of ecological niches. A fundamental, yet unresolved, question in craniofacial biology is about the origin of the premaxilla, the most distal bone present in the upper jaw of all amniotes. Recent reports have suggested that the mammalian premaxilla is derived from embryonic maxillary prominences rather than the frontonasal ectomesenchyme as previously shown in studies of chicken embryos. However, whether mammalian embryonic frontonasal ectomesenchyme contributes to the premaxillary bone has not been investigated and a tool to trace the contributions of the frontonasal ectomesenchyme to facial structures in mammals is lacking. The expression of the Alx3 gene is activated highly specifically in the frontonasal ectomesenchyme, but not in the maxillary mesenchyme, from the beginning of facial morphogenesis in mice. Here, we report the generation and characterization of a novel Alx3CreERT2 knock-in mouse line that express tamoxifen-inducible Cre DNA recombinase from the Alx3 locus. Tamoxifen treatment of Alx3CreERT2/+;Rosa26mTmG/+ embryos at E7.5, E8.5, E9.5, and E10.5, each induced specific labeling of the embryonic medial nasal and lateral nasal mesenchyme but not the maxillary mesenchyme. Lineage tracing of Alx3CreERT2-labeled frontonasal mesenchyme from E9.5 to E16.5 clearly showed that the frontonasal mesenchyme cells give rise to the osteoblasts generating the premaxillary bone. Furthermore, we characterize a Dlx1-Cre BAC transgenic mouse line that expresses Cre activity in the embryonic maxillary but not the frontonasal mesenchyme and show that the Dlx1-Cre labeled embryonic maxillary mesenchyme cells contribute to the maxillary bone as well as the soft tissues lateral to both the premaxillary and maxillary bones but not to the premaxillary bone. These results clearly demonstrate the developmental origin of the premaxillary bone from embryonic frontonasal ectomesenchyme cells in mice and confirm the evolutionary homology of the premaxilla across amniotes.

Effect of TCDD exposure in adult female and male mice on the expression of miRNA in the ovaries and testes and associated reproductive functions

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant found widely across the world. While animal and human studies have shown that exposure to TCDD may cause significant alterations in the reproductive tract, the effect of TCDD on the expression of miRNA in the reproductive organs has not been previously tested. In the current study, we exposed adult female or male mice to TCDD or vehicle and bred them to study the impact on reproduction. The data showed that while TCDD treatment of females caused no significant change in litter size, it did alter the survival of the pups. Also, TCDD exposure of either the male or female mice led to an increase in the gestational period. While TCDD did not alter the gross morphology of the ovaries and testes, it induced significant alterations in the miRNA expression. The ovaries showed the differential expression of 426 miRNAs, of which 315 miRNAs were upregulated and 111 miRNA that were downregulated after TCDD exposure when compared to the vehicle controls. In the testes, TCDD caused the differential expression of 433 miRNAs, with 247 miRNAs upregulated and 186 miRNAs downregulated. Pathway analysis showed that several of these dysregulated miRNAs targeted reproductive functions. The current study suggests that the reproductive toxicity of TCDD may result from alterations in the miRNA expression in the reproductive organs. Because miRNAs also represent one of the epigenetic pathways of gene expression, our studies suggest that the transgenerational toxicity of TCDD may also result from dysregulation in the miRNAs.

Linking Vertebrate Gene Duplications to the New Head Hypothesis

Vertebrates have diverse morphologies and various anatomical novelties that set them apart from their closest invertebrate relatives. A conspicuous head housing a large brain, paired sense organs, and protected by a skeleton of cartilage and bone is unique to vertebrates and is a defining feature of this taxon. Gans and Northcutt (1980s) proposed that the evolution of this "new head" was dependent on two key developmental innovations: neural crest cells (NCCs) and ectodermal placodes. NCCs are migratory embryonic cells that form bone, cartilage, and neurons in the new head. Based on genome size, Ohno (1970s) proposed a separate hypothesis, stating that vertebrate genome content was quadrupled via two rounds (2R) of whole genome duplications (WGDs), and the surplus of genetic material potentiated vertebrate morphological diversification. While both hypotheses offer explanations for vertebrate success, it is unclear if, and how, the "new head" and "2R" hypotheses are linked. Here, we consider both hypotheses and evaluate the experimental evidence connecting the two. Overall, evidence suggests that while the origin of the NC GRN predates the vertebrate WGDs, these genomic events may have potentiated the evolution of distinct genetic subnetworks in different neural crest subpopulations. We describe the general composition of the NC GRN and posit that its increased developmental modularity facilitated the independent evolution of NC derivatives and the diversification of the vertebrate head skeleton. Lastly, we discuss experimental strategies needed to test whether gene duplications drove the diversification of neural crest derivatives and the "new head".

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Estimates of Genetic Parameters for Shape Space Data in Franches-Montagnes Horses

Simple Summary

Equine breeding is often based on conformation traits, describing the proportions, shape and joint angles of a horse. These conformations traits are, however, mostly subjectively judged and not measured objectively, affecting the response of selection through lower heritabilities and precision. In this study, we measured joint angles, quantified the variation in shape of 608 Swiss Franches-Montagnes (FM) horses and estimated the heritability of these traits. We found that the poll angle had the highest heritability of all joint angles (h² = 0.37), and variation in shape describing the type (heavy–light) was also fairly heritable (h² = 0.36–0.37). Furthermore, the shape of the FM stallions has clearly evolved towards a lighter type from 1940 to 2018 without stabilisation in recent years, risking the loss of the light draught horse type. Phenotyping based on photographs allowed us to improve the accuracy of certain joint angle traits, and to monitor the conformational development of the FM breed.

Abstract

Conformation traits such as joint angles are important selection criteria in equine breeding, but mainly consist of subjective evaluation scores given by breeding judges, showing limited variation. The horse shape space model extracts shape data from 246 landmarks (LM) and objective joint angle measurements from triplets of LM on standardized horse photographs. The heritability was estimated for 10 joint angles (seven were measured twice with different LM placements), and relative warp components of the whole shape, in 608 Franches-Montagnes (FM) horses (480 stallions, 68 mares and 60 geldings born 1940–2018, 3–25 years old). The pedigree data comprised 6986 horses. Genetic variances and covariances were estimated by restricted maximum likelihood model (REML), including the fixed effects birth year, age (linear and quadratic), height at withers (linear and quadratic), as well as postural effects (head, neck, limb position and body alignment), together with a random additive genetic animal component and the residual effect. Estimated heritability varied from 0.08 (stifle joint) to 0.37 (poll). For the shape, the type was most heritable (0.36 to 0.37) and evolved from heavy to light over time. Image-based phenotyping can improve the selection of horses for conformation traits with moderate heritability (e.g., poll, shoulder and fetlock).

ALX1-transcribed LncRNA AC132217.4 promotes osteogenesis and bone healing via IGF-AKT signaling in mesenchymal stem cells

The osteogenic potential of bone marrow mesenchymal stem cells (BMSCs) is critical for bone formation and regeneration. A high non-/delayed-union rate of fracture healing still occurs in specific populations, implying an urgent need to discover novel targets for promoting osteogenesis and bone regeneration. Long non-coding (lnc)RNAs are emerging regulators of multiple physiological processes, including osteogenesis. Based on differential expression analysis of RNA sequencing data, we found that lncRNA AC132217.4, a 3'UTR-overlapping lncRNA of insulin growth factor 2 (IGF2), was highly induced during osteogenic differentiation of BMSCs. Afterward, both gain-of-function and loss-of-function experiments proved that AC132217.4 promotes osteoblast development from BMSCs. As for its molecular mechanism, we found that AC132217.4 binds with IGF2 mRNA to regulate its expression and downstream AKT activation to control osteoblast maturation and function. Furthermore, we identified two splicing factors, splicing component 35 KDa (SC35) and heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), which regulate the biogenesis of AC132217.4 at the post-transcriptional level. We also identified a transcription factor, ALX1, which regulates AC132217.7 expression at the transcriptional level to promote osteogenesis. Importantly, in-vivo over-expression of AC132217.4 essentially promotes the bone healing process in a murine tibial drill-hole model. Our study demonstrates that lncRNA AC132217.4 is a novel anabolic regulator of BMSC osteogenesis and could be a plausible therapeutic target for improving bone regeneration.