Environmental heterogeneity and balancing selection in the acorn barnacle Semibalanus balanoides

From tides to nucleotides: Genomic signatures of adaptation to environmental heterogeneity in barnacles

The northern acorn barnacle (Semibalanus balanoides) is a robust system to study the genetic basis of adaptations to highly heterogeneous environments. Adult barnacles may be exposed to highly dissimilar levels of thermal stress depending on where they settle in the intertidal (i.e., closer to the upper or lower tidal boundary). For instance, barnacles near the upper tidal limit experience episodic summer temperatures above recorded heat coma levels. This differential stress at the microhabitat level is also dependent on the aspect of sun exposure. In the present study, we used pool-seq approaches to conduct a genome wide screen for loci responding to intertidal zonation across the North Atlantic basin (Maine, Rhode Island, and Norway). Our analysis discovered 382 genomic regions containing SNPs which are consistently zonated (i.e., SNPs whose frequencies vary depending on their position in the rocky intertidal) across all surveyed habitats. Notably, most zonated SNPs are young and private to the North Atlantic. These regions show high levels of genetic differentiation across ecologically extreme microhabitats concomitant with elevated levels of genetic variation and Tajima's D, suggesting the action of non-neutral processes. Overall, these findings support the hypothesis that spatially heterogeneous selection is a general and repeatable feature for this species, and that natural selection can maintain functional genetic variation in heterogeneous environments.

Evolutionary Change in the Eastern Oyster, Crassostrea Virginica, Following Low Salinity Exposure

The presence of standing genetic variation will play a role in determining a population's capacity to adapt to environmentally relevant stressors. In the Gulf of Mexico, extreme climatic events and anthropogenic changes to local hydrology will expose productive oyster breeding grounds to stressful low salinity conditions. We identified genetic variation for performance under low salinity (due to the combined effects of low salinity and genetic load) using a single-generation selection experiment on larvae from two populations of the eastern oyster, Crassostrea virginica. We used pool-sequencing to test for allele frequency differences at 152 salinity-associated genes for larval families pre- and post-low salinity exposure. Our results have implications for how evolutionary change occurs during early life history stages at environmentally relevant salinities. Consistent with observations of high genetic load observed in oysters, we demonstrate evidence for purging of deleterious alleles at the larval stage in C. virginica. In addition, we observe increases in allele frequencies at multiple loci, suggesting that natural selection for low salinity performance at the larval stage can act as a filter for genotypes found in adult populations.

Temporally balanced selection during development of larval Pacific oysters (Crassostrea gigas) inherently preserves genetic diversity within offspring

Balancing selection is one of the mechanisms which has been proposed to explain the maintenance of genetic diversity in species across generations. For species with large populations and complex life histories, however, heterogeneous selection pressures may create a scenario in which the net effects of selection are balanced across developmental stages. With replicated cultures and a pooled sequencing approach, we show that genotype-dependent mortality in larvae of the Pacific oyster (Crassostrea gigas) is largely temporally dynamic and inconsistently in favour of a single genotype or allelic variant at each locus. Overall, the patterns of genetic change we observe to be taking place are more complex than what would be expected under classical examples of additive or dominant genetic interactions. They are also not easily explained by our current understanding of the effects of genetic load. Collectively, temporally heterogeneous selection pressures across different larval developmental stages may act to maintain genetic diversity, while also inherently sheltering genetic load within oyster populations.

Balanced polymorphism fuels rapid selection in an invasive crab despite high gene flow and low genetic diversity

Adaptation across environmental gradients has been demonstrated in numerous systems with extensive dispersal, despite high gene flow and consequently low genetic structure. The speed and mechanisms by which such adaptation occurs remain poorly resolved, but are critical to understanding species spread and persistence in a changing world. Here, we investigate these mechanisms in the European green crab Carcinus maenas, a globally distributed invader. We focus on a northwestern Pacific population that spread across >12 degrees of latitude in 10 years from a single source, following its introduction <35 years ago. Using six locations spanning >1500 km, we examine genetic structure using 9376 single nucleotide polymorphisms (SNPs). We find high connectivity among five locations, with significant structure between these locations and an enclosed lagoon with limited connectivity to the coast. Among the five highly connected locations, the only structure observed was a cline driven by a handful of SNPs strongly associated with latitude and winter temperature. These SNPs are almost exclusively found in a large cluster of genes in strong linkage disequilibrium that was previously identified as a candidate for cold tolerance adaptation in this species. This region may represent a balanced polymorphism that evolved to promote rapid adaptation in variable environments despite high gene flow, and which now contributes to successful invasion and spread in a novel environment. This research suggests an answer to the paradox of genetically depauperate yet successful invaders: populations may be able to adapt via a few variants of large effect despite low overall diversity.

Local Adaptation in Marine Foundation Species at Microgeographic Scales

AbstractNearshore foundation species in coastal and estuarine systems (e.g., salt marsh grasses, mangroves, seagrasses, corals) drive the ecological functions of ecosystems and entire biomes by creating physical structure that alters local abiotic conditions and influences species interactions and composition. The resilience of foundation species and the ecosystem functions they provide depends on their phenotypic and genetic responses to spatial and temporal shifts in environmental conditions. In this review, we explore what is known about the causes and consequences of adaptive genetic differentiation in marine foundation species over spatial scales shorter than dispersal capabilities (i.e., microgeographic scales). We describe the strength of coupling field and laboratory experiments with population genetic techniques to illuminate patterns of local adaptation, and we illustrate this approach by using several foundation species. Among the major themes that emerge from our review include (1) adaptive differentiation of marine foundation species repeatedly evolves along vertical (i.e., elevation or depth) gradients, and (2) mating system and phenology may facilitate this differentiation. Microgeographic adaptation is an understudied mechanism potentially underpinning the resilience of many sessile marine species, and this evolutionary mechanism likely has particularly important consequences for the ecosystem functions provided by foundation species.

