The Developmental Basis of Quantitative Craniofacial Variation in Humans and Mice

The human skull is a complex and highly integrated structure that has long held the fascination of anthropologists and evolutionary biologists. Recent studies of the genetics of craniofacial variation reveal a very complex and multifactorial picture. These findings contrast with older ideas that posit much simpler developmental bases for variation in cranial morphology such as the growth of the brain or the growth of the chondrocranium relative to the dermatocranium. Such processes have been shown to have major effects on cranial morphology in mice. It is not known, however, whether they are relevant to explaining normal phenotypic variation in humans. To answer this question, we obtained vectors of shape change from mutant mouse models in which the developmental basis for the craniofacial phenotype is known to varying degrees, and compared these to a homologous dataset constructed from human crania obtained from a single population with a known genealogy. Our results show that the shape vectors associated with perturbations to chondrocranial growth, brain growth, and body size in mice do largely correspond to axes of covariation in humans. This finding supports the view that the developmental basis for craniofacial variation funnels down to a relatively small number of key developmental processes that are similar across mice and humans. Understanding these processes and how they influence craniofacial shape provides fundamental insights into the developmental basis for evolutionary change in the human skull as well as the developmental-genetic basis for normal phenotypic variation in craniofacial form.     

Craniofacial variability and modularity in macaques and mice

Evolutionary developmental biology of primates will be driven largely by the developmental biology of the house mouse. Inferences from how known developmental perturbations produce phenotypic effects in model organisms, such as mice, to how the same perturbations would affect craniofacial form in primates must be informed by comparisons of phenotypic variation and variability in mice and the primate species of interest. We use morphometric methods to compare patterns of cranial variability in homologous datasets obtained for two strains of laboratory mice and rhesus macaques. C57BL/6J represents a common genetic background for transgenic models. A/WySnJ mice exhibit altered facial morphology which results from reduction in the growth of the maxillary process during formation of the face. This is relevant to evolutionary changes in facial prognathism in nonhuman primate and human evolution. Rhesus macaques represent a nonhuman primate about which a great deal of phenotypic and genetic information is available. We find significant similarities in covariation patterns between the C57BL/6J mice and macaques. Among-trait variation in genetic and phenotypic variances are fairly concordant among the three groups, but among-trait variation in developmental stability is not. Finally, analysis of modularity based on phenotypic and genetic correlations did not reveal a consistent pattern in the three groups. We discuss the implications of these results for the study of evolutionary developmental biology of primates and outline a research strategy for integrating mouse genomics and developmental biology into this emerging field.     

Mind the Gap: Genetic Manipulation of Basicranial Growth within Synchondroses Modulates Calvarial and Facial Shape in Mice through Epigenetic Interactions

Phenotypic integration patterns in the mammalian skull have long been a focus of intense interest as a result of their suspected influence on the trajectory of hominid evolution. Here we test the hypothesis that perturbation of cartilage growth, which directly affects only the chondrocranium during development, will produce coordinated shape changes in the adult calvarium and face regardless of mechanism. Using two murine models of cartilage undergrowth that target two very different mechanisms, we show that strong reduction in cartilage growth produces a short, wide, and more flexed cranial base. This in turn produces a short, wide face in both models. Cranial base and face are already correlated early in ontogeny, and the relationship between these modules gains structure through postnatal growth and development. These results provide further evidence that there exist physical interactions between developing parts of the phenotype that produce variation at a distance from the actual locus upon which a particular selective pressure is acting. Phenotypic changes observed over the course of evolution may not all require adaptationist explanations; rather, it is likely that a substantial portion of observed phenotypic variation over the history of a clade is not directly adaptive but rather a secondary consequence of some local response to selection.     

Experimental alteration of DNA methylation affects the phenotypic plasticity of ecologically relevant traits in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

