Comparative analysis of Mycobacterium tuberculosis pe and ppe genes reveals high sequence variation and an apparent absence of selective constraints

PPE59 antibodies in tuberculous patients and potential use for diagnosis when assayed with other rapid biomarkers

BACKGROUND

PPE 59, which is absent from bacillus Calmette Guérin (BCG) strains, seems to induce a humoral immune response in patients with tuberculosis (TB). Additional studies are needed to better evaluate this protein in immune response to tuberculosis.

OBJECTIVES

To evaluate the response of antibodies to PPE59 in TB individuals, its combination with IgG response to other, previously tested mycobacterial antigens (Ag) and with sputum smear microbiology (SM) results.

METHODS

We have cloned and expressed the rv3429 gene that encodes PPE59, then IgG, IgM, and IgA against PPE59 antigens measured by enzyme-linked immunosorbent assay (ELISA) in 212 sera samples obtained from the following subject cohorts: TB residents from Italy (79) and in Brazil (52); and an all-Brazilian cohort of 55 patients with other respiratory disorders; 10 patients infected with non-tuberculous mycobacteria, and 16 asymptomatic subjects. Drawing on results from a previous study(17) of serum samples from Brazilian subjects tested for IgG by ELISA against mycobacterial antigens ESAT-6, 16kDa, MT10.3, MPT-64 and 38kDa, the results were analysed in combination with those of the PPE59 and SM tests.

FINDINGS

Keeping the specificity rate at 97%, the overall PPE59 IgA sensitivity was 42.7%, while IgG and IgM showed lower performance (p < 0.0001). Combining PPE59 IgA/16kDa IgG results increased sensitivity to 71%, and even higher rates when the results were combined with SM results (86.5%, p = 0.001), at 88.9% specificity. Positive IgA was associated with pulmonary image alterations of high TB probability (p < 0.05).

MAIN CONCLUSIONS

Tests with TB patients found a moderate frequency of positivity for PPE59 IgA. However, the higher level of sensitivity attained in combination with PPE59 IgA/16kDa IgG/SM results unheard of before, although imperfect, suggests that this may be a potential additional tool for rapid detection of TB in low-resource areas.

PE/PPE Proteome and ESX-5 Substrate Spectrum in Mycobacterium marinum

PE/PPE proteins secreted by the ESX-5 type VII secretion system constitute a major protein repertoire in pathogenic mycobacteria and are essential for bacterial survival, pathogenicity, and host-pathogen interaction; however, little is known about their expression and secretion. The scarcity of arginine and lysine residues in PE/PPE protein sequences and the high homology of their N-terminal domains limit protein identification using classical trypsin-based proteomic methods. This study used endoproteinase AspN and trypsin to characterize the proteome of Mycobacterium marinum. Twenty-seven PE/PPE proteins were uniquely identified in AspN digests, especially PE_PGRS proteins. These treatments allowed the identification of approximately 50% of the PE/PPE pool encoded in the genome. Moreover, EspG5 pulldown assays retrieved 44 ESX-5-associated PPE proteins, covering 85% of the PPE pool in the identified proteome. The identification of PE/PE_PGRS proteins in the EspG5 interactome suggested the presence of PE-PPE pairs. The correlation analysis between protein abundance and phylogenetic relationships found potential PE/PPE pairs, indicating the presence of multiple PE/PE_PGRS partners in one PPE. We validated that EspG5 interacted with PPE31 and PPE32 and mapped critical residues for complex formation. The modified proteomic platform increases the coverage of PE/PPE proteins and elucidates the expression and localization of these proteins.

Multi-oligomeric and catalytically compromised serine acetyltransferase and cysteine regulatory complex of Mycobacterium tuberculosis

l-cysteine, a primary building block of mycothiol, plays an essential role in the defense mechanism of Mycobacterium tuberculosis (Mtb). However, it is unclear how Mtb regulates cysteine biosynthesis as no study has reported the cysteine regulatory complex (CRC) in Mtb. Serine acetyltransferase (SAT) and cysteine synthase (CS) interact to form CRC. Although MtCS has been characterized well, minimal information is available on MtSAT, which synthesizes, O-acetylserine (OAS), the precursor of cysteine. This study fills the gap and provides experimental evidence for the presence of MtCRC and a non-canonical multi-oligomeric MtSAT. We employed multiple analytical methods to characterize the oligomeric and kinetic properties of MtSAT and MtCRC. Results show that MtSAT, lacking >75 N-terminal amino acids exists in three different assembly states; trimer, hexamer, and dodecamer, compared to the single hexameric state of SAT of other bacteria. While hexamers display the highest catalytic turnover, the trimer is the least active. The predominance of trimers at low physiologically relevant concentrations suggests that MtSAT displays the lowest catalytic potential known. Further, the catalytic potential of MtSAT is also significantly reduced in CRC state, in contrast to enhanced activity of SAT in CRC of other organisms. Our study provides insights into multi-oligomeric MtSAT with reduced catalytic potential and demonstrates that both MtSAT and MtCS of Mycobacterium interact to form CRC, although with altered catalytic properties. We discuss our results in light of the altered biochemistry of the last step of canonical sulfate-dependent cysteine biosynthesis of Mycobacterium.

Advancing animal tuberculosis surveillance using culture-independent long-read whole-genome sequencing

Animal tuberculosis is a significant infectious disease affecting both livestock and wildlife populations worldwide. Effective disease surveillance and characterization of Mycobacterium bovis (M. bovis) strains are essential for understanding transmission dynamics and implementing control measures. Currently, sequencing of genomic information has relied on culture-based methods, which are time-consuming, resource-demanding, and concerning in terms of biosafety. This study explores the use of culture-independent long-read whole-genome sequencing (WGS) for a better understanding of M. bovis epidemiology in African buffaloes (Syncerus caffer). By comparing two sequencing approaches, we evaluated the efficacy of Illumina WGS performed on culture extracts and culture-independent Oxford Nanopore adaptive sampling (NAS). Our objective was to assess the potential of NAS to detect genomic variants without sample culture. In addition, culture-independent amplicon sequencing, targeting mycobacterial-specific housekeeping and full-length 16S rRNA genes, was applied to investigate the presence of microorganisms, including nontuberculous mycobacteria. The sequencing quality obtained from DNA extracted directly from tissues using NAS is comparable to the sequencing quality of reads generated from culture-derived DNA using both NAS and Illumina technologies. We present a new approach that provides complete and accurate genome sequence reconstruction, culture independently, and using an economically affordable technique.

The MAGMA pipeline for comprehensive genomic analyses of clinical Mycobacterium tuberculosis samples

BACKGROUND

Whole genome sequencing (WGS) holds great potential for the management and control of tuberculosis. Accurate analysis of samples with low mycobacterial burden, which are characterized by low (<20x) coverage and high (>40%) levels of contamination, is challenging. We created the MAGMA (Maximum Accessible Genome for Mtb Analysis) bioinformatics pipeline for analysis of clinical Mtb samples.

METHODS AND RESULTS

High accuracy variant calling is achieved by using a long seedlength during read mapping to filter out contaminants, variant quality score recalibration with machine learning to identify genuine genomic variants, and joint variant calling for low Mtb coverage genomes. MAGMA automatically generates a standardized and comprehensive output of drug resistance information and resistance classification based on the WHO catalogue of Mtb mutations. MAGMA automatically generates phylogenetic trees with drug resistance annotations and trees that visualize the presence of clusters. Drug resistance and phylogeny outputs from sequencing data of 79 primary liquid cultures were compared between the MAGMA and MTBseq pipelines. The MTBseq pipeline reported only a proportion of the variants in candidate drug resistance genes that were reported by MAGMA. Notable differences were in structural variants, variants in highly conserved rrs and rrl genes, and variants in candidate resistance genes for bedaquiline, clofazmine, and delamanid. Phylogeny results were similar between pipelines but only MAGMA visualized clusters.

CONCLUSION

The MAGMA pipeline could facilitate the integration of WGS into clinical care as it generates clinically relevant data on drug resistance and phylogeny in an automated, standardized, and reproducible manner.

Functional genetic variation in pe/ppe genes contributes to diversity in Mycobacterium tuberculosis lineages and potential interactions with the human host

Introduction

Around 10% of the coding potential of Mycobacterium tuberculosisis constituted by two poorly understood gene families, the pe and ppe loci, thought to be involved in host-pathogen interactions. Their repetitive nature and high GC content have hindered sequence analysis, leading to exclusion from whole-genome studies. Understanding the genetic diversity of pe/ppe families is essential to facilitate their potential translation into tools for tuberculosis prevention and treatment.

Methods

To investigate the genetic diversity of the 169 pe/ppe genes, we performed a sequence analysis across 73 long-read assemblies representing seven different lineages of M. tuberculosis and M. bovis BCG. Individual pe/ppe gene alignments were extracted and diversity and conservation across the different lineages studied.

Results

The pe/ppe genes were classified into three groups based on the level of protein sequence conservation relative to H37Rv, finding that >50% were conserved, with indels in pe_pgrs and ppe_mptr sub-families being major drivers of structural variation. Gene rearrangements, such as duplications and gene fusions, were observed between pe and pe_pgrs genes. Inter-lineage diversity revealed lineage-specific SNPs and indels.

Discussion

The high level of pe/ppe genes conservation, together with the lineage-specific findings, suggest their phylogenetic informativeness. However, structural variants and gene rearrangements differing from the reference were also identified, with potential implications for pathogenicity. Overall, improving our knowledge of these complex gene families may have insights into pathogenicity and inform the development of much-needed tools for tuberculosis control.

N-terminal PPE domain plays an integral role in extracellular transportation and stability of the immunomodulatory Rv3539 protein of the Mycobacterium tuberculosis

