Molecular tools reveal diets of insectivorous birds from predator fecal matter

DNA metabarcoding diet analysis in a generalist omnivore: feeding trials reveal the efficacy of extraction kits and a multi-locus approach for identifying diverse diets

Metabarcoding-based diet analysis is a valuable tool for understanding the feeding behavior of a wide range of species. However, many studies using these methods for wild animals assume accuracy and precision without experimental evaluation with known positive control food items. Here, we conducted a feeding trial experiment with a positive control community in pasture-raised chickens and assessed the efficacy of several commonly used DNA extraction kits and primer sets. We hand-fed 22 known food items, including insects and plants, to six backyard laying hens and collected their excreta for eight h. We evaluated the efficacy of three DNA extraction kits, three primer sets for plant identification (targeting rbcL, trnL, and internal transcribed spacer 2 [ITS2]), and three primer sets for arthropod identification (targeting cytochrome oxidase subunit I [COI]). The detection success rate of our positive control food items was highly variable, ranging from 2.04% to 93.88% for all kit/primer combinations and averaging 37.35% and 43.57% for the most effective kit/primer combination for plants and insects, respectively. Extraction kits using bead-based homogenization positively affected the recovery proportion of plant and insect DNA in excreta samples. The minimum time to detect known food items was 44 min post-feeding. Two COI primer sets significantly outperformed the third, and both recovery proportion and taxonomic resolution from ITS2 were significantly higher than those from rbcL and trnL. Taken together, these results display the potential variability that can be inherently present in DNA-based diet analyses and highlight the utility of experimental feeding trials in validating such approaches, particularly for omnivores with diverse diets.

16S rRNA gene-based microbiota profiles from diverse avian faeces are largely independent of DNA preservation and extraction method

The avian gut microbiota has been the subject of considerable recent attention, with potential implications for diverse fields such as the poultry industry, microbial ecology, and conservation. Faecal microbiotas are frequently used as a non-invasive proxy for the gut microbiota, however the extraction of high-quality microbial DNA from avian faeces has often proven challenging. Here we aimed to evaluate the performance of two DNA preservation methods (95% ethanol and RNAlater) and five extraction approaches (IndiSpin Pathogen Kit, QIAamp PowerFecal Pro DNA Kit, MicroGEM PrepGEM Bacteria Kit, ZymoBIOMICS DNA Miniprep Kit, and an in-house phase separation-based method) for studying the avian gut microbiota. Systematic testing of the efficacy of these approaches on faecal samples from an initial three avian species (chicken, ostrich, and the flightless parrot kākāpō) revealed substantial differences in the quality, quantity and integrity of extracted DNA, but negligible influence of applied method on 16S rRNA gene-based microbiota profiles. Subsequent testing with a selected combination of preservation and extraction method on 10 further phylogenetically and ecologically diverse avian species reiterated the efficacy of the chosen approach, with bacterial community structure clustering strongly by technical replicates for a given avian species. Our finding that marked differences in extraction efficacy do not appear to influence 16S rRNA gene-based bacterial community profiles provides an important foundation for ongoing research on the avian gut microbiota.

DNA metabarcoding quantifies the relative biomass of arthropod taxa in songbird diets: Validation with camera‐recorded diets

