Phenotypic clines, plasticity, and morphological trade-offs in an intertidal snail

A novel analytical framework to quantify co-gradient and countergradient variation

Spatial covariance between genotypic and environmental influences on phenotypes (Cov_GE ) can result in the nonrandom distribution of genotypes across environmental gradients and is a potentially important factor driving local adaptation. However, a framework to quantify the magnitude and significance of Cov_GE has been lacking. We develop a novel quantitative/analytical approach to estimate and test the significance of Cov_GE from reciprocal transplant or common garden experiments, which we validate using simulated data. We demonstrate how power to detect Cov_GE changes over a range of experimental designs. We confirm an inverse relationship between gene-by-environment interactions (GxE) and Cov_GE , as predicted by first principles, but show how phenotypes can be influenced by both. The metric provides a way to measure how phenotypic plasticity covaries with genetic differentiation and highlights the importance of understanding the dual influences of Cov_GE and GxE on phenotypes in studies of local adaptation and species' responses to environmental change.

Radular force performance of stylommatophoran gastropods (Mollusca) with distinct body masses

The forces exerted by the animal's food processing structures can be important parameters when studying trophic specializations to specific food spectra. Even though molluscs represent the second largest animal phylum, exhibiting an incredible biodiversity accompanied by the establishment of distinct ecological niches including the foraging on a variety of ingesta types, only few studies focused on the biomechanical performance of their feeding organs. To lay a keystone for future research in this direction, we investigated the in vivo forces exerted by the molluscan food gathering and processing structure, the radula, for five stylommatophoran species (Gastropoda). The chosen species and individuals have a similar radular morphology and motion, but as they represent different body mass classes, we were enabled to relate the forces to body mass. Radular forces were measured along two axes using force transducers which allowed us to correlate forces with the distinct phases of radular motion. A radular force quotient, AFQ = mean Absolute Force/bodymass^0.67, of 4.3 could be determined which can be used further for the prediction of forces generated in Gastropoda. Additionally, some specimens were dissected and the radular musculature mass as well as the radular mass and dimensions were documented. Our results depict the positive correlation between body mass, radular musculature mass, and exerted force. Additionally, it was clearly observed that the radular motion phases, exerting the highest forces during feeding, changed with regard to the ingesta size: all smaller gastropods rather approached the food by a horizontal, sawing-like radular motion leading to the consumption of rather small food particles, whereas larger gastropods rather pulled the ingesta in vertical direction by radula and jaw resulting in the tearing of larger pieces.

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

Extended phenotypes on coral reefs: cryptic phenotypes modulate coral-vermetid interactions

Phenotypic variation can lead to variation in the strength and outcome of species interactions. Variation in phenotypic traits can arise due to plastic responses to environmental stimuli, underlying genetic variation, or both, and may reflect differences in the focal organism or aspects of the extended phenotype (e.g., associated microbes). We used a reciprocal transplant experiment of Porites corals to evaluate the role of plasticity vs. heritable diversity on phenotypic traits and performance of corals that varied in their prior exposure to vermetid gastropods, an organism known to reduce coral growth and survival. We measured a suite of phenotypic traits associated with coral performance, many of which showed a plastic response to vermetid exposure. Vermetids decreased calcification of corals, increased microbial diversity, and shifted microbial composition. Most traits also showed a signature of previous exposure environment that persisted even when exposure was reversed: i.e., under the same conditions, corals naïve to vermetids had slower calcification rates, thicker tissues, higher Symbiodiniaceae densities, and different microbiomes than corals previously exposed to vermetids. We suggest the phenotypic differences are heritable, as reefs with and without vermetids were comprised of two different mitotypes, that revealed high, consistent genetic variation. Vermetids were only found on the fast-growing coral mitotype that was characterized by thin tissue, and that likely had a history of disturbance. As extended phenotypes can have community impacts, we suggest vermetid, in addition to microbes, are part of the extended community phenotype of these corals. Coral genotypes can establish different reef trajectories, with thin-tissue types more prone to disturbance and subsequent colonization by other species, like vermetids, which can further facilitate the degradation of coral reefs. The effects of the extended phenotype of species likely influence heterogeneity across landscapes as well as temporal differences in community composition.

