Impact Factor 2.029 | CiteScore 2.20
More on impact ›

Original Research ARTICLE

Front. Vet. Sci., 14 February 2019 | https://doi.org/10.3389/fvets.2019.00030

Is There a Relationship Between Bovine Tuberculosis (bTB) Herd Breakdown Risk and Mycobacterium avium subsp. paratuberculosis Status? An Investigation in bTB Chronically and Non-chronically Infected Herds

  • 1Veterinary Science Division, Agri-food and Biosciences Institute, Belfast, United Kingdom
  • 2School of Biological Sciences, Queen's University Belfast, Belfast, United Kingdom
  • 3Animal Health Ireland, Carrick on Shannon, Ireland
  • 4Animal Health and Welfare Northern Ireland, Dungannon, United Kingdom

Background: Bovine tuberculosis (bTB; Mycobacterium bovis) remains a significant problem in a number of countries, and is often found where M. avium subsp. paratuberculosis (MAP) is also present. In the United Kingdom, bTB has been difficult to eradicate despite long-term efforts. Co-infection has been proposed as one partial mechanism thwarting eradication.

Methods: A retrospective case-control study of 4,500 cattle herds in Northern Ireland, where serological testing of cattle for MAP, was undertaken (2004–2015). Blood samples were ELISA tested for MAP; infection of M. bovis was identified in herds by the comparative tuberculin test (CTT) and through post-mortem evidence of infection. Case-herds were those experiencing a confirmed bTB breakdown; control-herds were not experiencing a breakdown episode at the time of MAP testing. A second model included additional testing data of feces samples (culture and PCR results) to better inform herd MAP status. Multi-level hierarchical models were developed, controlling for selected confounders. A sensitivity analysis of the effect of MAP sample numbers per event and the prior timing of tuberculin-testing was undertaken.

Results: 45.2% (n = 250) of case observations and 36.0% (3,480) of control observations were positive to MAP by ELISA (45.8% and 36.4% when including ancillary fecal testing, respectively). Controlling for known confounders, the adjusted odds ratio (aOR) for this association was 1.339 (95%CI:1.085–1.652; including ancillary data aOR:1.356;95%CI:1.099–1.673). The size-effect of the association increased with the increasing number of samples per event used to assign herd MAP status (aOR:1.883 at >2 samples, to aOR:3.863 at >10 samples), however the estimated CI increased as N decreased. 41.7% of observations from chronic herds were MAP serology-positive and 32.2% from bTB free herds were MAP positive (aOR: 1.170; 95%ci: 0.481–2.849).

Discussion: Cattle herds experiencing a bTB breakdown were associated with increased risk of having a positive MAP status. Chronic herds tended to exhibit higher risk of a positive MAP status than bTB free herds, however there was less support for this association when controlling for repeated measures and confounding. MAP co-infection may be playing a role in the success of bTB eradiation schemes, however further studies are required to understand the mechanisms and to definitively establish causation.

Introduction