Balancing selection maintains ancient polymorphisms at conserved enhancers for the olfactory receptor genes of a Chinese marine fish

The study of balancing selection, as a selective force maintaining adaptive genetic variation in gene pools longer than expected by drift, is currently experiencing renewed interest due to the increased availability of new data, methods of analysis, and case studies. In this investigation, evidence of balancing selection operating on conserved enhancers of the olfactory receptor (OR) genes is presented for the Chinese sleeper (Bostrychus sinensis), a coastal marine fish that is emerging as a model species for evolutionary studies in the Northwest Pacific marginal seas. Coupled with tests for Gene Ontology enrichment and transcription factor binding, population genomic data allow for the identification of an OR cluster in the sleeper with a downstream flanking region containing three enhancers that are conserved with human and other fish species. Phylogenetic and population genetic analyses indicate that the enhancers are under balancing selection as evidenced by their translineage polymorphisms, excess common alleles, and increased within-group diversities. Age comparisons between the translineage polymorphisms and most recent common ancestors of neutral genealogies substantiate that the former are old, and thus, due to ancient balancing selection. The survival and reproduction of vertebrates depend on their sense of smell, and thereby, on their ORs. In addition to locus duplication and allelic variation of structural genes, this study highlights a third mechanism by which receptor diversity can be achieved for detecting and responding to the huge variety of environmental odorants (i.e., by balancing selection acting on OR gene expression through their enhancer variability).

Ecological Load and Balancing Selection in Circumboreal Barnacles

Acorn barnacle adults experience environmental heterogeneity at various spatial scales of their circumboreal habitat, raising the question of how adaptation to high environmental variability is maintained in the face of strong juvenile dispersal and mortality. Here, we show that 4% of genes in the barnacle genome experience balancing selection across the entire range of the species. Many of these genes harbor mutations maintained across 2 My of evolution between the Pacific and Atlantic oceans. These genes are involved in ion regulation, pain reception, and heat tolerance, functions which are essential in highly variable ecosystems. The data also reveal complex population structure within and between basins, driven by the trans-Arctic interchange and the last glaciation. Divergence between Atlantic and Pacific populations is high, foreshadowing the onset of allopatric speciation, and suggesting that balancing selection is strong enough to maintain functional variation for millions of years in the face of complex demography.

Interactive effects of large- and local-scale environmental gradients on phenotypic differentiation

Intraspecific differentiation across a steep environmental gradient depends on the relative influences of evolutionary, organismal, and environmental processes. But steep environmental gradients may be nested within larger-scale, regional conditions that could influence these processes at the local scale. Therefore, we hypothesized that phenotypic differentiation along a steep environmental gradient would vary among regions. To test this hypothesis, we conducted a reciprocal transplant experiment on rocky intertidal shores, a habitat characterized by gradients in abiotic and biotic stress, in three regions of the Gulf of Maine. We used the ubiquitous and ecologically important rockweed species Fucus vesiculosus to quantify differentiation in growth, tissue nitrogen, and nitrogen productivity between upper and lower intertidal individuals. We found that phenotypic differentiation between tide heights varied among traits and regions. Although tissue nitrogen did not vary among any treatment combinations, growth and nitrogen productivity response were region specific. A strong effect of transplant height was found in all regions; however, an effect of home (source) height was only detectable in the central Gulf of Maine. Our study reveals that intraspecific responses to steep environmental gradients vary among populations, but the mechanisms underlying these patterns remain unknown. Given the roles that rockweeds play as food and habitat, these in situ patterns of growth and nitrogen productivity could have important community- and ecosystem-level consequences.

Effects of balancing selection and microhabitat temperature variations on heat tolerance of the intertidal black mussel Septifer virgatus

Realistic assessments of the impacts of global warming on population extinction risk are likely to require an integrated analysis of the roles of standing genetic variation, microhabitat thermal complexity, and the inter-individual variation of heat tolerance due to both genetic differences and seasonal acclimatization effects. Here, we examine whether balancing selection and microhabitat temperature heterogeneity can interact to enhance the population persistence to thermal stress for the black mussel Septifer virgatus. We deployed biomimetic data loggers on the shore to measure the microhabitat-specific thermal variation from June 2014 to April 2016. Thermal tolerance of specimens was indexed by measuring effects of temperature on heart rate. Genotyping of specimens was performed using double digestion restriction association RADSeq (ddRADseq). Our results show that inter-individual variations in thermal tolerance correlate significantly with genetic differences at some specific gene loci, and that heterozygotes have higher thermal tolerances than homozygotes. The observed seasonal changes in genotype frequency suggest that these loci are under balancing selection. The ability of thermally resistant heterozygotes to survive in sun-exposed microhabitats acts to balance the loss of homozygotes during summer and enable the persistence of genetic polymorphisms. Population persistence of the mussel is also facilitated by the micro-scale variation in temperature, which provides refugia from thermal stress. Our results emphasize that inter-individual variation in thermal tolerance and in microhabitat heterogeneity in temperature are important for the persistence of populations in rocky shore habitats.

