Large‐Scale Diversification of Skull Shape in Domestic Dogs: Disparity and Modularity

Holm multiple correction for large-scale gene-shape association mapping

Background

Linkage Disequilibrium (LD) is a powerful approach for the identification and characterization of morphological shape, which usually involves multiple genetic markers. However, multiple testing corrections substantially reduce the power of the associated tests. In addition, the principle component analysis (PCA), used to quantify the shape variations into several principal phenotypes, further increases the number of tests. As a result, a powerful multiple testing correction for simultaneous large-scale gene-shape association tests is an essential part of determining statistical significance. Bonferroni adjustments and permutation tests are the most popular approaches to correcting for multiple tests within LD based Quantitative Trait Loci (QTL) models. However, permutations are extremely computationally expensive and may mislead in the presence of family structure. The Bonferroni correction, though simple and fast, is conservative and has low power for large-scale testing.

Results

We propose a new multiple testing approach, constructed by combining an Intersection Union Test (IUT) with the Holm correction, which strongly controls the family-wise error rate (FWER) without any additional assumptions on the joint distribution of the test statistics or dependence structure of the markers. The power improvement for the Holm correction, as compared to the standard Bonferroni correction, is examined through a simulation study. A consistent and moderate increase in power is found under the majority of simulated circumstances, including various sample sizes, Heritabilities, and numbers of markers. The power gains are further demonstrated on real leaf shape data from a natural population of poplar, Populus szechuanica var tietica, where more significant QTL associated with morphological shape are detected than under the previously applied Bonferroni adjustment.

Conclusion

The Holm correction is a valid and powerful method for assessing gene-shape association involving multiple markers, which not only controls the FWER in the strong sense but also improves statistical power.

Hedgehog signaling mediates adaptive variation in a dynamic functional system in the cichlid feeding apparatus

Adaptive variation in the craniofacial skeleton is a key component of resource specialization and habitat divergence in vertebrates, but the proximate genetic mechanisms that underlie complex patterns of craniofacial variation are largely unknown. Here we demonstrate that the Hedgehog (Hh) signaling pathway mediates widespread variation across a complex functional system that affects the kinematics of lower jaw depression--the opercular four-bar linkage apparatus--among Lake Malawi cichlids. By using a combined quantitative trait locus mapping and population genetics approach, we show that allelic variation in the Hh receptor, ptch1, affects the development of distinct bony elements in the head that represent two of three movable links in this functional system. The evolutionarily derived allele is found in species that feed from the water column, and is associated with shifts in anatomy that translate to a four-bar system capable of faster jaw rotation. Alternatively, the ancestral allele is found in species that feed on attached algae, and is associated with the development of a four-bar system that predicts slower jaw movement. Experimental manipulation of the Hh pathway during cichlid development recapitulates functionally salient natural variation in craniofacial geometry. In all, these results significantly extend our understanding of the mechanisms that fine-tune the craniofacial skeletal complex during adaptation to new foraging niches.

An analysis of sexual size dimorphism in goose

1. Sexual size dimorphism (SSD) is a common phenomenon in animals. Rensch's rule states that larger species generally exhibit a higher male to female body size ratio than smaller ones. 2. Domesticated animals offer excellent opportunities for testing predictions of the functional explanations of Rensch's rule and this was tested in a meta-analysis of SSD in 38 breeds of domestic geese compared among themselves and with their wild relatives (subfamily Anserinae, 35 species). 3. Domestic geese and wild Anser species taken together supported Rench's rule but the wild species did not. 4. The non-targeted sex selection hypothesis seems to provide the best intuitive explanation for the lack of SSD in geese.

Phenotypic integration of brain size and head morphology in Lake Tanganyika Cichlids

Background

Phenotypic integration among different anatomical parts of the head is a common phenomenon across vertebrates. Interestingly, despite centuries of research into the factors that contribute to the existing variation in brain size among vertebrates, little is known about the role of phenotypic integration in brain size diversification. Here we used geometric morphometrics on the morphologically diverse Tanganyikan cichlids to investigate phenotypic integration across key morphological aspects of the head. Then, while taking the effect of shared ancestry into account, we tested if head shape was associated with brain size while controlling for the potentially confounding effect of feeding strategy.

