Skip to main content
  • Research article
  • Open access
  • Published:

Runs of homozygosity and population history in cattle

Abstract

Background

Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an individual due to parents transmitting identical haplotypes to their offspring. The extent and frequency of ROHs may inform on the ancestry of an individual and its population. Here we use high density (n = 777,962) bi-allelic SNPs in a range of cattle breed samples to correlate ROH with the pedigree-based inbreeding coefficients and to validate subsequent analyses using 54,001 SNP genotypes. This study provides a first testing of the inference drawn from ROH through comparison with estimates of inbreeding from calculations based on the detailed pedigree data available for several breeds.

Results

All animals genotyped on the HD panel displayed at least one ROH that was between 1–5 Mb in length with certain regions of the genome more likely to be involved in a ROH than others. Strong correlations (r = 0.75, p < 0.0001) existed between the pedigree-based inbreeding coefficient and a statistic based on sum of ROH of length > 0.5 KB and suggests that in the absence of an animal’s pedigree data, the extent of a genome under ROH may be used to infer aspects of recent population history even from relatively few samples.

Conclusions

Our findings suggest that ROH are frequent across all breeds but differing patterns of ROH length and burden illustrate variations in breed origins and recent management.

Background

Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an animal due to parents transmitting identical haplotypes to their offspring. The extent and frequency of these may inform on the ancestry of an individual and its population. Particularly, consanguinity may be indicated from the presence of long ROH; the longer such segments are, the more likely that recent inbreeding occurred within a pedigree [1]. However, unusually long runs of homozygosity may also persist in out bred individuals, perhaps due to unusual mutation, linkage disequilibrium (LD), and recombination rates at certain genomic locations [2].

The distribution of shorter ROH may also inform on the presence of more ancient relatedness which is unaccounted for in an individual’s recorded pedigree due to the limitations in the recording process [3]. In the context of domestic animals these may result from breed or population founder effects or other restrictions.

The domestication of cattle, which occurred ~10,000 years ago, was a complex process, with evidence suggesting that it occurred in a minimum of two domestication events [4]. However, both natural and artificial selection of cattle, as well as regional variations due to drift has resulted in breeds that differ extensively in phenotypes. These processes and the extent of breeding control have differed greatly among populations and ROH may provide useful information on these disparate histories. Particularly, in recent times the practices of intense selection of sires, artificial insemination, and embryo transfer have featured heavily in some breeds, reducing effective population sizes, genetic diversity and affecting levels of homozygosity.

Runs of homozygosity have been extensively studied in human populations and are an established method of distinguishing a population history of consanguinity, and with homozygosity mapping analysis showing a relationship with susceptibility to recessive diseases [13, 5, 6]. Here we use high density (n = 777,962) bi-allelic SNP data in a range of cattle breeds to correlate ROH with the pedigree-based inbreeding coefficient and to validate further analysis using 54,001 SNP genotypes. This allows examination and interpretation of the level of ROH that exist in a wide range of cattle breeds samples.

Methods

Genotypes and quality control

Single nucleotide polymorphisms (SNPs) genotypes consisting of 777,972 bialleleic SNPs from the BovineHD BeadChip (Illumina Inc., San Diego, CA) were generated for 891 artificial insemination sires of multiple breeds. Breeds represented included Angus (n = 39), Belgian Blue (n = 38), Charolais (n = 117), Friesian (n = 98), Hereford (n = 40), Holstein (n = 262), Holstein-Friesian crosses (n = 111), Limousin (n = 128) and Simmental (n = 58). Additionally, the 48,734 SNPs common to both the HD and the Illumina BovineSNP50 Beadchip were retained in a reduced HD panel which was used to evaluate the lesser circa 50 k SNP density for identifying ROH; hereafter referred to as the reduced HD panel . In addition, three published genotype datasets at this density comprising of a total of 1166 animals from 42 different breeds (detailed in Additional file 1.) generated from the Illumina Bovine50 Beadchip, hereafter known as the SNP50 panel, were also available [79].

The impact of the following SNP edits on the number of SNPs remaining for each of the three density panels is summarised in Table 1. Only biallelic SNPs on the 29 autosomes were retained and both animals and SNPs with call rates <90% were discarded. Also, monomorphic SNPs and SNPs that deviated from the Hardy Weinberg equilibrium (p < 0.0001) within breed were discarded.

Table 1 Number of single nucleotide polymorphisms for each data edit implemented on each density panel

Definition of a run of homozygosity

Runs of homozygosity (ROH) were defined in the populations of animals for each of the three SNP density panels using PLINK v1.07, which slides a window of 50 SNPs, in one SNP intervals, across the genome estimating homozygosity. No more than two SNPs with missing genotypes were allowed per window and up to one possible heterozygous genotype was permitted. To ensure low SNP density did not affect ROH length, the minimum required density of SNPs per kb differed between the two genotyping densities, as well as the maximum distance length between two consecutive homozygous SNPs in a run. To minimize the number of ROH that occur by chance in the HD panel, the minimum number of SNPs that constituted a ROH (l) was calculated by a method similar to that proposed by Lencz et al., [6].

