Effects of environmental disturbance on phenotypic variation: an integrated assessment of canalization, developmental stability, modularity, and allometry in lizard head shape

Skull asymmetry in various sheep breeds: Directional asymmetry and fluctuating asymmetry

Sheep (Ovis aries) play an important role in the economy of Turkey and the Balkan Peninsula due to their use in farming. As a domesticated species, sheep's morphometric and morphological diversity is likely determined by selective breeding practices rather than geographic distribution. This study aimed to analyse four different sheep breed skulls and reveal skull asymmetry using geometric morphometric methods. For this purpose, 2D images of 52 sheep skulls from different breeds were analysed from the dorsal view of the skull, using 28 landmarks. In the comparison of sheep skulls from the dorsal view, the first principal components for directional asymmetry (DA) and fluctuating asymmetry (FA) were 32.98% and 39.62% of the total variation, respectively. Sharri and Ivesi (Awassi) sheep breeds had the broadest distribution of skull shapes among the breeds, while Lara e Polisit was the most conservative breed. DA was used as a measure of biomechanical constraints, and FA was used as an indicator of environmental stress. Consistent with DA, both differences in centroid size and shape between breeds were statistically significant. No differences between males and females related to asymmetry were revealed. Ivesi sheep revealed the highest fluctuating asymmetry. Geometric morphometric methods proved to be a useful tool for distinguishing differences in the shape of the skull of different sheep breeds and also can be useful for taxonomic purposes.

Rapid and biased evolution of canalization during adaptive divergence revealed by dominance in gene expression variability during Arctic charr early development

Adaptive evolution may be influenced by canalization, the buffering of developmental processes from environmental and genetic perturbations, but how this occurs is poorly understood. Here, we explore how gene expression variability evolves in diverging and hybridizing populations, by focusing on the Arctic charr (Salvelinus alpinus) of Thingvallavatn, a classic case of divergence between feeding habitats. We report distinct profiles of gene expression variance for both coding RNAs and microRNAs between the offspring of two contrasting morphs (benthic/limnetic) and their hybrids reared in common conditions and sampled at two key points of cranial development. Gene expression variance in the hybrids is substantially affected by maternal effects, and many genes show biased expression variance toward the limnetic morph. This suggests that canalization, as inferred by gene expression variance, can rapidly diverge in sympatry through multiple gene pathways, which are associated with dominance patterns possibly biasing evolutionary trajectories and mitigating the effects of hybridization on adaptive evolution.

Is Fluctuating Asymmetry a Sufficient Indicator of Stress Level in Two Lizard Species (Zootoca vivipara and Lacerta agilis) from Alpine Habitats?

Alpine habitats are exposed to increasing anthropogenic pressure and climate change. The negative impacts can lead to chronic stress that can affect the survival and reproductive success of individuals and even lead to population extinction. In this study, we analyse different morphological and ecological traits and indices of abiotic and biotic stressors (such as head size and shape, fluctuating asymmetry, body condition index, tail autotomy, and population abundance) in alpine and subalpine populations of two lacertid species (Zootoca vivipara and Lacerta agilis) from Serbia and North Macedonia. These lizards live under different conditions: allotopy/syntopy, different anthropogenic pressure, and different levels of habitat protection. We found differences between syntopic and allotopic populations in pileus size, body condition index (in both species), pileus shape, fluctuating asymmetry (in L. agilis), and abundance (in Z. vivipara). Differences between populations under anthropogenic pressure and populations without it were observed in pileus shape, body condition index (in both species), pileus size, fluctuating asymmetry, tail autotomy and abundance (in L. agilis). On the basis of our results, it is necessary to include other stress indicators in addition to fluctuating asymmetry to quickly observe and quantify the negative effects of threat factors and apply protective measures.

Effects of the egg incubation environment on turtle carapace development

Developing organisms are often exposed to fluctuating environments that destabilize tissue-scale processes and induce abnormal phenotypes. This might be common in species that lay eggs in the external environment and with little parental care, such as many reptiles. In turtles, morphological development has provided striking examples of abnormal phenotypic patterns, though the influence of the environment remains unclear. To this end, we compared fluctuating asymmetry, as a proxy for developmental instability, in turtle hatchlings incubated in controlled laboratory and unstable natural conditions. Wild and laboratory hatchlings featured similar proportions of supernumerary scales (scutes) on the dorsal shell (carapace). Such abnormal scutes likely elevated shape asymmetry, which was highest in natural nests. Moreover, we tested the hypothesis that hot and dry environments cause abnormal scute formation by subjecting eggs to a range of hydric and thermal laboratory incubation regimes. Shape asymmetry was similar in hatchlings incubated at five constant temperatures (26-30°C). A hot (30°C) and severely Dry substrate yielded smaller hatchlings but scutes were not overtly affected. Our study suggests that changing nest environments contribute to fluctuating asymmetry in egg-laying reptiles, while clarifying the conditions at which turtle shell development remains buffered from the external environment.

