Impaired reproduction, energy reserves and dysbiosis: The overlooked consequences of heatwaves in a bivalve mollusc.

Extreme events like Marine Heatwaves (MHWs) are becoming more intense, severe, and frequent, threatening benthic communities, specifically bivalves. However, the consequences of non-lethal MHWs on animals are still poorly understood. Here, we exposed the Manila clam Ruditapes philippinarum to non-lethal MHW for 30 days and provided an integrative view of its effects. Our result indicated that albeit non-lethal, MHW reduced clam's energy reserves (by reducing their hepato-somatic index), triggered antioxidant defenses (particularly in males), impaired reproduction (via the production of smaller oocytes in females), triggered dysbiosis in the digestive gland microbiota and altered animals' behaviour (by impacting their burying capacity) and filtration rate. Such effects were seen also at RNA-seq (i.e. many down-regulated genes belonged to reproduction) and metabolome level. Interestingly, negative effects were more pronounced in males than in females. Our results show that MHWs influence animal physiology at multiple levels, likely impacting its fitness and its ecosystem services.

Multidisciplinary long-term survey of Manila clam grown in farming sites subjected to different environmental conditions.

In recent years recurrent bivalve mass mortalities considerably increased around the world, causing the collapse of natural and farmed populations. Venice Lagoon has historically represented one of the major production areas of the Manila clam Ruditapes philippinarum in Europe. However, in the last 20 years a 75 % decrease in the annual production has been experienced. While climate change and anthropogenic interventions may have played a key role in natural and farmed stocks reductions, no studies investigated at multiple levels the environmental stressors affecting farmed Manila clam to date. In this work we carried out a long-term monitoring campaign on Manila clam reared in four farming sites located at different distances from the southern Venice Lagoon inlet, integrating (meta)genomic approaches (i.e. RNA-seq; microbiota characterization), biometric measurements and chemical-physical parameters. Our study allowed to characterize the molecular mechanisms adopted by this species to cope with the different environmental conditions characterizing farming sites and to propose hypotheses to explain mortality events observed in recent years. Among the most important findings, the disruption of clam's immune response, the spread of Vibrio spp., and the up-regulation of molecular pathways involved in xenobiotic metabolism suggested major environmental stressors affecting clams farmed in sites placed close to Chioggia's inlet, where highest mortality was also observed. Overall, our study provides knowledge-based tools for managing Manila clam farming on-growing areas. In addition, the collected data is a snapshot of the time immediately before the commissioning of MoSE, a system of mobile barriers aimed at protecting Venice from high tides, and will represent a baseline for future studies on the effects of MoSE on clams farming and more in general on the ecology of the Venice Lagoon.

Co-infection of two eukaryotic pathogens within clam populations in Arcachon Bay.

The parasitic species Perkinsus olseni (= atlanticus) (Perkinsea, Alveolata) infects a wide range of mollusc species and is responsible for mortality events and economic losses in the aquaculture industry and fisheries worldwide. Thus far, most studies conducted in this field have approached the problem from a "one parasite-one disease" perspective, notably with regards to commercially relevant clam species, while the impact of other Perkinsus species should also be considered as it could play a key role in the disease phenotype and dynamics. Co-infection of P. olseni and P. chesapeaki has already been sporadically described in Manila clam populations in Europe. Here, we describe for the first time the parasitic distribution of two Perkinsus species, P. olseni and P. chesapeaki, in individual clam organs and in five different locations across Arcachon Bay (France), using simultaneous in situ detection by quantitative PCR (qPCR) duplex methodology. We show that P. olseni single-infection largely dominated prevalence (46-84%) with high intensities of infection (7.2 to 8.5 log-nb of copies. g-1of wet tissue of Manila clam) depending on location, suggesting that infection is driven by the abiotic characteristics of stations and physiological states of the host. Conversely, single P. chesapeaki infections were observed in only two sampling stations, Ile aux Oiseaux and Gujan, with low prevalences 2 and 14%, respectively. Interestingly, the co-infection by both Perkinsus spp., ranging in prevalence from 12 to 34%, was distributed across four stations of Arcachon Bay, and was detected in one or two organs maximum. Within these co-infected organs, P. olseni largely dominated the global parasitic load. Hence, the co-infection dynamics between P. olseni and P. chesapeaki may rely on a facilitating role of P. olseni in developing a primary infection which in turn may help P. chesapeaki infect R. philippinarum as a reservoir for a preferred host. This ecological study demonstrates that the detection and quantification of both parasitic species, P. olseni and P. chesapeaki, is essential and timely in resolving cryptic infections and their consequences on individual hosts and clam populations.

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.

Long-lasting effects of chronic exposure to chemical pollution on the hologenome of the Manila clam.

Chronic exposure to pollutants affects natural populations, creating specific molecular and biochemical signatures. In the present study, we tested the hypothesis that chronic exposure to pollutants might have substantial effects on the Manila clam hologenome long after removal from contaminated sites. To reach this goal, a highly integrative approach was implemented, combining transcriptome, genetic and microbiota analyses with the evaluation of biochemical and histological profiles of the edible Manila clam Ruditapes philippinarum, as it was transplanted for 6 months from the polluted area of Porto Marghera (PM) to the clean area of Chioggia (Venice lagoon, Italy). One month post-transplantation, PM clams showed several modifications to its resident microbiota, including an overrepresentation of the opportunistic pathogen Arcobacter spp. This may be related to the upregulation of several immune genes in the PM clams, potentially representing a host response to the increased abundance of deleterious bacteria. Six months after transplantation, PM clams demonstrated a lower ability to respond to environmental/physiological stressors related to the summer season, and the hepatopancreas-associated microbiota still showed different compositions among PM and CH clams. This study confirms that different stressors have predictable effects in clams at different biological levels and demonstrates that chronic exposure to pollutants leads to long-lasting effects on the animal hologenome. In addition, no genetic differentiation between samples from the two areas was detected, confirming that PM and CH clams belong to a single population. Overall, the obtained responses were largely reversible and potentially related to phenotypic plasticity rather than genetic adaptation. The results here presented will be functional for the assessment of the environmental risk imposed by chemicals on an economically important bivalve species.

