Twenty-three years of PCR-based seafood authentication assay development: What have we learned?

Seafood is a prime target for fraudulent activities due to the complexity of its supply chain, high demand, and difficult discrimination among species once morphological characteristics are removed. Instances of seafood fraud are expected to increase due to growing demand. This manuscript reviews the application of DNA-based methods for commercial fish authentication and identification from 2000 to 2023. It explores (1) the most common types of commercial fish used in assay development, (2) the type of method used, (3) the gene region most often targeted, (4) provides a case study of currently published assays or primer-probe pairs used for DNA amplification, for specificity, and (5) makes recommendations for ensuring standardized assay-based reporting for future studies. A total of 313 original assays for the detection and authentication of commercial fish species from 191 primary articles published over the last 23 years were examined. The most explored DNA-based method was real-time polymerase chain reaction (qPCR), followed by DNA sequencing. The most targeted gene regions were cytb (cytochrome b) and COI (cytochrome c oxidase 1). Tuna was the most targeted commercial fish species. A case study of published tuna assays (n = 19) targeting the cytb region found that most assays were not species-specific through in silico testing. This was conducted by examining the primer mismatch for each assay using multiple sequence alignment. Therefore, there is need for more standardized DNA-based assay reporting in the literature to ensure specificity, reproducibility, and reliability of results. Factors, such as cost, sensitivity, quality of the DNA, and species, should be considered when designing assays.

A Measure of the DNA Barcode Gap for Applied and Basic Research.

DNA barcoding has largely established itself as a mainstay for rapid molecular taxonomic identification in both academic and applied research. The use of DNA barcoding as a molecular identification method depends on a "DNA barcode gap"-the separation between the maximum within-species difference and the minimum between-species difference. Previous work indicates the presence of a gap hinges on sampling effort for focal taxa and their close relatives. Furthermore, both theory and empirical work indicate a gap may not occur for related pairs of biological species. Here, we present a novel evaluation approach in the form of an easily calculated set of nonparametric metrics to quantify the extent of proportional overlap in inter- and intraspecific distributions of pairwise differences among target species and their conspecifics. The metrics are based on a simple count of the number of overlapping records for a species falling within the bounds of maximum intraspecific distance and minimum interspecific distance. Our approach takes advantage of the asymmetric directionality inherent in pairwise genetic distance distributions, which has not been previously done in the DNA barcoding literature. We apply the metrics to the predatory diving beetle genus Agabus as a case study because this group poses significant identification challenges due to its morphological uniformity despite both relative sampling ease and well-established taxonomy. Results herein show that target species and their nearest neighbor species were found to be tightly clustered and therefore difficult to distinguish. Such findings demonstrate that DNA barcoding can fail to fully resolve species in certain cases. Moving forward, we suggest the implementation of the proposed metrics be integrated into a common framework to be reported in any study that uses DNA barcoding for identification. In so doing, the importance of the DNA barcode gap and its components for the success of DNA-based identification using DNA barcodes can be better appreciated.

Fecal DNA metabarcoding helps characterize the Canada jay's diet and confirms its reliance on stored food for winter survival and breeding.

