Proportion of prey consumed can be determined from faecal DNA using real-time PCR

Combining vertebrate mitochondrial 12S rRNA gene sequencing and shotgun metagenomic sequencing to investigate the diet of the leopard cat (Prionailurus bengalensis) in Korea

The leopard cat (Prionailurus bengalensis), an endangered species in South Korea, is a small feline widely distributed in Asia. Here, we investigated the diet of leopard cats in the inland areas of Korea by examining their fecal contents using vertebrate mitochondrial 12S rRNA gene sequencing and shotgun metagenomic sequencing. Shotgun metagenomic sequencing revealed that the feces were rich in DNA not only of vertebrates but also of arthropods and plants, but care should be taken when using shotgun metagenomic sequencing to identify vertebrates at low taxonomic levels (e.g., genus level), as it was often erroneous. Meanwhile, vertebrate mitochondrial 12S rRNA gene sequencing was found to be accurate in the genus-level identification, as the genera identified were consistent with the Korean fauna. We found that small mammals such as murids were their main prey. By using these two sequencing methods in combination, this study demonstrated that accurate information about the overall dietary content and vertebrate prey of leopard cats could be obtained. We expect that the continued community efforts to expand the genome database of wildlife, including vertebrates, will alleviate the problem of erroneous identification of prey at low taxonomic levels by shotgun metagenomic sequencing in the near future.

Molecular diet analysis enables detection of diatom and cyanobacteria DNA in the gut of Macoma balthica

Detritivores are essential to nutrient cycling, but are often neglected in trophic networks, due to difficulties with determining their diet. DNA analysis of gut contents shows promise of trophic link discrimination, but many unknown factors limit its usefulness. For example, DNA can be rapidly broken down, especially by digestion processes, and DNA provides only a snapshot of the gut contents at a specific time. Few studies have been performed on the length of time that prey DNA can be detected in consumer guts, and none so far using benthic detritivores. Eutrophication, along with climate change, is altering the phytoplankton communities in aquatic ecosystems, on which benthic detritivores in aphotic soft sediments depend. Nutrient-poor cyanobacteria blooms are increasing in frequency, duration, and magnitude in many water bodies, while nutrient-rich diatom spring blooms are shrinking in duration and magnitude, creating potential changes in diet of benthic detritivores. We performed an experiment to identify the taxonomy and quantify the abundance of phytoplankton DNA fragments on bivalve gut contents, and how long these fragments can be detected after consumption in the Baltic Sea clam Macoma balthica. Two common species of phytoplankton (the cyanobacteria Nodularia spumigena or the diatom Skeletonema marinoi) were fed to M. balthica from two regions (from the northern and southern Stockholm archipelago). After removing the food source, M. balthica gut contents were sampled every 24 hours for seven days to determine the number of 23S rRNA phytoplankton DNA copies and when the phytoplankton DNA could no longer be detected by quantitative PCR. We found no differences in diatom 18S rRNA gene fragments of the clams by region, but the southern clams showed significantly more cyanobacteria 16S rRNA gene fragments in their guts than the northern clams. Interestingly, the cyanobacteria and diatom DNA fragments were still detectable by qPCR in the guts of M. balthica one week after removal from its food source. However, DNA metabarcoding of the 23S rRNA phytoplankton gene found in the clam guts showed that added food (i.e. N. spumigena and S. marinoi) did not make up a majority of the detected diet. Our results suggest that these detritivorous clams therefore do not react as quickly as previously thought to fresh organic matter inputs, with other phytoplankton than large diatoms and cyanobacteria constituting the majority of their diet. This experiment demonstrates the viability of using molecular methods to determine feeding of detritivores, but further studies investigating how prey DNA signals can change over time in benthic detritivores will be needed before this method can be widely applicable to both models of ecological functions and conservation policy.

Using DNA metabarcoding and a novel canid-specific blocking oligonucleotide to investigate the composition of animal diets of raccoon dogs (Nyctereutes procyonoides) inhabiting the waterside area in Korea

The raccoon dog (Nyctereutes procyonoides) is known to be an opportunistic generalist who feeds on a wide variety of foods. Historically, their diet has been investigated by morphological observation of undigested remains in feces, requiring specialized knowledge such as osteology, zoology, and phytology. Here, we used DNA metabarcoding of vertebrate 12S rRNA gene and invertebrate 16S rRNA gene to investigate their fecal contents. Additionally, we developed a blocking oligonucleotide that specifically inhibits the amplification of the canid 12S rRNA gene. We confirmed that the blocking oligonucleotide selectively inhibit the amplification of raccoon dog’s DNA without significantly changing the composition of the preys’ DNA. We found that the main foods of raccoon dogs in our study area, the waterside of paddy fields in Korea, were fishes such as Cyprinidae and insects such as mole crickets, which makes sense given the Korean fauna and their well-known opportunistic feeding behaviors. As a method to conveniently and objectively investigate feeding habits of raccoon dogs, this study provided baseline information on DNA metabarcoding. By using DNA metabarcoding, it is expected that the diet habits and ecology of raccoon dogs will be better understood by future research.

Diet composition of Korean water deer (Hydropotes inermis argyropus) from the Han River Estuary Wetland in Korea using fecal DNA

To investigate the foraging preferences of the Korean water deer (Hydropotes inermis argyropus) from the Janghang wetlands located in Han River estuary in South Korea, we elucidated their diet compositions using DNA extracted from their fecal samples. Samples were collected from the Janghang wetlands from June 2013 to January 2015. For dietary analysis, we amplified the internal transcribed spacer 2 (ITS2) and the large subunit of ribulose-1,5-bisphosphate carboxylase (rbcL) regions, followed by cloning and sequencing. We were able to detect 20 families and 35 genera of plants from 77 fecal samples, showing that water deer in the Janghang wetlands mainly preferred Salicaceae (27.5%), Fabaceae (17.5%), and Poaceae (14.2%). In winter, the relative frequency of woody plant detection (81.3%) was much higher than that in other seasons. Our results suggest that Korean water deer in Janghang wetlands browse or graze depending on the locally-available plants species in their habitat as intermediate feeders. Non-invasive approaches using fecal samples could be applied to study food webs of various species and ecosystems for conservation and management of wildlife.

