Linking Genes to Shape in Plants Using Morphometrics.

A transition from qualitative to quantitative descriptors of morphology has been facilitated through the growing field of morphometrics, representing the conversion of shapes and patterns into numbers. The analysis of plant form at the macromorphological scale using morphometric approaches quantifies what is commonly referred to as a phenotype. Quantitative phenotypic analysis of individuals with contrasting genotypes in turn provides a means to establish links between genes and shapes. The path from a gene to a morphological phenotype is, however, not direct, with instructive information progressing both across multiple scales of biological complexity and through nonintuitive feedback, such as mechanical signals. In this review, we explore morphometric approaches used to perform whole-plant phenotyping and quantitative approaches in capture processes in the mesoscales, which bridge the gaps between genes and shapes in plants. Quantitative frameworks involving both the computational simulation and the discretization of data into networks provide a putative path to predicting emergent shape from underlying genetic programs.

An interactive atlas of three-dimensional syndromic facial morphology.

Craniofacial phenotyping is critical for both syndrome delineation and diagnosis because craniofacial abnormalities occur in 30% of characterized genetic syndromes. Clinical reports, textbooks, and available software tools typically provide two-dimensional, static images and illustrations of the characteristic phenotypes of genetic syndromes. In this work, we provide an interactive web application that provides three-dimensional, dynamic visualizations for the characteristic craniofacial effects of 95 syndromes. Users can visualize syndrome facial appearance estimates quantified from data and easily compare craniofacial phenotypes of different syndromes. Our application also provides a map of morphological similarity between a target syndrome and other syndromes. Finally, users can upload 3D facial scans of individuals and compare them to our syndrome atlas estimates. In summary, we provide an interactive reference for the craniofacial phenotypes of syndromes that allows for precise, individual-specific comparisons of dysmorphology.

Craniofacial dysmorphology in Down syndrome is caused by increased dosage of Dyrk1a and at least three other genes.

Down syndrome (DS), trisomy of human chromosome 21 (Hsa21), occurs in 1 in 800 live births and is the most common human aneuploidy. DS results in multiple phenotypes, including craniofacial dysmorphology, which is characterised by midfacial hypoplasia, brachycephaly and micrognathia. The genetic and developmental causes of this are poorly understood. Using morphometric analysis of the Dp1Tyb mouse model of DS and an associated mouse genetic mapping panel, we demonstrate that four Hsa21-orthologous regions of mouse chromosome 16 contain dosage-sensitive genes that cause the DS craniofacial phenotype, and identify one of these causative genes as Dyrk1a. We show that the earliest and most severe defects in Dp1Tyb skulls are in bones of neural crest (NC) origin, and that mineralisation of the Dp1Tyb skull base synchondroses is aberrant. Furthermore, we show that increased dosage of Dyrk1a results in decreased NC cell proliferation and a decrease in size and cellularity of the NC-derived frontal bone primordia. Thus, DS craniofacial dysmorphology is caused by an increased dosage of Dyrk1a and at least three other genes.

The reconstructed cranium of Pierolapithecus and the evolution of the great ape face.

Pierolapithecus catalaunicus (~12 million years ago, northeastern Spain) is key to understanding the mosaic nature of hominid (great ape and human) evolution. Notably, its skeleton indicates that an orthograde (upright) body plan preceded suspensory adaptations in hominid evolution. However, there is ongoing debate about this species, partly because the sole known cranium, preserving a nearly complete face, suffers from taphonomic damage. We 1) carried out a micro computerized tomography (CT) based virtual reconstruction of the Pierolapithecus cranium, 2) assessed its morphological affinities using a series of two-dimensional (2D) and three-dimensional (3D) morphometric analyses, and 3) modeled the evolution of key aspects of ape face form. The reconstruction clarifies many aspects of the facial morphology of Pierolapithecus. Our results indicate that it is most similar to great apes (fossil and extant) in overall face shape and size and is morphologically distinct from other Middle Miocene apes. Crown great apes can be distinguished from other taxa in several facial metrics (e.g., low midfacial prognathism, relatively tall faces) and only some of these features are found in Pierolapithecus, which is most consistent with a stem (basal) hominid position. The inferred morphology at all ancestral nodes within the hominoid (ape and human) tree is closer to great apes than to hylobatids (gibbons and siamangs), which are convergent with other smaller anthropoids. Our analyses support a hominid ancestor that was distinct from all extant and fossil hominids in overall facial shape and shared many features with Pierolapithecus. This reconstructed ancestral morphotype represents a testable hypothesis that can be reevaluated as new fossils are discovered.

