High microsatellite and mitochondrial diversity in Anatolian native horse breeds shows Anatolia as a genetic conduit between Europe and Asia

Genetic Diversities and Historical Dynamics of Native Ethiopian Horse Populations (Equus caballus) Inferred from Mitochondrial DNA Polymorphisms

Matrilineal genetic diversity and relationship were investigated among eight morphologically identified native Ethiopian horse populations using polymorphisms in 46 mtDNA D-loop sequences (454 base pairs). The horse populations identified were Abyssinian, Bale, Borana, Horro, Kafa, Kundido feral horses, Ogaden and Selale. Mitochondrial DNA D-loop sequences were characterized by 15 variable sites that defined five different haplotypes. All genetic diversity estimates, including Reynolds' linearized genetic distance, genetic differentiation (F_ST) and nucleotide sequence divergence (D_A), revealed a low genetic differentiation in native Ethiopian horse populations. However, Kundido feral and Borana domestic horses were slightly diverged from the rest of the Ethiopian horse populations. We also tried to shed some light on the matrilineal genetic root of native Ethiopian horses from a network constructed by combining newly generated haplotypes and reference haplotypes deposited in the GenBank for Eurasian type Turkish Anatolian horses that were used as a genetic conduit between Eurasian and African horse populations. Ninety-two haplotypes were generated from the combined Ethio-Eurasian mtDNA D-loop sequences. A network reconstructed from the combined haplotypes using Median-Joining algorithm showed that haplotypes generated from native Ethiopian horses formed separate clusters. The present result encourages further investigation of the genetic origin of native African horses by retrieving additional mtDNA sequences deposited in the GenBank for African and Eurasian type horses.

Major inconsistencies of inferred population genetic structure estimated in a large set of domestic horse breeds using microsatellites

STRUCTURE remains the most applied software aimed at recovering the true, but unknown, population structure from microsatellite or other genetic markers. About 30% of structure-based studies could not be reproduced (Molecular Ecology, 21, 2012, 4925). Here we use a large set of data from 2,323 horses from 93 domestic breeds plus the Przewalski horse, typed at 15 microsatellites, to evaluate how program settings impact the estimation of the optimal number of population clusters K opt that best describe the observed data. Domestic horses are suited as a test case as there is extensive background knowledge on the history of many breeds and extensive phylogenetic analyses. Different methods based on different genetic assumptions and statistical procedures (dapc, flock, PCoA, and structure with different run scenarios) all revealed general, broad-scale breed relationships that largely reflect known breed histories but diverged how they characterized small-scale patterns. structure failed to consistently identify K opt using the most widespread approach, the ΔK method, despite very large numbers of MCMC iterations (3,000,000) and replicates (100). The interpretation of breed structure over increasing numbers of K, without assuming a K opt, was consistent with known breed histories. The over-reliance on K opt should be replaced by a qualitative description of clustering over increasing K, which is scientifically more honest and has the advantage of being much faster and less computer intensive as lower numbers of MCMC iterations and repetitions suffice for stable results. Very large data sets are highly challenging for cluster analyses, especially when populations with complex genetic histories are investigated.

Genetic variability of Akhal-Teke horses bred in Italy

Background

The Akhal-Teke horse (AKH) is native of the modern Turkmenistan area. It was introduced in Italy from 1991 to 2000 mainly as an endurance horse. This paper characterizes the genetic variability of the whole Italian AKH horse population and evaluates their inbreeding level by analyzing microsatellite markers and mitochondrial D-Loop sequences.

Methods

Seventeen microsatellite marker loci were genotyped on 95 DNA samples from almost all the AKH horses bred in Italy in the last 20 years. Standard genetic variability measures (H_o, H_e, F_IS) were compared against the same variables published on other eight AKH populations. In addition, 397 bp of mtDNA D-loop region were sequenced on a sub-group of 22 unrelated AKH out of the 95 sampled ones, and on 11 unrelated Arab horses. The haplotypes identified in the Italian population were aligned to sequences of AKH (56), Arab (five), Caspian Pony (13), Przewalskii (two) and Barb (15) horses available in GenBank. The Median Joining Network (MJN), Principal Component Analysis (PCA) and Neighbor-joining (NJ) tree were calculated on the total 126 sequences.

