SciELO - Scientific Electronic Library Online

 
vol.33 issue1Comparing RAS with and without biofloc: Transcriptional response of immune-related genes in Litopenaeus vannamei post-larvaeGenetic parameters between somatic cell score and production traits for Holstein cattle in Southern Brazil author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista Colombiana de Ciencias Pecuarias

Print version ISSN 0120-0690On-line version ISSN 2256-2958

Rev Colom Cienc Pecua vol.33 no.1 Medellín Jan./Mar. 2020

https://doi.org/10.17533/udea.rccp.v32n4a05 

Original research article

Evaluation of the Romosinuano cattle population structure in Mexico using pedigree analysis*

Evaluación de la estructura poblacional de bovinos Romosinuano en México mediante el análisis de pedigrí

Avaliação da estrutura populacional de bovinos Romosinuano no México usando análise de pedigree

Rafael Núñez-Domínguez1 

Ricardo E Martínez-Rocha1 

Jorge A Hidalgo-Moreno1 

Rodolfo Ramírez-Valverde1 

José G García-Muñiz1 

1Departamento de Zootecnia, Universidad Autónoma Chapingo, Chapingo, México.


Abstract

Background:

Romosinuano cattle breed in Mexico has endured isolation and it is necessary to characterize it in order to facilitate sustainable genetic management.

Objective:

To assess the evolution of the structure and genetic diversity of the Romosinuano breed in Mexico, through pedigree analysis.

Methods:

Pedigree data was obtained from Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical (AMCROLET). The ENDOG program (4.8 version) was used to analyze two datasets, one that includes upgrading from F1 animals (UP) and the other with only straight-bred cattle (SP). For both datasets, three reference populations were defined: 1998-2003 (RP1), 2004-2009 (RP2), and 2010-2017 (RP3). The pedigree included 3,432 animals in UP and 1,518 in SP. Demographic parameters were: Generation interval (GI), equivalent number of generations (EG), pedigree completeness index (PCI), and gene flow among herds. Genetic parameters were: Inbreeding (F) and average relatedness (AR) coefficients, effective population size (Nec), effective number of founders and ancestors, and number of founder genome equivalents.

Results:

The GI varied from 6.10 to 6.54 for UP, and from 6.47 to 7.16 yr for SP. The EG of the UP and SP improved >63% from RP1 to RP3. The PCI increased over time. No nucleus or isolated herds were found. For RP3, F and AR reached 2.08 and 5.12% in the UP, and 2.55 and 5.94% in the SP. For RP3, Nec was 57 in the UP and 45 in the SP. Genetic diversity losses were attributed mainly (>66%) to genetic drift, except for RP3 in the SP (44%).

Conclusions:

A reduction of the genetic diversity has been occurring after the Romosinuano breed association was established in Mexico, and this is mainly due to random loss of genes.

Keywords: effective population size; gene flow; genetic diversity; genetic drift; generation interval; inbreeding; pedigree; population structure; probability of gene origin; Romosinuano cattle

Resumen

Antecedentes:

La raza bovina Romosinuano ha estado prácticamente aislada en México y requiere ser caracterizada para un manejo genético sostenible.

Objetivo:

Evaluar la evolución de la estructura y diversidad genética de la raza Romosinuano en México, mediante el análisis del pedigrí.

Métodos:

Los datos genealógicos provinieron de la Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical (AMCROLET). Los análisis se realizaron con el programa ENDOG (versión 4.8) para dos bases de datos, una que incluyó animales en cruzamiento absorbente (UP) a partir de F1 y la otra con sólo animales puros (SP). Para ambas bases de datos se definieron tres poblaciones de referencia: 1998-2003 (RP1), 2004- 2009 (RP2), y 2010-2017 (RP3). El pedigrí incluyó 3.432 animales en la UP y 1.518 en la SP. Los parámetros demográficos fueron: intervalo generacional (GI), número de generaciones equivalentes (EG), índice de completitud del pedigrí (PCI), y flujo de genes entre hatos. Los parámetros genéticos fueron: coeficientes de consanguinidad (F) y de relación genética aditiva (AR), tamaño efectivo de la población (Nec), número efectivo de fundadores y ancestros, y número equivalente de genomas fundadores.

Resultados:

El GI varió de 6,10 a 6,54 para la UP, y de 6,47 a 7,16 años para la SP. El EG de la UP y la SP mejoró >63%, de RP1 a RP3. El PCI aumentó a través de los años, pero más para la SP que para la UP. No se encontraron hatos núcleo o aislados. Para RP3, F y AR alcanzaron 2,08 y 5,12% en la UP, y 2,55 y 5,94% en la SP. Para RP3, Nec fue 57 en la UP y 45 en la SP. Más de 66% de las pérdidas en diversidad genética se debieron a deriva genética, excepto para RP3 en la UP (44%).

Conclusiones:

una reducción de la diversidad genética ha estado ocurriendo después de que se formó la asociación de criadores de ganado Romosinuano en México, y es debida principalmente a pérdidas aleatorias de genes.

Palabras clave: consanguinidad; deriva genética; diversidad genética; estructura poblacional; flujo de genes; ganado Romosinuano; intervalo generacional; pedigrí; probabilidad de origen del gen; tamaño efectivo de población

Resumo

Antecedentes:

A raça bovina Romosinuano tem estado praticamente isolada no México e precisa ser caracterizada para um manejo genético sustentável.

Objetivo:

Avaliar a evolução da estrutura e diversidade genética da raça Romosinuano no México, através da análise de pedigree.

