- Open Access
Protective role of the vulture facial skin and gut microbiomes aid adaptation to scavenging
© The Author(s) 2018
- Received: 27 April 2018
- Accepted: 3 October 2018
- Published: 11 October 2018
Vultures have adapted the remarkable ability to feed on carcasses that may contain microorganisms that would be pathogenic to most other animals. The holobiont concept suggests that the genetic basis of such adaptation may not only lie within their genomes, but additionally in their associated microbes. To explore this, we generated shotgun DNA sequencing datasets of the facial skin and large intestine microbiomes of the black vulture (Coragyps atratus) and the turkey vulture (Cathartes aura). We characterized the functional potential and taxonomic diversity of their microbiomes, the potential pathogenic challenges confronted by vultures, and the microbial taxa and genes that could play a protective role on the facial skin and in the gut.
We found microbial taxa and genes involved in diseases, such as dermatitis and pneumonia (more abundant on the facial skin), and gas gangrene and food poisoning (more abundant in the gut). Interestingly, we found taxa and functions with potential for playing beneficial roles, such as antilisterial bacteria in the gut, and genes for the production of antiparasitics and insecticides on the facial skin. Based on the identified phages, we suggest that phages aid in the control and possibly elimination, as in phage therapy, of microbes reported as pathogenic to a variety of species. Interestingly, we identified Adineta vaga in the gut, an invertebrate that feeds on dead bacteria and protozoans, suggesting a defensive predatory mechanism. Finally, we suggest a colonization resistance role through biofilm formation played by Fusobacteria and Clostridia in the gut.
Our results highlight the importance of complementing genomic analyses with metagenomics in order to obtain a clearer understanding of the host-microbial alliance and show the importance of microbiome-mediated health protection for adaptation to extreme diets, such as scavenging.
- Colonization resistance
- Diet specialization
Vultures are composed of two clades of carrion-scavenging birds that diverged more than 60 million years ago , the New World vultures (Cathartidae) and Old World vultures (Accipitridae). Vultures global populations are under serious threats due to e.g. collisions with wind-energy turbines , their use for traditional medicine , or ingestion of lead bullets from deer carcasses . Vultures are known as “nature’s clean-up crew”, as they feed on tissues of animals that have died mainly from malnutrition, accidents, predation, and diseases [5–7]. Vultures are thus exposed to a variety of pathogens, including those that cause anthrax, tuberculosis, and brucellosis. A better understanding of various aspects for their biology are necessary, such as their susceptibility to the pathogens in their diet and their role in the transmission of infectious diseases [8, 9].
Vertebrate carcasses are very nutrient-rich resources. It has been speculated that the release of toxins and pathogenicity genes in the carcass microbiome are part of a microbial strategy for outcompeting other microbes [10, 11]. The main colonizers of a carcass are microbes originating from the microbiome of the animal when alive, some of which might become pathogenic in the carcass environment . Other components of the post-mortem microbiome include soil-dwelling bacteria, nematodes, fungi, and insects . In spite of the potentially serious health implications posed by their consumption, the pathogenic repertoire of the gut of these birds has not been fully characterized in relation to their possible implications in the environment. Thus, one of the most intriguing aspects of vulture biology is how they protect themselves against the health challenges posed by their dietary source. Physiologic, genetic and genomic analyses of different species of vultures have explored this aspect and identified genes associated with respiration, immunity, and gastric secretion as possible adaptations to its scavenging diet [14, 15]. For example, due to the previously reported very low stomach pH of a small sample of Old World vulture species , it has been suggested that the vulture stomach acidity serves as a filter of potential pathogens [14, 17].
With the genomic revolution, it has become apparent that besides genomic changes, host-associated microbiota plays an important role in diet specialization across vertebrates  and that the gut microbiome may play a highly relevant yet unexplored role in diet-driven speciation . The gut microbiome is intimately related to digestion functions, such as energy harvest, nutrient acquisition, and intestinal homeostasis . It has also been shown that the microbiome plays a health protective role to the host by interacting with the host’s immune system and mediating colonization-resistance against pathogens [21, 22]. Furthermore, disorders in the microbiome can lead to diseases such as irritable bowel syndrome, inflammatory bowel disease, obesity, and diabetes [23–25]. In light of the key roles that host-microbiome relationships can play in adaptation, it has been acknowledged that vulture genomic adaptations alone may not provide the full answer to the vulture adaptation to scavenging . However, neither the complete microbial taxonomic diversity (including non-bacterial microbes) nor the gene catalogue of the microbiome of any vulture species has been examined for their protective role against microbes that would normally pose serious health risks for other non-scavenging vertebrate species.
In order to evaluate the protective role of the vulture’s facial skin and gut microbiome, we generated metagenomic datasets from facial skin swabs and gut samples for two species of New World vultures, the black vulture (Coragyps atratus) and the turkey vulture (Cathartes aura), and performed taxonomic and functional metagenomic analyses.
Sampling method and DNA sequencing
We generated DNA shotgun metagenome datasets from a subset of the samples used by Roggenbuck et al. . Samples were collected over a period of several days in Tennessee, USA. Black vultures were live-trapped at deer carcasses and then transported to a central facility within a couple of hours of trapping. They were then euthanized with CO2, necropsied, and sampled within 30–45 min of death. Turkey vultures were shot at roosts, bagged individually, and transported to the processing facility where they were refrigerated 2–6 h before necropsy and sampling. To collect the gut samples, carcasses were opened to expose the entire gastrointestinal tract. A section of around 3–4 cm of the large intestine (hereafter called gut) located 2–3 cm above the cloaca was isolated with a pair of medical haemostats. Afterwards, 2–3 mL of sterile water was injected through the wall of the intestine with a sterile single-use syringe. The haemostat-blocked section of the intestine was gently massaged with the needle still inserted, and then the wash liquid was aspirated with the syringe. The aspirant was injected into a sterile vial containing RNAlater. Facial skin samples were taken by using sterile polyester swabs saturated with sterile water and wiped across the facial skin of the vultures. Swab tips were cut off and immersed in sterile vials filled with RNAlater. DNA was extracted and the shotgun libraries for HiSeq PE 100 were prepared using the Nextera library building kit following the manufacturer’s instructions, as in Roggenbuck et al. . From the total of 48 different sampled individuals (25 black vultures, and 23 turkey vultures), we used 33 facial skin samples (17 black vultures; 16 turkey vultures) and 47 intestinal samples (25 black vultures; 22 turkey vultures).
Two pipelines were used to process the raw reads. In the first approach, we removed adapter sequences and bases with quality < 15 using Trimmomatic v0.32 . Afterwards, in order to filter out non-bacterial reads derived from the vulture, human, and Phi phage (used as control required to increase library complexity on the Illumina HiSeq), the datasets were mapped against the bird genomes dataset of the avian phylogenomic project  (which includes the turkey vulture genome), the human (hg19), and the Phi phage genomes. Only the non-mapping reads were retained. The second approach was developed to take into account possible k-mer bias in the first bases of the reads that could have implications in the subsequent de novo assembly and gene prediction. To this end, we trimmed the first 16 bases of the reads with Trimmomatic v0.32. We then processed those reads with a Metagenomics Assembly and Gene Prediction Toolkit (MOCAT)  to clean them of low quality bases and adaptors and screen them versus the turkey vulture, human, and Phi phage genomes.
Relaxed filtering: In order to ensure the identification of low abundant taxa, we removed identifications with 90% of the abundance signal coming from only three samples.
Strict filtering: On top of the relaxed filtering, we removed the species with a coverage lower than the 1st quartile (Qu) value from the coverage distribution of the corresponding database.
With the filtered taxa, we classified the species as present in at least 90% and 50% of all the samples, thus defining a strict and a relaxed, respectively, taxonomic microbiome core.
For each database, we compared the taxa present only in the facial skin or gut datasets, those present in both, and those in significant differential abundance (P value < 0.05). To identify the differentially abundant species, we performed Wilcoxon and t-tests on the normalized abundance distribution of the identified species in the facial skin dataset versus the gut dataset. We also evaluated the taxonomic intra and inter sample variation between the facial skin and gut samples by calculating the Euclidean distances of their normalized abundances using the Ward.D method in R .
We then used the identifications of the relaxed filtering to test for microbial abundance correlations by calculating the Spearman correlation for each pairwise comparison of the microbes and calculated the P with a Bonferroni correction on those with a correlation value > 0.8 and < − 0.7. We also examined the enrichment and depletion of taxa within the facial skin and gut microbiome. To this end, we calculated their mean abundance (x̅) across the samples and compared them to the total distribution to calculate the Bonferroni corrected P. From these assignations, we also obtained a specific facial skin and gut core. As before, we defined two types of microbial cores: a strict one that retains those taxonomic identifications at the species level present in at least 80% of the samples of each sample type (facial skin or gut), and a relaxed one that retains those species present in at least 50% of the samples of each sample type.
We also identified the taxa of the top most abundant identified genes, which we defined as those genes with > 2000 mapping reads in the facial skin dataset and > 5000 in the gut dataset. We analysed the principal components (PCs) and rotation matrix of these taxa to identify the ones driving the variation within the facial skin and gut microbiomes. We defined “variation drivers” as those with an absolute rotation matrix value larger than the 3rd Qu value of the distributions from PC1, PC2 and PC3, and as “non-variation drivers” those with less than the 3rd Qu value of the distributions.
Besides the MGmapper identification, we used MOCAT as a complementary taxonomic identification method. For this approach, we used the taxonomic annotation given to the genes from the MOCAT strict non-redundant (NR) gene catalogue. This catalogue contains genes coding for proteins with a minimum length of 80 amino acids, not identified in low abundance, present only bacteria, fungi and virus, and that have an assigned Uniprot annotation (see “Methods”—“Functional profiling”). We analysed these identifications with MEGAN  having as input the search of the NR gene set catalogue against Uniprot using Ultra-Fast Sequence Search (USEARCH) .
We identified potential pathogens in the filtered bacterial and plasmid identifications. To this end, we obtained a list of the bacteria annotated with a disease from the database Pathosystems Resource Integration Center (PATRIC) . In PATRIC, bacteria are annotated as pathogenic if they have been reported with experimental data as causative of a disease. We further added the pathogenicity classification level of bacterial strain using the list from van Belkum , which was developed by the Commissie Genetische Modificatie (COGEM). Pathogenicity classes are defined as follows. Class 1 represents species that are commonly non-pathogenic, although there may be differences in virulence among the bacterial strains that should be taken into account. Class 2 contains species that can cause diseases in humans or animals but are unlikely to spread in the human population. Class 3 encompasses species that cause serious human diseases and can disseminate in the human population. We used the metadata of the pathogenic strains obtained from PATRIC to identify whether or not the identified bacteria are capable of sporulation and of antimicrobial resistance, together with the reported disease and host. For the identification of pathogenic plasmids, we used the list from Ho-Sui et al.  on the association of virulence factors with genomic islands of pathogenic bacteria.
We then used R v3.1.1  to examine the distribution of the total number of identified pathogenic bacteria. We grouped the samples by (i) vulture species (turkey and black vulture), and (ii) body sampling place (facial skin and gut). We then tested if the x̅ of the distributions were significantly different with a two-tailed and one-tailed (alternate greater) t-test. Next, we examined the number of samples in which each pathogenic bacterial strain, plasmid, resistance gene, and virulence factor was present. In order to get a potentially pathogenic core, we identified those taxa present in 50% (relaxed core) and 90% (strict core) of the samples. We also identified the taxa present only in the facial skin or the gut microbiome.
Abundance analyses of potential pathogens
In order to analyse the abundance of the potentially pathogenic microbes across the samples, we first rescaled the number of unique mapping reads by their percentage in the sample. We then removed taxa present in low abundance (those with 90% of their signal coming from < 4 samples). To determine if the identified potential pathogenic bacteria of facial skin and gut differed, we used the rescaled counts to build a dendrogram using a hierarchical clustering on the Euclidean distance. Afterwards, we examined which of the retained potentially pathogenic bacteria were present only on the facial skin or in the gut, and which ones were present in both. We then used a t-test to evaluate if the pathogenic bacterial abundance was statistically different by sample type and by vulture species.
16S bacterial taxonomic comparison
We compared the taxonomic bacterial identifications from both gut and facial datasets obtained with 16S analyses by Roggenbuck et al.  against the bacterial identifications from our metagenomics datasets using MGmapper, the taxonomic annotation of the de novo assembled genes with Uniprot, and the taxonomic identifications from the unmapped reads obtained with Double Index Alignment of Next-Generation Sequencing Data (DIAMOND) .
We performed de novo assembly with the De Bruijn Graph De Novo Assembler with Uneven Depth sequencing data (IDBA-UD)  and predicted genes with Prodigal . Afterwards, we generated a NR gene catalogue with USEARCH  by clustering the predicted genes with 90% identity and keeping the centroid sequences. We then searched the NR gene catalogue against Uniprot  with Ublast . The resulting identifications were functionally and taxonomically annotated with the use of a customized python script. Finally, we used DIAMOND v0.6.4  blastx to search the unmapped reads against Uniprot, keeping only the best hits for subsequent functional and taxonomic annotation.
For further functional assessment of the genes at the metabolic pathway level, we converted the Uniprot identifiers to the Kyoto Encyclopedia of Genes and Genomes (KEGG)  Enzyme Commission (EC) numbers and linked them to their corresponding metabolic pathway. Using the pathway classification information and the gene presence in each sample, we built a matrix for performing principal component analyses (PCA). Using the rotation matrix from the PCA we identified those pathways with an absolute rotation value within the minimum and 1st Qu values of the distributions of the PC1, PC2, and PC3. In order to distinguish the pathways driving most of the variation between the facial skin and gut microbiomes, we identified those pathways for which the absolute rotation value of their PC1, PC2, and PC3 was larger or equal to the 3rd Qu value of their corresponding distributions. We also obtained the Euclidean distances on the rescaled values of the matrices used for the PCAs.
As a second method, we used MOCAT with the Short Oligonucleotide Analysis Package for short-read de novo assembly (SOAPdenovo) v1.05  for assembling the reads cleaned with the approach that removed their first 16 bases. Subsequently, we corrected the assembly for indels and chimeric regions with SOAPdenovo. Using Prodigal, we then predicted the genes from all the samples, pooled them, and built an NR gene catalogue with Uclust  using a 90% identity threshold. Facial skin and gut datasets were treated separately. In the MGmapper core definition approach, the NR gene catalogue was obtained for each sample, then the catalogues were pooled and the unique genes were kept to compare their presence or absence across the samples. In contrast, in this approach using MOCAT, we built the cores based on the abundance of the reads mapping to the NR gene catalogue. To this end, we first mapped the reads of each sample against the NR gene catalogue and rescaled the counts values. Then, we removed those genes in low abundance (< 200 mapped reads), without a Uniprot annotation, and not derived from bacteria, archaea, virus, or fungi. We also removed genes coding for proteins with < 80 amino acids aligned to a hit from the Uniprot database. From these proteins, we also obtained a strict core (at least 80% of the samples) and a relaxed core (at least 50% of the samples). We also identified the top most abundant proteins (those with > 2000 mapping reads in the facial skin and > 5000 in the gut samples). On the relaxed functional core, we performed pathway functional analyses of their EC numbers with KEGG.
In order to search for antibiotic resistance genes, besides searching the ResFinder database with MGmapper, we used the Resfams v1.2 database  and the associated profile hidden Markov models. We searched the de novo assembled NR gene set of each sample against the Resfams profiles with a software for multiple alignment using hidden Markov models (HMMER v3.0) .
We produced a total of 342,279,763 raw read pairs from the facial skin samples and 512,803,778 from the gut samples. After cleaning and removing endogenous DNA by mapping against the bird genomes of the avian phylogenomic project , we obtained 79,938,910 read pairs from the facial skin samples (with a median of 1,378,000 read pairs per sample) and 144,877,366 from the gut samples (with a median of 1,118,000 per sample) (Additional file 1).
To prove the consistency of the taxonomic profiling between the two vulture species, we compared the number of identified microbial taxa in each species. We filtered the MGmapper  identifications of each whole-genome database by depth and breadth (percentage of covered reference sequence) of coverage and identified taxa differentially abundant in the facial skin and gut samples (Additional file 2). The number of identified bacteria was not significantly different between vulture species in the pooled datasets of facial skin and gut samples (P = 0.52, x̅black vulture = 366.97, x̅turkey vulture = 334.65). There were no significant differences between vulture species in the number of identified species of fungi (P = 0.43, x̅black vulture = 9, x̅turkey vulture = 7.5), viruses (P = 0.33, x̅black vulture = 21, x̅turkey vulture = 28.2), plasmids (P = 0.68, x̅black vulture = 186.85, x̅turkey vulture = 173.65), and protozoa (P = 0.21, x̅black vulture = 12, x̅turkey vulture = 9.62). Also, the number of identified proteins with resistance to antibiotics did not differ between vulture species (P = 0.64, x̅black vulture = 107.5, x̅turkey vulture = 100.78).
We compared our metagenomic bacterial identifications to those of Roggenbuck et al. . A total of 735 bacterial operational taxonomic units were identified analysing the 16S sequences, of which 93 were not found among our metagenomics identifications with strict filtering. When using the pre-filtering identifications from the whole-genome bacterial database and those identifications from the GreenGenes  and Silva  databases, only 14 genera were not identified in our analysis (Additional file 1).
Compared to the gut, the facial skin microbiome had higher microbial richness in terms of number of taxa and variation between individuals (Pprotozoa = 0.021, Pfungi = 0.029, Pbacteria = 0.0002). However, there was no significant difference in abundance (Pprotozoa = 0.514, Pfungi = 0.47, Pbacteria = 0.71). Although the number of identified virus was not significantly different between facial skin and gut (P = 0.58), viruses were statistically less abundant and variable in the facial skin than in the gut samples (P = 0.0002, Euclidean distancefacial skin = 9.04, Euclidean distancegut = 13.14). The most abundant bacterial genera in the facial skin microbiome were Pseudomonas, Bacteroides, and Prevotella, while the most abundant in the gut microbiome belonged to Escherichia, Campylobacter, and Clostridium (Additional file 3).
A total of 143 bacterial strains were significantly more abundant within the facial skin microbiome, 46 after the breadth filtering (mostly Pseudomonas). Within the gut microbiome, we identified 56 bacterial strains as the highest abundant, 33 after the breadth filtering (mostly Escherichia and Campylobacter). Bacterial strains most abundant in the facial skin dataset fell into three broad categories: (i) reported as potential human pathogens, (ii) associated with bioremediation (ionizing resistant, reducers of heavy metals, or oil degraders), and (iii) potentially beneficial (producers of antibiotics, insecticides and antifungals), usually intestinal bacteria, and related to water, plants, or soil. Those significantly more abundant in the gut dataset could be classified as: (i) reported as potential human pathogens, (ii) potentially beneficial, mostly intestinal or faecal bacteria from chicken, and (iii) fermenters and producers of intestinal metabolites.
Firmicutes was the most abundant phylum in the gut microbiome, followed by Proteobacteria, Fusobacteria in the third place, and Bacteroidetes in much less abundant in the fourth place. The most abundant class within Firmicutes was Clostridia (Fig. 2b). Among the Clostridiales, the most abundant families were Clostridiaceae, Peptostreptococcaceae and Lachnospiraceae (Fig. 2d). The most abundant taxa in the genus Clostridium were the potentially pathogenic C. perfringens and C. botulinum, followed by the beneficial C. carboxidovorans, C. sporogenes, and C. butyricum (Fig. 2e). The most abundant taxa from the Proteobacteria were Burkholderiales (from the Betaproteobacteria), Epsilonproteobacteria (from the delta/epsilon subdivision), and Enterobacteriales (mainly from Escherichia from the Gammaproteobacteria) (Fig. 2c). The most abundant taxa in the Fusobacteria were the potentially pathogenic Fusobacterium mortiferum, F. varium, and F. ulcerans (Fig. 2f).
Many Clostridia that are part of the normal human gut microbiome were also found in the vulture gut, likely playing roles in digestion. For instance, C. saccharolyticum was significantly more abundant in the gut. This bacterium, which is present in sewage sludge, ferments various carbohydrates into acetic acid, hydrogen, carbon dioxide, and ethanol , functions for which we identified related genes in the vulture gut microbiome. We also identified genes for cellulose degradation in the gut, along with the cellulose degraders C. cellulovorans and C. lentocellum , which were significantly more abundant in the gut microbiome. The gut microbiome also contained Bacteroides xylanisolvens, which breaks down xylan  and for which we identified a gene related to this function. Also, significantly more abundant in the gut microbiome than in the facial skin were the butanol producers C. beijerinckii  and C. saccharobutylicum .
We also identified protein coding genes involved in vitamin biosynthesis in the gut strict MOCAT NR functional core. For example, from the genera Hydrogenophaga, Herbaspirillum, and Gordonia, we identified the genes for d-threo-aldose 1-dehydrogenase, involved in ascorbate and aldarate metabolism, cobalamin biosynthesis, riboflavin biosynthesis, and vitamin B1 biosynthesis; and 2-ketopantoate reductase, involved in vitamin B5 production . Using the MGmapper gene core, which does not take taxa into account, we identified a larger abundance of genes in the gut than in the facial skin microbiome that belong to the metabolism of cofactors and vitamins (e.g. folate biosynthesis, vitamin B6 metabolism, riboflavin metabolism, and retinol metabolism). We also identified genes for the biosynthesis of various essential amino acids.
Comparison of the facial skin and gut microbiome variation
In the gut samples, we identified 604 species that showed significant variation in abundance among the samples and 348 that were relatively uniformly distributed. The variation drivers included 112 potentially pathogenic bacteria, such as species from the genera Listeria, Shigella, Yersinia, Bordetella, Shewanella, Erwinia, and Vibrio. Bacteria that were non-drivers included species from the genera Escherichia, Bacillus, Brucella, and Clostridium, among others. Non-variation driver phages included phages for Escherichia, Enterobacteria, and Shigella. Phages driving variation included phages for Clostridium, Yersinia, and Pseudomonas.
Of the 879 microbial species shared by the facial skin and gut samples (Fig. 3c), we identified 553 species as driving variation (62.9%), and 326 as non-variation drivers (37.1%). Among the most important variation drivers were Yersinia, Ralstonia, Rhizobium, Bifidobacterium, Bordetella, Listeria, and Burkholderia. The non-variation drivers included Brucella, Treponema, Clostridium, and Campylobacter. Looking at the phages, only Pseudomonas phages were variation drivers, while non-variation drivers included phages for Clostridium, Enterobacteria, Erwinia, and Shigella.
Functional potential characterization
Distances between and within the facial skin and gut samples
Facial skin vs facial skin
Gut vs facial skin
Gut vs gut
Amino acid metabolism
Biosynthesis of other secondary metabolites
Glycan biosynthesis and metabolism
Metabolism of cofactors and vitamins
Metabolism of other amino acids
Metabolism of terpenoids and polyketides
Xenobiotic biodegradation and metabolism
Core microbiome identification and attributes comparison
In the filtered MGmapper taxonomic profiling, we identified 1483 species in the facial skin samples, 638 of which occurred in the relaxed core, and only 184 species that occurred in the strict core. In the gut microbiome, we found 1419 microbial species, with 322 present in the relaxed core, and 129 in the strict core. In the functional characterization, we identified a total of 238,065 NR unique bacterial genes in the facial skin microbiome and 387,951 NR unique bacterial genes in the gut microbiome (Additional file 3).
We compared the taxonomic and functional composition of the facial skin and gut microbiomes and examined the microbial attributes of the taxa identified from the annotations of the assembled genes (Additional file 3). We found ~ 26× more habitat-specialized microbes in the facial skin than in the gut microbiome (facial skin = 8373, gut = 320). Consistent with the anaerobic gut environment, the gut microbiome had ~ 5× more anaerobic or microaerophilic bacteria than the facial microbiome (gut = 34,749, facial skin = 6699). The functional pathways clearly separated by sample type were the metabolism of other secondary metabolites, the glycan biosynthesis and metabolism, and the lipid metabolism, followed to a lesser extent by the metabolism of xenobiotics and the metabolism of other amino acids (Fig. 4). From the energy metabolism class, the methane metabolism was among the most abundant functions in the facial skin and gut microbiomes.
There was no significant difference (P = 0.44) in the number of identified potentially pathogenic plasmids between the facial skin (x̅ = 11.3) and gut (x̅ = 9.82) microbiomes. Likewise, there was no statistical difference in their abundance between the facial skin and gut samples (P = 0.78). Furthermore, no potentially pathogenic plasmid was present in 90% of the samples. Among those present in at least 50% of the samples were plasmids from Burkholderia vietnamiensis, Escherichia coli, and Ochrobactrum anthropic. Potentially pathogenic plasmids present only in the facial skin microbiome were from opportunistic pathogens such as Acinetobacter baumannii and Staphylococcus epidermidis, which is usually part of the normal skin microbiome . For example, the O. anthropi ATCC 49188 plasmid pOANT04 was present in 63 of the samples. O. anthropi is being increasingly recognized as a potentially problematic opportunistic and nosocomial human pathogen .
The Shiga toxin 1-converting phage BP-4795, which transmits virulence genes to its infected bacteria , was found in only five facial skin samples, and in 23 of the gut samples at various levels of abundance (max. coverage = 17.9%, max. mapping reads = 372). We also found the Shigella phage SfIV, which aids the virulence of Shigella flexneri , in 25 of the gut samples and 8 of the facial skin samples. We found that 75.2% of the identified pathogens are classified as level 2 pathogens in the facial skin samples (level 1 = 26 species, level 2 = 79 species), while 95.8% of those more abundant in the gut were classified as level 2 (level 1 = 2 species, level 2 = 46 species).
From the MGmapper strict core, the only protein identified in most of the facial skin samples (20 samples) and gut samples (36 samples) was an uncharacterized protein from Chlamydophila psittaci, an avian pathogen that causes avian chlamydiosis and epizootic outbreaks in mammals . Among the viruses from this functional strict core, we found the avian endogenous retrovirus EAV-HP and the avian leucosis virus. We also identified in higher abundance in the gut dataset Trichuris trichiura, causative of trichuriasis in humans [65, 66] (max. mapping readsfacial skin = 942, max. mapping readsgut = 565,870), and Eimeria brunetti, causative of haemorrhagic intestinal coccidiosis in poultry  (max. mapping readsfacial skin = 30, max. mapping readsgut = 1706). More abundant in the facial skin dataset we identified the fly Lucila cuprina (x̅gut = 11,910, x̅facial skin = 49,210), which causes sheep strike .
We identified only 49 bacteria in the gut dataset with potential for sporulation, including Fusobacterium necrophorum, Campylobacter jejuni, and L. monocytogenes (Additional file 4). Regarding the bacteria reported as zoonotic pathogens, we could only identify Streptococcus suis, a pathogen capable of transmission from swine to humans . Other identified bacteria with reported zoonotic capacity had very low abundance and were present in only one or two samples, so that they likely represent non-viable bacteria already dealt with by the vulture.
Microbiome composition and variability
Reduced core host microbiome
In order to differentiate the constant host microbiome from the one derived from variable and external influences (i.e. microbes derived from the carcass and the environment), we defined two types of microbiome cores. A relaxed core containing those elements (microbial species and genes) present in at least 50% of the samples, and a strict core with those present in at least 80% of the samples. We found that the relaxed core contained ~ 43% and ~ 22.7% of the facial skin and gut taxa, respectively, and the strict core only contained ~ 1% of the taxa in both gut and facial skin datasets. Notably, the distinction between carcass and established constant host-derived microbiome was complicated, even after the cores were defined. For example, the foodborne pathogen Salmonella enterica was present in the gut core (Fig. 6b). We additionally identified genes in the facial skin and gut microbiome functional strict cores that are related to putrescine, one of the main molecules produced in a carcass (Additional files 5 and 6). We found ~ 26× more habitat-specialized microbes in the facial skin than in the gut microbiome (facial skin = 8373, gut = 320), most likely due to the fact that a mammalian corpse is a disturbance habitat that selects for a specialized microbial community . Some of these community species likely derive from the carrion microbiota. For example, we identified phenol degrading bacteria in the vulture facial skin, such as Acinetobacter calcoaceticus  (max. mapping readsfacial skin = 478, max. mapping readsgut = 6), for which we identified its gene coding for phenol 2-monooxygenase in the NR gene catalogue. Phenolic compounds can act against foodborne pathogens and spoilage bacteria , suggesting that they derive from carrion dwellers, adapted to their competitive environment, instead of being part of the core vulture facial skin microbiome. Thus, we suggest that the vulture microbiome is a result of its scavenging diet, with part of the carcass microbiome leaving a profound footprint in the vulture microbiome.
Stomach acidity protection
The vulture stomach acidity has been suggested to prevent potential pathogens from reaching the gut [14, 17]. However, postprandial pH values observed in the stomach of black and turkey vultures appear to be no more acidic than those reported for domestic fowl and non-scavenging birds that consume large animal prey . Instrument readings of the pH are usually higher in a stomach with food contents, since the gastric acid is sparse and diluted by the water content in the lumpy food items. Sequential and independent probe values in the stomach can often provide different readings, even when the probe is used in the same location. This is reflected in the large standard deviations observed in the previously reported pH readings of black vultures (pH 3.8 ± 1.25) . The measurements were less acidic, and neutral in some occasions in the duodenum (pH 6.1 ± 0.48) and lower intestine (pH 6.0 ± 0.3) . Given that the carcass microbiome enters the vulture’s body mainly along with the ingested food items, the pH measurements suggest that the gastrointestinal acidity is not an efficient filter against all the potential pathogens present in a scavenging diet, rather it plays the general role of primary selection, which is not enough for all the potential pathogens in the carcass.
Top 10 potential disease-causing bacteria identified in the facial skin and gut
x̅ facial skin
Rank facial skin
Stenotrophomonas maltophilia K279a
Bacteraemia, bronchitis, pneumonia, urinary tract infection
Clostridium perfringens ATCC 13124
Plesiomonas shigelloides 302-73
Clostridium perfringens str. 13
Pseudomonas fluorescens A506
Clostridium perfringens SM101
Acinetobacter lwoffii SH145
Aeromonas salmonicida subsp. salmonicida A449
Acidovorax avenae subsp. avenae ATCC 19860
Bacterial leaf blight, brown stripe, red stripe
Propionibacterium propionicum F0230a
Campylobacter lari RM2100
Campylobacter jejuni subsp. doylei 269.97
Campylobacter jejuni RM1221
Fusobacterium nucleatum subsp. nucleatum ATCC 25586
Clostridium perfringens E str. JGS1987
Gastroenteritis, gas gangrene
Microbiome mediated protection
It has been shown that microbes provide protection to the host against pathogenic bacteria, thus we hypothesized that the vulture microbiome plays a protective role in terms of combating, preventing or maintaining in balance the abundance of potential pathogens. Accordingly, we identified functional and taxonomic protective elements that could be classified as related to (i) beneficial bacterial taxa and functions, (ii) phages, (iii) predatory eukaryotes, and (iv) colonization resistance (Fig. 7c, Additional file 5).
Beneficial bacterial taxa
Consistent with our microbiome-mediated protection hypothesis, we identified Hylemonella gracilis as part of the facial skin core, which has been shown to prevent long term colonization by Yersinia pestis . Other beneficial bacteria present in both gut and facial skin microbiomes include Lactobacillus sakei, an antilisterial bacterium . We also identified several genes for the biosynthesis of antibiotics such as carbapenem, tetracycline, macrolides, and ansamycins, as well as resistance genes towards them (Additional file 6). The identification of insecticide, fungicide, and antiparasite related taxa and genes in the facial skin microbiome suggests protective mechanisms against possible eukaryotic pathogens present in the carcass (Additional file 6). For example, Pseudomonas entomophila, which causes lethality in flies , and for which we identified a gene coding for an insecticidal toxin SepC/Tcc class in the facial skin NR gene set catalogue (Additional file 5). The production of antibiotics to outcompete for resources is known in soil microbiomes, and recently similar strategies have been reported in the human nasal microbiome from commensal bacteria against pathogens . Our results suggest that the vulture’s facial skin microbiome plays a similar defensive role. Regarding the gut microbiome, commensal Clostridia are known to play an important role in the production of butyrate that the colonocytes use . Notably, C. butyricum was among the most abundant Clostridia in the gut microbiome. Future studies would be needed to go beyond the description presented here to test for competitive exclusion among bacteria in the vulture gut microbiome [79, 80].
Beneficial bacterial functions
Besides containing potential pathogenic microbes, carcasses also contain toxic and carcinogenic compounds , which pose health risks to the vulture, particularly to its facial skin, which is in direct contact with such compounds. Among the bacteria identified in higher abundance in the facial microbiome was Arthrobacter phenanthrenivorans, which is able to degrade phenanthrene, a skin-irritating polycyclic aromatic hydrocarbon (PAH). PAH are xenobiotic pollutants with negative health-effects found to be emitted from animal carcass , and previously reported in high concentrations in other vulture species . Interestingly, the largest variation on the metabolism of xenobiotics biodegradation was in the facial skin dataset (Additional file 6), with PAHs degradation metabolism being the most abundant subclass from the xenobiotics degradation pathway in both facial skin and gut. These findings suggest a microbiome protective role for the vulture against such compounds. In regards to the gut microbiome, the second most abundant Fusobacteria in the gut from the NR gene set was the gut butyrate-producing F. varium , for which we also identified its gene formate C-acetyltransferase, which is involved in butanoate metabolism, as the most abundant gene in the gut dataset. Interestingly, besides the use of butyrate for the colonocytes, it has been shown that butyrate glycerides have antimicrobial activity against C. perfringens and Salmonella typhimurium .
Phage controlled pathogen abundance
Phages in the human gut microbiome have been shown to play a protective role and the increasing identification of antibiotic resistance genes in pathogenic bacteria has led to the proposition of using phages as alternative therapies . Given the identification of potential antibiotic resistance genes in the vulture facial skin and gut microbiomes (Additional file 5), we investigated the possible role of phages in eliminating or balancing the abundance of potential pathogens. In the facial skin microbiome, Clostridium phages positively correlated with C. perfringens and C. botulinum, whereas in the gut we observed enterobacteria phages correlating to Escherichia fergusonii (Additional file 7). Furthermore, in the taxonomic annotations of the gut functional core, we identified the Salmonella phage SPN3US (Additional file 7), which has shown effective inhibition of S. enterica . From the facial skin functional strict core, the most abundant virus was phage BPP-1 (Additional file 7), which infects pathogenic Bordetella bacteria . These findings show that the phage sets in both facial skin and gut microbiomes are related to the presence of the potential pathogens most abundant in the corresponding sample type. They also suggest that phages could represent an alternative defence mechanism for the control, and possibly elimination of potential pathogens, as in phage therapy  (Additional file 5).
Predatory defence mechanism
In spite of being important elements in the gut microbiome, gut microbial eukaryotes remain largely unexplored. Thus, we investigated whether vulture gut microbial eukaryotes played any protective role. We identified the invertebrate Adineta vaga, which feeds on dead bacteria and protozoans, to be ~ 6.8× more abundant in the gut core than in the facial skin core. This identification suggests that a predatory mechanism may be exploited for defence in the vulture’s gut.
Biofilm formation and colonization resistance
Biofilms are assemblages of microbes associated within a matrix composed of extracellular polymeric substances that facilitate their adhesion to the surface, protection against antimicrobials, and better nutrient acquisition. The abundance of Fusobacteria in the gut has been suggested to play a particularly relevant role in lumen biofilm formation in the gastrointestinal tract . To explore this hypothesis, we searched for proteins related to biofilm formation (Additional file 7). In the core of the gut functional potential we identified biofilm related proteins from F. mortiferum, such as rubrerythrin, as well as from C. perfringens, such as UDP-glucuronic acid epimerase. Given that bacteria form biofilms in which they can thrive under different patterns of gene expression , this suggests that the identified potential pathogenic Clostridia and Fusobacteria from the gut microbiome might not pose pathogenic risks and instead offer colonization resistance against other external pathogens (Additional file 5). To further explore the possible colonization resistance role of Clostridia and Fusobacteria, we examined the gut functional core for toxins with potential effects on the vulture. We identified only few potentially pathogenic toxin coding genes from Fusobacterium (Additional file 7). Considering that F. varium has been shown to affect its human host in a beneficial manner by antagonizing colonization by pathogenic agents , we suggest that an important role of the gut Fusobacteria could be the formation of biofilms and colonization resistance, without representing a serious pathogenic threat. We identified the pathogenicity genes perfringolysin O and phospholipase C in the gut microbiome from C. perfringens. However, we also identified genes for the biosynthesis of short chain fatty acids from Clostridia, which can provide protection against inflammatory responses . Thus, we could classify the observed Clostridia into two types, (i) the potentially pathogenic, mainly represented by C. perfringens and (ii) the non-pathogenic, which may contribute to biofilm formation and health defence (Additional file 4).
Our findings strongly suggest that the turkey and black vultures have adapted to their scavenging diet with the help of their facial skin and gut microbiomes. Surprisingly, most of their microbiome consists of a large variable pool of environmental and carcass-derived microbiota, with only a small set of constant inhabitants. In particular, the presence of a wide variety of microbes reported as pathogenic to non-scavengers (mammals and other birds) without an apparent or a reported pathogenic effect on the vultures calls for deeper study. Further studies would be required to determine whether the microbes reported here are pathogenic to the wild vultures or if they serve as reservoirs, and to determine what is their zoonotic potential. A better characterization of wild-life as potential pathogenic reservoirs (with microbiological, epidemiological and surveillance data) would allow for better informed wild-life protection programs, particularly for those species in endangered status, such as some species of vultures (e.g. the white-backed vulture ). We highlight the identification within the vultures’ facial skin and gut microbiomes of defence mechanisms that are alternative to the use of antimicrobials, such as the use of predatory microbes, and the protective nature of colonization resistance through biofilm formation by Fusobacteria and Clostridia. However, further microbiology studies would be needed to isolate the relevant microbes and validate the antimicrobial mechanisms reported here from the vultures’ microbiomes.
The establishment of these suggested protective mechanisms in the vulture microbiome unveiled by metagenomics analyses highlights the important role that vultures play in their ecosystem. This role is the essential but underrated service of cleaning up carcasses that otherwise would spread microbial elements pathogenic to species without a specialized microbiome like that of the vultures. In conclusion, our results show the importance of complementing genomic analyses with metagenomics on the host microbiome in order to obtain a clearer understanding of the host-microbial alliance that aids the evolution of extreme dietary adaptations.
MLZM performed the metagenomics analyses. KMV and MLZM performed the search for the resistance genes and potentially pathogenic microbes. MLZM, LHH, MTPG, MR, SR, and TSP interpreted the results. MLZM wrote the manuscript. LHH, MTPG, MR, SR, and TSP provided significant input on the manuscript drafting. All authors read and approved the final manuscript.
We thank the Danish National High-Throughput DNA Sequencing Centre for the generation of the sequencing data. We are also grateful to George Pacheco and Lillian Anne Petersen for DNA extractions and library preparations. We also gratefully acknowledge the Danish National Supercomputer for Life Sciences—Computerome (computerome.dtu.dk) for the computational resources to perform the sequence analyses. We thank Prof. Gary Graves (Smithsonian Institution) for providing the vulture facial skin and gut samples and valuable input on the manuscript drafting. We also thank Brian Schmidt and Christina Gebhard (both Smithsonian Institution) for necropsying vultures; Blaine Hyle, Talon Redding, William Simmons and J. D. Freye (all USDA) for collecting vultures; and Keith Wehner, Blaine Hyle and Brett Dunlap (all USDA) for providing critical logistic support in Nashville. The Alexander Wetmore Fund of the Smithsonian Institution provided funding for fieldwork. Søren Brunak thanks the Novo Nordisk Foundation (Grant NNF14CC0001).
The authors declare that they have no competing interests.
Availability of data and materials
The data reported in this paper are tabulated in the additional files and the sequencing datasets are available in the NCBI SRA database with the id SRP153398 under to the BioProject PRJNA479518.
Consent for publication
Ethics approval and consent to participate
Black and turkey vultures were captured and euthanized during scheduled population control activities conducted near Nashville, Tennessee (USA) on the 13th–17th February 2012 by the United States Department of Agriculture (USDA), Animal and Plant Health Inspection Service (APHIS), under the authority of the US Fish & Wildlife Service.
MLZM and MTPG thank Lundbeck Foundation Grant R52-A5062 for funding their research.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, et al. Whole genome analyses resolve early branches in the tree of life of modern birds. Science. 2014;346:1320–31.View ArticleGoogle Scholar
- de Lucas M, Ferrer M, Bechard MJ, Muñoz AR. Griffon vulture mortality at wind farms in southern Spain: distribution of fatalities and active mitigation measures. Biol Conserv. 2012;147:184–9.View ArticleGoogle Scholar
- Mckean S, Mander M, Diederichs N, Ntuli L, Mavundla K, Williams V, et al. The impact of traditional use on vultures in South Africa. Vulture News. 2013;65:15–36.View ArticleGoogle Scholar
- Snyder NFR, Snyder H. The California condor: a saga of natural history and conservation. Cambridge: Academic Press; 2000. p. 410.Google Scholar
- Houston DC. The ecology of Serengeti vultures. Oxford: Oxford Press; 1972.Google Scholar
- Houston DC. The role of griffon vultures Gyps africanus and Gyps ruppellii as scavengers. J Zool. 1974;172:35–46.View ArticleGoogle Scholar
- Grigg NP, Krilow JM, Gutierrez-Ibanez C, Wylie DR, Graves GR, Iwaniuk AN. Anatomical evidence for scent guided foraging in the turkey vulture. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-17794-0.View ArticlePubMedPubMed CentralGoogle Scholar
- Koenig R. Ornithology. Vulture research soars as the scavengers’ numbers decline. Science. 2006. https://doi.org/10.1126/science.312.5780.1591.View ArticlePubMedGoogle Scholar
- Margalida A, Colomer MA. Modelling the effects of sanitary policies on European vulture conservation. Sci Rep. 2012;2:753.View ArticleGoogle Scholar
- Anzen DH. Herbivores and the number of tree species in tropical forests. Am Nat. 1970;104(940):501–28.View ArticleGoogle Scholar
- Vollaard EJ, Clasener HA. Colonization resistance. Antimicrob Agents Chemother. 1994;38:409–14.View ArticleGoogle Scholar
- Loeffler AG, Hart MN. Chapter 26—Infectious diseases. In: Introduction to human disease: pathophysiology for health professionals. 6th ed. Burlington: Jones & Bartlett Learning; 2014. p. 396–7.Google Scholar
- Metcalf JL, Xu ZZ, Weiss S, Lax S, Van Treuren W, Hyde ER, et al. Microbial community assembly and metabolic function during mammalian corpse decomposition. Science. 2015;351:158–62.View ArticleGoogle Scholar
- Chung O, Jin S, Cho YS, Lim J, Kim H, Jho S, et al. The first whole genome and transcriptome of the cinereous vulture reveals adaptation in the gastric and immune defense systems and possible convergent evolution between the Old and New World vultures. Genome Biol. 2015;16:215.View ArticleGoogle Scholar
- Mateos-Hernndez L, Crespo E, Fuente JDL, de la Lastra JMP. Identification of key molecules involved in the protection of vultures against pathogens and toxins. In: Baptista GR, editor. An integrated view of the molecular recognition and toxinology—from analytical procedures to biomedical applications. Rijeka: InTech; 2013.Google Scholar
- Houston DC, Cooper JE. The digestive tract of the whiteback griffon vulture and its role in disease transmission among wild ungulates. J Wildl Dis. 1975;11:306–13.View ArticleGoogle Scholar
- Roggenbuck M, Bærholm Schnell I, Blom N, Bælum J, Bertelsen MF, Pontén TS, et al. The microbiome of New World vultures. Nat Commun. 2014. https://doi.org/10.1038/ncomms6498.View ArticlePubMedGoogle Scholar
- Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI. Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol. 2008;6:776–88.View ArticleGoogle Scholar
- Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320:1647–51.View ArticleGoogle Scholar
- Semova I, Carten JD, Stombaugh J, Mackey LC, Knight R, Farber SA, et al. Microbiota regulate intestinal absorption and metabolism of fatty acids in the zebrafish. Cell Host Microbe. 2012;12:277–88.View ArticleGoogle Scholar
- Hooper LV, Littman DR, Macpherson AJ. Interactions between the microbiota and the immune system. Science. 2012;336:1268–73.View ArticleGoogle Scholar
- Buffie CG, Pamer EG. Microbiota-mediated colonization resistance against intestinal pathogens. Nat Rev Immunol. 2013;13:790–801.View ArticleGoogle Scholar
- Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490:55–60.View ArticleGoogle Scholar
- Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444:1022–3.View ArticleGoogle Scholar
- Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermúdez-Humarán LG, Gratadoux J-J, et al. Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci. 2008;105:16731–6.View ArticleGoogle Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.View ArticleGoogle Scholar
- Zhang G, Li B, Li C, Gilbert M, Jarvis E, Consortium TAG, et al. Comparative genomic data of the avian phylogenomics project. GigaSci Database. 2014;3:26.View ArticleGoogle Scholar
- Kultima JR, Sunagawa S, Li J, Chen W, Chen H, Mende DR, et al. MOCAT: a metagenomics assembly and gene prediction toolkit. PLoS ONE. 2012;7:1–6.View ArticleGoogle Scholar
- Petersen TN, Lukjancenko O, Thomsen MCF, Maddalena Sperotto M, Lund O, Møller Aarestrup F, et al. MGmapper: reference based mapping and taxonomy annotation of metagenomics sequence reads. PLoS ONE. 2017;12:e0176469.View ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticleGoogle Scholar
- Nielsen HB, Almeida M, Juncker AS, Rasmussen S, Li J, Sunagawa S, et al. Identification and assembly of genomes and genetic elements in complex metagenomic samples without using reference genomes. Nat Biotechnol. 2014;32:822–8.View ArticleGoogle Scholar
- Peterson J, Garges S, Giovanni M, McInnes P, Wang L, Schloss JA, et al. The NIH human microbiome project. Genome Res. 2009;19:2317–23.View ArticleGoogle Scholar
- Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. J Antimicrob Chemother. 2012;67:2640–4.View ArticleGoogle Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.View ArticleGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.View ArticleGoogle Scholar
- R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.Google Scholar
- Huson DH, Weber N. Microbial community analysis using MEGAN. Methods Enzymol. 2013;531:465–85.View ArticleGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticleGoogle Scholar
- Gillespie JJ, Wattam AR, Cammer SA, Gabbard JL, Shukla MP, Dalay O, et al. PATRIC: the comprehensive bacterial bioinformatics resource with a focus on human pathogenic species. Infect Immun. 2011;79:4286–98.View ArticleGoogle Scholar
- Van Belkum A. Classification of bacterial pathogens. COGEM Res Rep CGM. 2011;2011-07:1–113.Google Scholar
- Ho Sui SJ, Fedynak A, Hsiao WWL, Langille MGI, Brinkman FSL. The association of virulence factors with genomic islands. PLoS ONE. 2009;4:e8094.View ArticleGoogle Scholar
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2014;12:59–60.View ArticleGoogle Scholar
- Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD: a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28:1420–8.View ArticleGoogle Scholar
- Hyatt D, Chen G-L, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:119.View ArticleGoogle Scholar
- UniProt Consortium. The universal protein resource (UniProt). Nucleic Acids Res. 2008;36:D190–5.View ArticleGoogle Scholar
- Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.View ArticleGoogle Scholar
- Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012. https://doi.org/10.1186/2047-217x-1-18.View ArticlePubMedPubMed CentralGoogle Scholar
- Gibson MK, Forsberg KJ, Dantas G. Improved annotation of antibiotic resistance determinants reveals microbial resistomes cluster by ecology. ISME J. 2015;9:207–16.View ArticleGoogle Scholar
- Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7:e1002195.View ArticleGoogle Scholar
- Martínez I, Lattimer JM, Hubach KL, Case JA, Yang J, Weber CG, et al. Gut microbiome composition is linked to whole grain-induced immunological improvements. ISME J. 2013;7:269–80.View ArticleGoogle Scholar
- Rey FE, Faith JJ, Bain J, Muehlbauer MJ, Stevens RD, Newgard CB, et al. Dissecting the in vivo metabolic potential of two human gut acetogens. J Biol Chem. 2010;285:22082–90.View ArticleGoogle Scholar
- Leisner JJ, Laursen BG, Prévost H, Drider D, Dalgaard P. Carnobacterium: positive and negative effects in the environment and in foods. FEMS Microbiol Rev. 2007;31:592–613.View ArticleGoogle Scholar
- Suzuki M, Yamamoto T, Kawai Y, Inoue N, Yamazaki K. Mode of action of piscicocin CS526 produced by Carnobacterium piscicola CS526. J Appl Microbiol. 2005;98:1146–51.View ArticleGoogle Scholar
- Mirray WD, Khan AW, van den Berg L. Clostridium saccharolyticum sp. nov., a saccharolytic species from sewage sludge. Int J Syst Bacteriol. 1982;32:132–5.View ArticleGoogle Scholar
- He YL, Ding YF, Long YQ. Two cellulolytic Clostridium species: Clostridium cellulosi sp. nov. and Clostridium cellulofermentans sp. nov. Int J Syst Bacteriol. 1991;41:306–9.View ArticleGoogle Scholar
- Chassard C, Delmas E, Lawson PA, Bernalier-Donadille A. Bacteroides xylanisolvens sp. nov., a xylan-degrading bacterium isolated from human faeces. Int J Syst Evol Microbiol. 2008;58:1008–13.View ArticleGoogle Scholar
- George HA, Johnson JL, Moore WE, Holdeman LV, Chen JS. Acetone, isopropanol, and butanol production by Clostridium beijerinckii (syn. Clostridium butylicum) and Clostridium aurantibutyricum. Appl Environ Microbiol. 1983;45:1160–3.PubMedPubMed CentralGoogle Scholar
- Madihah MS, Ariff AB, Khalil MS, Suraini AA, Karim MI. Anaerobic fermentation of gelatinized sago starch-derived sugars to acetone-1-butanol-ethanol solvent by Clostridium acetobutylicum. Folia Microbiol. 2001;46:197–204.View ArticleGoogle Scholar
- Frodyma ME, Downs D. The panE gene, encoding ketopantoate reductase, maps at 10 minutes and is allelic to apbA in Salmonella typhimurium. J Bacteriol. 1998;180:4757–9.PubMedPubMed CentralGoogle Scholar
- Grice EA, Kong HH, Renaud G, Young AC, Bouffard GG, Blakesley RW, et al. A diversity profile of the human skin microbiota. Genome Res. 2008;18:1043–50.View ArticleGoogle Scholar
- Chain PSG, Lang DM, Comerci DJ, Malfatti SA, Vergez LM, Shin M, et al. Genome of Ochrobactrum anthropi ATCC 49188T, a versatile opportunistic pathogen and symbiont of several eukaryotic hosts. J Bacteriol. 2011;193:4274–5.View ArticleGoogle Scholar
- Creuzburg K, Recktenwald J, Kuhle V, Herold S, Hensel M, Schmidt H. The Shiga toxin 1-converting bacteriophage BP-4795 encodes an NleA-like type III effector protein. J Bacteriol. 2005;187:8494–8.View ArticleGoogle Scholar
- Jakhetia R, Talukder KA, Verma NK. Isolation, characterization and comparative genomics of bacteriophage SfIV: a novel serotype converting phage from Shigella flexneri. BMC Genomics. 2013;14:677.View ArticleGoogle Scholar
- Harkinezhad T, Geens T, Vanrompay D. Chlamydophila psittaci infections in birds: a review with emphasis on zoonotic consequences. Vet Microbiol. 2009;135:68–77.View ArticleGoogle Scholar
- Klementowicz JE, Travis MA, Grencis RK. Trichuris muris: a model of gastrointestinal parasite infection. Semin Immunopathol. 2012;34:815–28.View ArticleGoogle Scholar
- Jung RC, Beaver PC. Clinical observations on Trichocephalus trichiurus (whipworm) infestation in children. Pediatrics. 1951;8:548–57.PubMedGoogle Scholar
- Shirley MW, Smith AL, Tomley FM. The biology of avian Eimeria with an emphasis on their control by vaccination. Adv Parasitol. 2005;60:285–330.View ArticleGoogle Scholar
- East IJ, Eisemann CH. Vaccination against Lucilia cuprina: the causative agent of sheep blowfly strike. Immunol Cell Biol. 1993;71:453–62.View ArticleGoogle Scholar
- Wertheim HFL, Nghia HDT, Taylor W, Schultsz C. Streptococcus suis: an emerging human pathogen. Clin Infect Dis. 2009;48:617–25.View ArticleGoogle Scholar
- Paller G, Hommel RK, Kleber HP. Phenol degradation by Acinetobacter calcoaceticus NCIB 8250. J Basic Microbiol. 1995;35:325–35.View ArticleGoogle Scholar
- Nazzaro F, Fratianni F, De Martino L, Coppola R, De Feo V. Effect of essential oils on pathogenic bacteria. Pharmaceuticals. 2013;6:1451–74.View ArticleGoogle Scholar
- Graves G. Field measurements of gastrointestinal pH of New World vultures in Guyana. J Raptor Res. 2017;51:465–9.View ArticleGoogle Scholar
- Hou PC. The pathology of Clonorchis sinensis infestation of the liver. J Pathol Bacteriol. 1955;70:53–64.View ArticleGoogle Scholar
- Pawlowski DR, Raslawsky A, Siebert G, Metzger DJ. Identification of Hylemonella gracilis as an antagonist of Yersinia pestis persistence. J Bioterror Biodef. 2011. https://doi.org/10.4172/2157-2526.s3-004.View ArticleGoogle Scholar
- Bredholt S, Nesbakken T, Holck A. Industrial application of an antilisterial strain of Lactobacillus sakei as a protective culture and its effect on the sensory acceptability of cooked, sliced, vacuum-packaged meats. Int J Food Microbiol. 2001;66:191–6.View ArticleGoogle Scholar
- Vodovar N, Vinals M, Liehl P, Basset A, Degrouard J, Spellman P, et al. Drosophila host defense after oral infection by an entomopathogenic Pseudomonas species. Proc Natl Acad Sci. 2005;102:11414–9.View ArticleGoogle Scholar
- Zipperer A, Konnerth MC, Laux C, Berscheid A, Janek D, Weidenmaier C, et al. Human commensals producing a novel antibiotic impair pathogen colonization. Nature. 2016;535:511–6.View ArticleGoogle Scholar
- Pryde SE, Duncan SH, Hold GL, Stewart CS, Flint HJ. The microbiology of butyrate formation in the human colon. FEMS Microbiol Lett. 2002;217:133–9.View ArticleGoogle Scholar
- Freilich S, Zarecki R, Eilam O, Segal ES, Henry CS, Kupiec M, et al. Competitive and cooperative metabolic interactions in bacterial communities. Nat Commun. 2011. https://doi.org/10.1038/ncomms1597.View ArticlePubMedGoogle Scholar
- Louca S, Doebeli M. Transient dynamics of competitive exclusion in microbial communities. Environ Microbiol. 2016;18:1863–74.View ArticleGoogle Scholar
- Phan NT, Kim KH, Jeon EC, Kim UH, Sohn JR, Pandey SK. Analysis of volatile organic compounds released during food decaying processes. Environ Monit Assess. 2012;184:1683–92.View ArticleGoogle Scholar
- Chen SJ, Hsieh LT, Chiu SC. Emission of polycyclic aromatic hydrocarbons from animal carcass incinerators. Sci Total Environ. 2003;313:61–76.View ArticleGoogle Scholar
- Dhananjayan V, Muralidharan S. Levels of polycyclic aromatic hydrocarbons, polychlorinated biphenyls, and organochlorine pesticides in various tissues of white-backed vulture in India. Biomed Res Int. 2013. https://doi.org/10.1155/2013/190353.View ArticlePubMedPubMed CentralGoogle Scholar
- Dzink JL, Socransky SS. Amino acid utilization by Fusobacterium nucleatum grown in a chemically defined medium. Oral Microbiol Immunol. 1990;5:172–4.View ArticleGoogle Scholar
- Namkung H, Yu H, Gong J, Leeson S. Antimicrobial activity of butyrate glycerides toward Salmonella typhimurium and Clostridium perfringens. Poult Sci. 2011;90:2217–22.View ArticleGoogle Scholar
- Sulakvelidze A, Alavidze Z, Morris JG. Bacteriophage therapy. Antimicrob Agents Chemother. 2001;45:649–59.View ArticleGoogle Scholar
- Kagawa H, Ono N, Enomoto M, Komeda Y. Bacteriophage chi sensitivity and motility of Escherichia coli K-12 and Salmonella typhimurium Fla-mutants possessing the hook structure. J Bacteriol. 1984;157:649–54.PubMedPubMed CentralGoogle Scholar
- Liu M, Deora R, Doulatov SR, Gingery M, Eiserling FA, Preston A, et al. Reverse transcriptase-mediated tropism switching in Bordetella bacteriophage. Science. 2002;295:2091–4.View ArticleGoogle Scholar
- Focà A, Liberto MC, Quirino A, Marascio N, Zicca E, Pavia G, et al. Gut inflammation and immunity: what is the role of the human gut virome? Mediat Inflamm. 2015. https://doi.org/10.1155/2015/326032.View ArticleGoogle Scholar
- Keenan SW, Engel AS, Elsey RM. The alligator gut microbiome and implications for archosaur symbioses. Sci Rep. 2013;3:2877.View ArticleGoogle Scholar
- Prigent-Combaret C, Vidal O, Dorel C, Lejeune P. Abiotic surface sensing and biofilm-dependent regulation of gene expression in Escherichia coli. J Bacteriol. 1999;181:5993–6002.PubMedPubMed CentralGoogle Scholar
- Ushijima T, Ozaki Y. Potent antagonism of Escherichia coli, Bacteroides ovatus, Fusobacterium varium, and Enterococcus faecalis, alone or in combination, for enteropathogens in anaerobic continuous flow cultures. J Med Microbiol. 1986;22:157–63.View ArticleGoogle Scholar
- Vinolo MAR, Rodrigues HG, Nachbar RT, Curi R. Regulation of inflammation by short chain fatty acids. Nutrients. 2011;3:858–76.View ArticleGoogle Scholar
- Birdlife International. Gyps africanus (amended version of 2016 assessment). The IUCN Red List of Threatened Species. 2017;2017:e.T22695189A118592149.Google Scholar