A Comprehensive Survey of Single Nucleotide Polymorphisms (SNPs) across Mycobacterium bovis Strains and M. bovis BCG Vaccine Strains Refines the Genealogy and Defines a Minimal Set of SNPs That Separate Virulent M. bovis Strains and M. bovis BCG Strains

ABSTRACT To further unravel the mechanisms responsible for attenuation of the tuberculosis vaccine Mycobacterium bovis BCG, comparative genomics was used to identify single nucleotide polymorphisms (SNPs) that differed between sequenced strains of Mycobacterium bovis and M. bovis BCG. SNPs were assayed in M. bovis isolates from France and the United Kingdom and from different BCG vaccines in order to identify those that arose during the attenuation process which gave rise to BCG. Informative data sets were obtained for 658 SNPs from 21 virulent M. bovis strains and 13 BCG strains; these SNPs showed phylogenetic clustering that was consistent with the geographical origin of the strains and previous schemes for BCG genealogies. The data revealed a closer relationship between BCG Tice and BCG Pasteur than was previously appreciated, while we were able to position BCG Beijing within a grouping of BCG Denmark-derived strains. Only 186 SNPs were identified between virulent M. bovis strains and all BCG strains, with 115 nonsynonymous SNPs affecting important functions such as global regulators, transcriptional factors, and central metabolism, which might impact on virulence. We therefore refine previous genealogies of BCG vaccines and define a minimal set of SNPs between virulent M. bovis strains and the attenuated BCG strain that will underpin future functional analyses.

[1]  Lynne Regan,et al.  TPR proteins: the versatile helix. , 2003, Trends in biochemical sciences.

[2]  N. Smith,et al.  The population structure of Mycobacterium bovis in Great Britain: Clonal expansion , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[3]  M. Behr,et al.  Reduced expression of antigenic proteins MPB70 and MPB83 in Mycobacterium bovis BCG strains due to a start codon mutation in sigK , 2005, Molecular microbiology.

[4]  P. Small,et al.  Unique Gene Expression Profiles in Infants Vaccinated with Different Strains of Mycobacterium bovis Bacille Calmette-Guérin , 2007, Infection and Immunity.

[5]  Tanya Parish,et al.  The senX3-regX3 two-component regulatory system of Mycobacterium tuberculosis is required for virulence. , 2003, Microbiology.

[6]  Julian Parkhill,et al.  The complete genome sequence of Mycobacterium bovis , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Martin C. J. Maiden,et al.  Bioinformatics Applications Note Sequence Type Analysis and Recombinational Tests (start) , 2022 .

[8]  R. Fleischmann,et al.  Modeling Bacterial Evolution with Comparative-Genome-Based Marker Systems: Application to Mycobacterium tuberculosis Evolution and Pathogenesis , 2003, Journal of bacteriology.

[9]  D van Soolingen,et al.  Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology , 1997, Journal of clinical microbiology.

[10]  R. Frothingham,et al.  Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. , 1998, Microbiology.

[11]  M. Behr,et al.  A Point Mutation in the mma3 Gene Is Responsible for Impaired Methoxymycolic Acid Production in Mycobacterium bovis BCG Strains Obtained after 1927 , 2000, Journal of bacteriology.

[12]  D. Collins,et al.  Different susceptibility of two animal species infected with isogenic mutants of Mycobacterium bovis identifies phoT as having roles in tuberculosis virulence and phosphate transport. , 2003, Microbiology.

[13]  M. Behr,et al.  Mycobacterium bovis BCG Vaccines Exhibit Defects in Alanine and Serine Catabolism , 2003, Infection and Immunity.

[14]  Anu Raghunathan,et al.  Comparative genome sequencing of Escherichia coli allows observation of bacterial evolution on a laboratory timescale , 2006, Nature Genetics.

[15]  M. Behr,et al.  The in vitro evolution of BCG vaccines. , 2003, Vaccine.

