Two I84 family protease inhibitors from Chinese razor clams Sinonovacula constricta expressed in response to environmental challenges

Ammonia-induced oxidative stress triggered apoptosis in the razor clam (Sinonovacula constricta)

As one of the most significant contaminants and stressors in aquaculture systems, ammonia adversely jeopardizes the health of aquatic animals. Ammonia exposure affects the development, metabolism, and survival of shellfish. However, the responses of the innate immune and antioxidant systems and apoptosis in shellfish under ammonia stress have rarely been reported. In this study, razor clams (Sinonovacula constricta) were exposed to different concentrations of non-ion ammonia (0.25 mg/L, 2.5 mg/L) for 72 h and then placed in ammonia-free seawater for 72 h for recovery. The immune responses induced by ammonia stress on razor clams were investigated by antioxidant enzyme activities and degree of apoptosis in digestive gland and gill tissues at different time points. The results showed that exposure to a high concentration of ammonia greatly disrupted the antioxidant system of the razor clam by exacerbating the accumulation of reactive oxygen species ( O 2 - , H2O2) and disordering the activities of antioxidant enzymes (superoxide dismutase, catalase, and glutathione peroxidase), and the level of activity remained at a significantly high level after recovering for 72 h (P < 0.05). In addition, there were significant differences (P < 0.05) in the expression of key genes (Caspase 7, Cyt-c, Bcl-2, and Bax) in the mitochondrial apoptotic pathway in the digestive glands and gills of razor clams as a result of ammonia stress and were unable to return to normal levels after 72 h of recovery. TUNEL staining indicated that apoptosis was more pronounced in gills, showing a dose and time-dependent pattern. As to the results, ammonia exposure leads to the activation of innate immunity in razor clams, disrupts the antioxidant system, and activates the mitochondrial pathway of apoptosis. This is important for comprehending the mechanism underlying the aquatic toxicity resulting from ammonia in shellfish.

Genome-Wide Analysis of Family I84 Protease Inhibitor Genes in Three Bivalves Reveals Important Information About the Protein Family's Evolution

Family I84 serine protease inhibitors are believed to be mollusk specific proteins involved in host defense. The molecular evolution of the family, however, remains to be understood. In this study, the genes of Family I84 protease inhibitors in 3 bivalves, Crassostrea gigas, Crassostrea virginica and Tegillarca granosa, were analyzed at the genomic level. A total of 66 Family I84 genes (22 in C. gigas, 28 in C. virginica and 16 in T. granosa) were identified from the 3 species. They distributed unevenly in the genomes involving 4 chromosomes in C. gigas and 5 chromosomes in C. virginica and T. granosa and some genes were tandemly duplicated. Most genes had 3 exons with 12 genes having 4 exons and 1 gene having 2 exons. All genes but 1 from C. gigas and 1 from T. granosa encoded peptides with a signal sequence at the N-terminus, and the properties of the predicted mature molecules were similar. Four conserved motifs were identified in the 66 amino acid sequences. Collinear analysis revealed higher collinearity between the 2 oyster species in general genes and in Family I84 genes. Phylogenetic analysis of the 66 genes with those previously reported from 3 other bivalves and 1 gastropod showed that Family I84 protease inhibitor genes from the same species tended to be grouped together in terminal branches of the constructed Maximum likelihood tree, but most internal nodes were poorly supported by the bootstrap values. In addition, differences in expression patterns between the genes of a same species were observed in the developmental stages and tissues of C. gigas and T. granosa. Moreover, the co-expression of genes within Family I84 and Family I84 genes with non-Family I84 were also detected in C. gigas and T. granosa. These results suggested that Family I84 protease inhibitor genes evolved by active duplications and structural and functional diversifications after the speciation of related mollusks, and the diversified protease inhibitor family was likely multifunctional.

Genome-wide identification and immune response analysis of serine protease inhibitor genes in the blood clam Tegillarca granosa

