Discovery of a novel Mycobacterium tuberculosis lineage that is a major cause of tuberculosis in Rio de Janeiro, Brazil

Molecular epidemiology of Mycobacterium tuberculosis in Ecuador: Recent advances and future challenges

Tuberculosis (TB) is one of the three leading causes of death from a single infectious agent, Mycobacterium tuberculosis (MTB), together with COVID-19 and HIV/AIDS. This disease places a heavy burden on countries with low socio-economic development and aggravates existing inequalities. For the year 2021, estimations for Ecuador were 8500 TB cases, of which 370 were associated to multiple drug resistance (TB-MDR), and 1160 deaths. In the same year, Ecuador notified 5973 total cases, 401 of them were TB-MDR, pointing out an under diagnosis problem. The few molecular epidemiology studies available conclude that L4 is the most prevalent MTB lineage in Ecuador (with LAM as the main L4 sublineage), but L2-Beijing family is also present at low prevalence. Nevertheless, with less than 1 % MTB isolates genetically characterized by either MIRU-VNTR, spolygotyping or WGS to date, molecular epidemiology research must me improved to assist the TB surveillance and control program in Ecuador.

Prevalence, drug resistance, and genotypic diversity of the RDRio subfamily of Mycobacterium tuberculosis in Ecuador: a retrospective analysis for years 2012-2016

Introduction

A major sublineage within the Mycobacterium tuberculosis (MTB) LAM family characterized by a new in-frame fusion gene Rv3346c/55c was discovered in Rio de Janeiro (Brazil) in 2007, called RDRio, associated to drug resistance. The few studies about prevalence of MTB RDRio strains in Latin America reported values ranging from 3% in Chile to 69.8% in Venezuela, although no information is available for countries like Ecuador.

Methods

A total of 814 MTB isolates from years 2012 to 2016 were screened by multiplex PCR for RDRio identification, followed by 24-loci MIRU-VNTR and spoligotyping.

Results

A total number of 17 MTB RDRio strains were identified, representing an overall prevalence of 2.09% among MTB strains in Ecuador. While 10.9% of the MTB isolates included in the study were multidrug resistance (MDR), 29.4% (5/17) of the RDRio strains were MDR.

Discussion

This is the first report of the prevalence of MTB RDRio in Ecuador, where a strong association with MDR was found, but also a very low prevalence compared to other countries in Latin America. It is important to improve molecular epidemiology tools as a part of MTB surveillance programs in Latin America to track the transmission of potentially dangerous MTB stains associated to MDR TB like MTB RDRio.

Single-tube protocol for culture-independent spoligotyping of Mycobacterium tuberculosis based on MeltArray

IMPORTANCE

Spacer oligonucleotide typing (spoligotyping), the first-line genotyping assay for Mycobacterium tuberculosis (MTB), plays a fundamental role in the investigation of its epidemiology and evolution. In this study, we established a single-tube spoligotyping assay using MeltArray, a highly multiplex polymerase chain reaction (PCR) approach that runs on a real-time PCR thermocycler. The MeltArray protocol included an internal positive control, gyrB, to indicate the abundance of MTB via the quantification cycle and 43 spacers to identify the spoligotype via melting curve analysis. The entire protocol was completed in a single step within 2.5 hours. The lowest detectable copy number for the tested strains was 20 copies/reaction and thus sufficient for analyzing both culture and sputum samples. We conclude that MeltArray-based spoligotyping could be used immediately in low- and middle-income countries with a high tuberculosis burden, given its easy access, improved throughput, and potential applicability to clinical samples.

BORC complex specific components and Kinesin-1 mediate autophagy evasion by the autophagy-resistant Mycobacterium tuberculosis Beijing strain

Autophagy induction by starvation has been shown to enhance lysosomal delivery to mycobacterial phagosomes, resulting in the restriction of the Mycobacterium tuberculosis reference strain H37Rv. In contrast to H37Rv, our previous study showed that strains belonging to the notorious M. tuberculosis Beijing genotype could evade autophagic elimination. Our recent RNA-Seq analysis also discovered that the autophagy-resistant M. tuberculosis Beijing strain (BJN) evaded autophagic control by upregulating the expression of Kxd1, a BORC complex component, and Plekhm2, both of which function in lysosome positioning towards the cell periphery in host macrophages, thereby suppressing enhanced lysosomal delivery to its phagosome and sparing the BJN from elimination as a result. In this work, we further characterised the other specific components of the BORC complex, BORC5-8, and Kinesin proteins in autophagy resistance by the BJN. Depletion of BORCS5-8 and Kinesin-1, but not Kinesin-3, reverted autophagy avoidance by the BJN, resulting in increased lysosomal delivery to the BJN phagosomes. In addition, the augmented lysosome relocation towards the perinuclear region could now be observed in the BJN-infected host cells depleted in BORCS5-8 and Kinesin-1 expressions. Taken together, the data uncovered new roles for BORCS5-8 and Kinesin-1 in autophagy evasion by the BJN.

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Genetic Diversity and Transmission of Multidrug Resistant Mycobacterium tuberculosis strains in Lusaka, Zambia

OBJECTIVE

Zambia is among the 30 high tuberculosis burden countries in the world. Despite increasing reports of multidrug resistant tuberculosis (MDR-TB) in routine surveillance, information on the transmission of MDR Mycobacterium tuberculosis strains is largely unknown. This study elucidated genetic diversity and transmission of MDR M. tuberculosis strains in Lusaka, Zambia.

METHODS

Eighty-five MDR M. tuberculosis samples collected from the year 2013 to 2017 at the University Teaching Hospital were used. Drug-resistance associated gene sequencing, spoligotyping, 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats, and multiplex PCR for RD-Rio sub-lineage identification were applied.

RESULTS

Clades identified were LAM (48%), CAS (29%), T (14%), X (6%) and Harlem (2%). Strains belonging to SITs 21/CAS1-Kili and 20/LAM1 formed the largest clonal complexes. Combined spoligotyping and 24 loci-MIRU-VNTR revealed 47 genotypic patterns with clustering rate of 63%. Ninety five percent of LAM strains belonged to RD-Rio sub-lineage.

CONCLUSION

The high clustering rate suggested that a large proportion of MDR-TB was due to recent transmission rather than independent acquisition of MDR. This spread was attributed to clonal expansion of SIT21/CAS1-Kili and SIT20/LAM1 strains. Therefore, TB control programs recommending genotyping coupled with conventional epidemiological methods can guide measures for stopping the spread of MDR-TB.

A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