Heritable phenotypic variation in plants can be caused not only by underlying genetic differences, but also by variation in epigenetic modifications such as DNA methylation. However, we still know very little about how relevant such epigenetic variation is to the ecology and evolution of natural populations. We conducted a greenhouse experiment in which we treated a set of natural genotypes of Arabidopsis thaliana with the demethylating agent 5-azacytidine and examined the consequences of this treatment for plant traits and their phenotypic plasticity. Experimental demethylation strongly reduced the growth and fitness of plants and delayed their flowering, but the degree of this response varied significantly among genotypes. Differences in genotypes’ responses to demethylation were only weakly related to their genetic relatedness, which is consistent with the idea that natural epigenetic variation is independent of genetic variation. Demethylation also altered patterns of phenotypic plasticity, as well as the amount of phenotypic variation observed among plant individuals and genotype means. We have demonstrated that epigenetic variation can have a dramatic impact on ecologically important plant traits and their variability, as well as on the fitness of plants and their ecological interactions. Epigenetic variation may thus be an overlooked factor in the evolutionary ecology of plant populations.     

How does mutation affect the distribution of phenotypes?

The potential for mutational processes to influence patterns of neutral or adaptive phenotypic evolution is not well understood. If mutations are directionally biased, shifting trait means in a particular direction, or if mutation generates more variance in some directions of multivariate trait space than others, mutation itself might be a source of bias in phenotypic evolution. Here, we use mutagenesis to investigate the affect of mutation on trait mean and (co)variances in zebrafish, Danio rerio. Mutation altered the relationship between age and both prolonged swimming speed and body shape. These observations suggest that mutational effects on ontogeny or aging have the potential to generate variance across the phenome. Mutations had a far greater effect in males than females, although whether this is a reflection of sex‐specific ontogeny or aging remains to be determined. In males, mutations generated positive covariance between swimming speed, size, and body shape suggesting the potential for mutation to affect the evolutionary covariation of these traits. Overall, our observations suggest that mutation does not generate equal variance in all directions of phenotypic space or in each sex, and that pervasive variation in ontogeny or aging within a cohort could affect the variation available to evolution.     

Maintenance of chondroitin sulfation balance by chondroitin-4-sulfotransferase 1 is required for chondrocyte development and growth factor signaling during cartilage morphogenesis

Glycosaminoglycans (GAGs) such as heparan sulfate and chondroitin sulfate are polysaccharide chains that are attached to core proteins to form proteoglycans. The biosynthesis of GAGs is a multistep process that includes the attachment of sulfate groups to specific positions of the polysaccharide chains by sulfotransferases. Heparan-sulfate and heparan sulfate-sulfotransferases play important roles in growth factor signaling and animal development. However, the biological importance of chondroitin sulfation during mammalian development and growth factor signaling is poorly understood. We show that a gene trap mutation in the BMP-induced chondroitin-4-sulfotransferase 1 (C4st1) gene (also called carbohydrate sulfotransferase 11 - Chst11), which encodes an enzyme specific for the transfer of sulfate groups to the 4-O-position in chondroitin, causes severe chondrodysplasia characterized by a disorganized cartilage growth plate as well as specific alterations in the orientation of chondrocyte columns. This phenotype is associated with a chondroitin sulfation imbalance, mislocalization of chondroitin sulfate in the growth plate and an imbalance of apoptotic signals. Analysis of several growth factor signaling pathways that are important in cartilage growth plate development showed that the C4st1(gt/gt) mutation led to strong upregulation of TGFbeta signaling with concomitant downregulation of BMP signaling, while Indian hedgehog (Ihh) signaling was unaffected. These results show that chondroitin 4-O-sulfation by C4st1 is required for proper chondroitin sulfate localization, modulation of distinct signaling pathways and cartilage growth plate morphogenesis. Our study demonstrates an important biological role of differential chondroitin sulfation in mammalian development.     

Increased deposition of glycosaminoglycans and altered structure of heparan sulfate in idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is characterized by aberrant deposition of extracellular matrix (ECM) constituents, including glycosaminoglycans (GAGs), that may play a role in remodelling processes by influencing critical mediators such as growth factors. We hypothesize that GAGs may be altered in IPF and that this contribute to create a pro-fibrotic environment. The aim of this study was therefore to examine the fine structure of heparan sulfate (HS), chondroitin/dermatan sulfate (CS/DS) and hyaluronan (HA) in lung samples from IPF patients and from control subjects. GAGs in lung samples from severe IPF patients and donor lungs were analyzed with HPLC. HS was assessed by immunohistochemistry and collagen was quantified as hydroxyproline content. The total amount of HS, CS/DS and HA was increased in IPF lungs but there was no significant difference in the total collagen content. We found a relative increase in total sulfation of HS due to increment of 2-O, 6-O and N-sulfation and a higher proportion of sulfation in CS/DS. Highly sulfated HS was located in the border zone between denser areas and more normal looking alveolar parenchyma in basement membranes of blood vessels and airways, that were immuno-positive for perlecan, as well as on the cell surface of spindle-shaped cells in the alveolar interstitium. These findings show for the first time that both the amount and structure of glycosaminoglycans are altered in IPF. These changes may contribute to the tissue remodelling in IPF by altering growth factor retention and activity, creating a pro-fibrotic ECM landscape.     