Serine protease inhibitors (SPIs) are the main regulators of serine protease activities. In this study, we present a genome-wide identification of SPI genes in T. granosa（TgSPI genes）and their expression characteristics in respond to Vibrio stress. A total of 102 TgSPI genes belonging to eight families, including Serpin, TIL (trypsin inhibitor like cysteine rich domain), Kunitz, Kazal, I84, Pacifastin, WAP (whey acidic protein) and A2M (Alpha-2-macroglobulin) were identified, while no genes belonging to Bowman-Birk, amfpi and Antistasin families were identified. The Kazal family has the most TgSPI genes with 38, and 11 TgSPI genes belong to the mollusc-specific I84 family. The TgSPI genes were found to be randomly distributed on 17 chromosomes with 12 tandem duplicate gene pairs. Expression profiles showed that most TgSPI genes were mainly expressed in immune-related tissues such as hepatopancreas, gill and mantle. In the hepatopancreas, most of TgSPI genes were sensitive to Vibrio stress, 28 and 29 TgSPI genes were up-regulated and down-regulated, respectively. Some up-regulated genes with signal peptides, such as the TgSPIs of I84 family, may act as a mechanism to directly prevent Vibrio from invasion. Six Kazal-type TgSPIs (TgSPI29, 45, 49, 50, 51 and 52) were intracellular proteins and their expression was down-regulated in hemocytes after Vibrio stress. This may have boosted protease activity in hemocytes to the point that more hemoglobin derived peptides were produced and secreted into the hemolymph to exert their anti-Vibrio effects. These findings may provide valuable information for further clarifying the roles of SPIs in the immune defense and will benefit future exploration of the immune function of SPIs in molluscs.

Genome Wide Identification and Expression Profiling Indicate Expansion of Family I84 Protease Inhibitor via Gene Tandem Duplication and Divergence in Razor Clam Sinonovacula constricta

Family I84 protease inhibitors represent a novel family in the MEROPS peptidase database and are likely unique for molluscan host defense. Two Family I84 members, scSI-1 and scSI-2, were reported from the razor clam Sinonovacula constricta in a previous research. In the present study, 12 additional genes, named scSI-3 to scSI-14, were identified via genome wide sequence analyses. Among them, 10 genes were predicted to have a signal sequence, but one (scSI-7) was not. Besides, one sequence (scSI-14) was likely to encode a prematurely terminated peptide. The predicted mature peptides shared characteristics including 12 conserved cysteine residues, isoelectric points of 4.98 to 6.11, and molecular weights of 7.1 to 9.3 kDa with previously reported family members. Four motifs were characterized in 13 predicted mature peptides (with exception of scSI-14), which shared two to four conserved cysteine residues, are possibly to form two functional domain comprised 6 cysteine residues, respectively. At genomic level, all the 14 razor clam Family I84 genes were organized into 3 exons and 2 introns; 13 of them clustered in 3 regions of 100 kb on 3 separate chromosomes, suggesting tandem duplications of related genes. The promoter region of all the 14 genes was predicted to share some transcription factor binding sites, in particular those responsive to pathological and physiological stimuli, but no shared motifs were identified. Analyses also revealed differences in expression patterns among the genes. One gene in a tandem duplicated gene pairs usually showed a higher expression level than the other whereas non-tandem duplicated genes exhibited a higher degree of correlation in expression level. In addition, 8 of the 14 genes demonstrated higher level of expression in Vibrio tolerant clams than in non-tolerant clams following challenges with Vibrio parahaemolyticus. These results generated important information about the evolution of Family I84 protease inhibitors in S. constricta.

Family I84 protease inhibitors likely constitute a Mollusca-specific protein family functioning in host defense

Protease inhibitors are proteins or small polypeptides functioning in numerous biological processes in all organisms. The I84 family of protease inhibitors in the MEROPS database represents a novel protease inhibitor family that has been reported in 2 bivalves, Crassostrea virginica and Sinonovacula constricta, and is believed to play a role in host defense. In the present study, 7 new members of Family I84 were identified in 2 bivalves, Meretrix meretrix and Mytilus galloprovincialis, and 1 gastropod, Haliotis discus hannai, at the mRNA level via cDNA cloning. The expression patterns of the newly identified genes varied in response to salinity stresses and pathogen-associated molecular pattern stimulations, suggesting their involvement in the host defense of related species. Additionally, analyses of sequence data in public databases did not reveal any Family I84 protease inhibitor molecules in non-molluscan animals. The results indicated that Family I84 protease inhibitors are likely mollusk specific, constituting a unique host defense mechanism in molluscan species.

Integrated analysis of transcriptomic and metabolomic data to evaluate responses to hypersalinity stress in the gill of the razor clam (Sinonovacula constricta)

