Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

Interspecific allometry for sexual shape dimorphism: Macroevolution of multivariate sexual phenotypes with application to Rensch's rule

Allometric trends in the degree of sexual dimorphism with body size have long fascinated evolutionary biologists. Many male-biased clades display more prominent sexual dimorphism in larger taxa (Rensch's rule), with most examples documenting this pattern for body size dimorphism. Although sexual dimorphism in traits other than body size is equally functionally relevant, characterizing allometric patterns of sexual dimorphism in such traits is hampered by lack of an analytical framework that can accommodate multivariate phenotypes. In this article, we derive a multivariate equivalency for investigating trends in sexual dimorphism—relative to overall body size—across taxa and provide a generalized test to determine whether such allometric patterns correspond with Rensch's rule. For univariate linear traits such as body size, our approach yields equivalent results to those from standard procedures, but our test is also capable of detecting trends in multivariate datasets such as shape. Computer simulations reveal that the method displays appropriate statistical properties, and an empirical example in Mediterranean lizards provides the first demonstration of Rensch's rule in a multivariate phenotype (head shape). Our generalized procedure substantially extends the analytical toolkit for investigating macroevolutionary patterns of sexual dimorphism and seeking a better understanding of the processes that underlie them.     

Relationships between head morphology, bite performance and ecology in two species of Podarcis wall lizards

Understanding the relationship between form and function is central to our comprehension of how phenotypic diversity evolves. Traits involved in multiple activities, such as social interactions and ecological resource use, are under the influence of different evolutionary forces potentially acting in opposite directions. Such systems provide the opportunity of understanding how potential constraints on morphological variation may influence whole-organism performance. In this study we examined morphology and bite performance in two closely related species of Podarcis wall lizards with divergent microhabitat preferences, to investigate how natural and sexual selection interact to shape the evolution of head traits. Our results show that although head morphology is markedly different between species and sexes, only sexes differ in bite force, indicating that the ecological differentiation between species is reflected in their morphology but does not constrain performance. Rather, the modification of the relative size of head components between species and a shift in the form-function relationship provide a potential explanation of how equal performance is attained by different morphological configurations. Geometric morphometrics provide a clear, biomechanically meaningful image of how this is achieved and show a bisexual pattern of head shape-bite force association in both species. This, together with a strong allometry of head size on body size and head shape on head size, provides indirect morphological evidence for the importance of sexual selection in shaping morphological and functional patterns. Finally, our findings suggest that the differences observed between species and sexes in head traits and bite performance are not reflected in their dietary ecology, implying that if trophic niche segregation between groups occurs, the reasons behind it are not primarily related to head morphology and functional variation.     

Ecomorphological variation in male and female wall lizards and the macroevolution of sexual dimorphism in relation to habitat use

Understanding how phenotypic diversity evolves is a major interest of evolutionary biology. Habitat use is an important factor in the evolution of phenotypic diversity of many animal species. Interestingly, male and female phenotypes have been frequently shown to respond differently to environmental variation. At the macroevolutionary level, this difference between the sexes is frequently analysed using phylogenetic comparative tools to assess variation in sexual dimorphism (SD) across taxa in relation to habitat. A shortcoming of such analyses is that they evaluate the degree of dimorphism itself and therefore they do not provide access to the evolutionary trajectories of each sex. As such, the relative contribution of male and female phenotypes on macroevolutionary patterns of sexual dimorphism cannot be directly assessed. Here, we investigate how habitat use shapes phenotypic diversity in wall lizards using phylogenetic comparative tools to simultaneously assess the tempo and mode of evolution in males, females and the degree of sexual dimorphism. We find that both sexes have globally diversified under similar, but not identical, processes, where habitat use seems to drive macroevolutionary variation in head shape, but not in body size or relative limb length. However, we also observe small differences in the evolutionary dynamics of male and female phenotypes that have a marked impact on macroevolutionary patterns of SD, with important implications for our interpretation of what drives phenotypic diversification within and between the sexes.     

Hybridization patterns between two marine snails,Littorina fabalis and L. obtusata

Characterizing the patterns of hybridization between closely related species is crucial to understand the role of gene flow in speciation. In particular, systems comprising multiple contacts between sister species offer an outstanding opportunity to investigate how reproductive isolation varies with environmental conditions, demography and geographic contexts of divergence. The flat periwinkles, Littorina obtusata and L. fabalis (Gastropoda), are two intertidal sister species with marked ecological differences compatible with late stages of speciation. Although hybridization between the two was previously suggested, its extent across the Atlantic shores of Europe remained largely unknown. Here, we combined genetic (microsatellites and mtDNA) and morphological data (shell and male genital morphology) from multiple populations of flat periwinkles in north‐western Iberia to assess the extent of current and past hybridization between L. obtusata and L. fabalis under two contrasting geographic settings of divergence (sympatry and allopatry). Hybridization signatures based on both mtDNA and microsatellites were stronger in sympatric sites, although evidence for recent extensive admixture was found in a single location. Misidentification of individuals into species based on shell morphology was higher in sympatric than in allopatric sites. However, despite hybridization, species distinctiveness based on this phenotypic trait together with male genital morphology remained relatively high. The observed variation in the extent of hybridization among locations provides a rare opportunity for future studies on the consequences of different levels of gene flow for reinforcement, thus informing about the mechanisms underlying the completion of speciation.     

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

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

Comparing Multiple Criteria for Species Identification in Two Recently Diverged Seabirds

Correct species identification is a crucial issue in systematics with key implications for prioritising conservation effort. However, it can be particularly challenging in recently diverged species due to their strong similarity and relatedness. In such cases, species identification requires multiple and integrative approaches. In this study we used multiple criteria, namely plumage colouration, biometric measurements, geometric morphometrics, stable isotopes analysis (SIA) and genetics (mtDNA), to identify the species of 107 bycatch birds from two closely related seabird species, the Balearic (Puffinus mauretanicus) and Yelkouan (P. yelkouan) shearwaters. Biometric measurements, stable isotopes and genetic data produced two stable clusters of bycatch birds matching the two study species, as indicated by reference birds of known origin. Geometric morphometrics was excluded as a species identification criterion since the two clusters were not stable. The combination of plumage colouration, linear biometrics, stable isotope and genetic criteria was crucial to infer the species of 103 of the bycatch specimens. In the present study, particularly SIA emerged as a powerful criterion for species identification, but temporal stability of the isotopic values is critical for this purpose. Indeed, we found some variability in stable isotope values over the years within each species, but species differences explained most of the variance in the isotopic data. Yet this result pinpoints the importance of examining sources of variability in the isotopic data in a case-by-case basis prior to the cross-application of the SIA approach to other species. Our findings illustrate how the integration of several methodological approaches can help to correctly identify individuals from recently diverged species, as each criterion measures different biological phenomena and species divergence is not expressed simultaneously in all biological traits.     

Locomotor Mode and the Evolution of the Hindlimb in Western Mediterranean Anurans

Evolutionary Biology - The evolutionary association between morphology, locomotor performance and habitat use is a central element of the ecomorphological paradigm, and it is known to underlie the...