Abstract
In this study, we genetically characterized the Uruguayan pig breed Pampa Rocha. Genetic variability was assessed by analyzing a panel of 25 microsatellite markers from a sample of 39 individuals. Pampa Rocha pigs showed high genetic variability with observed and expected heterozygosities of 0.583 and 0.603, respectively. The mean number of alleles was 5.72. Twenty-four markers were polymorphic, with 95.8% of them in Hardy Weinberg equilibrium. The level of endogamy was low (FIS = 0.0475). A factorial analysis of correspondence was used to assess the genetic differences between Pampa Rocha and other pig breeds; genetic distances were calculated, and a tree was designed to reflect the distance matrix. Individuals were also allocated into clusters. This analysis showed that the Pampa Rocha breed was separated from the other breeds along the first and second axes. The neighbour-joining tree generated by the genetic distances DA showed clustering of Pampa Rocha with the Meishan breed. The allocation of individuals to clusters showed a clear separation of Pampa Rocha pigs. These results provide insights into the genetic variability of Pampa Rocha pigs and indicate that this breed is a well-defined genetic entity.
genetic resources; microsatellites
Introduction
Pampa Rocha pigs are a genetic resource of the Oriental Republic of Uruguay (Castro, 2007Castro G (2007) Situación de los recursos genéticos porcinos locales en Uruguay. Arch Zootec 56:783–788.; Vadell, 2008Vadell A (2008) Una reseña corta sobre la raza criolla de cerdos Pampa Rocha y su utilización en Uruguay. RCPP 15:105–112.). These pigs are characterized by a black coat with six different regions of white hair. The ears are mostly of the celtic type, while the side of the face can be rectilinear or subconcave (Vadell, 2008Vadell A (2008) Una reseña corta sobre la raza criolla de cerdos Pampa Rocha y su utilización en Uruguay. RCPP 15:105–112.). Figure S1 (Supplementary Material) Figure S1 Pampa Rocha pigs. shows the phenotypic characteristics of these pigs. The origin of this breed dates back to the introduction of the first pigs into South America by Spanish and Portuguese colonists from Europe and Asia, with subsequent modern contributions from Berkshire and Poland China breeds (Urioste et al., 2002Urioste J, Vadell A and Barlocco N (2002) El cerdo Pampa Rocha como recurso zoogenético en Uruguay. Aspectos generales. In: III Simposio Iberoamericano sobre la Conservación de los Recursos Zoogenéticos Locales y el Desarrollo Rural Sostenible. Montevideo.). According to the classification of the Food and Agriculture Organization (FAO) of the United Nations, these pigs are not at risk of extinction, but preservation is required (over 1000 breeding animals) (Delgado, 2002Delgado J (2002) Gestión genética de las poblaciones. In: Jornadas Iberoamericanas sobre la Mejora y Conservación de Razas Ganaderas Locales para el Desarrollo Rural Sostenible. Antigua, Guatemala. 12–17.14.). Since 2002, the possibility of extinction of these animals has increased because of the decrease in the number of small farmers in the eastern region of Uruguay. This situation is compounded by the absence of policies from specialized organisations engaged in the conservation of this resource (Vadell, 2011Vadell A (2011) Los cerdos criollos en el ecosistema de los palmares de Rocha. In: Situación y conservación de recursos zoogenéticos porcinos (Curso-Taller de Posgrado). Facultad de Veterinaria - UdelaR. pp 27–30.). The number of animals in this breed is currently unknown.
The characterization of genetic resources is one of the steps in establishing possible conservation strategies. Currently, DNA microsatellite molecular markers are the most commonly used genetic markers for studying genetic variability in domestic species. In recent years, knowledge of the local pig resources of Uruguay has increased, following a similar trend in other countries. The productivity of Pampa Rocha pigs is already well-characterized. The productivity life-span, grazing adaptability, good maternal ability and high milk production of Pampa Rocha females is remarkable. In addition, the superiority of crossbreeding these pigs with Duroc or Large White breeds in terms of the growth stage and weight gain of offspring is also noteworthy (Barlocco and Vadell, 2005Barlocco N and Vadell A (2005) Experiencias en la caracterización del cerdo Pampa-Rocha de Uruguay. Agrociencia 9:495–503.; Vadell et al., 2010Vadell A, Barlocco N and Carballo C (2010) Prolificidad y longevidad productiva de cerdas Pampa Rocha en un sistema de producción al airelibre. RCPP 17:149–153.). However, the genetic characterization of Pampa Rocha is still insufficient. In Uruguay, the first studies using molecular markers in pigs were done by Kelly et al. (2004)Kelly L, Clop A, Vadell A, Nicolini P, Monteverde S, Amills M and Sanchez A (2004) El cerdo Pampa-Rocha como recurso zoogenético en Uruguay. Marcadores moleculares. Veterinaria (Montevideo) 39:15–16. who analyzed a panel of nine microsatellites and the haplotypes of the mitochondrial cytochrome B gene in ten Pampa Rocha pigs. This study revealed high genetic variability that suggested that the most probable origin of this breed was European breeds with the introgresion of Asian breeds.
The aim of this study was to analyze the genetic variability and population structure of Pampa Rocha using microsatellite markers and to determine the genetic relationships of Pampa Rocha with other pig breeds.
Materials and Methods
Collection of samples and DNA extraction
Thirty-nine DNA samples from Pampa Rocha pigs were analyzed. The animals belonged to the breeding group of the Pig Production Unit of the Southern Regional Center of the Experimental Station of the Faculty of Agronomy (University of the Republic) located in Progreso, Canelones, Uruguay. This breeding group was created in 1995 with animals from different farms from eastern Uruguay, where this resource originated. Some localities where Pampa Rocha is found include Castillos, Barrancas, Costas de Pelotas and Valizas. DNA was extracted from hair follicles using Chelex® resin (Walsh et al., 1991Walsh P, Metzger D and Higuchi R (1991) Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10:506–513.). The samples are currently in the DNA bank of the Genetics Area of the Faculty of Veterinary Medicine (University of the Republic).
Microsatellite analysis
A panel of 25 microsatellites recommended by the FAO/ISAG was amplified (FAO, 2011FAO (2011) Molecular Genetic Characterization of Animal Genetic Resources. FAO Animal Production and Health Guidelines. No. 9. FAO, Rome.) in the Laboratory of Molecular Genetics at the University of Cordoba, Spain, where the technique has been optimized based on the protocol used by Martínez et al. (2000)Martínez AM, Delgado JV, Rodero A and Vega Pla JL (2000) Genetic structure of Iberian pig breed using microsatellites. Anim Genet 31:295–301.. Markers were amplified in multiplex polymerase chain reactions (PCRs) using fluorescence-labelled primers. PCR products were separated by electrophoresis on ABI instruments (ABI PRISM 377XL, Applied Biosystems, Foster City, CA) according to manufacturers recommendations, and allele sizing was done by using the internal size standard of GeneScan-400HD ROX (Applied Biosystems, Warrington, UK). Genotyping was done with Genescan Analysis software v.3.1.2, which provides information about the sizes of the studied fragments, and Genotyper 2.5 software was used to identify the different alleles for each marker.
Statistical analysis
The within breed genetic diversity was investigated by calculating the observed (Ho) and expected (He) heterozygosities, polymorphism information content (PIC) and average number of alleles using Microsatellite Toolkit software (S Park, 2001, PhD thesis, University of Dublin, Ireland). The statistical parameter FIS and the Hardy-Weinberg equilibrium value (Fisher’s exact test, using the Markov chain Monte Carlo method; Guo and Thompson, 1992Guo S and Thompson E (1992) Performing the exact test of Hardy-Weinberg proportions for multiple alleles. Biometrics 48:361–372.) were calculated using Genetix 4.05.2 (Belkhir et al., 2004Belkhir K, Borsa P, Chikhi V, Raufaste N and Bonhomme F (2004) Genetix 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II, Montpellier, France.) and Genepop 4.0.10 (Raymond and Rousset, 1995Raymond M and Rousset F (1995) Genepop (Version 3.4): population genetics software for exact tests and ecumenicism. J Hered 86:248–249.), respectively.
To determine the genetic diversity among breeds, the following breeds were
included in the study: Berkshire, Celtic, Duroc, Iberian, Landrace, Large White,
Meishan and Pietrain, to a total of 415 individuals. Information about these
breeds was obtained from the BIOPIG Consortium’s database of the Iberoamerican
Porcine Biodiversity Project. This consortium is part of the Red CONBIAND. These
data have previously been published by Martínez
et al. (2000)Martínez AM, Delgado JV, Rodero A and Vega Pla JL (2000) Genetic
structure of Iberian pig breed using microsatellites. Anim Genet
31:295–301., Sancristobal et al. (2006)Sancristobal M, Chevalet C, Haley CS, Joosten R, Rattink AP,
Harlizius B, Groenen MAM, Amigues Y, Boscher M, Russell G, et
al. (2006) Genetic diversity within and between European pig breeds
using microsatellite markers. Anim Genet 37:189–198., Vicente et al. (2008)Vicente AA, Carolino MI, Sousa MCO, Ginja C, Silva FS, Martinez AM,
Vega-Pla JL, Carolino N and Gama L (2008) Genetic diversity in native and
commercial breeds of pigs in Portugal assessed by microsatellites. J Anim Sci
86:2496–2507. and Gama et al. (2013)Gama LT, Martínez A, Carolino I, Landi V, Delgado JV, Vicente AA,
Vega-Pla JL, Cortés O, Sousa CO and BIOPIG Consortium (2013) Genetic structure,
relationships and admixture with wild relatives in native pig breeds from Iberia
and its islands. Genet Sel Evol 45:e18.. BIOPIG
consortium analyses were done by A. Martínez at the University of Cordoba,
Spain, where the microsatellite markers used have been standardized. This
comparative analysis was done with 37 samples of Pampa Rocha pigs and 23
microsatellites (excluding the CGA and S0026 markers). Factorial correspondence
analysis was done using Genetix v.4.01 (Belkhir
et al., 2004Belkhir K, Borsa P, Chikhi V, Raufaste N and Bonhomme F (2004)
Genetix 4.05, logiciel sous Windows TM pour la génétique des populations.
Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de
Montpellier II, Montpellier, France.). Nei’s genetic distance
DA was calculated (Nei,
1983Nei M (1983) Genetic polymorphism and the role of mutation in
evolution. In: Nei M and Khoen R (eds) Evolution of Genes and Proteins. Sinauer
Press, Sunderland, MA, pp 165–190.) and a neighbour joining tree was built with Populations 1.2.28
software and viewed with Treeview (Page,
1996Page R (1996) TREEVIEW: an application to display phylogenetic trees
on personal computers. Comput Appl Biosci 12:357–358.). The population structure was assessed using Structure v.2.3.4
software, which uses a Bayesian algorithm (Markov chain Monte Carlo) to allocate
individuals to clusters (K) that are more genetically similar. For this
analysis, the markers that were in Hardy-Weinberg disequilibrium were removed.
Burn-in periods of 200,000 and 600,000 iterations were used, assuming K values
from 2 to 12, with five repetitions for each K and an admixture model and
correlated allele frequencies (Falush et
al., 2003Falush D, Stephens M and Pritchard JK (2003) Inference of population
structure using multilocus genotype data: linked loci and correlated allele
frequencies. Genetics 164:1567–1587.). To identify the most likely K modal
distribution, delta K was determined (Evanno
et al., 2005Evanno G, Regnaut S and Goudet J (2005) Detecting the number of
clusters of individuals using the software STRUCTURE: a simulation study. Mol
Ecol 14:2611–2620.) using Harvester software (Earl and vonHoldt, 2012Earl DA and vonHoldt BM (2012) Structure Harvester: A website and
program for visualizing STRUCTURE output and implementing the Evanno method.
Conserv Genet Resour 4:359–361. Core version: vA.1 http://taylor0.biology.ucla.edu/structureHarvester/ (August 27,
2013).
http://taylor0.biology.ucla.edu/structur...
). The corresponding
graphs were constructed with Distruct software (Rosenberg, 2004Rosenberg NA (2004) Distruct: a program for the graphical display of
population structure. Mol Ecol Notes 4:137–138.).
Results
Within-breed diversity
All the loci studied were polymorphic. Microatellites S0386 and S0002 showed the lowest number of alleles (two), whereas S0068 and SW240 had the highest number of alleles (10). The average number of alleles was 5.72 ± 2.46. PIC values ranged from 0.030 for S0386 to 0.819 for CGA. Seventeen of the 25 markers were highly informative (PIC > 0.5), seven were moderately informative (PIC between 0.25 and 0.5) and one was uninformative (S0386, PIC < 0.25). The expected heterozygosities ranged from 0.030 for S0386 to 0.837 for CGA, with a mean of 0.603 ± 0.199, while the observed heterozygosities ranged from 0.031 for S0386 to 0.875 for S0005, with a mean of 0.583 ± 0.208. The FIS statistic ranged from 0.115 for S0227 to 0.401 for S0002, with a confidence interval of 95%. Seventeen of the 25 markers had a positive sign, indicating an excess of homozygotes, while eight had a negative sign. An FIS multilocus of 0.0475 indicated a low value of inbreeding.
Of the 24 microsatellites analyzed using Fisher’s exact test, 95.8% were in Hardy-Weinberg equilibrium, and only the marker SW72 was out of equilibrium (p < 0.01). The microsatellite S0386 was analyzed by the χ2 test, and p = 1.0 was obtained, indicating that the marker was in equilibrium (Table 1).
Between breed diversity
Figure 1 shows the results of the factorial analysis of correspondence. The first three components accounted for nearly 60% of the total variability. In this analysis, Pampa Rocha was separated from the other breeds by the first and second axes. Figure 2 shows the neighbour-joining tree calculated based on Nei’s genetic distance matrix DA (Nei, 1983Nei M (1983) Genetic polymorphism and the role of mutation in evolution. In: Nei M and Khoen R (eds) Evolution of Genes and Proteins. Sinauer Press, Sunderland, MA, pp 165–190.). The tree showed a cluster of all European breeds and a cluster of the Pampa Rocha and Meishan breeds. This grouping has a 97% bootstrap, indicating that both races were far apart. This divergence most likely reflects a common ancestor influence rather than direct introgression.
Neighbour-joining tree based on Nei’s genetic distances (DA). Breed abbreviations: BK – Berkshire, CL – Celtic, DC – Duroc, IB – Iberian, LD – Landrace, LW – Large White, MS – Meishan, P – Pietrain, PR – Pampa Rocha.
Structure v.2.3.4 (Pritchard et al., 2000Pritchard J, Stephens M and Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959.) was used to allocate the individuals into clusters, with each individual being represented by a vertical bar (Figure 3). When the total number of individuals (415) was assigned to K = 2, the Pampa Rocha individuals were grouped together with the Asian breed Meishan. The remaining breeds used for the analysis (Berkshire, Celtic, Duroc, Iberian, Large White, Landrace and Pietrain) were grouped into the other cluster. For K = 5, the Pampa Rocha pigs were separated from Meishan pigs (q > 0.970). This separation was maintained for K values from 5 to 9; at K=9, all nine breeds were separated. An optimal K value of 9 was determined (Figure S2, Supplementary Material Figure S1 Pampa Rocha pigs. ).
Graphical representation of the results generated by Structure software for nine breeds with K = 2 to K = 9. Breed abbreviations: BK – Berkshire, CL – Celtic, DC – Duroc, IB – Iberian, LD – Landrace, LW – Large White, MS – Meishan, P – Pietrain, PR – Pampa Rocha.
Discussion
The average number of alleles found in this study (5.72 ± 2.46) was lower than that in Mexican Hairless pigs (7.07) (Canul et al., 2005Canul M, Sierra A, Martínez A, Ortiz O, Delgado JV, Vega-Pla J and Pérez G (2005) Caracterización genética del cerdo pelón mexicano mediante marcadores moleculares. Arch Zootec 54:267–272.), Cuban Creole pigs (8.2) (Martínez et al., 2005Martínez AM, Pérez Pineda E, Vega-Pla JL, Barba C, Velázquez FJ and Delgado JV (2005) Caracterización genética del cerdo criollo cubano con microsatélites. Arch Zootec 54:369–375.) and North East Argentina Creole pigs (9.25) (MA Revidatti, 2009, PhD thesis, University of Córdoba, Spain), but higher than the average number of alleles found by Martínez et al. (2000)Martínez AM, Delgado JV, Rodero A and Vega Pla JL (2000) Genetic structure of Iberian pig breed using microsatellites. Anim Genet 31:295–301. in some varieties of Iberian pigs (values of 3.44–5.44), except for the black Entrepelado pig variety (5.84). This value for Pampa Rocha was also higher than that found by Vicente et al. (2008)Vicente AA, Carolino MI, Sousa MCO, Ginja C, Silva FS, Martinez AM, Vega-Pla JL, Carolino N and Gama L (2008) Genetic diversity in native and commercial breeds of pigs in Portugal assessed by microsatellites. J Anim Sci 86:2496–2507. in local and commercial pig breeds from Portugal (5.36) and by Gama et al. (2013)Gama LT, Martínez A, Carolino I, Landi V, Delgado JV, Vicente AA, Vega-Pla JL, Cortés O, Sousa CO and BIOPIG Consortium (2013) Genetic structure, relationships and admixture with wild relatives in native pig breeds from Iberia and its islands. Genet Sel Evol 45:e18. in native pig breeds from Iberia and its island (5.14).
The average expected (0.603) and observed (0.583) heterozygosities of Pampa Rocha indicated a high degree of variability. The heterozygosity observed here resembled the values found for other breeds, such as Mexican Hairless pigs (He = 0.635) (Canul et al., 2005Canul M, Sierra A, Martínez A, Ortiz O, Delgado JV, Vega-Pla J and Pérez G (2005) Caracterización genética del cerdo pelón mexicano mediante marcadores moleculares. Arch Zootec 54:267–272.), Cuban Creole pigs (He = 0.653) (Martínez et al., 2005Martínez AM, Pérez Pineda E, Vega-Pla JL, Barba C, Velázquez FJ and Delgado JV (2005) Caracterización genética del cerdo criollo cubano con microsatélites. Arch Zootec 54:369–375.) and North East Argentina Creole pigs (He = 0.682) (MA Revidatti, 2009, PhD thesis, University of Córdoba, Spain). In Brazilian breeds, such as Monteiro, Moura and Piau, the reported He values were 0.573, 0.569 and 0.661, respectively (Sollero et al., 2009Sollero B, Paiva S, Faria D, Guimarães S, Castro S, Egito A, Albuquerque M, Piovezan U, Bertani G and Mariante A (2009) Genetic diversity of Brazilian pig breeds evidenced by microsatellite markers. Livest Sci 123:8–15.). Thus, the He for Pampa Rocha pigs was greater than the He values reported for the first two Brazilian breeds but lower than those for the Piau breed. For Iberian pigs, He ranges from 0.439 (Jabugo Spotted) to 0.567 (Entrepelado) (A Martínez, 2001, PhD thesis, University of Córdoba, Spain), while for Canary Black pigs, the reported He is 0.475 (Martínez et al., 2007Martínez A, Quiroz J, Marques J and Delgado JV (2007) Estudio de la diversidad genética del cerdo Negro Canario con microsatélites de ADN. Arch Zootec 56:425–428.). In local and commercial pig breeds from Portugal the reported He is 0.66 (Vicente et al., 2008Vicente AA, Carolino MI, Sousa MCO, Ginja C, Silva FS, Martinez AM, Vega-Pla JL, Carolino N and Gama L (2008) Genetic diversity in native and commercial breeds of pigs in Portugal assessed by microsatellites. J Anim Sci 86:2496–2507.), and in native pig breeds from Iberia and its island the reported He is 0.534 (Gama et al., 2013Gama LT, Martínez A, Carolino I, Landi V, Delgado JV, Vicente AA, Vega-Pla JL, Cortés O, Sousa CO and BIOPIG Consortium (2013) Genetic structure, relationships and admixture with wild relatives in native pig breeds from Iberia and its islands. Genet Sel Evol 45:e18.). For Pampa Rocha, the heterozygosity found here was similar to that reported by Kelly et al. (2004)Kelly L, Clop A, Vadell A, Nicolini P, Monteverde S, Amills M and Sanchez A (2004) El cerdo Pampa-Rocha como recurso zoogenético en Uruguay. Marcadores moleculares. Veterinaria (Montevideo) 39:15–16., who also observed high genetic variability with an He value of 0.653; however, the latter study was done with a sample of 10 pigs and nine microsatellite markers. In general, Pampa Rocha showed levels of genetic diversity in its geographical context that resulted from a low level of artificial selection in Creole breeds.
The FIS values in the sample of Pampa Rocha pigs ranged from −0.115 to 0.401. The degree of inbreeding determined by the FIS multilocus was low, in agreement with the genetic diversity results.
Of the markers analyzed, 95.8% were in Hardy-Weinberg equilibrium; only marker SW72 was not (p < 0.01). Microsatellite S0386 was analyzed with the χ2 test and p = 1.0 was obtained, indicating that this micro-satellite was in Hardy-Weinberg equilibrium. Overall, these results demonstrated a high genetic stability, and indicated that this breed was not undergoing genetic drift, selection or other genetic forces, such as mutation or migration. Factors that contribute to disequilibrium include selection and inbreeding in the case of closed populations. Selection forces can be ruled out in this case because the breed is not involved in any modern selection program.
The factorial analysis of correspondence (Figure 1) showed that Pampa Rocha was separated from the other breeds studied by the first and second axes. The Meishan breed was also separated from the remaining races. Structure software v.2.3.4 (Pritchard et al., 2000Pritchard J, Stephens M and Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959.) was used to allocate individuals into clusters based on their genetic similarity, with the assumption that gene frequencies were correlated and that populations were mixed and in Hardy-Weinberg equilibrium (the SW72 marker was excluded). The initial clustering of the races Pampa Rocha and Meishan indicated some level of ancestry, which has been shown before. However, upon further analysis, the Pampa Rocha breed proved to be completely different from the other breeds. The separation into a distinct single cluster indicated that Pampa Rocha individuals constituted a well-defined genetic entity.
The genetic distances between Pampa Rocha and the other breeds were high. The lowest estimated genetic distance was between Pampa Rocha and the Celtic breed (0.57), while the greatest estimated genetic distance was between Pampa Rocha and Meishan (0.73). The tree constructed from this genetic distance analysis consisted of a cluster of all European races and a very distant cluster containing Pampa Rocha and Meishan.
These analyses did not demonstrate an influence of Spanish breeds in Pampa Rocha,
which is unexpected because Spanish breeds have been assumed to have a strong
influence on the origin of Creole pig populations. The time gap and accumulation of
genetic differences may partially explain the separation observed in the tree, but
the successive bottleneck effects that occurred during pig colonization of America
could also explain this basic differentiation. Commercial breeds are thought to have
had less influence on the origin of Creole pig populations, except for the Berkshire
and Meishan breeds. The Berkshire breed was not grouped with Pampa Rocha in this
analysis. The tree analysis indicated some influence of the Asian breed on the
origin of Pampa Rocha because when the number of ancestral populations was small (K
= 2, K = 4), Pampa Rocha and Meishan were grouped in the same cluster. This is
consistent with the studies of Kelly et
al. (2004)Kelly L, Clop A, Vadell A, Nicolini P, Monteverde S, Amills M and
Sanchez A (2004) El cerdo Pampa-Rocha como recurso zoogenético en Uruguay.
Marcadores moleculares. Veterinaria (Montevideo) 39:15–16. who determined that the maternal origin of the
Pampa Rocha breed could be European and Asian. These authors identified the most
common haplotype of European wild and domestic pigs and the characteristic haplotype
of Asian breeds, such as the Japanese Wild Boar and Meishan, in Pampa Rocha. This
has also been described in other breeds, such as Berkshire, Canary Islands’ Black
Pig, Duroc, Landrace and Large White, and indicates the participation of Asian
breeds in their formation. These authors concluded that the origin of the Asian
haplotype in Pampa Rocha pigs likely derived from crosses with Poland China or the
Canary Islands’ Black Pig breeds (introduced during Iberian colonization).
Additionally, the most likely origin of Pampa Rocha pigs was European breeds, some
of which could have been introgressed with Chinese breeds, as indicated by
mitochondrial DNA haplotypes. Burgos-Paz et
al. (2013)Burgos-Paz W, Souza CA, Megens HJ, Ramayo-Caldas Y, Melo M,
Lemús-Flores C, Caal E, Soto HW, Martínez R, Alvarez LA, et al.
(2013) Porcine colonization of the Americas: a 60k SNP story. Heredity
110:321–330. affirmed that the pig breeds or races of the
Americas were the result of many independent colonization and introgression events,
possibly including direct Chinese into-gression. With regard to the origin of Pampa
Rocha, neither Asiatic nor European influences could be rejected because both were
strongly supported by the history of the movement of persons and, consequently,
domestic animals during colonization of the Americas (Rodero et al., 1992Rodero A, Delgado JV and Rodero E (1992) Primitive Andalusian
livestock and their implications in the discovery of America. Arch Zootec
42:383–400.; Gómez, 2007Gómez S (2007) El Galeón de Manila en el siglo XVIII. Navíos de la
Carrera de Filipinas. http://www.todoababor.es/articulos/nav_carrerafilipinas.htm
(August 20, 2013).
http://www.todoababor.es/articulos/nav_c...
; González
Fernández, 2011González Fernández M (2011) Proa a la Mar N162: 40–43. Revista de la
Real Liga Naval Española. http://rlne.info/proaenpdf/proa162pdfs/galeonmanila162.pdf
(August 20, 2013).
http://rlne.info/proaenpdf/proa162pdfs/g...
).
Conclusions
The various parameters and values used to quantify genetic variability, such as the high percentage of polymorphic loci and the expected and observed heterozygosities, indicated high genetic variability in the Pampa Rocha pig breed. A comparative analysis of Pampa Rocha and other breeds showed that Pampa Rocha pigs had high genetic differentiation and conformed to a well-defined and differentiated genetic group. The results described here did not relate Pampa Rocha with Iberian pigs, even though the latter are thought to have influenced the origin of Pampa Rocha because of the historical process that Creole breeds experienced. This lack of association may reflect the time gap between these resources in which genetic drift (bottleneck effects) may have generated differentiation. There have been additional contributions from other breeds. The Meishan breed appears to be the most closely related to Pampa Rocha because of its genetic introgression, a situation not observed for the Iberian or Berkshire breeds, even though they are thought to have participated in the formation of the Pampa Rocha breed. The influence of Asiatic genes on the origin of Pampa Rocha could not be rejected because of the introgression of these genes from the Philippines to the Americas through the so called “Navio de Manila” during the 16th century and later (Rodero et al., 1992Rodero A, Delgado JV and Rodero E (1992) Primitive Andalusian livestock and their implications in the discovery of America. Arch Zootec 42:383–400.). Other modern Asiatic introgressions also cannot be discarded.
Supplementary Material
Figure S1Pampa Rocha pigs.
Figure S2Graphical representation of delta K.
Figure S1 - Pampa Rocha pigs.
Figure S2 - Graphical representation of delta K.
This material is available as part of the online article from http://www.scielo.br/gmb.
The present work was funded by the National Agency of Investigation and Innovation (ANII) and the Sectorial Commission of Scientific Investigation (CSIC).
-
Associate Editor: Alexandre Rodrigues Caetano
References
- Barlocco N and Vadell A (2005) Experiencias en la caracterización del cerdo Pampa-Rocha de Uruguay. Agrociencia 9:495–503.
- Belkhir K, Borsa P, Chikhi V, Raufaste N and Bonhomme F (2004) Genetix 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II, Montpellier, France.
- Burgos-Paz W, Souza CA, Megens HJ, Ramayo-Caldas Y, Melo M, Lemús-Flores C, Caal E, Soto HW, Martínez R, Alvarez LA, et al. (2013) Porcine colonization of the Americas: a 60k SNP story. Heredity 110:321–330.
- Canul M, Sierra A, Martínez A, Ortiz O, Delgado JV, Vega-Pla J and Pérez G (2005) Caracterización genética del cerdo pelón mexicano mediante marcadores moleculares. Arch Zootec 54:267–272.
- Castro G (2007) Situación de los recursos genéticos porcinos locales en Uruguay. Arch Zootec 56:783–788.
- Delgado J (2002) Gestión genética de las poblaciones. In: Jornadas Iberoamericanas sobre la Mejora y Conservación de Razas Ganaderas Locales para el Desarrollo Rural Sostenible. Antigua, Guatemala. 12–17.14.
- Evanno G, Regnaut S and Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14:2611–2620.
- Falush D, Stephens M and Pritchard JK (2003) Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164:1567–1587.
- FAO (2011) Molecular Genetic Characterization of Animal Genetic Resources. FAO Animal Production and Health Guidelines. No. 9. FAO, Rome.
- Gama LT, Martínez A, Carolino I, Landi V, Delgado JV, Vicente AA, Vega-Pla JL, Cortés O, Sousa CO and BIOPIG Consortium (2013) Genetic structure, relationships and admixture with wild relatives in native pig breeds from Iberia and its islands. Genet Sel Evol 45:e18.
- Guo S and Thompson E (1992) Performing the exact test of Hardy-Weinberg proportions for multiple alleles. Biometrics 48:361–372.
- Kelly L, Clop A, Vadell A, Nicolini P, Monteverde S, Amills M and Sanchez A (2004) El cerdo Pampa-Rocha como recurso zoogenético en Uruguay. Marcadores moleculares. Veterinaria (Montevideo) 39:15–16.
- Martínez AM, Delgado JV, Rodero A and Vega Pla JL (2000) Genetic structure of Iberian pig breed using microsatellites. Anim Genet 31:295–301.
- Martínez AM, Pérez Pineda E, Vega-Pla JL, Barba C, Velázquez FJ and Delgado JV (2005) Caracterización genética del cerdo criollo cubano con microsatélites. Arch Zootec 54:369–375.
- Martínez A, Quiroz J, Marques J and Delgado JV (2007) Estudio de la diversidad genética del cerdo Negro Canario con microsatélites de ADN. Arch Zootec 56:425–428.
- Nei M (1983) Genetic polymorphism and the role of mutation in evolution. In: Nei M and Khoen R (eds) Evolution of Genes and Proteins. Sinauer Press, Sunderland, MA, pp 165–190.
- Page R (1996) TREEVIEW: an application to display phylogenetic trees on personal computers. Comput Appl Biosci 12:357–358.
- Pritchard J, Stephens M and Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959.
- Raymond M and Rousset F (1995) Genepop (Version 3.4): population genetics software for exact tests and ecumenicism. J Hered 86:248–249.
- Rodero A, Delgado JV and Rodero E (1992) Primitive Andalusian livestock and their implications in the discovery of America. Arch Zootec 42:383–400.
- Rosenberg NA (2004) Distruct: a program for the graphical display of population structure. Mol Ecol Notes 4:137–138.
- Sancristobal M, Chevalet C, Haley CS, Joosten R, Rattink AP, Harlizius B, Groenen MAM, Amigues Y, Boscher M, Russell G, et al. (2006) Genetic diversity within and between European pig breeds using microsatellite markers. Anim Genet 37:189–198.
- Sollero B, Paiva S, Faria D, Guimarães S, Castro S, Egito A, Albuquerque M, Piovezan U, Bertani G and Mariante A (2009) Genetic diversity of Brazilian pig breeds evidenced by microsatellite markers. Livest Sci 123:8–15.
- Urioste J, Vadell A and Barlocco N (2002) El cerdo Pampa Rocha como recurso zoogenético en Uruguay. Aspectos generales. In: III Simposio Iberoamericano sobre la Conservación de los Recursos Zoogenéticos Locales y el Desarrollo Rural Sostenible. Montevideo.
- Vadell A (2008) Una reseña corta sobre la raza criolla de cerdos Pampa Rocha y su utilización en Uruguay. RCPP 15:105–112.
- Vadell A (2011) Los cerdos criollos en el ecosistema de los palmares de Rocha. In: Situación y conservación de recursos zoogenéticos porcinos (Curso-Taller de Posgrado). Facultad de Veterinaria - UdelaR. pp 27–30.
- Vadell A, Barlocco N and Carballo C (2010) Prolificidad y longevidad productiva de cerdas Pampa Rocha en un sistema de producción al airelibre. RCPP 17:149–153.
- Vicente AA, Carolino MI, Sousa MCO, Ginja C, Silva FS, Martinez AM, Vega-Pla JL, Carolino N and Gama L (2008) Genetic diversity in native and commercial breeds of pigs in Portugal assessed by microsatellites. J Anim Sci 86:2496–2507.
- Walsh P, Metzger D and Higuchi R (1991) Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10:506–513.
Internet Resources
- BIOPIG Consortium’s database: Iberoamerican Porcine Biodiversity Project. http://biopig.jimdo.com/ (May 5, 2012).
» http://biopig.jimdo.com/ - Populations 1.2.28 software. http://www.cnrs-gif.fr/pge/bioinfo/populations/ (May 5, 2012).
» http://www.cnrs-gif.fr/pge/bioinfo/populations/ - Structure software v.2.3.4. http://pritch.bsd.uchicago.edu/structure.html (May 12, 2012).
» http://pritch.bsd.uchicago.edu/structure.html - Distruct software. http://www.stanford.edu/group/rosenberglab/distruct.html (May 12, 2012).
» http://www.stanford.edu/group/rosenberglab/distruct.html - Earl DA and vonHoldt BM (2012) Structure Harvester: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour 4:359–361. Core version: vA.1 http://taylor0.biology.ucla.edu/structureHarvester/ (August 27, 2013).
» http://taylor0.biology.ucla.edu/structureHarvester/ - Gómez S (2007) El Galeón de Manila en el siglo XVIII. Navíos de la Carrera de Filipinas. http://www.todoababor.es/articulos/nav_carrerafilipinas.htm (August 20, 2013).
» http://www.todoababor.es/articulos/nav_carrerafilipinas.htm - González Fernández M (2011) Proa a la Mar N162: 40–43. Revista de la Real Liga Naval Española. http://rlne.info/proaenpdf/proa162pdfs/galeonmanila162.pdf (August 20, 2013).
» http://rlne.info/proaenpdf/proa162pdfs/galeonmanila162.pdf
Publication Dates
-
Publication in this collection
Jan-Mar 2015
History
-
Received
09 May 2014 -
Accepted
25 July 2014