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

Increasing marine trophic web knowledge through DNA analyses of fish stomach content: a step towards an ecosystem-based approach to fisheries research

Multispecies and ecosystem models, which are key for the implementation of ecosystem-based approaches to fisheries management, require extensive data on the trophic interactions between marine organisms, including changes over time. DNA metabarcoding, by allowing the simultaneous taxonomic identification of the community present in hundreds of samples, could be used for speeding up large-scale stomach content data collection. Yet, for DNA metabarcoding to be routinely implemented, technical challenges should be addressed, such as the potentially complicated sampling logistics, the detection of a high proportion of predator DNA, and the inability to provide reliable abundance estimations. Here, we present a DNA metabarcoding assay developed to examine the diet of five commercially important fish, which can be feasibly incorporated into routinary samplings. The method is devised to speed up the analysis process by avoiding the stomach dissection and content extraction steps, while preventing the amplification of predator DNA by using blocking primers. Tested in mock samples and in real stomach samples, the method has proven effective and shows great effectiveness discerning diet variations due to predator ecology or prey availability. Additionally, by applying our protocol to mackerel stomachs previously analyzed by visual inspection, we showcase how DNA metabarcoding could complement visually based data by detecting overlooked prey by the visual approach. We finally discuss how DNA metabarcoding-based data can contribute to trophic data collection. Our work reinforces the potential of DNA metabarcoding for the study and monitoring of fish trophic interactions and provides a basis for its incorporation into routine monitoring programs, which will be critical for the implementation of ecosystem-based approaches to fisheries management.

Using eDNA to understand predator-prey interactions influenced by invasive species

Invasive predatory species may alter population dynamic processes of their prey and impact biological communities and ecosystem processes. Revealing biotic interactions, however, including the relationship between predator and prey, is a difficult task, in particular for species that are hard to monitor. Here, we present a case study that documents the utility of environmental DNA analysis (eDNA) to assess predator-prey interactions between two invasive fishes (Lepomis gibbosus, Pseudorasbora parva) and two potential amphibian prey species, (Triturus cristatus, Pelobates fuscus). We used species-specific TaqMan assays for quantitative assessment of eDNA concentrations from water samples collected from 89 sites across 31 ponds during three consecutive months from a local amphibian hotspot in Germany. We found a negative relationship between eDNA concentrations of the predators (fishes) and prey (amphibians) using Monte-Carlo tests. Our study highlights the potential of eDNA application to reveal predator-prey interactions and confirms the hypothesis that the observed local declines of amphibian species may be at least partly caused by recently introduced invasive fishes. Our findings have important consequences for local conservation management and highlight the usefulness of eDNA approaches to assess ecological interactions and guide targeted conservation action.

DNA metabarcoding highlights cyanobacteria as the main source of primary production in a pelagic food web model

Models that estimate rates of energy flow in complex food webs often fail to account for species-specific prey selectivity of diverse consumer guilds. While DNA metabarcoding is increasingly used for dietary studies, methodological biases have limited its application for food web modeling. Here, we used data from dietary metabarcoding studies of zooplankton to calculate prey selectivity indices and assess energy fluxes in a pelagic resource-consumer network. We show that food web dynamics are influenced by prey selectivity and temporal match-mismatch in growth cycles and that cyanobacteria are the main source of primary production in the investigated coastal pelagic food web. The latter challenges the common assumption that cyanobacteria are not supporting food web productivity, a result that is increasingly relevant as global warming promotes cyanobacteria dominance. While this study provides a method for how DNA metabarcoding can be used to quantify energy fluxes in a marine food web, the approach presented here can easily be extended to other ecosystems.

Food webs reveal coexistence mechanisms and community organization in carnivores

Globally, massive carnivore guild extirpations have led to trophic downgrading and compromised ecosystem services. However, the complexity of multi-carnivore food webs complicates accurate identification of species interactions and community organization. Here, we used fecal DNA metabarcoding to investigate three communities that together encompass eight large- and meso-carnivore species and their 44 prey taxa of the Qinghai-Tibet Plateau (QTP), one of the last places on Earth that still harbors intact carnivore assemblages. Quantitative food-web analyses revealed pronounced interspecific variations in the carnivores' prey compositions and dietary partitioning both between and within guilds. Additionally, body masses of the carnivores and their prey exhibited consistent hump-shaped correlations across communities. Overall, differences in prey diversity, size category, and proportional utilization among the carnivore species result in trophic niche segregation that likely promotes carnivore coexistence in the harsh QTP environment. Network structure analyses detected significant modularity in all food webs but nestedness in only one. Furthermore, network characterization identified pikas (Ochotona spp.), bharal (Pseudois nayaur), and domestic yak (Bos grunniens) as potential keystone prey across the areas. Our results paint a holistic and detailed picture of the QTP carnivore assemblages' trophic networks and demonstrate that the combined use of the molecular dietary approach and network analysis can generate structural insights into carnivore coexistence and can identify functionally important species in complex communities. Such knowledge can help safeguard carnivore guild integrity and enhance community resilience to environmental perturbations in the sensitive QTP ecosystems.

Domestic Cattle in a National Park Restricting the Sika Deer Due to Diet Overlap

