Genealogical analyses in open populations: the case of three Arab-derived Spanish horse breeds

Genetic Variability and Conservation Challenges in Lithuanian Dairy Cattle Populations

The purpose of the study was to investigate the genetic variability of open Lithuanian Red and Red-and-White (LRWP) and Lithuanian Black-and-White (LBWP) dairy cattle populations and indicate the differences from the old genotypes of Lithuanian Black-and-White (LBW) and Lithuanian Red cattle (LR), which are currently under a conservation program. In order to gain a better understanding of the populations under conservation and to minimize the potential influence of other breeds, a distinct subgroup was formed that comprised animals whose father and mother belonged to the same breed (LR_pure and LBW_pure). The genetic variability was estimated using the number of founders, pedigree completeness, number of males and females in reproduction and age distribution, generation interval (GI), inbreeding coefficient (F) and effective population size (Ne). The highest average pedigree completeness values in the second generations of the old genotype LR and LBW were 100%. Higher ages of females in the populations under conservation were related to a higher GI and their longer life expectancy. In 2021, the reproductive age of bulls used for insemination within these populations ranged from 5.1 to 27.8 years. The proportions of males producing offspring in their older age indicate that the semen was used from the national gene bank of commercial artificial insemination companies. The GI (>5) in LR and LBW females was higher than that in LRWP and LBWP. The analysis of the data over the 15-year period showed that the GI of males in LRWP and LBWP decreased equally by 38%, while in LR_pure population, it increased by 80%. A high (9.24%) average inbreeding coefficient (F) was found in inbred animals of LR_pure population, while in LBW_pure, it was 5.35% in 2021. The coefficient of inbreeding varied within the different cattle populations. In the open LR population, it ranged from 1.48% to 2.7%, while in the LRWP population, it fell between 2.12% and 3.72%. The lowest effective population size (Ne) concerning the rate of inbreeding was observed in LBW_pure (23) and LR_pure (59), with the highest Ne identified in the LBWP population (462). When considering Ne based on the number of parents, LR_pure displayed the lowest Ne (42), while the highest Ne was found in LBWP (4449). An analysis of local cattle populations reveals that LR faces the most critical situation. This particular population has been steadily declining for a number of years, necessitating additional measures and efforts to safeguard the LR's ancestral genetic makeup. The results of the LBWP analysis also highlight a concerning trend. Even in very large populations with open breeding programs, the effective population size per generation can experience a significant decrease.

Approaching autozygosity in a small pedigree of Gochu Asturcelta pigs

BACKGROUND

