The metabarcoding of Grubs: Traditional herbal medicine of Scarabaeidae larvae

Grubs, called Qicao in China, have a long tradition as herbal medicine in East Asia. These larvae belong to the diverse family Scarabaeidae and are typically harvested from the wild during their immature stage based on morphological characteristics. However, rapid and accurate identification becomes challenging when relying solely on external morphological features, as the lack of clarity on biological sources raises safety concerns for clinical applications. The application of DNA metabarcoding provides a solution by enabling the determination of the biological source of a large sample. In the current study, we collected 19 batches of Grubs, consisting of 11,539 individuals, from the market and analyzed their biological composition through metabarcoding. We identified 49 Amplicon Sequence Variants (ASVs), 21 of which were Grubs. The 21 ASVs were classified into seven Molecular Operational Taxonomic Units (MOTUs) through species delimitation, which revealed that commercially available Grubs are predominantly sourced from Protaetia brevitarsis seulensis, while species of Rutelinae, Anomala, and Holotrichia were also abundant in some commercial batches. Among the identified ASVs, 28 belonged to non-Grub species and indicated adulteration from different animal families; high abundances of these ASVs were detected for Bombycidae, Tabanidae, and Viviparidae. Our findings underscore the complexity of Grubs' species composition and advocate for a deeper understanding of the wildlife sources contributing to herbal products. This research contributes valuable insights into the molecular identification of Grubs, paving the way for enhanced quality assurance in traditional medicine applications to provide safe and effective medicines for humanity.

DNA metabarcoding of cloacal swabs provides insight into diets of highly migratory sharks in the Mid-Atlantic Bight

The abundances of migratory shark species observed throughout the Mid-Atlantic Bight (MAB) during productive summer months suggest that this region provides critical habitat and prey resources to these taxa. However, the principal prey assemblages sustaining migratory shark biomass in this region are poorly defined. We applied high-throughput DNA metabarcoding to shark feces derived from cloacal swabs across nine species of Carcharhinid and Lamnid sharks to (1) quantify the contribution of broad taxa (e.g., invertebrates, fishes) supporting shark biomass during seasonal residency in the MAB and (2) determine whether the species displayed distinct dietary preference indicative of resource partitioning. DNA metabarcoding resulted in high taxonomic (species-level) resolution of shark diets with actinopterygian and elasmobranch fishes as the dominant prey categories across the species. DNA metabarcoding identified several key prey groups consistent across shark taxa that are likely integral for sustaining their biomass in this region, including Atlantic menhaden (Brevoortia tyrannus), Atlantic mackerel (Scomber scombrus), and benthic elasmobranchs, including skates. Our results are consistent with previously published stomach content data for the shark species of similar size range in the Northwest Atlantic Ocean, supporting the efficacy of cloacal swab DNA metabarcoding as a minimally invasive diet reconstruction technique. The high reliance of several shark species on Atlantic menhaden could imply wasp-waist food-web conditions during the summer months, whereby high abundances of forage fishes sustain a diverse suite of migratory sharks within a complex, seasonal food web.

DNA High-Throughput Sequencing for Arthropod Gut Content Analysis to Evaluate Effectiveness and Safety of Biological Control Agents

The search for effective biological control agents without harmful non-target effects has been constrained by the use of impractical (field direct observation) or imprecise (cage experiments) methods. While advances in the DNA sequencing methods, more specifically the development of high-throughput sequencing (HTS), have been quickly incorporated in biodiversity surveys, they have been slow to be adopted to determine arthropod prey range, predation rate and food web structure, and critical information to evaluate the effectiveness and safety of a biological control agent candidate. The lack of knowledge on how HTS methods could be applied by ecological entomologists constitutes part of the problem, although the lack of expertise and the high cost of the analysis also are important limiting factors. In this review, we describe how the latest HTS methods of metabarcoding and Lazaro, a method to identify prey by mapping unassembled shotgun reads, can serve biological control research, showing both their power and limitations. We explain how they work to determine prey range and also how their data can be used to estimate predation rates and subsequently be translated into food webs of natural enemy and prey populations helping to elucidate their role in the community. We present a brief history of prey detection through molecular gut content analysis and also the attempts to develop a more precise formula to estimate predation rates, a problem that still remains. We focused on arthropods in agricultural ecosystems, but most of what is covered here can be applied to natural systems and non-arthropod biological control candidates as well.

