Phylogeography, Population History and Conservation Genetics of Jaguars (, Mammalia, Felidae)

Molecular Ecology (2001)
10

, 6579
2001 Blackwell Science Ltd

Blackwell Science, Ltd

Phylogeography, population history and conservation
genetics of jaguars (

Panthera onca

, Mammalia, Felidae)

EDUARDO EI ZI RI K, * J AE- HEUP KI M, * MARI LYN MENOTTI - RAYMOND, *
PETER G. CRAWSHAW J R. , STEPHEN J . O BRI EN* and WARREN E. J OHNSON*
*

Laboratory of Genomic Diversity, National Cancer Institute

FCRDC, Frederick, MD 217021201, USA,

Department of Biology,
University of Maryland at College Park, USA,

CENAP/IBAMA, C.P. 295, Araoiaba da Serra, SP 18190000, Brazil

Abstract
The jaguar (

Panthera onca

), the largest felid in the American Continent, is currently threatened
by habitat loss, fragmentation and human persecution. We have investigated the genetic
diversity, population structure and demographic history of jaguars across their geographical
range by analysing 715 base pairs of the mitochondrial DNA (mtDNA) control region and 29
microsatellite loci in

40 individuals sampled from Mexico to southern Brazil. Jaguars
display low to moderate levels of mtDNA diversity and medium to high levels of microsatellite
size variation, and show evidence of a recent demographic expansion. We estimate that extant
jaguar mtDNA lineages arose 280 000510 000 years ago (95% CI 137 000830 000 years ago),
a younger date than suggested by available fossil data. No strong geographical structure
was observed, in contrast to previously proposed subspecic partitions. However, major
geographical barriers such as the Amazon river and the Darien straits between northern
South America and Central America appear to have restricted historical gene ow in this
species, producing measurable genetic differentiation. Jaguars could be divided into four
incompletely isolated phylogeographic groups, and further sampling may reveal a ner
pattern of subdivision or isolation by distance on a regional level. Operational conservation
units for this species can be dened on a biome or ecosystem scale, but should take into
account the historical barriers to dispersal identied here. Conservation strategies for jaguars
should aim to maintain high levels of gene ow over broad geographical areas, possibly
through active management of disconnected populations on a regional scale.

Keywords

: conservation, evolution, microsatellites, mitochondrial DNA,

Panthera onca
Received 23 March 2000; revision received 16 August 2000; accepted 16 August 2000

Introduction

The jaguar (

Panthera onca

) is the largest Neotropical felid,
measuring up to 2 m in total length and weighing up to
120 kg (Kitchener 1991). It is closely related to the other
four big cats in the genus

Panthera

(Salles 1992; Johnson &
OBrien 1997), and is the only extant representative of this
genus in the Americas. It appears to have diverged from a
common ancestor with the other

Panthera

species at least
1.5 million years ago (Ma) and entered the American con-
tinent through Beringia in the early Pleistocene (Turner &
Antn 1997). Jaguar fossils have been reported in North
America as far north as Nebraska and Washington
(Seymour 1989; Turner & Antn 1997), the oldest denitive
specimens dating from approximately 850 000 years ago
(K. Seymour, personal communication).
In historical times jaguars ranged from the southern
USA throughout Central and South America to the
Argentinean Patagonia. Since the mid-1900s, however,
they have been extirpated from broad portions of their
original distribution and currently subsist in fragmented
populations of variable size (Nowell & Jackson 1996).
Jaguars occur in habitats ranging from dense rainforest
to more open scrubland and seasonally ooded marshes,
and prey mainly on large mammals such as peccary, deer,
tapir, and capybara, and in some areas also on reptiles
such as turtles and caiman (Oliveira 1994). Around human
agricultural lands they often prey on livestock, which
generates chronic conicts with ranchers and widespread

Correspondence: Warren E. Johnson. Fax: (1301) 846 6327; E-mail:
johnsonw@mail.ncifcrf.gov

MEC1144.fm Page 65 Friday, December 8, 2000 3:25 PM

66

E. EI ZI RI K

ET AL

.

2001 Blackwell Science Ltd,

Molecular Ecology

, 10, 6579

persecution of this species in many regions. Although
jaguars have been heavily hunted in the past as a trophy
or for their pelts, their survival is currently threatened
mostly by severe loss and fragmentation of habitat and
direct persecution by ranchers (Nowell & Jackson 1996).
The species as a whole is currently listed in Appendix I
(endangered species) of the Convention on International
Trade of Endangered Species of Wild Fauna and Flora
(CITES), and is considered to be vulnerable by IUCN
(Wozencraft 1993). Some remnant populations, however,
are known to be small and isolated, and can be considered
to be severely endangered (Medelln

et al.

, in press).
In order to devise adequate conservation and manage-
ment strategies for an endangered species, it is important
to incorporate a reliable understanding of its population
structure and history, the existence of demographic partitions
throughout its geographical range, and a characterization
of the distribution of its genetic diversity (Avise 1989;
OBrien 1994). Jaguars are currently divided into eight
subspecies based on classical studies of morphology
(reviewed by Seymour 1989), but a recent revision of skull
characters found no signicant geographical differenti-
ation in these traits (Larson 1997). However, there has been
no comparable attempt to describe the geographical
distribution of genetic diversity in jaguars, nor to use
molecular markers to clarify its population structure and
evolutionary history.
In the present paper, we report a study of genetic variation
in jaguars sampled throughout their geographical range,
applying mitochondrial DNA (mtDNA) control region
sequences and microsatellite loci size variation to address
issues of phylogeography and population history in
this species. Results from phylogenetic and population
genetic analyses of these data are used to draw inferences
on past demographic processes, clarify taxonomic issues,
and provide insight relevant for conservation efforts.
Table 1 Samples of Panthera onca analysed in the present study
ID Geographic Origin Institution/Contact
Pon16 Guatemala ARCAS
Pon20 Chepo, Panama Prov., Panama Summit Zoo
Pon21 Upala, Alajuela Prov., Costa Rica Las Pumas
Pon22 Puerto Limon Prov., Costa Rica Simon Bolivar Zoo
Pon23 Atlantico Sur Department, Nicaragua Juigalpa Zoo
Pon24 Nicaragua Radio Antenna, Managua
Pon25 Rio San Juan Department, Nicaragua Managua Zoo
Pon26 Guatemala Aurora Zoo
Pon28 Peru (Probably) Mendoza Zoo
Pon29 Southern South America (Probably) Montevideo Zoo
Pon30 Yucatan State, Mexico Proyecto Balam
Pon31 San Luis Potosi State, Mexico Leon Zoo
Pon34
, 35
Mexico, unknown State Cd.Victoria Zoo, Aragon Zoo

Pon37 Yucatan State, Mexico Zoologico Centenario
Pon40*, 41, 42, 44 Mato Grosso do Sul State, SW Brazil Proj. Carnvoros, CENAP/IBAMA
Pon43 P.N. Iguau, Paran State, S Brazil Proj. Carnvoros, CENAP/IBAMA
Pon46 Gois State, Brazil Zoolgico de Goinia
Pon47 Par State, N Brazil Zoolgico de Goinia
Pon48
Brazil, unknown State Zoolgico de Goinia

Pon49, 50
Chaco, Paraguay Itaipu, Paraguay

Pon51, 52 Jalisco State, Mexico Idaho State University
Pon54
Amazonas State, Venezuela Las Delicias

Pon56*, 58 Falcon State, Venezuela Las Delicias
Pon59 Barinas State, Venezuela Asociacion Ganaderas
Pon60 Bolivar State, Venezuela Merida
Pon61 Venezuela, unknown State Barquisimento
Pon62
, 63 Bolivar State, Venezuela Valencia Castellito

Pon64 Sta. Cruz Department, Bolivia Zoologico de Santa Cruz
Pon65, 66 N. Bolivia Zoologico de Santa Cruz
Pon67
Venezuela, unknown State PROFAUNA

Pon68 French Guiana Prog. Faune Sauvage/J.C. Vi
Pon71 Brazil, unknown State Zoolgico de Paranagu/R.N.Morais
Pon72 Brazil, unknown State CENAP/IBAMA
Pon73 Brazil, unknown State CENAP/IBAMA
Pon75 Sinaloa State, Mexico Brian Miller
, samples typed only for microsatellites; *, samples typed only for mtDNA.


