Nucleotide and phylogenetic analyses of the Chlamydia trachomatis ompA gene indicates it is a hotspot for mutation

BMC Research Notes

George F Sensabaugh¹  Brian W Brunelle² 

[1]Division of Infectious Diseases, School of Public Health, University of California, Berkeley, CA 94720, USA

[2]Food Safety and Enteric Pathogens Research Unit, USDA, ARS, National Animal Disease Center, Ames, IA 50010, USA

关键词: Phylogenetics;    ompA;    MOMP;    Evolution;    Chlamydia trachomatis;   

Others  :  1166738
DOI  :  10.1186/1756-0500-5-53

 received in 2011-09-29, accepted in 2012-01-20,  发布年份 2012

【 摘 要 】

Background

Serovars of the human pathogen Chlamydia trachomatis occupy one of three specific tissue niches. Genomic analyses indicate that the serovars have a phylogeny congruent with their pathobiology and have an average substitution rate of less than one nucleotide per kilobase. In contrast, the gene that determines serovar specificity, ompA, has a phylogenetic association that is not congruent with tissue tropism and has a degree of nucleotide variability much higher than other genomic loci. The ompA gene encodes the major surface-exposed antigenic determinant, and the observed nucleotide diversity at the ompA locus is thought to be due to recombination and host immune selection pressure. The possible contribution of a localized increase in mutation rate, however, has not been investigated.

Results

Nucleotide diversity and phylogenetic relationships of the five constant and four variable domains of the ompA gene, as well as several loci surrounding ompA, were examined for each serovar. The loci flanking the ompA gene demonstrated that nucleotide diversity increased monotonically as ompA is approached and that their gene trees are not congruent with either ompA or tissue tropism. The variable domains of the ompA gene had a very high level of non-synonymous change, which is expected as these regions encode the surface-exposed epitopes and are under positive selection. However, the synonymous changes are clustered in the variable regions compared to the constant domains; if hitchhiking were to account for the increase in synonymous changes, these substitutions should be more evenly distributed across the gene. Recombination also cannot entirely account for this increase as the phylogenetic relationships of the constant and variable domains are congruent with each other.

Conclusions

The high number of synonymous substitutions observed within the variable domains of ompA appears to be due to an increased mutation rate within this region of the genome, whereas the increase in nucleotide substitution rate and the lack of phylogenetic congruence in the regions flanking ompA are characteristic motifs of gene conversion. Together, the increased mutation rate in the ompA gene, in conjunction with gene conversion and positive selection, results in a high degree of variability that promotes host immune evasion.

【 授权许可】

   
2012 Brunelle and Sensabaugh; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

【 参考文献 】

