Fitting and validating the genomic evaluation model to Polish Holstein-Friesian cattle

The aim of the study was to fit the genomic evaluation model to Polish Holstein-Friesian dairy cattle. A training data set for the estimation of additive effects of single nucleotide polymorphisms (SNPs) consisted of 1227 Polish Holstein-Friesian bulls. Genotypes were obtained by the use of Illumina BovineSNP50 Genotyping BeadChip. Altogether 29 traits were considered: milk-, fat- and protein- yields, somatic cell score, four female fertility traits, and 21 traits describing conformation. The prediction of direct genomic values was based on a mixed model containing deregressed national proofs as a dependent variable and random SNP effects as independent variables. The correlations between direct genomic values and conventional estimated breeding values estimated for the whole data set were overall very high and varied between 0.98 for production traits and 0.78 for non return rates for cows. For the validation data set of 232 bulls the corresponding correlations were 0.38 for milk-, 0.37 for protein-, and 0.32 for fat yields, while the correlations between genomic enhanced breeding values and conventional estimated breeding values for the four traits were: 0.43, 0.44, 0.31, and 0.35. This model was able to pass the interbull validation criteria for genomic selection, which indicates that it is realistic to implement genomic selection in Polish Holstein-Friesian cattle.     

Genetic variability among Polish Red, Hereford and Holstein-Friesian cattle raised in Poland based on analysis of microsatellite DNA sequences

Polymorphism of 11 microsatellite DNA loci was analysed in Polish Red (PR), Hereford and Holstein-Friesian (HF) cattle raised in Poland and genetic distance among these breeds was determined. At the 11 loci (TGLA227, BM2113, TGLA53, ETH10, SPS115, TGLA126, TGLA122, INRA23, ETH3, ETH225 and BM1824) analysed with automated DNA sizing technology, a total of 213 alleles were identified: 76 in PR, 76 in HF, and 61 in Hereford. All the microsatellite DNA markers showed high polymorphism. Polymorphism information content (PIC) calculated for each marker exceeded 0.5, except for the ETH3 locus in Hereford cattle (PIC=0.475), and heterozygosity (H) ranged from 54.1% to as much as 85.2%. The coefficient of genetic distance was 0.354 between PR and Hereford, 0.414 between HF and Hereford, and 0.416 between PR and HF cattle.     

Identification of carriers of the mutation causing coagulation factor XI deficiency in Polish Holstein-Friesian cattle

Factor XI (FXI) deficiency is a hereditary coagulation disorder observed in various mammalian species. The molecular basis of coagulopathy has been recognized in Holstein cattle as a 76-bp insertion in the coding region of theFXI gene. Because the disorder seems to have an impact on reproductive traits and udder health in cattle, we tested 103 randomly selected cows, 28 cows with repeat breeding, and 9 cows with recurrent mastitis for the presence of an abnormalFXI allele. Three related cows were diagnosed as carriers.     

Effectiveness of a program aimed at the elimination of BLAD-carrier bulls from Polish Holstein-Friesian cattle

The molecular basis of BLAD is the D128G mutation of the gene coding for the CD18 subunit of beta-2 integrin. This mutation is lethal, since homozygous (BL/BL) animals die before they reach sexual maturity. In the 1990s, BLAD was the most widespread genetic disease in HF cattle worldwide. The aim of the present study was to determine the frequency of BLAD carriers among 4645 young breeding bulls in Poland in 1995-2006. The frequency of carriers of the mutated allele showed a clear decreasing trend. The highest frequency (7.9%) was recorded while implementing the BLAD control program (1995-1997). Regular monitoring has enabled a great reduction of this threat to the tested population. Today only sporadic cases of BL/TL heterozygotes are reported (ca. 0.8% in 2004-2006).     

Purging of inbreeding depression within the Irish Holstein-Friesian population

