Genetic diversity of Mycobacterium tuberculosis from Pará, Brazil, reveals a higher frequency of ancestral strains than previously reported in South America

Genomic Diversity of the Rarely Observed Genotype of the Mycobacterium tuberculosis Central Asian (CAS) Lineage 3 from North Brazil

Mycobacterium tuberculosis (Mtb) Central Asian Strain (CAS) Lineage 3 (L3) genotype is predominantly found in East-Africa, Central-Asia, Western-Asia, and South-Asia; however, a new spoligotyping CAS/SIT2545 was found in northern regions of Brazil. We aimed to characterize and describe the genetic diversity and perform a phylogenetic assessment of this novel genotype. We performed 24-MIRU-VNTR loci and Whole-genome sequencing (WGS) of six Brazilian isolates previously spoligotyped. The libraries were prepared using a Nextera-XT kit and sequenced in a NextSeq 550 Illumina instrument. We performed lineage assignment and genomic characterization. From publicly available genomes of Mtb L3 and other lineages, we created a robust dataset to run the MTBSeq pipeline and perform a phylogenetic analysis. MIRU-VNTR and WGS confirmed CAS/SIT2545 belongs to L3. Out of 1691 genomes, 1350 (79.83%) passed in quality control (genomic coverage > 95%). Strain 431 differed in 52 single nucleotide variants (SNV), confirming it does not belong to the same transmission chain. The eight genomes from a global dataset clustered closer to Brazilian strains differed in >52 SNVs. We hypothesized L3 and L1 were introduced in Brazilian Northern in the same historical event; however, there is a need for additional studies exploring the genetic diversity of Mtb Brazilian Northern.

Genome-wide estimation of recombination, mutation and positive selection enlightens diversification drivers of Mycobacterium bovis

Genome sequencing has reinvigorated the infectious disease research field, shedding light on disease epidemiology, pathogenesis, host-pathogen interactions and also evolutionary processes exerted upon pathogens. Mycobacterium tuberculosis complex (MTBC), enclosing M. bovis as one of its animal-adapted members causing tuberculosis (TB) in terrestrial mammals, is a paradigmatic model of bacterial evolution. As other MTBC members, M. bovis is postulated as a strictly clonal, slowly evolving pathogen, with apparently no signs of recombination or horizontal gene transfer. In this work, we applied comparative genomics to a whole genome sequence (WGS) dataset composed by 70 M. bovis from different lineages (European and African) to gain insights into the evolutionary forces that shape genetic diversification in M. bovis. Three distinct approaches were used to estimate signs of recombination. Globally, a small number of recombinant events was identified and confirmed by two independent methods with solid support. Still, recombination reveals a weaker effect on M. bovis diversity compared with mutation (overall r/m = 0.037). The differential r/m average values obtained across the clonal complexes of M. bovis in our dataset are consistent with the general notion that the extent of recombination may vary widely among lineages assigned to the same taxonomical species. Based on this work, recombination in M. bovis cannot be excluded and should thus be a topic of further effort in future comparative genomics studies for which WGS of large datasets from different epidemiological scenarios across the world is crucial. A smaller M. bovis dataset (n = 42) from a multi-host TB endemic scenario was then subjected to additional analyses, with the identification of more than 1,800 sites wherein at least one strain showed a single nucleotide polymorphism (SNP). The majority (87.1%) was located in coding regions, with the global ratio of non-synonymous upon synonymous alterations (dN/dS) exceeding 1.5, suggesting that positive selection is an important evolutionary force exerted upon M. bovis. A higher percentage of SNPs was detected in genes enriched into "lipid metabolism", "cell wall and cell processes" and "intermediary metabolism and respiration" functional categories, revealing their underlying importance in M. bovis biology and evolution. A closer look on genes prone to horizontal gene transfer in the MTBC ancestor and included in the 3R (DNA repair, replication and recombination) system revealed a global average negative value for Taijima's D neutrality test, suggesting that past selective sweeps and population expansion after a recent bottleneck remain as major evolutionary drivers of the obligatory pathogen M. bovis in its struggle with the host.

Molecular epidemiology of Mycobacterium tuberculosis in Brazil before the whole genome sequencing era: a literature review