Accurately determining the diet of wild animals can be challenging if food items are small, visible only briefly, or rendered visually unidentifiable in the digestive system. In some food caching species, an additional challenge is determining whether consumed diet items have been previously stored or are fresh. The Canada jay (Perisoreus canadensis) is a generalist resident of North American boreal and subalpine forests with anatomical and behavioural adaptations allowing it to make thousands of arboreal food caches in summer and fall that are presumably responsible for its high winter survival and late winter/early spring breeding. We used DNA fecal metabarcoding to obtain novel information on nestling diets and compiled a dataset of 662 published and unpublished direct observations or stomach contents identifications of natural foods consumed by Canada jays throughout the year. We then used detailed natural history information to make informed decisions on whether each item identified to species in the diets of winter adults and nestlings was best characterized as 'likely cached', 'likely fresh' (i.e., was available as a non-cached item when it appeared in a jay's feces or stomach), or 'either possible'. Of the 87 food items consumed by adults in the winter, 39% were classified as 'likely cached' and 6% were deemed to be 'likely fresh'. For nestlings, 29% of 125 food items identified to species were 'likely cached' and 38% were 'likely fresh'. Our results support both the indispensability of cached food for Canada jay winter survival and previous suggestions that cached food is important for late winter/early spring breeding. Our work highlights the value of combining metabarcoding, stomach contents analysis, and direct observations to determine the cached vs. non-cached origins of consumed food items and the identity of food caches, some of which could be especially vulnerable to degradation through climate change.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae).

BACKGROUND: Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species. AIM:  Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species. RESULTS: We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample. CONCLUSIONS: Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

Molecular Acquisition, Cleaning and Evaluation in R (MACER) - A tool to assemble molecular marker datasets from BOLD and GenBank.

Molecular sequence data is an essential component for many biological fields of study. The strength of these data is in their ability to be centralised and compared across research studies. There are many online repositories for molecular sequence data, some of which are very large accumulations of varying data types like NCBI's GenBank. Due to the size and the complexity of the data in these repositories, challenges arise in searching for data of interest. While data repositories exist for molecular markers, taxa and other specific research interests, repositories may not contain, or be suitable for, more specific applications. Manually accessing, searching, downloading, accumulating, dereplicating and cleaning data to construct project-specific datasets is time-consuming. In addition, the manual assembly of datasets presents challenges with reproducibility. Here, we present the MACER package to assist researchers in assembling molecular datasets and provide reproducibility in the process.

Validation of an Effective Protocol for Culicoides Latreille (Diptera: Ceratopogonidae) Detection Using eDNA Metabarcoding.

eDNA metabarcoding is an effective molecular-based identification method for the biosurveillance of flighted insects. An eDNA surveillance approach maintains specimens for secondary morphological identification useful for regulatory applications. This study identified Culicoides species using eDNA metabarcoding and compared these results to morphological identifications of trapped specimens. Insects were collected using ultraviolet (UV) lighted fan traps containing a saturated salt (NaCl) solution from two locations in Guelph, Ontario, Canada. There were forty-two Culicoides specimens collected in total. Molecular identification detected four species, C. biguttatus, C. stellifer, C. obsoletus, and C. mulrennani. Using morphological identification, two out of these four taxonomic ranks were confirmed at the species level (C. biguttatus and C. stellifer) and one was confirmed at the subgenus level (Avaritia [C. obsoletus]). No molecular detection of Culicoides species occurred in traps with an abundance of less than three individuals per taxon. The inconsistency in identifying Culicoides specimens to the species level punctuates the need for curated DNA reference libraries for Culicoides. In conclusion, the saturated salt (NaCl) solution preserved the Culicoides' morphological characteristics and the eDNA.

Biosurveillance for invasive insect pest species using an environmental DNA metabarcoding approach and a high salt trap collection fluid.

With the increase in global trade and warming patterns, the movement, introduction, and establishment of non-native insect species has increased. A rapid and effective early detection biosurveillance program to identify species of concern is needed to reduce future impacts and costs associated with introduced non-native species. One of the challenges facing insect surveillance trapping methods is the sheer volume of individual specimens in the collections. Although molecular identification methods are improving, they currently have limitations (e.g., destructive processing of specimens) and a protocol addressing these limitations can support regulatory applications that need morphological evidence to corroborate molecular data.The novel protocol presented here uses a metabarcoding approach to amplify environmental DNA from a saturated salt solution trap fluid, which retains trap specimens for downstream morphological identifications. The use of a saturated salt solution to preserve specimens in traps addresses issues with the high evaporation rate of ethanol in traps, and public safety concerns with other fluid preservation options with unattended traps in public settings.Using a metabarcoding approach, a 407-nucleotide segment of the cytochrome c oxidase subunit 1 (COI) animal barcode region was successfully amplified from Lindgren funnel trap collection fluids. These traps were placed in forested areas to survey for wood-boring beetles of regulatory concern. Our results displayed successful amplification of target taxa, including the molecular identification of the Japanese Beetle Popillia japonica, a species regulated in Canada. A second species, Anisandrus maiche, recently introduced to North America, was identified in every trap. The genus Lymantria, which contains numerous species of concern to North American woodlands, was also detected. Also, there were six other species identified of interest due to their potential impacts on native and crop flora and fauna.Our results show how this protocol can be used as an efficient method for the surveillance of insects using a trap with a saturated salt solution and eDNA metabarcoding to detect species of regulatory concern.

Molecular Detection Mapping and Analysis Platform for R (MDMAPR) facilitating the standardization, analysis, visualization, and sharing of qPCR data and metadata.

Quantitative polymerase chain reaction (qPCR) has been used as a standard molecular detection tool in many scientific fields. Unfortunately, there is no standard method for managing published qPCR data, and those currently used generally focus on only managing raw fluorescence data. However, associated with qPCR experiments are extensive sample and assay metadata, often under-examined and under-reported. Here, we present the Molecular Detection Mapping and Analysis Platform for R (MDMAPR), an open-source and fully scalable informatics tool for researchers to merge raw qPCR fluorescence data with associated metadata into a standard format, while geospatially visualizing the distribution of the data and relative intensity of the qPCR results. The advance of this approach is in the ability to use MDMAPR to store varied qPCR data. This includes pathogen and environmental qPCR species detection studies ideally suited to geographical visualization. However, it also goes beyond these and can be utilized with other qPCR data including gene expression studies, quantification studies used in identifying health dangers associated with food and water bacteria, and the identification of unknown samples. In addition, MDMAPR's novel centralized management and geospatial visualization of qPCR data can further enable cross-discipline large-scale qPCR data standardization and accessibility to support research spanning multiple fields of science and qPCR applications.

The Molecular Data Organization for Publication (MDOP) R package to aid the upload of data to shared databases.

Molecular identification methods, such as DNA barcoding, rely on centralized databases populated with morphologically identified individuals and their referential nucleotide sequence records. As molecular identification approaches have expanded in use to fields such as food fraud, environmental surveys, and border surveillance, there is a need for diverse international data sets. Although central data repositories, like the Barcode of Life Datasystems (BOLD), provided workarounds for formatting data for upload, these workarounds can be taxing on researchers with few resources and limited funding. To address these concerns, we present the Molecular Data Organization for Publication (MDOP) R package to assist researchers in uploading data to public databases. To illustrate the use of these scripts, we use the BOLD system as an example. The main intent of this writing is to assist in the movement of data, from academic, governmental, and other institutional computer systems, to public locations. The movement of these data can then better contribute to the global DNA barcoding initiative and other global molecular data efforts.

Using DNA barcoding to improve invasive pest identification at U.S. ports-of-entry.

Interception of potential invasive species at ports-of-entry is essential for effective biosecurity and biosurveillance programs. However, taxonomic assessment of the immature stages of most arthropods is challenging; characters for identification are often dependent on adult morphology and reproductive structures. This study aims to strengthen the identification of such specimens through DNA barcoding, with a focus on microlepidoptera. A sample of 241 primarily immature microlepidoptera specimens intercepted at U.S. ports-of-entry from 2007 to 2011 were selected for analysis. From this sample, 201 COI-5P sequences were generated and analyzed for concordance between morphology-based and DNA-based identifications. The retrospective analysis of the data over 10 years (2009 to 2019) using the Barcode of Life Data (BOLD) system demonstrates the importance of establishing and growing DNA barcode reference libraries for use in specimen identification. Additionally, analysis of specimen identification using public data (43.3% specimens identified) vs. non-public data (78.6% specimens identified) highlights the need to encourage researchers to make data publicly accessible. DNA barcoding surpassed morphological identification with 42.3% (public) and 66.7% (non-public) of the sampled specimens achieving a species-level identification, compared to 38.3% species-level identification by morphology. Whilst DNA barcoding was not able to identify all specimens in our dataset, its incorporation into border security programs as an adjunct to morphological identification can provide secondary lines of evidence and lower taxonomic resolution in many cases. Furthermore, with increased globalization, database records need to be clearly annotated for suspected specimen origin versus interception location.

Barcode-based species delimitation in the marine realm: a test using Hexanauplia (Multicrustacea: Thecostraca and Copepoda).

DNA barcoding has been used successfully for identifying specimens belonging to marine planktonic groups. However, the ability to delineate species within taxonomically diverse and widely distributed marine groups, such as the Copepoda and Thecostraca, remains largely untested. We investigate whether a cytochrome c oxidase subunit I (COI-5P) global pairwise sequence divergence threshold exists between intraspecific and interspecific divergences in the copepods plus the thecostracans (barnacles and allies). Using publicly accessible sequence data, we applied a graphical method to determine an optimal threshold value. With these thresholds, and using a newly generated planktonic marine data set, we quantify the degree of concordance using a bidirectional analysis and discuss different analytical methods for sequence-based species delimitation (e.g., BIN, ABGD, jMOTU, UPARSE, Mothur, PTP, and GMYC). Our results support a COI-5P threshold between 2.1% and 2.6% p-distance across methods for these crustacean taxa, yielding molecular groupings largely concordant with traditional, morphologically defined species. The adoption of internal methods for clustering verification enables rapid biodiversity studies and the exploration of unknown faunas using DNA barcoding. The approaches taken here for concordance assessment also provide a more quantitative comparison of clustering results (as contrasted with "success/failure" of barcoding), and we recommend their further consideration for barcoding studies.

The "Lillie transition": models of the onset of saltatory conduction in myelinating axons.

Almost 90 years ago, Lillie reported that rapid saltatory conduction arose in an iron wire model of nerve impulse propagation when he covered the wire with insulating sections of glass tubing equivalent to myelinated internodes. This led to his suggestion of a similar mechanism explaining rapid conduction in myelinated nerve. In both their evolution and their development, myelinating axons must make a similar transition between continuous and saltatory conduction. Achieving a smooth transition is a potential challenge that we examined in computer models simulating a segmented insulating sheath surrounding an axon having Hodgkin-Huxley squid parameters. With a wide gap under the sheath, conduction was continuous. As the gap was reduced, conduction initially slowed, owing to the increased extra-axonal resistance, then increased (the "rise") up to several times that of the unmyelinated fiber, as saltatory conduction set in. The conduction velocity slowdown was little affected by the number of myelin layers or modest changes in the size of the "node," but strongly affected by the size of the "internode" and axon diameter. The steepness of the rise of rapid conduction was greatly affected by the number of myelin layers and axon diameter, variably affected by internode length and little affected by node length. The transition to saltatory conduction occurred at surprisingly wide gaps and the improvement in conduction speed persisted to surprisingly small gaps. The study demonstrates that the specialized paranodal seals between myelin and axon, and indeed even the clustering of sodium channels at the nodes, are not necessary for saltatory conduction.