The objective of this study was to investigate whether inbreeding depression in milk production or fertility performance has been partially purged due to selection within the Irish Holstein-Friesian population. Classical, ancestral (i.e., the inbreeding of an individual''s ancestors according to two different formulae) and new inbreeding coefficients (i.e., part of the classical inbreeding coefficient that is not accounted for by ancestral inbreeding) were computed for all animals. The effect of each coefficient on 305-day milk, fat and protein yield as well as calving interval, age at first calving and survival to second lactation was investigated. Ancestral inbreeding accounting for all common ancestors in the pedigree had a positive effect on 305-day milk and protein yield, increasing yields by 4.85 kg and 0.12 kg, respectively. However, ancestral inbreeding accounting only for those common ancestors, which contribute to the classical inbreeding coefficient had a negative effect on all milk production traits decreasing 305-day milk, fat and protein yields by -8.85 kg, -0.53 kg and -0.33 kg, respectively. Classical, ancestral and new inbreeding generally had a detrimental effect on fertility and survival traits. From this study, it appears that Irish Holstein-Friesians have purged some of their genetic load for milk production through many years of selection based on production alone, while fertility, which has been less intensely selected for in the population demonstrates no evidence of purging.     

A mating method accounting for inbreeding and multi-trait selection in dairy cattle populations

Selection in dairy cattle populations usually takes into account both the breed profiles for many traits and their overall estimated breeding values (EBV). This can result in effective contributions of breeding animals departing substantially from contributions optimised for saving future genetic variability. In this work, we propose a mating method that considers not only inbreeding but also the detailed EBV of progeny or the EBV of sires in reference to acceptance thresholds. Penalties were defined for inbreeding and for inadequate EBV profiles. Relative reductions of penalties yielded by any mating design were expressed on a scale ranging from 0 to 1. A value of 0 represented the average performance of random matings and a value of 1 represented the maximal reduction allowed by a specialized, single-penalty, mating design. The core of the method was an adaptative simulated annealing, where the maximized function was the average of both ratios, under the constraints that both relative penalty reductions should be equal and that the within-herd concentration criterion should be equal to a predefined reasonable value. The method was tested on two French dairy cattle populations originating from the same AI organization. The optimised mating design allowed substantial reductions of penalty: 70% and 64% for the Holstein and the Normandy populations, respectively. Thus, this mating method decreased inbreeding and met various demands from breeders.     

A simple method to approximate gene content in large pedigree populations: application to the myostatin gene in dual-purpose Belgian Blue cattle

Gene content is the number of copies of a particular allele in a genotype of an animal. Gene content can be used to study additive gene action of candidate gene. Usually genotype data are available only for a part of population and for the rest gene contents have to be calculated based on typed relatives. Methods to calculate expected gene content for animals on large complex pedigrees are relatively complex. In this paper we proposed a practical method to calculate gene content using a linear regression. The method does not estimate genotype probabilities but these can be approximated from gene content assuming Hardy-Weinberg proportions. The approach was compared with other methods on multiple simulated data sets for real bovine pedigrees of 1 082 and 907 903 animals. Different allelic frequencies (0.4 and 0.2) and proportions of the missing genotypes (90, 70, and 50%) were considered in simulation. The simulation showed that the proposed method has similar capability to predict gene content as the iterative peeling method, however it requires less time and can be more practical for large pedigrees. The method was also applied to real data on the bovine myostatin locus on a large dual-purpose Belgian Blue pedigree of 235 133 animals. It was demonstrated that the proposed method can be easily adapted for particular pedigrees.     

PedNavigator: a pedigree drawing servlet for large and inbred populations

PedNavigator is a pedigree drawing application for large and complex pedigrees. It has been developed especially for genetic and epidemiological studies of isolated populations characterized by high inbreeding and multiple matrimonies. PedNavigator is written in Java and is intended as a server-side web application, allowing researchers to 'walk' through family ties by point-and-clicking on person's symbols. The application is able to enrich the pedigree drawings with genotypic and phenotypic information taken from the underlying relational database.     

Descriptive epidemiology of anestrus in Michigan Holstein-Friesian cattle

Twenty-two Michigan Holstein-Friesian herds were studied to determine the incidence and epidemiology of anestrus. In 3,309 lactations studied, 42% were classified as having exhibited preservice anestrus (no estrus detected by 70 d after calving). Organic reasons (pyometra, cystic follicles, static ovaries) were identified by palpation per rectum for 237 (23%) of the cows with preservice anestrus that were examined by a veterinarian. Postservice anestrus, defined as failure to show estrus within 35 d after an unsuccessful insemination, was identified in 790 (47%) of 1,691 lactations. Veterinary examination identified 104 (20%) of postservice anestrous cows as having this condition because of organic causes. The average cow with preservice anestrus had an increase of 30 d open compared to her herdmates, and the average cow with postservice anestrus had an increase of 37 d open. Anestrous cows produced more milk than their unaffected herdmates in both the current and the previous lactation. Analysis of composite lactation curves indicates that, for preservice anestrus, this additional milk production was obtained gradually over the entire lactation. Anestrous cows were culled at a significantly higher rate than their herdmates. Cows with preservice anestrus were more likely to have begun their lactation in the spring months.     

