Faecal genetic analysis to determine the presence and distribution of elusive carnivores: design and feasibility for the Iberian lynx

Noninvasive methods using genetic markers have been suggested as ways to overcome difficulties associated with documenting the presence of elusive species. We present and assess a novel, reliable and effective molecular genetic technique for the unequivocal genetic identification of faeces from the endangered Iberian lynx (Lynx pardinus). From mitochondrial DNA (mtDNA) cytochrome b and D-loop region sequences, we designed four species-specific primers (for products 130-161 bp long) that were considered to be likely to amplify degraded DNA. We compared two DNA extraction methods, various DNA amplification conditions and the robustness and specificity of the primer pairs with 87 lynx samples from 5 potentially different lynx populations and with 328 samples of other carnivore species. The utility of the identification technique was tested with faeces of different ages, with faeces from controlled field experiments, and with faeces collected from locales with possible lynx populations from throughout the state of Andalusia, Spain (8052 km2). Faecal mtDNA extraction was more efficient using PBS wash of the faeces instead of a faeces homogenate. Our assay increased from 92.6 to 99% efficiency with a second amplification and a reduction in template concentration to overcome polymerase chain reaction (PCR) inhibition. Our assay never produced false positives, and correctly identified all lynx faeces. Of 252 faeces samples of unknown species collected throughout Andalusia, 26.6% (from three different areas) were classified as Iberian lynx, 1.4% showed evidence of PCR inhibition and 1.2% were of uncertain origin. This method has proven to be a reliable technique that can be incorporated into large-scale surveys of Iberian lynx populations and exemplifies an approach that can easily be extended to other species.     

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.     

成都地区四种食尸性蝇类 mtDNA中 COⅠ基因序列检测

通过检测食尸性苍蝇线粒体DNA(mtDNA)上细胞色素氧化酶亚基Ⅰ(COⅠ)中278bp 基因序列,鉴定食尸性苍蝇的种类,解决依据形态学方法不能鉴定苍蝇卵的种类、很难鉴定幼虫种类的难题 ,作为法医鉴别食尸性苍蝇及其幼虫、卵种类依据。随机采集放置在成都地区室外草地兔尸体 上的4 种15个食尸性苍蝇。利用改进的小型昆虫DNA匀浆方法提取上述苍蝇mtDNA;通过Perkin Elmer 9600扩增仪进行PCR扩增;聚丙烯酰胺非变性凝胶连续缓冲体系垂直电泳和银染显色技 术进行扩增结果检测;PCR胶回收试剂盒纯化;ABI 377测序仪测序;MEGA2.1软件包进行序列 分析和构建系统发育树。在双翅目食尸性苍蝇的种内进化分歧均数小于1%,种间进化分歧均数 大于7%。mtDNA上COⅠ序列分析能有效地对主要的食尸性苍蝇进行种类鉴定。该检测方法快速、 简便和精确,能作为法医鉴别食尸性苍蝇种类的可靠依据。     

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.     

Differentiating between Steller sea lion (Eumetopias jubatus) and northern fur seal (Callorhinus ursinus) scats through analysis of faecal DNA

We describe a method to determine the species of pinniped from faeces collected from sympatric Steller sea lion (Eumetopias jubatus) and northern fur seal (Callorhinus ursinus) rookeries using newly developed species-specific primers that amplify a 667-669-base pair segment from the mitochondrial DNA (mtDNA) cytochrome B (cytB) gene region. The primers yielded the correct species in 100% of tissue samples from 10 known animals and 100% of faecal samples from 13 known animals. Species could be identified unequivocally for 87.7% of faecal samples from 122 unknown individuals. The ability to differentiate between scats of sympatrically breeding Steller sea lions and northern fur seals will contribute to the range-wide knowledge of the foraging strategies of both species as well as allow researchers to examine the niche partitioning and potential resource competition between the two predators.     

Simple identification of mitochondrial lineages in contact zones based on lineage-selective primers

