Evidence of genetic resistance of cattle to infection with Mycobacterium bovis

Integrative genomics sheds light on the immunogenetics of tuberculosis in cattle

Mycobacterium bovis causes bovine tuberculosis (bTB), an infectious disease of cattle that represents a zoonotic threat to humans. Research has shown that the peripheral blood (PB) transcriptome is perturbed during bTB disease but the genomic architecture underpinning this transcriptional response remains poorly understood. Here, we analyse PB transcriptomics data from 63 control and 60 confirmed M. bovis-infected animals and detect 2592 differently expressed genes perturbing multiple immune response pathways. Leveraging imputed genome-wide SNP data, we characterise thousands of cis-expression quantitative trait loci (eQTLs) and show that the PB transcriptome is substantially impacted by intrapopulation genomic variation during M. bovis infection. Integrating our cis-eQTL data with bTB susceptibility GWAS summary statistics, we perform a transcriptome-wide association study and identify 115 functionally relevant genes (including RGS10, GBP4, TREML2, and RELT) and provide important new omics data for understanding the host response to mycobacterial infections that cause tuberculosis in mammals.

Detection of genetic variability in dairy cattle infectivity for bovine tuberculosis

This study investigated the genetics of bovine tuberculosis (bTB) infectivity in Holstein-Friesian dairy cows using British national data. The analyses included cows with recorded sires from herds affected by bTB outbreaks between 2000 and 2022. Animals were considered bTB-positive if they reacted positively to the skin test and/or had positive post-mortem findings. We introduced the "index case approach," based on the assumption that once the initial positively tested animals (index cases) are detected in a herd, subsequent infections (secondary cases) in the early stages of the breakdown are likely to be attributed to these animals. Genetic analysis of the number of secondary cases (NSC) associated with a given index case was used to establish evidence of genetic variability in bTB infectivity of cattle, and derive Estimated Breeding Values (EBVs) for infectivity for the sires of the index cases. Data were analyzed by employing Markov Chain Monte Carlo techniques to fit Generalized Linear Mixed Models with either Poisson, Zero-Inflated Poisson (ZIP), Hurdle Poisson, or Geometric distributions. All 4 models demonstrated presence of genetic variance in cattle infectivity, with the strongest evidence provided by the ZIP and Hurdle Poisson models. The Hurdle Poisson model offered the most accurate and least biased predictions. Sire infectivity EBVs from the Poisson, ZIP, and Geometric models showed strong concordance, with pairwise correlations of 0.90 or higher. In contrast, correlations between EBVs from the Hurdle Poisson model and the other models ranged from 0.36 to 0.39. The association of the sire infectivity EBVs with the average observed NSC per sire and the proportion of infectious index case daughters per sire was generally moderate with correlations between 44 and 47% and 65-69%, respectively. Agreement among models for identifying the genetically most infectious sires was also reasonable, with 151 out of 285 sires appearing in the top 10% across models, and 122 (42.8%) also aligning with the top 10% based on observed average NSC. Results provide novel evidence for exploitable genetic variance in bTB infectivity allowing the derivation of meaningful EBVs. Based on the estimated posterior mean genetic variances obtained, reduction in infectivity by one genetic standard deviation would result in a 32% - 44% decrease in the expected NSC per index case. Further research is warranted to refine the phenotypic definition of infectivity and assess correlation with other dairy traits.

Genome-Wide Association Study Reveals Quantitative Trait Loci and Candidate Genes Associated with High Interferon-gamma Production in Holstein Cattle Naturally Infected with Mycobacterium Bovis

Mycobacterium bovis (Mb) is the causative agent of bovine tuberculosis (bTb). Genetic selection aiming to identify less susceptible animals has been proposed as a complementary measure in ongoing programs toward controlling Mb infection. However, individual animal phenotypes for bTb based on interferon-gamma (IFNɣ) and its use in bovine selective breeding programs have not been explored. In the current study, IFNɣ production was measured using a specific IFNɣ ELISA kit in bovine purified protein derivative (bPPD)-stimulated blood samples collected from Holstein cattle. DNA isolated from the peripheral blood samples collected from the animals included in the study was genotyped with the EuroG Medium Density bead Chip, and the genotypes were imputed to whole-genome sequences. A genome-wide association analysis (GWAS) revealed that the IFNɣ in response to bPPD was associated with a specific genetic profile (heritability = 0.23) and allowed the identification of 163 SNPs, 72 quantitative trait loci (QTLs), 197 candidate genes, and 8 microRNAs (miRNAs) associated with this phenotype. No negative correlations between this phenotype and other phenotypes and traits included in the Spanish breeding program were observed. Taken together, our results define a heritable and distinct immunogenetic profile associated with strong production of IFNɣ in response to Mb.

Divergent proinflammatory immune responses associated with the differential susceptibility of cattle breeds to tuberculosis

Tuberculosis (TB) in the bovine is one of the most predominant chronic debilitating infectious diseases primarily caused by Mycobacterium bovis. Besides, the incidence of TB in humans due to M. bovis, and that in bovines (bovine TB, bTB) due to M. tuberculosis- indicates cattle as a major reservoir of zoonotic TB. While India accounts for the highest global burden of both TB and multidrug-resistant TB in humans, systematic evaluation of bTB prevalence in India is largely lacking. Recent reports emphasized markedly greater bTB prevalence in exotic and crossbred cattle compared to indigenous cattle breeds that represent more than one-third of the total cattle population in India, which is the largest globally. This study aimed at elucidating the immune responses underlying the differential bTB incidence in prominent indigenous (Sahiwal), and crossbred (Sahiwal x Holstein Friesian) cattle reared in India. Employing the standard Single Intradermal Tuberculin Test (SITT), and mycobacterial gene-targeting single as well as multiplex-PCR-based screening revealed higher incidences of bovine tuberculin reactors as well as Mycobacterium tuberculosis Complex specific PCR positivity amongst the crossbred cattle. Further, ex vivo mycobacterial infection in cultures of bovine peripheral blood mononuclear cells (PBMC) from SITT, and myco-PCR negative healthy cattle exhibited significantly higher intracellular growth of M. bovis BCG, and M. tuberculosis H37Ra in the crossbred cattle PBMCs compared to native cattle. In addition, native cattle PBMCs induced higher pro-inflammatory cytokines and signaling pathways, such as interferon-gamma (IFN-γ), interleukin-17 (IL-17), tank binding kinase-1 (TBK-1), and nitric oxide (NO) upon exposure to live mycobacterial infection in comparison to PBMCs from crossbred cattle that exhibited higher expression of IL-1β transcripts. Together, these findings highlight that differences in the innate immune responses of these cattle breeds might be contributing to the differential susceptibility to bTB infection, and the resultant disparity in bTB incidence amongst indigenous, and crossbred cattle.

Selective breeding can contribute to bovine tuberculosis control and eradication

Bovine tuberculosis (bTB) persists in many countries having a significant impact on public health and livestock industry finances. The incidence and prevalence of new cases in parts of the UK and elsewhere over the past decades warrant intensified efforts towards achieving Officially Tuberculosis Free (OTF) status in the respective regions. Genetic selection aiming to identify and remove inherently susceptible animals from breeding has been proposed as an additional measure in ongoing programmes towards controlling the disease. The presence of genetic variation among individual animals in their capacity to respond to Mycobacterium bovis exposure has been documented and heritability estimates of 0.06-0.18 have been reported. Despite their moderate magnitude, these estimates suggest that host resistance to bTB is amenable to improvement with selective breeding. Although relatively slow, genetic progress can be constant, cumulative and permanent, thereby complementing ongoing disease control measures. Importantly, mostly no antagonistic genetic correlations have been found between bTB resistance and other animal traits suggesting that carefully incorporating the former in breeding decisions should not adversely affect bovine productivity. Simulation studies have demonstrated the potential impact of genetic selection on reducing the probability of a breakdown to occur or the duration and severity of a breakdown that has already been declared. Furthermore, research on the bovine genome has identified multiple genomic markers and genes associated with bTB resistance. Nevertheless, the combined outcomes of these studies suggest that host resistance to bTB is a complex, polygenic trait, with no single gene alone explaining the inherent differences between resistant and susceptible animals. Such results support the development of accurate genomic breeding values that duly capture the collective effect of multiple genes to underpin selective breeding programmes. In addition to improving host resistance to bTB, scientists and practitioners have considered the possibility of reducing host infectivity. Ongoing studies have suggested the presence of genetic variation for infectivity and confirmed that bTB eradication would be accelerated if selective breeding considered both host resistance and infectivity traits. In conclusion, research activity on bTB genetics has generated knowledge and insights to support selective breeding as an additional measure towards controlling and eradicating the disease.

