Population Genetics, Geographic Provenance, and Stable Isotopes
Population-level genetic assignment and stable isotope analysis can pinpoint where a seized wildlife product originated geographically, turning a tusk or horn into a map coordinate that prosecutors can place at a specific poaching hotspot.
Last updated:
Population genetics and stable isotope analysis are the two principal methods for tracing the geographic origin of seized wildlife products. Genetic assignment uses allele frequency differences that accumulate between isolated populations over generations, expressed as probability distributions across a reference map. Stable isotope ratios in bone, ivory, or horn reflect the local geology, water, and diet of the landscape where the animal lived. When both methods are applied to the same sample and their results converge, the combined evidence is robust enough to direct law enforcement to specific poaching hotspots and to withstand scrutiny in court.
A kilogram of ivory arriving in a Hong Kong seizure carries no reliable label of origin: the crate documentation is typically falsified and the trafficker is several steps removed from the source. The ivory itself, however, carries two independent geographic signals. The elephant's genome encodes population structure accumulated over thousands of years of limited gene flow, and the mineral chemistry of the tusk records the geology and diet of the landscape where the animal lived. Together, these genetic and isotopic signals have made geographic provenance tracing a practical investigative tool.
The genetic side runs through population-genetics software such as STRUCTURE and ADMIXTURE, which decompose a multi-locus genotype into fractional memberships across geographic clusters. Samuel Wasser's laboratory at the University of Washington calibrated this approach specifically for African elephants, building a continental map of SNP allele frequencies that can assign a tusk to a broad region with quantified probability. Rhino horn assignment works similarly, distinguishing South African from Zimbabwean populations at the SNP level.
The isotope side exploits the fact that strontium, carbon, and nitrogen ratios in biological tissue reflect local geology and diet. John Vogel and collaborators built Sr-87/86 isoscapes across sub-Saharan Africa; a ratio measured in ivory can be mapped against that isoscape to identify the most probable origin region. When genetic assignment and isotope provenance agree, the combined evidence is robust enough to hold up in court and direct anti-poaching effort to specific geographic corridors.
By the end of this topic you will be able to:
- Explain how limited gene flow between geographically isolated populations produces allele frequency gradients that can be used for provenance assignment.
- Describe how STRUCTURE and ADMIXTURE software express population assignment as membership coefficients rather than single-point conclusions, and why that matters for court presentation.
- Interpret strontium-87/86 isoscapes and explain how delta-C13 and delta-N15 values encode habitat type and aridity independently of genetic data.
- Apply the dual DNA-isotope framework to a seized specimen, identify what each signal contributes, and recognise the limitations imposed by uneven reference database coverage.
- Distinguish RhODIS individual matching from population-level SNP assignment for rhinoceros horn, and explain how keratin isotope profiling reconstructs an animal's movement history.
- Population structure
- The non-random distribution of genetic variation across subpopulations caused by limited gene flow between them. Different elephant herds separated by geography or habitat accumulate distinct allele frequencies over generations, creating a geographic genetic signal.
- SNP (single nucleotide polymorphism)
- A position in the genome where the nucleotide base differs between individuals. Genome-wide SNP panels are more powerful than STR panels for population assignment because thousands of informative loci can be assayed simultaneously.
- STRUCTURE / ADMIXTURE
- Software tools that infer population structure by estimating each individual's fractional ancestry (membership coefficient) across K genetic clusters. Used to assign unknown samples to source populations and to visualise admixture between groups.
- Isoscape
- A continental or global map of expected stable isotope ratios: typically Sr-87/86, delta-C13, or delta-N15: interpolated from measured values at known locations. Used to geolocate a biological sample by matching its measured ratio to the map.
- Strontium-87/86 ratio
- The ratio of two strontium isotopes in a sample, which mirrors the ratio in local bedrock and water. Because strontium substitutes for calcium in bone and dentine, the ratio measured in ivory reflects the geological substrate where the elephant lived.
- Delta-C13 / Delta-N15
- Stable carbon and nitrogen isotope ratios in tissue, expressed as deviation from a reference standard (delta notation). They reflect diet composition and trophic level, and can distinguish savanna from forest habitat and wet from dry ecosystems.
Genetic population structure as a geographic signal
Gene flow between elephant populations is constrained by rivers, mountains, farmland, and national park boundaries. Where populations are isolated long enough, their allele frequencies drift apart. The result is a gradient of genetic differentiation that correlates with geography. Populations in East Africa cluster separately from those in West Africa; Southern African populations form a third cluster; Central African forest elephants form yet another. This structure is the foundation of genetic provenance assignment.
Samuel Wasser's team at the University of Washington published a geographic assignment tool for elephant ivory using microsatellite and later SNP markers. The approach builds a map of allele frequencies from reference samples of known origin, then calculates the likelihood that an unknown sample's genotype was drawn from each map cell. The result is a probability surface across Africa rather than a single point, which is the correct way to present geographic uncertainty in court.
STRUCTURE and ADMIXTURE: the software backbone
STRUCTURE (Pritchard et al., 2000) and ADMIXTURE (Alexander et al., 2009) take as input a matrix of genotypes: alleles at multiple loci across multiple individuals: and output membership coefficients: the fraction of each individual's genome that traces to each of K hypothetical clusters. Running the analysis at increasing values of K and selecting the one that best fits the data (using delta-K or cross-validation error) reveals the natural number of genetically distinct groups in the dataset.

