Clam Genome Sequence Clarifies the Molecular Basis of Its Benthic Adaptation and Extraordinary Shell Color Diversity

Exploring the Heat-Responsive miRNAs and their Target Gene Regulation in Ruditapes philippinarum Under Acute Heat Stress

This study aimed to investigate the inherent molecular regulatory mechanisms of Ruditapes philippinarum in response to extremely high-temperature environments and to enhance the sustainable development of the R. philippinarum aquaculture industry. In this study, we established a differential expression profile of miRNA under acute heat stress and identified a total of 46 known miRNAs and 80 novel miRNAs, three of which were detected to be significantly differentially expressed. We analyzed the functions of target genes regulated by differentially expressed miRNAs (DEMs) of R. philippinarum. The findings of the KEGG enrichment analysis revealed that 29 enriched pathways in the group were subjected to acute heat stress. Notably, fatty acid metabolism, FoxO signaling pathway, TGF-β signaling pathway, and ubiquitin-mediated proteolysis were found to play significant roles in response to acute heat stress. We established a regulatory map of DEMs and their target genes in response to heat stress and constructed the miRNA-mRNA regulation network. This study provides valuable insights into the response of R. philippinarum to high temperature, helping to understand its underlying molecular regulatory mechanisms under high-temperature stress.

Physiological and transcriptomic analysis provides new insights into osmoregulation mechanism of Ruditapes philippinarum under low and high salinity stress

The Manila clam (Ruditapes philippinarum) is a commercially important marine bivalve, which inhabits the estuarine and mudflat areas. The osmoregulation is of great significance for molluscs adaptation to salinity fluctuations. In this study, we investigated the effects of low salinity (10 psu) and high salinity (40 psu) stress on survival and osmoregulation of the R. philippinarum. The results of physiological parameters showed that the ion (Na+, K+, Cl-) concentrations and Na+/K+-ATPase (NKA) activity of R. philippinarum decreased significantly under low salinity stress, but increased significantly under high salinity stress, indicating that there are differences in physiological adaptation of osmoregulation of R. philippinarum. In addition, we conducted the transcriptome analysis in the gills of R. philippinarum exposed to low (10 psu) and high (40 psu) salinity challenge for 48 h using RNA-seq technology. A total of 153 and 640 differentially expressed genes (DEGs) were identified in the low salinity (LS) group and high salinity (HS) group, respectively. The immune (IAP, TLR6, C1QL4, Ank3), ion transport (Slc34a2, SLC39A14), energy metabolism (PCK1, LDLRA, ACOX1) and DNA damage repair-related genes (Gadd45g, HSP70B2, GATA4) as well as FoxO, protein processing in endoplasmic reticulum and endocytosis pathways were involved in osmoregulation under low salinity stress of R. philippinarum. Conversely, the ion transport (SLC6A7, SLC6A9, SLC6A14, TRPM2), amino acid metabolism (GS, TauD, ABAT, ALDH4A1) and immune-related genes (MAP2K6, BIRC7A, CTSK, GVIN1), and amino acid metabolism pathways (beta-Alanine, Alanine, aspartate and glutamate, Glutathione) were involved in the process of osmoregulation under high salinity stress. The results obtained here revealed the difference of osmoregulation mechanism of R. philippinarum under low and high salinity stress through physiological and molecular levels. This study contributes to the assessment of salinity adaptation of bivalves in the context of climate change and provides useful information for marine resource conservation and aquaculture.

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

Genome-wide sequencing reveals geographical variations in the thermal adaptation of an aquaculture species with frequent seedling introductions

Climate change and human activity are essential factors affecting marine biodiversity and aquaculture, and understanding the impacts of human activities on the genetic structure to increasing high temperatures is crucial for sustainable aquaculture and marine biodiversity conservation. As a commercially important bivalve, the Manila clam Ruditapes philippinarum is widely distributed along the coast of China, and it has been frequently introduced from Fujian Province, China, to other regions for aquaculture. In this study, we collected four populations of Manila clams from different areas to evaluate their thermal tolerance by measuring cardiac performance and genetic variations using whole-genome resequencing. The upper thermal limits of the clams showed high variations within and among populations. Different populations displayed divergent genetic compositions, and the admixed population was partly derived from the Zhangzhou population in Fujian Province, implying a complex genomic landscape under the influence of local genetic sources and human introductions. Multiple single nucleotide polymorphisms (SNPs) were associated with the cardiac functional traits, and some of these SNPs can affect the codon usage and the structural stability of the resulting protein. This study shed light on the importance of establishing long-term ecological and genetic monitoring programs at the local level to enhance resilience to future climate change.

The lgi-miR-2d is Potentially Involved in Shell Melanin Synthesis by Targeting mitf in Manila Clam Ruditapes philippinarum

Shell color as an important economic trait is also the crucial target trait for breeding and production. MicroRNA (miRNA) is an endogenous small non-coding RNA that can post-transcriptionally regulate the expression of target genes, it plays important roles in many life activities and physiological processes, such as shell color, stress response, and disease traits. In this study, we investigated the function of lgi-miR-2d in shell melanin formation and the expression patterns of lgi-miR-2d and target gene Rpmitf in Manila clam Ruditapes philippinarum. We further explored and verified the relationship between Rpmitf and lgi-miR-2d and identified the expression level of shell color-related gene changes by RNAi and injecting the antagomir of lgi-miR-2d, respectively. Our results indicated that lgi-miR-2d antagomir affected the expression of its target gene Rpmitf. In addition, the dual-luciferase reporter assay was conducted to confirm the direct interaction between lgi-miR-2d and Rpmitf. The results showed that the expression levels of melanin-related genes such as Rpmitf and tyr were significantly decreased in the positive treatment group compared with the blank control group after the Rpmitf dsRNA injection, indicating Rpmitf plays a crucial role in the melanin synthesis pathway. Taken together, we speculated that lgi-miR-2d might be negatively modulating Rpmitf, which might regulate other shell color-related genes, thereby affecting melanin synthesis in R. philippinarum.

Genome-wide identification and expression analysis of Sox gene family in the Manila clam (Ruditapes philippinarum)

Sox transcription factors are vital in numerous fundamental biological processes. In this study, nine Sox gene family members were discovered in the Ruditapes philippinarum genome, classified into the SoxB1, SoxB2, SoxC, SoxD, SoxE, and SoxF groups, marking the first genome-wide identification of this gene family in R. philippinarum. Analyses of phylogeny, exon-intron structures, and domains bolster the support for their categorization and annotation. Furthermore, transcriptomic analyses across various developmental stages revealed that RpSox4, RpSox5, RpSox9, and RpSox11 were significantly expressed in the D-larval stage. Additionally, investigations into transcriptomes of clams with different shell colors indicated that most sox genes exhibited their highest expression levels in orange clams, followed by zebra, white zebra, and white clams, and the results of transcriptomes analysis in different tissues indicated that 8 Sox genes (except RpSox17) were highly expressed in the mantle tissue. Moreover, qPCR was used to detect the expression of Sox gene in R. philippinarum at different developmental periods, different shell colors and different tissues, and the results showed consistency with those of the transcriptomes. This study's findings lay the groundwork for additional exploration into the role of the Sox gene in melanin production in R. philippinarum shells.