A variety of research projects focus on genetic variation among and within maternal lineages as encompassed by mitochondrial DNA (mtDNA). While mtDNA often differs substantially between species, large differences may also be found within species. The evaluation of such divergent lineages, for example in intraspecific contact zones (hybrid zones), commonly involves sequencing numerous individuals. Large‐scale sequencing is both expensive and labour‐intensive. Based on sequences from 15 individuals, we devised a simple and quick polymerase chain reaction assay for identification of divergent mtDNA lineages in a secondary contact zone of the side‐blotched lizard (Uta stansburiana). The application uses lineage‐selective primers to amplify a lineage‐diagnostic product, and is based on each group of mtDNA haplotypes being a monophyletic assemblage of haplotypes sharing the same maternal ancestry, deeply divergent from the other group. The assay was tested on a larger sample (n = 147) of specimens from the contact zone, confirming its usefulness in quick and reliable identification of mtDNA lineages. This approach can be modified for other species, provided diagnostic lineage variation is available, and may also be performed in simple laboratory settings while conducting fieldwork.     

Molecular identification of crocodile species using novel primers for forensic analysis

All crocodilians are under varying degrees of threat due to over exploitation and these species have been listed in Appendix I or II of CITES. The lack of molecular techniques for the identification of confiscated samples makes it difficult to enforce the law. Conclusive forensic identification of species requires a complete gene sequence which is difficult in case of degraded samples. We have developed two novel sets of primers to amplify two partial cytochrome b gene sequences of six crocodile species i.e. Crocodylus palustris, Crocodylus porosus, Crocodylus siamensis, Crocodylus niloticus, Gavialis gangeticus and Caiman crocodilus. These partial sequences were edited to give a complete cyt b gene sequence, which can be used as an effective tool for forensic authentication of crocodile species. A phylogeny of crocodile species was reconstructed using these sequences. The described primers hold great promise in forensic identification of crocodile species, which can aid in the effective enforcement of law and conservation of these ancient species.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Sliding window analyses for optimal selection of mini-barcodes, and application to 454-pyrosequencing for specimen identification from degraded DNA

DNA barcoding remains a challenge when applied to diet analyses, ancient DNA studies, environmental DNA samples and, more generally, in any cases where DNA samples have not been adequately preserved. Because the size of the commonly used barcoding marker (COI) is over 600 base pairs (bp), amplification fails when the DNA molecule is degraded into smaller fragments. However, relevant information for specimen identification may not be evenly distributed along the barcoding region, and a shorter target can be sufficient for identification purposes. This study proposes a new, widely applicable, method to compare the performance of all potential 'mini-barcodes' for a given molecular marker and to objectively select the shortest and most informative one. Our method is based on a sliding window analysis implemented in the new R package SPIDER (Species IDentity and Evolution in R). This method is applicable to any taxon and any molecular marker. Here, it was tested on earthworm DNA that had been degraded through digestion by carnivorous landsnails. A 100 bp region of 16 S rDNA was selected as the shortest informative fragment (mini-barcode) required for accurate specimen identification. Corresponding primers were designed and used to amplify degraded earthworm (prey) DNA from 46 landsnail (predator) faeces using 454-pyrosequencing. This led to the detection of 18 earthworm species in the diet of the snail. We encourage molecular ecologists to use this method to objectively select the most informative region of the gene they aim to amplify from degraded DNA. The method and tools provided here, can be particularly useful (1) when dealing with degraded DNA for which only small fragments can be amplified, (2) for cases where no consensus has yet been reached on the appropriate barcode gene, or (3) to allow direct analysis of short reads derived from massively parallel sequencing without the need for bioinformatic consolidation.     

虎物种特异性鉴定的 PCR 方法研究