Developmental instability in domesticated mammals

Measures of fluctuating asymmetry (FA) have been adopted widely as an estimate of developmental instability. Arising from various sources of stress, developmental instability is associated with an organism's capacity to maintain fitness. The process of domestication has been framed as an environmental stress with human-specified parameters, suggesting that FA may manifest to a larger degree among domesticates compared to their wild relatives. This study used three-dimensional geometric morphometric landmark data to (a) quantify the amount of FA in the cranium of six domestic mammal species and their wild relatives and, (b) provide novel assessment of the commonalities and differences across domestic/wild pairs concerning the extent to which random variation arising from the developmental system assimilates into within-population variation. The majority of domestic mammals showed greater disparity for asymmetric shape, however, only two forms (Pig, Dog) showed significantly higher disparity as well as a higher degree of asymmetry compared to their wild counterparts (Wild Boar, Wolf). Contra to predictions, most domestic and wild forms did not show a statistically significant correspondence between symmetric shape variation and FA, however, a moderate correlation value was recorded for most pairs (r-partial least squares >0.5). Within pairs, domestic and wild forms showed similar correlation magnitudes for the relationship between the asymmetric and symmetric components. In domesticates, new variation may therefore retain a general, conserved pattern in the gross structuring of the cranium, whilst also being a source for response to selection on specific features.

Condition- and context-dependent variation of sexual dimorphism across lizard populations at different spatial scales

The evolution of sexual dimorphism (SD) is driven by intricate interplays between sexual and natural selection. When it comes to SD variation within populations, however, environmental factors play a major role. Sexually selected traits are expected to be strongly dependent on individual body condition, which is influenced by the local environment that individuals experience. As a consequence, the degree of SD may also depend on resource availability. Here, we investigated the potential drivers of SD expression at two sexually dimorphic morphometric traits, body size (snout vent length) and head shape (head geometric morphometrics), in the Italian wall lizard (Podarcis siculus). We assessed the existence of condition- and context-dependent SD across ten islands of the Aeolian archipelago (southern Italy), at within- and among-population scales. We observed strong geographical variation of SD among islands, and tested three potential SD predictors related to resource availability (individual body condition, ecosystem productivity, temperature). Body condition and ecosystem productivity were the main drivers of body size SD variation, and body condition was also the main driver for head shape SD. Our results highlight that the expression of SD in the Italian wall lizard is both condition- and context-dependent. These results are congruent at within- and among-populations scales highlighting that spatial multi-scale analysis represents a useful approach to understand patterns of SD expression.

Anatomical network analyses reveal evolutionary integration and modularity in the lizards skull

The morphology of lizard skulls is highly diverse, and it is crucial to understand the factors that constrain and promote their evolution to understand how lizards thrive. The results of interactions between cranial bones reflecting these factors can be detected as integration and modularity, and the analysis of integration and modularity allows us to explore the underlying factors. In this study, the integration and modularity of the skulls of lizards and the outgroup tuatara are analyzed using a new method, Anatomical Network Analysis (AnNA), and the factors causing lizards morphological diversity are investigated by comparing them. The comparison of modular structures shows that lizard skulls have high integration and anisomerism, some differences but basically common modular patterns. In contrast, the tuatara shows a different modular pattern from lizards. In addition, the presence of the postorbital bar by jugal and postorbital (postorbitofrontal) also reflect various functional factors by maintaining low integration. The maintenance of basic structures due to basic functional requirements and changes in integration within the modules play a significant role in increasing the morphological diversity of the lizard skull and in the prosperity of the lizards.

Native Lizards Living in Brazilian Cities: Effects of Developmental Environments on Thermal Sensitivity and Morpho-Functional Associations of Locomotion

Environmental conditions often affect developmental processes and consequently influence the range of phenotypic variation expressed at population level. Expansion of urban sites poses new challenges for native species, as urbanization usually affects the intensity of solar exposure and shade availability, determining the thermal regimes organisms are exposed to. In this study, we evaluate the effects of different developmental conditions in a Tropidurus lizard commonly found in Brazilian urban sites. After incubating embryos of Tropidurus catalanensis in two different thermal regimes (Developmental Environments [DE]: cold 24°C and warm 30°C), we measured morphological traits in the neonates and quantified locomotor performance in horizontal and vertical surfaces at three temperatures [Test Temperatures (TT) = 24°C, 30°C and 36°C]. Results indicate effects of developmental temperatures on morphological features, expressing functional implications that might be decisive for the viability of T. catalanensis in urbanized areas. Lizards ran similarly on horizontal and vertical surfaces, and isolated analyses did not identify significant effects of DE or TT on the sprint speeds measured. Absolute Vmax (i.e., the maximum sprint speed reached among all TTs) positively correlated with body size (SVL), and neonates from the warm DE (30°C) were larger than those from the cold DE (24°C). Morpho-functional associations of absolute Vmax also involved pelvic girdle width and forelimb, hindlimb, trunk, and tail lengths. Emerging discussions aim to understand how animals cope with abrupt environmental shifts, a likely common challenge in urbanized sites. Our findings add a new dimension to the topic, providing evidence that temperature, an environmental parameter often affected by urbanization, influences the thermal sensitivity of locomotion and the morphological profile of T. catalanensis neonates. Thermal sensitivity of specific developmental processes may influence the ability of these lizards to remain in habitats that change fast, as those suffering rapid urbanization due to city growth.

