Isolation and functional identification of immune cells in hemolymph of blood clams Tegillarca granosa.

Blood clam Tegillarca granosa is a type of economically cultivated bivalve mollusk with red blood, and it primarily relies on hemocytes in its hemolymph for immune defense. However, there are currently no reports on the isolation and identification of immune cells in T. granosa, which hinders our understanding of their immune defense. In this study, we employed single-cell transcriptome sequencing (scRNA-seq) to visualize the molecular profile of hemocytes in T. granosa. Based on differential expression of immune genes and hemoglobin genes, hemocytes can be molecularly classified into immune cells and erythrocytes. In addition, we separated immune cells using density gradient centrifugation and demonstrated their stronger phagocytic capacity compared to erythrocytes, as well as higher levels of ROS and NO. In summary, our experiments involved the isolation and functional identification of immune cells in hemolymph of T. granosa. This study will provide valuable insights into the innate immune system of red-blood mollusks and further deepen the immunological research of mollusks.

Genome-wide identification and characteristic analysis of ETS gene family in blood clam Tegillarca granosa.

BACKGROUND: ETS transcription factors, known as the E26 transformation-specific factors, assume a critical role in the regulation of various vital biological processes in animals, including cell differentiation, the cell cycle, and cell apoptosis. However, their characterization in mollusks is currently lacking. RESULTS: The current study focused on a comprehensive analysis of the ETS genes in blood clam Tegillarca granosa and other mollusk genomes. Our phylogenetic analysis revealed the absence of the SPI and ETV subfamilies in mollusks compared to humans. Additionally, several ETS genes in mollusks were found to lack the PNT domain, potentially resulting in a diminished ability of ETS proteins to bind target genes. Interestingly, the bivalve ETS1 genes exhibited significantly high expression levels during the multicellular proliferation stage and in gill tissues. Furthermore, qRT-PCR results showed that Tg-ETS-14 (ETS1) is upregulated in the high total hemocyte counts (THC) population of T. granosa, suggesting it plays a significant role in stimulating hemocyte proliferation. CONCLUSION: Our study significantly contributes to the comprehension of the evolutionary aspects concerning the ETS gene family, while also providing valuable insights into its role in fostering hemocyte proliferation across mollusks.

Genomic Insights into the Origin and Evolution of Molluscan Red-Bloodedness in the Blood Clam Tegillarca granosa.

Blood clams differ from their molluscan kins by exhibiting a unique red-blood (RB) phenotype; however, the genetic basis and biochemical machinery subserving this evolutionary innovation remain unclear. As a fundamental step toward resolving this mystery, we presented the first chromosome-level genome and comprehensive transcriptomes of the blood clam Tegillarca granosa for an integrated genomic, evolutionary, and functional analyses of clam RB phenotype. We identified blood clam-specific and expanded gene families, as well as gene pathways that are of RB relevant. Clam-specific RB-related hemoglobins (Hbs) showed close phylogenetic relationships with myoglobins (Mbs) of blood clam and other molluscs without the RB phenotype, indicating that clam-specific Hbs were likely evolutionarily derived from the Mb lineage. Strikingly, similar to vertebrate Hbs, blood clam Hbs were present in a form of gene cluster. Despite the convergent evolution of Hb clusters in blood clam and vertebrates, their Hb clusters may have originated from a single ancestral Mb-like gene as evidenced by gene phylogeny and synteny analysis. A full suite of enzyme-encoding genes for heme synthesis was identified in blood clam, with prominent expression in hemolymph and resembling those in vertebrates, suggesting a convergence of both RB-related Hb and heme functions in vertebrates and blood clam. RNA interference experiments confirmed the functional roles of Hbs and key enzyme of heme synthesis in the maintenance of clam RB phenotype. The high-quality genome assembly and comprehensive transcriptomes presented herein serve new genomic resources for the super-diverse phylum Mollusca, and provide deep insights into the origin and evolution of invertebrate RB.