Mycobacterium tuberculosis complex (MTBC) species are classic examples of genetically monomorphic microorganisms due to their low genetic variability. Whole-genome sequencing made it possible to describe both the main species within the complex and M. tuberculosis lineages and sublineages. This differentiation is based on single nucleotide polymorphisms (SNPs) and large sequence polymorphisms in the so-called regions of difference (RDs). Although a number of studies have been performed to elucidate RD localizations, their distribution among MTBC species, and their role in the bacterial life cycle, there are some inconsistencies and ambiguities in the localization of RDs in different members of the complex. To address this issue, we conducted a thorough search for all possible deletions in the WGS data collection comprising 721 samples representing the full MTBC diversity. Discovered deletions were compared with a list of all previously described RDs. As with the SNP-based analysis, we confirmed the specificities of 79 regions at the species, lineage, or sublineage level, 17 of which are described for the first time. We also present RDscan (https://github.com/dbespiatykh/RDscan), an open-source workflow, which detects deletions from short-read sequencing data and correlates the results with high-specificity RDs, curated in this study. Testing of the workflow on a collection comprising ∼7,000 samples showed a high specificity of the found RDs. This study provides novel details that can contribute to a better understanding of the species differentiation within the MTBC and can help to determine how individual clusters evolve within various MTBC species.

IMPORTANCE Reductive genome evolution is one of the most important and intriguing adaptation strategies of different living organisms to their environment. Mycobacterium offers several notorious examples of either naturally reduced (Mycobacterium leprae) or laboratory-reduced (Mycobacterium bovis BCG) genomes. Mycobacterium tuberculosis complex has its phylogeny unambiguously framed by large sequence polymorphisms that present unidirectional unique event changes. In the present study, we curated all known regions of difference and analyzed both Mycobacterium tuberculosis and animal-adapted MTBC species. For 79 loci, we have shown a relationship with phylogenetic units, which can serve as a marker for diagnosing or studying biological effects. Moreover, intersections were found for some loci, which may indicate the nonrandomness of these processes and the involvement of these regions in the adaptation of bacteria to external conditions.

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.

Lysosome repositioning as an autophagy escape mechanism by Mycobacterium tuberculosis Beijing strain

Induction of host cell autophagy by starvation was shown to enhance lysosomal delivery to mycobacterial phagosomes, resulting in the restriction of Mycobacterium tuberculosis reference strain H37Rv. Our previous study showed that strains belonging to M. tuberculosis Beijing genotype resisted starvation-induced autophagic elimination but the factors involved remained unclear. Here, we conducted RNA-Seq of macrophages infected with the autophagy-resistant Beijing strain (BJN) compared to macrophages infected with H37Rv upon autophagy induction by starvation. Results identified several genes uniquely upregulated in BJN-infected macrophages but not in H37Rv-infected cells, including those encoding Kxd1 and Plekhm2, which function in lysosome positioning towards the cell periphery. Unlike H37Rv, BJN suppressed enhanced lysosome positioning towards the perinuclear region and lysosomal delivery to its phagosome upon autophagy induction by starvation, while depletion of Kxd1 and Plekhm2 reverted such effects, resulting in restriction of BJN intracellular survival upon autophagy induction by starvation. Taken together, these data indicated that Kxd1 and Plekhm2 are important for the BJN strain to suppress lysosome positioning towards the perinuclear region and lysosomal delivery into its phagosome during autophagy induction by starvation to evade starvation-induced autophagic restriction.

Genomic analysis of Latin American-Mediterranean family of Mycobacterium tuberculosis clinical strains from Kazakhstan

The human-adapted strains of the Mycobacterium tuberculosis complex (MTBC) comprise seven phylogenetic lineages originally associated with their geographical distribution. Here, we report the genomes of three drug-resistant clinical isolates of the Latin American-Mediterranean (LAM) family collected in Kazakhstan. We utilised whole-genome sequencing to study the distribution and drug resistance of these isolates. Phylogenetic analysis grouped the genomes described in this study with the sequences from Russia, Uzbekistan, and Kazakhstan belonging to the LAM family. One isolate has acquired extensive drug resistance to seven antituberculosis drugs. Our results suggest at least two multi-drug resistant (MDR)/extensively drug-resistant (XDR)-associated genotypes of the LAM family circulate in Kazakhstan.

Most LAM Mycobacterium tuberculosis strains in Venezuela, but not SIT605, belong to the RDRio subfamily

Tuberculosis is a global public health problem that is resurgent in Venezuela, with 13 thousand estimated new cases in 2018. Strains of the Mycobacterium tuberculosis RD^Rio, subfamily belong to the Latín American Mediterranean (LAM) family and are a major cause of TB in Rio de Janeiro, Brazil. LAM strains predominate in Venezuela, where spoligotype SIT605 is common, but surprisingly not found elsewhere. We sought to assess the presence of RD^Rio strains in tuberculosis patients in different regions of Venezuela and determine whether SIT605 also belongs to the RD^Rio subfamily. Using spoligotyping and MIRU-VNTR 24 loci, we identified 86 clinical LAM and SIT605 isolates from the Venezuelan capital Caracas and several Venezuelan states. Region of difference deletion loci RD174 and RD^Rio, and also IS1561 were used to identify strains of the RD^Rio subfamily, while IS6110 at position 932,204 and the Ag85C¹⁰³ polymorphism were used to validate SIT 605 as a LAM family strain. We found that 69.8% of the isolates were RD^Río, including 94.3% of strains isolated in Caracas, 17.9% isolated in the state of Carabobo, the two strains analyzed from Delta Amacuro, and one each from Sucre, Apure and Aragua states. RD^Rio was in 100% of: SIT17 (LAM 2); SIT20 (LAM 1); SITs 93, 1694, 1696, 960, 1367 (LAM 5); and SITs 216 (LAM 9); but only 75% of SIT42 (LAM 9) strains. Thus, most of the LAM strains in Venezuela belong to the RD^Río subfamily. SIT 605 strains, although LAM, are not in the RD^Río subfamily.

Identification of gene fusion events in Mycobacterium tuberculosis that encode chimeric proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen responsible for causing tuberculosis. The harsh environment in which M. tuberculosis survives requires this pathogen to continuously adapt in order to maintain an evolutionary advantage. However, the apparent absence of horizontal gene transfer in M. tuberculosis imposes restrictions in the ways by which evolution can occur. Large-scale changes in the genome can be introduced through genome reduction, recombination events and structural variation. Here, we identify a functional chimeric protein in the ppe38-71 locus, the absence of which is known to have an impact on protein secretion and virulence. To examine whether this approach was used more often by this pathogen, we further develop software that detects potential gene fusion events from multigene deletions using whole genome sequencing data. With this software we could identify a number of other putative gene fusion events within the genomes of M. tuberculosis isolates. We were able to demonstrate the expression of one of these gene fusions at the protein level using mass spectrometry. Therefore, gene fusions may provide an additional means of evolution for M. tuberculosis in its natural environment whereby novel chimeric proteins and functions can arise.

Global Distribution and Evolution of Mycobacterium bovis Lineages

Mycobacterium bovis is the main causative agent of zoonotic tuberculosis in humans and frequently devastates livestock and wildlife worldwide. Previous studies suggested the existence of genetic groups of M. bovis strains based on limited DNA markers (a.k.a. clonal complexes), and the evolution and ecology of this pathogen has been only marginally explored at the global level. We have screened over 2,600 publicly available M. bovis genomes and newly sequenced four wildlife M. bovis strains, gathering 1,969 genomes from 23 countries and at least 24 host species, including humans, to complete a phylogenomic analyses. We propose the existence of four distinct global lineages of M. bovis (Lb1, Lb2, Lb3, and Lb4) underlying the current disease distribution. These lineages are not fully represented by clonal complexes and are dispersed based on geographic location rather than host species. Our data divergence analysis agreed with previous studies reporting independent archeological data of ancient M. bovis (South Siberian infected skeletons at ∼2,000 years before present) and indicates that extant M. bovis originated between 715 and 3,556 years BP, with later emergence in the New World and Oceania, likely influenced by trades among countries.

Molecular characterisation of multidrug-resistant Mycobacterium tuberculosis isolates from a high-burden tuberculosis state in Brazil

Tuberculosis (TB) is the leading cause of death among infectious diseases worldwide. Among the estimated cases of drug-resistant TB, approximately 60% occur in the BRICS countries (Brazil, Russia, India, China and South Africa). Among Brazilian states, primary and acquired multidrug-resistant TB (MDR-TB) rates were the highest in Rio Grande do Sul (RS). This study aimed to perform molecular characterisation of MDR-TB in the State of RS, a high-burden Brazilian state. We performed molecular characterisation of MDR-TB cases in RS, defined by drug susceptibility testing, using 131 Mycobacterium tuberculosis (M.tb) DNA samples from the Central Laboratory. We carried out MIRU-VNTR 24loci, spoligotyping, sequencing of the katG, inhA and rpoB genes and RD^Rio sublineage identification. The most frequent families found were LAM (65.6%) and Haarlem (22.1%). RD^Rio deletion was observed in 42 (32%) of the M.tb isolates. Among MDR-TB cases, eight (6.1%) did not present mutations in the studied genes. In 116 (88.5%) M.tb isolates, we found mutations associated with rifampicin (RIF) resistance in rpoB gene, and in 112 isolates (85.5%), we observed mutations related to isoniazid resistance in katG and inhA genes. An insertion of 12 nucleotides (CCAGAACAACCC) at the 516 codon in the rpoB gene, possibly responsible for a decreased interaction of RIF and RNA polymerase, was found in 19/131 of the isolates, belonging mostly to LAM and Haarlem families. These results enable a better understanding of the dynamics of transmission and evolution of MDR-TB in the region.

Circulation of M. tuberculosis Beijing genotype in Latin America and the Caribbean

Lineage 2 (East Asian), which includes the Beijing genotype, is one of the most prevalent lineages of Mycobacterium tuberculosis (Mtb) throughout the world. The Beijing family is associated to hypervirulence and drug-resistant tuberculosis. The study of this genotype's circulation in Latin America is crucial for achieving total control of TB, the goal established by the World Health Organization, for the American sub-continent, before 2035. In this sense, the present work presents an overview of the status of the Beijing genotype for this region, with a bibliographical review, and data analysis of MIRU-VNTRs for available Beijing isolates. Certain countries present a prevalent trend of <5%, suggesting low transmissibility for the region, with the exception of Cuba (17.2%), Perú (16%) and Colombia (5%). Minimum Spanning Tree analysis, obtained from MIRU-VNTR data, shows distribution of specific clonal complex strains in each country. From this data, in most countries, we found that molecular epidemiology has not been a tool used for the control of TB, suggesting that the Beijing genotype may be underestimated in Latin America. It is recommended that countries with the highest incidence of the Beijing genotype use effective control strategies and increased care, as a requirement for public health systems.

Usefulness of 3'- 5' IS6110-RFLP genotyping and spoligotyping of Mycobacterium tuberculosis isolated in a tertiary hospital: a retrospective study detecting unsuspected epidemiological events

A drug resistance survey involving Mycobacterium tuberculosis isolated from patients of a tertiary Hospital in the Rio de Janeiro city (RJ), Brazil, between the years 1996 and 1998 revealed a high frequency of isoniazid (H^R) resistance. These isolates were revisited and genotyped. Patients came from different RJ neighborhoods and municipalities, and 70% were outpatients. Applying the 3’ and 5’ IS 6110 -RFLP and the Spoligotype genotyping methods, the clonal structure of this population was investigated obtaining a snapshot of past epidemiological events. The 3’ clusters were subsequently 5’ IS 6110 -RFLP typed. Spoligotyping was analyzed in the SITVIT2 database. Epidemiological relationships were investigated. The major lineage was T (54.4%), and SIT 53/T1 and SIT 535/T1 were the most frequent. The T1 sublineage comprises 12.8% of resistant strains and SIT 535 were assigned for 31.8% of them. Orphan patterns corresponded to 12% and 73.3% and belonged to the T lineage. One pattern was unlisted in the SITVIT2. The 5’ IS 6110 -RFLP did not confirm 3/12 of the 3’ IS 6110 -RFLP clusters. A combination of all methods decreased the number of clusters to three. Nosocomial transmission was associated with one cluster involving a hospital cupbearer. This event was suspected in a multidrug resistant-TB inpatient caregiver who harbored a mixed infection. The 3’ IS 6110 clusters were associated with H^R (p=0.046). These genotypic retrospective data may reflect a fraction of more extensive recent transmission in different communities that may be corroborated by the concentration of H^R patients, and may serve as a database for further evolutionary and characterization evaluation of circulating strains and together with epidemiological data favors a more effective transmission control.

Mycobacterium tuberculosis RD-Rio Strain in Kazakhstan

Mycobacterium tuberculosis RD-Rio strains are still rare in the former Soviet Union countries and Asia. We describe a strain in Kazakhstan that belongs to the RD-Rio secondary branch, which is endemic to northwest Russia and eastern Europe. Although RD-Rio strains are frequently multidrug resistant, this heterogeneous branch included only drug-susceptible isolates.

Frequency of the Mycobacterium tuberculosis RDRio genotype and its association with multidrug-resistant tuberculosis

BACKGROUND

In recent decades, Mycobacterium tuberculosis with the RD^Rio genotype, frequently isolated from tuberculosis patients in Rio de Janeiro, has become part of the Latin American - Mediterranean (LAM) family and has been associated with multidrug-resistant tuberculosis (MDR-TB). The aim of this study was to investigate the frequency of M. tuberculosis RD^Rio in the state of Minas Gerais, Brazil, and its relationship with MDR-TB.

METHODS

For convenience, 172 susceptible and 63 MDR M. tuberculosis isolates were taken from pulmonary samples from patients diagnosed between January 2007 and December 2011. The DNA extracted from these isolates was analyzed by spoligotyping, PCR-RFLP to characterize fbpC¹⁰³/Ag85C103, multiplex PCR to detect RD^Rio and RD174, and MIRU-VNTR 24 loci.

RESULTS

Among the 235 isolates, the RD^Rio pattern was identified in 122 (51.9%) isolates (IC 0.45-0.58), with 100 (42.5%) wild-type and 13 (5.5%) mixed pattern isolates, whereas RD174 was identified in 93 of the 122 RD^Rio positive samples (76.3%). The LAM family and the LAM9 lineage were the most frequently identified among the isolates in this study. Among the 63 MDR isolates, 41 (65.1%) were RD^Rio and 28 (44.4%) RD174.

CONCLUSION

The association of both deletions with MDR proved to be statistically significant, corroborating the few reports that have associated RD^Rio with MDR.

Genetic diversity and drug susceptibility profiles of Mycobacterium tuberculosis obtained from Saint Peter's TB specialized Hospital, Ethiopia

Background

Tuberculosis (TB) is one of the major public health problems in Ethiopia. Data on genetic diversity and resistance profile of circulating TB strains is critical for informing the national TB control program.

Methods

A cross-sectional study was conducted on 213 smear positive pulmonary TB patients between 2015 and 2016. Sputum samples were cultured on LJ media following the Petroff’s method. Region of difference-9 (RD9)-deletion typing and spoligo-typing were performed for molecular analysis of M. tuberculosis at species and strain levels, respectively. Drug sensitivity and mutation patterns of the isolates were assessed by the conventional indirect proportion method and molecular line probe assays (LPAs), respectively. Data were analyzed using statistical package for social sciences (SPSS) software version 20.

Results

Spoligo-typing of 150 M. tuberculosis isolates led to 57 different patterns of which 25 were new strains. The majority (71.6%) of the isolates were grouped in to 17 clusters consisting 2 to 24 isolates. The majority of the strains belonged to Euro-American lineage and the predominant spoligotypes were SIT 37 and SIT 149. MDR-TB was detected in 5.2% and 20.3% of new and retreatment cases, respectively. Two MDR-TB isolates exhibited additional resistance to one of the second line anti-TB drugs. Common gene mutations including S531L, S315T1 and M306V were detected in RIF, INH and EMB resistant strains, respectively.

Conclusions

The identification of several new strains, higher proportion of MDR-TB and higher clustering rate in this study, warrants the need for re-enforcement of the national TB control program. The detection of common gene mutations in the majority drug resistant strains might suggest the feasibility of LPAs for rapid screening of drug resistant M. tuberculosis strains in Ethiopia.

Exploring the "Latin American Mediterranean" family and the RDRio lineage in Mycobacterium tuberculosis isolates from Paraguay, Argentina and Venezuela

Background

The Latin American & Mediterranean (LAM) spoligotype family is one of the most successful genotype of Mycobacterium tuberculosis worldwide and particularly prevalent in South-America. Within this family, a sublineage named Region of Difference Rio (RD^Rio) was reported initially in Brazil and is characterized by a genomic deletion of about 26.3 kb. This lineage seems to show a specific adaptation to the Euro-Latin American population. In this context, we sought to evaluate the LAM family and the presence of the RD^Rio genotype in samples from three Latin American countries including Paraguay, Venezuela and Argentina. To detect LAM strains reliably we applied a typing scheme using spoligotyping, 12 loci MIRU-VNTR, the Ag85C¹⁰³ SNP and the regions of difference RD^Rio and RD174. IS6110-RFLP results were also used when available.

Results

Genotyping of 413 M. tuberculosis isolates from three Latin-American countries detected LAM (46%) and the ill-defined T clade (16%) as the most frequent families. The highest clustering rate was detected in the sample population from the city of Caracas in Venezuela. We observed considerable differences in the presence of the RD^Rio lineage, with high frequency in Caracas-Venezuela (55%) and low frequency in Buenos Aires-Argentina (11%) and Paraguay (10%). The molecular markers (RD174, Ag85C¹⁰³, MIRU02-MIRU40 signature) of the RD^Rio lineage were essentially confirmed. For the LAM family, the most polymorphic loci were MIRU40, MIRU31, MIRU10, MIRU26, MIRU16 and the least polymorphic MIRU24, MIRU20, MIRU04, MIRU23.

Conclusions

Our results suggest a differential adaptation of LAM-sublineages in neighboring populations and that RD^Rio strains spread regionally with different rates of distribution. The Ag85C SNP and RDs (RD174, RD^Rio) tested in this study can in fact facilitate molecular epidemiological studies of LAM strains in endemic settings and low-income countries.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1479-6) contains supplementary material, which is available to authorized users.

