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.     

An improved method for utilizing high‐throughput amplicon sequencing to determine the diets of insectivorous animals

DNA analysis of predator faeces using high‐throughput amplicon sequencing (HTS) enhances our understanding of predator–prey interactions. However, conclusions drawn from this technique are constrained by biases that occur in multiple steps of the HTS workflow. To better characterize insectivorous animal diets, we used DNA from a diverse set of arthropods to assess PCR biases of commonly used and novel primer pairs for the mitochondrial gene, cytochrome oxidase C subunit 1 (COI). We compared diversity recovered from HTS of bat guano samples using a commonly used primer pair “ZBJ” to results using the novel primer pair “ANML.” To parameterize our bioinformatics pipeline, we created an arthropod mock community consisting of single‐copy (cloned) COI sequences. To examine biases associated with both PCR and HTS, mock community members were combined in equimolar amounts both pre‐ and post‐PCR. We validated our system using guano from bats fed known diets and using composite samples of morphologically identified insects collected in pitfall traps. In PCR tests, the ANML primer pair amplified 58 of 59 arthropod taxa (98%), whereas ZBJ amplified 24–40 of 59 taxa (41%–68%). Furthermore, in an HTS comparison of field‐collected samples, the ANML primers detected nearly fourfold more arthropod taxa than the ZBJ primers. The additional arthropods detected include medically and economically relevant insect groups such as mosquitoes. Results revealed biases at both the PCR and sequencing levels, demonstrating the pitfalls associated with using HTS read numbers as proxies for abundance. The use of an arthropod mock community allowed for improved bioinformatics pipeline parameterization.     

Blocking primers to enhance PCR amplification of rare sequences in mixed samples – a case study on prey DNA in Antarctic krill stomachs

Background

Identification of DNA sequence diversity is a powerful means for assessing the species present in environmental samples. The most common molecular strategies for estimating taxonomic composition depend upon PCR with universal primers that amplify an orthologous DNA region from a range of species. The diversity of sequences within a sample that can be detected by universal primers is often compromised by high concentrations of some DNA templates. If the DNA within the sample contains a small number of sequences in relatively high concentrations, then less concentrated sequences are often not amplified because the PCR favours the dominant DNA types. This is a particular problem in molecular diet studies, where predator DNA is often present in great excess of food-derived DNA.

Results

We have developed a strategy where a universal PCR simultaneously amplifies DNA from food items present in DNA purified from stomach samples, while the predator's own DNA is blocked from amplification by the addition of a modified predator-specific blocking primer. Three different types of modified primers were tested out; one annealing inhibiting primer overlapping with the 3' end of one of the universal primers, another annealing inhibiting primer also having an internal modification of five dI molecules making it a dual priming oligo, and a third elongation arrest primer located between the two universal primers. All blocking primers were modified with a C3 spacer. In artificial PCR mixtures, annealing inhibiting primers proved to be the most efficient ones and this method reduced predator amplicons to undetectable levels even when predator template was present in 1000 fold excess of the prey template. The prey template then showed strong PCR amplification where none was detectable without the addition of blocking primer. Our method was applied to identifying the winter food of one of the most abundant animals in the world, the Antarctic krill, Euphausia superba. Dietary item DNA was PCR amplified from a range of species in krill stomachs for which we had no prior sequence knowledge.

Conclusion

We present a simple, robust and cheap method that is easily adaptable to many situations where a rare DNA template is to be PCR amplified in the presence of a higher concentration template with identical PCR primer binding sites.     

Group-Specific Multiplex PCR Detection Systems for the Identification of Flying Insect Prey

The applicability of species-specific primers to study feeding interactions is restricted to those ecosystems where the targeted prey species occur. Therefore, group-specific primer pairs, targeting higher taxonomic levels, are often desired to investigate interactions in a range of habitats that do not share the same species but the same groups of prey. Such primers are also valuable to study the diet of generalist predators when next generation sequencing approaches cannot be applied beneficially. Moreover, due to the large range of prey consumed by generalists, it is impossible to investigate the breadth of their diet with species-specific primers, even if multiplexing them. However, only few group-specific primers are available to date and important groups of prey such as flying insects have rarely been targeted. Our aim was to fill this gap and develop group-specific primers suitable to detect and identify the DNA of common taxa of flying insects. The primers were combined in two multiplex PCR systems, which allow a time- and cost-effective screening of samples for DNA of the dipteran subsection Calyptratae (including Anthomyiidae, Calliphoridae, Muscidae), other common dipteran families (Phoridae, Syrphidae, Bibionidae, Chironomidae, Sciaridae, Tipulidae), three orders of flying insects (Hymenoptera, Lepidoptera, Plecoptera) and coniferous aphids within the genus Cinara. The two PCR assays were highly specific and sensitive and their suitability to detect prey was confirmed by testing field-collected dietary samples from arthropods and vertebrates. The PCR assays presented here allow targeting prey at higher taxonomic levels such as family or order and therefore improve our ability to assess (trophic) interactions with flying insects in terrestrial and aquatic habitats.     

