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

Gut microbial shifts in vampire bats linked to immunity due to changed diet in human disturbed landscapes

Anthropogenic land-use change alters wildlife habitats and modifies species composition, diversity, and contacts among wildlife, livestock, and humans. Such human-modified ecosystems have been associated with emerging infectious diseases, threatening human and animal health. However, human disturbance also creates new resources that some species can exploit. Common vampire bats (Desmodus rotundus) in Latin America constitute an important example, as their adaptation to human-modified habitats and livestock blood-feeding has implications for e.g., rabies transmission. Despite the well-known links between habitat degradation and disease emergence, few studies have explored how human-induced disturbance influences wildlife behavioural ecology and health, which can alter disease dynamics. To evaluate links among habitat disturbance, diet shifts, gut microbiota, and immunity, we quantified disturbance around roosting caves of common vampire bats in Costa Rica, measured their long-term diet preferences (livestock or wildlife blood) using stable isotopes of carbon and nitrogen, evaluated innate and adaptive immune markers, and characterized their gut microbiota. We observed that bats from roosting caves with more cattle farming nearby fed more on cattle blood. Moreover, gut microbial richness and the abundance of specific gut microbes differed according to feeding preferences. Interestingly, bats feeding primarily on wildlife blood harboured a higher abundance of the bacteria Edwardsiella sp., which tended to be associated with higher immunoglobulin G levels. Our results highlight how human land-use change may indirectly affect wildlife health and emerging infectious diseases through diet-induced shifts in microbiota, with implications for host immunity and potential consequences for susceptibility to pathogens.

High risk of bat bites in an indigenous village in Brazil: Warning of the re-emergence of rabies among the Maxakali People

Bat-mediated human rabies is a viral zoonotic disease that poses a serious threat to the public health of traditional peoples, especially indigenous populations that maintain primitive cultural and social habits, such as the Maxakali ethnic group, located in the southeastern region of Brazil. The sociocultural habit of this population led to the emergence between April and May 2022 of the viral spillover of rabies transmitted by bats, which decimated the lives of four children from this population who maintained contact with this animal as a recreational practice. Because the vampire bats Desmodus rotundus have exceptional ecology and social characteristics that can have important effects on the dynamics of viral dispersion in this indigenous population, I present the dynamics of contact between native children and the bat and the meaning of this relationship, which involves ritualistic and recreational significance. As important as knowing the reasons for this practice is discussing some intrinsic and extrinsic factors that imply risks that intensify the vulnerability of this population to the transmission of the rabies virus at any time. In view of this, I warn of the need to adopt efficient strategies to mitigate the risks of a new emergency in this region. Although emergency containment measures were carried out during the critical period of the outbreak, such animal and environmental control actions must become routine programmatic and structuring interventions. Essential for rabies surveillance in this population is to develop culturally adapted interethnic health education campaigns to guarantee the accessibility of the Maxakali indigenous people to the content taught, so that any attempt at domestication, captivity and recreational practices with bats of any species is discouraged, thus avoiding a possible re-emergence of this anthropozoonosis that has impacted not only the epidemiological scenario in this region, but throughout Brazil, and also throughout Latin America.

Qualitative identification of lonicerae japonicae flos in traditional chinese medicine using metabarcoding combined with specific mini-barcodes

BACKGROUND

Lonicerae japonicae flos, also known as Jinyinhua (JYH), is an important component of traditional Chinese patent medicine (TCPM) products. However, the potential for adulteration and substitution with low-quality materials highlights the need for a reliable and sensitive approach to identify the species composition of TCPM products for consumer safety.

METHODS AND RESULTS

We used universal ITS2 primers to amplify TCPMs containing JYH. However, the results were inconclusive, as only one operational taxonomic unit (OTU) was identified as Lonicera sp., which could not be identified at the species level. To confirm the species identification of Lonicera sp. in TCPM, we developed a short mini-barcode primer based on the psbA-trnH region, which, in combination with DNA metabarcoding technology, allowed for qualitative and quantitative analysis of artificially mixed samples. We applied the mini-barcode to distinguish TCPMs containing JYH and demonstrated its relatively accurate quantitative ability in identifying two Lonicera species.

CONCLUSIONS

Our study presents a method for qualitative and quantitative identification of JYH, providing a promising application of DNA metabarcoding technology in the quality control of TCPM products.

Metabarcoding of fecal DNA reveals the broad and flexible diet of a globally endangered bird

Knowing the diet of endangered wild animals is a prerequisite for species-specific conservation and habitat management. The Sichuan partridge Arborophila rufipectus is a globally endangered Galliformes species endemic to the mountains of southwest China. Existing information on the diet of this species is biased and fragmented owing to traditional observation methods. Little is known about their dietary composition or how they respond to temporal variations in food resources throughout the year. In this study, a dietary analysis was performed on 60 fecal samples using DNA Metabarcoding of invertebrates and plants to determine the primary animal and plant components of the diet across 3 critical periods of adult life history (breeding, postbreeding wandering, and overwintering). Preys from the dipteran order, followed by the lepidopteran and araneaen spp., were the predominant, animal-derived foods. Symplocos, Rubus, Celastrus, Holboellia, and Actinidia spp. supply a large abundance of fruits and seeds for this omnivorous bird. Substantial temporal dietary changes among the 3 periods and a general shift toward lower dietary diversity during the breeding season were observed, suggesting that the Sichuan partridge can adjust their diet according to the availability of food resources and their own needs. Characterizing the composition and seasonal changes in Sichuan partridge diets informs the habitat management of native flora (the plant taxa that can generate berries and seeds, such as Symplocos, Rubus, Celastrus, and Holboellia, which are likely of conservation interest) to achieve full life-cycle conservation.