Managers need to know the extent of the conflict between livestock and wild animals. Although many studies have reported the conflict between livestock and wild animals, few have checked the extent of the conflict. Cattle raising in the Northeast Tiger and Leopard National Park is considered one of the main driving forces behind the restricted distribution of sika deer. To understand whether foraging competition is contributing to avoidance patterns between sika deer and cattle, we investigated their feeding habits using DNA barcoding and high-throughput sequencing. Our study shows that although cattle are grazers in the traditional division of herbivores, their diet shifted to a predominance of dicotyledonous woody plants, and this diet shift resulted in a high degree of dietary overlap between sika deer and cattle. Moreover, compared to sika deer, cattle diets are more diverse at the species level with a wider ecological niche. Our results confirm that overlapping dietary niches and the superior competitive abilities of cattle contribute to the restricted distribution of the sika deer, which has critical implications for the conservation of their predators. Our study suggests that cattle grazing should be prohibited in the Park and effective measures should be taken for the benefit of sika deer.

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.

Fishing for fish environmental DNA: Ecological applications, methodological considerations, surveying designs, and ways forward

Vast global declines of freshwater and marine fish diversity and population abundance pose serious threats to both ecosystem sustainability and human livelihoods. Environmental DNA (eDNA)-based biomonitoring provides robust, efficient, and cost-effective assessment of species occurrences and population trends in diverse aquatic environments. Thus, it holds great potential for improving conventional surveillance frameworks to facilitate fish conservation and fisheries management. However, the many technical considerations and rapid developments underway in the eDNA arena can overwhelm researchers and practitioners new to the field. Here, we systematically analysed 416 fish eDNA studies to summarize research trends in terms of investigated targets, research aims, and study systems, and reviewed the applications, rationales, methodological considerations, and limitations of eDNA methods with an emphasis on fish and fisheries research. We highlighted how eDNA technology may advance our knowledge of fish behaviour, species distributions, population genetics, community structures, and ecological interactions. We also synthesized the current knowledge of several important methodological concerns, including the qualitative and quantitative power eDNA has to recover fish biodiversity and abundance, and the spatial and temporal representations of eDNA with respect to its sources. To facilitate ecological applications implementing fish eDNA techniques, recent literature was summarized to generate guidelines for effective sampling in lentic, lotic, and marine habitats. Finally, we identified current gaps and limitations, and pointed out newly emerging research avenues for fish eDNA. As methodological optimization and standardization improve, eDNA technology should revolutionize fish monitoring and promote biodiversity conservation and fisheries management that transcends geographic and temporal boundaries.

Temporal and spatial changes in benthic invertebrate trophic networks along a taxonomic richness gradient

Species interactions underlie most ecosystem functions and are important for understanding ecosystem changes. Representing one type of species interaction, trophic networks were constructed from biodiversity monitoring data and known trophic links to assess how ecosystems have changed over time. The Baltic Sea is subject to many anthropogenic pressures, and low species diversity makes it an ideal candidate for determining how pressures change food webs. In this study, we used benthic monitoring data for 20 years (1980-1989 and 2010-2019) from the Swedish coast of the Baltic Sea and Skagerrak to investigate changes in benthic invertebrate trophic interactions. We constructed food webs and calculated fundamental food web metrics evaluating network horizontal and vertical diversity, as well as stability that were compared over space and time. Our results show that the west coast of Sweden (Skagerrak) suffered a reduction in benthic invertebrate biodiversity by 32% between the 1980s and 2010s, and that the number of links, generality of predators, and vulnerability of prey have been significantly reduced. The other basins (Bothnian Sea, Baltic Proper, and Bornholm Basin) do not show any significant changes in species richness or consistent significant trends in any food web metrics investigated, demonstrating resilience at a lower species diversity. The decreased complexity of the Skagerrak food webs indicates vulnerability to further perturbations and pressures should be limited as much as possible to ensure continued ecosystem functions.

EcoDiet: A hierarchical Bayesian model to combine stomach, biotracer, and literature data into diet matrix estimation

Although quantifying trophic interactions is a critical path to understanding and forecasting ecosystem functioning, fitting trophic models to field data remains challenging. It requires flexible statistical tools to combine different sources of information from the literature and fieldwork samples. We present EcoDiet, a hierarchical Bayesian modeling framework to simultaneously estimate food-web topology and diet composition of all consumers in the food web, by combining (1) a priori knowledge from the literature on both food-web topology and diet proportions; (2) stomach content analyses, with frequencies of prey occurrence used as the primary source of data to update the prior knowledge on the topological food-web structure; (3) and biotracers data through a mixing model (MM). Inferences are derived in a Bayesian probabilistic rationale that provides a formal way to incorporate prior information and quantifies uncertainty around both the topological structure of the food web and the dietary proportions. EcoDiet was implemented as an open-source R package, providing a user-friendly interface to execute the model, as well as examples and guidelines to familiarize with its use. We used simulated data to demonstrate the benefits of EcoDiet and how the framework can improve inferences on diet matrix by comparison with classical network MM. We applied EcoDiet to the Celtic Sea ecosystem, and showed how combining multiple data types within an integrated approach provides a more robust and holistic picture of the food-web topology and diet matrices than the literature or classical MM approach alone. EcoDiet has the potential to become a reference method for building diet matrices as a preliminary step of ecosystem modeling and to improve our understanding of prey-predator interactions.

