Transmissible Mycobacterium tuberculosis Strains Share Genetic Markers and Immune Phenotypes

Genomic analysis of lineage-specific transmission of multidrug resistance tuberculosis in China

OBJECTIVES

We investigated the genetic diversities and lineage-specific transmission dynamics of multidrug-resistant tuberculosis (MDR-TB), with the goal of determining the potential factors driving the MDR epidemics in China.

METHODS

We curated a large nationwide Mycobacterium tuberculosis (M. tuberculosis) whole genome sequence data set, including 1313 MDR strains. We reconstructed the phylogeny and mapped the transmission networks of MDR-TB across China using Bayesian inference. To identify drug-resistance variants linked to enhanced transmissibility, we employed ordinary least-squares (OLS) regression analysis.

RESULT

The majority of MDR-TB strains in China belong to lineage 2.2.1. Transmission chain analysis has indicated that the repeated and frequent transmission of L2.2.1 plays a central role in the establishment of MDR epidemic in China, but no occurrence of a large predominant MDR outbreak was detected. Using OLS regression, the most common single nucleotide polymorphisms (SNPs) associated with resistance to isoniazid (katG_p.Ser315Thr and katG_p.Ser315Asn) and rifampicin (rpoB_p.Ser450Leu, rpoB_p.His445Tyr, rpoB_p.His445Arg, rpoB_p.His445Asp, and rpoB_p.His445Asn) were more likely to be found in L2 clustered strains. Several putative compensatory mutations in rpoA, rpoC, and katG were significantly associated with clustering. The eastern, central, and southern regions of China had a high level of connectivity for the migration of L2 MDR strains throughout the country. The skyline plot showed distinct population size expansion dynamics for MDR-TB lineages in China.

CONCLUSION

MDR-TB epidemic in China is predominantly driven by the spread of highly transmissible Beijing strains. A range of drug-resistance mutations of L2 MDR-TB strains displayed minimal fitness costs and may facilitate their transmission.

Iron-related gene mutations driving global Mycobacterium tuberculosis transmission revealed by whole-genome sequencing

BACKGROUND

Iron plays a crucial role in the growth of Mycobacterium tuberculosis (M. tuberculosis). However, the precise regulatory mechanism governing this system requires further elucidation. Additionally, limited studies have examined the impact of gene mutations related to iron on the transmission of M. tuberculosis globally. This research aims to investigate the correlation between mutations in iron-related genes and the worldwide transmission of M. tuberculosis.

RESULTS

A total of 13,532 isolates of M. tuberculosis were included in this study. Among them, 6,104 (45.11%) were identified as genomic clustered isolates, while 8,395 (62.04%) were classified as genomic clade isolates. Our results showed that a total of 12 single nucleotide polymorphisms (SNPs) showed a positive correlation with clustering, such as Rv1469 (ctpD, C758T), Rv3703c (etgB, G1122T), and Rv3743c (ctpJ, G676C). Additionally, seven SNPs, including Rv0104 (T167G, T478G), Rv0211 (pckA, A302C), Rv0283 (eccB3, C423T), Rv1436 (gap, G654T), ctpD C758T, and etgB C578A, demonstrated a positive correlation with transmission clades across different countries. Notably, our findings highlighted the positive association of Rv0104 T167G, pckA A302C, eccB3 C423T, ctpD C758T, and etgB C578A with transmission clades across diverse regions. Furthermore, our analysis identified 78 SNPs that exhibited significant associations with clade size.

CONCLUSIONS

Our study reveals the link between iron-related gene SNPs and M. tuberculosis transmission, offering insights into crucial factors influencing the pathogenicity of the disease. This research holds promise for targeted strategies in prevention and treatment, advancing research and interventions in this field.

Research progress of nanoparticle targeting delivery systems in bacterial infections

Bacterial infection is a huge threat to the health of human beings and animals. The abuse of antibiotics have led to the occurrence of bacterial multidrug resistance, which have become a difficult problem in the treatment of clinical infections. Given the outstanding advantages of nanodrug delivery systems in cancer treatment, many scholars have begun to pay attention to their application in bacterial infections. However, due to the similarity of the microenvironment between bacterial infection lesions and cancer sites, the targeting and accuracy of traditional microenvironment-responsive nanocarriers are questionable. Therefore, finding new specific targets has become a new development direction of nanocarriers in bacterial prevention and treatment. This article reviews the infectious microenvironment induced by bacteria and a series of virulence factors of common pathogenic bacteria and their physiological functions, which may be used as potential targets to improve the targeting accuracy of nanocarriers in lesions.

Design of a Multi-Epitope Vaccine against Tuberculosis from Mycobacterium tuberculosis PE_PGRS49 and PE_PGRS56 Proteins by Reverse Vaccinology

Tuberculosis is a disease caused by Mycobacterium tuberculosis, representing the second leading cause of death by an infectious agent worldwide. The available vaccine against this disease has insufficient coverage and variable efficacy, accounting for a high number of cases worldwide. In fact, an estimated third of the world's population has a latent infection. Therefore, developing new vaccines is crucial to preventing it. In this study, the highly antigenic PE_PGRS49 and PE_PGRS56 proteins were analyzed. These proteins were used for predicting T- and B-cell epitopes and for human leukocyte antigen (HLA) protein binding efficiency. Epitopes GGAGGNGSLSS, FAGAGGQGGLGG, GIGGGTQSATGLG (PE_PGRS49), and GTGWNGGKGDTG (PE_PGRS56) were selected based on their best physicochemical, antigenic, non-allergenic, and non-toxic properties and coupled to HLA I and HLA II structures for in silico assays. A construct with an adjuvant (RS09) plus each epitope joined by GPGPG linkers was designed, and the stability of the HLA-coupled construct was further evaluated by molecular dynamics simulations. Although experimental and in vivo studies are still necessary to ensure its protective effect against the disease, this study shows that the vaccine construct is dynamically stable and potentially effective against tuberculosis.

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.

Immune evasion and provocation by Mycobacterium tuberculosis

Mycobacterium tuberculosis, the causative agent of tuberculosis, has infected humans for millennia. M. tuberculosis is well adapted to establish infection, persist in the face of the host immune response and be transmitted to uninfected individuals. Its ability to complete this infection cycle depends on it both evading and taking advantage of host immune responses. The outcome of M. tuberculosis infection is often a state of equilibrium characterized by immunological control and bacterial persistence. Recent data have highlighted the diverse cell populations that respond to M. tuberculosis infection and the dynamic changes in the cellular and intracellular niches of M. tuberculosis during the course of infection. M. tuberculosis possesses an arsenal of protein and lipid effectors that influence macrophage functions and inflammatory responses; however, our understanding of the role that specific bacterial virulence factors play in the context of diverse cellular reservoirs and distinct infection stages is limited. In this Review, we discuss immune evasion and provocation by M. tuberculosis during its infection cycle and describe how a more detailed molecular understanding is crucial to enable the development of novel host-directed therapies, disease biomarkers and effective vaccines.

Mycobacterium tuberculosis transmission in Birmingham, UK, 2009-19: An observational study

Background

Over 10-years of whole-genome sequencing (WGS) of Mycobacterium tuberculosis in Birmingham presents an opportunity to explore epidemiological trends and risk factors for transmission in new detail.

Methods

Between 1st January 2009 and 15th June 2019, we obtained the first WGS isolate from every patient resident in a postcode district covered by Birmingham's centralised tuberculosis service. Data on patients' sex, country of birth, social risk-factors, anatomical locus of disease, and strain lineage were collected. Poisson harmonic regression was used to assess seasonal variation in case load and a mixed-effects multivariable Cox proportionate hazards model was used to assess risk factors for a future case arising in clusters defined by a 5 single nucleotide polymorphism (SNP) threshold, and by 12 SNPs in a sensitivity analysis.

Findings

