Metabarcoding of Fecal Samples to Determine Herbivore Diets: A Case Study of the Endangered Pacific Pocket Mouse

DNA metabarcoding analysis of the bare-nosed wombat (Vombatus ursinus) diet

Wombats are native herbivorous grazers that have adapted to Australia's low-quality forage. Studies on diet selection by bare-nosed wombats (Vombatus ursinus) are limited and are either observational or based on microhistological studies. The current study determined the diet of wombats through DNA metabarcoding across five study sites in New South Wales over a one-year period. Metabarcoding was chosen as it is non-invasive, less time consuming and more specific than traditional techniques. The list of 209 plant species identified as eaten by wombats in this study is much higher than previously reported, with grasses being the most common plant group identified in all samples. Most dietary items identified were introduced plant species. Seasonal differences in plants eaten occurred at four of the five study sites and may reflect dietary abundance and floristic composition at different times of year. Further studies are required to determine if the dietary items differ markedly across the entire range of wombats, and if nutrition influences dietary preferences.

Deer dietary responses to wildfire: Optimal foraging, individual specialization, or opportunism?

Increasing impacts of wildfire on arid regions of the world fuelled by climate change highlight the need to better understand how natural communities respond to fire. We took advantage of a large (1660-km2 ) wildfire that erupted in northern California during an in-progress study of black-tailed deer (Odocoileus hemionus columbianus) to investigate deer use of and diets within burned and unburned habitats before and after the fire. We compared deer diet breadth to predictions of optimal foraging theory, the niche variation hypothesis, and opportunistic (i.e., generalist) foraging expectations under the assumption that overall availability and diversity of forage in burned areas declined immediately after the fire and increased as the plant community recovered in the next 3 years after the fire. We used faecal pellet counts to document space use and metabarcoding to study diet during pre-fire, post-fire, and recovery periods. Pellet counts supported predictions that deer increased use of unburned sites and reduced use of burn sites after the fire and began to return to burned sites in subsequent sampling years. Diet diversity did not differ significantly between control and burn sites before the fire, but was lower in burn than control sites post-fire (p < .001), when and where diet was dominated by oak (Quercus spp). In contrast, during subsequent years, diet diversity was higher (including more herbaceous plants) in burn than control sites (p < .05). In contrast to predictions of optimal foraging and niche variation hypotheses, individual deer foraged as generalists for which changes in dietary niche breadth paralleled fire-induced changes in diversity of the plant community.

Unraveling the potential of environmental DNA for deciphering recent advances in plant-animal interactions: a systematic review

MAIN CONCLUSION

Environmental DNA-based monitoring provides critical insights for enhancing our understanding of plant-animal interactions in the context of worldwide biodiversity decrease for developing a global framework for effective plant biodiversity conservation. To understand the ecology and evolutionary patterns of plant-animal interactions (PAI) and their pivotal roles in ecosystem functioning advances in molecular ecology tools such as Environmental DNA (eDNA) provide unprecedented research avenues. These methods being non-destructive in comparison to traditional biodiversity monitoring methods, enhance the discernment of ecosystem health, integrity, and complex interactions. This review intends to offer a systematic and critical appraisal of the prospective of eDNA for investigating PAI. The review thoroughly discusses and analyzes the recent reports (2015-2022) employing preferred reporting items for systematic reviews and meta-analyses (PRISMA) to outline the recent progression in eDNA approaches for elucidating PAI. The current review envisages that eDNA has a significant potential to monitor both plants and associated cohort of prospective pollinators (avian species and flowering plants, bees and plants, arthropods and plants, bats and plants, etc.). Furthermore, a brief description of the factors that influence the utility and interpretation of PAI eDNA is also presented. The review establishes that factors such as biotic and abiotic, primer selection and taxonomic resolution, and indeterminate spatio-temporal scales impact the availability and longevity of eDNA. The study also identified the limitations that influence PAI detection and suggested possible solutions for better execution of these molecular approaches. Overcoming these research caveats will augment the assortment of PAI analysis through eDNA that could be vital for ecosystem health and integrity. This review forms a critical guide and offers prominent insights for ecologists, environmental managers and researchers to assess and evaluate plant-animal interaction through environmental DNA.

