Diets Alter the Gut Microbiome of Crocodile Lizards

The crocodile lizard is a critically endangered reptile, and serious diseases have been found in this species in recent years, especially in captive lizards. Whether these diseases are caused by changes in the gut microbiota and the effect of captivity on disease remains to be determined. Here, we examined the relationship between the gut microbiota and diet and disease by comparing the fecal microbiota of wild lizards with those of sick and healthy lizards in captivity. The gut microbiota in wild crocodile lizards was consistently dominated by Proteobacteria (∼56.4%) and Bacteroidetes (∼19.1%). However, the abundance of Firmicutes (∼2.6%) in the intestine of the wild crocodile lizards was distinctly lower than that in other vertebrates. In addition, the wild samples from Guangdong Luokeng Shinisaurus crocodilurus National Nature Reserve also had a high abundance of Deinococcus–Thermus while the wild samples from Guangxi Daguishan Crocodile Lizard National Nature Reserve had a high abundance of Tenericutes. The gut microbial community in loach-fed crocodile lizards was significantly different from the gut microbial community in the earthworm-fed and wild lizards. In addition, significant differences in specific bacteria were detected among groups. Notably, in the gut microbiota, the captive lizards fed earthworms resulted in enrichment of Fusobacterium, and the captive lizards fed loaches had higher abundances of Elizabethkingia, Halomonas, Morganella, and Salmonella, all of which are pathogens or opportunistic pathogens in human or other animals. However, there is no sufficient evidence that the gut microbiota contributes to either disease A or disease B. These results provide a reference for the conservation of endangered crocodile lizards and the first insight into the relationship between disease and the gut microbiota in lizards.

Introduction

The crocodile lizard (Shinisaurus crocodilurus Ahl, 1930) is the only species in the monotypic genus Shinisaurus and the monotypic family Shinisauridae. It is a relict reptile that now survives only in separated Pleistocene refugia. This species is distributed in southern China (Guangdong and Guangxi Provinces) and northern Vietnam (Qung Ninh and Bac Giang Provinces) with severely fragmented populations (van Schingen et al., 2014, 2016). This species is essential not only for taxonomical systematics but also for understanding the origin, adaptation, and evolution of reptiles. However, it faces extinction due to the pressure of being hunted; environmental changes; and habitat destruction (Huang et al., 2008; Nguyen and Ziegler, 2015). It has been listed as an endangered species on the IUCN Red List of Threatened Species (Nguyen et al., 2014), a class I protected species in China, and an appendix I species by the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES I). According to a recent survey, the total number of crocodile lizards in the wild has decreased from 6000 in 1978 to approximately 1200 in China (Huang et al., 2008; unpublished survey conducted by Wu et al., 2012). What’s worse, the population continues to decline sharply. Similarly, the wild population in Vietnam has decreased to fewer than 150 individuals in recent years (van Schingen et al., 2016). Some nature reserves, such as Guangdong Luokeng S. crocodilurus National Nature Reserve and Guangxi Daguishan Crocodile Lizard National Nature Reserve, are conducting captive breeding and release programs with the hope of restoring the wild populations. However, the captive individuals can become infected with serious diseases that cause many deaths each year (Figure 1). Crocodile lizards have two major types of diseases that can be distinguished by their symptoms. Disease A is characterized by one or more nodules in the underjaw or limbs covered by lesions (Figure 1A), while disease B is characterized by varying degrees of rot in the four limbs (Figure 1B). In the wild, individuals with disease B have been observed, while individuals with disease A have not been found.

FIGURE 1

FIGURE 1. Symptoms of disease (A) and (B). Arrows indicate the lesions.

Recent studies have revealed how variations and changes in the composition of gut microbial communities influence normal physiology and contribute to diseases (Clemente et al., 2012; Martin et al., 2014; Boursier et al., 2016). In addition, gut microbes affect host immunity, behavior, reproductive isolation, and metabolism (Cryan and Dinan, 2012; Brucker and Bordenstein, 2013; Ramakrishna, 2013; Thaiss et al., 2016; Sylvia et al., 2017). Conversely, many factors such as diet or host genetics can shape the microbial community (David et al., 2014; Goodrich et al., 2014).

Study on the gut microbiota has been conducted in a host of vertebrates, including mammals (Thaiss et al., 2016), birds (Hird et al., 2015; Waite and Taylor, 2015), fishes (Gajardo et al., 2016), amphibians (Bletz et al., 2016), and reptiles (Costello et al., 2010; Yuan et al., 2015; Kohl et al., 2016, 2017; Ren et al., 2016). However, the majority of these studies have been conducted in mammalian hosts. Surveys of the gut microbiota in reptiles, an ancient group with more than 10,000 extant species (Uetz et al., 2016), remain rare except in the case of some economically important species, such as snakes (Costello et al., 2010; Colston et al., 2015) and turtles (Huang and Zhang, 2013; Yuan et al., 2015). A few studies have been conducted on the gut microbiota of lizards, which represent about 60% of reptiles (Martin et al., 2010; Hong et al., 2011; Kohl et al., 2016, 2017; Ren et al., 2016). Nevertheless, the interaction between the gut microbiome and disease in reptiles remains unclear.

According to previous studies, the reptilian immune system differs from those of other vertebrates in several aspects (Zimmerman et al., 2010). Here, two questions are raised. Is the gut microbiota associated with disease susceptibility in crocodile lizards? How does cultivation shape the gut microbiome of crocodile lizards? To explore these questions and to facilitate the protection of this endangered species, we analyzed the amplicon-based microbiome of cloacal swab samples from crocodile lizards using 16S rRNA gene sequencing. We compared the gut microbiota of sick, healthy, captive, and wild lizards to identify the interactions among specific diseases, diets, and the gut microbiota in crocodile lizards.

