Significance of the identification in the Horn of Africa of an exceptionally deep branching Mycobacterium tuberculosis clade

Genetic diversity of Mycobacterium tuberculosis strains isolated from spiritual holy water site attendees in Northwest Ethiopia. A cross-sectional study

Background

The genetic diversity of Mycobacterium tuberculosis complex (MTBC) strains was characterized among isolates from individuals with pulmonary tuberculosis (PTB) symptoms attended holy water sites (HWSs) in the Amhara region, Ethiopia.

Methods

A cross-sectional study was done from June 2019 to March 2020 to describe the genetic diversity and drug-resistance profiles of MTBC isolates. Sputum specimens were collected and cultured in the Löwenstein-Jensen culture medium. Line Probe Assay, MTBDRplus VER 2.0, and MTBDRsl VER 2.0 were used to detect first-and second-line anti-TB drug-resistance patterns. A spoligotyping technique was utilized to characterize the genetic diversity. Statistical analysis was performed using STATA 15.

Results

Of 560 PTB-symptomatic participants, 122 (21.8%) were culture-positive cases. Spoligotyping of 116 isolates revealed diverse MTBC sublineages, with four major lineages: Euro-American (EA) (Lineage 4), East-African-Indian (EAI) (Lineage 3), Ethiopian (ETH) (Lineage 7), East Asian (EA) (Lineage 2). The majority (96.6%) of the isolates were EA (lineage 4) and EAI, with proportions of 54.3% and 42.2%, respectively. A total of 31 spoligotype patterns were identified, 26 of which were documented in the SITVIT2 database. Of these, there were 15 unique spoligotypes, while eleven were grouped with 2-17 isolates. SIT149/T3-ETH (n = 17), SIT26/CAS1-DELHI (n = 16), SIT25/CAS1-DELHI (n = 12), and SIT52/T2 (n = 11) spoligotypes were predominant. A rare spoligotype pattern: SIT41/Turkey and SIT1/Beijing, has also been identified in North Shewa. The overall clustering rate of sub-lineages with known SIT was 76.4%.Of the 122 culture-positive isolates tested, 16.4% were resistant to rifampicin (RIF) and/or isoniazid (INH). Multidrug-resistant TB (MDR-TB) was detected in 12.3% of isolates, five of which were fluoroquinolones (FLQs) resistant. SIT149/T3-ETH and SIT21/CAS1-KILI sublineages showed a higher proportion of drug resistance.

Conclusions

Diverse MTBC spoligotypes were identified, with the T and CAS families and EA (lineage 4) predominating. A high prevalence of drug-resistant TB, with SIT149/T3-ETH and CAS1-KILI sublineages comprising a greater share, was observed. A study with large sample size and a sequencing method with stronger discriminatory power is warranted to understand better the genetic diversity of circulating MTBC in this cohort of study, which would help to adopt targeted interventions.

Genomic transmission clusters and circulating lineages of Mycobacterium tuberculosis among refugees residing in refugee camps in Ethiopia

BACKGROUND

Understanding the transmission dynamics of Mycobacterium tuberculosis (Mtb) could benefit the design of tuberculosis (TB) prevention and control strategies for refugee populations. Whole Genome Sequencing (WGS) has not yet been used to document the Mtb transmission dynamics among refugees in Ethiopia. We applied WGS to accurately identify transmission clusters and Mtb lineages among TB cases in refugee camps in Ethiopia.

METHOD AND DESIGN

We conducted a cross-sectional study of 610 refugees in refugee camps in Ethiopia presenting with symptoms of TB. WGS data of 67 isolates was analyzed using the Maximum Accessible Genome for Mtb Analysis (MAGMA) pipeline; iTol and FigTree were used to visualize phylogenetic trees, lineages, and the presence of transmission clusters.

RESULTS

Mtb culture-positive refugees originated from South Sudan (52/67, 77.6%), Somalia (9/67, 13.4%). Eritrea (4/67, 6%), and Sudan (2/67, 3%). The majority (52, 77.6%) of the isolates belonged to Mtb lineage (L) 3, and one L9 was identified from a Somalian refugee. The vast majority (82%) of the isolates were pan-susceptible Mtb, and none were multi-drug-resistant (MDR)-TB. Based on the 5-single nucleotide polymorphisms cutoff, we identified eight potential transmission clusters containing 23.9% of the isolates. Contact investigation confirmed epidemiological links with either family or social interaction within the refugee camps or with neighboring refugee camps.

CONCLUSION

Four lineages (L1, L3, L4, and L9) were identified, with the majority of strains being L3, reflecting the Mtb L3 dominance in South Sudan, where the majority of refugees originated from. Recent transmission among refugees was relatively low (24%), likely due to the short study period. The improved understanding of the Mtb transmission dynamics using WGS in refugee camps could assist in designing effective TB control programs for refugees.

Genetic diversity of Mycobacterium tuberculosis isolates from the central, eastern and southeastern Ethiopia

Introduction

The population structure of Mycobacterium tuberculosis complex (MTBC) in Ethiopia is diverse but dominated by Euro-American (Lineage 4) and East-African-Indian (Lineage 3) lineages. The objective of this study was to describe the genetic diversity of MTBC isolates in Central, Eastern and Southeastern Ethiopia.

Methods

A total of 223 MTBC culture isolates obtained from patients referred to Adama and Harar TB reference laboratories were spoligotyped. Demographic and clinical characteristics were collected.

Results

Six major lineages: Euro-American (Lineage 4), East-African-Indian (Lineage 3), East Asian (Lineage 2), Indo-Oceanic (Lineage 1), Mycobacterium africanum (Lineage 5 and Lineage 6) and Ethiopian (Lineage 7) were identified. The majority (94.6 %) of the isolates were Euro-American and East-African-Indian, with proportions of 75.3 % and 19.3 %, respectively. Overall, 77 different spoligotype patterns were identified of which 42 were registered in the SITVIT2 database. Of these, 27 spoligotypes were unique, while 15 were clustered with 2-49 isolates. SIT149/T3_ETH (n = 49), SIT53/T1 (n = 33), SIT21/CAS1_Kili (n = 24) and SIT41/Turkey (n = 11) were the dominant spoligotypes. A rare Beijing spoligotype pattern, SIT541, has also been identified in Eastern Ethiopia. The overall clustering rate of sub-lineages with known SIT was 71.3 %. Age group (25-34) was significantly associated with clustering.

Conclusion

We found a heterogeneous population structure of MTBC dominated by T and CAS families, and the Euro-American lineage. The identification of the Beijing strain, particularly the rare SIT541 spoligotype in Eastern Ethiopia, warrants a heightened surveillance plan, as little is known about this genotype. A large-scale investigation utilizing a tool with superior discriminatory power, such as whole genome sequencing, is necessary to gain a thorough understanding of the genetic diversity of MTBC in the nation, which would help direct the overall control efforts.

Phylogenetic analysis of prospective M. bovis antigens with the aim of developing candidate vaccines for bovine tuberculosis

BACKGROUND

Bovine Tuberculosis is a respiratory disease caused by the pathogen Mycobacterium bovis (M. bovis) that infects cattle. Though rare, this disease can also affect humans, as well as domestic and wild animals, making it a serious concern. Therefore, searching for alternative and new vaccines with high efficiency and safety is the main goal in bovine tuberculosis prophylaxis. New vaccines, known as vector vaccines, have the potential to become safe and effective alternatives to the traditional BCG vaccine. In this study, two major immunodominant proteins of M. bovis Esat-6 and TB10.4 were utilized to create a vector vaccine for bovine tuberculosis.

METHODS

The Esat-6 and TB10.4 genes were amplified by PCR. The amplified and purified PCR products were sequenced by the Sanger method. Assembly and multiple alignments of amplicon nucleotides were carried out in the MEGA 11 software.

RESULT

Two genes of the local strain 0078-M. bovis-8/RIBSP were sequenced. The nucleotide sequences were deposited in the GenBank database. Comparative analysis of the nucleotide sequences of the ESAT-6 and TB10.4 genes established 100% identity of the compared strains of Mycobacterium.

CONCLUSION

Through the use of phylogenetic analysis, it has been confirmed that the amplified genes are related to the mycobacteria genus. This discovery allows the development of a vector vaccine against bovine tuberculosis utilising these genes.

