The role of animal hosts in shaping gut microbiome variation.

Millions of years of co-evolution between animals and their associated microbial communities have shaped and diversified the nature of their relationship. Studies continue to reveal new layers of complexity in host-microbe interactions, the fate of which depends on a variety of different factors, ranging from neutral processes and environmental factors to local dynamics. Research is increasingly integrating ecosystem-based approaches, metagenomics and mathematical modelling to disentangle the individual contribution of ecological factors to microbiome evolution. Within this framework, host factors are known to be among the dominant drivers of microbiome composition in different animal species. However, the extent to which they shape microbiome assembly and evolution remains unclear. In this review, we summarize our understanding of how host factors drive microbial communities and how these dynamics are conserved and vary across taxa. We conclude by outlining key avenues for research and highlight the need for implementation of and key modifications to existing theory to fully capture the dynamics of host-associated microbiomes. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

A multibiomarker approach in clams (Ruditapes philippinarum) for a toxicological evaluation of dredged sediments.

The Lagoon of Venice is often dredged for channel maintenance. To avoid harmful consequences to the ecosystem, a proper disposal of bottom sediments requires a preliminary evaluation of its potential toxicity before excavation. Here we evaluated the effects of polluted sediments on clams (Ruditapes philippinarum) using a multibiomarker approach. Bivalves were exposed for 3 and 14 days to five sediment samples collected along a navigation canal between Venice historical centre and the industrial area of Porto Marghera. Immunological, antioxidant, detoxification, and neurotoxicity biomarkers were analysed in haemolymph, gill, and digestive gland. As a control, sediment collected far from pollution sources was used. Two experiments were performed to assess potential seasonal/gametogenic influence in clam sensitivity. A different response of clam biomarkers was observed during the two experiments and among sampling sites. Clams' digestive gland resulted to be the most sensitive tissue analysed showing significant differences among sites in all biomarkers analysed. Greater differences were present due to seasonality rather than exposure. The concentrations of metals and organic pollutants increased from the city centre to the industrial area, highlighting the influence that industrial activities had on the lagoon ecosystem. However, bioaccumulation in clams did not follow the same clear pattern, suggesting low bioavailability of compounds due to relatively high organic matter content. Biomarkers modulation was mainly driven by metals, both present in sediments and bioaccumulated. In comparison, effects of organic pollutants on the biomarkers tested were negligible. Other sources of contamination not investigated (e.g. pesticides) were suggested by neurotoxicity biomarkers alteration.

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.

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.

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.

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.

Macrostructural Evolution of the Mitogenome of Butterflies (Lepidoptera, Papilionoidea).

The mitogenome of the species belonging to the Papilionodea (Lepidoptera) is a double stranded circular molecule containing the 37 genes shared by Metazoa. Eight mitochondrial gene orders are known in the Papilionoidea. MIQGO is the plesiomorphic gene order for this superfamily, while other mitochondrial arrangements have a very limited distribution. 2S1GO gene order is an exception and is present in several Lycaenidae and one species of Hesperiidae. We studied the macrostructural changes generating the gene orders of butterflies by analysing a large data set (611 taxa) containing 5 new mitochondrial sequences/assemblies and 87 de novo annotated mitogenomes. Our analysis supports a possible origin of the intergenic spacer trnQ-nad2, characterising MIQGO, from trnM. We showed that the homoplasious gene order IMQGO, shared by butterflies, species of ants, beetles and aphids, evolved through different transformational pathways. We identify a complicated evolutionary scenario for 2S1GO in Lycaenidae, characterised by multiple events of duplication/loss and change in anticodon of trnS1. We show that the gene orders ES1GO and S1NGO originated through a tandem duplication random loss mechanism. We describe two novel gene orders. Ampittia subvittatus (Hesperiidae) exhibits the gene order 2FFGO, characterised by two copies of trnF, one located in the canonical position and a second placed in the opposite strand between trnR and trnN. Bhutanitis thaidina (Papilionidae) exhibits the gene order 4QGO, characterised by the quadruplication of trnQ.

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.

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.

First Evidence of In Vitro Effects of C6O4-A Substitute of PFOA-On Haemocytes of the Clam Ruditapes philippinarum.

Alternative chemicals to per- and poly-fluoroalkyl substances have recently been introduced in various industrial processes. 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 replacement for perfluorooctanoic acid (PFOA). From an ecotoxicological point of view, in vitro assays are useful tools for assessing the negative effects and understanding the mechanisms of action of chemicals at the cellular level. Here, we present the results of an in vitro study in which the effects of C6O4 were evaluated-for the first time-on haemocytes of the clam Ruditapes philippinarum. Cells were exposed to three concentrations of C6O4 (0.05, 0.5, 5 µg/mL) and the effects on haemocyte viability, haemocyte morphology, differential haemocyte count, lysosomal membrane stability, superoxide anion production, acid phosphatase, and ß-glucuronidase activities, as well as on the percentage of micronuclei and chromosomal aberrations were evaluated. The results demonstrated that C6O4 significantly affected haemocyte morphology, lysosomal membrane stability, hydrolytic enzyme activity, and superoxide anion production, and promoted chromosomal aberrations. To the best of our knowledge, this is the first study revealing the in vitro effects of C6O4, a substitute for PFOA, on haemocytes from a bivalve species.

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.

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.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish.

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

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.