Species diversity analysis of commercial Mantidis Ootheca samples contaminated by store pests based on DNA metabarcoding

Mantidis Ootheca (Sangpiaoxiao, mantis egg case) is a typical multi-origin Chinese medicinal material. The Chinese Pharmacopoeia stipulates that the Mantidis Ootheca originates from three species of Mantis: Tenodera sinensis, Statilia maculate, and Hierodula patellifera. However, Mantidis Ootheca mainly relies on field collection, which leads to confusion of its actual origin in the market. As the clinical use of Mantidis Ootheca with unknown original mantis species will pose potential risks to drug safety, it is necessary to survey the commercially available Mantidis Ootheca origin species. However, as the egg case of Mantis, the morphological characters of Mantidis Ootheca are limited and usually cannot serve as accurate identification tool. DNA barcoding, which is widely used in taxonomic studies of animals, is severely affected by the impact of storage pests and DNA degradation. Thus, this study collected a total of 4580 Mantidis Ootheca and pooled separately Mantidis Ootheca samples according to 18 different sources as DNA samples to analyze the origin diversity of Mantidis Ootheca individuals contaminated by common store pests collected in in the market using DNA metabarcoding, and to provide a basis for quality control of Mantidis Ootheca. 37 Mantis ASVs and 9 Mantis MOTUs were identified through species delimitation, and the high-level intraspecific diversity was depicted as haplotype network plot. Besides Tenodera sinensis and Hierodula patellifera as genuine original mantis species defined in the Chinese Pharmacopoeia, Tenodera angustipennis was also the origin species of these Mantidis Ootheca samples.

A novel biological sources consistency evaluation method reveals high level of biodiversity within wild natural medicine: A case study of Amynthas earthworms as “Guang Dilong”

For wild natural medicine, unanticipated biodiversity as species or varieties with similar morphological characteristics and sympatric distribution may co-exist in a single batch of medical materials, which affects the efficacy and safety of clinical medication. DNA barcoding as an effective species identification tool is limited by its low sample throughput nature. In this study, combining DNA mini-barcode, DNA metabarcoding and species delimitation method, a novel biological sources consistency evaluation strategy was proposed, and high level of interspecific and intraspecific variations were observed and validated among 5376 Amynthas samples from 19 sampling points regarded as "Guang Dilong" and 25 batches of proprietary Chinese medicines. Besides Amynthas aspergillum as the authentic source, 8 other Molecular Operational Taxonomic Units (MOTUs) were elucidated. Significantly, even the subgroups within A. aspergillum revealed here differ significantly on chemical compositions and biological activity. Fortunately, this biodiversity could be controlled when the collection was limited to designated areas, as proved by 2796 "decoction pieces" samples. This batch biological identification method should be introduced as a novel concept regarding natural medicine quality control, and to offer guidelines for in-situ conservation and breeding bases construction of wild natural medicine.

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

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

Biological composition analysis of a natural medicine, Faeces Vespertilionis, with complex sources using DNA metabarcoding

Faeces Vespertilionis is a commonly used fecal traditional Chinese medicine. Traditionally, it is identified relying only on morphological characters. This poses a serious challenge to the composition analysis accuracy of this complex biological mixture. Thus, for quality control purposes, an accurate and effective method should be provided for taxonomic identification of Faeces Vespertilionis. In this study, 26 samples of Faeces Vespertilionis from ten provinces in China were tested using DNA metabarcoding. Seven operational taxonomic units (OTUs) were detected as belonging to bats. Among them, Hipposideros armiger (Hodgson, 1835) and Rhinolophus ferrumequinum (Schober and Grimmberger, 1997) were the main host sources of Faeces Vespertilionis samples, with average relative abundances of 59.3% and 24.1%, respectively. Biodiversity analysis showed that Diptera and Lepidoptera were the most frequently consumed insects. At the species level, 19 taxa were clearly identified. Overall, our study used DNA metabarcoding to analyze the biological composition of Faeces Vespertilionis, which provides a new idea for the quality control of this special traditional Chinese medicine.

Diet diversity and environment determine the intestinal microbiome and bacterial pathogen load of fire salamanders