The Mycobacterium tuberculosis genome at 25 years: lessons and lingering questions

First achieved in 1998 by Cole et al., the complete genome sequence of Mycobacterium tuberculosis continues to provide an invaluable resource to understand tuberculosis (TB), the leading cause of global infectious disease mortality. At the 25-year anniversary of this accomplishment, we describe how insights gleaned from the M. tuberculosis genome have led to vital tools for TB research, epidemiology, and clinical practice. The increasing accessibility of whole-genome sequencing across research and clinical settings has improved our ability to predict antibacterial susceptibility, to track epidemics at the level of individual outbreaks and wider historical trends, to query the efficacy of the bacille Calmette-Guérin (BCG) vaccine, and to uncover targets for novel antitubercular therapeutics. Likewise, we discuss several recent efforts to extract further discoveries from this powerful resource.

Lineage classification and antitubercular drug resistance surveillance of Mycobacterium tuberculosis by whole-genome sequencing in Southern India

Whole-genome sequencing has created a revolution in tuberculosis management by providing a comprehensive picture of the various genetic polymorphisms with unprecedented accuracy. Studies mapping genomic heterogeneity in clinical isolates of Mycobacterium tuberculosis using a whole-genome sequencing approach from high tuberculosis burden countries are underrepresented. We report whole-genome sequencing results of 242 clinical isolates of culture-confirmed M. tuberculosis isolates from tuberculosis patients referred to a tertiary care hospital in Southern India. Phylogenetic analysis revealed that the isolates in our study belonged to five different lineages, with Indo-Oceanic (lineage 1, n = 122) and East-African Indian (lineage 3, n = 80) being the most prevalent. We report several mutations in genes conferring resistance to first and second line antitubercular drugs including the genes rpoB, katG, ahpC, inhA, fabG1, embB, pncA, rpsL, rrs, and gyrA. The majority of these mutations were identified in relatively high proportions in lineage 1. Our study highlights the utility of whole-genome sequencing as a potential supplemental tool to the existing genotypic and phenotypic methods, in providing expedited comprehensive surveillance of mutations that may be associated with antitubercular drug resistance as well as lineage characterization of M. tuberculosis isolates. Further larger-scale whole-genome datasets with linked minimum inhibition concentration testing are imperative for resolving the discrepancies between whole-genome sequencing and phenotypic drug sensitivity testing results and quantifying the level of the resistance associated with the mutations for optimization of antitubercular drug and precise dose selection in clinics. IMPORTANCE Studies mapping genetic heterogeneity of clinical isolates of M. tuberculosis for determining their strain lineage and drug resistance by whole-genome sequencing are limited in high tuberculosis burden settings. We carried out whole-genome sequencing of 242 M. tuberculosis isolates from drug-sensitive and drug-resistant tuberculosis patients, identified and collected as part of the TB Portals Program, to have a comprehensive insight into the genetic diversity of M. tuberculosis in Southern India. We report several genetic variations in M. tuberculosis that may confer resistance to antitubercular drugs. Further wide-scale efforts are required to fully characterize M. tuberculosis genetic diversity at a population level in high tuberculosis burden settings for providing precise tuberculosis treatment.

Lack of methoxy-mycolates characterizes the geographically restricted lineage 7 of Mycobacterium tuberculosis complex

Lineage 7 (L7) emerged in the phylogeny of the Mycobacterium tuberculosis complex (MTBC) subsequent to the branching of 'ancient' lineage 1 and prior to the Eurasian dispersal of 'modern' lineages 2, 3 and 4. In contrast to the major MTBC lineages, the current epidemiology suggests that prevalence of L7 is highly confined to the Ethiopian population, or when identified outside of Ethiopia, it has mainly been in patients of Ethiopian origin. To search for microbiological factors that may contribute to its restricted distribution, we compared the genome of L7 to the genomes of globally dispersed MTBC lineages. The frequency of predicted functional mutations in L7 was similar to that documented in other lineages. These include mutations characteristic of modern lineages - such as constitutive expression of nitrate reductase - as well as mutations in the VirS locus that are commonly found in ancient lineages. We also identified and characterized multiple lineage-specific mutations in L7 in biosynthesis pathways of cell wall lipids, including confirmed deficiency of methoxy-mycolic acids due to a stop-gain mutation in the mmaA3 gene that encodes a methoxy-mycolic acid synthase. We show that the abolished biosynthesis of methoxy-mycolates of L7 alters the cell structure and colony morphology on selected growth media and impacts biofilm formation. The loss of these mycolic acid moieties may change the host-pathogen dynamic for L7 isolates, explaining the limited geographical distribution of L7 and contributing to further understanding the spread of MTBC lineages across the globe.

Zoonotic Tuberculosis: A Neglected Disease in the Middle East and North Africa (MENA) Region

Mycobacterium bovis is the etiologic agent of bovine tuberculosis (BTB), a serious infectious disease in both humans and animals. BTB is a zoonotic disease primarily affecting cattle and occasionally humans infected through close contact with infected hosts or the consumption of unpasteurized dairy products. Zoonotic tuberculosis is strongly associated with poverty and poor hygiene, and low- and middle-income countries bear the brunt of the disease. BTB has been increasingly recognized as a growing public health threat in developing countries. However, the lack of effective surveillance programs in many of these countries poses a barrier to accurately determining the true burden of this disease. Additionally, the control of BTB is threatened by the emergence of drug-resistant strains that affect the effectiveness of current treatment regimens. Here, we analyzed current trends in the epidemiology of the disease as well as the antimicrobial susceptibility patterns of M. bovis in the Middle East and North Africa (MENA) region, a region that includes several developing countries. Following PRISMA guidelines, a total of 90 studies conducted in the MENA region were selected. Our findings revealed that the prevalence of BTB among humans and cattle varied significantly according to the population size and country in the MENA region. Most of the available studies were based on culture and/or PCR strategies and were published without including data on antimicrobial resistance and molecular typing. Our findings highlighted the paramount need for the use of appropriate diagnostic tools and the implementation of sustainable control measures, especially at the human/animal interface, in the MENA region.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

GWAS and functional studies suggest a role for altered DNA repair in the evolution of drug resistance in Mycobacterium tuberculosis

The emergence of drug resistance in Mycobacterium tuberculosis (Mtb) is alarming and demands in-depth knowledge for timely diagnosis. We performed genome-wide association analysis using 2237 clinical strains of Mtb to identify novel genetic factors that evoke drug resistance. In addition to the known direct targets, we identified for the first time, a strong association between mutations in DNA repair genes and the multidrug-resistant phenotype. To evaluate the impact of variants identified in the clinical samples in the evolution of drug resistance, we utilized knockouts and complemented strains in Mycobacterium smegmatis and Mtb. Results show that variant mutations compromised the functions of MutY and UvrB. MutY variant showed enhanced survival compared with wild-type (Rv) when the Mtb strains were subjected to multiple rounds of ex vivo antibiotic stress. In an in vivo guinea pig infection model, the MutY variant outcompeted the wild-type strain. We show that novel variant mutations in the DNA repair genes collectively compromise their functions and contribute to better survival under antibiotic/host stress conditions.

Brucella abortus in Kazakhstan, population structure and comparison with worldwide genetic diversity

Brucella abortus is the main causative agent of brucellosis in cattle, leading to severe economic consequences in agriculture and affecting public health. The zoonotic nature of the infection increases the need to control the spread and dynamics of outbreaks in animals with the incorporation of high resolution genotyping techniques. Based on such methods, B. abortus is currently divided into three clades, A, B, and C. The latter includes subclades C1 and C2. This study presents the results of whole-genome sequencing of 49 B. abortus strains isolated in Kazakhstan between 1947 and 2015 and of 36 B. abortus strains of various geographic origins isolated from 1940 to 2004. In silico Multiple Locus Sequence Typing (MLST) allowed to assign strains from Kazakhstan to subclades C1 and to a much lower extend C2. Whole-genome Single-Nucleotide Polymorphism (wgSNP) analysis of the 46 strains of subclade C1 with strains of worldwide origins showed clustering with strains from neighboring countries, mostly North Caucasia, Western Russia, but also Siberia, China, and Mongolia. One of the three Kazakhstan strains assigned to subclade C2 matched the B. abortus S19 vaccine strain used in cattle, the other two were genetically close to the 104 M vaccine strain. Bayesian phylodynamic analysis dated the introduction of B. abortus subclade C1 into Kazakhstan to the 19th and early 20th centuries. We discuss this observation in view of the history of population migrations from Russia to the Kazakhstan steppes.

