Mammalian species identification by interspersed repeat PCR fingerprinting

Most DNA methods for species identification of animal tissues test the presence/absence of one species per assay, requiring several tests for a complete analysis and prior knowledge of the species that are potentially present in the sample. Here we demonstrate that PCR with fluorescently labeled MIR (mammalian-wide interspersed repeat) primers generate fingerprints that are suitable for rapid identification of known and unknown species on an automatic sequencing apparatus and with computer-assisted data processing. The method allows the analysis of processed meat samples and offers a convenient alternative to sequencing of mitochondrial DNA. Received 19 December 1997/ Accepted in revised form 15 June 1998     

降香黄檀的DNA条形码鉴定研究

DNA条形码技术是利用基因组中一段短的标准序列进行物种的鉴定并探索其亲缘进化关系。本研究对采自海南不同地区降香黄檀五个居群24份样品的psbA-trnH,rbcL,核ITS及ITS2序列进行PCR扩增和测序,比较各序列扩增和测序效率。种间和种内变异,采用BLAST1和邻接 (NJ) 法构建系统聚类树方法评价不同序列的鉴定能力。结果表明ITS2在所研究的材料中具有最高的扩增和测序效率,而ITS扩增效率较低。ITS2完整序列在区分黄檀属不同种间差异具有较大优势。因此可利用ITS2从分子水平区分降香黄檀与其他混伪种。     

Successful carnivore identification with faecal DNA across a fragmented Amazonian landscape

The use of scat surveys to obtain DNA has been well documented in temperate areas, where DNA preservation may be more effective than in tropical forests. Samples obtained in the tropics are often exposed to high humidity, warm temperatures, frequent rain and intense sunlight, all of which can rapidly degrade DNA. Despite these potential problems, we demonstrate successful mtDNA amplification and sequencing for faeces of carnivores collected in tropical conditions and quantify how sample condition and environmental variables influence the success of PCR amplification and species identification. Additionally, the feasibility of genotyping nuclear microsatellites from jaguar (Panthera onca) faeces was investigated. From October 2007 to December 2008, 93 faecal samples were collected in the southern Brazilian Amazon. A total of eight carnivore species was successfully identified from 71% of all samples obtained. Information theoretic analysis revealed that the number of PCR attempts before a successful sequence was an important negative predictor across all three responses (success of species identification, success of species identification from the first sequence and PCR amplification success), whereas the relative importance of the other three predictors (sample condition, season and distance from forest edge) varied between the three responses. Nuclear microsatellite amplification from jaguar faeces had lower success rates (15-44%) compared with those of the mtDNA marker. Our results show that DNA obtained from faecal samples works efficiently for carnivore species identification in the Amazon forest and also shows potential for nuclear DNA analysis, thus providing a valuable tool for genetic, ecological and conservation studies.     

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.     

Novel primers for complete mitochondrial cytochrome b gene sequencing in mammals

Sequence-based species identification relies on the extent and integrity of sequence data available in online databases such as GenBank. When identifying species from a sample of unknown origin, partial DNA sequences obtained from the sample are aligned against existing sequences in databases. When the sequence from the matching species is not present in the database, high-scoring alignments with closely related sequences might produce unreliable results on species identity. For species identification in mammals, the cytochrome b (cyt b) gene has been identified to be highly informative; thus, large amounts of reference sequence data from the cyt b gene are much needed. To enhance availability of cyt b gene sequence data on a large number of mammalian species in GenBank and other such publicly accessible online databases, we identified a primer pair for complete cyt b gene sequencing in mammals. Using this primer pair, we successfully PCR amplified and sequenced the complete cyt b gene from 40 of 44 mammalian species representing 10 orders of mammals. We submitted 40 complete, correctly annotated, cyt b protein coding sequences to GenBank. To our knowledge, this is the first single primer pair to amplify the complete cyt b gene in a broad range of mammalian species. This primer pair can be used for the addition of new cyt b gene sequences and to enhance data available on species represented in GenBank. The availability of novel and complete gene sequences as high-quality reference data can improve the reliability of sequence-based species identification.     

The impact of contaminants on the accuracy of genome skimming and the effectiveness of exclusion read filters