Potential for Genetic Improvement of Resistance to Perkinsus olseni in the Manila Clam, Ruditapes philippinarum, Using DNA Parentage Assignment and Mass Spawning.

The Manila clam Ruditapes philippinarum, a major cultured shellfish species, is threatened by infection with the microparasite Perkinsus olseni, whose prevalence increases with high water temperatures. Under the current trend of climate change, the already severe effects of this parasitic infection might rapidly increase the frequency of mass mortality events. Treating infectious diseases in bivalves is notoriously problematic, therefore selective breeding for resistance represents a key strategy for mitigating the negative impact of pathogens. A crucial step in initiating selective breeding is the estimation of genetic parameters for traits of interest, which relies on the ability to record parentage and accurate phenotypes in a large number of individuals. Here, to estimate the heritability of resistance against P. olseni, a field experiment mirroring conditions in industrial clam production was set up, a genomic tool was developed for parentage assignment, and parasite load was determined through quantitative PCR. A mixed-family cohort of potentially 1,479 clam families was produced in a hatchery by mass spawning of 53 dams and 57 sires. The progenies were seeded in a commercial clam production area in the Venice lagoon, Italy, where high prevalence of P. olseni had previously been reported. Growth and parasite load were monitored every month and, after 1 year, more than 1,000 individuals were collected for DNA samples and phenotype recording. A pooled sequencing approach was carried out using DNA samples from the hatchery broodstock and from a Venice lagoon clam population, providing candidate markers used to develop a 245-SNP panel. Parentage assignment for 246 F1 individuals showed sire and dam representation were high (75 and 85%, respectively), indicating a very limited risk of inbreeding. Moderate heritability (0.23 ± 0.11-0.35 ± 0.13) was estimated for growth traits (shell length, shell weight, total weight), while parasite load showed high heritability, estimated at 0.51 ± 0.20. No significant genetic correlations were found between growth-associated traits and parasite load. Overall, the preliminary results provided by this study show high potential for selecting clams resistant to parasite load. Breeding for resistance may help limit the negative effects of climate change on clam production, as the prevalence of the parasite is predicted to increase under a future scenario of higher temperatures. Finally, the limited genetic correlation between resistance and growth suggests that breeding programs could incorporate dual selection without negative interactions.

Depuration processes affect the Vibrio community in the microbiota of the Manila clam, Ruditapes philippinarum.

As filter-feeders, bivalve molluscs accumulate Vibrio into edible tissues. Consequently, an accurate assessment of depuration procedures and the characterization of the persistent Vibrio community in depurated shellfish represent a key issue to guarantee food safety in shellfish products. The present study investigated changes in the natural Vibrio community composition of the Ruditapes philippinarum microbiota with specific focus on human pathogenic species. For this purpose, the study proposed a MLSA-NGS approach (rRNA 16S, recA and pyrH) for the detection and identification of Vibrio species. Clam microbiota were analysed before and after depuration procedures performed in four depuration plants, using culture-dependent and independent approaches. Microbiological counts and NGS data revealed differences in terms of both contamination load and Vibrio community between depuration plants. The novel MLSA-NGS approach allowed for a clear definition of the Vibrio species specific to each depuration plant. Specifically, depurated clam microbiota showed presence of human pathogenic species. Ozone treatments and the density of clams in the depuration tank probably influenced the level of contamination and the Vibrio community composition. The composition of Vibrio community specific to each plant should be carefully evaluated during the risk assessment to guarantee a food-safe shellfish-product for the consumer.

Host-microbiota interactions shed light on mortality events in the striped venus clam Chamelea gallina.

Mass mortalities due to disease outbreaks have recently affected a number of major taxa in marine ecosystems. Climate- and pollution-induced stress may compromise host immune defenses, increasing the risk of opportunistic diseases. Despite growing evidence that mass mortality events affecting marine species worldwide are strongly influenced by the interplay of numerous environmental factors, the reductionist approaches most frequently used to investigate these factors hindered the interpretation of these multifactorial pathologies. In this study, we propose a broader approach based on the combination of RNA-sequencing and 16S microbiota analyses to decipher the factors underlying mass mortality in the striped venus clam, Chamelea gallina, along the Adriatic coast. On one hand, gene expression profiling and functional analyses of microbial communities showed the over-expression of several genes and molecular pathways involved in xenobiotic metabolism, suggesting potential chemical contamination in mortality sites. On the other hand, the down-regulation of several genes involved in immune and stress response, and the over-representation of opportunistic pathogens such as Vibrio and Photobacterium spp. indicates that these microbial species may take advantage of compromised host immune pathways and defense mechanisms that are potentially affected by chemical exposure, resulting in periodic mortality events. We propose the application of our approach to interpret and anticipate the risks inherent in the combined effects of pollutants and microbes on marine animals in today's rapidly changing environment.