Differences in gene expression between clonal variants of Mycobacterium tuberculosis emerging as a result of microevolution

The mutation rate of rpoB gene showed an upward trend with the increase of MIRU10, MIRU39 and QUB4156 repetitive number

BACKGROUND

Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) is a frequently used typing method for identifying the Beijing genotype of Mycobacterium tuberculosis (Mtb), which is easily transformed into rifampicin (RIF) resistance. The RIF resistance of Mtb is considered to be highly related with the mutation of rpoB gene. Therefore, this study aimed to analyze the relationship between the repetitive number of MIRU loci and the mutation of rpoB gene.

METHODS

An open-source whole-genome sequencing data of Mtb was used to detect the mutation of rpoB gene and the repetitive number of MIRU loci by bioinformatics methods. Cochran-Armitage analysis was performed to analyze the trend of the rpoB gene mutation rate and the repetitive number of MIRU loci.

RESULTS

Among 357 rifampicin-resistant tuberculosis (RR-TB), 304 strains with mutated rpoB genes were detected, and 6 of 67 rifampicin susceptible strains were detected mutations. The rpoB gene mutational rate showed an upward trend with the increase of MIRU10, MIRU39, QUB4156 and MIRU16 repetitive number, but only the repetitive number of MIRU10, MRIU39 and QUB4156 were risk factors for rpoB gene mutation. The Hunter-Gaston discriminatory index (HGDI) of MIRU10 (0.65) and QUB4156 (0.62) was high in the overall sample, while MIRU39 (0.39) and MIRU16 (0.43) showed a moderate discriminatory Power.

CONCLUSION

The mutation rate of rpoB gene increases with the addition of repetitive numbers of MIRU10, QUB4156 and MIRU39 loci.

Rapid adaptation of a complex trait during experimental evolution of Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (M. tb), is a leading cause of death due to infectious disease. TB is not traditionally associated with biofilms, but M. tb biofilms are linked with drug and immune tolerance and there is increasing recognition of their contribution to the recalcitrance of TB infections. Here, we used M. tb experimental evolution to investigate this complex phenotype and identify candidate loci controlling biofilm formation. We identified novel candidate loci, adding to our understanding of the genetic architecture underlying M. tb biofilm development. Under selective pressure to grow as a biofilm, regulatory mutations rapidly swept to fixation and were associated with changes in multiple traits, including extracellular matrix production, cell size, and growth rate. Genetic and phenotypic paths to enhanced biofilm growth varied according to the genetic background of the parent strain, suggesting that epistatic interactions are important in M. tb adaptation to changing environments.

Methods for Detecting Mycobacterial Mixed Strain Infections-A Systematic Review

Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with multiple strains of a single pathogenic species. Known to occur in humans and animals, MSIs deserve special consideration when studying transmission dynamics, evolution, and treatment of mycobacterial diseases, notably tuberculosis in humans and paratuberculosis (or Johne's disease) in ruminants. Therefore, a systematic review was conducted to examine how MSIs are defined in the literature, how widespread the phenomenon is across the host species spectrum, and to document common methods used to detect such infections. Our search strategy identified 121 articles reporting MSIs in both humans and animals, the majority (78.5%) of which involved members of the Mycobacterium tuberculosis complex, while only a few (21.5%) examined non-tuberculous mycobacteria (NTM). In addition, MSIs exist across various host species, but most reports focused on humans due to the extensive amount of work done on tuberculosis. We reviewed the strain typing methods that allowed for MSI detection and found a few that were commonly employed but were associated with specific challenges. Our review notes the need for standardization, as some highly discriminatory methods are not adapted to distinguish between microevolution of one strain and concurrent infection with multiple strains. Further research is also warranted to examine the prevalence of NTM MSIs in both humans and animals. In addition, it is envisioned that the accurate identification and a better understanding of the distribution of MSIs in the future will lead to important information on the epidemiology and pathophysiology of mycobacterial diseases.

Local adaptive evolution of two distinct clades of Beijing and T families of Mycobacterium tuberculosis in Chongqing: a Bayesian population structure and phylogenetic study

Background

Beijing sub-pedigree 2 (BSP2) and T sub-lineage 6 (TSL6) are two clades belonging to Beijing and T family of Mycobacterium tuberculosis (MTB), respectively, defined by Bayesian population structure analysis based on 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTR). Globally, over 99% of BSP2 and 89% of TSL6 isolates were distributed in Chongqing, suggesting their possible local adaptive evolution. The objective of this paper is to explore whether BSP2 and TSL6 originated by their local adaptive evolution from the specific isolates of Beijing and T families in Chongqing.

