Species‐level phylogeographical history of the endemic species Calligonum roborovskii and its close relatives in Calligonum section Medusa (Polygonaceae) in arid north‐western China

Quaternary climatic oscillations appear to have influenced the genetic diversity and evolutionary history of arid‐adapted plants. To understand the processes involved and reveal evolutionary relationships, haplotypes were examined from Calligonum roborovskii, an endemic species occurring in the arid zones across the desert regions of north‐western China, and seven other species also from Calligonum section Medusa, including C. gobicum, C. mongolicum and the narrow endemic species C. ebi‐nuricum, C. pumilum, C. taklimakanense, C. trifarium and C. yengisaricum. Forty‐three haplotypes were identified in 422 individuals from 51 natural populations, from variation of two plastid DNA intergenic spacers (rpl32‐trnL and ycf6‐psbM). A high level of total genetic diversity was found across species for which more than two populations were examined, including C. gobicum, C. mongolicum, C. pumilum and C. roborovskii. A distinct isolation‐by‐distance pattern in each of these species was suggested by the Mantel test, indicating that restricted gene flow caused high genetic differentiation among populations. Three haplotypes were shared by two or three species each, but the other 40 haplotypes were species‐specific. The 43 haplotypes split into three major clades, but not species‐specific lineages; most of the Calligonum species were not reciprocally monophyletic, probably due to incomplete lineage sorting or introgression. The identified haplotypes were dated to 1.97 Mya (95% highest posterior density: 2.95–0.99 Mya) and diverged until the late Pleistocene, possibly linked to aridification and enlargement of deserts caused by climate changes. Variation of desert habitats during the Pleistocene might play a key role in causing the divergence.     

Phylogeography of Lagochilus ilicifolius (Lamiaceae) in relation to Quaternary climatic oscillation and aridification in northern China

Using two chloroplast DNA intergenic spacers, a phylogeographical study of Lagochilus ilicifolius with 168 individuals from 14 populations was performed to assess geographical patterns and genetic variation in relation to past climate change. Population structure and history were inferred by nested clade analysis, neutrality tests, divergence time estimates, and mismatch distribution analysis. Twelve haplotypes were identified. Genetic differentiation among populations within groups was also elevated (F_SC = 0.772), suggesting restricted gene flow among populations within groups. Divergence time was at approximately 0.71 Ma, consistent with aridification and desert expansion during the middle Pleistocene. The Helan Mountains contained a large proportion of the haplotypes, which implies that the region may be the center of diversification for the species, whereas the Loess Plateau as a dispersal corridor for postglacial re-colonization northward. Climatic oscillation, aridification, and desert expansion within these regions have molded current distribution and biodiversity of L. ilicifolius.     

Chloroplast phylogeographic patterns of Calligonum sect. Pterococcus (Polygonaceae) in arid Northwest China

To understand the impacts of past climatic change and geological events on the evolutionary history of Calligonum sect. Pterococcus, including C. aphyllum, C. rubicundum and C. leucocladum, a total of 128 individuals from 14 populations, mainly from arid Northwest China, were sampled. Two cpDNA intergenic spacer regions (rpl32‐trnL and ycf6‐psbM) were sequenced and 11 haplotypes were identified. Levels of genetic differentiation between populations was low in C. rubicundum (F_ST = 0.54317, p < 0.001) and C. aphyllum (F_ST = 0.55795), while much higher in C. leucocladum (F_ST = 0.95800, p < 0.001), possibly as an effct of differences in geographic distributions and habitats. Analysis of molecular variance (AMOVA) revealed that most of the total genetic variations occurred among species (72.97%). Among eleven identified haplotypes, only H1 and H2 were shared between C. aphyllum and C. rubicundum, while nine were private for one of the three species. The eleven identified haplotypes were divided into two major clades, but they did not yield three species‐specific lineages. Calligonum sect. Pterococcus therefore not appeared reciprocally monophyletic, more likely due to incomplete lineage sorting than hybridization. Mismatch distribution analysis suggested that only C. aphyllum has experienced recent demographic expansion. Divergence time among the 11 haplotypes was estimated at between 2.84 Ma and 0.06 Ma. Within the two clades, haplotype divergence began in early Pleistocene and mainly occurred during the middle to late Pleistocene and was most likely triggered by Quaternary climatic oscillations and increasing aridity of the region.     