Ecological research is often hampered by the inability to quantify animal diets. Diet composition can be tracked through DNA metabarcoding of fecal samples, but whether (complex) diets can be quantitatively determined with metabarcoding is still debated and needs validation using free‐living animals. This study validates that DNA metabarcoding of feces can retrieve actual ingested taxa, and most importantly, that read numbers retrieved from sequencing can also be used to quantify the relative biomass of dietary taxa. Validation was done with the hole‐nesting insectivorous Pied Flycatcher whose diet was quantified using camera footage. Size‐adjusted counts of food items delivered to nestlings were used as a proxy for provided biomass of prey orders and families, and subsequently, nestling feces were assessed through DNA metabarcoding. To explore potential effects of digestion, gizzard and lower intestine samples of freshly collected birds were subjected to DNA metabarcoding. For metabarcoding with Cytochrome Oxidase subunit I (COI), we modified published invertebrate COI primers LCO1490 and HCO1777, which reduced host reads to 0.03%, and amplified Arachnida DNA without significant changing the recovery of other arthropod taxa. DNA metabarcoding retrieved all commonly camera‐recorded taxa. Overall, and in each replicate year (N = 3), the relative scaled biomass of prey taxa and COI read numbers correlated at R = .85 (95CI:0.68–0.94) at order level and at R = .75 (CI:0.67–0.82) at family level. Similarity in arthropod community composition between gizzard and intestines suggested limited digestive bias. This DNA metabarcoding validation demonstrates that quantitative analyses of arthropod diet is possible. We discuss the ecological applications for insectivorous birds.

Advances and Limitations of Next Generation Sequencing in Animal Diet Analysis

Diet analysis is a critical content of animal ecology and the diet analysis methods have been constantly improving and updating. Contrary to traditional methods of high labor intensity and low resolution, the next generation sequencing (NGS) approach has been suggested as a promising tool for dietary studies, which greatly improves the efficiency and broadens the application range. Here we present a framework of adopting NGS and DNA metabarcoding into diet analysis, and discuss the application in aspects of prey taxa composition and structure, intra-specific and inter-specific trophic links, and the effects of animal feeding on environmental changes. Yet, the generation of NGS-based diet data and subsequent analyses and interpretations are still challenging with several factors, making it possible still not as widely used as might be expected. We suggest that NGS-based diet methods must be furthered, analytical pipelines should be developed. More application perspectives, including nutrient geometry, metagenomics and nutrigenomics, need to be incorporated to encourage more ecologists to infer novel insights on they work.

High-throughput sequencing reveals dietary segregation in Malaysian babblers

The coexistence of numerous species within a community results from how those species use available resources. Babblers are one of the major groups of Malaysian insectivorous birds, which frequently forage in dense vegetation cover and have a high level of sympatry. Therefore, examining the diet, prey selection, and niche segregation of babblers can be challenging. In this study, we used high-throughput sequencing to investigate potential dietary overlap or segregation among 10 babbler species of the 4 genera of the family Pellorneidae and Timaliidae: Pellorneum, Malacopteron, Stachyris, and Cyanoderma in central peninsular Malaysia. We tested the hypothesis that trophically similar species may differ in resource use to avoid competitive exclusion. We identified 81 distinct arthropod taxa from fecal samples, belonging to 71 families representing 13 orders, which were predominantly from 16 dipteran, 13 lepidopteran, and 10 coleopteran families. Of all the prey taxa consumed, 45% were found to be distinct across the 10 babbler species, and ˂35% were shared simultaneously by ≥3 babbler species, indicating minimal dietary overlap. The black-throated babbler Stachyris nigricollis and moustached babbler Malacopteron magnirostre had the most generalist tendencies because they consumed a greater variety of prey taxa. Small dietary overlap values (Ojk) and a relatively wide range of food resources suggest that dietary segregation occurred among the studied babblers. The great diversity of prey consumed revealed the presence of dietary flexibility among the sympatric insectivorous birds, thus reducing any active dietary competition and facilitating the coexistence through niche partitioning.

The Rapid Rise of Next-Generation Natural History

