Adaptations to nitrogen availability drive ecological divergence of chemosynthetic symbionts.

Bacterial symbionts, with their shorter generation times and capacity for horizontal gene transfer (HGT), play a critical role in allowing marine organisms to cope with environmental change. The closure of the Isthmus of Panama created distinct environmental conditions in the Tropical Eastern Pacific (TEP) and Caribbean, offering a "natural experiment" for studying how closely related animals evolve and adapt under environmental change. However, the role of bacterial symbionts in this process is often overlooked. We sequenced the genomes of endosymbiotic bacteria in two sets of sister species of chemosymbiotic bivalves from the genera Codakia and Ctena (family Lucinidae) collected on either side of the Isthmus, to investigate how differing environmental conditions have influenced the selection of symbionts and their metabolic capabilities. The lucinid sister species hosted different Candidatus Thiodiazotropha symbionts and only those from the Caribbean had the genetic potential for nitrogen fixation, while those from the TEP did not. Interestingly, this nitrogen-fixing ability did not correspond to symbiont phylogeny, suggesting convergent evolution of nitrogen fixation potential under nutrient-poor conditions. Reconstructing the evolutionary history of the nifHDKT operon by including other lucinid symbiont genomes from around the world further revealed that the last common ancestor (LCA) of Ca. Thiodiazotropha lacked nif genes, and populations in oligotrophic habitats later re-acquired the nif operon through HGT from the Sedimenticola symbiont lineage. Our study suggests that HGT of the nif operon has facilitated niche diversification of the globally distributed Ca. Thiodiazotropha endolucinida species clade. It highlights the importance of nitrogen availability in driving the ecological diversification of chemosynthetic symbiont species and the role that bacterial symbionts may play in the adaptation of marine organisms to changing environmental conditions.

Eyeless razor clam Sinonovacula constricta discriminates light spectra through opsins to guide Ca2+ and cAMP signaling pathways.

Phototransduction is based on opsins that drive distinct types of Gα cascades. Although nonvisual photosensitivity has long been known in marine bivalves, the underlying molecular basis and phototransduction mechanism are poorly understood. Here, we introduced the eyeless razor clam Sinonovacula constricta as a model to clarify this issue. First, we showed that S. constricta was highly diverse in opsin family members, with a significant expansion in xenopsins. Second, the expression of putative S. constricta opsins was highly temporal-spatio specific, indicating their potential roles in S. constricta development and its peripheral photosensitivity. Third, by cloning four S. constricta opsins with relatively higher expression (Sc_opsin1, 5, 7, and 12), we found that they exhibited different expression levels in response to different light environments. Moreover, we demonstrated that these opsins (excluding Sc_opsin7) couple with Gαq and Gαi cascades to mediate the light-dependent Ca2+ (Sc_opsin1 and 5) and cAMP (Sc_opsin12) signaling pathways. The results indicated that Sc_opsin1 and 5 belonged to Gq-opsins, Sc_opsin12 belonged to Gi-opsins, while Sc_opsin7 might act as a photo-isomerase. Furthermore, we found that the phototransduction function of S. constricta Gq-opsins was dependent on the lysine at the seventh transmembrane domain, and greatly influenced by the external light spectra in a complementary way. Thus, a synergistic photosensitive system mediated by opsins might exist in S. constricta to rapidly respond to the transient or subtle changes of the external light environment. Collectively, our findings provide valuable insights into the evolution of opsins in marine bivalves and their potential functions in nonvisual photosensitivity.

The mud-dwelling clam Meretrix petechialis secretes endogenously synthesized erythromycin.

Although lacking an adaptive immune system and often living in habitats with dense and diverse bacterial populations, marine invertebrates thrive in the presence of potentially challenging microbial pathogens. However, the mechanisms underlying this resistance remain largely unexplored and promise to reveal novel strategies of microbial resistance. Here, we provide evidence that a mud-dwelling clam, Meretrix petechialis, synthesizes, stores, and secretes the antibiotic erythromycin. Liquid chromatography coupled with mass spectrometry, immunocytochemistry, fluorescence in situ hybridization, RNA interference, and enzyme-linked immunosorbent assay revealed that this potent macrolide antimicrobial, thought to be synthesized only by microorganisms, is produced by specific mucus-rich cells beneath the clam's mantle epithelium, which interfaces directly with the bacteria-rich environment. The antibacterial activity was confirmed by bacteriostatic assay. Genetic, ontogenetic, phylogenetic and genomic evidence, including genotypic segregation ratios in a family of full siblings, gene expression in clam larvae, phylogenetic tree, and synteny conservation in the related genome region further revealed that the genes responsible for erythromycin production are of animal origin. The detection of this antibiotic in another clam species showed that the production of this macrolide is not exclusive to M. petechialis and may be a common strategy among marine invertebrates. The finding of erythromycin production by a marine invertebrate offers a striking example of convergent evolution in secondary metabolite synthesis between the animal and bacterial domains. These findings open the possibility of engineering-animal tissues for the localized production of an antibacterial secondary metabolite.