The chloroplast triggers developmental reprogramming when mutS HOMOLOG1 is suppressed in plants

Multicellular eukaryotes demonstrate nongenetic, heritable phenotypic versatility in their adaptation to environmental changes. This inclusive inheritance is composed of interacting epigenetic, maternal, and environmental factors. Yet-unidentified maternal effects can have a pronounced influence on plant phenotypic adaptation to changing environmental conditions. To explore the control of phenotypy in higher plants, we examined the effect of a single plant nuclear gene on the expression and transmission of phenotypic variability in Arabidopsis (Arabidopsis thaliana). MutS HOMOLOG1 (MSH1) is a plant-specific nuclear gene product that functions in both mitochondria and plastids to maintain genome stability. RNA interference suppression of the gene elicits strikingly similar programmed changes in plant growth pattern in six different plant species, changes subsequently heritable independent of the RNA interference transgene. The altered phenotypes reflect multiple pathways that are known to participate in adaptation, including altered phytohormone effects for dwarfed growth and reduced internode elongation, enhanced branching, reduced stomatal density, altered leaf morphology, delayed flowering, and extended juvenility, with conversion to perennial growth pattern in short days. Some of these effects are partially reversed with the application of gibberellic acid. Genetic hemicomplementation experiments show that this phenotypic plasticity derives from changes in chloroplast state. Our results suggest that suppression of MSH1, which occurs under several forms of abiotic stress, triggers a plastidial response process that involves nongenetic inheritance.     

Analysis and visualization of growth-related and treatment-induced craniofacial changes

An analysis and visualization of craniofacial shape changes due to growth or orthodontic treatment is presented. The suggested method is based on an adapted Karhunen-Loève decomposition of time-discrete data based on landmarks in lateral X-rays of the skull. It allows for a reduction of the high-dimensional dynamic problem to a few spatial modes representing synchronous components of growth patterns with time-dependent mode coefficients. The growth-related shape changes as well as the orthodontic treatment effects are visualized by overdrawing the underlying shape changes. The results based on this technique give insight into the still controversially discussed question to which degree the craniofacial skeletal structures can be influenced by orthodontic appliances.     

Characterization of the interaction of interleukin-8 with hyaluronan, chondroitin sulfate, dermatan sulfate and their sulfated derivatives by spectroscopy and molecular modeling

The interactions between glycosaminoglycans (GAGs), important components of the extracellular matrix, and proteins such as growth factors and chemokines play critical roles in cellular regulation processes. Therefore, the design of GAG derivatives for the development of innovative materials with bio-like properties in terms of their interaction with regulatory proteins is of great interest for tissue engineering and regenerative medicine. Previous work on the chemokine interleukin-8 (IL-8) has focused on its interaction with heparin and heparan sulfate, which regulate chemokine function. However, the extracellular matrix contains other GAGs, such as hyaluronic acid (HA), dermatan sulfate (DS) and chondroitin sulfate (CS), which have so far not been characterized in terms of their distinct molecular recognition properties towards IL-8 in relation to their length and sulfation patterns. NMR and molecular modeling have been in great part the methods of choice to study the structural and recognition properties of GAGs and their protein complexes. However, separately these methods have challenges to cope with the high degree of similarity and flexibility that GAGs exhibit. In this work, we combine fluorescence spectroscopy, NMR experiments, docking and molecular dynamics simulations to study the configurational and recognition properties of IL-8 towards a series of HA and CS derivatives and DS. We analyze the effects of GAG length and sulfation patterns in binding strength and specificity, and the influence of GAG binding on IL-8 dimer formation. Our results highlight the importance of combining experimental and theoretical approaches to obtain a better understanding of the molecular recognition properties of GAG-protein systems.     

