Identification of a Serum amyloid A gene and the association of SNPs with Vibrio-resistance and growth traits in the clam Meretrix meretrix

Characterization of the lipidomic profile of clam Meretrix petechialis in response to Vibrio parahaemolyticus infection

Vibrio parahaemolyticus is a devastating pathogen of clam Meretrix petechialis, which brings about huge economic losses in aquaculture breeding industry. In our previous study, we have found that Vibrio infection is closely associated with lipid metabolism of clams. In this study, an untargeted lipidomics approach was used to explore the lipid profiling changes upon Vibrio infection. The results demonstrated that the hepatopancreas of clams was composed of five lipid categories including fatty acyls, glycerolipids, glycerophospholipids, sphingolipids and sterol lipids. And the content of lipid classes altered during Vibrio infection, implying that Vibrio infection altered intracellular lipid homeostasis in clams. Meanwhile, a total of 200 lipid species including 82 up-regulated and 118 down-regulated significantly were identified in response to Vibrio infection, of which ceramide (Cer), phosphatidylcholine (PC) and triglyceride (TG) accounted for the largest proportion. Notably, all Cers showed a significantly decreased trend while nearly all TG species were increased significantly during Vibrio infection, which suggested that Cer and TG could be determined as effective biomarkers. Furthermore, these differentially expressed lipid species were enriched in 20 metabolic pathways and sphingolipid metabolism was one of the most enriched pathways. These results evidenced how the lipid metabolism altered in the process of Vibrio infection and opened a new perspective on the response of marine bivalves to pathogen infection.

Characterization of a novel activating protein-1 (AP-1) gene and the association of its single nucleotide polymorphisms with vibrio resistance in Tegillarca granosa

The blood clam Tegillarca granosa is a commercial marine bivalve of economic value, accounting for approximately 50% of clam production in China. In recent years, the yield of blood clams has been threatened by bacterial infections caused by marine Vibrio species that thrive under a rising sea temperature. The transcription factor activating protein-1 (AP-1) is emerging as an important player in the innate immunity of marine bivalves against viral or bacterial infections. In this study, the full-length cDNA of a novel T. granosa AP-1 (TgAP-1) was cloned for the first time. The 1591-bp cDNA encoded a protein of 292 amino acid residues with a calculated molecular weight of 32.8 kDa. The TgAP-1 protein contained an N-terminal Jun domain and a C-terminal basic region leucine zipper domain typically found in Jun proteins (a subfamily of AP-1 proteins). TgAP-1 was ubiquitously expressed in T. granosa, with the highest expression detected in the gill and foot, followed by the mantle, hemolymph, and hepatopancreas. Exposure to Vibrio harveyi induced TgAP-1 expression in gill tissues and the expression levels of TgAP-1 of resistant blood clams were always lower than that of control population whether Vibro infection or not. A total of 18 single nucleotide polymorphisms (SNPs) of TgAP-1 were detected in T. granosa. SNP-typing and haplotyping of resistant and susceptible populations revealed that six SNPs (AG type of TgSNP-1, GA type of TgSNP-2, TG type of TgSNP-4, CT type of TgSNP-7, AG type of TgSNP-11, and GA type of TgSNP-12) and four haplotypes (fHap2, fHap3, fHap6, and fHap7) were significantly associated with V. harveyi resistance. Risk assessment showed that fHap2 (CG) and fHap7 (GA) were associated with an increased resistance, while fHap3 (CT) and fHap6 (AG) were associated with an increased susceptibility. The results from this study supported a potential role of TgAp-1 in the anti-Vibro immunity of T. granosa. The discovery of the genetic molecular markers and haplotypes related to Vibrio resistance can provide guidance for selective breeding of T. granosa in the future.

Identification and Expression Characterization of the Smad3 Gene and SNPs Associated with Growth Traits in the Hard Clam (Meretrix meretrix)

