Contaminants from dredged sediments alter the transcriptome of Manila clam and induce shifts in microbiota composition.

BACKGROUND: The reuse of dredged sediments in ports and lagoons is a big issue as it should not affect the quality and the equilibrium of ecosystems. In the lagoon of Venice, sediment management is of crucial importance as sediments are often utilized to built-up structures necessary to limit erosion. However, the impact of sediment reuse on organisms inhabiting this delicate area is poorly known. The Manila clam is a filter-feeding species of high economic and ecological value for the Venice lagoon experiencing a drastic decline in the last decades. In order to define the molecular mechanisms behind sediment toxicity, we exposed clams to sediments sampled from different sites within one of the Venice lagoon navigable canals close to the industrial area. Moreover, we investigated the impacts of dredged sediments on clam's microbial communities. RESULTS: Concentrations of the trace elements and organic chemicals showed increasing concentrations from the city of Venice to sites close to the industrial area of Porto Marghera, where PCDD/Fs and PCBs concentrations were up to 120 times higher than the southern lagoon. While bioaccumulation of organic contaminants of industrial origin reflected sediments' chemical concentrations, metal bioaccumulation was not consistent with metal concentrations measured in sediments probably due to the activation of ABC transporters. At the transcriptional level, we found a persistent activation of the mTORC1 signalling pathway, which is central in the coordination of cellular responses to chemical stress. Microbiota characterization showed the over-representation of potential opportunistic pathogens following exposure to the most contaminated sediments, leading to host immune response activation. Despite the limited acquisition of new microbial species from sediments, the latter play an important role in shaping Manila clam microbial communities. CONCLUSIONS: Sediment management in the Venice lagoon will increase in the next years to maintain and create new canals as well as to allow the operation of the new mobile gates at the three Venice lagoon inlets. Our data reveal important transcriptional and microbial changes of Manila clams after exposure to sediments, therefore reuse of dredged sediments represents a potential risk for the conservation of this species and possibly for other organisms inhabiting the Venice lagoon.

Evaluating Escherichia coli contamination in bivalve mollusks using the impedance method: a comparison with most probable number analyses and correlation with environmental parameters.

The application of an electrochemical (impedance) tool for monitoring Escherichia coli contamination in shellfish was evaluated after 13 months of observation. The primary aim of the present study was to compare the standard most probable number (MPN) and µ-trac 4200 (log imped/100 g) for the assessment of E. coli contamination (log MPN/100 g) in non-depurated bivalve mollusks (BM) from five sampling areas of the Veneto-Emilian coast (Italy) (118 samples). The secondary aim was to evaluate the correlation between E. coli concentrations in BM and environmental factors on a large data set (690). The methods showed a moderate, positive correlation (0.60 and 0.69 Pearson and Spearman coefficients, respectively; P<0.01) in Ruditapes philippinarum. The McNemar test indicated analogous sample classification between methods, and the impedance method overestimated the most contaminated class (P=0.03; >4,600 MPN/100 g). The results highlighted the suitability of the impedance method for a faster evaluation and routine use especially in clams, while in Mytilus it seemed less effective. Different models built by multivariate permutational variance analysis and multinomial logistic regression selected the suitable environmental features able to predict the E. coli load. Overall, salinity and season affected the E. coli contamination, whereas locally it was mainly influenced by hydrometry and salinity. The application of the impedance method coupled with environmental data analysis could help purification phase management to adhere to legal limits and could represent an advantage for local control authorities to define actions, considering extreme meteorological events' effects as a proactive reaction to climate change.

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.

Viral nervous necrosis resistance in gilthead sea bream (Sparus aurata) at the larval stage: heritability and accuracy of genomic prediction with different training and testing settings.

