Pelibuey is the most important hair sheep breed in Mexico. It entered the country between 1930 and 1940 from Cuba^1,2,3, and since then, no new genetic material has entered the country. Initially, Pelibuey sheep were distributed in tropical regions. Nowadays, it is spread throughout the country⁴. Pelibuey sheep are not highly productive; their biological importance lies in their ability to adapt to different environments and climates⁵ and reproduce throughout the year⁶.

• ¹
  Genética del borrego Tabasco o Peligüey

  Téc Pecu, 1975
  Genética del borrego Tabasco o Peligüey
• ²
  El ovino Pelibuey en el trópico mexicano.

  Biodiversidad ovina iberoamericana. Caracterización y uso sustentable, 2009
  2. Montalvo MP, Romualdo MJG, Sierra VA, Ortiz OJ, Hernández ZJ, Medrano HA. El ovino Pelibuey en el trópico mexicano. En: Delgado BJV, Nogales BS editores. Biodiversidad ovina iberoamericana. Caracterización y uso sustentable. España, Universidad de Córdoba. Servicio de Publicaciones 2009:363-375.
• ³
  Origen, historia y situación actual de la oveja Pelibuey en México

  Trop Subtrop Agroecosys, 2017
  Origen, historia y situación actual de la oveja Pelibuey en México
• ⁴
  Diagnóstico actual de los ovinos en México, 2007
  4. Arteaga CJD. Diagnóstico actual de los ovinos en México. En: Memorias del 8 Congreso Mundial de la Lana y el Cordero. Asociación Mexicana de Criadores de Ovinos. Querétaro, Qro.2007:9-12.
• ⁵
  A field study of reproductive performance and productivity of Pelibuey ewes in southeaster Mexico

  Trop Anim Health Prod, 2013
  A field study of reproductive performance and productivity of Pelibuey ewes in southeaster Mexico
• ⁶
  Actividad reproductiva de la oveja Pelibuey durante la época de anestro: influencia de la presencia del macho

  Rev Científica FCV-LUZ, 2006
  Actividad reproductiva de la oveja Pelibuey durante la época de anestro: influencia de la presencia del macho

The wide range of environments in which Pelibuey sheep have been bred has led to adaptive responses, which are part of their gene pool. However, in recent years, this breed has been subjected to indiscriminate crosses with breeds specialized in meat production to increase its productivity^3,7; this has put at risk the original genetic diversity of the breed. Additionally, artificial insemination and the intense flow of genetic material between flocks have exacerbated the problem⁸.

• ³
  Origen, historia y situación actual de la oveja Pelibuey en México

  Trop Subtrop Agroecosys, 2017
  Origen, historia y situación actual de la oveja Pelibuey en México
• ⁷
  Morphological characterization of Pelibuey sheep in Colima, México

  Trop Anim Health Prod, 2013
  Morphological characterization of Pelibuey sheep in Colima, México
• ⁸
  Are cattle, sheep and goats endangered species?

  Mol Ecol, 2008
  8. Taberlet P, Valentini A, Rezaei HR, Naderi S, Pompanon F, Negrini R, Ajmone-Marsan P. Are cattle, sheep and goats endangered species? Mol Ecol 2008;17(1):275-284.

An initial step for the conservation of animal genetic resources is breed characterization. The Food and Agriculture Organization of the United Nations (FAO) has proposed that the initial characterization should include a phenotypic description of the breed and its subsequent genetic characterization using molecular markers⁹. Molecular markers, such as microsatellites, have been used in population genetics to characterize several breeds⁹. More profound knowledge about the genetic diversity and variability and the population structure of different flocks in Mexico will allow us to determine the degree of risk of the Pelibuey sheep and, consequently, suggest viable strategies for its conservation. Therefore, this study aimed to characterize different Pelibuey subpopulations in Mexico using microsatellite markers.

• ⁹
  The second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture, 2015
  The second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture
• ⁹
  The second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture, 2015
  The second Report on the State of the World’s Animal Genetic Resources for Food and Agriculture