A convolutional neural network highlights mutations relevant to antimicrobial resistance in Mycobacterium tuberculosis

Long diagnostic wait times hinder international efforts to address antibiotic resistance in M. tuberculosis. Pathogen whole genome sequencing, coupled with statistical and machine learning models, offers a promising solution. However, generalizability and clinical adoption have been limited by a lack of interpretability, especially in deep learning methods. Here, we present two deep convolutional neural networks that predict antibiotic resistance phenotypes of M. tuberculosis isolates: a multi-drug CNN (MD-CNN), that predicts resistance to 13 antibiotics based on 18 genomic loci, with AUCs 82.6-99.5% and higher sensitivity than state-of-the-art methods; and a set of 13 single-drug CNNs (SD-CNN) with AUCs 80.1-97.1% and higher specificity than the previous state-of-the-art. Using saliency methods to evaluate the contribution of input sequence features to the SD-CNN predictions, we identify 18 sites in the genome not previously associated with resistance. The CNN models permit functional variant discovery, biologically meaningful interpretation, and clinical applicability.

The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis

Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.

Understanding the Genetic Diversity of Mycobacterium africanum Using Phylogenetics and Population Genomics Approaches

A total of two lineages of Mycobacterium tuberculosis var. africanum (Maf), L5 and L6, which are members of the Mycobacterium tuberculosis complex (MTBC), are responsible for causing tuberculosis in West Africa. Regions of difference (RDs) are usually used for delineation of MTBC. With increased data availability, single nucleotide polymorphisms (SNPs) promise to provide better resolution. Publicly available 380 Maf samples were analyzed for identification of “core-cluster-specific-SNPs,” while additional 270 samples were used for validation. RD-based methods were used for lineage-assignment, wherein 31 samples remained unidentified. The genetic diversity of Maf was estimated based on genome-wide SNPs using phylogeny and population genomics approaches. Lineage-based clustering (L5 and L6) was observed in the whole genome phylogeny with distinct sub-clusters. Population stratification using both model-based and de novo approaches supported the same observations. L6 was further delineated into three sub-lineages (L6.1–L6.3), whereas L5 was grouped as L5.1 and L5.2 based on the occurrence of RD711. L5.1 and L5.2 were further divided into two (L5.1.1 and L5.1.2) and four (L5.2.1–L5.2.4) sub-clusters, respectively. Unassigned samples could be assigned to definite lineages/sub-lineages based on clustering observed in phylogeny along with high-confidence posterior membership scores obtained during population stratification. Based on the (sub)-clusters delineated, “core-cluster-specific-SNPs” were derived. Synonymous SNPs (137 in L5 and 128 in L6) were identified as biomarkers and used for validation. Few of the cluster-specific missense variants in L5 and L6 belong to the central carbohydrate metabolism pathway which include His6Tyr (Rv0946c), Glu255Ala (Rv1131), Ala309Gly (Rv2454c), Val425Ala and Ser112Ala (Rv1127c), Gly198Ala (Rv3293) and Ile137Val (Rv0363c), Thr421Ala (Rv0896), Arg442His (Rv1248c), Thr218Ile (Rv1122), and Ser381Leu (Rv1449c), hinting at the differential growth attenuation. Genes harboring multiple (sub)-lineage-specific “core-cluster” SNPs such as Lys117Asn, Val447Met, and Ala455Val (Rv0066c; icd2) present across L6, L6.1, and L5, respectively, hinting at the association of these SNPs with selective advantage or host-adaptation. Cluster-specific SNPs serve as additional markers along with RD-regions for Maf delineation. The identified SNPs have the potential to provide insights into the genotype–phenotype correlation and clues for endemicity of Maf in the African population.

High-Resolution Melting PCR as Rapid Genotyping Tool for Brucella Species

Brucella sp. are the causative agents of brucellosis. One of the main characteristics of the Brucella genus concerns its very high genetic homogeneity. To date, classical bacteriology typing is still considered as the gold standard assay for direct diagnosis of Brucella. Molecular approaches are routinely used for the identification of Brucella at the genus level. However, genotyping is more complex, and to date, no method exists to quickly assign a strain into species and biovar levels, and new approaches are required. Next generation sequencing (NGS) opened a new era into the diagnosis of bacterial diseases. In this study, we designed a high-resolution melting (HRM) method for the rapid screening of DNA and direct assignment into one of the 12 species of the Brucella genus. This method is based on 17 relevant single nucleotide polymorphisms (SNPs), identified and selected from a whole genome SNP (wgSNP) analysis based on 988 genomes (complete and drafts). These markers were tested against the collection of the European Reference Laboratory (EU-RL) for brucellosis (1440 DNAs extracted from Brucella strains). The results confirmed the reliability of the panel of 17 SNP markers, allowing the differentiation of each species of Brucella together with biovars 1, 2, and 3 of B. suis and vaccine strain Rev1 (B. melitensis) within 3 h, which is a considerable gain of time for brucellosis diagnosis. Therefore, this genotyping tool provides a new and quick alternative for Brucella identification based on SNPs with the HRM-PCR assay.

Biodiversity of Mycobacterium tuberculosis in Bulgaria Related to Human Migrations or Ecological Adaptation

Bulgaria is among the 18 high-priority countries of the WHO European Region with high rates of tuberculosis. The causative agent of tuberculosis is thought to have emerged in Africa 70,000 years ago, or during the Neolithic age, and colonized the world through human migrations. The established main lineages of tuberculosis correlate highly with geography. The goal of our study was to investigate the biodiversity of Mycobacteriumtuberculosis in Bulgaria in association with human migration history during the last 10 centuries. We analyzed spoligotypes and MIRU-VNTR genotyping data of 655 drug-sensitive and 385 multidrug-resistant M. tuberculosis strains collected in Bulgaria from 2008 to 2018. We assigned the genotype of all isolates using SITVITWEB and MIRU-VNTRplus databases and software. We investigated the major well-documented historical events of immigration to Bulgaria that occurred during the last millennium. Genetic profiles demonstrated that, with the exceptions of 3 strains of Mycobacterium bovis and 18 strains of Lineage 2 (W/Beijing spoligotype), only Lineage 4 (Euro-American) was widely diffused in Bulgaria. Analysis of well-documented immigrations of Roma from the Indian subcontinent during the 10th to the 12th centuries, Turkic peoples from Central Asia in the medieval centuries, and more recently Armenians, Russians, and Africans in the 20th century influenced the biodiversity of M. tuberculosis in Bulgaria but only with genotypes of sublineages within the L4. We hypothesize that these sublineages were more virulent, or that ecological adaptation of imported M. tuberculosis genotypes was the main driver contributing to the current genetic biodiversity of M. tuberculosis in Bulgaria. We also hypothesize that some yet unknown local environmental factors may have been decisive in the success of imported genotypes. The ecological factors leading to local genetic biodiversity in M. tuberculosis are multifactorial and have not yet been fully clarified. The coevolution of long-lasting pathogen hosts should be studied, taking into account environmental and ecological changes.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

Mycobacterium tuberculosis strain lineage in mixed tribal population across India and Andaman Nicobar Island

