New mitochondrial primers for metabarcoding of insects,designed and evaluated using in silico methods

Insect metabarcoding has been mainly based on PCR amplification of short fragments within the “barcoding region” of the gene cytochrome oxidase I (COI). However, because of the variability of this gene, it has been difficult to design good universal PCR primers. Most primers used today are associated with gaps in the taxonomic coverage or amplification biases that make the results less reliable and impede the detection of species that are present in the sample. We identify new primers for insect metabarcoding using computational approaches (ecoprimers and degeprime ) applied to the most comprehensive reference databases of mitochondrial genomes of Hexapoda assembled to date. New primers are evaluated in silico against previously published primers in terms of taxonomic coverage and resolution of the corresponding amplicons. For the latter criterion, we propose a new index, exclusive taxonomic resolution, which is a more biologically meaningful measure than the standard index used today. Our results show that the best markers are found in the ribosomal RNA genes (12S and 16S); they resolve about 90% of the genetically distinct species in the reference database. Some markers in protein‐coding genes provide similar performance but only at much higher levels of primer degeneracy. Combining two of the best individual markers improves the effective taxonomic resolution with up to 10%. The resolution is strongly dependent on insect taxon: COI primers detect 40% of Hymenoptera, while 12S primers detect 12% of Collembola. Our results indicate that amplicon‐based metabarcoding of insect samples can be improved by choosing other primers than those commonly used today.     

Environmental DNA metabarcoding primers for freshwater fish detection and quantification: In silico and in tanks

Environmental DNA (eDNA) techniques refer to utilizing the organisms’ DNA extracted from environment samples to genetically identify target species without capturing actual organisms. eDNA metabarcoding via high‐throughput sequencing can simultaneously detect multiple fish species from a single water sample, which is a powerful tool for the qualitative detection and quantitative estimates of multiple fish species. However, sequence counts obtained from eDNA metabarcoding may be influenced by many factors, of which primer bias is one of the foremost causes of methodological error. The performance of 18 primer pairs for COI, cytb, 12S rRNA, and 16S rRNA mitochondrial genes, which are all frequently used in fish eDNA metabarcoding, were evaluated in the current study. The ribosomal gene markers performed better than the protein‐coding gene markers during in silico screening, resulting in higher taxonomic coverage and appropriate barcode lengths. Four primer pairs—AcMDB07, MiFish‐U, Ve16S1, and Ve16S3—designed for various regions of the 12S and 16S rRNA genes were screened for tank metabarcoding in a case study targeting six freshwater fish species. The four primer pairs were able to accurately detect all six species in different tanks, while only MiFish‐U, Ve16S1, and Ve16S3 revealed a significant positive relationship between species biomass and read count for the pooled tank data. The positive relationship could not be found in all species within the tanks. Additionally, primer efficiency differed depending on the species while primer preferential species varied in different fish assemblages. This case study supports the potential for eDNA metabarcoding to assess species diversity in natural ecosystems and provides an alternative strategy to evaluate the performance of candidate primers before application of eDNA metabarcoding in natural ecosystems.     

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.     

Quantitative DNA metabarcoding: improved estimates of species proportional biomass using correction factors derived from control material

DNA metabarcoding is a powerful new tool allowing characterization of species assemblages using high‐throughput amplicon sequencing. The utility of DNA metabarcoding for quantifying relative species abundances is currently limited by both biological and technical biases which influence sequence read counts. We tested the idea of sequencing 50/50 mixtures of target species and a control species in order to generate relative correction factors (RCFs) that account for multiple sources of bias and are applicable to field studies. RCFs will be most effective if they are not affected by input mass ratio or co‐occurring species. In a model experiment involving three target fish species and a fixed control, we found RCFs did vary with input ratio but in a consistent fashion, and that 50/50 RCFs applied to DNA sequence counts from various mixtures of the target species still greatly improved relative abundance estimates (e.g. average per species error of 19 ± 8% for uncorrected vs. 3 ± 1% for corrected estimates). To demonstrate the use of correction factors in a field setting, we calculated 50/50 RCFs for 18 harbour seal (Phoca vitulina) prey species (RCFs ranging from 0.68 to 3.68). Applying these corrections to field‐collected seal scats affected species percentages from individual samples (Δ 6.7 ± 6.6%) more than population‐level species estimates (Δ 1.7 ± 1.2%). Our results indicate that the 50/50 RCF approach is an effective tool for evaluating and correcting biases in DNA metabarcoding studies. The decision to apply correction factors will be influenced by the feasibility of creating tissue mixtures for the target species, and the level of accuracy needed to meet research objectives.     

