- Research Article
- Open Access
Does solar irradiation drive community assembly of vulture plumage microbiotas?
Animal Microbiome volume 2, Article number: 24 (2020)
Stereotyped sunning behaviour in birds has been hypothesized to inhibit keratin-degrading bacteria but there is little evidence that solar irradiation affects community assembly and abundance of plumage microbiota. The monophyletic New World vultures (Cathartiformes) are renowned for scavenging vertebrate carrion, spread-wing sunning at roosts, and thermal soaring. Few avian species experience greater exposure to solar irradiation. We used 16S rRNA sequencing to investigate the plumage microbiota of wild individuals of five sympatric species of vultures in Guyana.
The exceptionally diverse plumage microbiotas (631 genera of Bacteria and Archaea) were numerically dominated by bacterial genera resistant to ultraviolet (UV) light, desiccation, and high ambient temperatures, and genera known for forming desiccation-resistant endospores (phylum Firmicutes, order Clostridiales). The extremophile genera Deinococcus (phylum Deinococcus-Thermus) and Hymenobacter (phylum, Bacteroidetes), rare in vertebrate gut microbiotas, accounted for 9.1% of 2.7 million sequences (CSS normalized and log2 transformed). Five bacterial genera known to exhibit strong keratinolytic capacities in vitro (Bacillus, Enterococcus, Pseudomonas, Staphylococcus, and Streptomyces) were less abundant (totaling 4%) in vulture plumage.
Bacterial rank-abundance profiles from melanized vulture plumage have no known analog in the integumentary systems of terrestrial vertebrates. The prominence of UV-resistant extremophiles suggests that solar irradiation may play a significant role in the assembly of vulture plumage microbiotas. Our results highlight the need for controlled in vivo experiments to test the effects of UV on microbial communities of avian plumage.
Avian evolution has been shaped by millions of years of natural selection by bacteria, archaea, fungi, and metazoan ectoparasites [1,2,3,4,5,6] that digest or degrade the keratinous feathers, nails, and foot, leg, and bill coverings of birds [7,8,9]. Feathers are largely composed of polymeric filaments of β-keratin that are tightly bound to an amorphous polymer matrix. Once grown, feathers are metabolically inert and have no vascular connection to the body. The bulk of microbial research on avian plumage has concentrated on the industrial potential of keratinolytic bacteria isolated from soil and commercial poultry waste to convert feathers [8, 9] into animal feed [10,11,12,13,14,15,16]. Microbial keratinases target cross-linked structural peptides that make feather keratin insoluble . The continued focus on domestic poultry and keratinolytic bacteria has resulted in a substantial void in our knowledge of the taxonomic diversity, host specificity, and assembly of microbial communities in the plumage of the 10,135 wild bird species [5, 6, 17,18,19,20,21,22,23,24,25,26,27,28].
Birds have evolved a suite of behaviours such as sunning, dusting, and anting to maintain plumage quality [29,30,31]. Avian sunning is phylogenetically and geographically widespread and serves a multiplicity of purposes [29, 32, 33]. Sunning birds typically spread their wings and tail for maximal exposure to the sun [29, 33]. Heat absorption and thermoregulation are the most frequently hypothesized functions [29, 33,34,35], but sunning is often performed at thermoneutral temperatures in tropical and temperate latitudes. In fact, birds may voluntarily sunbathe to the brink of hyperthermia [29, 36], indicating some purpose other than regulation of body temperature. Other traditional explanations for sunning include plumage-drying [29, 35], inhibition of feather-chewing lice [30, 36], restoration of feather shape , moult facilitation , production of vitamin D , and preen gland stimulation . After 75 years of study, researchers have been unable to reach consensus on a primary driver that comfortably explains the occurrence of sunning in more than 50 taxonomic families and 21 orders of birds [29, 33].
The detection of keratinolytic bacteria in the plumage of a wide variety of wild bird species [5, 6] led to a new hypothesis for the function of avian sunning behavior . Saranathan and Burtt  studied the effect of sunlight (280–750 nm) on Bacillus licheniformis, a common soil-dwelling keratinolytic bacterium found widely on avian plumage . Based on in vitro experiments, they posited that sunning might regulate potentially harmful plumage microorganisms. Their hypothesis has failed to gain traction and there is little empirical evidence that solar irradiation alters the assembly and abundance of plumage microbiotas in wild birds. Recent studies, however, suggest that feather structure and pigmentation may affect ultraviolet (UV) reflectance and absorbance of solar radiation and may thus affect bacterial diversity and abundance [40, 41].
Bird plumage exposed to intense solar radiation constitutes one of the more extreme microbial environments found on or inside terrestrial and aquatic animals. For example, maximum plumage temperatures recorded in sunning brown-necked ravens (Corvus ruficollis) ranged from 49.2 °C at skin level to an extraordinary 83.9 °C at feather surfaces . Black plumage attains significantly higher temperatures than pale plumage in direct sun at low wind speeds [36, 42, 43]. The disproportionate number of bird species with black or dark gray plumage in arid regions  suggests that comparably high plumage temperatures are routine. Hundreds of diurnal bird species live in open-habitat formations or forage aerially in full sunlight, but few species experience as much solar irradiation as New World vultures , a monophyletic group (Cathartiformes) of obligate scavengers of vertebrate carrion [44,45,46]. Vultures engage in thermal soaring  and are renowned for conspicuous spread-wing sunning at roosts and loafing sites [35, 47]. Indeed, sunning appears to be a deeply-ingrained behaviour in the Cathartiformes. Sunning vultures perch upright and typically face directly toward or away from the sun to maximize the incidence of sunlight on their plumage (Fig. 1). Both dorsal and ventral surfaces of flight feathers are exposed to direct solar irradiation during spread-wing sunning as individuals shift positions. Ventral and dorsal contour (body) plumage also receives direct solar irradiation during sunning bouts.
We used 16S rRNA sequencing to investigate the taxonomic diversity and abundance of plumage microbiota in five species of sympatric vultures in Guyana: Sarcoramphus papa (king vulture), Coragyps atratus (black vulture), Cathartes aura (turkey vulture), Cathartes burrovianus (lesser yellow-headed vulture), and Cathartes melambrotos (greater yellow-headed vulture). We predicted that vulture plumage microbiotas would be exceptionally rich based on their natural history and the myriad potential transmission mechanisms. The downy plumage of nestling vultures is exposed to soil microbiota at terrestrial nest sites and carrion regurgitated by adults . Vertical transmission from adults to nestlings is enabled during brooding and feeding. In adults, horizontal transmission of plumage microbiota is facilitated by mating and allopreening behaviour [47, 48]. Adult plumage, facial skin, and legs and feet are fouled by the rich microbiotas of decomposing flesh and body fluids of vertebrate carrion . The habit of inserting their heads into body cavities undoubtedly explains the surprisingly diverse microbiotas of vulture facial skin . Intraspecific and interspecific transmission almost certainly takes place when multispecies feeding aggregations jostle for access to carcasses. Vulture plumage is also chronically exposed to soil and dust kicked up at feeding sites, transferred from feet to plumage during preening, and when vultures rub their plumage on soil . Adding to the long list of contamination routes, plumage microbiotas are inoculated by falling excrement at communal roosting sties and aerial particulates during flight. In sum, vulture plumage represents a model system for investigating the consequences of solar irradiation on rich microbial communities.
Vulture plumage hosts diverse microbial assemblages
We sequenced the hypervariable V4-V5 region of the 16S SSU rRNA gene isolated from dorsal and ventral feather tracts of 34 field-collected specimens (see Methods). We obtained 790 to 101,187 (median = 34,441) filtered sequences per plumage sample and a total of 2,722,291 sequences in the pooled dataset (Additional file 1). Only a single plumage sample was represented by fewer than 2000 sequences. We emphasize generic-level comparisons in our analyses because a majority of the amplicon sequence variants (ASVs) were unclassified at the species level in the SILVA database . We present cumulative sum scaled  and log2-transformed ASV normalized sequence counts (hereafter log2 CSS) to account for variation in sequencing depth among samples  and non-normal distribution of abundances among ASVs.
The filtered data contained 2670 ASVs, 631 genera, 239 families, 118 orders, 49 classes, and 24 phyla of Bacteria and Archaea (Additional files 2,3). Rarefaction curves for ASVs and microbial genera, families, and phyla reached asymptotes between 1000 and 2000 sequences per sample (Additional file 4). The majority of ASVs (85.4%) originated from four bacterial phyla: Firmicutes (29.7%), Proteobacteria (24.8%), Actinobacteria (15.6%), and Bacteroidetes (14.9%). Rank-abundance histograms of ASVs, genera, and families by phylum were strongly left-skewed (Fig. 2). The rarest 10 bacterial phyla contributed only 1.1% of the ASVs. Archaea (phyla Euryarchaeota and Thaumarchaeota) constituted a mere 0.2% of the ASVs. The plumage microbiota of the pooled sample of vulture species included 85 bacterial genera, 21 families, five orders, and two classes that were taxonomically unclassified in the SILVA  database (Additional file 1).
Vultures exhibit subtle species-specific differences in microbial assemblages
The observed number of microbial genera ranged from 9 to 152 per plumage sample, averaging 77.6 and 56.2 genera for dorsal and ventral plumage samples, respectively (Fig. 3a,c). Shannon-Wiener diversity indices exhibited a wide range of values (Fig. 3b,d). Principal coordinate analysis (PCoA) based on weighted UniFrac distances revealed only modest separation of dorsal plumage microbial communities of Coragyps atratus from those of the three species of Cathartes (Fig. 3e,f). The PCoA plot for ventral plumage displayed more spatial overlap, reflecting extensive sharing of dominant microbial ASVs. Between-group weighted UniFrac distances were significant for vulture species (PERMANOVA F4,56 = 3.78, R2 = 0.20, P < 0.001), marginally significant for plumage tracts (F1,56 = 2.51, R2 = 0.03, P = 0.04), but were non-significant for the interaction, vulture species × plumage tract (F4,56 = 0.90, R2 = 0.05, P = 0.57).
Extremophile bacteria are abundant in heavily melanized plumage
Rank-abundance plots of microbial genera were strongly left-skewed for dorsal and ventral plumage tracts (Fig. 4). Perhaps the most surprising finding of the 16S rRNA survey was the marked abundance of extremophiles in heavily-melanized plumage of Coragyps atratus and the three species of Cathartes. The genus Deinococcus (phylum Deinococcus-Thermus), an aerobic polyextremophile that is remarkably resistant to ionizing radiation, UV, and desiccation [54,55,56,57], constituted 6.3% of the log2 CSS normalized counts of filtered sequences in the pooled sample (Additional file 3). It was 1.7 × more abundant than the second most common genus, Clostridium sensu stricto 1 (phylum Firmicutes), which is commonly found in the vulture gut . The third most abundant genus, Hymenobacter (phylum Bacteriodetes) also exhibits extremophile traits .
Collectively, Deinococcus and Hymenobacter (of 631 genera) constituted 9.1% of the log2 CSS normalized counts of filtered sequences recovered from black plumage of Coragyps atratus and the three species of Cathartes. By comparison, Deinococcus was rare (rank-abundance <60th) and Hymenobacter was absent in the beige plumage of the sole specimen of Sarcoramphus papa. Neither Deinococcus nor Hymenobacter is known to be pathogenic in vertebrates [56, 58], nor have these extremophiles been identified as keratinolytic agents in wild bird populations [5, 18, 21, 39]. Deinococcus has been identified as a possible contaminant of laboratory reagents and DNA extraction kits . We ruled out contamination in our study because DNA was not detected in blank samples processed with the extraction and clean up kits. Additionally, only 5% of microbiota samples (n = 208) from other anatomical regions (e.g., facial swabs, esophagus, duodenum, colon) obtained from the same vulture specimens, and randomly interspersed with feather samples in the same Illumina MiSeq sequencing runs, yielded sequence classified as Deinococcus or Hymenobacter (unpublished data). Moreover, only 5 of 64 (non-plumage samples) represented by low DNA concentrations (< 1 ηg/μL) contained Deinococcus or Hymenobacter. High standardized sequence counts of Deinococcus and Hymenobacter were exclusive to melanized plumage samples of Coragpys atratus and the three Cathartes species.
Bacterial genera capable of forming desiccation-resistant endospores (phylum Firmicutes, order Clostridiales) were common in both black and beige vulture plumage. Genera which exhibit keratinolytic activity in vitro (Bacillus, Enterococcus, Pseudomonas, Staphylococcus, Streptomyces), and are known to occur in the plumage of wild birds [5, 18, 21, 39], collectively constituted only 4.0% of the CSS log2 normalized sequence counts in vulture plumage (Fig. 4, Additional file 3).
The five most abundant microbial genera comprised 17.5% of the CSS log2 normalized counts in the pooled sample of vultures (Fig. 5). Abundances of the extremophile genera Deinococcus and Hymenobacter (Fig. 5a,b) were strongly negatively correlated with ASV richness in dorsal and ventral plumage tracts. The abundance of three other common bacterial genera (Clostridium sensu stricto 1, Mobilicoccus, Clostridium sensu stricto 7) exhibited weak correlations with ASV richness (Fig. 5c,d).
Vulture plumage is exposed to a rich source pool of bacteria and archaea in carrion, soil, gut contents, aerial particulates, and other environmental sources . The majority of bacterial genera detected in vulture plumage are drawn from phyla that dominate avian gut microbiotas (Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes) [60,61,62,63]. The Chao1  estimate of generic richness in the combined dorsal and ventral plumage microbiotas was 956 genera (95% confidence interval, 858–1098 genera). A significant fraction (13.5%) of the 631 bacterial and archaeal genera detected in the 16S rRNA survey were taxonomically unclassified in SILVA database . More than 70% of the bacterial and archaeal genera detected in our analyses have not been previously reported in the associated tables and appendices of published papers that investigated plumage microbiota with 16S rRNA sequencing [5, 6, 17,18,19,20,21,22,23,24, 27]. Vulture plumage microbiotas can accurately be described as hyper-diverse and taxonomically cryptic.
Non-experimental data must be interpreted cautiously, but we infer the comparative abundance of UV-resistant genera Deinococcus and Hymenobacter in plumage microbiotas as a signature of solar irradiation. Deinococcus is environmentally widespread but relatively rare and patchily distributed in media ranging from soil, dung, and sewage to stratospheric dust [56, 57]. The prevalence of Deinococcus in gut and respiratory tract microbiotas of vertebrates is poorly understood but it appears to occur at low frequencies [50, 60,61,62,63, 65, 66], if it was detected at all. Higher abundances have been noted in the gut of bats  and crocodile lizards . Previous reports from avian plumage have been limited to trace occurrences in larks . Hymenobacter has been isolated from desert soil, airborne dust, glacial ice, irradiated meat, dung, and water . It is a rare component of vertebrate gut microbiotas  and the sole record in avian plumage was reported from black-plumaged swifts  that forage aerially and often fly nonstop for months .
The inverse correlations between the abundance of Deinococcus and Hymenobacter and ASV richness (Fig. 5) in vulture plumage microbiotas suggest that the combination of abiotic factors (UV radiation, heat, and desiccation) that suppress overall ASV richness, promotes the colonization, persistence, and growth of extremophiles. Agonistic interactions among ASVs waged through the production of bacteriocins [27, 70] or ASV competition for plumage-related resources such as preen-oil, dust, carrion smears, and feather keratin, may also contribute to the inverse relationships. We cannot directly address the underlying causes of the observed patterns because we lack data on the bacteriocins and nutritive constituents of plumage samples and the environmental conditions to which plumage was exposed in the months leading up to the collection date. The numerical dominance of Deinococcus and Hymenobacter in melanized vulture plumage has no known analog in the integumentary systems of terrestrial vertebrates. We note, however, that these bacterial genera are common constituents of biofilms recovered from solar panels [71, 72]. In any event, the unusual microbial rank-abundance profiles of vulture plumage appear to exhibit evidence of taxonomic filtering driven by solar irradiation.
All modern birds undergo periodic moult to replace worn, degraded, and damaged feathers [5, 7, 21] and the vast majority of species undergo at least one complete annual moult [73,74,75,76]. Bacterial degradation of feather keratin decreases the physical integrity of feathers, the functionality of plumage, and reduces the overall fitness of individuals. Moult reduces the load, at least temporarily, of bacteria and fungi that degrade plumage . Replacement of flight feathers (primaries and secondaries) in New World vultures requires as long as 2 years to complete [74, 77]. The prolonged staggered moult in vultures imposes strong selection on maintenance behaviours that preserve the aerodynamic integrity of flight feathers necessary for efficient soaring. We hypothesize that spread-wing sunning in the Cathartiformes [35, 47] evolved to maximize solar irradiation on ventral and dorsal surfaces of flight feathers—noting that dorsal and ventral contour plumage is also irradiated in the process. Frequent exposure to UV radiation, and associated desiccation and heat, is likely to kill, inhibit growth, or induce sporulation in mesophilic and UV-sensitive bacteria and fungi, including many that degrade keratin.
Spread-wing sunning in vultures is just one of many postural sunning variants observed in birds [29, 33]. The occurrence of stereotyped sunning in more than 50 avian families and 21 avian orders suggests an ancient origin dating to the early Palaeogene when most avian orders arose [45, 46]. This raises the question of when the arms race between keratinolytic bacteria and the bearers of pennaceous feathers arose, noting that well-feathered Archaeopteryx and maniraptoriform dinosaurs appeared in the fossil record in the late Jurassic . We forecast that intensive metagenomic surveys of plumage microbiotas of wild bird species will lead to novel insights on the evolution of pennaceous feathers in birds and dinosaurs.
Our results offers tangible support for the hypothesis  that the principal function of avian sunning behavior is the regulation of bacteria and other keratinolytic agents that degrade plumage quality. The prominence of UV-resistant extremophiles in generic rank-abundance curves suggests that solar irradiation may play a significant role in the assembly of vulture plumage microbiotas. Indeed, the high abundance of UV-resistant extremophiles in avian plumage may be a signature of taxonomic filtering via solar irradiation. In any event, bacterial rank-abundance profiles of melanized vulture plumage have no known analog in the integumentary systems of terrestrial vertebrates. Our results highlight the need for controlled in vivo experiments to test the effects of UV on microbial communities of avian plumage.
New World vultures (Cathartiformes)
The avian order Cathartiformes forms a monophylytic lineage composed of seven living species of obligate carrion scavengers that are now geographically restricted to the Americas [44, 47]. New World vultures are often confused with the ecologically convergent Old World vultures (Accipitriformes), which have been decimated by livestock pharmaceuticals ingested in carrion in Africa and Eurasia [79, 80]. The most recent common ancestor of the New World vultures and Old World vultures (Accipitriformes) lived 55–65 million year bp [44,45,46]. Five New World vulture species occur in Guyana in northern South America (Figs. 1,4): Sarcoramphus papa (body mass 3.1 kg, wingspan 167 cm), Coragyps atratus (body mass 1.4–1.9 kg, wingspan 123–136 cm), Cathartes aura  (body mass 1.2–1.5 kg, wingspan 153–159 cm), Cathartes melambrotos (body mass 1.3–1.8 kg, wingspan 153–172 cm), and Cathartes burrovianus (body mass 0.8–0.9 kg, wingspan 144–156 cm).
Vultures were sampled at two sites in southern Guyana in Upper Takuto-Upper Essequibo (Region 9), an area of exceptionally low human population density (0.42/km2).
Dadanawa Ranch (2° 49.28′ N, 59° 31.34′ W, 127 m above sea level). Thousands of free-range cattle graze a mosaic of unmanaged grassland and scrubby savannah on one of the largest ranches (> 4400 km2) in the Americas. Livestock mortality supports a large resident population of Coragyps atratus. More than 1500 vultures seasonly frequent a roost near the ranch headquarters. Smaller numbers of Cathartes burrovianus and Cathartes aura were observed daily. All three species were observed at baited net traps and carrion in multi-species feeding aggregations. Sarcoramphus papa is seasonally present but was not observed at the ranch headquarters during the brief sampling period. Cathartes melambrotos is unknown from open savannah grasslands on the ranch but occurs in adjacent blocks of forest.
Kanuku Mountains (3° 12.20′ N, 59° 24.20′ W, 109 m above sea level). The landcover consists of riverine forest and primary lower montane rainforest flanking sandbars of the Rupununi River. Little or no domestic livestock or agriculture occurs within 10 km of the study site. Cathartes melambrotos is common and outnumbers Coragyps atratus, Cathartes aura, and Sarcoramphus papa. Cathartes burrovianus is not known to occur in the forested interior of the Kanuku Mountains.
The collecting period (16–29 October 2015) was timed for the non-breeding season during the final weeks of the pronounced fall dry season in southern Guyana. No measurable rainfall was recorded during the 50 day period prior to collecting at Dadanawa. Daily low and high temperatures, respectively, averaged 25.6 C and 35.6 C during the collecting period. Vultures were netted or shot and necropsied soon after death (usually 15–30 min) by gloved preparators [81,82,83]. Clusters of contour plumage were plucked with sanitized forceps from the center of the mantle (dorsal feather tract) and center of the upper breast (ventral feather tract), placed in cryovials, and frozen immediately in liquid nitrogen. In contrast to the large flight feathers (primaries and secondaries) that are replaced in a staggered moult that may require up to 2 years to complete [74, 77], the smaller contour feathers of the mantle and breast are replaced annually . Contour plumage samples from Coragyps atratus and the three species of Cathartes were heavily melanized (black) whereas those from Sarcoramphus papa were pale beige (Figs. 1,4). All voucher specimens were deposited in the National Museum of Natural History (USNM), Smithsonian Institution (Additional file 1). A total of 34 vulture specimens were sampled: Sarcoramphus papa (n = 1), Coragyps atratus (n = 15), Cathartes aura (n = 3), Cathartes melambrotos (n = 10), and Cathartes burrovianus (n = 5).
DNA was extracted from plumage samples using the PowerSoil DNA Isolation Kit (Mo Bio Laboratories, Carlsbad, CA, USA) following manufacturer protocols. Extracted DNA was purified with the PowerClean Pro DNA Clean-Up Kit (Mo Bio Laboratories, Carlsbad, CA, USA) to remove PCR inhibitors. To check for potential contaminants from the extraction and clean-up kits, we processed an extraction control without any feather material alongside the plumage samples. We quantified DNA concentration with the 1X dsDNA High-Sensitivity Assay Kit on Qubit 4.0 (Thermo Fisher Scientific Waltham, MA, USA). For feather samples we used 2 μL of DNA extract. For the extraction control, we increased the volume of extract to 10 μL and observed no detectable contamination within the threshold of the assay (10 pg/μL). We amplified and sequenced the V4–V5 region of the 16S SSU rRNA gene (370 bp) from Bacteria and Archaea using primers 515FB and 926R , following procedures outlined in the Earth Microbiome Project 16S protocol [84, 85]. Detailed description of PCR conditions, library preparation, and sequencing on an Illumina MiSeq (Illumina, Inc., San Diego, CA, USA) are provided in Drovetski SV, O’Mahoney M, Ransome EJ, Matterson KO, Lim HC, Chesser RT and Graves GR . This protocol was modified to use AMPure XP beads (50 μL of library + 40 μL of beads; Beckman Coulter, Inc., Life Sciences Division, Indianapolis, IN, USA) to clean the final pooled library and the 500 or 600-cycle MiSeq Reagent Kits v2 or v3, (Illumina, Inc., San Diego, CA, USA), respectively, to accommodate the larger amplicon size. Each PCR reaction included a negative control.
Illumina data processing
Samples were distributed among three Illumina runs (two runs with MiSeq Reagent Kits v2 500 cycles and one with v3 600-cycles kit). Initially, raw FASTQ sequence reads for each Illumina run were demultiplexed in QIIME2 v. 2019.1 (qiime demux emp-paired) . Due to poor quality of the majority of reverse sequences, we processed only forward sequences (qiime demux emp-single). DADA2  was used for quality filtering and chimera removal (−-p-chimera-method consensus). All sequences were truncated to 250 bp. This fragment covers V4 region of the 16S SSU rRNA and was 2–4 bp shorter than amplicons produced with primers 515FB and 806R , which are frequently used in microbiota studies that follow the Earth Microbiome Project 16S protocol [84, 85]. We combined feature tables and representative sequences from all three runs. We assigned taxonomy using the SILVA -based 16S classifier for region V4 of the 16S SSU rRNA gene (SILVA 132 99% OTUs from 515F/806R region of sequences available at https://docs.qiime2.org/2019.1/data-resources/). Following taxonomic assignment, we filtered out sequences classified as Eukaryota, Mitochondria, Chloroplast, and those that were not classified to phylum. We aligned sequences and constructed phylogenetic trees using the align-to-tree-mafft-fasttree pipeline in QIIME2. The resulting feature table was exported into biom format and converted to the tab-delimited text in QIIME v. 1.9.1–20,150,604 . Taxonomy was downloaded in the comma-separated format, and the rooted tree was downloaded in the newik format. Taxonomy was added to the tab-delimited table in Excel 14.7.7 using VLOOKUP function and the table was converted back to the biom format in QIIME. Finally, we filtered out ASVs with abundance fractions < 0.001 in individual samples using script from https://gist.github.com/adamrp/7591573. Sequence abundance by ASV was cumulative sum scaled and log2-transformed ASV abundances in QIIME (normalize_table.py–a CSS). All statistical analyses were based on log2 CSS ASV normalized sequence counts.
We constructed rarefaction plots with R package Vegan 2.4–4  (Additional file 4) to ensure our sequencing coverage was sufficient. We used R version 3.3.3 (http://www.R-project.org) to generate linear regression, box, and column plots.
We used QIIME for α and β-diversity analyses. We calculated the number of observed ASVs and Shannon-Wiener diversity indices (alpha_diversity.py) to compare microbiota richness and evenness among vulture species and feather tracts. Microbial richness at different taxonomic levels was assessed in QIIME (summarize_taxa_through_plots.py or summarize_taxa.py). We calculated weighted UniFrac distances  among individual samples in QIIME (beta_diversity_through_plots.py). To compare microbiota composition between feather tracts and among vulture species, we generated principal coordinate analyses (PCoA) plots using the betadisper function in Vegan 2.4–4. Weighted UniFrac distances were used in the Permutational Multivariate Analysis of Variance (PERMANOVA)  implemented in Vegan 2.4–4  to test for the effect of species, plumage region (dorsal or ventral plumage tract), and their interaction on microbial variation among samples.
Availability of data and materials
The raw sequences archive has been deposited to GenBank: BioProject PRJNA521985: New World Vulture feather microbiota.
Amplicon sequence variant
Cumulative sum scaled
Principal coordinate analysis
Permutational multivariate analysis of variance
Smith VS, Ford T, Johnson KP, Johnson PCD, Yoshizawa K. Multiple lineages of lice pass through the K-Pg boundary. Biol Lett. 2011;7:782–5.
Marshall AG. The ecology of ectoparasitic insects. London: Academic Press Inc. (London) Ltd; 1981.
Xia Y, Massé DI, McAllister TA, Kong Y, Seviour R, Beaulieu C. Identity and diversity of archaeal communities during anaerobic co-digestion of chicken feathers and other animal wastes. Bioresour Technol. 2012;110:111–9.
Pugh GJF, Evans MD. Keratinophilic fungi associated with birds: I. Fungi isolated from feathers, nests and soils. Trans Br Mycol Soc. 1970;54:233–40.
Burtt EH, Ichida JM. Occurrence of feather-degrading bacilli in the plumage of birds. Auk. 1999;116:364–72.
Whitaker JM, Cristol DA, Forsyth MH. Prevalence and genetic diversity of Bacillus licheniformis in avian plumage. J Field Ornithol. 2005;76:264–70.
Stettenheim PR. The integumentary morphology of modern birds - an overview. Am Zool. 2000;40:461–77.
Bradbury JH. The structure and chemistry of keratin fibers. Adv Protein Chem. 1973;27:111–211. https://doi.org/10.1016/s0065-3233(08)60447-7.
Fraser RDB, MacRae TP, Parry DAD, Suzuki E. The structure of feather keratin. Polymer. 1971;12:35–56.
Noval JJ, Nickerson WJ. Decomposition of native keratin by Streptomycer fradiae. J Bacteriol. 1959;77:251–63.
Cheng S-W, Hu H-M, Shen S-W, Takagi H, Asano M, Tsai Y-C. Production and characterization of keratinase of a feather-degrading Bacillus licheniformis PWD-1. Biosci Biotechnol Biochem. 1995;59:2239–43.
Sangali S, Brandelli A. Feather keratin hydrolysis by a Vibrio sp strain kr2. J Appl Microbiol. 2000;89:735–43.
Williams CM, Richter CS, Mackenzie JM, Shih JCH. Isolation, identification, and characterization of a feather-degrading bacterium. Appl Environ Microbiol. 1990;56:1509–15.
Riffel A, Lucas F, Heeb P, Brandelli A. Characterization of a new keratinolytic bacterium that completely degrades native feather keratin. Arch Microbiol. 2003;179:258–65.
Gupta R, Ramnani P. Microbial keratinases and their prospective applications: an overview. Appl Microbiol Biotechnol. 2006;70:21.
Brandelli A. Bacterial keratinases: useful enzymes for bioprocessing agroindustrial wastes and beyond. Food Bioprocess Technol. 2008;1:105–16.
Shawkey MD, Mills KL, Dale C, Hill GE. Microbial diversity of wild bird feathers revealed through culture-based and culture-independent techniques. Microb Ecol. 2005;50:40–7.
Gunderson AR. Feather-degrading bacteria: a new frontier in avian and host-parasite research? Auk. 2008;125:972–9.
Dille JW, Rogers CM, Schneegurt MA. Isolation and characterization of bacteria from the feathers of wild dark-eyed juncos (Junco hyemalis). Auk. 2016;133:155–67.
Bisson IA, Marra PP, Burtt EH, Sikaroodi M, Gillevet PM. A molecular comparison of plumage and soil bacteria across biogeographic, ecological, and taxonomic scales. Microb Ecol. 2007;54:65–81.
Kent CM, Burtt EH. Feather-degrading bacilli in the plumage of wild birds: prevalence and relation to feather wear. Auk. 2016;133:583–92.
Musitelli F, Ambrosini R, Caffi M, Caprioli M, Rubolini D, Saino N, Franzetti A, Gandolfi I. Ecological features of feather microbiota in breeding common swifts. Ethol Ecol Evol. 2018;30:569–81.
Jacob S, Sallé L, Zinger L, Chaine AS, Ducamp C, Boutault L, Russell AF, Heeb P. Chemical regulation of body feather microbiota in a wild bird. Mol Ecol. 2018;27:1727–38.
van Veelen HPJ, Salles JF, Tieleman BI. Multi-level comparisons of cloacal, skin, feather and nest-associated microbiota suggest considerable influence of horizontal acquisition on the microbiota assembly of sympatric woodlarks and skylarks. Microbiome. 2017;5:156.
Dickinson EC, Christidis L. (eds). The Howard and Moore complete checklist of the birds of the world, 4th edition, Vol. 2: passerines. Eastbourne: Aves Press; 2014.
Dickinson EC, Remsen JV. (eds). The Howard and Moore complete checklist of the birds of the world, 4th edition, Vol. 1: non-passerines. Eastbourne: Aves Press; 2013.
Javůrková VG, Kreisinger J, Procházka P, Požgayová M, Ševčíková K, Brlík V, Adamík P, Heneberg P, Porkert J. Unveiled feather microcosm: feather microbiota of passerine birds is closely associated with host species identity and bacteriocin-producing bacteria. ISME J. 2019;13:2363–76.
Azcárate-García M, González-Braojos S, Díaz-Lora S, Ruiz-Rodríguez M, Martín-Vivaldi M, Martínez-Bueno M, Moreno J, Soler JJ. Interspecific variation in deterioration and degradability of avian feathers: the evolutionary role of microorganisms. J Avian Biol. 2020;51:e02320.
Simmons KEL. The sunning behaviour of birds: a guide for ornithologists. Bristol: Bristol Ornithological Club; 1986.
Clayton DH, Koop JAH, Harbison CW, Moyer BR, Bush SE. How birds combat ectoparasites. Open Ornithol J. 2010;3:41–71.
Ehrlich PR, Dobkin DS, Wheye D. The adaptive significance of anting. Auk. 1986;103:835.
Hauser DC. Some observations on sun-bathing in birds. Wilson Bull. 1957;69:78–90.
Kennedy RJ. Sunbathing behaviour of birds. British Birds. 1969;62:249–58.
Ohmart RD, Lasiewski RC. Roadrunners: energy conservation by hypothermia and absorption of sunlight. Science. 1971;172:67–9.
Clark RG, Ohmart RD. Spread-winged posture of Turkey vultures: single or multiple function? Condor. 1985;87:350–5.
Moyer BR, Wagenbach GE. Sunning by black Noddies (Anous minutus) may kill chewing lice (Quadraceps hopkinsi). Auk. 1995;112:1073–7.
Houston DC. A possible function of sunning behavior by griffon vultures, Gyps spp., and other large soaring birds. Ibis. 1980;122:366–9.
Potter EF, Hauser DC. Relationship of anting and sunbathing to molting in wild birds. Auk. 1974;91:537–63.
Saranathan V, Burtt EH. Sunlight on feathers inhibits feather-degrading bacteria. Wilson J Ornithol. 2007;119:239–45.
Javůrková VG, Enbody ED, Kreisinger J, Chmel K, Mrázek J, Karubian J. Plumage iridescence is associated with distinct feather microbiota in a tropical passerine. Sci Rep. 2019;9:12921.
Shawkey MD, Igic B, Rogalla S, Goldenberg J, Clusella-Trullas S, D’Alba L. Beyond colour: consistent variation in near infrared and solar reflectivity in sunbirds (Nectariniidae). Sci Nat. 2017;104:78.
Marder J. Body temperature regulation in the brown-necked raven (Corvus corax ruficollis)—II. Thermal changes in the plumage of ravens exposed to solar radiation. Comp Biochem Physiol A Physiol. 1973;45:431–40.
Ward JM, Blount JD, Ruxton GD, Houston DC. The adaptive significance of dark plumage for birds in desert environments. Ardea. 2002;90:311–23.
Johnson JA, Brown JW, Fuchs J, Mindell DP. Multi-locus phylogenetic inference among New World vultures (Aves: Cathartidae). Mol Phylogenet Evol. 2016;105:193–9.
Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, Ho SYW, Faircloth BC, Nabholz B, Howard JT, et al. Whole-genome analyses resolve early branches in the tree of life of modern birds. Science. 2014;346:1320–31.
Prum RO, Berv JS, Dornburg A, Field DJ, Townsend JP, Lemmon EM, Lemmon AR. A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing. Nature. 2015;526:569–73.
Houston DC. Cathartidae (New World vultures). In: del Hoyo J, Elliott A, Sargatal J, editors. Handbook of Birds of the World, New World vultures to guineafowl, vol. 2. Barcelona: Lynx Edicions; 1994. p. 24–41.
Buckley NJ. Black vulture (Coragyps atratus). In: Rodewald PG, editor. The birds of North America. Ithaca: Cornell Lab of Ornithology; 1999. https://birdsna.org/Species-Account/bna/species/blkvul.
Metcalf JL, Xu ZZ, Weiss S, Lax S, Van Treuren W, Hyde ER, Song SJ, Amir A, Larsen P, Sangwan N, et al. Microbial community assembly and metabolic function during mammalian corpse decomposition. Science. 2016;351:158–62.
Roggenbuck M, Schnell IB, Blom N, Baelum J, Bertelsen MF, Ponten TS, Sorensen SJ, Gilbert MTP, Graves GR, Hansen LH. The microbiome of New World vultures. Nat Commun. 2014;5:5498. https://doi.org/10.1038/ncomms6498.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Paulson JN, Stine OC, Bravo HC, Pop M. Robust methods for differential abundance analysis in marker gene surveys. Nat Methods. 2013;10:1200–2.
McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10:e1003531.
Battista JR. Against all odds: the survival strategies of Deinococcus radiodurans. Annu Rev Microbiol. 1997;51:203–24.
Daly MJ. A new perspective on radiation resistance based on Deinococcus radiodurans. Nat Rev Microbiol. 2009;7:237–45.
Battista JR, Rainey FA. Deinococcus. In: Boone DR, Castenholz RW, editors. Bergey's Manual of Systematic Bacteriology Second Edition Volume One The Archaea and the Deeply Branching and Phototrophic Bacteria. New York: Springer; 2001. p. 396–403.
Makarova KS, Aravind L, Wolf YI, Tatusov RL, Minton KW, Koonin EV, Daly MJ. Genome of the extremely radiation-resistant bacterium Deinococcus radiodurans viewed from the perspective of comparative genomics. Microbiol Mol Biol Rev. 2001;65:44.
Buczolits S, Busse HJ. Hymenobacter. In Bergey's Manual of Systematic Bacteriology Second Edition Volume Four The Bacteroidetes, Spirochaetes, Tenericutes (Mollicutes), Acidobacteria, Fibrobacteres, Fuscobacteria, Dictyoglomi, Gemmatimonadetes, Lentisphaerae, Verrucomicrobia, Chlamydiae, and Planctomycetes. Edited by Krieg NR, Staley JT, Brown DR, Hedlund BP, Paster BJ, Ward NL, Ludwig W, Whitman WB. New York: Springer; 2010: 397–404.
Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, Turner P, Parkhill J, Loman NJ, Walker AW. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87.
Waite D, Taylor M. Exploring the avian gut microbiota: current trends and future directions. Front Microbiol. 2015;6:673.
Hird SM, Sánchez C, Carstens BC, Brumfield RT. Comparative gut microbiota of 59 Neotropical bird species. Front Microbiol. 2015;6:1403.
Drovetski SV, O'Mahoney M, Ransome EJ, Matterson KO, Lim HC, Chesser RT, Graves GR. Spatial organization of the gastrointestinal microbiota of urban Canada geese. Sci Rep. 2018;8:3713.
Drovetski SV, O’Mahoney MJV, Matterson KO, Schmidt BK, Graves GR. Distinct microbiotas of anatomical gut regions display idiosyncratic seasonal variation in an avian folivore. Anim Microbiome. 2019;1:2.
Chao A. Non-parametric estimation of the number of classes in a population. Scand J Stat. 1984;11:265–70.
McLaughlin RW, Chen MM, Zheng JS, Zhao QZ, Wang D. Analysis of the bacterial diversity in the fecal material of the endangered Yangtze finless porpoise, Neophocaena phocaenoides asiaeorientalis. Mol Biol Rep. 2012;39:5669–76.
Kropáčková L, Pechmanová H, Vinkler M, Svobodová J, Velová H, Těšičký M, Martin J-F, Kreisinger J. Variation between the oral and faecal microbiota in a free-living passerine bird, the great tit (Parus major). PLoS One. 2017;12:e0179945.
Banskar S, Mourya DT, Shouche YS. Bacterial diversity indicates dietary overlap among bats of different feeding habits. Microbiol Res. 2016;182:99–108.
Jiang H-Y, Ma J-E, Li J, Zhang X-J, Li L-M, He N, Liu H-Y, Luo S-Y, Wu Z-J, Han R-C, Chen J-P. Diets alter the gut microbiome of crocodile lizards. Front Microbiol. 2017;8:2073.
Hedenström A, Norevik G, Warfvinge K, Andersson A, Bäckman J, Åkesson S. Annual 10-month aerial life phase in the common swift Apus apus. Curr Biol. 2016;26:3066–70.
Riley MA, Wertz JE. Bacteriocins: evolution, ecology, and application. Annu Rev Microbiol. 2002;56:117–37.
Porcar M, Louie KB, Kosina SM, Van Goethem MW, Bowen BP, Tanner K, Northen TR. Microbial ecology on solar panels in Berkeley, CA, United States. Front Microbiol. 2018;9:1–14.
Dorado-Morales P, Vilanova C, Peretó J, Codoñer FM, Ramón D, Porcar M. A highly diverse, desert-like microbial biocenosis on solar panels in a Mediterranean city. Sci Rep. 2016;6:29235.
Pyle P. Identification guide to north American birds. Part 1. Columbidae to Ploceidae. Bolinas: Slate Creek Press; 1997.
Pyle P. Identification guide to north American birds. Part II. Anatidae to Alcidae. Point Reyes Station: Slate Creek Press; 2008.
Jenni L, Winkler R. Moult and ageing of European passerines. London: Christopher Helm; 2011.
Gill FB, Prum RO. Ornithology. New York: W. H. Freeman, Macmillan Learning; 2019.
Chandler RM, Pyle P, Flannery ME, Long DJ, Howell SNG. Flight feather molt of Turkey vultures. Wilson J Ornithol. 2010;122:354–60.
Xu X, Zhou Z, Dudley R, Mackem S, Chuong C-M, Erickson GM, Varricchio DJ. An integrative approach to understanding bird origins. Science. 2014;346:1253293.
Shultz S, Baral HS, Charman S, Cunningham AA, Das D, Ghalsasi GR, Goudar MS, Green RE, Jones A, Nighot P, et al. Diclofenac poisoning is widespread in declining vulture populations across the Indian subcontinent. Proc R Soc Lond Ser B Biol Sci. 2004;271:S458–60.
Oaks JL, Gilbert M, Virani MZ, Watson RT, Meteyer CU, Rideout BA, Shivaprasad HL, Ahmed S, Iqbal Chaudhry MJ, Arshad M, et al. Diclofenac residues as the cause of vulture population decline in Pakistan. Nature. 2004;427:630–3.
Graves GR. Head color and caruncles of sympatric Cathartes vultures (Aves: Cathartidae) in Guyana and their possible function in intra- and interspecific signaling. Proc Biol Soc Wash. 2016;129:66–75.
Graves GR. Field measurements of gastrointestinal pH of New World vultures in Guyana. J Raptor Res. 2017;51:465–9.
Graves GR. Urohidrosis and tarsal color in Cathartes vultures (Aves: Cathartidae). Proc Biol Soc Wash. 2019;132:56–64.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.
Earth Microbiome Project: 16S Illumina amplicon protocol [http://press.igsb.anl.gov/earthmicrobiome/protocols-and-standards/16s/].
Drovetski SV, O’Mahoney M, Ransome EJ, Matterson KO, Lim HC, Chesser RT, Graves GR. Spatial organization of the gastrointestinal microbiota in urban Canada geese. Sci Rep. 2018;8:3713.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37:852–7.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
vegan: community ecology package. R package version 2.4–3. [https://CRAN.R-project.org/package=vegan].
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73:1576–85.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.
We thank Justin DeFreitas, the staff of the Dadanawa Ranch, and Jim Whatton for field assistance. Hillary Potter and International SOS Inc. provided medical support. All laboratory and computer work were conducted in and with the support of the Laboratories of Analytical Biology (LAB) facilities of the National Museum of Natural History, Smithsonian Institution. The Guyana Wildlife Division provided assistance with permits. Damian Fernandes of the Guyana Protected Areas Commission issued permission to work in the Kanuku Mountains. This is publication number 234 in the Smithsonian Biological Diversity of the Guianan Shield Publication Series.
Holt Thrasher, Kevin Kimberlin, Brad Palmer, and David Ford funded the Guyana field expedition. Additional financial support was provided by Mike and Carol Shealy (Basis Foundation) and the Jaideep and Rachel Khanna Foundation. GG acknowledges continuing support of the Wetmore Fund (Smithsonian Institution) and the Smoketree Trust.
Vulture specimens were collected under IACUC permit 2015–11 (National Museum of Natural History [NMNH], Smithsonian Institution) and Guyana EPA research permit number 101215 BR 023 in accordance with the approved regulations and protocols.
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.
Metadata. Laboratory identification numbers of feather samples, museum catalog numbers for voucher specimens, DNA concentration, number of filtered sequence reads, amplicon sequence variants (ASV), and number of microbial genera and families.
Amplicon sequence variants (ASV) identified in the plumage samples of New World vultures.
Log2 CSS abundance of bacterial genera identified in plumage samples of New World vultures.
Rarefaction curves for ASVs and microbial genera, families, and phyla.
About this article
Cite this article
Graves, G.R., Matterson, K.O., Milensky, C.M. et al. Does solar irradiation drive community assembly of vulture plumage microbiotas?. anim microbiome 2, 24 (2020). https://doi.org/10.1186/s42523-020-00043-7
- Keratin-degrading bacteria
- Melanized plumage
- New World vultures
- Plumage microbiota
- Solar irradiation
- Spread-wing sunning
- Ultraviolet (UV) light