Local adaptation in populations of Mycobacterium tuberculosis endemic to the Indian Ocean Rim

Population structure and spatial distribution of Mycobacterium tuberculosis in Ethiopia

Ethiopia is one of the countries with a high tuberculosis (TB) burden, yet little is known about the spatial distribution of Mycobacterium tuberculosis (Mtb) lineages. This study identifies the spoligotyping of 1735 archived Mtb isolates from the National Drug Resistance Survey, collected between November 2011 and June 2013, to investigate Mtb population structure and spatial distribution. Spoligotype International Types (SITs) and lineages were retrieved from online databases. The distribution of lineages was evaluated using Fisher's exact test and logistic regression models. The Global Moran's Index and Getis-Ord Gi statistic were utilized to identify hotspot areas. Our results showed that spoligotypes could be interpreted and led to 4 lineages and 283 spoligotype patterns in 91% of the isolates, including 4% of those with multidrug/rifampicin resistance (MDR/RR) TB. The identified Mtb lineages were lineage 1 (1.8%), lineage 3 (25.9%), lineage 4 (70.6%) and lineage 7 (1.6%). The proportion of lineages 3 and 4 varied by regions, with lineage 3 being significantly greater than lineage 4 in reports from Gambella (AOR = 4.37, P < 0.001) and Tigray (AOR = 3.44, P = 0.001) and lineage 4 being significantly higher in Southern Nations Nationalities and Peoples Region (AOR = 1.97, P = 0.026) than lineage 3. Hotspots for lineage 1 were located in eastern Ethiopia, while a lineage 7 hotspot was identified in northern and western Ethiopia. The five prevalent spoligotypes, which were SIT149, SIT53, SIT25, SIT37 and SIT26 account for 42.8% of all isolates under investigation, while SIT149, SIT53 and SIT21 account for 52-57.8% of drug-resistant TB cases. TB and drug resistant TB are mainly caused by lineages 3 and 4, and significant proportions of the prevalent spoligotypes also influence drug-resistant TB and the total TB burden. Regional variations in lineages may result from both local and cross-border spread.

HIV co-infection is associated with reduced Mycobacterium tuberculosis transmissibility in sub-Saharan Africa

Persons living with HIV are known to be at increased risk of developing tuberculosis (TB) disease upon infection with Mycobacterium tuberculosis (Mtb). However, it has remained unclear how HIV co-infection affects subsequent Mtb transmission from these patients. Here, we customized a Bayesian phylodynamic framework to estimate the effects of HIV co-infection on the Mtb transmission dynamics from sequence data. We applied our model to four Mtb genomic datasets collected in sub-Saharan African countries with a generalized HIV epidemic. Our results confirm that HIV co-infection is a strong risk factor for developing active TB. Additionally, we demonstrate that HIV co-infection is associated with a reduced effective reproductive number for TB. Stratifying the population by CD4+ T-cell count yielded similar results, suggesting that, in this context, CD4+ T-cell count is not a better predictor of Mtb transmissibility than HIV infection status alone. Together, our genome-based analyses complement observational household contact studies, and more firmly establish the negative association between HIV co-infection and Mtb transmissibility.

Bacterial diversity dominates variable macrophage responses of tuberculosis patients in Tanzania

The Mycobacterium tuberculosis complex (MTBC) comprises nine human-adapted lineages that differ in their geographical distribution. Local adaptation of specific MTBC genotypes to the respective human host population has been invoked in this context. We aimed to assess if bacterial genetics governs MTBC pathogenesis or if local co-adaptation translates into differential susceptibility of human macrophages to infection by different MTBC genotypes. We generated macrophages from cryopreserved blood mononuclear cells of Tanzanian tuberculosis patients, from which the infecting MTBC strains had previously been phylogenetically characterized. We infected these macrophages ex vivo with a phylogenetically similar MTBC strain ("matched infection") or with strains representative of other MTBC lineages ("mismatched infection"). We found that L1 infections resulted in a significantly lower bacterial burden and that the intra-cellular replication rate of L2 strains was significantly higher compared the other MTBC lineages, irrespective of the MTBC lineage originally infecting the patients. Moreover, L4-infected macrophages released significantly greater amounts of TNF-α, IL-6, IL-10, MIP-1β, and IL-1β compared to macrophages infected by all other strains. While our results revealed no measurable effect of local adaptation, they further highlight the strong impact of MTBC phylogenetic diversity on the variable outcome of the host-pathogen interaction in human tuberculosis.

