High-resolution MIRU-VNTRs typing reveals the unique nature of Mycobacterium tuberculosis Beijing genotype in Okinawa, Japan

Distribution and Clonality of drug-resistant tuberculosis in South Africa

Background

Studies have shown that drug-resistant tuberculosis (DR-TB) in South Africa (SA) is clonal and is caused mostly by transmission. Identifying transmission chains is important in controlling DR-TB. This study reports on the sentinel molecular surveillance data of Rifampicin-Resistant (RR) TB in SA, aiming to describe the RR-TB strain population and the estimated transmission of RR-TB cases.

Method

RR-TB isolates collected between 2014 and 2018 from eight provinces were genotyped using combination of spoligotyping and 24-loci mycobacterial interspersed repetitive-units-variable-number tandem repeats (MIRU-VNTR) typing.

Results

Of the 3007 isolates genotyped, 301 clusters were identified. Cluster size ranged between 2 and 270 cases. Most of the clusters (247/301; 82.0%) were small in size (< 5 cases), 12.0% (37/301) were medium sized (5–10 cases), 3.3% (10/301) were large (11–25 cases) and 2.3% (7/301) were very large with 26–270 cases. The Beijing genotype was responsible for majority of RR-TB cases in Western and Eastern Cape, while the East-African-Indian-Somalian (EAI1_SOM) genotype accounted for a third of RR-TB cases in Mpumalanga. The overall proportion of RR-TB cases estimated to be due to transmission was 42%, with the highest transmission-rate in Western Cape (64%) and the lowest in Northern Cape (9%).

Conclusion

Large clusters contribute to the burden of RR-TB in specific geographic areas such as Western Cape, Eastern Cape and Mpumalanga, highlighting the need for community-wide interventions. Most of the clusters identified in the study were small, suggesting close contact transmission events, emphasizing the importance of contact investigations and infection control as the primary interventions in SA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02232-z.

A landscape of genomic alterations at the root of a near-untreatable tuberculosis epidemic

BACKGROUND

Atypical Beijing genotype Mycobacterium tuberculosis strains are widespread in South Africa and have acquired resistance to up to 13 drugs on multiple occasions. It is puzzling that these strains have retained fitness and transmissibility despite the potential fitness cost associated with drug resistance mutations.

METHODS

We conducted Illumina sequencing of 211 Beijing genotype M. tuberculosis isolates to facilitate the detection of genomic features that may promote acquisition of drug resistance and restore fitness in highly resistant atypical Beijing forms. Phylogenetic and comparative genomic analysis was done to determine changes that are unique to the resistant strains that also transmit well. Minimum inhibitory concentration (MIC) determination for streptomycin and bedaquiline was done for a limited number of isolates to demonstrate a difference in MIC between isolates with and without certain variants.

RESULTS

Phylogenetic analysis confirmed that two clades of atypical Beijing strains have independently developed resistance to virtually all the potent drugs included in standard (pre-bedaquiline) drug-resistant TB treatment regimens. We show that undetected drug resistance in a progenitor strain was likely instrumental in this resistance acquisition. In this cohort, ethionamide (ethA A381P) resistance would be missed in first-line drug-susceptible isolates, and streptomycin (gidB L79S) resistance may be missed due to an MIC close to the critical concentration. Subsequent inadequate treatment historically led to amplification of resistance and facilitated spread of the strains. Bedaquiline resistance was found in a small number of isolates, despite lack of exposure to the drug. The highly resistant clades also carry inhA promoter mutations, which arose after ethA and katG mutations. In these isolates, inhA promoter mutations do not alter drug resistance, suggesting a possible alternative role.

CONCLUSION

The presence of the ethA mutation in otherwise susceptible isolates from ethionamide-naïve patients demonstrates that known exposure is not an adequate indicator of drug susceptibility. Similarly, it is demonstrated that bedaquiline resistance can occur without exposure to the drug. Inappropriate treatment regimens, due to missed resistance, leads to amplification of resistance, and transmission. We put these results into the context of current WHO treatment regimens, underscoring the risks of treatment without knowledge of the full drug resistance profile.

Evaluation of the 15 and 24-loci MIRU-VNTR genotyping tools with spoligotyping in the identification of Mycobacterium tuberculosis strains and their genetic diversity in molecular epidemiology studies