The regression of the bovine tuberculin reaction: Results from the Reactor Quality Assurance study in Northern Ireland

Tuberculin skin tests remain widely used in the control of bovine tuberculosis (bTB) in cattle. Little is known about the rate of regression of tuberculin reactions after the comparative intradermal cervical test (CICT) in cattle. This study aimed to collect data to describe tuberculin regression in reactors following the CICT at 72 ± 4 h post injection. Reactors were also tested using the interferon gamma (IFN-γ) assay to establish if any pattern existed between these results and the CICT reaction regression. The data were derived from 108 herds, 112 herd-level CICTs and 1008 animals. A multivariable linear mixed model was built to explore the regression of the bovine tuberculin reaction over time and the influence of potential predictors. The results confirmed a proportional decline in the bovine tuberculin reaction occurred over time. The predictors in the final model demonstrated that regression of the tuberculin reaction differed between reactors according to their IFN-γ test results and whether visible lesions were present at slaughter. Follow-up measurement of tuberculin reactions and the serial use of the IFN-γ assay in large breakdowns has the potential to provide both a mechanism for quality assurance of the current CICT bTB surveillance and the identification of atypical breakdowns or reactors requiring further investigation.

Risk factors and true prevalence of bovine tuberculosis in Bangladesh

Bovine tuberculosis (bTb) is endemic in Bangladesh but the true prevalence has not yet been reported. Our objectives for this study were to determine the true prevalence and identify risk factors for bTb at the animal- and herd-level in Bangladesh. A total of 510 cows were randomly selected during January 2018 to December 2018. Caudal fold (CFT) and comparative cervical tuberculin tests (CCT) were serially interpreted. Animal- and herd-level risk factor data were collected using a pre-tested questionnaire. The hierarchical true prevalence of bTb was estimated within a Bayesian framework. The herd- and animal-level risk factors were identified using mixed effects logistic regression. The apparent prevalence of bTb was 20.6% [95% Confidence Interval (CI): 17.3; 24.3] based on CFT. The animal-level true prevalence of bTb was 21.9 (13.0; 32.4). The herd-level true prevalence in different regions varied from 41.9% to 88.8%. The region-level true prevalence was 49.9 (13.8; 91.2). There is a 100% certainty that herds from Bhaluka and Mymensingh Sadar upazilas are not free from bTb. The odds of bTb were 3.9 times (1.2; 12.6) higher in herds having more than four cows than those with ≤ 4 cows. On the other hand, the risk of bTb was 3.3 times higher (1.0; 10.5) in non-grazing cows than grazing cows. Crossbred cows were 2.9 times (1.5; 5.9) more likely to be infected with bTb than indigenous cows. The risk of bTb in animals with cough was 2.3 times (1.2; 4.3) higher than those without cough. Crossbred, non-grazing cows with cough should be targeted for bTb surveillance. Herds of the Mymensingh, Sadar and Bhaluka regions should be emphasized for bTb control programs. Estimation of Bayesian hierarchical true prevalence facilitates identification of areas with higher prevalence and can be used to indicate regions that where true prevalence exceeds a pre-specified critical threshold.

Epidemiology of Bovine Tuberculosis and Its Zoonotic Implication in Addis Ababa Milkshed, Central Ethiopia

Bovine tuberculosis (bTB) continues to be one of the most widely distributed chronic infectious diseases of zoonotic importance, which causes a significant economic loss in animal production. A cross-sectional study was conducted to estimate the prevalence of bTB and its associated risk factors and type the Mycobacterium bovis isolated in central Ethiopia. A total of 65 dairy farms and 654 cattle were tested for bTB using a single intradermal comparative cervical tuberculin (SICCT) test. Data on farm management, animal-related characteristics, and the owner's knowledge of the zoonotic importance of bTB were collected using a structured questionnaire. In addition, a total of 16 animals from different farms were identified for postmortem examination. Lowenstein Jensen (LJ) culture was also conducted, and spoligotyping was used to type the M. bovis strains isolated. Chi-square test and logistic regression models were used to analyze the herd- and animal-level risk factors. Herd- and animal-level prevalence rates of bTB were 58.5% (95% CI: 46.2%-69.2%) and 39.3% (95% CI: 35.5%-43.5%), respectively. At the herd level, poor farm management was the predictor for bTB positivity (p < 0.05). Animal breed, poor BCS, farm type, and poor farm management conditions were significant predictors of bTB positivity (p < 0.05) at an individual animal level. All animals identified for postmortem examination were found to have gross TB-like lesions. A total of 14 M. bovis strains were identified from 12 animals that were positive for LJ culture. The strain with the largest number of clusters (five isolates) was SB1176, followed by SB0134 (three isolates), SB0192 (two isolates), and SB2233 (two isolates), and two new strains, each consisting of only one isolate. The majority (58.5%) of the respondents did not know the zoonotic importance of bTB. The result of this study showed a high prevalence of bTB in the Addis Ababa milkshed and a low level of consciousness of the owners on its transmission to humans. Therefore, the launching of acceptable control measures of bTB and the creation of public awareness about its zoonotic transmission and prevention measures are required.

Risk factors associated with testing positive for tuberculosis in high-yielding Holstein cows

Tuberculosis remains a significant infectious disease of farm animals and humans worldwide. The objective of this study was to assess various risk factors associated with testing positive for bovine tuberculosis (bTB) in high-yielding Holstein cows in an intensive dry-lot dairy operation. In a retrospective observational study, 9312 records from Holstein cows from a large dairy herd in northern Mexico were used. The incidence rate of lactating cows reactor to bTB was 7.3 cases/100 cow years (95% CI = 6.7-7.9%). Multiple logistic regression models indicated that cows with total milk yield during the first lactation >10,200 kg were 1.3 times (95% confidence intervals (CI) for odds ratio (OR) = 1.2-1.6) more likely to be detected as bTB reactors than cows with total milk yield <10,200 kg (8.2 vs. 6.4%; P < 0.01). The risk of being reactor to bTB increased with x̅ total milk yield/lactation >11,500 kg (OR=2.7, 95% CI = 2.3-3.3; 10.4 vs. 4.2%) compared with cows with x̅ total milk yield/lactation <11,500 kg. Cows with >10,200 kg 305-day milk yield in the last lactation was associated (P < 0.01) with 1.3 times higher odds (95% CI for OR = 1.1-1.5) of being classed as reactors than cows with <10,200 kg. Cows with peak milk yield >48 kg were 1.9 times (95% CI for OR = 1.6-2.2) more likely to be reactor to bTB than cows with peak milk yield <48 kg (9.2 vs. 5.1%; P < 0.01). Cows with either puerperal metritis (OR = 0.07, 95% CI = 0.5-0.9) or carrying twins (OR = 0.05, 95% CI = 0.01-0.19) had a protective role for being reactor to bTB. This study showed that increased milk production was associated with a higher risk of becoming positive to tuberculin skin test in high-yielding Holstein cows.

Genetic Diversity of Cameroon Cattle and a Putative Genomic Map for Resistance to Bovine Tuberculosis

Bovine Tuberculosis (bTB) caused by Mycobacterium bovis is a livestock disease of global economic and public health importance. There are currently no effective vaccines available for livestock and so control relies on animal level surveillance and pasteurization of dairy products. A new alternative control approach is to exploit the genetic variability of the host; recent studies have demonstrated that breeding European taurine cattle, such as Holsteins for increased resistance to bTB is feasible. The utility of such an approach is still unknown for African cattle populations. This study aims to assess genetic variation in bTB resistance and the underlying genomic architecture of cattle in Cameroon. We conducted a cross-sectional study of 2,346 slaughter cattle in Cameroon. Retropharyngeal lymph node samples were collected and cultured on Lowenstein Jensen media and the BACTEC MGIT 960 system, and M. bovis was identified using the Hain® Genotype kits. A total of 153 cattle were positive for M. bovis and were archived along with a random selection of negative samples. In this study, we genotyped archived samples from 212 cattle. Their genomic diversity was characterized using PCA, hierarchical clustering and admixture analysis. We assessed genetic variation in bTB resistance using heritability analysis and compared quantitative trait loci. Previous research on this study population have shown that Fulani cattle are more susceptible to bTB than mixed breeds. However, here we show that these apparent phenotypic differences in breeds are not reflected by clear genomic differences. At the genetic level, both the Fulani and mixed cattle show similar patterns of admixture with evidence of both taurine and indicine ancestry. There was little European taurine introgression within the studied population. Hierarchical clustering showed clusters of cattle that differed in their susceptibility to bTB. Our findings allude to bTB resistance being polygenic in nature. This study highlights the potential for genetic control of bTB in Africa and the need for further research into the genetics of bTB resistance within African cattle populations.