511/1653 (31%) patients were genomically clustered with another. A seasonal variation in diagnoses was observed, peaking in spring, but only among clustered cases. Risk-factors for a future clustered case included UK-birth (aHR=2·03 (95%CI 1·35-3·04), p < 0·001), infectious (pulmonary/laryngeal/miliary) tuberculosis (aHR=3·08 (95%CI 1·98-4·78), p < 0·001), and M. tuberculosis lineage 3 (aHR=1·91 (95%CI 1·03-3·56), p = 0·041) and 4 (aHR=2·27 (95%CI 1·21-4·26), p = 0·011), vs. lineage 1. Similar results pertained to 12 SNP clusters, for which social risk-factors were also significant (aHR 1·72 (95%CI 1·02-2·93), p = 0·044). There was marked heterogeneity in transmission patterns between postcode districts.

Interpretation

There is seasonal variation in the diagnosis of genomically clustered, but not non-clustered, cases. Risk factors for clustering include UK-birth, infectious forms of tuberculosis, and infection with lineage 3 or 4.

Funding

Wellcome Trust, MRC, UKHSA.

PowerBacGWAS: a computational pipeline to perform power calculations for bacterial genome-wide association studies

Genome-wide association studies (GWAS) are increasingly being applied to investigate the genetic basis of bacterial traits. However, approaches to perform power calculations for bacterial GWAS are limited. Here we implemented two alternative approaches to conduct power calculations using existing collections of bacterial genomes. First, a sub-sampling approach was undertaken to reduce the allele frequency and effect size of a known and detectable genotype-phenotype relationship by modifying phenotype labels. Second, a phenotype-simulation approach was conducted to simulate phenotypes from existing genetic variants. We implemented both approaches into a computational pipeline (PowerBacGWAS) that supports power calculations for burden testing, pan-genome and variant GWAS; and applied it to collections of Enterococcus faecium, Klebsiella pneumoniae and Mycobacterium tuberculosis. We used this pipeline to determine sample sizes required to detect causal variants of different minor allele frequencies (MAF), effect sizes and phenotype heritability, and studied the effect of homoplasy and population diversity on the power to detect causal variants. Our pipeline and user documentation are made available and can be applied to other bacterial populations. PowerBacGWAS can be used to determine sample sizes required to find statistically significant associations, or the associations detectable with a given sample size. We recommend to perform power calculations using existing genomes of the bacterial species and population of study.

Rv0180c contributes to Mycobacterium tuberculosis cell shape and to infectivity in mice and macrophages

Mycobacterium tuberculosis, the main causative agent of human tuberculosis, is transmitted from person to person via small droplets containing very few bacteria. Optimizing the chance to seed in the lungs is therefore a major adaptation to favor survival and dissemination in the human population. Here we used TnSeq to identify genes important for the early events leading to bacterial seeding in the lungs. Beside several genes encoding known virulence factors, we found three new candidates not previously described: rv0180c, rv1779c and rv1592c. We focused on the gene, rv0180c, of unknown function. First, we found that deletion of rv0180c in M. tuberculosis substantially reduced the initiation of infection in the lungs of mice. Next, we established that Rv0180c enhances entry into macrophages through the use of complement-receptor 3 (CR3), a major phagocytic receptor for M. tuberculosis. Silencing CR3 or blocking the CR3 lectin site abolished the difference in entry between the wild-type parental strain and the Δrv0180c::km mutant. However, we detected no difference in the production of both CR3-known carbohydrate ligands (glucan, arabinomannan, mannan), CR3-modulating lipids (phthiocerol dimycocerosate), or proteins in the capsule of the Δrv0180c::km mutant in comparison to the wild-type or complemented strains. By contrast, we established that Rv0180c contributes to the functionality of the bacterial cell envelope regarding resistance to toxic molecule attack and cell shape. This alteration of bacterial shape could impair the engagement of membrane receptors that M. tuberculosis uses to invade host cells, and open a new perspective on the modulation of bacterial infectivity.