An invasive appetite: Combining molecular and stable isotope analyses to reveal the diet of introduced house mice (Mus musculus) on a small, subtropical island

House mice (Mus musculus) pose a conservation threat on islands, where they adversely affect native species' distributions, densities, and persistence. On Sand Island of Kuaihelani, mice recently began to depredate nesting adult mōlī (Laysan Albatross, Phoebastria immutabilis). Efforts are underway to eradicate mice from Sand Island, but knowledge of mouse diet is needed to predict ecosystem response and recovery following mouse removal. We used next-generation sequencing to identify what mice eat on Sand Island, followed by stable isotope analysis to estimate the proportions contributed by taxa to mouse diet. We collected paired fecal and hair samples from 318 mice between April 2018 to May 2019; mice were trapped approximately every eight weeks among four distinct habitat types to provide insight into temporal and spatial variation. Sand Island's mice mainly consume arthropods, with nearly equal (but substantially smaller) contributions of C3 plants, C4 plants, and mōlī. Although seabird tissue is a small portion of mouse diet, mice consume many detrital-feeding arthropods in and around seabird carcasses, such as isopods, flesh flies, ants, and cockroaches. Additionally, most arthropods and plants eaten by mice are non-native. Mouse diet composition differs among habitat types but changes minimally throughout the year, indicating that mice are not necessarily limited by food source availability or accessibility. Eradication of house mice may benefit seabirds on Sand Island (by removing a terrestrial, non-native predator), but it is unclear how arthropod and plant communities may respond and change. Non-native and invasive arthropods and plants previously consumed (and possibly suppressed) by mice may be released post-eradication, which could prevent recovery of native taxa. Comprehensive knowledge of target species' diet is a critical component of eradication planning. Dietary information should be used both to identify and to monitor which taxa may respond most strongly to invasive species removal and to assess if proactive, pre-eradication management activities are warranted.

Fecal metabarcoding of the endangered Pacific pocket mouse (Perognathus longimembris pacificus) reveals a diverse and forb rich diet that reflects local habitat availability

Information on diet breadth and preference can assist in understanding links between food resources and population growth and inform habitat restoration for rare herbivores. We assessed the diet of the endangered Pacific pocket mouse using metabarcoding of fecal samples and compared it to plant community composition in long-term study plots in two populations on Marine Corps Base Camp Pendleton, San Diego County, CA. Fecal samples (n = 221) were collected between spring 2016 and fall 2017 during monthly live-trap surveys. Concurrently, percent cover and plant phenology were measured in plots centered on trap locations. Fecal samples were sequenced with paired-end reads of the internal transcribed spacer 2 region of the nuclear ribosomal gene, and the resulting amplicons were matched to a regionally specific database. Seventy-three plant taxa were detected, which were mostly forbs and perennial herbs (70-90%). Diet composition differed between populations, years, seasons, and plots. Overall, diet and local habitat composition in plots were significantly correlated. However, we detected some differences in above-ground seed availability and proportion in fecal samples that indicate diet preferences for some forbs, perennial herbs, and native bunch grasses over perennial shrubs and non-native grasses. This is the first study of PPM to pair plant phenology surveys with diet metabarcoding to estimate resource selection, and results suggest that managing habitat for diverse native forb communities and reducing non-native grass cover may be beneficial for this critically endangered species.

Potential of plant DNA information in determining the provenance and identify of unknown victims

Determination of the personal identity of victims is particularly important for the settlement of criminal cases. Unfortunately, useful information for identification is not always available. We here propose that the particles (pollens) of some plants with specific geographical distributions extracted from human lung tissues contribute to further determining the provenance or long-term residence of unknown victims, thereby considerably narrowing the search scope of the victims. We collected lung tissues from 155 victims with diverse causes of death, extracted DNA from lung tissues, sequenced the DNA fragments of plants on the Illumina Hiseq platform, and barcoded the plant species using phylogenetic methods. Finally, 108 unique plant sequences were detected in 55 samples and identified to belong to 36 species in 32 genera of 29 families. These plants were predominantly insect-pollinated crops and ornamental plants. No significant difference was observed between male and female samples, between urban and rural samples, or among samples of different ages and different sample sizes. There were 16 samples with 21 wild plant species. The original sources of 15 samples were overlapped with the distribution regions of detected plants; 2 samples narrowed the original sources to 2 provinces, which were quite coincident with their source places; 1 sample had no overlapping with its victim source region. Although plant information was only found in one-third of the samples, we further demonstrated the great potential of plant eDNA in identifying the source of unnamed corpses in a real-world case. We used plant eDNA from lung tissues to explore the provenance of an unknown female corpse found in Beijing. The source place of this victim was narrowed to Guangdong and Guangxi provinces, and finally, we confirmed her true identity in the list of missing persons in Guangxi Province. In the presence of a well-covered local reference library, the plant species detected in the lungs can be accurately identified. In difficult criminal cases where physical evidence is relatively weak, plant DNA information may provide new clues. In conclusion, the plant particles trapped in the lungs are promising to help forensic experts narrow the search scope for the identity of unknown victims.