Predicting bovine tuberculosis status of dairy cows from mid-infrared spectral data of milk using deep learning

Bovine tuberculosis (bTB) is a zoonotic disease in cattle that is transmissible to humans, distributed worldwide, and considered endemic throughout much of England and Wales. Mid-infrared (MIR) analysis of milk is used routinely to predict fat and protein concentration, and is also a robust predictor of several other economically important traits including individual fatty acids and body energy. This study predicted bTB status of UK dairy cows using their MIR spectral profiles collected as part of routine milk recording. Bovine tuberculosis data were collected as part of the national bTB testing program for Scotland, England, and Wales; these data provided information from over 40,500 bTB herd breakdowns. Corresponding individual cow life-history data were also available and provided information on births, movements, and deaths of all cows in the study. Data relating to single intradermal comparative cervical tuberculin (SICCT) skin-test results, culture, slaughter status, and presence of lesions were combined to create a binary bTB phenotype labeled 0 to represent nonresponders (i.e., healthy cows) and 1 to represent responders (i.e., bTB-affected cows). Contemporaneous individual milk MIR spectral data were collected as part of monthly routine milk recording and matched to bTB status of individual animals on the single intradermal comparative cervical tuberculin test date (±15 d). Deep learning, a sub-branch of machine learning, was used to train artificial neural networks and develop a prediction pipeline for subsequent use in national herds as part of routine milk recording. Spectra were first converted to 53 × 20-pixel PNG images, then used to train a deep convolutional neural network. Deep convolutional neural networks resulted in a bTB prediction accuracy (i.e., the number of correct predictions divided by the total number of predictions) of 71% after training for 278 epochs. This was accompanied by both a low validation loss (0.71) and moderate sensitivity and specificity (0.79 and 0.65, respectively). To balance data in each class, additional training data were synthesized using the synthetic minority over sampling technique. Accuracy was further increased to 95% (after 295 epochs), with corresponding validation loss minimized (0.26), when synthesized data were included during training of the network. Sensitivity and specificity also saw a 1.22- and 1.45-fold increase to 0.96 and 0.94, respectively, when synthesized data were included during training. We believe this study to be the first of its kind to predict bTB status from milk MIR spectral data. We also believe it to be the first study to use milk MIR spectral data to predict a disease phenotype, and posit that the automated prediction of bTB status at routine milk recording could provide farmers with a robust tool that enables them to make early management decisions on potential reactor cows, and thus help slow the spread of bTB.

Genome-Wide Association Study in Mexican Holstein Cattle Reveals Novel Quantitative Trait Loci Regions and Confirms Mapped Loci for Resistance to Bovine Tuberculosis

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

Enhancing genetic disease control by selecting for lower host infectivity and susceptibility

Infectious diseases have a huge impact on animal health, production and welfare, and human health. Understanding the role of host genetics in disease spread is important for developing disease control strategies that efficiently reduce infection incidence and risk of epidemics. While heritable variation in disease susceptibility has been targeted in livestock breeding, emerging evidence suggests that there is additional genetic variation in host infectivity, but the potential benefits of including infectivity into selection schemes are currently unknown. A Susceptible-Infected-Recovered epidemiological model incorporating polygenic genetic variation in both susceptibility and infectivity was combined with quantitative genetics selection theory to assess the non-linear impact of genetic selection on field measures of epidemic risk and severity. Response to 20 generations of selection was calculated in large simulated populations, exploring schemes differing in accuracy and intensity. Assuming moderate genetic variation in both traits, 50% selection on susceptibility required seven generations to reduce the basic reproductive number R₀ from 7.64 to the critical threshold of <1, below which epidemics die out. Adding infectivity in the selection objective accelerated the decline towards R₀ < 1, to 3 generations. Our results show that although genetic selection on susceptibility reduces disease risk and prevalence, the additional gain from selection on infectivity accelerates disease eradication and reduces more efficiently the risk of new outbreaks, while it alleviates delays generated by unfavourable correlations. In conclusion, host infectivity was found to be an important trait to target in future genetic studies and breeding schemes, to help reducing the occurrence and impact of epidemics.

Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data

Bovine tuberculosis (bTB) is an infectious disease of cattle generally caused by Mycobacterium bovis, a bacterium that can elicit disease humans. Since the 1950s, the objective of the national bTB eradication program in Republic of Ireland was the biological extinction of bTB; that purpose has yet to be achieved. Objectives of the present study were to develop the statistical methodology and variance components to undertake routine genetic evaluations for resistance to bTB; also of interest was the detection of regions of the bovine genome putatively associated with bTB infection in dairy and beef breeds. The novelty of the present study, in terms of research on bTB infection, was the use of beef breeds in the genome-wide association and the utilization of imputed whole genome sequence data. Phenotypic bTB data on 781,270 animals together with imputed whole genome sequence data on 7,346 of these animals' sires were available. Linear mixed models were used to quantify variance components for bTB and EBVs were validated. Within-breed and multi-breed genome-wide associations were undertaken using a single-SNP regression approach. The estimated genetic standard deviation (0.09), heritability (0.12), and repeatability (0.30) substantiate that genetic selection help to eradicate bTB. The multi-breed genome-wide association analysis identified 38 SNPs and 64 QTL regions associated with bTB infection; two QTL regions (both on BTA23) identified in the multi-breed analysis overlapped with the within-breed analyses of Charolais, Limousin, and Holstein-Friesian. Results from the association analysis, coupled with previous studies, suggest bTB is controlled by an infinitely large number of loci, each having a small effect. The methodology and results from the present study will be used to develop national genetic evaluations for bTB in the Republic of Ireland. In addition, results can also be used to help uncover the biological architecture underlying resistance to bTB infection in cattle.

Can We Breed Cattle for Lower Bovine TB Infectivity?

Host resistance and infectivity are genetic traits affecting infectious disease transmission. This Perspective discusses the potential exploitation of genetic variation in cattle infectivity, in addition to resistance, to reduce the risk, and prevalence of bovine tuberculosis (bTB). In bTB, variability in M. bovis shedding has been previously reported in cattle and wildlife hosts (badgers and wild boars), but the observed differences were attributed to dose and route of infection, rather than host genetics. This article addresses the extent to which cattle infectivity may play a role in bTB transmission, and discusses the feasibility, and potential benefits from incorporating infectivity into breeding programmes. The underlying hypothesis is that bTB infectivity, like resistance, is partly controlled by genetics. Identifying and reducing the number of cattle with high genetic infectivity, could reduce further a major risk factor for herds exposed to bTB. We outline evidence in support of this hypothesis and describe methodologies for detecting and estimating genetic parameters for infectivity. Using genetic-epidemiological prediction models we discuss the potential benefits of selection for reduced infectivity and increased resistance in terms of practical field measures of epidemic risk and severity. Simulations predict that adding infectivity to the breeding programme could enhance and accelerate the reduction in breakdown risk compared to selection on resistance alone. Therefore, given the recent launch of genetic evaluations for bTB resistance and the UK government's goal to eradicate bTB, it is timely to consider the potential of integrating infectivity into breeding schemes.

Impact of Genetic Selection for Increased Cattle Resistance to Bovine Tuberculosis on Disease Transmission Dynamics

Bovine tuberculosis (bTB) poses a challenge to animal health and welfare worldwide. Presence of genetic variation in host resistance to Mycobacterium bovis infection makes the trait amenable to improvement with genetic selection. Genetic evaluations for resistance to infection in dairy cattle are currently available in the United Kingdom (UK), enabling genetic selection of more resistant animals. However, the extent to which genetic selection could contribute to bTB eradication is unknown. The objective of this study was to quantify the impact of genetic selection for bTB resistance on cattle-to-cattle disease transmission dynamics and prevalence by developing a stochastic genetic epidemiological model. The model was used to implement genetic selection in a simulated cattle population. The model considered various levels of selection intensity over 20 generations assuming genetic heterogeneity in host resistance to infection. Our model attempted to represent the dairy cattle population structure and current bTB control strategies in the UK, and was informed by genetic and epidemiological parameters inferred from data collected from UK bTB infected dairy herds. The risk of a bTB breakdown was modeled as the percentage of herds where initially infected cows (index cases) generated secondary cases by infecting herd-mates. The model predicted that this risk would be reduced by half after 4, 6, 9, and 15 generations for selection intensities corresponding to genetic selection of the 10, 25, 50, and 70% most resistant sires, respectively. In herds undergoing bTB breakdowns, genetic selection reduced the severity of breakdowns over generations by reducing both the percentage of secondary cases and the duration over which new secondary cases were detected. Selection of the 10, 25, 50, and 70% most resistant sires reduced the percentage of secondary cases to <1% in 4, 5, 7, and 11 generations, respectively. Similarly, the proportion of long breakdowns (breakdowns in which secondary cases were detected for more than 365 days) was reduced by half in 2, 2, 3, and 4 generations, respectively. Collectively, results suggest that genetic selection could be a viable tool that can complement existing management and surveillance methods to control and ultimately eradicate bTB.

