Whole‐Genome association analysis of susceptibility to paratuberculosis in Holstein cattle

The objective of this study was to identify genetic markers and genomic regions associated with susceptibility to Mycobacterium avium ssp. paratuberculosis (MAP) infection in Holstein cattle. Associated single nucleotide polymorphisms (SNPs) were identified by genotyping 521 MAP‐infected Holstein cows and comparing SNP allele frequencies of these infected cows with allele frequencies estimated from specific reference populations. Reference population allele frequency estimates used Holstein sire genotype data and were weighted estimates based on sire usage within the population in question. The 521 infected cows were 233 and 288 cows from two resource populations of approximately 5000 cows each, collected independently. Population 1 was comprised primarily of daughters of twelve Holstein artificial insemination sires used heavily within the US dairy cattle population. Samples were obtained from 300 co‐operating commercial dairy herds throughout the US and were tested by both MAP faecal culture and blood‐enzyme‐linked immunosorbent assay (ELISA). Population 2 consisted of dairy cattle from six co‐operating dairy herds in Wisconsin, with all animals in the herds tested by blood enzyme‐linked immunosorbent assay (ELISA) for MAP infection. Genotyping was performed with the Illumina Bovine SNP50 Bead Chip, providing genotypes for 35 772 informative SNPs. Data from the two resource populations were analysed both in separate and combined analyses. The most significant autosomal markers from the individual and combined analyses (n = 197, nominal P < 0.001) were used in a stepwise logistic regression analysis to identify a set of 51 SNPs that could be used as a predictor of genetics for Holstein cattle susceptibility to MAP infection.     

Quantitative trait loci affecting lactose and total solids on chromosome 6 in Brazilian Gir dairy cattle

Fourteen Brazilian Gir sire families with 657 daughters were analyzed for quantitative trait loci (QTL) on chromosome 6 affecting lactose and total solids. Cows and sires were genotyped with 27 microsatellites with a mean spacing between markers of 4.9 cM. We used a 1% chromosome-wide threshold for QTL qualification. A QTL for lactose yield was found close to marker MNB66 in three families. A QTL for total solid yield was identified close to marker BMS2508 in three families. A QTL for lactose percentage, close to marker DIK1182, was identified in two families. A QTL for total solid percentage, close to marker MNB208, was identified in four families. These QTLs could be used for selection of animals in dairy production systems.     

Quantitative trait locus mapping in dairy cattle by means of selective milk DNA pooling using dinucleotide microsatellite markers: analysis of milk protein percentage.

"Selective DNA pooling" accomplishes quantitative trait locus (QTL) mapping through densitometric estimates of marker allele frequencies in pooled DNA samples of phenotypically extreme individuals. With poly(TG) microsatellites, such estimates are confounded by "shadow" ("stutter") bands. A correction procedure was developed on the basis of an observed linear regression between shadow band intensity and allele TG repeat number. Using this procedure, a selective DNA pooling study with respect to milk protein percentage was implemented in Israel-Holstein dairy cattle. Pools were prepared from milk samples of high and low daughters of each of seven sires and genotyped with respect to 11 markers. Highly significant associations with milk protein percentage were found for 5 of the markers; 4 of these markers confirmed previous reports. Selective DNA pooling accessed 80.6 and 48.3%, respectively, of the information that would have been available through individual selective genotyping or total population genotyping. In effect, the statistical power of 45,600 individual genotypings was obtained from 328 pool genotypings. This methodology can make genome-wide mapping of QTL accessible to moderately sized breeding organizations.     

Identification of an ovulation rate QTL in cattle on BTA14 using selective DNA pooling and interval mapping