Multigene PE/PPE family is exclusively present in mycobacterium species. Only few selected genes of this family have been characterized till date. Rv3539 was annotated as PPE63 with conserved PPE domain at N-terminal and PE-PPE at C-terminal. An α/β hydrolase structural fold, characteristic of lipase/esterase, was present in the PE-PPE domain. To assign the biochemical function to Rv3539, the corresponding gene was cloned in pET-32a (+) as full-length, PPE, and PE-PPE domains individually, followed by expression in E. Coli C41 (DE3). All three proteins demonstrated esterase activity. However, the enzyme activity in the N-terminal PPE domain was very low. The enzyme activity of Rv3539 and PE-PPE proteins was approximately same with the pNP-C4 as optimum substrate at 40 °C and pH 8.0. The loss of enzyme activity after mutating the predicted catalytic triad (Ser296Ala, Asp369Ala, and His395Ala) found only in the PE-PPE domain, confirmed the candidature of the bioinformatically predicted active site residue. The optimal activity and thermostability of the Rv3539 protein was altered by removing the PPE domain. CD-spectroscopy analysis confirmed the role of PPE domain to the thermostability of Rv3539 by maintaining the structural integrity at higher temperatures. The presence of the N-terminal PPE domain directed the Rv3539 protein to the cell membrane/wall and the extracellular compartment. The Rv3539 protein could generate humoral response in TB patients. Therefore, results demonstrated that Rv3539 demonstrated esterase activity. PE-PPE domain of Rv3539 is functionally automated, however, N-terminus domain played a role in protein stabilization and its transportation. Both domains participated in immunomodulation.

Digerati - A multipath parallel hybrid deep learning framework for the identification of mycobacterial PE/PPE proteins

The genome of Mycobacterium tuberculosis contains a relatively high percentage (10%) of genes that are poorly characterised because of their highly repetitive nature and high GC content. Some of these genes encode proteins of the PE/PPE family, which are thought to be involved in host-pathogen interactions, virulence, and disease pathogenicity. Members of this family are genetically divergent and challenging to both identify and classify using conventional computational tools. Thus, advanced in silico methods are needed to identify proteins of this family for subsequent functional annotation efficiently. In this study, we developed the first deep learning-based approach, termed Digerati, for the rapid and accurate identification of PE and PPE family proteins. Digerati was built upon a multipath parallel hybrid deep learning framework, which equips multi-layer convolutional neural networks with bidirectional, long short-term memory, equipped with a self-attention module to effectively learn the higher-order feature representations of PE/PPE proteins. Empirical studies demonstrated that Digerati achieved a significantly better performance (∼18-20%) than alignment-based approaches, including BLASTP, PHMMER, and HHsuite, in both prediction accuracy and speed. Digerati is anticipated to facilitate community-wide efforts to conduct high-throughput identification and analysis of PE/PPE family members. The webserver and source codes of Digerati are publicly available at http://web.unimelb-bioinfortools.cloud.edu.au/Digerati/.

Attenuated but immunostimulatory Mycobacterium tuberculosis variant bovis strain Ravenel shows variation in T cell epitopes

Tuberculosis, caused by Mycobacterium tuberculosis complex (MTBC) organisms, affects a range of humans and animals globally. Mycobacterial pathogenesis involves manipulation of the host immune system, partially through antigen presentation. Epitope sequences across the MTBC are evolutionarily hyperconserved, suggesting their recognition is advantageous for the bacterium. Mycobacterium tuberculosis var. bovis (MBO) strain Ravenel is an isolate known to provoke a robust immune response in cattle, but typically fails to produce lesions and persist. Unlike attenuated MBO BCG strains that lack the critical RD1 genomic region, Ravenel is classic-type MBO structurally, suggesting genetic variation is responsible for defective pathogenesis. This work explores variation in epitope sequences in MBO Ravenel by whole genome sequencing, and contrasts such variation against a fully virulent clinical isolate, MBO strain 10-7428. Validated MTBC epitopes (n = 4818) from the Immune Epitope Database were compared to their sequences in MBO Ravenel and MBO 10-7428. Ravenel yielded 3 modified T cell epitopes, in genes rpfB, argC, and rpoA. These modifications were predicted to have little effect on protein stability. In contrast, T cells epitopes in 10-7428 were all WT. Considering T cell epitope hyperconservation across MTBC variants, these altered MBO Ravenel epitopes support their potential contribution to overall strain attenuation. The affected genes may provide clues on basic pathogenesis, and if so, be feasible targets for reverse vaccinology.

Structural Basis of PE_PGRS Polymorphism, a Tool for Functional Modulation

BACKGROUND

The mycobacterial PE_PGRS protein family is present only in pathogenic strains of the genus mycobacterium, such as Mtb and members of the MTB complex, suggesting a likely important role of this family in pathogenesis. Their PGRS domains are highly polymorphic and have been suggested to cause antigenic variations and facilitate pathogen survival. The availability of AlphaFold2.0 offered us a unique opportunity to better understand structural and functional properties of these domains and a role of polymorphism in Mtb evolution and dissemination.

METHODS

We made extensive use of AlphaFold2.0 computations and coupled them with sequence distribution phylogenetic and frequency analyses, and antigenic predictions.

RESULTS

Modeling of several polymorphic forms of PE_PGRS33, the prototype of the PE_PGRS family and sequence analyses allowed us to predict the structural impact of mutations/deletions/insertions present in the most frequent variants. These analyses well correlate with the observed frequency and with the phenotypic features of the described variants.

CONCLUSIONS

Here, we provide a thorough description of structural impacts of the observed polymorphism of PE_PGRS33 protein and we correlate predicted structures to the known fitness of strains containing specific variants. Finally, we also identify protein variants associated with bacterial evolution, showing sophisticated modifications likely endowed with a gain-of-function role during bacterial evolution.

The PE-PPE Family of Mycobacterium tuberculosis: Proteins in Disguise

Mycobacterium tuberculosis has thrived in parallel with humans for millennia, and despite our efforts, M. tuberculosis continues to plague us, currently infecting a third of the world's population. The success of M. tuberculosis has recently been attributed, in part, to the PE-PPE family; a unique collection of 168 proteins fundamentally involved in the pathogenesis of M. tuberculosis. The PE-PPE family proteins have been at the forefront of intense research efforts since their discovery in 1998 and whilst our knowledge and understanding has significantly advanced over the last two decades, many important questions remain to be elucidated. This review consolidates and examines the vast body of existing literature regarding the PE-PPE family proteins, with respect to the latest developments in elucidating their evolution, structure, subcellular localisation, function, and immunogenicity. This review also highlights significant inconsistencies and contradictions within the field. Additionally, possible explanations for these knowledge gaps are explored. Lastly, this review poses many important questions, which need to be addressed to complete our understanding of the PE-PPE family, as well as highlighting the challenges associated with studying this enigmatic family of proteins. Further research into the PE-PPE family, together with technological advancements in genomics and proteomics, will undoubtedly improve our understanding of the pathogenesis of M. tuberculosis, as well as identify key targets/candidates for the development of novel drugs, diagnostics, and vaccines.

ppe51 Variants Enable Growth of Mycobacterium tuberculosis at Acidic pH by Selectively Promoting Glycerol Uptake

In defined media supplemented with single carbon sources, Mycobacterium tuberculosis (Mtb) exhibits carbon source specific growth restriction. When supplied with glycerol as the sole carbon source at pH 5.7, Mtb establishes a metabolically active state of nonreplicating persistence known as acid growth arrest. We hypothesized that acid growth arrest on glycerol is not a metabolic restriction, but rather an adaptive response. To test this hypothesis, we selected for and identified several Mtb mutants that could grow under these restrictive conditions. All mutations were mapped to the ppe51 gene and resulted in variants with 3 different amino acid substitutions- S211R, E215K, and A228D. Expression of the ppe51 variants in Mtb promoted growth at acidic pH showing that the mutant alleles are sufficient to cause the dominant gain-of-function, Enhanced Acid Growth (EAG) phenotype. Testing growth on other single carbon sources showed the PPE51 variants specifically enhanced growth on glycerol, suggesting PPE51 plays a role in glycerol uptake. Using radiolabeled glycerol, enhanced glycerol uptake was observed in Mtb expressing the PPE51 (S211R) variant, with glycerol overaccumulation in triacylglycerol. Notably, the EAG phenotype is deleterious for growth in macrophages, where the mutants have selectively faster replication and reduced survival in activated macrophages compared to resting macrophages. Recombinant PPE51 protein exhibited differential thermostability in the wild type (WT) or S211R variants in the presence of glycerol, supporting the model that EAG substitutions alter PPE51-glycerol interactions. Together, these findings support that PPE51 variants selectively promote glycerol uptake and that slowed growth at acidic pH is an important adaptive mechanism required for macrophage pathogenesis. IMPORTANCE It is puzzling why Mycobacterium tuberculosis (Mtb) cannot grow on glycerol at acidic pH, as it has a carbon source and oxygen, everything it needs to grow. In this study, we found that Mtb limits uptake of glycerol at acidic pH to restrict its growth and that mutations in ppe51 promote uptake of glycerol at acidic pH and enable growth. That is, Mtb can grow well at acidic pH on glycerol, but has adapted instead to stop growth. Notably, ppe51 variants exhibit enhanced replication and reduced survival in activated macrophages, supporting a role for pH-dependent slowed growth during macrophage pathogenesis.

MycoVarP: Mycobacterium Variant and Drug Resistance Prediction Pipeline for Whole-Genome Sequence Data Analysis

Whole-genome sequencing (WGS) provides a comprehensive tool to analyze the bacterial genomes for genotype–phenotype correlations, diversity of single-nucleotide variant (SNV), and their evolution and transmission. Several online pipelines and standalone tools are available for WGS analysis of Mycobacterium tuberculosis (Mtb) complex (MTBC). While they facilitate the processing of WGS data with minimal user expertise, they are either too general, providing little insights into bacterium-specific issues such as gene variations, INDEL/synonymous/PE-PPE (IDP family), and drug resistance from sample data, or are limited to specific objectives, such as drug resistance. It is understood that drug resistance and lineage-specific issues require an elaborate prioritization of identified variants to choose the best target for subsequent therapeutic intervention. Mycobacterium variant pipeline (MycoVarP) addresses these specific issues with a flexible battery of user-defined and default filters. It provides an end-to-end solution for WGS analysis of Mtb variants from the raw reads and performs two quality checks, viz, before trimming and after alignments of reads to the reference genome. MycoVarP maps the annotated variants to the drug-susceptible (DS) database and removes the false-positive variants, provides lineage identification, and predicts potential drug resistance. We have re-analyzed the WGS data reported by Advani et al. (2019) using MycoVarP and identified some additional variants not reported so far. We conclude that MycoVarP will help in identifying nonsynonymous, true-positive, drug resistance–associated variants more effectively and comprehensively, including those within the IDP of the PE-PPE/PGRS family, than possible from the currently available pipelines.

