Population Genomics of Mycobacterium tuberculosis in Ethiopia Contradicts the Virgin Soil Hypothesis for Human Tuberculosis in Sub-Saharan Africa

Whole-genome sequencing-based genetic diversity, transmission dynamics, and drug-resistant mutations in Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

Background

Extrapulmonary tuberculosis (EPTB) refers to a form of Tuberculosis (TB) where the infection occurs outside the lungs. Despite EPTB being a devastating disease of public health concern, it is frequently overlooked as a public health problem. This study aimed to investigate genetic diversity, identify drug-resistance mutations, and trace ongoing transmission chains.

Methods

A cross-sectional study was undertaken on individuals with EPTB in western Ethiopia. In this study, whole-genome sequencing (WGS) was employed to analyze Mycobacterium tuberculosis (MTB) samples obtained from EPTB patients. Out of the 96 genomes initially sequenced, 89 met the required quality standards for genetic diversity, and drug-resistant mutations analysis. The data were processed using robust bioinformatics tools.

Results

Our analysis reveals that the majority (87.64%) of the isolates can be attributed to Lineage-4 (L4), with L4.6.3 and L4.2.2.2 emerging as the predominant sub-lineages, constituting 34.62% and 26.92%, respectively. The overall clustering rate and recent transmission index (RTI) were 30 and 17.24%, respectively. Notably, 7.87% of the isolates demonstrated resistance to at least one anti-TB drug, although multi-drug resistance (MDR) was observed in only 1.12% of the isolates.

Conclusions

The genetic diversity of MTBC strains in western Ethiopia was found to have low inter-lineage diversity, with L4 predominating and exhibiting high intra-lineage diversity. The notably high clustering rate in the region implies a pressing need for enhanced TB infection control measures to effectively disrupt the transmission chain. It's noteworthy that 68.75% of resistance-conferring mutations went undetected by both GeneXpert MTB/RIF and the line probe assay (LPA) in western Ethiopia. The identification of resistance mutations undetected by both GeneXpert and LPA, along with the detection of mixed infections through WGS, emphasizes the value of adopting WGS as a high-resolution approach for TB diagnosis and molecular epidemiological surveillance.

Genomic Interactions Between Mycobacterium tuberculosis and Humans

Mycobacterium tuberculosis is considered by many to be the deadliest microbe, with the estimated annual cases numbering more than 10 million. The bacteria, including Mycobacterium africanum, are classified into nine major lineages and hundreds of sublineages, each with different geographical distributions and levels of virulence. The phylogeographic patterns can be a result of recent and early human migrations as well as coevolution between the bacteria and various human populations, which may explain why many studies on human genetic factors contributing to tuberculosis have not been replicable in different areas. Moreover, several studies have revealed the significance of interactions between human genetic variations and bacterial genotypes in determining the development of tuberculosis, suggesting coadaptation. The increased availability of whole-genome sequence data from both humans and bacteria has enabled a better understanding of these interactions, which can inform the development of vaccines and other control measures.

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.

Molecular epidemiology and drug sensitivity of Mycobacterium tuberculosis in homeless individuals in the Addis Ababa city, Ethiopia

Although homeless segment of the society could be the hotspots for tuberculosis (TB) transmission, there is little data on TB in homeless individuals in Ethiopia. The objective of this study was to investigate the molecular epidemiology and drug sensitivity of Mycobacterium tuberculosis (M. tuberculosis) isolated from homeless individuals in Addis Ababa, Ethiopia. The study was conducted on 59 M. tuberculosis isolates, which were recovered by the clinical screening of 5600 homeless individuals and bacteriological examination of 641 individuals with symptoms of pulmonary tuberculosis (PTB). Region of difference-9 (RD9) based polymerase-chain reaction (PCR), Spoligotyping and 24-loci Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) typing were used for genotyping of the isolates. In addition, drug sensitivity test was performed on the isolates using BD Bactec Mycobacterial Growth Inhibition Tube (MGIT) 960. Fifty-eight of the 59 isolates were positive by spoligotyping and spoligotyping International type (SIT) 53, SIT 37, and SIT 149 were the dominant spoligotypes; each consisting of 19%, 15.5%, and10.3% of the isolates, respectively. The majority of the isolates (89.7%) were members of the Euro-American (EA) major lineage. MIRU-VNTR identified Ethiopia_3, Delhi/CAS, Ethiopia_2, TUR, X-type, Ethiopia_H37Rv-like strain, Haarlem and Latin-American Mediterranean (LAM) sub lineages. The proportion of clustering was 77.6% (45/58) in spoligotyping while it was 39.7% (23/58) in 24-loci MIRU-VNTR typing. Furthermore, the proportion of clustering was significantly lowered to 10.3% (6/58) when a combination of spoligotyping and 24-loci MIRU-VNTRplus was used. The recent transmission index (RTI) recorded by spoligotyping, 24-loci MIRU-VNTR typing, and a combination of the two genotyping methods were 58.6%, 27.6% and 5.2%, respectively. Young age and living in groups were significantly associated with strain clustering (P < 0.05). The drug sensitivity test (DST) result showed 8.9% (4/58) of the isolates were resistant to one or more first line ant-TB drugs; but multidrug resistant isolate was not detected. Clustering and RTI could suggest the transmission of TB in the homeless individuals, which could suggest a similar pattern of transmission between homeless individuals and the general population. Hence, the TB control program should consider homeless individuals during the implementation of TB control program.

Sensitive lipid biomarker detection for tuberculosis in late Neanderthal skeletons from Subalyuk Cave, Hungary

Skeletal remains of two Neanderthal individuals, a 25-35 year-old woman and a 3-4 year-old child, were discovered in a Subalyuk Cave in North-Eastern Hungary. Radiocarbon dating of the female and child remains revealed an age of 39,732-39,076 and 36,117-35,387 cal BP, respectively. Paleopathological studies of these Neanderthal remains revealed probable evidence of skeletal mycobacterial infection, including in the sacrum of the adult specimen and the endocranial surface of the child's skull. Application of PCR amplification to the juvenile cranium and a vertebra gave a positive result (IS6110) for tuberculosis, backed up by spoligotyping. Lipid biomarker analyses of the same two specimens revealed definitive signals for C32 mycoserosates, a very characteristic component of the Mycobacterium tuberculosis complex (MTBC). A vertebra from the adult provided weak evidence for mycocerosate biomarkers. The correlation of probable skeletal lesions with characteristic amplified DNA fragments and a proven lipid biomarker points to the presence of tuberculosis in these Neanderthals. In particular, the closely similar biomarker profiles, for two distinct juvenile cranial and vertebral bones, strengthen this diagnosis.

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.