Phenotypic variation in urban environments: mechanisms and implications

In the past decade, numerous studies have explored how urbanisation affects the mean phenotypes of populations, but it remains unknown how urbanisation impacts phenotypic variation, a key target of selection that shapes, and is shaped by, eco-evolutionary processes. Our review suggests that urbanisation may often increase intraspecific phenotypic variation through several processes; a conclusion aligned with results from our illustrative analysis on tit morphology across 13 European city/forest population pairs. Urban-driven changes in phenotypic variation will have immense implications for urban populations and communities, particularly through urbanisation's effects on individual fitness, species interactions, and conservation. We call here for studies that incorporate phenotypic variation in urban eco-evolutionary research alongside advances in theory.

Effects of biotic and abiotic stressors on asymmetries and head size in two sympatric lizard species

Organisms face numerous environmental stressors, which can affect developmental precision, including symmetry of various physical characteristics. Fluctuating asymmetry (FA) has therefore been suggested as a simple and efficient tool for assessing sub-lethal stress levels. We analyzed FA in two sympatric lizard species (Iberolacerta horvathi and Podarcis muralis) to determine potential effects of interspecific competition and urbanization, as proxies of stress, taking into account sexual dimorphism and environmental conditions. We sampled 16 syntopic and allotopic populations and used geometric morphometrics of head morphology. We detected significant but mixed effects on the head asymmetry from the environment and the syntopic occurrence that differed between species. P. muralis lizards had more asymmetric heads at higher altitudes, while I. horvathi lizards did at mid altitudes, which may be explained by P. muralis experiencing environmental stress of colder conditions at higher altitudes. The mid-altitude effect on asymmetries in I. horvathi might be explained by a lower availability of stony walls and higher abundance of P. muralis, thus higher competition. The asymmetry of supraciliary granules was affected by the presence of other species. However, lizards from allotopic populations attained larger asymmetries compared to lizards from syntopic populations, which was the opposite from what was expected. There was no effect of urbanization in P. muralis, which could be due to relatively low pollution and habitat degradation in study locations. Overall, we highlighted the possibility of using lizards and FA for bioindication of environmental stressors and especially improved the knowledge gap in the research of biotic stressors.

A New Integrated Tool to Calculate and Map Bilateral Asymmetry on Three-Dimensional Digital Models

The observation and the quantification of asymmetry in biological structures are deeply investigated in geometric morphometrics. Patterns of asymmetry were explored in both living and fossil species. In living organisms, levels of directional and fluctuating asymmetry are informative about developmental processes and health status of the individuals. Paleontologists are primarily interested in asymmetric features introduced by the taphonomic process, as they may significantly alter the original shape of the biological remains, hampering the interpretation of morphological features which may have profound evolutionary significance. Here, we provide a new R tool that produces the numerical quantification of fluctuating and directional asymmetry and charts asymmetry directly on the specimens under study, allowing the visual inspection of the asymmetry pattern. We tested this show.asymmetry algorithm, written in the R language, on fossil and living cranial remains of the genus Homo. show.asymmetry proved successful in discriminating levels of asymmetry among sexes in Homo sapiens, to tell apart fossil from living Homo skulls, to map effectively taphonomic distortion directly on the fossil skulls, and to provide evidence that digital restoration obliterates natural asymmetry to unnaturally low levels.

Does stress mess with rodents' heads? Influence of habitat amount and genetic factors in mandible fluctuating asymmetry in South American water rats (Nectomys squamipes, Sigmodontinae) from Brazilian Atlantic rainforest remnants