Uncovering the shell game with barcodes: diversity of meiofaunal Caecidae snails (Truncatelloidea, Caenogastropoda) from Central America

Caecidae is a species-rich family of microsnails with a worldwide distribution. Typical for many groups of gastropods, caecid taxonomy is largely based on overt shell characters. However, identification of species using shell characteristics is problematic due to their rather uniform, tubular shells, the presence of different growth stages, and a high degree of intraspecific variability. In the present study, a first integrative approach to caecid taxonomy is provided using light-microscopic investigation with microsculptural analyses and multi-marker barcoding, in conjunction with molecular species delineation analyses (ABGD, haplotype networks, GMYC, and bPTP). In total 132 specimens of Caecum and Meioceras collected during several sampling trips to Central America were analyzed and delineated into a minimum of 19 species to discuss putative synonyms, and supplement the original descriptions. Molecular phylogenetic analyses suggest Meiocerasnitidum and M.cubitatum should be reclassified as Caecum, and the genus Meioceras might present a junior synonym of Caecum. Meiofaunal caecids morphologically resembling C.glabrum from the Northeast Atlantic are a complex of cryptic species with independent evolutionary origins, likely associated with multiple habitat shifts to the mesopsammic environment. Caecuminvisibile Egger & Jörger, sp. nov. is formally described based on molecular diagnostic characters. This first integrative approach towards the taxonomy of Caecidae increases the known diversity, reveals the need for a reclassification of the genus Caecum and serves as a starting point for a barcoding library of the family, thereby enabling further reliable identifications of these taxonomically challenging microsnails in future studies.

Island effect in the shell phenotypic plasticity of an intertidal gastropod

The pattern of shell shape variation in populations of the top shell, Phorcus sauciatus from Madeira Island (NE Atlantic) was analysed as a function of substrate type and wave exposure, using geometric morphometrics. We hypothesized that P. sauciatus shows morphological variations to inhabit contrasting environments. Highly significant differences in shell shape were found depending on both substrate type and coastal exposure. The most marked differences were found between exposed and sheltered environments. Rounded shells in exposed environments may be explained by physiological reasons, since larger muscles are needed to ensure attachment to substrate in areas subject to harsh conditions. On the other side, conically-shaped shells may accommodate a larger body but with a smaller foot, an adaptation to sheltered environments. Slight shape variations were also observed among substrates, mostly in the degree of differentiation of some whorls. Differences could be related to a particular use of habitat and/or to the degree of exposure to water current. These results suggest that P. sauciatus is locally adapted to varying coastal hydrodynamics and may be considered as a good model in studies on adaptations of fauna to certain climate change effects.

Detecting trends in body size: empirical and statistical requirements for intraspecific analyses

Attributing biological explanations to observed ecogeographical and ecological patterns require eliminating potential statistical and sampling artifacts as alternative explanations of the observed patterns. Here, we assess the role of sample size, statistical power, and geographic inclusivity on the general validity and statistical significance of relationships between body size and latitude for 3 well-studied species of turtles. We extend those analyses to emphasize the importance of using statistically robust data in determining macroecological patterns. We examined intraspecific trends in body size with latitude in Chelydra serpentina, Chrysemys picta, and Trachemys scripta using Pearson’s correlations, diagnostic tests for influential points, and resampling. Existing data were insufficient to ascertain a latitudinal trend in body size for C. serpentina or T. scripta. There was a significant relationship for C. picta, however, resampling analyses show that, on average, 16 of the 23 available independent populations were needed to demonstrate a significant relationship and that at least 20 of 23 populations were required to obtain a statistically powerful correlation between body size and latitude. Furthermore, restricting the latitudes of populations resampled shows that body size trends of C. picta were largely due to leveraging effects of populations at the edge of the species range. Our results suggest that broad inferences regarding ecological trends in body size should be made with caution until underlying (intraspecific) patterns in body size can be statistically and conclusively demonstrated.

