Heritability estimates for antibody response to Mycobacterium avium subspecies paratuberculosis in German Holstein cattle

Progress in Paratuberculosis Control Programmes for Dairy Herds

While paratuberculosis control has been studied for over a century, knowledge gaps still exist regarding the uptake and efficacy of control programmes. This narrative review aims to summarise studies on control programmes presented at the IDF ParaTB Fora in 2021 and 2022 and the International Colloquium on Paratuberculosis in 2022. Studies were grouped by topic as follows: successful control, field studies, education and extension, voluntary and compulsory control programmes, and surveillance. Various Map control programmes resulted in a decreasing animal and herd level Map prevalence. Long-term stakeholder commitment, stable funding, involvement of herd veterinarians and incentives for farmers to participate were shown to be pivotal for long-term success. Control measures focused on vertical and calf-to-calf transmission may improve Map control in infected herds. Easy-to-capture visualisation of surveillance test results to inform participants on the progress of Map control in their herds was developed. The probability of freedom from disease and estimated within-herd prevalence were identified as good candidates for categorisation of herds to support low-risk trade of cattle. Results of the surveillance schemes may inform genetic selection for resistance to Map infection. In conclusion, successful paratuberculosis control is feasible at both the herd and country level provided that crucial prerequisites are met.

Genetic Parameters and Genetic Trends for Johne's Disease in U.S. Holsteins: An Updated Study

Johne's disease (JD) is an infectious enteric disease in ruminants, causing substantial economic loss annually worldwide. This work aimed to estimate JD's genetic parameters and the phenotypic and genetic trends by incorporating recent data. It also explores the feasibility of a national genetic evaluation for JD susceptibility in Holstein cattle in the United States (US). The data were extracted from a Johne's disease data repository, maintained at the Council on Dairy Cattle Breeding (CDCB), and initially supplied by 2 dairy records processing centers. The data comprised 365,980 Holstein cows from 1,048 herds participating in a voluntary control program for JD. Two protocol kits, Idexx Paratuberculosis Screening Ab Test (IDX) and Parachek® 2 (PCK), were used to analyze milk samples with the Enzyme-Linked Immunosorbent Assay (ELISA) technique. Test results from the first 5 parities were considered. An animal was considered infected if it had at least one positive outcome. The overall average of JD incidence was 4.72% in these US Holstein cattle. Genotypes of 78,964 SNP markers were used for 25,000 animals randomly selected from the phenotyped population. Variance components and genetic parameters were estimated based on 3 models, namely, a pedigree-only threshold model (THR), a single-step threshold model (ssTHR), and a single-step linear model (ssLR). The posterior heritability estimates of JD susceptibility were low to moderate: 0.11 to 0.16 based on the 2 threshold models and 0.05 to 0.09 based on the linear model. The average reliability of estimated breeding values of JD susceptibility using single-step analysis for animals with or without phenotypes varied from 0.18 (THR) to 0.22 (ssLR) for IDX and from 0.14 (THR) to 0.18 (ssTHR and ssLR) for PCK. Despite no prior direct genetic selection against JD, the estimated genetic trends of JD susceptibility were negative and highly significant (P-value <0.01). The correlations of bulls' predicted transmitting abilities with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a non-overlapping genetic selection process with traits in current genetic evaluations. Our results suggest the feasibility of reducing the JD incidence rate by incorporating it into the national genetic evaluation programs.

Advances in understanding the genetic architecture of antibody response to paratuberculosis in sheep by heritability estimate and LDLA mapping analyses and investigation of candidate regions using sequence-based data

BACKGROUND

Paratuberculosis is a contagious and incurable disease that is caused by Mycobacterium avium subsp. paratuberculosis (MAP) with significant negative effects on animal welfare and farm profitability. Based on a large naturally infected flock over 12 years, we analyzed repeated enzyme-linked immunosorbent assay tests (ELISA), OvineSNP50 BeadChip genotypes and whole-genome sequences imputed from 56 influential animals. The main goals were to estimate the genetic parameters of proxy traits for resistance to MAP, identify genomic regions associated with the host's immune response against MAP and search for candidate genes and causative mutations through association and functional annotation analyses of polymorphisms identified by sequencing.

RESULTS

