Skip to main content

ORIGINAL RESEARCH article

Front. Vet. Sci., 13 August 2019
Sec. Veterinary Epidemiology and Economics

Modeling the Accuracy of Two in-vitro Bovine Tuberculosis Tests Using a Bayesian Approach

  • 1Department of Veterinary Population Medicine, University of Minnesota, Saint Paul, MN, United States
  • 2Facultad de Veterinaria, Universidad de la Republica, Montevideo, Uruguay
  • 3División Laboratorios Veterinarios “Miguel C. Rubino”, Ministerio de Ganadería, Agricultura y Pesca, Montevideo, Uruguay
  • 4VISAVET Health Surveillance Centre, Universidad Complutense, Madrid, Spain
  • 5Departamento de Sanidad Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain

Accuracy of new or alternative diagnostic tests is typically estimated in relation to a well-standardized reference test referred to as a gold standard. However, for bovine tuberculosis (bTB), a chronic disease of cattle, affecting animal and public health, no reliable gold standard is available. In this context, latent-class models implemented using a Bayesian approach can help to assess the accuracy of diagnostic tests incorporating previous knowledge on test performance and disease prevalence. In Uruguay, bTB-prevalence has increased in the past decades partially because of the limited accuracy of the diagnostic strategy in place, based on intradermal testing (caudal fold test, CFT, for screening and comparative cervical test, CCT, for confirmation) and slaughter of reactors. Here, we evaluated the performance of two alternative bTB-diagnostic tools, the interferon-gamma assay, IGRA, and the enzyme-linked immunosorbent assay (ELISA), which had never been used in Uruguay in the absence of a gold standard. In order to do so animals from two heavily infected dairy herds and tested with CFT-CCT were also analyzed with the IGRA using two antigens (study 1) and the ELISA (study 2). The accuracy of the IGRA and ELISA was assessed fitting two latent-class models: a two test-one population model (LCA-a) based on the analysis of CFT/CFT-CCT test results and one in-vitro test (IGRA/ELISA), and a one test-one population model (LCA-b) using the IGRA or ELISA information in which the prevalence was modeled using information from the skin tests. Posterior estimates for model LCA-a suggested that IGRA was as sensitive (75–78%) as the CFT and more sensitive than the serial use of CFT-CCT. Its specificity (90–96%) was superior to the one for the CFT and equivalent to the use of CFT-CCT. Estimates from LCA-b models consistently yielded lower posterior Se estimates for the IGRA but similar results for its Sp. Estimates for the Se (52% 95%PPI:44.41-71.28) and the Sp (92% 95%PPI:78.63–98.76) of the ELISA were however similar regardless of the model used. These results suggest that the incorporation of IGRA for detection of bTB in highly infected herds could be a useful tool to improve the sensitivity of the bTB-control in Uruguay.

Introduction

Accuracy of diagnostic tests has been traditionally estimated by comparing the test results with those of a reference test, sometimes referred to as the gold standard, which unequivocally indicates the true status of an individual (infected/not infected). In the absence of such a reference test, latent class analyses based on Bayesian methods provide an alternative strategy for evaluation of diagnostic tests when the true status of the individual is unknown. The use of this approach in the context of veterinary medicine has been described elsewhere (1). Briefly, the use of latent class analyses based on Bayesian methods involves the combination of previous knowledge on test performance (when available) with the evidence provided by newly collected data to obtain a posterior estimate on test performance and disease prevalence, often achieved through Monte Carlo simulations using Gibbs sampling (2). The prior knowledge on test performance is typically obtained through the review of the scientific literature and/or the elicitation of expert opinion (3). Methodologies to elicit expert opinion have been described elsewhere (3).

Use of latent class models in veterinary epidemiology has increased in the past decades, particularly for the assessment of diagnostic tests for chronic and complex diseases for which gold standard tests are not available, such as bovine tuberculosis (bTB) (47).

Bovine tuberculosis, mainly caused by infection with Mycobacterium bovis (M. bovis), is an important chronic disease of cattle that causes a substantial impact on animal and public health, and that imposes a significant economic burden associated to its control and international trade restrictions (8, 9).

Control programs worldwide are based on test and removal of positive animals or, in some cases, complete herds (10). In Uruguay, the bTB-national program involves serial intradermal testing (caudal fold test –CFT- followed by the comparative cervical test –CCT- for confirmation) of all dairy herds annually for the detection of infected animals and its posterior removal (11, 12). In the past decade, the number of bTB-positive dairy herds detected every year, the within-herd prevalence in infected farms, and the time from outbreak detection to control has increased in Uruguay despite measures implemented as part of the national bTB control program (12, 13). The evolution of the dairy industry in the country, characterized by an increase in herd sizes and production intensification, has been associated with the limited success of bTB-control in recent years (12, 14). Additionally, insufficient sensitivity of bTB diagnostic tests may also contribute to the persistence of potentially infectious individuals in the herd that can further spread the disease within and between herds (15).

In Europe the use of the interferon-gamma release assay (IGRA) in parallel with the skin test has been incorporated in many eradication programs to maximize diagnostic sensitivity (Council Directive 64/432/EEC, 1964) (16, 17). Other tests based in the detection of specific antibodies (such as the enzyme-linked immunosorbent assay, ELISA) have been developed and proven useful for detection of specific subpopulations of M. bovis-infected animals that may not react to the skin test, although their field use has been mostly limited so far to experimental purposes (1823). Differences in the performance of bTB diagnostic tests can be related with local factors related with the personnel conducting the tests (experience) or with the cattle population (frequency of testing, presence of other diseases compromising the immune response, breed, among others). Characterization of the performance of alternative diagnostic tools (IGRA and ELISA) previously never implemented in Uruguay may help to design strategies for the improvement of the diagnostic sensitivity in high bTB-prevalence infected dairy herds, currently a priority for the control and eradication of bTB in the country.