Genome-wide association analysis reveals the genetic basis of thermal tolerance in dwarf surf clam Mulinia lateralis.

Recently, elevated seawater temperatures have resulted numerous adverse effects, including significant mortality among bivalves. The dwarf surf clam, Mulinia lateralis, is considered a valuable model species for bivalve research due to its rapid growth and short generation time. The successful cultivation in laboratory setting throughout its entire life cycle makes it an ideal candidate for exploring the potential mechanisms underlying bivalve responses to thermal stress. In this study, a total of 600 clams were subjected to a 17-day thermal stress experiment at a temperature of 30 °C which is the semi-lethal temperature for this species. Ninety individuals who perished initially were classified as heat-sensitive populations (HSP), while 89 individuals who survived the experiment were classified as heat-tolerant populations (HTP). Subsequently, 179 individuals were then sequenced, and 21,292 single nucleotide polymorphisms (SNPs) were genotyped for downstream analysis. The heritability estimate for survival status was found to be 0.375 ± 0.127 suggesting a genetic basis for thermal tolerance trait. Furthermore, a genome-wide association study (GWAS) identified three SNPs and 10 candidate genes associated with thermal tolerance trait in M. lateralis. These candidate genes were involved in the ETHR/EHF signaling pathway and played pivotal role in signal sensory, cell adhesion, oxidative stress, DNA damage repair, etc. Additionally, qPCR results indicated that, excluding MGAT4A, ZAN, and RFC1 genes, all others exhibited significantly higher expression in the HTP (p < 0.05), underscoring the critical involvement of the ETHR/EHF signaling pathway in M. lateralis' thermal tolerance. These results unveil the presence of standing genetic variations associated with thermal tolerance in M. lateralis, highlighting the regulatory role of the ETHR/EHF signaling pathway in the bivalve's response to thermal stress, which contribute to comprehension of the genetic basis of thermal tolerance in bivalves.

SmithHunter: a workflow for the identification of candidate smithRNAs and their targets.

BACKGROUND: SmithRNAs (Small MITochondrial Highly-transcribed RNAs) are a novel class of small RNA molecules that are encoded in the mitochondrial genome and regulate the expression of nuclear transcripts. Initial evidence for their existence came from the Manila clam Ruditapes philippinarum, where they have been described and whose activity has been biologically validated through RNA injection experiments. Current evidence on the existence of these RNAs in other species is based only on small RNA sequencing. As a preliminary step to characterize smithRNAs across different metazoan lineages, a dedicated, unified, analytical workflow is needed. RESULTS: We propose a novel workflow specifically designed for smithRNAs. Sequence data (from small RNA sequencing) uniquely mapping to the mitochondrial genome are clustered into putative smithRNAs and prefiltered based on their abundance, presence in replicate libraries and 5' and 3' transcription boundary conservation. The surviving sequences are subsequently compared to the untranslated regions of nuclear transcripts based on seed pairing, overall match and thermodynamic stability to identify possible targets. Ample collateral information and graphics are produced to help characterize these molecules in the species of choice and guide the operator through the analysis. The workflow was tested on the original Manila clam data. Under basic settings, the results of the original study are largely replicated. The effect of additional parameter customization (clustering threshold, stringency, minimum number of replicates, seed matching) was further evaluated. CONCLUSIONS: The study of smithRNAs is still in its infancy and no dedicated analytical workflow is currently available. At its core, the SmithHunter workflow builds over the bioinformatic procedure originally applied to identify candidate smithRNAs in the Manila clam. In fact, this is currently the only evidence for smithRNAs that has been biologically validated and, therefore, the elective starting point for characterizing smithRNAs in other species. The original analysis was readapted using current software implementations and some minor issues were solved. Moreover, the workflow was improved by allowing the customization of different analytical parameters, mostly focusing on stringency and the possibility of accounting for a minimal level of genetic differentiation among samples.

The genome of a giant clam zooxanthella (Cladocopium infistulum) offers few clues to adaptation as an extracellular symbiont with high thermotolerance.