Intraspecific molecular phylogeny and phylogeography of the Meriones meridianus (Rodentia: Cricetidae) complex in northern China reflect the processes of desertification and the Tianshan Mountains uplift

The phylogeographical patterns and demographic history of mitochondrial DNA (cytochrome b, N = 327; D‐loop, N = 252) and nuclear DNA (IRBP gene, N = 235) haplotypes were studied for the Meriones meridianus complex in northern China, a desert‐dwelling gerbil species complex. The phylogenetic analyses, which were performed on the separate and combined (mitochondrial + nuclear) datasets, revealed two divergent clades (Clade A and Clade B) corresponding to distinct geographical regions. Clade A contained the haplotypes found mostly in individuals from the Tianshan Mountains area. Clade B contained haplotypes from populations located in other deserts in northern China. The divergence times indicated that the history of the M. meridianus complex was influenced by the uplift of the Tianshan Mountains and climate‐induced habitat fluctuations. In the Pleistocene, the expansion of forests and grasslands during interglacial period led to the isolation of the M. meridianus complex, which preferred to inhabit deserts. Hence, long geological isolation and the M. meridianus complex adaptation to local ecological conditions led to its genetic divergence. Clade A had long‐lasting demographic stability, most likely because the populations of this clade remained in a stable desert environment for a long time. However, the extension of other deserts and disappearance of palaeolakes during the last glacial period resulted in demographic expansion of Clade B. Furthermore, our genetic data indicated that two subspecies may exist within the M. meridianus complex. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110 , 362–383.     

Genetic structure of Salsola junatovii (Chenopodiaceae) in the northern edge of the Taklimakan desert and conservational implications

Salsola junatovii Botsch. (Chenopodiaceae), a desert shrub endemic to the northern edge of the Taklimakan desert in China, is listed as an endangered species because of its narrow habitat requirements and a decline in numbers of individuals. In order to provide information on conservation strategies, we investigated the genetic patterns and evolutionary history of this species. Two cpDNA spacers, (psbA-trnH and psbK-psbI), and the nrITS sequence were employed in 116 individuals from 15 populations across the whole range of the species. A total of 7 haplotypes and 12 ribotypes were detected. Based on separate BEAST analyses, we suggest that intraspecific differentiation of S. junatovii occurred at 2.284–0.264 Ma for cpDNA haplotypes, and 1.304–0.1 Ma for nrITS ribotypes, both in the Pleistocene. Only mismatch analysis and neutrality tests for nrITS sequences over all populations were indicative of demographic expansion. Inconsistency between results from chloroplast and nuclear sequences was detected, apparently due to a stronger seed-mediated gene flow among populations than pollen-mediated flow under the severe desert conditions. Based on genetic structure results, one population, KMS, with the highest number of unique genotypes, should be a conservation priority; another nine populations containing unique genotypes were suggested for ex situ conservation.     

A range wide geographic pattern of genetic diversity and population structure of Hexinia polydichotoma (Asteraceae) in Tarim Basin and adjacent areas