Genome-wide identification and characteristic analysis of PGRP gene family in Tegillarca granosa reveals distinct immune response of the invasive pathogen.

The peptidoglycan recognition proteins (PGRPs) are conserved innate immune molecular in invertebrates and vertebrates, which play important roles in immune system by recognize the peptidoglycans of bacterial cell walls. Although PGRPs have been extensively characterized in insects, a systematic analysis of PGRPs in bivalves is lacking. In the present study, the phylogenic relationships, gene structures and expression profiles of PGRPs in marine bivalves were analyzed. The results indicated that the most PGRPs of bivalves were predicted to degrade the peptidoglycans and prevent excessive immunostimulation of bacteria. In addition, the results of the present study showed that the protein diversity of PGRPs in most marine bivalves was mainly generated by the alternative splicing of genes, however the alternative splicing of PGRP gene family was absent in Tegillarca granosa. The differences of PGRPs might be related to the genetic and environmental differences of marine bivalves. Spatiotemporal expression profiling in T. granosa suggested that PGRPs play important roles in the immune response of invasive pathogens. The present study describes a comprehensive view of PGRPs in the blood clam T. granosa and provides a foundation for functional characterization of this gene family in innate immune of marine bivalves.

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.

Prediction and characterization of a novel hemoglobin-derived mutant peptide (mTgHbP7) from Tegillarca granosa.

The hemoglobin (Hb) is identified in Tegillarca granosa and its derived peptides have been proved to possess antibacterial activity against gram-positive and gram-negative bacteria. In this study, we identified a series of novel antimicrobial peptides (AMPs) and artificially mutated AMPs derived from subunits of T. granosa Hbs, among which, a mutant T. granosa hemoglobin peptide (mTgHbP) mTgHbP7, was proved to possess predominant antibacterial activity against three bacteria strains (Vibrio alginolyticus, V. parahaemolyticus and Escherichia coli). Besides, mTgHbP7 was predicted to form α-helical structure, which was known to be an important feature of bactericidal AMPs. Furthermore, upon contact with HEK293 cell line, we confirmed that mTgHbP7 had no cytotoxicity to mammalian cell even at a high concentration of 160 µM. Therefore, the findings reported here provide a rationalization for antimicrobial peptide prediction and optimization from mollusk hemoglobin, which will be useful for future development of antimicrobial agents.

Evolution of mitochondrial gene arrangements in Arcidae (Bivalvia: Arcida) and their phylogenetic implications.

Arcidae is a diverse group of ark shells with over 260 described species. The phylogenetic relationships and the evolution of the mitochondrial genomes in this family were poorly understood. Comparisons of mitogenomes have been widely used to explore the phylogenetic relationship among animal taxa. We described the complete mitogenomes of Arca navicularis, Scapharca gubernaculum and one nearly complete mitogenome of Anadara consociata. The mitogenome of A. navicularis (18,103 bp) is currently the smallest known Arcidae mitogenome, while the mitogenomes of S. gubernaculum (45,697 bp) and A. consociata (44,034 bp) are relatively large. The mitochondrial gene orders of the three taxa were substantially different from each other, as well as the patterns found in other ark shells. The relationships among Arcidae species recovered from different mitochondrial characters (nucleotide sequence versus gene order) were in disagreement. The phylogeny based on nucleotide sequences did not support the monophyly of Arcidae, as Cucullaea labiata (Cucullaeidae) appeared as a subgroup within Arcinae, rather than sister group to the family Arcidae. In addition, we presented the first time-calibrated evolutionary tree of Arcidae based on mitochondrial DNA (mtDNA) sequences, which placed the deepest divergence within Arcidae at 342.36 million years ago (Mya), around the Carboniferous (360-300 Mya).

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size.

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Functional characterization of a putative lipopolysaccharide-induced TNF-alpha factor (LITAF) from blood clam Tegillarca granosa in innate immunity.

Lipopolysaccharide-induced TNF-alpha factor (LITAF), as a transcription factor, activates the transcription of TNF and other cytokines in inflammatory response upon lipopolysaccharide (LPS) stimulation. In the present study, we cloned and identified the full-length cDNA of LITAF homolog from blood clam Tegillarca granosa for the first time. The full-length cDNA of TgLITAF was 1801 bp encoding a polypeptide of 147 amino acids with an estimated molecular mass of 16.13 kDa. TgLITAF contained a zf-LITAF-like zinc ribbon domain at the C-terminal of the protein and the TgLITAF domain showed 48-74% amino acid sequence identity with other known LITAFs from other species. Subcellular localization study showed that TgLITAF was mainly expressed in the nucleus. qRT-PCR analysis showed that the TgLITAF transcription expressed constitutively in all the examined tissues with the highest expression level in the gills. After LPS or V. alginolyticus treatment, expression of TgLITAF in hemocytes was both up-regulated significantly at 3-6 h. Furthermore, in vitro study indicated that overexpression of TgLITAF in HeLa cells resulted in the activation of TNFα, p53, and influenced the expression levels of apoptotic-related genes Bax, Bcl-2, Caspase-3, Caspase-6, and Caspase-7. The proliferation of HeLa cells was inhibited by overexpression of TgLITAF. Apoptotic fluorescence assay further revealed that TgLITAF participated in the apoptotic process of HeLa cells. Western blotting analysis showed that overexpression of TgLITAF increased endogenous level of cleaved Caspase-7. Taken together, these results revealed that TgLITAF participates in the innate immune response to the pathogen invasion in blood clams and induces apoptosis in HeLa cells.

The transcriptomic responses of the ark shell, Anadara broughtonii, to sulfide and hypoxia exposure.

Sulfide and hypoxia threaten marine organisms in various ways. Anadara broughtonii, a commercial marine bivalve in China which has great potential exposure to sulfide and hypoxia, was selected to test the responses to these stresses. Digital gene expression profile (DGE) analysis was performed on the juveniles' gills after exposed to normal condition (CG group), hypoxia (LO group), and low/high concentration of sulfide (LS/HS group, administered in hypoxia), respectively, using RNA-seq technology. A total of over 30 million clean reads were filtered from each DGE library and over 90% of them were annotated successfully. In total, 774 significant differentially expressed genes (DEGs) were detected and assigned to Gene ontology (GO) classification and KEGG Pathway enrichment analysis. The results show that many of the upregulated DEGs are related to hemoglobin, immunology, and stress responding. In the stressed A. broughtonii, cytochrome P450 and phosphoenolpyruvate carboxykinase may stimulate the glycolysis process to reduce oxygen consumption; Aminoacyl-tRNA synthetases, heat shock protein and protein disulfide isomerase probably help to maintain the genome integrity; Baculoviral IAP repeat-containing protein 2/3, mitogen-activated protein kinase and tumor necrosis factor pathways were probably responsible for protein repair, proteolysis, apoptosis and immune responses to high concentration of sulfide. Combined challenges also induced alternative oxidase and sushi repeat-containing protein, which have indistinct but probably indispensable function in invertebrates. For the first time, comprehensive transcriptome information on A. broughtonii in response to sulfide and hypoxia were provided. Our research offers new insights into the molecular mechanism behind the resistance of shellfish to sulfide and hypoxia.

Immunotoxicity of nanoparticle nTiO2 to a commercial marine bivalve species, Tegillarca granosa.

