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

Genetic parameters and trends for Johne's disease in US 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. The data were extracted from a JD data repository, maintained at the Council on Dairy Cattle Breeding, and initially supplied by 2 dairy record 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 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 EBVs 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. The correlations of bulls' PTA with economically important traits such as milk yield, milk protein, milk fat, somatic cell score, and mastitis were low, indicating a nonoverlapping 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.

The role of interleukin-10 receptor alpha (IL10Rα) in Mycobacterium avium subsp. paratuberculosis infection of a mammary epithelial cell line

BACKGROUND

Johne's disease is a chronic wasting disease caused by the bacterium Mycobacterium avium subspecies paratuberculosis (MAP). Johne's disease is highly contagious and MAP infection in dairy cattle can eventually lead to death. With no available treatment for Johne's disease, genetic selection and improvements in management practices could help reduce its prevalence. In a previous study, the gene coding interleukin-10 receptor subunit alpha (IL10Rα) was associated with Johne's disease in dairy cattle. Our objective was to determine how IL10Rα affects the pathogenesis of MAP by examining the effect of a live MAP challenge on a mammary epithelial cell line (MAC-T) that had IL10Rα knocked out using CRISPR/cas9. The wild type and the IL10Rα knockout MAC-T cell lines were exposed to live MAP bacteria for 72 h. Thereafter, mRNA was extracted from infected and uninfected cells. Differentially expressed genes were compared between the wild type and the IL10Rα knockout cell lines. Gene ontology was performed based on the differentially expressed genes to determine which biological pathways were involved.

RESULTS

Immune system processes pathways were targeted to determine the effect of IL10Rα on the response to MAP infection. There was a difference in immune response between the wild type and IL10Rα knockout MAC-T cell lines, and less difference in immune response between infected and not infected IL10Rα knockout MAC-T cells, indicating IL10Rα plays an important role in the progression of MAP infection. Additionally, these comparisons allowed us to identify other genes involved in inflammation-mediated chemokine and cytokine signalling, interleukin signalling and toll-like receptor pathways.

CONCLUSIONS

Identifying differentially expressed genes in wild type and ILR10α knockout MAC-T cells infected with live MAP bacteria provided further evidence that IL10Rα contributes to mounting an immune response to MAP infection and allowed us to identify additional potential candidate genes involved in this process. We found there was a complex immune response during MAP infection that is controlled by many genes.

A scoping review on associations between paratuberculosis and productivity in cattle

Paratuberculosis (PTB), or Johne's disease, is a disease with worldwide distribution caused by Mycobacterium avium subsp. paratuberculosis (MAP) that leads to chronic enteritis, primarily in ruminants. Even subclinical infection significantly reduces the animals' performance, and consequences of the disease lead to high economic losses for the cattle industry. To estimate the economic burden of bovine PTB and to evaluate the benefits of a potential control program, accurate estimates of the production effects associated with the disease are required. Therefore, the aim of this scoping review was to provide a comprehensive overview of associations between MAP infection and production parameters in cattle. The studies were collected from three electronic databases. Of the total 1,605 identified studies, 1,432 did not meet the set criteria in the title and abstract screening and a further 106 were excluded during full-text review. Finally, data on 34 different production parameters were extracted from 67 publications. Results show that the magnitude of reported performance losses varies depending on several factors, such as the type of diagnostic test applied, disease status or number of lactations. Studies reported a reduction in milk yield, changes in milk quality (e.g., higher somatic cell count, lower amount of produced milk fat and protein), reduced fertility (e.g., prolonged calving interval and service period, higher abortion rate and calving difficulties), reduced weaning weight, slaughter weight and slaughter value, or a higher risk for mastitis. Results from the studies included in our review show a median decrease of milk yield per infected cow of -452 kg/lactation for raw and -405 kg/lactation for modeled data. Similarly, the amount of produced milk protein fell by a median of -14.41 kg/lactation for modeled data and the amount of produced milk fat by a median of -13.13 kg/lactation. The reviewed studies revealed a prolonged calving interval by around 30 days and a 1.5 to 3 times higher likeliness of culling per lactation in PTB positive animals. Results from this scoping review provide evidence-based inputs for the development of economic models aiming at the estimation of the costs and benefits associated with different disease control scenarios for PTB.

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 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.

Association of paratuberculosis sero-status with milk production and somatic cell counts across 5 lactations, using multilevel mixed models, in dairy cows

