Finding the weakest link: mechanical sensitivity in a fish cranial linkage system

Understanding the physical mechanics behind morphological systems can offer insights into their evolution. Recent work on linkage systems in fish and crustaceans has suggested that the evolution of such systems may depend on mechanical sensitivity, where geometrical changes to different parts of a biomechanical system have variable influence on mechanical outputs. While examined at the evolutionary level, no study has directly explored this idea at the level of the mechanism. We analyse the mechanical sensitivity of a fish cranial linkage to identify the influence of linkage geometry on the kinematic transmission (KT) of the suspensorium, hyoid and lower jaw. Specifically, we answer two questions about the sensitivity of this linkage system: (i) What changes in linkage geometry affect one KT while keeping the other KTs constant? (ii) Which geometry changes result in the largest and smallest changes to KT? Our results show that there are ways to alter the morphology that change each KT individually, and that there are multiple ways to alter a single link that have variable influence on KT. These results provide insight into the morphological evolution of the fish skull and highlight which structural features in the system may have more freedom to evolve than others.

Analysis of modularity and integration suggests evolution of dragonfly wing venation mainly in response to functional demands

Insect wings show a high variability in wing venation. Selection for function, developmental pathways and phylogeny likely influenced wing vein diversification, however, quantitative data to estimate these influences and their interplay are missing. Here, it is tested how dragonfly wing vein configuration is influenced by functional demands, development, phylogeny and allometry using the concepts of modularity and integration. In an evolutionary context, modules are sets of characters that evolve in relative independence to other characters, while integration refers to a high degree of association between subparts of a structure. Results show allometric and phylogenetic signal in the wing shape variation, however, patterns of integration and modularity are not influenced by these two factors. Overall, dragonfly wings are highly integrated structures with almost no modular signal. Configuration changes in one wing vein or wing area thus influence wing shape as a whole. Moreover, the fore- and hindwings correlate with each other in their evolutionary shape variation supporting biomechanical data of wing interdependence. Despite the overall high degree of evolutionary integration, functional hypotheses of modularity could be confirmed for two wing areas, the arculus-triangle complex at the base of the wing which is responsible for passive wing folding especially during flapping flight and the location of the pterostigma, a coloured wing cell which is more heavy that other wing cells and passively regulates wing pitch as well as critical flight speeds during gliding. Although evolving as distinct modules, these specific vein regions also show high integration and evolve at the same rates like the whole wing which suggests an influence of these structures on the shape evolution of the rest of the wing. Their biomechanical role as passive regulators of wing corrugation and wing pitch suggests that these structures decisively influenced the evolution of advanced modern flight styles and explains their retention once they had evolved early within the lineage Odonatoptera.

Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae)

BACKGROUND

One of the most perplexing questions in evolutionary biology is why some lineages diversify into many species, and others do not. In many cases, ecological opportunity has played an important role, leading to diversification along trophic or habitat-based axes. The Goodeidae (Teleostomi: Cyprinodontiformes) are a family of freshwater fishes with two subfamilies: Goodeinae (42 species, viviparous, heterogeneous habitats, Mesa Central of Mexico) and Empetrichthyinae (4 species, oviparous, homogeneous habitats, Great Basin of the United States). These discrepant sets of characteristics and their sister-group relationship make the goodeids amenable to a comparative study of diversification. We gathered lateral body images from more than 1600 specimens of all extant species in the family. Geometric morphometric, and phylogenetic comparative analyses were used to address whether higher species diversity correlates with higher rates of morphological shape evolution and whether there are differences in functional/habitat modules between the two subfamilies.

RESULTS

This study recovered a higher rate of overall body shape evolution in the Goodeinae that is nearly double in magnitude compared to the Empetrichthyinae. A modularity test indicated that the Goodeinae displayed elevated rates of morphological evolution in comparison to the Empetrichthyinae when only trunk (locomotor) regions were compared between subfamilies. No significant differences in evolutionary shape rates were recovered when the trophic (head) regions were compared between subfamilies.

DISCUSSION

These results support the hypothesis that Mexican goodeids radiated via an ecological opportunity scenario into a wide-array of novel habitats in the island-like Mesa Central as evidenced by their high rate of shape evolution, relative to the Empetrichthyinae. This study quantitatively unraveled the drivers of evolution and eliminated trophic specialization as a driving force within the Goodeidae.

CONCLUSIONS

A combination of phylogenetic and morphometric data, and phylogenetic comparative analyses were used to examine body shape rate evolution within the Goodeidae. Results support the hypothesis that species in the subfamily Goodeinae on the central Mexican plateau had a higher rate of body shape evolution relative to its sister subfamily Empetrichthyinae in the Great Basin suggesting that the Goodeinae diversified via an ecological opportunity scenario along habitat, rather than trophic axes.

Modularity promotes morphological divergence in ray-finned fishes

Modularity is considered a prerequisite for the evolvability of biological systems. This is because in theory, individual modules can follow quasi-independent evolutionary trajectories or evolve at different rates compared to other aspects of the organism. This may influence the potential of some modules to diverge, leading to differences in disparity. Here, we investigated this relationship between modularity, rates of morphological evolution and disparity using a phylogenetically diverse sample of ray-finned fishes. We compared the support for multiple hypotheses of evolutionary modularity and asked if the partitions delimited by the best-fitting models were also characterized by the highest evolutionary rate differentials. We found that an evolutionary module incorporating the dorsal, anal and paired fins was well supported by the data, and that this module evolves more rapidly and consequently generates more disparity than other modules. This suggests that modularity may indeed promote morphological disparity through differences in evolutionary rates across modules.

Patterns of growth in the presacral vertebral column of the leopard gecko (Eublepharis macularius)