Materials and Methods

Sample Collection

All samples were collected from Guangdong Luokeng S. crocodilurus National Nature Reserve (referred to as “Luokeng Nature Reserve” in the following sections) and Guangxi Daguishan Crocodile Lizard National Nature Reserve (referred to as “Daguishan Nature Reserve” in the following sections). Thirty crocodile lizards were separated into six groups, namely, the wild group from Luokeng Nature Reserve (WLK, n = 7), the healthy earthworm-fed group (NLK, n = 4), the sick earthworm-fed group with disease A (SLK, n = 5), the wild group from Daguishan Nature Reserve (WDG, n = 8), the healthy loach-fed group (NDG, n = 3), and the sick loach-fed group with disease B (SDG, n = 3). In addition, because of the highly similarity, the sick and healthy groups that fed the same diet were merged and recalculated. The sick and healthy individuals that fed earthworm were merged as earthworm-fed group (CLK), and the sick and healthy individuals that fed loach were merged as loach-fed group (CDG). Detailed sample information is shown in Table 1. Cloacal swabs were used for nondestructive sampling of the gut microbiota (Colston et al., 2015). The cloacal swabs were collected and stored in absolute ethyl alcohol or liquid nitrogen and then transported to the lab for DNA extraction within 24 h.

TABLE 1

TABLE 1. Sample information used in this study.

All experimental animal procedures were approved by the Committee on the Ethics of Animal Experiments of the Guangdong Institute of Applied Biological Resources following basic principles.

DNA Extraction and Sequencing

Total DNA was extracted from the cloacal swabs using a PowerFecal^® DNA Isolation Kit (MOBIO Laboratories, Inc., United States). The V4 hypervariable region of the 16S rRNA gene was amplified with the primers 515F (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′), followed by library preparation using an NEB Next^® Ultra^TM DNA Library Prep Kit for Illumina (NEB, United States). Sequencing on an Illumina HiSeq platform (250 bp paired-end reads) was performed by Novogene Corporation (Beijing, China).

Data Analysis

Raw tags were filtered using the QIIME V1.7.0 package (Caporaso et al., 2010) in order to remove the low-quality sequences and chimeras. Then, sequences with ≥97% similarity were assigned to the same operational taxonomic units (OTUs) using UCLUST in QIIME V1.7.0 package (Caporaso et al., 2010). A representative sequence for each OTU was annotated with threshold 0.8 using RDP Classifier 2.2 by searching the SILVA database (Wang et al., 2007; Quast et al., 2013).

For comparisons between samples, the OTU abundances were normalized by the number obtained from the sample with the lowest counts.

For each sample, alpha diversity was estimated by calculating the Shannon and abundance-based coverage estimator (ACE) indices. These indices were calculated by QIIME 1.7.0 (Caporaso et al., 2010) and displayed using R software. Alpha diversity indices were compared among samples using the Tukey method (P = 0.05) with R software.

Beta diversity was measured by principal coordinate analysis (PCoA) on unweighted and weighted UniFrac distances and were displayed using R software. The unweighted and weighted UniFrac distances were calculated by QIIME 1.7.0 (Caporaso et al., 2010). In addition, unweighted pair-group method with arithmetic means (UPGMA) clustering was also performed using QIIME 1.7.0 (Caporaso et al., 2010). The unweighted UniFrac distance accounts for membership in a community whereas the weighted UniFrac distance considers both membership and the relative abundance. Permutational multivariate analysis of variance (PERMANOVA) statistical analyses were conducted based on unweighted and weighted UniFrac distances with 999 permutations using function adonis in R’s vegan package.

To identify microbes accounting for the effects of disease and diet, the linear discriminatory analysis (LDA) effect size (LEfSe) method was used to compare the differential abundances of bacteria among groups at family and genus levels. LEfSe analysis emphasizes statistical significance, biological consistency, and effect relevance. It first robustly identifies taxa that are statistically different among groups. Then it investigates biological consistent using a set of pairwise tests among subgroups. At last, it uses LDA to estimate the effect size of each selected taxon. LEfSe analysis was performed using LEfSe software (Segata et al., 2011). The threshold of P-value in the Kruskal–Wallis test among groups was 0.05. Only those taxa with a log LDA score >4 (more than four orders of magnitude) were considered in this study.

All raw sequences obtained in this study have been deposited in the Sequence Read Archive (SRA) under accession number SRP107074.

Results

We analyzed the bacterial composition of 30 crocodile lizard cloacal swab samples (Table 1). Each sample contained at least 30,000 effective sequences (Supplementary Figure 1). The rarefaction curves showed that these sequence depths were sufficient for capturing the major microbiota in each sample (Supplementary Figure 2). More than 99% of the OTUs could be well annotated at the family level in each sample (Supplementary Figure 3).

General Pattern of the Gut Microbiota of the Wild Crocodile Lizard

