The Isolation of Vibrio crassostreae and V. cyclitrophicus in Lesser-Spotted Dogfish (Scyliorhinus canicula) Juveniles Reared in a Public Aquarium

The genus Vibrio currently contains 147 recognized species widely distributed, including pathogens for aquatic organisms. Vibrio infections in elasmobranchs are poorly reported, often with identifications as Vibrio sp. and without detailed diagnostic insights. The purpose of this paper is the description of the isolation and identification process of Vibrio spp. following a mortality event of Scyliorhinus canicula juvenile reared in an Italian public aquarium. Following investigations aimed at excluding the presence of different pathogens of marine fish species (parasites, bacteria, Betanodavirus), several colonies were isolated and subjected to species identification using the available diagnostic techniques (a biochemical test, MALDI-TOF MS, and biomolecular analysis). Discrepancies were observed among the methods; the limits of biochemistry as a unique tool for Vibrio species determination were detected through statistical analysis. The use of the rpoB gene, as a diagnostic tool, allowed the identification of the isolates as V. crassostreae and V. cyclotrophicus. Although the pathogenic role of these microorganisms in lesser-spotted dogfish juveniles has not been demonstrated, and the presence of further pathogens cannot be excluded, this study allowed the isolation of two Vibrio species in less-studied aquatic organisms, highlighting the weaknesses and strengths of the different diagnostic methods applied.


Introduction
The bacteria of the Vibrionaceae family are Gram-negative curved rods that globally occur in marine, estuarine and freshwater ecosystems. They occupy habitats ranging from the deep-sea to shallow aquatic environments [1]. Among them, some species are important for natural systems, including carbon cycle and osmoregulation, and as free-living inhabitants in the water column or associated with particulate matter [2]. The genus Vibrio currently contains 147 recognized species widely distributed in aquatic environments [3]. Some vibrios cause water-and seafood-related outbreaks of gastrointestinal infections in humans [4], and many can be pathogenic for marine vertebrates [5,6] and invertebrates [7,8]. Paillard et al. recognized that the emergence of vibrios as etiological agents of diseases is likely to increase over the coming years due to ocean warming [9]. Classical vibriosis is generally characterized by lethargic movement of affected fish [10][11][12], the presence of various other pathogens floating in the water are impacted. Nevertheless, nitrifying bacteria (with key roles in the nitrogen cycle) are not carried to the UV sterilizer and therefore do not come in range and will not be lost. The following water physicochemical parameters were measured by a multiparametric field probe (HACH, HQd Field Case): temperature (17 ± 1 °C); salinity (36 ± 1 PSU); dissolved oxygen (8.5 ± 0.4 mg/L); pH (8.0 ± 0.3). The specimens ( Figure 1) showed the first symptoms a few months after hatching, i.e., anomalies in swimming, fast breathing, and lack of appetite. For the analysis, symptomatic and moribund fish were sent inside a double plastic bag (containing ⅓ water and ⅔ air/oxygen) and sent refrigerated to the fish diseases laboratory of the Istituto Zooprofilattico Sperimentale of Piemonte, Liguria and Valle d'Aosta, Turin, Italy.

Necropsy
In the laboratory, the fish were euthanatized using an overdose of ethyl 3aminobenzoate methanesulfonate (MS-222, Sigma Aldrich, St. Louis, USA) following the current legislation. Before necropsy, biometrical parameters of each animal were recorded (total length (TL) and total weight (TW)). The fish were then macroscopically examined to highlight external alterations. After dissection with sterile tools, the coelomic cavity and the internal organs were subjected to visual inspection for the evaluation of anatomopathological changes.

Parasitological Examination
A general approach was used for parasitological examination in order to detect esoand ectoparasites. Gills, skin, visceral cavity, and digestive tract were analyzed, both macroscopically and microscopically, with the use of a microscope. The microscope used for the analysis is an OLYMPUS BX40, with magnification ranging from 4× to 40×.

Virological Examination
For virological analysis, a sample of a portion of the brain was taken for the search for viral encephalopathy and retinopathy virus (Betanodavirus). The total RNA was extracted using the AllPrep DNA/RNA Micro Kit (QIAGEN, Hilden, Germany) following the manufacturerʹs protocol. Subsequently, cDNA was synthesized from extracted RNA using the QuantiTect Reverse Transcription Kit (QIAGEN, Hilden, Germany). A real-time RT-PCR targeting the viral encephalopathy and retinopathy virus RNA2 portion was used for viral diagnosis, following the protocol described by Panzarin et al. [40]. The reactions were carried out using the BioRad CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, USA), and the data were analyzed by Bio-Rad CFX Maestro Software version 1.1.