Mycobacterium tuberculosis Sub-Lineage 4.2.2/SIT149 as Dominant Drug-Resistant Clade in Northwest Ethiopia 2020-2022: In-silico Whole-Genome Sequence Analysis

Introduction

Drug resistance (DR) in Mycobacterium tuberculosis complex (MTBC) is mainly associated with certain lineages and varies across regions and countries. The Beijing genotype is the leading resistant lineage in Asia and western countries. M. tuberculosis (Mtb) (sub) lineages responsible for most drug resistance in Ethiopia are not well described. Hence, this study aimed to identify the leading drug resistance sub-lineages and characterize first-line anti-tuberculosis drug resistance-associated single nucleotide polymorphisms (SNPs).

Methods

A facility-based cross-sectional study was conducted in 2020-2022 among new and presumptive multidrug resistant-TB (MDR-TB) cases in Northwest Ethiopia. Whole-genome sequencing (WGS) was performed on 161 isolates using Illumina NovaSeq 6000 technology. The SNP mutations associated with drug resistance were identified using MtbSeq and TB profiler Bioinformatics softwares.

Results

Of the 146 Mtb isolates that were successfully genotyped, 20 (13.7%) harbored one or more resistance-associated SNPs. L4.2.2.ETH was the leading drug-resistant sub-lineage, accounting for 10/20 (50%) of the resistant Mtb. MDR-TB isolates showed extensive mutations against first-line anti-TB drugs. Ser450Leu/(tcg/tTg) for Rifampicin (RIF), Ser315Thr/(agc/aCc) for Isoniazid (INH), Met306Ile/(atg/atA(C)) for Ethambutol (EMB), and Gly69Asp for Streptomycin (STR) were the leading resistance associated mutations which accounted for 56.5%, 89.5%, 47%, and 29.4%, respectively. The presence of both clustered and non-clustered drug resistance (DR) isolates indicated that the epidemics is driven by both new DR development and acquired resistance.

Conclusion

The high prevalence of drug-resistant TB due to geographically restricted sub-lineages (L4.2.2.ETH) indicates the ongoing local micro epidemics. The Mtb drug resistance surveillance system must be improved. Further evolutionary analysis of L4.2.2.ETH strain is highly desirable to understand evolutionary forces that leads L4.2.2.ETH in to high level DR and transmissible sub-lineage.

Mycobacterium tuberculosis complex drug-resistance, phylogenetics, and evolution in Nigeria: Comparison with Ghana and Cameroon

In this article, we provide an in-depth analysis on the drug-resistance phenotypic characteristics of a cohort of 325 tuberculosis and characterize by Whole Genome Sequencing 24 isolates from Nigeria belonging to L4, L5 and L6. Our results suggest an alarming rate of drug-resistance of the L4.6.2.2 Mycobacterium tuberculosis complex (MTBC) lineage and a high diversity of L5. We compiled these new Sequence Read Archives (SRAs) to previously published ones from available Bioprojects run in Nigeria. We performed RAxML phylogenetic reconstructions of larger samples that include public NCBI SRAs from some neighboring countries (Cameroon, Ghana). To confront phylogenetic reconstruction to metadata, we used a new proprietary database named TB-Annotator. We show that L5 genomes in Northern Nigeria belong to new clades as the ones described until now and allow an update of the taxonomy of L5. In addition, we describe the L4.6.2.2 lineage in Nigeria, Cameroon and Ghana. We provide computations on the likely divergence time of L4.6.2.2 and suggest a new hypothesis concerning its origin. Finally we provide a short overview on M. bovis diversity in Nigeria. This study constitutes a baseline knowledge on the global genomic diversity, phylogeography and phylodynamics of MTBC in Nigeria, as well as on the natural history of this largely ignored but densely populated country of Africa. These results highlight the need of sequencing additional MTBC genomes in Nigeria and more generally in West-Africa, both for public health and for academic reasons. The likelihood of replacement of L5-L6 by L4.6.2.2 isolates, leave potentially little time to gather historical knowledge informative on the ancient history of tuberculosis in West-Africa.

Molecular epidemiology and transmission dynamics of multi-drug resistant tuberculosis strains using whole genome sequencing in the Amhara region, Ethiopia

BACKGROUND

Drug resistant Mycobacterium tuberculosis prevention and care is a major challenge in Ethiopia. The World health organization has designated Ethiopia as one of the 30 high burden multi-drug resistant tuberculosis (MDR-TB) countries. There is limited information regarding genetic diversity and transmission dynamics of MDR-TB in Ethiopia.

OBJECTIVE

To investigate the molecular epidemiology and transmission dynamics of MDR-TB strains using whole genome sequence (WGS) in the Amhara region.

METHODS

Forty-five MDR-TB clinical isolates from Amhara region were collected between 2016 and 2018, and characterized using WGS and 24-loci Mycobacterium Interspersed Repetitive Units Variable Number of Tandem Repeats (MIRU-VNTR) typing. Clusters were defined based on the maximum distance of 12 single nucleotide polymorphisms (SNPs) or alleles as the upper threshold of genomic relatedness. Five or less SNPs or alleles distance or identical 24-loci VNTR typing is denoted as surrogate marker for recent transmission.

RESULTS

Forty-one of the 45 isolates were analyzed by WGS and 44% (18/41) of the isolates were distributed into 4 clusters. Of the 41 MDR-TB isolates, 58.5% were classified as lineage 4, 36.5% lineage 3 and 5% lineage 1. Overall, TUR genotype (54%) was the predominant in MDR-TB strains. 41% (17/41) of the isolates were clustered into four WGS groups and the remaining isolates were unique strains. The predominant cluster (Cluster 1) was composed of nine isolates belonging to lineage 4 and of these, four isolates were in the recent transmission links.

CONCLUSIONS

Majority of MDR-TB strain cluster and predominance of TUR lineage in the Amhara region give rise to concerns for possible ongoing transmission. Efforts to strengthen TB laboratory to advance diagnosis, intensified active case finding, and expanded contact tracing activities are needed in order to improve rapid diagnosis and initiate early treatment. This would lead to the interruption of the transmission chain and stop the spread of MDR-TB in the Amhara region.

Comparative whole-genome sequence analysis of Mycobacterium tuberculosis isolated from pulmonary tuberculosis and tuberculous lymphadenitis patients in Northwest Ethiopia

Background

Tuberculosis (TB), caused by the Mycobacterium tuberculosis complex (MTBC), is a chronic infectious disease with both pulmonary and extrapulmonary forms. This study set out to investigate and compare the genomic diversity and transmission dynamics of Mycobacterium tuberculosis (Mtb) isolates obtained from tuberculous lymphadenitis (TBLN) and pulmonary TB (PTB) cases in Northwest Ethiopia.

