Large Sequence Polymorphisms of the Euro-American lineage of Mycobacterium tuberculosis: A phylogenetic reconstruction and evidence for convergent evolution in the DR locus

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Association between genotype and drug resistance profiles of Mycobacterium tuberculosis strains circulating in China in a national drug resistance survey

We describe the population structure of a representative collection of 3,133 Mycobacterium tuberculosis isolates, collected within the framework of a national resistance survey from 2007 in China. Genotyping data indicate that the epidemic strains in China can be divided into seven major complexes, of which 92% belonged to the East Asian (mainly Beijing strains) or the Euro-American lineage. The epidemic Beijing strains in China are closely related to the Beijing B0/W148 strain earlier described in Russia and a large cluster of these strains has spread national wide. The density of Beijing strains is high in the whole of China (average 70%), but the highest prevalence was found North of the Yellow river. The Euro-American lineage consists of three sublineages (sublineage_1, 2, and 3) and is more prevalent in the South. Beijing lineage showed the highest cluster rate of 48% and a significantly higher level of resistance to rifampicin (14%, p<0.001), ethambutol (9%, p = 0.001), and ofloxacin (5%, p = 0.011). Within the Euro-American Lineage, sublineage_3 revealed the highest cluster rate (28%) and presented a significantly elevated level of resistance to streptomycin (44%, p<0.001). Our findings suggest that standardised treatment in this region may have contributed to the successful spread of certain strains: sublineage_3 in the Euro-American lineage may have thrived when streptomycin was used without rifampicin for treatment, while later under DOTS based treatment, in which rifampicin plays a key role, Beijing lineage appears to be spreading.

Genetic Diversity of Mycobacterium tuberculosis Isolates from Assam, India: Dominance of Beijing Family and Discovery of Two New Clades Related to CAS1_Delhi and EAI Family Based on Spoligotyping and MIRU-VNTR Typing

Tuberculosis (TB) is one of the major public health concerns in Assam, a remote state located in the northeastern (NE) region of India. The present study was undertaken to explore the circulating genotypes of Mycobacterium tuberculosis complex (MTBC) in this region. A total of 189 MTBC strains were collected from smear positive pulmonary tuberculosis cases from different designated microscopy centres (DMC) from various localities of Assam. All MTBC isolates were cultured on Lowenstein-Jensen (LJ) media and subsequently genotyped using spoligotyping and 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing. Spoligotyping of MTBC isolates revealed 89 distinct spoligo patterns. The most dominant MTBC strain belonged to Beijing lineage and was represented by 35.45% (n = 67) of total isolates, followed by MTBC strains belonging to Central Asian-Delhi (CAS/Delhi) lineage and East African Indian (EAI5) lineage. In addition, in the present study 43 unknown spoligo patterns were detected. The discriminatory power of spoligotyping was found to be 0.8637 based on Hunter Gaston Discriminatory Index (HGDI). On the other hand, 24-loci MIRU-VNTR typing revealed that out of total 189 MTBC isolates from Assam 185 (97.9%) isolates had unique MIRU-VNTR profiles and 4 isolates grouped into 2 clusters. Phylogenetic analysis of 67 Beijing isolates based on 24-loci MIRU-VNTR typing revealed that Beijing isolates from Assam represent two major groups, each comprising of several subgroups. Neighbour-Joining (NJ) phylogenetic tree analysis based on combined spoligotyping and 24-loci MIRU-VNTR data of 78 Non-Beijing isolates was carried out for strain lineage identification as implemented by MIRU-VNTRplus database. The important lineages of MTBC identified were CAS/CAS1_Delhi (41.02%, n = 78) and East-African-Indian (EAI, 33.33%). Interestingly, phylogenetic analysis of orphan (23.28%) MTBC spoligotypes revealed that majority of these orphan isolates from Assam represent two new sub-clades Assam/EAI and Assam/CAS. The prevalence of multidrug resistance (MDR) in Beijing and Non-Beijing strains was found to be 10.44% and 9.01% respectively. In conclusion, the present study has shown the predominance of Beijing isolates in Assam which is a matter of great concern because Beijing strains are considered to be ecologically more fit enabling wider dissemination of M. tuberculosis. Other interesting finding of the present study is the discovery of two new clades of MTBC isolates circulating in Assam. More elaborate longitudinal studies are required to be undertaken in this region to understand the transmission dynamics of MTBC.

Mycobacterium tuberculosis Type II Toxin-Antitoxin Systems: Genetic Polymorphisms and Functional Properties and the Possibility of Their Use for Genotyping

