Hox genes specify vertebral types in the presomitic mesoderm

Ventral body wall closure: Mechanistic insights from mouse models and translation to human pathology

The ventral body wall (VBW) that encloses the thoracic and abdominal cavities arises by extensive cell movements and morphogenetic changes during embryonic development. These morphogenetic processes include embryonic folding generating the primary body wall; the initial ventral cover of the embryo, followed by directed mesodermal cell migrations, contributing to the secondary body wall. Clinical anomalies in VBW development affect approximately 1 in 3000 live births. However, the cell interactions and critical cellular behaviors that control VBW development remain little understood. Here, we describe the embryonic origins of the VBW, the cellular and morphogenetic processes, and key genes, that are essential for VBW development. We also provide a clinical overview of VBW anomalies, together with environmental and genetic influences, and discuss the insight gained from over 70 mouse models that exhibit VBW defects, and their relevance, with respect to human pathology. In doing so we propose a phenotypic framework for researchers in the field which takes into account the clinical picture. We also highlight cases where there is a current paucity of mouse models for particular clinical defects and key gaps in knowledge about embryonic VBW development that need to be addressed to further understand mechanisms of human VBW pathologies.

Prevalence and characteristics of lumbar ribs: a meta-analysis with anatomical and clinical considerations

BACKGROUND

Lumbar ribs (LR) are a rare and relatively unknown anatomical abnormality of the lumbar spine. The literature provides better understanding regarding other spinal congenital variations like cervical ribs or lumbosacral transitional vertebrae, which are rather commonly recognised conditions. Thus, this meta-analysis aimed to provide data on prevalence and key characteristics of LR.

METHODS

Relevant databases were systematically searched for studies reporting the prevalence, laterality and geographic distribution of LR. No exclusion criteria based on language and date of original articles were employed. The pooled prevalence estimates (PPE) were calculated using a random-effects model. To assess the between-study heterogeneity, the I2 statistic and Chi-square test were utilized. Throughout the investigation, the PRISMA guidelines were adhered to scrupulously. Evaluation of the included studies' reliability was made with the AQUA tool.

RESULTS

In total, 9 studies were included in this meta-analysis. The pooled prevalence estimate (PPE) of LR was 2.1% (95%CI: 1.0-4.6). In studies based on CT imaging, LR were found in 1.6% (95%CI: 0.6-4.3) of patients and in Xray based studies in 2.1% (95%CI: 0.4-11.1) of patients. Lumbar ribs were bilateral in majority of individuals (65.4%, 95%CI: 39.4-84.6) and could be most frequently encountered in Europe with PPE of 2.8% (95%CI: 3.0-20.0), then in East Asia with PPE of 1.5% (95%CI: 1.0-19.2) and Middle East with PPE of 1.1% (95%CI: 0.6-20.0).

CONCLUSIONS

The findings of our study indicate that LR are a common anatomical variation of lumbar spine, contrary to previous beliefs. In a clinical practice, when a patient presents with a non-specific low back pain, a possible occurrence of LR should be taken into consideration. The presence of LR may be misleading for surgeons and result in wrong-level spine surgeries.

Homeotic and nonhomeotic patterns in the tetrapod vertebral formula

Vertebrate development and phylogeny are intimately connected through the vertebral formula, the numerical distribution of vertebrae along the body axis into different categories such as neck and chest. A key window into this relationship is through the conserved Hox gene clusters. Hox gene expression boundaries align with vertebral boundaries, and their manipulation in model organisms often results in the transformation of one vertebral type into its neighbor, a homeotic transformation. If the variety in the vertebrate body plan is produced by homeotic shifts, then the number of adjacent vertebrae will be inversely related when making interspecies comparisons since the gain in one vertebra is due to the loss in its neighbor. To date, such a pattern across species consistent with homeotic transitions has only been found in the thoracolumbar vertebral count of mammals. To further investigate potential homeotic relationships in other vertebrate classes and along the entire body axis, we compiled a comprehensive dataset of complete tetrapod vertebral formulas and systematically searched for patterns by analyzing combinations of vertebrae. We uncovered mammalian homeotic patterns and found balances between distal vertebrae not anticipated by a Hox-vertebral homeotic relationship, including one that emerged during the progression from theropods to birds. We also identified correlations between vertebral counts and intergenic distances in the HoxB gene cluster which do not align with the common picture of a colinear relationship between Hox expression and vertebral categories. This quantitative approach revises our expectations for the diversity of a Hox-mediated vertebrate body plan.

Conservation of rib skeleton regionalization in the homoplastic evolution of the snake-like body form in squamates

Squamates have independently evolved an elongate, limb-reduced body form numerous times. This transition has been proposed to involve either changes to regulatory gene expression or downstream modification of target enhancers to produce a homogeneous, deregionalized axial skeleton. Analysis of vertebral morphology has suggested that regionalization is maintained in snake-like body forms, but morphological variation in the other primary component of the axial skeleton, the dorsal ribs, has not been previously examined. We quantified rib morphology along the anterior-posterior axis in limbed and snake-like squamates to test different regionalization models. We find that the relative position of regional boundaries remains consistent across taxa of differing body types, including in the homoplastic evolution of snake-like body forms. The consistent retention of regional boundaries in this primaxial domain is uncorrelated with more plastic abaxial region markers. Rather than loss of regions, rib shape at the anterior and posterior of the axis converges on those in the middle, resulting in axial regions being distinguishable by allometric shape changes rather than by discrete morphologies. This complexity challenges notions of deregionalization, revealing a nuanced evolutionary history shaped by shared functions.

40 years of the homeobox: mechanisms of Hox spatial-temporal collinearity in vertebrates

Animal body plans are established during embryonic development by the Hox genes. This patterning process relies on the differential expression of Hox genes along the head-to-tail axis. Hox spatial collinearity refers to the relationship between the organization of Hox genes in clusters and the differential Hox expression, whereby the relative order of the Hox genes within a cluster mirrors the spatial sequence of expression in the developing embryo. In vertebrates, the cluster organization is also associated with the timing of Hox activation, which harmonizes Hox expression with the progressive emergence of axial tissues. Thereby, in vertebrates, Hox temporal collinearity is intimately linked to Hox spatial collinearity. Understanding the mechanisms contributing to Hox temporal and spatial collinearity is thus key to the comprehension of vertebrate patterning. Here, we provide an overview of the main discoveries pertaining to the mechanisms of Hox spatial-temporal collinearity.

Conserved patterns and locomotor-related evolutionary constraints in the hominoid vertebral column