Evolutionary lineages of marine snails identified using molecular phylogenetics and geometric morphometric analysis of shells

The relationship between morphology and inheritance is of perennial interest in evolutionary biology and palaeontology. Using three marine snail genera Penion, Antarctoneptunea and Kelletia, we investigate whether systematics based on shell morphology accurately reflect evolutionary lineages indicated by molecular phylogenetics. Members of these gastropod genera have been a taxonomic challenge due to substantial variation in shell morphology, conservative radular and soft tissue morphology, few known ecological differences, and geographical overlap between numerous species. Sampling all sixteen putative taxa identified across the three genera, we infer mitochondrial and nuclear ribosomal DNA phylogenetic relationships within the group, and compare this to variation in adult shell shape and size. Results of phylogenetic analysis indicate that each genus is monophyletic, although the status of some phylogenetically derived and likely more recently evolved taxa within Penion is uncertain. The recently described species P. lineatus is supported by genetic evidence. Morphology, captured using geometric morphometric analysis, distinguishes the genera and matches the molecular phylogeny, although using the same dataset, species and phylogenetic subclades are not identified with high accuracy. Overall, despite abundant variation, we find that shell morphology accurately reflects genus-level classification and the corresponding deep phylogenetic splits identified in this group of marine snails.

The Fossil Record of Predator-Prey Arms Races: Coevolution and Escalation Hypotheses

Arms races between predators and prey may be driven by two related processes—escalation and coevolution. Escalation is enemy-driven evolution. In this top-down view of an arms race, the role of prey (with the exception of dangerous prey) is downplayed. In coevolution, two or more species change reciprocally in response to one another; prey are thought to drive the evolution of their predator, and vice versa. In the fossil record, the two processes are most reliably distinguished when the predator-prey system is viewed within the context of the other species that may influence the interaction, thus allowing for a relative ranking of the importance of selective agents. Detailed documentation of the natural history of living predator-prey systems is recommended in order to distinguish the processes in some fossil systems. A geographic view of species interactions and the processes driving their evolution may lead to a more diverse array of testable hypotheses on how predator-prey systems evolve and what constraints interactions impose on the evolution of organisms. Scale is important in evaluating the role of escalation and coevolution in the evolution of species interactions. If short-term reciprocal adaptation (via phenotypic plasticity or selection mosaics among populations) between predator and prey is a common process, then prey are likely to exert some selective pressure over their predators over the short term (on ecological time scales), but in the long run predators may still exert primary “top-down” control in directing evolution. On the scale of evolutionary time, predators of large effect likely control the overall directionality of evolution due to the inequalities of predator and prey in control of resources.

Disentangling the genetic and morphological structure of Patella candei complex in Macaronesia (NE Atlantic)

The uptake of natural living resources for human consumption has triggered serious changes in the balance of ecosystems. In the archipelagos of Macaronesia (NE Atlantic), limpets have been extensively exploited probably since islands were first colonized. This has led to profound consequences in the dynamics of rocky shore communities. The Patella candei complex includes various subspecies of limpets that are ascribed to a particular archipelago and has been the focus of several taxonomic surveys without much agreement. Under a conservational perspective, we apply morphometric and genetic analyses to test subspecies boundaries in P. candei and to evaluate its current population connectivity throughout Macaronesia (Azores, Madeira, and Canaries). A highly significant genetic break between archipelagos following isolation by distance was detected (F_ST = 0.369, p < .001). Contrastingly, significant genetic differentiation among islands (i.e., Azores) was absent possibly indicating ongoing gene flow via larval exchange between populations. Significant shell‐shape differences among archipelagos were also detected using both distance‐based and geometric morphometric analyses. Adaptive processes associated with niche differentiation and strong barriers to gene flow among archipelagos may be the mechanisms underlying P. candei diversification in Macaronesia. Under the very probable assumption that populations of P. candei from each archipelago are geographically and/or ecologically isolated populations, the various subspecies within the P. candei complex may be best thought of as true species using the denomination: P. candei in Selvagens, Patella gomesii in Azores, Patella ordinaria in Madeira, and Patella crenata for Canaries. This would be in agreement with stock delimitation and units of conservation of P. candei sensu latu along Macaronesia.

