The problem of omnivory: A synthesis on omnivory and DNA metabarcoding

Networking nutrients: How nutrition determines the structure of ecological networks

Nutrients can shape ecological interactions but remain poorly integrated into ecological networks. Concepts such as nutrient-specific foraging nevertheless have the potential to expose the mechanisms structuring complex ecological systems. Nutrients also present an opportunity to predict dynamic processes, such as interaction rewiring and extinction cascades, and increase the accuracy of network analyses. Here, we propose the concept of nutritional networks. By integrating nutritional data into ecological networks, we envisage significant advances to our understanding of ecological processes from individual to ecosystem scales. We show that networks can be constructed with nutritional data to illuminate how nutrients structure ecological interactions in natural systems through an empirical example. Throughout, we identify fundamental ecological hypotheses that can be explored in a nutritional network context, alongside methods for resolving those networks. Nutrients influence the structure and complexity of ecological networks through mechanistic processes and concepts including nutritional niche differentiation, functional responses, landscape diversity, ecological invasions and ecosystem robustness. Future research on ecological networks should consider nutrients when investigating the drivers of network structure and function.

A global synthesis of predation on bivalves

Predation is a dominant structuring force in ecological communities. In aquatic environments, predation on bivalves has long been an important focal interaction for ecological study because bivalves have central roles as ecosystem engineers, basal components of food webs, and commercial commodities. Studies of bivalves are common, not only because of bivalves' central roles, but also due to the relative ease of studying predatory effects on this taxonomic group. To understand patterns in the interactions of bivalves and their predators we synthesised data from 52 years of peer-reviewed studies on bivalve predation. Using a systematic search, we compiled 1334 studies from 75 countries, comprising 61 bivalve families (N = 2259), dominated by Mytilidae (29% of bivalves), Veneridae (14%), Ostreidae (8%), Unionidae (7%), and Dreissenidae and Tellinidae (6% each). A total of 2036 predators were studied, with crustaceans the most studied predator group (34% of predators), followed by fishes (24%), molluscs (17%), echinoderms (10%) and birds (6%). The majority of studies (86%) were conducted in marine systems, in part driven by the high commercial value of marine bivalves. Studies in freshwater ecosystems were dominated by non-native bivalves and non-native predator species, which probably reflects the important role of biological invasions affecting freshwater biodiversity. In fact, while 81% of the studied marine bivalve species were native, only 50% of the freshwater species were native to the system. In terms of approach, most studies used predation trials, visual analysis of digested contents and exclusion experiments to assess the effects of predation. These studies reflect that many factors influence bivalve predation depending on the species studied, including (i) species traits (e.g. behaviour, morphology, defence mechanisms), (ii) other biotic interactions (e.g. presence of competitors, parasites or diseases), and (iii) environmental context (e.g. temperature, current velocity, beach exposure, habitat complexity). There is a lack of research on the effects of bivalve predation at the population and community and ecosystem levels (only 7% and 0.5% of studies respectively examined impacts at these levels). At the population level, the available studies demonstrate that predation can decrease bivalve density through consumption or the reduction of recruitment. At the community and ecosystem level, predation can trigger effects that cascade through trophic levels or effects that alter the ecological functions bivalves perform. Given the conservation and commercial importance of many bivalve species, studies of predation should be pursued in the context of global change, particularly climate change, acidification and biological invasions.

DNA metabarcoding, diversity partitioning and null models reveal mechanisms of seasonal trophic specialization in a Mediterranean warbler

Optimal Foraging Theory (OFT) predicts that a population's trophic niche expansion should occur in periods of food scarcity as individuals begin to opportunistically exploit sub-optimal food items. However, the Niche Variation Hypothesis (NVH) posits that niche widening may result from increased among-individual differentiation due to food partitioning to avoid competition. We tested these hypotheses through a DNA metabarcoding study of the Sardinian Warbler (Curruca melanocephala) diet over a year. We used null models and the decomposition of beta diversity on among-individual dietary differentiation to infer the mechanisms driving the population's niche variation. Warblers fed frequently on berries, with a peak in late summer and, to a lesser extent, in autumn. Their diet also included a wide range of arthropods, with their prevalence varying among seasons. Consistent with OFT, the population's niche width was narrower in spring/summer when the population was strongly specialized in berries. In winter, the population's niche expanded, possibly reflecting seasonal declines in food abundance. As predicted by NVH, among-individual differentiation tended to be higher in winter, but this was mainly due to increased differences in dietary richness rather than to the partitioning of resources. Overall, our results suggest that within-individual niche does not increase in lean periods, and instead, individuals adopt either a more opportunistic or more specialized foraging strategy. Increased competition in periods of scarcity may help explain such patterns, but instead of showing increased food partitioning as expected from NVH, it may reflect OFT mechanisms on individuals with differential competitive ability to access better food resources.

DNA metabarcoding the diet of Podarcis lizards endemic to the Balearic Islands

Dietary studies are essential to unravel the functioning of ecosystems and ultimately to understand biodiversity. This task, which at first may seem simple, becomes especially complex in those cases of omnivorous species with highly variable diets. In this regard, the emergence of next-generation DNA sequencing methodologies represents a powerful tool to address the problem. Here we implement a high-throughput metabarcoding strategy based on the analysis of four molecular markers aimed at sequencing both mitochondrial (animal prey) and chloroplast (diet plants) genome fragments from fecal samples of two lizard species endemic to the Balearic Archipelago (Podarcis lilfordi and P. pityusensis) obtained through non-invasive methods. The results allowed for the characterization of their diets with a high degree of taxonomic detail and have contributed a large number of new trophic records. The reported diets are based mainly on the consumption of arthropods, mollusks and plants from a diversity of taxonomic orders, as well as carrion and marine subsidies. Our analyses also reveal inter- and intra-specific differences both in terms of seasonality and geographical distribution of the sampled lizard populations. These molecular findings provide new insights into the trophic interactions of these threatened endemic lizards in their unique and isolated ecosystems.

Otterly delicious: Spatiotemporal variation in the diet of a recovering population of Eurasian otters (Lutra lutra) revealed through DNA metabarcoding and morphological analysis of prey remains

Eurasian otters are apex predators of freshwater ecosystems and a recovering species across much of their European range; investigating the dietary variation of this predator over time and space, therefore, provides opportunities to identify changes in freshwater trophic interactions and factors influencing the conservation of otter populations. Here we sampled feces from 300 dead otters across England and Wales between 2007 and 2016, conducting both morphological analyses of prey remains and dietary DNA metabarcoding. Comparison of these methods showed that greater taxonomic resolution and breadth could be achieved using DNA metabarcoding but combining data from both methodologies gave the most comprehensive dietary description. All otter demographics exploited a broad range of taxa and variation likely reflected changes in prey distributions and availability across the landscape. This study provides novel insights into the trophic generalism and adaptability of otters across Britain, which is likely to have aided their recent population recovery, and may increase their resilience to future environmental changes.

Multi-marker DNA metabarcoding reveals spatial and sexual variation in the diet of a scarce woodland bird

Avian diet can be affected by site-specific variables, such as habitat, as well as intrinsic factors such as sex. This can lead to dietary niche separation, which reduces competition between individuals, as well as impacting how well avian species can adapt to environmental variation. Estimating dietary niche separation is challenging, due largely to difficulties in accurately identifying food taxa consumed. Consequently, there is limited knowledge of the diets of woodland bird species, many of which are undergoing serious population declines. Here, we show the effectiveness of multi-marker fecal metabarcoding to provide in-depth dietary analysis of a declining passerine in the UK, the Hawfinch (Coccothraustes coccothraustes). We collected fecal samples from (n = 262) UK Hawfinches prior to, and during, the breeding seasons in 2016-2019. We detected 49 and 90 plant and invertebrate taxa, respectively. We found Hawfinch diet varied spatially, as well as between sexes, indicating broad dietary plasticity and the ability of Hawfinches to utilize multiple resources within their foraging environments.

Dietary Variation Is Driven by Landscape Heterogeneity in an Insular Omnivorous Endemic Lizard, Revealed by DNA Metabarcoding