The evolution of the hominoid lineage is characterized by pervasive homoplasy, notably in regions such as the vertebral column, which plays a central role in body support and locomotion. Few isolated and fewer associated vertebrae are known for most fossil hominoid taxa, but identified specimens indicate potentially high levels of convergence in terms of both form and number. Homoplasy thus complicates attempts to identify the anatomy of the last common ancestor of hominins and other taxa and stymies reconstructions of evolutionary scenarios. One way to clarify the role of homoplasy is by investigating constraints via phenotypic integration, which assesses covariation among traits, shapes evolutionary pathways, and itself evolves in response to selection. We assessed phenotypic integration and evolvability across the subaxial (cervical, thoracic, lumbar, sacral) vertebral column of macaques (n = 96), gibbons (n = 77), chimpanzees (n = 92), and modern humans (n = 151). We found a mid-cervical cluster that may have shifted cranially in hominoids, a persistent thoracic cluster that is most marked in chimpanzees, and an expanded lumbosacral cluster in hominoids that is most expanded in gibbons. Our results highlight the highly conserved nature of the vertebral column. Taxa appear to exploit existing patterns of integration and ontogenetic processes to shift, expand, or reduce cluster boundaries. Gibbons appear to be the most highly derived taxon in our sample, possibly in response to their highly specialized locomotion.

Hox genes and patterning the vertebrate body

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

Regionalization of the vertebral column and its correlation with heart position in snakes: Implications for evolutionary pathways and morphological diversification

Spinal regionalization has important implications for the evolution of vertebrate body plans. We determined the variation in the number and morphology of vertebrae across the vertebral column (i.e., vertebral formula) for 63 snake species representing 13 families using intracolumnar variation in vertebral shape. Vertebral counts were used to determine the position of the heart, pylorus, and left kidney for each species. Across all species we observed a conspicuous midthoracic transition in vertebral shape, indicating four developmental domains of the precloacal vertebral column (cervical, anterior thoracic, posterior thoracic, and lumbar). Using phylogenetic analyses, the boundary between the anterior and posterior thoracic vertebrae was correlated with heart position. No associations were found between shifts in morphology of the vertebral column and either the pylorus or left kidney. We observed that among taxa, the number of preapex and postapex vertebrae could change independently from one another and from changes in the total number of precloacal vertebrae. Ancestral state reconstruction of the preapex and postapex vertebrae illustrated several evolutionary pathways by which diversity in the vertebral column and heart position have been attained. In addition, no conspicuous pattern was observed among the heart, pylorus, or kidney indicating that their relative positions to each other evolve independently. We conclude that snakes exhibit four morphologically distinct regions of the vertebral column. We discuss the implications of the forebody and hindbody vertebral formula on the morphological diversification of snakes.

The control of transitions along the main body axis

Although vertebrates display a large variety of forms and sizes, the mechanisms controlling the layout of the basic body plan are substantially conserved throughout the clade. Following gastrulation, head, trunk, and tail are sequentially generated through the continuous addition of tissue at the caudal embryonic end. Development of each of these major embryonic regions is regulated by a distinct genetic network. The transitions from head-to-trunk and from trunk-to-tail development thus involve major changes in regulatory mechanisms, requiring proper coordination to guarantee smooth progression of embryonic development. In this review, we will discuss the key cellular and embryological events associated with those transitions giving particular attention to their regulation, aiming to provide a cohesive outlook of this important component of vertebrate development.

Evaluation of spinopelvic parameters in patients with lumbosacral transitional vertebra: a cross sectional and comparative study

INTRODUCTION

Lumbosacral transitional vertebra (LSTV) is the most common congenital anomaly of the lumbosacral junction and is a frequent cause of back pain in young patients with a prevalence of 4.6% to 36% in different regions.

OBJECTIVE

The objective of this study was to evaluate spinopelvic parameters in patients with lumbosacral transitional vertebra and to compare them with the same parameters of low back ache patients without lumbosacral transitional vertebra.

METHODS

This was a cross-sectional and comparative study conducted among low back ache patients presenting to our tertiary care center. Low back ache patients presenting to the outpatient department of AIIMS Jodhpur were screened for LSTV using radiographs. The spinopelvic parameters of those with LSTV were measured using Surgimap software and compared with the parameters of low back ache patients without LSTV. An Independent sample t-test was done and p-values were calculated.

RESULTS

The spinopelvic parameters, pelvic incidence, pelvic tilt and lumbar lordosis differed significantly in the patients with LSTV. Pelvic incidence was higher in the group with LSTV (58.5+9.3) when compared to the group without LSTV (50+8.8) with a p-value (<0.001). Pelvic tilt was higher in the group with LSTV (19.4+8.8) when compared to the group without LSTV (13.6+7.8) with a p-value (0.001). Lumbar lordosis was significantly higher in the group with LSTV (57.6+13.2) when compared to the group without LSTV (50.7+12.2) with a p-value (0.007). No significant differences were obtained in sacral slope and Pelvic-incidence and lumbar lordosis mismatch.

CONCLUSION

LSTV alters the spinopelvic parameters. Altered spinopelvic parameters predispose to spondylolisthesis, degenerative disc disease, and facet joint arthritis and are important in preoperative planning in spine and pelvic surgeries.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Spatiotemporal Characterization of Human Early Intervertebral Disc Formation at Single-Cell Resolution

The intervertebral disc (IVD) acts as a fibrocartilaginous joint to anchor adjacent vertebrae. Although several studies have demonstrated the cellular heterogeneity of adult mature IVDs, a single-cell transcriptomic atlas mapping early IVD formation is still lacking. Here, the authors generate a spatiotemporal and single cell-based transcriptomic atlas of human IVD formation at the embryonic stage and a comparative mouse transcript landscape. They identify two novel human notochord (NC)/nucleus pulposus (NP) clusters, SRY-box transcription factor 10 (SOX10)⁺ and cathepsin K (CTSK)⁺ , that are distributed in the early and late stages of IVD formation and they are validated by lineage tracing experiments in mice. Matrisome NC/NP clusters, T-box transcription factor T (TBXT)⁺ and CTSK⁺ , are responsible for the extracellular matrix homeostasis. The IVD atlas suggests that a subcluster of the vertebral chondrocyte subcluster might give rise to an inner annulus fibrosus of chondrogenic origin, while the fibroblastic outer annulus fibrosus preferentially expresseds transgelin and fibromodulin . Through analyzing intercellular crosstalk, the authors further find that notochordal secreted phosphoprotein 1 (SPP1) is a novel cue in the IVD microenvironment, and it is associated with IVD development and degeneration. In conclusion, the single-cell transcriptomic atlas will be leveraged to develop preventative and regenerative strategies for IVD degeneration.

Evolution of vertebral numbers in primates, with a focus on hominoids and the last common ancestor of hominins and panins

