Genome-wide association study reveals the underlying regulatory mechanisms of red blood traits in Anadara granosa.

BACKGROUND: Anadara granosa, commonly known as the blood clam, exhibits the unusual characteristic of having red blood among invertebrates. There is significant individual variation in blood color intensity among blood clams; individuals with vibrant red blood are deemed healthier and exhibit stronger stress resistance. However, the molecular basis underlying these red blood traits (RBTs) remains poorly understood. RESULTS: In this study, we performed genome-wide association studies (GWAS) in a population of 300 A. granosa individuals, focusing on RBTs as measured by hemoglobin concentration (HC), total hemocyte count (THC), and heme concentration (HEME). Our analysis identified 18 single nucleotide polymorphisms (SNPs) correlated with RBTs, subsequently selected 117 candidate genes within a 100 kb flanking region of these SNPs, potentially involved in the RBTs of A. granosa. Moreover, we discovered two haplotype blocks specifically associated with THC and HEME. Further analysis revealed eight genes (Septin7, Hox5, Cbfa2t3, Avpr1b, Hhex, Eif2ak3, Glrk, and Rpl35a) that significantly influence RBTs. Notably, a heterozygous A/T mutation in the 3'UTR of Cbfa2t3 was found to promote blood cell proliferation. These genes suggest that the hematopoietic function plays a significant role in the variability of RBTs in A. granosa. CONCLUSIONS: Our findings reveal a conservation of the regulatory mechanisms of RBTs between blood clams and vertebrates. The results not only provide a scientific basis for selective breeding in blood clams, but also offer deeper insights into the evolutionary mechanisms of RBTs in invertebrates.

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.

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.

Expression of ABCA3 transporter gene in Tegillarca granosa and its association with cadmium accumulation.

Exposure to cadmium (Cd), a heavy metal, can cause strong and toxic side effects. Cd can enter the body of organisms in several ways, leading to various pathological reactions in the body. Tegillarca granosa is a kind of bivalve shellfish favored by people in the coastal areas of China. Bivalve shellfish can easily absorb heavy metal pollutants from water bodies while filter feeding. T. granosa is considered a hyper-accumulator of Cd, and the TgABCA3 gene is highly expressed in individuals with a high content of Cd-exposed blood clam. However, it is unclear whether TgABCA3 is involved in Cd ion transport in blood clam and the molecular mechanism for the mechanism of the Cd-induced responses for maintaining cell homeostasis. In this study, the complete cDNA of the TgABCA3 gene was analyzed to provide insights into the roles of TgABCA3 in resistance against Cd in blood clam. The complete sequence of TgABCA3 showed high identity to that of TgABCA3 from other bivalves and contained some classical motifs of ATP-binding cassette transport proteins. TgABCA3 expression in different tissues was measured using real-time quantitative polymerase chain reaction (qRT-PCR) and western blot analysis. The tissue-specific expression showed that TgABCA3 expression was highest in the gill tissue. The TgABCA3 expression in the gill tissue was silenced using the RNA interference technique. After TgABCA3 silencing, the TgABCA3 expression decreased, the Cd content increased, the oxygen consumption and ammonia excretion rates increased, and the ingestion rate decreased. These results showing that the extents of Cd accumulation and resulting toxic effects are related to expression levels and activity of TgABCA3 indicate that TgABCA3 has a protective function against Cd in the clam. This increase in Cd accumulation results in serious damage to the body, leading to the enhancement of its physiological metabolism. Therefore, the findings of the study demonstrated that TgABCA3 can participate in the transport of Cd ions in the blood clam through active transport and play a vital role in Cd detoxification.

Molecular cloning, characterization, and tissue distribution of c-Myc from blood clam Tegillarca granosa and its role in cadmium-induced stress response.