Loss of developmental stability can lead to deviations from bilateral symmetry (i.e. Fluctuating Asymmetry - FA), and is thought to be caused by environmental and genetic factors associated with habitat loss and stress. Therefore, levels of FA might be a valuable tool to monitor wild populations if FA serves as an indicator of exposure to stress due to impacts of habitat loss and fragmentation. In studies examining FA and habitat fragmentation, FA levels are often explained by loss of genetic variation, though few studies have addressed FA's use as indicator of environmental impact. Here, we investigated whether habitat loss, genetic variation, and/or inbreeding affect the developmental instability in Brazilian Atlantic forest populations of a Neotropical water rat (Nectomys squamipes). We sampled individuals from eight sites within Atlantic forest remnants with different amounts of available forest habitat and assessed FA levels with geometric morphometric techniques using adult mandibles. We used observed heterozygosity (Ho) and inbreeding coefficient (Fis), from seven microsatellite markers, as a proxy of genetic variation at individual and population levels. Populations were not significantly different for shape or size FA levels. Furthermore, interindividual variation in both shape and size FA levels and interpopulational differences in size FA levels were best explained by chance. However, habitat amount was negatively associated with both interpopulational variance and average shape FA levels. This association was stronger in populations living in areas with <28% of forest cover, which presented higher variance and higher average FA, suggesting that Nectomys squamipes might have a tolerance threshold to small availability of habitat. Our work is one of the first to use FA to address environmental stress caused by habitat loss in small mammal populations from a Neotropical biome. We suggest that shape FA might serve as a conservation tool to monitor human impact on natural animal populations.

Fluctuating asymmetry as biomarker of pesticides exposure in the Italian wall lizards (Podarcis siculus)

The extensive use of pesticides in agricultural environments produces drastic effects on wildlife, hence the need for less invasive indicators of environmental stress to monitoring the impact of agriculture treatments on biological systems. Fluctuating asymmetry (FA), as measure of developmental instability, has recently been proposed as reliable biomarker of populations stress due to environmental disturbance. We investigate femoral pores (FP) and dorsal head shape (HS) traits in populations of the Italian wall lizard inhabiting agricultural environments to examine whether different pesticide exposures (conventional, organic and control) can cause distinctive degree of FA. High-resolution photographs of FP and HS were taken in the field with a digital camera. The number of FP were counted twice on both sides and HS was analysed using geometric morphometrics with 25 landmarks and 12 semilandmarks. Individuals under conventional management showed higher levels of FA compared to control ones, and females exhibited higher FA levels than males for the FP. However, no significant difference was found for the HS trait. Our study provided evidence that FA may have a real potential as biomarker of population stress in wall lizards, highlighting the importance in the choice of the experimental design and the traits adopted for estimating DI.

Temporal instability of lake charr phenotypes: Synchronicity of growth rates and morphology linked to environmental variables?

Pathways through which phenotypic variation among individuals arise can be complex. One assumption often made in relation to intraspecific diversity is that the stability or predictability of the environment will interact with expression of the underlying phenotypic variation. To address biological complexity below the species level, we investigated variability across years in morphology and annual growth increments between and within two sympatric lake charr Salvelinus namaycush ecotypes in Rush Lake, USA. A rapid phenotypic shift in body and head shape was found within a decade. The magnitude and direction of the observed phenotypic change were consistent in both ecotypes, which suggests similar pathways caused the variation over time. Over the same time period, annual growth increments declined for both lake charr ecotypes and corresponded with a consistent phenotypic shift of each ecotype. Despite ecotype-specific annual growth changes in response to winter conditions, the observed annual growth shift for both ecotypes was linked, to some degree, with variation in the environment. Particularly, a declining trend in regional cloud cover was associated with an increase of early-stage (ages 1-3) annual growth for lake charr of Rush Lake. Underlying mechanisms causing changes in growth rates and constrained morphological modulation are not fully understood. An improved knowledge of the biology hidden within the expression of phenotypic variation promises to clarify our understanding of temporal morphological diversity and instability.

How do developmental and parental exposures to predation affect personality and immediate behavioural plasticity in the snail Physa acuta?

Individuals differ in personality and immediate behavioural plasticity. While developmental environment may explain this group diversity, the effect of parental environment is still unexplored-a surprising observation since parental environment influences mean behaviour. We tested whether developmental and parental environments impacted personality and immediate plasticity. We raised two generations of Physa acuta snails in the laboratory with or without developmental exposure to predator cues. Escape behaviour was repeatedly assessed on adult snails with or without predator cues in the immediate environment. On average, snails were slower to escape if they or their parents had been exposed to predator cues during development. Snails were also less plastic in response to immediate predation risk on average if they or their parents had been exposed to predator cues. Group diversity in personality was greater in predator-exposed snails than unexposed snails, while parental environment did not influence it. Group diversity in immediate plasticity was not significant. Our results suggest that only developmental environment plays a key role in the emergence of group diversity in personality, but that parental environment influences mean behavioural responses to the environmental change. Consequently, although different, both developmental and parental cues may have evolutionary implications on behavioural responses.

When European meets African honeybees (Apis mellifera L.) in the tropics: Morphological changes related to genetics in Mauritius Island (South-West Indian Ocean)

