Microsatellite data resolve phylogeographic patterns in European grayling, Thymallus thymallus, Salmonidae

The phylogeography of an endangered salmonid, European grayling (Thymallus thymallus), was studied based on analysis of 17 nuclear microsatellite DNA loci. In agreement with earlier mitochondrial DNA (mtDNA) studies, phylogenetic relationships of the populations suggested that northern Europe was colonized from two distinct Pleistocene refugia. Furthermore, microsatellites revealed highly supported grouping of mainland Swedish, Norwegian, Danish, German and Slovenian populations, suggesting that grayling from the northwestern and central Europe have descended from their southern conspecifics. The level of divergence between populations was substantial, even across short geographical distances. Although this was in part due to postglacial colonization patterns and contemporary barriers for gene flow, the high divergence estimates between hydrologically connected sampling locations implied efficient interpopulation reproductive isolation. Microsatellites revealed that the populations exhibited, on average, only 3.5 (+/-2.2) alleles per locus, indicating that T. thymallus has strikingly low levels of intrapopulation genetic diversity as compared with other freshwater fish species. Accordingly, as indicated by analysis of molecular variance (AMOVA), only 49.1-58.0% of the total grayling microsatellite diversity resided within populations. A latitudinal genetic diversity gradient, potentially resulting from glaciation-mediated founder events, was not evident. Alternatively, it is possible that grayling display limited dispersal behaviour/capability, leading to low long-term effective population sizes and, consequently, depauperate intrapopulation polymorphism. These findings have implications for conservation of T. thymallus. Importantly, they exemplify that microsatellites can be highly informative for intraspecific phylogeography studies dealing with substantial divergence scales.     

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.     

Historical introgression and the persistence of ghost alleles in the intermediate horseshoe bat (Rhinolophus affinis)

Phylogenetic conflicts between genetic markers can help to disentangle complex histories of phylogeography and introgression among taxa. We previously proposed that the Chinese mainland subspecies of the intermediate horseshoe bat Rhinolophus affinis himalayanus colonized Hainan Island to form the subspecies R. a. hainanus. Subsequent recolonization of the mainland formed a third taxon, R. a macrurus, and a secondary contact zone with the ancestral himalayanus. To test for historical and recurrent genetic exchange between these mainland subspecies, we sampled populations of each from two parapatric zones and undertook analyses using one mtDNA marker, three nuclear genes and 14 microsatellites. Nuclear DNA, echolocation call frequencies and morphological data all recovered two taxa; however, a mtDNA phylogeny revealed two himalayanus clades, of which one clustered with macrurus, as well as some shared or related mtDNA haplotypes in eastern populations. Isolation‐with‐migration (IM) models suggested some mtDNA gene flow from macrurus to himalayanus. However, strong population structure in himalayanus raises the possibility that macrurus captured mtDNA from a coastal population of himalayanus that has since become rare or extinct. To reconcile these two sets of results, we suggest that the IM estimates might reflect historical mtDNA gene flow among populations of himalayanus, before mtDNA was subsequently captured by macrurus. Finally, microsatellite‐based ABC analyses supported the island origin of macrurus; however, mtDNA‐based ABC analyses suggest this taxon might have evolved on the mainland. Our findings highlight the importance of understanding population history and structure for interpreting hybridization and introgression events.     

Population Structure of the Giant Garter Snake, Thamnophis gigas

