Transcriptomic insights into the dominance of two phototrophs throughout the water column of a tropical hypersaline-alkaline crater lake (Dziani Dzaha, Mayotte).

Saline-alkaline lakes often shelter high biomasses despite challenging conditions, owing to the occurrence of highly adapted phototrophs. Dziani Dzaha (Mayotte) is one such lake characterized by the stable co-dominance of the cyanobacterium Limnospira platensis and the picoeukaryote Picocystis salinarum throughout its water column. Despite light penetrating only into the uppermost meter, the prevailing co-dominance of these species persists even in light- and oxygen-deprived zones. Here, a depth profile of phototrophs metatranscriptomes, annotated using genomic data from isolated strains, is employed to identify expression patterns of genes related to carbon processing pathways including photosynthesis, transporters and fermentation. The findings indicate a prominence of gene expression associated with photosynthesis, with a peak of expression around 1 m below the surface, although the light intensity is very low and only red and dark red wavelengths can reach it, given the very high turbidity linked to the high biomass of L. platensis. Experiments on strains confirmed that both species do grow under these wavelengths, at rates comparable to those obtained under white light. A decrease in the expression of photosynthesis-related genes was observed in L. platensis with increasing depth, whereas P. salinarum maintained a very high pool of psbA transcripts down to the deepest point as a possible adaptation against photodamage, in the absence and/or very low levels of expression of genes involved in protection. In the aphotic/anoxic zone, expression of genes involved in fermentation pathways suggests active metabolism of reserve or available dissolved carbon compounds. Overall, L. platensis seems to be adapted to the uppermost water layer, where it is probably maintained thanks to gas vesicles, as evidenced by high expression of the gvpA gene. In contrast, P. salinarum occurs at similar densities throughout the water column, with a peak in abundance and gene expression levels which suggests a better adaptation to lower light intensities. These slight differences may contribute to limited inter-specific competition, favoring stable co-dominance of these two phototrophs.

Environmental Transmission of Symbionts in the Mangrove Crabs Aratus pisonii and Minuca rapax: Acquisition of the Bacterial Community through Larval Development to Juvenile Stage.

Aratus pisonii and Minuca rapax are two brachyuran crabs living with bacterial ectosymbionts located on gill lamellae. One previous study has shown that several rod-shaped bacterial morphotypes are present and the community is dominated by Alphaproteobacteria and Bacteroidota. This study aims to identify the mode of transmission of the symbionts to the new host generations and to identify the bacterial community colonizing the gills of juveniles. We tested for the presence of bacteria using PCR with universal primers targeting the 16S rRNA encoding gene from gonads, eggs, and different larval stages either obtained in laboratory conditions or from the field. The presence of bacteria on juvenile gills was also characterized by scanning electron microscopy, and subsequently identified by metabarcoding analysis. Gonads, eggs, and larvae were negative to PCR tests, suggesting that bacteria are not present at these stages in significant densities. On the other hand, juveniles of both species display three rod-shaped bacterial morphotypes on gill lamellae, and sequencing revealed that the community is dominated by Bacteroidota and Alphaproteobacteria on A. pisonii juveniles, and by Alphaprotobacteria, Bacteroidota, and Acidimicrobia on M. rapax juveniles. Despite the fact that juveniles of both species co-occur in the same biotope, no shared bacterial phylotype was identified. However, some of the most abundant bacteria present in adults are also present in juveniles of the same species, suggesting that juvenile-associated communities resemble those of adults. Because some of these bacteria were also found in crab burrow water, we hypothesize that the bacterial community is established gradually during the life of the crab starting from the megalopa stage and involves epibiosis-competent bacteria that occur in the environment.

Genomes of nine biofilm-forming filamentous strains of Cyanobacteria (genera Jaaginema, Scytonema, and Karukerafilum gen. nov.) isolated from mangrove habitats of Guadeloupe (Lesser Antilles).

