Using a metabolomics approach to investigate the sensitivity of a potential Arctic-invader and its Arctic sister-species to marine heatwaves and traditional harvesting disturbances

Coastal species are threatened by fishing practices and changing environmental conditions, such as marine heatwaves (MHW). The mechanisms that confer tolerance to such stressors in marine invertebrates are poorly understood. However, differences in tolerance among different species may be attributed to their geographical distribution. To test the tolerance of species occupying different thermal ranges, we used two closely related bivalves the softshell clam Mya arenaria (Linnaeus, 1758), a cold-temperate invader with demonstrated potential for establishment in the Arctic, and the blunt gaper Mya truncata (Linnaeus, 1758), a native polar species. Clams were subjected to a thermal stress, mimicking a MHW, and harvesting stress in a controlled environment. Seven acute temperature changes (2, 7, 12, 17, 22, 27, and 32 °C) were tested at two harvesting disturbance intensities (with, without). Survival was measured after 12 days and three tissues (gills, mantle, and posterior adductor muscle) collected from surviving individuals for targeted metabolomic profiling. MHW tolerance differed significantly between species: 26.9 °C for M. arenaria and 17.8 °C for M. truncata, with a negligeable effect of harvesting. At the upper thermal limit, M. arenaria displayed a more profound metabolomic remodelling when compared to M. truncata, and this varied greatly between tissue types. Network analysis revealed differences in pathway utilization at the upper MHW limit, with M. arenaria displaying a greater reliance on multiple DNA repair and expression and cell signalling pathways, while M. truncata was limited to fewer pathways. This suggests that M. truncata is ill equipped to cope with warming environments. MHW patterning in the Northwest Atlantic may be a strong predictor of population survival and future range shifts in these two clam species. As polar environments undergo faster rates of warming compared to the global average, M. truncata may be outcompeted by M. arenaria expanding into its native range.

Centuries of genome instability and evolution in soft-shell clam, Mya arenaria, bivalve transmissible neoplasia

Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. We investigated the evolutionary history of a cancer lineage that has spread though the soft-shell clam (Mya arenaria) population by assembling a chromosome-scale soft-shell clam reference genome and characterizing somatic mutations in transmissible cancer. We observe high mutation density, widespread copy-number gain, structural rearrangement, loss of heterozygosity, variable telomere lengths, mitochondrial genome expansion and transposable element activity, all indicative of an unstable cancer genome. We also discover a previously unreported mutational signature associated with overexpression of an error-prone polymerase and use this to estimate the lineage to be >200 years old. Our study reveals the ability for an invertebrate cancer lineage to survive for centuries while its genome continues to structurally mutate, likely contributing to the evolution of this lineage as a parasitic cancer.

The Genome of the Softshell Clam Mya arenaria and the Evolution of Apoptosis

Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

Population Genetic Structure is Unrelated to Shell Shape, Thickness and Organic Content in European Populations of the Soft-Shell Clam Mya Arenaria

The soft-shell clam Mya arenaria is one of the most ancient invaders of European coasts and is present in many coastal ecosystems, yet little is known about its genetic structure in Europe. We collected 266 samples spanning a latitudinal cline from the Mediterranean to the North Sea and genotyped them at 12 microsatellite loci. In parallel, geometric morphometric analysis of shell outlines was used to test for associations between shell shape, latitude and genotype, and for a selection of shells we measured the thickness and organic content of the granular prismatic (PR), the crossed-lamellar (CL) and the complex crossed-lamellar (CCL) layers. Strong population structure was detected, with Bayesian cluster analysis identifying four groups located in the Mediterranean, Celtic Sea, along the continental coast of the North Sea and in Scotland. Multivariate analysis of shell shape uncovered a significant effect of collection site but no associations with any other variables. Shell thickness did not vary significantly with either latitude or genotype, although PR thickness and calcification were positively associated with latitude, while CCL thickness showed a negative association. Our study provides new insights into the population structure of this species and sheds light on factors influencing shell shape, thickness and microstructure.

Biogeography of resistance to paralytic shellfish toxins in softshell clam, Mya arenaria (L.), populations along the Atlantic coast of North America

