CONTRASTING PATTERNS OF HERITABLE GEOGRAPHIC VARIATION IN SHELL MORPHOLOGY AND GROWTH POTENTIAL IN THE MARINE GASTROPODBEMBICIUM VITTATUM: EVIDENCE FROM FIELD EXPERIMENTS

Deciphering mollusc shell production: the roles of genetic mechanisms through to ecology, aquaculture and biomimetics

Most molluscs possess shells, constructed from a vast array of microstructures and architectures. The fully formed shell is composed of calcite or aragonite. These CaCO₃ crystals form complex biocomposites with proteins, which although typically less than 5% of total shell mass, play significant roles in determining shell microstructure. Despite much research effort, large knowledge gaps remain in how molluscs construct and maintain their shells, and how they produce such a great diversity of forms. Here we synthesize results on how shell shape, microstructure, composition and organic content vary among, and within, species in response to numerous biotic and abiotic factors. At the local level, temperature, food supply and predation cues significantly affect shell morphology, whilst salinity has a much stronger influence across latitudes. Moreover, we emphasize how advances in genomic technologies [e.g. restriction site-associated DNA sequencing (RAD-Seq) and epigenetics] allow detailed examinations of whether morphological changes result from phenotypic plasticity or genetic adaptation, or a combination of these. RAD-Seq has already identified single nucleotide polymorphisms associated with temperature and aquaculture practices, whilst epigenetic processes have been shown significantly to modify shell construction to local conditions in, for example, Antarctica and New Zealand. We also synthesize results on the costs of shell construction and explore how these affect energetic trade-offs in animal metabolism. The cellular costs are still debated, with CaCO₃ precipitation estimates ranging from 1-2 J/mg to 17-55 J/mg depending on experimental and environmental conditions. However, organic components are more expensive (~29 J/mg) and recent data indicate transmembrane calcium ion transporters can involve considerable costs. This review emphasizes the role that molecular analyses have played in demonstrating multiple evolutionary origins of biomineralization genes. Although these are characterized by lineage-specific proteins and unique combinations of co-opted genes, a small set of protein domains have been identified as a conserved biomineralization tool box. We further highlight the use of sequence data sets in providing candidate genes for in situ localization and protein function studies. The former has elucidated gene expression modularity in mantle tissue, improving understanding of the diversity of shell morphology synthesis. RNA interference (RNAi) and clustered regularly interspersed short palindromic repeats - CRISPR-associated protein 9 (CRISPR-Cas9) experiments have provided proof of concept for use in the functional investigation of mollusc gene sequences, showing for example that Pif (aragonite-binding) protein plays a significant role in structured nacre crystal growth and that the Lsdia1 gene sets shell chirality in Lymnaea stagnalis. Much research has focused on the impacts of ocean acidification on molluscs. Initial studies were predominantly pessimistic for future molluscan biodiversity. However, more sophisticated experiments incorporating selective breeding and multiple generations are identifying subtle effects and that variability within mollusc genomes has potential for adaption to future conditions. Furthermore, we highlight recent historical studies based on museum collections that demonstrate a greater resilience of molluscs to climate change compared with experimental data. The future of mollusc research lies not solely with ecological investigations into biodiversity, and this review synthesizes knowledge across disciplines to understand biomineralization. It spans research ranging from evolution and development, through predictions of biodiversity prospects and future-proofing of aquaculture to identifying new biomimetic opportunities and societal benefits from recycling shell products.

Estimating Trait Heritability in Highly Fecund Species

Increasingly, researchers are interested in estimating the heritability of traits for nonmodel organisms. However, estimating the heritability of these traits presents both experimental and statistical challenges, which typically arise from logistical difficulties associated with rearing large numbers of families independently in the field, a lack of known pedigree, the need to account for group or batch effects, etc. Here we develop both an empirical and computational methodology for estimating the narrow-sense heritability of traits for highly fecund species. Our experimental approach controls for undesirable culturing effects while minimizing culture numbers, increasing feasibility in the field. Our statistical approach accounts for known issues with model-selection by using a permutation test to calculate significance values and includes both fitting and power calculation methods. We further demonstrate that even with moderately high sample-sizes, the p-values derived from asymptotic properties of the likelihood ratio test are overly conservative, thus reducing statistical power. We illustrate our methodology by estimating the narrow-sense heritability for larval settlement, a key life-history trait, in the reef-building coral Orbicella faveolata. The experimental, statistical, and computational methods, along with all of the data from this study, are available in the R package multiDimBio.

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.