cAMP-Mediated CREM-MITF-TYR Axis Regulates Melanin Synthesis in Pacific Oysters

Colorful shells in bivalves are mostly caused by the presence of biological pigments, among which melanin is a key component in the formation of shell colours. Cyclic adenosine monophosphate (cAMP) is an important messenger in the regulation of pigmentation in some species. However, the role of cAMP in bivalve melanogenesis has not yet been reported. In this study, we performed in vitro and in vivo experiments to determine the role of cAMP in regulating melanogenesis in Pacific oysters. Besides, the function of cAMP-responsive element modulator (CREM) and the interactions between CREM and melanogenic genes were investigated. Our results showed that a high level of cAMP promotes the expression of melanogenic genes in Pacific oysters. CREM controls the expression of the MITF gene under cAMP regulation. In addition, CREM can regulate melanogenic gene expression, tyrosine metabolism, and melanin synthesis. These results indicate that cAMP plays an important role in the regulation of melanogenesis in Pacific oysters. CREM is a key transcription factor in the oyster melanin synthesis pathway, which plays a crucial role in oyster melanin synthesis through a cAMP-mediated CREM-MITF-TYR axis.

Hypoxia in the Blue Mussel Mytilus chilensis Induces a Transcriptome Shift Associated with Endoplasmic Reticulum Stress, Metabolism, and Immune Response

The increase in hypoxia events, a result of climate change in coastal and fjord ecosystems, impacts the health and survival of mussels. These organisms deploy physiological and molecular responses as an adaptive mechanism to maintain cellular homeostasis under environmental stress. However, the specific effects of hypoxia on mussels of socioeconomic interest, such as Mytilus chilensis, are unknown. Using RNA-seq, we investigated the transcriptomic profiles of the gills, digestive gland, and adductor muscle of M. chilensis under hypoxia (10 days at 2 mg L-1) and reoxygenation (10 days at 6 mg L-1). There were 15,056 differentially expressed transcripts identified in gills, 11,864 in the digestive gland, and 9862 in the adductor muscle. The response varied among tissues, showing chromosomal changes in Chr1, Chr9, and Chr10 during hypoxia. Hypoxia regulated signaling genes in the Toll-like, mTOR, citrate cycle, and apoptosis pathways in gills, indicating metabolic and immunological alterations. These changes suggest that hypoxia induced a metabolic shift in mussels, reducing reliance on aerobic respiration and increasing reliance on anaerobic metabolism. Furthermore, hypoxia appeared to suppress the immune response, potentially increasing disease susceptibility, with negative implications for the mussel culture industry and natural bed populations. This study provides pivotal insights into metabolic and immunological adaptations to hypoxia in M. chilensis, offering candidate genes for adaptive traits.

Transcriptome analysis provides new insights into the immune response of Ruditapes philippinarum infected with Vibrio alginolyticus

Manila clam (Ruditapes philippinarum) is a bivalve species with commercial value, but it is easily infected by pathogenic microorganisms in aquaculture, which restricts the shellfish industry. Notably, the impact of Vibrio alginolyticus on clam culture is obvious. In this study, RNA-seq was performed to analyze clam hepatopancreas tissue in 48 h (challenge group, G48h) and 96 h (challenge group, G96h) after infection with V. alginolyticus and 0 h after injection of PBS (control group, C). The results showed that a total of 1670 differentially expressed genes were detected in the G48h vs C group, and 1427 differentially expressed genes were detected in the G96h vs C group. In addition, KEGG analysis showed that DEGs were significantly enriched in pathways such as Lysosome and Mitophagy. Moreover, 15 immune related DEGs were selected for qRT-PCR analysis to verify the accuracy of RNA-seq, and the results showed that the expression level of DEGs was consistent with that of RNA-seq. Therefore, the results obtained in this study provides a preliminary understanding of the immune defense of R. philippinarum and molecular insights for genetic breeding of V. alginolyticus resistance in Manila clam.

Transcriptomic and Physiological Analysis Reveal Melanin Synthesis-Related Genes and Pathways in Pacific Oysters (Crassostrea gigas)

Shell color is one of the shell traits of molluscs, which has been regarded as an economic trait in some bivalves. Pacific oysters (Crassostrea gigas) are important aquaculture shellfish worldwide. In the past decade, several shell color strains of C. gigas were developed through selective breeding, which provides valuable materials for research on the inheritance pattern and regulation mechanisms of shell color. The inheritance patterns of different shell colors in C. gigas have been identified in certain research; however, the regulation mechanism of oyster pigmentation and shell color formation remains unclear. In this study, we performed transcriptomic and physiological analyses using black and white shell oysters to investigate the molecular mechanism of melanin synthesis in C. gigas. Several pigmentation-related pathways, such as cytochrome P450, melanogenesis, tyrosine metabolism, and the cAMP signaling pathway were found. The majority of differentially expressed genes and some signaling molecules from these pathways exhibited a higher level in the black shell oysters than in the white, especially after L-tyrosine feeding, suggesting that those differences may cause a variation of tyrosine metabolism and melanin synthesis. In addition, the in vitro assay using primary cells from mantle tissue showed that L-tyrosine incubation increased cAMP level, gene and protein expression, and melanin content. This study reveals the difference in tyrosine metabolism and melanin synthesis in black and white shell oysters and provides evidence for the potential regulatory mechanism of shell color in oysters.

Expansion of the HSP70 gene family in Tegillarca granosa and expression profiles in response to zinc toxicity

Zinc (Zn) is an essential micronutrient in organisms and an abundant element in the Earth's crust. Trace amounts of Zn released from natural sources can enter aquatic ecosystems through weathering and erosion. Zn accumulates in organisms, and when its intracellular concentration exceeds a certain level, it can induce oxidative stress and trigger oxidative stress-mediated heat shock protein (HSP) modulation. HSP70 is the most evolutionarily conserved among the HSP families. Despite extensive research on HSP70 genes in bivalves, the HSP70 gene family of Tegillarca granosa is still poorly characterized. We identified 65 HSP70 genes belonging to 6 families in the T. granosa genome, with 50 HSPa12 and 11 HSPa B2 genes highly expanded. On chromosome 11, 39 HSP70 (60%) genes were identified, and the HSPa12A genes were highly duplicated. A total of 527 and 538 differentially expressed genes were identified in the gills and mantle based on Zn exposure, respectively. The Gene Ontology of cellular anatomical entities was significantly enriched with upregulated differentially expressed genes in the gills and mantle. Eight of the 11 HSPa B2 genes were upregulated in both tissues. Most of the genes identified in both tissues were involved in "protein homeostasis" and "inhibition of apoptosis," which are associated with the HSP70 family's resistance to extrinsic and intrinsic stress. Hence, this study identified that the HSP70 gene family plays a vital role in the adaptation of aquatic organisms to heavy metal (e.g., Zn) stress in contaminated environments by compiling the different physiological responses to preserve homeostasis.