It has been demonstrated that the sekelsky mothers against decapentaplegic homolog 3 (Smad3) plays an important role in the growth and development of vertebrates. However, little is known about the association between the Smad3 gene and the growth traits of mollusks. In this study, Smad3 from the hard clam Meretrix meretrix (Mm-Smad3) was cloned, characterized, and screened for growth-related single nucleotide polymorphisms (SNPs) in its exons. The full-length cDNA of Mm-Smad3 was 1938 bp, encoding a protein with 428 amino acid residues. The protein sequence included an MH1 (27–135 aa) and MH2 domain (233–404 aa). Promoter analysis showed that the promoter sequence of Mm-Smad3 was 2548 bp, and the binding sites of Pit-1a, Antp, Hb, and other transcription factors are related to the growth and development of hard clams. The phylogenetic tree was divided into two major clusters, including mollusks and vertebrate. The expression level of Mm-Smad3 was predominantly detected in the mantle and foot, while extremely less expression was observed in the digestive gland. The low expression level of Mm-Smad3 was detected at the stages of unfertilized mature eggs, fertilized eggs, four-cell embryos, blastula, gastrulae, trochophore, and D-shaped larvae, whereas an opposite trend was observed regarding the highest expression at the umbo larvae stage (p < 0.05). In the mantle repair experiment, the time-course expression profiles showed that compared to the expression level at 0 h, Mm-Smad3 significantly decreased at 6 h (p < 0.05) but increased at 12 and 48 h. Further, the association analysis identified 11 SNPs in the exons of Mm-Smad3, of which three loci (c.597 C > T, c.660 C > T, c.792 A > T) were significantly related to the growth traits of clam (p < 0.05). Overall, our findings indicated that Mm-Smad3 is a growth-related gene and the detected SNP sites provide growth-related markers for molecular marker-assisted breeding of this species.

TAF5L functions as transcriptional coactivator of MITF involved in the immune response of the clam Meretrix petechialis

TAF5L is a component of the P300/CBP-associated factor (PCAF) histone acetylase complex, which serves as a coactivator and takes part in basal transcription such as promoter recognition, complex assembly and transcription initiation. In our study, the full-length sequence of MpTAF5L was identified and characterized in the clam M. petechialis. Sequence analysis showed that the predicted MpTAF5L protein had a N-terminal TAF5-NTD2 domain and a C-terminal WD40-repeats domain. The annotation and evolutionary analysis revealed MpTAF5L had close evolutionary relationship with other invertebrate species. Tissue distribution analysis of TAF5L claimed that it was highly expressed in the mantle, adductor muscle, foot and hepatopancreas. The mRNA expression of MpTAF5L was significantly up-regulated after Vibrio parahaemolyticus challenge, indicating its involvement in the immune response of clam. Yeast two-hybrid assays verified that MpTAF5L can interact with MpMITF (a critical immune-related transcription factor), and our further research clarified this interaction depended upon the N-terminal TAF5-NTD2 domain of MpTAF5L. Moreover, the mRNA expression of MpBcl-2 (a target gene of MITF) was significantly decreased but the mRNA expression of MpMITF was not significantly changed after knockdown of MpTAF5L, which indicated the reduction of MpMITF regulating activity at the same time. These results revealed that MpTAF5L interacted with MpMITF and enhanced the activation of MpMITF, which plays roles in the immune defense against V. parahaemolyticus.

MITF Regulates Downstream Genes in Response to Vibrio parahaemolyticus Infection in the Clam Meretrix Petechialis

The microphthalmia-associated transcription factor (MITF) is a basic helix-loop-helix-leucine zipper protein that plays a key role in cell proliferation, survival and immune defense through the direct transcriptional control of downstream genes. We have found that MITF participates in the immune response to Vibrio parahaemolyticus infection in the clam Meretrix petechialis. In this study, we focused on how MITF functions in immunity. First, PO, CTSK, and BCL-2 were identified as the target genes of MpMITF in the clam by RNAi. EMSAs showed direct binding between the MpMITF protein and the E-box of the MpPO, MpCTSK, and MpBCL-2 promoters. Yeast one-hybrid assays also suggested that MpMITF could activate the expression of these three downstream genes. These results demonstrated that the transcriptional expression of MpPO, MpCTSK, and MpBCL-2 is directly regulated by MpMITF. Second, we analyzed the roles of MpPO, MpCTSK, and MpBCL-2 in clam immunity. The mRNA expression of MpPO, MpCTSK, and MpBCL-2 increased significantly after V. parahaemolyticus challenge, which implied that these genes might take part in the immune defense against V. parahaemolyticus challenge in clams. The purified recombinant proteins, MpPO and MpCTSK, inhibited the growth of V. parahaemolyticus. Additionally, the apoptosis rate of clam haemocytes rose significantly when the activity of MpBCL-2 was suppressed. These results revealed that MpPO, MpCTSK, and MpBCL-2 are involved in the immune defense against V. parahaemolyticus. This study supports the idea that the MpMITF pathway plays a key role in immune defense through the direct regulation of the downstream genes MpPO, MpCTSK, and MpBCL-2 in the clam, M. petechialis.

