Transcriptional dynamics following freezing stress reveal selection for mechanisms of freeze tolerance at the poleward range margin in the cold water intertidal barnacle Semibalanus balanoides

The ability to survive freezing has repeatedly evolved across multiple phyla. This suggests that the mechanisms of freeze tolerance must be readily evolvable from basal physiological traits. While several biochemical correlates to freeze tolerance have been described, the mechanism that confers freeze tolerance is still not well understood. To understand both the basic biochemical mechanisms of freeze tolerance as well as their role in local adaptation at the poleward range edge, we conducted a transcriptomic study on two populations (one from the poleward range margin in the White Sea, Russia, and one from the central coast of British Columbia, Canada) of the cold water acorn barnacle Semibalanus balanoides on a time series following a freezing event. We found that the British Columbia population (at the equatorward range margin) was significantly less freeze tolerant than the White Sea population (at the poleward range margin). After assembling and annotating a de novo transcriptome for S. balanoides, we found that the patterns of differential transcript expression following freezing were almost entirely non-overlapping between the two populations, with the White Sea population expressing a series of heat shock proteins in response to freezing stress as well as several aquaporins, while the British Columbia population expressed a series of proteases instead, indicating severe protein damage. We found strong evidence of purifying selection on the significantly upregulated transcripts in the White Sea population, suggesting local adaptation to freezing threat. Taken together, this shows the importance of freeze tolerance to population survival at the poleward range margin, and highlights the central roles of aquaporins and heat shock proteins to the trait of freeze tolerance across taxa.