An analysis of effects of heterozygosity in dairy cattle for bovine tuberculosis resistance

Genetic selection of cattle more resistant to bovine tuberculosis (bTB) may offer a complementary control strategy. Hypothesising underlying non-additive genetic variation, we present an approach using genome-wide high density markers to identify genomic loci with dominance effects on bTB resistance and to test previously published regions with heterozygote advantage in bTB. Our data comprised 1151 Holstein-Friesian cows from Northern Ireland, confirmed bTB cases and controls, genotyped with the 700K Illumina BeadChip. Genome-wide markers were tested for associations between heterozygosity and bTB status using marker-based relationships. Results were tested for robustness against genetic structure, and the genotypic frequencies of a significant locus were tested for departures from Hardy-Weinberg equilibrium. Genomic regions identified in our study and in previous publications were tested for dominance effects. Genotypic effects were estimated through ASReml mixed models. A SNP (rs43032684) on chromosome 6 was significant at the chromosome-wide level, explaining 1.7% of the phenotypic variance. In the controls, there were fewer heterozygotes for rs43032684 (P < 0.01) with the genotypic values suggesting that heterozygosity confers a heterozygote disadvantage. The region surrounding rs43032684 had a significant dominance effect (P < 0.01). SNP rs43032684 resides within a pseudogene with a parental gene involved in macrophage response to infection and within a copy-number-variation region previously associated with nematode resistance. No dominance effect was found for the region on chromosome 11, as indicated by a previous candidate region bTB study. These findings require further validation with large-scale data.

Variation in the Early Host-Pathogen Interaction of Bovine Macrophages with Divergent Mycobacterium bovis Strains in the United Kingdom

Bovine tuberculosis has been an escalating animal health issue in the United Kingdom since the 1980s, even though control policies have been in place for over 60 years. The importance of the genetics of the etiological agent, Mycobacterium bovis, in the reemergence of the disease has been largely overlooked. We compared the interaction between bovine monocyte-derived macrophages (bMDM) and two M. bovis strains, AF2122/97 and G18, representing distinct genotypes currently circulating in the United Kingdom. These M. bovis strains exhibited differences in survival and growth in bMDM. Although uptake was similar, the number of viable intracellular AF2122/97 organisms increased rapidly, while G18 growth was constrained for the first 24 h. AF2122/97 infection induced a greater transcriptional response by bMDM than G18 infection with respect to the number of differentially expressed genes and the fold changes measured. AF2122/97 infection induced more bMDM cell death, with characteristics of necrosis and apoptosis, more inflammasome activation, and a greater type I interferon response than G18. In conclusion, the two investigated M. bovis strains interact in significantly different ways with the host macrophage. In contrast to the relatively silent infection by G18, AF2122/97 induces greater signaling to attract other immune cells and induces host cell death, which may promote secondary infections of naive macrophages. These differences may affect early events in the host-pathogen interaction, including granuloma development, which could in turn alter the progression of the disease. Therefore, the potential involvement of M. bovis genotypes in the reemergence of bovine tuberculosis in the United Kingdom warrants further investigation.

Bovine tuberculosis in Ethiopia: A systematic review and meta-analysis

Bovine tuberculosis (bTB) is a known endemic disease of cattle in Ethiopia; however, there is lack of a comprehensive information on the status and distribution of the disease in the country. The objectives of this systematic review and meta-analysis were to provide a pooled prevalence estimate of bTB at a national level, assess the level of in-between variance among study reports and illustrate the spatial distribution pattern in the country. Articles published on bTB from January 2000 to December, 2016 in English language were included in the review. Pubmed, CAB direct, AJOL and Web of Science were the databases used in electronic search. A total of 127 articles were retrieved from online sources, of which 56 articles were selected for data extraction based on the specified inclusion criteria. From these selected published articles, 114 animal level data were extracted for quantitative analysis. A pooled prevalence estimate of bovine tuberculosis in Ethiopia was found to be 5.8% (95% CI: 4.5, 7.5). In a multivariable meta-regression analysis, breed and production system explained 40.9% of the explainable proportion of the in-between study variance computed. The prevalence of bovine tuberculosis in Holstein-Friesians, 21.6% (95% CI: 14.7-30.7), was higher than the prevalence in local zebus 4.1 (95% CI: 3.4-4.9). Cattle kept under intensive and semi-intensive production systems had higher prevalence, 16.6% (95% CI: 12.4-21.6), of bTB than those kept in extensive livestock production system, 4.6 (95% CI: 3.4-6.2). Bovine tuberculosis is widely distributed across major livestock producing regions of Ethiopia. However, no valid data could be retrieved from Benishanul-Gumuz, Harari and Dire Dawa. Data obtained on bTB from Somali and Gambella regional states are also few and further studies are suggested in these regions. In conclusion, this review showed that bTB in cattle in Ethiopia is widespread with high prevalence in intensive and semi-intensive management systsems that keep exotic breeds and their crosses in urban and peri-urban areas. Thus, it is suggested that the design and implementation of bTB control strategies in Ethiopia should prioritize these hotspots in order to reduce the impact of the disease on the growing dairy sector.

Fine-mapping host genetic variation underlying outcomes to Mycobacterium bovis infection in dairy cows

BACKGROUND

Susceptibility to Mycobacterium bovis infection in cattle is governed in part by host genetics. However, cattle diagnosed as infected with M. bovis display varying signs of pathology. The variation in host response to infection could represent a continuum since time of exposure or distinct outcomes due to differing pathogen handling. The relationships between host genetics and variation in host response and pathological sequelae following M. bovis infection were explored by genotyping 1966 Holstein-Friesian dairy cows at 538,231 SNPs with three distinct phenotypes. These were: single intradermal cervical comparative tuberculin (SICCT) test positives with visible lesions (VLs), SICCT-positives with undetected visible lesions (NVLs) and matched controls SICCT-negative on multiple occasions.

RESULTS

Regional heritability mapping identified three loci associated with the NVL phenotype on chromosomes 17, 22 and 23, distinct to the region on chromosome 13 associated with the VL phenotype. The region on chromosome 23 was at genome-wide significance and candidate genes overlapping the mapped window included members of the bovine leukocyte antigen class IIb region, a complex known for its role in immunity and disease resistance. Chromosome heritability analysis attributed variance to six and thirteen chromosomes for the VL and NVL phenotypes, respectively, and four of these chromosomes were found to explain a proportion of the phenotypic variation for both the VL and NVL phenotype. By grouping the M. bovis outcomes (VLs and NVLs) variance was attributed to nine chromosomes. When contrasting the two M. bovis infection outcomes (VLs vs NVLs) nine chromosomes were found to harbour heritable variation. Regardless of the case phenotype under investigation, chromosome heritability did not exceed 8% indicating that the genetic control of bTB resistance consists of variants of small to moderate effect situated across many chromosomes of the bovine genome.

CONCLUSIONS

These findings suggest the host genetics of M. bovis infection outcomes is governed by distinct and overlapping genetic variants. Thus, variation in the pathology of M. bovis infected cattle may be partly genetically determined and indicative of different host responses or pathogen handling. There may be at least three distinct outcomes following M. bovis exposure in dairy cattle: resistance to infection, infection resulting in pathology or no detectable pathology.

Genomic regions underlying susceptibility to bovine tuberculosis in Holstein-Friesian cattle

BACKGROUND

The significant social and economic loss as a result of bovine tuberculosis (bTB) presents a continuous challenge to cattle industries in the UK and worldwide. However, host genetic variation in cattle susceptibility to bTB provides an opportunity to select for resistant animals and further understand the genetic mechanisms underlying disease dynamics.

METHODS

