Mycobacterium Tuberculosis Lineages Research Articles

Whole-genome sequencing-based epidemiological analysis of anti-tuberculosis drug resistance genes in Japan in 2007: Application of the Genome Research for Asian Tuberculosis (GReAT) database

We investigated the lineages of Mycobacterium tuberculosis (Mtb) isolates from the RYOKEN study in Japan in 2007 and the usefulness of genotypic drug susceptibility testing (DST) using the Genome Research for Asian Tuberculosis (GReAT) database. In total, 667 isolates were classified into lineage 1 (4.6%), lineage 2 (0.8%), lineage 2/Beijing (72.1%), lineage 3 (0.5%), and lineage 4 (22.0%). The nationality, gender, and age groups associated with the isolates assigned to lineage 1 were significantly different from those associated with other lineages. In particular, isolates of lineage 1.2.1 (EAI2) formed sub-clusters and included a 2,316-bp deletion in the genome. The proportion of the isolates resistant to at least one anti-tuberculosis (TB) drug was 10.8%, as determined by either the genotypic or phenotypic method of DST. However, the sensitivities to isoniazid, streptomycin, and ethambutol determined by the genotypic method were low. Thus, unidentified mutations in the genome responsible for drug resistance were explored, revealing previously unreported mutations in the katG, gid, and embB genes. This is the first nationwide report of whole-genome analysis of TB in Japan.

Indo-Oceanic Mycobacterium tuberculosis strains from Thailand associated with higher mortality.

SETTING: This study was conducted among tuberculosis (TB) patients in a highly endemic Thai province.OBJECTIVE: To evaluate the association between different Mycobacterium tuberculosis lineages and clinical characteristics, especially mortality.DESIGN: We enrolled 1,304 TB patients registered from 2002-2011 with culture isolates whose lineages were identified by specific regions of deletion. Data on mortality within 1 year of follow-up were extracted from the registration system and hospital records. Mortality-associated risk factors, including bacterial lineages, as independent variables were analysed using Cox regression models.RESULTS: Of 1,304 isolates, 521 (40.0%) and 582 (44.6%) belonged to Indo-Oceanic and East-Asian lineages, respectively. Indo-Oceanic strains significantly increased the mortality risk compared with East-Asian strains (adjusted hazard ratio [aHR] 1.42, 95%CI 1.02-1.99) or modern lineages (aHR 1.49, 95%CI 1.08-2.06) in the 172 patients who died within 1 year after TB diagnosis. The former also caused significantly higher mortality than modern lineages among patients who died within 6 months after TB diagnosis (aHR 1.62, 95%CI 1.12-2.35). No significant association was found between drug resistance and death.CONCLUSION: In Thailand, the Indo-Oceanic lineage of M. tuberculosis increased mortality risk compared with modern lineages or the East-Asian lineage, the latter being considered highly virulent in previous studies.

The role of IS6110 in micro- and macroevolution of Mycobacterium tuberculosis lineage 2

The insertion sequence 6110 (IS6110) is the most studied transposable element in the Mycobacterium tuberculosis complex species. The element plays a significant role in genome plasticity of this important human pathogen, but still many causes and consequences of its transposition have not been fully studied. Here, we analyzed insertion sites for 902 Mycobacterium tuberculosis lineage 2 strains using whole-genome sequencing data. In total, 17,972 insertions were found, corresponding to 827 independent positions in the genome of the reference strain H37Rv. To trace the history of IS6110 expansion since proto-Beijing strains up to modern sublineages, we looked at the distribution of IS6110 across the genome-wide SNP-based phylogenetic tree. This analysis demonstrated a stepwise transposition of IS6110 that occurs by «copy-and-paste» mechanism. Additionally, we detected evolutionary-scale and sublineage-specific integration sites, which can be used for typing and for understanding the reasons for the success of the lineage. A significant part of such insertions affected the genes that are essential for the pathogen. Finally, we identified and confirmed deletions that occurred between differently oriented elements, which is uncommon for this family of insertion elements and appears to be another mechanism of genome variability.

Severe Clinical Outcomes of Tuberculosis in Kharkiv Region, Ukraine, Are Associated with Beijing Strains of Mycobacterium tuberculosis.