Evolution and transmission of antibiotic resistance is driven by Beijing lineage Mycobacterium tuberculosis in Vietnam

IMPORTANCE

Drug-resistant tuberculosis (TB) infection is a growing and potent concern, and combating it will be necessary to achieve the WHO's goal of a 95% reduction in TB deaths by 2035. While prior studies have explored the evolution and spread of drug resistance, we still lack a clear understanding of the fitness costs (if any) imposed by resistance-conferring mutations and the role that Mtb genetic lineage plays in determining the likelihood of resistance evolution. This study offers insight into these questions by assessing the dynamics of resistance evolution in a high-burden Southeast Asian setting with a diverse lineage composition. It demonstrates that there are clear lineage-specific differences in the dynamics of resistance acquisition and transmission and shows that different lineages evolve resistance via characteristic mutational pathways.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Comparative Study on Tuberculosis Drug Resistance and Molecular Detection Methods Among Different Mycobacterium Tuberculosis Lineages

Purpose

This study aims to compare drug resistance and detection efficacy across different Mycobacterium tuberculosis lineages, offering insights for precise treatment and molecular diagnosis.

Methods

161 strains of Mycobacterium tuberculosis (M.tb) were tested for drug resistance using Phenotypic Drug Susceptibility Testing (pDST), High-Resolution Melting analysis (HRM), and Whole Genome Sequencing (WGS) methods. The main focus was on evaluating the accuracy of different methods for detecting resistance to rifampicin (RIF), isoniazid (INH), and streptomycin (SM).

Results

Among the 161 strains of M.tb, 83.85% (135/161) were fully sensitive to RIF, INH, and SM according to pDST, and the rate of multidrug resistance was 4.35% (7/161). The drug resistance rates of lineage 2 M.tb to the three drugs (26/219, 11.87%) were significantly higher than those of non-lineage 2 M.tb (12/264, 4.45%) (P<0.05). Compared with pDST, WGS had a sensitivity of 100%, 94.12%, and 92.31% and a specificity of 100%, 99.31%, and 98.65% for RIF, INH, and SM, respectively, with no significant difference. The sensitivity of HRM for RIF, INH, and SM was 87.50%, 52.94%, and 76.92%, respectively, while the specificity was 96.08%, 99.31%, and 99.32%, respectively. The sensitivity of HRM for detecting INH resistance was significantly lower than that of pDST (P=0.039). Compared with HRM, WGS increased the sensitivity of RIF, INH, and SM by 12.50%, 41.18%, and 15.38%, respectively.

Conclusion

There are significant differences in drug resistance rates among different lineages of M.tb, with lineage 2 having higher rates of RIF, INH, and SM resistance than lineages 3 and 4. The sensitivity of HRM is far lower than that of pDST, and currently, the accuracy of HRM is not sufficient to replace pDST. WGS has no significant difference in detecting drug resistance compared with pDST but can identify new anti-tuberculosis drug-resistant mutations, providing effective guidance for clinical decision-making.

Genomic changes underpinning the emergence of a successful Mycobacterium tuberculosis Latin American and Mediterranean clonal complex

Introduction

The Latin American and Mediterranean sublineage (L4.3/LAM) is the most common generalist sublineage of Mycobacterium tuberculosis lineage 4 (L4), yet certain L4.3/LAM genotypes appear to be confined to particular geographic regions. This is typically the case of a L4.3/LAM clonal complex (CC), TUN4.3_CC1, which is the most preponderant in Tunisia (61.5% of L4.3/LAM).

Methods

Here, we used whole-genome sequencing data of 346 globally distributed L4 clinical strains, including 278 L4.3/LAM isolates, to reconstruct the evolutionary history of TUN4.3_CC1 and delineate critical genomic changes underpinning its success.

Results and Discussion