The primate vertebral column has been extensively studied, with a particular focus on hominoid primates and the last common ancestor of humans and chimpanzees. The number of vertebrae in hominoids-up to and including the last common ancestor of humans and chimpanzees-is subject to considerable debate. However, few formal ancestral state reconstructions exist, and none include a broad sample of primates or account for the correlated evolution of the vertebral column. Here, we conduct an ancestral state reconstruction using a model of evolution that accounts for both homeotic (changes of one type of vertebra to another) and meristic (addition or loss of a vertebra) changes. Our results suggest that ancestral primates were characterized by 29 precaudal vertebrae, with the most common formula being seven cervical, 13 thoracic, six lumbar, and three sacral vertebrae. Extant hominoids evolved tail loss and a reduced lumbar column via sacralization (homeotic transition at the last lumbar vertebra). Our results also indicate that the ancestral hylobatid had seven cervical, 13 thoracic, five lumbar, and four sacral vertebrae, and the ancestral hominid had seven cervical, 13 thoracic, four lumbar, and five sacral vertebrae. The last common ancestor of humans and chimpanzees likely either retained this ancestral hominid formula or was characterized by an additional sacral vertebra, possibly acquired through a homeotic shift at the sacrococcygeal border. Our results support the 'short-back' model of hominin vertebral evolution, which postulates that hominins evolved from an ancestor with an African ape-like numerical composition of the vertebral column.

Prevalence of lumbosacral transitional vertebra among 4816 consecutive patients with low back pain: A computed tomography, magnetic resonance imaging, and plain radiographic study with novel classification schema

Study Design

A retrospective single-center study.

Background

The prevalence of the lumbosacral anomalies remains controversial. The existing classification to characterize these anomalies is more complex than necessary for clinical use.

Purpose

To assessment of the prevalence of lumbosacral transitional vertebra (LSTV) in patients with low back pain and the development of clinically relevant classification to describe these anomalies.

Materials and Methods

During the period from 2007 to 2017, all cases of LSTV were preoperatively verified, and classified according to Castellvi, as well as O'Driscoll. We then developed modifications of those classifications that are simpler, easier to remember, and clinically relevant. At the surgical level, this was assessed intervertebral disc and facet joint degeneration.

Results

The prevalence of the LSTV was 8.1% (389/4816). The most common L5 transverse process anomaly type was fused, unilaterally or bilaterally (48%), to the sacrum and were O'Driscoll's III (40.1%) and IV (35.8%). The most common type of S1-2 disc was a lumbarized disc (75.9%), where the disc's anterior-posterior diameter was equal to the L5-S1 disc diameter. In most cases, neurological compression symptoms (85.5%) were verified to be due to spinal stenosis (41.5%) or herniated disc (39.5%). In the majority of patients without neural compression, the clinical symptoms were due to mechanical back pain (58.8%).

Conclusions

LSTV is a fairly common pathology of the lumbosacral junction, occurring in 8.1% of the patients in our series (389 out of 4,816 cases). The most common types were Castellvi's type IIA (30.9%) and IIIA (34.9%) and were O'Driscoll's III (40.1%) and IV (35.8%).

Nr6a1 controls Hox expression dynamics and is a master regulator of vertebrate trunk development

The vertebrate main-body axis is laid down during embryonic stages in an anterior-to-posterior (head-to-tail) direction, driven and supplied by posteriorly located progenitors. Whilst posterior expansion and segmentation appears broadly uniform along the axis, there is developmental and evolutionary support for at least two discrete modules controlling processes within different axial regions: a trunk and a tail module. Here, we identify Nuclear receptor subfamily 6 group A member 1 (Nr6a1) as a master regulator of trunk development in the mouse. Specifically, Nr6a1 was found to control vertebral number and segmentation of the trunk region, autonomously from other axial regions. Moreover, Nr6a1 was essential for the timely progression of Hox signatures, and neural versus mesodermal cell fate choice, within axial progenitors. Collectively, Nr6a1 has an axially-restricted role in all major cellular and tissue-level events required for vertebral column formation, supporting the view that changes in Nr6a1 levels may underlie evolutionary changes in axial formulae.

Shaping Hox gene activity to generate morphological diversity across vertebrate phylogeny

The importance of Hox genes for the development and evolution of the vertebrate axial skeleton and paired appendages has been recognized for already several decades. The steady growth of genomic sequence data from an increasing number of vertebrate species, together with the improvement of methods to analyze genomic structure and interactions, as well as to control gene activity in various species has refined our understanding of Hox gene activity in development and evolution. Here, I will review recent data addressing the influence of Hox regulatory processes in the evolution of the fins and the emergence of the tetrapod limb. In addition, I will discuss the involvement of posterior Hox genes in the control of vertebrate axial extension, focusing on an apparently divergent activity that Hox13 paralog group genes have on the regulation of tail bud development in mouse and zebrafish embryos.

The pioneering function of the hox transcription factors

Ever since the discovery that the Hox family of transcription factors establish morphological diversity in the developing embryo, major efforts have been directed towards understanding Hox-dependent patterning. This has led to important discoveries, notably on the mechanisms underlying the collinear expression of Hox genes and Hox binding specificity. More recently, several studies have provided evidence that Hox factors have the capacity to bind their targets in an inaccessible chromatin context and trigger the switch to an accessible, transcriptional permissive, chromatin state. In this review, we provide an overview of the evidences supporting that Hox factors behave as pioneer factors and discuss the potential mechanisms implicated in Hox pioneer activity as well as the significance of this functional property in Hox-dependent patterning.

Seeking Sense in the Hox Gene Cluster

The Hox gene cluster, responsible for patterning of the head-tail axis, is an ancestral feature of all bilaterally symmetrical animals (the Bilateria) that remains intact in a wide range of species. We can say that the Hox cluster evolved successfully only once since it is commonly the same in all groups, with labial-like genes at one end of the cluster expressed in the anterior embryo, and Abd-B-like genes at the other end of the cluster expressed posteriorly. This review attempts to make sense of the Hox gene cluster and to address the following questions. How did the Hox cluster form in the protostome-deuterostome last common ancestor, and why was this with a particular head-tail polarity? Why is gene clustering usually maintained? Why is there collinearity between the order of genes along the cluster and the positions of their expressions along the embryo? Why do the Hox gene expression domains overlap along the embryo? Why have vertebrates duplicated the Hox cluster? Why do Hox gene knockouts typically result in anterior homeotic transformations? How do animals adapt their Hox clusters to evolve new structural patterns along the head-tail axis?

What snakes and caecilians have in common? Molecular interaction units and the independent origins of similar morphotypes in Tetrapoda

Developmental pathways encompass transcription factors and cis-regulatory elements that interact as transcription factor-regulatory element (TF-RE) units. Independent origins of similar phenotypes likely involve changes in different parts of these units, a hypothesis promisingly tested addressing the evolution of the rib-associated lumbar (RAL) morphotype that characterizes emblematic animals such as snakes and elephants. Previous investigation in these lineages identified a polymorphism in the Homology region 1 [H1] enhancer of the Myogenic factor-5 [Myf5], which interacts with HOX10 proteins to modulate rib development. Here we address the evolution of TF-RE units focusing on independent origins of RAL morphotypes. We compiled an extensive database for H1-Myf5 and HOX10 sequences with two goals: (i) evaluate if the enhancer polymorphism is present in amphibians exhibiting the RAL morphotype and (ii) test a hypothesis of enhanced evolutionary flexibility mediated by TF-RE units, according to which independent origins of the RAL morphotype might involve changes in either component of the interaction unit. We identified the H1-Myf5 polymorphism in lineages that diverged around 340 Ma, including Lissamphibia. Independent origins of the RAL morphotype in Tetrapoda involved sequence variation in either component of the TF-RE unit, confirming that different changes may similarly affect the phenotypic outcome of a given developmental pathway.