Comparative Analysis of Bivalve and Sea Urchin Genetics and Development: Investigating the Dichotomy in Bilateria

This comprehensive review presents a comparative analysis of early embryogenesis in Protostomia and Deuterostomia, the first of which exhibit a mosaic pattern of development, where cells are fated deterministically, while Deuterostomia display a regulatory pattern of development, where the fate of cells is indeterminate. Despite these fundamental differences, there are common transcriptional mechanisms that underline their evolutionary linkages, particularly in the field of functional genomics. By elucidating both conserved and unique regulatory strategies, this review provides essential insights into the comparative embryology and developmental dynamics of these groups. The objective of this review is to clarify the shared and distinctive characteristics of transcriptional regulatory mechanisms. This will contribute to the extensive areas of functional genomics, evolutionary biology and developmental biology, and possibly lay the foundation for future research and discussion on this seminal topic.

Genome-Wide Identification and Characterization of MITF Genes in Ruditapes philippinarum and Their Involvement in the Immune Response to Vibrio anguillarum Infection

Studies have shown that the shellfish have innate immune system, which is a very important immune form of shellfish, and they rely on the innate immune system to resist diseases. As a transcription factor, Microphthalmia-associated transcription factor (MITF) plays a regulatory role in immune response and the shell color is also an important index for the breeding of excellent varieties of R. philippinarum. The research on immune response mechanism of RPMITFs can provide important reference data for the breeding of excellent clam varieties. In the genome of R. philippinarum, the RPMITF genes family of shell color-related gene family was selected as the target gene of this experiment. There are 12 RpMITF genes named RpMITF1, RpMITF2, RpMITF3, RpMITF4, RpMITF5, RpMITF6, RpMITF7, RpMITF8, RpMITF9, RpMITF10, RpMITF11, and RpMITF12. The open reading frame length is 639, 1233, 996, 1239, 675, 624, 816, 1365, 612, 1614, 1122, and 486 bp, encoding 212, 410, 331, 412, 224, 207, 271, 454, 203, 537, 373, and 161 aa, respectively. The predicted molecular weight range of amino acids is 18.85-62.61 kda, and the isoelectric point range is 5.26-9.44. Real-time quantitative PCR was used to detect the gene expression of RpMITF gene family in hepatopancreas tissues of two populations of Manila clam at 6 time points (0, 3, 6, 12, 24, and 48 h) after Vibrio anguillarum stress. The results show that RpMITF gene family was significantly expressed in hepatopancreas of two clam populations after V. anguillarum stress (P < 0.05).

Bioactive Peptides from Ruditapes philippinarum Attenuate Hypertension and Cardiorenal Damage in Deoxycorticosterone Acetate-Salt Hypertensive Rats

Hypertension is a common disease that affects human health and can lead to damage to the heart, kidneys, and other important organs. In this study, we investigated the regulatory effects of bioactive peptides derived from Ruditapes philippinarum (RPP) on hypertension and organ protection in deoxycorticosterone acetate (DOCA)-salt hypertensive rats. We found that RPPs exhibited significant blood pressure-lowering properties. Furthermore, the results showed that RPPs positively influenced vascular remodeling and effectively maintained a balanced water-sodium equilibrium. Meanwhile, RPPs demonstrated anti-inflammatory potential by reducing the serum levels of inflammatory cytokines (TNF-α, IL-2, and IL-6). Moreover, we observed the strong antioxidant activity of RPPs, which played a critical role in reducing oxidative stress and alleviating hypertension-induced damage to the aorta, heart, and kidneys. Additionally, our study explored the regulatory effects of RPPs on the gut microbiota, suggesting a possible correlation between their antihypertensive effects and the modulation of gut microbiota. Our previous studies have demonstrated that RPPs can significantly reduce blood pressure in SHR rats. This suggests that RPPs can significantly improve both essential hypertension and DOAC-salt-induced secondary hypertension and can ameliorate cardiorenal damage caused by hypertension. These findings further support the possibility of RPPs as an active ingredient in functional anti-hypertensive foods.

Transcriptome Responses to Different Environments in Intertidal Zones in the Peanut Worm Sipunculus nudus

The peanut worm (Sipunculus nudus) is an important intertidal species worldwide. Species living in the same aquaculture area might suffer different environmental impacts. To increase knowledge of the molecular mechanisms underlying the response to environmental fluctuations, we performed a transcriptome analysis of S. nudus from different intertidal zones using a combination of the SMRT platform and the Illumina sequencing platform. (1) A total of 105,259 unigenes were assembled, and 23,063 unigenes were perfectly annotated. The results of the PacBio Iso-Seq and IIIumina RNA-Seq enriched the genetic database of S. nudus. (2) A total of 830 DEGs were detected in S. nudus from the different groups. In particular, 33 DEGs had differential expression in the top nine KEGG pathways related to pathogens, protein synthesis, and cellular immune response and signaling. The results indicate that S. nudus from different zones experience different environmental stresses. (3) Several DEGs (HSPA1, NFKBIA, eEF1A, etc.) in pathways related to pathogens (influenza A, legionellosis, measles, and toxoplasmosis) had higher expression in groups M and L. HSPA1 was clearly enriched in most of the pathways, followed by NFKBIA. The results show that the peanut worms from the M and L tidal flats might have suffered more severe environmental conditions. (4) Some DEGs (MKP, MRAS, and HSPB1) were upregulated in peanut worms from the H tidal flat, and these DEGs were mainly involved in the MAPK signaling pathway. These results indicate that the MAPK pathway may play a vital role in the immune response of the peanut worm to the effects of different intertidal flats. This study provides a valuable starting point for further studies to elucidate the molecular basis of the response to different environmental stresses in S. nudus.

Genome-wide identification and expression analyses of Toll-like receptors provide new insights on adaptation to intertidal benthic environments in Urechis unicinctus (Annelida, Echiura)

Toll-like receptors (TLR) are an important class of pattern recognition receptors involved in innate immunity that recognize pathogen-associated and damage-associated molecular patterns. Although the role of TLRs in immunity has been extensively studied, a systematic investigation of their function in environmental adaptation is still in its infancy. In this study, a genome-wide search was conducted to systematically investigate TLR family members of Urechis unicinctus, a typical benthic organism in intertidal mudflats. A total of 28 TLR genes were identified in the U. unicinctus genome, and their fundamental physiological and biochemical properties were characterized. Gene copy number analysis among species in different habitats indicated that TLR family gene expansion may be probably related with benthic environmental adaptation. To further investigate the expression patterns of TLR members under environmental stress, transcriptome data was analyzed from different developmental stages and the hindgut under sulfide stress. Transcriptome analysis of different developmental stages showed that most TLR genes were highly expressed during a key period of benthic environment adaptation (worm-shaped larva). Transcriptome analysis of the hindgut under sulfide stress showed that the expression of 12 TLR members was significantly induced under sulfide stress. These results indicate that the regulation of TLR gene expression may be probably involved in the adaptation of U. unicinctus to the benthic intertidal zone environment. Taken together, this study may lay the foundation for future functional analysis of the specific role of TLRs in host immune responses against sulfide exposure and benthic environmental stress in annelid.