Living on islands entails numerous challenges for animals, among which population density approaching the carrying capacity of trophic resources stands out. To overcome this limitation, many insular lizards can supplement their insectivorous diet with increasing portions of plant material. The Madeira wall lizard, Teira dugesii, is a medium-sized lacertid, endemic to the Madeira and Selvagens archipelagos. As common in this family, adults are sexually dimorphic with males being bigger than females. Previous dietary studies on morphological scatology identified a higher proportion of plant over animal prey items, changing according to the location and sex. Here, we used DNA metabarcoding to examine the diet of this lizard species quantifying it at a higher taxonomical resolution and enhancing the detection of soft-body prey that often go undetected in morphology-based studies. In a sample of 151 faecal samples from eight populations including different habitats and altitudes in Madeira, we identified 289 prey items belonging to eight animal and three plant Classes, encompassing 58 distinct orders and 140 families. Of these, 63 were identified up to the species level. The results support a strong trend towards herbivory in this species with plants representing almost 74% of the diet occurrences in contrast to the 26% of animal prey. Remarkably, the plant fraction of the diet remained stable across localities but varied with size and mass in males. As males grew bigger and heavier, they significantly increased their plant matter intake. Likely, larger bodies and abdomens allowed allocating longer and more complex digestive tracts harbouring intestinal flora to better decompose plant organic compounds. This allowed heavier animals to have a richer diet regime. However, diet richness and composition were not affected by either sex or size, while the locality had a significant effect on both diet components likely in response to local variation in prey availability. By including an increasing plant fraction into a primarily insectivorous diet, this insular lizard has not only enlarged its trophic niche but is also able to exploit more efficiently the highly variable resources provided by insular environments.

Using eDNA metabarcoding to understand the effect of fire on the diet of small mammals in a woodland ecosystem

Food acquisition is a fundamental process that drives animal distribution and abundance, influencing how species respond to changing environments. Disturbances such as fire create significant shifts in available dietary resources, yet, for many species, we lack basic information about what they eat, let alone how they respond to a changing resource base. In order to create effective management strategies, faunal conservation in flammable landscapes requires a greater understanding of what animals eat and how this change following a fire. What animals eat in postfire environments has received little attention due to the time‐consuming methodologies and low‐resolution identification of food taxa. Recently, molecular techniques have been developed to identify food DNA in scats, making it possible to identify animal diets with enhanced resolution. The primary aim of this study was to utilize eDNA metabarcoding to obtain an improved understanding of the diet of three native Australian small mammal species: yellow‐footed antechinus (Antechinus flavipes), heath mouse (Pseudomys shortridgei), and bush rat (Rattus fuscipes). Specifically, we sought to understand the difference in the overall diet of the three species and how diet changed over time after fire. Yellow‐footed antechinus diets mostly consisted of moths, and plants belonging to myrtles and legume families while bush rats consumed legumes, myrtles, rushes, and beetles. Heath mouse diet was dominated by rushes. All three species shifted their diets over time after fire, with most pronounced shifts in the bush rats and least for heath mice. Identifying critical food resources for native animals will allow conservation managers to consider the effect of fire management actions on these resources and help conserve the species that use them.

The complex epistemological challenge of data curation in dietary metabarcoding: Comment on "The precautionary principle and dietary DNA metabarcoding: Commonly used abundance thresholds change ecological interpretation" by Littleford-Colquhoun et al. (2022)

In their article, Littleford-Colquhoun et al. (2022) advise against using arbitrary relative read abundance (RRA) thresholds (i.e., minimum sequence copy thresholds) for removing low-abundance sequences since they can increase false negative rates in dietary DNA metabarcoding data sets. The main criticisms presented against these widespread methods are that they (i) are arbitrary, often existing as standard values or defined based on researcher-selected delineations, (ii) are subjective, varying between studies and contexts, and, most problematically, (iii) result in the exclusion of true positives, particularly rarely consumed taxa, to the detriment of ecological insight. We commend the authors for presenting a refreshing and timely perspective on this often neglected topic, which is certainly in need of greater discussion following over a decade of significant advances in dietary metabarcoding. In this complex epistemological problem of false positives versus false negatives, we feel that several of the points raised deserve additional discussion. We address these aspects below, including measured approaches to data filtration and consistent representation of RRAs, and we welcome any further discourse to solidify or refute the concepts therein.