BACKGROUND

The transmission dynamics of Mycobacterium tuberculosis (Mtb) using various genotyping tools has been studied globally and a particular tool for genotyping Mtb is the mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR). Tuberculosis (TB) remains an important public health problem worldwide and Sri Lanka being a country of tourist destination; because of major development projects undergoing, it has a high proportion of tourists and immigrants from Asia and Europe that are characterized with highest TB incidences and drug-resistant clinical isolates. Hence, in order to address the question of Mtb genetic diversity, we investigated the discriminatory power of both MIRU-VNTR typing of 15 and 24 loci with spoligotyping to differentiate Mtb isolates.

METHOD

Acid-fast bacilli positive sputum samples (n = 150) from first visit patients were collected. Decontamination of sputum and extraction of genomic DNA were carried out using standard techniques. The isolates were characterized by MIRU-VNTR for both the 15 and 24 loci and spoligotyping.

RESULTS

In our study population, MIRU-VNTR 15 and 24 loci did not show a significant difference among the identified M. tuberculosis strains. However, MIRU 24 loci yielded an additional strain LAM, which is of T1 origin. 15 loci strain grouping had more clusters of strains grouped together while 24 loci differentiated the same cluster of strains into distinct strain types.

CONCLUSION

We conclude that the use of 15-locus MIRU-VNTR typing is sufficient for a first-line epidemiological study to genotype M. tuberculosis, but the additional discriminatory power of 24 loci MIRU-VNTR has been able to differentiate samples within highly homologous groups.

Mycobacterium tuberculosis phylogeography in the context of human migration and pathogen's pathobiology: Insights from Beijing and Ural families

Here, I review the population structure and phylogeography of the two contrasting families of Mycobacterium tuberculosis, Beijing and Ural, in the context of strain pathobiology and human history and migration. Proprietary database (12-loci MIRU-VNTR profiles of 3067 Beijing genotype isolates) was subjected to phylogenetic and statistical analysis. The highest rate (90%) and diversity (HGI 0.80-0.95) of the Beijing genotype in North China suggest it to be its area of origin. Under VNTR-based MDS analysis the interpopulation genetic distances correlated with geography over uninterrupted landmasses. In contrast, large water distances together with long time generated remarkable outliers. Weak and less expected affinities of the distant M. tuberculosis populations may reflect hidden epidemiological links due to unknown migration. Association with drug-resistance or increased virulence/transmissibility along with particular human migration flows shape global dissemination of some Beijing clones. The paucity of data on the Ural genotype prevents from high-resolution analysis that was mainly based on the available spoligotyping data. The North/East Pontic area marked with the highest prevalence of the Ural family may have been the area of its origin and primary dispersal in Eurasia. Ural strains are not marked by increased pathogenic capacities, increased transmissibility and association with drug resistance (but most recent reports describe an alarming increase of MDR Ural strains in some parts of eastern Europe and northwestern Russia). Large-scale SNP or WGS population-based studies targeting strains from indigenous populations and, eventually, analysis of ancient DNA will better test these hypotheses. Host genetics factors likely play the most prominent role in differential dissemination of particular M. tuberculosis genotypes.

Evaluation of the discriminatory power of variable number of tandem repeat typing of Mycobacterium bovis isolates from southern Africa

The usefulness of variable number tandem repeat (VNTR) typing based on limited numbers of loci has previously proven inferior compared to IS6110-RFLP typing when applied to the study of the molecular epidemiology of bovine tuberculosis (BTB) in both livestock and wildlife in southern Africa. In this study, the discriminatory power of 29 published VNTR loci in the characterization of 131 Mycobacterium bovis strains isolated predominantly from wildlife and a smaller number from livestock in southern Africa was assessed. Allelic diversities calculated when loci were evaluated on a selected panel of 23 M. bovis isolates with identified varying degrees of genetic relatedness from different geographic origins as well as M. bovis BCG ranged from 0.00 to 0.63. Of the 29 loci tested, 13 were polymorphic (QUB 11a, QUB 11b, QUB 18, ETR-B and -C, Mtub 21, MIRU 16 and 26, ETR-E, QUB 26, MIRU 23, ETR-A, and Mtub 12). In addition, a comparative evaluation of the 13 loci on a panel of 65 isolates previously characterized by IS6110 restriction fragment length polymorphism (RFLP) typing and further evaluation on 41 isolates with no typing history from Kruger National Park (KNP) highlighted that M. bovis from epidemiologically unrelated cases of BTB in different geographic regions can be adequately distinguished. However, there is a need for improvement of the method to fully discriminate between the parental KNP strain and its clones to allow the detection of evolutionary events causing transmission between and within wildlife species.

