Variation among Genome Sequences of H37Rv Strains of Mycobacterium tuberculosis from Multiple Laboratories

Abstract

The publication of the complete genome sequence for Mycobacterium tuberculosis H37Rv in 1998 has had a great impact on the research community. Nonetheless, it is suspected that genetic differences have arisen in stocks of H37Rv that are maintained in different laboratories. In order to assess the consistency of the genome sequences among H37Rv strains in use and the extent to which they have diverged from the original strain sequenced, we carried out whole-genome sequencing on six strains of H37Rv from different laboratories. Polymorphisms at 73 sites were observed, which were shared among the lab strains, though 72 of these were also shared with H37Ra and are likely to be due to sequencing errors in the original H37Rv reference sequence. An updated H37Rv genome sequence should be valuable to the tuberculosis research community as well as the broader microbial research community. In addition, several polymorphisms unique to individual strains and several shared polymorphisms were identified and shown to be consistent with the known provenance of these strains. Aside from nucleotide substitutions and insertion/deletions, multiple IS6110 transposition events were observed, supporting the theory that they play a significant role in plasticity of the M. tuberculosis genome. This genome-wide catalog of genetic differences can help explain any phenotypic differences that might be found, including a frameshift mutation in the mycocerosic acid synthase gene which causes two of the strains to be deficient in biosynthesis of the surface glycolipid phthiocerol dimycocerosate (PDIM). The resequencing of these six lab strains represents a fortuitous "in vitro evolution" experiment that demonstrates how the M. tuberculosis genome continues to evolve even in a controlled environment.Publication of the whole genome sequence of the H37Rv strain of Mycobacterium tuberculosis by Stewart Cole and colleagues in 1998 provided a breakthrough in tuberculosis (TB) research (8), leading to insights into the biology, metabolism, and evolution of this infectious pathogen. Large protein families related to fatty acid and polyketide biosynthesis, regulation (e.g., sigma factors and two-component sensor systems), drug efflux pumps and transporters, and the PE_PGRS proteins (a large duplicated family unique to the M. tuberculosis group of mycobacteria) were identified. In addition, transposons, prophage-like elements, and other repetitive and/or mobile genetic elements were identified (18). This genomic information has played an essential role in interpreting gene expression studies, modeling persistence, and identifying essential proteins as putative targets for drug discovery. However, to date the functions of only half of the genes (1,756/ 4,066) have been determined or predicted, and the rest remain annotated as "hypothetical proteins" (6).The H37Rv strain was initially selected for sequencing because it is a widely used laboratory strain that has retained its virulence. H37Rv was initially derived from a clinical isolate, H37, obtained fr...

Analysis of Mutations in Streptomycin-Resistant Strains Reveals a Simple and Reliable Genetic Marker for Identification of the Mycobacterium tuberculosis Beijing Genotype

“…All strains analyzed in our study, including the reference H37Rv strain, had the mutations A363G in rpsL and C299T in gidB, which have been reported as sequencing errors in the sequence of H37Rv included in both the TubercuList database and GenBank accession no. AL123456.2 (15,16). The significant mutations found in the strains analyzed in our study are summarized in Table 3.…”

First Evaluation of GenoType MTBDR plus 2.0 Performed Directly on Respiratory Specimens in Central America

“…Base calling was performed using OLB 1.9.3 (Illumina, Inc.). Genome assembly was performed using custom software that implements a comparative assembly approach (12) by aligning reads with the genome sequence of M. tuberculosis H37Rv (GenBank accession no. NC_000962.2).…”

Spontaneous Phthiocerol Dimycocerosate-Deficient Variants of Mycobacterium tuberculosis Are Susceptible to Gamma Interferon-Mediated Immunity

“…The PDIM-negative H37Rv clone had 151 putative SNPs and 15 putative indels. These sequence alterations included 72 polymorphisms (57 SNPs and 15 indels) that were recently identified by whole-genome resequencing of six H37Rv isolates and that are likely to be errors within the H37Rv reference sequence (14). Analysis of the remaining sequence polymorphisms, excluding those in the highly repetitive GC-rich PE_PGRS coding regions, revealed 22 putative SNPs unique to the PDIM-negative H37Rv clone.…”