Physiological and molecular responses to hypoxia stress in Manila clam Ruditapes philippinarum

Hypoxia has become one of the major environmental problems in the aquaculture industry. As one of the most commercially important bivalves, Manila clam Ruditapes philippinarum may be suffering substantial mortality attributable to hypoxia. The physiological and molecular responses to hypoxia stress in Manila clam were evaluated at two levels of low dissolved oxygen: 0.5 mg/L (DO 0.5 mg/L) and 2.0 mg/L (DO 2.0 mg/L). With the prolongation of hypoxia stress, the mortality rate was 100% at 156 h under DO 0.5 mg/L. In contrast, 50% of clams survived after 240 h of stress at DO 2.0 mg/L. After the hypoxia stress, some severe structural damages were observed in gill, axe foot, hepatopancreas tissues, such as cell rupture and mitochondrial vacuolization. For the hypoxia-stressed clams, the significant rise and decline of enzyme activity (LDH and T-AOC) was observed in gills, in contrast to the reduction of glycogen content. Furthermore, the expression levels of genes related to energy metabolism (SDH, PK, Na⁺/K⁺-ATPase, NF-κB and HIF-1α) was significantly affected by the hypoxia stress. It is therefore suggested that the short-term survival of clams under hypoxia may be dependent on stress protection by antioxidants, energy allocation, and tissue energy reserves (such as glycogen stores). Despite this, the prolongation of hypoxia stress at DO 2.0 mg/L may cause the irreversible damages of cellular structures in clam tissues, eventually leading to the death of clams. We therefore support the hypothesis that the extent of hypoxia impacts on marine bivalves may be underestimated in the coastal areas.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves

BACKGROUND

Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution.

RESULTS

We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta.

CONCLUSIONS

Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.

Revealing the potential regulatory relationship between HSP70, HSP90 and HSF genes under temperature stress

Heat Shock Protein (HSPs) gene family members play fundamental roles in different environmental stress tolerances, protect the structure and function of cells, and perform a significant task in cellular homeostasis. In this study, we conducted a genome-wide identification, evolutionary relationship analysis and gene expression analysis of the HSP70, HSP90, and HSF gene families in Ruditapes philippinarum. We identified 83 RpHSP70, 6 RpHSP90, and 3 RpHSF genes in R. philippinarum. The structural characteristics, chromosomal localization, and the gene structure map were constructed to reveal the characteristics of protein structures. Furthermore, the expression profiling of transcriptome data showed the expression pattern of HSP70, HSP90 and HSF genes in Manila clam from different populations, and under high and low temperature stress. In addition, we performed protein-protein interaction network analysis between HSP70, HSP90, and HSF gene family which enabled us to recognize the regulatory relationship between the two HSP gene families and the HSF gene family. Furthermore, the predicted sub-cellular location revealed a diversified subcellular distribution of HSP70, HSP90, and HSF proteins, which may be directly or indirectly associated with functional diversification under heat stress condition.

Chromosome-level genome assembly of the European flat oyster (Ostrea edulis) provides insights into its evolution and adaptation

The European flat oyster (Ostrea edulis) is an endangered and economically important marine bivalve species that plays a critical role in the coastal ecosystem. Here, we report a high-quality chromosome-level genome assembly of O. edulis, generated using PacBio HiFi-CCS long reads and annotated with Nanopore full-length transcriptome. The O. edulis genome covers 946.06 Mb (scaffold N50 94.82 Mb) containing 34,495 protein-coding genes and a high proportion of repeat sequences (58.49 %). The reconstructed demographic histories show that O. edulis population might be shaped by breeding habit (embryo brooding) and historical climatic change. Comparative genomic analysis indicates that transposable elements may drive lineage-specific evolution in oysters. Notably, the O. edulis genome has a Hox gene cluster rearrangement that has never been reported in bivalves, making this species valuable for evolutionary studies of molluscan diversification. Moreover, genome expansion of O. edulis is probably central to its adaptation to filter-feeding and sessile lifestyles, as well as embryo brooding and pathogen resistance, in coastal ecosystems. This chromosome-level genome assembly provides new insights into the genome feature of oysters, and presents an important resource for genetic research, evolutionary studies, and biological conservation of O. edulis.

Transient receptor potential (TRP) channels in the Manila clam (Ruditapes philippinarum): Characterization and expression patterns of the TRP gene family under heat stress in Manila clams based on genome-wide identification

In this study, we identified a total of 40 transient receptor potential genes (RpTRP) in Manila clam by genome-wide identification and classified them into four categories (TRPV, TRPA, TRPM, TRPC) based on gene structure and subfamily relationships. The protein length of RpTRP genes ranges from 281 amino acids to 1601 amino acids. Molecular weight and theoretical PI values range from 182.82 kDa to 32.43 kDa, respectively, with PI values between 5.17 and 9.25. By comparing the expression profiles of TRP genes during heat stress in Manila clams at different latitudes, we found that most genes in the TRP gene family were up-regulated in expression during heat challenge. Therefore, we determined that TRP genes have an important role in the heat stress of Manila clams. This work provides a basis for further studies on the molecular mechanisms of TRP-mediated heat tolerance in Manila clam and for explaining differences in heat tolerance in Manila clam at different latitudes through key differential TRP genes at the molecular level.

Multi-tissue RNA-Seq Analysis and Long-read-based Genome Assembly Reveal Complex Sex-specific Gene Regulation and Molecular Evolution in the Manila Clam

The molecular factors and gene regulation involved in sex determination and gonad differentiation in bivalve molluscs are unknown. It has been suggested that doubly uniparental inheritance (DUI) of mitochondria may be involved in these processes in species such as the ubiquitous and commercially relevant Manila clam, Ruditapes philippinarum. We present the first long-read-based de novo genome assembly of a Manila clam, and a RNA-Seq multi-tissue analysis of 15 females and 15 males. The highly contiguous genome assembly was used as reference to investigate gene expression, alternative splicing, sequence evolution, tissue-specific co-expression networks, and sexual contrasting SNPs. Differential expression (DE) and differential splicing (DS) analyses revealed sex-specific transcriptional regulation in gonads, but not in somatic tissues. Co-expression networks revealed complex gene regulation in gonads, and genes in gonad-associated modules showed high tissue specificity. However, male gonad-associated modules showed contrasting patterns of sequence evolution and tissue specificity. One gene set was related to the structural organization of male gametes and presented slow sequence evolution but high pleiotropy, whereas another gene set was enriched in reproduction-related processes and characterized by fast sequence evolution and tissue specificity. Sexual contrasting SNPs were found in genes overrepresented in mitochondrial-related functions, providing new candidates for investigating the relationship between mitochondria and sex in DUI species. Together, these results increase our understanding of the role of DE, DS, and sequence evolution of sex-specific genes in an understudied taxon. We also provide resourceful genomic data for studies regarding sex diagnosis and breeding in bivalves.

Examination of the potential roles of insulin-like peptide receptor in regulating the growth of Manila clam Ruditapes philippinarum