Diverse communities of symbiotic microbes inhabit the digestive systems of vertebrates and play a crucial role in animal health, and host diet plays a major role in shaping the composition and diversity of these communities. Here, we characterized diet and gut microbiome of fire salamander populations from three Belgian forests. We carried out DNA metabarcoding on fecal samples, targeting eukaryotic 18S rRNA of potential dietary prey items, and bacterial 16S rRNA of the concomitant gut microbiome. Our results demonstrated an abundance of soft-bodied prey in the diet of fire salamanders, and a significant difference in the diet composition between males and females. This sex-dependent effect on diet was also reflected in the gut microbiome diversity, which is higher in males than female animals. Proximity to human activities was associated with increased intestinal pathogen loads. Collectively, the data supports a relationship between diet, environment and intestinal microbiome in fire salamanders, with potential health implications.

High-throughput sequencing reveals dietary segregation in Malaysian babblers

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

Molecular gut content analysis indicates the inter- and intra-guild predation patterns of spiders in conventionally managed vegetable fields

Inter- and intra-guild interactions are important in the coexistence of predators and their prey, especially in highly disturbed vegetable cropping systems with sporadic food resources. Assessing the dietary range of a predator taxon characterized by diverse foraging behavior using conventional approaches, such as visual observation and conventional molecular approaches for prey detection, has serious logistical problems. In this study, we assessed the prey compositions and compare the dietary spectrum of a functionally diverge group of predators-spiders-to characterize their trophic interactions and assess biological control potential in Brassica vegetable fields. We used high-throughput sequencing (HTS) and biotic interaction networks to precisely annotate the predation spectrum and highlight the predator-predator and predator-prey interactions. The prey taxa in the gut of all spider families were mainly enriched with insects (including dipterans, coleopterans, orthopterans, hemipterans, and lepidopterans) with lower proportions of arachnids (such as Araneae) along with a wide range of other prey factions. Despite the generalist foraging behavior of spiders, the community structure analysis and interaction networks highlighted the overrepresentation of particular prey taxa in the gut of each spider family, as well as showing the extent of interfamily predation by spiders. Identifying the diverse trophic niche proportions underpins the importance of spiders as predators of pests in highly disturbed agroecosystems. More specifically, combining HTS with advanced ecological community analysis reveals the preferences and biological control potential of particular spider taxa (such as Salticidae against lepidopterans and Pisauridae against dipterans), and so provides a valuable evidence base for targeted conservation biological control efforts in complex trophic networks.

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

Dietary analysis using DNA metabarcoding is a powerful tool that is increasingly being used to further our knowledge of trophic interactions in highly complex food webs but is not without limitations. Omnivores, the most generalist of consumers, pose unique challenges when using such methods. Here, we provide the rationale to understand the problems associated with analysing the complex diets of omnivores. By reviewing existing metabarcoding studies of omnivorous diet, and constructing hypothetical scenarios arising from each, we outline that great caution is required when interpreting sequencing data in such cases. In essence, the problems of accidental consumption and secondary ingestion are significant sources of error when investigating omnivorous diets. The integration of multiple high throughput sequencing markers increases the taxonomic breadth of taxa detected but we reveal how some detections may be misleading. Disentangling which taxa have been deliberately or accidentally consumed by the focal omnivore is challenging and can falsely emphasise those that were not intentionally consumed, obscuring biologically meaningful interactions. Although we suggest ways to disentangle these issues, we urge that the results of such analyses should be interpreted with caution and all possible scenarios for the presence of biota within omnivores given due consideration.

A powerful long metabarcoding method for the determination of complex diets from faecal analysis of the European pond turtle (Emys orbicularis, L. 1758)

High‐throughput sequencing has become an accurate method for the identification of species present in soil, water, faeces, gut or stomach contents. However, information at the species level is limited due to the choice of short barcodes and based on the idea that DNA is too degraded to allow longer sequences to be amplified. We have therefore developed a long DNA metabarcoding method based on the sequencing of short reads followed by de novo assembly, which can precisely identify the taxonomic groups of organisms associated with complex diets, such as omnivorous individuals. The procedure includes 11 different primer pairs targeting the COI gene, the large subunit of the ribulose‐1,5‐bisphosphate carboxylase gene, the maturase K gene, the 28S rRNA and the trnL‐trnF chloroplastic region. We validated this approach using 32 faeces samples from an omnivorous reptile, the European pond turtle (Emys orbicularis, L. 1758). This metabarcoding approach was assessed using controlled experiments including mock communities and faecal samples from captive feeding trials. The method allowed us to accurately identify prey DNA present in the diet of the European pond turtles to the species level in most of the cases (82.4%), based on the amplicon lengths of multiple markers (168–1,379 bp, average 546 bp), and produced by de novo assembly. The proposed approach can be adapted to analyse various diets, in numerous conservation and ecological applications. It is consequently appropriate for detecting fine dietary variations among individuals, populations and species as well as for the identification of rare food items.