Sacralization of Coccygeal Vertebra: A Descriptive Observational Study in Bangladesh

Background: In the sacrococcygeal region, anatomical variation is due to the sacralization of the coccygeal vertebra, which is the due union of/fusion of the fifth sacral with the first coccygeal vertebra of five couples of sacral foramina under-detected or asymptomatic beyond radiological assessment. That is why it is challenging to know the cause of coccydynia, caudal block failure, the difficult second stage of labor, and perineal tears. The present study aims to improve knowledge about the anatomical variation of sacralization of the coccygeal vertebra. Additionally, to find the prevalence of sacralization of coccygeal vertebra in Sylhet, Bangladesh.

Methods: This study was performed on 60 parched, totally calcified, typical sacra of mature-age individuals of undetermined sexes, fulfilling the inclusion criteria from the bone bank of the osteology museum of the Department of Anatomy, Sylhet MAG Osmani Medical College, Sylhet, Bangladesh, from July 2017 to June 2018. Sex determination of the collected unknown sacra was conducted using discriminant function analysis. It was found that 50% (30) were male and 50% (30%) were female. The unpaired t-tests and chi-square were utilized to determine the statistical significance.

Results: Out of 60 sacra, eight (13.33%) samples presented with sacralization. This study found that males had significantly higher straight (p=0.05) and curved (p=0.05) lengths of sacrococcygeal vertebrae. The sacrococcygeal curvature index (SCI) showed statistically significant (p=0.05) differences between the sexes.

Conclusion: Sacralization may exert an impact on the caudal block. It could extend the second stage of the labor process with perineal tears. Therefore, knowledge about the anatomical variation of the coccygeal vertebra is essential.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Lumbosacral Transitional Vertebra Contributed to Lumbar Spine Degeneration: An MR Study of Clinical Patients

We aimed to comprehensively characterize degenerative findings associated with various types of lumbosacral transitional vertebra (LSTV) on magnetic resonance images. Three hundred and fifty patients with LSTV (52.3 ± 10.9 years), including 182 Castellvi type I, 107 type II, 43 type III, and 18 type IV, and 179 controls without LSTV (50.6 ± 13.1 years), were studied. Discs, endplates, and posterior vertebral structures were assessed and compared to those of controls for the most caudal three discs on MRIs. There were no differences in degenerative findings between patients with type I LSTV and controls. For types III and IV, the transitional discs had smaller sizes, lower Pfirrmann scores, and lower rates of disc bulging (2.3% and 5.6% vs. 39.1%), osteophytes (2.3% vs. 15.1%), disc herniation (2.3% and 5.6% vs. 31.8%), and Modic changes (2.3% and 5.6% vs. 16.8%) than controls. However, the cranial discs had more severe Pfirrmann scores, disc narrowing and spinal canal narrowing, and greater rates of disc herniation (41.9% and 50.0% vs. 25.7%), endplate defects (27.9% and 33.3% vs. 14.4%) and spondylolisthesis (18.6% vs. 7.3%) than controls. Type II LSTV was associated with degenerative findings in the cranial segments but to a lesser degree, as compared with type III/IV LSTV. Thus, Castellvi type III/IV LSTV predisposed the adjacent spinal components to degeneration and protected the transitional discs. Type II LSTV had significant effects in promoting transitional and adjacent disc degeneration. Type I LSTV was not related to spinal degeneration.

Increased vertebral body area, disc and facet joint degeneration throughout the lumbar spine in patients with lumbosacral transitional vertebrae

OBJECTIVES

Lumbosacral transitional vertebrae (LSTV), developmental spinal canal stenosis, lumbar disc and facet joint degeneration are all common lumbar spine conditions. This study explores the relationship between these entities.

METHODS

Five hundred and five patients (254 females, 251 males), mean age 57 ± 13 years, who underwent abdominopelvic CT examinations were studied. Vertebral body cross-sectional area (CSA), spinal canal CSA, width and depth were measured from L1 to L5. The presence of LSTV was documented and the severity of disc/facet joint degeneration graded.

RESULTS

LSTV was present in 155 (31%) of 505 patients, being twice as common in males than in females. After controlling for age, height and weight, vertebral body CSAs from L1 to L5 in male LSTV patients were, on average, 9.8% larger than those in male non-LSTV patients (p < 0.0001 for all) while vertebral body areas at L3 and L4 were 4.5% larger in female LSTV patients than in non-LSTV female patients. Disc degeneration from L2/3 to L5/S1 and facet joint degeneration from L2/3 to L4/5 were more severe in LSTV patients. LSTV Castellvi Type II patients appear to be most prone to disc and facet degeneration.

CONCLUSIONS

LSTV patients, and particularly male patients, have an increased vertebral body CSA compared to non-LSTV patients throughout the lumbar spine. LSTV patients also have more severe disc and facet joint degeneration. The increase in vertebral body area may be contributing to the increased lumbar and facet joint degeneration seen in LSTV patients.

KEY POINTS

• LSTV patients have increased vertebral body cross-sectional area throughout their lumbar spine compared to non-LSTV patients. This vertebral body area increase was more pronounced in male patients and also apparent in younger patients with LSTV. • LSTV patients have increased disc and facet joint degeneration affecting most of the lumbar spine than non-LSTV patients. • Increased vertebral body area may contribute to the increased lumbar disc and facet joint degeneration seen in LSTV patients.

The morphological consequences of segmentation anomalies in the human sacrum

OBJECTIVES

Despite the high frequency of segmentation anomalies in the human sacrum, their evolutionary and clinical implications remain controversial. Specifically, inconsistencies involving the classification and counting methods obscure accurate assessment of lumbosacral transitional vertebrae. Therefore, we aim to establish more reliable morphological and morphometric methods for differentiating between sacralizations and lumbarizations in clinical and paleontological contexts.

MATERIALS AND METHODS

Using clinical CT data from 145 individuals aged 14-47 years, vertebral counts and the spatial relationship between the sacrum and adjoining bony structures were assessed, while the morphological variation of the sacrum was assessed using geometric morphometrics based on varied landmark configurations.

RESULTS

The prevalence of lumbosacral and sacrococcygeal segmentation anomalies was 40%. Lumbarizations and sacralizations were reliably distinguishable based on the spatial relationship between the iliac crest and the upward or downward trajectory of the linea terminalis on the sacrum. Different craniocaudal orientations of the alae relative to the corpus of the first sacral vertebra were also reflected in the geometric morphometric analyses. The fusion of the coccyx (32%) was frequently coupled with lumbarizations, suggesting that the six-element sacra more often incorporate the coccyx rather than the fifth lumbar vertebra.

CONCLUSIONS

