PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45470.
Published online 2017 April 6. doi:  10.1038/srep45470
PMCID: PMC5382539

Badger macrophages fail to produce nitric oxide, a key anti-mycobacterial effector molecule

Abstract

The European badger is recognised as a wildlife reservoir for bovine tuberculosis (bTB); the control of which is complex, costly and controversial. Despite the importance of badgers in bTB and the well-documented role for macrophages as anti-mycobacterial effector cells, badger macrophage (bdM[var phi]) responses remain uncharacterised. Here, we demonstrate that bdM[var phi] fail to produce nitric oxide (NO) or upregulate inducible nitric oxide synthase (iNOS) mRNA following Toll-like receptor (TLR) agonist treatment. BdM[var phi] also failed to make NO after stimulation with recombinant badger interferon gamma (bdIFNγ) or a combination of bdIFNγ and lipopolysaccharide. Exposure of bdM[var phi] to TLR agonists and/or bdIFNγ resulted in upregulated cytokine (IL1β, IL6, IL12 and TNFα) mRNA levels indicating that these critical pathways were otherwise intact. Although stimulation with most TLR agonists resulted in strong cytokine mRNA responses, weaker responses were evident after exposure to TLR9 agonists, potentially due to very low expression of TLR9 in bdM[var phi]. Both NO and TLR9 are important elements of innate immunity to mycobacteria, and these features of bdM[var phi] biology would impair their capacity to resist bTB infection. These findings have significant implications for the development of bTB management strategies, and support the use of vaccination to reduce bTB infection in badgers.

European badgers (Meles meles) are implicated as a major wildlife reservoir for bovine tuberculosis (bTB), the control of which poses substantial practical and political challenges1,2,3. In the United Kingdom (UK), attempts to control the spread of bTB cost the taxpayer c. £99 m/annum (2013–14)4, yet the problem continues to worsen and control measures have sparked great controversy5. Badger culling has been implemented in bTB endemic areas despite strong protests6 yet the efficacy of culling strategies has been questioned7,8. Vaccination of badgers with Bacillus Calmette–Guérin (BCG) is efficacious under laboratory and field conditions9,10 and a current initiative is examining how badger vaccination can control bTB11. Understanding of badger immune responses to bTB challenge is limited despite the potential for this research to better inform bTB control efforts as well as enhancing vaccine development and deployment strategies. Here we present a critical study of badger innate immunity focussing on macrophages, TLR-based pathogen recognition and outcomes of macrophage activation.

Macrophages are the principal target cell for the growth of mycobacteria and key effector cells in immunity. Efficient anti-mycobacterial responses require induction of strong nitric oxide (NO) responses and mice deficient in inducible nitric oxide synthase (iNOS) are highly susceptible to TB12. Macrophages produce NO by upregulation of the iNOS gene in response to stimulation by Toll-like receptor (TLR) agonists and/or cytokines such as interferon gamma (IFNγ)13. Mammals typically express 10–12 TLRs that recognise different pathogen-associated molecular patterns (PAMPs)14 and initiate a signalling cascade, triggering macrophages to produce cytokines (including IL1β, IL6, IL12 and TNFα) and enter an enhanced antimicrobial state (reviewed in ref. 15). In humans, single nucleotide polymorphisms in TLR1, 2 and 9 associate with increased susceptibility to TB16,17 and in murine models, TLR2 and TLR9 knockouts exhibit increased levels of TB replication18. We investigated NO production and cytokine mRNA upregulation in bdM[var phi] following activation of TLR and IFNγ pathways revealing two aspects of macrophage function in badgers relevant to bTB susceptibility and vaccination strategies.

Results

Badger macrophages do not produce nitric oxide and fail to effectively upregulate iNOS gene expression

Forty-eight hour cultured badger blood monocyte-derived adherent cells exhibited classical macrophage morphology and were phagocytic (Supplementary Fig. 1). Badger macrophages did not produce NO detectable by Griess assay following 48 hours exposure to lipopolysaccharide (LPS, Fig. 1a). In contrast, similar treatment of murine and chicken macrophages resulted in substantial NO production (Fig. 1a). Despite extensive analyses with macrophages isolated from wild badgers (n = 33) NO was never detected in response to a wide range of TLR agonists (Pam3CSK4, PAM; lipoarabinomannan, LAM; poly (I:C), PIC; flagellin, FGN; R848; ODN M362) or heat killed TB, BCG or Listeria monocytogenes (LM). To test whether this lack of NO phenotype was badger-specific or found in another Mustelid, we examined macrophages derived from ferret peripheral blood monocytes (Fig. 1a), or spleen and neither produced detectable NO after exposure to LPS, or other TLR agonists.

Figure 1
Badger macrophages do not produce NO.

The complete absence of NO may have been due to a disruption in the badger iNOS gene. To test for this possibility, we designed primers targeting conserved regions (using human, mouse, dog, ferret and giant panda iNOS genome sequence). A putative badger iNOS gene fragment was obtained by PCR, cloned and the sequence verified. This sequence was used to develop a QRTPCR, which revealed that the iNOS mRNA signal was very low in untreated and TLR agonist-treated M[var phi] or those treated with heat-killed bacterial preparations (<1000 copies/1.1 million copies of GAPDH) (Fig. 1b). In comparison, TLR agonist treated murine M[var phi] upregulate iNOS mRNA to levels equivalent to the levels of GAPDH signal19. Hence, the very low iNOS mRNA levels detected with bdM[var phi] are consistent with the total lack of NO response.

