Next Article in Journal
Deletion of the Neuronal Transcription Factor Satb1 Induced Disturbance of the Kinome and Mechanisms of Hypoxic Preconditioning
Previous Article in Journal
Interaction between Autophagy and Senescence in Pancreatic Beta Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biology
Volume 12
Issue 9
10.3390/biology12091206
biology-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparisons of the Oral Microbiota from Seven Species of Wild Venomous Snakes in Taiwan Using the High-Throughput Amplicon Sequencing of the Full-Length 16S rRNA Gene

by
Wen-Hao Lin
1 and
Tein-Shun Tsai
2,*
1
Institute of Wildlife Conservation, National Pingtung University of Science and Technology, Pingtung 912301, Taiwan
2
Department of Biological Science and Technology, National Pingtung University of Science and Technology, Pingtung 912301, Taiwan
*
Author to whom correspondence should be addressed.
Biology 2023, 12(9), 1206; https://doi.org/10.3390/biology12091206
Submission received: 30 June 2023 / Revised: 26 August 2023 / Accepted: 2 September 2023 / Published: 4 September 2023
(This article belongs to the Section Microbiology)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

There are millions of snake envenoming cases annually worldwide, which is an important but often neglected tropical disease. The morbidity and mortality from snake bites are not only caused by the venom toxins but also the secondary infection from bacteria from the snake itself or an external source. Therefore, if we do not understand the actual pathogens of infection, the expected treatment efficacy is reduced. In this study, we investigated the microbiota carried in the oral cavity of seven wild venomous snakes in Taiwan using full-length 16S rRNA next-generation sequencing. The results showed that non-pathogenic bacteria and pathogenic bacteria were similar in dominance in the oral cavity of snakes, Gram-negative bacteria were more common than Gram-positive bacteria, and the microbiota were significantly different between the Elapidae and Viperidae families of snakes. Our findings will assist in resolving the diversity of oral microbiomes in snakes and can be applied to the veterinary medicine and clinical therapy of venomous snakebites.

Abstract

A venomous snake’s oral cavity may harbor pathogenic microorganisms that cause secondary infection at the wound site after being bitten. We collected oral samples from 37 individuals belonging to seven species of wild venomous snakes in Taiwan, including Naja atra (Na), Bungarus multicinctus (Bm), Protobothrops mucrosquamatus (Pm), Trimeresurus stejnegeri (Ts), Daboia siamensis (Ds), Deinagkistrodon acutus (Da), and alpine Trimeresurus gracilis (Tg). Bacterial species were identified using full-length 16S rRNA amplicon sequencing analysis, and this is the first study using this technique to investigate the oral microbiota of multiple Taiwanese snake species. Up to 1064 bacterial species were identified from the snake’s oral cavities, with 24 pathogenic and 24 non-pathogenic species among the most abundant ones. The most abundant oral bacterial species detected in our study were different from those found in previous studies, which varied by snake species, collection sites, sampling tissues, culture dependence, and analysis methods. Multivariate analysis revealed that the oral bacterial species compositions in Na, Bm, and Pm each were significantly different from the other species, whereas those among Ts, Ds, Da, and Tg showed fewer differences. Herein, we reveal the microbial diversity in multiple species of wild snakes and provide potential therapeutic implications regarding empiric antibiotic selection for wildlife medicine and snakebite management.

1. Introduction

There are approximately 2.7 million cases of snake envenoming annually worldwide, resulting in 81,000–138,000 deaths [1]. The World Health Organization (WHO) placed it in the category A of the neglected tropical diseases in 2017 [2,3,4]. The morbidity and mortality from snakebites are caused by not only the venom toxins but also the secondary infection from pathogenic bacteria [5,6,7]. There may be abundant microorganisms within a snake’s oral cavity, including harmless environmental bacterial species and opportunistic, human-pathogenic species [8,9]. After a snakebite, the venom causes bruising around the bite site, and abscesses, pustules, or necrosis can develop, which creates a suitable environment for the bacteria from the snake’s oral cavity and saliva [4,7,10,11,12]. Many bacterial species could cause secondary infection, including Gram-positive, Gram-negative, and anaerobic bacteria. Anaerobic bacteria can cause gas gangrene after infecting the bite wound; Gram-positive bacteria can cause infection symptoms, such as an abscess, cellulitis, sepsis, meningitis, and urinary tract infections; and Gram-negative bacteria can cause several infection symptoms, such as gastroenteritis, sepsis, respiratory infection, meningitis, diarrhea, fever, and soft tissue infection [7,8,12,13,14,15,16,17,18]. Although anti-serum therapy can reduce venom toxicity to human tissue and the organ system, it cannot prevent a secondary infection caused by bacteria. The treatment for bacterial infections is typically empiric antibiotic administration, but the microbiota of the infected wound and oropharynx of the culprit snakes must be properly established. If the actual pathogens of infection are unknown, the expected treatment efficacy would decrease, and antibiotic resistance would likely be induced [4,9,19].
In addition to potentially infecting humans by snakebites, the oral bacteria of snakes also affect snakes themselves [16,20,21]. Snakes may experience oral wounds or infection when feeding prey or contacting foreign substances in the environment [22,23,24]. The bacteria may invade through the respiratory and gastrointestinal tracts, causing systemic infection and even death [25,26,27]. In reptiles, Gram-negative bacteria, such as Pseudomonas, Aeromonas, Proteus, and Escherichia, may be common taxa that cause infection [24,27]. Pseudomonas aeruginosa, Providencia rettgeri, and Stenotrophomonas maltophilia are often the dominant bacteria in the oral cavity of snakes suffering from stomatitis. In addition, Klebsiella, Acinetobacter, Citobacter, Chlamydia, Morganella, Staphylococcus, Bacteroides, Fusobacterium, Clostridium, Mycobacterium, and Peptostreptococcus are also common pathogenic bacteria in the oral cavity of snakes [24,25]. Previous studies showed that the oral cavity microbiomes may vary by the species, season, habitat, health situation, and predation strategy of snakes [7,9,12,13,16,17,20,28,29]. Even if the phylogenetic relationship among snake species is close or they are of the same species, they may have different oral cavity bacterial species if their habitat and activity range are not similar [9,28]. In terms of the development of wildlife medicine, only by better understanding the oral bacteria can the subsequent therapy and empiric antibiotic selection be precisely and efficiently made on snakes.
Owing to the lack of information on the complete microbiota in the oral cavity of Taiwanese snake species (except [30]) and their relationship with environmental factors, this study conducted high-throughput sequencing (next-generation sequencing; NGS) of the full-length 16S rRNA gene to identify the full compositions of microbiota within the oral cavity of seven Taiwanese snake species belonging to the Viperidae and Elapidae families. The nine hypervariable regions (V1–V9) of bacterial 16S ribosomal RNA genes can be targeted to identify bacterial taxa in 16S amplicon NGS studies [31,32], in contrast to the short-read NGS used in some previous studies, which used short fragments, such as V3–V4 regions, as target references [33]. Without the sequence information from other hypervariable regions, information at the species level usually cannot be fully acquired by short-read NGS, whereas long-read sequencing can fix the issues revealed with short-reads, such as genome-wide repeats and structural variant detection [34]. In this study, we conducted long-read sequencing by virtual PCR with primers covering the V1–V9 region of the 16S rRNA gene as references. We aimed to investigate not only the diversity of snake oral cavity microbiota within and among species, but also their correlations with altitude. We also discussed the results of previous studies compared to the present study, regarding different snake species, collection sites (in clinic, out-of-hospital, in animal house, or in wild), sampling tissues (oral cavity, gut, or bite wound), culture dependence, and analysis methods (CBtest, mass spectrometry, or 16S rRNA sequencing).

2. Materials and Methods

2.1. Sample Collection

A total of 37 snakes from seven venomous species were sampled in this study, including Naja atra (Na; n = 6), Bungarus multicinctus (Bm; n = 5), Protobothrops mucrosquamatus (Pm; n = 6), Trimeresurus stejnegeri (Ts; n = 5), Daboia siamensis (Ds; n = 5), Deinagkistrodon acutus (Da; n = 5) and alpine Trimeresurus gracilis [35] (Tg; n = 5) (Figure 1a–g), which were captured from the wild in Kaohsiung, Pingtung, Taitung, Yilan, and Nantou in Taiwan (Figure 2). To avoid the original microbiota within the oral cavity being affected by traffic and transporting procedures, we sampled the snakes in as soon as possible (2.4 ± 2.3 [mean ± SD] days) after collection from the wild. To prevent cross-contamination, we opened the mouth of the snake using a sterilized, hollowed out wooden stir stick and used sterilized cotton-tipped swab sticks to collect oral samples (Figure 1h) [36]. Samples were taken by rotating the cotton tip of the swab on the floor of the oral cavity between the larynx and mandibular teeth from each snake after securing the head [12]. Subsequently, sample swabs were placed separately into sterilized tubes and stored immediately at −80 °C for later analysis. Snakes were released back into the wild after sampling.

2.2. DNA Extraction, PCR Amplification, and Purification