Our approach allowed the consistent identification of segmentation anomalies even in isolated sacra. Additionally, our outcomes either suggest that homeotic border shifts often affect multiple spinal regions in a unidirectional way, or that sacrum length is highly conserved perhaps due to functional constraints. Our results elucidate the potential clinical, biomechanical, and evolutionary significance of lumbosacral transitional vertebrae.

Regulation of posterior Hox genes by sex steroids explains vertebral variation in inbred mouse strains

A series of elegant embryo transfer experiments in the 1950s demonstrated that the uterine environment could alter vertebral patterning in inbred mouse strains. In the intervening decades, attention has tended to focus on the technical achievements involved and neglected the underlying biological question: how can genetically homogenous individuals have a heterogenous number of vertebrae? Here I revisit these experiments and, with the benefit of knowledge of the molecular-level processes of vertebral patterning gained over the intervening decades, suggest a novel hypothesis for homeotic transformation of the last lumbar vertebra to the adjacent sacral type through regulation of Hox genes by sex steroids. Hox genes are involved in both axial patterning and development of male and female reproductive systems and have been shown to be sensitive to sex steroids in vitro and in vivo. Regulation of these genes by sex steroids and resulting alterations to vertebral patterning may hint at a deep evolutionary link between the ribless lumbar region of mammals and the switch from egg-laying to embryo implantation. An appreciation of the impact of sex steroids on Hox genes may explain some puzzling aspects of human disease, and highlights the spine as a neglected target for in utero exposure to endocrine disruptors.

Quantitative measurements at the lumbosacral junction are more reliable parameters for identifying and numbering lumbosacral transitional vertebrae

OBJECTIVES

To evaluate quantitative parameters to identify the anatomic variation lumbosacral transitional vertebrae (LSTV) and compare them with the landmarks commonly used at present.

METHODS

A total of 2,845 PET/CT scans were reviewed, and the patients with 23 and 25 presacral vertebrae were included. The quantitative parameters, including the anterior-edge vertebral angle (AVA) of the lowest lumbar-type vertebra, the ratio of the length of the inferior endplate to that of the superior endplate (RISE) of the uppermost sacral-type vertebra and the lumbosacral intervertebral disc angle (LSIVDA), and the anatomical landmarks, including the iliac crest tangent (ICT) level, the iliolumbar ligament (ILL) origin level and psoas proximal insertion, were all evaluated to determine their ability to identify LSTV.

RESULTS

The values of AVA and RISE were significantly different between the LSTV group and the control group, and between subgroups of LSTV. The cutoff value for AVA was 73.0°, with an accuracy, sensitivity, and specificity of 91.1%, 77.5%, and 88.3%, and that for RISE was 0.79, with an accuracy, sensitivity, and specificity of 90.3%, 77.5%, and 94.2%, while that for LSIVDA was 14.15°, with an accuracy, sensitivity, and specificity of 75.9%, 65.7%, and 78.3%, to differentiate L5 sacralization from S1 lumbarization. For differentiating the controls from LSTV, the accuracy, sensitivity, and specificity of the ICT level and proximal psoas insertion were 78.0%, 70.2%, and 95.0%, versus 71.7%, 61.7%, and 94.0%.

CONCLUSIONS

Compared with the anatomical landmarks, the quantitative measurements at the lumbosacral junction, including AVA and RISE, may be more helpful for differentiating subgroups of LSTV especially if only lumbar spine imaging is available.

KEY POINTS

• The quantitative parameters, the anterior-edge vertebral angle (AVA) of the lowest lumbar-type vertebra and the ratio of the length of the inferior endplate to that of the superior endplate (RISE) of the uppermost sacral-type vertebra, are more helpful for distinguishing L5 sacralization from S1 lumbarization than the previously proposed anatomic landmarks. • AVA and RISE represent relevant changes in the curvature at the lumbosacral region and the shape of the transitional vertebral body, respectively. • AVA and RISE are easily assessed, with high intra- and inter-reader reliability.

Breaking constraint of mammalian axial formulae

The vertebral column of individual mammalian species often exhibits remarkable robustness in the number and identity of vertebral elements that form (known as axial formulae). The genetic mechanism(s) underlying this constraint however remain ill-defined. Here, we reveal the interplay of three regulatory pathways (Gdf11, miR-196 and Retinoic acid) is essential in constraining total vertebral number and regional axial identity in the mouse, from cervical through to tail vertebrae. All three pathways have differing control over Hox cluster expression, with heterochronic and quantitative changes found to parallel changes in axial identity. However, our work reveals an additional role for Hox genes in supporting axial elongation within the tail region, providing important support for an emerging view that mammalian Hox function is not limited to imparting positional identity as the mammalian body plan is laid down. More broadly, this work provides a molecular framework to interrogate mechanisms of evolutionary change and congenital anomalies of the vertebral column.

Vertebral column length and shape exhibits remarkable robustness within a species but diversity across species. Here the authors reveal the molecular logic constraining vertebral number in mouse and a novel role for posterior Hox genes in this context.

[Analysis of unfavorable postoperative outcomes in patients with lumbosacral junction anomalies]

BACKGROUND

The absence of a unified diagnostic approach and heterogeneous treatment strategies in patients with lumbosacral transitional vertebra contribute to registration of contradictory data on postoperative outcomes in these patients.

OBJECTIVE

A retrospective analysis of postoperative outcomes in patients with degenerative lumbar spine diseases associated with lumbosacral junction anomalies and causes of unfavorable results.

MATERIAL AND METHODS

Lumbosacral transitional vertebra was preoperatively verified and classified according to Castellvi A.E. in 352 cases between 2007 and 2017. Full-time examination, phone and e-mail survey covered 314 (89.2%) patients throughout 5.2-year follow-up period. There were unfavorable clinical outcomes in 42 cases. We analyzed the causes of these results and relationship of surgical strategy with the type of lumbosacral junction anomaly.

RESULTS

Lumbosacral transitional vertebra Castellvi type IIA (30.9%) and IIIA (34.9%) prevailed. Unfavorable outcomes occurred in 7 cases after percutaneous interventions, in 12 patients after decompression surgery and in 23 patients after decompression and stabilization procedures. In all cases (n=42), redo surgery was performed: rigid fixation or extension of stabilization system in 16 cases, additional intracanal decompression in 10 patients, foraminal/extraforaminal decompression in 8 patients. Seven patients underwent repeated laser denervation. Of these, 5 patients underwent dorsal decompression without discectomy due to ineffective denervation. Correlation analysis showed the relationship between unfavorable outcomes after percutaneous procedures with L_V transverse process dysplasia and sacralization, decompressive interventions with additional articular joints and lumbarization, decompression/stabilization procedures with L_V transverse process dysplasia and sacralization.

CONCLUSION

Unfavorable outcomes occurred after percutaneous interventions for anomaly type IB and IIIB, decompressive procedures for anomaly type IIA, IIB and IV, decompression/ stabilization procedures for types IA and IIIA. It is the first report devoted to analysis of relationships between unsatisfactory postoperative outcomes, surgical technique and type of lumbosacral dysgenesis.