In India, the tribal population constitutes almost 8.6% of the nation's total population. Despite their large presence, there are only a few reports available on Mycobacterium tuberculosis (M. tb) strain prevalence in Indian tribal communities considering the mobile nature of this population and also the influence of the mainstream populations they coexist within many areas for their livelihood. This study attempts to provide critical information pertaining to the TB strain diversity, its public health implications, and distribution among the tribal population in eleven Indian states and Andaman & Nicobar (A&N) Island. The study employed a population-based molecular approach. Clinical isolates were received from 66 villages (10 states and Island) and these villages were selected by implying situation analysis. A total of 78 M. tb clinical isolates were received from 10 different states and A&N Island. Among these, 16 different strains were observed by spoligotyping technique. The major M. tb strains spoligotype belong to the Beijing, CAS1_DELHI, and EAI5 family of M. tb strains followed by EAI1_SOM, EAI6_BGD1, LAM3, LAM6, LAM9, T1, T2, U strains. Drug-susceptibility testing (DST) results showed almost 15.4% of clinical isolates found to be resistant to isoniazid (INH) or rifampicin (RMP) + INH. Predominant multidrug-resistant (MDR-TB) isolates seem to be Beijing strain. Beijing, CAS1_DELHI, EAI3_IND, and EAI5 were the principal strains infecting mixed tribal populations across India. Despite the small sample size, this study has demonstrated higher diversity among the TB strains with significant MDR-TB findings. Prevalence of Beijing MDR-TB strains in Central, Southern, Eastern India and A&N Island indicates the transmission of the TB strains.

Mycobacterial Lineages Associated with Drug Resistance in Patients with Extrapulmonary Tuberculosis in Addis Ababa, Ethiopia

BACKGROUND

In Ethiopia, tuberculosis (TB) is one of the most common causes of illness and death. However, there is limited information available on lineages associated with drug resistance among extrapulmonary tuberculosis patients in Ethiopia. In this study, researchers looked into Mycobacterium tuberculosis lineages linked to drug resistance in patients with extrapulmonary tuberculosis in Addis Ababa, Ethiopia.

METHODS

On 151 Mycobacterium tuberculosis isolates, a cross-sectional analysis was performed. Spoligotyping was used to characterize mycobacterial lineages, while a phenotypic drug susceptibility test was performed to determine the drug resistance pattern. Data were analyzed using SPSS version 23.

RESULTS

Among 151 Mycobacterium tuberculosis complex (MTBC) genotyped isolates, four lineages (L1-L4), and Mycobacterium bovis were identified. The predominantly identified lineage was Euro-American (73.5%) followed by East-African-Indian (19.2%). Any drug resistance (RR) and multidrug-resistant (MDR) tuberculosis was identified among 16.2% and 7.2% of the Euro-American lineage, respectively, while it was 30.8% and 15.4% among the East-African-Indian lineages. Among all three preextensively drug-resistance (pre-XDR) cases identified, two isolates belong to T3-ETH, and the other one strain was not defined by the database. There was no statistically significant association between any type of drug resistance and either lineage or sublineages of Mycobacterium tuberculosis.

CONCLUSION

A higher proportion of any type of drug resistance and MDR was detected among the East-African-Indian lineage compared to others. However, there was no statistically significant association between any type of drug resistance and either lineages or sublineages. Thus, the authors recommend a large-scale study.

Whole genome sequencing of clinical samples reveals extensively drug resistant tuberculosis (XDR TB) strains from the Beijing lineage in Nigeria, West Africa

Multi-drug (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) continues to be a global public health problem especially in high TB burden countries like Nigeria. Many of these cases are undetected and go on to infect high risk individuals. Clinical samples from positive rifampicin resistant Xpert®MTB/Rif assay were subjected to direct whole genome sequencing and bioinformatics analysis to identify the full antibiotics resistance and lineage profile. We report two (2) XDR TB samples also belonging to the East-Asian/Beijing family of lineage 2 Mycobacterium tuberculosis complex from clinical samples in Nigeria. Our findings further reveal the presence of mutations that confer resistance to first-line drugs (rifampicin, isoniazid, ethambutol and pyrazanimide), second-line injectables (capreomycin, streptomycin, kanamycin and/or amikacin) and at least one of the fluoroquinolones (ofloxacin, moxifloxacin, levofloxacin and/or ciprofloxacin) in both samples. The genomic sequence data from this study not only provide the first evidence of XDR TB in Nigeria and West Africa, but also emphasize the importance of WGS in accurately detecting MDR and XDR TB, to ensure adequate and proper management treatment regimens for affected individuals. This will greatly aid in preventing the spread of drug resistance TB in high burden countries.

The Relevance of Genomic Epidemiology for Control of Tuberculosis in West Africa

Tuberculosis (TB), an airborne infectious disease caused by Mycobacterium tuberculosis complex (MTBC), remains a global health problem. West Africa has a unique epidemiology of TB that is characterized by medium- to high-prevalence. Moreover, the geographical restriction of M. africanum to the sub-region makes West Africa have an extra burden to deal with a two-in-one pathogen. The region is also burdened with low case detection, late reporting, poor treatment adherence leading to development of drug resistance and relapse. Sporadic studies conducted within the subregion report higher burden of drug resistant TB (DRTB) than previously thought. The need for more sensitive and robust tools for routine surveillance as well as to understand the mechanisms of DRTB and transmission dynamics for the design of effective control tools, cannot be overemphasized. The advancement in molecular biology tools including traditional fingerprinting and next generation sequencing (NGS) technologies offer reliable tools for genomic epidemiology. Genomic epidemiology provides in-depth insight of the nature of pathogens, circulating strains and their spread as well as prompt detection of the emergence of new strains. It also offers the opportunity to monitor treatment and evaluate interventions. Furthermore, genomic epidemiology can be used to understand potential emergence and spread of drug resistant strains and resistance mechanisms allowing the design of simple but rapid tools. In this review, we will describe the local epidemiology of MTBC, highlight past and current investigations toward understanding their biology and spread as well as discuss the relevance of genomic epidemiology studies to TB control in West Africa.

Genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan

Tularemia is a highly dangerous zoonotic infection due to the bacteria Francisella tularensis. Low genetic diversity promoted the use of polymorphic tandem repeats (MLVA) as first-line assay for genetic description. Whole genome sequencing (WGS) is becoming increasingly accessible, opening the perspective of a time when WGS might become the universal genotyping assay. The main goal of this study was to describe F. tularensis strains circulating in Kazakhstan based on WGS data and develop a MLVA assay compatible with in vitro and in silico analysis. In vitro MLVA genotyping and WGS were performed for the vaccine strain and for 38 strains isolated in Kazakhstan from natural water bodies, ticks, rodents, carnivores, and from one migratory bird, an Isabellina wheatear captured in a rodent burrow. The two genotyping approaches were congruent and allowed to attribute all strains to two F. tularensis holarctica lineages, B.4 and B.12. The seven tandem repeats polymorphic in the investigated strain collection could be typed in a single multiplex PCR assay. Identical MLVA genotypes were produced by in vitro and in silico analysis, demonstrating full compatibility between the two approaches. The strains from Kazakhstan were compared to all publicly available WGS data of worldwide origin by whole genome SNP (wgSNP) analysis. Genotypes differing at a single SNP position were collected within a time interval of more than fifty years, from locations separated from each other by more than one thousand kilometers, supporting a role for migratory birds in the worldwide spread of the bacteria.

Genotyping of Francisella tularensis has become a routine practice in epidemiology. Despite rapidly accumulating knowledge, the phylogeography of the pathogen is still poorly understood and discussions about geographic and temporal origins continue. One important reason is the poor characterization of the pathogen in many tularemia-endemic countries. This article describes the genetic diversity of Francisella tularensis subsp. holarctica in Kazakhstan using tandem repeat polymorphisms as well as whole genome sequencing. Thirty-nine strains were analyzed and two lineages were identified, namely B.4 and B.12. The study demonstrates a wider distribution of genotype B.4 in Asia, and identified a more basal branching point in this subclade. The obtained data support the Asian origin hypothesis for F. tularensis. The finding of identical genotypes in strains separated in time by decades and a thousand-kilometers geographic distance, confirms the ability of the bacteria for long-term preservation and fast long distances spread. The isolation of F. tularensis subsp. holarctica from the bird species Isabellina wheatear allows speculating about a major contribution of birds to the phylogeography of the pathogen. A genotyping protocol was developed utilizing seven polymorphic tandem repeats, two of which were identified within the framework of this work. The in vitro and in silico results are identical when using sequencing reads of 300 base-pairs or more.

Drug Resistance Conferring Mutation and Genetic Diversity of Mycobacterium tuberculosis Isolates in Tuberculosis Lymphadenitis Patients; Ethiopia

Background

Tuberculosis lymphadenitis (TBLN) is a growing public health concern in Ethiopia. However, there is limited information available on gene mutations conferring drug resistance and genetic diversity of M. tuberculosis isolates from TBLN patients.