DNA was extracted by organic extraction (chloroform: isoamyl alcohol = 96:4) [37], and DNA concentration was determined using a Qubit 4.0 Fluorometer (Thermo Scientific, Waltham, MA, USA) and adjusted to 1 ng/μL for the following process. The full-length 16S genes (V1–V9 regions) were amplified by barcoded 16S gene-specific primers. According to the Amplification of Full-Length 16S Gene with Barcoded Primers for Multiplexed SMRTbell Library Preparation and Sequencing Procedure (PacBio, Menlo Park, CA, USA), each primer was designed to contain a 5′ buffer sequence (GCATC) with a 5′ phosphate modification, a 16-base barcode, and the degenerate 16S gene-specific forward or reverse primer sequences (Forward: 5′Phos/GCATC-16-base barcode-AGRGTTYGATYMTGGCTCAG-3′, Reverse: 5′Phos/GCATC-16-base barcode-RGYTACCTTGTTACGACTT-3′). Degenerate base identities are as follows: R = A, G; Y = C, T; M = A, C. Briefly, 2 ng of gDNA was used for the PCR reaction performed with KAPA HiFi HotStart ReadyMix (Roche, Basel, Switzerland) under the following PCR conditions: 95 °C for 3 min; 20–30 cycles (depending on the sample) at 95 °C for 30 s, 57 °C for 30 s, 72 °C for 60 s; 72 °C for 5 min, and hold at 4 °C. The PCR products were monitored on 1% agarose gel. Samples with bright main strip of approximately 1500 bp were chosen and purified using the AMPure PB beads for the following library preparation.

2.3. SMRTbell Library Construction and Sequencing

The SMRTbell library was prepared according to PacBio. Briefly, equal molars of each barcoded PCR product were pooled, and 500–1000 ng of the pooled amplicon sample was used for DNA damage repair followed by end repair/A-tailing and ligation steps to introduce the universal hairpin adapters onto double-stranded DNA fragments. After purification with AMPure PB beads to remove the adapter dimer, the SMRTbell library was incubated with Sequel II primer 3.1 and Sequel II Binding Kit 3.1 (PacBio, Menlo Park, CA, USA) for the primer annealing and polymerase binding. Finally, sequencing was performed in the circular consensus sequence (CCS) mode on a PacBio Sequel IIe instrument (PacBio, Menlo Park, CA, USA) to generate the HiFi reads with predicted accuracy (Phred Scale) = 30.

2.4. Data Analysis

The CCS reads were determined with a minimum predicted accuracy of 0.9 and the minimum number of passes set to 3 in the official workflow of PacBio through the SMRT Link software (version 11.1). Only CCS reads with a quality score greater than Q30, referred to as Q30 HiFi reads, were used in the next stage. After demultiplexing, the HiFi reads were further processed with DADA2 (version 1.20) to obtain amplicons with single-nucleotide resolution [38,39]. The trimming and filtering were performed with a maximum of two expected errors per read (maxEE = 2). DADA2 algorithm resolves exact amplicon sequence variants (ASVs) with single-nucleotide resolution from the full-length 16S rRNA gene with a near-zero error rate. For each representative sequence, the feature-classifier [40] and classify-consensus-vsearch [41] algorithm in QIIME2 (version 2022.11) [42] was employed to annotate the taxonomy classification based on the information retrieved from the NCBI 16S ribosomal RNA database. To analyze the sequence similarities among different ASVs, multiple sequence alignment was conducted using the QIIME2 alignment, MAFFT [43], against the NCBI 16S ribosomal RNA database [44,45,46]. In addition, the estimation of the abundance was deduced from the assigned read counts of the sequences (=non-chimeric read counts − unassigned read counts); the relative abundance was calculated from the assigned read counts for each snake or bacteria species divided by the total assigned read counts.
Several alpha-diversity indices were used to measure alpha diversity for each sample, including observed-species, Menhinick’s richness, Margalef’s richness, Shannon, Simpson, Pielou’s evenness, PD whole tree, and Good’s coverage [47]. All indices were calculated with QIIME. Kruskal–Wallis tests were conducted to test the differences on the alpha diversity of oral bacteria among snake species. The “ggplot2” package in R (version 3.6.0) was used to generate the plots, including rarefaction curves. Beta diversity analysis was also applied to analyze the difference in species complexity between samples. An UpSet plot was generated using the “UpSetR” package in R to present the number of common and uniquely identified bacterial species in the samples among snake species. Heatmap representation was generated using the “ggplot2” package in R to show the abundance distribution of the top 35 bacterial genera or species identified among the individual oral samples of different snake species. Dissimilarity matrices (Bray-Curtis and weighted Unifrac) [48,49] were generated using the “micro-eco” package in R. To visualize patterns in multidimensional data, the dimensionality reduction method of constrained principal co-ordinates analysis (CPCoA) [50,51] was conducted using the “micro-eco” and “vegan” packages in R. The p-value of CPCoA test was calculated using the “anova.cca” permutation package in R. In addition, we used UPGMA to construct a tree that reflects the similarities between samples presented in the UniFrac distance matrix. “ggplot2” and “plotly” were used to generate the plots. For statistical analysis, the significance of all species among groups at various taxonomic level was detected using differential abundance analysis with a zero-inflated Gaussian (ZIG) log-normal model as implemented in the “fitFeatureModel” function of the Bioconductor “metagenomeSeq” package in R [52]. Moreover, Welch’s t-test was performed using the “stat” package in R to compare metagenomic profiles of the oral bacterial communities among snake species.

3. Results

3.1. Sequence Analysis

Total non-chimeric reads of 71,848 (unassigned reads 670), 54,328 (unassigned reads 0), 68,820 (unassigned reads 35), 80,195 (unassigned reads 40), 61,205 (unassigned reads 105), 55,877 (unassigned reads 145), and 56,836 (unassigned reads 109) were obtained for Na, Bm, Pm, Ts, Ds, Da, and Tg samples, respectively (Table 1). After DADA2 denoising, full-length 16S rRNA sequences for each sample were generated. From a total of 3241 ASVs, 49 were removed because the species was unassigned, and 3192 were considered for further analysis. The rarefaction plot (Supplementary Figure S1) reveals that the number of observed bacterial species plateaus when the sequencing depths reach 10,000 in most samples.

3.2. Taxonomic Profiling of Metagenomic Sequences

We discovered 163, 207, 204, 192, 478, 253, and 471 bacterial species in Na, Bm, Pm, Ts, Ds, Da, and Tg, respectively (Table 2). At the phylum level, Proteobacteria, Bacteroidetes, and Firmicutes were commonly distributed among the samples (Figure 3a), while the Tenericutes, Actinobacteria, Fusobacteria, Cyanobacteria, Verrucomicrobia, and Planctomycetes did not present in some parts of samples (Figure 4a). Proteobacteria (Na: 65.3%, Bm: 62.7%, Pm: 37.6%, Ts: 40.1%, Ds: 59.5%, Da: 70.5%, Tg: 58.9%) and Bacteroidetes (Na: 15.3%, Bm: 18.0%, Pm: 18.9%, Ts: 58.6%, Ds: 15.9%, Da: 19.3%, Tg: 12.3%) were identified as the dominant phyla (Figure 3a). At the species level, Puia dinghuensis, Cupriavidus numazuensis, Mycoplasma fastidiosum, Erysipelatoclostridium innocuum, Bacteroides fragilis, Morganella morganii, Haemophilus felis, Citrobacter freundii, Stenotrophomonas maltophilia, and Bordetella trematum are the 10 most abundant species as a whole, which accounted for 51.71 ± 21.03% of the total (mean ± SD; Figure 3b), while the 20 most abundant species accounted for 67.28 ± 18.31% of the total. None of the 10 most abundant species presented in all individual samples (Figure 4b); the corresponding results for each snake species are shown in Table 3.
Taxonomy annotations and relative abundances of the bacteria are summarized in Supplementary Table S1. A total of 44 unique ASVs with 1104 reads are reported as unassigned. Of the 448,005 total reads, we identified 1064 unique species among the 37 samples. The number of common and unique bacterial species among seven snake species is shown in Figure 5. Seven bacterial species (eight ASVs) were shared by all snake species, including Methylibium petroleiphilum, Delftia lacustris, Caulobacter segnis, Pedobacter nutrimenti, Cupriavidus numazuensis, Sphingomonas paucimobilis, and Brevundimonas diminuta. In addition, 22 species (207 ASVs), 32 species (230 ASVs), 39 species (277 ASVs), 39 species (197 ASVs), 193 species (798 ASVs), 108 species (337 ASVs), and 189 species (591 ASVs) were unique to Na, Bm, Pm, Ts, Ds, Da, and Tg, respectively. The summary of the alpha diversity results is shown in Supplementary Table S2. The alpha diversity of oral bacteria in Ds was the highest (e.g., Kruskal–Wallis test: p = 0.037 for Menhinick’s richness index).
Furthermore, we also presented the abundance distribution of the top 35 bacterial genera and species among individual oral samples of the seven snake species in heatmaps (Figure 6). The abundance patterns of dominant bacterial genera or species varied in the samples among snake species and were more consistent within the same snake species. At the genus level, Bordetella was the most abundant, followed by Pseudomonas and Mycoplasma in Na; Haemophilus, Mycoplasma, and Stenotrophomonas were the most abundant in Bm; Erysipelatoclostridium, Klebsiella, and Puia in Pm; Puia, Citrobacter, and Bacteroides in Ts; Cupriavidus, Romboutsia, and Salmonella in Ds; Morganella, Puia, and Providencia in Da; Cupriavidus, Mycoplasma, and Pseudomonas in Tg (Figure 6a).
Among the top 10 abundant bacterial species in the oral cavity of seven snake species, totally, 24 species have been reported as pathogens or opportunistic pathogens to humans, while 24 species had no pathogenicity toward humans, although they might have been isolated from an infection in a human, and most of them are non-spore-forming bacteria (Table 4). For the top 10 most abundant bacterial species in all snake species, five were pathogenic, four were non-pathogenic, and one remains to be elucidated; one was Gram-positive, and nine were Gram-negative; and six were aerobic (two were facultative anaerobic), and two were anaerobic.