Phylogenomic coupled to phylogeographic analyses indicated that TUN4.3_CC1 has evolved locally, being confined mainly to North Africa. Maximum likelihood analyses using the site and branch-site models of the PAML package disclosed strong evidence of positive selection in the gene category "cell wall and cell processes" of TUN4.3_CC1. Collectively, the data indicate that TUN4.3_CC1 has inherited several mutations, which could have potentially contributed to its evolutionary success. Of particular interest are amino acid replacements at the esxK and eccC2 genes of the ESX/Type VII secretion system, which were found to be specific to TUN4.3_CC1, being common to almost all isolates. Because of its homoplastic nature, the esxK mutation could potentially have endowed TUN4.3_CC1 with a selective advantage. Moreover, we noticed the occurrence of additional, previously described homoplasic nonsense mutations in ponA1 and Rv0197. The mutation in the latter gene, a putative oxido-reductase, has previously been shown to be correlated with enhanced transmissibility in vivo. In sum, our findings unveiled several features underpinning the success of a locally evolved L4.3/LAM clonal complex, lending further support to the critical role of genes encoded by the ESX/type VII secretion system.

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

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

Fine-grain population structure and transmission patterns of Mycobacterium tuberculosis in southern Mozambique, a high TB/HIV burden area

Genomic studies of the Mycobacterium tuberculosis complex (MTBC) might shed light on the dynamics of its transmission, especially in high-burden settings, where recent outbreaks are embedded in the complex natural history of the disease. To this end, we conducted a 1 year prospective surveillance-based study in Mozambique. We applied whole-genome sequencing (WGS) to 295 positive cultures. We fully characterized MTBC isolates by phylogenetics and dating evaluation, and carried out a molecular epidemiology analysis to investigate further associations with pre-defined transmission risk factors. The majority of strains (49.5%, 136/275) belonged to lineage (L) 4; 57.8 % of them (159/275) were in genomic transmission clusters (cut-off 5 SNPs), and a strikingly high proportion (45.5%) shared an identical genotype (0 SNP pairwise distance). We found two ‘likely endemic’ clades, comprising 67 strains, belonging to L1.2, which dated back to the late 19th century and were associated with recent spread among people living with human immunodeficiency virus (PLHIV). We describe for the first time the population structure of MTBC in our region, a high tuberculosis (TB)/HIV burden area. Clustering analysis revealed an unforeseen pattern of spread and high rates of progression to active TB, suggesting weaknesses in TB control activities. The long-term presence of local strains in Mozambique, which were responsible for large transmission among HIV/TB-coinfected patients, calls into question the role of HIV in TB transmission.

Understanding drivers of phylogenetic clustering and terminal branch lengths distribution in epidemics of Mycobacterium tuberculosis

Detecting factors associated with transmission is important to understand disease epidemics, and to design effective public health measures. Clustering and terminal branch lengths (TBL) analyses are commonly applied to genomic data sets of Mycobacterium tuberculosis (MTB) to identify sub-populations with increased transmission. Here, I used a simulation-based approach to investigate what epidemiological processes influence the results of clustering and TBL analyses, and whether differences in transmission can be detected with these methods. I simulated MTB epidemics with different dynamics (latency, infectious period, transmission rate, basic reproductive number R0, sampling proportion, sampling period, and molecular clock), and found that all considered factors, except for the length of the infectious period, affect the results of clustering and TBL distributions. I show that standard interpretations of this type of analyses ignore two main caveats: (1) clustering results and TBL depend on many factors that have nothing to do with transmission, (2) clustering results and TBL do not tell anything about whether the epidemic is stable, growing, or shrinking, unless all the additional parameters that influence these metrics are known, or assumed identical between sub-populations. An important consequence is that the optimal SNP threshold for clustering depends on the epidemiological conditions, and that sub-populations with different epidemiological characteristics should not be analyzed with the same threshold. Finally, these results suggest that different clustering rates and TBL distributions, that are found consistently between different MTB lineages, are probably due to intrinsic bacterial factors, and do not indicate necessarily differences in transmission or evolutionary success.

Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection

SignificancePrevious attempts to identify the action of natural selection in the Mycobacterium tuberculosis complex (MTBC) were limited by sample size and averaging across time and lineages. We investigate changes in selective pressures across time for every single gene of the MTBC. We developed a methodology to analyze temporal signals of selection in a large dataset (∼5,000 complete genomes) and showed that 1) almost half of the genes seem to have been under positive selection at some point in time; 2) experimentally confirmed epitopes tend to accumulate more mutations in deeper branches than in external branches; and 3) temporal signals identify genes that were conserved in the past but under positive selection in the present, suggesting ongoing adaptation to the host.

