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

Toxic effects of nanoplastics on biological nitrogen removal in constructed wetlands: Evidence from iron utilization and metabolism

Nanoplastics (NPs) in wastewaters may present a potential threat to biological nitrogen removal in constructed wetlands (CWs). Iron ions are pivotal in microbially mediated nitrogen metabolism, however, explicit evidence demonstrating the impact of NPs on nitrogen removal regulated by iron utilization and metabolism remains unclear. Here, we investigated how NPs disturb intracellular iron homeostasis, consequently interfering with the coupling mechanism between iron utilization and nitrogen metabolism in CWs. Results indicated that microorganisms affected by NPs developed a siderophore-mediated iron acquisition mechanism to compensate for iron loss. This deficiency resulted from NPs internalization limited the activity of the electron transport system and key enzymes involved in nitrogen metabolism. Microbial network analysis further suggested that NPs exposure could potentially trigger destabilization in microbial networks and impair effective microbial communication, and ultimately inhibit nitrogen metabolism. These adverse effects, accompanied by the dominance of Fe3+ over certain electron acceptors engaged in nitrogen metabolism under NPs exposure, were potentially responsible for the observed significant deterioration in nitrogen removal (decreased by 30 %). This study sheds light on the potential impact of NPs on intracellular iron utilization and offers a substantial understanding of the iron-nitrogen coupling mechanisms in CWs.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

NO production of granulocytes associated with antibacterial immune response in Tegillarca granosa

Nitric oxide (NO) is a signaling molecule and immune effector produced by the nitric oxide synthases (NOS), which involved to various physiological processes of animals. In marine bivalves, hemocytes play important roles in antimicrobial innate immune response. Although hemocyte-derived NO has been detected in several bivalves, the immune function of hemocyte-derived NO is not well understood. Here, we investigated the antibacterial response of hemocyte-derived NO in the blood clam Tegillarca granosa. Two types of hemocytes including erythrocytes and granulocytes were isolated by Percoll density gradient centrifugation, their NO production and TgNOS expression level were analyzed. The results showed that NO was mainly produced in granulocytes and almost no detected in erythrocytes. The granulocytes showed significantly higher NO level and TgNOS expression level than the erythrocytes. And the TgNOS expression level was significantly increased in granulocytes after Vibro parahemolyticus challenge. In addition, the NO donor sodium nitroprusside (SNP) significantly increased the NO production of hemocytes to kill pathogenic bacteria. In summary, the results revealed that granulocytes-derived NO play vital roles in the antimicrobial immune response of the blood clam.

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.

Hemoglobin wonders: a fascinating gas transporter dive into molluscs

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Taxonomic Distribution and Molecular Evolution of Mytilectins

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a β-trefoil fold. Mytilectins represent a recently described subgroup of β-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the β-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

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.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

Chromosome-level genome assembly of sea scallop Placopecten magellanicus provides insights into the genetic characteristics and adaptive evolution of large scallops

Placopecten magellanicus (Gmelin, 1791), a deep-sea Atlantic scallop, holds significant commercial value as a benthic marine bivalve along the northwest Atlantic coast. Recognizing its economic importance, the need to reconstruct its genome assembly becomes apparent, fostering insights into natural resources and generic breeding potential. This study reports a high-quality chromosome-level genome of P. magellanicus, achieved through the integration of Illumina short read sequencing, PacBio HiFi sequencing, and Hi-C sequencing techniques. The resulting assembly spans 1778 Mb with a scaffold N50 of 86.71 Mb. An intriguing observation arises - the genome size of P. magellanicus surpasses that of its Pectinidae family peers by 1.80 to 2.46 times. Within this genome, 28,111 protein-coding genes were identified. Comparative genomic analysis involving five scallop species unveils the critical determinant of this expanded genome: the proliferation of repetitive sequences recently inserted, contributing to its enlarged size. The landscape of whole genome collinearity sheds light on the relationships among scallop species, enhancing our broader understanding of their genomic framework. This genome provides genomic resources for future molecular biology research on scallops and serves as a guide for the exploration of longevity-related genes in scallops.

A C-type lectin with a single carbohydrate-recognition domain (CRD) containing unique QPN/WDD motifs from Tegillarca granosa is involved in the innate immune defense

C-type lectins (CTLs), a superfamily of Ca2+-dependent carbohydrate-recognition proteins, serve as pattern recognition receptors (PRRs) in the immune response of many species. However, little is currently known about the CTLs of the commercially and ecologically important bivalve species, blood clam (Tegillarca granosa). In this study, a CTL (designated as TgCTL-1) with a single carbohydrate-recognition domain (CRD) containing unique QPN/WDD motifs was identified in the blood clam through transcriptome and whole-genome searching. Multiple alignment and phylogenetic analysis strongly suggested that TgCTL-1 was a new member of the CTL superfamily. Expression analysis demonstrated that TgCTL-1 was highly expressed in the hemocytes and visceral mass of the clam under normal condition. In addition, the expression of TgCTL-1 was shown to be significantly up-regulated upon pathogen challenge. Moreover, the recombinant TgCTL-1 (rTgCTL-1) displayed agglutinating and binding activities against both the gram-positive and gram-negative bacteria tested in a Ca2+-dependent manner. Furthermore, it was found that the in vitro phagocytic activity of hemocytes was significantly enhanced by rTgCTL-1. In general, our results showed that TgCTL-1 was an inducible acute-phase secretory protein, playing crucial roles in recognizing, agglutinating, and binding to pathogenic bacteria as well as modulating phagocytic activity of hemocytes in the innate immune defense of blood clam.

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 N¹⁰⁰ residue in known Tgr-Hbs with Q¹⁰⁰, 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 N¹⁰⁰ to Q¹⁰⁰ 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.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

Immune responses of hemocytes in the blood clam Tegillarca granosa in response to in vivo Vibrio harveyi infection

Aquaculture of the blood clam Tegillarca granosa accounts for approximately 50% of Arcidae (ark shell) production in China. Vibrio infection severely threatens the sustainability of the clam aquaculture industry. Exposure to Vibrio induces an immune response in blood clams. However, the underlying mechanism remains poorly understood. In this study, immune responses of hemocytes in blood clams were detected after Vibrio infection; the immersion method was used in vivo to mimic the clam's natural infection process. After 24 h of exposure to Vibrio infection, the Vibrio load in hemolymph fluid in both the treatment Ⅰ (25,033.33 ± 19,563.11 CFU/mL) and treatment Ⅱ (122,163.33 ± 194,409.49 CFU/mL) groups were significantly higher, than that in the control group (13.67 ± 37.73 CFU/mL) (P < 0.05). Correspondingly, the production of intracellular reactive oxygen species was approximately 1.40 (treatment Ⅰ) and 2.12 (treatment Ⅱ) fold higher than that in the control group (P < 0.05), and the induced DNA damage showed a similar trend (P < 0.05). Vibrio infection also significantly increased lysozyme content, adenosine triphosphate content, and peroxidase isozyme activity, in both the serum and hemocyte lysates (P < 0.05). The expression of immune-associated genes (ABCA3, c-Myc, Caspase 3, and HSP70) was upregulated under infection conditions. The phagocytic activity was approximately 1.99 (treatment Ⅰ) and 2.57 (treatment Ⅱ) fold that in control clams (P < 0.05). In addition, the total hemocyte count and red granulocyte percentage both significantly decreased by approximately 75-90% after Vibrio infection. These results provided novel insights into the mechanism of hemocyte immunity in T. granosa against Vibrio infection, which may aid in the future prevention and control of Vibrio infection in vivo.

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.

A high-quality chromosome-level genome assembly of the bivalve mollusk Mactra veneriformis

Mactra veneriformis (Bivalvia: Mactridae) is a bivalve mollusk of major economic importance in China. Decreased natural yields of M. veneriformis have led to an urgent need for genomic resources. To address this problem and the currently limited knowledge of molecular evolution in this genus, we here report a high-quality chromosome-level genome assembly of M. veneriformis. Our approach yielded a 939.32 Mb assembled genome with an N50 contig length of 7,977.84 kb. Hi-C scaffolding of the genome resulted in assembly of 19 pseudochromosomes. Repetitive elements made up ∼51.79% of the genome assembly. A total of 29,315 protein-coding genes (PCGs) were predicted in M. veneriformis. Construction of a genome-level phylogenetic tree demonstrated that M. veneriformis and Ruditapes philippinarum diverged around 231 million years ago (MYA). Inter-species comparisons revealed that 493 gene families have undergone expansion and 449 have undergone contraction in the M. veneriformis genome. Chromosome-based macrosynteny analysis revealed a high degree of synteny between the 19 chromosomes of M. veneriformis and those of Patinopecten yessoensis. These results suggested that M. veneriformis has a similar karyotype to that of P. yessoensis, and that a highly conserved 19-chromosome karyotype was formed in the early differentiation stages of bivalves. In summary, the genomic resources generated in this work serve as a valuable reference for investigating the molecular mechanisms underlying biological functions in M. veneriformis and will facilitate future genetic improvement and disease treatment in this economically important species. Furthermore, the assembled genome greatly improves our understanding of early genomic evolution of the Bivalvia.

