PubMed Health⌕ Search

Biomedical subjects

P J Boettcher

Publications and source records attributed to P J Boettcher.

At least 19 recordsLinked to original sources

Genetic analysis of somatic cell scores in US Holsteins with a Bayesian mixture model.

The objective of this study was to apply finite mixture models to field data for somatic cell scores (SCS) for estimation of genetic parameters. Data were approximately 170,000 test-day records for SCS from first-parity Holstein cows in Wisconsin. Five different models of increasing level of complexity were fitted. Model 1 was the standard single-component model, and the others were 2-component Gaussian mixtures consisting of similar but distinct linear models. All mixture models (i.e., 2 to 5) included separate means for the 2 components. Model 2 assumed entirely homogeneous variances for both components. Models 3 and 4 assumed heterogeneous variances for either residual (model 3) or genetic and permanent environmental variances (model 4). Model 5 was the most complex, in which variances of all random effects were allowed to vary across components. A Bayesian approach was applied and Gibbs sampling was used to obtain posterior estimates. Five chains of 205,000 cycles were generated for each model. Estimates of variance components were based on posterior means. Models were compared by use of the deviance information criterion. Based on the deviance information criterion, all mixture models were superior to the linear model for analysis of SCS. The best model was one in which genetic and PE variances were heterogeneous, but residual variances were homogeneous. The genetic analysis suggested that SCS in healthy and infected cattle are different traits, because the genetic correlation between SCS in the 2 components of 0.13 was significantly different from unity.

Animals↗

Molecular typing of Staphylococcus aureus isolated from cows, goats and sheep with intramammary infections on the basis of gene polymorphisms and toxins genes.

We investigated 116 Staphylococcus aureus isolates from cows, goats and sheep with intramammary infections (IMI) in Italy to provide information about the spread of enterotoxigenic strains and to compare strains isolated from different ruminant species. The isolates were typed by restriction fragment length polymorphism (RFLP) analysis of the coagulase (coa) gene, by analysis of polymorphisms of the X region of protein A (spa) gene and by detection of genes encoding enterotoxins (sea, seb, sec, sed, see, seg, seh, sei, sej and sel). Seven different coa types and 12 different spa types were distinguished. On the basis of polymerase chain reaction-RFLP, 29 different coa subtypes were identified. Two different coa subtypes accounted for 49% and 67% of bovine and ovine isolates respectively. Only seven coa subtypes were observed in isolates from more than one host species and no coa subtype was present in isolates from all three ruminant species. Furthermore, 85 of the isolates (73%) harboured at least one enterotoxin gene (se) with a predominance of sea, sed and sej among isolates from bovine IMI, and sec and sel among isolates from caprine and ovine IMI. Comparing the S. aureus isolates on the basis of gene polymorphisms and presence of se genes, significant differences were found in distributions of genotypes among isolates from cows, goats and sheep.

Animals↗

Relationships between somatic cell count and intramammary infection in buffaloes.

The objectives of this study were to evaluate the presence of intramammary infections (IMI) in dairy buffaloes and to examine the relationships among IMI, somatic cell counts (SCC), and milk production traits. Two farms in northern Italy were visited monthly for a complete milking season. Quarter-based milk samples were collected at each visit from 46 buffaloes. A total of 1,912 samples were assessed in this experiment. Samples were cultured for bacterial presence and were tested for SCC and percentages of milk protein and fat. In addition, daily milk yield was recorded from each buffalo. Prevalence of IMI was large; 63% of quarters were infected. No buffalo remained free from IMI throughout the course of the study. Coagulase-negative staphylococci were the most common pathogen (66% of positive samples). The SCC was distinctly greater in infected quarters; 100% of quarters with SCC >200,000 cell/mL had IMI, whereas 98% of quarters with SCC below this threshold were uninfected. The somatic cell scores (SCS) in these buffaloes were much lower than those commonly observed in dairy cattle. The mean SCS from quarters with IMI was only 2.93. The highest SCS was observed in quarters infected by streptococci. No drastic decrease in milk yield was observed among infected buffaloes relative to healthy contemporaries. The relatively low SCS and lack of a strong effect on milk yield provide evidence to discourage antibiotic treatment of buffaloes for subclinical IMI during lactation.

Animals↗

Risk factors for intramammary infections and relationship with somatic-cell counts in Italian dairy goats.

