Plastic and genomic change of a newly established lizard population following a founder event.

Understanding how phenotypic divergence arises among natural populations remains one of the major goals in evolutionary biology. As part of competitive exclusion experiment conducted in 1971, 10 individuals of Italian wall lizard (Podarcis siculus (Rafinesque-Schmaltz, 1810)) were transplanted from Pod Kopiste Island to the nearby island of Pod Mrcaru (Adriatic Sea). Merely 35 years after the introduction, the newly established population on Pod Mrcaru Island had shifted their diet from predominantly insectivorous towards omnivorous and changed significantly in a range of morphological, behavioural, physiological and ecological characteristics. Here, we combine genomic and quantitative genetic approaches to determine the relative roles of genetic adaptation and phenotypic plasticity in driving this rapid phenotypic shift. Our results show genome-wide genetic differentiation between ancestral and transplanted population, with weak genetic erosion on Pod Mrcaru Island. Adaptive processes following the founder event are indicated by highly differentiated genomic loci associating with ecologically relevant phenotypic traits, and/or having a putatively adaptive role across multiple lizard populations. Diverged traits related to head size and shape or bite force showed moderate heritability in a crossing experiment, but between-population differences in these traits did not persist in a common garden environment. Our results confirm the existence of sufficient additive genetic variance for traits to evolve under selection while also demonstrating that phenotypic plasticity and/or genotype by environment interactions are the main drivers of population differentiation at this early evolutionary stage.

Geometric growth of the normal human craniocervical junction from 0 to 18 years old.

The craniocervical junction (CCJ) forms the bridge between the skull and the spine, a highly mobile group of joints that allows the mobility of the head in every direction. The CCJ plays a major role in protecting the inferior brainstem (bulb) and spinal cord, therefore also requiring some stability. Children are subjected to multiple constitutive or acquired diseases involving the CCJ: primary bone diseases such as in FGFR-related craniosynostoses or acquired conditions such as congenital torticollis, cervical spine luxation, and neurological disorders. To design efficient treatment plans, it is crucial to understand the relationship between abnormalities of the craniofacial region and abnormalities of the CCJ. This can be approached by the study of control and abnormal growth patterns. Here we report a model of normal skull base growth by compiling a collection of geometric models in control children. Focused analyses highlighted specific developmental patterns for each CCJ bone, emphasizing rapid growth during infancy, followed by varying rates of growth and maturation during childhood and adolescence until reaching stability by 18 years of age. The focus was on the closure patterns of synchondroses and sutures in the occipital bone, revealing distinct closure trajectories for the anterior intra-occipital synchondroses and the occipitomastoid suture. The findings, although based on a limited dataset, showcased specific age-related changes in width and closure percentages, providing valuable insights into growth dynamics within the first 2 years of life. Integration analyses revealed intricate relationships between skull and neck structures, emphasizing coordinated growth at different stages. Specific bone covariation patterns, as found between the first and second cervical vertebrae (C1 and C2), indicated synchronized morphological changes. Our results provide initial data for designing inclusive CCJ geometric models to predict normal and abnormal growth dynamics.

Relationships between dietary breadth and flexibility in jaw movement: A case study of two recently diverged insular populations of Podarcis lizards.

The kinematics of lizard feeding are the result of complex interactions between the craniocervical, the hyolingual, and the locomotor systems. The coordinated movement of these elements is driven by sensory feedback from the tongue and jaws during intraoral transport. The kinematics of jaw movements have been suggested to be correlated with the functional characteristics of the prey consumed, such as prey mobility and hardness. However, whether and how dietary breadth correlates with the flexibility in the behavioral response has rarely been tested, especially at the intraspecific level. Here we tested whether an increase in dietary breadth was associated with a greater behavioral flexibility by comparing two recently diverged populations of insular Podarcis lizards differing in dietary breadth. To do so, we used a stereoscopic high-speed camera set-up to analyze the jaw kinematics while offering them different prey types. Our results show that prey type impacts kinematics, especially maximum gape, and maximum opening and closing speed. Furthermore, the behavioral flexibility was greater in the population with the greater dietary breadth, suggesting that populations which naturally encounter and feed on more diverse prey items show a greater ability to modulate their movements to deal with variation in functionally relevant prey properties. Finally, the more generalist population showed more stereotyped movements suggesting a finer motor control.

3D quantitative comparative analysis of long bone diaphysis variations in microanatomy and cross-sectional geometry.