The previous genetic characterization of the honeybee population of Mauritius Island (Indian Ocean) revealed an ongoing process of hybridization between the first established African subspecies Apis mellifera unicolor and recently imported European subspecies (A. m. ligustica, A. m. carnica and A. m. mellifera). This context offers the rare opportunity to explore the influence of hybridization between African and European honeybees on phenotypic traits out of the case largely studied of the Africanized honeybee (hybrid between A. m. scutellata from South Africa and European subspecies). We thus conducted geometric morphometric analyses on forewings of 283 workers genetically characterized at 14 microsatellite loci to evaluate (1) if the morphological variability coincides well with the neutral genetic variability, (2) if hybrids exhibited rather parental, intermediate or transgressive traits, and (3) to test if fluctuating asymmetry (FA) of size and shape, as a measure of developmental stability, was elevated in hybrids (due to genetic stress) and/or European bees (due to unsuitable environment) compared to African bees. A strong concordance was found between morphological variability and neutral genetic variability, especially for wing shape, based on partial least-square analyses (PLS). However, on average, the morphology of hybrids was more similar to the African bees, potentially reflecting the dynamics and direction of introgression. Significant FA for wing size as well as wing shape was detected, suggesting the overall presence of stress during the development of the studied individuals. In contrast, the asymmetry levels do not differ according to the ancestry (African, European or hybrid) of the individuals. Therefore, if ongoing hybridization contributed to increasing the genetic and phenotypic diversity of the populations and influences its adaptive potential, developmental stressors could not be identified and their evolutionary consequences remain uncertain.

Geographical differences in competitive hierarchy in a native–invasive system

Biological invasions can create novel competitive interactions and force ecological shifts in both native and invasive species. Anuran tadpoles are able to modify their behaviour, morphology, growth and development to cope with competitive pressure. This plasticity is a good target for natural selection and can drive rapid evolutionary changes in response to novel interactions. Here, we explore changes in plastic responses and fitness of competing invasive and native tadpoles by exposing tadpoles from different locations with contrasting evolutionary histories to the same set of varied competitive conditions. Eggs were collected from one site near the first introduction of the invasive frog (~110 years of coexistence) and from a second site that was invaded recently. We hypothesized less favourable outcomes for the invasive species in long-coexisting populations, where the native competitor might have developed adaptive responses. Most results support the hypothesis. Where the invasion was older, invasive tadpoles exposed to native competitors grew less, developed more slowly and displayed morphologies linked to competitive stress, whereas the developmental stability and canalization of native tadpoles increased. On the whole, the asymmetric competitive relationship thus appeared to approach symmetry after ~35 generations, highlighting a noteworthy example of rapid adaptation after an invasion.

Righting behaviour in the European pond turtle (Emys orbicularis): relations between behavioural and morphological lateralization

Lateralization represents a key property of many behavioural traits, with the right and left sides of the brain providing different and integrative functions. Common ecological contexts where lateralization can be observed are foraging and predatory ones, where both visual and auditory lateralization may provide advantages such as faster response and increasing neural processing capacity. This is crucial in selecting a safe refuge during a predatory attack and may strongly affect the outcome of predator-prey interactions. For animals like turtles, a critical condition may occur when they are overturned on their own shell, which causes difficulties in breathing and thermoregulation, making them more vulnerable to predators. Therefore, the ability to right is a critical adaptive component related to survival in aquatic turtles, which has been observed to be lateralized. However, an overlooked feature of behavioural lateralization is its possible relationship with asymmetry in external morphology. Here we investigated this topic in freshwater European pond turtles Emys orbicularis, looking at a possible relation between lateralization in righting behaviour response and asymmetry in the shape of turtles' plastron and carapace. Righting performance (total time needed to completely turn) appeared to depend on shell shape. We found that none of the morphometric variables was related to a lateralization index calculated as the first side from which turtles tried to right. However, a strong negative correlation between the asymmetry index of plastron and the turning direction emerged, with more symmetric animals tending to turn to the right side.

Genetic and morphological variation of metacercariae of Microphallus piriformes (Trematoda, Microphallidae): Effects of paraxenia and geographic location

Host organism offers an environment for a parasite, and this environment is heterogenous within the host, variable among individual as well as between the hosts, and changing during the host's lifetime. This heterogeneity may act as a prerequisite for parasite species divergence. Intraspecific variability related to a certain type of heterogeneity may indicate an initial stage of speciation, and thus poses an evolutionary importance. Here we analyzed genetic and morphologic variation of trematode metacercariae of Microphallus piriformes (Trematoda, Microphallidae). Genetic variability of trematodes was assessed from sequences of cytochrome c oxidase subunit 1 (COI) and internal transcribed spacer region (ITS-1). Morphological variation of metacercarial body shape was for the first time analyzed using geometric morphometrics. Parasites from the White Sea and the Barents Sea coasts demonstrated partial genetic divergence (according to COI sequence analysis) and had significantly different body shape. Neither genetic nor morphological variation of metacercariae was related to intermediate host species. We discuss possible causes of the observed genetic divergence of parasite populations in different geographic regions.

Town and Country Reptiles: A Review of Reptilian Responses to Urbanization