Routine examination of milk was performed on five herds of lactating goats in northern Italy as part of a milk quality-monitoring program in the year 2000. As part of the study, aseptic samples of foremilk were collected monthly from both half udders during the entire lactation for 305 goats, resulting in a total of 4571 samples. The samples were tested with cytological and bacteriological analyses to evaluate the relationship between mammary infections and somatic-cell count (SCC; Fossomatic (TM) method). Prevalence of intramammary infection (IMI) was 40.2% (n = 1837) of all udder-half samples examined. The most-prevalent mastitis agents were coagulase-negative Staphylococci (CNS), 80% (n = 1474 udder-half samples); within this group, Staphylococcus epidermidis was the most-prevalent species (38%). Other prevalence were Staphylococcus aureus 6% (n = 112 udder-half samples) and environmental pathogens 14% of infected udder-half samples (n = 251) with a diverse mixture of species, none of which had a frequency of > 4%. Enterococcus faecalis was the most-frequently isolated among this group. Neither Salmonella spp. nor Listeria monocytogenes were detected. The risk (sample level) of infection differed across herds, parities, and stage of lactation according to results from logistic multiple regression. Infection was more common among goats in third and fourth parities and during the later stages of lactation. Of the 2734 samples from uninfected udder halves, the mean log2 SCC was 3.9 cell/ml; of the 1837 bacteriological positive samples, the mean log2 SCC was 5.6 cell/ml. According to results from a linear mixed model, concentrations of somatic cells tended to increase with increasing age and days in milk and with the presence of bacteria. Infection with S. aureus was associated with the highest SCS.

Animals↗

Application of a finite mixture model to somatic cell scores of Italian goats.

The objectives of this study were to apply a finite mixture model (FMM) to data for somatic cell count in goats and to compare the fit of the FMM with that of a standard linear mixed effects model. Bacteriological information was used to assess the ability of the model to classify records from healthy or infected goats. Data were 4518 observations of somatic cell score (SCS) and bacterial infection from both udder halves of 310 goats from 5 herds in Northern Italy. The records were from a complete production season, and were taken monthly from February to November 2000. Explanatory factors in both models included a 3-parameter regression on days in milk (DIM); fixed class effects of herd-test-day, parity group, and udder side (left or right); and random effects of goat and udder half within goat. In addition, the 2-component FMM included a fixed mean for the second component of the model (theoretically corresponding to infected udder halves), as well as an unknown probability of membership to a given putative infection status. A Bayesian statistical approach was used for the analysis with Gibbs sampling used to obtain draws from posterior distributions of parameters of interest. Two sampling chains of 200,000 cycles each were generated for each model. The FMM yielded a much lower estimate of residual variance than the standard model (1.28 vs. 3.02 SCS2), and a slightly higher estimate for the between-goat variance (1.79 vs. 1.48). The deviance information criterion (DIC) was used to compare the fit of the 2 models. The DIC was much lower for the FMM, indicating a better fit to the data. The FMM was able to classify correctly 60 and 48% of the healthy and infected observations, respectively. This was slightly higher than what would be expected from random classification, but not high enough for useful mastitis diagnosis. Nevertheless, increased precision of genetic evaluation is the goal of applying the FMM, rather than timely and accurate mastitis diagnosis. The results suggest that more research on FMM for SCS is merited and necessary for proper application.

Analysis of Variance↗

Genotype x environment interaction for grazing versus confinement. I. Production traits.

The objective of this study was to investigate the possible existence of a genotype x environment interaction (GxE) for production traits of US Holsteins in grazing versus confinement herds. Grazing herds were defined as those that utilized grazing for at least 6 mo and were enrolled in dairy herd improvement (DHI). Control herds were confinement DHI herds of comparable size in similar regions. The performance of daughters in grazing herds and control herds was examined using linear regression of mature equivalent milk, fat, and protein yield on the November 2000 USDA-DHI predicted transmitting abilities (PTA) of their sires for those traits. Heritabilities and genetic correlations were estimated using restricted maximum likelihood in a bivariate animal model that considered the same trait in different environments as different traits. Product-moment and rank correlations were calculated between sires' estimated breeding values, estimated separately in both environments. For grazing herds, the coefficient of regression of milk, fat and protein on PTA were 0.78, 0.76, and 0.78, respectively. Corresponding coefficients in the control herds were 0.99, 0.96, and 0.98. Estimates of heritability for the traits ranged from 0.2 to 0.25, and differences between grazing and control environments were small. Estimates of the genetic correlations for the traits in both environments were 0.89, 0.88, and 0.91 for milk, fat, and protein, respectively. Within-quartile analyses revealed a lower correlation for milk and protein between the upper and lower grazing quartiles, while the same quartiles for the control herds did not differ from unity. Rank correlation coefficients between sire estimated breeding values from the 2 environments were 0.59, 0.63, and 0.66 for milk, fat, and protein, respectively. The mean rank change for the top 100 sires between the two environments was 27. The regression coefficients indicate that expected daughter differences may be overstated by current sire PTA in grazing herds. Genetic correlations less than unity suggests that there is, at least, some reranking among sires in both environments, while the rank correlations indicate the possibility of sire reranking when evaluations were performed within management system. However, differences are not so large as to justify separate genetic evaluations for each system.