l = log e a n s . n i log e ( 1 h e t ) ¯

where ns is the number of SNPs per individual, ni is the number of individuals, α is the percentage of false positive ROH (set to 0.05 in the present study), h e t ¯ is the mean SNP heterozygosity across all SNPs. For the exclusion of very short and common ROH that occur prevalently throughout the genome due to LD, a minimum ROH length of 500 kb was set.

For analysis of the HD panel genotypes, the minimum SNP density was 1 SNP every 50 kb to ensure low SNP density did not falsify ROH length, a minimum run length of 58 SNPs was needed to produce <5% randomly generated ROH and the maximum gap between two consecutive homozygous SNPs in a run was set at a 100 kb. In the analysis of the reduced HD panel and the Bovine SNP50 panel, the minimum SNP density was altered to 1 SNP every 120 kb and no restriction was placed on the minimum number of SNPs in a ROH and the maximum gap length between two consecutive homozygous SNPs in a run was kept at the default value of 1000 kb to account for the lesser genotype density.

Runs of homozygosity were identified for each animal separately on the HD and reduced panel and for comparative purposes only, the maximum gap length in the HD panel for ROH identification was altered to 1000 kb, SNP density to 1 SNP every 120 kb and there was no restriction on the minimum number of SNPs that constituted a ROH, in order to limit bias. To establish that the reduced density also predicts the correct ROH length category for ROH, the extent that reduced panel ROH were correctly assigned to the HD ROH length category was plotted.

Animals with overlapping ROH, and those ROH that were an allelic match, were also identified in the HD panel. The identification of overlapping regions was done by using the sliding window approach as mentioned above, and then for each SNP by calculating the proportion of homozygous windows in the population dataset that overlap that same position. The percentage of animals that had the region with the most overlapping ROH on each chromosome was plotted and the percentage of these overlapping ROH that were an allelic match ≥ 95 % were identified.

The percentage population of each breed with ROH present at different ROH length categories was calculated, as well as the mean overall ROH sum per animal for each breed. The mean sum of ROH within each ROH length category was also calculated by summing all ROH per animal in each ROH length category and averaging this per breed population. The percentage of SNP involvement in ROH was also calculated by counting the amount of times a SNP appeared in a ROH in the population dataset.

Inbreeding coefficient vs. runs of homozygosity

The measure of homozygosity per animal was calculated using the HD panel and reduced HD panel, by a method similar to that as proposed by McQuillan et al., (2008) except in the present study the centromeric region was also included in the calculation,

F ROH = L ROH L AUTO

in which LROH is the sum of ROH per animal above a certain criterion length and LAUTO is the total length of autosome covered by SNPs. LROH was calculated separately as the sum of all ROH >500 kb and the sum of all ROH >10,000 kb where SNP autosomal genome coverage was 2,510,611 kb in the HD panel and 2,500,265 kb in the reduced HD panel. Pedigree based inbreeding coefficients for all animals were calculated using the Meuwissen and Luo [10] algorithm. Depth of pedigree known was measured in complete generation equivalents (CGE) for all animals as described in [11] and correlations between all measures of inbreeding were calculated only on animals (n = 230) with a CGE value ≥6.

Results

ROH in animals with HD panel genotypes

We used a definition of ROH as tracts of homozygous genotypes that were >500 kb in length identified in a genome sliding window of 50 SNPs. No more than two missing genotypes and one possible heterozygous genotype were allowed in a window, and within our HD SNP panel data, all 867 animals tested, displayed at least one ROH that was between 1–5 Mb in length with almost all (i.e., 98%) of the population dataset also having at least one ROH between 5–10 Mb in length. Differences among breeds existed in their frequency in different ROH length categories (Figure 1). The two British breeds, the Angus and the Hereford, had a larger mean portion of their genome, 198.6 Mb and 198.7 Mb, respectively, covered in shorter length (i.e. 1–5 Mb) ROH; coverage ranged from 80.58 to 93.48 Mb in the remaining breeds. For all breeds, most ROH segment coverage was in the shorter length categories with the Holstein, Holstein-Friesian and Friesian samples showing the greatest frequency in the longer ROH categories which are more indicative of recent inbreeding. A smaller effective population size and its associated unfavourable effects due to inbreeding depression in these dairy cattle is a current concern in the industry [12]. The mean sum of the ROH per animal for all ROH length categories > 15-20 Mb was less than 20 Mb and no animal had more than five ROH that were >30 Mb in length. The three most homozygous animals present in our dataset were from the dairy breeds and had on average 700.3 Mb classified as ROH; this is equivalent to almost a quarter of their genome.

Figure 1
figure 1

The mean sum of Run of Homozygosity (ROH) per HD genotyped animal. The sum of ROH was calculated per animal, measured in megabases (Mb) within each ROH length category and averaged per breed. Breeds from left to right are Angus, Belgian Blue, Charolais, Friesian, Hereford, Holstein, Holstein-Friesian, Limousin and Simmental. Angus and Hereford had high amount of mean short ROH (<5 Mb) possibly due to ancestral relatedness whereas the Holsteins had the greater mean long ROH (>20 Mb) due to recent consanguineous matings.