Necropsy
In the laboratory, the fish were euthanatized using an overdose of ethyl 3-aminobenzoate methanesulfonate (MS-222, Sigma Aldrich, St. Louis, USA) following the current legislation. Before necropsy, biometrical parameters of each animal were recorded (total length (TL) and total weight (TW)). The fish were then macroscopically examined to highlight external alterations. After dissection with sterile tools, the coelomic cavity and the internal organs were subjected to visual inspection for the evaluation of anatomopathological changes.

Parasitological Examination
A general approach was used for parasitological examination in order to detect esoand ectoparasites. Gills, skin, visceral cavity, and digestive tract were analyzed, both macroscopically and microscopically, with the use of a microscope. The microscope used for the analysis is an OLYMPUS BX40, with magnification ranging from 4× to 40×.

Virological Examination
For virological analysis, a sample of a portion of the brain was taken for the search for viral encephalopathy and retinopathy virus (Betanodavirus). The total RNA was extracted using the AllPrep DNA/RNA Micro Kit (QIAGEN, Hilden, Germany) following the manufacturer's protocol. Subsequently, cDNA was synthesized from extracted RNA using the QuantiTect Reverse Transcription Kit (QIAGEN, Hilden, Germany). A real-time RT-PCR targeting the viral encephalopathy and retinopathy virus RNA2 portion was used for viral diagnosis, following the protocol described by Panzarin et al. [40]. The reactions were carried out using the BioRad CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, USA), and the data were analyzed by Bio-Rad CFX Maestro Software version 1.1.

Bacteriological and Biochemical Analyses
Bacteriological sampling was carried out from the head kidney, brain, and blood. The organs samples were taken using a 1 µL sterile calibrated loop, while blood sampling was performed immediately after the euthanasia procedure, taking it directly from the sharks' hearts with an insulin syringe, inoculated on plates, and spread using a 1 µL sterile loop. The first isolation was performed on Columbia blood agar (BA), tryptic soy agar (TSA) supplemented with 2% NaCl, and Monsur medium (thiosulfate citrate bile sucrose, TCBS); the latter medium was used for its selectivity toward Vibrionaceae. Plates were incubated at 22 ± 2 • C for a total of 72 h, and growth was checked daily. After growth, colonies were cloned on BA for identification analyses. After Gram staining, pure isolates were identified by API ® 20E tests, matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) using VITEK-MS (BioMerieux, Marcy-l'Étoile, France), and molecular analyses. API ® tests were also used for the biochemical characterization of strains. API ® galleries were incubated for 48 h at 22 ± 2 • C; tests were read every 24 h, and reagents were added according to the manufacturer's instructions. The results were issued using the specific API ® test reading software (APIWEB TM ) available on the BioMerieux website (https://apiweb.biomerieux.com, accessed on 30 September 2021).

Biomolecular Analysis for Bacterial Identification
DNA extraction was performed using the boiling and freeze-thawing protocol as described by Pastorino et al. [41]. RNA polymerase beta subunit (rpoB) gene was amplified to discriminate among different Vibrio species. Specific primers for the rpoB gene, developed by Ki et al. [42], and the protocol suggested by the same authors were used. Amplicons were run on 2% GelGreen (Biotium, Landing Pkwy, USA) stained agarose gel and then visualized under UV exposure. A 50-2000 kb ladder (Amplisize Molecular Ruler, Bio-Rad, Hercules, USA) was used as a molecular marker. The amplicons were purified with an ExtractMe DNA Clean-Up and Gel-Out Kit (Blirt, Gdańsk, Poland), according to the manufacturer's instructions. The purified PCR products were bidirectionally sequenced using Big Dye 3.1 (Applied Biosystems, Waltham, USA) chemistry and the same primers used for PCR amplification. Cycle-sequencing products were purified using Dye Ex 2.0 Spin Kit (QIAGEN, Hilden, Germany) and sequenced in an ABI3130xl Genetic analyzer (Applied Biosystems, Waltham, USA). Contig assembly of the DNA sequences was performed using the Lasergene Software package (DNASTAR, Madison, USA). The rpoB sequences were compared with nucleotide sequences in the GenBank database using the Basic Local Alignment Search Tool (BLAST) search algorithm. Subsequently, rpoB sequences were deposited on the GenBank database.

