Potential efficacy of mitochondrial genes for animal DNA barcoding: a case study using eutherian mammals

Enhancing the authenticity of animal by-products: harmonization of DNA extraction methods from novel ingredients

Introduction: The increasing global pressure to explore alternative protein sources derived from animal by-products has opened-up opportunities, but it has also created the need to assess their compliance with labelling statements, to ensure consumer's trust in the composition of both feed and food products. Assessing the authenticity of highly processed animal by-products, particularly within the rapidly expanding Halal food market, presents a significant challenge due to the lack of robust and standardized methodologies. However, the success of DNA based authenticity system is highly dependent on the extracted DNA quantity, quality, and purity ratios from heterogeneous matrices. Material and methods: In this work, nine DNA extraction methods were tested on selected processed animal by-products with high-value and interest for the feed industry: meals from poultry meat, blood and feather, and hydrolysates from swine meat and bone, fish, and black soldier fly. The proposed DNA extraction methods are developed to specifically target swine-specific mitochondrial region, as a case study. Results and discussion: Both the conventional CTAB method and the commercial kits, specifically Invisorb® Spin Tissue Mini and NucleoSpin™ Food, demonstrated superior extraction efficiency and quality ratios. Nevertheless, commercial kits enabled faster detection in comparison to the conventional methods. The absence of swine DNA was successfully validated and confirmed in all animal meals and hydrolysates that did not contain swine in their composition beforehand, demonstrating their compliance with the Halal market requirements.

Disclosing the hidden nucleotide sequences: a journey into DNA barcoding of raptor species in public repositories

BACKGROUND

In nucleotide public repositories, studies discovered data errors which resulted in incorrect species identification of several accipitrid raptors considered for conservation. Mislabeling, particularly in cases of cryptic species complexes and closely related species, which were identified based on morphological characteristics, was discovered. Prioritizing accurate species labeling, morphological taxonomy, and voucher documentation is crucial to rectify spurious data.

OBJECTIVE

Our study aimed to identify an effective DNA barcoding tool that accurately reflects the efficiency status of barcodes in raptor species (Accipitridae).

METHODS

Barcode sequences, including 889 sequences from the mitochondrial cytochrome c oxidase I (COI) gene and 1052 sequences from cytochrome b (Cytb), from 150 raptor species within the Accipitridae family were analyzed.

RESULTS

The highest percentage of intraspecific nearest neighbors from the nearest neighbor test was 88.05% for COI and 95.00% for Cytb, suggesting that the Cytb gene is a more suitable marker for accurately identifying raptor species and can serve as a standard region for DNA barcoding. In both datasets, a positive barcoding gap representing the difference between inter-and intra-specific sequence divergences was observed. For COI and Cytb, the cut-off score sequence divergences for species identification were 4.00% and 3.00%, respectively.

CONCLUSION

Greater accuracy was demonstrated for the Cytb gene, making it the preferred primary DNA barcoding marker for raptors.

DNA barcoding of indigenous fowl of Manipur, Kaunayen (Gallus gallus domesticus)

The jungle fowl (Gallus gallus) is a tropical bird with important hereditary and phenotypical traits like disease resistance and resistance to harsh conditions and can often survive with scanty diet. However, as commercial chicken breeds replace them, their population is dwindling, which poses a significant threat to fowl genetic resources. There is minimal information on the variety of Indian poultry, mainly native chicken from Northeast India. As a result, the record of the fowl's genetic diversity is essential for its preservation and formulation of conservation strategies. The current study sought to identify indigenous chicken, Kaunayen (Gallus gallus domesticus), from Manipur using barcoding based on DNA sequences. A total of 5 CO1 DNA barcodes from several indigenous chickens were sequenced and compared to the previous data of diverse taxa of Phasianidae using the conventional methodology and were recognized as Gallus gallus. The Phasianid birds that were researched were accurately classified into their appropriate species. There is a minuscule genomic difference between G. gallus and G. varius (1.2%) and the highest between Arborophila rufipectus and Tympanuchus pallidicinctus (22.5%). The phylogenetic relationship established on the NJ tree revealed a coherent gathering of indigenous fowl with G. gallus and unique to all other species studied, showing their taxonomic classification. Nonetheless, the investigation offered a genetic identity tag for indigenous chicken for the first time. It will be a potential guide for identifying distinctive and genetically unique poultry sequences for later requirements.

Comparative Mitogenomic Analysis of Two Snake Eels Reveals Irregular Gene Rearrangement and Phylogenetic Implications of Ophichthidae

The family Ophichthidae has the largest number and the most various species (about 359 valid species) in the order Anguilliformes worldwide. Both morphological and molecular characteristics have been used to assess their taxonomic status. However, due to the ambiguous morphological features, molecular data such as mitochondrial DNA sequences have been implemented for the correct identification and classification of these fishes. In this study, the gene arrangement and structure characteristics of two Ophichthidae mitochondrial genomes were investigated for the first time. The total mitogenome lengths of O. evermanni and O. erabo were 17,759 bp and 17,856 bp, respectively. Comparing with the ancestral mitochondrial gene order, the irregular gene rearrangement happened between ND6 and tRNA-Pro (P) genes with another similar control region emerging between tRNA-Thr (T) and ND6 genes, which could be explained by the tandem duplication and random loss (TDRL) model appropriately. ML phylogenetic tree demonstrated that the family Ophichthidae was monophyletic origin, but genus Ophichthus might be polyphyletic because of the confused cluster relationships among different species.

Mitochondrial genome diversity of Balamuthia mandrillaris revealed by a fatal case of granulomatous amoebic encephalitis

Introduction

Balamuthia (B.) mandrillaris is a free-living amoeba that can cause rare yet fatal granulomatous amoebic encephalitis (GAE). However, efficacious treatment for GAE is currently unavailable, especially when genomic studies on B. mandrillaris are limited.

