Whole genome SNP analysis suggests unique virulence factor differences of the Beijing and Manila families of Mycobacterium tuberculosis found in Hawaii

Rv0100: An essential acyl carrier protein from M. tuberculosis important in dormancy

We have identified an acyl-carrier protein, Rv0100, that is up-regulated in a dormancy model. This protein plays a critical role in the fatty acid biosynthesis pathway, which is important for energy storage and cell wall synthesis in Mycobacterium tuberculosis (MTB). Knocking out the Rv0100 gene resulted in a significant reduction of growth compared to wild-type MTB in the Wayne model of non-replicating persistence. We have also shown that Rv0100 is essential for the growth and survival of this pathogen during infection in mice and a macrophage model. Furthermore, knocking out Rv0100 disrupted the synthesis of phthiocerol dimycocerosates, the virulence-enhancing lipids produced by MTB and Mycobacterium bovis. We hypothesize that this essential gene contributes to MTB virulence in the state of latent infection. Therefore, inhibitors targeting this gene could prove to be potent antibacterial agents against this pathogen.

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

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

On the Identification of Clinically Relevant Bacterial Amino Acid Changes at the Whole Genome Level Using Auto-PSS-Genome

The identification of clinically relevant bacterial amino acid changes can be performed using different methods aimed at the identification of genes showing positively selected amino acid sites (PSS). Nevertheless, such analyses are time consuming, and the frequency of genes showing evidence for PSS can be low. Therefore, the development of a pipeline that allows the quick and efficient identification of the set of genes that show PSS is of interest. Here, we present Auto-PSS-Genome, a Compi-based pipeline distributed as a Docker image, that automates the process of identifying genes that show PSS using three different methods, namely codeML, FUBAR, and omegaMap. Auto-PSS-Genome accepts as input a set of FASTA files, one per genome, containing all coding sequences, thus minimizing the work needed to conduct positively selected sites analyses. The Auto-PSS-Genome pipeline identifies orthologous gene sets and corrects for multiple possible problems in input FASTA files that may prevent the automated identification of genes showing PSS. A FASTA file containing all coding sequences can also be given as an external global reference, thus easing the comparison of results across species, when gene names are different. In this work, we use Auto-PSS-Genome to analyse Mycobacterium leprae (that causes leprosy), and the closely related species M. haemophilum, that mainly causes ulcerating skin infections and arthritis in persons who are severely immunocompromised, and in children causes cervical and perihilar lymphadenitis. The genes identified in these two species as showing PSS may be those that are partially responsible for virulence and resistance to drugs.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

Development of small-molecule inhibitors of fatty acyl-AMP and fatty acyl-CoA ligases in Mycobacterium tuberculosis

Lipid metabolism in Mycobacterium tuberculosis (Mtb) relies on 34 fatty acid adenylating enzymes (FadDs) that can be grouped into two classes: fatty acyl-CoA ligases (FACLs) involved in lipid and cholesterol catabolism and long chain fatty acyl-AMP ligases (FAALs) involved in biosynthesis of the numerous essential and virulence-conferring lipids found in Mtb. The precise biochemical roles of many FACLs remain poorly characterized while the functionally non-redundant FAALs are much better understood. Here we describe the systematic investigation of 5'-O-[N-(alkanoyl)sulfamoyl]adenosine (alkanoyl adenosine monosulfamate, alkanoyl-AMS) analogs as potential multitarget FadD inhibitors for their antitubercular activity and biochemical selectivity towards representative FAAL and FACL enzymes. We identified several potent compounds including 12-azidododecanoyl-AMS 28, 11-phenoxyundecanoyl-AMS 32, and nonyloxyacetyl-AMS 36 with minimum inhibitory concentrations (MICs) against M. tuberculosis ranging from 0.098 to 3.13 μM. Compound 32 was notable for its impressive biochemical selectivity against FAAL28 (apparent Ki = 0.7 μM) versus FACL19 (Ki > 100 μM), and uniform activity against a panel of multidrug and extensively drug-resistant TB strains with MICs ranging from 3.13 to 12.5 μM in minimal (GAST) and rich (7H9) media. The SAR analysis provided valuable insights for further optimization of 32 and also identified limitations to overcome.