Genetic by environmental variation but no local adaptation in oysters (Crassostrea virginica)

Functional trait variation within and across populations can strongly influence population, community, and ecosystem processes, but the relative contributions of genetic vs. environmental factors to this variation are often not clear, potentially complicating conservation and restoration efforts. For example, local adaptation, a particular type of genetic by environmental (G*E) interaction in which the fitness of a population in its own habitat is greater than in other habitats, is often invoked in management practices, even in the absence of supporting evidence. Despite increasing attention to the potential for G*E interactions, few studies have tested multiple populations and environments simultaneously, limiting our understanding of the spatial consistency in patterns of adaptive genetic variation. In addition, few studies explicitly differentiate adaptation in response to predation from other biological and environmental factors. We conducted a reciprocal transplant experiment of first‐generation eastern oyster (Crassostrea virginica) juveniles from six populations across three field sites spanning 1000 km in the southeastern Atlantic Bight in both the presence and absence of predation to test for G*E variation in this economically valuable and ecologically important species. We documented significant G*E variation in survival and growth, yet there was no evidence for local adaptation. Condition varied across oyster cohorts: Offspring of northern populations had better condition than offspring from the center of our region. Oyster populations in the southeastern Atlantic Bight differ in juvenile survival, growth, and condition, yet offspring from local broodstock do not have higher survival or growth than those from farther away. In the absence of population‐specific performance information, oyster restoration and aquaculture may benefit from incorporating multiple populations into their practices.

Shaving a Shell: Effect of Manipulated Sculpture and Feeding on Shell Growth and Sculpture Development in Nucella lamellosa (Muricidae: Ocenebrinae)

Gastropod shell sculpture offers a novel tool for studying morphological patterning. Existing shell features may be manipulated experimentally to test how alteration affects subsequent shell growth and form. Axial sculpture occurs in many gastropod groups, and spacing of sculpture may be regular or irregular. But how gastropods control sculpture placement during shell growth is unknown. We studied the growth and positioning of axial lamellae in the muricid Nucella lamellosa, and compared these to the superficially similar axial varices seen in other muricids. First, we tested whether the feeding rate had any effect on the rate of addition or positioning of new lamellae. Second, we tested what effect previous shell sculpture had on lamellar placement, and shell growth in general, by removing all shell sculpture and allowing snails to grow over the "shaved" shell surface. Lamellar growth appeared to be relatively plastic; spacing was highly variable both within and among individual snails, and 1-2 weeks were required to complete the addition of a new lamella. Body growth rate was the primary determinant of lamellar growth; past lamellae had no effect on placement of new lamellae or rate of shell length increase. Feeding rate and body size affected only growth in shell length, and had no direct effect on spacing or on the rate of addition of new lamellae. The growth of axial lamellae in N. lamellosa differed from that of varices by exhibiting neither a) regular spacing nor b) a growth hiatus after completion of a lamella. Significantly, despite the obvious impediment of previous sculpture to future shell growth, removal of this sculpture had no observable effect on the rate of body growth or on any aspect of subsequent lamellar growth.

What can aquatic gastropods tell us about phenotypic plasticity? A review and meta-analysis