In the total dataset, most of the bacteria were identified as Proteobacteria (47.9%) and Bacteroidetes (32.1%). At the phylum level, majority of species in the wild crocodile lizard gut microbiota were classified as Proteobacteria (56.4%), Bacteroidetes (19.1%), and Firmicutes (2.6%). In addition, the wild samples from Luokeng Nature Reserve also had a high abundance of Deinococcus–Thermus (13.6%) while the samples from Daguishan Nature Reserve had a high abundance of Tenericutes (5.9%). At the family level, the most abundant taxa in the wild crocodile lizard gut microbiota from Luokeng Nature Reserve were Pasteurellaceae, Deinococcaceae, Comamonadaceae, and Flavobacteriaceae. However, in Daguishan Nature Reserve, the wild crocodile lizard gut microbiota was dominated by Helicobacteraceae, Mycoplasmataceae, Pseudomonadaceae, and Chitinophagaceae. In Luokeng Nature Reserve, the most frequently occurring genera in the wild crocodile lizard intestine were Niabella, Deinococcus, Alysiella, and Chryseobacterium. In Daguishan Nature Reserve, the most frequently occurring genera in the wild crocodile lizard intestine were Helicobacter, Mycoplasma, Pseudomonas, and Niabella. However, the wild samples from Daguishan Nature Reserve were obviously separated into two patterns. Four samples had an extremely high abundance of Mycoplasma and Helicobacter (Figure 2).

FIGURE 2

FIGURE 2. Composition of the gut microbiota of crocodile lizards at the phylum (A), family (B), and genus (C) levels.

The ACE and Shannon indices of the wild crocodile lizard gut microbiota from Luokeng Nature Reserve were significantly higher than that from Daguishan Nature Reserve (Figure 3 and Supplementary Tables 1, 2). Therefore, both community richness and community diversity of the gut microbiota in wild crocodile lizards were different between distinct locations.

FIGURE 3

FIGURE 3. The alpha diversity of the gut microbial composition. (A) ACE index and (B) Shannon index. The bottom and top of the box are the first and third quartiles, the band inside the box is the median, and the ends of the whisker present the minimum and maximum.

Comparison of Gut Microbial Community Diversity between Groups

A comparison of alpha diversity indices between groups is presented in Figure 3. The significance between groups was detected by turkey method. For the community richness estimator (the ACE index), there was significant difference between earthworm-fed group and loach-fed group (NLK versus NDG). In addition, the loach-fed group was notably dissimilar with the wild group from the same place (NDG versus WDG), while the earthworm-fed group had no significant difference with the wild group from the same place (NLK versus WLK). However, there were no detectable significant differences between healthy and sick groups with the same diet. For the community diversity estimator (the Shannon index), no significant differences were found between these groups (Figure 3 and Supplementary Tables 1, 2).

The PCoA results also showed that samples from earthworm-fed group and loach-fed group were distantly separated for both community membership (Figure 4A) and community diversity (Figure 4B). In addition, the wild samples were separated from the loach-fed samples in Daguishan Nature Reserve. Nevertheless, for community diversity, four wild samples (WDG02, WDG03, WDG05, and WDG10) from Daguishan Nature Reserve were distantly separated from the other samples. These differences were also revealed by UPMA clustering (Figure 5). The separation of the four samples (WDG02, WDG03, WDG05, and WDG10) was resulted by the abnormally high abundances of Mycoplasma and Helicobacter, which decrease the community diversity. The earthworm-fed individuals were similar with the wild samples (Figures 4, 5). The sick and healthy groups overlapped for both community membership and community diversity (Figures 4, 5). In addition, the observed clusters were also supported by PERMANOVA analyses based on weighted and unweighted UniFrac metrics (P = 0.001, Supplementary Table 3).

FIGURE 4

FIGURE 4. The beta diversity of the gut microbial composition. Principal coordinate analysis (PCoA) was conducted based on unweighted (A) and weighted (B) UniFrac distance matrices. The variation explained by the plotted principal coordinates is indicated in the axis labels.

FIGURE 5

FIGURE 5. The beta diversity of the gut microbial composition. The UPGMA tree was generated based on weighted UniFrac distances.

Together with the results of alpha and beta diversity analyses, the gut microbiota in the sick and healthy groups was highly similar. Therefore, to well interpret the effect of diet, the sick and healthy samples that fed same diet were merged. The sick and healthy individuals that fed earthworm were merged as earthworm-fed group (CLK), and the sick and healthy individuals that fed loach were merged as loach-fed group (CDG). Then we recalculated the difference of alpha diversity indices between earthworm-fed group and loach-fed group. The results confirmed that the differences between groups of different diets were significant in terms of community richness (Supplementary Table 1).

In conclusion, there were significant differences in community richness and membership between groups of different diets. But there were no significant differences in the community diversity between any groups. However, no significant differences were found between healthy and sick groups in both community richness and community diversity.

Differential Microbes among Groups

The LEfSe analysis was used to screen the differential microbes among groups. Fourteen genera and 13 families were significantly enriched in distinct groups (Figure 6). The relative abundance of each selected genus is presented in Figure 7.

FIGURE 6

FIGURE 6. Differences in bacterial taxa among groups determined by linear discriminative analysis effect size (LEfSe). The highlighted taxa were significantly enriched in the group that corresponds to each color. LDA scores can be interpreted as the degree of difference in relative abundance. Abbreviation: g_, genus and f_, family.

FIGURE 7

FIGURE 7. Relative abundances of differential microbes among groups. The taxa were selected by LEfSe analysis at the genus level. The straight line is the mean. The dot line is the median. (A) Dominant bacteria in the gut of wild crocodile lizards from Luokeng Nature Reserve. (B) Dominant bacteria in the gut of wild crocodile lizards from Daguishan Nature Reserve. (C) Dominant bacteria in the gut of captive crocodile lizards fed earthworms. (D) Dominant bacteria in the gut of captive crocodile lizards fed loaches.