Identification of covalent inhibitors that disrupt M. tuberculosis growth by targeting multiple serine hydrolases involved in lipid metabolism

The increasing incidence of antibiotic-resistant Mycobacterium tuberculosis infections is a global health threat necessitating the development of new antibiotics. Serine hydrolases (SHs) are a promising class of targets because of their importance for the synthesis of the mycobacterial cell envelope. We screen a library of small molecules containing serine-reactive electrophiles and identify narrow-spectrum inhibitors of M. tuberculosis growth. Using these lead molecules, we perform competitive activity-based protein profiling and identify multiple SH targets, including enzymes with uncharacterized functions. Lipidomic analyses of compound-treated cultures reveal an accumulation of free lipids and a substantial decrease in lipooligosaccharides, linking SH inhibition to defects in cell envelope biogenesis. Mutant analysis reveals a path to resistance via the synthesis of mycocerates, but not through mutations to SH targets. Our results suggest that simultaneous inhibition of multiple SH enzymes is likely to be an effective therapeutic strategy for the treatment of M. tuberculosis infections.

Early alveolar macrophage response and IL-1R-dependent T cell priming determine transmissibility of Mycobacterium tuberculosis strains

Mechanisms underlying variability in transmission of Mycobacterium tuberculosis strains remain undefined. By characterizing high and low transmission strains of M.tuberculosis in mice, we show here that high transmission M.tuberculosis strain induce rapid IL-1R-dependent alveolar macrophage migration from the alveolar space into the interstitium and that this action is key to subsequent temporal events of early dissemination of bacteria to the lymph nodes, Th1 priming, granulomatous response and bacterial control. In contrast, IL-1R-dependent alveolar macrophage migration and early dissemination of bacteria to lymph nodes is significantly impeded in infection with low transmission M.tuberculosis strain; these events promote the development of Th17 immunity, fostering neutrophilic inflammation and increased bacterial replication. Our results suggest that by inducing granulomas with the potential to develop into cavitary lesions that aids bacterial escape into the airways, high transmission M.tuberculosis strain is poised for greater transmissibility. These findings implicate bacterial heterogeneity as an important modifier of TB disease manifestations and transmission.

Structure-Aware Mycobacterium tuberculosis Functional Annotation Uncloaks Resistance, Metabolic, and Virulence Genes

Accurate and timely functional genome annotation is essential for translating basic pathogen research into clinically impactful advances. Here, through literature curation and structure-function inference, we systematically update the functional genome annotation of Mycobacterium tuberculosis virulent type strain H37Rv. First, we systematically curated annotations for 589 genes from 662 publications, including 282 gene products absent from leading databases. Second, we modeled 1,711 underannotated proteins and developed a semiautomated pipeline that captured shared function between 400 protein models and structural matches of known function on Protein Data Bank, including drug efflux proteins, metabolic enzymes, and virulence factors. In aggregate, these structure- and literature-derived annotations update 940/1,725 underannotated H37Rv genes and generate hundreds of functional hypotheses. Retrospectively applying the annotation to a recent whole-genome transposon mutant screen provided missing function for 48% (13/27) of underannotated genes altering antibiotic efficacy and 33% (23/69) required for persistence during mouse tuberculosis (TB) infection. Prospective application of the protein models enabled us to functionally interpret novel laboratory generated pyrazinamide (PZA)-resistant mutants of unknown function, which implicated the emerging coenzyme A depletion model of PZA action in the mutants’ PZA resistance. Our findings demonstrate the functional insight gained by integrating structural modeling and systematic literature curation, even for widely studied microorganisms. Functional annotations and protein structure models are available at https://tuberculosis.sdsu.edu/H37Rv in human- and machine-readable formats.

IMPORTANCEMycobacterium tuberculosis, the primary causative agent of tuberculosis, kills more humans than any other infectious bacterium. Yet 40% of its genome is functionally uncharacterized, leaving much about the genetic basis of its resistance to antibiotics, capacity to withstand host immunity, and basic metabolism yet undiscovered. Irregular literature curation for functional annotation contributes to this gap. We systematically curated functions from literature and structural similarity for over half of poorly characterized genes, expanding the functionally annotated Mycobacterium tuberculosis proteome. Applying this updated annotation to recent in vivo functional screens added functional information to dozens of clinically pertinent proteins described as having unknown function. Integrating the annotations with a prospective functional screen identified new mutants resistant to a first-line TB drug, supporting an emerging hypothesis for its mode of action. These improvements in functional interpretation of clinically informative studies underscore the translational value of this functional knowledge. Structure-derived annotations identify hundreds of high-confidence candidates for mechanisms of antibiotic resistance, virulence factors, and basic metabolism and other functions key in clinical and basic tuberculosis research. More broadly, they provide a systematic framework for improving prokaryotic reference annotations.

Proline-Glutamate/Proline-Proline-Glutamate (PE/PPE) proteins of Mycobacterium tuberculosis: The multifaceted immune-modulators

The PE/PPE proteins encoded by seven percent (7%) of Mycobacterium tuberculosis (Mtb) genome are the chief constituents to pathogen's virulence reservoir. The fact that these genes have evolved along ESX secretory system in pathogenic Mtb strains make their investigation very intriguing. There is lot of speculation about the prominent role of these proteins at host pathogen interface and in disease pathogenesis. Nevertheless, the exact function of PE/PPE proteins still remains a mystery which calls for further research targeting these proteins. This article is an effort to document all the facts known so far with regard to these unique proteins which involves their origin, evolution, transcriptional control, and most important their role as host immune-modulators. Our understanding strongly points towards the versatile nature of these PE/PPE proteins as Mtb's host immune sensors and as decisive factors in shaping the outcome of infection. Further investigation on these proteins will surely pave way for newer and effective vaccines and therapeutics to control Tuberculosis (TB).

Whole-Genome Sequencing of Mycobacterium tuberculosis Directly from Sputum Samples

Whole-genome sequencing is a powerful, high-resolution tool that can be used to generate accurate data on bacterial population structure, phylogeography, and mutations associated with antimicrobial resistance. The ability to sequence pathogen genomes directly from clinical specimens, without the requirement for in vitro culturing, is attractive in terms of time- and labor-saving, especially in the case of slow growing pathogens, such as Mycobacterium tuberculosis. However, clinical samples typically contain too low levels of pathogen nucleic acid, plus relatively high levels of human and natural microbiota DNA/RNA, to make this a viable option. Using a combination of whole-genome enrichment and deep sequencing, which has been proven to be a nonmutagenic approach, we can capture all known variations found within M. tuberculosis genomes. The method is a consistent and sensitive tool that enables rapid whole-genome sequencing of M. tuberculosis directly from clinical samples and has the potential to be adapted to other pathogens with a similar clonal nature.

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.

Correlation of drug resistance with single nucleotide variations through genome analysis and experimental validation in a multi-drug resistant clinical isolate of M. tuberculosis

BACKGROUND

Genome sequencing and genetic polymorphism analysis of clinical isolates of M. tuberculosis is carried out to gain further insight into molecular pathogenesis and host-pathogen interaction. Therefore the functional evaluation of the effect of single nucleotide variation (SNV) is essential. At the same time, the identification of invariant sequences unique to M. tuberculosis contributes to infection detection by sensitive methods. In the present study, genome analysis is accompanied by evaluation of the functional implication of the SNVs in a MDR clinical isolate VPCI591.

RESULT

By sequencing and comparative analysis of VPCI591 genome with 1553 global clinical isolates of M. tuberculosis (GMTV and tbVar databases), we identified 141 unique strain specific SNVs. A novel intergenic variation in VPCI591 in the putative promoter/regulatory region mapping between embC (Rv3793) and embA (Rv3794) genes was found to enhance the expression of embAB, which correlates with the high resistance of the VPCI591 to ethambutol. Similarly, the unique combination of three genic SNVs in RNA polymerase β gene (rpoB) in VPCI591 was evaluated for its effect on rifampicin resistance through molecular docking analysis. The comparative genomics also showed that along with variations, there are genes that remain invariant. 173 such genes were identified in our analysis.

CONCLUSION

The genetic variation in M. tuberculosis clinical isolate VPCI591 is found in almost all functional classes of genes. We have shown that SNV in rpoB gene mapping outside the drug binding site along with two SNVs in the binding site can contribute to quantitative change in MIC for rifampicin. Our results show the collective effect of SNVs on the structure of the protein, impacting the interaction between the target protein and the drug molecule in rpoB as an example. The study shows that intergenic variations bring about quantitative changes in transcription in embAB and in turn can lead to drug resistance.

Comparison of routine field epidemiology and whole genome sequencing to identify tuberculosis transmission in a remote setting

Yukon Territory (YT) is a remote region in northern Canada with ongoing spread of tuberculosis (TB). To explore the utility of whole genome sequencing (WGS) for TB surveillance and monitoring in a setting with detailed contact tracing and interview data, we used a mixed-methods approach. Our analysis included all culture-confirmed cases in YT (2005-2014) and incorporated data from 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping, WGS and contact tracing. We compared field-based (contact investigation (CI) data + MIRU-VNTR) and genomic-based (WGS + MIRU-VNTR + basic case data) investigations to identify the most likely source of each person's TB and assessed the knowledge, attitudes and practices of programme personnel around genotyping and genomics using online, multiple-choice surveys (n = 4) and an in-person group interview (n = 5). Field- and genomics-based approaches agreed for 26 of 32 (81%) cases on likely location of TB acquisition. There was less agreement in the identification of specific source cases (13/22 or 59% of cases). Single-locus MIRU-VNTR variants and limited genetic diversity complicated the analysis. Qualitative data indicated that participants viewed genomic epidemiology as a useful tool to streamline investigations, particularly in differentiating latent TB reactivation from the recent transmission. Based on this, genomic data could be used to enhance CIs, focus resources, target interventions and aid in TB programme evaluation.

