Temperature optima of a natural diatom population increases as global warming proceeds

Studies in laboratory-based experimental evolution have demonstrated that phytoplankton species can rapidly adapt to higher temperatures. However, adaptation processes and their pace remain largely unknown under natural conditions. Here, by comparing resurrected Skeletonema marinoi strains from the Baltic Sea during the past 60 years, we show that modern S. marinoi have increased their temperature optima by 1 °C. With the increasing ability to grow in higher temperatures, growth rates in cold water decreased. Modern S. marinoi modified their valve:girdle ratio under warmer temperatures, which probably increases nutrient uptake ability. This was supported by the upregulation of several genes related to nitrate metabolism in modern strains grown under high temperatures. Our approach using resurrected strains demonstrates the adaptation potential of naturally occurring marine diatoms to increasing temperatures as global warming proceeds and exemplifies a realistic pace of evolution, which is an order of magnitude slower than estimated by experimental evolution. The Anthropocene has moved the planet into a new human-mediated geological epoch 1 , causing a rapid loss of oceanic biodiversity on a global scale 2 . The increasing average temperature is one of the most evident human-induced changes. It has wide-reaching effects on many organisms 2 due to the temperature dependency of biological processes 3 . This is especially relevant for marine organisms, as the ocean is a sink for most surplus heat 4 . Global sea surface temperature (SST) has risen by 0.7 °C and is projected to increase substantially by the end of the century 5–8 . Our study focuses on adaptation to human-induced global warming in a keystone phytoplankton species in the Baltic Sea—a region considered a ‘time machine’ for future ecosystem change due to experiencing warming levels above the global

Studies in laboratory-based experimental evolution have demonstrated that phytoplankton species can rapidly adapt to higher temperatures.However, adaptation processes and their pace remain largely unknown under natural conditions.Here, by comparing resurrected Skeletonema marinoi strains from the Baltic Sea during the past 60 years, we show that modern S. marinoi have increased their temperature optima by 1 °C.With the increasing ability to grow in higher temperatures, growth rates in cold water decreased.Modern S. marinoi modified their valve:girdle ratio under warmer temperatures, which probably increases nutrient uptake ability.This was supported by the upregulation of several genes related to nitrate metabolism in modern strains grown under high temperatures.Our approach using resurrected strains demonstrates the adaptation potential of naturally occurring marine diatoms to increasing temperatures as global warming proceeds and exemplifies a realistic pace of evolution, which is an order of magnitude slower than estimated by experimental evolution.
The Anthropocene has moved the planet into a new human-mediated geological epoch 1 , causing a rapid loss of oceanic biodiversity on a global scale 2 .The increasing average temperature is one of the most evident human-induced changes.It has wide-reaching effects on many organisms 2 due to the temperature dependency of biological processes 3 .This is especially relevant for marine organisms, as the ocean is a sink for most surplus heat 4 .Global sea surface temperature (SST) has risen by 0.7 °C and is projected to increase substantially by the end of the century [5][6][7][8] .Our study focuses on adaptation to human-induced global warming in a keystone phytoplankton species in the Baltic Sea-a region considered a 'time machine' for future ecosystem change due to experiencing warming levels above the global average 9 .
Unicellular phytoplankton, accounting for 40% of global primary production 10 , are essential contributors to oxygen generation, carbon sequestration and biogeochemical cycles 11 , and constitute the foundation of marine food webs 12 .Distinct thermal responses of marine phytoplankton species could lead to alterations in community composition and geographical distribution with increasing SST 13 .Further, productivity and diversity are expected to decline under increasing SST if species cannot shift their distribution range or adapt to the novel environment 14 .This may be especially true for partly isolated populations as studied here 15 .However, the widespread correlation between temperature optima (T opt ) in phytoplankton and SST over a 150° latitudinal gradient shows that phytoplankton can adapt to different temperatures 14 .
The high adaptive capacity of phytoplankton has been confirmed by experimental evolution 16 , indicating rapid adaptation to new environments through selection on existing genetic diversity 17 or de novo mutations 18 .Temperature adaptation to a delta of ≥4 °C has resulted in a shift either in the T opt (refs.18-23) or the upper thermal limit 18,21,22 . https://doi.org/10.1038/s41558-024-01981-9 decreasing water density at high temperatures 33 .In addition, higher temperatures can indirectly favour smaller cells through increased resource competition and reproductive rate.Smaller cells have a larger surface:volume ratio compared with larger cells and thus a larger relative area supporting resource acquisition 34 , making them better competitors for resources.A shift towards smaller cells was observed in coccolithophores 18 and green algae 23 adapted to high temperatures.However, there is conflicting evidence for diatoms, with some species increasing and others decreasing in cell size in response to increasing temperature 25 .Overall, cell size decline is important to consider, as it entails the potential for far-reaching ecological consequences, including the observed productivity decline in open oceans 35 .
Here we study potential temperature adaptation in a natural diatom population.We compared thermal response curves, cell size and morphology, and gene expression of resurrected early-Anthropocene (1960s) strains of the key diatom species Skeletonema marinoi to strains subjected to increasing global warming (1990s and 2010s).The spring-blooming marine diatom S. marinoi is a key primary producer in the Baltic Sea 36 and may periodically constitute up to 60-80% of the total biomass 37 .The study area in the northern Baltic Sea has been subjected to temperature increase and eutrophication during the past decades.Therefore, we expected to observe adaptation to increasing temperatures manifested as a higher T opt in modern S. marinoi (2010s).Further, we expected to observe a shift in cell size in S. marinoi from the past 60 years when grown in higher temperatures, and differential gene expression of metabolic pathways that are linked to the observed shifts in thermal reaction norms and related cell-size shifts.