Postnatal growth patterns within the vertebral column may be informative about body proportions and regionalization. We measured femur length, lengths of all pre-sacral vertebrae, and lengths of intervertebral spaces, from radiographs of a series of 21 Eublepharis macularius, raised under standard conditions and covering most of the ontogenetic body size range. Vertebrae were grouped into cervical, sternal, and dorsal compartments, and lengths of adjacent pairs of vertebrae were summed before analysis. Femur length was included as an index of body size. Principal component analysis of the variance-covariance matrix of these data was used to investigate scaling among them. PC1 explained 94.19% of total variance, interpreted as the variance due to body size. PC1 differed significantly from the hypothetical isometric vector, indicating overall allometry. The atlas and axis vertebrae displayed strong negative allometry; the remainder of the vertebral pairs exhibited weak negative allometry, isometry or positive allometry. PC1 explained a markedly smaller amount of variance for the vertebral pairs of the cervical compartment than for the remainder of the vertebral pairs, with the exception of the final pair. The relative standard deviations of the eigenvalues from the PCAs of the three vertebral compartments indicated that the vertebrae of the cervical compartment were less strongly integrated by scaling than were the sternal or dorsal vertebrae, which did not differ greatly between themselves in their strong integration, suggesting that the growth of the cervical vertebrae is constrained by the mechanical requirements of the head. Regionalization of the remainder of the vertebral column is less clearly defined but may be associated with wave form propagation incident upon locomotion, and by locomotory changes occasioned by tail autotomy and regeneration. Femur length exhibits negative allometry relative to individual vertebral pairs and to vertebral column length, suggesting a change in locomotor requirements over the ontogenetic size range.

ShapeRotator: An R tool for standardized rigid rotations of articulated three-dimensional structures with application for geometric morphometrics

The quantification of complex morphological patterns typically involves comprehensive shape and size analyses, usually obtained by gathering morphological data from all the structures that capture the phenotypic diversity of an organism or object. Articulated structures are a critical component of overall phenotypic diversity, but data gathered from these structures are difficult to incorporate into modern analyses because of the complexities associated with jointly quantifying 3D shape in multiple structures. While there are existing methods for analyzing shape variation in articulated structures in two-dimensional (2D) space, these methods do not work in 3D, a rapidly growing area of capability and research. Here, we describe a simple geometric rigid rotation approach that removes the effect of random translation and rotation, enabling the morphological analysis of 3D articulated structures. Our method is based on Cartesian coordinates in 3D space, so it can be applied to any morphometric problem that also uses 3D coordinates (e.g., spherical harmonics). We demonstrate the method by applying it to a landmark-based dataset for analyzing shape variation using geometric morphometrics. We have developed an R tool (ShapeRotator) so that the method can be easily implemented in the commonly used R package geomorph and MorphoJ software. This method will be a valuable tool for 3D morphological analyses in articulated structures by allowing an exhaustive examination of shape and size diversity.

Study of morphological variation of northern Neotropical Ariidae reveals conservatism despite macrohabitat transitions

BACKGROUND

Morphological convergence triggered by trophic adaptations is a common pattern in adaptive radiations. The study of shape variation in an evolutionary context is usually restricted to well-studied fish models. We take advantage of the recently revised systematics of New World Ariidae and investigate skull shape evolution in six genera of northern Neotropical Ariidae. They constitute a lineage that diversified in the marine habitat but repeatedly adapted to freshwater habitats. 3D geometric morphometrics was applied for the first time in catfish skulls and phylogenetically informed statistical analyses were performed to test for the impact of habitat on skull diversification after habitat transition in this lineage.

RESULTS

We found that skull shape is conserved throughout phylogeny. A morphospace analysis revealed that freshwater and marine species occupy extreme ends of the first principal component axis and that they exhibit similar Procrustes variances. Yet freshwater species occupy the smallest shape space compared to marine and brackish species (based on partial disparity), and marine and freshwater species have the largest Procrustes distance to each other. We observed a single case of shape convergence as derived from 'C-metrics', which cannot be explained by the occupation of the same habitat.

CONCLUSIONS

Although Ariidae occupy such a broad spectrum of different habitats from sea to freshwater, the morphospace analysis and analyses of shape and co-variation with habitat in a phylogenetic context shows that conservatism dominates skull shape evolution among ariid genera.

Integration of the nasal complex: Implications for developmental and evolutionary change in modern humans

OBJECTIVES

Assessing the strength of integration among different regions of the modern human nasal complex is important for developing a more thorough understanding of the determinants of nasal morphology. Given the morphogenetic influence of cartilage on adjacent intramembranous growth sites, the interaction between chondrocranial- versus intramembranous-derived nasal structures may have a significant influence on patterns of nasal variation. The purpose of this study is to examine integration between the chondrocranial- and intramembranous-derived regions of the nasal complex.

MATERIALS AND METHODS

Using computed tomograph (CT) scans, we collected three-dimensional coordinate landmark data from a static adult sample (n = 62). First, using centroid size, and the symmetric and asymmetric components of shape variation, we examined the strength of integration between landmarks representing chondrocranial-derived structures (e.g., ethmoid, external nasal cartilages) and landmarks representing intramembranous-derived structures (nasal floor, anterior nasal aperture, etc.). Second, given that the strength of integration is a relative measure, we compared integration between chondrocranial- and intramembranous-derived structures to the more modularized external and internal regions of the nasal complex.

RESULTS

There was significant moderate morphological integration between chondrocranial- versus intramembranous-derived regions of the nasal complex. Moreover, integration between chondrocranial- versus intramembranous-derived structures was consistently stronger when compared to external versus internal regions for both the symmetric and asymmetric components of variation. Thus, more covariation within the nasal complex could be explained by the relationship between chondrocranial- and intramembranous-derived structures.

CONCLUSIONS

Our results suggest that the interaction between chondrocranial- and intramembranous-derived structures may be an important determinant in the patterning of nasal complex variation.