The giant garter snake, Thamnophis gigas, is a threatened species endemic to California’s Central Valley. We tested the hypothesis that current watershed boundaries have caused genetic differentiation among populations of T. gigas. We sampled 14 populations throughout the current geographic range of T. gigas and amplified 859 bp from the mitochondrial gene ND4 and one nuclear microsatellite locus. DNA sequence variation from the mitochondrial gene indicates there is some genetic structuring of the populations, with high F_ST values and unique haplotypes occurring at high frequency in several populations. We found that clustering populations by watershed boundary results in significant between-region genetic variance for mtDNA. However, analysis of allele frequencies at the microsatellite locus NSU3 reveals very low F_ST values and little between-region variation in allele frequencies. The discordance found between mitochondrial and microsatellite data may be explained by aspects of molecular evolution and/or T. gigas life history characteristics. Differences in effective population size between mitochondrial and nuclear DNA, or male-biased gene flow, result in a lower migration rate of mitochondrial haplotypes relative to nuclear alleles. However, we cannot exclude homoplasy as one explanation for homogeneity found for the single microsatellite locus. The mitochondrial nucleotide sequence data supports conservation practices that identify separate management units for T. gigas.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

Regional differentiation among populations of the Diamondback terrapin (Malaclemys terrapin)

The Diamondback terrapin (Malaclemys terrapin) is a brackish-water turtle species whose populations have been fragmented due to anthropogenic activity such as development of coastal habitat and entrapment in commercial blue crab (Callinectes sapidus) fishing gear. Genetic analyses can improve conservation efforts for the long-term protection of the species. We used microsatellite DNA analysis to investigate levels of gene flow among and genetic variability within 21 geographically separate collections of the species distributed from Massachusetts to Texas. Quantified levels of genetic variability (allelic diversity, genotypic frequencies, and heterozygosity) revealed three zones of genetic discontinuity, resulting in four discrete populations: Northeast Atlantic, Coastal Mid-Atlantic, Florida and Texas/Louisiana. The average number of alleles and expected heterozygosity for the four genetic clusters were N_A = 6.54 and H_E = 0.050, respectively. However, the geographic boundaries of the populations did not correspond to accepted terrapin subspecies limits. Our results illuminate not only the need to sample terrapins in additional sites, specifically in the southeast, but also the necessity for allowing uninterrupted gene flow among population groupings to preserve current levels of genetic diversity.     

Polymorphic microsatellite loci for tiger salamanders,Ambystoma tigrinum

Microsatellites have been developed for few amphibian species. However, developing genetic markers for population genetic studies in amphibians is critical because amphibians are declining globally. The tiger salamander, Ambystoma tigrinum, is widespread throughout the United States and includes the endangered subspecies, A. t. stebbinsi. We present primers and amplification conditions for 10 polymorphic microsatellite loci that have produced successful results in three subspecies of A. tigrinum. Number of alleles per locus ranged from one to 11 and heterozygosity ranged from 0 to 0.815 depending on the subspecies and locus analysed. These markers should prove useful for future studies of genetic diversity and population subdivision.     

Regional genetic structuring and evolutionary history of the impala Aepyceros melampus

Samples of 162 impala antelope (Aepyceros melampus) from throughout its distribution range in sub-Saharan Africa were surveyed using eight polymorphic microsatellite loci. Furthermore, 155 previously published mitochondrial DNA (mtDNA) sequences from the same localities were reanalyzed. Two subspecies of impala are presently recognized--the isolated black-faced impala (Aepyceros melampus petersi) in southwest Africa and the common impala (Aepyceros melampus melampus) abundant in southern and east Africa. All tests performed indicated significant genetic differentiation at the subspecific level. Furthermore, individual-based analyses split the common impala subspecies into two distinct genetic groups, conforming with regional geographic affiliation to southern or east Africa. This was supported by assignment tests, genetic distance measures, pairwise theta values, and analysis of molecular variance. We suggest that the presence of such previously unknown regional structuring within the subspecies reflects a pattern of colonization from a formerly large panmictic population in southern Africa toward east Africa. This scenario was supported by a progressive decline in population diversity indices toward east Africa and a significant increase in the quantity theta/(1 - theta). Both microsatellite and mtDNA data indicated a genetic distinctiveness of the Samburu population in Kenya.     

Positive relationships between genetic diversity and abundance in fishes

