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

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.

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.

Association between major histocompatibility complex haplotypes and susceptibility of unvaccinated and vaccinated cattle to paratuberculosis

Bovine Johne's disease (BJD) or paratuberculosis is caused by Mycobacterium avium spp. paratuberculosis (MAP) and is a worldwide problem among domestic and wild ruminants. While vaccines are available, natural differences in background immunity between breeds within species and between individuals within herds suggest that genetic differences may be able to be exploited in marker-assisted selection as an aid to disease control. The major histocompatibility complex (MHC) is an important component in immune recognition with considerable genetic variability. In this study, associations between the MHC and resistance to BJD were explored in dairy cattle across two herds in which some of the cattle had been vaccinated with Silirum® (n = 540 cows). A BJD susceptible animal was exposed to MAP and became infected, while a resistant animal was exposed but did not become infected. There are different ways to define both exposure and infection, with different levels of stringency, therefore many classifications of the same set of animals are possible and were included in the analysis. The polymorphic regions of major histocompatibility complex class I (MHC I) and class II (MHC II) genes were amplified from the genomic DNA by PCR and sequenced, targeting exons 2 and 3 of the classical and non-classical MHC I genes and exon 2 from the DRB3, DQA1, DQA2 + 3 and DQB MHC II genes. The frequencies of MHC I and MHC II haplotypes and alleles were determined in susceptible and resistant populations. In unvaccinated animals, seven MHC I haplotypes and seven MHC II haplotypes were associated with susceptibility while two MHC I and six MHC II haplotypes were associated with resistance (P < 0.05). In vaccinated animals, two MHC I and three MHC II haplotypes were associated with susceptibility, while one MHC I and two MHC II haplotypes were associated with resistance (P < 0.05). The alleles in significant haplotypes were also identified. Case definitions with higher stringency resulted in fewer animals being included in the analyses, but the power to detect an association was not reduced and there was an increase in strength and consistency of associations. Consistent use of stringent case definitions is likely to improve agreement in future association studies.

Identification of loci associated with pathological outcomes in Holstein cattle infected with Mycobacterium avium subsp. paratuberculosis using whole-genome sequence data

Bovine paratuberculosis (PTB), caused by Mycobacterium avium subsp. paratuberculosis (MAP), is a chronic granulomatous enteritis that affects cattle worldwide. According to their severity and extension, PTB-associated histological lesions have been classified into the following groups; focal, multifocal, and diffuse. It is unknown whether these lesions represent sequential stages or divergent outcomes. In the current study, the associations between host genetic and pathology were explored by genotyping 813 Spanish Holstein cows with no visible lesions (N = 373) and with focal (N = 371), multifocal (N = 33), and diffuse (N = 33) lesions in gut tissues and regional lymph nodes. DNA from peripheral blood samples of these animals was genotyped with the bovine EuroG MD Bead Chip, and the corresponding genotypes were imputed to whole-genome sequencing (WGS) data using the 1000 Bull genomes reference population. A genome-wide association study (GWAS) was performed using the WGS data and the presence or absence of each type of histological lesion in a case–control approach. A total of 192 and 92 single nucleotide polymorphisms (SNPs) defining 13 and 9 distinct quantitative trait loci (QTLs) were highly-associated (P ≤ 5 × 10⁻⁷) with the multifocal (heritability = 0.075) and the diffuse (heritability = 0.189) lesions, respectively. No overlap was seen in the SNPs controlling these distinct pathological outcomes. The identified QTLs overlapped with some QTLs previously associated with PTB susceptibility, bovine tuberculosis susceptibility, clinical mastitis, somatic cell score, bovine respiratory disease susceptibility, tick resistance, IgG level, and length of productive life. Pathway analysis with candidate genes overlapping the identified QTLs revealed a significant enrichment of the keratinization pathway and cholesterol metabolism in the animals with multifocal and diffuse lesions, respectively. To test whether the enrichment of SNP variants in candidate genes involved in the cholesterol metabolism was associated with the diffuse lesions; the levels of total cholesterol were measured in plasma samples of cattle with focal, multifocal, or diffuse lesions or with no visible lesions. Our results showed reduced levels of plasma cholesterol in cattle with diffuse lesions. Taken together, our findings suggested that the variation in MAP-associated pathological outcomes might be, in part, genetically determined and indicative of distinct host responses.

Johne's Disease in Dairy Cattle: An Immunogenetic Perspective

Johne's disease (JD), also known as paratuberculosis, is a severe production-limiting disease with significant economic and welfare implications for the global cattle industry. Caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP), JD manifests as chronic enteritis in infected cattle. In addition to the economic losses and animal welfare issues associated with JD, MAP has attracted public health concerns with potential association with Crohn's disease, a human inflammatory bowel disease. The lack of effective treatment options, such as a vaccine, has hampered JD control resulting in its increasing global prevalence. The disease was first reported in 1895, but in recognition of its growing economic impact, extensive recent research facilitated by a revolution in technological approaches has led to significantly enhanced understanding of the immunological, genetic, and pathogen factors influencing disease pathogenesis. This knowledge has been derived from a variety of diverse models to elucidate host-pathogen interactions including in vivo and in vitro experimental infection models, studies measuring immune parameters in naturally-infected animals, and by studies conducted at the population level to enable the estimation of genetic parameters, and the identification of genetic markers and quantitative trait loci (QTL) putatively associated with susceptibility or resistance to JD. The main objectives of this review are to summarize these recent developments from an immunogenetics perspective and attempt to extract the principal and common findings emerging from this wealth of recent information. Based on these analyses, and in light of emerging technologies such as gene-editing, we conclude by discussing potential future avenues for effectively mitigating JD in cattle.