Effectiveness of Annealing Blocking Primers versus Restriction Enzymes for Characterization of Generalist Diets: Unexpected Prey Revealed in the Gut Contents of Two Coral Reef Fish Species

Characterization of predator-prey interactions is challenging as researchers have to rely on indirect methods that can be costly, biased and too imprecise to elucidate the complexity of food webs. DNA amplification and sequencing techniques of gut and fecal contents are promising approaches, but their success largely depends on the ability to amplify the taxonomic array of prey consumed and then match prey amplicons with reference sequences. When little a priori information on diet is available or a generalist predator is targeted, versatile primer sets (also referred to as universal or general primers) as opposed to group- or species-specific primer sets are the most powerful to unveil the full range of prey consumed. However, versatile primers are likely to preferentially amplify the predominant, less degraded predator DNA if no manipulation is performed to exclude this confounding DNA template. In this study we compare two approaches that eliminate the confounding predator template: restriction digestion and the use of annealing blocking primers. First, we use a preliminary DNA barcode library provided by the Moorea BIOCODE project to 1) evaluate the cutting frequency of commercially available restriction enzymes and 2) design predator specific annealing blocking primers. We then compare the performance of the two predator removal strategies for the detection of prey templates using two versatile primer sets from the gut contents of two generalist coral reef fish species sampled in Moorea. Our study demonstrates that blocking primers should be preferentially used over restriction digestion for predator DNA removal as they recover greater prey diversity. We also emphasize that a combination of versatile primers may be required to best represent the breadth of a generalist''s diet.     

A framework for developing and validating taxon‐specific primers for specimen identification from environmental DNA

Taxon‐specific DNA tests are applied to many ecological and management questions, increasingly using environmental DNA (eDNA). eDNA facilitates noninvasive ecological studies but introduces additional risks of bias and error. For effective application, PCR primers must be developed for each taxon and validated in each system. We outline a nine step framework for the development and validation of taxon‐specific primers for eDNA analysis in ecological studies, involving reference database construction, phylogenetic evaluation of the target gene, primer design, primer evaluation in silico, and laboratory evaluation of primer specificity, sensitivity and utility. Our framework makes possible a rigorous evaluation of likely sources of error. The first five steps can be conducted relatively rapidly and (where reference DNA sequences are available) require minimal laboratory resources, enabling assessment of primer suitability before investing in further work. Steps six to eight require more costly laboratory analyses but are essential to evaluate risks of false‐positive and false‐negative results, while step 9 relates to field implementation. As an example, we have developed and evaluated primers to specifically amplify part of the mitochondrial ND2 gene from Australian bandicoots. If adopted during the early stages of primer development, our framework will facilitate large‐scale implementation of well‐designed DNA tests to detect specific wildlife from eDNA samples. This will provide researchers and managers with an understanding of the strengths and limitations of their data and the conclusions that can be drawn from them.     

Molecular identification of prey in the stomach contents of Harp Seals (Pagophilus groenlandicus) using species-specific oligonucleotides

All methods of diet analysis in marine mammals, including hard part analysis (HPA), have biases affecting the accuracy of prey-species identification and frequency in the estimated diet due to differential consumption, digestion and retention. Using PCR amplification of specific prey DNA with species-specific primers, we developed a DNA-based method that complements HPA and provides an alternative means to detect prey from stomach contents of Harp Seals (Pagophilus groenlandicus). The target size that could be reliably amplified was determined using a digestion time-series of Atlantic Cod (Gadus morhua) tissue in simulated seal stomachs. Various target lengths were trialed using general teleost primers; amplicons of approximately 800 bp or less were consistently obtained. Prey species-specific PCR primers for Atlantic Cod, Arctic Cod (Boreogadus saida) and Capelin (Mallotus villosus) were designed and tested with DNA from the stomach contents of 31 Harp Seals. Amplicons were obtained for all three species-specific primer sets. Amplification results compared with HPA revealed: (i) Atlantic Cod hard parts were found in five stomachs where no Atlantic Cod DNA amplified, suggesting that Atlantic Cod may be over-represented in the estimated diet, (ii) amplification of Arctic Cod DNA occurred for 17 stomachs, including all 12 stomachs with, and five stomachs without, Arctic Cod hard parts, and (iii) Capelin DNA amplified for four of five stomachs with Capelin hard parts and for one stomach without Capelin hard parts. We conclude that PCR amplification of specific prey DNA provides a viable means to complement Harp Seal diet analysis by HPA, but suggest that valuable information for quantitative diet analysis rests in a quantitative PCR approach.     