Genetic diversity of Mycobacterium tuberculosis isolates from Beijing, China assessed by Spoligotyping, LSPs and VNTR profiles

Background

Tuberculosis is one of the most infectious diseases in the world. Molecular typing methods such as spoligotyping, and VNTR (variable number tandem repeats), IS6110 in the NTF region and LSP (large sequence polymorphisms) analysis are generally useful tools for the resolution of various issues related to the classical epidemiology of Mycobacterium tuberculosis (M. tuberculosis).

Methods

To determine the transmission characteristics of M. tuberculosis strains isolated in Beijing, China, and their genetic relationships, especially those among Beijing family strains, 260 M. tuberculosis strains isolated from patients presenting pulmonary tuberculosis were analyzed by spoligotyping, and by examining 22 VNTR loci and the presence/absence of IS6110 in the NTF region, RD105 and RD181.

Results

81% (211 strains) of the isolates studied were Beijing family strains, 174 (82.5%) of which were identified as modern Beijing strains based on the presence of IS6110 upstream of the NTF region. RD181 was intact in 9 of the other 37 (17.5%) ancestral Beijing strains. The percentage of Beijing family strains in this study was consistent with previous reports. There are many differences, however, in allele diversity among VNTR loci between reports on strains from different areas.

Conclusions

The Beijing family is the most prevalent genotype in Beijing city and the predominance of Beijing family strains has not altered in almost twenty years. Differences in the alleles and discrimination ability of VNTR loci between different regions is likely due to population differences in the regions where these M. tuberculosis strains were isolated or to differences in sampling times.

Genetic diversity and transmission characteristics of Beijing family strains of Mycobacterium tuberculosis in Peru

Beijing family strains of Mycobacterium tuberculosis have attracted worldwide attention because of their wide geographical distribution and global emergence. Peru, which has a historical relationship with East Asia, is considered to be a hotspot for Beijing family strains in South America. We aimed to unveil the genetic diversity and transmission characteristics of the Beijing strains in Peru. A total of 200 Beijing family strains were identified from 2140 M. tuberculosis isolates obtained in Lima, Peru, between December 2008 and January 2010. Of them, 198 strains were classified into sublineages, on the basis of 10 sets of single nucleotide polymorphisms (SNPs). They were also subjected to variable number tandem-repeat (VNTR) typing using an international standard set of 15 loci (15-MIRU-VNTR) plus 9 additional loci optimized for Beijing strains. An additional 70 Beijing family strains, isolated between 1999 and 2006 in Lima, were also analyzed in order to make a longitudinal comparison. The Beijing family was the third largest spoligotyping clade in Peru. Its population structure, by SNP typing, was characterized by a high frequency of Sequence Type 10 (ST10), which belongs to a modern subfamily of Beijing strains (178/198, 89.9%). Twelve strains belonged to the ancient subfamily (ST3 [n=3], ST25 [n=1], ST19 [n=8]). Overall, the polymorphic information content for each of the 24 loci values was low. The 24 loci VNTR showed a high clustering rate (80.3%) and a high recent transmission index (RTI(n-1)=0.707). These strongly suggest the active and on-going transmission of Beijing family strains in the survey area. Notably, 1 VNTR genotype was found to account for 43.9% of the strains. Comparisons with data from East Asia suggested the genotype emerged as a uniquely endemic clone in Peru. A longitudinal comparison revealed the genotype was present in Lima by 1999.

The genotypic population structure of Mycobacterium tuberculosis complex from Moroccan patients reveals a predominance of Euro-American lineages

BACKGROUND

Tuberculosis (TB) remains a major health problem in Morocco. Characterization of circulating Mycobacterium tuberculosis genotypic lineages, important to understand the dynamic of the disease, was hereby addressed for the first time at a national level.

METHODOLOGY/PRINCIPAL FINDINGS