Untangling the relationship between developmental and evolutionary integration.

Patterns of integration and modularity among organismal traits are prevalent across the tree of life, and at multiple scales of biological organization. Over the past several decades, researchers have studied these patterns at the developmental, and evolutionary levels. While their work has identified the potential drivers of these patterns at different scales, there appears to be a lack of consensus on the relationship between developmental and evolutionary integration. Here, we review and summarize key studies and build a framework to describe the conceptual relationship between these patterns across organismal scales and illustrate how, and why some of these studies may have yielded seemingly conflicting outcomes. We find that among studies that analyze patterns of integration and modularity using morphological data, the lack of consensus may stem in part from the difficulty of fully disentangling the developmental and functional causes of integration. Nonetheless, in some empirical systems, patterns of evolutionary modularity have been found to coincide with expectations based on developmental processes, suggesting that in some circumstances, developmental modularity may translate to evolutionary modularity. We also advance an extension to Hallgrímsson et al.'s palimpsest model to describe how patterns of trait modularity may shift across different evolutionary scales. Finally, we also propose some directions for future research which will hopefully be useful for investigators interested in these issues.

Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy.

South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.

Dental data challenge the ubiquitous presence of Homo in the Cradle of Humankind.

The origins of Homo, as well as the diversity and biogeographic distribution of early Homo species, remain critical outstanding issues in paleoanthropology. Debates about the recognition of early Homo, first appearance dates, and taxonomic diversity within Homo are particularly important for determining the role that southern African taxa may have played in the origins of the genus. The correct identification of Homo remains also has implications for reconstructing phylogenetic relationships between species of Australopithecus and Paranthropus, and the links between early Homo species and Homo erectus. We use microcomputed tomography and landmark-free deformation-based three-dimensional geometric morphometrics to extract taxonomically informative data from the internal structure of postcanine teeth attributed to Early Pleistocene Homo in the southern African hominin-bearing sites of Sterkfontein, Swartkrans, Drimolen, and Kromdraai B. Our results indicate that, from our sample of 23 specimens, only 4 are unambiguously attributed to Homo, 3 of them coming from Swartkrans member 1 (SK 27, SK 847, and SKX 21204) and 1 from Sterkfontein (Sts 9). Three other specimens from Sterkfontein (StW 80 and 81, SE 1508, and StW 669) approximate the Homo condition in terms of overall enamel-dentine junction shape, but retain Australopithecus-like dental traits, and their generic status remains unclear. The other specimens, including SK 15, present a dominant australopith dental signature. In light of these results, previous dietary and ecological interpretations can be reevaluated, showing that the geochemical signal of one tooth from Kromdraai (KB 5223) and two from Swartkrans (SK 96 and SKX 268) is consistent with that of australopiths.

Factors associated with engraftment success of patient-derived xenografts of breast cancer.