The temperature optima of S. marinoi increased over the past 60 years
Strains of S. marinoi isolated across the past 60 years showed a shift in temperature-dependent growth and thermal performance curves (Fig. 1).The T opt of modern strains shifted by almost 1° compared with strains from the early Anthropocene (Fig. 1a,b; F 2,17 = 4.98, P = 0.019).The mean T opt of 14.99 °C in the 1960s strains shifted to 15.5 °C in the 1990s strains and increased significantly to 15.88 °C in the 2010s strains (Fig. 1b).A shift of ~1.50 °C was observed in the lower temperature limit (T min ; Fig. 1a,c; F 2,17 = 15.02,P < 0.001).The strains from the 1960s showed a higher growth rate at low temperatures (Fig. 1a).Their T min (2.93 °C) was significantly lower compared with the 1990s and 2010s strains (3.79 °C and 4.42 °C, respectively; Fig. 1c).No differences in the upper temperature limit (T max ) were observed between decades (Supplementary Fig. 5).The maximum growth rate (µ max ) was highest in the 1990s strains at 1.01 day −1 , which was 15% higher compared with the 1960s strains.
Modern S. marinoi display a shift in cell shape in response to increased temperature.Shifts in cell size (measured as volume) along the temperature gradient from 6-22 °C were significantly different between strains of S. marinoi from different decades (Fig. 2a; F 2,158 = 75.26,P < 0.001).While the size was rather similar across temperatures in strains from the 1960s and the 1990s, the modern population increased in size by 220% from low to high temperatures.Shifts in cell size affect the surface to volume ratio of the cells, which similarly showed an interaction between the temporal populations and temperature (Fig. 2b; F 2,158 = 94.86,P < 0.001).The response within strains from the 1960s and 1990s was similar across temperatures, while a strong decrease (40%) occurred with increasing temperatures within the strains from the 2010s.These shifts were driven by changes in cell width (Fig. 2c,d).The cell length decreased by ~30% across strains from all decades.In contrast, shifts in cell width across temperatures were dependent on the age of strains (Fig. 2b; F 2,158 = 77.00,P < 0.001).An increase of 120% was observed within the strains from the 2010s, while strains from the 1990s and 1960s showed no shift.Also, increased growth and photosynthetic rates are linked to temperature adaptation 23,24 .Further, analyses of differential gene expression have revealed changes in pathways related to photosynthesis and energy metabolism 25,26 .However, the constraints on temperature adaptation, whether due to thermodynamics 6 , generalist-specialist trade-offs 19 or resource allocation 22 , remain a subject of ongoing debate.For example, temperature adaptation seems limited under low nitrogen concentrations 20,27 .Growth at elevated temperatures leads to an increased demand for nitrogen to sustain energetically costly repair mechanisms.This demand cannot be met under nutrient limitation because of trade-offs between allocating resources to reproduction and nutrient uptake 22 .Overall, experimental evolution studies have enhanced our understanding of phytoplankton evolution, but it remains uncertain whether the adaptive potential observed in laboratory experiments translates to real-world conditions.
Most experimental evolution studies do not include realistic selection pressures, natural diversity, interactions among individuals and species 28 , or sexual reproduction, all of which can alter the adaptation potential and/or the rate of evolution.Therefore, other approaches are required for studying adaptation under real-world global warming.This is possible through a 'backward-in-time' method, by resurrecting phytoplankton resting stages from seafloor sediment archives [29][30][31] .A previous resurrection study using phytoplankton demonstrated a shift in life cycle processes with increasing temperatures 29 .Analysing resurrected strains provides insights into how evolutionary processes, under a natural pace of global warming, compare to evolution occurring under simulated laboratory conditions.
With increasing SST, the general expectation is that phytoplankton cells become smaller 32,33 .One postulated mechanism is that size shifts allow cells to maintain the same sinking velocity with https://doi.org/10.1038/s41558-024-01981-9 Gene expression analysis revealed evolutionary effects on nutrient metabolism.On average, 98.2% of the RNA-seq reads per sample mapped to the S. marinoi reference genome (total size ~55 Mb).We observed no differences in mapping success between different strains (Supplementary Table 1).We compared potential differences in gene expression and found that in total, 8,280 of the 22,438 predicted genes were differentially expressed (DE) using a 5% false discovery rate (FDR) level.Approximately 76% or 6,328 of these DE genes (DEGs) received a functional annotation.When visualizing the top 500 DEGs, we observed a clear difference in expression patterns between strains from different time points (Fig. 3).The number of DEGs was higher between decades, when grown under the same temperature.The most pronounced difference was observed between the 1960s and the 2010s (4,948 ± 55 DEGs).This was significantly more than what was observed in the 1960s versus 1990s (4,021 ± 158) and the 1990s versus 2010s contrasts (3,532 ± 106) (t-test, t = 3.18, P < 0.05) (Fig. 4 and Supplementary Fig. 1).
We also observed a difference when comparing strains within time points subjected to different temperatures (Fig. 3).Between temperatures, strains from the 1960s displayed a total of 1,062 ± 138 DEGs (per strain), which is significantly less compared with that of the 1990s (1,302 ± 189) (t-test, t = 2.44, P < 0.011).Also, strains from the 2010s had significantly fewer DEGs (913 ± 240) between temperatures compared with strains from the 1990s (t-test, t = 2.44, P < 0.001) (Fig. 4).There was no difference between the 1960s and 2010s in the number of DEGs between temperatures (t-test, t = 2.44, P = 0.643).Especially in the 14 °C versus 20 °C and the 8 °C versus 14 °C contrasts, a low number of DEGs was observed, while a higher number of DEGs was observed in the 8 °C versus 20 °C contrast (Supplementary Fig. 2).
We consider genes that were uniquely DE in 2010s strains when comparing 8 °C to 20 °C to represent a portion of evolved functionality in S. marinoi (Fig. 5).There was a total of 416 such genes when including genes with log fold change (FC) values >2 or <−2.These included, for example, upregulation of thioredoxin (Acht1, Acht4), trypsin (Loc5578510) and one heat-shock gene (Hsf1) (Fig. 5).In contrast, we observed several different upregulated heat-shock genes and transcription factors (Hsf1, Hsf4, Loc4342550) in strains from the 1960s (Supplementary Table 2).When grouping genes with related functions into Gene Ontology (GO) terms, we observed a downregulation of functions within, for example, nicotinamide adenine dinucleotide phosphate (NADP) biosynthetic energy metabolism processes in strains from the 1990s and 2010s with increasing temperature (Fig. 6).This was not observed in strains from the 1960s.Also, differential gene expression in photosynthesis-related processes was observed in strains from the 1990s and 2010s but not seen in the 1960s strains.Mitotic sister chromatid cohesion (GO:0007064) was significantly upregulated in the 1990s and 2010s strains, which was not seen in the strains from the 1960s.In addition to shifts in these biological processes, we observed changes in related energy processes at the molecular function level (Supplementary Fig. 3).Changes to the photosynthesis machinery in modern strains were further supported by differential gene expression affecting cellular components within photosystem I and II (Supplementary Fig. 4).