Spoligotyping was performed on a panel of 592 M. tuberculosis complex strains covering a 2-year period (2004-2006). It identified 129 patterns: 105 (n = 568 strains) corresponded to a SIT number in the SITVIT2 database, while 24 patterns were labeled as orphan. A total of 523 (88.3%) strains were clustered vs. 69 or 11.7% unclustered. Classification of strains within 3 large phylogenetical groups was as follows: group 1- ancestral/TbD1+/PGG1 (EAI, Bovis, Africanum), group 2- modern/TbD1-/PGG1 group (Beijing, CAS), group 3- evolutionary recent/TbD1-/PGG2/3 (Haarlem, X, S, T, LAM; alternatively designated as the Euro-American lineage). As opposed to group 3 strains (namely LAM, Haarlem, and T) that predominated (86.5% of all isolates), 6 strains belonged to group 2 (Beijing n = 5, CAS n = 1), and 3 strains (BOV_1 n = 2, BOV_4-CAPRAE) belonged to ancestral group 1 (EAI and AFRI lineage strains were absent). 12-loci MIRU-VNTR typing of the Casablanca subgroup (n = 114 strains) identified 71 patterns: 48 MITs and 23 orphan patterns; it allowed to reduce the clustering rate from 72.8% to 29.8% and the recent transmission rate from 64% to 20.2%.

CONCLUSION

The M. tuberculosis population structure in Morocco is highly homogeneous, and is characterized by the predominance of the Euro-American lineages, namely LAM, Haarlem, and T, which belong to the "evolutionary recent" TbD1-/PGG2/3 phylogenetic group. The combination of spoligotyping and MIRUs decreased the clustering rate significantly, and should now be systematically applied in larger studies. The methods used in this study appear well suited to monitor the M. tuberculosis population structure for an enhanced TB management program in Morocco.

Scanning of genetic diversity of evolutionarily sequential Mycobacterium tuberculosis Beijing family strains based on genome wide analysis

The Beijing family is an endemic lineage of Mycobacterium tuberculosis in eastern Asia. In Japan, five evolutionarily sequential sublineages composing the lineage have predominated. Comparative genomic sequencing based on a microarray technique was conducted for five representative strains of those respective sublineages. Results revealed approximately 200 point mutations specific to each strain. Subsequently, to investigate the genetic diversity of each sublineage, we analysed the phylogenetic divergence of 103 domestic strains belonging to them using genetic markers derived from the mutation information. Results show that the five sublineages have comprised smaller lineages which had diverged at various points. The smaller sub-sublineages have emerged with respective bottlenecks, which are reflected in the excessive monophyletic evolution of the species. Our data provide necessary information to grasp a comprehensive picture of genetic diversity of the lineage constructed in its evolution.

Investigation on Mycobacterium tuberculosis diversity in China and the origin of the Beijing clade

Background

Investigation of the genetic diversity of Mycobacterium tuberculosis in China has shown that Beijing genotype strains play a dominant role in the tuberculosis (TB) epidemic. In order to examine the strain diversity in the whole country, and to study the evolutionary development of Beijing strains, we sought to genotype a large collection of isolates using different methods.

Methodology/Principal Findings

We applied a 15-loci VNTR typing analysis on 1,586 isolates from the Beijing municipality and 12 Chinese provinces or autonomous regions. The data was compared to that of 900 isolates from various other worldwide geographic regions outside of China. A total of 1,162/1,586 (73.2%) of the isolates, distributed into 472 VNTR types, were found to belong to the Beijing genotype family and this represented 56 to 94% of the isolates in each of the localizations. VNTR typing revealed that the majority of the non-Beijing isolates fall into two genotype families, which represented 17% of the total number of isolates, and seem largely restricted to China. A small number of East African Indian genotype strains was also observed in this collection. Ancient Beijing strains with an intact region of difference (RD) 181, as well as strains presumably resembling ancestors of the whole Beijing genotype family, were mainly found in the Guangxi autonomous region.

Conclusions/Significance

This is the largest M. tuberculosis VNTR-based genotyping study performed in China to date. The high percentage of Beijing isolates in the whole country and the presence in the South of strains representing early branching points may be an indication that the Beijing lineage originated from China, probably in the Guangxi region. Two modern lineages are shown here to represent the majority of non-Beijing Chinese isolates. The observed geographic distribution of the different lineages within China suggests that natural frontiers are major factors in their diffusion.