BACKGROUND: The gilthead sea bream (Sparus aurata) has long been considered resistant to viral nervous necrosis (VNN), until recently, when significant mortalities caused by a reassortant nervous necrosis virus (NNV) strain were reported. Selective breeding to enhance resistance against NNV might be a preventive action. In this study, 972 sea bream larvae were subjected to a NNV challenge test and the symptomatology was recorded. All the experimental fish and their parents were genotyped using a genome-wide single nucleotide polymorphism (SNP) array consisting of over 26,000 markers. RESULTS: Estimates of pedigree-based and genomic heritabilities of VNN symptomatology were consistent with each other (0.21, highest posterior density interval at 95% (HPD95%): 0.1-0.4; 0.19, HPD95%: 0.1-0.3, respectively). The genome-wide association study suggested one genomic region, i.e., in linkage group (LG) 23 that might be involved in sea bream VNN resistance, although it was far from the genome-wide significance threshold. The accuracies (r) of the predicted estimated breeding values (EBV) provided by three Bayesian genomic regression models (Bayes B, Bayes C, and Ridge Regression) were consistent and on average were equal to 0.90 when assessed in a set of cross-validation (CV) procedures. When genomic relationships between training and testing sets were minimized, accuracy decreased greatly (r = 0.53 for a validation based on genomic clustering, r = 0.12 for a validation based on a leave-one-family-out approach focused on the parents of the challenged fish). Classification of the phenotype using the genomic predictions of the phenotype or using the genomic predictions of the pedigree-based, all data included, EBV as classifiers was moderately accurate (area under the ROC curve 0.60 and 0.66, respectively). CONCLUSIONS: The estimate of the heritability for VNN symptomatology indicates that it is feasible to implement selective breeding programs for increased resistance to VNN of sea bream larvae/juveniles. Exploiting genomic information offers the opportunity of developing prediction tools for VNN resistance, and genomic models can be trained on EBV using all data or phenotypes, with minimal differences in classification performance of the trait phenotype. In a long-term view, the weakening of the genomic ties between animals in the training and test sets leads to decreased genomic prediction accuracies, thus periodical update of the reference population with new data is mandatory.

An innovative index to incorporate transcriptomic data into weight of evidence approaches for environmental risk assessment.

The sharp decrease in the cost of RNA-sequencing and the rapid improvement in computational analysis of eco-toxicogenomic data have brought new insights into the adverse effects of chemicals on aquatic organisms. Yet, transcriptomics is generally applied qualitatively in environmental risk assessments, hampering more effective exploitation of this evidence through multidisciplinary studies. In view of this limitation, a methodology is here presented to quantitatively elaborate transcriptional data in support to environmental risk assessment. The proposed methodology makes use of results from the application of Gene Set Enrichment Analysis to recent studies investigating the response of Mytilus galloprovincialis and Ruditapes philippinarum exposed to contaminants of emerging concern. The degree of changes in gene sets and the relevance of physiological reactions are integrated in the calculation of a hazard index. The outcome is then classified according to five hazard classes (from absent to severe), providing an evaluation of whole-transcriptome effects of chemical exposure. The application to experimental and simulated datasets proved that the method can effectively discriminate different levels of altered transcriptomic responses when compared to expert judgement (Spearman correlation coefficient of 0.96). A further application to data collected in two independent studies of Salmo trutta and Xenopus tropicalis exposed to contaminants confirmed the potential extension of the methodology to other aquatic species. This methodology can serve as a proof of concept for the integration of "genomic tools" in environmental risk assessment based on multidisciplinary investigations. To this end, the proposed transcriptomic hazard index can now be incorporated into quantitative Weight of Evidence approaches and weighed, with results from other types of analysis, to elucidate the role of chemicals in adverse ecological effects.

Functional Evolution of Clustered Aquaporin Genes Reveals Insights into the Oceanic Success of Teleost Eggs.

Aquaporin-mediated oocyte hydration is considered important for the evolution of pelagic eggs and the radiative success of marine teleosts. However, the molecular regulatory mechanisms controlling this vital process are not fully understood. Here, we analyzed >400 piscine genomes to uncover a previously unknown teleost-specific aquaporin-1 cluster (TSA1C) comprised of tandemly arranged aqp1aa-aqp1ab2-aqp1ab1 genes. Functional evolutionary analysis of the TSA1C reveals a ∼300-million-year history of downstream aqp1ab-type gene loss, neofunctionalization, and subfunctionalization, but with marine species that spawn highly hydrated pelagic eggs almost exclusively retaining at least one of the downstream paralogs. Unexpectedly, one-third of the modern marine euacanthomorph teleosts selectively retain both aqp1ab-type channels and co-evolved protein kinase-mediated phosphorylation sites in the intracellular subdomains together with teleost-specific Ywhaz-like (14-3-3ζ-like) binding proteins for co-operative membrane trafficking regulation. To understand the selective evolutionary advantages of these mechanisms, we show that a two-step regulated channel shunt avoids competitive occupancy of the same plasma membrane space in the oocyte and accelerates hydration. These data suggest that the evolution of the adaptive molecular regulatory features of the TSA1C facilitated the rise of pelagic eggs and their subsequent geodispersal in the oceanic currents.

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.

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.

