1
|
Faggion S, Bonfatti V, Carnier P. Genome-Wide Association Study for Weight Loss at the End of Dry-Curing of Hams Produced from Purebred Heavy Pigs. Animals (Basel) 2024; 14:1983. [PMID: 38998095 PMCID: PMC11240668 DOI: 10.3390/ani14131983] [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/31/2024] [Revised: 06/24/2024] [Accepted: 07/03/2024] [Indexed: 07/14/2024] Open
Abstract
Dissecting the genetics of production traits in livestock is of outmost importance, both to understand biological mechanisms underlying those traits and to facilitate the design of selection programs incorporating that information. For the pig industry, traits related to curing are key for protected designation of origin productions. In particular, appropriate ham weight loss after dry-curing ensures high quality of the final product and avoids economic losses. In this study, we analyzed data (N = 410) of ham weight loss after approximately 20 months of dry-curing. The animals used for ham production were purebred pigs belonging to a commercial line. A genome-wide association study (GWAS) of 29,844 SNP markers revealed the polygenic nature of the trait: 221 loci explaining a small percentage of the variance (0.3-1.65%) were identified on almost all Sus scrofa chromosomes. Post-GWAS analyses revealed 32 windows located within regulatory regions and 94 windows located in intronic regions of specific genes. In total, 30 candidate genes encoding receptors and enzymes associated with ham weight loss (MTHFD1L, DUSP8), proteolysis (SPARCL1, MYH8), drip loss (TNNI2), growth (CDCA3, LSP1, CSMD1, AP2A2, TSPAN4), and fat metabolism (AGPAT4, IGF2R, PTDSS2, HRAS, TALDO1, BRSK2, TNNI2, SYT8, GTF2I, GTF2IRD1, LPCAT3, ATN1, GNB3, CMIP, SORCS2, CCSER1, SPP1) were detected.
Collapse
Affiliation(s)
- Sara Faggion
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università 16, 35020 Padova, Italy
| | - Valentina Bonfatti
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università 16, 35020 Padova, Italy
| | - Paolo Carnier
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell'Università 16, 35020 Padova, Italy
| |
Collapse
|
2
|
Ajasa AA, Boison SA, Gjøen HM, Lillehammer M. Accuracy of genomic prediction using multiple Atlantic salmon populations. Genet Sel Evol 2024; 56:38. [PMID: 38750427 PMCID: PMC11094890 DOI: 10.1186/s12711-024-00907-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 05/06/2024] [Indexed: 05/19/2024] Open
Abstract
BACKGROUND The accuracy of genomic prediction is partly determined by the size of the reference population. In Atlantic salmon breeding programs, four parallel populations often exist, thus offering the opportunity to increase the size of the reference set by combining these populations. By allowing a reduction in the number of records per population, multi-population prediction can potentially reduce cost and welfare issues related to the recording of traits, particularly for diseases. In this study, we evaluated the accuracy of multi- and across-population prediction of breeding values for resistance to amoebic gill disease (AGD) using all single nucleotide polymorphisms (SNPs) on a 55K chip or a selected subset of SNPs based on the signs of allele substitution effect estimates across populations, using both linear and nonlinear genomic prediction (GP) models in Atlantic salmon populations. In addition, we investigated genetic distance, genetic correlation estimated based on genomic relationships, and persistency of linkage disequilibrium (LD) phase across these populations. RESULTS The genetic distance between populations ranged from 0.03 to 0.07, while the genetic correlation ranged from 0.19 to 0.99. Nonetheless, compared to within-population prediction, there was limited or no impact of combining populations for multi-population prediction across the various models used or when using the selected subset of SNPs. The estimates of across-population prediction accuracy were low and to some extent proportional to the genetic correlation estimates. The persistency of LD phase between adjacent markers across populations using all SNP data ranged from 0.51 to 0.65, indicating that LD is poorly conserved across the studied populations. CONCLUSIONS Our results show that a high genetic correlation and a high genetic relationship between populations do not guarantee a higher prediction accuracy from multi-population genomic prediction in Atlantic salmon.
Collapse
Affiliation(s)
- Afees A Ajasa
- Nofima (Norwegian Institute of Food, Fisheries and Aquaculture Research), PO Box 210, 1431, Ås, Norway.
- Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, 1430, Ås, Norway.
| | | | - Hans M Gjøen
- Department of Animal and Aquacultural Sciences, Norwegian University of Life Sciences, 1430, Ås, Norway
| | - Marie Lillehammer
- Nofima (Norwegian Institute of Food, Fisheries and Aquaculture Research), PO Box 210, 1431, Ås, Norway
| |
Collapse
|
3
|
Liu P, Liang Y, Li L, Lv X, He Z, Gu Y. Identification of Selection Signatures and Candidate Genes Related to Environmental Adaptation and Economic Traits in Tibetan Pigs. Animals (Basel) 2024; 14:654. [PMID: 38396622 PMCID: PMC10886212 DOI: 10.3390/ani14040654] [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: 12/12/2023] [Revised: 02/01/2024] [Accepted: 02/13/2024] [Indexed: 02/25/2024] Open
Abstract
Tibetan pigs are indigenous to the Qinghai-Tibet Plateau and have been the subject of extensive genomic research primarily focused on their adaptation to high altitudes. However, genetic modifications associated with their response to low-altitude acclimation have not been thoroughly explored. To investigate the genetic basis underlying the low-altitude acclimation of Tibetan pigs, we generated and analyzed genotyping data of Tibetan pigs that inhabit high-altitude regions (average altitude 4000 m) and Tibetan pigs that have inhabited nearby low-altitude regions (average altitude 500 m) for approximately 20 generations. We found that the highland and lowland Tibetan pigs have distinguishable genotype and phenotype variations. We identified 46 and 126 potentially selected SNPs associated with 29 and 56 candidate genes in highland and lowland Tibetan pigs, respectively. Candidate genes in the highland Tibetan pigs were involved in immune response (NFYC and STAT1) and radiation (NABP1), whereas candidate genes in the lowland Tibetan pigs were related to reproduction (ESR2, DMRTA1, and ZNF366), growth and development (NTRK3, FGF18, and MAP1B), and blood pressure regulation (CARTPT). These findings will help to understand the mechanisms of environmental adaptation in Tibetan pigs and offer valuable information into the genetic improvement of Tibetan pigs pertaining to low-altitude acclimation and economic traits.
Collapse
Affiliation(s)
- Pengliang Liu
- Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization, Ministry of Education, Southwest Minzu University, Chengdu 610041, China;
| | - Yan Liang
- Animal Breeding and Genetics Key Laboratory of Sichuan Province, Sichuan Animal Science Academy, Chengdu 610066, China; (Y.L.)
| | - Li Li
- Renshou County Bureau of Agriculture and Rural Affairs, Meishan 620500, China
| | - Xuebin Lv
- Animal Breeding and Genetics Key Laboratory of Sichuan Province, Sichuan Animal Science Academy, Chengdu 610066, China; (Y.L.)
| | - Zhiping He
- Animal Breeding and Genetics Key Laboratory of Sichuan Province, Sichuan Animal Science Academy, Chengdu 610066, China; (Y.L.)
| | - Yiren Gu
- Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization, Ministry of Education, Southwest Minzu University, Chengdu 610041, China;
| |
Collapse
|
4
|
Zhao W, Zhang Z, Wang Z, Ma P, Pan Y, Wang Q, Zhang Z. Factors affecting the accuracy of genomic prediction in joint pig populations. Animal 2023; 17:100980. [PMID: 37797495 DOI: 10.1016/j.animal.2023.100980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2023] [Revised: 08/28/2023] [Accepted: 08/31/2023] [Indexed: 10/07/2023] Open
Abstract
Genomic prediction (GP) has greatly advanced animal and plant breeding over the past two decades. GP in joint populations is a feasible method to improve the accuracy of genomic estimated breeding values in small populations. However, there is still a need to understand the factors that influence GP in joint populations. This study used simulated data and real data from Duroc pig populations to examine the impact of linkage disequilibrium (LD), causal variants effect sizes (CVESs), and minor allele frequencies (MAF) of SNPs on the accuracy of genomic prediction in joint populations. Three prediction methods were used: genomic best linear unbiased prediction (GBLUP), single-step GBLUP and multi-trait GBLUP. Results from the simulated datasets showed that the accuracies of GP in joint populations were always higher than those in a single population when only LD inconsistencies existed. However, single-step GBLUP accuracy in joint populations decreased as the correlation of MAF between populations decreased, while the accuracy of GBLUP is consistently higher in joint populations than in a single population. As the correlation of CVES between populations decreased, the accuracy of both GBLUP and single-step GBLUP in joint populations declined. Analysis of real Duroc populations showed low genetic correlation, similar to the simulated relationship between the most distant populations. In most cases in Duroc populations, GP have higher accuracies in joint populations than in individual population. In conclusion, the consistency of CVES plays a more important role in multi-population GP. The genetic relatedness of the Duroc populations is so weak that the prediction accuracy of GP in joint populations is reduced in some traits. Multi-trait GBLUP is a competitive method for the joint breeding evaluation.
Collapse
Affiliation(s)
- Wei Zhao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiaotong University, 800# Dongchuan Road, Shang, East 200240, China
| | - Zhenyang Zhang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, 866# Yuhangtang Road, Hangzhou, East 310058, China
| | - Zhen Wang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, 866# Yuhangtang Road, Hangzhou, East 310058, China
| | - Peipei Ma
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiaotong University, 800# Dongchuan Road, Shang, East 200240, China
| | - Yuchun Pan
- Department of Animal Science, College of Animal Sciences, Zhejiang University, 866# Yuhangtang Road, Hangzhou, East 310058, China; Hainan Institute, Zhejiang University, Yongyou Industrial Park, Yazhou Bay Sci-Tech City, Sanya 572000, China
| | - Qishan Wang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, 866# Yuhangtang Road, Hangzhou, East 310058, China; Hainan Institute, Zhejiang University, Yongyou Industrial Park, Yazhou Bay Sci-Tech City, Sanya 572000, China
| | - Zhe Zhang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, 866# Yuhangtang Road, Hangzhou, East 310058, China.
| |
Collapse
|
5
|
Faggion S, Boschi E, Veroneze R, Carnier P, Bonfatti V. Genomic Prediction and Genome-Wide Association Study for Boar Taint Compounds. Animals (Basel) 2023; 13:2450. [PMID: 37570259 PMCID: PMC10417264 DOI: 10.3390/ani13152450] [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: 06/05/2023] [Revised: 07/13/2023] [Accepted: 07/27/2023] [Indexed: 08/13/2023] Open
Abstract
With a perspective future ban on surgical castration in Europe, selecting pigs with reduced ability to accumulate boar taint (BT) compounds (androstenone, indole, skatole) in their tissues seems a promising strategy. BT compound concentrations were quantified in the adipose tissue of 1075 boars genotyped at 29,844 SNPs. Traditional and SNP-based breeding values were estimated using pedigree-based BLUP (PBLUP) and genomic BLUP (GBLUP), respectively. Heritabilities for BT compounds were moderate (0.30-0.52). The accuracies of GBLUP and PBLUP were significantly different for androstenone (0.58 and 0.36, respectively), but comparable for indole and skatole (~0.43 and ~0.47, respectively). Several SNP windows, each explaining a small percentage of the variance of BT compound concentrations, were identified in a genome-wide association study (GWAS). A total of 18 candidate genes previously associated with BT (MX1), reproduction traits (TCF21, NME5, PTGFR, KCNQ1, UMODL1), and fat metabolism (CTSD, SYT8, TNNI2, CD81, EGR1, GIPC2, MIGA1, NEGR1, CCSER1, MTMR2, LPL, ERFE) were identified in the post-GWAS analysis. The large number of genes related to fat metabolism might be explained by the relationship between sexual steroid levels and fat deposition and be partially ascribed to the pig line investigated, which is selected for ham quality and not for lean growth.
Collapse
Affiliation(s)
- Sara Faggion
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell’Università 16, 35020 Padova, Italy; (E.B.); (P.C.); (V.B.)
| | - Elena Boschi
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell’Università 16, 35020 Padova, Italy; (E.B.); (P.C.); (V.B.)
| | - Renata Veroneze
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa 36570-999, Brazil;
| | - Paolo Carnier
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell’Università 16, 35020 Padova, Italy; (E.B.); (P.C.); (V.B.)
| | - Valentina Bonfatti
- Department of Comparative Biomedicine and Food Science, University of Padova, Viale dell’Università 16, 35020 Padova, Italy; (E.B.); (P.C.); (V.B.)
| |
Collapse
|
6
|
Clasen JB, Fikse WF, Su G, Karaman E. Multibreed genomic prediction using summary statistics and a breed-origin-of-alleles approach. Heredity (Edinb) 2023; 131:33-42. [PMID: 37231157 PMCID: PMC10313778 DOI: 10.1038/s41437-023-00619-4] [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: 12/14/2021] [Revised: 04/11/2023] [Accepted: 04/26/2023] [Indexed: 05/27/2023] Open
Abstract
Because of an increasing interest in crossbreeding between dairy breeds in dairy cattle herds, farmers are requesting breeding values for crossbred animals. However, genomically enhanced breeding values are difficult to predict in crossbred populations because the genetic make-up of crossbred individuals is unlikely to follow the same pattern as for purebreds. Furthermore, sharing genotype and phenotype information between breed populations are not always possible, which means that genetic merit (GM) for crossbred animals may be predicted without the information needed from some pure breeds, resulting in low prediction accuracy. This simulation study investigated the consequences of using summary statistics from single-breed genomic predictions for some or all pure breeds in two- and three-breed rotational crosses, rather than their raw data. A genomic prediction model taking into account the breed-origin of alleles (BOA) was considered. Because of a high genomic correlation between the breeds simulated (0.62-0.87), the prediction accuracies using the BOA approach were similar to a joint model, assuming homogeneous SNP effects for these breeds. Having a reference population with summary statistics available from all pure breeds and full phenotype and genotype information from crossbreds yielded almost as high prediction accuracies (0.720-0.768) as having a reference population with full information from all pure breeds and crossbreds (0.753-0.789). Lacking information from the pure breeds yielded much lower prediction accuracies (0.590-0.676). Furthermore, including crossbred animals in a combined reference population also benefitted prediction accuracies in the purebred animals, especially for the smallest breed population.
Collapse
Affiliation(s)
- J B Clasen
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Box 7023, 75007, Uppsala, Sweden.
- Center for Quantitative Genetics and Genomics, Aarhus University, C. F. Møllers Allé 8, DK-8000, Aarhus, Denmark.
| | - W F Fikse
- Växa Sverige, Swedish University of Agricultural Sciences, Ulls väg 26, 756 51, Uppsala, Sweden
| | - G Su
- Center for Quantitative Genetics and Genomics, Aarhus University, C. F. Møllers Allé 8, DK-8000, Aarhus, Denmark
| | - E Karaman
- Center for Quantitative Genetics and Genomics, Aarhus University, C. F. Møllers Allé 8, DK-8000, Aarhus, Denmark
| |
Collapse
|
7
|
Botelho ME, Lopes MS, Mathur PK, Knol EF, e Silva FF, Lopes PS, Gimarães SEF, Marques DB, Veroneze R. Weighted genome-wide association study reveals new candidate genes related to boar taint compounds 1. Livest Sci 2022. [DOI: 10.1016/j.livsci.2022.104845] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
8
|
Verardo LL, E Silva FF, Machado MA, do Carmo Panetto JC, de Lima Reis Faza DR, Otto PI, de Almeida Regitano LC, da Silva LOC, do Egito AA, do Socorro Maués Albuquerque M, Zanella R, da Silva MVGB. Genome-Wide Analyses Reveal the Genetic Architecture and Candidate Genes of Indicine, Taurine, Synthetic Crossbreds, and Locally Adapted Cattle in Brazil. Front Genet 2021; 12:702822. [PMID: 34386042 PMCID: PMC8353373 DOI: 10.3389/fgene.2021.702822] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Accepted: 06/30/2021] [Indexed: 11/24/2022] Open
Abstract
Cattle population history, breeding systems, and geographic subdivision may be reflected in runs of homozygosity (ROH), effective population size (Ne), and linkage disequilibrium (LD) patterns. Thus, the assessment of this information has become essential to the implementation of genomic selection on purebred and crossbred cattle breeding programs. In this way, we assessed the genotype of 19 cattle breeds raised in Brazil belonging to taurine, indicine, synthetic crossbreds, and Iberian-derived locally adapted ancestries to evaluate the overall LD decay patterns, Ne, ROH, and breed composition. We were able to obtain a general overview of the genomic architecture of cattle breeds currently raised in Brazil and other tropical countries. We found that, among the evaluated breeds, different marker densities should be used to improve the genomic prediction accuracy and power of genome-wide association studies. Breeds showing low Ne values indicate a recent inbreeding, also reflected by the occurrence of longer ROH, which demand special attention in the matting schemes to avoid extensive inbreeding. Candidate genes (e.g., ABCA7, PENK, SPP1, IFNAR1, IFNAR2, SPEF2, PRLR, LRRTM1, and LRRTM4) located in the identified ROH islands were evaluated, highlighting biological processes involved with milk production, behavior, rusticity, and fertility. Furthermore, we were successful in obtaining the breed composition regarding the taurine and indicine composition using single-nucleotide polymorphism (SNP) data. Our results were able to observe in detail the genomic backgrounds that are present in each breed and allowed to better understand the various contributions of ancestor breeds to the modern breed composition to the Brazilian cattle.
Collapse
Affiliation(s)
- Lucas Lima Verardo
- Animal Breeding Lab, Department of Animal Science, Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, Brazil
| | | | | | | | | | - Pamela Itajara Otto
- Department of Animal Science, Universidade Federal de Santa Maria, Santa Maria, Brazil
| | | | | | | | | | - Ricardo Zanella
- Department of Veterinary Medicine, Universidade de Passo Fundo, Passo Fundo, Brazil
| | | |
Collapse
|
9
|
Genetic Analyses and Genome-Wide Association Studies on Pathogen Resistance of Bos taurus and Bos indicus Cattle Breeds in Cameroon. Genes (Basel) 2021; 12:genes12070976. [PMID: 34206759 PMCID: PMC8307268 DOI: 10.3390/genes12070976] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Revised: 06/23/2021] [Accepted: 06/24/2021] [Indexed: 11/16/2022] Open
Abstract
Autochthonous taurine and later introduced zebu cattle from Cameroon differ considerably in their resistance to endemic pathogens with little to no reports of the underlying genetic make-up. Breed history and habitat variations are reported to contribute significantly to this diversity worldwide, presumably in Cameroon as well, where locations diverge in climate, pasture, and prevalence of infectious agents. In order to investigate the genetic background, the genotypes of 685 individuals of different Cameroonian breeds were analysed by using the BovineSNP50v3 BeadChip. The variance components including heritability were estimated and genome-wide association studies (GWAS) were performed. Phenotypes were obtained by parasitological screening and categorised in Tick-borne pathogens (TBP), gastrointestinal nematodes (GIN), and onchocercosis (ONC). Estimated heritabilities were low for GIN and TBP (0.079 (se = 0.084) and 0.109 (se = 0.103) respectively) and moderate for ONC (0.216 (se = 0.094)). Further than revealing the quantitative nature of the traits, GWAS identified putative trait-associated genomic regions on five chromosomes, including the chromosomes 11 and 18 for GIN, 20 and 24 for TBP, and 12 for ONC. The results imply that breeding for resistant animals in the cattle population from Northern Cameroon might be possible for the studied pathogens; however, further research in this field using larger datasets will be required to improve the resistance towards pathogen infections, propose candidate genes or to infer biological pathways, as well as the genetic structures of African multi-breed populations.
Collapse
|
10
|
Karaman E, Su G, Croue I, Lund MS. Genomic prediction using a reference population of multiple pure breeds and admixed individuals. Genet Sel Evol 2021; 53:46. [PMID: 34058971 PMCID: PMC8168010 DOI: 10.1186/s12711-021-00637-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 05/11/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND In dairy cattle populations in which crossbreeding has been used, animals show some level of diversity in their origins. In rotational crossbreeding, for instance, crossbred dams are mated with purebred sires from different pure breeds, and the genetic composition of crossbred animals is an admixture of the breeds included in the rotation. How to use the data of such individuals in genomic evaluations is still an open question. In this study, we aimed at providing methodologies for the use of data from crossbred individuals with an admixed genetic background together with data from multiple pure breeds, for the purpose of genomic evaluations for both purebred and crossbred animals. A three-breed rotational crossbreeding system was mimicked using simulations based on animals genotyped with the 50 K single nucleotide polymorphism (SNP) chip. RESULTS For purebred populations, within-breed genomic predictions generally led to higher accuracies than those from multi-breed predictions using combined data of pure breeds. Adding admixed population's (MIX) data to the combined pure breed data considering MIX as a different breed led to higher accuracies. When prediction models were able to account for breed origin of alleles, accuracies were generally higher than those from combining all available data, depending on the correlation of quantitative trait loci (QTL) effects between the breeds. Accuracies varied when using SNP effects from any of the pure breeds to predict the breeding values of MIX. Using those breed-specific SNP effects that were estimated separately in each pure breed, while accounting for breed origin of alleles for the selection candidates of MIX, generally improved the accuracies. Models that are able to accommodate MIX data with the breed origin of alleles approach generally led to higher accuracies than models without breed origin of alleles, depending on the correlation of QTL effects between the breeds. CONCLUSIONS Combining all available data, pure breeds' and admixed population's data, in a multi-breed reference population is beneficial for the estimation of breeding values for pure breeds with a small reference population. For MIX, such an approach can lead to higher accuracies than considering breed origin of alleles for the selection candidates, and using breed-specific SNP effects estimated separately in each pure breed. Including MIX data in the reference population of multiple breeds by considering the breed origin of alleles, accuracies can be further improved. Our findings are relevant for breeding programs in which crossbreeding is systematically applied, and also for populations that involve different subpopulations and between which exchange of genetic material is routine practice.
Collapse
Affiliation(s)
- Emre Karaman
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830, Tjele, Denmark.
| | - Guosheng Su
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830, Tjele, Denmark
| | | | - Mogens S Lund
- Center for Quantitative Genetics and Genomics, Aarhus University, 8830, Tjele, Denmark
| |
Collapse
|
11
|
Short communication: Genome wide association study for gastrointestinal nematodes resistance in Bos taurus x Bos indicus crossbred cattle. Livest Sci 2021. [DOI: 10.1016/j.livsci.2021.104403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
12
|
Botelho ME, Lopes MS, Mathur PK, Knol EF, Guimarães SEF, Marques DBD, Lopes PS, Silva FF, Veroneze R. Applying an association weight matrix in weighted genomic prediction of boar taint compounds. J Anim Breed Genet 2020; 138:442-453. [PMID: 33285013 DOI: 10.1111/jbg.12528] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2020] [Revised: 10/13/2020] [Accepted: 11/14/2020] [Indexed: 12/14/2022]
Abstract
Biological information regarding markers and gene association may be used to attribute different weights for single nucleotide polymorphism (SNP) in genome-wide selection. Therefore, we aimed to evaluate the predictive ability and the bias of genomic prediction using models that allow SNP weighting in the genomic relationship matrix (G) building, with and without incorporating biological information to obtain the weights. Firstly, we performed a genome-wide association studies (GWAS) in data set containing single- (SL) or a multi-line (ML) pig population for androstenone, skatole and indole levels. Secondly, 1%, 2%, 5%, 10%, 30% and 50% of the markers explaining the highest proportions of the genetic variance for each trait were selected to build gene networks through the association weight matrix (AWM) approach. The number of edges in the network was computed and used to derive weights for G (AWM-WssGBLUP). The single-step GBLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used as standard scenarios. All scenarios presented predictive abilities different from zero; however, the great overlap in their confidences interval suggests no differences among scenarios. Most of scenarios of based on AWM provide overestimations for skatole in both SL and ML populations. On the other hand, the skatole and indole prediction were no biased in the ssGBLUP (S1) in both SL and ML populations. Most of scenarios based on AWM provide no biased predictions for indole in both SL and ML populations. In summary, using biological information through AWM matrix and gene networks to derive weights for genomic prediction resulted in no increase in predictive ability for boar taint compounds. In addition, this approach increased the number of analyses steps. Thus, we can conclude that ssGBLUP is most appropriate for the analysis of boar taint compounds in comparison with the weighted strategies used in the present work.
Collapse
Affiliation(s)
- Margareth E Botelho
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Marcos S Lopes
- Topigs Norsvin, Curitiba, Brazil.,Topigs Norsvin Research Center, Beuningen, the Netherlands
| | | | - Egbert F Knol
- Topigs Norsvin Research Center, Beuningen, the Netherlands
| | | | - Daniele B D Marques
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Paulo S Lopes
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Fabyano F Silva
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| | - Renata Veroneze
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Brazil
| |
Collapse
|
13
|
Cai Z, Sarup P, Ostersen T, Nielsen B, Fredholm M, Karlskov-Mortensen P, Sørensen P, Jensen J, Guldbrandtsen B, Lund MS, Christensen OF, Sahana G. Genomic diversity revealed by whole-genome sequencing in three Danish commercial pig breeds. J Anim Sci 2020; 98:5873883. [PMID: 32687196 DOI: 10.1093/jas/skaa229] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Accepted: 07/14/2020] [Indexed: 01/04/2023] Open
Abstract
Whole-genome sequencing of 217 animals from three Danish commercial pig breeds (Duroc, Landrace [LL], and Yorkshire [YY]) was performed. Twenty-six million single-nucleotide polymorphisms (SNPs) and 8 million insertions or deletions (indels) were uncovered. Among the SNPs, 493,099 variants were located in coding sequences, and 29,430 were predicted to have a high functional impact such as gain or loss of stop codon. Using the whole-genome sequence dataset as the reference, the imputation accuracy for pigs genotyped with high-density SNP chips was examined. The overall average imputation accuracy for all biallelic variants (SNP and indel) was 0.69, while it was 0.83 for variants with minor allele frequency > 0.1. This study provides whole-genome reference data to impute SNP chip-genotyped animals for further studies to fine map quantitative trait loci as well as improving the prediction accuracy in genomic selection. Signatures of selection were identified both through analyses of fixation and differentiation to reveal selective sweeps that may have had prominent roles during breed development or subsequent divergent selection. However, the fixation indices did not indicate a strong divergence among these three breeds. In LL and YY, the integrated haplotype score identified genomic regions under recent selection. These regions contained genes for olfactory receptors and oxidoreductases. Olfactory receptor genes that might have played a major role in the domestication were previously reported to have been under selection in several species including cattle and swine.
Collapse
Affiliation(s)
- Zexi Cai
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Pernille Sarup
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Tage Ostersen
- SEGES Danish Pig Research Centre, Copenhagen, Denmark
| | | | - Merete Fredholm
- Department of Veterinary and Animal Sciences, University of Copenhagen, Copenhagen, Denmark
| | | | - Peter Sørensen
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Just Jensen
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Bernt Guldbrandtsen
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Mogens Sandø Lund
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Ole Fredslund Christensen
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| | - Goutam Sahana
- Center for Quantitative Genetics and Genomics, Faculty of Technical Sciences, Aarhus University, Tjele, Denmark
| |
Collapse
|
14
|
Thomasen JR, Liu H, Sørensen AC. Genotyping more cows increases genetic gain and reduces rate of true inbreeding in a dairy cattle breeding scheme using female reproductive technologies. J Dairy Sci 2019; 103:597-606. [PMID: 31733861 DOI: 10.3168/jds.2019-16974] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2019] [Accepted: 09/23/2019] [Indexed: 12/26/2022]
Abstract
Both small dairy cattle populations and dairy cattle populations with a low level of linkage disequilibrium (LD) suffer from low reliability of genomic prediction. In this study, we investigated whether adding more genotyped cows to the reference population influences the rate of genetic gain and rate of inbreeding by affecting the reliability. A standard breeding program with a large reference population and high LD, which mimicked a breeding program for Danish Holstein population, was simulated as a reference. A Danish Jersey population with a small reference population and high LD and a Red Dairy Cattle population with a large reference population and low LD were also simulated. Two additional breeding programs were simulated for Danish Jersey and Red Dairy Cattle populations, where 2,000 additional genotyped cows were included in the population for genomic selection. All 5 simulated breeding programs were initiated by a founder population to generate LD resembling the real LD pattern, followed by a 20-yr conventional progeny-testing scheme with 1,000 or 10,000 genotyped progeny-tested bulls and a 10-yr genomic selection scheme with or without 2,000 additional genotyped cows. Evaluation criteria were annual monetary genetic gain and rate of true inbreeding. Our results showed that adding more genotyped cows to the reference in dairy cattle populations has the potential to increase genetic gain and reduce the rate of inbreeding, regardless of reference population size and level of LD. However, it is still not possible to reach the same genetic gain as in the simulated Danish Holstein population with either a small reference population or low LD. Our results also showed that in a small reference population with high LD, it is difficult to manage inbreeding because of lower accuracy compared with the simulated Danish Holstein population and a smaller number of relevant families to select from. Therefore, breeding strategies need to be chosen to match population size and structure. The rate of true inbreeding is always underestimated by pedigree inbreeding and even more in genomic breeding programs, indicating that some forms of genome-wide inbreeding, instead of pedigree-based inbreeding, should be used to monitor inbreeding when genomic selection is implemented.
Collapse
Affiliation(s)
| | - H Liu
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Faculty of Science and Technology, Aarhus University, DK-8830, Tjele, Denmark.
| | - A C Sørensen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and Genetics, Faculty of Science and Technology, Aarhus University, DK-8830, Tjele, Denmark
| |
Collapse
|
15
|
Silva ÉF, Lopes MS, Lopes PS, Gasparino E. A genome-wide association study for feed efficiency-related traits in a crossbred pig population. Animal 2019; 13:2447-2456. [PMID: 31133085 DOI: 10.1017/s1751731119000910] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Feed efficiency (FE) is one of the most important traits in pig production. However, it is difficult and costly to measure it, limiting the collection of large amount of data for an accurate selection for better FE. Therefore, the identification of single-nucleotide polymorphisms (SNPs) associated with FE-related traits to be used in the genetic evaluation is of great interest of pig breeding programs for increasing the prediction accuracy and the genetic progress of these traits. The objective of this study was to identify SNPs significantly associated with FE-related traits: average daily gain (ADG), average daily feed intake (ADFI) and feed conversion ratio (FCR). We also aimed to identify potential candidate genes for these traits. Phenotypic information recorded on a population of 2386 three-way crossbreed pigs that were genotyped for 51 468 SNPs was used. We identified three loci of quantitative trait (QTL) regions associated with ADG and three QTL regions associated with ADFI; however, no significant association was found for FCR. A false discovery rate (FDR) ≤ 0.005 was used as the threshold for declaring an association as significant. The QTL regions associated with ADG on Sus scrofa chromosome (SSC) 1 were located between 177.01 and 185.47 Mb, which overlaps with the QTL regions for ADFI on SSC1 (173.26 and 185.47 Mb). The other QTL region for ADG was located on SSC12 (2.87 and 3.22 Mb). The most significant SNPs in these QTL regions explained up to 3.26% of the phenotypic variance of these traits. The non-identification of genomic regions associated with FCR can be explained by the complexity of this trait, which is a ratio between ADG and ADFI. Finally, the genes CDH19, CDH7, RNF152, MC4R, PMAIP1, FEM1B and GAA were the candidate genes found in the 1 Mb window around the QTL regions identified in this study. Among them, the MC4R gene (SSC1) has a well-known function related to ADG and ADFI. In this study, we identified three QTL regions for ADG (SSC1 and SSC12) and three for ADFI (SSC1). These regions were previously described in purebred pig populations; however, to our knowledge, this is the first study to confirm the relevance of these QTL regions in a crossbred pig population. The potential use of the SNPs and genes identified in this study in prediction models that combine genomic selection and marker-assisted selection should be evaluated for increasing the prediction accuracy of these traits in this population.
Collapse
Affiliation(s)
- É F Silva
- Departamento de Zootecnia, UEM - Universidade Estadual de Maringá, Av. Colombo, 5790, 87.020-900, Maringá, PR, Brazil
- Topigs Norsvin, Rua Visconde do Rio Branco, 1310 - Sala 52, 80.420-210, Curitiba, PR, Brazil
| | - M S Lopes
- Topigs Norsvin, Rua Visconde do Rio Branco, 1310 - Sala 52, 80.420-210, Curitiba, PR, Brazil
- Topigs Norsvin Research Center, Schoenaker 6, 6641 SZ, Beuningen, the Netherlands
| | - P S Lopes
- Departamento de Zootecnia, UFV - Universidade Federal de Viçosa, Campus Universitário, 36.570-000, Viçosa, MG, Brazil
| | - E Gasparino
- Departamento de Zootecnia, UEM - Universidade Estadual de Maringá, Av. Colombo, 5790, 87.020-900, Maringá, PR, Brazil
| |
Collapse
|
16
|
Duenk P, Calus MPL, Wientjes YCJ, Breen VP, Henshall JM, Hawken R, Bijma P. Validation of genomic predictions for body weight in broilers using crossbred information and considering breed-of-origin of alleles. Genet Sel Evol 2019; 51:38. [PMID: 31286857 PMCID: PMC6613268 DOI: 10.1186/s12711-019-0481-7] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Accepted: 06/25/2019] [Indexed: 02/01/2023] Open
Abstract
Background Pig and poultry breeding programs aim at improving crossbred (CB) performance. Selection response may be suboptimal if only purebred (PB) performance is used to compute genomic estimated breeding values (GEBV) because the genetic correlation between PB and CB performance (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc) is often lower than 1. Thus, it may be beneficial to use information on both PB and CB performance. In addition, the accuracy of GEBV of PB animals for CB performance may improve when the breed-of-origin of alleles (BOA) is considered in the genomic relationship matrix (GRM). Thus, our aim was to compare scenarios where GEBV are computed and validated by using (1) either CB offspring averages or individual CB records for validation, (2) either a PB or CB reference population, and (3) a GRM that either accounts for or ignores BOA in the CB individuals. For this purpose, we used data on body weight measured at around 7 (BW7) or 35 (BW35) days in PB and CB broiler chickens and evaluated the accuracy of GEBV based on the correlation GEBV with phenotypes in the validation population (validation correlation). Results With validation on CB offspring averages, the validation correlation of GEBV of PB animals for CB performance was lower with a CB reference population than with a PB reference population for BW35 (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc = 0.96), and about equal for BW7 (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc = 0.80) when BOA was ignored. However, with validation on individual CB records, the validation correlation was higher with a CB reference population for both traits. The use of a GRM that took BOA into account increased the validation correlation for BW7 but reduced it for BW35. Conclusions We argue that the benefit of using a CB reference population for genomic prediction of PB animals for CB performance should be assessed either by validation on CB offspring averages, or by validation on individual CB records while using a GRM that accounts for BOA in the CB individuals. With this recommendation in mind, our results show that the accuracy of GEBV of PB animals for CB performance was equal to or higher with a CB reference population than with a PB reference population for a trait with an \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc of 0.8, but lower for a trait with an \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc of 0.96. In addition, taking BOA into account was beneficial for a trait with an \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc of 0.8 but not for a trait with an \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$r_{pc}$$\end{document}rpc of 0.96. Electronic supplementary material The online version of this article (10.1186/s12711-019-0481-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Pascal Duenk
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands.
| | - Mario P L Calus
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Yvonne C J Wientjes
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | | | | | - Rachel Hawken
- Cobb-vantress Inc., Siloam Springs, AR, 72761-1030, USA
| | - Piter Bijma
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| |
Collapse
|
17
|
Wang L, Mu Y, Xu L, Li K, Han J, Wu T, Liu L, Gao Q, Xia Y, Hou G, Yang S, He X, Liu GE, Feng S. Genomic Analysis Reveals Specific Patterns of Homozygosity and Heterozygosity in Inbred Pigs. Animals (Basel) 2019; 9:E314. [PMID: 31159442 PMCID: PMC6617223 DOI: 10.3390/ani9060314] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 05/27/2019] [Accepted: 05/28/2019] [Indexed: 11/29/2022] Open
Abstract
The inbred strain of miniature pig is an ideal model for biomedical research due to its high level of homozygosity. In this study, we investigated genetic diversity, relatedness, homozygosity, and heterozygosity using the Porcine SNP60K BeadChip in both inbred and non-inbred Wuzhishan pigs (WZSPs). Our results from multidimensional scaling, admixture, and phylogenetic analyses indicated that the inbred WZSP, with its unique genetic properties, can be utilized as a novel genetic resource for pig genome studies. Inbreeding depression and run of homozygosity (ROH) analyses revealed an average of 61 and 12 ROH regions in the inbred and non-inbred genomes of WZSPs, respectively. By investigating ROH number, length, and distribution across generations, we further briefly studied the impacts of recombination and demography on ROH in these WZSPs. Finally, we explored the SNPs with higher heterozygosity across generations and their potential functional implications in the inbred WZSP. We detected 56 SNPs showing constant heterozygosity with He = 1 across six generations in inbred pigs, while only one was found in the non-inbred population. Among these SNPs, we observed nine SNPs located in swine RefSeq genes, which were found to be involved in signaling and immune processes. Together, our findings indicate that the inbred-specific pattern of homozygosity and heterozygosity in inbred pigs can offer valuable insights for elucidating the mechanisms of inbreeding in farm animals.
Collapse
Affiliation(s)
- Ligang Wang
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Yulian Mu
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Linyang Xu
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Kui Li
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Jianlin Han
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Tianwen Wu
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Lan Liu
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Qian Gao
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Ying Xia
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Guanyu Hou
- Institute of Tropical Crop Variety Resources, Chinese Academy of Tropical Agricultural Sciences, Haikou, Hainan 571101, China.
| | - Shulin Yang
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - Xiaohong He
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| | - George E Liu
- Animal Genomics and Improvement Laboratory, U.S. Department of Agriculture-Agricultural Research Services, Beltsville, MD 20705, USA.
| | - Shutang Feng
- Key Laboratory of Farm Animal Genetic Resources and Utilization of Ministry of Agriculture, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (CAAS), Beijing 100193, China.
| |
Collapse
|
18
|
Wientjes YCJ, Calus MPL, Duenk P, Bijma P. Required properties for markers used to calculate unbiased estimates of the genetic correlation between populations. Genet Sel Evol 2018; 50:65. [PMID: 30547748 PMCID: PMC6295113 DOI: 10.1186/s12711-018-0434-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Accepted: 11/28/2018] [Indexed: 11/27/2022] Open
Abstract
BACKGROUND Generally, populations differ in terms of environmental and genetic factors, which can create differences in allele substitution effects between populations. Therefore, a single genotype may have different additive genetic values in different populations. The correlation between the two additive genetic values of a single genotype in two populations is known as the additive genetic correlation between populations and thus, can differ from 1. Our objective was to investigate whether differences in linkage disequilibrium (LD) and allele frequencies of markers and causal loci between populations affect the bias of the estimated genetic correlation. We simulated two populations that were separated by 50 generations and differed in LD pattern between markers and causal loci, as measured by the LD-statistic r. We used a high marker density to represent a high consistency of LD between populations, and lower marker densities to represent situations with a lower consistency of LD between populations. Markers and causal loci were selected to have either similar or different allele frequencies in the two populations. RESULTS Our results show that genetic correlations were underestimated only slightly when the difference in allele frequencies between the two populations was similar for the markers and the causal loci. A lower marker density, representing a lower consistency of LD between populations, had only a minor effect on the underestimation of the genetic correlation. When the difference in allele frequencies between the two populations was not similar for markers and causal loci, genetic correlations were severely underestimated. This bias occurred because the markers did not predict accurately the relationships at causal loci. CONCLUSIONS For an unbiased estimation of the genetic correlation between populations, the markers should accurately predict the relationships at the causal loci. To achieve this, it is essential that the difference in allele frequencies between populations is similar for markers and causal loci. Our results show that differences in LD phase between causal loci and markers across populations have little effect on the estimated genetic correlation.
Collapse
Affiliation(s)
- Yvonne C. J. Wientjes
- Animal Breeding and Genomics, Wageningen University and Research, 6700 AH Wageningen, The Netherlands
| | - Mario P. L. Calus
- Animal Breeding and Genomics, Wageningen University and Research, 6700 AH Wageningen, The Netherlands
| | - Pascal Duenk
- Animal Breeding and Genomics, Wageningen University and Research, 6700 AH Wageningen, The Netherlands
| | - Piter Bijma
- Animal Breeding and Genomics, Wageningen University and Research, 6700 AH Wageningen, The Netherlands
| |
Collapse
|
19
|
Boré R, Brito LF, Jafarikia M, Bouquet A, Maignel L, Sullivan B, Schenkel FS. Genomic data reveals large similarities among Canadian and French maternal pig lines. CANADIAN JOURNAL OF ANIMAL SCIENCE 2018. [DOI: 10.1139/cjas-2017-0103] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
Combining reference populations from different countries and breeds could be an affordable way to enlarge the size of the reference populations for genomic prediction of breeding values. Therefore, the main objectives of this study were to assess the genetic diversity within and between two Canadian and French pig breeds (Landrace and Yorkshire) and the genomic relatedness among populations to evaluate the feasibility of an across-country reference population for pig genomic selection. A total of 14 756 pigs were genotyped on two single nucleotide polymorphism (SNP) chip panels (∼65K SNPs). A principal component analysis clearly discriminated Landrace and Yorkshire breeds, and also, but to a lesser extent, the Canadian and French purebred pigs of each breed. Linkage disequilibrium (LD) between adjacent SNPs was similar within Yorkshire populations. However, levels of LD were slightly different for Landrace populations. The consistency of gametic phase was very high between Yorkshire populations (0.96 at 0.05 Mb) and high for Landrace (0.88 at 0.05 Mb). Based on consistency of gametic phase, Canadian and French pig maternal lines are genetically close to each other. These results are promising, as they indicate that the accuracy of estimated genomic breeding values may increase by combining reference populations from the two countries.
Collapse
Affiliation(s)
- Raphael Boré
- Institut de la Filière Porcine, La Motte au Vicomte, BP 35104, Le Rheu, France
| | - Luiz F. Brito
- Centre for Genetic Improvement of Livestock, University of Guelph, 50 Stone Road East, Guelph, ON N1G 2W1, Canada
| | - Mohsen Jafarikia
- Centre for Genetic Improvement of Livestock, University of Guelph, 50 Stone Road East, Guelph, ON N1G 2W1, Canada
- Canadian Centre for Swine Improvement, Central Experimental Farm, Building No. 75, 960 Carling Avenue, Ottawa, ON K1A 0C6, Canada
| | - Alban Bouquet
- Institut de la Filière Porcine, La Motte au Vicomte, BP 35104, Le Rheu, France
| | - Laurence Maignel
- Canadian Centre for Swine Improvement, Central Experimental Farm, Building No. 75, 960 Carling Avenue, Ottawa, ON K1A 0C6, Canada
| | - Brian Sullivan
- Canadian Centre for Swine Improvement, Central Experimental Farm, Building No. 75, 960 Carling Avenue, Ottawa, ON K1A 0C6, Canada
| | - Flávio S. Schenkel
- Centre for Genetic Improvement of Livestock, University of Guelph, 50 Stone Road East, Guelph, ON N1G 2W1, Canada
| |
Collapse
|
20
|
Sevillano CA, ten Napel J, Guimarães SEF, Silva FF, Calus MPL. Effects of alleles in crossbred pigs estimated for genomic prediction depend on their breed-of-origin. BMC Genomics 2018; 19:740. [PMID: 30305017 PMCID: PMC6180412 DOI: 10.1186/s12864-018-5126-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 09/27/2018] [Indexed: 01/18/2023] Open
Abstract
BACKGROUND This study investigated if the allele effect of a given single nucleotide polymorphism (SNP) for crossbred performance in pigs estimated in a genomic prediction model differs depending on its breed-of-origin, and how these are related to estimated effects for purebred performance. RESULTS SNP-allele substitution effects were estimated for a commonly used SNP panel using a genomic best linear unbiased prediction model with breed-specific partial relationship matrices. Estimated breeding values for purebred and crossbred performance were converted to SNP-allele effects by breed-of-origin. Differences between purebred and crossbred, and between breeds-of-origin were evaluated by comparing percentage of variance explained by genomic regions for back fat thickness (BF), average daily gain (ADG), and residual feed intake (RFI). From ten regions explaining most additive genetic variance for crossbred performance, 1 to 5 regions also appeared in the top ten for purebred performance. The proportion of genetic variance explained by a genomic region and the estimated effect of a haplotype in such a region were different depending upon the breed-of-origin. To illustrate underlying mechanisms, we evaluated the estimated effects across breeds-of-origin for haplotypes associated to the melanocortin 4 receptor (MC4R) gene, and for the MC4Rsnp itself which is a missense mutation with a known effect on BF and ADG. Although estimated allele substitution effects of the MC4Rsnp mutation were very similar across breeds, explained genetic variance of haplotypes associated to the MC4R gene using a SNP panel that does not include the mutation, was considerably lower in one of the breeds where the allele frequency of the mutation was the lowest. CONCLUSIONS Similar regions explaining similar additive genetic variance were observed across purebred and crossbred performance. Moreover, there was some overlap across breeds-of-origin between regions that explained relatively large proportions of genetic variance for crossbred performance; albeit that the actual proportion of variance deviated across breeds-of-origin. Results based on a missense mutation in MC4R confirmed that even if a causal locus has similar effects across breeds-of-origin, estimated effects and explained variance in its region using a commonly used SNP panel can strongly depend on the allele frequency of the underlying causal mutation.
Collapse
Affiliation(s)
- Claudia A Sevillano
- Wageningen University & Research Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH 6700 The Netherlands
- Topigs Norsvin Research Center, P.O. Box 43, Beuningen, 6640 AA The Netherlands
| | - Jan ten Napel
- Wageningen University & Research Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH 6700 The Netherlands
| | - Simone E F Guimarães
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Minas 36570-000 Brazil
| | - Fabyano F Silva
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, Minas 36570-000 Brazil
| | - Mario P L Calus
- Wageningen University & Research Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH 6700 The Netherlands
| |
Collapse
|
21
|
Duarte DAS, Fortes MRS, Duarte MDS, Guimarães SEF, Verardo LL, Veroneze R, Ribeiro AMF, Lopes PS, de Resende MDV, Fonseca e Silva F. Genome-wide association studies, meta-analyses and derived gene network for meat quality and carcass traits in pigs. ANIMAL PRODUCTION SCIENCE 2018. [DOI: 10.1071/an16018] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
A large number of quantitative trait loci (QTL) for meat quality and carcass traits has been reported in pigs over the past 20 years. However, few QTL have been validated and the biological meaning of the genes associated to these QTL has been underexploited. In this context, a meta-analysis was performed to compare the significant markers with meta-QTL previously reported in literature. Genome association studies were performed for 12 traits, from which 144 SNPs were found out to be significant (P < 0.05). They were validated in the meta-analysis and used to build the Association Weight Matrix, a matrix framework employed to investigate co-association of pairwise SNP across phenotypes enabling to derive a gene network. A total of 45 genes were selected from the Association Weight Matrix analysis, from which 25 significant transcription factors were identified and used to construct the networks associated to meat quality and carcass traits. These networks allowed the identification of key transcription factors, such as SOX5 and NKX2–5, gene–gene interactions (e.g. ATP5A1, JPH1, DPT and NEDD4) and pathways related to the regulation of adipose tissue metabolism and skeletal muscle development. Validated SNPs and knowledge of key genes driving these important industry traits might assist future strategies in pig breeding.
Collapse
|
22
|
Sevillano CA, Vandenplas J, Bastiaansen JWM, Bergsma R, Calus MPL. Genomic evaluation for a three-way crossbreeding system considering breed-of-origin of alleles. Genet Sel Evol 2017; 49:75. [PMID: 29061123 PMCID: PMC5653471 DOI: 10.1186/s12711-017-0350-1] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2017] [Accepted: 10/10/2017] [Indexed: 12/31/2022] Open
Abstract
BACKGROUND Genomic prediction of purebred animals for crossbred performance can be based on a model that estimates effects of single nucleotide polymorphisms (SNPs) in purebreds on crossbred performance. For crossbred performance, SNP effects might be breed-specific due to differences between breeds in allele frequencies and linkage disequilibrium patterns between SNPs and quantitative trait loci. Accurately tracing the breed-of-origin of alleles (BOA) in three-way crosses is possible with a recently developed procedure called BOA. A model that accounts for breed-specific SNP effects (BOA model), has never been tested empirically on a three-way crossbreeding scheme. Therefore, the objectives of this study were to evaluate the estimates of variance components and the predictive accuracy of the BOA model compared to models in which SNP effects for crossbred performance were assumed to be the same across breeds, using either breed-specific allele frequencies ([Formula: see text] model) or allele frequencies averaged across breeds ([Formula: see text] model). In this study, we used data from purebred and three-way crossbred pigs on average daily gain (ADG), back fat thickness (BF), and loin depth (LD). RESULTS Estimates of variance components for crossbred performance from the BOA model were mostly similar to estimates from models [Formula: see text] and [Formula: see text]. Heritabilities for crossbred performance ranged from 0.24 to 0.46 between traits. Genetic correlations between purebred and crossbred performance ([Formula: see text]) across breeds ranged from 0.30 to 0.62 for ADG and from 0.53 to 0.74 for BF and LD. For ADG, prediction accuracies of the BOA model were higher than those of the [Formula: see text] and [Formula: see text] models, with significantly higher accuracies only for one maternal breed. For BF and LD, prediction accuracies of models [Formula: see text] and [Formula: see text] were higher than those of the BOA model, with no significant differences. Across all traits, models [Formula: see text] and [Formula: see text] yielded similar predictions. CONCLUSIONS The BOA model yielded a higher prediction accuracy for ADG in one maternal breed, which had the lowest [Formula: see text] (0.30). Using the BOA model was especially relevant for traits with a low [Formula: see text]. In all other cases, the use of crossbred information in models [Formula: see text] and [Formula: see text], does not jeopardize predictions and these models are more easily implemented than the BOA model.
Collapse
Affiliation(s)
- Claudia A Sevillano
- Wageningen University & Research Animal Breeding and Genomics, 6700 AH, Wageningen, The Netherlands. .,Topigs Norsvin Research Center, 6640 AA, Beuningen, The Netherlands.
| | - Jeremie Vandenplas
- Wageningen University & Research Animal Breeding and Genomics, 6700 AH, Wageningen, The Netherlands
| | - John W M Bastiaansen
- Wageningen University & Research Animal Breeding and Genomics, 6700 AH, Wageningen, The Netherlands
| | - Rob Bergsma
- Topigs Norsvin Research Center, 6640 AA, Beuningen, The Netherlands
| | - Mario P L Calus
- Wageningen University & Research Animal Breeding and Genomics, 6700 AH, Wageningen, The Netherlands
| |
Collapse
|
23
|
Lopes MS, Bovenhuis H, Hidalgo AM, van Arendonk JAM, Knol EF, Bastiaansen JWM. Genomic selection for crossbred performance accounting for breed-specific effects. Genet Sel Evol 2017. [PMID: 28651536 PMCID: PMC5485705 DOI: 10.1186/s12711-017-0328-z] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Background Breed-specific effects are observed when the same allele of a given genetic marker has a different effect depending on its breed origin, which results in different allele substitution effects across breeds. In such a case, single-breed breeding values may not be the most accurate predictors of crossbred performance. Our aim was to estimate the contribution of alleles from each parental breed to the genetic variance of traits that are measured in crossbred offspring, and to compare the prediction accuracies of estimated direct genomic values (DGV) from a traditional genomic selection model (GS) that are trained on purebred or crossbred data, with accuracies of DGV from a model that accounts for breed-specific effects (BS), trained on purebred or crossbred data. The final dataset was composed of 924 Large White, 924 Landrace and 924 two-way cross (F1) genotyped and phenotyped animals. The traits evaluated were litter size (LS) and gestation length (GL) in pigs. Results The genetic correlation between purebred and crossbred performance was higher than 0.88 for both LS and GL. For both traits, the additive genetic variance was larger for alleles inherited from the Large White breed compared to alleles inherited from the Landrace breed (0.74 and 0.56 for LS, and 0.42 and 0.40 for GL, respectively). The highest prediction accuracies of crossbred performance were obtained when training was done on crossbred data. For LS, prediction accuracies were the same for GS and BS DGV (0.23), while for GL, prediction accuracy for BS DGV was similar to the accuracy of GS DGV (0.53 and 0.52, respectively). Conclusions In this study, training on crossbred data resulted in higher prediction accuracy than training on purebred data and evidence of breed-specific effects for LS and GL was demonstrated. However, when training was done on crossbred data, both GS and BS models resulted in similar prediction accuracies. In future studies, traits with a lower genetic correlation between purebred and crossbred performance should be included to further assess the value of the BS model in genomic predictions.
Collapse
Affiliation(s)
- Marcos S Lopes
- Topigs Norsvin Research Center, P.O. Box 43, 6640 AA, Beuningen, The Netherlands. .,Animal Breeding and Genomics Centre, Wageningen University, 6708 PB, Wageningen, The Netherlands. .,Topigs Norsvin, Curitiba, PR, 80.420-210, Brazil.
| | - Henk Bovenhuis
- Animal Breeding and Genomics Centre, Wageningen University, 6708 PB, Wageningen, The Netherlands
| | - André M Hidalgo
- Animal Breeding and Genomics Centre, Wageningen University, 6708 PB, Wageningen, The Netherlands
| | - Johan A M van Arendonk
- Animal Breeding and Genomics Centre, Wageningen University, 6708 PB, Wageningen, The Netherlands
| | - Egbert F Knol
- Topigs Norsvin Research Center, P.O. Box 43, 6640 AA, Beuningen, The Netherlands
| | - John W M Bastiaansen
- Animal Breeding and Genomics Centre, Wageningen University, 6708 PB, Wageningen, The Netherlands
| |
Collapse
|
24
|
Gene networks for total number born in pigs across divergent environments. Mamm Genome 2017; 28:426-435. [PMID: 28577119 DOI: 10.1007/s00335-017-9696-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2017] [Accepted: 05/23/2017] [Indexed: 10/19/2022]
Abstract
For reproductive traits such as total number born (TNB), variance due to different environments is highly relevant in animal breeding. In this study, we aimed to perform a gene-network analysis for TNB in pigs across different environments using genomic reaction norm models. Thus, based on relevant single-nucleotide polymorphisms and linkage disequilibrium blocks across environments obtained from GWAS, different sets of candidate genes having biological roles linked to TNB were identified. Network analysis across environment levels resulted in gene interactions consistent with known mammal's fertility biology, captured relevant transcription factors for TNB biology and pointing out different sets of candidate genes for TNB in different environments. These findings may have important implication for animal production, as optimal breeding may vary depending on later environments. Based on these results, genomic diversity was identified and inferred across environments highlighting differential genetic control in each scenario.
Collapse
|
25
|
Verardo LL, Lopes MS, Wijga S, Madsen O, Silva FF, Groenen MAM, Knol EF, Lopes PS, Guimarães SEF. After genome-wide association studies: Gene networks elucidating candidate genes divergences for number of teats across two pig populations. J Anim Sci 2017; 94:1446-58. [PMID: 27136004 DOI: 10.2527/jas.2015-9917] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Number of teats (NT) is an important trait affecting both piglet's welfare and the production level of pig farms. Biologically, embryonic mammary gland development requires the coordination of many signaling pathways necessary for the proper development of teats. Several QTL for NT have been identified; however, further analysis is still lacking. Therefore, gene networks derived from genomewide association study (GWAS) results can be used to examine shared pathways and functions of putative candidate genes. Besides, such analyses may also be helpful to understand the genetic diversity between populations for the same trait or traits. In this study, we identified significant SNP for Landrace-based (line C) and Large White-based (line D) dam lines. Besides, gene-transcription factor (TF) networks were constructed aiming to obtain the most likely candidate genes for NT in each line followed by a comparative analysis between both lines to access similarities or dissimilarities at the marker and gene level. We identified 24 and 19 significant SNP (Bayes factor ≥ 100) for lines C and D, respectively. Only 1 significant SNP overlapped both lines. Network analysis illustrated gene interactions consistent with known mammal's breast biology and captured known TF. We observed different sets of putative candidate genes for NT in each line evaluated that may have common effects on the phenotype. Based on these results, we demonstrated the importance of post-GWAS analyses increasing the biological understanding of relevant genes for a complex trait. Moreover, we believe that this genomic diversity across lines should be taken into account, considering breed-specific reference populations for genomic selection.
Collapse
|
26
|
Grossi DA, Jafarikia M, Brito LF, Buzanskas ME, Sargolzaei M, Schenkel FS. Genetic diversity, extent of linkage disequilibrium and persistence of gametic phase in Canadian pigs. BMC Genet 2017; 18:6. [PMID: 28109261 PMCID: PMC5251314 DOI: 10.1186/s12863-017-0473-y] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2016] [Accepted: 01/13/2017] [Indexed: 01/12/2023] Open
Abstract
Background Knowledge on the levels of linkage disequilibrium (LD) across the genome, persistence of gametic phase between breed pairs, genetic diversity and population structure are important parameters for the successful implementation of genomic selection. Therefore, the objectives of this study were to investigate these parameters in order to assess the feasibility of a multi-herd and multi-breed training population for genomic selection in important purebred and crossbred pig populations in Canada. A total of 3,057 animals, representative of the national populations, were genotyped with the Illumina Porcine SNP60 BeadChip (62,163 markers). Results The overall LD (r2) between adjacent SNPs was 0.49, 0.38, 0.40 and 0.31 for Duroc, Landrace, Yorkshire and Crossbred (Landrace x Yorkshire) populations, respectively. The highest correlation of phase (r) across breeds was observed between Crossbred animals and either Landrace or Yorkshire breeds, in which r was approximately 0.80 at 1 Mbp of distance. Landrace and Yorkshire breeds presented r ≥ 0.80 in distances up to 0.1 Mbp, while Duroc breed showed r ≥ 0.80 for distances up to 0.03 Mbp with all other populations. The persistence of phase across herds were strong for all breeds, with r ≥ 0.80 up to 1.81 Mbp for Yorkshire, 1.20 Mbp for Duroc, and 0.70 Mbp for Landrace. The first two principal components clearly discriminate all the breeds. Similar levels of genetic diversity were observed among all breed groups. The current effective population size was equal to 75 for Duroc and 92 for both Landrace and Yorkshire. Conclusions An overview of population structure, LD decay, demographic history and inbreeding of important pig breeds in Canada was presented. The rate of LD decay for the three Canadian pig breeds indicates that genomic selection can be successfully implemented within breeds with the current 60 K SNP panel. The use of a multi-breed training population involving Landrace and Yorkshire to estimate the genomic breeding values of crossbred animals (Landrace × Yorkshire) should be further evaluated. The lower correlation of phase at short distances between Duroc and the other breeds indicates that a denser panel may be required for the use of a multi-breed training population including Duroc. Electronic supplementary material The online version of this article (doi:10.1186/s12863-017-0473-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Daniela A Grossi
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada
| | - Mohsen Jafarikia
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,Canadian Centre for Swine Improvement Inc, Ottawa, Ontario, Canada
| | - Luiz F Brito
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada
| | - Marcos E Buzanskas
- Departamento de Zootecnia, Centro de Ciências Agrárias - Campus II, Universidade Federal da Paraíba, Areia, Paraíba, Brazil
| | - Mehdi Sargolzaei
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,The Semex Alliance, Guelph, Ontario, Canada
| | - Flávio S Schenkel
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.
| |
Collapse
|
27
|
Veroneze R, Lopes PS, Lopes MS, Hidalgo AM, Guimarães SEF, Harlizius B, Knol EF, van Arendonk JAM, Silva FF, Bastiaansen JWM. Accounting for genetic architecture in single- and multipopulation genomic prediction using weights from genomewide association studies in pigs. J Anim Breed Genet 2016; 133:187-96. [PMID: 27174095 DOI: 10.1111/jbg.12202] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2015] [Accepted: 12/18/2015] [Indexed: 11/28/2022]
Abstract
We studied the effect of including GWAS results on the accuracy of single- and multipopulation genomic predictions. Phenotypes (backfat thickness) and genotypes of animals from two sire lines (SL1, n = 1146 and SL3, n = 1264) were used in the analyses. First, GWAS were conducted for each line and for a combined data set (both lines together) to estimate the genetic variance explained by each SNP. These estimates were used to build matrices of weights (D), which was incorporated into a GBLUP method. Single population evaluated with traditional GBLUP had accuracies of 0.30 for SL1 and 0.31 for SL3. When weights were employed in GBLUP, the accuracies for both lines increased (0.32 for SL1 and 0.34 for SL3). When a multipopulation reference set was used in GBLUP, the accuracies were higher (0.36 for SL1 and 0.32 for SL3) than in single-population prediction. In addition, putting together the multipopulation reference set and the weights from the combined GWAS provided even higher accuracies (0.37 for SL1, and 0.34 for SL3). The use of multipopulation predictions and weights estimated from a combined GWAS increased the accuracy of genomic predictions.
Collapse
Affiliation(s)
- R Veroneze
- Departamento de Zootecnia, Universidade Federal de Viçosa, Viçosa, Brazil.,Animal Breeding and Genomics Centre, Wageningen University, Wageningen, the Netherlands
| | - P S Lopes
- Departamento de Zootecnia, Universidade Federal de Viçosa, Viçosa, Brazil
| | - M S Lopes
- Animal Breeding and Genomics Centre, Wageningen University, Wageningen, the Netherlands.,Topigs Norsvin Research Center, Beuningen, the Netherlands
| | - A M Hidalgo
- Animal Breeding and Genomics Centre, Wageningen University, Wageningen, the Netherlands.,Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| | - S E F Guimarães
- Departamento de Zootecnia, Universidade Federal de Viçosa, Viçosa, Brazil
| | - B Harlizius
- Topigs Norsvin Research Center, Beuningen, the Netherlands
| | - E F Knol
- Topigs Norsvin Research Center, Beuningen, the Netherlands
| | - J A M van Arendonk
- Animal Breeding and Genomics Centre, Wageningen University, Wageningen, the Netherlands
| | - F F Silva
- Departamento de Zootecnia, Universidade Federal de Viçosa, Viçosa, Brazil
| | - J W M Bastiaansen
- Animal Breeding and Genomics Centre, Wageningen University, Wageningen, the Netherlands
| |
Collapse
|
28
|
Verardo LL, Silva FF, Lopes MS, Madsen O, Bastiaansen JWM, Knol EF, Kelly M, Varona L, Lopes PS, Guimarães SEF. Revealing new candidate genes for reproductive traits in pigs: combining Bayesian GWAS and functional pathways. Genet Sel Evol 2016; 48:9. [PMID: 26830357 PMCID: PMC4736284 DOI: 10.1186/s12711-016-0189-x] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2015] [Accepted: 01/20/2016] [Indexed: 12/18/2022] Open
Abstract
Background Reproductive traits such as number of stillborn piglets (SB) and number of teats (NT) have been evaluated in many genome-wide association studies (GWAS). Most of these GWAS were performed under the assumption that these traits were normally distributed. However, both SB and NT are discrete (e.g. count) variables. Therefore, it is necessary to test for better fit of other appropriate statistical models based on discrete distributions. In addition, although many GWAS have been performed, the biological meaning of the identified candidate genes, as well as their functional relationships still need to be better understood. Here, we performed and tested a Bayesian treatment of a GWAS model assuming a Poisson distribution for SB and NT in a commercial pig line. To explore the biological role of the genes that underlie SB and NT and identify the most likely candidate genes, we used the most significant single nucleotide polymorphisms (SNPs), to collect related genes and generated gene-transcription factor (TF) networks. Results Comparisons of the Poisson and Gaussian distributions showed that the Poisson model was appropriate for SB, while the Gaussian was appropriate for NT. The fitted GWAS models indicated 18 and 65 significant SNPs with one and nine quantitative trait locus (QTL) regions within which 18 and 57 related genes were identified for SB and NT, respectively. Based on the related TF, we selected the most representative TF for each trait and constructed a gene-TF network of gene-gene interactions and identified new candidate genes. Conclusions Our comparative analyses showed that the Poisson model presented the best fit for SB. Thus, to increase the accuracy of GWAS, counting models should be considered for this kind of trait. We identified multiple candidate genes (e.g. PTP4A2, NPHP1, and CYP24A1 for SB and YLPM1, SYNDIG1L, TGFB3, and VRTN for NT) and TF (e.g. NF-κB and KLF4 for SB and SOX9 and ELF5 for NT), which were consistent with known newborn survival traits (e.g. congenital heart disease in fetuses and kidney diseases and diabetes in the mother) and mammary gland biology (e.g. mammary gland development and body length). Electronic supplementary material The online version of this article (doi:10.1186/s12711-016-0189-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Lucas L Verardo
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, 36570000, Brazil. .,Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Fabyano F Silva
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, 36570000, Brazil.
| | - Marcos S Lopes
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands. .,Topigs Norsvin, Research Center, 6641 SZ, Beuningen, The Netherlands.
| | - Ole Madsen
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - John W M Bastiaansen
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Egbert F Knol
- Topigs Norsvin, Research Center, 6641 SZ, Beuningen, The Netherlands.
| | - Mathew Kelly
- Queensland Alliance for Agriculture and Food Innovation, The University of Queensland, St Lucia, QLD, 4072, Australia.
| | - Luis Varona
- Departamento de Anatomía, Embriología y Genética, Universidad de Zaragoza, 50013, Saragossa, Spain.
| | - Paulo S Lopes
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, 36570000, Brazil.
| | - Simone E F Guimarães
- Department of Animal Science, Universidade Federal de Viçosa, Viçosa, 36570000, Brazil.
| |
Collapse
|
29
|
|
30
|
Veroneze R, Lopes MS, Hidalgo AM, Guimarães SEF, Silva FF, Harlizius B, Lopes PS, Knol EF, M. van Arendonk JA, Bastiaansen JWM. Accuracy of genome-enabled prediction exploring purebred and crossbred pig populations1. J Anim Sci 2015; 93:4684-91. [DOI: 10.2527/jas.2015-9187] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
|
31
|
Accuracy of Predicted Genomic Breeding Values in Purebred and Crossbred Pigs. G3-GENES GENOMES GENETICS 2015; 5:1575-83. [PMID: 26019187 PMCID: PMC4528314 DOI: 10.1534/g3.115.018119] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Genomic selection has been widely implemented in dairy cattle breeding when the aim is to improve performance of purebred animals. In pigs, however, the final product is a crossbred animal. This may affect the efficiency of methods that are currently implemented for dairy cattle. Therefore, the objective of this study was to determine the accuracy of predicted breeding values in crossbred pigs using purebred genomic and phenotypic data. A second objective was to compare the predictive ability of SNPs when training is done in either single or multiple populations for four traits: age at first insemination (AFI); total number of piglets born (TNB); litter birth weight (LBW); and litter variation (LVR). We performed marker-based and pedigree-based predictions. Within-population predictions for the four traits ranged from 0.21 to 0.72. Multi-population prediction yielded accuracies ranging from 0.18 to 0.67. Predictions across purebred populations as well as predicting genetic merit of crossbreds from their purebred parental lines for AFI performed poorly (not significantly different from zero). In contrast, accuracies of across-population predictions and accuracies of purebred to crossbred predictions for LBW and LVR ranged from 0.08 to 0.31 and 0.11 to 0.31, respectively. Accuracy for TNB was zero for across-population prediction, whereas for purebred to crossbred prediction it ranged from 0.08 to 0.22. In general, marker-based outperformed pedigree-based prediction across populations and traits. However, in some cases pedigree-based prediction performed similarly or outperformed marker-based prediction. There was predictive ability when purebred populations were used to predict crossbred genetic merit using an additive model in the populations studied. AFI was the only exception, indicating that predictive ability depends largely on the genetic correlation between PB and CB performance, which was 0.31 for AFI. Multi-population prediction was no better than within-population prediction for the purebred validation set. Accuracy of prediction was very trait-dependent.
Collapse
|