Research article

Characterisation and comparison of the mucosa-associated bacterial communities across the gastrointestinal tract of stranded green turtles, Chelonia mydas

  • Received: 18 July 2020 Accepted: 08 October 2020 Published: 19 October 2020
  • Chelonia mydas are primarily herbivorous long-distance migratory sea turtles that contribute to marine ecosystems. Extensive research has been conducted to restore the populations of green turtles. Little is known about their gut microbiota which plays a vital role in their health. We investigated the mucosa-associated bacterial communities across the gastrointestinal (GI) tract of a total four (3, juvenile and 1, adult) stranded green turtles. Samples taken from four GI regions including oesophagus, stomach, small intestine and large intestine were analysed by high-throughput sequencing targeting hypervariable V1-V3 regions of the bacterial 16S rRNA gene. Bacterial diversity and richness decreased longitudinally along the GI tract from oesophagus to the small intestine of stranded turtles. The large intestine showed a higher bacterial diversity and richness compared to small intestine. The bacterial community of green turtles' GI tract was largely dominated by Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes and Fusobacteria. Aerobic and facultative anaerobic bacteria prevailed primarily in the oesophagus while anaerobes (Lachnospiraceae, Peptostreptococcaceae and Ruminococcaceae) constituted the bulk of large intestinal microbiota. Firmicutes dominated the GI tract except within the small intestine where Proteobacteria prevailed. At the OTU level, six percent of the total OTUs (>1% relative abundance) were common in all GI regions. This is a comprehensive characterisation of bacterial microbiota across the GI tract in green turtles which will provide a reference for future studies on turtle gut microbiome and their metabolism to improve their health and nutrition during rehabilitation.

    Citation: Mohammad Shamim Ahasan, Thomas B. Waltzek, Leigh Owens, Ellen Ariel. Characterisation and comparison of the mucosa-associated bacterial communities across the gastrointestinal tract of stranded green turtles, Chelonia mydas[J]. AIMS Microbiology, 2020, 6(4): 361-378. doi: 10.3934/microbiol.2020022

    Related Papers:

    [1] Javier Fernández, Sául Redondo-Blanco, Elisa M. Miguélez, Claudio J. Villar, Alfonso Clemente, Felipe Lombó . Healthy effects of prebiotics and their metabolites against intestinal diseases and colorectal cancer. AIMS Microbiology, 2015, 1(1): 48-71. doi: 10.3934/microbiol.2015.1.48
    [2] Makayla Schissel, Rebecca Best, Shelby Liesemeyer, Yuan-De Tan, Darby J. Carlson, Julie J. Shaffer, Nagavardhini Avuthu, Chittibabu Guda, Kimberly A. Carlson . Effect of Nora virus infection on native gut bacterial communities of Drosophila melanogaster. AIMS Microbiology, 2021, 7(2): 216-237. doi: 10.3934/microbiol.2021014
    [3] Souliphone Sivixay, Gaowa Bai, Takeshi Tsuruta, Naoki Nishino . Cecum microbiota in rats fed soy, milk, meat, fish, and egg proteins with prebiotic oligosaccharides. AIMS Microbiology, 2021, 7(1): 1-12. doi: 10.3934/microbiol.2021001
    [4] Alexandra S McHarg, Steven Leach . The role of the gut microbiome in paediatric irritable bowel syndrome. AIMS Microbiology, 2022, 8(4): 454-469. doi: 10.3934/microbiol.2022030
    [5] Esther Mezhibovsky, Yue Wu, Fiona G. Bawagan, Kevin M. Tveter, Samantha Szeto, Diana Roopchand . Impact of grape polyphenols on Akkermansia muciniphila and the gut barrier. AIMS Microbiology, 2022, 8(4): 544-565. doi: 10.3934/microbiol.2022035
    [6] Stephen T. Abedon . Active bacteriophage biocontrol and therapy on sub-millimeter scales towards removal of unwanted bacteria from foods and microbiomes. AIMS Microbiology, 2017, 3(3): 649-688. doi: 10.3934/microbiol.2017.3.649
    [7] Tanim Jabid Hossain, Mukta Das, Ferdausi Ali, Sumaiya Islam Chowdhury, Subrina Akter Zedny . Substrate preferences, phylogenetic and biochemical properties of proteolytic bacteria present in the digestive tract of Nile tilapia (Oreochromis niloticus). AIMS Microbiology, 2021, 7(4): 528-545. doi: 10.3934/microbiol.2021032
    [8] Paola Brun . The profiles of dysbiotic microbial communities. AIMS Microbiology, 2019, 5(1): 87-101. doi: 10.3934/microbiol.2019.1.87
    [9] Chunze Zhang, Mingqian Ma, Zhenying Zhao, Zhiqiang Feng, Tianhao Chu, Yijia Wang, Jun Liu, Xuehua Wan . Gut mucosal microbiota profiles linked to development of positional-specific human colorectal cancer. AIMS Microbiology, 2024, 10(4): 812-832. doi: 10.3934/microbiol.2024035
    [10] Dawn B. Goldsmith, Zoe A. Pratte, Christina A. Kellogg, Sara E. Snader, Koty H. Sharp . Stability of temperate coral Astrangia poculata microbiome is reflected across different sequencing methodologies. AIMS Microbiology, 2019, 5(1): 62-76. doi: 10.3934/microbiol.2019.1.62
  • Chelonia mydas are primarily herbivorous long-distance migratory sea turtles that contribute to marine ecosystems. Extensive research has been conducted to restore the populations of green turtles. Little is known about their gut microbiota which plays a vital role in their health. We investigated the mucosa-associated bacterial communities across the gastrointestinal (GI) tract of a total four (3, juvenile and 1, adult) stranded green turtles. Samples taken from four GI regions including oesophagus, stomach, small intestine and large intestine were analysed by high-throughput sequencing targeting hypervariable V1-V3 regions of the bacterial 16S rRNA gene. Bacterial diversity and richness decreased longitudinally along the GI tract from oesophagus to the small intestine of stranded turtles. The large intestine showed a higher bacterial diversity and richness compared to small intestine. The bacterial community of green turtles' GI tract was largely dominated by Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes and Fusobacteria. Aerobic and facultative anaerobic bacteria prevailed primarily in the oesophagus while anaerobes (Lachnospiraceae, Peptostreptococcaceae and Ruminococcaceae) constituted the bulk of large intestinal microbiota. Firmicutes dominated the GI tract except within the small intestine where Proteobacteria prevailed. At the OTU level, six percent of the total OTUs (>1% relative abundance) were common in all GI regions. This is a comprehensive characterisation of bacterial microbiota across the GI tract in green turtles which will provide a reference for future studies on turtle gut microbiome and their metabolism to improve their health and nutrition during rehabilitation.


    Microbial communities inhabiting the gastrointestinal (GI) tract comprise a complex ecosystem and play a vital role in maintaining host physiology, ranging from metabolic activity to host immune homeostasis [1],[2]. They are also involved in many functions that are not encoded in the host‘s DNA [3]. In recent decades, several studies on gut microbiota have confirmed that a balanced gut microbiome is essential for the host's ability to maintain a healthy state. Perturbations in the stability of gut microbial communities dispose the host to pathogenic invasions which may lead to several GI diseases and disorders [4][7]. Knowledge of microbial community and population dynamics in symbiosis, as well as in dysbiosis, are essential for developing management strategies for treating GI associated diseases and disorders.

    Recent advancement in the molecular biology and computational techniques have empowered the scientific community to explore several influencing factors on the establishment and maintenance of host gut microbiome. Studies revealed that diet is one of the primary drivers of functional capacity within the GI tract, resulting in a convergence of bacterial communities and the phylogenetically un-related host [8],[9]. Besides diet, several host factors such as physiology, gut structure and genetics can also shape gut microbial communities [10],[11]. Studies in terrestrial and aquatic animals revealed that several environmental factors such as habitat, temperature and salinity can influence the microbial community compositions within the GI tract of the host [12][14].

    To date, large scale research on host gut microbiome has been conducted on several aquatic and terrestrial mammals [9],[15],[16], fishes [17], birds [18][20] and to a lesser extent invertebrates [21][23]. Few studies have investigated marine reptiles, including the endangered green turtles [24][26]. Green sea turtles are a long-distance migrators, hind gut fermenters, which forages primarily on sea grasses [27]. The gut microbiota of green turtles are believed to play a crucial role in gaining energy from their food sources [28],[29]. Moreover, they can also contribute to several other aspects of their health and development of disease, as noted in other animals [15],[30],[31]. Therefore, detailed understanding of resident gut bacterial communities along the digestive tract of green turtles is very important to develop strategies to treat gut associated disorders and restore the host's normal gut microbiome during rehabilitation.

    To date, most investigations of the green turtle's gut microbiome have typically involved bacterial identification in cloacal or faecal contents. The faecal bacterial communities of green turtle mostly dominated by bacteria with in the phyla Firmicutes, Bacteroidetes, Proteobacteria and Actinobacteria while the relative abundance of these bacteria may vary at different health [32], age [24],[25], nutrition [33],[34] and environmental conditions [25],[35]. Additionally, bacteria present in faecal contents only represent the overall gut microbiota rather than bacterial communities of specific anatomical regions of the GI tract [36],[37]. Recently, an investigation of two Hawaiian green turtles' gut microbiota was carried out by in situ sampling of gut contents from hindguts (caecum, large intestine and rectum) of the GI tracts and reported more than half of the sequences belonged to the phylum Firmicutes, Bacteroidetes and Proteobacteria while no consistent pattern of bacterial diversity across hindgut was recorded [26]. It is hypothesized that microbiota present in faecal or gut contents represent a pool of both resident and transient microbial population in the GI tract [38]. The mucosal microbiota is more likely to be resident and have a higher potential to influence the host than transient microbial populations and is therefore of more interest to intestinal health [39]. In this study, we therefore, aimed to characterise and compare the mucosa-associated bacterial communities across different regions of the GI tract of stranded green turtles using high-throughput sequencing analysis.

    Stranded green turtles are occasionally cared for in rehabilitation centers until they recover. In this study, samples were taken from a total of four stranded green turtles that were collected from the Whitsunday beach of the central Great Barrier Reef (Table S1). These turtles died shortly (1 to 4 days) after arriving for rehabilitation. During rehabilitation, turtles were offered human grade squid (Loligo opalescens), but were not interested in feeding. The freshly dead turtles were kept frozen (below −20 °C) for subsequent post-mortem examination. The metadata, curved carapace length (CCL) and body weight were recorded in accordance to DEHP (Department of Environment and Heritage Protection) standard operating procedures upon the turtle's arrival to the rehabilitation center [41] (Table S1). Turtle 1 was underweight and emaciated. Turtle 2 was determined to have a healed boat-strike scar on the carapace. No visible traumatic injuries were recorded for turtle 3 and turtle 4 however, they were emaciated and infested by barnacles.

    During the necropsies, samples were collected from four different regions of the GI tract, including the oesophagus, stomach, small intestine (duodenum) and large intestine (caecum) of the animals. Using a sterile scalpel blade, the luminal cavity of each GI region was opened and the gut contents were separated from the mucosal wall. A deep mucosal scraping sample was collected using a sterile cotton swab rubbing thoroughly against the outer mucosal layer of the gut lumen with minimum contamination by gut contents. Swabs were placed in microfuge tubes, transported on ice within 4–6 hrs directly to the laboratory at James Cook University (Townsville, QLD), and frozen at −80 °C until processing. This study was conducted under permit from James Cook University Animal Ethic Committee (Permit no. A2101), Department of Environment and Heritage Protection Authority (permit no. WISP15015914) and Great Barrier Reef Marine Park Authority (Permit no. G14/37285.1).

    The PowerLyzer® PowerSoil® DNA Isolation Kit was used to extract bacterial nucleic acid from mucosal swab samples following the manufacturer's standard protocol (MO BIO Laboratories, Inc., Carlsbad, CA, USA). Once extracted, the quantity of the extracted DNA was determined using GeneQuant Prospectrophotometer (Amersham Pharmacia Biotech) and stored at −20 °C until used.

    Genomic DNA extracted from mucosal swab samples were used as templates for PCR amplification of the hypervariable V1-V3 regions of the bacterial 16S rRNA gene using the forward (27F 5′-AGAGTTTGATCMTGGCTCAG-3′) and reverse (519R 5′-GWATTACCGCGGCKGCTG-3′) primers, and required Illumina flowcell adaptor sequences [42]. The PCR reaction was carried out using an initial denaturation at 95 °C for 7 min, followed by 29 cycles of denaturation at 94 °C for 45 second, primer annealing at 50 °C for 1 minute and extension at 72 °C for 1 minute, with a final elongation at 72 °C for 7 minutes. The AmpliTaq Gold 360 mastermix (Life Technologies, Australia) was used for the primary PCR and a secondary PCR to index the amplicons was performed with TaKaRa Taq DNA Polymerase (Clontech). The resulting amplicons were measured by fluorometry (Invitrogen Picogreen) and normalized. The eqimolar pool was then estimated using qPCR (KAPA), followed by sequencing on the Illumina MiSeq platform (San Diego, CA, USA). Sequencing and library preparation were done by the Australian Genome Research Facility, Brisbane, Australia.

    All bioinformatics analysis were performed using the software package ‘The Quantitative Insights Into Microbial Ecology (QIIME v1.9.1)’ [42]. A unique barcode was used to label each sample and Illumina sequencing resulted in 17,22,641 paired-end reads from a total of 16 samples. The raw paired-ends Fastq files were quality trimmed (q < 20) to remove low-quality and single end reads using Sickle (version 1.33) (https://github.com/najoshi/sickle). The forward and reverse reads were assembled using PANDAseq Assembler [43]. Using USEARCH (version 8.0.1623), sequences were quality filtered to remove the chimeric reads and full length duplicate sequences [44]. Sequences were then clustered into operational taxonomic units (OTUs) using UCLUST taxonomy assigner in QIIME. To obtain the number of reads in each OTU, reads were mapped back to OTUs with a minimum identity of 97%. The reads that failed to hit the reference database were subsequently clustered as de novo. Using QIIME, taxonomy was assigned using Silva database (version 128, Sep 2016) (https://www.mothur.org/wiki/Silva_reference_files#Release_128). The sequences that were unable to assign to any known taxa were categorized as ‘Unclassified’.

    A new OTU table was constructed considering OTUs that showed a relative abundance > 1% and samples that retained > 1000 sequence reads. Bacterial count data were then normalised by total sum scaling (TSS) normalization method. The abundance based coverage estimator (ACE), Chao1 and Shannon diversity indices were used to determine bacterial richness and diversity of samples using R software (version 3.3.0) [45]. One-way ANOVA was used to test for differences in the bacterial diversity/richness among different GI regions. Individually based rarefaction curves were constructed to determine whether sampling yielded sufficient OTU coverage to delineate the gut bacterial composition of each turtle. Furthermore, conditional uncovered probability (CUP) and Good's coverage were estimated for each sample in QIIME. Distances between bacterial communities of different samples of different GI regions were calculated using weighted UniFrac distance matric [46] and further visualised by Principle Coordinate Analysis (PCoA) in Calypso [47]. The analysis of similarity (ANOSIM) was performed by applying Bray-Curtis [48] and weighted UniFrac [46] distance matrices to evaluate the association between bacterial communities of different GI regions. The taxa level differences among bacterial communities of different GI regions were assessed by ANOVA and Wilcoxon signed-rank test. A significant p-value (P < 0.05) indicated if the taxa were significantly different between GI regions. P-values were also adjusted for multiple testing by the false discovery rate (FDR) using Benjamini and Hochberg method [49]. Venn diagrams were constructed to visualise the amount of unique and mutually shared gut bacterial communities in different GI regions.

    A total of 703,799 high-quality sequences obtained after assembly and quality filtering were used for downstream analyses. The mean sequence length ranged from 434 to 500 bp (avg. 468.43 bp). Individual samples were covered by an average 43,987 reads (ranging from 15,125 to 75, 952). The highest number of reads was recorded for the oesophageal samples. All sequences were delineated into OTUs with 97% nucleotide sequence identity threshold, where 23,392 unique OTUs were identified in total and 405 OTUs were retained after filtering the low abundant OTUs (relative abundance < 1%) and OTUs that retained above 1,000 sequences (up to 57,928 sequences). The highest number of OTUs were recorded for large intestinal (LI) samples (123.00 ± 29.16) while a lower number of OTUs (64.25 ± 6.49) were in the small intestinal (SI) samples of green turtles (Table 1). The individually based rarefaction curves were constructed that tended to approach saturation plateau in our current sequencing depth (Figure S1). The mean conditional uncovered probability (PE) for different GI region samples ranged from 0.02 to 0.04 (Table 1) and the Good's coverage was > 92% (ranging from 92 to 99%) for all samples, suggesting the presence of sufficient phylotypes for the current analysis although more phylotypes may be recovered by increasing the sequencing depth (Table 1).

    Table 1.  Alpha diversity metrics for the bacterial communities across different regions of the gastrointestinal (GI) tract at an evolutionary distance D = 0.03. Values indicates means and standard errors derived from total samples for each GI regions. ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine.
    Indices ES ST SI LI
    Sequences 48535.50 ± 10864.05 37154.00 ± 2571.05 24939.25 ± 3098.35 41552.50 ± 10551.54
    OTUs 109.25 ± 19.01 114.00 ± 9.14 64.25 ± 6.49 123.00 ± 29.16
    Chao 1 index 144.80 ± 22.71 139.86 ± 12.77 87.07 ± 11.33a 157.41 ± 31.99
    ACE index 148.61 ± 25.40 138.05 ± 13.10 90.62 ± 11.93b 157.28 ± 28.99
    Shannon index 4.91 ± 0.47 4.33 ± 0.40 2.82 ± 0.63c 4.24 ± 1.01
    Good's coverage 0.97 ± 0.01 0.97 ± 0.01 0.97 ± 0.01 0.96 ± 0.01
    Conditional Uncovered Probability (CUP)
    PE 0.03 ± 0.01 0.02 ± 0.00 0.02 ± 0.01 0.04 ± 0.01
    Lower bound 0.004 ± 0.002 0.003 ± 0.001 0.004 ± 0.001 0.006 ± 0.002
    Upper bound 0.042 ± 0.016 0.028 ± 0.006 0.040 ± 0.014 0.056 ± 0.017

    Note: Superscript ‘a, b’ indicate Chao 1 and ACE index of SI were significantly different from ST indices (ANOVA, P < 0.05). ‘c’ indicates Shannon index of SI was significantly different (ANOVA, P < 0.05) from ES and ST.

     | Show Table
    DownLoad: CSV

    The Chao 1 and abundance-based coverage estimators (ACE) were used to evaluate the community richness of each sample from different regions of GI tract (Table 1). Bacterial richness was lower in the small intestinal samples compared to oesophageal, stomach and large intestinal samples. However, no significant difference (P > 0.2) was found between observed OTUs and actual OTUs, estimated by Chao 1 and ACE estimators for any GI region. Bacterial diversity as estimated by the Shannon index (SI) revealed significant difference between oesophageal and small intestinal samples (t test, P = 0.037). Additionally, the stomach and large intestinal mucosa also harbored a relatively higher level of diverse bacteria than the small intestinal mucosa. However, at this sequencing depth, oesophageal samples showed the highest bacterial diversity (mean SI 4.91 ± 0.47) compared to other samples from different GI regions.

    Although several bacterial communities exhibited a consistent presence across the length of the green turtles' GI tract, the overall composition of mucosa-associated bacterial microbiome was found to be significantly associated with GI regions (ANOSIM R: 0.28, P < 0.05, Bray-Curtis, Figure S2). The strongest association between bacterial microbiomes and GI regions was recorded for the small and large intestinal mucosa-associated bacterial communities (ANOSIM R: 0.469, P < 0.05, Bray-Curtis). PCoA plot and hierarchical dendrograms were constructed based on weighted UniFrac distance metric, clearly demonstrated that the majority of large intestinal samples were clustered together and distinctly separated from small intestinal samples although a higher variation was noticed for small intestinal samples (Figure 1). The oesophageal and stomach samples were clustered closely to one another and clearly separated from small intestinal samples.

    Figure 1.  Principle coordinates analysis (PCoA) of the dissimilarity among samples from different regions of gastrointestinal tract of green turtles. (a) PCoA plots and (b) hierarchical dendrograms were constructed using weighted UniFrac distance matric. T1: turtle 1, T2: turtle 2, T3: turtle 3, T4: turtle 4, ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine.

    Phylum-level analysis of the mucosa-associated bacterial communities of green turtle GI tract revealed a total of 29 bacterial phyla, however, majority (17/29) are lower in abundant (<1%). The majority of the bacterial community present within the GI tract belonged to Firmicutes (57.8%) and Proteobacteria (21.3%), while the rest were distributed amongst Actinobacteria (6.4%), Bacteroidetes (3.6%), Fusobacteria (2.4%), Spirochaetae (0.19%), Saccharibacteria (0.16%), and unclassified bacteria (7.9%) (Figure 2a). Other phyla, such as Thermotogae, Tenericutes, Synergistetes, Gracilibacteria, Chloroflexi, comprised <0.25% of total composition. Oesophageal and stomach mucosa harboured the majority of the identified phyla (24 and 22 of the 29), while the lowest number of bacterial phyla were observed in the small intestinal mucosa. Among the phyla identified in this study, only Firmicutes, Proteobacteria and Actinobacteria were present in samples of all GI regions. These phyla were ubiquitously distributed in the GI tract whilst the relative abundance varied between the gut sites. The phylum Firmicutes was significantly abundant in oesophageal (30.3%), stomach (85.1%) and large intestinal (86.5%) mucosa-associated microbial communities, while a drastically lower abundance was observed for the small intestinal bacterial population (Figure 2b). The relative abundance of Proteobacteria was significantly higher (P < 0.05) in small intestinal microbial communities than other GI regions (Figure 2b). In summation, Proteobacteria showed an inverse profile to Firmicutes and was the second most abundant phylum across the GI tract. The relative abundance of Actinobacteria, Bacteroidetes and Fusobacteria was higher in oesophageal mucosa-associated bacterial communities compared to other GI regions (Figure 2a).

    Figure 2.  Composition of green turtle gut microbiomes. (a) cumulative abundance (%) of the bacterial communities at phylum level across different regions of the gastrointestinal (GI) tract. (b) Relative abundance of two most dominant phyla across the GI tract. T1: turtle 1, T2: turtle 2, T3: turtle 3, T4: turtle 4, ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine. Significantly different regions are shown in barchart using Wilcoxon signed-rank test (* p < 0.05, **: p < 0.01, ***: p < 0.001). Standard error is depicted by error bars.

    A total of 190 families were identified from the complete data set. The 30 most abundant families (Figure 3) belonged to 11 different classes that include Clostridia, Gammaproteobacteria, Bacilli, Actinobacteria, Bacteroidia, Alphaproteobacteria, Fusobacteria, Epsilonproteobacteria, Spirochaetes, Coriobacteria and Erysipelotrichia (Figure 4). Clostridia was one of the most abundant bacterial classes in all anatomic regions, except the small intestine, where a significantly lower (P < 0.05) abundance of Clostridia was observed. Clostridia was mainly comprised of seven families, Peptostreptococcaceae, Lachnospiraceae, Clostridiaceae 1, Ruminococcaceae, Eubacteriaceae, Family XIII and Family XI that accounted for an average 45.6% of the total identified families. Bacteria within the family Peptostreptococcaceae were significantly abundant (P < 0.007) in stomach and large intestinal samples compared to oesophageal and small intestinal samples (Figure 5). A high abundance (P = 0.02) of Lachnospiraceae was observed only in large intestinal samples compared to oesophageal and small intestinal samples. Bacteria within the family Ruminococcaceae was significantly lower (P = 0.02) in the oesophageal mucosa compared to other GI regions. The relative abundance of Gammaproteobacteria was higher in small intestinal samples without any significant difference (P = 0.064) to other GI regions. Gammaproteobacteria includes Aeromonadaceae, Vibrionaceae, Enterobacteriaceae and Xanthomonadaceae, all of which were dominant in small intestinal samples except Xanthomonadaceae, found only in oesophageal samples. The prevalence of Bacilli was lower in the large intestinal mucosa compared to other GI regions. A higher abundance of Bacteroidia (Porphyromonadaceae and Marinilabiliaceae) was observed within the oesophageal and stomach mucosa-associated bacterial communities (Figure 5). Likewise, Alphaproteobacteria (Rhodobacteraceae) and Actinobacteria (Dietziaceae and Propionibacteriaceae) were absent in large intestinal samples, but were higher in abundance in oesophageal samples. The Campylobacteraceae, Helicobacteraceae, Acidaminococcaceae, Desulfobulbaceae, Flavobacteriaceae and Marinilabiaceae significantly prevailed in the oesophageal mucosa-associated bacterial communities. Our analysis revealed that one or two families represented the majority of OTUs and were highly abundant in the turtles. For example, Vibrionaceae 17162, Aeromonas 19906, uncultured Clostridium sp., Peptostreptococcaceae 17846 and Clostridiaceae 1 5902 comprised of 5.3%, 4.6%, 4.5%, 4.2%, 3.3%, and 2.0% of the total sequences respectively. At the genus level, 459 taxa were identified across the GI tract of green turtles while more than 35.8% of the total sequences were not identified at the genus level.

    Figure 3.  The heatmap indicating relative abundance of major bacteria families across different regions of the gastrointestinal tract of green turtles. ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine.
    Figure 4.  Composition (a) and intra-class variations (b) of most abundant bacteria classes across different regions of the gastrointestinal tract of green turtles. T1: turtle 1, T2: turtle 2, T3: turtle 3, T4: turtle 4, ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine.

    We constructed a Venn diagram to determine the unique and mutually shared bacterial taxa present across the GI tract of green turtles. Our results showed that only 11 (5.6%) of the total OTUs (relative abundance > 1%) were common in all gut regions, while 38 (19.2%), 6 (3%), 23 (11.6%) and 57 (28.6%) OTUs were strictly associated with oesophageal, stomach, small and large intestinal bacterial communities respectively (Figure 6a). Moreover, eight of the 11 common OTUs were within the phylum Firmicutes and three OTUs represented Proteobacteria, Bacteroidetes, and Fusobacteria respectively (Table S2). Sixteen (8.1%) of the total OTUs were common between small intestinal and large intestinal bacterial communities. At the family level, 19 (38%) families were shared among all GI regions and 32 (64%) families were exclusively shared between oesophageal and stomach mucosa-associated microbiomes (Figure 6a). Six of the 10 phyla (relative abundance >1%) were distributed across the GI tract of green turtles while three phyla were exclusively present in oesophageal and stomach mucosa-associated bacterial communities (Figure 6a). It is worth noting that inter-turtle variations were also observed from the phylum to OTU levels of the bacterial communities of green turtles (Figure 6b). Only six of the total (405) OTUs were found in all green turtles, where five represented the phylum Firmicutes.

    Figure 5.  The top 10 bacteria families that were significantly different across the gastrointestinal tract of green turtles. ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine. Significantly different regions are shown in barchart using Wilcoxon signed-rank test (* p < 0.05, **: p < 0.01, ***: p < 0.001). Standard error is depicted by error bars.
    Figure 6.  Venn diagrams represent the numbers of total shared and exclusive taxa. (a) Taxa at phylum (i), family (ii), and OTU (iii) level within the different gut regions of green turtles. (b) Taxa at OTU level in each green turtle. T1: turtle 1, T2: turtle 2, T3: turtle 3, T4: turtle 4, ES: Oesophagus, ST: Stomach, SI: Small intestine and LI: Large intestine.

    To date, little is known about the diverse gut microbiota of endangered green turtles [24],[25],[33] in contrast to the extensive knowledge of their ecology, as well as distribution [50],[51]. Knowledge of the microbial community is imperative to understand the host-microbial interaction in both healthy and diseased animals. In this study, we explored the community composition and structure of the mucosa-associated bacterial microbiomes across the length of gastrointestinal (GI) tract of stranded green turtles using high-throughput sequencing and identified an enormous diversity of microbes.

    As green turtles are endangered species, it is close to impossible to obtain permits to kill healthy turtles for research purposes in Australia. The samples for this study were therefore collected from turtles that died within one week of entering rehabilitation and without having consumed the food offered. Although the underlying causes of the strandings may have varied between the 4 turtles and thereby influenced the intestinal microbiota in different ways, there was little variation between samples collected in the four regions of the intestinal tract and we therefore feel that is it appropriate to proceed to make general conclusions based on these samples. Gut bacterial diversity decreased across the longitudinal axis of the green turtle GI tract from oesophagus to small intestine while the large intestine showed a higher bacterial diversity compared to the small intestine but less than stomach. This finding is in accordance with findings in terrestrial herbivores, such as cattle, where bacterial diversity was higher in the stomach compared to small intestinal mucosa [15]. The presence of highly diverse bacterial communities in the oesophageal and stomach mucosa could be due to the existence of transient microbiota. Highly diverse bacterial communities in sea water might influence the diversity of bacterial microbiomes in the oesophageal mucosa of green turtles because turtles are continuously ingesting bacteria associated with their food and water [52],[53]. These findings support the theory speculating that the mouth is an access point for gut microbes (highest diversity) with selection for a subetaoup of the cumulative diversity appearing as these microbes pass through the GI tract. However, studies in laboratory mice revealed that the microbial diversity in stomach and small intestinal mucosa were similar to the large intestinal mucosa [54]. The presence of a diverse bacterial community in the large intestinal mucosa of green turtles indicates the occurrence of a more complex micro-ecosystem that might be associated with nutrient absorption and assimilation [55] and possibly the availability of oxygen [56]. A similar finding was also observed in the colonic mucosal samples of rat and humans [36].

    The GI mucosa-associated microbiota of green turtles were largely dominated by Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, Fusobacteria and a considerable number of unclassified bacteria. However, the detailed composition of these phyla at a lower taxonomic level differed notably in different GI regions. Our results demonstrated that Firmicutes were consistently present across different regions of the GI tract mucosa, however, these bacteria were significantly dominant in stomach and large intestine compared to esophagus and small intestine. The higher prevalence of Firmicutes was supported by findings in other studies on green turtle using cloacal swab or faecal contents samples [24],[25],[57] as faecal swab represents the overall gut microbial population of green turtles. A similar finding was also reported by Mcdermid et. al., 2020 [26] in Hawaiian green turtles where in situ gut contents samples of the hindgut were used for analysis. Additionally, the study suggested to include samples across the GI tract to investigate whether these differences are significant. Members of Firmicutes have been found to be associated with harvesting energy and absorption of nutrients through microbial fermentation of feed components [8],[58]. Studies suggested that the caecum is the initial site for this fermentation activity while little digestion and absorption may take place in the stomach and small intestine [28],[29]. In green turtles, examination of the GI ingesta demonstrated that many of the plant cell materials remained undigested until it reached the caecum where microbial fermentation occurs [28]. Our current study revealed that several families within the phyla Firmicutes, such as Ruminococcaceae, Lachnospiraceae, and Peptostreptococcaceae, had a significant presence in the large intestinal mucosa while Peptostreptococcaceae and Ruminococcaceae also dominated the stomach mucosa. Lachnospiraceae was more abundant than Ruminococcaceae in the large intestine, as reported previously by Campos and her colleagues for the excrements of green turtles from Brazil [57] and by Hong and her colleagues for excrement of marine iguanas [59]. The opposite is true in terrestrial herbivores such as cattle [15], and pigs [39]. It is speculated that the composition of cell wall of macroalgae and terrestrial plants may explain the difference [57]. Juvenile green turtles in eastern Australia forage primarily on macroalgae, while adults preferentially consume seagrass [60]. Green turtles sampled in this study were all juveniles and a macroalgael based diet might explain the prevalence of Lachnospiraceae over Ruminococcaceae in these animals. Bacteria within these families are known to proliferate by hydrolyzing dietary fibre and complex carbohydrates [61] with Ruminococcaceae and Lachnospiraceae being major butyrate producers that reduce leeky gut syndrome, reduce available oxygen in the gut and prevent pathogenic Proteobacteria from dominating the gut [56]. Peptostreptococcaceae can also play a vital role in feed degradation and digestion [62],[63]. However, bacterial community structure and membership within the GI tract can be influenced by factors such as, physiochemical conditions along the GI tract including GI motility, pH and GI secretions. A recent study using humanized mice (germ-free mice colonized with human faecal microbiota) clearly revealed that GI motility can alter the gut microbial composition of the host although it is dependent on diet [64]. Studies in human colonic microbiota demonstrated that a lowering pH from substrate fermentation in the large intestine may increase butyrate production and populations of butyrate-producing bacteria while at the same time inhibiting the Bacteroides spp. [65]. The availability of different substrates and nutrients along the GI tract may also promote distinct bacterial populations in different GI regions of the host [66]. Moreover, the availability of oxygen in different GI regions can also play a vital role in gut bacterial compositions [67]. Oxygen availability in the oesophagus and stomach was higher, therefore, the aerobic and facultative bacteria including Propionibacteriaceae, Desulfobulbaceae, Marinilabiaceae, Campylobacteriaceae, Helicobacteriaceae, Flavobacteriaceae were significantly enriched in oesophageal and stomach mucosa. In contrast, in absence of oxygen, anaerobes or oxygen sensitive bacterial communities such as Lachnospiraceae, Clostridiaceae, Ruminococcaceae were enriched in the large intestinal mucosa of GI tract.

    Our results demonstrate that bacteria within the phylum Proteobacteria was significantly less abundant in oesophageal, stomach and large intestinal mucosa while its prevalence was higher in small intestinal mucosa of the GI tract. The present findings are in support to the findings of other study where bacteria within the phylum Proteobacteria recorded less abundant compared to members of Firmicutes and Bacteroidetes in hindguts (caecum, colon and rectum) of Hawaiian green turtles [26]. In this study, the abundance of Proteobacteria in the small intestinal mucosa was largely due to a higher prevalence of the bacteria belonging to the class Gammaproteobacteria, which includes the Vibrionaceae, Enterobacteriaceae and Aeromonadaceae. The majority of these bacteria are known to be facultative anaerobes and multiply at neutral to slightly basic pH (~7.4) [68]. However, previous studies in green turtles revealed a fairly acidic pH (5.5–6) in the large intestinal contents (different from the other reptiles) and a very acidic pH (3.85–4) in the stomach contents [28]. In contrast, the small intestine exhibited a fairly neutral pH (7.0) which provides a suitable environment for the occurrence of Gammaproteobacteria [28],[68]. Gammaproteobacteria, a physiologically and metabolically diverse class, can play an important role in preparing the gut for successive colonization by strict anaerobes through utilizing oxygen, changing the gut pH, producing CO2 and nutrients [69],[70]. However, an over colonization of these opportunistic bacteria (Vibrionaceae and Enterobacteriaceae) common markers of dysbiosis [56] within the small intestinal mucosa may competitively exclude the other normal flora and result in decreased bacterial diversity in the small intestinal mucosa [71],[72]. Further investigation is suggested to explore these phenomena.

    Actinobacteria was the third most abundant phylum in the oesophageal mucosa. However, they were very low in other GI regions. The higher prevalence of Actinobacteria in the oesophageal mucosa was mainly due to higher abundance of bacteria in the Dietziaceae and Propionibacteriaceae. The members of Dietziaceae are characterized by the presence of mycolic acids [73]. They are mainly chemoorganotrophic and exhibit an oxidation type of metabolism. However, the diversity of the metabolic activities of different members of this family are still unknown which demand better genetic tools to characterize and exploit their diverse functions within the GI tract of animals [73]. Bacteria within the Propionibacteriaceae are well recognized, normal inhabitants of the GI tract for both humans and animals [74],[75]. The pathogenic potential for many Propionibacteria, such as Propionibacterium or Eubacterium, have been claimed for many years and later on, recognized as opportunistic rather than obligate pathogens [74]. These bacteria are often encountered in multiple infection processes where it is difficult to determine which component of the mixed microbial community is a true pathogen and which is only a ‘hanger-on’ [74]. Furthermore, in this study, bacteria within the phylum Bacteroidetes were less abundant in the small and large intestinal mucosa while it was prominent in oesophageal and stomach mucosa and composed of mainly Porphyromonadaceae known members of the buccal cavity associated with the salivary biome [76]. Bacteria within the Porphyromonadaceae were also observed in previous studies on cloacal swabs of green turtles [25]. However, in this study, the reason for lower abundance of Bacteroidetes in the small and large intestinal mucosa is still unclear because the majority of food fermentation of green turtles usually takes place in the caecum [28] and Bacteroidetes species are well known to degrade a variety of plant polysaccharides including pectin, xylan, galactomannan, arabinogalactan, alginate and glucomannan [77][79] and are major butyrate producers. Therefore, further research is recommended to explore the exact reasons of the low level of Bacteroidetes prevalence in small and large intestinal mucosa.

    Furthermore, it is worth noting that the predominance of different phyla in study samples were largely due to the dominance of a small number of distinct OTUs. The top 10 most abundant OTUs, within each GI region, comprised 70.2%, 87.1% and 77.9% of the total resident microbiomes in stomach, small intestinal and large intestinal mucosa respectively. Additionally, four of the top 10 most abundant OTUs of stomach and large intestinal mucosa were common and accounted for 38.4% and 50.4% of the total OTUs. This finding indicates the possibility of similar mucosa-associated microbiota between the stomach and large intestine of green turtles. However, further investigation is suggested to determine the core mucosa-associated gut bacterial community of healthy green turtles.

    We have generated a comprehensive landscape of green turtle bacterial microbiomes across the GI tract in stranded turtles. Our results reveal that mucosa-associated bacterial communities of various compositions occupied different regions of the GI tract. The oesophagus and large intestine exhibit the highest bacterial diversity compared to stomach and small intestine of green turtles. This finding supports the recognized notion that different anatomic regions of the GI tract have their own physiochemical conditions such pH and oxygen, which exert selective pressures on the bacterial communities and play an important role in shaping the microbiota of the GI tract. Interestingly, our study identified several members of opportunistic bacteria across the GI tract which might be overrepresented in stranded turtles. Our results also indicate that faecal samples can only represent the overall bacterial communities rather than microbiota of specific regions of the GI tract. It is therefore advisable to be cautious when interpreting results of faecal samples when investigating gut-associated disease in turtles. Our study provides a helpful reference for further detailed investigation of healthy sea turtle gut microbiomes and their metabolic functions to improve their health and nutrition during rehabilitation.


    Acknowledgments



    The authors thank Eco Barge Clean Seas Inc., Airlie Beach, QLD and Bowen Sea Turtle Assessment and Rehabilitation (BSTAR), QLD for their help with sample collection and necessary information. We thank Dr Jacqueline Picard, JCU and Dr Roger Huerlimann, JCU for useful advice on sample and data processing.

    Conflict of interest



    The authors declare no conflict of interest.

    Author Contributions



    Conceptualization, MSA and EA; Formal analysis, MSA and LO; Funding acquisition, MSA and EA; Investigation, MSA and EA; Methodology, MSA, LO and EA; Software, MSA and TBW; Supervision, LO and EA; Visualization, MSA, TBW, LO and EA; Writing-original draft, MSA; Writing-review & editing, TBA, LO and EA.

    Funding



    The study was supported by the post-graduate research funds from the College of Public Health, Medical and Veterinary Sciences (CPHMVS) under James Cook University, QLD 4811 Australia.

    [1] Chung H, Kasper DL (2010) Microbiota-stimulated immune mechanisms to maintain gut homeostasis. Curr Opin Immunol 22: 455-460. doi: 10.1016/j.coi.2010.06.008
    [2] Dethlefsen L, McFall-Ngai M, Relman DA (2007) An ecological and evolutionary perspective on human-microbe mutualism and disease. Nature 449. doi: 10.1038/nature06245
    [3] Gill SR, Pop M, DeBoy RT, et al. (2006) Metagenomic analysis of the human distal gut microbiome. Science 312: 1355-1359. doi: 10.1126/science.1124234
    [4] Sartor RB, Mazmanian SK (2012) Intestinal microbes in inflammatory bowel diseases. Am J Gastroenterol Suppl 1: 15-21. doi: 10.1038/ajgsup.2012.4
    [5] Sobhani I, Tap J, Roudot-Thoraval F, et al. (2011) Microbial dysbiosis in colorectal cancer (CRC) patients. PloS one 6: e16393. doi: 10.1371/journal.pone.0016393
    [6] MacFarlane GT, Macfarlane L (2009) Acquisition, evolution and maintenance of the normal gut microbiota. Dig Dis 27: 90-98. doi: 10.1159/000268127
    [7] Vaarala O (2011) The gut as a regulator of early inflammation in type 1 diabetes. Curr Opin Endocrinol Diabetes Obes 18: 241-247. doi: 10.1097/MED.0b013e3283488218
    [8] Power SE, O'Toole PW, Stanton C, et al. (2014) Intestinal microbiota, diet and health. Br J Nutr 111: 387-402. doi: 10.1017/S0007114513002560
    [9] MM OD, Harris HM, Jeffery IB, et al. (2013) The core faecal bacterial microbiome of Irish Thoroughbred racehorses. Lett Appl Microbiol 57: 492-501. doi: 10.1111/lam.12137
    [10] Benson AK, Kelly SA, Legge R, et al. (2010) Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc Natl Acad Sci 107: 18933-18938. doi: 10.1073/pnas.1007028107
    [11] Sommer F, Backhed F (2013) The gut microbiota--masters of host development and physiology. Nat Rev Microbiol 11: 227-238. doi: 10.1038/nrmicro2974
    [12] Carey HV, Walters WA, Knight R (2013) Seasonal restructuring of the ground squirrel gut microbiota over the annual hibernation cycle. Am J Physiol Regul Integr Comp Physiol 304: R33-R42. doi: 10.1152/ajpregu.00387.2012
    [13] MacFariane RD, McLaughlin JJ, Bullock G (1986) Quantitative and qualitative studies of gut flora in striped bass from estuarine and coastal marine environments. J wildl dis 22: 344-348. doi: 10.7589/0090-3558-22.3.344
    [14] Sullam KE, Essinger SD, Lozupone CA, et al. (2012) Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol Ecol 21: 3363-3378. doi: 10.1111/j.1365-294X.2012.05552.x
    [15] Mao S, Zhang M, Liu J, et al. (2015) Characterising the bacterial microbiota across the gastrointestinal tracts of dairy cattle: membership and potential function. Sci Rep 5: 16116. doi: 10.1038/srep16116
    [16] Nelson TM, Rogers TL, Brown MV (2013) The gut bacterial community of mammals from marine and terrestrial habitats. PLoS One 8: e83655. doi: 10.1371/journal.pone.0083655
    [17] Li T, Long M, Gatesoupe F-J, et al. (2015) Comparative Analysis of the Intestinal Bacterial Communities in Different Species of Carp by Pyrosequencing. Microb Ecol 69: 25-36. doi: 10.1007/s00248-014-0480-8
    [18] Wang W, Cao J, Li J-R, et al. (2016) Comparative analysis of the gastrointestinal microbial communities of bar-headed goose (Anser indicus) in different breeding patterns by high-throughput sequencing. Microbiol Res 182: 59-67. doi: 10.1016/j.micres.2015.10.003
    [19] Wilkinson N, Hughes RJ, Aspden WJ, et al. (2016) The gastrointestinal tract microbiota of the Japanese quail, Coturnix japonicaAppl Microbiol Biotechnol 100: 4201-4209. doi: 10.1007/s00253-015-7280-z
    [20] Yang Y, Deng Y, Cao L (2016) Characterising the interspecific variations and convergence of gut microbiota in Anseriformes herbivores at wintering areas. Sci Rep 6: 32655. doi: 10.1038/srep32655
    [21] Li K, Guan W, Wei G, et al. (2007) Phylogenetic analysis of intestinal bacteria in the Chinese mitten crab (Eriocheir sinensis). J Appl Microbiol 103: 675-682. doi: 10.1111/j.1365-2672.2007.03295.x
    [22] Li L, Yan B, Li S, et al. (2016) A comparison of bacterial community structure in seawater pond with shrimp, crab, and shellfish cultures and in non-cultured pond in Ganyu, Eastern China. Ann Microbiol 66: 317-328. doi: 10.1007/s13213-015-1111-4
    [23] Souza DT, Genuário DB, Silva FSP, et al. (2016) Analysis of bacterial composition in marine sponges reveals the influence of host phylogeny and environment. FEMS Microbiol Ecol 93: fiw204-fiw204. doi: 10.1093/femsec/fiw204
    [24] Price JT, Paladino FV, Lamont MM, et al. (2017) Characterization of the juvenile green turtle (Chelonia mydas) microbiome throughout an ontogenetic shift from pelagic to neritic habitats. PloS one 12: e0177642. doi: 10.1371/journal.pone.0177642
    [25] Ahasan MS, Waltzek TB, Huerlimann R, et al. (2017) Fecal bacterial communities of wild-captured and stranded green turtles (Chelonia mydas) on the Great Barrier Reef. FEMS Microbiol Ecol 93: fix139. doi: 10.1093/femsec/fix139
    [26] McDermid KJ, Kittle III RP, Veillet A, et al. (2020) Identification of Gastrointestinal Microbiota in Hawaiian Green Turtles (Chelonia mydas). Evol Bioinf Online 16: 1176934320914603. doi: 10.1177/1176934320914603
    [27] Bjorndal KA, Lutz P, Musick J (1997) Foraging ecology and nutrition of sea turtles. Biol sea turtles 1: 199-231.
    [28] Bjorndal LA (1979) Cellulose digestion and volatile fatty-acid production in the green turtle Chelonia mydasComp Biochem Physiol 63: 127-134. doi: 10.1016/0300-9629(79)90638-8
    [29] Bjorndal KA, Suganuma H, Bolten AB (1991) Digestive fermentation in green turtles, Chelonia mydas, feeding on algae. B Mar Sci 48: 66-71.
    [30] Costa MC, Arroyo LG, Allen-Vercoe E, et al. (2012) Comparison of the fecal microbiota of healthy horses and horses with colitis by high throughput sequencing of the V3-V5 region of the 16S rRNA gene. PLoS One 7: e41484. doi: 10.1371/journal.pone.0041484
    [31] Guarner F, Malagelada JR (2003) Gut flora in health and disease. Lancet 361: 512-519. doi: 10.1016/S0140-6736(03)12489-0
    [32] Ahasan MS, Waltzek TB, Huerlimann R, et al. (2018) Comparative analysis of gut bacterial communities of green turtles (Chelonia mydas) pre-hospitalization and post-rehabilitation by high-throughput sequencing of bacterial 16S rRNA gene. Microbiol Res 207: 91-99. doi: 10.1016/j.micres.2017.11.010
    [33] Bloodgood JC, Hernandez SM, Isaiah A, et al. (2020) The effect of diet on the gastrointestinal microbiome of juvenile rehabilitating green turtles (Chelonia mydas). PLOS One 15: e0227060. doi: 10.1371/journal.pone.0227060
    [34] Morais RA, Dos Santos RG, Longo GO, et al. (2014) Direct evidence for gradual ontogenetic dietary shift in the green turtle, Chelonia mydas. Chelonian Conserv Biol 13: 260-266. doi: 10.2744/CCB-1058.1
    [35] McNally KL (2020) Characterization of Microbial Communities Across Disease States and Environmental Conditions in Kemp's Ridley (Lepidochelys Kempii) and Green Sea Turtles (Chelonia Mydas). Graduate Doctoral Dissertations 551.
    [36] Li D, Chen H, Mao B, et al. (2017) Microbial biogeography and core microbiota of the Rat digestive tract. Sci Rep 8.
    [37] Neumann LM, Dehority BA (2008) An investigation of the relationship between fecal and rumen bacterial concentrations in sheep. Zoo Biol 27: 100-108. doi: 10.1002/zoo.20166
    [38] Eckburg P, Bik E, Bernstein C, et al. (2005) Diversity of the human intestinal microbial flora. Science 308: 1635-1638. doi: 10.1126/science.1110591
    [39] Kelly J, Daly K, Moran AW, et al. (2016) Composition and diversity of mucosa–associated microbiota along the entire length of the pig gastrointestinal tract; dietary influences. Environ Microbiol 19: 1425-1438. doi: 10.1111/1462-2920.13619
    [40] Daly K, Stewart CS, Flint HJ, et al. (2001) Bacterial diversity within the equine large intestine as revealed by molecular analysis of cloned 16S rRNA genes. FEMS Microbiol Ecol 38: 141-151. doi: 10.1111/j.1574-6941.2001.tb00892.x
    [41] (2013) DEHPDepartment of Environment and Heritage Protection. Monitoring and Sampling Manual 2009, Version 2, July 2013, format edits ed.. Australia: Department of Environment and Heritage Protection, QLD.
    [42] Caporaso JG, Kuczynski J, Stombaugh J, et al. (2010) QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7: 335-336. doi: 10.1038/nmeth.f.303
    [43] Masella AP, Bartram AK, Truszkowski JM, et al. (2012) PANDAseq: paired-end assembler for illumina sequences. BMC Bioinformatics 13: 31. doi: 10.1186/1471-2105-13-31
    [44] Edgar RC, Haas BJ, Clemente JC, et al. (2011) UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27: 2194-2200. doi: 10.1093/bioinformatics/btr381
    [45] McMurdie PJ, Holmes S (2013) phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLOS One 8: e61217. doi: 10.1371/journal.pone.0061217
    [46] Lozupone C, Knight R (2005) UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 8228-8235. doi: 10.1128/AEM.71.12.8228-8235.2005
    [47] Zakrzewski M, Proietti C, Ellis JJ, et al. (2017) Calypso: a user-friendly web-server for mining and visualizing microbiome–environment interactions. Bioinformatics 33: 782-783.
    [48] Bray JR, Curtis JT (1957) An ordination of the upland forest communities of Southern Wisconsin. Ecol Monogr 27: 325-349. doi: 10.2307/1942268
    [49] Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B 289-300.
    [50] Bjorndal KA (1995)  Biology and conservation of sea turtles Washington: Smithsonian Institution Press.
    [51] Lutz PL, Musick JA, Wyneken J (2002)  The biology of sea turtles CRC press. doi: 10.1201/9781420040807
    [52] Pedrós-Alió C (2006) Marine microbial diversity: can it be determined? Trends Microbiol 14: 257-263. doi: 10.1016/j.tim.2006.04.007
    [53] Zinger L, Amaral-Zettler LA, Fuhrman JA, et al. (2011) Global patterns of bacterial beta-diversity in seafloor and seawater ecosystems. PLoS One 6: e24570. doi: 10.1371/journal.pone.0024570
    [54] Gu S, Chen D, Zhang JN, et al. (2013) Bacterial community mapping of the mouse gastrointestinal tract. PloS One 8: e74957. doi: 10.1371/journal.pone.0074957
    [55] Flint HJ, Scott KP, Louis P, et al. (2012) The role of the gut microbiota in nutrition and health. Nat Rev Gastroenterol Hepatol 9: 577-589. doi: 10.1038/nrgastro.2012.156
    [56] Byndloss MX, Olsan EE, Rivera-Chávez F, et al. (2017) Microbiota-activated PPAR-γ signaling inhibits dysbiotic Enterobacteriaceae expansion. Science 357: 570-575. doi: 10.1126/science.aam9949
    [57] Campos P, Guivernau M, Prenafeta-Boldú FX, et al. (2018) Fast acquisition of a polysaccharide fermenting gut microbiome by juvenile green turtles Chelonia mydas after settlement in coastal habitats. Microbiome 6: 69. doi: 10.1186/s40168-018-0454-z
    [58] Turnbaugh PJ, Bäckhed F, Fulton L, et al. (2008) Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome. Cell Host Microbe 3: 213-223. doi: 10.1016/j.chom.2008.02.015
    [59] Hong PY, Wheeler E, Cann IK, et al. (2011) Phylogenetic analysis of the fecal microbial community in herbivorous land and marine iguanas of the Galapagos Islands using 16S rRNA-based pyrosequencing. ISME J 5: 1461-1470. doi: 10.1038/ismej.2011.33
    [60] Arthur KE, Boyle MC, Limpus CJ (2008) Ontogenetic changes in diet and habitat use in green sea turtle (Chelonia mydas) life history. Marine Ecology Progress Series 362: 303-311. doi: 10.3354/meps07440
    [61] Kim M, Morrison M, Yu Z (2011) Status of the phylogenetic diversity census of ruminal microbiomes. FEMS Microbiol Ecol 76: 49-63. doi: 10.1111/j.1574-6941.2010.01029.x
    [62] Kim M, Kim J, Kuehn L, et al. (2014) Investigation of bacterial diversity in the feces of cattle fed different diets. J Anim Sci 92: 683-694. doi: 10.2527/jas.2013-6841
    [63] Mao S, Zhang R, Wang D, et al. (2012) The diversity of the fecal bacterial community and its relationship with the concentration of volatile fatty acids in the feces during subacute rumen acidosis in dairy cows. BMC Vet Res 8: 237. doi: 10.1186/1746-6148-8-237
    [64] Kashyap PC, Marcobal A, Ursell LK, et al. (2013) Complex interactions among diet, gastrointestinal transit, and gut microbiota in humanized mice. Gastroenterology 144: 967-977. doi: 10.1053/j.gastro.2013.01.047
    [65] Walker AW, Duncan SH, William Leitch EC, et al. (2005) pH and peptide supply can radically alter bacterial populations and short-chain fatty acid ratios within microbial communities from the human colon. Appl Environ Microbiol 71: 3692-3700. doi: 10.1128/AEM.71.7.3692-3700.2005
    [66] Bartram AK, Lynch MDJ, Stearns JC, et al. (2011) Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end Illumina reads. Appl Environ Microbiol 77: 3846-3852. doi: 10.1128/AEM.02772-10
    [67] Savage DC (1977) Microbial ecology of the gastrointestinal tract. Annu Rev Microbiol 31: 107-133. doi: 10.1146/annurev.mi.31.100177.000543
    [68] Bergey DH, Garrity GM (2004)  Bergey's manual of systematic bacteriology: the proteobacteria Springer.
    [69] Chow WL, Lee YK (2006) Mucosal interactions and gastrointestinal microbiota. Gastrointestinal Microbiology New York, NY: Taylor and Francis, 123-136. doi: 10.3109/9781420014952-7
    [70] Wilson M (2005)  Microbial inhabitants of humans: their ecology and role in health and disease Cambridge University Press.
    [71] Amaro C, Biosca EG (1996) Vibrio vulnificus biotype 2, pathogenic for eels, is also an opportunistic pathogen for humans. Appl Environ Microbiol 62: 1454-1457. doi: 10.1128/AEM.62.4.1454-1457.1996
    [72] Brenner DJ, Farmer J (1984) Enterobacteriaceae. Bergey's Manual of Systematics of Archaea and Bacteria .
    [73] Dworkin M, Falkow S, Rosenberg E, et al. (2006) Proteobacteria: Gamma Subclass. The Prokaryotes Springer Science & Business Media.
    [74] Schaal K, Schofield G, Pulverer G (1980) Taxonomy and clinical significance of Actinomycetaceae and Propionibacteriaceae. Infection 8: S122-S130. doi: 10.1007/BF01639868
    [75] Stackebrandt E, Cummins CS, Johnson JL (2006) Family propionibacteriaceae: the genus Propionibacterium. The prokaryotes Springer, 400-418. doi: 10.1007/0-387-30743-5_19
    [76] Wang K, Lu W, Tu Q, et al. (2016) Preliminary analysis of salivary microbiome and their potential roles in oral lichen planus. Sci Rep 6: 22943. doi: 10.1038/srep22943
    [77] Flint HJ, Scott KP, Duncan SH, et al. (2012) Microbial degradation of complex carbohydrates in the gut. Gut Microbes 3: 289-306. doi: 10.4161/gmic.19897
    [78] Martens EC, Lowe EC, Chiang H, et al. (2011) Recognition and degradation of plant cell wall polysaccharides by two human gut symbionts. PLoS Biol 9: e1001221. doi: 10.1371/journal.pbio.1001221
    [79] Salyers A, Vercellotti J, West S, et al. (1977) Fermentation of mucin and plant polysaccharides by strains of Bacteroides from the human colon. Appl Environ Microbiol 33: 319-322. doi: 10.1128/AEM.33.2.319-322.1977
  • microbiol-06-04-022-s001.pdf
  • This article has been cited by:

    1. KL McNally, CR Mott, JR Guertin, JL Bowen, Microbial communities of wild-captured Kemp’s ridley (Lepidochelys kempii) and green sea turtles (Chelonia mydas), 2021, 45, 1863-5407, 21, 10.3354/esr01116
    2. Lucía Díaz-Abad, Natassia Bacco-Mannina, Fernando Miguel Madeira, Ester A. Serrao, Aissa Regalla, Ana R. Patrício, Pedro R. Frade, Red, Gold and Green: Microbial Contribution of Rhodophyta and Other Algae to Green Turtle (Chelonia mydas) Gut Microbiome, 2022, 10, 2076-2607, 1988, 10.3390/microorganisms10101988
    3. Karina Jones, Ellen Ariel, 2023, 9781119839675, 136, 10.1002/9781119839729.ch4
    4. Klara Filek, Adriana Trotta, Romana Gračan, Antonio Di Bello, Marialaura Corrente, Sunčica Bosak, Characterization of oral and cloacal microbial communities of wild and rehabilitated loggerhead sea turtles (Caretta caretta), 2021, 3, 2524-4671, 10.1186/s42523-021-00120-5
    5. ZR Forbes, AK Scro, SH Patel, KM Dourdeville, RL Prescott, RM Smolowitz, Fecal and cloacal microbiomes of cold-stunned loggerhead Caretta caretta, Kemp’s ridley Lepidochelys kempii, and green sea turtles Chelonia mydas, 2023, 50, 1863-5407, 93, 10.3354/esr01220
    6. Frances Widjaja, Ivonne M. C. M. Rietjens, From-Toilet-to-Freezer: A Review on Requirements for an Automatic Protocol to Collect and Store Human Fecal Samples for Research Purposes, 2023, 11, 2227-9059, 2658, 10.3390/biomedicines11102658
    7. T. Franciscus Scheelings, Thi Thu Hao Van, Robert J. Moore, Lee F. Skerratt, Location Matters: Variations in Cloacal Microbiota Composition of Spatially Separated Freshwater Turtles, 2024, 87, 0095-3628, 10.1007/s00248-024-02452-6
    8. Deirdre Mikkelsen, Alexandra M McGowan, Justine S Gibson, Janet M Lanyon, Sara Horsman, Jennifer M Seddon, Faecal bacterial communities differ amongst discrete foraging populations of dugongs along the east Australian coast, 2024, 100, 1574-6941, 10.1093/femsec/fiae051
    9. Xin Li Ching, Syamsyahidah Samsol, Mohd Uzair Rusli, Mohamad Aqmal-Naser, Joseph Anak Bidai, Christian Sonne, Xin Wu, Nyuk Ling Ma, Blood and cloacal microbiome profile of captive green turtles (Chelonia mydas) and hawksbill turtles (Eretmochelys imbricata): Water quality and conservation implications, 2025, 375, 00456535, 144223, 10.1016/j.chemosphere.2025.144223
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5240) PDF downloads(118) Cited by(9)

Figures and Tables

Figures(6)  /  Tables(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog