Skip to main content

Spoligotype-specific risk of finding lesions in tissues from cattle infected by Mycobacterium bovis

Abstract

Background

Although the pathogenic effect of members of the Mycobacterium tuberculosis complex in susceptible hosts is well known, differences in clinical signs and pathological findings observed in infected animals have been reported, likely due to a combination of host and pathogen-related factors. Here, we investigated whether Mycobacterium bovis strains belonging to different spoligotypes were associated with a higher risk of occurrence of visible/more severe lesions in target organs (lungs and/or lymph nodes) from infected animals. A large collection of 8889 samples belonging to cattle were classified depending on the presence/absence of tuberculosis-like lesions and its degree of severity. All samples were subjected to culture irrespective of the presence of lesions, and isolates retrieved were identified and subjected to spoligotyping. The association between the presence/severity of the lesions and the isolation of strains from a given spoligotype was assessed using non-parametric tests and Bayesian mixed multivariable logistic regression models that accounted for origin (region and herd) effects.

Results

Results suggested a difference in severity in lesioned samples depending on the strain’s spoligotype. An association between specific spoligotypes and presence of lesions was observed, with a higher risk of finding lesions in animals infected with strains with spoligotypes SB0120, SB0295 and SB1142 compared with SB0121, and in those coming from certain regions in Spain.

Conclusions

Our results suggest that strains belonging to certain spoligotypes may be associated with a higher probability in the occurrence of gross/macroscopic lesions in infected cattle, although these observational findings should be confirmed in further studies that allow accounting for the effect of other possible confounders not considered here, and ultimately through experimental studies.

Peer Review reports

Background

Animal tuberculosis (TB) is a multi-host chronic-debilitating infectious disease caused by members of the Mycobacterium tuberculosis complex (MTC). Even though it is most commonly found in cattle, its causative agents (M. bovis and to a lesser extent M. caprae) are capable of infecting a wide spectrum of species [1, 2], including domestic species such as goat, pig and sheep, as well as wild species (e.g., wild boar, deer and badger) that can act as reservoirs and hamper efforts invested into the eradication of the disease.

Members of the MTC cause a granulomatous-caseous-necrotical lesion that mainly affects the lungs and regional lymph nodes but can also extend to other organs such as liver, spleen, kidneys, mammary glands, pericardium, uterus and brain [3]. Their pathogenicity is mainly based on four aspects: (i) their ability to multiply inside the macrophages of the host, (ii) their resistance to phagocytosis, (iii) the chronical inflammatory response elicited in the host leading to the development of granulomas, and (iv) a slow evolution of the infection. TB lesions can be easily confused with those produced by other pathogens such as Actinomyces spp., Rhodococcus equi or other mycobacteria belonging to the Mycobacterium avium complex [4, 5].

Previous studies have suggested the existence of differences in the presentation and evolution of the disease in infected animals, what could be due to the strain involved (characterized by DVR-spoligotyping, MIRU-VNTR typing or whole-genome sequence), owing to its ability to evoke a different immune response [6], and to factors related with host susceptibility [7,8,9]. Given that the virulence of a pathogen is typically related with its ability to spread [10], the existence of differences in the capacity of MTC strains to cause clinical lesions in infected hosts could have implications in disease eradication given that it could be hypothesized that more virulent strains may be transmitted faster (though could be also detected easier due to the induction of a stronger immune response in the host).

For this reason, we conducted a study to test the hypothesis that the recovery of M. bovis isolates belonging to specific spoligotypes could be associated with an increased probability in detecting lesions/severity of such lesions in an infected host, considering a large collection of cattle samples retrieved in the frame of the Spanish eradication program. This sample collection had been assessed for the presence of macroscopical TB-like lesion to evaluate the relationship between the presence and degree of lesions and the spoligotype involved accounting for other possible covariates.

Methods

Population of study

A total of 8889 samples from lung and lymph nodes coming from cattle collected from 2011 to 2019 were included in the study. Samples originated from 13 autonomous communities in Spain, and the most represented ones were Madrid (n = 3429), Castilla-La Mancha (2966), Aragon (1429) and Valencia (375).