3.3. Comparison of Bacterial Community Structure

The CPCoA analysis revealed that Ts, Ds, Da, and Tg clustered closely by sharing identical ASVs at the species level, whereas the bacterial species identified in Na, Bm, and Pm were uniquely distributed and significantly deviated from Ts, Ds, Da, and Tg (p = 0.001; Figure 7).
Comparison of a metagenomic profile using Welch’s t-test revealed that there are significant differences in abundances of at least one bacterial taxon between Na vs. Bm, Na vs. Pm, Na vs. Ts, Na vs. Ds, Na vs. Da, Na and Tg, Bm vs. Pm, Bm vs. Ts, Bm vs. Ds, Bm vs. Da, Bm vs. Tg, Pm vs. Da, Ts vs. Ds, and Ds vs. Da samples (Figure 8). Among the taxa with a mean proportion >1%, Proteiniphilum acetatigenes was more abundant in Na than in other snake species (all ps < 0.05; Figure 8a–f). Haemophilus felis was more abundant in Bm than in other snake species (all ps = 0.0275; Figure 8a,g–k). Mycoplasma fastidiosum was more abundant in Na than in Pm, Ts, and Ds (all ps < 0.03; Figure 8b–d). Stenotrophomonas maltophilia in Da was less abundant than in Na (p = 0.0345; Figure 8e) or Pm (p = 0.0424; Figure 8l).

4. Discussion

In Taiwan, previous studies on the bacterial species in a snake’s oral cavity were much fewer than those in clinical wound samples. The methods to identify the bacterial species usually included bacterial culture and biochemical testing (CBtest), Sanger sequencing, and Vitek2 system for culturing bacteria [15,33]. However, not all bacteria are easy or able to be cultured, and the bacterial compositions after being cultured may differ from the original ones. In addition, for the 16S rRNA NGS used to investigate the microbiota in a snake’s oral cavity, the precision and accuracy of bacterial identification would be reduced if full-length variable regions are not considered as a sequencing target [34]. Therefore, the bacterial identification results in this study are more robust as we investigated the oral bacterial community composition using full-length 16S rRNA amplicon sequencing analysis.
Our results showed that Proteobacteria and Bacteroidetes were the most dominant phyla among the snakes analyzed. In a previous study, the amplicon 16S rRNA gene V3–V4 hypervariable region showed that Proteobacteria and Actinobacteria are the dominant phyla in the oral bacterial community of Naja naja, Ophiophagus hannah, and Python molurus [36]. Another study using the V4 region of the 16S rRNA gene [53] showed the dominant phyla are Proteobacteria, Firmicutes, and Actinobacteria in saliva samples of Komodo dragons. In the present study, Actinobacteria were not listed in the top three phyla in our results. In addition, previous bacteriomic studies on a snake’s gastrointestinal tract using 16S rRNA gene high-throughput sequencing demonstrated that Proteobacteria and Bacteroidetes are the dominant bacterial phyla in Naja atra, Ptyas mucosa, Elaphe carinata, and Deinagkistrodon acutus [54,55], as well as Rhabdophis subminiatus [56]. Another gut bacteriomic study on sea snakes (Hydrophis curtus and Hydrophis cyanocinctus) using the amplicon V3–V4 regions of the 16S rRNA genes also showed that Proteobacteria and Bacteroidetes were the most abundant phyla [57]. Thus, the dominant bacterial phyla in the present study are similar to those in the studies on bacterial communities in the gastrointestinal tract of reptiles.
The most previous studies in Taiwan collected samples of infected tissue or necrosis sites from the wounds of snake-bite patients and identified the bacterial composition to the genus or species level using the CBtest [33,58]. More research was conducted on common venomous species, like Na, Pm, and Ts, and indicated that Enterococcus faecalis and Morganella morganii are the most dominant bacterial species [59,60,61,62], followed by Aeromonas hydrophila, Serratia marcescens, Proteus vulgaris, Bacteroides fragilis, and Pseudomonas aeruginosa [15,58]. Another study on Pm and Ts showed the dominant bacterial species were Enterococcus faecalis and Morganella morganii, followed by Staphylococcus spp., Corynebacterium spp., and Enterobacter spp. [63]. In addition, a recent study, which sampled the snake’s oral cavity and not bite tissue, used a cultivation and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry analysis, which showed that the most common bacterial species harbored in Na, Bm, Ts, and Pm of southern Taiwan were Enterococcus faecalis, Pseudomonas aeruginosa, Clostridium spp., Bacteroides fragilis, Proteus vulgaris, and Citrobacter freundii [64]. In contrast to the above studies, our results (Table 3) identified different dominant bacteria, which were Bordetella trematum, Mycoplasma fastidiosum, and Phocoenobacter uteri in Na (Morganella morganii ranked 9th); Haemophilus felis, Stenotrophomonas maltophilia, and Puia dinghuensis in Bm; Clostridium innocuum, Cupriavidus numazuensis, and Puia dinghuensis in Pm; and Puia dinghuensis, Bacteroides fragilis, and Citrobacter freundii in Ts (Morganella morganii ranked sixth). Rather, Morganella morganii was the most dominant in Da. Enterococcus faecalis was not in the top 10 species for all snake species but ranked 24th, 50th, 130th, and 200th in Na, Ds, Da, and Tg, respectively.
Furthermore, several reports have presented potential advantages of molecular diagnostics over microbial culture, namely, a shorter turnaround time, detection of bacteria difficult to culture, or detection after prior administration of antibiotics that inhibit bacterial growth [33,65]. The advantage of the NGS-method 16S metagenomics assay is not only the identification of species but also the quantification of the relative abundance of all bacteria in polymicrobial infections [33]. Previous results from a culture-based method in Na only detected 2% of the bacterial species identified using the 16S rRNA NGS method. The most abundant species in Na observed using the 16S rRNA method were different from those determined using culture-based methods in [33], which revealed that the most common bacterial species in the oropharynx of Na were Pseudomonas azotoformans, P. lundensis, Delftia tsuruhatensis, and Methylobacterium goesingense, where only P. azotoformans was also dominant in our study. The variations of bacterial species identified among the different studies may result from the differences in snake species, collection sites (in clinic, out-of-hospital, in animal house, or in wild), sampling tissues (oral cavity, gut, or bite wound), culture dependence, and analysis methods (CBtest, mass spectrometry, or 16S rRNA sequencing).
In our study, Puia dinghuensis was the most abundant bacteria species (ranked in the top 10) among all snake species except Ds (Table 3), but it is not a pathogenic or opportunistic-pathogenic bacterium to humans and is typically isolated from the monsoon evergreen broad-leaved forest soil [66]. Most previous studies focused on pathogenic bacterial species carried in a snake’s oral cavity. Among the seven bacterial species shared by all snake species in our study, only Delftia lacustris and Sphingomonas paucimobilis have been reported as pathogenic, and both are Gram-negative bacteria [67,68,69]. The five remaining bacterial species are all non-pathogenic, including Methylibium petroleiphilum, a Gram-negative bacterium, which prefers aerobic, warm, and close to neutral pH conditions [70,71]; Caulobacter segnis, a Gram-negative bacterium, which is typically isolated from soil [72,73]; Pedobacter nutrimenti, a Gram-negative bacterium and non-spore-forming bacterial strain that is typically isolated from soils, other environmental sources, or from food [74]; Cupriavidus numazuensis, a Gram-negative bacterium typically isolated from soil [75]; Brevundimonas diminuta, a Gram-negative bacterium with a worldwide distribution and is typically isolated from several sites, including water, soil, and plants [76,77].
Generally, bacterial isolates from the oral cavity of snakes are not all pathogenic: They may cause wound infections or abscesses in humans but so do normal environmental contaminants or soil pathogens [33]. Most of pathogenic bacterial species are Gram-negative [7]. Five bacterial species in the top 10 abundant bacterial species found in our study have been reported as pathogenic, including Erysipelatoclostridium innocuum, a Gram-positive bacterium common in intestinal flora [78]; Bacteroides fragilis, a Gram-negative bacterium commonly found as part of the normal microbiota of the human colon; Morganella morganii, a Gram-negative bacterium found in the feces and intestines of humans, dogs, and other mammals [79,80]; Citrobacter freundii, a Gram-negative bacterium found in water, soil, food, and the intestines of humans and other animals [81,82]; Stenotrophomonas maltophilia, a ubiquitous Gram-negative bacterium frequently isolated from soil, water, animals, plant matter, and hospital equipment [83,84]. In addition, the pathogenic potential and infection virulence of Bordetella trematum, commonly isolated from human wounds and ear infections [85], remains to be elucidated [86].
Our results showed that Elapidae snakes, including Na and Bm, had noticeably different oral cavity microbiota from each other and from the species of Viperidae, which may have similar microbiota among each other. It is known that a Na bite often causes severe tissue necrosis while a Bm bite rarely causes a significant local tissue problem [87,88]. Compared to Ts, victims of Pm have a higher rate (37% vs. 6%) of local necrosis or cellulitis [89]. Victims of Ds, Da, and Tg are rare in Taiwan [90,91], and a Ds bite seems to have less complication of local tissue necrosis [91,92]. Tissue necrosis may be initiated by the action of snake toxins [93,94]. How it is related to the differences of oral bacterial strains between Na and Bm or among Ds, Pm, and other pitvipers needs further studies to verify the ability of the microbiota identified in this study that cause infection and necrosis. The microorganisms harbored by reptiles and their abundances could be affected by their environment [53]. The environment in a snake’s oral cavity is likely affected by the snake’s habitat, home range, and predation strategy [7,16]. We also investigated the oral bacteriome of Tg living at a high altitude in Taiwan to test whether different microbiota could be detected by altitudinal factors; however, the oral bacterial community in this species was not significantly different from the low-altitude species on the existence of confounding factors, such as controversial phylogeny [95]. Further studies on the effects of habitat variations on the metagenomics of microorganisms in snakes are required.
There are several limitations in our research. First, there was usually a delay of several days between snake collection and oral sampling, although we sampled the snakes as quickly as possible after wild collection. Second, there were only five or six snake samples per species, which may not properly represent the wild populations. Third, most of our samples were collected from snakes in southern Taiwan, except for Trimeresurus gracilis, which was collected from a high-altitude area around the island; therefore, the collection sites can be expanded in future studies. In addition, although understanding the pathogens from the oral cavities of venomous snakes may assist in empiric antibiotic selection for veterinary medicine and human medicine, physicians still need to consider secondary skin and nosocomial infections after snake envenomation [64]. Our research only focused on venomous snake samples. Most studies focused on venomous snakes because most clinical records are from venomous snake bites. However, pathogens harbored by non-venomous snakes could still be threatening if their teeth are left in the wound or the bite creates a massive open wound. Therefore, assessing the bacterial species carried by non-venomous snakes would fulfill the diversity database of microbiota within a snake’s oral cavity and should be considered in future studies.