Tuberculosis Genetic Epidemiology: A Latin American Perspective

There are an estimated 10 million new cases of tuberculosis worldwide annually, with 282,000 new or relapsed cases each year reported from the Americas. With improvements in genome sequencing technology, it is now possible to study the genetic diversity of tuberculosis with much greater resolution. Although tuberculosis bacteria do not engage in horizontal gene transfer, the genome is far more variable than previously thought. The study of genome-wide variation in tuberculosis has improved our understanding of the evolutionary origins of tuberculosis, the arrival of tuberculosis in Latin America, the genetic determinants of drug resistance, and lineage-specific associations with important clinical phenotypes. This article reviews what is known about the arrival of tuberculosis in Latin America, the genetic diversity of tuberculosis in Latin America, and the genotypic determinants of clinical phenotypes.

Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Genotypic diversity of Mycobacterium tuberculosis in Buenos Aires, Argentina

Buenos Aires is an overpopulated port city historically inhabited by people of European descent. Together with its broader metropolitan area, the city exhibits medium tuberculosis rates, and receives migrants, mainly from tuberculosis highly endemic areas of Argentina and neighboring countries. This work was aimed to gain insight into the Mycobacterium tuberculosis population structure in two suburban districts of Buenos Aires which are illustrative of the overall situation of tuberculosis in Argentina. The Lineage 4 Euro-American accounted for >99% of the 816 isolates analyzed (one per patient). Frequencies of spoligotype families were T 35.9%, LAM 33.2%, Haarlem 19.5%, S 3.2%, X 1.5%, Ural 0.7%, BOV 0.2%, Beijing 0.2%, and Cameroon 0.2%. Unknown signatures accounted for 5.3% isolates. Of 55 spoligotypes not matching any extant shared international type (SIT) in SITVIT database, 22 fitted into 15 newly-issued SITs. Certain autochthonous South American genotypes were found to be actively evolving. LAM3, which is wild type for RD^rio, was the predominant LAM subfamily in both districts and the RD^rio signature was rare among autochthonous, newly created, SITs and orphan patterns. Two genotypes that are rarely observed in neighboring countries ̶ SIT2/H2 and SIT159/T1 Tuscany ̶ were conspicuously represented in Argentina. The infrequent Beijing patterns belonged to Peruvian patients. We conclude that the genotype diversity observed reflects the influence of the Hispanic colonization and more recent immigration waves from Mediterranean and neighboring countries. Unlike in Brazil, the RD^rio type does not play a major role in the tuberculosis epidemic in Buenos Aires.

