BackgroundBovine tuberculosis (bTB) poses serious risks to animal welfare and economy, as well as to public health as a zoonosis. Its etiological agent, Mycobacterium bovis, belongs to the Mycobacterium tuberculosis complex (MTBC), a group of genetically monomorphic organisms featured by a remarkably high overall nucleotide identity (99.9%). Indeed, this characteristic is of major concern for correct typing and determination of strain-specific traits based on sequence diversity. Due to its historical economic dependence on cattle production, Uruguay is deeply affected by the prevailing incidence of Mycobacterium bovis. With the world’s highest number of cattle per human, and its intensive cattle production, Uruguay represents a particularly suited setting to evaluate genomic variability among isolates, and the diversity traits associated to this pathogen.ResultsWe compared 186 genomes from MTBC strains isolated worldwide, and found a highly structured population in M. bovis. The analysis of 23 new M. bovis genomes, belonging to strains isolated in Uruguay evidenced three groups present in the country. Despite presenting an expected highly conserved genomic structure and sequence, these strains segregate into a clustered manner within the worldwide phylogeny. Analysis of the non-pe/ppe differential areas against a reference genome defined four main sources of variability, namely: regions of difference (RD), variable genes, duplications and novel genes. RDs and variant analysis segregated the strains into clusters that are concordant with their spoligotype identities. Due to its high homoplasy rate, spoligotyping failed to reflect the true genomic diversity among worldwide representative strains, however, it remains a good indicator for closely related populations.ConclusionsThis study introduces a comprehensive population structure analysis of worldwide M. bovis isolates. The incorporation and analysis of 23 novel Uruguayan M. bovis genomes, sheds light onto the genomic diversity of this pathogen, evidencing the existence of greater genetic variability among strains than previously contemplated.Electronic supplementary materialThe online version of this article (10.1186/s12864-017-4249-6) contains supplementary material, which is available to authorized users.

Highlights In Uruguay (Uy), bTB highly prevalent dairy herds have emerged. Johne's disease (JD) control programs lack at the national level in Uy. In high prevalence coinfected herds, bTB- and JD-diagnosis interact. An integrated control program is needed for high prevalent coinfected dairies in Uy.

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.