Species boundaries in Malagasy snakes of the genus Madagascarophis (Serpentes: Colubridae sensu lato) assessed by nuclear and mitochondrial markers

We studied mitochondrial divergence in 27 individuals of colubrid snakes of the genus Madagascarophis Mertens from most of its distribution area in Madagascar. Combined analyses of 16S rRNA and cytochrome b sequences identified six major clades which only partly agreed with previously proposed classifications. Analysis of nuclear DNA sequences of the c-mos gene as well as of ISSR fingerprints revealed consistent differences only among three clades which we consider as distinct species: a widespread Madagascarophis colubrinus (Schlegel), with M. citrinus (Boettger) as a junior synonym, a southern M. meridionalis Domergue, and a presumably undescribed species from the extreme north of Madagascar. The species M. ocellatus Domergue was not available for our study. Within M. colubrinus there are two populations from the north-west, each showing two divergent haplotypes with pairwise divergences of up to 5.2% in the cytochrome b gene. Maximum divergence in this gene within M. colubrinus was 7.1%. These high values emphasise that caution needs to be applied before genetic distance values are used for species delimitation. Phylogeographically, most of the genetic variation in M. colubrinus is found in northern Madagascar, indicating that the species might have originated in this region. Later one haplotype clade colonised western and eastern Madagascar, with a putative secondary introgression into north-western populations.     

Phylogeography of the introduced species Rattus rattus in the western Indian Ocean,with special emphasis on the colonization history of Madagascar

Aim To describe the phylogeographic patterns of the black rat, Rattus rattus, from islands in the western Indian Ocean where the species has been introduced (Madagascar and the neighbouring islands of Réunion, Mayotte and Grande Comore), in comparison with the postulated source area (India). Location Western Indian Ocean: India, Arabian Peninsula, East Africa and the islands of Madagascar, Réunion, Grande Comore and Mayotte. Methods Mitochondrial DNA (cytochrome b, tRNA and D‐loop, 1762 bp) was sequenced for 71 individuals from 11 countries in the western Indian Ocean. A partial D‐loop (419 bp) was also sequenced for eight populations from Madagascar (97 individuals), which were analysed in addition to six previously published populations from southern Madagascar. Results Haplotypes from India and the Arabian Peninsula occupied a basal position in the phylogenetic tree, whereas those from islands were distributed in different monophyletic clusters: Madagascar grouped with Mayotte, while Réunion and Grand Comore were present in two other separate groups. The only exception was one individual from Madagascar (out of 190) carrying a haplotype that clustered with those from Réunion and South Africa. ‘Isolation with migration’ simulations favoured a model with no recurrent migration between Oman and Madagascar. Mismatch distribution analyses dated the expansion of Malagasy populations on a time‐scale compatible with human colonization history. Higher haplotype diversity and older expansion times were found on the east coast of Madagascar compared with the central highlands. Main conclusions Phylogeographic patterns supported the hypothesis of human‐mediated colonization of R. rattus from source populations in either the native area (India) or anciently colonized regions (the Arabian Peninsula) to islands of the western Indian Ocean. Despite their proximity, each island has a distinct colonization history. Independent colonization events may have occurred simultaneously in Madagascar and Grande Comore, whereas Mayotte would have been colonized from Madagascar. Réunion was colonized independently, presumably from Europe. Malagasy populations may have originated from a single successful colonization event, followed by rapid expansion, first in coastal zones and then in the central highlands. The congruence of the observed phylogeographic pattern with human colonization events and pathways supports the potential relevance of the black rat in tracing human history.     

Invasion genetics of a human commensal rodent: the black rat Rattus rattus in Madagascar

Studies focusing on geographical genetic patterns of commensal species and on human history complement each other and provide proxies to trace common colonization events. On Madagascar, the unintentional introduction and spread of the commensal species Rattus rattus by people may have left a living clue of human colonization patterns and history. In this study, we addressed this question by characterizing the genetic structure of natural populations of R. rattus using both microsatellites and mitochondrial sequences, on an extensive sampling across the island. Such data sets were analysed by a combination of methods using population genetics, phylogeography and approximate Bayesian computation. Our results indicated two introduction events to Madagascar from the same ancestral source of R. rattus, one in the extreme north of the island and the other further south. The latter was the source of a large spatial expansion, which may have initially started from an original point located on the southern coast. The inferred timing of introduction events—several centuries ago—is temporally congruent with the Arabian trade network in the Indian Ocean, which was flourishing from the middle of the first millennium.     