Genome locations of ROH

The number of ROH per chromosome was greatest for chromosome 1 (11,513 runs across all 867 animals with the HD panel) with on average, 12.39% of the chromosome consisting of a ROH. The number of ROH per chromosome tended to decrease with chromosome length (see Additional file 2). The fraction of chromosome residing in ROH was greatest on chromosomes 14 and 16, with 13.71% and 14.16%, respectively; these had the largest extent of overlapping ROH positions among individuals (see Additional file 3). In particular, 756 animals (87.2% of the dataset) had a ROH segment on chromosome 14 involving the same genomic region centred around the 25 Mb position that consisted of 28 SNPs and was 127.3 kb in length (Figure 2; Additional file 4). When interrogated using Haploview [13] the majority of SNPs in this region were in high LD with each other ( Additional file 5). Also, regions on each of chromosomes 7, 8, 14, 16, 18 featured in a ROH in over 55 % of the sampled animals (Figure 2). Chromosome 12 displayed the longest run of contiguous SNPs (197) that were uninvolved in any ROH within the sample population and interrogation in Haploview [13] revealed that this genomic region had corresponding low levels of LD between SNPs ( Additional file 6). The overlapping ROH regions on chromosomes 8 and 23 were found to be the most allelically similar with over 80% of these overlaps having a >95% allelic match. As this was only examined in the HD population it was difficult to discern a geographical pattern affecting this but the Hereford breed displayed the strongest effect on both chromosomes.

Figure 2
figure 2

Incidence of each single nucleotide polymorphism (SNP) in a Run of Homozygosity (ROH) in the HD panel population. Particular genomic positions can be seen, on each of chromosomes 7, 8, 14, 16, 18 featured in ROH which were shared in over 55% of the sampled animals.

Correlation between ROH and inbreeding coefficient

Cattle offer an opportunity to assess ROH as an indicator of inbreeding by comparing pedigree-based estimates of inbreeding using the same individuals. Figure 3 illustrates a scatter plot of either individual animal FROH or FROH10000 (the sums of all ROH per animal that are >0.5 Mb and >10 Mb respectively) against the inbreeding coefficient based on pedigree recorded (McParland et al., 2007). Conservatively, we restricted our sample to the 230 animals that had complete generation equivalent values ≥6 of recorded pedigree ancestry. A clear linear relationship exists between the ROH and pedigree based estimates of inbreeding. The correlations between the pedigree inbreeding coefficient and individual sums of ROH were different from zero with values of r = 0.75 (P < 0.0001) for FROH and r = 0.71 (P < 0.0001) for FROH10000 using the HD panel and r = 0.73 (P < 0.0001) for FROH and r = 0.70 (P < 0.0001) for FROH10000 using the reduced HD panel. The intercept (0.088, s.e = 0.002) of the regression of FROH on the pedigree inbreeding coefficient was greater than zero suggesting that the pedigree-based inbreeding coefficient may underestimate the levels of ancestral genomic relatedness that may exist. The intercept (0.006; s.e = 0.002) of the regression of FROH10000 on the pedigree inbreeding coefficient was also greater than zero, but to a lesser extent when the 10 Mb ROH length limit is applied; consistent with longer ROH resulting from more recent inbreeding, such as that recorded in these pedigrees.

Figure 3
figure 3

Scatter plot of A) F ROH on pedigree inbreeding coefficient B) F ROH10000 on pedigree Inbreeding coefficient. FROH is the sum of all runs of homozygosity (ROH) per animal that are >0.5 Mb, FROH10000 is the sum of all ROH per animal that are >10 Mb and the pedigree inbreeding coefficient is based on pedigree recorded data calculated by the Meuwissen and Lou algorithm (1992).

Validation of Bovine SNP50 density

The most common SNP density platform for cattle is that from 54,000 SNP genotyping arrays (Illumina BovineSNP50 beadchip); an important question is whether this panel can reliably infer ROH that have been detected with the greater accuracy provided by the HD array. To facilitate this, a reduced genotype dataset comprising 48,734 SNPs that were common to both the HD and Bovine50 Beadchip genotype panel on the 867 genotyped sires was used for comparison purposes. The same ROH criteria were applied to each dataset, where no more than two missing genotypes and one possible heterozygote were permitted per window, no restriction was placed on the number of SNPs in a ROH in either density and the maximum gap between two consecutive SNPs was kept at the default value of 1000 kb to account for the lesser genotype density in the reduced HD panel and to limit bias for comparison of ROH. There were in total 157,600 ROH identified using the HD panel, whereas only 19,078 ROH were identified in the same individuals using the reduced HD panel. However, the vast majority of ROH that the reduced panel failed to recognise were between 0.5-1 Mb in length (Figure 4). Also, only 27.7% of the 1–5 Mb ROH length category in the HD panel were identified by the reduced panel but almost all ROH in length categories of 5 Mb and above were identified (Figure 4). In the larger ROH length categories, ROH from the HD panel that were recognised by the reduced panel were correctly assigned >70% of the time to the correct length category where the majority of incorrectly assigned ROH were due to inflation of length in the reduced HD panel. When the reduced panel density ROH estimates were compared to pedigree inbreeding values using the same individuals the correlations obtained were very similar (0.73;P < 0.0001) using FROH and 0.70 (P < 0.0001) using FROH10000 to those observed for the HD genotype panel.

