High mitochondrial differentiation levels between wild and domestic Bactrian camels: a basis for rapid detection of maternal hybridization

Hybridization between wild species and their domestic congeners often threatens the gene pool of the wild species. The last wild Bactrian camel (Camelus ferus) populations in Mongolia and China are examples of populations facing such a hybridization threat. To address this key issue in the conservation of wild camels, we analysed wild, hybrid and domestic Bactrian camels (Camelus bactrianus) originating from Mongolia, China and Austria. Through screening of an 804‐base‐pair mitochondrial fragment, we identified eight mitochondrial haplotypes and found high sequence divergence (1.9%) between C. ferus and C. bactrianus. On the basis of a mitochondrial DNA sequence fixed difference, we developed a diagnostic PCR restriction fragment length polymorphism (PCR‐RFLP) assay to differentiate between wild and domestic camel samples. We applied the assay to 81 individuals and confirmed the origin of all samples including five hybrids with known maternal ancestry. The PCR‐RFLP system was effective for both traditional (blood, skin) and non‐invasive samples (faeces, hair), as well as for museum specimens. Our results demonstrate high levels of mitochondrial differentiation between wild and domestic Bactrian camels and that maternal hybridization can be detected by a rapid and reliable PCR‐RFLP system.     

Old World camels in a modern world – a balancing act between conservation and genetic improvement

Old World camels have served humans in cross‐continental caravans, transporting people and goods, connecting different cultures and providing milk, meat, wool and draught since their domestication around 3000–6000 years ago. In a world of modern transport and fast connectivity, these beasts of burden seem to be out‐dated. However, a growing demand for sustainable milk and meat production, especially in countries affected by climate change and increasing desertification, brings dromedaries (Camelus dromedarius) and Bactrian camels (Camelus bactrianus) back onstage and into the focus of animal breeders and scientists. In this review on the molecular genetics of these economically important species we give an overview about the evolutionary history, domestication and dispersal of Old World camels, whereas highlighting the need for conservation of wild two‐humped camels (Camelus ferus) as an evolutionarily unique and highly endangered species. We provide cutting‐edge information on the current molecular resources and on‐going sequencing projects. We cannot emphasise enough the importance of balancing the need for improving camel production traits with maintaining the genetic diversity in two domestic species with specific physiological adaptation to a desert environment.     

SNP discovery and population structure analysis in Lassi and Marecha camel breeds using a genotyping by sequencing method

Pakistani camels have been classified socio-geographically into 20 breeds, but they have not yet been subjected to substantial selective pressures and the genetic basis for these breeds is not understood. However, it should be possible to distinguish them by use of molecular data. This study investigated the genetic diversity and population structure within and between two major Pakistani camel breeds, Marecha and Lassi. As no SNP array is currently available, we first identified 63 619 SNPs using a genotyping by sequencing approach. After quality control, a panel of 36 926 SNPs was used in the analysis. Population structure was investigated with a principal coordinate analysis as well as a cluster analysis using NetView , and multilocus heterozygosity analysis to explore between- and within-breed genetic variation. In addition, between-breed variation was explored using the fixation index, F_ST. We also compared relationship matrices computed using the VanRaden SNP-based method and a method developed specifically for genotyping by sequencing data. Among the two camel breeds, Lassi showed a lower level of genetic diversity whereas Marecha showed a higher level. As a genotyping platform has not yet been developed for the camel, the SNPs discovered in this study will be useful in future genetic studies in camels.     

Monophyletic origin of domestic bactrian camel (Camelus bactrianus) and its evolutionary relationship with the extant wild camel (Camelus bactrianus ferus)

The evolutionary relationship between the domestic bactrian camel and the extant wild two-humped camel and the factual origin of the domestic bactrian camel remain elusive. We determined the sequence of mitochondrial cytb gene from 21 camel samples, including 18 domestic camels (three Camelus bactrianus xinjiang , three Camelus bactrianus sunite , three Camelus bactrianus alashan , three Camelus bactrianus red , three Camelus bactrianus brown and three Camelus bactrianus normal ) and three wild camels ( Camelus bactrianus ferus ). Our phylogenetic analyses revealed that the extant wild two-humped camel may not share a common ancestor with the domestic bactrian camel and they are not the same subspecies at least in their maternal origins. Molecular clock analysis based on complete mitochondrial genome sequences indicated that the sub-speciation of the two lineages had begun in the early Pleistocene, about 0.7 million years ago. According to the archaeological dating of the earliest known two-humped camel domestication (5000–6000 years ago), we could conclude that the extant wild camel is a separate lineage but not the direct progenitor of the domestic bactrian camel. Further phylogenetic analysis suggested that the bactrian camel appeared monophyletic in evolutionary origin and that the domestic bactrian camel could originate from a single wild population. The data presented here show how conservation strategies should be implemented to protect the critically endangered wild camel, as it is the last extant form of the wild tribe Camelina.     

New Pliocene remains of Camelus grattardi (Mammalia,Camelidae) from the Shungura Formation,Lower Omo Valley,Ethiopia, and the evolution of African camels

Abstract

Camels are exceptionally rare in the Plio-Pleistocene fossil record of Africa, hindering attempts to understand the evolution of this family on the continent. Here we describe recently collected camel specimens from the Shungura Formation, Lower Omo Valley, Ethiopia, and attribute these remains to Camelus grattardi. The new specimens date to the late Pliocene (~3 to 2.6 Ma) and consist of three upper molars, one upper premolar, and two proximal metatarsals. The dental specimens confirm this species’ small P4 relative to its molars, a trait that differs significantly from all extant and fossil Old World camels. The metatarsals indicate that C. grattardi was similar in size to the living Bactrian camel, C. bactrianus. Phylogenetically, we find no suitable ancestor, sister, or descendant of the eastern African fossil camel, which suggests greater lineage diversity in Plio-Pleistocene Camelus than previously recognised. Microwear analyses suggest that C. grattardi was likely a mixed-feeder preferring browse, which is consistent with carbon isotopes of enamel from the Turkana Basin. A review of the fossil record of African camels suggests no clear paleoenvironmental association, as fossil camels occur in a range of environments from dry savannas with no permanent water bodies to closed woodlands along the paleo-Omo River.     

Genetic diversity and population structure of Mongolian domestic Bactrian camels (Camelus bactrianus)

The tradition of animal husbandry in the context of a nomadic lifestyle has been of great significance in the Mongolian society. Both Bactrian camels and horses have been invaluable for the survival and development of human activities in the harsh arid environment of the Mongolian steppe. As camels offer unique and sustainable opportunities for livestock production in marginal agro‐ecological zones, we investigated the current genetic diversity of three local Mongolian camel breeds and compared their levels of variation with common native Mongolian camels distributed throughout the country. Based on mitochondrial and nuclear markers, we found levels of genetic diversity in Mongolian populations similar to that reported for Chinese Bactrian camels and for dromedaries. Little differentiation was detected between single breeds, except for a small group originating from the northwestern Mongolian Altai. We found neither high inbreeding levels in the different breeds nor evidence for a population decline. Although the Mongolian camel census size has severely declined over the past 20 years, our analyses suggest that there still exists a stable population with adequate genetic variation for continued sustainable utilization.     

Results of <Emphasis Type="Italic">Camelus dromedarius</Emphasis> and <Emphasis Type="Italic">Camelus bactrianus</Emphasis> Genotyping by Alpha-S<Subscript>1</Subscript>-Casein,Kappa-Casein Loci,and DNA Fingerprinting