High-resolution multi-marker DNA metabarcoding reveals sexual dietary differentiation in a bird with minor dimorphism

Although sexual dietary differentiation is well known in birds, it is usually linked with significant morphological dimorphism between males and females, with lower differentiation reported in sexually monomorphic or only slightly dimorphic species. However, this may be an artifact of poor taxonomic resolution achieved in most conventional dietary studies, which may be unable to detect subtle intraspecific differentiation in prey consumption. Here, we show the power of multi‐marker metabarcoding to address these issues, focusing on a slightly dimorphic generalist passerine, the black wheatear Oenanthe leucura. Using markers from four genomic regions (18S, 16S, COI, and trnL), we analyzed fecal droppings collected from 93 adult black wheatears during the breeding season. We found that sexes were rather similar in bill and body features, though males had a slightly thicker bill and longer wings and tail than females. Diet was dominated in both sexes by a very wide range of arthropod species and a few fleshy fruits, but the overall diet diversity was higher for males than females, and there was a much higher frequency of occurrence of ants in female (58%) than male (29%) diets. We hypothesize that the observed sexual differentiation was likely related to females foraging closer to their offspring on abundant prey, while males consumed a wider variety of prey while foraging more widely. Overall, our results suggest that dietary sexual differentiation in birds may be more widespread than recognized at present and that multi‐marker DNA metabarcoding is a particularly powerful tool to unveiling such differences.

Questioning the proverb 'more haste, less speed': classic versus metabarcoding approaches for the diet study of a remote island endemic gecko

Dietary studies can reveal valuable information on how species exploit their habitats and are of particular importance for insular endemics conservation as these species present higher risk of extinction. Reptiles are often neglected in island systems, principally the ones inhabiting remote areas, therefore little is known on their ecological networks. The Selvagens gecko Tarentola (boettgeri) bischoffi, endemic to the remote and integral reserve of Selvagens Archipelago, is classified as Vulnerable by the Portuguese Red Data Book. Little is known about this gecko’s ecology and dietary habits, but it is assumed to be exclusively insectivorous. The diet of the continental Tarentola species was already studied using classical methods. Only two studies have used next-generation sequencing (NGS) techniques for this genus thus far, and very few NGS studies have been employed for reptiles in general. Considering the lack of information on its diet and the conservation interest of the Selvagens gecko, we used morphological and DNA metabarcoding approaches to characterize its diet. The traditional method of morphological identification of prey remains in faecal pellets collected over a longer period was compared with metabarcoding of samples collected during rapid surveys. Molecular results revealed that this species is a generalist, feeding on invertebrate, plant and vertebrate items, whereas the morphological approaches were unable to detect the latter two. These results opened up new questions on the ecological role of the Selvagens gecko that deserves to be further explored, such as the possible predation on seabirds, plant services or trophic competition with the sympatric Madeira lizard Teira dugesii. Metabarcoding identified a greater diversity of dietary items at higher taxonomic resolution, but morphological identification enabled calculation of relative abundances and biomasses of ingested arthropods, and detected a dietary shift on invertebrate preys between seasons. Results of this study highlight the global applicability of rapid metabarcoding surveys for understudied taxa on remote islands that are difficult to access. We recommend using the metabarcoding approach, even if ‘speedy’ sampling only is possible, but we must highlight that disregarding long-term ecological data may lead to ‘hasty’ conclusion.

What Is the Giant Wall Gecko Having for Dinner? Conservation Genetics for Guiding Reserve Management in Cabo Verde

Knowledge on diet composition of a species is an important step to unveil its ecology and guide conservation actions. This is especially important for species that inhabit remote areas within biodiversity hotspots, with little information about their ecological roles. The emblematic giant wall gecko of Cabo Verde, Tarentola gigas, is restricted to the uninhabited Branco and Raso islets, and presents two subspecies. It is classified as Endangered, and locally Extinct on Santa Luzia Island; however, little information is known about its diet and behaviour. In this study, we identified the main plant, arthropods, and vertebrates consumed by both gecko subspecies using next generation sequencing (NGS) (metabarcoding of faecal pellets), and compared them with the species known to occur on Santa Luzia. Results showed that plants have a significant role as diet items and identified vertebrate and invertebrate taxa with higher taxonomic resolution than traditional methods. With this study, we now have data on the diet of both subspecies for evaluating the reintroduction of this threatened gecko on Santa Luzia as potentially successful, considering the generalist character of both populations. The information revealed by these ecological networks is important for the development of conservation plans by governmental authorities, and reinforces the essential and commonly neglected role of reptiles on island systems.

