Genome-wide association analysis for susceptibility to infection by Mycobacterium avium ssp. paratuberculosis in US Holsteins

Paratuberculosis, or Johne's disease, is a chronic, granulomatous, gastrointestinal tract disease of cattle and other ruminants caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Control of Johne's disease is based on programs of testing and culling animals positive for infection with MAP and concurrently modifying management to reduce the likelihood of infection. The current study was motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. Two separate GWAS analyses were conducted, the first using 897 genotyped Holstein artificial insemination sires with phenotypes derived from incidence of MAP infection among daughters based on milk ELISA testing records. The second GWAS analysis was a case-control design using US Holstein cows phenotyped for MAP infection by serum ELISA or fecal culture tests. Cases included cows positive for either serum ELISA, fecal culture, or both. Controls consisted of animals negative for all tests conducted. A total of 376 samples (70 cases and 306 controls) from a University of Minnesota Johne's management demonstration project and 184 samples (76 cases and 108 controls) from a Michigan State University study were used. Medium-density (sires) and high-density (cows) genotype data were imputed to full genome sequence for the analyses. Marker-trait associations were analyzed using the single-step (ss)GWAS procedure implemented in the BLUPF90 suite of programs. Evidence of significant genomic contributions for susceptibility to MAP infection were observed on multiple chromosomes. Results were combined across studies in a meta-analysis, and increased support for genomic regions on BTA7 and BTA21 were observed. Gene set enrichment analysis suggested pathways for antigen processing and presentation, antimicrobial peptides and natural killer cell–mediated cytotoxicity are relevant to variation in host susceptibility to MAP infection, among others. Genomic prediction was evaluated using a 5-fold cross-validation, and moderate correlations were observed between genomic breeding value predictions and daughter averages (～0.43 to 0.53) for MAP infection in testing data sets. These results suggest that genomic selection against susceptibility to MAP infection is feasible in Holstein cattle.     

Short communication: Uncovering quantitative trait loci associated with resistance to Mycobacterium avium ssp. paratuberculosis infection in Holstein cattle using a high-density single nucleotide polymorphism panel

Mycobacterium avium ssp. paratuberculosis (MAP) is the etiological agent of Johne's disease in cattle. Johne's disease is a disease of significant economic, animal welfare, and public health concern around the globe. Therefore, understanding the genetic architecture of resistance to MAP infection has great relevance to advance genetic selection methods to breed more resistant animals. The objectives of this study were to perform a genome-wide association study of previously analyzed 50K genotypes now imputed to a high-density single nucleotide polymorphism panel (777K), aiming to validate previously reported associations and potentially identify additional single nucleotide polymorphisms associated with antibody response to MAP infection. A principal component regression-based genome-wide association study revealed 15 putative quantitative trait loci (QTL) associated with the MAP infection phenotype (serum or milk ELISA tests) on 9 different chromosomes (Bos taurus autosomes 5, 6, 7, 10, 14, 15, 16, 20, and 21). These results validated previous findings and identified new QTL on Bos taurus autosomes 15, 16, 20, and 21. The positional candidate genes NLRP3, IFi47, TRIM41, TNFRSF18, and TNFRSF4 lying within these QTL were identified. Further functional validation of these genes is now warranted to investigate their roles in regulating the immune response and, consequently, cattle resistance to MAP infection.     

Short communication: Heritability of susceptibility to infection by Mycobacterium avium ssp. paratuberculosis in Holstein cattle

Johne's disease in cattle is the result of infection of the small intestine by Mycobacterium avium ssp. paratuberculosis (MAP), leading to an incurable inflammatory bowel disease (Johne's disease or paratuberculosis). The disease is a concern both for its direct cost to dairy producers and for its zoonotic potential. The objective of this study was to estimate the heritability for susceptibility to infection of cattle by MAP using Johne's testing records (ELISA test for presence of antibodies to MAP in milk or blood) from US Holstein cattle from 2009 to 2016. Data sets were edited to include records from herds with 100 or more total records and sires with 50 or more daughters. Data sets were further edited to include (1) only herds with at least 1 positive test, (2) herds with at least 2.5% positive test results, and (3) herds with at least 5% positive test results to examine the effect of data from herds with higher proportions of positive tests, and presumably higher pathogen exposure, on heritability estimates. Two models were used in this study, a linear sire model and a binary threshold-probit sire model. Both were mixed models considering fixed effects of herd and age at test, the latter as a covariate accounting for linear and quadratic effects; random effects included sire and residual. Analyses were conducted using a restricted maximum likelihood method. Heritability estimates (±standard error) from the linear model were 0.041 ± 0.004, 0.050 ± 0.004, and 0.062 ± 0.007 for data from herds with at least 1 positive test, ≥2.5% positive tests, and ≥5% positive tests, respectively. Heritability estimates from the threshold model were 0.157 ± 0.014, 0.174 ± 0.016, and 0.186 ± 0.021 for data from herds with at least 1 positive test, ≥2.5% positive tests, and ≥5% positive tests, respectively. Heritability estimates from the linear model were affected by population incidence for positive tests, in contrast to estimates from the threshold model, likely accounting for the difference in magnitude of heritability estimates between models and suggesting that the threshold model analysis is the better choice. Heritability estimates increased as data were restricted to herds with presumed higher MAP exposure for both linear model and threshold model analyses. These estimates are similar to previous estimates in other dairy cattle populations and suggest the potential for selection to lessen susceptibility to MAP infection.     