BACKGROUND: Patient-derived xenograft (PDX) models serve as a valuable tool for the preclinical evaluation of novel therapies. They closely replicate the genetic, phenotypic, and histopathological characteristics of primary breast tumors. Despite their promise, the rate of successful PDX engraftment is various in the literature. This study aimed to identify the key factors associated with successful PDX engraftment of primary breast cancer. METHODS: We integrated clinicopathological data with morphological attributes quantified using a trained artificial intelligence (AI) model to identify the principal factors affecting PDX engraftment. RESULTS: Multivariate logistic regression analyses demonstrated that several factors, including a high Ki-67 labeling index (Ki-67LI) (p < 0.001), younger age at diagnosis (p = 0.032), post neoadjuvant chemotherapy (NAC) (p = 0.006), higher histologic grade (p = 0.039), larger tumor size (p = 0.029), and AI-assessed higher intratumoral necrosis (p = 0.027) and intratumoral invasive carcinoma (p = 0.040) proportions, were significant factors for successful PDX engraftment (area under the curve [AUC] 0.905). In the NAC group, a higher Ki-67LI (p < 0.001), lower Miller-Payne grade (p < 0.001), and reduced proportion of intratumoral normal breast glands as assessed by AI (p = 0.06) collectively provided excellent prediction accuracy for successful PDX engraftment (AUC 0.89). CONCLUSIONS: We found that high Ki-67LI, younger age, post-NAC status, higher histologic grade, larger tumor size, and specific morphological attributes were significant factors for predicting successful PDX engraftment of primary breast cancer.

Utilizing geometric morphometrics to investigate gene function during organ growth: Insights through the study of beetle horn shape allometry.

Static allometry is a major component of morphological variation. Much of the literature on the development of allometry investigates how functional perturbations of diverse pathways affect the relationship between trait size and body size. Often, this is done with the explicit objective to identify developmental mechanisms that enable the sensing of organ size and the regulation of relative growth. However, changes in relative trait size can also be brought about by a range of other distinctly different developmental processes, such as changes in patterning or tissue folding, yet standard univariate biometric approaches are usually unable to distinguish among alternative explanations. Here, we utilize geometric morphometrics to investigate the degree to which functional genetic manipulations known to affect the size of dung beetle horns also recapitulate the effect of horn shape allometry. We reasoned that the knockdown phenotypes of pathways governing relative growth should closely resemble shape variation induced by natural allometric variation. In contrast, we predicted that if genes primarily affect alternative developmental processes, knockdown effects should align poorly with shape allometry. We find that the knockdown effects of several genes (e.g., doublesex, Foxo) indeed closely aligned with shape allometry, indicating that their corresponding pathways may indeed function primarily in the regulation of relative trait growth. In contrast, other knockdown effects (e.g., Distal-less, dachs) failed to align with allometry, implicating these pathways in potentially scaling-independent processes. Our findings moderate the interpretation of studies focusing on trait length and highlight the usefulness of multivariate approaches to study allometry and phenotypic plasticity.

Got milkweed? Genetic assimilation as potential source for the evolution of nonmigratory monarch butterfly wing shape.

Monarch butterflies (Danaus plexippus) are well studied for their annual long-distance migration from as far north as Canada to their overwintering grounds in Central Mexico. At the end of the cold season, monarchs start to repopulate North America through short-distance migration over the course of multiple generations. Interestingly, some populations in various tropical and subtropical islands do not migrate and exhibit heritable differences in wing shape and size, most likely an adaptation to island life. Less is known about forewing differences between long- and short-distance migrants in relation to island populations. Given their different migratory behaviors, we hypothesized that these differences would be reflected in wing morphology. To test this, we analyzed forewing shape and size of three different groups: nonmigratory, lesser migratory (migrate short-distances), and migratory (migrate long-distances) individuals. Significant differences in shape appear in all groups using geometric morphometrics. As variation found between migratory and lesser migrants has been shown to be caused by phenotypic plasticity, and lesser migrants develop intermediate forewing shapes between migratory and nonmigratory individuals, we suggest that genetic assimilation might be an important mechanism to explain the heritable variation found between migratory and nonmigratory populations. Additionally, our research confirms previous studies which show that forewing size is significantly smaller in nonmigratory populations when compared to both migratory phenotypes. Finally, we found sexual dimorphism in forewing shape in all three groups, but for size in nonmigratory populations only. This might have been caused by reduced constraints on forewing size in nonmigratory populations.