Two variables were derived from ELISA tests. The first, a binary variable, assessed the infection status of each animal over the entire productive life, while the second considered the level of antibody recorded over time. Very similar results were obtained for both variables. Heritability estimates of about 0.20 were found and a significant region capturing 18% and 13% of the genetic variance was detected on ovine chromosome 20 by linkage disequilibrium and linkage analysis on OvineSNP50 positions. Functional annotation and association analyses on the imputed sequence polymorphisms that were identified in this region were carried out. No significant variants showed a functional effect on the genes that mapped to this region, most of which belong to the major histocompatibility complex class II (MHC II). However, the conditional analysis led to the identification of two significant polymorphisms that can explain the genetic variance associated with the investigated genomic region.

CONCLUSIONS

Our results confirm the involvement of the host's genetics in susceptibility to MAP in sheep and suggest that selective breeding may be an option to limit the infection. The estimated heritability is moderate with a relevant portion being due to a highly significant region on ovine chromosome 20. The results of the combined use of sequence-based data and functional analyses suggest several genes belonging to the MHC II as the most likely candidates, although no mutations in their coding regions showed a significant association. Nevertheless, information from genotypes of two highly significant polymorphisms in the region can enhance the efficiency of selective breeding programs.

New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation

Background

Bovine paratuberculosis, or Johne’s disease (JD), is a contagious and incurable disease caused by Mycobacterium avium subsp. paratuberculosis (MAP). It has adverse effects on animal welfare and is very difficult to control, leading to serious economic consequences. An important line of defense to this disease is host genetic resistance to MAP, which, when it will be more fully understood, could be improved through selective breeding. Using a large dataset of Holstein cows (161,253 animals including 56,766 cows with ELISA serological phenotypes and 12,431 animals with genotypes), we applied a single-step single nucleotide polymorphism (SNP) best linear unbiased prediction approach to investigate the genetic determinism underlying resistance to this disease (heritability estimate and identification of relevant genomic regions) and estimated genetic trends, reliability, and relative risk factors associated with genomic predictions.

Results

Resistance to JD was moderately heritable (0.14) and 16 genomic regions were detected that accounted for at least 0.05% of the breeding values variance (GV) in resistance to JD, and were located on chromosomes 1, 3, 5, 6, 7, 19, 20, 21, 23, 25, and 27, with the highest percentage of variance explained by regions on chromosomes 23 (0.36% GV), 5 (0.22% GV), 1 (0.14% GV), and 3 (0.13% GV). When estimated for the whole chromosomes, the autosomes with the largest overall contributions were chromosomes 3 (5.3% GV), 10 (4.8%), 23 (4.7%), 1 (3.6%), 7 (3.4%), 5 (2.9%), 12 (2.5%), 11 (2.2%), and 13 (2%). We estimated a slightly favorable genetic trend in resistance to JD over the last two decades, which can be explained by a low positive genetic correlation between resistance to JD and total merit index (+ 0.06). Finally, in a validation population of 907 cows, relatively reliable genomic predictions (reliability = 0.55) were obtained, which allowed the identification of cows at high risk of infection.

Conclusions

This study provides new insights into the genetic determinism of resistance to JD and shows that this trait can be predicted from SNP genotypes. It has led to the implementation of a single-step genomic evaluation that should rapidly become an effective tool for controlling paratuberculosis on French Holstein farms.

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.

Elucidating Transmission Patterns of Endemic Mycobacterium avium subsp. paratuberculosis Using Molecular Epidemiology

Mycobacterial diseases are persistent and characterized by lengthy latent periods. Thus, epidemiological models require careful delineation of transmission routes. Understanding transmission routes will improve the quality and success of control programs. We aimed to study the infection dynamics of Mycobacterium avium subsp. paratuberculosis (MAP), the causal agent of ruminant Johne’s disease, and to distinguish within-host mutation from individual transmission events in a longitudinally MAP-defined dairy herd in upstate New York. To this end, semi-annual fecal samples were obtained from a single dairy herd over the course of seven years, in addition to tissue samples from a selection of culled animals. All samples were cultured for MAP, and multi-locus short-sequence repeat (MLSSR) typing was used to determine MAP SSR types. We concluded from these precise MAP infection data that, when the tissue burden remains low, the majority of MAP infections are not detectable by routine fecal culture but will be identified when tissue culture is performed after slaughter. Additionally, we determined that in this herd vertical infection played only a minor role in MAP transmission. By means of extensive and precise longitudinal data from a single dairy herd, we have come to new insights regarding MAP co-infections and within-host evolution.