Blood samples were collected from 119 Pelibuey sheep from 23 domestic flocks and one Cuban flock. Domestic flocks were sampled in the agroecological zones that correspond to the tropical wet, tropical dry, and central mountainous regions. The geographical location of the domestic flocks is shown in Figure 1. Of the 24 flocks, 11 belong to universities or research institutes; the remaining 13 belong to private producers. All the flocks in this study were handled according to the Institutional Animal Care and Use Committee (SICUAE, as per its acronym in Spanish) of the School of Veterinary Medicine and Zootechnics of the UNAM. The flocks that belong to the Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP) in Yucatan (INIFAP Mocochá=IN-MOC) and Puebla (INIFAP Las Margaritas= IN-MAR) were considered as the domestic reference subpopulations because they were the first to settle in Mexico and their breeding practices have targeted the conservation of their breed. Furthermore, both have served as the basis for the formation of other domestic flocks. The inclusion criteria were the following: an external appearance characteristic of the Pelibuey sheep, female specimens, clinically healthy, and not related to each other.

Following FAO recommendations, the analysis included five individuals from each flock¹⁰. We only included four samples of the Rancho Jalapa and Rancho El Carrizal flocks because of the DNA degradation of some samples.

• ¹⁰
  Molecular genetic characterization of animal genetic resources, 2011
  Molecular genetic characterization of animal genetic resources

Blood samples (6 ml) were obtained by jugular venipuncture and collected in Vacutainer® tubes with EDTA as an anticoagulant; samples were labeled at the time of sampling. Samples were stored at -80 °C until further processing and analysis.

DNA isolation was carried out based on a simple and inexpensive protocol for extracting DNA from poultry blood samples¹¹. This protocol was subsequently adapted for sheep blood¹².

• ¹¹
  An inexpensive, simple protocol for DNA isolation from blood for high-throughput genotyping by polymerase chain reaction or restriction endonuclease digestion

  Poult Sci, 2007
  An inexpensive, simple protocol for DNA isolation from blood for high-throughput genotyping by polymerase chain reaction or restriction endonuclease digestion
• ¹²
  Identificación de ovejas Pelibuey con actividad reproductiva continua mediante polimorfismos de ADN, 2015
  Identificación de ovejas Pelibuey con actividad reproductiva continua mediante polimorfismos de ADN

The yield and purity of the extracted DNA were determined by spectrophotometry (NanoDrop Thermo Scientific, Wilmington, DE). These results were corroborated by a 0.8% agarose gel electrophoresis using 10 µL of the extracted DNA. The extracted DNA was stored at -80 °C until further used.

Nine microsatellite markers were amplified through PCR following the International Society for Animal Genetics (ISAG) and the FAO^10,13 recommendations. The PCR reactions were carried out in a final volume of 25 µL and consisted of 100 ng of genomic DNA (2 µL), 0.2 µM of each primer (Forward and Reverse; 0.5 µL), 200 µM of each dNTP (2.5 µL), 2.5 mM of MgCl²⁺ (0.5 µL), 1.25 U of DNA Taq polymerase (0.25 µL), 10 X PCR buffer (2.5 µL), and 16.25 µL of sterile water.

• ¹⁰
  Molecular genetic characterization of animal genetic resources, 2011
  Molecular genetic characterization of animal genetic resources
• ¹³
  Panel of markers for parentage verification tested 2001/02 ISAG comparison test, 2002
  Panel of markers for parentage verification tested 2001/02 ISAG comparison test

PCR reactions were carried out in a thermocycler (Axygen Scientific Inc.). The PCR protocol was the following: an initial denaturation step at 95 °C for 5 min; followed by 35 cycles of denaturation at 94 °C for 45 s, annealing for 1 min at variable temperature, and extension at 72 °C for 1 min; with a final extension at 72 °C for 10 min. Amplifications were stored at 4 °C until further analysis.

Amplifications were subjected to a non-denaturing polyacrylamide (12 %) gel electrophoresis. Electrophoresis was carried out with 0.5% TBE as a running buffer. The size of the DNA fragments was determined using a 25 bp molecular weight ladder (Invitrogen Life Technologies, Carlsbad, USA). Gels were stained with 0.5 µg/ml of ethidium bromide and visualized under ultraviolet light (Kodak Gel Logic 2200 Imaging System). The gel images were processed using the MyImageAnalysisSystem^TM (Fisher) software.

The number of alleles (NA), effective number of alleles (NE), mean number of alleles (MNA), Shannon index (I), observed heterozygosity (HO), and expected heterozygosity (HE) for each marker were estimated using the POPGENE v.1.31¹⁴ and FSTAT¹⁵ software. The polymorphic information content (PIC) was calculated with the CERVUS v 3.0.7¹⁶ software. For each marker, exact tests were carried out to determine deviations from Hardy-Weinberg equilibrium (HWE) through the Markov chain Monte Carlo simulation method using the GENEPOP v 4.7.0^17,18 software.

