A fecal sequel: Testing the limits of a genetic assay for bat species identification

Dietary niche partitioning of three Sky Island Sceloporus lizards as revealed through DNA metabarcoding

Lizard diets are highly diverse and have contributed to the diversification, biogeographical distributions, and evolution of novel traits across this global radiation. Many parts of a lizard's ecology-including habitat preferences, foraging modes, predation risks, interspecific competition, and thermal constraints, among others-interact to shape diets, and dietary niche partitioning simultaneously contributes to co-occurrence within communities. We used DNA metabarcoding of fecal samples to identify prey items in the diets of three sympatric Sceloporus lizards in the Madrean Sky Islands of Arizona, USA. We found evidence for dietary niche partitioning between interacting species concomitant with their respective ecologies. We also compared diet composition between populations to understand how conserved or plastic species' diets are between different environments. Our findings suggest that habitat generalists are also diet generalists in this system, while the same may be true for specialists. The identification of prey items to much lower taxonomic levels than previously documented further reveals hidden diversity in the diets of these species and underscores the utility of metabarcoding for understanding the full complexity of lizard diets.

THAPBI PICT-a fast, cautious, and accurate metabarcoding analysis pipeline

THAPBI PICT is an open source software pipeline for metabarcoding analysis of Illumina paired-end reads, including cases of multiplexing where more than one amplicon is amplified per DNA sample. Initially a Phytophthora ITS1 Classification Tool (PICT), we demonstrate using worked examples with our own and public data sets how, with appropriate primer settings and a custom database, it can be applied to other amplicons and organisms, and used for reanalysis of existing datasets. The core dataflow of the implementation is (i) data reduction to unique marker sequences, often called amplicon sequence variants (ASVs), (ii) dynamic thresholds for discarding low abundance sequences to remove noise and artifacts (rather than error correction by default), before (iii) classification using a curated reference database. The default classifier assigns a label to each query sequence based on a database match that is either perfect, or a single base pair edit away (substitution, deletion or insertion). Abundance thresholds for inclusion can be set by the user or automatically using per-batch negative or synthetic control samples. Output is designed for practical interpretation by non-specialists and includes a read report (ASVs with classification and counts per sample), sample report (samples with counts per species classification), and a topological graph of ASVs as nodes with short edit distances as edges. Source code available from https://github.com/peterjc/thapbi-pict/ with documentation including installation instructions.

Our good neighbors: Understanding ecosystem services provided by insectivorous bats in Rwanda

Bats are prodigious consumers of agricultural and forest pests, and are, therefore, a natural asset for agricultural productivity, suppressing populations of such pests. This study provides baseline information of diet of 143 bats belonging to eight insectivorous bat species from agricultural areas of Rwanda while evaluating the effectiveness of bats as pest suppressors. Using DNA metabarcoding to analyze bat fecal pellets, 85 different insect species were detected, with 60% (n = 65), 64% (n = 11) and 78% (n = 9) found to be agricultural pests from eastern, northern and western regions, respectively. Given the high percentages of agricultural pests detected, we submit that Rwandan insectivorous bats have the capacity for biocontrol of agricultural pests. Rwandan bat populations should be protected and promoted since they may foster higher crop yields and sustainable livelihoods.

Bat Use of Hollows in California’s Old-Growth Redwood Forests: From DNA to Ecology