Short communication: Heritability estimates for susceptibility to Mycobacterium avium ssp. paratuberculosis infection in Chinese Holstein cattle

Paratuberculosis in ruminants, which is caused by Mycobacterium avium ssp. paratuberculosis (MAP), is a contagious, chronic enteric disease associated with economic losses, animal welfare, and health implications in dairy cattle production. In this study, we estimated the variance components and heritability of susceptibility to MAP infection in Chinese Holstein cattle. We collected 4,937 serum samples from cows in 7 dairy herds in the Beijing region of China and used the ELISA test to detect antibodies to MAP. Three statistical models were implemented to estimate heritabilities: (1) a linear model (ELISA sample-to-positive ratios as a continuous trait); (2) a binary threshold model (positive/negative from ELISA results); and (3) an ordered threshold model (ELISA results as an ordered categorical model with categories 1 to 5 corresponding to negative, uncertain, mildly positive, intermediate positive, and strongly positive). The heritability estimates ranged from 0.0389 to 0.1069, indicating that genetic factors affect MAP infection susceptibility in Chinese Holstein cattle.     

Gene mapping,gene-set analysis,and genomic prediction of postpartum blood calcium in Holstein cows

The onset of lactation results in a sudden irreversible loss of Ca for colostrum and milk synthesis. Some cows are unable to quickly adapt to this demand and succumb to clinical hypocalcemia, whereas a larger proportion of cows develop subclinical hypocalcemia that predisposes them to other peripartum diseases. The objective of this study was to perform a comprehensive genomic analysis of blood total Ca concentration in periparturient Holstein cows. We first performed a genomic scan and a subsequent gene-set analysis to identify candidate genes, biological pathways, and molecular mechanisms affecting postpartum Ca concentration. Then, we assessed the prediction of postpartum Ca concentration using genomic information. Data consisted of 7,691 records of plasma or serum concentrations of Ca measured in the first, second, and third day after parturition of 959 primiparous and 1,615 multiparous cows that calved between December 2015 and June 2020 in 2 dairy herds. All cows were genotyped with 80k SNPs. The statistical model included lactation (1 to 5+), calf category (male, females, twins), and day as fixed effects, and season-treatment-experiment, animal, and permanent environmental as random effects. Model predictive ability was evaluated using 10-fold cross-validation. Heritability and repeatability estimates were 0.083 (standard error = 0.017) and 0.444 (standard error = 0.028). The association mapping identified 2 major regions located on Bos taurus autosome (BTA)6 and BTA16 that explained 1.2% and 0.7% of additive genetic variance of Ca concentration, respectively. Interestingly, the region on BTA6 harbors the GC gene, which encodes the vitamin D binding protein, and the region on BTA16 harbors LRRC38, which is actively involved in K transport. Other sizable peaks were identified on BTA5, BTA2, BTA7, BTA14, and BTA9. These regions harbor genes associated with Ca channels (CACNA1S, CRACR2A), K channels (KCNK9), bone remodeling (LRP6), and milk production (SOCS2). The gene-set analysis revealed terms related to vitamin transport, calcium ion transport, calcium ion binding, and calcium signaling. Genomic predictions of phenotypic and genomic estimated breeding values of Ca concentration yielded predictive correlations up to 0.50 and 0.15, respectively. Overall, the present study contributes to a better understanding of the genetic basis of postpartum blood Ca concentration in Holstein cows. In addition, the findings may contribute to the development of novel selection and management strategies for reducing periparturient hypocalcemia in dairy cattle.     

Longitudinal relationship between fecal culture,fecal quantitative PCR,and milk ELISA in Mycobacterium avium ssp. paratuberculosis-infected cows from low-prevalence dairy herds

