Genome-Wide Association and Expression Analysis Revealed the Candidate Variants and Molecular Underpinnings of Cold-Stress Response in Large Yellow Croaker

Large yellow croaker (Larimichthys crocea) is one of the most economically important fish in China. Recently, global climate change has caused more and more intense and extreme low temperature weathers, resulting in huge losses to the large yellow croaker industry. Therefore, it is essential to understand the mechanisms of low-temperature tolerance in large yellow croaker. Here, we conducted an integrative analysis of genome-wide association study (GWAS) and transcriptome analysis to identify candidate variants and reveal the molecular underpinning of cold-stress response in large yellow croaker. A total of 8 significant single nucleotide polymorphisms (SNPs) loci on 6 chromosomes were identified in the GWAS analysis, and 5764 (gill) and 3588 (liver) differentially expressed genes (DEGs) were detected in cold-stressed large yellow croaker, respectively. Further comparative and functional analysis of the candidate genes and DEGs highlighted the importance of pathways/genes related to immune response, cellular stress response, lipid transport, and metabolism in the cold-stress response of large yellow croaker. Our results provide insights into the cold tolerance of large yellow croaker and contribute to genomic-based selection for low-temperature-resistant large yellow croaker.

Low temperature-induced variation in plasma biochemical indices and aquaglyceroporin gene expression in the large yellow croaker Larimichthys crocea

Low temperature influences multiple physiological processes in fish. To explore the adaptability of the large yellow croaker (Larimichthys crocea) to low temperature, the concentrations of glycerol, blood urea nitrogen (BUN), and triglycerides (TG) in plasma, as well as the expression levels of metabolism-related genes aqp7 and aqp10, were measured after exposure to low temperature stress and during subsequent rewarming. In addition, tissue samples from the intestine and liver were histologically analyzed. We found that the concentrations of plasma glycerol, BUN, and TG, decreased under low temperature stress, suggesting the metabolism of fat and protein slowed at low temperature. The expression levels of aqp7 and aqp10 mRNA were also downregulated under exposure to low temperature. Interestingly, above plasma indices and gene expression returned to basic levels within 24 h after rewarming. Furthermore, the liver and the intestine were damaged under continuous low temperature stress, whereas they were repaired upon rewarming. From the above results, we concluded that aqp7 and aqp10 genes were sensitive to low temperature, and that the decrease in their expression levels at low temperature might reduce energy consumption by L. crocea. However, the adaptation to low temperature was limited because the key metabolic tissues were damaged under continuous exposure to low temperature. Interestingly, the metabolism of L. crocea was basically back to normal within 24 h of rewarming, showing that it has high capacity of self-recovery.

Molecular cloning and expression analysis of scd1 gene from large yellow croaker Larimichthys crocea under cold stress

Desaturation of fatty acids is an important adaptation mechanism to maintain membrane fluidity under cold stress. To comprehend the mechanism of adaptation to low temperatures in fish, we investigated stearoyl-CoA desaturase 1 (SCD1) endocrine expression in the process of cold acclimation from 15°C to 7°C in Larimichthys crocea. The cDNA and genomic sequences of scd1 were cloned and characterized and named as Lcscd1. The cDNA encoded an iron-containing protein of 337 amino acids with functional motifs. The full-length genome sequence of Lcscd1 was composed of 2556 nucleotides, including five exons and four introns. Tissue expression profiles by qPCR and western blot analysis revealed that Lcscd1 was highly expressed in the liver, followed by the brain. The expression of Lcscd1 mRNA in the liver was firstly down-regulated from 15°C to 11°C, and then up-regulated until the first day of 7°C, followed by a decline until the last day. In the brain, the expression showed no significant change from 15°C to 9°C, but then significantly increased until the last day of 7°C. SCD1 protein expression in the liver decreased from 15°C to the first day of 7°C, and then gradually recovered to the starting level. In the brain, SCD1 protein expression maintained rising trends in the whole process. Immunoelectron microscopic analysis showed that SCD1 was localized in fat granules, mitochondria and granular endoplasmic reticulum of hepatic cells, but only in mitochondria of encephalic cells. The results above suggested that SCD1 expression was responsive to both cold and starvation stresses in the liver, but only to cold stress in the brain. In conclusion, these findings suggested that SCD1 may be involved in fish adaptation to cold stress.