In order to investigate the genetic diversity and influence of climate oscillations on evolutionary processes of organisms in northwest China, we selected Hexinia polydichotoma, a species endemic to China, and examined the geographic pattern of genetic variation in its entire cover range, Tarim Basin and adjacent areas. In the study, 22 chloroplast (cp) haplotypes were identified based on two cpDNA sequences (trnH–psbA and ycf6–psbM), and five ITS sequence variations were found. Shown in the cp haplotype network, the two common cp haplotypes mainly distributed along the northern and southern rims of the basin respectively and intersected in the center of the basin, whereas in the ITS network, ITS genotype 1 was widespread across the whole distribution area, and rare genotypes were concentrated in the western rim of the basin. Genetic variation primarily occurred among populations and SAMOVA groups. Fragmented desert habitat may have caused gene flow barrier among populations or groups far from each other, leading to significant genetic differentiation at these levels. It appears that expansion and contraction cycles of river systems and oases during the middle Pleistocene was the source of the observed fragmentation. Geographic range expansion in H. polydichotoma was supported by the significant value for Tajima's D, Fu's F_S, and by a unimodal mismatch distribution. It was possible that the enlargement of the Taklimakan Desert during the middle Pleistocene may have provided appropriate conditions for the range dispersal. We identified western rim of the basin as the center of genetic diversity of H. polydichotoma based on the present dataset.     

Phylogeographic and demographic effects of Quaternary climate oscillations in Hexinia polydichotoma (Asteraceae) in Tarim Basin and adjacent areas

In order to investigate the genetic diversity and influence of climate oscillations on evolutionary processes of organisms in Northwest China, we selected Hexinia polydichotoma, a species endemic to China, and examined the phylogeographic structure and historical factors that influenced the evolutionary history of this species in its entire cover range, Tarim Basin and adjacent areas. In the study, 17 haplotypes were identified in H. polydichotoma on the basis of two chloroplast DNA sequences (trnH–psbA and ycf6–psbM). Shown in the network, the two common haplotypes, A and D, respectively, mainly distribute along the northern and southern rims of the basin. The analyses of molecular variance analysis suggest that genetic variation primarily occurs among populations, and all populations were subdivided into five groups by SAMOVA. Geographic range expansion along the southern and northern rims of the basin was supported by the significant value for Tajima’s D and by the unimodal mismatch distribution. It is possible that during the interglacial period of the middle Pleistocene, a large amount of snow and glacial ice melted from the mountains surrounding Tarim Basin. This increased water, the expanding desert, and the dispersal ability of H. polydichotoma were important factors driving not only geographic range expansion, but also the current phylogeographic structure of this species. It is possible that during the middle Pleistocene, the climatic fluctuations resulted in expansion and contraction cycles of river systems and oases, and may consequently have caused population fragmentation.     

Remarkable genetic divergence of Gymnodiptychus dybowskii between south and north of Tianshan Mountain in northwest China

Gymnodiptychus dybowskii is endemic to Xinjiang, China and has been locally listed as protected animals. To investigate its genetic diversity and structure, specimens were collected from six localities in Yili River system and Kaidu River. Fragments of 1092bp Cyt b gene were sequenced for 116 individuals. A total of 21 haplotypes were found in all samples, and no haplotype was shared between Yili River system and Kaidu River population. Sequence comparisons revealed 123 variable sites, with eight singleton sites and 115 parsimony informative sites. For all the populations examined, the haplotype diversity (h) was 0.8298 ± 0.0226, nucleotide diversity (π) was 0.2521 ± 0.1202, and average number of pairwise nucleotide differences (k) was 275.3369 ± 118.5660. AMOVA analysis showed that the differences were significant for total populations except for Yili River system populations. The pairwise F_st values revealed same conclusion with AMOVA analysis: Kaidu River population was divergent from Yili River system populations. The genetic distance between two groups was 0.108 and the divergence time was estimated at 5.4–6.6 Ma, the uplift of Tianshan Mountain might have separated them and resulted in the genetic differentiation. The neutrality test and mismatch analysis indicated that both two groups of G. dybowskii had went through population expansion, the expansion time of Yili River system and Kaidu River population was estimated at 0.5859–0.7146 Ma and 0.5151–0.6282 Ma, respectively. The climate changes of Qinghai-Tibetan Plateau might have influenced the demographic history of G. dybowskii.     

Phylogeography of the temperate tree species Quercus acutissima in China: Inferences from chloroplast DNA variations