Methods

A facility-based cross-sectional study was conducted using two groups of samples collected between February 2021 and June 2022 (Group 1) and between June 2020 and June 2022 (Group 2) in Northwest Ethiopia. Deoxyribonucleic acid (DNA) was extracted from 200 heat-inactivated Mtb isolates. Whole-genome sequencing (WGS) was performed from 161 isolates having ≥1 ng DNA/μl using Illumina NovaSeq 6000 technology.

Results

From the total 161 isolates sequenced, 146 Mtb isolates were successfully genotyped into three lineages (L) and 18 sub-lineages. The Euro-American (EA, L4) lineage was the prevailing (n = 100; 68.5%) followed by Central Asian (CAS, L3, n = 43; 25.3%) and then L7 (n = 3; 2.05%). The L4.2.2.ETH sub-lineage accounted for 19.9%, while Haarlem estimated at 13.7%. The phylogenetic tree revealed distinct Mtb clusters between PTB and TBLN isolates even though there was no difference at lineages and sub-lineages levels. The clustering rate (CR) and recent transmission index (RTI) for PTB were 30 and 15%, respectively. Similarly, the CR and RTI for TBLN were 31.1 and 18 %, respectively.

Conclusion and recommendations

PTB and TBLN isolates showed no Mtb lineages and sub-lineages difference. However, at the threshold of five allelic distances, Mtb isolates obtained from PTB and TBLN form distinct complexes in the phylogenetic tree, which indicates the presence of Mtb genomic variation among the two clinical forms. The high rate of clustering and RTI among TBLN implied that TBLN was likely the result of recent transmission and/or reactivation from short latency. Hence, the high incidence rate of TBLN in the Amhara region could be the result of Mtb genomic diversity and rapid clinical progression from primary infection and/or short latency. To validate this conclusion, a similar community-based study with a large sample size and better sampling technique is highly desirable. Additionally, analysis of genomic variants other than phylogenetic informative regions could give insightful information. Combined analysis of the host and the pathogen genome (GXG) together with environmental (GxGxE) factors could give comprehensive co-evolutionary information.

Rapid Bioethical Assessment for Mycobacterium tuberculosis and Host Genetic Study in Amhara Regional State, Ethiopia: Towards a Context-Specific Ethical Approach

Background

Rapid Ethical Assessment (REA) is a rapid qualitative study anticipated to understand the ethical sphere of the research setting prior to recruiting study subjects. This study assessed the communities' knowledge about tuberculosis (TB) and research, understand the social arrangements advisable for recruiting research participant and appraised the information provision and consent process.

Methods

The study was conducted in Amhara region, Ethiopia from 5th-30th January 2021. Google-based survey, face-to-face in-depth interview and focus group discussion were carried out to collect the data from researchers, data collectors, health professionals, TB program officers. A structured questionnaire was administered to assess the knowledge of TB patients and healthy controls about TB, research, gene, (co)evolution and consent process.

Results

Over 71% of researchers were not satisfied with the current consent process, and 82.7% of researchers agreed that the best interest of the research participants was not adequately addressed in the current research practices in ANRS. TB patients and healthy controls misunderstood research and its goals. Participants advised the researchers to approach the community with the assistance of health extension workers (HEW) or religious/local leaders. Combined use of verbal and written based information provision at individual participant level is the preferred way for information provision.

Conclusions

The adherence of researchers to standard information provision and consent process was very low. Healthy controls and TB patients have low level of knowledge and awareness about research, ethics and genomic research-related common terms. Hence, public education is required to strengthen the research ethics in the region.

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.

An ancestral mycobacterial effector promotes dissemination of infection

The human pathogen Mycobacterium tuberculosis typically causes lung disease but can also disseminate to other tissues. We identified a M. tuberculosis (Mtb) outbreak presenting with unusually high rates of extrapulmonary dissemination and bone disease. We found that the causal strain carried an ancestral full-length version of the type VII-secreted effector EsxM rather than the truncated version present in other modern Mtb lineages. The ancestral EsxM variant exacerbated dissemination through enhancement of macrophage motility, increased egress of macrophages from established granulomas, and alterations in macrophage actin dynamics. Reconstitution of the ancestral version of EsxM in an attenuated modern strain of Mtb altered the migratory mode of infected macrophages, enhancing their motility. In a zebrafish model, full-length EsxM promoted bone disease. The presence of a derived nonsense variant in EsxM throughout the major Mtb lineages 2, 3, and 4 is consistent with a role for EsxM in regulating the extent of dissemination.

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.

Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection

Genetic differences between different Mycobacterium tuberculosis complex (MTBC) strains determine their ability to transmit within different host populations, their latency times, and their drug resistance profiles. Said differences usually emerge through de novo mutations and are maintained or discarded by the balance of evolutionary forces. Using a dataset of ∼5,000 strains representing global MTBC diversity, we determined the past and present selective forces that have shaped the current variability observed in the pathogen population. We identified regions that have evolved under changing types of selection since the time of the MTBC common ancestor. Our approach highlighted striking differences in the genome regions relevant for host–pathogen interaction and, in particular, suggested an adaptive role for the sensor protein of two-component systems. In addition, we applied our approach to successfully identify potential determinants of resistance to drugs administered as second-line tuberculosis treatments.

Whole Genome Sequencing of Drug Resistant and Drug Susceptible Mycobacterium tuberculosis Isolates From Tigray Region, Ethiopia

Background: Tuberculosis, mainly caused by Mycobacterium tuberculosis (Mtb), is an ancient human disease that gravely affects millions of people annually. We wanted to explore the genetic diversity and lineage-specific association of Mtb with drug resistance among pulmonary tuberculosis patients. Methods: Sputum samples were collected from pulmonary tuberculosis patients at six different healthcare institutions in Tigray, Ethiopia, between July 2018 and August 2019. DNA was extracted from 74 Mtb complex isolates for whole-genome sequencing (WGS). All genomes were typed and screened for mutations with known associations with antimicrobial resistance using in silico methods, and results were cross-verified with wet lab methods. Results: Lineage (L) 4 (55.8%) was predominant, followed by L3 (41.2%); L1 (1.5%) and L2 (1.5%) occurred rarely. The most frequently detected sublineage was CAS (38.2%), followed by Ural (29.4%), and Haarlem (11.8%). The recent transmission index (RTI) was relatively low. L4 and Ural strains were more resistant than the other strains to any anti-TB drug (P < 0.05). The most frequent mutations to RIF, INH, EMB, SM, PZA, ETH, FLQs, and 2nd-line injectable drugs occurred at rpoB S450L, katG S315T, embB M306I/V, rpsL K43R, pncA V139A, ethA M1R, gyrA D94G, and rrs A1401G, respectively. Disputed rpoB mutations were also shown in four (16%) of RIF-resistant isolates. Conclusion: Our WGS analysis revealed the presence of diverse Mtb genotypes. The presence of a significant proportion of disputed rpoB mutations highlighted the need to establish a WGS facility at the regional level to monitor drug-resistant mutations. This will help control the transmission of DR-TB and ultimately contribute to the attainment of 100% DST coverage for TB patients as per the End TB strategy.

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.