BACKGROUND: Cladocopium infistulum (Symbiodiniaceae) is a dinoflagellate specialized to live in symbiosis with western Pacific giant clams (Tridacnidae). Unlike coral-associated symbionts, which reside within the host cells, C. infistulum inhabits the extracellular spaces of the clam's digestive diverticula. It is phylogenetically basal to a large species complex of stress-tolerant Cladocopium, many of which are associated with important reef-building corals in the genus Porites. This close phylogenetic relationship may explain why C. infistulum exhibits high thermotolerance relative to other tridacnid symbionts. Moreover, past analyses of microsatellite loci indicated that Cladocopium underwent whole-genome duplication prior to the adaptive radiations that led to its present diversity. RESULTS: A draft genome assembly of C. infistulum was produced using long- and short-read sequences to explore the genomic basis for adaptations underlying thermotolerance and extracellular symbiosis among dinoflagellates and to look for evidence of genome duplication. Comparison to three other Cladocopium genomes revealed no obvious over-representation of gene groups or families whose functions would be important for maintaining C. infistulum's unique physiological and ecological properties. Preliminary analyses support the existence of partial or whole-genome duplication among Cladocopium, but additional high-quality genomes are required to substantiate these findings. CONCLUSION: Although this investigation of Cladocopium infistulum revealed no patterns diagnostic of heat tolerance or extracellular symbiosis in terms of overrepresentation of gene functions or genes under selection, it provided a valuable genomic resource for comparative analyses. It also indicates that ecological divergence among Cladocopium species, and potentially among other dinoflagellates, is partially governed by mechanisms other than gene content. Thus, additional high-quality, multiomic data are needed to explore the molecular basis of key phenotypes among symbiotic microalgae.

Mitonuclear Sex Determination? Empirical Evidence from Bivalves.

Genetic elements encoded in nuclear DNA determine the sex of an individual in many animals. In certain bivalve lineages that possess doubly uniparental inheritance (DUI), mitochondrial DNA (mtDNA) has been hypothesized to contribute to sex determination. In these cases, females transmit a female mtDNA to all offspring, while male mtDNA (M mtDNA) is transmitted only from fathers to sons. Because M mtDNA is inherited in the same way as Y chromosomes, it has been hypothesized that mtDNA may be responsible for sex determination. However, the role of mitochondrial and nuclear genes in sex determination has yet to be validated in DUI bivalves. In this study, we used DNA, RNA, and mitochondrial short noncoding RNA (sncRNA) sequencing to explore the role of mitochondrial and nuclear elements in the sexual development pathway of the freshwater mussel Potamilus streckersoni (Bivalvia: Unionida). We found that the M mtDNA sheds a sncRNA partially within a male-specific mitochondrial gene that targets a pathway hypothesized to be involved in female development and mitophagy. RNA-seq confirmed the gene target was significantly upregulated in females, supporting a direct role of mitochondrial sncRNAs in gene silencing. These findings support the hypothesis that M mtDNA inhibits female development. Genome-wide patterns of genetic differentiation and heterozygosity did not support a nuclear sex-determining region, although we cannot reject that nuclear factors are involved with sex determination. Our results provide further evidence that mitochondrial loci contribute to diverse, nonrespiratory functions and additional insights into an unorthodox sex-determining system.

Bivalves Present the Largest and Most Diversified Repertoire of Toll-Like Receptors in the Animal Kingdom, Suggesting Broad-Spectrum Pathogen Recognition in Marine Waters.

Toll-like receptors (TLRs) are the most widespread class of membrane-bound innate immune receptors, responsible of specific pathogen recognition and production of immune effectors through the activation of intracellular signaling cascades. The repertoire of TLRs was analyzed in 85 metazoans, enriched on molluscan species, an underrepresented phylum in previous studies. Following an ancient evolutionary origin, suggested by the presence of TLR genes in Anthozoa (Cnidaria), these receptors underwent multiple independent gene family expansions, the most significant of which occurred in bivalve molluscs. Marine mussels (Mytilus spp.) had the largest TLR repertoire in the animal kingdom, with evidence of several lineage-specific expanded TLR subfamilies with different degrees of orthology conservation within bivalves. Phylogenetic analyses revealed that bivalve TLR repertoires were more diversified than their counterparts in deuterostomes or ecdysozoans. The complex evolutionary history of TLRs, characterized by lineage-specific expansions and losses, along with episodic positive selection acting on the extracellular recognition domains, suggests that functional diversification might be a leading evolutionary force. We analyzed a comprehensive transcriptomic data set from Mytilus galloprovincialis and built transcriptomic correlation clusters with the TLRs expressed in gills and in hemocytes. The implication of specific TLRs in different immune pathways was evidenced, as well as their specific modulation in response to different biotic and abiotic stimuli. We propose that, in a similar fashion to the remarkable functional specialization of vertebrate TLRs, the expansion of the TLR gene family in bivalves attends to a functional specification motivated by the biological particularities of these organisms and their living environment.