The role of insulin/insulin-like growth factor (IGF) signaling pathway in the growth regulation of marine invertebrates has not been fully studied. In this study, the economically important species Ruditapes philippinarum was sacrificed to clarify the role of IGF system in the growth regulation of R. philippinarum by real-time quantitative PCR. Systematic bioinformatics analysis can identify the major genes of IGF signaling pathway and insulin-like peptide receptor (ILPR) - mediated signaling pathway in R. philippinarum. The expression levels of IGF and its downstream signaling pathway genes in larger clams were significantly higher than those in small clams, indicating that they were involved in the growth regulation of R. philippinarum. These results suggest that IGF signaling pathway and ILPR mediated signaling pathway to regulate the growth of R. philippinarum.

A high-quality, haplotype-phased genome reconstruction reveals unexpected haplotype diversity in a pearl oyster

Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.

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.

Animal board invited review: Widespread adoption of genetic technologies is key to sustainable expansion of global aquaculture

•

The extent of application of genetic technologies to aquaculture production varies widely by species and geography.

•

Achieving a more universal application of seed derived from scientifically based breeding programmes is an important goal in order to meet increasing global demands for seafood production.

•

This article reviews the status of genetic technologies across the world’s top 10 highly produced species.

•

Opportunities and barriers to achieving broad-scale uptake of genetic technologies in global aquaculture are discussed.

•

A future outlook for potential disruptive genetic technologies and how they might affect global aquaculture production is given.

Aquaculture production comprises a diverse range of species, geographies, and farming systems. The application of genetics and breeding technologies towards improved production is highly variable, ranging from the use of wild-sourced seed through to advanced family breeding programmes augmented by genomic techniques. This technical variation exists across some of the most highly produced species globally, with several of the top ten global species by volume generally lacking well-managed breeding programmes. Given the well-documented incremental and cumulative benefits of genetic improvement on production, this is a major missed opportunity. This short review focusses on (i) the status of application of selective breeding in the world’s most produced aquaculture species, (ii) the range of genetic technologies available and the opportunities they present, and (iii) a future outlook towards realising the potential contribution of genetic technologies to aquaculture sustainability and global food security.

Comparative transcriptomic analysis revealed dynamic changes of distinct classes of genes during development of the Manila clam (Ruditapes philippinarum)

Background

The Manila clam Ruditapesphilippinarum is one of the most economically important marine shellfish. However, the molecular mechanisms of early development in Manila clams are largely unknown. In this study, we collected samples from 13 stages of early development in Manila clam and compared the mRNA expression pattern between samples by RNA-seq techniques.

Results

We applied RNA-seq technology to 13 embryonic and larval stages of the Manila clam to identify critical genes and pathways involved in their development and biological characteristics. Important genes associated with different morphologies during the early fertilized egg, cell division, cell differentiation, hatching, and metamorphosis stages were identified. We detected the highest number of differentially expressed genes in the comparison of the pediveliger and single pipe juvenile stages, which is a time when biological characteristics greatly change during metamorphosis. Gene Ontology (GO) enrichment analysis showed that expression levels of microtubule protein-related molecules and Rho genes were upregulated and that GO terms such as ribosome, translation, and organelle were enriched in the early development stages of the Manila clam. Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that the foxo, wnt, and transforming growth factor-beta pathways were significantly enriched during early development. These results provide insights into the molecular mechanisms at work during different periods of early development of Manila clams.

Conclusion

These transcriptomic data provide clues to the molecular mechanisms underlying the development of Manila clam larvae. These results will help to improve Manila clam reproduction and development.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08813-0.

Genome-wide identification and analysis of HECT E3 ubiquitin ligase gene family in Ruditapes philippinarum and their involvement in the response to heat stress and Vibrio anguillarum infection

E3 ubiquitin ligase (E3s) plays an important role in ubiquitin proteasome pathway, proteins containing homologous E6-AP carboxyl terminus (HECT) domains. However, the role of HECT E3 ubiquitin ligase in mollusk was rarely explored. In this study, we performed a genome-wide analysis of the HECT domain-containing gene in Ruditapes philippinarum to identify and predict the structural and functional characterization of HECT genes in response to abiotic and biotic stress. A total of sixteen members of HECT gene family were identified and analyzed for the gene structure, phylogenetic relation, three-dimensional structure, protein interaction network, and expression patterns. Experimental results demonstrated that Rph.HUWE1, Rph.HECTD1, Rph.Ubr5 were significantly up-regulated in response to heat stress and bacterial challenge. Taken together, our data provide insights into the potential function of HECT E3 ligase in heat stress and Vibrio anguillarum infection.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.

Apextrin from Ruditapes philippinarum functions as pattern recognition receptor and modulates NF-κB pathway

Apextrin belongs to ApeC-containing proteins (ACPs) and features a signal-peptide, an N-terminal membrane attack complex component/perforin (MACPF) domain, and a C-terminal ApeC domain. Recently, apextrin-like proteins were identified as pattern recognition receptor (PRR), which recognize the bacterial cell wall component and participate in innate immunity. Here, an apextrin (Rpape) was identified and characterized in Ruditapes philippinarum. Our results showed that Rpape mRNA was significantly induced under bacterial challenges. The Rpape recombinant protein exhibited a significant inhibitory effect on gram-positive bacteria (Bacillus subtilis and Staphylococcus aureus) and bound with Vibrio anguillarum, S. aureus and B. subtilis. We found Rpape protein positively activated the NF-κB signaling cascade and increased the activity of Nitric oxide (NO). This study revealed the immunity role of apextrin in R. philippinarum and provided a reference for further study on the role of apextrin in bivalves.

Domain-Dependent Evolution Explains Functional Homology of Protostome and Deuterostome Complement C3-Like Proteins

Complement proteins emerged early in evolution but outside the vertebrate clade they are poorly characterized. An evolutionary model of C3 family members revealed that in contrast to vertebrates the evolutionary trajectory of C3-like genes in cnidarian, protostomes and invertebrate deuterostomes was highly divergent due to independent lineage and species-specific duplications. The deduced C3-like and vertebrate C3, C4 and C5 proteins had low sequence conservation, but extraordinarily high structural conservation and 2-chain and 3-chain protein isoforms repeatedly emerged. Functional characterization of three C3-like isoforms in a bivalve representative revealed that in common with vertebrates complement proteins they were cleaved into two subunits, b and a, and the latter regulated inflammation-related genes, chemotaxis and phagocytosis. Changes within the thioester bond cleavage sites and the a-subunit protein (ANATO domain) explained the functional differentiation of bivalve C3-like. The emergence of domain-related functions early during evolution explains the overlapping functions of bivalve C3-like and vertebrate C3, C4 and C5, despite low sequence conservation and indicates that evolutionary pressure acted to conserve protein domain organization rather than the primary sequence.

Comparative analysis of the Mercenaria mercenaria genome provides insights into the diversity of transposable elements and immune molecules in bivalve mollusks

BACKGROUND