Molecular Epidemiology of Mycobacterium tuberculosis strains isolated from pulmonary tuberculosis patients in south Ethiopia

Introduction:

Understanding the epidemiology of tuberculosis is limited by lack of genotyping data. We sought to characterize the drug susceptibility testing patterns and genetic diversity of M. tuberculosis isolates in southern Ethiopia.

Methodology:

A cross-sectional study was conducted among newly diagnosed sputum smear positive patients with tuberculosis visiting nine health facilities in southern Ethiopia from June 2015 to May 2016. Three consecutive sputum samples (spot-morning-spot) per patient were examined using acid-fast bacilli smear microscopy with all smear positive specimens having acid-fast bacilli cultures performed. M. tuberculosis isolates had drug susceptibility testing performed using indirect proportion method and were genotyped with RD9 deletion analysis and spoligotyping. Mapping of strain was made using geographic information system.

Results:

Among 250 newly diagnosed patients with tuberculosis, 4% were HIV co-infected. All 230 isolates tested were M. tuberculosis strains belonging to three lineages: Euro-American, 187 (81%), East-African-Indian, 31 (14%), and Lineage 7 (Ethiopian lineage), 8 (4%); categorized into 63 different spoligotype patterns, of which 85% fell into 28 clusters. M. tuberculosis strains were clustered by geographic localities. The dominant spoligotypes were SIT149 (21%) and SIT53 (19%). Drug susceptibility testing found that 14% of isolates tested were resistant to ≥ 1 first line anti- tuberculosis drugs and 11% to INH. SIT 149 was dominant among drug resistant isolates.

Conclusions:

The study revealed several clusters and drug resistant strains of M. tuberculosis in the study area, suggesting recent transmission including of drug resistant tuberculosis. Wider monitoring of drug susceptibility testing and geospatial analysis of transmission trends is required to control tuberculosis in southern Ethiopia.

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.

Diversity of Mycobacterium tuberculosis Complex Lineages Associated with Pulmonary Tuberculosis in Southwestern, Uganda

Uganda is among the 22 countries in the world with a high burden of tuberculosis. The southwestern region of the country has consistently registered a high TB/HIV incidence rate. This study is aimed at characterizing the Mycobacterium tuberculosis complex (MTBC) genotypic diversity in southwestern Uganda. A total of 283 sputum samples from patients with pulmonary tuberculosis were genotyped using specific single nucleotide polymorphism markers for lineages 3 and 4. Most of the patients were males with a mean age of 34. The lineage 4 Ugandan family was found to be the most dominant strains accounting for 59.7% of all cases followed by lineage 3 at 15.2%. The lineage 4 non-Ugandan family accounted for 14.5% of all cases while 4.2% showed amplification for both lineage 4 and lineage 3. Eighteen samples (6.4%) of the strains remained unclassified since they could not be matched to any lineage based on the genotyping technique used. This study demonstrates that a wide diversity of strains is causing pulmonary tuberculosis in this region with those belonging to the lineage 4 Ugandan family being more predominant. However, to confirm this, further studies using more discriminative genotyping methods are necessary.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

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.

Evolutionary Genetics of Mycobacterium tuberculosis and HIV-1: "The Tortoise and the Hare"

The already enormous burden caused by Mycobacterium tuberculosis and Human Immunodeficiency Virus type 1 (HIV-1) alone is aggravated by co-infection. Despite obvious differences in the rate of evolution comparing these two human pathogens, genetic diversity plays an important role in the success of both. The extreme evolutionary dynamics of HIV-1 is in the basis of a robust capacity to evade immune responses, to generate drug-resistance and to diversify the population-level reservoir of M group viral subtypes. Compared to HIV-1 and other retroviruses, M. tuberculosis generates minute levels of genetic diversity within the host. However, emerging whole-genome sequencing data show that the M. tuberculosis complex contains at least nine human-adapted phylogenetic lineages. This level of genetic diversity results in differences in M. tuberculosis interactions with the host immune system, virulence and drug resistance propensity. In co-infected individuals, HIV-1 and M. tuberculosis are likely to co-colonize host cells. However, the evolutionary impact of the interaction between the host, the slowly evolving M. tuberculosis bacteria and the HIV-1 viral "mutant cloud" is poorly understood. These evolutionary dynamics, at the cellular niche of monocytes/macrophages, are also discussed and proposed as a relevant future research topic in the context of single-cell sequencing.

Genetic diversity of the Mycobacterium tuberculosis complex strains from newly diagnosed tuberculosis patients in Northwest Ethiopia reveals a predominance of East-African-Indian and Euro-American lineages

OBJECTIVES

This study described the population structure of M. tuberculosis complex (MTBc) strains among patients with pulmonary or lymph node tuberculosis (TB) in Northwest Ethiopia and tested the performance of culture isolation and MPT64-based speciation for Lineage 7 (L7).

METHODS

Patients were recruited between April 2017 and June 2019 in North Gondar, Ethiopia. The MPT64 assay was used to confirm MTBc, and spoligotyping was used to characterize mycobacterial lineages. Line probe assay (LPA) was used to detect resistance to rifampicin and isoniazid.

RESULTS

Among 274 MTBc genotyped isolates, there were five MTBc lineages: L1-L4 and L7 were identified, with predominant East-African-Indian (L3) (53.6%) and Euro-American (L4) (40.1%) strains, and low prevalence (2.6%) of Ethiopia L7. The genotypes were similarly distributed between pulmonary and lymph node TB, and all lineages were equally isolated by culture and recognized as MTBc by the MPT64 assay. Additionally, LPA showed that 259 (94.5%) MTBc were susceptible to both rifampicin and isoniazid, and one (0.4%) was multi-drug resistant (resistant to both rifampicin and isoniazid).

CONCLUSION

These findings show that TB in North Gondar, Ethiopia, is mainly caused by L3 and L4 strains, with low rates of L7, confirmed as MTBc by MPT64 assay and with limited resistance to rifampicin and isoniazid.

Whole genome sequencing-based analysis of tuberculosis (TB) in migrants: rapid tools for cross-border surveillance and to distinguish between recent transmission in the host country and new importations