Quantitative and qualitative assessment of pollen DNA metabarcoding using constructed species mixtures

Pollen DNA metabarcoding—marker‐based genetic identification of potentially mixed‐species pollen samples—has applications across a variety of fields. While basic species‐level pollen identification using standard DNA barcode markers is established, the extent to which metabarcoding (a) correctly assigns species identities to mixes (qualitative matching) and (b) generates sequence reads proportionally to their relative abundance in a sample (quantitative matching) is unclear, as these have not been assessed relative to known standards. We tested the quantitative and qualitative robustness of metabarcoding in constructed pollen mixtures varying in species richness (1–9 species), taxonomic relatedness (within genera to across class) and rarity (5%–100% of grains), using Illumina MiSeq with the markers rbcL and ITS2. Qualitatively, species composition determinations were largely correct, but false positives and negatives occurred. False negatives were typically driven by lack of a barcode gap or rarity in a sample. Species richness and taxonomic relatedness, however, did not strongly impact correct determinations. False positives were likely driven by contamination, chimeric sequences and/or misidentification by the bioinformatics pipeline. Quantitatively, the proportion of reads for each species was only weakly correlated with its relative abundance, in contrast to suggestions from some other studies. Quantitative mismatches are not correctable by consistent scaling factors, but instead are context‐dependent on the other species present in a sample. Together, our results show that metabarcoding is largely robust for determining pollen presence/absence but that sequence reads should not be used to infer relative abundance of pollen grains.     

Toward an ecoregion scale evaluation of eDNA metabarcoding primers: A case study for the freshwater fish biodiversity of the Murray–Darling Basin (Australia)

High‐throughput sequencing of environmental DNA (i.e., eDNA metabarcoding) has become an increasingly popular method for monitoring aquatic biodiversity. At present, such analyses require target‐specific primers to amplify DNA barcodes from co‐occurring species, and this initial amplification can introduce biases. Understanding the performance of different primers is thus recommended prior to undertaking any metabarcoding initiative. While multiple software programs are available to evaluate metabarcoding primers, all programs have their own strengths and weaknesses. Therefore, a robust in silico workflow for the evaluation of metabarcoding primers will benefit from the use of multiple programs. Furthermore, geographic differences in species biodiversity are likely to influence the performance of metabarcoding primers and further complicate the evaluation process. Here, an in silico workflow is presented that can be used to evaluate the performance of metabarcoding primers on an ecoregion scale. This workflow was used to evaluate the performance of published and newly developed eDNA metabarcoding primers for the freshwater fish biodiversity of the Murray–Darling Basin (Australia). To validate the in silico workflow, a subset of the primers, including one newly designed primer pair, were used in metabarcoding analyses of an artificial DNA community and eDNA samples. The results show that the in silico workflow allows for a robust evaluation of metabarcoding primers and can reveal important trade‐offs that need to be considered when selecting the most suitable primer. Additionally, a new primer pair was described and validated that allows for more robust taxonomic assignments and is less influenced by primer biases compared to commonly used fish metabarcoding primers.     

Comparison of destructive and nondestructive DNA extraction methods for the metabarcoding of arthropod bulk samples

