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.

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 million years 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.

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.