Microsatellite mutations and inferences about human demography

Microsatellites have been widely used as tools for population studies. However, inference about population processes relies on the specification of mutation parameters that are largely unknown and likely to differ across loci. Here, we use data on somatic mutations to investigate the mutation process at 14 tetranucleotide repeats and carry out an advanced multilocus analysis of different demographic scenarios on worldwide population samples. We use a method based on less restrictive assumptions about the mutation process, which is more powerful to detect departures from the null hypothesis of constant population size than other methods previously applied to similar data sets. We detect a signal of population expansion in all samples examined, except for one African sample. As part of this analysis, we identify an "anomalous" locus whose extreme pattern of variation cannot be explained by variability in mutation size. Exaggerated mutation rate is proposed as a possible cause for its unusual variation pattern. We evaluate the effect of using it to infer population histories and show that inferences about demographic histories are markedly affected by its inclusion. In fact, exclusion of the anomalous locus reduces interlocus variability of statistics summarizing population variation and strengthens the evidence in favor of demographic growth.     

Heparin--a unique stimulator of human colon cancer cells' growth

The cancer microenvironment and the interactions between cancer and surrounding tissue cells are thought to play a pivotal role in tumor development and progression. Glycosaminoglycans (GAGs)/proteoglycans (PGs) are major constituents of the extracellular matrix, the composition of which may affect various cellular functions. In the present study, the effects of GAGs on the proliferation of HT29, SW1116, and HCT116 human colon cancer cell lines were examined using exogenously added GAGs, an inhibitor of endogenous GAG sulfation and specific glycosidase digestions. Our results demonstrate that colon cancer cell growth was exclusively stimulated by exogenously added heparin and insensitive to endogenous GAGs/PGs production, in a sulfation pattern-related manner. Treatment of the tested cell lines with the FGF-2 neutralizing antibody showed that the stimulatory effect of heparin on the cells' growth was not FGF-2-dependent. Responsiveness of colon cancer cell lines to exogenous heparin/heparan sulfate may play a role in their growth and metastasis.     

Current state on the enzymatic synthesis of glycosaminoglycans

Glycosaminoglycans (GAGs) are linear anionic polysaccharides, and most of them show a specific sulfation pattern. GAGs have been studied for decades, and still, new biological functions are discovered. Hyaluronic acid and heparin are sold for medical or cosmetic applications. With increased market and applications, the production of GAGs stays in the focus of research groups and the industry. Common industrial GAG production relies on the extraction of animal tissue. Contamination, high dispersity, and uncontrolled sulfation pattern are still obstacles to this process. Tailored production strategies for the chemoenzymatic synthesis have been developed to address these obstacles. In recent years, enzyme cascades, including uridine-5′-diphosphate sugar syntheses, were established to obtain defined polymer size and dispersity, as well as defined sulfation patterns. Nevertheless, the complex synthesis of GAGs is still a challenging research field.     

How Changes in Extracellular Matrix Mechanics and Gene Expression Variability Might Combine to Drive Cancer Progression

Changes in extracellular matrix (ECM) structure or mechanics can actively drive cancer progression; however, the underlying mechanism remains unknown. Here we explore whether this process could be mediated by changes in cell shape that lead to increases in genetic noise, given that both factors have been independently shown to alter gene expression and induce cell fate switching. We do this using a computer simulation model that explores the impact of physical changes in the tissue microenvironment under conditions in which physical deformation of cells increases gene expression variability among genetically identical cells. The model reveals that cancerous tissue growth can be driven by physical changes in the microenvironment: when increases in cell shape variability due to growth-dependent increases in cell packing density enhance gene expression variation, heterogeneous autonomous growth and further structural disorganization can result, thereby driving cancer progression via positive feedback. The model parameters that led to this prediction are consistent with experimental measurements of mammary tissues that spontaneously undergo cancer progression in transgenic C3(1)-SV40Tag female mice, which exhibit enhanced stiffness of mammary ducts, as well as progressive increases in variability of cell-cell relations and associated cell shape changes. These results demonstrate the potential for physical changes in the tissue microenvironment (e.g., altered ECM mechanics) to induce a cancerous phenotype or accelerate cancer progression in a clonal population through local changes in cell geometry and increased phenotypic variability, even in the absence of gene mutation.     