Long bone inner structure and cross-sectional geometry display a strong functional signal, leading to convergences, and are widely analyzed in comparative anatomy at small and large taxonomic scales. Long bone microanatomical studies have essentially been conducted on transverse sections but also on a few longitudinal ones. Recent studies highlighted the interest in analyzing variations of the inner structure along the diaphysis using a qualitative as well as a quantitative approach. With the development of microtomography, it has become possible to study three-dimensional (3D) bone microanatomy and, in more detail, the form-function relationships of these features. This study focused on the selection of quantitative parameters to describe in detail the cross-sectional shape changes and distribution of the osseous tissue along the diaphysis. Two-dimensional (2D) virtual transverse sections were also performed in the two usual reference planes and results were compared with those obtained based on the whole diaphysis analysis. The sample consisted in 14 humeri and 14 femora of various mammalian taxa that are essentially terrestrial. Comparative quantitative analyses between different datasets made it possible to highlight the parameters that are strongly impacted by size and phylogeny and the redundant ones, and thus to estimate their relevance for use in form-function analyses. The analysis illustrated that results based on 2D transverse sections are similar for both sectional planes; thus if a strong bias exists when mixing sections from the two reference planes in the same analysis, it would not problematic to use either one plane or the other in comparative studies. However, this may no longer hold for taxa showing a much stronger variation in bone microstructure along the diaphysis. Finally, the analysis demonstrated the significant contribution of the parameters describing variations along the diaphysis, and thus the interest in performing 3D analyses; this should be even more fruitful for heterogeneous diaphyses. In addition, covariation analyses showed that there is a strong interest in removing the size effect to access the differences in the microstructure of the humerus and femur. This methodological study provides a reference for future quantitative analyses on long bone inner structure and should make it possible, through a detailed knowledge of each descriptive parameter, to better interpret results from the multivariate analyses associated with these studies. This will have direct implications for studies in vertebrate anatomy, but also in paleontology and anthropology.

New diagnostic criteria for metopic ridges and trigonocephaly: a 3D geometric approach.

BACKGROUND: Trigonocephaly occurs due to the premature fusion of the metopic suture, leading to a triangular forehead and hypotelorism. This condition often requires surgical correction for morphological and functional indications. Metopic ridges also originate from premature metopic closure but are only associated with mid-frontal bulging; their surgical correction is rarely required. Differential diagnosis between these two conditions can be challenging, especially in minor trigonocephaly. METHODS: Two hundred seven scans of patients with trigonocephaly (90), metopic rigdes (27), and controls (90) were collected. Geometric morphometrics were used to quantify skull and orbital morphology as well as the interfrontal angle and the cephalic index. An innovative method was developed to automatically compute the frontal curvature along the metopic suture. Different machine-learning algorithms were tested to assess the predictive power of morphological data in terms of classification. RESULTS: We showed that control patients, trigonocephaly and metopic rigdes have distinctive skull and orbital shapes. The 3D frontal curvature enabled a clear discrimination between groups (sensitivity and specificity > 92%). Furthermore, we reached an accuracy of 100% in group discrimination when combining 6 univariate measures. CONCLUSION: Two diagnostic tools were proposed and demonstrated to be successful in assisting differential diagnosis for patients with trigonocephaly or metopic ridges. Further clinical assessments are required to validate the practical clinical relevance of these tools.

Anatomy and mobility in the adult cadaveric craniocervical junction.

Genetic diseases with craniofacial malformations can be associated with anomalies of the craniocervical joint (CCJ). The functions of the CCJ are thus impaired, as mobility may be either limited by abnormal bone fusion causing headaches, or exaggerated in the case of hypermobility, which may cause irreparable damage to the spinal cord. Restoring the balance between mobility and stability requires surgical correction in children. The anatomy and biomechanics of the CCJ are quite unique, yet have been overlooked in the past decades. Pediatric evidence is so scarce, that investigating the adult CCJ is our best shot to disentangle the form-function relationships of this anatomical region. The motivation of the present study was to understand the morphological and functional basis of motion in the CCJ, in the hope to find morphological features accessible from medical imaging able to predict mobility. To do so, we have quantified the in-vitro kinematics of the CCJ in nine cadaveric asymptomatic adults, and estimated a wide range of mobility variables covering the complexity of spinal motion. We compared these variables with the shape of the occipital, the atlas and the axis, obtained using a dense geometric morphometric approach. Morphological joint congruence was also quantified. Our results suggest a strong relationship between bone shape and motion, with the overall geometry predicting best the primary movements, and the joint facets predicting best the secondary movements. We propose a functional hypothesis stating that the musculoligamental system determines movements of great amplitude, while the shape and congruence of joint facets determine the secondary and coupled movements, especially by varying the geometry of bone stops and the way ligaments are tensioned. We believe this work will provide valuable insights in understanding the biomechanics of the CCJ. Furthermore, it should help surgeons treating CCJ anomalies by enabling them to translate objectives of functional and clinical outcome into clear objectives of morphological outcome.