Footprints of natural selection at the mannose-6-phosphate isomerase locus in barnacles

The rocky intertidal is a natural laboratory to study how natural selection acts on the genes and proteins responsible for organismal survival and reproduction. Alternative forms of enzymes that differ across the intertidal have been known for decades and have provided examples of selection, but the genetic basis of such enzyme variation is known in only a few cases. In this paper, we present molecular evidence of natural selection at the Mpi gene, a key enzyme in energy metabolism that alters survival of barnacles living across the stress gradient imposed by the intertidal. Our study demonstrates how natural selection can facilitate survival in highly heterogeneous environments through the maintenance of multiple molecular solutions to ecological stresses.

The mannose-6-phosphate isomerase (Mpi) locus in Semibalanus balanoides has been studied as a candidate gene for balancing selection for more than two decades. Previous work has shown that Mpi allozyme genotypes (fast and slow) have different frequencies across Atlantic intertidal zones due to selection on postsettlement survival (i.e., allele zonation). We present the complete gene sequence of the Mpi locus and quantify nucleotide polymorphism in S. balanoides, as well as divergence to its sister taxon Semibalanus cariosus. We show that the slow allozyme contains a derived charge-altering amino acid polymorphism, and both allozyme classes correspond to two haplogroups with multiple internal haplotypes. The locus shows several footprints of balancing selection around the fast/slow site: an enrichment of positive Tajima’s D for nonsynonymous mutations, an excess of polymorphism, and a spike in the levels of silent polymorphism relative to silent divergence, as well as a site frequency spectrum enriched for midfrequency mutations. We observe other departures from neutrality across the locus in both coding and noncoding regions. These include a nonsynonymous trans-species polymorphism and a recent mutation under selection within the fast haplogroup. The latter suggests ongoing allelic replacement of functionally relevant amino acid variants. Moreover, predicted models of Mpi protein structure provide insight into the functional significance of the putatively selected amino acid polymorphisms. While footprints of selection are widespread across the range of S. balanoides, our data show that intertidal zonation patterns are variable across both spatial and temporal scales. These data provide further evidence for heterogeneous selection on Mpi.

Adaptive landscape genetics and malaria across divergent island bird populations

Environmental conditions play a major role in shaping the spatial distributions of pathogens, which in turn can drive local adaptation and divergence in host genetic diversity. Haemosporidians, such as Plasmodium (malaria), are a strong selective force, impacting survival and fitness of hosts, with geographic distributions largely determined by habitat suitability for their insect vectors. Here, we have tested whether patterns of fine-scale local adaptation to malaria are replicated across discrete, ecologically differing island populations of Berthelot's pipits Anthus berthelotii. We sequenced TLR4, an innate immunity gene that is potentially under positive selection in Berthelot's pipits, and two SNPs previously identified as being associated with malaria infection in a genome-wide association study (GWAS) in Berthelot's pipits in the Canary Islands. We determined the environmental predictors of malaria infection, using these to estimate variation in malaria risk on Porto Santo, and found some congruence with previously identified environmental risk factors on Tenerife. We also found a negative association between malaria infection and a TLR4 variant in Tenerife. In contrast, one of the GWAS SNPs showed an association with malaria risk in Porto Santo, but in the opposite direction to that found in the Canary Islands GWAS. Together, these findings suggest that disease-driven local adaptation may be an important factor in shaping variation among island populations.

Mapping physiology: biophysical mechanisms define scales of climate change impacts

The rocky intertidal zone is a highly dynamic and thermally variable ecosystem, where the combined influences of solar radiation, air temperature and topography can lead to differences greater than 15°C over the scale of centimetres during aerial exposure at low tide. For most intertidal organisms this small-scale heterogeneity in microclimates can have enormous influences on survival and physiological performance. However, the potential ecological importance of environmental heterogeneity in determining ecological responses to climate change remains poorly understood. We present a novel framework for generating spatially explicit models of microclimate heterogeneity and patterns of thermal physiology among interacting organisms. We used drone photogrammetry to create a topographic map (digital elevation model) at a resolution of 2 × 2 cm from an intertidal site in Massachusetts, which was then fed into to a model of incident solar radiation based on sky view factor and solar position. These data were in turn used to drive a heat budget model that estimated hourly surface temperatures over the course of a year (2017). Body temperature layers were then converted to thermal performance layers for organisms, using thermal performance curves, creating 'physiological landscapes' that display spatially and temporally explicit patterns of 'microrefugia'. Our framework shows how non-linear interactions between these layers lead to predictions about organismal performance and survivorship that are distinct from those made using any individual layer (e.g. topography, temperature) alone. We propose a new metric for quantifying the 'thermal roughness' of a site (RqT, the root mean square of spatial deviations in temperature), which can be used to quantify spatial and temporal variability in temperature and performance at the site level. These methods facilitate an exploration of the role of micro-topographic variability in driving organismal vulnerability to environmental change using both spatially explicit and frequency-based approaches.

Negative Frequency-Dependent Selection Is Frequently Confounding

Persistent genetic variation within populations presents an evolutionary problem, as natural selection and genetic drift tend to erode genetic diversity. Models of balancing selection were developed to account for the maintenance of genetic variation observed in natural populations. Negative frequency-dependent selection is a powerful type of balancing selection that maintains many natural polymorphisms, but it is also commonly misinterpreted. This review aims to clarify the processes underlying negative frequency-dependent selection, describe classes of polymorphisms that can and cannot result from these processes, and discuss the empirical data needed to accurately identify processes that generate or maintain diversity in nature. Finally, the importance of accurately describing the processes affecting genetic diversity within populations as it relates to research progress is considered.