Here, we aimed to estimate the accuracy of two commercial assays for in-vitro diagnosis of bTB that had never been used in Uruguay, namely an IGRA (using two alternative antigens –referred to as study1-) and an antibody-based ELISA (referred to as study 2), fitting two different latent-class models in a Bayesian framework. Results from this research will help to quantify the potential impact that alternative diagnostic strategies may have in improving the effectiveness of the bTB-control program in Uruguay.

Methods

We followed the STARD-BLCM guidelines to describe the materials and methods in our study (24).

Study Design and Source Population

Two cross-sectional diagnostic accuracy studies (referred to as study 1 and study 2 –Figure 1) using a “single-gate” diagnostic design were performed to evaluate the performance of the in-vitro bTB-diagnostic assays (25).

FIGURE 1
www.frontiersin.org

Figure 1. Schematic diagram showing the study design, with the diagnostic tests used for study 1 and study 2, and the Bayesian latent-class fitted models LCA-a and LCA-b.

Sampling for both studies was carried in 2016, and included 121 and 279 Holstein cows for studies 1 and 2, respectively. All animals were selected from two commercial dairy herds belonging to the same company (with similar management practices and that frequently and systematically mingle their animals) located in the Department of Florida. Both herds were bTB positive since 2013. The two herds were subjected to the intradermal test as regulated by the national bTB-control program in Uruguay for dairies based on the status of the herd (26). In addition, blood and serum samples were drawn from the selected animals (Figure 1). For logistics reasons blood samples from 158 animals were not collected reducing the sample size for study 1 in comparison to study 2. To avoid selection bias, the proportion of CFT-positive results in animals with missing blood samples in comparison to the ones that have blood samples for study 1 were assessed using a Pearson's chi-square test.

Sampling and Diagnostic Assays

All dairy >12 month animals were tested using the CFT as a screening test, involving the intradermal inoculation of a purified protein derivate from M. bovis (PPDb) in the caudal area. Animals with an increase in skin thickness and/or presence of in-situ clinical signs of inflammation 72 h post inoculation were considered reactors and subjected to the CCT for confirmation within the following seven days. In this test, two PPD inoculations from M. bovis (PPDb) and M. avium (PPDa) are performed in the cervical area. When the difference in skinfold thickness in the PPDb inoculation site was ≥4 mm than the PPDa animals were considered infected and culled.

Blood samples from the coccygeal vein of cows enrolled in studies 1 and 2 were collected in tubes with (IGRA) or without (antibody ELISA) heparin, respectively, after the results of the serial CCT test were assessed (if applicable). Blood samples were maintained in a container at monitored environmental temperatures (22 ± 3°C) until arrival to the official veterinary diagnostic laboratory (Miguel C. Rubino) within the first 8 h post extraction to perform the IGRA (Prionics, Schlieren-Zurich, Switzerland). Serum samples were refrigerated until the performance of the ELISA (IDEXX Laboratories, Westbrook, ME, respectively).

In study 1, blood samples were stimulated with specific antigens as described elsewhere (27). All samples were divided into five aliquots and incubated for 18 h with pokeweed mitogen, PBS (blank), PPDa, PPDb and an antigenic cocktail formed by the early secretory antigenic target-6 (ESAT-6) and the culture filtrate protein 10 (CFP-10), two highly specific M. bovis antigenic proteins (28). Samples were then centrifuged and the supernatant was analyzed using the Bovigam 2.G (Prionics, Schlieren-Zurich, Switzerland) according to the manufacturer's recommendations. Two criteria based on different sets of antigens were applied to classify animals as positive; for criteria A (IGRAb) animals were considered positive if the optical density (OD) obtained after stimulation with PPDb (ODPPDb) minus the OD of the aliquot stimulated with PBS (ODPBS) was ≥0.1 and ODPPDb-ODPPDa ≥0.1; in the case of criteria B (IGRAc), animals were classified as positive when ODcocktail-ODPBS was ≥0.1. For study 2 a commercial ELISA (IDEXX Laboratories, Westbrook, ME) was used to detect MPB83 and MPB70 bTB specific antibodies as described elsewhere (19). Animals with an S/P ratio ≥ 0.3 were considered positive and negative if else as recommended by the manufacturer.

Statistical Models

Latent-class models were used to estimate diagnostic test accuracy (sensitivity –Se-, and specificity -Sp-) of the IGRA using the different antigens (IGRAb and IGRAc) and the ELISA in the absence of a gold standard assay (1, 29). Samples collected were assumed to originate from a single population given they were drawn from herds belonging to the same company with similar animal health status regarding bTB and similar production management standards.

For each study (1 and 2) two different models were used alternatively: a two dependent tests-one population model (LCA-a) using the results from the skin test (CFT or CFT-CCT) and one of the in-vitro tests (IGRA or ELISA), and a one test-one population model (LCA-b) analyzing the results of the in-vitro tests separately (Figure 1).