SPECIES IN THE FECES: DNA METABARCODING TO DETECT POTENTIAL GASTROPOD HOSTS OF PARELAPHOSTRONGYLUS TENUIS Consumed By Moose (Alces Alces)

Our understanding of wildlife multihost pathogen transmission systems is often incomplete due to the difficulty of observing contact between hosts. Understanding these interactions can be critical for preventing disease-induced wildlife declines. The proliferation of high-throughput sequencing technologies provides new opportunities to better explore these cryptic interactions. Parelaphostrongylus tenuis, a multihost parasite, is a leading cause of death in some moose (Alces alces) populations threatened by local extinction in the midwestern and northeastern US and southeastern Canada. Moose contract P. tenuis by consuming infected gastropod intermediate hosts, but little is known about which gastropod species moose consume. To gain more insight, we used a genetic metabarcoding approach on 258 georeferenced and temporally stratified moose fecal samples collected May-October 2017-20 from a declining population in the north-central US. We detected moose consumption of three species of gastropods across five positive samples. Two of these (Punctum minutissimum and Helisoma sp.) have been minimally investigated for the ability to host P. tenuis, while one (Zonitoides arboreus) is a well-documented host. Moose consumption of gastropods documented herein occurred in June and September. Our findings prove that moose consume gastropod species known to become infected by P. tenuis and demonstrate that fecal metabarcoding can provide novel insight on interactions between hosts of a multispecies pathogen transmission system. After determining and improving test sensitivity, these methods may also be extended to document important interactions in other multihost disease systems.

Using eDNA to understand predator-prey interactions influenced by invasive species

Invasive predatory species may alter population dynamic processes of their prey and impact biological communities and ecosystem processes. Revealing biotic interactions, however, including the relationship between predator and prey, is a difficult task, in particular for species that are hard to monitor. Here, we present a case study that documents the utility of environmental DNA analysis (eDNA) to assess predator-prey interactions between two invasive fishes (Lepomis gibbosus, Pseudorasbora parva) and two potential amphibian prey species, (Triturus cristatus, Pelobates fuscus). We used species-specific TaqMan assays for quantitative assessment of eDNA concentrations from water samples collected from 89 sites across 31 ponds during three consecutive months from a local amphibian hotspot in Germany. We found a negative relationship between eDNA concentrations of the predators (fishes) and prey (amphibians) using Monte-Carlo tests. Our study highlights the potential of eDNA application to reveal predator-prey interactions and confirms the hypothesis that the observed local declines of amphibian species may be at least partly caused by recently introduced invasive fishes. Our findings have important consequences for local conservation management and highlight the usefulness of eDNA approaches to assess ecological interactions and guide targeted conservation action.

Haplotype-level metabarcoding of freshwater macroinvertebrate species: A prospective tool for population genetic analysis

Metabarcoding is a molecular-based tool capable of large quantity high-throughput species identification from bulk samples that is a faster and more cost-effective alternative to conventional DNA-sequencing approaches. Still, further exploration and assessment of the laboratory and bioinformatics strategies are required to unlock the potential of metabarcoding-based inference of haplotype information. In this study, we assessed the inference of freshwater macroinvertebrate haplotypes from metabarcoding data in a mock sample. We also examined the influence of DNA template concentration and PCR cycle on detecting true and spurious haplotypes. We tested this strategy on a mock sample containing twenty individuals from four species with known haplotypes based on the 658-bp Folmer region of the mitochondrial cytochrome c oxidase gene. We recovered fourteen zero-radius operational taxonomic units (zOTUs) of 421-bp length, with twelve zOTUs having a 100% match with the Sanger haplotype sequences. High-quality reads relatively increased with increasing PCR cycles, and the relative abundance of each zOTU was consistent for each cycle. This suggests that increasing the PCR cycles from 24 to 64 did not affect the relative abundance of each zOTU. As metabarcoding becomes more established and laboratory protocols and bioinformatic pipelines are continuously being developed, our study demonstrated the method's ability to infer intraspecific variability while highlighting the challenges that must be addressed before its eventual application for population genetic studies.

A review of the diet of the common vampire bat (Desmodus rotundus) in the context of anthropogenic change