Evaluating metabarcoding to analyse diet composition of species foraging in anthropogenic landscapes using Ion Torrent and Illumina sequencing

DNA metabarcoding of faecal samples is being successfully used to study the foraging niche of species. We assessed the ability of two benchtop high-throughput sequencing (HTS) platforms, to identify a large taxonomic array of food items from domestic cats Felis silvestris catus, including prey and human-related food taxa (pet food and leftovers leaving undetectable solid remains in faeces). Scats from a captive feeding trial (n = 41) and from free-ranging individuals (n = 326) were collected and analysed using a cytb mini-barcode in independent PCR replicates on the Ion PGM and the MiSeq platforms. Outputs from MiSeq were more sensitive and reproducible than those from Ion PGM due to a higher sequencing depth and sequence quality on MiSeq. DNA from intact prey taxa was detected more often (82% of the expected occurrences) than DNA from pet food (54%) and raw fish and meat (31%). We assumed that this variability was linked to different degree of DNA degradation: The Ion PGM detected significantly less human-linked food, birds, field voles, murids and shrews in the field-collected samples than the MiSeq platform. Pooling the replicates from both platforms and filtering the data allowed identification of at least one food item in 87.4% of the field-collected samples. Our DNA metabarcoding approach identified 29 prey taxa, of which 25 to species level (90% of items) including 9 rodents, 3 insectivores, 12 birds and 1 reptile and 33 human-related food taxa of which 23 were identified to genus level (75% of items). Our results demonstrate that using HTS platforms such as MiSeq, which provide reads of sufficiently high quantity and quality, with sufficient numbers of technical replicates, is a robust and non-invasive approach for further dietary studies on animals foraging on a wide range of food items in anthropogenic landscapes.

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

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

Molecular gut content analysis of different spider body parts

Molecular gut-content analysis has revolutionized the study of food webs and feeding interactions, allowing the detection of prey DNA within the gut of many organisms. However, successful prey detection is a challenging procedure in which many factors affect every step, starting from the DNA extraction process. Spiders are liquid feeders with branched gut diverticula extending into their legs and throughout the prosoma, thus digestion takes places in different parts of the body and simple gut dissection is not possible. In this study, we investigated differences in prey detectability in DNA extracts from different parts of the spider´s body: legs, prosoma and opisthosoma, using prey-specific PCR and metabarcoding approaches. We performed feeding trials with the woodlouse hunter spider Dysdera verneaui Simon, 1883 (Dysderidae) to estimate the time at which prey DNA is detectable within the predator after feeding. Although we found that all parts of the spider body are suitable for gut-content analysis when using prey-specific PCR approach, results based on metabarcoding suggested the opisthosoma is optimal for detection of predation in spiders because it contained the highest concentration of prey DNA for longer post feeding periods. Other spiders may show different results compared to D. verneaui, but given similarities in the physiology and digestion in different families, it is reasonable to assume this to be common across species and this approach having broad utility across spiders.

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.

DNA metabarcoding for diet analysis and biodiversity: A case study using the endangered Australian sea lion (Neophoca cinerea)

The analysis of apex predator diet has the ability to deliver valuable insights into ecosystem health, and the potential impacts a predator might have on commercially relevant species. The Australian sea lion (Neophoca cinerea) is an endemic apex predator and one of the world's most endangered pinnipeds. Given that prey availability is vital to the survival of top predators, this study set out to understand what dietary information DNA metabarcoding could yield from 36 sea lion scats collected across 1,500 km of its distribution in southwest Western Australia. A combination of PCR assays were designed to target a variety of potential sea lion prey, including mammals, fish, crustaceans, cephalopods, and birds. Over 1.2 million metabarcodes identified six classes from three phyla, together representing over 80 taxa. The results confirm that the Australian sea lion is a wide-ranging opportunistic predator that consumes an array of mainly demersal fauna. Further, the important commercial species Sepioteuthis australis (southern calamari squid) and Panulirus cygnus (western rock lobster) were detected, but were present in <25% of samples. Some of the taxa identified, such as fish, sharks and rays, clarify previous knowledge of sea lion prey, and some, such as eel taxa and two gastropod species, represent new dietary insights. Even with modest sample sizes, a spatial analysis of taxa and operational taxonomic units found within the scat shows significant differences in diet between many of the sample locations and identifies the primary taxa that are driving this variance. This study provides new insights into the diet of this endangered predator and confirms the efficacy of DNA metabarcoding of scat as a noninvasive tool to more broadly define regional biodiversity.