The majority of the world population is now inhabiting urban areas, and with staggering population growth, urbanization is also increasing. While the work studying the effects of changing landscapes and specific urban pressures on wildlife is beginning to amass, the majority of this work focuses on avian or mammalian species. However, the effects of urbanization likely vary substantially across taxonomic groups due to differences in habitat requirements and life history. The current article aims first to broaden the review of urban effects across reptilian species; second, to summarize the responses of reptilian fauna to specific urban features; and third, to assess the directionality of individual and population level responses to urbanization in reptile species. Based on our findings, urban research in reptilian taxa is lacking in the following areas: (1) investigating interactive or additive urban factors, (2) measuring multiple morphological, behavioral, and physiological endpoints within an animal, (3) linking individual to population-level responses, and (4) testing genetic/genomic differences across an urban environment as evidence for selective pressures.

Juvenile social experience generates differences in behavioral variation but not averages

Developmental plasticity is known to influence the mean behavioral phenotype of a population. Yet, studies on how developmental plasticity shapes patterns of variation within populations are comparatively rare and often focus on a subset of developmental cues (e.g., nutrition). One potentially important but understudied developmental experience is social experience, as it is explicitly hypothesized to increase variation among individuals as a way to promote "social niches." To test this, we exposed juvenile black widow spiders (Latrodectus hesperus) to the silk of conspecifics by transplanting them onto conspecific webs for 48 h once a week until adulthood. We also utilized an untouched control group as well as a disturbed group. This latter group was removed from their web at the same time points as the social treatment, but was immediately placed back on their own web. After repeatedly measuring adult behavior and web structure, we found that social rearing drove higher or significant levels of repeatability relative to the other treatments. Repeatability in the social treatment also decreased in some traits, paralleling the decreases observed in the disturbed treatments. Thus, repeated juvenile disturbance may decrease among-individual differences in adult spiders. Yet, social rearing appeared to override the effect of disturbance in some traits, suggesting a prioritization effect. The resulting individual differences were maintained over at least one-third of the adult lifespan and thus appear to represent stable, canalized developmental effects and not temporal state differences. These results provide proximate insight into how a broader range of developmental experiences shape trait variation.

Muscle-induced loading as an important source of variation in craniofacial skeletal shape

The shape of the craniofacial skeleton is constantly changing through ontogeny and reflects a balance between developmental patterning and mechanical-load-induced remodeling. Muscles are a major contributor to producing the mechanical environment that is crucial for "normal" skull development. Here, we use an F₅ hybrid population of Lake Malawi cichlids to characterize the strength and types of associations between craniofacial bones and muscles. We focus on four bones/bone complexes, with different developmental origins, alongside four muscles with distinct functions. We used micro-computed tomography to extract 3D information on bones and muscles. 3D geometric morphometrics and volumetric measurements were used to characterize bone and muscle shape, respectively. Linear regressions were performed to test for associations between bone shape and muscle volume. We identified three types of associations between muscles and bones: weak, strong direct (i.e., muscles insert directly onto bone), and strong indirect (i.e., bone is influenced by muscles without a direct connection). In addition, we show that although the shape of some bones is relatively robust to muscle-induced mechanical stimulus, others appear to be highly sensitive to muscular input. Our results imply that the roles for muscular input on skeletal shape extend beyond specific points of origin or insertion and hold significant potential to influence broader patterns of craniofacial geometry. Thus, changes in the loading environment, either as a normal course of ontogeny or if an organism is exposed to a novel environment, may have pronounced effects on skeletal shape via near and far-ranging effects of muscular loading.

Where does diversity come from? Linking geographical patterns of morphological, genetic, and environmental variation in wall lizards

BACKGROUND

Understanding how phenotypic variation scales from individuals, through populations, up to species, and how it relates to genetic and environmental factors, is essential for deciphering the evolutionary mechanisms that drive biodiversity. We used two species of Podarcis wall lizards to test whether phenotypic diversity within and divergence across populations follow concordant patterns, and to examine how phenotypic variation responds to genetic and environmental variability across different hierarchical levels of biological organization, in an explicit geographic framework.

RESULTS

We found a general concordance of phenotypic variation across hierarchical levels (i.e. individuals and populations). However, we also found that within-population diversity does not exhibit a coherent geographic structure for most traits, while among-population divergence does, suggesting that different mechanisms may underlie the generation of diversity at these two levels. Furthermore, the association of phenotypic variation with genetic and environmental factors varied extensively between hierarchical levels and across traits, hampering the identification of simple rules to explain what yields diversity.

CONCLUSIONS

Our results in some cases comply with general ecological and evolutionary predictions, but in others they are difficult to explain in the geographic framework used, suggesting that habitat characteristics and other regulatory mechanisms may have a more substantial contribution in shaping phenotypic diversity.

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.

Analysing morphological variation of appendages and labrum in 10 species of Heterocypris Claus, 1893 (Podocopida : Cyprididae) with additional description of Heterocypris exigua