Genotyping of Kazakh camels Camelus dromedarius (milk breed) (n = 18) and Camelus bactrianus (meat breed) (n = 18) by alpha-S₁-casein (αs₁-CN) and kappa-casein (κ-CN) loci was conducted using the PCR–RFLP analysis method. A new pair of primers was suggested for the amplification of the CSN3 gene fragment with subsequent cleavage of the reaction products by AluI restriction endonuclease in order to identify the gene genetic variants. DNA polymorphism was detected only for the kappa-casein locus; no genetic polymorphism for alpha-S₁-casein gene was found in the studied populations. Analysis of the results of DNA fingerprinting demonstrated that the band sharing (BS) coefficient between the groups was low enough (0.13), and the genetic distance (D) between Dromedary and Bactrian breeds was 0.305. The results of genotyping of Bactrian and Dromedary Kazakh camel breeds by alpha-S₁-casein, kappa-casein loci, and DNA fingerprinting indicate that the Dromedary breed female camels are more polymorphic as compared with Bactrian.     

Population genetic analysis of the domestic Bactrian camel in China by RAD‐seq

Restriction site‐associated DNA sequencing (RAD‐seq) is one of the most effective high‐throughput sequencing technologies for SNP development and utilization and has been applied to studying the origin and evolution of various species. The domestic Bactrian camels play an important role in economic trade and cultural construction. They are precious species resources and indispensable animals in China's agricultural production. Recently, the rapid development of modern transportation and agriculture, and the deterioration of the environment have led to a sharp decline in the number of camels. Although there have been some reports on the evolution history of the domestic Bactrian camel in China, the origin, evolutionary relationship, and genetic diversity of the camels are unclear due to the limitations of sample size and sequencing technology. Therefore, 47 samples of seven domestic Bactrian camel species from four regions (Inner Mongolia, Gansu, Qinghai, and Xinjiang) were prepared for RAD‐seq analysis to study the evolutionary relationship and genetic diversity. In addition, seven domestic Bactrian camel species are located in different ecological zones, forming different characteristics and having potential development value. A total of 6,487,849 SNPs were genotyped. On the one hand, the filtered SNP information was used to conduct polymorphism mapping construction, LD attenuation analysis, and nucleotide diversity analysis. The results showed that the number of SNPs in Dongjiang camel was the highest, the LD coefficient decayed the fastest, and the nucleotide diversity was the highest. It indicates that Dongjiang camel has the highest genetic diversity. On the other hand, the filtered SNPs information was used to construct the phylogenetic tree, and F_ST analysis, inbreeding coefficient analysis, principal component analysis, and population structure analysis were carried out. The results showed that Nanjiang camel and Beijiang camels grouped together, and the other five Bactrian camel populations gathered into another branch. It may be because the mountains in the northern part of Xinjiang and the desert in the middle isolate the two groups from the other five groups.     

Mixed ancestry from wild and domestic lineages contributes to the rapid expansion of invasive feral swine

Invasive alien species are a significant threat to both economic and ecological systems. Identifying the processes that give rise to invasive populations is essential for implementing effective control strategies. We conducted an ancestry analysis of invasive feral swine (Sus scrofa, Linnaeus, 1758), a highly destructive ungulate that is widely distributed throughout the contiguous United States, to describe introduction pathways, sources of newly emergent populations and processes contributing to an ongoing invasion. Comparisons of high‐density single nucleotide polymorphism genotypes for 6,566 invasive feral swine to a comprehensive reference set of S. scrofa revealed that the vast majority of feral swine were of mixed ancestry, with dominant genetic associations to Western heritage breeds of domestic pig and European populations of wild boar. Further, the rapid expansion of invasive feral swine over the past 30 years was attributable to secondary introductions from established populations of admixed ancestry as opposed to direct introductions of domestic breeds or wild boar. Spatially widespread genetic associations of invasive feral swine to European wild boar deviated strongly from historical S. scrofa introduction pressure, which was largely restricted to domestic pigs with infrequent, localized wild boar releases. The deviation between historical introduction pressure and contemporary genetic ancestry suggests wild boar‐hybridization may contribute to differential fitness in the environment and heightened invasive potential for individuals of admixed domestic pig–wild boar ancestry.     

Comparison of virokine from camel pseudocowpoxvirus (PCPV) with Interleukin 10 of the Dromedary camel (Camelus dromedarius)