PCR-CTPP: a rapid and reliable genotyping technique based on ZFX/ZFY alleles for sex identification of tiger (<Emphasis Type="Italic">Panthera tigris</Emphasis>) and four other endangered felids

An inexpensive, time-saving and reliable method, polymerase chain reaction with confronting two-pair primers (PCR-CTPP), was developed for sex identification in tiger (Panthera tigris) based on zinc finger alleles (ZFX/ZFY). A site of “C/G” transversion representing fixed differences that discriminated between ZFX and ZFY exons among felids was identified for primers designing. This primer set was successfully tested on samples including blood, shed hairs, dried skin, and stool which contained potential contamination caused by prey DNA. Cross species tests shown that this primer set was also useful for sex identification in four other endangered felids.     

Faecal DNA amplification in Pacific walruses (Odobenus rosmarus divergens)

Dietary information is critical for assessing the population status of seals, sea lions and walruses—and is determined for most species of pinnipeds using non-invasive methods. However, diets of walruses continue to be described from the stomach contents of dead individuals. Our goal was to assess whether DNA could be extracted from the faeces of Pacific walruses (Odobenus rosmarus divergens) collected at haulout sites, and whether potential prey species or taxa could be amplified from that DNA. We extracted DNA from 70 faecal samples collected from ice pans in the Bering Sea during the spring of 2008 and 2009 (with between 4.6 and 308.9 ng/μl of DNA in every sample). We also extracted DNA from 12 potential prey species or taxa collected by bottom-grabs in 2009 to identify positive controls for primers and to test the ability of previously published taxon-specific and species-specific primers to correctly identify the prey using conventional PCR. We tested primers that successfully amplified DNA from the tissue of at least one potential prey species or taxon on all 70 walrus faecal samples. We found that two sets of primers successfully amplified many of the potential prey species or taxa using DNA from their tissue, and that one of these primer sets produced positive amplification in 4 of the 70 faecal samples. The band size that was produced for prey organisms and in the faecal samples was consistent with expectations, although prey identities were not verified with sequencing. Our pilot study demonstrates that DNA can be successfully extracted and amplified from walrus faeces, providing a stepping stone towards describing the diets of walruses from faecal DNA.     

Optimization of RAPD‐PCR conditions in cattle

Abstract

Random Amplified Polymorphic DNA polymerase chain reaction (RAPD‐PCR) is a fast and easy way of identifying DNA polymorphisms generated from several regions of the genome. This could expedite the process of identifying informative polymorphic markers that may be linked to important genes controlling economic traits. In cattle, failure to obtain consistent amplification patterns in RAPD‐PCR has been a cause for concern. This has been attributed to the fact that decamer primers that are used in RAPD‐PCR reactions are likely to amplify regions of DNA where the primer‐template base pairing has some degree of mismatch and that these mismatches fail to repeat from reaction to reaction. This paper describes the use of tricine buffer along with changes in reaction components and thermal cycling conditions that has yielded consistent and reproducible RAPD‐PCR amplifications using single primers and double primer combinations on bovine DNA.     

扬子鳄鞣制皮革和鳞片的DNA提取方法

运用一种改进的提取方法 ,作者从鞣制皮革中成功地提取了总DNA ,同时还对尾尖皮、鳞片、盐腌生皮等皮质进行了DNA提取 ;用 12SrRNA基因扩增的通用引物、扬子鳄鉴别引物、微卫星引物及RAPD引物进行PCR扩增 ,并对部分扩增结果进行测序 ,以检验提取效果。结果证明 ,几种皮质标本都可提取出DNA ,其中尾尖皮和鳞片的提取效果较好 ,用四种引物都可扩增出明显亮带 ;盐腌生皮和鞣制皮提取结果也很好 ,并且用12SrRNA通用引物、扬子鳄鉴别引物扩增的亮带较明显 ,可进行扬子鳄皮质用品等的分子鉴定及部分序列的扩增和测序研究     