Molecular markers, such as mitochondrial DNA and microsatellite loci, are widely studied to assess population genetics and phylogeography; however, the selective neutrality of these markers is increasingly being questioned. Given the importance of molecular markers in fisheries science and conservation, we evaluated the neutrality of both mtDNA and microsatellite loci through their associations with population size. We surveyed mtDNA and microsatellite data from the primary literature and determined whether genetic diversity increased with abundance across a total of 105 marine and freshwater fishes, with both global fisheries catch data and body size as proxies for abundance (with an additional 57 species for which only body size data were assessed). We found that microsatellite data generally yielded higher associations with abundance than mtDNA data, and within mtDNA analyses, number of haplotypes and haplotype diversity were more strongly associated with abundance than nucleotide diversity, particularly for freshwater fishes. We compared genetic diversity between freshwater and marine fishes and found that marine fishes had higher values of all measures of genetic diversity than freshwater fishes. Results for both mtDNA and microsatellites generally conformed to neutral expectations, although weaker relationships were often found between mtDNA nucleotide diversity and ‘abundance’ compared to any other genetic statistic. We speculate that this is because of historical events unrelated to natural selection, although a role for selection cannot be ruled out.     

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.     

Population differentiation in the Pacific white-sided dolphin Lagenorhynchus obliquidens inferred from mitochondrial DNA and microsatellite analyses

We investigated genetic diversity and differentiation of the Pacific white-sided dolphin (Lagenorhynchus obliquidens) in Japanese coastal waters and offshore North Pacific by analyzing mitochondrial DNA and nuclear microsatellite variation. A total of 519 bp of the mitochondrial control region was sequenced and five microsatellite locus were genotyped for 59 individuals. A high level of haplotypic diversity (h=96.1%), moderate level of nucleotide diversity (pi=1.65%) and average expected heterozygosity (HE=0.66-0.76) were within an extent of those reported for other odontocetes. Consistent genetic difference between the samples from Japanese coastal Pacific-Sea of Japan and offshore North Pacific was indicated by analyses of molecular variance (AMOVAs) based on mtDNA and microsatellite variations, comparison of genetic variabilities, and geographical distributions of mtDNA haplotypes and microsatellite alleles. This result suggests that Pacific white-sided dolphins in each of the above two areas belong to different populations between which gene flow has been severely restricted. The low genetic diversity and mtDNA genealogy of the population in Japanese coastal waters suggest that it originated from a small population that colonized the Sea of Japan or that experienced population reduction when this Sea was isolated from the North Pacific during a glacial period in the Late Pleistocene.     

滇金丝猴线粒体D-loop 片段SNPs 位点的鉴别及初步分析

本文利用已测序的157 个滇金丝猴控制区(D-loop)片段,通过与参考序列比对,鉴别了线粒体D-loop 片段中的52 个SNP (Single Nucleotide Polymorphisms)位点,定义了30 种滇金丝猴单倍型,排除概率为0. 938。谱系及种群遗传结构分析结果与以前利用D-loop 片段的研究结果相似。同时表明基于粪便样品进行保护遗传学、谱系生物地理学、种群遗传学等研究时,与线粒体标记和微卫星标记相比,SNP 标记可能具有一定的优越性,并建议进一步分析滇金丝猴线粒体D-loop 全序列甚至线粒体全基因组上的SNPs 位点的信息,以促进滇金丝猴保护遗传学等研究的开展。     

Intraspecific genetic variation in the common midwife toad (Alytes obstetricans): subspecies assignment using mitochondrial and microsatellite markers

