Skip to content

Species Identification by Cytochrome-b, COI Barcoding and 16S rRNA

The mitochondrial and ribosomal markers that distinguish a tiger from a leopard, a beef adulterant from a horse-meat substitute, and one bacterial pathogen from another: cytochrome-b for vertebrates, the COI barcode of life, 16S rRNA for bacteria, ITS regions for fungi, and the BOLD Systems and SILVA reference databases.

Last updated:

Share

When a border seizure, food mislabelling case, or biothreat sample raises the question "what organism is this DNA from?", the answer since the late 1980s has come from mitochondrial and ribosomal molecular markers that are faster, more objective, and court-defensible compared to morphology or biochemical phenotyping.

Key takeaways

  • Cytochrome-b (~350 bp forensic fragment) reliably distinguishes vertebrate species at 5-9% inter-species divergence; the 2013 European horsemeat scandal was confirmed using a 117 bp cyt-b real-time PCR assay targeting horse-specific sequence.
  • The COI Folmer fragment (658 bp, LCO1490/HCO2198 primers) is the BOLD Systems barcode standard for animals; a hit at 98% or above identity is required for species-level identification, and 97-99% is reported as genus-level pending confirmatory analysis.
  • 16S rRNA cannot distinguish Bacillus anthracis from B. cereus (more than 99.5% shared identity); Amerithrax-level attribution requires targeted PCR for the pXO1 and pXO2 virulence plasmids, or whole-genome sequencing.
  • BOLD Systems requires every reference barcode to be linked to a physically vouchered, taxonomist-verified specimen, giving it stronger forensic evidentiary standing than an unchecked GenBank entry for CITES casework.
  • Fungal forensics uses the ITS1-5.8S-ITS2 region (~550-700 bp) curated in the UNITE database; the same ITS barcoding approach is used for forensic mycology and palynological species confirmation.

The logic is elegant. Mitochondria are inherited maternally and accumulate substitutions at a relatively predictable rate, producing regions that are conserved enough to be amplifiable with universal primers yet variable enough to discriminate between species. The cytochrome-b gene (cyt-b) and the cytochrome-c oxidase subunit I gene (COI, also written CO1) carry that combination for vertebrates and the wider animal kingdom respectively. Ribosomal RNA genes carry the same principle in a different genomic compartment: the 16S subunit in bacteria, and the internal transcribed spacer (ITS) regions flanking the 5.8S rRNA gene in fungi, accumulate variation between lineages while remaining recognisable at the universal-primer level.

Reference databases convert a sequence into a species name. BOLD Systems (Barcode of Life Data System, University of Guelph, Canada) holds more than 9 million barcodes for over 250,000 species and is the primary reference for COI identifications. SILVA (hosted at the Max Planck Institute for Marine Microbiology, Bremen, Germany) is the curated benchmark database for 16S, 18S, and 23S/28S rRNA sequences. GenBank at the National Center for Biotechnology Information (NCBI, US) accepts all sequence types and is the broadest reference, though curation quality varies by submission. Database querying relies on the BLAST and sequence-alignment tools described separately. Together these databases underpin a discipline that operates in food-authenticity laboratories in Frankfurt and Singapore, wildlife-seizure labs in Ashland, Oregon and Dehradun, India, and public-health genomics centres responding to biothreat events.

Cytochrome-b: The Vertebrate Workhorse

A fragment of a gene from the mitochondrial respiratory chain has become the standard molecular witness in any case where a vertebrate species must be identified from tissue, blood, or hair.

The cytochrome-b gene encodes a transmembrane protein in the mitochondrial electron transport chain. Its forensic value has nothing to do with its biochemical function. What matters is that it sits on the mitochondrial genome (present in hundreds to thousands of copies per cell, making it recoverable from trace material), it is 1,140 base pairs long in most vertebrates, its ends are flanked by conserved primer-binding sites, and its middle third accumulates substitutions at a rate sufficient to distinguish species while remaining stable enough to be sequenced reliably in degraded samples.

