Original Research ARTICLE
Comparative Analyses of the Digestive Tract Microbiota of New Guinean Passerine Birds
- 1Section for Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark
- 2Section for Biosystematics, Natural History Museum of Denmark, Copenhagen, Denmark
- 3Biology Centre AS CR v. v. i., Faculty of Science, Institute of Entomology and University of South Bohemia, Ceske Budejovice, Czechia
The digestive tract microbiota (DTM) plays a plethora of functions that enable hosts to exploit novel niches. However, our understanding of the DTM of birds, particularly passerines, and the turnover of microbial communities along the digestive tract are limited. To better understand how passerine DTMs are assembled, and how the composition changes along the digestive tract, we investigated the DTM of seven different compartments along the digestive tract of nine New Guinean passerine bird species using Illumina MiSeq sequencing of the V4 region of the 16S rRNA. Overall, passerine DTMs were dominated by the phyla Firmicutes and Proteobacteria. We found bird species-specific DTM assemblages and the DTM of different compartments from the same species tended to cluster together. We also found a notable relationship between gut community similarity and feeding guilds (insectivores vs. omnivores). The dominant bacterial genera tended to differ between insectivores and omnivores, with insectivores mainly having lactic acid bacteria that may contribute to the breakdown of carbohydrates. Omnivorous DTMs were more diverse than insectivores and dominated by the bacterial phyla Proteobacteria and Tenericutes. These bacteria may contribute to nitrogen metabolism, and the diverse omnivorous DTMs may allow for more flexibility with varying food availability as these species have wider feeding niches. In well-sampled omnivorous species, the dominant bacterial genera changed along the digestive tracts, which was less prominent for insectivores. In conclusion, the DTMs of New Guinean passerines seem to be species specific and, at least in part, be shaped by bird diet. The sampling of DTM along the digestive tract improved capturing of a more complete set of members, with implications for our understanding of the interactions between symbiont and gut compartment functions.
Animal bodies are sophisticated centers of symbiotic interactions (Ley et al., 2008; Huttenhower and The human microbiome project consortium, 2012; Macke et al., 2017) that have allowed for adaptations enabling hosts to exploit new niches (Russell et al., 2009; Godoy-Vitorino et al., 2012; Dietrich et al., 2014; Macke et al., 2017). Of the many microbial communities associated with hosts, the digestive tract microbiota (DTM) is likely the most diverse, fulfilling essential functions in nutrition absorption, breakdown of macromolecules, vitamin synthesis, detoxification, and defense (Kohl, 2012; Poulsen et al., 2014; Ceja-Navarro et al., 2015). The DTM is not only of importance in the digestive tract, but also impacts the growth, development, and behavior of hosts (Round and Mazmanian, 2009; Sommer and Bäckhed, 2013). Although extensive work on DTMs has been carried out on a number of hosts, a thorough understanding of the composition and role of the DTM is limited to relatively few taxa, including mammals (Ley et al., 2008; Huttenhower and The human microbiome project consortium, 2012; Kohl et al., 2018), insects (Macke et al., 2017), and some captive and livestock animals (Waite and Taylor, 2015; Hird, 2017).
Despite the diversity and global distribution of birds, studies of their microbial communities are limited (Lee, 2015; Waite and Taylor, 2015), especially in wild populations (Hird, 2017). The majority of avian DTM studies have been performed on commercial animals and charismatic species such as penguins (order: Sphenisciformes, family: Spheniscidae), hoatzin (order: Opisthocomiformes, family: Opisthocomidae), kakapo (order: Psittaciformes, family: Strigopidae), and vultures (Order: Accipitriformes, family: Cathartidae and Accipitridae) (Godoy-Vitorino et al., 2012; Kohl, 2012; Dewar et al., 2013; Roggenbuck et al., 2014; Vela et al., 2015; Waite and Taylor, 2015). Although Passeriformes comprises more than 50% of bird diversity (~6,000 species) (Dickinson and Christidis, 2014), their DTM composition has only been explored in fewer than 100 species (e.g., Hird et al., 2015; Kropáčková et al., 2017; García-Amado et al., 2018). Passerine birds occupy a multitude of habitats and feeding guilds, making them an ideal group for investigation of the evolution and roles of DTM associations. This is also likely to provide a better understanding of how the DTM influences passerine richness, distribution, and environmental adaptations.
The assemblage of DTM may vary according to the function of each region in the digestive tract (Kohl et al., 2018). For example, the acidic stomach environment is one of the major regions of digestion of macromolecules, and would thus be expected to house a community of acidophilus bacteria that play roles in digestion. In contrast, the lower regions of the small intestine and large intestine likely harbor specialized bacteria that contribute to nutrient absorption and processing (Stevens and Hume, 1995; Lu and Domingo, 2008; Svihus, 2014). However, most avian DTM studies have been based only on cloacal swabs or fecal samples (Waite and Taylor, 2015), and these alone are unlikely to fully represent the microbial communities across the entire digestive tract. This assumption is supported by a few studies that have compared different digestive tract compartments and found differences in microbial community compositions (Videvall et al., 2017; Zhang et al., 2017; Drovetski et al., 2018; García-Amado et al., 2018). These studies demonstrated that the DTMs of cloacal or fecal samples may qualitatively represent the microbiota of other sections of the digestive tract but not quantitatively (Videvall et al., 2017; Zhang et al., 2017). Lu and Domingo (2008) also showed that the microbial genes in feces are skewed toward genes related to nutrient absorption and processing rather than to digestion. This indicates that the fecal microbiota may be overrepresented by bacteria in the last regions of the digestive tract, where nutrient metabolism takes place (Stevens and Hume, 1995). More comparative studies on the DTM of digestive tract compartments of wild birds may improve our understanding of what determines the composition of microbial communities in natural conditions and identify the bacterial symbionts that play important roles along the digestive tract. In this study we investigate multiple digestive tract compartments of nine New Guinean passerine bird species belonging to two feeding guilds (insectivores and omnivores); Colluricincla megarhyncha, Ifrita kowaldi, Rhipidura atra, Crateroscelis robusta, Sericornis nouhuysi (insectivores), Melipotes fumigatus, Melanocharis nigra, Melanocharis versteri, Toxorhamphus poliopterus (omnivores), and assess the appropriateness of investigating single compartments of the digestive tract to capture the full DTM diversity.
Materials and Methods
A total of 33 alcohol bird specimens from the museum collections at the Natural History Museum of Denmark were used in this study. Twenty-nine specimens were collected in 2015 and 2016 in Madang Province in Papua New Guinea (PNG) at five elevations, ranging from 200 to 3,700 m a.s.l., along the Mount Wilhelm elevational transect. Three specimens were collected at Wasaunon (~2,950 m a.s.l.) in the Saruwaged Range on the Huon Peninsula in Morobe Province and one specimen was collected in the lowland rainforest of the Wanang Conservation Area (200 m a.s.l.) in Madang Province (Figure 1A, Table S1). Birds were captured using mist-nets, euthanized and immediately injected with 96% ethanol and preserved in 99% ethanol. An additional 51 regurgitated samples that had been collected between 2010 and 2012 from the same locations were included to analyse the crop microbiota. These samples were collected using the tartar emetic method to investigate the variation in diet of birds along the Mount Wilhelm elevational gradient, as described in detail by Sam et al. (2017). These samples were stored in 96% ethanol and select species were analyzed (Table 1 and Table S1).
Figure 1. (A) A map of the sampling localities in Papua New Guinea: the Wanang conservation area, the Mount Wilhelm transect in Madang province, and Wasaunon on the Huon Peninsula in the Morobe province. (B) The digestive tract of Melipotes fumigatus indicating six compartments dissected from the alcohol specimen. The regurgitated samples were used to analyse the crop (not shown) microbiota.
Table 1. The nine Papua New Guinean passerine (Order: Passeriformes) bird species and number of original samples per sample type (R, regurgitated samples; A, dissected alcohol specimens), feeding guilds (Sam et al., 2017), and habitat and elevation (Pratt and Beehler, 2015).
Digestive tracts were dissected from alcohol specimens of nine species (Table 1) to represent multiple passerine clades (Corvides, Passerides, and Meliphagides). One to six individuals of each species were included (Table 2) and the number of bird individuals per species was governed by the availability of alcohol specimens. Six segments (stomach, beginning of the small intestine, middle of the small intestine, ceca region and the ileum, large intestine, and cloacal region) (Figure 1B) of the digestive tract were used.
Table 2. The number of samples used for DNA extraction (S), samples that amplified during MiSeq sequencing (A) and samples that passed quality checks, filtering, and were deemed suitable based on visual inspection of the rarefaction curves (Figure S1) (Q).
DNA Extraction, Initial PCR, and MiSeq Amplicon Sequencing
DNA was extracted using DNeasy Blood and Tissue kits (Qiagen, Germany) following the protocol, except that only 75 μl AE buffer was used to elute the DNA. The amount of tissue varied from 15 to 30 mg depending on the compartment, and the extraction protocol was adjusted accordingly. Diagnostic PCR was conducted at the University of Copenhagen to confirm the presence of sufficient bacterial DNA. We used primers for the V4 region of the 16S rRNA gene (′V4.SA504 and ′V4.SB711) following Otani et al. (2016) and PCRs were conducted under conditions of 96°C for 4 min, 35 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 30 s, and an extension step of 74°C for 4 min (cf., Otani et al., 2016). Amplifications were assessed on a 1.5% agarose gel.
Successfully amplified samples from the diagnostic PCR were sequenced by the Microbial System Molecular Biology lab at the University of Michigan using Illumina MiSeq. PCR reactions were 20 μl and the master mix contained AccuPrime PCR buffer 2, AccuPrime HiFi polymerase and 4 μM of the same primers (′V4.SA504 and ′V4.SB711). Reactions were performed under the PCR conditions of denaturing of 96°C for 1 min, 30 cycles of annealing (95°C for 20 s, 55°C 15 s, and 72°C for 5 min), and extension at 72°C for 10 min.
MiSeq sequences were analyzed using Mothur 1.35.1 (Kozich et al., 2013). Chimeric sequences were filtered using the VSEARCH package and non-bacterial sequences (e.g., mitochondrial sequences) were removed. Operational Taxonomic Units (OTUs) (Cut-off 97%) of bacteria were aligned and identified using the SILVA bacterial reference library (Kozich et al., 2013). Sample coverage was evaluated by generating rarefaction curves with Mothur 1.35.1 (Figure S1). Samples with less than 2,500 reads and samples that did not display a flattening rarefaction curve based on visual inspection were excluded (Table 2). We also calculated the inverse Simpson's diversity index and the Bray-Curtis dissimilarity of bacterial communities in different compartments. We used vegan and Bidiversity R packages in R version 3.5.0 (R Core Team, 2013) to conduct statistical analysis. To visualize community similarities, we conducted a principal component analysis (PCA), and a scree plot analysis was conducted to evaluate the variation explained by the first ten principal components. We also conducted an adonis statistical test (a permutational multivariate analysis of variance using Bray-Curtis distance matrices) with 999 permutations to examine whether there were species-level and digestive tract compartment-level differences in DTMs. FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/) was used to generate trees.
After a chimera check and quality filtering, a total of 7,505,272 bacterial DNA sequences were retained from 149 intestinal compartment samples. The number of sequences per sample ranged from 2,904 to 261,625 (average 50,371), classified to 15,767 OTUs at the 97% sequence-similarity cut-off level (range 29–4,017 per sample; Table S2). Classified sequences belonged to 27 phyla, dominated by Firmicutes (49.2% of the OTUs), followed by Proteobacteria (25.6%) and Actinobacteria (12.7%). Only 5% of the sequences could not be classified at the phylum level (Table S3). The most abundant genera include an unclassified Carnobacteriaceae genus (Firmicutes) representing 31.4% of all sequences, followed by an unclassified Enterobacteriaceae genus (8%: Proteobacteria), Lactobacillus (4.9%: Firmicutes), an unclassified Gammaproteobacteria genus (4.7%: Proteobacteria), an unclassified Actinomycetales genus (2.6%: Actinobacteria), Rothia (2.5%: Actinobacteria), an unclassified Rhizobiales genus (1.9%: Proteobacteria), and an unclassified Mycoplasmataceae genus (1.6%: Tenericutes). We did not distinguish between intra- and extra-cellular bacteria in our approaches, but since we only found one well-represented genus of a (facultative) intracellular bacteria (Bartonella) in one individual of R. atra, we anticipate that the bacteria identified are most likely to primarily be extracellular.
The adonis tests showed that DTMs were significantly different between bird species (p = 0.001, R2 = 0.24) but there were no significant difference between digestive tract compartments (p = 0.627, R2 = 0.04). Loading value analyses (Table S4) indicated that OTU 01 (Firmicutes, unclassified Carnobacteriaceae) contributed the most to the separation of S. nouhuysi, C. robusta, and R. atra communities (insectivores) in PCA space, while OTU 04 (Firmicutes, Lactobacillus) contributed to the separation of most M. nigra samples from the remaining communities. Furthermore, OTU 08 (Tenericutes, unclassified Mycoplasmataceae) and OTU 10 (Proteobacteria, Helicobacter) accounted for the separation of M. versteri and T. poliopterus (omnivores) from the others (Figure 2). Overall, Bray-Curtis distances revealed that bacterial communities in insectivorous species were more similar to each other than bacterial communities were to each other in omnivorous species (Figure 3). However, clustering was not in two perfectly distinct groups; rather, there were multiple clusters comprised of either insectivorous or omnivorous DTMs, and the omnivores were sub-clustered according to nectar feeders vs. fruit feeders (Figure 3). Different gut compartments of insectivorous species were more similar than those of omnivorous species, and compartments from the same individual were more similar to each other than to those of other individuals (Figures 2, 3).
Figure 2. Principal component analysis (PCA) of DTM communities in different digestive tract compartments (shapes) in nine New Guinean passerine species (colors). The main dietary guilds of the bird species are indicated in front of their names in the figure key. The scree plot represents the variation in the data captured by the first 10 principal components. Arrows indicate the OTUs that mainly drive the microbial community separation: OTU 01 (unclassified Carnobacteriaceae genus, Firmicutes) and OTU 04 (Lactobacillus, Firmicutes), OTU 08 (unclassified Mycoplasmataceae genus, Tenericutes), and OTU 10 (Helicobacter, Proteobacteria).
Figure 3. Dendrogram based on Bray-Curtis dissimilarity of DTMs between different compartments (different shapes) of the digestive tract of nine New Guinean passerine species (different colors). Gray scale shading of clades represents the feeding guilds of the bird species.
The inverse Simpson's diversity index of the DTM between different compartments varied between bird species; however, variation, based on the standard errors, between individuals was low in well-sampled species (Figure 4). Generally, diversity increased as one moved along the digestive tract in most insectivorous species, which had lower diversity than omnivorous species, with the exception being I. kowaldi. This trend was less consistent for omnivorous species, but lower diversity was observed in the midgut of most species, except for M. versteri. Statistical tests for individual species was not possible due to the low sample sizes.
Figure 4. Average ± SE Inverse Simpson's diversity index of DTM in different digestive tract compartments. The infraorders of the species and the dietary guilds indicated above the figure. Numbers of replicates are given in parentheses.
Insectivores and omnivores displayed distinct differences in abundances of bacterial phyla. The DTM of insectivorous species were dominated by Firmicutes, except for I. kowaldi, which had a higher representation of Actinobacteria, whereas Proteobacteria dominated the DTM of omnivores. The nectar-feeding omnivore T. poliopterus had a higher number of OTUs belong to the Tenericutes. Furthermore, omnivores had a higher diversity of bacterial phyla compared to insectivores (Figure 5).
Figure 5. The proportion of major bacterial phyla present in the DTMs of different digestive tract compartments. The dendrogram represents the phylogenetic relationship between the bird species and the tips of the dendrogram represent the feeding guild of each species: black: insectivores, solid gray: omnivores (fruits and insects), and dashed gray: omnivores (nectar and insects). Red line indicates the separation of the bird infraorders. Photo credits: Andrew Hart Reeve.
Due to the different abundances of bacterial phyla, we determined the 20 most abundant genera separately for each feeding guild, where only nine of the most abundant genera were shared (Table 3). Nine of the 20 most abundant genera in insectivores belonged to the Firmicutes (Figure 6). On average 75.2 ± 8.6% (SE) [83.4 ± 3.3% (SE) if we exclude I. kowaldi] of the OTUs belonged to the 20 most abundant genera, with the unclassified Carnobacteriaceae genus typically dominating. I. kowaldi appears to be a special case within the insectivores, with a more diverse and more evenly distributed microbiota. Consequently, only 42.6% of the sequences were covered by the 20 most abundant insectivorous genera. In omnivorous birds, eight out of the 20 most abundant genera were Proteobacteria and they represented 59 ± 9.7% (SE) of the OTUs, with the dominant bacterial genera varying between species. M. nigra was dominated by the genus Lactobacillus, M. versteri by Helicobacter and multiple other proteobacterial genera, the T. poliopterus DTM was mainly dominated by Ureaplasma and an unclassified genus in the Mycoplasmataceae, and the twenty most abundant genera were more evenly distributed in M. fumigatus (Figure 7).
Figure 6. Relative abundance of the twenty most abundant bacterial genera in digestive tract compartments of insectivorous bird species. Bacterial genera belonging to the same phylum are represented by similar shades of colors.
Figure 7. Relative abundance of the twenty most abundant bacterial genera in different digestive tract compartments of omnivorous bird species. Bacterial genera belonging to the same phylum are represented by similar shades of colors.
In insectivorous species, the stomach and the beginning of the small intestine were typically dominated by 2–3 bacterial genera mainly belonging to Firmicutes (Figures 5, 6), while bacterial communities in the crop and the lower digestive tract were inhabited by more genera and phyla, such as Firmicutes, Proteobacteria, and Verrucomicrobia. In contrast, omnivores showed higher species-level and individual variation in the distribution of their 20 most common bacterial genera along the digestive tract (Figure 7). Moreover, the DTM of different compartments were typically dominated by multiple genera of Proteobacteria; although, the dominant genera tended to change along the digestive tract.
We investigated the DTM along the digestive tract of nine New Guinean passerine bird species and found that Firmicutes and Proteobacteria dominate the DTMs, which is consistent with previous studies of the passerine gut microbiota (Hird et al., 2014, 2015; Lewis et al., 2016; Kropáčková et al., 2017; García-Amado et al., 2018; Teyssier et al., 2018). DTM assemblages were species-specific and mainly shaped by diet, which agrees to some extent with recent studies on Neotropical (Hird et al., 2015; García-Amado et al., 2018) and European birds (Kropáčková et al., 2017). Even though we found species-specific DTMs, our microbial communities did not aggregate according to the phylogenetic placement of host birds; rather they clustered according to bird feeding guild. Previous studies on passerines have suggested that host phylogeny plays a main role in DTM assembly and ecological and life histories, such as feeding guilds, play secondary roles in shaping the DTM (Hird et al., 2015; Kropáčková et al., 2017; García-Amado et al., 2018). Our findings, however, are consistent with studies on mammals, for which host diet appears to be more important than host phylogeny (Ley et al., 2008; Muegge et al., 2011). Similar to previous DTM studies in passerines our study also demonstrates high individual variation in the DTM suggesting that the inclusion of more individuals is needed to assess the levels of intra-specific DTM variation.
The observed differences between the two feeding guilds (insectivores and omnivores) are most likely the result of differences in digestive needs (cf. Wang et al., 2016). Insectivorous diets are dominated by protein-rich arthropods, whereas omnivorous diets contain fruits, seeds, nectar, and insects (Sam et al., 2017). Lactic acid bacteria such as Lactobacillus, Enterococcus, and Carnobacteriaceae genera dominate the digestive tract compartments of insectivores. One of their main functions in the human digestive tract is carbohydrate metabolism (Hammes and Hertel, 2006; Pikuta, 2014), and a similar function is expected in birds. Many Lactobacillus and Carnobacteriaceae genera also produce antimicrobial substances that may play roles in the defense against antagonists, and they may play a role in detoxification through bile acid hydrolysis and removal of by-products of protein hydrolysis (Hammes and Hertel, 2006). Enterococcus bacteria can hydrolyse amino acids (LeBlanc, 2006), which has been shown to be important in carnivorous mammals to generate energy (Muegge et al., 2011). Thus, it is conceivable that insectivorous passerines with a protein-rich diet benefit in similar ways from lactic acid bacteria to improve energy uptake from protein-rich diets and by reducing harmful by-products of protein hydrolysis.
Omnivorous species, most of which had Proteobacteria-rich DTMs, were more variable in their DTMs between species than insectivorous birds. The majority of the OTUs driving these differences, e.g., Gammaproteobacteria (Enterobacteriaceae) and Alphaproteobacteria (Rhizobiaceae), include bacterial species that are capable of nitrogen fixation and amino acid synthesis (Panina et al., 2001; Zehr et al., 2003). The presence of these bacteria may thus contribute to the nitrogen budget of omnivorous birds, which feed on plant material that is generally low in nitrogen. Accordingly, Muegge et al. (2011) showed higher levels of amino acid synthesizing enzymes coded for by members of the microbiota in herbivorous mammals. A similar specialized DTM member function may apply for the nectar-feeding omnivore T. poliopterus. Nectar contains little nitrogen, so Mycoplasmataceae in the small and large intestines of T. poliopterus may help conserve nitrogen, as bacteria in this family are known to be able to hydrolyse urea (Glass et al., 2000). This is consistent with work by Preest et al. (2003), who suggested that uric acid may be broken down to re-cycle nitrogen in the lower intestine of nectar feeders, such as hummingbirds.
Some DTMs of the omnivorous M. nigra were more similar to those of insectivores, which is consistent with findings by Sam et al. (2017) that insects are a major component of the diet in some individuals of this species. Variable diets between individuals thus likely induce variation in DTMs. This is supported by the finding that the DTMs of some individuals of M. nigra are similar to insectivores, while others resemble omnivorous species. Comparable variability in DTM composition was recently documented in frugivorous Neotropical passerines (García-Amado et al., 2018), and a potentially more flexible and diverse DTM may enable omnivorous species to adjust their diet according to the availability of food resources. It appears a promising avenue of further research to determine how particular diets may alter the DTM composition for a given individual.
The DTM of I. kowaldi, which contains batrachotoxins in its tissue (possibly stemming from a diet of beetles, which also contain this toxin) (Dumbacher et al., 2000, 2008), was markedly different from all other insectivores in this study, with a higher proportion of Actinobacteria, including Conexibacter and Pseudonocardia. Bacteria from the genus Pseudonocardia inhabit toxic environments and typically have capacities to break down toxins (Veinberg et al., 2006), and members of the genus associate with fungus-farming ants for antimicrobial production (Cafaro et al., 2011). The genome of Conexibacter has a high proportion (11%) of genes related to amino acid metabolism (Pukall et al., 2010), suggesting a potential importance of Conexibacter for nitrogen recycling. Alternatively, the distinct composition of the DTM of I. kowaldi might be a result of its isolated phylogenetic position in the passerine phylogeny, within which it is the sole member of the family (Ifritidae) (Dickinson and Christidis, 2014). We were only able to include a single individual and were therefore unable to assess individual variation within this species. Future work comparing more individuals is clearly needed to corroborate our findings and to determine the reasons for the potentially relatively unique DTM composition of I. kowaldi.
We did not find common avian pathogenic bacterial genera, such as Salmonella, in our study. Some species belonging to common genera identified, such as Lactobacillus, Enterococcus, and Staphylococcus, can be pathogenic in birds (Benskin et al., 2009). However, these genera also contain a multitude of beneficial bacterial symbionts (Hammes and Hertel, 2006; LeBlanc, 2006).
Thus far, the majority of studies of avian DTM have focused on single compartments of the digestive tract to assess the composition of the entire digestive tract (Lee, 2015; Waite and Taylor, 2015) with evidence both for (Sekelja et al., 2012) and against (Lu and Domingo, 2008; Stanley et al., 2015) such an approach. We found notable species-specific differences in major bacterial groups along the digestive tract, with changes apparently depending on whether the birds were insectivores or omnivores. The crop region contained a higher DTM diversity and a different profile of bacteria compared to the mid gut regions. Some of this diversity might be accounted for by bacteria that entered with the food. We did not pursue this matter further but speculate that the reduced diversity and the changed composition of DTM in the stomach regions may be caused by the acidic stomach, which may act as a barrier to exogenous bacteria. Compared to the mid gut regions, the DTM in the large intestine and cloacal regions were more diverse and reflected the majority of bacterial OTUs from other regions of the digestive tract. The changes in DTM diversity along the passerine digestive tract and higher diversity in the hind regions, is consistent with a recent study on Deer mice (Kohl et al., 2018). Focusing on hind regions of the digestive tract thus likely provide sufficient qualitative insights into what microbes are present. However, they may not be sufficient to identify quantitative differences in major bacterial groups along the digestive tract, which align with previous comparative studies (Sekelja et al., 2012; Stanley et al., 2015; Videvall et al., 2017; Zhang et al., 2017). Moreover, the cloaca may also be affected by factors that mask the true endogenous DTM, such as the transfer of microbes during mating (White et al., 2010).
Our study was able to use museum alcohol specimens to investigate the DTM of multiple passerine species. This provides the first insights into New Guinean passerine DTMs and how DTMs change along the digestive tract of wild birds. We show that the DTM of insectivores is more consistent among species than omnivores and that the DTMs change along the digestive tract. The omnivorous DTMs were more variable but also show a change in major bacterial groups in different regions of the digestive tract. These differences in the DTM between feeding guilds appear to be driven by diet, where insectivores have a somewhat consistent diet based on insects, while omnivores have a more flexible diet including animal and plant material. Improved sampling of freshly dissected guts for an increased number of passerine species will enable a better understanding of the evolution of DTMs in this diverse bird order. It will also be relevant to explore different gut regions to improve our understanding of roles of these DTMs in different compartments based on the function of the digestive tract region.
MiSeq data is available from the NCBI SRA archive (SRP143623).
The fieldwork protocol, including ethical standards, was approved by the ethics committee at the Natural History Museum of Denmark. The fieldwork protocol follows Guidelines to the Use of Wild Birds in Research formulated by The Ornithological Council.
KB dissected the digestive tracts and conducted the molecular work and sequence analyses. KJ collected bird specimens. KS collected the regurgitated samples. KB, KJ, and MP wrote the manuscript.
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.
This work was funded by the Danida Foundation to KB (A31661), The Villum Kann Rasmussen Foundation to MP (VKR10101) and KJ (VKR15560), GACR to KS (14-36098G), and The Carlsberg Foundation to MP (CF17-0597) and KJ (CF17-0248). We would like to thank Vojtech Novotny and the staff at the Binatang Research Center in Madang, PNG for their help with collection of the samples. We also acknowledge the PNG Immigration & Citizen Service Authority and the Department of Environment and Conservation, Boroko, Papua New Guinea for granting research and export permits. We thank Andrew Hart Reeve for permission to use the photos included in Figure 5. Finally, we thank Jan Bolding Kristensen, Saria Otani, Sylvia Mathiasen, Haofu Hu, Panagiotis Sapountzis, Ana Barrientos, Lisa Brandenborg, and David R. Nash for help and guidance throughout this project.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01830/full#supplementary-material
Benskin, C. M. H., Wilson, K., Jones, K., and Hartley, I. R. (2009). Bacterial pathogens in wild birds: a review of the frequency and effects of infection. Biol. Rev. 84, 349–373. doi: 10.1111/j.1469-185X.2008.00076.x
Cafaro, M. J., Poulsen, M., Little, A. E. F., Price, S. L., Gerardo, N. M., Wong, B., et al. (2011). Specificity in the symbiotic association between fungus-growing ants and protective Pseudonocardia bacteria. Proc. Biol. Sci. 278, 1814–1822. doi: 10.1098/rspb
Ceja-Navarro, J. A., Vega, F. A., Karaoz, U., Hao, Z., Jenkins, S., Lim, H. C., et al. (2015). Gut microbiota mediate caffeine detoxification in the primary insect pest of coffee. Nat. Commun. 6:7618. doi: 10.1038/ncomms8618
Dewar, M. L., Arnold, J. P. Y., Dann, P., Trathan, P., Groscolas, R., and Smith, S. (2013). Interspecific variations in the gastrointestinal microbiota in penguins. Microbiologyopen 2, 195–204. doi: 10.1002/mbo3.66
Dietrich, C., Köhler, T., and Brune, A. (2014). The cockroach origin of the termite gut microbiota: patterns in bacterial community structure reflect major evolutionary events. Appl. Environ. Microbiol. 80, 2261–2269. doi: 10.1128/AEM.04206-13
Drovetski, S. V., O'Mahony, M., Ransome, E. J., Matterson, K. O., Lim, H. C., Chesser, R. T., et al. (2018). Spatial organization of the gastrointestinal microbiota in urban Canada geese. Sci. Rep. 8:3713. doi: 10.1038/s41598-018-21892-y
Dumbacher, J. P., Deiner, K., Thompson, L., and Fleischer, R. C. (2008). Phylogeny of the avian genus Pitohui and the evolution of toxicity in bird. Mol. Phylogenet. Evol. 49, 774–781. doi: 10.1016/j.ympev
Dumbacher, J. P., Spande, T. F., and Daly, J. W. (2000). Batrachotoxin alkaloids from passerine birds: a second toxic bird genus (Ifrita kowaldi) from New Guinea. Proc. Natl. Acad. Sci. U.S.A. 97, 12970–12975. doi: 10.1073/pnas.200346897
García-Amado, M. A., Shin, H., Sanz, V., Lentino, M., Martínez, L. M., Contreras, M., et al. (2018). Comparison of gizzard and intestinal microbiota of wild neotropical birds. PLoS ONE 13:e0194857. doi: 10.1371/journal.pone.0194857
Glass, J. I., Lefkowitz, E. J., Glass, J. S., Heiner, C. R., Chen, E. Y., and Cassell, G. H. (2000). The complete sequence of the mucosal pathogen Ureaplasma urealyticum. Nature 407, 757–762. doi: 10.1038/35037619
Godoy-Vitorino, F., Goldfarb, K. C., Karaoz, U., Leal, S., Garcia-Amado, M. A., Hugenholtz, P., et al. (2012). Comparative analyses of foregut and hindgut bacterial communities in hoatzins and cows. ISME J. 6, 531–541. doi: 10.1038/ismej.2011.131
Hammes, W. P., and Hertel, C. (2006). “The genera of Lactobacillus and Carnobacterium,” in The Prokaryotes: A Handbook on the Biology of Bacteria, eds M. Dworkin, S. Falkow, E. Rosenberg, K. Schleifer, and E. Stackebrandt (New York, NY: Springer-Verlag), 320–403.
Hird, S. M., Carstens, B. C., Cardiff, S. W., Dittmann, D. L., and Brumfield, R. B. (2014). Sampling locality is more detectable than taxonomy or ecology in the gut microbiota of the brood-parasitic Brown-headed Cowbird (Molothrus ater). PeerJ 2:e321. doi: 10.7717/peerj.321
Kohl, K. D., Dearing, M. D., and Bordenstein, S. R. (2018). Microbial communities exhibit host species distinguishability and phylosymbiosis along the length of the gastrointestinal tract. Mol. Ecol. 27, 1874–1883. doi: 10.1111/mec.14460
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Kropáčková, L., Tesicky, M., Albrecht, T., Kubovčiak, J., Cisková, D., Tomásek, O., et al. (2017). Codiversification of gastrointestinal microbiota and phylogeny in passerines is not explained by ecological divergence. Mol. Ecol. 26, 5292–5304. doi: 10.1111/mec.14144
LeBlanc, D. J. (2006). “Enterococcus,” in “The Prokaryotes: A Handbook on the Biology of Bacteria,” eds M. Dworkin, S. Falkow, E. Rosenberg, K. Schleifer, and E. Stackebrandt (New York, NY: Springer Science), 175–204.
Lee, W. Y. (2015). Avian gut microbiota and behavioral studies. Korean J. Ornithol. 22, 1–11. Available online at: https://www.researchgate.net/profile/Won_Young_Lee2/publication/281350575_Avian_gut_microbiota_and_behavioral_studies/links/55e3ac1908aede0b5733cae4/Avian-gut-microbiota-and-behavioral-studies.pdf
Lewis, W. B., Moore, F. R., and Wang, S. (2016). Characterization of the gut microbiota of migratory passerines during stopover along the northern coast of the Gulf of Mexico. J. Avian Biol. 47, 659–668. doi: 10.1111/jav.00954
Lu, J., and Domingo, J. S. (2008). Turkey fecal microbial community structure and functional gene diversity revealed by 16S rRNA and metagenomic sequencing. J. Microbiol. 46, 469–477. doi: 10.1007/s12275-008-0117-z
Muegge, B. D., Kuczynski, J., Knights, D., Clemente, J. C., González, A., Fontana, L., et al. (2011). Diet drives convergence in microbiome functions across mammalian phylogeny and within humans. Science 332, 970–974. doi: 10.1126/science.1198719
Otani, S., Hansen, L. H., Sørensen, S. J., and Poulsen, M. (2016). Bacterial communities in termite fungus combs are comprised of consistent gut deposits and contributions from the environment. Microb. Ecol. 71, 207–220. doi: 10.1007/s00248-015-0692-6
Pikuta, E. V. (2014). “Section 2: Family Carnobacteraceae” in The Lactic Acid Bacteria: Biodiversity and Taxonomy, 1st edn, eds W. H. Holzapfel and B. J. B. Wood (Hoboken, NJ: John Wiley & Sons, Ltd), 107–108.
Poulsen, M., Hu, H., Li, C., Chen, Z., Xu, L., Otani, S., et al. (2014). Complementary symbiont contributions to plant decomposition in a fungus-farming termite. Proc. Natl. Acad. Sci. U.S.A. 111, 14500–14505. doi: 10.1073/pnas.1319718111
Pukall, R., Lipidus, A., Rio, T. G. D., Copeland, A., Tice, H., Cheng, J., et al. (2010). Complete genome sequence of Conexibacter woesei type strain (ID131577T). Stand. Genomic Sci. 2, 212–219. doi: 10.4056/sigs.751339
R Core Team (2013). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org/
Russell, J. A., Moreau, C. S., Goldman-Huertas, B., Fujiwara, M., Lohman, D. J., and Pierce, N. E. (2009). Bacterial gut symbionts are tightly linked with the evolution of herbivory in ants. Proc. Natl. Acad. Sci. U.S.A. 106, 21236–21241. doi: 10.1073/pnas.0907926106
Sekelja, M., Rud, I., Knutsen, S. H., Denstadli, V., Westereng, B., Næs, T., et al. (2012). Abrupt temporal fluctuations in the chicken fecal microbiota are explained by its gastrointestinal origin. Appl. Environ. Microbiol. 78, 2941–2948. doi: 10.1128/aem.05391-11
Stanley, D., Geier, M. S., Chen, H., Hughes, R. J., and Moore, R. J. (2015). Comparison of fecal and cecal microbiotas reveals qualitative similarities but quantitative differences. BMC Microbiol. 15:51. doi: 10.1186/s12866-015-0388-6
Teyssier, A., Rouffaer, L. O., Hudin, N. S., Strubbe, D., Matthysen, E., Lens, L., et al. (2018). Inside the guts of the city: urban-induced alterations of the gut microbiota in a wild passerine. Sci. Total Environ. 612, 1276–1286. doi: 10.1016/j.scitotenv.2017.09.035
Veinberg, S., McClay, K., Masuda, H., Root, D., Condee, C., Zylstra, D. J., et al. (2006). Biodegradation of ether pollutants by Pseudonocardia sp. strain ENV478. Appl. Environ. Microbiol. 72, 5218–5224. doi: 10.1128/AEM.00160-06
Vela, A. I., Casas-Diaz, E., Fernandez-Garayzabal, J. F., Serrano, E., Agusti, S., Porrero, M. C., et al. (2015). Estimation of cultivable bacterial diversity in the cloacae and pharynx in eurasian griffon vultures (Gyps fulvus). Microb. Ecol. 69, 597–607. doi: 10.1007/s00248-014-0513-3
Videvall, E., Strandh, M., Engelbrecht, A., Cloete, S., and Cornwallis, C. K. (2017). Measuring the gut microbiome in birds: comparison of faecal and cloacal sampling. Mol. Ecol. Resour. 18, 424–434. doi: 10.1111/1755-0998.12744
Wang, W., Cao, J., Li, J. R., Yang, F., Li, Z., and Li, L. X. (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
White, J., Mirleau, P., Danchin, E., Mulard, H., Hatch, S. A., Heeb, P., et al. (2010). Sexually transmitted bacteria affect female cloacal assemblages in a wild bird. Ecol. Lett. 13, 1515–1524. doi: 10.1111/j.1461-0248.2010.01542.x
Zehr, J. P., Jenkins, B. D., Short, S. M., and Steward, G. F. (2003). Nitrogenase gene diversity and microbial community structure: a cross-system comparison. Environ. Microbiol. 5, 539–554. doi: 10.1046/j.1462-2920.2003.00451.x
Keywords: symbiosis, microbiome, MiSeq amplicon sequencing, insectivores, omnivores, nutrition
Citation: Bodawatta KH, Sam K, Jønsson KA and Poulsen M (2018) Comparative Analyses of the Digestive Tract Microbiota of New Guinean Passerine Birds. Front. Microbiol. 9:1830. doi: 10.3389/fmicb.2018.01830
Received: 15 May 2018; Accepted: 23 July 2018;
Published: 10 August 2018.
Edited by:David William Waite, University of Auckland, New Zealand
Reviewed by:Meagan Dewar, Federation University, Australia
Kirsten Grond, University of Connecticut, United States
Shahid Karim, University of Southern Mississippi, United States
Copyright © 2018 Bodawatta, Sam, Jønsson and Poulsen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kasun H. Bodawatta, email@example.com