The hard clam Mercenaria mercenaria is a major marine resource along the Atlantic coasts of North America and has been introduced to other continents for resource restoration or aquaculture activities. Significant mortality events have been reported in the species throughout its native range as a result of diseases (microbial infections, leukemia) and acute environmental stress. In this context, the characterization of the hard clam genome can provide highly needed resources to enable basic (e.g., oncogenesis and cancer transmission, adaptation biology) and applied (clam stock enhancement, genomic selection) sciences.

RESULTS

Using a combination of long and short-read sequencing technologies, a 1.86 Gb chromosome-level assembly of the clam genome was generated. The assembly was scaffolded into 19 chromosomes, with an N50 of 83 Mb. Genome annotation yielded 34,728 predicted protein-coding genes, markedly more than the few other members of the Venerida sequenced so far, with coding regions representing only 2% of the assembly. Indeed, more than half of the genome is composed of repeated elements, including transposable elements. Major chromosome rearrangements were detected between this assembly and another recent assembly derived from a genetically segregated clam stock. Comparative analysis of the clam genome allowed the identification of a marked diversification in immune-related proteins, particularly extensive tandem duplications and expansions in tumor necrosis factors (TNFs) and C1q domain-containing proteins, some of which were previously shown to play a role in clam interactions with infectious microbes. The study also generated a comparative repertoire highlighting the diversity and, in some instances, the specificity of LTR-retrotransposons elements, particularly Steamer elements in bivalves.

CONCLUSIONS

The diversity of immune molecules in M. mercenaria may allow this species to cope with varying and complex microbial and environmental landscapes. The repertoire of transposable elements identified in this study, particularly Steamer elements, should be a prime target for the investigation of cancer cell development and transmission among bivalve mollusks.

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.

A chromosome-level genome of Portunus trituberculatus provides insights into its evolution, salinity adaptation and sex determination

Portunus trituberculatus (Crustacea: Decapoda: Brachyura), commonly known as the swimming crab, is of major ecological importance, as well as being important to the fisheries industry. P. trituberculatus is also an important farmed species in China due to its rapid growth rate and high economic value. Here, we report the genome sequence of the swimming crab, which was assembled at the chromosome scale, covering ~1.2 Gb, with 79.99% of the scaffold sequences assembled into 53 chromosomes. The contig and scaffold N50 values were 108.7 kb and 15.6 Mb, respectively, with 19,981 protein-coding genes. Based on comparative genomic analyses of crabs and shrimps, the C2H2 zinc finger protein family was found to be the only gene family expanded in crab genomes, suggesting it was closely related to the evolution of crabs. The combination of transcriptome and bulked segregant analysis provided insights into the genetic basis of salinity adaptation and rapid growth in P. trituberculatus. In addition, the specific region of the Y chromosome was located for the first time in the genome of P. trituberculatus, and three genes were preliminarily identified as candidate genes for sex determination in this region. Decoding the swimming crab genome not only provides a valuable genomic resource for further biological and evolutionary studies, but is also useful for molecular breeding of swimming crabs.

MiRNA-mRNA Integration Analysis Reveals the Regulatory Roles of MiRNAs in Shell Pigmentation of the Manila clam (Ruditapes philippinarum)

The shell color of the Manila clam (Ruditapes philippinarum) is an economically important trait. We used high-throughput sequencing and transcriptome analysis to study the molecular mechanisms that underlie shell color formation and regulation in this species. We constructed small RNA libraries from mantle tissues from four shell color strains of Manila clam, subjected them to high-throughput sequencing. Notably, the results suggested that a number of pigment-associated genes including Mitf, HERC2, were negatively regulated by nvi-miR-2a, tgu-miR-133-3p, respectively. They might be involved in melanin formation via the activation of the melanogenesis pathway. And aae-miR-71-5p and dme-miR-7-5p linked to shell formation-related genes such as Calmodulin and IMSP3 were considered to participate in the calcium signaling pathway. We then used quantitative PCR to verify the candidate miRNAs and target genes in different shell color groups. Our results indicated that miR-7, miR-71, and miR-133 may regulate target mRNAs to participate in shell color pigmentation. These results provide the foundation to further characterize miRNA effects on the regulation of shell color and have significant implications for the breeding of new varieties of clams.

Phalloidin fluorescence and confocal microscopy reveal the musculature development of clam Ruditapes philippinarum

Knowledge of early development in bivalves is of great importance to understand the function of animal organ systems and the evolution of phenotypic diversity. Manila clam (Ruditapes philippinarum) is an economically important bivalve living in marine intertidal zones. To determine the pattern of muscle development in the clams, we investigate the characteristics of musculature development utilizing phalloidin staining and confocal microscopy. Myofilaments first appear at the early trochophore stage, and gradually become orderly arranged during the transition from trochophore to veliger. For veliger, larval muscle system is mainly composed of dorsal velum retractors, medio-dorsal velum retractors, ventral velum retractors, ventral larval retractors and anterior and posterior adductor muscles. After metamorphosis, the muscle system of late veliger has been reconstructed, showing the irreversible shrink of velum retractor muscles, the rapid growth of wedge-shaped foot and mantle margins. One of the most striking changes in settled spats is the development of sophisticated architecture of foot musculature, which consists of transverse pedal muscles, anterior foot retractor and posterior foot retractor. These findings will not only provide the basis to understand the dynamic pattern of myogenesis in the burrowing bivalves, but also provide valuable information for comparative analysis of muscle development among bivalves.

Dissecting the chromosome-level genome of the Asian Clam (Corbicula fluminea)

The Asian Clam (Corbicula fluminea) is a valuable commercial and medicinal bivalve, which is widely distributed in East and Southeast Asia. As a natural nutrient source, the clam is rich in protein, amino acids, and microelements. The genome of C. fluminea has not yet been characterized; therefore, genome-assisted breeding and improvements cannot yet be implemented. In this work, we present a de novo chromosome-scale genome assembly of C. fluminea using PacBio and Hi-C sequencing technologies. The assembled genome comprised 4728 contigs, with a contig N50 of 521.06 Kb, and 1,215 scaffolds with a scaffold N50 of 70.62 Mb. More than 1.51 Gb (99.17%) of genomic sequences were anchored to 18 chromosomes, of which 1.40 Gb (92.81%) of genomic sequences were ordered and oriented. The genome contains 38,841 coding genes, 32,591 (83.91%) of which were annotated in at least one functional database. Compared with related species, C. fluminea had 851 expanded gene families and 191 contracted gene families. The phylogenetic tree showed that C. fluminea diverged from Ruditapes philippinarum, ~ 228.89 million years ago (Mya), and the genomes of C. fluminea and R. philippinarum shared 244 syntenic blocks. Additionally, we identified 2 MITF members and 99 NLRP members in C. fluminea genome. The high-quality and chromosomal Asian Clam genome will be a valuable resource for a range of development and breeding studies of C. fluminea in future research.

Evolutionary History of DNA Methylation Related Genes in Bivalvia: New Insights From Mytilus galloprovincialis