The Evolution of Diet and Morphology in Insular Lizards: Insights from a Replicated Island Introduction Experiment.

Resource-limited environments may drive the rapid evolution of phenotypic traits and ecological preferences optimizing the exploitation of resources. Very small islands are often characterized by reduced food availability, seasonal fluctuations in resources and strong unpredictability. These features may drive the evolution of phenotypic traits such as high bite forces, allowing animals to exploit a wider variety of the available resources. They may also lead to more generalist dietary patterns in response to food scarcity. However, the lack of predators and competitors on such small islands often also leads to high densities and the evolution of strong sexual dimorphism, which may also drive the evolution of bite force. Here, we take advantage of a unique replicated introduction experiment to test whether lizards introduced into very small islands alter their feeding ecology and use different resources, resulting in the evolution of a large body size, large head size and large bite forces. Our results show that three years after their introduction, the island lizards were larger and had greater bite forces and more pronounced sexual dimorphism. However, the diets were only marginally different between animals from the source population on a very large nearby island and those on the islets. Moreover, distinct differences in diet between animals on the different islets were observed, suggesting that the local environment is a strong driver of resource use. Overall, lizards with absolutely and relatively (adjusted for body size) large bite forces did eat larger and harder prey. Taken together, our data suggest that intraspecific competition is an important driver of the rapid evolution of bite force, which may allow these lizards to exploit the scarce and fluctuating resources on the islets. Whether or not lizards will evolve to include other types of food such as plants in their diet, facilitated by their large bite forces, remains to be explored in future studies.

Characterization of labiomandibular movements induced after isolated LeFort I osteotomy in the surgical management of class III malocclusion.

INTRODUCTION: Maxillary surgery alone can be proposed for the surgical management of class III malocclusion, but anticipating outcomes for the labiomental muscle complex is challenging due to the mandibular autorotation phenomenon. The objective of this study was to quantify the mandibular and labiomental movements induced by maxillary osteotomy alone in the management of class III malocclusion according to different clinical and surgical variables. METHODS: The post-operative changes in mandibular and labiomental shapes were studied by geometric morphometry from the pre- and post-operative lateral cephalometric radiograph of 25 patients. The explanatory variables tested were maxillary advancement, maxillary rotation, and divergence. RESULTS: Soft tissues repositioning are different from postoperative mandibular repositioning after maxillary osteotomy. Neuromuscular adjustments of mandible depend on divergence and the maxillary rotation. Labiomental response only depends on divergence. CONCLUSION: The surgical procedure does not have the same bone-related and musculocutaneous effects on patients with the same class III malocclusion. It is therefore essential for surgeons to understand the effects of their procedure on musculocutaneous tissues in order to best anticipate post-operative outcomes.

Assessing craniofacial growth and form without landmarks: A new automatic approach based on spectral methods.

We present a novel method for the morphometric analysis of series of 3D shapes, and demonstrate its relevance for the detection and quantification of two craniofacial anomalies: trigonocephaly and metopic ridges, using CT-scans of young children. Our approach is fully automatic, and does not rely on manual landmark placement and annotations. Our approach furthermore allows to differentiate shape classes, enabling successful differential diagnosis between trigonocephaly and metopic ridges, two related conditions characterized by triangular foreheads. These results were obtained using recent developments in automatic nonrigid 3D shape correspondence methods and specifically spectral approaches based on the functional map framework. Our method can capture local changes in geometric structure, in contrast to methods based, for instance, on global shape descriptors. As such, our approach allows to perform automatic shape classification and provides visual feedback on shape regions associated with different classes of deformations. The flexibility and generality of our approach paves the way for the application of spectral methods in quantitative medicine.

Craniofacial growth and function in achondroplasia: a multimodal 3D study on 15 patients.