Quercus acutissima is one of the most widespread temperate deciduous tree species in China. To study its phylogeographical pattern and demographic history, three chloroplast DNA fragments (atpB-rbcL, psbA-trnH and trnS-trnG) from 401 individuals representing 30 populations were sequenced. A total of 19 haplotypes were identified, and these showed a weak phylogeographical structure (N_ST = 0.689 > G_ST = 0.630, P > 0.05) at the species level. The Q. acutissima population harboured a high level of genetic diversity (H_T = 0.791), and the genetic variation mainly resided among populations (59.54%). The unimodal mismatch distribution and significantly negative Fu's F_S value indicate that the Q. acutissima population experienced rapid range expansion, which probably occurred between 0.37 and 0.12 Ma. Molecular phylogeography and ecological niche modelling (ENM) data suggest the existence of multiple localized glacial refugia in central China (e.g., the Qinling, Dabashan and Dabieshan mountain ranges) and southwestern China (Yunnan-Guizhou Plateau and its adjacent regions) during the Quaternary glaciations. Our study showed that geographical heterogeneity and climate changes may have shaped the genetic structure and phylogeographical pattern of this tree species in China.     

Molecular evidence for Pleistocene glacial cycles driving diversification of a North American desert spider, Agelenopsis aperta

The influence of historical climatic vs. geological changes on species diversification patterns was investigated in a widely distributed North American desert spider, Agelenopsis aperta (Araneae: Agelenidae), with particular reference to Pleistocene glacial cycles and earlier patterns of mountain building. Levels of sequence divergence obtained from the mitochondrial gene, cytochrome oxidase I, dated to the Pleistocene, eliminating Rocky Mountain orogeny as a cause of diversification, as orogeny ended 4 million years ago. The results of phylogenetic and network analyses showed the presence of three geographically defined clades, which were consistent with the presence of at least three glacial refugia: (i) east of the Rocky Mountains; (ii) between the Rocky Mountains and Sierra Nevadas; and (iii) west of the Sierra Nevadas. In addition, populations within the Rocky Mountains exhibited significantly lower genetic diversity than populations east of the Rocky Mountains and the haplotypes found within the Rockies were a subset of eastern haplotypes. These patterns suggest that a post-Pleistocene range expansion occurred out of an eastern glacial refugium into the Rocky Mountains. Examination of phylogeographical studies of other North American desert taxa indicated that mountain building explained diversification patterns more effectively for some taxa but Pleistocene climate change was more important for others, including A. aperta.     

Genetic diversity and population differentiation of Capparis spinosa (Capparaceae) in Northwestern China

To explore the influence of vicariance on differentiation patterns of taxa in arid regions, we systematically assessed the genetic diversity and variability of Capparis spinosa, a typical xerophyte that is widely distributed in the Tianshan Mountains and adjacent areas. In total, 300 individuals from 25 populations were sampled, and 14 haplotypes were identified using two cpDNA sequences (rpS12-rpL20 and ndhF). A high level of total genetic diversity (H_T = 0.859) was detected, and this was attributed to the extensive distribution range, which included numerous large populations. The SAMOVA results suggested that the 25 populations were clustered into 4 major geographical groups; a similar divergence trend was found by constructing a BEAST phylogenetic tree and a network diagram. The AMOVA results revealed that significant genetic differentiation occurred among groups. Our results indicated a considerable correlation between genetic divergence and geographical distribution. Isolation due to complex mountain and desert geography might limit gene exchange between disjunct populations, resulting in high differentiation between geographical groups.     

Genetic evidence and species distribution modelling reveal the response of Larix sibirica and its related species to Quaternary climatic and ancient historical events