Unveiling the complexity and ecological function of aquatic macrophyte-animal networks in coastal ecosystems

Network theory offers innovative tools to explore the complex ecological mechanisms regulating species associations and interactions. Although interest in ecological networks has grown steadily during the last two decades, the application of network approaches has been unequally distributed across different study systems: while some kinds of interactions (e.g. plant-pollinator and host-parasite) have been extensively investigated, others remain relatively unexplored. Among the latter, aquatic macrophyte-animal associations in coastal environments have been largely neglected, despite their major role in littoral ecosystems. The ubiquity of macrophyte systems, their accessibility and multi-faceted ecological, economical and societal importance make macrophyte-animal systems an ideal subject for ecological network science. In fact, macrophyte-animal networks offer an aquatic counterpart to terrestrial plant-animal networks. In this review, we show how the application of network analysis to aquatic macrophyte-animal associations has the potential to broaden our understanding of how coastal ecosystems function. Network analysis can also provide a key to understanding how such ecosystems will respond to on-going and future threats from anthropogenic disturbance and environmental change. For this, we: (i) identify key issues that have limited the application of network theory and modelling to aquatic animal-macrophyte associations; (ii) illustrate through examples based on empirical data how network analysis can offer new insights on the complexity and functioning of coastal ecosystems; and (iii) provide suggestions for how to design future studies and establish this new research line into network ecology.

Using DNA metabarcoding to decipher the diet plant component of mammals from the Eastern Mediterranean region

Longevity of species populations depends largely on interactions among animals and plants in an ecosystem. Predation and seed dispersal are among the most important interactions necessary for species conservation and persistence, and diet analysis is a prerequisite tool to evaluate these interactions. Understanding these processes is crucial for identifying conservation targets and for executing efficient reforestation and ecological restoration. In this study, we applied a scat DNA metabarcoding technique using the P6-loop of the trnL (UAA) chloroplastic marker to describe the seasonal plant diet composition of 15 mammal species from a highly biodiverse Lebanese forest in the Eastern Mediterranean. We also recovered plant seeds, when present, from the scats for identification. The mammal species belong to 10 families from 5 different orders. More than 133 plant species from 54 plant families were detected and identified. Species from the Rosaceae, Poaceae, Apiaceae, Fabaceae, Fagaceae and Berberidaceae families were consumed by the majority of the mammals and should be taken into consideration in future reforestation and conservation projects. Our results showed that the DNA metabarcoding approach provides a promising method for tracking the dietary plant components of a wide diversity of mammals, yielding key insights into plant-animal interactions inside Lebanon’s forests.

Prey partitioning and livestock consumption in the world's richest large carnivore assemblage

Large mammalian carnivores have undergone catastrophic declines during the Anthropocene across the world. Despite their pivotal roles as apex predators in food webs and ecosystem dynamics, few detailed dietary datasets of large carnivores exist, prohibiting deep understanding of their coexistence and persistence in human-dominated landscapes. Here, we present fine-scaled, quantitative trophic interactions among sympatric carnivores from three assemblages in the Mountains of Southwest China, a global biodiversity hotspot harboring the world's richest large-carnivore diversity, derived from DNA metabarcoding of 1,097 fecal samples. These assemblages comprise a large-carnivore guild ranging from zero to five species along with two mesocarnivore species. We constructed predator-prey food webs for each assemblage and identified 95 vertebrate prey taxa and 260 feeding interactions in sum. Each carnivore species consumed 6-39 prey taxa, and dietary diversity decreased with increased carnivore body mass across guilds. Dietary partitioning was more evident between large-carnivore and mesocarnivore guilds, yet different large carnivores showed divergent proportional utilization of different-sized prey correlating with their own body masses. Large carnivores particularly selected livestock in Tibetan-dominated regions, where the indigenous people show high tolerance toward wild predators. Our results suggest that dietary niche partitioning and livestock subsidies facilitate large-carnivore sympatry and persistence and have key implications for sustainable conservation promoting human-carnivore coexistence.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

Altered structure of bat-prey interaction networks in logged tropical forests revealed by metabarcoding

Habitat degradation is pervasive across the tropics and is particularly acute in Southeast Asia, with major implications for biodiversity. Much research has addressed the impact of degradation on species diversity; however, little is known about how ecological interactions are altered, including those that constitute important ecosystem functions such as consumption of herbivores. To examine how rainforest degradation alters trophic interaction networks, we applied DNA metabarcoding to construct interaction networks linking forest-dwelling insectivorous bat species and their prey, comparing old-growth forest and forest degraded by logging in Sabah, Borneo. Individual bats in logged rainforest consumed a lower richness of prey than those in old-growth forest. As a result, interaction networks in logged forests had a less nested structure. These network structures were associated with reduced network redundancy and thus increased vulnerability to perturbations in logged forests. Our results show how ecological interactions change between old-growth and logged forests, with potentially negative implications for ecosystem function and network stability.

Cascading effects of moth outbreaks on subarctic soil food webs

