Perspective on sequence evolution of microsatellite locus (CCG)n in Rv0050 gene from Mycobacterium tuberculosis
© Qin et al; licensee BioMed Central Ltd. 2011
Received: 10 June 2011
Accepted: 31 August 2011
Published: 31 August 2011
The mycobacterial genome is inclined to polymerase slippage and a high mutation rate in microsatellite regions due to high GC content and absence of a mismatch repair system. However, the exact molecular mechanisms underlying microsatellite variation have not been fully elucidated. Here, we investigated mutation events in the hyper-variable trinucleotide microsatellite locus MML0050 located in the Rv0050 gene of W-Beijing and non-W-Beijing Mycobacterium tuberculosis strains in order to gain insight into the genomic structure and activity of repeated regions.
Size analysis indicated the presence of five alleles that differed in length by three base pairs. Moreover, nucleotide gains occurred more frequently than loses in this trinucleotide microsatellite. Mutation frequency was not completely related with the total length, though the relative frequency in the longest allele was remarkably higher than that in the shortest. Sequence analysis was able to detect seven alleles and revealed that point mutations enhanced the level of locus variation. Introduction of an interruptive motif correlated with the total allele length and genetic lineage, rather than the length of the longest stretch of perfect repeats. Finally, the level of locus variation was drastically different between the two genetic lineages.
The Rv0050 locus encodes the bifunctional penicillin-binding protein ponA1 and is essential to mycobacterial survival. Our investigations of this particularly dynamic genomic region provide insights into the overall mode of microsatellite evolution. Specifically, replication slippage was implicated in the mutational process of this microsatellite and a sequence-based genetic analysis was necessary to determine that point mutation events acted to maintain microsatellite size integrity while providing genomic diversity.
Simple sequence repeats (SSRs, also known as microsatellites) are repetitive sequence motifs of one to six base pairs (bps) scattered throughout all known genomes . The extensive length differences that may be achieved by microsatellites and their high rate of polymorphism have facilitated their use as molecular markers in epidemiological investigations. However, knowledge is lacking about the mutational mechanism(s) that lead to variations in microsatellite loci.
Gaining a detailed understanding of the features underlying microsatellite genomic structure will aid subsequent interpretations of data from these clinically useful genomic regions. Strand slippage during DNA replication can cause insertion and deletion of repeat units in newly-synthesized nucleotide chains. These events are the most common cause of expansion or contraction of microsatellites [2, 3]. Recent studies have reported differences in rates and patterns of mutation among distinct loci and species; thus, allele size, motif size, genetic lineage, G/C content, functional potential of the transcribed product, and effectiveness of mismatch repair enzymes might all act as mediators of the mutation patterns of such loci [4–7].
Although several models have been proposed to explain the mutation processes that effect microsatellite evolution, they have yet to be confirmed . Studies into genomic evolution exploit the hyper-variable nature of microsatellite sequences to observe mutation events directly . Specifically, pedigree analysis has provided substantial amounts of mutation data for broad ranges of chromosomal loci and organisms [8–10]. However, most of these studies on microsatellites have focused on size variation among alleles, and have not addressed potential sequence variations within otherwise similarly sized alleles.
Sequence analysis is an alternative empirical approach for studying microsatellite evolution. Elucidating the sequence structure of alleles allows for direct comparison with other alleles within a single species or with orthologous loci from different species, effectively allowing for the study of accumulated mutational effects over evolutionary timescales. Intraspecific comparisons that reveal the sequence structure of individual alleles may provide significant insights into the otherwise complex process of maintaining genomic integrity under selective pressure.
Mycobacterial genomes harbor a number of polymorphic microsatellites . Microsatellites in these genomes impart a certain degree of genome plasticity and probably account for many biological functions in the context of pathogen adaptability, virulence and survival . Usually, errors resulting from strand slippage are promptly repaired by a three-enzyme system composed of mutL, mutS, and mutH; however, mycobacterial genomes lack these enzymes . Thus, such genomes serve as interesting systems to investigate the rates of mutations in microsatellites and the existence of regulatory mechanisms that govern microsatellite mutations.
The repetitive CCG sequence located in the Rv0050 gene in Mycobacterium tuberculosis (MTB) is a trinucleotide microsatellite locus (MML0050) that exhibits high a polymorphism rate. The Rv0050 locus encodes a bifunctional penicillin-binding protein (ponA1) [14, 15]. The hyper-variable locus has gained popularity as a Variable Number Tandem Repeat (VNTR) biomarker in epidemiological investigations of MTB strains . However, the mutational mechanisms that are responsible for generating the high levels of variation in the MML0050 locus remain unclear.
In this study, we sought to explain the mutation tempo and mode of the MML0050 locus and its polymorphic nature using clinical isolates from two MTB genotypes: W-Beijing and non-W-Beijing. W-Beijing strains are genetically closely related, present with a characteristic spoligotype pattern, and have enjoyed wide global dispersion [17–19]. They account for 80-90% of the MTB strains isolated from the Beijing area since the 1950s and remain prevalent in other parts of China, including the Ningxia, Shanghai, and Guangdong provinces . Such strains have attracted much research interest due to their reported association with multiple drug resistance, relapses, treatment failure, hypervirulent phenotype in mice, and faster growth rates in human monocytes [21–25]. We sequenced a set of MML0050 locus alleles of different size classes from different MTB families to analyze the stabilizing effect of interrupting motifs in microsatellite regions, the effect of allele length, genetic lineage for the introduction of interruptive motif, and the relation between number of repeat units and mutation frequency.
Sampling and DNA extraction
The MTB reference strain H37Rv was studied along with 461 clinical strains that had been isolated between July 2009 and April 2010 from regions of Eastern China, including Shanghai, Jiangsu, Zhejiang, Shandong, Fujian, Anhui, and Jiangxi (collected by Laboratory of Tuberculosis, Shanghai Pulmonary Hospital, Medical School, Tongji University, Shanghai, China).
Strains were grown in Sauton culture medium (0.5 g/L KH2PO4, 0.5 g/L MgSO4.7H2O, 2 g/L citric acid, 0.05 g/L ferric ammonium citrate, 4.0 g/L L-asparagine, 6% glycerol and 0.02% Tween 80). Cells were sterilized at 80°C for 30 min, and were harvested by centrifugation (12,000 × g for 5 min). The bacterial pellet was washed three times with sterilized saline and re-collected by centrifugation (12,000 × g for 10 min each time). DNA Lysis Buffer (10 mmol/L NaCl, 1 mg/mL SDS, 15% Chelex-100, 1% Tween 20) was added and cells were incubated at 50°C for 1 h, followed by 100°C for 10 min. The mixture was centrifuged (5,000 × g for 10 min) to obtain the aqueous phase containing genomic DNA that was used for PCR amplification.
Identification of W-Beijing strains
W-Beijing strains were identified by deletion-targeted multiplex PCR (DTM-PCR) to detect the genomic deletion RD105, which defines the W-Beijing family as a separate lineage within MTB [26–28]. The DTM-PCR primers P1 (5'-GGAGTCGTTGAGGGTGTTCATCAGCTCAGTC-3') and P2 (5'-CGCCAAGGCCGCATAGTCACGGTCG-3') were designed to amplify a 1466 bp product from the non-W-Beijing strains, while P1 and P3 (5'- GGTTGCCCACTGGTCGATATGGTGGACTT-3') were designed to amplify a 761 bp fragment from the W-Beijing genotype. PCR was performed in a 20 μL reaction mixture containing 2× Long Taq mixture (Tiangen, Co., Beijing, China), 0.2 μM of each primer and 10 ng of template DNA, under the following conditions: denaturation at 94°C for 5 min; 30 amplification cycles of 94°C for 30 s, 68°C for 30 s and 72°C for 2 min; and a final extension at 72°C for 7 min. PCR products were separated on 0.8% agarose gels.
To detect polymorphisms within the MML0050 locus, PCR primer pairs were designed for the flanking sequences: HponAF (5'-TTGAAGGGCACGTCGAACGAG-3') and HponAR (5'-GGGACCGATCGGGATGGTAA-3').
PCR was performed in a 20 μL reaction mixture containing 2× Taq mixture (Tiangen), 0.2 μM of each primer and 10 ng of template DNA, under the following conditions: denaturation at 94°C for 5 min; 35 amplification cycles of 94°C for 30 s, 65°C for 30 s and 72°C for 2 min; and a final extension at 72°C for 7 min. PCR products were analyzed by the PCR-Single Strand Conformation Polymorphism (SSCP) technique to detect point mutations; briefly, DNA was separated on a 6% polyacrylamide gel, visualized by silver staining, excised and sequenced on an ABI 3730xl DNA Analyzer (procedure carried out by Sangon, Shanghai, China).
GraphPad Prism 5.0 statistical software package (http://www.graphpad.com/prism/prism.htm) was used to analyze differences in allele distributions detected among the W-Beijing and non-W-Beijing genotypes. The character of bias was analyzed by Student's t-test with statistical confidence limit of 95%. Allele sequences were compared to sequences previously deposited in the GenBank database (http://www.ncbi.nlm.nih.gov/BLAST/). Phylogenetic trees were constructed by the neighbor-joining method using MEGA 4.0 software to ascertain the evolutionary distance of the locus mutation rates .
Results and Discussion
MML0050 microsatellite features revealed by fragment length analysis
MML0050 microsatellite features revealed by sequence analyses
Alleles of MML0050 locus in the Rv0050 gene from M. tuberculosis
H37Ra, CDC1551, KZN1435
M. Bovis BCG
Variations in the genome may occur by other modes besides insertion/deletion. Nucleotide substitutions, such as transitions (A→G) or transversions (G→T), may alter the genome and its function without affecting the structure itself. Sequence analysis revealed that the Rv0050*01 allele was composed of seven perfect CCG repeats from the MTB reference strains H37Rv. All alleles from the clinical strains had C→T substitutions in the last repeat motif. Rv0050*07 was composed of 10 repeat units and had the largest fragment size; moreover, this allele had G→T substitutions in the fifth repeat motif, which interrupted the CCG consecutive repeat.
Microsatellite alleles can generally be divided into two types: alleles that are identical in both length and sequence and those identical in length but not in sequence, the latter being known as homoplastic microsatellite alleles. Rv0050*05 and Rv0050*06 had been determined by fragment length analyses to be identical in length (both having nine repeat units) but sequence analysis revealed significant differences in the underlying genomic sequences. Rv0050*05 had a perfect repeat with eight CCG units, while Rv0050*06 had a G→T substitution in the fourth repeat motif and the same mutation pattern of Rv0050*07. Overall, only 7.4% of the MML0050 mutations detected in our samples resulted in homoplastic character.
We also analyzed the relation between number of interruptive motifs and (1) the total length of allele or (2) genetic lineage. The Rv0050*06 allele from W-Beijing genotype with nine repeat motifs was determined to be the point of introduction of a second interruptive motif. In contrast, the Rv0050*05 allele belonging to the non-W-Beijing genotype had a perfect nine repeat. Alleles from the non-W-Beijing genotype had a greater amount of perfect repeat motifs than did those from the W-Beijing genotype. Introduction of an interruptive motif had the strongest statistical correlations with the total length of allele and genetic lineage (P < 0.0001, 95% CI: 6.482 to 8.946, t = 13.64). Surprisingly, the length of the longest stretch of perfectly repeated units was not associated with any interruptive motif, suggesting the existence of a threshold level for the maximum length of perfectly repeated trinucleotide motifs in stable MML0050 alleles. In addition, this finding suggests that genetic lineage and point mutation may be moderators of mutational events in those alleles with nine repeat motifs.
Evolution of the MML0050 locus
Undoubtedly, the high GC content (65.6% in MTB)  and absence of mismatch repair system (MMR) make the mycobacterial genome prone to polymerase slippage  and lead to a high mutation rate in simple sequence repeats . However, the hyper-variable trinucleotide microsatellite locus MML0050 locus does not contain conspicuously long simple sequence repeats. Our analysis of W-Beijing and non-W-Beijing strains suggested that 10 repeats represented the upper limit for allele size in the MTB genome. Wanner et al. proposed context-dependent codon choice as an alternative mechanism used by the bacteria to reduce the number of mutations ; however, this mechanism does not fit with trinucleotide-related events. It has also been suggested by several others that infinite microsatellite growth may be disabled by introductions of point mutations [6, 8, 33]. Specifically, if the balance between slippage and point mutation favors point mutations within the repeated region, the mutations may interrupt the feature of the microsatellite without eliciting large changes in length. Indeed, all of the clinical strains examined in our study contained point mutations. This finding has led us to speculate that the first point mutation in the last repeat motif of all sequences may represent one of the bacteria's mechanisms to confine the expansion or contraction of microsatellite sequences. In addition, the second point mutation in the middle of the sequence may act as an alternative way to maintain the stabilization of long alleles by breaking repeat arrays into smaller units.
Trinucleotide repeating sequences from the MML0050 locus are transcribed into a repeating series of prolines in the Rv005 protein sequence. It has been reported that repeating sequences of amino acids can affect the physical and chemical properties of proteins, and harbor the potential for producing gradual and predictable changes in protein action [34, 35]. Penicillin-binding protein encoded by Rv0050 is not only a key cell wall synthesizing enzyme [14, 15] but also plays an important role in regulating cell wall hydrolysis . Depletion of the penicillin-binding protein results in misshapen bacterial cells and impaired growth [15, 36]. However, it has been unclear whether repeating sequences of amino acids caused by the MML0050 locus can affect the protein function or bacterial phenotype and growth. Nevertheless, the results from our study indicated that nearly all samples in the W-Beijing group with unique phenotypes harbored the allele with the highest repeat numbers. It is possible that polymorphism of the MML0050 locus, to some extent, can affect functional properties of the proteins and may, subsequently, impact bacterial phenotype.
Size analysis indicated the presence of five MML0050 alleles that differed in length by three base pairs, implicating replication slippage in the mutational process of this microsatellite. Although the relative frequency in the longest allele was remarkably higher than that in the shortest for all of the tested strains, mutation frequency was not completely related with the total length in this trinucleotide microsatellite. Moreover, nucleotide gains were found to have occurred more frequently than loses. Sequence analysis supported the notion that point mutation events acted to maintain microsatellite size integrity while providing genomic diversity. Introduction of an interruptive motif correlated with the total allele length and genetic lineage, rather than the length of the longest stretch of perfect repeats. Finally, the comparison of W-Beijing vs. non-W-Beijing strains revealed that the levels of the locus variation were significantly different between the two groups.
This work was supported by Ministry of Science and Technology of the People's Republic of China (No.2009CB522605), and Science and Technology Commission of Shanghai Municipality, Shanghai, P. R. China (No.10411955100 and No.11ZR1430200).
- Schlotterer C: Evolutionary dynamics of microsatellite DNA. Chromosoma. 2000, 109: 365-371. 10.1007/s004120000089.View ArticlePubMed
- Levinson G, Gutman GA: Slipped-strand mispairing: a major mechanism for DNA sequence evolution. Mol Biol Evol. 1987, 4: 203-221.PubMed
- Schlotterer C, Tautz D: Slippage synthesis of simple sequence DNA. Nucleic Acids Res. 1992, 20: 211-215. 10.1093/nar/20.2.211.View ArticlePubMedPubMed Central
- Paun O, Horandl E: Evolution of Hypervariable Microsatellites in Apomictic Polyploid Lineages of Ranunculus carpaticola: Directional Bias at Dinucleotide Loci. Genetics. 2006, 174: 387-398. 10.1534/genetics.105.052761.View ArticlePubMedPubMed Central
- Calabrese P, Durrett R: Dinucleotide repeats in the Drosophila and human genomes have complex, length dependent mutation processes. Mol Biol Evol. 2003, 20: 715-725. 10.1093/molbev/msg084.View ArticlePubMed
- Microsatellites Ellegren: Simple sequences with complex evolution. Nat Rev Genet. 2004, 5: 435-445.View Article
- Xu X, Peng M, Fang Z, Xu X: The direction of microsatellite mutations is dependent upon allele length. Nat Genet. 2000, 24: 396-399. 10.1038/74238.View ArticlePubMed
- Brohede J, Primmer CR, Moller A, Ellegren H: Heterogeneity in the rate and pattern of germline mutation at individual microsatellite loci. Nucleic Acids Res. 2002, 30: 1997-2003. 10.1093/nar/30.9.1997.View ArticlePubMedPubMed Central
- Beck NR, Double MC, Cockburn A: Microsatellite evolution at two hypervariable loci revealed by extensive avian pedigrees. Mol Biol Evol. 2003, 20: 54-61. 10.1093/molbev/msg005.View ArticlePubMed
- Heyer E, Puymirat J, Dieltjes P, Bakker E, de Knijff P: Estimating Y chromosome specific microsatellite mutation frequencies using deep rooting pedigrees. Hum Mol Genet. 1997, 6: 799-803. 10.1093/hmg/6.5.799.View ArticlePubMed
- Sreenu VB, Kumar P, Nagaraju J, Nagarajaram HA: Simple sequence repeats in mycobacterial genomes. J Biosci. 2007, 32: 3-15. 10.1007/s12038-007-0002-7.View ArticlePubMed
- Sreenu VB, Kumar P, Nagaraju J, Nagarajaram HA: Microsatellite polymorphism across the M. tuberculosis and M. bovis genomes: implications on genome evolution and plasticity. BMC Genomics. 2006, 7: 78-10.1186/1471-2164-7-78.View ArticlePubMedPubMed Central
- Springer B, Sander P, Sedlacek L, Hardt W, Mizrahi V, Schär P, Böttger EC: Lack of mismatch correction facilitates genome evolution in mycobacteria. Mol Microbiol. 2004, 53: 1601-1609. 10.1111/j.1365-2958.2004.04231.x.View ArticlePubMed
- Bhakta S, Basu J: Overexpression, purification and biochemical characterization of a class A high-molecular-mass penicillin-binding protein (PBP), PBP1* and its soluble derivative from Mycobacterium tuberculosis. J Biochem. 2002, 361: 635-639. 10.1042/0264-6021:3610635.
- Hett EC, Rubin EJ: Bacterial Growth and Cell Division: a Mycobacterial Perspective. Microbiol. Mol Biol Rev. 2008, 72: 126-156. 10.1128/MMBR.00028-07.View Article
- Spurgiesz RS, Quitugua TN, Smith KL, Schupp J, Palmer EG, Cox RA, Keim P: Molecular Typing of Mycobacterium tuberculosis by Using Nine Novel Variable-Number Tandem Repeats across the Beijing Family and Low-Copy-Number IS6110 Isolates. J Clin Microbiol. 2003, 41: 4224-4230. 10.1128/JCM.41.9.4224-4230.2003.View ArticlePubMed
- Bifani PJ, Mathema B, Kurepina NE, Kreiswirth BN: Global dissemination of the Mycobacterium tuberculosis W-Beijing family strains. Trends Microbiol. 2002, 10: 45-52. 10.1016/S0966-842X(01)02277-6.View ArticlePubMed
- Park YK, Shin S, Ryu S, Cho SN, Koh WJ, Kwon OJ, Shim YS, Lew WJ, Bai GH: Comparison of drug resistance genotypes between Beijing and non-Beijing family strains of Mycobacterium tuberculosis in Korea. J Microbiol Methods. 2005, 63: 165-172. 10.1016/j.mimet.2005.03.002.View ArticlePubMed
- Kremer K, Glynn JR, Lillebaek T, Niemann S, Kurepina NE, Kreiswirth BN, Bifani PJ, van Soolingen D: Definition of the Beijing/W lineage of Mycobacterium tuberculosis on the basis of genetic markers. J Clin Microbiol. 2004, 42: 4040-4049. 10.1128/JCM.42.9.4040-4049.2004.View ArticlePubMedPubMed Central
- Li WM, Wang SM, Li CY, Liu YH, Shen GM, Zhang XX, Niu TG, Gao Q, van Soolingen D, Kremer K, Duanmu HJ: Molecular epidemiology of Mycobacterium tuberculosis in China: a nationwide random survey in 2000. Int J Tuberc Lung Dis. 2005, 9: 1314-1319.PubMed
- Hanekom M, van der Spuy GD, Streicher E, Ndabambi SL, McEvoy CR, Kidd M, Beyers N, Victor TC, van Helden PD, Warren RM: A recently evolved sublineage of the Mycobacterium tuberculosis Beijing strain family is associated with an increased ability to spread and cause disease. J Clin Microbiol. 2007, 45: 1483-1490. 10.1128/JCM.02191-06.View ArticlePubMedPubMed Central
- Aristimuño L, España M, Guilarte A, Ramírez C, Rojas J, Gavín P, López-Calleja AI, Lezcano MA, Revillo MJ, Cebollada A, Martín C, Samper S: Multidrug-resistant Mycobacterium tuberculosis Beijing/W genotype in Venezuela. J Med Microbiol. 2007, 56: 1707-1708. 10.1099/jmm.0.47310-0.View ArticlePubMed
- Kubica T, Agzamova R, Wright A, Aziz MA, Rakishev G, Bismilda V, Richter E, Rusch-Gerdes S, Niemann S: The Beijing genotype is a major cause of drug-resistant tuberculosis in Kazakhstan. Int J Tuberc Lung Dis. 2005, 9: 646-653.PubMed
- Richardson M, van Lill SW, van der Spuy GD, Munch Z, Booysen CN, Beyers N, van Helden PD, Warren RM: Historic and recent events contribute to the disease dynamics of Beijing-like Mycobacterium tuberculosis isolates in a high incidence region. Int J Tuberc Lung Dis. 2002, 6: 1001-1011.PubMed
- Kong Y, Cave MD, Zhang L, Foxman B, Marrs CF, Bates JH, Yang ZH: Association between Mycobacterium tuberculosis Beijing/W Lineage Strain Infection and Extrathoracic Tuberculosis: Insights from Epidemiologic and Clinical Characterization of the Three Principal Genetic Groups of M. tuberculosis Clinical Isolates. J Clin Microbiol. 2007, 45: 409-414. 10.1128/JCM.01459-06.View ArticlePubMedPubMed Central
- Chen J, Tsolaki AG, Shen X, Jiang X, Mei J, Gao Q: Deletion-targeted multiplex PCR (DTM-PCR) for identification of Beijing/W genotypes of Mycobacterium tuberculosis. Tuberculosis. 2007, 87: 446-449. 10.1016/j.tube.2007.05.014.View ArticlePubMed
- Rindi L, Lari N, Cuccu B, Garzelli C: Evolutionary pathway of the Beijing lineage of Mycobacterium tuberculosis based on genomic deletions and mutT genes polymorphisms. Infect Genet Evol. 2009, 9: 48-53. 10.1016/j.meegid.2008.09.006.View ArticlePubMed
- Tsolaki AG, Gagneux S, Pym AS, Goguet de la Salmoniere YO, Kreiswirth BN, Van Soolingen D, Small PM: Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis. J Clin Microbiol. 2005, 43: 3185-3191. 10.1128/JCM.43.7.3185-3191.2005.View ArticlePubMedPubMed Central
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4:Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMed
- Metzgar D, Liu L, Hansen C, Dybvig K, Wills C: Domain-level differences in microsatellite distribution and content result from different relative rates of insertion and deletion mutations. Genome Res. 2002, 12: 408-413. 10.1101/gr.198602.View ArticlePubMedPubMed Central
- Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, Gordon SV, Eiglmeier K, Gas S, Barry CE, Tekaia F, Badcock K, Basham D, Brown D, Chillingworth T, Connor R, Davies R, Devlin K, Feltwell T, Gentles S, Hamlin N, Holroyd S, Hornsby T, Jagels K, Krogh A, McLean J, Moule S, Murphy L, Oliver K, Osborne J, et al: Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998, 393: 537-544. 10.1038/31159.View ArticlePubMed
- Wanner RM, Güthlein C, Springer B, Böttger EC, Ackermann M: Stabilization of the genome of the mismatch repair deficient Mycobacterium tuberculosis by context-dependent codon choice. BMC Genomics. 2008, 9: 249-10.1186/1471-2164-9-249.View ArticlePubMedPubMed Central
- Anmarkrud JA, Kleven O, Bachmann L, Lifjeld JT: Microsatellite evolution: Mutations, sequence variation, and homoplasy in the hypervariable avian microsatellite locus HrU10. BMC Evol Biol. 2008, 8: 138-10.1186/1471-2148-8-138.View ArticlePubMedPubMed Central
- Hancock JM, Simon M: Simple sequence repeats in proteins and their significance for network evolution. Gene. 2005, 345: 113-118. 10.1016/j.gene.2004.11.023.View ArticlePubMed
- Moxon ER, Rainey PB, Nowak MA, Lenski RE: Adaptive evolution of highly mutable loci in pathogenic bacteria. Curr Bio. 1994, 4: 24-32. 10.1016/S0960-9822(00)00005-1.View Article
- Hett EC, Chao MC, Rubin EJ: Interaction and Modulation of Two Antagonistic Cell Wall Enzymes of Mycobacteria. PLoS Pathog. 2010, 6: e1001020-10.1371/journal.ppat.1001020.View ArticlePubMedPubMed Central
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.