Methods

In this study, B. mandrillaris strain KM-20 was isolated from the brain tissue of a GAE patient, and its mitochondrial genome was de novo assembled using high-coverage Nanopore long reads and Illumina short reads.

Results and Discussion

Phylogenetic and comparative analyses revealed a range of diversification in the mitochondrial genome of KM-20 and nine other B. mandrillaris strains. According to the mitochondrial genome alignment, one of the most variable regions was observed in the ribosomal protein S3 (rps3), which was caused by an array of novel protein tandem repeats. The repeating units in the rps3 protein tandem region present significant copy number variations (CNVs) among B. mandrillaris strains and suggest KM-20 as the most divergent strain for its highly variable sequence and highest copy number in rps3. Moreover, mitochondrial heteroplasmy was observed in strain V039, and two genotypes of rps3 are caused by the CNVs in the tandem repeats. Taken together, the copy number and sequence variations of the protein tandem repeats enable rps3 to be a perfect target for clinical genotyping assay for B. mandrillaris. The mitochondrial genome diversity of B. mandrillaris paves the way to investigate the phylogeny and diversification of pathogenic amoebae.

Survey of mitochondrial sequences integrated into the bovine nuclear genome

Nuclear copies of the mitochondrial DNA (NUMTs) have already been described in several species. In this context, we identified and analysed 166 bovine NUMT regions with a total length of 430 kbp, representing about 0.02% of the cattle nuclear genome. Copies of all mitochondrial regions were detected in the nuclear genome, with distinct degrees of sequence similarity to the mitogenome. Some NUMT regions include large mitogenome segments and show high similarity to the organelle genome sequence. NUMT regions are frequently modified by insertion of repetitive sequences and by sequence rearrangements. We confirmed the existence of 29 NUMT regions by PCR amplification using DNA from the cow (Dominette) which was used to generate the bovine genome reference sequence, ruling out the possibility that these NUMTs could be artifacts of the genome assembly. As there are NUMT regions with high similarity to the mitogenome, special care is needed when designing primers for mitochondrial DNA amplification. Our results can therefore be used to avoid co-amplification of bovine nuclear sequences similar to mitochondrial DNA.

Mitogenomic phylogenetic analyses of Leptogorgia virgulata and Leptogorgia hebes (Anthozoa: Octocorallia) from the Gulf of Mexico provides insight on Gorgoniidae divergence between Pacific and Atlantic lineages

The use of genetics in recent years has brought to light the need to reevaluate the classification of many gorgonian octocorals. This study focuses on two Leptogorgia species-Leptogorgia virgulata and Leptogorgia hebes-from the northwestern Gulf of Mexico (GOM). We target complete mitochondrial genomes and mtMutS sequences, and integrate this data with previous genetic research of gorgonian corals to resolve phylogenetic relationships and estimate divergence times. This study contributes the first complete mitochondrial genomes for L. ptogorgia virgulata and L. hebes. Our resulting phylogenies stress the need to redefine the taxonomy of the genus Leptogorgia in its entirety. The fossil-calibrated divergence times for Eastern Pacific and Western Atlantic Leptogorgia species based on complete mitochondrial genomes shows that the use of multiple genes results in estimates of more recent speciation events than previous research based on single genes. These more recent divergence times are in agreement with geologic data pertaining to the formation of the Isthmus of Panama.

Portunus pelagicus mtDNA heteroplasmy inheritance and its effect on the use of mtCR and mtCOI sequence data

Mitochondrial DNA (mtDNA) sequences, especially mitochondrial control region (mtCR) and mitochondrial cytochrome c oxidase subunit I (mtCOI), have been widely used in population and evolutionary genetic analyses of metazoan. The presence of mtDNA heteroplasmy - a mixture of mtDNA haplotypes - possibly affects these analyses. This study aimed to reveal mtDNA heteroplasmy in mtCR, mtCOI, and mtND2 (mitochondrial NADH dehydrogenase subunit 2) of Portunus pelagicus, and examine its effect on the use of mtCR and mtCOI sequences. The screening result showed that the probability of observing mtDNA heteroplasmy was approximately 8%. Across the three targeted regions, 92 heteroplasmic variants were observed from seven samples comprising three mothers and four offspring. Most inherited heteroplasmy presented transition and silence mutation. By comparing the proportion of shared variants among maternal relatives to that among non-relatives, the result suggested that most heteroplasmic variants observed in an individual are inherited. Statistical analyses carried out on the inter-generational differences suggested that random drift and purifying selection play roles in determining the offspring's heteroplasmy level. The size of the random shift varies according to the location of variants and the mothers. The phylogenetic analysis showed that the presence of mtDNA heteroplasmy in mtCR and mtCOI does not affect familial and species identification, respectively. This study firstly reported the mtDNA heteroplasmy in P. pelagicus, its inheritance pattern, and its effect on the use of mtDNA sequence data. This basic knowledge would be useful for the study based on mtDNA sequence data, especially in other invertebrates.

What does mitogenomics tell us about the evolutionary history of the Drosophila buzzatii cluster (repleta group)?

