Isolation of vibrionaceae groups and determination of diversity
As per the farmer’s information during sampling, the infected shrimps had shown abnormal swimming, loss of appetite, reddish body discoloration and necrosis of hepatopancreas which lead to the mortality of shrimps from both farms (Supplementary file ESM Fig. 1). More than 100 colonies each were segregated on TCBS agar plates from different parts of shrimp samples of two farms. Majority of the colonies were green colored and others exhibited yellow and black color. Based on the unique features like color, size, shape, margin etc. suspected vibrio colonies were selected from the pooled shrimp samples of both farms. While assessing the salt tolerance of the isolates, all the isolates have grown with NaCl up to 12% and the growth was arrested on further salinity. This indicates that the isolates are halotolerant which might be because they inhabit brackish water environment and however, not halophilic. All the recovered isolates were found as Gram negative and motile. On phenotypic characterization, the isolates have shown typical biochemical features and based on that dendrogram with 5 clusters for farm 1 with a distance coefficient of 0.50 to 1 (Supplementary file ESM Fig. 2) and 6 for farm 2 with 0.48 to 1 (Supplementary file ESM Fig. 3) were generated using NtSys software. Thus the representative 5 isolates from farm 1 and 6 from farm 2 were selected for further characterization.
According to Primer E software, the Shannon wiener diversity index of the recovered isolates was found to be H’ (loge) = 1.506 and 1.69 for farm 1 and 2 respectively. This in turn is determined by the factors such as species richness (d = 0.868 for farm 1 and 1.086 for farm 2), evenness (J’= 0.935 for farm 1 and 0.943 for farm 2) and dominance (1-lambda = 0.767 for farm 1 and 0.809 for farm 2). Thus microbial diversity is found to be higher in farm 1 compared to farm 2.
On sequence analysis of the 16SrRNA gene of the representative isolates, Vibrio alginolyticus, V. cholera, V. fluvialis, V. mimicus and Shewanella algae were identified from the samples of farm 1 and Vibrio parahaemolyticus, V. cholera, Aeromonas hydrophila, A. salmonicida, Citrobacter freundii and S. algae from that of farm 2. Phylogenetic tree (Fig. 1) consisting 5 operational taxonomic units (OTUs) for farm 1 and 6 for farm 2 was constructed on comparison with similar 16rRNA gene sequences from the Genbank databases. The distance coefficient seems to be low for farm 1 and 2 with 0.03 r to 0.0r and 0.04r to 0.0r respectively since the isolates are of related groups. The phylogenetic tree of farm 1 isolates comprised mainly two branches separated at the distance coefficient of 0.03r. All the Vibrio sp were clustered together in first branch with V. cholerae and V. mimicus under one sub-branch separated from V. fluvialis and V. alginolyticus at 0.02r under another one. While S. algae was clustered alone under second branch. While the distance coefficient of phylogenetic tree constructed with farm 2 isolates is 0.04r with main two branches. The aeromonads group and Shewanella were grouped under one branch with separation at 0.04r from the second branch consisting Citrobacter and vibrios with the similarity coefficient of 0.03r. The overall mean distance, determined using Mega X software, between the 16srRNA gene sequences of farm 1 and farm 2 isolates was found to be 0.04 and 0.07 respectively.
Even though the intention of the present study was to recover vibrios from the infected shrimps, other microbes belonging to aeromonads (A. hydrophila and A. salmonicida) and enterobacteraceae groups (Citrobacter freundii) along with Shewanella algae were also resolved on TCBS medium. As noticed in previous studies, although the medium is selective for vibrios, it may promote the growth of other organisms like Aeromonas, Shewanella, Pseudomonas etc. and can be easily misinterpreted as Vibrio by showing similar morphology and biochemical test results (Canellas et al. 2021). Shewanella was also reported to grow on TCBS medium as black colored colonies (Wang et al. 2009) while on sub culturing it was difficult to get differentiated from vibrios in our study. Hence molecular identification became essential for the genus and species level confirmation.
The most predominant organism isolated from all the infected parts of shrimps belonging to both farms were Vibrio cholera. It was followed by V. alginolyticus and V. fluvialis in the case of farm 1 and V. parahaemolyticus from farm 2. While the number of colonies of Aeromonas and S. algae recovered from the infected samples was less. It is reported that non-vibrio species like Aeromonas and Shewanella could grow on standard TCBS medium but at lower growth rate compared to vibrios (Valente and Wan 2021). The diversity of Vibrio sp. recovered from farm 1 was high compared to farm 2 while the latter showed increased overall diversity. The vibrios segregated from both farms such as V. cholerae, V. alginolyticus, V. parahaemolyticus, V. fluvialis and V. mimicus were noticed to cause significant diseases in decapod crustaceans (Valente and Wan 2021). Most of the vibriosis causing bacteria are found to be opportunistic in shrimps and prawns and could induce infections when environmental stresses enhanced. Even though V. cholera is ubiquitous in aquatic environments, it was observed to cause severe infections like yellowing of legs and white feces syndrome in fresh water cultured white leg shrimps (Penaeus vannamei) (Cao et al. 2015). Hepatopancreatic necrosis disease induced by V. alginolyticus and V. parahaemolyticus is severe and prevalent in all kinds of cultured shrimps (Morales-Covarrubias et al. 2018). In juveniles and adults of crustaceans, shell diseases caused by V. fluvialis, V. mimicus, V. alginolyticus and V. parahaemolyticus are also relevant in severe economic losses (Valente and Wan 2021). In addition to vibrios, A. hydrophila was also recently reported as an emerging pathogen in causing mass mortalities of Litopenaeus vannamei (Zhou et al. 2019). Hence the occurrence of aeromonads within the infected shrimps might also have contributed pathogenicity to the shrimps. The coinfection of pathogenic Aeromonas, Vibrio and Citrobacter in crustaceans were already reported in earlier studies (Thancharoen et al. 2019) and which is in close agreement with the present study. Another isolate which was resolved from both farms was S. algae and was reported for the first time by Cao et al. (2018) as emerging pathogen in causing black spot disease in white leg shrimps. Despite the ubiquitous occurrence of S. algae in high salinity aquatic environments, the reports of shrimps with this species infection is less; while a novel species, S. khirikhana was recently isolated from shrimps in causing early mortality syndrome by Prachumwat et al. (2020). Hence more studies are essential to recover novel species from crustaceans utilizing appropriate selective media for the pathogens. Since diverse species of vibrios and other related pathogens were recovered from the infected shrimps, their response towards different classes of antimicrobials need to be monitored for the effective treatment and management measures.
Antibiogram Profiling
The antimicrobial resistance pattern of the isolates from farm 1 revealed that they were resistant towards a maximum of 21 and a minimum of 9 antibiotics and those from farm 2 showed a maximum of 15 and minimum of 9. The resistance and sensitivity pattern of all the isolates resolved from both shrimp farms towards 35 tested antibiotics are represented as heat map in Fig. 2. Out of the 35 antibiotics tested, all the isolates were found to be resistant to tetracycline and ampicillin and all were sensitive to erythromycin, ciprofloxacin and cefazolin. Besides, amikacin, meropenem and tobramycin were found to ineffective against all the isolates from farm 1 and vancomycin and all the quinolones group were ineffective for any of the isolates from farm 2. The details showing the percentage of isolates from both farms exhibiting resistance towards 15 different antibiotic classes are given in Table.1 and figures are depicted in Fig. 3. In the case of isolates from farm 1, hundred percent of the isolates have shown AMR towards betalactams, first and third generation cephalosporins and tetracyclines and more than fifty percent towards carbapenems, aminoglycosides, quinolones and phenicols; while least percent to macrolides, glycopeptides, nitrofurans, macrolides and sulphonamides. While, those from farm 2 exhibited hundred percent resistance to beta lactams and tetracyclines and greater than fifty percent towards all the other classes except glycopeptides and first and third generation cephalosporins.
Amidst the 15 classes tested, the isolates from farm 1 possessed AMR towards at least one antibiotic of minimum 8 classes and a maximum of 12, while those from farm 2 have shown minimum of 7 and maximum of 12 classes. The graph showing the percentage of AMR shown by each isolate from the samples of both farms towards different classes of antibiotics is given in Fig. 3. The most resistant isolate recovered from the shrimp samples was V. alginolyticus, recovered from farm 1 sample, by showing resistance towards 21 antibiotics belonging to 11 different classes with the MAR index of 0.6. It is to be pointed out that they were found to be insusceptible towards all the new generation cephalosporins and is challenging. The second most resistant isolate, V. mimicus was also recovered from the farm 1 shrimps with the MAR index of 0.54 by exhibiting AMR towards 19 antibiotics of 12 different classes. This is followed by V. fluvialis from farm 1 and V. parahaemolyticus from farm 2 with the MAR index of 0.43 by expressing resistance towards 15 antibiotics of 8 and 9 various classes respectively. The least resistant isolates viz. Shewanella algae and Vibrio cholera segregated from farm 1 were found to have the MAR index of 0.34 and 0.25 respectively. Meanwhile, V. cholerae isolated from farm 2 displayed MAR index of 0.34 followed by 0.31 for A. hydrophila and 0.28 for A. salmonicida. However, C. freundii from farm 2 showed resistance towards 13 antibiotics of 9 different classes with the MAR index of 0.37. The most sensitive isolate compared to other tested ones was S. algae from farm 2 with the MAR index of 0.35. Thus all the isolates recovered from both farms exhibited MAR index of > 0.2, which highlighted the increased contamination of antibiotic residues in both systems (Preena et al. 2020b). The minimum inhibitory concentration of most of the resistant antibiotics was found to be higher than 256mcg/ml and the details are given in supplementary file ESM.Table.1 and Table. 2.
Since V. alginolyticus is one of the most significant pathogens in shrimp cultured systems, their resistance towards multiple drugs with higher MAR index raise major concerns. Multiple antibiotic resistant V. alginolyticus is already reported as a worldwide threatening species to human and animal health and food safety (Yin et al. 2022). Yasir et al. (2020) reported the draft genome sequence of V. alginolyticus from marine system, which was found to comprise multiple ARG linked with multidrug resistant efflux pumps and highlighted their major role in AMR dissemination. In addition to V. alginolyticus, V.parahaemolyticus is also a significant zoonotic pathogen found among the shrimps and hence their multidrug resistance also presents potential risks to human health (Ina-Salwany et al. 2019). The existence of AMR among V. parahaemolyticus isolates is evidenced in several studies and pointed out their role as reservoirs of ARG in aquaculture industry (Reyhanath et al. 2014; Letchumanan et al. 2015). Similar to the present study, Banerjee et al. (2012) demonstrated the coexistence of multidrug resistant V.alginolyticus, V.cholerae, V.parahaemolyticus and V.mimicus along with other pathogens in diseased cultured Litopenaeus vannamei. Most of the vibrios isolated in the study, especially V.cholerae, V.parahaemolyticus and V.fluvialis are usually considered as opportunistic pathogens and reported to cause disease outbreaks when shrimps are stressed and also found to be responsible for major food borne and zoonotic diseases (Kitiyodom et al. 2010). Hence the presence of enhanced antimicrobial resistance in such zoonotic pathogens pose threats even to human health.
Prevalence of antibiotic resistance in virulent aeromonads, especially A. hydrophila among the cultured fishes, ornamental fishes and crustaceans were noticed in several studies earlier and is in accordance with this study (Vivekanandhan et al. 2002; Preena et al. 2019a, b; Preena et al. 2020c; Preena et al. 2021). However, although A. salmonicida is found to cause profound infections, like furunculosis, in salmonid fish, it is least reported in fresh water or brackish water fish shrimps. Likewise, infection of shrimps with Citrobacter and Shewanella is also least reported. Hence it is surprising to notice the recovery of those multi drug resistant strains from shrimps under study and need further characterization. Nevertheless, their antimicrobial resistance pattern indicated that the prevalence of such least reported bacteria also play a key role in spreading AMR among the shrimp pathogens. Higher MAR index of putative pathogen, S. algae was recently reported from marine environments also corroborated with the present study (Ibrahim et al. 2021). Even though diverse bacteria were recovered and an overall increased AMR and uniform resistance pattern were observed among the samples from farm 2 compared to farm 1, the MAR index was found to be higher in the case of farm 1 with V. alginoyliticus and V. mimicus as the most resistant ones.
Plasmid Curing And Antimicrobial Susceptibility Studies
On analyzing the plasmids for determining the occurrence of ARG, all the isolates were found to have plasmids (> 10kb). Hence all the isolates were treated with various concentration of SDS and 0.2mg/ml was found to be optimum for plasmid curing for every isolate. The high frequency of plasmid elimination efficiency of SDS was reported in the earlier period onwards (Tomoeda et al. 1974) and because of the safe handling and easy disposal compared to ethidium bromide and acridine orange, SDS was selected for the present study (Letchumanan et al. 2015). After curing, plasmid elimination was confirmed by re-extraction followed by visualization on 1% agarose gel and assured that cells were devoid of plasmids. Altered resistance phenotype was observed following repeated antimicrobial sensitivity tests and a drastic change in MAR index was noticed for all organisms. V. cholerae and V. mimicus from farm 1 have lost antimicrobial resistance completely after plasmid elimination and sensitive towards all antibiotics. V. alginolyticus, V. fluvialis and S. algae were also found to be susceptible to all antibiotics except one or two antibiotics and exhibited a marked changed in MAR index. The case is not different with farm 2 also, where all aeromonads and Shewanella possessed the MAR index of 0 after plasmid curing. V. parahaemolyticus, V. cholerae and C. freundii also showed a great alteration in resistance phenotype. The antimicrobial resistance phenotype along with MAR index before and after plasmid curing is given in Table.2. The information generated denotes that most of the antimicrobial resistance genes in majority of the isolates from both farms are of plasmid borne. Even the samples were collected from two different geographic locations, an increased antimicrobial resistance pattern was observed ubiquitously. This might be due to the global indiscriminate usage of prophylactic drugs and possible rapid AMR spread through horizontal gene transfer.
It is well known reported that most of the vibrios harbor mobile genetic elements like plasmids, integrons, integrating conjugative elements, transposons etc which could acquire genetic determinants of AMR (Dutta et al. 2021). Aeromonas and Citrobacter were noticed to possess plasmids of various sizes and are well known for their capability to acquire multiple drug resistance genes within their mobile genetic elements (Fang et al. 2021).The horizontal movement of antibiotic resistance gene harboring plasmids between A. hydrophila and other strains like E.coli was demonstrated in previous studies and bring out their significance in rapid dissemination of AMR from aquatic to terrestrial environment through fish food (Stratev and Odeyemi 2016). The mobilizing capacity of multidrug resistance genes in S. algae and various vibrio strains and their active role as the reservoir of AMR in aquaculture farms and further spread to clinical pathogens was clearly determined through resistome analysis by Zago et al. (2020). It is also stated that antimicrobial resistant pathogens, especially Vibrio and Aeromonas released in to the shrimp cultured ponds have the capability to become predominant strains in the system, which raise the sudden chance of AMR circulation (Thornber et al. 2020). In addition, virulence of pathogens can be further enhanced by the integration of antimicrobial resistance gene harboring plasmids, which highlights the contribution of ARGs in causing diseases (Yu et al. 2012). Thus our study is in line with the previous studies proving the higher antimicrobial resistance index of the isolated vibrios and other bacteria and their potential in spreading diseases and AMR to the whole ecosystem.
As discussed earlier, due to the excess of use of antibiotics, shrimp aquaculture system could act as one of the major natural hub of multi drug resistance genes. It is evident in a recent study through advanced high-throughput sequencing techniques that there exists a strong correlation between bacterial diversity, antimicrobial resistance genes and mobile genetic elements in shrimp aquaculture environments (Fang et al. 2021). The prevalence of AMR among the shrimp pathogens raise major concerns to the human health also through the consumption of such under cooked shrimps. There is also a chance of exchange of ARG from aquatic bacteria to human gut microbiome and further to clinical pathogens through HGT thereby establishing AMR to a great extent in humans and make the treatment difficult (Thornber et al. 2020). This is exemplified by the laboratory conjugation experiments in which high transfer frequency of ARG from marine strains to human gut flora like E. coli has been demonstrated (Pepi and Focardi 2021). It can be concluded that the shrimp cultured system could serve as active “hotspots” of ARGs and shrimps as delivery vehicle in contributing AMR from aquatic pathogens to human. In this scenario, the attempt of this present study is of utmost importance to unravel the mode of transfer of AMR among the pathogens in shrimp cultured systems. Hence appropriate investigation systems and continuous surveillance programs need to be implemented for the timely detection of antibiotic resistant pathogens and adoption of preventive measures. Thus the fight against global AMR menace can be successfully overcome by considering the uniform fitness of all living beings under “One-Health” umbrella.