Many ecologists have lamented the demise of natural history and have attributed this decline to a misguided view that natural history is outdated and unscientific. Although there is a perception that the focus in ecology and conservation have shifted away from descriptive natural history research and training toward hypothetico-deductive research, we argue that natural history has entered a new phase that we call “next-generation natural history.” This renaissance of natural history is characterized by technological and statistical advances that aid in collecting detailed observations systematically over broad spatial and temporal extents. The technological advances that have increased exponentially in the last decade include electronic sensors such as camera-traps and acoustic recorders, aircraft- and satellite-based remote sensing, animal-borne biologgers, genetics and genomics methods, and community science programs. Advances in statistics and computation have aided in analyzing a growing quantity of observations to reveal patterns in nature. These robust next-generation natural history datasets have transformed the anecdotal perception of natural history observations into systematically collected observations that collectively constitute the foundation for hypothetico-deductive research and can be leveraged and applied to conservation and management. These advances are encouraging scientists to conduct and embrace detailed descriptions of nature that remain a critically important component of the scientific endeavor. Finally, these next-generation natural history observations are engaging scientists and non-scientists alike with new documentations of the wonders of nature. Thus, we celebrate next-generation natural history for encouraging people to experience nature directly.

Nutritional consequences of breeding away from riparian habitats in Bank Swallows: new evidence from multiple endogenous markers

The Bank Swallow (Riparia riparia), a threatened species in Canada, breeds primarily in banks at lakeshores and rivers and in artificial (typically inland) aggregate mining pits. Inland pits may be ecological traps for this species, but relative dietary trade-offs between these two nesting habitats have not been investigated. The availability of aquatic emergent insects at lakeshores may have associated nutritional benefits for growing nestlings due to increased omega-3 fatty acids (FAs) in prey. We compared the diets of juvenile swallows from lakeshore and inland pit sites using assays of stable isotope values (δ¹³C, δ¹⁵N, δ²H) of feathers, faecal DNA metabarcoding and blood plasma FAs. Colony proximity to Lake Erie influenced the use of aquatic versus terrestrial insects by Bank Swallow adults and juveniles. Feather δ²H was particularly useful as a tracer of aquatic emergent versus terrestrial prey, and inland juveniles had feathers enriched in ²H, reflective of diets composed of fewer aquatic emergent insects. DNA metabarcoding of juvenile and adult faecal material indicated that lakeshore birds consumed more aquatic-emergent chironomids than inland birds. Lakeshore juveniles had elevated plasma omega-3 eicosapentaenoic acid levels compared with inland pit-breeding birds. We discuss the need to consider 'nutritional landscapes' and the importance of this concept in conservation of declining species and populations.

Bitter fruits of hard labour: diet metabarcoding and telemetry reveal that urban songbirds travel further for lower-quality food

Rapidly increasing urbanisation requires mitigation against associated losses of biodiversity and species abundance. In urban-breeding birds, altered food availability for nestlings is thought to reduce reproductive success compared to forest populations. To compensate for shortages of preferred foods, urban parents could increase their search effort for optimal diets or provision other foods. Here, we used telemetry and faecal metabarcoding on blue tits from one urban and one forest populations to compare parental effort and comprehensively describe nestling diet. Urban parents travelled on average 30% further than those in the forest, likely to offset limited availability of high-quality nestling food (i.e. caterpillars) in cities. Metabarcoding, based on a mean number of 30 identified taxa per faeces, revealed that the diets of urban chicks were nonetheless substantially shifted to include alternative foods. While in the forest caterpillars comprised 82 ± 11% of taxa provisioned to nestlings, in the city they constituted just 44 ± 10%. Pre-fledging chick mass as well as offspring numbers were lower in urban than in forest-reared broods. Thus, at least in our comparison of two sites, the hard labour of urban parents did not fully pay off, suggesting that improved habitat management is required to support urban-breeding birds.

Gradients in richness and turnover of a forest passerine's diet prior to breeding: A mixed model approach applied to faecal metabarcoding data