When compared the effects of diet on the crocodile lizard gut microbiota, the earthworm-fed group (CLK) showed increased abundances of genera Petrimonas, Bacteroides, Fusobacterium, Microbacter, and Proteiniclasticum, and families Porphyromonadaceae, Bacteriodaceae, Fusobacteriaceae, and Clostridiaceae. The loach-fed lizards (CDG) had a significant higher abundance of genera Elizabethkingia, Halomonas, and Morganella, and families Flavobacteriaceae, Enterobacteriaceae, and Halomonadaceae (Figure 6). When checking the relative abundance of each genus, Salmonella appeared in all samples, and significantly enriched in the gut of loach-fed lizards (Figure 8).

FIGURE 8

FIGURE 8. Relative abundances of Salmonella.

Compared with captive crocodile lizards, the gut microbiota of wild lizards from Luokeng Nature Reserve showed increased significantly in the abundances of genera Deinococcus and Lysobacter, and families Pasteurellaceae, Deinococcaceae, and Comamonadaceae. However, the wild lizards from Daguishan Nature Reserve had higher abundances of genera Helicobacter, Mycoplasma, Pseudomonas, and Niabella, and families Helicobacteraceae, Mycoplasmataceae, and Pseudomonadaceae (Figure 6).

Discussion

Determining the role of the intestinal bacterial community in digestion and pathogenesis depends critically upon defining the “wild” and “normal” states. With high-throughput sequencing, it is now possible to comprehensively identify the bacteria in a given community, including fastidious and unculturable taxa. This study aimed to describe, for the first time, the wild state of the crocodile lizard gut microbiota as well as potential alterations in the composition of the gut microbiota of crocodile lizards with different diseases and diets compared with that in wild and healthy control subjects.

The composition of the gut microbiota in crocodile lizards is unique. Previous reports have indicated that the gut microbiota of lizards is dominated by the phyla Firmicutes (33.2–73%), Bacteroidetes (6.2–45.7%), and Proteobacteria (5.7–62.3%) (Martin et al., 2010; Nelson et al., 2010; Hong et al., 2011; Ren et al., 2016; Kohl et al., 2017). In other reptiles, the gut microbiota also appeared to be consistently dominated by Firmicutes, followed by Bacteroidetes, while Proteobacteria ranged from dominant to minor components (Costello et al., 2010; Colston et al., 2015; McLaughlin et al., 2015; Yuan et al., 2015). An exception is the alligator gut microbiota, which was dominated by Fusobacteria (Keenan et al., 2013). In other vertebrates, Firmicutes and Bacteroidetes also represent overwhelming majority of the gut microbiome (summarized by Keenan et al., 2013). Conversely, in the crocodile lizard gut microbiota, the proportion of Firmicutes was low, in all samples (0.13–14.56%), while Proteobacteria and Bacteroidetes were dominant (Figure 2). Moreover, in the wild samples from Luokeng Nature Reserve, the high prevalence of phylum Deinococcus–Thermus in the gut microbiome has not been found in other vertebrates. Deinococcus–Thermus spp. are usually found in extreme environments (Theodorakopoulos et al., 2013). They have also been detected in the feces of some animals with only a few clones (Lagier et al., 2012; McLaughlin et al., 2012). The role of Deinococcus–Thermus in the intestinal tract is not clear. However, Deinococcus–Thermus bacteria have been reported to show remarkable resistance to a range of stresses such as ionizing radiation, UV radiation, oxidizing agents, and desiccation (Theodorakopoulos et al., 2013). In addition, the relative abundance of Deinococcus–Thermus was significantly higher in the wild group than in the captive groups from the same place (Figures 2, 7). Therefore, the high proportion of Deinococcus–Thermus may help crocodile lizards adapt to the wild environment in Luokeng Nature Reserve. The wild samples from the Daguishan Nature Reserve had high abundances of Mycoplasma and Helicobacter. Particularly, the total abundances of Mycoplasma and Helicobacter reached up to 71.2–91.4% in samples WDG02, WDG03, WDG05, and WDG10. Because Mycoplasma and Helicobacter were detected in all samples, they are likely commensal inhabitants of crocodile lizards. Mycoplasma is a genus of bacteria that lack a cell wall and is primarily obligate commensals or parasites. Some species of Mycoplasma are significant pathogens of birds, mammals, fish, and reptiles, although many of them are harmless commensals to the hosts (Bano et al., 2007; Ossiboff et al., 2015). The overwhelming dominance of Mycoplasma have been reported in the fish gut or oyster stomach (Bano et al., 2007; King et al., 2012). However, the dominance of Mycoplasma in reptile guts was first reported. In the terms of Helicobacter, the most widely known member is H. pylori, which is strongly associated with peptic ulcers, chronic gastritis, duodenitis, and stomach cancer in humans (Khalifa et al., 2010). However, it is unclear whether the dominant Mycoplasma and Helicobacter would be harmful for the crocodile lizards.

Many studies showed that the gut microbiota had a relationship with a variety of diseases. In crocodile lizards, there was no significant correlation between the two diseases (A and B) and the gut microbiota. This may be because both disease A and B are located on the skin surface and do not have direct contact with the gut system.

Previous studies have shown that captivity can change the diversity of the gut microbiota (Nelson et al., 2013; Kohl and Dearing, 2014). In a study of Anolis sagrei insectivorous lizards also demonstrated that captivity led to a shift in microbial diversity (Ren et al., 2016). However, an investigation of omnivorous and herbivorous lizards showed that captivity had no significant effect on gut microbial diversity in terms of alpha diversity (Kohl et al., 2017). In this study, a comparison of alpha diversity indices, PCoA analysis, and UPGMA cluster suggested that the effect of captivity in shaping of gut microbiota was depend on diet. The diet of loach significantly changed the community richness of the gut microbiota of crocodile lizards but the diet of earthworm did not. However, captivity had no effect on the community diversity of gut microbiota of crocodile lizards according to the Shannon index.