Mycobacterium avium ssp. paratuberculosis (MAP), the causative agent of ruminant Johne's disease, presents a particular challenge with regard to infection mitigation on dairy farms. Diagnostic testing strategies to identify and quantify MAP and associated antibodies are imperfect, and certain facets of the relationship between diagnostic tests remain to be explored. Additional repeated-measures data from known infected animals are needed to complement the body of cross-sectional research on Johne's disease-testing methods. Statistical models that accurately account for multiple diagnostic results while adjusting for the effects of individual animals and herds over time can provide a more detailed understanding of the interplay between diagnostic outcomes. Further, test results may be considered as continuous wherever possible so as to avoid the information loss associated with dichotomization. To achieve a broader understanding of the relationship between diagnostic tests, we collected a large number of repeated fecal and milk samples from 14 infected cows, in addition to bulk milk samples, from 2 low-prevalence dairy herds in the northeast United States. Predominately through the use of mixed linear modeling, we identified strong associations between milk ELISA optical density, fecal quantitative PCR, and fecal culture in individual animals while concurrently adjusting for variables that could alter these relationships. Notably, we uncovered subtleties in the predictive abilities of fecal shedding level on milk ELISA results, with animals categorized as disease progressors reaching higher ELISA optical density levels. Moreover, we observed that spikes in fecal shedding could predict subsequent high ELISA values up to 2 mo later. We also investigated the presence of MAP in individual milk samples via PCR and noted an association between poor udder hygiene and MAP positivity in milk, suggesting some level of environmental contamination. The paucity of positive milk samples and the complete absence of detectable MAP in the bulk tank throughout the study period indicate that contamination of milk with MAP may not be of chief concern in low-prevalence herds. An enhanced understanding of the interrelationships between diagnostic tests can only benefit the development of testing strategies and objectives, which in turn may lessen MAP infection prevalence in dairy herds.     

Evaluation of environmental sampling and culture to determine Mycobacterium avium subspecies paratuberculosis distribution and herd infection status on US dairy operations

The objectives of this study were to determine the distribution of Mycobacterium avium subspecies paratuberculosis (MAP) in the environment and assess the relationship between the culture status of MAP in the farm environment and herd infection status. The National Animal Health Monitoring System's Dairy 2002 study surveyed dairy operations in 21 states. One component of the study involved collection and culturing of environmental samples for MAP from areas on farms where manure accumulated from a majority of a herd's cows. Operations were selected for inclusion based on perceived risk factors for MAP infection identified in a previously administered questionnaire. Individual animal and environmental samples were collected and used to determine the efficiency of environmental sampling for determination of herd infection status. Individual animal fecal, serum, and milk samples were used to classify herds as infected or not infected based on the presence of at least one test-positive animal in the herd. A total of 483 environmental samples (approximately 5 per farm) were collected, and 218 (45.1%) were culture-positive for MAP. A similar percentage of environmental cultures collected from all designated areas were positive [parlor exits (52.3%), floors of holding pens (49.1%), common alleyways (48.8%), lagoons (47.4%), manure spreaders (42.3%), and manure pits (41.5%)]. Of the 98 operations tested with the environmental sample culture, 97 had individual serum ELISA results, 60 had individual fecal culture results, and 34 had individual milk ELISA results. Sixty-nine of the 98 operations (70.4%) had at least one environmental sample that was culture-positive. Of the 50 herds classified as infected by fecal culture, 38 (76.0%) were identified by environmental culture. Two of the 10 operations classified as not infected based on individual animal fecal culture were environmental culture-positive. Of the 80 operations classified as infected based on serum ELISA-positive results, 61 (76.3%) were identified as environmental-positive, whereas 20 of the 28 (71.4%) operations identified as infected based on milk ELISA were detected by environmental sampling. Environmental sample culturing is less costly than individual animal sampling, does not require animal restraint, and identified more than 70% of infected operations. Environmental sampling is another diagnostic tool that veterinarians and dairy producers can use to determine herd infection status for MAP.     

Genome-wide association study for milk production traits in a Brazilian Holstein population