The epidemic efficiency of tuberculosis bacilli is determined by their capacity to transmit via aerosol. Currently, the bacterial functions that favor Mycobacterium tuberculosis seeding in the lung of naïve host remain mostly unknown. Here we implemented a genome-wide approach to identify M. tuberculosis mutants deficient for seeding and early replication in the lung of mice. In addition to genes known to encode virulence factors, we identified three genes not previously described. We used complementary approaches to characterize the phenotype of a M. tuberculosis mutant with insertion within the rv0180c gene. We found that this mutant is impaired for seeding in the lung of mice and for invasion and replication in human macrophages. In macrophages, the defect relies on a lack of engagement of CR3 receptor. Although we did not detect any difference between the wild type strain and the rv0180c mutant with regard to potential CR3-ligand, we found that the bacterial cell envelope is altered in the rv0180c mutant. Our study provides new insight into bacterial genes required for early interaction of M. tuberculosis with the host and perspective to understand the bacterial functions enhancing infectivity.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Globally diverse Mycobacterium tuberculosis resistance acquisition: a retrospective geographical and temporal analysis of whole genome sequences

Background

Mycobacterium tuberculosis whole genome sequencing (WGS) data can provide insights into temporal and geographical trends in resistance acquisition and inform public health interventions. We aimed to use a large clinical collection of M tuberculosis WGS and resistance phenotype data to study how, when, and where resistance was acquired on a global scale.

Methods

We did a retrospective analysis of WGS data. We curated a set of clinical M tuberculosis isolates with high-quality sequencing and culture-based drug susceptibility data (spanning four lineages and 52 countries in Africa, Asia, the Americas, and Europe) using public databases and literature curation. For inclusion, sequence quality criteria and country of origin data were required. We constructed geographical and lineage specific M tuberculosis phylogenies and used Bayesian molecular dating with BEAST, version 1.10.4, to infer the most recent common susceptible ancestor age for 4869 instances of resistance to ten drugs.

Findings

Between Jan 1, 1987, and Sept 12, 2014, of 10 299 M tuberculosis clinical isolates, 8550 were curated, of which 6099 (71%) from 15 countries met criteria for molecular dating. The number of independent resistance acquisition events was lower than the number of resistant isolates across all countries, suggesting ongoing transmission of drug resistance. Ancestral age distributions supported the presence of old resistance, 20 years or more before, in most countries. A consistent order of resistance acquisition was observed globally starting with resistance to isoniazid, but resistance ancestral age varied by country. We found a direct correlation between gross domestic product per capita and resistance age (r ²=0·47; p=0·014). Amplification of fluoroquinolone and second-line injectable resistance among multidrug-resistant isolates is estimated to have occurred very recently (median ancestral age 4·7 years [IQR 1·9-9·8] before sample collection). We found the sensitivity of commercial molecular diagnostics for second-line resistance to vary significantly by country (p<0·0003).

Interpretation

Our results highlight that both resistance transmission and amplification are contributing to disease burden globally but vary by country. The observation that wealthier nations are more likely to have old resistance (most recent common susceptible ancestor >20 years before isolation) suggests that programmatic improvements can reduce resistance amplification, but that fit resistant strains can circulate for decades subsequently implies the need for continued surveillance.

Improved Prediction of Bacterial Genotype-Phenotype Associations Using Interpretable Pangenome-Spanning Regressions