Animals↗

Genotype x environment interaction for grazing vs. confinement. II. Health and reproduction traits.

Continual selection for increased milk yield for more than 40 yr, combined with the antagonistic relationship between increasing yield, somatic cell count, and fertility, have resulted in sires that may not be optimal for producing daughters for grazing systems where seasonal calving is very important. The objective of this study was to investigate the possible existence of a genotype x environment interaction (G x E) in grazing vs. confinement herds within the United States for lactation average somatic cell score (LSCS), days open (DO), days to first service (DFS), and number of services per conception (SPC). Grazing herds were defined as those that utilized grazing for at least 6 mo each year and were enrolled in Dairy Herd Improvement (DHI). Control herds were confinement DHI herds of similar size in the same states. For LSCS, the performance of daughters in grazing and control herds was examined using linear regression of LSCS on the November 2000 USDA-DHIA sire predicted transmitting abilities (PTA) for SCS. Genetic parameters for all traits were estimated using REML in a bivariate animal model that treated the same trait in different environments as different traits. Rank correlations were calculated between sires' estimated breeding values for LSCS, calculated separately for sires in both environments. The coefficient of regression of daughter LSCS on sire PTA was less in grazing herds than in control herds. The coefficient of regression for control herds was closer to expectation. Estimates of heritability were approximately 0.12 for LSCS, and less than 0.05 for the reproduction traits. Heritabilities for DO, DFS, and SPC were slightly higher for control herds. Estimates of genetic correlation for each reproductive trait between the 2 environments were high and not significantly different from unity. Generally, these traits appear to be under similar genetic control, but a lower coefficient of regression of LSCS on sire PTA for SCS in grazing herds suggests differences in daughter performance in grazing herds may be overstated based on current PTA for SCS.

Animals↗

Phenotypic and genetic influences on test-day measures of acetone concentration in milk.

The objectives of this study were to estimate heritability of acetone concentration in milk, based on monthly samples of milk obtained as part of a routine milk testing program, and to evaluate the feasibility of using such data in a genetic evaluation program for selection against ketosis incidence. Milk samples were collected from January to December of 1999 in herds enrolled in the Ontario Dairy Herd Improvement Association, and acetone concentration was measured using an inline chemical procedure. The original data included more than 50,000 records. Because ketosis is generally a problem during early lactation, only the single test with the fewest days in milk was retained. In addition, data were retained only for cows with pedigree information. The final data set included 10,375 records. Among these data, only 6.56% had detectable levels of acetone. Acetone data were log-transformed prior to statistical analysis. Simple ANOVA indicated that herd, parity number, days in milk, and month of test had significant effects on acetone concentration. Acetone levels increased with lactation number and were higher in early lactation. Three approaches were applied for genetic analysis. First, REML was used with a simple linear animal model. Then, a separate procedure used data augmentation and Gibbs Sampling to obtain continuously distributed underlying values for records with zero acetone concentration, and these data were analyzed with both an animal and sire model. Heritability of acetone concentration was less than 1% for all 3 analyses. Herd effects accounted for about 5% of the phenotypic variance. Low estimates of heritability were due either to low actual levels of genetic variance or inability to detect all of the genetic variance present, due to infrequent recording during the early part of the lactation. Genetic evaluation based on recording of acetone concentration on a monthly basis seems of little use as a selection tool to decrease incidence of ketosis.

Acetone↗

A Monte Carlo approach for estimation of haplotype probabilities in half-sib families.