5. Conclusions

This study is the first to use full-length 16S rRNA next-generation sequencing (NGS) to investigate the complete oral microbiota carried by seven wild snakes in Taiwan. The results showed similar dominance between non-pathogenic bacteria and pathogenic bacteria, and Gram-negative bacteria were more common than Gram-positive ones. The oral microbiota were significantly different between the snakes belonging to Viperidae and Elapidae and between Naja atra and Bungarus multicinctus within Elapidae, while there was considerable overlap of bacterial species among the species within Viperidae. Our study identified different dominant bacteria compared to previous studies, possibly owing to differences in snake species, collection sites, sampling tissues, culture dependence, and analysis methods. Long-read microbiome analyses of bite wound samples from patients need to be carried out in future studies. This study provides potential therapeutic implications for wildlife medicine and snakebite management.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology12091206/s1, Figure S1: Rarefaction curves showing the relationships between the observed bacterial taxa and sequencing depth for the individual samples from seven Taiwanese snake species; Table S1: Summary of taxonomy annotation and abundances of the bacteria in the oral samples of seven Taiwanese snake species; Table S2: Summary of the alpha diversity results for individual oral samples of seven Taiwanese snake species.

Author Contributions

Conceptualization, T.-S.T.; methodology, W.-H.L. and T.-S.T.; software, W.-H.L. and T.-S.T.; validation, W.-H.L. and T.-S.T.; formal analysis, W.-H.L. and T.-S.T.; investigation, W.-H.L. and T.-S.T.; resources, T.-S.T.; data curation, T.-S.T.; writing—original draft preparation, W.-H.L. and T.-S.T.; writing—review and editing, T.-S.T.; visualization, W.-H.L. and T.-S.T.; supervision, T.-S.T.; project administration, T.-S.T.; funding acquisition, T.-S.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology (grant number 110-2621-B-020-001), the Council of Agriculture, and the Kaohsiung City Government.

Institutional Review Board Statement

The study was approved (approval numbers: NPUST-106-033, Date: 15 August 2019; NPUST-109-023, Date: 24 April 2020; NPUST-109-024, Date: 16 April 2020; NPUST-109-121, Date: 4 February 2021) by the Animal Care and Use Committee of the National Pingtung University of Science and Technology.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article or Supplementary Material.

Acknowledgments