为了建立可对虎DNA 进行特异性检测的PCR 方法, 应用在野外调查中获得的虎疑似样品和保护执法工作中难以检查辨认的虎产品进行物种鉴定, 从Genbank 数据库下载5 个虎亚种及其它6 种猫科动物和6 种鹿科动物的mtDNA 细胞色素b 基因序列, 并用Wdnasis (V2.5) 软件对上述不同动物的该基因碱基序列进行比较。在此基础上, 综合考虑了设计引物的基本原则, 选择了虎与其它动物碱基序列上差异位点较多的两个片段, 设计出PCR 引物(引物1、2) 。用该对引物分别对从东北虎、华南虎及8 种猫科动物和6 种非猫科动物的肌肉、脏器组织、皮或毛发中提取的DNA 进行PCR (聚合酶链反应) 扩增。结果表明, 所设计的引物对虎DNA 具有特异性,从而达到了对该物种进行特异性检测鉴定的目的。     

Species-specific mitochondrial DNA markers for identification of non-invasive samples from sympatric carnivores in the Iberian Peninsula

Genetic species identification of non-invasively collected samples has become an important tool in ecological research, management and conservation and wildlife forensics. This is especially true for carnivores, due to their elusive nature, and is crucial when several ecologically and phylogenetically close species, with similar faeces, hairs, bones and/or pelts, occur in sympatry. This is the case of the Iberian Peninsula, a region with a carnivore community of 16 species—about two-thirds of the European carnivore fauna. Here we present a simple, efficient and reliable PCR-based protocol, using a novel set of species-specific primers, for the unambiguous identification to species of non-invasively collected samples or forensic materials from Iberian carnivores. For each species, from the consensus of all cytochrome b haplotypes, found here and previously reported, we designed species-specific primer pairs for short fragments, the most likely to persist in low-quantity and degraded DNA samples. The predicted specificity of each primer pair was assessed through PCR of positive DNA extracts from the carnivore species, from an exhaustive array of potential prey and from humans. The robustness of PCR amplification for non-invasively sampled DNA was tested with scat samples. The primers did not produce false positives and correctly identified all carnivore samples to the species level. In comparison with sequencing and PCR-RFLP assays, our method is, respectively, cost- and time-effective, and is especially suited for monitoring surveys targeting multiple populations/species. It also introduces an approach that works for a whole community of carnivores living sympatrically over a large geographic area.     

A modified mini-primer set for analyzing mitochondrial DNA control region sequences from highly degraded forensic samples

To facilitate the analysis of mitochondrial DNA (mtDNA) control region sequences from highly degraded skeletal remains, a modified mini-primer set was designed to overcome the limitations of the Armed Forces DNA Identification Laboratory (AFDIL) mini-primer set. This modified mini-primer set is less affected by nucleotide variability and PCR amplification conditions than the AFDIL mini-primer set, and was able to amplify the mtDNA sequences of 55-year-old skeletal remains with high efficiency, indicating that it is a useful tool for analyzing mtDNA control region sequences from highly degraded forensic samples.     

Detection of species-specific long VNTRs in mitochondrial control region and their application to identifying sympatric Hong Kong grouper (Epinephelus akaara) and yellow grouper (Epinephelus awoara)

Length variation and heteroplasmy were observed in PCR products of the first half of mtDNA control region of both Hong Kong grouper (Epinephelus akaara) and yellow grouper (Epinephelus awoara). DNA sequencing unveiled the phenomena were caused by the presence of species-specific long variable number tandem repeats (VNTRs). This is the first report on the mtDNA VNTRs and their heteroplasmy in groupers. Moreover, these VNTRs are also the longest such structure found in teleost fish. Thereafter, we designed two species-specific PCR reverse primers according to the 3' end sequences of the VNTRs and successfully established assays for the identification of these two sympatric grouper species.     

Conserved PCR primer set designing for closely-related species to complete mitochondrial genome sequencing using a sliding window-based PSO algorithm

Background

Complete mitochondrial (mt) genome sequencing is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. For long template sequencing, i.e., like the entire mtDNA, it is essential to design primers for Polymerase Chain Reaction (PCR) amplicons which are partly overlapping each other. The presented chromosome walking strategy provides the overlapping design to solve the problem for unreliable sequencing data at the 5′ end and provides the effective sequencing. However, current algorithms and tools are mostly focused on the primer design for a local region in the genomic sequence. Accordingly, it is still challenging to provide the primer sets for the entire mtDNA.