The objective of this work was to propose an algorithm (HAPROB) to estimate haplotype probabilities for genotyped members of half-sib families for which parents lacked genotypic information. The algorithm had 2 basic steps. First, a Monte Carlo-based approach was used to estimate haplotype probabilities for sires conditional upon offspring genotypes and population allelic frequencies, and then offspring-haplotype probabilities were estimated conditional upon sire probabilities and population frequencies. The 2 steps were alternated iteratively until estimates of population frequencies were essentially unchanged. Simulation was used to evaluate effects of the number of Monte Carlo cycles on the accuracy of the reconstructed haplotypes. Fifty thousand cycles was found to be sufficient for the haplotype configurations considered. Accuracy of the algorithm was compared with that obtained by the public domain SIMWALK2 software. Predictions of the most likely haplotype configurations are produced by SIM-WALK2, but no estimates of probability are given. The accuracy of the current approach was comparable to that obtained from SIMWALK2. The proportions of times that haplotypes were reconstructed correctly were 87.0 and 92.4% (sires and offspring) for HAPROB vs. 87.5 and 91.5% for SIMWALK2. Effects of family size on accuracy of reconstruction were examined. Accuracy of reconstruction was only about 4% for sires with 2 offspring, but accuracy among the offspring themselves was 65%. Accuracy increased quickly as family size increased and reached 100% for sires with 30 offspring. Maximum accuracy for offspring was about 96%. Estimates of haplotype probabilities produced can be used in regression analyses to estimate effects of haplotypes on quantitative phenotypes.

Algorithms↗

Effects of casein haplotypes on milk production traits in Italian Holstein and Brown Swiss cattle.

The objective of this study was to estimate the effects of different haplotypes of the casein genes on milk production traits in Italian dairy cattle. Traits of interest were yields of milk, fat, and protein, and percentages of fat and protein in milk. The data included 728 multiparous records from 347 Holsteins and 773 records from 298 Brown Swiss cows. Records were preadjusted for effects of age and parity, season of calving, and region, and expressed as deviations from herdmate averages. Twenty half-sib families were represented in each breed. Haplotype probabilities were estimated for each animal and phenotypes were regressed on these probabilities. Nine haplotypes were observed in Holsteins and 17 were identified among the Brown Swiss. For Holsteins, significant effects were observed for protein percentage, with some indication of an effect for fat percentage. For the Brown Swiss, effects of haplotypes were significant for milk yield and fat and protein percentages. Effects were strongest for protein percentage. Correlation coefficients of solutions across breeds tended to be strong and positive, indicating that the same haplotypes had similar estimated effects in the 2 breeds. Although the data were limited (<350 cows in each study), this latter result may suggest that genes in the casein complex itself are responsible for the effects observed, rather than loci that are physically linked on either side of the casein cluster.

Animals↗

Estimation of variances for gametic effects on litter size in Yorkshire and Landrace swine.

The objective of this study was to test for effects of gametic imprinting on litter size in swine by estimating variances for parent-specific gametic effects. Data were 64,047 and 137,009 multiparous records of number born alive for the U.S. Landrace and Yorkshire breeds, respectively. The statistical model included fixed effects of parity number and herd, and random effects of herd-year-season, mate, permanent environment, animal (additive genetic), and either maternal or paternal gametes. A Bayesian approach that used Gibbs sampling to obtain posterior distributions was employed. To aid in the interpretation of results, the Landrace data structure was used to simulate data with and without effects of imprinting. Analyses of the simulated records indicated that the model applied was capable of detecting effects of imprinting when such effects were present. Small, but non-zero, estimates of gametic variances were obtained when no imprinting was simulated. Estimates of the proportion of total variance accounted for by paternally transmitted gametes were 0.8 and 0.9% for Landrace and Yorkshires, respectively. These estimates were different from zero, but were similar to the results observed for data simulated without an imprinting effect. Corresponding results for maternally transmitted gametes were 1.6% for Landrace and 0.8% for Yorkshires. The estimate for Landrace was significantly greater than that observed for Yorkshires and for the simulations without a true effect and suggested the presence of a non-Mendelian genetic influence on litter size. Paternally imprinted genes are a plausible reason for the observed results. Assuming that the effect observed was due to paternal imprinting at a single biallelic locus, the substitution effect of the superior allele could be greater than 0.7 piglets per litter. Identification of a genetic marker for such an allele would be useful in marker-assisted selection of females. Other possible explanations exist for the increased gametic variance in the Landrace breed, but these explanations (such as maternal or cytoplasmic effects) may be less likely than paternal imprinting.

