Distinct molecular responses to acute cold exposure revealed by comparative transcriptomic and metabolomic profiling in the bay scallop Argopecten irradians.

Abstract

Acute cold stress can elicit distinct molecular responses even when bay scallop populations show similar phenotypic outcomes. We compared a seventh-generation fast-growing bay scallop line (BS) with a commercial control population (CC) during a 72-h acute cold exposure at -1 ± 0.3 °C. RNA-seq was used as the discovery layer, representative BS cold-responsive genes were evaluated by qRT-PCR, and paired LC-MS profiles provided a comparative metabolic layer. At baseline, 138 genes differed between BS and CC; after cold exposure, 134 of these baseline differences disappeared and 61 of 65 cold-state differences newly emerged. BS showed a larger transcriptomic response magnitude than CC, with 1129 cold-responsive genes compared with 28 genes in CC, and this ordering remained robust across multiple sensitivity analyses. Survival after 72 h was identical in BS and CC (83/90, 92.2% in each population). Biochemical responses were time-dependent and marker-specific: CAT, LZM, T-SOD and T-AOC showed population-by-time interactions, whereas GSH-Px and MDA did not, and the 72-h differences were not consistently favourable to BS. Metabolomic cold effects were strongly concordant between populations, and no feature showed a significant population-by-cold interaction. Features putatively assigned to arachidonic acid metabolism were enriched, but this provider-annotated pathway signal remains exploratory because authentic-standard confirmation was not performed. These findings indicate population-specific differences in molecular responsiveness but do not establish superior cold tolerance in BS.