• ¹⁴
  Microsoft Windows based freeware for population genetic analysis, 1999
  Microsoft Windows based freeware for population genetic analysis
• ¹⁵
  FSTAT, a program to estimate and test gene diversities and fixation indices F-statistics, 2001
  FSTAT, a program to estimate and test gene diversities and fixation indices F-statistics
• ¹⁶
  Revising how the computer program CERVUS accommodates error increases success in paternity assignment

  Mol Ecol, 2007
  Revising how the computer program CERVUS accommodates error increases success in paternity assignment
• ¹⁷
  GENEPOP (version 1.2) population genetics software for exact test and ecumenicism

  J Hered, 1995
  GENEPOP (version 1.2) population genetics software for exact test and ecumenicism
• ¹⁸
  Genepop’007: a complete reimplementation of the GENEPOP software for Windows and Linux

  Mol Ecol Resour, 2008
  Genepop’007: a complete reimplementation of the GENEPOP software for Windows and Linux

F statistics measured genetic structure. The inbreeding index within a population (FIS), fixation index of the total population (FIT), and the genetic differentiation index between populations (FST) were calculated with the GENEPOP v 4.7.0 software, carrying out 1,000 permutations^17,18. Private alleles were analyzed using the GENALEX v 6.5 software^19,20. Shared alleles between subpopulations were analyzed by creating a binary matrix that recorded the presence (1) or absence (0) of each of the alleles in the subpopulations. Subsequently, the proportions of the subpopulations with the most and least common shared alleles were calculated.

• ¹⁷
  GENEPOP (version 1.2) population genetics software for exact test and ecumenicism

  J Hered, 1995
  GENEPOP (version 1.2) population genetics software for exact test and ecumenicism
• ¹⁸
  Genepop’007: a complete reimplementation of the GENEPOP software for Windows and Linux

  Mol Ecol Resour, 2008
  Genepop’007: a complete reimplementation of the GENEPOP software for Windows and Linux
• ¹⁹
  Population genetic software for teaching and research

  Mol Ecol Notes, 2006
  Population genetic software for teaching and research
• ²⁰
  , 2012
  20. Peakall R, Smouse PE. GenAlEx 6.5: genetics analysis in Excel. Population software for teaching and research-an update. Bioinformatics 2012;28(19):2537-2539.

The genetic flow between the subpopulations was measured through the number of migrants per generation (Nm) and calculated from the FST between pairs of subpopulations using the GENETIX v. 4.05 software²¹.

• ²¹
  GENETIX 4.05, logiciel sous WindowsTM pour la génetique des populations, 2004
  GENETIX 4.05, logiciel sous WindowsTM pour la génetique des populations

Isolation by distance between subpopulations was analyzed using a Mantel test²², which consisted of correlating the genetic distances (measured through Nei's unbiased minimum distances²³) and the geographic distances between the subpopulation pairs using the XLSTAT software. This test used a significance of 0.05 and 10,000 permutations.

• ²²
  The detection of disease clustering and a generalized regression approach

  Cancer Res, 1967
  The detection of disease clustering and a generalized regression approach
• ²³
  Estimation of average heterozygosity and genetic distance from a small number of individuals

  Genetics, 1978
  Estimation of average heterozygosity and genetic distance from a small number of individuals

The spatial relationships between the subpopulations were analyzed using a principal component analysis (PCA) implemented in the PCAGEN software²⁴. The population structure and degree of mixture were estimated using a Bayesian clustering model implemented in the STRUCTURE v 2.3²⁵ software. This analysis involved a mixture model with correlated allele frequencies. To choose the correct number of inferred clusters (K) needed to model the data, 2-24 clusters were inferred with 20 independent runs in each cluster. All runs used 100,000 interactions (burn-in), followed by 1,000,000 interactions (MCMC). The most probable number of K was calculated with the ΔK algorithm using the online software STRUCTURE HARVESTER²⁶. Finally, results were processed using the online software CLUMPAK²⁷ to interpret the obtained inferences.

• ²⁴
  PCAGEN. Principal component analysis of gene frequency data, 1999
  PCAGEN. Principal component analysis of gene frequency data
• ²⁵
  Inference of populations structure using multilocus genotype data

  Genetics, 2000
  Inference of populations structure using multilocus genotype data
• ²⁶
  STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method

  Conserv Genet Resour, 2012
  STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method
• ²⁷
  CLUMPAK: a program for identifying clustering modes and packaging populations structure inferences across K

  Mol Ecol Resour, 2015
  CLUMPAK: a program for identifying clustering modes and packaging populations structure inferences across K

Table 1 shows the global statistics obtained with the nine microsatellites. For the entire Pelibuey population, 99 alleles were detected in 119 animals genotyped with nine microsatellite primers. All loci were polymorphic in the analyzed subpopulations. The number of alleles per locus fluctuated from 9 (OarCP34, OarFCB304, and ILSTS5) to 14 (OarJMP58). The observed average number of alleles per locus was 11. The NE ranged from 5.31 (OarCP34) to 8.68 (OarJMP58).

The HO values ranged from 0.657 (SRCRSP9) to 0.812 (OarFCB304). The average HO for the entire population was 0.726. Moreover, HE fluctuated from 0.664 (ILSTS5) to 0.807 (OarJMP29), with an average of 0.732. The PIC for each marker ranged from 0.79 (OarCP34) to 0.87 (OarJMP58 and DYMS1), with a global average of 0.84. The I value was between 1.90 (OarCP34) and 2.35 (OarJMP58), with a global average of 2.1.

The microsatellite markers were tested for deviation from the HWE. Most loci were under HWE. However, OarFCB304, OarJMP29, and ILSTS5 showed a significant deviation from the HWE (P<0.05). The global values of FIS, FIT, and FST were 0.007, 0.158, and 0.151, respectively, with no significant differences (P>0.05). The highest FIS value (0.136) was observed in the OarJMP29 marker; the lowest value (-0.176) was observed in the ILSTS5 marker. The highest FIT value (0.255) was observed in the SRCRSP9 locus, the lowest (0.075) in the OarFCB304 locus. Furthermore, the highest FST value (0.228) was observed in the ILSTS5 locus, the lowest (0.056) in the OarJMP29 locus. The differences between subpopulations, evaluated by the multilocus values of FST, revealed that, for the most part, genetic variation corresponds to differences between individuals within the subpopulations (84.9 %); 15.1 % results from differences between subpopulations.

Table 2 shows the genetic diversity measures for each subpopulation. The subpopulation with the higher genetic diversity was UNPA (49 alleles); the subpopulation with the lowest genetic diversity was IN-MOC (28 alleles), with a global average of 36.6 alleles per subpopulation. The average MNA and NE were 4.08 and 3.25, respectively. The highest values were observed in the UNPA subpopulation (5.44 and 4.20, respectively); the IN-MOC subpopulation had the lowest values (3.11 and 2.46, respectively).

The average HO was 0.726. The lowest value was observed in the IN-MOC subpopulation (0.666), the highest in COL-TEX (0.800). The average HE was 0.731. The lowest value was observed in the IN-MOC subpopulation (0.627), the highest in UNPA (0.840).

The average FIS for each subpopulation, considering all loci, fluctuated from -0.171 (CEPIPSA) to 0.138 (Rancho Libertad). Nine of the 24 analyzed subpopulations (IN-MOC, ITC, UV, CEIEGT, IN-MAR, Rancho El Paraíso, CEPIPSA, UCOL, and Cuba) showed negative FIS values; the remaining 15 had positive values. Six subpopulations showed (Table 2) FIS values different from zero (P<0.05), with positive (Rancho Garrido, Rancho Libertad, and Rancho El Carrizal) and negative values (ITC, IN-MAR, and CEPIPSA).

The allele analysis (Table 3; end of manuscript) demonstrated the presence of nine private alleles (9.09% of the 99 alleles) distributed in seven subpopulations (29.1% of the 24 flocks). The Posta El Cuatro subpopulation had three private alleles; CEIEGT, COLPOS-CORD, Rancho San Alberto, Rancho Belbesah, Rancho Libertad, and Rancho Jalapa had one private allele. The loci that contributed to private alleles were OarFCB304 (2), OarJMP29 (2), SRCRSP5 (2), MAF33 (2), and SRCRSP9 (1). The frequencies of private alleles ranged from 0.083 in Posta El Cuatro (markers OarFCB304 and OarJMP29) to 0.400 in CEIEGT (marker SRCRSP9). The Cuban Pelibuey flock had no private alleles.

Although not all subpopulations shared alleles, it was identified two alleles ("110" and "150" of the OarCP34 and OarJMP29 loci, respectively) that were shared by 23 subpopulations (95.83 %). Table 4 (end of the manuscript) shows the analysis of the genetic flow (measured by Nm). The highest value of Nm was observed between the CEIEGT and CEPIPSA subpopulations. Of the 276 possible subpopulation comparisons, 222 (80.43 %) had a Nm higher than 1; 14 (5.07 %) had a Nm higher than 4.

The isolation by distance in the subpopulations was corroborated using Mantel's test (Figure 2). This analysis showed that the genetic differentiation between the Pelibuey sheep subpopulations (measured through Nei's unbiased minimum distances) is not related to their geographic distances (r= -0.062; P>0.05).