The common vampire bat (Desmodus rotundus) maintains a diverse, sanguivorous diet, utilizing a broad range of prey taxa. As anthropogenic change alters the distribution of this species, shifts in predator-prey interactions are expected. Understanding prey richness and patterns of prey selection is, thus, increasingly informative from ecological, epidemiological, and economic perspectives. We reviewed D. rotundus diet and assessed the geographical, taxonomical, and behavioral features to find 63 vertebrate species within 21 orders and 45 families constitute prey, including suitable host species in regions of invasion outside D. rotundus' range. Rodentia contained the largest number of species utilized by D. rotundus, though cattle were the most commonly reported prey source, likely linked to the high availability of livestock and visibility of bite wounds compared to wildlife. Additionally, there was tendency to predate upon species with diurnal activity and social behavior, potentially facilitating convenient and nocturnal predation. Our review highlights the dietary heterogeneity of D. rotundus across its distribution. We define D. rotundus as a generalist predator, or parasite, depending on the ecological definition of its symbiont roles in an ecosystem (i.e., lethal vs. non-lethal blood consumption). In view of the eminent role of D. rotundus in rabies virus transmission and its range expansion, an understanding of its ecology would benefit public health, wildlife management, and agriculture.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42991-023-00358-3.

Rabies in a Captive Lowland Tapir (Tapirus terrestris)

Rabies is a severe viral zoonosis of mammals and causes irreversible neurological damage. We describe the clinical presentation and anatomopathological lesions of rabies in a captive lowland tapir (Tapirus terrestris) in Bauru, São Paulo State, Brazil. The clinical course of the disease lasted 6 days and was characterized by progressive neurological deterioration and death. The main anatomopathological findings were non-suppurative encephalitis and presence of Negri bodies within neurons. Direct immunofluorescence and mouse inoculation tests were positive for rabies virus. This is the first report of rabies in a lowland tapir and highlights the importance of disease prevention under managed care and continuous control measures in urbanized environments.

Ecological determinants of rabies virus dynamics in vampire bats and spillover to livestock

The pathogen transmission dynamics in bat reservoirs underpin efforts to reduce risks to human health and enhance bat conservation, but are notoriously challenging to resolve. For vampire bat rabies, the geographical scale of enzootic cycles, whether environmental factors modulate baseline risk, and how within-host processes affect population-level dynamics remain unresolved. We studied patterns of rabies exposure using an 11-year, spatially replicated sero-survey of 3709 Peruvian vampire bats and co-occurring outbreaks in livestock. Seroprevalence was correlated among nearby sites but fluctuated asynchronously at larger distances. A generalized additive mixed model confirmed spatially compartmentalized transmission cycles, but no effects of bat demography or environmental context on seroprevalence. Among 427 recaptured bats, we observed long-term survival following rabies exposure and antibody waning, supporting hypotheses that immunological mechanisms influence viral maintenance. Finally, seroprevalence in bats was only weakly correlated with outbreaks in livestock, reinforcing the challenge of spillover prediction even with extensive data. Together our results suggest that rabies maintenance requires transmission among multiple, nearby bat colonies which may be facilitated by waning of protective immunity. However, the likelihood of incursions and dynamics of transmission within bat colonies appear largely independent of bat ecology. The implications of these results for spillover anticipation and controlling transmission at the source are discussed.

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.

Strategies for sample labelling and library preparation in DNA metabarcoding studies

Metabarcoding of DNA extracted from environmental or bulk specimen samples is increasingly used to profile biota in basic and applied biodiversity research because of its targeted nature that allows sequencing of genetic markers from many samples in parallel. To achieve this, PCR amplification is carried out with primers designed to target a taxonomically informative marker within a taxonomic group, and sample-specific nucleotide identifiers are added to the amplicons prior to sequencing. The latter enables assignment of the sequences back to the samples they originated from. Nucleotide identifiers can be added during the metabarcoding PCR and during "library preparation", that is, when amplicons are prepared for sequencing. Different strategies to achieve this labelling exist. All have advantages, challenges and limitations, some of which can lead to misleading results, and in the worst case compromise the fidelity of the metabarcoding data. Given the range of questions addressed using metabarcoding, ensuring that data generation is robust and fit for the chosen purpose is critically important for practitioners seeking to employ metabarcoding for biodiversity assessments. Here, we present an overview of the three main workflows for sample-specific labelling and library preparation in metabarcoding studies on Illumina sequencing platforms; one-step PCR, two-step PCR, and tagged PCR. Further, we distill the key considerations for researchers seeking to select an appropriate metabarcoding strategy for their specific study. Ultimately, by gaining insights into the consequences of different metabarcoding workflows, we hope to further consolidate the power of metabarcoding as a tool to assess biodiversity across a range of applications.

A two-step metagenomics approach for the identification and mitochondrial DNA contig assembly of vertebrate prey from the blood meals of common vampire bats (Desmodus rotundus)

The feeding behaviour of the sanguivorous common vampire bat (Desmodus rotundus) facilitates the transmission of pathogens that can impact both human and animal health. To formulate effective strategies in controlling the spread of diseases, there is a need to obtain information on which animals they feed on. One DNA-based approach, shotgun sequencing, can be used to obtain such information. Even though it is costly, shotgun sequencing can be used to simultaneously retrieve prey and vampire bat mitochondrial DNA for population studies within one round of sequencing. However, due to the challenges of analysing shotgun sequenced metagenomic data such as false negatives/positives and typically low proportion of reads mapped to diet items, shotgun sequencing has not been used for the identification of prey from common vampire bat blood meals. To overcome these challenges and generate longer mitochondrial contigs which could be useful for prey population studies, we shotgun sequenced common vampire bat blood meal samples (n = 8) and utilised a two-step metagenomic approach based on combining existing bioinformatic workflows (alignment and mtDNA contig assembly) to identify prey. After validating our results from detections made through metabarcoding, we accurately identified the common vampire bats’ prey in six out of eight samples without any false positives. We also generated prey mitochondrial contig lengths between 138 bp to 3231 bp (median = 770 bp, Q1 = 262 bp, Q3 = 1766 bp). This opens the potential to conduct phylogenetic and phylogeographic monitoring of elusive prey species in future studies, through the analyses of blood meal metagenomic data.