Machine Learning and Canine Chronic Enteropathies: A New Approach to Investigate FMT Effects.

Fecal microbiota transplantation (FMT) represents a very promising approach to decreasing disease activity in canine chronic enteropathies (CE). However, the relationship between remission mechanisms and microbiome changes has not been elucidated yet. The main objective of this study was to report the clinical effects of oral freeze-dried FMT in CE dogs, comparing the fecal microbiomes of three groups: pre-FMT CE-affected dogs, post-FMT dogs, and healthy dogs. Diversity analysis, differential abundance analysis, and machine learning algorithms were applied to investigate the differences in microbiome composition between healthy and pre-FMT samples, while Canine Chronic Enteropathy Clinical Activity Index (CCECAI) changes and microbial diversity metrics were used to evaluate FMT effects. In the healthy/pre-FMT comparison, significant differences were noted in alpha and beta diversity and a list of differentially abundant taxa was identified, while machine learning algorithms predicted sample categories with 0.97 (random forest) and 0.87 (sPLS-DA) accuracy. Clinical signs of improvement were observed in 74% (20/27) of CE-affected dogs, together with a statistically significant decrease in CCECAI (median value from 5 to 2 median). Alpha and beta diversity variations between pre- and post-FMT were observed for each receiver, with a high heterogeneity in the response. This highlighted the necessity for further research on a larger dataset that could identify different healing patterns of microbiome changes.

Genomic Predictions of Phenotypes and Pseudo-Phenotypes for Viral Nervous Necrosis Resistance, Cortisol Concentration, Antibody Titer and Body Weight in European Sea Bass.

In European sea bass (Dicentrarchus labrax L.), the viral nervous necrosis mortality (MORT), post-stress cortisol concentration (HC), antibody titer (AT) against nervous necrosis virus and body weight (BW) show significant heritability, which makes selective breeding a possible option for their improvement. An experimental population (N = 650) generated by a commercial broodstock was phenotyped for the aforementioned traits and genotyped with a genome-wide SNP panel (16,075 markers). We compared the predictive accuracies of three Bayesian models (Bayes B, Bayes C and Bayesian Ridge Regression) and a machine-learning method (Random Forest). The prediction accuracy of the EBV for MORT was approximately 0.90, whereas the prediction accuracies of the EBV and the phenotype were 0.86 and 0.21 for HC, 0.79 and 0.26 for AT and 0.71 and 0.38 for BW. The genomic prediction of the EBV for MORT used to classify the phenotype for the same trait showed moderate classification performance. Genome-wide association studies confirmed the polygenic nature of MORT and demonstrated a complex genetic structure for HC and AT. Genomic predictions of the EBV for MORT could potentially be used to classify the phenotype of the same trait, though further investigations on a larger experimental population are needed.

Molecular Monitoring of SARS-CoV-2 in Different Sewage Plants in Venice and the Implications for Genetic Surveillance.

Wastewater-based epidemiology is now widely used as an indirect tool to monitor the spread of SARS-CoV-2. In this study, five different sample matrices representing diverse phases of the wastewater treatment process were collected during the second wave of SARS-CoV-2 from two wastewater treatment plants (WWTPs) serving the Civil Hospital and Sacca Fisola island in Venice, Italy. Positive SARS-CoV-2 detections occurred at both WWTPs, and data on viral genome detection rate and quantification suggest that the pellet (i.e., the particulate resulting from the influent) is a sensitive matrix that permits reliable assessment of infection prevalence while reducing time to results. On the contrary, analysis of post-treatment matrices provides evidence of the decontamination efficacy of both WWTPs. Finally, direct sequencing of wastewater samples enabled us to identify B.1.177 and B.1.160 as the prevalent SARS-CoV-2 lineages circulating in Venice at the time of sampling. This study confirmed the suitability of wastewater testing for studying SARS-CoV-2 circulation and established a simplified workflow for the prompt detection and characterization of the virus.

Gene Expression Profiles of the Immuno-Transcriptome in Equine Asthma.