Comparative transcriptomics revealed parallel evolution and innovation of photosymbiosis molecular mechanisms in a marine bivalve.

Photosymbioses between heterotrophic hosts and autotrophic symbionts are evolutionarily prevalent and ecologically significant. However, the molecular mechanisms behind such symbioses remain less elucidated, which hinders our understanding of their origin and adaptive evolution. This study compared gene expression patterns in a photosymbiotic bivalve (Fragum sueziense) and a closely related non-symbiotic species (Trigoniocardia granifera) under different light conditions to detect potential molecular pathways involved in mollusc photosymbiosis. We discovered that the presence of algal symbionts greatly impacted host gene expression in symbiont-containing tissues. We found that the host immune functions were suppressed under normal light compared with those in the dark. In addition, we found that cilia in the symbiont-containing tissues play important roles in symbiont regulation or photoreception. Interestingly, many potential photosymbiosis genes could not be annotated or do not exhibit orthologues in T. granifera transcriptomes, indicating unique molecular functions in photosymbiotic bivalves. Overall, we found both novel and known molecular mechanisms involved in animal-algal photosymbiosis within bivalves. Given that many of the molecular pathways are shared among distantly related host lineages, such as molluscs and cnidarians, it indicates that parallel and/or convergent evolution is instrumental in shaping host-symbiont interactions and responses in these organisms.

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia).

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

eDNA-based detection reveals invasion risks of a biofouling bivalve in the world's largest water diversion project.

Environmental DNA (eDNA) has increasingly been used to detect rare species (e.g., newly introduced nonindigenous species) in both terrestrial and aquatic ecosystems, often with distinct advantages over traditional methods. However, whether water eDNA signals can be used to inform invasion risks remains debatable owing to inherent uncertainties associated with the methods used and the varying conditions among study systems. Here, we sampled eDNA from canals of the central route of the South-to-North Water Diversion Project (hereafter SNWDP) in China to investigate eDNA distribution and efficacy to inform invasion risks in a unique lotic system. We first conducted a total of 16 monthly surveys in this system (two sites in the source reservoir and four sites in the main canal) to test if eDNA could be applied to detect an invasive, biofouling bivalve, the golden mussel Limnoperna fortunei. Second, we initiated a one-time survey in a sub-canal of the SNWDP using refined sampling (12 sites in ~22 km canal) and considered a few environmental predictors. We found that detection of target eDNA in the main canal was achieved up to 1100 km from the putative source population but was restricted to the warmer months (May-November). Detection probability exhibited a significant positive relationship with average daily minimum air temperature and with water temperature, consistent with the expected spawning season. eDNA concentration in the main canal generally fluctuated across months and sites and was generally higher in warmer months. Golden mussel eDNA concentration in the sub-canal decreased significantly with distance from the source and with increasing water temperature and became almost undetectable at ~22 km distance. Given the enormity of the SNWDP, golden mussels may eventually expand their distribution in the main canal, with established "bridgehead" populations facilitating further spread. Our findings suggest an elevated invasion risk of golden mussels in the SNWDP in warm months, highlighting the critical period for spread and, possibly, management.

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

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

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.

Effects of acute hypoxia and reoxygenation on histological structure, antioxidant response, and apoptosis in razor clam Sinonovacula constricta.