Biofilm-forming cyanobacteria are abundant in mangrove ecosystems, colonizing various niches including sediment surface and periphyton where they can cover large areas, yet have received limited attention. Several filamentous isolates were recently isolated from Guadeloupe, illustrating the diversity and novelty present in these biofilms. In this study, nine strains belonging to three novel lineages found abundantly in Guadeloupe biofilms are characterized by genome sequencing, morphological and ultrastructural examination, metabolome fingerprinting and searched for secondary metabolites biosynthesis pathways. Assignation of two lineages to known genera is confirmed, namely Scytonema and Jaaginema. The third lineage corresponds to a new Coleofasciculales genus herein described as Karukerafilum gen. nov. The four strains belonging to this genus group into two subclades, one of which displays genes necessary for nitrogen fixation as well as the complete pathway for geosmin production. This study gives new insights into the diversity of mangrove biofilm-forming cyanobacteria, including genome-based description of a new genus and the first genome sequence available for the genus Jaaginema.

Multi-omics analyses from a single sample: prior metabolite extraction does not alter the 16S rRNA-based characterization of prokaryotic community in a diversity of sample types.

Massive sequencing of the 16S rRNA gene has become a standard first step to describe and compare microbial communities from various samples. Parallel analysis of high numbers of samples makes it relevant to the statistical testing of the influence of natural or experimental factors and variables. However, these descriptions fail to document changes in community or ecosystem functioning. Nontargeted metabolomics are a suitable tool to bridge this gap, yet extraction protocols are different. In this study, prokaryotic community compositions are documented by 16S rRNA gene sequencing after direct DNA extraction or after metabolites extraction followed by DNA extraction. Results obtained using the V3-V4 region on nonaxenic cultures of cyanobacteria, lake water column, biofilm, and gut of wild and lab-reared fish indicate that prior extraction of metabolites does not influence the obtained image of prokaryotic communities. This validates sequential extraction of metabolites followed by DNA as a way to combine 16S rRNA sequencing with metabolome characterization from a single sample. This approach has the potential to complement community structure characterization with a proxy of their functioning, without the uncertainties associated with the use of separate samples.

Phytoplanktonic species in the haloalkaline Lake Dziani Dzaha select their archaeal microbiome.

Microorganisms are key contributors of aquatic biogeochemical cycles but their microscale ecology remains largely unexplored, especially interactions occurring between phytoplankton and microorganisms in the phycosphere, that is the region immediately surrounding phytoplankton cells. The current study aimed to provide evidence of the phycosphere taking advantage of a unique hypersaline, hyperalkaline ecosystem, Lake Dziani Dzaha (Mayotte), where two phytoplanktonic species permanently co-dominate: a cyanobacterium, Arthrospira fusiformis, and a green microalga, Picocystis salinarum. To assay phycospheric microbial diversity from in situ sampling, we set up a flow cytometry cell-sorting methodology for both phytoplanktonic populations, coupled with metabarcoding and comparative microbiome diversity. We focused on archaeal communities as they represent a non-negligible part of the phycospheric diversity, however their role is poorly understood. This work is the first which successfully explores in situ archaeal diversity distribution showing contrasted phycospheric compositions, with P. salinarum phycosphere notably enriched in Woesearchaeales OTUs while A. fusiformis phycosphere was enriched in methanogenic lineages affiliated OTUs such as Methanomicrobiales or Methanofastidiosales. Most archaeal OTUs, including Woesearchaeales considered in literature as symbionts, were either ubiquitous or specific of the free-living microbiome (i.e. present in the 3-0.2 µm fraction). Seminally, several archaeal OTUs were enriched from the free-living microbiome to the phytoplankton phycospheres, suggesting (i) either the inhibition or decrease of other OTUs, or (ii) the selection of specific OTUs resulting from the physical influence of phytoplanktonic species on surrounding Archaea.

Intra-population genomic diversity of the bloom-forming cyanobacterium, Aphanizomenon gracile, at low spatial scale.

Cyanobacteria are oxygenic photosynthetic bacteria that perform a substantial part of the global primary production. Some species are responsible for catastrophic environmental events, called blooms, which have become increasingly common in lakes and freshwater bodies as a consequence of global changes. Genotypic diversity is considered essential for marine cyanobacterial population, allowing it to cope with spatio-temporal environmental variations and to adapt to specific micro-niches in the ecosystem. This aspect is underestimated in the study of bloom development, however, and given little notice in studies of the ecology of harmful cyanobacteria. Here we compared the genomes of four strains of Aphanizomenon gracile, a species of filamentous toxinogenic cyanobacteria (Nostocales) found worldwide in fresh and brackish water. Millimeter-sized fascicles were isolated from a single water sample and have been maintained in culture since 2010. A comparative study revealed extensive heterogeneity in gene contents, despite similar genome size and high similarity indices. These variations were mainly associated with mobile genetic elements and biosynthetic gene clusters. For some of the latter, metabolomic analysis confirmed the production of related secondary metabolites, such as cyanotoxins and carotenoids, which are thought to play a fundamental role in the cyanobacterial fitness. Altogether, these results demonstrated that an A. gracile bloom could be a highly diverse population at low spatial scale and raised questions about potential exchanges of essential metabolites between individuals.