The effect of captivity in shaping of gut microbiota was also reflected in the enrichment of specific microbes. The earthworm-fed group exhibited a gut microbiota enriched in Petrimonas, Bacteroides, Fusobacterium, Microbacter, and Proteiniclasticum compared with the wild and loach-fed groups (Figures 6, 7). Moreover, the captive lizards were consistently had a high abundance of Bacteroides compared with the wild samples (Figure 7). Petrimonas and Microbacter are fermenters isolated from environment (Grabowski et al., 2005; Sánchez-Andrea et al., 2014). The Proteiniclasticum, a genus of Clostridiaceae, is known as polysaccharide degrader in the gut (Zhang et al., 2010; Wust et al., 2011). A high relative abundance of Clostridiaceae was also observed in Anolis lizards (Hong et al., 2011). Proteiniclasticum may facilitate energy consumption in crocodile lizards. Notably, the diet of earthworm resulted in enrichment of Fusobacterium, which has been reported as human and animal pathogens (Signat et al., 2011). In addition, Fusobacterium has been commonly isolated from infected reptiles (Stewart, 1990). In the gut microbiota of other lizards, the abundance of Fusobacterium, if any, was very low, similar to that of the wild group of crocodile lizards (Martin et al., 2010; Hong et al., 2011; Kohl et al., 2016, 2017; Ren et al., 2016).

The loach-based diet resulted in the significant predominance of Elizabethkingia, Halomonas, Morganella, and Salmonella in the crocodile lizard gut microbiota. These genera are pathogens or opportunistic pathogens in animals and/or humans and have also been found in disease cases in reptiles (O’Hara et al., 2000; Bernardet et al., 2006; Stevens et al., 2009; Kim et al., 2013; Yeo et al., 2016). Particularly, attention should be paid to Salmonella, which causes enteritis and typhoid fever in mammalian and avian species. The Salmonella was detected in all samples, and its relative abundance in loach-fed crocodile lizards reached up to an average of 11.42% (Figure 8). Moreover, the zoonotic potential of Salmonella in reptiles has been widely documented (Murphy and Oshin, 2015). In addition, the Salmonella isolated from the focus of crocodile lizard with disease A caused the death of the Chinese skink in our lab (data have not been published). Moreover, Elizabethkingia was isolated from the focus of crocodile lizard with disease B (data have not been published). Therefore, hand washing should be recommended after contact with crocodile lizards, especially for those who contact these animals frequently. Although these genera have been detected as normal flora in many reptiles, there are many reports that associate these bacteria with reptile diseases such as bacterial pneumonia, osteomyelitis, septicemia, and hepatitis (Huchzermeyer, 1991; Ramsay et al., 2002; Grupka et al., 2006; Chinnadurai and Devoe, 2009). The high content of potential pathogenic bacteria in loach-fed group suggests that the diet of loach not only altered the structure of gut microbiota but also increased the risk of infection for crocodile lizards.

Conclusion

The composition of the crocodile lizard gut microbiota is unique compared with other animals. Diets altered the bacterial community richness and the relative abundance of certain bacteria in the intestine. The gut microbiota of loach-fed crocodile lizard was significantly different from the gut microbiota of the wild and the earthworm-fed crocodile lizards. The earthworm-fed crocodile lizards had a higher abundance of Fusobacterium in the gut compared with the wild lizards. The intestine of loach-fed crocodile lizard was enriched in Elizabethkingia, Halomonas, Morganella, and Salmonella. These bacteria were reported to be pathogens or opportunistic pathogens in human or other animals. This may be a consequence of unbalanced nutrition, as the crocodile lizards in this study were routinely fed with only earthworms or loaches at the Nature Reserves. It seems that the diet of loach was not suitable for crocodile lizards. However, there is no sufficient evidence that the gut microbiota contributes to either disease A or disease B.

This study provides an overview of the gut microbiota of the crocodile lizard, an extremely endangered lizard, at different states as well as the first examination of the relationship between disease and the gut microbiota in lizards. These findings have numerous implications for the practice of crocodile lizard conservation, from the perspectives of both captivity and disease prevention. For instance, the results emphasize that a more diverse diet could improve the care of crocodile lizards in wildlife rescue centers and nature reserves. In addition, the diet of earthworm is better for crocodile lizard than the diet of loach.

Ethics Statement

This study was carried out in accordance with the recommendations of Guidelines of Animal Experiments, the Committee on the Ethics of Animal Experiments of the Guangdong Institute of Applied Biological Resources. The protocol was approved by the Committee on the Ethics of Animal Experiments of the Guangdong Institute of Applied Biological Resources.

Author Contributions

JPC and RCH designed the research. HYJ, NH, HYL, SYL, and ZJW collected the samples. HYJ, LML, and XJZ conducted the research. HYJ, JEM, and JL analyzed the data. HYJ and JPC wrote the manuscript. All authors approved the final version of the manuscript.

Funding

This project was supported by the Training Fund of Guangdong Institute of Applied Biological Resources for PhDs, Masters, and Postdoctoral Researchers (No. GIABR-pyjj201604), the Planning Funds of Science and Technology of Guangdong Province (2016B070701016), the Funds for Environment Construction and Capacity Building of GDAS’ Research Platform (2016GDASPT-0107), GDAS Special Project of Science and Technology Development (2017GDASCX-0107), and Special Funds for Forestry Development and Protection of Guangdong Province (2017).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The reviewers AM and KK, and handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.02073/full#supplementary-material

References