Results

The shapes of the anterior and posterior parts of the head were strongly correlated, indicating that the head represents an integrated morphological unit in Lake Tanganyika cichlids. After controlling for phylogenetic non-independence, we also found evolutionary associations between head shape, brain size and feeding ecology.

Conclusions

Geometric morphometrics and phylogenetic comparative analyses revealed that the anterior and posterior parts of the head are integrated, and that head morphology is associated with brain size and feeding ecology in Tanganyikan cichlid fishes. In light of previous results on mammals, our results suggest that the influence of phenotypic integration on brain diversification is a general process.

Genetic and morphometric evidence on a Galápagos Island exposes founder effects and diversification in the first-known (truly) feral western dog population

Domesticated animals that revert to a wild state can become invasive and significantly impact native biodiversity. Although dogs can be problematic locally, only the Australasian dingo is known to occur in isolation from humans. Western dogs have experienced more intense artificial selection, which potentially limits their invasiveness. However, feral dogs eradicated from Isabela Island, Galápagos in the 1980s could be the first-known exception. We used DNA and morphometric data from 92 of these dogs to test the hypotheses that (i) these dogs persisted independently of humans for up to a century and a half since descending from a handful of dogs introduced in the early 1800s, vs. (ii) similarly to other western feral dog populations, they reflected continuous recruitment of strays from human settlements on a portion of the Island. We detected one dominant maternal lineage and one dominant paternal lineage shared by the three subpopulations, along with low autosomal genetic diversity, consistent with the hypothesized common origins from a small founder population. Genetic diversity patterns among the three island subpopulations were consistent with stepping-stone founder effects, while morphometric differentiation suggested rapid phenotypic divergence, possibly due to drift and reinforced by selection corresponding to distinct microclimates and habitats on Isabela. Despite the continued presence of free-ranging dogs in the vicinity of settlements on Isabela and other Galápagos Islands, feral populations have not reestablished in remote areas since the 1980s, emphasizing the rarity of conditions necessary for feralization of modern western dogs.

Modularity of a leaf moth-wing pattern and a versatile characteristic of the wing-pattern ground plan

BACKGROUND

One of the most intriguing questions in evolutionary developmental biology is how an insect acquires a mimicry pattern within its body parts. A striking example of pattern mimicry is found in the pattern diversity of moth and butterfly wings, which is thought to evolve from preexisting elements illustrated by the nymphalid ground plan (NGP). Previous studies demonstrated that individuality of the NGP facilitates the decoupling of associated common elements, leading to divergence. In contrast, recent studies on the concept of modularity have argued the importance of a combination of coupling and decoupling of the constituent elements. Here, we examine the modularity of a mimicry wing pattern in a moth and explore an evolvable characteristic of the NGP.

RESULTS

This study examined the wings of the noctuid moth Oraesia excavata, which closely resemble leaves with a leaf venation pattern. Based on a comparative morphological procedure, we found that this leaf pattern was formed by the NGP common elements. Using geometric morphometrics combined with network analysis, we found that each of the modules in the leaf pattern integrates the constituent components of the leaf venation pattern (i.e., the main and lateral veins). Moreover, the detected modules were established by coupling different common elements and decoupling even a single element into different modules. The modules of the O. excavata wing pattern were associated with leaf mimicry, not with the individuality of the NGP common elements. For comparison, we also investigated the modularity of a nonmimetic pattern in the noctuid moth Thyas juno. Quantitative analysis demonstrated that the modules of the T. juno wing pattern regularly corresponded to the individuality of the NGP common elements, unlike those in the O. excavata wing pattern.

CONCLUSIONS

This study provides the first evidence for modularity in a leaf mimicry pattern. The results suggest that the evolution of this pattern involves coupling and decoupling processes to originate these modules, free from the individuality of the NGP system. We propose that this evolution has been facilitated by a versatile characteristic of the NGP, allowing the association of freely modifiable subordinate common elements to make modules.

The genetics of canine skull shape variation