[16]  Jun Liu,et al.  Differential productions of lipid virulence factors among BCG vaccine strains and implications on BCG safety. , 2007, Vaccine.

[17]  W. Jacobs,et al.  The Two-Component Regulatory System senX3-regX3 Regulates Phosphate-Dependent Gene Expression in Mycobacterium smegmatis , 2007, Journal of bacteriology.

[18]  B. Barrell,et al.  Genomics of Mycobacterium bovis. , 2001, Tuberculosis.

[19]  G. Cook,et al.  The Phn system of Mycobacterium smegmatis: a second high-affinity ABC-transporter for phosphate. , 2006, Microbiology.

[20]  K. McDonough,et al.  The Mycobacterium bovis BCG Cyclic AMP Receptor-Like Protein Is a Functional DNA Binding Protein In Vitro and In Vivo, but Its Activity Differs from That of Its M. tuberculosis Ortholog, Rv3676 , 2007, Infection and Immunity.

[21]  Eduardo P C Rocha,et al.  Comparisons of dN/dS are time dependent for closely related bacterial genomes. , 2006, Journal of theoretical biology.

[22]  C. Hutchinson,et al.  Mapping the DNA‐binding domain and target sequences of the Streptomyces peucetius daunorubicin biosynthesis regulatory protein, DnrI , 2002, Molecular microbiology.

[23]  G. Schoolnik,et al.  Comparative genomics of BCG vaccines by whole-genome DNA microarray. , 1999, Science.

[24]  P. Wheeler,et al.  The pyruvate requirement of some members of the Mycobacterium tuberculosis complex is due to an inactive pyruvate kinase: implications for in vivo growth , 2005, Molecular microbiology.

[25]  C. Leclerc,et al.  Comparison of immune responses of mice immunized with five different Mycobacterium bovis BCG vaccine strains , 1996, Infection and immunity.

[26]  M. Worobey,et al.  Genomics: Anthrax and the art of war (against ascertainment bias) , 2005, Heredity.

[27]  G. Besra,et al.  The embAB genes of Mycobacterium avium encode an arabinosyl transferase involved in cell wall arabinan biosynthesis that is the target for the antimycobacterial drug ethambutol. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[28]  M. Nei,et al.  Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. , 1986, Molecular biology and evolution.

[29]  H. David.  Resistance to D-cycloserine in the tubercle bacilli: mutation rate and transport of alanine in parental cells and drug-resistant mutants. , 1971, Applied microbiology.

[30]  Jason Hinds,et al.  Comparative transcriptomics reveals key gene expression differences between the human and bovine pathogens of the Mycobacterium tuberculosis complex. , 2007, Microbiology.

[31]  G. Besra,et al.  Mycobacterium tuberculosis pks12 Produces a Novel Polyketide Presented by CD1c to T Cells , 2004, The Journal of experimental medicine.

[32]  E. Myers,et al.  Basic local alignment search tool. , 1990, Journal of molecular biology.

[33]  Influence of BCG vaccine strain on the immune response and protection against tuberculosis. , 2008, FEMS microbiology reviews.

[34]  Lippincott-Schwartz,et al.  Supporting Online Material Materials and Methods Som Text Figs. S1 to S8 Table S1 Movies S1 to S3 a " Silent " Polymorphism in the Mdr1 Gene Changes Substrate Specificity Corrected 30 November 2007; See Last Page , 2022 .

[35]  Brian D Robertson,et al.  The OtsAB Pathway Is Essential for Trehalose Biosynthesis in Mycobacterium tuberculosis* , 2005, Journal of Biological Chemistry.

[36]  E. Rubin,et al.  Genes required for mycobacterial growth defined by high density mutagenesis , 2003, Molecular microbiology.

[37]  B. Barrell,et al.  Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence , 1998, Nature.

[38]  G. Kaplan,et al.  The effect of bacille Calmette-Guérin vaccine strain and route of administration on induced immune responses in vaccinated infants. , 2006, The Journal of infectious diseases.