Introducing an rbcL and a trnL reference library to aid in the metabarcoding analysis of foraged plants from two semi-arid eastern South African savanna bioregions

Success of a metabarcoding study is determined by the extent of taxonomic coverage and the quality of records available in the DNA barcode reference database used. This study aimed to create an rbcL and a trnL (UAA) DNA barcode sequence reference database of plant species that are potential herbivore foraging targets and commonly found in semi-arid savannas of eastern South Africa. An area-specific species list of 765 species was compiled according to plant collection records available and areas comparable to an eastern semi-arid South African savanna. Thereafter, rbcL and trnL sequences of species from this list were mined from GenBank and BOLD sequence databases according to specific quality criteria to ensure accurate taxonomic coverage and resolution. These were supplemented with sequences of 24 species sequenced for this study. A phylogenetic approach, employing Neighbor-Joining, was used to verify the topology of the reference libraries to known angiosperm phylogeny. The taxonomic reliability of these reference libraries was evaluated by testing for the presence of a barcode gap, identifying a data-appropriate identification threshold, and determining the identification accuracy of reference sequences via primary distance-based criteria. The final rbcL reference dataset consisted of 1238 sequences representing 318 genera and 562 species. The final trnL dataset consisted of 921 sequences representing 270 genera and 461 species. Barcode gaps were found for 76% of the taxa in the rbcL barcode reference dataset and 68% of the taxa in the trnL barcode reference dataset. The identification success rate, calculated with the k-nn criterion was 85.86% for the rbcL dataset and 73.72% for the trnL dataset. The datasets for rbcL and trnL combined during this study are not presented as complete DNA reference libraries, but rather as two datasets that should be used in unison to identify plants present in the semi-arid eastern savannas of South Africa.

Best practice for wildlife gut microbiome research: A comprehensive review of methodology for 16S rRNA gene investigations

Extensive research in well-studied animal models underscores the importance of commensal gastrointestinal (gut) microbes to animal physiology. Gut microbes have been shown to impact dietary digestion, mediate infection, and even modify behavior and cognition. Given the large physiological and pathophysiological contribution microbes provide their host, it is reasonable to assume that the vertebrate gut microbiome may also impact the fitness, health and ecology of wildlife. In accordance with this expectation, an increasing number of investigations have considered the role of the gut microbiome in wildlife ecology, health, and conservation. To help promote the development of this nascent field, we need to dissolve the technical barriers prohibitive to performing wildlife microbiome research. The present review discusses the 16S rRNA gene microbiome research landscape, clarifying best practices in microbiome data generation and analysis, with particular emphasis on unique situations that arise during wildlife investigations. Special consideration is given to topics relevant for microbiome wildlife research from sample collection to molecular techniques for data generation, to data analysis strategies. Our hope is that this article not only calls for greater integration of microbiome analyses into wildlife ecology and health studies but provides researchers with the technical framework needed to successfully conduct such investigations.

Using high-throughput sequencing to investigate the dietary composition of the Korean water deer (Hydropotes inermis argyropus): a spatiotemporal comparison