Molecular assessment of predation by hoverflies (Diptera: Syrphidae) in Mediterranean lettuce crops

BACKGROUND

Hoverflies (Diptera: Syrphidae) are generalist predators of a great variety of pests. Nasonovia ribisnigri (Hemiptera: Aphididae) and Frankliniella occidentalis (Thysanoptera: Thripidae) are two common pests in Mediterranean lettuce crops, where they occur alongside alternative prey (e.g., Collembola). A semi-field experiment was conducted in an experimental lettuce plot where hoverfly predation on N. ribisnigri, F. occidentalis and Collembola was studied by conventional PCR and qPCR using specific primers, as well as by next-generation sequencing (NGS) in order to reveal other potential trophic interactions.

RESULTS

Trophic linkages between hoverflies and N. ribisnigri were the strongest both in spring and summer. F. occidentalis and Collembolans were also detected in both seasons, but with less frequency. qPCR detected a higher frequency of consumption than conventional PCR when both tests were run at optimal conditions. NGS analyses showed intraguild predation on other hoverfly species, as well as on anthocorids, spiders and even aphid parasitoids.

CONCLUSIONS

Conventional PCR and qPCR provided important insights into Mediterranean hoverfly species predation on target pest and non-pest prey. NGS gave a complementary approach revealing a broader diet of these predators within the studied ecosystem.

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

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

Sexual dimorphism of tropical green pit viper Trimeresurus (Cryptelytrops) macrops in Northeast Thailand

Total body size for arboreal vipers is normally biased toward small males and larger females. We evaluated sexual dimorphism in adult and sub-adultTrimeresurus macrops, a small, arboreal green pit viper commonly found in Southeast Asia, but severely understudied. We evaluated morphological characters of 139 adultT. macropsobtained by active and opportunistic searches at Sakaerat Biosphere Reserve in northeast Thailand, from May 2012 to October 2014. We compared 7 external characters including residual index (Ri) and scaled mass index (SMI) between male and femaleT. macrops. Body length measurements and mass differed between males and females. Females had greater snout-vent length, body mass and head length and width, while males had greater tail lengths. A postocular stripe was always present on males, but never on females. Male head size was negatively correlated with SMI, which may reveal intersexual competition inT. macrops. Sexually dimorphic characters inT. macropsmay have evolved through intraspecific resource partitioning.

Universal and blocking primer mismatches limit the use of high-throughput DNA sequencing for the quantitative metabarcoding of arthropods

The quantification of the biological diversity in environmental samples using high-throughput DNA sequencing is hindered by the PCR bias caused by variable primer-template mismatches of the individual species. In some dietary studies, there is the added problem that samples are enriched with predator DNA, so often a predator-specific blocking oligonucleotide is used to alleviate the problem. However, specific blocking oligonucleotides could coblock nontarget species to some degree. Here, we accurately estimate the extent of the PCR biases induced by universal and blocking primers on a mock community prepared with DNA of twelve species of terrestrial arthropods. We also compare universal and blocking primer biases with those induced by variable annealing temperature and number of PCR cycles. The results show that reads of all species were recovered after PCR enrichment at our control conditions (no blocking oligonucleotide, 45 °C annealing temperature and 40 cycles) and high-throughput sequencing. They also show that the four factors considered biased the final proportions of the species to some degree. Among these factors, the number of primer-template mismatches of each species had a disproportionate effect (up to five orders of magnitude) on the amplification efficiency. In particular, the number of primer-template mismatches explained most of the variation (~3/4) in the amplification efficiency of the species. The effect of blocking oligonucleotide concentration on nontarget species relative abundance was also significant, but less important (below one order of magnitude). Considering the results reported here, the quantitative potential of the technique is limited, and only qualitative results (the species list) are reliable, at least when targeting the barcoding COI region.

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

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

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

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

Pyrosequencing of prey DNA in faeces of carnivorous land snails to facilitate ecological restoration and relocation programmes