Genotypic variation in Beijing lineages of Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis (TB), has been associated with hyper virulence and the spread of extensively and multiple drug (X/MDR) resistant MTB strains in Eastern Europe, Central Asia, and East Asia. The clinical outcomes of 215 new cases of TB among the population of the Kharkiv region of Eastern Ukraine were analyzed to uncover factors associated with severe infection. Infecting MTB strains were profiled by 5 locus exact tandem repeats (ETRs) and 15 locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) genotyping. Among diverse MTB genotypes discovered in Ukraine, the Beijing genotype (MIRU-VNTR 42425) was significantly associated with risk factors for severe outcomes of disease in the study population, including TB/HIV co-infection and treatment failure. Strain replacement (superinfection) was observed in 10 patients, suggesting repeated exposure to novel MTB strains in hospital or community settings. Inclusion of MTB genotyping data may identify at-risk patients and improve treatment adherence to prevent X/MDR development for effective public health response against tuberculosis in Ukraine.

Mycobacterium tuberculosis lineage 1 genetic diversity in Pará, Brazil, suggests common ancestry with east-African isolates potentially linked to historical slave trade

Lineage 1 (L1) is one of seven Mycobacterium tuberculosis complex (MTBC) lineages. The objective of this study was to improve the complex taxonomy of L1 using phylogenetic SNPs, and to look for the origin of the main L1 sublineage prevalent in Para, Brazil. We developed a high-throughput SNPs-typing assay based on 12-L1-specific SNPs. This assay allowed us to experimentally retrieve SNP patterns on nine of these twelve SNPs in 277 isolates previously tentatively assigned to L1 spoligotyping-based sublineages. Three collections were used: Pará-Brazil (71); RIVM, the Netherlands (102), Madagascar (104). One-hundred more results were generated in Silico using the PolyTB database. Based on the final SNPs combination, the samples were classified into 11 clusters (C1-C11). Most isolates within a SNP-based cluster shared a mutual spoligotyping-defined lineage. However, L1/EAI1-SOM (SIT48) and L1/EAI6-BGD1 (SIT591) showed a poor correlation with SNP data and are not monophyletic. L1/EAI8-MDG and L1/EAI3-IND belonged to C5; this result suggests that they share a common ancestor. L1.1.3/SIT129, a spoligotype pattern found in SNPs-cluster C6, was found to be shared between Pará/Brazil and Malawi. SIT129 was independently found to be highly prevalent in Mozambique, which suggests a migration history from East-Africa to Brazil during the 16th–18th slave trade period to Northern Brazil.

Spoligotyping analysis of Mycobacterium tuberculosis in Khyber Pakhtunkhwa area, Pakistan.

Background: Spoligotyping is a reproducible, reverse hybridization approach for genotyping of Mycobacterium tuberculosis complex (MTBC). Molecular typing of MTBC is helpful for understanding and controlling tuberculosis epidemics.Methods: Spoligotyping was performed on 166 clinical isolates of Mycobacterium tuberculosis (MTB) collected from 25 districts of Khyber Pakhtunkhwa, Pakistan. Results were compared to SITVIT2, an online database developed by the Institut Pasteur de la Guadeloupe, France.Results: Spoligotyping results showed that 145 strains (88%) displayed known patterns while 21 (12%) were new. Lineage 3/Central Asian strain (L3/CAS) was the predominant family (73%, χ2=19.9, P=0.001), followed by L2/Beijing (5.4%) and L4 (4.2%). L3/CAS1-Delhi was the major sublineage (82%) among the L3/CAS family (χ2=664, P=0.0001). Analysis showed that the majority of the clinical isolates with an unknown pattern had an evolutionary link with the L3/CAS strain, and nine (5.4%) of the unknown strains were epidemiologically linked and were tentatively named L3/CAS-KP (Khyber Pakhtunkhwa).Conclusion: The present study demonstrated that L3/CAS is the predominant lineage of MTB, widely distributed in different areas of the Khyber Pakhtunkhwa province of Pakistan. Spoligotyping patterns of some clinical isolates could not be matched to other reported patterns in an international database. Other tools, such as mycobacterial interspersed repetitive unit–variable number tandem repeat (MIRU-VNTR), will be helpful in future investigations into the epidemiological characteristics of clinical isolates in the Khyber Pakhtunkhwa province.

An eight-year epidemiologic study of head and neck tuberculosis in Texas, USA