Figure 4
figure 4

The percentage of Runs of Homozygosity (ROH) in the HD panel that are recognised by the reduced density panel. Also the percentage of the reduced panel ROH that are correctly identified into the correct ROH length category is also plotted. Greater than 70 % of the ROH identified using the reduced panel were placed in the correct length category. The reduced density failed to identify ROHs between 0.5-1 Mb in length due to the low SNP density.

ROH in Bovine SNP50 genotypes across the 42 breeds

The success of SNP50 validation exercise allows the testing of a large previously published data collection of wide provenance [79] for ROH distributions using 1 Mb as a ROH minimum. Figure 5 describes the total (Mb) per individual genome estimated to be under ROH in different breed groupings. There were clear differences among breeds in both the levels and variation of ROH frequency. Also, there were discernable patterns with respect to breed origin. The British Isles breed samples (Angus, Guernsey, Hereford and Jersey) clearly displayed the highest individual sum of ROH per animal. Mainland European breed samples showed varying ranges of ROH levels and zebu breeds had intermediate levels. African breeds had a tendency towards low ROH but the West African taurine breeds (Baoule, Lagune, N’Dama and Somba) showed high variability in individual sums of ROH within breed samples (Figure 5).

Figure 5
figure 5

Individual sum of total Run of Homozygosity (ROH) per Bovine SNP50 animal split into origin of breed. The British Isles breeds had higher sums of ROH per animal in comparison to other geographic origins. Breeds originating in Africa have the most disparate range of sums a reflection of the current breeding practices.

Long ROH arise as a result of recent inbreeding, while shorter ROH can indicate more distant ancestral effects such as breed founder effects. In order to discriminate between these effects the average sum of the genome under ROH of size >20 Mb against the average sum <20 Mb for each breed were compared (Figure 6). A correlation of r = 0.68 (P < 0.0001) existed between the two measures but the relative contributions of long and short ROH differed sharply among breeds and breed types. The African Bos taurus breeds showed low levels of ROH for both measures, except for the Somba, Oulmès Zaer and Lagune breeds where elevated levels were clearly highly influenced by long ROH and hence recent consanguinity. The two Channel Islands breeds also showed this trend. In contrast, the zebu breeds showed intermediate values of ROH but these were more strongly influenced by shorter ROH, perhaps indicating a founder effect rather than recent inbreeding.

Figure 6
figure 6

Average sum of Runs of Homozygosity (ROH) per bovine SNP50 animal <20 Mb in length within breed VS average sum of ROH >20 Mb in length. Breeds that originate from the same geographical area tend to cluster together with African Bos taurus breeds showing low levels for both measures, except for the Somba, Oulmès Zaer and Lagune breeds where elevated levels were clearly highly influenced by long ROH whereas the zebu breeds showed intermediate values of ROH but these were more strongly influenced by shorter ROH.

Discussion

Our findings show that ROH are frequent across all breeds and that certain ROH length categories can be used as an indication of consanguinity. They can also inform on breed population history as the effects of population bottlenecks, selection pressure and breeding management on the bovine genome may potentially leave an imprint on ROH length.

The bovine HD SNP assay allows an analysis of ROH at similar density to that employed to generate genomic signatures of endogamy that differ markedly among human populations [13]. Moreover, cattle allow a first comprehensive testing of the inference drawn from ROH through comparison with estimates of inbreeding from calculations based on the detailed pedigree data available for many breeds. The strong correlation between the pedigree inbreeding coefficient and sum of ROH of length > 0.5 kb suggests that, in the absence of an animal’s pedigree data, the extent of a genome under ROH may be used to infer aspects of recent population history even from relatively few samples, as previously suggested by McQuillan et al. (2008). However, 44% of the variance in ROH distribution remains unexplained by pedigree inbreeding and may partly reflect the limitations of ancestry recording in cattle where founder animals are generally, and often inaccurately, assumed to be unrelated (McParland et al., 2007). Additionally, the propensity for multiple megabase scale ancestral haplotypes in certain genome regions to persist even in outbred animals, perhaps due to localised low levels of recombination and high levels of LD may contribute [14]. Lastly, we note that pedigree relatedness gives an expected, not actual, proportion of genomic identity by descent among individuals and it might be anticipated that genotype-based estimates provide greater accuracy on relatedness [15].

Whereas HD SNP panels facilitate more accurate detection of ROH, the vast majority of cattle SNP genotype data, and emerging data in other livestock, is available at ~50,000 SNP density. It is therefore of interest as to whether these sparser genotypic data can reliably inform on ROH and inbreeding history. We found that HD ROH were not accurately identified in such a reduced panel if between 0.5-1 Mb in length. However, the SNP50 density genotypes were sufficient to recognise almost all ROH >5 Mb but they also has the potential to inflate ROH length. Importantly, ROH levels at the lower SNP density correlate equally well with the pedigree estimates of inbreeding. We conclude that this prevalent marker density is appropriate in identifying ROH.

