A GENERAL FRAMEWORK FOR THE ANALYSIS OF PHENOTYPIC TRAJECTORIES IN EVOLUTIONARY STUDIES

Does mercury affect morphology, developmental stability and canalization of the skull in the Common wall lizard (Podarcis muralis)?

Environmental pollution by metals and metalloids can have a detrimental effect on the fitness and development of organisms. Studies on the influence of metals and metalloids as environmental stressors on developmentally and functionally complex morphological structures of reptiles are important, as this group of vertebrates is highly threatened and is an important component of food webs. To assess the effects of chronic mercury exposure on cranium morphology and post-natal development in a model species of lizards, we analysed the concentration of this metal in liver tissues in the population of the Common wall lizard (Podarcis muralis) from the mercury mine tailings and the control population, and possible differences in skull size, shape, developmental stability and canalization between the two groups. Patterns of variation and asymmetry of the cranium shape were analysed using geometric morphometrics. The mercury concentration was significantly higher in the population from the polluted locality, but had very little to no effect on the cranial morphology. Juveniles and females from both sites had the same size and shape of the dorsal and ventral cranium, while males showed small differences in ventral cranium shape, reflected in slightly longer maxillae and wider crania at the polluted site. The pattern of sexual dimorphism remained constant at both localities. Both static and ontogenetic allometry were significant in both groups, and allometric trajectories did not differ between the two sites. The differences in fluctuating asymmetry (FA) between localities were not statistically significant for the dorsal and ventral cranium. The results do not support the idea that FA can be used as an early indicator mercury exposure at the population level.

Ontogenetic trajectories and sexual dimorphism of a neotropical ground snake genus Erythrolamprus (Serpentes: Dipsadidae)

Allometry and sexual dimorphism are potential sources of intraspecific morphological variation. Advances in the methodological framework of geometric morphometrics allow in-depth analysis of these issues, enabling the descriptions and comparisons between groups of multivariate phenotypic attributes and the visualization of allometric trajectories. In the present study, we evaluated the presence of secondary sexual dimorphism in Erythrolamprus miliaris merremi, a semi-aquatic snake occurring in the Atlantic Forest of southeastern Brazil. We assessed linear body measurements and digitized landmarks in the dorsal view of the head of 107 specimens grouped into four categories combining ontogenetic stages and sex. We performed linear models and principal component analyses to estimate and visualize head shape variation and the allometric trajectories of these categories. Adults from both sexes and immature females showed significant static allometry, while immature males showed isometric growth. Although we recovered no sexual dimorphism for head shape, ontogenetic trajectories of both sexes are distinct, with females showing head shape variation throughout ontogenetic development, while males exhibited more conspicuous changes only after sexual maturity. Comparisons with literature data suggest a high degree of variation in ontogenetic allometry of snakes, depending on the phylogenetic group or structure (head/skull) analyzed, while the results for static allometry are more similar between the studied taxa. Only the collection of data on distinct groups of snakes will provide more clues as to whether there is any allometric pattern to these structures, and which factors (evolutionary, ecological or both) are prevalent over it, especially for ontogenetic allometry.

Divergence time and environmental similarity predict the strength of morphological convergence in stick and leaf insects

Independent evolution of similar traits in lineages inhabiting similar environments (convergent or repeated evolution) is often taken as evidence for adaptation by natural selection, and used to illustrate the predictability of evolution. Yet convergence is rarely perfect for two reasons. First, environments may not be as similar as they appear. Second, responses to selection are contingent upon available genetic variation and independent lineages may differ in the alleles, genetic backgrounds, and even the developmental mechanisms responsible for the phenotypes in question. Both impediments to convergence are predicted to increase as the length of time separating two lineages increases, making it difficult to discern their relative importance. We quantified environmental similarity and the extent of convergence to show how habitat and divergence time each contribute to observed patterns of morphological evolution in 212 species of stick and leaf insects (order Phasmatodea). Dozens of phasmid lineages independently colonized similar habitats, repeatedly evolving in parallel directions on a 23-trait morphospace, though the magnitude and direction of these shifts varied. Lineages converging toward more similar environments ended up closer on the morphospace, as did closely related lineages, and closely related lineages followed more parallel evolutionary trajectories to arrive there than more distantly related ones. Remarkably, after accounting for habitat similarity, we show that divergence time reduced the extent of convergence at a constant rate across more than 100 My of separation, suggesting even the magnitude of contingency can be predictable, given sufficient spans of time.

Combined toxicity of perfluoroalkyl substances and microplastics on the sentinel species Daphnia magna: Implications for freshwater ecosystems

Persistent chemicals from industrial processes, particularly perfluoroalkyl substances (PFAS), have become pervasive in the environment due to their persistence, long half-lives, and bioaccumulative properties. Used globally for their thermal resistance and repellence to water and oil, PFAS have led to widespread environmental contamination. These compounds pose significant health risks with exposure through food, water, and dermal contact. Aquatic wildlife is particularly vulnerable as water bodies act as major transport and transformation mediums for PFAS. Their co-occurrence with microplastics may intensify the impact on aquatic species by influencing PFAS sorption and transport. Despite progress in understanding the occurrence and fate of PFAS and microplastics in aquatic ecosystems, the toxicity of PFAS mixtures and their co-occurrence with other high-concern compounds remains poorly understood, especially over organisms' life cycles. Our study investigates the chronic toxicity of PFAS and microplastics on the sentinel species Daphnia, a species central to aquatic foodwebs and an ecotoxicology model. We examined the effects of perfluorooctane sulfonate (PFOS), perfluorooctanoic acid (PFOA), and polyethylene terephthalate microplastics (PET) both individually and in mixtures on Daphnia ecological endpoints. Unlike conventional studies, we used two Daphnia genotypes with distinct histories of chemical exposure. This approach revealed that PFAS and microplastics cause developmental failures, delayed sexual maturity and reduced somatic growth, with historical exposure to environmental pollution reducing tolerance to these persistent chemicals due to cumulative fitness costs. We also observed that the combined effect of the persistent chemicals analysed was 59% additive and 41% synergistic, whereas no antagonistic interactions were observed. The genotype-specific responses observed highlight the complex interplay between genetic background and pollutant exposure, emphasizing the importance of incorporating multiple genotypes in environmental risk assessments to more accurately predict the ecological impact of chemical pollutants.

Interactions between phenanthrene exposure and historical chemical stress: Implications for fitness and ecological resilience of the sentinel species Daphnia magna