Association and Functional Analyses Revealed That PPP1R3B Plays an Important Role in the Regulation of Glycogen Content in the Pacific Oyster Crassostrea gigas

The Pacific oyster (Crassostrea gigas) is one of the most important aquaculture species worldwide. Glycogen contributes greatly to the special taste and creamy white color of oysters. Previous genome-wide association studies (GWAS) identified several single nucleotide polymorphism (SNP) sites that were strongly related to glycogen content. Genes within 100 kb upstream and downstream of the associated SNPs were screened. One gene annotated as protein phosphatase 1 regulatory subunit 3B (PPP1R3B), which can promote glycogen synthesis together with protein phosphatase 1 catalytic subunit (PPP1C) in mammals, was selected as a candidate gene in this study. First, full-length CgPPP1R3B was cloned and its function was characterized. The gene expression profiles of CgPPP1R3B in different tissues and seasons showed a close relationship to glycogen content. RNA interference (RNAi) experiments of this gene in vivo showed that decreased CgPPP1R3B levels resulted in lower glycogen contents in the experimental group than in the control group. Co-immunoprecipitation (Co-IP) and yeast two-hybrid (Y2H) assays indicated that CgPPP1R3B can interact with CgPPP1C, glycogen synthase (CgGS) and glycogen phosphorylase (CgGP), thus participating in glycogen metabolism. Co-sedimentation analysis in vitro demonstrated that the CgPPP1R3B protein can bind to glycogen molecules directly, and these results indicated the conserved function of the CgPPP1R3B protein compared to that of mammals. In addition, thirteen SNPs were precisely mapped in this gene. Ten of the thirteen SNPs were confirmed to be significantly (p < 0.05) related to glycogen content in an independent wild population (n = 288). The CgPPP1R3B levels in oysters with high glycogen content were significantly higher than those of oysters with low glycogen content, and gene expression levels were significantly associated with various genotypes of four associated SNPs (p < 0.05). The data indicated that the associated SNPs may control glycogen content by regulating CgPPP1R3B expression. These results suggest that CgPPP1R3B is an important gene for glycogen metabolic regulation and that the associated SNPs of this gene are potential markers for oyster molecular breeding for increased glycogen content.

Identification of an MITF gene and its polymorphisms associated with the Vibrio resistance trait in the clam Meretrix petechialis

Microphthalmia-associated transcription factor (MITF) regulates the transcription of its target genes by binding to their promoters. In this study, an MITF gene, MpMITF was identified in the clam Meretrix petechialis. The full-length cDNA of MpMITF is 3564 bp with an ORF of 1365 bp. The deduced amino acid sequence consists of a conserved functional structure of bHLH-LZ, which could bind with E-box. The mRNA and protein expression levels of MpMITF were significantly up-regulated 6 h post-Vibrio injection. The mRNA expression of MpMITF increased on day 2 and peaked on day 10 post-Vibrio immersion. Furthermore, MpMITF expression was significantly up-regulated in most resistant families of clams (P < 0.05) but did not change significantly in most susceptive families of clams after the Vibrio immersion challenge. These results suggest that, in clams, MpMITF participates in the immune response against a Vibrio infection. Genotyping in two clam groups with different resistant levels to Vibrio parahaemolyticus (i.e., 11-R and 11-S), thirteen SNPs and five haplotypes were detected in the DNA sequence of MpMITF, of which five SNPs and two haplotypes were associated with Vibrio resistance. Four SNPs (SNP2, 5, 6 and 13) and one haplotype (Hap1) were further confirmed to be associated with Vibrio resistance in M. petechialis by association analysis in different clam families. This study deepens the understanding of MITF in marine bivalves and provides potential candidate markers for resistance selection in the clam M. petechialis.