Methods

The genotyping data of 16 090 MTB isolates were collected from laboratory collection, published literatures and SITVIT database before subjected to Bayesian population structure analysis based on 24-loci MIRU-VNTR. Spacer Oligonucleotide Forest (Spoligoforest) and 24-loci MIRU-VNTR-based minimum spanning tree (MST) were used to explore their phylogenetic pathways, with Bayesian demographic analysis for exploring the recent demographic change of TSL6.

Results

Phylogenetic analysis suggested that BSP2 and TSL6 in Chongqing may evolve from BSP4 and TSL5, respectively, which were locally predominant in Tibet and Jiangsu, respectively. Spoligoforest showed that Beijing and T families were genetically distant, while the convergence of the MIRU-VNTR pattern of BSP2 and TSL6 was revealed by WebLogo. The demographic analysis concluded that the recent demographic change of TSL6 might take 111.25 years.

Conclusions

BSP2 and TSL6 clades might originate from BSP4 and TSL5, respectively, by their local adaptive evolution in Chongqing. Our study suggests MIRU-VNTR be combined with other robust markers for a more comprehensive genotyping approach, especially for families of clades with the same MIRU-VNTR pattern.

Deciphering Within-Host Microevolution of Mycobacterium tuberculosis through Whole-Genome Sequencing: the Phenotypic Impact and Way Forward

The Mycobacterium tuberculosis genome is more heterogenous and less genetically stable within the host than previously thought. Currently, only limited data exist on the within-host microevolution, diversity, and genetic stability of M. tuberculosis As a direct consequence, our ability to infer M. tuberculosis transmission chains and to understand the full complexity of drug resistance profiles in individual patients is limited. Furthermore, apart from the acquisition of certain drug resistance-conferring mutations, our knowledge on the function of genetic variants that emerge within a host and their phenotypic impact remains scarce. We performed a systematic literature review of whole-genome sequencing studies of serial and parallel isolates to summarize the knowledge on genetic diversity and within-host microevolution of M. tuberculosis We identified genomic loci of within-host emerged variants found across multiple studies and determined their functional relevance. We discuss important remaining knowledge gaps and finally make suggestions on the way forward.

Mycobacterium tuberculosis Acquires Limited Genetic Diversity in Prolonged Infections, Reactivations and Transmissions Involving Multiple Hosts

Background:Mycobacterium tuberculosis (MTB) has limited ability to acquire variability. Analysis of its microevolution might help us to evaluate the pathways followed to acquire greater infective success. Whole-genome sequencing (WGS) in the analysis of the transmission of MTB has elucidated the magnitude of variability in MTB. Analysis of transmission currently depends on the identification of clusters, according to the threshold of variability (<5 SNPs) between isolates. Objective: We evaluated whether the acquisition of variability in MTB, was more frequent in situations which could favor it, namely intrapatient, prolonged infections or reactivations and interpatient transmissions involving multiple sequential hosts. Methods: We used WGS to analyze the accumulation of variability in sequential isolates from prolonged infections or translations from latency to reactivation. We then measured microevolution in transmission clusters with prolonged transmission time, high number of involved cases, simultaneous involvement of latency and active transmission. Results: Intrapatient and interpatient acquisition of variability was limited, within the ranges expected according to the thresholds of variability proposed, even though bursts of variability were observed. Conclusions: The thresholds of variability proposed for MTB seem to be valid in most circumstances, including those theoretically favoring acquisition of variability. Our data point to multifactorial modulation of microevolution, although further studies are necessary to elucidate the factors underlying this modulation.

In-Depth Characterization and Functional Analysis of Clonal Variants in a Mycobacterium tuberculosis Strain Prone to Microevolution