BackgroundHead and neck tuberculosis (HNTB) especially cervical lymphadenopathy are the most common extrapulmonary indications of TB, but remain a diagnostic challenge. In this study, we describe and analyze the epidemiologic characteristics of HNTB on a population-level. Materials and methodsWe retrospectively assessed 547 HNTB cases reported to the Centers for Disease Control and Prevention's TB Genotyping Information Management System in Texas from 2009 to 2016 and compare and contrast differences between diagnosed exclusively HNTB and HNTB with concurrent pulmonary tuberculosis (PTB). ResultsThe majority of patients with HNTB were diagnosed with cervical lymphadenopathy (96.9%), age 25–44 (47.3%) and female (52.7%). Co-infection with human immunodeficiency virus, being homeless, excessive alcohol use within the past 12 months and drug use were more frequently seen for HNTB with concurrent pulmonary involvement compared to reported patients with exclusively HNTB. The highest prevalence of Mycobacterium tuberculosis lineage in exclusively HNTB was Euro-American L4 (52.3%), followed by Indo-Oceanic L1 (21.5%) and East-Asian L2 (13.1%). One multidrug resistant TB case was identified. Seven deaths were reported during treatment. ConclusionOur findings provide a better understanding of the epidemiology of HNTB and characteristics associated with the disease at the population-level, which is important in managing HNTB patients.

Multiple Introductions of Mycobacterium tuberculosis Lineage 2–Beijing Into Africa Over Centuries

The Lineage 2–Beijing (L2–Beijing) sub-lineage of Mycobacterium tuberculosis has received much attention due to its high virulence, fast disease progression, and association with antibiotic resistance. Despite several reports of the recent emergence of L2–Beijing in Africa, no study has investigated the evolutionary history of this sub-lineage on the continent. In this study, we used whole genome sequences of 781 L2 clinical strains from 14 geographical regions globally distributed to investigate the origins and onward spread of this lineage in Africa. Our results reveal multiple introductions of L2–Beijing into Africa linked to independent bacterial populations from East- and Southeast Asia. Bayesian analyses further indicate that these introductions occurred during the past 300 years, with most of these events pre-dating the antibiotic era. Hence, the success of L2–Beijing in Africa is most likely due to its hypervirulence and high transmissibility rather than drug resistance.

Prevalence, Drug Resistance, and Genotypic Diversity of the Mycobacterium tuberculosis Beijing Family in Ecuador.

The Beijing family, the most successful Mycobacterium tuberculosis lineage, is considered hypervirulent, associated with clustering and has a strong association with multidrug-resistant tuberculosis. The Beijing strains have spread worldwide and also to Latin America. Genotyping of a countrywide collection of 380 M. tuberculosis strains from Ecuador, with 24-loci mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR), revealed only six Beijing strains, but four of these were MDR-TB. There was no clustering as all six strains had very distinct MIRU-VNTR profiles that have not been reported in the rest of Latin America. Although active transmission for Beijing has been described for the neighboring countries Peru and Colombia, there is no evidence that Beijing strains in Ecuador are more frequently transmitted than other strains. Moreover, the low prevalence (1.6%) of the Beijing sublineage in Ecuador challenges the concept of hyperadaptability and transmissibility of the Beijing strains in our country.

Macro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database

