Genetic Diversity of Lusitano Horse in Brazil Using Pedigree Information
Total Page:16
File Type:pdf, Size:1020Kb
Journal of Equine Veterinary Science 69 (2018) 149e158 Contents lists available at ScienceDirect Journal of Equine Veterinary Science journal homepage: www.j-evs.com Original Research Genetic Diversity of Lusitano Horse in Brazil Using Pedigree Information * Ricardo Antonio da Silva Faria a, , Antonio Pedro Andrade Vicente c, d, Rute Isabel Duarte Guedes dos Santos e, f, Amanda Marchi Maiorano a, Rogerio Abdallah Curi b, Luis Artur Loyola Chardulo b, Josineudson Augusto II Vasconcelos Silva b a Universidade Estadual Paulista (Unesp), Faculdade de Ci^encias Agrarias e Veterinarias, CEP 14.884-900, Jaboticabal, Sao~ Paulo, Brazil b Universidade Estadual Paulista (Unesp), Faculdade de Medicina Veterinaria e Zootecnia, CEP 18.618-307, Botucatu, Sao~ Paulo, Brasil c Escola Superior Agraria do Instituto Politecnico de Santarem, Apartado 310, 2001-904 Santarem, Portugal d CIISA - Faculdade de Medicina Veterinaria, Universidade de Lisboa, 1300-477 Lisboa, Portugal e Instituto Politecnico de Portalegre - Escola Superior Agraria de Elvas, 7350-903 Elvas, Portugal f VALORIZA - Centro de Investigaçao~ para a Valorizaçao~ de Recursos Endogenos, 7300- 555 Portalegre, Portugal article info abstract Article history: This study aimed to evaluate population parameters and to describe the genetic diversity of the Lusitano Received 21 March 2018 breed in Brazil using pedigree data. Two populations were evaluated: total population (TP) containing Received in revised form 18,922 animals, and reference population (RP) composed of a part of TP containing 8,329 animals, rep- 12 July 2018 resenting the last generation. The generation interval (10.1 ± 5.1 years) was in the range for horse pop- Accepted 14 July 2018 ulations. Pedigree completeness in RP shows almost 100% filling in the three most recent generations, Available online 3 August 2018 indicating improvement in the pedigree data and accuracy of the results, and the inbreeding coefficient (4.46%) and average relatedness (5.97%) for RP, indicating control on the part of breeders. The effective Keywords: Ancestors population size was 89 (TP) and 90 (RP). The effective number of founders (fe) were 33 and 29, effective Demographic characterization number of ancestors (fa) were 30 and 26, and effective number of founder genomes (fg) were 19 and 15 for Effective population size TP and RP, respectively, indicating a reduction of genetic variability in the last generations. The total Founders number of ancestors that explains 100% of the genetic diversity in the Lusitano breed in Brazil was 427 (TP) Inbreeding and 341 (RP). The reproductive parameters, probabilities of gene origin showing loss of variability in the last generations, and the genetic contributions of ancestors suggest the need to monitor genetic diversity over time in breeding programs to allow control of the next generations and to increase their variability. © 2018 Elsevier Inc. All rights reserved. 1. Introduction due to their use for working cattle and cutting bulls [1]. Although records date back to 1824, the studbook of the Lusitano breed was The Puro Sangue Lusitano (PSL), also known as the Lusitano horse, officially established in 1967 [2], a time when the main Iberian represents a breed with a rich and ancient history, originated in breeds, Portuguese and Spanish, were separated and received the southern Portugal on the Iberian Peninsula. Its evolution is essentially denomination PSL in Portugal and Pura Raza Espanola~ in Spain. The main studbook of the Lusitano breed is managed in Portugal by the Portuguese Lusitano Horse Breeder Association (APSL in the Animal welfare/ethical statement: No animal was used in the present research Portuguese acronym), which has the largest number of birth re- “Genetic diversity of Lusitano horse in Brazil using pedigree information”; only the cords of the breed in the world [2]. Brazil is the country with the pedigree data of the animals were used. second largest number of animal records of this breed. This statistic was determined after the reciprocity partnership between APSL Conflict of interest statement: The authors wish to confirm that there are no known conflicts of interest associated with this publication and there has been no signif- and the Brazilian Lusitano Horse Breeder Association (ABPSL in the icant financial support for this work that could have influenced its outcome. Portuguese acronym) was established in 1991. This agreement * Corresponding author at: Ricardo Antonio da Silva Faria, FMVZ - Unesp, DMNA - permitted all Lusitano horses registered in the Brazilian studbook Fazenda Experimental Lageado, Rua Jose Barbosa de Barros, n 1780 CEP: 18.618- to be equally registered in the Portuguese studbook, which contains 307, Botucatu, SP, Brasil. E-mail address: [email protected] (R.A. da Silva Faria). the birth records of all countries. https://doi.org/10.1016/j.jevs.2018.07.009 0737-0806/© 2018 Elsevier Inc. All rights reserved. 150 R.A. da Silva Faria et al. / Journal of Equine Veterinary Science 69 (2018) 149e158 Studies investigating population parameters of the Lusitano proportion of known ancestors in each ascending generation per breed have been conducted only with data from the APSL studbook descendant and was calculated as proposed by MacCluer (Portugal) for all Lusitano horses registered in the world [2,3]. Thus, et al [4], in which ancestors without progenitors in the pedigree data on the population statistics and genetic diversity of the Lusi- were considered founders [21];(2)basedonthemean tano breed in Brazil are lacking. number of generation equivalents (GE), computed over the sum Parameters such as pedigree completeness [4], generation inter- of all known ancestors by calculating (½)n,wheren is the val (GI), inbreeding coefficient [5], and probability of gene origin [6] number of generations between the animal and each known are important to design strategies for the selection and improvement ancestor [7]. of animals. The data generated permit to verify genetic diversity and its changes over time [7]. The complete or partial results of genetic fi diversity and population parameters are reported in the literature for 2.4. Parameters Related to the Inbreeding Coef cient various horse breeds, including Andalusian [8], Holstein [9], Paint fi fi Horse [10], Hanoverian [11], Spanish Arab Horse [12],Lipizzaner[13], The inbreeding coef cient (F) de nes the probability that an and Old Kladruber [14]. In Brazil, studies were conducted on the individual has two identical alleles by descent and was calculated Pantaneiro [15],Mangalarga[16], Campolina [17],BrazilianSport using the algorithm proposed by Meuwissen and Luo [23]. fi Horse [18], and Quarter Horse [19] breeds. The average relatedness (AR) coef cient of each animal was This study aimed to evaluate the population parameters of the described by Gutierrez and Goyache [21] as the representation of a Lusitano horse in Brazil to observe the genetic diversity and to given animal in the pedigree and was obtained based on the contribute to the knowledge and development of this breed in probability that a randomly selected allele of the population be- longs to a given animal. Brazil and to compare our results with the ones obtained by D different authors for the worldwide population [2,3], providing the The increase in inbreeding ( F) for each generation was ob- information necessary to implement a breeding program of the tained as follows: breed in Brazil. DF ¼ðFt À FtÀ1Þ=ð1 À FtÀ1Þ (1) 2. Materials and Methods where Ft and FtÀ1 are the average inbreeding in generations t and t À 1, respectively. 2.1. Data and Computer Programs According to Falconer and Mackay [24], the effective size of a population (N )isdefined as the number of individuals of a fi e The pedigree le of the Lusitano breed containing information population with a nonideal structure that would give origin to a about the animal's name, sire and dam, sex, date of birth, and origin certain rate of consanguinity if its structure were ideal (e.g., equal was provided by ABPSL. Animals born and registered in Brazil be- number of males and females, absence of selection, random tween 1912 and 2012 were considered and made up the total pop- matings, etc.). Using DF, Ne was estimated considering Ne ¼ 1/2DF, ulation (TP) of 18,922 animals (48.1% males, stallions and geldings) which represents the number of animals that equally contribute to and the reference population (RP) consisting of 8,329 animals (51.0% the next generation and that would promote a similar increase in males, stallions and geldings), corresponding to the birth records inbreeding in the population studied [21]. In addition, Ne was also from 2003 to 2012. RP was used as the reference of the active pop- calculated based on the individual increase in inbreeding as sug- ulation, representing the last generation, which was equal to one gested by Gutierrez et al [25] and was used for the calculation of mean GI (10.1 years) and contained animals (stallions and mares) genetic drift. that potentially could transmit their genes to the next generation. Data preparation and statistical analysis were performed with the MEAN and FREQ procedures of the SAS program [20]. The population 2.5. Probability of Gene Origin and Genetic Drift and reproductive parameters, probability of gene origin, and genetic diversity were obtained with the ENDOG V4.8 program [21]. The effective number of founders (fe) is obtained by measuring the contributions of the most influential founders. Lacy [26] defines 2.2. Reproductive Parameters and Generation Interval fe as the expected number of founders that contribute equally and produce the same genetic diversity in the population studied. This The existence of imported animals registered in the studbook of parameter was calculated using the formula ABPSL permitted to observe their use as breeding stock and to , compare it with that of sires and dams born in Brazil. The co- Xf ¼ 2 efficients of variation for the mean number of offspring, age at first fe 1 qk (2) ¼ and last progeny, and time in reproduction obtained for stallions k 1 and mares explain the high standard deviation between the where q is the probability of gene origin of founder k.