Role of the Lumbosacral Transition Vertebra and Vertebral Lamina in the Pathogenesis of Lumbar Disc Herniation

OBJECTIVE

To investigate the prevalence of lumbosacral transition vertebrae (LSTVs) in both the normal population and the lumbar disc herniation (LDH) population and to determine the risk factors for LDH.

METHODS

Between January 2019 and September 2020, all patients aged 18-39 years and underwent an anteroposterior (AP) X-ray of the lumbar vertebrae were retrospective reviewed in our institution. Those patients who were diagnosed with LDH were eligible for inclusion in the LDH group. During the same period, those patients admitted to our hospital who underwent an anteroposterior X-ray of the lumbar spine and had not been diagnosed with LDH were included in the control group. Those patients with disease that might affect the lumbar anatomy were excluded from both groups. The type of LSTV was classified according to the Castellvi classification. The height of the lumbar vertebral lamina was evaluated through the h/H index. The inter- and intra-observer reliability was evaluated by one senior radiologist and one senior orthopedist using intraclass correlation coefficient (ICC). The association between the LSTV and the herniation level was also investigated. Binary logistic regression was used to explore the association of different factors between the LDH group and the control group.

RESULTS

Two hundred LDH patients (115 male and 85 female) and 200 individuals (108 male and 92 female) were investigated retrospectively. The prevalence of LSTVs was 71.5% (n = 143) in the LDH group and 34.0% (n = 68) in the control group. The most frequent LSTV types were type Ib and type IIa. The inter- and intra-observer ICCs of the measurement of "h/H" index and the classification of LSTV were all "excellent" (ICC > 0.90). The median h/H index in the control group was significantly higher than that in the LDH group (0.28 (0.26, 0.31) vs 0.34 (0.31, 0.37), P = 0.000). The distribution of the Castellvi classification in the L4/5 and L5/S1 herniation patients was significantly different (P = 0.048). LSTVs, BMI and the h/H index were closely associated with LDH, with odds ratios of 3.06 (95% CI: 2.12-4.43), 1.23 (95% CI: 1.13-1.33) and 0.09 (95% CI: 0.05-0.15), respectively. The incidence of L4/5 disc herniation in patients with an LSTV was significantly more common than that in patients with L5/S1 disc herniation (P = 0.048).

CONCLUSION

The prevalence of LSTVs was 34.0% in the control group and 71.5% in the LDH group; LSTVs and BMI were positively correlated with LDH, and h/H was negatively correlated with LDH.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

Somite development and regionalisation of the vertebral axial skeleton

A critical stage in the development of all vertebrate embryos is the generation of the body plan and its subsequent patterning and regionalisation along the main anterior-posterior axis. This includes the formation of the vertebral axial skeleton. Its organisation begins during early embryonic development with the periodic formation of paired blocks of mesoderm tissue called somites. Here, we review axial patterning of somites, with a focus on studies using amniote model systems - avian and mouse. We summarise the molecular and cellular mechanisms that generate paraxial mesoderm and review how the different anatomical regions of the vertebral column acquire their specific identity and thus shape the body plan. We also discuss the generation of organoids and embryo-like structures from embryonic stem cells, which provide insights regarding axis formation and promise to be useful for disease modelling.

Interaction between retinoic acid and FGF/ERK signals are involved in Dexamethasone-induced abnormal myogenesis during embryonic development

Despite the common application in pregnancy at clinical practice, it remains ambiguous whether dexamethasone (Dex) exposure can affect embryonic myogenesis. In this study, firstly we showed that 10^-6 M Dex (Cheng et al., 2016; 2017) treatment resulted in abnormal myogenesis in chicken embryos. Secondly, we demonstrated that 10^-6 M Dex-induced abnormality of myogenesis resulted from aberrant cell proliferation, as well as from alteration of the differentiation process from the early stage of somitogenesis up to the late stage of myogenesis. The above-mentioned results caused by Dex exposure might be due to the aberrant gene expressions of somite formation (Raldh2, Fgf8, Wnt3a, β-catenin, Slug, Paraxis, N-cadherin) and differentiation (Pax3, MyoD, Wnt3a, Msx1, Shh). Thirdly, RNA sequencing implied the statistically significant differential gene expressions in regulating the myofibril and systemic development, as well as a dramatical alteration of retinoic acid (RA) signaling during somite development in the chicken embryos exposed to Dex. The subsequent validation experiments verified that Dex treatment indeed led to a metabolic change of RA signaling, which was up-regulated and principally mediated by FGF-ERK signaling revealed by means of the combination of chicken embryos and in vitro C2C12 cells. These findings highlight that 10^-6 M Dex exposure enhances the risk of abnormal myogenesis through interfering with RA signaling during development.

Homeotic change in segment identity derives the human vertebral formula from a chimpanzee-like one

OBJECTIVES

One of the most contentious issues in paleoanthropology is the nature of the last common ancestor of humans and our closest living relatives, chimpanzees and bonobos (panins). The numerical composition of the vertebral column has featured prominently, with multiple models predicting distinct patterns of evolution and contexts from which bipedalism evolved. Here, we study total numbers of vertebrae from a large sample of hominoids to quantify variation in and patterns of regional and total numbers of vertebrae in hominoids.

MATERIALS AND METHODS

We compile and study a large sample (N = 893) of hominoid vertebral formulae (numbers of cervical, thoracic, lumbar, sacral, caudal segments in each specimen) and analyze full vertebral formulae, total numbers of vertebrae, and super-regional numbers of vertebrae: presacral (cervical, thoracic, lumbar) vertebrae and sacrococcygeal vertebrae. We quantify within- and between-taxon variation using heterogeneity and similarity measures derived from population genetics.

RESULTS

We find that humans are most similar to African apes in total and super-regional numbers of vertebrae. Additionally, our analyses demonstrate that selection for bipedalism reduced variation in numbers of vertebrae relative to other hominoids.

DISCUSSION

The only proposed ancestral vertebral configuration for the last common ancestor of hominins and panins that is consistent with our results is the modal formula demonstrated by chimpanzees and bonobos (7 cervical-13 thoracic-4 lumbar-6 sacral-3 coccygeal). Hox gene expression boundaries suggest that a rostral shift in Hox10/Hox11-mediated complexes could produce the human modal formula from the proposal ancestral and panin modal formula.

Of Necks, Trunks and Tails: Axial Skeletal Diversity among Vertebrates

The axial skeleton of all vertebrates is composed of individual units known as vertebrae. Each vertebra has individual anatomical attributes, yet they can be classified in five different groups, namely cervical, thoracic, lumbar, sacral and caudal, according to shared characteristics and their association with specific body areas. Variations in vertebral number, size, morphological features and their distribution amongst the different regions of the vertebral column are a major source of the anatomical diversity observed among vertebrates. In this review I will discuss the impact of those variations on the anatomy of different vertebrate species and provide insights into the genetic origin of some remarkable morphological traits that often serve to classify phylogenetic branches or individual species, like the long trunks of snakes or the long necks of giraffes.