Samples were collected because they originated from animals subjected to slaughter in the frame of the Spanish bovine tuberculosis eradication program (i.e., they were positive in ante-mortem tests). Samples coming from passive surveillance (i.e., abattoir inspection) were discarded to avoid biases (since they were already selected due to the presence of visible lesions).

Laboratory analysis

The presence of lesions in the samples was established qualitatively (yes/no) and, for samples with lesions, semi quantitatively depending on the type of lesion at the laboratory as previously reported [11]: briefly, in the case of lung samples, a scale from 1 to 5 was used, with:

  • Grade 1: small visible lesion of variable size of 1–2 mm.

  • Grade 2: five or fewer lesions less than 1 cm.

  • Grade 3: more than five lesions less than 1 cm, or one bigger than 1 cm.

  • Grade 4: more than one lesion bigger than 1 cm.

  • Grade 5: granulomatous lesions along all lung parenchyma.

In the case of lymph nodes (mediastinal, retropharyngeal, tracheobronchial and prescapular), lesions were scored on a 1–3 scale as follows:

  • Grade 1: five or less visible lesions of variable size of 1–2 mm.

  • Grade 2: more than five lesions less than 1–2 mm, or one bigger than 0.5 cm.

  • Grade 3: widely spread lesion affecting all the lymph node.

Samples were then subjected to bacteriological analysis as previously described [12]. Briefly, tissues collected before 2018 coming from the same animal were pooled and homogenized, independently of the presence of lesions but always including lesioned tissues when found, for subsequent decontamination with a hexadecyl pyridinium chloride solution and finally inoculated either in solid Löwenstein-Jensen with 1% sodium pyruvate and Coletsos media (Difco, Madrid, Spain) at 37 °C for up to 12 weeks. From 2018 onwards pooled tissues were decontaminated with N-acetyl-L-cysteine-sodium hydroxide solution (BBL® MycoPrep™) and inoculated into an automated liquid medium (BD BACTEC™ MGIT™) and Löwenstein-Jensen with 1% sodium pyruvate medium at 37 °C for up to 56 days. The presence of M. bovis was then confirmed in samples in which bacterial growth was observed by conventional PCR [13]. Finally, the spoligotype and MTC species of the isolated strains were determined by DVR-spoligotyping as previously described [14, 15]. Spoligotype profiles were assigned according to the spoligotype database Mbovis.org [16].

Statistical analysis

The proportion of samples from that were lesion and/or culture positive was determined and for the latter the most common spoligotypes found was described. Differences in the observed proportion of lesion-positive samples depending on the spoligotype was then assessed using a chi-square test. The diversity of spoligotypes in culture-positive samples was determined through Simpson’s index of diversity (SI) using the package “vegan” [17] in R [18]. Confidence intervals (CI) for the Simpson’s indexes were estimated through 1000 bootstrap replicates using the package “boot” [19]. Then, the severity of the lesions in lesioned samples as determined by the semiquantitative scores given to lesioned lung [1,2,3,4,5] or lymph node [1,2,3] samples was compared between spoligotypes through a Kruskal-Wallis test followed by a Holm-Bonferroni post-hoc test. Only samples in which spoligotypes were presented in at least 3% of the total population were subjected to the analysis.

Finally, the association between the spoligotype and the presence of lesions in culture-positive samples was evaluated through multivariable mixed logistic regression models. Again, only spoligotypes present in at least 3% of the samples were considered. Models were fitted in a Bayesian framework as:

$$ {logit}^{-1}\left({p}_{ij}\right)={\alpha}_{j\left[i\right]}+{\beta}_1{X}_{i1}+{\beta}_2{X}_{i2}+\cdots +{\beta}_k{X}_{ik} $$

Where pij is the probability of finding TB-like lesions in samples from animal i from herd j, αj[i] is the random effect for the herd j, β1…βk are the coefficients of the covariables assessed, and X1…Xk the values for the covariables itself.

Covariables available included the autonomous region of origin of the samples and the reason of culling (i.e., positive to intradermotuberculinization [IDTB] or gamma-interferon [IFN-γ]). Herd was included as a random variable to account for the lack of independence of samples with the same origin.