Identifying and understanding predator diets is of high importance in biological conservation. This is particularly true for the introduction, establishment and maintenance of predator populations in newly created or modified ecological communities, such as translocation sites or restored habitats. Conservation status of predators may not permit captive feeding trials or intrusive gut-content methods, so non-intrusive diet assessment is required, such as faecal analysis. However, prey such as earthworms leave no morphological clues suitable for accurately discriminating between species consumed through visual faecal analysis. This study uses non-intrusive molecular methods on earthworm DNA extracted from the faeces of the carnivorous land snail Powelliphanta patrickensis to identify its earthworm diet and any seasonal trends. Data from 454-pyrosequencing revealed earthworm DNA in all samples (n = 60). Sequences were compared to a DNA library created from published and unpublished studies of New Zealand's endemic earthworms and online databases. Unidentified earthworm sequences were clustered into molecular operational taxonomic units (MOTUs). Twenty-six MOTUs were identified, 17 of which matched the library, whereas nine did not. Similarity indices indicate that there were seasonal differences (P < 0.05) in the earthworm communities represented in the summer and the winter diets. This study highlights the importance of utilising the vast body of data generated by pyrosequencing to investigate potential temporal diet shifts in protected species. The method described here is widely applicable to a wide range of predatory species of conservation interest and can further inform habitat restoration and relocation programmes to optimize the long-term survival of the target species.

From puffins to plankton: a DNA-based analysis of a seabird food chain in the northern Gulf of Maine

The predator-prey interactions within food chains are used to both characterize and understand ecosystems. Conventional methods of constructing food chains from visual identification of prey in predator diet can suffer from poor taxonomic resolution, misidentification, and bias against small or completely digestible prey. Next-generation sequencing (NGS) technology has become a powerful tool for diet reconstruction through barcoding of DNA in stomach content or fecal samples. Here we use multi-locus (16S and CO1) next-generation sequencing of DNA barcodes on the feces of Atlantic puffin (Fratercula arctica) chicks (n=65) and adults (n=64) and the stomach contents of their main prey, Atlantic herring (Clupea harengus, n=44) to investigate a previously studied food chain. We compared conventional and molecular-derived chick diet, tested the similarity between the diets of puffin adults and chicks, and determined whether herring prey can be detected in puffin diet samples. There was high variability in the coverage of prey groups between 16S and CO1 markers. We identified more unique prey with our 16S compared to CO1 barcoding markers (51 and 39 taxa respectively) with only 12 taxa identified by both genes. We found no significant difference between the 16S-identified diets of puffin adults (n=17) and chicks (n=41). Our molecular method is more taxonomically resolved and detected chick prey at higher frequencies than conventional field observations. Many likely planktonic prey of herring were detected in feces from puffin adults and chicks, highlighting the impact secondary consumption may have on the interpretation of molecular dietary analysis. This study represents the first simultaneous molecular investigation into the diet of multiple components of a food chain and highlights the utility of a multi-locus approach to diet reconstruction that is broadly applicable to food web analysis.

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

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

Dietary competition between the alien Asian Musk Shrew (Suncus murinus) and a re-introduced population of Telfair's Skink (Leiolopisma telfairii)

Re-introduction of rare species to parts of their historical range is becoming increasingly important as a conservation strategy. Telfair's Skinks (Leiolopisma telfairii), once widespread on Mauritius, were until recently found only on Round Island. There it is vulnerable to stochastic events, including the introduction of alien predators that may either prey upon it or compete for food resources. Consequently, skinks have been introduced to Ile aux Aigrettes, another small Mauritian island that has been cleared of rats. However, the island has been invaded by Asian Musk Shrews (Suncus murinus), a commensal species spread by man well beyond its natural Asian range. Our aim was to use next-generation sequencing to analyse the diets of the shrews and skinks to look for niche competition. DNA was extracted from skink faeces and from the stomach contents of shrews. Application of shrew- and skink-specific primers revealed no mutual predation. The DNA was then amplified using general invertebrate primers with tags to identify individual predators, and then sequenced by 454 pyrosequencing. 119 prey MOTUs (molecular taxonomic units) were isolated, although none could be identified to species. Seeding of cladograms with known sequences allowed higher taxonomic assignments in some cases. Although most MOTUs were not shared by shrews and skinks, Pianka's niche overlap test showed significant prey overlap, suggesting potentially strong competition where food resources are limited. These results suggest that removal of the shrews from the island should remain a priority.