Balancing selection maintains diversity in a cold tolerance gene in broadly distributed live oaks

Cold poses major physiological challenges to plants, especially long-lived trees. In trees occurring along variable temperature clines, the expected direction and consequences of selection on cold acclimation ability and freezing tolerance are not straightforward. Here we estimated selection in cold acclimation genes at two evolutionary timescales in all seven species of the American live oaks (Quercus subsection Virentes). Two cold response candidate genes were chosen: ICE1, a key gene in the cold acclimation pathway, and HOS1, which modulates cold response by negatively regulating ICE1. Two housekeeping genes, GAPDB and CHR11, were also analyzed. At the shallow evolutionary timescale, we demonstrate that HOS1 experienced recent balancing selection in the two most broadly distributed species, Q. virginiana and Q. oleoides. At a deeper evolutionary scale, a codon-based model of evolution revealed the signature of negative selection in ICE1. In contrast, three positively selected codons have been identified in HOS1, possibly a signature of the diversification of Virentes into warmer climates from a freezing adapted lineage of oaks. Our findings indicate that evolution has favored diversity in cold tolerance modulation through balancing selection in HOS1 while maintaining core cold acclimation ability, as evidenced by purifying selection in ICE1.

Thermal variation, thermal extremes and the physiological performance of individuals

In this review we consider how small-scale temporal and spatial variation in body temperature, and biochemical/physiological variation among individuals, affect the prediction of organisms' performance in nature. For ‘normal’ body temperatures – benign temperatures near the species' mean – thermal biology traditionally uses performance curves to describe how physiological capabilities vary with temperature. However, these curves, which are typically measured under static laboratory conditions, can yield incomplete or inaccurate predictions of how organisms respond to natural patterns of temperature variation. For example, scale transition theory predicts that, in a variable environment, peak average performance is lower and occurs at a lower mean temperature than the peak of statically measured performance. We also demonstrate that temporal variation in performance is minimized near this new ‘optimal’ temperature. These factors add complexity to predictions of the consequences of climate change. We then move beyond the performance curve approach to consider the effects of rare, extreme temperatures. A statistical procedure (the environmental bootstrap) allows for long-term simulations that capture the temporal pattern of extremes (a Poisson interval distribution), which is characterized by clusters of events interspersed with long intervals of benign conditions. The bootstrap can be combined with biophysical models to incorporate temporal, spatial and physiological variation into evolutionary models of thermal tolerance. We conclude with several challenges that must be overcome to more fully develop our understanding of thermal performance in the context of a changing climate by explicitly considering different forms of small-scale variation. These challenges highlight the need to empirically and rigorously test existing theories.

Worldwide genetic differentiation in the common fouling barnacle, Amphibalanus amphitrite

Amphibalanus amphitrite is a common fouling barnacle distributed globally in tropical and subtropical waters. In the present study, the genetic (mitochondrial cytochrome oxidase subunit I) and morphological differentiation in A. amphitrite from 25 localities around the world were investigated. The results revealed three clades within A. amphitrite with a genetic divergence of ~ 4% among clades, whereas there were no diagnostic morphological differences among clades. Clade 1 is widely distributed in both temperate and tropical waters, whereas Clade 3 is currently restricted to the tropical region. The deep divergence among clades suggests historical isolation within A. amphitrite; thus, the present geographical overlaps are possibly a result of the combined effects of rising sea level and human-mediated dispersals. This study highlights the genetic differentiation that exists in a common, widely distributed fouling organism with great dispersal potential; future antifouling research should take into account the choice of lineages.

Repeated landmass reformation limits diversification in the widespread littoral zone mosquito Anopheles sundaicus sensu lato in the Indo-Oriental Region

Southeast Asia harbours abundant biodiversity, hypothesized to have been generated by Pliocene and Pleistocene climatic and environmental change. Vicariance between the island of Borneo, the remaining Indonesian archipelago and mainland Southeast Asia caused by elevated sea levels during interglacial periods has been proposed to lead to diversification in the littoral zone mosquito Anopheles (Cellia) sundaicus (Rodenwaldt) sensu lato. To test this biogeographical hypothesis, we inferred the population history and assessed gene flow of A. sundaicus s.l. sampled from 18 populations across its pan-Asian species range, using sequences from mitochondrial cytochrome c oxidase subunit 1 (CO1), the internal transcribed spacer 2 (ITS2) and the mannose phosphate isomerase (Mpi) gene. A hypothesis of ecological speciation for A. sundaicus involving divergent adaptation to brackish and freshwater larval habitats was also previously proposed, based on a deficiency of heterozygotes for Mpi allozyme alleles in sympatry. This hypothesis was not supported by Mpi sequence data, which exhibited no fixed differences between brackish and freshwater larval habitats. Mpi and CO1 supported the presence of up to eight genetically distinct population groupings. Counter to the hypothesis of three allopatric species, divergence was often no greater between Borneo, Sumatra/Java and the Southeast Asian mainland than it was between genetic groupings within these landmasses. An isolation-with-migration (IM) model indicates recurrent gene flow between the current major landmasses. Such gene flow would have been possible during glacial periods when the current landmasses merged, presenting opportunities for dispersal along expanding and contracting coastlines. Consequently, Pleistocene climatic variation has proved a homogenizing, rather than diversifying, force for A. sundaicus diversity.