Animals↗

The relationship between milk production and antibody response to ovalbumin during the peripartum period.

Suboptimal innate and immune mechanisms of host resistance during the peripartum period may contribute to increased incidence of mastitis. To evaluate associations between antibody response to ovalbumin and milk production during the peripartum period, 136 Holstein cows and heifers from three herds with known antibody response profiles, were evaluated for projected 305-d milk, protein, and fat yield. Using a previously described index (Wagter et al., 2000), cows were quantitatively classified based on their profile of antibody response to ovalbumin into high, average, or low antibody response groups. The single-effect antibody response group contributed significantly to variation in fat and protein yield, but not milk yield. The interaction between antibody response and parity significantly contributed to the variation in milk, fat, and protein yields; therefore the effects of group were reported on a within-parity basis. Among first-parity cows, low responders had a higher fat and protein yield than high or average antibody responder animals. Among older cows (parity 3 or greater) milk yield was significantly higher for those in the high antibody response group compared with average and low response groups. However, no significant differences in fat or protein yields were observed between high and low antibody response groups. These results suggest the possibility to select cows for enhanced immune response with no adverse effects on yield. That first-parity cows with low antibody response produce more fat and protein may be offset by the fact that mastitis occurrence was highest in this group in two out of three herds investigated. Selection for high immune response may prove beneficial to herd life by maintaining optimal yield, yet minimizing occurrence of disease.

Animals↗

Genotype x environment interactions in conventional versus pasture-based dairies in Canada.

The objective of this study was to evaluate effects of genotype x environmental interaction (GxE) on dairy traits among herds in Canada that practice intensive rotational grazing versus conventional methods that rely on stored feeds. Based on responses to questionnaires, 22 herds were selected for the grazing group and the conventional group consisted of 34 large free-stall dairies. Production data consisted of 6,749 lactations of 2,817 cows in the grazing herds and 29,371 lactations of 12,774 cows in the control herds. Subsets of data for reproduction and conformation also were created. Per cow production was greater in the conventionally managed herds (9,947 kg of milk) than in the grazing herds (9,400 kg). Phenotypic data were regressed on sire predicted transmitting abilities (PTA) and genetic correlations for all traits were estimated between environments. A scaling effect was observed across environments for yield traits, as phenotypic differences among cows were greater than predicted by sire PTA in conventional herds and consistent with sire PTA in grazing herds. Heritabilities for yield also were greater in conventional herds. Other effects of GxE were minor, with genetic correlations of near 0.90 or greater for all traits. Among yield traits, genetic correlation was lowest for fat (0.88 with SE = 0.04). These results demonstrate that effects of GxE are small between the two types of management systems analyzed and that graziers can accurately select sires based on national rankings.

Animal Feed↗

Genetic parameters for feet and leg traits evaluated in different environments.

The objective of this study was to test for genotype x environmental interaction (GXE) for feet and leg traits scored in different environments. Genetic correlations of seven feet and leg traits were estimated across different management systems: free versus tie stalls, slatted versus solid flooring, and intact versus trimmed hooves. Data were records from first-lactation Holstein cattle. Traits were claw uniformity, depth of heel, rear leg rear view, foot angle, bone quality, rear leg side view, and overall feet and legs. Different subsets of data were used for each comparison, resulting in 147,400; 53,550; and 145,160 records for housing, flooring, and hoof trimming management systems, respectively. Genetic parameters were estimated using REML and two-trait models in which for each animal a given trait was observed in one environment and missing in the other. Phenotypic scores were lower with tie stalls, slatted floors, and no trimming. Heritabilities tended to be greater in herds with tie stalls and slatted floors. Trimming had little effect on genetic parameters. The genetic correlations of feet and leg traits across pairs of management systems were > or = 0.85, except for rear legs, rear view. Therefore, effects of GXE were assumed to be of little importance and modification of genetic evaluation procedures on the basis of housing, flooring, and hoof conditions seems unnecessary.

Animal Husbandry↗

Estimation of genetic parameters for concentrations of milk urea nitrogen.