DNA metabarcoding is routinely used for biodiversity assessment, in particular targeting highly diverse groups for which limited taxonomic expertise is available. Various protocols are currently in use, although standardization is key to its application in large-scale monitoring. DNA metabarcoding of arthropod bulk samples can be conducted either destructively from sample tissue, or nondestructively from sample fixative or lysis buffer. Nondestructive methods are highly desirable for the preservation of sample integrity but have yet to be experimentally evaluated in detail. Here, we compare diversity estimates from 14 size-sorted Malaise trap samples processed consecutively with three nondestructive approaches (one using fixative ethanol and two using lysis buffers) and one destructive approach (using homogenized tissue). Extraction from commercial lysis buffer yielded comparable species richness and high overlap in species composition to the ground tissue extracts. A significantly divergent community was detected from preservative ethanol-based DNA extraction. No consistent trend in species richness was found with increasing incubation time in lysis buffer. These results indicate that nondestructive DNA extraction from incubation in lysis buffer could provide a comparable alternative to destructive approaches with the added advantage of preserving the specimens for postmetabarcoding taxonomic work but at a higher cost per sample.     

A new approach to molecular biosurveillance of invasive species using DNA metabarcoding

Non‐indigenous species (NIS) reach every corner of the world, at times wreaking havoc on ecosystems and costing the global economy billions of dollars. A rapid and accurate biosurveillance tool tailored to a particular biogeographic region is needed to detect NIS when they are first introduced into an area as traditional detection methods are expensive and require specialized expertise. Metabarcoding of environmental and community DNA meets those biosurveillance requirements; a novel tool tailored to the Northwest Pacific Ocean is presented here using an approach that could revolutionize early detection of NIS. Eight newly designed genetic markers for multiple gene regions were implemented to meet the stringent taxonomic requirements for the detection of NIS across four major marine phyla. The tool was considered highly successful because it identified 12 known NIS in the study area and a further seven species representing potential new records. Overall community composition detected here was statistically different between substrate types; zooplankton sampling accounted for significantly higher species richness than filtered sea water in most cases, but this was dominated by mollusk and arthropod species. Both substrate types sampled were required to identify the wide taxonomic breadth of known NIS in the study area. Intensive sampling is known to be paramount for the detection of rare species, including new incursions of NIS, thus it is recommended to include diverse DNA sampling protocols based on species’ life‐history characteristics for broad detection capacity. Application of a metabarcoding‐based molecular biosurveillance tool optimized for biogeographic regions enables rapid and accurate early detection across a wide taxonomic range to allow quick implementation of eradication or control efforts and potentially mitigate some of the devastating effects of NIS worldwide.     

Characterization of RNA primers synthesized by the human breast cancer cell DNA synthesome

We previously reported on the purification and characterization of a functional multi‐protein DNA replication complex (the DNA synthesome) from human cells and tissues. The synthesome is fully competent to carry‐out all phases of the DNA replication process in vitro. In this study, DNA primase, a component of the synthesome, is examined to determine its activity and processivity in the in vitro synthesis and extension of RNA primers. Our results show that primase activity in the P4 fraction of the synthesome is 30‐fold higher than that of crude cell extracts. The synthesome synthesizes RNA primers that are 7–10 ribonucleotides long and DNA primers that are 20–40 deoxyribonucleotides long using a poly(dT) template of exogenous single‐stranded DNA. The synthesome‐catalyzed RNA primers can be elongated by E. coli DNA polymerase I to form the complementary DNA strands on the poly(dT) template. In addition, the synthesome also supports the synthesis of native RNA primers in vitro using an endogenous supercoiled double‐stranded DNA template. Gel analysis demonstrates that native RNA primers are oligoribonucleotides of 10–20 nt in length and the primers are covalently link to DNA to form RNA‐primed nascent DNA of 100–200 nt. Our study reveals that the synthesome model is capable of priming and continuing DNA replication. The ability of the synthesome to synthesize and extend RNA primers in vitro elucidates the organizational and functional properties of the synthesome as a potentially useful replication apparatus to study the function of primase and the interaction of primase with other replication proteins. J. Cell. Biochem. 106: 798–811, 2009. © 2009 Wiley‐Liss, Inc.     

Identification of the pathogenic Aspergillus species by nested PCR using a mixture of specific primers to DNA topoisomerase II gene