Results

Nucleic markers showed a high degree of polymorphism (H_o = 0.642; H_e = 0.649) and a low inbreeding level (F_IS = 0.016) in Italian horses, compared to other AKH populations (ranged from −0.103 AKH from Estonia to 0.114 AKH from Czech Republic). High variability was also recorded in the D-Loop region. 11 haplotypes were identified with haplotype diversity (hd), nucleotide diversity (π) and average number of nucleotide differences (k) of 0.938, 0.021 and 6.448, respectively. When all the 126 D-Loop sequences were compared, 51 haplotypes were found, and four were here found only in the Italian AKH horses. The 51 haplotypes were conformed to eight recognized mtDNA haplogroups (A, C, F, G, L, M, P and Q) and confirmed by MJN analysis, Italian horses being assigned to five haplogroups (A, C, G, L and M). Using a PCA approach to the same data, the total haplotypes were grouped into two clusters including A+C+M+P and G+F haplogroups, while L and Q haplogroups remained ungrouped. Finally, the NJ algorithm effectively discretizes only the L haplogroup. All the above data univocally indicate good genetic variability and accurate management of the Akhal-Teke population in Italy.

Genetic analysis of maternal and paternal lineages in Kabardian horses by uniparental molecular markers

Studies of mitochondrial DNA (mtDNA) as well as the non-recombining part of the Y chromosome help to understand the origin and distribution of maternal and paternal lineages. The Kabardian horse from Northern Caucasia which is well-known for strength, stamina and endurance in distance riding has a large gap in its breeding documentation especially in the recent past. A 309 bp fragment of the mitochondrial D-loop (156 Kabardian horses) and six mutations in Y chromosome (49 Kabardian stallions), respectively, were analyzed to get a better insight into breeding history, phylogenetic relationship to related breeds, maternal and paternal diversity and genetic structure. We found a high mitochondrial diversity represented by 64 D-loop haplotypes out of 14 haplogroups. The most frequent haplogroups were G (19.5%), L (12.3%), Q (11.7%), and B (11.0%). Although these four haplogroups are also frequently found in Asian riding horses (e.g. Buryat, Kirghiz, Mongolian, Transbaikalian, Tuvinian) the percentage of the particular haplogroups varies sometimes remarkable. In contrast, the obtained haplogroup pattern from Kabardian horse was more similar to that of breeds reared in the Middle East. No specific haplotype cluster was observed in the phylogenetic tree for Kabardian horses. On Kabardian Y chromosome, two mutations were found leading to three haplotypes with a percentage of 36.7% (haplotype HT1), 38.8% (haplotype HT2) and 24.5% (haplotype HT3), respectively. The high mitochondrial and also remarkable paternal diversity of the Kabardian horse is caused by its long history with a widely spread maternal origin and the introduction of Arabian as well as Thoroughbred influenced stallions for improvement. This high genetic diversity provides a good situation for the ongoing breed development and performance selection as well as avoiding inbreeding.

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

Mitochondrial diversity in mountain horse population from the South-Eastern Europe

It is believed that population structure of mountain horse breeds is preserved from any genetic introgression, because of their geographical isolation and specific semi-wild life style of husbandry. Till date there are no molecular data for the Balkan horses. In this study we try to give information about some autochthonous mountain horse breeds from Bulgaria. A total of 121 horses from three different mountain massifs are presented: Stara Planina (the Balkan mountain), the Rhodopes and Rila-Pirin massif were genotyped according to mitochondrial D-loop region. The results showed huge diversity of all known haplogroups with exception of C, F and R. West Eurasian haplogroups B, D, M and L were with the highest frequencies. Haplogroups A, J, I, O'P and Q were also observed with the highest frequencies, but not equally distributed among the three populations. Analyses of the horse breeds reveal preserved genetic profile of the Balkan and the Rhodopes mountains populations. In contrast, a Rila-Pirin breed unexpectedly showed mixed profile - a massive genetic introgression with an Asiatic-type haplogroups. A similar mixed Euro-Asiatic haplotype profile possessed the Carpathian mountain pony, although both populations are separated geographically and historically. The genetic pool of three Bulgarian mountain horse populations is highly heterogenic and because of that these breeds should be preserved.