DNA methylation is an essential epigenetic mechanism influencing gene expression in all organisms. In metazoans, the pattern of DNA methylation changes during embryogenesis and adult life. Consequently, differentiated cells develop a stable and unique DNA methylation pattern that finely regulates mRNA transcription during development and determines tissue-specific gene expression. Currently, DNA methylation remains poorly investigated in mollusks and completely unexplored in Mytilus galloprovincialis. To shed light on this process in this ecologically and economically important bivalve, we screened its genome, detecting sequences homologous to DNA methyltransferases (DNMTs), methyl-CpG-binding domain (MBD) proteins and Ten-eleven translocation methylcytosine dioxygenase (TET) previously described in other organisms. We characterized the gene architecture and protein domains of the mussel sequences and studied their phylogenetic relationships with the ortholog sequences from other bivalve species. We then comparatively investigated their expression levels across different adult tissues in mussel and other bivalves, using previously published transcriptome datasets. This study provides the first insights on DNA methylation regulators in M. galloprovincialis, which may provide fundamental information to better understand the complex role played by this mechanism in regulating genome activity in bivalves.

Chromosome-level genome assembly of Scapharca kagoshimensis reveals the expanded molecular basis of heme biosynthesis in ark shells

Ark shells are commercially important clam species that inhabit in muddy sediments of shallow coasts in East Asia. For a long time, the lack of genome resources has hindered scientific research of ark shells. Here, we report a high-quality chromosome-level genome assembly of Scapharca kagoshimensis, with an aim to unravel the molecular basis of heme biosynthesis, and develop genomic resources for genetic breeding and population genetics in ark shells. Nineteen scaffolds corresponding to 19 chromosomes were constructed from 938 contigs (contig N50 = 2.01 Mb) to produce a final high-quality assembly with a total length of 1.11 Gb and scaffold N50 around 60.64 Mb. The genome assembly represents 93.4% completeness via matching 303 eukaryota core conserved genes. A total of 24,908 protein-coding genes were predicted and 24,551 genes (98.56%) of which were functionally annotated. The enrichment analyses suggested that genes in heme biosynthesis pathways were expanded and positive selection of the haemoglobin genes was also found in the genome of S. kagoshimensis, which gives important insights into the molecular mechanisms and evolution of the heme biosynthesis in mollusca. The valuable genome assembly of S. kagoshimensis would provide a solid foundation for investigating the molecular mechanisms that underlie the diverse biological functions and evolutionary adaptations of S. kagoshimensis.

Integrated transcriptome and metabolome analysis reveals molecular responses of the clams to acute hypoxia

Mudflat shellfish have evolved well-adapted strategies for coping with dynamic environmental fluxes and stressful conditions, including oxygen availability. The Manila clams Ruditapes philippinarum are worldwide cultured shellfish in marine intertidal zone, which usually encounter great risk of acute hypoxia exposure in coastal habitats. To reveal the effects of acute hypoxia on metabolic changes of the clams, we performed the integrated analysis of transcriptomics and metabolomics to investigate the global changes of genes and metabolites during acute hypoxia stress at the whole-organism level. The comparative transcriptome analysis reveals that the clams show the remarkable depression in a variety of biological performance, such as metabolic rates, neuronal activity, biomineralization activity, and cell proliferation and differentiation at the hypoxic condition. The metabolomic analysis reveals that amino acid metabolism plays a critical role in the metabolic changes of the clams in response to acute hypoxia. A variety of free amino acids may not only be served as the potential osmolytes for osmotic regulation, but also may contribute to energy production during the acute hypoxia exposure. The metabolite analysis also reveals several important biomarkers for metabolic changes, and provides new insights into how clams deal with acute hypoxia. These findings suggest that clams may get through acute hypoxia stress by the adaptive metabolic strategy to survive short-period of acute hypoxia which is likely to occur in their typical habitat. The present findings will not only shed lights on the molecular and metabolic mechanisms of adaptive strategies under stressful conditions, but also provide the signaling metabolites to assess the physiological states of clams in aquaculture.

The Crown Pearl: a draft genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Since historical times, the inherent human fascination with pearls turned the freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758) into a highly valuable cultural and economic resource. Although pearl harvesting in M. margaritifera is nowadays residual, other human threats have aggravated the species conservation status, especially in Europe. This mussel presents a myriad of rare biological features, e.g. high longevity coupled with low senescence and Doubly Uniparental Inheritance of mitochondrial DNA, for which the underlying molecular mechanisms are poorly known. Here, the first draft genome assembly of M. margaritifera was produced using a combination of Illumina Paired-end and Mate-pair approaches. The genome assembly was 2.4 Gb long, possessing 105,185 scaffolds and a scaffold N50 length of 288,726 bp. The ab initio gene prediction allowed the identification of 35,119 protein-coding genes. This genome represents an essential resource for studying this species’ unique biological and evolutionary features and ultimately will help to develop new tools to promote its conservation.

Integrated transcriptomics and metabolomics analyses reveal benzo[a]pyrene enhances the toxicity of mercury to the Manila clam, Ruditapes philippinarum

Mercury (Hg²⁺) and benzo[a]pyrene (BaP) are ubiquitous and persistent pollutants with multiple toxicities in bivalve molluscs. Here, the toxicological responses in the gills of Manila clams, Ruditapes philippinarum, to Hg²⁺ (10 μg L^-1), BaP (3 μg L^-1), and their mixture were analysed using transcriptomics and metabolomics approaches. Comparisons of the transcriptomes and metabolomes of Hg²⁺-and/or BaP-treated clams with control animals revealed the involvement of the detoxification metabolism, immune defence, energy-related pathways, and osmotic regulation in the stress response of R. philippinarum. Exposure to Hg²⁺ alone primarily enhanced the detoxification and energy metabolic pathways by significantly increasing the expression of genes associated with heat-shock proteins and oxidative phosphorylation. However, co-exposure to Hg²⁺ and BaP caused greater immunotoxicity and disrupted detoxification metabolism, the TCA cycle, glycolysis, and ATP generation. The expression levels of cytochrome P450 1A1 (CYP1A1), multidrug resistance-associated protein 1 (MRP1), and myosin (MYO), and the activity of electron transport system (ETS) in gills were detected, supporting the underlying toxic mechanisms of Hg²⁺ and BaP. We suggest that the presence of BaP enhances the toxicity of Hg²⁺ by 1) hampering the detoxification of Hg²⁺, 2) increasing the immunotoxicity of Hg²⁺, and 3) constraining energy availability for clams.

Identification of shell-color-related microRNAs in the Manila clam Ruditapes philippinarum using high-throughput sequencing of small RNA transcriptomes

Shell-color polymorphism is a common phenomenon in several mollusk species and has been associated with thermal capacity, developmental stability, shell strength, and immunity. Shell-color polymorphism has been related to the differential expression of genes in several signal transduction pathways; however, the functions of micro-RNAs (miRNAs) in shell-color formation remain unclear. In the present study, we compared high-quality, small-RNA transcriptomes in three strains of the Manila clam Ruditapes philippinarum with specific shell-color patterns, artificially selected for six generations. Totals of 114 known and 208 novel miRNAs were identified by high-throughput sequencing, of which nine known and one novel miRNA were verified by stem-loop quantitative real time-polymerase chain reaction. Predicted miRNA targets were subjected to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. miR-137 and miR-216b and the Hedgehog signaling pathway and Wnt signaling pathway were identified as being potentially involved in pigment formation and regulation in R. philippinarum. These results may help to clarify the role of miRNAs in shell coloration and shed light on the mechanisms regulating color formation in bivalve shells.