Rather little is known about the dietary richness and variation of generalist insectivorous species, including birds, due primarily to difficulties in prey identification. Using faecal metabarcoding, we provide the most comprehensive analysis of a passerine's diet to date, identifying the relative magnitudes of biogeographic, habitat and temporal trends in the richness and turnover in diet of Cyanistes caeruleus (blue tit) along a 39 site and 2° latitudinal transect in Scotland. Faecal samples were collected in 2014-2015 from adult birds roosting in nestboxes prior to nest building. DNA was extracted from 793 samples and we amplified COI and 16S minibarcodes. We identified 432 molecular operational taxonomic units that correspond to putative dietary items. Most dietary items were rare, with Lepidoptera being the most abundant and taxon-rich prey order. Here, we present a statistical approach for estimation of gradients and intersample variation in taxonomic richness and turnover using a generalised linear mixed model. We discuss the merits of this approach over existing tools and present methods for model-based estimation of repeatability, taxon richness and Jaccard indices. We found that dietary richness increases significantly as spring advances, but changes little with elevation, latitude or local tree composition. In comparison, dietary composition exhibits significant turnover along temporal and spatial gradients and among sites. Our study shows the promise of faecal metabarcoding for inferring the macroecology of food webs, but we also highlight the challenge posed by contamination and make recommendations of laboratory and statistical practices to minimise its impact on inference.

Multi-locus DNA metabarcoding of zooplankton communities and scat reveal trophic interactions of a generalist predator

To understand the ecosystem dynamics that underpin the year-round presence of a large generalist consumer, the Bryde's whale (Balaenoptera edeni brydei), we use a DNA metabarcoding approach and systematic zooplankton surveys to investigate seasonal and regional changes in zooplankton communities and if whale diet reflects such changes. Twenty-four zooplankton community samples were collected from three regions throughout the Hauraki Gulf, New Zealand, over two temperature regimes (warm and cool seasons), as well as 20 samples of opportunistically collected Bryde's whale scat. Multi-locus DNA barcode libraries were constructed from 18S and COI gene fragments, representing a trade-off between identification and resolution of metazoan taxa. Zooplankton community OTU occurrence and relative read abundance showed regional and seasonal differences based on permutational analyses of variance in both DNA barcodes, with significant changes in biodiversity indices linked to season in COI only. In contrast, we did not find evidence that Bryde's whale diet shows seasonal or regional trends, but instead indicated clear prey preferences for krill-like crustaceans, copepods, salps and ray-finned fishes independent of prey availability. The year-round presence of Bryde's whales in the Hauraki Gulf is likely associated with the patterns of distribution and abundance of these key prey items.

From feces to data: A metabarcoding method for analyzing consumed and available prey in a bird-insect food web

Diets play a key role in understanding trophic interactions. Knowing the actual structure of food webs contributes greatly to our understanding of biodiversity and ecosystem functioning. The research of prey preferences of different predators requires knowledge not only of the prey consumed, but also of what is available. In this study, we applied DNA metabarcoding to analyze the diet of 4 bird species (willow tits Poecile montanus, Siberian tits Poecile cinctus, great tits Parus major and blue tits Cyanistes caeruleus) by using the feces of nestlings. The availability of their assumed prey (Lepidoptera) was determined from feces of larvae (frass) collected from the main foraging habitat, birch (Betula spp.) canopy. We identified 53 prey species from the nestling feces, of which 11 (21%) were also detected from the frass samples (eight lepidopterans). Approximately 80% of identified prey species in the nestling feces represented lepidopterans, which is in line with the earlier studies on the parids' diet. A subsequent laboratory experiment showed a threshold for fecal sample size and the barcoding success, suggesting that the smallest frass samples do not contain enough larval DNA to be detected by high-throughput sequencing. To summarize, we apply metabarcoding for the first time in a combined approach to identify available prey (through frass) and consumed prey (via nestling feces), expanding the scope and precision for future dietary studies on insectivorous birds.

Statistical assessment of DNA extraction methodology for culture-independent analysis of microbial community associated with diverse environmental samples