Discovery of genetic variants underlying bacterial phenotypes and the prediction of phenotypes such as antibiotic resistance are fundamental tasks in bacterial genomics. Genome-wide association study (GWAS) methods have been applied to study these relations, but the plastic nature of bacterial genomes and the clonal structure of bacterial populations creates challenges. We introduce an alignment-free method which finds sets of loci associated with bacterial phenotypes, quantifies the total effect of genetics on the phenotype, and allows accurate phenotype prediction, all within a single computationally scalable joint modeling framework. Genetic variants covering the entire pangenome are compactly represented by extended DNA sequence words known as unitigs, and model fitting is achieved using elastic net penalization, an extension of standard multiple regression. Using an extensive set of state-of-the-art bacterial population genomic data sets, we demonstrate that our approach performs accurate phenotype prediction, comparable to popular machine learning methods, while retaining both interpretability and computational efficiency. Compared to those of previous approaches, which test each genotype-phenotype association separately for each variant and apply a significance threshold, the variants selected by our joint modeling approach overlap substantially.IMPORTANCE Being able to identify the genetic variants responsible for specific bacterial phenotypes has been the goal of bacterial genetics since its inception and is fundamental to our current level of understanding of bacteria. This identification has been based primarily on painstaking experimentation, but the availability of large data sets of whole genomes with associated phenotype metadata promises to revolutionize this approach, not least for important clinical phenotypes that are not amenable to laboratory analysis. These models of phenotype-genotype association can in the future be used for rapid prediction of clinically important phenotypes such as antibiotic resistance and virulence by rapid-turnaround or point-of-care tests. However, despite much effort being put into adapting genome-wide association study (GWAS) approaches to cope with bacterium-specific problems, such as strong population structure and horizontal gene exchange, current approaches are not yet optimal. We describe a method that advances methodology for both association and generation of portable prediction models.

Homoplastic single nucleotide polymorphisms contributed to phenotypic diversity in Mycobacterium tuberculosis

Homoplastic mutations are mutations independently occurring in different clades of an organism. The homoplastic changes may be a result of convergence evolution due to selective pressures. Reports on the analysis of homoplastic mutations in Mycobacterium tuberculosis have been limited. Here we characterized the distribution of homoplastic single nucleotide polymorphisms (SNPs) among genomes of 1,170 clinical M. tuberculosis isolates. They were present in all functional categories of genes, with pe/ppe gene family having the highest ratio of homoplastic SNPs compared to the total SNPs identified in the same functional category. Among the pe/ppe genes, the homoplastic SNPs were common in a relatively small number of homologous genes, including ppe18, the protein of which is a component of a promising candidate vaccine, M72/AS01E. The homoplastic SNPs in ppe18 were particularly common among M. tuberculosis Lineage 1 isolates, suggesting the need for caution in extrapolating the results of the vaccine trial to the population where L1 is endemic in Asia. As expected, homoplastic SNPs strongly associated with drug resistance. Most of these mutations are already well known. However, a number of novel mutations associated with streptomycin resistance were identified, which warrants further investigation. A SNP in the intergenic region upstream of Rv0079 (DATIN) was experimentally shown to increase transcriptional activity of the downstream gene, suggesting that intergenic homoplastic SNPs should have effects on the physiology of the bacterial cells. Our study highlights the potential of homoplastic mutations to produce phenotypic changes. Under selective pressure and during interaction with the host, homoplastic mutations may confer advantages to M. tuberculosis and deserve further characterization.

Bayesian reconstruction of Mycobacterium tuberculosis transmission networks in a high incidence area over two decades in Malawi reveals associated risk factors and genomic variants

Understanding host and pathogen factors that influence tuberculosis (TB) transmission can inform strategies to eliminate the spread of Mycobacterium tuberculosis (Mtb). Determining transmission links between cases of TB is complicated by a long and variable latency period and undiagnosed cases, although methods are improving through the application of probabilistic modelling and whole-genome sequence analysis. Using a large dataset of 1857 whole-genome sequences and comprehensive metadata from Karonga District, Malawi, over 19 years, we reconstructed Mtb transmission networks using a two-step Bayesian approach that identified likely infector and recipient cases, whilst robustly allowing for incomplete case sampling. We investigated demographic and pathogen genomic variation associated with transmission and clustering in our networks. We found that whilst there was a significant decrease in the proportion of infectors over time, we found higher transmissibility and large transmission clusters for lineage 2 (Beijing) strains. By performing evolutionary convergence testing (phyC) and genome-wide association analysis (GWAS) on transmitting versus non-transmitting cases, we identified six loci, PPE54, accD2, PE_PGRS62, rplI, Rv3751 and Rv2077c, that were associated with transmission. This study provides a framework for reconstructing large-scale Mtb transmission networks. We have highlighted potential host and pathogen characteristics that were linked to increased transmission in a high-burden setting and identified genomic variants that, with validation, could inform further studies into transmissibility and TB eradication.

