Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2008, p. 5905-5912, Vol. 74, No. 19
0099-2240/08/$08.00+0 doi:10.1128/AEM.00574-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Biology, University of Puerto Rico, San Juan, Puerto Rico 00931-3360,1 Center for Genome Sciences, Washington University School of Medicine, St. Louis, Missouri 63108,2 Department of Microbiology, New York University School of Medicine, 550 First Avenue, New York, New York 10016,3 High Performance Computing Facility, University of Puerto Rico, P.O. Box 23334, San Juan, Puerto Rico 00931-3334,4 Department of Mathematics, College of Natural Sciences, University of Puerto Rico, San Juan, Puerto Rico 00931-3360,5 Venezuelan Institute of Scientific Research (IVIC), Caracas, Venezuela6
Received 10 March 2008/ Accepted 11 July 2008
|
|
|---|
75 known phyla in Bacteria were identified in this gut habitat; the Firmicutes were dominant (67% of sequences, belonging to the classes Clostridia, Mollicutes, and Bacilli), followed by the Bacteroidetes (30%, mostly in the order Bacteroidales), Proteobacteria (1.8%), and Lentisphaerae, Verrucomicrobia, TM7, Spirochaetes, Actinobacteria, and Aminanaerobia (all <0.1%). The novelty in this ecosystem is great; 94% of the phylotypes were unclassified at the "species" level and thus likely include novel cellulolytic lineages. |
|
|---|
![]() View larger version (38K): [in a new window] |
FIG. 1. Schematic representation of the hoatzin digestive tract. (A) Location of the crop and expanded esophagus in the hoatzin body. The anterior sternum is much reduced to make room for the large crop. (Reprinted from Natural History [[17] with permission of the publisher.) (B) Extended complete digestive tract of the hoatzin. (Courtesy of Alejandro Grajal.)
|
3 kg (6); a small body size restricts the luminal volume needed for fermentation, and energy demand scales to metabolic weight (w0.75) (26). In addition, a bird must generate enough energy to maintain its body temperature, which can be as high as 43.5°C (3). This paradox has been explained by considering that flying can facilitate extreme diet selectivity, which in turn could maximize fermentation rates and digestive efficiency (16). In this view, crop fermentation allows the hoatzin to be a unique browsing bird that exploits resources normally available only to specialized mammals and not to other avian species (16).
Foregut fermentation originated on several independent occasions during the evolution of herbivorous vertebrates (28, 29, 41). Species of the orders Artiodactyla, Edentata, and Primates have also evolved foregut fermentation (49). The Opisthocomiformes appeared in the Eocene, about 55 million years ago, according to fossil records from Argentina (5), at about the same time that primitive ungulates (Condylarths) migrated to South America (37). A browser, the hoatzin ancestor may have evolved during the Eocene at the same time as herbivorous mammalian browsers, such as extinct South American ungulates (37), arboreal primates (4), and sloths (48).
Almost two decades have passed since the first report describing foregut fermentation in the hoatzin appeared (16). Nonetheless, the microbial composition of the crop of this bird has remained elusive. In this study, we characterized the crop bacteria of the hoatzin using culture-independent molecular methods. The results reveal an unexpected degree of unclassified lineages that likely include novel cellulolytic bacteria.
|
|
|---|
1 month prior to DNA isolation.
Genomic DNA extraction.
DNA was extracted from crop contents using a modification of the method of Löffler et al. (32). The modification consisted of microbial disruption by bead beating (0.8 g of sterile 0.1-mm-diameter zirconium beads [Biospec Products, Bartlesville, OK] in 1 ml of TEN buffer [100 mM Tris, 10 mM EDTA; 2 M NaCl, pH 7.4] containing 200 mg of the contents from a single crop, with a bead beater [Biospec Products] set at 5,000 rpm for 1 min at room temperature). DNA samples were kept frozen until they were used.
PCR amplification of 16S rRNA genes.
PCR (two replicates per sample) were performed by using 50-µl mixtures containing 2 PuReTaq Ready-To-Go PCR beads (GE Healthcare Bio-Sciences Corp., Piscataway, NJ),
50 ng of DNA template, and 20 pmol of universal bacterial primers 8F and 1513R (38). PCR for two or three replicates per sample were performed using 30 cycles and an annealing temperature of 52°C, as described by Saitou and Nei (42).
Clone library construction and DNA sequencing.
PCR products were subsequently purified (PCR purification kit; Qiagen, Valencia, CA) and subcloned (pGEM-T Easy [Promega, Madison, WI], using Escherichia coli XL1-Blue [Stratagene, La Jolla, CA] or pCR4.0, a TOPO TA cloning kit, and TOP10 competent cells [Invitrogen, Carlsbad, CA]). Cloned amplicons were sequenced using vector-specific primers and an ABI 3730xl instrument (Applied Biosystems, Foster City, CA).
DNA sequence analyses and chimera detection.
Sequences were edited with Sequencher 4.6, aligned with the Greengenes database, and classified according to the Hugenholtz taxonomy (as of April 2007 [7]). Hoatzin crop 16S rRNA gene sequences and their closest neighbors in Greengenes were then imported into ARB (34). A BLAST script and the classification tool in Greengenes were used to determine similarity to GenBank sequences. Chimeras (n = 519) (checked with Bellerophon, version 3 [19]) and sequences corresponding to plant chloroplasts (n = 84) were excluded from further analysis.
Richness and coverage estimators.
To estimate richness, sequences were binned into "species-level" operational taxonomic units (OTUs). Since in species the level of similarity of 16S rRNA gene sequences can be as high as 98 to 99% (27), we chose a sequence identity level of 99% (with the hypervariable regions of the sequence masked using lanemaskPH [20] in DOTUR [43]). One representative sequence per OTU was randomly selected for calculating a phylogenetic tree using the neighbor-joining algorithm (42) in ARB (34). OTU abundance values were used to create a heat map for interindividual comparisons (TreeView) (http://rana.lbl.gov).
Collector's curves for observed and estimated (Chao1 and CACE) OTUs were computed using DOTUR (43) and EstimateS (R. K. Colwell, EstimateS: statistical estimation of species richness and shared species from samples, version 7, 2004 [http://viceroy.eeb.uconn.edu/estimates; www.purl.oclc.org/estimates]). Bacterial diversity was estimated with Shannon and Simpson indexes in DOTUR; the Shannon index takes into account the number and evenness of species (2), while the Simpson index estimates the probability that two randomly selected individuals belong to the same species (46). The Pielou evenness index (39) was also calculated for each animal's microbial community. Good's coverage index (12) was estimated using singleton sequences obtained from DOTUR.
Bacterial community comparisons.
Bacterial communities from the different hoatzin crops were compared using UniFrac (33). UniFrac analyses were based on an ARB neighbor-joining tree of 16S rRNA gene sequences, their frequencies, and their assigned environments (individual crops). Phylum-level comparisons between the hoatzin crop communities and the communities in other fermentative organs were also performed using previously described 16S rRNA gene data sets in which the sequences were more than 500 bp long and data sets containing more than 150 sequences from particular host species, including 425 sequences from the chicken intestine (53; P. T. Lan, M. Sakamoto, and Y. Benno, unpublished data from GenBank), 752 sequences from cows (25, 51; L. Cauquil, S. Combes, V. Montelis, J. Gordon, and T. Gidenne, unpublished data from GenBank; E. C. Shin, W. J. Lim, H. Kim, and H. D. Yun, unpublished data from GenBank), 176 sequences from yaks (1; S. Cai and X. Dong, unpublished data from GenBank), 180 sequences from deer (35, 50), 436 sequences from sheep (30, 40; S. Sawanon, S. Koike, and Y. Kobayashi, unpublished data from GenBank; T. Shinkai, N. Matsumoto, and Y. Kobayashi, unpublished data from GenBank), and 165 sequences from water buffalo (H. Mao, S. Ma, J. Chen, and W. Deng, unpublished data from GenBank).
Nucleotide sequence accession numbers.
Bacterial 16S rRNA gene sequences obtained in this study have been deposited in the GenBank database. The accession numbers for the 580 representative OTUs are EU344158 to EU344737, and the accession numbers for the remaining 655 cloned sequences are EU747884 to EU748538.
|
|
|---|
Sample representativeness and bacterial richness.
We obtained a total of 1,235 bacterial 16S rRNA gene sequences (166 to 332sequences/bird) (Table 1), which were binned into 580 OTUs (threshold cutoff for each OTU,
99% nucleotide sequence identity). Most sequences were full-length sequences (1,400 ± 50 bp); the exceptions were animal 4 sequences, whose average length was 800 bp. Microbial communities are complex, particularly those associated with the vertebrate gut; thus, 16S rRNA gene surveys seldom sample the complete diversity. Therefore, for each bird sample, we assessed richness (actual diversity) using the Chao1 and CACE estimators and coverage (how well a sample represents the population) using Good's coverage estimator. We chose 99% sequence identity (with hypervariable regions removed) as a cutoff for differentiation of "species" (10). The results (Table 1) indicate that, on average, 67% of the pooled species richness (55% for each individual) was sampled.
|
View this table: [in a new window] |
TABLE 1. Diversity indices for bacterial communities in individual hoatzin crops
|
![]() View larger version (27K): [in a new window] |
FIG. 2. Phylogenetic analysis of the hoatzin crop bacteria. (A) Tree of the 580 identified bacterial OTUs, showing relative clone abundance for six individual animals. In the neighbor-joining tree, the lines on the left indicate OTUs and the columns on the right show data for animals. The phyla are color coded as described in Fig. S1 in the supplemental material and (from top to bottom) are Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, Spirochaetes, Lentisphaerae, Verrucomicrobia, TM7, and Aminanaerobia. The relative clone abundance for each OTU in each animal is indicated by a shade of gray (white, 0%; black, 100%). (B) Distribution of the 1,235 16S rRNA gene sequences in nine bacterial phyla for individual hoatzins. Data for the two phyla represented most often, Bacteroidetes and Firmicutes, are shown in the top panel, and data for the other seven phyla are shown in the bottom panel. Proteobacteria were the predominant bacteria in animal 4 but were absent in animal 2, whereas Verrucomicrobia and Lentisphaerae were detected only in animal 6.
|
![]() View larger version (21K): [in a new window] |
FIG. 3. Collector's and rarefaction curves of OTU richness for the 1,235 16S rRNA gene sequences from the pooled hoatzin crops. (A) Collector's curves of observed and estimated species-level OTU richness using a 99% identity cutoff value. The richness was estimated to be 1,795 OTUs using the Chao1 estimator and 1,709 OTUs using the CACE estimator. (B) Rarefaction at identity cutoff values of 90 to 99%. There was a notable increase in the number of OTUs with only a small increase in the cutoff value (580 OTUs at a cutoff value of 99% and 470, 402, 226, and 152 OTUs at cutoff values of 98, 97, 93, and 90%, respectively), reflecting greater richness at lower phylogenetic branching depth.
|
75 known phyla in the domain Bacteria (see Fig. S1 in the supplemental material), similar to the results for the rumens of cow and sheep (12 and 11 phyla, respectively) (see Fig. S2 in the supplemental material). Bacteria belonged principally to the Firmicutes (67% of sequences) and the Bacteroidetes (30%) (see Fig. S1 in the supplemental material). Other phyla represented included the Proteobacteria (1.8%) and the Lentisphaerae, Verrucomicrobia, TM7, Spirochaetes, Actinobacteria, and Aminanaerobia (Synergystes plus Dethiosulfovibrio) (all <0.1%). Remarkably, most OTUs exhibited less than 95% identity with any previously described 16S rRNA gene sequence (Fig. 4); of the 580 OTUs, 545 (94%) were unclassified at the species level, while 327 (56%), 147 (25%), 70 (12%), and 46 (8%) were unclassified at the genus, family, order, and class levels, respectively (Table 2). Within the Firmicutes phylum, the class Clostridia was the most abundant (329/396 OTUs; 83%); it was represented by the order Clostridiales almost exclusively, including the family Lachnospiraceae (248 OTUs, 75%), candidate group RC6 (46 OTUs, 16%), the family Peptostreptococcaceae (31 OTUs, 9.5%), and the family Clostridiaceae (4 OTUs, 1%). The class Mollicutes was represented by 27 OTUs (7%), while there were only three Bacilli OTUs (<1%). The Bacteroidetes phylum contained 159 OTUs, all of which belonged to order Bacteroidales, including the family Prevotellaceae (82 OTUs, 51%), candidate group rf14 (36 OTUs, 23%), and a few other candidate groups in the Hugenholtz taxonomy (http://greengenes.lbl.gov/cgi-bin/nph-classify.cgi). The Proteobacteria contained 12 OTUs, 11 of which were Gammaproteobacteria or Betaproteobacteria (the remaining OTU belonged to the Epsilonproteobacteria). The OTUs whose best match was a rumen bacterium (49% of the OTUs) had 16S rRNA gene sequence identities ranging from 88 to 95% with bacteria such as Butyrivibrio fibrisolvens, Eubacterium cellulosolvens, or Bacteroides xylanolyticus (see Table S1 in the supplemental material).
![]() View larger version (17K): [in a new window] |
FIG. 4. Distribution of species-level OTU frequencies in relation to known sequences. Most of the hoatzin crop bacterial OTUs have 95% identity to previously reported sequences.
|
|
View this table: [in a new window] |
TABLE 2. Numbers of unclassified crop bacterial OTUs for each phylum
|
|
|
|---|
95% identity with the main rumen cellulolytic genera Ruminococcus, Bacteroides, Prevotella, and Butyrivibrio (22, 23). Furthermore, no members of the Fibrobacteria phylum, which includes cellulolytic bacteria common to rumens, were identified. The conspicuous absence of known cellulolytic bacteria and the discovery of many novel unclassified sequences related to rumen bacteria raise the possibility that many of the unclassified lineages may be cellulolytic. The hoatzin's body temperature is
39°C (13), while Artyodactyla have body temperatures below this value (3). The gut temperature is likely to increase with fermentative activity; however, the crop temperature is 39°C ± 1°C (n = 4 animals) (our unpublished observation).
Despite the partial coverage obtained in this survey, the 580 OTUs (or even the 402 OTUs that could be classified using a less stringent threshold cutoff,
97% 16S rRNA gene sequence identity) represent greater richness than that described in the more deeply sampled human colon (10). The majority of the OTUs were unique in each individual, perhaps due in part to the partial sample coverage. The high level of bacterial diversity is likely related to the high dietary diversity of the hoatzin, which feeds on leaves of members of the Fabaceae, Sterculiaceae, Rutaceae, and Vitaceae (9).
While cellulose digestion is a very important function of the rumen (52), cellulose is not a major digested polysaccharide in the crop (8). Hemicellulose degradation appears to be more important than cellulolysis in the hoatzin crop (24), as it is in other selective browsers, such as the howler monkey (36), which chooses young leaves and buds with thin cell walls and high cellular contents (9). However, the metabolic activities of the hoatzin crop microbiota remain elusive; the surprising diversity and novelty suggest that comparative metagenomic sequencing of whole-community microbial DNA prepared from the hoatzin and other foregut fermenters would be a worthwhile undertaking.
The results may provide important new insights into (i) how gut microbial communities are affected by diet and host phylogeny, (ii) the degree to which glycobiomes (genes involved in the acquisition and metabolism of carbohydrates) vary among communities that consist of seemingly divergent collections of organismal lineages, and (iii) how fermentation by the crop microbiota may provide sufficient energy from both cellular contents and cell wall polysaccharides to satisfy the hoatzin's energetic requirements.
We acknowledge Daniel Ayala for providing the BLAST script and Falk Warnecke, Todd DeSantis, Philip Hugenholtz, and Leslie Dethlefsen for helpful suggestions concerning the use of sequence analysis software. We also thank Alejandro Grajal for providing his drawing of the hoatzin digestive tract. The logistic field support provided by Hato Piñero personnel and José and Antonio González from Hato Mataclara and the collection permits provided by the Venezuelan Ministry of Environment are deeply appreciated.
Published ahead of print on 8 August 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»