Larix Mill., a boreal forest genus, distributes in Eurasia and North America. Larix sibirica and some related species of this genus have wide distributions in the northeastern, northwestern and southwestern regions of China. To elucidate the response of populations to different regional historical events, we investigated the spatial genetic structures and the divergence time of genotypes in these Larix species based on two cpDNA fragments. To approach this question, we have also modelled the change of glacial-interglacial distributions for populations in different regions using the species distribution modelling. The results of analyses of molecular variance indicate genetic variation mainly occurred between three isolated regions (74.4%). Distributions of Larix have possibly undergone contractions during the Last Glacial Maximum (LGM) period in comparison with their present-day distributions in all three regions. Recent rapid uplift of Qinghai-Tibet Plateau likely promoted the high level of genetic diversity (π = 0.69 × 10⁻³, Hd = 0.711) and Plio-Pleistocene (1.9–2.7 Ma) lineage divergence in the southwestern region. Founder effect and genetic drift accompanied by postglacial colonizations are expected to shape one dominated haplotype and southern LGM distribution in the northeastern region. Quaternary Glacial period is also speculated to trigger low level of genetic diversity (π = 0.09 × 10⁻³, Hd = 0.133) in the northwestern region. Disjunction of L. sibirica and lineage divergence at about 1.6 Ma in the Altai Mountains and the eastern Tianshan Mountains are related to the enhanced aridity in the Asian interior and the desert expansion in Dzungarian Basin during the early Pleistocene.     

Lack of mitochondrial DNA structure in Balkan donkey is consistent with a quick spread of the species after domestication

A total of 132 mtDNA sequences from 10 Balkan donkey populations were analysed to ascertain their regional genetic structure and to contribute to the knowledge of the spreading of the species after domestication. The Balkan donkey sequences were compared with those from 40 Burkina Faso donkeys as an African outgroup to account for possible local Balkan scenarios. The 172 sequences gave 62 different haplotypes (55 in Balkan donkey). Virtually all the analysed populations had haplotypes assigned to either Clade 1 or Clade 2 even though the relative proportion of Clade 1 or 2 haplotypes differed across populations. Geographical maps constructed using factors computed via principal component analysis showed that the Balkan donkey populations are not spatially structured. AMOVA confirmed a lack of genetic structure in Balkan donkey mtDNA. Balkan populations were poorly differentiated (Φ_ST = 0.071). Differentiation between the Balkan donkey and the African outgroup also was low. The lack of correspondence between geographical areas and maternal genetic structure is consistent with the hypothesis suggesting a very quick spread of the species after domestication. The current research illustrates the difficulties to trace routes of expansion in donkey, as the species has no geographical structure.     

基于线粒体DNA控制区的斑翅草螽不同地理种群遗传分化研究

采用PCR扩增结合DNA克隆测序技术,分析了斑翅草螽Conocephalus maculates 9个地理种群mtDNA控制区序列的变异及遗传多样性。切除侧翼RNA基因序列后,最终获得的斑翅草螽mtDNA控制区比对后全长为676 bp,平均碱基组成T(37.8%),C(11.7%),A(41.3%)和G(9.1%)。共检测到98个可变位点,占总位点数的14.5%,其中,9处碱基插入/缺失,74处转换(40个T/C,34个A/G),50处颠换(18个A/T,11个T/G,15个A/C,6个C/G)。共定义46个单倍型,其中,4个为种群间共享单倍型(H02、H05、H08和H10),其余42个为各种群独有单倍型,包括6个种群内共享单倍型(H09、H11、H15、H18、H26和H38)。单倍型总数占实验个体总数的69.7%,除四川峨眉山外,其余种群单倍型百分比均﹥50%。通过两两地理种群间的FST值差异显著性检验,将这些群体分为4组,分别为SC+CQ,GX+FLB+HN+YN,XZ和HB。以长瓣草螽C.gladiatus、峨眉草螽C.emeiensis、悦鸣草螽C.melaenus、竹草螽C.bambusanus为外群,构建的斑翅草螽mtDNA控制区单倍型NJ法系统树形成3个自举支持度较高的分支,其中,分支A由28种单倍体组成,包括本研究中除四川峨眉山(SC)和重庆万州(CQ)以外的7个种群;分支B由12种单倍体组成,包含除菲律宾拉乌尼翁(FLB)和江西南昌(JX)以外的7个种群;分支C由6种单倍型组成,全部来自西藏林芝(XZ)的单倍型。聚类结果表明,斑翅草螽不同地理种群间的遗传分化并不明显,即使是两两群体间FST值差异显著的群体,也未能形成完全独立的分支。     