A PCA was carried out using the frequencies of the 99 alleles observed in the population. Figure 3 shows the global PCA. The first two principal components explained 30.21 % of the total variation. The first axis contributed 17.04 % of the variance and separated the flocks into two clusters. Subpopulations IN-MOC, IN-MAR, CEIEGT, CEPIPSA, and ITC comprise the first cluster. The second axis contributed 13.17% of the variance and separated the UCOL, COLPOS-CORD, Rancho el Paraíso, Universidad Veracruzana, ITC, and IN-MAR subpopulations.

The analysis with the STRUCTURE software demonstrated the presence of two ancestral populations (K), similar to the PCA. The first group included the IN-MOC, IN-MAR, CEIEGT, CEPIPSA, and ITC subpopulations (Figure 4). Hypothetically, if K=3 or K=4, it would be possible to observe the formation of new clusters; however, the subpopulations of the first group remain constant.

This study is the first genetic characterization of Pelibuey sheep using microsatellites. The average NA was 11, which is lower than the 14.27 value reported in a previous study carried out in 13 Colombian sheep breeds using 11 microsatellite markers²⁸; however, it is very similar to the 10.96 value reported in nine Spanish, Cuban, Mexican, and African breeds using 26 microsatellite markers²⁹. The average NE was 3.25, which is lower than the 3.73 observed in five Chinese sheep breeds³⁰ and the 4.68 value reported in Colombian sheep breeds²⁸. When contemplating all loci, the NA was higher than the NE, which indicates that alleles are irregularly distributed in each of the loci; this could be because of the geographic isolation or the selection or flow of genetic material between subpopulations³¹.

• ²⁸
  Genetic diversity of Colombian sheep by microsatellite markers

  Chil J Agric Res, 2016
  Genetic diversity of Colombian sheep by microsatellite markers
• ²⁹
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism

  Small Ruminant Res, 2012
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism
• ³⁰
  Analysis of geographic and pairwise genetic distances among sheep populations

  Biochem Genet, 2010
  Analysis of geographic and pairwise genetic distances among sheep populations
• ²⁸
  Genetic diversity of Colombian sheep by microsatellite markers

  Chil J Agric Res, 2016
  Genetic diversity of Colombian sheep by microsatellite markers
• ³¹
  Study of polymorphisms of microsatellite DNA of six Chinese indigenous sheep and goat breeds

  Front Agr China, 2007
  Study of polymorphisms of microsatellite DNA of six Chinese indigenous sheep and goat breeds

Heterozygosity is a parameter that reflects the genetic diversity in the population. The average values of HO and HE (0.726 and 0.732, respectively) were very similar in all loci. Values are similar to those observed in several Spanish, Cuban, Mexican, and African sheep breeds, in which an average HE value of 0.731 was reported²⁹. A different study reported HO and HE values of 0.744 and 0.755, respectively, in several Hungarian sheep breeds³². In four of the nine analyzed markers (OarFCB304, ILSTS5, SRCRSP5, and MAF33), the HO was higher than the HE, suggesting an excess of heterozygotes in the analyzed subpopulations.

• ²⁹
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism

  Small Ruminant Res, 2012
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism
• ³²
  Genetic diversity and population structure of Zackel sheep and other Hungarian sheep breeds

  Arch Anim Breed, 2015
  Genetic diversity and population structure of Zackel sheep and other Hungarian sheep breeds

The PIC values in this study were greater than 0.7; this indicates that the markers used in this study correctly measure genetic diversity. The average value of the Shannon Index (2.1) was similar to the 2.2 value obtained in 14 types of Iranian sheep³³ and the 2.38 value reported in four Nigerian sheep breeds³⁴. This result reflects the high genetic variability of the analyzed subpopulations.

• ³³
  Using microsatellite markers to analyze genetic diversity in 14 sheep types in Iran

  Arch Anim Breed, 2017
  Using microsatellite markers to analyze genetic diversity in 14 sheep types in Iran
• ³⁴
  Morphological and microsatellite DNA diversity of Nigerian indigenous sheep

  J Anim Sci Biotechnol, 2012
  Morphological and microsatellite DNA diversity of Nigerian indigenous sheep