The present study identified genomic regions associated with susceptibility to bTB using genome-wide association (GWA), regional heritability mapping (RHM) and chromosome association approaches. Phenotypes comprised de-regressed estimated breeding values of 804 Holstein-Friesian sires and pertained to three bTB indicator traits: i) positive reactors to the skin test with positive post-mortem examination results (phenotype 1); ii) positive reactors to the skin test regardless of post-mortem examination results (phenotype 2) and iii) as in (ii) plus non-reactors and inconclusive reactors to the skin tests with positive post-mortem examination results (phenotype 3). Genotypes based on the 50 K SNP DNA array were available and a total of 34,874 SNPs remained per animal after quality control.

RESULTS

The estimated polygenic heritability for susceptibility to bTB was 0.26, 0.37 and 0.34 for phenotypes 1, 2 and 3, respectively. GWA analysis identified a putative SNP on Bos taurus autosomes (BTA) 2 associated with phenotype 1, and another on BTA 23 associated with phenotype 2. Genomic regions encompassing these SNPs were found to harbour potentially relevant annotated genes. RHM confirmed the effect of these genomic regions and identified new regions on BTA 18 for phenotype 1 and BTA 3 for phenotypes 2 and 3. Heritabilities of the genomic regions ranged between 0.05 and 0.08 across the three phenotypes. Chromosome association analysis indicated a major role of BTA 23 on susceptibility to bTB.

CONCLUSION

Genomic regions and candidate genes identified in the present study provide an opportunity to further understand pathways critical to cattle susceptibility to bTB and enhance genetic improvement programmes aiming at controlling and eradicating the disease.

Genetic evaluation for bovine tuberculosis resistance in dairy cattle

Genetic evaluations for resistance to bovine tuberculosis (bTB) were calculated based on British national data including individual animal tuberculin skin test results, postmortem examination (presence of bTB lesions and bacteriological culture for Mycobacterium bovis), animal movement and location information, production history, and pedigree records. Holstein cows with identified sires in herds with bTB breakdowns (new herd incidents) occurring between the years 2000 and 2014 were considered. In the first instance, cows with a positive reaction to the skin test and a positive postmortem examination were defined as infected. Values of 0 and 1 were assigned to healthy and infected animal records, respectively. Data were analyzed with mixed models. Linear and logit function heritability estimates were 0.092 and 0.172, respectively. In subsequent analyses, breakdowns were split into 2-mo intervals to better model time of exposure and infection in the contemporary group. Intervals with at least one infected individual were retained and multiple intervals within the same breakdown were included. Healthy animal records were assigned values of 0, and infected records a value of 1 in the interval of infection and values reflecting a diminishing probability of infection in the preceding intervals. Heritability and repeatability estimates were 0.115 and 0.699, respectively. Reliabilities and across time stability of the genetic evaluation were improved with the interval model. Subsequently, 2 more definitions of "infected" were analyzed with the interval model: (1) all positive skin test reactors regardless of postmortem examination, and (2) all positive skin test reactors plus nonreactors with positive postmortem examination. Estimated heritability was 0.085 and 0.089, respectively; corresponding repeatability estimates were 0.701 and 0.697. Genetic evaluation reliabilities and across time stability did not change. Correlations of genetic evaluations for bTB with other traits in the current breeding goal were mostly not different from zero. Correlation with the UK Profitable Lifetime Index was moderate, significant, and favorable. Results demonstrated the feasibility of a national genetic evaluation for bTB resistance. Selection for enhanced resistance will have a positive effect on profitability and no antagonistic effects on current breeding goal traits. Official genetic evaluations are now based on the interval model and the last bTB trait definition.

Quantitative genetic analysis of the bTB diagnostic single intradermal comparative cervical test (SICCT)

BACKGROUND

Bovine tuberculosis (bTB) is a disease of significant economic importance and is a persistent animal health problem with implications for public health worldwide. Control of bTB in the UK has relied on diagnosis through the single intradermal comparative cervical test (SICCT). However, limitations in the sensitivity of this test hinder successful eradication and the control of bTB remains a major challenge. Genetic selection for cattle that are more resistant to bTB infection can assist in bTB control. The aim of this study was to conduct a quantitative genetic analysis of SICCT measurements collected during bTB herd testing. Genetic selection for bTB resistance will be partially informed by SICCT-based diagnosis; therefore it is important to know whether, in addition to increasing bTB resistance, this might also alter genetically the epidemiological characteristics of SICCT.

RESULTS

Our main findings are that: (1) the SICCT test is robust at the genetic level, since its hierarchy and comparative nature provide substantial protection against random genetic changes that arise from genetic drift and from correlated responses among its components due to either natural or artificial selection; (2) the comparative nature of SICCT provides effective control for initial skin thickness and age-dependent differences; and (3) continuous variation in SICCT is only lowly heritable and has a weak correlation with SICCT positivity among healthy animals which was not significantly different from zero (P > 0.05). These emerging results demonstrate that genetic selection for bTB resistance is unlikely to change the probability of correctly identifying non-infected animals, i.e. the test's specificity, while reducing the overall number of cases.

CONCLUSIONS

This study cannot exclude all theoretical risks from selection on resistance to bTB infection but the role of SICCT in disease control is unlikely to be rapidly undermined, with any adverse correlated responses expected to be weak and slow, which allow them to be monitored and managed.

A review of risk factors for bovine tuberculosis infection in cattle in the UK and Ireland

Bovine tuberculosis (bTB) is an important disease of cattle caused by infection withMycobacterium bovis, a pathogen that may be extremely difficult to eradicate in the presence of a true wildlife reservoir. Our objective was to identify and review relevant literature and provide a succinct summary of current knowledge of risk factors for transmission of infection of cattle. Search strings were developed to identify publications from electronic databases to February 2015. Abstracts of 4255 papers identified were reviewed by three reviewers to determine whether the entire article was likely to contain relevant information. Risk factors could be broadly grouped as follows: animal (including nutrition and genetics), herd (including bTB and testing history), environment, wildlife and social factors. Many risk factors are inter-related and study designs often do not enable differentiation between cause and consequence of infection. Despite differences in study design and location, some risk factors are consistently identified, e.g. herd size, bTB history, presence of infected wildlife, whereas the evidence for others is less consistent and coherent, e.g. nutrition, local cattle movements. We have identified knowledge gaps where further research may result in an improved understanding of bTB transmission dynamics. The application of targeted, multifactorial disease control regimens that address a range of risk factors simultaneously is likely to be a key to effective, evidence-informed control strategies.

Single nucleotide polymorphisms in toll-like receptor genes and case-control association studies with bovine tuberculosis

Aim:

Toll-like receptor 2 (TLR2) and TLR4 genes play critical roles in host recognition of Mycobacterium bovis infection and initiation of innate and adaptive immune response. The present study was aimed at exploring the association of seven single nucleotide polymorphisms (SNPs) in TLR2 and TLR4 genes with susceptibility/resistance against bovine tuberculosis (bTB) infection in cattle.

Materials and Methods:

A case-control resource population of 35 positive and 45 negative animals was developed after screening with single intradermal tuberculin test for bTB. Resource population was screened for SNPs in TLR2 and TLR4 genes using polymerase chain reaction-restriction fragment length polymorphism. The PROC LOGISTIC procedure of SAS 9.3 was used to find an association of allelic and genotypic frequencies with bTB.

Results:

In TLR2 gene, two of SNPs under study (rs55617172 and rs68268253) revealed polymorphism while in the case of TLR4 gene all four SNPs under investigation (rs8193041, rs207836014, rs8193060, and rs8193069) were found to be polymorphic in case-control population. SNP locus rs55617172 in TLR2 gene was found significantly (p<0.01) associated with susceptibility/resistance to TB in cattle.

Conclusion:

These findings indicate the presence of SNPs in TLR2 and TLR4 genes in our resource population. Upon validation in independent, large resource population and following biological characterization, SNP rs55617172 can be incorporated in marker panel for selection of animals with greater resistance to bTB.

A genome-wide association study for genetic susceptibility to Mycobacterium bovis infection in dairy cattle identifies a susceptibility QTL on chromosome 23

Background

Bovine tuberculosis (bTB) infection in cattle is a significant economic concern in many countries, with annual costs to the UK and Irish governments of approximately €190 million and €63 million, respectively, for bTB control. The existence of host additive and non-additive genetic components to bTB susceptibility has been established.

Methods

Two approaches i.e. single-SNP (single nucleotide polymorphism) regression and a Bayesian method were applied to genome-wide association studies (GWAS) using high-density SNP genotypes (n = 597,144 SNPs) from 841 dairy artificial insemination (AI) sires. Deregressed estimated breeding values for bTB susceptibility were used as the quantitative dependent variable. Network analysis was performed using the quantitative trait loci (QTL) that were identified as significant in the single-SNP regression and Bayesian analyses separately. In addition, an identity-by-descent analysis was performed on a subset of the most prolific sires in the dataset that showed contrasting prevalences of bTB infection in daughters.