• [1]CDC: Sexually Transmitted Disease Surveillance. Atlanta: US Department of Health and Human Services, Centers for Disease Control and Prevention 2010.
• [2]Mabey DC, Solomon AW, Foster A: Trachoma. Lancet 2003, 362(9379):223-229.
• [3]Caldwell HD, Kromhout J, Schachter J: Purification and partial characterization of the major outer membrane protein of Chlamydia trachomatis. Infect Immun 1981, 31(3):1161-1176.
• [4]Bavoil P, Ohlin A, Schachter J: Role of disulfide bonding in outer membrane structure and permeability in Chlamydia trachomatis. Infect Immun 1984, 44(2):479-485.
• [5]Stephens RS: Genomic autobiographies of Chlamydiae. In Chlamydia: Intracellular Biology, Pathogenesis, and Immunity. Edited by Stephens RS. Washington, DC: American Society for Microbiology; 1999:6-26.
• [6]Brunelle BW, Nicholson TL, Stephens RS: Microarray-based genomic surveying of gene polymorphisms in Chlamydia trachomatis. Genome Biol 2004, 5(6):R42. BioMed Central Full Text
• [7]Jeffrey BM, Suchland RJ, Quinn KL, Davidson JR, Stamm WE, Rockey DD: Genome sequencing of recent clinical Chlamydia trachomatis strains identifies loci associated with tissue tropism and regions of apparent recombination. Infect Immun 2010, 78(6):2544-2553.
• [8]Brunelle BW, Sensabaugh GF: The ompA gene in Chlamydia trachomatis differs in phylogeny and rate of evolution from other regions of the genome. Infect Immun 2006, 74(1):578-585.
• [9]Gomes JP, Nunes A, Bruno WJ, Borrego MJ, Florindo C, Dean D: Polymorphisms in the nine polymorphic membrane proteins of Chlamydia trachomatis across all serovars: evidence for serovar Da recombination and correlation with tissue tropism. J Bacteriol 2006, 188(1):275-286.
• [10]Stothard DR, Toth GA, Batteiger BE: Polymorphic membrane protein H has evolved in parallel with the three disease-causing groups of Chlamydia trachomatis. Infect Immun 2003, 71(3):1200-1208.
• [11]Stothard DR, Boguslawski G, Jones RB: Phylogenetic analysis of the Chlamydia trachomatis major outer membrane protein and examination of potential pathogenic determinants. Infect Immun 1998, 66(8):3618-3625.
• [12]Kubo A, Stephens RS: Characterization and functional analysis of PorB, a Chlamydia porin and neutralizing target. Mol Microbiol 2000, 38(4):772-780.
• [13]Kawa DE, Stephens RS: Antigenic topology of chlamydial PorB protein and identification of targets for immune neutralization of infectivity. J Immunol 2002, 168(10):5184-5191.
• [14]Fitch WM, Peterson EM, de la Maza LM: Phylogenetic analysis of the outer-membrane-protein genes of Chlamydiae, and its implication for vaccine development. Mol Biol Evol 1993, 10(4):892-913.
• [15]Feil E, Zhou J, Maynard Smith J, Spratt BG: A comparison of the nucleotide sequences of the adk and recA genes of pathogenic and commensal Neisseria species: evidence for extensive interspecies recombination within adk. J Mol Evol 1996, 43(6):631-640.
• [16]Gomes JP, Bruno WJ, Nunes A, Santos N, Florindo C, Borrego MJ, Dean D: Evolution of Chlamydia trachomatis diversity occurs by widespread interstrain recombination involving hotspots. Genome Res 2007, 17(1):50-60.
• [17]Joseph SJ, Didelot X, Gandhi K, Dean D, Read TD: Interplay of recombination and selection in the genomes of Chlamydia trachomatis. Biol Direct 2011, 6:28. BioMed Central Full Text
• [18]Stephens RS: Chlamydial evolution: a billion years and counting. In Chlamydial Infections, Proceedings of the Tenth International Symposium on Human Chlamydial Infections: 16-21 June 2002: 2002. Antalya, Turkey: International Chlamydia Symposium; 2002:3-16.
• [19]Lampe MF, Suchland RJ, Stamm WE: Nucleotide sequence of the variable domains within the major outer membrane protein gene from serovariants of Chlamydia trachomatis. Infect Immun 1993, 61(1):213-219.
• [20]Brunham R, Yang C, Maclean I, Kimani J, Maitha G, Plummer F: Chlamydia trachomatis from individuals in a sexually transmitted disease core group exhibit frequent sequence variation in the major outer membrane protein (omp1) gene. J Clin Invest 1994, 94(1):458-463.
• [21]Millman KL, Tavare S, Dean D: Recombination in the ompA gene but not the omcB gene of Chlamydia contributes to serovar-specific differences in tissue tropism, immune surveillance, and persistence of the organism. J Bacteriol 2001, 183(20):5997-6008.
• [22]Millman K, Black CM, Johnson RE, Stamm WE, Jones RB, Hook EW, Martin DH, Bolan G, Tavare S, Dean D: Population-based genetic and evolutionary analysis of Chlamydia trachomatis urogenital strain variation in the United States. J Bacteriol 2004, 186(8):2457-2465.
• [23]Dugan J, Rockey DD, Jones L, Andersen AA: Tetracycline resistance in Chlamydia suis mediated by genomic islands inserted into the chlamydial inv-like gene. Antimicrob Agents Chemother 2004, 48(10):3989-3995.
• [24]Yang Z, Yoder AD: Estimation of the transition/transversion rate bias and species sampling. J Mol Evol 1999, 48(3):274-283.
• [25]Stephens RS, Sanchez-Pescador R, Wagar EA, Inouye C, Urdea MS: Diversity of Chlamydia trachomatis major outer membrane protein genes. J Bacteriol 1987, 169(9):3879-3885.
• [26]Nunes A, Nogueira PJ, Borrego MJ, Gomes JP: Chlamydia trachomatis diversity viewed as a tissue-specific coevolutionary arms race. Genome Biol 2008, 9(10):R153. BioMed Central Full Text
• [27]Santoyo G, Romero D: Gene conversion and concerted evolution in bacterial genomes. FEMS Microbiol Rev 2005, 29(2):169-183.
• [28]Somboonna N, Wan R, Ojcius DM, Pettengill MA, Joseph SJ, Chang A, Hsu R, Read TD, Dean D: Hypervirulent Chlamydia trachomatis clinical strain is a recombinant between lymphogranuloma venereum (L(2)) and D lineages. MBio 2011, 2(3):e00045-e00011.
• [29]Jordan IK, Makarova KS, Wolf YI, Koonin EV: Gene conversions in genes encoding outer-membrane proteins in H. pylori and C. pneumoniae. Trends Genet 2001, 17(1):7-10.
• [30]Futse JE, Brayton KA, Knowles DP Jr, Palmer GH: Structural basis for segmental gene conversion in generation of Anaplasma marginale outer membrane protein variants. Mol Microbiol 2005, 57(1):212-221.
• [31]Nystedt B, Frank AC, Thollesson M, Andersson SG: Diversifying selection and concerted evolution of a type IV secretion system in Bartonella. Mol Biol Evol 2008, 25(2):287-300.
• [32]Rich SM, Sawyer SA, Barbour AG: Antigen polymorphism in Borrelia hermsii, a clonal pathogenic bacterium. Proc Natl Acad Sci USA 2001, 98(26):15038-15043.
• [33]Zhang JR, Norris SJ: Genetic variation of the Borrelia burgdorferi gene vlsE involves cassette-specific, segmental gene conversion. Infect Immun 1998, 66(8):3698-3704.
• [34]Noormohammadi AH, Markham PF, Kanci A, Whithear KG, Browning GF: A novel mechanism for control of antigenic variation in the haemagglutinin gene family of Mycoplasma synoviae. Mol Microbiol 2000, 35(4):911-923.
• [35]Kline KA, Sechman EV, Skaar EP, Seifert HS: Recombination, repair and replication in the pathogenic Neisseriae: the 3 R's of molecular genetics of two human-specific bacterial pathogens. Mol Microbiol 2003, 50(1):3-13.
• [36]Haas R, Meyer TF: The repertoire of silent pilus genes in Neisseria gonorrhoeae: evidence for gene conversion. Cell 1986, 44(1):107-115.
• [37]Liao D: Concerted evolution: molecular mechanism and biological implications. Am J Hum Genet 1999, 64(1):24-30.
• [38]Liao D: Gene conversion drives within genic sequences: concerted evolution of ribosomal RNA genes in bacteria and archaea. J Mol Evol 2000, 51(4):305-317.
• [39]Santoyo G, Martinez-Salazar JM, Rodriguez C, Romero D: Gene conversion tracts associated with crossovers in Rhizobium etli. J Bacteriol 2005, 187(12):4116-4126.
• [40]Fehlner-Gardiner C, Roshick C, Carlson JH, Hughes S, Belland RJ, Caldwell HD, McClarty G: Molecular basis defining human Chlamydia trachomatis tissue tropism. A possible role for tryptophan synthase. J Biol Chem 2002, 277(30):26893-26903.
• [41]Caldwell HD, Wood H, Crane D, Bailey R, Jones RB, Mabey D, Maclean I, Mohammed Z, Peeling R, Roshick C, et al.: Polymorphisms in Chlamydia trachomatis tryptophan synthase genes differentiate between genital and ocular isolates. J Clin Invest 2003, 111(11):1757-1769.
• [42]Sharp PM, Li WH: The rate of synonymous substitution in enterobacterial genes is inversely related to codon usage bias. Mol Biol Evol 1987, 4(3):222-230.
• [43]Mira A, Ochman H: Gene location and bacterial sequence divergence. Mol Biol Evol 2002, 19(8):1350-1358.
• [44]Nunes A, Gomes JP, Mead S, Florindo C, Correia H, Borrego MJ, Dean D: Comparative expression profiling of the Chlamydia trachomatis pmp gene family for clinical and reference strains. PLoS One 2007, 2(9):e878.
• [45]Hudson RE, Bergthorsson U, Ochman H: Transcription increases multiple spontaneous point mutations in Salmonella enterica. Nucleic Acids Res 2003, 31(15):4517-4522.
• [46]Maki H: Origins of spontaneous mutations: specificity and directionality of base-substitution, frameshift, and sequence-substitution mutageneses. Annu Rev Genet 2002, 36:279-303.
• [47]Wright BE: A biochemical mechanism for nonrandom mutations and evolution. J Bacteriol 2000, 182(11):2993-3001.
• [48]Wright BE: Stress-directed adaptive mutations and evolution. Mol Microbiol 2004, 52(3):643-650.
• [49]Aguilera A: The connection between transcription and genomic instability. EMBO J 2002, 21(3):195-201.
• [50]Beletskii A, Bhagwat AS: Transcription-induced mutations: increase in C to T mutations in the nontranscribed strand during transcription in Escherichia coli. Proc Natl Acad Sci USA 1996, 93(24):13919-13924.
• [51]Liu X, Hou J, Liu J: Chlamydial DNA polymerase I can bypass lesions in vitro. Biochem Biophys Res Commun 2006, 345(3):1083-1091.
• [52]Davis LG, Kuehl WM, Battey JF: Basic Methods in Molecular Biology. 2nd edition. Edited by Davis LG, Kuehl WM, Battey JF. Norwalk, CT: Appleton and Lange; 1994:16-21.
• [53]Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, White O, Hickey EK, Peterson J, Utterback T, Berry K, et al.: Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res 2000, 28(6):1397-1406.
• [54]Thompson JD, Gibson TJ, Higgins DG: Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics 2002., Chapter 2(Unit 2 3)
• [55]Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24(8):1596-1599.
• [56]Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25(11):1451-1452.