Phylogenetic Analysis
A neighbor-joining analysis was performed [43] using Molecular Evolutionary Genetics Analysis software (MEGAX) [44]. The rpoB gene sequences of 19 isolated strains (sequences with 100% identity with the reference sequence on BLAST were not used in the analysis), 10 rpoB sequences of different Vibrio species, and 4 genomic sequences for Vibrio species for which rpoB gene sequences were not available alone (V. crassostreae, V. cyclitrophicus, V. celticus, and V. gigantis) were considered for the phylogenetic analysis. Evolutionary distances were ascertained via the maximum composite likelihood method [45]. A bootstrap test of 1000 replicates was performed; a cut-off of 50% was set for the computation of the bootstrap condensed tree.

Statistical Analysis
Principal component analysis (PCA) was performed to reduce the dimensionality of the biochemical dataset, increasing interpretability but, at the same time, minimizing information loss. PCA was performed to illustrate the clustering of bacteria species based on biochemical features. This was accomplished considering both phenotypic (API ® 20E) and molecular (PCR) species identification. Statistical analysis was performed using R software (version 1.1.463, RStudio, Inc., Boston, MA, USA).

Necropsy, Parasitological, and Virological Examinations
A total of 20 lesser-spotted dogfish with an average weight of 2.16 ± 0.29 g (minimum weight: 1.70 g; maximum weight: 2.78 g) and an average length of 94.60 ± 4.03 mm (minimum length: 87 mm; maximum length: 103 mm) were analyzed. No external lesions were found, and the necropsy of the sharks did not show macroscopic lesions of the viscera. Parasitological (internal and external examination) and virological (Betanodavirus RT-PCR) tests yielded negative results.

Bacteriological Analysis
Out of 20 fish tested, 12 (12/20; 60%) tested positive from at least one matrix (head kidney, brain, or blood). A first visual examination of the plates was conducted, leading to the identification of 24 morphologically different isolates (Table 1). Each selected colony was subcloned to BA for subsequent tests. Following Gram staining, all isolates were identified as Gram-negative curved rods. Based on morphological and coltural characteristics, the identification by API ® test envisaged the use of the API ® 20E galleries. The results of the biochemical analyses carried out in a micro-method are summarized in Figure 2. On the other hand, identification by VITEK-MS did not produce any valid results according to the instrument cut-offs, so it was not possible to identify the isolates with this method. Table 1. Characteristics of the isolates. The table shows the number of sharks, the letter linked to the different isolated colonies (A, B, C; only in specimens showing more than one colony), the matrix from which the sampling was made, and the first isolation medium. Columbia blood agar (BA); tryptic soy agar supplemented with 2% NaCl (TSA2); thiosulfate citrate bile sucrose (TCBS).

Specimen
Isolates

Biomolecular Analysis for Bacterial Identification
All the 24 strains isolated showed a fragment of 730 bp from the rpoB gene amplification. Of those, 11 isolates (11/24; 45.8%) were identified at species level as Vibrio cyclitrophicus and 13 (13/24; 54.2%) as V. crassostreae (BLASTn nucleotide sequence identity value ranging from 98 to 100%). Results are reported in Table 2.
All obtained sequences, except for sequences with 100% identity value with the reference sequences in BLAST, were deposited on GenBank, a freely available online database (Accession Numbers: OM158211-OM158229).

Biomolecular Analysis for Bacterial Identification
All the 24 strains isolated showed a fragment of 730 bp from the rpoB gene amplification. Of those, 11 isolates (11/24; 45.8%) were identified at species level as Vibrio cyclitrophicus and 13 (13/24; 54.2%) as V. crassostreae (BLASTn nucleotide sequence identity value ranging from 98 to 100%). Results are reported in Table 2.
All obtained sequences, except for sequences with 100% identity value with the reference sequences in BLAST, were deposited on GenBank, a freely available online database (Accession Numbers: OM158211-OM158229).

Phylogenetic Analysis
The bootstrap condensed tree is shown in Figure 3. Two main clusters are present: one containing the 10 rpoB sequences of different Vibrio species subdivided into three different subclusters, the second containing V. cyclitrophicus sequences with a high bootstrap value (97%). V. crassostreae isolates are not shown to be clustering with other Vibrio species, considering a bootstrap value cut-off of 50%.