The standard forensic amplicon targets approximately 350-400 base pairs within the 3' portion of the gene, using universal vertebrate primers such as L14724 and H15915 (Kocher et al., 1989) or their shorter derivatives. This fragment reliably discriminates between closely related species. Tiger (Panthera tigris) and leopard (Panthera pardus) differ at cyt-b by roughly 5-7%, enough for unambiguous identification even from degraded trophy material. Elephant species (Loxodonta africana vs Elephas maximus) differ at cyt-b by approximately 9%, which matters when a tusk is seized and its provenance needs to be placed under the Appendix I listing of CITES. Domestic cattle (Bos taurus), water buffalo (Bubalus bubalis), and horse (Equus caballus) can all be distinguished using 359 bp cyt-b fragments with a simple PCR-RFLP approach, which was the basis of the first widely reported horse-meat adulteration testing.

In 2013, a Food Safety Authority of Ireland investigation found horse-meat DNA in 37% of beef-burger products sampled from major UK and Irish retailers. Subsequent European Food Safety Authority-coordinated testing across 27 member states used cyt-b PCR-RFLP and sequencing to confirm that horse meat originating from Romania and Poland had entered the beef supply chain at multiple points. The UK Food Standards Agency used the same marker set. In India, the Food Safety and Standards Authority of India (FSSAI) has adopted species-identification methods for meat products, with cyt-b cited in its 2019 food testing protocols.

Trace sample(tissue / hair /blood)DNA extractioncyt-b PCR(350-400 bp)Sanger or NGSsequencingBLAST vsGenBank
Cytochrome-b forensic workflow: universal primers bracket a 350-400 bp vertebrate amplicon; the sequence is queried against GenBank to return a species hit with percentage identity. A hit above 99% is regarded as species-level identification in most validation studies.

The COI Barcode of Life: 658 Base Pairs That Index the Animal Kingdom

The idea that a single gene fragment could serve as a global species identifier for all multicellular animals seemed extravagant when Paul Hebert proposed it in 2003, two decades of validation later, it underpins wildlife seizure casework on five continents.

The COI barcode is the 658 base-pair Folmer fragment of the cytochrome-c oxidase subunit I gene, amplified with the LCO1490/HCO2198 primer pair (Folmer et al., 1994) or one of several derivative primer sets developed for difficult groups. COI evolves faster than cyt-b in many invertebrate lineages, making it better for species discrimination across the animal kingdom as a whole, though cyt-b remains preferred for mammals and birds where the primer sets are more extensively validated.

The Barcode of Life Data System (BOLD Systems, boldsystems.org) was established at the University of Guelph in 2003 with the explicit purpose of hosting COI barcodes linked to voucher specimens, collection locality, and taxonomic authority. As of 2024, BOLD holds reference barcodes for over 250,000 species spanning insects, fish, birds, reptiles, and mammals. A forensic query against BOLD uses the Species Level Barcode Records (SLBR) dataset: a hit sharing 98% or more sequence identity with a BOLD reference barcode is treated as probable species-level identification in most published validation studies; below 95% indicates genus or higher-level assignment.

Forensic laboratories using COI include the US National Fish and Wildlife Forensics Laboratory (USFWFL) in Ashland, Oregon, which validates sequences against both BOLD and GenBank before reporting a species identification to the agent handling a CITES casework file. In Australia, the Department of Agriculture uses COI for rapid species screening at mail parcels suspected of containing protected wildlife parts. In India, the Wildlife Institute of India (WII, Dehradun) has a molecular systematic laboratory that uses both cyt-b and COI for wildlife-crime casework, and its protocols have been admitted in cases heard before High Courts and the National Green Tribunal. The admissibility standards for sequence-based identification evidence overlap with the broader Daubert framework examined in admissibility and ethics: Daubert, Frye and ELSI.