The role of clonal complexity has gradually been accepted in infection by Mycobacterium tuberculosis (MTB), although analyses of this issue are limited. We performed an in-depth study of a case of recurrent MTB infection by integrating genotyping, whole genome sequencing, analysis of gene expression and infectivity in in vitro and in vivo models. Four different clonal variants were identified from independent intrapatient evolutionary branches. One of the single-nucleotide polymorphisms in the variants mapped in mce3R, which encodes a repressor of an operon involved in virulence, and affected expression of the operon. Competitive in vivo and in vitro co-infection assays revealed higher infective efficiency for one of the clonal variants. A new clonal variant, which had not been observed in the clinical isolates, emerged in the infection assays and showed higher fitness than its parental strain. The analysis of other patients involved in the same transmission cluster revealed new clonal variants acquired through novel evolutionary routes, indicating a high tendency toward microevolution in some strains that is not host-dependent. Our study highlights the need for integration of various approaches to advance our knowledge of the role and significance of microevolution in tuberculosis.

Biological and Epidemiological Consequences of MTBC Diversity

Tuberculosis is caused by different groups of bacteria belonging to the Mycobacterium tuberculosis complex (MTBC). The combined action of human factors, environmental conditions and bacterial virulence determine the extent and form of human disease. MTBC virulence is a composite of different clinical phenotypes such as transmission rate and disease severity among others. Clinical phenotypes are also influenced by cellular and immunological phenotypes. MTBC phenotypes are determined by the genotype, therefore finding genotypes responsible for clinical phenotypes would allow discovering MTBC virulence factors. Different MTBC strains display different cellular and clinical phenotypes. Strains from Lineage 5 and Lineage 6 are metabolically different, grow slower, and are less virulent. Also, at least certain groups of Lineage 2 and Lineage 4 strains are more virulent in terms of disease severity and human-to-human transmission. Because phenotypic differences are ultimately caused by genotypic differences, different genomic loci have been related to various cellular and clinical phenotypes. However, defining the impact of specific bacterial genomic loci on virulence when other bacterial determinants, human and environmental factors are also impacting the phenotype would contribute to a better knowledge of tuberculosis virulence and ultimately benefit tuberculosis control.

Whole genome sequencing reveals mycobacterial microevolution among concurrent isolates from sputum and blood in HIV infected TB patients

Background

In the context of advanced immunosuppression, M. tuberculosis is known to cause detectable mycobacteremia. However, little is known about the intra-patient mycobacterial microevolution and the direction of seeding between the sputum and blood compartments.

Methods

From a diagnostic study of HIV-infected TB patients, 51 pairs of concurrent blood and sputum M. tuberculosis isolates from the same patient were available. In a previous analysis, we identified a subset with genotypic concordance, based on spoligotyping and 24 locus MIRU-VNTR. These paired isolates with identical genotypes were analyzed by whole genome sequencing and phylogenetic analysis.

Results

Of the 25 concordant pairs (49 % of the 51 paired isolates), 15 (60 %) remained viable for extraction of high quality DNA for whole genome sequencing. Two patient pairs were excluded due to poor quality sequence reads. The median CD4 cell count was 32 (IQR; 16–101)/mm³ and ten (77 %) patients were on ART. No drug resistance mutations were identified in any of the sequences analyzed. Three (23.1 %) of 13 patients had SNPs separating paired isolates from blood and sputum compartments, indicating evidence of microevolution.

Using a phylogenetic approach to identify the ancestral compartment, in two (15 %) patients the blood isolate was ancestral to the sputum isolate, in one (8 %) it was the opposite, and ten (77 %) of the pairs were identical.

Conclusions

Among HIV-infected patients with poor cellular immunity, infection with multiple strains of M. tuberculosis was found in half of the patients. In those patients with identical strains, whole genome sequencing indicated that M. tuberculosis intra-patient microevolution does occur in a few patients, yet did not reveal a consistent direction of spread between sputum and blood. This suggests that these compartments are highly connected and potentially seed each other repeatedly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1737-2) contains supplementary material, which is available to authorized users.

Mycobacterial Interspersed Repetitive Unit Can Predict Drug Resistance of Mycobacterium tuberculosis in China

Background: Recently, Mycobacterial Interspersed Repetitive Unit (MIRU) was supposed to be associated with drug resistance in Mycobacterium tuberculosis (M. tuberculosis), but whether the association exists actually in local strains in China was still unknown. This research was conducted to explore that association and the predictability of MIRU to drug resistance of Tuberculosis (TB).

Methods: The clinical isolates were collected and the susceptibility test were conducted with Lowenstein–Jensen (LJ) medium for five anti-TB drug. Based on PCR of MIRU-VNTR (Variable Number of Tandem Repeat) genotyping, we tested the number of the repeat unite of MIRU. Then, we used logistic regression to evaluate the association between 15 MIRU and drug resistance. In addition, we explored the most suitable MIRU locus of identified MIRU loci for drug resistance by multivariate logistic regression.