There have been few attempts to synthesise the growing body of literature on phenotypic plasticity to reveal patterns and generalities about the extent and magnitude of plastic responses. Here, we conduct a review and meta-analysis of published literature on phenotypic plasticity in aquatic (marine and freshwater) gastropods, a common system for studying plasticity. We identified 96 studies, using pre-determined search terms, published between 1985 and November 2013. The literature was dominated by studies of predator-induced shell form, snail growth rates and life history parameters of a few model taxa, accounting for 67% of all studies reviewed. Meta-analyses indicated average plastic responses in shell thickness, shell shape, and growth and fecundity of freshwater species was at least three times larger than in marine species. Within marine gastropods, species with planktonic development had similar average plastic responses to species with benthic development. We discuss these findings in the context of the role of costs and limits of phenotypic plasticity and environmental heterogeneity as important constraints on the evolution of plasticity. We also consider potential publication biases and discuss areas for future research, indicating well-studied areas and important knowledge gaps.

Evolutionary tradeoffs, Pareto optimality and the morphology of ammonite shells

Background

Organisms that need to perform multiple tasks face a fundamental tradeoff: no design can be optimal at all tasks at once. Recent theory based on Pareto optimality showed that such tradeoffs lead to a highly defined range of phenotypes, which lie in low-dimensional polyhedra in the space of traits. The vertices of these polyhedra are called archetypes- the phenotypes that are optimal at a single task. To rigorously test this theory requires measurements of thousands of species over hundreds of millions of years of evolution. Ammonoid fossil shells provide an excellent model system for this purpose. Ammonoids have a well-defined geometry that can be parameterized using three dimensionless features of their logarithmic-spiral-shaped shells. Their evolutionary history includes repeated mass extinctions.

Results

We find that ammonoids fill out a pyramid in morphospace, suggesting five specific tasks - one for each vertex of the pyramid. After mass extinctions, surviving species evolve to refill essentially the same pyramid, suggesting that the tasks are unchanging. We infer putative tasks for each archetype, related to economy of shell material, rapid shell growth, hydrodynamics and compactness.

Conclusions

These results support Pareto optimality theory as an approach to study evolutionary tradeoffs, and demonstrate how this approach can be used to infer the putative tasks that may shape the natural selection of phenotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0149-z) contains supplementary material, which is available to authorized users.

Geometric morphometric analysis of the head of Microlophus atacamensis (Tropiduridae) in a latitudinal gradient

BACKGROUND

Clinal variation is defined as gradual variation in a character associated with geographic distance among sites. Microlophus atacamensis is a medium large lizard species which inhabits the intertidal zone of northern Chile, distributed from Antofagasta (23° 39' S) to Arrayán, La Serena (29° 41' S) in a gradient which could show clinal variation. Geometric morphometrics analyzes differences in shape independent of size; information about shape allows a more complete biological interpretation than information on morphological variation. This studyevaluated clinal variation in the head shape of M.atacamensis fromfive localities (Antofagasta, Paposo, Pan de Azúcar, Caldera, and Tres Playitas) using dorsal and lateral views, comparing form variation with latitudinal distribution. The heads of adults collected were photographed in lateral and dorsal views.

RESULTS

The analysis did not find significant differences in form among the five localities, in contrast to the proposal of earlier studies, and no differences were recorded between the sexes. Possible reasons why these populations are not differentiated in the latitudinal gradient are discussed.

CONCLUSIONS

Our results show that there are no differences between the studied, among the explanations are that populationsof this species are subjected to similar environments that promote convergence of the structures analyzed.

AFLPs and mitochondrial haplotypes reveal local adaptation to extreme thermal environments in a freshwater gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.

Developmental thermal plasticity among Drosophila melanogaster populations