The aim of this work was to investigate associations between individual cow Mycobacterium avium ssp. paratuberculosis (MAP) seropositivity, 305-d corrected milk production, and somatic cell count during 5 lactations lifespan in Portuguese dairy herds using multilevel mixed models. We used MAP serum ELISA (Idexx MAP Ac, Idexx Laboratories Inc., Westbrook, ME) results (n = 23,960) from all the 20,221 adult cows present in 329 farms and corresponding 47,586 lactation records from the National Dairy Improvement Association. Cows and farms were classified as positive or negative. Multilevel mixed models were used to investigate the association of cow MAP status with variation in milk production and somatic cell count. Cow MAP status, farm status, and lactation number were considered as independent variables. A quadratic function of lactation number was used to mimic the effect of lactation order on milk production. The models considered 3 levels: measurement occasion (level 1) within cow (level 2) and cow within farm (level 3). Four final models were produced, including all herds and cows, to address the effect of farm status (models 1 and 2) or the effect of cow status (models 3 and 4) on the outcome variables. Our results show that MAP status affects milk production. Losses are detectable from third lactation onward. During the first 5 lactations, positive cows accumulated an average loss of 1,284.8 kg of milk when compared with the negative cows. We also observed that somatic cell counts were higher in positive cows and a positive interaction occurs between cow status and lactation number, suggesting a positive association between MAP infection and increased somatic cell counts. Our results are in line with previous studies, suggesting a possible positive relation between cow milk production and susceptibility to MAP infection.

Knowledge gaps that hamper prevention and control of Mycobacterium avium subspecies paratuberculosis infection

In the last decades, many regional and country-wide control programmes for Johne's disease (JD) were developed due to associated economic losses, or because of a possible association with Crohn's disease. These control programmes were often not successful, partly because management protocols were not followed, including the introduction of infected replacement cattle, because tests to identify infected animals were unreliable, and uptake by farmers was not high enough because of a perceived low return on investment. In the absence of a cure or effective commercial vaccines, control of JD is currently primarily based on herd management strategies to avoid infection of cattle and restrict within-farm and farm-to-farm transmission. Although JD control programmes have been implemented in most developed countries, lessons learned from JD prevention and control programmes are underreported. Also, JD control programmes are typically evaluated in a limited number of herds and the duration of the study is less than 5 year, making it difficult to adequately assess the efficacy of control programmes. In this manuscript, we identify the most important gaps in knowledge hampering JD prevention and control programmes, including vaccination and diagnostics. Secondly, we discuss directions that research should take to address those knowledge gaps.

The effect of paratuberculosis on milk yield--A systematic review and meta-analysis

Bovine paratuberculosis is a disease characterized by chronic granulomatous enteritis causing protein-losing enteropathy. Adverse effects on animal productivity are key drivers in the attempt to control paratuberculosis at the farm level. Economic models require an accurate estimation of the production effects associated with paratuberculosis. The aim of this study was to conduct a systematic review and meta-analysis to investigate the effect of paratuberculosis on milk production. A total of 20 effect estimates from 15 studies were included in the final meta-analysis. Substantial between-study heterogeneity was observed. Subgroup analysis by case definition and study design was carried out to investigate heterogeneity. The majority of between-study variation was attributed to studies that defined cases on serology. Calculation of a pooled effect estimate was only appropriate for studies that defined cases by organism detection. A reduction in milk yield, corrected for lactation number and herd of origin of 1.87 kg/d, equivalent to 5.9% of yield, was associated with fecal culture or PCR positivity in individual cows.

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.

Short communication: Heritability estimates for susceptibility to Mycobacterium avium subspecies paratuberculosis infection defined by ELISA and fecal culture test results in Jersey cattle

Paratuberculosis (Johne's disease), an enteric disorder in ruminants caused by Mycobacterium avium ssp. paratuberculosis, causes economic losses in excess of $200 million annually to the US dairy industry. Costly diagnostic testing, cumbersome control programs, incurability, and ineffective vaccination all make M. avium ssp. paratuberculosis susceptibility a good candidate for genetic studies and genetic selection a potentially useful adjunct to management-based control programs. No report has been published for heritability of susceptibility to M. avium ssp. paratuberculosis infection in Jersey cattle. The objective of this study was to estimate variance components and heritability for susceptibility to M. avium ssp. paratuberculosis infection in US Jersey cattle. Data consisted of complete serum ELISA and partial fecal culture results on a total of 2,861 Jersey cows from 23 commercial herds throughout the United States after editing. Four M. avium ssp. paratuberculosis susceptibility phenotypes were defined using (1) ELISA sample-to-positive ratios as a continuous trait, (2) ELISA results as a binary trait (positive=1, negative=0), (3) ELISA results as an ordered categorical trait, and (4) a combined test in which ELISA and fecal culture results were both taken into account in a binary analysis. Three statistical models, including linear, binary threshold, and ordered threshold sire models, were used to analyze the data. All analyses were executed using the restricted maximum likelihood method in ASReml 3 software. The heritability estimates were low to moderate and ranged from 0.08 (±0.03) to 0.27 (±0.11) based on different trait definitions. The nonzero heritability indicates that susceptibility to M. avium ssp. paratuberculosis infection in Jersey cattle is influenced by genetic factors. Therefore, selection of the least susceptible animals could decrease genetic predisposition to M. avium ssp. paratuberculosis infection in Jersey populations in future generations.

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).