- Open Access
The gut of the finch: uniqueness of the gut microbiome of the Galápagos vampire finch
- Alice J. Michel†1,
- Lewis M. Ward†1,
- Shana K. Goffredi2,
- Katherine S. Dawson1, 8,
- Daniel T. Baldassarre3,
- Alec Brenner1,
- Kiyoko M. Gotanda4,
- John E. McCormack2,
- Sean W. Mullin1,
- Ariel O’Neill1,
- Gabrielle S. Tender1,
- J. Albert C. Uy5,
- Kristie Yu1,
- Victoria J. Orphan1Email author and
- Jaime A. Chaves6, 7Email author
© The Author(s). 2018
- Received: 29 May 2018
- Accepted: 5 September 2018
- Published: 19 September 2018
Darwin’s finches are a clade of 19 species of passerine birds native to the Galápagos Islands, whose biogeography, specialized beak morphologies, and dietary choices—ranging from seeds to blood—make them a classic example of adaptive radiation. While these iconic birds have been intensely studied, the composition of their gut microbiome and the factors influencing it, including host species, diet, and biogeography, has not yet been explored.
We characterized the microbial community associated with 12 species of Darwin’s finches using high-throughput 16S rRNA sequencing of fecal samples from 114 individuals across nine islands, including the unusual blood-feeding vampire finch (Geospiza septentrionalis) from Darwin and Wolf Islands. The phylum-level core gut microbiome for Darwin’s finches included the Firmicutes, Gammaproteobacteria, and Actinobacteria, with members of the Bacteroidetes at conspicuously low abundance. The gut microbiome was surprisingly well conserved across the diversity of finch species, with one exception—the vampire finch—which harbored bacteria that were either absent or extremely rare in other finches, including Fusobacterium, Cetobacterium, Ureaplasma, Mucispirillum, Campylobacter, and various members of the Clostridia—bacteria known from the guts of carnivorous birds and reptiles. Complementary stable isotope analysis of feathers revealed exceptionally high δ15N isotope values in the vampire finch, resembling top marine predators. The Galápagos archipelago is also known for extreme wet and dry seasons, and we observed a significant seasonal shift in the gut microbial community of five additional finch species sampled during both seasons.
This study demonstrates the overall conservatism of the finch gut microbiome over short (< 1 Ma) divergence timescales, except in the most extreme case of dietary specialization, and elevates the evolutionary importance of seasonal shifts in driving not only species adaptation, but also gut microbiome composition.
- Galápagos Islands
- Darwin’s finches
First introduced to science by Charles Darwin in his diaries from the voyage of the HMS Beagle , Darwin’s finches are a classic example of adaptive radiation, the phenomenon by which species diverge from a common ancestor as they adapt to different ecological niches. Divergence of Darwin’s finches from relatives in South America has taken place in the last ~ 1.5 My, when ancestral finches first colonized the islands . Since this time, 19 formally recognized finch species have evolved within this clade—18 in the Galápagos and 1 in the Cocos Islands [3, 4], with major radiations occurring in the last 300 ka [5, 6].
The ecological drivers of adaptation in Darwin’s finches have been extensively studied [3, 7–9]. Diet has long been recognized as a major factor in the adaptive radiation of finches in the Galápagos, as lineages on different islands developed beaks specialized for food sources available on their islands, with selective pressure especially high during the dry season when food sources are most limited . Darwin’s finches, as a group, include plant-matter consumers—the seed-eating ground finches (Geospiza fortis, G. fuliginosa, G. magnirostris), herbivorous vegetarian finch (Platyspiza crassirostris), Opuntia nectar- and pollen-specializing cactus finches (Geospiza scandens, G. conirostris)—and insectivorous finches—the sharp-beaked ground finches (G. difficilis, G. acutirostris), woodpecker (Camarhynchus pallidus), tree finch (C. parvulus), warbler finches (Certhidia olivacea, C. fusca) —and the curious blood-feeding vampire finches of Darwin and Wolf Islands (Geospiza septentrionalis) [4, 11, 12]. The vampire finches are so called because they supplement their diet with blood harvested from Nazca and red-footed boobies (Sula granti and S. sula) during the dry season when resources are scarce on these remote islands, a unique strategy not used by any other species [13–15].
Growing evidence suggests that gut microbes exert a major influence on animal nutrition, health, immunity, and behavior [16–18]. Reciprocally, large-scale sequencing of 16S rRNA of associated gut bacteria and archaea has revealed that gut microbial diversity can be influenced by diet, host morphology, host phylogeny, or environment [16, 19, 20]. Comparisons between gut communities of closely related vertebrates with different diets have yielded a range of patterns. Diet has been shown to be a strong predictor of gut microbial composition in phyllostomid bats [21, 22], dolphins , humans , and phylogenetically diverse ant-eating mammals . On the other hand, links between host phylogeny and gut microbiota have been reported from bird species  and in selected mammals, such as bears, where despite their restricted bamboo diet, pandas maintain a microbiome similar to other bears . Other studies have identified a combination of determining factors, including phylogeny and diet (e.g., in some birds and baleen whales [28–30]) and biogeography and diet (e.g., in recently radiated African cichlids and Galápagos iguanas [31–33]). Crucially, studies have focused on clades that diverged tens of millions of years ago, whereas less is known about the degree of gut microbiome divergence in younger, relatively recent species radiations.
Recent studies in birds have generated important data regarding the composition of the gut microbiome, yet there remains an enormous gap in knowledge of most bird lineages , particularly non-domesticated species . The gut microbiome of the domestic chicken [17, 35], folivorous hoatzin , various seabirds including penguin  and petrel , and Passeriformes (the order that includes Darwin’s finches)  are all dominated by members of the Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria, as well as other, low-abundance phyla. A variety of factors influence the gut microbiome, but the role of diet and environment are thought to be of paramount importance in birds , even when controlling for other factors, as shown in fledglings of the brood parasite (Molothrus ater), a cowbird that lays its eggs in the nests of heterospecific hosts .
Here, we analyzed the gut microbiota of 12 species of Darwin’s finches, sampled across 9 of the 18 Galápagos Islands during both the wet and dry seasons. During the wet season, from January to early June, food is abundant, while food becomes limiting and the birds resort to their “adaptive” diet in the dry season, from July to December . Fecal samples were collected from 114 individual finches, including the two remote populations of vampire finch (Geospiza septentrionalis) on Darwin and Wolf Islands. Microbial diversity was screened using Illumina high-throughput sequencing of the 16S rRNA gene. A complementary analysis of stable carbon and nitrogen isotopes of feathers from vampire finches and a subset of others was also performed as an independent means to identify dietary signatures. The diversity in diet and lifestyle of Darwin’s finches and the extensive ecological and genetic knowledge of this avian group make them an attractive study subject to examine how host diet, biogeography, phylogeny, and other environmental factors affect the vertebrate gut microbial community over relatively short evolutionary timescales.
The core gut microbiome of Darwin’s finches
Finch species sampled in this study, including information about general diet category, sampling island, and number of specimens collected during each of the dry and wet season
No. samples (season)
Carnivorous2, inc. blood
Sharp-beaked ground finch
Sharp-beaked ground finch
Medium ground finch
Small ground finch
Large ground finch
Large cactus finch
Herbivorous (Opuntia cactus)
Herbivorous (Opuntia cactus)
Green warbler finch
Small tree finch
The analysis of the core microbiome, consisting of OTUs that were universally shared among the 113 finches (after removal of the single C. pallida), revealed a dominance of Firmicutes (50% average relative abundance, by 16S rRNA recovery), Proteobacteria (40%), and Actinobacteria (8%; Fig. 2; Additional file 3: Figure S1). This pattern has striking resemblance to other avian microbiomes, where surveys of diverse Neotropical birds recovered a similar number of phylotypes (~ 201 bacterial OTUs) belonging to the same three phyla (Proteobacteria, Firmicutes, and Actinobacteria, average relative abundances 46%, 37%, and 1%, respectively) . Two specific OTUs present in every finch included representatives of Enterobacter (17% average abundance; Gammaproteobacteria) and Enterococcus (13% average abundance; Bacilli). A Shigella-affiliated OTU (8% average abundance; Gammaproteobacteria) was also detected in all but one finch from Darwin Island.
Representatives from the Bacteroidetes and five additional phyla (Chloroflexi, Cyanobacteria, Deferribacteres, Fusobacteria, and Tenericutes) were detected, but occurred in only a subset of the finches (typically at < 1% average abundance; Fig. 2; Additional file 1: Table S1). While polysaccharide-degrading Bacteroidetes represent a significant component of the gut microbiome in many vertebrates [16, 38, 41, 42], they were rare, or not detected, in Darwins’ finches (< 0.5% on average in 23/113 finches, and absent in 88 finches), with the exception of 2 finch individuals (C. parvulus, San Cristóbal, wet season and G. conirostris, Genovesa, dry season; Bacteroidetes at 12% and 26%, respectively). Decreased Bacteroidetes presence has been previously observed in birds; however, the previously reported relative abundance of 3–20% is comparably higher than the vast majority of the finches analyzed in this study [20, 39, 40]. Even finches specializing in a carbohydrate-rich diet (e.g., the vegetarian and cactus finches, P. crassirostris and G. scandens) did not show a positive relationship with this bacterial phylum, with only one P. crassirostris and one G. scandens harboring Bacteroidetes OTUs, < 0.02% total microbiota. It is currently unknown which bacterial lineages may occupy a similar niche as Bacteroidetes in the finch microbiome.
Uniqueness of the vampire finch microbiome
Additional taxa that were shown to be statistically different in the vampire finch, relative to other finches, included an OTU affiliated with the genus Campylobacter (Epsilonproteobacteria), at 0.9% average relative abundance versus 0.0005% from other Darwin’s finches (Mann-Whitney p = 0.0002); two Ureaplasma (Tenericutes; Mycoplasmataceae) OTUs at 0.54% (10/31 samples) vs. 0.0001% in the other finches (3/82 samples; Mann-Whitney p = 0.004); and one Mucispirillum (Deferribacteraceae) OTU, accounting for 0.7% of the average vampire finch microbiome (in 17/31 specimens) versus 0.0002% in others (2/82 samples; Mann-Whitney p = 0.000009). The recovered Campylobacter OTU was closely related to Campylobacter volucris, isolated from gulls . Campylobacter have been reported in association with the gut microbiome of some Passeriformes species, as well as other birds, reptiles, and mammals [29, 41], and their distribution and inter-species transmission as pathogens has been well studied . Ureaplasma and Mucispirillum have also been previously described in association with shorebirds , but little is known about their overall distribution and functional role in the host. Given their common occurrence in the vampire finches, it is possible that the acquisition of these lower abundance genera may be linked to the close association of the vampire finch with the co-occurring Nazca and red-footed boobies on Darwin and Wolf Islands.
While other exclusively sanguivorous animals, such as vampire bats and leeches, harbor Aeromonas as an obligate gut symbiont [53, 54], this bacterial taxon was not recovered from the vampire finches, possibly due to the non-obligate nature of their blood feeding. However, a recent investigation of the microbiome of the common vampire bat (Desmodus rotundus) found only very low abundances of Aeromonas sp. (< 0.2%), indicating that Aeromonas might not be essential for obligate sanguivory, as previously presumed .
Blood feeding in Darwin’s finches was highly distinct from the other diet categories of insectivory and plant-based diets (seeds, Opuntia cactus nectar, and leaves) (ANOSIM R = 0.23, p = 0.001; Fig. 4b). The statistical significance of this difference was maintained after accounting for the larger vampire finch sample size (ANOSIM R = 0.38, p = 0.003), sampling only in the dry season (ANOSIM R = 0.34, p = 0.001), as well as focusing specifically on a taxonomically narrow group of finches (previously named G. difficilis), of which the vampires belong (ANOSIM R = 0.57, p = 0.001). In the full dataset, the distinctiveness of the vampire finch species (ANOSIM R = 0.23, p = 0.001) accounts for nearly all of the variation attributable to species differences (R = 0.27, p = 0.001; Fig. 4a; Additional file 5: Table S3), with species differences absent when the vampire finches are omitted (ANOSIM R = 0, p = 0.5). By contrast, no discernable difference between insectivory versus plant-based diets was detected in either the dry season (ANOSIM R = 0.02, p = 0.3) or the wet season (R = 0.01, p = 0.4), or when considering the full dataset, minus vampire finches (R = 0, p = 0.9; Additional file 5: Table S3).
Geographic isolation, host phylogeny, and unusual diet are all possible factors influencing the composition of the gut microbial community in the vampire finch. Darwin and Wolf Islands are geographically isolated from the other islands of the Galápagos (~ 300 km from Santa Cruz, more than double the distance between all other Galápagos Islands), thus limiting genetic exchange with finches within the archipelago . Targeted comparison of the gut microbiome of vampire finches from either Darwin or Wolf Island, separated by a mere 40 km, found them to differ significantly from each other (ANOSIM R = 0.16, p = 0.001; Additional file 6: Figure S2C, Additional file 5: Table S3), due to differences in the abundance of minor taxa (Figs. 2 and 5). This suggests that geographic separation may be a factor influencing the gut microbiome. However, preliminary analysis of a more limited sample set across the remaining islands did not identify a trend. For example, an analysis of G. fuliginosa and G. fortis did not reveal a strong divergence across the islands of San Cristóbal and Santa Cruz (~ 65 km apart; ANOSIM R = 0.06, p = 0.15). Deeper sampling of individual species across islands will be necessary to more confidently examine the potential role of island biogeography.
Additionally, finch phylogeny appeared insignificant for the deepest-sampled non-vampire finches (G. fuliginosa and G. fortis, n = 11 each; ANOSIM R = 0, p = 0.48), suggesting that host phylogeny is not likely a major influence on the gut microbiome of Darwin’s finches. Additional sampling from more distantly related and under-sampled species is necessary to thoroughly clarify the role of phylogeny. Nevertheless, of these potential variables examined, a blood-feeding diet of Geospiza septentrionalis, independently supported by isotopic evidence, appears to be the significant contributor to the unique diversity of the gut microbiome.
Season has a significant influence on the gut microbiome of Darwin’s finches
Investigations of Darwin’s finches provided a unique opportunity to compare the influence of extreme seasonal shifts on the gut microbiome. The Galápagos Islands are well known for strong seasonal cycles driving the availability of vegetation and other food sources utilized by Darwin’s finches, and it is during the dry season that natural selection for unique feeding adaptations is most intense [3, 55]. It is thought that food limitation during the dry season induces specialization and speciation of the various finch species, as species adapt to recalcitrant resources such as Tribulus seeds, a process documented in Darwin’s ground finches especially during extreme El Niño events [3, 10]. Plant growth and seed production occurs during the wet season (January to June), which is also associated with an increase in arthropod abundance . During the dry season, however, finches undergo significant dietary stress as plant-based food sources become limiting  and it is then that speciation becomes observable . Season has been identified as a significant driver of gut microbiome divergence in other animals, including alligators, which engage in protracted fasting during the winter . Whether the microbiome of Darwin’s finches responds to seasonal extremes on the Galápagos Islands is an outstanding question.
Using our dataset consisting of 72 individuals from the dry season and 41 from the wet season, we found season to be significantly associated with differences in microbial composition. Statistical differences in the finch gut microbiome between the wet and dry seasons occurred across the entire dataset (ANOSIM R = 0.27, p = 0.001; Fig. 4d), as well as when vampire finches were excluded (ANOSIM R = 0.17, p = 0.001). Seasonal differentiation was also detected strictly within the two finch species sampled during both seasons (G. fuliginosa and G. fortis) (ANOSIM R = 0.28, p = 0.001; Additional file 6: Figure S2D, Additional file 5: Table S3). Members of the Gammaproteobacteria were the primary driver of the variation in finch microbiome during the dry season (compositional abundance 48% in dry vs. 21% in wet, Mann-Whitney U test p = 0.01), while Bacilli were dominant during the wet season (35% in dry vs. 65% in wet, Mann-Whitney p = 0.00008) (Fig. 2). Higher relative abundances of Actinobacteria were also documented in the wet season (6% vs. 12%), although this difference was not significant (Mann-Whitney p = 0.5).
A significant shift in the finch gut microbial community with season may reflect temporal shifts in response to changes in weather, food resource availability, or other factors. The seasonal variation in the gut microbiome among all granivorous/herbivorous finches (G. magnirostris, G. fortis, G. fuliginosa, G. scandens, G. conirostris, P. crassirostris) supports the ecological observations of finches optimizing their feeding strategies by season (ANOSIM R = 0.21, p = 0.001; n = 50). Omnivorous/insectivorous finches (G. difficilis, G. acutirostris, C. parvulus, and Ce. olivacea), on the other hand, were not significantly different by season (ANOSIM R = 0.08, p = 0.2; n = 20), consistent with lesser shifts in types of foodstuffs between seasons. It is conceivable, however, that either hormonal or ecological changes coincident with season cause the microbiome to shift in composition, irrespective of diet. For instance, when the finches nest during the wet season, they may come into contact with different microorganisms, such as those carried by the invasive nest parasite, Philornis downsii . Firmicutes and Proteobacteria were recently reported as the dominant taxa in the microbiome of P. downsii on Santa Cruz Island .
A preliminary assessment of diet specialization in a limited number of these finches (broadly categorized as granivory, omnivory, insectivory, nectarivory, and herbivory) revealed no obvious relationship with the gut microbiome, despite well-documented shifts in finch beak morphology and behavior. Other bird species have shown similar patterns, where only the most extreme diets (e.g., vultures, and other scavengers) resulted in a marked difference in the microbiome . Most of Darwin’s finch species have remained omnivorous, with an emphasis on seeds; thus, the influence of diet specialization on the gut microbiome in these species is expected to be minor compared to other animal groups where more extreme dietary specialization has led to shifts in composition (e.g., phyllostomid bats [21, 22] and cichlid fishes ). Additional sampling—particularly of the seed-specializing ground finches, G. fuliginosa, G. fortis, and G. magnirostris, which are distributed across multiple islands, as well as dietary specialists such as the vegetarian P. crassirostris and insectivorous Ce. olivacea—may reveal more subtle trends driving gut microbiome composition.
Darwin’s finches have long been recognized as a model system for investigating interactions between biogeography, dietary specialization, morphology, and other aspects of evolutionary radiation. To this list can now be added the interaction of these factors with the composition of the gut microbial community. The radiation of Darwin’s finches has occurred relatively rapidly and recently (~ 300 ka). During this time, it appears that the gut microbial composition of most finch species has remained conserved among the majority of species examined in this study. A conspicuous exception to this conservation occurs in the case of the vampire finch, Geospiza septentrionalis, with a microbial community significantly distinct from the other finches. Because of the barrenness of Darwin and Wolf Islands, the vampire finches undergo extreme dietary limitations during the dry season, thus leading to blood feeding . Their unique and carnivore-like diet likely contributes to the distinct gut microbial composition, suggesting that cases of extreme dietary divergence can overpower phylogenetic inertia to drive shifts in microbiome composition, even over short evolutionary timescales. Future studies in other recent avian radiations that include dietary extremes, such as the blood-feeding Galápagos Hood mockingbird (Mimus macdonaldi), should help with investigating this hypothesis.
A robust seasonal difference was observed between finches collected during the wet and dry seasons and may reflect temporal shifts in response to changes in food resource availability, weather, mating or nesting, or other factors. Further sampling—particularly resampling individual finches through time, across seasons—may help to resolve these uncertainties. It will be interesting to see whether the trends observed in this year-long wet-to-dry cycle is reflected in longer term datasets or with more extreme weather periods experienced, for example, during El Niño/La Niña conditions. This goal will be challenging given their remote nature, enhanced diversity, and numerous overlapping influences. Regardless, Darwin’s finches continue to capture our attention and reveal new secrets after more than a century and a half of study.
Individuals (n = 114) from 12 species of Darwin’s finches, from nine islands over both the wet and dry seasons, were captured in mist nets and their fecal material sampled (Table 1, Fig. 1). Dry season samples were collected during November–December 2015, while wet season samples were collected between March and June 2016. Fecal samples were collected by placing the bird in a paper bag on top of a metal mesh grate above a sterilized weigh boat for 3-5 min . Samples were immediately transferred into ~ 2 ml of LifeGuard Preservation Solution (MoBio Laboratories, Carlsbad, CA, USA). Feather samples from three species (G. fortis, G. fuliginosa, and G. septentrionalis) were collected for isotopic analysis in Jan 2016 (Fig. 3a).
DNA extraction and microbial community analysis
Prior to extracting DNA, samples were centrifuged to remove the LifeGuard preservation solution. Genomic DNA was then recovered from the pelleted fecal material using a PowerSoil DNA isolation kit (MoBio Laboratories, Carlsbad, CA, USA) following the manufacturer’s instructions, with the exception of the addition of a bead-beating step using a Fast Prep 120 instrument; Thermo Electron Corporation). The V4 region of the 16S rRNA gene was PCR-amplified from each extract using the archaeal and bacterial targeted primer set 515F and 806R , following the protocol outlined by Case et al. . Successful PCR amplifications were pooled, in duplicate, and barcodes were added according to the Earth Microbiome Project protocol [59, 61, 62]; 5 μl of the amplicon product from the first PCR was used as template in a 5-cycle, 25-μl reconditioning reaction with the same EMP-recommended conditions and the full EMP primers. Samples were mixed together in equimolar amounts and purified in bulk through a Qiagen PCR Purification kit. At all PCR steps, amplification success and purity was checked by gel electrophoresis. Paired-end sequences (2x 250 basepair) were generated from barcoded amplicon products at Laragen, Inc. on an Illumina MiSeq platform. At Laragen, the raw data was passed through a filter which demultiplexed the library into individual samples and removed any sequences which had > 1 basepair mismatch on the 12-basepair barcode sequence, and assigned quality scores to each basepair call on every sequence. At the same time, adapter, barcode, and primer sequences were removed. Raw reads were deposited and are available through the Sequence Read Archive under accession number SRP130314.
Sequence data was processed with both the DADA2 pipeline for unclustered high-resolution exact sequence variants (SVs)  and QIIME version 1.8.0  for operational taxonomic units (OTUs) at the 97% similarity level. In QIIME, raw sequence pairs were joined and quality-trimmed using the default parameters. Sequences were clustered into de novo OTUs using UCLUST open reference clustering protocol . Then, the most abundant sequence was chosen as representative for each de novo OTU . Taxonomic identification for each representative sequence was assigned using the Silva-119 database . For downstream analyses, the QIIME 97% OTU dataset was trimmed to only include those OTUs representing at least 1% of the total gut community of at least one finch (Additional file 1: Table S1). DADA2 exact sequence variants were calculated using the published pipeline tutorial (v1.4) in R and were analyzed the same as QIIME results, without the 1% abundance cutoff (errors are already accounted for in the DADA2 error model ).
Sequences from the 16S rRNA region were aligned using MAFFT  and a phylogeny constructed using FastTree . Alpha diversity was estimated using the Shannon index, Simpson index, Fisher alpha, and Chao1 richness and by rarefaction. Diversity statistics were calculated in “phyloseq” v-1.25.2  and are reported at the 97% OTU level in Additional file 2: Table S2. Since results across phylotype extraction methods were highly similar, all in-text figures, values, and statistics are derived from the 97% OTU dataset.
To assess the influence of various factors on the microbiome, the 114-finch dataset was divided 19 different ways, each subdivision accounting for one or more of six different ecological and biological factors. One individual (C. pallida) was removed from all analyses because it was the sole representative of its species, leaving 113 samples for beta diversity analysis. Each division of the dataset was balanced in sampling within the factor in question, unless otherwise noted. Where category sample sizes differed, the data subdivision was randomly subset by each factor category down to the sample size of the smallest category, with statistics taken from the average of 10 random subsets, as described in Additional file 7: Table S4. The influence of sampling season, island, finch diet, finch species, and the extreme diet of sanguivory were assessed using these subdivisions (Additional file 7: Table S4).
Statistical analysis of the role of each factor in shaping finch microbiome beta diversity was done within each data subdivision on relative (compositional) abundance of each OTU in each finch using non-metric multidimensional scaling (NMDS) ordination plots, analysis of similarities (ANOSIM) of both Bray-Curtis and UniFrac dissimilarities, and PERMANOVA (Adonis) of both Bray-Curtis and UniFrac dissimilarities. Dissimilarity matrices derived from both unweighted (presence-absence) and weighted taxonomic relative abundance (scaled to the fourth-root) were calculated for data from each finch. Summary statistics for all data subdivisions, dissimilarities, and approaches listed above can be found in Additional file 5: Table S3. For data subdivisions where statistical balance was achieved via sub-setting one or more categories within a factor, random subsets were run 10 times and the summary statistics reported represent the average over each run.
The contribution of individual OTUs to statistically significant grouping by factor was queried using similarity percentage (SIMPER), and the significance of differential abundance in individual OTUs by factor grouping were statistically tested using two-tailed Mann-Whitney U tests. In an approach similar to beta diversity tests, for data subdivisions where random sub-setting was necessary to achieve equal sample sizes within each factor category, Mann-Whitney tests were run over 100 random subsets of the larger category size, and the average p value was reported. Summary statistics and plots were generated in R using the packages “phyloseq,” “vegan,” “ggplot2,” and “RColorBrewer” [70–74].
Carbon, nitrogen, and sulfur isotope analysis of feathers
The isotopic signature of an organism—including that of δ13C, δ15N, and δ34S of proteins in animal tissues such as hair or feathers—is primarily influenced by food source and is frequently used to discern differences in diet (e.g., ). Feathers were collected from a subset of finch species during the transition between dry and wet seasons (Jan 2016). Prior to isotopic analysis, feathers were immersed in 2:1 dichloro-methane:methanol to remove surface oils as described by Blight et al. . After air-drying, feather length and weight were recorded and then individual feathers were split along the rachis and transferred to tin capsules (0.3–1.0 mg), for parallel carbon, nitrogen, and sulfur isotope analysis. Feather δ13C and δ15N and weight percent C and N (wt.% TOC and TON) were determined via continuous flow (He; 100 ml/min) on a Costech Instruments Elemental Combustion System model 4010 by oxidation at 980 °C over chromium (III) oxide and silvered cobalt (II, III) oxide followed by reduction over elemental copper at 650 °C. CO2 was subsequently passed through a water trap and then a 5-Å molecular sieve GC column at 50 °C to separate N2 from CO2, which was diluted with helium in a Conflo IV interface/open split prior to analysis. Fast jump was calibrated and applied to measure both CO2 and N2 in the same run. δ13C and δ15N values were measured on a Thermo Scientific Delta V Plus irMS. δ13C and δ15N values were corrected for sample size dependency and then normalized to the VPDB scale with a two-point calibration . Error was determined by analyzing sucrose (NIST 8542), acetanilide (Costech Analytical Technologies Inc.), and nitrate (IAEA-NO-3) in combination with in-house standards (from − 45.93‰ to − 10.45‰ for C and − 3.02‰ to 4.70‰ for N). For C measurements, accuracy was ± 0.18‰ (n = 50) and precision was ± 0.23‰ (n = 50; 1σ). For N measurements, accuracy was ± 0.22‰ (n = 36) and precision was ± 0.43‰ (n = 36; 1σ). Feather δ34S and weight percent S (wt.% TOS) was measured using continuous flow (He; 120 ml/min) on a ThermoQuest NC2500 mass spectrometer (ThermoQuest Italia, Milan, Italy) by oxidation at 1000 °C over tungsten (VI) oxide and elemental copper. SO2 was subsequently passed through a water trap and then a 5-Å molecular sieve GC column at 85 °C. A Conflo III interface/open split was used to introduce SO2 to a Thermo Scientific Delta plus XL irMS. The δ34S values were corrected for sample size dependency and the normalized to the VCDT scale with a two-point calibration. Error was determined by analyzing two silver sulfides (IAEA-S-2 and IAEA-S-3) in combination with in-house standards including pyrite and sulfanilamide. Accuracy on the δ34S measurement was ± 1.84‰ (n = 20) and precision was ± 0.35‰ (n = 20; 1σ).
We are grateful to Sofia Carvajal, Kurt Gielow, Sarah Knutie, Andres Leon, Simón Villamar, Angela Hansen, Sabrina McNew, Ashley Saulsberry, Carlos Vera, Ruben Heleno, Manuel Nogales, and Sandra Hervías for the invaluable assistance in the field; Alexis Pasulka, David Case, and Stephanie Connon for the assistance in lab and advice with the data analysis, as well as other members of the Caltech GE/Bi/ESE_246 molecular geobiology and Bi/GE_105 Evolution course including Rob Phillips, Courtney Chen, Anne Christian, Bianca Lepe, Kristin Anderson, Tristan Murphy, Matthew Smalley, and Tiffany Zhou; and to Parque Nacional Galápagos and Galápagos Science Center for their logistical support. All samples were collected under the Contrato Marco de Acceso a los Recursos Genéticos MAE-DNB-CM-2016-0041 and Galápagos National Park permit PC-0615.
Funds were provided by GAIAS-USFQ Grant to JAC, the Gordon and Betty Moore Foundation through grant GMBF3780 to VJO, NSF-Postdoctoral Grant to DTB, Instrumental Crowdfunding, Le Fonds Québécois de la Recherche sur la Nature et les Technologies Postdoctoral Fellowship, a British Ornithologists’ Union Research Grant to KMG, and a NASA Earth and Space Science Fellowship to LMW. We also acknowledge support from Terence Barr and Caltech's Center for Environmental Microbial Interactions (CEMI).
Availability of data and materials
Raw reads were deposited and are available through the Sequence Read Archive under accession number SRP130314.
AJM, LMW, SKG, and VJO analyzed the data and wrote the manuscript. KSD collected the stable isotope data and assisted with the interpretation. AJM, LMW, SKG, VJO, SM extracted DNA and prepared samples for sequencing. AB, JEM, SM, AO, GT, and KY helped analyze the data. DTB, KMG, ACU, and JAC collected fecal samples and observation data on finches on the Galápagos and assisted with the writing of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
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.
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.
- Darwin C. Journal of Researches into the Geology and Natural History of the various countries visited by H.M.S. Beagle, under the command of Captain Fitzroy, R.N. from 1832 to 1836. London: Henry Colburn; 1839.Google Scholar
- Petren K, Grant PR, Grant BR, Keller LF. Comparative landscape genetics and the adaptive radiation of Darwin’s finches: the role of peripheral isolation. Mol Ecol. 2005;14:2943–57.View ArticleGoogle Scholar
- Grant PR and Grant RB. 40 Years of Evolution: Darwin's Finches on Daphne Major Island. Princeton: Princeton University Press; 2014.Google Scholar
- Remsen JV, Jr, Areta JI, Cadena CD, Claramunt S, Jaramillo A, Pacheco JF, Robbins MB, Stiles FG, Stotz DF, and Zimmer KJ. A classification of the bird species of South America. American Ornithologists' Union. Version January 2015. http://www.museum.lsu.edu/~Remsen/SACCBaseline.htm.
- Lamichhaney S, Berglund J, Almén MS, Maqbool K, Grabherr M, Martinez-Barrio A, et al. Evolution of Darwin’s finches and their beaks revealed by genome sequencing. Nature. 2015;518:371–5.View ArticleGoogle Scholar
- Lamichhaney S, Han F, Webster MT, Andersson L, Grant BR, Grant PR. Rapid hybrid speciation in Darwin’s finches. Science. 2018;359:224.View ArticleGoogle Scholar
- Grant PR. Ecology and evolution of Darwin’s finches. Princeton: Princeton University Press; 1999.Google Scholar
- Lack D. Subspecies and sympatry in Darwin’s finches. Evolution. 1969;23:252–63.View ArticleGoogle Scholar
- Lack D. Evolution of the Galápagos finches. Nature. 1940;146:324–7.View ArticleGoogle Scholar
- De León LF, Podos J, Gardezi T, Herrel A, Hendry AP. Darwin’s finches and their diet niches: the sympatric coexistence of imperfect generalists. J Evol Biol. 2014;27:1093–104.View ArticleGoogle Scholar
- Rothschild W and Hartert E. A Review of the Ornithology of the Galapagos Islands. With Notes on the Webster-Harris Expedition. Novitates Zoologicae. 1899; Vol.VI(No.2):pp.85-205, 2 plates.Google Scholar
- Farrington HL, Lawson LP, Clark CM, Petren K. The evolutionary history of Darwin’s finches: speciation, gene flow, and introgression in a fragmented landscape. Evolution. 2014;68:2932–44.View ArticleGoogle Scholar
- Bowman RI, Billeb SI. Blood-eating in a Galápagos finch. Living Bird. 1968;4:29–44.Google Scholar
- Koster F, Koster H. Twelve days among the “vampire finches” of Wolf Island. Not Galápagos. 1983;38:4–10.Google Scholar
- Schluter D, Grant PR. Ecological correlates of morphological evolution in a Darwin’s finch, Geospiza Difficilis. Evolution. 1984;38:856–69.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
- Stanley D, Hughes RJ, Moore RJ. Microbiota of the chicken gastrointestinal tract: influence on health, productivity and disease. Appl Microbiol Biotechnol. 2014;98:4301–10.View ArticleGoogle Scholar
- Heijtz RD, Wang S, Anuar F, Qian Y, Björkholm B, Samuelsson A, et al. Normal gut microbiota modulates brain development and behavior. Proc Natl Acad Sci. 2011;108:3047–52.View ArticleGoogle Scholar
- Spor A, Koren O, Ley R. Unravelling the effects of the environment and host genotype on the gut microbiome. Nat Rev Microbiol. 2011;9:279.View ArticleGoogle Scholar
- Waite D, Taylor M. Exploring the avian gut microbiota: current trends and future directions. Front Microbiol. 2015;6:673.View ArticleGoogle Scholar
- Carrillo-Araujo M, Taş N, Alcántara-Hernández RJ, Gaona O, Schondube JE, Medellín RA, et al. Phyllostomid bat microbiome composition is associated to host phylogeny and feeding strategies. Front Microbiol. 2015;6:447.View ArticleGoogle Scholar
- Mendoza MLZ, Xiong Z, Escalera-Zamudio M, Runge AK, Thézé J, Streicker D, et al. Hologenomic adaptations underlying the evolution of sanguivory in the common vampire bat. Nat Ecol Evol. 2018;2:659–68.View ArticleGoogle Scholar
- Soverini M, Quercia S, Biancani B, Furlati S, Turroni S, Biagi E, et al. The bottlenose dolphin (Tursiops truncatus) faecal microbiota. FEMS Microbiol Ecol. 2016;92:fiw055.View ArticleGoogle Scholar
- Muegge BD, Kuczynski J, Knights D, Clemente JC, González A, Fontana L, et al. Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science. 2011;332:970–4.View ArticleGoogle Scholar
- Delsuc F, Metcalf JL, Wegener Parfrey L, Song SJ, González A, Knight R. Convergence of gut microbiomes in myrmecophagous mammals. Mol Ecol. 2014;23:1301–17.View ArticleGoogle Scholar
- Waite DW, Taylor MW. Characterizing the avian gut microbiota: membership, driving influences, and potential function. Front Microbiol. 2014;5:223.View ArticleGoogle Scholar
- Xue Z, Zhang W, Wang L, Hou R, Zhang M, Fei L, et al. The bamboo-eating giant panda harbors a carnivore-like gut microbiota, with excessive seasonal variations. mBio. 2015;6:e00022–15.PubMedPubMed CentralGoogle Scholar
- Hird SM, Sánchez C, Carstens BC, Brumfield RT. Comparative gut microbiota of 59 neotropical bird species. Front Microbiol. 2015;6:1403.View ArticleGoogle Scholar
- Maul JD, Gandhi JP, Farris JL. Community-level physiological profiles of cloacal microbes in songbirds (Order: Passeriformes): variation due to host species, host diet, and habitat. Microb Ecol. 2005;50:19–28.View ArticleGoogle Scholar
- Sanders JG, Beichman AC, Roman J, Scott JJ, Emerson D, JJ MC, et al. Baleen whales host a unique gut microbiome with similarities to both carnivores and herbivores. Nat Commun. 2015;6:ncomms9285.View ArticleGoogle Scholar
- Baldo L, Riera JL, Tooming-Klunderud A, Albà MM, Salzburger W. Gut microbiota dynamics during dietary shift in eastern African cichlid fishes. PLOS ONE. 2015;10:e0127462.View ArticleGoogle Scholar
- Hong P-Y, Wheeler E, Cann IKO, Mackie RI. Phylogenetic analysis of the fecal microbial community in herbivorous land and marine iguanas of the Galápagos Islands using 16S rRNA-based pyrosequencing. ISME J. 2011;5:1461–70.View ArticleGoogle Scholar
- Lankau EW, Hong P-Y, Mackie RI. Ecological drift and local exposures drive enteric bacterial community differences within species of Galápagos iguanas. Mol Ecol. 2012;21:1779–88.View ArticleGoogle Scholar
- Hird SM. Evolutionary biology needs wild microbiomes. Front Microbiol. 2017;8:725.View ArticleGoogle Scholar
- Oakley BB, Lillehoj HS, Kogut MH, Kim WK, Maurer JJ, Pedroso A, et al. The chicken gastrointestinal microbiome. FEMS Microbiol Lett. 2014;360:100–12.View ArticleGoogle Scholar
- Godoy-Vitorino F, Goldfarb KC, Karaoz U, Leal S, Garcia-Amado MA, Hugenholtz P, et al. Comparative analyses of foregut and hindgut bacterial communities in hoatzins and cows. ISME J. 2012;6:531.View ArticleGoogle Scholar
- Dewar ML, Arnould JPY, Dann P, Trathan P, Groscolas R, Smith S. Interspecific variations in the gastrointestinal microbiota in penguins. Microbiology Open. 2013;2:195–204.View ArticleGoogle Scholar
- Pearce DS, Hoover BA, Jennings S, Nevitt GA, Docherty KM. Morphological and genetic factors shape the microbiome of a seabird species (Oceanodroma leucorhoa) more than environmental and social factors. Microbiome. 2017;5:146.View ArticleGoogle Scholar
- Kropáčková L, Těšický M, Albrecht T, Kubovčiak J, Čížková D, Tomášek O, et al. Codiversification of gastrointestinal microbiota and phylogeny in passerines is not explained by ecological divergence. Mol Ecol. 2017;26:5292–304.View ArticleGoogle Scholar
- Hird SM, Carstens BC, Cardiff SW, Dittmann DL, Brumfield RT. Sampling locality is more detectable than taxonomy or ecology in the gut microbiota of the brood-parasitic Brown-headed Cowbird (Molothrus ater). PeerJ. 2014;2:e321.View ArticleGoogle Scholar
- Colston TJ, Jackson CR. Microbiome evolution along divergent branches of the vertebrate tree of life: what is known and unknown. Mol Ecol. 2016;25:3776–800.View ArticleGoogle Scholar
- Benskin CMH, Rhodes G, Pickup RW, Wilson K, Hartley IR. Diversity and temporal stability of bacterial communities in a model passerine bird, the zebra finch. Mol Ecol. 2010;19:5531–44.View ArticleGoogle Scholar
- Kelly JF. Stable isotopes of carbon and nitrogen in the study of avian and mammalian trophic ecology. Can J Zool. 2000;78:1–27.View ArticleGoogle Scholar
- Polischuk SC, Hobson KA, Ramsay MA. Use of stable-carbon and -nitrogen isotopes to assess weaning and fasting in female polar bears and their cubs. Can J Zool. 2001;79:499–511.View ArticleGoogle Scholar
- Schoeninger MJ, DeNiro MJ. Nitrogen and carbon isotopic composition of bone collagen from marine and terrestrial animals. Geochim Cosmochim Acta. 1984;48:625–39.View ArticleGoogle Scholar
- Hobson KA, Clark RG. Assessing avian diets using stable isotopes I: turnover of 13C in tissues. The Condor. 1992;94:181–8.View ArticleGoogle Scholar
- Kohl KD, Dearing MD. Experience matters: prior exposure to plant toxins enhances diversity of gut microbes in herbivores. Ecol Lett. 2012;15:1008–15.View ArticleGoogle Scholar
- Keenan S, Engel A, M Elsey R. The alligator gut microbiome and implications for archosaur symbioses. Sci Rep. 2013;3:2877.View ArticleGoogle Scholar
- Roggenbuck M, Schnell IB, Blom N, Bælum J, Bertelsen MF, Sicheritz-Pontén T, et al. The microbiome of New World vultures. Nat Commun. 2014;5:ncomms6498.View ArticleGoogle Scholar
- Debruyne L, Broman T, Bergström S, Olsen B, On SLW, Vandamme P. Campylobacter volucris sp. nov., isolated from black-headed gulls (Larus ridibundus). Int J Syst Evol Microbiol. 2010;60:1870–5.View ArticleGoogle Scholar
- Wassenaar TM, Newell DG. The genus Campylobacter. In: Dworkin M, Falkow S, Rosenberg E, Schleifer K-H, Stackebrandt E, editors. Prokaryotes Vol 7 Proteobacteria Delta Epsil Subclass. New York, NY: Springer New York; 2006. p. 119–38.Google Scholar
- Ryu H, Grond K, Verheijen B, Elk M, Buehler DM, Domingo JWS. Intestinal microbiota and species diversity of Campylobacter and Helicobacter spp. in migrating shorebirds in Delaware Bay. Appl Environ Microbiol. 2014;80:1838–47.View ArticleGoogle Scholar
- Graf J. Symbiosis of Aeromonas veronii Biovar sobria and Hirudo medicinalis, the medicinal leech: a novel model for digestive tract associations. Infect Immun. 1999;67:1–7.PubMedPubMed CentralGoogle Scholar
- Muller HEM, Pinus M, Schmidt U. Aeromonas hydrophila as a normal intestinal bacterium of the vampire bat (Desmodus rotundus). Zentralblatt Vet. 1980;5:419–24.Google Scholar
- Grant PR, Boag PT. Rainfall on the Galápagos and the demography of Darwin’s finches. Auk. 1980;97:227–44.Google Scholar
- Kleindorfer S, Dudaniec RY. Host-parasite ecology, behavior and genetics: a review of the introduced fly parasite Philornis downsi and its Darwin’s finch hosts. BMC Zool. 2016;1:1.View ArticleGoogle Scholar
- Ben-Yosef M, Zaada DSY, Dudaniec RY, Pasternak Z, Jurkevitch E, Smith RJ, et al. Host-specific associations affect the microbiome of Philornis downsi, an introduced parasite to the Galápagos Islands. Mol Ecol. 2017;26:4644–56.View ArticleGoogle Scholar
- Knutie SA, Gotanda KM. A Non-invasive Method to Collect Fecal Samples from Wild Birds for Microbiome Studies. Microb Ecol. 2018;5:1–5.Google Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.View ArticleGoogle Scholar
- Case DH, Pasulka AL, Marlow JJ, Grupe BM, Levin LA, Orphan VJ. Methane seep carbonates host distinct, diverse, and dynamic microbial assemblages. mBio. 2015;6:e01348–15.View ArticleGoogle Scholar
- Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci. 2011;108:4516.View ArticleGoogle Scholar
- Gilbert JA, Meyer F, Jansson J, Gordon J, Pace N, Tiedje J, et al. The Earth Microbiome Project: meeting report of the “1st EMP meeting on sample selection and acquisition” at Argonne National Laboratory October 6th 2010. Stand Genomic Sci. 2010;3:249–53.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.View ArticleGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticleGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.View ArticleGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.View ArticleGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.View ArticleGoogle Scholar
- Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLOS ONE. 2010;5:e9490.View ArticleGoogle Scholar
- McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.View ArticleGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara R, et al. Vegan: community ecology package. R package version. 2016;2:3–5.Google Scholar
- Wickham, H. ggplot2: elegant graphics for data analysis. J Stat Softw, 2010;35:65–88.Google Scholar
- Neuwirth, E. and Neuwirth, M.E. Package ‘RColorBrewer’. 2011. CRAN 2011-06-17 08: 34: 00. Apache License 2.0Google Scholar
- R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2014. R Foundation for Statistical Computing. 2016Google Scholar
- Post DM. Using stable isotopes to estimate trophic position: models, methods, and assumptions. Ecology. 2002;83:703–18.View ArticleGoogle Scholar
- Blight LK, Hobson KA, Kyser KT, Arcese P. Changing gull diet in a changing world: A 150-year stable isotope (C, N) record from feathers collected in the Pacific Northwest of North America. Glob Chang Biol. 2015;21(4):1497–507.View ArticleGoogle Scholar
- Coplen TB, Brand WA, Gehre M, Gröning M, Meijer HAJ, Toman B, et al. New guidelines for δ13C measurements. Anal Chem. 2006;78:2439–41.View ArticleGoogle Scholar
- Hobson KA, Piatt JF, Pitocchelli J. Using stable isotopes to determine seabird trophic relationships. J Anim Ecol. 1994;63:786–98.View ArticleGoogle Scholar
- Hobson KA. Stable-carbon and nitrogen isotope ratios of songbird feathers grown in two terrestrial biomes: implications for evaluating trophic relationships and breeding origins. The Condor. 1999;101:799–805.View ArticleGoogle Scholar
- Roth JD, Hobson KA. Stable carbon and nitrogen isotopic fractionation between diet and tissue of captive red fox: implications for dietary reconstruction. Can J Zool. 2000;78:848–52.View ArticleGoogle Scholar
- Steinfurth A. Marine ecology and conservation of the Galápagos penguin, Spheniscus mendiculus. CAU Kiel: Christian-Albrechts Universität Kiel. 2007.Google Scholar
- Seminoff JA, Jones TT, Eguchi T, Hastings M, Jones DR. Stable carbon and nitrogen isotope discrimination in soft tissues of the leatherback turtle (Dermochelys coriacea): insights for trophic studies of marine turtles. J Exp Mar Biol Ecol. 2009;381:33–41.View ArticleGoogle Scholar
- Matich P, Heithaus MR, Layman CA. Size-based variation in intertissue comparisons of stable carbon and nitrogen isotopic signatures of bull sharks (Carcharhinus leucas) and tiger sharks (Galeocerdo cuvier). Can J Fish Aquat Sci. 2010;67:877–85.View ArticleGoogle Scholar
- Valenzuela LO, Chesson LA, O’Grady SP, Cerling TE, Ehleringer JR. Spatial distributions of carbon, nitrogen and sulfur isotope ratios in human hair across the central United States. Rapid Commun Mass Spectrom. 2011;25:861–8.View ArticleGoogle Scholar
- Jeglinski JWE, Werner C, Robinson PW, Costa DP, Trillmich F. Age, body mass and environmental variation shape the foraging ontogeny of Galápagos sea lions. Mar Ecol Prog Ser. 2012;453:279–96.View ArticleGoogle Scholar
- Zarate P, Bjorndal KA, Seminoff J, Bolten A. Understanding migratory and foraging behavior of green turtles Chelonia mydas in the Galápagos Islands through stable isotopes. 2012. Conference: Thirty-first Annual Symposium on Sea Turtle Biology and Conservation. San Diego: Volume: NOAA Technical Memorandum NOAA NMFS-SEFSC-631.Google Scholar
- Connan M, Hofmeyr GJG, Pistorius PA. Reappraisal of the trophic ecology of one of the world’s most threatened spheniscids, the African penguin. PLOS ONE. 2016;11:e0159402.View ArticleGoogle Scholar