New insights into the mycobacterial PE and PPE proteins provide a framework for future research

The PE and PPE proteins of Mycobacterium tuberculosis have been studied with great interest since their discovery. Named after the conserved proline (P) and glutamic acid (E) residues in their N-terminal domains, these proteins are postulated to perform wide-ranging roles in virulence and immune modulation. However, technical challenges in studying these proteins and their encoding genes have hampered the elucidation of molecular mechanisms and leave many open questions regarding the biological functions mediated by these proteins. Here, I review the shared and unique characteristics of PE and PPE proteins from a molecular perspective linking this information to their functions in mycobacterial virulence. I discuss how the different subgroups (PE_PGRS, PPE-PPW, PPE-SVP and PPE-MPTR) are defined and why this classification of paramount importance to understand the PE and PPE proteins as individuals and or groups. The goal of this MicroReview is to summarize and structure the existing information on this gene family into a simplified framework of thinking about PE and PPE proteins and genes. Thereby, I hope to provide helpful starting points in studying these genes and proteins for researchers with different backgrounds. This has particular implications for the design and monitoring of novel vaccine candidates and in understanding the evolution of the M. tuberculosis complex.

Genotyping of Mycobacterium tuberculosis spreading in Hanoi, Vietnam using conventional and whole genome sequencing methods

Hanoi is the capital of Vietnam, one of the 30 countries with a high tuberculosis (TB) burden. Fundamental data on the molecular epidemiology of the disease is required for future TB management. To identify lineages and genotypes of Mycobacterium tuberculosis (Mtb), conventional genotyping data from clinical isolates of the Hanoi area was compared with whole genome sequencing (WGS) analysis from 332 of 470 samples. It was obtained from lineage-specific single nucleotide variants (SNVs), large sequence polymorphisms, spoligotyping, and variable number of tandem repeats (VNTR) analysis using mycobacterial interspersed repetitive unit (MIRU) and Japan anti-tuberculosis association (JATA) locus sets. This information was directly compared with results obtained from WGS. Mini-satellite repeat unit variants were identified using BLAST search against concatenated short read sequences, the RepUnitTyping tool. WGS analysis revealed that the Mtb strains tested are diverse and classified into lineage (L) 1, 2 and 4 (24.7, 57.2 and 18.1% respectively). The majority of the L2 strains were further divided into ancient and modern Beijing genotypes, and most of the L1 group were EAI4_VNM strains. Although conventional PCR-based genotyping results were mostly consistent with information obtained through WGS analysis, in-depth analysis identified aberrant deletions and spacers that may cause discordance. JATA-VNTR sets, including hypervariable loci, separated large Beijing genotypic clusters generated by MIRU15 into smaller groups. The distribution of repeat unit variants observed within 33 VNTR loci showed clear variation depending on the three lineages. WGS-based pairwise-SNV differences within VNTR-defined genotypic clusters were greater in L1 than in L2 and L4 (P = .001). Direct comparisons between results of PCR-based genotyping and in silico analysis of WGS data would bridge a gap between classical and modern technologies during this transition period, and provide further information on Mtb genotypes in specific geographical areas.

Type VII Secretion Substrates of Pathogenic Mycobacteria Are Processed by a Surface Protease

Tuberculosis, one of the world's most severe infectious diseases, is caused by Mycobacterium tuberculosis A major weapon of this pathogen is a unique cell wall that protects the pathogen from eradication by the immune system. Mycobacteria have specialized secretion systems, e.g., type VII secretion or ESX systems, to transport substrates across this cell wall. The largest group of proteins that are secreted by these ESX systems are the PE proteins. Previously, it was shown that the N-terminal PE domain of about 100 amino acids is required for secretion. Here, we describe the identification of an aspartic protease, designated PecA, that removes (part of) this PE domain at the cell surface. Nearly all of the observed PE_PGRS proteins are processed by PecA. Interestingly, the protease itself is also a secreted PE protein and subject to self-cleavage. Furthermore, a defect in surface processing has no effect on the activity of the PE lipase protein LipY but does seem to affect the functioning of other virulence factors, as a pecA mutant strain of Mycobacterium marinum shows moderate attenuation in zebrafish larvae. In conclusion, our results reveal the presence of a functional aspartic acid protease in M. marinum that cleaves LipY, itself as well as other members of the PE_PGRS family. Finally, mutants lacking PecA show growth attenuation in vivo, suggesting that PecA plays a role during infection.IMPORTANCE Aspartic proteases are common in eukaryotes and retroviruses but are relatively rare among bacteria (N. D. Rawlings and A. Bateman, BMC Genomics 10:437, 2009, https://doi.org/10.1186/1471-2164-10-437). In contrast to eukaryotic aspartic proteases, bacterial aspartic proteases are generally located in the cytoplasm. We have identified a surface-associated mycobacterial aspartic protease, PecA, which cleaves itself and many other type VII secretion substrates of the PE_PGRS family. PecA is present in most pathogenic mycobacterial species, including M. tuberculosis In addition, pathogenicity of M. marinum is reduced in the ΔpecA mutant, indicating that PecA contributes to virulence.

Whole genome sequencing, analyses of drug resistance-conferring mutations, and correlation with transmission of Mycobacterium tuberculosis carrying katG-S315T in Hanoi, Vietnam

Drug-resistant tuberculosis (TB) is a serious global problem, and pathogen factors involved in the transmission of isoniazid (INH)-resistant TB have not been fully investigated. We performed whole genome sequencing of 332 clinical Mycobacterium tuberculosis (Mtb) isolates collected from patients newly diagnosed with smear-positive pulmonary TB in Hanoi, Vietnam. Using a bacterial genome-wide approach based on linear mixed models, we investigated the associations between 31-bp k-mers and clustered strains harboring katG-S315T, a major INH-resistance mutation in the present cohort and in the second panel previously published in South Africa. Five statistically significant genes, namely, PPE18/19, gid, emrB, Rv1588c, and pncA, were shared by the two panels. We further identified variants of the genes responsible for these k-mers, which are relevant to the spread of INH-resistant strains. Phylogenetic convergence test showed that variants relevant to PPE46/47-like chimeric genes were significantly associated with the same phenotype in Hanoi. The associations were further confirmed after adjustment for the confounders. These findings suggest that genomic variations of the pathogen facilitate the expansion of INH-resistance TB, at least in part, and our study provides a new insight into the mechanisms by which drug-resistant Mtb maintains fitness and spreads in Asia and Africa.

Evidence of nitrite acting as a stable and robust inducer of non-cultivability in Mycobacterium tuberculosis with physiological relevance

Tuberculosis (TB) is the ninth leading cause of death worldwide, ranking above human immunodeficiency virus. Latency is the major obstacle in the eradication of this disease. How the physiology of the pathogen changes in transition to the latent stage needs to be understood. The latent bacteria extracted from animal hosts exist in a nonculturable (NC) phase, whereas bacteria extracted from most in vitro models are culture-positive. In the present study, we observed that nitrite, up to a concentration of 5 mM, shows the growth of Mycobacterium tuberculosis (MTB) in liquid media, but this effect starts reversing at higher concentrations. At a concentration of 10 mM, nitrite induces rapid nonculturability of MTB at the aerobic stage. This noncultivable dormancy was confirmed by analyzing the characteristics of NC bacteria. Further differential gene expression analyses clearly supported the formation of a dormancy phenotype. This study will be helpful for the use of this bacillus as a dormancy model in future studies on TB latency.

Approximate, simultaneous comparison of microbial genome architectures via syntenic anchoring of quiver representations

Motivation

A long-standing limitation in comparative genomic studies is the dependency on a reference genome, which hinders the spectrum of genetic diversity that can be identified across a population of organisms. This is especially true in the microbial world where genome architectures can significantly vary. There is therefore a need for computational methods that can simultaneously analyze the architectures of multiple genomes without introducing bias from a reference.

Results

In this article, we present Ptolemy: a novel method for studying the diversity of genome architectures-such as structural variation and pan-genomes-across a collection of microbial assemblies without the need of a reference. Ptolemy is a 'top-down' approach to compare whole genome assemblies. Genomes are represented as labeled multi-directed graphs-known as quivers-which are then merged into a single, canonical quiver by identifying 'gene anchors' via synteny analysis. The canonical quiver represents an approximate, structural alignment of all genomes in a given collection encoding structural variation across (sub-) populations within the collection. We highlight various applications of Ptolemy by analyzing structural variation and the pan-genomes of different datasets composing of Mycobacterium, Saccharomyces, Escherichia and Shigella species. Our results show that Ptolemy is flexible and can handle both conserved and highly dynamic genome architectures. Ptolemy is user-friendly-requires only FASTA-formatted assembly along with a corresponding GFF-formatted file-and resource-friendly-can align 24 genomes in ∼10 mins with four CPUs and <2 GB of RAM.

Availability and implementation

Github: https://github.com/AbeelLab/ptolemy.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification and Evaluation of Novel Protective Antigens for the Development of a Candidate Tuberculosis Subunit Vaccine

The development of a vaccine against tuberculosis (TB), a disease caused by Mycobacterium tuberculosis, is urgently needed. The only currently available vaccine, M. bovis BCG, has variable efficacy. One approach in the global vaccine development effort is focused on boosting BCG using subunit vaccines. The identification of novel antigens for inclusion in subunit vaccines is a critical step in the TB vaccine development pathway. We selected four novel mycobacterial antigens recognized during the course of human infection. A replication-deficient chimpanzee adenovirus (ChAdOx1) was constructed to express each antigen individually, and these vectors were evaluated for protective efficacy in murine M. tuberculosis challenge experiments. One antigen, PPE15 (Rv1039c), conferred significant and reproducible protection when administered alone and as a boost to BCG vaccination. We identified immunodominant epitopes to define the protective immune responses using tetramers and intravascular staining. Lung parenchymal CD4⁺ and CD8⁺ CXCR3⁺ KLRG1⁻ T cells, previously associated with protection against M. tuberculosis, were enriched in the vaccinated groups compared to the control groups. Further work to evaluate the protective efficacy of PPE15 in more stringent preclinical animal models, together with the identification of further novel protective antigens using this selection strategy, is now merited.