The increasing production and extensive application of nanoparticles (NPs) inevitably leads to increased release of NPs into the marine environment and therefore poses a potential threat to marine organisms, especially the sessile benthic bivalves. However, the impacts of NPs on the immunity of commercial and ecological important bivalve species, Tegillarca granosa, still remain unknown to date. In addition, the molecular mechanism of the immunotoxicity of NPs still remains unclear in marine invertebrates. Therefore, the immunotoxicity of nTiO2 exposure to T. granosa at environmental realistic concentrations was investigated in the present study. Results obtained showed that the total number, phagocytic activity, and red granulocytes ratio of the haemocytes were significantly reduced after 30 days nTiO2 exposures at the concentrations of 10 and 100 µg/L. Furthermore, the expressions of genes encoding Pattern Recognition Receptors (PPRs) and downstream immune-related molecules were significantly down-regulated by nTiO2 exposures, indicating a reduced sensitivity to pathogen challenges. In conclusion, evident immunotoxicity of nTiO2 to T. granosa at environmental realistic concentrations was detected by the present study. In addition, the gene expression analysis suggests that the PRRs (both TLRs and RIG1 investigated) may be the molecules for NPs recognition in marine invertebrates.

Molecular cloning, characterization and functional analysis of a putative mitogen-activated protein kinase kinase kinase 4 (MEKK4) from blood clam Tegillarca granosa.

The mitogen-activated protein kinase (MAPK) cascades stand for one of the most important signaling mechanisms in response to environmental stimuli. In the present study, we cloned and identified for the first time the full-length cDNA of MAPK kinase kinase 4 (TgMEKK4) from Blood clam Tegillarca granosa using rapid amplification of cDNA ends method. The full-length cDNA of TgMEKK4 was of 1605 bp in length, encoding a polypeptide of 364 amino acids with a predicted molecular mass of 41.22 kDa and theoretical isoelectric point of 6.29. The conserved MEKK4-domain was identified in TgMEKK4 by SMART program analysis. Homology analysis of the deduced amino acid sequence of TgMEKK4 with other known sequences revealed that TgMEKK4 shared 58%-80% identity to MEKK4s from other species. TgMEKK4 mRNA transcripts could be detected in all tissues examined with the highest expression level in the gill by qRT-PCR. The mRNA expression of TgMEKK4 was up-regulated significantly in hemocytes after Vibrio parahaemolyticus, Vibrio alginolyticus and Lipopolysaccharide (LPS) challenges. Overexpression of TgMEKK4 in HEK 293T cells resulted in the activation of JNK and ERK, but not p38. Consistently, In vivo study indicated that LPS stimulation enhanced JNK, ERK and p38 phosphorylation in blood clams. These results suggest that TgMEKK4 is a powerful factor in the regulation of genes that may be involved in innate immune response of blood clam.

Clarifying phylogenetic relationships and the evolutionary history of the bivalve order Arcida (Mollusca: Bivalvia: Pteriomorphia).

The systematics of the bivalve order Arcida constitutes an unresolved conundrum in bivalve systematics. The current definition of Arcida encompasses two superfamilies: Limopsoidea, which includes the recent families Philobryidae and Limopsidae, and Arcoidea, which encompasses the families Arcidae, Cucullaeidae, Noetiidae, Glycymerididae and Parallelodontidae. This classification, however, is controversial particularly with respect to the position and taxonomic status of Glycymerididae. Previous molecular phylogenies were limited either by the use of only a single molecular marker or by including only a few limopsoid and glycymeridid taxa. The challenging nature of Arcida taxonomy and the controversial results of some of the previous studies, prompted us to use a broad range of taxa (55 species), three nuclear markers (18S rRNA, 28S rRNA and histone H3) and a wide range of algorithmic approaches. This broad but stringent approach led to a number of results that differ significantly from previous studies. We provide the first molecular evidence that supports the separation of Arcoidea from Limopsoidea, although the exact position of Glycymerididae remains unresolved, and the monophyly of Limopsoidea is algorithm-dependent. In addition, we present the first time-calibrated evolutionary tree of Arcida relationships, indicating a significant increase in the diversification of arcidan lineages at the beginning of the Cretaceous, around 140Ma. The monophyly of Arcida, which has been supported previously, was confirmed in all our analyses. Although relationships among families remain somehow unresolved we found support for the monophyly of most arcidan families, at least under some analytical conditions (i.e., Glycymerididae, Noetiidae, Philobryidae, and Limopsidae). However, Arcidae, and particularly Arcinae, remain a major source of inconsistency in the current system of Arcida classification and are in dire need of taxonomic revision.