Importance of plasticity and local adaptation for coping with changing salinity in coastal areas: a test case with barnacles in the Baltic Sea

Background

Salinity plays an important role in shaping coastal marine communities. Near-future climate predictions indicate that salinity will decrease in many shallow coastal areas due to increased precipitation; however, few studies have addressed this issue. The ability of ecosystems to cope with future changes will depend on species’ capacities to acclimatise or adapt to new environmental conditions. Here, we investigated the effects of a strong salinity gradient (the Baltic Sea system – Baltic, Kattegat, Skagerrak) on plasticity and adaptations in the euryhaline barnacle Balanus improvisus. We used a common-garden approach, where multiple batches of newly settled barnacles from each of three different geographical areas along the Skagerrak-Baltic salinity gradient were exposed to corresponding native salinities (6, 15 and 30 PSU), and phenotypic traits including mortality, growth, shell strength, condition index and reproductive maturity were recorded.

Results

We found that B. improvisus was highly euryhaline, but had highest growth and reproductive maturity at intermediate salinities. We also found that low salinity had negative effects on other fitness-related traits including initial growth and shell strength, although mortality was also lowest in low salinity. Overall, differences between populations in most measured traits were weak, indicating little local adaptation to salinity. Nonetheless, we observed some population-specific responses – notably that populations from high salinity grew stronger shells in their native salinity compared to the other populations, possibly indicating adaptation to differences in local predation pressure.

Conclusions

Our study shows that B. improvisus is an example of a true brackish-water species, and that plastic responses are more likely than evolutionary tracking in coping with future changes in coastal salinity.

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.

Adaptation and acclimatization to ocean acidification in marine ectotherms: an in situ transplant experiment with polychaetes at a shallow CO2 vent system

Metabolic rate determines the physiological and life-history performances of ectotherms. Thus, the extent to which such rates are sensitive and plastic to environmental perturbation is central to an organism's ability to function in a changing environment. Little is known of long-term metabolic plasticity and potential for metabolic adaptation in marine ectotherms exposed to elevated pCO2. Consequently, we carried out a series of in situ transplant experiments using a number of tolerant and sensitive polychaete species living around a natural CO2 vent system. Here, we show that a marine metazoan (i.e. Platynereis dumerilii) was able to adapt to chronic and elevated levels of pCO2. The vent population of P. dumerilii was physiologically and genetically different from nearby populations that experience low pCO2, as well as smaller in body size. By contrast, different populations of Amphiglena mediterranea showed marked physiological plasticity indicating that adaptation or acclimatization are both viable strategies for the successful colonization of elevated pCO2 environments. In addition, sensitive species showed either a reduced or increased metabolism when exposed acutely to elevated pCO2. Our findings may help explain, from a metabolic perspective, the occurrence of past mass extinction, as well as shed light on alternative pathways of resilience in species facing ongoing ocean acidification.

Host-specific phenotypic plasticity of the turtle barnacle Chelonibia testudinaria: a widespread generalist rather than a specialist

Turtle barnacles are common epibionts on marine organisms. Chelonibia testudinaria is specific on marine turtles whereas C. patula is a host generalist, but rarely found on turtles. It has been questioned why C. patula, being abundant on a variety of live substrata, is almost absent from turtles. We evaluated the genetic (mitochondrial COI, 16S and 12S rRNA, and amplified fragment length polymorphism (AFLP)) and morphological differentiation of C. testudinaia and C. patula from different hosts, to determine the mode of adaptation exhibited by Chelonibia species on different hosts. The two taxa demonstrate clear differences in shell morphology and length of 4-6(th) cirri, but very similar in arthropodal characters. Moreover, we detected no genetic differentiation in mitochondrial DNA and AFLP analyses. Outlier detection infers insignificant selection across loci investigated. Based on combined morphological and molecular evidence, we proposed that C. testudinaria and C. patula are conspecific, and the two morphs with contrasting shell morphologies and cirral length found on different host are predominantly shaped by developmental plasticity in response to environmental setting on different hosts. Chelonibia testudinaria is, thus, a successful general epibiotic fouler and the phenotypic responses postulated can increase the fitness of the animals when they attach on hosts with contrasting life-styles.

Longitudinal trends in climate drive flowering time clines in North American Arabidopsis thaliana

Introduced species frequently show geographic differentiation, and when differentiation mirrors the ancestral range, it is often taken as evidence of adaptive evolution. The mouse-ear cress (Arabidopsis thaliana) was introduced to North America from Eurasia 150–200 years ago, providing an opportunity to study parallel adaptation in a genetic model organism. Here, we test for clinal variation in flowering time using 199 North American (NA) accessions of A. thaliana, and evaluate the contributions of major flowering time genes FRI, FLC, and PHYC as well as potential ecological mechanisms underlying differentiation. We find evidence for substantial within population genetic variation in quantitative traits and flowering time, and putatively adaptive longitudinal differentiation, despite low levels of variation at FRI, FLC, and PHYC and genome-wide reductions in population structure relative to Eurasian (EA) samples. The observed longitudinal cline in flowering time in North America is parallel to an EA cline, robust to the effects of population structure, and associated with geographic variation in winter precipitation and temperature. We detected major effects of FRI on quantitative traits associated with reproductive fitness, although the haplotype associated with higher fitness remains rare in North America. Collectively, our results suggest the evolution of parallel flowering time clines through novel genetic mechanisms.