Coordinate regulation of virulence and metabolic genes by the transcription factor HbhR in Mycobacterium marinum

The heparin-binding hemagglutinin (HBHA) is a multifunctional protein involved in adherence of Mycobacterium tuberculosis to non-phagocytic cells and in the formation of intracytosolic lipid inclusions. We demonstrate that the expression of hbhA is regulated by a transcriptional repressor, named HbhR, in Mycobacterium marinum. The hbhR gene, located upstream of hbhA, was identified by screening a transposon insertion library and detailed analysis of a mutant overproducing HBHA. HbhR was found to repress both hbhA and hbhR transcription by binding to the promoter regions of both genes. Complementation restored production of HBHA. RNA-seq analyses comparing the mutant and parental strains uncovered 27 genes, including hbhA, that were repressed and 20 genes activated by HbhR. Among the former, the entire locus of genes coding for a type-VII secretion system, including esxA, esxB and pe-ppe paralogs, as well as the gene coding for PspA, present in intracellular lipid vesicles, was identified, as was katG, a gene involved in the sensitivity to isoniazid. The latter category contains genes that play a role in diverse functions, such as metabolism and resistance to oxidative conditions. Thus, HbhR appears to be a master regulator, linking the transcriptional regulation of virulence, metabolic and antibiotic sensitivity genes in M. marinum.

IL-32 and its splice variants are associated with protection against Mycobacterium tuberculosis infection and skewing of Th1/Th17 cytokines

Studies in IL‐32 transgenic mice and in vitro suggest that IL‐32 may have protective effects against Mycobacterium tuberculosis, but so far there are barely any studies in humans. We studied the role of IL‐32 and its splice variants in tuberculosis (TB) in vivo and in vitro. Blood transcriptional analysis showed lower total IL‐32 mRNA levels in pulmonary TB patients compared to patients with latent TB infection and healthy controls. Also, among Indonesian household contacts who were heavily exposed to an infectious TB patient, IL‐32 mRNA levels were higher among those who remained uninfected compared to those who became infected with M. tuberculosis. In peripheral blood mononuclear cells from healthy donors, we found that IL‐32γ, the most potent isoform, was down‐regulated upon M. tuberculosis stimulation. This decrease in IL‐32γ was mirrored by an increase of another splice variant, IL‐32β. Also, a higher IL‐32γ/IL‐32β ratio correlated with IFN‐γ production, whereas a lower ratio correlated with production of IL‐1Ra, IL‐6, and IL‐17. These data suggest that IL‐32 contributes to protection against M. tuberculosis infection, and that this effect may depend on the relative abundance of different IL‐32 isoforms.

Whole genome sequencing identifies bacterial factors affecting transmission of multidrug-resistant tuberculosis in a high-prevalence setting

Whole genome sequencing (WGS) can elucidate Mycobacterium tuberculosis (Mtb) transmission patterns but more data is needed to guide its use in high-burden settings. In a household-based TB transmissibility study in Peru, we identified a large MIRU-VNTR Mtb cluster (148 isolates) with a range of resistance phenotypes, and studied host and bacterial factors contributing to its spread. WGS was performed on 61 of the 148 isolates. We compared transmission link inference using epidemiological or genomic data and estimated the dates of emergence of the cluster and antimicrobial drug resistance (DR) acquisition events by generating a time-calibrated phylogeny. Using a set of 12,032 public Mtb genomes, we determined bacterial factors characterizing this cluster and under positive selection in other Mtb lineages. Four of the 61 isolates were distantly related and the remaining 57 isolates diverged ca. 1968 (95%HPD: 1945-1985). Isoniazid resistance arose once and rifampin resistance emerged subsequently at least three times. Emergence of other DR types occurred as recently as within the last year of sampling. We identified five cluster-defining SNPs potentially contributing to transmissibility. In conclusion, clusters (as defined by MIRU-VNTR typing) may be circulating for decades in a high-burden setting. WGS allows for an enhanced understanding of transmission, drug resistance, and bacterial fitness factors.

