Genome-wide characterization of the cytosolic sulfotransferase 1B member 1 (SULT1B1) family and its expression responses to sulfide stress in the razor clam Sinonovacula constricta

The razor clam (Sinonovacula constricta), a typical burrowing organism in the intertidal zones, is often exposed to sulfide environment and shows strong sulfide tolerance. Located downstream of the sulfur metabolism pathway, cytosolic sulfotransferase family 1B member 1 (SULT1B1) is a key enzyme catalysing the sulfonation reaction, and plays an important role in the biotransformation of endogenous substances such as thyroid hormones (THs). To investigate their roles in sulfide resistance, a systematic analysis of S. constricta SULT1B1s (ScSULT1B1s), including genomic distribution, phylogenetic relationships, gene structure, conserved motifs, and expression profiles under sulfide stress, was performed. A total of 10 ScSULT1B1 genes were found in the S. constricta genome. Sequence analysis showed that ScSULT1B1 gene family encoded 155-425 amino acids, containing four catalytic active sites (K, N, H, and S), one PAPS binding domain at the N-terminus, and one PAPS binding and dimerization domain at the C-terminus. The spatial-temporal expression patterns of ScSULT1B1s were further estimated by quantitative real-time PCR (qRT-PCR). Among them, partial ScSULT1B1s showed significantly high expression in the gill, hepatopancreas, and siphon. Furthermore, the response expression of certain ScSULT1B1s significantly fluctuated under sulfide stress. Together, our results suggest that ScSULT1B1s, by mediating the sulfonation reaction, may regulate THs levels to maintain basic metabolic and immune functions, making S. constricta highly sulfide tolerant.