Climate change in the oceans: evolutionary versus phenotypically plastic responses of marine animals and plants

I summarize marine studies on plastic versus adaptive responses to global change. Due to the lack of time series, this review focuses largely on the potential for adaptive evolution in marine animals and plants. The approaches were mainly synchronic comparisons of phenotypically divergent populations, substituting spatial contrasts in temperature or CO2 environments for temporal changes, or in assessments of adaptive genetic diversity within populations for traits important under global change. The available literature is biased towards gastropods, crustaceans, cnidarians and macroalgae. Focal traits were mostly environmental tolerances, which correspond to phenotypic buffering, a plasticity type that maintains a functional phenotype despite external disturbance. Almost all studies address coastal species that are already today exposed to fluctuations in temperature, pH and oxygen levels. Recommendations for future research include (i) initiation and analyses of observational and experimental temporal studies encompassing diverse phenotypic traits (including diapausing cues, dispersal traits, reproductive timing, morphology) (ii) quantification of nongenetic trans-generational effects along with components of additive genetic variance (iii) adaptive changes in microbe-host associations under the holobiont model in response to global change (iv) evolution of plasticity patterns under increasingly fluctuating environments and extreme conditions and (v) joint consideration of demography and evolutionary adaptation in evolutionary rescue approaches.

Differential shell strength of Cepaea nemoralis colour morphs--implications for their anti-predator defence

One of the most spectacular evolutionary forces is predation, evidenced to stimulate polymorphism in many prey species. Shell colour polymorphism of the land snail Cepaea nemoralis is a well-known model in evolutionary research. Nevertheless, the knowledge on the ecological causes driving its evolution remains incomplete and proximal factors shaping predatory pressure on C. nemoralis morphs are unknown. We evaluated shell crushing resistance and thickness, constituting crucial snail anti-predator defences in two shell areas (the apex and labium) of eight C. nemoralis morphotypes differing in shell colour and banding pattern. A GLM showed a significant effect of shell colour, banding pattern and shell thickness on shell strength. Pink shells were stronger than yellow ones, and banded forms had stronger shells than unbanded snails. The labium (usually attacked by mice) was generally thicker and more resistant than the apex (usually crushed by birds). Thicker shells were more resistant to crushing, and the rate of shell strength increase per unit of shell thickness was greater in pink and banded individuals compared to yellow and unbanded ones. Yellow and unbanded morphs have been found to be preferred by mice in the previous studies, which suggests that shell strength may be an important trait used in prey selection by these shell-crushing predators. The differences in potential anti-predator defences among snail morphs, found in the present study, justify future research on direct effect of C. nemoralis morphs shell strength on predator selectivity.

Decadal Changes in the World's Coastal Latitudinal Temperature Gradients

Most of the world's living marine resources inhabit coastal environments, where average thermal conditions change predictably with latitude. These coastal latitudinal temperature gradients (CLTG) coincide with important ecological clines,e.g., in marine species diversity or adaptive genetic variations, but how tightly thermal and ecological gradients are linked remains unclear. A first step is to consistently characterize the world's CLTGs. We extracted coastal cells from a global 1°×1° dataset of weekly sea surface temperatures (SST, 1982–2012) to quantify spatial and temporal variability of the world's 11 major CLTGs. Gradient strength, i.e., the slope of the linear mean-SST/latitude relationship, varied 3-fold between the steepest (North-American Atlantic and Asian Pacific gradients: −0.91°C and −0.68°C lat⁻¹, respectively) and weakest CLTGs (African Indian Ocean and the South- and North-American Pacific gradients: −0.28, −0.29, −0.32°C lat⁻¹, respectively). Analyzing CLTG strength by year revealed that seven gradients have weakened by 3–10% over the past three decades due to increased warming at high compared to low latitudes. Almost the entire South-American Pacific gradient (6–47°S), however, has considerably cooled over the study period (−0.3 to −1.7°C, 31 years), and the substantial weakening of the North-American Atlantic gradient (−10%) was due to warming at high latitudes (42–60°N, +0.8 to +1.6°C,31 years) and significant mid-latitude cooling (Florida to Cape Hatteras 26–35°N, −0.5 to −2.2°C, 31 years). Average SST trends rarely resulted from uniform shifts throughout the year; instead individual seasonal warming or cooling patterns elicited the observed changes in annual means. This is consistent with our finding of increased seasonality (i.e., summer-winter SST amplitude) in three quarters of all coastal cells (331 of 433). Our study highlights the regionally variable footprint of global climate change, while emphasizing ecological implications of changing CLTGs, which are likely driving observed spatial and temporal clines in coastal marine life.

