The precautionary principle and dietary DNA metabarcoding: Commonly used abundance thresholds change ecological interpretation

Linking diet switching to reproductive performance across populations of two critically endangered mammalian herbivores

Optimal foraging theory predicts that animals maximise energy intake by consuming the most valuable foods available. When resources are limited, they may include lower-quality fallback foods in their diets. As seasonal herbivore diet switching is understudied, we evaluate its extent and effects across three Kenyan reserves each for Critically Endangered eastern black rhino (Diceros bicornis michaeli) and Grevy's zebra (Equus grevyi), and its associations with habitat quality, microbiome variation, and reproductive performance. Black rhino diet breadth increases with vegetation productivity (NDVI), whereas zebra diet breadth peaks at intermediate NDVI. Black rhino diets associated with higher vegetation productivity have less acacia (Fabaceae: Vachellia and Senegalia spp.) and more grass suggesting that acacia are fallback foods, upending conventional assumptions. Larger dietary shifts are associated with longer calving intervals. Grevy's zebra diets in high rainfall areas are consistently grass-dominated, whereas in arid areas they primarily consume legumes during low vegetation productivity periods. Whilst microbiome composition between individuals is affected by the environment, and diet composition in black rhino, seasonal dietary shifts do not drive commensurate microbiome shifts. Documenting diet shifts across ecological gradients can increase the effectiveness of conservation by informing habitat suitability models and improving understanding of responses to resource limitation.

Simultaneous species detection and discovery with environmental DNA metabarcoding: A freshwater mollusk case study

Environmental DNA (eDNA) sampling is a powerful tool for rapidly characterizing biodiversity patterns for specious, cryptic taxa with incomplete taxonomies. One such group that are also of high conservation concern are North American freshwater gastropods. In particular, springsnails of the genus Pyrgulopsis (Family: Hydrobiidae) are prevalent throughout the western United States where >140 species have been described. Many of the described species are narrow endemics known from a single spring or locality, and it is believed that there are likely many additional species which have yet to be described. The distribution of these species across the landscape is of interest because habitat loss and degradation, climate change, groundwater mining, and pollution have resulted in springsnail imperilment rates as high as 92%. Determining distributions with conventional sampling methods is limited by the fact that these snails are often <5 mm in length with few distinguishing morphological characters, making them both difficult to detect and to identify. We developed an eDNA metabarcoding protocol that is both inexpensive and capable of rapid, accurate detection of all known Pyrgulopsis species. When compared with conventional collection techniques, our pipeline consistently resulted in detection at sites previously known to contain Pyrgulopsis springsnails and at a cost per site that is likely to be substantially less than the conventional sampling and individual barcoding that has been done historically. Additionally, because our method uses eDNA extracted from filtered water, it is non-destructive and suitable for the detection of endangered species where "no take" restrictions may be in effect. This effort represents both a tool which is immediately applicable to taxa of high conservation concern across western North America and a case study in the broader application of eDNA sampling for landscape assessments of cryptic taxa of conservation concern.

Foraging history of individual elephants using DNA metabarcoding

Individual animals should adjust diets according to food availability. We used DNA metabarcoding to construct individual-level dietary timeseries for elephants from two family groups in Kenya varying in habitat use, social position and reproductive status. We detected at least 367 dietary plant taxa, with up to 137 unique plant sequences in one fecal sample. Results matched well-established trends: elephants tended to eat more grass when it rained and other plants when dry. Nested within these switches from 'grazing' to 'browsing' strategies, dietary DNA revealed seasonal shifts in food richness, composition and overlap between individuals. Elephants of both families converged on relatively cohesive diets in dry seasons but varied in their maintenance of cohesion during wet seasons. Dietary cohesion throughout the timeseries of the subdominant 'Artists' family was stronger and more consistently positive compared to the dominant 'Royals' family. The greater degree of individuality within the dominant family's timeseries could reflect more divergent nutritional requirements associated with calf dependency and/or priority access to preferred habitats. Whereas theory predicts that individuals should specialize on different foods under resource scarcity, our data suggest family bonds may promote cohesion and foster the emergence of diverse feeding cultures reflecting links between social behaviour and nutrition.

Multi-factorial examination of amplicon sequencing workflows from sample preparation to bioinformatic analysis

BACKGROUND

The development of sequencing technologies to evaluate bacterial microbiota composition has allowed new insights into the importance of microbial ecology. However, the variety of methodologies used among amplicon sequencing workflows leads to uncertainty about best practices as well as reproducibility and replicability among microbiome studies. Using a bacterial mock community composed of 37 soil isolates, we performed a comprehensive methodological evaluation of workflows, each with a different combination of methodological factors spanning sample preparation to bioinformatic analysis to define sources of artifacts that affect coverage, accuracy, and biases in the resulting compositional profiles.

RESULTS

Of the workflows examined, those using the V4-V4 primer set enabled the highest level of concordance between the original mock community and resulting microbiome sequence composition. Use of a high-fidelity polymerase, or a lower-fidelity polymerase with an increased PCR elongation time, limited chimera formation. Bioinformatic pipelines presented a trade-off between the fraction of distinct community members identified (coverage) and fraction of correct sequences (accuracy). DADA2 and QIIME2 assembled V4-V4 reads amplified by Taq polymerase resulted in the highest accuracy (100%) but had a coverage of only 52%. Using mothur to assemble and denoise V4-V4 reads resulted in a coverage of 75%, albeit with marginally lower accuracy (99.5%).