Weakly informative normal N (0, 100) prior distributions were used for the β coefficients. The α random effect was assumed to follow a normal N(μ, σ− 2) distribution, with μ ~ N (0, 1) and σ ~ uniform (0, 10). The best model was selected based on the lowest DIC (Deviance Information Criteria) [20].

Models were fitted in WinBUGS [20] through the R package “R2WinBUGS” [21]. Three Markov Chain Monte Carlo chains were run for 100,000 iterations, with a “burn-in” of the first 25,000 iterations, and posterior distributions were obtained after thinning every 10 iteration. Convergence was assessed visually and more formally using the Gelman-Rubin statistic [22].

Results

Descriptive analysis

Macroscopical lesions were detected in approximately one-fourth of the samples included in the study (n = 2187/8889; 24.6%). Among the 8889 samples, 2080 (23.4%) were culture-positive and subsequently spoligotyped, leading to the identification of 88 different spoligotypes. Among the culture-positive samples, most presented TB-like lesions (1835/2080; 88.2%). The most common spoligotype retrieved from cattle samples in the collection was SB0121, the most prevalent in the Iberian Peninsula [21]. Overall, spoligotype diversity in the isolate collection was high (SI = 0.871, 95% CI 0.805–0.944), and if only one isolate per herd was considered, the estimated SI did not vary significantly (0.875, 95% CI 0.818–0.936), although the relative proportion of certain spoligotypes was different (data not shown). Overall, seven of the spoligotypes (SB0121, SB0339, SB0134, SB0265, SB0295, SB0120 and SB1142) were isolated in more than 3% of all the samples (and represented 73.8% of all culture-positive samples).

The proportion of lesions in culture-positive samples between spoligotypes was significantly different (p ≤ 0.001). SB0120 was the spoligotype in which the highest proportion of lesioned samples was observed (95/98; 96.9%), closely followed by SB1142 (70/74; 94.6%) and SB0295 (98/104; 94.2%), while lesions were less frequent in samples that yielded isolates with the SB0265 (112/129; 86.8%) and SB0339 (351/404; 86.9%) profiles. A total of 1986 lesioned TB target-organ samples from 1600 cattle were evaluated semi-quantitatively depending on the type of lesion. Lesions assessed were found in lungs (n = 311), mediastinal (834), retropharyngeal (639), tracheobronchial (150), and prescapular lymph nodes (52) (Additional file 1). The most prevalent type of lesion in lungs was grade 3 (85/311; 27.3%). In the case of the lymph nodes, grade 3 was the most common score (845/1675; 51.0%) (Table 1). Significant differences between the type of lesion observed and the spoligotype present were only found in lungs (p = 0.029) and mediastinal lymph nodes (p = 0.027), although differences between specific spoligotypes were not detected using post-hoc tests.

Table 1 Scores observed in lesioned samples (n = 1986) from lung and lymph nodes from cattle

Multivariable analysis

A total of 1534 samples were included in the multivariable analysis. Animals originated from 423 herds (mean = 3.63 animals per herd, median = 2, range = 1–38). All isolates belonged to seven spoligotypes that had been retrieved in more than 3% of the samples (SB0121, SB0339, SB0134, SB0265, SB0295, SB0120, and SB1142) (1534/2080; 73.8% of all cattle samples). Nearly all (93.9%, 1441/1534) samples came from 4 regions (Madrid, Castilla-La Mancha, Aragon and Valencia), with the remaining regions (Balearic Islands, Andalucia, La Rioja, Extremadura, Navarra, Murcia, Canary Islands and Cataluña) combined in a miscellaneous category (Table 2) (Additional file 2).

Table 2 Number of lesioned cattle samples and its association with the spoligotype and region of origin

The final model for cattle included the spoligotype and the region of origin as covariables. Two spoligotypes, SB0120 and SB0295, were associated with an increased probability of finding lesions (Table 2): animals infected with strains belonging to spoligotypes SB0120 and SB0295 had 6.1 (Posterior Probability Interval (PPI) 95% [1.96–27.9]) and 2.7 (PPI 95% [1.1–7.1]) higher odds, respectively, of having a detected lesion compared to the reference (SB0121) (Fig. 1). In addition, an increased risk was observed in samples from animals infected with SB1142 while the opposite was true (i.e., there was a decreased risk) for samples infected with SB0339, although in those cases the 95% PPI included 1 (Table 2).