RD5-mediated lack of PE_PGRS and PPE-MPTR export in BCG vaccine strains results in strong reduction of antigenic repertoire but little impact on protection

Tuberculosis is the deadliest infectious disease worldwide. Although the BCG vaccine is widely used, it does not efficiently protect against pulmonary tuberculosis and an improved tuberculosis vaccine is therefore urgently needed. Mycobacterium tuberculosis uses different ESX/Type VII secretion (T7S) systems to transport proteins important for virulence and host immune responses. We recently reported that secretion of T7S substrates belonging to the mycobacteria-specific Pro-Glu (PE) and Pro-Pro-Glu (PPE) proteins of the PGRS (polymorphic GC-rich sequences) and MPTR (major polymorphic tandem repeat) subfamilies required both a functional ESX-5 system and a functional PPE38/71 protein for secretion. Inactivation of ppe38/71 and the resulting loss of PE_PGRS/PPE-MPTR secretion were linked to increased virulence of M. tuberculosis strains. Here, we show that a predicted total of 89 PE_PGRS/PPE-MPTR surface proteins are not exported by certain animal-adapted strains of the M. tuberculosis complex including M. bovis. This Δppe38/71-associated secretion defect therefore also occurs in the M. bovis-derived tuberculosis vaccine BCG and could be partially restored by introduction of the M. tuberculosis ppe38-locus. Epitope mapping of the PPE-MPTR protein PPE10, further allowed us to monitor T-cell responses in splenocytes from BCG/M. tuberculosis immunized mice, confirming the dependence of PPE10-specific immune-induction on ESX-5/PPE38-mediated secretion. Restoration of PE_PGRS/PPE-MPTR secretion in recombinant BCG neither altered global antigenic presentation or activation of innate immune cells, nor protective efficacy in two different mouse vaccination-infection models. This unexpected finding stimulates a reassessment of the immunomodulatory properties of PE_PGRS/PPE-MPTR proteins, some of which are contained in vaccine formulations currently in clinical evaluation.

One of the major findings of the pioneering Mycobacterium tuberculosis H37Rv genome sequencing project was the identification of the highly abundant PE and PPE proteins, named after their N-terminal motifs Pro–Glu (PE) or Pro–Pro–Glu (PPE). Within the 20 years of research since then, many claims were made that PE/PPE proteins, including the two large subgroups encoded by repetitive sequences with very high GC content (PE_PGRS and PPE-MPTR families), are exported to the bacterial surface or beyond, and show broad immunomodulatory impact on host-pathogen interaction. We thus screened strains from different branches of the M. tuberculosis complex, including the attenuated Mycobacterium bovis BCG vaccine strains, for their capacity to export PE_PGRS/PPE-MPTR proteins. Strikingly, we found that BCG strains were unable to export the plethora of PE_PGRS/PPE-MPTR proteins due to the absence of the region of difference RD5, which in M. tuberculosis encodes PPE38, required for PE_PGRS/PPE-MPTR export. Surprisingly, the restoration of PE_PGRS/PPE-MPTR export by genetic complementation in recombinant BCG did not result in immunomodulatory changes or altered protection in mouse models. Our results thus put into perspective the numerous reports on virulence-associated immunomodulatory impact of individual PE_PGRS and PPE-MPTR proteins and open novel questions on their biological function(s).

A complete high-quality MinION nanopore assembly of an extensively drug-resistant Mycobacterium tuberculosis Beijing lineage strain identifies novel variation in repetitive PE/PPE gene regions

A better understanding of the genomic changes that facilitate the emergence and spread of drug-resistant Mycobacterium tuberculosis strains is currently required. Here, we report the use of the MinION nanopore sequencer (Oxford Nanopore Technologies) to sequence and assemble an extensively drug-resistant (XDR) isolate, which is part of a modern Beijing sub-lineage strain, prevalent in Western Province, Papua New Guinea. Using 238-fold coverage obtained from a single flow-cell, de novo assembly of nanopore reads resulted into one contiguous assembly with 99.92 % assembly accuracy. Incorporation of complementary short read sequences (Illumina) as part of consensus error correction resulted in a 4 404 064 bp genome with 99.98 % assembly accuracy. This assembly had an average nucleotide identity of 99.7 % relative to the reference genome, H37Rv. We assembled nearly all GC-rich repetitive PE/PPE family genes (166/168) and identified variants within these genes. With an estimated genotypic error rate of 5.3 % from MinION data, we demonstrated identification of variants to include the conventional drug resistance mutations, and those that contribute to the resistance phenotype (efflux pumps/transporter) and virulence. Reference-based alignment of the assembly allowed detection of deletions and insertions. MinION sequencing provided a fully annotated assembly of a transmissible XDR strain from an endemic setting and showed its utility to provide further understanding of genomic processes within Mycobacterium tuberculosis.

Comparative Genomic Analysis of Two Clonally Related Multidrug Resistant Mycobacterium tuberculosis by Single Molecule Real Time Sequencing

Background: Multidrug-resistant tuberculosis (MDR-TB) is posing a major threat to global TB control. In this study, we focused on two consecutive MDR-TB isolated from the same patient before and after the initiation of anti-TB treatment. To better understand the genomic characteristics of MDR-TB, Single Molecule Real-Time (SMRT) Sequencing and comparative genomic analyses was performed to identify mutations that contributed to the stepwise development of drug resistance and growth fitness in MDR-TB under in vivo challenge of anti-TB drugs. Result: Both pre-treatment and post-treatment strain demonstrated concordant phenotypic and genotypic susceptibility profiles toward rifampicin, pyrazinamide, streptomycin, fluoroquinolones, aminoglycosides, cycloserine, ethionamide, and para-aminosalicylic acid. However, although both strains carried identical missense mutations at rpoB S531L, inhA C-15T, and embB M306V, MYCOTB Sensititre assay showed that the post-treatment strain had 16-, 8-, and 4-fold elevation in the minimum inhibitory concentrations (MICs) toward rifabutin, isoniazid, and ethambutol respectively. The results have indicated the presence of additional resistant-related mutations governing the stepwise development of MDR-TB. Further comparative genomic analyses have identified three additional polymorphisms between the clinical isolates. These include a single nucleotide deletion at nucleotide position 360 of rv0888 in pre-treatment strain, and a missense mutation at rv3303c (lpdA) V44I and a 6-bp inframe deletion at codon 67-68 in rv2071c (cobM) in the post-treatment strain. Multiple sequence alignment showed that these mutations were occurring at highly conserved regions among pathogenic mycobacteria. Using structural-based and sequence-based algorithms, we further predicted that the mutations potentially have deleterious effect on protein function. Conclusion: This is the first study that compared the full genomes of two clonally-related MDR-TB clinical isolates during the course of anti-TB treatment. Our work has demonstrated the robustness of SMRT Sequencing in identifying mutations among MDR-TB clinical isolates. Comparative genome analysis also suggested novel mutations at rv0888, lpdA, and cobM that might explain the difference in antibiotic resistance and growth pattern between the two MDR-TB strains.

Clinical implication of novel drug resistance-conferring mutations in resistant tuberculosis

Evolving novel and/or unfamiliar mutations are revolutionizing the pathways of antibiotic resistance of clinical tuberculosis. The accumulation and interaction of these poorly characterized mutations augment the complexity of resistant pathogenic strains and raise public health concerns. This article reviews our current understanding of the genetic changes that characterize drug resistance in tuberculosis and highlights the imperative for further investigations focusing on the effects of an individual mutation and interacting mutations with detailed strain epidemiology, particularly as these pertain to technology-limited countries with high tuberculosis incidence rates. Concomitantly, there is a need for the development, testing, and uptake of new tools for studying the effects of these mutations in drug resistance and fitness cost of the pathogen. Such genetic data are critical for effective localized and global tuberculosis control interventions and for accurate epidemiological predictions.

Impact of pe_pgrs33 Gene Polymorphisms on Mycobacterium tuberculosis Infection and Pathogenesis

PE_PGRS33 is a surface-exposed protein of Mycobacterium tuberculosis (Mtb) which exerts its role in macrophages entry and immunomodulation. In this study, we aimed to investigate the polymorphisms in the pe_pgrs33 gene of Mtb clinical isolates and evaluate their impact on protein functions. We sequenced pe_pgrs33 in a collection of 135 clinical strains, genotyped by 15-loci MIRU-VNTR and spoligotyping and belonging to the Mtb complex (MTBC). Overall, an association between pe_pgrs33 alleles and MTBC genotypes was observed and a dN/dS ratio of 0.64 was obtained, suggesting that a purifying selective pressure is acting on pe_pgrs33 against deleterious SNPs. Among a total of 19 pe_pgrs33 alleles identified in this study, 5 were cloned and used to complement the pe_pgrs33 knock-out mutant strain of Mtb H37Rv (MtbΔ33) to assess the functional impact of the respective polymorphisms in in vitro infections of primary macrophages. In human monocyte-derived macrophages (MDMs) infection, large in-frame and frameshift mutations were unable to restore the phenotype of Mtb H37Rv, impairing the cell entry capacity of Mtb, but neither its intracellular replication rate nor its immunomodulatory properties. In vivo studies performed in the murine model of tuberculosis (TB) demonstrated that the MtbΔ33 mutant strain was not impaired in the ability to infect and replicate in the lung tissue compared to the parental strain. Interestingly, MtbΔ33 showed an enhanced virulence during the chronic steps of infection compared to Mtb H37Rv. Similarly, the complementation of MtbΔ33 with a frameshift allele also resulted in a Mtb strain capable of causing a surprisingly enhanced tissue damage in murine lungs, during the chronic steps of infection. Together, these results further support the role of PE_PGRS33 in the pathogenesis and virulence of Mtb.