In order to provide a global overview of genotypic, epidemiologic, demographic, phylogeographical, and drug resistance characteristics related to the prevailing tuberculosis (TB) epidemic, we hereby report an update of the 6th version of the international genotyping database SITVIT2. We also make all the available information accessible through a dedicated website (available at http://www.pasteur-guadeloupe.fr:8081/SITVIT2). Thanks to the public release of SITVIT2 which is currently the largest international multimarker genotyping database with a compilation of 111,635 clinical isolates from 169 countries of patient origin (131 countries of isolation, representing 1032 cities), our major aim is to highlight macro- and micro-geographical cleavages and phylogeographical specificities of circulating Mycobacterium tuberculosis complex (MTBC) clones worldwide. For this purpose, we retained strains typed by the most commonly used PCR-based methodology for TB genotyping, i.e., spoligotyping based on the polymorphism of the direct repeat (DR) locus, 5-loci Exact Tandem Repeats (ETRs), and MIRU-VNTR minisatellites used in 12–, 15–, or 24–loci formats. We describe the SITVIT2 database and integrated online applications that permit to interrogate the database using easy drop-down menus to draw maps, graphics and tables versus a long list of parameters and variables available for individual clinical isolates (year and place of isolation, origin, sex, and age of patient, drug-resistance, etc.). Available tools further allow to generate phylogenetical snapshot of circulating strains as Lineage-specific WebLogos, as well as minimum spanning trees of their genotypes in conjunction with their geographical distribution, drug-resistance, demographic, and epidemiologic characteristics instantaneously; whereas online statistical analyses let a user to pinpoint phylogeographical specificities of circulating MTBC lineages and conclude on actual demographic trends. Available associated information on gender (n = 18,944), age (n = 16,968), drug resistance (n = 19,606), and HIV serology (n = 2673), allowed to draw some important conclusions on TB geo-epidemiology; e.g. a positive correlation exists between certain Mycobacterium tuberculosis lineages (such as CAS and Beijing) and drug resistance (p-value<.001), while other lineages (such as LAM, X, and BOV) are more frequently associated with HIV-positive serology (p-value<.001). Besides, availability of information on the year of isolation of strains (range 1759–2012), also allowed to make tentative correlations between drug resistance information and lineages – portraying probable evolution trends over time and space. To conclude, the present approach of geographical mapping of predominant clinical isolates of tubercle bacilli causing the bulk of the disease both at country and regional level in conjunction with epidemiologic and demographic characteristics allows to shed new light on TB geo-epidemiology in relation with the continued waves of peopling and human migration.

Mycobacterium africanum (Lineage 6) shows slower sputum smear conversion on tuberculosis treatment than Mycobacterium tuberculosis (Lineage 4) in Bamako, Mali.

ObjectiveAncestral M. tuberculosis complex lineages such as M. africanum are underrepresented among retreatment patients and those with drug resistance. To test the hypothesis that they respond faster to TB treatment, we determined the rate of smear conversion of new pulmonary tuberculosis patients in Bamako, Mali by the main MTBc lineages.MethodsBetween 2015 and 2017, we conducted a prospective cohort study of new smear positive pulmonary tuberculosis patients in Bamako. Confirmed MTBc isolates underwent genotyping by spoligotyping for lineage classification. Patients were followed at 1 month (M), 2M and 5M to measure smear conversion in auramine (AR) and Fluorescein DiAcetate (FDA) vital stain microscopy.ResultAll the first six human MTBc lineages were represented in the population, plus M. bovis in 0.8% of the patients. The most widely represented lineage was the modern Euro-American lineage (L) 4, 57%, predominantly the T family, followed by L6 (M. africanum type 2) in 22.9%. Ancestral lineages 1, 5, 6 and M. bovis combined amounted to 28.8%. Excluding 25 patients with rifampicin resistance, smear conversion, both by AR and FDA, occurred later in L6 compared to L4 (HR 0.80 (95% CI 0.66–0.97) for AR, and HR 0.81 (95%CI 0.68–0.97) for FDA). In addition we found that HIV negative status, higher BMI at day 0, and patients with smear grade at baseline ≤ 1+ were associated with earlier smear conversion.ConclusionThe six major human lineages of the MTBc all circulate in Bamako. Counter to our hypothesis, we found that patients diseased with modern M. tuberculosis complex L4 respond faster to TB treatment than those with M. africanum L6.

A Closer Look at the Genomic Variation of Geographically Diverse Mycobacterium abscessus Clones That Cause Human Infection and Disease.

Mycobacterium abscessus is a multidrug resistant bacterium that causes pulmonary and extrapulmonary disease. The reported prevalence of pulmonary M. abscessus infections appears to be increasing in the United States (US) and around the world. In the last five years, multiple studies have utilized whole genome sequencing to investigate the genetic epidemiology of two clinically relevant subspecies, M. abscessus subsp. abscessus (MAB) and M. abscessus subsp. massiliense (MMAS). Phylogenomic comparisons of clinical isolates revealed that substantial proportions of patients have MAB and MMAS isolates that belong to genetically similar clusters also known as ‘dominant clones’. Unlike the genetic lineages of Mycobacterium tuberculosis that tend to be geographically clustered, the MAB and MMAS clones have been found in clinical populations from the US, Europe, Australia and South America. Moreover, the clones have been associated with worse clinical outcomes and show increased pathogenicity in macrophage and mouse models. While some have suggested that they may have spread locally and then globally through ‘indirect transmission’ within cystic fibrosis (CF) clinics, isolates of these clones have also been associated with sporadic pulmonary infections in non-CF patients and unrelated hospital-acquired soft tissue infections. M. abscessus has long been thought to be acquired from the environment, but the prevalence, exposure risk and environmental reservoirs of the dominant clones are currently not known. This review summarizes the genomic studies of M. abscessus and synthesizes the current knowledge surrounding the geographically diverse dominant clones identified from patient samples. Furthermore, it discusses the limitations of core genome comparisons for studying these genetically similar isolates and explores the breadth of accessory genome variation that has been observed to date. The combination of both core and accessory genome variation among these isolates may be the key to elucidating the origin, spread and evolution of these frequent genotypes.

Toxin-antitoxin systems shows variability among Mycobacterium tuberculosis lineages.

The toxin-antitoxin (TA) systems are operons involved in the formation of persistent cells and in stress situations in microorganism. TA systems are widely distributed in Mycobacterium tuberculosis (MTB). The objective of this study was to determine the distribution and variability of protein sequences of TA systems in seven MTB lineages. Protein prediction on 73 genomes of different lineage was made using Prodigal, and profile hidden Markov models (PHMMs) of 225 reference TA proteins were constructed with HMMER. An homology search was made comparing the predicted proteins to PHMMs. Consecutive proteins that matched PHMMs (forming an operon) were selected. Thereafter, clustering of orthologous genes was made for further mutation scanning through multiple alignments. All proteins found belong to TA types II and IV, and 45 proteins were found completely conserved. Six uncharacterized antitoxins homologous to VapB11, VapB5, VapB45, VapB13, ParD1 and RelB were found. Multiple alignments revealed differences among lineages with specific mutations suitable for genotyping. Significant changes in amino acid sequences caused by frameshift mutations were found in 46 TA proteins.

760. The Diversity of Mycobacterium tuberculosis in India Is Underestimated and Underreported

BackgroundThe seven lineages of Mycobacterium tuberculosis differ in virulence, transmissibility, drug resistance, and immune responses. We reviewed spoligotype data published from India to determine the distribution and diversity of TB lineages in India.MethodsSpoligoytpe patterns were extracted from studies published from India and compared against the SpolSimilarity database to determine Spoligotype International Type (SIT) number, sub-lineage and lineage. Minimum spanning trees (MSTs) were created with Phyloviz.ResultsThirty-four studies list 1,473 spoligotypes containing 8,242 isolates (Table 1), with 6,733 isolates belonging to 472 SITs. Of the remaining 1,001 patterns (1,509 isolates), 351 patterns (491 isolates) match orphans, while 650 patterns (1,018 isolates) are unreported; these 1,509 isolates could potentially form 532 new SITs and 469 orphans. The most common SITs are SIT11 (Lineage 1; 1,218 isolates, 14.78%), SIT26 (Lineage 3; 982,11.91%) and SIT1 (Lineage 2; 618, 7.5%). Figure 1 demonstrates regional differences in lineage distribution. Figure 2 is an MST of all isolates by lineage (“not assigned” are known SITs with no lineage and “unknowns” are previously unreported patterns). Figure 3 is an MST by superfamily, “others” referring to patterns with no known superfamily.ConclusionLineages 1 and 3 dominate in India, in contrast to lineages 2 and 4 prevalent worldwide. The 1,473 SITs reflect the significant diversity of M. tuberculosis in India; the 1,001 orphan and previously unreported patterns suggest this diversity is underestimated. Understanding M. tuberculosis diversity is necessary for clinical and public health interventions to control India’s TB epidemic.Table 1.LineageSuperfamilyNumber of SITsNumber of IsolatesPercentage (%)Lineage 1EAI2343,21439.00Manu704815.84Total3043,69544.83Lineage 2Beijing116557.95Lineage 3CAS1821,86522.63Lineage 4H26891.08LAM20550.67S5130.16T1114765.78Ural11370.45X17560.68Total1907268.81UnknownUnknown7861,30115.79Total1,4738,242100.00Figure 1.Figure 2.Figure 3.Disclosures All authors: No reported disclosures.

The role of Mycobacterium tuberculosis lineages on lung tissue damage and TNF-α level among tuberculosis patients, Indonesia

BackgroundMycobacterium tuberculosis (MTB) is divided into ancient and modern lineage. Both lineages have different virulence property and therefore associated with different immune responses and disease severity in tuberculosis (TB) patients. The aim of this study was to determine the association between MTB lineages and the level of tumor necrosis factor alpha (TNF-α), and lung tissue damage among Indonesian TB patients. MethodsThirty new and active TB patients were recruited randomly from Dr. Soewandhie Hospital, Surabaya, Indonesia. MTB isolates were isolated from bronchoalveolar lavage fluid (BALF) and their lineages were determined by primer-specific PCR targeting TbD1 and RD9 region gene. The degree of lung tissue damage was assessed and classified using NICE Scoring System and TNF-α level of BALF was measured by ELISA. ResultsMTB were detected in all patients in which 17 and 13 isolates were classified as modern and ancient lineage, respectively. Modern lineage was associated with the degree of lung damage (odds ratio [OR]: 3.97; 95% CI: 1.57–6.37), p < 0.001). This lineage was also associated with high TNF-α level (OR: 2.06; 95% CI: 1.31–47.00, p=0.026). In addition, the level of TNF-α was higher in severe compared to mild lung damage cases (p = 0.001). ConclusionModern lineage is likely associated with high TNF-α level and increases the risk of having severe lung damage three times among TB patients from Indonesia. Therefore, laboratory test to determine MTB lineages might crucial to be integrated in National Tuberculosis Control Program in the country.