Replicated Functional Evolution in Cichlid Adaptive Radiations.

AbstractAdaptive radiations highlight the mechanisms by which species and traits diversify and the extent to which these patterns are predictable. We used 1,110 high-speed videos of suction feeding to study functional and morphological diversification in 300 cichlid species from three African Great Lake radiations of varying ages (Victoria, Malawi, and Tanganyika) and an older, spatially dispersed continental radiation in the Neotropics. Among African radiations, standing diversity was reflective of time. Morphological and functional variance in Lake Victoria, the youngest radiation, was a subset of that within Lake Malawi, which itself was nested within the older Tanganyikan radiation. However, functional diversity in Neotropical cichlids was often lower than that in Lake Tanganyika, despite being much older. These two radiations broadly overlapped, but each diversified into novel trait spaces not found in the youngest lake radiations. Evolutionary rates across radiations were inversely related to age, suggesting extremely rapid trait evolution at early stages, particularly in lake radiations. Despite this support for early bursts, other patterns of trait diversity were inconsistent with expectations of adaptive radiations. This work suggests that cichlid functional evolution has played out in strikingly similar fashion in different radiations, with contingencies eventually resulting in lineage-specific novelties.

MRI-based axis-referenced morphometric model corresponding to lamellar organization for assessing hippocampal atrophy in dementia.

Research on the local hippocampal atrophy for early detection of dementia has gained considerable attention. However, accurately quantifying subtle atrophy remains challenging in existing morphological methods due to the lack of consistent biological correspondence with the complex curving regions like the hippocampal head. Thereby, this article presents an innovative axis-referenced morphometric model (ARMM) that follows the anatomical lamellar organization of the hippocampus, which capture its precise and consistent longitudinal curving trajectory. Specifically, we establish an "axis-referenced coordinate system" based on a 7 T ex vivo hippocampal atlas following its entire curving longitudinal axis and orthogonal distributed lamellae. We then align individual hippocampi by deforming this template coordinate system to target spaces using boundary-guided diffeomorphic transformation, while ensuring that the lamellar vectors adhere to the constraint of medial-axis geometry. Finally, we measure local thickness and curvatures based on the coordinate system and boundary surface reconstructed from vector tips. The morphometric accuracy is evaluated by comparing reconstructed surfaces with those directly extracted from 7 T and 3 T MRI hippocampi. The results demonstrate that ARMM achieves the best performance, particularly in the curving head, surpassing the state-of-the-art morphological models. Additionally, morphological measurements from ARMM exhibit higher discriminatory power in distinguishing early Alzheimer's disease from mild cognitive impairment compared to volume-based measurements. Overall, the ARMM offers a precise morphometric assessment of hippocampal morphology on MR images, and sheds light on discovering potential image markers for neurodegeneration associated with hippocampal impairment.

A landmark-free morphometrics pipeline for high-resolution phenotyping: application to a mouse model of Down syndrome.

Characterising phenotypes often requires quantification of anatomical shape. Quantitative shape comparison (morphometrics) traditionally uses manually located landmarks and is limited by landmark number and operator accuracy. Here, we apply a landmark-free method to characterise the craniofacial skeletal phenotype of the Dp1Tyb mouse model of Down syndrome and a population of the Diversity Outbred (DO) mouse model, comparing it with a landmark-based approach. We identified cranial dysmorphologies in Dp1Tyb mice, especially smaller size and brachycephaly (front-back shortening), homologous to the human phenotype. Shape variation in the DO mice was partly attributable to allometry (size-dependent shape variation) and sexual dimorphism. The landmark-free method performed as well as, or better than, the landmark-based method but was less labour-intensive, required less user training and, uniquely, enabled fine mapping of local differences as planar expansion or shrinkage. Its higher resolution pinpointed reductions in interior mid-snout structures and occipital bones in both the models that were not otherwise apparent. We propose that this landmark-free pipeline could make morphometrics widely accessible beyond its traditional niches in zoology and palaeontology, especially in characterising developmental mutant phenotypes.