The ability to detect the identity of a sample obtained from its environment is a cornerstone of molecular ecological research. Thanks to the falling price of shotgun sequencing, genome skimming, the acquisition of short reads spread across the genome at low coverage, is emerging as an alternative to traditional barcoding. By obtaining far more data across the whole genome, skimming has the promise to increase the precision of sample identification beyond traditional barcoding while keeping the costs manageable. While methods for assembly‐free sample identification based on genome skims are now available, little is known about how these methods react to the presence of DNA from organisms other than the target species. In this paper, we show that the accuracy of distances computed between a pair of genome skims based on k‐mer similarity can degrade dramatically if the skims include contaminant reads; i.e., any reads originating from other organisms. We establish a theoretical model of the impact of contamination. We then suggest and evaluate a solution to the contamination problem: Query reads in a genome skim against an extensive database of possible contaminants (e.g., all microbial organisms) and filter out any read that matches. We evaluate the effectiveness of this strategy when implemented using Kraken‐II, in detailed analyses. Our results show substantial improvements in accuracy as a result of filtering but also point to limitations, including a need for relatively close matches in the contaminant database.     

SPInDel: a multifunctional workbench for species identification using insertion/deletion variants

The majority of the available methods for the molecular identification of species use pairwise sequence divergences between the query and reference sequences (DNA barcoding). The presence of multiple insertions and deletions (indels) in the target genomic regions is generally regarded as a problem, as it introduces ambiguities in sequence alignments. However, we have recently shown that a high level of species discrimination is attainable in all taxa of life simply by considering the length of hypervariable regions defined by indel variants. Each species is tagged with a numeric profile of fragment lengths—a true numeric barcode. In this study, we describe a multifunctional computational workbench (named SPInDel for SPecies Identification by Insertions/Deletions) to assist researchers using variable‐length DNA sequences, and we demonstrate its applicability in molecular ecology. The SPInDel workbench provides a step‐by‐step environment for the alignment of target sequences, selection of informative hypervariable regions, design of PCR primers and the statistical validation of the species‐identification process. In our test data sets, we were able to discriminate all species from two genera of frogs (Ansonia and Leptobrachium) inhabiting lowland rainforests and mountain regions of South‐East Asia and species from the most common genus of coral reef fishes (Apogon). Our method can complement conventional DNA barcoding systems when indels are common (e.g. in rRNA genes) without the required step of DNA sequencing. The executable files, source code, documentation and test data sets are freely available at http://www.portugene.com/SPInDel/SPInDel_webworkbench.html .     

Phylogenomic relationships and species identification of the olive genus Olea (Oleaceae)

The olive genus Olea includes c. 30–40 taxa in three subgenera (Olea, Tetrapilus, and Paniculatae) within the family Oleaceae. Historically, the Olea genus was classified into four groups that were overall well supported by reconstructed phylogenies, despite incomplete sampling of subgenus Tetrapilus and poor resolution within clades. These analyses also showed that the genus was not monophyletic. Reliable identification of Olea species is important for both their conservation and utilization of this economically important genus. In this study, we used phylogenomic data from genome skimming to resolve relationships within Olea and to identify molecular markers for species identification. We assembled the complete plastomes, and nrDNA of 26 individuals representing 13 species using next-generation sequencing and added 18 publicly available accessions of Olea. We also developed nuclear SNPs using the genome skimming data to infer the phylogenetic relationships of Olea. Large-scale phylogenomic analyses of 138 samples of tribe Oleeae supported the polyphyly of Olea, with Olea caudatilimba and Olea subgenus Tetrapilus not sharing their most recent common ancestor with the main Olea clade (subgenus Paniculatae and subgenus Olea). The interspecific phylogenetic resolution was poor owing to a possible rapid radiation. By comparing with the plastome data, we identified the markers ycf1b and psbE-petL as the best Olea-specific chloroplast DNA barcodes. Compared with universal barcodes, specific DNA barcodes and super-barcode exhibited higher discriminatory power. Our results demonstrated the power of phylogenomics to improve phylogenetic relationships of intricate groups and provided new insights into barcodes that allow for accurate identification of Olea species.     

DNA barcoding of oomycetes with cytochrome c oxidase subunit I and internal transcribed spacer

Oomycete species occupy many different environments and many ecological niches. The genera Phytophthora and Pythium for example, contain many plant pathogens which cause enormous damage to a wide range of plant species. Proper identification to the species level is a critical first step in any investigation of oomycetes, whether it is research driven or compelled by the need for rapid and accurate diagnostics during a pathogen outbreak. The use of DNA for oomycete species identification is well established, but DNA barcoding with cytochrome c oxidase subunit I (COI) is a relatively new approach that has yet to be assessed over a significant sample of oomycete genera. In this study we have sequenced COI, from 1205 isolates representing 23 genera. A comparison to internal transcribed spacer (ITS) sequences from the same isolates showed that COI identification is a practical option; complementary because it uses the mitochondrial genome instead of nuclear DNA. In some cases COI was more discriminative than ITS at the species level. This is in contrast to the large ribosomal subunit, which showed poor species resolution when sequenced from a subset of the isolates used in this study. The results described in this paper indicate that COI sequencing and the dataset generated are a valuable addition to the currently available oomycete taxonomy resources, and that both COI, the default DNA barcode supported by GenBank, and ITS, the de facto barcode accepted by the oomycete and mycology community, are acceptable and complementary DNA barcodes to be used for identification of oomycetes.     