Candidate genes associated with the heritable humoral response to Mycobacterium avium ssp. paratuberculosis in dairy cows have factors in common with gastrointestinal diseases in humans

Infection of cattle with bovine paratuberculosis (i.e., Johne's disease) is caused by Mycobacterium avium ssp. paratuberculosis (MAP) and results in a chronic incurable gastroenteritis. This disease, which has economic ramifications for the cattle industry, is increasing in detected prevalence globally; subclinically infected animals can silently shed the bacterium into the environment for years, exposing contemporaries and hampering disease-control programs. The objective of the present study was to first quantify the genetic parameters for humoral response to MAP in dairy cattle. This was followed by a genome-based association analysis and subsequent downstream bioinformatic analyses from imputed whole genome sequence SNP data. After edits, ELISA test records were available on 136,767 cows; analyses were also undertaken on a subset of 33,818 of these animals from herds with at least 5 MAP ELISA-positive cows, with at least 1 of those positive cows being homebred. Variance components were estimated using univariate animal and sire linear mixed models. The heritability calculated from the animal model for humoral response to MAP using alternative phenotype definitions varied from 0.02 (standard error = 0.003) to 0.05 (standard error = 0.008). The genome-based associations were undertaken within a mixed model framework using weighted deregressed estimated breeding values as a dependent variable on 1,883 phenotyped animals that were ≥87.5% Holstein-Friesian. Putative susceptibility quantitative trait loci (QTL) were identified on Bos taurus autosome 1, 3, 5, 6, 8, 9, 10, 11, 13, 14, 18, 21, 23, 25, 26, 27, and 29; mapping the most significant SNP to genes within and overlapping these QTL revealed that the most significant associations were with the 10 functional candidate genes KALRN, ZBTB20, LPP, SLA2, FI3A1, LRCH3, DNAJC6, ZDHHC14, SNX1, and HAS2. Pathway analysis failed to reveal significantly enriched biological pathways, when both bovine-specific pathway data and human ortholog data were taken into account. The existence of genetic variation for MAP susceptibility in a large data set of dairy cows signifies the potential of breeding programs for reducing MAP susceptibility. Furthermore, the identification of susceptible QTL facilitates greater biological understanding of bovine paratuberculosis and potential therapeutic targets for future investigation. The novel molecular similarities identified between bovine paratuberculosis and human inflammatory bowel disease suggest potential for human therapeutic interventions to be translated to veterinary medicine and vice versa.

Genome-wide association study of Mycobacterium avium subspecies Paratuberculosis infection in Chinese Holstein

Background

Paratuberculosis is a contagious, chronic and enteric disease in ruminants, which is caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection, resulting in enormous economic losses worldwide. There is currently no effective cure for MAP infection or a vaccine, it is thus important to explore the genetic variants that contribute to host susceptibility to infection by MAP, which may provide a better understanding of the mechanisms of paratuberculosis and benefit animal genetic improvement. Herein we performed a genome-wide association study (GWAS) to identify genomic regions and candidate genes associated with susceptibility to MAP infection in dairy cattle.

Results

Using Illumina Bovine 50 K (54,609 SNPs) and GeneSeek HD (138,893 SNPs) chips, two analytical approaches were performed, GRAMMAR-GC and ROADTRIPS in 937 Chinese Holstein cows, among which individuals genotyped by the 50 K chip were imputed to HD SNPs with Beagle software. Consequently, 15 and 11 significant SNPs (P < 5 × 10^− 5) were identified with GRAMMAR-GC and ROADTDRIPS, respectively. A total of 10 functional genes were in proximity to (i.e., within 1 Mb) these SNPs, including IL4, IL5, IL13, IRF1, MyD88, PACSIN1, DEF6, TDP2, ZAP70 and CSF2. Functional enrichment analysis showed that these genes were involved in immune related pathways, such as interleukin, T cell receptor signaling pathways and inflammatory bowel disease (IBD), implying their potential associations with susceptibility to MAP infection. In addition, by examining the publicly available cattle QTLdb, a previous QTL for MAP was found to be overlapped with one of regions detected currently at 32.5 Mb on BTA23, where the TDP2 gene was anchored.

Conclusions

In conclusion, we identified 26 SNPs located on 15 chromosomes in the Chinese Holstein population using two GWAS strategies with high density SNPs. Integrated analysis of GWAS, biological functions and the reported QTL information helps to detect positional candidate genes and the identification of regions associated with susceptibility to MAP traits in dairy cattle.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5385-3) contains supplementary material, which is available to authorized users.