Blooms of Alexandrium spp., the causative agent of paralytic shellfish poisoning (PSP), recur with varying frequency and intensity on the Northwest Atlantic coast of North America, from New York, USA, to northern Canadian waters. Along this latitudinal range blooms co-occur with abundant, intertidal populations of softshell clams, Mya arenaria. Prior work identified a naturally-occurring genetic mutation in Domain II α-subunit of the clams' voltage-gated sodium channels (Na_V), which significantly reduces the binding affinity of the paralytic shellfish toxin, saxitoxin (STX). This mutation provides clams with resistance to the deleterious effects of STX, allowing them to continue feeding during Alexandrium spp. blooms and attain very high tissue toxicities. This study used genetic sequencing of the Na_V mutation locus in clams from four coastal regions of the Bay of Fundy-Gulf of Maine and the mid-Atlantic to determine the percentage of clams in each region that possess the resistant Na_V mutation. The genotype composition was related to the occurrence and magnitude of PSP outbreaks based on shellfish toxicity, primarily that of mussels, Mytilus edulis, used as a proxy for the prevalence and severity of Alexandrium blooms in each region. As hypothesized, the proportion of clams bearing the resistant mutation generally matched up well with the historical incidence and intensity of Alexandrium spp. blooms. The highest percentage of homozygote resistant clams (RR = 70.0%), and the lowest percentage of sensitive clams (SS = 4.5%) were found in eastern Gulf of Maine populations. Exceptions at a few sites where anomalously high numbers of M. arenaria with the resistant mutation were found despite the absence of blooms, may be attributable to larval gene flow. There was no evidence that Alexandrium blooms occurring in Northport Harbor, Long Island, have resulted in a shift in genotypic composition of the local clam population, presumably due to their low cell toxicity. Seasonal mismatch of highly vulnerable M. arenaria postset with toxic blooms at this latitude may also partly explain this result. This study provides strong supporting evidence that Alexandrium blooms can select for resistance to PSP-toxins in M. arenaria populations and proposes a mechanism for the persistence of the sensitive allele throughout the region. Implications for clam aquaculture (seeding) efforts, as well as for shellfish toxicity monitoring are discussed.

Cellular stress responses to chronic heat shock and shell damage in temperate Mya truncata

Acclimation, via phenotypic flexibility, is a potential means for a fast response to climate change. Understanding the molecular mechanisms underpinning phenotypic flexibility can provide a fine-scale cellular understanding of how organisms acclimate. In the last 30 years, Mya truncata populations around the UK have faced an average increase in sea surface temperature of 0.7 °C and further warming of between 1.5 and 4 °C, in all marine regions adjacent to the UK, is predicted by the end of the century. Hence, data are required on the ability of M. truncata to acclimate to physiological stresses, and most notably, chronic increases in temperature. Animals in the present study were exposed to chronic heat-stress for 2 months prior to shell damage and subsequently, only 3, out of 20 damaged individuals, were able to repair their shells within 2 weeks. Differentially expressed genes (between control and damaged animals) were functionally enriched with processes relating to cellular stress, the immune response and biomineralisation. Comparative transcriptomics highlighted genes, and more broadly molecular mechanisms, that are likely to be pivotal in this lack of acclimation. This study demonstrates that discovery-led transcriptomic profiling of animals during stress-response experiments can shed light on the complexity of biological processes and changes within organisms that can be more difficult to detect at higher levels of biological organisation.

A mu class glutathione S-transferase from Manila clam Ruditapes philippinarum (RpGSTμ): cloning, mRNA expression, and conjugation assays

Glutathione S-transferases (GSTs) are enzymes that catalyze xenobiotic metabolism in the phase II detoxification process. GSTs have a potential for use as indicators or biomarkers to assess the presence of organic and inorganic contaminants in aquatic environments. In this study, a full-length cDNA of a mu (μ) class GST (RpGSTμ) was identified from Manila clam (Ruditapes philippinarum) and biochemically characterized. The 1356 bp of the cDNA included an open reading frame of 651 bp encoding a polypeptide of 217 amino acid residues with a molecular mass of 25.04 kDa and an estimated pI of 6.34. Sequence analysis revealed that the RpGSTμ possessed several characteristic features of μ class GSTs, such as a thioredoxin-like N-terminal domain containing binding sites for glutathione (GSH), a C-terminal domain containing substrate binding sites, and a μ loop. The recombinant RpGSTμ (rRpGSTμ) protein exhibited GSH-conjugating catalytic activity towards several substrates, and significantly strong activity was detected against 4-nitrophenethyl bromide (5.77 ± 0.55) and 1-chloro-2,4-dinitrobenzene (CDNB, 3.19 ± 0.05). Kinetic analysis as a function of GSH and CDNB concentrations revealed relatively low Km values of 1.03 ± 0.46 mM and 0.56 ± 0.20 mM, respectively, thereby indicating a GSH-conjugation attributed with high rates. The optimum pH and temperature for the catalytic activity of the rRpGSTμ protein were 7.7 and 37°C, respectively. The effect of two inhibitors, Cibacron blue and hematin, on the activity of rRpGSTμ was evaluated and the IC50 values of 0.65 μM and 9 μM, respectively, were obtained. While RpGSTμ transcripts were highly expressed in gills and hemocytes, a significant elevation in mRNA levels was detected in these tissues after lipopolysaccharide (LPS), polyinosinic-polycytidylic acid (poly I:C) and live bacterial (Vibrio tapetis) challenges. These findings collectively suggest that RpGSTμ functions as a potent detoxifier of xenobiotic toxicants present in the aquatic environment, and that its mRNA expression could be modulated by pathogenic stress signal(s).