In spite of the availability of single nucleotide polymorphism (SNP) array data, differentiation between observed homozygosity and that caused by mating between relatives (autozygosity) introduces major difficulties. Homozygosity estimators show large variation due to different causes, namely, Mendelian sampling, population structure, and differences among chromosomes. Therefore, the ascertainment of how inbreeding is reflected in the genome is still an issue. The aim of this research was to study the usefulness of genomic information for the assessment of genetic diversity in the highly endangered Gochu Asturcelta pig breed. Pedigree depth varied from 0 (founders) to 4 equivalent discrete generations (t). Four homozygosity parameters (runs of homozygosity, FROH; heterozygosity-rich regions, FHRR; Li and Horvitz's, FLH; and Yang and colleague's FYAN) were computed for each individual, adjusted for the variability in the base population (BP; six individuals) and further jackknifed over autosomes. Individual increases in homozygosity (depending on t) and increases in pairwise homozygosity (i.e., increase in the parents' mean) were computed for each individual in the pedigree, and effective population size (Ne) was computed for five subpopulations (cohorts). Genealogical parameters (individual inbreeding, individual increase in inbreeding, and Ne) were used for comparisons.

RESULTS

The mean F was 0.120 ± 0.074 and the mean BP-adjusted homozygosity ranged from 0.099 ± 0.081 (FLH) to 0.152 ± 0.075 (FYAN). After jackknifing, the mean values were slightly lower. The increase in pairwise homozygosity tended to be twofold higher than the corresponding individual increase in homozygosity values. When compared with genealogical estimates, estimates of Ne obtained using FYAN tended to have low root-mean-squared errors. However, Ne estimates based on increases in pairwise homozygosity using both FROH and FHRR estimates of genomic inbreeding had lower root-mean-squared errors.

CONCLUSIONS

Parameters characterizing homozygosity may not accurately depict losses of variability in small populations in which breeding policy prohibits matings between close relatives. After BP adjustment, the performance of FROH and FHRR was highly consistent. Assuming that an increase in homozygosity depends only on pedigree depth can lead to underestimating it in populations with shallow pedigrees. An increase in pairwise homozygosity computed from either FROH or FHRR is a promising approach for characterizing autozygosity.

The Relevance of the Expected Value of the Proportion of Arabian Genes in Genetic Evaluations for Eventing Competitions

The Arabian horse is a generally reliable sport horse, and continues to be a remarkable endurance horse, so the relevance of the expected value of the proportion of Arabian genes (EV%AG) in horses participating in eventing could be a relevant factor. A total of 1089 horses participating in eventing (8866 records) were used. A GLM revealed that the EV%AG was significant in dressage, show jumping and cross-country. A BLUP genetic evaluation was computed with five genetic models (without the EV%AG (0) using as a covariate (A), as a fixed effect (B), with variance heterogeneity, and in genetic groups without (C) and with (D)). Dressage heritability ranged from 0.103 to 0.210, show jumping ranged from 0.117 to 0.203 and cross-country ranged from 0.070 to 0.099. The lowest DIC value was used as a criterion of fitness. The best fits (those which included variance heterogeneity) showed fewer than two points of difference in DIC values. The highest average estimated breeding value in dressage, show jumping and cross-country was found for horses with an expected value of the proportion of Arabian genes of 0%, ≥1% to <25%, and 100%, respectively. Therefore, the best way to model the EV%AG effect seems to be by considering the variance heterogeneity.

Effects of Selection on Breed Contribution in the Caballo de Deporte Español

The equine breeding industry for sport’s performance has evolved into a fairly profitable economic activity. In particular, the Caballo de Deporte Español (CDE) is bred for different disciplines with a special focus on Show Jumping. The main aim of this study was to determine the effects of 15 years of selection and to find out whether it has been effective and adequate regarding the CDE main breeding objectives. The whole pedigree of 19,045 horses registered as CDE was used, comprising 47,884 animals (18,799 males and 29,085 females). An analysis performed to check for the pedigree completeness level yielded a number of equivalent complete generations (t) equal to 1.95, an average generation interval (GI) of 10.87 years, mean inbreeding coefficient (F) of 0.32%, an average relatedness coefficient (AR) of 0.09% and an effective population size (Ne) of 204. For the analyses, animals were divided into fourteen breed groups. Additionally, in order to study the evolution of these breeds over time and their influence on CDE pedigree, five different periods were considered according to the year of birth of the animals. Performance data used in the analyses were the estimated breeding values (EBV) of the Show Jumping sport discipline of 12,197 horses in the CDE pedigree, available from the 2020 routine genetic evaluations of the CDE breeding program (starting in 2004). Dressage and Eventing EBV values were also assessed. Results showed values of F higher than expected under random mating; this pointed to some degree of inbred matings. With regard to the evolution of breeding values, we found that, in general, EBVs of offspring were higher than the EBVs of parents. Notwithstanding, there is still a need for improvement in population management and the coordination of the breeders to get higher responses but controlling the loss of genetic diversity in the CDE breed.

One Hundred Years of Coat Colour Influences on Genetic Diversity in the Process of Development of a Composite Horse Breed

Genetic diversity and demographic parameters were computed to evaluate the historic effects of coat colour segregation in the process of configuration of the Hispano-Arabian horse (Há). Pedigree records from 207,100 individuals born between 1884 and 2019 were used. Although coat colour is not a determinant for the admission of Hispano-Arabian individuals as apt for breeding, it may provide a representative visual insight into the gene contribution of Spanish Purebred horses (PRE), given many of the dilution genes described in Há are not present in the Arabian Purebred breed (PRá). The lack of consideration of coat colour inheritance patterns by the entities in charge of individual registration and the dodging behaviour of breeders towards the historic banning policies, may have acted as a buffer for diversity loss (lower than 8%). Inbreeding levels ranged from 1.81% in smokey cream horses to 8.80 for white horses. Contextually, crossbred breeding may increase the likelihood for double dilute combinations to occur as denoted by the increased number of Há horses displaying Pearl coats (53 Há against 3 PRE and 0 PRá). Bans against certain coat colours and patterns may have prevented an appropriate registration of genealogical information from the 4th generation onwards for decades. This may have brought about the elongation of generation intervals. Breeder tastes may have returned to the formerly officially-recognised coat colours (Grey and Bay) and Chestnut/Sorrel. However, coat colour conditioning effects must be evaluated timely for relatively short specific periods, as these may describe cyclic patterns already described in owners’ and breeders’ tastes over the centuries.

Discriminant Canonical Analysis of the Contribution of Spanish and Arabian Purebred Horses to the Genetic Diversity and Population Structure of Hispano-Arabian Horses

Simple Summary

The demographic and genetic diversity structure and the contributions of Spanish (PRE) and Arabian Purebred (PRá) horses to the process of conformation of the Hispano-Arabian (Há) horse breed were evaluated. Genetic diversity parameters (inbreeding coefficient, genetic conservation index, coancestry coefficient, non-random mating degree, relatedness coefficient, maximum, complete, and equivalent generations, and number of offspring) were evaluated using a discriminant canonical analysis to determine the partial contributions of each founder breed to the development of Há horse breed. The calculation of Nei genetic distances suggests the overlapping could be estimated in 29.55% of the gene pool of the Há having a PRE origin while 70.45% of the gene pool of the breed may derive from a PRá origin. Although a progressive loss of founder representation may have occurred, breeding strategies implemented considering mating between animals with the highest genetic conservation indices (GCI) may compensate for its effects.

Abstract

Genetic diversity and population structure were analyzed using the historical and current pedigree information of the Arabian (PRá), Spanish Purebred (PRE), and Hispano-Arabian (Há) horse breeds. Genetic diversity parameters were computed and a canonical discriminant analysis was used to determine the contributions of ancestor breeds to the genetic diversity of the Há horse. Pedigree records were available for 207,100 animals born between 1884 and 2019. Nei’s distances and the equivalent subpopulations number indicated the existence of a highly structured, integrated population for the Há breed, which is more closely genetically related to PRá than PRE horses. An increase in the length of the generation interval might be an effective solution to reduce the increase in inbreeding found in the studied breeds (8.44%, 8.50%, and 2.89%, for PRá, PRE, and Há, respectively). Wright’s fixation statistics indicated slight interherd inbreeding. Pedigree completeness suggested genetic parameters were highly reliable. High GCI levels found for number of founders and non-founders and their relationship to the evolution of inbreeding permit controlling potential deleterious negative effects from excessively frequent mating between interrelated individuals. For instance, the use of individuals presenting high GCI may balance founders’ gene contributions and consequently preserve genetic diversity levels (current genetic diversity loss in PRá, PRE, and Há is 6%, 7%, and 4%, respectively).

Relationships between Thoroughbreds’ Contribution in the Pedigree and the Level of Fearfulness and Performance in Warmblood Stallions

The aim of this study was to test the hypothesis that Thoroughbred (TB) ancestors blood percentage in a pedigree is a potential factor differentiating the level of fearfulness (based on behavior and heart rate (HR)) and performance scores in warmblood stallions. A total of 267 three-year-old warmblood stallions were subjected to the novel-object based test during their participation in the performance tests’ program. The effect of ancestors percentage contribution in the pedigree was analyzed by taking into account horses of four TB (<25.01, 25.01-50.00, 50.01-75.00, >75.00%) ancestor groups. It was found that the stallions with a higher proportion (>75.00%) of TB blood revealed higher values of HR measured at the riding hall just before the start of the fearfulness test. The results showed that the higher proportions of TB ancestors blood in the stallions’ pedigree were correlated with the higher values for total time to pass novel objects when led by a handler during the fearfulness test. Negative correlations were found between the TB percentage and the scores for character and trot evaluated by the trainer as well as rideability evaluated by the test rider. The fi ndings imply that a high proportion of Thoroughbreds in the pedigree may negatively infl uence some performance traits and increase fear reactivity in warmblood stallions.

Population study of the Pura Raza Español Horse regarding its coat colour

Coat colour has always been a valuable trait for horse breeders. However, preferences for this feature have changed over the years. In this research, the Pura Raza Español horse (PRE) population was divided into four subpopulations (Grey, Bay, Black and Others), according to the most frequent coat colours and those of their ancestors. The purpose was to analyse genetic variability, reproductive parameters and distances among subpopulations during three key periods in the history of the breed: before 1960, from 1960 to 2000 and after 2000. The subpopulations composed of animals with ancestors with the same coat colour showed higher values of recent inbreeding (ranging from 7.13% to 10.44%) and a greater Nei’s minimum distance between them, as a result of more inbred matings than those carried out in families with members with different coat colours. Non-pure subpopulations also showed more similar recent inbreeding values (between 6.63% and 6.74%). Strikingly, the productive life of Pure bay, Pure black and other subpopulations with minority coat colours was considerably longer (10.79, 10.08 and 9.11 years, respectively) compared to the values of grey PRE horses (6.01 and 7.98 years), which is the subpopulation with the highest census. These results, together with shorter generation intervals of black stallion-offspring (5.51 years via father-son and 6.39 years via father-daughter) and the fact that this coat colour was not present in the breed until two decades ago, highlight the recent trend towards the breeding of black animals.

Assessing performance of single-sample molecular genetic methods to estimate effective population size: empirical evidence from the endangered Gochu Asturcelta pig breed

Estimating effective population size (N_e) using linkage disequilibrium (LD) information (N_e(LD)) has the operational advantage of using a single sample. However, N_e(LD) estimates assume discrete generations and its performance are constrained by demographic issues. However, such concerns have received little empirical attention so far. The pedigree of the endangered Gochu Asturcelta pig breed includes individuals classified into discrete filial generations and individuals with generations overlap. Up to 780 individuals were typed with a set of 17 microsatellites. Performance of N_e(LD) was compared with N_e estimates obtained using genealogical information, molecular coancestry (N_e(M)) and a temporal (two‐sample) method (N_e(JR)). Molecular‐based estimates of N_e exceeded those obtained using pedigree data. Estimates of N_e(LD) for filial generations F3 and F4 (17.0 and 17.3, respectively) were lower and steadier than those obtained using yearly or biannual samplings. N_e(LD) estimated for samples including generations overlap could only be compared with those obtained for the discrete filial generations when sampling span approached a generation interval and demographic correction for bias was applied. Single‐sample N_e(M) estimates were lower than their N_e(LD) counterparts. N_e(M) estimates are likely to partially reflect the number of founders rather than population size. In any case, estimates of LD and molecular coancestry tend to covary and, therefore, N_e(M) and N_e(LD) can hardly be considered independent. Demographically adjusted estimates of N_e(JR) and N_e(LD) took comparable values when: (1) the two samples used for the former were separated by one equivalent to discrete generations in the pedigree and (2) sampling span used for the latter approached a generation interval. Overall, the empirical evidence given in this study suggested that the advantage of using single‐sample methods to obtain molecular‐based estimates of N_e is not clear in operational terms. Estimates of N_e obtained using methods based in molecular information should be interpreted with caution.

Crossbreed genetic performance study in the eventing horse competition

Eventing is an equestrian discipline combining dressage, show jumping and cross-country exercises. The Spanish Sport Horse (SSH) was used in this study as an example to develop a method to ascertain the influence of parental breeds on particular performances by linking their parental genetic contribution with the individual phenotype value computing the optimal breed contribution for each trait evaluated in eventing. Data included 1220 eventing records from 210 SSH animals. The genetic contribution of six main founder populations were computed for each SSH: the Spanish Purebreed (SPB), Arab Horse (A), Thoroughbred (TB), Selle Français (SF), German breeds (G) and North(-west) European (N) breeds. For this analysis we used BLUP (best linear unbiased predictor) animal models for five defined traits and for the total score. The genetic contribution as linear and quadratic adjustment and the interaction between genetic contributions were included as covariates in the genetic model to separate the possible combining ability from the breeding values. The optimal breed genetic contribution for each trait was maximised by using a simplex method. The best combination for conformation was 51% SPB and 49% N, for cross-country aptitude 23% SPB and 77% N, for dressage 48% SF and 52% N, and for show jumping 64% A and 36% N. For cross-country the best performance was predicted using 100% N genetic contribution and for the total score of 24% SF and 76% N. The combination of two breeds (group of breeds) seemed to be the best option for most of the traits. Although the results should be interpreted with caution, the importance of this paper is that it can be considered as a starting point of the analysis. The methodology applied here performed nicely in searching for the best contribution of several breeds to find the best combination for particular interests and could, therefore, be useful for other species/populations.

Genetic structure and gene flows within horses: a genealogical study at the french population scale

Since horse breeds constitute populations submitted to variable and multiple outcrossing events, we analyzed the genetic structure and gene flows considering horses raised in France. We used genealogical data, with a reference population of 547,620 horses born in France between 2002 and 2011, grouped according to 55 breed origins. On average, individuals had 6.3 equivalent generations known. Considering different population levels, fixation index decreased from an overall species F_IT of 1.37%, to an average of −0.07% when considering the 55 origins, showing that most horse breeds constitute populations without genetic structure. We illustrate the complexity of gene flows existing among horse breeds, a few populations being closed to foreign influence, most, however, being submitted to various levels of introgression. In particular, Thoroughbred and Arab breeds are largely used as introgression sources, since those two populations explain together 26% of founder origins within the overall horse population. When compared with molecular data, breeds with a small level of coancestry also showed low genetic distance; the gene pool of the breeds was probably impacted by their reproducer exchanges.

Founder and present maternal diversity in two endangered Spanish horse breeds assessed via pedigree and mitochondrial DNA information

Pedigree information and 179 mtDNA sequences from two endangered Spanish horse breeds, the Asturcón pony (143) and the Mallorquí horse (36), were analysed to asses: (i) the pedigree and molecular maternal genetic diversity of the two breeds; (ii) the concordance between the dam lines recorded in the corresponding studbooks and the mtDNA haplotypes identified; and (iii) to assess the losses of maternal genetic variability occurred from the foundation of the studbooks to present. Up to 50 Asturcón and 18 Mallorquí founder dam lines were identified in the studbooks analysed. Up to 315 Asturcón mares and 51 Mallorquí mares that foaled in the last 5 years of recording formed a reference population. Only 35 Asturcón and 13 Mallorquí founder dam lines were represented in their reference populations. Sequences from a total of 38 Asturcón and 12 Mallorquí dam lines could be obtained. The 179 sequences obtained gave 15 different haplotypes, 11 and 9 of them being identified, respectively, in the Asturcón pony and in the Mallorquí horse. Five different haplotypes (roughly two-thirds of the sequences) were shared by the two horse breeds. Most dam lines analysed had a single mtDNA haplotype. However, more than one haplotype was detected within eight of the dam lines in Asturcón pony. The found inconsistencies are likely to result from deficiencies in genebank management. The maternal N(e) (mN(e)) computed using the dam line information was higher in the Asturcón pony (20.5) than in the Mallorquí horse (15.9), while these figures were on the opposite direction for the haplotypic line information (6.4 and 9.4, respectively). The ratio of the computed mN(e) values to the actual number of founder dam lines were always higher in the Mallorquí horse probably due to a more balanced distribution of individuals kept for reproduction among studs. Consequences for the conservation programmes of the analysed breeds are discussed.

Estimation of effective population size from the rate of coancestry in pedigreed populations

We introduce a simple method to estimate effective population size from increase in coancestry (Δc(jk)) for all pairs of individuals j and k in a reference subpopulation. An increase in pairwise coancestry for any pair of individuals j and k can be defined assuming that a hypothetical mating between them would give an individual with an inbreeding coefficient equal to c(jk), where c(jk) is the coancestry coefficient between the individuals j and k. The equivalent measure to discrete generations value (g(jk)) corresponding to the individual jk can be computed by averaging discrete equivalents generations of its parents (g(j) and g(k)). The mean increase in coancestry for all pairs of individuals in a reference subpopulation can be used to estimate a realized effective population size based on coancestries that would provide information on the effective size of a population under random mating. Performance of the new parameter was tested on simulated and empirical (horse) populations with different mating strategies and population structures. The routines needed to compute the introduced parameters have been included in a new version of the program ENDOG.