A new compartmental model of Mycobacterium avium subsp. paratuberculosis infection dynamics in cattle

Models of Mycobacterium avium subsp. paratuberculosis (MAP), a chronic infectious agent of cattle, are used to identify effective control programs. However, new biological findings show that adult infections occur and that infected animals can be separated into 2 paths: animals that will become high-shedding and, eventually, experience clinical disease (high-path); and animals that will shed only small quantities of MAP and will remain subclinical (low-path). Longitudinal data analysis found that high-path animals progress more quickly than previously believed. A standard model of MAP transmission in dairy herds was modified to include adult low-path infections and 2 infection pathways for infected calves. Analysis of this model showed that adult infection may play an important role in MAP dynamics on a dairy farm, and that the increased rate of progression for high-path animals influences both the prevalence and the persistence of MAP on a dairy farm. This new model will be able to determine the effectiveness of MAP control programs more accurately than previous models.

Susceptibility to and diagnosis of Mycobacterium avium subspecies paratuberculosis infection in dairy calves: A review

The primary objectives of paratuberculosis control programs are reducing exposure of calves to Mycobacterium avium subspecies paratuberculosis (MAP), reducing herd infection pressure and regular testing of cattle >36 months of age. Although control programs based on these principles have reduced prevalence of MAP infection in dairy herds, they have generally not eliminated the infection. Recent infection trial(s) have yielded new knowledge regarding diagnostic testing and age- and dose-dependent susceptibility to MAP infection. Calves up to 1 year of age are still susceptible to MAP infection; therefore, control programs should refrain from referring to specific ages with respect to susceptibility and prevention of new infections. Notwithstanding, lesions were more severe when calves were inoculated at 2 weeks versus 1 year of age. Furthermore, a high inoculation dose resulted in more pronounced lesions than a low inoculation dose, especially in young calves. Consequently, keeping infection pressure low should decrease the incidence of new MAP infections and severity of JD in cattle that do acquire the infection. It was also evident that early diagnosis of MAP infection was possible and could improve efficacy of control programs. Although its use will still need to be validated in the field, a combination of antibody ELISA and fecal culture in young stock, in addition to testing cattle >36 months of age when screening a herd for paratuberculosis, was expected to improve detection of dairy cattle infected with MAP. Although calves were inoculated using a standardized method in a controlled environment, there were substantial differences among calves with regards to immune response, shedding and pathology. Therefore, we inferred there were genetic differences in susceptibility. Important insights were derived from experimental infection trials. Therefore, it was expected that these could improve paratuberculosis control programs by reducing severity and incidence of JD by lowering infection pressure on-farm, and reducing exposure of young calves and older cattle. Furthermore, an earlier diagnosis could be achieved by combining ELISA and fecal shedding in young stock, in addition to testing cattle >36 months of age.

Phenotype definition is a main point in genome-wide association studies for bovine Mycobacterium avium ssp. paratuberculosis infection status

Paratuberculosis caused by Mycobacterium avium ssp. paratuberculosis (MAP) causes economic losses and is present in dairy herds worldwide. Different studies used different diagnostic tests to detect infection status and are the basis of genome-wide association (GWA) studies with inconsistent results. Therefore, the aim of this study was to identify and compare genomic regions associated with MAP susceptibility in the same cohort of cattle using different diagnostic tests. The GWA study was performed in German Holsteins within a case-control assay using 305 cows tested for MAP by fecal culture and additional with four different commercial ELISA-tests. Genotyping was performed with the Illumina Bovine SNP50 BeadChip. The results using fecal culture or ELISA test led to the identification of different genetic loci. Two single-nucleotide polymorphisms showed significant association with the ELISA-status. However, no significant association for MAP infection could be confirmed. Our results show that the definition of the MAP-phenotype has an important impact on the outcome of GWA studies for paratuberculosis.