CONSERVATI ON GENETI CS OF J AGUARS

67


Molecular Ecology

, 10, 6579

Materials and methods

Sample collection and laboratory procedures

Samples of 44 jaguar individuals were obtained from
across most of the species range (Table 1, Fig. 1). Whole
blood preserved in a saturated salt solution (100 m

m

Tris,
100 m

m

EDTA, 2% SDS) was collected from wild-caught
individuals captured in eld ecology projects; blood
(separated into plasma, red blood cell and leucocyte
phases) and skin biopsies were obtained from animals
kept in captive breeding Institutions, most of which had
known geographical origin; and skin fragments were
collected from a museum pelt (sample Pon67). DNA
extraction was performed using standard techniques: salt
precipitation (Montgomery & Sise 1990) was used for
leucocyte pellets; phenol/chloroform (Sambrook

et al

. 1989)
was used for whole blood and primary broblast cultures;
and a silica-guanidium technique (Pbo

et al

. 1988; Boom

et al

. 1990) was employed with the museum material.
The 5

portion (including the First Hypervariable
Segment, HVS-I) and the central conserved portion of
the mtDNA control region (the structure of the felid
mtDNA control region is described elsewhere: Lopez

et al

.
1996; Eizirik

et al

. 1998; J.-H. Kim

et al

. unpublished) were
amplied by polymerase chain reaction (PCR) (Saiki

et al

.
1985) from genomic DNA. For the 5

segment, human
primers L15997 [modied to match the domestic cat ref-
erence sequence (Lopez

et al

. 1996)] and H16498 (Ward

et al

. 1991) were used, and for the central portion primers
PAN-CCR-F (CTCAACTGTCCGAAAGTGCTT) and PAN-
CCR-R (CCTGTGGAAGCAATAGGAATT) ( J.-H. Kim

et al

.
unpublished data) were employed. PCR reactions were
performed using standard buffer conditions, 1.5 m

m

MgCl

2

, and 3336 cycles of the following steps: 94

C
for 45 s, 53

C for 45 s, and 72

C for 1 min, preceded by
3 min of initial denaturing, and followed by 10 min of nal
extension. PCR products were puried using

centricon

-
100 (Amicon), and cycle-sequenced from both ends, using
either ABI FS Dye Primer kits (in which case PCR primers
Fig. 1 Map showing historical and current
geographical distribution of jaguars (modied
from Oliveira 1994), approximate sample
collection sites (black circles), and the location
of the Amazon river. Numbers next to circles
are sample identication labels of jaguar
individuals (number after Pon in Table 1)
from each area. Boxes indicate individuals
from the same region. Only individuals with
know regional origin are shown.


68

E. EI ZI RI K

ET AL

.


Molecular Ecology

, 10, 6579

had been engineered to contain M13 tails) or ABI Big Dye
Terminator kits; resulting fragments were analysed in an
ABI 373 A automated sequencer.
The 5

portion of the mtDNA control region of felids
contains a rather long tandem repetitive series, which leads
to size variation and heteroplasmy in several species,
including jaguars (Lopez

et al

. 1996; Eizirik

et al

. 1998;
J.-H. Kim

et al

. unpublished data). Due to the difculty
in sequencing through this array, only one strand was
sequenced on each ank of the repeats. Sequences were
visually inspected and corrected using

sequencher

(Gene
Codes Co.); sites or segments from which sequence could
not be unambiguously scored after three attempts (inde-
pendent PCR and sequencing reactions) were treated
as missing information and excluded from the analysis.
Similar care was taken with the central segment, although
in this case both strands were available for analysis. As part
of the mtDNA of cats of the

Panthera

genus is included in
a nuclear transposition ( Johnson

et al

. 1996), our data were
compared to reference sequences obtained from isolated
mitochondrial and nuclear extracts (which exhibited several
nucleotide differences; J.-H. Kim

et al

. unpublished data)
to ascertain that the haplotypes included in this study
were indeed of mitochondrial origin.
Thirty-ve microsatellite loci developed originally
for the domestic cat (Menotti-Raymond

et al

. 1999) were
amplied by PCR from jaguar genomic DNA using uore-
scently labelled primers and the standard conditions
described by Menotti-Raymond

et al

. (1999). The products
were analysed in an ABI 310 automated sequencer. To
assess amplication efciency and yield, and to establish
the size range for each microsatellite in the jaguar, an initial
run was performed with three individuals for each locus
separately. Twenty-nine loci (FCA005, FCA008, FCA026,
FCA043, FCA075, FCA077, FCA90, FCA091, FCA094,
FCA453 and FCA678 Menotti-Raymond

et al

. 1999) were
selected for use in this study. PCR reactions were carried
out for each locus separately, and products from 3 to 5
loci were diluted and pooled together based on yield, size
range and uorescent dye, and subsequently analysed in
an ABI 310 or ABI 377 automated sequencer. Resulting
patterns were scored and analysed using the ABI computer
programs

genescan

2.1 and

genotyper

2.1.

Data analysis

mtDNA control region sequences were aligned using

clustalx

(Thompson

et al

. 1997) and visually checked. Initial
sequence comparisons and measures of variability were
performed using

mega

(Kumar

et al

. 1993). Transition/
transversion ratios and the

parameter of the gamma
distribution of rate variation among sites (method of Yang
& Kumar 1996) were estimated using

pamp

(included in
the package

paml

2.0

Yang 1999). Phylogenetic analyses
of the identied sequences were performed using three
approaches: (i) minimum evolution (ME) heuristic search
as implemented in

paup

*4.0b2 (Swofford 1998), consisting
of a distance-based (TamuraNei model with gamma
correction) neighbour-joining (Saitou & Nei 1987) tree
followed by branch-swapping; (ii) maximum parsimony
(MP) using

paup

*, with a heuristic search and random
addition of taxa; and (iii) maximum likelihood (ML) using

paup

* and

puzzle

4.0 (Strimmer & Von Haeseler 1996),
incorporating a gamma-corrected HKY (HasegawaKishino
Yano) model with parameters estimated from the data set.
Reliability of nodes dened by the phylogenetic trees was
assessed using 100 bootstrap replications (Felsenstein
1985; Hillis & Bull 1993) in the ME and MP analyses,
and with the quartet puzzling method in the

puzzle

ML
analysis. Minimum-spanning networks (Excofer & Smouse
1994) were constructed using

minspnet

