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

Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica)

Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.

A common garden experiment supports a genetic component underlying the increased resilience of common cockle (Cerastoderma edule) to the parasite Marteilia cochillia

The common cockle is a valuable bivalve species inhabiting the Atlantic European coasts. The parasite Marteilia cochillia has devastated cockle beds in the southern Galician (NW Spain) rias since 2012. Previous data suggested that cockles from Ría de Arousa acquired some resilience to this parasite through natural selection after consecutive annual marteiliosis outbreaks and candidate markers associated with marteiliosis resilience were identified using population genomics and transcriptomics approaches. Here, a common garden experiment was performed using a naïve stock (from Ría de Muros-Noia) and an affected stock (from Ría de Arousa) to test this hypothesis. Breeders from both stocks were used to produce seed cohorts at hatchery, which were pre-grown in a raft (outdoor nursery stage) and deployed in two shellfish beds affected by marteiliosis in Ría de Arousa (growing-out stage). In both beds, the naïve stock showed high marteiliosis prevalence and was fully depleted in a short period, while the affected stock barely showed evidence of marteiliosis. A set of 45 SNPs putatively associated with marteiliosis resilience were fitted for MassARRAY genotyping to check their role in the differential resilience detected between both stocks. Though no significant differentiation was found between the naïve and the affected stocks with neutral markers, 28 SNPs showed significant divergence between them, suggesting that these SNPs were involved in directional selection during eight generations (to the most) of marteiliosis pressure (long-term selection). Furthermore, signals of selection were also detected in the naïve stock along the marteiliosis outbreak in the growing-out stage (short-term selection) and six SNPs, all shared with the long-term evaluation, showed consistent signals of differentiation according to the infection severity. Some of these SNPs were located within immune genes pertaining to families such as proteasome, ubiquitin, tumor necrosis factor, and glutathione S-transferase. These resilience-associated markers will be useful to recover cockle production in Galicia.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Molecular Features Associated with Resilience to Ocean Acidification in the Northern Quahog, Mercenaria mercenaria

The increasing concentration of CO₂ in the atmosphere and resulting flux into the oceans will further exacerbate acidification already threatening coastal marine ecosystems. The subsequent alterations in carbonate chemistry can have deleterious impacts on many economically and ecologically important species including the northern quahog (Mercenaria mercenaria). The accelerated pace of these changes requires an understanding of how or if species and populations will be able to acclimate or adapt to such swift environmental alterations. Thus far, studies have primarily focused on the physiological effects of ocean acidification (OA) on M. mercenaria, including reductions in growth and survival. However, the molecular mechanisms of resilience to OA in this species remains unclear. Clam gametes were fertilized under normal pCO₂ and reared under acidified (pH ~ 7.5, pCO₂ ~ 1200 ppm) or control (pH ~ 7.9, pCO₂ ~ 600 ppm) conditions before sampled at 2 days (larvae), 32 days (postsets), 5 and 10 months (juveniles) and submitted to RNA and DNA sequencing to evaluate alterations in gene expression and genetic variations. Results showed significant shift in gene expression profiles among clams reared in acidified conditions as compared to their respective controls. At 10 months of exposure, significant shifts in allele frequency of single nucleotide polymorphisms (SNPs) were identified. Both approaches highlighted genes coding for proteins related to shell formation, bicarbonate transport, cytoskeleton, immunity/stress, and metabolism, illustrating the role these pathways play in resilience to OA.

Transcriptomic, Proteomic, and Functional Assays Underline the Dual Role of Extrapallial Hemocytes in Immunity and Biomineralization in the Hard Clam Mercenaria mercenaria

Circulating hemocytes in the hemolymph represent the backbone of innate immunity in bivalves. Hemocytes are also found in the extrapallial fluid (EPF), the space delimited between the shell and the mantle, which is the site of shell biomineralization. This study investigated the transcriptome, proteome, and function of EPF and hemolymph in the hard clam Mercenaria mercenaria. Total and differential hemocyte counts were similar between EPF and hemolymph. Overexpressed genes in the EPF were found to have domains previously identified as being part of the "biomineralization toolkit" and involved in bivalve shell formation. Biomineralization related genes included chitin-metabolism genes, carbonic anhydrase, perlucin, and insoluble shell matrix protein genes. Overexpressed genes in the EPF encoded proteins present at higher abundances in the EPF proteome, specifically those related to shell formation such as carbonic anhydrase and insoluble shell matrix proteins. Genes coding for bicarbonate and ion transporters were also overexpressed, suggesting that EPF hemocytes are involved in regulating the availability of ions critical for biomineralization. Functional assays also showed that Ca²⁺ content of hemocytes in the EPF were significantly higher than those in hemolymph, supporting the idea that hemocytes serve as a source of Ca²⁺ during biomineralization. Overexpressed genes and proteins also contained domains such as C1q that have dual functions in biomineralization and immune response. The percent of phagocytic granulocytes was not significantly different between EPF and hemolymph. Together, these findings suggest that hemocytes in EPF play a central role in both biomineralization and immunity.

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.

An Overview of Transcriptional Responses of Schistosome-Susceptible (M line) or -Resistant (BS-90) Biomphalaria glabrata Exposed or Not to Schistosoma mansoni Infection

Background

We seek to provide a comprehensive overview of transcriptomics responses of immune-related features of the gastropod Biomphalaria glabrata (Bg) following exposure to Schistosoma mansoni (Sm), a trematode causing human schistosomiasis. Responses of schistosome-susceptible (M line, or SUS) and -resistant (BS-90, or RES) Bg strains were characterized following exposure to Sm for 0.5, 2, 8 or 40 days post-exposure (dpe).

Methods

RNA-Seq and differential expression analysis were undertaken on 56 snails from 14 groups. We considered 7 response categories: 1) constitutive resistance factors; 2) constitutive susceptibility factors; 3) generalized stress responses; 4) induced resistance factors; 5) resistance factors suppressed in SUS snails; 6) suppressed/manipulated factors in SUS snails; and 7) tolerance responses in SUS snails. We also undertook a gene co-expression network analysis. Results from prior studies identifying schistosome resistance/susceptibility factors were examined relative to our findings.

Results

A total of 792 million paired-end reads representing 91.2% of the estimated 31,985 genes in the Bg genome were detected and results for the 7 categories compiled and highlighted. For both RES and SUS snails, a single most supported network of genes with highly correlated expression was found.

Conclusions

1) Several constitutive differences in gene expression between SUS and RES snails were noted, the majority over-represented in RES; 2) There was little indication of a generalized stress response shared by SUS and RES snails at 0.5 or 2 dpe; 3) RES snails mounted a strong, multi-faceted response by 0.5 dpe that carried over to 2 dpe; 4) The most notable SUS responses were at 40 dpe, in snails shedding cercariae, when numerous features were either strongly down-regulated indicative of physiological distress or parasite manipulation, or up-regulated, suggestive of tolerance or survival-promoting effects; 5) Of 55 genes previously identified in genome wide mapping studies, 29 (52.7%) were responsive to Sm, as were many familiar resistance-associated genes (41.0%) identified by other means; 6) Both network analysis and remarkably specific patterns of expression of lectins and G protein-coupled receptors in categories 4, 6 and 7 were indicative of orchestrated responses of different suites of genes in SUS or RES snails following exposure to Sm.