Bano, N., deRae Smith, A., Bennett, W., Vasquez, L., and Hollibaugh, J. T. (2007). Dominance of Mycoplasma in the guts of the long-jawed mudsucker, Gillichthys mirabilis, from five California salt marshes. Environ. Microbiol. 9, 2636–2641. doi: 10.1111/j.1462-2920.2007.01381.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardet, J. F., Hugo, C., and Bruun, B. (2006). The Genera Chryseobacterium and Elizabethkingia. New York, NY: Springer.

Google Scholar

Bletz, M. C., Goedbloed, D. J., Sanchez, E., Reinhardt, T., Tebbe, C. C., Bhuju, S., et al. (2016). Amphibian gut microbiota shifts differentially in community structure but converges on habitat-specific predicted functions. Nat. Commun. 7:13699. doi: 10.1038/ncomms13699

PubMed Abstract | CrossRef Full Text | Google Scholar

Boursier, J., Mueller, O., Barret, M., Machado, M., Fizanne, L., Araujo-Perez, F., et al. (2016). The severity of nonalcoholic fatty liver disease is associated with gut dysbiosis and shift in the metabolic function of the gut microbiota. Hepatology 63, 764–775. doi: 10.1002/hep.28356

PubMed Abstract | CrossRef Full Text | Google Scholar

Brucker, R. M., and Bordenstein, S. R. (2013). The hologenomic basis of speciation: gut bacteria cause hybrid lethality in the genus Nasonia. Science 341, 667–669. doi: 10.1126/science.1240659

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Chinnadurai, S. K., and Devoe, R. S. (2009). Selected infectious diseases of reptiles. Vet. Clin. North Am. Exot. Anim. Pract. 12, 583–596. doi: 10.1016/j.cvex.2009.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Clemente, J. C., Ursell, L. K., Parfrey, L. W., and Knight, R. (2012). The impact of the gut microbiota on human health: an integrative view. Cell 148, 1258–1270. doi: 10.1016/j.cell.2012.01.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Colston, T. J., Noonan, B. P., and Jackson, C. R. (2015). Phylogenetic analysis of bacterial communities in different regions of the gastrointestinal tract of Agkistrodon piscivorus, the cottonmouth snake. PLOS ONE 10:e0128793. doi: 10.1371/journal.pone.0128793

PubMed Abstract | CrossRef Full Text | Google Scholar

Costello, E. K., Gordon, J. I., Secor, S. M., and Knight, R. (2010). Postprandial remodeling of the gut microbiota in Burmese pythons. ISME J. 4, 1375–1385. doi: 10.1038/ismej.2010.71

PubMed Abstract | CrossRef Full Text | Google Scholar

Cryan, J. F., and Dinan, T. G. (2012). Mind-altering microorganisms: the impact of the gut microbiota on brain and behaviour. Nat. Rev. Neurosci. 13, 701–712. doi: 10.1038/nrn3346

PubMed Abstract | CrossRef Full Text | Google Scholar

David, L. A., Maurice, C. F., Carmody, R. N., Gootenberg, D. B., Button, J. E., Wolfe, B. E., et al. (2014). Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559–563. doi: 10.1038/nature12820

PubMed Abstract | CrossRef Full Text | Google Scholar

Gajardo, K., Rodiles, A., Kortner, T. M., Krogdahl, A., Bakke, A. M., Merrifield, D. L., et al. (2016). A high-resolution map of the gut microbiota in Atlantic salmon (Salmo salar): a basis for comparative gut microbial research. Sci. Rep. 6:30893. doi: 10.1038/srep30893

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodrich, J. K., Waters, J. L., Poole, A. C., Sutter, J. L., Koren, O., Blekhman, R., et al. (2014). Human genetics shape the gut microbiome. Cell 159, 789–799. doi: 10.1016/j.cell.2014.09.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabowski, A., Tindall, B. J., Bardin, V., Blanchet, D., and Jeanthon, C. (2005). Petrimonas sulfuriphila gen. nov., sp. nov., a mesophilic fermentative bacterium isolated from a biodegraded oil reservoir. Int. J. Syst. Evol. Microbiol. 55, 1113–1121. doi: 10.1099/ijs.0.63426-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Grupka, L. M., Ramsay, E. C., and Bemis, D. A. (2006). Salmonella surveillance in a collection of rattlesnakes (Crotalus spp.). J. Zoo Wildl. Med. 37, 306–312. doi: 10.1638/05-059.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hird, S. M., Sanchez, C., Carstens, B. C., and Brumfield, R. T. (2015). Comparative gut microbiota of 59 neotropical bird Species. Front. Microbiol. 6:1403. doi: 10.3389/fmicb.2015.01403

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, P.-Y., Wheeler, E., Cann, I. K., and Mackie, R. I. (2011). Phylogenetic analysis of the fecal microbial community in herbivorous land and marine iguanas of the Galápagos Islands using 16S rRNA-based pyrosequencing. ISME J. 5, 1461–1470. doi: 10.1038/ismej.2011.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Yu, H., Wu, Z., Li, Y., Wei, F., and Gong, M. (2008). Population and conservation strategies for the Chinese crocodile lizard (Shinisaurus crocodilurus) in China. Anim. Biodivers. Conserv. 31, 63–70. doi: 10.1371/journal.pone.0091570

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S., and Zhang, H. (2013). The impact of environmental heterogeneity and life stage on the hindgut microbiota of Holotrichia parallela larvae (Coleoptera: Scarabaeidae). PLOS ONE 8:e57169. doi: 10.1371/journal.pone.0057169

PubMed Abstract | CrossRef Full Text | Google Scholar