The increasing severity and frequency of natural disturbances requires a better understanding of their effects on all compartments of biodiversity. In Northern Fennoscandia, recent large-scale moth outbreaks have led to an abrupt change in plant communities from birch forests dominated by dwarf shrubs to grass-dominated systems. However, the indirect effects on the belowground compartment remained unclear. Here, we combined eDNA surveys of multiple trophic groups with network analyses to demonstrate that moth defoliation has far-reaching consequences on soil food webs. Following this disturbance, diversity and relative abundance of certain trophic groups declined (e.g., ectomycorrhizal fungi), while many others expanded (e.g., bacterivores and omnivores) making soil food webs more diverse and structurally different. Overall, the direct and indirect consequences of moth outbreaks increased belowground diversity at different trophic levels. Our results highlight that a holistic view of ecosystems improves our understanding of cascading effects of major disturbances on soil food webs.

DNA metabarcoding reveals trophic niche diversity of micro and mesozooplankton species

Alternative pathways of energy transfer guarantee the functionality and productivity in marine food webs that experience strong seasonality. Nevertheless, the complexity of zooplankton interactions is rarely considered in trophic studies because of the lack of detailed information about feeding interactions in nature. In this study, we used DNA metabarcoding to highlight the diversity of trophic niches in a wide range of micro- and mesozooplankton, including ciliates, rotifers, cladocerans, copepods and their prey, by sequencing 16- and 18S rRNA genes. Our study demonstrates that the zooplankton trophic niche partitioning goes beyond both phylogeny and size and reinforces the importance of diversity in resource use for stabilizing food web efficiency by allowing for several different pathways of energy transfer. We further highlight that small, rarely studied zooplankton (rotifers and ciliates) fill an important role in the Baltic Sea pelagic primary production pathways and the potential of ciliates, rotifers and crustaceans in the utilization of filamentous and picocyanobacteria within the pelagic food web. The approach used in this study is a suitable entry point to ecosystem-wide food web modelling considering species-specific resource use of key consumers.

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.

Resolving Food-Web Structure

Food webs are a major focus and organizing theme of ecology, but the data used to assemble them are deficient. Early debates over food-web data focused on taxonomic resolution and completeness, lack of which had produced spurious inferences. Recent data are widely believed to be much better and are used extensively in theoretical and meta-analytic research on network ecology. Confidence in these data rests on the assumptions ( a) that empiricists correctly identified consumers and their foods and ( b) that sampling methods were adequate to detect a near-comprehensive fraction of the trophic interactions between species. Abundant evidence indicates that these assumptions are often invalid, suggesting that most topological food-web data may remain unreliable for inferences about network structure and underlying ecological and evolutionary processes. Morphologically cryptic species are ubiquitous across taxa and regions, and many trophic interactions routinely evade detection by conventional methods. Molecular methods have diagnosed the severity of these problems and are a necessary part of the cure.

Can the carbon and nitrogen isotope values of offspring be used as a proxy for their mother's diet? Using foetal physiology to interpret bulk tissue and amino acid δ15N values

The measurement of bulk tissue nitrogen (δ¹⁵N) and carbon isotope values (δ¹³C) chronologically along biologically inert tissues sampled from offspring can provide a longitudinal record of their mothers' foraging habits. This study tested the important assumption that mother-offspring stable isotope values are positively and linearly correlated. In addition, any change in the mother-offspring bulk tissues and individual amino acids that occurred during gestation was investigated. Whiskers sampled from southern elephant seal pups (Mirounga leonina) and temporally overlapping whiskers from their mothers were analyzed. This included n = 1895 chronologically subsampled whisker segments for bulk tissue δ¹⁵N and δ¹³C in total and n = 20 whisker segments for amino acid δ¹⁵N values, sampled from recently weaned pups (n = 17), juvenile southern elephant seals (SES) < 2 years old (n = 23) and adult female SES (n = 17), which included nine mother-offspring pairs. In contrast to previous studies, the mother-offspring pairs were not in isotopic equilibrium or linearly correlated during gestation: the Δ¹⁵N and Δ¹³C mother-offspring offsets increased by 0.8 and 1.2‰, respectively, during gestation. The foetal bulk δ¹⁵N values were 1.7 ± 0.5‰ (0.9-2.7‰) higher than mothers' δ¹⁵N values before birth, while the foetal δ¹³C increased by ~1.7‰ during gestation and were 1.0 ± 0.5‰ (0.0-1.9‰) higher than their mothers' δ¹³C at the end of pregnancy. The mother-offspring serine and glycine Δ¹⁵N differed by ~4.3‰, while the foetal alanine δ¹⁵N values were 1.4‰ lower than that of their mothers during the third trimester of pregnancy. The observed mother-offspring δ¹⁵N differences are likely explained by shuttling of glutamate-glutamine and glycine-serine amongst skeletal muscle, liver, placenta and foetal tissue. Foetal development relies primarily on remobilized endogenous maternal proteinaceous sources. Researchers should consider foetal physiology when using offspring bulk tissue isotope values as biomarkers for the mother's isotopic composition as part of monitoring programmes.

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.

Wing morphology predicts individual niche specialization in Pteronotus mesoamericanus (Mammalia: Chiroptera)