Fig. 1
figure1

Posterior probability interval of finding a lesion in cattle samples depending on the spoligotype isolated. Spoligotypes were compared with the reference category (SB0121) according to a Bayesian mixed logistic regression model

In addition, samples coming from the regions of Castilla-La Mancha, Aragon and Valencia had lower odds of having lesions than those coming from the reference category (Madrid), although the 95% PPI for samples coming from Castilla-La Mancha included 1 (Table 2).

Discussion

Animal tuberculosis is one of the most important diseases impacting wildlife and livestock worldwide due to the problems associated to control it, related to its complex epidemiology [23, 24]. Molecular characterization of its causative agent is now performed routinely in many countries (usually through spoligotyping, although whole genome sequencing is increasingly applied) in order to help to identify sources of infection and control disease spread [25, 26]. Differences in the pathogenicity/virulence of MTC strains could have important consequences in disease control, since a more virulent strain may induce the formation of larger lesions in a shorter period (thus increasing transmissibility) [27]. Furthermore, an increased survival in the environment could also provide a competitive advantage to specific strains, since indirect transmission (i.e., due to contact with contaminated environments instead of direct contact with an infected animal) can lead to a significant proportion of new infections in certain areas of Spain [28,29,30]. However, an increased virulence may also induce a stronger immune response that would be more easily detected in host species subjected to routine diagnostics tests [31], in agreement with previous reports describing a positive correlation between the reaction size of the single intradermal comparative cervical test (SICCT) and the number and size of lesions found [32]. Nevertheless, other studies have found no association between genotype of strains and the outcome of SICCT and therefore this relationship may not (always) hold [33]. Here, we assessed the presence of an association between a certain strain type (as defined by the spoligotype) and the presence and severity of the lesions in a large collection of samples coming from cattle retrieved in the frame of the Spanish TB eradication program [34].

Because of the observational nature of this study, several limitations must be considered when interpreting the results. First, we used the spoligotyping profile as our definition of strain, even though it is determined considering only a small fraction of the genome of M. bovis (differences in the Direct Repeat locus, a non-coding region), and thus, distantly related strains may present the same spoligotype while more closely related isolates may have different profiles [35, 36], limiting its discriminatory power, and even more with predominant spoligotypes as in this study. Further studies involving WGS are therefore needed to identify if differences in the spoligotype may be associated with genetic traits that could lead to increased/decreased pathogenicity.

In addition, several variables that could be associated with an increased risk of lesion in infected animals (not related with the specific M. bovis strain) and that were not formally included in the analysis could have influenced our results. For example, given the relatively long period (2011–2019) and regions included in the study, there may have been differences in the post-mortem inspection at the abattoir leading to the samples collection analysed here. Nevertheless, given that the presence of lesions/lesion severity was assessed at the same laboratory by the same technicians over time once the samples were collected there should not be differences in the criterium used to score the tissues.

Unfortunately, we had no information of other host-related factors (such as age, breeding type, sex, presence of coinfections, history of tuberculosis in the herd) that could be related with the early/late development of lesions. Also, the fact that only differences in the severity of lesions were seen in organs with a high number of samples could be mainly due to a lack of sample size, impairing our ability to detect significant differences. Nevertheless, accounting for the unit and the region of origin could partially control for the confounding effect of some of these unknown variables. In any case our results must therefore be interpreted carefully, and should be confirmed in animals for which more information is available and, ideally, in experimental conditions.

According to our results, the proportion of culture-positive samples was relatively low (23.4%). This was not unexpected even if all samples originated from animals that were positive to the official TB diagnostic tests (i.e., IDTB and IFN-γ) and were thus presumably infected: bacteriological culture has a limited sensitivity [37], particularly in the early stages of the disease [38]. Since the majority of the cattle population over 6 weeks of age are subjected to regular routine testing, infected animals are typically detected shortly after exposure to pathogen, and in this context ante-mortem tests are expected to outperform bacteriology, which should not be used as a gold-standard to determine the true infectious status of an animal [39].

When only the subset of culture-positive samples was considered, a much higher proportion of lesioned samples was found as expected (88.2%). In turn, most of the lesioned samples were also culture-positive (1755/2187; 80.2%), thus demonstrating the superior performance that culture has to confirm infection in animals in a more advance stage of disease [35, 36].