A dog's craniofacial diversity is the result of continual human intervention in natural selection, a process that began tens of thousands of years ago. To date, we know little of the genetic underpinnings and developmental mechanisms that make dog skulls so morphologically plastic. In this Perspectives, we discuss the origins of dog skull shapes in terms of history and biology and highlight recent advances in understanding the genetics of canine skull shapes. Of particular interest are those molecular genetic changes that are associated with the development of distinct breeds.

Comparative Anatomy of the Bony Labyrinth (Inner Ear) of Placental Mammals

BACKGROUND

Variation is a naturally occurring phenomenon that is observable at all levels of morphology, from anatomical variations of DNA molecules to gross variations between whole organisms. The structure of the otic region is no exception. The present paper documents the broad morphological diversity exhibited by the inner ear region of placental mammals using digital endocasts constructed from high-resolution X-ray computed tomography (CT). Descriptions cover the major placental clades, and linear, angular, and volumetric dimensions are reported.

PRINCIPAL FINDINGS

The size of the labyrinth is correlated to the overall body mass of individuals, such that large bodied mammals have absolutely larger labyrinths. The ratio between the average arc radius of curvature of the three semicircular canals and body mass of aquatic species is substantially lower than the ratios of related terrestrial taxa, and the volume percentage of the vestibular apparatus of aquatic mammals tends to be less than that calculated for terrestrial species. Aspects of the bony labyrinth are phylogenetically informative, including vestibular reduction in Cetacea, a tall cochlear spiral in caviomorph rodents, a low position of the plane of the lateral semicircular canal compared to the posterior canal in Cetacea and Carnivora, and a low cochlear aspect ratio in Primatomorpha.

SIGNIFICANCE

The morphological descriptions that are presented add a broad baseline of anatomy of the inner ear across many placental mammal clades, for many of which the structure of the bony labyrinth is largely unknown. The data included here complement the growing body of literature on the physiological and phylogenetic significance of bony labyrinth structures in mammals, and they serve as a source of data for future studies on the evolution and function of the vertebrate ear.

Shape variation in the skull and lower carnassial in a wild population of raccoon dog (Nyctereutes procyonoides)

Individual variations in skull and lower carnassial morphology within a wild population of raccoon dog were examined using geometric morphometric techniques. We compared individual morphological variations by using relative warp analysis, and then tested morphological integration between the skull and carnassial by using partial least square (PLS) analysis. The most marked variation in skull shape was the dorsoventral flexion; i.e., deformation from klinorhynchy to airorhynchy. Two remarkable variations were observed, including tilting between the trigonid (or carnassial blade) and the talonid in the lower carnassial, and the relative sizes of the trigonid and the talonid. This observed variation in skull shape was similar to previous reports of variations among dog breeds that correlate with a polymorphism of the Runx2 gene. This polymorphism has also been reported to correlate with snout length, which is strongly related to carnivorous or omnivorous dietary adaptations, across the entire order Carnivora. Our results in the lower carnassial were also similar to previously reported patterns observed for carnivorous or omnivorous dietary adaptations among Carnivora. However, in our PLS analysis between skull and carnassial shapes, we only found a significant correlation in a lower dimension, suggesting a lower degree of integration. These results indicate that shape variations, which could be sources of natural selection in the skull and carnassial, were present in a wild population, suggesting high evolvability of these variations in the raccoon dog and the order Carnivora in general.

Patterns and implications of extensive heterochrony in carnivoran cranial suture closure