Métodos:

Os dados genealógicos vieram da Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical (AMCROLET). As análises foram feitas com o programa ENDOG (versão 4.8) para duas bases de dados, uma que incluiu animais em cruzamento absorvente (UP) a partir da F1 e a outra base de dados somente com animais puros (SP). Para ambas bases de dados foram definidas três populações de referência: 1998-2003 (RP1), 2004-2009 (RP2) e 2010-2017 (RP3). O pedigree incluiu 3.432 animais na UP e 1.518 na SP. Os parâmetros demográficos foram: intervalo entre gerações (GI), número de gerações equivalentes (EG), índice de completude do pedigree (PCI), e fluxo de genes entre rebanhos. Os parâmetros genéticos foram: coeficiente de consanguinidade (F) e da relação genética aditiva (AR), tamanho efetivo da população (Nec), número efetivo de fundadores e ancestrais, e número equivalente de genomas fundadores.

Resultados:

O GI variou de 6,10 a 6,54 para a UP, e de 6,47 a 7,16 anos para a SP. EG da UP e a SP melhorou >63%, de RP1 a RP3. O PCI aumentou ao longo dos anos, mas mais para a SP do que para o UP. Não se encontraram rebanhos núcleo ou isolados. Para RP3, F e AR alcançaram 2,08 e 5,12% na UP, e 2,55 e 5,94% na SP. Para RP3, Nec foi 57 na UP e 45 na SP. Mais de 66% das perdas em diversidade genética foram ocasionadas pela deriva genética, exceto para RP3 no UP (44%).

Conclusões:

Depois que a associação da raça Romosinuano foi estabelecida no México, tem ocorrido uma redução da diversidade genética, principalmente devido a perdas aleatórias de genes.

Palavras-chave: consanguinidade; deriva genética; diversidade genética, estrutura populacional; fluxo de genes; intervalo entre gerações; pedigree; probabilidade de origem do gene; Romosinuano; tamanho efetivo da população

Introduction

Romosinuano, a Creole breed developed in Colombia, is based on cattle that arrived from the Iberian Peninsula five centuries ago. According to De Alba (2011), Romosinuano germplasm entered into México as semen. Semen from Colombia was first exported to North Carolina, USA, in the middle of the 20th century. Then, semen was exported to Venezuela in 1982. A research herd was founded in Turrialba, Costa Rica, in 1955 using sires and heifers imported from North Caroline. Subsequently Romosinuano embryos were imported by Florida State University, from Venezuela. The Mexican herd was developed from germplasm of Turrialba, and later from genetics of Florida State University.

Several breeders in the Mexican tropics have used Romosinuano due to ease of management, longevity (Carroll et al., 2011; 2012), fertility (Riley et al., 2007), tolerance to heat, humidity (Scharf et al., 2010), and to the toxic effects of endophytes commonly present in tall fescue (Riley et al., 2016). The Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical (AMCROLET) was established in 1998 (De Alba, 2011), and it is responsible for the herd-book and the performance records. Recently, the first Romosinuano genetic evaluation in Mexico was carried out (AMCROLET, 2016), which would allow breeders to identify superior breeding stock for genetic improvement. However, the selection of breeding stock may result in a smaller population size and in higher inbreeding (Bernardes et al., 2016). This is even more important considering that the Mexican herd of Romosinuano is relatively isolated due to sanitary restrictions that would not allow import genetic material directly from Colombia (SAGARPA, 2017).

The main goal of conservation programs is to control the effects of genetic drift in order to maintain the genetic variability of the population and to keep a low rate of inbreeding (Jamieson and Allendorf et al., 2012; Hammerly et al., 2013; Toro et al., 2014). Genetic variability works simmilar to an ´insurance´ for the population to face potential environmental and marketing changes; at the same time, it reduces the unfavorable effects on performance due to inbreeding (Burrow et al., 1993; Santana et al., 2012). The studies on genetic variability of populations have been conducted using pedigree analysis or genetic markers. Pedigree information provides a useful tool for conservation programs to maintain genetic variability and minimize inbreeding (Pinheiro et al., 2013; Santana et al., 2016; Sheikhlou and Abbasi et al., 2016).

Considering the potential of Romosinuano cattle to improve beef production under tropical conditions, the availability of predicted breeding values, and the fact that it is a breed with a small population size and relatively isolated, it is imperative to characterize the structure of this population and evaluate changes in its genetic variability. Therefore, the objective of this study was to assess the evolution of the structure and genetic diversity of Romosinuano breed in Mexico through pedigree analysis.

Materials and methods

Description of the data

Pedigree data (identification of animal, sire and dam, sex, birth date and herd) was obtained from AMCROLET. The database included 4,809 animals born between 1953 and 2017, belong ing to 39 herds located mainly in the Southeast region of Mexico.

Given that the Romosinuano population has been developing from upgrading, pedigree analyses were performed for two datasets: one that includes upgrading from F1 animals (UP; n = 3,432) and the other with only straightbred cattle (SP; n = 1,518); thus, SP is included in UP. For both datasets, three reference populations were defined: 1998-2003 (RP1), 2004-2009 (RP2), and 2010-2017 (RP3) to describe the time-trend of pedigree parameters for approximately three generation intervals after AMCROLET was created.

Pedigree analyses

The analyses of genealogical informa tion were carried out using ENDOG ver sion 4.8 (Gutiérrez and Goyache, 2005).

Demographic parameters

