- Open Access
Identification of genomic regions associated with inbreeding depression in Holstein and Jersey dairy cattle
Genetics Selection Evolutionvolume 46, Article number: 71 (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.
In this study, we used data on Holstein and Jersey dairy cows with pedigree data and phenotypes for milk production and fertility traits and for which a subset had SNP genotypes to:
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 .
Only cows with at least two generations of complete pedigree and born after 1994 were retained. The proportion of cows with complete pedigree over the first four generations was approximately 0.75 in both the Holstein and Jersey populations and was derived using the following equation ,:
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.
The following measures of inbreeding were calculated using the genomic data:
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):
For the Holstein breed, 887 561 lactation records for milk, fat and protein yields were available for 338 696 cows (i.e. each cow had on average 2.6 lactation records) and 755 618 calving interval records (a measure of fertility that records the interval between two consecutive calving dates in days) were available for 299 590 cows. In Australia, calving interval is available on more cows than measures of fertility calculated from insemination and pregnancy diagnosis data , which is why calving interval was chosen as the fertility trait for analysis. For the Jersey breed, 179 108 milk, fat and protein yield records (on average 2.8 lactation records per cow) were available for 64 049 animals and 153 347 calving interval records were available for 57 049 animals. Phenotypic means of the traits included in the analysis are in Table 1. The average number of lactation records for cows with genotypes was 3.84 and 3.35 for the Holstein and Jersey data, respectively.
All analyses were conducted separately for the Holstein and Jersey datasets. Following previous research using a very similar dataset by Haile-Mariam et al. , we used the following model to estimate the effect of inbreeding on the phenotype being analysed:
where y ijklmn is the nth record of the phenotype (milk, fat or protein yield or calving interval), HYS i the ith herd-year-season of calving, parity j the jth parity, month k the kth 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 lth random permanent environmental effect to account for multiple records on each individual, and cow m was the random genetic effect for the mth 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.
The genome-wide association study (GWAS) was designed to investigate the effect of overlapping ROH at consecutive positions across the genome using sliding windows, so the number of analyses was equal to the number of SNPs. Since the effect of inbreeding estimated using ROH, was of interest in this analysis, the model also included the SNP present at the start of the ROH to correct for its additive effect. In this analysis, milk yield was used to represent milk production traits, and calving interval to represent fertility traits. The model used for the GWAS was:
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).
For the Holstein breed, the mean inbreeding coefficients, calculated using pedigree, were greater for genotyped cows than for the population that had only pedigree-based inbreeding coefficients. However for the Jersey breed, the opposite was true: inbreeding coefficients of genotyped cows were smaller than those of the pedigree recorded population (Table 2; Figure 1). Genomic measures of inbreeding were all higher for the Jersey than for the Holstein population (Table 2). For example, the average proportion of the genome that was homozygous was 0.66 and 0.72 for the Holstein and Jersey populations, respectively. Jersey animals also had longer ROH than Holstein animals (Figure 2). The rate of inbreeding, calculated as the regression of pedigree inbreeding on birth year (1997 to 2007), was 0.13%/year in Holstein cows with pedigree information and 0.18%/year in the genotyped Holstein population. The equivalent rates for the Jersey animals were both 0.13% (Figure 1).
Correlations of pedigree-based inbreeding coefficients with each measure of genomic inbreeding (i.e., GRM_F, proportion homozygous, ROH_F) were lower than correlations among the genomic measures of inbreeding (Table 3). The correlation of the proportion of the genome that was homozygous with inbreeding coefficient GRM_F was equal to 0.72 for Holstein and 0.66 for Jersey animals. The correlation of ROH_F with the pedigree-based inbreeding coefficient was largest when the number of SNPs in the ROH_F was equal to 10 (Figure 3a) for both breeds. However, for the Holstein population, the correlation between GRM_F and ROH_F was greater for ROH equal to 1 (0.66), i.e., the proportion of SNPs that are homozygous, than for ROH equal to 10 (0.54) (Figure 3b). For the Jersey population, the correlations between ROH_F and GRM_F were more consistent and ranged from 0.66 for ROH equal to 1 to 0.62 for ROH equal to 10 (Figure 3b).
For Holstein cows, a 1% increase in pedigree-based inbreeding or in GRM_F was associated with a reduction in milk yield of 21 and 28 L/lactation, respectively (Table 4), which represents 0.3% to 0.4% of the phenotypic mean, respectively. For Jersey cows, a 1% increase in pedigree-based inbreeding or in GRM_F was associated with a reduction in milk yield of 12 and 27 L/lactation, respectively (Table 4), which represents 0.2% and 0.5% of the phenotypic mean, respectively. The effect of inbreeding on fertility was only significant (P <0.05) for the Holstein breed and when using pedigree-based inbreeding coefficients, where a 1% increase in inbreeding was associated with an extension of calving interval of +0.18 d, which represents 0.04% of the phenotypic mean (Table 4). Compared to the estimates of inbreeding depression based on the population that had only pedigree information, inbreeding depression effects were generally lower in the genotyped population. Furthermore, in the genotyped population, the effect of inbreeding depression was 2 to 3 times greater based on GRM_F compared to pedigree-based inbreeding for both breeds.
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%).
A 1% increase in homozygous SNPs resulted on average in a 1% and 1.5% reduction in milk, fat and protein yields for Holstein and Jersey cows, respectively (Table 5). The association with fertility was significant for the Holstein breed (P <0.05) and was equal to 0.93 d per 1% increase in homozygous SNPs (Table 5). Increasing GRM_F by 1% increases homozygosity by on average 0.33%, so the effect of a 1% increase in homozygosity is expected to be about 3 times the effect of a 1% increase in GRM_F, which is approximately what we observed.
Longer ROH, which indicates a more recent common ancestor, had a stronger unfavourable effect on milk yield than shorter ROH; a ROH of five SNPs led to a much smaller reduction in milk yield than a ROH of 100 SNPs (Figure 4a). However, general homozygosity may explain part of this relationship. To test this, the model was extended to include the overall level of homozygosity as an additional covariate. At the same level of overall homozygosity, shorter ROH (<60 SNPs) had no statistically significant association with milk yield in both breeds. However, longer ROH were associated with a reduction in milk yield that was independent of the proportion of the genome that was homozygous in the Holstein breed (Figure 4b).
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.
Negative associations of ROH with milk yield were found on chromosomes 11, 14, 16, 17, 20, 26, 28 for the Holstein breed and on chromosomes 8, 17, 20 and 24 for the Jersey breed (Table 6). Unfavourable associations of ROH with calving interval were found on chromosomes 2, 5, 8, 9, 15 and 24 for the Holstein breed and on chromosomes 24 and X for the Jersey breed.
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.
Keller LF, Waller DM: Inbreeding effects in wild populations. Trends Ecol Evol. 2002, 17: 230-241. 10.1016/S0169-5347(02)02489-8.
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.
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.
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.
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.
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.
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.
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.
VanRaden PM: Efficient methods to compute genomic predictions. J Dairy Sci. 2008, 91: 4414-4423. 10.3168/jds.2007-0980.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Meuwissen T, Luo Z: Computing inbreeding coefficients in large populations. Genet Sel Evol. 1992, 24: 305-313. 10.1186/1297-9686-24-4-305.
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.
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.
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.
Gilmour AR, Gogel B, Cullis B, Thompson R: ASReml user guide release 3.0. 2009, VSN International Ltd, Hemel Hempstead
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.
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.
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.
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.
Gulisija D, Crow JF: Inferring purging from pedigree data. Evolution. 2007, 61: 1043-1051. 10.1111/j.1558-5646.2007.00088.x.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The authors declare they have no competing interests.