Discussion
Despite mounting evidence that thermal adaptation is possible under controlled conditions [19][20][21]23 , it remains uncertain how adaptation plays out under natural conditions. Mutiple drivers including species interactions and abiotic environmental changes 27,38 can alter the selection pressure.Adaptation may also be affected by the rate of temperature change over decades and by diurnal and annual fluctuations in temperature conditions 39 .This plethora of contributing factors is essential for understanding 'real-world' selection.Here we demonstrate that a natural phytoplankton population can increase its T opt apace with natural global warming over six decades.This provides evidence for a realistic pace of evolution of phytoplankton to global warming and several other metrics that have been under selection from global warming.We demonstrate that modern S. marinoi decreases its surface to volume ratio with increasing temperatures. This as not observed in the strains from the 1960s and 1990s.We observed shifts in how the population has altered its gene expression in relation to the increase in SST in the study area.We suggest that strains from the 1960s experienced higher temperatures as more stressful than the recent strains. Ths is supported by their lower growth rates in high temperatures and by a higher number of upregulated genes coding for heat-shock proteins.Heat-shock proteins are known to repair cell damage under high-temperature conditions 40 .The less-clear stress response in strains from 2010s exposed to above-optimum In experimental evolution, temperature increase is frequently applied at the upper temperature limit of the species, posing a strong selective pressure.Here we have investigated the effect of a gradual increase in SST of ~1.5 °C in the study area since the 1960s.This entails a more subtle selection pressure compared with most experimental evolution studies that use a delta of ≥4 °C (refs.19-21,23).The upper temperature tolerance limit of S. marinoi (>27 °C; Supplementary Fig. 5 and ref. 41) is never exceeded during the main growth season in this study area (Supplementary Fig. 6).This diverging selection pressure between experimental and natural evolution might partly explain why a comparable increase in T opt required more generations in our natural population compared with populations exposed to experimental evolution.For example, Chlamydomonas reinhardtii exposed for a decade to +4 °C above ambient temperature in a mesocosm study showed a 1.6 °C increase in the T opt (ref.23).A comparable shift of 1-2 °C in T opt was observed in two diatom species grown for 200-600 generations at 4 °C above ambient temperature 19 .Using the climatology of the sampling area and the temperature response curves, we estimate that the observed 1 °C increase in T opt of S. marinoi required ~7,000 mitotic generations (Supplementary Fig. 6 and Supplementary Table 3).Skeletonema marinoi undergoes sexual reproduction and meiosis can be induced experimentally.However, it is unknown how often sexual events occur in nature 42 .Therefore, we are adhering to asexual generations in line with most laboratory evolution experiments.Also, other changing drivers in the environment (for example, light conditions and grazing), for which we have not accounted here, may have affected the population dynamics.Overall, our study demonstrates that the rate of temperature adaptation, while slower than initially estimated through experimental evolution, enables a population to adapt to ongoing global warming.
In addition, it is important to consider diurnal and seasonal temperature fluctuations, which may have contradictory effects on temperature adaptation.Diurnal fluctuations have been described to accelerate the molecular evolution of thermal tolerance in the diatom Thalassiosira pseudonana compared with constant exposure 39 .Further, diurnal fluctuations can lead to the evolution of plasticity 43 .Using our data in the framework of a reaction norm (as in ref. 44) suggests that the evolution of plasticity does not play a notable role in the temperature adaptation of modern strains (Supplementary Fig. 7).Seasonal temperature fluctuations may have a contradictory role by slowing down evolution.An increase in T opt is often accompanied by a performance trade-off, including a reduced growth rate at low temperatures 19,22 .We observed a higher growth capacity at low temperatures in early-Anthropocene strains and a decrease in this capacity as the optimum increased.Consequently, strains with high temperature optimum favoured in late spring can have a disadvantage during cold winter seasons.The maintenance of high phenotypic diversity in isolated strains across 60 years of selection under global warming suggests temporal fluctuation in selection pressure, which favours the maintenance of high diversity 45 .Moreover, strong coupling between the pelagic population and benthic resting stages could mitigate adaptation to short-term environmental fluctuations.Benthic-pelagic coupling has been suggested to result in a homogeneous population structure across the seasons in dinoflagellates in the Baltic Sea 46 .Strong benthic-pelagic coupling has also been observed in S. marinoi in the study area 47,48 .Thus, it is highly unlikely that the temperature adaptation we observed was driven by seasonal or short-term fluctuations.However, increasing winter temperatures could shift the balance in favour of strains with higher T opt .The selection for higher T opt and reduced growth capacity at low temperatures might be one of the reasons underlying the later onset of the S. marinoi spring bloom observed elsewhere 49,50 .
The strains isolated from the 1990s displayed the highest maximum growth rates.This may be explained by a trade-off between https://doi.org/10.1038/s41558-024-01981-9 growth and nutrient affinity.Generally, this trade-off arises as cellular resources can be allocated to reproduction or nutrient uptake 51 .
The nitrogen requirement of repair mechanisms such as heat-shock proteins is expected to increase during temperature adaptation.Nitrogen limitation can thus constrain adaptation to increasing temperatures 20 , making phytoplankton in areas with low nutrient conditions more vulnerable to global warming 52 .Hence, even though not directly tested here, adaptation to global warming in the population of S. marinoi in this study might have been enabled by the strong and documented eutrophication in the area 53 .When the limiting nutrient is present in excess, it gives individuals with higher growth a competitive advantage 54 .A higher growth rate is positively correlated with a higher maximum uptake velocity (V max ), which is generally observed in the velocity-adapted diatoms.These traits are advantageous under high nutrient environments 55 .In our study area, eutrophication peaked in the 1990s 9,56 , which may explain the high growth rates in strains from this time.In agreement, a recent study from the Baltic Sea comparing populations of S. marinoi growing https://doi.org/10.1038/s41558-024-01981-9 under different trophic states showed that eutrophication probably drives selection for faster growth 57 .This is further supported by the differential gene expression relating to energy metabolism and cell division (for example, GO:0046496, GO:0006739, GO:0007064).These processes are related to NADP, the central electron carrier during the light-dependent photosynthesis reactions which were most differentially expressed in strains from the 1990s and 2010s.Overall, this prompts the question of whether adaptation with ongoing global warming will be of broad applicability or if it is limited to regions subject to strong eutrophication.Elevated temperatures tend to decrease individual cell sizes 32,33 due to resource constraints and enhanced reproductive rates.However, the silica frustules of diatoms may contribute to conflicting evidence of size changes in higher temperatures 25 .Here we show opposing responses in size and surface:volume ratio within the same species with different T opt .While both early-Anthropocene and modern strains showed a reduced cell length under increasing temperatures, modern strains with higher T opt showed an increase in cell width.This is in line with shifts in cell size associated with temperature evolution in the diatom Thalassiosira pseudonana 58 .Diatoms respond to resource limitation by adjusting the valve:girdle ratio.This ratio is relevant for nutrient acquisition, as the silica frustules restrict nutrient uptake at the elongated girdle band (length) while facilitating it at the circular valves (width) that are equipped with punctae.With increasing valve:girdle ratio, the ratio of the surface area of the valve to the volume of the cell also increased 58 (Supplementary Fig. 8) and the cell shape changed (Supplementary Fig. 9).Overall, the counterintuitive increase in size with increasing temperature adaptation is probably driven by an increased need for resources.We also found support for this higher demand for nutrients in the gene expression of modern S. marinoi.
Among the upregulated genes in high temperature, we observed several with functions related to trypsin metabolism, which is a known regulator of N:P stoichiometric homoeostasis in phytoplankton 59 .The expression of trypsin in diatoms is especially responsive to shifts in the environment 60 .Here, it may be the key to fulfil the higher demand for nitrate when temperatures have increased and the competition for this resource simultaneously intensifies.Further, two thioredoxins, which are central regulators of CO 2 fixation and nitrogen in chloroplasts 61 , were highly upregulated in strains from the 2010s.This suggests some modification of the nutrient acquisition and photosynthesis machinery in modern S. marinoi strains.

