Age-related variation in the oral microbiome of urban Cooper’s hawks (Accipiter cooperii)
BMC Microbiology volume 19, Article number: 47 (2019)
Bird species worldwide are affected by trichomoniasis caused by the protist Trichomonas gallinae. In avivorous raptors such as Cooper’s hawks (Accipiter cooperii), nestlings are more susceptible than fledglings and adults. Previous research suggested a link between oral pH and susceptibility: the oral pH of fledgling and adult hawks is more than seven times more acidic than that of nestlings. We speculated that this age-specific difference in pH would correspond to age-specific differences in the oral microbiota of Cooper’s hawks. We examined the oral microbiomes of 31 healthy, wild Cooper’s hawks in Tucson, Arizona (USA). Individuals represented three age classes (nestlings, fledglings, and adults). We designed our study with multiple controls, replicated sampling, mock communities, and stringent quality-controls to address challenges that can limit the inferential quality of microbiome data sets.
Richness of bacterial communities in oral cavities of Cooper’s hawks differed as a function of age but not as a function of sex, sampling date, or sampling location. Bacterial communities in oral cavities of nestlings differed from those of fledglings and adults, whereas communities in fledglings and adults did not differ from each other. Communities were similar in males and females and did not differ over the sampling season. Prevalence of acid-producing bacteria in fledgling and adults vs. nestlings is consistent with previous reports of age-specific variation in oral pH, but further research is needed to establish a causal link to pH levels or susceptibility to disease. Analyses of mock communities demonstrated high repeatability and showed that operon number and read abundance were highly correlated.
The oral microbiota of wild Cooper’s hawks differs between nestlings and older birds. Variation in the oral microbiome is consistent with differences in oral pH between nestlings and older individuals. Overall our study provides a first perspective on bacterial communities associated with oral cavities of a wild raptor.
Bird species worldwide are affected by diverse diseases including trichomoniasis, a disease caused by the flagellated protist Trichomonas gallinae [1, 2]. Birds of prey that consume members of the family Columbidae are especially susceptible to trichomoniasis, as pigeons and doves are primary hosts for T. gallinae [3, 4]. Consequently, infectious diseases such as trichomoniasis rank among the most common causes of morbidity in some birds of prey . Susceptibility to trichomoniasis varies among age classes in several species of predatory birds. In general, young birds are more likely to contract the disease than older birds, often leading to high rates of nestling mortality [2, 6, 7]. For example, trichomoniasis was identified as the cause of mortality in ca. 40% of nestlings produced annually in a population of urban-nesting Cooper’s hawks in Tucson, Arizona . For susceptible Cooper’s hawks, infection by T. gallinae results in lesions in the mouth and oropharynx, with subsequent spread of lesions to the esophagus and crop . Reduced food consumption and intermittent regurgitation resulting from severe lesions leads to weight loss and ultimately can be fatal .
Previous research suggests that susceptibility of Cooper’s hawks to trichomoniasis is related strongly to the age-specific pH of the oral cavity . The mean pH of fluid in the oral cavity of nestling Cooper’s hawks is 6.8, whereas that of fledglings and adults is 6.0–6.1 (i.e., is at least seven times more acidic) . Trichomonads generally are sensitive to environmental pH (e.g., T. gallinae , T. gallinae and T. vaginalis , T. vaginalis , T. foetus ). Trichomonas gallinae thrives when pH is between 6.5 and 7.5 (optimum 7.2 [11, 12]). When the pH is near optimum for trichomonads, some species (e.g., T. vaginalis) secrete a substance that facilitates invasion of host epithelial cells , but the ability to do so declines rapidly when pH is above or below the optimum.
Fluid in the oral cavity of Cooper’s hawks becomes more acidic after birds have fledged and are nearing independence (≥50 days of age; ), but the reason for the change is unknown. Many animals undergo similar changes in body chemistry (i.e., a change in acidity) during maturation, and those changes often are associated with changes in their bacterial communities [16,17,18,19,20,21,22]. We speculated that the community of bacteria in the oral cavity of nestling Cooper’s hawks differs from that in older hawks. Our objectives in this study were to provide a first perspective on the bacterial communities associated with the oral cavity of a wild raptor, and to compare the oral microbiomes of apparently healthy Cooper’s hawks in three age classes, sampled contemporaneously from a shared geographic area.
We conducted the field portion of the study in Tucson, Arizona from April to July 2015. We sampled the oral microbiome of hawks from three age groups: nestlings (14–28 days old), fledglings (55–120 days old), and breeding adults (> 1 year old). We used a bal-chatri trap  to capture fledglings and adult males, and a dho-gaza trap  to capture adult females. Nestlings were taken by hand from their nests. We differentiated males from females by diameter of the tarsometatarsus . All hawks were captured from nests monitored in a long-term study of population dynamics . The University of Arizona’s IACUC reviewed and approved all of the activities in our study that involved the capture and handling of adult, fledgling, and nestling Cooper’s hawks (Protocol 12–329). Sampling locations are shown in Table 1 and Additional file 1: Figure S1.
Sampling the oral microbiome
We used a sterile foam-tipped applicator (Whatman, model WB100032) to swab the tongue and palate of each hawk. Prior to swabbing, we used scissors to trim the foam tip so it fit easily in the oral cavity. Scissors were soaked thoroughly with 95% ethanol prior to each use. We wore latex gloves and used a wooden tongue depressor to open the mouths of adult hawks. We used a new depressor for each hawk. Immediately after swabbing, we stored each foam tip in a sterile, 2 mL Eppendorf tube. We stored all tubes at − 80 °C within 2 h of sampling. In total we sampled 34 hawks; of these, 31 yielded high-quality DNA extractions and were included in our analyses (Table 1).
We used two types of controls in the field that allowed us to distinguish between bacterial associates of the oral cavity of hawks and those that might be captured incidentally during our sampling process. First, at each sampling event we retained a foam applicator that had been trimmed with scissors and exposed to the air as above, but not exposed to a hawk (field, air control). Second, we used an applicator that was trimmed as described above to swab a section of each tongue depressor that had not contacted a hawk or a person (field, wood control). These field-control swabs were processed contemporaneously with oral swabs as described below.
We used half of each foam tip for DNA extraction, retaining the second half for archival purposes. We used the PowerSoil DNA Isolation Kit (Mo Bio Laboratories: Carlsbad, California, USA) to extract total genomic DNA from each half-tip. We followed the manufacturer’s instructions, except that we added an initial incubation for 10 min at 65 °C in lysis buffer. We used the Qubit 2.0 Fluorometer (Invitrogen: Carlsbad, California, USA) with the Qubit dsDNA HS (High Sensitivity) Assay Kit to quantify DNA concentrations following the manufacturer’s instructions.
Library preparation and Illumina sequencing
We used a dual-barcoded, two-step PCR approach for high throughput amplicon sequencing on the Illumina MiSeq platform (Illumina, Inc.) . In the first PCR (PCR1) we amplified the V4 region of the 16S rRNA with primer pair 515F/806R . Each PCR1 primer contained a universal 22 basepair (bp) consensus sequence tag (i.e., CS1 forward and CS2 reverse), 0–5 bp for phase-shifting, a 2 bp linker, and the locus-specific primer (515F or 806R). We pooled forward or reverse PCR1 primers with different phase-shifting lengths in equimolar concentrations prior to PCR, such that each amplification consisted of a random mixture of different phase-shifting lengths. This allowed for the maximum amount of sequence diversity in the first four bases, which is critical for accurate cluster identification and color matrix estimation on the Illumina MiSeq . We chose linkers and phase-shifting bases with low identity between primers and the target sequences.
In the second amplification step (PCR2) we used forward and reverse primers that each contained the complement of the CS tag, a 12 bp barcode, and the corresponding Illumina sequencing primer. The addition of barcodes in PCR2 allowed for maximum flexibility in the locus of interest or to include multiple targets in the same sequencing reaction, without needing to purchase a large number of barcoded target-specific primers.
We performed PCR1 in triplicate in 15 μl reaction volumes that each contained 0.5 μl DNA template, 7.5 μl 1X Phusion Flash High-Fidelity PCR Master Mix (ThermoFisher Scientific, Austin, Texas, USA), 0.2 μl of 50 μM forward and reverse primers (505F and 806R respectively) , 1 mg/mL of molecular grade bovine serum albumin (BSA; New England Biolabs, Ipswich, Massachusetts, USA), and 6.0 μl of molecular grade water. We used the following cycling protocol: an initial denaturing step at 98 °C for 10 s; 25 cycles consisting of denaturation at 98 °C for 1 s, annealing at 57 °C for 5 s, and extension at 72 °C for 20 s; and a final extension at 72 °C for 1 min. We used sterile, molecular grade water instead of template for negative controls.
We used SYBR Green 1 (Molecular Probes, Invitrogen) to visualize amplification on a 2% agarose gel after electrophoresis. We then pooled the three PCR1 products for each sample. From this, we diluted 5 μl of the pooled amplicons with molecular grade water to a final concentration of 1:15. We then used 1 μl of the pooled, diluted PCR1 product as the template for PCR2.
Each PCR2 reaction contained a final concentration of 1X Phusion Flash High-Fidelity PCR Master Mix, 0.075 μM barcoded primers (forward and reverse previously pooled at a concentration of 2 μM), and 0.24 mg/mL of BSA, for a final volume of 20 μL. We used the following cycling protocol: an initial denaturing step at 98 °C for 10 s; five cycles consisting of denaturation at 98 °C for 1 s, annealing at 51 °C for 5 s, and extension at 72 °C for 15 s; and a final extension step at 72 °C for 1 min.
We visualized products of PCR2 on a 2% agarose gel to verify amplification with minimal primer dimers. We quantified PCR2 products with PicoGreen and the Biotek Synergy H1 Multi-Mode Reader (Winooski, Vermont, USA). We normalized amplicons to 1 ng/μL and pooled 2 μL of each for Illumina sequencing. We purified the final amplicon pool with Agencourt AMPure XP beads (Beckman Coulter, Inc., Brea, California, USA) at a ratio of 1:1 to remove excess primers, nucleotides, salts, enzymes, and primer dimers, following the manufacturer’s instructions.
We evaluated the amplicon library on a BioAnalyzer 2100 (Agilent Technologies) to determine concentration and fragment size distribution. Paired-end sequencing was performed on an Illumina MiSeq with Reagent Kit v3 (2 × 300 bp) according to protocols at the IBEST Genomics Core at the University of Idaho.
We complemented our field controls (i.e., air and wood controls) with several precautionary measures in the laboratory to limit contamination. We prepared PCR mixes in a sterile, dedicated “pre-PCR” hood, which never was exposed to amplified DNA. We decontaminated pipettes and all surfaces in the pre-PCR hood with DNA Away (Molecular Bioproducts, Inc., San Diego, California, USA) prior to each use. We used a dedicated “post-PCR” hood for all pre-sequencing steps after PCR1 (i.e., PCR1 pooling, dilutions, addition of diluted products to PCR2 master mix, and PCR2 amplicon pooling). We used DNA Away and treated the hood and all equipment with ultra-violet light for a minimum of 30 min prior to each use. We used sterile, aerosol-resistant pipette tips at all steps to minimize cross-contamination of samples, and we used separate reagents, pipettes, tips, and consumables for pre- and post-PCR setup. We pooled negative controls from PCR1 and used them as template for reactions in PCR2 to ensure that no contamination occurred during pooling or PCR2 setup. Finally, although we detected no contamination, we combined 5 μL of each PCR negative control, extraction blank, and field control in a separate pool and subjected them to the same pre-sequencing treatment as positive amplicon pools. We sequenced these negative controls in parallel with our samples.
To measure error rate, determine the consistency of quality among replicates, assess the correlation of operon number and read count, and filter sequence data, we amplified and sequenced bacterial mock communities in parallel with our samples. The mock communities were obtained through BEI resources (American Type Culture Collection, Manassas, Virginia, USA), the National Institute of Allergy and Infectious Diseases, and the National Institutes of Health as part of the Human Microbiome Project. We used two preparations: Genomic DNA from Microbial Mock Community B (Even, Low Concentration), v5.1 L, for 16S rRNA Gene Sequencing, HM-782D; and Genomic DNA from Microbial Mock Community B (Staggered, Low Concentration), v5.2 L, for 16S rRNA Gene Sequencing, HM-783D. Each preparation contained 20 phylogenetically diverse bacterial taxa either at even concentrations of 100,000 copies per organism per μL (even mock community) or staggered concentrations from 1000 to 1000,000 copies per organism per μL (staggered mock community). We amplified 1 μL of each mock community for PCR1 as above. To assess within-run variation, we used 1 μL of diluted PCR1 product (even and staggered) for five separate PCR2 amplifications, resulting in amplicons with five different barcodes for each mock community. Overall, we sequenced 10 replicate mock samples in parallel with oral swabs, field controls, and laboratory controls.
We demultiplexed raw Illumina reads with custom scripts at the IBEST Genomics Core (https://github.com/msettles/dbcAmplicons) (one mismatch allowed in barcode, four mismatches allowed in primers with ends matching). Overall, we obtained 1,011,997 paired-end reads after demultiplexing. We used FastQC to assess the quality of demultiplexed reads (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/).
We used the UPARSE pipeline with USEARCH v.8.1.1861_i86linux32  for all subsequent analyses. We used the fastq_eestats2 command to create a summary report showing the number of reads that would pass an expected error filter (maxEE 0.25, 0.50, 1.0) at different length thresholds (length_cutoffs 200, 250, 10). Based on this report we chose a length cutoff that would yield a sufficient number of high quality reads per sample while also retaining maximum read length. We used the fastq_filter command to trim forward (R1) reads to 200 bp (fastq_trunclen = 200) and to remove reads with greater than one error rate (max_ee = 1). Overall, 952,214 reads passed this quality filtering (94%).
We used the command derep_fulllength to dereplicate quality-trimmed reads and then removed singletons (parameters -sizeout -minuniquesize 1). We clustered dereplicated sequences (48,221 sequences) into operational taxonomic units (OTU) at 97% sequence similarity. In addition to de novo chimera checking performed during clustering, we used the RDP classifier for reference-based chimera checking of representative sequences for each OTU (16S rRNA reference database v.9). We used the command–usearch_global (id = 0.97) to map raw reads back to 478 chimera-checked OTU. From the resulting OTU table we removed 48 OTU that represented potential contaminants (i.e., OTU present in the PCR negative controls (16), extraction blanks (4), or field controls (11), or spurious OTU found in the mock community (17)). We retained 430 OTU for analysis (Additional file 2: Table S1). Following quality control we retained high-quality data from 31 hawks (Table 1).
We used the RDP naïve Bayesian 16S rRNA classifier to estimate taxonomy of each OTU. We used analyses of variance (ANOVA) to compare read number among age classes (nestlings, fledglings, and adults), sampling dates (Table 1), and sampling locations (Table 1). We used a t-test to compare read number between male and female hawks. We used a Kruskal-Wallis test to compare richness among age classes. We used ANOVA to compare read number for replicate analyses of the even mock community. We used linear regression to compare read number and operon number for the staggered mock community. Normality was confirmed before all parametric tests. We used analyses of similarity (ANOSIM) based on presence-absence data (Jaccard’s Index) and abundance data (Simpson’s Index) to evaluate community composition, implementing two sets of analyses: one that included only the OTU represented in the data set by ≥1000 reads (55 OTU), and one that included all nonsingleton OTU (236 OTU). We used a Bonferroni correction to correct for multiple comparisons in each analysis and non-metric multidimensional scaling implemented in PAST  to visualize results.
We obtained an average of 20,345 sequences from each hawk (95% CI: 15,252–25,437 sequences). Read number did not differ among nestlings, fledglings and adults (n = 31 hawks; log-transformed, F2,28 = 0.5154, P = 0.6028). Read number was consistent between males and females (t29 = − 0.2844, P = 0.7782), among sample-collection dates (n = 10 dates, log-transformed, F9,21 = 1.3857, P = 0.2559), and among sampling locations (n = 15 sampling locations, log-transformed, F14,16 = 0.8039, P = 0.6560). This consistency in read number permitted direct comparisons of richness among sample classes.
We observed a mean of 84.5 bacterial OTU per individual hawk (95% CI: 79.3–89.7 OTU). One adult hawk was excluded from further analyses because richness from that sample was ≥2 standard deviations lower than the mean. Analysis of data from the remaining 30 hawks showed that bacterial richness differed as a function of age class (Fig. 1) (χ2 = 6.51, df = 2, P = 0.0386). Mean richness ranged from 92.3 OTU for adults (95% CI: 83.6–101.0) to lower values in nestlings and fledglings (respectively: 83.4 OTU, 95% CI: 77.4–90.5; 72.0 OTU, 95% CI: 54.6–89.4) (Fig. 1).
We analyzed residuals from the previous analysis to determine whether richness differed as a function of sex, sampling date, or sampling location once differences due to age class were taken into account. We observed no difference in richness between male and female hawks (t28 = − 0.3657, P = 0.7173), among sampling dates (F9,20 = 1.3524, P = 0.2730), or among sampling locations (F14,15 = 0.5887, P = 0.8356).
Analyses of mock communities demonstrated the consistency of results among replicates when operon counts were equal (even mock community, Additional file 3: Table S2), and a strong correlation of read number with operon number when inputs were unequal (staggered mock community, R2 = 0.88; Additional file 1: Figure S2). The latter observation led us to use read count as a proxy for abundance in subsequent analyses, permitting analyses based on both presence/absence data and abundance.
Overall, 55 OTU that passed the quality-control process were represented by > 1000 reads. None was present in any controls from the field or laboratory. These OTU were primarily members of Firmicutes (31.5% of OTU), Actinobacteria (22.2%), Proteobacteria (20.4%), Bacterioidetes (18.5%), and Tenericutes (3.7%), with one representative each of Fusobacteria and the SR1 clade (Table 2). Clostridia, Gammaproteobacteria, Actinobacteria, and Bacteroidia were numerically dominant, together accounting for 54% of OTU (Table 2). Porphyromonas was prevalent along with Micrococcineae, Corynebacterineae, Suttonella, Mycoplasma, and Actinomycineae (Table 2). Analyses of these OTU revealed that sampling was statistically complete for all age classes (Additional file 1: Figure S3), providing the basis for comparisons among age classes (below).
When analyzed in terms of bacterial OTU, bacterial communities from the oral cavities of nestlings differed from those of fledglings and adults, whereas communities in fledglings and adults did not differ from each other (Fig. 2). Results were consistent when based on presence-absence data (Jaccard’s Index: ANOSIM R = 0.2059, P = 0.0168; Bonferroni-corrected pairwise comparisons: fledglings vs. nestlings, P = 0.0414; adults vs. nestlings, P = 0.0128; adults vs. fledglings, P = 0.9399) and abundance data (Simpson’s Index: ANOSIM R = 0.1925, P = 0.0015; Bonferroni-corrected pairwise comparisons: fledglings vs. nestlings, P = 0.0031; adults vs. nestlings, P = 0.0083; adults vs. fledglings, P = 0.6955). Analyses based on the full data set without singletons yielded the same results (Additional file 1: Figure S4).
To assess the quality of our inference, we evaluated whether differences in community structure among age classes could instead reflect differences in the sex of individuals (i.e., males vs. females) or the timing of sampling events. However, analyses restricted to fledglings and adults revealed that the oral microbiota of female and male hawks was similar (analysis of 55 OTU, as above; Jaccard’s index, ANOSIM R = 0.0583, P = 0.1743; Simpson’s index, R = 0.0306, P = 0.3171). Bacterial communities in these age classes did not differ as a function of sampling month (Table 1; Jaccard’s Index, ANOSIM R = 0.1181, P = 0.1363; Simpson’s Index, R = 0.0880, P = 0.1931). Analyses based on the full data set without singletons yielded the same results (Additional file 1: Figure S4).
Read counts for 13 of the 55 most common OTU (i.e., 24.1%) were at least 2.5 times as common in fledglings than in nestlings (Table 2). Twelve of these were also at least 2.5 times as common in adults than in nestlings (Table 2). These 12 OTU included putative members of the genera Porphyromonas, Proteiniphilum, Veillonella, Peptoniphilus, Kingella, Parvimonas, and Suttonella, as well as representatives of the Corynebacteriaceae and Dermabacteraceae (Table 2). Seven additional OTU were at least 2.5 times as common in adults vs. nestlings but were not as common in the transitional fledgling stage, although sampling of the fledgling age class was relatively limited (Table 2).
The oral microbiota of urban Cooper’s hawks differs between nestlings and older individuals. Differences in the oral microbiota could not be attributed to sampling month, location, or the sex of birds. These differences correspond to previously detected differences in the oral pH of Cooper’s hawks . It is unclear whether the previously documented change in oral pH between nestlings and older birds could be attributed to an increase in the abundance of acid-producing bacteria as birds mature, or whether the shift is endogenous and selects for the establishment of acidophilic and/or acid-producing bacteria in older individuals. However, we observed that communities in the oral cavities of fledgling and adult hawks have high abundances of acid-producing bacteria, providing a basis for further study.
Comparison of the oral microbiome of Cooper’s hawks relative to that of other vertebrates
The oral microbiome of Cooper’s hawks included bacterial phyla that are prevalent in the oral cavity of other vertebrates. The most common OTU in the present study represented Firmicutes, Actinobacteria, Proteobacteria, and Bacteroidetes. Firmicutes, Proteobacteria, and Bacteroidetes are especially common in the oral microbiome of healthy dogs . These phyla also predominate in the oral microbiome of domestic cats [32, 33] and humans , wherein Spirochaetes, Actinobacteria, and Synergisetes also are common. As in those vertebrates, Porphyromonas was particularly prevalent in the oral microbiome of Cooper’s hawks, especially in older age classes (Table 2).
Previous studies on bird-associated microbiota have not focused on the oral cavity
Previous work on the microbiomes of birds has focused primarily on cloacal , gut [36,37,38,39,40], and facial microbiomes . Studies of vultures suggest an oral-gastrointestinal-fecal route for Clostridia and other microbial taxa, consistent with our observations of these bacteria in the oral cavity of Cooper’s hawks . Many of the dominant genera observed here are more consistently found in oral cavities of other vertebrates than in studies of other bodily microbiomes in birds [31,32,33,34,35,36,37,38,39,40].
Age-specific variation in the oral microbiota
Several studies focusing on microbial communities associated with birds have explored bacterial community composition as a function of age. For example, the gut microbiota of juvenile kakapos (Strigops habroptilus) differed from that of adults in the prevalence of Lactobacillus, but communities did not differ otherwise among individuals of different ages . Cloacal microbiomes differed significantly between chicks and adults of kittiwakes (Rissa tridactyla) . In poultry such as chickens (Gallus gallus domesticus), major changes in the cecal microbiome are observed several weeks after hatching relative to birds on the day of hatch . Among 59 neotropical bird species, the age of individuals appears to be less important than taxonomic affiliation in predicting the composition of the gut microbiome, suggesting that age-related structure is more easily observed within bird species, rather than among them .
We found that nestlings differed markedly in the composition of their oral microbiomes relative to adults and fledglings. This difference is striking given the feeding behavior of this species. Early in the nestling phase, Cooper’s hawks are fed by the adult female. She tears pieces of flesh from prey items delivered by the adult male and individually feeds each nestling. Late in the nestling phase, chicks are strong enough to dismantle prey delivered to the nest on their own. After fledging, prey items are delivered whole to fledglings who usually are perched near the nest. Prey species delivered to nestlings in Tucson are dominated by birds, 57% of which are members of the Columbidae . Breeding adults consume the same prey items they feed their nestlings, so it is unlikely that changes in diet among age classes play a role in the changes in oral microbiomes we observed.
Taxonomic composition and potential relevance to the pH of the oral cavity in Cooper’s hawks
Taxa such as Porphyromonas, Proteiniphilum, Parvimonas, Kingella, Suttonella, Peptoniphilus, and Veillonella were especially common in samples from fledglings and adult hawks but were less commonly observed in nestlings. Porphyromonas species are known to produce acidic products (e.g., butyric acid) . Proteiniphilum strains grow at a pH similar to that of the oral microbiome of fledgling and adult hawks (ca. 6.0) and can produce acetic acid . Parvimonas strains metabolize peptone and amino acids to form acetic acid . Acid production is known in Kingella, Suttonella, and Peptoniphilus species [46,47,48], and Veillonella strains can produce acetic and proprionic acids  (also characteristic of some Corynebacteriaceae, Table 2). Although it is unclear whether these bacteria are responsible for the age-related change in pH in the oral microbiome of Cooper’s hawks, the results are consistent with observed differences in pH reported previously . In future work, isolating culturable members of the oral microbiota may be useful to support experiments evaluating whether it is pH per se, or microbial community composition, that correlates with the age-specific differences in susceptibility to T. gallinae observed previously in Cooper’s hawks.
We designed our study to maximize quality control in the sampling, molecular analysis, and bioinformatics steps. Analyses of mock communities demonstrated the consistency of runs and a strong correlation of read number with operon number when inputs were unequal, as is anticipated for natural communities. Use of multiple field controls, replicated sampling within each age class, phylogenetically diverse mock communities prepared in two different ways, diverse measures to limit contamination in the laboratory, stringent quality-control thresholds, analyses that considered the most thoroughly sampled subset of the microbial community as well as the community as a whole, and analyses based on both presence/absence and read abundance allowed us to overcome some of the challenges that can restrict the inferential quality of microbiome data sets. Overall, our detection of robust differences in community structure among bacteria of the oral cavity in nestling hawks vs. fledglings and adults, and our testing of alternative explanations for patterns observed here (e.g., sex-related differences or differences as a function of sampling month) provide support for our conclusions.
Our observations of the age-structured oral microbiome of Cooper’s hawks are consistent with established, age-related differences in oral pH between nestlings and older individuals. Previous work has shown that distinctive microbial communities are associated with diseases caused by trichomonads (e.g., in the human vaginal microbiome ). Although further study is needed to evaluate causal relationships among microbiome composition, oral pH, and trichomoniasis, our study provides a first perspective on the bacterial communities associated with the oral cavity of a wild raptor and sets the stage for further research to explore how susceptibility to parasites such as T. gallinae may be placed in a microbial community context.
- 16S rRNA:
16S ribosomal ribonucleic acid
analysis of similarity
analysis of variance
operational taxonomic unit
polymerase chain reaction
United States of America
- USFWS ID:
United States Fish and Wildlife Service identification number
Cooper JE, Petty SJ. Trichomoniasis in free-living goshawks (Accipiter gentilis gentilis) from Great Britain. J Wildlife Dis. 1988;24:80–7.
Boal CW, Mannan RW, Hudelson KS. Trichomoniasis in Cooper’s hawks from Arizona. J Wildlife Dis. 1998;34:590–3.
Stabler RM. Variations in virulence of strains of Trichomonas gallinae in pigeons. J Parasit. 1948;34:147–9.
Erwin KC, Kloss C, Lyles J, Felderhoff J, Fedynich AM, Henke SE, Roberson JA. Survival of Trichomonas gallinae in white-winged dove carcasses. J Wildlife Dis. 2000;36:551–4.
Naldo JL, Samour JH. Causes of morbidity and mortality in falcons in Saudi Arabia. J Avian Med Surg. 2004;18:229–41.
Real J, Manosa S, Munoz E. Trichomoniasis in a Bonelli’s eagle population in Spain. J Wildlife Dis. 2000;36:64–70.
Krone O, Altenkamp R, Kenntner N. Prevalence of Trichomonas gallinae in northern goshawks from the Berlin area of northeastern Germany. J Wildlife Dis. 2005;41:3304–9.
Hedlund CA. Trichomonas gallinae in avian populations in urban Tucson, Arizona. In: Masters thesis; 1998.
Samour JH, Naldo JL. Diagnosis and therapeutic management of trichomoniasis in falcons in Saudi Arabia. J Avian Med Surg. 2003;17:136–43.
Urban EH, Mannan RW. The potential role of oral pH in the persistence of Trichomonas gallinae in Cooper’s hawks (Accipiter cooperii). J Wildlife Dis. 2014;50:50–5.
Cailleau R. La nutrition de Trichomonas columbae en culture. Comp Rendus Societe de Biologie. 1935.
Read CP. Comparative studies on the physiology of trichomonad protozoa. J Parasitol. 1957;43:385–94.
Garber GE, Bowie WR. The effect of Trichomonas vaginalis and the role of pH on cell culture monolayer viability. Clin Invest Med. 1990;13:71–6.
Daniel GE. Observations on the pH of a culture medium and on cultures of Trichomonas foetus. J Parasitol. 1948:496–504.
Fiori PL, Addis MF, Sechi A, Cappuccinelli. Trichomonas vaginalis haemolysis: pH regulates a contact-independent mechanism based on pore-forming proteins. Microb pathogenesis 1996;109–118.
Grice EA, Segre JA. The skin microbiome. Nat Rev Microbiol. 2011;9:244–53.
Belkaid YB, Segre JA. Dialogue between skin microbiota and immunity. Science. 2014;346:954–9.
Beasely DE, Koltz AM, Lambert JE, Fierer N, Dunn RR. The evolution of stomach acidity and its relevance to the human microbiome. PLoS One. 2015;10:e0134116. https://doi.org/10.1371/journal.pone.0134116.
Hickey RJ, Zhou X, Settles ML, Erb J, Malone K, Hansmann MA, Shew ML, Van Der Pol B, Fortenberry JD, Forney LJ. Vaginal microbiota of adolescent girls prior to menarche resemble those of reproductive-age women. MBio. 2015. https://doi.org/10.1128/mBio.00097-15.
Uchihashi M, Bergin IL, Bassis CM, Hashway SA, Chai D, Bell JD. Influence of age, reproductive cycling status, and the menstruation on the vaginal microbiome in baboons (Papio anubis). Am J Primatol. 2015;77:563–78.
Cundell AM. Microbial ecology of the human skin. Microb Ecol. 2016. https://doi.org/10.1007/s00248-016-0789-6.
Jasarevic E, Morrison KE, Bale TL. Sex differences in the gut microbiome-brain axis across the lifespan. Philos T Roy Soc B. 2016;doi:https://doi.org/10.1098/rsto.2015.0122.
Bloom PH, Clark WS, Kidd JW. Capture techniques. In: Bird DM, Bildstein KL, editors. Raptor research and management techniques. Surrey, British Columbia, Hancock house: Publishers; 2007. p. 193–219.
Hill NP. Sexual dimorphism in the Falconiformes. Auk. 1944;61:228–34.
Mannan RW, Steidl RJ, Boal CW. Identifying habitat sinks: a case study of Cooper’s hawks in an urban environment. Urban Ecosyst. 2008;38:26–34.
U'Ren JM, Arnold AE. Illumina MiSeq dual-barcoded two-step PCR amplicon sequencing protocol. Protocols.io. 2016; https://doi.org/10.17504/protocols.io.fs9bnh6.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. P Natl Acad Sci USA. 2010;108:4516–22.
Wu L, Wen C, Quin Y, Yin H, Tu O, Van Nostrand JD, Tong Y, Yuan M, Deng Y, Zhou J. Phasing amplicon sequencing on Illumina Miseq for robust environmental microbial community analysis. BMC Microbiol. 2015;15:125.
Edgar RC. UPARSE: highly accurate OUT sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.
Hammer O, Harper DAT, Ryan PD. PAST: paleontological statistics software package for education and data analysis. Paleontol Electron. 2001;4.
Sturgeon A, Stull JW, Costa MC, Weese JS. Metagenomic analysis of the canine oral cavity as revealed by high-throughput pyrosequencing of the 16S rRNA gene. Vet Microb. 2013;162:182–8.
Dewhirst FE, Klein EA, Bennett M-L, Croft JM, Harris SJ, Marshall-Jones ZV. The feline oral microbiome: a provisional 16S rRNA gene based taxonomy with full-length reference sequences. Vet Microb 215. 175:294–303.
Adler JC, Malik R, Browne GV, Norris JM. Diet may influence the oral microbiome composition in cats. Microbiome. 2016;4:23.
Dewhirst FE, Chen T, Izard J, Paster BJ, Tanner ACR, Yu W-H, Lakshamanan A, Wade WG. The human oral microbiome. J Bacteriol. 2010;192:5002–17.
Van Dongen WFD, White J, Brandi HB, Moodley Y, Markling T, Leclaire S, Blanchard D, Danchin E, Hatch SA, Wagner RA. Age-related differences in the cloacal microbiota of a wild bird species. BMC Ecol. 2013;13:11.
Hird SH, Sanchez C, Carstens BC, Brumfield RT. Comparative gut microbiota of 59 neotropical bird species. Front Microbiol. 2015;6:1403.
Oakley BB, Buhr J, Ritz CW, Kiepper BH, Berrang ME, Seal BS, Cox NA. Successional changes in the chicken cecal microbiome over 42 days of growth are independent of organic acid feed. BMC Vet Res 2014;10:282.
Waite DW, Taylor MW. Exploring the avian gut microbiota: current trends and future directions. Front Microbiol. 2015;6:673.
Pan D, Yu Z. Intestinal microbiome of poultry and its interaction with host and diet. Gut Microbes. 2013;5:108–19.
Roggenbuck M, Schnell IB, Blom N, Baelum J, Bertelsen MF, Sicheriz-Ponten T, Sørensen SJ, Gilbert MTP, Graves GR, Hansen LH. The microbiome of New World vultures. Nature Comm. 2014;5. https://doi.org/10.1038/ncomms6498.
Waite DW, Eason DK, Taylor MW. Influence of hand rearing and bird age on the fecal microbiome of the critically endangered kakapo. Appl Environ Microbiol. 2014;80:4650–8.
Estes WA, Mannan RW. Feeding behavior of Cooper’s hawks at urban and rural nests in southeastern Arizona. Condor. 2003;105:107–16.
Shah H. The genus Porphyromonas. In: The prokaryotes: a handbook on the biology of bacteria: ecophysiology, isolation, identification, applications, vol. IV. 2nd ed. New York: Balows A, Trüper HG, Dworkin M, Harder W, Schleifer K-H. Springer Science+Business Media; 1992.
Chen S, Dong X. Proteiniphilum acetatigenes gen. Nov., sp. nov., from a UASB reactor treating brewery wastewater. Int J Syst Microbiol. 2005;55:2257–61.
The Editorial Board. Parvimonas Bergey’s Manual of Systematics of Archaea and Bacteria. 2015:1–4.
Bonacorsi S, Bidet P, Yagupsky P. Microbiology, genomics, and population structure. In: St Geme J, editor. Advances in understanding Kingella kingae. New York: Springer; 2016.
Foster G, Malnick H, Lawson PA, Kirkwood J, MacGregor SK, Collins MD. Suttonella ornitholoca sp. nov., from birds of the tit families, and emended description of the genus Suttonella. Int J Syst Evol Microbiol. 2005;55:2269–72.
Song Y, Liu C, Finegold SM. Peptoniphilus gorbachii sp. nov., Peptoniphilus olsenii sp. nov., and Anaerococcus murdochii sp. nov. isolated from clinical specimens of human origin. J Clin Microbiol. 2007;45:1746–52.
Rogosa M. The genus Veillonella I. General cultural, ecological, and biochemical considerations. J Bacteriol. 1964;87:162–70.
Brotman RM, Bradford LL, Conrad M, Gajer P, Ault K, Peralta L, Forney LJ, Carlton JM, Abdo Z, Ravel J. Association between Trichomonas vaginalis and vaginal bacterial community composition among reproductive-age women. Sex Trans Dis. 2012;39:807–12.
We thank the Undergraduate Biology Research Program at The University of Arizona for supporting a research opportunity for MJT. We thank M-M. Lee, T.K. O’Connor, and N.B. Zimmerman for technical assistance in the laboratory, K. Youens-Clark, B.L. Hurwitz, and A. Gerritsen for computational assistance, IBEST Genomics Core at the University of Idaho for technical support and sequencing, and P.M. Laskus and L. McHugh for assistance in the field.+
Availability of data and materials
The molecular data generated in the current study are available in the Short-read Archive at GenBank, https://www.ncbi.nlm.nih.gov/sra (SRA BioProject PRJNA507286). All other data are in the supplementary files accompanying this manuscript.
Ethics approval and consent to participate
The University of Arizona’s IACUC reviewed and approved all of the activities in our study that involved the capture and handling of adult, fledgling, and nestling Cooper’s hawks (Protocol 12–329). All specimens were collected locally in Tucson, Arizona and followed institutional and local standards for research. It was required that the number of hawks handled be reported to the IACUC at the end of each year.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Sampling locations in Tucson, Arizona, USA. Points are numbered by representative hawks, corresponding to hawk ID numbers shown in Table 1. Scale, 0.6 cm = 1 km. Figure S2. Strong positive correlation between read count and operon number for staggered mock community, used to confirm that read abundance could be used as a proxy for OTU abundance. Figure S3. Species-accumulation curve for the 55 most common OTU indicates thorough sampling of bacterial communities, providing a basis for statistical analyses presented in the text. Figure S4. Communities of bacteria in the oral cavity differed as a function of age class in Cooper’s hawks. Non-metric multi-dimensional scaling analyses of all nonsingleton OTU that passed quality control for nestlings, fledglings, and adults. Panel A, Jaccard’s Index; panel B, Simpson’s Index. (DOCX 1850 kb)
Table S1. Summary of operational taxonomic units for oral microbiome of urban Cooper’s hawks. Excel file containing operational taxonomic units, taxonomy, and occurrence of each in individual hawks. (XLSX 115 kb)
Table S2. Mock community, number of operons, and average reads observed for the even and staggered mock communities. Table containing details for each mock community. (XLSX 56 kb)