Wingless-type MMTV integration site family member 2 (WNT2) gene is associated with resistance to MAP in faecal culture and antibody response in Holstein cattle

Mycobacterium avium subspecies paratuberculosis (MAP) is a pathogenic bacterium responsible for the lethal Johne's disease in cattle. So far, several genome-wide association studies (GWAS) have been carried out to identify chromosomal regions highly associated with Johne's disease. The aim of this study was to investigate the genetic variability within a pool of seven genes (LAMB1, DLD, WNT2, PRDM1, SOCS5, PTGER4 and IL10) indicated by former GWAS/RNA-Seq studies as putatively associated with MAP infections and to achieve a confirmation study of association with paratuberculosis susceptibility in a population of 324 German Holstein cattle (162 cases MAP positive and 162 controls MAP negative) using ELISA and fecal cultural tests. SNP validation and genotyping information are provided, quick methods for allelic discrimination were set up and transcription factor binding analyses were performed. The rs43390642:G>TSNP in the WNT2 promoter region is associated with paratuberculosis susceptibility (P = 0.013), suggesting a protective role of the T allele (P = 0.043; odds ratio 0.50 [0.25-0.97]). The linkage disequilibrium with the DLD rs134692583:A>T might suggest a combined mechanism of action of these neighboring genes in resistance to MAP infection, which is also supported by a significant effect shown by the haplotype DLD(T) /WNT2(T) (P = 0.047). In silico analysis predicted rs43390642:G>T and rs134692583:A>T as essential parts of binding sites for the transcription factors GR, C/EBPβ and GATA-1, hence suggesting a potential influence on WNT2 and DLD gene expression. This study confirmed the region on BTA 4 (UMD 3.1: 50639460-51397892) as involved in tolerance/resistance to Johne's disease. In addition, this study clarifies the involvement of the investigated genes in MAP infection and contributes to the understanding of genetic variability involved in Johne's disease susceptibility.

Genome-wide association study of susceptibility to infection by Mycobacterium avium subspecies paratuberculosis in Holstein cattle

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 while concurrently modifying management to reduce the likelihood of infection. The current study is motivated by the hypothesis that genetic variation in host susceptibility to MAP infection can be dissected and quantifiable associations with genetic markers identified. For this purpose, a case-control, genome-wide association study was conducted using US Holstein cattle phenotyped for MAP infection using a serum ELISA and/or fecal culture test. Cases included cows positive for either serum ELISA, fecal culture or both. Controls consisted of animals negative for the serum ELISA test or both serum ELISA and fecal culture when both were available. Controls were matched by herd and proximal birth date with cases. A total of 856 cows (451 cases and 405 controls) were used in initial discovery analyses, and an additional 263 cows (159 cases and 104 controls) from the same herds were used as a validation data set. Data were analyzed in a single marker analysis controlling for relatedness of individuals (GRAMMAR-GC) and also in a Bayesian analysis in which multiple marker effects were estimated simultaneously (GenSel). For the latter, effects of non-overlapping 1 Mb marker windows across the genome were estimated. Results from the two discovery analyses were generally concordant; however, discovery results were generally not well supported in analysis of the validation data set. A combined analysis of discovery and validation data sets provided strongest support for SNPs and 1 Mb windows on chromosomes 1, 2, 6, 7, 17 and 29.

Genome-wide association analysis and genomic prediction of Mycobacterium avium subspecies paratuberculosis infection in US Jersey cattle