Conclusions
We conclude that S. marinoi in the Baltic Sea has adapted to ongoing global warming in the past 60 years.Our gene expression data and the observed shifts in cell morphology support earlier experimental evolution studies showing that nutrient conditions have the potential to affect adaptation to temperature increases in phytoplankton.In agreement, we also observed a trade-off between high growth at warm versus cold temperatures.However, the number of generations required by the natural population to reach the same evolutionary change was an order of magnitude greater.Overall, the underlying mechanisms of evolutionary change can be well understood in experimental studies, but the estimation of the rate of evolution requires the study of natural communities.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions Fig. 6 | Highly significant GO terms in high temperatures.Only biological processes that were significantly up-or downregulated in a higher temperature (8 °C→14 °C, 8 °C→20 °C, 14 °C→20 °C) are represented by a bar showing the proportion of DEGs per biological process category (gene ratio).The grey scale represents the P value from one-tailed Fisher's exact test.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.© The Author(s) 2024 https://doi.org/10.1038/s41558-024-01981-9

Model organism
The spring bloom in the northeast Baltic Sea is often dominated by the marine diatom S. marinoi which serves as an important food source for zooplankton 37 .This model species has been extensively studied in terms of biogeography, physiology and genetic variation across space and time in temperate areas of the world 62,63 .It is known to form resting stages that sink to the seafloor when the bloom phase is over towards the end of the spring bloom 64 .This 'biological archive' may contain up to ~57,000 S. marinoi cells per gram sediment 64 , with the potential to reveal the diversity and ecophysiology of past populations.