We are grateful to Jun-Wei Zhang, Jui-Hsiang Fan, Tsz-Chun Tse, Fu-Cheng Mi, and Sih-Sheng Cai for helping in the collection of snake specimens and oral samples. We also thank Yuen-Ying Chan for helping in figure preparations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Noutsos, T.; Currie, B.J.; Wijewickrama, E.S.; Isbister, G.K. Snakebite associated thrombotic microangiopathy and recommendations for clinical practice. Toxins 2022, 14, 57. [Google Scholar] [CrossRef]
  2. Chippaux, J.P. Snake-bites: Appraisal of the global situation. Bull. World Health Organ. 1998, 76, 515–524. [Google Scholar]
  3. Sminkey, L. World report on child injury prevention. Inj. Prev. 2008, 14, 69. [Google Scholar] [CrossRef]
  4. Ulima, A.T.; Fawzy, A. Complication in neglected non-venomous snake bite in rural area; what we need to anticipate. Int. J. Med. Sci. Clin. Res. 2022, 2, 1038–1043. [Google Scholar] [CrossRef]
  5. Gutiérrez, J.M.; Theakston, R.D.G.; Warrell, D.A. Confronting the neglected problem of snake bite envenoming: The need for a global partnership. PLoS Med. 2006, 3, e150. [Google Scholar] [CrossRef]
  6. Suchithra, N.; Pappachan, J.; Sujathan, P. Snakebite envenoming in Kerala, South India: Clinical profile and factors involved in adverse outcomes. Emerg. Med. J. 2008, 25, 200–204. [Google Scholar] [CrossRef]
  7. Shaikh, I.K.; Dixit, P.P.; Pawade, B.S.; Potnis-Lele, M.; Kurhe, B.P. Assessment of cultivable oral bacterial flora from important venomous snakes of India and their antibiotic susceptibilities. Curr. Microbiol. 2017, 74, 1278–1286. [Google Scholar] [CrossRef]
  8. Goldstein, E.; Agyare, E.; Vagvolgyi, A.; Halpern, M. Aerobic bacterial oral flora of garter snakes: Development of normal flora and pathogenic potential for snakes and humans. J. Clin. Microbiol. 1981, 13, 954–956. [Google Scholar] [CrossRef]
  9. Padhi, L.; Panda, S.K.; Mohapatra, P.P.; Sahoo, G. Antibiotic susceptibility of cultivable aerobic microbiota from the oral cavity of Echis carinatus from Odisha (India). Microb. Pathog. 2020, 143, 104121. [Google Scholar] [CrossRef]
  10. Gold, B.S.; Dart, R.C.; Barish, R.A. Bites of venomous snakes. N. Engl. J. Med. 2002, 347, 347–356. [Google Scholar] [CrossRef]
  11. Abubakar, S.; Habib, A.; Mathew, J. Amputation and disability following snakebite in Nigeria. Trop. Dr. 2010, 40, 114–116. [Google Scholar] [CrossRef]
  12. Lam, K.K.; Crow, P.; Ng, K.H.; Shek, K.C.; Fung, H.T.; Ades, G.; Grioni, A.; Tan, K.S.; Yip, K.T.; Lung, D.C.; et al. A cross-sectional survey of snake oral bacterial flora from Hong Kong, SAR, China. Emerg. Med. J. 2011, 28, 107–114. [Google Scholar] [CrossRef] [PubMed]
  13. Theakston, R.; Phillips, R.; Looareesuwan, S.; Echeverria, P.; Makin, T.; Warrell, D. Bacteriological studies of the venom and mouth cavities of wild Malayan pit vipers (Calloselasma rhodostoma) in southern Thailand. Trans. R. Soc. Trop. Med. Hyg. 1990, 84, 875–879. [Google Scholar] [CrossRef]
  14. Gilmore, M.S.; Ferretti, J.J. The thin line between gut commensal and pathogen. Science 2003, 299, 1999–2002. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, C.-M.; Wu, K.-G.; Chen, C.-J.; Wang, C.-M. Bacterial infection in association with snakebite: A 10-year experience in a northern Taiwan medical center. J. Microbiol. Immunol. Infect. 2011, 44, 456–460. [Google Scholar] [CrossRef]
  16. Jho, Y.-S.; Park, D.-H.; Lee, J.-H.; Cha, S.-Y.; Han, J.S. Identification of bacteria from the oral cavity and cloaca of snakes imported from Vietnam. Lab. Anim. Res. 2011, 27, 213–217. [Google Scholar] [CrossRef] [PubMed]
  17. Artavia-Leon, A.; Romero-Guerrero, A.; Sancho-Blanco, C.; Rojas, N.; Umana-Castro, R. Diversity of aerobic bacteria isolated from oral and cloacal cavities from free-living snakes species in Costa Rica rainforest. Int. Sch. Res. Notices 2017, 2017, 8934285. [Google Scholar] [CrossRef]
  18. Liu, P.Y.; Weng, L.L.; Tseng, S.Y.; Huang, C.C.; Cheng, C.C.; Mao, Y.C.; Tung, K.C. Colistin resistance of Pseudomonas aeruginosa isolated from snakes in Taiwan. Can. J. Infect. Dis. Med. Microbiol. 2017, 2017, 7058396. [Google Scholar] [CrossRef]
  19. Sadeghi, M.; Barazandeh, M.; Zakariaei, Z.; Davoodi, L.; Tabaripour, R.; Fakhar, M.; Zakariaei, A. Massive cutaneous complications due to snakebite: A case report and literature review. Clin. Case Rep. 2021, 9, e04129. [Google Scholar] [CrossRef]
  20. Blaylock, R. Normal oral bacterial flora from some southern African snakes. Onderstepoort J. Vet. Res. 2001, 68, 175–182. [Google Scholar]
  21. Janda, J.M.; Abbott, S.L. The genus Aeromonas: Taxonomy, pathogenicity, and infection. Clin. Microbiol. Rev. 2010, 23, 35–73. [Google Scholar] [CrossRef] [PubMed]
  22. Mehler, S.J.; Bennett, R.A. Oral, dental, and beak disorders of reptiles. Vet. Clin. N. Am. Exot. Anim. Pract. 2003, 6, 477–503. [Google Scholar] [CrossRef]
  23. Mehler, S.; Bennett, R. Chapter 72—Upper Alimentary Tract Disease. In Reptile Medicine and Surgery, 2nd ed.; Mader, D.R., Ed.; Saunders Elsevier: St. Louis, MI, USA, 2006; pp. 924–930. [Google Scholar]
  24. Hedley, J. Anatomy and disorders of the oral cavity of reptiles and amphibians. Vet. Clin. N. Am. Exot. Anim. Pract. 2016, 19, 689–706. [Google Scholar] [CrossRef]
  25. Schumacher, J. Selected infectious diseases of wild reptiles and amphibians. J. Exot. Pet Med. 2006, 15, 18–24. [Google Scholar] [CrossRef]
  26. Mitchell, M.; Tully, T.N. Manual of Exotic Pet Practice; Elsevier Health Sciences: Saint Louis, MI, USA, 2008; pp. 136–163. [Google Scholar]
  27. Abba, Y.; Ilyasu, Y.M.; Noordin, M.M. Isolation and identification of bacterial populations of zoonotic importance from captive non-venomous snakes in Malaysia. Microb. Pathog. 2017, 108, 49–54. [Google Scholar] [CrossRef]
  28. Blaylock, R. The identification and syndromic management of snakebite in South Africa. S. Afr. Fam. Pract. 2005, 47, 46–53. [Google Scholar] [CrossRef]
  29. Fonseca, M.G.; Moreira, W.; Cunha, K.; Ribeiro, A.; Almeida, M. Oral microbiota of Brazilian captive snakes. J. Venom. Anim. Toxins Incl. Trop. Dis. 2009, 15, 54–60. [Google Scholar] [CrossRef]
  30. Su, H.Y.; Hussain, B.; Hsu, B.M.; Lee, K.H.; Mao, Y.C.; Chiang, L.C.; Chen, J.S. Bacterial community analysis identifies Klebsiella pneumoniae as a native symbiotic bacterium in the newborn Protobothrops mucrosquamatus. BMC Microbiol. 2023, 23, 213. [Google Scholar] [CrossRef]
  31. Greay, T.L.; Gofton, A.W.; Paparini, A.; Ryan, U.M.; Oskam, C.L.; Irwin, P.J. Recent insights into the tick microbiome gained through next-generation sequencing. Parasit. Vectors 2018, 11, 12. [Google Scholar] [CrossRef] [PubMed]
  32. Martinez-Porchas, M.; Villalpando-Canchola, E.; Vargas-Albores, F. Significant loss of sensitivity and specificity in the taxonomic classification occurs when short 16S rRNA gene sequences are used. Heliyon 2016, 2, e00170. [Google Scholar] [CrossRef] [PubMed]
  33. Mao, Y.C.; Chuang, H.N.; Shih, C.H.; Hsieh, H.H.; Jiang, Y.H.; Chiang, L.C.; Lin, W.L.; Hsiao, T.H.; Liu, P.Y. 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, e0009331. [Google Scholar] [CrossRef]
  34. Hu, T.; Chitnis, N.; Monos, D.; Dinh, A. Next-generation sequencing technologies: An overview. Hum. Immunol. 2021, 82, 801–811. [Google Scholar] [CrossRef]
  35. Tse, T.-C.; Tsai, I.-H.; Chan, Y.-Y.; Tsai, T.-S. Venom proteomics of Trimeresurus gracilis, a Taiwan-endemic pitviper, and comparison of its venom proteome and VEGF and CRISP sequences with those of the most related species. Toxins 2023, 15, 408. [Google Scholar] [CrossRef]
  36. Krishnankutty, S.P.; Muraleedharan, M.; Perumal, R.C.; Michael, S.; Benny, J.; Balan, B.; Kumar, P.; Manazhi, J.; Kumar, B.D.; Santhosh, S.; 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. [Google Scholar] [CrossRef]
  37. McKiernan, H.; Danielson, P. Chapter 21—Molecular diagnostic applications in forensic science. In Molecular Diagnostics, 3rd ed.; Patrinos, G.P., Ansorge, W., Danielson, P.B., Eds.; Elsevier: San Diego, CA, USA, 2017; pp. 371–394. [Google Scholar]
  38. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [PubMed]
  39. Quin, C.; Estaki, M.; Vollman, D.; Barnett, J.; Gill, S.; Gibson, D. Probiotic supplementation and associated infant gut microbiome and health: A cautionary retrospective clinical comparison. Sci. Rep. 2018, 8, 8283. [Google Scholar] [CrossRef] [PubMed]
  40. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2′s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef] [PubMed]
  41. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahé, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ 2016, 4, e2584. [Google Scholar] [CrossRef]
  42. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  43. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  44. Balvočiūtė, M.; Huson, D.H. SILVA, RDP, Greengenes, NCBI and OTT—How do these taxonomies compare? BMC Genom. 2017, 18, 114. [Google Scholar] [CrossRef] [PubMed]
  45. Gyarmati, P.; Kjellander, C.; Aust, C.; Song, Y.; Öhrmalm, L.; Giske, C. Metagenomic analysis of bloodstream infections in patients with acute leukemia and therapy-induced neutropenia. Sci. Rep. 2016, 6, 23532. [Google Scholar] [CrossRef] [PubMed]
  46. Hong, X.; Chen, J.; Liu, L.; Wu, H.; Tan, H.; Xie, G.; Xu, Q.; Zou, H.; Yu, W.; Wang, L. Metagenomic sequencing reveals the relationship between microbiota composition and quality of Chinese Rice Wine. Sci. Rep. 2016, 6, 26621. [Google Scholar] [CrossRef] [PubMed]
  47. Whittaker, R.H. Evolution and measurement of species diversity. Taxon 1972, 21, 213–251. [Google Scholar] [CrossRef]
  48. Lozupone, C.; Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Appl. Environ. Microbiol. 2005, 71, 8228–8235. [Google Scholar] [CrossRef] [PubMed]
  49. Lozupone, C.; Lladser, M.E.; Knights, D.; Stombaugh, J.; Knight, R. UniFrac: An effective distance metric for microbial community comparison. ISME J. 2011, 5, 169–172. [Google Scholar] [CrossRef] [PubMed]
  50. Legendre, P.; Anderson, M.J. Distance-based redundancy analysis: Testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 1999, 69, 1–24. [Google Scholar] [CrossRef]
  51. Zgadzaj, R.; Garrido-Oter, R.; Jensen, D.B.; Koprivova, A.; Schulze-Lefert, P.; Radutoiu, S. Root nodule symbiosis in Lotus japonicus drives the establishment of distinctive rhizosphere, root, and nodule bacterial communities. Proc. Natl. Acad. Sci. USA 2016, 113, E7996–E8005. [Google Scholar] [CrossRef]
  52. Paulson, J.N.; Stine, O.C.; Bravo, H.C.; Pop, M. Differential abundance analysis for microbial marker-gene surveys. Nat. Methods 2013, 10, 1200–1202. [Google Scholar] [CrossRef]
  53. Hyde, E.R.; Navas-Molina, J.A.; Song, S.J.; Kueneman, J.G.; Ackermann, G.; Cardona, C.; Humphrey, G.; Boyer, D.; Weaver, T.; Mendelson, J.R., 3rd; et al. The oral and skin microbiomes of captive Komodo dragons are significantly shared with their habitat. mSystems 2016, 1, e00046-16. [Google Scholar] [CrossRef]
  54. Qin, Z.; Wang, S.; Guo, D.; Zhu, J.; Chen, H.; Bai, L.; Luo, X.; Yin, Y. Comparative analysis of intestinal bacteria among venom secretion and non-secrection snakes. Sci. Rep. 2019, 9, 6335. [Google Scholar] [CrossRef] [PubMed]
  55. Zhang, B.; Ren, J.; Yang, D.; Liu, S.; Gong, X. Comparative analysis and characterization of the gut microbiota of four farmed snakes from southern China. PeerJ 2019, 7, e6658. [Google Scholar] [CrossRef] [PubMed]
  56. Tang, W.; Zhu, G.; Shi, Q.; Yang, S.; Ma, T.; Mishra, S.K.; Wen, A.; Xu, H.; Wang, Q.; Jiang, Y.; et al. Characterizing the microbiota in gastrointestinal tract segments of Rhabdophis subminiatus: Dynamic changes and functional predictions. Microbiologyopen 2019, 8, e00789. [Google Scholar] [CrossRef]
  57. Zhong, J.; Guo, K.; Liao, Z.-L.; Hu, S.-C.; Du, Y.; Ji, X. Comparative analysis of the gut microbiota composition between two sea snakes, Hydrophis curtus, and Hydrophis cyanocinctus. Coral Reefs 2022, 41, 53–62. [Google Scholar] [CrossRef]
  58. Huang, L.-W.; Wang, J.-D.; Huang, J.-A.; Hu, S.-Y.; Wang, L.-M.; Tsan, Y.-T. Wound infections secondary to snakebite in central Taiwan. J. Venom. Anim. Toxins Incl. Trop. Dis. 2011, 18, 272–276. [Google Scholar] [CrossRef]
  59. Yeh, H.; Gao, S.-Y.; Lin, C.-C. Wound infections from Taiwan Cobra (Naja atra) bites: Determining bacteriology, antibiotic susceptibility, and the use of antibiotics-a cobra BITE study. Toxins 2021, 13, 183. [Google Scholar] [CrossRef]
  60. Mao, Y.-C.; Liu, P.-Y.; Hung, D.-Z.; Lai, W.-C.; Huang, S.-T.; Hung, Y.-M.; Yang, C.-C. Bacteriology of Naja atra snakebite wound and its implications for antibiotic therapy. Am. J. Trop. Med. Hyg. 2016, 94, 1129. [Google Scholar] [CrossRef]
  61. Mao, Y.-C.; Liu, P.-Y.; Chiang, L.-C.; Lai, C.-S.; Lai, K.-L.; Ho, C.-H.; Wang, T.-H.; Yang, C.-C. Naja atra snakebite in Taiwan. Clin. Toxicol. 2018, 56, 273–280. [Google Scholar] [CrossRef]
  62. Yeh, H.; Gao, S.Y.; Lin, C.C. Wound infection of snakebite from venomous Protobothrops mucrosquamatus, Viridovipera stejnegeri and Naja atra in Taiwan: Validation of BITE and cobra BITE scoring systems and their bacteriological differences in wound cultures. Toxins 2023, 15, 78. [Google Scholar] [CrossRef]
  63. Lin, C.-C.; Chen, Y.-C.; Goh, Z.N.L.; Seak, C.-K.; Seak, J.C.-Y.; Shi-Ying, G.; Seak, C.-J.; Investigators, S. 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 2020, 12, 575. [Google Scholar] [CrossRef]
  64. 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, 951. [Google Scholar] [CrossRef] [PubMed]
  65. Rampini, S.K.; Bloemberg, G.V.; Keller, P.M.; Büchler, A.C.; Dollenmaier, G.; Speck, R.F.; Böttger, E.C. Broad-range 16S rRNA gene polymerase chain reaction for diagnosis of culture-negative bacterial infections. Clin. Infect. Dis. 2011, 53, 1245–1251. [Google Scholar] [CrossRef]
  66. Lv, Y.Y.; Gao, Z.H.; Xia, F.; Chen, M.H.; Qiu, L.H. Puia dinghuensis gen. nov., sp. nov., isolated from monsoon evergreen broad-leaved forest soil. Int. J. Syst. Evol. Microbiol. 2017, 67, 4639–4645. [Google Scholar] [CrossRef]
  67. Shin, S.; Choi, J.; Ko, K. Four cases of possible human infections with Delftia lacustris. Infection 2012, 40, 709–712. [Google Scholar] [CrossRef] [PubMed]
  68. Sohn, K.M.; Baek, J.-Y.; Cheon, S.; Kim, Y.-S.; Koo, S.H. Ocular infection associated with Delftia lacustris: First report. Braz. J. Infect. Dis. 2015, 19, 449–450. [Google Scholar] [CrossRef]
  69. Jørgensen, N.O.; Brandt, K.K.; Nybroe, O.; Hansen, M. Delftia lacustris sp. nov., a peptidoglycan-degrading bacterium from fresh water, and emended description of Delftia tsuruhatensis as a peptidoglycan-degrading bacterium. Int. J. Syst. Evol. Microbiol. 2009, 59, 2195–2199. [Google Scholar] [CrossRef]
  70. Kane, S.R.; Chakicherla, A.Y.; Chain, P.S.; Schmidt, R.; Shin, M.W.; Legler, T.C.; Scow, K.M.; Larimer, F.W.; Lucas, S.M.; Richardson, P.M. Whole-genome analysis of the methyl tert-butyl ether-degrading beta-proteobacterium Methylibium petroleiphilum PM1. J. Bacteriol. 2007, 189, 1931–1945. [Google Scholar] [CrossRef]
  71. Nakatsu, C.H.; Hristova, K.; Hanada, S.; Meng, X.-Y.; Hanson, J.R.; Scow, K.M.; Kamagata, Y. Methylibium petroleiphilum gen. nov., sp. nov., a novel methyl tert-butyl ether-degrading methylotroph of the Betaproteobacteria. Int. J. Syst. Evol. Microbiol. 2006, 56, 983–989. [Google Scholar] [CrossRef]
  72. Moore, G.M.; Gitai, Z. Both clinical and environmental Caulobacter species are virulent in the Galleria mellonella infection model. PLoS ONE 2020, 15, e0230006. [Google Scholar] [CrossRef] [PubMed]
  73. Patel, S.; Fletcher, B.; Scott, D.C.; Ely, B. Genome sequence and phenotypic characterization of Caulobacter segnis. Curr. Microbiol. 2015, 70, 355–363. [Google Scholar] [CrossRef]
  74. Derichs, J.; Kämpfer, P.; Lipski, A. Pedobacter nutrimenti sp. nov., isolated from chilled food. Int. J. Syst. Evol. Microbiol. 2014, 64, 1310–1316. [Google Scholar] [CrossRef] [PubMed]
  75. Martínez-Aguilar, L.; Caballero-Mellado, J.; Estrada-De Los Santos, P. Transfer of Wautersia numazuensis to the genus Cupriavidus as Cupriavidus numazuensis comb. nov. Int. J. Syst. Evol. Microbiol. 2013, 63, 208–211. [Google Scholar] [CrossRef] [PubMed]
  76. Almuzara, M.N.; Barberis, C.M.; Rodríguez, C.H.; Famiglietti, A.M.; Ramirez, M.S.; Vay, C.A. First report of an extensively drug-resistant VIM-2 metallo-β-lactamase-producing Brevundimonas diminuta clinical isolate. J. Clin. Microbiol. 2012, 50, 2830–2832. [Google Scholar] [CrossRef] [PubMed]
  77. Lupande-Mwenebitu, D.; Tshiyongo, R.K.; Lunguya-Metila, O.; Lavigne, J.-P.; Rolain, J.-M.; Diene, S.M. First isolation and clinical case of Brevundimonas diminuta in a newborn with low birth weight, in Democratic Republic of Congo: A case report. Medicina 2021, 57, 1227. [Google Scholar] [CrossRef]
  78. Cherny, K.E.; Muscat, E.B.; Reyna, M.E.; Kociolek, L.K. Clostridium innocuum: Microbiological and clinical characteristics of a potential emerging pathogen. Anaerobe 2021, 71, 102418. [Google Scholar] [CrossRef]
  79. Kim, J.H.; Cho, C.R.; Um, T.H.; Rhu, J.Y.; Kim, E.S.; Jeong, J.W.; Lee, H.R. Morganella morganii sepsis with massive hemolysis. J. Korean Med. Sci. 2007, 22, 1082–1084. [Google Scholar] [CrossRef]
  80. Liu, H.; Zhu, J.; Hu, Q.; Rao, X. Morganella morganii, a non-negligent opportunistic pathogen. Int. J. Infect. Dis. 2016, 50, 10–17. [Google Scholar] [CrossRef] [PubMed]
  81. Wang, J.-T.; Chang, S.-C.; Chen, Y.-C.; Luh, K.-T. Comparison of antimicrobial susceptibility of Citrobacter freundii isolates in two different time periods. J. Microbiol. Immunol. Infect. 2000, 33, 258–262. [Google Scholar]
  82. Liu, H.; Zhao, Z.; Xue, Y.; Ding, K.; Xue, Q. Fatal cases of Citrobacter freundii septicemia and encephalitis in sheep. J. Vet. Diagn. Invest. 2018, 30, 245–248. [Google Scholar] [CrossRef]
  83. Denton, M.; Kerr, K.G. Microbiological and clinical aspects of infection associated with Stenotrophomonas maltophilia. Clin. Microbiol. Rev. 1998, 11, 57–80. [Google Scholar] [CrossRef]
  84. Garcia Paez, J.; Tengan, F.; Barone, A.; Levin, A.; Costa, S. Factors associated with mortality in patients with bloodstream infection and pneumonia due to Stenotrophomonas maltophilia. Eur. J. Clin. Microbiol. 2008, 27, 901–906. [Google Scholar] [CrossRef]
  85. Vandamme, P.; Heyndrickx, M.; Vancanneyt, M.; Hoste, B.; De Vos, P.; Falsen, E.; Kersters, K.; Hinz, K.-H. Bordetella trematum sp. nov., isolated from wounds and ear infections in humans, and reassessment of Alcaligenes denitrificans Rüger and Tan 1983. Int. J. Syst. Evol. Microbiol. 1996, 46, 849–858. [Google Scholar] [CrossRef] [PubMed]
  86. Kukla, R.; Svarc, M.; Bolehovska, R.; Ryskova, L.; Paterova, P.; Fajfr, M.; Malisova, L.; Zemlickova, H. Isolation of Bordetella trematum from the respiratory tract of a patient with lung cancer: A case report. Folia Microbiol. 2020, 65, 623–627. [Google Scholar] [CrossRef]
  87. Mao, Y.C.; Liu, P.Y.; Chiang, L.C.; Liao, S.C.; Su, H.Y.; Hsieh, S.Y.; Yang, C.C. Bungarus multicinctus multicinctus snakebite in Taiwan. Am. J. Trop. Med. Hyg. 2017, 96, 1497–1504. [Google Scholar] [CrossRef]
  88. Lai, C.S.; Liu, P.Y.; Lee, C.H.; Ho, C.H.; Chen, W.L.; Lai, K.L.; Su, H.Y.; Lin, W.L.; Chung, K.C.; Yang, Y.Y.; et al. The development of surgical risk score and evaluation of necrotizing soft tissue infection in 161 Naja atra envenomed patients. PLoS Negl. Trop. Dis. 2022, 16, e0010066. [Google Scholar] [CrossRef]
  89. Chen, Y.-W.; Chen, M.-H.; Chen, C.-K.; Lee, C.-H.; Hung, D.-Z.; Huang, C.-I.; Yang, C.-C.; Chen, Y.-C.; Yen, D.H.-T.; Wang, L.-M. Differences in clinical profiles of patients with Protobothrops mucrosquamatus and Viridovipera stejnegeri envenoming in Taiwan. Am. J. Trop. Med. Hyg. 2009, 80, 28–32. [Google Scholar] [CrossRef] [PubMed]
  90. Tsai, T.S.; Chan, Y.Y.; Huang, S.M.; Chuang, P.C. Case report: Symptoms and prognosis of Trimeresurus gracilis envenomation. Am. J. Trop. Med. Hyg. 2022, 106, 1281–1284. [Google Scholar] [CrossRef] [PubMed]
  91. Su, H.Y.; Huang, S.W.; Mao, Y.C.; Liu, M.W.; Lee, K.H.; Lai, P.F.; Tsai, M.J. Clinical and laboratory features distinguishing between Deinagkistrodon acutus and Daboia siamensis envenomation. J. Venom. Anim. Toxins Incl. Trop. Dis. 2018, 24, 43. [Google Scholar] [CrossRef]
  92. Tsai, T.S.; Liu, C.C.; Chuang, P.C. Personal experience of Daboia siamensis envenomation. Case Rep. Med. 2021, 2021, 3396373. [Google Scholar] [CrossRef]
  93. Ho, C.H.; Chiang, L.C.; Mao, Y.C.; Lan, K.C.; Tsai, S.H.; Shih, Y.J.; Tzeng, Y.S.; Lin, C.S.; Lin, W.L.; Fang, W.H.; et al. Analysis of the necrosis-inducing components of the venom of Naja atra and assessment of the neutralization ability of freeze-dried antivenom. Toxins 2021, 13, 619. [Google Scholar] [CrossRef]
  94. Chiang, L.C.; Chien, K.Y.; Su, H.Y.; Chen, Y.C.; Mao, Y.C.; Wu, W.G. Comparison of protein variation in Protobothrops mucrosquamatus venom between northern and southeast Taiwan and association with human envenoming effects. Toxins 2022, 14, 643. [Google Scholar] [CrossRef] [PubMed]
  95. Zaher, H.; Murphy, R.W.; Arredondo, J.C.; Graboski, R.; Machado-Filho, P.R.; Mahlow, K.; Montingelli, G.G.; Quadros, A.B.; Orlov, N.L.; Wilkinson, M. Large-scale molecular phylogeny, morphology, divergence-time estimation, and the fossil record of advanced caenophidian snakes (Squamata: Serpentes). PLoS ONE 2019, 14, e0216148. [Google Scholar] [CrossRef]