Gal/GalNAc specific multiple lectins in marine bivalve Anadara granosa.

Complete lectin mapping of molluscs with their diversified recognition pattern and possible role in lectin-carbohydrate interaction based immune response triggering need much attention. In this communication, Gal/GalNAc specific three lectins AGL-IA (Anadara granosa lectin-IA), AGL-IB (A. granosa lectin-IB) and AGL-IV (A. granosa lectin-IV) and a lectin having hemolytic activity AGL-III (A. granosa lectin-III) were purified from the plasma of A. granosa bivalve by a combination of gel filtration and affinity chromatography. AGL-IA and IB were oligomeric lectins whereas, AGL-III and IV were monomeric. The molecular weight of AGL-IA, IB, III and IV were 375, 260, 45 and 33 kDa respectively. AGL-IA and IV agglutinated both rabbit and pronase treated human erythrocytes, whereas AGL-IB agglutinated only rabbit erythrocytes. AGL-III was found to agglutinate rabbit erythrocytes, however, it caused hemolysis of pronase treated human erythrocytes. The activity of all four lectins was calcium dependent and maximum at a pH range 7-8. Apart from Gal/GalNAc specific, the four lectins showed substantial differences in their carbohydrate recognition pattern. Moreover, there was a difference in the carbohydrate specificity between AGL-III and other three lectins (AGL-IA, AGL-IB and AGL-IV) towards polyvalent glycotope. On the one hand, 'cluster glycoside effect' i.e., an enhancement of the activity of a multivalent ligand, was observed for carbohydrate specificities of AGL-IA, AGL-IB, AGL-IV. On the other hand, the effect of multivalent ligands on the carbohydrate specificity of AGL-III was opposite of cluster glycoside effect. The affinity of AGL-IA, AGL-IB and AGL-IV for ligands can be ranked as follows: glycoproteins >> polysaccharide > oligosaccharides and monosaccharides. However, Gal related monosaccharides were the best inhibitors of AGL-III and the inhibitory activity decreased gradually in the following order: monosaccharide > disaccharide > polysaccharide. Thus, the diverse specificity of multiple lectins in A. granosa plasma possibly enables to recognize a wide range of microorganisms.

The transcriptomic and biochemical responses of blood clams (Tegillarca granosa) to prolonged intermittent hypoxia.

The blood clam (Tegillarca granosa), a marine bivalve of ecological and economic significance, often encounters intermittent hypoxia in mudflats and aquatic environments. To study the response of blood clam foot to prolonged intermittent hypoxia, the clams were exposed to intermittent hypoxia conditions (0.5 mg/L dissolved oxygen, with a 12-h interval) for 31 days. Initially, transcriptomic analysis was performed, uncovering a total of 698 differentially expressed genes (DEGs), with 236 upregulated and 462 downregulated. These genes show enrichments in signaling pathways related to glucose metabolism, sugar synthesis and responses to oxidative stress. Furthermore, the activity of the enzyme glutathione peroxidase (GPx) and the levels of gpx1 mRNA showed gradual increases, reaching their peak on the 13th day of intermittent hypoxia exposure. This observation suggests an indirect protective role of GPx against oxidative stress. The results of this study make a significantly contribute to our broader comprehensive of the physiological, biochemical responses, and molecular reactions governing the organization of foot muscle tissue in marine bivalves exposed to prolonged intermittent hypoxic conditions.

Polymorphism of the multiple hemoglobins in blood clam Tegillarca granosa and its association with disease resistance to Vibrio parahaemolyticus.

