First Whole-Genome Sequence of a Clinical Isolate of Multidrug-Resistant Mycobacterium bovis BCG

The attenuated BCG strain of Mycobacterium bovis is widely used as a vaccine against tuberculosis. However, in rare cases, it can be pathogenic to humans. Here, we report the first draft of a whole-genome sequence of a multidrug-resistant clinical isolate of M. bovis BCG.

rivative of Mycobacterium bovis, has been used to immunize against tuberculosis, as it is currently the only available vaccine. As of 1974, it has been part of the WHO's Expanded Program on Immunization because of its proven efficacy in preventing extrapulmonary tuberculosis in children. However, in adults, its efficacy against pulmonary disease varies, possibly due to the use of different derived daughter strains (1)(2)(3).
The attenuated strain BCG was derived by Albert Calmette and Camille Guérin at the Institut Pasteur from an M. bovis isolate after 230 serial passages between 1908 and 1921 on potato slices soaked with glycerol. Once the safety of the vaccine was confirmed, BCG was disseminated around 1924 to different laboratories that maintained their own daughter strains through passaging, until the introduction of archival seed lots in the 1960s. Consequently, the original strain of BCG has produced many heterogenic daughter strains (substrains), which are the progenitors of the commonly used vaccines (1)(2)(3)(4)(5)(6).
Here, we report the first draft genome sequence of a multidrugresistant (MDR) clinical isolate of M. bovis BCG, BCG-MDR, which was isolated from a patient sputum sample and was found to be phenotypically resistant to rifampin, isoniazid, ethionamide, para-aminosalicylic acid, and cycloserine. Very few BCG strains resistant to rifampin or to both rifampin and isoniazid have been identified so far, and none have been sequenced (12)(13)(14)(15). The genome sequence was obtained using Illumina sequencing (MiSeq technology), producing a total of 2,894,644 paired-end reads. These sequences were assembled using SPAdes version 3.0.0 (16), resulting in 62 contigs. The draft genome sequence is 4,333,148 bp in length, with a GϩC content of 65% and an average coverage of 133ϫ. The genome sequence was annotated using the NCBI Prokaryotic Genomes Automatic Annotation Pipeline (htt p://www.ncbi.nlm.nih.gov/genomes/static/Pipeline.html). Using this automated pipeline, we identified 3,943 predicted coding sequences; out of them, 3,825 had a functional assignment. It also predicted 65 pseudogenes, 3 rRNAs, and 45 tRNAs.
Comparison with available M. bovis BCG genomes indicates that the clinical isolate is related to BCG Tokyo and BCG Russia. Nonsynonymous mutations were identified in genes known to be involved in the development of resistant forms of the M. tuberculosis complex. The availability of this genome will allow comparative genomic analysis with other M. bovis BCG daughter strains and might provide insights into the emergence of MDR clinical strains of M. bovis BCG.
Nucleotide sequence accession numbers. This whole-genome shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession no. JNAF00000000. The version described in this paper is version JNAF01000000.

ACKNOWLEDGMENTS
This work was funded by ESCMID research grant 2013 to A. Renvoisé and by DIM MalInf grant 120085 to S. Pang.