Salinity is an important ecological factor that affects physiological metabolism, survival, and distribution of marine organisms. Despite changes in the osmolarity and composition of the cytosol during salinity shifts, marine mollusks are able to maintain their metabolic function. The razor clam (Sinonovacula constricta) survives the wide range of salinity in the intertidal zone via changes in behavior and physiology. To explore the stress responses and mechanisms of salinity tolerance in razor clams, we collected transcriptomic and metabolomic data from a control group (salinity 20‰, S20) and a salinity-stress group (salinity 35‰, S35). The transcriptome data showed that genes related to the immune system, cytoskeleton remodeling, and signal transduction pathways dominated in the S35 group to counteract hypersalinity stress in the gill. The metabolomic analysis showed that 142 metabolites were significantly different between the S35 and S20 groups and that amino acid and carbohydrate metabolism were affected by hypersalinity stress. Levels of amino acids and energy substances, such as l-proline, isoleucine, and fructose, were higher in the gill of the S35 group. The combination of transcriptomic and metabolomic data indicated that metabolism of amino acids, carbohydrates, and lipids was enhanced in the gill during adaptation to high salinity. These results clarified the complex physiological processes involved in the response to hyperosmotic stress and maintenance of metabolism in the gill of razor clams. These findings provide a reference for further study of the biological responses of euryhaline shellfish to hyperosmotic stress and a molecular basis for the search for populations with high salinity tolerance.

The razor clam Sinonovacula constricta uses the strategy of conversion of toxic ammonia to glutamine in response to high environmental ammonia exposure

High ammonia can inhibit the survival and growth, and even cause mortality of razor clam (S. constricta). The accumulation of ammonia to lethal concentrations in some invertebrates may be partially prevented by converting some of the ammonia into glutamine (Gln). Glutamine dehydrogenase (GDH) and glutamine synthetase (GS) have been widely implicated a central role in response to ammonia stress. However, the molecular and physiological response of GDH and GS to ammonia alterations has not yet been determined in clams. To investigate the possible participatory role of GDH and GS genes in altered ammonia conditions, we have cloned their gene sequences and examined the mRNA expression and western blotting under ammonia exposure in S. constricta (ScGDH and ScGS), and detected the levels of GS and GDH, and the content of glutamate (Glu) and Gln. The full-length cDNA of ScGDH was 3924 bp, with a 1629 bp open reading frame (ORF) encoding a 542 amino-acid polypeptide. The complete cDNA sequence for ScGS had 2739 bp with an ORF of 1110 bp encoding 369 amino acids. To investigate ammonia detoxification strategies, the clams were exposed to ammonia for 96 h at four different concentrations (0, 100, 140, and 180 mg/L). Exposure to ammonia resulted in a significant increase of glutamate concentration and Gln in the haemocytes. GDH activity, GDH relative mRNA and protein expression, GS activity, GS relative mRNA and protein expression increased significantly and showed a pronounced time and dosage interaction in the liver. The results suggested that the protective strategies of Gln formation existed in S. constricta, which could convert ammonia to non- or less toxic nitrogenous compounds on the exposure of ammonia. Glutamate content in the haemocytes increased significantly, which is to ensure sufficient Glu to meet the needs for GS to catalyze the conversion of ammonia to Gln. We proposed that the induction of Glu synthesis-related genes and the subsequent formation of the active protein occurred in preparation for the increased capacity of the body to convert ammonia, into Gln. The results of this study suggested that GDH and GS play an important role in the synthesis of Gln, emphasizing, the protective strategies of Gln formation in S. constricta convert ammonia to nontoxic or less toxic nitrogenous compounds upon exposure to ammonia.

Pathogen proteases and host protease inhibitors in molluscan infectious diseases

The development of infectious diseases represents an outcome of dynamic interactions between the disease-producing agent's pathogenicity and the host's self-defense mechanism. Proteases secreted by pathogenic microorganisms and protease inhibitors produced by host species play an important role in the process. This review aimed at summarizing major findings in research on pathogen proteases and host protease inhibitors that had been proposed to be related to the development of mollusk diseases. Metalloproteases and serine proteases respectively belonging to Family M4 and Family S8 of the MEROPS system are among the most studied proteases that may function as virulence factors in mollusk pathogens. On the other hand, a mollusk-specific family (Family I84) of novel serine protease inhibitors and homologues of the tissue inhibitor of metalloprotease have been studied for their potential in the molluscan host defense. In addition, research at the genomic and transcriptomic levels showed that more proteases of pathogens and protease inhibitor of hosts are likely involved in mollusk disease processes. Therefore, the pathological significance of interactions between pathogen proteases and host protease inhibitors in the development of molluscan infectious diseases deserves more research efforts.