Bacterial Microbiota of Ostreobium, the Coral-Isolated Chlorophyte Ectosymbiont, at Contrasted Salinities.

Microscopic filaments of the siphonous green algae Ostreobium (Ulvophyceae, Bryopsidales) colonize and dissolve the calcium carbonate skeletons of coral colonies in reefs of contrasted salinities. Here, we analyzed their bacterial community's composition and plasticity in response to salinity. Multiple cultures of Pocillopora coral-isolated Ostreobium strains from two distinct rbcL lineages representative of IndoPacific environmental phylotypes were pre-acclimatized (>9 months) to three ecologically relevant reef salinities: 32.9, 35.1, and 40.2 psu. Bacterial phylotypes were visualized for the first time at filament scale by CARD-FISH in algal tissue sections, within siphons, at their surface or in their mucilage. Ostreobium-associated microbiota, characterized by bacterial 16S rDNA metabarcoding of cultured thalli and their corresponding supernatants, were structured by host genotype (Ostreobium strain lineage), with dominant Kiloniellaceae or Rhodospirillaceae (Alphaproteobacteria, Rhodospirillales) depending on Ostreobium lineage, and shifted Rhizobiales' abundances in response to the salinity increase. A small core microbiota composed of seven ASVs (~1.5% of thalli ASVs, 19-36% cumulated proportions) was persistent across three salinities in both genotypes, with putative intracellular Amoebophilaceae and Rickettsiales_AB1, as well as Hyphomonadaceae and Rhodospirillaceae also detected within environmental (Ostreobium-colonized) Pocillopora coral skeletons. This novel knowledge on the taxonomic diversity of Ostreobium bacteria paves the way to functional interaction studies within the coral holobiont.

Disruption of fish gut microbiota composition and holobiont's metabolome during a simulated Microcystis aeruginosa (Cyanobacteria) bloom.

BACKGROUND: Cyanobacterial blooms are one of the most common stressors encountered by metazoans living in freshwater lentic systems such as lakes and ponds. Blooms reportedly impair fish health, notably through oxygen depletion and production of bioactive compounds including cyanotoxins. However, in the times of the "microbiome revolution", it is surprising that so little is still known regarding the influence of blooms on fish microbiota. In this study, an experimental approach is used to demonstrate that blooms affect fish microbiome composition and functions, as well as the metabolome of holobionts. To this end, the model teleost Oryzias latipes is exposed to simulated Microcystis aeruginosa blooms of various intensities in a microcosm setting, and the response of bacterial gut communities is evaluated in terms of composition and metabolome profiling. Metagenome-encoded functions are compared after 28 days between control individuals and those exposed to highest bloom level. RESULTS: The gut bacterial community of O. latipes exhibits marked responses to the presence of M. aeruginosa blooms in a dose-dependent manner. Notably, abundant gut-associated Firmicutes almost disappear, while potential opportunists increase. The holobiont's gut metabolome displays major changes, while functions encoded in the metagenome of bacterial partners are more marginally affected. Bacterial communities tend to return to original composition after the end of the bloom and remain sensitive in case of a second bloom, reflecting a highly reactive gut community. CONCLUSION: Gut-associated bacterial communities and holobiont functioning are affected by both short and long exposure to M. aeruginosa, and show evidence of post-bloom resilience. These findings point to the significance of bloom events to fish health and fitness, including survival and reproduction, through microbiome-related effects. In the context of increasingly frequent and intense blooms worldwide, potential outcomes relevant to conservation biology as well as aquaculture warrant further investigation. Video Abstract.

Metagenome-Based Exploration of Bacterial Communities Associated with Cyanobacteria Strains Isolated from Thermal Muds.