The Drosophila repleta group is an array of more than 100 species endemic to the “New World”, many of which are cactophilic. The ability to utilize decaying cactus tissues as breeding and feeding sites is a key aspect that allowed the successful diversification of the repleta group in American deserts and arid lands. Within this group, the Drosophila buzzatii cluster is a South American clade of seven closely related species in different stages of divergence, making them a valuable model system for evolutionary research. Substantial effort has been devoted to elucidating the phylogenetic relationships among members of the D. buzzatii cluster, including molecular phylogenetic studies that have generated ambiguous results where different tree topologies have resulted dependent on the kinds of molecular marker used. Even though mitochondrial DNA regions have become useful markers in evolutionary biology and population genetics, none of the more than twenty Drosophila mitogenomes assembled so far includes this cluster. Here, we report the assembly of six complete mitogenomes of five species: D. antonietae, D. borborema, D. buzzatii, two strains of D. koepferae and D. seriema, with the aim of revisiting phylogenetic relationships and divergence times by means of mitogenomic analyses. Our recovered topology using complete mitogenomes supports the hypothesis of monophyly of the D. buzzatii cluster and shows two main clades, one including D. buzzatii and D. koepferae (both strains), and the other containing the remaining species. These results are in agreement with previous reports based on a few mitochondrial and/or nuclear genes, but conflict with the results of a recent large-scale nuclear phylogeny, indicating that nuclear and mitochondrial genomes depict different evolutionary histories.

Spatial and temporal dynamics of Antarctic shallow soft-bottom benthic communities: ecological drivers under climate change

Background

Marine soft sediments are some of the most widespread habitats in the ocean, playing a vital role in global carbon cycling, but are amongst the least studied with regard to species composition and ecosystem functioning. This is particularly true of the Polar Regions, which are currently undergoing rapid climate change, the impacts of which are poorly understood. Compared to other latitudes, Polar sediment habitats also experience additional environmental drivers of strong seasonality and intense disturbance from iceberg scouring, which are major structural forces for hard substratum communities. This study compared sediment assemblages from two coves, near Rothera Point, Antarctic Peninsula, 67°S in order to understand the principal drivers of community structure, for the first time, evaluating composition across all size classes from mega- to micro-fauna.

Results

Morpho-taxonomy identified 77 macrofaunal species with densities of 464–16,084 individuals m⁻². eDNA metabarcoding of microfauna, in summer only, identified a higher diversity, 189 metazoan amplicon sequence variants (ASVs) using the 18S ribosomal RNA and 249 metazoan ASVs using the mitochondrial COI gene. Both techniques recorded a greater taxonomic diversity in South Cove than Hangar Cove, with differences in communities between the coves, although the main taxonomic drivers varied between techniques. Morphotaxonomy identified the main differences between coves as the mollusc, Altenaeum charcoti, the cnidarian Edwardsia sp. and the polychaetes from the family cirratulidae. Metabarcoding identified greater numbers of species of nematodes, crustaceans and Platyhelminthes in South Cove, but more bivalve species in Hangar Cove. There were no detectable differences in community composition, measured through morphotaxonomy, between seasons, years or due to iceberg disturbance.

Conclusions

This study found that unlike hard substratum communities the diversity of Antarctic soft sediment communities is correlated with the same factors as other latitudes. Diversity was significantly correlated with grain size and organic content, not iceberg scour. The increase in glacial sediment input as glaciers melt, may therefore be more important than increased iceberg disturbance.

Electronic supplementary material

The online version of this article (10.1186/s12898-019-0244-x) contains supplementary material, which is available to authorized users.

Reconciling the Mitonuclear Compatibility Species Concept with Rampant Mitochondrial Introgression

The mitonuclear compatibility species concept defines a species as a population that is genetically isolated from other populations by uniquely coadapted mitochondrial (mt) and nuclear genes. A key prediction of this hypothesis is that the mt genotype of each species will be functionally distinct and that introgression of mt genomes will be prevented by mitonuclear incompatibilities that arise when heterospecific mt and nuclear genes attempt to cofunction to enable aerobic respiration. It has been proposed, therefore, that the observation of rampant introgression of mt genotypes from one species to another constitutes a strong refutation of the mitonuclear speciation. The displacement of a mt genotype from a nuclear background with which it co-evolved to a foreign nuclear background will necessarily lead to fitness loss due to mitonuclear incompatibilities. Here I consider two potential benefits of mt introgression between species that may, in some cases, overcome fitness losses arising from mitonuclear incompatibilities. First, the introgressed mt genotype may be better adapted to the local environment than the native mt genotype such that higher fitness is achieved through improved adaptation via introgression. Second, if the mitochondria of the recipient taxa carry a high mutational load, then introgression of a foreign, less corrupt mt genome may enable the recipient taxa to escape its mutational load and gain a fitness advantage. Under both scenarios, fitness gains from novel mt genotypes could theoretically compensate for the fitness that is lost via mitonuclear incompatibility. I also consider the role of endosymbionts in non-adaptive rampant introgression of mt genomes. I conclude that rampant introgression is not necessarily evidence against the idea of tight mitonuclear coadaptation or the mitonuclear compatibility species concept. Rampant mt introgression will typically lead to erasure of species but in some cases could lead to hybrid speciation.

Genetic status of indigenous poultry (red jungle fowl) from India

The global biodiversity of domesticated red jungle fowl (Gallus gallus) is gradually eroding by replacement with commercial poultry breeds and results loss of valuable genetic and physical traits like resistance to disease, extreme environment, etc. posing a threat to the poultry genetic resources. Very fewer reports exist on Indian poultry diversity, especially native chicken of India. Therefore, species identification and inventorying of the poultry genetic resource is indispensable. Thus, the present study aimed to characterize indigenous chicken from bio-diversity hotspot of Sunderban and Northeast India using DNA sequence based barcoding approach. A total of 15 CO1 (Cytochrome c Oxidase-I) DNA barcode of different indigenous chicken were newly sequenced along with 6 previously published sequences from our laboratory and compared with the available data of distinctive genera of Phasianidae as per the standard protocol and are identified as Gallus gallus. About 98.96% of the Phasianid birds were successfully delimitated into the respective species except for 12 congeneric pairs whose minimum interspecific K2P (Kimura 2-parameter) distance overlaps with the maximum intraspecific distance (3.9%). The least genetic divergence is observed between G. gallus and G. varius (0.013%) and highest between G. gallus and G. lafayettei (0.059%). The NJ tree showed a cohesive clustering of indigenous chicken with G. gallus and distinct with respect to all the different species under study, thereby revealing their taxonomic position except for few G. sonneratti that showed mixed clustering with G. gallus. This may be due to the genetic introgression between the species. Nevertheless, the study for the first time provided the molecular identification tag of indigenous poultry from biodiversity hotspot of East and Northeast India and will remain as a potential guide to recognize inimitable and valuable poultry genetic resources for future needs.