The common midwife toad (Alytes obstetricans), widely distributed in the northern half of the Iberian Peninsula and part of Western Europe, is currently subdivided into four subspecies: A. o. obstetricans, A. o. boscai, A. o. pertinax and A. o. almogavarii. However, the delimitation of these subspecies and their ranges are still under discussion because strong discordances have been found between morphological and molecular data, and especially among different genetic markers. Here, we screen a set of novel microsatellite loci and mtDNA sequences of A. obstetricans populations representative of all currently recognized subspecies to investigate the correspondence between genetic groupings inferred from clustering analysis of microsatellite genotypes and the described subspecies and test whether patterns of mtDNA variation are concordant with those genetic clusters. Our results confirm previous expectations of extremely high intraspecific diversity in A. obstetricans in Iberia. Analyses of microsatellite and mtDNA data were concordant in recovering five well‐defined groups, of which three correspond to previously defined subspecies, while the two additional clusters correspond to populations of subspecies A. o. boscai separated by the Douro River. Our results suggest the occurrence of two distinct genetic units within A. o. boscai that likely result from a long independent evolutionary history, thus deserving special attention from a conservation point of view.     

Low support for separate species within the redpoll complex (Carduelis flammea-hornemanni-cabaret) from analyses of mtDNA and microsatellite markers

The redpoll complex, consisting of three currently recognized species (Carduelis flammea, C. hornemanni and C. cabaret), is polytypic in biometry, morphology, physiology and behaviour. However, previous genetic work has not revealed any indications of genetic differentiation. We analysed sequence variation in the mtDNA control region, and allele frequencies of supposedly faster evolving microsatellites (n=10), in an attempt to detect molecular genetic support for the three species, as well as two subspecies of C. flammea (ssp. flammea and rostrata), within this complex. We used samples from two subspecies of the twite (Carduelis flavirostris, ssp. flavirostris and rufostrigata) as outgroup. We found no structure among redpoll individuals in mtDNA haplotypes or microsatellite allele frequencies, and only marginal differences between redpoll taxa in analyses of molecular variance (AMOVAs) of predefined groups. In contrast, the two twite subspecies constituted two well-supported monophyletic groups. Our study thus strengthens previous indications of low genetic support for current redpoll taxa. Two major alternative interpretations exist. Either redpolls form a single gene pool with geographical polymorphisms possibly explained by Bergmann's and Gloger's rules, or there are separate gene pools of recent origin but with too little time elapsed for genetic differentiation to have evolved in the investigated markers. Future studies should therefore examine whether reproductive isolation mechanisms and barriers to gene flow exist in areas with sympatric breeding.     

Ancestral polymorphisms in genetic markers obscure detection of evolutionarily distinct populations in the endangered Florida grasshopper sparrow (Ammodramus savannarum floridanus)

Genetic analyses of bird subspecies designated as conservation units can address whether they represent units with independent evolutionary histories and provide insights into the evolutionary processes that determine the degree to which they are genetically distinct. Here we use mitochondrial DNA control region sequence and six microsatellite DNA loci to examine phylogeographical structure and genetic differentiation among five North American grasshopper sparrow (Ammodramus savannarum) populations representing three subspecies, including a population of the endangered Florida subspecies (A. s. floridanus). This federally listed taxon is of particular interest because it differs phenotypically from other subspecies in plumage and behaviour and has also undergone a drastic decline in population size over the past century. Despite this designation, we observed no phylogeographical structure among populations in either marker: mtDNA haplotypes and microsatellite genotypes from floridanus samples did not form clades that were phylogenetically distinct from variants found in other subspecies. However, there was low but significant differentiation between Florida and all other populations combined in both mtDNA (FST = 0.069) and in one measure of microsatellite differentiation (theta = 0.016), while the non-Florida populations were not different from each other. Based on analyses of mtDNA variation using a coalescent-based model, the effective sizes of these populations are large (approximately 80,000 females) and they have only recently diverged from each other (< 26,000 ybp). These populations are probably far from genetic equilibrium and therefore the lack of phylogenetic distinctiveness of the floridanus subspecies and minimal genetic differentiation is due most probably to retained ancestral polymorphism. Finally, levels of variation in Florida were similar to other populations supporting the idea that the drastic reduction in population size which has occurred within the last 100 years has not yet had an impact on levels of variation in floridanus. We argue that despite the lack of phylogenetic distinctiveness of floridanus genotypes the observed genetic differentiation and previously documented phenotypic differences justify continued designation of this subspecies as a protected population segment.     