DNA integrity assessment in hemocytes of soft-shell clams (Mya arenaria) in the Saguenay Fjord (Québec, Canada)

The purpose of this study was to examine the effects of pollution on DNA integrity in the feral soft-shell clam (Mya arenaria) in the Saguenay Fjord. Intertidal clams were collected downstream and upstream of the fjord at sites under anthropogenic pollution. DNA integrity was assessed by following changes in single- and double-stranded breaks, variation in DNA content and micro-nuclei (MN) incidence in hemocytes. The results revealed that clams collected at polluted sites had reduced DNA strand breaks (lower DNA repair activity), increased DNA content variation and MN frequency in hemocytes. The data revealed that DNA content variation was closely related to MN frequency and negatively with DNA strand breaks formation. Water conductivity was also related to reduced MN frequency and DNA content variation, indicating that, in addition to the effects of pollution, the gradual dilution of saltwater could compromise mussel health.

Isolation and characterization of new microsatellite markers for the invasive softshell clam, Mya arenaria (L.) (Bivalvia: Myidae)

The invasive softshell clam (Mya arenaria Linnaeus, 1758) is native to the northwestern region of the Atlantic Ocean. This species has been introduced in the northeast Pacific and along the European coasts, due to intense naval transports and aquaculture, and it is now present in all the European seas. In this paper we describe seven new microsatellite loci for Mya arenaria. The isolated loci are polymorphic with a number of alleles per locus between 6 and 14. The observed and expected heterozygosities ranged from 0.417 to 0.951, and from 0.643 to 0.895, with an average of 0.716 and 0.775, respectively. These microsatellite markers should be useful in analyzing this species’ genetic diversity, which could explain various processes of its invasion history.

Effects of pesticide compounds (chlorothalonil and mancozeb) and benzo[a]pyrene mixture on aryl hydrocarbon receptor, p53 and ubiquitin gene expression levels in haemocytes of soft-shell clams (Mya arenaria)

The aim of this study is to investigate the effects of the pesticides/polycyclic aromatic hydrocarbon mixture on aryl hydrocarbon receptor (AhR), p53 and ubiquitin mRNA level in haemocytes of Mya arenaria exposed to a mixture of chlorothalonil, mancozeb and benzo[a]pyrene (BaP) for 48 and 72 h. AhR, p53 and ubiquitin gene expression levels were quantified using quantitative Real-time PCR. For robust and accurate quantification of transcripts, suitable housekeeping genes were selected from four sets of ribosomal and elongation factors transcripts previously sequenced from Mya arenaria using geNorm open source software. Quantitative Real-time PCR data exhibited a significantly high expression of AhR after 72 h of exposure (P ≤ 0.05). p53 gene expression seems to be up-regulated by the mixture after 48 h, however not significantly; but the level of p53 mRNA is down-regulated by the xenobiotics between 48 and 72 h after exposure. This study postulates that AhR mRNA levels could be used as an indicator of the exposure of clams' haemocytes to a mixture of xenobiotics such as chlorothalonil, mancozeb and BaP. However, further studies have to be pursued in order to unravel the molecular mechanisms involved in the p53 signaling pathway.

Induction of transposase and polyprotein RNA levels in disseminated neoplastic hemocytes of soft-shell clams: Mya arenaria

In Prince Edward Island, a high mortality of soft-shell clams Mya arenaria was found to be related to the disease known as disseminated neoplasia (DN). However, the molecular mechanisms by which hemocytes of clams are transformed in the course of DN remain by far unknown. This study aims at identifying the transcripts involved in the development of the disease. Four subtractive cDNA sequence libraries were generated and more than 200,000 reads were obtained. Following similarity searches in genome databases, the transcripts were assigned to cellular functions including mitochondrial respiration, structural proteins, cytoskeleton, nucleic acid regulation, general metabolism, signal transduction, apoptosis, cell cycle regulation, as well as virus transcripts. The expression levels of transposase and polyprotein genes were evaluated in clams with various percentages of tetraploid hemocytes. Data have shown that expression levels were significantly higher in clams with a high percentage of tetraploid hemocytes. These results reinforce the hypothesis of endogenous retrotransposon involvement in the etiology of the disease. Further investigations are needed, however, to elucidate the role of transposase and polyprotein in the disease development.

Identification and expression of immune-related genes in hemocytes of soft-shell clams, Mya arenaria, challenged with Vibrio splendidus