DNA metabarcoding reveals that coyotes in New York City consume wide variety of native prey species and human food

Carnivores are currently colonizing cities where they were previously absent. These urban environments are novel ecosystems characterized by habitat degradation and fragmentation, availability of human food, and different prey assemblages than surrounding areas. Coyotes (Canis latrans) established a breeding population in New York City (NYC) over the last few decades, but their ecology within NYC is poorly understood. In this study, we used non-invasive scat sampling and DNA metabarcoding to profile vertebrate, invertebrate, and plant dietary items with the goal to compare the diets of urban coyotes to those inhabiting non-urban areas. We found that both urban and non-urban coyotes consumed a variety of plants and animals as well as human food. Raccoons (Procyon lotor) were an important food item for coyotes within and outside NYC. In contrast, white-tailed deer (Odocoileus virginianus) were mainly eaten by coyotes inhabiting non-urban areas. Domestic chicken (Gallus gallus) was the human food item found in most scats from both urban and non-urban coyotes. Domestic cats (Felis catus) were consumed by urban coyotes but were detected in only a small proportion of the scats (<5%), which differs markedly from high rates of cat depredation in some other cities. In addition, we compared our genetic metabarcoding analysis to a morphological analysis of the same scat samples. We found that the detection similarity between the two methods was low and it varied depending on the type of diet item.

The predator problem and PCR primers in molecular dietary analysis: swamped or silenced; depth or breadth?

Dietary metabarcoding has vastly improved our ability to analyse the diets of animals, but it is hampered by a plethora of technical limitations including potentially reduced data output due to the disproportionate amplification of the DNA of the focal predator, here termed 'the predator problem'. We review the various methods commonly used to overcome this problem, from deeper sequencing to exclusion of predator DNA during PCR, and how they may interfere with increasingly common multi-predator-taxon studies. We suggest that multi-primer approaches with an emphasis on achieving both depth and breadth of prey detections may overcome the issue to some extent, although multi-taxon studies require further consideration, as highlighted by an empirical example. We also review several alternative methods for reducing the prevalence of predator DNA that are conceptually promising but require additional empirical examination. The predator problem is a key constraint on molecular dietary analyses but, through this synthesis, we hope to guide researchers in overcoming this in an effective and pragmatic way.

Dietary preferences and feeding strategies of Colombian highland woolly monkeys

Primates are very selective in the foods they include in their diets with foraging strategies that respond to spatial and temporal changes in resource availability, distribution and quality. Colombian woolly monkeys (Lagothrix lagotricha lugens), one of the largest primate species in the Americas, feed mainly on fruits, but they also eat a high percentage of arthropods. This differs from closely related Atelid species that supplement their diet with leaves. In an 11 month study, we investigated the foraging strategies of this endemic monkey and assessed how resource availability affects dietary selection. Using behavioural, phenological, arthropod sampling and metabarcoding methods, we recorded respectively foraging time, forest productivity, arthropod availability in the forest and arthropod consumption. Scat samples and capturing canopy substrates (i.e. moss, bromeliads, aerial insects) were used for assigning arthropod taxonomy. The most important resource in the diet was fruits (54%), followed by arthropods (28%). Resource availability predicted feeding time for arthropods but not for fruits. Further, there was a positive relationship between feeding time on fruits and arthropods, suggesting that eating both resources during the same periods might work as an optimal strategy to maximize nutrient intake. Woolly monkeys preferred and avoided some fruit and arthropod items available in their home range, choosing a wide variety of arthropods. Geometrid moths (Lepidoptera) were the most important and consistent insects eaten over time. We found no differences in the type of arthropods adults and juveniles ate, but adults invested more time foraging for this resource, especially in moss. Although woolly monkeys are generalist foragers, they do not select their food items randomly or opportunistically.

Interspecific coprophagia by wild red foxes: DNA metabarcoding reveals a potentially widespread form of commensalism among animals