Morphological variation between individuals can increase niche segregation and decrease intraspecific competition when heterogeneous individuals explore their environment in different ways. Among bat species, wing shape correlates with flight maneuverability and habitat use, with species that possess broader wings typically foraging in more cluttered habitats. However, few studies have investigated the role of morphological variation in bats for niche partitioning at the individual level. To determine the relationship between wing shape and diet, we studied a population of the insectivorous bat species Pteronotus mesoamericanus in the dry forest of Costa Rica. Individual diet was resolved using DNA metabarcoding, and bat wing shape was assessed using geometric morphometric analysis. Inter-individual variation in wing shape showed a significant relationship with both dietary dissimilarity based on Bray-Curtis estimates, and nestedness derived from an ecological network. Individual bats with broader and more rounded wings were found to feed on a greater diversity of arthropods (less nested) in comparison to individuals with triangular and pointed wings (more nested). We conclude that individual variation in bat wing morphology can impact foraging efficiency leading to the observed overall patterns of diet specialization and differentiation within the population.

Towards a unified study of multiple stressors: divisions and common goals across research disciplines

Anthropogenic environmental changes, or 'stressors', increasingly threaten biodiversity and ecosystem functioning worldwide. Multiple-stressor research is a rapidly expanding field of science that seeks to understand and ultimately predict the interactions between stressors. Reviews and meta-analyses of the primary scientific literature have largely been specific to either freshwater, marine or terrestrial ecology, or ecotoxicology. In this cross-disciplinary study, we review the state of knowledge within and among these disciplines to highlight commonality and division in multiple-stressor research. Our review goes beyond a description of previous research by using quantitative bibliometric analysis to identify the division between disciplines and link previously disconnected research communities. Towards a unified research framework, we discuss the shared goal of increased realism through both ecological and temporal complexity, with the overarching aim of improving predictive power. In a rapidly changing world, advancing our understanding of the cumulative ecological impacts of multiple stressors is critical for biodiversity conservation and ecosystem management. Identifying and overcoming the barriers to interdisciplinary knowledge exchange is necessary in rising to this challenge. Division between ecosystem types and disciplines is largely a human creation. Species and stressors cross these borders and so should the scientists who study them.

Methylmercury exposure in wildlife: A review of the ecological and physiological processes affecting contaminant concentrations and their interpretation

Exposure to methylmercury (MeHg) can result in detrimental health effects in wildlife. With advances in ecological indicators and analytical techniques for measurement of MeHg in a variety of tissues, numerous processes have been identified that can influence MeHg concentrations in wildlife. This review presents a synthesis of theoretical principals and applied information for measuring MeHg exposure and interpreting MeHg concentrations in wildlife. Mercury concentrations in wildlife are the net result of ecological processes influencing dietary exposure combined with physiological processes that regulate assimilation, transformation, and elimination. Therefore, consideration of both physiological and ecological processes should be integrated when formulating biomonitoring strategies. Ecological indicators, particularly stable isotopes of carbon, nitrogen, and sulfur, compound-specific stable isotopes, and fatty acids, can be effective tools to evaluate dietary MeHg exposure. Animal species differ in their physiological capacity for MeHg elimination, and animal tissues can be inert or physiologically active, act as sites of storage, transformation, or excretion of MeHg, and vary in the timing of MeHg exposure they represent. Biological influences such as age, sex, maternal transfer, and growth or fasting are also relevant for interpretation of tissue MeHg concentrations. Wildlife tissues that represent current or near-term bioaccumulation and in which MeHg is the predominant mercury species (such as blood and eggs) are most effective for biomonitoring ecosystems and understanding landscape drivers of MeHg exposure. Further research is suggested to critically evaluate the use of keratinized external tissues to measure MeHg bioaccumulation, particularly for less-well studied wildlife such as reptiles and terrestrial mammals. Suggested methods are provided to effectively use wildlife for quantifying patterns and drivers of MeHg bioaccumulation over time and space, as well as for assessing the potential risk and toxicological effects of MeHg on wildlife.

Characterizing ecological interaction networks to support risk assessment in classical biological control of weeds

A key element in weed biological control is the selection of a biological control agent that minimizes the risks of non-target attack and indirect effects on the recipient community. Network ecology is a promising approach that could help decipher tritrophic interactions in both the native and the invaded ranges, to complement quarantine-based host-specificity tests and gain insights on potential interactions of biological control agents. This review highlights practical questions addressed by networks, including 1) biological control agent selection, based on specialization indices, 2) risk assessment of biological control agent release into a novel environment, via particular patterns of association such as apparent competition between agent(s) and native herbivore(s), 3) network comparisons through structural metrics, 4) potential of network modelling and 5) limits of network construction methods.

Threats from the air: Damselfly predation on diverse prey taxa