Advances in the molecular area of selection have expanded knowledge of the genetic architecture of complex traits through genome-wide association studies (GWAS). Several GWAS have been performed so far, but confirming these results is not always possible due to several factors, including environmental conditions. Thus, our objective was to identify genomic regions associated with traditional milk production traits, including milk yield, somatic cell score, fat, protein and lactose percentages, and fatty acid composition in a Holstein cattle population producing under tropical conditions. For this, 75,228 phenotypic records from 5,981 cows and genotypic data of 56,256 SNP from 1,067 cows were used in a weighted single-step GWAS. A total of 46 windows of 10 SNP explaining more than 1% of the genetic variance across 10 Bos taurus autosomes (BTA) harbored well-known and novel genes. The MGST1 (BTA5), ABCG2 (BTA6), DGAT1 (BTA14), and PAEP (BTA11) genes were confirmed within some of the regions identified in our study. Potential novel genes involved in tissue damage and repair of the mammary gland (COL18A1), immune response (LTTC19), glucose homeostasis (SLC37A1), synthesis of unsaturated fatty acids (LTBP1), and sugar transport (SLC37A1 and MFSD4A) were found for milk yield, somatic cell score, fat percentage, and fatty acid composition. Our findings may assist genomic selection by using these regions to design a customized SNP array to improve milk production traits on farms with similar environmental conditions.     

A scoping review of the testing of bulk tank milk to detect nonbacterial pathogens or herd exposure to nonbacterial pathogens in dairy cattle

In this scoping review, we characterized the literature reporting on the testing of bulk milk samples to detect microorganisms other than bacteria that can cause diseases in dairy cattle, including viruses, helminths, algae, and protozoa. A search strategy was completed by screening databases, conference proceedings, animal health agency websites, disease surveillance program websites, and handbooks of cattle-related diagnostic tests for potentially relevant articles. Two reviewers independently screened articles in English, Portuguese, or Spanish; original studies reporting on the testing of farm-level, unprocessed bulk milk samples for presence of pathogens or specific antibodies against agents other than bacteria that can cause diseases in cows were retained. From all studies, we used spreadsheets to extract relevant information, including pathogen screened, test used, and country of origin of bulk milk samples. Additionally, for studies reporting sufficient data to estimate test characteristics, we extracted detailed information about herd eligibility, testing protocol, and herd-level infection definition. A total of 8,829 records were identified, from which 1,592 were retained and assessed for eligibility, and 306 were included. Bovine viral diarrhea virus, Fasciola hepatica, Ostertagia ostertagi, and bovine herpesvirus 1 were the most frequently screened agents, reported from 107, 45, 45, and 33 studies, respectively. Sensitivity of bulk milk ELISA to detect herds with animals infected by bovine herpesvirus 1 ranged from 2 to 100%, and was affected mostly by antigen selection, cut-off adopted, herd vaccination status, and seroprevalence of lactating cows. Bulk milk ELISA had very high specificity to detect herds free of bovine leukemia virus, and varying sensitivity to detect herds with infected animals, which depended on the within-herd seroprevalence of lactating cattle. As for bovine viral diarrhea virus, in general, the sensitivity of bulk milk ELISA was moderate to high (>80%) when infection status was defined based on presence of persistently infected cattle or a high proportion of seropositive lactating cattle. Nevertheless, bulk milk ELISA was not able to distinguish infected and noninfected herds based on presence of seropositive unvaccinated weanlings. The PCR or quantitative PCR protocols employed had very low sensitivities (<40%) and very high specificities (>95%) to classify bovine viral diarrhea virus infection status of dairy herds. Sensitivity and specificity of bulk milk ELISA to classify herds with regards to presence of F. hepatica- or O. ostertagi-parasitized cattle were generally high and driven mostly by the definition of herd infection status. Conversely, bulk milk ELISA demonstrated varying characteristics to detect herds with or without Dictyocaulus viviparus-parasitized cattle, depending primarily on the antigen selected and presence of cattle with clinical signs of lungworm infection.     

Estimation of genetic parameters and single-step genome-wide association studies for milk urea nitrogen in Holstein cattle

The main objectives of this study were to estimate genetic parameters for milk urea nitrogen (MUN) in Holstein cattle and to conduct a single-step (ss)GWAS to identify candidate genes associated with MUN. Phenotypic measurements from 24,435 Holstein cows were collected from March 2013 to July 2019 in 9 dairy farms located in the Beijing area, China. A total of 2,029 cows were genotyped using the Illumina 150K Bovine Bead Chip, containing 121,188 SNP. A single-trait repeatability model was used to evaluate the genetic background of MUN. We found that MUN is a trait with low heritability (0.06 ± 0.004) and repeatability (0.12). Considering similar milk production levels, a lower MUN concentration indicates higher nitrogen digestibility. The genetic correlations between MUN and milk yield, net energy concentration, fat percentage, protein percentage, and lactose percentage were positive and ranged from 0.02 to 0.26. The genetic correlation between MUN and somatic cell score (SCS) was negative (?0.18), indicating that animals with higher MUN levels tend to have lower SCS. Both ssGWAS and pathway enrichment analyses were used to explore the genetic mechanisms underlying MUN. A total of 18 SNP (located on BTA11, BTA12, BTA14, BTA17, and BTA18) were found to be significantly associated with MUN. The genes CFAP77, CAMSAP1, CACNA1B, ADGRB1, FARP1, and INTU are considered to be candidate genes for MUN. These candidate genes are associated with important biological processes such as protein and lipid metabolism and binding to specific proteins. This set of candidate genes, metabolic pathways, and their functions provide a better understanding of the genomic architecture and physiological mechanisms underlying MUN in Holstein cattle.     

