Evolution and Diversity of Clonal Bacteria: The Paradigm of Mycobacterium tuberculosis

International audienceBACKGROUND: Mycobacterium tuberculosis complex species display relatively static genomes and 99.9% nucleotide sequence identity. Studying the evolutionary history of such monomorphic bacteria is a difficult and challenging task. PRINCIPAL FINDINGS: We found that single-nucleotide polymorphism (SNP) analysis of DNA repair, recombination and replication (3R) genes in a comprehensive selection of M. tuberculosis complex strains from across the world, yielded surprisingly high levels of polymorphisms as compared to house-keeping genes, making it possible to distinguish between 80% of clinical isolates analyzed in this study. Bioinformatics analysis suggests that a large number of these polymorphisms are potentially deleterious. Site frequency spectrum comparison of synonymous and non-synonymous variants and Ka/Ks ratio analysis suggest a general negative/purifying selection acting on these sets of genes that may lead to suboptimal 3R system activity. In turn, the relaxed fidelity of 3R genes may allow the occurrence of adaptive variants, some of which will survive. Furthermore, 3R-based phylogenetic trees are a new tool for distinguishing between M. tuberculosis complex strains. CONCLUSIONS/SIGNIFICANCE: This situation, and the consequent lack of fidelity in genome maintenance, may serve as a starting point for the evolution of antibiotic resistance, fitness for survival and pathogenicity, possibly conferring a selective advantage in certain stressful situations. These findings suggest that 3R genes may play an important role in the evolution of highly clonal bacteria, such as M. tuberculosis. They also facilitate further epidemiological studies of these bacteria, through the development of high-resolution tools. With many more microbial genomes being sequenced, our results open the door to 3R gene-based studies of adaptation and evolution of other, highly clonal bacteria

Elucidation of the Dual Role of Mycobacterial MoeZR in Molybdenum Cofactor Biosynthesis and Cysteine Biosynthesis

The pathway of molybdenum cofactor biosynthesis has been studied in detail by using proteins from Mycobacterium species, which contain several homologs associated with the first steps of Moco biosynthesis. While all Mycobacteria species contain a MoeZR, only some strains have acquired an additional homolog, MoeBR, by horizontal gene transfer. The role of MoeBR and MoeZR was studied in detail for the interaction with the two MoaD-homologs involved in Moco biosynthesis, MoaD1 and MoaD2, in addition to the CysO protein involved in cysteine biosynthesis. We show that both proteins have a role in Moco biosynthesis, while only MoeZR, but not MoeBR, has an additional role in cysteine biosynthesis. MoeZR and MoeBR were able to complement an E. coli moeB mutant strain, but only in conjunction with the Mycobacterial MoaD1 or MoaD2 proteins. Both proteins were able to sulfurate MoaD1 and MoaD2 in vivo, while only MoeZR additionally transferred the sulfur to CysO. Our in vivo studies show that Mycobacteria have acquired several homologs to maintain Moco biosynthesis. MoeZR has a dual role in Moco- and cysteine biosynthesis and is involved in the sulfuration of MoaD and CysO, whereas MoeBR only has a role in Moco biosynthesis, which is not an essential function for Mycobacteria

First Description of Natural and Experimental Conjugation between Mycobacteria Mediated by a Linear Plasmid

