Long non-coding RNAs (lncRNAs) play crucial roles in a variety of biological processes, including stress response. However, the number, characteristics and stress-related expression of lncRNAs in turbot are still largely unknown. In this study, a total of 12,999 lncRNAs were identified at the genome-wide level of turbot for the first time using 24 RNA-seq datasets. Sequence characteristic analyses of transcripts showed that lncRNA transcripts were shorter in average length, lower in average GC content and in average expression level as compared to the coding genes. Expression pattern analyses of lncRNAs in 12 distinct tissues showed that lncRNAs, especially lincRNA, exhibited stronger tissue-specific expression than coding genes. Moreover, 612, 1351, 1060, 875, 420 and 1689 differentially expressed (DE) lncRNAs under Vibrio anguillarum, Enteromyxum scophthalmi, and Megalocytivirus infection and heat, oxygen, and salinity stress conditions were identified, respectively. Among them, 151 and 62 lncRNAs showed differential expression under various abiotic and biotic stresses, respectively, and 11 lncRNAs differentially expressed under both abiotic and biotic stresses were selected as comprehensive stress-responsive lncRNA candidates. Furthermore, expression pattern analysis and qPCR validation both verified the comprehensive stress-responsive functions of these 11 lncRNAs. In addition, 497 significantly co-expressed target genes (correlation coefficient (R) > 0.7 and q-value < 0.05) for these 11 comprehensive stress-responsive lncRNA candidates were identified. Finally, GO and KEGG enrichment analyses indicated that these target genes were enriched mainly in molecular function, such as cytokine activity and active transmembrane transporter activity, in biological processes, such as response to stimulus and immune response, and in pathways, such as protein families: signaling and cellular processes, transporters and metabolism. These findings not only provide valuable reference resources for further research on the molecular basis and function of lncRNAs in turbot but also help to accelerate the progress of molecularly selective breeding of stress-resistant turbot strains or varieties.

There are obvious individual differences in the growth and development of Cynoglossus semilaevis, mainly due to female bias. We selected 500 female Cynoglossus semilaevis of different sizes for GWAS and transcriptome analysis to screen for differential genes. qPCR was performed to detect the expression of the genes in various tissues, and RNAi experiments were performed in testicular cells to knock down the grid1 and grid2 genes and transcriptome sequencing was performed to check the changes of the downstream genes. Grid gene was screened for the common genes by GWAS sequencing and transcriptome sequencing. In the QPCR results, the expression of the grid gene family was negatively correlated with fish size, and was slightly higher in males than in females; in the transcriptome results, the expression of shcbp1, sass6, cdca7, and gh was up-regulated, and the expression of igf1 was down-regulated. It is speculated that igf1 has an antagonistic effect on gh, which is deregulated when the grid gene family is knocked down. The grid gene family may affect the growth of individual Cynoglossus semilaevis through the gh-igf1 axis, which provides a basis for the study of the differences in the growth size of Cynoglossus semilaevis.