New method of estimating inbreeding in large semi-isolated populations with application to historic Britain

The purpose of this paper is to introduce a new method of estimating inbreeding in large, relatively isolated, populations over historic times, to demonstrate its application, and indicate some of its limitations and future developments. The method is based on the "paradox" of genealogy, and requires only that the variation of the population size be known, at least reasonably well, over an extended historic period. In this study a method has been developed to model this "paradox" which allows an estimation of the minimum level of inbreeding necessary for a given population curve in terms of values of Pearl's coefficients for each generation. As an example, the method is applied to the population of Britain. It is found that the frequency of siblings occurring in the same generation of a pedigree varies with the population size according to the Fermi-Dirac equation of statistical physics. The effect of introducing a single known estimate of inbreeding into the model is to make the otherwise diverse results for both the actual numbers of ancestors in a generation and the corresponding coefficients of inbreeding to converge.     

Evidence of adoption,monozygotic twinning,and low inbreeding rates in a large genetic pedigree of polar bears

Multigenerational pedigrees have been developed for free-ranging populations of many species, are frequently used to describe mating systems, and are used in studies of quantitative genetics. Here, we document the development of a 4449-individual pedigree for the Western Hudson Bay subpopulation of polar bears (Ursus maritimus), created from relationships inferred from field and genetic data collected over six generations of bears sampled between 1966 and 2011. Microsatellite genotypes for 22–25 loci were obtained for 2945 individuals, and parentage analysis was performed using the program FRANz, including additional offspring–dam associations known only from capture data. Parentage assignments for a subset of 859 individuals were confirmed using an independent medium-density set of single nucleotide polymorphisms. To account for unsampled males in our population, we performed half-sib–full-sib analysis to reconstruct males using the program COLONY, resulting in a final pedigree containing 2957 assigned maternities and 1861 assigned paternities with only one observed case of inbreeding between close relatives. During genotyping, we identified two independently captured 2-year-old males with identical genotypes at all 25 loci, showing—for the first time—a case of monozygotic twinning among polar bears. In addition, we documented six new cases of cub adoption, which we attribute to cub misidentification or misdirected maternal care by a female bereaved of her young. Importantly, none of these adoptions could be attributed to reduced female vigilance caused by immobilization to facilitate scientific handling, as has previously been suggested.     

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.     

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.     

Prediction of rates of inbreeding in selected populations

A method is presented for the prediction of rate of inbreeding for populations with discrete generations. The matrix of Wright's numerator relationships is partitioned into 'contribution' matrices which describe the contribution of the Mendelian sampling of genes of ancestors in a given generation to the relationship between individuals in later generations. These contributions stabilize with time and the value to which they stabilize is shown to be related to the asymptotic rate of inbreeding and therefore also the effective population size, Ne approximately 2N/(mu 2r + sigma 2r), where N is the number of individuals per generation and mu r and sigma 2r are the mean and variance of long-term relationships or long-term contributions. These stabilized values are then predicted using a recursive equation via the concept of selective advantage for populations with hierarchical mating structures undergoing mass selection. Account is taken of the change in genetic parameters as a consequence of selection and also the increasing 'competitiveness' of contemporaries as selection proceeds. Examples are given and predicted rates of inbreeding are compared to those calculated in simulations. For populations of 20 males and 20, 40, 100 or 200 females the rate of inbreeding was found to increase by as much as 75% over the rate of inbreeding in an unselected population depending on mating ratio, selection intensity and heritability of the selected trait. The prediction presented here estimated the rate of inbreeding usually within 5% of that calculated from simulation.     

Effects of inbreeding on the genetic diversity of populations

The study of variability within species is important to all biologists who use genetic markers. Since the discovery of molecular variability among normal individuals, data have been collected from a wide range of organisms, and it is important to understand the major factors affecting diversity levels and patterns. Comparisons of inbreeding and outcrossing populations can contribute to this understanding, and therefore studying plant populations is important, because related species often have different breeding systems. DNA sequence data are now starting to become available from suitable plant and animal populations, to measure and compare variability levels and test predictions.