Background: in a previous study, we detected the presence of a Mycobacterium avium species-specific insertion sequence, IS1245, in Mycobacterium kansasii. Both species were isolated from a mixed M. avium-M. kansasii bone marrow culture from an HIV-positive patient. the transfer mechanism of this insertion sequence to M. kansasii was investigated here.Methodology/Principal Findings: A linear plasmid (pMA100) was identified in all colonies isolated from the M. avium-M. kansasii mixed culture carrying the IS1245 element. the linearity of pMA100 was confirmed. Other analyses suggested that pMA100 contained a covalently bound protein in the terminal regions, a characteristic of invertron linear replicons. Partial sequencing of pMA100 showed that it bears one intact copy of IS1245 inserted in a region rich in transposase-related sequences. These types of sequences have been described in other linear mycobacterial plasmids. Mating experiments were performed to confirm that pMA100 could be transferred in vitro from M. avium to M. kansasii. pMA100 was transferred by in vitro conjugation not only to the M. kansasii strain from the mixed culture, but also to two other unrelated M. kansasii clinical isolates, as well as to Mycobacterium bovis BCG Moreau.Conclusions/Significance: Horizontal gene transfer (HGT) is one of most important mechanisms leading to the evolution and diversity of bacteria. This work provides evidence for the first time on the natural occurrence of HGT between different species of mycobacteria. Gene transfer, mediated by a novel conjugative plasmid, was detected and experimentally reproduced.Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)Cooperacion Interuniversitaria UAM-Banco Santander con America Latina (CEAL), UAM, SpainConselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)Universidade Federal de São Paulo, Dept Microbiol Imunol & Parasitol, Escola Paulista Med, São Paulo, BrazilLab Nacl Comp Cient, Petropolis, BrazilUniv Autonoma Madrid, Fac Med, Dept Prevent Med, Madrid, SpainInst Adolfo Lutz Registro, Nucleo TB & Micobacterioses, São Paulo, BrazilUniversidade Federal de São Paulo, Dept Microbiol Imunol & Parasitol, Escola Paulista Med, São Paulo, BrazilFAPESP: FAPESP - 06/01533-9Web of Scienc

Simultaneous Analysis of Multiple Mycobacterium tuberculosis Knockdown Mutants In Vitro and In Vivo

Mycobacterium tuberculosis (Mtb) represents one of the most persistent bacterial threats to human health and new drugs are needed to limit its impact. Conditional knockdown mutants can help validate new drug targets, but the analysis of individual mutants is laborious and time consuming. Here, we describe quantitative DNA tags (qTags) and their use to simultaneously analyze conditional Mtb knockdown mutants that allowed silencing the glyoxylate and methylcitrate cycles (via depletion of isocitrate lyase, ICL), the serine protease Rv3671c, and the core subunits of the mycobacterial proteasome, PrcB and PrcA. The impact of gene silencing in multi-strain cultures was determined by measuring the relative abundance of mutant-specific qTags with real-time PCR. This achieved accurate quantification over a broad range of qTag abundances and depletion of ICL, Rv3671c, or PrcBA resulted in the expected impairment of growth of Mtb with butyrate as the primary carbon source, survival during oxidative stress, acid stress and starvation. The impact of depleting ICL, Rv3671c, or PrcBA in multi-strain mouse infections was analyzed with two approaches. We first measured the relative abundance of mutant-specific qTags in total chromosomal DNA isolated from bacteria that were recovered from infected lungs on agar plates. We then developed a two-step amplification procedure, which allowed us to measure the abundances of individual mutants directly in infected lung tissue. Both strategies confirmed that inactivation of Rv3671c and PrcBA severely reduced persistence of Mtb in mice. The multi-strain infections furthermore suggested that silencing ICL not only prevented growth of Mtb during acute infections but also prevented survival of Mtb during chronic infections. Analyses of the ICL knockdown mutant in single-strain infections confirmed this and demonstrated that silencing of ICL during chronic infections impaired persistence of Mtb to the extent that the pathogen was cleared from the lungs of most mice

Uncovering the Prevalence and Diversity of Integrating Conjugative Elements in Actinobacteria