Background

The analysis of transmission of tuberculosis (TB) is challenging in areas with a large migrant population. Standard genotyping may fail to differentiate transmission within the host country from new importations, which is key from an epidemiological perspective.

Aim

To propose a new strategy to simplify and optimise cross-border surveillance of tuberculosis and to distinguish between recent transmission in the host country and new importations

Methods

We selected 10 clusters, defined by 24-locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR), from a population in Spain rich in migrants from eastern Europe, north Africa and west Africa and reanalysed 66 isolates by whole-genome sequencing (WGS). A multiplex-allele-specific PCR was designed to target strain-specific marker single nucleotide polymorphisms (SNPs), identified from WGS data, to optimise the surveillance of the most complex cluster.

Results

In five of 10 clusters not all isolates showed the short genetic distances expected for recent transmission and revealed a higher number of SNPs, thus suggesting independent importations of prevalent strains in the country of origin. In the most complex cluster, rich in Moroccan cases, a multiplex allele-specific oligonucleotide-PCR (ASO-PCR) targeting the marker SNPs for the transmission subcluster enabled us to prospectively identify new secondary cases. The ASO-PCR-based strategy was transferred and applied in Morocco, demonstrating that the strain was prevalent in the country.

Conclusion

We provide a new model for optimising the analysis of cross-border surveillance of TB transmission in the scenario of global migration.

Dissemination of VIM-producing Pseudomonas aeruginosa associated with high-risk clone ST654 in a tertiary and quaternary hospital in Makkah, Saudi Arabia

To the best of our knowledge, no molecular surveillance that has been conducted to identify the most common clones of carbapenem-resistant Pseudomonas aeruginosa (CRPA) in western Saudi Arabia. Therefore, this study aimed to identify genetic diversity and the most common CRPA clones in this region. Thirty-five CRPA isolates were collected from a tertiary and quaternary hospital in Makkah. bla _VIM was the most common carbapenemase-encoding gene (11 CRPA isolates), while bla_GES was reported in only three isolates. CRPA isolates were subjected to multi- locus sequence typing and showed relatively high genetic diversity with 20 sequence types. Approximately one-third (31.4%) of the CRPA isolates belonged to two high-risk clones (ST235 and ST654). This troublesome finding raises serious concerns about the emergence and further dissemination of CRPA high-risk clones in local hospitals and suggest that surveillance programs should be established in this region to monitor and control clonal dissemination of all multidrug resistant bacteria, including CRPA.

A sister lineage of the Mycobacterium tuberculosis complex discovered in the African Great Lakes region

The human- and animal-adapted lineages of the Mycobacterium tuberculosis complex (MTBC) are thought to have expanded from a common progenitor in Africa. However, the molecular events that accompanied this emergence remain largely unknown. Here, we describe two MTBC strains isolated from patients with multidrug resistant tuberculosis, representing an as-yet-unknown lineage, named Lineage 8 (L8), seemingly restricted to the African Great Lakes region. Using genome-based phylogenetic reconstruction, we show that L8 is a sister clade to the known MTBC lineages. Comparison with other complete mycobacterial genomes indicate that the divergence of L8 preceded the loss of the cobF genome region - involved in the cobalamin/vitamin B12 synthesis - and gene interruptions in a subsequent common ancestor shared by all other known MTBC lineages. This discovery further supports an East African origin for the MTBC and provides additional molecular clues on the ancestral genome reduction associated with adaptation to a pathogenic lifestyle.

The human- and animal-adapted lineages of the Mycobacterium tuberculosis complex (MTBC) are thought to be evolved from a common progenitor in Africa. Here, the authors identify two MTBC strains isolated from patients with multidrug-resistant tuberculosis, representing an as-yet-unknown lineage further supporting an East African origin for the MTBC.

Mycobacterium bovis: From Genotyping to Genome Sequencing

Mycobacterium bovis is the main pathogen of bovine, zoonotic, and wildlife tuberculosis. Despite the existence of programs for bovine tuberculosis (bTB) control in many regions, the disease remains a challenge for the veterinary and public health sectors, especially in developing countries and in high-income nations with wildlife reservoirs. Current bTB control programs are mostly based on test-and-slaughter, movement restrictions, and post-mortem inspection measures. In certain settings, contact tracing and surveillance has benefited from M. bovis genotyping techniques. More recently, whole-genome sequencing (WGS) has become the preferential technique to inform outbreak response through contact tracing and source identification for many infectious diseases. As the cost per genome decreases, the application of WGS to bTB control programs is inevitable moving forward. However, there are technical challenges in data analyses and interpretation that hinder the implementation of M. bovis WGS as a molecular epidemiology tool. Therefore, the aim of this review is to describe M. bovis genotyping techniques and discuss current standards and challenges of the use of M. bovis WGS for transmission investigation, surveillance, and global lineages distribution. We compiled a series of associated research gaps to be explored with the ultimate goal of implementing M. bovis WGS in a standardized manner in bTB control programs.

Mycobacterium tuberculosis associated with severe tuberculosis evades cytosolic surveillance systems and modulates IL-1β production

Genetic diversity of Mycobacterium tuberculosis affects immune responses and clinical outcomes of tuberculosis (TB). However, how bacterial diversity orchestrates immune responses to direct distinct TB severities is unknown. Here we study 681 patients with pulmonary TB and show that M. tuberculosis isolates from cases with mild disease consistently induce robust cytokine responses in macrophages across multiple donors. By contrast, bacteria from patients with severe TB do not do so. Secretion of IL-1β is a good surrogate of the differences observed, and thus to classify strains as probable drivers of different TB severities. Furthermore, we demonstrate that M. tuberculosis isolates that induce low levels of IL-1β production can evade macrophage cytosolic surveillance systems, including cGAS and the inflammasome. Isolates exhibiting this evasion strategy carry candidate mutations, generating sigA recognition boxes or affecting components of the ESX-1 secretion system. Therefore, we provide evidence that M. tuberculosis strains manipulate host-pathogen interactions to drive variable TB severities.

Genetic diversity and drug resistance pattern of Mycobacterium tuberculosis strains isolated from pulmonary tuberculosis patients in the Benishangul Gumuz region and its surroundings, Northwest Ethiopia

Introduction

Tuberculosis (TB) remains a major global public health problem and is the leading cause of death from a single bacterium, Mycobacterium tuberculosis (MTB) complex. The emergence and spread of drug-resistant strains aggravate the problem, especially in tuberculosis high burden countries such as Ethiopia. The supposedly high initial cost of laboratory diagnosis coupled with scarce financial resources has limited collection of information about drug resistance patterns and circulating strains in peripheral and emerging regions of Ethiopia. Here, we investigated drug susceptibility and genetic diversity of mycobacterial isolates among pulmonary tuberculosis patients in the Benishangul Gumuz region and its surroundings in northwest Ethiopia.

