Can We Breed Cattle for Lower Bovine TB Infectivity?

Integrative and comparative genomic analyses of mammalian macrophage responses to intracellular mycobacterial pathogens

Mycobacterium tuberculosis, the causative agent of human tuberculosis (hTB), is a close evolutionary relative of Mycobacterium bovis, which causes bovine tuberculosis (bTB), one of the most damaging infectious diseases to livestock agriculture. Previous studies have shown that the pathogenesis of bTB disease is comparable to hTB disease, and that the bovine and human alveolar macrophage (bAM and hAM, respectively) transcriptomes are extensively reprogrammed in response to infection with these intracellular mycobacterial pathogens. In this study, a multi-omics integrative approach was applied with functional genomics and GWAS data sets across the two primary hosts (Bos taurus and Homo sapiens) and both pathogens (M. bovis and M. tuberculosis). Four different experimental infection groups were used: 1) bAM infected with M. bovis, 2) bAM infected with M. tuberculosis, 3) hAM infected with M. tuberculosis, and 4) human monocyte-derived macrophages (hMDM) infected with M. tuberculosis. RNA-seq data from these experiments 24 h post-infection (24 hpi) was analysed using three computational pipelines: 1) differentially expressed genes, 2) differential gene expression interaction networks, and 3) combined pathway analysis. The results were integrated with high-resolution bovine and human GWAS data sets to detect novel quantitative trait loci (QTLs) for resistance to mycobacterial infection and resilience to disease. This revealed common and unique response macrophage pathways for both pathogens and identified 32 genes (12 bovine and 20 human) significantly enriched for SNPs associated with disease resistance, the majority of which encode key components of the NF-κB signalling pathway and that also drive formation of the granuloma.

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.

Disturbance Ecology Meets Bovine Tuberculosis (bTB) Epidemiology: A Before-and-After Study on the Association between Forest Clearfelling and bTB Herd Risk in Cattle Herds

Disturbance ecology refers to the study of discrete processes that disrupt the structure or dynamics of an ecosystem. Such processes can, therefore, affect wildlife species ecology, including those that are important pathogen hosts. We report on an observational before-and-after study on the association between forest clearfelling and bovine tuberculosis (bTB) herd risk in cattle herds, an episystem where badgers (Meles meles) are the primary wildlife spillover host. The study design compared herd bTB breakdown risk for a period of 1 year prior to and after exposure to clearfelling across Ireland at sites cut in 2015–2017. The percent of herds positive rose from 3.47% prior to clearfelling to 4.08% after exposure. After controlling for confounders (e.g., herd size, herd type), we found that cattle herds significantly increased their odds of experiencing a bTB breakdown by 1.2-times (95%CIs: 1.07–1.36) up to 1 year after a clearfell risk period. Disturbance ecology of wildlife reservoirs is an understudied area with regards to shared endemic pathogens. Epidemiological observational studies are the first step in building an evidence base to assess the impact of such disturbance events; however, such studies are limited in inferring the mechanism for any changes in risk observed. The current cohort study suggested an association between clearfelling and bTB risk, which we speculate could relate to wildlife disturbance affecting pathogen spillback to cattle, though the study design precludes causal inference. Further studies are required. However, ultimately, integration of epidemiology with wildlife ecology will be important for understanding the underlying mechanisms involved, and to derive suitable effective management proposals, if required.

Review: Livestock disease resilience: from individual to herd level

Infectious diseases are a major threat to the sustainable production of high-producing animals. Control efforts, such as vaccination or breeding approaches often target improvements to individual resilience to infections, i.e., they strengthen an animal's ability to cope with infection, rather than preventing infection per se. There is increasing evidence for the contribution of non-clinical carriers (animals that become infected and are infectious but do not develop clinical signs) to the overall health and production of livestock populations for a wide range of infectious diseases. Therefore, we strongly advocate a shift of focus from increasing the disease resilience of individual animals to herd disease resilience as the appropriate target for sustainable disease control in livestock. Herd disease resilience not only captures the direct effects of vaccination or host genetics on the health and production performance of individuals but also the indirect effects on the environmental pathogen load that herd members are exposed to. For diseases primarily caused by infectious pathogens shed by herd members, these indirect effects on herd resilience are mediated both by individual susceptibility to infection and by characteristics (magnitude of infectiousness, duration of infectious period) that influence pathogen shedding from infected individuals. We review what is currently known about how vaccination and selective breeding affect herd disease resilience and its underlying components, and outline the changes required for improvement. To this purpose, we also seek to clarify and harmonise the terminology used in the different animal science disciplines to facilitate future collaborative approaches to infectious disease control in livestock.