(by L. Excofer,
http://anthropologie.unige.ch/LGB/software/win/) to
depict phylogenetic, geographical, and potential ancestor
descendent relationships among the identied sequences.
Measures of population genetic parameters such as
gene diversity (the probability that two randomly chosen
mtDNA sequences are different in the sample) and nucle-
otide diversity (

per nucleotide site, i.e. the probability
that two randomly chosen homologous nucleotides are
different in the sample) (Nei 1987) were estimated from
the mtDNA data set using

arlequin

1.1 (Schneider

et al

.
1997) and

sendbs

(by N. Takezaki, methods described in
Nei & Jin 1989; Jin & Nei 1990; http://www.bio.psu.edu/
People/Faculty/Nei/Lab/Programs.html.

sendbs

was
applied to derive measures of nucleotide diversity (

) and
their condence intervals (using a bootstrap approach),
which were used to date the origin of extant mtDNA line-
ages of jaguars. Dating of the origin of jaguar mtDNA
haplotypes was also performed using the linearized tree
method with the program

lintree

(Takezaki

et al

. 1995),
that incorporates two tests for the assumption of a molecular
clock. Inference of past population expansion events was
performed using mismatch distribution analyses (Rogers
& Harpending 1992) as implemented in

arlequin

, and
lineages-through-time plots (Nee

et al

. 1994) as imple-
mented in

end-epi

(Rambaut

et al

. 1997). To assess the
extent of differentiation among populations an Analysis
of Molecular Variance (

amova

) (Excofer

et al

. 1992) was
used to estimate

F

ST

values, whose statistical signicance
was tested using 10 000 permutations as implemented in

arlequin

. We also applied the nested cladistic analysis
proposed by Templeton

et al

. (1995) to make phylogeo-
graphic inferences with our mtDNA data set, using

geodis

2.0 (Posada

et al

. 2000) and 10 000 permutations to test the
signicance of alternative historical scenarios.


CONSERVATI ON GENETI CS OF J AGUARS

69


Molecular Ecology

, 10, 6579

Microsatellite data was analysed with

microsat

(Minch
1997) for general estimates of diversity and calculation
of several measures of genetic distance among indi-
viduals. Resulting distance matrices were used to con-
struct neighbour-joining phylogenetic trees with the
program

neighbor

[included in the package

phylip

3.5.
(Felsenstein 1993)].

arlequin

1.1 was used to test for
deviations from HardyWeinberg equilibrium (using an
exact test based on the procedure described by Guo &
Thompson 1992) and to derive estimates of population
subdivision (using

F

ST

and

R

ST

analogs; Weir & Cockerham
1984; Slatkin 1995; Michalakis & Excofer 1996). The
statistical signicance of

F

ST

and

R

ST

values was tested
using 10 000 permutations as implemented in

arlequin

.
To evaluate if there is support from microsatellite data to
the inference of a population expansion in jaguars, the
shape of the distribution of allele frequencies in each
locus was inspected, and two statistical approaches
were applied: (i) the P

k

distribution method of pairwise
comparisons suggested by Shriver

et al

. (1997); and
(ii) the interlocus

g

-test proposed by Reich & Goldstein
(1998) and Goldstein

et al

. (1999). The former is based
on the shape of the distribution of pairwise differences in
repeat number among all alleles of each locus, averaged
across loci. We have performed a weighted average in
which loci with the most complete data sets are pro-
portionally more represented. The second approach (

g

-test)
compares the observed variance across loci of the vari-
ance in repeat scores to that expected in a population of
constant size. We obtained the observed variance of
the variance (V

l

[V

r

]) empirically, based on the output
of

microsat

(using the total variance in repeat number
per locus as V

r

). The expected value under constant
population size (V

g

[V

r

] ) was derived using the method
suggested by Reich & Goldstein (1998) and Goldstein

et al

. (1999), with the formula V

g

[V

r

] = 4/3 E[V

r

]

2

+ 1/6 E[V

r

],
in which the average variance across loci is substituted
for E[Vr].

Results

mtDNA sequences

mtDNA control region PCR products from 37 jaguar
individuals (Table 1) were sequenced in three separate
units. Fragment A was sequenced with primer L15997, and
consisted of 183 bp at the 5

ank of the repetitive array.
Fragment B was sequenced with primer H16498 and
contained

240 bp at the 3

ank of the repetitive region,
including the rst (3

side) 80- bp repeat (all repeats but
the rst were excluded from the data set). Fragment C
spanned the central conserved portion of the mtDNA
control region, containing

390 bp sequenced for both
strands using primers PAN-CCR-F and PAN-CCR-R, and
overlapping with

100 bp of fragment B. Fragment A
sequences were concatenated with the contigs generated
from fragments B and C, producing a complete aligned
data set of 715 bp [positions 1631500218 in the domestic
cat reference sequence (Lopez

et al

. 1996) ]. Sequences
of the same control region segment from a lion (

Panthera
leo

) and a leopard (

P. pardus

) were used as outgroups
( J.-H. Kim

et al

. unpublished data). A 66- bp segment of
the leopard sequence on the 5

ank of the repetitive
array could not be aligned with the other sequences,
and was removed from the analyses. The alignment of
the jaguar sequences (GenBank accession nos AF244814
AF244887) is available at http://lgd.nci.nih.gov (as supple-
mentary information linked to this paper).
Sequences from the 37 jaguar individuals comprised
22 different haplotypes dened by 33 variable sites
(Table 2). The estimated transition:transversion ratios
varied from 1:6 to 1:16 depending on the inclusion of
outgroup sequences. Among-site rate variation was high
(

= 0.00230.3945 depending on inclusion of outgroups).
Moderate to high levels of gene diversity and considerably
low nucleotide diversity were observed among jaguar
individuals (Table 3). Sequence identity was observed only
among individuals originating from the same or adjacent
geographical regions, with the exception of one haplotype
which was widespread in the northern portion of the
range, from Venezuela to Mexico (Pon21, 25, 31, 51, 58, 59,
60, and 75, Figs 2 and 3). Phylogenetic analyses of jaguar
mtDNA sequences did not strongly support the existence of
any major phylogeographic partitions (Fig. 2). The overall
pattern approximated a star-phylogeny, with generally short
branches and little robust structure, a pattern reminiscent
of a relatively recent population expansion (Lavery

et al

.
1996; Avise 2000). The only phylogeographic partition which
received consistent support in all analyses was the grouping
of southern South American sequences (all the sequences
from Brazil and Paraguay, one sequence probably from
Peru, and most from Bolivia) in two monophyletic groups
(Fig. 2). The rooting of the trees using other

Panthera

genus
species as outgroups suggested that the northern South
American/Central American sequences were more basal
than the two southern South American clusters.
The minimum-spanning network (Fig. 3) of the jaguar
mtDNA sequences provided additional support for the
two southern South American clusters, and for their inde-
pendent connection to a more basal northern sequence.
This diagram is also indicative of a relatively recent popu-
lation expansion in the group of northern sequences, in
which several localized lineages are connected by short
branches to the most common, potentially ancestral,
widespread haplotype. With the exception of the basal
Pon37, all Mexican individuals either carried the widespread
Pon21 haplotype or a sequence which had been recently
derived (being one step away) from it.


70

E. EI ZI RI K

ET AL

.


Molecular Ecology

, 10, 6579
Table 2 Mitochondrial DNA control region haplotypes identied for jaguars (Panthera onca) in this study. Nucleotide positions showing
variation among jaguar haplotypes are depicted; numbers (vertical) refer to the aligned site in our 715 bp data set [which corresponds to
positions 1631500218 in the domestic cat reference sequence (Lopez et al. 1996)]. Nucleotides present at these sites in the outgroups
(Panthera leo and P. pardus) are shown in the bottom. The geographical origin of jaguar individuals bearing each haplotype is also
indicated
Nucleotide positions
111112 2222222334 4444445666 677
2356335783 4445689280 0135783137 801
Individuals 4231497368 1582601077 9182658185 300 Geographic origin
Pon21, 25, 31, 51, 58, 59, 60, 75 ATTTAGCTCA AAACTCTGTA TATATGGTGA TAT V, CR, Nic, Mex
Pon16, 26, 30 .......C.. .......... .......... ... Guat, Mexico
Pon52 .......... .......... .......... .G. Mexico
Pon37 .......... .....?.... ?G...A..A. ... Mexico
Pon20 .......... .........G .....A..A. ... Panama
Pon22, 24 .....A.... ........C. ........A. ... CR, Nic
Pon23 .....A...T G?TG?...C. ........A. ... Nicaragua
Pon56 ........T. .......... .......... C.. Venezuela
Pon61 .......... .......A.. .........G ... Venezuela
Pon63 .CC?....T. .....TC... ...G....A. ... Venezuela
Pon68 .......... .......... .........G ... Fr. Guiana
Pon29, 41, 42, 44 G...G.T.T. .G........ C?...A..A. ... Brazil
Pon40, 46 G...G.T.T. .G........ .....A..A. ... Brazil
Pon71 G..CG.T... .G.?C..... ...G....A. ... Brazil
Pon73 G...G.T.T. .G........ ...G....A. ... Brazil
Pon72 ....G...T. .......A.. ??...A..A. ... Brazil
Pon28, 47 ....G..... .......A.. ...GC.A..G .G. Peru (P), Bra
Pon43 ....G..... .......A.. ..CGC.A..G ... Brazil
Pon64 ....G..... .......... ...GC.A..G .GC Bolivia
Pon65 ...?...... .......... .......C.. ... Bolivia
Pon66 ....G..... .......A.. .G.GC.A..G .G. Bolivia
Pon49 ?...G..... .......A.. ...G?.A..G ... Paraguay
Ple181 .CC...T... ....C-CA.. ...G.A.... ...
Ppa21 G.C.----.. .......... .G.G.A.CA. ...
Abbreviations: V, Venezuela; CR, Costa Rica; Nic, Nicaragua; Mex, Mexico; Guat, Guatemala; Fr. Guiana, French Guiana, Bra, Brazil; (P),
Probable origin.
Table 3 Measures of mitochondrial DNA diversity observed in the two major phylogeographic groups of Panthera onca identied in this
study, Southern and Northern (see text for denition), as well as separate estimates for Central America (C.A.) and Northern South
America (N.S.A.), which are included in the Northern group. Estimates for the two separate subgroups in Central America (four-
population scenario) are not provided due the small sample size available for each of them
n
No. of
haplotypes
No. of bp
differences
% Pairwise
divergence* Gene diversity
Nucleotide
diversity ()*
Total 37 22 115 0.142.35 0.9399 0.0262 0.00771 0.00010
Southern 16 11 113 0.141.94 0.9333 0.0477 0.00815 0.00000
Northern 21 11 112 0.141.87 0.8476 0.0697 0.00448 0.00015
C.A. 14 7 18 0.141.25 0.8462 0.0742 0.00411 0.00024
N.S.A. 7 5 19 0.141.29 0.8571 0.1371 0.00475 0.00000
*Estimated using Kimura 2-parameter distances (Kimura 1980).
Sites with any missing information were completely excluded.
CONSERVATI ON GENETI CS OF J AGUARS 71
2001 Blackwell Science Ltd, Molecular Ecology, 10, 6579
The mismatch distribution of jaguar mtDNA haplotypes
(Fig. 4A) was approximately unimodal and wave-like,
which is compatible with a population bottleneck and/or
expansion in the recent past (Rogers & Harpending 1992).
The lineages-through-time plots (Fig. 4B) produced a
clearly convex-shaped curve, particularly when restricted
to the earliest haplotype coalescence events, which seems
to be also consistent with an initial period of exponential
population growth (Nee et al. 1994).
The measures of nucleotide diversity () and divergence
(d
xy
) (Nei 1987) were used to estimate a coalescence date
for jaguar mtDNA sequences (derived as the average
pairwise coalescent time of the sampled jaguar haplo-
types) and its 95% condence interval (CI: 2 SE). For
this analysis, all sites which contained gaps or missing
information in any haplotype, including outgroups,
were excluded, and all estimates were derived from the
remaining 556 nucleotides. This minimized the effects of
saturation in interspecic comparisons as some of the
most variable sites were excluded. To estimate the
substitution rate of this mtDNA segment, the average
distance (d
xy
) between the jaguar sequences and the
outgroups (lion and leopard) was computed (d
xy
= 0.07499
0.002839), based on a Kimura 2-parameter model of
nucleotide substitution, corrected with a gamma distri-
bution of rate variation among sites, with an (shape)
parameter derived specically for this data set ( = 0.13).
The standard equation d
xy
= 2T (where is the substitu-
tion rate and T is time; Nei 1987) and a conservative
fossil record calibration (3 Ma) (Turner & Antn 1997)
for the separation of jaguars from lions/leopards pro-
duced a substitution rate of 0.012498 per site per million
years (Myr), or approximately 1.2%/Myr. Applying the
same equation and this substitution rate to the overall
Fig. 2 Phylogenetic relationships among the jaguar mitochondrial DNA (mtDNA) sequences identied in this study, based on 715 bp
of the mtDNA control region. Labels are individual identication numbers. The two southern clades are identied by brackets.
(A) Maximum parsimony tree constructed with paup* using a heuristic search with 10 replicates of random taxon addition and the
accelerated transformation option for character-state optimization; shown is the strict consensus of eight trees (length = 122, CI = 0.818).
Values above selected branches indicate number of steps/number of homoplasies; values below branches represent support from 100
bootstrap replications (heuristic searches with random taxon-addition, constrained to generate a maximum of 1000 trees); black circles
indicate other branches with bootstrap support 50%. Asterisks below branches indicate support from 1000 quartet puzzling steps in the
maximum likelihood analysis with puzzle (* = support for branch >50%; ** = support for branch >80%). (B) Minimum evolution tree,
constructed with paup* using TamuraNei distances corrected for the gamma distribution of rate variation among sites ( = 0.39) and the
neighbour-joining algorithm followed by a branch-swapping procedure. Values above selected branches indicate branch lengths as
percentages; values below branches represent support from 100 bootstrap replications (heuristic searches constrained to a maximum of
500 trees); black circles indicate other branches with bootstrap support 50%.
72 E. EI ZI RI K ET AL.
jaguar nucleotide diversity in this segment ( = 0.00699
0.000733) produced an estimated age of 279 645 years (95%
CI = 220 990338 290 years) for the coalescence of jaguar
mtDNA haplotypes. If the 95% CI of d
xy
is also considered
when estimating the substitution rate (Bonatto & Salzano
1997), an upper bound of 366 000 years ago is obtained.
Using the 95% CI for both and the substitution rate,
and a more recent calibration date for the separation of
jaguars from lions/leopards (2 Myr), provided an overall
lower bound of 136 956 years ago for the coalescence of
jaguar mtDNA sequences. The molecular clock tests
implemented in lintree showed that the set of jaguar
sequences used for the analysis above did not depart
signicantly from rate constancy (P > 0.05), but indicated
that three haplotypes exhibited a marginally signicant
reduction in rate. Removal of haplotypes that failed either
of the rate-constancy tests in separate analyses with each
of the two outgroups produced a range of estimates of
407 905510 401 years ago for the basal divergence of
jaguar mtDNA lineages (applying the substitution rate
derived above with the 3 Myr calibration). The 95% CI of
these estimates ranged from 165 627 to 830 293 years ago.
Fig. 3 Minimum-spanning network depicting phylogenetic and geographical relationships among jaguar mitochondrial DNA (mtDNA)
haplotypes (circles) based on sequences of 715 bp of the control region. The area of the circle is roughly proportional to the haplotype
frequency, and the length of connecting lines is roughly proportional to the corrected distance (TamuraNei gamma, = 0.39) between
haplotypes. Use of the raw number of nucleotide differences among haplotypes produced a concordant result. Labels inside each circle
indicate individual jaguars bearing each haplotype. Dark grey circles (with no upper horizontal line) are haplotypes occurring only in
Mexico and Guatemala; dark grey circles with an upper horizontal line occurred only in southern Central America (see text for
denition); light grey circles are those found in northern South America (Venezuela and French Guiana); the large graded grey circle is a
widespread haplotype found from Mexico to Venezuela; white dotted circles are haplotypes found only in southern South America
(Brazil, Bolivia, Paraguay and potentially Peru). Bars placed on connecting lines indicate the exact number of nucleotide differences
between haplotypes. White circles bordered by heavy lines are outgroups; the number of nucleotide differences in their connecting lines
(not drawn to scale) is indicated in brackets. Polygons demarcate selected clades for which signicant results were observed in the nested
cladistic analyses (no outgroup was used in this approach). The dotted line internal polygon indicates a geographically restricted 3-step
clade that showed a signicantly small (P < 0.001) within-clade distance, and a signicantly large (P = 0.0001) nested-clade distance (see
Templeton 1998 for denition of statistics). The solid line polygon demarcates a 4-step clade that showed a signicant value (P < 0.0001)
for nonrandom geographical distribution of subclades in the exact contingency test implemented in geodis, and a signicantly large
(P < 0.0001) interior tip distance. This 4-step clade also showed signicantly small within-clade and nested-clade distances (P < 0.05)
when compared to the whole network.
To investigate the extent of population differentiation
in jaguars based on mtDNA haplotypes, we compared
different scenarios of geographical subdivision. In the rst
one, individuals were divided into two broad geographical
groups, which had been suggested by the topology of the
mtDNA phylogeny. One group encompassed the northern
part of the jaguar distribution (Northern group: Mexico,
Central America, Venezuela and French Guiana), and the
other contained the southern part (Southern group: Brazil,
Bolivia, Paraguay, and two samples probably from Peru and
southern South America, respectively). This subdivision
scheme implied that the Southern group contained all
sequences present in the two independent southern clades
(Figs 2 and 3), as well as the haplotype Pon65. In the second
hypothesized scenario of geographical subdivision, we
divided jaguar haplotypes into three groups. The Southern
group was maintained, but the Northern group was
separated into Central America (including Mexico) and
northern South America, aiming to test the signicance of
the Darien straits in Panama as a historical barrier to gene
ow in this species. Finally, a third scenario was tested in
which Central America was subdivided into two separate
groups (Mexico + Guatemala; Nicaragua + Costa Rica +
Panama), producing a total of four groups. This pattern of
geographical subdivision in Central America has been
observed in a previous study of Neotropical felids (Eizirik
et al. 1998). Other possible scenarios of population sub-
division were considered, but could not be thoroughly
assessed at the present time due to limitations in sample
size for particular areas and lack of precision in the
geographical origin of most individuals.
Using the two-group scheme, the estimated F
ST
value was
0.342 (P < 0.0001) (Table 4). Using the three- and four-group
scenarios, the F
ST
estimates decreased to 0.295 and 0.300,
respectively (P < 0.0001 for both). Pairwise F
ST
estimates
between the southern group (south of the Amazon River)
and all northern groups remained signicant, but the F
ST
between northern South America and part of Central
America was not (P > 0.05) (see Table 4). This is likely due
Fig. 4 (A) Mismatch distribution analysis of jaguar mitochondrial
DNA (mtDNA) control region sequences. (B) Lineages-through-time
plot based on a phylogeny of jaguar mtDNA haplotypes constructed
from Kimura 2-parameter distances. The number of lineages is
shown as a semi-logarithmic transform.
Table 4 Measures of geographical population differentiation in
Panthera onca, based on an Analysis of Molecular Variance approach
with mitochondrial DNA (mtDNA) and microsatellite data
Subdivision
mtDNA
F
ST
Microsatellites
F
ST
R
ST
2 populations 0.342** 0.045** 0.154**

3 populations 0.295** 0.065** 0.171**
S. vs. N.S.A. 0.261* 0.020* 0.067*
S. vs. C.A. 0.354** 0.083** 0.237**
N.S.A. vs. C.A. 0.079 0.077** 0.140**
4 populations 0.300** 0.063** 0.183**
S. vs. N.S.A. 0.261* 0.020* 0.067*
S. vs. S.C.A. 0.333* 0.085** 0.223**
S. vs. M + G 0.358** 0.084** 0.266**
N.S.A. vs. S.C.A. 0.185* 0.075* 0.131*
N.S.A. vs. M + G 0.082 0.083** 0.185**
S.C.A. vs. M + G 0.257* 0.025 0.140*
*signicant (P < 0.05); **signicant (P < 0.001).
Population subdivision scenarios are described in the text; the
two-population scheme depicts south (S) vs. north (N)
comparisons; the three-population scheme retains the southern
group (S), and divides the northern group into northern South
America (N.S.A.) and Central America (C.A.); the four-
population scheme further divides C.A. into separate S.C.A.
(southern Central America: Panama, Costa Rica, Nicaragua) and
M + G (Mexico + Guatemala).
Calculated using Tamura-Nei distances (Tamura & Nei 1993).
Calculated with the F
ST
option (number of different alleles) in
arlequin.
Calculated with the R
ST
option (sum of squared size difference)
in arlequin (Michalakis & Excofer 1996).
to the predominance in Mexico and Guatemala of the
widespread Pon21 haplotype and others closely related
to it. These mtDNA results support the recognition of two
main phylogeographic groups of jaguars, Southern and
Northern (probably separated by the Amazon River),
but also suggest a pattern of regional subdivision in the
northern part of the range. F
ST
estimates based on gamma-
corrected distances with extremely low values of (not
shown) also supported a southnorth division of jaguar
haplotypes.
The nested cladistic analysis proposed by Templeton
et al. (1995) was performed on the basis of the haplotype
tree shown in Fig. 3, and considering four broad geo-
graphical areas (four-population scenario described above)
as operational units. Signicant association of clades and
geography was observed only at the 3-step and 4-step
levels of the nested analysis. In one of the two 4-step
clades (demarcated within a solid-line polygon in Fig. 3)
the exact contingency test signicantly rejected (P < 0.0001)
the null hypothesis of no geographical association. In this
clade, signicantly small (P < 0.05) within-clade distances
(D
c
) were detected for all nested tip 3-step subclades, and
a signicantly large (P < 0.0001) interior tip D
c
was also
observed. One of these nested 3-step subclades (shown as
a dotted-line polygon in Fig. 3) exhibited a signicant
reversal of small within-clade distance vs. large nested-
clade distances (D
n
), which led to the inference of long-
distance colonization (from north to south) as the cause
for the observed pattern (see Templeton 1998 for a review
on the interpretation of this analysis). When the whole
tree was examined (a dichotomous comparison of the
indicated 4-step clade vs. the remaining jaguar haplo-
types), the null hypothesis was again rejected (P < 0.05)
in the exact contingency test, signicantly small (P < 0.05)
D
c
and D
n
were detected for the demarcated clade, and
a signicantly large D
n
was detected for the second major
clade. These results suggest a historical scenario of
restricted gene ow with isolation by distance over broad
geographical areas.
Microsatellite variation
Allele size variation at 29 microsatellite loci was recorded
for 42 jaguars (Table 1), as well as for two lions and two
leopards to be used as outgroups. Twenty-seven of these
loci were dinucleotide repeats, and two (FCA441 and
FCA453) were tetranucleotide repeats (Menotti-Raymond
et al. 1999). One of the dinucleotide loci (FCA293) was
monomorphic in our jaguar sample. The genotypes
used in the present analyses are available at http://
lgd.nci.nih.gov (as supplementary information linked to
this paper). All loci, except FCA098, have been mapped in
the domestic cat (Menotti-Raymond et al. 1999). Twenty-
four were at least 15 cm apart implying independent
allele assortment; two pairs of loci [FCA075, FCA096] and
[FCA161, FCA224] were linked at distances of 9.0 cm and
4.0 cm, respectively. Deviations from HardyWeinberg
(HW) equilibrium in the microsatellite loci were tested
separately for three geographical populations: southern
South America, northern South America, and Central
America (see below). Five loci were out of HW equilibrium
in the former, six in the second population, and nine in
the latter ( = 0.05). After applying the sequential Bonferroni
correction (Rice 1989), only one locus in the southern
group and another in Central America remained signicantly
out of HW equilibrium. Comparative analyses excluding
the 14 loci that were out of equilibrium in any population
(uncorrected assessment) indicated no signicant change
in the results relative to our complete data set, other than
slight reductions in estimates of genetic diversity and
geographical differentiation. We have, therefore, used our
complete data set of 29 loci throughout the analyses
presented here, except where otherwise stated.
Overall levels of genetic diversity were moderate to high,
with an average expected heterozygosity of 0.739 (Table 5).
The phylogenetic analysis of microsatellite variation across
individuals showed no evidence of strong geographical
structure, nor provided reliable support for any major
partition, but suggested a clustering of individuals from
Table 5 Measures of diversity at 29 microsatellite loci in the two major phylogeographic groups of Panthera onca identied in this study,
Southern and Northern (see text for denition), as well as separate estimates for Central America (C.A.) and Northern South America
(N.S.A.), which are included in the Northern group. Estimates for the two separate subgroups in Central America (four-population
scenario) are not provided due the small sample size available for each of them
n
Average
heterozygosity
Average
No. of alleles
Average
variance
Average allele
size range
No. of private
alleles
Total 42 0.739 8.310 6.946 9.862
Southern 17 0.724 6.655 6.520 8.034 43
Northern 25 0.695 6.828 6.352 8.690 48
C.A. 16 0.622 5.172 6.118 7.310 22
N.S.A. 9 0.690 5.138 6.514 6.931 12
nearby regions (Fig. 5). Phylogenetic trees generated
using different distance measures (Dps, Dkf Bowcock
et al. 1994) led to very similar conclusions.
The shape of the distribution of allele sizes varied across
microsatellite loci. Of the 28 variable loci, 13 had clearly
unimodal distributions, which according to Reich & Goldstein
(1998) can be expected for populations that underwent a
demographic expansion in the past. Five other loci had a
distribution that was approximately unimodal, while 10
loci clearly departed from this pattern. The two statistical
approaches to estimate demographic history with micro-
satellite data were performed for the jaguar sample as a
whole, using only the 26 polymorphic dinucleotide loci,
and provided results that are also compatible with a
population expansion. The P
k
method of Shriver et al. (1997)
produced an average distribution of pairwise differences
in repeat number which showed the highest point at a
plateau between values 0 and 1 (Fig. 6). The g-test of
Reich & Goldstein (1998) produced a value of 0.592 for
the ratio V
l
[V
r
]/V
g
[V
r
]. Inclusion of only 24 unlinked loci
in the analysis produced a similar result.
As with the mtDNA data, the statistical analysis of
population differentiation was performed for different
hypothesized geographical partitions. Measures of differ-
entiation (F
ST
and R
ST
analogs) were signicant (P < 0.05)
in all comparisons except for the subdivision of Central
America into two separate groups (four-population
scenario) in the F
ST
estimate (Table 4). R
ST
estimates, which
correct for allele size variation, indicated higher levels of
differentiation than F
ST
estimates, which only consider
the number of different alleles. However, both of these
methods produced higher overall levels of differentiation
with three geographical groups compared with two, and
the R
ST
estimate was maximized when four separate groups
were hypothesized (see Table 4). These microsatellite
results support the recognition of at least three or four
incompletely isolated phylogeographic groups in jaguars,
but may also reect, at least in part, a pattern of isolation
by distance over most of the species range. Measures of
microsatellite diversity were similar in different geographical
regions (Table 5).
Discussion
Population history
Observed levels of genetic diversity in jaguars were low
to moderate for the mtDNA data, and medium to high in
terms of microsatellite variability. The mtDNA diversity
estimates presented here cannot be directly compared to
other felids because data on the equivalent segment is
unavailable from other species. However, if only the
453-bp long segment used by Eizirik et al. (1998) in two
smaller Neotropical cats, ocelots (Leopardus pardalis) and
margays (L. wiedii), is considered, jaguars show a three-
to tenfold reduction in several measures of diversity
(e.g. per cent sequence divergence among haplotypes,
nucleotide diversity) relative to these two species (Eizirik
Fig. 5 Phylogenetic relationships among jaguar individuals based
on variation at 29 microsatellite loci. The tree was constructed
using the Neighbour-Joining algorithm based on the Proportion
of Shared Alleles distance measure (Bowcock et al. 1994), using
the option (1-ps) in microsat. Values above branches are percent
support in 100 bootstrap replications (only values above 50% are
shown). The geographical origin of individuals is indicated next
to each sample label.
Fig. 6 P
k
distribution of differences in repeat number among
alleles at a microsatellite locus, averaged in a weighted fashion
across 26 dinucleotide loci in our complete Panthera onca data set.
et al. 1998; this study). By contrast, the microsatellite
variation observed in jaguars is comparable to, or higher
than, that of other felid species analysed so far, including
the domestic cat, for which these markers were originally
designed (Driscoll 1998; Johnson et al. 1999; Menotti-
Raymond et al. 1999; Culver et al. 2000). Taken together,
these results suggest a younger age for jaguars as a species
compared to other Neotropical felids, particularly ocelots
and margays (Eizirik et al. 1998; Johnson et al. 1998, 1999).
The pattern observed in the jaguar intraspecic mtDNA
phylogeny is consistent with a fairly recent population
expansion, with little divergence among haplotypes (Table 3,
Fig. 2) and insufcient time and/or isolation to generate
regional differentiation. A similar inference can be made
from the microsatellite phylogeny (Fig. 5). More individual
variation is observed, but there is still no strong support
for major geographical structure, suggesting a recent origin
and high levels of gene ow in this species. The inference
of a population expansion in the relatively recent past is
supported by the mtDNA mismatch distribution analysis
and the lineages-through-time-plot (Fig. 4). The mismatch
distribution has a prominent major peak suggestive of a
population expansion event (Rogers & Harpending 1992).
However, the curve is not perfectly smooth, which could
be due to our small sample size or to a more complex
population history (e.g. subdivision) after the inferred
demographic growth. The convex shape of the lineages-
through-time plot is also consistent with exponential
population growth (Nee et al. 1994). But it also portrays
a more recent period when this pattern is less visible,
indicating that the demographic expansion may have
been restricted to an earlier phase.
The microsatellite allele frequency distributions are
also concordant with the inference of an expansion due to
the fact that more than half of the variable loci show a
unimodal or approximately unimodal pattern. The absence
of such indications in 10 of the loci may reect our small
sample size or alternatively derive from demographic
processes occurring after the hypothesized expansion.
One of the 10 loci (FCA026) has a bimodal distribution
that can be largely explained by dividing the jaguar sample
in northsouth geographical groups (not shown).
The P
k
distribution of our jaguar microsatellite data set
is similar to the simulation results of Shriver et al. (1997)
of either a population at equilibrium with an effective
population size (N
e
) around 1000, or a population result-
ing from a 1000-fold expansion 500 generations ago.
Although no strict comparisons should be attempted, it is
almost certain that the N
e
of jaguars as a species was,
until very recently, much larger than 1000. Therefore, this P
k
pattern is suggestive of a population which has still not
reached equilibrium after an expansion event. The inter-
locus g-test result of less than one is also suggestive of
population growth (Goldstein et al. 1999), although it seems
to fall short of signicance based on the critical value sug-
gested by Reich & Goldstein (1998) for a similar data set.
The estimated coalescence date of jaguar mtDNA
haplotypes puts a lower bound on the origin of extant
lineages around 137 000 years ago, which indicates that
the inferred population expansion happened around
or before this time. This estimation suggested an origin
280 000510 000 years ago, with an overall upper bound
of 830 000 years ago, still slightly younger than the con-
rmed fossil record of jaguars in north America (over 850 000
years old; Seymour 1989). The phylogenetic relationships
of jaguar mtDNA sequences (particularly Fig. 3) suggest
that most extant lineages arose in northern South America,
and that the colonization of Central America and Mexico
was largely a subsequent process. This suggests that
the postulated extinction of jaguars from most of North
America towards the end of the Pleistocene left little residual
genetic diversity, and that subsequent re-colonization
occurred by immigration of South American founders.
The occurrence of two independent clades of jaguar
mtDNA sequences occurring south of the Amazon River
does not seem to indicate geographically segregated clusters,
at least with the samples currently available. Rather, these
clades seem to result from a recent colonization process of
southern South America in which at least two unrelated
founder haplotypes were present. This view is supported
by the observation that these two clades have similar
estimates of nucleotide diversity (not shown), and there-
fore can be interpreted as being of a similar age.
Population structure
Analyses of the genetic structure of jaguar populations
showed no evidence of major geographical partitions, old
subdivision events, or complete barriers to historical gene
ow. This is consistent with the ndings of Larson (1997),
that were based on skull morphology, and in contrast with
the traditional division of jaguar subspecies (reviewed by
Seymour 1989). We could recognize up to four incompletely
isolated phylogeographic partitions: Mexico + Guatemala,
southern Central America, northern South America (north
of the Amazon river), and southern South America (south
of the river). These geographical groups are rather recent,
however, and none demonstrated complete reciprocal
monophyly in relation to the others (Figs 2, 3 and 5).
Furthermore, the mtDNA data did not support the
separation of Central America from northern South
America, because the F
ST
value derived from this marker
was maximized when only two major phylogeographic
groups (Southern and Northern) were considered (Table 4).
On the other hand, the microsatellite data indicated
higher levels of differentiation between Central America
and northern South America than between northern and
southern South America. In most cases, estimated levels of
population differentiation were lower with the microsatellite
data than with mtDNA (Table 4), which could be due to
differences in the effective population size assessed by
these two types of markers, and/or in the patterns of
dispersal by male and female jaguars. Estimates of the
number of migrants per generation (N
m
) derived from
mtDNA F
ST
and microsatellite F
ST
and R
ST
values were
high (around or above 1), but cannot be interpreted as a
direct measure of current gene ow because they are also
likely to be inuenced by the recent common origin that
these populations share.
The observed patterns can be interpreted as being
caused by a combination of several demographic phe-
nomena, starting with the recent expansion of the species
as a whole, preventing sufcient time for extensive genetic
differentiation. In addition, jaguars are able to disperse
efciently over broad areas and diverse habitats (Oliveira
1994), which is consistent with our inference of high levels
of gene ow on wide geographical scales. The Amazon river
may represent a historical barrier to gene ow predomin-
antly in females, producing the only visible structure in
the mtDNA phylogenies. However, the Amazon appears
to have been less of an impediment for male dispersal, as
inferred from our microsatellite data (Table 4, Fig. 5). This
is in agreement with known aspects of the dispersal biology
of felids in general, which indicate that females are
largely philopatric and males are more likely to disperse
longer distances (Kitchener 1991; Oliveira 1994). Our data
also provide evidence for more recent barriers to gene ow
between northern South America and southern Central
America, and between Nicaragua and Guatemala, producing
a measurable impact on the frequencies of microsatellite
alleles and mtDNA haplotypes. Similar phylogeographic
partitions that are likely caused by the same historical
barriers have been observed for other Neotropical cat
species, such as the puma (Culver et al. 2000), ocelot and
margay (Eizirik et al. 1998). In addition to the inferred
impact of localized historical barriers to gene ow, an
overall pattern of isolation by distance over broad regions
may also account for the observed patterns, and should
be further investigated with additional sampling and
rened geographical coverage.
Implications for jaguar conservation
These results have implications for the conservation and
management of jaguars in the wild and in captivity.
Genetic variability as estimated with mtDNA sequence
diversity and microsatellite size variation are indicative of
an outbred species with a somewhat recent origin, with
no appreciable evidence for population- or species-level
reduction in genetic variation.
The absence of deep geographical subdivision, the
evidence for recent and incomplete isolation among major
regions, and the inference of historically high levels of gene
ow over broad areas argue for a revision of the traditionally
held view of jaguar subspecies, which are usually the basis
for conservation and management plans. Our data pro-
vide no support for the existence of major geographical
partitions dening old, extensively isolated groups which
could be viewed as subspecies or Evolutionary Signicant
Units (Ryder 1986; Moritz 1994). The only evidence for
statistically signicant geographical differentiation seems
to support reduced levels of gene ow across barriers such
as the Amazon river and the Darien Straits, and isolation
by distance on a broad scale. Although the isolation gen-
erated by these localized barriers does not seem to be
complete, management options which increase migration
across them (e.g. eld translocations or captive breeding
of Amazonian individuals from the two sides of the river)
should be avoided or minimized to approximate natural
historical processes. Our results are compatible with
currently favoured strategies for jaguar conservation on a
broad regional basis (Medelln et al., in press), having large-
scale biomes or ecosystems as operational management
units. This strategy is advisable both in terms of main-
taining jaguars as an important component of functional
ecosystems and as a agship species for overall regional
conservation. This management approach is also appro-
priate considering the possibility of local adaptation to
different habitats, which would probably not be detected
by our molecular markers (Lynch 1996), particularly given
the number of samples analysed.
Finally, because of continued reduction, fragmentation
and isolation of remaining jaguar populations, it is import-
ant to consider the inference presented here for high
historical levels of gene ow over broad areas. Our results
support management alternatives that permit the main-
tenance of this original and likely adaptive population
dynamics, possibly including active intervention to mediate
gene ow among isolated remnant populations.
Acknowledgements
We would like to express our thanks to all the people and Institu-
tions that helped collect samples of biological material from jaguars,
including those mentioned in Table 1, as well as C.B. Indrusiak,
J.K. F. Mhler Jr., D.A. Sana, R.G. Morato, R. Gasparini-Morato,
Companhia Energtica de So Paulo (CESP) and Associao
Pr-Carnvoros. We would like to thank L.R.P. Utz, S.L. Bonatto,
N. Takezaki, C.B. Indrusiak, D.A. Sana, L. Cullen, R.J. Sarno,
C.A. Driscoll, O. Uphyrkina, G. Bar-Gal, V.A. David, S. Cevario and
A. Snyder for support, helpful information, discussions, sugges-
tions and technical assistance. We also thank three anonymous
reviewers for helpful suggestions and comments on a previous
version of this paper. The American Association of Zoological Parks
and Aquariums (AZA) and the Purina Conservation Endowment
Fund (CEF) provided nancial support for this research. E. Eizirik
is supported by a fellowship from CNPq, Brazil. Samples of bio-
logical material were collected with the appropriate permits, and
imported with a permit issued to the National Cancer Institute,
National Institutes of Health (principal ofcer S.J. OBrien).
References
Avise JC (1989) Role of molecular genetics in recognition and
conservation of endangered species. Trends in Ecology and Evo-
lution, 4, 279281.
Avise JC (2000) Phylogeography: the History and Formation of Species.
Harvard University Press, Cambridge.
Bonatto SL, Salzano FM (1997) A single and early migration for the
peopling of the Americas supported by mitochondrial DNA
sequence data. Proceedings of the National Academy of Sciences of
the USA, 94, 18661871.
Boom R, Sol CJA, Salimans MMM et al. (1990) Rapid and simple
method for purication of nucleic acids. Journal of Clinical
Microbiology, 28, 495503.
Bowcock AM, Ruiz-Linares A, Tomfohrde J et al. (1994) High
resolution of human evolutionary trees with polymorphic
microsatellites. Nature, 368, 455457.
Culver M, Johnson WE, Pecon-Slattery J, OBrien SJ (2000) Genomic
ancestry of the American puma (Puma concolor). Journal of Heredity,
91, 186197.
Driscoll CA (1998) A characterization of microsatellite loci variation
in Panthera leo, Acinonyx jubatus and Felis concolor. MSc Thesis.
Hood College, Frederick, MD, USA.
Eizirik E, Bonatto SL, Johnson WE et al. (1998) Phylogeographic
patterns and mitochondrial DNA control region evolution in
two Neotropical cats (Mammalia, Felidae). Journal of Molecular
Evolution, 47, 613624.
Excofer L, Smouse PE (1994) Using allele frequencies and geo-
graphic subdivision to reconstruct gene trees within a species:
molecular variance parsimony. Genetics, 136, 343359.
Excofer L, Smouse P, Quattro J (1992) Analysis of molecular variance
inferred from metric distances among DNA haplotypes: applica-
tion to human mitochondrial DNA restriction data. Genetics,
131, 479491.
Felsenstein J (1985) Condence limits on phylogenies: an approach
using the bootstrap. Evolution, 39, 783791.
Felsenstein J (1993) PHYLIP 3.5. Washington University, Seattle.
Goldstein DB, Roemer GW, Smith DA et al. (1999) The use of micro-
satellite variation to infer population structure and demographic
history in a natural model system. Genetics, 151, 797801.
Guo S, Thompson E (1992) Performing the exact test of Hardy
Weinberg proportion for multiple alleles. Biometrics, 48, 361372.
Hillis DM, Bull JJ (1993) An empirical test of bootstrapping as a
method for assessing condence in phylogenetic analysis.
Systematic Biology, 42, 182192.
Jin L, Nei M (1990) Limitations of the evolutionary parsimony
method of phylogenetic analysis. Molecular Biology and Evolution,
7, 82102 (Published erratum appears in Molecular Biology and
Evolution, 7, 201).
Johnson WE, Culver M, Iriarte A et al. (1998) Tracking the evolution
of the elusive Andean Mountain Cat (Oreailurus jacobita) from
mitochondrial DNA. Journal of Heredity, 89, 227232.
Johnson WE, Dratch PA, Martenson JS, OBrien SJ (1996) Resolu-
tion of recent radiations within three evolutionary lineages of
Felidae using mitochondrial restriction fragment length poly-
morphism variation. Journal of Mammalian Evolution, 3, 97120.
Johnson WE, OBrien SJ (1997) Phylogenetic reconstruction of the
Felidae using 16S rRNA and NADH-5 mitochondrial genes.
Journal of Molecular Evolution, 44 (Suppl. 1), S98S116.
Johnson WE, Pecon-Slattery J, Eizirik E et al. (1999) Disparate
phylogeographic patterns of mitochondrial DNA variation in
four closely related South American small cat species. Molecular
Ecology, 8, S79S94.
Kimura M (1980) A simple method for estimating evolutionary rate
of base substitution through comparative studies of nucleotide
sequences. Journal of Molecular Evolution, 16, 111120.
Kitchener A (1991) The Natural History of the Wild Cats. Comstock/
Cornell, Ithaca.
Kumar S, Tamura K, Nei M (1993) MEGA, Version 1.1. The Pennsy-
lvania State University, University Park, PA.
Larson SE (1997) Taxonomic re-evaluation of the jaguar. Zoo Biology,
16, 107120.
Lavery S, Moritz C, Fielder DR (1996) Genetic patterns suggest
exponential population growth in a declining species. Molecular
Biology and Evolution, 13, 11061113.
Lopez JV, Cevario S, OBrien SJ (1996) Complete nucleotide
sequences of the domestic cat (Felis catus) mitochondrial genome
and a transposed mtDNA tandem repeat (Numt) in the nuclear
genome. Genomics, 33, 229246.
Lynch M (1996) A quantitative-genetic perspective in conservation
issues. In: Conservation Genetics: Case Histories from Nature (eds
Avise JC, Hamrick JL), pp. 471501. Chapman & Hall, New York.
Medelln RA, Chetkiewicz C, Rabinowitz A et al. (eds) (2001)
Jaguars in the new millennium: A status assessment, priority detec-
tion, and recommendations for the conservation of jaguars in the
Americas. Universidad Nacional Autonoma de Mexico/Wildlife
Conservation Society, Mexico D. F (in press).
Menotti-Raymond M, David VA, Lyons LA et al. (1999) A genetic
linkage map of microsatellites in the domestic cat (Felis catus).
Genomics, 57, 923.
Michalakis Y, Excofer L (1996) A generic estimation of popula-
tion subdivision using distances between alleles with special
reference to microsatellite loci. Genetics, 142, 10611064.
Minch E (1997) MICROSAT 1.5d. Stanford University, http://
lotka.stanford.edu/microsat.html.
Montgomery GW, Sise JA (1990) Extraction of DNA from sheep
white blood cells. New Zealand Journal of Agricultural Research,
33, 437441.
Moritz C (1994) Dening Evolutionarily Signicant Units for
conservation. Trends in Ecology and Evolution, 9, 373375.
Nee S, May RM, Harvey PH (1994) The reconstructed evolutionary
process. Philosophical Transactions of the Royal Society (London) B,
344, 305311.
Nei M (1987) Molecular Evolutionary Genetics. Columbia University
Press, New York.
Nei M, Jin L (1989) Variances of the average numbers of nucle-
otide substitutions within and between populations. Molecular
Nowell K, Jackson P (1996) Wild Cats: status survey and conservation
action plan. IUCN/SSC Cat Specialist Group, Gland, Switzerland.
OBrien SJ (1994) A role for molecular genetics in biological con-
servation. Proceedings of the National Academy of Sciences of the USA,
91, 57485755.
Oliveira TG (1994) Neotropical cats: ecology and conservation. EDUFMA,
So Lus, Brazil.
Pbo S, Gifford JA, Wilson AC (1988) Mitochondrial DNA
sequences from a 7000-year-old brain. Nucleic Acids Research,
16, 12451255.
Posada D, Crandall KA, Templeton AR (2000) GeoDis: a program
for the cladistic nested analysis of the geographical distribution
of genetic haplotypes. Molecular Ecology, 9, 487488.
Rambaut A, Harvey PH, Nee S (1997) End-Epi: An application for
inferring phylogenetic and population dynamical processes from
molecular sequences. Computer Applications in the Biosciences,
13, 303306.
Reich DE, Goldstein DB (1998) Genetic evidence for a Paleolithic
human population exansion in Africa. Proceedings of the National
Academy of Sciences of the USA, 95, 81198123.
Rice WR (1989) Analyzing tables of statistical tests. Evolution, 43,
223225.
Rogers AR, Harpending H (1992) Population growth makes waves
in the distribution of pairwise genetic differences. Molecular
Ryder OA (1986) Species conservation and systematics: the dilemma
of subspecies. Trends in Ecology and Evolution, 1, 910.
Saiki RK, Scharf S, Faloona F et al. (1985) Enzymatic amplication
of beta-globin genomic sequences and restriction site analysis
for diagnosis of sickle cell anemia. Science, 230, 13501354.
Saitou N, Nei M (1987) The neighbor-joining method: a new method
for reconstructing phylogenetic trees. Molecular Biology and
Salles LO (1992) Felid phylogenetics: extant taxa and skull
morphology (Felidae, Aeluroidea). American Museum Novitates,
3047, 167.
Sambrook J, Fritsch EF, Maniatis T (1989) Molecular Cloning: a
Laboratory Manual. 2nd edn. Cold Spring Harbor Laboratory
Press, New York.
Schneider S, Kueffer JM, Roessli D, Excofer L (1997) ARLEQUIN,
Version 1.1. University of Geneva, Switzerland. (http://lgb.unige.ch/
arlequin).
Seymour KL (1989) Panthera onca. Mammalian Species, 340, 19.
Shriver MD, Jin L, Ferrel RE, Deka R (1997) Microsatellite data
support an early population expansion in Africa. Genome Research,
7, 586591.
Slatkin M (1995) A measure of population subdivision based on
microsatellite allele frequencies. Genetics, 139, 457462.
Strimmer K, Von Haeseler A (1996) Quartet puzzling: a quartet
maximum likelihood method for reconstructing tree topologies.
Molecular Biology and Evolution, 13, 964969.
Swofford DL (1998) PAUP* Phylogenetic Analysis Using Parsimony
and Other Methods, Computer Program. Sinauer, Sunderland,
MA.
Takezaki N, Rzhetsky A, Nei M (1995) Phylogenetic test of the
molecular clock and linearized trees. Molecular Biology and
Tamura K, Nei M (1993) Estimation of the number of nucleotide
substitutions in the control region of mitochondrial DNA in
humans and chimpanzees. Molecular Biology and Evolution, 10,
512526.
Templeton AR (1998) Nested clade analyses of phylogeographic
data: testing hypotheses about gene ow and population history.
Molecular Ecology, 7, 381397.
Templeton AR, Routman E, Phillips CA (1995) Separating popu-
lation structure from population history: a cladistic analysis of
the geographical distribution of mitochondrial DNA haplotypes
in the tiger salamander, Ambystoma tigrinum. Genetics, 140,
767782.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG
(1997) The ClustalX windows interface: exible strategies for
multiple sequence alignment aided by quality analysis tools.
Nucleic Acids Research, 25, 48764882.
Turner A, Antn M (1997) The Big Cats and Their Fossil Relatives.
Columbia University Press, New York.
Ward RH, Frazier BL, Dew-Jager K, Pbo S (1991) Extensive
mitochondrial diversity within a single Amerindian tribe.
Proceedings of the National Academy of Sciences of the USA, 88,
87208724.
Weir BS, Cockerham CC (1984) Estimating F-statistics for the
analysis of population structure. Evolution, 38, 13581370.
Wozencraft WC (1993) Felidae. In: Mammal Species of the World: a
Taxonomic and Geographic Reference 2nd edn. (eds Wilson DE,
Reeder DM), pp. 288300. Smithsonian Institution Press/
American Society of Mammalogists, Washington.
Yang Z (1999) Phylogenetic Analysis by Maximum Likelihood (PAML),
Version 2.0. University College London, London.
Yang Z, Kumar S (1996) Approximate methods for estimating
the patterns of nucleotide substitution and the variation of
substitution rates among sites. Molecular Biology and Evolution,
13, 650659.
This study was part of a broad effort to characterize the genetic
diversity, phylogenetic relationships and phylogeographic patterns
of Neotropical felid species using molecular markers. Parallel
studies on these species were performed addressing aspects of
their reproductive physiology and exposure to infectious disease.
The authors and several of their colleagues at the Laboratory of
Genomic Diversity (LGD), National Cancer Institute, are studying
the genetics of felids and other mammals to address issues such
as phylogenetic reconstruction, genome evolution, hereditary and
infectious disease, origin of phenotypic variation, and conservation
of endangered species. Eduardo Eizirik is currently developing
his Ph.D. at the LGD and the University of Maryland at College
Park, on the genetics of coat colour variation in felids, which
relates to this study in the context of the evolution of melanism
in jaguar populations.

Phylogeography, Population History and Conservation Genetics of Jaguars (, Mammalia, Felidae)

Uploaded by

Copyright:

Available Formats

Phylogeography, Population History and Conservation Genetics of Jaguars (, Mammalia, Felidae)

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Phylogeography, Population History and Conservation Genetics of Jaguars (, Mammalia, Felidae)

Uploaded by

Copyright:

Available Formats

Molecular Ecology (2001)

Mexico, unknown State Cd.Victoria Zoo, Aragon Zoo

Brazil, unknown State Zoolgico de Goinia

Chaco, Paraguay Itaipu, Paraguay

Amazonas State, Venezuela Las Delicias

, 63 Bolivar State, Venezuela Valencia Castellito

Venezuela, unknown State PROFAUNA

2 populations 0.342 0.045 0.154**

You might also like