Methodology/Principal Findings

The purpose of this study is to develop an integrated primer design algorithm for entire mt genome in general, and for the common primer sets for closely-related species in particular. We introduce ClustalW to generate the multiple sequence alignment needed to find the conserved sequences in closely-related species. These conserved sequences are suitable for designing the common primers for the entire mtDNA. Using a heuristic algorithm particle swarm optimization (PSO), all the designed primers were computationally validated to fit the common primer design constraints, such as the melting temperature, primer length and GC content, PCR product length, secondary structure, specificity, and terminal limitation. The overlap requirement for PCR amplicons in the entire mtDNA is satisfied by defining the overlapping region with the sliding window technology. Finally, primer sets were designed within the overlapping region. The primer sets for the entire mtDNA sequences were successfully demonstrated in the example of two closely-related fish species. The pseudo code for the primer design algorithm is provided.

Conclusions/Significance

In conclusion, it can be said that our proposed sliding window-based PSO algorithm provides the necessary primer sets for the entire mt genome amplification and sequencing.     

核糖体转录间隔子2应用于鱼类种属的鉴别

为了防止珍稀鱼类的非法捕捞和销售, 鱼类种属的鉴别就成为非常关键的问题, 特别是形态学方法无法区分的样品(如鱼苗、鱼鳞、鱼卵、鱼肉及其加工产品等)。为了帮助珍稀鱼类资源的管理和保护, 文章报道了一种利用核糖体基因的转录间隔子2鉴别鱼类种属的分子遗传学方法: (1) 利用同一目鱼类5.8S rRNA和28S rRNA基因的保守性, 设计出扩增鲤形目鱼类这两个基因间转录间隔子2 DNA片段, 测序获得它们的碱基排列顺序; (2) 再根据不同鱼类转录间隔子2序列的差异, 设计出每种鱼的种属特异引物、种属鉴别标准物, 构建鱼类分子分类图谱, 利用PCR复合扩增技术鉴别鱼类种属。通过对国内不同地方采集的5种鲤形目鱼类的210个单一品种样本和40个混合样本的鉴别检验, 该方法能够准确、灵敏和快速鉴别这5种鱼, 可用于鱼类资源保护和评估、管理和开发, 特别是在渔业管理人员渔业执法、海关打击珍稀鱼类走私、防止商业欺诈和外来有害生物入侵等方面非常有用     

Cross-species amplification of mitochondrial DNA sequence-tagged-site markers in conifers: the nature of polymorphism and variation within and among species in Picea

Primers previously developed to amplify specific non-coding regions of the mitochondrial genome in Angiosperms, and new primers for additional non-coding mtDNA regions, were tested for their ability to direct DNA amplification in 12 conifer taxa and to detect sequence-tagged-site (STS) polymorphisms within and among eight species in Picea. Out of 12 primer pairs, nine were successful at amplifying mtDNA in most of the taxa surveyed. In conifers, indels and substitutions were observed for several loci, allowing them to distinguish between families, genera and, in some cases, between species within genera. In Picea, interspecific polymorphism was detected for four loci, while intraspecific variation was observed for three of the mtDNA regions studied. One of these (SSU rRNA V1 region) exhibited indel polymorphisms, and the two others ( nad1 intron b/c and nad5 intron1) revealed restriction differences after digestion with Sau3AI (PCR-RFLP). A fourth locus, the nad4L- orf25 intergenic region, showed a multibanding pattern for most of the spruce species, suggesting a possible gene duplication. Maternal inheritance, expected for mtDNA in conifers, was observed for all polymorphic markers except the intergenic region nad4L- orf25. Pooling of the variation observed with the remaining three markers resulted in two to six different mtDNA haplotypes within the different species of Picea. Evidence for intra-genomic recombination was observed in at least two taxa. Thus, these mitotypes are likely to be more informative than single-locus haplotypes. They should be particularly useful for the study of biogeography and the dynamics of hybrid zones.     

Molecular identification of cryptic bumblebee species from degraded samples using PCR–RFLP approach