Quantitative molecular detection of larval Atlantic herring (Clupea harengus) in stomach contents of Atlantic mackerel (Scomber scombrus) marks regions of predation pressure

Mortality rates in the early life-history stages of fishes are generally high yet identifying the causes remain unclear. Faltering recruitment rates of Atlantic herring (Clupea harengus) in the Norwegian Sea indicate a need to identify which mortality factors influence larval herring survival. Previous research suggests that increased predation pressure by Atlantic mackerel (Scomber scombrus) may contribute to the disconnect between spawning stock biomass and recruitment. To quantify the contribution of predation pressure by Atlantic mackerel to herring larval mortality, two research cruises were conducted within a probable “hot spot” (67–72° N) for intensified mackerel predation based on particle drift simulations. Mackerel stomach contents were analysed for herring larvae content using droplet digital polymerase chain reaction (ddPCR) with a quantitative molecular detection assay specific for herring. The ddPCR results demonstrate clear predation by mackerel on herring larvae and also suggest that the alternative use of visual examination may give misleading results. Our results show that mackerel should be considered a potentially important predator on herring larvae. The quantitative molecular assay presented here shows great promise as an efficient and specific tool to correctly identify and quantify predation pressure on early life-history stages of fishes.

Fish as predators and prey: DNA-based assessment of their role in food webs

Fish are both consumers and prey, and as such part of a dynamic trophic network. Measuring how they are trophically linked, both directly and indirectly, to other species is vital to comprehend the mechanisms driving alterations in fish communities in space and time. Moreover, this knowledge also helps to understand how fish communities respond to environmental change and delivers important information for implementing management of fish stocks. DNA-based methods have significantly widened our ability to assess trophic interactions in both marine and freshwater systems and they possess a range of advantages over other approaches in diet analysis. In this review we provide an overview of different DNA-based methods that have been used to assess trophic interactions of fish as consumers and prey. We consider the practicalities and limitations, and emphasize critical aspects when analysing molecular derived trophic data. We exemplify how molecular techniques have been employed to unravel food web interactions involving fish as consumers and prey. In addition to the exciting opportunities DNA-based approaches offer, we identify current challenges and future prospects for assessing fish food webs where DNA-based approaches will play an important role.

Potential causative factors of noble scallop Chlamys nobilis mass mortality in Nan'ao Island, Shantou, China in 2017

Noble scallop Chlamys nobilis is an important marine bivalve that has been extensively cultured in the south coast of China since the 1980s. Unfortunately, since the late 1990s, the farmed scallops often suffered from regional mass mortality, which results in enormous economic losses to farmers and industries. In 2017, another mass mortality event occurred in Nan'ao Island, Shantou, China. In this study, the cause of C. nobilis mass mortality in 2017 was first investigated in the field, and then validated in a laboratory experiment. In the field, three sampling sites were selected according to the scallop mortality rate: Hunter Bay (90% mortality), Baisha Bay (67% mortality) and Longhai (6% mortality). Meanwhile, environmental parameters (temperature, salinity, DO, pH and chlorophyll a) of each site were also measured in situ. Then, water and scallop samples were collected randomly for the analysis of phytoplankton diversity and algal toxin activity using 18S rDNA and PP2A inhibition assay, respectively. In laboratory, healthy scallops were challenged with Karenia mikimotoi (1 × 10³ cells/mL) for 30 h. The field results showed that no significant difference in those environmental parameters existed among the three sites, but the relative abundance of K. mikimotoi in seawater and scallops' intestines in Hunter Bay and Baisha Bay was significantly higher than that in Longhai, and sick scallops contained significantly higher algal toxin activity than healthy ones. Laboratory results revealed that challenged scallops with K. mikimotoi showed significantly higher mortality rate and algal toxin activity than healthy ones, and low density of K. mikimotoi (1 × 10³ cells/mL) was sufficient to cause >50% scallops' mortality within 26 h. This study provides the first evidence that low K. mikimotoi cell density can cause massive mortality in C. nobilis, and provides useful information as guide to prevent scallop mass mortality in the future.

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.

DNA metabarcoding-based diet survey for the Eurasian otter (Lutra lutra): Development of a Eurasian otter-specific blocking oligonucleotide for 12S rRNA gene sequencing for vertebrates

The Eurasian otter (Lutra lutra) is an endangered species for which diet analyses are needed as part of its conservation efforts. Eurasian otters feed on vertebrates, such as fishes, and invertebrates, such as crustaceans, but their detailed taxonomies are not fully understood in part due to limited resolving power of traditional morphological identification methods. Here, we used high-throughput sequencing (HTS)-based DNA metabarcoding approaches to analyze diet profiles of Eurasian otters inhabiting a marshy estuary area in Korea. We investigated their diet profiles based on spraint sampling followed by DNA metabarcoding analyses targeting 12S rRNA gene region for vertebrates, 16S rRNA gene region for invertebrates, and cytochrome c oxidase 1 (COI) gene region for fishes. For the vertebrate analysis, a blocking oligonucleotide (OBS1) was designed to suppress amplification of DNA fragments derived from the otters. The 12S rRNA gene sequencing assay detected species belonging to fishes (95%) and amphibians (3.3%). Fishes detected by 12S rRNA gene sequencing included crucian carp (Carassius auratus), mullets (Mugil spp.), bluegill (Lepomis macrochirus), and northern snakehead (Channa argus), which were also detected by COI gene sequencing. Among invertebrates, mud flat crabs (Helicana spp.) and shrimps (Palaemon spp.) were abundant. The designed blocking oligonucleotide OBS1 effectively inhibited amplification of the otter’s DNA, with only up to 0.21% of vertebrate sequence reads assigned to the otter. This study demonstrated that HTS-based DNA metabarcoding methods were useful to provide in-depth information regarding diet profiles of the otters at our sampling site. By using HTS-based DNA metabarcoding approaches, future research will explore detailed taxonomies of their diets across locations and seasons.