TERAD: Extraction of transposable element composition from RADseq data

Transposable elements (TEs) – selfish DNA sequences that can move within the genome – comprise a large proportion of the genomes of many organisms. Although low‐coverage whole‐genome sequencing can be used to survey TE composition, it is noneconomical for species with large quantities of DNA. Here, we utilize restriction‐site associated DNA sequencing (RADSeq) as an alternative method to survey TE composition. First, we demonstrate in silico that double digest restriction‐site associated DNA sequencing (ddRADseq) markers contain the same TE compositions as whole genome assemblies across arthropods. Next, we show empirically using eight Synalpheus snapping shrimp species with large genomes that TE compositions from ddRADseq and low‐coverage whole‐genome sequencing are comparable within and across species. Finally, we develop a new bioinformatic pipeline, TERAD, to extract TE compositions from RADseq data. Our study expands the utility of RADseq to study the repeatome, making comparative studies of genome structure for species with large genomes more tractable and affordable.     

Chloroplast genome sequences from total DNA for plant identification

Chloroplast DNA sequence data are a versatile tool for plant identification or barcoding and establishing genetic relationships among plant species. Different chloroplast loci have been utilized for use at close and distant evolutionary distances in plants, and no single locus has been identified that can distinguish between all plant species. Advances in DNA sequencing technology are providing new cost‐effective options for genome comparisons on a much larger scale. Universal PCR amplification of chloroplast sequences or isolation of pure chloroplast fractions, however, are non‐trivial. We now propose the analysis of chloroplast genome sequences from massively parallel sequencing (MPS) of total DNA as a simple and cost‐effective option for plant barcoding, and analysis of plant relationships to guide gene discovery for biotechnology. We present chloroplast genome sequences of five grass species derived from MPS of total DNA. These data accurately established the phylogenetic relationships between the species, correcting an apparent error in the published rice sequence. The chloroplast genome may be the elusive single‐locus DNA barcode for plants.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.     

Assessing DNA for fish identifications from reference collections: the good,bad and ugly shed light on formalin fixation and sequencing approaches

Natural history collections are repositories of biodiversity and are potentially used by molecular ecologists for comparative taxonomic, phylogenetic, biogeographic and forensic purposes. Specimens in fish collections are preserved using a combination of methods with many fixed in formalin and then preserved in ethanol for long-term storage. Formalin fixation damages DNA, thereby limiting genetic analyses. In this study, the authors compared the DNA barcoding and identification success for frozen and formalin-fixed tissues obtained from specimens in the CSIRO Australian National Fish Collection. They studied 230 samples from fishes (consisting of >160 fish species). An optimized formalin-fixed, paraffin-embedded DNA extraction method resulted in usable DNA from degraded tissues. Four mini barcoding assays of the mitochondrial DNA (mtDNA) were characterized with Sanger and Illumina amplicon sequencing. In the good quality DNA (without exposure to formalin), up to 88% of the specimens were correctly matched at the species level using the cytochrome oxidase subunit 1 (COI) mini barcodes, whereas up to 58% of the specimens exposed to formalin for less than 8 weeks were correctly identified to species. In contrast, 16S primers provided higher amplification success with formalin-exposed tissues, although the COI gene was more successful for identification. Importantly, the authors found that DNA of a certain size and quality can be amplified and sequenced despite exposure to formalin, and Illumina sequencing provided them with greater power of resolution for taxa identification even when there was little DNA present. Overall, within parameter constraints, this study highlights the possibilities of recovering DNA barcodes for identification from formalin-fixed fish specimens, and the authors provide guidelines for when successful identification could be expected.     

Rapid multiplex PCR based species identification of wild tigers using non-invasive samples

Conservation and management of rare and elusive species requires accurate data on presence or absence. In such cases, molecular genetics based species identification approaches can prove invaluable, especially in conjuncture with non-invasive DNA sampling. However, non-invasive sources yield DNA in low concentration that is degraded, which could result in false negatives for species identification. In this paper, we developed a set of primers for PCR-based species identification of tigers. Our results reveal high rates (upto 90%) of species identification for both fresh (less than 48 h) and old (between 7 days and 3 months) fecal samples from the field. Experiments reveal that multiplex PCR (amplifying more than one genomic region) results in an increase in conclusive species identification (and a consequent decrease in the number of false negatives) from 55% to 89% for old fecal samples. We demonstrate that this increased success is because we experimentally overcome the problems of low DNA template quantity (using the multiplex PCR kit, increases species identification from 55% to 72%) and low template DNA quality (two sets of primers increase the species identification success from 72% to 89%). We recommend that multiplex PCR based methods be used (in conjuncture with species specific primers) for other rare and elusive species since such methods will potentially significantly decrease error in species identification.     

