Identification of genomic regions associated with inbreeding depression in Holstein and Jersey dairy cattle
© Pryce et al.; licensee BioMed Central Ltd. 2014
Received: 2 May 2014
Accepted: 6 October 2014
Published: 18 November 2014
Inbreeding reduces the fitness of individuals by increasing the frequency of homozygous deleterious recessive alleles. Some insight into the genetic architecture of fitness, and other complex traits, can be gained by using single nucleotide polymorphism (SNP) data to identify regions of the genome which lead to reduction in performance when identical by descent (IBD). Here, we compared the effect of genome-wide and location-specific homozygosity on fertility and milk production traits in dairy cattle.
Genotype data from more than 43 000 SNPs were available for 8853 Holstein and 4138 Jersey dairy cows that were part of a much larger dataset that had pedigree records (338 696 Holstein and 64 049 Jersey animals). Measures of inbreeding were based on: (1) pedigree data; (2) genotypes to determine the realised proportion of the genome that is IBD; (3) the proportion of the total genome that is homozygous and (4) runs of homozygosity (ROH) which are stretches of the genome that are homozygous.
A 1% increase in inbreeding based either on pedigree or genomic data was associated with a decrease in milk, fat and protein yields of around 0.4 to 0.6% of the phenotypic mean, and an increase in calving interval (i.e. a deterioration in fertility) of 0.02 to 0.05% of the phenotypic mean. A genome-wide association study using ROH of more than 50 SNPs revealed genomic regions that resulted in depression of up to 12.5 d and 260 L for calving interval and milk yield, respectively, when completely homozygous.
Genomic measures can be used instead of pedigree-based inbreeding to estimate inbreeding depression. Both the diagonal elements of the genomic relationship matrix and the proportion of homozygous SNPs can be used to measure inbreeding. Longer ROH (>3 Mb) were found to be associated with a reduction in milk yield and captured recent inbreeding independently and in addition to overall homozygosity. Inbreeding depression can be reduced by minimizing overall inbreeding but maybe also by avoiding the production of offspring that are homozygous for deleterious alleles at specific genomic regions that are associated with inbreeding depression.
Inbreeding depression is the reduction in fitness of offspring that result from the mating between individuals that share at least one common ancestor. Inbreeding reduces fitness by increasing the number of homozygous deleterious recessive alleles that affect traits related to fitness, such as survival, disease resistance, predation and birth weight in avian and mammalian populations . Some published examples include a dramatic effect on offspring survival and reproduction in both mice  and freshwater snails . In livestock and poultry species, preventing inbreeding is important because it is associated with a reduction in animal performance and consequently profitability - and with an increased frequency of genetic defects . Ideally, livestock breeding programs increase the frequency of favourable alleles and possibly their fixation, while minimizing homozygosity for deleterious recessive alleles. With the availability of genotyping data based on high-density SNP (single nucleotide polymorphism) chips, these objectives become possible.
Genetic markers can be used to estimate inbreeding as the realised proportion of the genome that is identical by descent (IBD), for example by calculating a genomic relationship matrix (GRM) among individuals -. The expected value of a diagonal element of the GRM, conditional on the pedigree of the animal, is 1 + F where F is the inbreeding coefficient of the animal. Unlike pedigree-based inbreeding, a genomic estimate of inbreeding does not suffer from lack of depth of pedigree data and pedigree errors and it measures realised inbreeding that can vary between animals that have the same pedigree, i.e. full sibs.
The diagonal elements of the GRM and the proportion of the genome where SNPs are homozygous are both examples of genomic measures that can be used to estimate the level of inbreeding across the genome. A limitation of both of these measures is that they do not distinguish between identity by state and identity by descent. One possibility to resolve this is to use runs of homozygosity (ROH), which are regions of the genome where the copies inherited from the two parents are IBD . ROH can be used to investigate genome-wide inbreeding and to identify the localization of specific regions of the genome that are IBD. Compared to the diagonal or the GRM or proportion of homozygous SNPs, long ROH are unlikely to have arisen by chance, and are more likely to be stretches of homologous chromosomes within the same individual that are IBD . In fact, the length of ROH inversely correlates to the distance in the pedigree where a common ancestor is present , which makes it a powerful method to detect inbreeding effects . ROH are shorter when the common ancestor that gave rise to the inbreeding arose many generations back (distant inbreeding), while they are longer when inbreeding occurred more recently . Inbreeding arising from a distant common ancestor should have less effect on fitness compared with inbreeding from a recent common relative because natural selection over long periods of time should act to purge deleterious alleles from the population . For example, in mice, the effect of inbreeding depression on litter size has been observed to be greater for inbreeding from recent compared with more distant common ancestors .
Dairy cattle are potentially a good model to investigate the effects of genome-wide and genome location-specific inbreeding depression, since phenotypes of traits associated with fitness (e.g., fertility) and quantitative performance traits (such as milk yield) are measured on a large scale. In addition, intense selection of males to be the sires of the next generation has made it almost impossible to find dairy cattle without multiple common ancestors , so there is substantial variation in the level of inbreeding among individuals in dairy populations.
Estimate inbreeding depression for milk production and fitness traits using pedigree-based inbreeding coefficients.
Compare these estimates to inbreeding depression estimated from genome-wide estimates that were derived from SNP data.
Investigate whether there are specific regions of the genome that are associated with inbreeding depression in fertility and milk production.
Pedigree-based inbreeding coefficients
Pedigree data that traced back to the 1950s were available on 2 070 219 Holstein and 293 588 Jersey cows . Animals with unknown parents were assigned to genetic groups based on their birth year, country of origin and sex, following standard Australian Dairy Herd Improvement Scheme (ADHIS; Melbourne, Australia) genetic evaluation procedures .
where C is the pedigree completeness, d is the number of generations and a i is the proportion of known ancestors in generation i. Selection of cows was restricted to animals born in or after 1994 to be able to compare the dataset with only pedigree information with the genotyped population. After editing, pedigree-based inbreeding coefficients were calculated using the algorithm of Meuwissen and Luo  for 338 696 Holstein and 64 049 Jersey individuals.
Genomic inbreeding coefficients
The genotyped cows that were included in this study were part of an initiative of the Dairy Future’s Cooperative Research Centre (Melbourne, Australia). Briefly, cows were selected to become part of the Australian genomic reference population on the basis of completeness of their phenotypic data. They were also selected to have sires that were genotyped and part of the national Australian genomic reference population. For each cow in the dataset, the validity of its sire was verified by checking for excessive “opposing homozygotes”, which is when an individual’s alleles at a given locus are homozygous, but for different alleles than its sire . The pedigree was corrected in cases for which a parent could be identified . Offspring that had more than 20 genotypes that were incompatible with those of their parents were removed. After quality control, genotypes were available for 8853 Holstein and 4138 Jersey females, which were genotyped with the Illumina BovineSNP50 BeadChip (Illumina, San Diego, CA; ) using a previously described quality control method for SNP editing . After editing, 45 753 and 43 737 SNPs remained for the Holstein and Jersey datasets, respectively. The SNPs were ordered by chromosome position using Bos taurus build UMD 3.1 (Center for Bioinformatics and Computational Biology, University of Maryland, MD). The average spacing between SNPs retained for this study was 58 kb. Homozygous SNPs were coded as 2 or 0 and heterozygous SNPs as 1.
- (1)The genomic relationship of an individual with itself relative to a base population was calculated as the diagonal of the genomic relationship matrix minus 1 (GRM_F) using ,:
where N is the number of SNPs, p m is the allele frequency of SNP m, and x m i s the genotype code at SNP m (0, 1 or 2). Allele frequencies p m were calculated separately for the Jersey and Holstein populations, which included genotyped males of the same breed, which amounted to a total of 12 649 and 5240 Holstein and Jersey individuals, respectively. Thus, GRM_F was relative to the current population, as represented by the sample of genotyped females and bulls that are part of the Australian genomic reference population.
The proportion of the genome that was homozygous (i.e. the proportion of total SNP genotypes that were either AA or BB).
The proportion of the genome that consists of runs of homozygosity (ROH_F) that were at least n SNPs in length, where n ranged from 5 to 100. If the ROH at SNP position i exceeded n, then ROH i was coded as 1 or else 0. ROH_F was the sum of ROH i divided by the total number of SNPs (nSNP):
Means and standard deviations (SD) of milk, fat and protein yields and calving interval
Calving interval (d)
where y ijklmn is the n th record of the phenotype (milk, fat or protein yield or calving interval), HYS i the i th herd-year-season of calving, parity j the j th parity, month k the k th month of calving, and b1 the regression coefficient on ageijklmn , which was the age at first calving, b2 is the regression coefficient on F ijklmn , which was the measure of inbreeding (pedigree, GRM_F, homozygosity, or ROH_F), perm l is the l th random permanent environmental effect to account for multiple records on each individual, and cow m was the random genetic effect for the m th animal, assumed to follow the distribution N(0, Aσ2), where A is the numerator relationship matrix, and e ijklmn is the random error associated with the observation. The quadratic effects of inbreeding were also tested in preliminary analyses but are not presented since none were significant. ASReml version 3.0 was used for the analysis . Finally, the genomic measures of inbreeding were tested after correcting for pedigree inbreeding by simultaneously fitting pedigree inbreeding and each genomic measure of inbreeding in turn.
This model was similar to that defined above, except that F was replaced by two new covariates: (1) the regression coefficient (b3) on SNP ijklmn , which was to correct for the additive effect of the presence of additional copies of the allele, and (2) the regression coefficient (b4) on ROH ijklmn , where, at each SNP position, ROH ijklmn was coded as 1 when a run of homozygosity of at least 50 SNPs was present at SNP ijklmn for that animal, and as 0 otherwise. A ROH was selected as having a statistically significant association with the trait analyzed if the regression coefficient had a p-value less than 0.001.
To identify possible candidate genes, QTL regions were defined based on SNP positions with significant associations (P <0.001) of ROH with milk yield or fertility. A list of annotated genes between the positions of the first and last SNP in the ROH or cluster of ROH was obtained from Ensembl BioMart MartView (http://asia.ensembl.org/biomart/martview/). QTL regions were also compared to QTL previously identified for the analysed trait through the Cattle QTLdb  (http://www.animalgenome.org/cgi-bin/QTLdb/BT/index).
Means and standard deviations in parentheses of measures of inbreeding in Holstein and Jersey populations
Year of birth3
Correlations between pedigree and genomic inbreeding coefficients in Holstein and Jersey populations
Effect on milk production traits and fertility in Holstein and Jersey populations of increasing inbreeding by 1%
Calving interval (d)
Calving interval (d)
For most other traits, the effects of increasing GRM_F were similar for the Holstein and Jersey populations (Table 4). However, as a percentage of the phenotypic mean, the average effects of GRM_F on milk, fat and protein yields were greater for the Jersey breed (0.57%) than for the Holstein breed (0.41%).
Regression coefficients of milk production traits and fertility on homozygosity in Holstein and Jersey breeds
Calving interval (d)
The association between ROH_F (the proportion of the genome that consists of ROH) and fertility was not significant for either breed (results not shown). In fact, for the Holstein breed, the effect of ROH_F on fertility was almost entirely captured by the overall homozygosity.
Genome-wide association studies
After correcting for the additive effect of SNPs, some regions of the genome contained stretches of ROH of 50 or more SNPs that were significantly associated with both milk yield and fertility for the Holstein and Jersey breeds. Significant ROH frequently formed clusters at overlapping positions on the genome, which may indicate that they were in linkage disequilibrium with the same QTL.
Genomic regions associated with a reduction in fertility and milk yield identified using runs of homozygosity
Number of ROH P <0.0013
Frequency of ROH4(%)
For calving interval in Holstein cows, seven clusters of ROH were identified on different chromosomes. These clusters comprised a total of 63 ROH that were significant, of which 52 had an unfavourable effect. The false discovery rate (FDR) of the 63 ROH was estimated at 72.5% using the method of Bolormaa et al. . The most probable real associations are likely to be the largest clusters; for example, there was a region on chromosome 9 in which 26 ROH had an unfavourable association with fertility in the Holstein population (Table 6). The size of this region was nearly 3 Mb (i.e. between 6.4 and 9.2 Mb) and included all the ROH (of 50 SNPs or more) present and led to an average increase in calving interval of 7.6 d (ranging from 6.7 to 8.4 d) in the Holstein population. For milk yield in the Holstein population, 117 ROH were significant (P <0.001), which corresponds to a FDR of 39%. Of these 117 ROH, 65 had an unfavourable association with milk yield and formed eight clusters. A group of 22 significant ROH (at P <0.001) was found on chromosome 16, which spanned a region between 64.9 and 66.2 Mb that contained 18 candidate genes (Table 6). On average, these ROH were associated with a reduction in milk yield of ~228 L/lactation. A region on chromosome 2 between 12.8 and 14.3 Mb that contained 24 ROH showed a significant (P <0.01) and positive (i.e. favourable) association with milk yield and an unfavourable association with calving interval. Although a more stringent significance level was necessary to detect this relationship, it is possible that there are many other regions of the genome with such opposing effects, since correlations between fertility and milk yield are generally unfavourable.
Since the Jersey population was smaller, there was less power to detect ROH that had an association with either milk yield or fertility. For this breed, only 44 and 17 ROH were significantly associated with milk yield and fertility, respectively, with FDR of 100% or more indicating that all these associations were probably identified by chance. However, 14 ROH on chromosome X had an effect on fertility in this breed (Table 6), with an FDR of only 6%, which suggests that this region is associated with inbreeding depression. Three regions on chromosomes 8, 20 and 24 were associated with milk yield in the Jersey breed, although none of these overlapped with those found for the Holstein breed. Detecting identical ROH with associations in two breeds is useful for validation, since it is unlikely that the same association is found in two breeds by chance. In our data, such cases were rare and only occurred when the significance threshold for the Jersey breed was relaxed to P <0.01, which is justified by its smaller population size. Of particular interest was a region found on chromosome 24 at 60 Mb that had an unfavourable association with fertility in the Holstein dataset (Table 6), which was validated in the Jersey dataset at P <0.01.
Correlations of 0.26 and 0.29 were found between inbreeding coefficients based on pedigree and GRM_F for the Jersey and Holstein datasets, respectively (Table 3). These values are lower than those observed for USA Holstein bulls (which ranged from 0.5 to 0.56; ) and for Australian Holstein bulls (which ranged from 0.67 to 0.87 for bulls with two to eight generations of recorded ancestry; ). There are three possible reasons why the correlations found in our study are lower than those reported in other studies: (1) individuals from sub-populations for which allele frequencies diverge from those of the entire population are estimated to have high inbreeding based on GRM_F; (2) pedigree completeness and (3) a possible bias introduced by the selection of the cows that were genotyped, since cows that had multiple lactation records were preferentially selected, which may have purged strongly deleterious effects out of this older population. In fact, the genotyped populations were on average 1.2 lactations older for the Holstein dataset and 0.6 lactations older for the Jersey dataset. So, purging may explain why pedigree inbreeding depression (Table 4) was on average 59% and 12.8% lower in the genotyped Holstein and Jersey populations, respectively, compared with the population that had only pedigree information. This comparison suggests that the effect of purging in the genotyped population could be greater than the value of 12.6% predicted by Gulisija and Crow .
Correlations of the proportion of homozygous SNPs with pedigree-based inbreeding coefficients were stronger than the equivalent correlations using GRM_F for both Holstein and Jersey populations. The lower correlations obtained with GRM_F are probably due to the effect of allele frequencies in the sub-population on the diagonal elements of the GRM, which results in artificially elevated “inbreeding coefficients” in minority sub-populations. To overcome this problem, several authors ,, have proposed calculating the GRM with allele frequencies fixed at 0.5. In fact, the resulting GRM is the same as the proportion of homozygous SNPs, since the correlation between these was equal to 0.99 . Therefore, using either the proportion of homozygous SNPs or the equivalent measure of the diagonal of the GRM calculated with an allele frequency of 0.5 is preferred for the estimation of a base population.
In accordance with other studies, the correlations between pedigree and genomic measures of inbreeding in cows were lower than the equivalent estimates in bulls ,, which implies that pedigree recording was worse for females than for males. This is expected for commercial dairy populations and is exacerbated in countries like Australia where calving occurs in unsupervised outside conditions and mismothering of calves is a problem. Furthermore, partial or incomplete pedigrees also reduce pedigree inbreeding estimates . Pedigree completeness was lower for the Jersey dataset than for the Holstein dataset, which probably explains why the correlation between proportion of homozygous SNPs and pedigree-based inbreeding coefficients was weaker for the Jersey than for the Holstein population. This is consistent with the fact that pedigree inbreeding tends to be lower for genotyped Jersey cows than for all Jersey and Holstein animals (Figure 1). However, inbreeding for all measures derived from genomic data was higher for the Jersey than the Holstein population. Furthermore, Jersey cows had on average a higher proportion of homozygous SNPs than Holstein cows (0.72 and 0.66, respectively) which is probably because (1) Jerseys were more inbred or (2) many of the SNPs on the genotyping panel were discovered using Holstein animals .
The correlation between ROH and pedigree-based inbreeding was strongest for ROH that consisted of 10 SNPs (0.58 Mb) (Figure 4a) and stronger than the correlation of any other genomic measure with pedigree-based inbreeding (Table 3). In the absence of pedigree data, the proportion of the genome that comprises ROH of 10 SNPs or more can be used to infer recent common ancestors ,, which makes ROH an attractive tool to monitor and evaluate inbreeding depression.
A 1% increase in inbreeding in the population with pedigree information only, or in the population with GRM_F resulted in a reduction in milk production traits of around 0.3% of the phenotypic mean and a deterioration in fertility for both the Holstein (Table 4) and the Jersey breeds. Comparable results have previously been published using pedigree-based inbreeding; a 1% increase in inbreeding has been estimated to be associated with a reduction in milk yield of around 20 to 30 L/lactation -,, and a lengthening of calving interval of up to 0.7 d ,. The small effect of inbreeding on calving interval may be due to the tendency of Australian farmers to cull cows that are expected to calve late, causing a downward bias in the effect on calving interval.
There have been relatively few studies on inbreeding depression using genomic measures of inbreeding. Recently, Bjelland et al.  estimated inbreeding depression using a measure similar to GRM_F and found a reduction in milk yield of 47 L/lactation and an increase of 1.06 d for days open (which is a trait very similar to calving interval) per % increase in inbreeding. The results presented in our study and those of Bjelland et al.  show that genomic estimates of inbreeding can be used instead of pedigree estimates to calculate the effects of inbreeding on performance and fitness traits. Under ideal conditions, the effect of pedigree inbreeding and GRM_F would be the same. However, errors in pedigree records are expected to reduce pedigree-based estimates of inbreeding depression. Estimates based on GRM_F should have a lower sampling error than those based on pedigree because it is not affected by incomplete pedigrees and it uses observed rather than expected IBD. However, this benefit had less effect on standard errors obtained in this study than the large difference in sample size used for the pedigree and genomic analyses. The magnitude of the difference in inbreeding depression estimated based on GRM_F and based on pedigree information for the genotyped cows implies that estimates of inbreeding depression on traits of economic importance may in fact be underestimated using pedigree information.
The proportion of homozygous SNPs was associated with inbreeding depression in all traits analysed, with the exception of fertility in the Jersey breed, which could be due to a lack of statistical power and bias in calving interval caused by culling. However, the association between homozygosity and milk production was not stronger than that between GRM_F and milk production (Tables 4 and 5). The allele frequencies of SNPs contain information on the probability that chromosome segments are IBD because segments that share a rare allele are more likely to be IBD than segments that share a common allele. Thus, including the base population in the calculation of inbreeding (i.e. GRM_F) can be advantageous. However, the correlation between pedigree-based inbreeding and SNP homozygosity was higher and as previously discussed, sub-populations that have allele frequencies that differ markedly from those of the base population appear to be more inbred than they actually are. Therefore, GRM_F has both advantages and disadvantages compared to SNP homozygosity as a measure of inbreeding. These advantages and disadvantages appear to approximately compensate each other, so while homozygosity may be a better indicator of pedigree-based inbreeding, it appears that GRM_F and homozygosity are equally good indicators of inbreeding depression.
Long ROH (>60 SNP or 3.5 Mb) were associated with a decrease in milk yield after correcting for average homozygosity, but short ROH were not. This could be due to selection eliminating deleterious mutations before inbreeding occurs if the inbreeding is due to an ancient common ancestor. Inbreeding arising from recent common ancestors has also been found to be associated with reproductive performance and body weight in mice  and with a large decrease in height in humans .
To account for recent and older inbreeding simultaneously in breeding programs that aim at reducing inbreeding, longer ROH (>50 SNPs or 3 Mb) should be explicitly considered in addition to other measures of genomic inbreeding. In this study, the Illumina BovineSNP50 BeadChip (Illumina, San Diego, CA; ) was used because of its popularity to estimate genomic breeding values . However, while this chip may be suitable for genomic evaluations, it has been reported to overestimate the number of ROH present compared to the Illumina HD panel . However, both SNP panels are equally effective at detecting segments that are longer than 4 Mb, which corresponds to approximately 70 SNPs . Thus, a negative aspect of our study could be that short ROH are not necessarily IBD. However, relying on longer ROH has the disadvantage that they have a lower frequency (Figure 2) and thus fail to identify small IBD regions. An area for future research could be to compare inbreeding effects of ROH obtained with high-density chips, or even whole-genome sequences with those obtained with lower density SNP panels.
Genome-wide association study
For the GWAS analysis, a ROH length of at least 50 SNPs (~3 Mb) was chosen because, in this case, the ROH frequency was large enough to detect statistically significant differences and the ROH length was sufficient to detect recent inbreeding. Other studies  have reported that using the Illumina BovineSNP50 BeadChip with ROH length set at 50 SNPs allowed the detection of fragments of more variable size (especially those less than 5 Mb) when compared to ROH of 50 SNPs. After correcting for the additive effect of the SNP at a given position, using a GWAS analysis, we still detected an effect of ROH on either calving interval or milk yield of up to 12.5 d and 260 L, respectively. This implies that there is considerable value in selecting for heterozygosity in regions where homozygosity has an unfavourable effect on valuable traits, such as fertility and milk yield. However, although our analysis identified several regions of interest, the associated FDR were high, which means that further validation is required.
The regions that were identified here were not close to the recessive lethal haplotypes reported by VanRaden et al.  and Fritz et al. , for which no homozygotes were observed. This is not surprising, since the homozygous regions detected in our study were present in the population and therefore, were unlikely to have a lethal effect.
However, the ROH that were identified as significant could be associated with genes that, if impaired, can lead to loss-of-function. In addition, when the presence of a ROH has a positive effect on a trait under selection, such as milk yield or fertility, it can be indicative of selection signatures .
For the Holstein population, a region on chromosome 16 contained 22 ROH that had a negative effect on milk volume. This region included 18 genes (see Additional file 1). The most promising candidate gene in this region is glutamate-ammonia ligase (GLUL), which is involved in alanine, aspartate, and glutamate metabolism, arginine and proline metabolism and nitrogen metabolism. This gene has previously been reported to control fat to protein ratio . Another region on chromosome 28 that affected milk yield contained eight genes. Bouwman et al  reported the identification of SNPs associated with milk fatty acid composition in this region.
On chromosome 20, two ROH were detected that, when homozygous, had a negative effect on milk yield in both the Holstein and Jersey breeds. Selection signatures close to this region have been identified in the genomes of US and German Holsteins , and the growth hormone receptor (GHR) gene has been proposed as the likely candidate gene , although some publications suggest other possible candidates, such as prolactin receptor (PRLR) ,. Neither PRLR nor GHR aligned with the genes identified in our scan. One hypothesis is that intense selection for GHR or PRLR has made this genomic region globally more homozygous and while homozygosity of specific alleles may be beneficial, overall homozygosity is detrimental to production because of unfavourable alleles that are in linkage disequilibrium with the favourable alleles of GHR/PRLR.
For the Jersey breed, there was a cluster of 14 ROH on the X chromosome, located between 60.4 and 62.1 Mb (Table 6), that were significantly (P <0.001) associated with fertility. A candidate gene in this region (insulin receptor substrate 4; IRS4) has been reported to be associated with mild defects in reproduction and growth in mice . Insulin receptor substrates mediate the actions of insulin and control energy balance and glucose homeostasis .
One of the strengths of this study was that two breeds, Holstein and Jersey, were used. Thus, any effects that were detected in one breed and validated in the other are more likely to be true rather than false discoveries. Unfortunately, there was only one region (on chromosome 24 at around 60 Mb), with an unfavourable association with fertility, that was validated in both breeds. No candidate genes were identified in this region (see Additional file 1), although the ROH could have been in linkage disequilibrium with QTL that were located further away. Our understanding of genomic inbreeding will improve as datasets increase in size and the ability to find causative mutations will also increase as a result of high-density genotyping and sequencing.
One potential use of non-additive effects is to predict the future performance of a cow or of the potential heifer resulting from a planned mating. This prediction would use both the additive and dominance effects and should reduce the production of calves that are homozygous for deleterious mutations. In the past, customised mating programs have focused on genome-wide inbreeding , but, with the availability of genomic data, information on the genomic locations for which deleterious recessive alleles segregate, could be used.
Genomic measures of inbreeding can be used instead of pedigree-based inbreeding to estimate the effects of inbreeding depression and give similar estimates. Milk, fat and protein yields and fertility were unfavourably impacted by increasing genomic and pedigree measures of inbreeding for the Holstein and Jersey breeds of dairy cattle. Although inbreeding was higher for the Jersey than the Holstein breed, the effect of increasing inbreeding on the traits studied was similar for both breeds. In situations where sub-populations are suspected, it may be preferable to consider a measure of inbreeding that does not depend on allele frequencies, such as the proportion of homozygous SNPs. However, in our data, homozygosity did not have a larger effect on milk production and fertility than the diagonal element of the GRM. ROH can provide additional benefits, since longer ROH (>3 Mb) capture recent inbreeding, which had an effect independent of overall homozygosity on milk production. Using a GWAS with ROH of 50 or more SNPs, genomic regions were detected that had an effect of up to 12.5 d and 260 L for calving interval and milk yield, respectively, when completely homozygous. Breeding programs could exploit whole-genome data, as well as these site-specific regions.
JP did the analyses and wrote the paper with BJH, MHM and MEG. MHM prepared the phenotype data and calculated pedigree inbreeding coefficients. All authors read and approved the final manuscript.
The authors acknowledge the Department of Environment and Primary Industries, Victoria, Australia and Dairy Futures Co-operative Research Council, Melbourne, Australia for funding this research and the Australian Dairy Herd Improvement Scheme (ADHIS, Melbourne, Australia) for providing the data.
- Keller LF, Waller DM: Inbreeding effects in wild populations. Trends Ecol Evol. 2002, 17: 230-241. 10.1016/S0169-5347(02)02489-8.View ArticleGoogle Scholar
- Meagher S, Penn DJ, Potts WK: Male-male competition magnifies inbreeding depression in wild house mice. Proc Natl Acad Sci U S A. 2000, 97: 3324-3329. 10.1073/pnas.97.7.3324.PubMed CentralView ArticlePubMedGoogle Scholar
- Jarne P, Perdieu M, Pernot A, Delay B, David P: The influence of self-fertilization and grouping on fitness attributes in the freshwater snail Physa acuta: population and individual inbreeding depression. J Evol Biol. 2000, 13: 645-655. 10.1046/j.1420-9101.2000.00204.x.View ArticleGoogle Scholar
- Sewalem A, Johansson K, Wilhelmson M, Lillpers K: Inbreeding and inbreeding depression on reproduction and production traits of White Leghorn lines selected for egg production traits. Br Poult Sci. 1999, 40: 203-208. 10.1080/00071669987601.View ArticlePubMedGoogle Scholar
- Haile-Mariam M, Bowman PJ, Goddard ME: A practical approach for minimising inbreeding and maximising genetic gain in dairy cattle. Genet Sel Evol. 2007, 39: 369-389. 10.1186/1297-9686-39-4-369.PubMed CentralView ArticlePubMedGoogle Scholar
- Mc Parland S, Kearney JF, Rath M, Berry DP: Inbreeding effects on milk production, calving performance, fertility, and conformation in Irish Holstein-Friesians. J Dairy Sci. 2007, 90: 4411-4419. 10.3168/jds.2007-0227.View ArticlePubMedGoogle Scholar
- Smith LA, Cassell BG, Pearson RE: The effects of inbreeding on the lifetime performance of dairy cattle. J Dairy Sci. 1998, 81: 2729-2737. 10.3168/jds.S0022-0302(98)75830-8.View ArticlePubMedGoogle Scholar
- Charlier C, Coppieters W, Rollin F, Desmecht D, Agerholm JS, Cambisano N, Carta E, Dardano S, Dive M, Fasquelle C, Frennet JC, Hanset R, Hubin X, Jorgensen C, Karim L, Kent M, Harvey K, Pearce BR, Simon P, Tama N, Nie H, Vandeputte S, Lien S, Longeri M, Fredholm M, Harvey RJ, Georges M: Highly effective SNP-based association mapping and management of recessive defects in livestock. Nat Genet. 2008, 40: 449-454. 10.1038/ng.96.View ArticlePubMedGoogle Scholar
- VanRaden PM: Efficient methods to compute genomic predictions. J Dairy Sci. 2008, 91: 4414-4423. 10.3168/jds.2007-0980.View ArticlePubMedGoogle Scholar
- Hayes BJ, Visscher PM, Goddard ME: Increased accuracy of artificial selection by using the realized relationship matrix. Genet Res (Camb). 2009, 91: 47-60. 10.1017/S0016672308009981.View ArticleGoogle Scholar
- Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, Madden PA, Heath AC, Martin NG, Montgomery GW, Goddard ME, Visscher PM: Common SNPs explain a large proportion of the heritability for human height. Nat Genet. 2010, 42: 565-569. 10.1038/ng.608.PubMed CentralView ArticlePubMedGoogle Scholar
- MacLeod IM, Meuwissen TH, Hayes BJ, Goddard ME: A novel predictor of multilocus haplotype homozygosity: comparison with existing predictors. Genet Res (Camb). 2009, 91: 413-426. 10.1017/S0016672309990358.View ArticleGoogle Scholar
- Keller MC, Simonson MA, Ripke S, Neale BM, Gejman PV, Howrigan DP, Lee SH, Lencz T, Levinson DF, Sullivan PF: Runs of homozygosity implicate autozygosity as a schizophrenia risk factor. PLoS Genet. 2012, 8: e1002656-10.1371/journal.pgen.1002656.PubMed CentralView ArticlePubMedGoogle Scholar
- Bjelland DW, Weigel KA, Vukasinovic N, Nkrumah JD: Evaluation of inbreeding depression in Holstein cattle using whole-genome SNP markers and alternative measures of genomic inbreeding. J Dairy Sci. 2013, 96: 4697-4706. 10.3168/jds.2012-6435.View ArticlePubMedGoogle Scholar
- Keller MC, Visscher PM, Goddard ME: Quantification of inbreeding due to distant ancestors and its detection using dense single nucleotide polymorphism data. Genetics. 2011, 189: 237-249. 10.1534/genetics.111.130922.PubMed CentralView ArticlePubMedGoogle Scholar
- Purfield DC, Berry DP, McParland S, Bradley DG: Runs of homozygosity and population history in cattle. BMC Genet. 2012, 13: 70-10.1186/1471-2156-13-70.PubMed CentralView ArticlePubMedGoogle Scholar
- Holt M, Meuwissen T, Vangen O: The effect of fast created inbreeding on litter size and body weights in mice. Genet Sel Evol. 2005, 37: 523-537. 10.1186/1297-9686-37-6-523.PubMed CentralView ArticlePubMedGoogle Scholar
- Hinrichs D, Meuwissen TH, Odegard J, Holt M, Vangen O, Woolliams JA: Analysis of inbreeding depression in the first litter size of mice in a long-term selection experiment with respect to the age of the inbreeding. Heredity (Edinb). 2007, 99: 81-88. 10.1038/sj.hdy.6800968.View ArticleGoogle Scholar
- Kim ES, Cole JB, Huson H, Wiggans GR, Van Tassell CP, Crooker BA, Liu G, Da Y, Sonstegard TS: Effect of artificial selection on runs of homozygosity in US Holstein cattle. PLoS One. 2013, 8: e80813-10.1371/journal.pone.0080813.PubMed CentralView ArticlePubMedGoogle Scholar
- Haile-Mariam M, Bowman PJ, Pryce JE: Genetic analyses of fertility and predictor traits in Holstein herds with low and high mean calving intervals and in Jersey herds. J Dairy Sci. 2013, 96: 655-667. 10.3168/jds.2012-5671.View ArticlePubMedGoogle Scholar
- Cassell BG, Adamec V, Pearson RE: Effect of incomplete pedigrees on estimates of inbreeding and inbreeding depression for days to first service and summit milk yield in Holsteins and Jerseys. J Dairy Sci. 2003, 86: 2967-2976. 10.3168/jds.S0022-0302(03)73894-6.View ArticlePubMedGoogle Scholar
- Sørensen AC, Sørensen MK, Berg P: Inbreeding in Danish dairy cattle breeds. J Dairy Sci. 2005, 88: 1865-1872. 10.3168/jds.S0022-0302(05)72861-7.View ArticlePubMedGoogle Scholar
- Meuwissen T, Luo Z: Computing inbreeding coefficients in large populations. Genet Sel Evol. 1992, 24: 305-313. 10.1186/1297-9686-24-4-305.PubMed CentralView ArticleGoogle Scholar
- Hayes B: Efficient parentage assignment and pedigree reconstruction with dense single nucleotide polymorphism data. J Dairy Sci. 2011, 94: 2114-2117. 10.3168/jds.2010-3896.View ArticlePubMedGoogle Scholar
- Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O’Connell J, Moore SS, Smith TP, Sonstegard TS, Van Tassell CP: Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009, 4: e5350-10.1371/journal.pone.0005350.PubMed CentralView ArticlePubMedGoogle Scholar
- Erbe M, Hayes BJ, Matukumalli LK, Goswami S, Bowman PJ, Reich CM, Mason BA, Goddard ME: Improving accuracy of genomic predictions within and between dairy cattle breeds with imputed high-density single nucleotide polymorphism panels. J Dairy Sci. 2012, 95: 4114-4129. 10.3168/jds.2011-5019.View ArticlePubMedGoogle Scholar
- Gilmour AR, Gogel B, Cullis B, Thompson R: ASReml user guide release 3.0. 2009, VSN International Ltd, Hemel HempsteadGoogle Scholar
- Hu ZL, Park CA, Wu XL, Reecy JM: Animal QTLdb: an improved database tool for livestock animal QTL/association data dissemination in the post-genome era. Nucleic Acids Res. 2013, 41: D871-D879. 10.1093/nar/gks1150.PubMed CentralView ArticlePubMedGoogle Scholar
- Bolormaa S, Pryce JE, Hayes BJ, Goddard ME: Multivariate analysis of a genome-wide association study in dairy cattle. J Dairy Sci. 2010, 93: 3818-3833. 10.3168/jds.2009-2980.View ArticlePubMedGoogle Scholar
- VanRaden PM, Olson KM, Wiggans GR, Cole JB, Tooker ME: Genomic inbreeding and relationships among Holsteins, Jerseys, and Brown Swiss. J Dairy Sci. 2011, 94: 5673-5682. 10.3168/jds.2011-4500.View ArticlePubMedGoogle Scholar
- Pryce JE, Hayes BJ, Goddard ME: Novel strategies to minimize progeny inbreeding while maximizing genetic gain using genomic information. J Dairy Sci. 2012, 95: 377-388. 10.3168/jds.2011-4254.View ArticlePubMedGoogle Scholar
- Gulisija D, Crow JF: Inferring purging from pedigree data. Evolution. 2007, 61: 1043-1051. 10.1111/j.1558-5646.2007.00088.x.View ArticlePubMedGoogle Scholar
- Forni S, Aguilar I, Misztal I: Different genomic relationship matrices for single-step analysis using phenotypic, pedigree and genomic information. Genet Sel Evol. 2011, 43: 1-10.1186/1297-9686-43-1.PubMed CentralView ArticlePubMedGoogle Scholar
- McQuillan R, Eklund N, Pirastu N, Kuningas M, McEvoy BP, Esko T, Corre T, Davies G, Kaakinen M, Lyytikäinen LP, Kristiansson K, Havulinna AS, Gögele M, Vitart V, Tenesa A, Aulchenko Y, Hayward C, Johansson A, Boban M, Ulivi S, Robino A, Boraska V, Igl W, Wild SH, Zgaga L, Amin N, Theodoratou E, Polašek O, Girotto G, Lopez LM: Evidence of inbreeding depression on human height. PLoS Genet. 2012, 8: e1002655-10.1371/journal.pgen.1002655.PubMed CentralView ArticlePubMedGoogle Scholar
- Miglior F, Burnside EB, Kennedy BW: Production traits of Holstein cattle: estimation of nonadditive genetic variance components and inbreeding depression. J Dairy Sci. 1995, 78: 1174-1180. 10.3168/jds.S0022-0302(95)76735-2.View ArticlePubMedGoogle Scholar
- Wiggans GR, VanRaden PM, Zuurbier J: Calculation and use of inbreeding coefficients for genetic evaluation of United States dairy cattle. J Dairy Sci. 1995, 78: 1584-1590. 10.3168/jds.S0022-0302(95)76782-0.View ArticlePubMedGoogle Scholar
- Croquet C, Mayeres P, Gillon A, Vanderick S, Gengler N: Inbreeding depression for global and partial economic indexes, production, type, and functional traits. J Dairy Sci. 2006, 89: 2257-2267. 10.3168/jds.S0022-0302(06)72297-4.View ArticlePubMedGoogle Scholar
- Fuerst C, Solkner J: Additive and nonadditive genetic variances for milk yield, fertility, and lifetime performance traits of dairy cattle. J Dairy Sci. 1994, 77: 1114-1125. 10.3168/jds.S0022-0302(94)77047-8.View ArticlePubMedGoogle Scholar
- Hayes B, Bowman P, Chamberlain A, Goddard M: Invited review: genomic selection in dairy cattle: progress and challenges. J Dairy Sci. 2009, 92: 433-443. 10.3168/jds.2008-1646.View ArticlePubMedGoogle Scholar
- Ferenčaković M, S Lkner J, Curik I: Estimating autozygosity from high-throughput information: effects of SNP density and genotyping errors. Genet Sel Evol. 2013, 45: 42-10.1186/1297-9686-45-42.PubMed CentralView ArticlePubMedGoogle Scholar
- VanRaden PM, Olson KM, Null DJ, Hutchison JL: Harmful recessive effects on fertility detected by absence of homozygous haplotypes. J Dairy Sci. 2011, 94: 6153-6161. 10.3168/jds.2011-4624.View ArticlePubMedGoogle Scholar
- Fritz S, Capitan A, Djari A, Rodriguez SC, Barbat A, Baur A, Grohs C, Weiss B, Boussaha M, Esquerre D, Klopp C, Rocha D, Boichard D: Detection of haplotypes associated with prenatal death in dairy cattle and identification of deleterious mutations in GART, SHBG and SLC37A2. PLoS One. 2013, 8: e65550-10.1371/journal.pone.0065550.PubMed CentralView ArticlePubMedGoogle Scholar
- Tetens J, Seidenspinner T, Buttchereit N, Thaller G: Whole-genome association study for energy balance and fat/protein ratio in German Holstein bull dams. Anim Genet. 2013, 44: 1-8. 10.1111/j.1365-2052.2012.02357.x.View ArticlePubMedGoogle Scholar
- Bouwman AC, Bovenhuis H, Visker MHPW, van Arendonk JAM: Genome-wide association of milk fatty acids in Dutch dairy cattle. BMC Genet. 2011, 12: 43-10.1186/1471-2156-12-43.PubMed CentralView ArticlePubMedGoogle Scholar
- Qanbari S, Pimentel E, Tetens J, Thaller G, Lichtner P, Sharifi AR, Simianer H: A genome-wide scan for signatures of recent selection in Holstein cattle. Anim Genet. 2010, 41: 377-389.PubMedGoogle Scholar
- Viitala S, Szyda J, Blott S, Schulman N, Lidauer M, Mäki-Tanila A, Georges M, Vilkki J: The role of the bovine growth hormone receptor and prolactin receptor genes in milk, fat and protein production in Finnish Ayrshire dairy cattle. Genetics. 2006, 173: 2151-2164. 10.1534/genetics.105.046730.PubMed CentralView ArticlePubMedGoogle Scholar
- Chamberlain A, Hayes B, Savin K, Bolormaa S, McPartlan HC, Bowman PJ, Van Der Jagt C, MacEachern S, Goddard ME: Validation of single nucleotide polymorphisms associated with milk production traits in dairy cattle. J Dairy Sci. 2012, 95: 864-875. 10.3168/jds.2010-3786.View ArticlePubMedGoogle Scholar
- Fantin VR, Wang Q, Lienhard GE, Keller SR: Mice lacking insulin receptor substrate 4 exhibit mild defects in growth, reproduction, and glucose homeostasis. Am J Physiol-Endocrinol Metab. 2000, 278: E127-E133.PubMedGoogle Scholar
- Sadagurski M, Dong XC, Myers MG, White MF: Irs2 and Irs4 synergize in non-LepRb neurons to control energy balance and glucose homeostasis. Mol Metab. 2014, 3: 55-63. 10.1016/j.molmet.2013.10.004.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun C, VanRaden PM, O’Connell JR, Weigel KA, Gianola D: Mating programs including genomic relationships and dominance effects. J Dairy Sci. 2013, 96: 8014-8023. 10.3168/jds.2013-6969.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.