Heterochronic changes in the rate or timing of development underpin many evolutionary transformations. In particular, the onset and rate of bone development have been the focus of many studies across large clades. In contrast, the termination of bone growth, as estimated by suture closure, has been studied far less frequently, although a few recent studies have shown this to represent a variable, although poorly understood, aspect of developmental evolution. Here, we examine suture closure patterns across 25 species of carnivoran mammals, ranging from social-insectivores to hypercarnivores, to assess variation in suture closure across taxa, identify heterochronic shifts in a phylogenetic framework and elucidate the relationship between suture closure timing and ecology. Our results show that heterochronic shifts in suture closure are widespread across Carnivora, with several shifts identified for most major clades. Carnivorans differ from patterns identified for other mammalian clades in showing high variability of palatal suture closure, no correlation between size and level of suture closure, and little phylogenetic signal outside of musteloids. Results further suggest a strong influence of feeding ecology on suture closure pattern. Most of the species with high numbers of heterochronic shifts, such as the walrus and the aardwolf, feed on invertebrates, and these taxa also showed high frequency of closure of the mandibular symphysis, a state that is relatively rare among mammals. Overall, caniforms displayed more heterochronic shifts than feliforms, suggesting that evolutionary changes in suture closure may reflect the lower diversity of cranial morphology in feliforms.

Quantitative genetic analysis of subspecific differences in body shape in the snail-feeding carabid beetle Damaster blaptoides

A dimorphic pattern of macrocephalic (wide, short) and stenocephalic (narrow, long) body shapes is observed in snail-feeding carabid beetles globally. The former exhibits high performance in crushing snail shells with powerful jaws, whereas the latter specializes in eating snails' soft body directly by inserting the head into the shell. In the snail-feeding species Damaster blaptoides, the subspecies D. b. capito has a wide, short forebody, and D. b. fortunei has a narrow, long forebody. They exhibit distinct morphologies despite their geographic and phylogenetic proximity. To examine the genetic basis of the morphological differences between these two subspecies, we conducted quantitative genetic analyses by crossing these subspecies and producing F(1) and backcross hybrids. The hybrids had body shapes intermediate between the parental subspecies. The variation between wide, short and narrow, long forebodies was based on negative genetic correlations between width and length of the head and thorax. Between one and eight genetic factors were involved in the morphological differences between subspecies. We suggest that the morphological integration of forebody parts in a small number of loci has facilitated the marked morphological diversification between subspecies of D. blaptoides.

Allometry of sexual size dimorphism in domestic dog

BACKGROUND

The tendency for male-larger sexual size dimorphism (SSD) to scale with body size - a pattern termed Rensch's rule - has been empirically supported in many animal lineages. Nevertheless, its theoretical elucidation is a subject of debate. Here, we exploited the extreme morphological variability of domestic dog (Canis familiaris) to gain insights into evolutionary causes of this rule.

METHODOLOGY/PRINCIPAL FINDINGS

We studied SSD and its allometry among 74 breeds ranging in height from less than 19 cm in Chihuahua to about 84 cm in Irish wolfhound. In total, the dataset included 6,221 individuals. We demonstrate that most dog breeds are male-larger, and SSD in large breeds is comparable to SSD of their wolf ancestor. Among breeds, SSD becomes smaller with decreasing body size. The smallest breeds are nearly monomorphic.

CONCLUSIONS/SIGNIFICANCE

SSD among dog breeds follows the pattern consistent with Rensch's rule. The variability of body size and corresponding changes in SSD among breeds of a domestic animal shaped by artificial selection can help to better understand processes leading to emergence of Rensch's rule.

The generation of variation and the developmental basis for evolutionary novelty

Organisms exhibit an incredible diversity of form, a fact that makes the evolution of novelty seemingly self-evident. However, despite the "obvious" case for novelty, defining this concept in evolutionary terms is highly problematic, so much so that some have suggested discarding it altogether. Approaches to this problem tend to take either an adaptation- or development-based perspective, but we argue here that an exclusive focus on either of these misses the original intent of the novelty concept and undermines its practical utility. We propose instead that for a feature to be novel, it must have evolved both by a transition between adaptive peaks on the fitness landscape and that this transition must have overcome a previous developmental constraint. This definition focuses novelty on the explanation of apparently difficult or low-probability evolutionary transitions and highlights how the integration of developmental and functional considerations are necessary to evolutionary explanation. It further reinforces that novelty is a central concern not just of evolutionary developmental biology (i.e., "evo-devo") but of evolutionary biology more generally. We explore this definition of novelty in light of four examples that range from the obvious to subtle.

Rethinking dog domestication by integrating genetics, archeology, and biogeography