For PCR-based identification of Aspergillus species, a common primer of the DNA topoisomerase II genes of Candida, Aspergillus and Penicillium, and species-specific primers of the genomic sequences of DNA topoisomerase II of A. fumigatus, A. niger, A. flavus (A. oryzae), A. nidulans and A. terreus were tested for their specificities in PCR amplifications. The method consisted of amplification of the genomic DNA topoisomerase II gene by a common primer set, followed by a second PCR with a primer mix consisting of 5 species-specific primer pairs for each Aspergillus species. By using the common primer pair, a DNA fragment of approximately 1,200 bp was amplified from the Aspergillus and Penicillium genomic DNAs. Using each species-specific primer pair, unique sizes of PCR products were amplified, all of which corresponded to a species of Aspergillus even in the presence of DNAs of several fungal species. The sensitivity of A. fumigatus to the nested PCR was found to be 100 fg of DNA in the reaction mixture. In the nested PCR obtained by using the primer mix (PsIV), the specific DNA fragment of A. fumigatus was amplified from clinical specimens. These results suggest that this nested PCR method is rapid, simple and available as a tool for identification of pathogenic Aspergillus to a species level.     

Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.     

Quantification of mesocosm fish and amphibian species diversity via environmental DNA metabarcoding

Freshwater fauna are particularly sensitive to environmental change and disturbance. Management agencies frequently use fish and amphibian biodiversity as indicators of ecosystem health and a way to prioritize and assess management strategies. Traditional aquatic bioassessment that relies on capture of organisms via nets, traps and electrofishing gear typically has low detection probabilities for rare species and can injure individuals of protected species. Our objective was to determine whether environmental DNA (eDNA) sampling and metabarcoding analysis can be used to accurately measure species diversity in aquatic assemblages with differing structures. We manipulated the density and relative abundance of eight fish and one amphibian species in replicated 206‐L mesocosms. Environmental DNA was filtered from water samples, and six mitochondrial gene fragments were Illumina‐sequenced to measure species diversity in each mesocosm. Metabarcoding detected all nine species in all treatment replicates. Additionally, we found a modest, but positive relationship between species abundance and sequencing read abundance. Our results illustrate the potential for eDNA sampling and metabarcoding approaches to improve quantification of aquatic species diversity in natural environments and point the way towards using eDNA metabarcoding as an index of macrofaunal species abundance.     

Design of phylum-specific hybrid primers for DNA barcoding: addressing the need for efficient COI amplification in the Echinodermata

Recent research has shown the usefulness of the Folmer region of the cytochrome oxidase I (COI) as a genetic barcode to assist in species delimitation of echinoderms. However, amplification of COI is often challenging in echinoderms (low success or pseudogenes). We present a method that allows the design of phylum-specific hybrid primers, and use this to develop COI primers for the Echinodermata. We aligned COI sequences from 310 echinoderm species and designed all possible primers along the consensus sequence with two methods (standard degenerate and hybrid). We found much lower degeneracy for hybrid primers (4-fold degeneracy) than for standard degenerate primers (≥48-fold degeneracy). We then designed the most conserved hybrid primers to amplify a >500-bp region within COI. These primers successfully amplified this gene region in all tested taxa (123 species across all echinoderm classes). Sequencing of 30 species among these confirmed both the quality of the sequences (>500 bp, no pseudogenes) and their utility as a DNA barcode. This method should be useful for developing primers for other mitochondrial genes and other phyla. The method will also be of interest for the development of future projects involving both community-based genetic assessments on macroorganisms and biodiversity assessment of environmental samples using high-throughput sequencing.     

Meeting the challenge of DNA barcoding Neotropical amphibians: polymerase chain reaction optimization and new COI primers