To understand the diversity and strength of predation in natural communities, researchers must quantify the total amount of prey species in the diet of predators. Metabarcoding approaches have allowed widespread characterization of predator diets with high taxonomic resolution. To determine the wider impacts of predators, researchers should combine DNA techniques with estimates of population size of predators using mark-release-recapture (MRR) methods, and with accurate metrics of food consumption by individuals. Herein, we estimate the scale of predation exerted by four damselfly species on diverse prey taxa within a well-defined 12-ha study area, resolving the prey species of individual damselflies, to what extent the diets of predatory species overlap, and which fraction of the main prey populations are consumed. We identify the taxonomic composition of diets using DNA metabarcoding and quantify damselfly population sizes by MRR. We also use predator-specific estimates of consumption rates, and independent data on prey emergence rates to estimate the collective predation pressure summed over all prey taxa and specific to their main prey (non-biting midges or chironomids) of the four damselfly species. The four damselfly species collectively consumed a prey mass equivalent to roughly 870 (95% CL 410-1,800) g, over 2 months. Each individual consumed 29%-66% (95% CL 9.4-123) of its body weight during its relatively short life span (2.1-4.7 days; 95% CL 0.74-7.9) in the focal population. This predation pressure was widely distributed across the local invertebrate prey community, including 4 classes, 19 orders and c. 140 genera. Different predator species showed extensive overlap in diets, with an average of 30% of prey shared by at least two predator species. Of the available prey individuals in the widely consumed family Chironomidae, only a relatively small proportion (0.76%; 95% CL 0.35%-1.61%) were consumed. Our synthesis of population sizes, per-capita consumption rates and taxonomic distribution of diets identifies damselflies as a comparatively minor predator group of aerial insects. As the next step, we should add estimates of predation by larger odonate species, and experimental removal of odonates, thereby establishing the full impact of odonate predation on prey communities.

Consuming alternative prey does not influence the DNA detectability half-life of pest prey in spider gut contents

BACKGROUND

Key natural enemy-pest interactions can be mapped in agricultural food webs by analysing predator gut content for the presence of a focal pest species. For this, PCR-based approaches are the most widely used methods providing the incidence of consumption of a focal pest in field sampled predators. To interpret such data the rate of prey DNA decay in the predators' gut, described by DNA detectability half-life (t 1/2), is needed. DNA decay may depend on the presence of alternative prey in the gut of generalist predators, but this effect has not been investigated in one of the major predatory arthropod groups, spiders.

METHODS

In a laboratory feeding experiment, we determined t 1/2 of the key cereal pest virus vector leafhopper Psammotettix alienus in the digestive tracts of its natural enemy, the spider Tibellus oblongus. We followed the fate of prey DNA in spiders which received only the focal prey as food, or as an alternative prey treatment they also received a meal of fruit flies after leafhopper consumption. After these feeding treatments, spiders were starved for variable time intervals prior to testing for leafhopper DNA in order to establish t 1/2.

RESULTS

We created a PCR protocol that detects P. alienus DNA in its spider predator. The protocol was further calibrated to the digestion speed of the spider by establishing DNA decay rate. Detectability limit was reached at 14 days, where c. 10% of the animals tested positive. The calculated t 1/2 = 5 days value of P. alienus DNA did not differ statistically between the treatment groups which received only the leafhopper prey or which also received fruit fly. The PCR protocol was validated in a field with known P. alienus infestation. In this applicability trial, we showed that 12.5% of field collected spiders were positive for the leafhopper DNA. We conclude that in our model system the presence of alternative prey did not influence the t 1/2 estimate of a pest species, which makes laboratory protocols more straightforward for the calibration of future field data.

The marine fish food web is globally connected

The productivity of marine ecosystems and the services they provide to humans are largely dependent on complex interactions between prey and predators. These are embedded in a diverse network of trophic interactions, resulting in a cascade of events following perturbations such as species extinction. The sheer scale of oceans, however, precludes the characterization of marine feeding networks through de novo sampling. This effort ought instead to rely on a combination of extensive data and inference. Here we investigate how the distribution of trophic interactions at the global scale shapes the marine fish food web structure. We hypothesize that the heterogeneous distribution of species ranges in biogeographic regions should concentrate interactions in the warmest areas and within species groups. We find that the inferred global metaweb of marine fish-that is, all possible potential feeding links between co-occurring species-is highly connected geographically with a low degree of spatial modularity. Metrics of network structure correlate with sea surface temperature and tend to peak towards the tropics. In contrast to open-water communities, coastal food webs have greater interaction redundancy, which may confer robustness to species extinction. Our results suggest that marine ecosystems are connected yet display some resistance to perturbations because of high robustness at most locations.

Review of Recent DNA-Based Methods for Main Food-Authentication Topics

Adulteration and mislabeling of food products and the commercial fraud derived, either intentionally or not, is a global source of economic fraud to consumers but also to all stakeholders involved in food production and distribution. Legislation has been enforced all over the world aimed at guaranteeing the authenticity of the food products all along the distribution chain, thereby avoiding food fraud and adulteration. Accordingly, there is a growing need for new analytical methods able to verify that all the ingredients included in a foodstuff match the qualities claimed by the manufacturer or distributor. In this sense, the improved performance of most recent DNA-based tools in term of sensitivity, multiplexing ability, high-throughput, and relatively low-cost give them a game-changing role in food-authenticity-related topics. Here, we provide a thorough and updated vision on the recently reported approaches that are applying these DNA-based tools to assess the authenticity of food components and products.

Intricate trophic links between threatened vertebrates confined to a small island in the Atlantic Ocean