Cadmium (Cd) pollution threatens the cultivation of the blood clam Tegillarca granosa (T. granosa) in coastal regions of the East China Sea. The molecular mechanisms regulating Cd stress response and detoxification in blood clams are largely unclear. In the present study, the full-length T. granosa c-Myc (Tgc-Myc) cDNA was cloned for the first time. The 3063-bp cDNA consisted of a 129-bp 5' untranslated region (UTR), a 1746-bp 3' UTR, and a 1188-bp open reading frame encoding a predicted protein of 395 amino acid residues. The predicted protein had a calculated molecular weight of 44.9 kDa and an estimated isoelectric point of 6.82. The predicted protein contained an N-terminal transactivation domain and a C-terminal basic helix-loop-helix leucine zipper domain, which are conserved functional domains of c-Myc proteins. Tgc-Myc showed broad tissue distribution in blood clams, with the highest expression detected in the gill and hepatopancreas. Exposure to Cd, a major heavy metal pollutant in coastal regions of the East China Sea, induced Tgc-Myc expression in gill tissues. Tgc-Myc knockdown led to reduced expression of a variety of stress response/detoxification genes in blood clams cultivated in Cd-contaminated seawater. Tgc-Myc knockdown also led to decreased expression of IGF1R, a proto-oncogene that promotes cell proliferation. These findings indicated that Tgc-Myc regulates Cd-induced stress response and detoxification in blood clams. The upregulation of Tgc-Myc may serve as an approach to generate strains with an enhanced detoxification response and consequently a low heavy metal buildup.

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.

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.

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.

Hemoglobin of the bloody clam Tegillarca granosa (Tg-HbI) is involved in the immune response against bacterial infection.

Hemoglobins (Hb) are the major protein components of erythrocytes circulating in the red blood, but can serve additional functions besides the transport of oxygen. Here, the cDNA of the bloody clam (Tegillarca granosa) Hb dimer (designated Tg-HbI) was cloned and was found to be 748 bp in length, consisting of an open reading frame of 441 bp encoding a polypeptide of 147 amino acids. The deduced amino acid sequence of Tg-HbI shared 81.6% similarity with HbI from two species of the genus Scapharca and 46-51% similarity with the Hb proteins from other mollusks. The 3D structure of bloody clam Tg-HbI was predicted by the SWISS-MODEL Protein Modelling Server and compared with that of Scapharca kagoshimensis. The mRNA transcript of Tg-HbI was detected in all of the clam cells/tissues examined, including haemocytes, the adductor muscle, foot, hepatopancreas, gill and mantle. The mRNA expression of Tg-HbI was significantly up-regulated after Vibrio parahaemolyticus, lipopolysaccharide and peptidoglycan challenge, indicating that Tg-HbI was involved in the immune defence responses against bacterial infection and exposure to bacterial pathogenic factors. As the first functional research on the Hb protein in bloody clam, our findings provide new insight into the innate immune defence mechanisms of T. granosa and other mollusks.

Identification and characterization of a Tegillarca granosa ferritin regulated by iron ion exposure and thermal stress.

Ferritin, a conserved iron storage protein of most living organisms, plays a crucial role in iron metabolism. In this study, the ferritin gene from Tegillarca granosa (denoted as TgFER) was identified by expressed sequence tag (EST) and PCR approaches. The full-length cDNA of TgFER was of 895bp, consisting of a 5'-UTR of 163 bp with a putative iron regulatory element (IRE), a 3'-UTR of 213 bp, and a complete open reading frame of 519 bp encoding a polypeptide with 172 amino acid residues. The predicted molecular mass of deduced amino acid of TgFER was 20.00 kDa and the theoretical pI was 4.89. The deduced amino acid of TgFER shared high identities to ferritin from abalone, oyster, clam and human. The tissue distribution of TgFER in the tissues of mantle, foot, gill, haemocytes and hepatopancreas was examined by quantitative real-time PCR (q-PCR) and mRNA transcripts of TgFER were found to be dominately expressed in haemocytes, hepatopancreas and gill and weakly in foot and mantle. The temporal expression of TgFER in haemocytes or hepatopancreases after challenged by metals ion (FeCl2 and FeCl3) exposure and thermal stress were also analyzed with q-PCR. The diverse expression patterns of TgFER were detected depending upon the types of stimulators and tissues. The ployconal antibodies generated from the recombinant product of TgFER could be specifically identified not only the recombinant product, but also the native protein from haemocytes. All these results strongly suggested that TgFER was involved in the iron metabolism and thermal stress regulation in T. granosa.