Many biotic and abiotic variables influence the dispersal and distribution of organisms. Temperature has a major role in determining these patterns because it changes daily, seasonally and spatially, and these fluctuations have a significant impact on an organism's behaviour and fitness. Most ecologically relevant phenotypes that are adaptive are also complex and thus they are influenced by many underlying loci that interact with the environment. In this study, we quantified the degree of thermal phenotypic plasticity within and among populations by measuring chill-coma recovery times of lines reared from egg to adult at two different environmental temperatures. We used sixty genotypes from six natural populations of Drosophila melanogaster sampled along a latitudinal gradient in South America. We found significant variation in thermal plasticity both within and among populations. All populations exhibit a cold acclimation response, with flies reared at lower temperatures having increased resistance to cold. We tested a series of environmental parameters against the variation in population mean thermal plasticity and discovered the mean thermal plasticity was significantly correlated with altitude of origin of the population. Pairing our data with previous experiments on viability fitness assays in the same populations in fixed and variable environments suggests an adaptive role of this thermal plasticity in variable laboratory environments. Altogether, these data demonstrate abundant variation in adaptive thermal plasticity within and among populations.

The cold-water connection: Bergmann's rule in North American freshwater fishes

Understanding general rules governing macroecological body size variations is one of the oldest pursuits in ecology. However, this science has been dominated by studies of terrestrial vertebrates, spurring debate over the validity of such rules in other taxonomic groups. Here, relationships between maximum body size and latitude, temperature, and elevation were evaluated for 29 North American freshwater fish species. Bergmann's rule (i.e., that body size correlates positively with latitude and negatively with temperature) was observed in 38% of species, converse Bergmann's rule (that body size correlates negatively with latitude and positively with temperature) was observed in 34% of species, and 28% of species showed no macroecological body size relationships. Most notably, every species that expressed Bergmann's rule was a cool- or cold-water species while every species that expressed converse Bergmann's rule was a warm-water species, highlighting how these patterns are likely connected to species thermal niches. This study contradicts previous research suggesting Bergmann's rule does not apply to freshwater fishes, and is congruent with an emerging paradigm of variable macroecological body size patterns in poikilotherms.

The geometry of the Pareto front in biological phenotype space

When organisms perform a single task, selection leads to phenotypes that maximize performance at that task. When organisms need to perform multiple tasks, a trade-off arises because no phenotype can optimize all tasks. Recent work addressed this question, and assumed that the performance at each task decays with distance in trait space from the best phenotype at that task. Under this assumption, the best-fitness solutions (termed the Pareto front) lie on simple low-dimensional shapes in trait space: line segments, triangles and other polygons. The vertices of these polygons are specialists at a single task. Here, we generalize this finding, by considering performance functions of general form, not necessarily functions that decay monotonically with distance from their peak. We find that, except for performance functions with highly eccentric contours, simple shapes in phenotype space are still found, but with mildly curving edges instead of straight ones. In a wide range of systems, complex data on multiple quantitative traits, which might be expected to fill a high-dimensional phenotype space, is predicted instead to collapse onto low-dimensional shapes; phenotypes near the vertices of these shapes are predicted to be specialists, and can thus suggest which tasks may be at play.

Marine invertebrate skeleton size varies with latitude, temperature and carbonate saturation: implications for global change and ocean acidification

There is great concern over the future effects of ocean acidification on marine organisms, especially for skeletal calcification, yet little is known of natural variation in skeleton size and composition across the globe, and this is a prerequisite for identifying factors currently controlling skeleton mass and thickness. Here, taxonomically controlled latitudinal variations in shell morphology and composition were investigated in bivalve and gastropod molluscs, brachiopods, and echinoids. Total inorganic content, a proxy for skeletal CaCO₃ , decreased with latitude, decreasing seawater temperature, and decreasing seawater carbonate saturation state (for CaCO₃ as calcite (Ω_cal )) in all taxa. Shell mass decreased with latitude in molluscs and shell inorganic content decreased with latitude in buccinid gastropods. Shell thickness decreased with latitude in buccinid gastropods (excepting the Australian temperate buccinid) and echinoids, but not brachiopods and laternulid clams. In the latter, the polar species had the thickest shell. There was no latitudinal trend in shell thickness within brachiopods. The variation in trends in shell thickness by taxon suggests that in some circumstances ecological factors may override latitudinal trends. Latitudinal gradients may produce effects similar to those of future CO₂ -driven ocean acidification on CaCO₃ saturation state. Responses to latitudinal trends in temperature and saturation state may therefore be useful in informing predictions of organism responses to ocean acidification over long-term adaptive timescales.