Three of the nine markers showed a deviation from the HWE due to either high (OarFCB304 and ILSTS5) or low heterozygosity (OarJMP29). Although the deviation from the HWE can have biological implications, such as intensive selection, inbreeding, migration, mutation, or an insufficient number of samples³⁵, it is also possible that this deviation results from errors during genotyping³⁶. In this study, it was only observed deviation from the HWE in 20 of the 216 Fisher's exact tests; therefore, it was decided to keep all markers in subsequent analyses.

• ³⁵
  Probability of random sire exclusion using microsatellite markers for parentage verification

  Anim Genetic, 1995
  Probability of random sire exclusion using microsatellite markers for parentage verification
• ³⁶
  Significant deviations from Hardy-Weimberg equilibrium caused by low levels of microsatellite genotyping errors

  Mol Ecol Res, 2009
  Significant deviations from Hardy-Weimberg equilibrium caused by low levels of microsatellite genotyping errors

The observed FIT value (0.158) indicates a global deficiency of 15.8 % of heterozygotes in the population. Furthermore, the FST value (0.151) indicates that 15.1 % of the genetic variation corresponds to differences between subpopulations. Similar values in both indexes were found in several sheep breeds, where FIT and FST values of 14.2 and 13.4 %, respectively, were observed²⁹. Thus, both values indicate that there is moderate genetic differentiation between the studied flocks.

• ²⁹
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism

  Small Ruminant Res, 2012
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism

The highest NA, MNA, and NE values corresponded to the UNPA, COL-TEX, Rancho Libertad, and Posta El Cuatro subpopulations. The genetic diversity in these subpopulations can be attributed to the introduction of new individuals, which was corroborated with the flocks' record books. Furthermore, Posta El Cuatro corresponds to a unit intended for the production of record Pelibuey specimens. In recent years, pure Pelibuey breeders have substantially improved the productive yield of this breed by resorting to the crosses between different Pelibuey lines. In Rancho Libertad, the high values of FIS (0.138) do not support the latter. However, this result could be attributed to the reduced sample size of the subpopulation. The lowest NA, MNA, and NE were observed in the IN-MOC, IN-MAR, CEIEGT, and CEPIPSA subpopulations; this is probably because these subpopulations belong to universities and research institutes, which have remained as closed centers and have not allowed the introduction of new animals. Thus, the breed's original genetic material is expected to be more conserved in these subpopulations. The negative FIS values confirm that these flocks follow a minimum inbreeding policy.

The HO (0.726) and HE (0.731) values obtained in this study were higher than those observed in a Pelibuey flock in Queretaro, Mexico (HO= 0.652 and HE= 0.659)²⁹. Furthermore, an HO of 0.72 and an HE of 0.71 was observed in a Colombian Pelibuey flock²⁸. The high global values of HO and HE could be attributed to the analysis of several subpopulations since the other studies only evaluated one flock.

• ²⁹
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism

  Small Ruminant Res, 2012
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism
• ²⁸
  Genetic diversity of Colombian sheep by microsatellite markers

  Chil J Agric Res, 2016
  Genetic diversity of Colombian sheep by microsatellite markers

The total value of FIS was 0.0072. However, in nine subpopulations (IN-MOC, UV, Rancho El Paraíso, IN-MAR, ITC, UCOL, CEIEGT, CEPIPSA, and Cuba), FIS was negative, suggesting an excess of heterozygotes. In the remaining subpopulations, FIS was positive, indicating a deficiency of heterozygotes; this coincides with the heterozygosity values since, in the same subpopulations, the HO was higher than the HE. Previous studies in Pelibuey sheep have reported FIS values of 0.023²⁹ and 0.02²⁸, indicating a slight deficiency of heterozygotes. The negative FIS values could be explained by the absence of inbreeding, high heterozygosity values, and low selection pressure, similar to what has been observed in other studies³³. High FIS values (Rancho Garrido, Rancho Libertad, and Rancho El Carrizal) indicate a deficiency of heterozygotes commonly associated with selection and inbreeding³⁷. In farm animal populations, inbreeding is a common finding due to flaws in reproductive programs and the fact that they are relatively small populations. Furthermore, the small number of sampled individuals per subpopulation could also contribute to the estimations obtained in these subpopulations.

• ²⁹
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism

  Small Ruminant Res, 2012
  Genetic relationships of the Cuban hair sheep inferred from microsatellite polymorphism