Seroprevalence to Rabies Virus in Wildlife in Brazil

Serum samples of 638 free-ranging wild mammals from São Paulo state, Brazil, were tested for neutralizing antibodies against rabies virus by the rapid fluorescent focus inhibition test. Overall seroprevalence was 1.7% among 24 species surveyed, with individuals of six species having positive results indicating exposure to rabies virus.

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

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

Harnessing the Power of Metabarcoding in the Ecological Interpretation of Plant-Pollinator DNA Data: Strategies and Consequences of Filtering Approaches

Although DNA metabarcoding of pollen mixtures has been increasingly used in the field of pollination biology, methodological and interpretation issues arise due to its high sensitivity. Filtering or maintaining false positives, contaminants, and rare taxa or molecular features could lead to different ecological results. Here, we reviewed how this choice has been addressed in 43 studies featuring pollen DNA metabarcoding, which highlighted a very high heterogeneity of filtering methods. We assessed how these strategies shaped pollen assemblage composition, species richness, and interaction networks. To do so, we compared four processing methods: unfiltering, filtering with a proportional 1% of sample reads, a fixed threshold of 100 reads, and the ROC approach (Receiver Operator Characteristic). The results indicated that filtering impacted species composition and reduced species richness, with ROC emerging as a conservative approach. Moreover, in contrast to unfiltered networks, filtering decreased network Connectance and Entropy, and it increased Modularity and Connectivity, indicating that using cut-off thresholds better describes interactions. Overall, unfiltering might compromise reliable ecological interpretations, unless a study targets rare species. We discuss the suitability of each filtering type, plead for justifying filtering strategies on biological or methodological bases and for developing shared approaches to make future studies more comparable.

Dung beetles as samplers of mammals in Malaysian Borneo-a test of high throughput metabarcoding of iDNA

Invertebrate-derived DNA (iDNA) sampling in biodiversity surveys is becoming increasingly widespread, with most terrestrial studies relying on DNA derived from the gut contents of blood-feeding invertebrates, such as leeches and mosquitoes. Dung beetles (superfamily Scarabaeoidea) primarily feed on the faecal matter of terrestrial vertebrates and offer several potential benefits over blood-feeding invertebrates as samplers of vertebrate DNA. Importantly, these beetles can be easily captured in large numbers using simple, inexpensive baited traps, are globally distributed, and occur in a wide range of habitats. To build on the few existing studies demonstrating the potential of dung beetles as sources of mammalian DNA, we subjected the large-bodied, Bornean dung beetle (Catharsius renaudpauliani) to a controlled feeding experiment. We analysed DNA from gut contents at different times after feeding using qPCR techniques. Here, we first describe the window of DNA persistence within a dung beetle digestive tract. We found that the ability to successfully amplify cattle DNA decayed over relatively short time periods, with DNA copy number decreasing by two orders of magnitude in just 6 h. In addition, we sampled communities of dung beetles from a lowland tropical rainforest in Sabah, Malaysia, in order to test whether it is possible to identify vertebrate sequences from dung beetle iDNA. We sequenced both the gut contents from large dung beetle species, as well as whole communities of smaller beetles. We successfully identified six mammalian species from our samples, including the bearded pig (Sus barbatus) and the sambar deer (Rusa unicolor)—both vulnerable species on the IUCN red list. Our results represent the first use of dung beetle iDNA to sample Southeast Asian vertebrate fauna, and highlight the potential for dung beetle iDNA to be used in future biodiversity monitoring surveys.

The potential of aquatic bloodfeeding and nonbloodfeeding leeches as a tool for iDNA characterisation

Leeches play important roles in food webs due to their abundance, diversity and feeding habits. Studies using invertebrate‐derived DNA (iDNA) extracted from leech gut contents to target vertebrate DNA have focused on the Indo‐Pacific region and mainly leveraged the leech family Haemadipsidae, composed of bloodfeeding terrestrial leeches, while predatory, fluid/tissue‐feeding and aquatic bloodfeeding species have been largely disregarded. While there is some general knowledge regarding the taxonomic groups that leeches prefer to feed on, detailed taxonomic resolution is missing and, therefore, their potential use for monitoring animals is unknown. In this study, 116 leeches from 12 species (six families) and spanning the three feeding habits were collected in Mexico and Canada. We used DNA metabarcoding to investigate their diet and assess their potential use for biodiversity monitoring. We detected vertebrates from five orders including fish, turtles and birds in the diet of aquatic bloodfeeding leeches; eight invertebrate orders of annelids, arthropods and molluscs in leeches that feed on body fluids and tissues; and 10 orders of invertebrates belonging to Arthropoda and Annelida, as well as one vertebrate and one parasitic nematode, in predatory leeches. These results show the potential use of iDNA from aquatic bloodfeeding leeches for retrieving vertebrate taxa, and from predatory and fluid‐feeding leeches for invertebrates. Our study provides information about the dietary range of freshwater leeches and one terrestrial leech and contributes proof‐of‐concept for the use of these leeches for animal monitoring, expanding our knowledge of the use of iDNA from leech gut contents to North America.

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.