According to the Bayesian mixed multivariable logistic regression model, a strong effect of the spoligotype was observed, with profiles associated with both high (SB0120, SB0295, SB1142) and low (SB0339) risk of finding lesions in infected animals compared to the reference profile (SB0121) (Table 2). Interestingly, the SB0120 profile (BCG-like spoligotype) had the strongest association with the presence of lesions (Odds = 6.1 [95% PPI 1.96–27.9]). This spoligotype was also considered in a previous study comparing prevalent spoligotypes in cattle in Argentina [40], but in that case it was not associated with increased risk of lesions. However, in that study only SB0120 out of the seven spoligotypes evaluated here was considered. Furthermore, samples considered in this study originated from cattle identified as infected during abattoir inspection, a population discarded here.

Regarding the findings related to the severity of TB lesions according to a semiquantitative score, samples infected with SB0120 did not have higher scores than those with other spoligotypes in which higher proportions of mediastinal lymph nodes and lung samples reached the maximum degree severity (Additional file 3): most of the lesions due to SB0120 in these organs were of grade one out of three (mediastinal lymph node, 23/53; 43.4%) and grade three out of five (lungs, 7/17; 41.2%). Interestingly, the other spoligotype with an increased risk of detection of lesions (SB0295) was associated with an increased severity in lesions found in both mediastinal lymph nodes and lung (Additional file 3), although significant differences between specific spoligotypes were not confirmed in the post-hoc analysis.

The relative frequency of the spoligotypes found in samples was determined including more than one isolate from a given epidemiological unit (i.e., herd) and thus could be affected by oversampling of certain units infected with specific spoligotypes. Therefore, it should not be considered as a representative sample of the true underlying diversity in the target populations sampled. Still, when only one isolate per unit was considered a similar diversity and proportion of spoligotypes was observed, suggesting that this was not a major source of bias in our results. Furthermore, the inclusion of epidemiological unit in the multivariable model corrected for the lack of independence between observations coming from the same herd.

Conclusions

This preliminary analysis suggests that the infection by M. bovis strains may lead to the occurrence of lesions preferentially when certain strains (spoligotypes) are involved. Therefore, in-depth studies to determine how the strain can influence disease progression accounting for host-related factors are needed to further understand the epidemiology of the infection, especially with regards to the immune response triggered by it, a key factor for early diagnosis.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

IDTB :

Intradermotuberculinization

IFN-γ:

Gamma-interferon

MTC:

Mycobacterium tuberculosis complex

PPI:

Posterior probability interval

SI:

Simpson’s diversity index

SICCT:

Single intradermal comparative cervical test

TB:

Animal tuberculosis