Evolution of sex-specific pace-of-life syndromes: genetic architecture and physiological mechanisms

Sex differences in life history, physiology, and behavior are nearly ubiquitous across taxa, owing to sex-specific selection that arises from different reproductive strategies of the sexes. The pace-of-life syndrome (POLS) hypothesis predicts that most variation in such traits among individuals, populations, and species falls along a slow-fast pace-of-life continuum. As a result of their different reproductive roles and environment, the sexes also commonly differ in pace-of-life, with important consequences for the evolution of POLS. Here, we outline mechanisms for how males and females can evolve differences in POLS traits and in how such traits can covary differently despite constraints resulting from a shared genome. We review the current knowledge of the genetic basis of POLS traits and suggest candidate genes and pathways for future studies. Pleiotropic effects may govern many of the genetic correlations, but little is still known about the mechanisms involved in trade-offs between current and future reproduction and their integration with behavioral variation. We highlight the importance of metabolic and hormonal pathways in mediating sex differences in POLS traits; however, there is still a shortage of studies that test for sex specificity in molecular effects and their evolutionary causes. Considering whether and how sexual dimorphism evolves in POLS traits provides a more holistic framework to understand how behavioral variation is integrated with life histories and physiology, and we call for studies that focus on examining the sex-specific genetic architecture of this integration.

Contrasting lengths of Pelargonium floral nectar tubes result from late differences in rate and duration of growth

BACKGROUND AND AIMS

Much of morphological evolution in flowers has arisen from pollinator-mediated selection, often manifest as a match between the length of the pollinator's proboscis and the depth of tubular corollas or spurs. We investigate development, growth and homology of the unique nectar tube of Pelargonium, frequently described as 'a spur adnate to the pedicel'.

METHODS

We focused on two species. The nectar tube of P. ionidiflorum is three times longer than that of P. odoratissimum. Light and scanning electron microscopy were carried out, and daily growth measurements were used to compare nectar tube development and vascular patterns.

KEY RESULTS

Nectar tubes in both species are initiated centripetally to the dorsal sepal in a space created by lateral displacement of two antepetalous stamens. The cavity deepens through subsequent intercalary growth of the receptacle that proceeds at the same rate in both species until tubes reach approx. 10 mm in length. Differences in final nectar tube lengths arise via an increase in the rate and duration of growth of the receptacle that begins just before anthesis (floral opening) and continues for several days past anthesis in P. ionidiflorum but does not occur in P. odoratissimum. Epidermal cells of the dorsal surface of the nectar tube in P. ionidiflorum are approx. 1.6 times longer than those in P. odoratissimum. Histological sections show no evidence that the nectar tube is a spur that became evolutionarily fused to the pedicel.

CONCLUSIONS

Nectar tubes in Pelargonium are localized cavities that form in the receptacle via intercalary growth. Differences in the rate and duration of growth just prior to and following anthesis underlie differences in final tube lengths. Because differences in cell lengths do not fully account for differences in nectar tube lengths, evolutionary diversification must involve changes in both cell cycle and cell expansion.

Genome-wide mapping of global-to-local genetic effects on human facial shape

Genome-wide association scans of complex multipartite traits like the human face typically use preselected phenotypic measures. Here we report a data-driven approach to phenotyping facial shape at multiple levels of organization, allowing for an open-ended description of facial variation while preserving statistical power. In a sample of 2,329 persons of European ancestry, we identified 38 loci, 15 of which replicated in an independent European sample (n = 1,719). Four loci were completely new. For the others, additional support (n = 9) or pleiotropic effects (n = 2) were found in the literature, but the results reported here were further refined. All 15 replicated loci highlighted distinctive patterns of global-to-local genetic effects on facial shape and showed enrichment for active chromatin elements in human cranial neural crest cells, suggesting an early developmental origin of the facial variation captured. These results have implications for studies of facial genetics and other complex morphological traits.

Phenotypic integration and the evolution of signal repertoires: A case study of treefrog acoustic communication

Animal signals are inherently complex phenotypes with many interacting parts combining to elicit responses from receivers. The pattern of interrelationships between signal components reflects the extent to which each component is expressed, and responds to selection, either in concert with or independently of others. Furthermore, many species have complex repertoires consisting of multiple signal types used in different contexts, and common morphological and physiological constraints may result in interrelationships extending across the multiple signals in species' repertoires. The evolutionary significance of interrelationships between signal traits can be explored within the framework of phenotypic integration, which offers a suite of quantitative techniques to characterize complex phenotypes. In particular, these techniques allow for the assessment of modularity and integration, which describe, respectively, the extent to which sets of traits covary either independently or jointly. Although signal and repertoire complexity are thought to be major drivers of diversification and social evolution, few studies have explicitly measured the phenotypic integration of signals to investigate the evolution of diverse communication systems. We applied methods from phenotypic integration studies to quantify integration in the two primary vocalization types (advertisement and aggressive calls) in the treefrogs Hyla versicolor, Hyla cinerea, and Dendropsophus ebraccatus. We recorded male calls and calculated standardized phenotypic variance-covariance (P) matrices for characteristics within and across call types. We found significant integration across call types, but the strength of integration varied by species and corresponded with the acoustic similarity of the call types within each species. H. versicolor had the most modular advertisement and aggressive calls and the least acoustically similar call types. Additionally, P was robust to changing social competition levels in H. versicolor. Our findings suggest new directions in animal communication research in which the complex relationships among the traits of multiple signals are a key consideration for understanding signal evolution.

Developmental constraint through negative pleiotropy in the zygomatic arch

BACKGROUND

