Bacterial community analysis identifies Klebsiella pneumoniae as a native symbiotic bacterium in the newborn Protobothrops mucrosquamatus
BMC Microbiology volume 23, Article number: 213 (2023)
The study of the native microbiome of organisms is crucial. The connection between the native microbiome and the host affects the formation of the innate immune system and the organism’s growth. However, the native microbiome of newborn venomous snakes has not been reported. Therefore, we aimed to determine the oral and skin microbiomes of newborn Protobothrops mucrosquamatus.
We performed 16 S full-length sequencing on 14 samples collected from 7 newborn P. mucrosquamatus individuals, specifically targeting their oral and skin microbiomes. In terms of the oral and skin microbiome, the main species were Klebsiella pneumoniae lineages. According to subspecies/species analysis, the proportion from highest to lowest was K. quasipneumoniae subsp. similipneumoniae, K. pneumoniae subsp. pneumoniae, and K. pneumoniae subsp. rhinoscleromatis. These three bacteria accounted for 62.5% and 85% of the skin and oral activity, respectively. The oral microbiome of newborn P. mucrosquamatus did not comprise common bacteria found in snakebite wounds or oral cultures in adult snakes. Therefore, the source of other microbiomes in the oral cavities of adult snakes may be the environment or prey. Functional Annotation of the Prokaryotic Taxa analysis showed that the skin/oral native microbiome metabolism was related to fermentation and human infection owing to the dominance of K. pneumoniae lineages. The characteristics of K. pneumoniae may impact the development of venom in venomous snakes.
The results of the native microbiome in the oral cavity and skin of newborn P. mucrosquamatus demonstrated that the habitat environment and prey capture may affect the composition of bacteria in adult snakes. We hypothesized that the native microbiome influences newborn venomous snakes and that K. pneumoniae lineages related to citrate fermentation may play a role in venom growth. However, further verification of this is required.
With the development of biotechnology, identifying and understanding microbiomes has become a famous field worldwide. In addition to emphasising the human microbiota, other animal or environmental microbiota have also been investigated [1,2,3,4]. The interaction between microbiomes and the host is an important research direction, as microbiomes play critical roles in nutrient acquisition, immunity, and disease processes that benefit or harm their hosts [3, 5,6,7,8]. Hence, studies on understanding host microbiomes could address many questions associated with the cross-talk between the host and microbial communities, thereby providing many novel findings. The limitations of microbiome studies include time-consuming culture methods and the low efficiency of sequencing technology before the application of next-generation sequencing (NGS) technology [9, 10]. The 16 S ribosomal RNA (rRNA) gene is a powerful tool for identifying bacteria; thus, it is essential to use NGS to harvest 16 S rRNA sequences from samples and classify them with 16 S rRNA microbial inventory approaches for microbiome studies . Third-generation sequencing technologies, such as nanopore or Pacbio single-molecule real-time (SMRT) sequencing, have recently provided high-resolution taxonomic and species-level classification [11, 12]. With increasing studies on microbiomes based on NGS or third-generation sequencing, there is increasing knowledge of the wide microbial taxonomic diversity of environments, humans, mammalian organisms, fish, etc. [6, 13,14,15]. However, studies on wild non-mammalian vertebrate microbiomes are limited .
Squamates are diverse and species-rich animals found in various habitats, ranging from tropical oceans to temperate mountaintops . They have played a critical role in the food chain and ecological recycling; hence, they can be used together with their derivatives to study medical innovations or understand evolutionary biology [3, 17,18,19]. Most research on the microbiomes of squamate reptiles has concentrated on the gastrointestinal tract of the same host species. These studies revealed gut microbial communities’ evolutionary factors, diversity, and composition [20,21,22]. However, few studies have investigated squamate reptiles’ oral or skin microbiomes [3, 20,21,22].
Snakes are one of the most abundant and charismatic squamates worldwide. According to venom characteristics, the snake can be classified as non-venomous or venomous. Snake envenomation causes 20,000 deaths annually in tropical areas. Some patients may develop cellulitis, tissue necrosis, or necrotizing fasciitis complicated by bacterial infection [23,24,25,26]. Whether the snakes are venomous or not, bacterial infections after snake envenomation are a significant global public health issue in neglected tropical diseases [27,28,29]. Six venomous snake species are common in Taiwan, including Naja atra, P. mucrosquamatus, Deinagkistrodon acutus, Trimeresurus stejnegeri, Bungarus multicinctus, and Daboia siamensis, posing a clinical management challenge [24, 30]. Wound infection has been frequently investigated in patients with snakebites, especially venomous snakebites. Further treatment, including antivenom, antibiotics, or surgical intervention, should be performed according to the patient’s condition. Some bacteria, such as Morganella morganii, Enterococcus faecalis, Aeromonas spp., Enterobacter spp., and Pseudomonas spp., have been frequently identified in surgical wounds. In contrast, Shewanella spp., Bacteroides fragilis, Klebsiella pneumoniae, Proteus spp., Providentia spp., Serratia marcescens, and Salmonella spp. were occasionally identified [24, 26, 31,32,33]. Antibiotic resistance is another critical issue in snakebite wound infections. For example, N. atra envenomation causes the most severe wound infection and the highest surgical rate. The infectious bacterial species abundance was ranked as M. morganii, followed by A. hydrophila and Enterococcus spp. [30, 34]. M. morganii isolated from the wound of N. atra envenomation was resistant to cefazolin and Augmentin (amoxicillin/clavulanate); therefore, advanced antibiotics (e.g., fluoroquinolone) are required as an effective first-line treatment for patients [34, 35]. Thus, the characteristics and metagenomics of bacteria in snakes provides excellent research and clinical value.
There are diverse bacterial symbioses in a variety of organs in snakes; for example, Aeromonas spp., Salmonella spp., Pseudomonas aeruginosa, and K. pneumoniae were the most frequently isolated bacteria in the visceral organs of snakes . Similar to the microbiome in the visceral organs of snakes, oral microbiomes such as Aeromonas hydrophila, M. morganii, and K. pneumoniae were also commonly found in culture methods . However, using 16 S Sanger sequencing with the culture method, the bacterial species identified only represented 2% of those identified under 16 S metagenomics using NGS . Using 16 S V3-4 metagenomics revealed that the opportunistic pathogens, Escherichia coli, Aeromonas spp., Propionibacterium acnes, M. morganii, Brevibacterium aureum, B. fragilis, and Shigella spp., showed high relative abundances in the oral microbiomes of N. naja, Ophiophagus hannah, Python molurus, Laticauda laticaudata, Trimeresurus flavomaculatus, and Boiga dendrophila from the India and Babuyan Island Group [3, 22]. In contrast, another study revealed Pseudomonas (62%), Delftia (19%), and Methylobacterium (9%) as the dominant bacterial genera from oral swabs of N. atra in Taiwan . The above-mentioned studies used 16 S V3-4 metagenomic analysis, which is limited to the genus level. Native symbiotic bacterial communities, especially skin and gut bacteria, play a crucial role in sustaining host-microbe symbiosis, neurodevelopment, and immune development [37, 38]. Gut microbiome development during infancy and early childhood affects mental and physical health . With individual growth, the bacterial communities of the individual are changed by the environment and diet. The individual may then face some challenges associated with microbe-host cross-talk, such as cancer . Many studies have shown that probiotics are native symbiotic bacteria from infants or juvenile animals . The skin and oral cavity are the organs where animals first encounter external substances and microorganisms. Therefore, understanding the native microbial communities of the skin and oral cavity can help elucidate the interaction between the immune system and microorganisms, and further explore the effects of native and acquired microbial communities on individuals. This study is the first to investigate the microbiome of newborn snakes. The study aimed to determine the oral and skin microbiomes of newborn P. mucrosquamatus snakes and provide a comparison between their oral and skin native symbiotic bacterial communities using 16 S full-length sequencing. By investigating the oral and skin symbiotic bacterial microbiomes of newborn snakes, the study will provide an understanding of the interactions between symbiotic bacteria, applied microbes, and snakes, as well as the ecology or habitat impact on snakes, characteristics of symbiotic bacteria, and snakebite wound infections.
Materials and methods
Seven hatchlings (newborn snakes) from a live female P. mucrosquamatus captured from the wild in the Taipei Basin, with coordinates (25.104949502832916, 121.6283376179647), Taiwan, were used in this study. Once the eggs hatched, the snakelets were temporarily kept in a tank until there was no observable fluid, blood, or other tissue on their bodies. They were then sampled by a professional with sterile water and the animals were immediately released to the Xizhi area of Taipei Basin, where the female P. mucrosquamatus was discovered.
Before sampling, commercially available sterile swabs were soaked in a standard saline solution. We then applied the swab to the inside of the hatchling’s cheek and slowly rotated the swab clockwise for 10 s to collect the oral sample. For skin samples, the same method was repeated on the skin’s ventral, lateral, and dorsal sides for 3 s each. These swabs were then placed in a centrifuge tube and kept in a refrigerator at 4 °C before being analysed within 24 h. Each sample was numbered O.AP.003 to O.AP.009 and the corresponding skin sample were numbered S.AP.010 to S.AP.016. The Institutional Animal Ethics Committee of Come Win Biotechnology Ltd. (IACUC22010) approved the experimental protocols. The Taipei City Government approved the use of the snakes.
Genomic DNA extraction
Genomic DNA (gDNA) extraction from oral and skin swabs was performed using the QIAamp DNA Microbiome Kit (Qiagen, USA). The gDNA extraction was performed according to the manufacturer’s instructions. The concentration of extracted gDNA was determined spectrophotometrically using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc., Wilmington, DE, USA). The quality of the extracted gDNA was evaluated by electrophoresis separation (1.5% gel in Tris-acetate ethylenediaminetetraacetic acid buffer), and the purified gDNA was stored at − 20 °C for the 16 S rRNA gene sequencing analysis.
16 S rRNA gene sequencing, library construction, and microbial community analysis
Purified gDNA was used to amplify the bacterial full-length 16 S rRNA using a universal primer set (27 F: AGRGTTYGATYMTGGCTCAG and 1492R: RGYTACCTTGTTACGACTT), and the amplicons were sequenced using the PacBio SMRT sequencing platform (Pacific Biosciences Inc., San Diego, CA, USA). The above conditions were used according to the manufacturer’s instructions. Circular consensus sequence reads were obtained from raw PacBio sequencing data using the standard software tools provided by the manufacturer (Pacific Biosciences). The obtained sequences were further analysed using the Quantitative Insights into Microbial Ecology (QIIME2) software package  in R (v 4.2.1). The sequence data were trimmed to remove chimeric sequences, marginal sequence errors, and noisy sequences while picking amplicon sequence variants (ASVs) using DADA2  which is also an R package. The taxonomy classification was performed with a 97% threshold limit of similarity against the SILVA database. Additionally, beta diversity was measured based on the Bray–Curtis index, followed by a permutational multivariate analysis of variance (PERMANOVA) using phyloseq , and the obtained results were visualized with ggplot2 . The relative abundance of microbes at the genus and species levels was determined using QIIME2 and the obtained results were visualized using ggplot2 and eulerr packages . The distribution pattern of the microbial community based on similarity and dissimilarity was further assessed using statistical analysis of taxonomic and functional profiles . Additionally, linear discriminant analysis effect size (LEfSe) was performed using Galaxy software (http://huttenhower.sph.harvard.edu/lefse/) following Linear discriminant analysis (LDA) score > 2 and p-value < 0.05 to identify the differential abundance in microbiota among the experimental groups.
Microbial functional prediction
The metabolic functions of the microbial communities were predicted using representative sequences and denoised ASV abundance table data. For this purpose, we used the Functional Annotation of the Prokaryotic Taxa (FAPROTAX) pipeline by transforming the ASV table into relevant ecological and metabolic functions associated with the microbial community based on the culture representative of a strain, species, or genus [46, 47]. Furthermore, the significance of the differences in the relative abundance of annotated microbial functions among the experimental groups was analysed using statistical analysis of taxonomic and functional profiles based on a two-tailed Welch t-test (p < 0.05), followed by the Benjamini-Hochberg procedure . Pearson correlation analysis was then applied to explore the relationship between the significant enriched taxa and annotated microbial functions, with a p-value range of 0.01–0.05. The analysis was performed using IBM SPSS Statistics 24 (IBM, Armonk, North Castle, NY, USA), and the obtained result was visualized using an online software (https://software.broadinstitute.org/morpheus/).
Sequence and clustering analysis based on full-length 16 S rRNA sequencing
This study obtained a total of 821,031 sequence reads associated with the skin and oral microbiota from seven newborn P. mucrosquamatus snakes based on full-length 16 S rRNA sequencing. Among these sequences, 569,106 reads were retained following quality filtration using DADA2, which were later rarefied at a sequence depth of 25,000 to compare the microbial diversity (Fig. 1A), and assigned to 2859 ASVs based on a 97% similarity index using the SILVA database. The rarefaction analysis revealed that the rarefaction curves of all samples approached a saturation plateau, indicating that the sequencing effort was sufficient to estimate microbial richness and diversity at the 97% similarity threshold in the skin and oral microbiota of newborn P. mucrosquamatus. According to ASV analysis, skin microbiota showed higher unique features (56.6%) as compared to those of oral microbiota (39.6%), as shown in the Venn diagram (Fig. 1B). However, a small portion of ASVs was shared between the skin and oral microbiota of newborn P. mucrosquamatus, indicating a diverse microbial community for the body part. Although the Venn diagram shows a significant difference in microbial communities between the skin and oral cavity, upon inspection of the raw data, it was found that the ASVs with a proportion greater than 1% are similar, while the main difference lies in the ASVs with a proportion less than 1%, where the microbial species on the skin and in the oral cavity are different.
Oral and skin microbial diversity analysis
The alpha diversity of newborn P. mucrosquamatus associated with the oral cavity and skin was evaluated using Shannon (Fig. 2A) and Simpson (Fig. 2B) diversity indices. The alpha diversity indices showed that the skin had higher alpha diversity than that of the oral cavity, which is consistent with the ASVs analysis. The difference in alpha diversity between the oral and skin microbiota of juvenile snakes was statistically significant. Similarly, beta diversity analysis based on Bray-Curtis showed distinct microbiota clustering between the oral cavity and skin, indicating a higher degree of dissimilarity in beta diversity, as shown in Fig. 2C. Furthermore, the adonis function based on PERMANOVA revealed that the difference in beta diversity was statistically significant between the oral cavity and skin of newborn P. mucrosquamatus (p < 0.05).
Oral and skin microbial community composition and abundance analysis
Taxonomic analysis based on full-length 16 S rRNA revealed striking differences in the composition and abundance of the microbial community at the taxonomic level between the oral cavity and skin of newborn P. mucrosquamatus. A total of 26 phyla, 415 genera, and 642 species were associated with the oral cavity and skin of newborn P. mucrosquamatus. Among the 26 phyla, five unique phyla were associated with the skin (19.2%). In contrast, only one phylum was associated with the oral cavity (3.8%), as shown in Fig. 3A. However, a large proportion of phyla were common (76.9%) between the oral cavity and skin, indicating a higher similarity at the phylum level. The relative abundance of the top ten phyla associated with the oral cavity and skin of newborn P. mucrosquamatus are shown in Fig. 3B). Proteobacteria was the most common phylum in the skin and oral cavity. The second-most abundant phylum in the skin was Bacteroidetes, whereas Firmicutes was the second-most abundant phylum in the oral microbiota. At the genus level, 415 genera and 642 species were obtained, with striking differences in composition and abundance between the oral and skin microbiota. Among them, 29.9% of genera and 21.8% of species were common in the oral cavity and skin (Fig. 3C and E). However, a large proportion of the genera (60.7%) was observed in the skin, whereas the remaining 9.4% of genera were associated with the oral microbiota, indicating that the skin microbiota is more diverse compared to that of the oral microbiota at the genus level, which is consistent with ASVs and diversity analysis. Klebsiella was the most dominant genus in the oral cavity and skin (Fig. 3D). K. pneumoniae lineages were discussed as critical species in the oral cavity and skin by species-level analysis (Fig. 3F). The second-most abundant genus in the oral cavity group was Neochroococcus, whereas Enterococcus was present in the oral cavity. Taxonomic analysis revealed a total of 642 species with a striking difference in composition and abundance, similar to the genus-level comparison between the oral cavity and skin of newborn P. mucrosquamatus. The relative abundance of each sample is shown in Fig. 3. For oral samples, the bacterial composition of O.AP.008 was different from others; for skin samples, the bacterial composition of S.AP.011, S.AP.014, and S.AP.016 were different from others, especially Neochroococcus gongqingensis, which was rich in S.AP.011. In brief, bacterial diversity was higher in the skin samples than that in the oral samples.
Furthermore, we performed Principal Component Analysis (PCA) to determine the variation in the bacterial community based on composition and abundance at the species level. PCA analysis showed a distinct clustering of these species based on their uniqueness and distribution in the oral cavity and skin of newborn P. mucrosquamatus, as shown in Fig. S1. At the genus level, only 21.8% of genera were common between the oral cavity and skin (Fig. 3E). However, 64.8% of genera were associated with the skin. In contrast, only 13.4% of genera were associated with oral bacteria, indicating a more diverse microbial community associated with skin than that associated with the oral cavity.
Differential analysis at taxonomic levels associated with the oral cavity and skin
LEfSe analysis coupled with a logarithmic LDA score analysis was performed to explore the statistically differential bacterial taxa between the oral and skin microbiota of newborn P. mucrosquamatus as shown in Fig. 4A. LEfSe analysis revealed 30 statistically differential bacteria at different taxonomic levels. Among them, 21 significant bacteria were enriched in the skin, whereas only nine were statistically enriched in the oral environment. Additionally, these significantly enriched bacteria associated with the skin were higher in abundance than oral bacteria, as shown in Fig. 4B, which shows species with significant differences between groups.
Microbial functional analysis associated with the oral cavity and skin
Functional annotation analysis based on the FAPROTAX pipeline revealed 73 ecological and metabolic functions associated with the oral cavity and skin microbiota of newborn P. mucrosquamatus. We performed PCA analysis to further determine the variation in these annotated functions based on composition and abundance associated with the oral cavity and skin, as shown in Fig. S2. PCA analysis showed distinct clustering of these annotated functions based on their uniqueness and distribution in the oral cavity and skin of newborn P. mucrosquamatus.
The abundance and composition of these annotated microbial functions are presented in the form of a bubble plot in Fig. 5. The majority of these annotated microbial functions, including chemoheterotrophy, nitrate reduction, fermentation, aerobic chemoheterotrophy, human pathogens, and animal symbionts, were more abundant in the oral cavity than those in the skin.
Furthermore, we performed statistical analysis using STAMP software to explore the significant enrichment of these annotated microbial functions in the oral cavity and skin of newborn P. mucrosquamatus. Comparative statistical analysis revealed that three annotated microbial functions, including animal parasites or symbionts, fermentation, and human pathogens, were significantly enriched (p < 0.05) in the oral cavity with a higher abundance than that in the skin. However, the remaining five annotated functions: cellulolysis, nitrate respiration, nitrogen fixation, nitrogen respiration, and ureolysis, were significantly enriched (p < 0.05) in the skin with a higher abundance than that in the oral environment.
Association between the annotated functions and microbial community
Pearson correlation analysis was applied to further determine the contribution of the significantly enriched bacteria at the genus and species levels to the significantly enriched annotated microbial functions associated with the oral cavity and skin (Fig. 6A and B). The correlation analysis revealed that Klebsiella at the genus, species, and subspecies levels had a significant positive correlation (p < 0.01) with human pathogens, animal parasites or symbionts, and fermentation. In contrast, it had a significant negative correlation (p < 0.01) with cellulolysis, nitrite respiration, nitrogen fixation, nitrogen respiration, and ureolysis. Conversely, Epilithonimonas at the genus and species levels had a significant positive relationship with cellulolysis (p < 0.05), nitrite respiration (p < 0.05), nitrogen fixation (p < 0.01), and ureolysis (p < 0.01). In contrast, it was significantly negatively correlated with animal parasites or symbionts (p < 0.05) and fermentation (p < 0.05). Additionally, Dyadobacter at the genus level was significantly positively correlated with cellulolysis (p < 0.01), nitrogen fixation (p < 0.05), and nitrogen respiration (p < 0.05). However, Cellvibrio at the genus and species levels had an insignificant correlation with any annotated microbial functions.
The native microbiome plays a critical role in the innate immune response after birth. Few native microbiomes are considered probiotics and have been used in various applications [48,49,50]. Investigating the native microbiome of various organisms is important for understanding the role of microorganisms in their early growth stages. Food ingestion and contact with wildlife can provide newborn organisms with various microorganisms that can act as commensal microorganisms. This study describes the native oral and skin microbiomes of newborn P. mucrosquamatus without food ingestion or environmental contamination. This is the first study to investigate the microbiota in newborn snakes.
Most snake oral bacterial flora studies have relied on culture methods [51, 52]. A local study revealed that the dominant isolated bacterial species from the oral cavity of adult P. mucrosquamatus in Taiwan were P. aeruginosa, B. fragilis, Clostridium spp., Proteus vulgaris, and E. faecalis, whereas K. pneumoniae accounted for only 16.7% . Bacteroides spp., Clostridium spp., Enterococcus spp., Pseudomonas spp., and Acinetobacter sp. were the dominant bacterial species found in P. mucrosquamatus bite wound infections [53, 54]. The oral microbiomes of P. mucrosquamatus newborns in our study differed from those of adults and biting wound infections, indicating that the acquired environment and diet significantly impact microbiome changes.
Skin microbiomes are critical for protection against infection and maintaining health . However, few studies have focused on the reptilian skin microbiome, and most have focused on the pathogenic microorganisms of reptile skin diseases [55, 56]. A study indicated that the skin microbiome of komodo dragons (Varanus komodoensis) has a higher microbial diversity than the oral microbiome, which agrees with our findings . Fungal dermatitis is a common snake disease that shifts the skin microbiome [57, 58]. The protective microbiome is a critical concern, and the native microbiome is a key factor [59, 60]. Our results serve as a reference for the prevention of skin diseases in snakes.
Previous studies have shown significant differences between oral and skin microbiomes in reptiles and snakes . However, no such differences were found in the newborn P. mucrosquamatus; the predominant bacterial communities found in both the skin and oral cavity of newborn snakes in this study were similar. However, there are greater relative differences in the ASVs that comprise less than 1% of the skin microbiome when compared to the oral cavity. Additionally, these ASVs have more homology with environmental microbes, indicating that the skin microbiome of newborn snakes is more susceptible to environmental influences. The similarity of the dominant bacterial communities suggests that these microbes are native to P. mucrosquamatus. Consequently, further analysing the interior or exterior of the eggshell can help understand the sources of the native microbial species in the current study. Additionally, the diverse bacterial communities found in adult snakes from different countries, combined with our results, suggests an adaptive microbiome dependent on habitat, thereby aligning with the findings from a study on Komodo dragons . Unfortunately, to date, no studies have characterised newborn snake skin or oral microbial communities; hence, there is no existing knowledge to compare.
In this study, K. pneumoniae lineages (K. pneumoniae and K. quasipneumoniae) were dominant in the oral and skin microbiomes of newborn P. mucrosquamatus, suggesting a key role in the snake’s innate immune system. Therefore, it will be helpful to understand the interaction between native bacteria and the host by analysing the genomic and proteomic features of Klebsiella spp. isolated from newborn snakes. Moreover, K. pneumoniae is responsible for many hospital-acquired infections and is one of the most critical multidrug-resistant microorganisms worldwide [61, 62]. Because K. pneumoniae can be found in various environments and animals, it has been considered a public health issue and requires urgent to monitoring . Klebsiella pneumoniae has been previously subdivided into three taxonomic subspecies, K. pneumoniae subsp. pneumoniae, K. pneumoniae subsp. rhinoscleromatis and K. pneumoniae subsp. ozaenae based on the KpI phylogenetic tree . With the development of molecular technology, an increasing number of subgroups have been reported, such as K. quasipneumoniae, which is a novel lineage of K. pneumoniae . Klebsiella pneumoniae and K. quasipneumoniae are zoonotic pathogens with severe antibiotic resistance .
Klebsiella spp. play a significant role in the skin and oral microbiome of newborn P. mucrosquamatus snakes. Results from the microbial function predictions indicated high fermentation activities of K. pneumoniae and K. quasipneumoniae. Snakes use their bodies to capture prey and potentially transfer K. pneumoniae lineages into their digestive systems after swallowing. Many studies have indicated that K. pneumoniae possesses diverse citrate fermentation genes and participates in citrate fermentation [67,68,69]. Citrate is responsible for 2.3–12.9% of various snake venoms and inhibits snake venom proteases . This study raises the hypothesis that native K. pneumoniae or its lineage may play a crucial role in venom development in newborn venomous snakes; however, further research is needed to confirm this hypothesis. Therefore, isolating K. pneumoniae lineages from newborn venomous snakes is an important task to gain further insight into the interactions between the host and microbiome in snakes.
Klebsiella pneumoniae lineages play a crucial role in the skin and oral microbiome of newborn P. mucrosquamatus, accounting for more than 60% of the microbiome. A limitation of this study is the lack of isolation of K. pneumoniae from the samples. Although the oral microbiome is consistent among individuals, the skin microbiome varies. This study proposed two essential points. First, adult snake oral and skin microbiome composition is affected by the external environment or diet, resulting in soil-associated bacteria and Enterobacter being dominant. Thus, snakebite-related infections may be associated with the habitat and environment. Second, our study hypothesises that K. pneumoniae lineages may play a role in venom production and maturation during venomous snake development. This study serves as an essential reference for future research on venomous snakes.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the [Figshare] and [NIBI] repositories, [https://doi.org/10.6084/m9.figshare.23551629] and [PRJNA990549]. The samples were assigned consecutive accession numbers ranging from SAMN36268686 to SAMN36268699, while the raw sequencing data were assigned consecutive numbers from SRR25114049 to SRR25114062 in NCBI.
Functional Annotation of the Prokaryotic Taxa
Next Generation Sequencing
Permutational multivariate analysis of variance
Quantitative Insights into Microbial Ecology
Amplicon sequence variants
Linear discriminant analysis effect size
Chen JS, Hussain B, Tsai HC, Nagarajan V, Koner S, Hsu BM. Analysis and interpretation of hot springs water, biofilms, and sediment bacterial community profiling and their metabolic potential in the area of Taiwan geothermal ecosystem. Sci Total Environ. 2022;856:159115. https://doi.org/10.1016/j.scitotenv.2022.159115.
Egan SL, Loh SM, Banks PB, Gillett A, Ahlstrom L, Ryan UM, et al. Bacterial community profiling highlights complex diversity and novel organisms in wildlife ticks. Ticks Tick Borne Dis. 2020;11(3):101407. https://doi.org/10.1016/j.ttbdis.2020.101407.
Smith SN, Colston TJ, Siler CD. Venomous snakes reveal ecological and phylogenetic factors influencing variation in gut and oral microbiomes. Front Microbiol. 2021;12: 657754. https://doi.org/10.3389/fmicb.2021.657754.
Ursell LK, Metcalf JL, Parfrey LW, Knight R. Defining the human microbiome. Nutr Rev. 2012;70 Suppl 1:S38-44. https://doi.org/10.1111/j.1753-4887.2012.00493.x.
Arizza V, Vecchioni L, Caracappa S, Sciurba G, Berlinghieri F, Gentile A, et al. New insights into the gut microbiome in loggerhead sea turtles Caretta caretta stranded on the Mediterranean coast. PLoS ONE. 2019;14(8): e0220329. https://doi.org/10.1371/journal.pone.0220329.
Fang CY, Chen JS, Hsu BM, Hussain B, Rathod J, Lee KH. Colorectal cancer stage-specific fecal bacterial community fingerprinting of the Taiwanese population and underpinning of potential taxonomic biomarkers. Microorganisms. 2021;9(8): 1548. https://doi.org/10.3390/microorganisms9081548.
Qin Z, Wang S, Guo D, Zhu J, Chen H, Bai L, et al. Comparative analysis of intestinal bacteria among venom secretion and non-secrection snakes. Sci Rep. 2019;9(1):6335. https://doi.org/10.1038/s41598-019-42787-6.
Zhang J, Ge Z, Ma Z, Huang D, Zhang J. Seasonal changes driving shifts of aquatic rhizosphere microbial community structure and the functional properties. J Environ Manage. 2022;322:116124. https://doi.org/10.1016/j.jenvman.2022.116124.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7. https://doi.org/10.1038/s41587-019-0209-9.
Wensel CR, Pluznick JL, Salzberg SL, Sears CL. Next-generation sequencing: insights to advance clinical investigations of the microbiome. J Clin Invest. 2022;132(7): e154944. https://doi.org/10.1172/JCI154944.
Matsuo Y, Komiya S, Yasumizu Y, Yasuoka Y, Mizushima K, Takagi T, et al. Full-length 16S rRNA gene amplicon analysis of human gut microbiota using MinION nanopore sequencing confers species-level resolution. BMC Microbiol. 2021;21(1):35. https://doi.org/10.1186/s12866-021-02094-5.
Wagner J, Coupland P, Browne HP, Lawley TD, Francis SC, Parkhill J. Evaluation of PacBio sequencing for full-length bacterial 16S rRNA gene classification. BMC Microbiol. 2016;16(1):274. https://doi.org/10.1186/s12866-016-0891-4.
Bletz MC, Goedbloed DJ, Sanchez E, Reinhardt T, Tebbe CC, Bhuju S, et al. Amphibian gut microbiota shifts differentially in community structure but converges on habitat-specific predicted functions. Nat Commun. 2016;7: 13699. https://doi.org/10.1038/ncomms13699.
Gajardo K, Rodiles A, Kortner TM, Krogdahl A, Bakke AM, Merrifield DL, et al. A high-resolution map of the gut microbiota in Atlantic salmon (Salmo salar): a basis for comparative gut microbial research. Sci Rep. 2016;6: 30893. https://doi.org/10.1038/srep30893.
Hird SM, Sanchez C, Carstens BC, Brumfield RT. Comparative gut microbiota of 59 neotropical bird species. Front Microbiol. 2015;6: 1403. https://doi.org/10.3389/fmicb.2015.01403.
Vitt LJ, Pianka ER, Cooper WE Jr, Schwenk K. History and the global ecology of squamate reptiles. Am Nat. 2003;162(1):44–60. https://doi.org/10.1086/375172.
Calvete JJ. Venomics: integrative venom proteomics and beyond. Biochem J. 2017;474(5):611–34. https://doi.org/10.1042/BCJ20160577.
Sunagar K, Morgenstern D, Reitzel AM, Moran Y. Ecological venomics: how genomics, transcriptomics and proteomics can shed new light on the ecology and evolution of venom. J Proteom. 2016;135:62–72. https://doi.org/10.1016/j.jprot.2015.09.015.
Watanabe A, Fabre AC, Felice RN, Maisano JA, Muller J, Herrel A, et al. Ecomorphological diversification in squamates from conserved pattern of cranial integration. Proc Natl Acad Sci U S A. 2019;116(29):14688–97. https://doi.org/10.1073/pnas.1820967116.
Goldstein EJ, Tyrrell KL, Citron DM, Cox CR, Recchio IM, Okimoto B, et al. Anaerobic and aerobic bacteriology of the saliva and gingiva from 16 captive Komodo dragons (Varanus komodoensis): new implications for the “bacteria as venom” model. J Zoo Wildl Med. 2013;44(2):262–72. https://doi.org/10.1638/2012-0022R.1.
Hyde ER, Navas-Molina JA, Song SJ, Kueneman JG, Ackermann G, Cardona C, et al. The oral and skin microbiomes of captive komodo dragons are significantly shared with their habitat. mSystems. 2016;1(4):10. https://doi.org/10.1128/mSystems.00046-16.
Krishnankutty SP, Muraleedharan M, Perumal RC, Michael S, Benny J, Balan B, et al. Next-generation sequencing analysis reveals high bacterial diversity in wild venomous and non-venomous snakes from India. J Venom Anim Toxins Incl Trop Dis. 2018;24:41. https://doi.org/10.1186/s40409-018-0181-8.
Alirol E, Sharma SK, Bawaskar HS, Kuch U, Chappuis F. Snake bite in South Asia: a review. PLoS Negl Trop Dis. 2010;4(1): e603. https://doi.org/10.1371/journal.pntd.0000603.
Mao YC, Liu PY, Chiang LC, Lai CS, Lai KL, Ho CH, et al. Naja atra snakebite in Taiwan. Clin Toxicol (Phila). 2018;56(4):273–80. https://doi.org/10.1080/15563650.2017.1366502.
Sharma SK, Chappuis F, Jha N, Bovier PA, Loutan L, Koirala S. Impact of snake bites and determinants of fatal outcomes in southeastern Nepal. Am J Trop Med Hyg. 2004;71(2):234–8.
Warrell DA. Snake bite. Lancet. 2010;375(9708):77–88. https://doi.org/10.1016/S0140-6736(09)61754-2.
Abba Y, Ilyasu Y, Yusoff M, Noordin M. Bacterial co-infections in a captive Python bivittatus with septicemia. Sokoto J Veterinary Sci. 2016;14(2):67–71.
Abba Y, Ilyasu YM, Noordin MM. Isolation and identification of bacterial populations of zoonotic importance from captive non-venomous snakes in Malaysia. Microb Pathog. 2017;108:49–54. https://doi.org/10.1016/j.micpath.2017.04.038.
Weinstein SA. Non-front-fanged colubroid snakes. In: Brent J, et al. Crit Care Toxicol. Cham: Springer; 2017. p. 1–41. https://doi.org/10.1007/978-3-319-17900-1_93.
Huang LW, Wang JD, Huang JA, Hu SY, Wang LM, Tsan YT. Wound infections secondary to snakebite in central Taiwan. J Venom Anim Toxins including Trop Dis. 2011;18:272–6. https://doi.org/10.1590/S1678-91992012000300004.
Liu PY, Shi ZY, Lin CF, Huang JA, Liu JW, Chan KW, et al. Shewanella infection of snake bites: a twelve-year retrospective study. Clin (Sao Paulo). 2012;67(5):431–5. https://doi.org/10.6061/clinics/2012(05)05.
Mao YC, Chuang HN, Shih CH, Hsieh HH, Jiang YH, Chiang LC, et al. An investigation of conventional microbial culture for the Naja atra bite wound, and the comparison between culture-based 16S Sanger sequencing and 16S metagenomics of the snake oropharyngeal bacterial microbiota. PLoS Negl Trop Dis. 2021;15(4): e0009331. https://doi.org/10.1371/journal.pntd.0009331.
Mitchell MA, Shane SM. Salmonella in reptiles, Seminars in Avian and Exotic Pet Medicine. 2001;10(1):25–35. ISSN 1055-937X. https://doi.org/10.1053/saep.2001.19798.
Su HY, Wang MJ, Li YH, Tang CN, Tsai MJ. Can surgical need in patients with Naja atra (Taiwan or Chinese cobra) envenomation be predicted in the emergency department? Hong Kong Med J. 2016;22(5):435–44. https://doi.org/10.12809/hkmj154739.
Ngo ND, Le QX, Pham AQ, Nguyen NT, Ha HT, Dinh MMQ, et al. Clinical features, bacteriology, and antibiotic treatment among patients with presumed Naja bites in Vietnam. Wilderness Environ Med. 2020;31(2):151–6. https://doi.org/10.1016/j.wem.2020.01.002.
Resiere D, Olive C, Kallel H, Cabie A, Neviere R, Megarbane B, et al. Oral microbiota of the snake Bothrops lanceolatus in Martinique. Int J Environ Res Public Health. 2018;15(10): 2122. https://doi.org/10.3390/ijerph15102122.
Dhariwala MO, Scharschmidt TC. Baby’s skin bacteria: first impressions are long-lasting. Trends Immunol. 2021;42(12):1088–99. https://doi.org/10.1016/j.it.2021.10.005.
Sarkar A, Prescott SM, Dutra S, Yoo JY, Gordon J, Shaffer E, et al. Relationships of the very low birth weight infant microbiome with neurodevelopment at 2 and 4 years of age. Dev Psychobiol. 2022;64(7): e22317. https://doi.org/10.1002/dev.22317.
Alving-Jessep E, Botchway E, Wood AG, Hilton AC, Blissett JM. The development of the gut microbiome and temperament during infancy and early childhood: a systematic review. Dev Psychobiol. 2022;64(7):e22306. https://doi.org/10.1002/dev.22306.
Jomehzadeh N, Javaherizadeh H, Amin M, Saki M, Al-Ouqaili MTS, Hamidi H, et al. Isolation and identification of potential probiotic Lactobacillus species from feces of infants in southwest Iran. Int J Infect Dis. 2020;96:524–30. https://doi.org/10.1016/j.ijid.2020.05.034.
Callahan BJ, Sankaran K, Fukuyama JA, McMurdie PJ, Holmes SP. Bioconductor workflow for microbiome data analysis: from raw reads to community analyses. F1000Research. 2016;5(1492):1492.
McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4): e61217.
Wickham H. ggplot2: elegant graphics for data analysis New York. NY: Springer; 2009.
Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12(1):1–7.
Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4.
Sansupa C, Wahdan SFM, Hossen S, Disayathanoowat T, Wubet T, Purahong W. Can we use functional annotation of prokaryotic taxa (FAPROTAX) to assign the ecological functions of soil bacteria? Appl Sci. 2021;11(2): 688.
Louca S, Parfrey LW, Doebeli M. Decoupling function and taxonomy in the global ocean microbiome. Science. 2016;353(6305):1272–7.
Xiang Q, Yan X, Shi W, Li H, Zhou K. Early gut microbiota intervention in premature infants: application perspectives. J Adv Res. 2022. https://doi.org/10.1016/j.jare.2022.11.004.
Campos-Cerda F, Bohannan BJM. The nidobiome: a framework for understanding microbiome assembly in neonates. Trends Ecol Evol. 2020;35(7):573–82. https://doi.org/10.1016/j.tree.2020.03.007.
Thaiss CA, Zmora N, Levy M, Elinav E. The microbiome and innate immunity. Nature. 2016;535(7610):65–74. https://doi.org/10.1038/nature18847.
Chuang P-C, Lin W-H, Chen Y-C, Chien C-C, Chiu I-M, Tsai T-S. Oral bacteria and their antibiotic susceptibilities in Taiwanese venomous snakes. Microorganisms. 2022;10(5): 951.
Shaikh IK, Dixit PP, Pawade BS, Potnis-Lele M, Kurhe BP. Assessment of cultivable oral bacterial flora from important venomous snakes of India and their antibiotic susceptibilities. Curr Microbiol. 2017;74(11):1278–86. https://doi.org/10.1007/s00284-017-1313-z.
Lin CC, Chen YC, Goh ZNL, Seak CK, Seak JC, Shi-Ying G, et al. Wound infections of snakebites from the Venomous Protobothrops mucrosquamatus and Viridovipera stejnegeri in Taiwan: Bacteriology, antibiotic susceptibility, and predicting the need for antibiotics-A bite study. Toxins (Basel). 2020;12(9):575. https://doi.org/10.3390/toxins12090575.
Mao YC, Liu PY, Chiang LC, Lee CH, Lai CS, Lai KL, et al. Clinical manifestations and treatments of Protobothrops mucrosquamatus bite and associated factors for wound necrosis and subsequent debridement and finger or toe amputation surgery. Clin Toxicol (Phila). 2021;59(1):28–37. https://doi.org/10.1080/15563650.2020.1762892.
Ross AA, Rodrigues Hoffmann A, Neufeld JD. The skin microbiome of vertebrates. Microbiome. 2019;7(1):79. https://doi.org/10.1186/s40168-019-0694-6.
Schumacher J. Selected infectious diseases of wild reptiles and amphibians. J Exotic Pet Med. 2006;15(1):18–24.
Cheatwood JL, Jacobson ER, May PG, Farrell TM, Homer BL, Samuelson DA, et al. An outbreak of fungal dermatitis and stomatitis in a free-ranging population of pigmy rattlesnakes (Sistrurus miliarius barbouri) in Florida. J Wildl Dis. 2003;39(2):329–37.
Allender MC, Baker S, Britton M, Kent AD. Snake fungal disease alters skin bacterial and fungal diversity in an endangered rattlesnake. Sci Rep. 2018;8(1):12147. https://doi.org/10.1038/s41598-018-30709-x.
Chiu L, Bazin T, Truchetet ME, Schaeverbeke T, Delhaes L, Pradeu T. Protective microbiota: from localized to long-reaching co-immunity. Front Immunol. 2017;8: 1678. https://doi.org/10.3389/fimmu.2017.01678.
Weiland-Brauer N, Pinnow N, Langfeldt D, Roik A, Gullert S, Chibani CM, et al. The native microbiome is crucial for offspring generation and fitness of Aurelia aurita. mBio. 2020;11(6):e02336. https://doi.org/10.1128/mBio.02336-20.
Marisa BDJ, Marthie ME, Ricardo FDS, Marleen MK. Review - understanding β-lactamase producing Klebsiella pneumoniae. In: Maria Cristina O, editor. Antimicrobial Resistance. Rijeka: IntechOpen; 2015. p. Ch. 4.
Podschun R, Ullmann U. Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors. Clin Microbiol Rev. 1998;11(4):589–603. https://doi.org/10.1128/CMR.11.4.589.
Ludden C, Moradigaravand D, Jamrozy D, Gouliouris T, Blane B, Naydenova P, et al. A One Health study of the genetic relatedness of Klebsiella pneumoniae and their mobile elements in the East of England. Clin Infect Dis. 2020;70(2):219–26. https://doi.org/10.1093/cid/ciz174.
Brisse S, Verhoef J. Phylogenetic diversity of Klebsiella pneumoniae and Klebsiella oxytoca clinical isolates revealed by randomly amplified polymorphic DNA, gyrA and parC genes sequencing and automated ribotyping. Int J Syst Evol Microbiol. 2001;51(Pt 3):915–24. https://doi.org/10.1099/00207713-51-3-915.
Brisse S, Passet V, Grimont PAD. Description of Klebsiella quasipneumoniae sp. nov., isolated from human infections, with two subspecies, Klebsiella quasipneumoniae subsp. quasipneumoniae subsp. nov. and Klebsiella quasipneumoniae subsp. similipneumoniae subsp. nov., and demonstration that Klebsiella singaporensis is a junior heterotypic synonym of Klebsiella variicola. Int J Syst Evol Microbiol. 2014;64(Pt 9):3146–52. https://doi.org/10.1099/ijs.0.062737-0.
Chew KL, Octavia S, Lai D, Lin RTP, Teo JWP. Genomic characterization of Klebsiella quasipneumoniae from clinical specimens in Singapore. Antimicrob Agents Chemother. 2021;65(8):e0041221. https://doi.org/10.1128/AAC.00412-21.
Bott M, Meyer M, Dimroth P. Regulation of anaerobic citrate metabolism in Klebsiella pneumoniae. Mol Microbiol. 1995;18(3):533–46. https://doi.org/10.1111/j.1365-2958.1995.mmi_18030533.x.
Chen YT, Liao TL, Wu KM, Lauderdale TL, Yan JJ, Huang IW, et al. Genomic diversity of citrate fermentation in Klebsiella pneumoniae. BMC Microbiol. 2009;9: 168. https://doi.org/10.1186/1471-2180-9-168.
Liu M, Hao G, Li Z, Zhou Y, Garcia-Sillas R, Li J, et al. CitAB two-component system-regulated citrate utilization contributes to Vibrio cholerae competitiveness with the gut microbiota. Infect Immun. 2019;87(3): e00746-18. https://doi.org/10.1128/IAI.00746-18.
Odell GV, Ferry PC, Vick LM, Fenton AW, Decker LS, Cowell RL, et al. Citrate inhibition of snake venom proteases. Toxicon. 1998;36(12):1801–6. https://doi.org/10.1016/s0041-0101(98)00084-1.
The authors are grateful to Kun-Tzu Li and Yi-Ting Lu for their technical support.
This work was supported by the E-DA Hospital (grant number EDAHP110044 and EDAHP110022) and National Science and Technology Council, Taiwan (grant number 110-2314-B-650 -011 -MY2, EDPJ111040). The funding bodies played no role in the design of the study and collection, analysis, interpretation of data, and in writing the manuscript.
Ethics approval and consent to participate
Animal studies were approved by the Institutional Animal Care and Use Committee (IACUC) of Come Win Biotechnology Ltd. in Taiwan (Approval No.: IACUC22010). All experiments in this study were performed in accordance with relevant guidelines and regulations (ARRIVE guidelines).
Consent for publication
The authors declare they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
PCA plot showing the distinct clustering in the bacterial diversity at the genus level in the oral cavity and skin. Fig. S2. PCA plot showing the distinct clustering in the annotated microbial functions using FAPROTAX in the oral cavity and skin of newborn Protobothrops mucrosquamatus.
About this article
Cite this article
Su, HY., Hussain, B., Hsu, BM. et al. Bacterial community analysis identifies Klebsiella pneumoniae as a native symbiotic bacterium in the newborn Protobothrops mucrosquamatus. BMC Microbiol 23, 213 (2023). https://doi.org/10.1186/s12866-023-02936-4