References

  1. 1.

    Aranaz A, de Juan L, Montero N, Sánchez C, Galka M, Delso C, et al. Bovine tuberculosis (Mycobacterium bovis) in wildlife in Spain. J Clin Microbiol. 2004;42(6):2602–8. https://doi.org/10.1128/JCM.42.6.2602-2608.2004.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Liébana E, Aranaz A, Francis B, Cousins D. Assessment of genetic markers for species differentiation within the Mycobacterium tuberculosis complex. J Clin Microbiol. 1996;34(4):933–8. https://doi.org/10.1128/JCM.34.4.933-938.1996.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Zachary JF. Pathologic basis of veterinary disease. Sixth. St. Louis, Missouri: Elsevier, Inc.; 2017.

    Google Scholar 

  4. 4.

    Vicente J, Höfle U, Garrido JM, Fernandez de Mera IG, Juste R, Barral M, et al. Wild boar and red deer display high prevalences of tuberculosis-like lesions in Spain. Vet Res. 2005;37:107–19.

    Article  Google Scholar 

  5. 5.

    Savini V, Fazii P, Favaro M, Astolfi D, Polilli E, Pompilio A, et al. Tuberculosis-like pneumonias by the aerobic actinomycetes Rhodococcus, Tsukamurella and Gordonia. Microbes and Infection. 2012;14(5):401–10. https://doi.org/10.1016/j.micinf.2011.11.014.

    Article  PubMed  Google Scholar 

  6. 6.

    Aguilar León D, Zumárraga MJ, Jiménez Oropeza R, Gioffré AK, Bernardelli A, Orozco Estévez H, et al. Mycobacterium bovis with different genotypes and from different hosts induce dissimilar immunopathological lesions in a mouse model of tuberculosis. Clin Exp Immunol. 2009;157(1):139–47. https://doi.org/10.1111/j.1365-2249.2009.03923.x.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Zanella G, Duvauchelle A, Hars J, Moutou F, Boschiroli ML, Durand B. Patterns of lesions of bovine tuberculosis in wild red deer and wild boar. Vet Rec. 2010;163(2):43–7.

    Article  Google Scholar 

  8. 8.

    Bezos J, Álvarez J, de Juan L, Romero B, Rodríguez S, Fernández-de-Mera IG, et al. Assessment of in vivo and in vitro tuberculosis diagnostic tests in Mycobacterium caprae naturally infected caprine flocks. Prev Vet Med. 2011;100(3–4):187–92. https://doi.org/10.1016/j.prevetmed.2011.03.012.

    Article  PubMed  Google Scholar 

  9. 9.

    Vicente J, Höfle U, Garrido JM, Fernández-de-Mera IG, Acevedo P, Juste R, et al. Risk factors associated with the prevalence of tuberculosis-like lesions in fenced wild boar and red deer in south Central Spain. Vet Res. 2007;38(3):451–64. https://doi.org/10.1051/vetres:2007002.

    Article  PubMed  Google Scholar 

  10. 10.

    André J-B, Hochberg ME. Virulence evolution in emerging infectious diseases. Evolution. 2005;59(7):1406–12. https://doi.org/10.1111/j.0014-3820.2005.tb01791.x.

    Article  PubMed  Google Scholar 

  11. 11.

    Vordermeier HM, Chambers MA, Cockle PJ, Whelan AO, Simmons J, Hewinson RG. Correlation of ESAT-6-specific gamma interferon production with pathology in cattle following Mycobacterium bovis BCG vaccination against experimental bovine tuberculosis. Infect Immun. 2002;70(6):3026–32. https://doi.org/10.1128/IAI.70.6.3026-3032.2002.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Lorente-Leal V, Liandris E, Castellanos E, Bezos J, Domínguez L, de Juan L, et al. Validation of a real-time PCR for the detection of Mycobacterium tuberculosis complex members in bovine tissue samples. Front Vet Sci. 2019;6(61):1–9.

    Google Scholar 

  13. 13.

    Wilton S, Cousins D. Detection and identification of multiple mycobacterial pathogens by DNA amplification in a single tube. Genome Res. 1992;1(4):269–73. https://doi.org/10.1101/gr.1.4.269.

    CAS  Article  Google Scholar 

  14. 14.

    Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997;35(4):907–14. https://doi.org/10.1128/JCM.35.4.907-914.1997.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Aranaz A, Liébana E, Mateos A, Dominguez L, Vidal D, Domingo M, et al. Spacer oligonucleotide typing of Mycobacterium bovis strains from cattle and other animals: a tool for studying epidemiology of tuberculosis. J Clin Microbiol. 1996;34(11):2734–40. https://doi.org/10.1128/JCM.34.11.2734-2740.1996.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Smith NH, Upton P. Naming spoligotype patterns for the RD9-deleted lineage of the Mycobacterium tuberculosis complex; https://www.Mbovis.org. Infect Genet Evol. 2012;12(4):873-6. https://doi.org/10.1016/j.meegid.2011.08.002.

  17. 17.

    Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, Mcglinn D, et al. vegan. Community ecology package. R package version 2.5–6; 2019. p. 1–296. Available from: https://github.com/vegandevs/vegan

    Google Scholar 

  18. 18.

    R Core Team. R: A language and environment for statistical computing. Viena: R Foundation for statistical Computing; 2019. Available from: https://www.R-project.org/

    Google Scholar 

  19. 19.

    Canty A, Ripley BD. boot: Bootstrap R (S-Plus) Functions. R package version 1.3–26; 2021.

    Google Scholar 

  20. 20.

    Lunn DJ, Thomas A, Best N, Spiegelhalter D. WinBugs - a Bayesian modelling framework: concepts, structure and extensibility. Stat Comput. 2000;10(4):325–37. https://doi.org/10.1023/A:1008929526011.

    Article  Google Scholar 

  21. 21.

    Sturtz S, Ligges U, Gelman A. R2WinBUGS: a package for running WinBUGS from R. J Stat Softw. 2005;12(3):1–16.

    Article  Google Scholar 

  22. 22.

    Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7(4):457–511.

    Google Scholar 

  23. 23.

    Pesciaroli M, Alvarez J, Boniotti MB, Cagiola M, di Marco V, Marianelli C, et al. Tuberculosis in domestic animal species. Res Vet Sci. 2014;97:78–85.

    Article  Google Scholar 

  24. 24.

    Pozo P, VanderWaal K, Grau A, de la Cruz ML, Nacar J, Bezos J, et al. Analysis of the cattle movement network and its association with the risk of bovine tuberculosis at the farm level in Castilla y Leon, Spain. Transbound Emerg Dis. 2019;66(1):327–40. https://doi.org/10.1111/tbed.13025.

    Article  PubMed  Google Scholar 

  25. 25.

    Crispell J, Zadoks RN, Harris SR, Paterson B, Collins DM, de-Lisle GW, et al. Using whole genome sequencing to investigate transmission in a multi-host system: Bovine tuberculosis in New Zealand. BMC Genomics. 2017;18(1):180. https://doi.org/10.1186/s12864-017-3569-x.

  26. 26.

    Trewby H, Wright D, Breadon EL, Lycett SJ, Mallon TR, McCormick C, et al. Use of bacterial whole-genome sequencing to investigate local persistence and spread in bovine tuberculosis. Epidemics. 2016;14:26–35. https://doi.org/10.1016/j.epidem.2015.08.003.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Acevedo MA, Dillemuth FP, Flick AJ, Faldyn MJ, Elderd BD. Virulence-driven trade-offs in disease transmission: a meta-analysis. Evolution. 2019;73(4):636–47. https://doi.org/10.1111/evo.13692.

    Article  PubMed  Google Scholar 

  28. 28.

    Cowie CE, Hutchings MR, Barasona JA, Gortázar C, Vicente J, White PCL. Interactions between four species in a complex wildlife: livestock disease community: implications for Mycobacterium bovis maintenance and transmission. Eur J Wildl Res. 2016;62(1):51–64. https://doi.org/10.1007/s10344-015-0973-x.

    Article  Google Scholar 

  29. 29.

    Triguero-Ocaña R, Barasona JA, Carro F, Soriguer RC, Vicente J, Acevedo P. Spatio-temporal trends in the frequency of interspecific interactions between domestic and wild ungulates from Mediterranean Spain. PLoS One. 2019;14(1):e0211216. https://doi.org/10.1371/journal.pone.0211216.

  30. 30.

    Allen AR, Ford T, Skuce RA. Does Mycobacterium tuberculosis var. bovis survival in the environment confound bovine tuberculosis control and eradication? A literature review. Morck D, editor. Vet Med Int. 2021;2021:1–19 Available from: https://www.hindawi.com/journals/vmi/2021/8812898/.

    Article  Google Scholar 

  31. 31.

    Bezos J, Casal C, Díez-Delgado I, Romero B, Liandris E, Álvarez J, Sevilla IA, Juan L, Domínguez L, Gortázar C Goats challenged with different members of the Mycobacterium tuberculosis complex display different clinical pictures. Vet Immunol Immunopathol 2015;167:185–189, 3-4, doi: https://doi.org/10.1016/j.vetimm.2015.07.009.

  32. 32.

    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(3):844–58. https://doi.org/10.1111/tbed.12814.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Wright DM, Allen AR, Mallon TR, McDowell SWJ, Bishop SC, Glass EJ, et al. Detectability of bovine TB using the tuberculin skin test does not vary significantly according to pathogen genotype within Northern Ireland. Infect Genet Evol. 2013;19:15–22. https://doi.org/10.1016/j.meegid.2013.05.011.

    Article  PubMed  Google Scholar 

  34. 34.

    Rodríguez S, Romero B, Bezos J, de Juan L, Álvarez J, Castellanos E, et al. High spoligotype diversity within a Mycobacterium bovis population: clues to understanding the demography of the pathogen in Europe. Vet Microbiol. 2010;141(1-2):89–95. https://doi.org/10.1016/j.vetmic.2009.08.007.

    Article  PubMed  Google Scholar 

  35. 35.

    Smith NH. The global distribution and phylogeography of Mycobacterium tuberculosis clonal complexes. Infect Genet Evol. 2012;12(4):857–65. https://doi.org/10.1016/j.meegid.2011.09.007.

    Article  PubMed  Google Scholar 

  36. 36.

    Rodriguez-Campos S, Navarro Y, Romero B, de Juan L, Bezos J, Mateos A, et al. Splitting of a prevalent Mycobacterium bovis spoligotype by variable-number tandem-repeat typing reveals high heterogeneity in an evolving clonal group. J Clin Microbiol. 2013;51(11):3658–65. https://doi.org/10.1128/JCM.01271-13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    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. Prev Vet Med. 2018;153:94–107. https://doi.org/10.1016/j.prevetmed.2017.02.017.

    Article  PubMed  Google Scholar 

  38. 38.

    Gormley E, Corner LAL, Costello E, Rodriguez-Campos S. Bacteriological diagnosis and molecular strain typing of Mycobacterium bovis and Mycobacterium caprae. Res Vet Sci. 2014;97:S30–40. https://doi.org/10.1016/j.rvsc.2014.04.010.

    Article  PubMed  Google Scholar 

  39. 39.

    de la Cruz ML, Pozo P, Grau A, Nacar J, Bezos J, Perez A, et al. Assessment of the sensitivity of the bovine tuberculosis eradication program in a high prevalence region of Spain using scenario tree modeling. Prev Vet Med. 2019;173:104800. https://doi.org/10.1016/j.prevetmed.2019.104800.

  40. 40.

    Garbaccio S, Macias A, Shimizu E, Paolicchi F, Pezzone N, Magnano G, et al. Association between spoligotype-VNTR types and virulence of Mycobacterium bovis in cattle. Virulence. 2014;5(2):297–302. https://doi.org/10.4161/viru.27193.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