Identification of loci associated with susceptibility to Mycobacterium avium subsp. paratuberculosis infection in Holstein cattle using combinations of diagnostic tests and imputed whole-genome sequence data

Bovine paratuberculosis (PTB) is a chronic inflammatory disease caused by Mycobacterium avium susbp. paratuberculosis (MAP). Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) significantly associated with susceptibility to bovine PTB. The main objective of this study was to identify quantitative trait loci (QTLs) associated with MAP infection in Spanish Holstein cows (N = 983) using combinations of diagnostic tests and imputed whole-genome sequence (WGS) data. The infection status of these animals was defined by three diagnostic methods including ELISA for MAP-antibodies detection, and tissue culture and PCR for MAP detection. The 983 cows included in this study were genotyped with the Bovine MD SNP50 Bead Chip, and the corresponding genotypes were imputed to WGS using the 1,000 Bull genomes reference population. In total, 33.77 million SNP variants per animal were identified across the genome. Linear mixed models were used to calculate the heritability (h²) estimates for each diagnostic test and test combinations. Next, we performed a case-control GWAS using the imputed WGS datasets and the phenotypes and combinations of phenotypes with h² estimates > 0.080. After performing the GWAS, the test combinations that showed SNPs with a significant association (P_FDR ≤ 0.05), were the ELISA-tissue PCR-tissue culture, ELISA-tissue culture, and ELISA-tissue PCR. A total of twelve quantitative trait loci (QTLs) highly associated with MAP infection status were identified on the Bos taurus autosomes (BTA) 4, BTA5, BTA11, BTA12, BTA14, BTA23, BTA24, and BTA28, and some of these QTLs were linked to immune-modulating genes. The identified QTLs on BTA23 spanning from 18.81 to 22.95 Mb of the Bos taurus genome overlapped with several QTLs previously found to be associated with PTB susceptibility, bovine tuberculosis susceptibility, and clinical mastitis. The results from this study provide more clues regarding the molecular mechanisms underlying susceptibility to PTB infection in cattle and might be used to develop national genetic evaluations for PTB in Spain.

Integrative genomics of the mammalian alveolar macrophage response to intracellular mycobacteria

Background

Bovine TB (bTB), caused by infection with Mycobacterium bovis, is a major endemic disease affecting global cattle production. The key innate immune cell that first encounters the pathogen is the alveolar macrophage, previously shown to be substantially reprogrammed during intracellular infection by the pathogen. Here we use differential expression, and correlation- and interaction-based network approaches to analyse the host response to infection with M. bovis at the transcriptome level to identify core infection response pathways and gene modules. These outputs were then integrated with genome-wide association study (GWAS) data sets to enhance detection of genomic variants for susceptibility/resistance to M. bovis infection.

Results

The host gene expression data consisted of RNA-seq data from bovine alveolar macrophages (bAM) infected with M. bovis at 24 and 48 h post-infection (hpi) compared to non-infected control bAM. These RNA-seq data were analysed using three distinct computational pipelines to produce six separate gene sets: 1) DE genes filtered using stringent fold-change and P-value thresholds (DEG-24: 378 genes, DEG-48: 390 genes); 2) genes obtained from expression correlation networks (CON-24: 460 genes, CON-48: 416 genes); and 3) genes obtained from differential expression networks (DEN-24: 339 genes, DEN-48: 495 genes). These six gene sets were integrated with three bTB breed GWAS data sets by employing a new genomics data integration tool—gwinteR. Using GWAS summary statistics, this methodology enabled detection of 36, 102 and 921 prioritised SNPs for Charolais, Limousin and Holstein-Friesian, respectively.

Conclusions

The results from the three parallel analyses showed that the three computational approaches could identify genes significantly enriched for SNPs associated with susceptibility/resistance to M. bovis infection. Results indicate distinct and significant overlap in SNP discovery, demonstrating that network-based integration of biologically relevant transcriptomics data can leverage substantial additional information from GWAS data sets. These analyses also demonstrated significant differences among breeds, with the Holstein-Friesian breed GWAS proving most useful for prioritising SNPS through data integration. Because the functional genomics data were generated using bAM from this population, this suggests that the genomic architecture of bTB resilience traits may be more breed-specific than previously assumed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07643-w.