Phylogenetic status of brown trout Salmo trutta populations in five rivers from the southern Caspian Sea and two inland lake basins,Iran: a morphogenetic approach

Interrelationships, origin and phylogenetic affinities of brown trout Salmo trutta populations from the southern Caspian Sea basin, Orumieh and Namak Lake basins in Iran were analysed from complete mtDNA control region sequences, 12 microsatellite loci and morphological characters. Among 129 specimens from six populations, seven haplotypes were observed. Based on mtDNA haplotype data, the Orumieh and southern Caspian populations did not differ significantly, but the Namak basin–Karaj population presented a unique haplotype closely related to the haplotypes of the other populations (0·1% Kimura two‐parameter, K2P divergence). All Iranian haplotypes clustered as a distinct group within the Danube phylogenetic grouping, with an average K2P distance of 0·41% relative to other Danubian haplotypes. The Karaj haplotype in the Namak basin was related to a haplotype (Da26) formerly identified in the Tigris basin in Turkey, to a Salmo trutta oxianus haplotype from the Aral Sea basin, and to haplotype Da1a with two mutational steps, as well as to other Iranian haplotypes with one to two mutational steps, which may indicate a centre of origin in the Caspian basin. In contrast to results of the mtDNA analysis, more pronounced differentiation was observed among the populations studied in the morphological and microsatellite DNA data, except for the two populations from the Orumieh basin, which were similar, possibly due to anthropogenic causes.     

Genetic characterization of common carp (Cyprinus carpio) populations from Greece using mitochondrial DNA sequences

Wild common carp from two lakes and two rivers in Greece were genetically characterized with sequencing analysis of two mitochondrial DNA segments: cytochrome b (1119 bp) and D-loop (646 bp). A total of 9 variable singleton sites and 7 unique haplotypes were detected. A common haplotype was found in three out of the four populations examined, which seems to be the ancestral one and represents the European origin of common carp from Greece. This haplotype could be also justified by the introductions reported with individuals belonging to the Central European race, into many natural habitats in Greece. Limited genetic variation — in Evros and Aliakmonas populations — could be due to bottleneck effects and small effective population sizes, whereas the different haplotypes found in Lake Volvi could represent different common carp stocks. Values of sequence divergence among Greek haplotypes ranged from 0.0006 to 0.0023. The Neighbour-Joining (NJ) phylogenetic tree constructed based on the combined sequences, reveals that the populations of common carp from Greece belong to the European group of populations — which is highly divergent from the South East-Asia cluster — and to the subspecies Cyprinus carpio carpio.     

基于叶绿体DNA非编码序列的蒙古扁桃谱系地理学研究

该研究选取中国西北干旱区第三纪孑遗植物蒙古扁桃（Amygdalus mongolica）,基于叶绿体DNA非编码trnH psbA序列对蒙古扁桃17个居群324个个体进行了谱系地理学研究。结果表明：（1）蒙古扁桃trnH psbA序列长度350 bp,变异位点63个,共有9种单倍型,居群间总遗传多样性为（H_t）为0.758,居群内平均遗传多样性为（H_s）为0.203,贺兰山东麓及阴山南麓边缘的居群具有较高的单倍型多样性及核苷酸多样性并固定较多特有单倍型,推测这2个地区是蒙古扁桃在第四纪冰期时的重要避难所。（2）AMOVA分析表明,居群间的遗传变异为83.84%,居群内的遗传变异为16.16%,居群间遗传分化系数N_st>G_st（N_st=0.733, G_st=0.655, P>0.05）,表明蒙古扁桃不存在明显的谱系地理结构;根据单倍型地理分布及网络关系图,把蒙古扁桃自然地理居群分为东、西两大地理组群,而且东、西地理组群没有共享单倍型;居群遗传结构分析表明,两大地理组群遗传分化较大。（3）蒙古扁桃居群在间冰期或冰期后经历了近期的居群扩张,由于奠基者效应使得多数居群只固定了单一的单倍型。     

Phylogeny and genetic structure of the goitered gazelle (Artiodactyla,Bovidae) in north-western China indicated by the hypervariable mitochondrial control region

