Comparative genomics of Esx genes from clinical isolates of Mycobacterium tuberculosis provides evidence for gene conversion and epitope variation

The 23-membered Esx protein family is involved in the host-pathogen interactions of Mycobacterium tuberculosis. These secreted proteins are among the most immunodominant antigens recognized by the human immune system and have thus been used to develop vaccines and immunodiagnostic tests for tuberculosis (TB). Gene pairs for 10 Esx proteins are contained in the ESX-1 to ESX-5 loci, encoding type VII secretion systems. A subset of Esx proteins can be further classified into the Mtb9.9, QILSS, and TB10.4 subfamilies. To survey genetic diversity in the Esx family and its potential for antigenic variation, we sequenced all esx genes from 108 clinical isolates of M. tuberculosis from different clades by using a targeted approach. A total of 109 unique single nucleotide polymorphisms (SNPs) were observed, and 59 of these were nonsynonymous. Some of the resultant amino acid substitutions affect known Esx epitopes, including two in the EsxB (CFP-10) and EsxH (TB10.4) antigens. Assessment of the SNP distribution across the Esx proteins revealed high genetic variability, especially in the Mtb9.9 and QILSS subfamilies, and more conservation in the ESX-1 to ESX-4 loci. Comparison of the DNA sequences of variable esx genes provided clear evidence for recombination events between different genes in the same strain, some of which are predicted to truncate the corresponding protein. Many of these polymorphisms escape detection by ultrahigh-throughput sequencing using short sequence reads, as such approaches cannot distinguish between closely related genes. The esx gene family is dynamic, and sequence changes likely lead to immune variation.

The PhoP-Dependent ncRNA Mcr7 Modulates the TAT Secretion System in Mycobacterium tuberculosis

Andrej Benjak, Stewart Cole, Jacques Rougemont, Claudia Sala, Swapna Uplekar

The PhoPR two-component system is essential for virulence in Mycobacterium tuberculosis where it controls expression of approximately 2% of the genes, including those for the ESX-1 secretion apparatus, a major virulence determinant. Mutations in phoP lead to compromised production of pathogen-specific cell wall components and attenuation both ex vivo and in vivo. Using antibodies against the native protein in ChIP-seq experiments (chromatin immunoprecipitation followed by high-throughput sequencing) we demonstrated that PhoP binds to at least 35 loci on the M. tuberculosis genome. The PhoP regulon comprises several transcriptional regulators as well as genes for polyketide synthases and PE/PPE proteins. Integration of ChIP-seq results with high-resolution transcriptomic analysis (RNA-seq) revealed that PhoP controls 30 genes directly, whilst regulatory cascades are responsible for signal amplification and downstream effects through proteins like EspR, which controls Esx1 function, via regulation of the espACD operon. The most prominent site of PhoP regulation was located in the intergenic region between rv2395 and PE_PGRS41, where the mcr7 gene codes for a small non-coding RNA (ncRNA). Northern blot experiments confirmed the absence of Mcr7 in an M. tuberculosis phoP mutant as well as low-level expression of the ncRNA in M. tuberculosis complex members other than M. tuberculosis. By means of genetic and proteomic analyses we demonstrated that Mcr7 modulates translation of the tatC mRNA thereby impacting the activity of the Twin Arginine Translocation (Tat) protein secretion apparatus. As a result, secretion of the immunodominant Ag85 complex and the beta-lactamase BlaC is affected, among others. Mcr7, the first ncRNA of M. tuberculosis whose function has been established, therefore represents a missing link between the PhoPR two-component system and the downstream functions necessary for successful infection of the host.

Public Library of Science2014

Functional Characterization of Nucleoid Associated Proteins Acting as Global Transcription Factors in Mycobacterium tuberculosis

Nina Theres Odermatt