The Korean water deer (Hydropotes inermis argyropus) is considered a vermin in Korea because it damages crops, but also listed as a vulnerable species on the IUCN's red list. Therefore, it is indispensable to manage them appropriately by understanding the ecology such as food habits. Here, we aimed to apply high-throughput sequencing (HTS), a sensitive and objective method, to investigate the dietary composition of the Korean water deer inhabiting the lowland and forest areas in summer and winter. We targeted the internal transcribed spacer 2 (ITS2) region for plant identification. From a total of 40 fecal samples analyzed, 63 plant genera were identified, with Morus being the most abundant, and some of the plant taxa identified by HTS were detected for the first time as the diets of Korean water deer. By type, woody plants (68.6%) were the most predominant, followed by forbs (7.0%) and graminoids (0.7%). We found that the deer in the forest area ate more woody plants (84.6%) than those in the lowland area (52.7%). It was also found that the type of woody plants that the deer ate changed by season. Overall, our results indicate that the Korean water deer is a browser that is seasonally adaptable and feeds on a wide variety of woody plants. We expect that the results and genetics methods reported here, by parallelly investigating their habitat range and reproductive behavior in the future, will help the management and conservation of the Korean water deer, which is in contradictory situations.

Correlation and Influence of Seasonal Variation of Diet with Gut Microbiota Diversity and Metabolism Profile of Chipmunk

Tamias Sibiricus is the only member of the genus Tamias, a significant and vigorous seed distributor and vital food for their predators. No information is known about the strict diet, gut microbiota structure, and metabolism profile of chipmunks and how they diversify seasonally. The above factors, as well as flexibility toward seasonal shifts, are critical in defining its growth rates, health, survivorship, and population stability. This study explored the diet, gut microbiota composition, and chipmunk metabolism. Additionally, the influence of different seasons was also investigated by using next-generation sequencing. Results revealed that seasons strongly affected a diet: streptophyte accounted for 37% in spring, which was lower than in summer (34.3%) and autumn (31.4%). Further, Ascomycota was observed at 43.8% in spring, which reduced to 36.6% in summer and the lowest (31.3%) in autumn. Whereas, nematodes showed maximum abundance from spring (15.8%) to summer (20.6%) and autumn (24.1%). These results signify the insectivorous nature of the chipmunk in summer and autumn. While herbivorous and fungivorous nature in spring. The DNA analysis revealed that chipmunk mainly feeds on fungi, including Aspergillus and Penicillium genus. Similar to diet composition, the microbiome also exhibited highly significant dissimilarity (p < 0.001, R = 0.235) between spring/autumn and spring/summer seasons. Proteobacteria (35.45%), Firmicutes (26.7%), and Bacteroidetes (23.59%) were shown to be the better discriminators as they contributed the most to causing differences between seasons. Moreover, PICRUSt showed that the assimilation of nutrients were also varied seasonally. The abundance of carbohydrates, lipids, nucleotides, xenobiotics, energy, terpenoids, and polyketides metabolism was higher in spring than in other seasons. Our study illustrates that seasonal reconstruction in the chipmunk diet has a significant role in shaping temporal variations in gut microbial community structure and metabolism profile.

Dietary characterization of the endangered salt marsh harvest mouse and sympatric rodents using DNA metabarcoding

The salt marsh harvest mouse (Reithrodontomys raviventris; RERA) is an endangered species endemic to the coastal wetlands of the San Francisco Estuary, California. RERA are specialized to saline coastal wetlands, and their historical range has been severely impacted by landscape conversion and the introduction of non-native plant and rodent species. A better understanding of their diet is needed to assess habitat quality, particularly in relation to potential competitors. We investigated three questions using DNA metabarcoding with ITS2 and trnL markers: (1) Do RERA specialize on the native plant, pickleweed (Salicornia pacifica), (2) Do RERA consume non-native plants, and (3) What is the dietary niche breadth and overlap with three sympatric native and non-native rodents? RERA diet was dominated by two plants, native Salicornia and non-native salt bush (Atriplex spp.), but included 48 plant genera. RERA diet breadth was narrowest in fall, when they consumed the highest frequencies of Salicornia and Atriplex, and broadest in spring, when the frequencies of these two plants were lowest. Diet breadth was slightly lower for RERA than for co-occurring species in pairwise comparisons. All four species consumed similarly high frequencies of wetland plants, but RERA consumed fewer grasses and upland plants, suggesting that it may be less suited to fragmented habitat than sympatric rodents. Diet overlap was lowest between RERA and the native California vole (Microtis californicus). In contrast, RERA diet overlapped substantially with the native western harvest mouse (R. megalotis) and non-native house mouse (Mus musculus), suggesting potential for competition if these species become sufficiently abundant.