Methods and material

In a cross-sectional study, 107 consecutive sputum smear-positive pulmonary tuberculosis (PTB) patients diagnosed at two hospitals and seven health centers were enrolled between October 2013 and June 2014. Sputum samples were cultured at Armauer Hansen Research Institute (AHRI) TB laboratory, and drug susceptibility testing (DST) was performed against Isoniazid, Rifampicin, Ethambutol, and Streptomycin using the indirect proportion method. Isolates were characterized using polymerase chain reaction (PCR)based Region of Difference 9 (RD9) testing and spoligotyping. Statistical analysis was performed using Statistical Package for the Social Sciences (SPSS) for Windows version 24.0.

Results

Of 107 acid-fast-bacilli (AFB) smear-positive sputum samples collected, 81.3% (87/107) were culture positive. A PCR based RD9 testing revealed that all the 87 isolates were M. tuberculosis. Of these isolates, 16.1% (14/87) resistance to one or more drugs was observed. Isoniazid monoresistance occurred in 6.9% (6/87). Multidrug resistance (MDR) was observed in two isolates (2.3%), one of which was resistant to all the four drugs tested. Spoligotyping revealed that the majority, 61.3% (46/75) of strains could be grouped into ten spoligotype patterns containing two to 11 isolates each while the remaining 38.7% (29/75) were unique. SIT289 (11 isolates) and SIT53 (nine isolates) constituted 43.5% (20/46) among clustered isolates while 29.3% (22/75) were ‘‘New” to the database. The dominant families were T, 37% (28/75), CAS, 16.0% (12/75), and H, 8% (6/75), adding up to 51.3% (46/75) of all isolates identified.

Conclusion and recommendations

The current study indicates a moderate prevalence of MDR TB. However, the observed high monoresistance to Isoniazid, one of the two proxy drugs for MDR-TB, reveals the hidden potential threat fora sudden increase in MDR-TB if resistance to Rifampicin would increase. Clustered spoligotype patterns suggest ongoing active tuberculosis transmission in the area. The results underscore the need for enhanced monitoring of TB drug resistance and epidemiological studies in this and other peripheral regions of the country using robust molecular tools with high discriminatory power such as the Mycobacterial Interspersed Repetitive Units -Variable Number of Tandem Repeats (MIRU-VNTR) typing and whole-genome sequencing (WGS).

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.

Tuberculosis in Liberia: high multidrug-resistance burden, transmission and diversity modelled by multiple importation events

Tuberculosis (TB) surveillance is scarce in most African countries, even though it is the continent with the greatest disease incidence according to the World Health Organization. Liberia is within the 30 countries with the highest TB burden, probably as a consequence of the long civil war and the recent Ebola outbreak, both crippling the health system and depreciating the TB prevention and control programmes. Due to difficulties working in the country, there is a lack of resistance surveys and bacillus characterization. Here, we use genome sequencing of Mycobacteriumtuberculosis clinical isolates to fill this gap. Our results highlight that the bacillus population structure is dominated by lineage 4 strains that harbour an outstanding genetic diversity, higher than in the rest of Africa as a whole. Coalescent analyses demonstrate that strains currently circulating in Liberia were introduced several times beginning in the early year 600 CE until very recently coinciding with migratory movements associated with the civil war and Ebola epidemics. A higher multidrug-resistant (MDR)-TB frequency (23.5 %) than current estimates was obtained together with non-catalogued drug-resistance mutations. Additionally, 39 % of strains were in genomic clusters revealing that ongoing transmission is a major contribution to the TB burden in the country. Our report emphasizes the importance of TB surveillance and control in African countries where bacillus diversity, MDR-TB prevalence and transmission are coalescing to jeopardize TB control programmes.

Translocation of Mycobacterium tuberculosis after experimental ingestion

Human tuberculosis is a life-threatening infection following the inhalation of Mycobacterium tuberculosis, while the closely related bacteria Mycobacterium bovis and Mycobacterium canettii are thought to be transmitted by ingestion. To explore whether M. tuberculosis could also infect individuals by ingestion, male BALBc mice were fed 2 x 10⁶ CFUs of M. tuberculosis Beijing or phosphate-buffered saline as a negative control, over a 28-day experiment. While eight negative control mice remained disease-free, M. tuberculosis was identified in the lymph nodes and lungs of 8/14 mice and in the spleens of 4/14 mice by microscopy, PCR-based detection and culture. Whole-genome sequencing confirmed the identity of the inoculum and the tissue isolates. In these genetically identical mice, the dissemination of M. tuberculosis correlated with the results of the culture detection of four intestinal bacteria. These observations indicate that ingested M. tuberculosis mycobacteria can translocate, notably provoking lymphatic tuberculosis.

Spoligotyping and drug sensitivity of Mycobacterium tuberculosis isolated from pulmonary tuberculosis patients in the Arsi Zone of southeastern Ethiopia

Tuberculosis (TB) is one of the leading causes of morbidity and mortality in different zones of Ethiopia. This study was undertaken to identify the strains of Mycobacterium tuberculosis and evaluate their drug sensitivity profiles in the Arsi Zone. A total of 111 isolates of M. tuberculosis from individuals with pulmonary TB were included and speciation and strain identification were performed using Region of difference 9 and spoligotyping, respectively. The drug sensitivity patterns were assessed using Bactec MGIT 960 SIRE and GenoType MTBDRplus line probe assays. Of 111 isolates, 83% were interpretable and 56 different spoligotype patterns were identified. From these, 22 patterns were shared types while the remaining 34 were orphans. The predominant shared types were spoligotype international type (SIT) 149 and SIT53, comprising 12 and 11 isolates, respectively. Euro-American lineage was the dominant lineage followed by East-African-Indian. Phenotypically, 17.2% of tested isolates were resistant to any first-line drugs and 3.1% were multidrug-resistant. Higher (6.2%) mono-resistance was observed to streptomycin, and no resistance was observed to rifampicin or ethambutol. Genotypically, five (5.4%) isolates were resistant to isoniazid and mutated at codon S315T1 of katG. In contrast, only 1.1% of the isolates were resistant to rifampicin and were mutated at codon S531L of rpoB gene. In this study, a high proportion of orphan strains were isolated, which could suggest the presence of new strains and a high percentage of mono-resistance, warranting the need to strengthen control efforts.

Development and application of affordable SNP typing approaches to genotype Mycobacterium tuberculosis complex strains in low and high burden countries