Storage and shipping of tissue samples for DNA analyses: A case study on earthworms

Nowadays, molecular analyses play an important role in studies of soil dwelling animals, for example in taxonomy, phylogeography or food web analyses. The quality of the DNA, used for later molecular analyses, is an important factor and depends on collection and preservation of samples prior to DNA extraction. Ideally, DNA samples are frozen immediately upon collection, but if samples are collected in the field, suitable preservation methods might be limited due to unavailability of resources or remote field sites. Moreover, shipping samples over long distances can cause loss of DNA quality e.g. by thawing or leaking of preservation liquid. In this study we use earthworms, a key organism in soil research, to compare three different DNA preservation methods - freezing at -20 °C, storing in 75% ethanol, and freeze drying. Samples were shipped from the United States of America to Austria. The DNA of the samples was extracted using two different extraction methods, peqGOLD™ and Chelex^® 100. The DNA amplification success was determined by amplifying four DNA fragments of different length. The PCR amplification success is significantly influenced by preservation method and extraction method and differed significantly depending on the length of the DNA fragment. Freeze drying samples was the best preservation method when samples were extracted using the silica based extraction method peqGOLD™. For samples that were extracted with Chelex^® 100, storage in ethanol was the best preservation method. However, the overall amplification success was significantly lower for the extraction procedure based on Chelex^® 100. The detection of the small DNA fragments was higher and independent from the extraction method, while the amplification success was significantly reduced for the longer DNA fragments. We recommend freeze drying of DNA samples, especially when they have to be shipped for longer distances. No special packaging or declaration is needed for freeze dried samples, and the risk of thawing is excluded. Storage of freeze dried samples also reduces costs because samples can be kept at room temperature in a desiccator. It should be noted, that the extraction methods showed significant differences in DNA amplification success. Thus, the extraction method should be taken into account when choosing the preservation method.

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

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

Ultra-deep sequencing enables high-fidelity recovery of biodiversity for bulk arthropod samples without PCR amplification

BACKGROUND

Next-generation-sequencing (NGS) technologies combined with a classic DNA barcoding approach have enabled fast and credible measurement for biodiversity of mixed environmental samples. However, the PCR amplification involved in nearly all existing NGS protocols inevitably introduces taxonomic biases. In the present study, we developed new Illumina pipelines without PCR amplifications to analyze terrestrial arthropod communities.

RESULTS

Mitochondrial enrichment directly followed by Illumina shotgun sequencing, at an ultra-high sequence volume, enabled the recovery of Cytochrome c Oxidase subunit 1 (COI) barcode sequences, which allowed for the estimation of species composition at high fidelity for a terrestrial insect community. With 15.5 Gbp Illumina data, approximately 97% and 92% were detected out of the 37 input Operational Taxonomic Units (OTUs), whether the reference barcode library was used or not, respectively, while only 1 novel OTU was found for the latter. Additionally, relatively strong correlation between the sequencing volume and the total biomass was observed for species from the bulk sample, suggesting a potential solution to reveal relative abundance.

CONCLUSIONS

The ability of the new Illumina PCR-free pipeline for DNA metabarcoding to detect small arthropod specimens and its tendency to avoid most, if not all, false positives suggests its great potential in biodiversity-related surveillance, such as in biomonitoring programs. However, further improvement for mitochondrial enrichment is likely needed for the application of the new pipeline in analyzing arthropod communities at higher diversity.

|SE|S|AM|E| Barcode: NGS-oriented software for amplicon characterization--application to species and environmental barcoding

Progress in NGS technologies has opened up new opportunities for characterizing biodiversity, both for individual specimen identification and for environmental barcoding. Although the amount of data available to biologist is increasing, user-friendly tools to facilitate data analysis have yet to be developed. Our aim, with |SE|S|AM|E| Barcode, is to provide such support through a unified platform. The sequences are analysed through a pipeline that (i) processes NGS amplicon runs, filtering markers and samples, (ii) builds reference libraries and finally (iii) identifies (barcodes) the sequences in each amplicon from the reference library. We use a simulated data set for specimen identification and a recently published data set for environmental barcoding to validate the method. The results obtained are consistent with the expected characterizations (in silico and previously published, respectively). |SE|S|AM|E| Barcode and its documentation are freely available under the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported Licence for Windows and Linux from http://www1.montpellier.inra.fr/CBGP/NGS/.

Who is eating what: diet assessment using next generation sequencing

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