Polycyclic aromatic hydrocarbons (PAHs) arise from incomplete combustion of oil, coal, and gasoline, with lipophilic properties facilitating their widespread distribution and persistence. Due to their biochemical attributes, PAHs can accumulate in animal tissues, potentially causing mutagenic and carcinogenic effects. Since the industrial revolution, PAH concentrations in the environment have risen, with lakes showing levels from 0.159 to 33,090 μg/kg sediment. Despite acute toxicity studies showing adverse effects on freshwater organisms, the long-term impacts and synergistic interactions with other pollutants remain largely unexplored. This study investigates the impact of phenanthrene (PHE), a prominent PAH found in aquatic environments, on Daphnia magna, a species of significant ecological importance in freshwater ecosystems globally, being both a sentinel species for chemical pollution and a keystone organism in freshwater aquatic ecosystems. Leveraging the dormancy of D. magna, which spans decades or even centuries, we exposed strains with diverse histories of chemical contaminant exposure to environmentally relevant concentrations of PHE. Initially, acute exposure experiments were conducted in accordance with OECD guidelines across 16 Daphnia strains, revealing substantial variation in acute toxic responses, with strains naïve to chemical pollutants showing the lowest toxicity. Utilizing the median effect concentration EC10 derived from acute exposures, we assessed the impacts of chronic PHE exposure on life history traits and ecological endpoints of the 16 strains. To elucidate how historical exposure to other environmental stressors may modulate the toxicity of PHE, temporal populations of D. magna resurrected from a lake with a well-documented century-spanning history of environmental impact were utilized. Our findings demonstrate that PHE exposure induces developmental failure, delays sexual maturation, and reduces adult size in Daphnia. Populations of Daphnia historically exposed to chemical stress exhibited significantly greater fitness impacts compared to naïve populations. This study provides crucial insights into the augmented effects of PAHs interacting with other environmental stressors.

Replicated Functional Evolution in Cichlid Adaptive Radiations

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

The influence of the land-to-sea macroevolutionary transition on vertebral column disparification in Pinnipedia

The repeated returns of vertebrates to the marine ecosystems since the Triassic serve as an evolutionary model to understand macroevolutionary change. Here we investigate the effects of the land-to-sea transition on disparity and constraint of the vertebral column in aquatic carnivorans (Carnivora; Pinnipedia) to assess how their functional diversity and evolutionary innovations influenced major radiations of crown pinnipeds. We use three-dimensional geometric morphometrics and multivariate analysis for high-dimensional data under a phylogenetic framework to quantify vertebral size and shape in living and extinct pinnipeds. Our analysis demonstrates an important shift in vertebral column evolution by 10-12 million years ago, from an unconstrained to a constrained evolutionary scenario, a point of time that coincides with the major radiation of crown pinnipeds. Moreover, we also demonstrate that the axial skeleton of phocids and otariids followed a different path of morphological evolution that was probably driven by their specialized locomotor strategies. Despite this, we found a significant effect of habitat preference (coastal versus pelagic) on vertebral morphology of crown taxa regardless of the family they belong. In summary, our analysis provides insights into how the land-to-sea transition influenced the complex evolutionary history of pinniped vertebral morphology.

Functionally mediated cranial allometry evidenced in a genus of rock-wallabies

In assessments of skeletal variation, allometry (disproportionate change of shape with size) is often corrected to examine size-independent variation for hypotheses relating to function. However, size-related trade-offs in functional demands may themselves be an underestimated driver of mammalian cranial diversity. Here, we use geometric morphometrics alongside dental measurements to assess craniodental allometry in the rock-wallaby genus Petrogale (all 17 species, 370 individuals). We identified functional aspects of evolutionary allometry that can be both extensions of, and correlated negatively with, static or ontogenetic allometric patterns. Regarding constraints, larger species tended to have relatively smaller braincases and more posterior orbits, the former of which might represent a constraint on jaw muscle anatomy. However, they also tended to have more anterior dentition and smaller posterior zygomatic arches, both of which support the hypothesis of relaxed bite force demands and accommodation of different selective pressures that favour facial elongation. By contrast, two dwarf species had stouter crania with divergent dental adaptations that together suggest increased relative bite force capacity. This likely allows them to feed on forage that is mechanically similar to that consumed by larger relatives. Our results highlight a need for nuanced considerations of allometric patterns in future research of mammalian cranial diversity.

Multiple morphophysiological responses of a tropical frog to urbanization conform to the pace-of-life syndrome

The Pace-of-Life syndrome proposes that behavioural, physiological and immune characteristics vary along a slow-fast gradient. Urbanization poses several physiological challenges to organisms. However, little is known about how the health status of frogs is affected by urbanization in the Tropics, which have a faster and more recent urbanization than the northern hemisphere. Here, we analysed a suite of physiological variables that reflect whole organism health, reproduction, metabolic and circulatory physiology and leukocyte responses in Leptodactylus podicipinus. Specifically, we tested how leukocyte profile, erythrocyte morphometrics and germ cell density, as well as somatic indices and erythrocyte nuclear abnormalities differ throughout the adult life span between urban and rural populations. We used Phenotypic Trajectory Analysis to test the effect of age and site on each of the multivariate data sets; and a Generalised Linear Model to test the effect of site and age on nuclear abnormalities. Somatic indices, erythrocyte nuclear abnormalities, erythrocyte morphometrics and leukocyte profile differed between populations, but less so for germ cell density. We found a large effect of site on nuclear abnormalities, with urban frogs having twice as many abnormalities as rural frogs. Our results suggest that urban frogs have a faster pace of life, but the response of phenotypic compartments is not fully concerted.

The effects of host ecology and phylogeny on gut microbiota (non)parallelism across birds and mammals

IMPORTANCE

What are the roles of determinism and contingency in evolution? The paleontologist and evolutionary biologist Stephen J. Gould raised this question in his famous thought experiment of "replaying life's tape." Settings where independent lineages have repeatedly adapted to similar ecological niches (i.e., parallel evolution) are well suited to address this question. Here, we quantified whether repeated ecological shifts across 53 mammalian and 50 avian host species are associated with parallel gut microbiota changes. Our results indicate that parallel shifts in host diet are associated with greater gut microbiota parallelism (i.e., more deterministic). While further research will be necessary to obtain a comprehensive picture of the circumstances under which deterministic gut microbiota changes might be expected, our study can be instrumental in motivating the use of more quantitative methods in microbiota research. This, in turn, can help us better understand microbiota dynamics during adaptive evolution of their hosts.

Role of phenotypic and transcriptomic plasticity in alpine adaptation of Arabidopsis arenosa