Genome-wide analysis of Mycobacterium tuberculosis polymorphisms reveals lineage-specific associations with drug resistance

BACKGROUND

Continuing evolution of the Mycobacterium tuberculosis (Mtb) complex genomes associated with resistance to anti-tuberculosis drugs is threatening tuberculosis disease control efforts. Both multi- and extensively drug resistant Mtb (MDR and XDR, respectively) are increasing in prevalence, but the full set of Mtb genes involved are not known. There is a need for increased sensitivity of genome-wide approaches in order to elucidate the genetic basis of anti-microbial drug resistance and gain a more detailed understanding of Mtb genome evolution in a context of widespread antimicrobial therapy. Population structure within the Mtb complex, due to clonal expansion, lack of lateral gene transfer and low levels of recombination between lineages, may be reducing statistical power to detect drug resistance associated variants.

RESULTS

To investigate the effect of lineage-specific effects on the identification of drug resistance associations, we applied genome-wide association study (GWAS) and convergence-based (PhyC) methods to multiple drug resistance phenotypes of a global dataset of Mtb lineages 2 and 4, using both lineage-wise and combined approaches. We identify both well-established drug resistance variants and novel associations; uniquely identifying associations for both lineage-specific and -combined GWAS analyses. We report 17 potential novel associations between antimicrobial resistance phenotypes and Mtb genomic variants.

CONCLUSIONS

For GWAS, both lineage-specific and -combined analyses are useful, whereas PhyC may perform better in contexts of greater diversity. Unique associations with XDR in lineage-specific analyses provide evidence of diverging evolutionary trajectories between lineages 2 and 4 in response to antimicrobial drug therapy.

Transmission phenotype of Mycobacterium tuberculosis strains is mechanistically linked to induction of distinct pulmonary pathology

In a study of household contacts (HHC), households were categorized into High (HT) and Low (LT) transmission groups based on the proportion of HHC with a positive tuberculin skin test. The Mycobacterium tuberculosis (Mtb) strains from HT and LT index cases of the households were designated Mtb-HT and Mtb-LT, respectively. We found that C3HeB/FeJ mice infected with Mtb-LT strains exhibited significantly higher bacterial burden compared to Mtb-HT strains and also developed diffused inflammatory lung pathology. In stark contrast, a significant number of mice infected with Mtb-HT strains developed caseating granulomas, a lesion type with high potential to cavitate. None of the Mtb-HT infected animals developed diffused inflammatory lung pathology. A link was observed between increased in vitro replication of Mtb-LT strains and their ability to induce significantly high lipid droplet formation in macrophages. These results support that distinct early interactions of Mtb-HT and Mtb-LT strains with macrophages and subsequent differential trajectories in pathological disease may be the mechanism underlying their transmission potential.

Global variation in bacterial strains that cause tuberculosis disease: a systematic review and meta-analysis

BACKGROUND

The host, microbial, and environmental factors that contribute to variation in tuberculosis (TB) disease are incompletely understood. Accumulating evidence suggests that one driver of geographic variation in TB disease is the local ecology of mycobacterial genotypes or strains, and there is a need for a comprehensive and systematic synthesis of these data. The objectives of this study were to (1) map the global distribution of genotypes that cause TB disease and (2) examine whether any epidemiologically relevant clinical characteristics were associated with those genotypes.

METHODS

We performed a systematic review of PubMed and Scopus to create a comprehensive dataset of human TB molecular epidemiology studies that used representative sampling techniques. The methods were developed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). We extracted and synthesized data from studies that reported prevalence of bacterial genotypes and from studies that reported clinical characteristics associated with those genotypes.

RESULTS

