26Cultural eutrophication is one of the largest threats to aquatic ecosystems with 27 devastating ecological and economic consequences. Although Daphnia, a keystone 28 grazer crustacean, can mitigate some negative effects of eutrophication in freshwater 29 habitats, it is itself affected by changes in nutrient composition. Previous studies have 30 shown evolutionary adaptation of Daphnia to environmental change, including high 31 phosphorus levels. While phosphorus is an essential macronutrient, it negatively 32affects Daphnia life history when present in excessive concentrations. Here, we map 33 weighted gene co-expression networks to phosphorus-related phenotypic traits in 34 centuries-old ancestral Daphnia and their modern descendants, contrasting pre-and 35 post-eutrophication environments. We find that evolutionary fine-tuning of 36 transcriptional responses is manifested at a basic (cellular) phenotypic level, and is 37 strongly correlated to an evolved plasticity of gene expression. At a higher phenotypic 38 level, this contributes to the maintenance of physiological homeostasis, and a robust 39 preservation of somatic phosphorus concentration. 40 41 3
Introduction 42Cultural eutrophication is one of the largest threats to aquatic ecosystems with 43 devastating ecological and economic effects 1 . Rising phosphorus (P) concentrations in 44 surface water result in high algal biomasses that pose a serious ecological threat to 45 aquatic habitats, impacting biodiversity, drinking water quality, recreational 46 resources, and fisheries 1-2 . One of the major P-sources of eutrophication is the use of 47 phosphate fertilizers, for which the global consumption has increased from 1 to 15 Mt 48 P year -1 in the 20th century 3 . Vital ecosystem services provided by Daphnia, a 49 keystone herbivore crustacean in freshwater habitats, can effectively mitigate some of 50 these negative impacts by decimating phytoplankton biomass 4 . Although P is a vital 51 element for many cellular components, including RNA, DNA and ATP, excessive P-52 concentrations typical of eutrophication can negatively affect somatic and population 53 growth of Daphnia 5-6 . In addition, algal blooms often involve toxin-producing 54 cyanobacteria that can negatively affect the growth rate of Daphnia, although local 55 adaptation to these toxins has been observed 7-9 . Insights from resurrection ecology 10 56 suggest that Daphnia can adapt to changing environmental conditions, such as climate 57 change 11 , anthropogenic stressors 12 , toxic algae 13 , and phosphorus supply 14 . Although 58 previous analyses have examined the functional genome underlying Daphnia 59 phenotypes (e.g. [15][16][17][18][19] ), the link between phenotypic adaptation and the evolution of 60 gene expression reaction norms in this keystone species is poorly understood. 61Adaptation to environmental change has been associated with divergent gene 62 expression patterns 20 and the maintenance or evolution of gene expression reaction 63 norms [21][22] . However, directly linking...