Generation interval (GI). The GI (average age of the parents when their replacements were born) was obtained for the four gametic paths (sire-son, sire-daughter, dam-son, dam-daughter).

Pedigree completeness level. To evaluate the level of integrity of the pedigree, the number of complete (CG), maximum (MG) and equivalent (EG) generations were calculated. The CG implies that the 2n ancestors of an animal are known, where n is the number of fully traced generations, whereas MG counts the number of generations that separate the animal from the oldest ancestor (Gutiérrez y Goyache et al., 2005). For a given animal, EG was calculated as the sum of (1/2)n coefficients of all known ancestors, where n is the number of generations separating the individual to each known ancestor (Maignel et al., 1996). The percentage of known ancestors during the last three generations was calculated, a period that encompasses the time from AMCROLET formation. Finally, the pedigree completeness index (PCI; MacCluer et al., 1983) was generated for each generation.

Gene flow among herds. The role of each herd was evaluated according to the contribution of the sires to the population (Vassallo et al., 1986). Herds were grouped as: a) nucleus, if breeders only use their own sires, which may also be sent to other herds; b) multiplier, when breeders use sires from other herds and provide them to other herds; c) commercial, if breeders use sires from other herds, but do not deliver sires; and d) isolated, when breeders only use own sires and do not provide them to other herds.

Genetic parameters

Inbreeding (F) and average relatedness coefficients (AR). The F was calculated using the algorithm of Meuwissen and Luo et al.(1992), assuming unrelated and noninbred founders. For each animal, AR was defined as the probability that an allele randomly chosen in the population belongs to a particular animal (Gutiérrez and Goyache et al., 2005). The AR for a particular animal is equivalent to the average of the coefficients in the row of the numerator relationship matrix.

Effective population size (Ne). The Ne was estimated using two methods. The first was based on the individual increase in inbreeding (ΔFi) and was computed as follows (Gutiérrez et al., 2009):

where, t is the number of equivalent generations. The Ne was calculated from the average of the inbreeding rate of the n animals included, as .

The second method was based on increases in coancestry (Δcjk) between any pair of individuals j and k (Cervantes et al., 2011), as:

where, cjk is the inbreeding of a progeny from j and k, and gj and gk are the discrete equivalent generation for the parents. Then, Ne was estimated by averaging the rate of coancestry (Nec) for all pairs of the individuals as.

Probability of gene origin. To assess the amount and kind of genetic diversity losses in the populations, the number of founders (f), effective number of founders (fe), effective number of ancestors (fa), and founder genome equivalents (fg) were obtained. The fe counts founders with unknown genetic relationships with other animals in the pedigree, except to its descendant (Lacy, 1989). Preservation of genetic diversity of founders toward the studied population may be evaluated by its contributions, that must add up to one (Boichard et al., 1997).

The fe is the number of equally contributing founders that would be expected to produce the same amount of genetic diversity as observed in the population under study (Lacy, 1989), and was calculated as:

where, pk is the probability of gene origin of founder k to the reference population. The fa is the minimum number of ancestors (founders or not) needed to explain the full genetic diversity of the population (Boichard et al., 1997), and it was obtained as:

where, pi is the marginal contribution of the ith ancestor (genetic contribution not yet explained by other ancestor chosen before), and a is the number of influential ancestors (i.e. ancestors with non-zero marginal contributions). Additionally, based on the marginal contributions, the number of ancestors explaining 50% of the genes in the population and the contribution of the main ancestor were obtained.

The fg was calculated as the inverse of twice the average coancestry of individuals in the predefined reference population. This parameter accounts for all causes of loss of genetic diversity in the population (Caballero and Toro et al., 2000).

The effective number of non-founders (fn) accounts only for the genetic drift effects in non-founder generations (Caballero and Toro et al., 2000):

Thus, the fg has two components, fe that accounts for the contributions of founders to the actual population, and fn which accounts for the contribution of non-founders, and is accumulated over generations. Based on these parameters, the genetic diversity (GD) in the reference population relative to the base population was obtained as (Lacy, 1989, 1995):

Therefore, 1 - GD represents the loss of genetic diversity in the population since the founder generation. The loss of genetic diversity because of unequal contributions of founders was estimated as 1 - GD*, where (Caballero and Toro et al., 2000):

Finally, the difference between GD and GD* accounts for the loss of genetic diversity due to genetic drift accumulated over the non-founder generations (Caballero and Toro, 2000), and it can be expressed as the inverse of twice the number of non-founders (Lacy, 1995).

Results

Demographic parametersTable 6

For the UP, average GI was similar from RP1 to RP3 (6.51 to 6.54 yr); however, for the SP these values tended to decrease from 7.16 in RP1 to 6.57 yr in RP3 (Table 1). Among parent-offspring pathways, the GI for path sire-son were larger and tended to decrease from RP1 to RP3 in both the UP and the SP, and were half to one year longer for the SP than the UP.

Table 1 Generation intervals (yr) for the upgraded and straightbred Romosinuano and their respective reference populations (RP1, RP2, and RP3). 

Generation intervals = average age of the parents when their replacements are born; RP1 = 1998-2003; RP2 = 2004-2009; RP3 = 2010-2017.

Table 2 Number of maximum (MG), complete (CG) and equivalent (EG) generations traced for the upgraded and straightbred Romosinuano and their respective reference populations (RP1, RP2, and RP3). 

RP1 = 1998-2003; RP2 = 2004-2009; RP3 = 2010-2017.

The level of completeness of the UP and the SP pedigrees improved from RP1 to RP3 more than 69, 43 and 63% for MG, CG and EG, respectively, after herd-book recording started (Table 2).