Postnatal Dynamics of Developmental Stability and Canalization of Lizard Head Shape Under Different Environmental Conditions

Developmental stability (DS) and canalization are key determinants of phenotypic variation. To provide a better understanding of how postnatal growth is involved in determining the effects of DS and canalization on phenotypic variation, we studied within- and among-individual variation in head shape in ontogenetic series of lizards inhabiting urban and rural environments. Urban lizards exhibited increased fluctuating asymmetry during the early postnatal stages, but asymmetry levels decreased during growth. By contrast, asymmetry remained constant across the investigated size range in the rural population. In addition, urban juveniles were more variable for symmetric shape and deviated more from the group shape-size allometric trajectory, but both indices declined across ontogeny. Congruent patterns of within- and among-individual variation suggest that both DS and canalization may rely on similar underlying mechanisms. Further, the ontogenetic reduction of variation in the urban population suggests that compensatory growth may aid in buffering phenotypic variation and correcting deviances from the established developmental path. Alternatively, passive mechanisms and population dynamics may also explain the decrease of phenodeviants in urban populations. Significant correlations between symmetric and asymmetric shape, as well as similar integration patterns between the two populations, suggest that similar developmental mechanisms regulate head shape in both environments. Overall, these results highlight the relevance of both pre- and post-natal dynamics in determining levels of phenotypic variation, enhancing our understanding of how organisms respond to perturbations to DS and canalization under stressful conditions.     

VARIATION IN ACOUSTIC SIGNALLING TRAITS EXHIBITS FOOTPRINTS OF SEXUAL SELECTION

Phenotypic variation is ubiquitous in nature and a precondition for adaptive evolution. However, theory predicts that the extent of phenotypic variation should decrease with increasing strength of selection on a trait. Comparative analyses of trait variability have repeatedly used this expectation to infer the type or strength of selection. Yet, the suggested influence of selection on trait variability has rarely been tested empirically. In the present study, I compare estimates of sexual selection strength and trait variability from published data. I constricted the analysis to acoustic courtship traits in amphibians and insects with known variability and corresponding results of female binary choice experiments on these traits. Trait variability and strength of sexual selection were significantly correlated, and both were correlated with signal duration. Because traits under stronger selection had lower variation even after the effect of signal duration was eliminated, I conclude that traces of the strength of selection can be observed with respect to variation of acoustic signaling traits in insects and amphibians. The analysis also shows that traits under stabilizing selection have significantly lower phenotypic variability than traits under directional selection.     

CONTRASTING PATTERNS OF HERITABLE GEOGRAPHIC VARIATION IN SHELL MORPHOLOGY AND GROWTH POTENTIAL IN THE MARINE GASTROPOD BEMBICIUM VITTATUM: EVIDENCE FROM FIELD EXPERIMENTS

Similar phenotypes do not always imply similar genotypes. In species distributed over a broad latitudinal range, geographical variation in morphological and life-history traits may reflect very different relations between genotypic and environmental effects on these traits. Patterns of selection among latitudinally separated sites may minimize phenotypic differences in life-history traits but promote phenotypic differences in form. Thus, for example, latitudinal variation in temperature often leads to genetically based metabolic differences that minimize differences in growth rate among populations at different latitudes (countergradient variation). However, variation in habitat experienced by the same populations may promote genetically based differences in shell form (cogradient variation). Few attempts have been made to assess simultaneously such mosaic effects of natural selection on the genetic basis of variation in both morphological and life-history traits among geographically separated populations. I quantified the extent to which widely separated populations of the rocky shore marine gastropod Bembicium vittatum exhibited genetic differences in shell shape, shell pattern, and growth rate. Bembicium vittatum occurs naturally at only three widely separated locations on the Western Australian coast. Individuals were transplanted from all three locations to a latitudinally intermediate site, where they were released in different pairwise combinations and allowed to reproduce. F₁ offspring from crosses between same- or different-source parents were identified using allozyme markers. When grown in a common environment, offspring from same-source parents exhibited similar differences in shell shape and pattern, but dramatic differences in growth rates, compared to native populations. Genetic variation therefore exists for all three traits. Growth rates in the common environment were positively correlated with latitude of the source population, confirming the existence of countergradient variation for growth associated with metabolic compensation. In addition, for both shell shape and growth rate, hybrids exhibited phenotypes roughly midway between the same-source parents, suggesting that genetic differences have a large additive component. In contrast, when one parent had pigmented spots, the offspring also had spots, suggesting a strong dominance component to the genetic basis of shell pattern. Genetic differences therefore yield different morphological phenotypes but similar life-history phenotypes, among latitudinally distant populations, and confirm a pattern of mosaic evolution in B. vittatum.     