Natural selection can lead to speciation, but its effects depend on amount of morphological variation within populations. In Ostracoda, the appendages enclosed within the calcified carapace are less influenced by environmental conditions in comparison with carapace valves. Here, we explored morphological variation of appendages and labrum in cosmopolitan genus Heterocypris, investigating the species occurring in Europe. With the aim to overcome taxonomical problems, both traditional and geometric morphometrics were applied. The relative lengths were analysed by one-way ANOVA and post-hoc tests to determine interspecific variation and sexual dimorphism. Significant traits were further analysed by Canonical Variate Analysis separately for males and females as segment width/length ratios and relative length of claws and setae ratios. Landmark-based geometric morphometrics were applied to analyse shape variation of the labrum and traits of the maxillular palp. From this, a new perspective on morphological variation within this genus emerged. Relative measurements of the second segment of the maxillular palp, and relative lengths of the natatory setae were found to be the characters that contribute most to discrimination among Heterocypris species. The most important discriminant traits in females are the relative lengths of the G2 and G1 claws. Heterocypris exigua is additionally described, as all findings demonstrated that this is the most distinctive species, the most similar to H. gevgelica. In conclusion, we found that in comparison with maxillular palp and other analysed characters, labrum traits showed as more reliable and deserve more attention in descriptions of new species.

The ontogeny of developmental buffering in lizard head shape

Canalization and developmental stability (DS) are important organismal properties involved in determining the level of phenotypic variation. Ontogenetic patterns of phenotypic variance components can shed light on the mechanistic basis of developmental buffering (DB). Here, we analyze how individual FA and among-individual variation in head shape change in ontogenetic series of three lizard species raised in laboratory. The degree of asymmetry increased slightly with size, suggesting that developmental mechanisms hypothesized to correct for deviations either do not exist, or that their efficiency is truncated with increasing size. Alternatively, they may need the disturbance as a trigger. The relationship between asymmetry and age was complex, with asymmetry being stable across the age range in two species but increased with age in the third. Lack of congruence in ontogenetic patterns of asymmetry might be due to intrinsic differences in buffering mechanisms or a result of species-specific growth patterns. Head shape was shown to be equally canalized across both size and age range in all species, probably as a result of a balance between the buffering mechanisms and mechanisms generating variance. The patterns of symmetric and asymmetric head shape variation were highly correlated across species meaning that DS and canalization may rely on similar mechanisms.

The causes and ecological correlates of head scale asymmetry and fragmentation in a tropical snake

The challenge of identifying the proximate causes and ecological consequences of phenotypic variation can be facilitated by studying traits that are usually but not always bilaterally symmetrical; deviations from symmetry likely reflect disrupted embryogenesis. Based on a 19-year mark-recapture study of >1300 slatey-grey snakes (Stegonotus cucullatus) in tropical Australia, and incubation of >700 eggs, we document developmental and ecological correlates of two morphological traits: asymmetry and fragmentation of head scales. Asymmetry was directional (more scales on the left side) and was higher in individuals with lower heterozygosity, but was not heritable. In contrast, fragmentation was heritable and was higher in females than males. Both scale asymmetry and fragmentation were increased by rapid embryogenesis but were not affected by hydric conditions during incubation. Snakes with asymmetry and fragmentation exhibited slightly lower survival and increased (sex-specific) movements, and females with more scale fragmentation produced smaller eggs. Counterintuitively, snakes with more asymmetry had higher growth rates (possibly reflecting trade-offs with other traits), and snakes with more fragmentation had fewer parasites (possibly due to lower feeding rates). Our data paint an unusually detailed picture of the complex genetic and environmental factors that, by disrupting early embryonic development, generate variations in morphology that have detectable correlations with ecological performance.

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.

Development and evolution of segmentation assessed by geometric morphometrics: The centipede Strigamia maritima as a case study

Using the centipede model species Strigamia maritima as a subject of study, we illustrate the potential of geometric morphometrics for investigating the development and evolution of segmentation, with a specific focus on post-embryonic segmental patterning. We show how these techniques can contribute detailed descriptive data for comparative purposes, but also precious information on some features of the developmental system that are considered relevant for the evolvability of a segmented body architecture, such as developmental stability and canalization. Morphometric analyses allow to separately investigate several sources of phenotypic variation along a segmented body axis, like constitutive and random segment heteronomy, both within and among individuals. Specifically, in S. maritima, the segmental pattern of ventral sclerite shapes mirrors that of their bilateral fluctuating asymmetry and among-individual variation in associating the most anterior and most posterior segments in diverging from the central ones. Also, among segments, there seems to be a correlation between fluctuating asymmetry and shape variation among individuals, suggesting that canalization and developmental stability are somehow associated. Overall, these associations might stem from a joint influence of the segmental position on the two processes of developmental buffering.

Urban-Rural Differences in Eye, Bill, and Skull Allometry in House Finches (Haemorhous mexicanus)