Improved method for genotyping the causative agent of crayfish plague (Aphanomyces astaci) based on mitochondrial DNA

Aphanomyces astaci causes crayfish plague, which is a devastating disease of European freshwater crayfish. The likely first introduction of A. astaci into Europe was in the mid-19th century in Italy, presumably with the introduction of North American crayfish. These crayfish can carry A. astaci in their cuticle as a benign infection. Aphanomyces astaci rapidly spread across Europe causing the decline of the highly susceptible indigenous crayfish species. Random amplified polymorphic DNA-PCR analysis of A. astaci pure cultures characterized five genotype groups (A, B, C, D and E). Current A. astaci genotyping techniques (microsatellites and genotype-specific regions, both targeting nuclear DNA) can be applied directly to DNA extracted from infected cuticles but require high infection levels. Therefore, they are not suitable for genotyping benign infections in North American crayfish (carriers). In the present study, we combine bioinformatics and molecular biology techniques to develop A. astaci genotyping molecular markers that target the mitochondrial DNA, increasing the sensitivity of the genotyping tools. The assays were validated on DNA extracts of A. astaci pure cultures, crayfish tissue extractions from crayfish plague outbreaks and tissue extractions from North American carriers. We demonstrate the presence of A. astaci genotype groups A and B in UK waters.

An Extension of the Kimura Two-Parameter Model to the Natural Evolutionary Process

Accurate estimates of genetic difference are required for research in evolutionary biology. Here we extend the Kimura two-parameter (K2P) model by considering gaps (insertions and/or deletions) and introduce a new measure for estimating genetic difference between two nucleotide sequences in terms of nucleotide changes that have occurred during the evolutionary process. Using the nuclear ribosomal DNA internal transcribed spacer 2 region from the genus Physalis, we demonstrate that species identification and phylogenetic studies strongly depend on evolutionary models. It is especially noteworthy that the use of different models affects the degree of overlap between intraspecific and interspecific genetic differences. We observe that the percentage of interspecific sequence pairs with values less than the maximum intraspecific genetic difference is 43.2% for the K2P model which is calculated by removing gap sites across all sequences, 22.7% for the K2P model which is calculated by removing gap sites for sequence pairs, and 16.9% for our model which is calculated without removing gap sites. Additionally, the numbers of sequence pairs with interspecific genetic differences of zero are 50 for the K2P model and 29 for our model. The genetic difference measure based on the K2P model, compared to our model, overestimates 21 sequence pairs that are not originally identical. These results indicate the importance of estimating genetic differences under the model of sequence evolution that includes insertions and deletions in addition to substitutions.

A new efficient method for analyzing fungi species using correlations between nucleotides

BACKGROUND

In recent years, DNA barcoding has become an important tool for biologists to identify species and understand their natural biodiversity. The complexity of barcode data makes it difficult to analyze quickly and effectively. Manual classification of this data cannot keep up to the rate of increase of available data.

RESULTS

In this study, we propose a new method for DNA barcode classification based on the distribution of nucleotides within the sequence. By adding the covariance of nucleotides to the original natural vector, this augmented 18-dimensional natural vector makes good use of the available information in the DNA sequence. The accurate classification results we obtained demonstrate that this new 18-dimensional natural vector method, together with the random forest classifier algorthm, can serve as a computationally efficient identification tool for DNA barcodes. We performed phylogenetic analysis on the genus Megacollybia to validate our method. We also studied how effective our method was in determining the genetic distance within and between species in our barcoding dataset.

CONCLUSIONS

The classification performs well on the fungi barcode dataset with high and robust accuracy. The reasonable phylogenetic trees we obtained further validate our methods. This method is alignment-free and does not depend on any model assumption, and it will become a powerful tool for classification and evolutionary analysis.

Is molecular evolution faster in the tropics?

The evolutionary speed hypothesis (ESH) suggests that molecular evolutionary rates are higher among species inhabiting warmer environments. Previously, the ESH has been investigated using small numbers of latitudinally-separated sister lineages; in animals, these studies typically focused on subsets of Chordata and yielded mixed support for the ESH. This study analyzed public DNA barcode sequences from the cytochrome c oxidase subunit I (COI) gene for six of the largest animal phyla (Arthropoda, Chordata, Mollusca, Annelida, Echinodermata, and Cnidaria) and paired latitudinally-separated taxa together informatically. Of 8037 lineage pairs, just over half (51.6%) displayed a higher molecular rate in the lineage inhabiting latitudes closer to the equator, while the remainder (48.4%) displayed a higher rate in the higher-latitude lineage. To date, this study represents the most comprehensive analysis of latitude-related molecular rate differences across animals. While a statistically-significant pattern was detected from our large sample size, our findings suggest that the EHS may not serve as a strong universal mechanism underlying the latitudinal diversity gradient and that COI molecular clocks may generally be applied across latitudes. This study also highlights the merits of using automation to analyze large DNA barcode datasets.

Over 2.5 million COI sequences in GenBank and growing