The results of this study are twofold. First, we identified 206 studies for inclusion in the study, representing over 200,000 bacterial isolates collected over 27 years in 85 countries. We mapped the genotypes and found that, consistent with previously published maps, Euro-American lineage 4 and East Asian lineage 2 strains are widespread, and West African lineages 5 and 6 strains are geographically restricted. Second, 30 studies also reported transmission chains and 4 reported treatment failure associated with genotypes. We performed a meta-analysis and found substantial heterogeneity across studies. However, based on the data available, we found that lineage 2 strains may be associated with increased risk of transmission chains, while lineages 5 and 6 strains may be associated with reduced risk, compared with lineage 4 strains.

CONCLUSIONS

This study provides the most comprehensive systematic analysis of the evidence for diversity in bacterial strains that cause TB disease. The results show both geographic and epidemiological differences between strains, which could inform our understanding of the global burden of TB. Our findings also highlight the challenges of collecting the clinical data required to inform TB diagnosis and treatment. We urge future national TB programs and research efforts to prioritize and reinforce clinical data collection in study designs and results dissemination.

Large-scale genomic analysis shows association between homoplastic genetic variation in Mycobacterium tuberculosis genes and meningeal or pulmonary tuberculosis

Background

Meningitis is the most severe manifestation of tuberculosis. It is largely unknown why some people develop pulmonary TB (PTB) and others TB meningitis (TBM); we examined if the genetic background of infecting M. tuberculosis strains may be relevant.

Methods

We whole-genome sequenced M. tuberculosis strains isolated from 322 HIV-negative tuberculosis patients from Indonesia and compared isolates from patients with TBM (n = 106) and PTB (n = 216). Using a phylogeny-adjusted genome-wide association method to count homoplasy events we examined phenotype-related changes at specific loci or genes in parallel branches of the phylogenetic tree. Enrichment scores for the TB phenotype were calculated on single nucleotide polymorphism (SNP), gene, and pathway level. Genetic associations were validated in an independent set of isolates.

Results

Strains belonged to the East-Asian lineage (36.0%), Euro-American lineage (61.5%), and Indo-Oceanic lineage (2.5%). We found no association between lineage and phenotype (Chi-square = 4.556; p = 0.207). Large genomic differences were observed between isolates; the minimum pairwise genetic distance varied from 17 to 689 SNPs. Using the phylogenetic tree, based on 28,544 common variable positions, we selected 54 TBM and 54 PTB isolates in terminal branch sets with distinct phenotypes. Genetic variation in Rv0218, and absence of Rv3343c, and nanK were significantly associated with disease phenotype in these terminal branch sets, and confirmed in the validation set of 214 unpaired isolates.

Conclusions

Using homoplasy counting we identified genetic variation in three separate genes to be associated with the TB phenotype, including one (Rv0218) which encodes a secreted protein that could play a role in host-pathogen interaction by altering pathogen recognition or acting as virulence effector.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4498-z) contains supplementary material, which is available to authorized users.

The Cording Phenotype of Mycobacterium tuberculosis Induces the Formation of Extracellular Traps in Human Macrophages

The causative agent of tuberculosis, Mycobacterium tuberculosis, shares several characteristics with organisms that produce biofilms during infections. One of these is the ability to form tight bundles also known as cords. However, little is known of the physiological relevance of the cording phenotype. In this study, we investigated whether cord-forming M. tuberculosis induce the formation of macrophage extracellular traps (METs) in human monocyte-derived macrophages. Macrophages have previously been shown to produce extracellular traps in response to various stimuli. We optimized bacterial culturing conditions that favored the formation of the cord-forming phenotype as verified by scanning electron microscopy. Microscopy analysis of METs formation during experimental infection of macrophages with M. tuberculosis revealed that cord-forming M. tuberculosis induced significantly more METs compared to the non-cording phenotype. Deletion of early secreted antigenic target-6 which is an important virulence factor of M. tuberculosis, abrogated the ability of the bacteria to induce METs. The release of extracellular DNA from host cells during infection may represent a defense mechanism against pathogens that are difficult to internalize, including cord-forming M. tuberculosis.