Previous analysis suggested that the relative contribution of individual bones to regional skull lengths differ between inbred mouse strains. If the negative correlation of adjacent bone lengths is associated with genetic variation in a heterogeneous population, it would be an example of negative pleiotropy, which occurs when a genetic factor leads to opposite effects in two phenotypes. Confirming negative pleiotropy and determining its basis may reveal important information about the maintenance of overall skull integration and developmental constraint on skull morphology.

RESULTS

We identified negative correlations between the lengths of the frontal and parietal bones in the midline cranial vault as well as the zygomatic bone and zygomatic process of the maxilla, which contribute to the zygomatic arch. Through gene association mapping of a large heterogeneous population of Diversity Outbred (DO) mice, we identified a quantitative trait locus on chromosome 17 driving the antagonistic contribution of these two zygomatic arch bones to total zygomatic arch length. Candidate genes in this region were identified and real-time PCR of the maxillary processes of DO founder strain embryos indicated differences in the RNA expression levels for two of the candidate genes, Camkmt and Six2.

CONCLUSIONS

A genomic region underlying negative pleiotropy of two zygomatic arch bones was identified, which provides a mechanism for antagonism in component bone lengths while constraining overall zygomatic arch length. This type of mechanism may have led to variation in the contribution of individual bones to the zygomatic arch noted across mammals. Given that similar genetic and developmental mechanisms may underlie negative correlations in other parts of the skull, these results provide an important step toward understanding the developmental basis of evolutionary variation and constraint in skull morphology.

Mechanical sensitivity and the dynamics of evolutionary rate shifts in biomechanical systems

The influence of biophysical relationships on rates of morphological evolution is a cornerstone of evolutionary theory. Mechanical sensitivity-the correlation strength between mechanical output and the system's underlying morphological components-is thought to impact the evolutionary dynamics of form-function relationships, yet has rarely been examined. Here, we compare the evolutionary rates of the mechanical components of the four-bar linkage system in the raptorial appendage of mantis shrimp (Order Stomatopoda). This system's mechanical output (kinematic transmission (KT)) is highly sensitive to variation in its output link, and less sensitive to its input and coupler links. We found that differential mechanical sensitivity is associated with variation in evolutionary rate: KT and the output link exhibit faster rates of evolution than the input and coupler links to which KT is less sensitive. Furthermore, for KT and, to a lesser extent, the output link, rates of evolution were faster in 'spearing' stomatopods than 'smashers', indicating that mechanical sensitivity may influence trait-dependent diversification. Our results suggest that mechanical sensitivity can impact morphological evolution and guide the process of phenotypic diversification. The connection between mechanical sensitivity and evolutionary rates provides a window into the interaction between physical rules and the evolutionary dynamics of morphological diversification.

Developmental pathways inferred from modularity, morphological integration and fluctuating asymmetry patterns in the human face

Facial asymmetries are usually measured and interpreted as proxies to developmental noise. However, analyses focused on its developmental and genetic architecture are scarce. To advance on this topic, studies based on a comprehensive and simultaneous analysis of modularity, morphological integration and facial asymmetries including both phenotypic and genomic information are needed. Here we explore several modularity hypotheses on a sample of Latin American mestizos, in order to test if modularity and integration patterns differ across several genomic ancestry backgrounds. To do so, 4104 individuals were analyzed using 3D photogrammetry reconstructions and a set of 34 facial landmarks placed on each individual. We found a pattern of modularity and integration that is conserved across sub-samples differing in their genomic ancestry background. Specifically, a signal of modularity based on functional demands and organization of the face is regularly observed across the whole sample. Our results shed more light on previous evidence obtained from Genome Wide Association Studies performed on the same samples, indicating the action of different genomic regions contributing to the expression of the nose and mouth facial phenotypes. Our results also indicate that large samples including phenotypic and genomic metadata enable a better understanding of the developmental and genetic architecture of craniofacial phenotypes.

Musculoskeletal networks reveal topological disparity in mammalian neck evolution

BACKGROUND

