Metamorphosis is a captivating process of change during which the morphology of the larva is completely reshaped to face the new challenges of adult life. In the case of fish, this process initiated in the brain has traditionally been considered to be a critical rearing point and despite the pioneering molecular work carried out in other flatfishes, the underlying molecular basis is still relatively poorly characterized. Turbot brain transcriptome of three developmental stages (pre-metamorphic, climax of metamorphosis and post-metamorphic) were analyzed to study the gene expression dynamics throughout the metamorphic process. A total of 1570 genes were differentially expressed in the three developmental stages and we found a specific pattern of gene expression at each stage. Unexpectedly, at the climax stage of metamorphosis, we found highly expressed genes related to the immune response, while the biological pathway enrichment analysis in pre-metamorphic and post-metamorphic were related to cell differentiation and oxygen carrier activity, respectively. In addition, our results confirm the importance of thyroid stimulating hormone, increasing its expression during metamorphosis. Based on our findings, we assume that immune system activation during the climax of metamorphosis stage could be related to processes of larval tissue inflammation, resorption and replacement, as occurs in other vertebrates.

Metamorphosis is a widely studied post-embryonic process in which many tissues undergo dramatic modifications to adapt to the new adult lifestyle. Flatfishes represent a good example of metamorphosis in teleost fishes. During metamorphosis of flatfish, organ regression and neoformation occur, with one of the most notable changes being the migration of one of the eyes to the other side of the body. In order to create a useful and reliable tool to advance the molecular study of metamorphosis in flatfish, we generated a chromatin accessible atlas as well as gene expression profile during four developmental stages ranging from a phylotypic to a post-metamorphic stage. We identified 29,019 differentially accessible chromatin regions and 3,253 differentially expressed genes. We found stage-specific regulatory regions and gene expression profiles, supporting the quality of the results. Our work provides strongly reproducible data for further studies to elucidate the regulatory elements that ensure successful metamorphosis in flatfish species.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).