The EG reached 3.43 and 3.34 yr for the UP and the SP in RP3. The percentage of known ancestors improved during the recent 20 yr and for the UP averaged 76, 73 and 66% for the generations of parents, grandparents and great-grandparents, respectively, and those values for the SP were 93, 85 and 65% (Table 3); thus, the improvement was higher for the SP.

Table 3 Percentage of known ancestors by generation for the upgraded and straightbred Romosinuano and their respective reference populations (RP1, RP2, and RP3). 

RP1 = 1998-2003; RP2 = 2004-2009; RP3 = 2010-2017.

Considering various generations, there was a steady improvement of the PCI during the recent 20 yr (Figure 1). Including three generations back, the PCI of RP3 were similar, 65.3 and 65.2% for the UP and the SP, but it increased to 75.4 and 93.1%, respectively, when only one generation was considered.

Figure 1 Pedigree completeness index for the upgraded (UP) and straightbred (SP) Romosinuano and their respective reference populations (RP1, 1998-2003; RP2, 2004-2009; and RP3, 2010-2017) for generations 1 to 6. 

Mostly multiplier herds were found in this Romosinuano population (67 and 77% in the UP and the SP, respectively). Commercial herds in the UP were about twice as many as those in the SP, while no nucleus was found (Table 4).

Table 4 Classification of herds according to the origin and use of sires for the upgraded (UP) and straightbred (SP) Romosinuano populations. 

UFS = use foreign sires; UOS = use own sires; SS= sell sires; NoH = number of herds; PFS = percent of foreign sires.

As previously indicated, development of Mexican Romosinuano population was based mainly on germplasm importation from a herd in Turrialba, Costa Rica (herd 75), which is depicted in Figure 2. Regarding the UP, some herds (e.g. herd 12) included a large number of dams in upgrading to Romosinuano. In the SP, the genetic representation of herd 75 is declining, whereas that of Mexican herds are increasing (e.g. herd 40; 9.7 to 23.5% from RP1 to RP3, respectively).

Figure 2 Main founder herds genetically represented (%) in RP1 (1998-2003), RP2 (2004-2009), and RP3 (2010-2017) in upgraded (UP) and straightbred (SP) Romosinuano populations. 

Generally, all inbreeding statistics increased from RP1 to RP3 (Table 5). Considering RP3, all inbreeding statistics were higher in the SP than the UP, except for the percentage of aver age inbreeding in inbred animals.

Table 5 Inbreeding statistics (%) for upgraded and straightbred Romosinuano and their respective reference population (RP1, RP2, and RP3). 

RP1 = 1998-2003; RP2 = 2004-2009; RP3 = 2010-2017; IA = inbred animals; AR = average relatedness; IP = average inbreeding in the population; IF = average inbreeding in inbred animals.

Table 6 Average inbreeding coefficients (F), inbred animals and relatedness coefficients (AR) for upgraded and straightbred Romosinuano populations, by complete generations. 

FAR

Even when both Ne and Nec were larger for the UP than the SP, for Ne there was a reduction of these sizes in RP2 and a subsequent increase to reach sizes of 79 and 48 in RP3 (Table 7).

Estimates of Nec were smaller than Ne and increased over the recent 20 yr for the UP and the SP, reaching in RP3 57 and 45, respectively

Table 7 Realized effective population size (Ne) and effective population size from an increase in coancestry (Nec) for the upgraded and straightbred Romosinuano and their respective reference populations (RP1, RP2, and RP3). 

RP1 = 1998-2003; RP2 = 2004-2009; RP3 = 2010-2017.

The parameters describing the probabilities of gene origin for the UP and the SP and their respective reference populations (RP1, RP2, and RP3) are shown in Table 8. The f in the UP were larger than in the SP and increased from RP1 (273) to RP3 (827). Although fe remained approximately constant across reference populations the ratio fe/f for the UP decreased from RP1 (22.0%) to RP3 (7.4%) whereas those ratios for the SP were larger and varied from 20.3 to 32.2%. During the last 20 years, fa in the SP decreased 16.7% from RP1 to RP3 and fg had a reduction of 17.6 and 30.1% for the UP and the SP, respectively. The fa represented about 50% of the fe. The fg represented 27 and 33% of fe for the UP and the SP, respectively. From RP1 to RP3 there was a reduction of the proportion of founder genomes that have been retained in the reference populations (fg/f), from 7.3 to 2.0% in the UP and from 10.6 to 7.4% in the SP.

Table 8 Probability of gene origin parameters for the upgraded and straightbred Romosinuano and their respective reference populations (RP1, RP2, and RP3). 

RP1 = 1998-2003; RP2 = 2004-2009; RP3 = 2010-2017; f = number of founders; fe = effective number of founders; fa = effective number of ancestors; fg = equivalent number of founder genomes; Cmax = contribution of the main ancestor; N50= number of ancestors explaining 50% of the gene pool.

Total losses of genetic variability from RP1 to RP3 increased from 2.5 to 3.0% in the UP and from 2.6 to 3.7% in the SP (Figure 3).

More than 66% of the losses on genetic diversity were accounted for genetic drift, except for RP3 in the SP (44%).

Figure 3 Losses (%) of genetic diversity due to genetic drift (GD*-GD), unequal founder contribution (1-GD*), and both reasons (1-GD) since the founder generation in upgraded (UP) and straightbred (SP) Romosinuano reference populations (RP1, 1998-2003; RP2, 2004-2009; and RP3, 2010-2017). 