RDRioMycobacterium tuberculosis lineage in the Brazil/Paraguay/Argentina triple border

The high tuberculosis (TB) incidence rates, the closeness of the cities and the high migration flux on the Brazil/Paraguay/Argentina border deserves an in-depth study, using Mycobacterial Interspersed Repetitive Unit (MIRU) and Spoligotyping genetic markers to explore the impact of the Mycobacterium tuberculosis RD^Rio lineage on disease transmission and resistance to anti-TB drugs in this setting. Although without the totality of M. tuberculosis isolates causing TB in this studied setting, a number of 97 isolates obtained from sputa samples culture of patients with confirmed TB, from 2013 to 2015, were submitted to 24 loci MIRU, Spoligotyping, detection of RD^Rio lineage and detection of mutation related to isoniazid and rifampicin resistance by MTBDRplus/DNA STRIP. In this sample, it was observed high clonal variability of circulating M. tuberculosis isolates causing TB in Brazilian cities bordering Paraguay and Argentina. The percentage of RD^Rio lineage causing TB in this setting was 15.46%, and lower than the detected in different areas of Brazil. According to 24 loci MIRU, the major MIRU International Type (MIT) related with RD^Rio lineage were MIT 26, MIT 738, MIT 601 with four, two and one isolates, respectively. Eight isolates with RD^Rio marker were classified as orphans. The mainly Spoligofamily related with RD^Rio lineage was LAM1 and LAM9 and no relationship between RD^Rio lineage and resistance in M. tuberculosis isolates circulating in this setting could be established. This work is pioneer in studying the dynamics of RD^Rio lineage transmission on the Brazil/Paraguay/Argentina border and deserves further studies to analyze the real contribution of the RD^Rio lineage in outbreaks and the risk of significant development of MDR-TB in the setting studied.

Genetic diversity of drug and multidrug-resistant Mycobacterium tuberculosis circulating in Veracruz, Mexico

Background