Horizontal gene transfer greatly facilitates rapid genetic adaptation of bacteria to shifts in environmental conditions and colonization of new niches by allowing one-step acquisition of novel functions. Conjugation is a major mechanism of horizontal gene transfer mediated by conjugative plasmids and integrating conjugative elements (ICEs). While in most bacterial conjugative systems DNA translocation requires the assembly of a complex type IV secretion system (T4SS), in Actinobacteria a single DNA FtsK/SpoIIIE-like translocation protein is required. To date, the role and diversity of ICEs in Actinobacteria have received little attention. Putative ICEs were searched for in 275 genomes of Actinobacteria using HMM-profiles of proteins involved in ICE maintenance and transfer. These exhaustive analyses revealed 144 putative FtsK/SpoIIIE-type ICEs and 17 putative T4SS-type ICEs. Grouping of the ICEs based on the phylogenetic analyses of maintenance and transfer proteins revealed extensive exchanges between different sub-families of ICEs. 17 ICEs were found in Actinobacteria from the genus Frankia, globally important nitrogen-fixing microorganisms that establish root nodule symbioses with actinorhizal plants. Structural analysis of ICEs from Frankia revealed their unexpected diversity and a vast array of predicted adaptive functions. Frankia ICEs were found to excise by site-specific recombination from their host's chromosome in vitro and in planta suggesting that they are functional mobile elements whether Frankiae live as soil saprophytes or plant endosymbionts. Phylogenetic analyses of proteins involved in ICEs maintenance and transfer suggests that active exchange between ICEs cargo-borne and chromosomal genes took place within the Actinomycetales order. Functionality of Frankia ICEs in vitro as well as in planta lets us anticipate that conjugation and ICEs could allow the development of genetic manipulation tools for this challenging microorganism and for many other Actinobacteria

High Content Phenotypic Cell-Based Visual Screen Identifies Mycobacterium tuberculosis Acyltrehalose-Containing Glycolipids Involved in Phagosome Remodeling

The ability of the tubercle bacillus to arrest phagosome maturation is considered one major mechanism that allows its survival within host macrophages. To identify mycobacterial genes involved in this process, we developed a high throughput phenotypic cell-based assay enabling individual sub-cellular analysis of over 11,000 Mycobacterium tuberculosis mutants. This very stringent assay makes use of fluorescent staining for intracellular acidic compartments, and automated confocal microscopy to quantitatively determine the intracellular localization of M. tuberculosis. We characterised the ten mutants that traffic most frequently into acidified compartments early after phagocytosis, suggesting that they had lost their ability to arrest phagosomal maturation. Molecular analysis of these mutants revealed mainly disruptions in genes involved in cell envelope biogenesis (fadD28), the ESX-1 secretion system (espL/Rv3880), molybdopterin biosynthesis (moaC1 and moaD1), as well as in genes from a novel locus, Rv1503c-Rv1506c. Most interestingly, the mutants in Rv1503c and Rv1506c were perturbed in the biosynthesis of acyltrehalose-containing glycolipids. Our results suggest that such glycolipids indeed play a critical role in the early intracellular fate of the tubercle bacillus. The unbiased approach developed here can be easily adapted for functional genomics study of intracellular pathogens, together with focused discovery of new anti-microbials

The use of genomic signature distance between bacteriophages and their hosts displays evolutionary relationships and phage growth cycle determination

<p>Abstract</p> <p>Background</p> <p>Bacteriophage classification is mainly based on morphological traits and genome characteristics combined with host information and in some cases on phage growth lifestyle. A lack of molecular tools can impede more precise studies on phylogenetic relationships or even a taxonomic classification. The use of methods to analyze genome sequences without the requirement for homology has allowed advances in classification.</p> <p>Results</p> <p>Here, we proposed to use genome sequence signature to characterize bacteriophages and to compare them to their host genome signature in order to obtain host-phage relationships and information on their lifestyle. We analyze the host-phage relationships in the four most representative groups of Caudoviridae, the dsDNA group of phages. We demonstrate that the use of phage genomic signature and its comparison with that of the host allows a grouping of phages and is also able to predict the host-phage relationships (lytic <it>vs</it>. temperate).</p> <p>Conclusions</p> <p>We can thus condense, in relatively simple figures, this phage information dispersed over many publications.</p

Non Mycobacterial Virulence Genes in the Genome of the Emerging Pathogen Mycobacterium abscessus

