Genome-wide diversity, population structure and demographic history of dromedaries in the central desert of Iran

The development of camel husbandry for good production in a desert climate is very important, thus we need to understand the genetic basis of camels and give attention to genomic analysis. We assessed genome-wide diversity, linkage disequilibrium (LD), effective population size (Ne) and relatedness in 96 dromedaries originating from five different regions of the central desert of Iran using genotyping-by-sequencing (GBS). A total of 14,522 Single Nucleotide Polymorphisms (SNPs) with an average minor allele frequency (MAF) of 0.19 passed quality control and filtering steps. The average observed heterozygosity in the population was estimated at 0.25 ± 0.03. The mean of LD at distances shorter than 40 kb was low (r² = 0.089 ± 0.234). The camels sampled from the central desert of Iran exhibited higher relatedness than Sudanese and lower than Arabian Peninsula dromedaries. Recent Ne of Iran’s camels was estimated to be 89. Predicted Tajima’s D (1.28) suggested a bottleneck or balancing selection in dromedary camels in the central desert of Iran. A general decrease in effective and census population size poses a threat for Iran’s dromedaries. This report is the first SNP calling report on nearly the chromosome level and a first step towards understanding genomic diversity, population structure and demography in Iranian dromedaries.