Time to the occurrence of a decline in milk production in cows with various paratuberculosis antibody profiles

Infection with Mycobacterium avium ssp. paratuberculosis (MAP) in dairy cattle often results in reduced milk production and premature culling. Some test-positive animals can live for years without being affected by infection, whereas others are test negative when they die from the infection. Our objective was to describe the deviation in milk production of cows with various MAP antibody profiles compared with their repeatedly test-negative herdmates in the same parity. Data were obtained from herds participating in the Danish control program on paratuberculosis, for which 4 annual MAP antibody ELISA of individual cows were performed per herd per year. A total of 136,489 ELISA results from 38,998 dairy cows in 64 herds were used along with 484,285 test-day records on energy-corrected milk (ECM) yield. Cows were divided into 6 antibody groups based on their repeated milk ELISA results: A0) repeated ELISA negative; A1) ELISA negative, but only once; A2) ELISA negative on the last 3 tests, but with 1 previous positive result; A3) ELISA negative on the last test, but with 1 or more previous positive results; A4) last sample was ELISA positive, but all previous were negative; and A5) at least the last 2 samples were ELISA positive. The expected test-day kilograms of ECM by herd and parity were estimated for cows in antibody group A0. Deviations from expected milk production were then assessed for cows in the other antibody groups relative to the time of the first test-positive ELISA result (D 0). Cows in groups A2, A3, and A5 produced approximately 0.5 kg of ECM/d more than cows in group A0 at 300 d before D 0. Cows in group A4 had a decline in milk production from d 300 before D 0, with daily milk production reduced by 5 kg of ECM at 200 d after D 0. Milk production of cows in group A5 was reduced by 2.5 kg of ECM at 300 d after D 0 compared with 300 d before D 0, whereas cows in groups A2 and A3 produced 0.5 kg of ECM more than cows in group A0. The conclusions of the study were that 1) increasing the number ELISA tests increases the predictive value of ELISA for inference on milk production losses, 2) a combination of ELISA with assessment of observed milk production may be a valuable tool for decisions on culling, and 3) the declines in milk production attributable to MAP occurred over a long time period, and may not be realized by the herd manager without more advanced management tools such as the model proposed here.     

Evaluation of the Voluntary Johne's Disease Herd Status Program as a source of replacement cattle

The objectives of this study were to evaluate the perceived value that enrolled producers gained from participation in the Voluntary Johne's Disease Herd Status Program (VJDHSP), and to evaluate the risk of infection with Mycobacterium avium ssp. paratuberculosis (MAP) of cattle reared in a presumed Johne's free environment vs. cattle raised in an environment of unknown Johne's status. Producers enrolled in level 3 or 4 of the VJDHSP (98 and 99% probability, respectively, of being free from MAP infection) were interviewed via telephone. Producers were asked questions pertaining to their participation in the VJDHSP, and asked to identify herds to which they had sold replacement heifers. These cattle were presumed to be uninfected before sale. Fifty-nine cattle (identified as having been purchased into infected herds as heifers) were identified as having been raised in uninfected herds and sold to producers with herds of unknown Johne's disease status. On the purchasing farm, fecal and blood samples were collected from each cow of VJDHSP origin and 3 randomly selected home-reared cows per VJDHSP cow matched by lactation as controls. Samples were tested using commercial ELISA (serum) and bacterial culture (feces). Results indicated that enrolled producers saw value in VJDHSP participation, and that cattle reared in VJDHSP herds were less likely to be infected with MAP than herdmates as measured by serum ELISA MAP antibody and by fecal culture. This study provides evidence of the value of the VJDHSP in providing economic value to participants and supports the promotion of VJDHSP herds as a source of replacement cattle of low infection risk for MAP.     

Evaluation of national surveillance methods for detection of Irish dairy herds infected with Mycobacterium avium ssp. paratuberculosis