Intraspecific trait cospecialization of constitutive and inducible morphological defences in a marine snail from habitats with different predation risk

1. Studies examining the integration of constitutive and inducible aspects of multivariate defensive phenotypes are rare. 2. I asked whether marine snails (Nucella lamellosa) from habitats with and without abundant predatory crabs differed in constitutive and inducible aspects of defensive shell morphology. 3. I examined multivariate shell shape development of snails from each habitat in the presence and absence of waterborne cues from feeding crabs (Cancer productus). I also examined the influence of constitutive and inducible shell morphology on resistance to crushing. 4. Regardless of the presence of crabs, snails from high-risk (HR) habitats developed rotund, short-spired shells, while snails from low-risk habitats developed elongate shells, tall-spired shells, indicating among-habitat divergence in constitutive shell shape. Moreover, allometry analyses indicated that constitutive developmental patterns underlying this variation also differed between habitats. However, snails from HR habitats showed greater plasticity for apertural lip thickness and apertural area in the presence of crab cues, indicating among-habitat variation in defence inducibility. 5. Both shell shape and apertural lip thickness contributed to shell strength suggesting that constitutive shell shape development and inducible lip thickening have evolved jointly to form an effective defence in habitats where predation risk is high.

Life-history variation with respect to experienced thermal environments in the lizard, Eremias multiocellata (Lacertidae)

We compared adult size, female reproductive traits, and offspring phenotypes between multiocellated racerunners (Eremias multiocellata) from two thermally different sites (populations) in Inner Mongolia (North China): the colder one in Wulatehouqi (WQ) and the warmer one in Dalateqi (DQ). Both adults and neonates were smaller in the colder site. Females from the two sites both produced a single litter of 2-5 young per season, and did not differ in allocation of energy to reproduction after accounting for differences in body size. Female neonates had more ventral scales than did males, and the WQ neonates had fewer ventral scales than did the DQ neonates. The WQ neonates were slower than the DQ neonates. When body length was normalized across populations, we found that (1) hindlimb length correlated positively with sprint speed in both WQ and DQ neonates, (2) forelimb length correlated positively with sprint speed only in the DQ neonates, and (3) tail length correlated positively with sprint speed only in the WQ neonates. Hindlimb length played a more important role in locomotion than did tail length or forelimb length. Though differing in size and morphology, neonates from the two sites did not differ in early growth and survival under identical laboratory conditions. Our data are consistent with many studies that have shown countergradient variation in physiological traits (growth rate and reproductive output) and cogradient variation in morphological traits.

Spatial geographic mosaic in an aquatic predator-prey network

The geographic mosaic theory of coevolution predicts 1) spatial variation in predatory structures as well as prey defensive traits, and 2) trait matching in some areas and trait mismatching in others mediated by gene flow. We examined gene flow and documented spatial variation in crushing resistance in the freshwater snails Mexipyrgus churinceanus, Mexithauma quadripaludium, Nymphophilus minckleyi, and its relationship to the relative frequency of the crushing morphotype in the trophically polymorphic fish Herichthys minckleyi. Crushing resistance and the frequency of the crushing morphotype did show spatial variation among 11 naturally replicated communities in the Cuatro Ciénegas valley in Mexico where these species are all endemic. The variation in crushing resistance among populations was not explained by geographic proximity or by genetic similarity in any species. We detected clear phylogeographic patterns and limited gene flow for the snails but not for the fish. Gene flow among snail populations in Cuatro Ciénegas could explain the mosaic of local divergence in shell strength and be preventing the fixation of the crushing morphotype in Herichthys minckleyi. Finally, consistent with trait matching across the mosaic, the frequency of the fish morphotype was negatively correlated with shell crushing resistance likely reflecting the relative disadvantage of the crushing morphotype in communities where the snails exhibit relatively high crushing resistance.