Allometry, the proportional scaling of log trait size with log body size, evolves to optimize allocation to growth of separate structures and is a major constraint on the functional limits of animal traits. While there are many studies demonstrating the rigidity of allometry across traits and taxa, comparatively less work has been done on allometric variation across environments within species. Rapidly changing environments, such as cities, may be prime systems for studying the flexibility of allometry because they uniquely alter many environmental parameters (e.g., habitat, light, noise). We studied size variation, allometry, and allometric dispersion of craniofacial traits in both sexes of urban and rural house finches (Haemorhous mexicanus) because many traits in the head are ecologically critical to the survival and acclimation of birds to their environment (e.g., brain: response to cognitive challenges; bill: foraging modes). We found that urban finches had shorter eye axial lengths and skull widths, but longer (but not wider or deeper) bills, than rural finches. Also, allometric slopes of eye, skull, and bill traits differed based on sex and environment. In the rural environment, females had a far steeper allometric slope for eye axial length than males, but such slopes were similar between males and females in the city. Skull allometry was similar for males and females in both environments, but urban birds had a shallower slope for skull length (but not width) than rural birds. Other traits only differed by sex (males had a steeper slope for bill width), and one trait did not differ based on either sex or environment (bill depth). The dispersion of points around the allometric line did not differ by sex or environment for any craniofacial variable. Due to the significant but low genetic divergence between urban and rural finch populations, allometric differences are probably largely driven by plastic forces. We suggest that differences in diet and cognitive demand of urban environments may drive these allometric patterns. Overall, these results indicate that allometry may shift due to rapid environmental change and differentially so between the sexes.

Size, shape, and form: concepts of allometry in geometric morphometrics

Allometry refers to the size-related changes of morphological traits and remains an essential concept for the study of evolution and development. This review is the first systematic comparison of allometric methods in the context of geometric morphometrics that considers the structure of morphological spaces and their implications for characterizing allometry and performing size correction. The distinction of two main schools of thought is useful for understanding the differences and relationships between alternative methods for studying allometry. The Gould-Mosimann school defines allometry as the covariation of shape with size. This concept of allometry is implemented in geometric morphometrics through the multivariate regression of shape variables on a measure of size. In the Huxley-Jolicoeur school, allometry is the covariation among morphological features that all contain size information. In this framework, allometric trajectories are characterized by the first principal component, which is a line of best fit to the data points. In geometric morphometrics, this concept is implemented in analyses using either Procrustes form space or conformation space (the latter also known as size-and-shape space). Whereas these spaces differ substantially in their global structure, there are also close connections in their localized geometry. For the model of small isotropic variation of landmark positions, they are equivalent up to scaling. The methods differ in their emphasis and thus provide investigators with flexible tools to address specific questions concerning evolution and development, but all frameworks are logically compatible with each other and therefore unlikely to yield contradictory results.

Size, shape, and form: concepts of allometry in geometric morphometrics

Allometry refers to the size-related changes of morphological traits and remains an essential concept for the study of evolution and development. This review is the first systematic comparison of allometric methods in the context of geometric morphometrics that considers the structure of morphological spaces and their implications for characterizing allometry and performing size correction. The distinction of two main schools of thought is useful for understanding the differences and relationships between alternative methods for studying allometry. The Gould–Mosimann school defines allometry as the covariation of shape with size. This concept of allometry is implemented in geometric morphometrics through the multivariate regression of shape variables on a measure of size. In the Huxley–Jolicoeur school, allometry is the covariation among morphological features that all contain size information. In this framework, allometric trajectories are characterized by the first principal component, which is a line of best fit to the data points. In geometric morphometrics, this concept is implemented in analyses using either Procrustes form space or conformation space (the latter also known as size-and-shape space). Whereas these spaces differ substantially in their global structure, there are also close connections in their localized geometry. For the model of small isotropic variation of landmark positions, they are equivalent up to scaling. The methods differ in their emphasis and thus provide investigators with flexible tools to address specific questions concerning evolution and development, but all frameworks are logically compatible with each other and therefore unlikely to yield contradictory results.

Genetic and phenotypic component in head shape of common wall lizard Podarcis muralis

Head shape in lizards correlates with a wide range of environmental pressures, supporting the hypothesis that patterns of phenotypic change represent adaptive responses to selective processes. However, natural selection promotes evolutionary adaptation only if the trait under selection has enough heritable variation. In this study we used geometric morphometrics and quantitative genetics to assess the heritability patterns of the head shape and size of common wall lizards (Podarcis muralis). Genetic and phenotypic components were estimated using animal models, which showed that more than half of the variation in head morphology is inheritable. Furthermore, at least five independent patterns of genetically determined phenotypic change were detected. These outcomes confirm that morphological differentiation in common wall lizards may reliably be regarded as the result of adaptive processes driven by natural selection.