Trophic networks in small isolated islands are in a fragile balance, and their disturbance can easily contribute toward the extinction vortex of species. Here, we show, in a small Atlantic island (Raso) in the Cabo Verde Archipelago, using DNA metabarcoding, the extent of trophic dependence of the Endangered giant wall gecko Tarentola gigas on endemic populations of vertebrates, including one of the rarest bird species of the world, the Critically Endangered Raso lark Alauda razae. We found that the Raso lark (27%), Iago sparrow Passer iagoensis (12%), Bulwer's petrel Bulweria bulwerii (15%), and the Cabo Verde shearwater Calonectris edwardsii (10%) are the most frequent vertebrate signatures found in the feces of the giant wall gecko. This work provides the first integrative assessment of their trophic links, an important issue to be considered for the long-term conservation of these small and isolated island ecosystems.

One-locus-several-primers: A strategy to improve the taxonomic and haplotypic coverage in diet metabarcoding studies

In diet metabarcoding analyses, insufficient taxonomic coverage of PCR primer sets generates false negatives that may dramatically distort biodiversity estimates. In this paper, we investigated the taxonomic coverage and complementarity of three cytochrome c oxidase subunit I gene (COI) primer sets based on in silico analyses and we conducted an in vivo evaluation using fecal and spider web samples from different invertivores, environments, and geographic locations. Our results underline the lack of predictability of both the coverage and complementarity of individual primer sets: (a) sharp discrepancies exist observed between in silico and in vivo analyses (to the detriment of in silico analyses); (b) both coverage and complementarity depend greatly on the predator and on the taxonomic level at which preys are considered; (c) primer sets' complementarity is the greatest at fine taxonomic levels (molecular operational taxonomic units [MOTUs] and variants). We then formalized the "one-locus-several-primer-sets" (OLSP) strategy, that is, the use of several primer sets that target the same locus (here the first part of the COI gene) and the same group of taxa (here invertebrates). The proximal aim of the OLSP strategy is to minimize false negatives by increasing total coverage through multiple primer sets. We illustrate that the OLSP strategy is especially relevant from this perspective since distinct variants within the same MOTUs were not equally detected across all primer sets. Furthermore, the OLSP strategy produces largely overlapping and comparable sequences, which cannot be achieved when targeting different loci. This facilitates the use of haplotypic diversity information contained within metabarcoding datasets, for example, for phylogeography and finer analyses of prey-predator interactions.

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.

Table for five, please: Dietary partitioning in boreal bats

Differences in diet can explain resource partitioning in apparently similar, sympatric species. Here, we analyzed 1,252 fecal droppings from five species (Eptesicus nilssonii, Myotis brandtii, M. daubentonii, M. mystacinus, and Plecotus auritus) to reveal their dietary niches using fecal DNA metabarcoding. We identified nearly 550 prey species in 13 arthropod orders. Two main orders (Diptera and Lepidoptera) formed the majority of the diet for all species, constituting roughly 80%-90% of the diet. All five species had different dietary assemblages. We also found significant differences in the size of prey species between the bat species. Our results on diet composition remain mostly unchanged when using either read counts as a proxy for quantitative diet or presence-absence data, indicating a strong biological pattern. We conclude that although bats share major components in their ecology (nocturnal life style, insectivory, and echolocation), species differ in feeding behavior, suggesting bats may have distinctive evolutionary strategies. Diet analysis helps illuminate life history traits of various species, adding to sparse ecological knowledge, which can be utilized in conservation planning.

Finding flies in the mushroom soup: Host specificity of fungus-associated communities revisited with a novel molecular method

Fruiting bodies of fungi constitute an important resource for thousands of other taxa. The structure of these diverse assemblages has traditionally been studied with labour-intensive methods involving cultivation and morphology-based species identification, to which molecular information might offer convenient complements. To overcome challenges in DNA extraction and PCR associated with the complex chemical properties of fruiting bodies, we developed a pipeline applicable for extracting amplifiable total DNA from soft fungal samples of any size. Our protocol purifies DNA in two sequential steps: (a) initial salt-isopropanol extraction of all nucleic acids in the sample is followed by (b) an extra clean-up step using solid-phase reversible immobilization (SPRI) magnetic beads. The protocol proved highly efficient, with practically all of our samples-regardless of biomass or other properties-being successfully PCR-amplified using metabarcoding primers and subsequently sequenced. As a proof of concept, we apply our methods to address a topical ecological question: is host specificity a major characteristic of fungus-associated communities, that is, do different fungus species harbour different communities of associated organisms? Based on an analysis of 312 fungal fruiting bodies representing 10 species in five genera from three orders, we show that molecular methods are suitable for studying this rich natural microcosm. Comparing to previous knowledge based on rearing and morphology-based identifications, we find a species-rich assemblage characterized by a low degree of host specialization. Our method opens up new horizons for molecular analyses of fungus-associated interaction webs and communities. Fruiting bodies of fungi constitute an important resource for thousands of other taxa. The structure of these diverse assemblages has traditionally been studied with labour-intensive methods involving cultivation and morphology-based species identification, to which molecular information might offer convenient complements. To overcome challenges in DNA extraction and PCR associated with the complex chemical properties of fruiting bodies, we developed a pipeline applicable for extracting amplifiable total DNA from soft fungal samples of any size. Our protocol purifies DNA in two sequential steps: (a) initial salt-isopropanol extraction of all nucleic acids in the sample is followed by (b) an extra clean-up step using solid-phase reversible immobilization (SPRI) magnetic beads. The protocol proved highly efficient, with practically all of our samples-regardless of biomass or other properties-being successfully PCR-amplified using metabarcoding primers and subsequently sequenced. As a proof of concept, we apply our methods to address a topical ecological question: is host specificity a major characteristic of fungus-associated communities, that is, do different fungus species harbour different communities of associated organisms? Based on an analysis of 312 fungal fruiting bodies representing 10 species in five genera from three orders, we show that molecular methods are suitable for studying this rich natural microcosm. Comparing to previous knowledge based on rearing and morphology-based identifications, we find a species-rich assemblage characterized by a low degree of host specialization. Our method opens up new horizons for molecular analyses of fungus-associated interaction webs and communities.

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.

