The HAB-forming, toxic dinoflagellate Karenia mikimotoi, previously found to benefit from ocean acidification (OA), was cultivated to investigate its transcriptional response to simulated OA for 30 generations. Batch cultures were grown under two CO2 concentrations, 450 (control) and 1100 (simulated OA) μatm, and physiological parameters [growth, pigments, catalase (CAT), glutathione reductase (GR), and superoxide dismutase (SOD) activity], as well as transcriptomes (obtained via RNA-seq), were compared. Chlorophyll a (Chl a) and carotenoid (Caro) contents, as well as CAT and GR activities, were significantly increased under OA conditions. Transcriptomic analysis revealed 2,490 differentially expressed unigenes in response to OA, which comprised 1.54% of all unigenes. A total of 1,121 unigenes were upregulated, and 1,369 unigenes were downregulated in OA compared to control conditions. The downregulated expression of bicarbonate transporter and carbonic anhydrase genes was a landmark of OA acclimation. Key genes involved in energy metabolism, e.g., photosynthesis, tricarboxylic acid cycle, oxidative phosphorylation, and nitrogen metabolism, were highly upregulated under OA, contributing to increases in the Chl a (55.05%) and Caro (28.37%). The enhanced antioxidant enzyme activities (i.e. CAT, GR) and upregulated genes (i.e. glutathione peroxidase, ascorbate peroxidase, heat shock protein, 20S proteasome, aldehyde dehydrogenase, and apolipoprotein) benefit cells against the potential lower pH stress condition under OA. In addition, the downregulation of four genes associated with motility suggested that the preserved energy could further boost growth. In conclusion, the present study suggests that K. mikimotoi exhibits efficient gene expression regulation for the utilization of energy and resistance to OA-induced stress. Taken together, K. mikimotoi appeared as a tolerant species in response to OA. Thus, more extensive algal blooms that threaten marine organisms are likely in the future. These findings expand current knowledge on the gene expression of HAB-forming species in response to future OA.

In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2021) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2022-06-08.