The Mycobacterium tuberculosis complex (MTBC) comprises the species that causes tuberculosis (TB) which affects 10 million people every year. A robust classification of species, lineages, and sub-lineages is important to explore associations with drug resistance, epidemiological patterns or clinical outcomes. We present a rapid and easy-to-follow methodology to classify clinical TB samples into the main MTBC clades. Approaches are based on the identification of lineage and sub-lineage diagnostic SNP using a real-time PCR high resolution melting assay and classic Sanger sequencing from low-concentrated, low quality DNA. Thus, suitable for implementation in middle and low-income countries. Once we validated our molecular procedures, we characterized a total of 491 biological samples from human and cattle hosts, representing countries with different TB burden. Overall, we managed to genotype ~95% of all samples despite coming from unpurified and low-concentrated DNA. Our approach also allowed us to detect zoonotic cases in eight human samples from Nigeria. To conclude, the molecular techniques we have developed, are accurate, discriminative and reproducible. Furthermore, it costs less than other classic typing methods, resulting in an affordable alternative method in TB laboratories.

Molecular epidemiology of M. tuberculosis in Ethiopia: A systematic review and meta-analysis

The molecular epidemiology of Mycobacterium tuberculosis (M. tuberculosis, Mtb) is poorly documented in Ethiopia. The data that exists has not yet been collected in an overview metadata form. Thus, this review summarizes available literature on the genomic diversity, geospatial distribution and transmission patterns of Mtb lineages (L) and sublineages in Ethiopia. Spoligotyping and Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats (MIRU-VNTR) based articles were identified from MEDLINE via PubMed and Scopus. The last date of article search was done on 12th February 2019. Articles were selected following the PRISMA flow diagram. The proportion of (sub)lineages was summarized at national level and further disaggregated by region. Clustering and recent transmission index (RTI) were determined using metan command and random effect meta-analysis model. The meta-analysis was computed using Stata 14 (Stata Corp. College Station, TX, USA). Among 4371 clinical isolates, 99.5% were Mtb and 0.5% were M. bovis. Proportionally, L4, L3, L1 and L7 made up 62.3%, 21.7%, 7.9% and 3.4% of the total isolates, respectively. Among sublineages, L4.2. ETH/SIT149, L4.10/SIT53, L3. ETH1/SIT25 and L4.6/SIT37 were the leading clustered isolates accounting for 14.4%, 9.7%, 7.2% and 5.5%, respectively. Based on MIRU-VNTR, the rate of clustering was 41% and the secondary case rate from a single source case was estimated at 29%. Clustering and recent transmission index was higher in eastern and southwestern Ethiopia compared with the northwestern part of the country. High level of genetic diversity with a high rate of clustering was noted which collectively mirrored the phenomena of micro-epidemics and super-spreading. The largest set of clustered strains deserves special attention and further characterization using whole genome sequencing (WGS) to better understand the evolution, genomic diversity and transmission dynamics of Mtb.

Nicotinamide and Demographic and Disease transitions: Moderation is Best

Good health and rapid progress depend on an optimal dose of nicotinamide. Too little meat triggers the neurodegenerative condition pellagra and tolerance of symbionts such as tuberculosis (TB), risking dysbioses and impaired resistance to acute infections. Nicotinamide deficiency is an overlooked diagnosis in poor cereal-dependant economies masquerading as 'environmental enteropathy' or physical and cognitive stunting. Too much meat (and supplements) may precipitate immune intolerance and autoimmune and allergic disease, with relative infertility and longevity, via the tryptophan-nicotinamide pathway. This switch favours a dearth of regulatory T (Treg) and an excess of T helper cells. High nicotinamide intake is implicated in cancer and Parkinson's disease. Pro-fertility genes, evolved to counteract high-nicotinamide-induced infertility, may now be risk factors for degenerative disease. Moderation of the dose of nicotinamide could prevent some common diseases and personalised doses at times of stress or, depending on genetic background or age, may treat some other conditions.

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.

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.

Genomic determinants of speciation and spread of the Mycobacterium tuberculosis complex

Models on how bacterial lineages differentiate increase our understanding of early bacterial speciation events and the genetic loci involved. Here, we analyze the population genomics events leading to the emergence of the tuberculosis pathogen. The emergence is characterized by a combination of recombination events involving core pathogenesis functions and purifying selection on early diverging loci. We identify the phoR gene, the sensor kinase of a two-component system involved in virulence, as a key functional player subject to pervasive positive selection after the divergence of the Mycobacterium tuberculosis complex from its ancestor. Previous evidence showed that phoR mutations played a central role in the adaptation of the pathogen to different host species. Now, we show that phoR mutations have been under selection during the early spread of human tuberculosis, during later expansions, and in ongoing transmission events. Our results show that linking pathogen evolution across evolutionary and epidemiological time scales points to past and present virulence determinants.

Insights into the genetic diversity of Mycobacterium tuberculosis in Tanzania

BACKGROUND

Human tuberculosis (TB) is caused by seven phylogenetic lineages of the Mycobacterium tuberculosis complex (MTBC), Lineage 1-7. Recent advances in rapid genotyping of MTBC based on single nucleotide polymorphisms (SNP), allow for phylogenetically robust strain classification, paving the way for defining genotype-phenotype relationships in clinical settings. Such studies have revealed that, in addition to host and environmental factors, strain variation in the MTBC influences the outcome of TB infection and disease. In Tanzania, such molecular epidemiological studies of TB however are scarce in spite of a high TB burden.

METHODS AND FINDINGS

Here we used SNP-typing to characterize a nationwide collection of 2,039 MTBC clinical isolates representative of 1.6% of all new and retreatment TB cases notified in Tanzania during 2012 and 2013. Four lineages, namely Lineage 1-4 were identified within the study population. The distribution and frequency of these lineages varied across regions but overall, Lineage 4 was the most frequent (n = 866, 42.5%), followed by Lineage 3 (n = 681, 33.4%) and 1 (n = 336, 16.5%), with Lineage 2 being the least frequent (n = 92, 4.5%). We found Lineage 2 to be independently associated with female sex (adjusted odds ratio [aOR] 2.14; 95% confidence interval [95% CI] 1.31 - 3.50, p = 0.002) and retreatment cases (aOR 1.67; 95% CI 0.95 - 2.84, p = 0. 065) in the study population. We found no associations between MTBC lineage and patient age or HIV status. Our sublineage typing based on spacer oligotyping on a subset of Lineage 1, 3 and 4 strains revealed the presence of mainly EAI, CAS and LAM families. Finally, we detected low levels of multidrug resistant isolates among a subset of 144 retreatment cases.

CONCLUSIONS

This study provides novel insights into the MTBC lineages and the possible influence of pathogen-related factors on the TB epidemic in Tanzania.