Biology 12 01206 g001
Figure 1. Venomous snakes used for oral sample collection: (a) Naja atra; (b) Bungarus multicinctus; (c) Deinagkistrodon acutus; (d) Daboia siamensis; (e) Protobothrops mucrosquamatus; (f) Trimeresurus stejnegeri; (g) Trimeresurus gracilis; (h) oral sample collection from P. mucrosquamatus using sterilized cotton-tipped swab and sterilized, hollowed out wooden stir stick. Photos courtesy of Jun-Wei Zhang (a,b,e), Jui-Hsiang Fan (c,f), and Tsz-Chun Tse (d,g).
Figure 1. Venomous snakes used for oral sample collection: (a) Naja atra; (b) Bungarus multicinctus; (c) Deinagkistrodon acutus; (d) Daboia siamensis; (e) Protobothrops mucrosquamatus; (f) Trimeresurus stejnegeri; (g) Trimeresurus gracilis; (h) oral sample collection from P. mucrosquamatus using sterilized cotton-tipped swab and sterilized, hollowed out wooden stir stick. Photos courtesy of Jun-Wei Zhang (a,b,e), Jui-Hsiang Fan (c,f), and Tsz-Chun Tse (d,g).
Biology 12 01206 g001
Biology 12 01206 g002
Figure 2. The collection sites and sample sizes of two elapid and five viperid species in five administrative districts of Taiwan.
Figure 2. The collection sites and sample sizes of two elapid and five viperid species in five administrative districts of Taiwan.
Biology 12 01206 g002
Biology 12 01206 g003
Figure 3. Relative abundances of the dominant bacterial phyla (a) and species (b) in the oral samples of seven Taiwanese snake species. Those with sequence reads that did not associate with any known reference taxon or the relative abundance is lower than that of the top 10 species were classified as Others. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Figure 3. Relative abundances of the dominant bacterial phyla (a) and species (b) in the oral samples of seven Taiwanese snake species. Those with sequence reads that did not associate with any known reference taxon or the relative abundance is lower than that of the top 10 species were classified as Others. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Biology 12 01206 g003
Biology 12 01206 g004
Figure 4. Relative abundances of the dominant bacterial phyla (a) and species (b) in the individual oral samples of seven Taiwanese snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Figure 4. Relative abundances of the dominant bacterial phyla (a) and species (b) in the individual oral samples of seven Taiwanese snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Biology 12 01206 g004
Biology 12 01206 g005
Figure 5. The UpSet plot showing the number of common and uniquely identified bacterial species among the oral samples of seven Taiwanese snake species. For example, for all seven snake species, only seven bacterial species (eight ASVs) were common in all 1064 bacterial species found. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Figure 5. The UpSet plot showing the number of common and uniquely identified bacterial species among the oral samples of seven Taiwanese snake species. For example, for all seven snake species, only seven bacterial species (eight ASVs) were common in all 1064 bacterial species found. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Biology 12 01206 g005
Biology 12 01206 g006
Figure 6. Heatmap showing of the abundance distribution of the top 35 bacterial genera and species identified among the individual oral samples of seven Taiwanese snake species: (a) The distribution of the top 35 abundant bacterial genera and (b) the top 35 abundant bacterial species among individual snake samples. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Figure 6. Heatmap showing of the abundance distribution of the top 35 bacterial genera and species identified among the individual oral samples of seven Taiwanese snake species: (a) The distribution of the top 35 abundant bacterial genera and (b) the top 35 abundant bacterial species among individual snake samples. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Biology 12 01206 g006
Biology 12 01206 g007
Figure 7. The coordinate plot for the constrained principal coordinates analysis (CPCoA) of the oral bacterial communities among the samples of seven Taiwanese snake species. The percentages in the brackets represent the contribution to sample variation. The significance of the CPCoA calculated from the “anova.cca” processed permutation in R was p = 0.001. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Figure 7. The coordinate plot for the constrained principal coordinates analysis (CPCoA) of the oral bacterial communities among the samples of seven Taiwanese snake species. The percentages in the brackets represent the contribution to sample variation. The significance of the CPCoA calculated from the “anova.cca” processed permutation in R was p = 0.001. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Biology 12 01206 g007
Biology 12 01206 g008
Figure 8. Metagenomic profile comparisons of the oral bacterial communities among the seven Taiwanese snake species, determined using Welch’s t-test. Only the taxa abundance comparisons with a significant difference (p < 0.05) are shown in (a) Na vs. Bm, (b) Na vs. Pm, (c) Na vs. Ts, (d) Na vs. Ds, (e) Na vs. Da, (f) Na and Tg, (g) Bm vs. Pm, (h) Bm vs. Ts, (i) Bm vs. Ds, (j) Bm vs. Da, (k) Bm vs. Tg, (l) Pm vs. Da, (m) Ts vs. Ds, and (n) Ds vs. Da samples. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Figure 8. Metagenomic profile comparisons of the oral bacterial communities among the seven Taiwanese snake species, determined using Welch’s t-test. Only the taxa abundance comparisons with a significant difference (p < 0.05) are shown in (a) Na vs. Bm, (b) Na vs. Pm, (c) Na vs. Ts, (d) Na vs. Ds, (e) Na vs. Da, (f) Na and Tg, (g) Bm vs. Pm, (h) Bm vs. Ts, (i) Bm vs. Ds, (j) Bm vs. Da, (k) Bm vs. Tg, (l) Pm vs. Da, (m) Ts vs. Ds, and (n) Ds vs. Da samples. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Biology 12 01206 g008
Table 1. Data summary of sequence reads for the individual oral samples of seven Taiwanese snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Table 1. Data summary of sequence reads for the individual oral samples of seven Taiwanese snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Sample NameRaw HiFi ReadsRemove PrimersFiltered ReadsDenoised ReadsNon-Chimeric Reads
Na 120,24119,17918,34517,94017,910
Na 212,77911,29310,98310,85010,850
Na 313,11211,00310,55910,2859955
Na 411,74011,07010,80010,70910,709
Na 511,53510,93610,70710,60710,607
Na 613,05712,33311,97511,81711,817
Bm 113,65212,61112,22611,95111,940
Bm 210,80610,38810,13510,02910,029
Bm 312,30611,21310,88210,60710,607
Bm 415,04313,33713,04812,91512,915
Bm 510,4369241899088378837
Pm 115,81314,10813,73213,69013,690
Pm 214,11312,00011,40511,25211,039
Pm 312,35010,178980497639763
Pm 413,03912,31511,89311,64311,643
Pm 513,77111,94111,61111,48811,488
Pm 613,37211,61011,32711,19711,197
Ts 124,25623,65922,68422,58922,589
Ts 213,49012,72312,34912,26912,269
Ts 315,99915,14714,73914,60514,605
Ts 412,48411,83911,56811,50711,507
Ts 521,20519,91219,37919,22519,225
Ds 113,79312,31811,98311,85411,852
Ds 215,17613,64913,22913,17513,175
Ds 313,42213,16612,73912,43312,298
Ds 414,30813,08412,73412,55812,546
Ds 512,14111,74011,42911,33411,334
Da 113,65112,42611,96411,81911,816
Da 213,39513,03212,65712,61912,606
Da 313,20012,48912,05811,65511,655
Da 411,58310,217999699109910
Da 511,51710,151994998909890
Tg 113,63712,41912,04411,88611,886
Tg 212,07110,96710,69010,53510,535
Tg 312,85011,45311,14910,82710,827
Tg 413,42411,99311,71211,54911,549
Tg 513,91312,60912,32012,03912,039
Table 2. Numbers of discovered non-chimeric reads and taxa for the oral samples from seven Taiwanese snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Table 2. Numbers of discovered non-chimeric reads and taxa for the oral samples from seven Taiwanese snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
NaBmPmTsDsDaTg
Reads71,84854,32868,82080,19561,20555,87756,836
Phylum88678913
Class22202121232430
Order39373742454555
Family697072699581103
Genus115115124121226155246
Species163207204192478253471
Table 3. Characteristics (highlighted in blue: aerobic Gram-positive bacteria; red: aerobic Gram-negative bacteria; green: anaerobic bacteria) and relative abundances (in parentheses) of the top ten bacterial taxa present in oral cavities of each snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
Table 3. Characteristics (highlighted in blue: aerobic Gram-positive bacteria; red: aerobic Gram-negative bacteria; green: anaerobic bacteria) and relative abundances (in parentheses) of the top ten bacterial taxa present in oral cavities of each snake species. Na: Naja atra, Bm: Bungarus multicinctus, Pm: Protobothrops mucrosquamatus, Ts: Trimeresurus stejnegeri, Ds: Daboia siamensis, Da: Deinagkistrodon acutus, and Tg: Trimeresurus gracilis.
OrderNaBmPmTsDsDaTg
1Bordetella
trematum
(16.11%) 1		Haemophilus
felis
(29.60%)		Clostridium
innocuum
(31.73%)		Puia
dinghuensis
(36.04%)		Clostridium
dakarense
(10.71%)		Morganella
morganii
(29.12%)		Cupriavidus
numazuensis
(22.37%)
2Mycoplasma
fastidiosum
(12.43%)		Stenotrophomonas
maltophilia
(15.23%)		Cupriavidus
numazuensis
(11.77%)		Bacteroides
fragilis
(21.40%)		Cupriavidus
numazuensis
(9.47%)		Providencia
rustigianii
(14.96%)		Mycoplasma
fastidiosum
(14.90%)
3Phocoenobacter
uteri
(12.21%)		Puia
dinghuensis
(10.74%)		Puia
dinghuensis
(11.58%)		Citrobacter
freundii
(17.14%)		Salmonella
enterica enterica
serovar
Typhimurium
(4.63%)		Puia
dinghuensis
(12.63%)		Pseudomonas
aeruginosa
(9.39%)
4Pseudomonas
helleri
(11.83%)		Mycoplasma
fastidiosum
(7.93%)		Klebsiella
pneumoniae
(9.89%)		Yokenella
regensburgei
(6.13%)		Sphingomonas
paucimobilis
(4.34%)		Mycoplasma
fastidiosum
(6.68%)		Puia
dinghuensis
(6.69%)
5Bisgaardia
hudsonensis
(11.21%)		Mycoplasma
iguana
(4.82%)		Clostridium
cylindrosporum
(8.14%)		Stenotrophomonas
maltophilia
(5.46%)		Bacteroides
fragilis
(4.20%)		Rheinheimera
pacifica
(5.33%)		Sphingomonas
paucimobilis
(4.50%)
6Puia
dinghuensis
(6.58%)		Bergeyella
zoohelcum
(4.56%)		Sphingobacterium
detergens
(2.88%)		Morganella
morganii
(2.29%)		Paracoccus
chinensis
(3.90%)		Cupriavidus
numazuensis
(4.14%)		Sediminibacterium
lactis
(2.45%)
7Proteiniphilum
acetatigenes
(4.63%)		Achromobacter
insuavis
(4.00%)		Luteibacter
anthropic
(2.39%)		Phyllobacterium
myrsinacearum
(1.86%)		Stenotrophomonas
pavanii
(3.86%)		Clostridium
innocuum
(2.09%)		Latilactobacillus
curvatus
(1.50%)
8Serratia
marcescens
(2.33%)		Acinetobacter
bereziniae
(3.13%)		Desulfovibrio
desulfuricans
(2.18%)		Beijerinckia
fluminensis
(1.71%)		Chryseobacterium
salipaludis
(3.24%)		Rheinheimera
hassiensis
(1.90%)		Comamonas
testosterone
(1.45%)
9Morganella
morganii
(1.74%)		Clostridium
cylindrosporum
(2.78%)		Epilithonimonas
lactis
(2.01%)		Cupriavidus
numazuensis
(0.80%)		Mammaliicoccus
sciuri
(2.70%)		Acinetobacter
proteolyticus
(1.85%)		Paracoccus
suum
(1.35%)
10Haoranjiania
flava
(1.64%)		Comamonas
testosterone
(2.12%)		Oceanisphaera
ostreae
(0.94%)		Tepidimonas
aquatica
(0.77%)		Acinetobacter
variabilis
(2.35%)		Rheinheimera
aquimaris
(1.61%)		Paenarthrobacter
nitroguajacolicus
(1.35%)
1 The relative abundance was calculated from the assigned read counts for each bacteria species divided by the total assigned read counts for a single snake species.
Table 4. The 24 pathogenic (to humans) bacterial species and 24 non-pathogenic bacterial species among the top 10 abundant bacterial species of seven Taiwanese snake species. G+: Gram-positive, G−: Gram-negative, An: anaerobic bacteria.
Table 4. The 24 pathogenic (to humans) bacterial species and 24 non-pathogenic bacterial species among the top 10 abundant bacterial species of seven Taiwanese snake species. G+: Gram-positive, G−: Gram-negative, An: anaerobic bacteria.
Pathogen Bacterial SpeciesNon-Pathogen Bacterial Species
Achromobacter insuavis (G−)Beijerinckia fluminensis (G−)
Acinetobacter bereziniae (G−)Chryseobacterium salipaludis (G−)
Acinetobacter proteolyticus (G−)Clostridium cylindrosporum (G+, An)
Acinetobacter variabilis (G−)Clostridium dakarense (G+, An)
Bacteroides fragilis (G−, An) 1Cupriavidus numazuensis (G−) 1
Bergeyella zoohelcum (G−)Epilithonimonas lactis (G−)
Bisgaardia hudsonensis (G−)Haemophilus felis (G−) 1
Citrobacter freundii (G−) 1Haoranjiania flava (G−)
Clostridium innocuum (G+, An) 1Latilactobacillus curvatus (G+, An)
Comamonas testosterone (G−)Mycoplasma fastidiosum (G−) 1
Desulfovibrio desulfuricans (G−, An)Mycoplasma iguana (G−)
Klebsiella pneumoniae (G−)Oceanisphaera ostreae (G−)
Morganella morganii (G−) 1Paracoccus chinensis (G−)
Pseudomonas aeruginosa (G−)Phocoenobacter uter (G−)
Serratia marcescens (G−)Proteiniphilum acetatigenes (G−, An)
Luteibacter anthropi (G−)Pseudomonas helleri (G−)
Mammaliicoccus sciuri (G+)Puia dinghuensis (G−) 1
Paenarthrobacter nitroguajacolicus (G+)Rheinheimera aquimaris (G−)
Phyllobacterium myrsinacearum (G−)Rheinheimera hassiensis (G−)
Providencia rustigianii (G−)Rheinheimera pacifica (G−)
Salmonella enterica enterica serovar Typhimurium (G−)Sediminibacterium lactis (G−)
Sphingomonas paucimobilis (G−)Sphingobacterium detergens (G−)
Stenotrophomonas maltophilia (G−) 1Stenotrophomonas pavanii (G−)
Yokenella regensburgei (G−)Tepidimonas aquatica (G−)
1 This species belongs to the top 10 abundant bacterial species in all snake species.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lin, W.-H.; Tsai, T.-S. Comparisons of the Oral Microbiota from Seven Species of Wild Venomous Snakes in Taiwan Using the High-Throughput Amplicon Sequencing of the Full-Length 16S rRNA Gene. Biology 2023, 12, 1206. https://doi.org/10.3390/biology12091206

AMA Style

Lin W-H, Tsai T-S. Comparisons of the Oral Microbiota from Seven Species of Wild Venomous Snakes in Taiwan Using the High-Throughput Amplicon Sequencing of the Full-Length 16S rRNA Gene. Biology. 2023; 12(9):1206. https://doi.org/10.3390/biology12091206

Chicago/Turabian Style

Lin, Wen-Hao, and Tein-Shun Tsai. 2023. "Comparisons of the Oral Microbiota from Seven Species of Wild Venomous Snakes in Taiwan Using the High-Throughput Amplicon Sequencing of the Full-Length 16S rRNA Gene" Biology 12, no. 9: 1206. https://doi.org/10.3390/biology12091206

APA Style

Lin, W. -H., & Tsai, T. -S. (2023). Comparisons of the Oral Microbiota from Seven Species of Wild Venomous Snakes in Taiwan Using the High-Throughput Amplicon Sequencing of the Full-Length 16S rRNA Gene. Biology, 12(9), 1206. https://doi.org/10.3390/biology12091206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop