Using metabarcoding to compare the suitability of two blood-feeding leech species for sampling mammalian diversity in North Borneo

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.

Measuring protected-area effectiveness using vertebrate distributions from leech iDNA

Protected areas are key to meeting biodiversity conservation goals, but direct measures of effectiveness have proven difficult to obtain. We address this challenge by using environmental DNA from leech-ingested bloodmeals to estimate spatially-resolved vertebrate occupancies across the 677 km² Ailaoshan reserve in Yunnan, China. From 30,468 leeches collected by 163 park rangers across 172 patrol areas, we identify 86 vertebrate species, including amphibians, mammals, birds and squamates. Multi-species occupancy modelling shows that species richness increases with elevation and distance to reserve edge. Most large mammals (e.g. sambar, black bear, serow, tufted deer) follow this pattern; the exceptions are the three domestic mammal species (cows, sheep, goats) and muntjak deer, which are more common at lower elevations. Vertebrate occupancies are a direct measure of conservation outcomes that can help guide protected-area management and improve the contributions that protected areas make towards global biodiversity goals. Here, we show the feasibility of using invertebrate-derived DNA to estimate spatially-resolved vertebrate occupancies across entire protected areas.

Invertebrate-derived eDNA (iDNA) is an emerging tool for taxonomic and spatial biodiversity monitoring. Here, the authors use metabarcoding of leech-derived iDNA to estimate vertebrate occupancy over an entire protected area, the Ailaoshan Nature Reserve, China.

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.

Characterization of the complete mitochondrial genome of Haemadipsa tianmushana Song 1977 (Hirudiniformes, Haemadipsidae) and its phylogenetic analysis

The complete mitochondrial genome of Haemadipsa tianmushana Song 1977 from China has been determined and reported for the first time in this study. It was 14,625 bp in length and consisted of 22 tRNA genes, 2 rRNA genes, 13 protein-coding genes (PCGs), and 3 control regions. The nucleotide base content of the complete mitogenome for this species was 35.1% A, 10.5% C, 11.6% G, and 42.8% T. The tRNA genes were ranged from 57 bp (Ser^TCT) to 66 bp (Gln^TTG) in length. The phylogenetic analyses indicated that Hirudinea is a mono-phyletic clade. And it includes Whitmania acranulata, Whitmania pigra, Whitmania laevis, Zeylanicobdella arugamensis, Ozobranchus jantseanus and Placobdella lamothei. In Hirudiniformes, H. tianmushana and three species of Haemopidae were obviously clustered into two independent branches. This result is consistent with a taxonomy that they all belong to the same suborder. This study adds to the genetic resources currently available for the species.

Measuring biodiversity from DNA in the air

The crisis of declining biodiversity¹ exceeds our current ability to monitor changes in ecosystems. Rapid terrestrial biomonitoring approaches are essential to quantify the causes and consequences of global change. Environmental DNA² has revolutionized aquatic ecology,³ permitting population monitoring⁴ and remote diversity assessments matching or outperforming conventional methods of community sampling.^3-5 Despite this model, similar methods have not been widely adopted in terrestrial ecosystems. Here, we demonstrate that DNA from terrestrial animals can be filtered, amplified, and then sequenced from air samples collected in natural settings representing a powerful tool for terrestrial ecology. We collected air samples at a zoological park, where spatially confined non-native species allowed us to track DNA sources. We show that DNA can be collected from air and used to identify species and their ecological interactions. Air samples contained DNA from 25 species of mammals and birds, including 17 known terrestrial resident zoo species. We also identified food items from air sampled in enclosures and detected taxa native to the local area, including the Eurasian hedgehog, endangered in the United Kingdom. Our data demonstrate that airborne eDNA concentrates around recently inhabited areas but disperses away from sources, suggesting an ecology to airborne eDNA and the potential for sampling at a distance. Our findings demonstrate the profound potential of air as a source of DNA for global terrestrial biomonitoring.

Invertebrates for vertebrate biodiversity monitoring: Comparisons using three insect taxa as iDNA samplers

Metabarcoding of environmental DNA (eDNA) is now widely used to build diversity profiles from DNA that has been shed by species into the environment. There is substantial interest in the expansion of eDNA approaches for improved detection of terrestrial vertebrates using invertebrate-derived DNA (iDNA) in which hematophagous, sarcophagous, and coprophagous invertebrates sample vertebrate blood, carrion, or faeces. Here, we used metabarcoding and multiple iDNA samplers (carrion flies, sandflies, and mosquitos) collected from 39 forested sites in the southern Amazon to profile gamma and alpha diversity. Our main objectives were to (1) compare diversity found with iDNA to camera trapping, which is the conventional method of vertebrate diversity surveillance; and (2) compare each of the iDNA samplers to assess the effectiveness, efficiency, and potential biases associated with each sampler. In total, we collected and analysed 1759 carrion flies, 48,686 sandflies, and 4776 mosquitos. Carrion flies revealed the greatest total vertebrate species richness at the landscape level, despite the least amount of sampling effort and the fewest number of individuals captured for metabarcoding, followed by sandflies. Camera traps had the highest median species richness at the site-level but showed strong bias towards carnivore and ungulate species and missed much of the diversity described by iDNA methods. Mosquitos showed a strong feeding preference for humans as did sandflies for armadillos, thus presenting potential utility to further study related to host-vector interactions.

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.

eDNAir: proof of concept that animal DNA can be collected from air sampling

Environmental DNA (eDNA) is one of the fastest developing tools for species biomonitoring and ecological research. However, despite substantial interest from research, commercial and regulatory sectors, it has remained primarily a tool for aquatic systems with a small amount of work in substances such as soil, snow and rain. Here we demonstrate that eDNA can be collected from air and used to identify mammals. Our proof of concept successfully demonstrated that eDNA sampled from air contained mixed templates which reflect the species known to be present within a confined space and that this material can be accessed using existing sampling methods. We anticipate this demonstration will initiate a much larger research programme in terrestrial airDNA sampling and that this may rapidly advance biomonitoring approaches. Lastly, we outline these and potential related applications we expect to benefit from this development.

Leech blood-meal invertebrate-derived DNA reveals differences in Bornean mammal diversity across habitats

The application of metabarcoding to environmental and invertebrate‐derived DNA (eDNA and iDNA) is a new and increasingly applied method for monitoring biodiversity across a diverse range of habitats. This approach is particularly promising for sampling in the biodiverse humid tropics, where rapid land‐use change for agriculture means there is a growing need to understand the conservation value of the remaining mosaic and degraded landscapes. Here we use iDNA from blood‐feeding leeches (Haemadipsa picta) to assess differences in mammalian diversity across a gradient of forest degradation in Sabah, Malaysian Borneo. We screened 557 individual leeches for mammal DNA by targeting fragments of the 16S rRNA gene and detected 14 mammalian genera. We recorded lower mammal diversity in the most heavily degraded forest compared to higher quality twice logged forest. Although the accumulation curves of diversity estimates were comparable across these habitat types, diversity was higher in twice logged forest, with more taxa of conservation concern. In addition, our analysis revealed differences between the community recorded in the heavily logged forest and that of the twice logged forest. By revealing differences in mammal diversity across a human‐modified tropical landscape, our study demonstrates the value of iDNA as a noninvasive biomonitoring approach in conservation assessments.

Sika deer presence affects the host-parasite interface of a Japanese land leech

Since the 1990s, increasing populations of a blood feeding land leech (Haemadipsa japonica) have become a serious issue in several Japanese prefectures, and it may be caused by the increases in sika deer (Cervus nippon) populations seen over the last quarter of the century. Therefore, this study aimed to reveal the host animal species of H. japonica using iDNA (vertebrate DNA isolated from invertebrates) and to test the hypothesis that the increasingly widespread distribution of sika deer results in increased H. japonica populations through changes to the host-parasite interface. We amplified mitochondrial DNA 16S ribosome RNA fragments from iDNA isolated from the blood clots of H. japonica collected across Japan. We identified 17 host animal species, including four orders of Mammalia (Carnivora, Artiodactyla, Rodentia, and Lagomorpha) and two orders of Amphibia (Caudata and Anura). The sika deer was the dominant host species of H. japonica. Additionally, the host animal species composition of H. japonica differed according to the presence or absence of sika deer. In the sites where sika deer were not found, Anura (frog) species were the most commonly identified hosts of H. japonica. These results suggest that the increases in H. japonica populations might have occurred via a change in host preference to sika deer. This change might be driven by the increases in sika deer populations and subsequent increase in the frequency that H. japonica uses the sika deer as easy prey, as well as by sika deer providing more reproductive energy per blood meal than blood meal from frog species. The present study suggests that a more widespread distribution of sika deer resulted in an increase in H. japonica through a change in the host-parasite interface. Therefore, management that focuses on decreasing sika deer populations would likely be an effective method for the reduction of H. japonica populations.

Using metabarcoding to compare the suitability of two blood-feeding leech species for sampling mammalian diversity in North Borneo

The application of high-throughput sequencing (HTS) for metabarcoding of mixed samples offers new opportunities in conservation biology. Recently, the successful detection of prey DNA from the guts of leeches has raised the possibility that these, and other blood-feeding invertebrates, might serve as useful samplers of mammals. Yet little is known about whether sympatric leech species differ in their feeding preferences, and whether this has a bearing on their relative suitability for monitoring local mammalian diversity. To address these questions, we collected spatially matched samples of two congeneric leech species Haemadipsa picta and Haemadipsa sumatrana from lowland rainforest in Borneo. For each species, we pooled ~500 leeches into batches of 10 individuals, performed PCR to target a section of the mammalian 16S rRNA locus and undertook sequencing of amplicon libraries using an Illumina MiSeq. In total, we identified sequences from 14 mammalian genera, spanning nine families and five orders. We found greater numbers of detections, and higher diversity of OTUs, in H. picta compared with H. sumatrana, with rodents only present in the former leech species. However, comparison of samples from across the landscape revealed no significant difference in mammal community composition between the leech species. We therefore suggest that H. picta is the more suitable iDNA sampler in this degraded Bornean forest. We conclude that the choice of invertebrate sampler can influence the detectability of different mammal groups and that this should be accounted for when designing iDNA studies.

Approaches to integrating genetic data into ecological networks

As molecular tools for assessing trophic interactions become common, research is increasingly focused on the construction of interaction networks. Here, we demonstrate three key methods for incorporating DNA data into network ecology and discuss analytical considerations using a model consisting of plants, insects, bats and their parasites from the Costa Rica dry forest. The simplest method involves the use of Sanger sequencing to acquire long sequences to validate or refine field identifications, for example of bats and their parasites, where one specimen yields one sequence and one identification. This method can be fully quantified and resolved and these data resemble traditional ecological networks. For more complex taxonomic identifications, we target multiple DNA loci, for example from a seed or fruit pulp sample in faeces. These networks are also well resolved but gene targets vary in resolution and quantification is difficult. Finally, for mixed templates such as faecal contents of insectivorous bats, we use DNA metabarcoding targeting two sequence lengths (157 and 407 bp) of one gene region and a MOTU, BLAST and BIN association approach to resolve nodes. This network type is complex to generate and analyse, and we discuss the implications of this type of resolution on network analysis. Using these data, we construct the first molecular-based network of networks containing 3,304 interactions between 762 nodes of eight trophic functions and involving parasitic, mutualistic and predatory interactions. We provide a comparison of the relative strengths and weaknesses of these data types in network ecology.