• ²⁸
  Genetic diversity of Colombian sheep by microsatellite markers

  Chil J Agric Res, 2016
  Genetic diversity of Colombian sheep by microsatellite markers
• ³³
  Using microsatellite markers to analyze genetic diversity in 14 sheep types in Iran

  Arch Anim Breed, 2017
  Using microsatellite markers to analyze genetic diversity in 14 sheep types in Iran
• ³⁷
  Genetic diversity in Kheri-A pastarolist developed Indian sheep using microsatellite markers

  Indian J Biotechnol, 2008
  Genetic diversity in Kheri-A pastarolist developed Indian sheep using microsatellite markers

Private alleles are defined as those that occur in a single population or subpopulation. Nine private alleles were detected in six of the evaluated subpopulations. The subpopulation with the highest number of private alleles (3) was Posta El Cuatro. The CEIEGT, COLPOS-CORD, Rancho San Alberto, Rancho Belbesah, Rancho Libertad, and Rancho Jalapa had one private allele. Except for allele "95" (CEIEGT subpopulation, locus SRCRSP9), which had a 0.400 frequency, private alleles showed relatively low frequencies (0.083-0.125). Thus, these private alleles are more at risk of disappearing if not considered in breeding programs. Private alleles with higher frequencies indicate that a population is unique and has a certain degree of isolation³⁸, which is compatible with the recent history of the CEIEGT flock; this flock has behaved as a closed core, avoiding the introduction of new genetic material. Additionally, the negative FIS (-0.064) rules out the possibility that the high frequency of the allele mentioned above is due to an inbreeding process.

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

  J Anim Sci, 2009
  Genetic variation and population structure of Italian native sheep breeds undergoing in situ conservation