The dog was the first domesticated animal but it remains uncertain when the domestication process began and whether it occurred just once or multiple times across the Northern Hemisphere. To ascertain the value of modern genetic data to elucidate the origins of dog domestication, we analyzed 49,024 autosomal SNPs in 1,375 dogs (representing 35 breeds) and 19 wolves. After combining our data with previously published data, we contrasted the genetic signatures of 121 breeds with a worldwide archeological assessment of the earliest dog remains. Correlating the earliest archeological dogs with the geographic locations of 14 so-called "ancient" breeds (defined by their genetic differentiation) resulted in a counterintuitive pattern. First, none of the ancient breeds derive from regions where the oldest archeological remains have been found. Second, three of the ancient breeds (Basenjis, Dingoes, and New Guinea Singing Dogs) come from regions outside the natural range of Canis lupus (the dog's wild ancestor) and where dogs were introduced more than 10,000 y after domestication. These results demonstrate that the unifying characteristic among all genetically distinct so-called ancient breeds is a lack of recent admixture with other breeds likely facilitated by geographic and cultural isolation. Furthermore, these genetically distinct ancient breeds only appear so because of their relative isolation, suggesting that studies of modern breeds have yet to shed light on dog origins. We conclude by assessing the limitations of past studies and how next-generation sequencing of modern and ancient individuals may unravel the history of dog domestication.

The developmental renaissance in adaptationism

From an adaptation perspective, unoccupied patches of morphological space are inferred to be empty because they are of low fitness and selected against. These inferences hinge on venturesome assumptions, because emptiness is explained by low fitness and low fitness is inferred from emptiness. Moreover, non-adaptive factors, such as developmental constraint, could also plausibly account for empty morphospace. In response, biologists increasingly study ontogeny to test the assumption that unobserved phenotypes could be produced if selection were to favor them; finding that empty space morphologies can be readily produced in development helps reject constraint and lends support to adaptive hypotheses. This developmental approach to adaptation calls on manifold techniques, including embryology, artificial selection and comparative methods. Belying their diversity, all of these methods examine the causes of empty morphospace and mark a return of development, long excluded from traditional evolutionary biology, to adaptationist practice.

Pervasive genetic integration directs the evolution of human skull shape

It has long been unclear whether the different derived cranial traits of modern humans evolved independently in response to separate selection pressures or whether they resulted from the inherent morphological integration throughout the skull. In a novel approach to this issue, we combine evolutionary quantitative genetics and geometric morphometrics to analyze genetic and phenotypic integration in human skull shape. We measured human skulls in the ossuary of Hallstatt (Austria), which offer a unique opportunity because they are associated with genealogical data. Our results indicate pronounced covariation of traits throughout the skull. Separate simulations of selection for localized shape changes corresponding to some of the principal derived characters of modern human skulls produced outcomes that were similar to each other and involved a joint response in all of these traits. The data for both genetic and phenotypic shape variation were not consistent with the hypothesis that the face, cranial base, and cranial vault are completely independent modules but relatively strongly integrated structures. These results indicate pervasive integration in the human skull and suggest a reinterpretation of the selective scenario for human evolution where the origin of any one of the derived characters may have facilitated the evolution of the others.

Developmental plasticity, morphological variation and evolvability: a multilevel analysis of morphometric integration in the shape of compound leaves

The structure of compound leaves provides flexibility for morphological change by variation in the shapes, sizes and arrangement of leaflets. Here, we conduct a multilevel analysis of shape variation in compound leaves to explore the developmental plasticity and evolutionary potential that are the basis of diversification in leaf shape. We use the methods of geometric morphometrics to study the shapes of individual leaflets and whole leaves in 20 taxa of Potentilla (sensu lato). A newly developed test based on the bootstrap approach suggests that uncertainty in the molecular phylogeny precludes firm conclusions whether there is a phylogenetic signal in the data on leaf shape. For variation among taxa, variation within taxa, as well as fluctuating asymmetry, there is evidence of strong morphological integration. The patterns of variation are similar across all three levels, suggesting that integration within taxa may act as a constraint on evolutionary change.

Beyond bilateral symmetry: geometric morphometric methods for any type of symmetry

BACKGROUND

Studies of symmetric structures have made important contributions to evolutionary biology, for example, by using fluctuating asymmetry as a measure of developmental instability or for investigating the mechanisms of morphological integration. Most analyses of symmetry and asymmetry have focused on organisms or parts with bilateral symmetry. This is not the only type of symmetry in biological shapes, however, because a multitude of other types of symmetry exists in plants and animals. For instance, some organisms have two axes of reflection symmetry (biradial symmetry; e.g. many algae, corals and flowers) or rotational symmetry (e.g. sea urchins and many flowers). So far, there is no general method for the shape analysis of these types of symmetry.

RESULTS

We generalize the morphometric methods currently used for the shape analysis of bilaterally symmetric objects so that they can be used for analyzing any type of symmetry. Our framework uses a mathematical definition of symmetry based on the theory of symmetry groups. This approach can be used to divide shape variation into a component of symmetric variation among individuals and one or more components of asymmetry. We illustrate this approach with data from a colonial coral that has ambiguous symmetry and thus can be analyzed in multiple ways. Our results demonstrate that asymmetric variation predominates in this dataset and that its amount depends on the type of symmetry considered in the analysis.

CONCLUSIONS

The framework for analyzing symmetry and asymmetry is suitable for studying structures with any type of symmetry in two or three dimensions. Studies of complex symmetries are promising for many contexts in evolutionary biology, such as fluctuating asymmetry, because these structures can potentially provide more information than structures with bilateral symmetry.

Evolutionary patterns and processes in the radiation of phyllostomid bats

Background

The phyllostomid bats present the most extensive ecological and phenotypic radiation known among mammal families. This group is an important model system for studies of cranial ecomorphology and functional optimisation because of the constraints imposed by the requirements of flight. A number of studies supporting phyllostomid adaptation have focused on qualitative descriptions or correlating functional variables and diet, but explicit tests of possible evolutionary mechanisms and scenarios for phenotypic diversification have not been performed. We used a combination of morphometric and comparative methods to test hypotheses regarding the evolutionary processes behind the diversification of phenotype (mandible shape and size) and diet during the phyllostomid radiation.

Results

The different phyllostomid lineages radiate in mandible shape space, with each feeding specialisation evolving towards different axes. Size and shape evolve quite independently, as the main directions of shape variation are associated with mandible elongation (nectarivores) or the relative size of tooth rows and mandibular processes (sanguivores and frugivores), which are not associated with size changes in the mandible. The early period of phyllostomid diversification is marked by a burst of shape, size, and diet disparity (before 20 Mya), larger than expected by neutral evolution models, settling later to a period of relative phenotypic and ecological stasis. The best fitting evolutionary model for both mandible shape and size divergence was an Ornstein-Uhlenbeck process with five adaptive peaks (insectivory, carnivory, sanguivory, nectarivory and frugivory).

Conclusions

The radiation of phyllostomid bats presented adaptive and non-adaptive components nested together through the time frame of the family's evolution. The first 10 My of the radiation were marked by strong phenotypic and ecological divergence among ancestors of modern lineages, whereas the remaining 20 My were marked by stasis around a number of probable adaptive peaks. A considerable amount of cladogenesis and speciation in this period is likely to be the result of non-adaptive allopatric divergence or adaptations to peaks within major dietary categories.

MorphoJ: an integrated software package for geometric morphometrics

Increasingly, data on shape are analysed in combination with molecular genetic or ecological information, so that tools for geometric morphometric analysis are required. Morphometric studies most often use the arrangements of morphological landmarks as the data source and extract shape information from them by Procrustes superimposition. The MorphoJ software combines this approach with a wide range of methods for shape analysis in different biological contexts. The program offers an integrated and user-friendly environment for standard multivariate analyses such as principal components, discriminant analysis and multivariate regression as well as specialized applications including phylogenetics, quantitative genetics and analyses of modularity in shape data. MorphoJ is written in Java and versions for the Windows, Macintosh and Unix/Linux platforms are freely available from http://www.flywings.org.uk/MorphoJ_page.htm.