The objective of this study was to use field data collected by dairy herd improvement programs to estimate genetic parameters for concentrations of milk urea nitrogen (MUN). Edited data were 36,074 test-day records of MUN and yields of milk, fat, and protein obtained from 6102 cows in Holstein herds in Ontario, Canada. Data were divided into three sets, for the first three lactations. Two analyses were performed on data from each lactation. The first procedure used ANOVA to estimate the significance of the effects of several environmental factors on MUN. Herd-test-day effects had the most significant impact on MUN. Effects of stage of lactation were also important, and MUN levels tended to increase from the time of peak yield until the end of lactation. The second analysis used a random regression model to estimate heritabilities and genetic correlations of MUN and the yield traits. Heritability estimates for MUN in lactations one, two, and three were 0.44, 0.59, and 0.48, respectively. Heritabilities for the yield traits were of a similar magnitude. Little relationship was observed between MUN and yield. Raw phenotypic correlations were all <0.10 (absolute value). Genetic correlations with production traits were close to zero in lactations one and three and only slightly positive in lactation two. The results indicate that selection on MUN is possible, but relationships between MUN and other economically important traits such as metabolic disease and fertility are needed.

Animals↗

Evaluation of sire predicted transmitting abilities for evidence of X-chromosomal inheritance in north american sire families.

This study tested for differences between genetic merits of sons and daughters of sires and for evidence of segregating quantitative trait loci on the X chromosomes of North American Holsteins. Son PTA adjusted for sire PTA was used as the dependent variable to test for biases and for genes that were passed from sire to daughter but not to son. The test of variability across sires of sons merely indicated an unaccounted source of variation, for which genes on X chromosomes might be responsible. Critical values for this test and power were determined by simulation for a variety of populations and traits differing in heritability, size of the X chromosome effect, and allelic frequency. Simulated genes on the X chromosome were detected with high power at intermediate frequencies of the favorable allele. The power of the test increased as the size of the effect increased and as genetic variance attributed to autosomes decreased. The test was then applied to recently evaluated data from US and Canadian Holstein populations. Genetic evaluations for >17,000 bulls from the US and >9000 from Canada were included. Results suggested that little extra variation was present for some traits formally evaluated in North America, but that genes on the X chromosome were unlikely to be the cause.

Animals↗

Bayesian inference for categorical traits with an application to variance component estimation.

We implemented statistical models of Bayesian inference that included direct and maternal genetic effects for genetic parameter estimation of categorical traits by Gibbs sampling. The estimation errors and variances of estimates of animal versus sire and maternal grandsire models, of linear versus threshold models, of single-trait versus multiple-trait models, and of treating herd-year-season as fixed versus random effects in the model were compared. The results indicated that linear models yielded biased estimates of genetic parameters for categorical traits. The animal model was improper for analysis of categorical traits using a threshold model and the Gibbs sampler. Moreover, linear versus threshold models and animal versus sire-maternal grandsire models resulted in larger Monte Carlo errors and increased auto-correlations among posterior samples. Treating herd-year-seasons as random effects in the threshold models decreased the Monte Carlo error, auto-correlations, and the variances of estimates. Efficiency of the single-trait threshold sire model, as measured by the variance of the estimates, was lower than for a multiple-trait model that included a correlated continuous trait, but both estimates were unbiased. Therefore, the threshold single-trait sire and maternal grandsire model is a feasible alternative to the multiple-trait model for analysis of variance components of categorical traits affected by direct and maternal genetic factors.

Animals↗

Bayesian segregation analysis of somatic cell scores of Ontario Holstein cattle.

Bayesian segregation analysis using a Gibbs sampling approach was applied to four sets of simulated data and one set of field data to detect evidence of major genes affecting the evaluated trait. The substitution effect of a major gene and its allelic frequency were estimated for each set of data. For two datasets simulated with a model with no major gene effect, the resulting estimates of polygenic variance and heritability agreed with the simulated values and tests for the presence of a major gene were not significant. Analyses of two sets of data simulated with a major gene produced posterior distributions that gave significant evidence of major gene effects but underestimated the substitution values of the major gene. The segregation analysis of field data suggested that a major gene significantly affected somatic cell score (SCS) in the population of Ontario Holstein cattle. The estimated heritability of SCS was approximately 0.16. The major gene variance accounted for about 17% of the total genetic variance and the point estimate of the frequency of the allele having a positive effect on SCS was 0.30. However, the precision of these estimates is questionable based on the simulation results. The effect of the major gene may be underestimated.

Animals↗