INTRODUCTION
The Old World camels are represented by the Tylopoda suborder which comprises domesticated Bactrian camel (Camelus Bactrianus) and dromedary (Camelus dromedarius), and the only representative wild species remaining, wild Bactrian camel (Camelus ferus), following the extinction of the wild dromedary [1]. Traditionally, the Old World camels have provided a wide range of useful products to Gobi desert rural communities, including meat, milk, and wool, since their domestication around 3,000 to 6,000 years ago. Currently, the high nutritional and medical value of their milk has increased the number of breeders, who have extended the exploitation of these species to other regions outside the Gobi Desert, especially Bactrian camels. The natural habitat of the domestic Bactrian camel are the cold desert areas of Northeast and Central Asia, in contrast to the dromedaries, which are adapted to the semi-arid and hot desert regions of North Africa and East Africa, the Arabian Peninsula and Southwest Asia [2]. The wild Bactrian camel is restricted to few remaining refuge areas, with their range extending to only three locations in China (Taklamakan desert, Gashun Gobi Desert and Arjin Mountains in the Lop Nur Lake region) and one in Mongolia (Great Gobi Strictly Protected Area ‘A’) [3]. Nowadays, the wild Bactrian camel is listed as Critically Endangered [4] and its population is estimated to number from a few hundred to 2,000 individuals [5,6]. The Old World camels show several biological and physiological traits that may be connected with adaptation to such extreme heat and harsh and dry environments, including resistance to hunger and thirst, fluctuating body temperatures [7], tolerance of a high dietary intake of salt, and an immune system producing unique immunoglobulin [8]. Recent genomic studies have revealed positive selection [9] and immune gene loci, which may be related to desert environmental adaptation in this species [10].
Mitochondrial DNA (mtDNA), as a molecular marker, is widely used in the identification of wild ancestors and maternal lines that have contributed to a breed or population [11–14]. Over recent years, nuclear and mtDNA markers have been used intensively to unravel the genetic diversity and phylogeographic structure in the Old World camels [15–17]. Within mtDNA, the control region (CR) has been revealed to evolve five times faster than the coding region [18]. Therefore, the CR is considered suitable for assessing population genetic structure, variability and phylogeny.
With the modernization of society and increasing desertification, the genomic resources from the two domesticated Old World camel species, dromedary (Camelus dromedarius) and Bactrian camel (Camelus bactrianus), and one wild species, wild Bactrian camel (Camelus Ferus), have been seriously affected; this is especially the case for the wild Bactrian camel species, which is on the verge of extinction [4]. Unfortunately, little is known about the genetic diversity and phylogeography of Old World camels. Therefore, we selected mtDNA CR sequences of 182 individuals from domestic and wild Bactrian camels, and dromedaries, to investigate the genetic diversity, phylogeography and demographic dynamics of the different camel populations.
MATERIALS AND METHODS
Data sources
The present study was conducted on 182 sequences which were downloaded from GenBank. These included wild Bactrian camel from Mongolia (n = 27); domestic Bactrian camel from China (n = 57), Mongolia (n = 43), Russia (n = 10), Kazakhstan (n = 6), and Iran (n = 13); and dromedary (n = 26) from Africa (Supplementary Table S1). In accordance with references [15,17,19], we restricted the 809 bp mtDNA fragment nt 15120–15928 from each sequence, and multiple sequence alignments were performed using Mega version 6 [20] against the wild Bactrian camel reference sequence (Genbank accession number NC_009628) for Bactrian camel populations, and the dromedary reference sequence (Genbank accession number NC_009849) for dromedaries.
Genetic diversity estimation
DnaSP v5 was used to determine the level of genetic diversity, including polymorphic sites (S), number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (π), the average number of nucleotide differences (K) and the standard deviations (SD). These parameters were identified for each population and across all populations [21].
Phylogenetic analysis
In order to assess the possible genetic relationship between individuals and populations, a phylogenetic tree was constructed with the neighbor-joining (NJ) algorithm implemented in MEGA6 [20]. The confidence level for each branch was evaluated with 1,000 bootstrap replications. To further infer the genetic relationship between the haplotypes from different camel populations, a median-joining (MJ) network was built using the NETWORK 5.0.0 with equal weight mutations and character states [22].
Neutrality test and demographic dynamics
We calculated the demographic profiles for each geographic region from mismatch distribution patterns [23]. At the same time, the Harpending’s raggedness index “r” [24] statistic was used to evaluate the significance of the deviations of the observed sum of squares differences (SSD) from the simulated model of expansion based on 1,000 coalescent simulations. Furthermore, Fu’s Fs [25] and Tajima’s D [26] statistics were also considered using the infinite sites model in Arlequin v3.5 [27]. Tajima’s D is negative, indicating that there is an excess of low-frequency polymorphisms relative to expectations, and the population size is expanding (both bottlenecks or selective sweep) or purification selection is occurring [28]. Positive Tajima’s D signifies low levels of high frequency polymorphisms, explaining a decrease in population size and/or balancing selection. The same holds for Fu’s Fs test, which is more powerful than Tajima’s D [29].
RESULTS
mtDNA sequence haplotype variation and genetic diversity
Thirty-two haplotypes confirmed by 54 polymorphic sites were identified in the 156 sequences, which included 129 domestic Bactrian camels and 27 wild Bactrian camels. There were 30 haplotypes determined in the domestic Bactrian camel population and two in the wild Bactrian camel population. Meanwhile, 14 haplotypes were defined by 47 polymorphic sites from 26 sequences in the dromedaries (Table 1, Supplementary Figure S1, S2). In the Bactrian camel populations, six haplotypes were shared by at least two populations (Supplementary Table S2), and the most frequent haplotype was H_4 (29%; 45 sequences out of 156), followed by haplotype H_3 (13%; 21 sequences out of 156). Several haplotypes included H_9, which was present only in samples from Russia (n = 5), and H_18, which was present only in samples from China (n = 2). Two specific haplotypes (W_1 and W_2) were determined in the extant wild Bactrian camel populations. In the dromedaries, the commonest haplotypes were D_1 and D_9, representing 27% and 15% of all sequences, respectively (Supplementary Table S2; Figure S2).
At the country level, excepting the wild Bactrian camel populations, all other populations showed a high level of maternal genetic diversity, including dromedaries. In the domestic Bactrian camel populations, the highest level of Hd (0.876±0.033) was observed in Mongolia, while the lowest (0.600±0.0215) was found in Kazakhstan; the lowest level of nucleotide diversity (0.00188) was observed in China, while the highest (0.066) was observed in Iran’s domestic Bactrian camel population.
Furthermore, we divided domestic Bactrian camel populations into two groups to investigate genetic diversity: East Asia (China and Mongolia’s domestic Bactrian camel populations) and Central Asia (Iran, Kazakhstan and Russia’s domestic Bactrian camel populations). The number of haplotypes from east Asia (H = 21) was higher than central Asia (H = 13); however, the genetic diversity for populations from Central Asia (Hd = 0.896±0.033 and π = 0.01043) was slightly higher than East Asia (Hd = 0.801±0.030 and π = 0.00238). Interestingly, the wild Bactrian camel population showed the lowest Hd (0.462±0.065) and nucleotide diversity (0.00115), while dromedaries revealed the highest Hd (0.908±0.041) in this study (Table 1).
Population phylogenetic analysis
In order to assess genetic relationships between Bactrian camels and dromedaries, the NJ tree was constructed for the 809 bp CR (Figure 1) and haplotype sequences (Supplementary Figure S3). The NJ tree revealed three well-resolved clusters, namely, domestic Bactrian camels, wild Bactrian camels and dromedaries (Figure 1). Interestingly, two individuals from Kazakhstan (H_2, MH109985) and Russia (H_8, MH109974) were clustered with the dromedaries, which indicated that introgression of mtDNA CR existed between Bactrian camels and dromedaries.
Furthermore, the phylogenetic relationships among domestic Bactrian camel population haplotypes, with H_2 and H_8 excluded, are shown in a MJ network (Figure 2). Haplotypes of each population from different geographic regions (China, Mongolia, Iran, and Kazakhstan) did not cluster together according to their geographic regions. Several common haplotypes, such as H_1, H_3, H_4 and H_20, were shared by individuals from different geographic regions. These results indicated that there was no correspondence between the geographic regions of origin and relationships among populations. All haplotypes networks are shown in Supplementary Figure S4. This result clearly showed that all camel haplotypes were divided into three distinct lineages, similar to the findings from the NJ tree. Meanwhile, for haplotypes H_2 (MH109985) and H_8 (MH109974) from the Bactrian camels, rather than clustering with the Bactrian camel groups, they were closer to the Iranian dromedary. This indicated that there was introgression between Bactrian camels from Russia and Kazakhstan and Iranian dromedary.
Population history and demographic dynamics
We predicted mismatch distribution patterns for domestic Bactrian camel populations from different geographic regions, the wild Bactrian camel population and the dromedary population. To clarify the demographic dynamics of Old World camels, the mismatch distribution graphs revealed uni-, bi- and multi-modal patterns (Supplementary Figure S5), which suggest that there has been a demographic expansion, bottleneck or purifying selection in different camel populations. The above results were compatible with the calculation of population demographic parameters from Table 2. Negative non-significant Tajima’s D values were observed for almost all domestic Bactrian camel populations, with the exception of the Russian Bactrian camel, for which Tajima’s D was negative but significant (p<0.01). The wild Bactrian camel population and dromedaries showed positive and non-significant Tajima’s D values. The test of Tajima’s D is used to distinguish between a DNA sequence evolving randomly (neutrally) and one evolving under a non-random process, including directional selection or balancing selection, demographic expansion or contraction, genetic hitchhiking, or introgression [26]. The Central Asian domestic Bactrian camel populations, the wild Bactrian camel population and dromedaries exhibited positive and non-significant Fu’s Fs values. These results suggest a complex demographic history for the different camel populations. While no significant differences were found for the Harpending raggedness index (r) and SSD in Old World camel populations, the results demonstrated that our total dataset showed relatively good fit to the population expansion model [30].
DISCUSSION
In this study, we assessed the mtDNA genetic diversity of 182 Old World camels including domestic Bactrian camels from China, Mongolia, Russia, Iran and Kazakhstan, wild Bactrian camel from Mongolia; and dromedaries. The aim of our research would reveal Old World camel populations history and demographic dynamics using the most polymorphic CR in mtDNA. This highly polymorphic CR has been used in previous research to assess the genetic diversity and phylogeny of camels, especially in relation to the maternal origin of Old World Camels [15,17,19], including the Bactrian camels.
We observed higher haplotype and nucleotide diversity in the dromedary CRs compared to the domestic and wild Bactrian camel populations (Table 1). This corresponds to previous estimates on whole mitochondrial genome sequencing, and the haplotype and nucleotide diversity of domestic (0.952±0.096) and wild Bactrian camels (0.822±0.097) are lower than dromedary (1.0±0.045) [12]. The wider geographical distribution of the domestic Bactrian camel, which come from China, Mongolia, Kazakhstan, Russia and Iran, and the dromedaries, which come from Iran, Saudi Arabia, United Arab Emirates, Kenya, Sudan, Pakistan and Morocco, could in principle explain these differences. In contrast, the wild Bactrian camel population originated from only one region (the Strictly Protected Area “A” in the Mongolian Gobi Desert), which may restrict its genetic diversity. In addition, researchers believe that the small size of the wild Bactrian camel population, with a census indicating as few as 2,000 individuals remaining [6], is also an important reason for their low genetic diversity [12,15,31]. Similar to previous results [19,15,31], only two haplotypes were discovered (W_1 and W_2) in the wild Bactrian camel population, which again showed its low genetic diversity. A possible caveat of our study is that CR (809 bp) sequence variation alone might not be sufficient to recover population genetic divergence in the wild Bactrian camel population.
In the domestic Bactrian camel populations, there was no distinguishing geographic structuring in the NJ tree and median-joining network among populations from China, Mongolia, Russia, Kazakhstan, and Iran. This result indicated the populations’ single maternal lineage and that strong gene introgression existed in different domestic Bactrian camel populations, especially individuals MH109985 (H_2) and MH109974 (H_8) from Kazakhstan and Russia, which were clustered with the dromedaries; this finding has also been confirmed in whole genome sequencing research [32]. The Silk Road could, in principle, account for this introgression circumstance.
During the process of domestication, population growth and dispersion of animals across a wider geographic range can be determined from molecular signals of sudden expansion [8]. In our study, we obtained positive values of Tajima’s D and Fu’s Fs test in the wild Bactrian camel population, which demonstrates a decrease in population size and/or balancing selection. At the same time, positive values of Fu’s Fs test were found in the Central Asian (i.e., from Russia, Kazakhstan, and Iran) domestic Bactrian camels, which may also be related to a decreasing population size, especially considering that an Iranian population census indicated there were just several hundred individuals remaining. These results are in agreement with reduced genetic diversity and small effective population size in wild and domestic Iranian Bactrian camels. Meanwhile, negative values of Tajima’s D and Fu’s Fs test in East Asian populations (i.e., from China and Mongolia) of domestic Bactrian camels can be explained by demographic expansion and/or positive selection.
CONCLUSION
In summary, we report novel information regarding genetic diversity, population structure and demographic dynamics of Old World camels. The findings obtained from the present study reveal that abundant genetic diversity occurs in domestic Bactrian camel populations and dromedaries, while low levels of haplotype and nucleotide diversity were found in the wild Bactrian camel population.
PDF Links
PubReader
ePub Link
Full text via DOI
Download Citation
Supplement
Print