Although the mollusc immune system has been studied at the cellular level, the response to pathogens at gene expression level has not been thoroughly investigated. This study aimed to investigate the early molecular response of hemocytes of soft-shell clams, Mya arenaria, to Vibrio splendidus strain LGP32 by identification of transcripts involved in immune defense. The Suppression Subtractive Hybridization (SSH) was used to selectively identify differentially expressed genes in hemocytes exposed to V. splendidus at a ratio 1:1 for 2 h. Both forward and reverse subtracted cDNA were constructed and a total of 16,000 reads were obtained and analyzed. Identity searches in genome databases were performed using BlastX program and transcripts were clustered to cellular functions including structural proteins, immunity, stress proteins, apoptosis, cell process, metabolism and signal transduction. Among the differentially expressed immune associated genes were ficolin, killer cell lectin-like receptor, natural resistance-associated macrophage protein 1 (Nramp-1), mitogen-activated protein kinases (MAPK), ferritin, heat shock proteins 90 (HSP90) and cathepsin and their expressions were quantified using Reverse Transcription quantitative Polymerase Chain Reaction (RT-qPCR) at 1, 2 and 3 h post-Vibrio challenge. These genes showed similar expression patterns, up-regulation at 1 h, followed by a down-regulation at 2 and 3 h. These data corroborates our previous observations of cell rounding, reduced phagocytosis and respiratory burst activity. To our knowledge, this is the first study to demonstrate an effect of V. splendidus on expression of genes related to immune system in soft-shell clams M. arenaria. However, further investigations are needed to unravel the molecular mechanisms of hemocytes subjected to V. splendidus.

Differential gene expression of gamma-actin, Toll-like receptor 2 (TLR-2) and interleukin-1 receptor-associated kinase 4 (IRAK-4) in Mya arenaria haemocytes induced by in vivo infections with two Vibrio splendidus strains

Immune function gene expression in Mya arenaria haemocytes was evaluated following in vivo infection with Vibrio splendidus LGP32-GFP and 7SHRW. Elongation factor 1alpha (EF-1alpha) with 2 (EF-2), after challenge with LGP32-GFP, and EF-1alpha with the ribosomal protein S-18, after challenge with 7SHRW, were found to be the most stable housekeeping genes. Using these internal controls and comparing the regulation induced by both strains, up-regulation of gamma-actin, down-regulation of TLR-2 and up-regulation of IRAK-4 was significantly higher after challenge with LGP32-GFP (p<0.001, p=0.001 and p<0.05, respectively). These results suggest specific responses at a molecular level modulated by the bacterial strains. LGP32-GFP induced marked responses which coincide with a similar trend previously found on phenotypic responses under our experimental model.

Selecting a set of housekeeping genes for quantitative real-time PCR in normal and tetraploid haemocytes of soft-shell clams, Mya arenaria

The transcripts involved in the molecular mechanisms of haemic neoplasia in relation to the haemocyte ploidy status of the soft-shell clam, Mya arenaria, have yet to be identified. For this purpose, real-time quantitative RT-PCR constitutes a sensitive and efficient technique, which can help determine the gene expression involved in haemocyte tetraploid status in clams affected by haemic neoplasia. One of the critical steps in comparing transcription profiles is the stability of selected housekeeping genes, as well as an accurate normalization. In this study, we selected five reference genes, S18, L37, EF1, EF2 and actin, generally used as single control genes. Their expression was analyzed by real-time quantitative RT-PCR at different levels of haemocyte ploidy status in order to select the most stable genes. Using the geNorm software, our results showed that L37, EF1 and S18 represent the most stable gene expressions related to various ploidy status ranging from 0 to 78% of tetraploid haemocytes in clams sampled in North River (Prince Edward Island, Canada). However, actin gene expression appeared to be highly regulated. Hence, using it as a housekeeping gene in tetraploid haemocytes can result in inaccurate data. To compare gene expression levels related to haemocyte ploidy status in Mya arenaria, using L37, EF1 and S18 as housekeeping genes for accurate normalization is therefore recommended.

Adaptive evolution of tetrodotoxin resistance in animals

Tetrodotoxin (TTX), first isolated from pufferfish (tetraodontids), is a highly potent neurotoxin that selectively binds to voltage-gated sodium channels (Na(v)) in muscle and nerve tissues causing paralysis and death. Saxitoxin (STX) is a TTX-related neurotoxin produced by dinoflagellates. Recent investigations have implicated diverse substitutions in the P-loop regions of skeletal muscle and neuronal Na(v) channels in the convergent evolution of neurotoxin resistance in pufferfish, garter snakes and softshell clams, which has enabled them to feed on TTX- and STX-bearing organisms.