A tale of two fluids: does storing specimens together in liquid preservative cause DNA cross‐contamination in molecular gut‐content studies?

The study of food webs and trophic interactions increasingly relies on PCR‐based molecular gut‐content analysis. However, this technique may be prone to error from contamination of minute quantities of DNA; i.e., simply storing specimens together in a liquid medium may lead to cross‐contamination. In this study, we used PCR to determine the contamination rate when (1) specimens were stored together in 95% ethanol for various time periods, and (2) predators fall into ethylene glycol‐filled pitfall traps where the dying predator may inadvertently consume prey DNA‐contaminated liquid. We designed experiments and PCR primers to quantify the risk of contamination for both situations and found no contamination by storing specimens together in 95% ethanol. Furthermore, zero predators contained prey DNA in their gut contents from imbibing prey DNA‐contaminated ethylene glycol. These results support the use of mass sampling techniques, like wet pitfall traps, for molecular gut‐content analysis.     

OpenADAM: an open source genome-wide association data management system for Affymetrix SNP arrays

Background

The goal of DNA barcoding is to develop a species-specific sequence library for all eukaryotes. A 650 bp fragment of the cytochrome c oxidase 1 (CO1) gene has been used successfully for species-level identification in several animal groups. It may be difficult in practice, however, to retrieve a 650 bp fragment from archival specimens, (because of DNA degradation) or from environmental samples (where universal primers are needed).

Results

We used a bioinformatics analysis using all CO1 barcode sequences from GenBank and calculated the probability of having species-specific barcodes for varied size fragments. This analysis established the potential of much smaller fragments, mini-barcodes, for identifying unknown specimens. We then developed a universal primer set for the amplification of mini-barcodes. We further successfully tested the utility of this primer set on a comprehensive set of taxa from all major eukaryotic groups as well as archival specimens.

Conclusion

In this study we address the important issue of minimum amount of sequence information required for identifying species in DNA barcoding. We establish a novel approach based on a much shorter barcode sequence and demonstrate its effectiveness in archival specimens. This approach will significantly broaden the application of DNA barcoding in biodiversity studies.     

Diagnostic PCR can be used to illuminate meiofaunal diets and trophic relationships

Analysis of the meiofaunal food web is hampered because few prey have features that persist long enough in a predator's digestive tract to allow identification to species. Hence, at least for platyhelminth predators, direct observations of prey preference are almost nonexistent, and where they occur, prey identification is often limited to the phylum level. Studies using an in vitro approach are rare because they are extremely time‐consuming and are subject to the criticism that predators removed from their natural environment may exhibit altered behaviors. Although PCR‐based approaches have achieved wide application in food‐web analysis, their application to meiofaunal flatworms suffers from a number of limitations. Most importantly, the microscopic size of both the predator and prey does not allow for removal of prey material from the digestive tract of the predator, and thus the challenge is to amplify prey sequences in the presence of large quantities of predator sequence. Here, we report on the successful use of prey‐taxon‐specific primers in diagnostic PCR to identify, to species level, specific prey items of 13 species of meiofaunal flatworms. Extension of this method will allow, for the first time, the development of a species‐level understanding of trophic interactions among the meiofauna.     

General-use polymerase chain reaction primers for amplification and direct sequencing of enolase, a single-copy nuclear gene, from different animal phyla

In contrast to mitochondrial DNA, remarkably few general-use primer sets are available for single-copy nuclear genes across animal phyla. Here, we present a primer set that yields a c. 364-bp coding fragment of the metabolic gene enolase, which includes an intron in some taxa. In species where introns are absent or have few insertions/deletions, the amplified fragment can be sequenced directly for phylogenetic or population analysis. Between-species variation in the coding region occurs widely at third codon positions, even between closely related taxa, making the fragment useful for species-level systematics. In low gene-flow species, the primers may also be of use for population genetics, as intraspecific polymorphisms occur at several silent positions in the taxa examined.     

Group-specific polymerase chain reaction for DNA-based analysis of species diversity and identity in dietary samples

Unique DNA sequences are present in all species and can be used as biomarkers for the detection of cells from that species. These DNA sequences can most easily be detected using the polymerase chain reaction (PCR), which allows very small quantities of target DNA sequence to be amplified even when the target is mixed with large amounts of nontarget DNA. PCR amplification of DNA markers that are present in a wide range of species has proven very useful for studies of species diversity in environmental samples. The taxonomic range of species to be identified from environmental samples may often need to be restricted to simplify downstream analyses and to ensure that less abundant sequences are amplified. Group-specific PCR primer sets are one means of specifying the range of taxa that produce an amplicon in a PCR. We have developed a range of group-specific PCR primers for studying the prey diversity found in predator stomach contents and scats. These primers, their design and their application to studying prey diversity and identity in predator diet are described.     