Hypoxia-mediated immunotoxicity in the blood clam Tegillarca granosa

In marine ecosystems, dissolved oxygen (DO) is essential for maintaining intracellular energy balance during aerobic metabolism. Bivalve mollusks are frequently exposed to hypoxia environments due to tides, temperature changes, and anthropogenic activities. The blood clam, Tegillarca granosa, mainly inhabits intertidal mudflats and is more susceptible to low oxygen events. In this study, we investigated the effect of hypoxia on immune responses in clams, and showed that hypoxia exposure reduced total hemocyte counts (THC), hemoglobin concentrations, and intracellular reactive oxygen species (ROS) levels. Also, phagocytic and cell activities of hemocyte were significantly inhibited. Furthermore, immune-related gene expression was also down-regulated. In conclusion, hypoxia greatly affected immune functions in blood clams, and our research provided the foundation for further mechanistic studies on hypoxia tolerance in clams.

Identification, characterization, and antimicrobial activity of a novel big defensin discovered in a commercial bivalve mollusc, Tegillarca granosa

Molluscs, the second largest animal phylum on earth, primarily rely on cellular and humoral immune responses to fight against pathogen infection. Although antimicrobial peptides (AMPs) such as big defensin play crucial roles in the humoral immune response, it remains largely unknown in the ecological and economic important blood clam (Tegillarca granosa). In this study, a novel big defensin gene (TgBD) was identified in T. granosa through transcripts and whole genome searching. Bioinformatic analyses were conducted to explore the molecular characteristics of TgBD, and comparisons of TgBD with those reported in other molluscs were performed by multiple alignments and phylogenetic analysis. In addition, the expression patterns of TgBD in various tissues and upon bacterial challenge were investigated while the antimicrobial activity of synthetic N-terminal domain of TgBD was confirmed in vitro by radial diffusion experiment. Results obtained showed TgBD had an open reading frame (ORF) of 369 bp, encoding a prepropeptide containing a signal peptide and a propeptide. Similar to big defensins reported in other species, TgBD consists of a hydrophobic N-terminal domain containing β1-α1-α2-β2 folds and a cysteine-rich cationic C-terminal domain with three disulfide bonds between C1-C5, C2-C4, and C3-C6. Phylogenetic analysis showed that TgBD shared 76.80% similarity to its close relative ark shell (Scapharca broughtoni). In addition, TgBD expression was observed in all tissues investigated under normal conditions and was significantly induced by injection of Vibrio parahaemolyticus. Furthermore, synthetic N-terminal peptide of TgBD exhibited strong antimicrobial activity against Gram-positive bacteria tested. Our results indicated that TgBD is a constitutive and inducible acute phase AMP, which provides a universal and prompt protection for T. granosa.

The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

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 of the estuarine oyster provides insights into climate impact and adaptive plasticity

Understanding the roles of genetic divergence and phenotypic plasticity in adaptation is central to evolutionary biology and important for assessing adaptive potential of species under climate change. Analysis of a chromosome-level assembly and resequencing of individuals across wide latitude distribution in the estuarine oyster (Crassostrea ariakensis) revealed unexpectedly low genomic diversity and population structures shaped by historical glaciation, geological events and oceanographic forces. Strong selection signals were detected in genes responding to temperature and salinity stress, especially of the expanded solute carrier families, highlighting the importance of gene expansion in environmental adaptation. Genes exhibiting high plasticity showed strong selection in upstream regulatory regions that modulate transcription, indicating selection favoring plasticity. Our findings suggest that genomic variation and population structure in marine bivalves are heavily influenced by climate history and physical forces, and gene expansion and selection may enhance phenotypic plasticity that is critical for the adaptation to rapidly changing environments.