Genome-Wide Adaptation to a Complex Environmental Gradient in a Keystone Phytoplankton Species

Marine phytoplankton play essential roles in global primary production and biogeochemical cycles. Yet, the evolutionary genetic underpinnings of phytoplankton adaptation to complex marine and coastal environments, where many environmental variables fluctuate and interact, remain unclear. We combined population genomics with experimental transcriptomics to investigate the genomic basis underlying a natural evolutionary experiment that has played out over the past 8000 years in one of the world's largest brackish water bodies: the colonisation of the Baltic Sea by the ancestrally marine diatom Skeletonema marinoi. To this end, we combined target capture of the entire nuclear genome with pooled shotgun sequencing, and showed that the method performs well on both cultures and single cells. Genotype-environment association analyses identified > 1000 genes with signals of selection in response to major environmental gradients in the Baltic Sea, which apart from salinity, include marked differences in temperature and nutrient supply. Locally adapted genes were related to diverse metabolic processes, including signal transduction, cell cycle, DNA methylation and maintenance of homeostasis. The locally adapted genes showed significant overlap with salinity-responsive genes identified in a laboratory common garden experiment, suggesting the Baltic salinity gradient contributes to local adaptation of S. marinoi. Taken together, our data show that local adaptation of phytoplankton to complex coastal environments, which are characterised by a multitude of environmental gradients, is driven by widespread changes in diverse metabolic pathways and functions.