Environmental pharmaceuticals and climate change: The case study of carbamazepine in M. galloprovincialis under ocean acidification scenario.

Contaminants of emerging concern and ocean changes are key environmental stressors for marine species with possibly synergistic, but still unexplored, deleterious effects. In the present study the influence of a simulated ocean acidification scenario (pH = 7.6) was investigated on metabolism and sub-lethal effects of carbamazepine, CBZ (1 µg/L), chosen as one of the most widely diffused pharmaceuticals in marine organisms. A multidisciplinary approach was applied on mussels, M. galloprovincialis, integrating measurement of drug bioaccumulation with changes in the whole transcriptome, responsiveness of various biochemical and cellular biomarkers including immunological parameters, lipid and oxidative metabolism, onset of genotoxic effects. Chemical analyses revealed a limited influence of hypercapnia on accumulation and excretion of CBZ, while a complex network of biological responses was observed in gene expression profile and functional changes at cellular level. The modulation of gamma-aminobutyric acid (GABA) pathway suggested similarities with the Mechanism of Action known for vertebrates: immune responses, cellular homeostasis and oxidative system represented the processes targeted by combined stressors. The overall elaboration of results through a quantitative Weight of Evidence model, revealed clearly increased cellular hazard due to interactions of CBZ with acidification compared to single stressors.

Draft genome assembly and transcriptome data of the icefish Chionodraco myersi reveal the key role of mitochondria for a life without hemoglobin at subzero temperatures.

Antarctic fish belonging to Notothenioidei represent an extraordinary example of radiation in the cold. In addition to the absence of hemoglobin, icefish show a number of other striking peculiarities including large-diameter blood vessels, high vascular densities, mitochondria-rich muscle cells, and unusual mitochondrial architecture. In order to investigate the bases of icefish adaptation to the extreme Southern Ocean conditions we sequenced the complete genome of the icefish Chionodraco myersi. Comparative analyses of the icefish genome with those of other teleost species, including two additional white-blooded and five red-blooded notothenioids, provided a new perspective on the evolutionary loss of globin genes. Muscle transcriptome comparative analyses against red-blooded notothenioids as well as temperate fish revealed the peculiar regulation of genes involved in mitochondrial function in icefish. Gene duplication and promoter sequence divergence were identified as genome-wide patterns that likely contributed to the broad transcriptional program underlying the unique features of icefish mitochondria.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Genome of the Komodo dragon reveals adaptations in the cardiovascular and chemosensory systems of monitor lizards.

Monitor lizards are unique among ectothermic reptiles in that they have high aerobic capacity and distinctive cardiovascular physiology resembling that of endothermic mammals. Here, we sequence the genome of the Komodo dragon Varanus komodoensis, the largest extant monitor lizard, and generate a high-resolution de novo chromosome-assigned genome assembly for V. komodoensis using a hybrid approach of long-range sequencing and single-molecule optical mapping. Comparing the genome of V. komodoensis with those of related species, we find evidence of positive selection in pathways related to energy metabolism, cardiovascular homoeostasis, and haemostasis. We also show species-specific expansions of a chemoreceptor gene family related to pheromone and kairomone sensing in V. komodoensis and other lizard lineages. Together, these evolutionary signatures of adaptation reveal the genetic underpinnings of the unique Komodo dragon sensory and cardiovascular systems, and suggest that selective pressure altered haemostasis genes to help Komodo dragons evade the anticoagulant effects of their own saliva. The Komodo dragon genome is an important resource for understanding the biology of monitor lizards and reptiles worldwide.

Tracing seafood at high spatial resolution using NGS-generated data and machine learning: Comparing microbiome versus SNPs.

Developing reliable tools to trace food origin represents a major goal for producers and control authorities. Here, we test the hypothesis whether NGS-generated data could provide a reliable tool to ensure seafood traceability. As a test case, we used the Manila clam Ruditapes philippinarum, a bivalve mollusk of high commercial interest with worldwide distribution, collected in the Venice lagoon sites subjected to prohibition of clam harvesting because of chemical contamination as well as in authorized clam harvesting areas. The results obtained demonstrated that the geographic origin of Manila clam may be more accurately determined basing on microbiome data than single nucleotide polymorphisms. In particular, combining microbiome data with machine-learning techniques, we provide the experimental evidence that it is possible to trace the clam place of origin at high spatial resolution. Considering its low cost and portability, NGS-analysis of microbiome data might represent a cost-effective, high-resolution tool for reliable food traceability.