Transcriptome profiles of the clam Meretrix petechialis hepatopancreas in response to Vibrio infection

Microbial diseases have received much attention due to their enormous destruction of aquaculture, and Vibrio parahaemolyticus is one of the main pathogens that cause bacterial disease in the clam Meretrix petechialis. To better understand the molecular mechanisms of the immune response to Vibrio in M. petechialis, RNA-Seq was applied to explore global expression changes of hepatopancreas from this clam after Vibrio challenge. There were 199,318,966 clean reads obtained by Illumina sequencing, which were further assembled into 214,577 transcripts, and then 147,255 unigenes with an N50 of 1393 bp were identified. Gene ontology (GO) analysis revealed 21 biological process subcategories, 15 cellular component subcategories and 12 molecular function subcategories. A total of 8358 unigenes were mapped onto 267 biological signaling pathways by KEGG, among which there were 16 pathways related to the immune system. In total, 206 differentially expressed genes (DEGs) were identified, including 113 up-regulated unigenes and 93 down-regulated unigenes. In these DEGs, 96 DEGs were annotated in at least one database, accounting for 46.60% of all significant DEGs. To validate the transcriptome dataset, 15 DEGs were selected for real-time qPCR confirmation and the results showed that expression patterns of 13 genes (86.7%) agreed well with the RNA-Seq analysis. Fourteen of the 206 DEGs were annotated to be immune-related genes, and we examined the expression patterns of four immune-related DEGs using clams post immersion challenge. This study enriched the M. petechialis transcriptome database and provided insight into the immune response of M. petechialis against Vibrio infection.

The polymorphisms of a MIF gene and their association with Vibrio resistance in the clam Meretrix meretrix

Macrophage migration inhibitory factor (MIF) is an important proinflammatory cytokine that mediates both innate and adaptive immune responses. In this study, a homolog of MIF was identified in the clam Meretrix meretrix. Ten SNPs in the DNA partial sequence of MmMIF were found to be significantly associated with Vibrio resistance (P < 0.05). Distinct expression patterns of MmMIF among different haplotypes were observed after Vibrio challenge. The results showed that haplotypes did not affect MmMIF expression in the negative control group, while the expression of MmMIF in clams with Hap1 and Hap1/Hap2 was significantly lower than that with Hap2 at 24 h in the PBS-injected group but significantly higher than that with Hap2 in the Vibrio-injected group. The results also indicate that Hap1 and Hap1/Hap2 can specifically respond to mechanical stimulation while Hap2 can specifically respond to Vibrio infection. The effect of a missense mutation was detected by site-directed mutagenesis using fusion expression of the protein, which showed that the SNP g.737 (I > V) has no effect on redox activity and tautomerase activity. These studies identified a potential marker that is enriched in Vibrio-resistant clams that can be used for the genetic breeding of Meretrix meretrix.

Serum amyloid A in marine bivalves: An acute phase and innate immunity protein

Serum amyloid A (SAA) is among the most potent acute phase proteins (APP) in vertebrates. After injury, its early expression can dramatically increase to promote the recruitment of immuno-competent cells, expression of pro-inflammatory proteins and the activation of the innate immune defences. Although APP have been studied in many vertebrates, only recently their search was extended to invertebrates and the finding of SAA-like molecules has opened new questions on the immune-regulatory functions of these soluble proteins in the animal kingdom. Taking advantage of the considerable amount of genomic and transcriptomic data currently available, we retrieved 51 SAA-like proteins in several protostome taxa comprising 21 marine bivalve species and basal metazoans. In addition to vertebrate-like SAAs, we identified a second protein type with peculiar features. In the bivalves Crassostrea gigas and Mytilus galloprovincialis, both digital expression analysis and qPCR data indicated an induction of the classical SAA after bacterial challenge.