The Future of DNA Barcoding: Reflections from Early Career Researchers

Over the last two decades, the use of DNA barcodes has transformed our ability to identify and assess life on our planet. Both strengths and weaknesses of the method have been exemplified through thousands of peer-reviewed scientific articles. Given the novel sequencing approaches, currently capable of generating millions of reads at low cost, we reflect on the questions: What will the future bring for DNA barcoding? Will identification of species using short, standardized fragments of DNA stand the test of time? We present reflected opinions of early career biodiversity researchers in the form of a SWOT analysis and discuss answers to these questions.

Molecular diet analysis of neotropical bats based on fecal DNA metabarcoding

Bat communities in the Neotropics are some of the most speciose assemblages of mammals on Earth, with regions supporting more than 100 sympatric species with diverse feeding ecologies. Because bats are small, nocturnal, and volant, it is difficult to directly observe their feeding habits, which has resulted in their classification into broadly defined dietary guilds (e.g., insectivores, carnivores, and frugivores). Apart from these broad guilds, we lack detailed dietary information for many species and therefore have only a limited understanding of interaction networks linking bats and their diet items. In this study, we used DNA metabarcoding of plants, arthropods, and vertebrates to investigate the diets of 25 bat species from the tropical dry forests of Lamanai, Belize. Our results report some of the first detection of diet items for the focal bat taxa, adding rich and novel natural history information to the field of bat ecology. This study represents a comprehensive first effort to apply DNA metabarcoding to bat diets at Lamanai and provides a useful methodological framework for future studies testing hypotheses about coexistence and niche differentiation in the context of modern high-throughput molecular data.

DNA metabarcoding provides insights into seasonal diet variations in Chinese mole shrew (Anourosorex squamipes) with potential implications for evaluating crop impacts

Diet analysis of potential small mammals pest species is important for understanding feeding ecology and evaluating their impact on crops and stored foods. Chinese mole shrew (Anourosorex squamipes), distributed in Southwest China, has previously been reported as a farmland pest. Effective population management of this species requires a better understanding of its diet, which can be difficult to determine with high taxonomic resolution using conventional microhistological methods. In this study, we used two DNA metabarcoding assays to identify 38 animal species and 65 plant genera from shrew stomach contents, which suggest that A. squamipes is an omnivorous generalist. Earthworms are the most prevalent (>90%) and abundant (>80%) food items in the diverse diet of A. squamipes. Species of the Fabaceae (frequency of occurrence [FO]: 88%; such as peanuts) and Poaceae (FO: 71%; such as rice) families were the most common plant foods identified in the diet of A. squamipes. Additionally, we found a seasonal decrease in the diversity and abundance of invertebrate foods from spring and summer to winter. Chinese mole shrew has a diverse and flexible diet throughout the year to adapt to seasonal variations in food availability, contributing to its survival even when food resources are limited. This study provides a higher resolution identification of the diet of A. squamipes than has been previously described and is valuable for understanding shrew feeding ecology as well as evaluating possible species impacts on crops.

DNA-Based Arthropod Diversity Assessment in Amazonian Iron Mine Lands Show Ecological Succession Towards Undisturbed Reference Sites

Human activities change natural landscapes, and in doing so endanger biodiversity and associated ecosystem services. To reduce the net impacts of these activities, such as mining, disturbed areas are rehabilitated and restored. During this process, monitoring is important to ensure that desired trajectories are maintained. In the Carajás region of the Brazilian Amazon, exploration for iron ores has transformed the original ecosystem; natural forest and a savanna formation with lateritic iron duricrust outcrops named canga. Here, native vegetation is logged and topsoil removed and deposited in waste piles along with mine waste. During rehabilitation, these waste piles are hydroseeded with non-native plant species to achieve rapid revegetation. Further, seeds of native canga and forest plant species are planted to point ecological succession towards natural ecosystems. In this study, we investigate diversity and composition of the arthropod community along a post-mining rehabilitation and restoration gradient, taking seasonality and primer bias into account. We use DNA metabarcoding of bulk arthropod samples collected in both the dry and rainy seasons from waste-pile benches at various stages of revegetation: non-revegetated exposed soils, initial stage with one-to-three-year-old stands, intermediate stage with four-to-five-year-old stands, and advanced stage with six-to-seven-year-old stands. We use samples from undisturbed cangas and forests as reference sites. In addition, we vegetation diversity and structure were measured to investigate relations between arthropod community and vegetation structure. Our results show that, over time, the arthropod community composition of the waste piles becomes more similar to the reference forests, but not to the reference cangas. Nevertheless, even the communities in the advanced-stage waste piles are different from the reference forests, and full restoration in these highly diverse ecosystems is not achieved, even after 6 to 7 years. Finally, our results show seasonal variation in arthropod communities and primer bias.

Sensory ecology of the frog-eating bat, Trachops cirrhosus, from DNA metabarcoding and behavior