Huchzermeyer, K. D. (1991). Treatment and control of an outbreak of salmonellosis in hatchling Nile crocodiles (Crocodylus niloticus). J. S. Afr. Vet. Assoc. 62, 23–25.

PubMed Abstract | Google Scholar

Keenan, S. W., Engel, A. S., and Elsey, R. M. (2013). The alligator gut microbiome and implications for archosaur symbioses. Sci. Rep. 3:2877. doi: 10.1038/srep02877

PubMed Abstract | CrossRef Full Text | Google Scholar

Khalifa, M. M., Sharaf, R. R., and Aziz, R. K. (2010). Helicobacter pylori: a poor man’s gut pathogen? Gut. Pathog. 2:2. doi: 10.1186/1757-4749-2-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, K. K., Lee, J. S., and Stevens, D. A. (2013). Microbiology and epidemiology of Halomonas species. Future Microbiol. 8, 1559–1573. doi: 10.2217/fmb.13.108

PubMed Abstract | CrossRef Full Text | Google Scholar

King, G. M., Judd, C., Kuske, C. R., and Smith, C. (2012). Analysis of stomach and gut microbiomes of the eastern oyster (Crassostrea virginica) from coastal Louisiana, USA. PLOS ONE 7:e51475. doi: 10.1371/journal.pone.0051475

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohl, K. D., Brun, A., Magallanes, M., Brinkerhoff, J., Laspiur, A., Acosta, J. C., et al. (2016). Physiological and microbial adjustments to diet quality permit facultative herbivory in an omnivorous lizard. J. Exp. Biol. 219(Pt 12), 1903–1912. doi: 10.1242/jeb.138370

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohl, K. D., Brun, A., Magallanes, M., Brinkerhoff, J., Laspiur, A., Acosta, J. C., et al. (2017). Gut microbial ecology of lizards: insights into diversity in the wild, effects of captivity, variation across gut regions and transmission. Mol. Ecol. 26, 1175–1189. doi: 10.1111/mec.13921

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohl, K. D., and Dearing, M. D. (2014). Wild- caught rodents retain a majority of their natural gut microbiota upon entrance into captivity. Environ. Microbiol. Rep. 6, 191–195. doi: 10.1111/1758-2229.12118

PubMed Abstract | CrossRef Full Text | Google Scholar

Lagier, J. C., Armougom, F., Million, M., Hugon, P., Pagnier, I., Robert, C., et al. (2012). Microbial culturomics: paradigm shift in the human gut microbiome study. Clin. Microbiol. Infect. 18, 1185–1193. doi: 10.1111/1469-0691.12023

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. O., Gilman, F. R., and Weiss, S. L. (2010). Sex-specific asymmetry within the cloacal microbiota of the striped plateau lizard, Sceloporus virgatus. Symbiosis 51, 97–105. doi: 10.1007/s13199-010-0078-y

CrossRef Full Text | Google Scholar

Martin, R., Miquel, S., Langella, P., and Bermudez-Humaran, L. G. (2014). The role of metagenomics in understanding the human microbiome in health and disease. Virulence 5, 413–423. doi: 10.4161/viru.27864

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaughlin, R. W., Chen, M., Zheng, J., Zhao, Q., and Wang, D. (2012). Analysis of the bacterial diversity in the fecal material of the endangered Yangtze finless porpoise, Neophocaena phocaenoides asiaeorientalis. Mol. Biol. Rep. 39, 5669–5676. doi: 10.1007/s11033-011-1375-0

PubMed Abstract | CrossRef Full Text | Google Scholar

McLaughlin, R. W., Cochran, P. A., and Dowd, S. E. (2015). Metagenomic analysis of the gut microbiota of the Timber Rattlesnake, Crotalus horridus. Mol. Biol. Rep. 42, 1187–1195. doi: 10.1007/s11033-015-3854-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, D., and Oshin, F. (2015). Reptile-associated salmonellosis in children aged under 5 years in South West England. Arch. Dis. Child. 100, 364–365. doi: 10.1136/archdischild-2014-306134

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, D. M., Cann, I. K., Altermann, E., and Mackie, R. I. (2010). Phylogenetic evidence for lateral gene transfer in the intestine of marine iguanas. PLOS ONE 5:e10785. doi: 10.1371/journal.pone.0010785

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, T. M., Rogers, T. L., Carlini, A. R., and Brown, M. V. (2013). Diet and phylogeny shape the gut microbiota of Antarctic seals: a comparison of wild and captive animals. Environ. Microbiol. 15, 1132–1145. doi: 10.1111/1462-2920.12022

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, T. Q., Hamilton, P., and Ziegler, T. (2014). Shinisaurus crocodilurus. The IUCN Red List of Threatened Species 2014. Available at: http://dx.doi.org/10.2305/IUCN.UK.2014-1.RLTS.T57287221A57287235.en. doi: 10.2305/IUCN.UK.2014-1.RLTS.T57287221A57287235.en

CrossRef Full Text | Google Scholar

Nguyen, T. Q., and Ziegler, T. (2015). Last chance to See? A review of the threats to and use of the crocodile lizard. Seizure Prosecution 27:19.

Google Scholar

O’Hara, C. M., Brenner, F. W., and Miller, J. M. (2000). Classification, identification, and clinical significance of Proteus, Providencia, and Morganella. Clin. Microbiol. Rev. 13, 534–546. doi: 10.1128/CMR.13.4.534-546.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Ossiboff, R. J., Raphael, B. L., Ammazzalorso, A. D., Seimon, T. A., Niederriter, H., Zarate, B., et al. (2015). A Mycoplasma species of Emydidae turtles in the Northeastern USA. J. Wildl. Dis. 51, 466–470. doi: 10.7589/2014-04-086