Effects of developmental and functional interactions on mouse cranial variability through late ontogeny

The mammalian skull performs a variety of functions and its growth and development mirrors this complexity. Cranial growth and development have been actively studied for many years. Despite this interest, the variation in the patterns and processes of skull growth has attracted little attention. An important and unanswered question is the extent to which patterns of cranial covariation and variation are dynamically reworked throughout postnatal growth. To address this question, we examine patterns of variability in random-bred mouse skulls aged 35, 90, and 150 days. Using a battery of both Procrustes coordinate and Euclidean distance-based methods, we measure mean shape, canalization, developmental stability, and morphological integration in these skulls. We predict that the patterns of variability are dynamic, particularly between the youngest and the two oldest age groups due to the influence of functional effects such as postweaning mastication. We also hypothesize that patterns of variability are structured by the same functional and developmental factors that have been shown to influence cranial growth in primates. Our results indicate that contrary to our predictions, patterns of canalization, developmental stability, and morphological integration are stabilized before 35 days. The mean shape, however, changed significantly with growth. We found that only the facial region showed significant integration as predicted by the functional matrix model used in other studies of integration. These results indicate that phenotypic integration in these mice does not closely match those found for primate species, suggesting that comparisons between species should be made with care.     

Sulfated glycosaminoglycans in protein aggregation diseases

Protein aggregation diseases are characterized by intracellular or extracellular deposition of misfolded and aggregated proteins. These aggregated deposits contain multiple proteinaceous and non-protein components that are thought to play critical roles in the etiology and pathogenesis of protein aggregation diseases in vivo. One of these components, the sulfated glycosaminoglycans (GAGs), includes heparan sulfate, chondroitin sulfate, and keratan sulfate. The sulfated GAGs are negatively charged heteropolysaccharides expressed in all mammalian tissues. Enzymatically generated structural patterns and the degree of sulfation in GAGs determine GAGs’ specific interactions with their protein ligands. Here, we review the potential roles of the sulfated GAGs in the pathogenesis and progression of protein aggregation diseases from a perspective of their sulfation modification. We also discuss the possibility of sulfated GAGs as therapeutic targets for protein aggregation diseases.     

Effect of bite force and diet composition on craniofacial diversification of Southern South American human populations

Ecological factors can be important to shape the patterns of morphological variation among human populations. Particularly, diet plays a fundamental role in craniofacial variation due to both the effect of the nutritional status—mostly dependent on the type and amount of nutrients consumed—on skeletal growth and the localized effects of masticatory forces. We examine these two dimensions of diet and evaluate their influence on morphological diversification of human populations from southern South America during the late Holocene. Cranial morphology was measured as 3D coordinates defining the face, base and vault. Size, form, and shape variables were obtained for 474 adult individuals coming from 12 samples. Diet composition was inferred from carious lesions and δ¹³C data, whereas bite forces were estimated using traits of main jaw muscles. The spatial structure of the morphological and ecological variables was measured using correlograms. The influence of diet composition and bite force on morphometric variation was estimated by a spatial regression model. Cranial variation and diet composition display a geographical structure, while no geographical pattern was observed in bite forces. Cranial variation in size and form is significantly associated with diet composition, suggesting a strong effect of systemic factors on cranial growth. Conversely, bite forces do not contribute significantly to the pattern of morphological variation among the samples analyzed. Overall, these results show that an association between diet composition and hardness cannot be assumed, and highlight the complex relationship between morphological diversification and diet in human populations. Am J Phys Anthropol 155:114–127, 2014. © 2014 Wiley Periodicals, Inc.