Mycobacterium abscessus is an emerging rapidly growing mycobacterium (RGM) causing a pseudotuberculous lung disease to which patients with cystic fibrosis (CF) are particularly susceptible. We report here its complete genome sequence. The genome of M. abscessus (CIP 104536T) consists of a 5,067,172-bp circular chromosome including 4920 predicted coding sequences (CDS), an 81-kb full-length prophage and 5 IS elements, and a 23-kb mercury resistance plasmid almost identical to pMM23 from Mycobacterium marinum. The chromosome encodes many virulence proteins and virulence protein families absent or present in only small numbers in the model RGM species Mycobacterium smegmatis. Many of these proteins are encoded by genes belonging to a “mycobacterial” gene pool (e.g. PE and PPE proteins, MCE and YrbE proteins, lipoprotein LpqH precursors). However, many others (e.g. phospholipase C, MgtC, MsrA, ABC Fe(3+) transporter) appear to have been horizontally acquired from distantly related environmental bacteria with a high G+C content, mostly actinobacteria (e.g. Rhodococcus sp., Streptomyces sp.) and pseudomonads. We also identified several metabolic regions acquired from actinobacteria and pseudomonads (relating to phenazine biosynthesis, homogentisate catabolism, phenylacetic acid degradation, DNA degradation) not present in the M. smegmatis genome. Many of the “non mycobacterial” factors detected in M. abscessus are also present in two of the pathogens most frequently isolated from CF patients, Pseudomonas aeruginosa and Burkholderia cepacia. This study elucidates the genetic basis of the unique pathogenicity of M. abscessus among RGM, and raises the question of similar mechanisms of pathogenicity shared by unrelated organisms in CF patients

Whole genome identification of Mycobacterium tuberculosis vaccine candidates by comprehensive data mining and bioinformatic analyses

<p>Abstract</p> <p>Background</p> <p><it>Mycobacterium tuberculosis</it>, the causative agent of tuberculosis (TB), infects ~8 million annually culminating in ~2 million deaths. Moreover, about one third of the population is latently infected, 10% of which develop disease during lifetime. Current approved prophylactic TB vaccines (BCG and derivatives thereof) are of variable efficiency in adult protection against pulmonary TB (0%–80%), and directed essentially against early phase infection.</p> <p>Methods</p> <p>A genome-scale dataset was constructed by analyzing published data of: (1) global gene expression studies under conditions which simulate intra-macrophage stress, dormancy, persistence and/or reactivation; (2) cellular and humoral immunity, and vaccine potential. This information was compiled along with revised annotation/bioinformatic characterization of selected gene products and <it>in silico </it>mapping of T-cell epitopes. Protocols for scoring, ranking and prioritization of the antigens were developed and applied.</p> <p>Results</p> <p>Cross-matching of literature and <it>in silico</it>-derived data, in conjunction with the prioritization scheme and biological rationale, allowed for selection of 189 putative vaccine candidates from the entire genome. Within the 189 set, the relative distribution of antigens in 3 functional categories differs significantly from their distribution in the whole genome, with reduction in the Conserved hypothetical category (due to improved annotation) and enrichment in Lipid and in Virulence categories. Other prominent representatives in the 189 set are the PE/PPE proteins; iron sequestration, nitroreductases and proteases, all within the Intermediary metabolism and respiration category; ESX secretion systems, resuscitation promoting factors and lipoproteins, all within the Cell wall category. Application of a ranking scheme based on qualitative and quantitative scores, resulted in a list of 45 best-scoring antigens, of which: 74% belong to the dormancy/reactivation/resuscitation classes; 30% belong to the Cell wall category; 13% are classical vaccine candidates; 9% are categorized Conserved hypotheticals, all potentially very potent T-cell antigens.</p> <p>Conclusion</p> <p>The comprehensive literature and <it>in silico</it>-based analyses allowed for the selection of a repertoire of 189 vaccine candidates, out of the whole-genome 3989 ORF products. This repertoire, which was ranked to generate a list of 45 top-hits antigens, is a platform for selection of genes covering all stages of <it>M. tuberculosis </it>infection, to be incorporated in rBCG or subunit-based vaccines.</p