The increasing popularity of cytochrome c oxidase subunit 1 (COI) DNA metabarcoding warrants a careful look at the underlying reference databases used to make high-throughput taxonomic assignments. The objectives of this study are to document trends and assess the future usability of COI records for metabarcode identification. The number of COI records deposited to the NCBI nucleotide database has increased by a geometric average of 51% per year, from 8,137 records deposited in 2003 to a cumulative total of ~ 2.5 million by the end of 2017. About half of these records are fully identified to the species rank, 92% are at least 500 bp in length, 74% have a country annotation, and 51% have latitude-longitude annotations. To ensure the future usability of COI records in GenBank we suggest: 1) Improving the geographic representation of COI records, 2) Improving the cross-referencing of COI records in the Barcode of Life Data System and GenBank to facilitate consolidation and incorporation into existing bioinformatic pipelines, 3) Adherence to the minimum information about a marker gene sequence guidelines, and 4) Integrating metabarcodes from eDNA and mixed community studies with existing reference sequences. The growth of COI reference records over the past 15 years has been substantial and is likely to be a resource across many fields for years to come.

Gelatin controversies in food, pharmaceuticals, and personal care products: Authentication methods, current status, and future challenges

Gelatin is a highly purified animal protein of pig, cow, and fish origins and is extensively used in food, pharmaceuticals, and personal care products. However, the acceptability of gelatin products greatly depends on the animal sources of the gelatin. Porcine and bovine gelatins have attractive features but limited acceptance because of religious prohibitions and potential zoonotic threats, whereas fish gelatin is welcomed in all religions and cultures. Thus, source authentication is a must for gelatin products but it is greatly challenging due to the breakdown of both protein and DNA biomarkers in processed gelatins. Therefore, several methods have been proposed for gelatin identification, but a comprehensive and systematic document that includes all of the techniques does not exist. This up-to-date review addresses this research gap and presents, in an accessible format, the major gelatin source authentication techniques, which are primarily nucleic acid and protein based. Instead of presenting these methods in paragraph form which needs much attention in reading, the major methods are schematically depicted, and their comparative features are tabulated. Future technologies are forecasted, and challenges are outlined. Overall, this review paper has the merit to serve as a reference guide for the production and application of gelatin in academia and industry and will act as a platform for the development of improved methods for gelatin authentication.

Utility of mtDNA-COI Barcode Region for Phylogenetic Relationship and Diagnosis of Five Common Pest Cockroaches

BACKGROUND

Cockroaches are of vital importance medically and hygienically as they can disperse human pathogenic agents and are especially responsible for food contamination and spreading of food borne pathogens. In this study, part of mtDNA-COI gene of five common pest cockroaches was tested for diagnostic and phylogenetic purposes.

METHODS

We have described barcode region of mtDNA-COI gene of five cockroach species: Blattella germanica, Blatta orientalis, Periplaneta americana, Shelfordella lateralis, and Supella longipalpa, along with the development of a PCR-RFLP method for rapid detection and differentiation of these health pest species.

RESULTS

The PCR generates a single 710 bp-sized amplicon in all cockroach specimens, followed by direct sequencing. AluI predicted from the sequencing data provided different RFLP profiles among five species. There was a significant intra-species variation within the American cockroach populations, but no genetic variation within other species. Accordingly, phylogenetic analysis demonstrates common monophyly for cockroach families in agreement with conventional taxonomy. However S. longipalpa (Ectobiidae) diverged as an early ancestor of other cockroaches and was not associated with other Ectobiidae.

CONCLUSION

The PCR-RFLP protocol might be useful when the conventional taxonomic methods are not able to identify specimens, particularly when only small body parts of specimens are available or they are in a decaying condition. mtDNA-COI gene shows potentially useful for studying phylogenetic relationships of Blattodea order.

Advances in DNA metabarcoding for food and wildlife forensic species identification

Species identification using DNA barcodes has been widely adopted by forensic scientists as an effective molecular tool for tracking adulterations in food and for analysing samples from alleged wildlife crime incidents. DNA barcoding is an approach that involves sequencing of short DNA sequences from standardized regions and comparison to a reference database as a molecular diagnostic tool in species identification. In recent years, remarkable progress has been made towards developing DNA metabarcoding strategies, which involves next-generation sequencing of DNA barcodes for the simultaneous detection of multiple species in complex samples. Metabarcoding strategies can be used in processed materials containing highly degraded DNA e.g. for the identification of endangered and hazardous species in traditional medicine. This review aims to provide insight into advances of plant and animal DNA barcoding and highlights current practices and recent developments for DNA metabarcoding of food and wildlife forensic samples from a practical point of view. Special emphasis is placed on new developments for identifying species listed in the Convention on International Trade of Endangered Species (CITES) appendices for which reliable methods for species identification may signal and/or prevent illegal trade. Current technological developments and challenges of DNA metabarcoding for forensic scientists will be assessed in the light of stakeholders' needs.

The establishment of species-specific primers for the molecular identification of ten stored-product psocids based on ITS2 rDNA

Psocids are important stored product pests found worldwide that can be spread through grain trade. Most stored-product psocids, including eggs, nymphs, and adults, are very small (~1 mm) and difficult to identify morphologically. Here, we collected 10 economically important stored-product Liposcelis spp. psocids (L. bostrychophila, L. entomophila, L. decolor, L. paeta, L. brunnea, L. corrodens, L. mendax, L. rufa, L. pearmani, and L. tricolor) from 35 geographical locations in 5 countries (China, Czech Republic, Denmark, Germany, and the United States). The ITS2 rDNA gene was extracted and sequenced. The interspecific genetic distance of the stored-product psocids was significantly higher than the intraspecific genetic distance according to the barcoding gap analysis. Ten pairs of species-specific primers based on the ITS2 rDNA were developed for psocid identification. The sensitivity estimation indicated that the species-specific primers could correctly amplify the target ITS2 gene and successfully identify psocids at 1.0 ng/mL. Additionally, these species-specific primers could quantify specificity and identify 10 stored-product psocids; this approach could also be used to accurately identify other stored-product psocids. This work provides a practical approach for the precise examination of 10 stored-product psocid species and also contributes to the development of an identification method using ITS2 rDNA.