Results: Of the 102 strains, one isolate was resistant to rifampicin and one isolate was resistant to streptomycin. Among these fifteen MIRU, there was a association between MIRU loci polymorphism and anti-tuberculosis drug resistance, ETRB (P = 0.03, OR = 0.19, 95% CI 0.05–0.81) and ETRC (P = 0.01, OR = 0.14, 95% CI 0.03–0.64) were negatively related to isoniazid resistance; MIRU20 (P = 0.05, OR = 2.87, 95% CI 1.01–8.12) was positively associated with ethambutol resistance; and QUB11a (P = 0.02, OR = 0.79, 95% CI 0.65–0.96) was a negative association factor of p-aminosalicylic acid resistance.

Conclusion: Our research showed that MIRU loci may predict drug resistance of tuberculosis in China. However, the mechanism still needs further exploration.

Drug-resistant tuberculosis can be predicted by Mycobacterial interspersed repetitive unit locus

It is unknown whether MIRU-VNTR (Mycobacterial Interspersed Repetitive Unit-Variable Number of Tandem Repeat) is associated with drug resistance of Mycobacterium tuberculosis. The purpose of this study was to explore the ability of 24 MIRU loci to predict the drug resistance of Isoniazid (INH), Rifampicin (RFP), Streptomycin (SM), Ethambutol (EMB) and Pyrazinamide (PZA). We collected the drug resistance and MIRU loci information of 109 strains of M. tuberculosis from an open database. The results of multivariate logistic regression showed that the VNTR polymorphism of MTUB04 was related to INH resistance [odds ratio (OR) = 2.82, P = 0.00], RFP resistance (OR = 1.91, P = 0.02), SM resistance (OR = 1.98, P = 0.01) and EMB resistance (OR = 1.95, P = 0.03). MIRU40 was associated with INH resistance (OR = 2.22, P = 0.00). MTUB21 was connected with INH resistance (OR = 1.63, P = 0.02) and SM resistance (OR = 1.69, P = 0.01). MIRU26 was correlated with SM resistance (OR = 1.52, P = 0.04). MIRU39 was associated with EMB resistance (OR = 4.07, P = 0.02). The prediction power of MIRU loci were 0.84, 0.70, 0.85, and 0.74 respectively for INH (predicted by MTUB04, MIRU20, and MTUB21), RFP (predicted by MTUB04), SM (predicted by MTUB21 and MIRU26) and EMB (MTUB04 and MIRU39) through ROC analysis. Our results showed that MIRU loci were related to anti-tuberculosis drug and could predict the drug resistance of tuberculosis.

Consequences of genomic diversity in Mycobacterium tuberculosis

The causative agent of human tuberculosis, Mycobacterium tuberculosis complex (MTBC), comprises seven phylogenetically distinct lineages associated with different geographical regions. Here we review the latest findings on the nature and amount of genomic diversity within and between MTBC lineages. We then review recent evidence for the effect of this genomic diversity on mycobacterial phenotypes measured experimentally and in clinical settings. We conclude that overall, the most geographically widespread Lineage 2 (includes Beijing) and Lineage 4 (also known as Euro-American) are more virulent than other lineages that are more geographically restricted. This increased virulence is associated with delayed or reduced pro-inflammatory host immune responses, greater severity of disease, and enhanced transmission. Future work should focus on the interaction between MTBC and human genetic diversity, as well as on the environmental factors that modulate these interactions.

Multiple sampling and discriminatory fingerprinting reveals clonally complex and compartmentalized infections by M. bovis in cattle