Variation in size and growth of the great scallop Pecten maximus along a latitudinal gradient

Understanding the relationship between growth and temperature will aid in the evaluation of thermal stress and threats to ectotherms in the context of anticipated climate changes. Most Pecten maximus scallops living at high latitudes in the northern hemisphere have a larger maximum body size than individuals further south, a common pattern among many ectotherms. We investigated differences in daily shell growth among scallop populations along the Northeast Atlantic coast from Spain to Norway. This study design allowed us to address precisely whether the asymptotic size observed along a latitudinal gradient, mainly defined by a temperature gradient, results from differences in annual or daily growth rates, or a difference in the length of the growing season. We found that low annual growth rates in northern populations are not due to low daily growth values, but to the smaller number of days available each year to achieve growth compared to the south. We documented a decrease in the annual number of growth days with age regardless of latitude. However, despite initially lower annual growth performances in terms of growing season length and growth rate, differences in asymptotic size as a function of latitude resulted from persistent annual growth performances in the north and sharp declines in the south. Our measurements of daily growth rates throughout life in a long-lived ectothermic species provide new insight into spatio-temporal variations in growth dynamics and growing season length that cannot be accounted for by classical growth models that only address asymptotic size and annual growth rate.

Adaptation to climate change: contrasting patterns of thermal-reaction-norm evolution in Pacific versus Atlantic silversides

How organisms may adapt to rising global temperatures is uncertain, but concepts can emerge from studying adaptive physiological trait variations across existing spatial climate gradients. Many ectotherms, particularly fish, have evolved increasing genetic growth capacities with latitude (i.e. countergradient variation (CnGV) in growth), which are thought to be an adaptation primarily to strong gradients in seasonality. In contrast, evolutionary responses to gradients in mean temperature are often assumed to involve an alternative mode, 'thermal adaptation'. We measured thermal growth reaction norms in Pacific silverside populations (Atherinops affinis) occurring across a weak latitudinal temperature gradient with invariant seasonality along the North American Pacific coast. Instead of thermal adaptation, we found novel evidence for CnGV in growth, suggesting that CnGV is a ubiquitous mode of reaction-norm evolution in ectotherms even in response to weak spatial and, by inference, temporal climate gradients. A novel, large-scale comparison between ecologically equivalent Pacific versus Atlantic silversides (Menidia menidia) revealed how closely growth CnGV patterns reflect their respective climate gradients. While steep growth reaction norms and increasing growth plasticity with latitude in M. menidia mimicked the strong, highly seasonal Atlantic coastal gradient, shallow reaction norms and much smaller, latitude-independent growth plasticity in A. affinis resembled the weak Pacific latitudinal temperature gradient.

Countergradient variation in body shape between two populations of Atlantic cod (Gadus morhua)

Variation in morphological traits is generally thought to be cogradient, with environmental effects on phenotypic expression reinforcing genetic differences between populations. We compared body shape between two populations of Atlantic cod (Gadus morhua). Striking shape differences occurred between juveniles from the two populations when reared in a common laboratory environment. However, no difference in body shape occurred between wild-reared juveniles from the two populations, suggesting that the genetic differences between populations were obscured by opposing effects of the environmental differences experienced in the wild. We suggest that much of the genetic diversity in body shape of fishes may be cryptic, with stabilizing selection for the same optimal phenotype resulting in genetic divergence between populations subject to contrasting environmental influences.