Mexico is one of the most important contributors of drug and multidrug-resistant tuberculosis in Latin America; however, knowledge of the genetic diversity of drug-resistant tuberculosis isolates is limited.

Methods

In this study, the genetic structure of 112 Mycobacterium tuberculosis strains from the southeastern Mexico was determined by spoligotyping and 24-loci MIRU-VNTRs.

Findings

The results show eight major lineages, the most of which was T1 (24%), followed by LAM (16%) and H (15%). A total of 29 (25%) isolates were identified as orphan. The most abundant SITs were SIT53/T1 and SIT42/LAM9 with 10 isolates each and SIT50/H3 with eight isolates. Fifty-two spoligotype patterns, twenty-seven clusters and ten clonal complexes were observed, demonstrating an important genetic diversity of drug and multidrug-resistant tuberculosis isolates in circulation and transmission level of these aggravated forms of tuberculosis. Being defined as orphan or as part of an orphan cluster, was a risk factor for multidrug resistant-tuberculosis (OR 2.5, IC 1.05–5.86 and OR 3.3, IC 1–11.03, respectively). Multiple correspondence analyses showed association of some clusters and SITs with specific geographical locations.

Conclusions

Our study provides one of the most detailed description of the genetic structure of drug and multidrug-resistant tuberculosis strains in southeast Mexico, establishing for the first time a baseline of the genotypes observed in resistant isolates circulating, however further studies are required to better elucidate the genetic structure of tuberculosis in region and the factors that could be participating in their dispersion.

Risk factors associated with cluster size of Mycobacterium tuberculosis (Mtb) of different RFLP lineages in Brazil

BACKGROUND

Tuberculosis (TB) transmission is influenced by patient-related risk, environment and bacteriological factors. We determined the risk factors associated with cluster size of IS6110 RFLP based genotypes of Mycobacterium tuberculosis (Mtb) isolates from Vitoria, Espirito Santo, Brazil.

METHODS

Cross-sectional study of new TB cases identified in the metropolitan area of Vitoria, Brazil between 2000 and 2010. Mtb isolates were genotyped by the IS6110 RFLP, spoligotyping and RDRio. The isolates were classified according to genotype cluster sizes by three genotyping methods and associated patient epidemiologic characteristics. Regression Model was performed to identify factors associated with cluster size.

RESULTS

Among 959 Mtb isolates, 461 (48%) cases had an isolate that belonged to an RFLP cluster, and six clusters with ten or more isolates were identified. Of the isolates spoligotyped, 448 (52%) were classified as LAM and 412 (48%) as non-LAM. Our regression model found that 6-9 isolates/RFLP cluster were more likely belong to the LAM family, having the RDRio genotype and to be smear-positive (adjusted OR = 1.17, 95% CI 1.08-1.26; adjusted OR = 1.25, 95% CI 1.14-1.37; crude OR = 2.68, 95% IC 1.13-6.34; respectively) and living in a Serra city neighborhood decrease the risk of being in the 6-9 isolates/RFLP cluster (adjusted OR = 0.29, 95% CI, 0.10-0.84), than in the others groups. Individuals aged 21 to 30, 31 to 40 and > 50 years were less likely of belonging the 2-5 isolates/RFLP cluster than unique patterns compared to individuals < 20 years of age (adjusted OR = 0.49, 95% CI 0.28-0.85, OR = 0.43 95% CI 0.24-0.77and OR = 0. 49, 95% CI 0.26-0.91), respectively. The extrapulmonary disease was less likely to occur in those infected with strains in the 2-5 isolates/cluster group (adjustment OR = 0.45, 95% CI 0.24-0.85) than unique patterns.

CONCLUSIONS

We found that a large proportion of new TB infections in Vitoria is caused by prevalent Mtb genotypes belonging to the LAM family and RDRio genotypes. Such information demonstrates that some genotypes are more likely to cause recent transmission. Targeting interventions such as screening in specific areas and social risk groups, should be a priority for reducing transmission.

Molecular epidemiology of Mycobacterium tuberculosis strains from prison populations in Santa Catarina, Southern Brazil

The Tuberculosis (TB) notification rates are 5 to 81 times higher in prisons worldwide when compared to the general population. The state of Santa Catarina (SC) has few epidemiological data regarding TB in prisons. The aim of this study was to evaluate the molecular epidemiology of circulating strains in prisons of SC. The study comprised 95 clinical samples from six prisons. Among the cases included, all subjects were male, predominantly caucasians, and young adults, with low education level. The positive smear in the TB diagnosis comprised 62.0% of cases. About 50% of subjects had some condition associated with TB. The Spoligotyping results showed that the most frequent lineages were LAM (50.7%), T (22.2%) and S (11.6%). The 12-loci MIRU generated 62 different genotypes. The MSTs showed evolutionary relationships between Mycobacterium tuberculosis spoligotypes from SC and evolutionary relationships between the prison isolates and studied parameters. This first study on TB in prison units of SC highlighted the predominance of SIT216/LAM5, and SIT34/S. Interestingly, his profile was found to be different from that observed in a previous study performed with the state's general population. This data shows the need for continued surveillance of episodes of TB occurring among prison inmates in an emerging country like Brazil.

Proteogenomic Investigation of Strain Variation in Clinical Mycobacterium tuberculosis Isolates

Mycobacterium tuberculosis consists of a large number of different strains that display unique virulence characteristics. Whole-genome sequencing has revealed substantial genetic diversity among clinical M. tuberculosis isolates, and elucidating the phenotypic variation encoded by this genetic diversity will be of the utmost importance to fully understand M. tuberculosis biology and pathogenicity. In this study, we integrated whole-genome sequencing and mass spectrometry (GeLC-MS/MS) to reveal strain-specific characteristics in the proteomes of two clinical M. tuberculosis Latin American-Mediterranean isolates. Using this approach, we identified 59 peptides containing single amino acid variants, which covered ∼9% of all coding nonsynonymous single nucleotide variants detected by whole-genome sequencing. Furthermore, we identified 29 distinct peptides that mapped to a hypothetical protein not present in the M. tuberculosis H37Rv reference proteome. Here, we provide evidence for the expression of this protein in the clinical M. tuberculosis SAWC3651 isolate. The strain-specific databases enabled confirmation of genomic differences (i.e., large genomic regions of difference and nonsynonymous single nucleotide variants) in these two clinical M. tuberculosis isolates and allowed strain differentiation at the proteome level. Our results contribute to the growing field of clinical microbial proteogenomics and can improve our understanding of phenotypic variation in clinical M. tuberculosis isolates.

Evolution of Mycobacterium tuberculosis: New Insights into Pathogenicity and Drug Resistance

The tuberculosis agent Mycobacterium tuberculosis has undergone a long and selective evolution toward human infection and represents one of the most widely spread pathogens due to its efficient aerosol-mediated human-to-human transmission. With the availability of more and more genome sequences, the evolutionary trajectory of this obligate pathogen becomes visible, which provides us with new insights into the molecular events governing evolution of the bacterium and its ability to accumulate drug-resistance mutations. In this review, we summarize recent developments in mycobacterial research related to this matter that are important for a better understanding of the current situation and future trends and developments in the global epidemiology of tuberculosis, as well as for possible public health intervention possibilities.

Relation of Mycobacterium tuberculosis mutations at katG315 and inhA-15 with drug resistance profile, genetic background, and clustering in Argentina

