Genome-wide identification and characterization of the HSP gene superfamily in apple snails (Gastropoda: Ampullariidae) and expression analysis under temperature stress

Apple snails from the family Ampullariidae have become economically important due to several species mainly from the genus Pomacea being invasive. The heat shock proteins (HSPs), which are important molecular chaperones for species responding to various stresses, have been proved to play critical roles in adapting harsh environments in the invasive apple snails. The recent release of the genomes of Pomacea canaliculata, Pomacea maculata, Lanistes nyassanus, and Marisa cornuarietis has opened the opportunity for a comprehensive analysis of HSP superfamily in the ampullariids. We identified the number of HSP from P. canaliculata (PcaHSPs) was greater than that from the other three species. A total of 42 PcaHSPs were distributed on 12 chromosomes and were classified into the families of HSP90, HSP70, HSP60, HSP40, HSP20, and HSP10. Each family formed a monoophyletic clade on the phylogenetic tree with strong support values, except for the HSP90 and HSP70 families. The RNA-seq data shows that most the PcaHSPs were of tissue-specific expression levels. Moreover, we identified more HSP genes with stronger transcription levels in response to heat than cold stress. Our findings are inforomative for future studies on stress adaptation and developing effective management strategies focusing on HSPs in invasive apple snails.