title: Single nucleotide polymorphism in the genes of mce1 and mce4 operons of Mycobacterium tuberculosis: analysis of clinical isolates and standard reference strains.
abstract: BACKGROUND
  The presence of four mammalian cell entry (mce) operons in Mycobacterium tuberculosis suggests the essentiality of the functions of the genes in these operons. The differential expression of the four mce operons in different phases of in vitro growth and in infected animals reported earlier from our laboratory further justifies the apparent redundancy for these genes in the genome.Here we investigate the extent of polymorphism in eight genes in the mce1 and mce4 operons of M. tuberculosis from four standard reference strains (H37Rv, H37Ra, LVS (Low Virulent Strain) and BCG) and 112 clinical isolates varying in their drug susceptibility profile, analysed by direct sequencing and Sequenom MassARRAY platform.
  
  
  RESULTS
  We discovered 20 single nucleotide polymorphisms (SNPs) in the two operons. The comparative analysis of the genes of mce1 and mce4 operons revealed that yrbE1A [Rv0167] was most polymorphic in mce1 operon while yrbE4A [Rv3501c] and lprN [Rv3495c] had the highest number of SNPs in the mce4 operon. Of 20 SNPs, 12 were found to be nonsynonymous and were further analysed for their pathological relevance to M. tuberculosis using web servers PolyPhen and PMut, which predicted five deleterious nonsynonymous SNPs. A mutation from proline to serine at position 359 of the native Mce1A protein was most deleterious as predicted by both PolyPhen and PMut servers. Energy minimization of the structure of native Mce1A protein and mutated protein was performed using InsightII. The mutated Mce1A protein showed structural changes that could account for the effects of this mutation.
  
  
  CONCLUSIONS
  Our results show that SNPs in the coding sequences of mce1 and mce4 operons in clinical isolates can be significantly high. Moreover, mce4 operon is significantly more polymorphic than mce1 operon (p < 0.001). However, the frequency of nonsynonymous substitutions is higher in mce1 operon and synonymous substitutions are more in mce4 operon. In silico modeling predict that nonsynonymous SNP at mce1A [Rv0169], a virulence gene could play a pivotal role in causing functional changes in M. tuberculosis that may reflect upon the biology of the bacteria.
authors.affiliation: Department of Microbiology, Vallabhbhai Patel Chest Institute, University of Delhi, Delhi-110007, India.
authors.forename: Rashmi
authors.initials: R
authors.lastname: Pasricha
authors.forename: Amita
authors.initials: A
authors.lastname: Chandolia
authors.forename: Prija
authors.initials: P
authors.lastname: Ponnan
authors.forename: Neeraj Kumar
authors.initials: NK
authors.lastname: Saini
authors.forename: Sangeeta
authors.initials: S
authors.lastname: Sharma
authors.forename: Madhu
authors.initials: M
authors.lastname: Chopra
authors.forename: Mandira Varma
authors.initials: MV
authors.lastname: Basil
authors.forename: Vani
authors.initials: V
authors.lastname: Brahmachari
authors.forename: Mridula
authors.initials: M
authors.lastname: Bose
mesh_terms: D000942:Antigens, Bacterial / Q000235:genetics*; D001426:Bacterial Proteins / Q000235:genetics*; D003062:Codon; D019295:Computational Biology / Q000379:methods; D004269:DNA, Bacterial / Q000235:genetics; D008958:Models, Molecular; D009154:Mutation; D009169:Mycobacterium tuberculosis / Q000235:genetics*; D016571:Neural Networks, Computer; D009876:Operon; D020641:Polymorphism, Single Nucleotide; D017434:Protein Structure, Tertiary; D017422:Sequence Analysis, DNA
chemical_list: D000942:Antigens, Bacterial; D001426:Bacterial Proteins; D003062:Codon; D004269:DNA, Bacterial; C510664:Mce4A protein, Mycobacterium tuberculosis; C443482:mce2 protein, Mycobacterium tuberculosis