Cyanobacteria constitute a pioneer colonizer of specific environments for whom settlement in new biotopes precedes the establishment of composite microbial consortia. Some heterotrophic bacteria constitute cyanobacterial partners that are considered as their cyanosphere, being potentially involved in mutualistic relationships through the exchange and recycling of key nutrients and the sharing of common goods. Several non-axenic cyanobacterial strains have been recently isolated, along with their associated cyanospheres, from the thermal mud of Balaruc-les-Bains (France) and the biofilms of the retention basin where they develop. The community structure and relationships among the members of the isolated cyanobacterial strains were characterized using a metagenomic approach combined with taxonomic and microscopic descriptions of the microbial consortia. The results provided insights into the potential role and metabolic capabilities of the microorganisms of thermal mud-associated cyanobacterial biofilms. Thus, the physical proximity, host-specificity, and genetic potential functions advocate for their complementarity between cyanobacteria and their associated microbiota. Besides these findings, our results also highlighted the great influence of the reference protein database chosen for performing functional annotation of the metagenomes from organisms of the cyanosphere and the difficulty of selecting one unique database that appropriately covers both autotroph and heterotroph metabolic specificities.

Computational Metabolomics Tools Reveal Subarmigerides, Unprecedented Linear Peptides from the Marine Sponge Holobiont Callyspongia subarmigera.

A detailed examination of a unique molecular family, restricted to the Callyspongia genus, in a molecular network obtained from an in-house Haplosclerida marine sponge collection (including Haliclona, Callyspongia, Xestospongia, and Petrosia species) led to the discovery of subarmigerides, a series of rare linear peptides from Callyspongia subarmigera, a genus mainly known for polyacetylenes and lipids. The structure of the sole isolated peptide, subarmigeride A (1) was elucidated through extensive 1D and 2D NMR spectroscopy, HRMS/MS, and Marfey's method to assign its absolute configuration. The putative structures of seven additional linear peptides were proposed by an analysis of their respective MS/MS spectra and a comparison of their fragmentation patterns with the heptapeptide 1. Surprisingly, several structurally related analogues of subarmigeride A (1) occurred in one distinct cluster from the molecular network of the cyanobacteria strains of the Guadeloupe mangroves, suggesting that the true producer of this peptide family might be the microbial sponge-associated community, i.e., the sponge-associated cyanobacteria.

Establishment of the Bacterial Microbiota in a Lab-Reared Model Teleost Fish, the Medaka Oryzias latipes.

Oryzias latipes is an important model organism for physiology, genetics, and developmental studies, and has also emerged as a relevant vertebrate model for aquatic ecotoxicology. Knowledge regarding its associated microbiota on the other hand is still scarce and limited to adults, despite the relevance of the associated microbiome to the host's biology. This study provides the first insights into the establishment of bacterial microbiota during early developmental stages of laboratory-reared medaka using a 16S-rRNA-sequencing-based approach. Major shifts in community compositions are observed, from a Proteobacteria-dominated community in larvae and juveniles to a more phylum-diverse community towards adulthood, with no obvious difference between female and male specimens. Major bacterial taxa found in adults, including genera Cetobacterium and ZOR0006, establish progressively and are rare during early stages. Dominance shifts are comparable to those documented in another major model teleost, the zebrafish. Results from this study provide a basis for future work investigating the influence of medaka-associated bacteria during host development.

Fish gut-associated bacterial communities in a tropical lagoon (Aghien lagoon, Ivory Coast).

Aghien lagoon (Ivory Coast) is a eutrophic freshwater lagoon that harbors high biomasses of phytoplankton. Despite Increasing interest in fish gut microbiomes diversity and functions, little data is currently available regarding wild species from tropical west African lakes. Here, gut-associated bacterial communities are investigated in four fish species that are consumed by locale populations, namely the Cichlidae Hemichromis fasciatus, Tilapia guineensis and Sarotherodon melanotheron, and the Claroteidae Chrysichthys nigrodigitatus. Species-related differences are identified, that can be attributed to host phylogeny and diet. Important variations throughout the year are observed in T. guineensis and C. nigrodigitatus. This result emphasized the importance of time-series sampling and comparison with environmental variables even in tropical regions, that are not often conducted in wild populations. Effects of environmental factors (anthropogenic or not) on the microbiota and potential outcomes for fish health and populations sustainability need to be further explored. Interestingly, fish appear as major reservoirs of bacterial diversity, suggesting that they could contribute to the overall stability and resilience of bacterial communities present in the Aghien lagoon.

Gut microbiota and holobiont metabolome composition of the medaka fish (Oryzias latipes) are affected by a short exposure to the cyanobacterium Microcystis aeruginosa.