Hemoglobin (Hb) is the major protein component of erythrocytes in animals with red blood, but it can serve additional functions beyond the transport of oxygen. In this study, we identified polymorphism in the blood clam Tegillarca granosa Hb (Tg-Hb) genes and investigated the association of this polymorphism with resistance/susceptibility to Vibrio parahaemolyticus. Analysis of the 540 sequences revealed 28 SNPs in the coding region of three Tg-Hbs, corresponding to about one SNP per 48 bp. Three SNPS: HbIIA-E2-146, HbIIB-E2-23, HbIIB-E2-121 showed a significant association with resistance/susceptibility to V. parahaemolyticus (P < 0.05). To further demonstrate that three significant SNPs of Tg-Hbs is associated with resistance of clams to V. parahaemolyticus, SNPs were genotyped in V. parahaemolyticus resistant strain clams and the wild base population from which this strain was derived. The results indicated that the nonsynonymous mutation T allele at HbIIA-E2-146 and A allele at HbIIB-E2-23 are associated with V. parahaemolyticus resistance in the blood clam, and its association with disease resistance may be due to its cause changes in amino acid sequences to a functional polymorphism. Together with previous bacterial challenge study, these results provides direct evidence that variation at HbIIA-E2-146 and HbIIB-E2-23 are associated with disease resistance in the blood clam, and these two polymorphic loci could be potential gene markers for the future molecular selection of strains that are resistant to diseases caused by V. parahaemolyticus.

Population structure of the blood clam (Tegillarca granosa) in China based on microsatellite markers.

The blood clam, Tegillarca granosa, is widely cultivated in China. We isolated 6 microsatellite loci from T. granosa and used them to investigate genetic diversity and population structure of 5 widely distributed populations of blood clam collected from eastern and southeastern China. The allele number per locus varied from 4 to 9, and the polymorphism information content value was 0.301 to 0.830. The mean observed and expected heterozygosities varied from 0.304 to 0.460 and 0.556 to 0.621, respectively; the population from Yueqing had the smallest observed heterozygosity. In the neighbor-joining tree, Shandong, Fenghua and Yueqing populations clustered together, and there was geographic divergence between Shandong and Guangxi populations. Some microsatellite loci that were isolated from these mainland China samples were not found in blood clams collected from Malaysia.

Expression of glutamine synthetase in Tegillarca granosa (Bivalvia, Arcidae) hemocytes stimulated by Vibrio parahaemolyticus and lipopolysaccharides.

The blood cockle, Tegillarca granosa, is a widely consumed clam in the Indo-Pacific region. Glutamine synthetase (GS) is an enzyme that plays an essential role in the metabolism of nitrogen by catalyzing the condensation of glutamate and ammonia to form glutamine. We identified the GS of T. granosa (Tg-GS) from hemocytes by 3'- and 5'-rapid amplification of cDNA ends (RACE)-PCR. The full-length cDNA consisted of 1762 bp, with a 1104-bp open reading frame encoding 367 amino acids. Sequence comparison showed that Tg-GS has homology to GS of other organisms, with 79.78% identity with GS from the Pacific oyster Crassostrea gigas, 71.98% identity with GS from the zebrafish Danio rerio, and 68.96% identity with human Homo sapiens GS. A C-beta-Grasp domain and an N-catalytic domain were identified in Tg-GS, indicating that Tg-GS should be classified as a new member of the GS family. A quantitative RT-PCR assay was used to detect mRNA expression of Tg-GS in five different tissues. Higher levels of mRNA expression of GS were detected in the tissues of hemocytes and the mantle. Up-regulation of GS by challenge with the bacteria Vibrio parahaemolyticus and with bacterial wall lipopolysaccharides showed that GS plays a role in anti-bacterial immunity. We conclude that pathogen infection significantly induces expression level of Tg- GS, and that activation of GS influences the immune response of T. granosa by increasing glutamine concentration.