The role of the Ord Arid Intrusion in the historical and contemporary genetic division of long-tailed finch subspecies in northern Australia

The effect of separation by biogeographic features followed by secondary contact can blur taxonomic boundaries and produce complex genetic signatures. We analyzed population structure and gene flow across the range of the long-tailed finch (Poephila acuticauda) in northern Australia (1) to test the hypothesis that Ord Arid Intrusion acted as the causative barrier that led to divergence of P. acuticauda subspecies, (2) to determine whether genetic data support the presence of a gradual cline across the range or a sudden shift, both of which have been suggested based on morphological data, and (3) to estimate levels of contemporary gene flow within this species complex. We collected samples from 302 individuals from 10 localities. Analyses of 12 microsatellite loci and sequence data from 333 base pairs of the mitochondrial control region were used to estimate population structure and gene flow, using analysis of molecular variance (AMOVA), haplotype network analysis, frequency statistics, and clustering methods. Mitochondrial sequence data indicated the presence of three genetic groups (regions) across the range of P. acuticauda. Genetic diversity was highest in the east and lowest in the west. The Ord Arid Intrusion appears to have functioned as a biogeographic barrier in the past, according to mtDNA evidence presented here and evidence from previous studies. The absence of isolation by distance between adjacent regions and the lack of population genetic structure of mtDNA within regions indicates that genetic changes across the range of P. acuticauda subspecies are characterized by discrete breaks between regions. While microsatellite data indicate a complete absence of genetic structure across this species' range, it appears unlikely that this results from high levels of gene flow. Mitochondrial data do not support the presence of contemporary gene flow across the range of this species.     

Reduced genetic variation and strong genetic population structure in the freshwater killifish Valencia letourneuxi (Valenciidae) based on nuclear and mitochondrial markers

The genetic variation of the critically endangered Corfu killifish (Valencia letourneuxi), an endemic freshwater fish species of the western Balkans, was assessed for nine populations sampled in eight water systems in western continental Greece, the Peloponnese and the Ionian Island of Corfu, using mitochondrial and microsatellite markers. The analyses were based on data from three mtDNA regions (D‐loop, COI and 16S rRNA sequences) and 14 microsatellite loci. Samples from the congeneric species Valencia hispanica and the phylogenetically closely related species Aphanius fasciatus were also used in the study as outgroups. Both the mitochondrial and the microsatellite analyses revealed three distinct population groupings associated with the geographical distribution of the populations: one southern group occupying rivers draining to the Patraikos Gulf, the second one including the populations flowing into the Amvrakikos Gulf and the third, more northern group, including the other populations from rivers in Corfu Island and Epirus flowing into the Ionian Sea. Within these groupings there is limited genetic differentiation between populations; in addition, there is reduced intrapopulation genetic variation, evidenced by low heterozygosity values, number of alleles and haplotype diversity. In terms of taxonomic implications and appropriate management actions for conservation, our data suggest that the major population groups should be regarded at least as three distinct conservation units (CUs), with translocation and restocking actions to take place only within the geographical range of the CU concerned. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 334–349.     

Genetic population structure in the yellow mongoose, Cynictis penicillata

Phylogeographic structure was determined for the yellow mongoose, Cynictis penicillata , using mtDNA RFLPs and control region sequences. The RFLP analysis revealed 13 haplotypes which showed weak geographical patterning consistent with a recent range expansion from a refugial population(s). An analysis of molecular variance ( AMOVA ) revealed no correspondence between mtDNA phylogeography and subspecies delimitation, nor between matrilines and areas characterized by a high incidence of the viverrid-type rabies, of which the yellow mongoose is the principal vector. The lack of structure was also shown by control region sequences although four of the maternal lineages shared a near-perfect 81 bp repeat. We speculate that regional hot spots of the viverrid rabies biotype reflect population density differences in the yellow mongoose that are not underscored by genetic partitioning, at least at the level of resolution provided by our analyses.     