Genomic sequencing is required for identification of tuberculosis transmission in Hawaii

BACKGROUND

Tuberculosis (TB) caused an estimated 1.4 million deaths and 10.4 million new cases globally in 2015. TB rates in the United States continue to steadily decline, yet rates in the State of Hawaii are perennially among the highest in the nation due to a continuous influx of immigrants from the Western Pacific and Asia. TB in Hawaii is composed of a unique distribution of genetic lineages, with the Beijing and Manila families of Mycobacterium tuberculosis (Mtb) comprising over two-thirds of TB cases. Standard fingerprinting methods (spoligotyping plus 24-loci Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats [MIRU-VNTR] fingerprinting) perform poorly when used to identify actual transmission clusters composed of isolates from these two families. Those typing methods typically group isolates from these families into large clusters of non-linked isolates with identical fingerprints. Next-generation whole-genome sequencing (WGS) provides a new tool for molecular epidemiology that can resolve clusters of isolates with identical spoligotyping and MIRU-VNTR fingerprints.

METHODS

We performed WGS and SNP analysis and evaluated epidemiological data to investigate 19 apparent TB transmission clusters in Hawaii from 2003 to 2017 in order to assess WGS' ability to resolve putative Mtb clusters from the Beijing and Manila families. This project additionally investigated MIRU-VNTR allele prevalence to determine why standard Mtb fingerprinting fails to usefully distinguish actual transmission clusters from these two Mtb families.

RESULTS

WGS excluded transmission events in seven of these putative clusters, confirmed transmission in eight, and identified both transmission-linked and non-linked isolates in four. For epidemiologically identified clusters, while the sensitivity of MIRU-VNTR fingerprinting for identifying actual transmission clusters was found to be 100%, its specificity was only 28.6% relative to WGS. We identified that the Beijing and Manila families' significantly lower Shannon evenness of MIRU-VNTR allele distributions than lineage 4 was the cause of standard fingerprinting's poor performance when identifying transmission in Beijing and Manila family clusters.

CONCLUSIONS

This study demonstrated that WGS is necessary for epidemiological investigation of TB in Hawaii and the Pacific.

Virulence of Mycobacterium tuberculosis Clinical Isolates Is Associated With Sputum Pre-treatment Bacterial Load, Lineage, Survival in Macrophages, and Cytokine Response

It is uncertain whether differences in Mycobacterium tuberculosis (Mtb) virulence defined in vitro influence clinical tuberculosis pathogenesis, transmission, and mortality. We primarily used a macrophage lysis model to characterize the virulence of Mtb isolates collected from 153 Vietnamese adults with pulmonary tuberculosis. The virulence phenotypes were then investigated for their relationship with sputum bacterial load, bacterial lineages, bacterial growth, and cytokine responses in macrophages. Over 6 days of infection, 34 isolates (22.2%) showed low virulence (< 5% macrophages lysed), 46 isolates (30.1%) showed high virulence (≥90% lysis of macrophages), and 73 isolates (47.7%) were of intermediate virulence (5–90% macrophages lysed). Highly virulent isolates were associated with an increased bacterial load in patients' sputum before anti-tuberculosis therapy (P = 0.02). Isolate-dependent virulence phenotype was consistent in both THP-1 and human monocyte-derived macrophages. High virulence isolates survived better and replicated in macrophages one hundred fold faster than those with low virulence. Macrophages infected with high virulence isolates produced lower concentrations of TNF-α and IL-6 (P = 0.002 and 0.0005, respectively), but higher concentration of IL-1β (P = 5.1 × 10⁻⁵) compared to those infected with low virulence isolates. High virulence was strongly associated with East Asian/Beijing lineage [P = 0.002, Odd ratio (OR) = 4.32, 95% confident intervals (CI) 1.68–11.13]. The association between virulence phenotypes, bacterial growth, and proinflammatory cytokines in macrophages suggest the suppression of certain proinflammatory cytokines (TNF-α and IL-6) but not IL-1β allows better intracellular survival of highly virulent Mtb. This could result in rapid macrophage lysis and higher bacterial load in sputum of patients infected with high virulence isolates, which may contribute to the pathogenesis and success of the Beijing lineage.