The circadian transcriptome of marine fish (Sparus aurata) larvae reveals highly synchronized biological processes at the whole organism level

The regulation of circadian gene expression remains largely unknown in farmed fish larvae. In this study, a high-density oligonucleotide microarray was used to examine the daily expression of 13,939 unique genes in whole gilthead sea bream (Sparus aurata) larvae with fast growth potentiality. Up to 2,229 genes were differentially expressed, and the first two components of Principal Component Analysis explained more than 81% of the total variance. Clustering analysis of differentially expressed genes identified 4 major clusters that were triggered sequentially, with a maximum expression at 0 h, 3 h, 9–15 h and 18-21 h zeitgeber time. Various core clock genes (per1, per2, per3, bmal1, cry1, cry2, clock) were identified in clusters 1–3, and their expression was significantly correlated with several genes in each cluster. Functional analysis revealed a daily consecutive activation of canonical pathways related to phototransduction, intermediary metabolism, development, chromatin remodeling, and cell cycle regulation. This daily transcriptome of whole larvae resembles a cell cycle (G1/S, G2/M, and M/G1 transitions) in synchronization with multicellular processes, such as neuromuscular development. This study supports that the actively feeding fish larval transcriptome is temporally organized in a 24-h cycle, likely for maximizing growth and development.