1
|
Arias KD, Lee H, Bozzi R, Álvarez I, Gutiérrez JP, Fernandez I, Menéndez J, Beja-Pereira A, Goyache F. Ascertaining the genetic background of the Celtic-Iberian pig strain: A signatures of selection approach. J Anim Breed Genet 2024; 141:96-112. [PMID: 37807719 DOI: 10.1111/jbg.12829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2023] [Revised: 09/22/2023] [Accepted: 09/26/2023] [Indexed: 10/10/2023]
Abstract
Celtic-Iberian pig breeds were majority in Spain and Portugal until the first half of the 20th century. In the 1990s, they were nearly extinct as a result of the introduction of foreign improved pig breeds. Despite its historical importance, the genetic background of the Celtic-Iberian pig strain is poorly documented. In this study, we have identified genomic regions that might contain signatures of selection peculiar of the Celtic-Iberian genetic lineage. A total of 153 DNA samples of Celtic-Iberian pigs (Spanish Gochu Asturcelta and Portuguese Bísara breeds), Iberian pigs (Spanish Iberian and Portuguese Alentejano breeds), Cinta Senese pig, Korean local pig and Cosmopolitan pig (Hampshire, Landrace and Large White individuals) were analysed. A pairwise-comparison approach was applied: the Gochu Asturcelta and the Bísara samples as test populations and the five other pig populations as reference populations. Three different statistics (XP-EHH, FST and ΔDAF) were computed on each comparison. Strict criteria were used to identify selection sweeps in order to reduce the noise brought on by the Gochu Asturcelta and Bísara breeds' severe population bottlenecks. Within test population, SNPs used to construct potential candidate genomic areas under selection were only considered if they were identified in four of ten two-by-two pairwise comparisons and in at least two of three statistics. Genomic regions under selection constructed within test population were subsequently overlapped to construct candidate regions under selection putatively unique to the Celtic-Iberian pig strain. These genomic regions were finally used for enrichment analyses. A total of 39 candidate regions, mainly located on SSC5 and SSC9 and covering 3130.5 kb, were identified and could be considered representative of the ancient genomic background of the Celtic-Iberian strain. Enrichment analysis allowed to identify a total of seven candidate genes (NOL12, LGALS1, PDXP, SH3BP1, GGA1, WIF1, and LYPD6). Other studies reported that the WIF1 gene is associated with ear size, one of the characteristic traits of the Celtic-Iberian pig strain. The function of the other candidate genes could be related to reproduction, adaptation and immunity traits, indirectly fitting with the rusticity of a non-improved pig strain traditionally exploited under semi-extensive conditions.
Collapse
Affiliation(s)
| | | | - Riccardo Bozzi
- DAGRI, Università degli Studi di Firenze, Firenze, Italy
| | | | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Universidad Complutense de Madrid, Madrid, Spain
| | | | - Juan Menéndez
- ACGA, C/ Párroco José Fernández Teral 5A, Avilés, Asturias, Spain
| | - Albano Beja-Pereira
- CIBIO-InBio, Universidade do Porto, Vairão, Portugal
- BIOPOLIS Program in Genomics, Biodiversity and Land Planning, CIBIO, Vairão, Portugal
- DGAOT, Faculty of Sciences, Universidade do Porto, Porto, Portugal
| | | |
Collapse
|
2
|
Arias KD, Gutiérrez JP, Fernández I, Álvarez I, Goyache F. Approaching autozygosity in a small pedigree of Gochu Asturcelta pigs. Genet Sel Evol 2023; 55:74. [PMID: 37880572 PMCID: PMC10601182 DOI: 10.1186/s12711-023-00846-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2023] [Accepted: 10/03/2023] [Indexed: 10/27/2023] Open
Abstract
BACKGROUND In spite of the availability of single nucleotide polymorphism (SNP) array data, differentiation between observed homozygosity and that caused by mating between relatives (autozygosity) introduces major difficulties. Homozygosity estimators show large variation due to different causes, namely, Mendelian sampling, population structure, and differences among chromosomes. Therefore, the ascertainment of how inbreeding is reflected in the genome is still an issue. The aim of this research was to study the usefulness of genomic information for the assessment of genetic diversity in the highly endangered Gochu Asturcelta pig breed. Pedigree depth varied from 0 (founders) to 4 equivalent discrete generations (t). Four homozygosity parameters (runs of homozygosity, FROH; heterozygosity-rich regions, FHRR; Li and Horvitz's, FLH; and Yang and colleague's FYAN) were computed for each individual, adjusted for the variability in the base population (BP; six individuals) and further jackknifed over autosomes. Individual increases in homozygosity (depending on t) and increases in pairwise homozygosity (i.e., increase in the parents' mean) were computed for each individual in the pedigree, and effective population size (Ne) was computed for five subpopulations (cohorts). Genealogical parameters (individual inbreeding, individual increase in inbreeding, and Ne) were used for comparisons. RESULTS The mean F was 0.120 ± 0.074 and the mean BP-adjusted homozygosity ranged from 0.099 ± 0.081 (FLH) to 0.152 ± 0.075 (FYAN). After jackknifing, the mean values were slightly lower. The increase in pairwise homozygosity tended to be twofold higher than the corresponding individual increase in homozygosity values. When compared with genealogical estimates, estimates of Ne obtained using FYAN tended to have low root-mean-squared errors. However, Ne estimates based on increases in pairwise homozygosity using both FROH and FHRR estimates of genomic inbreeding had lower root-mean-squared errors. CONCLUSIONS Parameters characterizing homozygosity may not accurately depict losses of variability in small populations in which breeding policy prohibits matings between close relatives. After BP adjustment, the performance of FROH and FHRR was highly consistent. Assuming that an increase in homozygosity depends only on pedigree depth can lead to underestimating it in populations with shallow pedigrees. An increase in pairwise homozygosity computed from either FROH or FHRR is a promising approach for characterizing autozygosity.
Collapse
Affiliation(s)
- Katherine D Arias
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Universidad Complutense de Madrid, Avda. Puerta de Hierro S/N, 28040, Madrid, Spain
| | - Iván Fernández
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | - Isabel Álvarez
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | - Félix Goyache
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain.
| |
Collapse
|
3
|
Giovannini S, Strillacci MG, Bagnato A, Albertini E, Sarti FM. Genetic and Phenotypic Characteristics of Belted Pig Breeds: A Review. Animals (Basel) 2023; 13:3072. [PMID: 37835678 PMCID: PMC10571877 DOI: 10.3390/ani13193072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2023] [Revised: 09/26/2023] [Accepted: 09/27/2023] [Indexed: 10/15/2023] Open
Abstract
Belted pig breeds have unique, distinguishing phenotypic characteristics. This review summarises the current knowledge on pig breeds displaying a belted coat pattern. Belts of different widths and positions around the animal's trunk characterise specific pig breeds from all around the world. All the breeds included in the present paper have been searched through the FAO domestic animal diversity information system (DAD-IS), Every country was checked to identify all breeds described as having black or red piebald coat pattern variations. Advances in genomic technologies have made it possible to identify the specific genes and genetic markers associated with the belted phenotype and explore the genetic relationships between different local breeds. Thus, the origin, history, and production traits of these breeds, together with all the genomic information related to the mechanism of skin pigmentation, are discussed. By increasing our understanding of these breeds, we can appreciate the richness of our biological and cultural heritage and work to preserve the biodiversity of the world's animals.
Collapse
Affiliation(s)
- Samira Giovannini
- Department of Agricultural, Food and Environmental Sciences, Università degli Studi di Perugia, Borgo XX Giugno 74, 06121 Perugia, Italy; (E.A.); (F.M.S.)
| | - Maria Giuseppina Strillacci
- Department of Veterinary and Animal Science, Università degli Studi di Milano, Via Dell’Università 6, 26900 Lodi, Italy; (M.G.S.); (A.B.)
| | - Alessandro Bagnato
- Department of Veterinary and Animal Science, Università degli Studi di Milano, Via Dell’Università 6, 26900 Lodi, Italy; (M.G.S.); (A.B.)
| | - Emidio Albertini
- Department of Agricultural, Food and Environmental Sciences, Università degli Studi di Perugia, Borgo XX Giugno 74, 06121 Perugia, Italy; (E.A.); (F.M.S.)
| | - Francesca Maria Sarti
- Department of Agricultural, Food and Environmental Sciences, Università degli Studi di Perugia, Borgo XX Giugno 74, 06121 Perugia, Italy; (E.A.); (F.M.S.)
| |
Collapse
|
4
|
Arias KD, Gutiérrez JP, Fernández I, Álvarez I, Goyache F. Copy Number Variation Regions Differing in Segregation Patterns Span Different Sets of Genes. Animals (Basel) 2023; 13:2351. [PMID: 37508128 PMCID: PMC10376189 DOI: 10.3390/ani13142351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 07/17/2023] [Accepted: 07/18/2023] [Indexed: 07/30/2023] Open
Abstract
Copy number variations regions (CNVRs) can be classified either as segregating, when found in both parents, and offspring, or non-segregating. A total of 65 segregating and 31 non-segregating CNVRs identified in at least 10 individuals within a dense pedigree of the Gochu Asturcelta pig breed was subjected to enrichment and functional annotation analyses to ascertain their functional independence and importance. Enrichment analyses allowed us to annotate 1018 and 351 candidate genes within the bounds of the segregating and non-segregating CNVRs, respectively. The information retrieved suggested that the candidate genes spanned by segregating and non-segregating CNVRs were functionally independent. Functional annotation analyses allowed us to identify nine different significantly enriched functional annotation clusters (ACs) in segregating CNVR candidate genes mainly involved in immunity and regulation of the cell cycle. Up to five significantly enriched ACs, mainly involved in reproduction and meat quality, were identified in non-segregating CNVRs. The current analysis fits with previous reports suggesting that segregating CNVRs would explain performance at the population level, whereas non-segregating CNVRs could explain between-individuals differences in performance.
Collapse
Affiliation(s)
- Katherine D Arias
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394 Gijón, Spain
| | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Iván Fernández
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394 Gijón, Spain
| | - Isabel Álvarez
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394 Gijón, Spain
| | - Félix Goyache
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394 Gijón, Spain
| |
Collapse
|
5
|
Arias KD, Pablo Gutiérrez J, Fernandez I, Menéndez-Arias NA, Álvarez I, Goyache F. Segregation patterns and inheritance rate of copy number variations regions assessed in a Gochu Asturcelta pig pedigree. Gene X 2023; 854:147111. [PMID: 36509293 DOI: 10.1016/j.gene.2022.147111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Revised: 11/30/2022] [Accepted: 12/06/2022] [Indexed: 12/14/2022] Open
Abstract
Copy Number Variation Regions (CNVR) were subjected to pedigree analysis to contribute to the understanding of their segregation patterns. Up to 492 Gochu Asturcelta pig individuals forming 478 different parents-offspring trios (61 different families) were genotyped using the Axiom_PigHDv1 Array (658,692 SNPs). CNVR calling, performed using two different platforms (PennCNV and QuantiSNP), allowed to identify a total of 344 candidate CNVR on the 18 porcine autosomes covering about 106.8 Mb of the pig genome. Sixty-nine CNVR were identified, to some extent, in both the parents and the offspring and were classified as segregating CNVR. The other candidate CNVR were called in one or more progeny but in neither parent and classified either as singleton or recurrent de novo CNVR. Segregating CNVR were, on average, larger and more frequent than the recurrent de novo CNVR (444.8 kb vs 287.9 kb long and 34 vs 5 individuals, respectively). In any case, segregating CNVR did not conform to strict Mendelian inheritance patterns: estimates of average paternal and maternal transmission rates ranged from 11.0 % to 13.4 % and mean inheritance rate was below 21 %. This issue should be carefully considered when interpreting the results of CNV studies. Segregating CNVR, present across generations, are unlikely to be artifacts or false positives and can be hypothesized to be important at the population level.
Collapse
Affiliation(s)
| | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Universidad Complutense de Madrid, Avda. Puerta de Hierro s/n, 28040 Madrid, Spain
| | | | | | | | - Félix Goyache
- SERIDA-Deva, Camino de Rioseco 1225, 33394-Gijón, Spain.
| |
Collapse
|
6
|
Understanding Mendelian errors in SNP arrays data using a Gochu Asturcelta pig pedigree: genomic alterations, family size and calling errors. Sci Rep 2022; 12:19686. [PMID: 36385499 PMCID: PMC9668983 DOI: 10.1038/s41598-022-24340-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 11/14/2022] [Indexed: 11/17/2022] Open
Abstract
Up to 478 Gochu Asturcelta pig parents-offspring trios (61 different families) were genotyped using the Axiom_PigHDv1 Array to identify the causes of Mendelian errors (ME). Up to 545,364 SNPs were retained. Up to 40,540 SNPs gathering 292,297 allelic mismatches were identified and were overlapped with SINEs and LINEs (Sscrofa genome 11.1). Copy number variations (CNV) were called using PennCNV. ME were classified into eight different classes according to the trio member ("Trio" meaning no assignment) and the allele on which ME was identified: TrioA/B, FatherA/B, MotherA/B, OffspringA/B. Most ME occurred due to systematic causes: (a) those assigned to the Father, Mother or Offspring occurred by null or partial null alleles characterized by heterozygote deficiency, varied with family size, involved a low number of loci (6506), and gathered most mismatches (228,145); (b) TrioB errors varied with family size, covaried with SINEs, LINEs and CNV, and involved most ME loci (33,483) and mismatches (65,682); and (c) TrioA errors were non-systematic ME with no sampling bias involving 1.2% of mismatches only and a low number of loci (1939). The influence of TrioB errors on the overall genotyping quality may be low and, since CNV vary among populations, their removal should be considered in each particular dataset. ME assignable to the Father, Mother or Offspring may be consistent within technological platforms and may bias severely linkage or association studies. Most ME caused by null or partial null alleles can be removed using heterozygote deficiency without affecting the size of the datasets.
Collapse
|
7
|
Genetic Differentiation among Livestock Breeds—Values for Fst. Animals (Basel) 2022; 12:ani12091115. [PMID: 35565543 PMCID: PMC9103131 DOI: 10.3390/ani12091115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 04/19/2022] [Accepted: 04/25/2022] [Indexed: 12/02/2022] Open
Abstract
Simple Summary The degree of relationship among livestock breeds can be quantified by the Fst statistic, which measures the extent of genetic differentiation between them. An Fst value of 0.1 has often been taken as indicating that two breeds are indeed genetically distinct, but this concept has not been evaluated critically. Here, Fst values have been collated for the six major livestock species: cattle, sheep, goats, pigs, horses, and chickens. These values are remarkably variable both within and between species, demonstrating that Fst > 0.1 is not a reliable criterion for breed distinctiveness. However, the large body of Fst data accumulated in the last 20–30 years represents an untapped database that could contribute to the development of interdisciplinary research involving livestock breeds. Abstract (1) Background: The Fst statistic is widely used to characterize between-breed relationships. Fst = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding Fst values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed Fst were as follows, cattle: microsatellite 0.06–0.12, SNP 0.08–0.15; sheep: microsatellite 0.06–0.10, SNP 0.06–0.17; horses: microsatellite 0.04–0.11, SNP 0.08–0.12; goats: microsatellite 0.04–0.14, SNP 0.08–0.16; pigs: microsatellite 0.06–0.27, SNP 0.15–0.22; chickens: microsatellite 0.05–0.28, SNP 0.08–0.26. (4) Conclusions: (1) Large amounts of Fst data are available for a substantial proportion of the world’s livestock breeds, (2) the value for between-breed Fst of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated Fst data may have value for interdisciplinary research.
Collapse
|
8
|
An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds. Animals (Basel) 2021; 11:ani11072016. [PMID: 34359144 PMCID: PMC8300386 DOI: 10.3390/ani11072016] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Revised: 06/30/2021] [Accepted: 07/01/2021] [Indexed: 12/18/2022] Open
Abstract
Simple Summary The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity. Abstract Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.
Collapse
|
9
|
Menéndez J, Álvarez I, Fernandez I, Menéndez-Arias NA, Goyache F. Assessing performance of single-sample molecular genetic methods to estimate effective population size: empirical evidence from the endangered Gochu Asturcelta pig breed. Ecol Evol 2016; 6:4971-80. [PMID: 27547327 PMCID: PMC4979721 DOI: 10.1002/ece3.2240] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2016] [Revised: 05/11/2016] [Accepted: 05/16/2016] [Indexed: 11/14/2022] Open
Abstract
Estimating effective population size (Ne) using linkage disequilibrium (LD) information (Ne(LD)) has the operational advantage of using a single sample. However, Ne(LD) estimates assume discrete generations and its performance are constrained by demographic issues. However, such concerns have received little empirical attention so far. The pedigree of the endangered Gochu Asturcelta pig breed includes individuals classified into discrete filial generations and individuals with generations overlap. Up to 780 individuals were typed with a set of 17 microsatellites. Performance of Ne(LD) was compared with Ne estimates obtained using genealogical information, molecular coancestry (Ne(M)) and a temporal (two‐sample) method (Ne(JR)). Molecular‐based estimates of Ne exceeded those obtained using pedigree data. Estimates of Ne(LD) for filial generations F3 and F4 (17.0 and 17.3, respectively) were lower and steadier than those obtained using yearly or biannual samplings. Ne(LD) estimated for samples including generations overlap could only be compared with those obtained for the discrete filial generations when sampling span approached a generation interval and demographic correction for bias was applied. Single‐sample Ne(M) estimates were lower than their Ne(LD) counterparts. Ne(M) estimates are likely to partially reflect the number of founders rather than population size. In any case, estimates of LD and molecular coancestry tend to covary and, therefore, Ne(M) and Ne(LD) can hardly be considered independent. Demographically adjusted estimates of Ne(JR) and Ne(LD) took comparable values when: (1) the two samples used for the former were separated by one equivalent to discrete generations in the pedigree and (2) sampling span used for the latter approached a generation interval. Overall, the empirical evidence given in this study suggested that the advantage of using single‐sample methods to obtain molecular‐based estimates of Ne is not clear in operational terms. Estimates of Ne obtained using methods based in molecular information should be interpreted with caution.
Collapse
Affiliation(s)
- Juan Menéndez
- ACGA C/Párroco José Fernández Teral nº 5 A Avilés Asturias 33403 Spain
| | - Isabel Álvarez
- SERIDA-Deva Camino de Rioseco 1225 Gijón Asturias 33394 Spain
| | - Iván Fernandez
- SERIDA-Deva Camino de Rioseco 1225 Gijón Asturias 33394 Spain
| | | | - Félix Goyache
- SERIDA-Deva Camino de Rioseco 1225 Gijón Asturias 33394 Spain
| |
Collapse
|