Low impact of different SNP panels from two building-loci pipelines on RAD-Seq population genomic metrics: case study on five diverse aquatic species

Background

The irruption of Next-generation sequencing (NGS) and restriction site-associated DNA sequencing (RAD-seq) in the last decade has led to the identification of thousands of molecular markers and their genotyping for refined genomic screening. This approach has been especially useful for non-model organisms with limited genomic resources. Many building-loci pipelines have been developed to obtain robust single nucleotide polymorphism (SNPs) genotyping datasets using a de novo RAD-seq approach, i.e. without reference genomes. Here, the performances of two building-loci pipelines, STACKS 2 and Meyer’s 2b-RAD v2.1 pipeline, were compared using a diverse set of aquatic species representing different genomic and/or population structure scenarios. Two bivalve species (Manila clam and common edible cockle) and three fish species (brown trout, silver catfish and small-spotted catshark) were studied. Four SNP panels were evaluated in each species to test both different building-loci pipelines and criteria for SNP selection. Furthermore, for Manila clam and brown trout, a reference genome approach was used as control.

Results

Despite different outcomes were observed between pipelines and species with the diverse SNP calling and filtering steps tested, no remarkable differences were found on genetic diversity and differentiation within species with the SNP panels obtained with a de novo approach. The main differences were found in brown trout between the de novo and reference genome approaches. Genotyped vs missing data mismatches were the main genotyping difference detected between the two building-loci pipelines or between the de novo and reference genome comparisons.

Conclusions

Tested building-loci pipelines for selection of SNP panels seem to have low influence on population genetics inference across the diverse case-study scenarios here studied. However, preliminary trials with different bioinformatic pipelines are suggested to evaluate their influence on population parameters according with the specific goals of each study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07465-w.

A High-Quality Reference Genome for a Parasitic Bivalve with Doubly Uniparental Inheritance (Bivalvia: Unionida)

From a genomics perspective, bivalves (Mollusca: Bivalvia) have been poorly explored with the exception for those of high economic value. The bivalve order Unionida, or freshwater mussels, has been of interest in recent genomic studies due to their unique mitochondrial biology and peculiar life cycle. However, genomic studies have been hindered by the lack of a high-quality reference genome. Here, I present a genome assembly of Potamilus streckersoni using Pacific Bioscience single-molecule real-time long reads and 10X Genomics-linked read sequencing. Further, I use RNA sequencing from multiple tissue types and life stages to annotate the reference genome. The final assembly was far superior to any previously published freshwater mussel genome and was represented by 2,368 scaffolds (2,472 contigs) and 1,776,755,624 bp, with a scaffold N50 of 2,051,244 bp. A high proportion of the assembly was comprised of repetitive elements (51.03%), aligning with genomic characteristics of other bivalves. The functional annotation returned 52,407 gene models (41,065 protein, 11,342 tRNAs), which was concordant with the estimated number of genes in other freshwater mussel species. This genetic resource, along with future studies developing high-quality genome assemblies and annotations, will be integral toward unraveling the genomic bases of ecologically and evolutionarily important traits in this hyper-diverse group.

Genome-wide investigation and expression analysis of MACPF gene family reveals its immune role in response to bacterial challenge of Manila clam

In this study, 18 MACPF genes (RpMACPF) were identified and classed into three types (Macrophage-expressed gene 1, Apextrin, and MACPF domain contain protein) based on gene structure and phylogenetic relationship in R. philippinarum. The length of RpMACPF proteins varied from 287 to 785 amino acids. The molecular weights and Theoretical PI values ranged from 3.2 kDa to 8.7 kDa and 4.7 to 8.6, respectively. RNA-seq data analysis revealed that 14 of 18 RpMACPF genes were highly expressed at the pediveliger larvae stage indicate RpMACPF might contribute to the early development and metamorphosis processes of the R. philippinarum. Besides, we found RpMACPF genes were significantly regulated by pathogen-associated molecular patterns (PAMPs) and Vibrio parahemolyticus, which indicates RpMACPF genes may play significant roles in response to invading pathogens. The results obtained in this work will provide valuable insight into the immune function of MACPF gene in R. philippinarum.

Molecular mechanisms of wound healing and regeneration of siphon in the Manila clam Ruditapes philippinarum revealed by transcriptomic analysis

Ruditapes philippinarum is an economically important marine shellfish aquaculture species, and it has the ability to regenerate its siphons. To gain a greater understanding of the molecular mechanisms at work during siphon regeneration and to provide evidence for morphological regeneration, we examined transcriptome responses of siphon tissue of R. philippinarum during regeneration and observed regenerative siphons under the stereomicroscope. The overall process of siphon regeneration was dissected based on the morphological changes of siphon and the identification of up-regulated key differentially expressed genes (DEGs). The protein biosynthesis and metabolism played important roles in wound healing and siphon regeneration of R. philippinarum. Transcriptomic analysis identified the Wnt and TGF-β signaling pathways by focusing on the function and expression pattern of genes in these pathways during siphon regeneration. In addition, we carried out a genome-wide identification and phylogenetic analysis of TGF-β superfamily in R. philippinarum. The expression profiles of the TGF-β superfamily genes were analyzed in eight adult tissues (adductor muscle, mantle, foot, gill, siphon, digestive gland, gonad, and labial palp) and regenerative siphon. This study shed new light on the process of morphological regeneration and regenerative mechanism of siphon of R. philippinarum.

Genome-wide identification and expression profiling of TYR gene family in Ruditapes philippinarum under the challenge of Vibrio anguillarum

Tyrosinase (EC1.14.18.1, TYR) is also called phenol oxidase, is not only involved in pigmentation but also plays an important role in modulating innate immunity in invertebrates. Tyrosinase is a copper containing metalloenzyme. The tyrosinase protein has two copper binding sites and three conserved histidines. In this study, 21 tyrosinase genes (RpTYR) were obtained from the whole genome of Ruditapes philippinarum. Their open reading frames were from 951 to 5424 aa, the range of predicted relative molecular weight from 36.72 to 203.81 kDa, and the range of isoelectric point from 4.72 to 9.88. Transcriptome analysis showed that RpTYR gene was expressed specifically in different developmental stages, adult tissues, four strains and two groups with different shell colors. Besides, the expression profiles of 21 RpTYRs were investigated against the immune response of R. philippinarum to a Vibrio challenge. The qPCR results showed that RpTYRs were involved in the immune response of R. philippinarum after Vibrio anguillarum challenge. This study provides preliminary evidence that the tyrosinases genes are involved in the immune defense and the potential immune function of R. philippinarum. Overall, these findings suggested that the expansion of TYR genes may play vital roles in larval development, the formation of shell color pattern, and immune response in R. philippinarum.

The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia

BACKGROUND

Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa.

RESULTS

Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors.

CONCLUSIONS

Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.