Paratuberculosis (Johne’s disease), an enteric disorder in ruminants caused by Mycobacterium avium subspecies paratuberculosis (MAP), causes economic losses in excess of $200 million annually to the US dairy industry. To identify genomic regions underlying susceptibility to MAP infection in Jersey cattle, a case-control genome-wide association study (GWAS) was performed. Blood and fecal samples were collected from ∼5,000 mature cows in 30 commercial Jersey herds from across the US. Discovery data consisted of 450 cases and 439 controls genotyped with the Illumina BovineSNP50 BeadChip. Cases were animals with positive ELISA and fecal culture (FC) results. Controls were animals negative to both ELISA and FC tests that matched cases on birth date and herd. Validation data consisted of 180 animals including 90 cases (positive to FC) and 90 controls (negative to ELISA and FC), selected from discovery herds and genotyped by Illumina BovineLD BeadChip (∼7K SNPs). Two analytical approaches were used: single-marker GWAS using the GRAMMAR-GC method and Bayesian variable selection (Bayes C) using GenSel software. GRAMMAR-GC identified one SNP on BTA7 at 68 megabases (Mb) surpassing a significance threshold of 5×10⁻⁵. ARS-BFGL-NGS-11887 on BTA23 (27.7 Mb) accounted for the highest percentage of genetic variance (3.3%) in the Bayes C analysis. SNPs identified in common by GRAMMAR-GC and Bayes C in both discovery and combined data were mapped to BTA23 (27, 29 and 44 Mb), 3 (100, 101, 106 and 107 Mb) and 17 (57 Mb). Correspondence between results of GRAMMAR-GC and Bayes C was high (70–80% of most significant SNPs in common). These SNPs could potentially be associated with causal variants underlying susceptibility to MAP infection in Jersey cattle. Predictive performance of the model developed by Bayes C for prediction of infection status of animals in validation set was low (55% probability of correct ranking of paired case and control samples).

Heritability estimates for Mycobacterium avium subspecies paratuberculosis status of German Holstein cows tested by fecal culture

The objective of this study was to estimate genetic manifestation of Mycobacterium avium ssp. paratuberculosis (MAP) infection in German Holstein cows. Incorporated into this study were 11,285 German Holstein herd book cows classified as MAP-positive and MAP-negative animals using fecal culture results and originating from 15 farms in Thuringia, Germany involved in a paratuberculosis voluntary control program from 2008 to 2009. The frequency of MAP-positive animals per farm ranged from 2.7 to 67.6%. The fixed effects of farm and lactation number had a highly significant effect on MAP status. An increase in the frequency of positive animals from the first to the third lactation could be observed. Threshold animal and sire models with sire relationship were used as statistical models to estimate genetic parameters. Heritability estimates of fecal culture varied from 0.157 to 0.228. To analyze the effect of prevalence on genetic parameter estimates, the total data set was divided into 2 subsets of data into farms with prevalence rates below 10% and those above 10%. The data set with prevalence above 10% show higher heritability estimates in both models compared with the data set with prevalence below 10%. For all data sets, the sire model shows higher heritabilities than the equivalent animal model. This study demonstrates that genetic variation exists in dairy cattle for paratuberculosis infection susceptibility and furthermore, leads to the conclusion that MAP detection by fecal culture shows a higher genetic background than ELISA test results. In conclusion, fecal culture seems to be a better trait to control the disease, as well as an appropriate feature for further genomic analyses to detect MAP-associated chromosome regions.

Genetic susceptibility to paratuberculosis

Multiple studies indicate that host animal genetics play a role in susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection. However, due to differences in methods used to define MAP-infected animals and controls and differences in methods of genetic analysis, there is as yet no clear consensus on the genes or markers to reliably define the MAP infection susceptibility of any animal species. Meta-analysis of combined studies and larger studies will help resolve the situation in the coming years.

Immunological and pathological responses of red deer resistant or susceptible genotypes, to experimental challenge with Mycobacterium avium subsp. paratuberculosis

This study aimed to monitor the clinical, immunological and pathological changes in red deer for 49 weeks after experimental oral challenge with Mycobacterium avium subsp. paratuberculosis (MAP) and to assess the heritability of resistance in the offspring of two red stags. Eighteen young deer, which were bred from unselected hinds and sired by two stags resistant (R) or susceptible (S) to paratuberculosis, were challenged with MAP and monitored for 49 weeks. Biopsy samples of the jejunal lymph node were collected at Weeks 4 and 13 and at necropsy after euthanasia of clinically affected animals or when electively killed at Week 49. Three animals (two S and one R) developed clinical disease and were euthanised. The nine S offspring had significantly more severe lesions than the nine R offspring (Mantel-Haenszel Chi-square P=0.017). The average Lesion Severity Score (LSS) of R offspring was 5.9 (mild), and 7/9 had no or very mild lesions. In contrast, the LSS of S offspring averaged 11.7 (severe), and 7/9 had severe lesions. Most of the resistant, but not the susceptible, animals showed evidence of resolving lesions and a reduction in the number of MAP between 13 and 49 weeks after challenge. One R offspring appeared to completely cure itself, and progressed from mild culture-positive paratuberculosis lesions at Week 13 to having no signs of disease or infection 36 weeks later. This study showed significant heritable resistance/susceptibility to paratuberculosis and key differences in immunological responses in the first 3 months after challenge, indicating different paths to relative success or failure to control MAP. In general, R deer had higher IFN-γ levels, low antibody titres and fewer MAP, while S deer had lower IFN-γ levels, higher antibody and more MAP.