Metabarcoding of prey DNA from fecal samples can be used to design behavioral experiments to study the foraging behavior and sensory ecology of predators. The frog-eating bat, Trachops cirrhosus, eavesdrops on the mating calls of its anuran prey. We captured wild T. cirrhosus and identified prey remains in the bats’ fecal samples using DNA metabarcoding of two gene regions (CO1 and 16S). Bats were preying on frogs previously unknown in their diet, such as species in the genus Pristimantis, which occurred in 29% of T. cirrhosus samples. Twenty-three percent of samples also contained DNA of Anolis lizards. We additionally report apparently rare predation events on hummingbirds and heterospecific bats. We used results from metabarcoding to design acoustic and 3D model stimuli to present to bats in behavioral experiments. We show predatory responses by T. cirrhosus to the calls of the frog Pristimantis taeniatus and to the rustling sounds of anoles moving through leaf-litter, as well as attacks on a stuffed hummingbird and a plastic anole model. The combination of species-specific dietary information from metabarcoding analyses with behavioral responses to prey cues provides a unique window into the foraging ecology of predators that are difficult to observe in the wild.

Defining New Pathways to Manage the Ongoing Emergence of Bat Rabies in Latin America

Rabies transmitted by common vampire bats (Desmodus rotundus) has been known since the early 1900s but continues to expand geographically and in the range of species and environments affected. In this review, we present current knowledge of the epidemiology and management of rabies in D. rotundus and argue that it can be reasonably considered an emerging public health threat. We identify knowledge gaps related to the landscape determinants of the bat reservoir, reduction in bites on humans and livestock, and social barriers to prevention. We discuss how new technologies including autonomously-spreading vaccines and reproductive suppressants targeting bats might manage both rabies and undesirable growth of D. rotundus populations. Finally, we highlight widespread under-reporting of human and animal mortality and the scarcity of studies that quantify the efficacy of control measures such as bat culling. Collaborations between researchers and managers will be crucial to implement the next generation of rabies management in Latin America.

A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

Tagsteady: A metabarcoding library preparation protocol to avoid false assignment of sequences to samples

Metabarcoding of environmental DNA (eDNA) and DNA extracted from bulk specimen samples is a powerful tool in studies of biodiversity, diet and ecological interactions as its inherent labelling of amplicons allows sequencing of taxonomically informative genetic markers from many samples in parallel. However, the occurrence of so-called 'tag-jumps' can cause incorrect assignment of sequences to samples and artificially inflate diversity. Two steps during library preparation of pools of 5' nucleotide-tagged amplicons have been suggested to cause tag-jumps: (a) T4 DNA polymerase blunt-ending in the end-repair step and (b) postligation PCR amplification of amplicon libraries. The discovery of tag-jumps has led to recommendations to only carry out metabarcoding PCR amplifications with primers carrying twin-tags to ensure that tag-jumps cannot result in false assignments of sequences to samples. As this increases both cost and workload, a metabarcoding library preparation protocol which circumvents the two steps that causes tag-jumps is needed. Here, we demonstrate Tagsteady, a PCR-free metabarcoding Illumina library preparation protocol for pools of nucleotide-tagged amplicons that enables efficient and cost-effective generation of metabarcoding data with virtually no tag-jumps. We use pools of twin-tagged amplicons to investigate the effect of T4 DNA polymerase blunt-ending and postligation PCR on the occurrence of tag-jumps and demonstrate that both blunt-ending and postligation PCR, alone or together, can result in detrimental amounts of tag-jumps (here, up to ca. 49% of total sequences), while leaving both steps out (the Tagsteady protocol) results in amounts of sequences carrying new combinations of used tags (tag-jumps) comparable to background contamination.

In silico and empirical evaluation of twelve metabarcoding primer sets for insectivorous diet analyses

During the most recent decade, environmental DNA metabarcoding approaches have been both developed and improved to minimize the biological and technical biases in these protocols. However, challenges remain, notably those relating to primer design. In the current study, we comprehensively assessed the performance of ten COI and two 16S primer pairs for eDNA metabarcoding, including novel and previously published primers. We used a combined approach of in silico, in vivo-mock community (33 arthropod taxa from 16 orders), and guano-based analyses to identify primer sets that would maximize arthropod detection and taxonomic identification, successfully identify the predator (bat) species, and minimize the time and financial costs of the experiment. We focused on two insectivorous bat species that live together in mixed colonies: the greater horseshoe bat (Rhinolophus ferrumequinum) and Geoffroy's bat (Myotis emarginatus). We found that primer degeneracy is the main factor that influences arthropod detection in silico and mock community analyses, while amplicon length is critical for the detection of arthropods from degraded DNA samples. Our guano-based results highlight the importance of detecting and identifying both predator and prey, as guano samples can be contaminated by other insectivorous species. Moreover, we demonstrate that amplifying bat DNA does not reduce the primers' capacity to detect arthropods. We therefore recommend the simultaneous identification of predator and prey. Finally, our results suggest that up to one-third of prey occurrences may be unreliable and are probably not of primary interest in diet studies, which may decrease the relevance of combining several primer sets instead of using a single efficient one. In conclusion, this study provides a pragmatic framework for eDNA primer selection with respect to scientific and methodological constraints.