A simulation study of sample size for DNA barcoding

For some groups of organisms, DNA barcoding can provide a useful tool in taxonomy, evolutionary biology, and biodiversity assessment. However, the efficacy of DNA barcoding depends on the degree of sampling per species, because a large enough sample size is needed to provide a reliable estimate of genetic polymorphism and for delimiting species. We used a simulation approach to examine the effects of sample size on four estimators of genetic polymorphism related to DNA barcoding: mismatch distribution, nucleotide diversity, the number of haplotypes, and maximum pairwise distance. Our results showed that mismatch distributions derived from subsamples of ≥20 individuals usually bore a close resemblance to that of the full dataset. Estimates of nucleotide diversity from subsamples of ≥20 individuals tended to be bell‐shaped around that of the full dataset, whereas estimates from smaller subsamples were not. As expected, greater sampling generally led to an increase in the number of haplotypes. We also found that subsamples of ≥20 individuals allowed a good estimate of the maximum pairwise distance of the full dataset, while smaller ones were associated with a high probability of underestimation. Overall, our study confirms the expectation that larger samples are beneficial for the efficacy of DNA barcoding and suggests that a minimum sample size of 20 individuals is needed in practice for each population.

Genetic diversity in cytochrome c oxidase I gene of Anopheles mosquitoes

Genetic diversity in cytochrome c oxidase I (coxI) among 7 species of Anopheles mosquitoes from Pakistan, and 37 species from different geographical regions of the world, was recorded. Automatic Barcode Gap Discovery (ABGD) analysis revealed a clear gap between intraspecific and interspecific distances of 7 species from Pakistan. However, genetic distances of 37 Anopheles species failed to adequately differentiate species in a global context. Intraspecific and interspecific divergences for 7 Anopheles species of Pakistan varied from 0.0% to 2.5% (mean = 0.49%) and 8% to 22.3% (mean = 12.77%), respectively. Similarly, intraspecific distances for 37 species from different parts of world ranged from 0.0% to 11.2% (mean = 0.65%) while values of interspecific divergences ranged from 3.4% to 35% (mean = 11.75%). Although phylogenetic tree revealed separate clades for 7 Anopheles species of Pakistan, it failed to produce separate clades for 37 species of the world. It is concluded that although standard barcode region is helpful for identifying Anopheles mosquitoes, combination of multi-locus approaches and morphology may be required to accurately identify species in this genus.

New tools (DNA barcoding), old hypothesis: the case of the taxonomic identity of the Argentine hakes (Actinopterygii: Merluccius)

The present study evaluated the possible occurrence of cryptic species among Merluccidae from Argentina by examining sequences of cytochrome c oxidase subunit I (coI) mtDNA. This approach can discriminate Merluccius hubbsi and Merluccius australis; specimens with morphological diagnostic characters of Merluccius patagonicus formed a cohesive cluster with M. hubbsi specimens. BIN analysis confirmed the effectiveness of barcoding within a global context.

Taxonomic resolutions based on 18S rRNA genes: a case study of subclass copepoda

Biodiversity studies are commonly conducted using 18S rRNA genes. In this study, we compared the inter-species divergence of variable regions (V1-9) within the copepod 18S rRNA gene, and tested their taxonomic resolutions at different taxonomic levels. Our results indicate that the 18S rRNA gene is a good molecular marker for the study of copepod biodiversity, and our conclusions are as follows: 1) 18S rRNA genes are highly conserved intra-species (intra-species similarities are close to 100%); and could aid in species-level analyses, but with some limitations; 2) nearly-whole-length sequences and some partial regions (around V2, V4, and V9) of the 18S rRNA gene can be used to discriminate between samples at both the family and order levels (with a success rate of about 80%); 3) compared with other regions, V9 has a higher resolution at the genus level (with an identification success rate of about 80%); and 4) V7 is most divergent in length, and would be a good candidate marker for the phylogenetic study of Acartia species. This study also evaluated the correlation between similarity thresholds and the accuracy of using nuclear 18S rRNA genes for the classification of organisms in the subclass Copepoda. We suggest that sample identification accuracy should be considered when a molecular sequence divergence threshold is used for taxonomic identification, and that the lowest similarity threshold should be determined based on a pre-designated level of acceptable accuracy.

Identification of ungulates used in a traditional Chinese medicine with DNA barcoding technology

Horns of Saiga antelope (Saiga tatarica) have always been an ingredient of “Lingyangjiao”, a traditional Chinese medicine (TCM). Persistent hunting for Saiga antelope has already threatened the survival of critical endangered populations in wild. To control the growing pressure, CITES and Chinese government have legislated for monitoring the trade of Saiga horns. However, similar ungulate horns are difficult to identify by their morphological characteristics, which has impeded the law enforcement. Besides Saiga antelope, other seven ungulate species which have similar horns are also sold and marked as “Lingyangjiao” in TCM markets to offset shortage of Saiga antelope horns. Such species are Gazella subgutturosa, Pantholops hodgsonii, Procapra picticaudata, Procapra gutturosa, Procapra przewalskii, Capra hircus, and Ovis aries. Our study aimed at implementing DNA barcoding technology to diagnose Saiga horns and the substitutes. We successfully extracted genomic DNA from horn samples. We recovered COI sequences of 644 bp with specific primers and 349 bp with nested PCR primers designed for degraded horn samples. The mean interspecific genetic distance of data set of the 644-bp full barcodes and the 349-bp mini-barcodes was 14.96% and 15.38%, respectively, and the mean intraspecific distance was 0.24% and 0.20%, respectively. Each species formed independent clades in neighbor-joining (NJ) phylogenetic tree of the two data sets with >99% supporting values, except P. gutturosa and P. przewalskii. The deep genetic distances gap and clear species clades in NJ tree of either full barcodes or mini-barcodes suggest that barcoding technology is an effective tool to diagnose Saiga horns and their substitutes. Barcoding diagnosis protocol developed here will simplify diagnosis of “Lingyangjiao” species and will facilitate conservation of endangered ungulates involved in TCM “Lingyangjiao” markets, especially the Saiga antelope.