Hypoxia is one of the major environmental problems limiting the healthy development of intensive aquaculture. Marine benthic shellfish are encountering heightened problems related to hypoxic stress as a result of ongoing human activities and aquaculture operations. Razor clam Sinonovacula constricta, a commercially valuable shellfish, has not yet been reported in studies on physiological changes caused by hypoxia and reoxygenation. To understand the negative effects of hypoxia and reoxygenation on the clams, we set up two low-oxygen concentration groups (DO 2.0 mg/L and DO 0.5 mg/L) and assessed multiple aspects of oxidative damage to their hepatopancreas and gills. After the hypoxic stress, the two tissues of the razor clam suffered varying degrees of damage, including cell degeneration and disruption of mitochondrial cristae. After reoxygenation, the 2.0 mg/L group recovered substantially, but the clams in the 0.5 mg/L group still unrecovered. The activities of antioxidant enzymes (MDA, T-AOC, SOD, GPX, and CAT) in clams were considerably altered by acute hypoxia and reoxygenation. Briefly, there was a growing and then declining trend in MDA, T-AOC, and SOD activities in the hepatopancreas, whereas GPX and CAT activities showed the converse trend. In the hepatopancreas and gills, the level of anti-apoptotic gene Bcl-2 transcripts gradually decreased with the duration of hypoxia and increased following reoxygenation. However, changes in the transcript level of the pro-apoptotic gene Bax were in contrast to that of Bcl-2. The TUNEL assay revealed that hypoxia caused apoptosis. Furthermore, at DO 0.5 mg/L, the degree of apoptosis was more significant than at DO 2.0 mg/L, and hepatopancreatic apoptosis was more severe than gill apoptosis. Collectively, our findings imply that hypoxia induces oxidative stress, histological damage, and apoptosis in razor clams in a concentration-dependent and tissue-specific manner. These consequences serve as a reminder that prolonged recovery periods may be required for razor clams to fully recover from oxidative damage resulting from hypoxia-reoxygenation episodes.

The combination of high temperature and Vibrio infection worsens summer mortality in the clam Meretrix petechialis by increasing apoptosis and oxidative stress.

The interaction between environmental factors and Vibrio in bivalves is not well understood, despite the widely held belief that pathogen infection and seawater temperature significantly impact summer mortality. In the present study, we conducted simulated experiments to explore the effects of high temperature and Vibrio infection on the clam Meretrix petechialis. The survival curve analysis revealed that the combined challenge of high temperature and Vibrio infection (31°C-vibrio) led to significantly higher clam mortality compared to the groups exposed solely to Vibrio (27°C-vibrio), high temperature (31°C-control), and the control condition (27°C-control). Furthermore, PCoA analysis of 11 immune genes indicated that Vibrio infection predominated during the incubation period, with a gradual equilibrium between these factors emerging during the course of the infection. Additionally, our investigations into apoptosis and autophagy processes exhibited significant induction of mTOR and Bcl2 of the 31°C-vibrio group in the early challenge stage, followed by inhibition in the later stage. Oxidative stress analysis demonstrated a substantial additive effect on malondialdehyde (MDA) and glutathione (GSH) content in the combined challenge group compared to the control group. Comparative transcriptome analysis revealed a significant increase in differentially expressed genes related to immunity, such as complement C1q-like protein, C-type lectin, big defensin, and lysozyme, in the 31°C-vibrio group, suggesting that the synergistic effect of high temperature and Vibrio infection triggers more robust antibacterial immune responses. These findings provide critical insights for understanding the infection process and uncovering the causes of summer mortality.

Molecular identification and functional characterization of a C-type lectin gene in Meretrix meretrix.

C-type lectins (CTLs) are a kind of Ca2+-dependent immunoreactive factors, which participated in pathogens recognition and defense. The present study identified a new CTL from hard clam Meretrix meretrix (designated as MmCTL4). The full-length of MmCTL4 cDNA was 608 bp, encoding a presumed signal peptide of 19 bp and a carbohydrate recognition domain (CRD) of 131 bp. The tertiary structure of recombinant MmCTL4 protein (rMmCTL4) was the typical long double-ring structure with three conserved disulfide bonds, and the motifs in Ca2+-binding sites of MmCTL4 were QPN and WSD. The SYBR Green real-time PCR analysis indicated that MmCTL4 was widely expressed in the hemocytes, hepatopancreas and mantle of healthy clams. After Vibrio splendidus stimulation, the temporal expression profile of MmCTL4 mRNA in hemocytes and hepatopancreas increased by 7.8-fold at 6 hpi and 3.9-fold at 12 hpi, respectively. The cDNA fragments encoding MmCTL4 were recombined into pET-32a (+) vectors, and transformed into Escherichia coli BL21 (DE3). The rMmCTL4 with the presence of Ca2+ performed obvious hemagglutination activity, and could agglutinate E. coli, Bacillus subtilis, and Staphylococcus aureus, while it only weakly agglutinate Vibrio parahaemolyticus and fungi P. pastoris. The agglutination activity of rMmCTL4 were significantly inhibited by D-mannose, D-xylose, D-lactose, maltose and lipopolysaccharides. These results indicated that MmCTL4, as a class of typical pattern recognition receptors (PRRs), could protect the host against pathogen invasion in the innate immunity of clams.