IFNγ initiates a TLR-independent pathway of NO production, which enhances LPS-induced NO production20. We therefore cloned badger IFNγ (bdIFNγ) and expressed it in HEK293T cells. Exposure of bdM[var phi] to culture medium containing 50 ng/ml bdIFNγ led to upregulation of TNFα mRNA but not iNOS mRNA (Fig. 1c), or release of NO. Similarly, we did not detect NO after exposing bdM[var phi] to a mixture of LPS and bdIFNγ, or to supernatants from Concanavalin A-stimulated badger peripheral blood lymphocytes (in which upregulation of bdIFNγ mRNA was detected by QRTPCR). Furthermore, Ficoll-purified leucocytes (a mixed population of lymphocytes, monocytes and other cells) did not produce NO or upregulate iNOS mRNA after stimulation with Concanavalin A for 48 hours.

Badgers have an intact iNOS gene but express an unusual mRNA isoform at low levels

Using a combination of RTPCR, 5′ RACE and genomic sequencing, we identified an in-frame coding sequence for an iNOS isoform with high homology to iNOS transcripts in other species (Fig. 2). Interestingly, the 5′ end of the transcript (corresponding to the first exon) was not homologous with the isoform generally considered the canonical 5′ iNOS sequence, but has high homology to a variant identified in RNA from human, mouse, dog and cow. Genomic sequencing revealed the potential for a transcript homologous to the canonical iNOS sequence, but this did not amplify using RTPCR or 5′ RACE. The predicted protein sequence of the observed transcript was highly conserved to rodent and human iNOS in structurally important areas and in the active site21 (Supplementary Fig. 2).

Figure 2
The cloned badger NOS transcript is an iNOS orthologue.

Badger macrophages upregulate cytokine mRNA in response to TLR agonists and bacterial lysates

To exclude the possibility that bdM[var phi] were simply failing to respond to TLR stimulation we examined the mRNA levels of various cytokine-encoding genes. Using the cross-species homology-based approach we cloned badger genes and developed badger-specific QRTPCR assays for GAPDH, β-actin, IL1β, TNFα, IL6, IFNγ and IL12 (Supplementary Tables 1 and 2), employing dilutions of plasmids containing the relevant gene to calculate copy numbers. Agonist-driven cytokine mRNA (IL1β, TNFα, IL6 and IL12) upregulation was evident with a wide range of TLR agonists, as with M[var phi] derived from other vertebrates13. For example, treatment of bdM[var phi] with the TLR4 agonist LPS led to large increases in target gene copy number (mean ± Standard Error for IL1β: 292635 ± 30286; TNFα: 261197 ± 46767; IL6: 331286 ± 83210; IL12: 10827 ± 2762; Fig. 3a–d). The strongest responses occurred with LPS, FGN, R848 (agonists recognised by TLR4, 5 and 7 and/or 8) and heat killed LM. The responses generated to TLR2 agonists, PAM (a synthetic triacylated lipopeptide) and LAM were readily detectable (mean increase in IL1β copy number 102882 ± 13932 and 109799 ± 61799 respectively, although lower than for LPS (292635 ± 30286)), despite using agonist concentrations that induce comparable IL1β responses in other species22,23. Responses to a heat-killed TB preparation were similar to TLR2 agonists. The TLR3 agonist, PIC, caused low but statistically significant (p < 0.01) increases in TNFα, IL6 and IL12 mRNA. BdM[var phi] did not respond to the TLR9 agonist ODN M362 (a synthetic unmethylated CpG) at 5 μg/ml, which stimulates responses in many vertebrates including humans24. Exposure of bdM[var phi] to heat-killed M. bovis BCG induced upregulation of IL1β mRNA to levels comparable with those induced by exposure to heat-killed TB (Supplementary Fig. 3).

Figure 3
Production of cytokines by badger macrophages in response to TLR agonists.

Badger macrophages express low levels of TLR9 and weakly upregulate cytokine mRNA in response to TLR9 agonists

Since TLR9 recognition is important during TB infections18 we considered it important to explore why bdM[var phi] failed to respond effectively to ODN M362. We first sequenced badger TLR9 (and TLR2 for comparison). The sequences obtained for bdTLR9 and bdTLR2 were clear orthologues of these TLR in other species and represented complete open reading frames for both molecules (Fig. 4 and Supplementary Table 1). Secondly, we compared the level of TLR2 and TLR9 mRNA in bdM[var phi] by QRTPCR (Fig. 5a). Notably the levels of TLR9 mRNA were much lower than detected for TLR2. Thirdly, we exposed bdM[var phi] to a dose titration of E. coli DNA (5–100 μg/ml) which induced a dose dependent upregulation of IL1β (Fig. 5b) although the amounts of DNA required to elicit a response were much higher than required for induction of TLR9 dependent responses in rodents and humans24,25. NO production was not detected after stimulation of bdM[var phi] with any dose of E. coli DNA.

Figure 4
The cloned badger TLR2 and TLR9 transcripts cluster with their orthologues from other species.
Figure 5
Badger macrophage TLR9 expression and response.