Issues of under-representation in quantitative DNA metabarcoding weaken the inference about diet of the tundra vole Microtus oeconomus

During the last decade, methods based on high-throughput sequencing such as DNA metabarcoding have opened up for a range of new questions in animal dietary studies. One of the major advantages of dietary metabarcoding resides in the potential to infer a quantitative relationship between sequence read proportions and biomass of ingested food. However, this relationship's robustness is highly dependent on the system under study, calling for case-specific assessments. Herbivorous small rodents often play important roles in the ecosystem, and the use of DNA metabarcoding for analyses of rodent diets is increasing. However, there has been no direct validation of the quantitative reliability of DNA metabarcoding for small rodents. Therefore, we used an experimental approach to assess the relationship between input plant biomass and sequence reads proportions from DNA metabarcoding in the tundra vole Microtus oeconomus. We found a weakly positive relationship between the number of high-throughput DNA sequences and the expected biomass proportions of food plants. The weak relationship was possibly caused by a systematic under-amplification of one of the three plant taxa fed. Generally, our results add to the growing evidence that case-specific validation studies are required to reliably make use of sequence read abundance as a proxy of relative food proportions in the diet.

Using Fecal DNA Metabarcoding to Investigate Foraging Reveals the Effects of Specific Herbage on the Improved n-3 Fatty Acid (PUFA) Composition in the Longissimus Dorsi Muscle of Grazing Tan Sheep

Understanding the natural diet of grazing sheep can help fulfill their nutritional requirements and positively affect the quality of their meat. Emerging fecal DNA (fDNA) metabarcoding technology can provide more accurate estimates for the dietary composition of free-ranging animals. This study has shown that pasture feeding can promote deposition of n-3 polyunsaturated fatty acids (PUFAs) in Tan lambs' muscle and decrease the ratio of n-6/n-3 fatty acids (FAs), and thus, we investigated the dietary composition of grazing lambs using fDNA metabarcoding to assess the prevalence of medicinal herbage plants in their diet. Herein, based on the full-time natural pasture grazing and 4-h natural pasture grazing with indoor feeding patterns, the herbage taxa (Bassia scoparia, Euphorbia humifusa, Arnebia euchroma, and Salsola sp.) most correlated to n-3 PUFAs were highlighted to elucidate how diversification in dietary components was associated with the muscle FA profile of lambs. Our findings provide experimental evidence for future feeding research.

The Rapid Rise of Next-Generation Natural History

Many ecologists have lamented the demise of natural history and have attributed this decline to a misguided view that natural history is outdated and unscientific. Although there is a perception that the focus in ecology and conservation have shifted away from descriptive natural history research and training toward hypothetico-deductive research, we argue that natural history has entered a new phase that we call “next-generation natural history.” This renaissance of natural history is characterized by technological and statistical advances that aid in collecting detailed observations systematically over broad spatial and temporal extents. The technological advances that have increased exponentially in the last decade include electronic sensors such as camera-traps and acoustic recorders, aircraft- and satellite-based remote sensing, animal-borne biologgers, genetics and genomics methods, and community science programs. Advances in statistics and computation have aided in analyzing a growing quantity of observations to reveal patterns in nature. These robust next-generation natural history datasets have transformed the anecdotal perception of natural history observations into systematically collected observations that collectively constitute the foundation for hypothetico-deductive research and can be leveraged and applied to conservation and management. These advances are encouraging scientists to conduct and embrace detailed descriptions of nature that remain a critically important component of the scientific endeavor. Finally, these next-generation natural history observations are engaging scientists and non-scientists alike with new documentations of the wonders of nature. Thus, we celebrate next-generation natural history for encouraging people to experience nature directly.

Composition and distribution of fish environmental DNA in an Adirondack watershed

Background

Environmental DNA (eDNA) surveys are appealing options for monitoring aquatic biodiversity. While factors affecting eDNA persistence, capture and amplification have been heavily studied, watershed-scale surveys of fish communities and our confidence in such need further exploration.

Methods