Genetic variation in the acorn barnacle from allozymes to population genomics

Understanding the patterns of genetic variation within and among populations is a central problem in population and evolutionary genetics. We examine this question in the acorn barnacle, Semibalanus balanoides, in which the allozyme loci Mpi and Gpi have been implicated in balancing selection due to varying selective pressures at different spatial scales. We review the patterns of genetic variation at the Mpi locus, compare this to levels of population differentiation at mtDNA and microsatellites, and place these data in the context of genome-wide variation from high-throughput sequencing of population samples spanning the North Atlantic. Despite considerable geographic variation in the patterns of selection at the Mpi allozyme, this locus shows rather low levels of population differentiation at ecological and trans-oceanic scales (F(ST) ~ 5%). Pooled population sequencing was performed on samples from Rhode Island (RI), Maine (ME), and Southwold, England (UK). Analysis of more than 650 million reads identified approximately 335,000 high-quality SNPs in 19 million base pairs of the S. balanoides genome. Much variation is shared across the Atlantic, but there are significant examples of strong population differentiation among samples from RI, ME, and UK. An F(ST) outlier screen of more than 22,000 contigs provided a genome-wide context for interpretation of earlier studies on allozymes, mtDNA, and microsatellites. F(ST) values for allozymes, mtDNA and microsatellites are close to the genome-wide average for random SNPs, with the exception of the trans-Atlantic F(ST) for mtDNA. The majority of F(ST) outliers were unique between individual pairs of populations, but some genes show shared patterns of excess differentiation. These data indicate that gene flow is high, that selection is strong on a subset of genes, and that a variety of genes are experiencing diversifying selection at large spatial scales. This survey of polymorphism in S. balanoides provides a number of genomic tools that promise to make this a powerful model for ecological genomics of the rocky intertidal.

Inference on the strength of balancing selection for epistatically interacting loci

Existing inference methods for estimating the strength of balancing selection in multi-locus genotypes rely on the assumption that there are no epistatic interactions between loci. Complex systems in which balancing selection is prevalent, such as sets of human immune system genes, are known to contain components that interact epistatically. Therefore, current methods may not produce reliable inference on the strength of selection at these loci. In this paper, we address this problem by presenting statistical methods that can account for epistatic interactions in making inference about balancing selection. A theoretical result due to Fearnhead (2006) is used to build a multi-locus Wright-Fisher model of balancing selection, allowing for epistatic interactions among loci. Antagonistic and synergistic types of interactions are examined. The joint posterior distribution of the selection and mutation parameters is sampled by Markov chain Monte Carlo methods, and the plausibility of models is assessed via Bayes factors. As a component of the inference process, an algorithm to generate multi-locus allele frequencies under balancing selection models with epistasis is also presented. Recent evidence on interactions among a set of human immune system genes is introduced as a motivating biological system for the epistatic model, and data on these genes are used to demonstrate the methods.

Local adaptation in marine invertebrates

Local adaptation in the sea was regarded historically as a rare phenomenon that was limited to a handful of species with exceptionally low dispersal potential. However, a growing body of experimental studies indicates that adaptive differentiation occurs in numerous marine invertebrates in response to selection imposed by strong gradients (and more complex mosaics) of abiotic and biotic conditions. Moreover, a surprisingly high proportion of the marine invertebrates known or suspected of exhibiting local adaptation are species with planktonic dispersal. Adaptive divergence among populations can occur over a range of spatial scales, including those that are fine-grained (i.e., meters to kilometers), reflecting a balance between scales of gene flow and selection. Addressing the causes and consequences of adaptive genetic differentiation among invertebrate populations promises to advance community ecology, climate change research, and the effective management of marine ecosystems.

Organismal climatology: analyzing environmental variability at scales relevant to physiological stress

Predicting when, where and with what magnitude climate change is likely to affect the fitness, abundance and distribution of organisms and the functioning of ecosystems has emerged as a high priority for scientists and resource managers. However, even in cases where we have detailed knowledge of current species' range boundaries, we often do not understand what, if any, aspects of weather and climate act to set these limits. This shortcoming significantly curtails our capacity to predict potential future range shifts in response to climate change, especially since the factors that set range boundaries under those novel conditions may be different from those that set limits today. We quantitatively examine a nine-year time series of temperature records relevant to the body temperatures of intertidal mussels as measured using biomimetic sensors. Specifically, we explore how a 'climatology' of body temperatures, as opposed to long-term records of habitat-level parameters such as air and water temperatures, can be used to extrapolate meaningful spatial and temporal patterns of physiological stress. Using different metrics that correspond to various aspects of physiological stress (seasonal means, cumulative temperature and the return time of extremes) we show that these potential environmental stressors do not always occur in synchrony with one another. Our analysis also shows that patterns of animal temperature are not well correlated with simple, commonly used metrics such as air temperature. Detailed physiological studies can provide guidance to predicting the effects of global climate change on natural ecosystems but only if we concomitantly record, archive and model environmental signals at appropriate scales.

From cells to coastlines: how can we use physiology to forecast the impacts of climate change?

