Challenges in inbreeding estimation of large populations based on Polish Holstein-Friesian cattle pedigree

The aim of this study was to evaluate observed and future inbreeding level in Polish Holstein-Friesian cattle population. In total, over 9.8 mln animals were used in the analysis coming from the pedigree of Polish Federation of Cattle Breeders and Dairy Farmers. Inbreeding level, as an average per birth year, was estimated with the method accounting for missing parent information with the assumption of year 1950 as the base year of the population. If an animal had no ancestral records, an average inbreeding level from its birth year was assigned. Twice the average inbreeding level served as relatedness of the animal to the population, which enabled estimation of inbreeding in its offspring. The future inbreeding of potential offspring was estimated as an average of animals (bulls and cows) available for mating in a certain year. It was observed that 30–50% of animals born between 1985 and 2015 had no relevant ancestral information, which is caused by a high number of new animals and/or entire farms entering the national milk recordings. For the year 2015, the observed inbreeding level was 3.30%, which was more than twice the inbreeding with the classical approach (without missing parent information) and higher by 0.4% than the future inbreeding. The average increase of inbreeding in years 2010–2015 was 0.10%, which is similar to other countries monitored by World Holstein-Friesian Federation. However, the values might be underestimated due to low pedigree completeness. The estimates of future inbreeding suggested that observed inbreeding could be even lower and also increase slower, which indicates a constant need to monitor rate of increase in inbreeding over time. The most important aspect of presented results is the necessity to advise individual farmers to keep precise recordings of the matings on their farm in order to improve the pedigree completeness of Polish Holstein-Friesian and to use suitable mating programs to avoid too rapid growth of inbreeding.     

Molecular and pedigree measures of relatedness provide similar estimates of inbreeding depression in a bottlenecked population

Individual‐based estimates of the degree of inbreeding or parental relatedness from pedigrees provide a critical starting point for studies of inbreeding depression, but in practice wild pedigrees are difficult to obtain. Because inbreeding increases the proportion of genomewide loci that are identical by descent, inbreeding variation within populations has the potential to generate observable correlations between heterozygosity measured using molecular markers and a variety of fitness related traits. Termed heterozygosity‐fitness correlations (HFCs), these correlations have been observed in a wide variety of taxa. The difficulty of obtaining wild pedigree data, however, means that empirical investigations of how pedigree inbreeding influences HFCs are rare. Here, we assess evidence for inbreeding depression in three life‐history traits (hatching and fledging success and juvenile survival) in an isolated population of Stewart Island robins using both pedigree‐ and molecular‐derived measures of relatedness. We found results from the two measures were highly correlated and supported evidence for significant but weak inbreeding depression. However, standardized effect sizes for inbreeding depression based on the pedigree‐based kin coefficients (k) were greater and had smaller standard errors than those based on molecular genetic measures of relatedness (RI), particularly for hatching and fledging success. Nevertheless, the results presented here support the use of molecular‐based measures of relatedness in bottlenecked populations when information regarding inbreeding depression is desired but pedigree data on relatedness are unavailable.     

Detection of selection signatures in Limousin cattle using whole-genome resequencing

Limousin, a renowned beef breed originating from central France, has been selectively bred over the last 100 years to improve economically important traits. We used whole-genome sequencing data from 10 unrelated Limousin bull calves to detect polymorphisms and identify regions under selection. A total of 13 943 766 variants were identified. Moreover, 311 852 bi-allelic SNPs and 92 229 indels located on autosomes were fixed for the alternative allele in all sequenced animals, including the previously reported missense deleterious F94L mutation in MSTN. We performed a whole-genome screen to discover genomic regions with excess homozygosity, using the pooled heterozygosity score and identified 171 different candidate selective sweeps. In total, 68 candidate genes were found in only 57 of these regions, indicating that a large fraction of the genome under selection might lie in non-coding regions and suggesting that a majority of adaptive mutations might be regulatory in nature. Many QTL were found within candidate selective sweep regions, including QTL associated with shear force or carcass weight. Among the putative selective sweeps, we located genes (MSTN, NCKAP5, RUNX2) that potentially contribute to important phenotypes in Limousin. Several candidate regions and genes under selection were also found in previous genome-wide selection scans performed in Limousin. In addition, we were able to pinpoint candidate causative regulatory polymorphisms in GRIK3 and RUNX2 that might have been under selection. Our results will contribute to improved understanding of the mechanisms and targets of artificial selection and will facilitate the interpretation of GWASs performed in Limousin.     

Patchwork: allele-specific copy number analysis of whole-genome sequenced tumor tissue

Whole-genome sequencing of tumor tissue has the potential to provide comprehensive characterization of genomic alterations in tumor samples. We present Patchwork, a new bioinformatic tool for allele-specific copy number analysis using whole-genome sequencing data. Patchwork can be used to determine the copy number of homologous sequences throughout the genome, even in aneuploid samples with moderate sequence coverage and tumor cell content. No prior knowledge of average ploidy or tumor cell content is required. Patchwork is freely available as an R package, installable via R-Forge (http://patchwork.r-forge.r-project.org/).     

Microsatellite measures of inbreeding: a meta-analysis

Abstract Meta-analyses of published and unpublished correlations between phenotypic variation and two measures of genetic variation at microsatellite loci, multilocus heterozygosity (MLH) and mean d², revealed that the strength of these associations are generally weak (mean r < 0.10). Effects on life-history trait variation were significantly greater than zero for both measures over all reported effect sizes ( r = 0. 0856 and 0.0479 for MLH and mean d ², respectively), whereas effects on morphometric traits were not ( r = 0.0052 and r = 0.0038), which is consistent with the prediction that life-history traits exhibit greater inbreeding depression than morphometric traits. Effect sizes reported using mean d ² were smaller and more variable than those reported using MLH, suggesting that MLH may be a better metric for capturing inbreeding depression most of the time. However, analyses of paired effect sizes reported using both measures from the same data did not differ significantly. Several lines of evidence suggest that published effects sizes are upwardly biased. First, effect sizes from published studies were significantly higher than those reported in unpublished studies. Second, fail-safe numbers for reported effect sizes were generally quite low, with the exception of correlations between MLH and life-history traits. Finally, the slope of the regression of effect size on sample size was negative for most sets of traits. Taken together, these results suggest that studies designed to detect inbreeding depression on a life-history trait using microsatellites will need to sample in excess of 600 individuals to detect an average effect size ( r = 0.10) with reasonable statistical power (0.80). Very few published studies have used samples sizes approaching this value.     

Genomic prediction using imputed whole-genome sequence data in Holstein Friesian cattle

Background

In contrast to currently used single nucleotide polymorphism (SNP) panels, the use of whole-genome sequence data is expected to enable the direct estimation of the effects of causal mutations on a given trait. This could lead to higher reliabilities of genomic predictions compared to those based on SNP genotypes. Also, at each generation of selection, recombination events between a SNP and a mutation can cause decay in reliability of genomic predictions based on markers rather than on the causal variants. Our objective was to investigate the use of imputed whole-genome sequence genotypes versus high-density SNP genotypes on (the persistency of) the reliability of genomic predictions using real cattle data.

Methods

Highly accurate phenotypes based on daughter performance and Illumina BovineHD Beadchip genotypes were available for 5503 Holstein Friesian bulls. The BovineHD genotypes (631,428 SNPs) of each bull were used to impute whole-genome sequence genotypes (12,590,056 SNPs) using the Beagle software. Imputation was done using a multi-breed reference panel of 429 sequenced individuals. Genomic estimated breeding values for three traits were predicted using a Bayesian stochastic search variable selection (BSSVS) model and a genome-enabled best linear unbiased prediction model (GBLUP). Reliabilities of predictions were based on 2087 validation bulls, while the other 3416 bulls were used for training.

Results

Prediction reliabilities ranged from 0.37 to 0.52. BSSVS performed better than GBLUP in all cases. Reliabilities of genomic predictions were slightly lower with imputed sequence data than with BovineHD chip data. Also, the reliabilities tended to be lower for both sequence data and BovineHD chip data when relationships between training animals were low. No increase in persistency of prediction reliability using imputed sequence data was observed.

Conclusions

Compared to BovineHD genotype data, using imputed sequence data for genomic prediction produced no advantage. To investigate the putative advantage of genomic prediction using (imputed) sequence data, a training set with a larger number of individuals that are distantly related to each other and genomic prediction models that incorporate biological information on the SNPs or that apply stricter SNP pre-selection should be considered.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0149-x) contains supplementary material, which is available to authorized users.     