The private alleles observed with low frequencies in the analyzed subpopulations could be due to relatively recent mutation events³⁹. Microsatellites have high mutation rates that range from 10^-6 to 1x10^-2(40; this allows the appearance of a mutation in 1/1,000,000-1/100 individuals per generation. Additionally, private alleles could also be due to the high gene flow between subpopulations⁴¹; this is corroborated by the global FST (0.151) observed in this study, which indicates a moderate gene flow between subpopulations.

• ³⁹
  Contrasting levels of genetic diversity among populations of the endangered tropical palm Euterpe edulis

  CERNE, 2017
  Contrasting levels of genetic diversity among populations of the endangered tropical palm Euterpe edulis
• ⁴⁰
  volutionary dinamics of microsatellite DNA

  Chromosoma, 2000
  40. Schlötterer C. Evolutionary dinamics of microsatellite DNA. Chromosoma 2000;109(6):365-371.
• ⁴¹
  Gene flow in natural populations

  Annu Rev Ecol Syst, 1985
  41. Slatkin M. Gene flow in natural populations. Annu Rev Ecol Syst 1985;16(1):393-430.

The private allele analysis could not differentiate the Cuban subpopulation from the Mexican flocks. It is important to consider the founder effect that resulted from the introduction of the Pelibuey breed into Mexico, in which only a few specimens of the total Cuban population were selected. The genetic drift of the Pelibuey specimens introduced into Mexico played an important role in the extinction or fixation of available alleles. Furthermore, mutations could significantly contribute to the emergence of new alleles within the newly formed flock. Therefore, a substantial difference between the genetic diversity of the Mexican and Cuban flocks would be expected. However, in this study, it was observed that the Cuban flock shares alleles with the Mexican flock; this could be because the Cuban subpopulation consisted of only five individuals, representing a minimum sample of the actual genetic diversity. Therefore, it is impossible to carry out an objective comparison of both subpopulations; results should only be considered indicative.

The analysis of shared alleles did not identify alleles shared between all subpopulations. However, alleles "110" and "150", corresponding to the OarCP34 and OarJMP29 loci, were found in 23 of the 24 analyzed subpopulations; only UV and Rancho Belbesah were missing these alleles. Future research should determine if these alleles could be used as markers of the Pelibuey breed.

The gene flow analysis, measured by the Nm, showed that CEIEGT and CEPIPSA (Nm= 16.39) had the highest gene flow, which corresponds with their record books. These flocks belong to the Universidad Nacional Autónoma de México (UNAM), where rams from CEIEGT have been regularly introduced to CEPIPSA for some years. Furthermore, the Nm analysis confirmed that the value was greater than one in 80.43 % of the subpopulation pairs. Meanwhile, 5.07 % of the subpopulation pairs reached values higher than four. Nm values lower than one indicate that the gene flow is not enough to counteract the differentiation caused by the genetic drift between subpopulations; therefore, subpopulations tend to differentiate. Values greater than one prevent differentiation. Moreover, if the Nm is higher than four, subpopulations behave as a panmictic population⁴²; this is corroborated by the global FST (0.151), which indicates a moderate gene flow between subpopulations.

• ⁴²
  La variabilidad genética de las especies: aspectos conceptuales y sus aplicaciones y perspectivas en México

  Capital natural de México, 2009
  La variabilidad genética de las especies: aspectos conceptuales y sus aplicaciones y perspectivas en México

The correlation between the genetic and geographic distances was nonsignificant; this means that the genetic structure of the Pelibuey flocks does not correspond to isolation by distance model between subpopulations. This could be attributed to the flow of Pelibuey specimens across the country.

The PCA showed that the first component separated the IN-MOC, IN-MAR, CEIEGT, CEPIPSA, and ITC subpopulations from the rest. This result was expected since these populations belong to research institutes and universities, where the main focus is research and breed conservation. The latter was previously corroborated by the negative FIS values that indicate high heterozygosity. The second group was conformed of flocks from other universities and research institutes (BUAP, COL-TEX, COL-COR, UCOL, and UNPA), all the private producers, and the Cuban flock. The divergence of the second group is probably due to the current flow of specimens between Pelibuey subpopulations in Mexico; this is supported by the moderate FST value of 0.151 and the significant proportion of subpopulation pairs with an Nm greater than one. Recently, reproductive technologies, such as embryo transfer and artificial insemination, have enabled faster genetic improvement⁸. Such technologies and the specimen flow (especially rams) between production units have substantially modified the original Pelibuey sheep genetic diversity. Additionally, the crossings with other breeds have significantly increased, resulting in greater divergence. Exhaustive studies should verify the breed's genetic erosion.

• ⁸
  Are cattle, sheep and goats endangered species?

  Mol Ecol, 2008
  8. Taberlet P, Valentini A, Rezaei HR, Naderi S, Pompanon F, Negrini R, Ajmone-Marsan P. Are cattle, sheep and goats endangered species? Mol Ecol 2008;17(1):275-284.

The second component of the PCA separated the UCOL, COLPOS-CORD, Rancho el Paraíso, Universidad Veracruzana, ITC, and IN-MAR subpopulations. These subpopulations belong to research institutes, universities, and private producers. Although it is difficult to explain, this clustering could be due to the practices to which the flocks have been subjected in recent years. For example, the UCOL and IN-MAR populations have drastically decreased; furthermore, Rancho El Paraíso has started introducing new genetic material.

The analysis with the STRUCTURE software suggests that all the Pelibuey subpopulations in Mexico originated from two ancestral populations that have diverged after several years of adaptation to different environments and practices. Furthermore, each of the subpopulations has a mixture of both clusters to a greater or lesser degree. The IN-MOC, IN-MAR, CEIEGT, CEPIPSA, and ITC subpopulations were grouped into one cluster, as observed in the PCA (first component). The relationship between the Mexican reference subpopulations and the other flocks was confirmed with the shared allele analysis previously mentioned. Furthermore, the analysis with STRUCTURE shows low mixture levels between the subpopulations in this cluster and the second one.

Additionally, it has been confirmed that the main objective of those populations is research and breed conservation. If hypothetically is considered a K= 3 or K= 4, it would be possible to observe that the IN-MOC, IN-MAR, CEIEGT, CEPIPSA, and ITC subpopulations continue clustering together; this allows to conclude that these subpopulations have the most preserved genetic material.

The Pelibuey sheep subpopulations in this study showed high genetic diversity and were genetically different from each other. The principal component analysis and the population structure study conclude that the IN-MOC, IN-MAR, CEIEGT, CEPIPSA, and ITC subpopulations have the most preserved genetic material. These results suggest that animals from these subpopulations could be used to implement national conservation programs for the Pelibuey breed. During this study, was observed that despite the significant importance of the Pelibuey breed in sheep farming and research, there are flocks, such as IN-MAR and UCOL, in which the number of individuals has decreased and could potentially disappear. Currently, there are research groups and organizations in Mexico, such as the Unidad Nacional de Ovinocultores [National Unit of Sheep Farmers], that have recognized the importance of the Pelibuey breed in Mexican sheep farming. Through individual efforts, they have tried to contribute to its conservation. A comprehensive plan is required for the conservation of this breed and the preservation of their genetic material. Further studies are essential to determine if there is genetic erosion of the Pelibuey sheep in Mexico.