We used three published [79] SNP50 genotype collections to examine patterns in ROH distribution and compare aspects of population history among a range of cattle breeds. The domestication process itself featured a limited sampling from the wild with a more recent bottleneck detectable 50–100 generations ago, presumably corresponding to breed formation [16]. However, this traditional breed formation is largely a European phenomenon and its absence is most apparent in the data from African cattle. These samples, including B. taurus breeds, humped B. Indicus breeds and indicine/taurine hybrids, tended toward low levels of ROH per genome, reflecting traditional management practices in Africa, characterised by less controlled mating [17].

An open village breeding system may also predispose to random consanguineous matings and many African breeds show outlying highly inbred individuals (Figure 5). The length distribution of ROH can help to distinguish different types of parental relatedness. Samples from human populations where cousin marriage is common show an excess of long ROH; whereas for example, Papuan and Melanesian human populations show an excess of shorter ROH, consistent with effects of reduced population and isolation rather than first degree relative unions (Kirin et al. 2010). Figure 6 compares the contributions of long (> 20 Mb) and short ROH to breed homozygosity in order to differentiate the effects of ancient and more recent relatedness among ancestors. Here, the three African taurine breeds (Oulmes zaer, Somba and Lagune) with higher homozygosity clearly show a strong influence of ROH of length greater than 20 Mb and hence of recent inbreeding. We note that Gautier et al. (2009) reported a high FIS value, as well as extensive linkage disequilibrium within the Lagune breed.

Zebu-taurus hybridisation is also a dynamic and contemporary process within Africa [1820]. This acts to increase genetic diversity and contributes to the interruption of stretches of homozygous genotypes within individuals. The effects of this process are evident in three hybrid breeds which show the lowest extent of ROH in African breeds. These include the Kuri; where previous work has shown a near 50:50 genetic admixture between surrounding zebu and themselves [21], and the Sheko breed, where the original taurine African Y chromosome is in danger of disappearing from this breed due to the use of zebu bulls [20].

Within European breeds, British breeds tended toward higher quantities of ROH, reflecting results of previous microsatellite research, where British Isles breeds had lower levels of observed heterozygosities and gene diversities in comparison to other Mediterranean and Northern European breeds analysed [22]. The Channel Island breeds showed strong influence from long ROH reflecting their unusually closed population histories due to strict importation restrictions on both the Jersey and Guernsey Islands implemented during the 1800 s [23].

The zebu breeds represented in this study have contrasting histories. The mainland African zebu breeds (Bororo and Fulani) which are products of ancient introductions from South Asia and are all hybrids to some extent, had much lower quantities of average ROH in comparison to the American zebu breeds analysed (Figure 6). Within the American and Madagascan zebu populations a stronger homozygosity signal, with a weighting toward smaller length ROH, suggests that these breeds were initially established by small founding populations but were not particularly affected by recent inbreeding (Figure 6). The initial introduction of the now prolific zebu animals in the Americas featured very limited numbers during the 19th and 20th century and Madagascan zebu were founded by ancient importations from Asia and East Africa which were probably limited in scope due to the isolation of the island [20, 24].

The ascertainment bias [8] towards European Bos taurus breeds that is associated with the Bovine SNP50 genotyping chip, does not seem to invalidate the trends in ROH levels observed here, as ROH levels were in fact higher in those breeds with a higher number of polymorphic SNPs, as validated by Illumina ( Additional file 7). Also, the existence of long ROH (>20 Mb) for example in many of the less polymorphic African village breeds (Oulmes zaer, Somba and Lagune) are unlikely to be artefactual due to the vanishingly small probability of long contiguous homozygous SNPs occurring by chance. However, some bias may exist in the Bos Indicus ROH levels, as an over estimation of ROH amounts is possible due to low amount of polymorphic markers found in these breeds due to the design of the genotyping chip [8], as a result some caution must be taken when inferring ROH levels within these breeds. The bovine HD genotyping chip was designed from a more comprehensive range of breeds comprising several temperate and tropically adapted Bos taurus, Bos indicus and hybrid breeds and thus does not exhibit the same level of ascertainment bias [25].

The Hapmap population data also allow comparison with an alternative inference of past population size. Linkage disequilibrium may be used to infer past population size where higher r2 indicates lower effective population size with LD at longer genetic distances corresponding to younger time depths [26, 27]. Interestingly, the Hapmap breed samples analysed here show a strikingly similar ranking in LD at distances of >200 kb to that which they show in average ROH [28].

Analyses of human ROH have previously established a correlation between extensive LD, locally low rates of recombination and high incidence of homozygous runs [2]. Intensive selection intensity in cattle has possibly acted to maintain long lengths of homozygous tracts. Previous work carried out in over 500 animals from 8 breeds noted that high levels of LD, particularly in the Holstein breed, existed on chromosomes 14 and 16, the two chromosomes with highest proportions of ROH in our study [14]. Conversely, chromosome 12 was found to have higher than average recombination rates and lower levels of LD (r2 <0.2) than the majority of chromosomes [14] and, interestingly, showed the highest proportion of SNPs uninvolved in a ROH within our sample population. The existence of recombination hotspots throughout the genome also can impact ROH, with multiple genomic regions that remained uninvolved in any ROH such as those on chromosomes 12 and 23 found to be well documented human and cattle recombination hotspots [2931].