For a seized specimen the workflow is: extract DNA, generate the same marker panel as the reference dataset, and project the unknown individual into the ADMIXTURE or STRUCTURE framework. The membership coefficient for each cluster is the assignment probability. A tusk that shows 87% membership in the East African cluster and 13% in the Southern African cluster is assigned to East Africa with 87% confidence: a quantified, auditable statistic that a court can evaluate.
| Feature | STRUCTURE | ADMIXTURE |
|---|---|---|
| Algorithm | Bayesian MCMC | Maximum likelihood |
| Speed | Slow for large datasets | Fast; scales to thousands of loci |
| Output | Membership coefficients + posterior probability | Membership coefficients + cross-validation error |
| Best for | Small panels, exploratory analysis | Genome-wide SNP datasets |
Stable strontium isotopes and the Vogel isoscapes
Strontium has four naturally occurring isotopes. The ratio of Sr-87 to Sr-86 varies across Earth's surface because Sr-87 is a decay product of rubidium-87, and rubidium concentrations differ by rock type. Old granitic cratons carry high Sr-87/86 ratios; younger volcanic basalts carry lower ones. Plants take up strontium from soil water in proportion to local ratios; elephants drink local water and eat local vegetation, so their ivory records the ratio of the landscapes they inhabited.
John Vogel and collaborators mapped Sr-87/86 ratios across sub-Saharan Africa by measuring bones and teeth of animals of known location, creating an isoscape: a spatial map of expected values. Measuring the ratio in a section of ivory and overlaying it on the isoscape produces a probability map of origin locations. Vogel's work showed that strontium alone could exclude large portions of Africa and in some cases narrow provenance to a region of a few hundred kilometres.

Delta-C13 and delta-N15: diet and habitat tracing
Stable carbon isotopes (delta-C13) in organic tissue reflect the photosynthetic pathway of the plants the animal ate. C4 grasses, which dominate open savanna, carry less negative delta-C13 values than C3 forest vegetation. An elephant that lived primarily in open savanna will show less negative ivory delta-C13 than one from dense forest. This separates savanna from forest habitats with a clean biochemical signal independent of the genetic data.
Nitrogen isotopes (delta-N15) increase up the food chain and also reflect aridity: arid-zone plants and soils have higher N-15 enrichment than humid ones, partly through ammonia volatilisation from dry soils. For elephants, which are herbivores, N-15 mainly tracks aridity and soil fertility rather than trophic level. Combining C13 and N15 with Sr-87/86 produces a three-isotope signature that, when mapped against an isoscape calibrated for all three, can constrain origin more tightly than any one isotope alone.
Rhino horn population assignment
White rhino (Ceratotherium simum simum) populations in South Africa and Zimbabwe show clear SNP-level genetic differentiation, reflecting their management history as semi-isolated populations within national parks. This makes population assignment feasible at the country level, which is directly relevant to legal jurisdiction. A horn that assigns to the South African population implicates South African protected areas; one that assigns to Zimbabwe points to a different national authority and different trafficking route.
Rhino horn is keratin, not bone, and it grows continuously. This means isotope sampling along the horn length records a time-series just as tusk sampling does. Stable isotopes in rhino horn keratin have been used to reconstruct the animal's movement between landscapes and to estimate how recently it was killed: important when prosecutors need to show that a horn is from a recently poached animal rather than a legal stockpile.
- RhODIS (Rhino DNA Index System): South Africa's national rhino DNA database, holding STR profiles from living animals, dehorned animals, and carcasses. It has been used to match seized horn to specific animals in court proceedings.
- Population assignment limitation: rhino translocation between countries for conservation purposes can blur the population genetic signal. An animal translocated from South Africa to Kenya will carry South African allele frequencies despite living in Kenya: a complication that should be addressed in the report.
The dual DNA-isotope approach
The most powerful provenance reports combine genetic assignment with isotope data on the same sample. The two methods are orthogonal: population genetics reads inherited sequence variation accumulated over thousands of generations, while stable isotopes read dietary and environmental chemistry incorporated during the individual's lifetime. When they converge on the same origin region, the joint probability of both patterns arising by coincidence from a different location is very small.
Researchers working on large ivory seizures, particularly the cargo analyses published by Wasser's group between 2007 and 2015 examining hundreds of tusks from major seizures, found that genetic assignment alone could identify the primary source populations and that isotope data corroborated those assignments while providing additional within-region discrimination. The combined approach has influenced wildlife trafficking investigations by identifying poaching hotspots in specific national parks rather than just broad geographic regions, directing law enforcement resources more precisely.
What biological basis allows genetic population assignment of elephant ivory to geographic regions?
Key Takeaways
- Geographic population assignment exploits the genetic differentiation that develops when elephant (or rhino) populations are isolated by landscape barriers, producing allele frequency gradients that correlate with geography.
- The Wasser et al. assignment tool and the STRUCTURE/ADMIXTURE software framework express assignments as probability distributions across geographic clusters, not single-point conclusions, which is the correct form for court presentation.
- Sr-87/86 ratios in ivory reflect local bedrock geology; delta-C13 separates savanna from forest habitat; delta-N15 tracks aridity. Together they form a multi-isotope signature that can constrain geographic origin independently of genetic data.
- Rhino horn assignment using SNPs can distinguish South African from Zimbabwean populations; RhODIS provides a reference STR database for individual matching; isotope profiling along the horn axis reconstructs the animal's movement and living history.
- A dual genetic-isotope approach is stronger than either method alone because the two signals are biologically independent; gaps in reference database coverage must be reported honestly rather than hidden behind a confident-sounding assignment.
How does geographic assignment work for elephant ivory?
What is the difference between STRUCTURE and ADMIXTURE?
How does strontium isotope analysis help in ivory provenance?
Why is a dual DNA-isotope approach stronger than either alone?
Can these methods be applied to rhino horn?
Test yourself on Wildlife Forensics with free, timed mocks.
Practice Wildlife Forensics questionsSpotted an error in this page? Report a correction or read our editorial standards.