We characterized fish eDNA compositions using rapid, low-volume filtering with replicate and control samples scaled for a single Illumina MiSeq flow cell, using the mitochondrial 12S ribosomal RNA locus for taxonomic profiling. Our goals were to determine: (1) spatiotemporal variation in eDNA abundance, (2) the filtrate needed to achieve strong sequencing libraries, (3) the taxonomic resolution of 12S ribosomal sequences in the study environment, (4) the portion of the expected fish community detectable by 12S sequencing, (5) biases in species recovery, (6) correlations between eDNA compositions and catch per unit effort (CPUE) and (7) the extent that eDNA profiles reflect major watershed features. Our bioinformatic approach included (1) estimation of sequencing error from unambiguous mappings and simulation of taxonomic assignment error under various mapping criteria; (2) binning of species based on inferred assignment error rather than by taxonomic rank; and (3) visualization of mismatch distributions to facilitate discovery of distinct haplotypes attributed to the same reference. Our approach was implemented within the St. Regis River, NY, USA, which supports tribal and recreational fisheries and has been a target of restoration activities. We used a large record of St. Regis-specific observations to validate our assignments.

Results

We found that 300 mL drawn through 25-mm cellulose nitrate filters yielded greater than 5 ng/µL DNA at most sites in summer, which was an approximate threshold for generating strong sequencing libraries in our hands. Using inferred sequence error rates, we binned 12S references for 110 species on a state checklist into 85 single-species bins and seven multispecies bins. Of 48 bins observed by capture survey in the St. Regis, we detected eDNA consistent with 40, with an additional four detections flagged as potential contaminants. Sixteen unobserved species detected by eDNA ranged from plausible to implausible based on distributional data, whereas six observed species had no 12S reference sequence. Summed log-ratio compositions of eDNA-detected taxa correlated with log(CPUE) (Pearson’s R = 0.655, P < 0.001). Shifts in eDNA composition of several taxa and a genotypic shift in channel catfish (Ictalurus punctatus) coincided with the Hogansburg Dam, NY, USA. In summary, a simple filtering apparatus operated by field crews without prior expertise gave useful summaries of eDNA composition with minimal evidence of field contamination. 12S sequencing achieved useful taxonomic resolution despite the short marker length, and data exploration with standard bioinformatic tools clarified taxonomic uncertainty and sources of error.

Diet composition analysis provides new management insights for a highly specialized endangered small mammal

The critically endangered Amargosa vole (Microtus californicus scirpensis) is found only in rare marsh habitat near Tecopa, California in a plant community dominated by three-square bulrush (Schoenoplectus americanus). Since the earliest research on the Amargosa vole, the existing paradigm has been that these voles are obligatorily dependent on bulrush as their only food source and for the three-dimensional canopy and litter structure it provides for predator avoidance. However, no prior research has confirmed the diet of the Amargosa vole. In this study we characterized the Amargosa vole' nutritional needs, analyzed the quality of bulrush by forage analysis, and performed microhistological and metabarcoding analyses of vole feces to determine what foods were consumed in the wild. All bulrush plant tissues analyzed were low in fat (from 0.9% of dry matter in roots to 3.6% in seeds), high in neutral detergent fiber (from 5.9% in rhizomes to 33.6% in seeds), and low in protein (7.3-8.4%). These findings support the conclusion that bulrush alone is unlikely to support vole survival and reproduction. Fecal microhistology and DNA metabarcoding revealed relatively diverse diets including plants in 14 families, with rushes (Juncaceae), bulrushes (Cyperaceae), and grasses (Poaceae) being the most common diet items. On microhistology, all analyzed samples contained bulrush, sedges (Carex sp.), rushes (Juncus sp.), and beaked spikerush (Eleocharis rostrellata) even from marshes where non-bulrush plants were uncommon. There was evidence of insects at <1% in two marshes but none in the remaining marshes. Metabarcoding detected ten genera of plants. When considering non-Schoenoplectus targets, for which metabarcoding had poor sensitivity, saltgrass (Distichlis spicata) was the most commonly detected species, with prominent contributions from seaside arrowgrass (Triglochin concinna) and yerba mansa (Anemopsis californica) as well. Diversity of vole diets generally increased with increasing site plant diversity, but differences were not statistically significant. Confirming details about dietary behaviors is critical for informing appropriate conservation planning including habitat management and reintroduction of voles into new sites.