Discussion

Animals can adopt different strategies to resist the effects of pathogens26,27. For example, an animal might mount very strong immune responses to kill the pathogen but then suffer the effects of immune-mediated pathology. Alternatively, an animal might mount a less potent anti-pathogen immune response and tolerate the continued presence of the pathogen but avoid excessive immunopathology. One example of this strategy might involve modulation of TLR-induced responses which are critical for induction of strong immune responses but can also be responsible for life threatening immune pathology (e.g. endotoxic shock). Herein we examined the TLR response profiles of M[var phi] derived from wild European badgers. Badgers represent one of the most intensively studied wild mammals and are subject to considerable scrutiny in relation to their potential role in bovine tuberculosis.

Exposure of bdM[var phi] to a wide range of TLR agonists led to upregulation of cytokine mRNA (including IL1β, IL6, IL12 and TNFα) in a similar manner to that seen with other mammals. However, TLR activation and/or exposure to IFNγ did not result in the induction of nitric oxide (NO) or an effective upregulation of iNOS mRNA. Production of NO is a feature of the immune response of many vertebrates following treatment with TLR agonists alone or in combination with cytokines, including IFNγ28,29,30,31,32,33,34. Human monocyte-derived macrophages have been studied intensively with different results reported35. Often, low levels of NO production were reported which have been attributed to peculiarities in human iNOS regulation (reviewed in ref. 36), including the lack of the LPS/IFNγ-responsive enhancer sequence present in the rodent iNOS promoter37, or to the methylation status of promoter elements38. Whether similar factors underpin the lack of effective iNOS upregulation or NO production in badgers deserves attention in future studies. A similar lack of NO production was also evident with ferret blood-derived or splenic M[var phi] which suggests that this NO-negative phenotype may be common amongst mustelids. High levels of NO can be harmful and the lack of an NO response might reflect a “tolerance to infection” strategy adopted by mustelids. However, it is important to note that the phenotype is targeted towards NO production since exposure of bdM[var phi] to TLR agonists or bdIFNγ led to substantial upregulation of mRNA encoding a variety of cytokines, including TNFα (which can also be damaging). The badger iNOS gene contains an intact coding sequence, consequently the lack of iNOS upregulation or NO production could relate to transcriptional regulation, which varies between species38,39. Importantly, NO has been implicated as the major anti-M. bovis (bTB) killing mechanism with bovine M[var phi]40 and with M. tuberculosis in M[var phi] from other mammals41,42.

Polymorphisms in the human and bovine iNOS genes have also been linked to increased rates of M. tuberculosis and bTB respectively suggesting that variation in the NO response is field relevant43,44,45. Moreover, iNOS-deficient mice are highly susceptible to TB infection, presenting with atypical granulomas that can facilitate mycobacterial reactivation, dissemination and transmission12,46. A lack of NO production would compromise badgers’ ability to resist bTB infection and contribute to the development of atypical granulomas which have a reduced capacity to limit bTB spread.

In addition to a lack of NO, we demonstrate that bdM[var phi] express only low levels of TLR9 and exhibit low responses to TLR9 agonists. Both TLR2 and TLR9 knockout mice have enhanced susceptibility to mycobacterial infections18 indicating the involvement of TLR9 recognition in protective immunity to infection with mycobacteria. TLR9 is also involved in the recognition of the anti-TB vaccine, BCG, by murine dendritic cells (DC)47, indeed the first stimulatory CpG motifs were isolated from BCG48. TLR9 mediated responses have been shown to play a role at multiple levels during TB infections including regulating granuloma formation, production of early protective type I IFN and promoting DC cross-presentation of exogenous antigen via the endogenous MHC class I pathway49,50,51. The very low level of TLR9 in bdM[var phi] limits their capacity to respond to unmethylated CpG agonists and may contribute to the susceptibility of badgers to bTB, although many other factors are likely to affect bTB-badger interactions. It is important that future studies consider whether the low TLR9 phenotype we report for blood monocyte derived bdM[var phi] is evident in other badger immune cell populations such as DC or other M[var phi] populations. In terms of vaccination with BCG it is noteworthy that the effective dose of BCG in badgers is 10 times larger than the human dose10 and this may be due to reduced recognition via TLR9. It may be possible to enhance the efficacy of BCG in badgers by supplementation with other TLR agonists. Interestingly, whereas TLR9 mediated signals enhanced DC cross-presentation, TLR2-mediated signals reduced cross-presentation and this could reduce the magnitude of protective CD8+ T cell responses51. Hence, the balance of TLR-mediated recognition events can be critical in the outcome of vaccination and enhanced vaccine efficacy may be achieved by employing alternative adjuvant supplements that increase cross-presentation. T cell responses can be detected after vaccination of badgers with BCG52 and future studies may explore the impact of alternative vaccine formulations on the magnitude and longevity of the protective CD8+ T cell response.

It is important to reiterate that bdM[var phi] upregulate a full range of cytokines via TLR-dependent pathways, including TNFα, IL12 and IL1β53,54,55, which contribute to protective immunity during infection, or immunisation with BCG. Moreover, vaccination of badgers with BCG is effective against bTB under laboratory and field conditions9,10.