The interdisciplinary fields of conservation physiology, macrophysiology, and mechanistic ecological forecasting have recently emerged as means of integrating detailed physiological responses to the broader questions of ecological and evolutionary responses to global climate change. Bridging the gap between large-scale records of weather and climate (as measured by remote sensing platforms, buoys and ground-based weather stations) and the physical world as experienced by organisms (niche-level measurements) requires a mechanistic understanding of how ;environmental signals' (parameters such as air, surface and water temperature, food availability, water flow) are translated into signals at the scale of the organism or cell (e.g. body temperature, food capture, hydrodynamic force, aerobic capacity). Predicting the impacts of how changing environments affect populations and ecosystems further mandates an understanding of how organisms ;filter' these signals via their physiological response (e.g. whether they respond to high or low frequencies, whether there is a time lag in response, etc.) and must be placed within the context of adult movement and the dispersal of larvae and gametes. Recent studies have shown that patterns of physiological stress in nature are far more complex in space and time than previously assumed and challenge the long-held paradigm that patterns of biogeographic distribution can be based on simple environmental gradients. An integrative, systems-based approach can provide an understanding of the roles of environmental and physiological variability in driving ecological responses and can offer considerable insight and predictive capacity to researchers, resource managers and policy makers involved in planning for the current and future effects of climate change.

Geographic distributions of Idh-1 alleles in a cricket are linked to differential enzyme kinetic performance across thermal environments

Background

Geographic clines within species are often interpreted as evidence of adaptation to varying environmental conditions. However, clines can also result from genetic drift, and these competing hypotheses must therefore be tested empirically. The striped ground cricket, Allonemobius socius, is widely-distributed in the eastern United States, and clines have been documented in both life-history traits and genetic alleles. One clinally-distributed locus, isocitrate dehydrogenase (Idh-1), has been shown previously to exhibit significant correlations between allele frequencies and environmental conditions (temperature and rainfall). Further, an empirical study revealed a significant genotype-by-environmental interaction (GxE) between Idh-1 genotype and temperature which affected fitness. Here, we use enzyme kinetics to further explore GxE between Idh-1 genotype and temperature, and test the predictions of kinetic activity expected under drift or selection.

Results

We found significant GxE between temperature and three enzyme kinetic parameters, providing further evidence that the natural distributions of Idh-1 allele frequencies in A. socius are maintained by natural selection. Differences in enzyme kinetic activity across temperatures also mirror many of the geographic patterns observed in allele frequencies.

Conclusion

This study further supports the hypothesis that the natural distribution of Idh-1 alleles in A. socius is driven by natural selection on differential enzymatic performance. This example is one of several which clearly document a functional basis for both the maintenance of common alleles and observed clines in allele frequencies, and provides further evidence for the non-neutrality of some allozyme alleles.

Ecological genetics in the North Atlantic: environmental gradients and adaptation at specific loci

The North Atlantic intertidal community provides a rich set of organismal and environmental material for the study of ecological genetics. Clearly defined environmental gradients exist at multiple spatial scales: there are broad latitudinal trends in temperature, meso-scale changes in salinity along estuaries, and smaller scale gradients in desiccation and temperature spanning the intertidal range. The geology and geography of the American and European coasts provide natural replication of these gradients, allowing for population genetic analyses of parallel adaptation to environmental stress and heterogeneity. Statistical methods have been developed that provide genomic neutrality tests of population differentiation and aid in the process of candidate gene identification. In this paper, we review studies of marine organisms that illustrate associations between an environmental gradient and specific genetic markers. Such highly differentiated markers become candidate genes for adaptation to the environmental factors in question, but the functional significance of genetic variants must be comprehensively evaluated. We present a set of predictions about locus-specific selection across latitudinal, estuarine, and intertidal gradients that are likely to exist in the North Atlantic. We further present new data and analyses that support and contradict these simple selection models. Some taxa show pronounced clinal variation at certain loci against a background of mild clinal variation at many loci. These cases illustrate the procedures necessary for distinguishing selection driven by internal genomic vs. external environmental factors. We suggest that the North Atlantic intertidal community provides a model system for identifying genes that matter in ecology due to the clarity of the environmental stresses and an extensive experimental literature on ecological function. While these organisms are typically poor genetic and genomic models, advances in comparative genomics have provided access to molecular tools that can now be applied to taxa with well-defined ecologies. As many of the organisms we discuss have tight physiological limits driven by climatic factors, this synthesis of molecular population genetics with marine ecology could provide a sensitive means of assessing evolutionary responses to climate change.

Highly divergent duplicate mannose-6-phosphate isomerase (Mpi) genes in the blue mussel, Mytilus edulis

Gene duplication has been hypothesized to play a major role in the evolution of genes and genomes and in generating phenotypic diversity among the proteins those genes encode. We have identified duplicate genes for the glycolytic enzyme mannose-6-phosphate isomerase (MPI: EC 5.3.1.8) in marine mussels in the genus Mytilus. Overall, there was only 52% sequence identity (72% similarity) between the proteins encoded by these two genes, designated as Mpi-A and Mpi-B. Based on a comparison of the rate of non-synonymous substitution between orthologous and paralogous Mpi coding sequences obtained from Mytilus edulis and the congener M. trossulus we estimate that the duplication of Mpi in mussels occurred ~170 MYA. We detected paralog-specific differences in the ratio of non-synonymous to synonymous substitutions (k(a)/k(s)) and in the predicted net charge for MPI-A and MPI-B. Using a real-time quantitative RT-PCR assay we observed substantial changes in Mpi-A and Mpi-B transcript levels between tissue types; the strongest expression of Mpi-A was observed in mantle tissue while Mpi-B expression exceeded that of Mpi-A in gill and hepatopancreas tissues. Taken together, these observations suggest that different functional roles have evolved for these two Mytilus Mpi genes subsequent to gene duplication.