Vertebrate animals are known to consume other species' faeces, yet the role of such coprophagy in species dynamics remains unknown, not least due to the methodological challenges of documenting it. In a large‐scale metabarcoding study of red fox and pine marten scats, we document a high occurrence of domestic dog DNA in red fox scats and investigate if it can be attributed to interspecific coprophagia. We tested whether experimental artifacts or other sources of DNA could account for dog DNA, regressed dog occurrence in the diet of fox against that of the fox’s main prey, short‐tailed field voles, and consider whether predation or scavenging could explain the presence of dog DNA. Additionally, we determined the calorific value of dog faeces through calorimetric explosion. The high occurrence of dog DNA in the diet of fox, the timing of its increase, and the negative relationship between dog and the fox's main prey, point to dog faeces as the source of DNA in fox scats. Dog faeces being highly calorific, we found that foxes, but not pine martens, regularly exploit them, seemingly as an alternative resource to fluctuating prey. Scattered accounts from the literature may suggest that interspecific coprophagia is a potentially frequent and widespread form of interaction among vertebrates. However, further work should address its prevalence in other systems and the implications for ecological communities. Tools such as metabarcoding offer a way forward.

Successes and limitations of quantitative diet metabarcoding in a small, herbivorous mammal

DNA metabarcoding is widely used to determine wild animal diets, however whether this technique provides accurate, quantitative measurements is still under debate. To test our ability to accurately estimate abundance of dietary items using metabarcoding, we fed wild caught desert woodrats (Neotoma lepida) diets comprised of constant amounts of juniper (Juniperus osteosperma, 15%) and varying amounts of creosote (Larrea tridentata, 1-60%), or cactus (Opuntia sp., 0-100%), and commercial chow (0-85%). Using metabarcoding, we compared the representation of items in the original diet samples to that in the fecal samples to test the sensitivity and accuracy of diet metabarcoding, the performance of different bioinformatic pipelines, and our ability to correct sequence counts. Metabarcoding, using standard trnL primers, detected creosote, juniper, and chow. Different pipelines for assigning taxonomy performed similarly. While creosote was detectable at dietary proportions as low as 1%, we failed to detect cactus in most samples, likely due to a primer mismatch. Creosote read counts increased as its proportion in the diet increased, and we could differentiate when creosote was a minor and major component of the diet. However, we found that estimates of juniper and creosote varied. Using previously suggested methods to correct these errors did not improve accuracy estimates of creosote, but did reduce error for juniper and chow. Our results indicate that metabarcoding can provide quantitative information on dietary composition, but may be limited. We suggest that researchers use caution in quantitatively interpreting diet metabarcoding results unless they first experimentally determine the extent of possible biases.

Specific gut bacterial responses to natural diets of tropical birds

The composition of gut bacterial communities is strongly influenced by the host diet in many animal taxa. For birds, the effect of diet on the microbiomes has been documented through diet manipulation studies. However, for wild birds, most studies have drawn on literature-based information to decipher the dietary effects, thereby, overlooking individual variation in dietary intake. Here we examine how naturally consumed diets influence the composition of the crop and cloacal microbiomes of twenty-one tropical bird species, using visual and metabarcoding-based identification of consumed diets and bacterial 16S rRNA microbiome sequencing. We show that diet intakes vary markedly between individuals of the same species and that literature-based dietary guilds grossly underestimate intraspecific diet variability. Furthermore, despite an effect of literature-based dietary guild assignment of host taxa, the composition of natural diets does not align with crop and cloacal microbiome similarity. However, host-taxon specific gut bacterial lineages are positively correlated with specific diet items, indicating that certain microbes associate with different diet components in specific avian hosts. Consequently, microbiome composition is not congruent with the overall consumed diet composition of species, but specific components of a consumed diet lead to host-specific effects on gut bacterial taxa.

DNA metabarcoding reveals introduced species predominate in the diet of a threatened endemic omnivore, Telfair's skink (Leiolopisma telfairii)

