|Genomic analysis of smooth tubercle bacilli provides insights into ancestry and pathoadaptation of Mycobacterium tuberculosis.|
|Jump to Full Text|
|PMID: 23291586 Owner: NLM Status: Publisher|
|Global spread and limited genetic variation are hallmarks of M. tuberculosis, the agent of human tuberculosis. In contrast, Mycobacterium canettii and related tubercle bacilli that also cause human tuberculosis and exhibit unusual smooth colony morphology are restricted to East Africa. Here, we sequenced and analyzed the whole genomes of five representative strains of smooth tubercle bacilli (STB) using Sanger (4-5× coverage), 454/Roche (13-18× coverage) and/or Illumina DNA sequencing (45-105× coverage). We show that STB isolates are highly recombinogenic and evolutionarily early branching, with larger genome sizes, higher rates of genetic variation, fewer molecular scars and distinct CRISPR-Cas systems relative to M. tuberculosis. Despite the differences, all tuberculosis-causing mycobacteria share a highly conserved core genome. Mouse infection experiments showed that STB strains are less persistent and virulent than M. tuberculosis. We conclude that M. tuberculosis emerged from an ancestral STB-like pool of mycobacteria by gain of persistence and virulence mechanisms, and we provide insights into the molecular events involved.|
|Philip Supply; Michael Marceau; Sophie Mangenot; David Roche; Carine Rouanet; Varun Khanna; Laleh Majlessi; Alexis Criscuolo; Julien Tap; Alexandre Pawlik; Laurence Fiette; Mickael Orgeur; Michel Fabre; Cécile Parmentier; Wafa Frigui; Roxane Simeone; Eva C Boritsch; Anne-Sophie Debrie; Eve Willery; Danielle Walker; Michael A Quail; Laurence Ma; Christiane Bouchier; Grégory Salvignol; Fadel Sayes; Alessandro Cascioferro; Torsten Seemann; Valérie Barbe; Camille Locht; Maria-Cristina Gutierrez; Claude Leclerc; Stephen D Bentley; Timothy P Stinear; Sylvain Brisse; Claudine Médigue; Julian Parkhill; Stéphane Cruveiller; Roland Brosch|
Related Documents :
|23324426 - Mycobacterium tuberculosis beijing outbreak in a school in marseille, france, 2012.
22924146 - Tuberculosis of the parotid gland.
23176186 - Complex sputum microbial composition in patients with pulmonary tuberculosis.
22679206 - Associations between awareness of beyondblue and mental health literacy in australian y...
19408076 - Surgical management of colorectal cancer in patients with psychiatric disorders.
22453166 - Being men with eating disorders: perspectives of male eating disorder service-users.
|Type: JOURNAL ARTICLE Date: 2013-1-06|
|Title: Nature genetics Volume: - ISSN: 1546-1718 ISO Abbreviation: Nat. Genet. Publication Date: 2013 Jan|
|Created Date: 2013-1-7 Completed Date: - Revised Date: -|
Medline Journal Info:
|Nlm Unique ID: 9216904 Medline TA: Nat Genet Country: -|
|Languages: ENG Pagination: - Citation Subset: -|
|1] Institut National de la Santé et de la Recherche Médicale (INSERM) U1019, Center for Infection and Immunity of Lille, Lille, France.  Centre National de la Recherche Scientifique (CNRS), Unité Mixte de Recherche (UMR) 8204, Center for Infection and Immunity of Lille, Lille, France.  Université Lille Nord de France, Center for Infection and Immunity of Lille, Lille, France.  Institut Pasteur de Lille, Center for Infection and Immunity of Lille, Lille, France.|
|APA/MLA Format Download EndNote Download BibTex|
Journal ID (nlm-journal-id): 9216904
Journal ID (pubmed-jr-id): 2419
Journal ID (nlm-ta): Nat Genet
Journal ID (iso-abbrev): Nat. Genet.
nihms-submitted publication date: Day: 30 Month: 11 Year: 2013
Electronic publication date: Day: 06 Month: 1 Year: 2013
Print publication date: Month: 2 Year: 2013
pmc-release publication date: Day: 09 Month: 12 Year: 2013
Volume: 45 Issue: 2
E-location ID: 10.1038/ng.2517
PubMed Id: 23291586
|Genome analysis of smooth tubercle bacilli provides insights into ancestry and pathoadaptation of the etiologic agent of tuberculosis|
|Eva C. Boritsch7|
|Michael A. Quail14|
|Timothy P. Stinear18|
1Institut National de la Santé et de la Recherche Médicale (INSERM), U1019, Center for Infection and Immunity of Lille, Lille, France
2Centre National de la Recherche Scientifique (CNRS), Unite mixte de recherche (UMR) 8204, Center for Infection and Immunity of Lille, Lille, France
3Univ Lille Nord de France, Center for Infection and Immunity of Lille, Lille, France
4Institut Pasteur de Lille, Center for Infection and Immunity of Lille, Lille, France
5CNRS-UMR 8030 , Evry, France
6Commissariat à l’Energie Atomique et aux Energies Alternatives CEA/DSV/IG/Genoscope, LABGeM, Evry, France
7Institut Pasteur, Unit for Integrated Mycobacterial Pathogenomics, Paris, France
8Institut Pasteur, Unité de Régulation Immunitaire et Vaccinologie, Paris, France
9INSERM U1041, Paris, France
10Institut Pasteur, Genotyping of Pathogens and Public Health (PF8), Paris, France
11Institut Pasteur, Unité d’Histopathologie Humaine et Modèles Animaux, Paris, France
12Université Versailles-Saint Quentin en Yvelines, Faculté de Médecine, DER Histologie, Versailles, France
13Laboratoire de Biologie Clinique, HIA Percy, Clamart, France
14Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
15Institut Pasteur, Genopole, Platform Genomics PF1, Paris, France
16Victorian Bioinformatics Consortium, Monash University, Clayton, Australia
17Institut Pasteur, Department d’Infection et d’Epidemiologie, Paris, France
18Department of Microbiology and Immunology, University of Melbourne, Parkville, Australia
|*To whom correspondence should be addressed: email@example.com and firstname.lastname@example.org
Author contributions P.S., M.-C.G. and R.B. designed the study; P.S. and R.B. analyzed data and wrote the paper with comments from all authors; M.M., D.W., J.P., C.M. and S.Be. annotated the genomes; S.M., E.C.B., M.A.Q. L.Ma, C.B. and V.B. performed and/or verified the finishing and assembly of sequences; D.R. and S.C. performed SNP analyses and database management with support of G.S.; C.R., A.P., W.F., R.S., A.-S.D., E.W., A.Ca. and M.-C.G. performed mouse infection experiments, in vivo data analyses and/or mycobacterial growth assays; L.Maj., F.S., C.Lo. and C.Le. conducted and/or analyzed immune assays; J.T., A.Cr. and S.Br. conducted MLST, recombination and/or phylogenetic analyses; L.F. conducted histopathological analyses; V.K, M.O. and C.P. created bioinformatics tools and analyzed data; M.F. isolated STB strains; T.S. and T.P.S conducted core genome and NeighborNet analyses.
Mycobacterium tuberculosis is a pervasive human pathogen, currently estimated to infect two billion people throughout the world1. The bacterial population size resulting from this massive spread is very large, yet the genetic diversity within the classical members of the M. tuberculosis complex (MTBC), comprising Mycobacterium africanum, Mycobacterium bovis, Mycobacterium microti, Mycobacterium pinnipedii, and M. tuberculosis is very limited. Tuberculosis is therefore assumed to be a recent human disease2,3 linked to clonal expansion of its causative organism4-6.
In contrast to MTBC, smooth tubercle bacilli (STB), defined as clinical isolates displaying a distinctive smooth colony phenotype on culture media, named Mycobacterium canettii and/or Mycobacterium prototuberculosis7-10, are less genetically restricted. Initial genotyping analysis suggested that these isolates possess a higher diversity with traces of intraspecies horizontal gene transfer (HGT) and might therefore represent early-branching lineages of tuberculosis-causing mycobacteria. Since their first isolation by Georges Canetti in 1969, less than one hundred strains of STB have been identified. All STB have been obtained from human tuberculosis patients, mostly from (or with connection to) East-Africa8,11. Thus, a collection of a few tens of STB strains from a geographically restricted region appears to contain greater genetic diversity than the worldwide population of MTBC strains. This observation raises intriguing questions about the origin of tuberculosis and provided an opportunity to examine the molecular and evolutionary events involved in the emergence of M. tuberculosis. Herein, we describe and compare complete genome sequences of five diverse STB isolates and the physiopathological properties of these mycobacteria relative to M. tuberculosis as well as whole genome shotgun (WGS) sequences of four additional STB strains for secondary screening and confirmation purposes.
We applied multilocus sequence typing (MLST) based on 12 house-keeping genes to a panel of 55 available STB isolates and identified a total of 13 sequence types among them (Fig. 1 and Supplementary Tables 1 and 2). From analyses of the concatenated sequences we inferred a highly reticulated phylogeny, suggestive of conflicting phylogenetic signals and possible HGT among the target genes. We then selected five representative isolates for the principal comprehensive genomic analysis. This selection included the original strain isolated by George Canetti, of sequence type A and an isolate from the most prevalent group of sequence type D (both belonging to the M. canettii cluster), as well as strains from the most distant sequence types L, J and K (Fig. 1 and Supplementary Fig. 1)9.
Comparison of these five STB genomes with those of M. tuberculosis H37Rv12 and other MTBC members13 revealed a very similar overall organization between STB and MTBC, with a high percentage of syntenic genes (from 96% for strain A and 93% for strain K, compared to only 77% between M. tuberculosis H37Rv and Mycobacterium marinum, one of the phylogenetically closest non-tuberculous mycobacterial species14). No major chromosomal rearrangements or plasmids were detected (Fig. 1). Pairwise analyses between the conserved STB and MTBC genome sequences showed that all combinations had average nucleotide identities of at least 97.3%, above the 95% threshold proposed for classification into the same species15. However, the genomes of STB are 10-115 kb larger than those of the MTBC members and thus represent the largest genomes known for tubercle bacilli, although they are still much smaller than those of M. marinum (6.6 Mb)16 and the other most closely related, non-tuberculous species M. kansasii (6.4 Mb)17. Excluding repetitive sequences such as PE_PGRS- and PPE_MPTR-encoding regions, which account for ~ 8% of the coding capacity of M. tuberculosis12, STB and MTBC share > 89.3% of their genomes, representing a core genome for tubercle bacilli of > 3.938 Mb. This core comprises 96.3% of the 774 M. tuberculosis H37Rv genes predicted as being essential for in vitro growth, and all 194 genes required for mycobacterial survival during mouse infection18-20, further reflecting the close affiliations of STB and M. tuberculosis. The accessory genomes of individual STB strains harbor from 124 (strain A) to 366 genes (strains K and J) not present in MTBC members that enlarge the known pan-genome of tubercle bacilli by 890 predicted coding sequences (CDS), representing a supplement of more than 20% relative to the gene pool of M. tuberculosis (Supplementary Table 3 and Supplementary Figs 2a, b). Interestingly, only nine of these CDS were common to all five STB genomes analyzed (Supplementary Table 3 and Supplementary Fig. 2c). Conversely, 51 genes partially overlapping with genomic islands21 present in MTBC were not found in any of the STB strains (Supplementary Table 4). These genes encode derivatives of mobile elements, such as the phiRv1 and phiRv2 prophage-like regions (24 CDS), 3 transposases, 5 unique members of a glycine-rich protein family (e.g. PE_PGRS33; Supplementary Fig. 3a) and 19 other hypothetical proteins (Supplementary Fig. 3b). It is noteworthy that Rv1989c-Rv1990c from one such MTBC specific region showed around 90% identity with proteins encoded on a plasmid from Mycobacterium gilvum and Mycobacterium sp. KMS, raising intriguing questions about possible transmission routes of the corresponding genes into the MTBC. Several other MTBC-specific hypothetical proteins had no or only weak amino-acid similarity with other mycobacterial proteins (Supplementary Table 4), suggesting HGT into MTBC from distant donors after the separation from the STB lineages.
We also identified prominent, HGT-related differences in clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins (CRISPR-Cas) systems between STB and MTBC. These systems may confer adaptive immunity against phages and plasmids in bacteria and archea via repeat/spacer-derived short RNAs22. The genomes of STB strains A and D contain a single CRISPR-Cas locus encoding a system of major type III-A that is similar to that of MTBC genomes, but with a few crispr spacers in common7,10 and substantially lower sequence similarities of their Cas proteins (down to 75%) than those of the core proteins (98%-100 %) (Fig. 2.). The same genomic region in the more distant K, L and J strains is occupied by a completely different CRISPR-Cas system of a rare type-Ic variant (Fig. 2), most closely related to those of environmental actinobacteria such as Gordonia amarae or purple sulfur bacteria Thioalkalivibrio sp.. Furthermore, in strain K, the presence of a second CRISPR-Cas module of a different type Ic was identified 260 kb upstream of the other locus (Fig. 2), whose Cas proteins were most similar to those of Moorella sp. or Thiorhodovibrio sp. Finally, screening of WGS-derived sequences from STB strains E, G, H, and I, located at well-distributed intermediate positions of the STB MLST-based network (Fig. 1), revealed the existence of yet another type I-E module in strains G and I that was most closely related to those of environmental actinobacteria such as Saccharomonospora sp., while in the two remaining E and H strains, a type Ic variant similar to those of STB-J, -K and -L was found (Fig. 2a). As CRISPR-Cas systems have not been identified in non-tuberculous mycobacterial species, these different systems were most likely acquired by independent HGT events that occurred after the divergence of STB and MTBC. While it is not known whether the CRISPR systems in tubercle bacilli are functional, their disparate origins suggest that the distinct, respective crispr spacer sets might not necessarily reflect genetic records of recent encounters of tubercle bacilli with distinct phage transgressors but, instead, older traces of interaction of the respective CRISPR-Cas donor organisms with non-mycobacterial phages. The identification of a 55 kb prophage region in the WGS-derived sequence of STB-I that is large enough to encode a potentially complete virion23 (Fig. 2), which to our knowledge represents the first such finding in tubercle bacilli, provides a promising future model for testing the functionality of mycobacterial CRISPR-Cas systems on adaptive immunity against phages.
Progressive genome downsizing is a hallmark of mycobacterial pathogen evolution17,24. Therefore, the larger genome sizes of STB compared with MTBC argue for their ancestral status. Further evidence for ancestrality of STB genome structures comes from inspection of interrupted coding sequence (ICDS) orthologs, thought to reflect molecular scars inherited during pseudogenization of the MTBC genomes25,26. Among the 81 reported ICDS in MTBC, most were found to be also interrupted both in STB and more distantly related mycobacteria, suggesting evolutionary ancient mycobacterial scars (Supplementary Table 5). However, we identified four ICDS, e.g. pks8 belonging to the pks multi-gene family encoding polyketide synthases that are involved in the biosynthesis of important cell envelope lipids16,27, which were intact in the genomes of STB (in one case - rv3741/42, the region was absent from STB-J) and from the M. marinum and/or M. kansasii outgroup genomes (Supplementary Fig. 4). Thus, these scars occurred in the most recent common ancestor of the MTBC after divergence from STB-like progenitors. The opposite situation, i.e. ICDS shared by the STB genomes corresponding to intact CDS in MTBC, was never observed, further supporting the ancestral status of the STB genome structure. In addition, we detected four independent loci (narX, pks5, pknH, lppV), where a likely ancestral gene organization present both in the mycobacterial outgroups M. marinum and/or M. kansasii and in STB, was rearranged to result in a single hybrid gene and loss of intervening gene(s) in MTBC genomes (Supplementary Fig. 5), similar to what has been observed for pknH in M. africanum28.
Ancient branching of STB lineages is also consistent with the much higher numbers of single nucleotide polymorphisms (SNPs) detected among STB genomes compared to MTBC. Pairwise comparisons of the STB-D, -A, -L, -J, and -K genome sequences with the M. tuberculosis H37Rv reference uncovered 16,168 - 61,228 SNPs (Fig. 1). This amount is within the 9,525-65,744 SNP range observed among the group of STB strains alone, and up to 25-fold higher than the 741-2437 SNPs previously observed among members of the MTBC13,29,30. Consistent with MLST data (here and ref. 9), a NeighborNet analysis based on pairwise comparisons of the genome-wide SNP data showed that MTBC forms a single compact group within a much larger, reticulated network of the STB genotypes (Fig. 1). Consistently, this reticulation was even increased when WGS-derived sequence data from four additional STB strains were included (Supplementary Fig. 6), further confirming the MLST-derived phylogeny at the genome level. The Phi test for recombination was highly significant (p=10−6). Importantly, the relative compactness of the MTBC branch is additionally confirmed by the structure of the phylogenetic tree, obtained after exclusion of the genome portions affected by recombination/HGT (see further and Fig. 3a). These results thus firmly demonstrate that the worldwide MTBC population only represents a genetically homogeneous subset branching from the larger diversity of recombinogenic STB isolates. Taken together with independent lines of evidence pointing to an earlier branching, they suggest that STB lineages diverged from the common ancestor of all tubercle bacilli well before the successful clonal radiation of MTBC began.
In order to compare the impact of selection on the evolution of the STB and MTBC genomes, we calculated global ratios of non-synonymous vs synonymous SNPs (dN/dS). The genome-wide dN/dS ratio is unusually high in MTBC, which has been suggested to reflect relaxed purifying selection against non-synonymous changes that are in general slightly deleterious31. The dN/dS ratios in the different gene categories among the STB strains were only about a third of those found in the MTBC (Table 1), and are thus compatible with a much longer time of exposure of STB to purifying selection, given the time dependence of dN/dS for closely related bacteria32-34 and assuming that purifying selection pressures were the same for STB as for MTBC.
As an important exception, protective human CD4+ and CD8+ T-cell antigens and epitopes of M. tuberculosis have been described to be under purifying selection, suggesting that MTBC members do not use T-cell antigen variation to escape human immune responses but, instead, might benefit from recognition by T-cells30. Similarly, we found that the dN/dS ratios based on pairwise, concatenated codon alignments35,36 of the 65 T-cell antigen-encoding STB genes conserved across all STB genomes were on average lower than those of the 2,300 genes classified as non essential and similar to slightly lower than the 710 essential genes18 conserved among all STB (Table 1). Overall similar results were also obtained when only the epitope regions of the T-cell antigens were considered. Thus, like the subset of essential proteins, human T-cell antigens tend to be more conserved in STB relative to the rest of the proteome. Following the argument of Comas and colleagues30, this sequence conservation suggests that STB and MTBC might have inherited a common strategy of immune subversion of the human host that predates the clonal emergence of the MTBC. However, there may be alternative explanations, as most of these low dN/dS antigens are also highly conserved in the environmental, facultative pathogens M. marinum and/or M. kansasii and/or other mycobacteria. For example, the 6-kD early secreted antigenic target (ESAT-6, Rv3875) and the 34.6-kD secreted antigen 85B (Ag85B; Rv1886c), that both show 100% amino-acid conservation in MTBC and STB, have orthologues in M. marinum that show 91% (ESAT-6) and 89% (Ag85B) amino-acid identity, which is above the average overall pairwise identity of 85.2%16. The conservation of these proteins might thus also be explained by their role in host-pathogen interaction such as phagosomal rupture37, cell envelope stability38 or other functions that are not necessarily linked to interactions with human T-cells.
Extensive recombination among STB, revealed by our comparative genome analysis, might also have played a role in the discrepancy in dN/dS between the MTBC and STB groups, as it could more efficiently oppose fixation of slightly deleterious mutations than in the more clonal MTBC population39. Consistent with this contention, strong variations in the local distribution of SNPs were observed throughout the aligned STB and MTBC genomes, suggestive of numerous recombination events. Approximately one-third of the core genome alignment consists of zones with significantly lower or higher SNP density compared to expectations for predicted recombination-free nucleotide differences between each pair of genomes. A stringent selection of informative regions among the predicted recombination-free blocks led to a minimal clonal backbone of 1,794,643 characters (~33% of the core genome), which was used to infer a phylogenetic tree (Fig. 3a). Inspection of the genomic regions with unexpected SNP densities allowed us to identify > 110 blocks, of up to 14 kb and including each from 1 to 5 complete genes (Supplementary Table 6), with homoplasic SNP distributions (relative to the tree), indicative of likely inter-strain recombination events among STB and/or between STB and MTBC strains (Fig. 3b and Supplementary Fig. 7a, b). The extensive impact of recombination was independently confirmed by the finding that ~8% to ~15% of the protein coding sequence alignments from the core genome has mosaic structures indicative of inter-strain intragenic recombination events. In contrast, the influence of exogenous importation from more distant mycobacterial species on the core genome sequence diversity is apparently minimal, as inferred by the detection of only few regions with unexpectedly high SNP densities in STB strains, yielding BLAST best hits closer to non-tuberculous mycobacteria than to STB and MTBC (Supplementary Fig. 7c).
Remarkably, the gene blocks in M. tuberculosis whose sequences perfectly match those of one or more STB strains, showed SNPs in the orthologous region in M. bovis and/or other MTBC strains (Fig. 3b), suggesting that gene fluxes between M. tuberculosis and the STB strain pool existed even well after the divergence of the MTBC, and perhaps still exist. We also found intermediate situations, where the SNP distribution clearly suggests recombination events that were more ancient and likely followed by accumulation of a few mutations in the recipient or the donor strains (Supplementary Fig. 7a). These data provide new, solid evidence to the question of inter-strain gene flux in M. tuberculosis40,41. Our findings also raise puzzling questions on the (micro-) environments and mechanisms favoring or having favored such extensive DNA exchanges. The high number of apparently recent recombination episodes, as suggested by numerous perfect large sequence matches detected among sequences from different STB lineages together with the almost exclusive isolation of STB strains from patients around the Horn of Africa strongly suggests a common local source. Aquatic environments rich in mycobacteria, potentially residing in protozoan hosts24,42, are one possible opportunity for genetic exchange to occur, as suggested by a recent report on detection of MTBC DNA in rural water sources in Ethiopia (E. Wellington, personal communication, Abstract, 16th International Symposium on the Biology of Actinomycetes). The presence of a 55 kb genomic segment corresponding to a putative complete phage-encoding region inserted into the Lys tRNA gene of strain STB-I (Fig. 2) suggests a possible mediation by phages, although alternative mechanisms such as DNA transfer by conjugation, reported for Mycobacterium smegmatis under biofilm conditions43, could also be involved.
To determine whether the genome differences between the STB and MTBC strains impact on host-pathogen interactions, we first measured their growth in in vitro cultures. Most STB grew 2 to 3 times faster than M. tuberculosis both in liquid (Supplementary Fig. 8) and on solid media (data not shown) at 30°c and 37°c, in line with previous observations10,11. Upon infection of BALB/c mice (Fig. 4) and C57BL6 mice (Supplementary Fig. 9) by aerosol, the STB strains effectively multiplied in lungs and disseminated to the spleens during the acute infection phase, but consistently persisted less well during the chronic infection phase compared to M. tuberculosis. While the latter was able to persist in the lungs for up to 30 weeks at levels close to those of the acute phase (peaking at 3 weeks with around 107.7 colony-forming units (CFUs)), the infection levels of all STB strains dropped by at least 1 log at all (and by 2 to 3 logs at most) later time points in these organs (p=0.05 by Mann-Whitney test, except for day 130 for strains D, L and K). The strongest difference with M. tuberculosis was observed for strain K, the strain phylogenomically most distant from MTBC and for which bacterial counts were undetectable after 30 weeks in BALB/c mice (Fig. 4c, d). Similar trends were observed in spleens, with strain K also almost completely cleared at day 210. In parallel, histopathological analyses revealed less intense lung lesions and inflammation 128 days after infection with the STB strains compared to M. tuberculosis infection, with strain K showing the least damages (Fig. 4 and Supplementary Table 7). Furthermore, C57BL/6 mice intravenously infected with high doses of STB survived in contrast to controls infected with M. tuberculosis strains of different lineages (data not shown), confirming decreased virulence of STB.
Finally, we determined whether these variations could be correlated to differences in innate or adaptive immune responses elicited by infection. The STB and M. tuberculosis strains were similarly able to induce maturation of innate immunity cells in vitro, such as dentritic cells derived from C57BL/6, tlr2°/°, tlr4°/°, or double KO mice (data not shown), suggesting shared major Pathogen-Associated Molecular Patterns (PAMPs)44. Consistently, substantial recruitment of activated innate immune cells, i.e., CD11b+ BST-2+ (Bone Marrow Stromal Cell Antigen-2)+ and CD11c+ MHC-IIhi, was observed in vivo in the lung parenchyma of SCID mice after 3 weeks of infection by STB, but to a lower extent as compared to M. tuberculosis infection (data not shown). Concerning adaptive responses, massive recruitment of activated CD4+ and CD8+ T-cells, displaying CD44 modulation and CD45RB, CD27, CD62L downregulation, was detected in the lungs of C57BL/6 mice after 13 weeks of infection by smooth strains. Again, the responses were overall quantitatively lower for STB compared to M. tuberculosis strains, especially for STB-K (Supplementary Fig. 10), in line with the lower virulence and persistence of STB.
With the larger pan genome reflecting the ancestral, wider gene pool of tubercle bacilli, their, lower virulence and faster growth especially at temperatures below 37°C, plausibly reflecting broader environmental adaptability, STB strains might thus come nearer to the as yet unknown missing link between the obligate pathogen M. tuberculosis and environmental mycobacteria. We propose that M. tuberculosis has evolved its so successful widespread, pathogenic lifestyle starting from a pool of STB-like mycobacteria by gaining additional virulence and persistence mechanisms through a potential combination of i) loss of gene function, ii) acquisition of novel genes via HGT, iii) inter-strain recombination of gene clusters and (iv) fixation of SNPs. From the data presented here, a rational experimental design to elucidate which of these genetic events were involved can now be undertaken. Primary candidates are MTBC-specific genes (Supplementary Table 4), including prophage-like phiRv1 / phiRv2 encoding regions reported to be important for late infection45, genes encoding PE_PGRS33 or other MTBC-acquired PE/PPE proteins known to enhance cellular toxicity46, polyketide synthase Pks8/17, the large prophage region in STB-I and/or CRISPR-cas systems. The insights gained through our analysis thus open novel perspectives to identify new targets to combat tuberculosis infection and disease.
The 55 STB and 10 reference MTBC isolates are described in Supplementary Table 1. Twelve house-keeping genes were selected for MLST47 (Supplementary Table 2). Phylogenetic groupings were identified by split decomposition analysis48 on the concatenated target sequences .
Genomic DNA was extracted from cultured single bacterial colonies12. For genome sequencing of STB-D, -J, -K and -L, Sanger reads from 10-kb fragment shotgun libraries at 4- 4.9 fold coverages were assembled with contigs obtained from Newbler assemblies of 454/Roche reads at 13-18.1 fold coverages, using Arachne49. Scaffolds were validated using Mekano interface (Genoscope). Primer walking, PCRs and in-vitro transposition were used for finishing. The assembled consensus sequences were validated using Illumina reads at 45-105 fold coverages and consed functionalities, and by mapping of termini-sequences from bacterial artificial chromosome libraries50. High quality, contiguous genome sequences of 4420 kb (STB-L, 9 contigs), 4432 kb (STB-D, 12 contigs), 4524 kb (STB-J, 11 contigs), and 4525 kb (STB-K, 9 contigs) were generated. Remaining gaps estimated not to exceed 2 kb correspond to GC-rich and repetitive regions coding for PE_PGRS proteins, and/or the pks5 region (STB-J). For STB-A, a fully finished, contiguous sequence of 4,482,059 bp was obtained by using ~ 80,000 shotgun Sanger reads, Illumina-generated reads and finishing12,28. WGS data from STB strains E, G, H, and I were generated using Illumina HiSseq technology and single lanes. Resulting reads that covered the genomes of these STB up to 900x were assembled using the Velvet software51 and contigs were ordered using M. canettii CIPT 140010059 (STB-A) and M. tuberculosis H37Rv as reference genomes.
Annotation and genome comparisons were performed with the Microscope platform52, Artemis and Artemis comparison tool (ACT)53. When applicable, annotations were transferred from those of M. tuberculosis orthologs in the TubercuList/Mycobrowser database, using BLAST matches of > 90% protein sequence identity, an alignable region of > 80% of the shortest protein length in pairwise comparisons and visual inspection of the gene synteny. Pairwise average nucleotide identities were calculated using JSpecies54. The core/accessory genomes of STB and M. tuberculosis were determined as described16.
SNiPer pipeline (Genoscope) based on the SSAHA2 package55 was used to map Illumina reads and detect SNPs and indels of STB strains against a corrected version56 of the M. tuberculosis H37Rv reference sequence (NC_000962)12. After exclusion of ambiguous maps on repeat regions, an average of 4.7 million split paired-end reads of 36 bp (STB-A, -D, -L, - J) or trimmed at 50 bp (STB-K) were mapped at a resulting genome coverage > 40x. SNPs with base coverage < 10, base quality < 25, or heterozygosity > 0.2 were removed. ACT53 comparison files were created by using MUMmer and NUCmer softwares57 to visualize the SNP distribution in local genome regions.
dN/dS ratios were calculated on orthologs conserved in all STB and M. tuberculosis H37Rv, as identified by bidirectional best hits, alignable region of >80% and sequence identity >=30%. Pairwise, concatenated codon alignments between M. tuberculosis H37Rv and each STB strain were generated using PAL2NAL58, after respective protein alignments obtained with MUSCLE59. Synonymous and non-synonymous substitutions were defined using Nei-Gojobori method-based SNAP35 or maximum likelihood-based PAML36. STB T-cell antigen, essential and non-essential gene categories, as well as T-cell epitope codon concatenates were constructed as described30.
The genomes of M. tuberculosis H37Rv, M. bovis AF2122/97 and the five STB strains were aligned using progressive Mauve60. Given a pair ij of aligned genomes, the number of SNPs xij observed between i and j within a region of length l follows a binomial distribution B(l, pij), where pij is the expected proportion of recombination-free nucleotide differences between taxa i and j. Regions containing at least one pair of sequences ij with an unexpectedly large or low number xij of SNPs, i.e. min [Pr (X ≥ xij ), Pr (X ≤ xij )] < 0.05 where X ~ B(l, pij), were identified by using a 200 character-long sliding window along the conserved (core) portions of the multiple genome alignment. The value pij inside each window was estimated as the proportion of SNPs between i and j within the 10000 aligned characters flanking the sliding window on both sides. To obtain a reference phylogeny, all regions of length ≥ 500 characters (excluding gaps) that did not contain an unexpected number of SNPs were concatenated. The derived supermatrix was used to infer a phylogenetic tree by using the Neighbor-Joining algorithm on the pairwise nucleotide p-distances61. All regions of length ≥ 500 characters with a significantly high or low number of SNPs were inspected visually for detection of concentration of homoplasic characters using ACT53, leading to similarities between strains incongruent with the phylogenetic tree. The proportion of protein coding sequences within the core genome likely affected by inter-strain recombination was assessed with the Pairwise Homoplasy Index62, the Maximum x2 test63, and the Neighbour Similarity Score64.
Growth rates of STB and reference MTBC strains in liquid media were measured by using a BACTEC 460 system (Beckton-Dickinson) as recommended by the manufacturer.
Mice were maintained according to the Institut Pasteur de Lille and Paris guidelines for laboratory animal husbandry. Animal experiments were approved by the Nord-Pas-De-Calais ethical committee (CEEA 15/2009) and the Institut Pasteur Hygiene Committee (authorization number 75-1469), in accordance with European and French guidelines (Directive 86/609/CEE and Decree 87–848). Eight-week-old female BALB/c mice were infected by the intranasal route with 103 CFUs of either STB or M. tuberculosis H37Rv strains, respectively. At indicated times, 4 mice per group were sacrificed, and colony counting was performed from homogenized individual lungs and spleens as described65. For histopathological evaluation, whole lungs were harvested from 3 BALB/c mice per group 128 days post-infection, fixed in 4% formalin, and embedded in paraffin. Four mm-thick sections were stained with hematoxylin-eosin. Virulence and cell-analysis-based immunological assays using C57BL/6 and/or SCID mice were performed as described66,67. Adaptive immune cells from infected mice were prepared, incubated with conjugated mAbs (Beckton-Dickinson), fixed, and analyzed using a CyAn system and Summit (Beckman Coulter) and FlowJo (Treestar) softwares.
FN2URLs Magnifying Genome (MaGe) server, https://www.genoscope.cns.fr/agc/microscope/about/collabprojects.php?P_id=44#ancreLogin; MycoBrowser database, http://mycobrowser.epfl.ch/; The Mycobacteriophage database, http://phagesdb.org/.
FN3Accession Codes The complete genome sequence for strain STB-A (CIPT 140010059) was deposited under Accession No. HE572590. Genome sequences of strains STB-D (CIPT 140060008), STB-J (CIPT 140070017), STB-K (CIPT 140070010), and STB-L (CIPT 140070008) were deposited in the EMBL database under accession numbers PRJEB94 - BN44 - FO203507, PRJEB93 - BN43 - FO203508, PRJEB92 - BN42 - FO203509, and PRJEB95 - BN45 - FO203510, respectively. Illumina-derived WGS sequences for strains STB-E (CIPT 140070002), STB-G (CIPT 140070005), STB-H (CIPT 140070013), and STB-I (CIPT 140070007) were deposited in the EMBL WGS repository under project numbers PRJEB584, PRJEB585, PRJEB586, and PRJEB587, respectively.
We are grateful to Stewart Cole for help in initiating the M. canettii CIPT 140010059 genome sequencing project and advice, to Thierry Garnier and Adamandia Kapopoulou for help in data management and to Zoe Rouy for help with sequence deposal. The work was supported in part by the Institut Pasteur (PTR 314, PTR 383), the European Community’s FP7 program Grant no. 260872, the Wellcome trust grant 098051, and a Genoscope collaborative grant no. 114.
|1.||Dye C,Williams BG. The population dynamics and control of tuberculosisScienceYear: 201032885686120466923|
|2.||Sreevatsan S,et al. Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global disseminationProc Natl Acad Sci U S AYear: 199794986998749275218|
|3.||Wirth T,et al. Origin, spread and demography of the Mycobacterium tuberculosis complexPLoS PathogYear: 20084e100016018802459|
|4.||Brosch R,et al. A new evolutionary scenario for the Mycobacterium tuberculosis complexProc Natl Acad Sci U S AYear: 2002993684368911891304|
|5.||Supply P,et al. Linkage disequilibrium between minisatellite loci supports clonal evolution of Mycobacterium tuberculosis in a high tuberculosis incidence areaMol MicrobiolYear: 20034752953812519202|
|6.||Hirsh AE,Tsolaki AG,DeRiemer K,Feldman MW,Small PM. Stable association between strains of Mycobacterium tuberculosis and their human host populationsProc Natl Acad Sci U S AYear: 20041014871487615041743|
|7.||van Soolingen D,et al. A novel pathogenic taxon of the Mycobacterium tuberculosis complex, Canetti: characterization of an exceptional isolate from AfricaInt J Syst BacteriolYear: 199747123612459336935|
|8.||Fabre M,et al. High genetic diversity revealed by variable-number tandem repeat genotyping and analysis of hsp65 gene polymorphism in a large collection of “Mycobacterium canettii” strains indicates that the M. tuberculosis complex is a recently emerged clone of “M. canettii”J Clin MicrobiolYear: 2004423248325515243089|
|9.||Gutierrez MC,et al. Ancient origin and gene mosaicism of the progenitor of Mycobacterium tuberculosisPLoS PathogYear: 20051e516201017|
|10.||Fabre M,et al. Molecular characteristics of “Mycobacterium canettii” the smooth Mycobacterium tuberculosis bacilliInfect Genet EvolYear: 2010101165117320692377|
|11.||Koeck JL,et al. Clinical characteristics of the smooth tubercle bacilli ‘Mycobacterium canettii’ infection suggest the existence of an environmental reservoirClin Microbiol InfectYear: 2011171013101920831613|
|12.||Cole ST,et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequenceNatureYear: 19983935375449634230|
|13.||Garnier T,et al. The complete genome sequence of Mycobacterium bovisProc Natl Acad Sci U S AYear: 20031007877788212788972|
|14.||Springer B,Stockman L,Teschner K,Roberts GD,Bottger EC. Two-laboratory collaborative study on identification of mycobacteria: molecular versus phenotypic methodsJ. Clin. MicrobiolYear: 1996342963038789004|
|15.||Goris J,et al. DNA-DNA hybridization values and their relationship to whole-genome sequence similaritiesInt J Syst Evol MicrobiolYear: 200757819117220447|
|16.||Stinear TP,et al. Insights from the complete genome sequence of Mycobacterium marinum on the evolution of Mycobacterium tuberculosisGenome ResYear: 20081872974118403782|
|17.||Veyrier FJ,Dufort A,Behr MA. The rise and fall of the Mycobacterium tuberculosis genomeTrends MicrobiolYear: 20111915616121277778|
|18.||Sassetti CM,Boyd DH,Rubin EJ. Genes required for mycobacterial growth defined by high density mutagenesisMol MicrobiolYear: 200348778412657046|
|19.||Sassetti CM,Rubin EJ. Genetic requirements for mycobacterial survival during infectionProc Natl Acad Sci U S AYear: 2003100129891299414569030|
|20.||Griffin JE,et al. High-resolution phenotypic profiling defines genes essential for mycobacterial growth and cholesterol catabolismPLoS PathogYear: 20117e100225121980284|
|21.||Becq J,et al. Contribution of horizontally acquired genomic islands to the evolution of the tubercle bacilliMol Biol EvolYear: 2007241861187117545187|
|22.||Makarova KS,et al. Evolution and classification of the CRISPR-Cas systemsNat Rev MicrobiolYear: 2011946747721552286|
|23.||Hendrix RW,Smith MC,Burns RN,Ford ME,Hatfull GF. Evolutionary relationships among diverse bacteriophages and prophages: all the world’s a phageProc Natl Acad Sci U S AYear: 1999962192219710051617|
|24.||Gordon SV,Bottai D,Simeone R,Stinear TP,Brosch R. Pathogenicity in the tubercle bacillus: molecular and evolutionary determinantsBioessaysYear: 20093137838819274661|
|25.||Deshayes C,et al. Detecting the molecular scars of evolution in the Mycobacterium tuberculosis complex by analyzing interrupted coding sequencesBMC Evol BiolYear: 200887818325090|
|26.||Smith NH,Hewinson RG,Kremer K,Brosch R,Gordon SV. Myths and misconceptions: the origin and evolution of Mycobacterium tuberculosisNat Rev MicrobiolYear: 2009753754419483712|
|27.||Reed MB,et al. A glycolipid of hypervirulent tuberculosis strains that inhibits the innate immune responseNatureYear: 2004431848715343336|
|28.||Bentley SD,et al. The genome of Mycobacterium africanum West African 2 reveals a lineage-specific locus and genome erosion common to the M. tuberculosis complexPLoS Negl Trop DisYear: 20126e155222389744|
|29.||Brosch R,et al. Genome plasticity of BCG and impact on vaccine efficacyProc Natl Acad Sci U S AYear: 20071045596560117372194|
|30.||Comas I,et al. Human T cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconservedNat GenetYear: 20104249850320495566|
|31.||Hershberg R,et al. High functional diversity in M. tuberculosis driven by genetic drift and human demographyPLoS BiolYear: 20086e31119090620|
|32.||Rocha EP,et al. Comparisons of dN/dS are time dependent for closely related bacterial genomesJ Theor BiolYear: 200623922623516239014|
|33.||Castillo-Ramirez S,et al. The impact of recombination on dN/dS within recently emerged bacterial clonesPLoS PathogYear: 20117e100212921779170|
|34.||Croucher NJ,et al. Rapid pneumococcal evolution in response to clinical interventionsScienceYear: 201133143043421273480|
|35.||Korber B. Rodrigo AG,Learn GHComputational Analysis of HIV Molecular SequencesYear: 20005572Kluwer Academic Publishers Ch. 4.|
|36.||Goldman N,Yang Z. A codon-based model of nucleotide substitution for protein-coding DNA sequencesMol Biol EvolYear: 1994117257367968486|
|37.||Simeone R,et al. Phagosomal rupture by Mycobacterium tuberculosis results in toxicity and host cell deathPLoS PathogYear: 20128e100250722319448|
|38.||Kalscheuer R,Weinrick B,Veeraraghavan U,Besra GS,Jacobs WR Jr.. Trehalose-recycling ABC transporter LpqY-SugA-SugB-SugC is essential for virulence of Mycobacterium tuberculosisProc Natl Acad Sci U S AYear: 2010107217612176621118978|
|39.||Felsenstein J. The evolutionary advantage of recombinationGeneticsYear: 1974787377564448362|
|40.||Achtman M. Insights from genomic comparisons of genetically monomorphic bacterial pathogensPhilos Trans R Soc Lond B Biol SciYear: 201236786086722312053|
|41.||Namouchi A,Didelot X,Schock U,Gicquel B,Rocha EP. After the bottleneck: Genome-wide diversification of the Mycobacterium tuberculosis complex by mutation, recombination, and natural selectionGenome ResYear: 20122272173422377718|
|42.||Mba Medie F,Ben Salah I,Henrissat B,Raoult D,Drancourt M. Mycobacterium tuberculosis complex mycobacteria as amoeba-resistant organismsPLoS OneYear: 20116e2049921673985|
|43.||Nguyen KT,Piastro K,Gray TA,Derbyshire KM. Mycobacterial biofilms facilitate horizontal DNA transfer between strains of Mycobacterium smegmatisJ BacteriolYear: 20101925134514220675473|
|44.||Medzhitov R,Janeway CA. Innate Immunity cecognition and control of adaptive immune responsesSemin ImmunolYear: 1998103513539799709|
|45.||Aagaard C,et al. A multistage tuberculosis vaccine that confers efficient protection before and after exposureNat MedYear: 20111718919421258338|
|46.||Cadieux N,et al. Induction of cell death after localization to the host cell mitochondria by the Mycobacterium tuberculosis PE_PGRS33 proteinMicrobiologyYear: 201115779380421081760|
|47.||Maiden MC,et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganismsProc Natl Acad Sci U S AYear: 199895314031459501229|
|48.||Huson DH,Bryant D. Application of phylogenetic networks in evolutionary studiesMol Biol EvolYear: 20062325426716221896|
|49.||Batzoglou S,et al. ARACHNE: a whole-genome shotgun assemblerGenome ResYear: 20021217718911779843|
|50.||Brosch R,et al. Use of a Mycobacterium tuberculosis H37Rv bacterial artificial chromosome library for genome mapping, sequencing, and comparative genomicsInfect ImmunYear: 199866222122299573111|
|51.||Zerbino DR,Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphsGenome ResYear: 20081882182918349386|
|52.||Vallenet D,et al. MicroScope: a platform for microbial genome annotation and comparative genomicsDatabase (Oxford)Year: 2009|
|53.||Carver T,et al. Artemis and ACT: viewing, annotating and comparing sequences stored in a relational databaseBioinformaticsYear: 2008242672267618845581|
|54.||Richter M,Rossello-Mora R. Shifting the genomic gold standard for the prokaryotic species definitionProc Natl Acad Sci U S AYear: 2009106191261913119855009|
|55.||Ning Z,Cox AJ,Mullikin JC. SSAHA: a fast search method for large DNA databasesGenome ResYear: 2001111725172911591649|
|56.||Niemann S,et al. Genomic diversity among drug sensitive and multidrug resistant isolates of Mycobacterium tuberculosis with identical DNA fingerprintsPLoS OneYear: 20094e740719823582|
|57.||Delcher AL,Phillippy A,Carlton J,Salzberg SL. Fast algorithms for large-scale genome alignment and comparisonNucleic Acids ResYear: 2002302478248312034836|
|58.||Suyama M,Torrents D,Bork P. PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignmentsNucleic Acids ResYear: 200634W60961216845082|
|59.||Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughputNucleic Acids ResYear: 2004321792179715034147|
|60.||Darling AE,Mau B,Perna NT. progressiveMauve: multiple genome alignment with gene gain, loss and rearrangementPLoS OneYear: 20105e1114720593022|
|61.||Saitou N,Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic treesMol Biol EvolYear: 198744064253447015|
|62.||Bruen TC,Philippe H,Bryant D. A simple and robust statistical test for detecting the presence of recombinationGeneticsYear: 20061722665268116489234|
|63.||Smith JM. Analyzing the mosaic structure of genesJ Mol EvolYear: 1992341261291556748|
|64.||Jakobsen IB,Easteal S. A program for calculating and displaying compatibility matrices as an aid in determining reticulate evolution in molecular sequencesComput Appl BiosciYear: 1996122912958902355|
|65.||Pethe K,et al. The heparin-binding haemagglutinin of M. tuberculosis is required for extrapulmonary disseminationNatureYear: 200141219019411449276|
|66.||Majlessi L,et al. Influence of ESAT-6 secretion system 1 (RD1) of Mycobacterium tuberculosis on the interaction between mycobacteria and the host immune systemJ ImmunolYear: 20051743570357915749894|
|67.||Bottai D,et al. ESAT-6 secretion-independent impact of ESX-1 genes espF and espG1 on virulence of Mycobacterium tuberculosisJ Infect DisYear: 20112031155116421196469|
Previous Document: Common variants in the HLA-DRB1-HLA-DQA1 HLA class II region are associated with susceptibility to v...
Next Document: Genome-wide association analysis identifies new susceptibility loci for Behçet's disease and epista...