The main ancestor contributed 8.6 to 11.0, and 9.8 to 10.3% of the genes in the UP and the SP, respectively. The numbers of ancestors that explained half of the genetic variability of the populations ranged from 10 to 12 in the UP, and 9 to 11 in the SP, with a decreasing trend from RP1 to RP3.

Discussion

The average GI observed in the present study for the different reference populations (6.41 to 6.54 yr for UP and 6.47 to 7.16 yr for SP; Table 1) were larger than those of Colombian creole breeds such as Costeño con Cuernos (CCC, 5.4 yr), Blanco Orejinegro (BON, 4.7 yr) and Romosinuano (ROM, 5.7 yr), but similar to Sanmartinero (SM, 6.8 yr) (Martínez et al., 2008), and to Brazilian Lageana Crioula cattle breed (6.41 yr; Pezzini et al., 2018). Gutiérrez et al. (2003) estimated an interval of GI in eight autochthonous Spanish cattle breeds (Alistana, Asturiana de la Montaña, Asturiana de los Valles, Avileña-Negra Ibérica, Bruna dels Pirineus, Morucha, Pirenaica and Sayaguesa) of 3.7 to 6.1 yr. However, compared to the present study, larger GI were reported in Kumamoto (9.4 yr) and Kouchi (10.4 yr) Japanese Brown cattle breeds (Honda et al., 2006). Differences of GI among breeds may be the result of variations in calving age, fluctuating environmental conditions throught the years in the tropic (Rosendo et al., 2018), number of years that sires and dams are used, and breeding objectives; for example, large GI are convenient in conservation breeding programs.

Regarding the integrity of the pedigree, the EG of the present study increased during the recent 20 yr and reached in RP3, 3.43 for UP and 3.34 for SP (Table 2). Those estimates were smaller than those for Romosinuano (4.8), SM (3.8) and CCC (3.7) Colombian creole breeds (Martínez et al., 2008) but higher than the interval 0.81 to 2.97 in eight Spanish autochthonous breeds (Gutiérrez et al., 2003), the EG in BON (3.1; Martínez et al., 2008), and the estimate in Lageana Crioula cattle breed (0.77, Pezzini et al., 2018). The PCI is an important measure of the quality of inbreeding estimates from the pedigree since PCI represents a harmonic mean of parental genetic contributions and equals zero if any parent is unknown regardless of how deep and complete is the pedigree for the other parent. Similarly, inbreeding can be estimated only if information for both parents is known (Battagin et al., 2010; Onogi et al., 2017). Estimates of PCI in RP3 from the present study (65.3 for UP and 65.2% for SP, and considering three generations back, and 75.4 and 93.1% considering one generation back; Figure 1) were lower than the PCI (90 and 69% in 2009 considering 2 and 5 generations, respectively) observed by Battagin et al. (2010) in Burlina native Italian cattle, and those reported by Santana et al. (2014) for Brazilian Gir cattle (99.8, 96.3 and 92.8% for generations 1, 2 and 3, respectively). The observed differences of EG and PCI among studies could be the result of variations in the number of generations traced and the quality of pedigree recording.

The three different ways to evaluate the degree of pedigree completeness (number of MG, CG and EG; percentage of known ancestors up to three generations; and PCI) indicated that for the UP and the SP the integrity of the pedigree was improved during the latest 20 yr. It is clear that an important progress has occurred after AMCROLET was established; however, these alternative measures suggest that greater efforts are needed in the recording of genealogical data to improve the accuracy of inbreeding parameters in this population (Boichard et al., 1997).

Given that Romosinuano breed was developed in Colombia and later introduced into Mexico, no nucleus herds were found. The high proportion of multiplier herds found in the present study indicates that the population is growing. Similar results were reported by Silva et al. (2016). Since the demand for Romosinuano germplasm in Mexico is increasing, commercial herds in the UP were about twice as many as those in the SP. Mexican herds (like herd 40, Figure 3) are having an increasing genetic representation in the population. Thus, it is beneficial if breeders share germplasm from various herds to avoid short-term significant increases in inbreeding, with the consequent reduction in animal performance (Burrow, 1993; Santana et al., 2012). Germplasm used in the UP matings mainly comes from the SP; thus, it is important to compare results of the SP and RP3 with other studies. The percentage of inbred animals in SP (65.8%) was higher than the range 5.8 to 30.8% in Colombian creole breeds (Martínez et al., 2008), and the interval 2 to 32% in Spanish autochthonous breeds (Gutiérrez et al., 2003), but lower than in Alentejana Portuguese (nearly 80%; Carolino and Gama, 2007), and in Burlina Italian cattle (>80%; Battagin et al., 2010). Possible causes of these differences among studies could be the pedigree completeness, monitoring of matings to avoid inbreeding, and gene flow among herds. Even though averages of F for the Romosinuano breed in Mexico have been increasing during the last 20 yr, they are still low (2.5 and 3.9% for all and inbred animals, respectively; Table 5). Lower values (<1.41% for the reference population, and <3.21% for inbred animals) were estimated for various breeds by Martínez et al. (2008), except for Romosinuano that had 3.1 and 4.0%, respectively. These differences may be due mainly to variations in effective population sizes, integrity level of the pedigree, and efficacy to prevent inbreeding. The AR was high (5.9%) as compared to the interval 2.0 to 3.9% reported by Martínez et al. (2008) in four Colombian breeds, and the range 0.10 to 1.70% by Gutiérrez et al. (2003) in eight Spanish breeds. The size of AR is important to predict the F of following generations since, under random mating, F of the progeny will be half of the AR between parents. For RP3, in the SP the expected F in the following generation would be still low at about 3%. However, these estimates may be biased downward because of incomplete pedigree data since pedigree recording in the Romosinuano Herd-Book was formally carried out after AMCROLET was created (1998).