The aim of this study was to evaluate the utility and cost-effectiveness of a range of national surveillance methods for paratuberculosis in Irish dairy herds. We simulated alternative surveillance strategies applied to dairy cattle herds for the detection of Mycobacterium avium ssp. paratuberculosis (MAP)-infected herds (case-detection) or for estimation of confidence of herd freedom from infection (assurance testing). Strategies simulated included whole-herd milk or serum serology, serology on cull cows at slaughter, bulk milk tank serology, environmental testing, and pooled fecal testing. None of the strategies evaluated were ideal for widespread national case-detection surveillance. Herd testing with milk or serum ELISA or pooled fecal testing were the most effective methods currently available for detection of MAP-infected herds, with median herd sensitivity >60% and 100% herd specificity, although they are relatively expensive for widespread use. Environmental sampling shows promise as an alternative, with median herd sensitivity of 69%, but is also expensive unless samples can be pooled and requires further validation under Irish conditions. Bulk tank milk testing is the lowest cost option and may be useful for detecting high-prevalence herds but had median herd sensitivity <10% and positive predictive value of 85%. Cull cow sampling strategies were also lower cost but had median herd sensitivity <40% and herd positive predictive values of <50%, resulting in an increased number of test-positive herds, each of which requires follow-up herd testing to clarify status. Possible false-positive herd testing results associated with prior tuberculosis testing also presented logistical issues for both cull cow and bulk milk testing. Whole-herd milk or serum ELISA testing are currently the preferred testing strategies to estimate confidence of herd freedom from MAP in dairy herds due to the good technical performance and moderate cost of these strategies for individual herd testing. Cull cow serology and bulk tank milk sampling provide only minimal assurance value, with confidence of herd freedom increasing only minimally above the prior estimate. Different testing strategies should be considered when deciding on cost-effective approaches for case-detection compared with those used for building confidence of herd freedom (assurance testing) as part of a national program.     

Associations between paratuberculosis ELISA results and test-day records of cows enrolled in the Irish Johne's Disease Control Program

The effect of the Mycobacterium avium ssp. paratuberculosis (MAP) ELISA status on test-day milk performance of cows from Irish herds enrolled in the pilot national voluntary Johne's disease control program during 2013 to 2015 was estimated. A data set comprising 92,854 cows and 592,623 complete test-day records distributed across 1,700 herds was used in this study. The resulting ELISA outcome (negative, inconclusive, and positive) of each cow within each year of the program was used to allocate the cow into different scenarios representing the MAP status. At MAPscenario₁, all cows testing ELISA nonnegative (i.e., inconclusive and positive) were assigned a MAP-positive status; at MAPscenario₂ only cows testing ELISA-positive were assigned a MAP-positive status; at MAPscenario₃ only cows testing ELISA nonnegative (inconclusive or positive) and gathered exclusively from herds where at least 2 further ELISA nonnegative (inconclusive or positive) cows were found were assigned a MAP-positive status; at MAPscenario₄ only cows testing ELISA-positive that were gathered exclusively from herds where at least 2 further ELISA-positive cows were found were assigned a MAP-positive status. Milk outputs based on test-day records were standardized for fat and protein contents (SMY) and the effect of MAP ELISA status on the SMY was estimated by a linear mixed effects model structure. The SMY mean difference recorded at test day between cows with a MAP-positive status and those with a MAP-negative status within MAPscenario₁ was estimated at ?0.182 kg/test day; the mean difference was ?0.297 kg/test day for MAPscenario₂; for MAPscenario₃ mean difference between MAP-positive status and MAP test-negative cows was ?0.209 kg/test day, and for MAPscenario₄, the difference was ?0.326 kg/test day.     

Single-step genome-wide association for selected milk fatty acids in Dual-Purpose Belgian Blue cows

The aim of this study was to estimate genetic parameters and identify genomic regions associated with selected individual and groups of milk fatty acids (FA) predicted by milk mid-infrared spectrometry in Dual-Purpose Belgian Blue cows. The used data were 69,349 test-day records of milk yield, fat percentage, and protein percentage along with selected individual and groups FA of milk (g/dL milk) collected from 2007 to 2020 on 7,392 first-parity (40,903 test-day records), and 5,185 second-parity (28,446 test-day records) cows distributed in 104 herds in the Walloon Region of Belgium. Data of 28,466 SNPs, located on 29 Bos taurus autosomes (BTA), of 1,699 animals (639 males and 1,060 females) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic best linear unbiased prediction approach. The proportion of genetic variance explained by each 25-SNP sliding window (with an average size of ~2 Mb) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for candidate genes. Average daily heritability estimated for the included milk FA traits ranged from 0.01 (C4:0) to 0.48 (C12:0) and 0.01 (C4:0) to 0.42 (C12:0) in the first and second parities, respectively. Genetic correlations found between milk yield and the studied individual milk FA, except for C18:0, C18:1 trans, C18:1 cis-9, were positive. The results showed that fat percentage and protein percentage were positively genetically correlated with all studied individual milk FA. Genome-wide association analyses identified 11 genomic regions distributed over 8 chromosomes [BTA1, BTA4, BTA10, BTA14 (4 regions), BTA19, BTA22, BTA24, and BTA26] associated with the studied FA traits, though those found on BTA14 partly overlapped. The genomic regions identified differed between parities and lactation stages. Although these differences in genomic regions detected may be due to the power of quantitative trait locus detection, it also suggests that candidate genes underlie the phenotypic expression of the studied traits may vary between parities and lactation stages. These findings increase our understanding about the genetic background of milk FA and can be used for the future implementation of genomic evaluation to improve milk FA profile in Dual-Purpose Belgian Blue cows.     