Estimation of individual level of inbreeding using relatedness measures in haplodiploids

Although male haploidy in haplodiploid species aids purging of deleterious alleles, haplodiploid animals may nevertheless suffer significant negative effects of inbreeding. The effects may even be stronger in social Hymenoptera because the negative fitness consequences may be expressed at two levels: the individual level (inbred queens) and colony level (inbred workers). Surprisingly, in natural populations the impact of inbreeding on fitness has been studied in very few insects, and even fewer haplodiploid ones. Hence there is currently little understanding of the potential effects of inbreeding. One reason may be the difficulties in estimating inbreeding especially at the individual level, apart from the additional problems posed by haplodiploidy. In order to study the impact of inbreeding, its individual level must be estimated as precisely as possible. When the population pedigree is unknown, relatedness-based estimates of the individual inbreeding coefficient can be used to estimate inbreeding. Here we examine the relationship between inbreeding coefficients and relatedness in diploid and haplodiploid organisms, and provide guidelines for estimating inbreeding both at the individual and the colony level. Received 7 March 2005; revised 18 April 2005, accepted 20 April 2005. An erratum to this article is available at .     

Population structure and inbreeding from pedigree analysis of purebred dogs

Dogs are of increasing interest as models for human diseases, and many canine population-association studies are beginning to emerge. The choice of breeds for such studies should be informed by a knowledge of factors such as inbreeding, genetic diversity, and population structure, which are likely to depend on breed-specific selective breeding patterns. To address the lack of such studies we have exploited one of the world's most extensive resources for canine population-genetics studies: the United Kingdom (UK) Kennel Club registration database. We chose 10 representative breeds and analyzed their pedigrees since electronic records were established around 1970, corresponding to about eight generations before present. We find extremely inbred dogs in each breed except the greyhound and estimate an inbreeding effective population size between 40 and 80 for all but 2 breeds. For all but 3 breeds, >90% of unique genetic variants are lost over six generations, indicating a dramatic effect of breeding patterns on genetic diversity. We introduce a novel index Psi for measuring population structure directly from the pedigree and use it to identify subpopulations in several breeds. As well as informing the design of canine population genetics studies, our results have implications for breeding practices to enhance canine welfare.     

A comparison of approaches to estimate the inbreeding coefficient and pairwise relatedness using genomic and pedigree data in a sheep population

Genome-wide SNP data provide a powerful tool to estimate pairwise relatedness among individuals and individual inbreeding coefficient. The aim of this study was to compare methods for estimating the two parameters in a Finnsheep population based on genome-wide SNPs and genealogies, separately. This study included ninety-nine Finnsheep in Finland that differed in coat colours (white, black, brown, grey, and black/white spotted) and were from a large pedigree comprising 319 119 animals. All the individuals were genotyped with the Illumina Ovine SNP50K BeadChip by the International Sheep Genomics Consortium. We identified three genetic subpopulations that corresponded approximately with the coat colours (grey, white, and black and brown) of the sheep. We detected a significant subdivision among the colour types (F _ST = 5.4%, P<0.05). We applied robust algorithms for the genomic estimation of individual inbreeding (F _SNP) and pairwise relatedness (Φ _SNP) as implemented in the programs KING and PLINK, respectively. Estimates of the two parameters from pedigrees (F _PED and Φ _PED) were computed using the RelaX2 program. Values of the two parameters estimated from genomic and genealogical data were mostly consistent, in particular for the highly inbred animals (e.g. inbreeding coefficient F>0.0625) and pairs of closely related animals (e.g. the full- or half-sibs). Nevertheless, we also detected differences in the two parameters between the approaches, particularly with respect to the grey Finnsheep. This could be due to the smaller sample size and relative incompleteness of the pedigree for them.We conclude that the genome-wide genomic data will provide useful information on a per sample or pairwise-samples basis in cases of complex genealogies or in the absence of genealogical data.     