Phylogenetic Analysis
The bootstrap condensed tree is shown in Figure 3. Two main clusters are present: one containing the 10 rpoB sequences of different Vibrio species subdivided into three different subclusters, the second containing V. cyclitrophicus sequences with a high bootstrap value (97%). V. crassostreae isolates are not shown to be clustering with other Vibrio species, considering a bootstrap value cut-off of 50%. Figure 3. Phylogenetic relationship among Vibrio species isolated from lesser-spotted dogfish and different Vibrio species. Phylogenetic tree was constructed using MEGAX and neighbor-joining method. A bootstrap test of 1000 replicates was performed; the bootstrap condensed tree, using a cut-off value of 50% is shown. The sequences with "G" displayed after the scientific name derived from a deposited genome.

Biochemical Features of Isolates
Morphological and biochemical analysis showed that all isolates were Gram-negative bacteria, positive for cytochrome oxidase, produced indole from tryptophan, and liquefied gelatin from enzyme gelatinase. None of the isolated strains produced lysine decarboxylase, ornithine decarboxylase, hydrogen sulfide, urease, nor were they able to utilize citrate. The enzyme tryptophan deaminase was negative in almost all strains. Generally, the isolates produced acid from glucose, arabinose, amygdalin, and mannitol. No acid was produced from amygdalin and arabinose. It is important to highlight a certain variability in biochemical features for the following test: decarboxylation of the amino acid arginine by arginine dihydrolase, detection of acetoin (acetyl methylcarbinol) produced by fermentation of glucose by bacteria utilizing the butylene glycol pathway (Voges-Proskauer test), and fermentation of sucrose and melibiose ( Figure 2). The first PCA (Figure 4) based on phenotypic identification (API ® 20E) showed that the first (PC1) and second (PC2) components accounted for meaningful amounts of the total variance (87.96%). PC1 explained 78.1% of the total variance, and PC2, 9.86%. PCA of bacteria strains yielded three main clusters. The blue cluster contains most of the bacteria species (Vibrio celticus, V. gigantis, and V. gallium) that were similar in biochemical features. The red cluster contains only V. fluvialis. Finally, the green cluster contains only Aeromonas hydrophila and V. splendidus. The second PCA cluster analysis ( Figure 5) based on molecular species identification showed that the first (PC1) and second (PC2) components accounted for meaningful amounts of the total variance (58.93%). PC1 explained 30.86% of the total variance, and PC2, 28.07%. PCA of bacteria strains yielded two main clusters. The red cluster contains most of the strains that belonged to V. crassostreae, whereas the blue cluster contains V. cyclitrophicus. It is possible to observe a partial overlap of the two bacteria clusters, confirmed by the overlap of the confidence ellipses (95%) of each one, indicating similar biochemical features for certain tests. . Phylogenetic relationship among Vibrio species isolated from lesser-spotted dogfish and different Vibrio species. Phylogenetic tree was constructed using MEGAX and neighbor-joining method. A bootstrap test of 1000 replicates was performed; the bootstrap condensed tree, using a cut-off value of 50% is shown. The sequences with "G" displayed after the scientific name derived from a deposited genome.

Biochemical Features of Isolates
Morphological and biochemical analysis showed that all isolates were Gramnegative bacteria, positive for cytochrome oxidase, produced indole from tryptophan, and liquefied gelatin from enzyme gelatinase. None of the isolated strains produced lysine decarboxylase, ornithine decarboxylase, hydrogen sulfide, urease, nor were they able to utilize citrate. The enzyme tryptophan deaminase was negative in almost all strains. Generally, the isolates produced acid from glucose, arabinose, amygdalin, and mannitol. No acid was produced from amygdalin and arabinose. It is important to highlight a certain variability in biochemical features for the following test: decarboxylation of the amino acid arginine by arginine dihydrolase, detection of acetoin (acetyl methylcarbinol) produced by fermentation of glucose by bacteria utilizing the butylene glycol pathway (Voges-Proskauer test), and fermentation of sucrose and melibiose (Figure 2). The first PCA (Figure 4) based on phenotypic identification (API ® 20E) showed that the first (PC1) and second (PC2) components accounted for meaningful amounts of the total variance (87.96%). PC1 explained 78.1% of the total variance, and PC2, 9.86%. PCA of bacteria strains yielded three main clusters. The blue cluster contains most of the bacteria species (Vibrio celticus, V. gigantis, and V. gallium) that were similar in biochemical features. The red cluster contains only V. fluvialis. Finally, the green cluster contains only Aeromonas hydrophila and V. splendidus. The second PCA cluster analysis ( Figure 5) based on molecular species identification showed that the first (PC1) and second (PC2) components accounted for meaningful amounts of the total variance (58.93%). PC1 explained 30.86% of the total variance, and PC2, 28.07%. PCA of bacteria strains yielded two main clusters. The red cluster contains most of the strains that belonged to V. crassostreae, whereas the blue cluster contains V. cyclitrophicus. It is possible to observe a partial overlap of the two bacteria clusters, confirmed by the overlap of the confidence ellipses (95%) of each one, indicating similar biochemical features for certain tests.

