Full-Length Transcriptome Sequencing and hsp Gene Family Analysis Provide New Insights into the Stress Response Mechanisms of Mystus guttatus.

Abstract

Mystus guttatus, a second-class protected species in China, has undergone severe population decline due to anthropogenic and environmental pressures, yet conservation efforts are hindered by limited genomic resources and a lack of mechanistic insights into its stress response systems. Here, the first full-length transcriptome of M. guttatus was generated via SMRT sequencing. A total of 32,647 full-length transcripts were obtained, with an average length of 1783 bp. After structure and function annotation of full-length transcripts, 30,977 genes, 1670 transcription factors (TF), 918 alternative splicing (AS), and 11,830 simple sequence repeats (SSR) were identified. In order to further explore the stress resistance of M. guttatus, 93 genes belonging to the heat shock protein (HSP) family were identified and categorized into HSP70 and HSP90 subgroups. After phylogenetic analysis and selective stress analysis, it was discovered that the hsp family has suffered purifying selection and gene loss, potentially contributing to a decrease in the stress resilience and population of M. guttatus. Using protein interaction network and molecular docking tools, we observed the intricate interplay among HSPs and discovered HSP70-HOP-HSP90 interaction, which is an essential stress response mechanism. Our study sequenced the first full-length transcriptome of M. guttatus to enhance its genomic resources for its conservation and breeding and provide new insights into the future study of stress response mechanisms on M. guttatus.