The evaluation of the usefulness of pedigree verification-dedicated SNPs for breed assignment in three polish cattle populations

The breed assignment in cattle is one of the issues of molecular genetics which needs further testing and development. Although several statistical approaches have been developed to enable such application, the obtained results strongly depend on specific populations differentiation and power of markers discrimination or their informativeness. Currently, all breeding animals are being tested for parentage with the use of panel of 12 microsatellite markers, which in near future probably will be replaced by about 100 single nucleotide polymorphisms (SNPs). Despite the fact that SNPs are mainly bi-allelic, the multilocus genotypes can reach the level of polymorphism of a panel of microsatellite markers. In this study we attempted to determine the breed of origin of 741 cattle by using 120 SNPs dedicated for parentage testing and included in the BovineSNP50 BeadChip genotyping assay (Illumina). The applied Bayesian and frequency-based methods allowed such differentiation, however, the reliability of the results was not completely satisfying, suggesting that the studied markers are not the best tool for breed assignment.     

Shade material evaluation using a cattle response model and meteorological instrumentation

Shade structures are often considered as one method of reducing stress in feedlot cattle. Selection of a suitable shade material can be difficult without data that quantify material effectiveness for stress reduction. A summer study was conducted during 2007 using instrumented shade structures in conjunction with meteorological measurements to estimate relative effectiveness of various shade materials. Shade structures were 3.6 m by 6.0 m by 3.0 m high at the peak and 2.0 m high at the sides. Polyethylene shade cloth was used in three of the comparisons and consisted of effective coverings of 100%, 60% with a silver reflective coating, and 60% black material with no reflective coating. Additionally, one of the structures was fitted with a poly snow fence with an effective shade of about 30%. Each shade structure contained a solar radiation meter and a black globe thermometer to measure radiant energy received under the shade material. Additionally, meteorological data were collected as a non-shaded treatment and included temperature, humidity, wind speed, and solar radiation. Data analyses was conducted using a physiological model based on temperature, humidity, solar radiation and wind speed; a second model using black globe temperatures, relative humidity, and wind speed was used as well. Analyses of the data revealed that time spent in the highest stress category was reduced by all shade materials. Moreover, significant differences (P < 0.05) existed between all shade materials (compared to no-shade) for hourly summaries during peak daylight hours and for ‘full sun’ days.     

Shade material evaluation using a cattle response model and meteorological instrumentation

Shade structures are often considered as one method of reducing stress in feedlot cattle. Selection of a suitable shade material can be difficult without data that quantify material effectiveness for stress reduction. A summer study was conducted during 2007 using instrumented shade structures in conjunction with meteorological measurements to estimate relative effectiveness of various shade materials. Shade structures were 3.6 m by 6.0 m by 3.0 m high at the peak and 2.0 m high at the sides. Polyethylene shade cloth was used in three of the comparisons and consisted of effective coverings of 100%, 60% with a silver reflective coating, and 60% black material with no reflective coating. Additionally, one of the structures was fitted with a poly snow fence with an effective shade of about 30%. Each shade structure contained a solar radiation meter and a black globe thermometer to measure radiant energy received under the shade material. Additionally, meteorological data were collected as a non-shaded treatment and included temperature, humidity, wind speed, and solar radiation. Data analyses was conducted using a physiological model based on temperature, humidity, solar radiation and wind speed; a second model using black globe temperatures, relative humidity, and wind speed was used as well. Analyses of the data revealed that time spent in the highest stress category was reduced by all shade materials. Moreover, significant differences (P < 0.05) existed between all shade materials (compared to no-shade) for hourly summaries during peak daylight hours and for ‘full sun’ days.     

Evaluation of d2, a microsatellite measure of inbreeding and outbreeding, in wolves with a known pedigree