We analyzed 362 isoniazid-resistant clinical isolates of Mycobacterium tuberculosis obtained countrywide for the presence of mutation at katG315 and inhA-15 in relation to genotype, pattern of phenotypic resistance to other drugs, and ability to spread. We found the following mutation frequencies: katG315MUT/inhA-15wt 53.0%, katG315wt/inhA-15MUT 27.4%, katG315wt/inhA-15wt 19.3%, and katG315MUT/inhA-15MUT only 0.3%. Mutation at katG315 associated with the LAM superfamily; mutation at inhA-15 associated with the S family and the T1 Tuscany genotype; the combination katG315wt/inhA-15wt associated with the T1 Ghana genotype. Isolates harboring katG315MUT/inhA-15wt tended to accumulate resistance to other drugs and were more frequently found in cluster; isolates harboring katG315wt/inhA-15wt were more frequently found as orphan isolates. Although epidemiological and host factors could also be modulating the events observed, in Argentina, the systematic genotyping of drug resistant clinical isolates could help to predict an enhanced risk of transmission and a propensity to develop resistance to increasing numbers of drugs.

A systematic review of East African-Indian family of Mycobacterium tuberculosis in Brazil

Introduction

The Mycobacterium tuberculosis East African-Indian (EAI) spoligotyping family (belonging to lineage 1, Indo-Oceanic, defined by the region of deletion RD239) is distributed worldwide, but is more prevalent in Southeast Asia, India, and East Africa. Studies in Latin America have rarely identified EAI. In this study, we describe the occurrence of the EAI family in Brazil.

Methods

EAI was identified in a systematic literature review of genetic diversity studies pertaining to M. tuberculosis in Brazil, as well as in a survey conducted in Salvador, Bahia, located in the northeastern region of this country.

Results

The EAI6-BGD1 spoligotyping family and the EAI5 Spoligotype International Type (SIT) 1983 clade were the most frequently reported, with wide distribution of this particular clade described in Brazil. The distribution of other EAI spoligotyping patterns with broader worldwide distribution was restricted to the southeastern region of the country.

Conclusions

EAI may be endemic at a low frequency in Brazil, with some clades indicating increased fitness with respect to this population.

The Evolution of Strain Typing in the Mycobacterium tuberculosis Complex

Tuberculosis (TB) is a contagious disease with a complex epidemiology. Therefore, molecular typing (genotyping) of Mycobacterium tuberculosis complex (MTBC) strains is of primary importance to effectively guide outbreak investigations, define transmission dynamics and assist global epidemiological surveillance of the disease. Large-scale genotyping is also needed to get better insights into the biological diversity and the evolution of the pathogen. Thanks to its shorter turnaround and simple numerical nomenclature system, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) typing, based on 24 standardized plus 4 hypervariable loci, optionally combined with spoligotyping, has replaced IS6110 DNA fingerprinting over the last decade as a gold standard among classical strain typing methods for many applications. With the continuous progress and decreasing costs of next-generation sequencing (NGS) technologies, typing based on whole genome sequencing (WGS) is now increasingly performed for near complete exploitation of the available genetic information. However, some important challenges remain such as the lack of standardization of WGS analysis pipelines, the need of databases for sharing WGS data at a global level, and a better understanding of the relevant genomic distances for defining clusters of recent TB transmission in different epidemiological contexts. This chapter provides an overview of the evolution of genotyping methods over the last three decades, which culminated with the development of WGS-based methods. It addresses the relative advantages and limitations of these techniques, indicates current challenges and potential directions for facilitating standardization of WGS-based typing, and provides suggestions on what method to use depending on the specific research question.

Mycobacterium tuberculosis lineage 4 comprises globally distributed and geographically restricted sublineages

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration.

Multidrug-Resistant Mycobacterium tuberculosis of the Latin American Mediterranean Lineage, Wrongly Identified as Mycobacterium pinnipedii (Spoligotype International Type 863 [SIT863]), Causing Active Tuberculosis in South Brazil

