- Research article
- Open Access
Local and global patterns of admixture and population structure in Iranian native cattle
- Karim Karimi†1Email author,
- Eva M. Strucken†2,
- Nasir Moghaddar2,
- Mohammad H. Ferdosi3,
- Ali Esmailizadeh1, 4 and
- Cedric Gondro2
© The Author(s). 2016
- Received: 13 May 2016
- Accepted: 8 July 2016
- Published: 15 July 2016
Two separate domestication events gave rise to humped zebu cattle in India and humpless taurine cattle in the Fertile Crescent of the Near and Middle East. Iran covers the Eastern side of the Fertile Crescent and exhibits a variety of native cattle breeds, however, only little is known about the admixture patterns of Iranian cattle and their contribution to the formation of modern cattle breeds.
Genome-wide data (700 k chip) of eight Iranian cattle breeds (Sarabi N = 19, Kurdi N = 7, Taleshi N = 7, Mazandarani N = 10, Najdi N = 7, Pars N = 7, Kermani N = 9, and Sistani N = 9) were collected from across Iran. For a local assessment, taurine (Holstein and Jersey) and indicine (Brahman) outgroup samples were used. For the global perspective, 134 world-wide cattle breeds were included. Between breed variation amongst Iranian cattle explained 60 % (p < 0.001) of the total molecular variation and 82.88 % (p < 0.001) when outgroups were included. Several migration edges were observed within the Iranian cattle breeds. The highest indicine proportion was found in Sistani. All Iranian breeds with higher indicine ancestry were more admixed with a complex migration pattern. Nineteen founder populations most accurately explained the admixture of 44 selected representative cattle breeds (standard error 0.4617). Low levels of African ancestry were identified in Iranian cattle breeds (on average 7.5 %); however, the signal did not persist through all analyses. Admixture and migration analyses revealed minimal introgression from Iranian cattle into other taurine cattle (Holstein, Hanwoo, Anatolian breeds).
The eight Iranian cattle breeds feature a discrete genetic composition which should be considered in conservation programs aimed at preserving unique species and genetic diversity. Despite a complex admixture pattern among Iranian cattle breeds, there was no strong introgression from other world-wide cattle breeds into Iranian cattle and vice versa. Considering Iran’s central location of cattle domestication, Iranian cattle might represent a local domestication event that remained contained and did not contribute to the formation of modern breeds, or genetics of the ancestral population that gave rise to modern cattle is too diluted to be linked directly to any current cattle breeds.
- Bos taurus
- Bos indicus
- Fertile crescent
The general consensus about the origin of domesticated cattle is that two separate domestication events took place and gave rise to the variety of cattle breeds we see nowadays . India is the origin of humped zebu cattle (Bos indicus) , and the Fertile Crescent of the Near East is the region of origin of humpless taurine cattle (Bos taurus) . Iran covers the Eastern side of the Fertile Crescent bordering to the West with Turkey and Saudi Arabia, which link to Europe and Africa, and to the East to Afghanistan and Pakistan, which links to India. First agricultural remnants date back 10,000 years when also first cattle domestication is believed to have started [1, 4]. Iran is home to a variety of cattle breeds, however, only little is known about the genetic diversity of Iranian native cattle.
Globalization of breeding programs has become more important and maintaining of local genetic resources is required to facilitate rapid adaptation to changing environment. Indigenous breeds have developed unique characteristics as a response to environmental pressures such as disease and parasite tolerance, heat tolerance, and adaptation to local feed resources [5, 6]. The loss of these breeds or their genetic diversity, which is the ultimate source for the ability to adapt to a changing environment, will significantly limit future breeding programs . Some of the Iranian breeds, such as Golpaigani, have already become extinct while other indigenous breeds have been shown to be on the brink of losing genetic diversity due to small effective population sizes and inbreeding [8, 9]. Middle Eastern cattle breeds represent the main links to the ancient history of taurine domestication and may also be relevant as a future source of currently untapped genetic material. Characterization of the genetic variability and breed composition of these populations will assist in guiding preservation programs and may provide additional insights into the domestication process of taurine cattle.
The availability of high density genome-wide SNP arrays has given researchers a powerful tool to characterize genetic diversity and breed composition [10–13]. Genome-wide information provides a fine-grained raster, compared to for example microsatellites, to trace even small differences between animal populations. Thus, the history of migration and mating events should be possible to be reconstructed more accurately .
We used high-density SNP data to investigate genetic diversity, admixture and population structures in eight Iranian native cattle breeds. Further, we incorporated our data set with information from 134 world-wide cattle breeds  based on 18,892 common SNPs in order to achieve an assessment of the genetic history and structure of Iranian cattle populations. These results are the first comprehensive evaluation of the valuable resource on native Iranian cattle diversity in a historically important geographic region for cattle domestication.
Genetic structure within Iranian cattle breeds
Analysis of molecular variance in 8 Iranian cattle populations and 3 outgroup cattle breeds
Source of variance
Variance components (%)
Iran + Outgroups
Population genetic estimates of 8 Iranian and 3 outgroup cattle breeds based on autosomal chromosomes
Ne [ 9 ]
% SNPs MAF ≤ 0.01
Pairwise FST values among 8 Iranian cattle breeds and 3 outgroup cattle breeds
The semiarid and arid South of the country is inhabited by the remaining four breeds. Sistani, Pars, and Najdi are linked via the Kermani breed of the central Kerman region, which is genetically the most closely related to all three breeds (Table 3). The Kermani population has some ancestry that can be traced to a population that separated from the Sarabi into Kurdi and Taleshi. The Najdi and Pars populations showed some influx from a population that was already on the verge to separating into the branch that developed into Taleshi, Mazandarani, and the southern breeds (Fig. 1).
Adding a third potential founder population separated the Holstein, Jersey, and Brahman (Fig. 3b), and revealed that the Iranian breeds are an admixture of taurine breeds, mainly corresponding to Holstein, and indicine breeds as represented by the Brahmans. To a lesser degree, a Jersey component was present in the Iranian breeds (Fig. 3b). It is often assumed that both Holstein and Jersey were used for crossbreeding with native cattle in Iran . However, the low influence of Jersey in all Iranian breeds suggest that they were not widely used for crossbreeding purposes (e.g., with Najdi cattle) and the breeding history should be reconsidered.
Assuming four ancestral populations provided the lowest cross-validation standard error (0.557) indicating that this is the most likely number of ancestral breeds based on the data of this study. With four founder breeds, the Sarabi separated into its own distinct breed (an Iranian taurine). Based on word-of-mouth from local farmers, Sarabi had no introgression from other breeds for the last 50 years and represents the fourth ancestral breed of the Iranian cattle populations in our study (Fig. 3b). With a fifth ancestor population, Taleshi and Mazandarani separate from the other breeds which could be associated with their relatively isolated location by the Caspian Sea. The most admixed breed is the Kurdi, made up of a high proportion of Iranian taurine followed by a Holstein and Jersey background (Fig. 3b).
From this Iranian focused perspective, we can summarize that there is a strong West to East distribution of taurine to indicine breed proportions. Predicting a breed (taurine or indicine) based on phenotypic appearance is highly misleading and should always backed-up by genetic analyses. There has been some migration within the Iranian breeds. The Sarabi appear to be most homogenous representing an independent taurine population whilst the Sistani represent a relatively homogenous indicine population. The Kurdi were the most heterogeneous population. To anchor the Iranian breeds in a global perspective and to infer breed proportions and migration events from populations outside Iran, we included a further 134 breeds in the following section.
Admixture patterns between Iranian cattle and other world-wide cattle breeds
We used the ADMIXTURE program in an unsupervised analysis with all 142 breeds as well as with a reduced data set of 44 breeds (highlighted breeds in Additional file 2: Table S1, note that Brown Swiss BSW were treated as two populations BRU and BSW for the analysis as per Decker et al. 2014 ) to infer general patterns of admixture and genetic structure. We will describe here only the analysis with 44 breeds; the detailed output with 2–5 assumed founder populations can be found in Additional files 3: Figure S2, 4: Figure S3, 5: Figure S4, and 6: Figure S5. We further calculated f 3 statistics for all possible triplet groups of the 44 breeds, and f 4 statistics with the Iranian cattle breeds as the sister, and African taurine (Somba, Lagune, Baole, N’Dama and Oulmès Zaer), Anatolian taurine (Turkish Grey, Anatolian Black, East Anatolian Red, Anatolian Southern Yellow, and South Anatolian Red), Asian indicine (Hissar, Gabrali, Dajal, Bhagnari, Rojhan, Sahiwal, Gir, Cholistani, Tharparkar, Red Sindhi, and Achai), Asian taurine (Hanwoo, Wagyu, and Mongolian), and European taurine (Angus, Hereford, Lincoln Red, Holstein, Jersey, Guernsey, and Brown Swiss) as opposing sister groups. The f 3 and f 4 statistics are similar to the fixation index (FST statistic), with the exception that the genetic relationship between three or four populations are considered simultaneously rather than just between two populations. Both f 3 and f 4 tests were designed to detect admixture in a population based on the other populations submitted to the test.
With 19 assumed ancestral populations (lowest cross-validation error), the Sarabi formed a distinct taurine genetic cluster of which larger ancestries can be found in Kurdi, Taleshi, Mazandarani, and Pars. The Sistani as well as the Kermani formed a distinct indicine population with minor traces of Indian Gir cattle (2.1 and 3.2 %, respectively). Proportions of Sistani/Kermani were also present in large percentages in the other Iranian cattle breeds apart from Sarabi (Fig. 5). For Iranian cattle breeds, we found only 26 significant f 3 tests that all belonged to the Kermani breed. The Sistani breed was always included as a sister group confirming the strong influence of Sistani on the formation of the Kermani breed. The most extreme test had Sistani and Jersey breeds as the sister groups with a Z-score of −9.95, which opposes our findings from the previous section in regards to the use of Jersey cattle for crossbreeding purposes with the Iranian native breeds. However, it is more likely that the entire admixture signal of the f 3 statistic stems purely from the Sistani breed. Furthermore, f 4 tests with Asian indicine breeds showed 43 significant tests among 420 possible tests. The most significant test included Sistani and Kurdi as sister and Achai and Gir as opposing sister groups (f 4 = 0.0014, Z-score = 6.97). Of 43 significant tests, 35 and 11 contained Achai and Gir, respectively. Sistani had the most significant tests (26) among Iranian cattle. These results of the f 4 statistic are in accordance with the ADMIXTURE findings of an introgression of Gir into Sistani and Kermani (Fig. 5). This introgression might be explained by the closer geographic location of these Iranian breeds to the Indian sub-continent where Gir originally stem from. The Persian Empire (550 BC-651 AD) occupied land from Africa to Eastern Europe and the Indus Valley . Again, whilst migration of cattle within the Persian Empire and across its borders can historically be assumed, further proof or recordings of such a migration is lacking. No significant introgression was found with Hissar, Gabrali, Dajal, Bhagnari, or Rojhan as opposing sister groups.
As in the previous section, the Kurdi were the most heterogeneous breed and some traces of Brown Swiss were detected (max 23.3 %, min 3.4 %). Even though some Kurdi animals also showed traces of Jersey, an active crossing can still be rejected based on these results. Out of 588 f 4 tests with European taurine cattle, 100 tests were significant, and Brown Swiss were included in most of the significant tests (46 tests), followed by Holstein (42 tests) and Jersey (27 tests).
Modern Iranian cattle breeds are located in a region associated with first cattle domestication, therefore allowing for the hypothesis that these breeds should represent some form of ancestor, we could not find substantial proof for this hypothesis. Breed proportions of Iranian populations in other world-wide cattle breeds are limited to some traces of Sarabi in Korean Hanwoo cattle, and Anatolian breeds (Mongolian, Turkish Grey, Anatolian Southern Yellow, South Anatolian Red, east Anatolian Red; Fig. 5). The geographical location of the Anatolian breeds close to Iran and the extend of historic Empires might explain migration events from Iranian breeds across borders. The breed content in the Korean Hanwoo cattle, however, might have been a larger migration such as the Silk Road traffic. Out of 85 possible f 4 tests with Iranian cattle as sister and Wagyu, Hanwoo, and Mongolian cattle as opposing sister populations, only one significant test was found (f 4 (Kurdi, Kermani; Hanwoo, Mongolian) = −0.0009 (Z-score = −3.55, alternative trees had Z-scores of 46.5 and 50.7).
Phylogeny and migration events
Ten migration events based on allele fractions between the populations were sequentially added with the final model explaining 99.32 % of the variance in relatedness between populations (Fig. 6). As in the ADMIXTURE and f 4 analyses, a larger influx of Brown Swiss into the Kurdi population was observed, underpinning that there was no active crossing with Jersey cattle. The Iranian breeds with larger indicine breed proportions (Najdi, Pars, Kermani, and Sistani) were influenced by a common ancestor with indicine Dajal cattle (Punjab, Pakistan). Dajal cattle have more than 50 % of Gir content according to our ADMIXTURE analysis which might explain the Gir proportion that we previously observed in the Sistani and Kermani breeds (Fig. 5).
Several migration events of Iranian cattle breeds were observed. A strong indication is given for a crossing of Sistani with Pars cattle, which was not observed in the previous section (Fig. 6). Further, migration from the Kurdi population into an ancestral population of the European Holstein, as well as from a common ancestor of the Kurdi to an ancestral population of African cattle breeds (Somba, Lagune, and Baoule) was observed and in concordance with our f 4 statistics. The Sarabi breed appears to be excluded from migration events confirming their unique position within the Iranian breeds.
Our results provide novel information about the genetic structure and admixture of present day cattle breeds inhabiting a center of a historical domestication event. We showed that the eight Iranian cattle breeds feature discrete genetic characteristics which have to be considered in conservation programs aiming at preserving unique species and genetic diversity. Despite a complex admixture pattern among Iranian cattle breeds, we did not find strong introgression from other world-wide cattle breeds into Iranian cattle. A clear geographical distribution of taurine influences in the North-West and indicine influence in the South-East was found with Sarabi forming a unique taurine breed and Sistani and Kermani forming unique indicine breeds. Other Iranian breeds are, with minor exceptions, composed of these unique Iranian breeds. Minor introgression from Romagnola, Brown Swiss, and Gir were found, however, in such low quantities that a prolonged or breed changing admixture could be excluded. Further, minor contents of the Iranian Sarabi breed were found in Korean Hanwoo and Anatolian cattle breeds. Despite the geographic location of the Iranian breeds in an ancient centre of cattle domestication, we did not find more evidence of introgression of these breeds into other world-wide cattle breeds. Possible interpretations of these results are: 1. Iranian breeds remained fairly unchanged since their formation and main migration events occurred within the country. Thus, the Iranian cattle breeds might be seen as a living link to ancient domestication events; however, their contribution to the formation of modern cattle breeds seems to be marginal. Whilst there might have been only two major domestication events leading to the differentiation of Bos taurus and Bos indicus breeds, many smaller and local events might have taken place at the same time, with the Iranian breeds representing one of these local domestication events that survived into the present. 2. Historical invasions and migration routes between Europe, Africa, and Asia resulted in a highly admixed cattle population inhabiting the ancient Fertile Crescent. Traces of this ancestral population might be so diluted in modern cattle breeds that it is difficult to pin-point the contribution or genetic link between cattle breeds of Iran and other global breeds.
Hair samples were collected from a total of 90 individuals representing eight different breeds of Iranian native cattle. As far as possible, unrelated animals were selected either based on pedigree recordings or information provided by farmers. To ensure a good representation of breeds across Iran, we considered characteristics of the area from which the cattle were sampled such as ethnic community, production system (pastoral or crop and livestock), and ecological zones (Fig. 3). Samples were collected from the following breeds: Sarabi (N = 20), Kurdi (N = 10), Taleshi (N = 10), Mazandarani (N = 10), Najdi (N = 10), Pars (N = 10), Kermani (N = 10), and Sistani (N = 10). We further included Holstein, Jersey, and Brahman cattle samples sourced from the Bovine HapMap Consortium to be taurine and indicine outgroups to the Iranian samples. Regarding the prominent role of Holstein and Jersey breeds in crossbreeding programs with Iranian cattle, this information was used to check whether animals correctly represented their predefined populations. In a global analysis, 1,557 animals from 134 cattle breeds sourced from Decker et al. (2014)  were used instead of the Bovine HapMap data.
Genotyping and quality control
Genomic DNA of the Iranian samples was extracted from hair roots and SNP genotyping was performed using the Illumina high-density Bovine BeadChip (Illumina, Inc, San Diego, CA, USA) designed to genotype 777,962 SNPs. Quality control of the autosomal genotypes was performed using the program snpQC  across all Iranian breeds and per outgroup breed. Filtering parameters included a GC score >0.9 (418,571 SNPs failed), call rates per marker >90 % (295,607 SNPs failed) and per animal >70 % (14 animals failed). Further, markers that deviated in their heterozygosity by more than 3 standard deviations from the mean heterozygosity (3,243 SNPs failed), and that deviated from Hardy-Weinberg equilibrium at P-value < 10e–16 (7,550 SNPs failed) were excluded. After quality control, 283,028 SNPs that passed the filter criteria in all breeds remained for further analysis. The number of Iranian samples used after quality control was 19 Sarabi, 7 Kurdi, 7 Taleshi, 10 Mazandarani, 7 Najdi, 7 Pars, 9 Kermani, and 9 Sistani animals. From the quality controlled outgroup breeds, 15 animals per breed were selected to create a balanced data set.
Genetic differentiation and population structure within Iranian cattle breeds
Inbreeding coefficients (FIS) and pairwise FST values, which describe the genetic differentiation between two populations, were estimated according to Weir and Cockerham (1984) . Analysis of molecular variances (AMOVA)  were carried out with the StAMPP package in R . The AMOVA analysis provides further information on the genetic differentiation within a population and between populations. Further, a principal components analysis based on the genetic data (PCA) was carried out with the SNPRelate package in R  describing patterns of population differentiation and overlap. Clustering of breeds into genetic groups and estimation of breed proportions of ancestral breeds was carried out with an unsupervised analysis in ADMIXTURE 1.23 . The best number of ancestral populations (K) was inferred via the program’s cross-validation procedure for 1 to 11 assumed populations. The number of ancestral populations with the best predictive accuracy was based on the lowest standard error of the cross-validation error estimate. These analyses were carried out with three outgroup breeds (Holstein, Hersey, and Brahman) to anchor the eight Iranian cattle breeds. The indicine breed proportion was calculated as the proportion of the Brahman breed at K = 2. The TreeMix software  was used to model gene flow between the Iranian cattle populations and create maximum likelihood trees. TreeMix computes a covariance matrix of all populations based on allele frequencies, comparing two populations (Xi and Xj) with respect to a common ancestral population (xA) (Cov (Xi,Xj) = E[(Xi-xA)(Xj-xA)]). Migration events are modelled by allowing ancestry from multiple ancestral populations and weighting the contribution of each ancestral population by its fraction of alleles in the descendent population. Four independent runs with different seeds for five migration edges were carried out and Sarabi were set as the root of the tree due to their clear separation from other breeds in the previously described analyses.
Placement of Iranian cattle breeds in world-wide patterns of admixture
Genotyping data from Iranian cattle samples were merged with the data set from Decker et al. (2014) . The final data set comprised 1,632 individuals from 142 cattle breeds genotyped for 18,892 common autosomal SNPs. As before, the SNPRelate package  was used for a PCA analysis.
To decrease computing time, 584 individuals from 44 cattle breeds were selected from Decker et al. (2014) . Selection was based on purity (based on PCA results), number of samples, and best representation of major genetic groups: European taurine, African taurine, Asian taurine, Anatolian cattle, Iranian cattle, and indicine populations. An unsupervised ADMIXTURE analyses was run for K values from 1 to 20 with a cross-validation procedure for the complete (142 breeds) and reduced dataset (44 breeds). Maximum likelihood phylogeny of the 44 breeds was created in TreeMix  with and without migration events. Balinese cattle were set as the graph root and blocks of 1,000 SNPs were used to account for potential linkage disequilibrium between nearby SNPs. Ten migration edges were sequentially added. Four independent runs with different seeds were carried out to examine the consistency of the migration edges. Furthermore, the THREEPOP program implemented in TreeMix was applied to evaluate the admixture history among populations based on f 3  and f 4  statistics. The f 3 statistic was carried out for all possible triplets among the 142 breeds [f 3 (X; A, B)]. The f 4 statistic was calculated for several subsets of the populations using the FOURPOP program from TreeMix [f 4 (X, Y; A, B)]. Both the f 3 and f 4 tests are designed to detect admixture in a population based on the other populations submitted to the test.
AD, after Christ; AMOVA, analysis of molecular variance; BC, before Christ; DF, degrees of freedom; FIS, inbreeding coefficient; FST, fixation index; GC score, GenCall score, confidence measure for genotype calling; He, heterozygosity; K, number of assumed ancestral population in an ADMIXTURE analysis; MAF, minor allele frequency; MSD, mean square error; N, sample size; NA, data not available; Ne, effective population size; PCA, principal component analysis; SNP, single nucleotide polymorphism; SSD, sums of square deviation
We would like to thank the farmers for their participation and cooperation in the data sampling of the Iranian cattle.
This project was partially funded by a grant from the Next-Generation BioGreen 21 Program (No. PJ01134903) and the Cooperative Research Program for Agriculture Science & Technology Development (PJ00640505), Rural Development Administration, Republic of Korea.
Availability of data and materials
The datasets analyzed during the current study are available in the Dryad repository, provisional DOI: doi:10.5061/dryad.nq189.
KK and EMS carried out the analyses and drafted the manuscript. CG and AEK conceived of the study, and critically revised the manuscript. KK, EMS, CG, AEK, NM and MF were involved in the design of the study and data sampling. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
No ethics approval was required to obtain hair samples from Iranian cattle. Permission for hair sampling was obtained from the farmers on site.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Loftus RT, MacHugh DE, Bradley DG, Sharp PM, Cunningham P. Evidence for two independent domestications of cattle. Proc Natl Acad Sci U S A. 1994;91(7):2757–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen S, Lin BZ, Baig M, Mitra B, Lopes RJ, Santos AM, Magee DA, Azevedo M, Tarroso P, Sasazaki S, et al. Zebu cattle are an exclusive legacy of the South Asia neolithic. Mol Biol Evol. 2010;27(1):1–6.View ArticlePubMedGoogle Scholar
- Bollongino R, Burger J, Powell A, Mashkour M, Vigne J-D, Thomas MG. Modern taurine cattle descended from small number of near-eastern founders. Mol Biol Evol. 2012;29(9):2101–4.View ArticlePubMedGoogle Scholar
- Hole F. Neolithic Age in Iran. In: Encyclopaedia Iranica. onlineth ed. New York: Columbia University; 2004.Google Scholar
- Soma P, Kotze A, Grobler JP, van Wyk JB. South African sheep breeds: population genetic structure and conservation implications. Small Rumin Res. 2012;103(2–3):112–9.View ArticleGoogle Scholar
- Hoffmann I. Climate change and the characterization, breeding and conservation of animal genetic resources. Anim Genet. 2010;41 Suppl 1:32–46.View ArticlePubMedGoogle Scholar
- Herrero-Medrano JM, Megens HJ, Groenen MA, Ramis G, Bosse M, Perez-Enciso M, Crooijmans RP. Conservation genomic analysis of domestic and wild pig populations from the Iberian Peninsula. BMC Genet. 2013;14:106.View ArticlePubMedPubMed CentralGoogle Scholar
- Strucken EM, Karimi K, Esmailizadeh Koshkoiyeh A, Moghaddar N, Gondro C. Genetic Diversity and Effective Population Sizes of Eight Iranian Cattle Breeds. In: Association for the Advancement of Animal Breeding and Genetics (AAABG). Lorne: Association for the Advancement of Animal Breeding and Genetics (AAABG); 2015.Google Scholar
- Karimi K, Esmailizadeh Koshkoiyeh A, Fozi MA, Porto Neto LR, Gondro C. Prioritization for conservation of Iranian native cattle breeds based on genome-wide SNP data. Conserv Genet. 2016;17:77–89.View ArticleGoogle Scholar
- Decker JE, McKay SD, Rolf MM, Kim J, Molina Alcala A, Sonstegard TS, Hanotte O, Gotherstrom A, Seabury CM, Praharani L, et al. Worldwide patterns of ancestry, divergence, and admixture in domesticated cattle. PLoS Genet. 2014;10(3), e1004254.View ArticlePubMedPubMed CentralGoogle Scholar
- Gautier M, Faraut T, Moazami-Goudarzi K, Navratil V, Foglio M, Grohs C, Boland A, Garnier JG, Boichard D, Lathrop GM, et al. Genetic and haplotypic structure in 14 European and African cattle breeds. Genetics. 2007;177(2):1059–70.View ArticlePubMedPubMed CentralGoogle Scholar
- McCue ME, Bannasch DL, Petersen JL, Gurr J, Bailey E, Binns MM, Distl O, Guerin G, Hasegawa T, Hill EW, et al. A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies. PLoS Genet. 2012;8(1), e1002451.View ArticlePubMedPubMed CentralGoogle Scholar
- Vonholdt BM, Pollinger JP, Lohmueller KE, Han E, Parker HG, Quignon P, Degenhardt JD, Boyko AR, Earl DA, Auton A, et al. Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication. Nature. 2010;464(7290):898–902.View ArticlePubMedPubMed CentralGoogle Scholar
- Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, Rogel-Gaillard C, Park C, Milan D, Megens HJ, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491(7424):393–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Strucken EM, Lee SH, Jang GW, Porto-Neto LR, Gondro C. Towards breed formation by island model divergence in Korean cattle. BMC Evol Biol. 2015;15(1):284.View ArticlePubMedPubMed CentralGoogle Scholar
- Amos W, Balmford A. When does conservation genetics matter? Heredity (Edinb). 2001;87(Pt 3):257–65.View ArticleGoogle Scholar
- Lachance J, Tishkoff SA. SNP ascertainment bias in population genetic analyses: why it is important, and how to correct it. Bioessays. 2013;35(9):780–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19(9):1655–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Pickrell JK, Pritchard JK. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012;8(11), e1002967.View ArticlePubMedPubMed CentralGoogle Scholar
- Bovine HapMap C, Gibbs RA, Taylor JF, Van Tassell CP, Barendse W, Eversole KA, Gill CA, Green RD, Hamernik DL, Kappes SM, et al. Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science. 2009;324(5926):528–32.View ArticleGoogle Scholar
- Lin BZ, Sasazaki S, Mannen H. Genetic diversity and structure in Bos taurus and Bos indicus populations analyzed by SNP markers. Anim Sci J. 2010;81(3):281–9.View ArticlePubMedGoogle Scholar
- McKay SD, Schnabel RD, Murdoch BM, Matukumalli LK, Aerts J, Coppieters W, Crews D, Dias Neto E, Gill CA, Gao C, et al. An assessment of population structure in eight breeds of cattle using a whole genome SNP panel. BMC Genet. 2008;9:37.View ArticlePubMedPubMed CentralGoogle Scholar
- O’Brien AM, Holler D, Boison SA, Milanesi M, Bomba L, Utsunomiya YT, Carvalheiro R, Neves HH, da Silva MV, VanTassell CP, et al. Low levels of taurine introgression in the current Brazilian Nelore and Gir indicine cattle populations. Genet Sel Evol. 2015;47:31.View ArticlePubMedPubMed CentralGoogle Scholar
- Domestic Animal Diversity Information Service (DAD-IS) [http://dad.fao.org/] accessed 4 Apr 2016.
- Bradley DG, MacHugh DE, Loftus RT, Sow RS, Hoste CH, Cunningham EP. Zebu-taurine variation in Y chromosomal DNA: a sensitive assay for genetic introgression in west African trypanotolerant cattle populations. Anim Genet. 1994;25(1):7–12.View ArticlePubMedGoogle Scholar
- Mason I. Evolution of Domesticated Animals. London: Longman; 1984.Google Scholar
- Ben Jemaa S, Boussaha M, Ben Mehdi M, Lee JH, Lee S-H. Genome-wide insights into population structure and genetic history of tunisian local cattle using the illumina bovinesnp50 beadchip. BMC Genomics. 2015;16(1):1–12.View ArticleGoogle Scholar
- Roux G. Ancient Iraq. 3rd ed. London, England: Penguin Books; 1992.Google Scholar
- Elisseeff V. The Silk Roads: Highways of Culture and Commerce. New York: Berghahn Books; 2001.Google Scholar
- Curtis VS, Stewart S. The Idea of Iran. 3rd ed. London: I.B. Tauris & Co Ltd; 2010.Google Scholar
- Gondro C, Porto-Neto LR, Lee SH. SNPQC--an R pipeline for quality control of Illumina SNP genotyping array data. Anim Genet. 2014;45(5):758–61.View ArticlePubMedGoogle Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38(6):1358–70.View ArticleGoogle Scholar
- Excoffier L, Smouse PE, Quattro JM. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992;131(2):479–91.PubMedPubMed CentralGoogle Scholar
- Pembleton LW, Cogan NOI, Forster JW. StAMPP: an R package for calculation of genetic differentiation and structure of mixed-ploidy level populations. Mol Ecol Resour. 2013;13(5):946–52.View ArticlePubMedGoogle Scholar
- Zheng X, Levine D, Shen J, Gogarten SM, Laurie C, Weir BS. A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics. 2012;28(24):3326–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19.Google Scholar
- Reich D, Thangaraj K, Patterson N, Price AL, Singh L. Reconstructing Indian population history. Nature. 2009;461(7263):489–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Keinan A, Mullikin JC, Patterson N, Reich D. Measurement of the human allele frequency spectrum demonstrates greater genetic drift in East Asians than in Europeans. Nat Genet. 2007;39(10):1251–5.View ArticlePubMedPubMed CentralGoogle Scholar