We have evaluated a microsatellite measure proposed as an indicator of inbreeding and outbreeding using a captive wolf population with known inbreeding levels and founder sources. The measure, which is based on the difference in the repeat number for microsatellite alleles within an individual, was not more predictive of the known inbreeding coefficient than microsatellite heterozygosity (it was actually less predictive). We also found no support that the measure was predictive of the level of outbreeding. However, we could not determine if the measure was predictive of very low levels of inbreeding due to matings between remote relatives. Overall, it appears that the usefulness of this measure to identify individuals on the inbred-outbred continuum beyond that of heterozygosity and identify biologically important associations with fitness-related traits may be limited. We suggest that the measure be examined theoretically to determine when (and how much) the predictive value of the measure is different from that of heterozygosity for inbreeding or outbreeding levels in a variety of different scenarios.     

Shade material evaluation using a cattle response model and meteorological instrumentation

Shade structures are often considered as one method of reducing stress in feedlot cattle. Selection of a suitable shade material can be difficult without data that quantify material effectiveness for stress reduction. A summer study was conducted during 2007 using instrumented shade structures in conjunction with meteorological measurements to estimate relative effectiveness of various shade materials. Shade structures were 3.6 m × 6.0 m × 3.0 m high at the peak and 2.0 m high at the sides. Polyethylene shade cloth was used in three of the comparisons and consisted of effective coverings of 100%, 60% with a silver reflective coating, and 60% black material with no reflective coating. Additionally, one of the structures was fitted with a poly snow fence with an effective shade of about 30%. Each shade structure contained a solar radiation meter and a black globe thermometer to measure radiant energy received under the shade material. Additionally, meteorological data were collected as a non-shaded treatment and included temperature, humidity, wind speed, and solar radiation. Data analyses was conducted using a physiological model based on temperature, humidity, solar radiation and wind speed; a second model using black globe temperatures, relative humidity, and wind speed was used as well. Analyses of the data revealed that time spent in the highest stress category was reduced by all shade materials. Moreover, significant differences (P < 0.05) existed between all shade materials (compared to no-shade) for hourly summaries during peak daylight hours and for ‘full sun’ days.     

Potential bias in inbreeding depression estimates when using pedigree relationships to assess the degree of homozygosity for loci under selection

A potential bias in estimation of inbreeding depression when using pedigree relationships to assess the degree of homozygosity for loci under selection is indicated. A comparison of inbreeding coefficients based on either pedigree or genotypic frequencies indicated that, as a result of selection, the inbreeding coefficient based on pedigree might not correspond with the random drift of allelic frequencies. Apparent differences in average levels of both inbreeding coefficients were obtained depending on the genetic model (additive versus dominance, initial allelic frequencies, heritability) and the selection system assumed (no versus mass selection). In the absence of selection, allelic frequencies within a small population change over generations due to random drift, and the pedigree-based inbreeding coefficient gives a proper assessment of the accompanying probability of increased homozygosity within a replicate by indicating the variance of allelic frequencies over replicates. With selection, in addition to random drift, directional change in allelic frequencies is not accounted for by the pedigree-based inbreeding coefficient. This result implies that estimation of inbreeding depression for traits under either direct or indirect selection, estimated by a regression of performance on pedigree-based coefficients, should be carefully interpreted.Deceased     

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F_ROH) in three local dairy cattle breeds. F_ROH values were compared with F estimated from the genomic relationship matrix (F_GRM), based on the difference between observed v. expected number of homozygous genotypes (F_HOM) and the genomic homozygosity of individual i (F_{MOL i}). The molecular coancestry coefficient (f_{MOL ij}) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F_ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F_HOM and F_ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F_ROH and other estimated F coefficients were generally lower ranged from 0.45 (F_{MOL i}−F_ROH) in Cinisara to 0.17 (F_GRM−F_ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.     

An ABC estimate of pedigree error rate: application in dog, sheep and cattle breeds

On the basis of correlations between pairwise individual genealogical kinship coefficients and allele sharing distances computed from genotyping data, we propose an approximate Bayesian computation (ABC) approach to assess pedigree file reliability through gene-dropping simulations. We explore the features of the method using simulated data sets and show precision increases with the number of markers. An application is further made with five dog breeds, four sheep breeds and one cattle breed raised in France and displaying various characteristics and population sizes, using microsatellite or SNP markers. Depending on the breeds, pedigree error estimations range between 1% and 9% in dog breeds, 1% and 10% in sheep breeds and 4% in cattle breeds.