BACKGROUND: Mild equine asthma (MEA) and severe equine asthma (SEA) are two of the most frequent equine airway inflammatory diseases, but knowledge about their pathogenesis is limited. The goal of this study was to investigate gene expression differences in the respiratory tract of MEA- and SEA-affected horses and their relationship with clinical signs. METHODS: Clinical examination and endoscopy were performed in 8 SEA- and 10 MEA-affected horses and 7 healthy controls. Cytological and microbiological analyses of bronchoalveolar lavage (BAL) fluid were performed. Gene expression profiling of BAL fluid was performed by means of a custom oligo-DNA microarray. RESULTS: In both MEA and SEA, genes involved in the genesis, length, and motility of respiratory epithelium cilia were downregulated. In MEA, a significant overexpression for genes encoding inflammatory mediators was observed. In SEA, transcripts involved in bronchoconstriction, apoptosis, and hypoxia pathways were significantly upregulated, while genes involved in the formation of the protective muco-protein film were underexpressed. The SEA group also showed enrichment of gene networks activated during human asthma. CONCLUSIONS: The present study provides new insight into equine asthma pathogenesis, representing the first step in transcriptomic analysis to improve diagnostic and therapeutic approaches for this respiratory disease.

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.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes.

BACKGROUND: The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations. METHODS: Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time. RESULTS: Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe. CONCLUSIONS: Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

New compounds, old problems. The case of C6O4 - a substitute of PFOA - and its effects to the clam Ruditapes philippinarum.

C6O4 (difluoro{[2,2,4,5-tetrafluoro-5-(trifluoromethoxy)-1,3-dioxolan-4-yl]oxy}acetic acid) is a new surfactant and emulsifier used as a substitute of perfluorooctanoic acid (PFOA). Recently, C6O4 has been detected in aquatic environments, but, at present, no information concerning the effects of C6O4 on aquatic species, such as bivalves, are available in the literature. Therefore, in this study we evaluated for the first time the effects of C6O4 (0.1 and 1 µg/L) and PFOA (1 µg/L) to the clam Ruditapes philippinarum. Short-term (7 days) and long-term (21 days) exposures of clams to the two compounds were carried out and numerous biomarkers were measured in haemocytes/haemolymph, as well as in gills and digestive gland. The MANOVA analysis demonstrated statistically significant effects of the independent variables "treatment", "time" and "treatment-time interaction" on the whole dataset of biomarker responses. The two-way ANOVA analysis performed for each biomarker response indicated that the two compounds affected most of the cellular and tissue parameters measured. Despite preliminary, the results obtained suggested that C6O4 - similarly to PFOA - can affect both cellular and biochemical parameters of clams.

Resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax L.): heritability and relationships with body weight, cortisol concentration, and antibody titer.

BACKGROUND: Susceptibility of European sea bass (Dicentrarchus labrax L.) to viral nervous necrosis (VNN) is well-known. Interest towards selective breeding as a tool to enhance genetic resistance in this species has increased sharply due to the major threat represented by VNN for farmed sea bass and limitations concerning specific therapeutical measures. A sea bass experimental population (N = 650) was challenged with nervous necrosis virus (NNV) to investigate genetic variation in VNN mortality. In addition, relationships of this trait with serum cortisol concentration after stress exposure, antibody titer against NNV antigens, and body weight at a fixed age were studied to identify potential indicator traits of VNN resistance. RESULTS: The estimate of heritability for VNN mortality was moderate and ranged from 0.15 (HPD95%, 95% highest posterior density interval: 0.02, 0.31) to 0.23 (HPD95%: 0.06, 0.47). Heritability estimates for cortisol concentration, antibody titer, and body weight were 0.19 (HPD95%: 0.07, 0.34), 0.36 (HPD95%: 0.16, 0.59) and 0.57 (HPD95%: 0.33, 0.84), respectively. Phenotypic relationships between traits were trivial and not statistically significant, except for the estimated correlation between antibody titer and body weight (0.24). Genetic correlations of mortality with body weight or antibody titer (- 0.39) exhibited a 0.89 probability of being negative. A negligible genetic correlation between mortality and cortisol concentration was detected. Antibody titer was estimated to be positively correlated with body weight (0.49). CONCLUSIONS: Antibody titer against NNV offers the opportunity to use indirect selection to enhance resistance, while the use of cortisol concentration as an indicator trait in breeding programs for VNN resistance is questionable. The estimate of heritability for VNN mortality indicates the feasibility of selective breeding to enhance resistance to NNV and raises attention to the development of genomic prediction tools to simplify testing procedures for selection candidates.