The increase in locomotor and metabolic performance during mammalian evolution was accompanied by the limitation of the number of cervical vertebrae to only seven. In turn, nuchal muscles underwent a reorganization while forelimb muscles expanded into the neck region. As variation in the cervical spine is low, the variation in the arrangement of the neck muscles and their attachment sites (i.e., the variability of the neck's musculoskeletal organization) is thus proposed to be an important source of neck disparity across mammals. Anatomical network analysis provides a novel framework to study the organization of the anatomical arrangement, or connectivity pattern, of the bones and muscles that constitute the mammalian neck in an evolutionary context.

RESULTS

Neck organization in mammals is characterized by a combination of conserved and highly variable network properties. We uncovered a conserved regionalization of the musculoskeletal organization of the neck into upper, mid and lower cervical modules. In contrast, there is a varying degree of complexity or specialization and of the integration of the pectoral elements. The musculoskeletal organization of the monotreme neck is distinctively different from that of therian mammals.

CONCLUSIONS

Our findings reveal that the limited number of vertebrae in the mammalian neck does not result in a low musculoskeletal disparity when examined in an evolutionary context. However, this disparity evolved late in mammalian history in parallel with the radiation of certain lineages (e.g., cetartiodactyls, xenarthrans). Disparity is further facilitated by the enhanced incorporation of forelimb muscles into the neck and their variability in attachment sites.

Evolution of a complex phenotype with biphasic ontogeny: Contribution of development versus function and climatic variation to skull modularity in toads

The theory of morphological integration and modularity predicts that if functional correlations among traits are relevant to mean population fitness, the genetic basis of development will be molded by stabilizing selection to match functional patterns. Yet, how much functional interactions actually shape the fitness landscape is still an open question. We used the anuran skull as a model of a complex phenotype for which we can separate developmental and functional modularity. We hypothesized that functional modularity associated to functional demands of the adult skull would overcome developmental modularity associated to bone origin at the larval phase because metamorphosis would erase the developmental signal. We tested this hypothesis in toad species of the Rhinella granulosa complex using species phenotypic correlation pattern (P-matrices). Given that the toad species are distributed in very distinct habitats and the skull has important functions related to climatic conditions, we also hypothesized that differences in skull trait covariance pattern are associated to differences in climatic variables among species. Functional and hormonal-regulated modules are more conspicuous than developmental modules only when size variation is retained on species P-matrices. Without size variation, there is a clear modularity signal of developmental units, but most species have the functional model as the best supported by empirical data without allometric size variation. Closely related toad species have more similar climatic niches and P-matrices than distantly related species, suggesting phylogenetic niche conservatism. We infer that the modularity signal due to embryonic origin of bones, which happens early in ontogeny, is blurred by the process of growth that occurs later in ontogeny. We suggest that the species differing in the preferred modularity model have different demands on the orbital functional unit and that species contrasting in climate are subjected to divergent patterns of natural selection associated to neurocranial allometry and T3 hormone regulation.

The genetics and development of mandibles and hypopharyngeal sclerite and cornua in larvae of Drosophila gaucha

The genetics and epigenetic processes associated with morphological organization are a principal aim of biology, ranging from cohesion between cells to shape and size of organisms. We investigate the post-embryonic development of Hypopharyngeal sclerite and cornua HPC and mandibles M of Drosophila gaucha larva. Integrated functioning of these Cephalopharyngeal skeleton parts of D. gaucha larva is essential for food acquisition, participating in locomotion and microhabitat selection. We examined two isolates by recording the growth of the HPC and M every 24 h for 8 days in parental, F₁, F₂ and backcross larvae. In F₁ larvae, the HPC and M growth was similar to the parental. In F₂ and backcross larvae, the growth was slower. Epistasis and dominance are the principal sources upon which the growth of HPC and M are based. Pleiotropic genes seem also to be involved in integrating the development of M and HPC. Our data suggest that hybridization of the isolates modified epigenetic processes involved in the development of those morphological structures of D. gaucha larva.

Sexual dimorphism, phenotypic integration, and the evolution of head structure in casque-headed lizards

Sexes can differ in features associated with differential reproduction, which can be used during courtship or aggressive encounters. Some traits tend to evolve independently between sexes and emerge as sexually dimorphic within the organismal phenotype. We characterize such a relationship by estimating the phenotypic integration of the head morphology and modularity of the crest in the casque‐headed lizards (Corytophanidae). In this clade, some species show extreme sexual dimorphism (e.g., head crests in the genus Basiliscus) while in others, both sexes are monomorphic. To characterize these patterns, we define phenotypic integration at the interspecific level as a pattern or network of traits evidenced by phylogenetically adjusted correlations that persist among species. At this level, modularity is an increased connectedness (e.g., higher correlation) among sections of these networks that persist in a lineage during the evolution of complex phenotypes. To test both concepts, we used phylogenetic geomorphometrics to characterize the head structure of corytophanid lizards, based on a time‐calibrated phylogeny that includes candidate fossil ancestors. We found evidence of an older diversification of corytophanids than previously reported (~67 vs. ~23.5 MYA) and show that this clade includes two morphological head architectures: (1) Sexually dimorphic crests present in males that are evolving independently from the rest of the head structure, and (2) full integration of the head morphology in monomorphic species. We propose that both architectures are optimal evolutionary trajectories of the parietal crest bones in the head of these lizards. In sexually dimorphic species, these bones are elongated and thinner, and gave rise to the extended crest used in male courtship displays. In monomorphic species, the parietal crest grew thicker in both sexes to allow for a better insertion of muscles associated with a stronger bite.

Functional modularity in lake-dwelling characin fishes of Mexico

Modular evolution promotes evolutionary change, allowing independent variation across morphological units. Recent studies have shown that under contrasting ecological pressures, patterns of modularity could be related to divergent evolution. The main goal of the present study was to evaluate the presence of modular evolution in two sister lacustrine species, Astyanax aeneus and A. caballeroi, which are differentiated by their trophic habits. Two different datasets were analyzed: (1) skull X-rays from 73 specimens (35 A. aeneus and 38 A. caballeroi) to characterize skull variation patterns, considering both species and sex effects. For this dataset, three different modularity hypotheses were tested, previously supported in other lacustrine divergent species; (2) a complete body shape dataset was also tested for four modularity hypotheses, which included a total of 196 individuals (110 Astyanax aeneus and 86 A. caballeroi). Skull shape showed significant differences among species and sex (P < 0.001), where Astyanax caballeroi species showed an upwardly projected mandible and larger preorbital region. For the skull dataset, the modularity hypothesis ranked first included three partitioning modules. While for the complete body dataset the best ranked hypothesis included two modules (head vs the rest of the body), being significant only for A. caballeroi.

The skull roof tracks the brain during the evolution and development of reptiles including birds

Major transformations in brain size and proportions, such as the enlargement of the brain during the evolution of birds, are accompanied by profound modifications to the skull roof. However, the hypothesis of concerted evolution of shape between brain and skull roof over major phylogenetic transitions, and in particular of an ontogenetic relationship between specific regions of the brain and the skull roof, has never been formally tested. We performed 3D morphometric analyses to examine the deep history of brain and skull-roof morphology in Reptilia, focusing on changes during the well-documented transition from early reptiles through archosauromorphs, including nonavian dinosaurs, to birds. Non-avialan taxa cluster tightly together in morphospace, whereas Archaeopteryx and crown birds occupy a separate region. There is a one-to-one correspondence between the forebrain and frontal bone and the midbrain and parietal bone. Furthermore, the position of the forebrain-midbrain boundary correlates significantly with the position of the frontoparietal suture across the phylogenetic breadth of Reptilia and during the ontogeny of individual taxa. Conservation of position and identity in the skull roof is apparent, and there is no support for previous hypotheses that the avian parietal is a transformed postparietal. The correlation and apparent developmental link between regions of the brain and bony skull elements are likely to be ancestral to Tetrapoda and may be fundamental to all of Osteichthyes, coeval with the origin of the dermatocranium.

Comparison of musculoskeletal networks of the primate forelimb

Anatomical network analysis is a framework for quantitatively characterizing the topological organization of anatomical structures, thus providing a way to compare structural integration and modularity among species. Here we apply this approach to study the macroevolution of the forelimb in primates, a structure whose proportions and functions vary widely within this group. We analyzed musculoskeletal network models in 22 genera, including members of all major extant primate groups and three outgroup taxa, after an extensive literature survey and dissections. The modules of the proximal limb are largely similar among taxa, but those of the distal limb show substantial variation. Some network parameters are similar within phylogenetic groups (e.g., non-primates, strepsirrhines, New World monkeys, and hominoids). Reorganization of the modules in the hominoid hand compared to other primates may relate to functional changes such as coordination of individual digit movements, increased pronation/supination, and knuckle-walking. Surprisingly, humans are one of the few taxa we studied in which the thumb musculoskeletal structures do not form an independent anatomical module. This difference may be caused by the loss in humans of some intrinsic muscles associated with the digits or the acquisition of additional muscles that integrate the thumb more closely with surrounding structures.

A population level atlas of Mus musculus craniofacial skeleton and automated image-based shape analysis

Laboratory mice are staples for evo/devo and genetics studies. Inbred strains provide a uniform genetic background to manipulate and understand gene-environment interactions, while their crosses have been instrumental in studies of genetic architecture, integration and modularity, and mapping of complex biological traits. Recently, there have been multiple large-scale studies of laboratory mice to further our understanding of the developmental basis, evolution, and genetic control of shape variation in the craniofacial skeleton (i.e. skull and mandible). These experiments typically use micro-computed tomography (micro-CT) to capture the craniofacial phenotype in 3D and rely on manually annotated anatomical landmarks to conduct statistical shape analysis. Although the common choice for imaging modality and phenotyping provides the potential for collaborative research for even larger studies with more statistical power, the investigator (or lab-specific) nature of the data collection hampers these efforts. Investigators are rightly concerned that subtle differences in how anatomical landmarks were recorded will create systematic bias between studies that will eventually influence scientific findings. Even if researchers are willing to repeat landmark annotation on a combined dataset, different lab practices and software choices may create obstacles for standardization beyond the underlying imaging data. Here, we propose a freely available analysis system that could assist in the standardization of micro-CT studies in the mouse. Our proposal uses best practices developed in biomedical imaging and takes advantage of existing open-source software and imaging formats. Our first contribution is the creation of a synthetic template for the adult mouse craniofacial skeleton from 25 inbred strains and five F1 crosses that are widely used in biological research. The template contains a fully segmented cranium, left and right hemi-mandibles, endocranial space, and the first few cervical vertebrae. We have been using this template in our lab to segment and isolate cranial structures in an automated fashion from a mixed population of mice, including craniofacial mutants, aged 4-12.5 weeks. As a secondary contribution, we demonstrate an application of nearly automated shape analysis, using symmetric diffeomorphic image registration. This approach, which we call diGPA, closely approximates the popular generalized Procrustes analysis (GPA) but negates the collection of anatomical landmarks. We achieve our goals by using the open-source advanced normalization tools (ANT) image quantification library, as well as its associated R library (ANTsR) for statistical image analysis. Finally, we make a plea to investigators to commit to using open imaging standards and software in their labs to the extent possible to increase the potential for data exchange and improve the reproducibility of findings. Future work will incorporate more anatomical detail (such as individual cranial bones, turbinals, dentition, middle ear ossicles) and more diversity into the template.

Genome-wide association studies dissect the genetic networks underlying agronomical traits in soybean

BACKGROUND

Soybean (Glycine max [L.] Merr.) is one of the most important oil and protein crops. Ever-increasing soybean consumption necessitates the improvement of varieties for more efficient production. However, both correlations among different traits and genetic interactions among genes that affect a single trait pose a challenge to soybean breeding.

RESULTS

To understand the genetic networks underlying phenotypic correlations, we collected 809 soybean accessions worldwide and phenotyped them for two years at three locations for 84 agronomic traits. Genome-wide association studies identified 245 significant genetic loci, among which 95 genetically interacted with other loci. We determined that 14 oil synthesis-related genes are responsible for fatty acid accumulation in soybean and function in line with an additive model. Network analyses demonstrated that 51 traits could be linked through the linkage disequilibrium of 115 associated loci and these links reflect phenotypic correlations. We revealed that 23 loci, including the known Dt1, E2, E1, Ln, Dt2, Fan, and Fap loci, as well as 16 undefined associated loci, have pleiotropic effects on different traits.

CONCLUSIONS

This study provides insights into the genetic correlation among complex traits and will facilitate future soybean functional studies and breeding through molecular design.

Comparative Analyses of Phenotypic Trait Covariation within and among Populations

Many morphological, behavioral, physiological, and life-history traits covary across the biological scales of individuals, populations, and species. However, the processes that cause traits to covary also change over these scales, challenging our ability to use patterns of trait covariance to infer process. Trait relationships are also widely assumed to have generic functional relationships with similar evolutionary potentials, and even though many different trait relationships are now identified, there is little appreciation that these may influence trait covariation and evolution in unique ways. We use a trait-performance-fitness framework to classify and organize trait relationships into three general classes, address which ones more likely generate trait covariation among individuals in a population, and review how selection shapes phenotypic covariation. We generate predictions about how trait covariance changes within and among populations as a result of trait relationships and in response to selection and consider how these can be tested with comparative data. Careful comparisons of covariation patterns can narrow the set of hypothesized processes that cause trait covariation when the form of the trait relationship and how it responds to selection yield clear predictions about patterns of trait covariation. We discuss the opportunities and limitations of comparative approaches to evaluate hypotheses about the evolutionary causes and consequences of trait covariation and highlight the importance of evaluating patterns within populations replicated in the same and in different selective environments. Explicit hypotheses about trait relationships are key to generating effective predictions about phenotype and its evolution using covariance data.

Integration of Brain and Skull in Prenatal Mouse Models of Apert and Crouzon Syndromes

The brain and skull represent a complex arrangement of integrated anatomical structures composed of various cell and tissue types that maintain structural and functional association throughout development. Morphological integration, a concept developed in vertebrate morphology and evolutionary biology, describes the coordinated variation of functionally and developmentally related traits of organisms. Syndromic craniosynostosis is characterized by distinctive changes in skull morphology and perceptible, though less well studied, changes in brain structure and morphology. Using mouse models for craniosynostosis conditions, our group has precisely defined how unique craniosynostosis causing mutations in fibroblast growth factor receptors affect brain and skull morphology and dysgenesis involving coordinated tissue-specific effects of these mutations. Here we examine integration of brain and skull in two mouse models for craniosynostosis: one carrying the FGFR2c C342Y mutation associated with Pfeiffer and Crouzon syndromes and a mouse model carrying the FGFR2 S252W mutation, one of two mutations responsible for two-thirds of Apert syndrome cases. Using linear distances estimated from three-dimensional coordinates of landmarks acquired from dual modality imaging of skull (high resolution micro-computed tomography and magnetic resonance microscopy) of mice at embryonic day 17.5, we confirm variation in brain and skull morphology in Fgfr2c^C342Y/+ mice, Fgfr2^+/S252W mice, and their unaffected littermates. Mutation-specific variation in neural and cranial tissue notwithstanding, patterns of integration of brain and skull differed only subtly between mice carrying either the FGFR2c C342Y or the FGFR2 S252W mutation and their unaffected littermates. However, statistically significant and substantial differences in morphological integration of brain and skull were revealed between the two mutant mouse models, each maintained on a different strain. Relative to the effects of disease-associated mutations, our results reveal a stronger influence of the background genome on patterns of brain-skull integration and suggest robust genetic, developmental, and evolutionary relationships between neural and skeletal tissues of the head.

Developmental and Evolutionary Significance of the Zygomatic Bone

The zygomatic bone is derived evolutionarily from the orbital series. In most modern mammals the zygomatic bone forms a large part of the face and usually serves as a bridge that connects the facial skeleton to the neurocranium. Our aim is to provide information on the contribution of the zygomatic bone to variation in midfacial protrusion using three samples; humans, domesticated dogs, and monkeys. In each case, variation in midface protrusion is a heritable trait produced by one of three classes of transmission: localized dysmorphology associated with single gene dysfunction, selective breeding, or long‐term evolution from a common ancestor. We hypothesize that the shape of the zygomatic bone reflects its role in stabilizing the connection between facial skeleton and neurocranium and consequently, changes in facial protrusion are more strongly reflected by the maxilla and premaxilla. Our geometric morphometric analyses support our hypothesis suggesting that the shape of the zygomatic bone has less to do with facial protrusion. By morphometrically dissecting the zygomatic bone we have determined a degree of modularity among parts of the midfacial skeleton suggesting that these components have the ability to vary independently and thus can evolve differentially. From these purely morphometric data, we propose that the neural crest cells that are fated to contribute to the zygomatic bone experience developmental cues that distinguish them from the maxilla and premaxilla. The spatiotemporal and molecular identity of the cues that impart zygoma progenitors with their identity remains an open question that will require alternative data sets. Anat Rec, 299:1616–1630, 2016. © 2016 The Authors The Anatomical Record Published by Wiley Periodicals, Inc.

Methods for Studying Morphological Integration and Modularity

Morphological integration and modularity are closely related concepts about how different traits of an organism are correlated. Integration is the overall pattern of intercorrelation; modularity is the partitioning of integration into evolutionarily or developmentally independent blocks of traits. Modularity and integration are usually studied using quantitative phenotypic data, which can be obtained either from extant or fossil organisms. Many methods are now available to study integration and modularity, all of which involve the analysis of patterns found in trait correlation or covariance matrices. We review matrix correlation, random skewers, fluctuating asymmetry, cluster analysis, Euclidean distance matrix analysis (EDMA), graphical modelling, two-block partial least squares, RV coefficients, and theoretical matrix modelling and discuss their similarities and differences. We also review different coefficients that are used to measure correlations. We apply all the methods to cranial landmark data from and ontogenetic series of Japanese macaques,Macaca fuscatato illustrate the methods and their individual strengths and weaknesses. We conclude that the exploratory approaches (cluster analyses of various sorts) were less informative and less consistent with one another than were the results of model testing or comparative approaches. Nevertheless, we found that competing models of modularity and integration are often similar enough that they are not statistically distinguishable; we expect, therefore, that several models will often be significantly correlated with observed data.

Geometric morphometrics analysis of the hind wing of leaf beetles: proximal and distal parts are separate modules

The success of beetles is mainly attributed to the possibility to hide the hindwings under the sclerotised elytra. The acquisition of the transverse folding function of the hind wing is an important event in the evolutionary history of beetles. In this study, the morphological and functional variances in the hind wings of 94 leaf beetle species (Coleoptera: Chrysomelinae) is explored using geometric morphometrics based on 36 landmarks. Principal component analysis and Canonical variate analysis indicate that changes of apical area, anal area, and middle area are three useful phylogenetic features at a subtribe level of leaf beetles. Variances of the apical area are the most obvious, which strongly influence the entire venation variance. Partial least squares analysis indicates that the proximal and distal parts of hind wings are weakly associated. Modularity tests confirm that the proximal and distal compartments of hind wings are separate modules. It is deduced that for leaf beetles, or even other beetles, the hind wing possibly exhibits significant functional divergences that occurred during the evolution of transverse folding that resulted in the proximal and distal compartments of hind wings evolving into separate functional modules.

Non-decoupled morphological evolution of the fore- and hindlimb of sabretooth predators

Specialized organisms are useful for exploring the combined effects of selection of functional traits and developmental constraints on patterns of phenotypic integration. Sabretooth predators are one of the most interesting examples of specialization among mammals. Their hypertrophied, sabre-shaped upper canines and their powerfully built forelimbs have been interpreted as adaptations to a highly specialized predatory behaviour. Given that the elongated and laterally compressed canines of sabretooths were more vulnerable to fracture than the shorter canines of conical-tooth cats, it has been long hypothesized that the heavily muscled forelimbs of sabretooths were used for immobilizing prey before developing a quick and precise killing bite. However, the effect of this unique adaptation on the covariation between the fore- and the hindlimb has not been explored in a quantitative fashion. In this paper, we investigate if the specialization of sabretooth predators decoupled the morphological variation of their forelimb with respect to their hindlimb or, in contrast, both limbs vary in the same fashion as in conical-tooth cats, which do not show such extreme adaptations in their forelimb. We use 3D geometric morphometrics and different morphological indices to compare the fore- and hindlimb of conical- and sabretooth predators. Our results indicate that the limb bones of sabretooth predators covary following the same trend of conical-tooth cats. Therefore, we show that the predatory specialization of sabretooth predators did not result in a decoupling of the morphological evolution of their fore- and hindlimbs. The role of developmental constraints and natural selection on this coordinate variation between the fore- and the hindlimb is discussed in the light of this new evidence.

Phylogenetic conservatism in skulls and evolutionary lability in limbs - morphological evolution across an ancient frog radiation is shaped by diet, locomotion and burrowing

Background

Quantifying morphological diversity across taxa can provide valuable insight into evolutionary processes, yet its complexities can make it difficult to identify appropriate units for evaluation. One of the challenges in this field is identifying the processes that drive morphological evolution, especially when accounting for shape diversification across multiple structures. Differential levels of co-varying phenotypic diversification can conceal selective pressures on traits due to morphological integration or modular shape evolution of different structures, where morphological evolution of different modules is explained either by co-variation between them or by independent evolution, respectively.

Methods

Here we used a 3D geometric morphometric approach with x-ray micro CT scan data of the skull and bones of forelimbs and hindlimbs of representative species from all 21 genera of the ancient Australo-Papuan myobatrachid frogs and analysed their shape both as a set of distinct modules and as a multi-modular integrative structure. We then tested three main questions: (i) are evolutionary patterns and the amount and direction of morphological changes similar in different structures and subfamilies?, (ii) do skulls and limbs show different levels of integration?, and (iii) is morphological diversity of skulls and limbs shaped by diet, locomotion, burrowing behavior, and ecology?.

Results

Our results in both skulls and limbs support a complex evolutionary pattern typical of an adaptive radiation with an early burst of phenotypic variation followed by slower rates of morphological change. Skull shape diversity was phylogenetically conserved and correlated with diet whereas limb shape was more labile and associated with diet, locomotion, and burrowing behaviour. Morphological changes between different limb bones were highly correlated, depicting high morphological integration. In contrast, overall limb and skull shape displayed semi-independence in morphological evolution, indicating modularity.

Conclusions

Our results illustrate how morphological diversification in animal clades can follow complex processes, entailing selective pressures from the environment as well as multiple trait covariance with varying degrees of independence across different structures. We suggest that accurately quantifying shape diversity across multiple structures is crucial in order to understand complex evolutionary processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0993-0) contains supplementary material, which is available to authorized users.