Assessment of seasonal variation of diet composition in rodents using DNA barcoding and Real-Time PCR

The study of animal diet and feeding behaviour is a fundamental tool for the illustration of the ecological role of species in the ecosystem. However, size and quality of food intake samples make it hard for researchers to describe the diet composition of many small species. In our study, we exploited genomic tools for the analysis of the diet composition of the Savi’s pine vole (Microtus savii) using DNA barcoding and qPCR techniques for the identification of ingested plant species retrieved from stomach contents. In contrast with previous studies, we found that, despite being a fossorial species, the Savi’s pine vole is a selective feeder that undergoes intense superficial activity in search for food. In addition, our study shows that with a a priori knowledge of the candidate plant species included in animal diet, qPCR is a powerful tool to assess presence/absence, frequency of occurrence and electivity of ingested species. We conclude that this approach offers new opportunities to implement the analysis of food selection in small animals, thereby revealing a detailed picture of plant-animal interactions.

Using molecular diet analysis to inform invasive species management: A case study of introduced rats consuming endemic New Zealand frogs

The decline of amphibians has been of international concern for more than two decades, and the global spread of introduced fauna is a major factor in this decline. Conservation management decisions to implement control of introduced fauna are often based on diet studies. One of the most common metrics to report in diet studies is Frequency of Occurrence (FO), but this can be difficult to interpret, as it does not include a temporal perspective. Here, we examine the potential for FO data derived from molecular diet analysis to inform invasive species management, using invasive ship rats (Rattus rattus) and endemic frogs (Leiopelma spp.) in New Zealand as a case study. Only two endemic frog species persist on the mainland. One of these, Leiopelma archeyi, is Critically Endangered (IUCN 2017) and ranked as the world's most evolutionarily distinct and globally endangered amphibian (EDGE, 2018). Ship rat stomach contents were collected by kill-trapping and subjected to three methods of diet analysis (one morphological and two DNA-based). A new primer pair was developed targeting all anuran species that exhibits good coverage, high taxonomic resolution, and reasonable specificity. Incorporating a temporal parameter allowed us to calculate the minimum number of ingestion events per rat per night, providing a more intuitive metric than the more commonly reported FO. We are not aware of other DNA-based diet studies that have incorporated a temporal parameter into FO data. The usefulness of such a metric will depend on the study system, in particular the feeding ecology of the predator. Ship rats are consuming both species of native frogs present on mainland New Zealand, and this study provides the first detections of remains of these species in mammalian stomach contents.

Approaches to integrating genetic data into ecological networks

As molecular tools for assessing trophic interactions become common, research is increasingly focused on the construction of interaction networks. Here, we demonstrate three key methods for incorporating DNA data into network ecology and discuss analytical considerations using a model consisting of plants, insects, bats and their parasites from the Costa Rica dry forest. The simplest method involves the use of Sanger sequencing to acquire long sequences to validate or refine field identifications, for example of bats and their parasites, where one specimen yields one sequence and one identification. This method can be fully quantified and resolved and these data resemble traditional ecological networks. For more complex taxonomic identifications, we target multiple DNA loci, for example from a seed or fruit pulp sample in faeces. These networks are also well resolved but gene targets vary in resolution and quantification is difficult. Finally, for mixed templates such as faecal contents of insectivorous bats, we use DNA metabarcoding targeting two sequence lengths (157 and 407 bp) of one gene region and a MOTU, BLAST and BIN association approach to resolve nodes. This network type is complex to generate and analyse, and we discuss the implications of this type of resolution on network analysis. Using these data, we construct the first molecular-based network of networks containing 3,304 interactions between 762 nodes of eight trophic functions and involving parasitic, mutualistic and predatory interactions. We provide a comparison of the relative strengths and weaknesses of these data types in network ecology.

Herbivore-Diet Analysis Based on Illumina MiSeq Sequencing: The Potential Use of an ITS2-Barcoding Approach to Establish Qualitative and Quantitative Predictions of Diet Composition of Mongolian Sheep