Introduced species can exert disproportionately negative effects on island ecosystems, but their potential role as food for native consumers is poorly studied. Telfair's skinks are endemic omnivores living on Round Island, Mauritius, a globally significant site of biodiversity conservation. We aimed to determine the dietary diversity and key trophic interactions of Telfair's skinks, whether introduced species are frequently consumed, and if diet composition changes seasonally between male and female skinks. We used DNA metabarcoding of skink fecal samples to identify animals (COI) and plants (ITS2) consumed by skinks. There were 389 dietary presence counts belonging to 77 dietary taxa found across the 73 Telfair's skink fecal samples. Introduced taxa were cumulatively consumed more frequently than other categories, accounting for 49.4% of all detections, compared to cryptogenic (20.6%), native (20.6%), and endemic taxa (9.5%). The most frequently consumed introduced species was the ant, Pheidole megacephala, present in 40% of samples. Blue latan palm, Latania loddigesii, was the most frequently consumed endemic species, present in 33% of samples but was only detected in the dry season, when fruits are produced. We found a strong seasonal difference in diet composition explained by the presence of certain plant species solely or primarily in one season and a marked increase in the consumption of animal prey in the dry season. Male and female skinks consumed several taxa at different frequencies. These results present a valuable perspective on the role of introduced species in the trophic network of their invaded ecosystem. Both native and introduced species provide nutritional resources for skinks, and this may have management implications in the context of species conservation and island restoration.

DNA-metabarcoding reveals the importance of gelatinous zooplankton in the diet of Pandalus borealis, a keystone species in the Arctic

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly molecular approaches such as DNA-metabarcoding are used in studying trophic relations, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA-metabarcoding with universal primers for cytochrome c oxidase I (COI), to study the diet composition of the Northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40-60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on SCA. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underscores the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

Handling of targeted amplicon sequencing data focusing on index hopping and demultiplexing using a nested metabarcoding approach in ecology

High-throughput sequencing platforms are increasingly being used for targeted amplicon sequencing because they enable cost-effective sequencing of large sample sets. For meaningful interpretation of targeted amplicon sequencing data and comparison between studies, it is critical that bioinformatic analyses do not introduce artefacts and rely on detailed protocols to ensure that all methods are properly performed and documented. The analysis of large sample sets and the use of predefined indexes create challenges, such as adjusting the sequencing depth across samples and taking sequencing errors or index hopping into account. However, the potential biases these factors introduce to high-throughput amplicon sequencing data sets and how they may be overcome have rarely been addressed. On the example of a nested metabarcoding analysis of 1920 carabid beetle regurgitates to assess plant feeding, we investigated: (i) the variation in sequencing depth of individually tagged samples and the effect of library preparation on the data output; (ii) the influence of sequencing errors within index regions and its consequences for demultiplexing; and (iii) the effect of index hopping. Our results demonstrate that despite library quantification, large variation in read counts and sequencing depth occurred among samples and that the sequencing error rate in bioinformatic software is essential for accurate adapter/primer trimming and demultiplexing. Moreover, setting an index hopping threshold to avoid incorrect assignment of samples is highly recommended.

Multiple species-specific molecular markers using nanofluidic array as a tool to detect prey DNA from carnivore scats

Large carnivore feeding ecology plays a crucial role for management and conservation for predators and their prey. One of the keys to this kind of research is to identify the species composition in the predator diet, for example, prey determination from scat content. DNA-based methods applied to detect prey in predators' scats are viable alternatives to traditional macroscopic approaches, showing an increased reliability and higher prey detection rate. Here, we developed a molecular method for prey species identification in wolf (Canis lupus) scats using multiple species-specific marker loci on the cytochrome b gene for 18 target species. The final panel consisted of 80 assays, with a minimum of four markers per target species, and that amplified specifically when using a high-throughput Nanofluidic array technology (Fluidigm Inc.). As a practical example, we applied the method to identify target prey species DNA in 80 wolf scats collected in Sweden. Depending on the number of amplifying markers required to obtain a positive species call in a scat, the success in determining at least one prey species from the scats ranged from 44% to 92%. Although we highlight the need to evaluate the optimal number of markers for sensitive target species detection, the developed method is a fast and cost-efficient tool for prey identification in wolf scats and it also has the potential to be further developed and applied to other areas and large carnivores as well.