Methods

Drug resistance and genetic diversity analysis were done on 91 M. tuberculosis isolates from culture positive TBLN patients collected between 2016 and 2017. Detection of mutations conferring resistance was carried out using GenoType MTBDRplus VER 2.0. Thereafter, isolates were typed using spoligotyping.

Results

Out of the 91 strains, mutations conferring resistance to rifampicin (RIF) and isoniazid (INH) were observed in two (2.2%) and six (6.6%) isolates, respectively. The two RIF resistant isolates displayed a mutation at codon 531 in the rpoB gene with amino acid change of S531L. Among the six INH resistant strains, four isolates had shown mutation at the KatG gene at codon 315 with amino acid change of S315T, one isolate had a mutation at the inhA gene at codon 15 with amino acid change of C15T and one isolate had a mutation at the inhA gene with unknown amino acid change. All drug resistant isolates were from treatment naive TBLN patients. The dominantly identified Spoligo International Types (SITs) were SIT25, SIT149, and SIT53, respectively; these accounted for 43% of the total number of strains. The isolates were grouped into four main lineages; Lineage 1 (2, 2.2%), Lineage 3 (38, 41.7%), Lineage 4 (49, 53.8%) and Lineage 7 (2, 2.2%). Four out of six (66.7%) isolates with drug resistance conferring mutations belonged to clustered strains (strains with shared SIT).

Conclusion

The detection of drug resistant conferring mutation in treatment naïve TBLN patients together with detection of drug resistant isolates among clustered strains might suggest resistant strains' transmission in the community. This needs to be carefully considered to prevent the spread of drug resistant clones in the country.

Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum. Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum, and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.

The Epidemiology of first and second-line drug-resistance Mycobacterium tuberculosis complex common species: Evidence from selected TB treatment initiating centers in Ethiopia

Background

Drug-resistance in Mycobacterium tuberculosis complex remains a major health burden in human history and still is a major leading cause of death in developing countries including Ethiopia. Early detection of all forms of drug-resistant Tuberculosis(TB) is a key factor to reduce and contain the spread of these resistant strains.

Methods

A health facility-based cross-sectional study was employed, based on demographic, clinical, and laboratory data collected from 204 patients with bacteriological confirmed TB. Sputum samples were analyzed using conventional TB culture and identification test followed by molecular species identification, and then phenotypic drug susceptibility tests. Data were entered using an excel spreadsheet and exported to SPSS version 20 for analysis. Descriptive analysis; frequencies, and proportions were computed.

Results

Among the 204 sputum samples inoculated in culture media, Mycobacterium species were recovered from 165 specimens, with 160 Mycobacterium tuberculosis complex and five Non- Tuberculosis Mycobacterium(NTM) species. All Mycobacterium tuberculosis complex was found to be M. tuberculosis. Of the five NTM species, 2 M.fortuitum, 2 M.intracellulare, and 1 M.gordonae were identified. Among 160 species of M. tuberculosis isolates, 110(68.8%) were resistant to any of the anti-TB drugs. The resistance pattern was; INH (109, 68.1%), RIF (99, 61.9%), STM (73,45.6%), and EMB (32,20.0%). Mono-resistance was found for INH (7,4.3%) and STM (1,0.6%). Ninety-nine (61.9%) isolates become MDR, while resistance to any of the second-line anti-TB drugs was detected in 9 (5.6%) strains, with 8(5%) Pre-XDR and one (0.6%) XDR cases.

Conclusion

Our findings highlight high frequencies of drug resistance to first and second-line anti-TB drugs.Determining the drug-resistance pattern of MTB is important for programmatic management of drug-resistant TB in Ethiopia. The circulating Pre-XDR and XDR case identified in the current study is alarming to the tuberculosis control program in the country.

Retrospective Analysis of the Relationship between Two Anthrax Outbreaks in Kazakhstan Based on Genomic Data

We present a retrospective analysis of strains from two anthrax outbreaks in western Kazakhstan in 2009. The outbreaks occurred during the same period and in the same area located close to main roads, favoring a single source of infection. However, multilocus variable-number tandem-repeat analysis (MLVA), canonical single-nucleotide polymorphism (CanSNP) analysis, and genome-wide analysis demonstrated that the outbreaks were not connected.

We present a retrospective analysis of strains from two anthrax outbreaks in western Kazakhstan in 2009. The outbreaks occurred during the same period and in the same area located close to main roads, favoring a single source of infection. However, multilocus variable-number tandem-repeat analysis (MLVA), canonical single-nucleotide polymorphism (CanSNP) analysis, and genome-wide analysis demonstrated that the outbreaks were not connected.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

BACKGROUND

Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes.

RESULTS

We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats.

CONCLUSION

This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

Population structure of multidrug-resistant Mycobacterium tuberculosis clinical isolates in Colombia

Emergence of multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) isolates is a major public health problem that threatens progress made in tuberculosis (TB) care and control worldwide. In Colombia, the prevalence of MDR tuberculosis (MDR-TB) has increased slowly but steadily since 2001. However, the population structure of the MDR-TB strains circulating in Colombia is sparsely known. In this work, 203 MDR isolates isolated in 2012-2013 were collected, and characterized by spoligotyping, followed by 24-loci MIRU-VNTR (data available for 190 isolates). The most prevalent genotypes corresponded to SIT42/LAM9 (12.81%), SIT62/H1 (10.34%), and SIT190/Beijing (10.34%). A fine analysis showed that although the MDR strains came from 29 of the 33 departments of Colombia, the distribution of these main lineages was not at random and depended on the city of isolation (p-value <0.000001). Both LAM and Beijing lineage strains were significantly associated with MDR-TB (p-value <0.0001): LAM lineage was associated with 2 patterns of MDR, namely combined resistance to INH + Rifampin (HR), and to SHRE (Streptomycin + INH + Rifampin + Ethambutol), while the Beijing lineage strains were essentially associated with MDR (SHRE). Interestingly, distribution of genotypic lineages in function of drug resistance information (e.g. pansusceptible vs. MDR) was different in our setting as compared to other countries in Latin America. However, MIRU-VNTR patterns were unique for all strains, an observation that did not support active transmission of circulating MDR clones.

The whole genome sequence data analyses of a Mycobacterium tuberculosis strain SBH321 isolated in Sabah, Malaysia, belongs to Ural family of Lineage 4

In 2019, 10 million new cases of tuberculosis have been reported worldwide. Our data reports genetic analyses of a Mycobacterium tuberculosis strain SBH321 isolated from a 31-year-old female with pulmonary tuberculosis. The genomic DNA of the strain was extracted from pure culture and subjected to sequencing using Illumina platform. M. tuberculosis strain SBH321 consists of 4,374,895 bp with G+C content of 65.59%. The comparative analysis by SNP-based phylogenetic analysis using maximum-likelihood method showed that our strain belonging to sublineage of the Ural family of Europe-America-Africa lineage (Lineage 4) and clustered with M. tuberculosis strain OFXR-4 from Taiwan. The whole genome sequence is deposited at DDBJ/ENA/GenBank under the accession WCJH00000000 (SRR10230353).

Whole Genome Sequencing and Spatial Analysis Identifies Recent Tuberculosis Transmission Hotspots in Ghana

Whole genome sequencing (WGS) is progressively being used to investigate the transmission dynamics of Mycobacterium tuberculosis complex (MTBC). We used WGS analysis to resolve traditional genotype clusters and explored the spatial distribution of confirmed recent transmission clusters. Bacterial genomes from a total of 452 MTBC isolates belonging to large traditional clusters from a population-based study spanning July 2012 and December 2015 were obtained through short read next-generation sequencing using the illumina HiSeq2500 platform. We performed clustering and spatial analysis using specified R packages and ArcGIS. Of the 452 traditional genotype clustered genomes, 314 (69.5%) were confirmed clusters with a median cluster size of 7.5 genomes and an interquartile range of 4–12. Recent tuberculosis (TB) transmission was estimated as 24.7%. We confirmed the wide spread of a Cameroon sub-lineage clone with a cluster size of 78 genomes predominantly from the Ablekuma sub-district of Accra metropolis. More importantly, we identified a recent transmission cluster associated with isoniazid resistance belonging to the Ghana sub-lineage of lineage 4. WGS was useful in detecting unsuspected outbreaks; hence, we recommend its use not only as a research tool but as a surveillance tool to aid in providing the necessary guided steps to track, monitor, and control TB.

