Bovine 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. We 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. This 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.

We sequenced maxicircles fromT.cruzistrains representative of the species evolutionary diversity by using long-read sequencing, which allowed us to uncollapse their repetitive regions, finding that their real lengths range from 35 to 50 kb.T.cruzimaxicircles have a common architecture composed of four regions: coding region (CR), AT-rich region, short (SR) and long repeats (LR). Distribution of genes, both in order and in strand orientation are conserved, being the main differences the presence of deletions affecting genes coding for NADH dehydrogenase subunits, reinforcing biochemical findings that indicate that complex I is not functional inT.cruzi. Moreover, the presence of complete minicircles into maxicircles of some strains lead us to think about the origin of minicircles. Finally, a careful phylogenetic analysis was conducted using coding regions of maxicircles from up to 29 strains, and 1108 single copy nuclear genes from all of the DTUs, clearly establishing that taxonomicallyT.cruziis a complex of species composed by group 1 that contains clades A (TcI), B (TcIII) and D (TcIV), and group 2 (1 and 2 do not coincide with groups I and II described decades ago) containing clade C (TcII), being all hybrid strains of the BC type. Three variants of maxicircles exist inT.cruzi: a, b and c, in correspondence with clades A, B, and C from mitochondrial phylogenies. While A and C carry maxicircles a and c respectively, both clades B and D carry b maxicircle variant; hybrid strains also carry the b- variant. We then propose a new nomenclature that is self-descriptive and makes use of both the phylogenetic relationships and the maxicircle variants present inT.cruzi.

Background: Tuberculosis is a serious infectious disease affecting millions of people. In spite of efforts to reduce the disease, increasing antibiotic resistance has contributed to persist in the top 10 causes of death worldwide. In fact, the increased cases of multi (MDR) and extreme drug resistance (XDR) worldwide remains the main challenge for tuberculosis control. Whole genome sequencing is a powerful tool for predicting drug resistance-related variants, studying lineages, tracking transmission, and defining outbreaks. This study presents the identification and characterization of resistant clinical isolates of Mycobacterium tuberculosis including a phylogenetic and molecular resistance profile study by sequencing the complete genome of 24 strains from different provinces of Ecuador. Results: Genomic sequencing was used to identify the variants causing resistance. A total of 15/21 isolates were identified as MDR, 4/21 as pre-XDR and 2/21 as XDR, with three isolates discarded due to low quality; the main sub-lineage was LAM (61.9%) and Haarlem (19%) but clades X, T and S were identified. Of the six pre-XDR and XDR strains, it is noteworthy that five come from females; four come from the LAM sub-lineage and two correspond to the X-class sub-lineage. A core genome of 3,750 genes, distributed in 295 subsystems, was determined. Among these, 64 proteins related to virulence and implicated in the pathogenicity of M. tuberculosis and 66 possible pharmacological targets stand out. Most variants result in nonsynonymous amino acid changes and the most frequent genotypes were identified as conferring resistance to rifampicin, isoniazid, ethambutol, para-aminosalicylic acid and streptomycin. However, an increase in the resistance to fluoroquinolones was detected. Conclusion: This work shows for the first time the variability of circulating resistant strains between men and women in Ecuador, highlighting the usefulness of genomic sequencing for the identification of emerging resistance. In this regard, we found an increase in fluoroquinolone resistance. Further sampling effort is needed to determine the total variability and associations with the metadata obtained to generate better health policies.

Introduction: Tuberculosis (TB) is among the deadliest diseases worldwide, and its impact is mainly due to the continuous emergence of resistant isolates during treatment due to the laborious process of resistance diagnosis, nonadherence to treatment and circulation of previously resistant isolates of Mycobacterium tuberculosis. In this study, we evaluated the performance and functionalities of web-based tools, including Mykrobe, TB-profiler, PhyResSE, KvarQ, and SAM-TB, for detecting resistance in 88 Ecuadorian isolates of Mycobacterium tuberculosis drug susceptibility tested previously. Statistical analysis was used to determine the correlation between genomic and phenotypic analysis. Our results showed that with the exception of KvarQ, all tools had the highest correlation with the conventional drug susceptibility test (DST) for global resistance detection (98% agreement and 0.941 Cohen’s kappa), while SAM-TB, PhyResSE, TB-profiler and Mykrobe had better correlations with DST for first-line drug analysis individually. We also identified that in our study, only 50% of mutations characterized by the web-based tools in the rpoB, katG, embB, pncA, gyrA and rrs regions were canonical and included in the World Health Organization (WHO) catalogue. Our findings suggest that SAM-TB, PhyResSE, TB-profiler and Mykrobe were efficient in determining canonical resistance-related mutations, but more analysis is needed to improve second-line detection. Improving surveillance programs using whole-genome sequencing tools for first-line drugs, MDR-TB and XDR-TB is essential to understand the molecular epidemiology of TB in Ecuador. Importance: Tuberculosis, an infectious disease caused by Mycobacterium tuberculosis, most commonly affects the lungs and is often spread through the air when infected people cough, sneeze, or spit. However, despite the existence of effective drug treatment, patient adherence, long duration of treatment, and late diagnosis have reduced the effectiveness of therapy and increased drug resistance. The increase in resistant cases, added to the impact of the COVID-19 pandemic, has highlighted the importance of implementing efficient and timely diagnostic methodologies worldwide. The significance of our research is in evaluating and identifying a more efficient and user-friendly web-based tool to characterize resistance in Mycobacterium tuberculosis by whole-genome sequencing, which will allow more routine application to improve TB strain surveillance programs locally.

Tuberculosis is a global public health concern, and understanding Mycobacterium tuberculosis transmission routes and genetic diversity of M. tuberculosis is crucial for outbreak control. This study aimed to explore the genomic epidemiology and genetic diversity of M. tuberculosis in Ecuador by analyzing 88 local isolates and 415 public genomes from 19 countries within the Euro-American lineage (L4). Our results revealed significant genomic diversity among the isolates, particularly in the genes related to protein processing, carbohydrate metabolism, lipid metabolism, and xenobiotic biodegradation and metabolism. The population structure analysis showed that sub-lineages 4.3.2/3 (35.4%), 4.1.2.1 (22.7%), 4.4.1 (12.7%), and 4.1.1. (10.7%) were the most prevalent. Phylogenetic and transmission network analyses suggest that these isolates circulating within Ecuador share genetic ties with isolates from other continents, implying historical and ongoing intercontinental transmission events. Our findings underscore the importance of integrating genomic data into public health strategies for tuberculosis control and suggest that enhanced genomic surveillance is essential for understanding and mitigating the global spread of M. tuberculosis. This study provides a comprehensive genomic framework for future epidemiological investigations and control measures targeting M. tuberculosis L4 in Ecuador.