From Asia to Europe across Siberia: phylogeography of the Siberian spined loach (Teleostei,Cobitidae)

We assessed the genetic variability of the Siberian spined loach Cobitis melanoleuca across its unusually broad distribution that encompasses areas greatly affected by Pleistocene glaciations. Due to extensive morphological variation among their populations, the taxonomic status of C. melanoleuca is complicated. It is unclear whether C. melanoleuca represents a single taxonomic unit or contains several species or subspecies. Our analyses showed low genetic variability in all populations without any phylogenetic structure. The absence of molecular distinctiveness indicates the conspecificity of all C. melanoleuca populations. Only a few common haplotypes shared by East Asian, Siberian and European populations were found at high frequency in the nuclear genes analysed. At the mitochondrial level, Siberian populations shared haplotypes with populations located at both extremes of the species’ range suggesting central populations as a source of current mitochondrial variability. Unimodal mismatch distributions and significant values from neutrality tests support a recent expansion of C. melanoleuca. Our time estimates suggest a postglacial colonisation of European waters around 1.0 MYA, indicating that C. melanoleuca may represent the last cobitid immigrant in Europe that used the northern route across Siberia to expand its range.     

The effect of past climatic oscillations on spatial genetic structure of Atraphaxis manshurica (Polygonoideae) in the Horqin sandlands,northern China

To investigate the effect of past climate oscillations on the genetic diversity and demographic history of organisms in the Horqin sandlands, at the margin of the monsoonal zone in northern China. We selected Atraphaxis manshurica as a model and examined the spatial genetic structure in populations of this species across most of its covered range. Five haplotypes were identified on the basis of two chloroplast DNA sequences (psbK-psbI and psbB-psbH) in 127 individuals from 11 natural populations. The total genetic diversity (H_T = 0.437) of A. manshurica is low compared with species to the west in arid northwestern China, outside of monsoonal influence. Analysis of molecular variance (AMOVA) indicates that genetic differentiation primarily occurred among two geographical groups distinguished by spatial analysis of molecular variance (SAMOVA). The BEAST molecular clock approach revealed that genetic divergence of the species mostly occurred in middle-late Pleistocene, in accordance with the periodic glacial periods and accompanying monsoonal changes. From species distribution modeling (SDM), we found that the species experienced range contraction and southward retreat during the Last Glacial Maximum (LGM). The current genetic pattern and demographic history may have thus been shaped by glacial–interglacial cycles and changes of habitat since the middle Pleistocene.     

Genetic diversity and maternal origin of Northeast African and South American donkey populations

To investigate the mtDNA variation and origin of maternal lineages in South American donkeys and to reassess the domestication of donkeys in northeast Africa, we analyzed sequences (489 bp of the D‐loop) from 323 domestic donkeys sampled from Peru, Brazil, Ethiopia and Egypt. Altogether, the 323 sequences displayed 53 different haplotypes (45 in Ethiopia, 14 in Egypt, eight in Peru and six in Brazil). Among the four populations, Egyptian donkeys possessed the highest haplotype diversity (0.910 ± 0.032), followed by Brazilian donkeys (0.879 ± 0.060). The Clade I haplotypes dominated in Peruvian donkeys (65%), whereas Clade II haplotypes dominated in Brazilian donkeys (67%). Estimates of F_ST values showed a high genetic differentiation between Peruvian and Brazilian donkey populations (F_ST = 0.4066), which could be explained by the complex introduction history of South American donkeys. Phylogeographic analysis indicates that northeast Africa could be the most probable domestication center for Clade I donkeys. Analysis of molecular variance confirmed a weak genetic structure in domestic donkey populations among four continents (Europe, Asia, Africa and South America).     

Range expansions in the flightless longhorn cactus beetles, Moneilema gigas and Moneilema armatum, in response to Pleistocene climate changes