Molecular epidemiology of drug resistant Mycobacterium tuberculosis in Africa: a systematic review

Background

The burden of drug resistant tuberculosis in Africa is largely driven by the emergence and spread of multidrug resistant (MDR) and extensively drug resistant (XDR) Mycobacterium tuberculosis strains. MDR-TB is defined as resistance to isoniazid and rifampicin, while XDR-TB is defined as MDR-TB with added resistance to any of the second line injectable drugs and any fluoroquinolone.

The highest burden of drug resistant TB is seen in countries further experiencing an HIV epidemic. The molecular mechanisms of drug resistance as well as the evolution of drug resistant TB strains have been widely studied using various genotyping tools. The study aimed to analyse the drug resistant lineages in circulation and transmission dynamics of these lineages in Africa by describing outbreaks, nosocomial transmission and migration. Viewed as a whole, this can give a better insight into the transmission dynamics of drug resistant TB in Africa.

Methods

A systematic review was performed on peer reviewed original research extracted from PubMed reporting on the lineages associated with drug resistant TB from African countries, and their association with outbreaks, nosocomial transmission and migration. The search terms “Tuberculosis AND drug resistance AND Africa AND (spoligotyping OR molecular epidemiology OR IS6110 OR MIRU OR DNA fingerprinting OR RFLP OR VNTR OR WGS)” were used to identify relevant articles reporting the molecular epidemiology of drug resistant TB in Africa.

Results

Diverse genotypes are associated with drug resistant TB in Africa, with variations in strain predominance within the continent. Lineage 4 predominates across Africa demonstrating the ability of “modern strains” to adapt and spread easily. Most studies under review reported primary drug resistance as the predominant type of transmission. Drug resistant TB strains are associated with community and nosocomial outbreaks involving MDR- and XDR-TB strains. The under-use of molecular epidemiological tools is of concern, resulting in gaps in knowledge of the transmission dynamics of drug resistant TB on the continent.

Conclusions

Genetic diversity of M. tuberculosis strains has been demonstrated across Africa implying that diverse genotypes are driving the epidemiology of drug resistant TB across the continent.

TbD1 deletion as a driver of the evolutionary success of modern epidemic Mycobacterium tuberculosis lineages

Mycobacterium tuberculosis (Mtb) strains are classified into different phylogenetic lineages (L), three of which (L2/L3/L4) emerged from a common progenitor after the loss of the MmpS6/MmpL6-encoding Mtb-specific deletion 1 region (TbD1). These TbD1-deleted “modern” lineages are responsible for globally-spread tuberculosis epidemics, whereas TbD1-intact “ancestral” lineages tend to be restricted to specific geographical areas, such as South India and South East Asia (L1) or East Africa (L7). By constructing and characterizing a panel of recombinant TbD1-knock-in and knock-out strains and comparison with clinical isolates, here we show that deletion of TbD1 confers to Mtb a significant increase in resistance to oxidative stress and hypoxia, which correlates with enhanced virulence in selected cellular, guinea pig and C3HeB/FeJ mouse infection models, the latter two mirroring in part the development of hypoxic granulomas in human disease progression. Our results suggest that loss of TbD1 at the origin of the L2/L3/L4 Mtb lineages was a key driver for their global epidemic spread and outstanding evolutionary success.

Mycobacterium tuberculosis (Mtb) modern strains emerged from a common progenitor after the loss of Mtb-specific deletion 1 region (TbD1). Here, the authors show that deletion of TbD1 correlates with enhanced Mtb virulence in animal models, mirroring the development of hypoxic granulomas in human disease progression.

Evolution of Mycobacterium tuberculosis complex lineages and their role in an emerging threat of multidrug resistant tuberculosis in Bamako, Mali

In recent years Bamako has been faced with an emerging threat from multidrug resistant TB (MDR-TB). Whole genome sequence analysis was performed on a subset of 76 isolates from a total of 208 isolates recovered from tuberculosis patients in Bamako, Mali between 2006 and 2012. Among the 76 patients, 61(80.3%) new cases and 15(19.7%) retreatment cases, 12 (16%) were infected by MDR-TB. The dominant lineage was the Euro-American lineage, Lineage 4. Within Lineage 4, the Cameroon genotype was the most prevalent genotype (n = 20, 26%), followed by the Ghana genotype (n = 16, 21%). A sub-clade of the Cameroon genotype, which emerged ~22 years ago was likely to be involved in community transmission. A sub-clade of the Ghana genotype that arose approximately 30 years ago was an important cause of MDR-TB in Bamako. The Ghana genotype isolates appeared more likely to be MDR than other genotypes after controlling for treatment history. We identified a clade of four related Beijing isolates that included one MDR-TB isolate. It is a major concern to find the Cameroon and Ghana genotypes involved in community transmission and MDR-TB respectively. The presence of the Beijing genotype in Bamako remains worrying, given its high transmissibility and virulence.

Identification and Characterization of Mycobacterium tuberculosis Beijing Genotype Strain SBH163, Isolated in Sabah, Malaysia

This is a report on the whole-genome sequence of Mycobacterium tuberculosis strain SBH163, which was isolated from a patient in the Malaysian Borneo state of Sabah. This report provides insight into the molecular characteristics of an M. tuberculosis Beijing genotype strain related to strains from Russia and South Africa.

An ancestral genomic locus in Mycobacterium tuberculosis clinical isolates from India hints the genetic link with Mycobacterium canettii

BACKGROUND

Tuberculosis remains a worldwide public health emergency. To better understand M. tuberculosis and to identify genomic variations characteristic to the Indian clinical isolates by a low-cost method, a genomic subtractive hybridization between M. tuberculosis H37Rv and a clinical isolate from South India was performed.

RESULTS

This revealed a novel 0.4-kb subtractive fragment which was used as a handle to pull out a 4.5-kb genomic region characteristic to the clinical isolate and was absent in H37Rv. On further studies, this 4.5-kb region was found to be present in 91% of the M. tuberculosis clinical isolates screened from Kerala, a state in South India. Interestingly, this novel region has 99% identity (with 100% query coverage) with genomic regions of M. canettii.

DISCUSSION

The present study hypothesizes that this locus was present in the recent common environmental ancestor of mycobacteria, retained to the maximum extent in M. canettii and ancestral isolates of M. tuberculosis, and later deleted in other modern lineages of M. tuberculosis. Thus, this region may serve as one of the links between the pathogenic mycobacteria and the environmental species. We also propose that the Indian isolates of M. tuberculosis might be closely related to the putative progenitor M. prototuberculosis with respect to this locus. More studies on other genomic loci from different strains of M. tuberculosis are required to establish more links in this direction.

Whole genome sequencing data and analysis of a rifampicin-resistant Mycobacterium tuberculosis strain SBH162 from Sabah, Malaysia

A Mycobacterium tuberculosis strain SBH162 was isolated from a 49-year-old male with pulmonary tuberculosis. GeneXpert MDR/RIF identified the strain as rifampicin-resistant M. tuberculosis. The whole genome sequencing was performed using Illumina HiSeq 4000 system to further investigate and verify the mutation sites of the strain through genetic analyses namely variant calling using bioinformatics tools. The de novo assembly of genome generated 100 contigs with N50 of 156,381bp. The whole genome size was 4,343,911 bp with G + C content of 65.58% and consisted of 4,306 predicted genes. The mutation site, S450L, for rifampicin resistance was detected in the rpoB gene. Based on the phylogenetic analysis using the Maximum Likelihood method, the strain was identified as belonging to the Europe America Africa lineage (Lineage 4). The genome dataset has been deposited at DDBJ/ENA/GenBank under the accession number SMOE00000000.

Towards better prediction of Mycobacterium tuberculosis lineages from MIRU-VNTR data