Results

A significant QTL region was identified on BTA23 (P value >1 × 10⁻⁵, Bayes factor >10) across all analyses. Sires with the minor allele (minor allele frequency = 0.136) for this QTL on BTA23 had estimated breeding values that conferred a greater susceptibility to bTB infection than those that were homozygous for the major allele. Imputation of the regions that flank this QTL on BTA23 to full sequence indicated that the most significant associations were located within introns of the FKBP5 gene.

Conclusions

A genomic region on BTA23 that is strongly associated with host susceptibility to bTB infection was identified. This region contained FKBP5, a gene involved in the TNFα/NFκ-B signalling pathway, which is a major biological pathway associated with immune response. Although there is no study that validates this region in the literature, our approach represents one of the most powerful studies for the analysis of bTB susceptibility to date.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0197-x) contains supplementary material, which is available to authorized users.

Use of Genomic Tools to Improve Cattle Health in the Context of Infectious Diseases

Although infectious diseases impose a heavy economic burden on the cattle industry, the etiology of many disorders that affect livestock is not fully elucidated, and effective countermeasures are often lacking. The main tools available until now have been vaccines, antibiotics and antiparasitic drugs. Although these have been very successful in some cases, the appearance of parasite and microbial resistance to these treatments is a cause of concern. Next-generation sequencing provides important opportunities to tackle problems associated with pathogenic illnesses. This review describes the rapid gains achieved to track disease progression, identify the pathogens involved, and map pathogen interactions with the host. Use of novel genomic tools subsequently aids in treatment development, as well as successful creation of breeding programs aimed toward less susceptible livestock. These may be important tools for mitigating the long term effects of combating infection and helping reduce the reliance on antibiotic treatment.

Relationship of bovine TNF-α gene polymorphisms with the risk of bovine tuberculosis in Holstein cattle

Many studies suggest significant genetic variation in the resistance of cattle and humans to infection with Mycobacterium bovis (M. bovis), the causative agent of zoonotic tuberculosis. TNF-α promotes inflammation and induces apoptosis in response to mycobacterial infection. The aim of the present study was to investigate the influence of single nucleotide polymorphisms of the TNF-α gene on bovine tuberculosis (bTB) susceptibility. We genotyped the TNF-α gene in 74 bTB-infected Holstein cows and 90 healthy control animals. The influence in the exon 3 region of TNF-α polymorphisms on bTB susceptibility was subsequently investigated by association analysis. Our finding demonstrated that the g.27534932A>C polymorphism of the TNF-α is associated with bTB in Holstein cattle. The susceptibility of cattle with the g.27534932A>C genotype compared with the CC genotype was 4.11-fold (95% CI, 1.27–13.36; P=0.02) higher. The g.27534932A>C polymorphism located in exon 3 of the TNF-α gene, and the functional consequence was missense. The deduced amino acid sequence for the protein product revealed an arginine to serine conversion at position 159, which may affect initiation of protein synthesis and disrupt normal TNF-α function that protects animals against mycobacterial infection. A significant association was observed with the A allele as a risk factor for bTB susceptibility (OR, 3.84; 95% CI, 1.21–12.17; P=0.02). In conclusion, this is the first report showing that the g.27534932A>C polymorphism may contribute to TNF-α-mediated bTB susceptibility.

Molecular Epidemiology of Mycobacterium bovis in Humans and Cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), is a serious re-emerging disease in both animals and humans. The evolution of the Multi- and Extensively drug-resistant M. bovis strains (MDR-TB and XDR-TB) represents a global threat to public health. Worldwide, the disease is responsible for great economic losses in the veterinary field, serious threat to the ecosystem, and about 3.1% of human TB cases, up to 16% in Tanzania. Only thorough investigation to understand the pathogen's epidemiology can help in controlling the disease and minimizing its threat. For this purpose, various tools have been developed for use in advanced molecular epidemiological studies of bTB, either alone or in combination with standard conventional epidemiological approaches. These techniques enable the analysis of the intra- and inter-species transmission dynamics of bTB. The delivered data can reveal detailed insights into the source of infection, correlations among human and bovine isolates, strain diversity and evolution, spread, geographical localization, host preference, tracing of certain virulence factors such as antibiotic resistance genes, and finally the risk factors for the maintenance and spread of M. bovis. They also allow for the determination of epidemic and endemic strains. This, in turn, has a significant diagnostic impact and helps in vaccine development for bTB eradication programs. The present review discusses many topics including the aetiology, epidemiology and importance of M. bovis, the prevalence of bTB in humans and animals in various countries, the molecular epidemiology of M. bovis, and finally applied molecular epidemiological techniques.

Relationship of bovine NOS2 gene polymorphisms to the risk of bovine tuberculosis in Holstein cattle

Many studies suggest significant genetic variation in the resistance of cattle and humans to infection with Mycobacterium bovis, the causative agent of zoonotic tuberculosis. The inducible nitric oxide synthase (iNOS which is encoded by the NOS2 gene) plays a key role in the immunological control of a broad spectrum of infectious agents. This study aimed to investigate the influence of genetic variations in the promoter of the NOS2 gene on bovine tuberculosis (bTB) susceptibility. In this study, the NOS2 genes of 74 bTB-infected Holstein cows and 90 healthy controls were genotyped using PCR followed by nucleotide sequencing. Polymorphisms at rs207692718, rs109279434, rs209895548, rs385993919, rs433717754, rs383366213, rs466730386, rs715225976, rs525673647, rs720757654 and g.19958101T>G in the promoter region of the NOS2 gene were detected. The g.19958101T>G SNP produced two different conformation patterns (TT and TG) and the TG genotype was over-represented in the bTB group (20.27%) compared with the control group (2.22%). The TG genotype frequency of the g.19958101T>G variant was significantly higher in bTB cattle than in healthy controls (OR, 11.19; 95% CI, 2.47–50.73; P=0.0002). The G allele of the g.19958101T>G polymorphism was more frequent in bTB group when compared to control group (10.14% versus 1.11%). Furthermore, the G allele was a risk factor for bTB susceptibility (OR, 10.04; 95% CI, 2.26–44.65; P=0.0002). In conclusion, the g.19958101T>G polymorphism of the NOS2 gene may contribute to the susceptibility of Holstein cattle to bTB.

Selective breeding: the future of TB management in African buffalo?

The high prevalence of bovine tuberculosis (BTB) in African buffalo (Syncerus caffer) in regions of southern African has a negative economic impact on the trade of animals and animal products, represents an ecological threat to biodiversity, and poses a health risk to local communities through the wildlife-cattle-human interface. Test and cull methods may not be logistically feasible in many free-range wildlife systems, and with the presence of co-existing BTB hosts and the limited effectiveness of the BCG vaccine in buffalo, there is a need for alternative methods of BTB management. Selective breeding for increased resistance to BTB in buffalo may be a viable method of BTB management in the future, particularly if genetic information can be incorporated into these schemes. To explore this possibility, we discuss the different strategies that can be employed in selective breeding programmes, and consider the implementation of genetic improvement schemes. We reflect on the suitability of applying this strategy for enhanced BTB resistance in African buffalo, and address the challenges of this approach that must be taken into account. Conclusions and the implications for management are presented.

Admixture mapping of tuberculosis and pigmentation-related traits in an African-European hybrid cattle population

Admixture mapping affords a powerful approach to genetic mapping of complex traits and may be particularly suited to investigation in cattle where many breeds and populations are hybrids of the two divergent ancestral genomes, derived from Bos taurus and Bos indicus. Here we design a minimal genome wide SNP panel for tracking ancestry in recent hybrids of Holstein–Friesian and local Arsi zebu in a field sample from a region of high bovine tuberculosis (BTB) endemicity in the central Ethiopian highlands. We first demonstrate the utility of this approach by mapping the red coat color phenotype, uncovering a highly significant peak over the MC1R gene and a second peak with no previously known candidate gene. Secondly, we exploit the described differential susceptibility to BTB between the ancestral strains to identify a region in which Bos taurus ancestry associates, at suggestive significance, with skin test positivity. Interestingly, this association peak contains the toll-like receptor gene cluster on chromosome 6. With this work we have shown the potential of admixture mapping in hybrid domestic animals with divergent ancestral genomes, a recurring condition in domesticated species.

Variance components for susceptibility to Mycobacterium bovis infection in dairy and beef cattle

BACKGROUND

