RNA-seq analysis reveals divergent adaptive response to hyper- and hypo-salinity in cobia, Rachycentron canadum

Salinity is an important abiotic stress that affects metabolic and physiological activities, breed, development, and growth of marine fish. Studies have shown that cobia (Rachycentron canadum), a euryhaline marine teleost fish, possesses the ability of rapid and effective hyper/hypo iono- and osmoregulation. However, genomic studies on this species are lacking and it has not been studied at the transcriptome level to identify the genes responsible for salinity regulation, which affects the understanding of the fundamental mechanism underlying adaptation to fluctuations in salinity. To describe the molecular response of cobia to different salinity levels, we used RNA-seq analysis to identify genes and biological processes involved in response to salinity changes. In the present study, 395,080,114 clean reads were generated and then assembled into 65,318 unigenes with an N50 size of 2758 bp. There were 20,671 significantly differentially expressed genes (DEGs) including 8805 genes adapted to hypo-salinity and 11,866 genes adapted to hyper-salinity. These DEGs were highly represented in steroid biosynthesis, unsaturated fatty acid metabolism, glutathione metabolism, energy metabolism, osmoregulation, and immune response. The candidate genes identified in cobia provide valuable information for studying the molecular mechanism of salinity adaptation in marine fish. Furthermore, the transcriptomic sequencing data acts not only as an important resource for the identification of novel genes but also for further investigations regarding cobia biology.