Study area, sediment sampling and age modelling
The Småholmen station (60.24°N, 22.04° E) close to Haverö in the Archipelago Sea is 20 m deep and has been seasonally anoxic or hypoxic at the seafloor for at least seven decades 65 .The low oxygen conditions and lack of bioturbation have resulted in the formation of laminated sediment at this site.Such conditions are optimal for preserving the chronology of sedimented cells.Historical data show that the average spring temperature (air) in the Archipelago Sea was stable at ~2 °C from the 1880s until the 1970s.Since the mid-1970s, it has increased by ~0.5 °C per decade 66 .An increase in the SST has been observed close to our study site (Seili monitoring station located <3 km from Småholmen).In April, when the spring bloom reaches its peak, SST has increased by ~2.5 °C since the early 1980s (Supplementary Fig. 10).
Sediment cores were retrieved in April 2020 using an HTH Kajak surface sediment gravity corer 67 at the Småholmen station.After core retrieval, the core tube containing undisturbed sediment profiles was attached to a stand on-site and carefully sliced into 2 cm subsamples.The outer edge (3 cm) of the entire sediment core was discarded to avoid cross-contamination between different depth layers due to smearing along the outer edge of the core.The subsamples were stored in the dark at 8 °C until further processing.An age-depth model for the sediment core was constructed with the Undatable software 68 using an xfactor of 0.1 and 15% bootstrapping (Supplementary Fig. 11).The xfactor value determines the sediment accumulation rate uncertainty between consecutive age-depth constraints and the bootstrapping function randomly removes a selected percentage of the age-depth constraints in each run of the age-depth simulation 68 .Age constraints for the modelling procedure were obtained through loss on ignition (LOI) correlation against previous studies in the same location 47 .LOI was determined at 2 cm intervals by drying the subsamples at 105 °C for 16 h and ashing at 550 °C for 2 h.The LOI correlation was further verified by 137 Cs dating of a replicate sediment core retrieved for this study (Supplementary Fig. 12).For this, 137 Cs activity of the untreated 2-cm-thick sediment slices was determined by gamma spectrometry using a BrightSpec bMCA-USB pulse height analyser coupled to a well-type NaI(Tl) detector 69 .In the northern hemisphere, 137 Cs contamination in sediments is mainly derived from the Chernobyl nuclear power plant accident in 1986 and atmospheric nuclear weapons testing in the early 1960s 70 .Therefore, the 137 Cs dating of Baltic Sea sediments is based on the recognition of these two horizons in the sediment profiles.