Abortive vampire bat rabies infections in Peruvian peridomestic livestock

Rabies virus infections normally cause universally lethal encephalitis across mammals. However, 'abortive infections' which are resolved prior to the onset of lethal disease have been described in bats and a variety of non-reservoir species. Here, we surveyed rabies virus neutralizing antibody titers in 332 unvaccinated livestock of 5 species from a vampire bat rabies endemic region of southern Peru where livestock are the main food source for bats. We detected rabies virus neutralizing antibody titers in 11, 5 and 3.6% of cows, goats and sheep respectively and seropositive animals did not die from rabies within two years after sampling. Seroprevalence was correlated with the number of local livestock rabies mortalities reported one year prior but also one year after sample collection. This suggests that serological status of livestock can indicate the past and future levels of rabies risk to non-reservoir hosts. To our knowledge, this is the first report of anti-rabies antibodies among goats and sheep, suggesting widespread abortive infections among livestock in vampire bat rabies endemic areas. Future research should resolve the within-host biology underlying clearance of rabies infections. Cost-effectiveness analyses are also needed to evaluate whether serological monitoring of livestock can be a viable complement to current monitoring of vampire bat rabies risk based on animal mortalities alone.

Synergistic Effects of Grassland Fragmentation and Temperature on Bovine Rabies Emergence

In 2007, common vampire bats were the source of the first outbreak of paralytic bovine rabies in Uruguay. The outbreak coincided in space and time with the fragmentation of native grasslands for monospecific forestry for wood and cellulose production. Using spatial analyses, we show that the increase in grassland fragmentation, together with the minimum temperature in the winter, accounts for the spatial pattern of outbreaks in the country. We propose that fragmentation may increase the connectivity of vampire bat colonies by promoting the sharing of feeding areas, while temperature modulates their home range plasticity. While a recent introduction of the virus from neighboring Brazil could have had an effect on outbreak occurrence, we show here that the distribution of rabies cases is unlikely to be explained by only an invasion process from Brazil. In accordance with previous modeling efforts, an increase in connectivity may promote spatial persistence of rabies virus within vampire bat populations. Our results suggest that land use planning might help to reduce grassland fragmentation and thus reduce risk of rabies transmission to livestock. This will be especially important in the context of climatic changes and increasing minimum temperatures in the winter.

An appetite for pests: Synanthropic insectivorous bats exploit cotton pest irruptions and consume various deleterious arthropods

Conservation biological control (CBC) seeks to minimize the deleterious effects of agricultural pests by enhancing the efficiency of natural enemies. Despite the documented potential of insectivorous bats to consume pests, many synanthropic bat species are still underappreciated as beneficial species. We investigated the diet of Kuhl's pipistrelle (Pipistrellus kuhlii), a common synanthropic insectivorous bat that forages in urban and agricultural areas, to determine whether it may function as a natural enemy in CBC. Faecal samples of P. kuhlii were collected throughout the cotton-growing season from five roost sites near cotton fields located in a Mediterranean agroecosystem, Israel, and analyzed using DNA metabarcoding. Additionally, data on estimated abundance of major cotton pests were collected. We found that the diet of P. kuhlii significantly varied according to sites and dates and comprised 27 species of agricultural pests that were found in 77.2% of the samples, including pests of key economic concern. The dominant prey was the widespread cotton pest, the pink bollworm, Pectinophora gossypiella, found in 31% of the samples and in all the roosts. Pink bollworm abundance was positively correlated with its occurrence in the bat diet. Furthermore, the bats' dietary breadth narrowed, while temporal dietary overlap increased, in relation to increasing frequencies of pink bollworms in the diet. This suggests that P. kuhlii exploits pink bollworm irruptions by opportunistic feeding. We suggest that synanthropic bats provide important pest suppression services, may function as CBC agents of cotton pests and potentially contribute to suppress additional deleterious arthropods found in their diet in high frequencies.

Relations between topography, feeding sites, and foraging behavior of the vampire bat, Desmodus rotundus

The vampire bat (Desmodus rotundus) plays a crucial role in the maintenance and transmission of the rabies virus to humans and livestock, impacting public health and economic production. Its importance lies not only in its capacity to transmit the virus but also in its ability to adapt to anthropic changes, as expressed in its wide geographic distribution in Latin America. Deforestation, livestock intensification, and other human activities have reduced the abundance of its wild prey but have also provided new and abundant shelter and foraging resources for the vampire bats. We used radiotelemetry to evaluate relations between topography, feeding site choice, and foraging behavior in southeastern Brazil, where three occupied D. rotundus roosts, out of 11 possible, were systematically monitored throughout a 1-year period once every 2 months. Sixty-two vampire bats were captured; biometric data were collected and 44 radiotransmitters were installed in adult individuals, producing telemetry data that were recorded in VHF receptors installed in the farms. Elevation of the roosts was related to the farms attacked by the vampire bats. Understanding the use of the environment and resources by vampire bats is critical to improving rabies control aiming at the reduction of disease impacts. From the perspective of the official veterinary service, telemetry would be ineffective as a rabies surveillance system due to the costs and limitations of the available technology. However, livestock rabies control measures would be greatly improved if ecological characteristics of the vampire bat were considered.