Evolution: Spinal Innovation Enabled by Genome Duplication

Solid vertebrae evolved multiple times across vertebrates, but the origins and relationships of different spine forms remain unclear. A new study reveals teleost fishes evolved their solid vertebrae following genome duplication, when a novel gene repressed ancestral spine programming.

Discovering a sparse set of pairwise discriminating features in high-dimensional data

MOTIVATION

Recent technological advances produce a wealth of high-dimensional descriptions of biological processes, yet extracting meaningful insight and mechanistic understanding from these data remains challenging. For example, in developmental biology, the dynamics of differentiation can now be mapped quantitatively using single-cell RNA sequencing, yet it is difficult to infer molecular regulators of developmental transitions. Here, we show that discovering informative features in the data is crucial for statistical analysis as well as making experimental predictions.

RESULTS

We identify features based on their ability to discriminate between clusters of the data points. We define a class of problems in which linear separability of clusters is hidden in a low-dimensional space. We propose an unsupervised method to identify the subset of features that define a low-dimensional subspace in which clustering can be conducted. This is achieved by averaging over discriminators trained on an ensemble of proposed cluster configurations. We then apply our method to single-cell RNA-seq data from mouse gastrulation, and identify 27 key transcription factors (out of 409 total), 18 of which are known to define cell states through their expression levels. In this inferred subspace, we find clear signatures of known cell types that eluded classification prior to discovery of the correct low-dimensional subspace.

AVAILABILITY AND IMPLEMENTATION

https://github.com/smelton/SMD.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Fossoriality and evolutionary development in two Cretaceous mammaliamorphs

Mammaliamorpha comprises the last common ancestor of Tritylodontidae and Mammalia plus all its descendants¹. Tritylodontids are nonmammaliaform herbivorous cynodonts that originated in the Late Triassic epoch, diversified in the Jurassic period^2-5 and survived into the Early Cretaceous epoch^6,7. Eutriconodontans have generally been considered to be an extinct mammalian group, although different views exist⁸. Here we report a newly discovered tritylodontid and eutriconodontan from the Early Cretaceous Jehol Biota of China. Eutriconodontans are common in this biota⁹, but it was not previously known to contain tritylodontids. The two distantly related species show convergent features that are adapted for fossorial life, and are the first 'scratch-diggers' known from this biota. Both species also show an increased number of presacral vertebrae, relative to the ancestral state in synapsids or mammals^10,11, that display meristic and homeotic changes. These fossils shed light on the evolutionary development of the axial skeleton in mammaliamorphs, which has been the focus of numerous studies in vertebrate evolution^12-17 and developmental biology^18-28. The phenotypes recorded by these fossils indicate that developmental plasticity in somitogenesis and HOX gene expression in the axial skeleton-similar to that observed in extant mammals-was already in place in stem mammaliamorphs. The interaction of these developmental mechanisms with natural selection may have underpinned the diverse phenotypes of body plan that evolved independently in various clades of mammaliamorph.

Retinoic Acid Signaling Plays a Crucial Role in Excessive Caffeine Intake-Disturbed Apoptosis and Differentiation of Myogenic Progenitors

Whether or not the process of somitogenesis and myogenesis is affected by excessive caffeine intake still remains ambiguous. In this study, we first showed that caffeine treatment results in chest wall deformities and simultaneously reduced mRNA expressions of genes involved in myogenesis in the developing chicken embryos. We then used embryo cultures to assess in further detail how caffeine exposure affects the earliest steps of myogenesis, and we demonstrated that the caffeine treatment suppressed somitogenesis of chicken embryos by interfering with the expressions of crucial genes modulating apoptosis, proliferation, and differentiation of myogenic progenitors in differentiating somites. These phenotypes were abrogated by a retinoic acid (RA) antagonist in embryo cultures, even at low caffeine doses in C2C12 cells, implying that excess RA levels are responsible for these phenotypes in cells and possibly in vivo. These findings highlight that excessive caffeine exposure is negatively involved in regulating the development of myogenic progenitors through interfering with RA signaling. The RA somitogenesis/myogenesis pathway might be directly impacted by caffeine signaling rather than reflecting an indirect effect of the toxicity of excess caffeine dosage.

Making and shaping endochondral and intramembranous bones

Skeletal elements have a diverse range of shapes and sizes specialized to their various roles including protecting internal organs, locomotion, feeding, hearing, and vocalization. The precise positioning, size, and shape of skeletal elements is therefore critical for their function. During embryonic development, bone forms by endochondral or intramembranous ossification and can arise from the paraxial and lateral plate mesoderm or neural crest. This review describes inductive mechanisms to position and pattern bones within the developing embryo, compares and contrasts the intrinsic vs extrinsic mechanisms of endochondral and intramembranous skeletal development, and details known cellular processes that precisely determine skeletal shape and size. Key cellular mechanisms are employed at distinct stages of ossification, many of which occur in response to mechanical cues (eg, joint formation) or preempting future load‐bearing requirements. Rapid shape changes occur during cellular condensation and template establishment. Specialized cellular behaviors, such as chondrocyte hypertrophy in endochondral bone and secondary cartilage on intramembranous bones, also dramatically change template shape. Once ossification is complete, bone shape undergoes functional adaptation through (re)modeling. We also highlight how alterations in these cellular processes contribute to evolutionary change and how differences in the embryonic origin of bones can influence postnatal bone repair.

Compares and contrasts Endochondral and intramembranous bone development

Reviews embryonic origins of different bones

Describes the cellular and molecular mechanisms of positioning skeletal elements.

Describes mechanisms of skeletal growth with a focus on the generation of skeletal shape

Development and patterning of rib primordia are dependent on associated musculature

The importance of skeletal muscle for rib development and patterning in the mouse embryo has not been resolved, largely because different experimental approaches have yielded disparate results. In this study, we utilize both gene knockouts and muscle cell ablation approaches to re-visit the extent to which rib growth and patterning are dependent on developing musculature. Consistent with previous studies, we show that rib formation is highly dependent on the MYOD family of myogenic regulatory factors (MRFs), and demonstrate that the extent of rib formation is gene-, allele-, and dosage-dependent. In the absence of Myf5 and MyoD, one allele of Mrf4 is sufficient for extensive rib growth, although patterning is abnormal. Under conditions of limiting MRF dosage, MyoD is identified as a positive regulator of rib patterning, presumably due to improved intercostal muscle development. In contrast to previous muscle ablation studies, we show that diphtheria toxin subunit A (DTA)-mediated ablation of muscle progenitors or differentiated muscle, using MyoDiCre or HSA-Cre drivers, respectively, profoundly disrupts rib development. Further, a comparison of three independently derived Rosa26-based DTA knockin alleles demonstrates that the degree of rib perturbations in MyoDiCre/+/DTA embryos is markedly dependent on the DTA allele used, and may in part explain discrepancies with previous findings. The results support the conclusion that the extent and quality of rib formation is largely dependent on the dosage of Myf5 and Mrf4, and that both early myotome-sclerotome interactions, as well as later muscle-rib interactions, are important for proper rib growth and patterning.