Infection of livestock with bovine tuberculosis (bTB; Mycobacterium bovis) is of major economical concern in many countries; approximately 15 000 to 20 000 cattle are infected per year in Ireland. The objective of this study was to quantify the genetic variation for bTB susceptibility in Irish dairy and beef cattle.

METHODS

A total of 105 914 cow, 56 904 heifer and 21 872 steer single intra-dermal comparative tuberculin test records (i.e., binary trait) collected from the years 2001 to 2010 from dairy and beef herds were included in the analysis. Only animal level data pertaining to periods of herd bTB infection were retained. Variance components for bTB were estimated using animal linear and threshold mixed models and co-variances were estimated using sire linear mixed models.

RESULTS

Using a linear model, the heritability for susceptibility to bTB in the entire dataset was 0.11 and ranged from 0.08 (heifers in dairy herds) to 0.19 (heifers in beef herds) among the sub-populations investigated. Differences in susceptibility to bTB between breeds were clearly evident. Estimates of genetic correlations for bTB susceptibility between animal types (i.e., cows, heifers, steers) were all positive (0.10 to 0.64), yet different from one. Furthermore, genetic correlations for bTB susceptibility between environments that differed in herd prevalence of bTB ranged from 0.06 to 0.86 and were all different from one.

CONCLUSIONS

Genetic trends for bTB susceptibility observed in this study suggest a slight increase in genetic susceptibility to bTB in recent years. Since bTB is of economic importance and because all animals are routinely tested at least once annually in Ireland and some other countries, the presence of genetic variation for bTB susceptibility suggests that bTB susceptibility should be included in a national breeding program to halt possible deterioration in genetic susceptibility to bTB infection.

Genomics and disease resistance studies in livestock

This paper considers the application of genetic and genomic techniques to disease resistance, the interpretation of data arising from such studies and the utilisation of the research outcomes to breed animals for enhanced resistance. Resistance and tolerance are defined and contrasted, factors affecting the analysis and interpretation of field data presented, and appropriate experimental designs discussed. These general principles are then applied to two detailed case studies, infectious pancreatic necrosis in Atlantic salmon and bovine tuberculosis in dairy cattle, and the lessons learnt are considered in detail. It is concluded that the rate limiting step in disease genetic studies will generally be provision of adequate phenotypic data, and its interpretation, rather than the genomic resources. Lastly, the importance of cross-disciplinary dialogue between the animal health and animal genetics communities is stressed.

Genomic prediction for tuberculosis resistance in dairy cattle

Background

The increasing prevalence of bovine tuberculosis (bTB) in the UK and the limitations of the currently available diagnostic and control methods require the development of complementary approaches to assist in the sustainable control of the disease. One potential approach is the identification of animals that are genetically more resistant to bTB, to enable breeding of animals with enhanced resistance. This paper focuses on prediction of resistance to bTB. We explore estimation of direct genomic estimated breeding values (DGVs) for bTB resistance in UK dairy cattle, using dense SNP chip data, and test these genomic predictions for situations when disease phenotypes are not available on selection candidates.

Methodology/Principal Findings

We estimated DGVs using genomic best linear unbiased prediction methodology, and assessed their predictive accuracies with a cross validation procedure and receiver operator characteristic (ROC) curves. Furthermore, these results were compared with theoretical expectations for prediction accuracy and area-under-the-ROC-curve (AUC). The dataset comprised 1151 Holstein-Friesian cows (bTB cases or controls). All individuals (592 cases and 559 controls) were genotyped for 727,252 loci (Illumina Bead Chip). The estimated observed heritability of bTB resistance was 0.23±0.06 (0.34 on the liability scale) and five-fold cross validation, replicated six times, provided a prediction accuracy of 0.33 (95% C.I.: 0.26, 0.40). ROC curves, and the resulting AUC, gave a probability of 0.58, averaged across six replicates, of correctly classifying cows as diseased or as healthy based on SNP chip genotype alone using these data.

Conclusions/Significance

These results provide a first step in the investigation of the potential feasibility of genomic selection for bTB resistance using SNP data. Specifically, they demonstrate that genomic selection is possible, even in populations with no pedigree data and on animals lacking bTB phenotypes. However, a larger training population will be required to improve prediction accuracies.

Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis

Tuberculosis (TB) caused by Mycobacterium bovis is a re-emerging disease of livestock that is of major economic importance worldwide, as well as being a zoonotic risk. There is significant heritability for host resistance to bovine TB (bTB) in dairy cattle. To identify resistance loci for bTB, we undertook a genome-wide association study in female Holstein-Friesian cattle with 592 cases and 559 age-matched controls from case herds. Cases and controls were categorised into distinct phenotypes: skin test and lesion positive vs skin test negative on multiple occasions, respectively. These animals were genotyped with the Illumina BovineHD 700K BeadChip. Genome-wide rapid association using linear and logistic mixed models and regression (GRAMMAR), regional heritability mapping (RHM) and haplotype-sharing analysis identified two novel resistance loci that attained chromosome-wise significance, protein tyrosine phosphatase receptor T (PTPRT; P=4.8 × 10(-7)) and myosin IIIB (MYO3B; P=5.4 × 10(-6)). We estimated that 21% of the phenotypic variance in TB resistance could be explained by all of the informative single-nucleotide polymorphisms, of which the region encompassing the PTPRT gene accounted for 6.2% of the variance and a further 3.6% was associated with a putative copy number variant in MYO3B. The results from this study add to our understanding of variation in host control of infection and suggest that genetic marker-based selection for resistance to bTB has the potential to make a significant contribution to bTB control.

Bovine TB in livestock and wildlife: what's in the genes?

Bovine tuberculosis (BTB) is a chronic, infectious disease found in domestic livestock and wildlife. It is caused predominantly by Mycobacterium bovis, which forms part of the Mycobacterium tuberculosis complex. BTB has serious implications for the movement of animals and animal products, biodiversity, and public health and is of significant economic concern. The existence of wildlife maintenance hosts makes it extremely difficult to eradicate BTB, even when established control strategies are in place, creating the need for alternative methods for controlling this disease. There are multiple factors that influence the outcome of infection by a pathogen, one of which is the host's genome. The identification of genetic variants involved in the susceptibility to BTB would supply a new selection of potential drug targets as well as the possibility for the breeding of animals with greater disease resistance. In this review, we collate the results of the BTB heritability and association studies performed in cattle and wildlife, discuss considerations and other methodologies (such as gene expression work) to be taken into account when performing genetic studies, and make some recommendations for future work in this area.

Gene polymorphisms in African buffalo associated with susceptibility to bovine tuberculosis infection

Bovine tuberculosis (BTB) is a chronic, highly infectious disease that affects humans, cattle and numerous species of wildlife. In developing countries such as South Africa, the existence of extensive wildlife-human-livestock interfaces poses a significant risk of Mycobacterium bovis transmission between these groups, and has far-reaching ecological, economic and public health impacts. The African buffalo (Syncerus caffer), acts as a maintenance host for Mycobacterium bovis, and maintains and transmits the disease within the buffalo and to other species. In this study we aimed to investigate genetic susceptibility of buffalo for Mycobacterium bovis infection. Samples from 868 African buffalo of the Cape buffalo subspecies were used in this study. SNPs (n = 69), with predicted functional consequences in genes related to the immune system, were genotyped in this buffalo population by competitive allele-specific SNP genotyping. Case-control association testing and statistical analyses identified three SNPs associated with BTB status in buffalo. These SNPs, SNP41, SNP137 and SNP144, are located in the SLC7A13, DMBT1 and IL1α genes, respectively. SNP137 remained significantly associated after permutation testing. The three genetic polymorphisms identified are located in promising candidate genes for further exploration into genetic susceptibility to BTB in buffalo and other bovids, such as the domestic cow. These polymorphisms/genes may also hold potential for marker-assisted breeding programmes, with the aim of breeding more BTB-resistant animals and herds within both the national parks and the private sector.

Differential responses of macrophages from bovines naturally resistant or susceptible to Mycobacterium bovis after classical and alternative activation