The worldwide decline and local extinctions of bumblebees have raised a need for fast and accurate tools for species identification. Morphological characters are often not sufficient, and molecular methods have been increasingly used for reliable identification of bumblebee species. Molecular methods often require high‐quality DNA which makes them less suitable for analysis of low‐quality or older samples. We modified the PCR–RFLP protocol for an efficient and cost‐effective identification of four bumblebee species in the subgenus Bombus s. str. (B. lucorum, B. terrestris, B. magnus and B. cryptarum). We used a short partial mitochondrial COI fragment (446 bp) and three diagnostic restriction enzymes (Hinf I, Hinc II and Hae III) to identify species from degraded DNA material. This approach allowed us to efficiently determine the correct species from all degraded DNA samples, while only a subset of samples 64.6% (31 of 48) resulted in successful amplification of a longer COI fragment (1064 bp) using the previously described method. This protocol can be applied for conservation and management of bumblebees within this subgenus and is especially useful for fast species identification from degraded samples.     

Enterobius vermicularis: ancient DNA from North and South American human coprolites

A molecular paleoparasitological diagnostic approach was developed for Enterobius vermicularis. Ancient DNA was extracted from 27 coprolites from archaeological sites in Chile and USA. Enzymatic amplification of human mtDNA sequences confirmed the human origin. We designed primers specific to the E. vermicularis 5S ribosomal RNA spacer region and they allowed reproducible polymerase chain reaction identification of ancient material. We suggested that the paleoparasitological microscopic identification could accompany molecular diagnosis, which also opens the possibility of sequence analysis to understand parasite-host evolution.     

Phylogenetic relationships within Hevea brasiliensis as deduced from a polymorphic mitochondrial DNA region

We have cloned a 4.5-kb mtDNA fragment showing a high RFLP polymorphism between various Hevea genotypes. Subcloning and sequencing of a 1.4-kb segment of this clone allowed us to design PCR amplification primers to isolate homologous mtDNA segments of about 0.9 kb from 23 representative genotypes of Hevea. Complete sequences from 4 genotypes showed between 6.7% and 20.2% of nucleotide diversity, suggesting the presence of a hypervariable, or hotspot, region. A sequence of 345 nucleotides within this region was determined for the 23 genotypes. The phylogenetic relationships inferred from the sequence comparison are in general agreement with the results obtained from mtDNA RFLP analysis, indicating that this polymorphic mtDNA region is a useful molecular marker for phylogenetic analysis within Hevea.     

Cloning and characterization of pejerrey mitochondrial DNA and its application for RFLP analysis

Probes were cloned, characterized, and developed for all regions of the mitochondrial DNA (mtDNA) of pejerrey Odontesthes bonariensis to provide the basis for the study of genetic diversity of South American atherinopsinii and to enable species identification from small amounts of tissue. The mtDNA was extracted from liver and cleaved with Eco RI, producing four fragments (7.4, 3.4, 3.1 and 2.9 kb) which were cloned using pUC118 plasmid vectors. Sequence analysis from both ends of the fragments showed that they encode tRNA (Asp, Phe, and Ser-TGA), 12 S rRNA, cytochrome oxidase (CO) II, NADH 4, 5, and 6, and the D-loop, and that the relative positions of these genes are identical to those in the mtDNA of other teleosts. A comparison of homology with carp mtDNA nucleotide sequences revealed that tRNA (Phe and Ser-TGA) and CO II were relatively conserved, whereas the D-loop region was highly divergent. The cloned mtDNA probes detected mtDNA fragments from about 800 ng of total DNA extracted from liver, muscle, and single embryos of O. bonariensis , and were effective for restriction length fragment polymorphism (RFLP) analysis of Patagonina hatcheri , the most distant atherinopsine relative of pejerrey. The cloned mtDNA probes may be useful for the analysis of genetic diversity and non-destructive species identification, including the examination of eggs, larvae and juveniles. The mtDNA sequences reported here provide the basis for the design of primers for PCR-based RFLP analysis.