Occupational exposure and challenges in tackling M. bovis at human-animal interface: a narrative review

Zoonotic tuberculosis caused by Mycobacterium bovis (M. bovis), a member of Mycobacterium tuberculosis complex (MTBC) has increasingly gathered attention as a public health risk, particularly in developing countries with higher disease prevalence. M. bovis is capable of infecting multiple hosts encompassing a number of domestic animals, in particular cattle as well as a broad range of wildlife reservoirs. Humans are the incidental hosts of M. bovis whereby its transmission to humans is primarily through the consumption of cattle products such as unpasteurized milk or raw meat products that have been contaminated with M. bovis or the transmission could be due to close contact with infected cattle. Also, the transmission could occur through aerosol inhalation of infective droplets or infected body fluids or tissues in the presence of wound from infected animals. The zoonotic risk of M. bovis in humans exemplified by miscellaneous studies across different countries suggested the risk of occupational exposure towards M. bovis infection, especially those animal handlers that have close and unreserved contact with cattle and wildlife populations These animal handlers comprising of livestock farmers, abattoir workers, veterinarians and their assistants, hunters, wildlife workers as well as other animal handlers are at different risk of contracting M. bovis infection, depending on the nature of their jobs and how close is their interaction with infected animals. It is crucial to identify the underlying transmission risk factors and probable transmission pathways involved in the zoonotic transmission of M. bovis from animals to humans for better designation and development of specific preventive measures and guidelines that could reduce the risk of transmission and to protect these different occupational-related/populations at risk. Effective control and disease management of zoonotic tuberculosis caused by M. bovis in humans are also hindered by various challenges and factors involved at animal–human interface. A closer look into factors affecting proper disease control and management of M. bovis are therefore warranted. Hence, in this narrative review, we have gathered a number of different studies to highlight the risk of occupational exposure to M. bovis infection and addressed the limitations and challenges underlying this context. This review also shed lights on various components and approaches in tackling M. bovis infection at animal–human interface.

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.

Characterization of potential superspreader farms for bovine tuberculosis: A review

Background

Variation in host attributes that influence their contact rates and infectiousness can lead some individuals to make disproportionate contributions to the spread of infections. Understanding the roles of such ‘superspreaders’ can be crucial in deciding where to direct disease surveillance and controls to greatest effect. In the epidemiology of bovine tuberculosis (bTB) in Great Britain, it has been suggested that a minority of cattle farms or herds might make disproportionate contributions to the spread of Mycobacterium bovis, and hence might be considered ‘superspreader farms’.

Objectives and Methods

We review the literature to identify the characteristics of farms that have the potential to contribute to exceptional values in the three main components of the farm reproductive number ‐ R_f: contact rate, infectiousness and duration of infectiousness, and therefore might characterize potential superspreader farms for bovine tuberculosis in Great Britain.

Results

Farms exhibit marked heterogeneity in contact rates arising from between‐farm trading of cattle. A minority of farms act as trading hubs that greatly augment connections within cattle trading networks. Herd infectiousness might be increased by high within‐herd transmission or the presence of supershedding individuals, or infectiousness might be prolonged due to undetected infections or by repeated local transmission, via wildlife or fomites.

Conclusions

Targeting control methods on putative superspreader farms might yield disproportionate benefits in controlling endemic bovine tuberculosis in Great Britain. However, real‐time identification of any such farms, and integration of controls with industry practices, present analytical, operational and policy challenges.

Correlations between single nucleotide polymorphisms in bovine CD209, SLC11A1, SP110 and TLR2 genes and estimated breeding values for several traits in Spanish Holstein cattle