Our study makes two important contributions; firstly, it identifies molecular mechanisms that are likely to contribute to the susceptibility of badgers to infection with bTB: specifically, no NO and low TLR9-responses. Secondly, our study identifies TLR pathways that could be exploited to improve the design of adjuvants to enhance vaccine efficacy in badgers (e.g. supplementing with non-TLR9 agonists). In the absence of a suitable vaccine for cattle, vaccination of badgers represents a proven, sustainable and humane approach to the control of this disease in the field. Our study also highlights the importance and potential benefits of studying immune function in wildlife species threatened by zoonotic disease.

Methods

Ex vivo sampling

Badgers were blood sampled under sedation with ketamine hydrochloride (Zoetis) as part of an on-going study of their socio-ecology at Wytham Woods, Oxford, United Kingdom56. All badgers were released at their site of capture after sampling and full recovery. Ferrets were obtained from a recognised supplier and housed in accredited facilities before termination by Schedule 1 methods; all samples were obtained post mortem. All sampling was performed under Animals (Scientific Procedures) Act, 1986 licence and in accordance with guidance from the University of Oxford’s Animal Welfare and Ethical Review Board.

Cell isolation and culture

Blood samples were collected into lithium heparin vacutainers (BD Biosciences) and centrifuged (1500 g, 10 min, 4 °C) to isolate plasma. Peripheral blood mononuclear cells were isolated following established procedures57. Briefly, cells were isolated using Ficoll (Life Technologies), washed and resuspended in culture medium comprising phenol red-free Dulbecco’s Modified Eagle Medium (Life Technologies) containing 10% FCS (Foetal Calf Serum, Life Technologies), 2 mM Penicillin and Streptomycin (Life Technologies), and 2 mM L-Glutamine (Sigma-Aldrich). Cells were plated at 5 × 105 cells/well on 96 well plates and incubated for 24 hours prior to removal of non-adherent cells by washing with medium. Following a further 24 hours incubation cells were either mock-treated or treated with defined TLR agonists, microbial lysates (all obtained from Invivogen, concentrations given in Supplementary Table (3) or recombinant IFNγ (50 ng/ml). For NO assays and IFNγ bioactivity assays, cells were cultured under treatment conditions as described for 48 hours; medium was removed and frozen for Griess assay (Promega). For TLR agonist induced cytokine analyses, cells were stimulated with agonist for 4 hours, following initial time course studies (Supplementary Fig. 4). Points shown in figures represent replicate wells where cell lines are analysed (Fig. 1a, HD11 and RAW), and cells derived from individual badgers in all other cases.

RNA extraction and cDNA synthesis

RNA was extracted from cells frozen in RLT buffer (Qiagen) using either single columns or 96 well format columns (both Qiagen) following the manufacturer’s spin protocol. For QRT-PCR, the input RNA amount was normalised for cDNA reactions across a sample set and cDNA was generated (Life Techologies).

Badger gene sequences

Primers were designed based upon regions of target genes conserved in multiple mammalian genomes. Appropriate primer sets were used to amplify fragments of badger cytokine or TLR genes from PBMC cDNA using Myfi PCR mix (Bioline), according to manufacturer’s instructions. Products were cloned into pGEM TEasy (Promega) and sequenced using plasmid targeting (M13) or gene specific primers with BigDye chemistry (Life Technologies), according to manufacturer’s instructions. Sequences were analysed using Bioedit58, and badger specific sequence used to design primers for 5′RACE PCR. 5′RACE cDNA was generated using a SMARTer RACE kit (Clontech) according to manufacturer’s instructions. PCR products generated from RACE PCR products were cloned into pGEM TEasy, as described, or into pTarget (Promega) if longer than 1 kb, and sequenced as above (European Nucleotide Archive accession numbers are given in Supplementary Table 1). Sequences were analysed using BioEdit and MEGA59 and compared to human, mouse dog and ferret sequences to identify target regions for intron spanning QRT-PCR assays.

QRT-PCR assays

Intron spanning SYBR green QRT-PCR assays were designed using badger gene sequences (Supplementary Table 2). Assays were tested for specificity by melt curve analysis, and cloning and sequencing of amplicons. Validated assays were then used to analyse cytokine cDNA levels in agonist treated PBMC. Dilution series of linearised plasmids containing cytokine or TLR gene sequences were used to calculate copy numbers. Ct values for TLR or cytokine target amplification for each sample were adjusted using the GAPDH Ct value for the same sample, to account for variation in sampling and RNA preparation. The slopes of a plot of Ct against log10 of the plasmid dilution series present on each plate were calculated. The slopes of the respective gene of interest (GOI) and GAPDH dilution series were used to calculate GOI Ct values and adjusted for differences in input total RNA as follows: corrected Ct value = Ct + (Nt−Ct’) 3 S/S’, where Ct is the GOI Ct value, Nt is the median GAPDH Ct for all samples within an experiment, Ct’ is the GAPDH value of the individual sample, S is the GOI slope, and S’ is the GAPDH slope60. A 10-fold dilution series of linearised plasmid DNA was prepared for each target, starting at 20 million copies. Ct values for the dilution series were calculated as described for samples. Plots of 40 – Ct against copy number of the serial dilution were used to derive sample copy numbers per well using the following equation: sample copy number = y*e(ab), where “y” is the y intercept, “a” represents the slope of the dilution series and “b” is the corrected 40-Ct of the sample.

Recombinant IFNγ

The published sequence61 for full-length badger IFNγ including Kozak and signal sequences was modified to encode 2 glycine linker residues and a 6 histidine residue tag at the C terminus, and synthesised (Geneart). The sequence was subcloned into pTarget (Invivogen) using standard molecular biology techniques, and HEK293T cells were transiently transfected using TransIT-2020 (Mirus Bio) and cultured in Freestyle protein-free medium (ThermoFisher Scientific). The presence of a protein of the expected molecular weight (~18 kDa) in cell culture supernatants was confirmed by Western Blot, using Penta-His antibody (Qiagen), and quantified using a dot blot with the same reagent, by comparison with a polyhistidine-tagged protein of known concentration. It was used to treat cells at 50 ng/ml.

Statistics

Statistical analyses were performed in R 3.1.262. Responses to agonist panels (Figs 1b and and3)3) were analysed using linear mixed effects models with individual badger as a random effect, and the dose response titration of E. coli DNA (Fig. 5b) was analysed using a linear model, because a random effect on individual badger was not significant. Where necessary, data were transformed in order to meet assumptions of normality and homoscedasticity of residuals. Probabilities were calculated from linear mixed effects models using the Satterthwaite approximation for degrees of freedom63, as implemented in the lmerTest package64. Where more than one measurement was involved (Figs 1b and and3),3), a Holm–Bonferroni65 family-wise error rate correction was applied using function p.adjust(). Mixed models were built using the using the lme4 package66, and the linear model using function l m(). TLR2 and TLR9 expression levels (Fig. 5a) were compared using a two-tailed paired t test. Nitric oxide production data (Fig. 1a) was compared using a two-tailed paired t test (function t.test()) with a Holm–Bonferroni correction. IFNγ response (Fig. 1c) data could not be conformed to assumptions for parametric tests, so distributions were compared using one-tailed Kolmogorov–Smirnov tests67, with a Bonferroni correction for multiple comparisons.