Chronic mycobacterial infections are some of the most challenging pathogens to eradicate at national and international scales, requiring large-scale long-term coordinated eradication efforts (13). For the United Kingdom, bovine tuberculosis (caused by infection with the pathogen Mycobacterium bovis) remains one of the highest priority pathogens affecting cattle farming industries, and one of the most costly pathogens to control to the state (4, 5). Indeed, in Northern Ireland alone, the annual cost of the eradication program in recent years has exceeded £30 million (http://www.tbhub.co.uk/tb-policy/northern-ireland/). Eradication in the UK is primarily predicated on the accurate identification of the pathogen within herds using a screening comparative tuberculin test (using bovine and avian tuberculins) to allow the slaughter of test positive animals. Subsequent repeated tuberculin testing is undertaken, and where there is evidence of within herd transmission, ancillary testing using the more sensitive interferon gamma test, and removal of additional test positive animals may be undertaken (6). The sensitivity of the screening test (single intradermal cervical comparative tuberculin (SICCT) test) is known to be moderate and can be variable in performance (710). A further complication is the presence of a known wildlife reservoir species that can harbor widespread and high levels of infection (11, 12), and which cull trials have provided evidence to suggest are involved in the risk of infection within cattle herds (13, 14). Despite these issues, modeling studies suggest that the mean reproductive number of the system (R0), which determines whether infection will be maintained in populations in perpetuity, is only slightly above 1 (1517). Therefore, only small positive changes to the eradication program could bring the trajectory below 1, which would in the long-term lead to eradication (5).

M. avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease, another significant mycobacterial pathogen of cattle, which can impact production and thus lead to economic loses to farmers (18). Unlike bTB, Johne's disease is not a statutory pathogen in Britain and Ireland, however, it is considered a priority non-regulatory pathogen of concern to farmers (19). There is evidence that co-infection of MAP can impact on the immunological response to tuberculins during statutory bTB skin testing and interferon gamma tests (2024). This is due to cross-reactivity, primarily on the avian tuberculin (which is derived from M. avium), but also impacts the reaction size on the bovine injection site (21, 25). Such cross-reactivity effects can lead to misdiagnosis (7, 26), allowing for TB positive animals to be mistakenly retained within herds. One could hypothesize that these processes may have impacts on the ability to clear infection once detected within herds (21), and subsequently impacts the duration of breakdown episodes. Unsurprisingly, the exposure to M. bovis can also have an impact on the performance of diagnostics for MAP (27). Therefore, it can be challenging to assign mechanistic inference, as the timing of infections and the testing regimes will have an impact on patterns that emerge from co-infected herds. For example, evidence suggests that tuberculin testing may impact on serum and milk ELISA testing for MAP, causing an increased risk of false positive disclosure (25, 28). Despite this, there is some epidemiological evidence that co-infection can hamper the ability to accurately identify infected animals and to clear infection from herds when the primary aim is either to eradicate bTB [e.g., (23)] or MAP [e.g., (27)].

In this study, we undertook a retrospective case-control study to investigate the hypothesis that there was an association between MAP herd status and the probability of a herd experiencing a confirmed bTB breakdown episode. We also wanted to investigate whether there were any patterns of co-infection and the longer-term bTB status of herds. We hypothesized that herds with “chronic” (i.e., long-term and recurrent problem herds) TB infection may be more likely to have evidence of co-infection than herds that did not experience chronic infections.

Methods

A retrospective, case-control study of cattle herds in Northern Ireland (United Kingdom) was undertaken (2004–2015). The primary dataset utilized during the study included serological testing data generated from surveillance activities, but primarily from part of a cattle health scheme in Northern Ireland undertaken by the Agri-food and Biosciences Institute (AFBI; www.afbi.gov.uk). These data included all blood (serological) tests undertaken during voluntary submissions, as well as those made by members of AFBI Cattle Health Scheme (CHS). Therefore, inclusion criteria restricted herds to those sampled and which voluntarily participated in MAP control schemes in Northern Ireland. It should be noted that all samples were handled and tested in the same way within AFBI irrespective of the scheme from which samples were received.

The AFBI CHS is a voluntary scheme which provides a structured approach to the control of five endemic diseases of cattle including Johne's disease. The scheme is licensed by Cattle Health Certification Standards (CHeCS) which provides the standards for licensed health schemes within the UK and Ireland. Herds enrolled in the Johne's programme have to carry out an annual herd screen consisting on blood testing all animals 2 years of age or older within the herd. In low seroprevalence herds, seropositive animals could be re-tested by fecal culture or PCR to determine the status (below for more). Introduced animals had to be blood/fecal sampled in isolation prior to joining the herd independently of the age.

Bloods were primarily sampled during concurrent statutory test sampling [i.e., TB testing or brucellosis monitoring; (29)]. All veterinarians and herd owners involved in CHS and other voluntary MAP testing with AFBI were given advice to avoid cross-reactions following intradermal tuberculin testing by either blood sampling on the first day of the test or waiting 3 months. These samples were date stamped within the dataset with the submission to the laboratory, with the statutory key performance indicator (KPI) ensuring that samples were processed within a week or less (S. Verner, pers. comm.). It should be noted that not all samples were taken at the same time as those for ancillary tests, with a risk that the time between sampling and sample submission within the field may have varied (e.g., the fecal samples, see below). However, we have no evidence to support the hypothesis that there would be any systematic bias in submission times. We assumed this random error contributed to the uncertainty in the error in our models, but potentially erring the models toward the null.

The testing was undertaken in the Diagnostic Surveillance and Investigation Branch (DSIB) at AFBI, Stormont. Blood samples were tested using the IDEXX Paratuberculosis Screening AB ELISA Test (formerly Pourquier) as per the manufacturer's recommendations. This test was accredited to ISO17025 standard. Additional information on kits used and laboratory procedures are presented in the Supplementary Material. There were also a small number of feces samples taken from animals within herds where MAP was suspected (after a positive ELISA result or clinical suspicion) or from animals added to AFBI CHS herds in the Johne's programme. These samples were tested for MAP via culturing of the pathogen in the laboratory (30), or using PCR methodologies (see Supplementary Material). These ancillary data were added to the serological dataset to improve the ability to designate a MAP herd status.

For the purposes of the present study, for each herd, a MAP status was attributed to the herd at the date of MAP testing. If any positive result for MAP was disclosed from the dataset, then the herd was attributed a “positive status” for that test event (note, this is not necessarily the case for the program, where a herd is not considered positive with a single seropositive result due to imperfect test specificity). Herds could be represented in the dataset on multiple occasions as repeated MAP tests occurred over the time series.

The Animal and Public Health and Information System (APHIS) maintained by the Department of Agriculture, Environment and Rural Affairs (DAERA, Northern Ireland) was used to gain information on herd bovine tuberculosis status, as well as deriving data on herd type, whether the herd held a milk license, which DVO the herd was assigned to, the number of animals moving into a herd per annum (buying in), and for estimating herd size.

Study Design and Modeling

The outcome variable of interest was the binary bTB status of herds at the time of the MAP testing event. The bTB status of the herd was based on breakdown episodes (31, 32). An episode was a period of time a herd was undergoing a breakdown, and was defined as the period from the first disclosing test to the last test at which the herd was released from restriction [for more on this classification see (31, 32)]. Only episodes whereby TB was post-mortem laboratory confirmed during an episode were included in this episode file (these are either multiple reactor breakdowns and/or breakdowns triggered by the presence of a TB lesion). Laboratory confirmation in Northern Ireland includes histopathology, culture, and/or strain/genotyping [spoligotyping and Variable Number Tandem Repeat genotyping; (33)]. This episode file then allowed for each herd to be classified as either “within episode,” which was our cases, and “not experiencing a breakdown” which was our control observations at the time of the MAP test event.

The primary explanatory variable was herd MAP status at the test event level based on the serological ELISA test classification. A secondary model also included ancillary data (feces) to classify herds with regards MAP. It is important to note that MAP and bTB are currently not co-managed in Northern Ireland as part of the bTB eradication scheme.

Bovine TB has a complex epidemiology, and we wanted to control for the most significant confounders of bTB risk in Northern Ireland but limiting the number of parameters in our model to ensure the model was reasonably parsimonious. Given previous experience in Northern Ireland [e.g., (6)] and published reviews of the literature [e.g., (34)], we made the a-priori decision to control for herd type [measured in two ways (Table 1)], herd size, number of animals bought in to the herd during the test year, and region (based on 10 District Veterinary Office administrative areas), while controlling for the non-independence of repeated observations within the same herd as a baseline model. We also tested whether there was any evidence of significant 2-way interactions between MAP status and (i) herd type (ii) herd size (iii) buying in, and also whether there was an interaction between herd type and (i) herd size (ii) buying in. Year was additionally tested for inclusion as a linear predictor, as a quadratic function (year2), as the natural log [log(year)], and a triennial categorical variable (dummy variable with four categories). Details of the variable are presented in Table 1. It should be noted that two herd-type variables were considered, HERD_TYPE and MILK_LICENSE. HERD_TYPE contained information about the enterprise type classification used in APHIS but also had some missing classifications, whereas MILK_LICENSE only allowed for binary classification of herds to those with/without a milk license irrespective of main enterprise type recorded within the APHIS dataset.

TABLE 1
www.frontiersin.org

Table 1. Cross-tabulations and univariable analysis of the relationship between explanatory variables and bovine TB episode cases (positive) and controls (negative) in cattle herds in Northern Ireland.

Simple unadjusted associations between variables was assessed using logistic regression. For the multivariable models, the data were fitted to a binary logit random effect regression model (random effect = herd_id). All independent variables were included in all models, with the exception of HERD_TYPE and MILK_LICENSE, which were entered into competing models and the model with the lowest Akaike's Information Criteria (AIC) value was considered the better fitting model. Where the difference in AIC was <2, models were considered equivalent (35). The random effect model was compared with a logistic regression to assess if there was significant clustering within the dataset, with a likelihood-ratio test of the intraclass correlation (ICC; ρ) being statistically different from zero. The classification ability of the model was assessed by calculating the Area under the ROC (AUC), with a model with >0.7 considered exhibiting reasonable performance.

Sensitivity Analysis

The mean number of animals sampled for Johne's Disease per herd per event was 10.02 (SD: 34.091), however there was a large variance in the distribution with the median sample of 1. The sensitivity of ELISA tests for MAP can be low [e.g., 26%; (36)], and can depend on the infectious stage of the host (37), and thus can also be highly variable in terms of performance. Recent work from the UK has suggested that moving toward a more probabilistic interpretation of MAP ELISA tests based on multiple (repeat) samples could improve the ability to assign infection status (38). Therefore, we undertook a sensitivity analysis, where separate models were fitted with the minimum number of samples per herd per event varying from ≥2 to ≥10. In each model, covariates were included [herd-type, herd-size, numbers of animals bought into the herd during the year of testing, region (DVO)], and a herd-level random effect was included to control for repeated measures within herds. The adjusted odds ratio (aOR; 95%ci), the number of observations, and the number of herds represented in the dataset for each model was reported. We also assessed models in terms of discrimination via the AUC and calibration via the Hosmer-Lemeshow test (for fixed effects only).

Despite guidelines that recommend that MAP ELISA testing should occur at least 3 months following an animal's last TB skin test (AHWNI 2018), there was a potential for some samples to have been taken within this “risk window”. To investigate this potential effect at the herd event-level, we used available animal-level data on the timing of animal's last bTB test prior to the current sampling period (the week of the MAP/bTB test, as generally MAP samples were taken on the first day of a SCITT test). In total, there were approximately 54,000 animals and 85,000 MAP tests that we could associate with a bTB skin test prior to the current skin test/MAP event. Recent research from Ireland and the UK, suggests this risk window may be < 71 days (28, 39). Any MAP test events where one or more animals had a tuberculin test within a 71 day risk window were then excluded, and the final multivariable random effect model was refit and parameter estimates compared as part of the sensitivity analysis.

Chronic Status

Milne et al. (31) and Milne et al. (32) used the episode file for bTB breakdowns to assign a status to herds based on herds experiencing problematic prolonged or recurrent breakdowns. Chronic cases were herds which experienced a prolonged (>365 days) or recurrent (1 breakdown of duration < 1year, followed by two or more bTB breakdowns within 2 years) over a time-series 2004–2015 (32). We used this classification to stratify herds recruited to the present study. We also included a classification “Non-chronic episodic herd” for herds that experienced at least one episode during the time series that were non-chronic in nature. We investigated whether there was any relationship between MAP risk and chronic herd status.

Results

Summary Statistics

There were 10,219 observations (events) where one or more MAP tests occurred constituting 553 (5.4%) TB cases and 9,666 (94.6%) controls. MAP positive results were disclosed at 3,730/10,219 (36.5%) events. Cross-tabulations at the event level between TB cases/controls and predictor variables are presented in Table 1.

Overall, 43.0% of herds (1,933/4,500) exhibited at least 1 MAP serum positive sample over the study period. 9.2% of herds (414/4,500) experienced a bTB breakdown episode during a MAP testing event. Aggregating across the whole period of the study, 138 (3.1%) herds were deemed ‘chronic' breakdown herds (i.e., experienced at least one chronic breakdown over the period), 2,237 (49.7%) were non-chronic breakdown herds (i.e., experienced at least one breakdown over the period, but this breakdown did not constitute a chronic breakdown), and 2,125 (47.2%) did not experience an episode over the period of the study (i.e., clear of infection over the whole study period).

Univariable Modeling

The unadjusted univariable associations between explanatory variables and the outcome is presented in Table 1. There was evidence of increased risk of a herd experiencing a bTB episode when disclosing with a MAP serum ELISA positive result (OR: 1.467 95%ci: 1.234–1.743), having an overall positive MAP status (serum plus additional feces data; OR: 1.476 95%ci: 1.243–1.754), being a dairy herd relative to a beef herd (OR: 1.957; 95%ci: 1.623–2.361), being in Armagh DVO relative to most other DVOs (referent in Table 1) with the exception of Newry (OR:1.116; 95%ci: 0.775–1.607) or Newtownards (OR: 1.356; 95%ci: 0.942–1.953). Increasing herd size was a risk factor (modeled as a linear predictor scaled to per 100 animals: OR: 1.308; 95%ci: 1.263–1.354), as was buying-in of cattle during the test year (linear predictor per 100 animals entering the farm: OR: 1.145; 95%ci: 1.087–1.207). There was little evidence of inter-annual variation in risk (overall categorical: p = 0.651), nor trends across years (modeled as linear: OR: 0.974; 95%ci: 0.940–1.009; quadratic: OR 0.998; 95%ci: 0.996–1.000; log(year): OR 0.891; 95%ci: 0.732–1.083).

Multivariable Model

The final multivariable model (Table 2) did not include any interaction terms as there was no evidence that these terms improved the fit of the model (i.e., aORs straddled 0, AIC increased with term addition). The model exhibited an adequate discriminatory ability with an estimated AUC of 70.29 (95%ci: 68.10–72.47). Controlling for confounders, there remained evidence that the MAP status of a herd as assessed using serological testing was positively associated with bTB episode risk at the time of testing (aOR: 1.339; 95%CI:1.085–1.652; Table 2; including ancillary feces sample data: aOR: 1.356; 95%CI:1.099–1.673).

TABLE 2
www.frontiersin.org

Table 2. Multivariable random effect binary logit model of the association between bovine TB (bTB) breakdown episode risk (outcome) and the M. avium paratuberculosis status assessed using ELISA serum testing (primary explanatory variable), while controlling for known bTB risk factors.

Sensitivity Analysis

The sensitivity analysis (Table 3; Figure 1) suggested that the effect size of this association increased as the number of samples used to inform the MAP status of herds increased. For example, a model fitted where a minimum of 2 samples was used to inform MAP herd status, suggested that bTB-MAP association had an aOR of 1.883. Increasing the minimum sample to 5 resulted in an aOR of 2.500, at ≥10 samples the aOR was 3.863. However, the overall sample (data set) size decreased with the increasing minimum number of samples required to assign a MAP status, and also the uncertainty in the parameter estimates increased (Figure 1). Despite this, the model's discriminatory ability generally improved with increasing minimum number of samples (minimum 1 sample AUC: 70.29; minimum of 10 samples AUC 75.80; Table 3), and calibration appeared not to be a large problem (based on HL on fixed effects; Hosmer-Lemeshow test: p > 0.2). It should be noted, that across all 10 models where minimum sample number was varied, the lower 95%CI aOR was always above 0, suggesting that the association was robust across the different subsamples of the population.

TABLE 3
www.frontiersin.org

Table 3. Adjusted odds ratio (aOR) for a MAP positive herd experiencing a bTB herd breakdown episode, depending on the number of samples used per herd per event to categorize the herd's MAP status.

FIGURE 1
www.frontiersin.org

Figure 1. The relationship between the number of serum samples used (2+ to 10+ per sampling event) to assign a herd-level status for Mycobacterium avium paratuberculosis (MAP) and the estimated adjusted odds ratio (aOR) for increasing bovine TB risk (left hand y-axis). The corresponding number of observations and herds included in each model declined with increasing minimum samples taken (right hand y-axis).

92.7% of MAP test events (9,477/10,219) were undertaken outside of the risk window (i.e., not < 71 days prior to samples tested). The sensitivity analysis (Table 3) suggested that excluding observations during the risk period increased the uncertainty in the parameter estimates (e.g., 2+ samples, excluding risk window obs: aOR: 1.597; 95%ci: 0.995–2.561). The magnitude of point estimate for association between MAP and bTB episode risk increased with increasing MAP samples used to assign herd MAP-status. However, this too resulted in an overall reduction in sample size. The highest AUC for any of the models during the sensitivity analysis was 78.83 for a dataset excluding observations within the risk window and having a minimum of 7 samples to inform the herd MAP status (aOR: 4.036; 95%ci: 1.451–11.226; n = 1160).

Chronic TB Herd—MAP Associations

Table 4 presents data on the cross-tabulation of the herd's bTB classification (31, 32) and their MAP status based on serology ELISA testing. There was a marginally higher proportion of serum positive MAP observations for chronic herds (41.67%), than non-chronic episodic herds (39.85%), or bTB free herds (32.22%). The association at the univariable level was moderate and above 1 for the 95%ci (OR: 1.502; 95%ci: 1.185–1.906), however the evidence in support of this association was found to be weak when controlling for repeated measures (herd random effect; OR: 1.403; 95%ci: 0.717–2.747) and when controlling for confounders (RE plus confounders: aOR 1.170; 95%ci: 0.481–2.849). There were 34 herds which had MAP tests undertaken during chronic episodes. Removing these herds diminished the difference between the proportion MAP positive from chronic herds from 41.7% to 36.2%.

TABLE 4
www.frontiersin.org

Table 4. Cross-tabulation of bTB herd classification and herd status for MAP using a serum ELISA test.

Discussion

This study identified a robust association between confirmed bovine TB breakdown episodic risk and herd MAP status, based primarily on MAP serum ELISA testing with some supplementary fecal testing. We found weaker evidence that there was an association between the long-term categorization of herds in terms of bTB risk (so called “chronic herds”). These data form part of an emerging evidence base to suggest that Johne's disease, caused by MAP infection, could be impacting on the dynamics of bTB in cattle herds, and consequently has a negative impact on bTB disease eradication (2024).

Previous research from Spain has suggested that co-infection of MAP and bTB can lead to problems in clearing infection (23). In a bullfighting herd where there was a mixed infection, Alvarez et al. (23) showed that co-infection affected the sensitivity of the skin test. Indeed, other researchers have found during experimental co-infections that the performance of both the skin test (based on Purified Protein Derivative (PPD) tuberculin injection) and the interferon gamma test [e.g., (20, 30)] can be affected. However, the direction of effects are not always consistent across reported studies [e.g., (27, 40, 41)], but it appears that overall the current consensus is that co-infection reduces the sensitivity of tests for bTB (21), but may also negatively impact on the specificity (indeed, there is still some debate as to the dominant effect). Decreasing sensitivity of the CTT could result in more bTB false-negative animals disclosed in exposed herds, lowering the probability of truly infected herds being disclosed. Also for herds that do breakdown, a lowering of sensitivity increases the likelihood of infected animals being missed during whole herd tests, resulting in either lengthening of breakdowns or increasing the probability of recrudescence (31). If MAP decreases the specificity of the CTT, then more false positives would be disclosed. This could increase bTB disclosure rates (i.e., could result in uninfected herds being classed as experiencing a breakdown), however it should not be a major factor when post-mortem/laboratory confirmation of M. bovis infection is used to assign herd breakdown status (as in the case during the present study). Specificity interference could also affect the duration of breakdowns, as herds in the UK only get released from bTB breakdowns when the herd achieves clear CTT tests. In terms of epidemiological impact, specificity interactions have financial and resource attribution costs for control programs (i.e., the cost of culling uninfected herds, and keeping trading restrictions in truly uninfected herds), while sensitivity interactions have more direct consequences for disease eradication programs as truly infected animals may not be disclosed. It should be noted though, for bTB, multiple studies have demonstrated that the CTT achieves very high specificity [e.g., 100% (95% Bayes'CrI 99%, 100%) (9)], but moderate and variable sensitivity performance [e.g., 50% (95% Bayes' CrI: 26%, 78%; (9)].

TB in a herd can also affect the performance of the MAP ELISA serum tests [e.g., (27)], potentially increasing the number of false positives. Therefore, in our study, there was the risk that herd's MAP status could be affected by bTB presence. At the animal level, leaving a minimum period between skin testing with avian tuberculin can reduce interference of the ELISA blood test (AWHNI), but if the tested animal has been exposed to M. bovis, there is the potential for cross reactivity with the MAP ELISA test. We investigated the potential impact of MAP testing within a “risk window” prior to MAP testing during the sensitivity analysis. This work supported the idea that there may have been some impact of prior tuberculin testing on the association between TB episode risk and MAP status—pushing results toward the null. However, in models with greater MAP sampling (≥3 samples) and excluding observations within the risk window, there remained a strong positive association.

The effect of prior M. bovis exposure within herds (separate from just the tuberculin use alone) could still impact on the MAP test. While we believe this to be a risk for establishing causation in our study design, there are lines of evidence that suggest this risk to be low. There is 1. an annual testing regime of all animals in Northern Ireland for bTB, where all skin test reactors are culled (29), 2. follow-up stringent use of severe interpretation and interferon gamma testing is undertaken (6), 3. in some instances additional serological tests are deployed (42), 4. animals are managed so to reduce susceptible exposure to at-risk groups, and 5. is coupled with post-mortem surveillance. With this system, it is rare that disseminated infection is found in reactor animals [~45% have lesions, and most have very few lesions; (43)], and therefore transmission within herds is low, leading to small breakdown size [median < 5 animals; (44)]. Therefore, the risk of MAP tested animals being exposed to M. bovis is low, even within breakdown herds. We utilized additional fecal data (culture and PCR) in an attempt to make our definition of MAP infection more robust, however we had limited fecal testing data available to inform our study (e.g., only 61 positive test result herds in our study).

As the impact of co-infection can impact on both diagnostics of bTB and MAP, and as MAP tests have low diagnostic sensitivity, we performed a sensitivity analysis on the number of samples used to infer a MAP status for herds. Interestingly, with increasing the numbers of samples (presumably increasing the sensitivity of the test (21), but also potential reducing the specificity of the MAP test) we found the association between the two pathogens increased, despite lowering the power to detect an effect (by reducing the overall n). As the serum ELISA MAP test has been shown to have reasonably high specificity for a screening test across studies [0.95–1.00; (37)], and as we sampled typically small numbers of animals per herd (mean 10), we do not anticipate that MAP false positive results would be the major factor influencing our inferences.

This work does add some evidence that exposure to non-tuberculosis bacteria could be having an effect in Northern Ireland (26, 30, 43). Barry et al. (30) showed using experimental animal models that co-infection could affect the gamma interferon test (but not necessarily the skin test). Using field data, our group previously found that animals with large avian tuberculin reactions (but CTT negative) had an increased risk of being TB confirmed and having greater number of lesions at slaughter, relative to non-reactor animals (26, 43). Such animals were potentially exposed to non-tuberculosis mycobacteria (including environmental exposure or MAP). Therefore, there was a risk that some co-infected animals were missed during antemortem TB testing. Coupled with the present study, there is a strong argument for further investigation into the mechanism by which MAP-bTB-coinfection is impacting on the bTB statutory eradication program.

Our models used a parsimonious base model from which we tested our co-infection hypothesis. This base model confirmed that herd size, the act of buying in of cattle, regional variation (likely representing spatial variation in underlying infection pressure) are all very important explanatory factors for bovine TB risk (34, 45). However, despite our previous work [e.g., (6, 31, 46)], there was less evidence that major enterprise type was very important in our final model (though invariably associated with the outcome). This could be related to our ability to consistently categorize herd types in Northern Ireland. In the present study we used two pieces of data available to the study, herd type and having a milk license. Both are recorded within the APHIS dataset, and both are known to have limitations. For example, where farms change their enterprise type, the recording of this change can be delayed within the dataset. Furthermore, where there are mixed enterprises, it is difficult to assign a single category. Despite this, we have previously reported differential risk at both the herd [e.g., (44)] and animal level [e.g., (43)] in Northern Ireland using similar datasets. This covariate was not our primary interest during this study, however further more refined “real-time” classification of herds is required to improve epidemiological modeling in Northern Ireland.

Limitations

Our models were based on retrospective observational datasets that were not collected for this specific purpose, and therefore the population they represent may not be fully representative of herds in the broader population in Northern Ireland. Indeed, inclusion criteria for the study was for voluntary participation in a MAP control/surveillance scheme in Northern Ireland, which may have introduced some selection bias. Furthermore, the median number of samples taken per event was 1, and therefore these animals are more likely representative of ‘suspect' cases, sampled as part of a veterinary investigations (though this is not the case for CHS herds). Indeed, therefore, there may have been some directed sampling toward the highest risk animals within the cohort. A further problem with the low number of animals tested per event is that, because of the sensitivity of the MAP test, there was the risk of misclassifying the herds. However, we addressed this with our sensitivity analysis. There was the overall problem of assessing co-infections at the herd level, though they have been undertaken previously [e.g., (47)]. Obviously, herds are a dynamic collection of animals, which can move in (buying in, birth) or out (market, abattoir, death) of herds. Therefore, finding associations at this level of organization is probably more biased toward the null than animal-level experimental analyses where there is little ambiguity as to exposure to either or both pathogens. Despite the limitations of this study, we believe there are useful associations to build further inferences about the dynamics of two difficult chronic pathogens of cattle herds.

Conclusion

We have found robust associations between bTB herd breakdown episode risk and concurrent MAP infections at the herd level. This association was positive, indicating that breakdown herds exhibited a higher risk of co-infection than sampled herds which were not experiencing a breakdown. Furthermore, there was a marginally higher proportion of chronically infected herds with evidence of MAP infection than bTB free herds over the study period. During the study we investigated whether the sampling frequency of herds to determine MAP status, and the time-period between a previous tuberculin test and MAP test, had an impact on these associations. Increasing sampling frequency increased the magnitude of the point estimate associations and the uncertainty around these estimates, while decreasing N. The sampling of animals within a “risk window” where cross-reactivity may occur after a TB tuberculin test did appear to impact the association by increasing uncertainty and lowering the magnitude of the association. The exact mechanism for our findings cannot be readily determined with these data, however we hypothesize given results from elsewhere that MAP may be masking bTB infection in herds experiencing TB breakdowns. This masking effect may impact on the ability to identify exposed cattle, cull them, and thus clear infection. Further research is required to better understand how these patterns emerge, and whether (or which) co-management tools would yield optimal control benefits.

Author Contributions

AB conceived and designed the study, project managed, undertook data analysis, and interpretation of data for the paper. AB drafted the initial manuscript and revised it critically. JG managed all research databases contributing to the work, and was involved in data merging, validation, and cleaning, and contributed to data interpretation. GM was involved in generating the chronic herds data set and episode file, contributed to the data interpretation, and critically revised the manuscript. MG-G led and coordinated the AFBI cattle health scheme (2010–2018) which contributed data, managed the diagnostic laboratories undertaking MAP testing, contributed to data interpretation and study design, and critically revised the manuscript. SS was involved in the initial concept for the work, managed a Johne's control program (AHWNI) in Northern Ireland which contributed data, contributed to data interpretation and study design, and critically revised the manuscript. All authors read and made comments and/or edits to the final manuscript.

Funding

Funding was provided as part of the Evidence and Innovation Strategy from the Department of Agriculture, Environment, and Rural Affairs (DAERA; Grant numbers: 15/03/08 Tuberculosis in chronic herds in Northern Ireland; PI: AB; 15/03/10 Concurrent endemic diseases in Northern Ireland; PI: AB). Open access publication fees were institutionally funded (AFBI).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We would like to acknowledge the laboratory staff at DSIB, AFBI. We also wish to acknowledge colleagues within Epidemiology, Molecular biology, and Immunology (EMbI; https://twitter.com/AFBI_EMbI) research cluster within Bacteriology branch AFBI for stimulating discussions about this topic and help in the thinking around this confection problem. Data were provided by DAERA (data controllers) under a data sharing agreement to AFBI.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2019.00030/full#supplementary-material

References

1. Benedictus G, Verhoeff J, Schukken YH, Hesselink JW. Dutch paratuberculosis programme history, principles and development. Vet Microbiol. (2000) 77:399–413. doi: 10.1016/S0378-1135(00)00325-4

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Olmstead AL, Rhode PW. An impossible undertaking: the eradication of bovine tuberculosis in the United States. J Econ History (2004) 64:734–72. doi: 10.1017/S0022050704002955

CrossRef Full Text | Google Scholar

3. Gordejo FR, Vermeersch JP. Towards eradication of bovine tuberculosis in the European Union. Vet Microbiol. (2006) 112:101–9. doi: 10.1016/j.vetmic.2005.11.034

CrossRef Full Text | Google Scholar

4. Reynolds D. A review of tuberculosis science and policy in Great Britain. Vet Microbiol. (2006) 112:119–26. doi: 10.1016/j.vetmic.2005.11.042

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Allen A, Skuce R, Byrne A. Bovine tuberculosis in Britain and Ireland–A Perfect Storm? The confluence of potential ecological and epidemiological impediments to controlling a chronic infectious disease Front. Veterin. Sci. (2018) 5:109. doi: 10.3389/fvets.2018.00109

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Lahuerta-Marin A, Gallagher M, McBride S, Skuce R, Menzies F, McNair J, et al. Should they stay, or should they go? Relative future risk of bovine tuberculosis for interferon-gamma test-positive cattle left on farms. Vet Res. (2015) 46:90. doi: 10.1186/s13567-015-0242-8

CrossRef Full Text | Google Scholar

7. De la Rua-Domenech R, Goodchild AT, Vordermeier HM, Hewinson RG, Christiansen KH, Clifton-Hadley RS. Ante mortem diagnosis of tuberculosis in cattle: a review of the tuberculin tests, γ-interferon assay and other ancillary diagnostic techniques. Res Vet Sci. (2006) 81:190–210. doi: 10.1016/j.rvsc.2005.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Clegg TA, Duignan A, Whelan C, Gormley E, Good M, Clarke J, et al. Using latent class analysis to estimate the test characteristics of the γ-interferon test, the single intradermal comparative tuberculin test and a multiplex immunoassay under Irish conditions. Vet Microbiol. (2011) 151:68–76. doi: 10.1016/j.vetmic.2011.02.027

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Nuñez-Garcia JSH, Downs JE, Parry DA, Abernethy JM, Broughan AR, Cameron AJ, Cook et al. Meta-analyses of the sensitivity and specificity of ante-mortem and post-mortem diagnostic tests for bovine tuberculosis in the, UK and Ireland. Prev Vet Med. (2018) 153:94. doi: 10.1016/j.prevetmed.2017.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lahuerta-Marin A, Milne MG, McNair J, Skuce R, McBride S, Menzies F, et al. Bayesian latent class estimation of sensitivity and specificity parameters of diagnostic tests for bovine tuberculosis in chronic herds-Northern Ireland. Vet J. (2018) 238:15–21. doi: 10.1016/j.tvjl.2018.04.019

CrossRef Full Text | Google Scholar

11. Murphy D, Gormley E, Costello E, O'Meara D, Corner LAL. The prevalence and distribution of Mycobacterium bovis infection in European badgers (Meles meles) as determined by enhanced post mortem examination and bacteriological culture. Res Vet Sci. (2010) 88:1–5. doi: 10.1016/j.rvsc.2009.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Byrne AW, Kenny K, Fogarty U, O'Keeffe JJ, More SJ, McGrath G, et al. Spatial and temporal analyses of metrics of tuberculosis infection in badgers (Meles meles) from the Republic of Ireland: trends in apparent prevalence. Prev Vet Med. (2015) 122:345–54. doi: 10.1016/j.prevetmed.2015.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Griffin JM, Williams DH, Kelly GE, Clegg TA, O'boyle I, Collins JD, et al. The impact of badger removal on the control of tuberculosis in cattle herds in Ireland. Prev. Vet. Med. (2005) 67:237–266. doi: 10.1016/j.prevetmed.2004.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Donnelly CA, Nouvellet P. The contribution of badgers to confirmed tuberculosis in cattle in high-incidence areas in England. PLoS Curr. (2013) 5. doi: 10.1371/currents.outbreaks.097a904d3f3619db2fe78d24bc776098

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Cox DR, Donnelly CA, Bourne FJ, Gettinby G, McInerney JP, Morrison WI, et al. Simple model for tuberculosis in cattle and badgers. Proc Nat Acad Sci USA. (2005) 102:17588–93. doi: 10.1073/pnas.0509003102

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Aznar I, Frankena K, Byrne AW, More SJ, De Jong MCM. Infection dynamics and effective control options of tuberculosis in cattle and badgers. In: 6th International M. bovis Conference. (Cardiff) (2014).

Google Scholar

17. Brooks-Pollock E, Wood JL. Eliminating bovine tuberculosis in cattle and badgers: insight from a dynamic model. Proc. R. Soc. B. (2015) 282:20150374. doi: 10.1098/rspb.2015.0374

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ott SL, Wells SJ, Wagner BA. Herd-level economic losses associated with Johne's disease on US dairy operations. Prev Vet Med. (1999) 40:179–92. doi: 10.1016/S0167-5877(99)00037-9

PubMed Abstract | CrossRef Full Text | Google Scholar

19. More SJ, McKenzie K, O'Flaherty J, Doherty ML, Cromie AR, Magan MJ. Setting priorities for non-regulatory animal health in Ireland: results from an expert Policy Delphi study and a farmer priority identification survey. Prev Vet Med. (2010) 95:198–207.

PubMed Abstract | Google Scholar

20. Dunn JR, Kaneene JB, Grooms DL, Bolin SR, Bolin CA, Bruning-Fann CS. Effects of positive results for Mycobacterium avium subsp paratuberculosis as determined by microbial culture of feces or antibody ELISA on results of caudal fold tuberculin test and interferon-γ assay for tuberculosis in cattle. J Am Vet Med Assoc. (2005) 226;429–435. doi: 10.2460/javma.2005.226.429

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Hope JC, Thom ML, Villarreal-Ramos B, Vordermeier H M, Hewinson RG, Howard CJ. Exposure to Mycobacterium avium induces low-level protection from Mycobacterium bovis infection but compromises diagnosis of disease in cattle. Clin Exp Immunol. (2005) 141:432–439. doi: 10.1111/j.1365-2249.2005.02882.x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Aranaz A, Bezos J, Álvarez J, Romero B, Lozano F, Paramio JL., et al. Assessment of diagnostic tools for eradication of bovine tuberculosis in cattle co-infected with Mycobacterium bovis and M. avium subsp. paratuberculosis. Vet. Res. (2006) 37:593–606.

PubMed Abstract | Google Scholar

23. Álvarez J, De Juan L, Bezos J, Romero B, Sáez JL, Marqués S, et al. Effect of paratuberculosis on the diagnosis of bovine tuberculosis in a cattle herd with a mixed infection using interferon-gamma detection assay. Vet. Microbiol. (2009) 135:389–393.

PubMed Abstract | Google Scholar

24. Roupie V, Alonso-Velasco E, Van Der Heyden S, Holbert S, Duytschaever L, Berthon P, et al. Evaluation of mycobacteria-specific gamma interferon and antibody responses before and after a single intradermal skin test in cattle naturally exposed to M. avium subsp paratuberculosis and experimentally infected with M. bovis. Vet Immunol Immunopathol. (2018) 196:35–47. doi: 10.1016/j.vetimm.2017.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kennedy AE, Byrne N, O'Mahony J, Sayers RG. Investigations and implications of associations between mycobacterial purified protein derivative hypersensitivity and MAP-antibody ELISA in Irish dairy cows. Res Vet Sci. (2017) 115:13–6. doi: 10.1016/j.rvsc.2017.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Byrne AW, Graham J, Brown C, Donaghy A, Guelbenzu-Gonzalo M, McNair J, et al. Modelling the variation in skin-test tuberculin reactions, post-mortem lesion counts and case pathology in tuberculosis-exposed cattle: Effects of animal characteristics, histories and co-infection. Transbound Emerg Dis. (2018) 65:844–58. doi: 10.1111/tbed.12814

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lilenbaum W, Marassi CD, Varges R, Medeiros L, Oelemann WMR, Fonseca LS. Occurrence of false-positive results in three Paratuberculosis-ELISAs performed in a tuberculous herd. Vet Res Comm. (2009) 33:693–9. doi: 10.1007/s11259-009-9218-3

PubMed Abstract | CrossRef Full Text | Google Scholar

28. May H, Orpin P, Pearse H, Jones P, Clough H. The Influence of Tuberculosis Testing in Dairy Cattle on Milk ELISA Tests for Johne's Disease (2008). Available online at: https://www.researchgate.net/publication/237445560_The_influence_of_tuberculosis_testing_in_dairy_cattle_on_milk_ELISA_tests_for_Johne's_disease

29. Abernethy DA, Denny GO, Menzies FD, McGuckian P, Honhold N, Roberts AR. The Northern Ireland programme for the control and eradication of Mycobacterium bovis. Vet Microbiol. (2006) 112:231–7. doi: 10.1016/j.vetmic.2005.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Barry C, Corbett D, Bakker D, Andersen P, McNair J, Strain S. The effect of Mycobacterium avium complex infections on routine Mycobacterium bovis diagnostic tests. Vet. Med. Int. (2011) 2011:145092 doi: 10.4061/2011/145092

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Milne GM, Graham J, Allen A, Lahuerta-Marin A, McCormick C, Presho E, et al. Spatiotemporal analysis of prolonged and recurrent bovine tuberculosis breakdowns in Northern Irish cattle herds reveals a new infection hotspot. Spatial Spatio-Temporal Epidemiol. (2019) 28:33–42. doi: 10.1016/j.sste.2018.11.002

CrossRef Full Text | Google Scholar

32. Milne G, Graham J, Allen A, Lahuerta-Marin A, McCormick C, Presho E, et al. Herd characteristics, wildlife risk and bacterial strain genotypes in persistent breakdowns of bovine tuberculosis in Northern Irish cattle herds. In: Society for Veterinary Epidemiology and Preventive Medicine–Proceedings. Tallinn, Estonia. (2018). p. 21–3.

Google Scholar

33. Skuce RA, Byrne AW, Lahuerta-Marin A, Allen A. Chapter 5. Mycobacterium bovis molecular typing and surveillance. In: Chambers M, Gordon S, Olea-Popelka F, Barrow P, editors. Bovine Tuberculosis. Wallingford, VT: CABI (2018). p. 58–79. doi: 10.1079/9781786391520.0058

CrossRef Full Text | Google Scholar

34. Skuce RA, Allen AR, McDowell SW. Herd-level risk factors for bovine tuberculosis: a literature review. Vet Med Int. (2012) 2102:621210. doi: 10.1155/2012/621210.

CrossRef Full Text | Google Scholar

35. Dohoo IR, Martin SW, Stryhn H. Veterinary Epidemiologic Research. Charlottetown, PE: AVC Inc (2009).

Google Scholar

36. Alinovi CA, Ward MP, Lin TL, Moore GE, Wu CC. Real-time PCR, compared to liquid and solid culture media and ELISA, for the detection of Mycobacterium avium ssp. paratuberculosis Vet Microbiol. (2009) 136:177–9. doi: 10.1016/j.vetmic.2008.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Nielsen SS, Toft N. Ante mortem diagnosis of paratuberculosis: a review of accuracies of ELISA, interferon-γ assay and faecal culture techniques. Vet Microbiol. (2008) 129:217–35. doi: 10.1016/j.vetmic.2007.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Meyer A, Bond K, Van Winden S, Green M, Guitian J. A probabilistic approach to the interpretation of milk antibody results for diagnosis of Johne's disease in dairy cattle. Prev Vet Med. (2018) 150:30–7. doi: 10.1016/j.prevetmed.2017.11.016

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Kennedy AE, Da Silva AT, Byrne N, Govender R, MacSharry J, O'Mahony J, et al. The single intradermal cervical comparative test interferes with Johne's disease ELISA diagnostics. Front Immunol. (2014) 5:564. doi: 10.3389/fimmu.2014.00564

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Seva J, Sanes JM, Ramis G, Mas A, Quereda JJ, Villarreal-Ramos B, et al. Evaluation of the single cervical skin test and interferon gamma responses to detect Mycobacterium bovis infected cattle in a herd co-infected with Mycobacterium avium subsp. paratuberculosis. Vet. Microbiol. (2014) 171:139–46. doi: 10.1016/j.vetmic.2014.03.035

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Raffo E, Steuer P, Monti G, Salgado M. Effect of Mycobacterium avium subsp. paratuberculosis (MAP) infection on the diagnostic accuracy for Mycobacterium bovis (M bovis) infection under field conditions in cattle belonging to low M bovis prevalence herds Trop Anim Health Prod. (2017) 49:771–5. doi: 10.1007/s11250-017-1259-y

PubMed Abstract | CrossRef Full Text | Google Scholar

42. McCallan L, Brooks C, Couzens C, Young F, McNair J, Byrne AW. Assessment of serological tests for diagnosis of bovine tuberculosis. Vet Record (2017) 181:90. doi: 10.1136/vr.104272

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Byrne AW, Guelbenzu-Gonzalo M, Strain SAJ, McBride S, Graham J, Lahuerta-Marin A, et al. Assessment of concurrent infection with bovine viral diarrhoea virus (BVDV) and Mycobacterium bovis: a herd-level risk factor analysis from Northern Ireland. Prev Vet Med. (2017) 141:38–47. doi: 10.1016/j.prevetmed.2017.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Byrne AW, Graham J, Brown C, Donaghy A, Guelbenzu-Gonzalo M, McNair J, et al. Bovine tuberculosis visible lesions in cattle culled during herd breakdowns: the effects of individual characteristics, trade movement and co-infection. BMC Vet. Res. (2017) 13:400. doi: 10.1186/s12917-017-1321-z

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Broughan JM, Judge J, Ely E, Delahay RJ, Wilson G, Clifton-Hadley RS, et al. A review of risk factors for bovine tuberculosis infection in cattle in the UK and Ireland. Epidemiol Infect. (2016) 144:2899–926. doi: 10.1017/S095026881600131X

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Lahuerta-Marin A, McNair J, Skuce R, McBride S, Allen M, Strain SA, et al. Risk factors for failure to detect bovine tuberculosis in cattle from infected herds across Northern Ireland (2004–2010). Res Vet Sci. (2016) 107:233–9. doi: 10.1016/j.rvsc.2016.06.014

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Vanleeuwen JA, Haddad JP, Dohoo IR, Keefe GP, Tiwari A, Scott HM. Risk factors associated with Neospora caninum seropositivity in randomly sampled Canadian dairy cows and herds. Prev Vet Med. (2010) 93:129–138. doi: 10.1016/j.prevetmed.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: bovine TB, Johne's disease, veterinary epidemiology, co-infection, infectious disease control, mycobacteria

Citation: Byrne AW, Graham J, Milne G, Guelbenzu-Gonzalo M and Strain S (2019) Is There a Relationship Between Bovine Tuberculosis (bTB) Herd Breakdown Risk and Mycobacterium avium subsp. paratuberculosis Status? An Investigation in bTB Chronically and Non-chronically Infected Herds. Front. Vet. Sci. 6:30. doi: 10.3389/fvets.2019.00030

Received: 30 October 2018; Accepted: 23 January 2019;
Published: 14 February 2019.

Edited by:

Javier Bezos, Complutense University of Madrid, Spain

Reviewed by:

Francisco Ruiz-Fons, Spanish National Research Council (CSIC), Spain
Ana Cláudia Coelho, University of Trás-os-Montes and Alto Douro, Portugal

Copyright © 2019 Byrne, Graham, Milne, Guelbenzu-Gonzalo and Strain. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Andrew W. Byrne, ecologicalepidemiology@gmail.com; andrew.byrne@afbini.gov.uk