The determination of lineages from strain-based molecular genotyping information is an important problem in tuberculosis. Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing is a commonly used molecular genotyping approach that uses counts of the number of times pre-specified loci repeat in a strain. There are three main approaches for determining lineage based on MIRU-VNTR data - one based on a direct comparison to the strains in a curated database, and two others, on machine learning algorithms trained on a large collection of labeled data. All existing methods have limitations. The direct approach imposes an arbitrary threshold on how much a database strain can differ from a given one to be informative. On the other hand, the machine learning-based approaches require a substantial amount of labeled data. Notably, all three methods exhibit suboptimal classification accuracy without additional data. We explore several computational approaches to address these limitations. First, we show that eliminating the arbitrary threshold improves the performance of the direct approach. Second, we introduce RuleTB, an alternative direct method that proposes a concise set of rules for determining lineages. Lastly, we propose StackTB, a machine learning approach that requires only a fraction of the training data to outperform the accuracy of both existing machine learning methods. Our approaches demonstrate superior performance on a training dataset collected in New York City over 10 years, and the improvement in performance translates to a held-out testing set. We conclude that our methods provide opportunities for improving the determination of pathogenic lineages based on MIRU-VNTR data.

Lineage specific histories of Mycobacterium tuberculosis dispersal in Africa and Eurasia

Mycobacterium tuberculosis (M.tb) is a globally distributed, obligate pathogen of humans that can be divided into seven clearly defined lineages. An emerging consensus places the origin and global dispersal of M.tb within the past 6,000 years: identifying how the ancestral clone of M.tb spread and differentiated within this timeframe is important for identifying the ecological drivers of the current pandemic. We used Bayesian phylogeographic inference to reconstruct the migratory history of M.tb in Africa and Eurasia and to investigate lineage specific patterns of spread from a geographically diverse sample of 552 M.tb genomes. Applying evolutionary rates inferred with ancient M.tb genome calibration, we estimated the timing of major events in the migratory history of the pathogen. Inferred timings contextualize M.tb dispersal within historical phenomena that altered patterns of connectivity throughout Africa and Eurasia: trans-Indian Ocean trade in spices and other goods, the Silk Road and its predecessors, the expansion of the Roman Empire, and the European Age of Exploration. We found that Eastern Africa and Southeast Asia have been critical in the dispersal of M.tb. Our results further reveal that M.tb populations have grown through range expansion, as well as in situ, and delineate the independent evolutionary trajectories of bacterial subpopulations underlying the current pandemic.

Insights from Bacillus anthracis strains isolated from permafrost in the tundra zone of Russia

This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common feature of these strains is their conservation in permafrost, from which they were extracted either due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In contrast, during the excavations in Yakutia, three genetically different strains were recovered from a single pit. One strain belongs to B.Br.001/002, and whole genome sequence analysis showed that it is most closely related to the Yamal strains in spite of the remoteness of Yamal from Yakutia. The two other strains contribute to two different branches of A.Br.008/011, one of the remarkable polytomies described so far in the B. anthracis species. The geographic distribution of the strains belonging to A.Br.008/011 is suggesting that the polytomy emerged in the thirteenth century, in combination with the constitution of a unified Mongol empire extending from China to Eastern Europe. We propose an evolutionary model for B. anthracis recent evolution in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage except in some geographically isolated areas.

Insights from Bacillus anthracis strains isolated from permafrost in the tundra zone of Russia

This article describes Bacillus anthracis strains isolated during an outbreak of anthrax on the Yamal Peninsula in the summer of 2016 and independently in Yakutia in 2015. A common feature of these strains is their conservation in permafrost, from which they were extracted either due to the thawing of permafrost (Yamal strains) or as the result of paleontological excavations (Yakut strains). All strains isolated on the Yamal share an identical genotype belonging to lineage B.Br.001/002, pointing to a common source of infection in a territory over 250 km in length. In contrast, during the excavations in Yakutia, three genetically different strains were recovered from a single pit. One strain belongs to B.Br.001/002, and whole genome sequence analysis showed that it is most closely related to the Yamal strains in spite of the remoteness of Yamal from Yakutia. The two other strains contribute to two different branches of A.Br.008/011, one of the remarkable polytomies described so far in the B. anthracis species. The geographic distribution of the strains belonging to A.Br.008/011 is suggesting that the polytomy emerged in the thirteenth century, in combination with the constitution of a unified Mongol empire extending from China to Eastern Europe. We propose an evolutionary model for B. anthracis recent evolution in which the B lineage spread throughout Eurasia and was subsequently replaced by the A lineage except in some geographically isolated areas.

Accurate Phylogenetic Relationships Among Mycobacterium bovis Strains Circulating in France Based on Whole Genome Sequencing and Single Nucleotide Polymorphism Analysis

In recent years the diversity of the French Mycobacterium bovis population responsible for bovine tuberculosis (bTB) outbreaks since 1970 has been described in detail. To further understand bTB evolution in France, we used single nucleotide polymorphisms (SNPs) based on whole genome sequence versus classical genotyping methods in order to identify accurate phylogenetic relationships between M. bovis strains. Whole genome sequencing was carried out on a selection of 87 strains which reflect the French M. bovis population’s genetic diversity. Sequences were compared to the M. bovis reference genome AF2122/97. Comparison among the 87 genomes revealed 9,170 sites where at least one strain shows a SNP with respect to the reference genome; 1,172 are intergenic and 7,998 in coding sequences, of which 2,880 are synonymous and 5,118 non-synonymous. SNP-based phylogenetic analysis using these 9,170 SNP is congruent with the cluster defined by spoligotyping and multilocus variable number of tandem repeat analysis typing. In addition, some SNPs were identified as specific to genotypic groups. These findings suggest new SNP targets that can be used for the development of high-resolving methods for genotyping as well as for studying M. bovis evolution and transmission patterns. The detection of non-synonymous SNPs on virulence genes enabled us to distinguish different clusters. Our results seem to indicate that genetically differentiated clusters could also display distinctive phenotypic traits.

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.

Pulmonary Isolation of Multidrug resistant "Mycobacterium simulans" and Mycobacterium tuberculosis from a patient in the Horn of Africa

In low-income countries of the Horn of Africa, pulmonary infections are usually considered as tuberculosis, which diagnosis relies on clinical data and positive microscopic observation. This strategy allows non-tuberculous mycobacteria to escape detection, facilitating their emergence in populations. A non-tuberculous mycobacterium strain FB-527 was unexpectedly cultured from the sputum of a Djiboutian patient otherwise diagnosed with multi-drug resistant (MDR) tuberculosis. The sequencing of the rpoB and 16S rRNA genes showed that the isolate was identical to strain FI-09026 previously named "Mycobacterium simulans" and reported only once from a Somali patient. Strain FB-527 mimicked Mycobacterium tuberculosis colonies and enzymatic profile using API ZYM strip and was in vitro resistant to rifampicin and isoniazid. Isolation of two MDR mycobacteria complicated the diagnosis and therapeutic management of the patient. We here report on the complete description of strain FB-527 and strain FI-09026 including genome sequencing, finalizing the description of the proposed new species "Mycobacterium simulans".

Unexpected Genomic and Phenotypic Diversity of Mycobacterium africanum Lineage 5 Affects Drug Resistance, Protein Secretion, and Immunogenicity

Mycobacterium africanum consists of Lineages L5 and L6 of the Mycobacterium tuberculosis complex (MTBC) and causes human tuberculosis in specific regions of Western Africa, but is generally not transmitted in other parts of the world. Since M. africanum is evolutionarily closely placed between the globally dispersed Mycobacterium tuberculosis and animal-adapted MTBC-members, these lineages provide valuable insight into M. tuberculosis evolution. Here, we have collected 15 M. africanum L5 strains isolated in France over 4 decades. Illumina sequencing and phylogenomic analysis revealed a previously underappreciated diversity within L5, which consists of distinct sublineages. L5 strains caused relatively high levels of extrapulmonary tuberculosis and included multi- and extensively drug-resistant isolates, especially in the newly defined sublineage L5.2. The specific L5 sublineages also exhibit distinct phenotypic characteristics related to in vitro growth, protein secretion and in vivo immunogenicity. In particular, we identified a PE_PGRS and PPE-MPTR secretion defect specific for sublineage L5.2, which was independent of PPE38. Furthermore, L5 isolates were able to efficiently secrete and induce immune responses against ESX-1 substrates contrary to previous predictions. These phenotypes of Type VII protein secretion and immunogenicity provide valuable information to better link genome sequences to phenotypic traits and thereby understand the evolution of the MTBC.

Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa

Objective

To investigate the distribution of Mycobacterium tuberculosis genotypes across Africa.

Methods

The SITVIT2 global repository and PUBMED were searched for spoligotype and published genotype data respectively, of M. tuberculosis from Africa. M. tuberculosis lineages in Africa were described and compared across regions and with those from 7 European and 6 South-Asian countries. Further analysis of the major lineages and sub-lineages using Principal Component analysis (PCA) and hierarchical cluster analysis were done to describe clustering by geographical regions. Evolutionary relationships were assessed using phylogenetic tree analysis.

Results

A total of 14727 isolates from 35 African countries were included in the analysis and of these 13607 were assigned to one of 10 major lineages, whilst 1120 were unknown. There were differences in geographical distribution of major lineages and their sub-lineages with regional clustering. Southern African countries were grouped based on high prevalence of LAM11-ZWE strains; strains which have an origin in Portugal. The grouping of North African countries was due to the high percentage of LAM9 strains, which have an origin in the Eastern Mediterranean region. East African countries were grouped based on Central Asian (CAS) and East-African Indian (EAI) strain lineage possibly reflecting historic sea trade with Asia, while West African Countries were grouped based on Cameroon lineage of unknown origin. A high percentage of the Haarlem lineage isolates were observed in the Central African Republic, Guinea, Gambia and Tunisia, however, a mixed distribution prevented close clustering.

Conclusions

This study highlighted that the TB epidemic in Africa is driven by regional epidemics characterized by genetically distinct lineages of M. tuberculosis. M. tuberculosis in these regions may have been introduced from either Europe or Asia and has spread through pastoralism, mining and war. The vast array of genotypes and their associated phenotypes should be considered when designing future vaccines, diagnostics and anti-TB drugs.

Reduced transmission of Mycobacterium africanum compared to Mycobacterium tuberculosis in urban West Africa

OBJECTIVE

Understanding transmission dynamics is useful for tuberculosis (TB) control. A population-based molecular epidemiological study was conducted to determine TB transmission in Ghana.

METHODS

Mycobacterium tuberculosis complex (MTBC) isolates obtained from prospectively sampled pulmonary TB patients between July 2012 and December 2015 were characterized using spoligotyping and standard 15-locus mycobacterial interspersed repetitive unit variable number tandem repeat (MIRU-VNTR) typing for transmission studies.

RESULTS

Out of 2309 MTBC isolates, 1082 (46.9%) unique cases were identified, with 1227 (53.1%) isolates belonging to one of 276 clusters. The recent TB transmission rate was estimated to be 41.2%. Whereas TB strains of lineage 4 belonging to M. tuberculosis showed a high recent transmission rate (44.9%), reduced recent transmission rates were found for lineages of Mycobacterium africanum (lineage 5, 31.8%; lineage 6, 24.7%).

CONCLUSIONS

The study findings indicate high recent TB transmission, suggesting the occurrence of unsuspected outbreaks in Ghana. The observed reduced transmission rate of M. africanum suggests other factor(s) (host/environmental) may be responsible for its continuous presence in West Africa.

Microbead-based spoligotyping of Mycobacterium tuberculosis from Ziehl-Neelsen-stained microscopy preparations in Ethiopia

The worldwide dissemination of Mycobacterium tuberculosis strains has led to the study of their genetic diversity. One of the most used genotyping methods is spoligotyping, based on the detection of spacers in the clustered regularly interspaced short palindromic repeats (CRISPR) locus. This study assessed the performance of a microbead-based spoligotyping assay using samples extracted from Ziehl-Neelsen-stained smear-microscopy preparations and described the genetic diversity of Mycobacterium tuberculosis among new TB patients in Southern Nations, Nationalities and Peoples’ Region (SNNPR) in Ethiopia. Among the 91 samples analysed, 59 (64.8%) generated spoligotyping patterns. Fifty (84.7%) samples were classified into 12 clusters (mostly Lineage 4 or 3) comprising 2–11 samples and nine had unique spoligotyping patterns. Among the 59 spoligotyping patterns, 25 belonged to the T1 sublineage, 11 to the T3-ETH, 5 to the URAL, 4 to the H3 and 14 to other L4 sublineages. There was a remarkable variation in genetic distribution in SNNPR compared to other regions of the country. Microbead-based spoligotyping is an easy-to-perform, high-throughput assay that can generate genotyping information using material obtained from smear microscopy preparations. The method provides an opportunity to obtain data of the M. tuberculosis genetic epidemiology in settings with limited laboratory resources.

Mycobacterium tuberculosis and whole-genome sequencing: how close are we to unleashing its full potential?

BACKGROUND

Nearly two decades after completion of the genome sequence of Mycobacterium tuberculosis (MTB), and with the advent of next generation sequencing technologies (NGS), whole-genome sequencing (WGS) has been applied to a wide range of clinical scenarios. Starting in 2017, England is the first country in the world to pioneer its use on a national scale for the diagnosis of tuberculosis, detection of drug resistance, and typing of MTB.

AIMS

This narrative review critically analyses the current applications of WGS for MTB and explains how close we are to realizing its full potential as a diagnostic, epidemiologic, and research tool.

SOURCES

We searched for reports (both original articles and reviews) published in English up to 31 May 2017, with combinations of the following keywords: whole-genome sequencing, Mycobacterium, and tuberculosis. MEDLINE, Embase, and Scopus were used as search engines. We included articles that covered different aspects of whole-genome sequencing in relation to MTB.

CONTENT

This review focuses on three main themes: the role of WGS for the prediction of drug susceptibility, MTB outbreak investigation and genetic diversity, and research applications of NGS.

IMPLICATIONS

Many of the original expectations have been accomplished, and we believe that with its unprecedented sensitivity and power, WGS has the potential to address many unanswered questions in the near future. However, caution is still needed when interpreting WGS data as there are some important limitations to be aware of, from correct interpretation of drug susceptibilities to the bioinformatic support needed.

The predominance of Ethiopian specific Mycobacterium tuberculosis families and minimal contribution of Mycobacterium bovis in tuberculous lymphadenitis patients in Southwest Ethiopia

BACKGROUND

Ethiopia has an extremely high rate of extrapulmonary tuberculosis, dominated by tuberculous lymphadenitis (TBLN). However, little is known about Mycobacterium tuberculosis complex (MTBc) lineages responsible for TBLN in Southwest Ethiopia.

METHODS

A total of 304 MTBc isolates from TBLN patients in Southwest Ethiopia were genotyped primarily by spoligotyping. Isolates of selected spoligotypes were further analyzed by 15-loci mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) (n=167) and qPCR-based single nucleotide polymorphism (n=38). Isolates were classified into main phylogenetic lineages and families by using the reference strain collections and identification tools available at MIRU-VNTRplus data base. Resistance to rifampicin was determined by Xpert MTB/RIF.

RESULTS

The majority of isolates (248; 81.6%) belonged to the Euro-American lineage (Lineage 4), with the ill-defined T and Haarlem as largest families comprising 116 (38.2%) and 43 (14.1%) isolates respectively. Of the T family, 108 isolates were classified as being part of the newly described Ethiopian families, namely Ethiopia_2 (n=44), Ethiopia_3 (n=34) and Ethiopia_H₃₇Rv-like (n=30). Other sub-lineages included URAL (n=18), S (n=17), Uganda I (n=16), LAM (n=13), X (n=5), TUR (n=5), Uganda II (n=4) and unknown (n=19). Lineage 3 (Delhi/CAS) was the second most common lineage comprising 44 (14.5%) isolates. Interestingly, six isolates (2%) were belonged to Lineage 7, unique to Ethiopia. Lineage 1 (East-African Indian) and Lineage 2 (Beijing) were represented by 3 and 1 isolates respectively. M. bovis was identified in only two (0.7%) TBLN cases. The cluster rate was highest for Ethiopia_3 isolates showing clonal similarity with isolates from North Ethiopia. Lineage 3 was significantly associated with rifampicin resistance.

CONCLUSIONS

In TBLN in Southwest Ethiopia, the recently described Ethiopia specific Lineage 4 families were predominant, followed by Lineage 3 and Lineage 4-Haarlem. The contribution of M. bovis in TBLN infection is minimal.