Genetic heterogeneity among intertidal habitats in the flat periwinkle, Littorina obtusata

Comparisons among patterns exhibited by functionally distinct genetic markers have been widely used to infer the impacts of demography and selection in structuring genetic variation in natural populations. However, such multilocus comparisons remain an indirect evaluation of selection at particular candidate loci; ideally, the identification of a candidate gene by comparative genetic methodologies should be complemented by functional analyses and experimental manipulations of genotypes in the laboratory or field. We examined genotype frequency variation among replicated intertidal habitats at two spatial scales in the grazing snail Littorina obtusata. Both of the candidate allozyme markers varied predictably with environment, and these patterns were consistent at both spatial scales. Three of four reference loci were spatially homogeneous, but one microsatellite exhibited significant structure at both geographical and mesoscales. To initiate a direct examination of whether the observed genotype frequency variation at one of the candidate markers, mannose-6-phosphate isomerase (MPI), was impacted by differential survivorship of genotypes, we conducted a series of laboratory-based thermal stress assays using snails from two geographically disparate source populations. When snails were exposed to bouts of thermal/desiccation stress, patterns of mortality were nonrandom with respect to MPI genotype. Furthermore, patterns of mortality in the laboratory manipulation coincided with the observed distribution of genotypes in the field. The data suggest the operation of selection at the Mpi or a linked locus, but functional studies and further experimentation are required to establish the relationship between MPI genotype and fitness across heterogeneous intertidal environments.

Single-locus polymorphism in a heterogeneous two-deme model

Environmental heterogeneity has long been considered a likely explanation for the high levels of genetic variation found in most natural populations: selection in a spatially heterogeneous environment can maintain more variation. While this theoretical result has been extensively studied in models with limited parameters (e.g., two alleles, fixed gene flow, and particular selection schemes), the effect of spatial heterogeneity is poorly understood for models with a wider range of parameters (e.g., multiple alleles, different levels of gene flow, and more general selection schemes). We have compared the volume of fitness space that maintains variation in a single-deme model to the volume in a two-deme model for multiple alleles, random selection schemes, and various levels of migration. Furthermore, equilibrium allele-frequency vectors were examined to see if particular patterns of variation are more prevalent than first expected. The two-deme model maintains variation for substantially larger volumes of fitness space with lower heterozygote fitness than the single-deme model. This result implies that selection schemes in the two-deme model can have a wider range of fitness patterns while still maintaining variation. The equilibrium allele-frequency patterns emerging from the two-deme model are more variable and strongly influenced by gene flow.

GENETIC IDENTITY DETERMINES RISK OF POST-SETTLEMENT MORTALITY OF A MARINE FISH

Longitudinal sampling of four cohorts of Neopomacentrus filamentosus, a common tropical damselfish from Dampier Archipelago, Western Australia, revealed the evolution of size structure after settlement. Light traps collected premetamorphic individuals from the water column ("settlers") to establish a baseline for each cohort. Subsequently, divers collected benthic juveniles ("recruits") at 1-3-month intervals to determine the relative impacts of post-settlement mortality during the first three months. Growth trajectories for individual fish were back-calculated from otolith records and compared with nonlinear mixed-effects models. Size-selective mortality was detected in all cohorts with the loss of smaller, slower growing individuals. Three months after settlement, recruits showed significantly faster growth as juveniles, faster growth as larvae, and larger sizes as hatchlings. The timing and intensity of post-settlement selection differed among cohorts and was correlated with density at settlement. The cohort with the greatest initial abundance experienced the strongest selective mortality, with most of this mortality occurring between one and two months after settlement when juveniles began foraging at higher positions in the water column. Significant genetic structure was found between settlers and three-month-old recruits in this cohort as a result of natural selection that changed the frequency of mtDNA haplotypes measured at the control region. The extent of this genetic difference was enlarged or reduced by artificially manipulating the intensity of size-based selection, thus establishing a link between phenotype and haplotype. Sequence variation in the control region of the mitochondrial genome has been linked to mitochondrial efficiency and weight gain in other studies, which provides a plausible explanation for the patterns observed here.

Adaptive phenotypic differentiation across the intertidal gradient in the alga Silvetia compressa

Populations of intertidal species span a steep environmental gradient driven by differences in emersion time. In spite of strong differential selection on traits related to this gradient, the small spatial scale over which differences occur may prevent local adaptation, and instead may favor a single intermediate phenotype, or nongenetic mechanisms of differentiation. Here I examine whether a common macroalga, Silvetia compressa, exhibits phenotypic differentiation across the intertidal gradient and evaluate how local adaptation, developmental plasticity, and maternal effects may interact to shape individual phenotypes. Reciprocal transplants of both adults and embryos showed a "home-height advantage" in two of the three populations tested. In laboratory trials, the progeny of upper-limit individuals survived exposure to air significantly better than lower-limit progeny from the same population. I compared the emersion tolerance of full-sib families generated from gametes produced in the field to those produced under common garden conditions. The relative advantage of upper-limit lineages was robust to maternal environment during gametogenesis; this pattern is consistent with genetic differentiation. The possible role of local adaptation has historically been ignored in studies of intertidal zonation. In S. compressa, phenotypic differentiation may have important consequences for vertical range, both within and among sites.