Sequence alignments

NOS sequences were aligned using Clustal Omega68 and manually refined, with trimming to remove unaligned or poorly aligned stretches at the N and C termini. The sub-alignment of iNOS sequences was formatted for publication using BioEdit 7.258 and CorelDRAW X5. Badger iNOS transcript as described in the main text. Other sequence data was drawn from publically available databases (see Supplementary methods): RefSeq Release 7269; Genbank70; Ensembl release 8271 Dog genome assembly CanFam3.1 and Ferret genome assembly MusPutFur1.0. Mouse and human iNOS transcripts proved the forms most closely corresponding to the badger transcript which are present in publically available databases. Each is supported by RNASeq/EST evidence. TIR domains were identified in TLR peptide sequences using SMART72, and aligned using Clustal Omega. TIR domains were aligned in preference to whole-sequence alignments as functional constraints mean that TIR domains can be aligned more reliably than other regions of TLRs73,74.

Phylogenetic trees

MEGA 6.0659 was used for model selection and maximum likelihood tree construction, with 100 bootstrap replications. For the NOS tree (Fig. 2), a JTT + G model (Jones-Taylor-Thornton, with rates among sites gamma distributed) was employed. For the TIR domain tree (Fig. 4), an LG + G + I model (Le and Gascuel, with rates among sites gamma distributed, and with invariant sites) was used. In both cases, five discrete gamma categories were used, gaps were partially deleted (95% coverage cutoff), and a moderate branch swap filter was used. Default tree inference options were used. Altering the number of gamma categories or strength branch swap filter, or completely deleting sites with gaps, did not affect the branching of the trees.

Additional Information

How to cite this article: Bilham, K. et al. Badger macrophages fail to produce nitric oxide, a key anti-mycobacterial effector molecule. Sci. Rep. 7, 45470; doi: 10.1038/srep45470 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:

Acknowledgments

We thank The WildCRU badger project, Mike Noonan, Alison Poole, Greg Albery, Tom Hart, Kerry Walker and Steven Fiddaman for their assistance with various aspects of this study, and Professor Daniel Anthony for provision of heat killed BCG.

Footnotes

The authors declare no competing financial interests.

Author Contributions K.M., C.B. and C.N. collected blood samples. K.M., A.B. and S.G.P. performed the experimental work and data analysis. D.M. led field programme. A.L.S. conceived and led design of all work. All authors contributed to writing and reviewing of the manuscript.

