1
|
Caballero A, Fernández A, Villanueva B, Toro MA. A comparison of marker-based estimators of inbreeding and inbreeding depression. Genet Sel Evol 2022; 54:82. [PMID: 36575379 PMCID: PMC9793638 DOI: 10.1186/s12711-022-00772-0] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Accepted: 12/14/2022] [Indexed: 12/28/2022] Open
Abstract
BACKGROUND The availability of genome-wide marker data allows estimation of inbreeding coefficients (F, the probability of identity-by-descent, IBD) and, in turn, estimation of the rate of inbreeding depression (ΔID). We investigated, by computer simulations, the accuracy of the most popular estimators of inbreeding based on molecular markers when computing F and ΔID in populations under random mating, equalization of parental contributions, and artificially selected populations. We assessed estimators described by Li and Horvitz (FLH1 and FLH2), VanRaden (FVR1 and FVR2), Yang and colleagues (FYA1 and FYA2), marker homozygosity (FHOM), runs of homozygosity (FROH) and estimates based on pedigree (FPED) in comparison with estimates obtained from IBD measures (FIBD). RESULTS If the allele frequencies of a base population taken as a reference for the computation of inbreeding are known, all estimators based on marker allele frequencies are highly correlated with FIBD and provide accurate estimates of the mean ΔID. If base population allele frequencies are unknown and current frequencies are used in the estimations, the largest correlation with FIBD is generally obtained by FLH1 and the best estimator of ΔID is FYA2. The estimators FVR2 and FLH2 have the poorest performance in most scenarios. The assumption that base population allele frequencies are equal to 0.5 results in very biased estimates of the average inbreeding coefficient but they are highly correlated with FIBD and give relatively good estimates of ΔID. Estimates obtained directly from marker homozygosity (FHOM) substantially overestimated ΔID. Estimates based on runs of homozygosity (FROH) provide accurate estimates of inbreeding and ΔID. Finally, estimates based on pedigree (FPED) show a lower correlation with FIBD than molecular estimators but provide rather accurate estimates of ΔID. An analysis of data from a pig population supports the main findings of the simulations. CONCLUSIONS When base population allele frequencies are known, all marker-allele frequency-based estimators of inbreeding coefficients generally show a high correlation with FIBD and provide good estimates of ΔID. When base population allele frequencies are unknown, FLH1 is the marker frequency-based estimator that is most correlated with FIBD, and FYA2 provides the most accurate estimates of ΔID. Estimates from FROH are also very precise in most scenarios. The estimators FVR2 and FLH2 have the poorest performances.
Collapse
Affiliation(s)
- Armando Caballero
- grid.6312.60000 0001 2097 6738Centro de Investigación Mariña, Universidade de Vigo, Facultade de Bioloxía, 36310 Vigo, Spain
| | - Almudena Fernández
- Departamento de Mejora Genética Animal, INIA-CSIC, Ctra. de La Coruña, Km 7.5, 28040 Madrid, Spain
| | - Beatriz Villanueva
- Departamento de Mejora Genética Animal, INIA-CSIC, Ctra. de La Coruña, Km 7.5, 28040 Madrid, Spain
| | - Miguel A. Toro
- grid.5690.a0000 0001 2151 2978Departamento de Producción Agraria, ETSI Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, 28040 Madrid, Spain
| |
Collapse
|
2
|
Villanueva B, Fernández A, Saura M, Caballero A, Fernández J, Morales-González E, Toro MA, Pong-Wong R. The value of genomic relationship matrices to estimate levels of inbreeding. Genet Sel Evol 2021; 53:42. [PMID: 33933002 PMCID: PMC8088726 DOI: 10.1186/s12711-021-00635-0] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 04/19/2021] [Indexed: 12/12/2022] Open
Abstract
BACKGROUND Genomic relationship matrices are used to obtain genomic inbreeding coefficients. However, there are several methodologies to compute these matrices and there is still an unresolved debate on which one provides the best estimate of inbreeding. In this study, we investigated measures of inbreeding obtained from five genomic matrices, including the Nejati-Javaremi allelic relationship matrix (FNEJ), the Li and Horvitz matrix based on excess of homozygosity (FL&H), and the VanRaden (methods 1, FVR1, and 2, FVR2) and Yang (FYAN) genomic relationship matrices. We derived expectations for each inbreeding coefficient, assuming a single locus model, and used these expectations to explain the patterns of the coefficients that were computed from thousands of single nucleotide polymorphism genotypes in a population of Iberian pigs. RESULTS Except for FNEJ, the evaluated measures of inbreeding do not match with the original definitions of inbreeding coefficient of Wright (correlation) or Malécot (probability). When inbreeding coefficients are interpreted as indicators of variability (heterozygosity) that was gained or lost relative to a base population, both FNEJ and FL&H led to sensible results but this was not the case for FVR1, FVR2 and FYAN. When variability has increased relative to the base, FVR1, FVR2 and FYAN can indicate that it decreased. In fact, based on FYAN, variability is not expected to increase. When variability has decreased, FVR1 and FVR2 can indicate that it has increased. Finally, these three coefficients can indicate that more variability than that present in the base population can be lost, which is also unreasonable. The patterns for these coefficients observed in the pig population were very different, following the derived expectations. As a consequence, the rate of inbreeding depression estimated based on these inbreeding coefficients differed not only in magnitude but also in sign. CONCLUSIONS Genomic inbreeding coefficients obtained from the diagonal elements of genomic matrices can lead to inconsistent results in terms of gain and loss of genetic variability and inbreeding depression estimates, and thus to misleading interpretations. Although these matrices have proven to be very efficient in increasing the accuracy of genomic predictions, they do not always provide a useful measure of inbreeding.
Collapse
Affiliation(s)
- Beatriz Villanueva
- Departamento de Mejora Genética Animal, INIA, Ctra. de La Coruña, km 7.5, 28040 Madrid, Spain
| | - Almudena Fernández
- Departamento de Mejora Genética Animal, INIA, Ctra. de La Coruña, km 7.5, 28040 Madrid, Spain
| | - María Saura
- Departamento de Mejora Genética Animal, INIA, Ctra. de La Coruña, km 7.5, 28040 Madrid, Spain
| | - Armando Caballero
- Centro de Investigación Mariña, Universidade de Vigo, Departamento de Bioquímica, Genética E Inmunología, Campus de Vigo, 36310 Vigo, Spain
| | - Jesús Fernández
- Departamento de Mejora Genética Animal, INIA, Ctra. de La Coruña, km 7.5, 28040 Madrid, Spain
| | | | - Miguel A. Toro
- Departamento de Producción Agraria, ETSI Agrónomos, Universidad Politécnica de Madrid, 28040 Madrid, Spain
| | - Ricardo Pong-Wong
- Genetics and Genomics, The Roslin Institute and the R(D)SVS, University of Edinburgh, Easter Bush, Midlothian, EH25 9RG UK
| |
Collapse
|
3
|
Martikainen K, Koivula M, Uimari P. Identification of runs of homozygosity affecting female fertility and milk production traits in Finnish Ayrshire cattle. Sci Rep 2020; 10:3804. [PMID: 32123255 PMCID: PMC7052207 DOI: 10.1038/s41598-020-60830-9] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 02/17/2020] [Indexed: 11/09/2022] Open
Abstract
Inbreeding gives rise to continuous lengths of homozygous genotypes called runs of homozygosity (ROH) that occur when identical haplotypes are inherited from both parents. ROHs are enriched for deleterious recessive alleles and can therefore be linked to inbreeding depression, defined as decreased phenotypic performance of the animals. However, not all ROHs within a region are expected to have harmful effects on the trait of interest. We aimed to identify ROHs that unfavourably affect female fertility and milk production traits in the Finnish Ayrshire population. The estimated effect of ROHs with the highest statistical significance varied between parities from 9 to 17 days longer intervals from calving to first insemination, from 13 to 38 days longer intervals from first to last insemination and from 0.3 to 1.0 more insemination per conception. Similarly, for milk production traits ROHs were associated with a reduction of 208 kg for milk yield, 7 kg for protein yield and 16 kg for fat yield. We also found regions where ROHs displayed unfavourable effects across multiple traits. Our findings can be exploited for more efficient control of inbreeding depression, for example by minimizing the occurrence of unfavourable haplotypes as homozygous state in breeding programmes.
Collapse
Affiliation(s)
- K Martikainen
- Department of Agricultural Sciences, University of Helsinki, P.O. Box 28, FI-00014, Helsinki, Finland.
| | - M Koivula
- Natural Resources Institute Finland (Luke), Green Technology, FI-31600, Jokioinen, Finland
| | - P Uimari
- Department of Agricultural Sciences, University of Helsinki, P.O. Box 28, FI-00014, Helsinki, Finland
| |
Collapse
|
4
|
Martikainen K, Sironen A, Uimari P. Estimation of intrachromosomal inbreeding depression on female fertility using runs of homozygosity in Finnish Ayrshire cattle. J Dairy Sci 2018; 101:11097-11107. [DOI: 10.3168/jds.2018-14805] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2018] [Accepted: 08/22/2018] [Indexed: 12/31/2022]
|
5
|
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
|
6
|
Mastrangelo S, Tolone M, Sardina MT, Sottile G, Sutera AM, Di Gerlando R, Portolano B. Genome-wide scan for runs of homozygosity identifies potential candidate genes associated with local adaptation in Valle del Belice sheep. Genet Sel Evol 2017; 49:84. [PMID: 29137622 PMCID: PMC5684758 DOI: 10.1186/s12711-017-0360-z] [Citation(s) in RCA: 79] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2017] [Accepted: 11/07/2017] [Indexed: 11/23/2022] Open
Abstract
Background Because very large numbers of single nucleotide polymorphisms (SNPs) are now available throughout the genome, they are particularly suitable for the detection of genomic regions where a reduction in heterozygosity has occurred and they offer new opportunities to improve the accuracy of inbreeding (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$F$$\end{document}F) estimates. Runs of homozygosity (ROH) are contiguous lengths of homozygous segments of the genome where the two haplotypes inherited from the parents are identical. Here, we investigated the occurrence and distribution of ROH using a medium-dense SNP panel to characterize autozygosity in 516 Valle del Belice sheep and to identify the genomic regions with high ROH frequencies. Results We identified 11,629 ROH and all individuals displayed at least one ROH longer than 1 Mb. The mean value of \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$F$$\end{document}F estimated from ROH longer than1 Mb was 0.084 ± 0.061. ROH that were shorter than 10 Mb predominated. The highest and lowest coverages of Ovis aries chromosomes (OAR) by ROH were on OAR24 and OAR1, respectively. The number of ROH per chromosome length displayed a specific pattern, with higher values for the first three chromosomes. Both number of ROH and length of the genome covered by ROH varied considerably between animals. Two hundred and thirty-nine SNPs were considered as candidate markers that may be under directional selection and we identified 107 potential candidate genes. Six genomic regions located on six chromosomes, corresponding to ROH islands, are presented as hotspots of autozygosity, which frequently coincided with regions of medium recombination rate. According to the KEGG database, most of these genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. A genome scan revealed the presence of ROH islands in genomic regions that harbor candidate genes for selection in response to environmental stress and which underlie local adaptation. Conclusions These results suggest that natural selection has, at least partially, a role in shaping the genome of Valle del Belice sheep and that ROH in the ovine genome may help to detect genomic regions involved in the determinism of traits under selection. Electronic supplementary material The online version of this article (10.1186/s12711-017-0360-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università degli Studi di Palermo, 90128, Palermo, Italy.
| | - Marco Tolone
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università degli Studi di Palermo, 90128, Palermo, Italy
| | - Maria T Sardina
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università degli Studi di Palermo, 90128, Palermo, Italy
| | - Gianluca Sottile
- Dipartimento di Scienze Economiche, Aziendali e Statistiche, Università degli Studi di Palermo, 90128, Palermo, Italy
| | - Anna M Sutera
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università degli Studi di Palermo, 90128, Palermo, Italy
| | - Rosalia Di Gerlando
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università degli Studi di Palermo, 90128, Palermo, Italy
| | - Baldassare Portolano
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università degli Studi di Palermo, 90128, Palermo, Italy
| |
Collapse
|
7
|
Martikainen K, Tyrisevä AM, Matilainen K, Pösö J, Uimari P. Estimation of inbreeding depression on female fertility in the Finnish Ayrshire population. J Anim Breed Genet 2017; 134:383-392. [PMID: 28748554 DOI: 10.1111/jbg.12285] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Accepted: 06/14/2017] [Indexed: 01/01/2023]
Abstract
Single nucleotide polymorphism (SNP) data enable the estimation of inbreeding at the genome level. In this study, we estimated inbreeding levels for 19,075 Finnish Ayrshire cows genotyped with a low-density SNP panel (8K). The genotypes were imputed to 50K density, and after quality control, 39,144 SNPs remained for the analysis. Inbreeding coefficients were estimated for each animal based on the percentage of homozygous SNPs (FPH ), runs of homozygosity (FROH ) and pedigree (FPED ). Phenotypic records were available for 13,712 animals including non-return rate (NRR), number of inseminations (AIS) and interval from first to last insemination (IFL) for heifers and up to three parities for cows, as well as interval from calving to first insemination (ICF) for cows. Average FPED was 0.02, FROH 0.06 and FPH 0.63. A correlation of 0.71 was found between FPED and FROH , 0.66 between FPED and FPH and 0.94 between FROH and FPH . Pedigree-based inbreeding coefficients did not show inbreeding depression in any of the traits. However, when FROH or FPH was used as a covariate, significant inbreeding depression was observed; a 10% increase in FROH was associated with 5 days longer IFL0 and IFL1, 2 weeks longer IFL3 and 3 days longer ICF2 compared to non-inbred cows.
Collapse
Affiliation(s)
- K Martikainen
- Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland
| | - A M Tyrisevä
- Natural Resources Institute Finland (Luke), Jokioinen, Finland
| | - K Matilainen
- Natural Resources Institute Finland (Luke), Jokioinen, Finland
| | - J Pösö
- Finnish Animal Breeding Association, Vantaa, Finland
| | - P Uimari
- Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland
| |
Collapse
|