Genetic structuring and recent demographic history of red pandas (Ailurus fulgens) inferred from microsatellite and mitochondrial DNA

Clarification of the genetic structure and population history of a species can shed light on the impacts of landscapes, historical climate change and contemporary human activities and thus enables evidence‐based conservation decisions for endangered organisms. The red panda (Ailurus fulgens) is an endangered species distributing at the edge of the Qinghai‐Tibetan Plateau and is currently subject to habitat loss, fragmentation and population decline, thus representing a good model to test the influences of the above‐mentioned factors on a plateau edge species. We combined nine microsatellite loci and 551 bp of mitochondrial control region (mtDNA CR) to explore the genetic structure and demographic history of this species. A total of 123 individuals were sampled from 23 locations across five populations. High levels of genetic variation were identified for both mtDNA and microsatellites. Phylogeographic analyses indicated little geographic structure, suggesting historically wide gene flow. However, microsatellite‐based Bayesian clustering clearly identified three groups (Qionglai‐Liangshan, Xiaoxiangling and Gaoligong‐Tibet). A significant isolation‐by‐distance pattern was detected only after removing Xiaoxiangling. For mtDNA data, there was no statistical support for a historical population expansion or contraction for the whole sample or any population except Xiaoxiangling where a signal of contraction was detected. However, Bayesian simulations of population history using microsatellite data did pinpoint population declines for Qionglai, Xiaoxiangling and Gaoligong, demonstrating significant influences of human activity on demography. The unique history of the Xiaoxiangling population plays a critical role in shaping the genetic structure of this species, and large‐scale habitat loss and fragmentation is hampering gene flow among populations. The implications of our findings for the biogeography of the Qinghai‐Tibetan Plateau, subspecies classification and conservation of red pandas are discussed.     

Phylogeography and genetic ancestry of tigers (Panthera tigris)

Eight traditional subspecies of tiger (Panthera tigris),of which three recently became extinct, are commonly recognized on the basis of geographic isolation and morphological characteristics. To investigate the species' evolutionary history and to establish objective methods for subspecies recognition, voucher specimens of blood, skin, hair, and/or skin biopsies from 134 tigers with verified geographic origins or heritage across the whole distribution range were examined for three molecular markers: (1) 4.0 kb of mitochondrial DNA (mtDNA) sequence; (2) allele variation in the nuclear major histocompatibility complex class II DRB gene; and (3) composite nuclear microsatellite genotypes based on 30 loci. Relatively low genetic variation with mtDNA,DRB,and microsatellite loci was found, but significant population subdivision was nonetheless apparent among five living subspecies. In addition, a distinct partition of the Indochinese subspecies P. t. corbetti in to northern Indochinese and Malayan Peninsula populations was discovered. Population genetic structure would suggest recognition of six taxonomic units or subspecies: (1) Amur tiger P. t. altaica; (2) northern Indochinese tiger P. t. corbetti; (3) South China tiger P. t. amoyensis; (4) Malayan tiger P. t. jacksoni, named for the tiger conservationist Peter Jackson; (5) Sumatran tiger P. t. sumatrae; and (6) Bengal tiger P. t. tigris. The proposed South China tiger lineage is tentative due to limited sampling. The age of the most recent common ancestor for tiger mtDNA was estimated to be 72,000-108,000 y, relatively younger than some other Panthera species. A combination of population expansions, reduced gene flow, and genetic drift following the last genetic diminution, and the recent anthropogenic range contraction, have led to the distinct genetic partitions. These results provide an explicit basis for subspecies recognition and will lead to the improved management and conservation of these recently isolated but distinct geographic populations of tigers.