Cost-effectiveness, quality, time-effectiveness and ease of the methodology are the most crucial factors in isolating quality DNA from wide variety of samples. Thus, research efforts focusing on the development of an efficient DNA extraction protocol is the need of the hour. The present study therefore, focuses on development of an efficient, rapid and free of inhibitory substances based methodology for extracting metagenomic DNA from diverse environmental samples viz. anaerobic biogas digesta, ruminant stomach, human feces, soil, and microbial starter cultures used for preparation of fermented food. PCR-DGGE based analysis and quality metagenomic library preparation, using DNA extraction methodology, validates the developed protocol. The developed protocol is cost effective, capable of isolating DNA from small sample size (100-1000 µl), time efficient (1.5-2.0 h protocol) and results in significantly higher DNA yield (4-8 times increased yield) when compared to previously available DNA extraction method and a commercial DNA extraction kit. The DNA extracted from the samples using different protocols was evaluated based on its ability to identify diverse microbial species using PCR-DGGE profiles targeting variable region within the 16S rRNA gene. The results of microbial community analysis revealed comparability of the developed protocol to commercial kits, in effectively identifying dominant representatives of the microbial community in different samples. Using the DNA extracted from the presented methodology, metagenomic libraries were prepared, which were found suitable for sequencing on Illumina platform.

The influence of meal size on prey DNA detectability in piscivorous birds

Molecular methods allow noninvasive assessment of vertebrate predator-prey systems at high taxonomic resolution by examining dietary samples such as faeces and pellets. To facilitate the interpretation of field-derived data, feeding trials, investigating the impacts of biological, methodological and environmental factors on prey DNA detection, have been conducted. The effect of meal size, however, has not yet been explicitly considered for vertebrate consumers. Moreover, different noninvasively obtained sample types remain to be compared in such experiments. Here, we present a feeding trial on abundant piscivorous birds, Great Cormorants (Phalacrocorax carbo), to assess meal size effects on postfeeding prey DNA detection success. Faeces and pellets were sampled twice a day after the feed of large (350-540 g), medium (190-345 g) and small (15-170 g) fish meals contributing either a large (>79%) or small (<38%) share to the daily consumption. Samples were examined for prey DNA and fish hard parts. Molecular analysis of faeces revealed that both large meal size and share had a significantly positive effect on prey DNA detection rate postfeeding. Furthermore, large meals were detectable for a significantly longer time span with a detection limit at ~76 hr and a 50% detection probability at ~32 hr postfeeding. In pellets, molecular methods reliably identified the meal consumed the previous day, which was not possible via morphological analysis or when examining individual faeces. The less reliable prey DNA detection of small meals or meal shares in faeces signifies the importance of large numbers of dietary samples to obtain reliable trophic data.

Molecular analysis of faecal samples from birds to identify potential crop pests and useful biocontrol agents in natural areas

Wild habitats adjoining farmland are potentially valuable sources of natural enemies, but also of pests. Here we tested the utility of birds as 'sampling devices', to identify the diversity of prey available to predators and particularly to screen for pests and natural enemies using natural ecosystems as refugia. Here we used PCR to amplify prey DNA from three sympatric songbirds foraging on small invertebrates in Phragmites reedbed ecosystems, namely the Reed Warbler (Acrocephalus scirpaceus), Sedge Warbler (Acrocephalus schoenobaenus) and Cetti's Warbler (Cettia cetti). A recently described general invertebrate primer pair was used for the first time to analyse diets. Amplicons were cloned and sequenced, then identified by reference to the Barcoding of Life Database and to our own sequences obtained from fresh invertebrates. Forty-five distinct prey DNA sequences were obtained from 11 faecal samples, of which 39 could be identified to species or genus. Targeting three warbler species ensured that species-specific differences in prey choice broadened the range of prey taken. Amongst the prey found in reedbeds were major pests (including the tomato moth Lacanobia oleracea) as well as many potentially valuable natural enemies including aphidophagous hoverflies and braconid wasps. Given the mobility of birds, this approach provides a practical way of sampling a whole habitat at once, providing growers with information on possible invasion by locally resident pests and the colonization potential of natural enemies from local natural habitats.

Molecular diagnosis of bird-mediated pest consumption in tropical farmland