Genetics of animal health and disease in cattle

There have been considerable recent advancements in animal breeding and genetics relevant to disease control in cattle, which can now be utilised as part of an overall programme for improved cattle health. This review summarises the contribution of genetic makeup to differences in resistance to many diseases affecting cattle. Significant genetic variation in susceptibility to disease does exist among cattle suggesting that genetic selection for improved resistance to disease will be fruitful. Deficiencies in accurately recorded data on individual animal susceptibility to disease are, however, currently hindering the inclusion of health and disease resistance traits in national breeding goals. Developments in 'omics' technologies, such as genomic selection, may help overcome some of the limitations of traditional breeding programmes and will be especially beneficial in breeding for lowly heritable disease traits that only manifest themselves following exposure to pathogens or environmental stressors in adulthood. However, access to large databases of phenotypes on health and disease will still be necessary. This review clearly shows that genetics make a significant contribution to the overall health and resistance to disease in cattle. Therefore, breeding programmes for improved animal health and disease resistance should be seen as an integral part of any overall national disease control strategy.

Genetic loci involved in antibody response to Mycobacterium avium ssp. paratuberculosis in cattle

BACKGROUND

Mycobacterium avium subsp. paratuberculosis (MAP) causes chronic enteritis in a wide range of animal species. In cattle, MAP causes a chronic disease called Johne's disease, or paratuberculosis, that is not treatable and the efficacy of vaccine control is controversial. The clinical phase of the disease is characterised by diarrhoea, weight loss, drop in milk production and eventually death. Susceptibility to MAP infection is heritable with heritability estimates ranging from 0.06 to 0.10. There have been several studies over the last few years that have identified genetic loci putatively associated with MAP susceptibility, however, with the availability of genome-wide high density SNP maker panels it is now possible to carry out association studies that have higher precision.

METHODOLOGY/PRINCIPAL FINDINGS

The objective of the current study was to localize genes having an impact on Johne's disease susceptibility using the latest bovine genome information and a high density SNP panel (Illumina BovineSNP50 BeadChip) to perform a case/control, genome-wide association analysis. Samples from MAP case and negative controls were selected from field samples collected in 2007 and 2008 in the province of Lombardy, Italy. Cases were defined as animals serologically positive for MAP by ELISA. In total 966 samples were genotyped: 483 MAP ELISA positive and 483 ELISA negative. Samples were selected randomly among those collected from 119 farms which had at least one positive animal.

CONCLUSION/SIGNIFICANCE

THE ANALYSIS OF THE GENOTYPE DATA IDENTIFIED SEVERAL CHROMOSOMAL REGIONS ASSOCIATED WITH DISEASE STATUS: a region on chromosome 12 with high significance (P<5x10(-6)), while regions on chromosome 9, 11, and 12 had moderate significance (P<5x10(-5)). These results provide evidence for genetic loci involved in the humoral response to MAP. Knowledge of genetic variations related to susceptibility will facilitate the incorporation of this information into breeding programmes for the improvement of health status.

Evidence of genetic resistance of cattle to infection with Mycobacterium bovis

Anecdotal evidence points to genetic variation in resistance of cattle to infection with Mycobacterium bovis, the causative agent of bovine tuberculosis (BTB), and published experimental evidence in deer and cattle suggests significant genetic variation in resistance and reactivity to diagnostic tests. However, such genetic variation has not been properly quantified in the United Kingdom dairy cattle population; it is possible that it exists and may be a factor influencing the occurrence of BTB. Using models based on the outcome of the process of diagnosis (ultimate fate models) and on the outcome of a single stage of diagnosis (continuation ratio models, herd test-date models), this study shows that there is heritable variation in individual cow susceptibility to BTB, and that selection for milk yield is unlikely to have contributed to the current epidemic. Results demonstrate that genetics could play an important role in controlling BTB by reducing both the incidence and the severity of herd breakdowns.

