Chromosome-Level Clam Genome Helps Elucidate the Molecular Basis of Adaptation to a Buried Lifestyle

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Comprehensive investigation of differentially expressed ncRNAs, mRNAs, and their ceRNA networks in the regulation of shell color formation in clam, Cyclina sinensis

Noncoding RNAs (ncRNAs) have gained significant attention in recent years due to their crucial roles in various biological processes. However, our understanding of the expression and functions of ncRNAs in Cyclina sinensis, an economically important marine bivalve, remains limited. This study aimed to address this knowledge gap by systematically identifying ncRNAs in the mantles of C. sinensis with purple and white shells. Through our analysis, we identified a differential expression of 1244 mRNAs, 196 lncRNAs, 49 circRNAs, and 23 miRNAs between purple- and white-shell clams. Functional enrichment analysis revealed the involvement of these differentially expressed ncRNAs in biomineralization and pigmentation processes. To gain further insights into the regulatory mechanisms underlying shell color formation, we established competitive endogenous RNA (ceRNA) networks. These networks allowed us to identify targeted differentially expressed miRNAs (DEMis) and genes associated with shell color formation. Based on the ceRNA networks, we obtained an up-down-up lncRNA-miRNA-mRNA network consisting of 13 upregulated lncRNAs and a circRNA-miRNA-mRNA network comprising three upregulated circRNAs (novel_circ_0004787, novel_circ_0001165, novel_circ_0000245). Through these networks, we identified and selected an upregulated novel gene (evm.TU.Hic_asm_7.988) and a downregulated sponge miRNA (hru-miR-1985) as potential contributors to shell color regulation. In summary, the present investigation offers a comprehensive analysis of ncRNA catalogs expressed in the clam mantle of C. sinensis. The findings enhance our comprehension of the molecular mechanisms governing shell coloration and offer new perspectives for selective breeding of C. sinensis in the future.

Comparative transcriptome elucidates key genes and pathways related to golden phenotype of Crassostrea gigas

Marine bivalves are economically important and exhibit a remarkable diversity in shell color. The Pacific oyster Crassostrea gigas stands out as an important economic species, with the successful development of four distinct color strains through selective breeding. While previous studies have shed light on the genetic mechanism underlying color segregation, the precise molecular regulatory mechanisms responsible for shell coloration in oysters remains elusive. In this study, we confirmed that the golden phenotype is primarily attributed to pheomelanin by histological and ultrastructural observations. Additionally, we conducted a comparative transcriptome analysis of the black and golden shell color oysters to explore the potential genes and pathways contributing to the golden phenotype in C. gigas. Our results revealed a significant increase in differentially expressed genes in the golden phenotype associated with pathways such as glutathione metabolism, and calcium signaling pathway, suggesting a potential role in the synthesis of pheomelanin. Of particular note, we highlighted the potential role of two-pore channel 2 (TPC2) in modulating tyrosinase activity and melanosomal pH, ultimately determining the shade of pigmentation. Our study in this work provided a preliminary exploration of the mechanism, shedding light on the melanosome microenvironment and shell color.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

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.

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.

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.

A chromosome-level genome of the kuruma shrimp (Marsupenaeus japonicus) provides insights into its evolution and cold-resistance mechanism

Marsupenaeus japonicus is an important marine crustacean species. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for genome-assisted breeding. Consequently, we determined the chromosome-level genome of M. japonicus. Here we determine the chromosome-level genome assembly for M. japonicus with a total of 665.19 Gb genomic sequencing data, yielding an approximately1.54 Gb assembly with a contig N50 size of 229.97 kb and a scaffold N50 size of 38.27 Mb. With the high-throughput chromosome conformation capture (Hi-C) technology, we anchored 18,019 contigs onto 42 pseudo-chromosomes, accounting for 99.40% of the total genome assembly. Analysis of the present M. japonicus genome revealed 24,317 protein-coding genes and a high proportion of repetitive sequences (61.56%). The high-quality genome assembly enabled the identification of genes associated with cold-stress and cold tolerance in kuruma shrimp through the comparison of eyestalk transcriptomes between the low temperature-stressed shrimp (10 °C) and normal temperature shrimp (28 °C). The genome assembly presented here could be useful in future studies to reveal the molecular mechanisms of M. japonicus in response to low temperature stress and the molecular assisted breeding of M. japonicus in low temperature.

Dynamic evolutionary history of spiralian-specific TALE homeobox genes in mollusks

Homeobox genes play essential roles in the early development of many animals. Although the repertoire of most homeobox genes, including three amino acid loop extension (TALE)-type homeobox genes, is conserved in animals, spiralian-TALE (SPILE) genes are a notable exception. In this study, SPILE genes were extracted from the genomic data of 22 mollusk species and classified into four clades (-A/C, -B, -D, and -E) to determine which SPILE genes exhibit dynamic repertoire changes. While SPILE-D and -E duplications were rarely observed, SPILE-B duplication was observed in the bivalve lineage and SPILE-A/C duplication was observed in multiple clades. Conversely, most or all SPILE genes were lost in cephalopods and in some gastropod lineages. SPILE gene expression patterns were also analyzed in multiple mollusk species using publicly available RNA-seq data. The majority of SPILE genes examined, particularly those in the A/C- and B-clades, were specifically expressed during early development, suggesting that most SPILE genes exert specific roles in early development. This comprehensive cataloging and characterization revealed a dynamic evolutionary history, including SPILE-A/C and -B gene duplications and the loss of SPILE genes in several lineages. Furthermore, this study provides a useful resource for studying the molecular mechanism of spiralian early development and the evolution of young and lineage-specific transcription factors.