The authors would like to thank the official veterinary services from the central government (Ministerio de Agricultura, Pesca y Alimentación) and from the Autonomous Regions for their collaboration and support.

Funding

This work was funded by the project Integrated Strategies for tuberculosis control and eradication in Spain (ERATUB) (RTI2018–096010-B-C22, Ministerio de Ciencia, Innovación y Universidades, MICINN) and is a contribution to projects RTA2015–00043-C02–02 (Ministerio de Economía y Competitividad (MINECO) and PID2019-105155RB-C31 (MICINN). JA is the recipient of a Ramón y Cajal contract from the Spanish Ministry of Economy, Industry and Competitiveness (MINECO, RYC-2016-20422).

Author information

Affiliations

Authors

Contributions

BR, LdJ and JA designed the study. JB, FL and AGB participated on laboratory analysis and data collection. AGB and JA performed the statistical analyses and drafted the first version of the manuscript. All authors provided written feedback and read and approved the final version of the manuscript.

Corresponding author

Correspondence to Julio Alvarez.

Ethics declarations

Ethics approval and consent to participate

No ethical approval and or permissions were required since all samples used in this study were coming from routine activities in the eradication program for bovine tuberculosis and no specific sampling for this study was carried out. Handling of animals and sample collection was performed according to Spanish and European Legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Scored lesioned lymph node samples.

Additional file 2.

Distribution of the most isolated spoligotypes in cattle within different regions from Spain.

Additional file 3.

Frequency of isolates belonging to the most frequent spoligotypes cultured from lesioned mediastinal lymph nodes and lungs from cattle depending on the type of lesion found.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gomez-Buendia, A., Romero, B., Bezos, J. et al. Spoligotype-specific risk of finding lesions in tissues from cattle infected by Mycobacterium bovis. BMC Vet Res 17, 148 (2021). https://doi.org/10.1186/s12917-021-02848-3

Download citation

Keywords

  • Animal tuberculosis
  • Mycobacterium bovis
  • Spoligotype
  • Lesion