Discussion
This study reports the finding of Vibrio crassostreae and V. cyclotrophicus in Scyliorhinus canicula isolated from different matrices (brain, blood, and kidney) and provides a description of the isolated strains. Notably, there are only a few reports in the literature on the isolation of pathogens from sharks. Therefore, our findings could clarify the pathogenic role of these bacteria despite the absence of gross pathology lesions.
The isolation process of Vibrio spp. from sharks started following the observation of external clinical signs common to many infectious diseases. For this reason, the presence of parasites, bacteria, and viruses was investigated. Parasitological and virological analyses tested negative, although most cases of diseases in sharks reported in the literature were of parasitic etiology [37]. In relation to the symptoms highlighted in the tank, a possible viral cause different from Betanodavirus cannot be excluded. However, it was not possible to carry out in-depth analyses due to the low number of samples, the limited amount of biological material available (in relation to shark size), and the scarce presence of information regarding elasmobranch viruses.
On the contrary, colonial growth was observed on agar media. Once pure cultures of the morphologically different colonies were obtained, the diagnostic techniques available in the laboratory clearly identified the isolates. Therefore, the isolates were tested by a biochemical test (API ® 20E), MALDI-TOF MS, and molecular biology techniques. Upon Gram staining, the isolates were all found to be Gram-negative rods. Since the colonies grew on TCBS and had the typical morphology of Vibrio spp., it was decided to use the API ® 20E galleries for species identification and characterization based on a biochemical test.
Biochemical tests led to the identification of bacteria of the genus Vibrio, as presumed, but also of Aeromonas hydrophila isolates. Therefore, the identification process continued by VITEK-MS, which did not lead to the identification of any of the isolates. The possible explanation of this was probably due to the absence in the database used by the instrument of reference strains comparable to those of the study. The implementation of