Resurrection of clonal cultures
Resurrection was initiated by mixing 0.5 g of sediment from different layers with each 50 ml of filtered (0.22 μm) seawater with f/2+Si medium 71 in June 2020.The sediment was from the 2-4 cm, 28-30 cm and 46-48 cm layers corresponding to median age ± 1σ error estimates of 2017-2019-2020 (referred to as ~2010s/modern), 1985-1992-2001 (~1990s) and 1957-1963-1969 (~1960s/early Anthropocene), respectively (Supplementary Fig. 8).The sediment 'slurries' were distributed on 24-well NUNC plates and incubated at 8 °C and 40 μmol photons m −2 s −1 (12 h:12 h light:dark cycle).When vegetative growth emerged (after 2-4 weeks), single chains of S. marinoi were isolated using a micropipette under an inverted light microscope (Nikon Diaphot 300).One chain per well was isolated to a new well to minimize the probability of isolating the same clone twice.As previous population genetic studies on S. marinoi have shown, the genotypic diversity of S. marinoi is extremely high 72 .Thus, it is highly unlikely to ever find the same genotype (using, for example, microsatellite markers) of this species in natural conditions when sampled over time or across space.Therefore, we considered that the strains included in this study were all different genotypes and consisted of a small percentage of the population's total intraspecific diversity.The isolated strains can theoretically stem from different seasons, as each layer contains resting stages across seasons from 2-5 years.We confirmed that the strains belonged to the species S. marinoi by sequencing the V4 region of the 18S ribosomal RNA gene (for more details see Supplementary Information, 'Genetic identification of S. marinoi').After 1-2 weeks, the entire volume in each well was transferred to 40 ml culture flasks containing 10 ml of f/2+Si media.After further growth for 7 days, the total volume was increased to 40 ml, and clonal cultures were maintained in the same conditions as described above.The resurrection rate varied between 44 and 58% (isolated cells that made it to a stable culture).The rates for the 1960s, 1990s and 2010s were 52%, 58% and 44%, respectively.