Patterns and processes of Mycobacterium bovis evolution revealed by phylogenomic analyses

Mycobacterium bovis is an important animal pathogen worldwide that parasitizes wild and domesticated vertebrate livestock as well as humans. A comparison of the five M. bovis complete genomes from the United Kingdom, South Korea, Brazil, and the United States revealed four novel large-scale structural variations of at least 2,000 bp. A comparative phylogenomic study including 2,483 core genes of 38 taxa from eight countries showed conflicting phylogenetic signal among sites. By minimizing this effect, we obtained a tree that better agrees with sampling locality. Results supported a relatively basal position of African strains (all isolated from Homo sapiens ), confirming that Africa was an important region for early diversification and that humans were one of the earliest hosts. Selection analyses revealed that functional categories such as “Lipid transport and metabolism,” “Cell cycle control, cell division, chromosome partitioning” and “Cell motility” were significant for the evolution of the group, besides other categories previously described, showing importance of genes associated with virulence and cholesterol metabolism in the evolution of M. bovis . PE/PPE genes, many of which are known to be associated with virulence, were major targets for large-scale polymorphisms, homologous recombination, and positive selection, evincing for the first time a plethora of evolutionary forces possibly contributing to differential adaptability in M. bovis . By assuming different priors, US strains originated and started to diversify around 150–5,210 ya. By further analyzing the largest set of US genomes to date (76 in total), obtained from 14 host species, we detected that hosts were not clustered in clades (except for a few cases), with some faster-evolving strains being detected, suggesting fast and ongoing reinfections across host species, and therefore, the possibility of new bovine tuberculosis outbreaks.

PPE Surface Proteins Are Required for Heme Utilization by Mycobacterium tuberculosis

Iron is essential for replication of Mycobacterium tuberculosis, but iron is efficiently sequestered in the human host during infection. Heme constitutes the largest iron reservoir in the human body and is utilized by many bacterial pathogens as an iron source. While heme acquisition is well studied in other bacterial pathogens, little is known in M. tuberculosis. To identify proteins involved in heme utilization by M. tuberculosis, a transposon mutant library was screened for resistance to the toxic heme analog gallium(III)-porphyrin (Ga-PIX). Inactivation of the ppe36, ppe62, and rv0265c genes resulted in resistance to Ga-PIX. Growth experiments using isogenic M. tuberculosis deletion mutants showed that PPE36 is essential for heme utilization by M. tuberculosis, while the functions of PPE62 and Rv0265c are partially redundant. None of the genes restored growth of the heterologous M. tuberculosis mutants, indicating that the proteins encoded by the genes have separate functions. PPE36, PPE62, and Rv0265c bind heme as shown by surface plasmon resonance spectroscopy and are associated with membranes. Both PPE36 and PPE62 proteins are cell surface accessible, while the Rv0265c protein is probably located in the periplasm. PPE36 and PPE62 are, to our knowledge, the first proline-proline-glutamate (PPE) proteins of M. tuberculosis that bind small molecules and are involved in nutrient acquisition. The absence of a virulence defect of the ppe36 deletion mutant indicates that the different iron acquisition pathways of M. tuberculosis may substitute for each other during growth and persistence in mice. The emerging model of heme utilization by M. tuberculosis as derived from this study is substantially different from those of other bacteria.

Tuberculosis is caused by Mycobacterium tuberculosis and is a devastating disease affecting eight million people each year. Iron is an essential nutrient for replication of M. tuberculosis in the human host. More than 70% of iron in the human body is bound in heme. Not surprisingly, many bacterial pathogens, including M. tuberculosis, are able to acquire iron from heme. However, the mechanism of heme uptake by M. tuberculosis is poorly understood. We have identified two novel surface proteins that bind heme and are required for heme utilization by M. tuberculosis. These findings constitute a major advancement of our understanding of iron acquisition by M. tuberculosis and show that M. tuberculosis has evolved heme uptake systems different from the paradigms established by other bacteria.

Comprehensive Essentiality Analysis of the Mycobacterium tuberculosis Genome via Saturating Transposon Mutagenesis

For decades, identifying the regions of a bacterial chromosome that are necessary for viability has relied on mapping integration sites in libraries of random transposon mutants to find loci that are unable to sustain insertion. To date, these studies have analyzed subsaturated libraries, necessitating the application of statistical methods to estimate the likelihood that a gap in transposon coverage is the result of biological selection and not the stochasticity of insertion. As a result, the essentiality of many genomic features, particularly small ones, could not be reliably assessed. We sought to overcome this limitation by creating a completely saturated transposon library in Mycobacterium tuberculosis In assessing the composition of this highly saturated library by deep sequencing, we discovered that a previously unknown sequence bias of the Himar1 element rendered approximately 9% of potential TA dinucleotide insertion sites less permissible for insertion. We used a hidden Markov model of essentiality that accounted for this unanticipated bias, allowing us to confidently evaluate the essentiality of features that contained as few as 2 TA sites, including open reading frames (ORF), experimentally identified noncoding RNAs, methylation sites, and promoters. In addition, several essential regions that did not correspond to known features were identified, suggesting uncharacterized functions that are necessary for growth. This work provides an authoritative catalog of essential regions of the M. tuberculosis genome and a statistical framework for applying saturating mutagenesis to other bacteria.

IMPORTANCE

Sequencing of transposon-insertion mutant libraries has become a widely used tool for probing the functions of genes under various conditions. The Himar1 transposon is generally believed to insert with equal probabilities at all TA dinucleotides, and therefore its absence in a mutant library is taken to indicate biological selection against the corresponding mutant. Through sequencing of a saturated Himar1 library, we found evidence that TA dinucleotides are not equally permissive for insertion. The insertion bias was observed in multiple prokaryotes and influences the statistical interpretation of transposon insertion (TnSeq) data and characterization of essential genomic regions. Using these insights, we analyzed a fully saturated TnSeq library for M. tuberculosis, enabling us to generate a comprehensive catalog of in vitro essentiality, including ORFs smaller than those found in any previous study, small (noncoding) RNAs (sRNAs), promoters, and other genomic features.

PE and PPE Genes: A Tale of Conservation and Diversity

PE and PPE are two large families of proteins typical of mycobacteria whose structural genes in the Mycobacterium tuberculosis complex (MTBC) occupy about 7% of the total genome. The most ancestral PE and PPE proteins are expressed by genes that belong to the same operon and in most cases are found inserted in the esx clusters, encoding a type VII secretion system. Duplication and expansion of pe and ppe genes, coupled with intragenomic and intergenomic recombination events, led to the emergence of the polymorphic pe_pgrs and ppe_mptr genes in the MTBC genome. The role and function of these proteins, and particularly of the polymorphic subfamilies, remains elusive, although it is widely accepted that PE and PPE proteins may represent a specialized collection used by MTBC to interact with the complex host immune system of mammals. In this chapter, we summarize what has been discovered since the identification of these genes in 1998, focusing on M. tuberculosis genetic variability, host-pathogen interaction and TB pathogenesis.

Type VII Secretion: A Highly Versatile Secretion System

Type VII secretion (T7S) systems of mycobacteria secrete substrates over the unusual diderm cell envelope. Furthermore, T7S gene clusters are present throughout the phylum Actinobacteria, and functional T7S-like systems have been identified in Firmicutes. Most of the T7S substrates can be divided into two families: the Esx proteins, which are found in both Firmicutes and Actinobacteria, and the PE and PPE proteins, which are more mycobacterium-specific. Members of both families have been shown to be secreted as folded heterodimers, suggesting that this is a conserved feature of T7S substrates. Most knowledge of the mechanism of T7S and the roles of T7S systems in virulence comes from studies of pathogenic mycobacteria. These bacteria can contain up to five T7S systems, called ESX-1 to ESX-5, each having its own role in bacterial physiology and virulence. In this article, we discuss the general composition of T7S systems and the role of the individual components in secretion. These conserved components include two membrane proteins with (predicted) enzymatic activities: a predicted ATPase (EccC), likely to be required for energy provision of T7S, and a subtilisin-like protease (MycP) involved in processing of specific substrates. Additionally, we describe the role of a conserved intracellular chaperone in T7S substrate recognition, based on recently published crystal structures and molecular analysis. Finally, we discuss system-specific features of the different T7S systems in mycobacteria and their role in pathogenesis and provide an overview of the role of T7S in virulence of other pathogenic bacteria.

Turkish and Japanese Mycobacterium tuberculosis sublineages share a remote common ancestor

Two geographically distant M. tuberculosis sublineages, Tur from Turkey and T3-Osaka from Japan, exhibit partially identical genotypic signatures (identical 12-loci MIRU-VNTR profiles, distinct spoligotyping patterns). We investigated T3-Osaka and Tur sublineages characteristics and potential genetic relatedness, first using MIRU-VNTR locus analysis on 21 and 25 samples of each sublineage respectively, and second comparing Whole Genome Sequences of 8 new samples to public data from 45 samples uncovering human tuberculosis diversity. We then tried to date their Most Recent Common Ancestor (MRCA) using three calibrations of SNP accumulation rate (long-term=0.03SNP/genome/year, derived from a tuberculosis ancestor of around 70,000years old; intermediate=0.2SNP/genome/year derived from a Peruvian mummy; short-term=0.5SNP/genome/year). To disentangle between these scenarios, we confronted the corresponding divergence times with major human history events and knowledge on human genetic divergence. We identified relatively high intrasublineage diversity for both T3-Osaka and Tur. We definitively proved their monophyly; the corresponding super-sublineage (referred to as "T3-Osa-Tur") shares a common ancestor with T3-Ethiopia and Ural sublineages but is only remotely related to other Euro-American sublineages such as X, LAM, Haarlem and S. The evolutionary scenario based on long-term evolution rate being valid until T3-Osa-Tur MRCA was not supported by Japanese fossil data. The evolutionary scenario relying on short-term evolution rate since T3-Osa-Tur MRCA was contradicted by human history and potential traces of past epidemics. T3-Osaka and Tur sublineages were found likely to have diverged between 800y and 2000years ago, potentially at the time of Mongol Empire. Altogether, this study definitively proves a strong genetic link between Turkish and Japanese tuberculosis. It provides a first hypothesis for calibrating TB Euro-American lineage molecular clock; additional studies are needed to reliably date events corresponding to intermediate depths in tuberculosis phylogeny.