Little is known about the correlations between the genetic susceptibility/resistance to Mycobacterium avium subsp. paratuberculosis (Map) infection and the estimated breeding values for type, production and functional traits. Previously, we identified 70 combinations of five single nucleotide polymorphisms (SNPs) in four bovine innate immune genes (SLC11A1, SP110, TLR2, CD209) that are associated with the genetic risk of paratuberculosis (PTB) or Johne's disease progression, which can be graded as low (LOWIN), latent (LATIN), or patent (PATIN) risk. Other possible combinations of these 5 SNPs were grouped in the average group (AVERIN). In the current study, differences in estimated breeding values (EBVs) for several traits were analyzed using linear models in a large cohort of Holstein cows (N = 15656) genotyped across Spain in 2016 or 2017. After the assignment of each genotyped cow to a risk group, cows within the PATIN risk group (N = 1448) had a superior combined genetic index (2797.57), type genetic index (524.62), milk yield (653.92 kg), protein yield (21.77 kg), fat yield (24.82 kg) and economic merit index (125 Euros) compared with the other three risk groups. Statistically significant differences in the longevity scores between the cows that were included in the PATIN risk group (108.85) and the LOWIN (107.82) and AVERIN (107.92) groups were also observed. The associations between the genetic risk groups and PTB diagnostic results were validated in a population of 99 cows from a Spanish farm with a high prevalence of PTB. Significant differences in ELISA readings between the PATIN (65.49 %) and the AVERIN (15.97 %), LATIN (2.11 %), and LOWIN (3.27 %) groups were observed. In addition, significant differences in Map DNA copies/gram of feces were observed between the PATIN and the other three risk groups. These results together with the substantial economic impact of PTB in dairy cattle support the selection of the animals with less susceptibility to PTB in the Spanish breeding program.

Triple Immunochromatographic System for Simultaneous Serodiagnosis of Bovine Brucellosis, Tuberculosis, and Leukemia

An immunochromatographic test system has been developed for the simultaneous rapid multiplex serodiagnostics of bovine brucellosis, tuberculosis, and leukemia. The test system is based on the use of a conjugate of gold nanoparticles with the chimeric protein Cysteine-A/G and three analytical zones with immobilized pathogen antigens: Brucella abortus lipolysaccharide, recombinant proteins MPB64 and MPB83-MPB63 of Mycobacterium bovis, and recombinant protein p24 of the bovine leukemia virus. Prototypes of the test system were tested on 98 samples of sera from healthy and infected animals. The diagnostic sensitivity of the developed test system was 92% for brucellosis, 92% for tuberculosis, and 96% for leukemia. False positive test results were not observed.

Predicting vaccine effectiveness in livestock populations: A theoretical framework applied to PRRS virus infections in pigs

Vaccines remain one of the main tools to control infectious diseases in domestic livestock. Although a plethora of veterinary vaccines are on the market and routinely applied to protect animals against infection with particular pathogens, the disease in question often continues to persist, sometimes at high prevalence. The limited effectiveness of certain vaccines in the field leaves open questions regarding the required properties that an effective vaccine should have, as well as the most efficient vaccination strategy for achieving the intended goal of vaccination programmes. To date a systematic approach for studying the combined effects of different types of vaccines and vaccination strategies is lacking. In this paper, we develop a theoretical framework for modelling the epidemiological consequences of vaccination with imperfect vaccines of various types, administered using different strategies to herds with different replacement rates and heterogeneity in vaccine responsiveness. Applying the model to the Porcine Reproductive and Respiratory Syndrome (PRRS), which despite routine vaccination remains one of the most significant endemic swine diseases worldwide, we then examine the influence of these diverse factors alone and in combination, on within-herd virus transmission. We derive threshold conditions for preventing infection invasion in the case of imperfect vaccines inducing limited sterilizing immunity. The model developed in this study has practical implications for the development of vaccines and vaccination programmes in livestock populations not only for PRRS, but also for other viral infections primarily transmitted by direct contact.

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.

Genetic differences in host infectivity affect disease spread and survival in epidemics

Survival during an epidemic is partly determined by host genetics. While quantitative genetic studies typically consider survival as an indicator for disease resistance (an individual's propensity to avoid becoming infected or diseased), mortality rates of populations undergoing an epidemic are also affected by endurance (the propensity of diseased individual to survive the infection) and infectivity (i.e. the propensity of an infected individual to transmit disease). Few studies have demonstrated genetic variation in disease endurance, and no study has demonstrated genetic variation in host infectivity, despite strong evidence for considerable phenotypic variation in this trait. Here we propose an experimental design and statistical models for estimating genetic diversity in all three host traits. Using an infection model in fish we provide, for the first time, direct evidence for genetic variation in host infectivity, in addition to variation in resistance and endurance. We also demonstrate how genetic differences in these three traits contribute to survival. Our results imply that animals can evolve different disease response types affecting epidemic survival rates, with important implications for understanding and controlling epidemics.