All-Food-Seq (AFS): a quantifiable screen for species in biological samples by deep DNA sequencing

Background

DNA-based methods like PCR efficiently identify and quantify the taxon composition of complex biological materials, but are limited to detecting species targeted by the choice of the primer assay. We show here how untargeted deep sequencing of foodstuff total genomic DNA, followed by bioinformatic analysis of sequence reads, facilitates highly accurate identification of species from all kingdoms of life, at the same time enabling quantitative measurement of the main ingredients and detection of unanticipated food components.

Results

Sequence data simulation and real-case Illumina sequencing of DNA from reference sausages composed of mammalian (pig, cow, horse, sheep) and avian (chicken, turkey) species are able to quantify material correctly at the 1% discrimination level via a read counting approach. An additional metagenomic step facilitates identification of traces from animal, plant and microbial DNA including unexpected species, which is prospectively important for the detection of allergens and pathogens.

Conclusions

Our data suggest that deep sequencing of total genomic DNA from samples of heterogeneous taxon composition promises to be a valuable screening tool for reference species identification and quantification in biosurveillance applications like food testing, potentially alleviating some of the problems in taxon representation and quantification associated with targeted PCR-based approaches.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-639) contains supplementary material, which is available to authorized users.     

Detecting host–parasitoid interactions in an invasive Lepidopteran using nested tagging DNA metabarcoding

Determining the host–parasitoid interactions and parasitism rates for invasive species entering novel environments is an important first step in assessing potential routes for biocontrol and integrated pest management. Conventional insect rearing techniques followed by taxonomic identification are widely used to obtain such data, but this can be time‐consuming and prone to biases. Here, we present a next‐generation sequencing approach for use in ecological studies which allows for individual‐level metadata tracking of large numbers of invertebrate samples through the use of hierarchically organised molecular identification tags. We demonstrate its utility using a sample data set examining both species identity and levels of parasitism in late larval stages of the oak processionary moth (Thaumetopoea processionea—Linn. 1758), an invasive species recently established in the United Kingdom. Overall, we find that there are two main species exploiting the late larval stages of oak processionary moth in the United Kingdom with the main parasitoid (Carcelia iliaca—Ratzeburg, 1840) parasitising 45.7% of caterpillars, while a rare secondary parasitoid (Compsilura concinnata—Meigen, 1824) was also detected in 0.4% of caterpillars. Using this approach on all life stages of the oak processionary moth may demonstrate additional parasitoid diversity. We discuss the wider potential of nested tagging DNA metabarcoding for constructing large, highly resolved species interaction networks.     

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build‐up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double‐digest restriction site‐associated DNA sequencing (ddRADseq) to obtain and analyse genomic‐scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.     

DNA barcoding of Pedicularis L. (Orobanchaceae): Evaluating four universal barcode loci in a large and hemiparasitic genus

Abstract One application of DNA barcoding is species identification based on sequences of a short and standardized DNA region. In plants, various DNA regions, alone or in combination, have been proposed and investigated, but consensus on a universal plant barcode remains elusive. In this study, we tested the utility of four candidate barcoding regions (rbcL, matK, trnH‐psbA, and internal transcribed spacer (ITS)) as DNA barcodes for discriminating species in a large and hemiparasitic genus Pedicularis (Orobanchaceae). Amplification and sequencing was successful using single primer pairs for rbcL, trnH‐psbA, and ITS, whereas two primer pairs were required for matK. Patterns of sequence divergence commonly showed a “barcoding gap”, that is, a bimodal frequency distribution of pairwise distances representing genetic diversity within and between species, respectively. Considering primer universality, ease of amplification and sequencing, and performance in discriminating species, we found the most effective single‐region barcode for Pedicularis to be ITS, and the most effective two‐region barcode to be rbcL + ITS. Both discriminated at least 78% of the 88 species and correctly identified at least 89% of the sequences in our sample, and were effective in placing unidentified samples in known species groups. Our results suggest that DNA barcoding has the potential to aid taxonomic research in Pedicularis, a species‐rich cosmopolitan clade much in need of revision, as well as ecological studies in its center of diversity, the Hengduan Mountains region of China.