Abstract
To further investigate the bacterial community and identify the bacterial biomarkers between venom secretion and non-venom secretion snakes, 50 intestinal samples (25 large intestine, 25 small intestine) were obtained from 29 snakes (13 gut samples from Deinagkistrodon, 26 from Naja and 11 from Ptyas mucosa). 16S rDNA high-throughput sequencing results showed that 29 bacterial phyla, 545 bacterial genera, and 1,725 OTUs (operational taxonomic units) were identified in these samples. OTU numbers and the Ace, Chao, Shannon, and Simpson indexes were very similar among the three breeds of snakes included in this study. The Bacteroidetes, Firmicutes, Fusobacteria and Proteobacteria were predominant bacterial phyla. The relative abundance at the phylum level among these samples was similar, and the difference between small and large intestinal samples was not obvious. However, at the genus level, venom secretion snakes Deinagkistrodon and Naja clustered together according to different breeds. 27, 24, and 16 genera were identified as core microbes for Deinagkistrodon, Naja, and Ptyas mucosa, respectively. Interestingly, the relative abundances of genera Hafnia_Obesumbacterium, Providencia, and Ureaplasma were found to be significantly higher in non-venom secretion snakes, and the genera Achromobacter, Cetobacterium, Clostridium innocuum group, Fusobacterium, Lachnoclostridium, Parabacteroides, and Romboutsia were only detected in venom secretion snakes. The function of these bacteria in venom secretion needs to be further studied, and these venom secretion related genera may be the promising target to improve venom production.
Similar content being viewed by others
Introduction
Gut microbiota refer to a large number of microorganisms existing in the intestine of human and animals. Gut microbitoa take part in host metabolization and regulate host immune systems development, which plays a very important role in maintaining the host health balance1,2. Many diseases have been reported to be associated with changes in diversity or richness of gut microbiota3. So far, most of the studies have focused on mammalian vertebrates’ gut microbiota, with fewer investigations on reptile vertebrates such as snake. As one of the special economic animals, snake has high commercial values as its whole body has value. Snake gall is regarded as a precious Chinese medicine, which can be used to treat damp and fire, improve eyesight, relieve a cough, and reduce phlegm4. Antithrombotic enzyme obtained from Agkistrodon can be utilized to treat thrombus5,6. In particular, known as “liquid gold” in the international market, and costing dozens of times more than gold, snake venom is insufficiently available as animal medicine in the international market. Due to its proven antivirus7, anti-inflammatory8, and immune regulatory8 functions, snake venom not only has been in use as tranquilizer9, and for colorectal cancer treatment10, but also can be developed into health care products such as snake wines11,12.
Although diet elements or chemical compounds can be metabolized by microbes or co-metabolized by the host and microbes13,14, the relationship between venom secretion and gut microbiota requires further study. To investigate the flora difference between venom-secretion and non-venom secretion snakes, the gut samples of three kinds of snakes (e.g., Deinagkistrodon, Naja, and Ptyas mucosa) were collected for high-throughput 16 s rDNA sequencing analysis in this study. Our results revealed Bacteroidetes, Firmicutes, Fusobacteria, Proteobacteria, Actinobacteria, and Saccharibacteria to be the major bacteria present in these three snakes. Some core microbes, which can be detected in all samples, were identified in this study. Some bacteria may be related to the venom secretion by comparison analysis between two venom secretion snakes (Deinagkistrodon and Naja) and a non-venom secretion snake (Ptyas mucosa).
Materials and Methods
Sample collection
A total of five small intestinal samples and six large intestinal samples were collected from seven Ptyas mucosa; seven small intestinal samples and six large intestinal samples were collected from eight Deinagkistrodon; 13 small intestinal samples and 13 large intestinal samples were collected from 14 Naja. All of these samples were used for bacterial genome DNA extraction and high-throughput 16S rDNA sequencing. Detailed information of these snakes, like gender, weight and length et al., is listed in Supplementary Table 1. As in routine procedures, snakes were anesthetized using diethyl ether and then were killed by severing the spinal cord immediately posterior to the head. After the snakes were sacrificed, the snake skin was cleaned and wiped with sterile alcohol, and the gastrointestinal tract of these snakes was exposed via a mid-abdomen incision. The contents of the small intestine and large intestine were drained into different sterile vials. The collected samples were quickly frozen in liquid nitrogen, and then stored at −80 °C for further analysis. The study was approved by the Ethics Committee of the Hunan University of Science and Engineering, and all procedures were performed according to the relevant guidelines and regulations.
Bacterial genomic DNA extraction
Bacterial genomic DNA was extracted using a QIAamp DNA Stool Mini Kit following the manufacturer’s instructions (QIAGEN, Germany). The concentration of the extracted DNA was determined by using a NanoDrop ND-2000 (NanoDrop Technologies, USA), and its integrity and size were confirmed by performing 1% agar gel electrophoresis.
16S rRNA gene high-throughput sequencing and analysis
According to previous work, bacterial genome DNA samples extracted from snake gut samples were amplified using barcoded 16S rRNA gene primers 338F (5′-ACT CCT ACG GGA GGC AGC A-3′) with 806R (5′-GGA CTA CHV GGG TWT CTA AT-3′)15. After adding sequence adapts and constructing sequencing libraries, next-generation sequencing was performed on an Illumina MiSeq 300PE system that was operated by XY Biotechnology Co., Ltd., Shanghai, China. In total, 50 snake gut samples were sent to high-throughput 16S rRNA gene sequencing.
The original image data from the second-generation high-throughput sequencing instruments were converted into sequence data by base calling and saved as in a FASTQ format. Next-generation sequencing reads were identified by barcodes using a QIIME pipeline. Clean and high-quality sequences were then used for downstream analysis. A 97% similarity cutoff was employed in defining OTUs using Mothur16. One sequence was picked out from each OTU as representative. The representative sequences were classified using the RDP classifier method17,18 and SILVA database (https://www.arb-silva.de/documentation/release-123/). Good’s coverage, alpha diversities including Simpson and Shannon index, and richness (observed number of OTUs) values were calculated using Mothur. The data have been deposited in the sequence read archive (SRA) of NCBI as GenBank Accession Number SRP132421. Basic data on 16S rRNA gene high-throughput sequencing are listed in Supplementary Table 2.
Statistical Analysis
SPSS software (version 20.0; SPSS Inc., USA) and the ANOVA LSD(L) test were employed in this study. P < 0.05 was considered to be statistically significant. Plot cladograms and significantly different bacterial taxa were analyzed using a LDA Effect Size (LEfSe) algorithm19.
Results
Barcode 16s pyrosequencing and diversity analysis
In this study, 50 gut samples collected from three snake species, Deinagkistrodon, Naja, and Ptyas mucosa, were subjected to the 16s rDNA high-throughput sequencing analysis. In total, 1,301,357 valid reads were obtained from 13 Deinagkistrodon gut samples, 2,461,474 valid reads were obtained from 26 Naja gut samples, and 971,049 valid reads were obtained from 11 Ptyas mucosa gut samples (Supplementary Table 2). Of these sequences, 94.4%, 94.6%, and 99.9% reads were classified as bacterial phyla for Deinagkistrodon, Naja, and Ptyas mucosa samples, respectively. The reads that could be classified as bacteria were used for further analysis. As shown in Fig. 1A and Supplementary Figure 1, the sequencing coverage was close to 100% and the rarefaction curves exhibits a tendency to saturation, indicating that most bacterial species were captured by sequencing in all samples.
After classification, 29 bacterial phyla, 545 bacterial genera, and 1,725 OTUs were identified from these snake gut samples. In this study, 97% similarity was set as the OTU cut-off point. For each kind of snakes, 104–495 OTUs, 150–517 OTUs, and 91–335 OTUs were identified in Deinagkistrodon, Naja, and Ptyas mucosa samples, respectively (Supplementary Table 2). The distribution of OTU numbers in each sample is shown in Fig. 1B. The number of OTUs was slightly higher in venom secretion snakes Deinagkistrodon and Naja than non-venom secretion snake Ptyas mucosa. Although the difference of small intestinal samples was not statistically significant among these groups, OTU numbers were greatly higher in Naja than in Ptyas mucosa for the large intestinal samples (Table 1).
To analyze the alpha diversity of these samples, the indexes of Ace, Chao, Shannon, and Simpson were calculated. This process allowed us to fully characterize the bacterial community diversity in samples. Detailed information on the estimators for each sample is presented in Supplementary Table 2. As shown in Fig. 1C–F, while the indexes of Ace, Chao, and Shannon were slightly higher in Deinagkistrodon and Naja samples than Ptyas mucosa, the index of Simpson appears to be higher in Ptyas mucosa samples. Statistical analysis results showed that Ace, Chao, and Shannon indexes were significantly different between Naja and Ptyas mucosa small intestinal samples; Shannon and Simpson were significantly different between Deinagkistrodon and Naja small intestinal samples; Chao index was significantly different both between Deinagkistrodon to Ptyas mucosa, and Deinagkistrodon to Naja large intestinal samples (Table 1).
To analyze the beta diversity of the snake gut samples in this study, principal coordinate analysis (PCoA) was employed. As illustrated in Fig. 2, PCoA results indicated that most of the samples from Deinagkistrodon and Naja groups were clustered together with those of Ptyas mucosa groups being more dispersive.
Microbial composition and core gut microbes of snake gut samples
Ten bacterial phyla with a relative abundance higher than 1% were identified from sequence samples. Among them, six phyla Bacteroidetes, Firmicutes, Fusobacteria, Proteobacteria, Actinobacteria, and Saccharibacteria were detected in all samples. As shown in Fig. 3, the relative abundances of the top 10 bacterial phyla accounted for most of the sequenced reads. The phyla Bacteroidetes, Firmicutes, Fusobacteria, and Proteobacteria were the predominant bacteria, which accounted for 80.5–99.9% of the reads of all sequenced samples, with Actinobacteria and Saccharibacteria very low (Supplementary Figure 2). In 46 of 50 samples, more than 90% of sequenced reads were identified as bacterial phyla Bacteroidetes, Firmicutes, Fusobacteria, and Proteobacteria. The average percentages of these phyla in Deinagkistrodon, Naja, and Ptyas mucosa samples were 30.5%, 29.2%, and 23.2% for Bacteroidetes; 22.7%, 19.1%, and 27.7% for Firmicutes; 18.2%, 5.5%, and 4.4% for Fusobacteria; 24.4%, 40.1%, and 43.3% for Proteobacteria. Although the difference of the phyla Bacteroidetes and Firmicutes among Deinagkistrodon, Naja, and Ptyas mucosa samples was not significant, the relative abundance of phylum Fusobacteria was significantly higher in Deinagkistrodon than in Naja and Ptyas mucosa samples, and the relative abundance of the phylum Proteobacteria was significantly lower in Deinagkistrodon samples (Supplementary Figure 3A and B). Difference in some other minor bacterial phyla were also found from the pairwise comparison. Synergistetes was significantly higher in Deinagkistrodon than in Naja samples (Supplementary Figure 3A); Actinobacteria was significantly higher in Deinagkistrodon than in Ptyas mucosa samples (Supplementary Figure 3B); Actinobacteria, Synergistetes and Chlamydiae were significantly higher in Naja than in Ptyas mucosa samples (Supplementary Figure 3C).
At the genus level, 108 genera with relative abundance higher than 1% were identified from sequence samples. The first top 10 genera, Bacteroides, Cetobacterium, Aeromonas, Providencia, Morganella, Clostridium innocuum group, Achromobacter, Clostridium sensu stricto 1, Salmonella and Fusobacterium (Fig. 4), accounted for 57.0%, 52.2% and 64.3% of sequenced reads in Deinagkistrodon, Naja, and Ptyas mucosa samples, respectively.
The genera that were detected in each sample were considered as core genera. As summarized in Table 2, 27, 24, and 16 genera were identified as core microbes for Deinagkistrodon, Naja, and Ptyas mucosa, respectively. The relative abundance of the core genera accounted for more than 50% of the sequence reads in most of the sequenced samples (Supplementary Figure 4A–C). 12 genera were detected in all sequenced snake samples, independent of breed (Table 2). The distribution of the 12 genera in all samples is presented in Fig. 5. The average abundance of the 12 core genera was upto 43.9%, and most of the genera like Acinetobacter20, Aeromonas21, Citrobacter22, Morganella23, Proteus24, Providencia25, Pseudomonas26, and Salmonella27 were human and animal pathogens or opportunistic pathogens (Table 2). Interestingly, genera Achromobacter, Cetobacterium, Clostridium innocuum group, Fusobacterium, Lachnoclostridium, Parabacteroides, Romboutsia can only be detected in venom secretion snakes Deinagkistrodon and Naja (Table 2).
Although, some genera can be identified as core microbes, the difference among different breeds was obvious. Between Deinagkistrodon and Naja, 42 genera were significantly higher in Deinagkistrodon samples, and another 56 genera were significantly higher in Naja samples (Supplementary Figure 5A). Between Deinagkistrodon and Ptyas mucosa, 49 genera were significantly higher in Deinagkistrodon samples, and another 10 genera were significantly higher in Ptyas mucosa samples (Supplementary Figure 5B). Between Naja and Ptyas mucosa, 68 genera were significantly higher in Naja samples, and another 11 genera were significantly higher in Ptyas mucosa samples (Supplementary Figure 5C).
Intestinal flora difference among studied snakes
Based on LEfSe analyses, the significantly different bacterial taxa (P < 0.05) among the three communities were identified and listed in the right side of the cladogram (Fig. 6). At the Class level, Epsilonproteobacteria was higher in Naja than other two breeds, and Gammaproteobacteria was higher in Ptyas mucosa than other two breeds. At the Order level, Vibrionales, Campylobacterales, and Pseudonocardiales were significantly higher in Naja samples. At the Family level, Rikenellaceae was significantly higher in Deinagkistrodon samples, with Helicobacteraceae, Family XI under Order Clostridiales, Veillonellaceae, and Pseudonocardiaceae significantly higher in Naja samples. At the genus level, uncultured bacteria under Family Flavobacteriaceae, Oxalobacter, other unidentified bacteria under Family Lachnospiraceae, Porphyromonas and Desulfovibrio were significantly higher in Deinagkistrodon samples; Helicobacter, Leucobacter, Arachidicoccus, Peptostreptococcus, Erysipelothrix, Saccharopolyspora, and Mobiluncus were significantly higher in Naja samples; Providencia, Hafnia_Obesumbacterium, and Ureaplasma were significantly higher in Ptyas mucosa samples (Fig. 6).
Discussion
Although a few studies have used high-thoughput sequencing method to investigate the bacterial communities of Agkistrodon piscivorus28, Crotalus horridus29 and Burmese pythons30, no bacterial taxa data of gut microbiota from Deinagkistrodon, Naja and Ptyas mucosa have been reported so far. In this study, 29 bacterial phyla, 545 bacterial genera, and 1,725 OTUs were identified in the gut samples of these snakes. The average OTU number of these samples was about 300 to 400, which was similar to the OTU numbers detected in Burmese pythons samples30, but a little bit higher than Agkistrodon piscivorus samples28 for which about 100 OTUs were detected. The OTU number difference among our study with previously published reports maybe due to their different life habits, or dependence on the sequence depth.
In this study, the phyla Bacteroidetes, Firmicutes, Fusobacteria, and Proteobacteria were the predominant bacteria, which were consistent with results obtained from Agkistrodon piscivorus28 and Crotalus horridus29. Bacteroidetes, Firmicutes, and Proteobacteria were reported as the predominant bacteria detected in Agkistrodon piscivorus samples28, and Proteobacteria and Firmicutes were reported as the predominant bacteria detected in Crotalus horridus samples29. However, Proteobacteria was the major bacterial phylum, accounting for more than 50% in Agkistrodon piscivorus small intestinal samples28, and Bacteroidetes was the predominant bacterium in large intestinal samples. In this study, bacterial composition patterns of small and large intestines at the phylum level were very similar (Fig. 3), the similarity coefficients of these patterns between small and large intestines for Deinagkistrodon, Naja Ptyas and Ptyas mucosa were 0.98, 0.96 and 0.78, respectively. The bacterial communities were not clustered together according to gender, weight, and length (Fig. 2). We know from Supplementary Table 1, that snakes D1, D2, D3 and D8 were female Deinagkistrodon, and snakes D4, D5, D6 and D7 were male Deinagkistrodon. PCoA analyses did not separate these samples according to different sex (Fig. 2). Similar results were also found for different weights and lengths (Fig. 2). In addition, the wild snake didn’t have an obvious different bacterial community when compared to bred snakes (Fig. 2). However, Deinagkistrodon and Naja groups were clustered together with those of the Ptyas mucosa group being more disperse. The character and strength of venom toxicity may be one of the major factors that affected their intestinal bacterial community. Both snakes of Deinagkistrodon and Naja can secret venom (https://en.wikipedia.org/wiki/), thus the venom may have the ability to affect the growth of intestinal bacteria. However, the snake of Ptyas mucosa didn’t have the ability to secret venom, so there was no venom to restrict their intestinal bacterial patterns. Of course, our hypothesis requires further investigation.
Despite the fact that core microbiomes have been identified for human and mouse gut samples31,32,33, snake core gut microbes are still missing due to the small number of samples from previous snake gut microbiota studies. In this study, we sequenced 50 gut samples, and found that the core microbes existed in these snake samples, and that different breeds had different core microbiomes. As shown in Table 2, 27, 24, and 16 genera were identified as core microbes for Deinagkistrodon, Naja, and Ptyas mucosa, respectively. Moreover, some bacterial genera were found existing in all detected snakes (Fig. 5). In consideration that most genera listed in Table 2, like Acinetobacter20, Aeromonas21, Citrobacter22, Morganella23, Proteus24, Providencia25, Pseudomonas26, and Salmonella27 were human and animal pathogens or opportunistic pathogens, they should be defined as one of the focuses of disease prevention and control. Since Ptyas mucosa and Naja species (mostly Naja naja) frequent human dwellings, their feces may contaminate human and animals, and become a bacterial infection source. A total of 17 bacterial genera identified from intestines of study snakes (Table 2) were found to be pathogens or opportunistic pathogens to animals and humans (https://en.wikipedia.org/wiki/Main_Page). Therefore, this study provides an opportunity for further investigations on how snakes contribute to zoonosis in their distribution ranges.
However, some significant differences in some bacteria were identified among different breeds. More interestingly, we found that some bacteria could be related to the venom secretion. As illustrated in Fig. 6, Hafnia_Obesumbacterium, Providencia, and Ureaplasma were significantly higher in non-venom secretion snake Ptyas mucosa than the other two venom secretion snakes Deinagkistrodon and Naja. The other 7 genera, Achromobacter, Cetobacterium, Clostridium innocuum group, Fusobacterium, Lachnoclostridium, Parabacteroides, and Romboutsia were only detected in the venom secretion snakes Deinagkistrodon and Naja. Certainly, the function of these bacteria in venom secretion needs to be further studied. In addition, the study of the varied buccal cavity bacterial community of snakes is needed. Some pathogenic bacteria has been cultured from the venom and oropharynx of snakes, and they have been implicated in secondary infections of snake bite wounds34,35.
In summary, the bacterial community of 50 snake gut samples was examined using 16S rDNA sequencing analysis, and 29 bacterial phyla, 545 bacterial genera, and 1,725 OTUs were identified in these samples. Among them, 27, 24, and 16 genera were identified as core microbes for Deinagkistrodon, Naja, and Ptyas mucosa, respectively. In addition, the relative abundance of 3 genera (Providencia, Hafnia_Obesumbacterium, and Ureaplasma) were found to be significantly higher in non-venom secretion snakes, and 7 genera (Achromobacter, Cetobacterium, Clostridium innocuum group, Fusobacterium, Lachnoclostridium, Parabacteroides, Romboutsia) were only detected in venom secretion snakes. There is a need for additional sophisticated study to understand whether these bacteria have targets to improve venom secretion in those snakes.
References
Makki, K., Deehan, E. C., Walter, J. & Backhed, F. The Impact of Dietary Fiber on Gut Microbiota in Host Health and Disease. Cell Host Microbe. 23, 705–715 (2018).
Jandhyala, S. M. et al. Role of the normal gut microbiota. World J Gastroenterol. 21, 8787–8803 (2015).
Xu, W. T., Nie, Y. Z., Yang, Z. & Lu, N. H. The crosstalk between gut microbiota and obesity and related metabolic disorders. Future Microbiol. 11, 825–836 (2016).
Sun, H. L. & Li, C. Q. The Pharmacological Action and Application of Snake’s Gallbladde. r. Northwest Pharmaceutical Journal. 19, 285–287 (2004).
Hao, W. X. Effect of an antithrombotic enzyme of Agkistrodon halys venom in the treatment of cerebral thrombosis–report of 322 cases. Zhonghua Shen Jing Jing Shen Ke Za Zhi. 17, 223–225 (1984).
Ding, B. et al. Antiplatelet Aggregation and Antithrombosis Efficiency of Peptides in the Snake Venom of Deinagkistrodon acutus: Isolation, Identification, and Evaluation. Evid Based Complement Alternat Med. 2015, 412841 (2015).
Yu, M. Y. Ozagrel sodium combined with snake anti-thrombosis enzyme therapy for cerebral hemorrhage clinical analysis. Chinese Medicine Guide. 15, 167–168 (2017).
Wang, S. Z. & Qin, Z. H. Anti-Inflammatory and Immune Regulatory Actions of Naja naja atra Venom. Toxins (Basel). 10 (2018).
Osipov, A. & Utkin, Y. Effects of snake venom polypeptides on central nervous system. Cent Nerv Syst Agents Med Chem. 12, 315–328 (2012).
Uzair, B. et al. Snake venom as an effective tool against colorectal cancer. Protein Pept Lett. 25, 626–632 (2018).
Tang, S. D. Analysis about the process and function of yongzhou Elapidae healthful liquor. Food Research and Development. 28, 149–152 (2007).
Bin, D. M., Zhou, J. T., Zhong, F. S. & Liu, J. Study on the Formula Design and Technology of Yongzhou Different snake wine. Food and Machinery. 22, 130–132 (2006).
De Angelis, M. et al. The food-gut human axis: the effects of diet on gut microbiota and metabolome. Curr Med Chem (2017).
Koppel, N., Maini Rekdal, V. & Balskus, E. P. Chemical transformation of xenobiotics by the human gut microbiota. Science. 356 (2017).
Dennis, K. L. et al. Adenomatous polyps are driven by microbe-instigated focal inflammation and are controlled by IL-10-producing T cells. Cancer Res. 73, 5905–5913 (2013).
Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 75, 7537–7541 (2009).
Cole, J. R. et al. The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, D141–145 (2009).
Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 73, 5261–5267 (2007).
Segata, N. et al. Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60 (2011).
Amorim, A. M. & Nascimento, J. D. Acinetobacter: an underrated foodborne pathogen? J Infect Dev Ctries. 11, 111–114 (2017).
Batra, P., Mathur, P. & Misra, M. C. Aeromonas spp.: An Emerging Nosocomial Pathogen. J Lab Physicians. 8, 1–4 (2016).
Collins, J. W. et al. Citrobacter rodentium: infection, inflammation and the microbiota. Nat Rev Microbiol. 12, 612–623 (2014).
Liu, H., Zhu, J., Hu, Q. & Rao, X. Morganella morganii, a non-negligent opportunistic pathogen. Int J Infect Dis. 50, 10–17 (2016).
Schaffer, J. N. & Pearson, M. M. Proteus mirabilis and Urinary Tract Infections. Microbiol Spectr. 3 (2015).
Sagar, S., Narasimhaswamy, N. & D’Souza, J. Providencia Rettgeri: An Emerging Nosocomial Uropathogen in an Indwelling Urinary Catheterised Patient. J Clin Diagn Res. 11, DD01–DD02 (2017).
Juan, C., Pena, C. & Oliver, A. Host and Pathogen Biomarkers for Severe Pseudomonas aeruginosa Infections. J Infect Dis. 215, S44–S51 (2017).
Crowley, S. M., Knodler, L. A. & Vallance, B. A. Salmonella and the Inflammasome: Battle for Intracellular Dominance. Curr Top Microbiol Immunol. 397, 43–67 (2016).
Leser, T. D. et al. Culture-independent analysis of gut bacteria: the pig gastrointestinal tract microbiota revisited. Appl Environ Microbiol. 68, 673–690 (2002).
McLaughlin, R. W., Cochran, P. A. & Dowd, S. E. Metagenomic analysis of the gut microbiota of the Timber Rattlesnake, Crotalus horridus. Mol Biol Rep. 42, 1187–1195 (2015).
Costello, E. K., Gordon, J. I., Secor, S. M. & Knight, R. Postprandial remodeling of the gut microbiota in Burmese pythons. ISME J. 4, 1375–1385 (2010).
Barzegari, A., Saeedi, N. & Saei, A. A. Shrinkage of the human core microbiome and a proposal for launching microbiome biobanks. Future Microbiol. 9, 639–656 (2014).
Li, J. et al. An integrated catalog of reference genes in the human gut microbiome. Nat Biotechnol. 32, 834–841 (2014).
Xiao, L. et al. A catalog of the mouse gut metagenome. Nat Biotechnol. 33, 1103–1108 (2015).
Suankratay, C., Wilde, H., Nunthapisud, P. & Khantipong, M. Tetanus after white-lipped green pit viper (Trimeresurus albolabris) bite. Wilderness Environ Med. 13, 256–261 (2002).
Theakston, R. D. et al. Bacteriological studies of the venom and mouth cavities of wild Malayan pit vipers (Calloselasma rhodostoma) in southern Thailand. Trans R Soc Trop Med Hyg. 84, 875–879 (1990).
Acknowledgements
This study was funded by the National Nature Science Foundation (NSFC, No. 21606079 and No. 31741109), the key research and development plan of Hunan Province (No. 2018WK2093, No. 2015CK3002 and 2019JJ3011), and the construct program of applied characteristic discipline in Hunan University of Science and Engineering. We thanks to Dr. Huaqiang Zeng (Institute of Bioengineering and Nanotechnology, Singapore) for linguistic assistance.
Author information
Authors and Affiliations
Contributions
Y.Y., Z.Q. designed the experiments; S.W., D.G., and H.C. carried out the research; J. Z., L.B., and X.L. contributed reagents/materials/analysis tools; Y.Y., Z.Q., and S.W. analyzed the data; Y.Y. and S.W. wrote the manuscript. All authors discussed the results and commented on the manuscript.
Corresponding author
Ethics declarations
Competing Interests
Authors Dezhi Guo and Le Bai are members of Yongzhou City Strange Snake Science and Technology Industrial Co., Ltd. Those authors provided the snake samples. The other authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Qin, Z., Wang, S., Guo, D. et al. Comparative analysis of intestinal bacteria among venom secretion and non-secrection snakes. Sci Rep 9, 6335 (2019). https://doi.org/10.1038/s41598-019-42787-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-42787-6
- Springer Nature Limited