The average Nec and Ne reached 45 and 48 in RP3 SP, respectively (Table 7). The Ne estimates from the present study are within the range (21 to 123) published by Gutiérrez et al. (2003). Larger values of Ne were reported in Romosinuano (66), BON (260), and CCC (143) Colombian breeds, except for SM (27) (Martínez et al., 2008), in Tabapuã (411; Bernardes et al., 2016), and in Brazilian Gir cattle (94; Santana et al., 2014). Estimates of Nec were smaller than Ne and increased over the last 20 yr. Smaller values of Nec than Ne were also reported by Santana et al. (2012) for Marchigiana (98 vs 140) and Bonsmara (55 vs 325) cattle breeds. However, Santana et al. (2014) found higher Nec than Ne in Brazilian Gir (166 vs 94). As indicated by Cervantes et al. (2011), the differences between Ne and Nec are mainly due to population substructures (caused by mating policies, breeding goals or geographical distances) that affect the increase in F, while the increase in coancestry practically does not change. For RP3 in the SP, the size of Nec (45) was lower than the limit (>50) recommended for conservation of genetic resources to have rates of inbreeding per generation smaller than 1.0% (FAO, 1998).

In the present study, there was an uncommon relationship between Ne and F, since there was a trend on both parameters to increase over time (Tables 5 and 7). There is a direct relationship between Ne and . Additionally, we expect that as Ne increases, F decreases. However, sometimes this does not happen, since the depends not only on the F, but also on the number of equivalent generations (). Other authors have reported this situation; for example, Rosendo et al. (2018) published Ne of 68.1 for the total and 64.6 for the selection nucleus Tropical Milking Criollo cattle, but corresponding F were 1.07 and 2.14.

Boichard et al. (1997) indicated that F trend, although commolnly used to quantify the rate of genetic drift, has some limitations; particularly when the pedigree completeness is low. The analysis of probabilities of gene origin is suggested as a complementary approach. In the present study, the fe/f ratio indicated that part of the initial genetic variability was lost due to unbalanced gene contributions of founders to the reference populations (Boichard et al., 1997). Additionaly, the fa represented about 50% of the fe, which results from bottlenecks during the development of the population. Estimates of fg for RP3 in the UP and the SP (Table 8) indicate that the amount of genetic diversity can be generated by 17 and 14 unrelated founders, respectively (Lacy, 1989; Ballou and Lacy, 1995).

Both, genetic drift and unbalanced founder contributions of genes to the reference populations were responsible for the reduction of genetic diversity, and the former accounted for >66% of the total losses, except for the SP most recent reference population (44%). Honda et al. (2006) also reported that the decay in genetic diversity in Kumamoto and Kouchi sub-breeds of Japanese Brown cattle was due mainly to accumulated genetic drift (>78%, except for Kumamoto, 16%, in 2000 reference population; estimated by the authors of the present study).

The numbers of ancestors that explained 50% of the genetic variability of Romosinuano in the present study (9 to 12) were similar to Lageana Crioula breed (10; Pezzini et al., 2018), Romosinuano (8) and CCC (12), but smaller than SM (22) and BON (22) Colombian breeds (Martínez et al., 2008), Marchigiana (13) and Bonsmara (41) (Santana et al., 2012), and Gir (38; Santana et al., 2014).

The estimates of demographic and genetic parameters suggest that the Romosinuano population in Mexico is on the borderline of risk status. One of the purposes of conservation is to maintain, as much as possible, the genetic diversity present in the founder population. Among the alternatives to maintain high levels of genetic diversity in the population, the most efficient has been to minimize the average pairwise coancestries every generation and by definition to maximize fg (Ballou and Lacy, 1995; Lacy 1995; Caballero and Toro, 2000). Therefore, it is important that Romosinuano breeders contribute to reduce losses of genetic variability through matings designed to minimize average pairwise coancestries and a more balanced use of sires.

In conclusion, reduction of the genetic diversity has been occurring since the Romosinuano breed association was established in Mexico, mainly due to random losses of genes. Even though no nucleus herds were found in this study, there was a tendency to concentrate the use of germplasm from a reduced number of herds; therefore, it is important to diversify the sources of genetic material for future matings. Monitoring genetic parameters and defining matings aimed to minimize pairwise coancestries are suggested to increase the effective population size and thus, to reduce the inbreeding rate per generation. Efforts are needed to improve the level of completeness of the pedigree and the accuracy of estimates of demographic and genetic parameters.

Acknowledgements

The authors want to thank the people at Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical (AMCROLET) for allowing us to use their database.

References

AMCROLET. Resumen de la evaluación genética para sementales Romosinuano 2016. Chapingo, México. Asociación Mexicana de Criadores de Ganado Romosinuano y Lechero Tropical, AC; 2017. [ Links ]

Ballou JD, Lacy RC. Identifying genetically important individuals for management of genetic diversity in pedigreed populations. In: Ballou JD, Gilpin M, Foose TJ, editors. Population Management for Survival & Recovery. Analytical Methods and Strategies in Small Population Conservation. New York: Columbia University Press; 1995. p.76-111. [ Links ]