Blooms of toxic cyanobacteria are a common stress encountered by aquatic fauna. Evidence indicates that long-lasting blooms affect fauna-associated microbiota. Because of their multiple roles, host-associated microbes are nowadays considered relevant to ecotoxicology, yet the respective timing of microbiota versus functional changes in holobionts response needs to be clarified. The response of gut microbiota and holobiont's metabolome to exposure to a dense culture of Microcystis aeruginosa was investigated as a microcosm-simulated bloom in the model fish species Oryzias latipes (medaka). Both gut microbiota and gut metabolome displayed significant composition changes after only 2 days of exposure. A dominant symbiont, member of the Firmicutes, plummeted whereas various genera of Proteobacteria and Actinobacteriota increased in relative abundance. Changes in microbiota composition occurred earlier and faster compared to metabolome composition. Liver and muscle metabolome were much less affected than guts, supporting that the gut and associated microbiota are in the front row upon exposure. This study highlights that even short cyanobacterial blooms, that are increasingly frequent, trigger changes in microbiota composition and holobiont metabolome. It emphasizes the relevance of multi-omics approaches to explore organism's response to an ecotoxicological stress.

Investigation of Capitella spp. symbionts in the context of varying anthropic pressures: First occurrence of a transient advantageous epibiosis with the giant bacteria Thiomargarita sp. to survive seasonal increases of sulfides in sediments.

Capitella spp. is considered as an important ecological indicator of eutrophication due to its high densities in organic-rich, reduced, and sometimes polluted coastal ecosystems. We investigated whether such ability to cope with adverse ecological contexts might be a response to the microorganisms these worms are associated with. In populations from the French Atlantic coast (Roscoff, Brittany), we observed an epibiotic association covering the tegument of 20-30% specimens from an anthropized site while individuals from a reference, non-anthropized site were devoid of any visible epibionts. Using RNAseq, molecular and microscopic analyses, we described and compared the microbial communities associated with the epibiotic versus the non-epibiotic specimens at both locations. Interestingly, data showed that the epibiosis is characterized by sulfur-oxidizing bacteria among which the giant bacterium Thiomargarita sp., to date only described in deep sea habitats. Survey of Capitella combined with the geochemical analysis of their sediment revealed that epibiotic specimens are always found in muds with the highest concentration of sulfides, mostly during the summer. Concomitantly, tolerance tests demonstrated that the acquisition of epibionts increased survival against toxic level of sulfides. Overall, the present data highlight for the first time a peculiar plastic adaptation to seasonal variations of the habitat based on a transcient epibiosis allowing a coastal species to survive temporary harsher conditions.

Endozoicomonadaceae symbiont in gills of Acesta clam encodes genes for essential nutrients and polysaccharide degradation.

Gammaproteobacteria from the family Endozoicomonadaceae have emerged as widespread associates of dense marine animal communities. Their abundance in coral reefs involves symbiotic relationships and possibly host nutrition. We explored functions encoded in the genome of an uncultured Endozoicomonadaceae 'Candidatus Acestibacter aggregatus' that lives inside gill cells of large Acesta excavata clams in deep-water coral reefs off mid-Norway. The dominance and deep branching lineage of this symbiont was confirmed using 16S rRNA gene sequencing and phylogenomic analysis from shotgun sequencing data. The 4.5 Mb genome binned in this study has a low GC content of 35% and is enriched in transposon and chaperone gene annotations indicating ongoing adaptation. Genes encoding functions potentially involved with the symbiosis include ankyrins, repeat in toxins, secretion and nutritional systems. Complete pathways were identified for the synthesis of eleven amino acids and six B-vitamins. A minimal chitinolytic machinery was indicated from a glycosyl hydrolase GH18 and a lytic polysaccharide monooxygenase LPMO10. Expression of the latter was confirmed using proteomics. Signal peptides for secretion were identified for six polysaccharide degrading enzymes, ten proteases and three lipases. Our results suggest a nutritional symbiosis fuelled by enzymatic products from extracellular degradation processes.

Draft Genome Sequence of the Toxic Freshwater Microcystis aeruginosa Strain PMC 728.11 (Cyanobacteria, Chroococcales).