Identificaiton and characterization of a novel hemoglobin gene (Tgr-HbIII) from blood clam Tegillarca granosa.

Hemoglobin (Hb) is the major protein component of red blood cells (hemocytes) of the blood clam Tegillarca granosa. Three T. granosa hemoglobin genes have been mentioned in the literature, designated Tgr-HbI, Tgr-HbIIA and Tgr-HbIIB. Previously, our group identified another novel gene, Tgr-HbIII, in the Hb cluster of the chromosome-level genome but the issue of whether this Hb gene expresses functional protein remains unclear. In the current study, phylogenetic analysis revealed that Tgr-HbIII resembles an ancient Hb gene. Sequence alignment and three-dimensional structural modeling results showed that Tgr-HbIII does not bind heme due to the completely different structure at amino acid position 96-100 and replacement of the N100 residue in known Tgr-Hbs with Q100, what causes loss of a single hydrogen bond linking heme with the globin fold. Interface prediction data suggest that Tgr-HbIII forms a homodimer (ΔG = -5.6 kcal/mol) with a similar conformation to the Tgr-HbI homodimer (ΔG = -3.5 kcal/mol). In adult T. granosa, mRNA expression of Tgr-HbIII was lower than that of Tgr-HbIIA and Tgr-HbIIB (up to 100 × ), but comparable to that of Tgr-HbI. Notably, protein expression of Tgr-HbIII was extremely low. Single-cell RNA sequencing analysis of Hb expression showed that all adult hemocytes expressed Tgr-HbI, Tgr-HbIIA and Tgr-HbIIB, while only 43 % (3872 of 8978) expressed Tgr-HbIII. Based on the collective data, we speculate that Tgr-HbIII carried oxygen prior to mutation of N100 to Q100 and subsequently evolved into a known functional remnant of Hb with an adequate mRNA/low protein expression profile. The current study provides a foundation for further research on the origin, evolution and function of molluscan Hbs.

A novel tissue inhibitor of metalloproteinase in blood clam Tegillarca granosa: molecular cloning, tissue distribution and expression analysis.

Tissue inhibitor of metalloproteinases (TIMPs) were originally characterized as inhibitors of matrix metalloproteinases (MMPs), but their range of activities has been found to be broader as it includes the inhibition of several of the MMPs, etc. The cDNA encoding TIMP-4-like gene from blood clam Tegillarca granosa (designated as Tg-TIMP-4-like) which is the first tissue inhibitor of metalloproteinase identified in blood clams, was cloned and characterized. It was of 1164 bp, and an open reading frame (ORF) of 666 bp encoding a putative protein of 222 amino acids. The predicted amino acid sequence comprised all recognized functional domains found in other TIMP homologues and showed the highest (30.56%) identity to the TIMP-1.3 from Crassostrea gigas. Several highly conserved motifs including several TIMP signatures, amino acid residue Cys³° responsible for coordinating the metal ions, the Cys-X-Cys motif and the putative NTR (netrin) domain were almost completely conserved in the deduced amino acid of Tg-TIMP-4 like, which indicated that Tg-TIMP-4-like should be a member of the TIMP family. The mRNA expression of Tg-TIMP-4-like in the tissues of mantle, adductor muscle, foot, gill, hemocyte and hepatopancreas was examined by quantitative real-time PCR (qT-PCR) and mRNA transcripts of Tg-TIMP-4-like were mainly detected in hemocyte, and weakly detected in the other tissues. We also observed that Tg-TIMP-4 like mRNA accumulated significantly during Vibrio parahaemolyticus, Peptidogylcan (PGN) and Lipopolysaccharide (LPS) challenge, whereas the timing and quantitative differences of mRNA expression against different challenge indicated that Tg-TIMP-4-like may play a pivotal role in mollusc defense mechanisms.