We recently detected the spoligotype patterns of strains of Mycobacterium pinnipedii, a species of the Mycobacterium tuberculosis complex, in sputum samples from nine cases with pulmonary tuberculosis residing in Porto Alegre, South Brazil. Because this species is rarely encountered in humans, we further characterized these nine isolates by additional genotyping techniques, including 24-locus mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing, verification of the loci TbD1, RD9, pks15/1, RD(Rio), and fbpC, the insertion of IS6110 at a site specific to the M. tuberculosis Latin American Mediterranean (LAM) lineage, and whole-genome sequencing. The combined analysis of these markers revealed that the isolates are in fact M. tuberculosis and more specifically belong to the LAM genotype. Most of these isolates (n8) were shown to be multidrug resistant (MDR), which prompted us to perform partial sequencing of the rpoA, rpoB, rpoC, katG, and inhA genes. Seven isolates (77.8%) carried the S315T mutation in katG, and one of these (11%) also presented the C((-17)T single-nucleotide polymorphism (SNP) in inhA. Interestingly, six of the MDR isolates also presented an undescribed insertion of 12 nucleotides (CCA GAA CAA CCC) in codon 516 of rpoB. No putative compensatory mutation was found in either rpoA or rpoC. This is the first report of an M. tuberculosis LAM family strain with a convergent M. pinnipedii spoligotype. These spoligotypes are observed in genotype databases at a modest frequency, highlighting that care must be taken when identifying isolates in the M. tuberculosis complex on the basis of single genetic markers.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

First insight into the molecular epidemiology of Mycobacterium tuberculosis in Santa Catarina, southern Brazil

Molecular epidemiology of Mycobacterium tuberculosis is useful for understanding disease transmission dynamics, and to establish strategic measures for TB control and prevention. The aim of this study was to analyze clinical, epidemiological and molecular characteristics of MTBC clinical isolates from Santa Catarina state, southern Brazil. During one-year period, 406 clinical isolates of MTBC were collected from Central Laboratory of Public Health and typed by spoligotyping. Demographic and clinical data were collected from the Brazilian National Mandatory Disease Reporting System. The majority of cases occurred in highest population densities regions and about 50% had some condition associated with TB. Among all isolates, 5.7% were MDR, which showed association with drug addiction. LAM was the most predominant lineage with 47.5%, followed by the T superfamily with 25.9% and Haarlem with 12.3%. The MST showed two major groups: the first was formed mainly by the LAM lineage and the second was mainly formed by the T and Haarlem lineages. Others lineages were distributed in peripheral positions. This study provides the first insight into the population structure of M. tuberculosis in SC State. Spoligotyping and other genotyping analyses are important to establish strategic measures for TB control and prevention.

Determinants of the Sympatric Host-Pathogen Relationship in Tuberculosis

Major contributions from pathogen genome analysis and host genetics have equated the possibility of Mycobacterium tuberculosis co-evolution with its human host leading to more stable sympatric host–pathogen relationships. However, the attribution to either sympatric or allopatric categories depends on the resolution or grain of genotypic characterization. We explored the influence on the sympatric host-pathogen relationship of clinical (HIV infection and multidrug-resistant tuberculosis [MDRTB]) and demographic (gender and age) factors in regards to the genotypic grain by using spacer oligonucleotide typing (spoligotyping) for classification of M. tuberculosis strains within the Euro-American lineage. We analyzed a total of 547 tuberculosis (TB) cases, from six year consecutive sampling in a setting with high TB-HIV coinfection (32.0%). Of these, 62.0% were caused by major circulating pathogen genotypes. The sympatric relationship was defined according to spoligotype in comparison to the international spoligotype database SpolDB4. While no significant association with Euro-American lineage was observed with any of the factors analyzed, increasing the resolution with spoligotyping evidenced a significant association of MDRTB with sympatric strains, regardless of the HIV status. Furthermore, distribution curves of the prevalence of sympatric and allopatric TB in relation to patients’ age showed an accentuation of the relevance of the age of onset in the allopatric relationship, as reflected in the trimodal distribution. On the contrary, sympatric TB was characterized by the tendency towards a typical (standard) distribution curve. Our results suggest that within the Euro-American lineage a greater degree of genotyping fine-tuning is necessary in modeling the biological processes behind the host-pathogen interplay. Furthermore, prevalence distribution of sympatric TB to age was suggestive of host genetic determinisms driven by more common variants.

Genetic Structuration, Demography and Evolutionary History of Mycobacterium tuberculosis LAM9 Sublineage in the Americas as Two Distinct Subpopulations Revealed by Bayesian Analyses

Tuberculosis (TB) remains broadly present in the Americas despite intense global efforts for its control and elimination. Starting from a large dataset comprising spoligotyping (n = 21183 isolates) and 12-loci MIRU-VNTRs data (n = 4022 isolates) from a total of 31 countries of the Americas (data extracted from the SITVIT2 database), this study aimed to get an overview of lineages circulating in the Americas. A total of 17119 (80.8%) strains belonged to the Euro-American lineage 4, among which the most predominant genotypic family belonged to the Latin American and Mediterranean (LAM) lineage (n = 6386, 30.1% of strains). By combining classical phylogenetic analyses and Bayesian approaches, this study revealed for the first time a clear genetic structuration of LAM9 sublineage into two subpopulations named LAM9C1 and LAM9C2, with distinct genetic characteristics. LAM9C1 was predominant in Chile, Colombia and USA, while LAM9C2 was predominant in Brazil, Dominican Republic, Guadeloupe and French Guiana. Globally, LAM9C2 was characterized by higher allelic richness as compared to LAM9C1 isolates. Moreover, LAM9C2 sublineage appeared to expand close to twenty times more than LAM9C1 and showed older traces of expansion. Interestingly, a significant proportion of LAM9C2 isolates presented typical signature of ancestral LAM-RD^Rio MIRU-VNTR type (224226153321). Further studies based on Whole Genome Sequencing of LAM strains will provide the needed resolution to decipher the biogeographical structure and evolutionary history of this successful family.

Clustured regularly interspersed short palindromic repeats (CRISPR) genetic diversity studies as a mean to reconstruct the evolution of the Mycobacterium tuberculosis complex

The natural history of tuberculosis may be tackled by various means, among which the record of molecular scars that have been registered by the Mycobacterium tuberculosis complex (MTBC) genomes transmitted from patient to patient for tens of thousands years and possibly more. Recently discovered polymorphic loci, the CRISPR sequences, are indirect witnesses of the historical phage-bacteria struggle, and may be related to the time when the ancestor of today's tubercle bacilli were environmental bacteria, i.e. before becoming intracellular parasites. In this article, we present what are CRISPRs and try to summarize almost 20 years of research results obtained using the genetic diversity of the CRISPR loci in MTBC as a perspective for studying new models. We show that the study of the diversity of CRISPR sequences, thanks to «spoligotyping», has played a great role in our global understanding of the population structure of MTBC.

A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile

Tuberculosis (TB) remains a significant public health problem worldwide, but the ecology of the prevalent mycobacterial strains, and their transmission, can vary depending on country and region. Chile is a country with low incidence of TB, that has a geographically isolated location in relation to the rest of South American countries due to the Andes Mountains, but recent migration from neighboring countries has changed this situation. We aimed to assess the genotypic diversity of Mycobacterium tuberculosis complex (MTBC) strains in Santiago, Chile, and compare with reports from other Latin-American countries. We analyzed MTBC isolates from pulmonary tuberculosis cases collected between years 2008 and 2013 in Central Santiago, using two genotyping methods: spoligotyping and 12-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTRs). Data obtained were analyzed and compared to the SITVIT2 database. Mean age of the patients was 47.5 years and 61% were male; 11.6% were migrants. Of 103 strains (1 isolate/patient) included, there were 56 distinct spoligotype patterns. Of these, 16 strains (15.5%) corresponded to orphan strains in the SITVIT2 database, not previously reported. Latin American and Mediterranean (LAM) (34%) and T (33%) lineages were the most prevalent strains, followed by Haarlem lineage (16.5%). Beijing family was scarcely represented with only two cases (1.9%), one of them isolated from a Peruvian migrant. The most frequent clustered spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%), and SIT37/T3 (6.8%). We conclude that LAM and T genotypes are the most prevalent genotypes of MTBC in Santiago, Chile, and together correspond to almost two thirds of analyzed strains, which is similar to strain distribution reported from other countries of Latin America. Nevertheless, the high proportion of SIT37/T3, which was rarely found in other Latin American countries, may underline a specific history or demographics of Chile related to probable human migrations and evolutions.

Unraveling Mycobacterium tuberculosis genomic diversity and evolution in Lisbon, Portugal, a highly drug resistant setting

Background

Multidrug- (MDR) and extensively drug resistant (XDR) tuberculosis (TB) presents a challenge to disease control and elimination goals. In Lisbon, Portugal, specific and successful XDR-TB strains have been found in circulation for almost two decades.

Results

In the present study we have genotyped and sequenced the genomes of 56 Mycobacterium tuberculosis isolates recovered mostly from Lisbon. The genotyping data revealed three major clusters associated with MDR-TB, two of which are associated with XDR-TB. Whilst the genomic data contributed to elucidate the phylogenetic positioning of circulating MDR-TB strains, showing a high predominance of a single SNP cluster group 5. Furthermore, a genome-wide phylogeny analysis from these strains, together with 19 publicly available genomes of Mycobacterium tuberculosis clinical isolates, revealed two major clades responsible for M/XDR-TB in the region: Lisboa3 and Q1 (LAM).

The data presented by this study yielded insights on microevolution and identification of novel compensatory mutations associated with rifampicin resistance in rpoB and rpoC. The screening for other structural variations revealed putative clade-defining variants. One deletion in PPE41, found among Lisboa3 isolates, is proposed to contribute to immune evasion and as a selective advantage. Insertion sequence (IS) mapping has also demonstrated the role of IS6110 as a major driver in mycobacterial evolution by affecting gene integrity and regulation.

Conclusions

Globally, this study contributes with novel genome-wide phylogenetic data and has led to the identification of new genomic variants that support the notion of a growing genomic diversity facing both setting and host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-991) contains supplementary material, which is available to authorized users.

Strain classification of Mycobacterium tuberculosis isolates in Brazil based on genotypes obtained by spoligotyping, mycobacterial interspersed repetitive unit typing and the presence of large sequence and single nucleotide polymorphism