PubMed Abstract | CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramakrishna, B. S. (2013). Role of the gut microbiota in human nutrition and metabolism. J. Gastroenterol. Hepatol. 28(Suppl. 4), 9–17. doi: 10.1111/jgh.12294

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramsay, E. C., Daniel, G. B., Tryon, B. W., Merryman, J. I., Morris, P. J., and Bemis, D. A. (2002). Osteomyelitis associated with Salmonella enterica SS arizonae in a colony of ridgenose rattlesnakes (Crotalus willardi). J. Zoo Wildl. Med. 33, 301–310.

PubMed Abstract | Google Scholar

Ren, T., Kahrl, A. F., Wu, M., and Cox, R. M. (2016). Does adaptive radiation of a host lineage promote ecological diversity of its bacterial communities? A test using gut microbiota of Anolis lizards. Mol. Ecol. 25, 4793–4804. doi: 10.1111/mec.13796

PubMed Abstract | CrossRef Full Text | Google Scholar

Sánchez-Andrea, I., Sanz, J. L., and Stams, A. J. (2014). Microbacter margulisiae gen. nov., sp. nov., a propionigenic bacterium isolated from sediments of an acid rock drainage pond. Int. J. Syst. Evol. Microbiol. 64, 3936–3942. doi: 10.1099/ijs.0.066241-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | CrossRef Full Text | Google Scholar

Signat, B., Roques, C., Poulet, P., and Duffaut, D. (2011). Role of Fusobacterium nucleatum in periodontal health and disease. Curr. Issues Mol. Biol. 13, 25–35.

Google Scholar

Stevens, D. A., Hamilton, J. R., Johnson, N., Kim, K. K., and Lee, J. S. (2009). Halomonas, a newly recognized human pathogen causing infections and contamination in a dialysis center: three new species. Medicine 88, 244–249. doi: 10.1097/MD.0b013e3181aede29

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, J. S. (1990). Anaerobic bacterial-infections in reptiles. J. Zoo Wildl. Med. 21, 180–184.

Google Scholar

Sylvia, K. E., Jewell, C. P., Rendon, N. M., St John, E. A., and Demas, G. E. (2017). Sex-specific modulation of the gut microbiome and behavior in Siberian hamsters. Brain Behav. Immun. 60, 51–62. doi: 10.1016/j.bbi.2016.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Thaiss, C. A., Zmora, N., Levy, M., and Elinav, E. (2016). The microbiome and innate immunity. Nature 535, 65–74. doi: 10.1038/nature18847

PubMed Abstract | CrossRef Full Text | Google Scholar

Theodorakopoulos, N., Bachar, D., Christen, R., Alain, K., and Chapon, V. (2013). Exploration of deinococcus-Thermus molecular diversity by novel group-specific PCR primers. Microbiologyopen 2, 862–872. doi: 10.1002/mbo3.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Uetz, P., Hošek, J., and Hallermann, J. (2016). The Reptile Database. Available at: http://www.reptile-database.org

Google Scholar

van Schingen, M., Duc Le, M., Thi Ngo, H., The Pham, C., Quy Ha, Q., Quang Nguyen, T., et al. (2016). Is there more than one Crocodile Lizard? An integrative taxonomic approach reveals Vietnamese and Chinese Shinisaurus crocodilurus represent separate conservation and taxonomic units. Der Zoologische Garten 85, 240–260. doi: 10.1016/j.zoolgart.2016.06.001

CrossRef Full Text | Google Scholar

van Schingen, M., Ihlow, F., Nguyen, T. Q., Ziegler, T., Bonkowski, M., Wu, Z., et al. (2014). Potential distribution and effectiveness of the protected area network for the crocodile lizard, Shinisaurus crocodilurus (Reptilia: Squamata: Sauria). Salamandra 50, 71–76.

Google Scholar

Waite, D. W., and Taylor, M. W. (2015). Exploring the avian gut microbiota: current trends and future directions. Front. Microbiol. 6:673. doi: 10.3389/fmicb.2015.00673

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Wust, P. K., Horn, M. A., and Drake, H. L. (2011). Clostridiaceae and Enterobacteriaceae as active fermenters in earthworm gut content. ISME J. 5, 92–106. doi: 10.1038/ismej.2010.99

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeo, S. H., Kwak, J. H., Kim, Y. U., Lee, J. S., Kim, H. J., Park, K. H., et al. (2016). Peritoneal dialysis-related peritonitis due to Halomonas hamiltonii: a first case report. Medicine 95:e5424. doi: 10.1097/MD.0000000000005424

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, M. L., Dean, S. H., Longo, A. V., Rothermel, B. B., Tuberville, T. D., and Zamudio, K. R. (2015). Kinship, inbreeding and fine-scale spatial structure influence gut microbiota in a hindgut-fermenting tortoise. Mol. Ecol. 24, 2521–2536. doi: 10.1111/mec.13169

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, K., Song, L., and Dong, X. (2010). Proteiniclasticum ruminis gen. nov., sp. nov., a strictly anaerobic proteolytic bacterium isolated from yak rumen. Int. J. Syst. Evol. Microbiol. 60, 2221–2225. doi: 10.1099/ijs.0.011759-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zimmerman, L. M., Vogel, L. A., and Bowden, R. M. (2010). Understanding the vertebrate immune system: insights from the reptilian perspective. J. Exp. Biol. 213, 661–671. doi: 10.1242/jeb.038315

PubMed Abstract | CrossRef Full Text | Google Scholar