Plasticity is an important component of the response of organism to environmental changes, but whether plasticity facilitates adaptation is still largely debated. Using transcriptomic and phenotypic data, we explored the evolution of ancestral plasticity during alpine colonization in Arabidopsis arenosa. We leveraged naturally replicated adaptation in four distinct mountain regions in Central Europe. We sampled seeds from ancestral foothill and independently formed alpine populations in each region and raised them in growth chambers under conditions approximating their natural environments. We gathered RNA-seq and genetic data of 48 and 63 plants and scored vegetative and flowering traits in 203 and 272 plants respectively. Then, we compared gene expression and trait values over two treatments differing in temperature and irradiance and elevations of origin and quantified the extent of ancestral and derived plasticity. At the transcriptomic level, initial plastic changes tended to be more reinforced than reversed in adapted alpine populations. Genes showing reinforcement were involved in the stress response, developmental processes and morphogenesis and those undergoing reversion were related to the stress response (light and biotic stress). At the phenotypic level, initial plastic changes in all but one trait were also reinforced supporting a facilitating role of phenotypic plasticity during colonization of an alpine environment. Our results contrasted with previous studies that showed generally higher reversion than reinforcement and supported the idea that ancestral plasticity tends to be reinforced in the context of alpine adaptation. However, plasticity may also be the source of potential maladaptation, especially at the transcriptomic level.

Eco-evolutionary dynamics of experimental Pseudomonas aeruginosa populations under oxidative stress

Within-host environments are likely to present a challenging and stressful environment for opportunistic pathogenic bacteria colonizing from the external environment. How populations of pathogenic bacteria respond to such environmental challenges and how this varies between strains is not well understood. Oxidative stress is one of the defences adopted by the human immune system to confront invading bacteria. In this study, we show that strains of the opportunistic pathogenic bacterium Pseudomonas aeruginosa vary in their eco-evolutionary responses to hydrogen peroxide stress. By quantifying their 24 h growth kinetics across hydrogen peroxide gradients we show that a transmissible epidemic strain isolated from a chronic airway infection of a cystic fibrosis patient, LESB58, is much more susceptible to hydrogen peroxide than either of the reference strains, PA14 or PAO1, with PAO1 showing the lowest susceptibility. Using a 12 day serial passaging experiment combined with a mathematical model, we then show that short-term susceptibility controls the longer-term survival of populations exposed to subinhibitory levels of hydrogen peroxide, but that phenotypic evolutionary responses can delay population extinction. Our model further suggests that hydrogen peroxide driven extinctions are more likely with higher rates of population turnover. Together, these findings suggest that hydrogen peroxide is likely to be an effective defence in host niches where there is high population turnover, which may explain the counter-intuitively high susceptibility of a strain isolated from chronic lung infection, where such ecological dynamics may be slower.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

Latitude-specific urbanization effects on life history traits in the damselfly Ischnura elegans

Many species are currently adapting to cities at different latitudes. Adaptation to urbanization may require eco-evolutionary changes in response to temperature and invasive species that may differ between latitudes. Here, we studied single and combined effects of increased temperatures and an invasive alien predator on the phenotypic response of replicated urban and rural populations of the damselfly Ischnura elegans and contrasted these between central and high latitudes. Adult females were collected in rural and urban ponds at central and high latitudes. Their larvae were exposed to temperature treatments (current [20°C], mild warming [24°C], and heat wave [28°C; for high latitude only]) crossed with the presence or absence of chemical cues released by the spiny-cheek crayfish (Faxonius limosus), only present at the central latitude. We measured treatment effects on larval development time, mass, and growth rate. Urbanization type affected all life history traits, yet these responses were often dependent on latitude, temperature, and sex. Mild warming decreased mass in rural and increased growth rate in urban populations. The effects of urbanization type on mass were latitude-dependent, with central-latitude populations having a greater phenotypic difference. Urbanization type effects were sex-specific with urban males being lighter and having a lower growth rate than rural males. At the current temperature and mild warming, the predator cue reduced the growth rate, and this independently of urbanization type and latitude of origin. This pattern was reversed during a heat wave in high-latitude damselflies. Our results highlight the context-dependency of evolutionary and plastic responses to urbanization, and caution for generalizing how populations respond to cities based on populations at a single latitude.

The Difficulty of Predicting Evolutionary Change in Response to Novel Ecological Interactions: A Field Experiment with Anolis Lizards

AbstractDetermining whether and how evolution is predictable is an important goal, particularly as anthropogenic disturbances lead to novel species interactions that could modify selective pressures. Here, we use a multigeneration field experiment with brown anole lizards (Anolis sagrei) to test hypotheses about the predictability of evolution. We manipulated the presence/absence of predators and competitors of A. sagrei across 16 islands in the Bahamas that had preexisting brown anole populations. Before the experiment and again after roughly five generations, we measured traits related to locomotor performance and habitat use by brown anoles and used double-digest restriction enzyme-associated DNA sequencing to estimate genome-wide changes in allele frequencies. Although previous work showed that predators and competitors had characteristic effects on brown anole behavior, diet, and population sizes, we found that evolutionary change at both phenotypic and genomic levels was difficult to forecast. Phenotypic changes were contingent on sex and habitat use, whereas genetic change was unpredictable and not measurably correlated with phenotypic changes, experimental treatments, or other environmental factors. Our work shows how differences in ecological context can alter evolutionary outcomes over short timescales and underscores the difficulty of forecasting evolutionary responses to multispecies interactions in natural conditions, even in a well-studied system with ample supporting ecological information.

Species-Specific Spatial Patterns of Variation in Sexual Dimorphism by Two Lizards Settled in the Same Geographic Context

The evolution of sexual dimorphism (SD) results from intricate interactions between sexual and natural selections. Sexually selected traits are expected to depend on individual condition, while natural selected traits should not be. Islands offer an ideal context to test how these drivers interact with one another, as the size is a reliable proxy for resource availability. Here, we analysed SD in body size (snout-vent length) and head shape (assessed by geometric morphometric) in two species of lizards (Podarcis muralis and P. siculus) inhabiting the Tuscan archipelago (Central Italy). We found a strong SD variation among islands in both species. Furthermore, in P. muralis emerged some significant correlations between SD and island size, supporting the occurrence of possible effects of individual condition on SD. By contrast, SD in P. siculus followed opposite trajectories than in P. muralis, suggesting that in this species, natural selection could play a major role as a driver of SD. Our findings show that natural and sexual selection can interact in complex ways, and the responses are species-specific. Therefore, spatial patterns of variation in SD may strongly differ among species, even when they settle in the same geographic contest.

The effect of captivity on craniomandibular and calcaneal ontogenetic trajectories in wild boar

Deciphering the plastic (i.e., nonheritable) changes induced by human control over wild animals in the archeological record is challenging. Previous studies detected morphological markers associated with captivity in the cranium, mandible, and calcaneus of adult wild boar (Sus scrofa) but the developmental trajectories leading up to these changes during ontogeny remain unknown. To assess the impact of growth in a captive environment on morphological structures during postnatal ontogeny, we used an experimental approach focusing on the same three structures and taxon. We investigated the form and size differences of captive-reared and wild-caught wild boar during growth using three-dimensional landmark-based geometric morphometrics. Our results provide evidence of an influence of captivity on the morphology of craniomandibular structures, as wild specimens are smaller than captive individuals at similar ages. The food resources inherent to anthropogenic environments may explain some of the observed differences between captive-reared and wild specimens. The calcaneus presents a different contrasted pattern of plasticity as captive and wild individuals differ in terms of form but not in terms of size. The physically more constrained nature of the calcaneus and the direct influence of mobility reduction on this bone may explain these discrepancies. These results provide new methodological perspectives for bioarchaeological approaches as they imply that the plastic mark of captivity can be observed in juvenile specimens in the same way it has been previously described in adults.

Quantifying (non)parallelism of gut microbial community change using multivariate vector analysis

Parallel evolution of phenotypic traits is regarded as strong evidence for natural selection and has been studied extensively in a variety of taxa. However, we have limited knowledge of whether parallel evolution of host organisms is accompanied by parallel changes of their associated microbial communities (i.e., microbiotas), which are crucial for their hosts' ecology and evolution. Determining the extent of microbiota parallelism in nature can improve our ability to identify the factors that are associated with (putatively adaptive) shifts in microbial communities. While it has been emphasized that (non)parallel evolution is better considered as a quantitative continuum rather than a binary phenomenon, quantitative approaches have rarely been used to study microbiota parallelism. We advocate using multivariate vector analysis (i.e., phenotypic change vector analysis) to quantify direction and magnitude of microbiota changes and discuss the applicability of this approach for studying parallelism, and we compiled an R package for multivariate vector analysis of microbial communities ('multivarvector'). We exemplify its use by reanalyzing gut microbiota data from multiple fish species that exhibit parallel shifts in trophic ecology. We found that multivariate vector analysis results were largely consistent with other statistical methods, parallelism estimates were not affected by the taxonomic level at which the microbiota is studied, and parallelism might be stronger for gut microbiota function compared to taxonomic composition. This approach provides an analytical framework for quantitative comparisons across host lineages, thereby providing the potential to advance our capacity to predict microbiota changes. Hence, we emphasize that the development and application of quantitative measures, such as multivariate vector analysis, should be further explored in microbiota research in order to better understand the role of microbiota dynamics during their hosts' adaptive evolution, particularly in settings of parallel evolution.

Parallel evolution in human populations: A biocultural perspective

Parallel evolution-where different populations evolve similar traits in response to similar environments-has been a topic of growing interest to biologists and biological anthropologists for decades. Parallel evolution occurs in human populations thanks to myriad biological and cultural mechanisms that permit humans to survive and thrive in diverse environments worldwide. Because humans shape and are shaped by their environments, biocultural approaches that emphasize the interconnections between biology and culture are key to understanding parallel evolution in human populations as well as the nuances of human biological variation and adaptation. In this review, we discuss how biocultural theory has been and can be applied to studies of parallel evolution and adaptation more broadly. We illustrate this through four examples of parallel evolution in humans: malaria resistance, lactase persistence, cold tolerance, and high-altitude adaptation.

Functional constraints channel mandible shape ontogenies in rodents

In mammals, postnatal growth plays an essential role in the acquisition of the adult shape. During this period, the mandible undergoes many changing functional constraints, leading to spatialization of bone formation and remodelling to accommodate various dietary and behavioural changes. The interactions between the bone, muscles and teeth drive this developmental plasticity, which, in turn, could lead to convergences in the developmental processes constraining the directionality of ontogenies, their evolution and thus the adult shape variation. To test the importance of the interactions between tissues in shaping the ontogenetic trajectories, we compared the mandible shape at five postnatal stages on three rodents: the house mouse, the Mongolian gerbil and the golden hamster, using geometric morphometrics. After an early shape differentiation, by both longer gestation and allometric scaling in gerbils or early divergence of postnatal ontogeny in hamsters in comparison with the mouse, the ontogenetic trajectories appear more similar around weaning. The changes in muscle load associated with new food processing and new behaviours at weaning seem to impose similar physical constraints on the mandible, driving the convergences of the ontogeny at that stage despite an early anatomical differentiation. Nonetheless, mice present a rather different timing compared with gerbils or hamsters.

The terroir of the finch: How spatial and temporal variation shapes phenotypic traits in DARWIN'S finches

The term terroir is used in viticulture to emphasize how the biotic and abiotic characteristics of a local site influence grape physiology and thus the properties of wine. In ecology and evolution, such terroir (i.e., the effect of space or "site") is expected to play an important role in shaping phenotypic traits. Just how important is the pure spatial effect of terroir (e.g., differences between sites that persist across years) in comparison to temporal variation (e.g., differences between years that persist across sites), and the interaction between space and time (e.g., differences between sites change across years)? We answer this question by analyzing beak and body traits of 4388 medium ground finches (Geospiza fortis) collected across 10 years at three locations in Galápagos. Analyses of variance indicated that phenotypic variation was mostly explained by site for beak size (η 2 = 0.42) and body size (η 2 = 0.43), with a smaller contribution for beak shape (η 2 = 0.05) and body shape (η 2 = 0.12), but still higher compared to year and site-by-year effects. As such, the effect of terroir seems to be very strong in Darwin's finches, notwithstanding the oft-emphasized interannual variation. However, these results changed dramatically when we excluded data from Daphne Major, indicating that the strong effect of terroir was mostly driven by that particular population. These phenotypic results were largely paralleled in analyses of environmental variables (rainfall and vegetation indices) expected to shape terroir in this system. These findings affirm the evolutionary importance of terroir, while also revealing its dependence on other factors, such as geographical isolation.

An effect size for comparing the strength of morphological integration across studies

Understanding how and why phenotypic traits covary is a major interest in evolutionary biology. Biologists have long sought to characterize the extent of morphological integration in organisms, but comparing levels of integration for a set of traits across taxa has been hampered by the lack of a reliable summary measure and testing procedure. Here, we propose a standardized effect size for this purpose, calculated from the relative eigenvalue variance,V r e l $V_{rel}$ . First, we evaluate several eigenvalue dispersion indices under various conditions, and show that onlyV r e l $V_{rel}$ remains stable across samples size and the number of variables. We then demonstrate thatV r e l $V_{rel}$ accurately characterizes input patterns of covariation, so long as redundant dimensions are excluded from the calculations. However, we also show that the variance of the sampling distribution ofV r e l $V_{rel}$ depends on input levels of trait covariation, makingV r e l $V_{rel}$ unsuitable for direct comparisons. As a solution, we propose transformingV r e l $V_{rel}$ to a standardized effect size (Z-score) for representing the magnitude of integration for a set of traits. We also propose a two-sample test for comparing the strength of integration between taxa, and show that this test displays appropriate statistical properties. We provide software for implementing the procedure, and an empirical example illustrates its use.

Mosaic adaptive peak shifts underlie body shape diversification in pelagiarian fishes (Acanthomorpha: Percomorpha)

Extreme body elongation in fishes is a major evolutionary transformation that extends the boundaries of morphological diversity and alters aspects of function, behaviour and ecology. Prior studies have identified features of the cranial and axial skeleton that characterize elongate fishes, but a lack of detailed reconstructions of anatomical evolution has limited inferences about factors that underlie major shifts in body shape. In this study, we fitted multi-peak adaptive (Ornstein–Uhlenbeck) evolutionary models to species body shape and anatomical dimensions in Pelagiaria, a radiation of open-ocean fishes whose species span a continuum from deep bodied to highly elongate. We inferred an ancestral fusiform adaptive peak that is retained by several major pelagiarian lineages (e.g. Scombridae) and found robust support for multiple transitions to deep-bodied optima (in the families Stromateidae, Bramidae and Caristiidae) and elongate-bodied optima (within Trichiuroidei), including two instances of sequential shifts towards increasingly elongate optima that followed distinct paths of anatomical evolution. Within Trichiuridae, initial increases in head length and the number of vertebrae were followed by changes in head and vertebral shape. Within an elongate-bodied subclade of taxa traditionally identified as ‘gempylids’, changes in head and vertebral shape and in the number of precaudal vertebrae preceded an increase in the number of caudal vertebrae. Altogether, this mosaic of anatomical peak shifts suggests that body shape transformations were associated with differing selective demands and developmental changes.

The evolutionary trajectory of drosophilid walking

Neural circuits must both execute the behavioral repertoire of individuals and account for behavioral variation across species. Understanding how this variation emerges over evolutionary time requires large-scale phylogenetic comparisons of behavioral repertoires. Here, we describe the evolution of walking in fruit flies by capturing high-resolution, unconstrained movement from 13 species and 15 strains of drosophilids. We find that walking can be captured in a universal behavior space, the structure of which is evolutionarily conserved. However, the occurrence of and transitions between specific movements have evolved rapidly, resulting in repeated convergent evolution in the temporal structure of locomotion. Moreover, a meta-analysis demonstrates that many behaviors evolve more rapidly than other traits. Thus, the architecture and physiology of locomotor circuits can execute precise individual movements in one species and simultaneously support rapid evolutionary changes in the temporal ordering of these modular elements across clades.

Variation in the genomic basis of parallel phenotypic and ecological divergence in benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Sympatric adaptive phenotypic divergence should be underlain by genomic differentiation between sub-populations. When divergence drives similar patterns of phenotypic and ecological variation within species we expect evolution to draw on common allelic variation. We investigated divergence histories and genomic signatures of adaptive divergence between benthic and pelagic morphs of Icelandic Arctic charr. Divergence histories for each of four populations were reconstructed using coalescent modelling and 14,187 single nucleotide polymorphisms. Sympatric divergence with continuous gene flow was supported in two populations while allopatric divergence with secondary contact was supported in one population; we could not differentiate between demographic models in the fourth population. We detected parallel patterns of phenotypic divergence along benthic-pelagic evolutionary trajectories among populations. Patterns of genomic differentiation between benthic and pelagic morphs were characterized by outlier loci in many narrow peaks of differentiation throughout the genome, which may reflect the eroding effects of gene flow on nearby neutral loci. We then used genome-wide association analyses to relate both phenotypic (body shape and size) and ecological (carbon and nitrogen stable isotopes) variation to patterns of genomic differentiation. Many peaks of genomic differentiation were associated with phenotypic and ecological variation in the three highly divergent populations, suggesting a genomic basis for adaptive divergence. We detected little evidence for a parallel genomic basis of differentiation as most regions and outlier loci were not shared among populations. Our results show that adaptive divergence can have varied genomic consequences in populations with relatively recent common origins, similar divergence histories, and parallel phenotypic divergence.

Hybrid evolution repeats itself across environmental contexts in Texas sunflowers (Helianthus)

To what extent is evolution repeatable? Little is known about whether the evolution of hybrids is more (or less) repeatable than that of nonhybrids. We used field experimental evolution in annual sunflowers (Helianthus) in Texas to ask the extent to which hybrid evolution is repeatable across environments compared to nonhybrid controls. We created hybrids between Helianthus annuus (L.) and H. debilis (Nutt.) and grew plots of both hybrids and nonhybrid controls through eight generations at three sites in Texas. We collected seeds from each generation and grew each generation × treatment × home site combination at two final common gardens. We estimated the strength and direction of evolution in terms of fitness and 24 traits, tested for repeated versus nonrepeated evolution, and assessed overall phenotypic evolution across lineages and in relation to a locally adapted phenotype. Hybrids consistently evolved higher fitness over time, while controls did not, although trait evolution varied in strength across home sites. Repeated evolution was more evident in hybrids versus nonhybrid controls, and hybrid evolution was often in the direction of the locally adapted phenotype. Our findings have implications for both the nature of repeatability in evolution and the contribution of hybridization to evolution across environmental contexts.

Historical exposure to chemicals reduces tolerance to novel chemical stress in Daphnia (waterflea)

Until the last few decades, anthropogenic chemicals used in most production processes have not been comprehensively assessed for their risk and impact on wildlife and humans. They are transported globally and usually end up in the environment as unintentional pollutants, causing long‐term adverse effects. Modern toxicology practices typically use acute toxicity tests of unrealistic concentrations of chemicals to determine their safe use, missing pathological effects arising from long‐term exposures to environmentally relevant concentrations. Here, we study the transgenerational effect of environmentally relevant concentrations of five chemicals on the priority list of international regulatory frameworks on the keystone species Daphnia magna. We expose Daphnia genotypes resurrected from the sedimentary archive of a lake with a known history of chemical pollution to the five chemicals to understand how historical exposure to chemicals influences adaptive responses to novel chemical stress. We measure within‐ and transgenerational plasticity in fitness‐linked life history traits following exposure of “experienced” and “naive” genotypes to novel chemical stress. As the revived Daphnia originate from the same genetic pool sampled at different times in the past, we are able to quantify the long‐term evolutionary impact of chemical pollution by studying genome‐wide diversity and identifying functional pathways affected by historical chemical stress. Our results suggest that historical exposure to chemical stress causes reduced genome‐wide diversity, leading to lower cross‐generational tolerance to novel chemical stress. Lower tolerance is underpinned by reduced gene diversity at detoxification, catabolism and endocrine genes in experienced genotypes. We show that these genes sit within pathways that are conserved and potential chemical targets in other species, including humans.

Variation in vertebrae shape across small-bodied newts reveals functional and developmental constraints acting upon the trunk region

The salamander vertebral column is largely undifferentiated with a series of more or less uniform rib-bearing presacral vertebrae traditionally designated as the trunk region. We explored regionalization of the salamander trunk in seven species and two subspecies of the salamander genus Lissotriton by the combination of microcomputed tomography scanning and geometric morphometrics. The detailed information on trunk vertebral shape was subjected to a multidimensional cluster analysis and a phenotypic trajectory analysis. With these complementary approaches, we observed a clear morphological regionalization. Clustering analysis showed that the anterior trunk vertebrae (T1 and T2) have distinct morphologies that are shared by all taxa, whereas the subsequent, more posterior vertebrae show significant disparity between species. The phenotypic trajectory analysis revealed that all taxa share a common pattern and amount of shape change along the trunk region. Altogether, our results support the hypothesis of a conserved anterior-posterior developmental patterning which can be associated with different functional demands, reflecting (sub)species' and, possibly, regional ecological divergences within species.

Warmer temperatures limit the effects of antidepressant pollution on life-history traits

Pharmaceutical pollutants pose a threat to aquatic ecosystems worldwide. Yet, few studies have considered the interaction between pharmaceuticals and other chronic stressors contemporaneously, even though the environmental challenges confronting animals in the wild seldom, if ever, occur in isolation. Thermal stress is one such environmental challenge that may modify the threat of pharmaceutical pollutants. Accordingly, we investigated how fluoxetine (Prozac), a common psychotherapeutic and widespread pollutant, interacts with temperature to affect life-history traits in the water flea, Daphnia magna. We chronically exposed two genotypes of Daphnia to two ecological relevant concentrations of fluoxetine (30 ng l^-1 and 300 ng l^-1) and a concentration representing levels used in acute toxicity tests (3000 ng l^-1) and quantified the change in phenotypic trajectories at two temperatures (20°C and 25°C). Across multiple life-history traits, we found that fluoxetine exposure impacted the fecundity, body size and intrinsic growth rate of Daphnia in a non-monotonic manner at 20°C, and often in genotypic-specific ways. At 25°C, however, the life-history phenotypes of individuals converged under the widely varying levels of fluoxetine, irrespective of genotype. Our study underscores the importance of considering the complexity of interactions that can occur in the wild when assessing the effects of chemical pollutants on life-history traits.

Detecting (non)parallel evolution in multidimensional spaces: angles, correlations and eigenanalysis

Parallelism between evolutionary trajectories in a trait space is often seen as evidence for repeatability of phenotypic evolution, and angles between trajectories play a pivotal role in the analysis of parallelism. However, properties of angles in multidimensional spaces have not been widely appreciated by biologists. To remedy this situation, this study provides a brief overview on geometric and statistical aspects of angles in multidimensional spaces. Under the null hypothesis that trajectory vectors have no preferred directions (i.e. uniform distribution on hypersphere), the angle between two independent vectors is concentrated around the right angle, with a more pronounced peak in a higher-dimensional space. This probability distribution is closely related to t- and beta distributions, which can be used for testing the null hypothesis concerning a pair of trajectories. A recently proposed method with eigenanalysis of a vector correlation matrix can be connected to the test of no correlation or concentration of multiple vectors, for which simple test procedures are available in the statistical literature. Concentration of vectors can also be examined by tools of directional statistics such as the Rayleigh test. These frameworks provide biologists with baselines to make statistically justified inferences for (non)parallel evolution.

Parallel shape divergence between ecotypes of the limpet Nacella concinna along the Antarctic Peninsula: a new model species for parallel evolution?

Parallel phenotypic divergence is the independent differentiation between phenotypes of the same lineage or species occupying ecologically similar environments in different populations. We tested in the Antarctic limpet Nacella concinna the extent of parallel morphological divergence in littoral and sublittoral ecotypes throughout its distribution range. These ecotypes differ in morphological, behavioural and physiological characteristics. We studied the lateral and dorsal outlines of shells and the genetic variation of the mitochondrial gene Cytochrome Oxidase subunit I from both ecotypes in 17 sample sites along more than 2,000 km. The genetic data indicate that both ecotypes belong to a single evolutionary lineage. The magnitude and direction of phenotypic variation differ between ecotypes across sample sites; completely parallel ecotype-pairs (i.e., they diverge in the same magnitude and in the same direction) were detected in 84.85% of lateral and 65.15% in dorsal view comparisons. Besides, specific traits (relative shell height, position of shell apex, and elliptical/pear-shape outline variation) showed high parallelism. We observed weak morphological covariation between the two shape shell views, indicating that distinct evolutionary forces and environmental pressures could be acting on this limpet shell shape. Our results demonstrate there is a strong parallel morphological divergence pattern in N. concinna along its distribution, making this Antarctic species a suitable model for the study of different evolutionary forces shaping the shell evolution of this limpet.

Size and shape regional differentiation during the development of the spine in the nine-banded armadillo (Dasypus novemcinctus)

Xenarthrans (armadillos, anteaters, sloths, and their extinct relatives) are unique among mammals in displaying a distinctive specialization of the posterior trunk vertebrae-supernumerary vertebral xenarthrous articulations. This study seeks to understand how xenarthry develops through ontogeny and if it may be constrained to appear within pre-existing vertebral regions. Using three-dimensional geometric morphometrics on the neural arches of vertebrae, we explore phenotypic, allometric, and disparity patterns of the different axial morphotypes during the ontogeny of nine-banded armadillos. Shape-based regionalization analyses showed that the adult thoracolumbar column is divided into three regions according to the presence or absence of ribs and the presence or absence of xenarthrous articulations. A three-region division was retrieved in almost all specimens through development, although younger stages (e.g., fetuses, neonates) have more region boundary variability. In size-based regionalization analyses, thoracolumbar vertebrae are separated into two regions: a prediaphragmatic, prexenarthrous region, and a postdiaphragmatic xenarthrous region. We show that posterior thoracic vertebrae grow at a slower rate, while anterior thoracics and lumbars grow at a faster rate relatively, with rates decreasing anteroposteriorly in the former and increasing anteroposteriorly in the latter. We propose that different proportions between vertebrae and vertebral regions might result from differences in growth pattern and timing of ossification.

Speciation dynamics and extent of parallel evolution along a lake-stream environmental contrast in African cichlid fishes

Understanding the dynamics of speciation is a central topic in evolutionary biology. Here, we investigated how morphological and genomic differentiation accumulated along the speciation continuum in the African cichlid fish Astatotilapia burtoni. While morphological differentiation was continuously distributed across different lake-stream population pairs, we found that there were two categories with respect to genomic differentiation, suggesting a “gray zone” of speciation at ~0.1% net nucleotide divergence. Genomic differentiation was increased in the presence of divergent selection and drift compared to drift alone. The quantification of phenotypic and genetic parallelism in four cichlid species occurring along a lake-stream environmental contrast revealed parallel and antiparallel components in rapid adaptive divergence, and morphological convergence in species replicates inhabiting the same environments. Furthermore, we show that the extent of parallelism was higher when ancestral populations were more similar. Our study highlights the complementary roles of divergent selection and drift on speciation and parallel evolution.

Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio

The independent and repeated adaptation of populations to similar environments often results in the evolution of similar forms. This phenomenon creates a strong correlation between phenotype and environment and is referred to as parallel evolution. However, we are still largely unaware of the dynamics of parallel evolution, as well as the interplay between phenotype and genotype within natural systems. Here, we examined phenotypic and genotypic parallel evolution in multiple parapatric Dune‐Headland coastal ecotypes of an Australian wildflower, Senecio lautus. We observed a clear trait‐environment association in the system, with all replicate populations having evolved along the same phenotypic evolutionary trajectory. Similar phenotypes have arisen via mutational changes occurring in different genes, although many share the same biological functions. Our results shed light on how replicated adaptation manifests at the phenotypic and genotypic levels within populations, and highlight S. lautus as one of the most striking cases of phenotypic parallel evolution in nature.

Divergent processes drive parallel evolution in marine and freshwater fishes

Evolutionary comparisons between major environmental divides, such as between marine and freshwater systems, can reveal the fundamental processes governing diversification dynamics. Although processes may differ due to the different scales of their biogeographic barriers, freshwater and marine environments nevertheless offer similar opportunities for diversification in benthic, demersal, and pelagic habitats. Here, we compare the evolutionary patterns and processes shaping teleost diversity both in each of these three habitats and between marine and freshwater systems. Using specimens from the National Museum of Natural History, we developed a dataset of linear measurements capturing body shape in 2,266 freshwater and 3,344 marine teleost species. With a novel comparative approach, we contrast the primary axis of morphological diversification in each habitat with the major axis defined by phylogenetic signal. By comparing angles between these axes, we find that fish in corresponding habitats have more similar primary axes of morphological diversity than would be expected by chance, but that different historical processes underlie these parallel patterns in freshwater and marine environments. Marine diversification is more strongly aligned with phylogenetic signal and shows a trend toward lineages occupying separate regions of morphospace. In contrast, ecological signal appears to be a strong driver of diversification in freshwater lineages through repeated morphological evolution in densely packed regions of morphospace. In spite of these divergent histories, our findings reveal that habitat has driven convergent patterns of evolutionary diversification on a global scale.

Postnatal ontogeny of the femur in fossorial and semiaquatic water voles in the 3D-shape space

Water voles of the genus Arvicola constitute an excellent subject to investigate to which extent function affects postnatal developmental growth of limb structures in phylogenetically close species. We performed a comparative analysis of postweaning femur form changes between Arvicola sapidus (semiaquatic) and Arvicola scherman (fossorial) using three-dimensional landmark-based geometric morphometrics. In both species, we observed greater femur robustness in juvenile individuals than in adult ones, probably due to the accommodation of high loads on the bone during initial locomotor efforts. Significant interspecific differences were also found in the femur size and shape of adult specimens, as well as in the postnatal allometric and phenotypic trajectories. In terms of phenotypic variation, fossorial water voles show relatively wider third and lesser trochanters, and greater femur robustness than A. sapidus, characters associated to the digging activity. In contrast, A. sapidus displays a slight increase of the greater trochanter in comparison with A. scherman, which is seemingly an adaptive response for enhancing propulsion through the water. Results evidence that certain morphological traits and differences between A. sapidus and A. scherman in the allometric and phenotypic trajectories of the femur are associated with their different locomotor mode.

Doublesex mediates species-, sex-, environment- and trait-specific exaggeration of size and shape

Context-dependent trait exaggeration is a major contributor to phenotypic diversity. However, the genetic modifiers instructing development across multiple contexts remain largely unknown. We use the arthropod tibia, a hotspot for segmental differentiation, as a paradigm to assess the developmental mechanisms underlying the context-dependent structural exaggeration of size and shape through nutritional plasticity, sexual dimorphism and segmental differentiation. Using an RNAseq approach in the sexually dimorphic and male-polyphenic dung beetle Digitonthophagus gazella, we find that only a small portion (3.7%) of all transcripts covary positively in expression level with trait size across contexts. However, RNAi-mediated knockdown of the conserved sex-determination gene doublesex suggests that it functions as a context-dependent master mediator of trait exaggeration in D. gazella as well as the closely related dung beetle Onthophagus taurus. Taken together, our findings suggest (i) that the gene networks associated with trait exaggeration are highly dependent on the precise developmental context, (ii) that doublesex differentially shapes morphological exaggeration depending on developmental contexts and (iii) that this context-specificity of dsx-mediated trait exaggeration may diversify rapidly. This mechanism may contribute to the resolution of conflict arising from environment-dependent antagonistic selection among sexes and divergent developmental contexts in a wide range of animals.

Temperature and pathogen exposure act independently to drive host phenotypic trajectories

Natural populations are experiencing an increase in the occurrence of both thermal stress and disease outbreaks. How these two common stressors interact to determine host phenotypic shifts will be important for population persistence, yet a myriad of different traits and pathways are a target of both stressors, making generalizable predictions difficult to obtain. Here, using the host Daphnia magna and its bacterial pathogen Pasteuria ramosa, we tested how temperature and pathogen exposure interact to drive shifts in multivariate host phenotypes. We found that these two stressors acted mostly independently to shape host phenotypic trajectories, with temperature driving a faster pace of life by favouring early development and increased intrinsic population growth rates, while pathogen exposure impacted reproductive potential through reductions in lifetime fecundity. Studies focussed on extreme thermal stress are increasingly showing how pathogen exposure can severely hamper the thermal tolerance of a host. However, our results suggest that under milder thermal stress, and in terms of life-history traits, increases in temperature might not exacerbate the impact of pathogen exposure on host performance, and vice versa.

Connexin 43 contributes to phenotypic robustness of the mouse skull

BACKGROUND

We compared skull shape and variation among genetically modified mice that exhibit different levels of connexin43 (Cx43) channel function, to determine whether Cx43 contributes to craniofacial phenotypic robustness. Specifically, we used two heterozygous mutant mouse models (G60S/+ and I130T/+) that, when compared to their wildtype counterparts, have an ~80% and ~50% reduction in Cx43 function, respectively.

RESULTS

Both mutant strains showed significant differences in skull shape compared to wildtype littermates and while these differences were more severe in the G60S/+ mouse, shape differences were localized to similar regions of the skull in both mutants. However, increased skull shape variation was observed in G60S/+ mutants only. Additionally, covariation of skull structures was disrupted in the G60S/+ mutants only, indicating that while a 50% reduction in Cx43 function is sufficient to cause a shift in mean skull shape, the threshold for Cx43 function for disrupting craniofacial phenotypic robustness is lower.

CONCLUSIONS

Collectively, our results indicate Cx43 can contribute to phenotypic robustness of the skull through a nonlinear relationship between Cx43 gap junctional function and phenotypic outcomes.

Effects of advanced maternal age and acute prenatal alcohol exposure on mouse offspring growth and craniofacial phenotype

BACKGROUND

Prenatal alcohol exposure (PAE) can result in developmental defects that include growth restriction, craniofacial anomalies, and cognitive behavioral deficits, though the presence and severity of these adverse outcomes can vary dramatically among exposed individuals. Preclinical animal models have demonstrated that the dose and timing of PAE account for much, but not all, of this phenotypic variation, suggesting that additional factors mitigate the effects of PAE. Here, we used a mouse model to investigate whether maternal age modulates the effects of PAE on the severity and variation in offspring growth and craniofacial outcomes.

METHODS

Nulliparous C57BL/6N dams received either an intraperitoneal injection of ethanol (EtOH) or vehicle solution on gestational day 7.5. Dams were divided into four groups: (1) EtOH-treated young dams (6 to 10 weeks); (2) control young dams; (3) EtOH-treated old dams (6 to 7 months); and (4) old control dams. Neonate offspring growth restriction was measured through body mass and organ-to-body mass ratios, while skeletal craniofacial features were imaged using micro-CT and analyzed for size, shape, and variation.

RESULTS

PAE and advanced maternal age each increased the risk of low birthweight and growth restriction in offspring, but these factors in combination changed the nature of the growth restriction. Similarly, both PAE and advanced maternal age individually caused changes to craniofacial morphology such as smaller skull size, dysmorphic skull shape, and greater skull shape variation and asymmetry. Interestingly, while the combination of PAE and advanced maternal age did not affect mean skull shape or size, it significantly increased the variation and asymmetry of those measures.

CONCLUSION

Our results indicate that maternal age modulates the effects of PAE, but that the effects of this combination on offspring outcomes are more complex than simply scaling the effects of either factor.

Body shape transformations by alternate anatomical adaptive peak shifts in blenniiform fishes

Extreme body elongation has occurred repeatedly in the evolutionary history of ray-finned fishes. Lengthening of the anterior-posterior body axis relative to depth and width can involve changes in the cranial skeleton and vertebral column, but to what extent is anatomical evolution determined by selective factors and intrinsic constraints that are shared broadly among closely related lineages? In this study, we fit adaptive (Ornstein-Uhlenbeck) evolutionary models to body shape and its anatomical determinants and identified two instances of extreme elongation by divergent anatomical peak shifts in the Blenniiformes, a radiation of small-bodied substrate-associated marine teleost fishes. Species in the genus Xiphasia (hairtail blennies) evolved toward a peak defined by a highly elongated caudal vertebral region but ancestral cranial and precaudal vertebral morphology. In contrast, a clade that includes the genera Chaenopsis and Lucayablennius (pike and arrow blennies) evolved toward a peak with a long slender skull but ancestral axial skeletal anatomy. Neither set of anatomical peak shifts aligns closely with the major axis of anatomical diversification in other blenniiform fishes. These results provide little evidence that ancestral constraints have affected body shape transformation, and instead suggest that extreme elongation arose with distinct shifts in selective factors and development.

Repeatability of Adaptive Radiation Depends on Spatial Scale: Regional Versus Global Replicates of Stickleback in Lake Versus Stream Habitats

The repeatability of adaptive radiation is expected to be scale-dependent, with determinism decreasing as greater spatial separation among "replicates" leads to their increased genetic and ecological independence. Threespine stickleback (Gasterosteus aculeatus) provide an opportunity to test whether this expectation holds for the early stages of adaptive radiation-their diversification in freshwater ecosystems has been replicated many times. To better understand the repeatability of that adaptive radiation, we examined the influence of geographic scale on levels of parallel evolution by quantifying phenotypic and genetic divergence between lake and stream stickleback pairs sampled at regional (Vancouver Island) and global (North America and Europe) scales. We measured phenotypes known to show lake-stream divergence and used reduced representation genome-wide sequencing to estimate genetic divergence. We assessed the scale dependence of parallel evolution by comparing effect sizes from multivariate models and also the direction and magnitude of lake-stream divergence vectors. At the phenotypic level, parallelism was greater at the regional than the global scale. At the genetic level, putative selected loci showed greater lake-stream parallelism at the regional than the global scale. Generally, the level of parallel evolution was low at both scales, except for some key univariate traits. Divergence vectors were often orthogonal, highlighting possible ecological and genetic constraints on parallel evolution at both scales. Overall, our results confirm that the repeatability of adaptive radiation decreases at increasing spatial scales. We suggest that greater environmental heterogeneity at larger scales imposes different selection regimes, thus generating lower repeatability of adaptive radiation at larger spatial scales.