Experiments
To test whether S. marinoi has adapted to an increase in temperature during the past ~60 years, we assayed seven strains per time point (1960s, 1990s, 2010s) in +6 °C, +8 °C, +10 °C, +14 °C, +18 °C, +20 °C, +22 °C and +26 °C in natural seawater (6 PSU, 0.22 μm filtered, with f/4+Si) and 100 μmol photons m −2 s −1 (12 h:12 h light:dark cycle).Due to technical limitations, each temperature was assessed in a separate experiment in a randomized order (Supplementary Table 4).Randomization minimizes the influence of measurement timing on different temperature treatments by ensuring that the timing effect is evenly distributed across the temperature range.Before each experimental start, strains were acclimated to f/4+Si media for 1 week.The experiment was started by transferring cells to new media, reaching a start concentration of 10,000 cells per ml.To avoid dilution of nutrients, the maximum inoculum was 5 ml (10% of the total volume).The growth of each strain (four replicates) was monitored daily by measuring the in vivo fluorescence of chlorophyll a in a 300 μl subsample on a 96-well plate (PerkinElmer, IsoPlate 96F) using a spectrophotometer (Tecan, Infinite 200 Pro with the software Magellan for Tecan Infinite Pro v.1) until the stationary phase was reached.The excitation wavelength was set to 425 nm and the emission wavelength to 680 nm.To account for uneven distribution of cells, nine positions were measured in each well.The cells were not dark-adapted, but to mitigate potential confounding effects of light adaptation during the diurnal light-dark cycle, measurements were consistently taken at the same time.At the end of the experiment, 1 ml subsamples were fixed with acidic Lugol's solution for later estimation of cell size.Three randomly chosen strains from each time point were consistently sampled for RNA when the stationary phase was reached (Supplementary Table 5).About 25-35 ml of cell culture was centrifuged for 30 min at 4 °C at 3,900 × g (Eppendorf, 5810R).The cell pellet was resuspended in 400 μl of TRIzol (Invitrogen) and incubated for 2 min at 60 °C until completely dissolved.The cells were immediately stored at −80 °C until RNA extraction.