Fish species identification provides a textbook COI casework context. Sea bass (Dicentrarchus labrax) and Nile perch (Lates niloticus) cannot be distinguished once filleted and skinned. In 2020, a UK Food Standards Agency investigation used COI barcoding to confirm that 35 of 116 fish samples sold as "sea bass" in UK restaurants were misidentified, including cheaper species such as basa (Pangasius bocourti). In a parallel investigation by the Food Safety Authority of Ireland, the same marker confirmed substitution of Atlantic cod (Gadus morhua) with Pacific pollock (Theragra chalcogramma) in fish-and-chip products.

16S rRNA for Bacteria and ITS for Fungi: The Ribosomal Barcode Strategy

Bacteria have no COI gene, but they carry a ribosomal RNA machinery whose conserved architecture and variable loops make the 16S subunit the most widely sequenced single gene in all of biology.

The bacterial ribosome is a 70S particle composed of a 30S small subunit (carrying the 16S rRNA gene, approximately 1,550 base pairs) and a 50S large subunit (carrying 23S and 5S rRNA genes). The 16S rRNA gene is present in one to fifteen copies per bacterial genome, is essential for cell survival, and therefore evolves slowly at its core while accumulating species-diagnostic variation in nine hypervariable regions designated V1 through V9. For forensic species identification, the V3-V4 region (approximately 470 base pairs) amplified with the 341F/806R primer pair provides sufficient resolution to distinguish bacterial species in mixed environmental samples using Illumina MiSeq.

SILVA (silva-project.eu) curates aligned 16S, 18S, and 23S/28S rRNA sequences with quality filtering, chimera checking, and taxonomic annotation updated through a consortium release cycle. The Ribosomal Database Project (RDP, US) and Greengenes (now Greengenes2, 2022, Australia/US) are the other major reference sets. For forensic biothreat attribution, the NCBI RefSeq collection of fully sequenced bacterial genomes provides the highest-resolution reference because it allows whole-genome or multi-locus comparison beyond the 16S region.