Experimental evidence for the rapid evolution of behavioral canalization in natural populations

Canalization-the evolutionary loss of the capacity of organisms to develop different phenotypes in different environments-is an evolutionary phenomenon suspected to occur widely, although examples in natural populations are elusive. Because behavior is typically a highly flexible component of an individual's phenotype, it provides fertile ground for studying the evolution of canalization. Here we report how snail populations exposed for different lengths of time to a predatory crab introduced from Europe to America exhibit different degrees of canalization of an adaptive antipredator behavior: soft tissue withdrawal, measured as angular retraction depth. Where crab-snail contact is shortest (60 years), snails showed the highest behavioral flexibility. Where crabs invaded 110 years ago, snails showed significantly less behavioral flexibility, and where the interaction is ancient (Europe), snails exhibited highly canalized behavior. Selection therefore appears to have acted rapidly to increase canalization in wild snail populations, leading ultimately to the hard-wired behavior seen in European conspecifics.

The covariance between genetic and environmental influences across ecological gradients: reassessing the evolutionary significance of countergradient and cogradient variation

Patterns of phenotypic change across environmental gradients (e.g., latitude, altitude) have long captivated the interest of evolutionary ecologists. The pattern and magnitude of phenotypic change is determined by the covariance between genetic and environmental influences across a gradient. Cogradient variation (CoGV) occurs when covariance is positive: that is, genetic and environmental influences on phenotypic expression are aligned and their joint influence accentuates the change in mean trait value across the gradient. Conversely, countergradient variation (CnGV) occurs when covariance is negative: that is, genetic and environmental influences on phenotypes oppose one another, thereby diminishing the change in mean trait expression across the gradient. CnGV has so far been found in at least 60 species, with most examples coming from fishes, amphibians, and insects across latitudinal or altitudinal gradients. Traits that display CnGV most often involve metabolic compensation, that is, the elevation of various physiological rates processes (development, growth, feeding, metabolism, activity) to counteract the dampening effect of reduced temperature, growing season length, or food supply. Far fewer examples of CoGV have been identified (11 species), and these most often involve morphological characters. Increased knowledge of spatial covariance patterns has furthered our understanding of Bergmann size clines, phenotypic plasticity, species range limits, tradeoffs in juvenile growth rate, and the design of conservation strategies for wild species. Moreover, temporal CnGV explains some cases of an apparent lack of phenotypic response to directional selection and provides a framework for predicting evolutionary responses to climate change.

Divergent shell shape as an antipredator adaptation in tropical land snails

Although many land snails exhibit amazingly divergent shell shapes in the tropics, the functions of these remain obscure. Here we show that a modified aperture shape acts as an impediment specifically to predation by a snail-eating snake. Pareas iwasakii (Colubridae: Pareatinae) uses a unique method to feed on land snails: the snake extracts the soft body from the shell through the aperture by alternately retracting its mandibles. The snail Satsuma caliginosa (Camaenidae: Camaeninae) has apertural variation in regard to the presence of snail-eating snakes. Our experiments demonstrated that the distorted aperture mechanically impeded predation by this gape-limited predator, interrupting the mandibular movements. In contrast, congeneric snails with round apertures did not escape predation by snakes. The paleobiogeography of the focal area indicates that the subspecies Satsuma caliginosa picta, which does not have apertural modification, was derived from a defensive ancestor after the extinction of snail-eating snakes. Our study suggests a possibility that snail-eating snakes are responsible for divergent shell shapes in a variety of tropical land snails.