Diet of a rare herbivore based on DNA metabarcoding of feces: Selection, seasonality, and survival

In herbivores, survival and reproduction are influenced by quality and quantity of forage, and hence, diet and foraging behavior are the foundation of an herbivore's life history strategy. Given the importance of diet to most herbivores, it is imperative that we know the species of plants they prefer, especially for herbivorous species that are at risk for extinction. However, it is often difficult to identify the diet of small herbivores because: (a) They are difficult to observe, (b) collecting stomach contents requires sacrificing animals, and (c) microhistology requires accurately identifying taxa from partially digested plant fragments and likely overemphasizes less-digestible taxa. The northern Idaho ground squirrel (Urocitellus brunneus) is federally threatened in the United States under the Endangered Species Act. We used DNA metabarcoding techniques to identify the diet of 188 squirrels at 11 study sites from fecal samples. We identified 42 families, 126 genera, and 120 species of plants in the squirrel's diet. Our use of three gene regions was beneficial because reliance on only one gene region (e.g., only trnL) would have caused us to miss >30% of the taxa in their diet. Northern Idaho ground squirrel diet differed between spring and summer, frequency of many plants in the diet differed from their frequency within their foraging areas (evidence of selective foraging), and several plant genera in their diet were associated with survival. Our results suggest that while these squirrels are generalists (they consume a wide variety of plant species), they are also selective and do not eat plants relative to availability. Consumption of particular genera such as Perideridia may be associated with higher overwinter survival.

Dual-locus DNA metabarcoding reveals southern hairy-nosed wombats (Lasiorhinus latifrons Owen) have a summer diet dominated by toxic invasive plants

Habitat degradation and summer droughts severely restrict feeding options for the endangered southern hairy-nosed wombat (SHNW; Lasiorhinus latifrons). We reconstructed SHNW summer diets by DNA metabarcoding from feces. We initially validated rbcL and ndhJ diet reconstructions using autopsied and captive animals. Subsequent diet reconstructions of wild wombats broadly reflected vegetative ground cover, implying local rather than long-range foraging. Diets were all dominated by alien invasives. Chemical analysis of alien food revealed Carrichtera annua contains high levels of glucosinolates. Clinical examination (7 animals) and autopsy (12 animals) revealed that the most degraded site also contained most individuals showing signs of glucosinolate poisoning. We infer that dietary poisoning through the ingestion of alien invasives may have contributed to the recent population crashes in the region. In floristically diverse sites, individuals appear to be able to manage glucosinolate intake by avoidance or episodic feeding but this strategy is less tractable in the most degraded sites. We conclude that recovery of the most affected populations may require effective Carrichtera management and interim supplementary feeding. More generally, we argue that protection against population decline by poisoning in territorial herbivores requires knowledge of their diet and of those food plants containing toxic principles.

Biostatistical Options for Quantitative Diet Analysis

Sufficient statistics knowledge is crucial for the correct design of a research plan. The elucidations of results are interpretive only if appropriate statistical methods are applied. Statistical strategies are a particular approach to demonstrate complicated information in broad and explicable conclusions. The emergence of biostatistical approaches for diet evaluation has improved the accuracy of diet estimation, and different methodologies of data integration promise to magnify our understanding of ecological communities. The present study aimed to compile multiple statistical methods used for diet analysis. More specifically, the significant analysis used in diet assessment, central expectations, and preferences related to each measure was conceptualized. In addition, the ability of each test to evaluate diversity, richness, differentiation, fluctuation, similarity, and quantification of multiple diet items was summarized. Moreover, different options were proposed for researchers to select the appropriate statistical tests. This study covers a framework, aim, and understanding of the statistical test methods of diet analysis.

How quantitative is metabarcoding: A meta-analytical approach