It is known that macrophages from naturally resistant animals possess a strong immune response against bovine tuberculosis to control mycobacterial infections. In the present study, the macrophage phagocytic activity, intracellular bacterial survival, and cytokine gene expression induced by classical and alternative activators against Mycobacterium bovis in naturally resistant or susceptible bovines, were evaluated. Animals were classified as naturally resistant or susceptible based on the capacity of their macrophages to allow M. bovis (BCG) growth. Peripheral blood macrophages from naturally resistant and susceptible animals were activated by classical and alternative stimuli and challenged with either non-pathogenic M. bovis BCG strain or pathogenic 9926 strain. Naturally resistant animals showed the highest phagocytosis index and microbial control after classical and alternative stimuli, being this response higher against the strain 9926 than the non-virulent strain. In addition, the response of macrophages activated by the classical pathway was higher than that under the alternative activation against both types of strains. Furthermore, classical pathway-activated macrophages derived from naturally resistant animals expressed higher levels of the pro-inflammatory markers iNOS, IL-1β, TNF-α, MIP-1 and MIP-3, and the anti-inflammatory markers ARGII and TGF-b, particularly to BCG. The results of this study showed that macrophages from naturally resistant animals produced stronger pro-inflammatory responses than those from susceptible ones to signals provided by classical pathway activators. Its role in innate immunity against M. bovis is yet to be determined.

Genetic predisposition to pass the standard SICCT test for bovine tuberculosis in British cattle

Bovine tuberculosis (bTB) imposes an important financial burden on the British cattle industry, yet despite intense efforts to control its spread, incidence is currently rising. Surveillance for bTB is based on a skin test that measures an immunological response to tuberculin. Cattle that fail the test are classified as "reactors" and slaughtered. Recent studies have identified genetic markers associated with the reaction of cattle to the tuberculin test. At marker INRA111 a relatively common '22' genotype occurs significantly more frequently in non-reactor cattle. Here we test the possibility that the putative protective '22' genotype does not confer resistance but instead causes cattle that carry it to react less strongly to the prescribed test, and hence avoid slaughter, potentially even though they are infected. We show that, after controlling for age and breed, '22' cattle react less strongly to the immunological challenge and may therefore be less likely to be classified as a reactor. These results highlight the potential discrepancy between infection and test status and imply that the effectiveness of the test-and-slaughter policy may be being compromised by selection for cattle that are genetically predisposed to react less strongly to tuberculin.

Innate immune markers that distinguish red deer (Cervus elaphus) selected for resistant or susceptible genotypes for Johne's disease

While many factors contribute to resistance and susceptibility to infectious disease, a major component is the genotype of the host and the way in which it is expressed. Johne’s disease is a chronic inflammatory bowel disease affecting ruminants and is caused by infection with Mycobacterium avium subspecies paratuberculosis (MAP). We have previously identified red deer breeds (Cervus elaphus) that are resistant; have a low rate of MAP infection and do not progress to develop Johne’s disease. In contrast, susceptible breeds have a high rate of MAP infection as seen by seroconversion and progress to develop clinical Johne’s disease. The aim of this study was to determine if immunological differences exist between animals of resistant or susceptible breeds. Macrophage cultures were derived from the monocytes of deer genotypically defined as resistant or susceptible to the development of Johne’s disease. Following in vitro infection of the cells with MAP, the expression of candidate genes was assessed by quantitative PCR as well as infection rate and cell death rate. The results indicate that macrophages from susceptible animals show a significantly higher upregulation of inflammatory genes (iNOS, IL-1α, TNF-α and IL-23p19) than the macrophages from resistant animals. Cells from resistant animals had a higher rate of apoptosis at 24 hours post infection (hpi) compared to macrophages from susceptible animals. The excessive expression of inflammatory mRNA transcripts in susceptible animals could cause inefficient clearing of the mycobacterial organism and the establishment of disease. Controlled upregulation of inflammatory pathways coupled with programmed cell death in the macrophages of resistant animals may predispose the host to a protective immune response against this mycobacterial pathogen.

Asymptomatic cattle naturally infected with Mycobacterium bovis present exacerbated tissue pathology and bacterial dissemination

Rational discovery of novel immunodiagnostic and vaccine candidate antigens to control bovine tuberculosis (bTB) requires knowledge of disease immunopathogenesis. However, there remains a paucity of information on the Mycobacterium bovis-host immune interactions during the natural infection. Analysis of 247 naturally PPD+ M. bovis-infected cattle revealed that 92% (n = 228) of these animals were found to display no clinical signs, but presented severe as well as disseminated bTB-lesions at post-mortem examination. Moreover, dissemination of bTB-lesions positively correlated with both pathology severity score (Spearman r = 0.48; p<0.0001) and viable tissue bacterial loads (Spearman r = 0.58; p = 0.0001). Additionally, granuloma encapsulation negatively correlated with M. bovis growth as well as pathology severity, suggesting that encapsulation is an effective mechanism to control bacterial proliferation during natural infection. Moreover, multinucleated giant cell numbers were found to negatively correlate with bacterial counts (Spearman r = 0.25; p = 0.03) in lung granulomas. In contrast, neutrophil numbers in the granuloma were associated with increased M. bovis proliferation (Spearman r = 0.27; p = 0.021). Together, our findings suggest that encapsulation and multinucleated giant cells control M. bovis viability, whereas neutrophils may serve as a cellular biomarker of bacterial proliferation during natural infection. These data integrate host granuloma responses with mycobacterial dissemination and could provide useful immunopathological-based biomarkers of disease severity in natural infection with M. bovis, an important cattle pathogen.

Understanding and predicting complex traits: knowledge from cattle

The genetic architecture of complex traits in cattle includes very large numbers of loci affecting any given trait. Most of these loci have small effects but occasionally there are loci with moderate-to-large effects segregating due to recent selection for the mutant allele. Genomic markers capture most but not all of the additive genetic variance for traits, probably because there are causal mutations with low allele frequency and therefore in incomplete linkage disequilibrium with the markers. The prediction of genetic value from genomic markers can achieve high accuracy by using statistical models that include all markers and assuming that marker effects are random variables drawn from a specified prior distribution. Recent effective population size is in the order of 100 within cattle breeds and ≈ 2500 animals with genotypes and phenotypes are sufficient to predict the genetic value of animals with an accuracy of 0.65. Recent effective population size for humans is much larger, in the order of 10,000-15,000, and more than 145,000 records would be required to reach a similar accuracy for people. However, our calculations assume that genomic markers capture all the genetic variance. This may be possible in the future as causal polymorphisms are genotyped using genome sequence data.

Genes controlling vaccine responses and disease resistance to respiratory viral pathogens in cattle

Farm animals remain at risk of endemic, exotic and newly emerging viruses. Vaccination is often promoted as the best possible solution, and yet for many pathogens, either there are no appropriate vaccines or those that are available are far from ideal. A complementary approach to disease control may be to identify genes and chromosomal regions that underlie genetic variation in disease resistance and response to vaccination. However, identification of the causal polymorphisms is not straightforward as it generally requires large numbers of animals with linked phenotypes and genotypes. Investigation of genes underlying complex traits such as resistance or response to viral pathogens requires several genetic approaches including candidate genes deduced from knowledge about the cellular pathways leading to protection or pathology, or unbiased whole genome scans using markers spread across the genome. Evidence for host genetic variation exists for a number of viral diseases in cattle including bovine respiratory disease and anecdotally, foot and mouth disease virus (FMDV). We immunised and vaccinated a cattle cross herd with a 40-mer peptide derived from FMDV and a vaccine against bovine respiratory syncytial virus (BRSV). Genetic variation has been quantified. A candidate gene approach has grouped high and low antibody and T cell responders by common motifs in the peptide binding pockets of the bovine major histocompatibility complex (BoLA) DRB3 gene. This suggests that vaccines with a minimal number of epitopes that are recognised by most cattle could be designed. Whole genome scans using microsatellite and single nucleotide polymorphism (SNP) markers has revealed many novel quantitative trait loci (QTL) and SNP markers controlling both humoral and cell-mediated immunity, some of which are in genes of known immunological relevance including the toll-like receptors (TLRs). The sequencing, assembly and annotation of livestock genomes and is continuing apace. In addition, provision of high-density SNP chips should make it possible to link phenotypes with genotypes in field populations without the need for structured populations or pedigree information. This will hopefully enable fine mapping of QTL and ultimate identification of the causal gene(s). The research could lead to selection of animals that are more resistant to disease and new ways to improve vaccine efficacy.

Uses and implications of field disease data for livestock genomic and genetics studies

This paper identifies issues associated with field disease data and their implications on the interpretation of estimated genetic parameters and experimental designs. The main focus is on concepts relating to the impacts of diagnostic test properties and exposure to infection, and how exposure to infection is intricately related to within-herd epidemic dynamics. The following are raised challenges: (i) to more fully understand and describe the dynamic impacts of disease epidemics on genetic interpretations; (ii) to develop statistical methods to jointly estimate epidemiological and genetic parameters from complex epidemiological data; (iii) to develop and explore optimal experimental designs for case-control studies, exploiting field disease data. Solving these problems would add insight to both disease genetic and epidemiological studies, as well as enabling us to better select animals for increased disease resistance.