Development and validation of genome-wide SSR molecular markers of Tapes dorsatus.

BACKGROUND: Tapes dorsatus is an economically important benthic animal in the Beibu Gulf of China. However, the deficiency of microsatellite markers has hindered the study of its genetics. The development of microsatellite markers will provide useful tools for genetic improvement, variety identification, phylogenetic analysis and resource conservation. METHODS AND RESULTS: Within the genome sequence, 145,008 simple sequence repeats (SSRs) were identified, and 29,691 primer pairs were designed successfully. A total of 100 primer pairs were randomly synthesized for testing, and 93 primers yielded products. Sixty highly polymorphic primers were used to reveal the genetic diversity of 50 T. dorsatus individuals. The average number of alleles (Na) of the population was 10.40; the average number of effective alleles was 6.16, the average expected heterozygosity (He) was 0.82, and the average polymorphic information content was 0.80. The genetic structure of the population was detected, by which the population could be divided into three subpopulations. CONCLUSION: We identified 145,008 SSRs in the genome of T. dorsatus and designed 29,691 primer pairs in this study. Of 100 synthesized primers, 60 were highly polymorphic and used to reveal the genetic diversity and structure of the population. The SSR markers identified here will provide useful tools and a foundation for genetic diversity, linkage mapping and molecular marker-aided breeding in T. dorsatus.

Glacial vicariance and secondary contact shape demographic histories in a freshwater mussel species complex.

Characterizing the mechanisms influencing the distribution of genetic variation in aquatic species can be difficult due to the dynamic nature of hydrological landscapes. In North America's Central Highlands, a complex history of glacial dynamics, long-term isolation, and secondary contact have shaped genetic variation in aquatic species. Although the effects of glacial history have been demonstrated in many taxa, responses are often lineage- or species-specific and driven by organismal ecology. In this study, we reconstruct the evolutionary history of a freshwater mussel species complex using a suite of mitochondrial and nuclear loci to resolve taxonomic and demographic uncertainties. Our findings do not support Pleurobema rubrum as a valid species, which is proposed for listing as threatened under the U.S. Endangered Species Act. We synonymize P. rubrum under Pleurobema sintoxia-a common and widespread species found throughout the Mississippi River Basin. Further investigation of patterns of genetic variation in P. sintoxia identified a complex demographic history, including ancestral vicariance and secondary contact, within the Eastern Highlands. We hypothesize these patterns were shaped by ancestral vicariance driven by the formation of Lake Green and subsequent secondary contact after the last glacial maximum. Our inference aligns with demographic histories observed in other aquatic taxa in the region and mirrors patterns of genetic variation of a freshwater fish species (Erimystax dissimilis) confirmed to serve as a parasitic larval host for P. sintoxia. Our findings directly link species ecology to observed patterns of genetic variation and may have significant implications for future conservation and recovery actions of freshwater mussels.

Transcriptomic responses of Antarctic clam Laternula elliptica to nanoparticles, at single and combined exposures reveal ecologically relevant biomarkers.

In recent years micro- and nanoplastics and metal-oxide nanomaterials have been found in several environmental compartments. The Antarctic soft clam Laternula elliptica is an endemic Antarctic species having a wide distribution in the Southern Ocean. Being a filter-feeder, it could act as suitable bioindicator of pollution from nanoparticles also considering its sensitivity to various sources of stress. The present study aims to assess the impact of polystyrene nanoparticles (PS-NP) and the nanometal titanium-dioxide (n-TiO2) on genome-wide transcript expression of L. elliptica either alone and in combination and at two toxicological relevant concentrations (5 and 50 µg/L) during 96 h exposure. Transcript-target qRT-PCR was performed with the aim to identify suitable biomarkers of exposure and effects. As expected, at the highest concentration tested, the clustering was clearer between control and exposed clams. A total of 221 genes resulted differentially expressed in exposed clams and control ones, and 21 of them had functional annotation such as ribosomal proteins, antioxidant, ion transport (osmoregulation), acid-base balance, immunity, lipid metabolism, cell adhesion, cytoskeleton, apoptosis, chromatin condensation and cell signaling. At functional level, relevant transcripts were shared among some treatments and could be considered as general stress due to nanoparticle exposure. After applying transcript-target approach duplicating the number of clam samples, four ecologically relevant transcripts were revealed as biomarkers for PS-NP, n-TiO2 and their combination at 50 µg/L, that could be used for monitoring clams' health status in different Antarctic localities.

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.