Battagin M, Penasa M, Pretto D, Cassandro M. Pedigree analysis of Burlina cattle population. Acta Agrar Kaposv 2010; 14:161-5. [ Links ]

Bernardes PA, Grossi DA, Savegnago RP, Buzanskas ME, Ramos SB, Romanzini EP, Guidolin DGF, Bezerra LAF, Lôbo RB, Munari DP. Population structure of Tabapuã beef cattle using pedigree analysis. Livest Sci 2016; 187:96-101. https://doi.org/10.1016/j.livsci.2016.03.002Links ]

Boichard D, Maignel L, Verrier E. The value of using probabilities of gene origin to measure genetic variability in a population. Genet Sel Evol 1997; 29:5-23. https://doi.org/10.1186/1297-9686-29-1-5Links ]

Burrow HM. The effects of inbreeding in beef cattle. Anim Breed Abstr 1993; 61:737-51. [ Links ]

Caballero A, Toro M. Interrelations between effective population size and other pedigree tools for the management of conserved populations. Genet Res 2000; 75:331-43. https://doi.org/10.1017/S0016672399004449Links ]

Carolino N, Gama LT. Indicators of genetic erosion in an endangered population: The Alentejana cattle breed in Portugal. J Anim Sci 2007; 86:47-56. https://doi.org/10.2527/jas.2007-0148Links ]

Carroll JA, Burdick NC, Reutera RR, Chase CCJr, Spiers DE, Arthington JD, Coleman SW. Differential acute phase immune responses by Angus and Romosinuano steers following an endotoxin challenge. Domest Anim Endocrin 2011; 41:163-73. https://doi.org/10.1016/j.domaniend.2011.06.002Links ]

Carroll JA, Burdick NC, Chase CCJr, Coleman SW, Spiers DE. Influence of environmental temperature on the physiological, endocrine, and immune responses in livestock exposed to a provocative immune challenge. Domest Anim Endocrin 2012; 43:146-53. https://doi.org/10.1016/j.domaniend.2011.12.008Links ]

Cervantes I, Goyache F, Molina A, Valera M, Gutiérrez JP. Estimation of effective population size from the rate of coancestry in pedigreed populations. J Anim Breed Genet 2011; 128:56-63. https://doi.org/10.1111/j.1439-0388.2010.00881.xLinks ]

De Alba J. El Libro de los Bovinos Criollos de América. Biblioteca Básica de Agricultura.. Montecillo (Estado de México): Ed. Colegio de Postgraduados. 2011. [ Links ]

FAO (Food and Agriculture Organization for the United Nations). Secondary guidelines for development of national farm animal genetic resources management plans: Management of small populations at risk. Rome, Italy; 1998. [ Links ]

Gutiérrez JP, Altarriba J, Díaz C, Quintanilla R, Cañón J, Piedrafita J. Pedigree analysis of eight Spanish beef cattle breeds. Genet Sel Evol 2003; 35:43-63. https://doi.org/10.1186/1297-9686-35-1-43Links ]

Gutiérrez JP, Goyache F. A note on Endog: a computer program for analysing pedigree information. J Anim Breed Genet 2005; 122:172-76. https://doi.org/10.1111/j.1439-0388.2005.00512.xLinks ]

Gutiérrez JP, Cervantes I, Goyache F. Improving the estimation of realized effective population sizes in farm animals. J Anim Breed Genet 2009; 126:327-32. https://doi.org/10.1111/j.1439-0388.2009.00810.xLinks ]

Hammerly SC, Morrow ME, Johnson JA. A comparison of pedigree- and DNA-based measures for identifying inbreeding depression in the critically endangered Attwater’s Prairie-chicken. Mol Ecol 2013; 22:5313-28. https://doi.org/10.1111/mec.12482Links ]

Honda T, Fujii T, Nomura T, Mukai F. Evaluation of genetic diversity in Japanese Brown cattle population by pedigree analysis. J Anim Breed Genet 2006; 123:172-9. https://doi.org/10.1111/j.1439-0388.2006.00586.xLinks ]

Jamieson IG, Allendorf FW. How does the 50/500 rule apply to MVPs?. Trends Ecol Evol 2012; 27:578-84. https://doi.org/10.1016/j.tree.2012.07.001Links ]

Lacy RC. Analysis of founder representation in pedigrees: founder equivalents and founder genome equivalents. Zoo Biol 1989; 8:111-23. https://doi.org/10.1002/zoo.1430080203Links ]

Lacy RC. Clarification of genetic terms and their use in the management of captive populations. Zoo Biol 1995; 14:565-78. https://doi.org/10.1002/zoo.1430140609Links ]

MacCluer JW, Boyce AJ, Dyke B, Weitkamp LR, Pfenning DW, Parsons CJ. Inbreeding and pedigree structure in Standardbred horses. J Hered 1983; 74:394-9. https://doi.org/10.1093/oxfordjournals.jhered.a109824Links ]

Maignel L, Boichard D, Verrier E. Genetic variability of French dairy breeds estimated from pedigree information. Interbull Bull 1996; 14:49-54. [ Links ]

Martínez RA, García D, Gallego JL, Onofre G, Pérez J, Cañón J. Genetic variability in Colombian Creole cattle populations estimated by pedigree information. J Anim Sci 2008; 86:545-52. https://doi.org/10.2527/jas.2007-0175Links ]

Meuwissen THE, Luo Z. Computing inbreeding coefficients in large populations. Genet Sel Evol 1992; 24:305-13. doi:10.1186/1297-9686-24-4-305 [ Links ]