Genome-wide association analysis identified both RNA-seq and DNA variants associated to paratuberculosis in Canadian Holstein cattle 'in vitro' experimentally infected macrophages

Background

Mycobacterium avium ssp. paratuberculosis (MAP) is the causative agent of paratuberculosis, or Johne’s disease (JD), an incurable bovine disease. The evidence for susceptibility to MAP disease points to multiple interacting factors, including the genetic predisposition to a dysregulation of the immune system. The endemic situation in cattle populations can be in part explained by a genetic susceptibility to MAP infection. In order to identify the best genetic improvement strategy that will lead to a significant reduction of JD in the population, we need to understand the link between genetic variability and the biological systems that MAP targets in its assault to dominate macrophages. MAP survives in macrophages where it disseminates. We used next-generation RNA (RNA-Seq) sequencing to study of the transcriptome in response to MAP infection of the macrophages from cows that have been naturally infected and identified as positive for JD (JD (+); n = 22) or negative for JD (healthy/resistant, JD (−); n = 28). In addition to identifying genetic variants from RNA-seq data, SNP variants were also identified using the Bovine SNP50 DNA chip.

Results

The complementary strategy allowed the identification of 1,356,248 genetic variants, including 814,168 RNA-seq and 591,220 DNA chip variants. Annotation using SnpEff predicted that the 2435 RNA-seq genetic variants would produce high functional effect on known genes in comparison to the 33 DNA chip variants. Significant variants from JD(+/−) macrophages were identified by genome-wide association study and revealed two quantitative traits loci: BTA4 and 11 at (P < 5 × 10^− 7). Using BovineMine, gene expression levels together with significant genomic variants revealed pathways that potentially influence JD susceptibility, notably the energy-dependent regulation of mTOR by LKB1-AMPK and the metabolism of lipids.

Conclusion

In the present study, we succeeded in identifying genetic variants in regulatory pathways of the macrophages that may affect the susceptibility of cows that are healthy/resistant to MAP infection. RNA-seq provides an unprecedented opportunity to investigate gene expression and to link the genetic variations to biological pathways that MAP normally manipulate during the process of killing macrophages. A strategy incorporating functional markers into genetic selection may have a considerable impact in improving resistance to an incurable disease. Integrating the findings of this research into the conventional genetic selection program may allow faster and more lasting improvement in resistance to bovine paratuberculosis in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07487-4.

Association of genetic polymorphisms related to Johne's disease with estimated breeding values of Holstein sires for milk ELISA test scores

BACKGROUND

Johne's disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls' estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle.

RESULTS

Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls' EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively.

CONCLUSIONS

The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs.

A review of paratuberculosis in dairy herds - Part 1: Epidemiology

Bovine paratuberculosis is a chronic infectious disease of cattle caused by Mycobacterium avium subspecies paratuberculosis (MAP). This is the first in a two-part review of the epidemiology and control of paratuberculosis in dairy herds. Paratuberculosis was originally described in 1895 and is now considered endemic among farmed cattle worldwide. MAP has been isolated from a wide range of non-ruminant wildlife as well as humans and non-human primates. In dairy herds, MAP is assumed to be introduced predominantly through the purchase of infected stock with additional factors modulating the risk of persistence or fade-out once an infected animal is introduced. Faecal shedding may vary widely between individuals and recent modelling work has shed some light on the role of super-shedding animals in the transmission of MAP within herds. Recent experimental work has revisited many of the assumptions around age susceptibility, faecal shedding in calves and calf-to-calf transmission. Further efforts to elucidate the relative contributions of different transmission routes to the dissemination of infection in endemic herds will aid in the prioritisation of efforts for control on farm.

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.

Case definition terminology for paratuberculosis (Johne's disease)

Paratuberculosis (Johne's disease) is an economically significant condition caused by Mycobacterium avium subsp. paratuberculosis. However, difficulties in diagnosis and classification of individual animals with the condition have hampered research and impeded efforts to halt its progressive spread in the global livestock industry. Descriptive terms applied to individual animals and herds such as exposed, infected, diseased, clinical, sub-clinical, infectious and resistant need to be defined so that they can be incorporated consistently into well-understood and reproducible case definitions. These allow for consistent classification of individuals in a population for the purposes of analysis based on accurate counts. The outputs might include the incidence of cases, frequency distributions of the number of cases by age class or more sophisticated analyses involving statistical comparisons of immune responses in vaccine development studies, or gene frequencies or expression data from cases and controls in genomic investigations. It is necessary to have agreed definitions in order to be able to make valid comparisons and meta-analyses of experiments conducted over time by a given researcher, in different laboratories, by different researchers, and in different countries. In this paper, terms are applied systematically in an hierarchical flow chart to enable classification of individual animals. We propose descriptive terms for different stages in the pathogenesis of paratuberculosis to enable their use in different types of studies and to enable an independent assessment of the extent to which accepted definitions for stages of disease have been applied consistently in any given study. This will assist in the general interpretation of data between studies, and will facilitate future meta-analyses.