Genetic identification of endangered North African ungulates using noninvasive sampling

North African ungulates include several threatened and emblematic species, yet are poorly studied mainly due to their remoteness and elusiveness. Noninvasive sampling provides a useful approach to obtain ecological and genetic information essential to guide conservation actions. The very first and most important step in conservation planning is to accurately identify species, and molecular genetics has been proved to be a useful tool. Several molecular genetics protocols are available for species identification, even for samples with poor quality DNA, such as faeces, hairs or bones. Most of these protocols use mitochondrial DNA for barcoding despite this marker being especially prone to problems, including mtDNA introgression, nuclear insert copies, high intraspecific diversity or heteroplasmy. In this work, we developed a molecular method based on polymorphisms in small fragments of the mitochondrial cytochrome b (cytb, mtDNA) and the nuclear kappa casein genes (KCAS, nDNA) for identifying endangered North African ungulates. These fragments revealed polymorphisms, including species-specific variation, which allowed species identification of nine ungulate species that co-occur in North Africa. The method was validated across more than 400 samples, including different types of noninvasive samples collected in the field. The simplicity, high reliability and relative low cost of the described method make it a promising tool to improve ecological studies of the North African ungulates and consequently, the implementation of more efficient management and conservation plans for these endangered ungulates.

Mitochondrial data are not suitable for resolving placental mammal phylogeny

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

Evaluation of the efficacy of twelve mitochondrial protein-coding genes as barcodes for mollusk DNA barcoding

In this study, we evaluated the efficacy of 12 mitochondrial protein-coding genes from 238 mitochondrial genomes of 140 molluscan species as potential DNA barcodes for mollusks. Three barcoding methods (distance, monophyly and character-based methods) were used in species identification. The species recovery rates based on genetic distances for the 12 genes ranged from 70.83 to 83.33%. There were no significant differences in intra- or interspecific variability among the 12 genes. The monophyly and character-based methods provided higher resolution than the distance-based method in species delimitation. Especially in closely related taxa, the character-based method showed some advantages. The results suggested that besides the standard COI barcode, other 11 mitochondrial protein-coding genes could also be potentially used as a molecular diagnostic for molluscan species discrimination. Our results also showed that the combination of mitochondrial genes did not enhance the efficacy for species identification and a single mitochondrial gene would be fully competent.

Using environmental DNA to census marine fishes in a large mesocosm

The ocean is a soup of its resident species' genetic material, cast off in the forms of metabolic waste, shed skin cells, or damaged tissue. Sampling this environmental DNA (eDNA) is a potentially powerful means of assessing whole biological communities, a significant advance over the manual methods of environmental sampling that have historically dominated marine ecology and related fields. Here, we estimate the vertebrate fauna in a 4.5-million-liter mesocosm aquarium tank at the Monterey Bay Aquarium of known species composition by sequencing the eDNA from its constituent seawater. We find that it is generally possible to detect mitochondrial DNA of bony fishes sufficient to identify organisms to taxonomic family- or genus-level using a 106 bp fragment of the 12S ribosomal gene. Within bony fishes, we observe a low false-negative detection rate, although we did not detect the cartilaginous fishes or sea turtles present with this fragment. We find that the rank abundance of recovered eDNA sequences correlates with the abundance of corresponding species' biomass in the mesocosm, but the data in hand do not allow us to develop a quantitative relationship between biomass and eDNA abundance. Finally, we find a low false-positive rate for detection of exogenous eDNA, and we were able to diagnose non-native species' tissue in the food used to maintain the mesocosm, underscoring the sensitivity of eDNA as a technique for community-level ecological surveys. We conclude that eDNA has substantial potential to become a core tool for environmental monitoring, but that a variety of challenges remain before reliable quantitative assessments of ecological communities in the field become possible.

DNA barcoding of Dutch birds

The mitochondrial cytochrome c oxidase subunit I (COI) can serve as a fast and accurate marker for the identification of animal species, and has been applied in a number of studies on birds. We here sequenced the COI gene for 387 individuals of 147 species of birds from the Netherlands, with 83 species being represented by > 2 sequences. The Netherlands occupies a small geographic area and 95% of all samples were collected within a 50 km radius from one another. The intraspecific divergences averaged 0.29% among this assemblage, but most values were lower; the interspecific divergences averaged 9.54%. In all, 95% of species were represented by a unique barcode, with 6 species of gulls and skua (Larus and Stercorarius) having at least one shared barcode. This is best explained by these species representing recent radiations with ongoing hybridization. In contrast, one species, the Lesser Whitethroat Sylvia curruca showed deep divergences, averaging 5.76% and up to 8.68% between individuals. These possibly represent two distinct taxa, S. curruca and S. blythi, both clearly separated in a haplotype network analysis. Our study adds to a growing body of DNA barcodes that have become available for birds, and shows that a DNA barcoding approach enables to identify known Dutch bird species with a very high resolution. In addition some species were flagged up for further detailed taxonomic investigation, illustrating that even in ornithologically well-known areas such as the Netherlands, more is to be learned about the birds that are present.

DNA barcoding and species boundary delimitation of selected species of Chinese Acridoidea (Orthoptera: Caelifera)