Molecular-typing can help in unraveling epidemiological scenarios and improvement for disease control strategies. A literature review of Mycobacterium tuberculosis transmission in Brazil through genotyping on 56 studies published from 1996-2019 was performed. The clustering rate for mycobacterial interspersed repetitive units - variable tandem repeats (MIRU-VNTR) of 1,613 isolates were: 73%, 33% and 28% based on 12, 15 and 24-loci, respectively; while for RFLP-IS6110 were: 84% among prison population in Rio de Janeiro, 69% among multidrug-resistant isolates in Rio Grande do Sul, and 56.2% in general population in São Paulo. These findings could improve tuberculosis (TB) surveillance and set up a solid basis to build a database of Mycobacterium genomes.

Evaluation of drug susceptibility profile of Mycobacterium tuberculosis Lineage 1 from Brazil based on whole genome sequencing and phenotypic methods

BACKGROUND

The evaluation of procedures for drug susceptibility prediction of Mycobacterium tuberculosis based on genomic data against the conventional reference method test based on culture is realistic considering the scenario of growing number of tools proposals based on whole-genome sequences (WGS).

OBJECTIVES

This study aimed to evaluate drug susceptibility testing (DST) outcome based on WGS tools and the phenotypic methods performed on isolates of M. tuberculosis Lineage 1 from the state of Pará, Brazil, generally associated with low levels of drug resistance.

METHODOLOGY

Culture based DST was performed using the Proportion Method in Löwenstein-Jensen medium on 71 isolates that had been submitted to WGS. We analysed the seven main genome sequence-based tools for resistance and lineage prediction applied to M. tuberculosis and for comparison evaluation we have used the Kappa concordance test.

FINDINGS

When comparing the WGS-based tools against the DST, we observed the highest level of agreement using TB-profiler. Among the tools, TB-profiler, KvarQ and Mykrobe were those which identified the largest number of TB-MDR cases. Comparing the four most sensitive tools regarding resistance prediction, agreement was observed for 43 genomes.

MAIN CONCLUSIONS

Drug resistance profiling using next-generation sequencing offers rapid assessment of resistance-associated mutations, therefore facilitating rapid access to effective treatment.

Characterization of Mycobacterium tuberculosis var. africanum isolated from a patient with pulmonary tuberculosis in Brazil

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC), including Mycobacterium tuberculosis var. tuberculosis (MTB) and Mycobacterium tuberculosis var. africanum (MAF). While MTB is isolated worldwide, MAF is almost completely restricted to the African continent, and despite the historical proximity between Brazil and Africa during the slave trade, no case of TB being caused by MAF has been reported in Brazil to date. We hereby describe the first case of TB caused by MAF in Brazil comparing its genome against the published ones. A female patient who had never visited Africa presented with clinical symptoms typical of pulmonary TB. Based on 16S rRNA gene sequencing, the cultured isolate was identified as belonging to MTBC and partial sequence of the hsp65 gene was identical to that of MAF. This was confirmed by genotyping based on detection of Single Nucleotide Polymorphism (SNP), Region of Difference (RD) and spoligotyping. The isolate presented the Shared International Typing (SIT) 181. In the whole-genome comparison against MAF genomes available on published EMBL-EBI European Nucleotide Archive (ENA), the Brazilian genome (MAFBRA00707) was identified as belonging to Lineage 6 and clustered with isolates from The Gambia. This is the first report of the isolation of MAF from a patient from Brazil, without evidence of having any contact with an African index case.

Laboratory tools for tuberculosis control in a setting with a high burden of HIV/AIDS

Introduction. Nosocomial transmission of Mycobacterium tuberculosis is an important health issue and the detection of tuberculosis (TB) cases is the main tool for controlling this disease.Aim. We aimed to assess the possible occurrence of nosocomial transmission of M. tuberculosis in a reference hospital for HIV/AIDS patients and evaluate both the performance of the Xpert MTB/RIF (Xpert) platform and drug resistance profiles.Methodology. We evaluated the performance of the Xpert platform. Samples that tested positive on the BACTEC MGIT 320 (MGIT320) platform were submitted for genotyping and drug susceptibility testing.Results. In this study, pulmonary and extrapulmonary samples from 407 patients were evaluated, and among these, 15.5 % were diagnosed with TB by the MGIT320 platform, with a TB/HIV coinfection rate of 52.4 %. The Xpert platform gave positive results for TB for 11 samples with negative results on the MGIT320 platform. In the genotyping results, 53.3 % of the strains clustered; of these strains, half were in two of the four clusters formed, and the patients had visited the hospital on the same day. Drug resistance was observed in 11.7 % of the strains.Conclusion. Putative nosocomial transmission of M. tuberculosis was detected, showing that genotyping is a powerful approach for understanding the dynamics of M. tuberculosis transmission, especially in a high-burden TB and HIV landscape.