Genetic analysis of leukosis incidence in United States Holstein and Jersey populations

Bovine leukosis (BL) is a retroviral disease caused by the bovine leukosis virus that affects only cattle. It is associated with decreased milk production and increased cull rates due to development of lymphosarcoma. The virus also affects the immune system. Infected cows display a weak response to some vaccinations. It is important to determine if the heritability of BL susceptibility is greater than zero, or if the environment is the only factor that can be used to reduce the transmission and incidence of the disease. Accordingly, the aim of this study was to estimate the heritability for BL incidence and the genetic merit of sires for leukosis resistance in Holstein and Jersey cattle. Continuous scores and binary milk ELISA results for 13,217 Holstein cows from 114 dairy herds across 16 states and 642 Jersey cows from 8 dairy herds were considered. Data were obtained from commercial testing records at Antel BioSystems (Lansing, MI). Out of the 13,859 animals tested, 38% were found to be infected with the disease. Linear and threshold animal models were used to analyze the continuous and binary data, respectively. Results from both models were similar in terms of estimated breeding values and variance components in their respective scales. Estimates of heritability obtained with the 2 approaches were approximately 8% for both breeds, indicating a considerable genetic component underlying BL disease incidence. The correlation between the estimated breeding values from the 2 models was larger than 0.90, and the lists of top 10% bulls selected from each model had about 80% overlap for both breeds. In summary, results indicate that a simple linear model using the continuous ELISA scores as the response variable was a reasonable approach for the genetic analysis of BL incidence in cattle. In addition, the levels of heritability found indicate that genetic selection could also be used to decrease susceptibility to bovine leukosis virus infection in Holstein and Jersey cattle. Further research is necessary to investigate the genetic correlations of BL with other production and reproduction traits, and to search for potential genomic regions harboring major genes affecting BL susceptibility.     

Genome-wide association study identifies loci influencing natural antibody titers in milk of Dutch Holstein-Friesian cattle

Natural antibodies (NAb) are produced without any antigenic stimulation as a part of the innate immune system and provide a first line of defense against pathogens. Hence, they may be a useful trait when estimating an animal's potential immune competence and in selection for disease resistance. The aim of this study was to identify genomic regions associated with different NAb traits in milk and potentially describe candidate genes. Milk samples from 1,695 first-lactation Holstein Friesian cows with titer measurements for keyhole limpet hemocyanin, lipopolysaccharide, lipoteichoic acid, and peptidoglycan-binding total NAb and isotypes IgG₁, IgM, and IgA were used. Genome-wide association study analyses were performed using imputed 777K SNP genotypes, accounting for relationships using pedigree information. Functional enrichment analysis was performed on the significantly associated genomic regions to look for candidate genes. For IgM NAb, significant associations (false discovery rate <0.05) were found on Bos taurus autosome (BTA) 17, 18, and 21 with candidate genes related to immunoglobulin structure and early B cell development. For IgG₁, associations were found on BTA3, and we confirmed a quantitative trait loci on BTA21 previously reported for IgG NAb in serum. Our results provide new insights into the regulation of milk NAb that will help unravel the complex relationship between milk immunoglobulins and disease resistance in dairy cattle.     

Risk factors for changing test classification in the Danish surveillance program for Salmonella in dairy herds