Thermal performance curves based on growth
Growth was calculated by fitting a linear model ('easy_linear') to blank-corrected fluorescence values using the 'growthrate' package 73 in R (v.4.3.0) 74and R Studio (v.2023.9.0.463) 75 .We excluded all replicates that showed negative growth (all replicates in 26 °C, 1960_05/06 in 22 °C) or unstable growth (maximum fluorescence <2,000) before fitting the growth model.For the samples with negative growth, the growth rate was set to zero for further analysis.Strain 2010_13 was potentially contaminated and was excluded from all downstream analyses.The calculated growth rates (total n of 216, 222 and 184 for

Fig. 1 |
Fig. 1 | Temperature response of modern and early-Anthropocene S. marinoi.a, Mean thermal performance curves of strains from the 1960s, 1990s and 2010s (green, orange and purple, respectively).Underlying thin lines show individual performance curves of 7 strains per population.Thermal performance curves were fitted as a quadratic model.Maximum growth rate denoted by µ. b-d,Mean and confidence interval (CI) of optimum temperature (T opt ) (b), lower temperature limit (T min ) (c) and maximum growth rate (µ max ) (d) were estimated for each strain by bootstrapping (n = 4 replicates per strain).Overlying strainspecific responses, the mean and CI of all strains per time point and differences between time points (lowercase grey letters) are shown.

Fig. 2 |Fig. 3 |
Fig. 2 | Size of modern and early-Anthropocene S. marinoi.a-d, Volume (a), surface to volume ratio (b; S:V), width (c) and length (d) of 3 strains from the 1960s (green), 1990s (orange) and 2 strains from the 2010s (purple).Points are the mean of 50 individual cells measured per replicate.A small jitter of data points was added around the temperature values on the x axis to enhance the visibility of individual measurements.Model predictions are shown as lines.

Fig. 4 |
Fig. 4 | The number of DEGs (up-and downregulated) within and between decades.The strains originate from the 1960s, 1990s and 2010s across all temperature contrasts.The number of DEGs is given as mean per strain.

Fig. 5 |
Fig. 5 | Unique DEGs in strains from the 2010s.The heatmap includes genes in strains from the 2010s that were uniquely up-or downregulated when compared at 8 °C vs 20 °C.Only the top 78 genes with logFC values >3 or <−3 are displayed.