Genome-wide analysis of the hard clam mitogen-activated protein kinase kinase gene family and their transcriptional profiles under abiotic stress

Mitogen-activated protein kinase kinase (MAPKK) was the hub component of the Mitogen-activated protein kinase (MAPK) signaling pathway and played an important role in the cellular response to environmental stress. In this study, we identified five MmMAPKK genes in hard clam Mercenaria mercenaria and found that all MmMAPKK genes contain a conserved protein kinase domain. The MmMAPKK genes derived from dispersed duplication were unevenly distributed in three chromosomes. Although the genome size was highly variable among different bivalve mollusks, the number of MAPKK genes was relatively stable. Phylogenetic analysis showed that bivalve MAPKK was divided into five clades, and amino acid sequences of MAPKK from the same clade consisted of similar conserved motifs. The syntenic analysis demonstrated that MmMAPKKs had the highest number of homologous gene pairs with Cyclina sinensis. MmMAPKKs were ubiquitously expressed in all examined tissues, and all MmMAPKK genes were highly expressed in the ovary. MmMAPKK genes showed stress-specific expression under envirionmental stress. MmMAPKK7 showed an upregulated in heat and heat plus hypoxia stress while MmMAPKK1 showed an upregulated in hypoxic stress groups. Dynamic changes of MmMAPKK7, MmMAPKK6 and MmMAPKK1 in hemocytes were observed in response to air exposure. MmMAPKK4 significantly downregulated after air exposure for five days. MmMAPKK7 and MmMAPKK6 might participate in adaptation to low salinity stress. Our results provided useful information about MAPKK and laid a foundation for further studies on MAPKK evolution in the bivalve.

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.

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.

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.

Evo-devo of shell colour in gastropods and bivalves

Recent technical innovations are revealing surprising patterns in mollusc shell pigmentation, such as an unexpectedly modest role for melanins and rapid divergences in the mix of pigments used to achieve similar colour patterns. The elucidation of the molecular genetic basis of shell pigmentation has been slow, probably because of the high genome complexity of gastropods and bivalves. Recent work within the old field of evolutionary ecology of shell pigmentation allows a greater role for the analysis of large-geographic-scale patterns (sometimes employing citizen-science data), as well as experimental field studies. However, the field remains dominated by land snails as model organisms, while colour pattern evolution in marine gastropods and bivalves, particularly those not exposed to visual predators, remains mysterious.

Reconstruction of ancient homeobox gene linkages inferred from a new high-quality assembly of the Hong Kong oyster (Magallana hongkongensis) genome

Background

Homeobox-containing genes encode crucial transcription factors involved in animal, plant and fungal development, and changes to homeobox genes have been linked to the evolution of novel body plans and morphologies. In animals, some homeobox genes are clustered together in the genome, either as remnants from ancestral genomic arrangements, or due to coordinated gene regulation. Consequently, analyses of homeobox gene organization across animal phylogeny provide important insights into the evolution of genome organization and developmental gene control, and their interaction. However, homeobox gene organization remains to be fully elucidated in several key animal ancestors, including those of molluscs, lophotrochozoans and bilaterians.

Results

Here, we present a high-quality chromosome-level genome assembly of the Hong Kong oyster, Magallana hongkongensis (2n = 20), for which 93.2% of the genomic sequences are contained on 10 pseudomolecules (~ 758 Mb, scaffold N50 = 72.3 Mb). Our genome assembly was scaffolded using Hi-C reads, facilitating a larger scaffold size compared to the recently published M. hongkongensis genome of Peng et al. (Mol Ecol Resources, 2020), which was scaffolded using the Crassostrea gigas assembly. A total of 46,963 predicted gene models (45,308 protein coding genes) were incorporated in our genome, and genome completeness estimated by BUSCO was 94.6%. Homeobox gene linkages were analysed in detail relative to available data for other mollusc lineages.

Conclusions

The analyses performed in this study and the accompanying genome sequence provide important genetic resources for this economically and culturally valuable oyster species, and offer a platform to improve understanding of animal biology and evolution more generally. Transposable element content is comparable to that found in other mollusc species, contrary to the conclusion of another recent analysis. Also, our chromosome-level assembly allows the inference of ancient gene linkages (synteny) for the homeobox-containing genes, even though a number of the homeobox gene clusters, like the Hox/ParaHox clusters, are undergoing dispersal in molluscs such as this oyster.

Identification of variants associated with hard clam, Mercenaria mercenaria, resistance to Quahog Parasite Unknown disease

Severe losses in aquacultured and wild hard clam (Mercenaria mercenaria) stocks have been previously reported in the northeastern United States due to a protistan parasite called QPX (Quahog Parasite Unknown). Previous work demonstrated that clam resistance to QPX is under genetic control. This study identifies single nucleotide polymorphism (SNP) associated with clam survivorship from two geographically segregated populations, both deployed in an enzootic site. The analysis contrasted samples collected before and after undergoing QPX-related mortalities and relied on a robust draft clam genome assembly. ~200 genes displayed significant variant enrichment at each sampling point in both populations, including 18 genes shared between both populations. Markers from both populations were identified in genes related to apoptosis pathways, protein-protein interaction, receptors, and signaling. This research begins to identify genetic markers associated with clam resistance to QPX disease, leading the way for the development of resistant clam stocks through marker-assisted selection.