Conditional correlation coefficients for the Se (rhoD) and Sp (rhoDc) were included in the LCA-a models as described elsewhere (29). We assumed results from the tests were conditionally dependent because results from diagnostic tests targeting a similar biological phenomenon, such as the intradermal tests and the IGRA (30), are likely dependent (29, 31). Similarly, and although the ELISA is based on the detection of the humoral immune response in the infected animals, there is a relationship between the initial predominant cellular-mediated immunity and the posterior humoral immunity observed as disease progresses in the animal (32), so results from the skin test and the ELISA were also assumed to be conditionally dependent.

Beta prior distributions for the Se and Sp of the CFT, CFT-CCT, IGRAb, IGRAc, and ELISA were chosen according to previous reports (Table 1, Supplementary Table 1). Distributions were fitted using Beta buster version 1.0 (downloadable at https://cadms.vetmed.ucdavis.edu/diagnostic/software). More informative distributions were used for the Se and Sp of the CFT-CCT due to the availability of Uruguay-specific information (40) compared with those used for the in-vitro assays, since most references for those originated from other countries with a different experience in the use of these techniques (Table 1).

TABLE 1
www.frontiersin.org

Table 1. Prior estimates (Mode and 5th percentiles) for sensitivity, specificity of the intradermal tests (CFT, CFT-CCT) and in-vitro (IGRAb, IGRAc, and ELISA) bTB tests, and prevalence for the two models implemented.

For the LCA-a (two-dependent-test) models, prevalence priors were formulated from expert opinion following procedures described elsewhere (3). For the LCA-b (one-test) models, prior distributions for prevalence were formulated using the results from the CFT-CCT as described previously (48). Briefly, we simulated the true prevalence distribution using the Rogan-Gladen estimation method to correct for the imperfect Se and Sp of the CFT-CCT (assumed to follow beta distributions as mentioned before) (Table 1) through 5,000 iterations in an Excel spreadsheet (Microsoft Office Professional Edition, 2016) using @Risk software version 7.0.0 (Palisade Corporation 2015). The outputs from the simulations were used to fit a beta distribution that was used as the prevalence prior for LCA-b models.

Three Markov chain Monte Carlo runs were implemented per model to visually assess convergence (also tested using the Gelman-RubinR statistic) (49). Models were run for 7,500 iterations for computing posterior estimates after an initial burn-in of 2,500 samples. To eliminate potential autocorrelation we applied thinning and selected one every 10 consecutive samples. Latent-class models were fitted using OpenBUGS 3.2.2 (50) via the R2OpenBUGS package (51) from the R 3.2.4 software. The influence of the selected priors on the posteriors distributions was evaluated by comparing the initial models with a model fitted using non-informative uniform (0,1) distributions for each parameter under evaluation. The possible independence between the results of the two tests being assessed was also evaluated by fitting models that did not include correlation terms. Model fit was assessed using the deviance information criterion (DIC), and the model selection (LCA-a or LCA-b) was based on lower DIC (52) and narrower posterior credibility intervals.

Results

Cross-tabulated dichotomous results for the combination of the intradermal tests (CFT or CFT-CCT), and the in-vitro assays (IGRAb, IGRAc, or ELISA) are presented in Table 2. Animals with missing IGRAs diagnostic results followed similar proportion of CFT results than those used for the analyses (Pearson's chi-square 2.91, P > 0.05).

TABLE 2
www.frontiersin.org

Table 2. Cross-tabulated dichotomous diagnostic results for intradermal test (CFT, CFT-CCT) and in-vitro (IGRAb, IGRAc, ELISA) bTB- diagnostic tests.

The estimated posterior estimates for the Se and Sp of the diagnostic tests and the prevalence in Study 1 and 2 are shown in Table 3.

TABLE 3
www.frontiersin.org

Table 3. Posterior estimates (median and 95% posterior probability interval) for CFT, CFT-CCT and in-vitro assays (IGRAb, IGRAc, ELISA) sensitivities, specificities, prevalence, and, when applicable, correlation terms (rhoD, rhoDc) distributions obtained for study 1 (121 animals) and study 2 (279 animals), applying the model “a,” or the model “b” in chronic naturally infected dairy herds in Uruguay.

Study 1

Median posterior estimates for the prevalence, Se and Sp of the intradermal tests (CFT, and CFT-CCT) using the LCA-a model were similar regardless the antigen used in the IGRA (IGRAb or IGRAc) (Table 3). The median posterior IGRAb Sp estimates were slightly lower than those obtained for the IGRAc, whit higher median Sp values for the models integrating CFT-CCT as second test as well, but with the overlapping of the PPIs (Table 3).

LCA-b models consistently yielded lower Se values for both IGRAs and higher prevalence estimates compared with LCA-a models, but with similar Sp posterior estimates.

Study 2

The LCA-a model yielded higher posterior estimates for the prevalence and Se of the intradermal tests, and a markedly lower Sp posterior values for CFT compared to those observed in study 1 using the same model. ELISA Se and Sp estimates obtained using the two models (LCA-a and b) were consistent.

Conditional correlation between intradermal and in-vitro test results in infected (rhoD) and non-infected animals (rhoDc) was low, with a 95% Posterior Probability Interval (95%PPI) including 0 in all LCA-models for study 1 (Table 3). However, no significant improvement was observed in the DIC when test independence was assumed for models using IGRAs and CFT (study 1: 19.4 vs. 19.4, 19.4 vs. 19) or IGRAs and CFT-CCT (study 1: 18.1 vs. 17.5, 17.7 vs. 19.7) respectively. Interestingly, the LCA-a model from study 2 showed the highest median correlation terms for infected animals (rhoD = 11.7 and 17.05), showing a poorer fit of the model when independent-tests models were assessed (ELISA and CFT DIC:19.4 vs. 24.7, ELISA and CFT-CCT DIC:24.4 vs. 31.9), although 95% PPI included 0.

The sensitivity analysis revealed that results obtained using LCA-a models for study 1 were not affected (changes <10.5%) by the use of weakly informative priors (Supplementary Table 2). However, various parameters were severely affected (changes>10.5% when weakly informative priors were used) by the choice of priors in the remaining models/studies. Results were most affected when LCA-a models were applied in study 2. The use of uniform distributions for the Sp of the in-vitro assays in both studies resulted in 17.8 to 38.2% decreased posterior median Sp values. Similarly, use of uniform priors for the prevalence resulted in a >20% reduction in posterior estimates of study 1 using the LCA-b model (76.5 to 54.3 and 65.4 to 52.2), and an increase in the Se estimates for IGRAb and IGRAc.

All models reached convergence as indicated by the visual inspection of the Markov chains and the Gelman-Rubin∧R statistic (<1.002) for all parameters.

Discussion

Due to the increasing number of bTB- infected herds in Uruguay (Animal Health Bureau, Uruguay -DSA MGAP-), the need for early and accurate detection, isolation and removal of infected animals from a herd is crucial when whole herd-culling is not an economically or socially sustainable option. Here, we aimed to assess the performance of bTB-in-vitro assays under field conditions for the first time in Uruguay with the ultimate goal of improving current bTB diagnostic strategies for chronic and high prevalence infected dairy herds.

In order to estimate the performance of the in-vitro assays evaluated here we used LCA, a suitable analytical approach when no reliable gold standard is available (1, 53, 54), as it is the case for bTB (5, 6, 55). We fitted two different latent-class models using prevalence priors based on expert opinion or diagnostic test results in order to evaluate the potential impact of a given methodological approach. Based on DIC models with three (Se, Sp, Prev) parameters were preferred above those with seven (Se1, Se2, rhoD, Sp1, Sp2, rhoDc, Prev). Correlation between test results were very low in all models/test pairs, what had been already described for the IGRA and single skin test (5, 6) such result is expected because the diagnostic tests evaluated here have high Sp (1). However, the comparatively higher correlation between the ELISA and CFT or CFT-CCT estimates in bTB-infected animals (rohD) was surprising, given that the ELISA and the skin tests target different immune responses and therefore a larger degree of independence is often assumed (32, 56).

Prevalence priors elicited from expert opinion were considerably lower than those based on the intradermal test results (median of 0.35 vs. 0.85). That finding could explain, at least in part, the lower posterior estimates for prevalence obtained in LCA-a models compared with those from LCA-b models. The higher posterior prevalence estimates obtained using all models in both studies; along with the fact that the two sampled herds remained infected with high rates of reactors 2 years after this study was completed (data not shown) suggest that bTB-infection was higher than what was estimated using expert opinion in this population. Comparison of results from the two modeling approaches illustrates the potential negative consequences of basing prior distributions exclusively on expert opinion.

Interestingly, estimates for the Se and Sp of the CFT test were lower than those described for the US (33, 34, 37), and more in line with Se values reported in field studies in Australia (35). Likewise, posterior estimates for the serial use of CFT-CCT, remained in the lower end of previous estimates (33, 36, 40, 41). This relatively low accuracy of the intradermal tests in Uruguay suggests that the bTB-control program may suffer from limited Se in heavily infected herds (~53%), what could lead to the persistence of infected animals in the dairy cattle population over time. Possible explanations for this finding include: presence of other infections that could compromise the reliability of bTB diagnostic tests such as Johne's disease, that is prevalent in dairy herds in Uruguay and whose impact on bTB diagnosis has been already suggested in the country (14), a high proportion of animals in an advanced stage of disease (state of anergy) (30, 57), what could be plausible given the high prevalence of infection in the tested herds, along with other factors associated with the performance of the technique itself or the animals tested, which, with the consolidation and intensification of the industry, may have contributed to the re-emergence of bTB observed in the last decade (14).

Posterior estimates for the Se of IGRAb and IGRAc obtained using LCA-a models (Table 3) are in agreement with previous reports suggesting IGRAs are at least as sensitive as intradermal assays (32, 58). IGRAs have two major advantages over intradermal tests, namely, the potential for detecting false negative animals in the skin test (41, 59, 5961), and the opportunity to maximize their sensitivity thanks to the anamnestic effect induced by the inoculation of PPDs when used in combination with intradermal tests (18, 41, 62). The population under study was sampled a post intradermal inoculation of the PPDb, while this time was variable, it could have contributed to an enhanced Se in agreement with previous studies in which IGRAs performance was assessed following intradermal tuberculin testing (63).

The low Se estimated for the intradermal testing protocol currently used in the bTB-control program in Uruguay (i.e., CFT-CCT) suggests that in-vitro tests might be advantageous if used in parallel to improve the Se of the program in these heavily infected herds. Assuming independence among diagnostic tests, as indicated by the correlation coefficients (rhoD and rhoDc) including zero in the dependent models, an estimate of the potential overall Se of the combination in parallel can be computed using the median estimates values as 1–[1–Se(in-vivo)] * [1-Se(in-vitro)] (64). This approximation showed a 19–25% improved Se in the different in-vivo and in-vitro combinations (Supplementary Table 3) which could vastly improve the detection of bTB-infected animals. A slightly higher Sp was obtained for the IGRAc compared with the IGRAb, what could be due to the use of more specific antigens (peptide-cocktail with ESAT-6 and CFP-10) (28, 42, 45, 46) although could be also a product of the different priors used for each test based on available knowledge. Interestingly, Sp of the IGRAc was equivalent to that of CFT-CCT, suggesting that the use of a single assay (IGRAc) could potentially replace serial testing with CFT-CCT for bTB screening in heavily infected dairy herds without increasing the rate of false positives compared with the current strategy in the project, thus avoiding the unnecessary culling of non-infected cattle.

Care should be taken when interpreting our results, since these were obtained in herds with a very high bTB prevalence and therefore may not be easily extrapolated to other epidemiological situations in Uruguay. The low sensitivity of the skin tests found here could be related with the presence of a larger proportion of infected animals in an advanced stage that is expected in herds with a high infectious pressure, as the ones evaluated here. Therefore, the gain in sensitivity that could be expected from the ancillary application of in-vitro tests (Supplementary Table 3) may not be as high in other situations (smaller dairy herds with lower infection levels, beef herds, etc.). Furthermore, other herd-level factors (such as the presence of Johne's disease) could have also contributed to the results found here. Still, highly infected large dairy herds are one of the main issues faced by the bTB eradication program and therefore our priority here was to assess the usefulness of additional diagnostic approaches to complement the current diagnostic strategy. Additional measures to control the spread of bTB in the Uruguayan cattle population include ban of animal movement, more frequent bTB-testing, the pasteurization of milk use for pre-weaned calves.

In conclusion, results found here, irrespective of the modeling approach followed, suggest that the use of IGRAs in Uruguay can dramatically improve the limited Se of the currently used diagnostic strategies based on skin tests, which would require numerous herd tests to eliminate disease from heavily infected dairy herds as the ones analyzed here. The ELISA could also have some potential for detection of bTB-infected animals if used as an ancillary test to skin test in these populations.

Data Availability

All datasets generated for this study are included in the manuscript and/or the Supplementary Files.

Ethics Statement

This study was performed in accordance with guidelines of the animal ethics committee of the University of Uruguay (Comision Honoraria de Experimentacion Animal).

Author Contributions

CP-R, JA, and AP designed the study and drafted the manuscript. AS, XS, and AN coordinated sample collection, field testing, and laboratory analysis. AG assembled the databases and helped in interpretation of the results. CP-R run the statistical analyses. All authors critically revised the manuscript.

Funding

This research was partially supported by Uruguayan Agency of Research and Innovation (ANII in Spanish language), Plan Nacional de Investigación en Salud Animal & Comité de Investigación Científica de la Universidad de la República (PLANISA–CSIC). JA was the recipient of a Ramón y Cajal postdoctoral contract from the Spanish Ministry of Economy, Industry and Competitiveness (MINECO) (RYC-2016-20422).

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

The authors thank the experts who provided insight for the prior distributions, and the veterinarians from the participant company for their willingness and time.

Supplementary Material

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

References

1. Branscum AJ, Gardner IA, Johnson WO. Estimation of diagnostic-test sensitivity and specificity through Bayesian modeling. Prevent Vet Med. (2005) 68:145–63. doi: 10.1016/j.prevetmed.2004.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Spiegelhalter D, Thomas A, Best NG, Gilks W. BUGS: Bayesian Inference Using Gibbs Sampling, version 0.50. MRC Biostatistics Unit, Cambridge (1996).

Google Scholar

3. Suess E, Gardner I, Johnson W. Hierarchical Bayesian model for prevalence inferences and determination of a country's status for an animal pathogen. Prevent Vet Med. (2002) 55:155–71. doi: 10.1016/S0167-5877(02)00092-2

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Al-Mouqatea S, Alkhamis M, Akbar B, Ali A, Al-Aqeel H, Bin-Heji A, et al. Bayesian estimation of ELISA and gamma interferon test accuracy for the detection of bovine tuberculosis in caudal fold test–negative dairy cattle in Kuwait. J Vet Diagn Invest. (2018) 30:468–70. doi: 10.1177/1040638718759574

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Alvarez J, Perez A, Bezos J, Marqués S, Grau A, Saez JL, et al. Evaluation of the sensitivity and specificity of bovine tuberculosis diagnostic tests in naturally infected cattle herds using a Bayesian approach. Vet Microbiol. (2012) 155:38–43. doi: 10.1016/j.vetmic.2011.07.034

PubMed Abstract | CrossRef Full Text | Google Scholar

6. de la Cruz ML, Branscum AJ, Nacar J, Pages E, Pozo P, Perez A, et al. Evaluation of the performance of the IDvet IFN-gamma test for diagnosis of bovine tuberculosis in Spain. Front Vet Sci. (2018) 5:229. doi: 10.3389/fvets.2018.00229

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Valerie-Beau Pucken, Knubben-Schweizer G, Döpfer D, Groll A, Hafner-Marx A, Hörmansdorfer S, et al. Evaluating diagnostic tests for bovine tuberculosis in the southern part of Germany: a latent class analysis. PLoS ONE. (2017) 12:0179847. doi: 10.1371/journal.pone.0179847

CrossRef Full Text | Google Scholar

8. Zinsstag J, Schelling E, Roth F, Kazwala R. Economics of Bovine tuberculosis. In: Thoen CO, Steele JH, Gilsdorf MJ, editors. Mycobacterium Bovis Infection in Animals and Humans. In: Charles O. Thoen, James H. Steele, Michael J. Gilsdorf, editors. (Ames, IA: John Wiley & Sons, Ltd.) (2008). p. 68–83. Available online at: https://onlinelibrary.wiley.com/doi/abs/10.1002/9780470344538.ch9

Google Scholar

9. Collins JD. Tuberculosis in cattle: Strategic planning for the future. Vet Microbiol. (2006) 112:369–81. doi: 10.1016/j.vetmic.2005.11.041

PubMed Abstract | CrossRef Full Text | Google Scholar

10. OIE. Bovine tuberculosis. In: Manual of Diagnostic Tests and Vaccines for Terrestrial Animals 2014. World Organisation for Animal Health. [Internet]. (2009). Available online at: http://www.oie.int/international-standard-setting/terrestrial-manual/access-online/

Google Scholar

11. Casas Olascoaga R. Antecedentes de la tuberculosis bovina en uruguay, periodo marzo 1888 – enero 1998. Servicio Medico Vet Uruguay. (2013) 49:14–30.

Google Scholar

12. Picasso C, Alvarez J, VanderWaal KL, Fernandez F, Gil A, Wells SJ, et al. Epidemiological investigation of bovine tuberculosis outbreaks in Uruguay (2011–2013). Prevent Vet Med. (2017) 138:156–61. doi: 10.1016/j.prevetmed.2017.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

13. WAHIS_OIE. Animal Health Situation. (2014). Available online at: http://www.oie.int/wahis_2/public/wahid.php/Countryinformation/Animalsituation

Google Scholar

14. Picasso-Risso C, Gil A, Nunez A, Suanes A, Macchi V, Salaberry X, et al. Diagnostic interaction between bovine tuberculosis (bTB) and Johne's disease in bTB highly prevalent dairy farms of Uruguay. Vet Animal Sci. (2019) 7:100052. doi: 10.1016/j.vas.2019.100052

CrossRef Full Text | Google Scholar

15. Alvarez J, Bezos J, de la Cruz ML, Casal C, Romero B, Domínguez L, et al. Bovine tuberculosis: within-herd transmission models to support and direct the decision-making process. Res Vet Sci. (2014) 97:S61–8. doi: 10.1016/j.rvsc.2014.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

16. EFSA Panel on Animal Health and Welfare (AHAW), More S, Bøtner A, Butterworth A, Calistri P, Depner K, et al. Ad hoc method for the assessment on listing and categorisation of animal diseases within the framework of the Animal Health Law. EFSA J. (2017) 15:e04783. doi: 10.2903/j.efsa.2017.4783

CrossRef Full Text | Google Scholar

17. EFSA (European Food Safety Authority). Scientific Opinion on the Use of a Gamma Interferon Test for the Diagnosis of Bovine tuberculosis. (2012). Available online at: https://efsa-onlinelibrary-wiley-com.ezp1.lib.umn.edu/doi/pdf/10.2903/j.efsa.2012.2975

Google Scholar

18. Casal, Díez-Guerrier A, Álvarez J, Rodriguez-Campos S, Mateos A, Linscott R, et al. Strategic use of serology for the diagnosis of bovine tuberculosis after intradermal skin testing. Vet Microbiol. (2014) 170:342–51. doi: 10.1016/j.vetmic.2014.02.036

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Waters WR, Buddle BM, Vordermeier HM, Gormley E, Palmer MV, Thacker TC, et al. Development and evaluation of an enzyme-linked immunosorbent assay for use in the detection of Bovine tuberculosis in Cattle. Clin Vaccine Immunol. (2011) 18:1882–8. doi: 10.1128/CVI.05343-11

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Waters WR, Palmer MV, Thacker TC, Bannantine JP, Vordermeier HM, Hewinson RG, et al. Early Antibody Responses to Experimental Mycobacterium bovis Infection of Cattle. Clin Vaccine Immunol. (2006) 13:648–54. doi: 10.1128/CVI.00061-06

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lyashchenko KP, Pollock JM, Colangeli R, Gennaro ML. Diversity of antigen recognition by serum antibodies in experimental Bovine tuberculosis. Infect Immun. (1998) 66:5344.

PubMed Abstract | Google Scholar

22. Lyashchenko K, Colangeli R, Houde M, Al Jahdali H, Menzies D, Gennaro ML. Heterogeneous antibody responses in tuberculosis. Infect Immun. (1998) 66:3936.

PubMed Abstract | Google Scholar

23. Radunz B, Lepper AWD. Suppression of skin reactivity to bovine tuberculin in repeat tests. Austr Vet J. (1985) 62:191–4. doi: 10.1111/j.1751-0813.1985.tb07294.x

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kostoulas P, Nielsen SS, Branscum AJ, Johnson WO, Dendukuri N, Dhand NK, et al. STARD-BLCM: Standards for the Reporting of Diagnostic accuracy studies that use Bayesian Latent Class Models. Prevent Vet Med. (2017) 138:37–47. doi: 10.1016/j.prevetmed.2017.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Rutjes AWS, Reitsma JB, Vandenbroucke JP, Glas AS, Bossuyt PMM. Case–control and two-gate designs in diagnostic accuracy studies. Clin Chem. (2005) 51:1335. doi: 10.1373/clinchem.2005.048595

PubMed Abstract | CrossRef Full Text | Google Scholar

26. MGAP M de GA y P. Legislacion Sanitaria Animal. [Internet]. (1989). Available online at: http://www2.mgap.gub.uy/portal/page.aspx?2,dgsg,dgsg-legislacion-sanitaria

Google Scholar

27. Wood PR, Corner LA, Plackett P. Development of a simple, rapid in vitro cellular assay for bovine tuberculosis based on the production of gamma interferon. Res Vet Sci. (1990) 49:46–9. doi: 10.1016/S0034-5288(18)31044-0

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Vordermeier HM, Whelan A, Cockle PJ, Farrant L, Palmer N, Hewinson RG. Use of synthetic peptides derived from the antigens ESAT-6 and CFP-10 for differential diagnosis of Bovine tuberculosis in Cattle. Clin Diagn Lab Immunol. (2001) 8:571–8. doi: 10.1128/CDLI.8.3.571-578.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Gardner IA, Stryhn H, Lind P, Collins MT. Conditional dependence between tests affects the diagnosis and surveillance of animal diseases. Prevent Vet Med. (2000) 45:107–22. doi: 10.1016/S0167-5877(00)00119-7

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Pollock JM, Welsh MD, McNair J. Immune responses in bovine tuberculosis: towards new strategies for the diagnosis and control of disease. Vet Immunol Immunopathol. (2005) 108:37–43. doi: 10.1016/j.vetimm.2005.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Georgiadis MP, Johnson WO, Gardner IA, Singh R. Correlation-adjusted estimation of sensitivity and specificity of two diagnostic tests. J R Stat Soc Series C (Applied Statistics). (2003) 52:63–76. doi: 10.1111/1467-9876.00389

CrossRef Full Text | Google Scholar

32. 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

33. Norby B, Bartlett PC, Fitzgerald SD, Granger LM, Bruning-Fann CS, Whipple DL, et al. The sensitivity of gross necropsy, caudal fold and comparative cervical tests for the diagnosis of bovine tuberculosis. J Vet Diagn Invest. (2004) 16:126–31. doi: 10.1177/104063870401600206

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Whipple DL, Bolin CA, Davis AJ, Jarnagin JL, Johnson DC, Nabors RS, et al. Comparison of the sensitivity of the caudal fold skin test and a commercial gamma-interferon assay for diagnosis of bovine tuberculosis. Am J Vet Res. (1995) 56:415–9.

PubMed Abstract | Google Scholar

35. Wood P, Corner LA, Rothel JS, Baldock C, Jones S, Cousins D, et al. Field comparison of the interferon-gamma assay and the intradermal tuberculin test for the diagnosis of bovine tuberculosis. Austr Vet J. (1991) 68:286–90. doi: 10.1111/j.1751-0813.1991.tb03254.x

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Nuñez-Garcia J, Downs SH, Parry JE, Abernethy DA, Broughan JM, Cameron AR, 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. Prevent Vet Med. (2018) 153:94–107. doi: 10.1016/j.prevetmed.2017.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Farnham MW, Norby B, Goldsmith TJ, Wells SJ. Meta-analysis of field studies on bovine tuberculosis skin tests in United States cattle herds. Prevent Vet Med. (2012) 103:234–42. doi: 10.1016/j.prevetmed.2011.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

38. USDA:APHIS. Assessment of Risk Factors for Mycobacterium bovis in the United States. Centers for Epidemiology and Animal Health, Fort Collins, CO. (1992).

Google Scholar

39. Francis J, Seiler RJ, Wilkie IW, O'Boyle D, Lumsden MJ, Frost AJ. The sensitivity and specificity of various tuberculin tests using bovine PPD and other tuberculins. Vet Rec. (1978) 103:420–5.

PubMed Abstract | Google Scholar

40. VanderWaal K, Enns EA, Picasso C, Alvarez J, Perez A, Fernandez F, et al. Optimal surveillance strategies for bovine tuberculosis in a low-prevalence country. Sci Rep. (2017) 7:4140. doi: 10.1038/s41598-017-04466-2

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Vordemeier M, Goodchild AV, Whelan A, Ewer K, Clifton-Hadley RS, Williams J, et al. The Bovigam assay as ancillary test to the tuberculin skin test. Government Vet J. (2006) 16:72–80.

Google Scholar

42. Aagaard C, Govaerts M, Meikle V, Vallecillo AJ, Gutierrez-Pabello JA, Suarez-Güemes F, et al. Optimizing antigen cocktails for detection of Mycobacterium bovis in herds with different prevalences of Bovine tuberculosis: ESAT6-CFP10 mixture shows optimal sensitivity and specificity. J Clin Microbiol. (2006) 44:4326. doi: 10.1128/JCM.01184-06

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Gormley E, Doyle MB, Fitzsimons T, McGill K, Collins JD. Diagnosis of Mycobacterium bovis infection in cattle by use of the gamma-interferon (Bovigam®) assay. Vet Microbiol. (2006) 112:171–9.

Google Scholar

44. Ryan T, Buddle B, de Lisle G. An evaluation of the gamma interferon test for detecting bovine tuberculosis in cattle 8 to 28 days after tuberculin skin testing. Res Vet Sci. (2000) 69:57–61.

PubMed Abstract | Google Scholar

45. Flores-Villalva S, Suárez-Güemes F, Espitia C, Whelan AO, Vordermeier M, Gutiérrez-Pabello JA. Specificity of the tuberculin skin test is modified by use of a protein cocktail containing ESAT-6 and CFP-10 in Cattle naturally infected with Mycobacterium bovis. Clin Vaccine Immunol. (2012) 19:797. doi: 10.1128/CVI.05668-11

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Casal C, Bezos J, Díez-Guerrier A, Álvarez J, Romero B, de Juan L, et al. Evaluation of two cocktails containing ESAT-6, CFP-10 and Rv-3615c in the intradermal test and the interferon-γ assay for diagnosis of bovine tuberculosis. Prevent Vet Med. (2012) 105:149–54. doi: 10.1016/j.prevetmed.2012.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Wood PR, Corner LA, Rothel JS, Ripper JL, Fifis T, McCormick BS, et al. A field evaluation of serological and cellular diagnostic tests for bovine tuberculosis. Vet Microbio. (1992) 31:71–9.

PubMed Abstract | Google Scholar

48. Rogan W, Gladen B. Estimating prevalence from the results of screening test. Am J Epidemiol. (1978) 107:71–6. doi: 10.1093/oxfordjournals.aje.a112510

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. (1992) 7:457–72. doi: 10.1214/ss/1177011136

CrossRef Full Text | Google Scholar

50. Lunn D, Spiegelhalter D, Thomas A, Best N. The BUGS project: Evolution, critique and future directions. Stat Med. (2009) 28:3049–67. doi: 10.1002/sim.3680

CrossRef Full Text | Google Scholar

51. Sturtz S, Ligges U, Gelman A. A Package for Running WinBUGS from R. J Statistical Software. (2005).

Google Scholar

52. Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A. Bayesian measures of model complexity and fit. J R Stat Soc Series B. (2002) 64:583–639. doi: 10.1111/1467-9868.00353

CrossRef Full Text | Google Scholar

53. Enøe C, Georgiadis MP, Johnson WO. Estimation of sensitivity and specificity of diagnostic tests and disease prevalence when the true disease state is unknown. Prevent Vet Med. (2000) 45:61–81. doi: 10.1016/S0167-5877(00)00117-3

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Gardner IA. The utility of Bayes' theorem and Bayesian inference in veterinary clinical practice and research. Austr Vet J. (2002) 80:758–61. doi: 10.1111/j.1751-0813.2002.tb11347.x

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Praud A, Boschiroli ML, Meyer L, Garin-Bastuji B, Dufour B. Assessment of the sensitivity of the gamma-interferon test and the single intradermal comparative cervical test for the diagnosis of bovine tuberculosis under field conditions. Epidemiol Infection. (2015) 143:157–66. doi: 10.1017/S0950268814000338

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Drewe JA, Tomlinson AJ, Walker NJ, Delahay RJ. Diagnostic accuracy and optimal use of three tests for tuberculosis in live badgers. PLoS ONE. (2010) 5:e11196. doi: 10.1371/journal.pone.0011196

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Lepper AWD, Pearson CW, Corner LA. Anergy to tuberculin in beef cattle. Austr Vet J. (1977) 53:214–6. doi: 10.1111/j.1751-0813.1977.tb00188.x

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Bezos J, Casal C, Romero B, Schroeder B, Hardegger R, Raeber AJ, et al. Current ante-mortem techniques for diagnosis of bovine tuberculosis. Res Vet Sci. (2014) 97(Suppl.) S44–52. doi: 10.1016/j.rvsc.2014.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Coad M, Clifford D, Rhodes SG, Hewinson RG, Vordermeier HM, Whelan AO. Repeat tuberculin skin testing leads to desensitisation in naturally infected tuberculous cattle which is associated with elevated interleukin-10 and decreased interleukin-1 beta responses. Vet Res. (2010) 41:14. doi: 10.1051/vetres/2009062

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Gormley E, Doyle MB, McGill K, Costello E, Good M, Collins JD. The effect of the tuberculin test and the consequences of a delay in blood culture on the sensitivity of a gamma-interferon assay for the detection of Mycobacterium bovis infection in cattle. Vet Immunol Immunopathol. (2004) 102:413–20. doi: 10.1016/j.vetimm.2004.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Monaghan ML, Doherty ML, Collins JD, Kazda JF, Quinn PJ. The tuberculin test. Vet Microbiol. (1994) 40:111–24. doi: 10.1016/0378-1135(94)90050-7

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Palmer MV, Waters WR, Thacker TC, Greenwald R, Esfandiari J, Lyashchenko KP. Effects of different tuberculin skin-testing regimenson gamma interferon and antibody responses in cattle experimentally infected with Mycobacterium bovis. Clin Vaccine Immunol. (2006) 13:387–94. doi: 10.1128/CVI.13.3.387-394.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Waters WR, Palmer MV, Stafne MR, Bass KE, Maggioli MF, Thacker TC, et al. Effects of serial skin testing with purified protein derivative on the level and quality of antibodies to complex and defined antigens in Mycobacterium bovis-infected cattle. Alexander TS, editor. Clin Vaccine Immunol. (2015) 22:641–9. doi: 10.1128/CVI.00119-15

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Thrusfield M. Veterinary Epidemiology. 3rd ed. University of Edimburg: Blackwell Science Ltd. (2005).

Google Scholar

Keywords: latent class analysis, diagnosis, interferon-gamma release assay, elisa, chronically infected, Uruguay

Citation: Picasso-Risso C, Perez A, Gil A, Nunez A, Salaberry X, Suanes A and Alvarez J (2019) Modeling the Accuracy of Two in-vitro Bovine Tuberculosis Tests Using a Bayesian Approach. Front. Vet. Sci. 6:261. doi: 10.3389/fvets.2019.00261

Received: 02 April 2019; Accepted: 25 July 2019;
Published: 13 August 2019.

Edited by:

Alejandra Victoria Capozzo, National Council for Scientific and Technical Research (CONICET), Argentina

Reviewed by:

Wendy Beauvais, Cornell University, United States
Ignacio De Blas, University of Zaragoza, Spain

Copyright © 2019 Picasso-Risso, Perez, Gil, Nunez, Salaberry, Suanes and Alvarez. 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: Catalina Picasso-Risso, picas001@umn.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.