Increased twinning incidence in beef cattle has the potential to improve production efficiency. However, phenotypic selection for twinning rate is difficult because of the trait's low heritability and the long time interval necessary to collect phenotypic records. Therefore, this trait and the correlated trait of ovulation rate are ideal candidates for marker-assisted selection. The objective of this study was to complete a genome-wide search for ovulation rate quantitative trait loci (QTL) in two related sire families. The families (paternal halfsib sires 839802 and 839803) were from a population of cattle selected for ovulation rate at the USDA Meat Animal Research Center, Clay Center, Nebraska. Putative ovulation rate QTL have previously been identified in the 839802 family on chromosomes 7 and 19; however, marker coverage in the original scan was not complete. This study fills the gaps in marker coverage of the earlier study by adding approximately 60 informative microsatellites to each sire family. Each family was genotyped using selective DNA pooling. Sons and daughters were included in either the high or low pool based on their estimated breeding value deviations from the mid-parent average (EBVMD) for ovulation rate. Approximately 40% (839802) and 26% (839803) of available progeny comprised the high and low pools combined. Pooled typing revealed possible associations (nominal P < 0.05) between ovulation rate and marker genotype for 11 and 15 microsatellites in the 839802 and 839803 families, respectively. Subsequent interval mapping strengthened support for the presence of an ovulation rate QTL on BTA14 (chromosome-wise P < 0.02).     

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design

We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.     

A whole genome scan for quantitative trait loci affecting milk protein percentage in Israeli-Holstein cattle, by means of selective milk DNA pooling in a daughter design, using an adjusted false discovery rate criterion

Selective DNA pooling was employed in a daughter design to screen all bovine autosomes for quantitative trait loci (QTL) affecting estimated breeding value for milk protein percentage (EBVP%). Milk pools prepared from high and low daughters of each of seven sires were genotyped for 138 dinucleotide microsatellites. Shadow-corrected estimates of sire allele frequencies were compared between high and low pools. An adjusted false discovery rate (FDR) method was employed to calculate experimentwise significance levels and empirical power. Significant associations with milk protein percentage were found for 61 of the markers (adjusted FDR = 0.10; estimated power, 0.68). The significant markers appear to be linked to 19--28 QTL. Mean allele substitution effects of the putative QTL averaged 0.016 (0.009--0.028) in units of the within-sire family standard deviation of EBVP% and summed to 0.460 EBVP%. Overall QTL heterozygosity was 0.40. The identified QTL appear to account for all of the variation in EBVP% in the population. Through use of selective DNA pooling, 4400 pool data points provided the statistical power of 600,000 individual data points.     

Optimal designs for linkage disequilibrium mapping and candidate gene association tests in livestock populations

Simulation was used to evaluate the performance of different selective genotyping strategies when using linkage disequilibrium across large half-sib families to position a QTL within a previously defined genomic region. Strategies examined included standard selective genotyping and different approaches of discordant and concordant sib selection applied to arbitrary or selected families. Strategies were compared as a function of effect and frequency of QTL alleles, heritability, and phenotypic expression of the trait. Large half-sib families were simulated for 100 generations and 2% of the population was genotyped in the final generation. Simple ANOVA was applied and the marker with the greatest F-value was considered the most likely QTL position. For traits with continuous phenotypes, genotyping the most divergent pairs of half-sibs from all families was the best strategy in general, but standard selective genotyping was somewhat more precise when heritability was low. When the phenotype was distributed in ordered categories, discordant sib selection was the optimal approach for positioning QTL for traits with high heritability and concordant sib selection was the best approach when genetic effects were small. Genotyping of a few selected sibs from many families was generally more efficient than genotyping many individuals from a few highly selected sires.     

Genetic parameters for paratuberculosis infection and effect of infection on production traits in Israeli Holsteins