References

  • Macdonald D. W., Woodroffe R. & Riordan P. In Wildlife Conservation on Farmland Volume 2: Conflict in the Countryside. (eds Macdonald D. W. & Feber R. E.) (Oxford University Press, 2015).
  • Zinsstag J., Schelling E., Roth F. & Kazwala R. In Mycobacterium bovis Infection in Animals and Humans (eds Thoen C. O., Steele J. H. & Gilsdorf M. J.) 68–83 (Blackwell, 2006).
  • Delahay R. J., Cheeseman C. L. & Clifton-Hadley R. S. Wildlife disease reservoirs: the epidemiology of Mycobacterium bovis infection in the European badger (Meles meles) and other British mammals. Tuberculosis 81, 43–49 (2001). [PubMed]
  • Department of Food and Rural Affairs. Request for information: Bovine TB control costs. www.gov.uk/government/uploads/system/uploads/attachment_data/file/323911/RFI_6505.pdf (2014).
  • Donnelly C. A. et al. . Positive and negative effects of widespread badger culling on tuberculosis in cattle. Nature 439, 843–846 (2006). [PubMed]
  • Donnelly C. A. & Woodroffe R. Epidemiology: Reduce uncertainty in UK badger culling. Nature 485, 582–582 (2012). [PubMed]
  • Donnelly C. A. & Woodroffe R. Bovine tuberculosis: Badger-cull targets unlikely to reduce TB. Nature 526, 640–640 (2015). [PubMed]
  • Carter S. P. et al. . Culling-induced social perturbation in Eurasian badgers Meles meles and the management of TB in cattle: an analysis of a critical problem in applied ecology. Proceedings of the Royal Society B 274, 2769–2777 (2007). [PMC free article] [PubMed]
  • Carter S. P. et al. . BCG Vaccination Reduces Risk of Tuberculosis Infection in Vaccinated Badgers and Unvaccinated Badger Cubs. Plos One 7, e49833 (2012). [PMC free article] [PubMed]
  • Chambers M. A. et al. . Bacillus Calmette-Guerin vaccination reduces the severity and progression of tuberculosis in badgers. Proceedings of the Royal Society B-Biological Sciences 278, 1913–1920 (2011). [PMC free article] [PubMed]
  • Wilson G. J., Carter S. P. & Delahay R. J. Advances and prospects for management of TB transmission between badgers and cattle. Veterinary Microbiology 151, 43–50 (2011). [PubMed]
  • MacMicking J. D. et al. . Identification of nitric oxide synthase as a protective locus against tuberculosis. Proc Natl Acad Sci USA 94, 5243–5248 (1997). [PubMed]
  • Bogdan C. Nitric oxide synthase in innate and adaptive immunity: an update. Trends in immunology 36, 161–178, doi: (2015).10.1016/j.it.2015.01.003 [PubMed] [Cross Ref]
  • Akira S., Uematsu S. & Takeuchi O. Pathogen Recognition and Innate Immunity. Cell 124, 783–801 (2006). [PubMed]
  • O’Garra A. et al. . The Immune Response in Tuberculosis. Annual Review of Immunology 31, 475–527 (2013). [PubMed]
  • Schurz H., Daya M., Moller M., Hoal E. G. & Salie M. TLR1, 2, 4, 6 and 9 Variants Associated with Tuberculosis Susceptibility: A Systematic Review and Meta-Analysis. Plos One 10 (2015). [PMC free article] [PubMed]
  • Zhang Y. et al. . Toll-Like Receptor -1, -2, and -6 Polymorphisms and Pulmonary Tuberculosis Susceptibility: A Systematic Review and Meta-Analysis. Plos One 8, doi: e63357 (2013). [PMC free article] [PubMed]
  • Bafica A. et al. . TLR9 regulates Th1 responses and cooperates with TLR2 in mediating optimal resistance to Mycobacterium tuberculosis. Journal of Experimental Medicine 202, 1715–1724 (2005). [PMC free article] [PubMed]
  • Tajima T. et al. . EP2 and EP4 receptors on muscularis resident macrophages mediate LPS-induced intestinal dysmotility via iNOS upregulation through cAMP/ERK signals. Am. J. Physiol. Gastrointest. Liver Physiol. 302, G524–G534 (2012). [PMC free article] [PubMed]
  • Ding A. H., Nathan C. F. & Stuehr D. J. Release of Reactive Nitrogen Intermediates and Reactive Oxygen Intermediates from Mouse Peritoneal-Macrophages - Comparison of Activating Cytokines and Evidence for Independent Production. Journal of Immunology 141, 2407–2412 (1988). [PubMed]
  • Fischmann T. O. et al. . Structural characterization of nitric oxide synthase isoforms reveals striking active-site conservation. Nature Structural Biology 6, 233–242 (1999). [PubMed]
  • Peroval M. Y., Boyd A. C., Young J. R. & Smith A. L. A critical role for MAPK signalling pathways in the transcriptional regulation of toll like receptors. PLoS ONE 8, e51243 (2013). [PMC free article] [PubMed]
  • Toshchakov V. et al. . TLR2 and TLR4 agonists stimulate unique repertoires of host resistance genes in murine macrophages: interferon-β-dependent signaling in TLR4-mediated responses. Journal of Endotoxin Research 9, 169–175 (2003). [PubMed]
  • Mutwiri G. et al. . Biological activity of immunostimulatory CpG DNA motifs in domestic animals. Veterinary Immunology and Immunopathology 91, 89–103 (2003). [PubMed]
  • Bauer S. et al. . Human TLR9 confers responsiveness to bacterial DNA via species-specific CpG motif recognition. Proceedings of the National Academy of Sciences of the United States of America 98, 9237–9242, doi: (2001).10.1073/pnas.161293498 [PubMed] [Cross Ref]
  • Babayan S. A. et al. . Wild immunology: converging on the real world. Ann NY Acad Sci 1236, 17–29, doi: (2011).10.1111/j.1749-6632.2011.06251.x [PubMed] [Cross Ref]
  • Raberg L., Sim D. & Read A. F. Disentangling genetic variation for resistance and tolerance to infectious diseases in animals. Science 318, 812–814, doi: (2007).10.1126/science.1148526 [PubMed] [Cross Ref]
  • Zhao B. Y., Collins M. T. & Czuprynski C. J. Induction of L-arginine-dependent production of nitric oxide in bovine monocytes by interferon gamma and lipopolysaccharide. Research in Veterinary Science 60, 190–192 (1996). [PubMed]
  • Piercy J., Werling D. & Coffey T. J. Differential responses of bovine macrophages to infection with bovine-specific and non-bovine specific mycobacteria. Tuberculosis 87, 415–420 (2007). [PubMed]
  • Drake M. G. et al. . Toll-like Receptor 7 Rapidly Relaxes Human Airways. American Journal of Respiratory and Critical Care Medicine 188, 664–672 (2013). [PMC free article] [PubMed]
  • Jesch N. K. et al. . Expression of inducible nitric oxide synthase and formation of nitric oxide by alveolar macrophages: An interspecies comparison. Environmental Health Perspectives 105, 1297–1300 (1997). [PMC free article] [PubMed]
  • Geller D. A. et al. . Molecular cloning and expression of inducible nitric oxide synthase from human hepatocytes. Proc Natl Acad Sci USA 90, 3491–3495 (1993). [PubMed]
  • Tunctan B. et al. . Comparison of nitric oxide production by monocyte/macrophages in healthy subjects and patients with active pulmonary tuberculosis. Pharmacological Research 37, 219–226 (1998). [PubMed]
  • Weinberg J. B. Nitric oxide production and nitric oxide synthase type 2 expression by human mononuclear phagocytes: a review. Mol Med 4, 557–591 (1998). [PMC free article] [PubMed]
  • Thomas A. C. & Mattila J. T. “Of mice and men”: arginine metabolism in macrophages. Frontiers in immunology 5, 479, doi: (2014).10.3389/fimmu.2014.00479 [PMC free article] [PubMed] [Cross Ref]
  • Pautz A. et al. . Regulation of the expression of inducible nitric oxide synthase. Nitric oxide: biology and chemistry 23, 75–93, doi: (2010).10.1016/j.niox.2010.04.007 [PubMed] [Cross Ref]
  • Alley E. W., Murphy W. J. & Russell S. W. A classical enhancer element responsive to both lipopolysaccharide and interferon-gamma augments induction of the iNOS gene in mouse macrophages. Gene 158, 247–251 (1995). [PubMed]
  • Gross T. J. et al. . Epigenetic silencing of the human NOS2 gene: rethinking the role of nitric oxide in human macrophage inflammatory responses. Journal of Immunology 192, 2326–2338, doi: (2014).10.4049/jimmunol.1301758 [PMC free article] [PubMed] [Cross Ref]
  • Rao K. M. K. Molecular mechanisms regulating iNOS expression in various cell types. Journal of Toxicology and Environmental Health-Part B-Critical Reviews 3, 27–58 (2000). [PubMed]
  • Esquivel-Solís H. et al. . Nitric Oxide Not Apoptosis Mediates Differential Killing of Mycobacterium bovis in Bovine Macrophages. PLoS ONE 8, e63464 (2013). [PMC free article] [PubMed]
  • Chan J., Xing Y., Magliozzo R. S. & Bloom B. R. Killing of Virulent Mycobacterium-Tuberculosis by Reactive Nitrogen Intermediates Produced by Activated Murine Macrophages. Journal of Experimental Medicine 175, 1111–1122, doi: (1992).10.1084/jem.175.4.1111 [PMC free article] [PubMed] [Cross Ref]
  • Thoma-Uszynski S. et al. . Induction of direct antimicrobial activity through mammalian toll-like receptors. Science 291, 1544–1547, doi: (2001).10.1126/science.291.5508.1544 [PubMed] [Cross Ref]
  • Cheng Y., Huang C. & Tsai H. J. Relationship of bovine NOS2 gene polymorphisms to the risk of bovine tuberculosis in Holstein cattle. Journal of Veterinary Medical Science 78, 281–286, doi: (2016).10.1292/jvms.15-0295 [PMC free article] [PubMed] [Cross Ref]
  • Gomez L. M. et al. . A polymorphism in the inducible nitric oxide synthase gene is associated with tuberculosis. Tuberculosis (Edinb) 87, 288–294, doi: (2007).10.1016/j.tube.2007.03.002 [PubMed] [Cross Ref]
  • Velez D. R. et al. . NOS2A, TLR4, and IFNGR1 interactions influence pulmonary tuberculosis susceptibility in African-Americans. Human Genetics 126, 643–653, doi: (2009).10.1007/s00439-009-0713-y [PMC free article] [PubMed] [Cross Ref]
  • Reece S. T. et al. . Serine protease activity contributes to control of Mycobacterium tuberculosis in hypoxic lung granulomas in mice. J Clin Invest 120, 3365–3376 (2010). [PMC free article] [PubMed]
  • von Meyenn F. et al. . Toll-like receptor 9 contributes to recognition of Mycobacterium bovis Bacillus Calmette-Guerin by Flt3-ligand generated dendritic cells. Immunobiology 211, 557–565, doi: (2006).10.1016/j.imbio.2006.05.004 [PubMed] [Cross Ref]
  • Krieg A. M. et al. . CpG motifs in bacterial DNA trigger direct B-cell activation. Nature 374, 546–549, doi: (1995).10.1038/374546a0 [PubMed] [Cross Ref]
  • Ito T. et al. . TLR9 regulates the mycobacteria-elicited pulmonary granulomatous immune response in mice through DC-derived Notch ligand delta-like 4. Journal of Clinical Investigation 119, 33–46, doi: (2009).10.1172/JCI35647 [PMC free article] [PubMed] [Cross Ref]
  • Kuchtey J., Fulton S. A., Reba S. M., Harding C. V. & Boom W. H. Interferon-alphabeta mediates partial control of early pulmonary Mycobacterium bovis bacillus Calmette-Guerin infection. Immunology 118, 39–49, doi: (2006).10.1111/j.1365-2567.2006.02337.x [PubMed] [Cross Ref]
  • Simmons D. P. et al. . Mycobacterium tuberculosis and TLR2 agonists inhibit induction of type I IFN and class I MHC antigen cross processing by TLR9. Journal of Immunology 185, 2405–2415, doi: (2010).10.4049/jimmunol.0904005 [PMC free article] [PubMed] [Cross Ref]
  • Lesellier S. et al. . The safety and immunogenicity of Bacillus Calmette-Guerin (BCG) vaccine in European badgers (Meles meles). Vet Immunol Immunopathol 112, 24–37, doi: (2006).10.1016/j.vetimm.2006.03.009 [PubMed] [Cross Ref]
  • Cooper A. M., Magram J., Ferrante J. & Orme I. M. Interleukin 12 (IL-12) is crucial to the development of protective immunity in mice intravenously infected with Mycobacterium tuberculosis. Journal of Experimental Medicine 186, 39–45 (1997). [PMC free article] [PubMed]
  • Flynn J. L. et al. . Tumor-Necrosis-Factor-Alpha Is Required in the Protective Immune-Response against Mycobacterium-Tuberculosis in Mice. Immunity 2, 561–572 (1995). [PubMed]
  • Kleinnijenhuis J., Oosting M., Joosten L. A. B., Netea M. G. & Van Crevel R. Innate Immune Recognition of Mycobacterium tuberculosis. Clinical & Developmental Immunology (2011). [PMC free article] [PubMed]
  • Macdonald D. W., Newman C. & Buesching C. D. In Wildlife Conservation on Farmland Volume 2: Conflict in the Coutryside. (eds Macdonald D. W. & Feber R. E.) 65 (Oxford University Press, 2015).
  • Wigley P., Hulme S. D., Bumstead N. & Barrow P. A. In vivo and in vitro studies of genetic resistance to systemic salmonellosis in the chicken encoded by the SAL1 locus. Microbes and Infection 4, 1111–1120 (2002). [PubMed]
  • Hall T. A. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic acids symposium series 41, 95–98 (1999).
  • Tamura K., Stecher G., Peterson D., Filipski A. & Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Molecular Biology and Evolution 30, 2725–2729 (2013). [PMC free article] [PubMed]
  • Wu Z. et al. . Generation and characterization of chicken bone marrow‐derived dendritic cells. Immunology 129, 133–145 (2010). [PubMed]
  • Dalley D. J., Hogarth P. J., Hughes S., Hewinson R. G. & Chambers M. A. Cloning and sequencing of badger (Meles meles) interferon gamma and its detection in badger lymphocytes. Veterinary Immunology and Immunopathology 101, 19–30 (2004). [PubMed]
  • R Core Team. R: A language and enviroment for statistical computing. R Foundation for Statistical Computing (2014).
  • Satterthwaite F. E. An approximate distribution of estimates of variance components. Biometrics bulletin, 110–114 (1946). [PubMed]
  • Kuznetsova A., Brockhoff P. & Christensen R. lmerTest: tests in linear mixed effects models. R package version 2.0-29 (2015).
  • Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat 65–70 (1979).
  • Bates D., Maechler M., Bolker B. & Walker S. lme4: Linear mixed-effects models using Eigen and S4. R package version 1.4 (2013).
  • Smirnov N. Table for estimating the goodness of fit of empirical distributions. Ann Math Stat 279–281 (1948).
  • Sievers F. et al. . Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Molecular Systems Biology 7 (2011). [PMC free article] [PubMed]
  • Pruitt K. D. et al. . RefSeq: an update on mammalian reference sequences. Nucleic Acids Research 42, D756–D763 (2014). [PMC free article] [PubMed]
  • Benson D. A. et al. . GenBank. Nucleic Acids Research 41, D36–D42 (2013). [PMC free article] [PubMed]
  • Cunningham F. et al. . Ensembl 2015. Nucleic Acids Research 43, D662–D669 (2015). [PMC free article] [PubMed]
  • Letunic I., Doerks T. & Bork P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Research 43, D257–D260 (2015). [PMC free article] [PubMed]
  • Roach J. C. et al. . The evolution of vertebrate Toll-like receptors. Proceedings of the National Academy of Sciences of the United States of America 102, 9577–9582 (2005). [PubMed]
  • Bateman A. et al. . UniProt: a hub for protein information. Nucleic Acids Research 43, D204–D212 (2015). [PMC free article] [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group