Cellular interleukin-10 (IL-10) gene from the peripheral blood mononuclear cells of the healthy Dromedary camel (Camelus dromedarius) and viral IL-10 (vIL-10) from the skin scabs of the Dromedary camels infected with contagious ecthyma (a parapoxviral infection in the camels) were amplified by polymerase chain reaction, cloned and characterized. Sequence analysis revealed that the open reading frame (ORF) of dromedarian camel IL-10 is 537 bp in length, encoding 178 amino acid polypeptide while open reading frame of vIL-10 from camel is 561 bp, encoding 187 amino acid polypeptide. The Dromedary camel IL-10 exhibited 62.6% and 68.5% sequence identity at the nucleotide and amino acid level, respectively, with vIL-10 from camel. Sequence analysis also revealed that the Dromedary camel IL-10 shared 99.4% and 98.3% identity at the nucleotide and amino acid level, respectively, with the Bactrian camel (Camelus bactrianus). But vIL-10 from camel shared 84.7% and 83.4% sequence identity at the nucleotide and amino acid level, respectively, with vIL-10 from reindeer (Rangifer tarandus), which is a ruminant species belonging to the order Artiodactyla. The present study was conducted to evaluate the evolutionary origin of the camel parapoxvirus with parapoxviruses of cattle and sheep and the resultant sequence analysis revealed that camel parapoxvirus is closely related to cattle parapoxvirus than sheep parapoxvirus (Orf virus).     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

野双峰驼各分布区的生存环境差异及评价

世界上野双峰驼(Camelus bactrianus ferus)目前仅分布于中国新疆境内及中蒙边境,即：中国新疆的塔克拉玛干沙漠、阿尔金山北麓、戛顺戈壁、中蒙边境及蒙古国西部的外阿尔泰戈壁。总数共计730～880头。目前已为极度濒危物种。本文对野驼在不同分布区由于生境及人类活动影响所造成的差异加以对比研究并进行评价,为野骆驼保护提供依据。     

Identification and management of a single large population of wild dromedary camels

The dromedary camel (Camelus dromedarius) is a significant invasive species in Australia. It is an unusual pest species that is of large body size with relatively low fecundity compared with other pest species. Camels are highly adapted to the arid regions that characterize a large proportion of Australia and occupy an almost completely undisturbed area of ≥3 million km². They have no history of invasion elsewhere in the world. Despite this, their population has been expanding at approximately 80,000 camels per annum, with the most recent estimate of population size around 1,000,000 individuals. We employed a landscape-genetic approach to evaluate the population structure and molecular ecology of Australian camels. We combined mitochondrial control region sequence (n = 209 animals) with 18 microsatellite markers to profile over 800 adult camels to identify the presence of a single panmictic population. We showed that demographically defined neighborhoods for wild camels are about 200 km; this value was supported by home range estimates. Distances greater than this display no pattern of isolation by distance across the Australian continent. The result is the largest single geographical population so far recorded for an invasive species in Australia. This pattern may be explained by the impressive and near-nomadic dispersal pattern of camels, in combination with an unpredictable environment virtually devoid of barriers to movement and predatory suppression. Although it is technically feasible, the reality is that it would not be economically or politically viable to have continental eradication of wild camels in Australia because of the vast size and movement dynamics of the camel population. As such, we advocate a change away from an expensive solution to an intractable reduction program (that is almost entirely focused on protection of biological refugia) and moves to include cultural, economic, and biodiversity asset protection for the management of this most unorthodox of invasive species. © 2012 The Wildlife Society.     

Signatures of de‐domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution

Autochthonous pig breeds are usually reared in extensive or semi‐extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo‐Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south‐eastern European countries (Kr?kopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro‐geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild‐type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of ‘de‐domestication’ process, and wild resources are challenged by a ‘domestication’ drift. Both need to be further investigated and managed.     

Cloning and Molecular Characterization of HSL and Its Expression Pattern in HPG Axis and Testis during Different Stages in Bactrian Camel

Hormone-sensitive lipase (HSL) is a key enzyme in animal fat metabolism and is involved in the rate-limiting step of catalyzing the decomposition of fat and cholesterol. It also plays an important regulatory role in maintaining seminiferous epithelial structure, androgen synthesis and primordial germ cell differentiation. We previously reported that HSL is involved the synthesis of steroids in Bactrian camels, although it is unclear what role it plays in testicular development. The present study was conducted to characterize the biological function and expression pattern of the HSL gene in the hypothalamic pituitary gonadal (HPG) axis and the development of testis in Bactrian camels. We analyzed cloning of the cDNA sequence of the HSL gene of Bactrian camels by RT-PCR, as well as the structural features of HSL proteins, using bioinformatics software, such as ProtParam, TMHMM, Signal P 4.1, SOPMA and MEGA 7.0. We used qRT-PCR, Western blotting and immunofluorescence staining to clarify the expression pattern of HSL in the HPG axis and testis of two-week-old (2W), two-year-old (2Y), four-year-old (4Y) and six-year-old (6Y) Bactrian camels. According to sequence analysis, the coding sequence (CDS) region of the HSL gene is 648 bp in length and encodes 204 amino acids. According to bioinformatics analysis, the nucleotide and amino acid sequence of Bactrian camel HSL are most similar to those of Camelus pacos and Camelus dromedarius, with the lowest sequence similarity with Mus musculus. In adult Bactrian camel HPG axis tissues, both HSL mRNA and protein expression were significantly higher in the testis than in other tissues (hypothalamus, pituitary and pineal tissues) (p < 0.05). The expression of mRNA in the testis increased with age and was the highest in six-year-old testis (p < 0.01). The protein expression levels of HSL in 2Y and 6Y testis were clearly higher than in 2W and 4Y testis tissues (p < 0.01). Immunofluorescence results indicate that the HSL protein was mainly localized in the germ cells, Sertoli cells and Leydig cells from Bactrian camel testis, and strong positive signals were detected in epididymal epithelial cells, basal cells, spermatocytes and smooth muscle cells, with partially expression in hypothalamic glial cells, pituitary suspensory cells and pineal cells. According to the results of gene ontology (GO) analysis enrichment, HSL indirectly regulates the anabolism of steroid hormones through interactions with various targets. Therefore, we conclude that the HSL gene may be associated with the development and reproduction of Bactrian camels in different stages of maturity, and these results will contribute to further understanding of the regulatory mechanisms of HSL in Bactrian camel reproduction.     

Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle

Complete mitochondrial DNA D‐loop sequences of 1105 individuals were used to assess the diversity of maternal lineages of cattle populations in China. In total, 250 taurine and 88 zebu haplotypes were identified. Five main haplogroups—T1a, T2, T3, T4 and T5—were identified in Bos taurus, whereas Bos indicus harbored two haplogroups—I1 and I2. Our results suggest that the distribution of T1a in Asia was concentrated mainly in the northeast region (northeast China, Korea and Japan); haplogroups T2, T3 and T4 were predominant in Chinese cattle; and T5 was sporadically detected in Mongolian and Pingwu cattle. In contrast to the widespread presence of I1, I2 was distributed only in southwestern China (Yunnan‐Guizhou Plateau and the Tibet Autonomous Region) and Xinjiang Uygur Autonomous Region. This is the first time that all five taurine haplogroups and two zebu haplogroups have been found in Mongolian cattle. In addition, eight individuals in Tibetan cattle carried the Bos grunniens mtDNA type. The high mtDNA diversity (H = 0.904 ± 0.008) and the weak genetic structure among the 57 Chinese cattle breeds/populations are consistent with their complex historical background, migration route and ecological environment.     

Detecting population structure and recent demographic history in endangered livestock breeds: the case of the Italian autochthonous donkeys

Since its domestication, about 5000 years ago, the donkey (Equus asinus) has been extensively used as a work or draft animal in agricultural activities and for the transportation of people and goods. In the last century, technology improvement and growing mechanization strongly affected agriculture and the management and use of this livestock species in the industrialized countries. Nowadays, the use of donkeys for work or transport has almost disappeared, together with the need for mules or hinny breeding. During the last five decades, Italian autochthonous donkey populations suffered from a severe reduction in population size, which led to the extinction of several breeds. At present, eight breeds remain, all classified by FAO as critically endangered or endangered: Asinara, Pantesco, Grigio Siciliano, Romagnolo, Amiatino, Sardo Grigio, Martina Franca, and Ragusano. To evaluate the extant genetic variability of Italian donkeys, we typed 16 microsatellite loci in 258 individuals from these breeds. The results highlighted moderate levels of inbreeding ( F _IS = 0.127) and a significant partition of genetic variation into breeds, as suggested by fixation index ( F _ST = 0.109) and analysis of molecular variance (10.86% of total variation assigned to the between‐breeds level) analyses. This was confirmed by a Bayesian clustering procedure that also highlighted a further partitioning at lower hierarchical levels corresponding to the farms of origin. This evidence suggests that an effective management strategy for Italian donkey populations should focus on breeds as conservation units. However, this requires a synergic management strategy at the farm level to maintain diversity and avoid inbreeding.     

Microscopic and molecular detection of Deraiophoronema evansi (Lewis, 1882) in domestic Bactrian camels (Camelus bactrianus) of Mongolia

Cameline filarosis is an important parasitic disease having an economic impact on the camel industry around the world. However, there has been no study on filarosis in Bactrian camels of Mongolia. Therefore, the aim of the present study was to detect and identify microfilariae of Deraiophoronema evansi (D. evansi) in Bactrian camels from three provinces, located in southern and southwestern Mongolia. Blood samples were obtained from 400 healthy two-humped camels of different ages and both sexes. All blood samples were analysed using a variety of diagnostic techniques. Microfilariae were detected in 30 Bactrian camels (7.5%) by the Knott technique, while 13 Bactrian camels (3.3%) tested positive in a direct smear test. D. evansi was detected in 18 Bactrian camels (4.5%) by PCR assay. Prevalence was shown to be high among Bactrian camels in the age group up to 5 years, while the lowest positive results were obtained for Bactrian camels in the 5–10-year age group and the over 10-year age group. To confirm the morphological identification, D. evansi–COI gene sequences were subjected to phylogenetic analyses. The D. evansi–COI gene sequences from Mongolian two-humped camels were identical to sequences from Iranian one-humped camels and were clustered together with these sequences in the phylogeny. This is the first report of molecular detection and identification of microfilariae of D. evansi in Bactrian camels of Mongolia.     

Population genetic structure and adaptive differentiation of iron walnut Juglans regia subsp. sigillata in southwestern China

Southwestern (SW) China is an area of active tectonism and erosion, yielding a dynamic, deeply eroded landscape that influences the genetic structure of the resident populations of plants and animals. Iron walnut (Juglans regia subsp. sigillata) is a deciduous tree species endemic to this region of China and cultivated there for its edible nuts. We sampled 36 iron walnut populations from locations throughout the species' range in SW China and genotyped a total of 765 individuals at five chloroplast DNA regions and 22 nuclear microsatellite loci. Species distribution models were produced to predict the evolution and historical biogeography of iron walnut and to estimate the impacts of climate oscillations and orographic environments on the species' demography. Our results indicated that J. regia subsp. sigillata had relatively low genetic diversity, high interpopulation genetic differentiation, and asymmetric interpopulation gene flow. Based on DIYABC analysis, we identified two lineages of J. sigillata in southwestern China. The lineages (subpopulations) diverge during the last glacial period (~1.34 Ma). Southwestern China was a glacial refuge during the last glacial period, but increasingly colder and arid climates might have fostered the fragmentation of J. regia subsp. sigillata within this refugium. Finally, we found that recent habitat fragmentation has led to a reduction in population connectivity and increased genetic differentiation by genetic drift in isolated populations. Our results support a conclusion that geological and climatic factors since the Miocene triggered the differentiation, evolutionary origin, and range shifts of J. sigillata in the studied region.