Sub-Lineage Specific Phenolic Glycolipid Patterns in the Mycobacterium tuberculosis Complex Lineage 1

“Ancestral” Mycobacterium tuberculosis complex (MTBC) strains of Lineage 1 (L1, East African Indian) are a prominent tuberculosis (TB) cause in countries around the Indian Ocean. However, the pathobiology of L1 strains is insufficiently characterized. Here, we used whole genome sequencing (WGS) of 312 L1 strains from 43 countries to perform a characterization of the global L1 population structure and correlate this to the analysis of the synthesis of phenolic glycolipids (PGL) – known MTBC polyketide-derived virulence factors. Our results reveal the presence of eight major L1 sub-lineages, whose members have specific mutation signatures in PGL biosynthesis genes, e.g., pks15/1 or glycosyltransferases Rv2962c and/or Rv2958c. Sub-lineage specific PGL production was studied by NMR-based lipid profiling and strains with a completely abolished phenolphthiocerol dimycoserosate biosynthesis showed in average a more prominent growth in human macrophages. In conclusion, our results show a diverse population structure of L1 strains that is associated with the presence of specific PGL types. This includes the occurrence of mycoside B in one sub-lineage, representing the first description of a PGL in an M. tuberculosis lineage other than L2. Such differences may be important for the evolution of L1 strains, e.g., allowing adaption to different human populations.

Whole-genome single nucleotide variant phylogenetic analysis of Mycobacterium tuberculosis Lineage 1 in endemic regions of Asia and Africa

Mycobacterium tuberculosis (Mtb) lineage 1 (L1) contributes considerably to the disease morbidity. While whole genome sequencing (WGS) is increasingly used for studying Mtb, our understanding of genetic diversity of L1 remains limited. Using phylogenetic analysis of WGS data from endemic range in Asia and Africa, we provide an improved genotyping scheme for L1. Mapping deletion patterns of the 68 direct variable repeats (DVRs) in the CRISPR region of the genome onto the phylogeny provided supporting evidence that the CRISPR region evolves primarily by deletion, and hinted at a possible Southeast Asian origin of L1. Both phylogeny and DVR patterns clarified some relationships between different spoligotypes, and highlighted the limited resolution of spoligotyping. We identified a diverse repertoire of drug resistance mutations. Altogether, this study demonstrates the usefulness of WGS data for understanding the genetic diversity of L1, with implications for public health surveillance and TB control. It also highlights the need for more WGS studies in high-burden but underexplored regions.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

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.

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.

Tuberculosis in Australia's tropical north: a population-based genomic epidemiological study

BACKGROUND

The Northern Territory (NT) has the highest tuberculosis (TB) rate of all Australian jurisdictions. We combined TB public health surveillance data with genomic sequencing of Mycobacterium tuberculosis isolates in the tropical 'Top End' of the NT to investigate trends in TB incidence and transmission.

METHODS

This retrospective observational study included all 741 culture-confirmed cases of TB in the Top End over three decades from 1989-2020. All 497 available M. tuberculosis isolates were sequenced. We used contact tracing data to define a threshold pairwise SNP distance for hierarchical single linkage clustering, and examined putative transmission clusters in the context of epidemiologic information.

FINDINGS

There were 359 (48%) cases born overseas, 329 (44%) cases among Australian First Nations peoples, and 52 (7%) cases were Australian-born and non-Indigenous. The annual incidence in First Nations peoples from 1989-2019 fell from average 50.4 to 11.0 per 100,000 (P<0·001). First Nations cases were more likely to die from TB (41/329, 12·5%) than overseas-born cases (11/359, 3·1%; P<0·001). Using a threshold of ≤12 SNPs, 28 clusters of between 2-64 individuals were identified, totalling 250 cases; 214 (86%) were First Nations cases and 189 (76%) were from a remote region. The time between cases and past epidemiologically- and genomically-linked contacts ranged from 4·5 months to 24 years.

INTERPRETATION

Our findings support prioritisation of timely case detection, contact tracing augmented by genomic sequencing, and latent TB treatment to break transmission chains in Top End remote hotspot regions.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.