Two tales: Worldwide distribution of Central Asian (CAS) versus ancestral East-African Indian (EAI) lineages of Mycobacterium tuberculosis underlines a remarkable cleavage for phylogeographical, epidemiological and demographical characteristics

The East African Indian (EAI) and Central Asian (CAS) lineages of Mycobacterium tuberculosis complex (MTBC) mainly infect tuberculosis (TB) patients in the eastern hemisphere which contains many of the 22 high TB burden countries including China and India. We investigated if phylogeographical, epidemiological and demographical characteristics for these 2 lineages differed in SITVIT2 database. Genotyping results and associated data (age, sex, HIV serology, drug resistance) on EAI and CAS lineages (n = 10,974 strains) were extracted. Phylogenetic and Bayesian, and other statistical analyses were used to compare isolates. The male/female sex ratio was 907/433 (2.09) for the EAI group vs. 881/544 (1.62) for CAS (p-value<0.002). The proportion of younger patients aged 0-20 yrs. with CAS lineage was significantly higher than for EAI lineage (18.07% vs. 10.85%, p-value<0.0001). The proportion of multidrug resistant and extensively drug resistant TB among CAS group (30.63% and 1.03%, respectively) was significantly higher than in the EAI group (12.14% and 0.29%, respectively; p-value<0.0001). Lastly, the proportion of HIV+ patients was 20.34% among the EAI group vs. 3.46% in the CAS group (p-value<0.0001). This remarkable split observed between various parameters for these 2 lineages was further corroborated by their geographic distribution profile (EAI being predominantly found in Eastern-Coast of Africa, South-India and Southeast Asia, while CAS was predominantly found in Afghanistan, Pakistan, North India, Nepal, Middle-east, Libya, Sudan, Ethiopia, Kenya and Tanzania). Some geo-specificities were highlighted. This study demonstrated a remarkable cleavage for aforementioned characteristics of EAI and CAS lineages, showing a North-South divide along the tropic of cancer in Eastern hemisphere-mainly in Asia, and partly prolonged along the horn of Africa. Such studies would be helpful to better comprehend prevailing TB epidemic in context of its historical spread and evolutionary features, and provide clues to better treatment and patient-care in countries and regions concerned by these lineages.

Comprehensive Analysis and Comparison on the Codon Usage Pattern of Whole Mycobacterium tuberculosis Coding Genome from Different Area

Phenomenon of unequal use of synonymous codons in Mycobacterium tuberculosis is common. Codon usage bias not only plays an important regulatory role at the level of gene expression, but also helps in improving the accuracy and efficiency of translation. Meanwhile, codon usage pattern of Mycobacterium tuberculosis genome is important for interpreting evolutionary characteristics in species. In order to investigate the codon usage pattern of the Mycobacterium tuberculosis genome, 12 Mycobacterium tuberculosis genomes from different area are downloaded from the GeneBank. The correlations between G₃, GC₁₂, whole GC content, codon adaptation index, codon bias index, and so on of Mycobacterium tuberculosis genomes are calculated. The ENC-plot, relationship between A₃/(A₃ + T₃) and G₃/(G₃ + C₃), GC₁₂ versus GC₃ plot, and the RSCU of overall/separated genomes all show that the codon usage bias exists in all 12 Mycobacterium tuberculosis genomes. Lastly, relationship between CBI and the equalization of ENC shows a strong negative correlation between them. The relationship between protein length and GC content (GC₃ and GC₁₂) shows that more obvious differences in the GC content may be in shorter protein. These results show that codon usage bias existing in the Mycobacterium tuberculosis genomes could be used for further study on their evolutionary phenomenon.