Whole Genome Sequence Analysis of a Large Isoniazid-Resistant Tuberculosis Outbreak in London: A Retrospective Observational Study

BACKGROUND

A large isoniazid-resistant tuberculosis outbreak centred on London, United Kingdom, has been ongoing since 1995. The aim of this study was to investigate the power and value of whole genome sequencing (WGS) to resolve the transmission network compared to current molecular strain typing approaches, including analysis of intra-host diversity within a specimen, across body sites, and over time, with identification of genetic factors underlying the epidemiological success of this cluster.

METHODS AND FINDINGS

We sequenced 344 outbreak isolates from individual patients collected over 14 y (2 February 1998-22 June 2012). This demonstrated that 96 (27.9%) were indistinguishable, and only one differed from this major clone by more than five single nucleotide polymorphisms (SNPs). The maximum number of SNPs between any pair of isolates was nine SNPs, and the modal distance between isolates was two SNPs. WGS was able to reveal the direction of transmission of tuberculosis in 16 cases within the outbreak (4.7%), including within a multidrug-resistant cluster that carried a rare rpoB mutation associated with rifampicin resistance. Eleven longitudinal pairs of patient pulmonary isolates collected up to 48 mo apart differed from each other by between zero and four SNPs. Extrapulmonary dissemination resulted in acquisition of a SNP in two of five cases. WGS analysis of 27 individual colonies cultured from a single patient specimen revealed ten loci differed amongst them, with a maximum distance between any pair of six SNPs. A limitation of this study, as in previous studies, is that indels and SNPs in repetitive regions were not assessed due to the difficulty in reliably determining this variation.

CONCLUSIONS

Our study suggests that (1) certain paradigms need to be revised, such as the 12 SNP distance as the gold standard upper threshold to identify plausible transmissions; (2) WGS technology is helpful to rule out the possibility of direct transmission when isolates are separated by a substantial number of SNPs; (3) the concept of a transmission chain or network may not be useful in institutional or household settings; (4) the practice of isolating single colonies prior to sequencing is likely to lead to an overestimation of the number of SNPs between cases resulting from direct transmission; and (5) despite appreciable genomic diversity within a host, transmission of tuberculosis rarely results in minority variants becoming dominant. Thus, whilst WGS provided some increased resolution over variable number tandem repeat (VNTR)-based clustering, it was insufficient for inferring transmission in the majority of cases.

Identification of a Non-Pentapeptide Region Associated with Rapid Mycobacterial Evolution

A large portion of the coding capacity of Mycobacterium tuberculosis is devoted to the production of proteins containing several copies of the pentapeptide-2 repeat, namely the PE/PPE_MPTR proteins. Protein domain repeats have a variety of binding properties and are involved in protein-protein interactions as well as binding to other ligands such as DNA and RNA. They are not as common in prokaryotes, compared to eukaryotes, but the enrichment of pentapeptide-2 repeats in Mycobacteria constitutes an exception to that rule. The genes encoding the PE/PPE_MPTR proteins have undergone many rearrangements and here we have identified the expansion patterns across the Mycobacteria. We have performed a reclassification of the PE/PPE_MPTR proteins using cohesive regions rather than sparse domain architectures. It is clear that these proteins have undergone large insertions of several pentapeptide-2 domains appearing adjacent to one another in a repetitive pattern. Further, we have identified a non-pentapeptide motif associated with rapid mycobacterial evolution. The sequence composition of this region suggests a different structure compared to pentapeptide-2 repeats. By studying the evolution of the PE/PPE_MPTR proteins, we have distinguished features pertaining to tuberculosis-inducing species. Further studies of the non-pentapeptide region associated with repeat expansions promises to shed light on the pathogenicity of Mycobacterium tuberculosis.

Genetic signatures of Mycobacterium tuberculosis Nonthaburi genotype revealed by whole genome analysis of isolates from tuberculous meningitis patients in Thailand

Genome sequencing plays a key role in understanding the genetic diversity of Mycobacterium tuberculosis (M.tb). The genotype-specific character of M. tb contributes to tuberculosis severity and emergence of drug resistance. Strains of M. tb complex can be classified into seven lineages. The Nonthaburi (NB) genotype, belonging to the Indo-Oceanic lineage (lineage 1), has a unique spoligotype and IS6110-RFLP pattern but has not previously undergone a detailed whole genome analysis. In addition, there is not much information available on the whole genome analysis of M. tb isolates from tuberculous meningitis (TBM) patients in public databases. Isolates CSF3053, 46-5069 and 43-13838 of NB genotype were obtained from the cerebrospinal fluids of TBM Thai patients in Siriraj Hospital, Bangkok. The whole genomes were subjected to high throughput sequencing. The sequence data of each isolate were assembled into draft genome. The sequences were also aligned to reference genome, to determine genomic variations. Single nucleotide polymorphisms (SNPs) were obtained and grouped according to the functions of the genes containing them. They were compared with SNPs from 1,601 genomes, representing the seven lineages of M. tb complex, to determine the uniqueness of NB genotype. Susceptibility to first-line, second-line and other antituberculosis drugs were determined and related to the SNPs previously reported in drug-resistant related genes. The assembled genomes have an average size of 4,364,461 bp, 4,154 genes, 48 RNAs and 64 pseudogenes. A 500 base pairs deletion, which includes ppe50, was found in all isolates. RD239, specific for members of Indo Oceanic lineage, and RD147c were identified. A total of 2,202 SNPs were common to the isolates and used to classify the NB strains as members of sublineage 1.2.1. Compared with 1,601 genomes from the seven lineages of M. tb complex, mutation G2342203C was found novel to the isolates in this study. Three mutations (T28910C, C1180580T and C152178T) were found only in Thai NB isolates, including isolates from previous study. Although drug susceptibility tests indicated pan-susceptibility, non-synonymous SNPs previously reported to be associated with resistance to anti-tuberculous drugs; isoniazid, ethambutol, and ethionamide were identified in all the isolates. Non-synonymous SNPs were found in virulence genes such as the genes playing roles in apoptosis inhibition and phagosome arrest. We also report polymorphisms in essential genes, efflux pumps associated genes and genes with known epitopes. The analysis of the TBM isolates and the availability of the variations obtained will provide additional resources for global comparison of isolates from pulmonary tuberculosis and TBM. It will also contribute to the richness of genomic databases towards the prediction of antibiotic resistance, level of virulence and of origin of infection.

Recombination in pe/ppe genes contributes to genetic variation in Mycobacterium tuberculosis lineages

Background

Approximately 10 % of the Mycobacterium tuberculosis genome is made up of two families of genes that are poorly characterized due to their high GC content and highly repetitive nature. The PE and PPE families are typified by their highly conserved N-terminal domains that incorporate proline-glutamate (PE) and proline-proline-glutamate (PPE) signature motifs. They are hypothesised to be important virulence factors involved with host-pathogen interactions, but their high genetic variability and complexity of analysis means they are typically disregarded in genome studies.

Results

To elucidate the structure of these genes, 518 genomes from a diverse international collection of clinical isolates were de novo assembled. A further 21 reference M. tuberculosis complex genomes and long read sequence data were used to validate the approach. SNP analysis revealed that variation in the majority of the 168 pe/ppe genes studied was consistent with lineage. Several recombination hotspots were identified, notably pe_pgrs3 and pe_pgrs17. Evidence of positive selection was revealed in 65 pe/ppe genes, including epitopes potentially binding to major histocompatibility complex molecules.

Conclusions

This, the first comprehensive study of the pe and ppe genes, provides important insight into M. tuberculosis diversity and has significant implications for vaccine development.

Electronic supplementary material

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

Prediction of Certain Well-Characterized Domains of Known Functions within the PE and PPE Proteins of Mycobacteria

The PE and PPE protein family are unique to mycobacteria. Though the complete genome sequences for over 500 M. tuberculosis strains and mycobacterial species are available, few PE and PPE proteins have been structurally and functionally characterized. We have therefore used bioinformatics tools to characterize the structure and function of these proteins. We selected representative members of the PE and PPE protein family by phylogeny analysis and using structure-based sequence annotation identified ten well-characterized protein domains of known function. Some of these domains were observed to be common to all mycobacterial species and some were species specific.

Double-Strand DNA Break Repair in Mycobacteria

Discontinuity of both strands of the chromosome is a lethal event in all living organisms because it compromises chromosome replication. As such, a diversity of DNA repair systems has evolved to repair double-strand DNA breaks (DSBs). In part, this diversity of DSB repair systems has evolved to repair breaks that arise in diverse physiologic circumstances or sequence contexts, including cellular states of nonreplication or breaks that arise between repeats. Mycobacteria elaborate a set of three genetically distinct DNA repair pathways: homologous recombination, nonhomologous end joining, and single-strand annealing. As such, mycobacterial DSB repair diverges substantially from the standard model of prokaryotic DSB repair and represents an attractive new model system. In addition, the presence in mycobacteria of a DSB repair system that can repair DSBs in nonreplicating cells (nonhomologous end joining) or when DSBs arise between repeats (single-strand annealing) has clear potential relevance to Mycobacterium tuberculosis pathogenesis, although the exact role of these systems in M. tuberculosis pathogenesis is still being elucidated. In this article we will review the genetics of mycobacterial DSB repair systems, focusing on recent insights.

Mycobacterium tuberculosis Resists Stress by Regulating PE19 Expression