DNA-barcoding approaches to estimate the diet compositions of grazing animals have received significant attention, and particularly when combined with next-generation sequencing, these techniques have substantially improved in recent years. In this study, the identity and species composition of plant material ingested by Mongolian sheep were estimated through the use of 350 bp ITS2 gene sequences of the vegetation found in fecal samples. Four diets were formulated using varying amounts of eight plant species that are common in the grasslands of northern China. Sixteen Mongolian sheep were taken from pastures and randomly assigned to four groups, and each group received one of four diets. Each sheep was randomly assigned to one of 16 confinement pens and fed its respective diet for 12 consecutive days. Fecal samples were removed from each pen from days 7-12, preserved, and composited for each pen. All herbage species included in the daily diets were detected in each fecal sample, with the exception of Phragmites australis. Moreover, 12 additional different plant species were retrieved from feces of the experimental sheep. The obtained data provided preliminary support for the use of the ITS2 barcode to determine which plants were consumed. Moreover, the proportions of the herbage DNA sequences recovered from sheep feces and those of the herbage masses in the daily diets did not completely match. These results indicate that the non-Gramineae DNA sequences amplified with ITS2 primers (including those of Chenopodium album, Artemisia scoparia, Artemisia tanacetifolia, and Medicago sativa) far exceeded those of the Gramineae species (including Leymus chinensis and Puccinellia distans), which constitute the largest share of the experimental diets. A significant positive correlation (Spearman's ρ = 0.376, P = 0.003) between the actual herbage mass proportions in the experimental diets and the herbage-DNA-sequence proportions provided sufficiently favorable support for the further investigation of DNA barcoding for the quantification of plants in feces. A significant regression coefficient was found between the relative DNA-sequence proportions of L. chinensis ( R² = 0.82, P < 0.0001), P. distans ( R² = 0.64, P = 0.0017), and C. album ( R² = 0.98, P < 0.0001) and their respective herbage mass proportions. The quantitative relationship can be expressed by the linear-regression equations y = 0.90 x - 0.22, y = 0.98 x - 0.03, and y = 5.00 x - 0.25, respectively. Thus, these results demonstrate that dietary-DNA-barcoding methods exhibited potential in providing valuable quantitative information regarding food-item components. However, it should be noted that this explorative data needs to be further improved by using additional genes and by creating a sophisticated reference database, thus enhancing both quality and accuracy of the obtained results.

Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data?

Advances in DNA sequencing technology have revolutionized the field of molecular analysis of trophic interactions, and it is now possible to recover counts of food DNA sequences from a wide range of dietary samples. But what do these counts mean? To obtain an accurate estimate of a consumer's diet should we work strictly with data sets summarizing frequency of occurrence of different food taxa, or is it possible to use relative number of sequences? Both approaches are applied to obtain semi-quantitative diet summaries, but occurrence data are often promoted as a more conservative and reliable option due to taxa-specific biases in recovery of sequences. We explore representative dietary metabarcoding data sets and point out that diet summaries based on occurrence data often overestimate the importance of food consumed in small quantities (potentially including low-level contaminants) and are sensitive to the count threshold used to define an occurrence. Our simulations indicate that using relative read abundance (RRA) information often provides a more accurate view of population-level diet even with moderate recovery biases incorporated; however, RRA summaries are sensitive to recovery biases impacting common diet taxa. Both approaches are more accurate when the mean number of food taxa in samples is small. The ideas presented here highlight the need to consider all sources of bias and to justify the methods used to interpret count data in dietary metabarcoding studies. We encourage researchers to continue addressing methodological challenges and acknowledge unanswered questions to help spur future investigations in this rapidly developing area of research.

DNA barcoding reveals seasonal shifts in diet and consumption of deep-sea fishes in wedge-tailed shearwaters

The foraging ecology of pelagic seabirds is difficult to characterize because of their large foraging areas. In the face of this difficulty, DNA metabarcoding may be a useful approach to analyze diet compositions and foraging behaviors. Using this approach, we investigated the diet composition and its seasonal variation of a common seabird species on the Ogasawara Islands, Japan: the wedge-tailed shearwater Ardenna pacifica. We collected fecal samples during the prebreeding (N = 73) and rearing (N = 96) periods. The diet composition of wedge-tailed shearwater was analyzed by Ion Torrent sequencing using two universal polymerase chain reaction primers for the 12S and 16S mitochondrial DNA regions that targeted vertebrates and mollusks, respectively. The results of a BLAST search of obtained sequences detected 31 and 1 vertebrate and mollusk taxa, respectively. The results of the diet composition analysis showed that wedge-tailed shearwaters frequently consumed deep-sea fishes throughout the sampling season, indicating the importance of these fishes as a stable food resource. However, there was a marked seasonal shift in diet, which may reflect seasonal changes in food resource availability and wedge-tailed shearwater foraging behavior. The collected data regarding the shearwater diet may be useful for in situ conservation efforts. Future research that combines DNA metabarcoding with other tools, such as data logging, may provide further insight into the foraging ecology of pelagic seabirds.

Main predators of insect pests: screening and evaluation through comprehensive indices

BACKGROUND

Predatory natural enemies play key functional roles in integrated pest management. However, the screening and evaluation of the main predators of insect pests has seldom been reported in the field. Here, we employed comprehensive indices for evaluating the predation of a common pest (Ectropis obliqua) by nine common spider species in Chinese tea plantations.

RESULTS

We established the relative dominance of the spider species and their phenological overlap with the pest species, and analyzed DNA from the nine spider species using targeted real-time quantitative polymerase chain reaction to identify the residual DNA of E. obliqua. The predation rates and predation numbers per predator were estimated by the positive rates of target fragments and the residual minimum number of E. obliqua in predators' guts, respectively. The results showed that only four spider species preyed on E. obliqua, and the order of potential of the spiders to control E. obliqua from greatest to smallest was Neoscona mellotteei, Xysticus ephippiatus, Evarcha albaria and Coleosoma octomaculatum by the Z-score method.

CONCLUSION

The orb-weaving spider N. mellotteei has the maximum potential as a biological control agent of E. obliqua in an integrated pest management strategy. An approach of screening and evaluating main predators of insect pests through comprehensive indices was preliminarily established. © 2017 Society of Chemical Industry.

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.

Simultaneous estimation of diet composition and calibration coefficients with fatty acid signature data

Knowledge of animal diets provides essential insights into their life history and ecology, although diet estimation is challenging and remains an active area of research. Quantitative fatty acid signature analysis (QFASA) has become a popular method of estimating diet composition, especially for marine species. A primary assumption of QFASA is that constants called calibration coefficients, which account for the differential metabolism of individual fatty acids, are known. In practice, however, calibration coefficients are not known, but rather have been estimated in feeding trials with captive animals of a limited number of model species. The impossibility of verifying the accuracy of feeding trial derived calibration coefficients to estimate the diets of wild animals is a foundational problem with QFASA that has generated considerable criticism. We present a new model that allows simultaneous estimation of diet composition and calibration coefficients based only on fatty acid signature samples from wild predators and potential prey. Our model performed almost flawlessly in four tests with constructed examples, estimating both diet proportions and calibration coefficients with essentially no error. We also applied the model to data from Chukchi Sea polar bears, obtaining diet estimates that were more diverse than estimates conditioned on feeding trial calibration coefficients. Our model avoids bias in diet estimates caused by conditioning on inaccurate calibration coefficients, invalidates the primary criticism of QFASA, eliminates the need to conduct feeding trials solely for diet estimation, and consequently expands the utility of fatty acid data to investigate aspects of ecology linked to animal diets.

The use of DNA barcodes in food web construction-terrestrial and aquatic ecologists unite!

By depicting who eats whom, food webs offer descriptions of how groupings in nature (typically species or populations) are linked to each other. For asking questions on how food webs are built and work, we need descriptions of food webs at different levels of resolution. DNA techniques provide opportunities for highly resolved webs. In this paper, we offer an exposé of how DNA-based techniques, and DNA barcodes in particular, have recently been used to construct food web structure in both terrestrial and aquatic systems. We highlight how such techniques can be applied to simultaneously improve the taxonomic resolution of the nodes of the web (i.e., the species), and the links between them (i.e., who eats whom). We end by proposing how DNA barcodes and DNA information may allow new approaches to the construction of larger interaction webs, and overcome some hurdles to achieving adequate sample size. Most importantly, we propose that the joint adoption and development of these techniques may serve to unite approaches to food web studies in aquatic and terrestrial systems-revealing the extent to which food webs in these environments are structured similarly to or differently from each other, and how they are linked by dispersal.

Molecular trophic markers in marine food webs and their potential use for coral ecology

Notable advances in ecological genomics have been driven by high-throughput sequencing technology and taxonomically broad sequence repositories that allow us to accurately assess species interactions with great taxonomic resolution. The use of DNA as a marker for ingested food is particularly relevant to address predator-prey interactions and disentangle complex marine food webs. DNA-based methods benefit from reductionist molecular approaches to address ecosystem scale processes, such as community structure and energy flow across trophic levels, among others. Here we review how molecular trophic markers have been used to better understand trophic interactions in the marine environment and their advantages and limitations. We focus on animal groups where research has been focused, such as marine mammals, seabirds, fishes, pelagic invertebrates and benthic invertebrates, and use case studies to illustrate how DNA-based methods unraveled food-web interactions. The potential of molecular trophic markers for disentangling the complex trophic ecology of corals is also discussed.

Omnivory of an Insular Lizard: Sources of Variation in the Diet of Podarcis lilfordi (Squamata, Lacertidae)

Through 17 years and from a sample of 7,790 faecal pellets and 26,346 prey items, we studied the diet of the Balearic lizard Podarcis lilfordi in Aire Island (Menorca, Balearic Islands, Spain). We analysed the diet in terms of prey frequencies, as well as by their volume and biomass contributions. The diet of the Balearic lizard was extremely variable through the years, months and areas under study. The dominance of small clumped prey, particularly ants, was confirmed. However, the main contribution by volume corresponded to beetles, with a relevant role for Diplopoda and terrestrial Isopoda during some months and at particular areas of the island. Several prey items were probably captured at the base of shrubs, under stones or inside rock crevices. Therefore, our estimations of electivity would only be reliable for epigeal and flying prey. The capacity of the Balearic lizard to include marine subsidies in its diet, such as coastal crustaceans, is noteworthy. Also, its consumption of carrion from carcasses of gulls and rabbits and leftovers from human visitors is remarkable. Juvenile conspecifics can also be a sporadic food resource, especially during the second half of summer, whereas the consumption of vegetal matter is constant for each whole year. The shifts of vegetal exploitation among areas of the island and months take place according to availability of different plant species at each area or during a given period. Thus, lizards are able to conduct a thorough monitoring of plant phenology, exploiting a large variety of plant species. Omnivory does not imply the indiscriminate inclusion of any edible food in its diet. Rather, the inclusion of several food items means the adoption of a wide range of foraging behaviours adapted to the exploitation of each food resource.

Maximizing dietary information retrievable from carcasses of Great Cormorants Phalacrocorax carbo using a combined morphological and molecular analytical approach

Avian carcasses can provide important information on the trophic ecology of birds. Usually, the number of carcasses available for examination is limited and therefore it is important to gain as much dietary information per specimen as possible. In piscivorous birds and raptors, the stomach has been the primary source of dietary information, whereas the gut (intestine) has so far been neglected as it usually contains only a few morphologically identifiable hard parts of prey. Molecular approaches have the potential to retrieve dietary information from the gut, although this has not yet been verified. As well as identifying the prey, it is important to estimate any secondary predation to avoid food web errors in dietary analyses. The assignment of accidentally consumed prey is notoriously difficult regardless of the prey identification approach used. In the present study, morphological and molecular analyses were, for the first time, combined to maximize the dietary information retrievable from the complete digestive tract of Great Cormorants Phalacrocorax carbo sinensis. Moreover, a novel approach based on predator-prey size ratios was applied to these piscivorous birds to minimize the number of samples that might contain secondarily predated prey. The stomach contents of the examined birds were found to provide the most dietary information when morphological and molecular analyses were used in combination. However, compared with the morphological approach, the molecular analysis increased the number of fish species detected by 39%. The molecular approach also permitted the identification of fish DNA in the Cormorant guts. Predator-prey size ratios derived from morphological analysis of fish hard parts can reduce the incidence of potential confounding influence of secondarily predated prey by 80%. Our findings demonstrate that a combination of morphological and molecular approaches maximizes the trophic information retrievable from bird carcasses.

A comparison of morphological and molecular diet analyses of predator scats

An understanding of a species’ diet is required to make sound conservation and management decisions. Traditionally, morphological analyses of undigested hard parts from food items remaining in scats have been used to assess diets. More recently, molecular analyses of scats have been used to identify plant and prey species’ DNA, but no studies have compared morphological and molecular diet analyses for large, terrestrial carnivores. We used molecular tools to determine the percentage of black bear and coyote scats that contained 3 common prey species (caribou, moose, and snowshoe hares) in Newfoundland and compared the results to a traditional morphological analysis. We found that a ranking of relative prey frequencies was consistent between the 2 methods, but molecular methods tended to detect prey species in a greater percentage of scats for all prey species. However, there were individual scats in which a prey species was detected by morphological methods only, and we provide evidence that molecular methods could result in false negatives if prey DNA is not uniformly distributed throughout a scat or as a result of PCR inconsistency. We also found that the per sample cost comparison between morphological and molecular analyses was dependent upon whether or not a molecular test was needed to identify scats to the predator species, the cost of developing molecular methods, and the number of samples being processed. We recommend that controlled feeding studies be performed to validate molecular methods and investigate the utility of molecular techniques to estimate the proportions of food items consumed.

DNA metabarcoding diet analysis for species with parapatric vs sympatric distribution: a case study on subterranean rodents

Closely related sympatric species commonly develop different ecological strategies to avoid competition. Ctenomys minutus and C. flamarioni are subterranean rodents parapatrically distributed in the southern Brazilian coastal plain, showing a narrow sympatric zone. To gain understanding on food preferences and possible competition for food resources, we evaluated their diet composition performing DNA metabarcoding analyzes of 67 C. minutus and 100 C. flamarioni scat samples, collected along the species geographical ranges. Thirteen plant families, mainly represented by Poaceae, Araliaceae, Asteraceae and Fabaceae, were identified in the diet of C. minutus. For C. flamarioni, 10 families were recovered, with a predominance of Poaceae, Araliaceae and Asteraceae. A significant correlation between diet composition and geographical distance was detected in C. minutus, whereas the diet of C. flamarioni was quite homogeneous throughout its geographical distribution. No significant differences were observed between males and females of each species. However, differences in diet composition between species were evident according to multivariate analysis. Our results suggest some level of diet partitioning between C. flamarioni and C. minutus in the sympatric region. While the first species is more specialized on few plant items, the second showed a more varied and heterogeneous diet pattern among individuals. These differences might have been developed to avoid competition in the region of co-occurrence. Resource availability in the environment also seems to influence food choices. Our data indicate that C. minutus and C. flamarioni are generalist species, but that some preference for Poaceae, Asteraceae and Araliaceae families can be suggested for both rodents.

Molecular detection of invertebrate prey in vertebrate diets: trophic ecology of Caribbean island lizards

Understanding community assembly and population dynamics frequently requires detailed knowledge of food web structure. For many consumers, obtaining precise information about diet composition has traditionally required sacrificing animals or other highly invasive procedures, generating tension between maintaining intact study populations and knowing what they eat. We developed 16S mitochondrial DNA sequencing methods to identify arthropods in the diets of generalist vertebrate predators without requiring a blocking primer. We demonstrate the utility of these methods for a common Caribbean lizard that has been intensively studied in the context of small island food webs: Anolis sagrei (a semi-arboreal 'trunk-ground' anole ecomorph). Novel PCR primers were identified in silico and tested in vitro. Illumina sequencing successfully characterized the arthropod component of 168 faecal DNA samples collected during three field trips spanning 12 months, revealing 217 molecular operational taxonomic units (mOTUs) from at least nine arthropod orders (including Araneae, Blattodea, Coleoptera, Hemiptera, Hymenoptera, Isoptera, Lepidoptera and Orthoptera). Three mOTUs (one beetle, one cockroach and one ant) were particularly frequent, occurring in ≥50% of samples, but the majority of mOTUs were infrequent (180, or 83%, occurred in ≤5% of samples). Species accumulation curves showed that dietary richness and composition were similar between size-dimorphic sexes; however, female lizards had greater per-sample dietary richness than males. Overall diet composition (but not richness) was significantly different across seasons, and we found more pronounced interindividual variation in December than in May. These methods will be generally useful in characterizing the diets of diverse insectivorous vertebrates.

Dining local: the microbial diet of a snail that grazes microbial communities is geographically structured

Achatinella mustelina is a critically endangered tree snail that subsists entirely by grazing microbes from leaf surfaces of native trees. Little is known about the fundamental aspects of these microbe assemblages: not taxonomic composition, how this varies with host plant or location, nor whether snails selectively consume microbes. To address these questions, we collected 102 snail faecal samples as a proxy for diet, and 102 matched-leaf samples from four locations. We used Illumina amplicon sequencing to determine bacterial and fungal community composition. Microbial community structure was significantly distinct between snail faeces and leaf samples, but the same microbes occurred in both. We conclude that snails are not 'picky' eaters at the microbial level, but graze the surface of whatever plant they are on. In a second experiment, the gut was dissected from non-endangered native tree snails in the same family as Achatinella to confirm that faecal samples reflect gut contents. Over 60% of fungal reads were shared between faeces, gut and leaf samples. Overall, location, sample type (faeces or leaf) and host plant identity all significantly explained the community composition and variation among samples. Understanding the microbial ecology of microbes grazed by tree snails enables effective management when conservation requires captive breeding or field relocation.

Molecular detection of trophic interactions: emerging trends, distinct advantages, significant considerations and conservation applications

The emerging field of ecological genomics contains several broad research areas. Comparative genomic and conservation genetic analyses are providing great insight into adaptive processes, species bottlenecks, population dynamics and areas of conservation priority. Now the same technological advances in high-throughput sequencing, coupled with taxonomically broad sequence repositories, are providing greater resolution and fundamentally new insights into functional ecology. In particular, we now have the capacity in some systems to rapidly identify thousands of species-level interactions using non-invasive methods based on the detection of trace DNA. This represents a powerful tool for conservation biology, for example allowing the identification of species with particularly inflexible niches and the investigation of food-webs or interaction networks with unusual or vulnerable dynamics. As they develop, these analyses will no doubt provide significant advances in the field of restoration ecology and the identification of appropriate locations for species reintroduction, as well as highlighting species at ecological risk. Here, I describe emerging patterns that have come from the various initial model systems, the advantages and limitations of the technique and key areas where these methods may significantly advance our empirical and applied conservation practices.

Utilizing stomach content and faecal DNA analysis techniques to assess the feeding behaviour of largemouth bass Micropterus salmoides and bluegill Lepomis macrochirus

In this study, the feeding behaviour of the non-native invasive predatory fishes largemouth bass Micropterus salmoides and bluegill Lepomis macrochirus was studied in the Ezura River, a northern tributary of Lake Biwa, Japan. Prey composition was estimated based on visual examination of stomach contents and faecal DNA analysis to determine feeding habits of these predatory fishes. Stomach content analysis showed that native fishes (e.g. ayu Plecoglossus altivelis and gobies Rhinogobius spp.) and shrimps (e.g. Palaemon paucidens) were the major prey items for M. salmoides, while snails, larval Chironomidae and submerged macrophytes were the dominant prey items of L. macrochirus. Micropterus salmoides tended to select larger fish in the case of crucian carp Carassius spp., but smaller fishes in the case of P. altivelis and Rhinogobius spp. Faecal DNA analyses revealed prey compositions similar to those identified in predator stomach contents, and identified additional prey species not detected in stomach content inspection. This study demonstrated that both stomach content inspection and DNA-based analysis bear several inherent shortcomings and advantages. The former method is straightforward, although identification of species can be inaccurate or impossible, whereas the latter method allows for accurate species identification, but cannot distinguish prey size or stage. Hence, integration of morphology-based and DNA-based methods can provide more reliable estimates of foraging habits of predatory fishes.

Improving accuracy of DNA diet estimates using food tissue control materials and an evaluation of proxies for digestion bias

Ecologists are increasingly interested in quantifying consumer diets based on food DNA in dietary samples and high-throughput sequencing of marker genes. It is tempting to assume that food DNA sequence proportions recovered from diet samples are representative of consumer's diet proportions, despite the fact that captive feeding studies do not support that assumption. Here, we examine the idea of sequencing control materials of known composition along with dietary samples in order to correct for technical biases introduced during amplicon sequencing and biological biases such as variable gene copy number. Using the Ion Torrent PGM(©) , we sequenced prey DNA amplified from scats of captive harbour seals (Phoca vitulina) fed a constant diet including three fish species in known proportions. Alongside, we sequenced a prey tissue mix matching the seals' diet to generate tissue correction factors (TCFs). TCFs improved the diet estimates (based on sequence proportions) for all species and reduced the average estimate error from 28 ± 15% (uncorrected) to 14 ± 9% (TCF-corrected). The experimental design also allowed us to infer the magnitude of prey-specific digestion biases and calculate digestion correction factors (DCFs). The DCFs were compared with possible proxies for differential digestion (e.g. fish protein%, fish lipid%) revealing a strong relationship between the DCFs and percent lipid of the fish prey, suggesting prey-specific corrections based on lipid content would produce accurate diet estimates in this study system. These findings demonstrate the value of parallel sequencing of food tissue mixtures in diet studies and offer new directions for future research in quantitative DNA diet analysis.

Using next-generation sequencing to analyse the diet of a highly endangered land snail (Powelliphanta augusta) feeding on endemic earthworms

Predation is often difficult to observe or quantify for species that are rare, very small, aquatic or nocturnal. The assessment of such species' diet can be conducted using molecular methods that target prey DNA remaining in predators' guts and faeces. These techniques do not require high taxonomic expertise, are applicable to soft-bodied prey and allow for identification at the species level. However, for generalist predators, the presence of mixed prey DNA in guts and faeces can be a major impediment as it requires development of specific primers for each potential prey species for standard (Sanger) sequencing. Therefore, next generation sequencing methods have recently been applied to such situations. In this study, we used 454-pyrosequencing to analyse the diet of Powelliphantaaugusta, a carnivorous landsnail endemic to New Zealand and critically endangered after most of its natural habitat has been lost to opencast mining. This species was suspected to feed mainly on earthworms. Although earthworm tissue was not detectable in snail faeces, earthworm DNA was still present in sufficient quantity to conduct molecular analyses. Based on faecal samples collected from 46 landsnails, our analysis provided a complete map of the earthworm-based diet of P. augusta. Predated species appear to be earthworms that live in the leaf litter or earthworms that come to the soil surface at night to feed on the leaf litter. This indicates that P. augusta may not be selective and probably predates any earthworm encountered in the leaf litter. These findings are crucial for selecting future translocation areas for this highly endangered species. The molecular diet analysis protocol used here is particularly appropriate to study the diet of generalist predators that feed on liquid or soft-bodied prey. Because it is non-harmful and non-disturbing for the studied animals, it is also applicable to any species of conservation interest.

Stranded dolphin stomach contents represent the free-ranging population's diet

Diet is a fundamental aspect of animal ecology. Cetacean prey species are generally identified by examining stomach contents of stranded individuals. Critical uncertainty in these studies is whether samples from stranded animals are representative of the diet of free-ranging animals. Over two summers, we collected faecal and gastric samples from healthy free-ranging individuals of an extensively studied bottlenose dolphin population. These samples were analysed by molecular prey detection and these data compared with stomach contents data derived from stranded dolphins from the same population collected over 22 years. There was a remarkable consistency in the prey species composition and relative amounts between the two datasets. The conclusions of past stomach contents studies regarding dolphin habitat associations, prey selection and proposed foraging mechanisms are supported by molecular data from live animals and the combined dataset. This is the first explicit test of the validity of stomach contents analysis for accurate population-scale diet determination of an inshore cetacean.

Quantifying sequence proportions in a DNA-based diet study using Ion Torrent amplicon sequencing: which counts count?

A goal of many environmental DNA barcoding studies is to infer quantitative information about relative abundances of different taxa based on sequence read proportions generated by high-throughput sequencing. However, potential biases associated with this approach are only beginning to be examined. We sequenced DNA amplified from faeces (scats) of captive harbour seals (Phoca vitulina) to investigate whether sequence counts could be used to quantify the seals' diet. Seals were fed fish in fixed proportions, a chordate-specific mitochondrial 16S marker was amplified from scat DNA and amplicons sequenced using an Ion Torrent PGM™. For a given set of bioinformatic parameters, there was generally low variability between scat samples in proportions of prey species sequences recovered. However, proportions varied substantially depending on sequencing direction, level of quality filtering (due to differences in sequence quality between species) and minimum read length considered. Short primer tags used to identify individual samples also influenced species proportions. In addition, there were complex interactions between factors; for example, the effect of quality filtering was influenced by the primer tag and sequencing direction. Resequencing of a subset of samples revealed some, but not all, biases were consistent between runs. Less stringent data filtering (based on quality scores or read length) generally produced more consistent proportional data, but overall proportions of sequences were very different than dietary mass proportions, indicating additional technical or biological biases are present. Our findings highlight that quantitative interpretations of sequence proportions generated via high-throughput sequencing will require careful experimental design and thoughtful data analysis.

Sliding window analyses for optimal selection of mini-barcodes, and application to 454-pyrosequencing for specimen identification from degraded DNA

DNA barcoding remains a challenge when applied to diet analyses, ancient DNA studies, environmental DNA samples and, more generally, in any cases where DNA samples have not been adequately preserved. Because the size of the commonly used barcoding marker (COI) is over 600 base pairs (bp), amplification fails when the DNA molecule is degraded into smaller fragments. However, relevant information for specimen identification may not be evenly distributed along the barcoding region, and a shorter target can be sufficient for identification purposes. This study proposes a new, widely applicable, method to compare the performance of all potential 'mini-barcodes' for a given molecular marker and to objectively select the shortest and most informative one. Our method is based on a sliding window analysis implemented in the new R package SPIDER (Species IDentity and Evolution in R). This method is applicable to any taxon and any molecular marker. Here, it was tested on earthworm DNA that had been degraded through digestion by carnivorous landsnails. A 100 bp region of 16 S rDNA was selected as the shortest informative fragment (mini-barcode) required for accurate specimen identification. Corresponding primers were designed and used to amplify degraded earthworm (prey) DNA from 46 landsnail (predator) faeces using 454-pyrosequencing. This led to the detection of 18 earthworm species in the diet of the snail. We encourage molecular ecologists to use this method to objectively select the most informative region of the gene they aim to amplify from degraded DNA. The method and tools provided here, can be particularly useful (1) when dealing with degraded DNA for which only small fragments can be amplified, (2) for cases where no consensus has yet been reached on the appropriate barcode gene, or (3) to allow direct analysis of short reads derived from massively parallel sequencing without the need for bioinformatic consolidation.

Who is eating what: diet assessment using next generation sequencing

The analysis of food webs and their dynamics facilitates understanding of the mechanistic processes behind community ecology and ecosystem functions. Having accurate techniques for determining dietary ranges and components is critical for this endeavour. While visual analyses and early molecular approaches are highly labour intensive and often lack resolution, recent DNA-based approaches potentially provide more accurate methods for dietary studies. A suite of approaches have been used based on the identification of consumed species by characterization of DNA present in gut or faecal samples. In one approach, a standardized DNA region (DNA barcode) is PCR amplified, amplicons are sequenced and then compared to a reference database for identification. Initially, this involved sequencing clones from PCR products, and studies were limited in scale because of the costs and effort required. The recent development of next generation sequencing (NGS) has made this approach much more powerful, by allowing the direct characterization of dozens of samples with several thousand sequences per PCR product, and has the potential to reveal many consumed species simultaneously (DNA metabarcoding). Continual improvement of NGS technologies, on-going decreases in costs and current massive expansion of reference databases make this approach promising. Here we review the power and pitfalls of NGS diet methods. We present the critical factors to take into account when choosing or designing a suitable barcode. Then, we consider both technical and analytical aspects of NGS diet studies. Finally, we discuss the validation of data accuracy including the viability of producing quantitative data.

DNA-based faecal dietary analysis: a comparison of qPCR and high throughput sequencing approaches

The genetic analysis of faecal material represents a relatively non-invasive way to study animal diet and has been widely adopted in ecological research. Due to the heterogeneous nature of faecal material the primary obstacle, common to all genetic approaches, is a means to dissect the constituent DNA sequences. Traditionally, bacterial cloning of PCR amplified products was employed; less common has been the use of species-specific quantitative PCR (qPCR) assays. Currently, with the advent of High-Throughput Sequencing (HTS) technologies and indexed primers it has become possible to conduct genetic audits of faecal material to a much greater depth than previously possible. To date, no studies have systematically compared the estimates obtained by HTS with that of qPCR. What are the relative strengths and weaknesses of each technique and how quantitative are deep-sequencing approaches that employ universal primers? Using the locally threatened Little Penguin (Eudyptula minor) as a model organism, it is shown here that both qPCR and HTS techniques are highly correlated and produce strikingly similar quantitative estimates of fish DNA in faecal material, with no statistical difference. By designing four species-specific fish qPCR assays and comparing the data to the same four fish in the HTS data it was possible to directly compare the strengths and weaknesses of both techniques. To obtain reproducible quantitative data one of the key, and often overlooked, steps common to both approaches is ensuring that efficient DNA isolation methods are employed and that extracts are free of inhibitors. Taken together, the methodology chosen for long-term faecal monitoring programs is largely dependent on the complexity of the prey species present and the level of accuracy that is desired. Importantly, these methods should not be thought of as mutually exclusive, as the use of both HTS and qPCR in tandem will generate datasets with the highest fidelity.