Highly effective sequencing whole chloroplast genomes of angiosperms by nine novel universal primer pairs

Chloroplast genomes supply indispensable information that helps improve the phylogenetic resolution and even as organelle‐scale barcodes. Next‐generation sequencing technologies have helped promote sequencing of complete chloroplast genomes, but compared with the number of angiosperms, relatively few chloroplast genomes have been sequenced. There are two major reasons for the paucity of completely sequenced chloroplast genomes: (i) massive amounts of fresh leaves are needed for chloroplast sequencing and (ii) there are considerable gaps in the sequenced chloroplast genomes of many plants because of the difficulty of isolating high‐quality chloroplast DNA, preventing complete chloroplast genomes from being assembled. To overcome these obstacles, all known angiosperm chloroplast genomes available to date were analysed, and then we designed nine universal primer pairs corresponding to the highly conserved regions. Using these primers, angiosperm whole chloroplast genomes can be amplified using long‐range PCR and sequenced using next‐generation sequencing methods. The primers showed high universality, which was tested using 24 species representing major clades of angiosperms. To validate the functionality of the primers, eight species representing major groups of angiosperms, that is, early‐diverging angiosperms, magnoliids, monocots, Saxifragales, fabids, malvids and asterids, were sequenced and assembled their complete chloroplast genomes. In our trials, only 100 mg of fresh leaves was used. The results show that the universal primer set provided an easy, effective and feasible approach for sequencing whole chloroplast genomes in angiosperms. The designed universal primer pairs provide a possibility to accelerate genome‐scale data acquisition and will therefore magnify the phylogenetic resolution and species identification in angiosperms.     

穿山甲标本和甲片的DNA提取及PCR扩增

为验证经处理后的穿山甲(Manis spp.)标本和甲片是否可以用于种间分子鉴定标记的开发及个体识别工作,本文在样品的预处理、消化、提取后纯化等方面对传统提取方法进行了改进,分别从穿山甲剥制标本、干皮标本及甲片中提取总DNA;然后用Cyt b基因扩增通用引物、12S rRNA基因全序列扩增引物、RAPD引物及微卫星引物进行了PCR扩增,并对部分扩增结果进行了序列测定.结果表明,除剥制标本的脚底皮张组织外,其他样品基本都可以提取出DNA.以此为模板的PCR扩增中,2种线粒体基因引物扩增出明显目的条带,RAPD引物扩增出种间特异条带,测序结果可用于种间特异性引物及SCAR引物的开发;微卫星引物在甲片样品中扩增稳定,可用于个体识别工作.     

Improving PCR detection of prey in molecular diet studies: importance of group‐specific primer set selection and extraction protocol performances

While the morphological identification of prey remains in predators' faeces is the most commonly used method to study trophic interactions, many studies indicate that this method does not detect all consumed prey. Polymerase chain reaction–based methods are increasingly used to detect prey DNA in the predator food bolus and have proven efficient, delivering highly accurate results. When studying complex diet samples, the extraction of total DNA is a critical step, as polymerase chain reaction (PCR) inhibitors may be co‐extracted. Another critical step involves a careful selection of suitable group‐specific primer sets that should only amplify DNA from the targeted prey taxon. In this study, the food boluses of five Rattus rattus and seven Rattus exulans were analysed using both morphological and molecular methods. We tested a panel of 31 PCR primer pairs targeting bird, invertebrate and plant sequences; four of them were selected to be used as group‐specific primer pairs in PCR protocols. The performances of four DNA extraction protocols (QIAamp^® DNA stool mini kit, DNeasy^® mericon food kit and two of cetyltrimethylammonium bromide‐based methods) were compared using four variables: DNA concentration, A₂₆₀/A₂₈₀ absorbance ratio, food compartment analysed (stomach or faecal contents) and total number of prey‐specific PCR amplification per sample. Our results clearly indicate that the A₂₆₀/A₂₈₀ absorbance ratio, which varies between extraction protocols, is positively correlated to the number of PCR amplifications of each prey taxon. We recommend using the DNeasy^® mericon food kit (QIAGEN), which yielded results very similar to those achieved with the morphological approach.     

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.