Single-cell morphometrics reveals ancestral principles of notochord development.

Embryonic tissues are shaped by the dynamic behaviours of their constituent cells. To understand such cell behaviours and how they evolved, new approaches are needed to map out morphogenesis across different organisms. Here, we apply a quantitative approach to learn how the notochord forms during the development of amphioxus: a basally branching chordate. Using a single-cell morphometrics pipeline, we quantify the geometries of thousands of amphioxus notochord cells, and project them into a common mathematical space, termed morphospace. In morphospace, notochord cells disperse into branching trajectories of cell shape change, revealing a dynamic interplay between cell shape change and growth that collectively contributes to tissue elongation. By spatially mapping these trajectories, we identify conspicuous regional variation, both in developmental timing and trajectory topology. Finally, we show experimentally that, unlike ascidians but like vertebrates, posterior cell division is required in amphioxus to generate full notochord length, thereby suggesting this might be an ancestral chordate trait that is secondarily lost in ascidians. Altogether, our novel approach reveals that an unexpectedly complex scheme of notochord morphogenesis might have been present in the first chordates. This article has an associated 'The people behind the papers' interview.

A shared pattern of midfacial bone modelling in hominids suggests deep evolutionary roots for human facial morphogenesis.

Midfacial morphology varies between hominoids, in particular between great apes and humans for which the face is small and retracted. The underlying developmental processes for these morphological differences are still largely unknown. Here, we investigate the cellular mechanism of maxillary development (bone modelling, BM), and how potential changes in this process may have shaped facial evolution. We analysed cross-sectional developmental series of gibbons, orangutans, gorillas, chimpanzees and present-day humans (n = 183). Individuals were organized into five age groups according to their dental development. To visualize each species's BM pattern and corresponding morphology during ontogeny, maps based on microscopic data were mapped onto species-specific age group average shapes obtained using geometric morphometrics. The amount of bone resorption was quantified and compared between species. Great apes share a highly similar BM pattern, whereas gibbons have a distinctive resorption pattern. This suggests a change in cellular activity on the hominid branch. Humans possess most of the great ape pattern, but bone resorption is high in the canine area from birth on, suggesting a key role of canine reduction in facial evolution. We also observed that humans have high levels of bone resorption during childhood, a feature not shared with other apes.

Rapid divergent evolution of internal female genitalia and the coevolution of male genital morphology revealed by micro-computed tomography.

Animal genitalia are thought to evolve rapidly and divergently in response to sexual selection. Studies of genital evolution have focused largely on male genitalia. The paucity of work on female genital morphology is probably due to problems faced in quantifying shape variation, due to their composition and accessibility. Here we use a combination of micro-computed tomography, landmark free shape quantification and phylogenetic analysis to quantify the rate of female genital shape evolution among 29 species of Antichiropus millipedes, and their coevolution with male genitalia. We found significant variation in female and male genital shape among species. Male genital shape showed a stronger phylogenetic signal than female genital shape, although the phylogenetic signal effect sizes did not differ significantly. Male genital shape was found to be evolving 1.2 times faster than female genital shape. Female and male genital shape exhibited strong correlated evolution, indicating that genital shape changes in one sex are associated with corresponding changes in the genital shape of the other sex. This study adds novel insight into our growing understanding of how female genitalia can evolve rapidly and divergently, and highlights the advantages of three-dimensional techniques and multivariate analyses in studies of female genital evolution.

Intra-leaf modeling of Cannabis leaflet shape produces leaf models that predict genetic and developmental identities.