Demographic and environmental drivers of metagenomic viral diversity in vampire bats

Viruses infect all forms of life and play critical roles as agents of disease, drivers of biochemical cycles and sources of genetic diversity for their hosts. Our understanding of viral diversity derives primarily from comparisons among host species, precluding insight into how intraspecific variation in host ecology affects viral communities or how predictable viral communities are across populations. Here we test spatial, demographic and environmental hypotheses explaining viral richness and community composition across populations of common vampire bats, which occur in diverse habitats of North, Central and South America. We demonstrate marked variation in viral communities that was not consistently predicted by a null model of declining community similarity with increasing spatial or genetic distances separating populations. We also find no evidence that larger bat colonies host greater viral diversity. Instead, viral diversity follows an elevational gradient, is enriched by juvenile-biased age structure, and declines with local anthropogenic food resources as measured by livestock density. Our results establish the value of linking the modern influx of metagenomic sequence data with comparative ecology, reveal that snapshot views of viral diversity are unlikely to be representative at the species level, and affirm existing ecological theories that link host ecology not only to single pathogen dynamics but also to viral communities.

Validation of COI metabarcoding primers for terrestrial arthropods

Metabarcoding can rapidly determine the species composition of bulk samples and thus aids biodiversity and ecosystem assessment. However, it is essential to use primer sets that minimize amplification bias among taxa to maximize species recovery. Despite this fact, the performance of primer sets employed for metabarcoding terrestrial arthropods has not been sufficiently evaluated. This study tests the performance of 36 primer sets on a mock community containing 374 insect species. Amplification success was assessed with gradient PCRs and the 21 most promising primer sets selected for metabarcoding. These 21 primer sets were also tested by metabarcoding a Malaise trap sample. We identified eight primer sets, mainly those including inosine and/or high degeneracy, that recovered more than 95% of the species in the mock community. Results from the Malaise trap sample were congruent with the mock community, but primer sets generating short amplicons produced potential false positives. Taxon recovery from both mock community and Malaise trap sample metabarcoding were used to select four primer sets for additional evaluation at different annealing temperatures (40-60 °C) using the mock community. The effect of temperature varied by primer pair but overall it only had a minor effect on taxon recovery. This study reveals the weak performance of some primer sets employed in past studies. It also demonstrates that certain primer sets can recover most taxa in a diverse species assemblage. Thus, based our experimental set up, there is no need to employ several primer sets targeting the same gene region. We identify several suitable primer sets for arthropod metabarcoding, and specifically recommend BF3 + BR2, as it is not affected by primer slippage and provides maximal taxonomic resolution. The fwhF2 + fwhR2n primer set amplifies a shorter fragment and is therefore ideal when targeting degraded DNA (e.g., from gut contents).

Assessment of seasonal variation of diet composition in rodents using DNA barcoding and Real-Time PCR

The study of animal diet and feeding behaviour is a fundamental tool for the illustration of the ecological role of species in the ecosystem. However, size and quality of food intake samples make it hard for researchers to describe the diet composition of many small species. In our study, we exploited genomic tools for the analysis of the diet composition of the Savi’s pine vole (Microtus savii) using DNA barcoding and qPCR techniques for the identification of ingested plant species retrieved from stomach contents. In contrast with previous studies, we found that, despite being a fossorial species, the Savi’s pine vole is a selective feeder that undergoes intense superficial activity in search for food. In addition, our study shows that with a a priori knowledge of the candidate plant species included in animal diet, qPCR is a powerful tool to assess presence/absence, frequency of occurrence and electivity of ingested species. We conclude that this approach offers new opportunities to implement the analysis of food selection in small animals, thereby revealing a detailed picture of plant-animal interactions.

The menu varies with metabarcoding practices: A case study with the bat Plecotus auritus

Metabarcoding of feces has revolutionized the knowledge of animal diets by providing unprecedented resolution of consumed resources. However, it is still unclear how different methodological approaches influence the ecological conclusions that can be drawn from such data. Here, we propose a critical evaluation of several data treatments on the inferred diet of the bat Plecotus auritus using guano regularly collected from various colonies throughout the entire active season. First and unlike previous claims, our data indicates that DNA extracted from large amounts of fecal material issued from guano accumulates yield broader taxonomic diversity of prey than smaller numbers of pellets would do, provided that extraction buffer volumes are adapted to such increased amounts of material. Second, trophic niche analyses based on prey occurrence data uncover strong seasonality in the bat’s diet and major differences among neighboring maternity colonies. Third, while the removal of rare prey items is not always warranted as it introduces biases affecting particularly samples with greater prey species richness. Fourth, examination of distinct taxonomic depths in diet analyses highlights different aspects of food consumption providing a better understanding of the consumer’s diet. Finally, the biologically meaningful patterns recovered with presence-absence approaches are virtually lost when attempting to quantify prey consumed using relative read abundances. Even in an ideal situation where reference barcodes are available for most potential prey species, inferring realistic patterns of prey consumption remains relatively challenging. Although best practice in metabarcoding analyses will depend on the aims of the study, several previous methodological recommendations seem unwarranted for studying such diverse diets as that of brown long-eared bats.

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.