Alpha diversity and data analysis
Rarefaction curves reached asymptotic levels for all samples, indicating sufficient depth of sequencing possible (Fig 1). A healthy shrimp gut sample (CKG) revealed the most shift towards higher species richness with 1673 OTUs and a Shannon alpha diversity value of 6.40, indicating highly diverse species. In the case of the infected shrimp gut sample (IAG), 668 OTUs were obtained with a lower Shannon alpha diversity value of 4.50, indicating less diverse OTUs. Similarly, the infected shrimp gut sample (IEG) showed a lesser number of OTUs when compared to the control at 1398 OTUs with more diverse species (Table 2). Our work indicated no significant difference in the shrimp gut microbiota between healthy and EHP-infected samples (P= 0.667). Intriguingly, we found higher representation of certain OTUs in CKG when compared to IAG and IEG, throughout the class, order, family and genus levels, respectively. A similar case of no significant difference in microbiota was reported in WSSV-infected shrimp samples
(Wang et al., 2019). Holt et al., (2021) also showed that diseased shrimp larvae and blue shell syndrome-affected shrimps indicated no significant difference between the healthy and infected groups.
Taxonomic classification of gut bacterial communities
Proteobacteria and Firmicutes maintained the most abundant bacterial status irrespective of the healthy or diseased state of the shrimp. The microbial community analysis revealed that the taxa Firmicutes accounted for the highest abundance for CKG with 39.62% abundance, IAG and IEG accounting for 39.56% and 22.64%, respectively. Proteobacteria were the most abundant taxa for the sample IAG at 44.75%, followed by CKG and IEG at 28.48% and 27.82% (Fig 2). These results were consistent with an earlier study where post-larval shrimps treated with
Bacillus subtilis also exhibited Proteobacteria as the most abundant phylum, similar to their control group
(Cao et al., 2020).
The phylum Planctomycetes has been shown to be one of the dominant OTUs in CKG with an abundance of 12.05% and the same bacterial abundance level was lower in IAG and IEG, showing an abundance percentage of 2.05% and 1.79%, respectively. Planctomycetes also showed a relatively lower abundance level in the shrimp gut during an ongoing Acute Hepatopancreatic Necrosis Disease (AHPND) infection
(Chen et al., 2017). Planctomycetes also form a superphylum with Verrucomicrobia and Chlamydaie, with a unique compartmentalized cell plan for their prokaryotic organization
(Lee et al., 2009). In accordance with this, Verrucomicrobia also showed a lesser abundance (0.57% in IAG, 1.25% in IEG) compared to the healthy sample (3.26% in CKG). Verrucomicrobia, which resides in the intestinal mucosa, has been reported to aid in processing complex polysaccharides and enhance the activity of the immune response in the gut
(Cardman et al., 2014; Martinez-Garcia et al., 2012). Another bacterial phylum, Actinobacteria, which has a role in detoxification and protection against pathogens through its biofilm production
(Anandan et al., 2016), also showed variation in its richness, with a higher abundance of 11.79% in CKG and a lesser abundance of 2.54% in IAG and 6.60% in IEG.
Bacterial phyla such as Fusobacteria and Spirochaetes showed a reverse trend of higher abundance in IAG and IEG. Spirochaetes showed an extremely low abundance of 0.007% in CKG, with a higher abundance in IAG and IEG at 3.57% and 6.61%, respectively. Spirochaete potential to cause infection in brine shrimp and artemia has been demonstrated (
Tyson, 1975). Similarly, Fusobacterium showed a lower abundance of 0.17% in CKG and a higher abundance of 2.82% and 1.63% in IAG and IEG, respectively.
Wang et al., (2019) showed that the phylum Fusobacteria had a higher abundance response with respect to WSSV infection, which is similar to our study and the data shows an increased representation of potential pathogenic bacterial phyla during EHP infection. The relative bacterial abundance observed at class and order levels is listed in Table 3.
At the family level, one of the surprising findings in our study is that,
Lactobacillaceae, a major probiotic bacterium
(Oscarsson et al., 2021), had a high level of abundance in CKG, when compared to the infected samples.
Lactobacillaceae ranked number one at the family level with an abundance of 30.20% in CKG, whereas IEG and IAG showed very low abundance levels of 1.70% and 0.29%, respectively. This is a clear indication of the involvement of EHP in bringing down beneficial bacteria in the host.
Hjelm et al., (2004) showed
Rhodobacteraceae members limit the growth of
Vibrio spp. with a higher abundance level in the midgut of
P.vannamei (Pilotto et al., 2018). In this study,
Rhobacteraceae ranked higher next to
Lactobacillaceae in CKG, with an abundance of 10.50%. In contrast, IAG and IEG showed lower abundance levels of 2.47% and 1.83%, respectively (Fig 3). This further strengthens the fact that the
Rhobacteraceae bacterial population is downregulated by EHP infection, resulting in adverse effects on the health status of the host.
Verrucomicrobiaceae, Bacillaceae, Planctomycetaceae, Clostridiceae, Microbacteriaceae, Pirellulaceae, and
Pseudoalteromonadaceae showed patterns of lower abundance in IEG and IAG (Table 4).
Microbacteriaceae and
Pirellulaceae can be associated as an indicator for the healthy state of the shrimp, since
Microbacteriaceae showed fluctuating abundance during the development of a shrimp’s intestinal microbiota
(Huang et al., 2016), and
Liu et al., (2018) showed a lower abundance of
Pirellulaceae when the shrimps were treated with microbial agents.
Clostidiceae, which is reported to aid in protein digestibility
(Bermingham et al., 2017), exhibited a lower abundance in IAG and IEG, and it could be speculated that a lower abundance of such beneficial bacteria may affect the shrimp’s ability to digest complex proteins, in turn affecting its overall growth.
In relation to genus level, three unclassified bacteria with considerable abundance in CKG also followed the lesser abundance trend. An unclassified genus from the
Rhodobacteraceae family showed an abundance level of 9.40% for CKG, 2.08% for IAG and 1.30% for IEG. Another unclassified genus from the
Pirellulaceae family showed an abundance level of 9.40% for CKG and 2.08% for IAG (Fig 4). Similarly, unclassified bacteria from the Actinomycetales (anaerobic bacteria) order displayed an abundance level of 6.06% in the healthy sample (CKG), 1.52% in IAG, and 0.33% in IEG.
Arcobacter,
Fusibacter,
Enterococcus and
Spirochaeta were also present at an abundance of 21.3%, 16.4%, 14.8%, and 3.5%, respectively, in the sample, IAG. It is worth noting that certain species within the
Arcobacter genus are known to be zoonotic, causing bacteremia and gastroenteritis in humans
(Uljanovas et al., 2021). IAG and IEG also showed a higher representation of an unclassified bacterium from the
Enterococcus genus, with abundances of 11.8% and 0.1%, respectively, compared to a lower abundance of 0.01% in the healthy gut sample, CKG. Similar to
Arcobacter,
Enterococcus also has the potential to cause urinary tract infection and other diseases in humans
(Said et al., 2021).
In terms of species-level,
Pediococcus acidilactici, being an important lactic acid bacteria known for its probiotic activity that also exerts antagonism against microorganisms such as enteric pathogens through its bacteriocins and lactic acid secretion (
Daeschel and Klaenhammer 1985;
Porto et al., 2017), ranked first at species level with an abundance level of 28.66% in CKG. An extremely low abundance was observed in IAG (0.004%) and IEG didn’t show any signs of
P.acidilactici in its microbial composition.
Heat map with hierarchical clustering
Lactobacillaceae formed two subclusters with
Vibrionaceae and an unassigned bacterium, along with
Rhodobacteraceae and
Pirellulaceae, which are prevalent in CKG and less abundant in IAG and IEG. Similarly,
Verrucomicrobiaceae and
Moraxellaceae formed subclusters with
Bacillaceae, Microbacteriaceae, and
Pseudoalteromonadaceae. From the dendrogram, we can see that the samples CKG and IAG are clustered together, indicating a closer relationship when compared to IEG, even though the relative abundance level for the microbiota in infected samples is on the lower side when compared to the healthy sample (Fig 5). There is a distinct difference between the EHP-infected samples themselves, where the sample IEG showed other OTUs (not shown here) that were not seen in the samples CKG and IAG, indicating uniqueness among the EHP-infected samples.