Identifying Mixed Mycobacterium tuberculosis Infection and Laboratory Cross-Contamination during Mycobacterial Sequencing Programs.

The detection of laboratory cross-contamination and mixed tuberculosis infections is an important goal of clinical mycobacteriology laboratories. The objective of this study was to develop a method to detect mixtures of different Mycobacterium tuberculosis lineages in laboratories performing mycobacterial next-generation sequencing (NGS). The setting was the Public Health England National Mycobacteriology Laboratory Birmingham, which performs Illumina sequencing on DNA extracted from positive mycobacterial growth indicator tubes. We analyzed 4,156 samples yielding M. tuberculosis from 663 MiSeq runs, which were obtained during development and production use of a diagnostic process using NGS. The counts of the most common (major) variant and all other variants (nonmajor variants) were determined from reads mapping to positions defining M. tuberculosis lineages. Expected variation was estimated during process development. For each sample, we determined the nonmajor variant proportions at 55 sets of lineage-defining positions. The nonmajor variant proportion in the two most mixed lineage-defining sets (F2 metric) was compared with that of the 47 least-mixed lineage-defining sets (F47 metric). The following three patterns were observed: (i) not mixed by either metric; (ii) high F47 metric, suggesting mixtures of multiple lineages; and (iii) samples compatible with mixtures of two lineages, detected by differential F2 metric elevations relative to F47. Pattern ii was observed in batches, with similar patterns in the M. tuberculosis H37Rv control present in each run, and is likely to reflect cross-contamination. During production, the proportions of samples in the patterns were 97%, 2.8%, and 0.001%, respectively. The F2 and F47 metrics described could be used for laboratory process control in laboratories sequencing M. tuberculosis genomes.