Discussion
This study reports the finding of Vibrio crassostreae and V. cyclotrophicus in Scyliorhinus canicula isolated from different matrices (brain, blood, and kidney) and provides a description of the isolated strains. Notably, there are only a few reports in the literature on the isolation of pathogens from sharks. Therefore, our findings could clarify the pathogenic role of these bacteria despite the absence of gross pathology lesions.
The isolation process of Vibrio spp. from sharks started following the observation of external clinical signs common to many infectious diseases. For this reason, the presence of parasites, bacteria, and viruses was investigated. Parasitological and virological analyses tested negative, although most cases of diseases in sharks reported in the literature were of parasitic etiology [37]. In relation to the symptoms highlighted in the tank, a possible viral cause different from Betanodavirus cannot be excluded. However, it was not possible to carry out in-depth analyses due to the low number of samples, the limited amount of biological material available (in relation to shark size), and the scarce presence of information regarding elasmobranch viruses.
On the contrary, colonial growth was observed on agar media. Once pure cultures of the morphologically different colonies were obtained, the diagnostic techniques available in the laboratory clearly identified the isolates. Therefore, the isolates were tested by a biochemical test (API ® 20E), MALDI-TOF MS, and molecular biology techniques. Upon Gram staining, the isolates were all found to be Gram-negative rods. Since the colonies grew on TCBS and had the typical morphology of Vibrio spp., it was decided to use the API ® 20E galleries for species identification and characterization based on a biochemical test.
Biochemical tests led to the identification of bacteria of the genus Vibrio, as presumed, but also of Aeromonas hydrophila isolates. Therefore, the identification process continued by VITEK-MS, which did not lead to the identification of any of the isolates. The possible explanation of this was probably due to the absence in the database used by the instrument of reference strains comparable to those of the study. The implementation of these databases with Vibrio spp. strains not belonging from food or clinical isolates could make MALDI-TOF MS a rapid and efficient method for diagnosing these bacteria. Moreover, the identification of strains identified as A. hydrophila by biochemical methods was not confirmed, although it is reported that VITEK-MS is able to distinguish Aeromonas spp. from some species of Vibrio (V. cholerae) [46]. For this reason, it was decided to use molecular biology techniques to obtain a more accurate identification.
For genetic identification, the 16S subunit of ribosomal RNA (16S rRNA) gene is commonly used for the study of phylogenetic relationship among bacterial taxa and identification at the species level; however, the 16S rRNA genes of vibrios are not sufficiently polymorphic to ensure a reliable identification [47][48][49]. For this reason, species-specific PCRs have been developed, and different genetic markers have been considered. Some examples are the gene encoding the α-subunit of bacterial ATP synthase (atpA) [50] and rpoB [42]; considering that the rpoB gene led to more suitable results than the 16S rRNA gene for Vibrio species identification; thus, this tool was chosen in this study for the molecular identification. The mentioned analysis allowed the identification of the isolates, respectively, as V. crassostreae (13/24; 54.2%) and V. cyclitrophicus (11/24; 45.8%). The first species was reported as a pathogen of mollusks and fish [51], while V. cyclitrophicus was reported in mollusks and in the microbiome of marine copepods [52,53]. There are no reports of these species in sharks, although bacteria of the genus Vibrio have already been reported in elasmobranchs but not specifically for S. canicula.
As further support of these results, a statistical analysis was carried out on the results of biochemical and molecular tests. The biochemical test (API ® 20E) allowed the identification of three main clusters of Vibrio species, based on biochemical features. In particular, the first PCA revealed that strains belonging to different Vibrio species (V. celticus, V. gigantis, and V. gallicus) had the same biochemical features. Considering this, it is clear how such identification is not a good diagnostic tool in fish pathology. On the contrary, the second PCA performed on biochemical features of Vibrio species identified throughout molecular tools allowed the discrimination of the biochemical features of V. crassostreae and V. cyclotrophicus, even though a partial overlap of the two main clusters was observed, indicating similar features for a certain aspect of the biochemical test (e.g., negative for decarboxylations of the amino acid ornithine, decarboxylation of the amino acid lysine urease, H 2 S production). These results support what has already been indicated in the description of the two Vibrio species isolated in this survey [54,55].
The phylogenetic analysis showed a clear separation between V. cyclitrophicus, characterized by a cluster with a high bootstrap value (97%), and V. crassostreae, even if a partial overlap was detected in the PCA analysis.
Similarly, Aeromonas hydrophila, V. splendidus, and V. fluvialis in the phylogenetic tree were clearly separated from V. crassostreae and V. cyclitrophicus even if sometimes identified as them by the API system.
Vibrio species identified in the study are generally linked to mortality in bivalve mollusks or found at the environmental level. A study by Petton et al. [56] shows how Vibrio infection is related to water temperature: In this article, it is shown that temperatures higher or lower than 16 • C affect oyster mortality due to the growth of Vibrio spp. As indicated by Sims [57], S. canicula fits perfectly in a thermal range between 14.9 • C and 17.7 • C. Thus, the tank's water temperature is ideal both for the growth of isolated Vibrio species and for S. canicula. Nevertheless, it has been found that in natural environments, the lesser-spotted dogfish prefers temperatures lower than 17 • C, for optimization of its metabolism [57].
Currently, few studies are available in the literature on elasmobranchs infectious diseases. Probably, these shortcomings are due to the lack of general elasmobranch studies, as supported by Abdul Malak et al. [32]. Among bacterial diseases, Vibrio spp. infections are reported, although in a few reports and with several shortcomings in the identification process. The present study, therefore, has the function of a report of the presence of vibrios, similar to previous research; in-depth studies will be required for details on the pathogenic mechanism of these bacteria in sharks. Nevertheless, our study made it possible to identify specifically V. crassostreae and V. cyclitrophicus, underlining the criticalities in the identification process in relation to the techniques adopted (API ® 20E, MALDI-TOF MS, molecular biology) and the need for developing specific guidelines for the identification of nonmajor pathogenic Vibrio species. Institutional Review Board Statement: The Ethics Committee was not included, as all samples were derived from a diagnostic service that the laboratory in which the authors work offers to users (http://www.izsplv.it/) Therefore, the study was derived from a routine activity for which it was not necessary to establish an Ethics Committee.