1
|
Smaragdov MG. Identification of homozygosity-rich regions in the Holstein genome. Vavilovskii Zhurnal Genet Selektsii 2023; 27:471-479. [PMID: 37808215 PMCID: PMC10556852 DOI: 10.18699/vjgb-23-57] [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: 12/04/2022] [Revised: 02/03/2023] [Accepted: 02/27/2023] [Indexed: 10/10/2023] Open
Abstract
In this study, 371 Holstein cows from six herds and 26 Holstein bulls, which were used in these herds, were genotyped by the Illumina BovineSNP50 array. For runs of homozygosity (ROH) identification, consecutive and sliding runs were performed by the detectRUNS and Plink software. The missing calls did not significantly affect the ROH data. The mean number of ROH identified by consecutive runs was 95.4 ± 2.7, and that by sliding runs was 86.0 ± 2.6 in cows, while this number for Holstein bulls was lower 58.9 ± 1.9. The length of the ROH segments varied from 1 Mb to over 16 Mb, with the largest number of ROH having a length of 1-2 Mb. Of the 29 chromosomes, BTA 14, BTA 16, and BTA 7 were the most covered by ROH. The mean coefficient of inbreeding across the herds was 0.111 ± 0.003 and 0.104 ± 0.004 based on consecutive and sliding runs, respectively, and 0.078 ± 0.005 for bulls based on consecutive runs. These values do not exceed those for Holstein cattle in North America. The results of this study confirmed the more accurate identification of ROH by consecutive runs, and also that the number of allowed heterozygous SNPs may have a significant effect on ROH data.
Collapse
Affiliation(s)
- M G Smaragdov
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Science Center for Animal Husbandry, St. Petersburg, Pushkin, Russia
| |
Collapse
|
2
|
Lozada-Soto EA, Tiezzi F, Jiang J, Cole JB, VanRaden PM, Maltecca C. Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle. J Dairy Sci 2022; 105:8956-8971. [PMID: 36153159 DOI: 10.3168/jds.2022-22116] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 06/20/2022] [Indexed: 11/19/2022]
Abstract
Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (FPED), marker-based using a genomic relationship matrix (FGRM), and segment-based using ROH (FROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using FPED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using FGRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using FROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.
Collapse
Affiliation(s)
| | - Francesco Tiezzi
- Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, 50144 Florence, Italy
| | - Jicai Jiang
- Department of Animal Science, North Carolina State University, Raleigh 27607
| | | | - Paul M VanRaden
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Service, USDA, Beltsville, MD 20705
| | - Christian Maltecca
- Department of Animal Science, North Carolina State University, Raleigh 27607
| |
Collapse
|
3
|
Rezende FM, Haile-Mariam M, Pryce JE, Peñagaricano F. Across-country genomic prediction of bull fertility in Jersey dairy cattle. J Dairy Sci 2020; 103:11618-11627. [PMID: 32981736 DOI: 10.3168/jds.2020-18910] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2020] [Accepted: 07/15/2020] [Indexed: 12/11/2022]
Abstract
The use of information across populations is an attractive approach to increase the accuracy of genomic predictions for numerically small breeds and traits that are time-consuming and difficult to measure, such as male fertility in cattle. This study was conducted to evaluate genomic prediction of Jersey bull fertility using an across-country reference population combining records from the United States and Australia. The data set consisted of 1,570 US Jersey bulls with sire conception rate (SCR) records, 603 Australian Jersey bulls with semen fertility value (SFV) records and SNP genotypes for roughly 90,000 loci. Both SCR and SFV are evaluations of service sire fertility based on cow field data, and both are intended as phenotypic evaluations because the estimates include genetic and nongenetic effects. Within- and across-country genomic predictions were evaluated using univariate and bivariate genomic best linear unbiased prediction models. Predictive ability was assessed in 5-fold cross-validation using the correlation between observed and predicted fertility values and mean squared error of prediction. Within-country genomic predictions exhibited predictive correlations of around 0.28 and 0.02 for the United States and Australia, respectively. The Australian Jersey population is genetically diverse and small in size, so careful selection of the reference population by including only closely related animals (e.g., excluding New Zealand bulls, which is a less-related population) increased the predictive correlations up to 0.20. Notably, the use of bivariate models fitting all US Jersey records and the optimized Australian population resulted in predictive correlations around of 0.24 for SFV values, which is a relative increase in predictive ability of 20%. Conversely, for predicting SCR values, the use of an across-country reference population did not outperform the standard approach using pure US Jersey reference data set. Our findings indicate that genomic prediction of male fertility in dairy cattle is feasible, and the use of an across-country reference population would be beneficial when local populations are small and genetically diverse.
Collapse
Affiliation(s)
- Fernanda M Rezende
- Department of Animal Sciences, University of Florida, Gainesville 32611; Faculdade de Medicina Veterinária, Universidade Federal de Uberlândia, Uberlândia MG 38410-337, Brazil
| | - Mekonnen Haile-Mariam
- Agriculture Victoria Research, AgriBio, Centre for AgriBioscience, Bundoora, Victoria 3083, Australia
| | - Jennie E Pryce
- Agriculture Victoria Research, AgriBio, Centre for AgriBioscience, Bundoora, Victoria 3083, Australia; School of Applied Systems Biology, La Trobe University, Bundoora, Victoria 3083, Australia
| | - Francisco Peñagaricano
- Department of Animal Sciences, University of Florida, Gainesville 32611; Department of Animal and Dairy Sciences, University of Wisconsin-Madison, 53706.
| |
Collapse
|
4
|
Cardoso DF, Fernandes Júnior GA, Scalez DCB, Alves AAC, Magalhães AFB, Bresolin T, Ventura RV, Li C, de Sena Oliveira MC, Porto-Neto LR, Carvalheiro R, de Oliveira HN, Tonhati H, Albuquerque LG. Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle. Sci Rep 2020; 10:8770. [PMID: 32471998 PMCID: PMC7260210 DOI: 10.1038/s41598-020-65565-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2019] [Accepted: 05/04/2020] [Indexed: 11/09/2022] Open
Abstract
Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77-65.25 Mb) and 14 (22.81-23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.
Collapse
Affiliation(s)
- Diercles Francisco Cardoso
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil.
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada.
| | - Gerardo Alves Fernandes Júnior
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
| | - Daiane Cristina Becker Scalez
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
- Centre for Genetic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, Canada
| | - Anderson Antonio Carvalho Alves
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
| | - Ana Fabrícia Braga Magalhães
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
| | - Tiago Bresolin
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
| | - Ricardo Vieira Ventura
- Department of Animal Nutrition and Production, School of Veterinary Medicine and Animal Science (FMVZ), University of Sao Paulo (USP), Pirassununga, SP, Brazil
| | - Changxi Li
- Department of Agricultural Food and Nutritional Science, Faculty of Agricultural, Life & Environmental Sciences, University of Alberta, Edmonton, AB, Canada
- Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
| | | | | | - Roberto Carvalheiro
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
- National Council for Science and Technological Development, Brasília, Distrito Federal, Brazil
| | - Henrique Nunes de Oliveira
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
- National Council for Science and Technological Development, Brasília, Distrito Federal, Brazil
| | - Humberto Tonhati
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil
- National Council for Science and Technological Development, Brasília, Distrito Federal, Brazil
| | - Lucia Galvão Albuquerque
- Department of Animal Science, School of Agricultural and Veterinarian Science, São Paulo State University (UNESP), Jaboticabal, SP, Brazil.
- National Council for Science and Technological Development, Brasília, Distrito Federal, Brazil.
| |
Collapse
|
5
|
Smaragdov MG, Kudinov AA. Assessing the power of principal components and wright's fixation index analyzes applied to reveal the genome-wide genetic differences between herds of Holstein cows. BMC Genet 2020; 21:47. [PMID: 32345235 PMCID: PMC7189535 DOI: 10.1186/s12863-020-00848-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 03/27/2020] [Indexed: 11/30/2022] Open
Abstract
Background Due to the advent of SNP array technology, a genome-wide analysis of genetic differences between populations and breeds has become possible at a previously unattainable level. The Wright’s fixation index (Fst) and the principal component analysis (PCA) are widely used methods in animal genetics studies. In paper we compared the power of these methods, their complementing each other and which of them is the most powerful. Results Comparative analysis of the power Principal Components Analysis (PCA) and Fst were carried out to reveal genetic differences between herds of Holsteinized cows. Totally, 803 BovineSNP50 genotypes of cows from 13 herds were used in current study. Obtained Fst values were in the range of 0.002–0.012 (mean 0.0049) while for rare SNPs with MAF 0.0001–0.005 they were even smaller in the range of 0.001–0.01 (mean 0.0027). Genetic relatedness of the cows in the herds was the cause of such small Fst values. The contribution of rare alleles with MAF 0.0001–0.01 to the Fst values was much less than common alleles and this effect depends on linkage disequilibrium (LD). Despite of substantial change in the MAF spectrum and the number of SNPs we observed small effect size of LD - based pruning on Fst data. PCA analysis confirmed the mutual admixture and small genetic difference between herds. Moreover, PCA analysis of the herds based on the visualization the results of a single eigenvector cannot be used to significantly differentiate herds. Only summed eigenvectors should be used to realize full power of PCA to differentiate small between herds genetic difference. Finally, we presented evidences that the significance of Fst data far exceeds the significance of PCA data when these methods are used to reveal genetic differences between herds. Conclusions LD - based pruning had a small effect on findings of Fst and PCA analyzes. Therefore, for weakly structured populations the LD - based pruning is not effective. In addition, our results show that the significance of genetic differences between herds obtained by Fst analysis exceeds the values of PCA. Proposed, to differentiate herds or low structured populations we recommend primarily using the Fst approach and only then PCA.
Collapse
Affiliation(s)
- M G Smaragdov
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the l.K. Ernst Federal Science Center for Animal Husbandry, St. Petersburg, Pushkin, Russia. .,, St. Petersburg, Russian Federation.
| | - A A Kudinov
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the l.K. Ernst Federal Science Center for Animal Husbandry, St. Petersburg, Pushkin, Russia.,Department of Agricultural Science, University of Helsinki, FI-00014, Helsinki, Finland
| |
Collapse
|
6
|
Haile-Mariam M, MacLeod IM, Bolormaa S, Schrooten C, O'Connor E, de Jong G, Daetwyler HD, Pryce JE. Value of sharing cow reference population between countries on reliability of genomic prediction for milk yield traits. J Dairy Sci 2019; 103:1711-1728. [PMID: 31864746 DOI: 10.3168/jds.2019-17170] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Accepted: 10/24/2019] [Indexed: 01/08/2023]
Abstract
Increasing the reliability of genomic prediction (GP) of economic traits in the pasture-based dairy production systems of New Zealand (NZ) and Australia (AU) is important to both countries. This study assessed if sharing cow phenotype and genotype data of NZ and AU improves the reliability of GP for NZ bulls. Data from approximately 32,000 NZ genotyped cows and their contemporaries were included in the May 2018 routine genetic evaluation of the Australian Dairy cattle in an attempt to provide consistent phenotypes for both countries. After the genetic evaluation, deregressed proofs of cows were calculated for milk yield traits. The April 2018 multiple across-country evaluation of Interbull was also used to calculate deregressed proofs for bulls on the NZ scale. Approximately 1,178 Jersey (Jer) and 6,422 Holstein (Hol) bulls had genotype and phenotype data. In addition to NZ cows, phenotype data of close to 60,000 genotyped Australian (AU) cows from the same genetic evaluation run as NZ cows were used. All AU and NZ females were genotyped using low-density SNP chips (<10K SNP) and were imputed first to 50K and then to ∼600K (referred to as high density; HD). We used up to 98,000 animals in the reference populations, both by expanding the NZ reference set (cow, bull, single breed to multi-breed set) and by adding AU cows. Reliabilities of GP were calculated for 508 Jer and 1,251 Hol bulls whose sires are not included in the reference set (RS) to ensure that real differences are not masked by close relationships. The GP was tested using 50K or high-density SNP chip using genomic BLUP in bivariate (considering country as a trait) or single trait models. The RS that gave the highest reliability for each breed were also tested using a hybrid GP method that combines expectation maximization with Bayes R. The addition of the AU cows to an NZ RS that included either NZ cows only, or cows and bulls, improved the reliability of GP for both NZ Hol and Jer validation bulls for all traits. Using single breed reference populations also increased reliability when NZ crossbred cows were added to reference populations that included only purebred NZ bulls and cows and AU cows. The full multi-breed RS (all NZ cows and bulls and AU cows) provided similar reliabilities in NZ Hol bulls, when compared with the single breed reference with crossbred NZ cows. For Jer validation bulls, the RS that included Jer cows and bulls and crossbred cows from NZ and Jer cows from AU was marginally better than the all-breed, all-country RS. In terms of reliability, the advantage of the HD SNP chip was small but captured more of the genomic variance than the 50K, particularly for Hol. The expectation maximization Bayes R GP method was slightly (up to 3 percentage points) better than genomic BLUP. We conclude that GP of milk production traits in NZ bulls improves by up to 7 percentage points in reliability by expanding the NZ reference population to include AU cows.
Collapse
Affiliation(s)
- M Haile-Mariam
- Agriculture Victoria, Department of Jobs, Precincts and Regions, Bundoora, VIC 3083, Australia.
| | - I M MacLeod
- Agriculture Victoria, Department of Jobs, Precincts and Regions, Bundoora, VIC 3083, Australia
| | - S Bolormaa
- Agriculture Victoria, Department of Jobs, Precincts and Regions, Bundoora, VIC 3083, Australia
| | | | | | - G de Jong
- CRV, 6800 AL Arnhem, the Netherlands
| | - H D Daetwyler
- Agriculture Victoria, Department of Jobs, Precincts and Regions, Bundoora, VIC 3083, Australia; School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
| | - J E Pryce
- Agriculture Victoria, Department of Jobs, Precincts and Regions, Bundoora, VIC 3083, Australia; School of Applied Systems Biology, La Trobe University, Bundoora, VIC 3083, Australia
| |
Collapse
|
7
|
Doekes HP, Veerkamp RF, Bijma P, Hiemstra SJ, Windig JJ. Trends in genome-wide and region-specific genetic diversity in the Dutch-Flemish Holstein-Friesian breeding program from 1986 to 2015. Genet Sel Evol 2018; 50:15. [PMID: 29642838 PMCID: PMC5896142 DOI: 10.1186/s12711-018-0385-y] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2017] [Accepted: 03/27/2018] [Indexed: 12/19/2022] Open
Abstract
Background In recent decades, Holstein–Friesian (HF) selection schemes have undergone profound changes, including the introduction of optimal contribution selection (OCS; around 2000), a major shift in breeding goal composition (around 2000) and the implementation of genomic selection (GS; around 2010). These changes are expected to have influenced genetic diversity trends. Our aim was to evaluate genome-wide and region-specific diversity in HF artificial insemination (AI) bulls in the Dutch-Flemish breeding program from 1986 to 2015. Methods Pedigree and genotype data (~ 75.5 k) of 6280 AI-bulls were used to estimate rates of genome-wide inbreeding and kinship and corresponding effective population sizes. Region-specific inbreeding trends were evaluated using regions of homozygosity (ROH). Changes in observed allele frequencies were compared to those expected under pure drift to identify putative regions under selection. We also investigated the direction of changes in allele frequency over time. Results Effective population size estimates for the 1986–2015 period ranged from 69 to 102. Two major breakpoints were observed in genome-wide inbreeding and kinship trends. Around 2000, inbreeding and kinship levels temporarily dropped. From 2010 onwards, they steeply increased, with pedigree-based, ROH-based and marker-based inbreeding rates as high as 1.8, 2.1 and 2.8% per generation, respectively. Accumulation of inbreeding varied substantially across the genome. A considerable fraction of markers showed changes in allele frequency that were greater than expected under pure drift. Putative selected regions harboured many quantitative trait loci (QTL) associated to a wide range of traits. In consecutive 5-year periods, allele frequencies changed more often in the same direction than in opposite directions, except when comparing the 1996–2000 and 2001–2005 periods. Conclusions Genome-wide and region-specific diversity trends reflect major changes in the Dutch-Flemish HF breeding program. Introduction of OCS and the shift in breeding goal were followed by a drop in inbreeding and kinship and a shift in the direction of changes in allele frequency. After introduction of GS, rates of inbreeding and kinship increased substantially while allele frequencies continued to change in the same direction as before GS. These results provide insight in the effect of breeding practices on genomic diversity and emphasize the need for efficient management of genetic diversity in GS schemes. Electronic supplementary material The online version of this article (10.1186/s12711-018-0385-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Harmen P Doekes
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands. .,Centre for Genetic Resources the Netherlands, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands.
| | - Roel F Veerkamp
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Piter Bijma
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Sipke J Hiemstra
- Centre for Genetic Resources the Netherlands, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Jack J Windig
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands.,Centre for Genetic Resources the Netherlands, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| |
Collapse
|
8
|
Blant A, Kwong M, Szpiech ZA, Pemberton TJ. Weighted likelihood inference of genomic autozygosity patterns in dense genotype data. BMC Genomics 2017; 18:928. [PMID: 29191164 PMCID: PMC5709839 DOI: 10.1186/s12864-017-4312-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2017] [Accepted: 11/16/2017] [Indexed: 12/14/2022] Open
Abstract
Background Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns. Methods We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data. Results Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies. Conclusions This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Alexandra Blant
- Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB, Canada
| | - Michelle Kwong
- Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB, Canada
| | - Zachary A Szpiech
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Trevor J Pemberton
- Department of Biochemistry and Medical Genetics, University of Manitoba, Winnipeg, MB, Canada.
| |
Collapse
|
9
|
Signer-Hasler H, Burren A, Neuditschko M, Frischknecht M, Garrick D, Stricker C, Gredler B, Bapst B, Flury C. Population structure and genomic inbreeding in nine Swiss dairy cattle populations. Genet Sel Evol 2017; 49:83. [PMID: 29115934 PMCID: PMC5674839 DOI: 10.1186/s12711-017-0358-6] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2016] [Accepted: 10/26/2017] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF). RESULTS The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26. CONCLUSIONS The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.
Collapse
Affiliation(s)
- Heidi Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | - Alexander Burren
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | | | - Mirjam Frischknecht
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
- Qualitas AG, Zug, Switzerland
| | | | | | | | | | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| |
Collapse
|
10
|
Howard JT, Pryce JE, Baes C, Maltecca C. Invited review: Inbreeding in the genomics era: Inbreeding, inbreeding depression, and management of genomic variability. J Dairy Sci 2017; 100:6009-6024. [DOI: 10.3168/jds.2017-12787] [Citation(s) in RCA: 84] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2017] [Accepted: 04/25/2017] [Indexed: 11/19/2022]
|
11
|
Brito LF, Kijas JW, Ventura RV, Sargolzaei M, Porto-Neto LR, Cánovas A, Feng Z, Jafarikia M, Schenkel FS. Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers. BMC Genomics 2017; 18:229. [PMID: 28288562 PMCID: PMC5348779 DOI: 10.1186/s12864-017-3610-0] [Citation(s) in RCA: 91] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Accepted: 03/07/2017] [Indexed: 01/08/2023] Open
Abstract
Background The detection of signatures of selection has the potential to elucidate the identities of genes and mutations associated with phenotypic traits important for livestock species. It is also very relevant to investigate the levels of genetic diversity of a population, as genetic diversity represents the raw material essential for breeding and has practical implications for implementation of genomic selection. A total of 1151 animals from nine goat populations selected for different breeding goals and genotyped with the Illumina Goat 50K single nucleotide polymorphisms (SNP) Beadchip were included in this investigation. Results The proportion of polymorphic SNPs ranged from 0.902 (Nubian) to 0.995 (Rangeland). The overall mean HO and HE was 0.374 ± 0.021 and 0.369 ± 0.023, respectively. The average pairwise genetic distance (D) ranged from 0.263 (Toggenburg) to 0.323 (Rangeland). The overall average for the inbreeding measures FEH, FVR, FLEUT, FROH and FPED was 0.129, −0.012, −0.010, 0.038 and 0.030, respectively. Several regions located on 19 chromosomes were potentially under selection in at least one of the goat breeds. The genomic population tree constructed using all SNPs differentiated breeds based on selection purpose, while genomic population tree built using only SNPs in the most significant region showed a great differentiation between LaMancha and the other breeds. We hypothesized that this region is related to ear morphogenesis. Furthermore, we identified genes potentially related to reproduction traits, adult body mass, efficiency of food conversion, abdominal fat deposition, conformation traits, liver fat metabolism, milk fatty acids, somatic cells score, milk protein, thermo-tolerance and ear morphogenesis. Conclusions In general, moderate to high levels of genetic variability were observed for all the breeds and a characterization of runs of homozygosity gave insights into the breeds’ development history. The information reported here will be useful for the implementation of genomic selection and other genomic studies in goats. We also identified various genome regions under positive selection using smoothed FST and hapFLK statistics and suggested genes, which are potentially under selection. These results can now provide a foundation to formulate biological hypotheses related to selection processes in goats. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3610-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Luiz F Brito
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.
| | - James W Kijas
- CSIRO Agriculture & Food, Brisbane, Queensland, Australia
| | - Ricardo V Ventura
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,Beef Improvement Opportunities, Guelph, Ontario, Canada
| | - Mehdi Sargolzaei
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,The Semex Alliance, Guelph, Ontario, Canada
| | | | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada
| | - Zeny Feng
- Department of Mathematics and Statistics, 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
| | - Flávio S Schenkel
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada
| |
Collapse
|
12
|
Howard JT, Tiezzi F, Huang Y, Gray KA, Maltecca C. Characterization and management of long runs of homozygosity in parental nucleus lines and their associated crossbred progeny. Genet Sel Evol 2016; 48:91. [PMID: 27884108 PMCID: PMC5123398 DOI: 10.1186/s12711-016-0269-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2016] [Accepted: 11/10/2016] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND In nucleus populations, regions of the genome that have a high frequency of runs of homozygosity (ROH) occur and are associated with a reduction in genetic diversity, as well as adverse effects on fitness. It is currently unclear whether, and to what extent, ROH stretches persist in the crossbred genome and how genomic management in the nucleus population might impact low diversity regions and its implications on the crossbred genome. METHODS We calculated a ROH statistic based on lengths of 5 (ROH5) or 10 (ROH10) Mb across the genome for genotyped Landrace (LA), Large White (LW) and Duroc (DU) dams. We simulated crossbred dam (LA × LW) and market [DU × (LA × LW)] animal genotypes based on observed parental genotypes and the ROH frequency was tabulated. We conducted a simulation using observed genotypes to determine the impact of minimizing parental relationships on multiple diversity metrics within nucleus herds, i.e. pedigree-(A), SNP-by-SNP relationship matrix or ROH relationship matrix. Genome-wide metrics included, pedigree inbreeding, heterozygosity and proportion of the genome in ROH of at least 5 Mb. Lastly, the genome was split into bins of increasing ROH5 frequency and, within each bin, heterozygosity, ROH5 and length (Mb) of ROH were evaluated. RESULTS We detected regions showing high frequencies of either ROH5 and/or ROH10 across both LW and LA on SSC1, SSC4, and SSC14, and across all breeds on SSC9. Long haplotypes were shared across parental breeds and thus, regions of ROH persisted in crossbred animals. Averaged across replicates and breeds, progeny had higher levels of heterozygosity (0.0056 ± 0.002%) and lower proportion of the genome in a ROH of at least 5 Mb (-0.015 ± 0.003%) than their parental genomes when genomic relationships were constrained, while pedigree relationships resulted in negligible differences at the genomic level. Across all breeds, only genomic data was able to target low diversity regions. CONCLUSIONS We show that long stretches of ROH present in the parents persist in crossbred animals. Furthermore, compared to using pedigree relationships, using genomic information to constrain parental relationships resulted in maintaining more genetic diversity and more effectively targeted low diversity regions.
Collapse
Affiliation(s)
- Jeremy T Howard
- Department of Animal Science, North Carolina State University, Raleigh, NC, 27695-7627, USA.
| | - Francesco Tiezzi
- Department of Animal Science, North Carolina State University, Raleigh, NC, 27695-7627, USA
| | - Yijian Huang
- Smithfield Premium Genetics, Rose Hill, NC, 28458, USA
| | - Kent A Gray
- Smithfield Premium Genetics, Rose Hill, NC, 28458, USA
| | - Christian Maltecca
- Department of Animal Science, North Carolina State University, Raleigh, NC, 27695-7627, USA.,Genetics Program, North Carolina State University, Raleigh, NC, 27695-7627, USA
| |
Collapse
|
13
|
Gómez-Romano F, Villanueva B, Fernández J, Woolliams JA, Pong-Wong R. The use of genomic coancestry matrices in the optimisation of contributions to maintain genetic diversity at specific regions of the genome. Genet Sel Evol 2016; 48:2. [PMID: 26763889 PMCID: PMC5518098 DOI: 10.1186/s12711-015-0172-y] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Accepted: 11/19/2015] [Indexed: 11/10/2022] Open
Abstract
Background Optimal contribution methods have proved to be very efficient for controlling the rates at which coancestry and inbreeding increase and therefore, for maintaining genetic diversity. These methods have usually relied on pedigree information for estimating genetic relationships between animals. However, with the large amount of genomic information now available such as high-density single nucleotide polymorphism (SNP) chips that contain thousands of SNPs, it becomes possible to calculate more accurate estimates of relationships and to target specific regions in the genome where there is a particular interest in maximising genetic diversity. The objective of this study was to investigate the effectiveness of using genomic coancestry matrices for: (1) minimising the loss of genetic variability at specific genomic regions while restricting the overall loss in the rest of the genome; or (2) maximising the overall genetic diversity while restricting the loss of diversity at specific genomic regions. Results Our study shows that the use of genomic coancestry was very successful at minimising the loss of diversity and outperformed the use of pedigree-based coancestry (genetic diversity even increased in some scenarios). The results also show that genomic information allows a targeted optimisation to maintain diversity at specific genomic regions, whether they are linked or not. The level of variability maintained increased when the targeted regions were closely linked. However, such targeted management leads to an important loss of diversity in the rest of the genome and, thus, it is necessary to take further actions to constrain this loss. Optimal contribution methods also proved to be effective at restricting the loss of diversity in the rest of the genome, although the resulting rate of coancestry was higher than the constraint imposed. Conclusions The use of genomic matrices when optimising contributions permits the control of genetic diversity and inbreeding at specific regions of the genome through the minimisation of partial genomic coancestry matrices. The formula used to predict coancestry in the next generation produces biased results and therefore it is necessary to refine the theory of genetic contributions when genomic matrices are used to optimise contributions.
Collapse
Affiliation(s)
| | | | - Jesús Fernández
- Departamento de Mejora Genética Animal, INIA, Madrid, Spain.
| | - John A Woolliams
- The Roslin Institute and the R(D)SVS, University of Edinburgh, Edinburgh, UK
| | - Ricardo Pong-Wong
- The Roslin Institute and the R(D)SVS, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
14
|
Moore SG, Pryce JE, Hayes BJ, Chamberlain AJ, Kemper KE, Berry DP, McCabe M, Cormican P, Lonergan P, Fair T, Butler ST. Differentially Expressed Genes in Endometrium and Corpus Luteum of Holstein Cows Selected for High and Low Fertility Are Enriched for Sequence Variants Associated with Fertility1. Biol Reprod 2016; 94:19. [DOI: 10.1095/biolreprod.115.132951] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2015] [Accepted: 11/24/2015] [Indexed: 11/01/2022] Open
|
15
|
Howard JT, Haile-Mariam M, Pryce JE, Maltecca C. Investigation of regions impacting inbreeding depression and their association with the additive genetic effect for United States and Australia Jersey dairy cattle. BMC Genomics 2015; 16:813. [PMID: 26481110 PMCID: PMC4612420 DOI: 10.1186/s12864-015-2001-7] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2015] [Accepted: 10/03/2015] [Indexed: 01/25/2023] Open
Abstract
Background Variation in environment, management practices, nutrition or selection objectives has led to a variety of different choices being made in the use of genetic material between countries. Differences in genome-level homozygosity between countries may give rise to regions that result in inbreeding depression to differ. The objective of this study was to characterize regions that have an impact on a runs of homozygosity (ROH) metric and estimate their association with the additive genetic effect of milk (MY), fat (FY) and protein yield (PY) and calving interval (CI) using Australia (AU) and United States (US) Jersey cows. Methods Genotyped cows with phenotypes on MY, FY and PY (n = 6751 US; n = 3974 AU) and CI (n = 5816 US; n = 3905 AU) were used in a two-stage analysis. A ROH statistic (ROH4Mb), which counts the frequency of a SNP being in a ROH of at least 4 Mb was calculated across the genome. In the first stage, residuals were obtained from a model that accounted for the portion explained by the estimated breeding value. In the second stage, these residuals were regressed on ROH4Mb using a single marker regression model and a gradient boosted machine (GBM) algorithm. The relationship between the additive and ROH4Mb of a region was characterized based on the (co)variance of 500 kb estimated genomic breeding values derived from a Bayesian LASSO analysis. Phenotypes to determine ROH4Mb and additive effects were residuals from the two-stage approach and yield deviations, respectively. Results Associations between yield traits and ROH4Mb were found for regions on BTA13, BTA23 and BTA25 for the US population and BTA3, BTA7, BTA17 for the AU population. Only one association (BTA7) was found for CI and ROH4Mb for the US population. Multiple potential epistatic interactions were characterized based on the GBM analysis. Lastly, the covariance sign between ROH4Mb and additive SNP effect of a region was heterogeneous across the genome. Conclusion We identified multiple genomic regions associated with ROH4Mb in US and AU Jersey females. The covariance of regions impacting ROH4Mb and the additive genetic effect were positive and negative, which provides evidence that the homozygosity effect is location dependent. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2001-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Jeremy T Howard
- Department of Animal Science and Genetics Program, North Carolina State University, Raleigh, NC, 27695-7627, USA.
| | - Mekonnen Haile-Mariam
- Department of Economic Development, Jobs, Transport and Resources and Dairy Futures Cooperative Research Centre, 5 Ring Road, Bundoora, VIC, 3083, Australia.
| | - Jennie E Pryce
- Department of Economic Development, Jobs, Transport and Resources and Dairy Futures Cooperative Research Centre, 5 Ring Road, Bundoora, VIC, 3083, Australia. .,La Trobe University, Bundoora, VIC, 3086, Australia.
| | - Christian Maltecca
- Department of Animal Science and Genetics Program, North Carolina State University, Raleigh, NC, 27695-7627, USA.
| |
Collapse
|
16
|
Kim ES, Sonstegard TS, Van Tassell CP, Wiggans G, Rothschild MF. The Relationship between Runs of Homozygosity and Inbreeding in Jersey Cattle under Selection. PLoS One 2015; 10:e0129967. [PMID: 26154171 PMCID: PMC4496098 DOI: 10.1371/journal.pone.0129967] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2014] [Accepted: 05/14/2015] [Indexed: 11/19/2022] Open
Abstract
Inbreeding is often an inevitable outcome of strong directional artificial selection but on average it reduces population fitness with increased frequency of recessive deleterious alleles. Runs of homozygosity (ROH) representing genomic autozygosity that occur from mating between selected and genomically related individuals may be able to reveal the regions affecting fitness. To examine the influence of genomic autozygosity on fitness, we used a genome-wide association test to evaluate potential negative correlations between ROH and daughter pregnancy rate (DPR) or somatic cell score (SCS) in US Jersey cattle. In addition, relationships between changes of local ROH and inbreeding coefficients (F) were assessed to locate genomic regions with increased inbreeding. Despite finding some decreases in fertility associated with incremental increases in F, most emerging local ROH were not significantly associated with DPR or SCS. Furthermore, the analyses of ROH could be approximated with the most frequent haplotype(s), including the associations of ROH and F or traits. The analysis of the most frequent haplotype revealed that associations of ROH and fertility could be accounted for by the additive genetic effect on the trait. Thus, we suggest that a change of autozygosity is more likely to demonstrate footprints of selected haplotypes for production rather than highlight the possible increased local autozygosity of a recessive detrimental allele resulting from the mating between closely related animals in Jersey cattle.
Collapse
Affiliation(s)
- Eui-Soo Kim
- Animal Genomics & Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, United States Department of Agriculture, Beltsville, Maryland, United States of America
- Department of Animal Science, Iowa State University, Ames, Iowa, United States of America
| | - Tad S. Sonstegard
- Animal Genomics & Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, United States Department of Agriculture, Beltsville, Maryland, United States of America
| | - Curtis P. Van Tassell
- Animal Genomics & Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, United States Department of Agriculture, Beltsville, Maryland, United States of America
| | - George Wiggans
- Animal Genomics & Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, United States Department of Agriculture, Beltsville, Maryland, United States of America
| | - Max F. Rothschild
- Department of Animal Science, Iowa State University, Ames, Iowa, United States of America
- * E-mail:
| |
Collapse
|