Biodiversity loss will likely have surprising and dramatic consequences for human wellbeing. Identifying species that benefit society represents a critical first step towards predicting the consequences of biodiversity loss. Though natural predators prevent billions of dollars in agricultural pest damage annually, characterizing which predators consume pests has proven challenging. Emerging molecular techniques may illuminate these interactions. In the countryside of Costa Rica, we identified avian predators of coffee's most damaging insect pest, the coffee berry borer beetle (Coleoptera:Scolytidae Hypothenemus hampeii), by assaying 1430 fecal samples of 108 bird species for borer DNA. While feeding trials confirmed the efficacy of our approach, detection rates were low. Nevertheless, we identified six species that consume the borer. These species had narrow diet breadths, thin bills, and short wings; traits shared with borer predators in other systems. Borer predators were not threatened; therefore, safeguarding pest control necessitates managing species beyond those at risk of regional extinction by maintaining populations in farmland habitats. Generally, our results demonstrate potential for pairing molecular methods with ecological analyses to yield novel insights into species interactions.

Protocols for metagenomic DNA extraction and Illumina amplicon library preparation for faecal and swab samples

Next-generation sequencing (NGS) technology has extraordinarily enhanced the scope of research in the life sciences. To broaden the application of NGS to systems that were previously difficult to study, we present protocols for processing faecal and swab samples into amplicon libraries amenable to Illumina sequencing. We developed and tested a novel metagenomic DNA extraction approach using solid phase reversible immobilization (SPRI) beads on Western Bluebird (Sialia mexicana) samples stored in RNAlater. Compared with the MO BIO PowerSoil Kit, the current standard for the Human and Earth Microbiome Projects, the SPRI-based method produced comparable 16S rRNA gene PCR amplification from faecal extractions but significantly greater DNA quality, quantity and PCR success for both cloacal and oral swab samples. We furthermore modified published protocols for preparing highly multiplexed Illumina libraries with minimal sample loss and without post-adapter ligation amplification. Our library preparation protocol was successfully validated on three sets of heterogeneous amplicons (16S rRNA gene amplicons from SPRI and PowerSoil extractions as well as control arthropod COI gene amplicons) that were sequenced across three independent, 250-bp, paired-end runs on Illumina's MiSeq platform. Sequence analyses revealed largely equivalent results from the SPRI and PowerSoil extractions. Our comprehensive strategies focus on maximizing efficiency and minimizing costs. In addition to increasing the feasibility of using minimally invasive sampling and NGS capabilities in avian research, our methods are notably not avian-specific and thus applicable to many research programmes that involve DNA extraction and amplicon sequencing.

An holistic ecological analysis of the diet of Cory's shearwaters using prey morphological characters and DNA barcoding

Knowledge of the dietary choices and trophic niches of organisms is the key to understanding their roles in ecosystems. In seabird diet studies, prey identification is a difficult challenge, often yielding results with technique-specific biases. Additionally, sampling efforts are often not extensive enough to reveal intrapopulational variation. Immature animals, which may constitute up to 50% of a population, may occupy a significantly different trophic niche to more experienced birds, but this remains largely unexplored. We investigated the diet of Cory's shearwater (Calonectris diomedea) from Selvagem Grande, an island located off the northwest African coast, collecting a total of 698 regurgitate samples over three consecutive breeding seasons. The diet was assessed using two complementary approaches for prey identification: conventional morphological analysis (using fish vertebrae, otoliths and cephalopod beaks) and DNA barcoding of the 16S rRNA mitochondrial gene, in cases where a positive identification could not be retrieved. Species assignments employed BLAST and distance-based methods, as well as direct optimization of the tree length based on unaligned sequences in POY. This method resulted in robust tree estimates and species assignments, showing its potential for DNA barcoding of stomach contents using hypervariable markers such as the 16S. The molecular approach increased taxonomic resolution and revealed an additional 17 taxa. Diet differed significantly according to breeding status, sex, breeding phase (prelaying and chick rearing) and year. Such direct evidence of trophic segregation within the same population has rarely been shown in seabirds and highlights the importance of including such variables in ecosystem-based management approaches.