Genetic diversity and admixture analysis of Sanfratellano and three other Italian horse breeds assessed by microsatellite markers

Structure and genetic variability of the Costa Rican Paso horse

The Costa Rican Paso Horse (CPC) is a breed developed in Costa Rica. The objectives were to estimate the genetic structure and evaluate the levels of genetic variability of the population. The genotypes of 14 microsatellites in 3654 records (2052 females and 1602 males) were analyzed. Expected (He) and observed (Ho) heterozygosity, polymorphic information content (CIP), fixation index (FIS), Shannon index, as well as Hardy-Weinberg disequilibrium (DHW) were evaluated. Kinship relationships (Rij) were estimated throughout the entire population. The effective population size (Ne) was calculated, alternating allele frequencies less than 0.05, 0.02 and 0.01. The Bayesian clustering study was carried out to infer how many lines are appropriate from the analysis of genotypes using multiple loci. The number of alleles per locus ranged from 7 to 17, with an average value of 9.6; nine loci presented DHW (P < 0.05); two loci presented negative FIS values, the same as Ho > He; the average of CIP, Ho and He was 0.254, 0.756 and 0.785, respectively. At the 12 loci where He > Ho, the differences ranged from 0.002 to 0.341 (0.036 on average). For Ne, the estimates were 201.9, 230.1, and 241.5. In the Rij, 54.86% of the estimates were in the interval of 0.01 to 77.7%. The number of lines that define the population corresponds to three, with an approximate composition of 33.1%, 32.4% and 34.5%, respectively. The CPC, as a subdivided population with DHW and a reduction in heterozygotes may be associated with possible Wahlund effects. Keywords: Wahlund effect, equines, genetic markers, synthetic breed, Hardy Weinberg.

Unravelling of the genetic diversity, population structure, and phylogenetic relationships of the indigenous horse breeds of Pakistan

The current study aimed to evaluate the genetic diversity, level of admixture, and phylogenetic relationship, of the Pakistani horse breeds, along with their morphological characterization. Data for the body measurements showed that Morna horses had the highest values of body height, body length, chest girth, leg length, and head length, whereas the Baluchi horses had the lowest values for these traits. For the genetic diversity 64 animals, 15/breed except for Baluchi(14) and Topras(05), were genotyped by using the 17-plex equine genotyping kit. The AMOVA results showed that 13% of genetic diversity was explained by breed differences, whereas 27% and 60% came from among and within individuals respectively. The highest values of genetic diversity parameters including Na(7.29±0.29), Ne(5.73±0.28), Ho(0.74±0.05) and He(0.82±0.01) were observed for Morna, whereas their lowest values were found for Topras. However, the highest value of inbreeding coefficient (Fis) was found for Baluchi and the lowest for KB horses. Among the markers, CA425 was found as the most polymorphic and ASB23 as the least polymorphic and highly fixed marker. Results of structure analysis revealed that, except Topras, each local horse breed had very admixed genetic structure perhaps due to their continuous crossing with other breeds in the past. Moreover, the structure analysis also showed that Morna and Shiean breeds had mixing of each other which was also confirmed by the lowest value of their pairwise Fst values and likewise the phylogenetic analysis also showed their close genetic relationship with each other. The phylogenetic analysis also revealed that Shiean breed had close genetic relationship with Arabian horses. Collectively, our data showed that Morna is the largest and genetically most diverse horse breed of Pakistan, whereas Baluchi horses are the smallest in size and have the highest values of inbreeding coefficient. And the phylogeny analysis showed that Shiean breed had close genetic relationship with the Arabian.

Exploring the Italian equine gene pool via high-throughput genotyping