Amphibians are one of the most threatened vertebrate classes, yet at the same time new species are being described every year, demonstrating that the number of existing species is grossly underestimated. In groups such as amphibians, with high extinction rates and poorly known species boundaries, DNA barcoding is a tool that can rapidly assess genetic diversity and estimate species richness for prioritizing conservation decisions. However, reliable recovery of the 5′ region of the cytochrome c oxidase subunit 1 (COI) gene is critical for the ongoing effort to gather DNA barcodes for all amphibian species. Here, we provide new PCR conditions and tested new primers that increase the efficiency of barcode recovery in amphibians. We found that a low extension temperature for PCR cycles significantly improves the efficiency of amplification for all combinations of primers. Combining low PCR extension temperature and primers AnF1 + AnR1, we were able to recover COI sequences for 100% of the species analysed (N = 161), encompassing ~15% of the species known from Brazil (representing 77 genera and 23 families), which is an important improvement over previous studies. The preliminary assessment of species diversity suggested that number of species might be underestimated by about 25%. We conclude that DNA barcoding is an efficient, simple, and standardized protocol for identifying cryptic diversity in amphibians and advocate for its use in biodiversity inventories and across widespread populations within known species.     

Focal vs. fecal: Seasonal variation in the diet of wild vervet monkeys from observational and DNA metabarcoding data

Assessing the diet of wild animals reveals valuable information about their ecology and trophic relationships that may help elucidate dynamic interactions in ecosystems and forecast responses to environmental changes. Advances in molecular biology provide valuable research tools in this field. However, comparative empirical research is still required to highlight strengths and potential biases of different approaches. Therefore, this study compares environmental DNA and observational methods for the same study population and sampling duration. We employed DNA metabarcoding assays targeting plant and arthropod diet items in 823 fecal samples collected over 12 months in a wild population of an omnivorous primate, the vervet monkey (Chlorocebus pygerythrus). DNA metabarcoding data were subsequently compared to direct observations. We observed the same seasonal patterns of plant consumption with both methods; however, DNA metabarcoding showed considerably greater taxonomic coverage and resolution compared to observations, mostly due to the construction of a local plant DNA database. We found a strong effect of season on variation in plant consumption largely shaped by the dry and wet seasons. The seasonal effect on arthropod consumption was weaker, but feeding on arthropods was more frequent in spring and summer, showing overall that vervets adapt their diet according to available resources. The DNA metabarcoding assay outperformed also direct observations of arthropod consumption in both taxonomic coverage and resolution. Combining traditional techniques and DNA metabarcoding data can therefore not only provide enhanced assessments of complex diets and trophic interactions to the benefit of wildlife conservationists and managers but also opens new perspectives for behavioral ecologists studying whether diet variation in social species is induced by environmental differences or might reflect selective foraging behaviors.     

Comparing whole‐genome shotgun sequencing and DNA metabarcoding approaches for species identification and quantification of pollen species mixtures

Molecular identification of mixed‐species pollen samples has a range of applications in various fields of research. To date, such molecular identification has primarily been carried out via amplicon sequencing, but whole‐genome shotgun (WGS) sequencing of pollen DNA has potential advantages, including (1) more genetic information per sample and (2) the potential for better quantitative matching. In this study, we tested the performance of WGS sequencing methodology and publicly available reference sequences in identifying species and quantifying their relative abundance in pollen mock communities. Using mock communities previously analyzed with DNA metabarcoding, we sequenced approximately 200Mbp for each sample using Illumina HiSeq and MiSeq. Taxonomic identifications were based on the Kraken k‐mer identification method with reference libraries constructed from full‐genome and short read archive data from the NCBI database. We found WGS to be a reliable method for taxonomic identification of pollen with near 100% identification of species in mixtures but generating higher rates of false positives (reads not identified to the correct taxon at the required taxonomic level) relative to rbcL and ITS2 amplicon sequencing. For quantification of relative species abundance, WGS data provided a stronger correlation between pollen grain proportion and sequence read proportion, but diverged more from a 1:1 relationship, likely due to the higher rate of false positives. Currently, a limitation of WGS‐based pollen identification is the lack of representation of plant diversity in publicly available genome databases. As databases improve and costs drop, we expect that eventually genomics methods will become the methods of choice for species identification and quantification of mixed‐species pollen samples.     

Polymorphism of the SOD1-DNA aggregation species can be modulated by DNA