The loss of roosting resources, either through disturbance or removal, negatively affects bats. Identifying sensitive species and determining roost requirements are critical components in conserving their habitat. Cavity-roosting bats on the North Coast of California are known to use hollows in large redwood trees. In this study, we examined the factors determining the use of basal tree hollows by different bat species at eight redwood forest sites in Del Norte, Humboldt, and Mendocino Counties, California. Bat guano was collected from 179 basal hollow roosts from 2017 to 2018, and guano mass was used as an index of roosting activity. Nine bat species and one species group were identified by analysis of DNA in guano. We made a total of 253 identifications from 83 hollows into the 10 species categories. The most prevalent species were Myotis californicus (California myotis; 28.5% of all identifications), the Myotis evotis-Myotis thysanodes group (17.4%), Corynorhinus townsendii (17.0%), and Myotis volans (15.0%). We evaluated the extent to which habitat variables at the scales of the hollow, vicinity, and site influenced the level of roost use. The correlations between guano mass and habitat variables were examined using generalized additive mixed models. At the hollow scale, guano mass increased with ceiling height above the opening. At the vicinity scale, guano mass increased with less cover of small trees. At the site scale, there was no association between guano mass and distance to foraging areas, elevation, or the number of nearby hollows. These tree hollow roost preferences can inform land managers when planning the management and conservation of redwood forests.

Multifaceted DNA Metabarcoding of Guano to Uncover Multiple Classes of Ecological Data in Two Different Bat Communities

DNA contained in animal scat provides a wealth of information about the animal, and DNA metabarcoding of scat collections can provide key information about animal populations and communities. Next‐generation DNA sequencing technologies and DNA metabarcoding provide an efficient means for obtaining information available in scat samples. We used multifaceted DNA metabarcoding (MDM) of noninvasively collected bat guano pellets from a Myotis lucifugus colony on Fort Drum Military Installation, New York, USA, and from two mixed‐species bat roosts on Fort Huachuca Military Installation, Arizona, USA, to identify attributes such as bat species composition, sex ratios, diet, and the presence of pathogens and parasites. We successfully identified bat species for nearly 98% of samples from Fort Drum and 90% of samples from Fort Huachuca, and identified the sex for 84% and 67% of samples from these same locations, respectively. Species and sex identification matched expectations based on prior censuses of bat populations utilizing those roosts, though samples from some species were more or less common than anticipated within Fort Huachuca roosts. Nearly 62% of guano samples from Fort Drum contained DNA from Pseudogymnoascus destructans, where bats with wing damage from White‐nose Syndrome were commonly observed. Putative dietary items were detected in a majority of samples from insectivorous bats on Fort Drum (81%) and Fort Huachuca (63%). A minority of guano samples identified as the nectarivorous Leptonycteris yerbabuenae (28%) provided DNA sequences from putative forage plant species. Finally, DNA sequences from both putative ecto‐ and endoparasite taxa were detected in 35% and 56% of samples from Fort Drum and Fort Huachuca, respectively. This study demonstrates that the combination of noninvasive sampling, DNA metabarcoding, and sample and locus multiplexing provide a wide array of data that are otherwise difficult to obtain.

Vampyrum spectrum (Phyllostomidae) movement and prey revealed by radio-telemetry and DNA metabarcoding

The spectral bat (Vampyrum spectrum), the largest bat species in the Americas, is considered Near Threatened by the International Union for Conservation of Nature and is listed as a species of special concern or endangered in several countries throughout its range. Although the species is known as carnivorous, data on basic ecology, including habitat selection and primary diet items, are limited owing to its relative rarity and difficulty in capturing the species. Leveraging advances in DNA metabarcoding and using radio-telemetry, we present novel information on the diet and movement of V. spectrum based on locations of a radio-collared individual and fecal samples collected from its communal roost (three individuals) in the Lowland Dry Forest of southern Nicaragua. Using a non-invasive approach, we explored the diet of the species with genetic markers designed to capture a range of arthropods and vertebrate targets from fecal samples. We identified 27 species of vertebrate prey which included birds, rodents, and other bat species. Our evidence suggested that V. spectrum can forage on a variety of species, from those associated with mature forests to forest edge-dwellers. Characteristics of the roost and our telemetry data underscore the importance of large trees for roosting in mature forest patches for the species. These data can inform conservation efforts for preserving both the habitat and the prey items in remnants of mature forest required by Vampyrum spectrum to survive in landscape mosaics.