Pollen cores and plant and animal fossils suggest that global climate changes at the end of the last glacial period caused range expansions in organisms indigenous to the North American desert regions, but this suggestion has rarely been investigated from a population genetic perspective. In order to investigate the impact of Pleistocene climate changes and glacial/interglacial cycling on the distribution and population structure of animals in North American desert communities, biogeographical patterns in the flightless, warm-desert cactus beetles, Moneilema gigas and Moneilema armatum, were examined using mitochondrial DNA (mtDNA) sequence data from the cytochrome oxidase I (COI) gene. Gene tree relationships between haplotypes were inferred using parsimony, maximum-likelihood, and Bayesian analysis. Nested clade analysis and coalescent modelling using the programs mdiv and fluctuate were used to identify demographically independent populations, and to test the hypothesis that Pleistocene climate changes caused recent range expansions in these species. A sign test was used to evaluate the probability of observing concerted population growth across multiple, independent populations. The phylogeographical and nested clade analyses reveal a history of northward expansion in both of these species, as well as a history of past range fragmentation, followed by expansion from refugia. The coalescent analyses provide highly significant evidence for independent range expansions from multiple refugia, but also identify biogeographical patterns that predate the most recent glacial period. The results indicate that widespread desert environments are more ancient than has been suggested in the past.     

Genetic diversity and population structure of Schizopygopsis younghusbandi Regan in the Yarlung Tsangpo River inferred from mitochondrial DNA sequence analysis

The genetic diversity and population structure of Schizopygopsis younghusbandi from six sites in the middle reach of the Yarlung Tsangpo River were examined based on mitochondrial DNA cytochrome b (Cyt b) gene and control region (D-loop) sequences. 1141 bp of complete Cyt b sequences and 737 bp of partial D-loop sequences for 153 individuals of S. younghusbandi were obtained by using PCR amplification and sequencing. The results showed that S. younghusbandi populations had high haplotype diversity and low nucleotide diversity, and the population genetic diversity of this species was at a low level. Analysis of molecular variance revealed that most genetic variation occurred within populations, and that genetic differentiation was at low or moderate levels. The network of haplotypes based on Cyt b gene showed that there were two groups within the examined populations. Neutrality tests and mismatch distribution analysis suggested that this species might have undergone a population expansion and the expansion time was estimated to be 0.25–0.46 Ma BP. All these results would be crucial to establish scientific strategies for effective conservation and sustainable exploitation of wild resources of S. younghusbandi.     

Low genetic diversity in the vulnerable Goitred Gazelle,Gazella subgutturosa (Cetartiodactyla: Bovidae), in Iran: potential genetic consequence of recent population declines

The Goitered Gazelle, Gazella subgutturosa, is the most widespread gazelle species in the Middle East and central Asia inhabiting desert and semi-desert habitats. Today it is threatened and its geographic range and population size have experienced significant decline in the last decades. In Iran, the remnant populations are confined to fragmented habitats. We aimed to characterise genetic diversity and phylogenetic status of the populations of Goitered Gazelle in Central Iran and to evaluate the potential effect of a historic population bottleneck on the genetic variation of today’s population. We used noninvasive sampling to uncover structure and level of genetic variation in a fragment of the cytochrome-b gene from 170 samples. Genealogical analyses were performed using HKY+I model and phylogenetic trees reconstructed using Bayesian inference and maximum likelihood. We found extremely low levels of genetic variation, with altogether only five haplotypes in samples from different populations. Overall haplotype diversity was 0.081 and nucleotide diversity 0.0003. The mean observed mismatch between any two sequences was 0.093 with the largest peak for small numbers. The mismatch distribution fit the model of population expansion and suggested that gazelles had experienced a sudden expansion. An unrooted median-joining network analysis of mtDNA haplotypes showed a star-like structure which few mutations steps separating the haplotypes from other regions. Our findings strengthen the urgency of preserving the species’ genetic diversity to prevent local extinction.