The effects of pastoral intensification on the feeding interactions of generalist predators in streams

Land‐use change can alter trophic interactions with wide‐ranging functional consequences, yet the consequences for aquatic food webs have been little studied. In part, this may reflect the challenges of resolving the diets of aquatic organisms using classical gut contents analysis, especially for soft‐bodied prey. We used next‐generation sequencing to resolve prey use in nearly 400 individuals of two predatory invertebrates (the Caddisfly, Rhyacophila dorsalis, and the Stonefly Dinocras cephalotes) in streams draining land with increasingly intensive livestock farming. Rhyacophila dorsalis occurred in all streams, whereas D. cephalotes was restricted to low intensities, allowing us to test whether: (i) apparent sensitivity to agriculture in the latter species reflects a more specialized diet and (ii) diet in R. dorsalis varied between sites with and without D. cephalotes. DNA was extracted from dissected gut contents, amplified without blocking probes and sequenced using Ion Torrent technology. Both predators were generalists, consuming 30 prey taxa with a preference for taxa that were abundant in all streams or that increased with intensification. Where both predators were present, their diets were nearly identical, and R. dorsalis's diet was virtually unchanged in the absence of D. cephalotes. The loss of D. cephalotes from more intensive sites was probably due to physicochemical stressors, such as sedimentation, rather than to dietary specialization, although wider biotic factors (e.g., competition with other predatory taxa) could not be excluded. This study provides a uniquely detailed description of predator diets along a land‐use intensity gradient, offering new insights into how anthropogenic stressors affect stream communities.

Performance of DNA metabarcoding, standard barcoding, and morphological approach in the identification of host-parasitoid interactions

Understanding interactions between herbivores and parasitoids is essential for successful biodiversity protection and monitoring and for biological pest control. Morphological identifications employ insect rearing and are complicated by insects’ high diversity and crypsis. DNA barcoding has been successfully used in studies of host–parasitoid interactions as it can substantially increase the recovered real host–parasitoid diversity distorted by overlooked species complexes, or by species with slight morphological differences. However, this approach does not allow the simultaneous detection and identification of host(s) and parasitoid(s). Recently, high-throughput sequencing has shown high potential for surveying ecological communities and trophic interactions. Using mock samples comprising insect larvae and their parasitoids, we tested the potential of DNA metabarcoding for identifying individuals involved in host–parasitoid interactions to different taxonomic levels, and compared it to standard DNA barcoding and morphological approaches. For DNA metabarcoding, we targeted the standard barcoding marker cytochrome oxidase subunit I using highly degenerate primers, 2*300 bp sequencing on a MiSeq platform, and RTAX classification using paired-end reads. Additionally, using a large host–parasitoid dataset from a Central European floodplain forest, we assess the completeness and usability of a local reference library by confronting the number of Barcoding Index Numbers obtained by standard barcoding with the number of morphotypes. Overall, metabarcoding recovery was high, identifying 92.8% of the taxa present in mock samples, and identification success within individual taxonomic levels did not significantly differ among metabarcoding, standard barcoding, and morphology. Based on the current local reference library, 39.4% parasitoid and 90.7% host taxa were identified to the species level. DNA barcoding estimated higher parasitoid diversity than morphotyping, especially in groups with high level of crypsis. This study suggests the potential of metabarcoding for effectively recovering host–parasitoid diversity, together with more accurate identifications obtained from building reliable and comprehensive reference libraries, especially for parasitoids.

A framework for evaluating food-web responses to hydrological manipulations in riverine systems

Environmental flows are used to restore elements of the hydrological regime altered by human use of water. One of the primary justifications and purposes for environmental flows is the maintenance of target species populations but, paradoxically, there has been little emphasis on incorporating the food-web and trophic dynamics that determine population-level responses into the monitoring and evaluation of environmental flow programs. We develop a generic framework for incorporating trophic dynamics into monitoring programs to identify the food-web linkages between hydrological regimes and population-level objectives of environmental flows. These linkages form the basis for objective setting, ecological targets and indicator selection that are necessary for planning monitoring programs with a rigorous scientific basis. Because there are multiple facets of trophic dynamics that influence energy production and transfer through food webs, the specific objectives of environmental flows need to be defined during the development of monitoring programs. A multitude of analytical methods exist that each quantify distinct aspects of food webs (e.g. energy production, prey selection, energy assimilation), but no single method can provide a basis for holistic understanding of food webs. Our paper critiques a range of analytical methods for quantifying attributes of food webs to inform the setting, monitoring and evaluation of trophic outcomes of environmental flows and advance the conceptual understanding of trophic dynamics in river-floodplain systems.