A surveillance program in which all cattle herds in Denmark are classified into Salmonella infection categories has been in place since 2002. Dairy herds were considered test negative and thus most likely free of infection if Salmonella antibody measurements were consistently low in bulk tank milk samples collected every 3 mo. Herds were considered test positive and thus most likely infected if the 4-quarter moving average bulk tank milk antibody concentration was high or if there was a large increase in the most recent measurement compared with the average value from the previous 3 samples. The objective of this study was to evaluate risk factors for changing from test negative to positive, which was indicative of herds becoming infected from one quarter of the year to the next, and risk factors for changing from test positive to negative, which was indicative of herds recovering from infection between 2 consecutive quarters of the year. The Salmonella serotypes in question were Salmonella Dublin or other serotypes that cross-react with the Salmonella Dublin antigen in the ELISA (e.g., some Salmonella Typhimurium types). Two logistic regression models that accounted for repeated measurements at the herd level and controlled for herd size and regional effects were used. Data from 2003 was used for the analyses. A change from test negative to positive occurred in 2.0% of the quarterly observations (n = 21,007) from test negative dairy herds. A change from test positive to negative occurred in 10.0% of quarterly observations (n = 6,168) available from test positive dairy herds. The higher the number of test-positive neighbor herds in the previous year-quarter, the more likely herds were to become test positive for Salmonella. The number of purchased cattle from test-positive herds was also associated with changing from test negative to positive. The bigger the herd, the more likely it was to change from negative to test positive. The effect of herd size on recovery was less clear. Large herds consisting mainly of large breeds or having test-positive neighbors in a 2-km radius were less likely to change from test positive to negative, whereas the breed and neighbor factors were not found to be important for small herds. Organic production was associated with remaining test positive, but not with becoming test positive. The results emphasize the importance of external and internal biosecurity measures to control Salmonella infections.     

The genetics of antibody response to paratuberculosis in dairy cattle

Genetic parameters were estimated for antibody response to paratuberculosis (Mycobacterium avium ssp. paratuberculosis) using milk ELISA test results, collected and analyzed by National Milk Records, from Holstein Friesian cows on UK dairy farms in their first 3 lactations. Milk ELISA test results were obtained from 2007 to 2012 and combined with milk recording data and pedigree information. The reduced data set edited for the purposes of genetic parameter estimation consisted of 148,054 milk ELISA records from 64,645 lactations in 40,142 cows of 908 sires, recorded in 641 herds. Milk ELISA test results were log_e-transformed and univariate analysis of 3 alternative animal models and equivalent sire models were considered. The most appropriate model included additive genetic and permanent environmental random effects, whereas maternal effects were significant according to likelihood ratio test and Akaike's information criterion but not for Bayesian information criterion. Heritability and repeatability estimates were 0.06 and 0.37, respectively, for the chosen animal model and its equivalent sire model. A subset of the data including herds with greater than 10% positive tests gave a slightly higher heritability of 0.08. Favorable but generally low significant genetic correlations were obtained between antibody response with 305-d milk yield (−0.16), 305-d protein yield (−0.16), log_e-transformed lactation-average somatic cell count (0.15), and the number of mastitis episodes (0.22). Thus, selection on the antibody response to paratuberculosis, should not be detrimental to production or udder health traits. Testing cattle for paratuberculosis is important for its use in control programs and although the heritability of antibody response was low, breeding against the disease might be a good prospect as a preventative measure to assist together with other approaches in an overall control strategy.     

Changes in the dynamics of Coxiella burnetii infection in dairy cattle: An approach to match field data with the epidemiological cycle of C. burnetii in endemic herds

This study aimed to evaluate changes in the epidemiological status of Coxiella burnetii in dairy cattle herds to better understand the epidemiology of the infection and to predict its evolution. Bulk-tank milk (BTM) and serum samples were collected from 94 dairy cattle herds and analyzed by ELISA (BTM and sera) and PCR (BTM) in study 1 (S₁). Two years later (study 2; S₂), the same farms were visited with a similar sampling approach. To estimate seroconversion during this period, blood samples were collected from the maximum possible number of animals surveyed in S₁. Environmental samples were collected in S₂ to identify active shedding. Farms were allocated into 3 different categories in each study according to PCR and ELISA results: category A, with BTM ELISA and PCR positive herds and at least 1 seropositive animal; category B, with BTM ELISA or PCR positive herds or individual sera positive; and category C, with all negative results among herds. Changes in herd category between S₁ and S₂ were grouped in 9 classes. Two statistical models, one to search for drives of within-herd changes in C. burnetii infection status and another to look for variables modulating individual changes in C. burnetii antibody level, were built. Several herds in category A in S₁ remained in that category 2 yr later, indicating that C. burnetii can remain within a herd for a long time. Most of the herds with seroconversion and detection of the bacterium in the environment belonged to category A, suggesting active and recent infections. Changes in the epidemiological status of herds were driven by local densities of domestic ruminants, showing the implication of neighbor reservoirs; whereas individual changes in antibody levels were modulated by variation in the epidemiological status of herds. Observed changes in epidemiological status allowed depiction of the hypothesized life cycle of C. burnetii within dairy cattle herds, which should be tested by future long-term series studies on C. burnetii infection to help fitting control measures (e.g., vaccination) to within-herd C. burnetii status.