The fatal lung disease tuberculosis is caused by the airborne Mycobacterium tuberculosis, a versatile pathogen adapted to rapidly changing environments. Instead of being eradicated by phagocytic cells of its human host, bacilli tune macrophages to support their own growth and even mask their presence from the immune system for several decades. Rapid adjustment of gene expression is critical for bacterial survival and heavily relies on nucleoid-associated proteins (NAPs). NAPs contribute to active DNA management by altering the chromosomal topology through bending, bridging and looping the DNA. These conformational changes can bring distant genetic loci into close spatial proximity or influence DNA supercoiling and therefore accessibility of the transcription machinery. Apart from their architectural role, NAPs moreover act as global transcription factors by direct regulation of numerous genes. In M. tuberculosis, five proteins were assigned a role as NAP. Among these, EspR, HupB and Lsr2 are crucial not only for virulence but also for cellular metabolism. This thesis focuses on the NAPs mIHF and H-NS with the objective of determining their function and target regulon. Both proteins were investigated by means of genetic manipulation, phenotype assessment and structural studies, and the efficiency of a potential new anti-tuberculosis drug acting on Lsr2 was assessed. Chrysomycin, described as specific inhibitor of Lsr2-DNA complex formation, was found to intercalate into the DNA. The resulting toxic effect on both its target M. tuberculosis as well as on eukaryotic cells rendered further development of the compound as an anti-tuberculosis drug futile. We demonstrated that Rv3852, formerly annotated as H-NS, does not act as a NAP. Deletion of the rv3852 gene had no effect on the in vitro phenotype of M. tuberculosis, did not alter nucleoid spread nor position and had no influence on virulence in mice. The mIHF protein on the other hand is not only essential for active bacterial growth, but also indispensable for survival. Generation of a conditional knockdown mutant showed that depletion of mIHF led to elongated cells devoid of septa with abnormal DNA localization and finally to cell death. The target regulon of mIHF was thoroughly studied by mapping its binding sites on the bacterial genome and by identifying genes that were differentially expressed upon depletion of the protein. We found that mIHF has a strong effect on virulence gene expression and, similar to EspR, possesses a major binding site upstream of one of the main virulence factor operons espACD. Analysis of the transcriptional response revealed that mIHF is further involved in the bacterial response to the hostâs immune system, including control of nutrient pathways as well as global protein and nucleic acid synthesis. To define how mIHF interacts with DNA and influences its 3D organisation, the protein structure of mIHF was determined by nuclear magnetic resonance spectroscopy. Binding of mIHF introduced left-hand loops into linear as well as supercoiled DNA substrates, therefore unwinding condensed DNA. We identified two DNA binding domains in mIHF and showed that its stability increased substantially upon DNA binding. All together, the findings of this thesis contribute to a better understanding of the complex gene regulatory network of M. tuberculosis, advancing the knowledge necessary to eventually defeat tuberculosis, a disease that has plagued humanity for millennia.

EPFL2017

Genomics and transcriptomics of the Mycobacterium tuberculosis complex

Swapna Uplekar

The goal of eliminating tuberculosis (TB) by 2050 depends on the development of improved TB diagnostics, drugs and vaccines. Advances in these areas require a deep understanding of the disease and its causative agent, Mycobacterium tuberculosis (M. tb). Mycobacterial species that cause TB in humans and other mammalian hosts are grouped within the M. tb complex. Development of powerful technologies such as next-generation sequencing and microarrays opened up new avenues for comparative and functional genomics of the M. tb complex. Due to the large and increasingly complex datasets generated from these technologies, the bottleneck in biological investigation has shifted from data generation to analysis. The objectives of this thesis were to establish and employ strategies for the analysis, integration, and interpretation of high-throughput sequencing and microarray datasets using a range of bioinformatics and statistical tools. In the area of comparative genomics, we assessed the genetic diversity in the M. tb complex using various methods, such as SNP (single nucleotide polymorphism) genotyping, automated Sanger sequencing and next-generation sequencing. In a study comparing the genomes of the virulent M. bovis and M. bovis BCG vaccine strains, we identified a set of SNPs that were common to all BCG strains, and could provide novel insights on the molecular basis of BCG attenuation. In another study, we surveyed the genetic variation in the highly immunodominant esx gene family among clinical isolates of M. tb and identified sequence polymorphisms in known T- cell epitopes on Esx proteins that could affect their immunogenicity. We exploited the power of next-generation sequencing to detect sequence variation among M. tb strains that could result in phenotypic differences. By comparing the genomes of drug-resistant mutants with the sensitive wild-type strain we were able to identify the target of the anti-TB drug, pyridomycin. Using a similar approach we identified a mutation that makes M. tb strains incapable of producing PDIMs (phthiocerol dimycocerosates), which are cell wall associated lipids involved in M. tb virulence. In the area of functional genomics, we mapped genome-wide binding sites for transcription factors using chromatin immunoprecipitation followed by hybridization to microarrays (ChIP-on-chip) or sequencing (ChIP-seq), and performed transcription profiling by means of high-throughput cDNA sequencing (RNA-seq). We carried out a comprehensive study to characterize the whole transcriptome of M. tb in exponential and stationary phases of growth, and understand the genome-wide dynamics of two key components of the transcription machinery, namely, RNA polymerase and NusA. By systematic integration of the ChIP-seq and RNA-seq data, we identified a set of transcription units (TU) in the M. tb genome, and mapped their putative promoters. Analysis of RNAP and NusA binding across the promoter and body of TUs and their correlation with transcription uncovered new functional aspects of the transcriptional complex in M. tb. We also exploited the ChIP-on-chip and ChIP-seq technologies to define the regulon of the M. tb sigma factor F, and gain a better understanding of the regulatory role of the nucleoid associated protein, EspR. Altogether, this thesis has improved our knowledge of the evolution, physiology and virulence of the M. tb complex. In addition, we have established next generation sequencing as a powerful tool for comparative and functional studies, with potential applications in the clinical setting.

EPFL2012