Proteinaceous aggregates rich in copper, zinc superoxide dismutase (SOD1) have been found in both in vivo and in vitro models. We have shown that double-stranded DNA that acts as a template accelerates the in vitro formation of wild-type SOD1 aggregates. Here, we examined the polymorphism of templated-SOD1 aggregates generated in vitro upon association with DNA under different conditions. Electron microscopy imaging indicates that this polymorphism is capable of being manipulated by the shapes, structures, and doses of the DNAs tested. The nanometer- and micrometer-scale aggregates formed under acidic conditions and under neutral conditions containing ascorbate fall into three classes: aggregate monomers, oligomeric aggregates, and macroaggregates. The aggregate monomers observed at given DNA doses exhibit a polymorphism that is markedly corresponded to the coiled shapes of linear DNA and structures of plasmid DNA. On the other hand, the regularly branched structures observed under both atomic force microscopy and optical microscope indicate that the DNAs tested are simultaneously condensed into a nanoparticle with a specific morphology during SOD1 aggregation, revealing that SOD1 aggregation and DNA condensation are two concurrent phenomena. The results might provide the basis of therapeutic approaches to suppress the formation of toxic protein oligomers or aggregates by screening the toxicity of the protein aggregates with various sizes and morphologies.     

Novel polymerase chain reaction primers for the specific detection of bacterial copper P‐type ATPases gene sequences in environmental isolates and metagenomic DNA

Aims: In the last decades, the worldwide increase in copper wastes release by industrial activities like mining has driven environmental metal contents to toxic levels. For this reason, the study of the biological copper‐resistance mechanisms in natural environments is important. Therefore, an appropriate molecular tool for the detection and tracking of copper‐resistance genes was developed. Methods and Results: In this work, we designed a PCR primer pair to specifically detect copper P‐type ATPases gene sequences. These PCR primers were tested in bacterial isolates and metagenomic DNA from intertidal marine environments impacted by copper pollution. As well, T‐RFLP fingerprinting of these gene sequences was used to compare the genetic composition of such genes in microbial communities, in normal and copper‐polluted coastal environments. New copper P‐type ATPases gene sequences were found, and a high degree of change in the genetic composition because of copper exposure was also determined. Conclusions: This PCR based method is useful to track bacterial copper‐resistance gene sequences in the environment. Significance and Impact of the Study: This study is the first to report the design and use of a PCR primer pair as a molecular marker to track bacterial copper‐resistance determinants, providing an excellent tool for long‐term analysis of environmental communities exposed to metal pollution.     

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.     

Using DNA metabarcoding for simultaneous inference of common vampire bat diet and population structure

Metabarcoding diet analysis has become a valuable tool in animal ecology; however, co‐amplified predator sequences are not generally used for anything other than to validate predator identity. Exemplified by the common vampire bat, we demonstrate the use of metabarcoding to infer predator population structure alongside diet assessments. Growing populations of common vampire bats impact human, livestock and wildlife health in Latin America through transmission of pathogens, such as lethal rabies viruses. Techniques to determine large‐scale variation in vampire bat diet and bat population structure would empower locality‐ and species‐specific projections of disease transmission risks. However, previously used methods are not cost‐effective and efficient for large‐scale applications. Using bloodmeal and faecal samples from common vampire bats from coastal, Andean and Amazonian regions of Peru, we showcase metabarcoding as a scalable tool to assess vampire bat population structure and feeding preferences. Dietary metabarcoding was highly effective, detecting vertebrate prey in 93.2% of the samples. Bats predominantly preyed on domestic animals, but fed on tapirs at one Amazonian site. In addition, we identified arthropods in 9.3% of samples, likely reflecting consumption of ectoparasites. Using the same data, we document mitochondrial geographic population structure in the common vampire bat in Peru. Such simultaneous inference of vampire bat diet and population structure can enable new insights into the interplay between vampire bat ecology and disease transmission risks. Importantly, the methodology can be incorporated into metabarcoding diet studies of other animals to couple information on diet and population structure.