BACKGROUND: Achondroplasia is the most frequent FGFR3-related chondrodysplasia, leading to rhizomelic dwarfism, craniofacial anomalies, stenosis of the foramen magnum, and sleep apnea. Craniofacial growth and its correlation with obstructive sleep apnea syndrome has not been assessed in achondroplasia. In this study, we provide a multimodal analysis of craniofacial growth and anatomo-functional correlations between craniofacial features and the severity of obstructive sleep apnea syndrome. METHODS: A multimodal study was performed based on a paediatric cohort of 15 achondroplasia patients (mean age, 7.8 ± 3.3 years), including clinical and sleep study data, 2D cephalometrics, and 3D geometric morphometry analyses, based on CT-scans (mean age at CT-scan: patients, 4.9 ± 4.9 years; controls, 3.7 ± 4.2 years). RESULTS: Craniofacial phenotype was characterized by maxillo-zygomatic retrusion, deep nasal root, and prominent forehead. 2D cephalometric studies showed constant maxillo-mandibular retrusion, with excessive vertical dimensions of the lower third of the face, and modifications of cranial base angles. All patients with available CT-scan had premature fusion of skull base synchondroses. 3D morphometric analyses showed more severe craniofacial phenotypes associated with increasing patient age, predominantly regarding the midface-with increased maxillary retrusion in older patients-and the skull base-with closure of the spheno-occipital angle. At the mandibular level, both the corpus and ramus showed shape modifications with age, with shortened anteroposterior mandibular length, as well as ramus and condylar region lengths. We report a significant correlation between the severity of maxillo-mandibular retrusion and obstructive sleep apnea syndrome (p < 0.01). CONCLUSIONS: Our study shows more severe craniofacial phenotypes at older ages, with increased maxillomandibular retrusion, and demonstrates a significant anatomo-functional correlation between the severity of midface and mandible craniofacial features and obstructive sleep apnea syndrome.

Mobility Assessment Using Multi-Positional MRI in Children with Cranio-Vertebral Junction Anomalies.

OBJECTIVE: This study aimed to assess the relevance of using multi-positional MRI (mMRI) to identify cranio-vertebral junction (CVJ) instability in pediatric patients with CVJ anomalies while determining objective mMRI criteria to detect this condition. MATERIAL AND METHODS: Data from children with CVJ anomalies who underwent a mMRI between 2017 and 2021 were retrospectively reviewed. Mobility assessment using mMRI involved: (1) morphometric analysis using hierarchical clustering on principal component analysis (HCPCA) to identify clusters of patients by considering their mobility similarities, assessed through delta (Δ) values of occipito-cervical parameters measured on mMRI; and (2) morphological analysis based on dynamic geometric CVJ models and analysis of displacement vectors between flexion and extension. Receiver operating characteristics (ROC) curves were generated for occipito-cervical parameters to establish instability cut-off values. (3) Additionally, an anatomical qualitative analysis of the CVJ was performed to identify morphological criteria of instability. RESULTS: Forty-seven patients with CVJ anomalies were included (26 females, 21 males; mean age: 10.2 years [3-18]). HCPCA identified 2 clusters: cluster №1 (stable patients, n = 39) and cluster №2 (unstable patients, n = 8). ΔpB-C2 (pB-C2 line delta) at ≥2.5 mm (AUC 0.98) and ΔBAI (Basion-axis Interval delta) ≥ 3 mm (AUC 0.97) predicted instability with 88% sensibility and 95% specificity and 88% sensitivity and 85% specificity, respectively. Geometric CVJ shape analysis differentiated patients along a continuum, from a low to a high CVJ motion that was characterized by a subluxation of C1 in the anterior direction. Qualitative analysis found correlations between instability and C2 anomalies, including fusions with C3 (body p = 0.032; posterior arch p = 0.045; inferior articular facets p = 0.012; lateral mass p = 0.029). CONCLUSIONS: We identified a cluster of pediatric patients with CVJ instability among a cohort of CVJ anomalies that were characterized by morphometric parameters with corresponding cut-off values that could serve as objective mMRI criteria. These findings warrant further validation through prospective case-control studies.

Hemifacial myohyperplasia is due to somatic muscular PIK3CA gain-of-function mutations and responds to pharmacological inhibition.

Hemifacial myohyperplasia (HFMH) is a rare cause of facial asymmetry exclusively involving facial muscles. The underlying cause and the mechanism of disease progression are unknown. Here, we identified a somatic gain-of-function mutation of PIK3CA in five pediatric patients with HFMH. To understand the physiopathology of muscle hypertrophy in this context, we created a mouse model carrying specifically a PIK3CA mutation in skeletal muscles. PIK3CA gain-of-function mutation led to striated muscle cell hypertrophy, mitochondria dysfunction, and hypoglycemia with low circulating insulin levels. Alpelisib treatment, an approved PIK3CA inhibitor, was able to prevent and reduce muscle hypertrophy in the mouse model with correction of endocrine anomalies. Based on these findings, we treated the five HFMH patients. All patients demonstrated clinical, esthetical, and radiological improvement with proof of target engagement. In conclusion, we show that HFMH is due to somatic alteration of PIK3CA and is accessible to pharmacological intervention.