A novel mutation alters the stability of PapA2 resulting in the complete abrogation of sulfolipids in clinical mycobacterial strains

The analysis of whole genomes has revealed specific geographical distribution of Mycobacterium tuberculosis (Mtb) strains across the globe suggestive of unique niche dependent adaptive mechanisms. We provide an important correlation of a genome-based mutation to a molecular phenotype across two predominant clinical Mtb lineages of the Indian subcontinent. We have identified a distinct lineage specific mutation-G247C, translating into an alanine-proline conversion in the papA2 gene of Indo-oceanic lineage 1 (L1) Mtb strains, and restoration of cell wall sulfolipids by simple genetic complementation of papA2 from lineage 3 (L3) or from H37Rv (lineage 4-L4) attributed the loss of this glycolipid to this specific mutation in Indo-Oceanic L1 Mtb. The investigation of structure of Mtb PapA2 revealed a distinct nonribosomal peptide synthetase (NRPS) C domain conformation with an unconventional presence of a zinc binding motif. Surprisingly, the A83P mutation did not map to either the catalytic center in the N-terminal subdomain or any of the substrate-binding region of the protein. On the contrary, the inherent ability of mutant PapA2 to form insoluble aggregates and molecular simulations with the wild-type/mutant (Wt/mut) PapA2 purports an important role for the surface associated 83rd residue in protein conformation. This study demonstrates the importance of a critical structural residue in the papA2 protein of Mtb and helps establish a link between observed genomic alteration and its molecular consequence in the successful human pathogen Mtb. Significance We demonstrate the effect of a unique SNP in PapA2 gene of Indo-oceanic Mycobacterium tuberculosis (Mtb) strains leading to the loss of sulfolipid from these strains. By X-ray crystallographic analysis and molecular dynamics (MD) simulations, we show the importance of this residue in the global PapA2 structure. The presence of a Zn atom has not been reported before for this class of proteins. Here, we provide an important link between genomic alteration and its molecular consequence in Mtb highlighting one of the many adaptive mechanisms that have contributed to its success as a human pathogen. A high degree of identity with PapA1, 3, or 4 would help in interpreting the structure of these PapA proteins and other acyl transferases of other biological systems.

Transmission dynamics study of tuberculosis isolates with whole genome sequencing in southern Sweden

Epidemiological contact tracing complemented with genotyping of clinical Mycobacterium tuberculosis isolates is important for understanding disease transmission. In Sweden, tuberculosis (TB) is mostly reported in migrant and homeless where epidemiologic contact tracing could pose a problem. This study compared epidemiologic linking with genotyping in a low burden country. Mycobacterium tuberculosis isolates (n = 93) collected at Scania University Hospital in Southern Sweden were analysed with the standard genotyping method mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR) and the results were compared with whole genome sequencing (WGS). Using a maximum of twelve single nucleotide polymorphisms (SNPs) as the upper threshold of genomic relatedness noted among hosts, we identified 18 clusters with WGS comprising 52 patients with overall pairwise genetic maximum distances ranging from zero to nine SNPs. MIRU-VNTR and WGS clustered the same isolates, although the distribution differed depending on MIRU-VNTR limitations. Both genotyping techniques identified clusters where epidemiologic linking was insufficient, although WGS had higher correlation with epidemiologic data. To summarize, WGS provided better resolution of transmission than MIRU-VNTR in a setting with low TB incidence. WGS predicted epidemiologic links better which could consolidate and correct the epidemiologically linked cases, avoiding thus false clustering.

Line-probe assay and molecular typing reveal a potential drug resistant clone of Mycobacterium tuberculosis in Ethiopia

Background

Antimicrobial resistance is a global concern of increasing significance. Multidrug resistant tuberculosis (MDR-TB) is spreading worldwide. It is important to monitor trends of antimycobacterial resistance. This is particularly true for high TB burden countries such as Ethiopia where disproportionally less drug sensitivity data are reported from.

Methods

The prevalence of drug resistance was assessed with the line probe assay GenoType MTBDRplus in a set of 161 M. tuberculosis strains that were selected from four common lineages and sub-lineages previously identified in Ethiopia. Most of the tested M. tuberculosis isolates had been genotyped by established Spoligotyping and MIRU-VNTR typing methods.

Results

The proportion of MDR-TB among the isolates was 3.1%. Mono-resistance was 1.2% to rifampicin and 4.3% to isoniazid, and resistance to either of the two first line drugs was 8.7%. Strains of Lineage 4 had the highest resistance rate (13.6%) followed by Lineage 3 (4.9%). None of the isolates representing Lineages 1 and Lineage 7 were drug resistant. Multidrug resistance among pulmonary TB and TB lymphadenitis clinical isolates was 2.8 and 3.7%, respectively. Drug resistance of strains carrying the most prevalent spoligotype in Ethiopia - SIT149 - was further explored. Stratification by MIRU-VNTR identified one genotype with a high rate of drug resistance against Rifampicin and Isoniazid and circulation of a potential MDR-TB clone is proposed.

Conclusion

Although the strain selection was not fully randomized, the overall M. tuberculosis drug resistance rate in this strain set was 8.7% while the rate of MDR was 3.1%. In parallel, we identified a sub-lineage that showed a high rate of resistance to both rifampicin and isoniazid. These resistant strains may belong to a clone of M. tuberculosis that is circulating in the highlands of Ethiopia.

Electronic supplementary material

The online version of this article (10.1186/s40794-018-0075-3) contains supplementary material, which is available to authorized users.

China's tuberculosis epidemic stems from historical expansion of four strains of Mycobacterium tuberculosis

A small number of high-burden countries account for the majority of tuberculosis cases worldwide. Detailed data are lacking from these regions. To explore the evolutionary history of Mycobacterium tuberculosis in China-the country with the third highest tuberculosis burden-we analysed a countrywide collection of 4,578 isolates. Little genetic diversity was detected, with 99.4% of the bacterial population belonging to lineage 2 and three sublineages of lineage 4. The deeply rooted phylogenetic positions and geographic restriction of these four genotypes indicate that their populations expanded in situ following a small number of introductions to China. Coalescent analyses suggest that these bacterial subpopulations emerged in China around 1,000 years ago, and expanded in parallel from the twelfth century onwards, and that the whole population peaked in the late eighteenth century. More recently, sublineage L2.3, which is indigenous to China and exhibited relatively high transmissibility and extensive global dissemination, came to dominate the population dynamics of M. tuberculosis in China. Our results indicate that historical expansion of four M. tuberculosis strains shaped the current tuberculosis epidemic in China, and highlight the long-term genetic continuity of the indigenous M. tuberculosis population.