Morphological modularity in the vertebral column of Felidae (Mammalia, Carnivora)

BACKGROUND

Previous studies have demonstrated that the clear morphological differences among vertebrae across the presacral column are accompanied by heterogeneous functional signals in vertebral shape. Further, several lines of evidence suggest that the mammalian axial skeleton is a highly modular structure. These include its composition of serial units, a trade-off between high shape variance and strong conservation of vertebral count, and direct association of regions with anterior expression sites of Hox genes. Here we investigate the modular organisation of the presacral vertebral column of modern cats (Felidae, Carnivora, Mammalia) with pairwise comparisons of vertebral shape covariation (i.e. integration) and evaluate our results against hypotheses of developmental and functional modularity. We used three-dimensional geometric morphometrics to quantify vertebral shape and then assessed integration between pairs of vertebrae with phylogenetic two-block partial least square analysis (PLS).

RESULTS

Six modules were identified in the pairwise analyses (vertebrae included are designated as 'C' for cervical, 'T' for thoracic, and 'L' for lumbar): an anterior module (C1 to T1); a transitional module situated between the last cervicals and first thoracics (C6 to T2); an anterior to middle thoracic set (T4 to T8); an anticlinal module (T10 and T11); a posterior set composed of the last two thoracics and lumbars (T12 to L7); and a module showing covariation between the cervicals and the posterior set (T12 to L7). These modules reflect shared developmental pathways, ossification timing, and observed ecological shape diversification in living species of felids.

CONCLUSIONS

We show here that patterns of shape integration reflect modular organisation of the vertebral column of felids. Whereas this pattern corresponds with hypotheses of developmental and functional regionalisation in the axial skeleton, it does not simply reflect major vertebral regions. This modularity may also have permitted vertebral partitions, specifically in the posterior vertebral column, to be more responsive to selection and achieve higher morphological disparity than other vertebral regions.