The Bacillus anthracis identification problem illustrates the limitation of 16S rRNA at intraspecies level. Bacillus anthracis is genetically near-identical to Bacillus cereus by 16S rRNA sequence. The two species share more than 99.5% 16S sequence identity and cannot be distinguished by 16S alone. The Amerithrax case (Section 4 of this module's microbial forensics topic) required whole-genome sequencing and morphotype analysis to go beyond species and reach the strain-level attribution that built the case. For identification purposes in a clinical or border-control context, however, 16S rRNA reliably places an unknown isolate at the genus level and often at the species level for the approximately 95% of bacterial species that are well-separated in 16S sequence space.

The ITS region, consisting of ITS1, 5.8S rRNA, and ITS2 flanked by the 18S and 28S subunits, is the primary barcode target for fungi. The UNITE database (unite.ut.ee, hosted at the University of Tartu, Estonia) curates fungal ITS sequences linked to species hypotheses. ITS1-5.8S-ITS2 is amplified with the ITS1/ITS4 primer pair (White et al., 1990). Forensic mycology uses ITS barcoding to identify fungal species from environmental samples in cases where mould growth patterns, food contamination, or spore profiles contribute to a forensic question. Palynology, covered in the plant forensics section of this module, uses a different set of nuclear markers (Section 5).

MarkerTarget organism groupStandard ampliconPrimary reference databaseForensic application
Cytochrome-bVertebrates (mammals, birds, reptiles, fish)~350 bp (3' end of gene)GenBank / NCBIWildlife ID, meat adulteration, food fraud
COI (Folmer fragment)Animals (all multicellular)658 bpBOLD SystemsSpecies ID, fish mislabelling, CITES casework
16S rRNA (V3-V4)Bacteria~470 bpSILVA / Greengenes2 / RDPPathogen ID, biothreat attribution, microbiome
ITS1-5.8S-ITS2Fungi~550-700 bp (variable)UNITEForensic mycology, food contamination, palynology
rbcL + matKPlants (vascular)~550 bp + ~870 bpGenBank / DNA barcode of lifePlant species ID, cannabis, paloverde-seed cases
MarkerOrganism scopeAmplicon (bp)Primary databaseSpecies-levelthresholdCytochrome-b(cyt-b)Vertebrates:mammals, birds,reptiles, fish350 to 400 bp (3'end of gene)GenBank / NCBI99% identity(species); 95 to 99%genusCOI (Folmerfragment)All animals(invertebrates,fish, birds,mammals)658 bp (LCO1490 /HCO2198 primers)BOLD Systems(voucher-linked)98% or above(species); below 95%= family16S rRNA (V3-V4)Bacteria (andarchaea)~470 bp (341F /806R primers)SILVA / Greengenes2/ RDPGenus only for B.anthracis vs B.cereus (99.5% shared)ITS(ITS1-5.8S-ITS2)Fungi (and vascularplants withrbcL+matK)550 to 700 bp (ITS1/ ITS4 primers)UNITE (Universityof Tartu)Species hypothesismatch at 97% or above
Four forensic barcoding markers compared: organism scope, diagnostic amplicon length, primary curated database, and the sequence-identity threshold required for a court-grade species-level call. The 16S rRNA row flags genus-only resolution for Bacillus anthracis vs B. cereus.

Reference Databases and the Quality of an Identification

A sequence match is only as reliable as the reference it is matched against, and the gap between a well-curated barcode database and an uncurated GenBank submission can be the difference between a correct species call and a misidentification in a prosecution.

BOLD Systems, SILVA, UNITE, and GenBank differ in their curation strategies, and those differences matter in a court setting where a defence expert can challenge the reference entry. BOLD requires that a submitted barcode be linked to a physical voucher specimen with a determined taxonomy from a specialist. SILVA applies automated and manual quality filters to rRNA sequences and releases versioned reference datasets that can be cited by accession. GenBank accepts submissions from any registered user with minimal quality checks beyond format compliance, meaning that a GenBank BLAST hit carries less inferential weight than a BOLD or SILVA hit unless the reference entry itself is scrutinised.

In a forensic context, the standard practice at accredited laboratories is to query unknown sequences against at least two reference databases and to report the taxonomy only when both return consistent results above a defined identity threshold. The USFWFL in Ashland requires a minimum of 98% identity with a reference sequence in BOLD or GenBank, supported by a phylogenetic tree showing the unknown clustering within a species-level clade, before reporting a species identification in writing to law enforcement. The UK Centre for Environment, Fisheries and Aquaculture Science (CEFAS) follows a similar dual-database requirement for fish species identification.

The CITES framework places molecular identification results into a prosecution context. Under CITES (the Convention on International Trade in Endangered Species of Wild Fauna and Flora, 1973), species are listed in Appendix I (trade prohibited), Appendix II (controlled trade), or Appendix III (national protection). A laboratory report identifying seized material as Panthera tigris (Appendix I) or Elephas maximus (Appendix I) is a material finding in any prosecution under CITES-implementing national legislation: the US Endangered Species Act, the UK Control of Trade in Endangered Species (Enforcement) Regulations, or India's Wildlife Protection Act 1972 (Schedule I). The molecular species identification must therefore meet an evidentiary standard, including documentation of the method, the database version consulted, and the identity score returned.

Food Fraud, Meat Adulteration and the Horsemeat Scandal

The most economically consequential application of molecular species identification has been in the global food supply chain, where the gap between label and contents has proved to be vast and systematic.

The 2013 European horsemeat scandal exposed a supply-chain fraud stretching across at least eight countries. The trigger was a Food Safety Authority of Ireland (FSAI) study that detected horse-specific DNA in beef products using a real-time PCR assay targeting a 117 base-pair fragment of horse cyt-b, run alongside a reference real-time PCR assay for cattle cyt-b to confirm that beef DNA was also present. The quantitative PCR approach (horse-to-beef DNA ratio) showed that some products contained horse meat at proportions exceeding 80% of total meat content. Subsequent investigations using Sanger sequencing and BLAST against GenBank confirmed species identification in regulatory submissions to the European Food Safety Authority (EFSA).

The regulatory framework that governs these investigations differs across jurisdictions. In the European Union, Regulation EC No 1169/2011 requires that the species composition of meat products be accurately labelled. Enforcement is the responsibility of national food safety authorities operating under the Official Controls Regulation (EU) 2017/625. In the UK after Brexit, the Food Safety Act 1990 and the Food Information (England) Regulations 2014 cover the same ground. In India, the Food Safety and Standards (Labelling and Display) Regulations 2020, enforced by FSSAI, require species declaration on meat product labels. In the US, the Federal Meat Inspection Act and USDA Food Safety and Inspection Service regulations require accurate species labelling, and FSIS maintains a molecular methods database for species authentication.

The forensic methodology for food fraud cases parallels wildlife casework. DNA is extracted from a processed food sample (heat treatment, grinding, and pH extremes reduce DNA quality, so short amplicons are preferred). A species-specific PCR or multiplex real-time PCR is used to detect and quantify target species. Sequencing is performed if the real-time PCR result is ambiguous or contested. The result is reported against a validated method, and the validation documentation must accompany any regulatory or prosecutorial submission.

In India, cyt-b and COI have been used in cases involving illegal trade in beef from protected bovine species (gaur, Bos gaurus) and camel (Camelus dromedarius), which is restricted in several states. The Central Forensic Science Laboratory (CFSL, New Delhi) and several state FSLs have validated molecular species identification protocols for meat casework. In Singapore, the Agri-Food and Veterinary Authority has adopted COI barcoding for import control of fish products.

Key terms
Cytochrome-b (cyt-b)
A gene on the mitochondrial genome encoding a component of the electron transport chain. Its 350-400 bp diagnostic fragment is the standard molecular marker for vertebrate species identification in forensic and food-fraud casework.
COI (COX1)
Cytochrome-c oxidase subunit I gene. The 658-base-pair Folmer fragment is the globally adopted barcode-of-life marker for animals, curated in BOLD Systems and used for CITES casework species identification.
16S rRNA
The small-subunit ribosomal RNA gene in bacteria (~1,550 bp). The V3-V4 hypervariable region (~470 bp) is amplified with universal bacterial primers for species identification, microbiome profiling, and biothreat attribution.
ITS region
Internal transcribed spacer region, flanking the 5.8S rRNA gene between the 18S and 28S subunits in fungi and plants. The ITS1-5.8S-ITS2 fragment (~550-700 bp) is the standard fungal barcode, curated in the UNITE database.
BOLD Systems
Barcode of Life Data System (University of Guelph, Canada). Holds over 9 million COI reference barcodes linked to voucher specimens, with mandatory taxonomic verification. Primary reference for animal species identification.
SILVA
Curated ribosomal RNA sequence database (Max Planck Institute for Marine Microbiology, Bremen). Provides quality-filtered 16S, 18S, 23S, and 28S reference alignments for microbial identification; versioned releases allow precise citation in forensic reports.
PCR-RFLP
PCR followed by restriction fragment length polymorphism analysis. A restriction enzyme cuts the amplicon at species-specific recognition sites; the resulting fragment pattern on a gel identifies the species. Rapid and cheap for known species pairs (e.g., cattle vs horse) but requires prior knowledge of the restriction site distribution.
Practice
Question 1 of 1· 0 answered

A forensic laboratory receives a 3 mm piece of dried meat from a border seizure and must confirm whether it is from Panthera tigris (tiger). Which marker and database combination provides the most defensible identification in a CITES prosecution?

Why is 99% COI identity the species-level threshold, and what should an examiner do when a result falls between 97% and 99%?
The 99% threshold is derived from empirical studies of the COI 'barcoding gap': within a species, COI typically varies by less than 2-3%; between species, divergence typically exceeds 3%. The 99% threshold builds in a conservative margin. A BOLD or BLAST result between 97% and 99% sits in the ambiguous zone, possibly an intraspecific variant, an undescribed species, or a poorly represented reference. BOLD's reporting standard for this range is genus-level identification, with a recommendation for confirmatory analysis using additional markers such as cytochrome-b or nuclear microsatellites. Reporting an ambiguous result as a species-level identification would overstate the evidence in a CITES prosecution. See [wildlife forensics and CITES](/topics/forensic-biotechnology/wildlife-forensics-ivory-tiger-pangolin-and-cites) for how this threshold is applied in pangolin and ivory casework.
Why cannot 16S rRNA alone distinguish Bacillus anthracis from Bacillus cereus?
The two species share over 99.5% 16S rRNA sequence identity, too conserved to distinguish them. In a biothreat investigation, 16S is used only as genus-level triage. Confirmatory identification requires targeted PCR for the pXO1 plasmid (encoding protective antigen and lethal factor) and the pXO2 plasmid (encoding the capsule biosynthesis operon), or whole-genome sequencing to place the isolate within the B. anthracis clade. This was the workflow in the Amerithrax investigation, where whole-genome SNP analysis of four morphotype variants distinguished the attack strain from other Ames-strain preparations. The full attribution chain is covered in [microbial forensics and biothreat attribution](/topics/forensic-biotechnology/microbial-forensics-anthrax-letters-and-biothreat-attribution).
What makes BOLD Systems more defensible than a raw GenBank query for COI species identification in court?
BOLD requires every reference barcode to be linked to a physically vouchered specimen confirmed by a recognised taxonomic authority, accessible for re-examination if challenged. GenBank accepts user submissions with minimal quality control, so the database contains mislabelled entries from sequencing errors, contamination, or taxonomic misidentification. A BLAST query against GenBank can return a high-identity hit to an incorrectly annotated entry. In court, a BOLD identification has a traceable evidentiary chain from the query sequence back to a physically verified voucher, whereas an unchecked GenBank hit lacks that chain. Validated wildlife forensics workflows use BOLD as primary reference and GenBank as a supplementary cross-check. The [sequence alignment and BLAST](/topics/forensic-biotechnology/sequence-alignment-blast-genbank-mitomap-empop) topic explains how GenBank QC limitations affect BLAST results more broadly.
Can cytochrome-b typing identify the geographic origin of seized ivory?
Cytochrome-b can distinguish African elephant (Loxodonta africana) from Asian elephant (Elephas maximus) at species level, determining which CITES appendix applies. Geographic origin within Africa requires mtDNA control-region haplotypes or autosomal SNP panels calibrated against population samples from specific geographic clusters. Samuel Wasser and colleagues at the University of Washington built a genotypic assignment map from 16 microsatellite loci that places seized ivory within broad African geographic zones. Cytochrome-b alone cannot achieve this resolution; population-calibrated reference data are required. The distinction matters forensically because it can link seized ivory to specific poaching hotspots and trafficking networks, as demonstrated in the 2014 analysis placing two large seizures in the DRC-Republic of Congo border zone.
Practice
Question 1 of 4· 0 answered

The 2013 European horsemeat scandal was first detected by a real-time PCR assay targeting a 117-base-pair fragment of horse cytochrome-b. Why was a short amplicon preferred over the standard 350-400 bp forensic fragment?

Test yourself on Forensic Biotechnology with free, timed mocks.

Practice Forensic Biotechnology questions

Found this useful? Pass it along.

Share

Spotted an error in this page? Report a correction or read our editorial standards.

Your journey to becoming a forensic professional starts here.

Practice with mock tests, learn from structured notes, and get your questions answered by a global forensic community, all in one place.