INTRODUCTION
In Costa Rica, the populations of Zebu cattle (Bos indicus) trace back to 1911. Starting in 1946, the Costa Rican government initiated the genealogical recording (GR) of domestic cattle populations. In 1989, the government transferred the GR to the registration and specialized breeders' associations, which later gave birth to the Costa Rican Zebu (ASOCEBU) Cattle Breeders Association (Quirós, 2006; Vásquez-Loaiza and Molina-Castro, 2020). The ASOCEBU coordinates and controls the GR and production data of Zebu cattle populations under its administration. Additionally, this association oversees the development of breeding programs in accordance with the selection criteria and objectives for each breed of cattle (ASOCEBU, 2007; Camacho et al., 2019).
Pedigree records include the genealogical information of the entire population. Their main functions are to ensure breed purity, support conservation, and breeding programs, and provide a unique and permanent identification to each of the animals in the system. Pedigree analysis is used to: identify ancestors and founders, and their contributions (%) to variability; evaluate the levels of genetic variability; characterize the genealogical structure and demographics of a population; analyze the levels, trends, and possible effects of inbreeding. The population genetic parameters generated by pedigree analysis are of importance to implement conservation and / or genetic improvement programs (Melka et al., 2013; Cañas-Álvarez et al., 2014; Pereira et al., 2016; Santana et al., 2016; Chiara et al., 2019).
The genetic evaluations of Costa Rican Brahman cattle started in 2004 (Camacho and Cruz, 2004). The Best Linear Unbiased Prediction (BLUP), which results from the genetic evaluations, favors the selection of related animals, thus increasing inbreeding. Furthermore, inbreeding and kinship levels are involved in the process of genetic evaluations using BLUP (Santana et al., 2010; Pereira et al., 2016). Therefore, the objective of this study was to analyze the population structure of the Costa Rican Brahman cattle populations using pedigree analysis.
MATERIALS AND METHODS
This study analyzed the GR of Brahman cattle of Costa Rica. These records included information from 98,749 animals (40.4% males and 59.6% females) born between 1965 and 2020, from 5,414 bulls and 29,880 mother cows, with an average of 16.2 and 2.8 calves, respectively. The analysis was carried out using the ENDOG software (Gutiérrez and Goyache, 2005). The population parameters used to analyze the structure and genetic variability were the following:
Pedigree integrity (Faria et al., 2009; Santana et al., 2016), as the proportion of known ancestors up to the fourth generation; the number of complete generations (NCG), maximum generations traced (NMGT), and equivalent complete generations (NECG).
Inbreeding coefficient (F). The inbreeding coefficient of each individual was estimated, and its trend over birth years was generated (Falconer and Mackay, 1996; Parland et al., 2007).
Generation interval (GI). The GI was determined using the mean age of a reproductive animal, at which it is replaced by one of its offspring. The GI was calculated using the four selection paths: father - son, father - daughter, mother - son, and mother - daughter (Ramírez-Valverde et al., 2018; Hagan and Cue, 2020).
Additive genetic relationship (AGR). The AGR was determined from the matrix of additive genetic relationships among all the individuals included in the pedigree. The calculation corresponds to the average value of the coefficients of each individual with the rest of the pedigree (Gutiérrez et al., 2003; Ocampo et al., 2020).
Effective number of founders (fe). Founder, any animal with unknown parents. The fe was defined as the founders that, by contributing in equal degree, would produce the existing genetic diversity in the population (Gutiérrez et al., 2009; Cañas-Álvarez et al., 2014).
Effective number of ancestors (fa). Ancestor, any animal, founder or not, that has contributed to the genetic variability of the population. The fa was defined as the number of ancestors required to explain the total genetic variability of the population; it considers the genetic variability, which is not explained by the contribution of any of its offspring, provided by an animal (Gutiérrez et al., 2009; Cañas-Álvarez et al., 2014).
Effective population size (Ne). The Ne defines the number of breeding stock that could generate the calculated inbreeding and change rate in the genetic variance under the ideal population scheme; it was estimated through three different approaches: NCG, NMGT, and NECG (Gutiérrez et al., 2008; Leroy et al., 2013).
RESULTS AND DISCUSSION
Fig. 1 represents the proportion of known ancestors up to the fourth generation, where the maximum proportions of unknown parents, grandparents, and great-grandparents were 14.1%, 38.1%, and 61.2%, respectively. A total of 15,496 animals were identified as founders and 13,391 as ancestors, with fa and fe equal to 172.0 and 617.8, respectively. The highest contributing ancestor explained 4.1% of the pedigree variability; 155 and 6,151 animals were responsible for 50% and 95% of the pedigree variability, respectively. The GR of Brahman cattle in Costa Rica requires at least three generations of known ancestors to issue a pure breed certificate (ASOCEBU, 2007); the proportion of individuals with unknown parents is attributed to the animals identified as founders. The precision and accuracy of a population structure analysis depend on pedigree integrity and the amount of genealogical information across generations.
Incomplete information results in an approximate generation assignment of individuals, and the calculations of F and Ne are imprecise.
Gutierrez et al.(2003) evaluated eight Spanish breeds and reported lower pedigree integrity than that observed in the present study. In fact, the authors described that pedigree integrity was substantially lower (lower than 20%) for the known ancestors of the second and third generations. However, our results are in agreement with Santana et al. (2012) and Abin et al. (2016), who evaluated two Brazilian and five outh African breeds, respectively. Similarly, Ramírez-Valverde et al. (2018) conducted a study on seven Bos taurus breeds in Mexico and estimated pedigree integrity of similar magnitude. I i Bos indicus breeds, Faria et al. (2009), Faria et al. (2010), Borges et al. (2013), Santana et al. (2014), and Costa et al. (2019) defined an interval for the proportion of known genealogical in the range of ur results.
The average F for the total population was .86%; 27.4% of the animals in the population were inbred, with an average of 3.1%. Fig. 2 resents the proportion of inbred individuals and their inbreeding levels over time. This study observed contrasting trends; the proportion of inbred animals increased as the mean inbreeding coefficient decreased. Previous studies on Brahman cattle have reported fa and fe values ranging from 5.4 to 50 and from 5.4 to 718, respectively. Additionally, average F values ranged from 0.69% to 11.9% (Faria et al., 2010; Santana et al., 2016; Cavani et al., 2018).
The AGR was 0.64, with a maximum value of 4.1 and a positive trend over time (Fig. 3), which correlates with the behavior of the average inbreeding percentage shown in Fig. 2. The knowledge of the AGR allows the designing of breeding schemes while maintaining certain levels of F in the offspring. The results observed in Figures 2 and 3 can be attributed to the use of sires related to a large number of animals in the population; selection base don BLUP that increases the probability of selecting related animals; and advances in reproductive technologies that reduce the number of parents required to breed the next generation. In the relation of fe and fa, the AGR indicates the percentage of the population that originated from a founder (Gutiérrez and Goyache, 2005).
The amount and trend of F, which depends on AGR and Ne, represent the genetic variability evolution over time. Pedigree integrity is related to the NCG, NMGT, and NECG (Table 1), and the population structure results from the selection and breeding schemes adopted by the breeders. The genetic relationships between the founders and the fe expose the initial genetic variability; the founders' contribution to pedigree variability exhibits the gene pool that has been maintained across generations (Biochard et al., 1997); a small number of ancestors that explain the pedigree variability is associated with increases in the levels of AGR and F. The fa comprises the possible causes of loss in genetic variability. Generally, fe > fa, as this difference increases, the participation of founders decreases across generations; the fe / fa ratio shows differential breeding schemes and considers the possible bottlenecks experienced by the population. When this ratio is extended, it indicates that most of the ancestors were founders
with the absence of bottlenecks (Bernardes et al.,2016).
The Ne values in Table 1 indicate lower than 1% increases in F. Ne, as a concept, was developed based on ideal population guidelines. It reflects the accumulation of genetic relationships between individuals and allows predicting changes in the levels of F. A reduced Ne is associated with decreased genetic variability, increased mating between related individuals, allele fixation, and reduced selection response (Gutiérrez et al., 2009; Leroy et al., 2013). For the populations under breeding schemes, the recommended Ne values optimize the selection response with a minimum increase in F (Meuwissen et al., 1998).
The GI estimated using the four selection paths were (years ± standar error): father - son 7.4 ± 0.10; father - daughter 7.1 ± 0.09; mother - son
Number of complete generations (NCG), maximum generations traced (NMGT), equivalent complete generations (NECG). Maximum (Max) and average (Mean) values. AF, inbreeding rate change. Ne, effective population size based on the number of known generations.
6.5 ± 0.08; and, mother - daughter 6.4 ± 0.07. The magnitude of the GI can be associated with the production system, selection path, and market. In dairy cattle, Melka et al. (2013) reported, for five specialized breeds in Canada, GIs ranging from 3.64 to 6.63 years, with an average of 5.03 years. Hagan and Cue (2020) evaluated the behavior of the GI overtime in four specialized breeds. At the beginning, they observed GIs ranging from 4.5 to 12 years; after 25 years, the GIs ranged from 3.5 to 6.0 years. Previous studies in Brahman cattle have reported average GIs of 4.0 years (Cavani et al., 2018), ranging from 5.0 to 9.8 years (Santana et al., 2016). In dairy Gyr cattle, reported GIs are greater than 7 years, with maximum values of 12 years (Reis et al., 2010; Santana et al., 2014; Santana et al., 2016). In Nelore populations, the GI ranges from 6.5 to 10 years (Borges et al., 2013; Costa et al., 2019). In Bos taurus beef cattle breeds, previous studies have reported GIs ranging from 5.4 to 7.7 years (Parland et al., 2007), 4.9 to 8.2 years (Ramírez-Valverde et al., 2018), and average values between 6.0 and 6.4 years (Piccoli et al., 2014).
CONCLUSIONS
Based on the effective number, number of generations, and large proportion of known ancestors in the third and fourth generations, the evaluated pedigree presented an integral, complete, and deep structure. The results show a summary of the genetic and reproductive management that breeders have carried out, including levels of additive genetic relationship and indicating that the percentage of inbred animals has increased. However, the mean inbreeding value presented a negative trend and the forecasted increment, given the effective population size, is less than 1%.