Short communication: Genetic and nongenetic factors influencing Ostertagia ostertagi antibodies in UK Holstein-Friesian cattle

The purpose of this study was to estimate and discuss the genetic variation, heritability, and effects of nongenetic factors on the ability of Holstein-Friesian cows to produce an immune response by producing IgG antibodies to Ostertagia ostertagi. Total IgG (IgG(1) and IgG(2)) antibody levels were determined using an ELISA and measured using optical density ratio (ODR=OD(sample) - OD(negative control)/OD(positive control) - OD(negative control)) from milk samples collected from 1,276 Holstein-Friesian cattle in 229 commercial dairy farms from 2002 to 2004 during their first (82%) and other (2 to 12) lactations. A sire (n=461) model was fitted to the ODR data using ASREML software, and variance components were estimated. The ability to produce O. ostertagi antibodies as measured by ODR had a heritability of 0.13+/-0.12, and both season of sample and herd had a significant effect on total IgG levels. To conclude, this study has ascertained that genetic variation is present in the ability of dairy cows to mount an immune response to the parasite O. ostertagi. Inasmuch as evidence exists that IgG is linked to protective immunity against the parasite via a reduction in its reproductive ability, this trait may be of potential interest to genetic selection programs as an aid to reduce the effect of O. ostertagi in dairy herds.

Gene expression profiling of PBMCs from Holstein and Jersey cows sub-clinically infected with Mycobacterium avium ssp. paratuberculosis

Infection of calves with intracellular Mycobacterium avium ssp. paratuberculosis (MAP) commonly results in a granulomatous, chronic inflammatory bowel disease known as Johne's disease. The asymptomatic stage of this infection can persist for the entire production life of an adult cow, resulting in reduced performance and premature culling, as well as transmission of MAP to progeny and herd-mates. It has been previously shown that the gene expression profiles of peripheral blood mononuclear cells (PBMCs) of healthy cows, and those chronically infected with MAP are inherently different, and that these changes may be indicative of disease progression. Since resistance to MAP infection is a heritable trait, and has been proposed to differ amongst domestic dairy cattle breeds, the objective of the present study was to compare gene expression profiles of PBMCs from healthy adult Holstein and Jersey cows to those considered to be sub-clinically infected with MAP, as indicated by serum ELISA. Microarray analysis using a platform containing more than 10,000 probes and ontological analysis identified differences in gene expression between a) healthy and infected cows, including genes involved in the inflammatory response, and calcium binding, and b) infected Holsteins and Jerseys, including genes involved in the immune response, and antigen processing and presentation. These results suggest a mixed pro- and anti-inflammatory phenotype of PBMCs from MAP-infected as compared to healthy control animals, and inherently different levels of immune and inflammatory-related gene expression between MAP-infected Holsteins and Jerseys.

A whole genome association analysis identifies loci associated with Mycobacterium avium subsp. paratuberculosis infection status in US holstein cattle

The purpose of this study was to identify loci associated with Mycobacterium avium subspecies paratuberculosis (Map) infection status in US Holsteins using the Illumina BovineSNP50 BeadChip whole genome single nucleotide polymorphism (SNP) assay. Two hundred forty-five cows from dairies in New York, Pennsylvania and Vermont enrolled in longitudinal herd studies between January 1999 and November 2007 were assessed for the presence of Map in both faecal and tissue samples. An animal was considered tissue infected if any sample contained at least one colony forming unit of Map per gram of tissue (CFU/g) and the same definition was employed for faecal samples. Each animal was genotyped with the Illumina BovineSNP50 BeadChip and after quality assurance filtering, 218 animals and 45 683 SNPs remained. We sought to identify loci associated with four different case/control classifications: presence of Map in the tissue, presence of Map in faeces, presence of Map in both tissue and faeces and presence of Map in tissue but not faeces. A case-control genome wide association study was conducted to test the four different classifications of Map infection status (cases) when compared with a Map-negative control group (control). Regions on chromosomes 1, 5, 7, 8, 16, 21 and 23 were identified with moderate significance (P < 5 x 10(-5)). Two regions, one on chromosome 3 (near EDN2) and another on chromosome 9 (no positional gene candidates), were identified with a high level of association to the presence of Map in tissue and both tissue and faeces respectively (P < 5 x 10(-7), genome-wide Bonferonni P < 0.05).