The goitered gazelle, Gazella subgutturosa, is a medium-sized ungulate inhabiting arid and semi-arid regions in the Middle East and central Asia. The intraspecific classification of the species remains unclear. We analysed the genetic diversity in mitochondrial DNA control region (CR) sequences (976?bp) from 104 wild samples from the Xinjiang Uyghur Autonomous Region (XUAR) in north-west China, and reconstructed phylogeny with additional sequences from across the species’ range. We detected 58 haplotypes in XUAR populations, all but three of which were specific to single sampling sites. The phylogenetic analysis displayed two obvious clades of mtDNA haplotypes and the other haplotypes differed from the two clades. A median-joining network showed three groups of haplotypes were to a high extent concordant with the phylogenetic tree. The haplotype clustering was consistent with their geographic distribution. Nei’s net sequence divergences amongst the three groups ranged from 0.010 to 0.018 and indicated three subspecies, two of which inhabit XUAR. We detected strong differentiation between northern (NX) and southern (SX) XUAR populations overall (F_ST?=?0.4448, P?相似文献   

No genetic differentiation between geographically isolated populations of Clarias macrocephalus Günther in Malaysia revealed by sequences of mtDNA Cytochrome b and D‐loop gene regions

In the present study, we assessed the genetic variation of three Clarias macrocephalus Günther populations collected from Kedah, Perlis and Kelantan (Peninsular Malaysia) using sequences of partial mitochondrial cytochrome b (Cyt b) and D‐loop genes. A total of 57 individuals were sequenced and 1470 bp were obtained (1053 bp Cyt‐b; 417 bp D‐loop). The analysis revealed 21 haplotypes based on 81 polymorphic sites. Nucleotide diversity (π) was 0.003 in all populations while haplotype diversity ranged from 0.657 to 0.765. No significant genetic differentiation among the three populations was observed. Nevertheless, a number of private haplotypes was discovered, providing valuable information for selective breeding programs.     

Mitochondrial haplotypes in populations of the plant-infecting fungus Sclerotinia sclerotiorum: wide distribution in agriculture, local distribution in the wild

Analysis of mitochondrial DNA (mtDNA) haplotypes of Sclerotinia sclerotiorum points to a common origin of some genotypes from agricultural populations, especially when compared with two wild populations that are sharply distinguished from the agricultural sample and from each other. Five agricultural population samples from canola (Alberta, Canada and Norway), cabbage (North Carolina, USA), sunflower (Manitoba, Canada and Queensland, Australia) and two Norwegian populations from a wild plant, Ranunculus ficaria were compared. Haplotypes were determined by Southern hybridization of purified organelle DNA from S. sclerotiorum and Neurospora crassa to total genomic DNA of S. sclerotiorum. Each isolate had one haplotype. Haplotypes of S. sclerotiorum from R. ficaria were different between the two wild populations and also from all haplotypes observed in the agricultural populations. Among the wild isolates, DNA fingerprint, mtDNA haplotype and location in the sampling transect were all associated. Among the agricultural isolates, four haplotypes were observed in at least two agricultural populations and one haplotype was observed in all agricultural populations. In the Canadian canola sample some clones had one mtDNA haplotype, indicating association with DNA fingerprint, some clones had more than one haplotype, and some groups of clones shared haplotypes. Some of the haplotype diversity may be due to the presence of extra-chromosomal elements associated with the mitochondria of S. sclerotiorum.     

Revealing the geographic origin of an invasive lizard: the problem of native population genetic diversity

The brown anole, Anolis sagrei, is one of the most widespread and successful colonisers of the diverse Anolis genus, which comprises c. 400 species occurring naturally in Central and South America and the Caribbean. Based on extensive between and within population sampling from a previously published study (334 mitochondrial DNA sequences) and sampling for this study (37 mtDNA sequences), we reconstruct a phylogeny and produce a haplotype network to assign a recently introduced population in St Vincent, Lesser Antilles to its geographic origin. A single haplotype was present in the St Vincent population, which was identical to a haplotype from Tampa, FL. We show that genetic diversity within native range populations, combined with low frequencies of introduced haplotypes in native ranges, may impair attempts to identify source populations, even despite intensive sampling effort. The absence of mtDNA haplotype diversity suggests a significant genetic founder effect within the St Vincent population.     

Genetic structure of <Emphasis Type="Italic">Varroa destructor</Emphasis> populations infesting <Emphasis Type="Italic">Apis mellifera</Emphasis> colonies in Argentina

Although mitochondrial DNA mapping of Varroa destructor revealed the presence of several haplotypes, only two of them (Korean and Japanese haplotypes) were capable to infest Apis mellifera populations. Even though the Korean haplotype is the only one that has been reported in Argentina, these conclusions were based on mites sampled in apiaries from a specific geographical place (Buenos Aires province). To study mites from several sites of Argentina could reveal the presence of the Japanese genotype, especially considering sites near to Brazil, where Japanese haplotype was already detected. The aim of this work was to study the genetic structure of V. destructor populations from apiaries located in various provinces of Argentina, in order to determine the presence of different haplotypes. The study was carried out between January 2006 and December 2009. Phoretic adult Varroa mites were collected from honey bee workers sampled from colonies of A. mellifera located in Entre Ríos, Buenos Aires, Corrientes, Río Negro, Santa Cruz and Neuquén provinces. Twenty female mites from each sampling site were used to carry out the genetic analysis. For DNA extraction a nondestructive method was used. DNA sequences were compared to Korean haplotype (AF106899) and Japanese haplotype (AF106897). All DNA sequences obtained from mite populations sampled in Argentina, share 98% of similitude with Korean Haplotype (AF106899). Taking into account these results, we are able to conclude that Korean haplotype is cosmopolite in Argentina.     

Population genetic structure of the seed pathogen Pyrenophora semeniperda on Bromus tectorum in western North America

We examined genetic variation in the ascomycete pathogen Pyrenophora semeniperda cultured from seeds of the invasive grass Bromus tectorum in the Intermountain West of North America. We sequenced the internal transcribed spacer (ITS) region of the nuclear ribosomal RNA genome in 417 monoconidial cultures collected from 20 sites in Washington, Idaho, Utah and Colorado, USA. ITS sequence diversity was surprisingly high; 12 unique haplotypes were identified, averaging 1.3% pairwise sequence divergence. All sites had at least two haplotypes present, and three sites had seven or more. One haplotype composed 60% of the isolates and occurred at all 20 locations; the remaining haplotypes generally occurred at low frequencies within sites but at multiple sites throughout the region. Sites in Washington and Idaho were more diverse than those in Utah and Colorado, averaging two more haplotypes and 67% more pairwise differences among haplotypes at a site. Analysis of molecular variance (AMOVA) indicated that more than 80% of the genetic variation was found within sampling locations, while 7-11% of the variation can be attributed to differences between northern (Washington and Idaho) and southern (Utah and Colorado) populations. The wide distribution of even uncommon haplotypes among sampling sites and weak correlations between genetic and geographic distances among populations (< 0.2) suggested that these populations recently were established from a common source. We hypothesize that the strains of P. semeniperda infecting B. tectorum in western North America probably arrived with the invasive grass from its native Eurasian range.     

The origin of Rattus rattus on the Îles Éparses,Western Indian Ocean

Tollenaere et al. (Journal of Biogeography, 2010, 37 , 398–410) present a phylogeographic analysis of Rattus rattus for the Western Indian Ocean, with particular emphasis on Madagascar, but do not include samples from three island groups centrally located in the Mozambique Channel. Haplotypes from these islands provide additional information on the colonization pathways of R. rattus in the Western Indian Ocean region. For each of the three Îles Éparses groups in the Mozambique Channel, we test the competing hypotheses that colonization by R. rattus was most likely: (1) from the Arabian Peninsula, (2) from East Africa, (3) from Madagascar, or (4) from independent shipping. These results are combined with historical observations of the presence of R. rattus on these islands to give stronger inference on the colonization pathways. Additionally, more accurate colonization dates provide guidance for contemporary conservation management.     

Invasion by Rattus rattus into native coastal forests of south‐eastern Australia: are native small mammals at risk?

The black rat, Rattus rattus, is an alien rodent in Australian ecosystems where niche overlap with native small mammals may lead to competition for resources and displacement of native species. In coastal habitats surrounding Jervis Bay in south‐eastern Australia, R. rattus co‐occurs with the native bush rat, Rattus fuscipes, and brown antechinus, Antechinus stuartii. Relative distributions and abundances, and fine‐scale space use suggest invasive and native rodents compete for use of space and habitat. Such competitive interactions were not evident between R. rattus and native A. stuartii, which was negatively influenced more by disturbance to habitat. Differences in rodent communities between spatially separate forests forming the northern and southern peninsulas of Jervis Bay potentially reflect symmetrical competition and differences in competitive outcomes. In southern forests, R. rattus was largely restricted to patches of disturbed forest associated with campgrounds. Competitive interference by native rodent populations inhabiting surrounding intact forests may have so far limited R. rattus colonization of these areas. In northern forests, R. rattus was the predominant rodent irrespective of disturbance, while populations of R. fuscipes were unusually low seemingly due to poor juvenile recruitment. Native individuals avoided areas frequented by adult R. rattus and given that species did not partition use of microhabitats, R. rattus most likely precluded R. fuscipes from suitable habitat and in doing so limited native populations. We discuss how natural disturbance of habitat and human activity have potentially facilitated successful invasion by R. rattus of the northern forests. Studies that manipulate rodent populations are required to support these interpretations of observed patterns.     

Molecular diversity and phylogeography of Andricus curtisii (Hymenoptera,Cynipidae) in Turkey

A 433-bp fragment of the mitochondrial cytochrome b gene and the entire nuclear ITS2 region were sequenced in an oak gall wasp species, Andricus curtisii, which represents 21 populations across its range in Turkey. Forty cytochrome b haplotypes and 63 ITS2 alleles were discovered among the 161 individuals studied. Remarkable genetic variation was detected in the species, with an average haplotype diversity (h) of 0.72 and 0.76 and an average nucleotide diversity (π) of 0.06 and 0.01 for the cyt b gene and the ITS region, respectively. Most populations were genetically distinct, possessing unique cytochrome b haplotypes. Maximum likelihood and Bayesian phylogenetic analyses and application of the haplotype network to both datasets revealed that the most basally located haplotypes/alleles were from the eastern part of Turkey. The detection of older sequences in the phylogenetic trees of the easterly located populations implied an eastern origin of the species. The application of a conventional insect mitochondrial DNA clock to the dataset suggested a splitting of ingroup haplotypes from the outgroup lineages predating the Pleistocene epoch and the formation of two major haplogroups that coincided with the oscillation of the early Pleistocene glaciations.     

Genetic Diversity of Moose (Alces alces L.) in Eurasia

Polymorphism of nucleotide sequence of D-loop fragment of the mitochondrial DNA was studied in 20 moose from several local populations on the territory of Eurasia. Three main haplotype variants of D-loop were detected by molecular phylogenetic method, which formed three clusters named European, Asian, and American. Intraspecies variation in the length of HVSI of D-loop of the mitochondrial DNA of moose was revealed. In the Far Eastern and Yakutian moose, haplotypes with a 75-bp deletion were found, which were most similar with haplotypes (also with the deletion), earlier observed in North American moose [1]. The highest diversity of the haplotypes of mitochondrial DNA is characteristic of Yakutia and the Far East (where three haplotype variants were found), which demonstrates the probable role of the region as the center of the species origin or as the region of ancient population mixture. The geographic region might be considered as a probable source of ancient moose migrations from Asia to America, basing on the data of distribution of mitochondrial haplotypes of D-loop and alleles of MhcAlal-DRB1. Divergence of nucleotide sequences of haplotypes with the 75-bp deletion (forming the American cluster on the phylogenetic tree) was the lowest (0.4%), which evidences respectively recent origin of the group of haplotypes. In Europe, only haplotypes of mitochondrial DNA referred to European variant were observed. Basing on analysis of variation of nucleotide sequences of D-loop, exon 4 of -Casein and exon 2 of MhcAlal-DRB1, we demonstrated that Eurasian moose studied belong to the unique species, which has probably passed through a bottle neck. The time of the origin of modern diversity of D-loop haplotypes of the species was estimated as 0.075–0.15 Myr ago.     

Disrupted phylogeographical microsatellite and chloroplast DNA patterns indicate a vicariance rather than long‐distance dispersal origin for the disjunct distribution of the Chilean endemic Dioscorea biloba (Dioscoreaceae) around the Atacama Desert

Aim The Chilean endemic Dioscorea biloba (Dioscoreaceae) is a dioecious geophyte that shows a remarkable 600 km north–south disjunction in the peripheral arid area of the Atacama Desert. Its restricted present‐day distribution and probable Neogene origin indicate that its populations have a history linked to that of the Atacama Desert, making this an ideal model species with which to investigate the biogeography of the region. Location Chile, Atacama Desert and peripheral arid area. Methods Two hundred and seventy‐five individuals from nine populations were genotyped for seven nuclear microsatellite loci, and plastid trnL–F and trnT–L sequences were obtained for a representative subset of these. Analyses included the estimation of genetic diversity and population structure through clustering, Bayesian and analysis of molecular variance analyses, and statistical parsimony networks of chloroplast haplotypes. Isolation by distance was tested against alternative dispersal hypotheses. Results Microsatellite markers revealed moderate to high levels of genetic diversity within populations, with those from the southern Limarí Valley showing the highest values and northern populations showing less exclusive alleles. Bayesian analysis of microsatellite data identified three genetic groups that corresponded to geographical ranges. Chloroplast phylogeography revealed no haplotypes shared between northern and southern ranges, and little haplotype sharing between the two neighbouring southern valleys. Dispersal models suggested the presence of extinct hypothetical populations between the southern and northern ranges. Main conclusions Our results are consistent with prolonged isolation of the northern and southern groups, mediated by the life‐history traits of the species. Significant isolation was revealed at both large and moderate distances as gene flow was not evident even between neighbouring valleys. Bayesian analyses of microsatellite and chloroplast haplotype diversity identified the southern area of Limarí as the probable area of origin of the species. Our data do not support recent dispersal of D. biloba from the southern range into Antofagasta, but indicate the fragmentation of an earlier wider range, concomitant with the Pliocene–Pleistocene climatic oscillations, with subsequent extinctions of the Atacama Desert populations and the divergence of the peripheral ones as a consequence of genetic drift.     

中国湛江市和哈尔滨市褐家鼠种群遗传结构及其年度变化特征

为探索褐家鼠Rattus norvegicus地理种群的遗传结构及其年度变化特点,本研究以广东省湛江市的褐家鼠指名亚种和黑龙江省哈尔滨市的褐家鼠东北亚种为主要研究对象,结合我国及世界其他褐家鼠种群的D-loop序列分析这2个褐家鼠地理种群间D-loop序列的遗传分化情况及系统进化关系,重点分析2008—2015年褐家鼠湛江种群和哈尔滨种群D-loop单倍型的年度频率变化特点。结果表明,褐家鼠湛江种群和哈尔滨种群共有32种不同的单倍型,其中有11种单倍型是2个种群共有的,有4种单倍型仅在湛江种群中出现,有17种单倍型仅在哈尔滨种群中出现。褐家鼠湛江种群D-loop区的核苷酸多态性为0.005,有27个变异位点,单倍型多态性为0.695,褐家鼠哈尔滨种群D-loop区的核苷酸多态性比湛江种群略高,为0.008,有35个变异位点,单倍型多态性为 0.793。褐家鼠湛江种群和哈尔滨种群没有经历过暴发性的扩增。褐家鼠湛江、哈尔滨和湖北3个地理种群的D-loop序列之间发生了明显的遗传分化,其中湛江种群和哈尔滨种群之间的分化程度最高,遗传分化系数F_st为0.245。褐家鼠湛江种群和哈尔滨种群的单倍型数目和主单倍型频率都发生明显波动,推测主要原因可能是由于灭鼠剂的大量使用或其他灭鼠活动导致种群出现瓶颈或更替的现象。     

From middens to molecules: phylogeography of the piñon pine,Pinus edulis

Aim Data from packrat middens have established a hypothesized historical biogeography of piñon pine, Pinus edulis, including locations of glacial refugia in the south‐western USA and subsequent migration out of the refugia. In this study, we used molecular techniques to test the glacial refugial hypotheses inferred from packrat (Neotoma) midden data for P. edulis. Location South‐western USA. Methods Two fragments of chloroplast DNA (a portion of the matK gene and a portion of the rbcL gene) for a total of 1045 base pairs were amplified and sequenced for 100 individuals. Thirty‐one populations were sampled throughout the range of P. edulis. Phylogenetic analyses included maximum parsimony and maximum likelihood. Results Very little variation existed among the individuals sampled. Four haplotypes were identified. The inferred ancestral haplotype was the most widespread; it was most common in Texas and New Mexico where, with the exception of one individual, it was the only haplotype found. Arizona and Utah populations were more diverse, with almost half of the populations containing two or more haplotypes. The most derived haplotype was most abundant in Arizona. Main conclusions The distribution of haplotypes is geographically informative. Only one haplotype exists in the south‐eastern portion of the range of P. edulis whereas up to four haplotypes are found in other populations, suggesting one of two hypotheses: either all modern populations are descended from a refugial population in central Arizona, or modern populations are descended from two refugial populations, one in central Arizona and another in Texas–southern New Mexico. Interpreting these data in the light of packrat midden data gives more support for the latter hypothesis.