The iconic, palmately compound leaves of Cannabis have attracted significant attention in the past. However, investigations into the genetic basis of leaf shape or its connections to phytochemical composition have yielded inconclusive results. This is partly due to prominent changes in leaflet number within a single plant during development, which has so far prevented the proper use of common morphometric techniques. Here, we present a new method that overcomes the challenge of nonhomologous landmarks in palmate, pinnate, and lobed leaves, using Cannabis as an example. We model corresponding pseudo-landmarks for each leaflet as angle-radius coordinates and model them as a function of leaflet to create continuous polynomial models, bypassing the problems associated with variable number of leaflets between leaves. We analyze 341 leaves from 24 individuals from nine Cannabis accessions. Using 3591 pseudo-landmarks in modeled leaves, we accurately predict accession identity, leaflet number, and relative node number. Intra-leaf modeling offers a rapid, cost-effective means of identifying Cannabis accessions, making it a valuable tool for future taxonomic studies, cultivar recognition, and possibly chemical content analysis and sex identification, in addition to permitting the morphometric analysis of leaves in any species with variable numbers of leaflets or lobes.

Phenotypic selection patterns in a hybrid zone between two Calceolaria species with contrasting pollinators: insights from field surveys and fitness assessments.

Hybrid zones provide natural experimental settings to test hypotheses about species divergence. We concentrated on a hybrid swarm in which oil-collecting bees and flower-pecking birds act as pollinators of two Calceolaria species. We asked whether both pollinators contributed to flower divergence by differentially promoting prezygotic fitness at the phenotypic extremes that represent parentals. We studied pollinator-mediated selection on phenotypic traits critical in plant-pollinator mechanical interaction, namely plant height, reward-to-stigma distance, and flower shape. We utilised the quantity and quality of pollen deposited as fitness measures and distinguished between the contribution of the two pollinator types. Results showed uni- and bivariate disruptive selection for most traits through pollen grains deposited by both pollinators. Bird-mediated fitness favoured low plants with a long reward-to-stigma distance and a straight corolla, while bee-mediated fitness favoured tall plants with a short reward-to-stigma distance and curved corolla. In addition, stabilising selection at one end of the phenotypic range showed a bird-mediated reproductive asymmetry within the swarm. The disruptive pattern was countered, albeit weakly, by hybrids receiving higher-quality pollen on the stigmas. Results suggest that pollinator-mediated selection promotes divergence of integrated flower phenotypes mechanically adjusted either to bees or birds underscoring the importance of pollinator specialisation in diversification.

Charting facial growth and development for Bantu Africans: Central tendencies, variational properties and sexual dimorphisms.

Current studies on facial growth and development have been largely based on European populations. Less studied are African populations, who because of their distinct genetic makeup and environmental conditions, provide deeper insights into patterns of facial development. Patterns of facial shape development in African populations remain largely uncharacterised. Our study aimed to establish facial growth and development trajectories based on a cohort of 2874 Bantu Africans from Tanzania aged 6-18 years, with particular focus on identifying morphogenetic processes that lead to observed developmental shape changes. Procrustes ANCOVA suggested sexually dimorphic patterns of facial shape development (p = 0.0036). The forehead was relatively contracted during development in both sexes. The glabella region was more anteriorly displaced in females due to expansion in the region laterosuperior to the eyes. Nasal protrusion increased with development, which was found to arise from local expansion in the nasal alae and columella. Local expansion in the upper and lower labial regions resulted in forward displaced lips in both sexes, with the effect more pronounced in males. The mentum was displaced more anteriorly in females due to comparatively more expanded mental regions with development. The lateral facial region corresponding to the underlying body of the mandible were developmentally expanded but were posteriorly positioned due to protrusive growth of surrounding structures. Generalised additive modelling of Procrustes variance suggested that facial variation decreased non-linearly with age (p < 0.05). Relative principal component analysis suggested that variations in facial outline shape were developmentally constrained, whereas nasolabial and mental regions, where developmental changes were significant, became morphologically diversified with development. In contrast to simple descriptive illustration of facial shape development, we gained transformative insights into patterns of facial shape development by analysing morphogenetic processes and variational properties. Our analytical framework is broadly applicable to morphometric studies on ontogenetic shape changes.

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.