CONCLUSIONS

Optimization of microbiome workflows is critical for accuracy and to support reproducibility and replicability among microbiome studies. These considerations will help reveal the guiding principles of microbial ecology and impact the translation of microbiome research to human and environmental health.

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.

VTAM: A robust pipeline for validating metabarcoding data using controls

To obtain accurate estimates for biodiversity and ecological studies, metabarcoding studies should be carefully designed to minimize both false positive (FP) and false negative (FN) occurrences. Internal controls (mock samples and negative controls), replicates, and overlapping markers allow controlling metabarcoding errors but current metabarcoding software packages do not explicitly integrate these additional experimental data to optimize filtering. We have developed the metabarcoding analysis software VTAM, which uses explicitly these elements of the experimental design to find optimal parameter settings that minimize FP and FN occurrences. VTAM showed similar sensitivity, but a higher precision compared to two other pipelines using three datasets and two different markers (COI, 16S). The stringent filtering procedure implemented in VTAM aims to produce robust metabarcoding data to obtain accurate ecological estimates and represents an important step towards a non-arbitrary and standardized validation of metabarcoding data for conducting ecological studies. VTAM is implemented in Python and available from: https://github.com/aitgon/vtam. The VTAM benchmark code is available from: https://github.com/aitgon/vtam_benchmark.

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 suggests adaptive seasonal variation of individual trophic traits in a critically endangered fish

Dietary studies are critical for understanding foraging strategies and have important applications in conservation and habitat management. We applied a robust metabarcoding protocol to characterize the diet of the critically endangered freshwater fish Zingel asper (the Rhone streber). We conducted modelling and simulation analyses to identify and characterize some of the drivers of individual trophic trait variation in this species. We found that population density and ontogeny had minor effects on the trophic niche of Z. asper. Instead, our results suggest that the majority of trophic niche variation was driven by seasonal variation in ecological opportunity. The total trophic niche width of Z. asper seasonally expanded to include a broader range of prey. Furthermore, null model simulations revealed that the increase of between-individual variation in autumn indicates that Z. asper become more opportunistic relative to summer and spring, rather than being associated with a seasonal specialization of individuals. Overall, our results suggest an adaptive variation of individual trophic traits in Z. asper: the species mainly consumes a few ephemeropteran taxa (Baetis fuscatus and Ecdyonurus) but seems to be capable of adapting its foraging strategy to maintain its body condition. This study illustrates how metabarcoding data obtained from faeces can be validated and combined with individual-based modelling and simulation approaches to explore inter- and intrapopulational individual trophic traits variation and to test hypotheses in the conventional analytic framework of trophic ecology.

Evidence-based strategies to navigate complexity in dietary DNA metabarcoding: A reply

It is clearly beneficial to eliminate low-abundance sequences that arise in error during dietary DNA metabarcoding studies, but to purge all low-abundance sequences is to risk eliminating real sequences and complicating ecological analyses. Our prior literature review noted that DNA sequence relative read abundance (RRA) thresholds can help ameliorate false-positive taxon occurrences, but that historical emphasis on this utility has fostered uncertainty about the associated risk of inflating the false-negative rate (Littleford-Colquhoun et al., 2022). To address this, we combined a simulation study and an empirical data set to both illustrate the issue and provide blueprints for simulation studies and sensitivity analyses that can help investigators avoid overcorrecting and thereby bolster confidence in ecological inferences. Awareness of both the costs and the benefits of abundance-filtering is needed because accurately characterizing dietary distributions can be critically important for understanding animal diets, nutrition and trophic networks. Highlighting the need to raise awareness, a critique of our paper emphasized the misleading notion that "false positive interactions between species can present fundamentally incorrect network structures in network ecology, whereas false negatives will provide a correct but incomplete version of the network" (Tercel & Cuff, 2022). Asserting that the reliability of results will be eroded by false positives but resilient to the omission of true positives is risky and runs counter to evidence. Unfortunately, abundance-filtering methods can introduce false negatives at higher rates than they eliminate false positives and thereby undermine the analysis of otherwise reliable sequencing data. Overcorrecting can qualitatively alter and ultimately undermine ecological interpretations.

The precautionary principle and dietary DNA metabarcoding: Commonly used abundance thresholds change ecological interpretation

Dietary DNA metabarcoding enables researchers to identify and characterize trophic interactions with a high degree of taxonomic precision. It is also sensitive to sources of bias and contamination in the field and laboratory. One of the earliest and most common strategies for dealing with such sensitivities has been to remove all low-abundance sequences and conduct ecological analyses based on the presence or absence of food taxa. Although this step is now often perceived to be necessary, evidence of its sufficiency is lacking and more attention to the risk of introducing other errors is needed. Using computer simulations, we demonstrate that common strategies to remove low-abundance sequences can erroneously eliminate true dietary sequences in ways that impact downstream inferences. Using real data from well-studied wildlife populations in Yellowstone National Park, we further show how these strategies can markedly alter the composition of dietary profiles in ways that scale-up to obscure ecological interpretations about dietary generalism, specialism, and composition. Although the practice of removing low-abundance sequences may continue to be a useful strategy to address research questions that focus on a subset of relatively abundant foods, its continued widespread use risks generating misleading perceptions about the structure of trophic networks. Researchers working with dietary DNA metabarcoding data-or similar data such as environmental DNA, microbiomes, or pathobiomes-should be aware of drawbacks and consider alternative bioinformatic, experimental, and statistical solutions.