Various genetic markers such as IS-elements, DR-elements, variable number tandem repeats (VNTR), single nucleotide polymorphisms (SNPs) in housekeeping genes and other groups of genes are being used for genotyping. We propose a different approach. We suggest the type II toxin-antitoxin (TA) systems, which play a significant role in the formation of pathogenicity, tolerance and persistence phenotypes, and thus in the survival of Mycobacterium tuberculosis in the host organism at various developmental stages (colonization, infection of macrophages, etc.), as the marker genes. Most genes of TA systems function together, forming a single network: an antitoxin from one pair may interact with toxins from other pairs and even from other families. In this work a bioinformatics analysis of genes of the type II TA systems from 173 sequenced genomes of M. tuberculosis was performed. A number of genes of type II TA systems were found to carry SNPs that correlate with specific genotypes. We propose a minimally sufficient set of genes of TA systems for separation of M. tuberculosis strains at nine basic genotype and for further division into subtypes. Using this set of genes, we genotyped a collection consisting of 62 clinical isolates of M. tuberculosis. The possibility of using our set of genes for genotyping using PCR is also demonstrated.

First insight into the genetic population structure of Mycobacterium tuberculosis isolated from pulmonary tuberculosis patients in Egypt

The present study aimed to assess the population structure of Mycobacterium tuberculosis (MTB) isolates from Egypt. A total of 230 MTB isolates were analysed using spoligotyping, large sequence polymorphism (LSPs), mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing and multi-locus sequence typing (MLST). The majority of isolates (93.0%) belonged to lineage 4, including 44.3, 13.4 and 10.8% of the ill-defined T clade, LAM and Haarlem families, respectively, and lineage 3 was identified in 7.0% of the isolates. MIRU-VNTRs typing allowed efficient discrimination of the spoligotype-defined clusters, including spoligo-international types (SIT) 53, 34, and 4, into 56 patterns, including 13 clusters and 43 unique patterns. A new SNP at position 311614 was identified in all six isolates to form the biggest MIRU-VNTR cluster, which suggested a recent clonal expansion. This SNP could possibly be used as a genetic marker for robust discriminations of Egyptian MTB isolates belonging to SIT53. The combination of spoligotyping, 12 MIRU-VNTRs loci and MLST provided insight into the genetic diversity and transmission dynamics of the Egyptian MTB genotypes and could be a key to implementation of effective control measures by public health authorities.

Predominance of Uganda genotype of Mycobacterium tuberculosis isolated from Ugandan patients with tuberculous lymphadenitis

Background

In Uganda, the emerging Uganda genotype of Mycobacterium tuberculosis is the most common cause of pulmonary tuberculosis (PTB), and accounts for up to 70 % of isolates. Extrapulmonary TB (EPTB) is less studied in Uganda.

Methods

Molecular characterization using deletion analysis and spoligotyping was performed on 121 M. tuberculosis isolates from lymph node fine needle biopsy aspirates of Ugandan patients with tuberculous lymphadenitis. The evolutionary relationships and worldwide distribution of the spoligotypes were analyzed.

Results

Mycobacterium tuberculosis was the only cause of EPTB in this study. The T2 sublineage was the most predominant lineage and the Uganda genotype was the dominant genotype. There were 54 spoligotype patterns among the 121 study isolates. The dominant spoligotypes were shared international types (SIT) SIT420, SIT53, SIT 135, SIT 128 and SIT590 in descending order. All but SIT420 were previously reported in pulmonary TB in this setting. The phylogenetic analysis showed a long descendant branch of spoligotypes belonging to the T2-Uganda sublineage containing specifically SITs 135, 128 and 420.

Conclusion

In most cases, the spoligotypes were similar to those causing PTB, but the Uganda genotype was found to be less common in EPTB than previously reported for PTB in Uganda. The phylogenetic analysis and the study of the worldwide distribution of clustered spoligotypes indicate an ongoing evolution of the Uganda genotype, with the country of Uganda at the center of this evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1362-y) contains supplementary material, which is available to authorized users.

Development and application of single-tube multiplex real-time PCR for lineage classification of Mycobacterium tuberculosis based on large sequence polymorphism in Northeast Thailand