[39]  Richard S. Clifton-Hadley,et al.  Bottlenecks and broomsticks: the molecular evolution of Mycobacterium bovis , 2006, Nature Reviews Microbiology.

[40]  G. Gao,et al.  Novel genome polymorphisms in BCG vaccine strains and impact on efficacy , 2008, BMC Genomics.

[41]  Christopher M. Sassetti,et al.  Genetic requirements for mycobacterial survival during infection , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[42]  R. Barletta,et al.  Roles of Mycobacterium smegmatisd-Alanine:d-Alanine Ligase and d-Alanine Racemase in the Mechanisms of Action of and Resistance to the Peptidoglycan Inhibitor d-Cycloserine , 2003, Antimicrobial Agents and Chemotherapy.

[43]  Julian Parkhill,et al.  Genome plasticity of BCG and impact on vaccine efficacy , 2007, Proceedings of the National Academy of Sciences.

[44]  M. Behr,et al.  Has BCG attenuated to impotence? , 1997, Nature.

[45]  R. Dubos,et al.  Differential characteristics in vitro and in vivo of several substrains of BCG. I. Multiplication and survival in vitro. , 1956, American review of tuberculosis.

[46]  V. Wendisch,et al.  Two-Component Systems of Corynebacterium glutamicum: Deletion Analysis and Involvement of the PhoS-PhoR System in the Phosphate Starvation Response , 2006, Journal of bacteriology.

[47]  Royston Goodacre,et al.  Metabolic fingerprints of Mycobacterium bovis cluster with molecular type: implications for genotype-phenotype links. , 2006, Microbiology.

[48]  Tim W. Overton,et al.  Point mutations in the DNA- and cNMP-binding domains of the homologue of the cAMP receptor protein (CRP) in Mycobacterium bovis BCG: implications for the inactivation of a global regulator and strain attenuation. , 2005, Microbiology.

[49]  Priscille Brodin,et al.  Loss of RD1 contributed to the attenuation of the live tuberculosis vaccines Mycobacterium bovis BCG and Mycobacterium microti , 2002, Molecular microbiology.

[50]  R. Heermann,et al.  The N-terminal Input Domain of the Sensor Kinase KdpD of Escherichia coli Stabilizes the Interaction between the Cognate Response Regulator KdpE and the Corresponding DNA-binding Site* , 2003, Journal of Biological Chemistry.

[51]  Molly Strom,et al.  Single Nucleotide Polymorphisms That Cause Structural Changes in the Cyclic AMP Receptor Protein Transcriptional Regulator of the Tuberculosis Vaccine Strain Mycobacterium bovis BCG Alter Global Gene Expression without Attenuating Growth , 2008, Infection and Immunity.

[52]  Lorenz Wernisch,et al.  Dissection of the heat-shock response in Mycobacterium tuberculosis using mutants and microarrays. , 2002, Microbiology.

[53]  D. Sherman,et al.  Deletion of RD1 from Mycobacterium tuberculosis mimics bacille Calmette-Guérin attenuation. , 2003, The Journal of infectious diseases.

[54]  Stefan Niemann,et al.  High Functional Diversity in Mycobacterium tuberculosis Driven by Genetic Drift and Human Demography , 2008, PLoS biology.

[55]  D. Smith,et al.  The protective and allergenic potency of four BCG substrains in use in China determined in two animal models. , 1988, Tubercle.

[56]  G. Besra,et al.  Molecular structure of EmbR, a response element of Ser/Thr kinase signaling in Mycobacterium tuberculosis. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[57]  G. Mahairas,et al.  Molecular analysis of genetic differences between Mycobacterium bovis BCG and virulent M. bovis , 1996, Journal of bacteriology.

[58]  A. Paterson.  The production of bovine tuberculoprotein. , 1948, The Journal of comparative pathology and therapeutics.