Metabarcoding has been used in a range of ecological applications such as taxonomic assignment, dietary analysis and the analysis of environmental DNA. However, after a decade of use in these applications there is little consensus on the extent to which proportions of reads generated corresponds to the original proportions of species in a community. To quantify our current understanding, we conducted a structured review and meta‐analysis. The analysis suggests that a weak quantitative relationship may exist between the biomass and sequences produced (slope = 0.52 ± 0.34, p < 0.01), albeit with a large degree of uncertainty. None of the tested moderators, sequencing platform type, the number of species used in a trial or the source of DNA, were able to explain the variance. Our current understanding of the factors affecting the quantitative performance of metabarcoding is still limited: additional research is required before metabarcoding can be confidently utilized for quantitative applications. Until then, we advocate the inclusion of mock communities when metabarcoding as this facilitates direct assessment of the quantitative ability of any given study.

New universal ITS2 primers for high-resolution herbivory analyses using DNA metabarcoding in both tropical and temperate zones

DNA metabarcoding is a rapidly growing technique for obtaining detailed dietary information. Current metabarcoding methods for herbivory, using a single locus, can lack taxonomic resolution for some applications. We present novel primers for the second internal transcribed spacer of nuclear ribosomal DNA (ITS2) designed for dietary studies in Mauritius and the UK, which have the potential to give unrivalled taxonomic coverage and resolution from a short-amplicon barcode. In silico testing used three databases of plant ITS2 sequences from UK and Mauritian floras (native and introduced) totalling 6561 sequences from 1790 species across 174 families. Our primers were well-matched in silico to 88% of species, providing taxonomic resolution of 86.1%, 99.4% and 99.9% at the species, genus and family levels, respectively. In vitro, the primers amplified 99% of Mauritian (n = 169) and 100% of UK (n = 33) species, and co-amplified multiple plant species from degraded faecal DNA from reptiles and birds in two case studies. For the ITS2 region, we advocate taxonomic assignment based on best sequence match instead of a clustering approach. With short amplicons of 187-387 bp, these primers are suitable for metabarcoding plant DNA from faecal samples, across a broad geographic range, whilst delivering unparalleled taxonomic resolution.

Dropping Hints: Estimating the diets of livestock in rangelands using DNA metabarcoding of faeces

The introduction of domesticated animals into new environments can lead to considerable ecological disruption, and it can be difficult to predict their impact on the new ecosystem. In this study, we use faecal metabarcoding to characterize the diets of three ruminant taxa in the rangelands of south-western New South Wales, Australia. Our study organisms included goats (Capraaegagrushircus) and two breeds of sheep (Ovisaries): Merinos, which have been present in Australia for over two hundred years, and Dorpers, which were introduced in the 1990s. We used High-Throughput Sequencing methods to sequence therbcLand ITS2 genes of plants in the faecal samples, and identified the samples using the GenBank and BOLD online databases, as well as a reference collection of sequences from plants collected in the study area. We found that the diets of all three taxa were dominated by the family Malvaceae, and that the Dorper diet was more diverse than the Merino diet at both the family and the species level. We conclude that Dorpers, like Merinos, are potentially a threat to some vulnerable species in the rangelands of New South Wales.

Metabarcoding for the parallel identification of several hundred predators and their prey: Application to bat species diet analysis

Assessing diet variability is of main importance to better understand the biology of bats and design conservation strategies. Although the advent of metabarcoding has facilitated such analyses, this approach does not come without challenges. Biases may occur throughout the whole experiment, from fieldwork to biostatistics, resulting in the detection of false negatives, false positives or low taxonomic resolution. We detail a rigorous metabarcoding approach based on a short COI minibarcode and two-step PCR protocol enabling the "all at once" taxonomic identification of bats and their arthropod prey for several hundreds of samples. Our study includes faecal pellets collected in France from 357 bats representing 16 species, as well as insect mock communities that mimic bat meals of known composition, negative and positive controls. All samples were analysed using three replicates. We compare the efficiency of DNA extraction methods, and we evaluate the effectiveness of our protocol using identification success, taxonomic resolution, sensitivity and amplification biases. Our parallel identification strategy of predators and prey reduces the risk of mis-assigning prey to wrong predators and decreases the number of molecular steps. Controls and replicates enable to filter the data and limit the risk of false positives, hence guaranteeing high confidence results for both prey occurrence and bat species identification. We validate 551 COI variants from arthropod including 18 orders, 117 family, 282 genus and 290 species. Our method therefore provides a rapid, resolutive and cost-effective screening tool for addressing evolutionary ecological issues or developing "chirosurveillance" and conservation strategies.