An appreciation of the genetic diversity of Mycobacterium tuberculosis (Mtb) is needed for effective planning of strategies in tuberculosis (TB) control. Large sequence polymorphisms (LSPs) are the molecular epidemiological and evolutionary markers for classification of Mtb into East Asian (EA) or Beijing, Indo-Oceanic (IO), Euro-American (EuA) and East African-Indian (EAI) lineages. We aimed to develop a single-tube multiplex real-time PCR assay using melting curve analysis for lineage classification of Mtb based on LSPs. The technique was optimized and tested with well-characterized strains (n = 89). The developed technique was then applied to classify Mtb isolates from TB patients (n = 256) randomly recruited from 19 provinces covering Northeast Thailand in 2013-2014. The technique demonstrated 100% sensitivity and specificity based on well-characterized strains compared to conventional techniques. The detection limit of the technique is 0.05 ng of genomic DNA of Mtb. The 256 Mtb isolates represented IO (n = 178, 70%), Beijing (n = 60, 23%) and EuA (n = 18, 7%) lineages. Significant associations of the Beijing lineage with drug resistance (p < 0.001) and younger average age of TB patients (p < 0.001) compared to other lineages were shown. The single-tube multiplex real-time PCR technique provides a simple, rapid and high performance tool for characterizing Mtb based on LSPs.

Strain classification of Mycobacterium tuberculosis isolates in Brazil based on genotypes obtained by spoligotyping, mycobacterial interspersed repetitive unit typing and the presence of large sequence and single nucleotide polymorphism

Rio de Janeiro is endemic for tuberculosis (TB) and presents the second largest prevalence of the disease in Brazil. Here, we present the bacterial population structure of 218 isolates of Mycobacterium tuberculosis, derived from 186 patients that were diagnosed between January 2008 and December 2009. Genotypes were generated by means of spoligotyping, 24 MIRU-VNTR typing and presence of fbpC¹⁰³, RD^Rio and RD174. The results confirmed earlier data that predominant genotypes in Rio de Janeiro are those of the Euro American Lineages (99%). However, we observed differences between the classification by spoligotyping when comparing to that of 24 MIRU-VNTR typing, being respectively 43.6% vs. 62.4% of LAM, 34.9% vs. 9.6% of T and 18.3% vs. 21.5% of Haarlem. Among isolates classified as LAM by MIRU typing, 28.0% did not present the characteristic spoligotype profile with absence of spacers 21 to 24 and 32 to 36 and we designated these conveniently as “LAM-like”, 79.3% of these presenting the LAM-specific SNP fbpC¹⁰³. The frequency of RD^Rio and RD174 in the LAM strains, as defined both by spoligotyping and 24 MIRU-VNTR loci, were respectively 11% and 15.4%, demonstrating that RD174 is not always a marker for LAM/RD^Rio strains. We conclude that, although spoligotyping alone is a tool for classification of strains of the Euro-American lineage, when combined with MIRU-VNTRs, SNPs and RD typing, it leads to a much better understanding of the bacterial population structure and phylogenetic relationships among strains of M. tuberculosis in regions with high incidence of TB.

Genomic variability of Mycobacterium tuberculosis strains of the Euro-American lineage based on large sequence deletions and 15-locus MIRU-VNTR polymorphism

A sample of 260 Mycobacterium tuberculosis strains assigned to the Euro-American family was studied to identify phylogenetically informative genomic regions of difference (RD). Mutually exclusive deletions of regions RD115, RD122, RD174, RD182, RD183, RD193, RD219, RD726 and RD761 were found in 202 strains; the RD^Rio deletion was detected exclusively among the RD174-deleted strains. Although certain deletions were found more frequently in certain spoligotype families (i.e., deletion RD115 in T and LAM, RD174 in LAM, RD182 in Haarlem, RD219 in T and RD726 in the “Cameroon” family), the RD-defined sublineages did not specifically match with spoligotype-defined families, thus arguing against the use of spoligotyping for establishing exact phylogenetic relationships between strains. Notably, when tested for katG463/gyrA95 polymorphism, all the RD-defined sublineages belonged to Principal Genotypic Group (PGG) 2, except sublineage RD219 exclusively belonging to PGG3; the 58 Euro-American strains with no deletion were of either PGG2 or 3. A representative sample of 197 isolates was then analyzed by standard 15-locus MIRU-VNTR typing, a suitable approach to independently assess genetic relationships among the strains. Analysis of the MIRU-VNTR typing results by using a minimum spanning tree (MST) and a classical dendrogram showed groupings that were largely concordant with those obtained by RD-based analysis. Isolates of a given RD profile show, in addition to closely related MIRU-VNTR profiles, related spoligotype profiles that can serve as a basis for better spoligotype-based classification.