Thierry Garnier

Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence.

Countless millions of people have died from tuberculosis, a chronic infectious disease caused by the tubercle bacillus. The complete genome sequence of the best-characterized strain of Mycobacterium tuberculosis, H37Rv, has been determined and analysed in order to improve our understanding of the biology of this slow-growing pathogen and to help the conception of new prophylactic and therapeutic interventions. The genome comprises 4,411,529 base pairs, contains around 4,000 genes, and has a very high guanine + cytosine content that is reflected in the biased amino-acid content of the proteins. M. tuberculosis differs radically from other bacteria in that a very large portion of its coding capacity is devoted to the production of enzymes involved in lipogenesis and lipolysis, and to two new families of glycine-rich proteins with a repetitive structure that may represent a source of antigenic variation.

Massive gene decay in the leprosy bacillus.

Leprosy, a chronic human neurological disease, results from infection with the obligate intracellular pathogen Mycobacterium leprae, a close relative of the tubercle bacillus. Mycobacterium leprae has the longest doubling time of all known bacteria and has thwarted every effort at culture in the laboratory. Comparing the 3.27-megabase (Mb) genome sequence of an armadillo-derived Indian isolate of the leprosy bacillus with that of Mycobacterium tuberculosis (4.41 Mb) provides clear explanations for these properties and reveals an extreme case of reductive evolution. Less than half of the genome contains functional genes but pseudogenes, with intact counterparts in M. tuberculosis, abound. Genome downsizing and the current mosaic arrangement appear to have resulted from extensive recombination events between dispersed repetitive sequences. Gene deletion and decay have eliminated many important metabolic activities including siderophore production, part of the oxidative and most of the microaerophilic and anaerobic respiratory chains, and numerous catabolic systems and their regulatory circuits.

A new evolutionary scenario for the Mycobacterium tuberculosis complex

The distribution of 20 variable regions resulting from insertion-deletion events in the genomes of the tubercle bacilli has been evaluated in a total of 100 strains of Mycobacterium tuberculosis, Mycobacterium africanum, Mycobacterium canettii, Mycobacterium microti, and Mycobacterium bovis. This approach showed that the majority of these polymorphisms did not occur independently in the different strains of the M. tuberculosis complex but, rather, resulted from ancient, irreversible genetic events in common progenitor strains. Based on the presence or absence of an M. tuberculosis specific deletion (TbD1), M. tuberculosis strains can be divided into ancestral and “modern” strains, the latter comprising representatives of major epidemics like the Beijing, Haarlem, and African M. tuberculosis clusters. Furthermore, successive loss of DNA, reflected by region of difference 9 and other subsequent deletions, was identified for an evolutionary lineage represented by M. africanum, M. microti, and M. bovis that diverged from the progenitor of the present M. tuberculosis strains before TbD1 occurred. These findings contradict the often-presented hypothesis that M. tuberculosis, the etiological agent of human tuberculosis evolved from M. bovis, the agent of bovine disease. M. canettii and ancestral M. tuberculosis strains lack none of these deleted regions, and, therefore, seem to be direct descendants of tubercle bacilli that existed before the M. africanum→M. bovis lineage separated from the M. tuberculosis lineage. This observation suggests that the common ancestor of the tubercle bacilli resembled M. tuberculosis or M. canettii and could well have been a human pathogen already.

The complete genome sequence of Mycobacterium bovis

Mycobacterium bovis is the causative agent of tuberculosis in a range of animal species and man, with worldwide annual losses to agriculture of

Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays

3 billion. The human burden of tuberculosis caused by the bovine tubercle bacillus is still largely unknown. M. bovis was also the progenitor for the M. bovis bacillus Calmette–Guérin vaccine strain, the most widely used human vaccine. Here we describe the 4,345,492-bp genome sequence of M. bovis AF2122/97 and its comparison with the genomes of Mycobacterium tuberculosis and Mycobacterium leprae. Strikingly, the genome sequence of M. bovis is >99.95% identical to that of M. tuberculosis, but deletion of genetic information has led to a reduced genome size. Comparison with M. leprae reveals a number of common gene losses, suggesting the removal of functional redundancy. Cell wall components and secreted proteins show the greatest variation, indicating their potential role in host–bacillus interactions or immune evasion. Furthermore, there are no genes unique to M. bovis, implying that differential gene expression may be the key to the host tropisms of human and bovine bacilli. The genome sequence therefore offers major insight on the evolution, host preference, and pathobiology of M. bovis.

Genome plasticity of BCG and impact on vaccine efficacy.

Whole‐genome comparisons of the tubercle bacilli were undertaken using ordered bacterial artificial chromosome (BAC) libraries of Mycobacterium tuberculosis and the vaccine strain, Mycobacterium bovis BCG‐Pasteur, together with the complete genome sequence of M. tuberculosis H37Rv. Restriction‐digested BAC arrays of M. tuberculosis H37Rv were used in hybridization experiments with radiolabelled M. bovis BCG genomic DNA to reveal the presence of 10 deletions (RD1–RD10) relative to M. tuberculosis. Seven of these regions, RD4–RD10, were also found to be deleted from M. bovis, with the three M. bovis BCG‐specific deletions being identical to the RD1–RD3 loci described previously. The distribution of RD4–RD10 in Mycobacterium africanum resembles that of M. tuberculosis more closely than that of M. bovis, whereas an intermediate arrangement was found in Mycobacterium microti, suggesting that the corresponding genes may affect host range and virulence of the various tubercle bacilli. Among the known products encoded by these loci are a copy of the proposed mycobacterial invasin Mce, three phospholipases, several PE, PPE and ESAT‐6 proteins, epoxide hydrolase and an insertion sequence. In a complementary approach, direct comparison of BACs uncovered a third class of deletions consisting of two M. tuberculosis H37Rv loci, RvD1 and RvD2, deleted from the genome relative to M. bovis BCG and M. bovis. These deletions affect a further seven genes, including a fourth phospholipase, plcD. In summary, the insertions and deletions described here have important implications for our understanding of the evolution of the tubercle complex.