Chiari malformation type I surgery in children: French multicenter 10-year cohort.

OBJECTIVE: Chiari malformation type I (CM-I) is frequent in children and remains a surgical challenge. Several techniques have been described for posterior fossa decompression. No decision algorithm has been validated, and strategies are highly variable between institutions. The goal of this study was to define therapeutic guidelines that take into consideration patient specificities. METHODS: The authors retrospectively collected data from patients who were < 18 years of age, were diagnosed with CM-I, and were treated surgically between 2008 and 2018 in 8 French pediatric neurosurgical centers. Data on clinical features, morphological parameters, and surgical techniques were collected. Clinical outcomes at 3 and 12 months after surgery were assessed by the Chicago Chiari Outcome Scale. The authors used a hierarchical clustering method to define clusters of patients by considering their anatomical similarities, and then compared outcomes between surgical strategies in each of these clusters. RESULTS: Data from 255 patients were collected. The mean age at surgery was 9.6 ± 5.0 years, syringomyelia was reported in 60.2% of patients, the dura mater was opened in 65.0% of patients, and 17.3% of patients underwent a redo surgery for additional treatment. The mean Chicago Chiari Outcome Scale score was 14.4 ± 1.5 at 3 months (n = 211) and 14.6 ± 1.9 at 12 months (n = 157). The hierarchical clustering method identified three subgroups with potentially distinct mechanisms underlying tonsillar herniation: bony compression, basilar invagination, and foramen magnum obstruction. Each cluster matched with specific outcomes. CONCLUSIONS: This French multicenter retrospective cohort study enabled the identification of three subgroups among pediatric patients who underwent surgery for CM-I, each of which was associated with specific outcomes. This morphological classification of patients might help in understanding the underlying mechanisms and providing personalized treatment.

From micro to macroevolution: drivers of shape variation in an island radiation of Podarcis lizards.

Phenotypictraits have been shown to evolve in response to variation in the environment. However, the evolutionary processes underlying the emergence of phenotypic diversity can typically only be understood at the population level. Consequently, how subtle phenotypic differences at the intraspecific level can give rise to larger-scale changes in performance and ecology remains poorly understood. We here tested for the covariation between ecology, bite force, jaw muscle architecture, and the three-dimensional shape of the cranium and mandible in 16 insular populations of the lizards Podarcis melisellensis and P. sicula. We then compared the patterns observed at the among-population level with those observed at the interspecific level. We found that three-dimensional head shape as well as jaw musculature evolve similarly under similar ecological circumstances. Depending on the type of food consumed or on the level of sexual competition, different muscle groups were more developed and appeared to underlie changes in cranium and mandible shape. Our findings show that the local selective regimes are primary drivers of phenotypic variation resulting in predictable patterns of form and function. Moreover, intraspecific patterns of variation were generally consistent with those at the interspecific level, suggesting that microevolutionary variation may translate into macroevolutionary patterns of ecomorphological diversity.

Diet variability among insular populations of Podarcis lizards reveals diverse strategies to face resource-limited environments.

Access to resources is a dynamic and multicausal process that determines the success and survival of a population. It is therefore often challenging to disentangle the factors affecting ecological traits like diet. Insular habitats provide a good opportunity to study how variation in diet originates, in particular in populations of mesopredators such as lizards. Indeed, high levels of population density associated with low food abundance and low predation are selection pressures typically observed on islands. In the present study, the diet of eighteen insular populations of two closely related species of lacertid lizards (Podarcis sicula and Podarcis melisellensis) was assessed. Our results reveal that despite dietary variability among populations, diet taxonomic diversity is not impacted by island area. In contrast, however, diet disparity metrics, based on the variability in the physical (hardness) and behavioral (evasiveness) properties of ingested food items, are correlated with island size. These findings suggest that an increase in intraspecific competition for access to resources may induce shifts in functional components of the diet. Additionally, the two species differed in the relation between diet disparity and island area suggesting that different strategies exist to deal with low food abundance in these two species. Finally, sexual dimorphism in diet and head dimensions is not greater on smaller islands, in contrast to our predictions.