Microcystis aeruginosa is one of the major species that cause toxic cyanobacterial blooms in freshwater systems worldwide. Here, we report the draft genome sequence of M. aeruginosa PMC 728.11, a microcystin-producing cyanobacterium isolated from the freshwater reservoir of Juanon in Valence, France. The genome sequence contains 276 contigs, consisting of 5,536,025 bp and 5,594 putative protein-coding genes, among which are several biosynthetic gene clusters encoding enzyme complexes involved in the production of various bioactive and toxic metabolites.

Microbiome-Aware Ecotoxicology of Organisms: Relevance, Pitfalls, and Challenges.

Over the last 15 years, the advent of high-throughput "omics" techniques has revealed the multiple roles and interactions occurring among hosts, their microbial partners and their environment. This microbiome revolution has radically changed our views of biology, evolution, and individuality. Sitting at the interface between a host and its environment, the microbiome is a relevant yet understudied compartment for ecotoxicology research. Various recent works confirm that the microbiome reacts to and interacts with contaminants, with consequences for hosts and ecosystems. In this paper, we thus advocate for the development of a "microbiome-aware ecotoxicology" of organisms. We emphasize its relevance and discuss important conceptual and technical pitfalls associated with study design and interpretation. We identify topics such as functionality, quantification, temporality, resilience, interactions, and prediction as major challenges and promising venues for microbiome research applied to ecotoxicology.

New Benthic Cyanobacteria from Guadeloupe Mangroves as Producers of Antimicrobials.

Benthic cyanobacteria strains from Guadeloupe have been investigated for the first time by combining phylogenetic, chemical and biological studies in order to better understand the taxonomic and chemical diversity as well as the biological activities of these cyanobacteria through the effect of their specialized metabolites. Therefore, in addition to the construction of the phylogenetic tree, indicating the presence of 12 potentially new species, an LC-MS/MS data analysis workflow was applied to provide an overview on chemical diversity of 20 cyanobacterial extracts, which was linked to antimicrobial activities evaluation against human pathogenic and ichtyopathogenic environmental strains.

Low-diversity bacterial microbiota in Southern Ocean representatives of lanternfish genera Electrona, Protomyctophum and Gymnoscopelus (family Myctophidae).

Myctophids are among the most abundant mesopelagic teleost fishes worldwide. They are dominant in the Southern Ocean, an extreme environment where they are important both as consumers of zooplankton as well as food items for larger predators. Various studies have investigated myctophids diet, but no data is yet available regarding their associated microbiota, despite that the significance of bacterial communities to fish health and adaptation is increasingly acknowledged. In order to document microbiota in key fish groups from the Southern Ocean, the bacterial communities associated with the gut, fin, gills and light organs of members of six species within the three myctophid genera Electrona, Protomyctophum and Gymnoscopelus were characterized using a 16S rRNA-based metabarcoding approach. Gut communities display limited diversity of mostly fish-specific lineages likely involved in food processing. Fin and skin communities display diversity levels and compositions resembling more those found in surrounding seawater. Community compositions are similar between genera Electrona and Protomyctophum, that differ from those found in Gymnoscopelus and in water. Low abundances of potentially light-emitting bacteria in light organs support the hypothesis of host production of light. This first description of myctophid-associated microbiota, and among the first on fish from the Southern Ocean, emphasizes the need to extend microbiome research beyond economically-important species, and start addressing ecologically-relevant species.

Characterization of ecto- and endoparasite communities of wild Mediterranean teleosts by a metabarcoding approach.

Next-generation sequencing methods are increasingly used to identify eukaryotic, unicellular and multicellular symbiont communities within hosts. In this study, we analyzed the non-specific reads obtained during a metabarcoding survey of the bacterial communities associated to three different tissues collected from 13 wild Mediterranean teleost fish species. In total, 30 eukaryotic genera were identified as putative parasites of teleosts, associated to skin mucus, gills mucus and intestine: 2 ascomycetes, 4 arthropods, 2 cnidarians, 7 nematodes, 10 platyhelminthes, 4 apicomplexans, 1 ciliate as well as one order in dinoflagellates (Syndiniales). These results highlighted that (1) the metabarcoding approach was able to uncover a large spectrum of symbiotic organisms associated to the fish species studied, (2) symbionts not yet identified in several teleost species were putatively present, (3) the parasitic diversity differed markedly across host species and (4) in most cases, the distribution of known parasitic genera within tissues is in accordance with the literature. The current work illustrates the large insights that can be gained by making maximum use of data from a metabarcoding approach.