New epidemic cluster of pre-extensively drug resistant isolates of Mycobacterium tuberculosis Ural family emerging in Eastern Europe

BackgroundUral genetic family is a part of the Euro-American lineage of Mycobacterium tuberculosis and is endemic in Northern Eurasia (former Soviet Union [FSU]). These strains were long described as drug susceptible and of low virulence, but recent studies reported an increasing circulation of the multidrug-resistant (MDR) and extensively drug-resistant (XDR) Ural strains. Here, we analyzed all publicly available whole genome sequence data of Ural genotype isolates, in order to elucidate their phylogenomic diversity with a special focus on MDR and potentially epidemic clones.ResultsA total of 149 M. tuberculosis genomes of Ural isolates from FSU countries were mined from the GMTV database and TB-ARC project. We identified 6002 variable amino acid positions that were assessed for functional significance and used to build ML, NJ trees and for Bayesian TMRCA estimation. Three robust monophyletic clades were identified: Clade A (31 isolates from Russia, Belarus, Moldova), Clade B (52 isolates from Russia), and Clade C (37 isolates from Moldova, 2 from Belarus). Clade C was significantly associated with XDR or pre-XDR status compared to the pooled Clades A and B (33/39 versus 5/83, P < 0.0001). Time of origin was estimated for Clade A at 77.7–137 years ago and for Clade B at 56.3–99.2 years ago compared to the significantly more recent origin for Clade C. in silico spoligotyping identified signatures specific of the Clade A (spoligotype SIT35), and Clades B and C (both SIT262).ConclusionsA genetically compact and evolutionarily young Ural Clade C, likely originated after collapse of the Soviet Union, and reached epidemic proportions in Moldova in the last 20 years. This epidemic pre-XDR clone (mostly rifampin, isoniazid and kanamycin resistant) is characterized by a specific combination of mutations: KatG Ser315Thr, fabG1 -15C > T, RpoB Ser450Leu, RpsL Lys88Arg, eis -12G > A and EmbB Ser297Ala/T > G. Its further dissemination may occur towards both Russia and European Union and should be taken into consideration by health authorities. The identified spoligotyping signatures can serve for rapid preliminary detection and surveillance of the more hazardous pre-XDR associated strains of the Ural family, both in populations from countries of their endemic circulation and migrant communities.

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.

Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates

Tuberculosis presents a global health challenge. Mycobacterium tuberculosis is divided into several lineages, each with a different geographical distribution. M. tuberculosis lineage 1 (L1) is common in the high-burden areas in East Africa and Southeast Asia. Although the founder effect contributes significantly to the phylogeographic profile, co-evolution between the host and M. tuberculosis may also play a role. Here, we reported the genomic analysis of 480 L1 isolates from patients in northern Thailand. The studied bacterial population was genetically diverse, allowing the identification of a total of 18 sublineages distributed into three major clades. The majority of isolates belonged to L1.1 followed by L1.2.1 and L1.2.2. Comparison of the single nucleotide variant (SNV) phylogenetic tree and the clades defined by spoligotyping revealed some monophyletic clades representing EAI2_MNL, EAI2_NTM and EAI6_BGD1 spoligotypes. Our work demonstrates that ambiguity in spoligotype assignment could be partially resolved if the entire DR region is investigated. Using the information to map L1 diversity across Southeast Asia highlighted differences in the dominant strain-types in each individual country, despite extensive interactions between populations over time. This finding supported the hypothesis that there is co-evolution between the bacteria and the host, and have implications for tuberculosis disease control.

Spoligotyping based genetic diversity of Mycobacterium tuberculosis in Ethiopia: A systematic review

Background: Understanding the types of strains and lineages of Mycobacterium tuberculosis (M. tuberculosis) circulating in a country is of paramount importance for tuberculosis (TB) control program of that country. Thus, the main aim of this research is to review and compile the results of studies conducted on strains and lineages of M. tuberculosis in Ethiopia. Methods & Materials: A systematic search and review of articles published on M. tuberculosis strains and lineages in Ethiopia were made. PubMed and Google Scholar databases were considered for the search while the key words used were M. tuberculosis, molecular epidemiology, molecular typing spoligotyping and Ethiopia. Results: Twenty-one studies were considered in this review and a total of 3,071 M. tuberculosis isolates and 3,067 strains were included. These studies used spoligotyping and identified five lineages including Indo-Ocean, East Asian/Beijing, East African-Indian, Euro-American and Ethiopian in a proportion of 7.1%, 0.2%, 23.0%, 64.8%, and 4.1%, respectively. Thus, Euro-American was the most frequently (64.8%) occurring Lineage while East Asian was the least (0.2%) frequently occurring Lineage in the country. Surprisingly, the Ethiopian Lineage seemed to be localized to northeastern Ethiopia. In addition, the top five clades identified by this review were T, CAS, H, Manu and Ethiopian comprising of 48.0%, 23.0%, 11.0%, 6.0% and 4.1% of the strains, respectively. Furthermore, predominant shared types (spoligotype patterns) identified were SIT149, SIT53, SIT25, SIT37, and SIT21, each consisting of 420, 343, 266, 162 and 102 isolates, respectively, while, on the other hand, 15% of the strains were orphan. Conclusion: According to the summary of the results of this review, diversified strains and lineages of M. tuberculosis were found in Ethiopia, and the frequencies of occurrence of these strains and lineages were variable in different regions of the country. This systematic review is registered in the PRISMA with the registration number of 42017059263.