The combination of new genotyping tools and a more exhaustive sampling policy in the analysis of infection by Mycobacterium tuberculosis has shown that infection by this pathogen is more complex than initially expected. Mixed infections, coexistence of clonal variants from a parental strain, and compartmentalized infections are all different modalities of this clonal complexity. Until recently, genotyping of Mycobacterium bovis in animal populations was based on spoligotyping and analysis of a single isolate per infection; therefore, clonal complexity is probably underdetected. We used multiple sampling combined with highly discriminatory MIRU-VNTR to study compartmentalized infections by M. bovis in a low-tuberculosis prevalence setting. We spoligotyped the M. bovis isolates from two or more anatomic locations sampled from 55 animals on 39 independent farms. Compartmentalized infections, with two different strains infecting independent lymph nodes in the same animal, were found in six cases (10.9%). MIRU-VNTR analysis confirmed that the compartmentalization was strict and that only one strain was present in each infected node. MIRU-VNTR analysis of additional infected animals on one of the farms confirmed that the compartmentalized infection was a consequence of superinfection, since the two strains were independently infecting other animals. This same analysis revealed the emergence of a microevolved clonal variant in one of the lymph nodes of the compartmentalized animal. Clonal complexity must also be taken into consideration in M. bovis infection, even in low-prevalence settings, and analyses must be adapted to detect it and increase the accuracy of molecular epidemiology studies.

Revealing hidden clonal complexity in Mycobacterium tuberculosis infection by qualitative and quantitative improvement of sampling

The analysis of microevolution events, its functional relevance and impact on molecular epidemiology strategies, constitutes one of the most challenging aspects of the study of clonal complexity in infection by Mycobacterium tuberculosis. In this study, we retrospectively evaluated whether two improved sampling schemes could provide access to the clonal complexity that is undetected by the current standards (analysis of one isolate from one sputum). We evaluated in 48 patients the analysis by mycobacterial interspersed repetitive unit-variable number tandem repeat of M. tuberculosis isolates cultured from bronchial aspirate (BAS) or bronchoalveolar lavage (BAL) and, in another 16 cases, the analysis of a higher number of isolates from independent sputum samples. Analysis of the isolates from BAS/BAL specimens revealed clonal complexity in a very high proportion of cases (5/48); in most of these cases, complexity was not detected when the isolates from sputum samples were analysed. Systematic analysis of isolates from multiple sputum samples also improved the detection of clonal complexity. We found coexisting clonal variants in two of 16 cases that would have gone undetected in the analysis of the isolate from a single sputum specimen. Our results suggest that analysis of isolates from BAS/BAL specimens is highly efficient for recording the true clonal composition of M. tuberculosis in the lungs. When these samples are not available, we recommend increasing the number of isolates from independent sputum specimens, because they might not harbour the same pool of bacteria. Our data suggest that the degree of clonal complexity in tuberculosis has been underestimated because of the deficiencies inherent in a simplified procedure.

Relaxed selection drives a noisy noncoding transcriptome in members of the Mycobacterium tuberculosis complex

Related species are often used to understand the molecular underpinning of virulence through examination of a shared set of biological features attributable to a core genome of orthologous genes. An important but insufficiently studied issue, however, is the extent to which the regulatory architectures are similarly conserved. A small number of studies have compared the primary transcriptomes of different bacterial species, but few have compared closely related species with clearly divergent evolutionary histories. We addressed the impact of differing modes of evolution within the genus Mycobacterium through comparison of the primary transcriptome of M. marinum with that of a closely related lineage, M. bovis. Both are thought to have evolved from an ancestral generalist species, with M. bovis and other members of the M. tuberculosis complex having subsequently undergone downsizing of their genomes during the transition to obligate pathogenicity. M. marinum, in contrast, has retained a large genome, appropriate for an environmental organism, and is a broad-host-range pathogen. We also examined changes over a shorter evolutionary time period through comparison of the primary transcriptome of M. bovis with that of another member of the M. tuberculosis complex (M. tuberculosis) which possesses an almost identical genome but maintains a distinct host preference.

Our comparison of the transcriptional start site (TSS) maps of M. marinum and M. bovis uncovers a pillar of conserved promoters, noncoding RNA (NCRNA), and a genome-wide signal in the −35 promoter regions of both species. We identify evolutionarily conserved transcriptional attenuation and highlight its potential contribution to multidrug resistance mediated through the transcriptional regulator whiB7. We show that a species population history is reflected in its transcriptome and posit relaxed selection as the main driver of an abundance of canonical −10 promoter sites in M. bovis relative to M. marinum. It appears that transcriptome composition in mycobacteria is driven primarily by the availability of such sites and that their frequencies diverge significantly across the mycobacterial clade. Finally, through comparison of M. bovis and M. tuberculosis, we illustrate that single nucleotide polymorphism (SNP)-driven promoter differences likely underpin many of the transcriptional differences between M. tuberculosis complex lineages.