Onogi A, Shirai K, Amano T. Investigation of genetic diversity and inbreeding in a Japanese native horse breed for suggestions on its conservation. Anim Sci J 2017; 88:1902-10. https://doi.org/10.1111/asj.12867Links ]

Pezzini T, Mariante AS, Martins E, Paiva S, Seixas L, Costas Jr JBG, Rolo J, McManus C. Population structure of Brazilian Crioula Lageana cattle (Bos taurus) breed. Rev Colomb Cienc Pecu 2018; 31:93-102. [ Links ]

Pinheiro M, Kjöllerström HJ, Oom MM. Genetic diversity and demographic structure of the endangered Sorraia horse breed assessed through pedigree analysis. Livest Sci 2013; 152:1-10. https://doi.org/10.1016/j.livsci.2012.11.017Links ]

Riley DG, Chase CCJr, Coleman SW, Olson TA. Evaluation of birth and weaning traits of Romosinuano calves as purebreds and crosses with Brahman and Angus. J Anim Sci 2007; 85:289-98. https://doi.org/10.2527/jas.2006-416Links ]

Riley DG, Burke JM, Chase CCJr, Coleman SW. Heterosis and direct effects for Charolais-sired calf weight and growth, cow weight and weight change, and rations of cow and calf weights and weight changes across warm season lactation in Romosinuano, Angus, and F1 cows in Arkansas. J Anim Sci 2016; 94:1-12. https://doi.org/10.2527/jas.2015-9484Links ]

Rosendo PA, Palacios JAL, Rosales MF, Torres HG, Ramírez VR, Becerril PCM. Genetic variability of Tropical Milking Criollo cattle of Mexico estimated from genealogical information. Rev Colomb Cienc Pecu 2018; 31:196-203. http://dx.doi.org/10.17533/udea.rccp.v31n3a04Links ]

Santana MLJr, Oliveira PS, Eler JP, Gutiérrez JP, Ferraz JBS. Pedigree analysis and inbreeding depression on growth traits in Brazilian Marchigiana and Bonsmara breeds. J Anim Sci 2012; 90:99-108. https://doi.org/10.2527/jas.2011-4079Links ]

Santana MLJr, Pereira RJ, Bignardi AB, El Faro L, Tonhati H, Albuquerque LG. History, structure, and genetic diversity of Brazilian Gir cattle. Livest Sci 2014; 163:26-33. https://doi.org/10.1016/j.livsci.2014.02.007Links ]

Santana MLJr, Pereira RJ, Bignardi AB, Ayres DR, Menezes GRO, Silva LOC, Leroy G, Machado CHC, Josahkian LA, Albuquerque LG. Structure and genetic diversity of Brazilian Zebu cattle breeds assessed by pedigree analysis. Livest Sci 2016; 187:6-15. https://doi.org/10.1016/j.livsci.2016.02.002Links ]

Scharf B, Carroll JA, Riley DG, Chase CCJr, Coleman SW, Keisler DH, Weaber RL, Spiers DE. Evaluation of physiological and blood serum differences in heat-tolerant (Romosinuano) and heat suceptible (Angus) Bos taurus cattle during controlled heat challenge. J Anim Sci 2010; 88:2321-36. https://doi.org/10.2527/jas.2009-2551Links ]

SAGARPA (Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación). Comunicado de Prensa BCS/226/2017. 2017; access date: 27/07/2017 2017; access date: 27/07/2017 http://sagarpa.gob.mx/Delegaciones/bajacaliforniasur/boletines/2017/junio/Documents/2017BS226.PDFLinks ]

Sheikhlou M, Abbasi MA. Genetic diversity of Iranian Lori-Bakhtiari sheep assessed by pedigree analysis. Small Ruminant Res 2016; 141:99-105. https://doi.org/10.1016/j.smallrumres.2016.07.009Links ]

Silva MHMA da, Malhado CHM, Costa JLJr, Cobuci JA, Costa CN, Carneiro PLS. Population genetic structure in the Holstein breed in Brazil. Trop Anim Health Prod 2016; 48:331-6. DOI: 10.1007/s11250-015-0956-7 [ Links ]

Toro MA, Villanueva B, Fernández J. Genomics applied to management strategies in conservation programmes. Livest Sci 2014; 166:48-53. https://doi.org/10.1016/j.livsci.2014.04.020Links ]

Vassallo JM, Diaz C, Garcia-Medina JR. A note on the population structure of the Avileña breed of cattle in Spain. Livest Prod Sci 1986; 15:285-8. https://doi.org/10.1016/0301-6226(86)90035-7Links ]

* To cite this article: Nuñez-Domínguez R, Martínez-Rocha RE, Hidalgo-Moreno JÁ, Ramírez-Valverde R, García-Muñiz JG. Evaluation of the Romosinuano cattle population structure in Mexico using pedigree analysis. Rev Colomb Cienc Pecu 2020; 33(1): 44-59 DOI: https://doi.org/10.17533/udea.rccp.v32n4a05

Received: February 14, 2018; Accepted: March 18, 2019

Conflicts of interest

The authors declare they have no conflicts of interest with regard to the work presented in this report.

Author contributions

This study was designed and directed by Rafael Núñez Domínguez, Ph.D. and Rodolfo Ramírez Valverde, Ph.D. Part of the study was presented as a Professional Thesis by Ricardo E. Martínez Rocha, M.C. All co-authors were closely involved in the collection, analysis and interpretation of data as well as in the writing and decision to submit the work for publication.

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License