Development of the amniote ventrolateral body wall

In vertebrates, the trunk consists of the musculoskeletal structures of the back and the ventrolateral body wall, which together enclose the internal organs of the circulatory, digestive, respiratory and urogenital systems. This review gives an overview on the development of the thoracic and abdominal wall during amniote embryogenesis. Specifically, I briefly summarize relevant historical concepts and the present knowledge on the early embryonic development of ribs, sternum, intercostal muscles and abdominal muscles with respect to anatomical bauplan, origin and specification of precursor cells, initial steps of pattern formation, and cellular and molecular regulation of morphogenesis.

Further evidence for paternal DNA transmission in gynogenetic grass carp

Gynogenesis is an important breeding method in aquaculture and has been widely applied to many fish species. If gynogenetic progenies are to inherit paternal partial genomic DNA, this will increase genetic variation and will provide a useful outcome for breeding. In this study, we investigated the genetic variation in homeobox (Hox) gene clusters (HoxA4a, HoxA9a, HoxA11b, HoxB1b, HoxC4a, HoxC6b, and HoxD10a) among koi carp (Cyprinus carpio haematopterus, KOC; the stimulation sperm source), grass carp (Ctenopharyngodon idellus), and gynogenetic grass carp (GGC). We found paternal DNA (a special DNA fragment and HoxC6b) derived from KOC and a recombinant gene belonging to HoxC6b in GGC. We are the first to report the recombinant HoxC6b in GGC. Our study provides further evidence for paternal DNA transmission to gynogenetic progenies, which is a finding with great significance for the genetic breeding of fish.

Dysbacteriosis-induced LPS elevation disturbs the development of muscle progenitor cells by interfering with retinoic acid signaling

Whether myogenesis is affected by the maternal gut dysbacteriosis still remains ambiguous. In this study, first we show the elevated level of lipopolysaccharides (LPS) in a gut microbiota dysbiosis mouse model. Second, we demonstrate that the diameter of muscle fibers, limb development, and somitogenesis were inhibited in both the gut microbiota dysbiosis and LPS exposed mice and chicken embryos. These might be due to LPS disturbed the cell survival and key genes which regulate the somitogenesis and myogenesis. RNA sequencing and subsequent validation experiments verified that retinoic acid (RA) signaling perturbation was mainly responsible for the aberrant somite formation and differentiation. Subsequently, we found that LPS-induced reactive oxygen species (ROS generation and antioxidant genes such as Nrf2, AKR1B10) contributed to the above -mentioned interference with RA signaling. These findings highlight that the gut microbiota homeostasis is also involved in regulating the development of muscle progenitor cells during pregnancy.

Insight into unique somitogenesis of yak (Bos grunniens) with one additional thoracic vertebra

Background

The yak is a species of livestock which is crucial for local communities of the Qinghai-Tibet Plateau and adjacent regions and naturally owns one more thoracic vertebra than cattle. Recently, a sub-population of yak termed as the Jinchuan yak has been identified with over half its members own a thoracolumbar vertebral formula of T15L5 instead of the natural T14L5 arrangement. The novel T15L5 positioning is a preferred genetic trait leading to enhanced meat and milk production. Selective breeding of this trait would have great agricultural value and exploration of the molecular mechanisms underlying this trait would both accelerate this process and provide us insight into the development and regulation of somitogenesis.

Results

Here we investigated the genetic background of the Jinchuan yak through resequencing fifteen individuals, comprising five T15L5 individuals and ten T14L5 individuals with an average sequencing depth of > 10X, whose thoracolumbar vertebral formulae were confirmed by anatomical observation. Principal component analysis, linkage disequilibrium analysis, phylogenetic analysis, and selective sweep analysis were carried out to explore Jinchuan yak’s genetic background. Three hundred and thirty candidate markers were identified as associated with the additional thoracic vertebrae and target sequencing was used to validate seven carefully selected markers in an additional 51 Jinchuan yaks. The accuracies of predicting 15 thoracic vertebrae and 20 thoracolumbar vertebrae with these 7 markers were 100.00 and 33.33% despite they both could only represent 20% of all possible genetic diversity. Two genes, PPP2R2B and TBLR1, were found to harbour the most candidate markers associated with the trait and likely contribute to the unique somitic number and identity according to their reported roles in the mechanism of somitogenesis.

Conclusions

Our findings provide a clear depiction of the Jinchuan yak’s genetic background and a solid foundation for marker-assistant selection. Further exploitation of this unique population and trait could be promoted with the aid of our genomic resource.

Sacral Dysmorphism and Lumbosacral Transitional Vertebrae (LSTV) Review

Background

Anatomic variation in the relationship between the lumbar spine and sacrum was first described in the literature nearly a century ago and continues to play an important role in spine deformity, low back pain (LBP), and pelvic trauma. This review will focus on the clinical and surgical implications of abnormal lumbosacral anatomy in the context of sacroiliac joint (SIJ) disease, spine deformity, and pelvic trauma.

Methods

A PubMed search using the keywords "lumbosacral transitional vertebrae," "LSTV," "transitional lumbosacral vertebrae," "TLSV," and "sacral dysmorphism" was performed. The articles presented here were evaluated by the authors.

Clinical Significance

The prevalence of LSTV varies widely in the literature from 3.9-% to 35.6% in the spine literature, and sacral dysmorphism is described in upwards of 50% of the population in the trauma literature. The relationship between LSTV and LBP is well established. While there is no agreed-on etiology, the source of pain is multifactorial and may be related to abnormal biomechanics and alignment, disc degeneration, and arthritic changes.

Surgical Implications

Understanding abnormal lumbosacral anatomy is crucial for preoperative planning of SIJ fusion, spine deformity, and pelvic trauma surgery. LSTV can alter spinopelvic parameters crucial in planning spine deformity correction. Traditional safe zones for sacroiliac screw placement do not apply in the first sacral segment in sacral dysmorphism and risk iatrogenic nerve injury.

Conclusions

LSTV and sacral dysmorphism are common anatomic variants found in the general population. Abnormal lumbosacral anatomy plays a significant role in clinical evaluation of LBP and surgical planning in SIJ fusion, spine deformity, and pelvic trauma. Further studies evaluating the influence of abnormal lumbosacral anatomy on LBP and surgical technique would help guide treatment for these patients.

Some Questions and Answers About the Role of Hox Temporal Collinearity in Vertebrate Axial Patterning

The vertebrate anterior-posterior (A-P = craniocaudal) axis is evidently made by a timing mechanism. Evidence has accumulated that tentatively identifies the A-P timer as being or involving Hox temporal collinearity (TC). Here, I focus on the two current competing models based on this premise. Common features and points of dissent are examined and a common model is distilled from what remains. This is an attempt to make sense of the literature.