Comparative genomic and phylogeographic analysis of Mycobacterium leprae.

To understand the evolution, attenuation, and variable protective efficacy of bacillus Calmette–Guérin (BCG) vaccines, Mycobacterium bovis BCG Pasteur 1173P2 has been subjected to comparative genome and transcriptome analysis. The 4,374,522-bp genome contains 3,954 protein-coding genes, 58 of which are present in two copies as a result of two independent tandem duplications, DU1 and DU2. DU1 is restricted to BCG Pasteur, although four forms of DU2 exist; DU2-I is confined to early BCG vaccines, like BCG Japan, whereas DU2-III and DU2-IV occur in the late vaccines. The glycerol-3-phosphate dehydrogenase gene, glpD2, is one of only three genes common to all four DU2 variants, implying that BCG requires higher levels of this enzyme to grow on glycerol. Further amplification of the DU2 region is ongoing, even within vaccine preparations used to immunize humans. An evolutionary scheme for BCG vaccines was established by analyzing DU2 and other markers. Lesions in genes encoding σ-factors and pleiotropic transcriptional regulators, like PhoR and Crp, were also uncovered in various BCG strains; together with gene amplification, these affect gene expression levels, immunogenicity, and, possibly, protection against tuberculosis. Furthermore, the combined findings suggest that early BCG vaccines may even be superior to the later ones that are more widely used.

Control of M. tuberculosis ESAT-6 Secretion and Specific T Cell Recognition by PhoP

Reductive evolution and massive pseudogene formation have shaped the 3.31-Mb genome of Mycobacterium leprae, an unculturable obligate pathogen that causes leprosy in humans. The complete genome sequence of M. leprae strain Br4923 from Brazil was obtained by conventional methods (6× coverage), and Illumina resequencing technology was used to obtain the sequences of strains Thai53 (38× coverage) and NHDP63 (46× coverage) from Thailand and the United States, respectively. Whole-genome comparisons with the previously sequenced TN strain from India revealed that the four strains share 99.995% sequence identity and differ only in 215 polymorphic sites, mainly SNPs, and by 5 pseudogenes. Sixteen interrelated SNP subtypes were defined by genotyping both extant and extinct strains of M. leprae from around the world. The 16 SNP subtypes showed a strong geographical association that reflects the migration patterns of early humans and trade routes, with the Silk Road linking Europe to China having contributed to the spread of leprosy.

Bacterial Artificial Chromosome-Based Comparative Genomic Analysis Identifies Mycobacterium microti as a Natural ESAT-6 Deletion Mutant

Analysis of mycobacterial strains that have lost their ability to cause disease is a powerful approach to identify yet unknown virulence determinants and pathways involved in tuberculosis pathogenesis. Two of the most widely used attenuated strains in the history of tuberculosis research are Mycobacterium bovis BCG (BCG) and Mycobacterium tuberculosis H37Ra (H37Ra), which both lost their virulence during in vitro serial passage. Whereas the attenuation of BCG is due mainly to loss of the ESAT-6 secretion system, ESX-1, the reason why H37Ra is attenuated remained unknown. However, here we show that a point mutation (S219L) in the predicted DNA binding region of the regulator PhoP is involved in the attenuation of H37Ra via a mechanism that impacts on the secretion of the major T cell antigen ESAT-6. Only H37Ra “knock-ins” that carried an integrated cosmid with the wild-type phoP gene from M. tuberculosis H37Rv showed changes in colony morphology, increased virulence, ESAT-6 secretion, and induction of specific T cell responses, whereas other H37Ra constructs did not. This finding established a link between the PhoP regulator and ESAT-6 secretion that opens exciting new perspectives for elucidating virulence regulation in M. tuberculosis.

Comparative genomics uncovers large tandem chromosomal duplications in Mycobacterium bovis BCG Pasteur

ABSTRACT Mycobacterium microti is a member of the Mycobacterium tuberculosis complex that causes tuberculosis in voles. Most strains of M. microti are harmless for humans, and some have been successfully used as live tuberculosis vaccines. In an attempt to identify putative virulence factors of the tubercle bacilli, genes that are absent from the avirulent M. microti but present in human pathogen M. tuberculosis or Mycobacterium bovis were searched for. A minimal set of 50 bacterial artificial chromosome (BAC) clones that covers almost all of the genome of M. microti OV254 was constructed, and individual BACs were compared to the corresponding BACs from M. bovis AF2122/97 and M. tuberculosis H37Rv. Comparison of pulsed-field gel-separated DNA digests of BAC clones led to the identification of 10 regions of difference (RD) between M. microti OV254 and M. tuberculosis. A 14-kb chromosomal region (RD1mic) that partly overlaps the RD1 deletion in the BCG vaccine strain was missing from the genomes of all nine tested M. microti strains. This region covers 13 genes, Rv3864 to Rv3876, in M. tuberculosis, including those encoding the potent ESAT-6 and CFP-10 antigens. In contrast, RD5mic, a region that contains three phospholipase C genes (plcA to -C), was missing from only the vole isolates and was present in M. microti strains isolated from humans. Apart from RD1mic and RD5mic other M. microti-specific deleted regions have been identified (MiD1 to MiD3). Deletion of MiD1 has removed parts of the direct repeat region in M. microti and thus contributes to the characteristic spoligotype of M. microti strains.