Paratuberculosis (Johne's disease) is an infectious enteric disease in dairy cattle and other species that causes substantial economic loss worldwide. In this study, two recursive Gaussian-threshold models were employed in order to infer the effects of Johne's disease on milk yield, fat yield, and protein yield while simultaneously estimating genetic parameters (i.e. heritability and genetic correlation) in an Israeli Holstein population. Disease diagnosis was based on ELISA serum antibody tests. Data were available for 4694 daughters of 361 sires; 3.5% were positive; and 1.6% were suspect for the disease test. Disease status was coded either as a binary character (negative vs. positive) or as an ordered categorical trait (negative, suspect, and positive) in the two recursive models and as a binary trait in a linear model. Among sires with ≥ 50 daughters, predicted probability of Mycobacterium avium ssp. paratuberculosis (MAP) infection in future daughters ranged from <1% to 16.5%. Heritability estimates for Johne's disease were near 0.15, confirming a genetic contribution to disease susceptibility. Genetic correlation estimates for Johne's disease with the three yield traits were 0.15-0.22. Residual correlations for Johne's disease with the yield traits were between -0.01 and -0.10. For the linear regression model, yield losses associated with a positive disease diagnosis during 305 days of lactation were 300 kg milk and around 10 kg for fat and protein. Yield loss estimates from the recursive models were 25-50% less than linear model estimates. Recursive modeling has theoretical advantages over linear models for these phenotypes, but the estimated genetic parameters in this study did not differ significantly between the two types of models.     

IS900 PCR to detect Mycobacterium paratuberculosis in retail supplies of whole pasteurized cows' milk in England and Wales.

IS900 PCR for Mycobacterium paratuberculosis was applied to cream, whey, and pellet fractions of centrifuged whole cows' milk. The test and simultaneous control reactions gave correct results for spiked milk and for native milk samples obtained directly from M. paratuberculosis-free, subclinically infected, and clinically infected cows. The test was then applied to units of whole pasteurized cows' milk widely obtained from retail outlets throughout central and southern England from September 1991 to March 1993. With peak periods in January to March and in September to November, when up to 25% of units were affected, an overall 22 of 312 samples (7%) tested positive for M. paratuberculosis. In 18 of the 22 positive samples (81%), the PCR signal segregated to the cream or pellet fractions or both, consistent with the presence of intact mycobacteria. Nine of 18 PCR-positive milk samples (50%) and 6 of 36 PCR-negative milk samples (16%) yielded long-term liquid cultures which tested positive for M. paratuberculosis after 13 to 40 months of incubation, despite overgrowth by other organisms. Taken together with data on the prevalence of M. paratuberculosis infection in herds in the United Kingdom, the known secretion of M. paratuberculosis in milk from subclinically infected animals, and the inability of laboratory conditions simulating pasteurization to ensure the killing of all these slowly growing or unculturable organisms, there is a high risk, particularly at peak times, that residual M. paratuberculosis will be present in retail pasteurized cows' milk in England.     

Modelling QTL effect on BTA06 using random regression test day models

In statistical models, a quantitative trait locus (QTL) effect has been incorporated either as a fixed or as a random term, but, up to now, it has been mainly considered as a time-independent variable. However, for traits recorded repeatedly, it is very interesting to investigate the variation of QTL over time. The major goal of this study was to estimate the position and effect of QTL for milk, fat, protein yields and for somatic cell score based on test day records, while testing whether the effects are constant or variable throughout lactation. The analysed data consisted of 23 paternal half-sib families (716 daughters of 23 sires) of Chinese Holstein-Friesian cattle genotyped at 14 microsatellites located in the area of the casein loci on BTA6. A sequence of three models was used: (i) a lactation model, (ii) a random regression model with a QTL constant in time and (iii) a random regression model with a QTL variable in time. The results showed that, for each production trait, at least one significant QTL exists. For milk and protein yields, the QTL effect was variable in time, while for fat yield, each of the three models resulted in a significant QTL effect. When a QTL is incorporated into a model as a constant over time, its effect is averaged over lactation stages and may, thereby, be difficult or even impossible to be detected. Our results showed that, in such a situation, only a longitudinal model is able to identify loci significantly influencing trait variation.     

A genome scan for quantitative trait loci affecting resistance to Trichostrongylus colubriformis in sheep

A genome linkage scan was carried out using a resource flock of 1029 sheep in six half-sib families. The families were offspring of sires derived by crossing divergent lines of sheep selected for response to challenge with the intestinal parasitic nematode Trichostrongylus colubriformis. All animals in the resource flock were phenotypically assessed for worm resistance soon after weaning using a vaccination/challenge regime. After correcting for fixed effects using a least squares linear model the faecal egg count data obtained following the first challenge and the faecal egg count data obtained after the second challenge were designated Trait 1 and Trait 2, respectively. A total of 472 lambs drawn from the phenotypic extremes of the Trait 2 faecal egg count distribution were genotyped with a panel of 133 microsatellite markers covering all 26 sheep autosomes. Detection of quantitative trait loci (QTL) for each of the faecal egg count traits was determined using interval analysis with the Animap program with recombination rates between markers derived from an existing marker map. No chromosomal regions attained genome-wide significance for QTL influencing either of the traits. However, one region attained chromosome-wide significance and five other regions attained point-wise significance for the presence of QTL affecting parasite resistance.     

Mycobacterium avium subspecies paratuberculosis cultured from locally and commercially pasteurized cow's milk in the Czech Republic

Between November 2002 and April 2003, 244 bottles and cartons of commercially pasteurized cow's milk were obtained at random from retail outlets throughout the Czech Republic. During the same period, samples of raw milk and of milk that was subsequently subjected to a minimum of 71.7 degrees C for 15 s in a local pasteurization unit were also obtained from two dairy herds, designated herds A and B, with low and high levels, respectively, of subclinical Mycobacterium avium subsp. paratuberculosis infection, and from one herd, herd C, without infection. Infection in individual cows in each herd was tested by fecal culturing. Milk samples were brought to the Veterinary Research Institute in Brno, Czech Republic, processed, inoculated onto Herrold's egg yolk slants, and incubated for 32 weeks. Colonies were characterized by morphology, Ziehl-Neelsen staining, mycobactin J dependency, and IS900 PCR results. M. avium subsp. paratuberculosis was cultured from 4 of 244 units (1.6%) of commercially pasteurized retail milk. M. avium subsp. paratuberculosis was also cultured from 2 of 100 (2%) cartons of locally pasteurized milk derived from infected herds A and B and from 0 of 100 cartons of milk from uninfected herd C. Raw milk from 1 of 10 (10%) fecal culture-positive cows in herd A and from 13 of 66 (19.7%) fecal culture-positive cows in herd B was culture positive for M. avium subsp. paratuberculosis. These findings confirm that M. avium subsp. paratuberculosis is present in raw milk from subclinically infected dairy cows. The culture of M. avium subsp. paratuberculosis in the Czech Republic from retail milk that had been pasteurized locally or commercially to the required national and European Union standards is in agreement with similar research on milk destined for consumers in the United Kingdom and the United States and shows that humans are being exposed to this chronic enteric pathogen by this route.     

Mapping Quantitative Trait Loci Controlling Milk Production in Dairy Cattle by Exploiting Progeny Testing

We have exploited ``progeny testing' to map quantitative trait loci (QTL) underlying the genetic variation of milk production in a selected dairy cattle population. A total of 1,518 sires, with progeny tests based on the milking performances of > 150,000 daughters jointly, was genotyped for 159 autosomal microsatellites bracketing 1645 centimorgan or approximately two thirds of the bovine genome. Using a maximum likelihood multilocus linkage analysis accounting for variance heterogeneity of the phenotypes, we identified five chromosomes giving very strong evidence (LOD score >/= 3) for the presence of a QTL controlling milk production: chromosomes 1, 6, 9, 10 and 20. These findings demonstrate that loci with considerable effects on milk production are still segregating in highly selected populations and pave the way toward marker-assisted selection in dairy cattle breeding.     

Chromosomal regions underlying noncoagulation of milk in Finnish Ayrshire cows

About 10% of Finnish Ayrshire cows produce noncoagulating milk, i.e., milk that does not form a curd in a standard 30-min testing time and is thus a poor raw material for cheese dairies. This phenomenon is associated with peak and midlactation, but some cows produce noncoagulating milk persistently. A genomewide scan under a selective DNA pooling method was carried out to locate genomic regions associated with the noncoagulation of milk. On the basis of the hypothesis of the same historical mutation, we pooled the data across sires. Before testing pools for homogeneity, allele intensities were corrected for PCR artifacts, i.e., shadow bands and differential amplification. Results indicating association were verified using daughter design and selective genotyping within families. Data consisted of 18 sire families with 477 genotyped daughters in total, i.e., 12% of each tail of the milk coagulation ability. Data were analyzed using interval mapping under maximum-likelihood and nonparametric methods. BMS1126 on chromosome 2 and BMS1355 on chromosome 18 were associated with noncoagulation of milk across families on an experimentwise 0.1% significance level. By scanning gene databases, we found two potential candidate genes: LOC538897, a nonspecific serine/threonine kinase on chromosome 2, and SIAT4B, a sialyltransferase catalyzing the last step of glycosylation of κ-casein on chromosome 18. Further studies to determine the role of the candidates in the noncoagulation of milk are clearly needed.     

Genome screen for twinning rate QTL in four North American Holstein families

The objective of this study was to identify twinning rate quantitative trait loci (QTL) by typing pooled samples in a preliminary screening followed by interval mapping to test QTL effects. Four elite North American Holstein half-sib sire families with high twinning rate predicted transmitting abilities (PTA) were used in this study. Chromosomes 5, 7, 19 and 23 were not genotyped as these chromosomes were scanned for QTL in these families in a previous study. DNA was extracted from phenotypically extreme sons in each sire family. Two pools were prepared from sons of sires in each phenotypic tail, two each for high and low PTA levels for twinning rates. Each pool contained DNA from 4 to 15% of all sons of the sire depending on family. A total of 268 fluorescently labelled microsatellite markers were tested for heterozygosity in sires. About 135--170 informative markers per family were genotyped using pooled DNA samples. Based on the preliminary evidence for potential twinning rate QTL from pooled typing, interval mapping was performed subsequently on 12 chromosomal regions by family combinations. Evidence of QTL for twinning rate was found in one family on BTA 21 and 29 at a chromosome-wide P<0.05 and on BTA 8, 10 and 14 with a chromosome-wide P<0.01.     

Low Fertility in Daughters of Bulls with 1/29 Translocation

Non-return rates at different stages after insemination were compared in daughters (n = 21,212) of five translocation heterozygous sires with all cows and heifers (n = 610,714) represented in the Norwegian Red Cattle (NRF) breeding statistics during the period October 1970 to September 1972. The daughters of the carrier bulls showed statistically significant lower non-return rates for 0–30, 30–60, 60–90, 90–120, 120–150, 150–180 and 180–270 days as compared with the control animals. The decline in non-return percentage for the daughters of translocation sires was greater than in the control animals during the period from 0–30 to 60–90 days indicating higher losses of fetuses during the early period after insemination. Daughter groups of carrier bulls consisting of adult cows only showed higher non-return rates than groups comprising both cows and heifers.     

Impact of bovine somatotropin on ranking for genetic value of dairy sires for milk yield traits and somatic cell score

Records of Holstein cows were used to examine how different models account for the effect of bovine somatotropin (bST) treatment on genetic evaluation of dairy sires for yield traits and somatic cell score. Data set 1 included 65,720 first-lactation records. Set 2 included 50,644 second-lactation records. Set 3 included 45,505 records for lactations three, four and five. Estimated breeding values (EBV) of sires were with three different animal models. With Model 1, bST administration was ignored. With Model 2, bST administration was used as a fixed effect. With Model 3, administration of bST was used to define the contemporary group (herd-year-month of calving-bST). Correlations for EBV of 1,366 sires with treated daughters between pairs of the three models were calculated for milk, fat and protein yields and somatic cell score for the three data sets. Correlations for EBV of sires between pairs of models for all traits ranged from 0.971 to 0.999. Fractions of sires with bST-treated progeny selected in common (top 10 to 15%) were 0.94 and usually greater for all pairs of models for all traits and parities. For this study, the method of statistical adjustment for bST treatment resulted in a negligible effect on genetic evaluations of sires when some daughters were treated with bST and suggests that selection of sires to produce the next generation of sires and cows might not be significantly affected by how the effect of bST is modeled for prediction of breeding values for milk, fat and protein yields and somatic cell score.     

Detection of Mycobacterium avium subspecies paratuberculosis in Swiss dairy cattle by real-time PCR and culture: a comparison of the two assays

AIMS: To compare the two different diagnostic assays for the detection of Mycobacterium avium ssp. paratuberculosis, the aetiological agent of paratuberculosis. METHODS AND RESULTS: Faecal samples were derived from 310 cows, representing 13 commercial dairy herds in various locations in Switzerland with expected increased risk because of a past history of disease. Detection assays for M. avium ssp. paratuberculosis were culture (gold standard) and a newly designed real-time PCR. Real-time PCR identified 31 of 310 animals as positive within this risk population whereas culture identified 20 positive animals. The specificity of real-time PCR was confirmed by DNA sequencing of the PCR product. Depending on the test used, the paratuberculosis prevalence in our tested risk population ranged from 6.5 to 10%. CONCLUSIONS: Real-time PCR and culture data were in good agreement, and real-time PCR generates data in a short time in contrast to culture. SIGNIFICANCE AND IMPACT OF THE STUDY: We consider real-time PCR as a suitable alternative method to culture for the detection of M. avium ssp. paratuberculosis in a national surveillance programme.     

Studies on changes of estimated breeding values of U.S. Holstein bulls for final score from the first to second crop of daughters

The purpose of this study was to find ways of reducing changes of sire predicted transmitting ability for type’s final scores (PTATs) from the first to second crop of daughters. The PTATs were estimated from two datasets: D01 (scores recorded up to 2001) and D05 (scores recorded up to 2005). The PTAT changes were calculated as the difference between the evaluations based on D01 and D05. The PTATs were adjusted to a common genetic base of all evaluated cows born in 1995. The single-trait (ST) animal model included the fixed effects of the herd–year–season–classifier, age by year group at classification, stage of lactation at classification, registry status of animals, and additive genetic and permanent environment random effects. Unknown parent groups (UPGs) were defined based on every other birth year starting from 1972. Modifications to the ST model included the usage of a single record per cow, separate UPGs for first and second crop daughters, separate UPGs for sires and dams, and deepened pedigrees for dams with missing phenotypic records. Also, the multiple-trait (MT) model treated records of registered and grade cows as correlated traits. The mean PTAT change, for all of the sires, was close to zero in all of the models analyzed. The estimated mean PTAT change for 145 sires with 40 to 100 first crop and ≥200 second crop daughters was −0.33, −0.20, −0.13, −0.28, and −0.12 with ST, only first records, only last records, updated pedigrees, and allowing separate parent groups (PGs) for sires and dams after updating the pedigrees, respectively. The percentages of sires showing PTAT decline were reduced from 74.5 (with ST) to 57.3 by using only the last records of cows, and to 56.4 by allowing separate UPGs for sires and dams after updating the pedigrees. Though updating of the pedigrees alone was not effective, separate UPGs for sires together with additional pedigree was helpful in reducing the bias.     

A chromosome-wide QTL study on BTA29 affecting temperament traits in German Angus beef cattle and mapping of DRD4

The behaviour of beef cattle is important for the safety and welfare of stockmen and animals. Ten microsatellites spanning BTA29 and, in addition, the candidate gene, dopamine receptor D4 gene, were analysed in 545 German Angus calves of six sires and included in a quantitative trait locus (QTL) study on the basis of three different behaviour tests. A putative QTL for the score while entering the scale (ScE) was detected at BMS764. The DRD4 fragment was mapped in the distal region of BTA29 15.3 cM distal of ILSTS081. The results clearly indicate that BTA29 with a putative QTL in the proximal part and the candidate gene, DRD4, in the distal part plays an important role in the regulation of temperament. During the study one of the sires was detected to be a blood chimera.