Rio de Janeiro is endemic for tuberculosis (TB) and presents the second largest prevalence of the disease in Brazil. Here, we present the bacterial population structure of 218 isolates of Mycobacterium tuberculosis, derived from 186 patients that were diagnosed between January 2008 and December 2009. Genotypes were generated by means of spoligotyping, 24 MIRU-VNTR typing and presence of fbpC¹⁰³, RD^Rio and RD174. The results confirmed earlier data that predominant genotypes in Rio de Janeiro are those of the Euro American Lineages (99%). However, we observed differences between the classification by spoligotyping when comparing to that of 24 MIRU-VNTR typing, being respectively 43.6% vs. 62.4% of LAM, 34.9% vs. 9.6% of T and 18.3% vs. 21.5% of Haarlem. Among isolates classified as LAM by MIRU typing, 28.0% did not present the characteristic spoligotype profile with absence of spacers 21 to 24 and 32 to 36 and we designated these conveniently as “LAM-like”, 79.3% of these presenting the LAM-specific SNP fbpC¹⁰³. The frequency of RD^Rio and RD174 in the LAM strains, as defined both by spoligotyping and 24 MIRU-VNTR loci, were respectively 11% and 15.4%, demonstrating that RD174 is not always a marker for LAM/RD^Rio strains. We conclude that, although spoligotyping alone is a tool for classification of strains of the Euro-American lineage, when combined with MIRU-VNTRs, SNPs and RD typing, it leads to a much better understanding of the bacterial population structure and phylogenetic relationships among strains of M. tuberculosis in regions with high incidence of TB.

Insights into the population structure of Mycobacterium tuberculosis using spoligotyping and RDRio in a southeastern Brazilian prison unit

Tuberculosis (TB) is still a serious public health problem, continuing to be an important threat for confined populations. We used spoligotyping to estimate the genotypic clades of Mycobacterium tuberculosis isolates from inmates in two blocks in a southeastern Brazilian prison unit, with TB incidence rate of 8185/100.000. The Latin American Mediterranean (LAM) clade is well represented in the country, and the LAM specific molecular markers, RD(Rio) large sequence polymorphism and the SNP on the Rv3062 [ligB(1212)], were used to characterize spoligotype signatures from prison isolates. Typing of RD(Rio) and ligB increase LAM clade from 66.7% (n=72/108) to 69.4% (n=75). The LAM2 SIT17 (n=23) and SIT179 (n=12) signatures comprised one third of all isolates, followed by Haarlem (11.5%, n=12), T (8.7%, n=9) and X (5.7%, n=6) clades. Strains with unknown signatures represented 5.5% (n=6), and four (3.7%) did not match any lineage. We observed RD(Rio) among 64 (59.2%) isolates, and 54 (50%) were of the LAM clade. In particular, the LAM2/RD(Rio) sub-lineage was significantly associated with clustering (p=0.02) and its frequency was higher (32%) when compared to that of the previous general TB cases in RJ (4.29%). Overall cluster frequency defined by spoligotyping/IS6110-RFLP was 62%. The two evolutionary markers helped to evaluate some LAM signature misconceptions and demonstrate that LAM2/RD(Rio) was found with high frequency, hitherto being unnoticed. All these data, allied to high clustering, imply that public health measures to minimize the escalation of TB in prison is essential, and both spoligotyping as well as RD(Rio) would be useful tools to monitor the effects of the measures with respect to M. tuberculosis lineage variation.

Mycobacterium tuberculosis Latin American-Mediterranean family and its sublineages in the light of robust evolutionary markers

Mycobacterium tuberculosis has a clonal population structure, and the Latin American-Mediterranean (LAM) family is one of the largest and most widespread within this species, showing evidence for remarkable pathobiology and a confusing phylogeny. Here, we applied robust phylogenetic markers to study the evolution of the LAM family and its major sublineages circulating in Russia and neighboring countries. A total of 250 M. tuberculosis isolates were confirmed to belong to the LAM family based on the analysis of the LAM-specific single-nucleotide polymorphisms (SNPs) in the Rv3062 and Rv0129c genes. At this stage, the family status was rectified for 121 isolates misleadingly assigned by CRISPR spoligotyping to non-LAM families (T1- or T5-RUS1). Consequently, the reestimated LAM prevalence rate increased 2-fold in Russia and Kazakhstan and 4-fold in Belarus. The majority (91.8 to 98.7%) of the LAM isolates from all three countries belonged to the LAM-RUS sublineage. In contrast, the Ibero-American LAM RD-Rio sublineage was identified in only 7 Russian isolates. Taken together, our findings and further analyses suggest a monophyletic origin of LAM-RUS: at a historically distant time, in Russia, in a small founding bacterial/human population. Its dissemination pattern and high prevalence rate in Northern Eurasia may indicate a long-term coexistence of the LAM-RUS sublineage and local human populations hypothetically leading to coadaptation and reduced pathogenicity of the relatively more ancient clones, such as spoligotype international type 254 (SIT254), compared to the more recent SIT252 and SIT266 clones. In contrast, rare LAM RD-Rio isolates were likely brought to Russia through occasional human contact. The spread of RD-Rio strains is not as global as commonly claimed and is determined largely by human migration flows (rather than by pathobiological properties of these strains). Consequently, a host population factor appears to play a major role in shaping the in situ dissemination pattern of the imported strains in an autochthonous population.

Single nucleotide polymorphism (SNP) analysis used for the phylogeny of the Mycobacterium tuberculosis complex based on a pyrosequencing assay

Background

Different polymorphisms have been described as markers to classify the lineages of the Mycobacterium tuberculosis complex. The analysis of nine single nucleotide polymorphisms (SNPs) was used to describe seven SNPs cluster groups (SCGs). We attempted to classify those strains that could not been categorized into lineages by the genotyping methods used in the routine testing.

Results

The M. tuberculosis complex isolates collected in 2010 in our region were analysed. A new method based on multiplex-PCRs and pyrosequencing to analyse these SNPs was designed. For the pyrosequencing assay nine SNPs that defined the seven SCGs were selected from the literature: 1977, 74092, 105139, 232574, 311613, 913274, 2460626, 3352929 and gyrA95. In addition, SNPs in katG⁴⁶³, mgtC¹⁸², Ag85C¹⁰³ and RD^Rio deletion were detected.

Conclusions

This work has permitted to achieve a better classification of Aragonian strains into SCGs and in some cases, to assign strains to its certain lineage. Besides, the description of a new pattern shared by two isolates “SCG-6c” reinforces the interest of SNPs to follow the evolution of M. tuberculosis complex.

Mycobacterium tuberculosis belonging to family LAM and sublineage RD(Rio): common strains in Southern Brazil for over 10 years

A sublineage of Mycobacterium tuberculosis called RD(Rio) was described in 2007. Although only recently described, this strain may have been present previously in the population, and its identification in clinical isolates will elucidate bacterial transmission dynamics and host-pathogen interactions. This study evaluated the clonal diversity of the RD(Rio) sublineage in clinical isolates from Rio Grande-RS obtained between 1998 and 2001. Among the 45 samples analyzed by the MIRU-VNTR method, there were six clusters with two samples each and 33 orphan strains with unique pattern. The strains were distributed across several different lineages including LAM (34.04%), × (14.89%), Haarlem (12.77%), UgandaI (10.64%), S (4.26%), NEW-1 (2.13%) and Cameroon (2.13%); 14.89% of the strains matched to multiple lineages. RD(Rio) strains were present in 28.9% of the samples and 81.25% of the identified strains belonged to the LAM family. The high clonal diversity observed in this study is a constant feature in this region. The RD(Rio) sublineage has been in Rio Grande-RS since 1998. The continued monitoring of RD(Rio) in clinical isolates will enhance the understanding of its epidemiological significance.