Introduction: The Italian peninsula is in the center of the Mediterranean area, and historically it has been a hub for numerous human populations, cultures, and also animal species that enriched the hosted biodiversity. Horses are no exception to this phenomenon, with the peculiarity that the gene pool has been impacted by warfare and subsequent "colonization". In this study, using a comprehensive dataset for almost the entire Italian equine population, in addition to the most influential cosmopolitan breeds, we describe the current status of the modern Italian gene pool. Materials and Methods: The Italian dataset comprised 1,308 individuals and 22 breeds genotyped at a 70 k density that was merged with publicly available data to facilitate comparison with the global equine diversity. After quality control and supervised subsampling to ensure consistency among breeds, the merged dataset with the global equine diversity contained data for 1,333 individuals from 54 populations. Multidimensional scaling, admixture, gene flow, and effective population size were analyzed. Results and Discussion: The results show that some of the native Italian breeds preserve distinct gene pools, potentially because of adaptation to the different geographical contexts of the peninsula. Nevertheless, the comparison with international breeds highlights the presence of strong gene flow from renowned breeds into several Italian breeds, probably due to historical introgression. Coldblood breeds with stronger genetic identity were indeed well differentiated from warmblood breeds, which are highly admixed. Other breeds showed further peculiarities due to their breeding history. Finally, we observed some breeds that exist more on cultural, traditional, and geographical point of view than due to actual genetic distinctiveness.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome‐wide single‐nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model‐based clustering and Neighbour‐Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (H_e) and contemporary effective population size (cNe) were reported in Siciliano (H_e = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (H_e = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome‐wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

Genetic diversity and population structure of Mongolian regional horses with 14 Microsatellite markers

Objective

This study aimed to identify the genetic diversity and population structure of Mongolian horse populations according to the province of residence (Khentii, KTP; Uvs, USP; Omnogovi and Dundgovi, GOP; Khovsgol, KGP) using 14 microsatellite (MS) markers.

Methods

A total of 269 whole blood samples were obtained from the four populations (KTP, USP, GOP, KGP) geographically distinct provinces. Multiplex polymerase chain reaction (PCR) was conducted using 14 MS markers (AHT4, ASB2, ASB17, ASB23, CA425, HMS1, HMS2, HMS3, HMS6, HMS7, HTG4, HTG6, HTG7, and VHL20), as recommended by the International Society for Animal Genetics. Capillary electrophoresis was conducted using the amplified PCR products, alleles were determined. Alleles were used for statistical analysis of genetic variability, Nei’s DA genetic distance, principal coordinate analysis (PCoA), factorial corresponding analysis (FCA), and population structure.

Results

On average, the number of alleles, expected heterozygosity (H_Exp), observed heterozygosity (H_Obs), and polymorphic information content among all populations were 11.43, 0.772, 0.757, and 0.737, respectively. In the PCoA and FCA, GOP, and KGP were genetically distinct from other populations, and the KTP and USP showed a close relationship. The two clusters identified using Nei’s DA genetic distance analysis and population structure highlighted the presence of structurally clear genetic separation.

Conclusion

Overall, the results of this study suggest that genetic diversity between KTP and USP was low, and that between GOP and KGP was high. It is thought that these results will help in the effective preservation and improvement of Mongolian horses through genetic diversity analysis and phylogenetic relationships.

Genetic Variability Trend of Lusitano Horse Breed Reared in Italy

The Lusitano Horse (LH) originates from Portugal, but is reared worldwide. Since 1994, the University of Milan has routinely tested the LHs bred in Italy for parentage control. This study aims to assess the genetic variability of the LH reared in Italy using 16 microsatellites markers. Moreover, the genetic variability changes over the years in the total population (n.384) and in unrelated horses (n.47) were evaluated. Horses were grouped according to their date of birth (1975–1990, 1991–2000, 2001–2010, 2010–2019). Standard genetic diversity parameters, including observed (Ho) and expected (He) heterozygosity, Hardy-Weinberg equilibrium (HWE; P-Val), allelic richness, and inbreeding coefficient (Fis) were estimated. In the whole period, the total population showed Ho as high as 0.69, low Fis (0.057), and imbalance for HWE. When considering the unrelated horses, Ho was seen to increase over time (from 0.594 in 1975–1990 to 0.68 in 2010–2019) and frequencies were in HWE, again having low and decreasing values of Fis (from 0.208 in 1975–1990 to 0.019 in 2010–2019). Bottleneck analysis excluded a recent population decline. Principal Coordinate Analysis at the individual level defined two clusters, the major cluster including all the most recent horses. An increasing number of dams (156% more from 2001–2010 to 2011–2019) supports the good variability recorded in the population so far. However, the high number of foals (77.2%) sired by only four stallions in recent years suggests caution in the choice of the sires for the future.

Assessment of genetic diversity and population structure of Colombian Creole cattle using microsatellites

To establish the diversity, structure, and phylogenetic relationships among Colombian Creole cattle, six native breeds and one introduced breed were genotyped for 20 microsatellite loci. The average number of alleles per breed ranged from 7050 (Romosinuano) to 10,100 (Casanareño), and the expected heterozygosity ranged from 0.691 (San martinero) to 0.785 (Casanareño). The deviation from the Hardy-Weinberg equilibrium (HWE) was statistically significant (p < 0.05) in 59 out of 120 tests carried out in the six breeds for the 20 microsatellite loci analyzed. Colombian Creole bovine breeds have maintained a high level of genetic differentiation within the same populations (93%), and the rest is explained by differences between breeds (7%). The differentiation pattern and the genetic relationships between the Colombian Creole bovine breeds showed high consistency with the evolutionary history of each. Both the Bayesian grouping analysis and the neighbor-joining tree exhibited a reliable grouping pattern, which revealed two main groups: one comprised by the breeds Blanco Orejinegro, Hartón del Valle, Costeño Con Cuernos, Romosinuano, and San Martinero, and the other one by the Creole breed Casanareño and Zebu. These were probably caused by different historical, reproductive, and geographic isolation precedents, as well as by different levels of inbreeding. This study will help understand the genetic characteristics of Colombian Creole cattle and will benefit future conservation programs.

Diversity and effective population size of four horse breeds from microsatellite DNA markers in South-Central Mexico

Abstract. The South-Central region of Mexico has experienced a sizeable introduction of purebred horses for recreational aims. A study was designed to assess effective population sizes and genetic diversity and to verify the genetic integrity of four horse breeds. Using a 12-microsatellite panel, Quarter Horse, Azteca, Thoroughbred and Creole (CRL) horses were sampled and analysed for diversity and genetic structure. Genetic diversity parameters showed high numbers of heterozygous horses but small effective population sizes in all breeds. Population structure results suggested some degree of admixture of CRL with the other reference breeds. The highly informative microsatellite panel allowed the verification of diversity in introduced horse populations and the confirmation of small effective population sizes, which suggests a risk for future breed integrity.

Genetic diversity in the Maremmano horse and its relationship with other European horse breeds

The Maremmano is an Italian warmblood horse breed from central Italy. We characterized the genetic diversity and the degree of admixture in Maremmano in comparison to 14 other European horse breeds using 30 microsatellites. Between-breed diversity explained about 9 per cent of the total genetic diversity. Cluster analysis, genetic distances and genetic differentiation coefficients showed a close relationship of Maremmano with Hanoverian and Lusitano in accordance with breed history.

Genetic diversity and variability in endangered Pantesco and two other Sicilian donkey breeds assessed by microsatellite markers

The genetic variability of Pantesco and other two Sicilian autochthonous donkey breeds (Ragusano and Grigio Siciliano) was assessed using a set of 14 microsatellites. The main goals were to describe the current differentiation among the breeds and to provide genetic information useful to safeguard the Pantesco breed as well as to manage Ragusano and Grigio Siciliano. In the whole sample, that included 108 donkeys representative of the three populations, a total of 85 alleles were detected. The mean number of alleles was lower in Pantesco (3.7), than in Grigio Siciliano and Ragusano (4.4 and 5.9, resp.). The three breeds showed a quite low level of gene diversity (He) ranging from 0.471 in Pantesco to 0.589 in Grigio. The overall genetic differentiation index (Fst) was quite high; more than 10% of the diversity was found among breeds. Reynolds' (D_R) genetic distances, correspondence, and population structure analysis reproduced the same picture, revealing that, (a) Pantesco breed is the most differentiated in the context of the Sicilian indigenous breeds, (b) within Ragusano breed, two well-defined subgroups were observed. This information is worth of further investigation in order to provide suitable data for conservation strategies.

Genetic structure and differentiation of the Italian catria horse

Catria is 1 of the 22 native Italian horse breeds that now survive from a larger number. Thirty individuals, representative of the Catria horse, were analyzed for 11 microsatellites and compared with data of 10 breeds reared in Italy. Three different approaches, genetic distances, correspondence analysis, and clustering methods, were considered to study genetic relationships among Catria and the other horse populations. Genetic differentiation among breeds was highly significant (P < 0.01) for all loci. Average F(ST) values indicate that around 10% of the total genetic variation was explained by the between-breed differences and the 3 approaches utilized gave similar results. Italian native breeds are clearly separated from the other examined breeds. However, by the correspondence analysis, the Catria appears closer to Maremmano and Murgese. The results of Bayesian approaches give further information showing for Catria a common origin with Maremmano and Italian Heavy Draught. Genetic relationships among Catria and the other breeds are consistent with the breed's documented history. The data and information found here can be utilized in the organization of conservation programmes planned to reduce inbreeding and to minimize loss of genetic variability.

Genetic analysis of Sicilian autochthonous horse breeds using nuclear and mitochondrial DNA markers

Genetic diversity and relationship among 3 Sicilian horse breeds were investigated using 16 microsatellite markers and a 397-bp length mitochondrial D-loop sequence. The analysis of autosomal DNA was performed on 191 horses (80 Siciliano [SIC], 61 Sanfratellano [SAN], and 50 Sicilian Oriental Purebred [SOP]). SIC and SAN breeds were notably higher in genetic variability than the SOP. Genetic distances and cluster analysis showed a close relationship between SIC and SAN breeds, as expected according to the breeds' history. Sequencing of hypervariable mitochondrial DNA region was performed on a subset of 60 mares (20 for each breed). Overall, 20 haplotypes with 31 polymorphic sites were identified: A higher haplotype diversity was detected in SIC and SAN breeds, with 13 and 11 haplotypes respectively, whereas only one haplotype was found in SOP. These were compared with 118 sequences from GenBank. BLAST showed that 17 of the 20 haplotypes had been reported previously in other breeds. One haplotype, found in SIC, traces back to a Bronze Age archaeological site (Inner Mongolia). The 3 Sicilian breeds are now rare, and 2 of them are officially endangered. Our results represent a valuable tool for management strategies as well as for conservation purposes.

Genetic variation and population structure of Italian native sheep breeds undergoing in situ conservation

The genetic variability and presence of population substructures in 4 native Northern Italian sheep breeds, Alpagota, Brogna, Foza, and Lamon, undergoing in situ conservation, and 1 widespread Italian breed, Bergamasca, were studied by investigating 19 microsatellite markers. The breeds showed considerable genetic variability in terms of number of alleles and heterozygosity, with the exception of Alpagota, which was the least variable (0.607). Nevertheless, a significant deficit of heterozygotes was observed in each breed due to rather increased levels of inbreeding or to the presence of population substructures, probably caused by increased genetic variation in the founder populations. The analyses evidenced clear genetic differentiation (F(ST) = 0.085), reduced levels of admixture, and presence of private alleles among the breeds, confirming their genetic uniqueness. In particular, according to Reynolds genetic distances, Alpagota was the most differentiated, perhaps because it had been bred mostly in a rather isolated area. Loss of any of the investigated breeds would result in a loss of genetic diversity ranging from 0.5 to 1.6% of the total observed gene diversity. Results supported the decision to safeguard these breeds as important reservoirs of genetic diversity and suggested breeding and mating practices to maintain variability and to overcome within-breed substructures.

Multiscale analysis reveals restricted gene flow and a linear gradient in heterozygosity for an island population of feral horses

We studied the genetic (microsatellite) diversity of a feral population of horses ( Equus caballus L., 1758) on Sable Island, Nova Scotia, Canada (1983–2003), at two spatial scales: (1) for the island as a whole and (2) at the level of four equally sized subdivisions along the length of Sable Island, which is a long (42 km) and narrow (1.5 km) vegetated sand bar. At the island scale (n = 264 horses), observed heterozygosity over 10 loci was 0.647 ± 0.035 (mean ± 1 SE), while expected heterozygosity was 0.696 ± 0.029; we observed significant heterozygote deficiency with all loci considered (P < 0.0001). At the subdivision scale, observed heterozygosity ranged from 0.589 to 0.694 in a gradient from west to east. We observed a corresponding gradient in effective number of alleles and allelic richness. Pairwise values of FST were significant for most subdivision pairs, ranging as high as 0.067 from west to east. Western areas showed highest levels of inbreeding (FIS = 0.113) with outbreeding indicated in the east (FIS = –0.008). Our results suggest that for a large mammal that lives in polygynous social groups, like the feral horse, gene flow along linear habitats (corridors) may be restricted (relative to the dispersal capabilities of the species), even over short distances.