The existence of QTL in ROH have been well documented in human studies [5, 6, 32]. Here, several of the highly involved genomic regions located on chromosomes 7, 14, 16 and 18 (Figure 2) all potentially contain genes of importance in cattle with associations ranging from immunity through to carcass and dystocia related traits [3335] when explored using three QTL databases available online (http://genomes.sapac.edu.au/bovineqtl/index.html, ahttp://www.animalgenome.org/QTLdb/cattle.html, http://www.ncbi.nlm.nih.gov/). In particular, Chromosomes 9 and 5, which had the highest amount of long ROHs (>20 Mb), are well documented to contain QTL pertaining to milk fat yield and weight related traits respectively [3639].

Conclusion

ROH analyses quantifies a feature of genomic variation that may be used in inference of population history and to associate with important production and disease traits and perhaps signatures of selection. We show that ROH analysis in cattle provides a sufficient predictor of the pedigree inbreeding coefficient and the prevalent SNP50 genotyping array may be a sufficient tool to predict ROH in the genome. Patterns of ROH may be decomposed to highlight the effects of recent and more ancient ancestral relatedness and match known aspects of breed history in a sample set of wide provenance.

References

  1. Kirin M, McQuillan R, Franklin CS, Campbell H, McKeigue PM, Wilson JF: Genomic runs of homozygosity record population history and consanguinity. PLoS One. 2010, 5 (11): e13996-10.1371/journal.pone.0013996.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Gibson J, Morton NE, Collins A: Extended tracts of homozygosity in outbred human populations. Hum Mol Genet. 2006, 15 (5): 789-795. 10.1093/hmg/ddi493.

    Article  CAS  PubMed  Google Scholar 

  3. McQuillan R, Leutenegger AL, Abdel-Rahman R, Franklin CS, Pericic M, Barac-Lauc L, Smolej-Narancic N, Janicijevic B, Polasek O, Tenesa A, et al: Runs of homozygosity in European populations. Am J Hum Genet. 2008, 83 (3): 359-372. 10.1016/j.ajhg.2008.08.007.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Loftus RT, MacHugh DE, Bradley DG, Sharp PM, Cunningham P: Evidence for two independent domestications of cattle. Proc Natl Acad Sci U S A. 1994, 91 (7): 2757-2761. 10.1073/pnas.91.7.2757.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Hildebrandt F, Heeringa SF, Rüschendorf F, Attanasio M, Nürnberg G, Becker C, Seelow D, Huebner N, Chernin G, Vlangos CN, et al: A systematic approach to mapping recessive disease genes in individuals from outbred populations. PLoS Genet. 2009, 5 (1): e1000353-10.1371/journal.pgen.1000353.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Lencz T, Lambert C, DeRosse P, Burdick KE, Morgan TV, Kane JM, Kucherlapati R, Malhotra AK: Runs of homozygosity reveal highly penetrant recessive loci in schizophrenia. Proc Natl Acad Sci U S A. 2007, 104 (50): 19942-19947. 10.1073/pnas.0710021104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Gautier M, Laloë D, Moazami-Goudarzi K: Insights into the genetic history of French cattle from dense SNP data on 47 worldwide breeds. PLoS One. 2010, 5 (9): e13038-10.1371/journal.pone.0013038.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, Heaton MP, O'Connell J, Moore SS, Smith TP, Sonstegard TS, et al: Development and characterization of a high density SNP genotyping assay for cattle. PLoS One. 2009, 4 (4): e5350-10.1371/journal.pone.0005350.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Gautier M, Flori L, Riebler A, Jaffrézic F, Laloé D, Gut I, Moazami-Goudarzi K, Foulley JL: A whole genome Bayesian scan for adaptive genetic divergence in West African cattle. BMC Genomics. 2009, 10: 550-10.1186/1471-2164-10-550.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Meuwissen T, Luo Z: Computing inbreeding coefficients in large populations. Genet Sel Evol. 1992, 24 (4): 305-313. 10.1186/1297-9686-24-4-305.

    Article  PubMed Central  Google Scholar 

  11. Mc Parland S, Kearney JF, Rath M, Berry DP: Inbreeding trends and pedigree analysis of Irish dairy and beef cattle populations. J Anim Sci. 2007, 85 (2): 322-331. 10.2527/jas.2006-367.

    Article  CAS  PubMed  Google Scholar 

  12. 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 (9): 4411-4419. 10.3168/jds.2007-0227.

    Article  CAS  PubMed  Google Scholar 

  13. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.

    Article  CAS  PubMed  Google Scholar 

  14. McKay SD, Schnabel RD, Murdoch BM, Matukumalli LK, Aerts J, Coppieters W, Crews D, Dias Neto E, Gill CA, Gao C, et al: Whole genome linkage disequilibrium maps in cattle. BMC Genet. 2007, 8: 74-

    Article  PubMed Central  PubMed  Google Scholar 

  15. Visscher PM, Medland SE, Ferreira MA, Morley KI, Zhu G, Cornes BK, Montgomery GW, Martin NG: Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings. PLoS Genet. 2006, 2 (3): e41-10.1371/journal.pgen.0020041.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Gautier M, Faraut T, Moazami-Goudarzi K, Navratil V, Foglio M, Grohs C, Boland A, Garnier JG, Boichard D, Lathrop GM, et al: Genetic and haplotypic structure in 14 European and African cattle breeds. Genetics. 2007, 177 (2): 1059-1070. 10.1534/genetics.107.075804.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Ibeagha-Awemu EM, Erhardt G: Genetic structure and differentiation of 12 African Bos indicus and Bos taurus cattle breeds, inferred from protein and microsatellite polymorphisms. J Anim Breed Genet. 2005, 122 (1): 12-20. 10.1111/j.1439-0388.2004.00478.x.

    Article  CAS  PubMed  Google Scholar 

  18. Kolbe JJ, Larson A, Losos JB, de Queiroz K: Admixture determines genetic diversity and population differentiation in the biological invasion of a lizard species. Biol Lett. 2008, 4 (4): 434-437. 10.1098/rsbl.2008.0205.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Lavergne S, Molofsky J: Increased genetic variation and evolutionary potential drive the success of an invasive grass. Proc Natl Acad Sci U S A. 2007, 104 (10): 3883-3888. 10.1073/pnas.0607324104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Hanotte O, Tawah CL, Bradley DG, Okomo M, Verjee Y, Ochieng J, Rege JE: Geographic distribution and frequency of a taurine Bos taurus and an indicine Bos indicus Y specific allele amongst sub-saharan African cattle breeds. Mol Ecol. 2000, 9 (4): 387-396. 10.1046/j.1365-294x.2000.00858.x.

    Article  CAS  PubMed  Google Scholar 

  21. Meghen C, MacHugh D, Sauveroche B, Kana G, Bradley DG: Characterization of the Kuri cattle of Lake Chad using molecular genetic techniques. The origin and development of African livestock. Edited by: MacDonald RMBKC. 2000, University College Press, London

    Google Scholar 

  22. Cymbron T, Freeman AR, Isabel Malheiro M, Vigne JD, Bradley DG: Microsatellite diversity suggests different histories for Mediterranean and Northern European cattle populations. Proc Biol Sci. 2005, 272 (1574): 1837-1843. 10.1098/rspb.2005.3138.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Wiener P, Burton D, Williams JL: Breed relationships and definition in British cattle: a genetic analysis. Heredity. 2004, 93 (6): 597-602. 10.1038/sj.hdy.6800566.

    Article  CAS  PubMed  Google Scholar 

  24. MacHugh DE, Shriver MD, Loftus RT, Cunningham P, Bradley DG: Microsatellite DNA variation and the evolution, domestication and phylogeography of taurine and zebu cattle (Bos taurus and Bos indicus). Genetics. 1997, 146 (3): 1071-1086.

    PubMed Central  CAS  PubMed  Google Scholar 

  25. Bovine HD Genotyping BeadChip Datasheet: http://www.illumina.com/documents/products/datasheets/datasheet_bovineHD.pdf,

  26. Hayes BJ, Visscher PM, McPartlan HC, Goddard ME: Novel multilocus measure of linkage disequilibrium to estimate past effective population size. Genome Res. 2003, 13 (4): 635-643. 10.1101/gr.387103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Tenesa A, Navarro P, Hayes BJ, Duffy DL, Clarke GM, Goddard ME, Visscher PM: Recent human effective population size estimated from linkage disequilibrium. Genome Res. 2007, 17 (4): 520-526. 10.1101/gr.6023607.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Gibbs RA, Taylor JF, Van Tassell CP, Barendse W, Eversole KA, Gill CA, Green RD, Hamernik DL, Kappes SM, Lien S, et al: Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science. 2009, 324 (5926): 528-532.

    Article  CAS  PubMed  Google Scholar 

  29. Wang K, Sun D, Xu R, Zhang Y: Identification of 19 new BoLA-DQB alleles. Anim Genet. 2005, 36 (2): 166-167. 10.1111/j.1365-2052.2005.01238.x.

    Article  CAS  PubMed  Google Scholar 

  30. Band M, Larson JH, Womack JE, Lewin HA: A radiation hybrid map of BTA23: identification of a chromosomal rearrangement leading to separation of the cattle MHC class II subregions. Genomics. 1998, 53 (3): 269-275. 10.1006/geno.1998.5507.

    Article  CAS  PubMed  Google Scholar 

  31. Abla N, Chinn LW, Nakamura T, Liu L, Huang CC, Johns SJ, Kawamoto M, Stryke D, Taylor TR, Ferrin TE, et al: The human multidrug resistance protein 4 (MRP4, ABCC4): functional analysis of a highly polymorphic gene. J Pharmacol Exp Ther. 2008, 325 (3): 859-868. 10.1124/jpet.108.136523.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Nalls MA, Guerreiro RJ, Simon-Sanchez J, Bras JT, Traynor BJ, Gibbs JR, Launer L, Hardy J, Singleton AB: Extended tracts of homozygosity identify novel candidate genes associated with late-onset Alzheimer's disease. Neurogenetics. 2009, 10 (3): 183-190. 10.1007/s10048-009-0182-4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Yun S, Lee SH, Yoon SR, Kim MS, Piao ZH, Myung PK, Kim TD, Jung H, Choi I: TOX regulates the differentiation of human natural killer cells from hematopoietic stem cells in vitro. Immunol Lett. 2011, 136 (1): 29-36. 10.1016/j.imlet.2010.11.008.

    Article  CAS  PubMed  Google Scholar 

  34. Casas E, Shackelford SD, Keele JW, Koohmaraie M, Smith TP, Stone RT: Detection of quantitative trait loci for growth and carcass composition in cattle. J Anim Sci. 2003, 81 (12): 2976-2983.

    CAS  PubMed  Google Scholar 

  35. Kühn C, Bennewitz J, Reinsch N, Xu N, Thomsen H, Looft C, Brockmann GA, Schwerin M, Weimann C, Hiendleder S, et al: Quantitative trait loci mapping of functional traits in the German Holstein cattle population. J Dairy Sci. 2003, 86 (1): 360-368. 10.3168/jds.S0022-0302(03)73614-5.

    Article  PubMed  Google Scholar 

  36. McClure MC, Morsci NS, Schnabel RD, Kim JW, Yao P, Rolf MM, McKay SD, Gregg SJ, Chapple RH, Northcutt SL, et al: A genome scan for quantitative trait loci influencing carcass, post-natal growth and reproductive traits in commercial Angus cattle. Anim Genet. 2010, 41 (6): 597-607. 10.1111/j.1365-2052.2010.02063.x.

    Article  CAS  PubMed  Google Scholar 

  37. Li C, Basarab J, Snelling WM, Benkel B, Murdoch B, Moore SS: The identification of common haplotypes on bovine chromosome 5 within commercial lines of Bos taurus and their associations with growth traits. J Anim Sci. 2002, 80 (5): 1187-1194.

    CAS  PubMed  Google Scholar 

  38. Lund MS, Guldbrandtsen B, Buitenhuis AJ, Thomsen B, Bendixen C: Detection of quantitative trait loci in Danish Holstein cattle affecting clinical mastitis, somatic cell score, udder conformation traits, and assessment of associated effects on milk yield. J Dairy Sci. 2008, 91 (10): 4028-4036. 10.3168/jds.2007-0290.

    Article  CAS  PubMed  Google Scholar 

  39. Wiener P, Maclean I, Williams JL, Woolliams JA: Testing for the presence of previously identified QTL for milk production traits in new populations. Anim Genet. 2000, 31 (6): 385-395. 10.1046/j.1365-2052.2000.00689.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by Science Foundation Ireland principal investigator award grant number 09/IN.1/B2642. We are grateful to those researchers who have made SNP50 genotypes used in this work publicly available. We are also grateful to the Irish Cattle Breeding Federation for access to the pedigree data.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Daniel G Bradley.

Additional information

Competing interests

The author(s) declares that there are no competing interests.

Authors’ contributions

DCP participated in the design of the study, carried out the ROH identification, the data analysis and drafted the manuscript. DPB participated in the design of the study, the data analysis and helped to draft the manuscript. SMP calculated all pedigree based inbreeding coefficients. DGB conceived of the study, and participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Details of animals genotyped with the Bovine SNP50 density panel.(PDF 31 KB)

12863_2012_1037_MOESM2_ESM.pdf

Additional file 2: The number of runs of homozygosity (ROH), in the HD panel population, per chromosome (bars) and the mean per animal percentage coverage of the chromosome covered by ROH graph.(PDF 82 KB)

12863_2012_1037_MOESM3_ESM.pdf

Additional file 3: The maximum amount of overlapping runs of homozygosity (ROH) per chromosome and the percentage of overlapping ROHs that are a >95% allelic match graph.(PDF 82 KB)

12863_2012_1037_MOESM4_ESM.pdf

Additional file 4: Proportion of appearances of each SNP on chromosome 14 in a ROH across dairy, beef and British animals graph.(PDF 153 KB)

12863_2012_1037_MOESM5_ESM.pdf

Additional file 5: Linkage disequilibrium plot using Haploview of the consensus overlap segment located on chromosome in 86% of the ROH.(PDF 241 KB)

12863_2012_1037_MOESM6_ESM.pdf

Additional file 6: Linkage disequilibrium plot using Haploview of SNPs on chromosome 12 that did not appear in a ROH.(PDF 241 KB)

12863_2012_1037_MOESM7_ESM.pdf

Additional file 7: Average total sum of Runs of Homozygosity (ROH) per animal within breed vs the number of validated polymorphic loci in the Bovine SNP50 genotyping beadchip graph.(PDF 82 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Purfield, D.C., Berry, D.P., McParland, S. et al. Runs of homozygosity and population history in cattle. BMC Genet 13, 70 (2012). https://doi.org/10.1186/1471-2156-13-70

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2156-13-70

Keywords