The new PFAS C6O4 and its effects on marine invertebrates: First evidence of transcriptional and microbiota changes in the Manila clam Ruditapes philippinarum.

There is growing concern for the wide use ofperfluorooctanoic acid (PFOA) because of its toxic effects on the environment and on human health. A new compound - the so called C6O4 (perfluoro ([5-methoxy-1,3-dioxolan-4-yl]oxy) acetic acid) - was recently introduced as one of the alternative to traditional PFOA, however this was done without any scientific evidence of the effects of C6O4 when dispersed into the environment. Recently, the Regional Agency for the Protection of the Environment of Veneto (Italy) detected high levels of C6O4 in groundwater and in the Po river, increasing the alarm for the potential effects of this chemical into the natural environment. The present study investigates for the first time the effects of C6O4 on the Manila clam Ruditapes philippinarum exposed to environmental realistic concentrations of C6O4 (0.1 µg/L and 1 µg/L) for 7 and 21 days. Furthermore, in order to better understand if C6O4 is a valid and less hazardous alternative to its substitute, microbial and transcriptomic alterations were also investigated in clams exposed to 1 µg/L ofPFOA. Results indicate that C6O4 may cause significant perturbations to the digestive gland microbiota, likely determining the impairment of host physiological homeostasis. Despite chemical analyses suggest a 5 times lower accumulation potential of C604 as compared to PFOA in clam soft tissues, transcriptional analyses reveal several alterations of gene expression profile. A large part of the altered pathways, including immune response, apoptosis regulation, nervous system development, lipid metabolism and cell membrane is the same in C6O4 and PFOA exposed clams. In addition, clams exposed to C6O4 showed dose-dependent responses as well as possible narcotic or neurotoxic effects and reduced activation of genes involved in xenobiotic metabolism. Overall, the present study suggests that the potential risks for marine organism following environmental contamination are not reduced by replacing PFOA with C6O4. In addition, the detection of both C6O4 and PFOA into tissues of clams inhabiting the Lagoon of Venice - where there are no point sources of either compounds - recommends a similar capacity to spread throughout the environment. These results prompt the urgent need to re-evaluate the use of C6O4 as it may represent not only an environmental hazard but also a potential risk for human health.

Genome-wide analysis clarifies the population genetic structure of wild gilthead sea bream (Sparus aurata).

Gilthead sea bream is an important target for both recreational and commercial fishing in Europe, where it is also one of the most important cultured fish. Its distribution ranges from the Mediterranean to the African and European coasts of the North-East Atlantic. Until now, the population genetic structure of this species in the wild has largely been studied using microsatellite DNA markers, with minimal genetic differentiation being detected. In this geographically widespread study, 958 wild gilthead sea bream from 23 locations within the Mediterranean Sea and Atlantic Ocean were genotyped at 1159 genome-wide SNP markers by RAD sequencing. Outlier analyses identified 18 loci potentially under selection. Neutral marker analyses identified weak subdivision into three genetic clusters: Atlantic, West, and East Mediterranean. The latter group could be further subdivided into an Ionian/Adriatic and an Aegean group using the outlier markers alone. Seascape analysis suggested that this differentiation was mainly due to difference in salinity, this being also supported by preliminary genomic functional analysis. These results are of fundamental importance for the development of proper management of this species in the wild and are a first step toward the study of the potential genetic impact of the sea bream aquaculture industry.

fshr: a fish sex-determining locus shows variable incomplete penetrance across flathead grey mullet populations.

Whole-genome sequencing data were produced from a single flathead grey mullet female and assembled into a draft genome sequence, whereas publicly available sequence data were used to obtain a male draft sequence. Two pools, each consisting of 60 unrelated individuals, respectively, of male and female fish were analyzed using Pool-Sequencing. Mapping and analysis of Pool-Seq data against the draft genome(s) revealed >30 loci potentially associated with sex, the most promising locus of which, encoding the follicle-stimulating hormone receptor (fshr) and harboring two missense variants, was genotyped on 245 fish from four Mediterranean populations. Genotype data showed that fshr represents a previously unknown sex-determining locus, although the incomplete association pattern between fshr genotype and sex-phenotype, the variability of such pattern across different populations, and the presence of other candidate loci reveal that a greater complexity underlies sex determination in the flathead grey mullet.

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.