Mycobacterium tuberculosis requires the phosphate-sensing signal transduction system Pst/SenX3-RegX3 to resist host immune responses. A ΔpstA1 mutant lacking a Pst phosphate uptake system component is hypersensitive to diverse stress conditions in vitro and is attenuated in vivo due to constitutive expression of the phosphate starvation-responsive RegX3 regulon. Transcriptional profiling of the ΔpstA1 mutant revealed aberrant expression of certain pe and ppe genes. PE and PPE proteins, defined by conserved N-terminal domains containing Pro-Glu (PE) or Pro-Pro-Glu (PPE) motifs, account for a substantial fraction of the M. tuberculosis genome coding capacity, but their functions are largely uncharacterized. Because some PE and PPE proteins localize to the cell wall, we hypothesized that overexpression of these proteins sensitizes M. tuberculosis to stress by altering cell wall integrity. To test this idea, we deleted pe and ppe genes that were overexpressed by ΔpstA1 bacteria. Deletion of a single pe gene, pe19, suppressed hypersensitivity of the ΔpstA1 mutant to both detergent and reactive oxygen species. Ethidium bromide uptake assays revealed increased envelope permeability of the ΔpstA1 mutant that was dependent on PE19. The replication defect of the ΔpstA1 mutant in NOS2(-/-) mice was partially reversed by deletion of pe19, suggesting that increased membrane permeability due to PE19 overexpression sensitizes M. tuberculosis to host immunity. Our data indicate that PE19, which comprises only a 99-amino-acid PE domain, has a unique role in the permeability of the M. tuberculosis envelope that is regulated to resist stresses encountered in the host.

Co-evolution of Mycobacterium tuberculosis and Homo sapiens

The causative agent of human tuberculosis (TB), Mycobacterium tuberculosis, is an obligate pathogen that evolved to exclusively persist in human populations. For M. tuberculosis to transmit from person to person, it has to cause pulmonary disease. Therefore, M. tuberculosis virulence has likely been a significant determinant of the association between M. tuberculosis and humans. Indeed, the evolutionary success of some M. tuberculosis genotypes seems at least partially attributable to their increased virulence. The latter possibly evolved as a consequence of human demographic expansions. If co-evolution occurred, humans would have counteracted to minimize the deleterious effects of M. tuberculosis virulence. The fact that human resistance to infection has a strong genetic basis is a likely consequence of such a counter-response. The genetic architecture underlying human resistance to M. tuberculosis remains largely elusive. However, interactions between human genetic polymorphisms and M. tuberculosis genotypes have been reported. Such interactions are consistent with local adaptation and allow for a better understanding of protective immunity in TB. Future 'genome-to-genome' studies, in which locally associated human and M. tuberculosis genotypes are interrogated in conjunction, will help identify new protective antigens for the development of better TB vaccines.

Polymorphisms in the PE35 and PPE68 antigens in Mycobacterium tuberculosis strains may affect strain virulence and reflect ongoing immune evasion

Previous studies have demonstrated that the Pro‑Glu/Pro‑Pro‑Glu (PE/PPE) genes in strains of Mycobacterium tuberculosis exhibit high sequence variation and may be involved in antigenic variation and immune evasion. Region of Difference 1 (RD1), encoding genes from Rv3871 to Rv3879, was observed to be lost during the original derivation of Bacillus Calmette‑Guérin between 1908 and 1921. It has been previously demonstrated that two PE/PPE proteins, PE35 (Rv3872) and PPE68 (Rv3873), are encoded by RD1 and exhibit immunodominance. To explore the genetic diversity of PE35 and PPE68, and to evaluate the impact of sequence variation on the immune recognition of these proteins, 161 clinical M. tuberculosis strains were selected from China and comparative sequence analysis of PE35 and PPE68 was performed. The results indicated that polymorphisms in PE35 and PPE68 may lead to alterations in the function of these proteins, which may potentially affect strain virulence. In addition, the human T‑cell epitopes of PE35 and PPE68 were highly variable, suggesting that the two antigens may be involved in diversifying selection to evade host immunity. The prevalence of strains with PE35 mutations in the non‑Beijing family was significantly greater compared with the Beijing family, indicating that Beijing strains may be more conservative than non‑Beijing strains in this gene.

Sequence conservation analysis and in silico human leukocyte antigen-peptide binding predictions for the Mtb72F and M72 tuberculosis candidate vaccine antigens

Background

Requisites for an efficacious tuberculosis (TB) vaccine are a minimal genomic diversity among infectious Mycobacterium tuberculosis strains for the selected antigen, and the capability to induce robust T-cell responses in the majority of human populations. A tool in the identification of putative T-cell epitopes is in silico prediction of major histocompatibility complex (MHC)-peptide binding. Candidate TB vaccine antigen Mtb72F and its successor M72 are recombinant fusion proteins derived from Mtb32A and Mtb39A (encoded by Rv0125 and Rv1196, respectively). Adjuvanted Mtb72F and M72 candidate vaccines were shown to induce CD4⁺ T-cell responses in European, US, African and Asian populations.

Methods

Sequence conservation of Mtb32A, Mtb39A, Mtb72F and M72 among 46 strains (prevalent Mycobacterium strains causing human TB disease, and H37Ra) was assessed by multiple alignments using ClustalX. For Mtb32A, Mtb39A and Mtb72F, 15-mer human leukocyte antigen (HLA)-class II-binding peptides were predicted for 158 DRB1 alleles prevailing in populations with high TB burden, 6 DRB3/4/5, 8 DQ and 6 DP alleles, using NetMHCII-pan-3.0. Results for 3 DRB1 alleles were compared with previously published allele-matched in vitro binding data. Additional analyses were done for M72. Nonameric MHC class I-binding peptides in Mtb72F were predicted for three alleles representative of class I supertypes A02, A03 and B07, using seven prediction algorithms.

Results

Sequence identity among strains was ≥98 % for each protein. Residue changes in Mtb39A comprised primarily single residue or nucleotide insertions and/or deletions in repeat regions, and were observed in 67 % of strains. For Mtb72F, 156 DRB1, 6 DRB3/4/5, 7 DQ and 5 DP alleles were predicted to contain at least one MHC class II-binding peptide, and class I-binding peptides were predicted for each HLA-A/B allele. Comparison of predicted MHC-II-binding peptides with experimental data indicated that the algorithm’s sensitivity and specificity were variable among alleles.

Conclusions

The sequences from which Mtb72F and M72 are derived are highly conserved among representative Mycobacterium strains. Predicted putative T-cell epitopes in M72 and/or Mtb72F covered a wide array of HLA alleles. In silico binding predictions for class I- and II-binding putative epitopes can be complemented with biochemical verification of HLA binding capacity, processing and immunogenicity of the predicted peptides.

Electronic supplementary material

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

Evolution of Extensively Drug-Resistant Tuberculosis over Four Decades: Whole Genome Sequencing and Dating Analysis of Mycobacterium tuberculosis Isolates from KwaZulu-Natal

BACKGROUND

The continued advance of antibiotic resistance threatens the treatment and control of many infectious diseases. This is exemplified by the largest global outbreak of extensively drug-resistant (XDR) tuberculosis (TB) identified in Tugela Ferry, KwaZulu-Natal, South Africa, in 2005 that continues today. It is unclear whether the emergence of XDR-TB in KwaZulu-Natal was due to recent inadequacies in TB control in conjunction with HIV or other factors. Understanding the origins of drug resistance in this fatal outbreak of XDR will inform the control and prevention of drug-resistant TB in other settings. In this study, we used whole genome sequencing and dating analysis to determine if XDR-TB had emerged recently or had ancient antecedents.

METHODS AND FINDINGS

We performed whole genome sequencing and drug susceptibility testing on 337 clinical isolates of Mycobacterium tuberculosis collected in KwaZulu-Natal from 2008 to 2013, in addition to three historical isolates, collected from patients in the same province and including an isolate from the 2005 Tugela Ferry XDR outbreak, a multidrug-resistant (MDR) isolate from 1994, and a pansusceptible isolate from 1995. We utilized an array of whole genome comparative techniques to assess the relatedness among strains, to establish the order of acquisition of drug resistance mutations, including the timing of acquisitions leading to XDR-TB in the LAM4 spoligotype, and to calculate the number of independent evolutionary emergences of MDR and XDR. Our sequencing and analysis revealed a 50-member clone of XDR M. tuberculosis that was highly related to the Tugela Ferry XDR outbreak strain. We estimated that mutations conferring isoniazid and streptomycin resistance in this clone were acquired 50 y prior to the Tugela Ferry outbreak (katG S315T [isoniazid]; gidB 130 bp deletion [streptomycin]; 1957 [95% highest posterior density (HPD): 1937-1971]), with the subsequent emergence of MDR and XDR occurring 20 y (rpoB L452P [rifampicin]; pncA 1 bp insertion [pyrazinamide]; 1984 [95% HPD: 1974-1992]) and 10 y (rpoB D435G [rifampicin]; rrs 1400 [kanamycin]; gyrA A90V [ofloxacin]; 1995 [95% HPD: 1988-1999]) prior to the outbreak, respectively. We observed frequent de novo evolution of MDR and XDR, with 56 and nine independent evolutionary events, respectively. Isoniazid resistance evolved before rifampicin resistance 46 times, whereas rifampicin resistance evolved prior to isoniazid only twice. We identified additional putative compensatory mutations to rifampicin in this dataset. One major limitation of this study is that the conclusions with respect to ordering and timing of acquisition of mutations may not represent universal patterns of drug resistance emergence in other areas of the globe.

CONCLUSIONS

In the first whole genome-based analysis of the emergence of drug resistance among clinical isolates of M. tuberculosis, we show that the ancestral precursor of the LAM4 XDR outbreak strain in Tugela Ferry gained mutations to first-line drugs at the beginning of the antibiotic era. Subsequent accumulation of stepwise resistance mutations, occurring over decades and prior to the explosion of HIV in this region, yielded MDR and XDR, permitting the emergence of compensatory mutations. Our results suggest that drug-resistant strains circulating today reflect not only vulnerabilities of current TB control efforts but also those that date back 50 y. In drug-resistant TB, isoniazid resistance was overwhelmingly the initial resistance mutation to be acquired, which would not be detected by current rapid molecular diagnostics employed in South Africa that assess only rifampicin resistance.