We tested the performance of DNA barcoding in Acridoidea and attempted to solve species boundary delimitation problems in selected groups using COI barcodes. Three analysis methods were applied to reconstruct the phylogeny. K2P distances were used to assess the overlap range between intraspecific variation and interspecific divergence. "Best match (BM)", "best close match (BCM)", "all species barcodes (ASB)" and "back-propagation neural networks (BP-based method)" were utilized to test the success rate of species identification. Phylogenetic species concept and network analysis were employed to delimitate the species boundary in eight selected species groups. The results demonstrated that the COI barcode region performed better in phylogenetic reconstruction at genus and species levels than at higher-levels, but showed a little improvement in resolving the higher-level relationships when the third base data or both first and third base data were excluded. Most overlaps and incorrect identifications may be due to imperfect taxonomy, indicating the critical role of taxonomic revision in DNA barcoding study. Species boundary delimitation confirmed the presence of oversplitting in six species groups and suggested that each group should be treated as a single species.

A mitochondrial phylogeny and biogeographical scenario for Asiatic water shrews of the genus Chimarrogale: implications for taxonomy and low-latitude migration routes

The six species and three subspecies in the genus Chimarrogale (Soricomorpha: Soricidae) are commonly referred to as Asiatic water shrews. The Chimarrogale are the most widely distributed group of Nectogaline shrews, extending throughout the Oriental region and Japan. Because of the limited numbers of specimens available for study, the phylogenetic relationships and biogeographical history of this genus have not been comprehensively discussed. We used mitochondrial cytochrome b gene sequences to estimate phylogenetic relationships and divergence times among four Chimarrogale species, including all three subspecies of Chimarrogale himalayica. We also conducted a species delimitation analysis and tested two alternative migration scenarios in Asia through species distribution modeling and a reconstruction of the ancestral distribution. Here, we present the first proposed hypothesis regarding the Asiatic water shrew phylogeny and reveal ten putative species within the four recognized species. Distinct phylogenetic statuses of Chimarrogale phaeura, Chimarrogale platycephala, and Chimarrogale styani were confirmed. Chimarrogale himalayica was strongly supported as paraphyletic. We suggest that three subspecies of Chimarrogale himalayica should be reconsidered as distinct species. However, these suggestions must be considered with caution because only a single locus of a mtDNA gene was used. Four additional putative species, possibly distributed in central southwestern China and Taiwan, are currently undescribed; therefore, comprehensive morphological analyses are warranted to test their taxonomic statuses. The estimated molecular divergence times indicated that rapid speciation occurred during the early Pliocene, and current distribution patterns may have been affected by global cooling during the Pliocene/Pleistocene boundary. Reconstruction of the ancestral distribution and species distribution modeling for Asiatic water shrews revealed a low-latitude migration route over which ancestral Chimarrogale migrated from Europe via Central Asia to their current distribution. Our results demonstrated that Asiatic water shrews could have evolved throughout the low-latitude migration route from Europe to East and Southeast Asia.

Description of new mitochondrial genomes (Spodoptera litura, Noctuoidea and Cnaphalocrocis medinalis, Pyraloidea) and phylogenetic reconstruction of Lepidoptera with the comment on optimization schemes

We newly sequenced mitochondrial genomes of Spodoptera litura and Cnaphalocrocis medinalis belonging to Lepidoptera to obtain further insight into mitochondrial genome evolution in this group and investigated the influence of optimal strategies on phylogenetic reconstruction of Lepidoptera. Estimation of p-distances of each mitochondrial gene for available taxonomic levels has shown the highest value in ND6, whereas the lowest values in COI and COII at the nucleotide level, suggesting different utility of each gene for different hierarchical group when individual genes are utilized for phylogenetic analysis. Phylogenetic analyses mainly yielded the relationships (((((Bombycoidea + Geometroidea) + Noctuoidea) + Pyraloidea) + Papilionoidea) + Tortricoidea), evidencing the polyphyly of Macrolepidoptera. The Noctuoidea concordantly recovered the familial relationships (((Arctiidae + Lymantriidae) + Noctuidae) + Notodontidae). The tests of optimality strategies, such as exclusion of third codon positions, inclusion of rRNA and tRNA genes, data partitioning, RY recoding approach, and recoding nucleotides into amino acids suggested that the majority of the strategies did not substantially alter phylogenetic topologies or nodal supports, except for the sister relationship between Lycaenidae and Pieridae only in the amino acid dataset, which was in contrast to the sister relationship between Lycaenidae and Nymphalidae in Papilionoidea in the remaining datasets.

Delimiting species using single-locus data and the Generalized Mixed Yule Coalescent approach: a revised method and evaluation on simulated data sets

DNA barcoding-type studies assemble single-locus data from large samples of individuals and species, and have provided new kinds of data for evolutionary surveys of diversity. An important goal of many such studies is to delimit evolutionarily significant species units, especially in biodiversity surveys from environmental DNA samples. The Generalized Mixed Yule Coalescent (GMYC) method is a likelihood method for delimiting species by fitting within- and between-species branching models to reconstructed gene trees. Although the method has been widely used, it has not previously been described in detail or evaluated fully against simulations of alternative scenarios of true patterns of population variation and divergence between species. Here, we present important reformulations to the GMYC method as originally specified, and demonstrate its robustness to a range of departures from its simplifying assumptions. The main factor affecting the accuracy of delimitation is the mean population size of species relative to divergence times between them. Other departures from the model assumptions, such as varying population sizes among species, alternative scenarios for speciation and extinction, and population growth or subdivision within species, have relatively smaller effects. Our simulations demonstrate that support measures derived from the likelihood function provide a robust indication of when the model performs well and when it leads to inaccurate delimitations. Finally, the so-called single-threshold version of the method outperforms the multiple-threshold version of the method on simulated data: we argue that this might represent a fundamental limit due to the nature of evidence used to delimit species in this approach. Together with other studies comparing its performance relative to other methods, our findings support the robustness of GMYC as a tool for delimiting species when only single-locus information is available.